Ecosyste.ms: Awesome
An open API service indexing awesome lists of open source software.
https://github.com/scivision/signal_subspace
vector and ensemble signal analysis w/ subspace methods such as Esprit and RootMusic
https://github.com/scivision/signal_subspace
esprit f2py frequency-estimation modern-fortran music root-music
Last synced: about 2 months ago
JSON representation
vector and ensemble signal analysis w/ subspace methods such as Esprit and RootMusic
- Host: GitHub
- URL: https://github.com/scivision/signal_subspace
- Owner: scivision
- License: mit
- Created: 2016-03-13T08:30:35.000Z (almost 9 years ago)
- Default Branch: main
- Last Pushed: 2024-05-10T03:20:00.000Z (8 months ago)
- Last Synced: 2024-09-29T23:12:42.919Z (3 months ago)
- Topics: esprit, f2py, frequency-estimation, modern-fortran, music, root-music
- Language: Python
- Homepage:
- Size: 360 KB
- Stars: 24
- Watchers: 8
- Forks: 6
- Open Issues: 0
-
Metadata Files:
- Readme: README.md
- License: LICENSE.txt
Awesome Lists containing this project
README
# Signal Subspace
![ci_python](https://github.com/scivision/signal_subspace/workflows/ci_python/badge.svg)
![ci_cmake](https://github.com/scivision/signal_subspace/workflows/ci_cmake/badge.svg)[![PyPi Download stats](http://pepy.tech/badge/signal_subspace)](http://pepy.tech/project/signal_subspace)
1-D and ensemble signal subspace analysis with methods such as Esprit
and RootMusic in Fortran, C, and Pythonbased in part upon the
[Spectral Analysis Lib public domain code](https://github.com/vincentchoqueuse/spectral_analysis_project)The core subspace code is written in Fortran 2008 and is called from other languages (Python, C).
Since the programs are Fortran / Python based, they should compile and run in virtually any platform from embedded to supercomputer.```sh
cmake -B buildcmake --build build
```If Lapack is not available, it is built automatically.
In 2023-2024, F2PY and Numpy are going through a transition of build systems.
[CMake script build the f2py bindings](https://numpy.org/doc/stable/f2py/buildtools/cmake.html)
are used to build the f2py targets that allow Python use of this library.```sh
cmake -B build -Dpython=yescmake --build build
```In the examples below, observe the frequency estimates printed along with their corresponding eigenvalues.
A larger eigenvalue is increased confidence in that particular frequency estimate.## ESPRIT examples
All example use a noisy sinusoid.
Some are "complex" using complex numbers, while others use real numbers.There are two versions of the complex program, one a full accuracy using `double complex` numbers, and the other using `single real` numbers as input.
The single real (4 bytes/number) runs about 4 times faster than the double complex (16 bytes/number) program.```sh
build/f_esprit_cmplbuild/f_esprit_real
```Call Fortran Esprit from C using real single precision float:
```sh
build/c_esprit
```Call Fortran Esprit from C++ using real single precision float:
```sh
build/cpp_esprit
```---
Plots comparing real vs. complex ESPRIT:
```sh
python BasicEspritExample.py
```