Source: http://ab-initio.mit.edu/wiki/index.php?title=Harminv&diff=881&oldid=876
Timestamp: 2019-04-19 07:05:24+00:00

Document:
* [http://www.chem.uci.edu/people/faculty/mandelsh/ V. A. Mandelshtam] and H. S. Taylor, "Harmonic inversion of time signals," ''J. Chem. Phys.'' '''107''' (17), 6756-6769 (1997). Erratum, ''ibid.'' '''109''' (10), 4128 (1998). * [http://www.chem.uci.edu/people/faculty/mandelsh/ V. A. Mandelshtam] and H. S. Taylor, "Harmonic inversion of time signals," ''J. Chem. Phys.'' '''107''' (17), 6756-6769 (1997). Erratum, ''ibid.'' '''109''' (10), 4128 (1998).
- This kind of spectral analysis has wide applications in many areas of physics and engineering, as well as other fields. For example, it could be used to extract the vibrational or "eigen" modes of a system from its response to some stimulus, and also their rates of decay in dissipative systems. FDM has been applied to analyze, e.g., NMR experimental data. It is especially appropriate for analyzing numerical simulations, e.g. of quantum mechanics or [[Meep|classical electromagnetism]]. In general, it is useful when you know on physical grounds that your system consists of a small number of decaying &amp; oscillating modes in the bandwidth of interest, and is '''not appropriate to analyze more arbitrary waveforms'''. + This kind of spectral analysis has wide applications in many areas of physics and engineering, as well as other fields. For example, it could be used to extract the vibrational or "eigen" modes of a system from its response to some stimulus, and also their rates of decay in dissipative systems. FDM has been applied to analyze, e.g., NMR experimental data. It is especially appropriate for analyzing numerical simulations, e.g. of quantum mechanics or [[Meep|classical electromagnetism]]. In general, it is useful when you know on physical grounds that your system consists of a small number of decaying &amp; oscillating modes in the bandwidth of interest, plus a limited amount of noise, and is '''not appropriate to analyze completely arbitrary waveforms'''.
Harminv was developed by [http://math.mit.edu/~stevenj/ Steven G. Johnson] ([mailto:stevenj@alum.mit.edu stevenj@alum.mit.edu]), and is [http://www.gnu.org/philosophy/free-sw.html free software] that should [[Harminv installation|install]] under any Unix-like operating system (''e.g.'' GNU/Linux). Harminv was developed by [http://math.mit.edu/~stevenj/ Steven G. Johnson] ([mailto:stevenj@alum.mit.edu stevenj@alum.mit.edu]), and is [http://www.gnu.org/philosophy/free-sw.html free software] that should [[Harminv installation|install]] under any Unix-like operating system (''e.g.'' GNU/Linux).
Harminv is a free program (and accompanying library) to solve the problem of harmonic inversion — given a discrete-time, finite-length signal that consists of a sum of finitely-many sinusoids (possibly exponentially decaying) in a given bandwidth, it determines the frequencies, decay constants, amplitudes, and phases of those sinusoids.
It can, in principle, provide much better accuracy than straightforwardly extracting FFT peaks, essentially because it assumes a specific form for the signal. (Fourier transforms, in contrast, attempt to represent any data as a sum of sinusoidal components, and are thus limited by the uncertainty principle.) It is also often more robust than directly least-squares fitting the data (which can have problematic convergence), since it re-expresses the problem in terms of simply finding the eigenvalues of a small matrix.
V. A. Mandelshtam and H. S. Taylor, "Harmonic inversion of time signals," J. Chem. Phys. 107 (17), 6756-6769 (1997). Erratum, ibid. 109 (10), 4128 (1998).
This kind of spectral analysis has wide applications in many areas of physics and engineering, as well as other fields. For example, it could be used to extract the vibrational or "eigen" modes of a system from its response to some stimulus, and also their rates of decay in dissipative systems. FDM has been applied to analyze, e.g., NMR experimental data. It is especially appropriate for analyzing numerical simulations, e.g. of quantum mechanics or classical electromagnetism. In general, it is useful when you know on physical grounds that your system consists of a small number of decaying & oscillating modes in the bandwidth of interest, plus a limited amount of noise, and is not appropriate to analyze completely arbitrary waveforms.
Harminv was developed by Steven G. Johnson (stevenj@alum.mit.edu), and is free software that should install under any Unix-like operating system (e.g. GNU/Linux).
What's new in each version is described in the Harminv release notes. Harminv is distributed under the GNU GPL and comes with NO WARRANTY (see the license for more details).
It would be courteous of you to cite Harminv and its author in any publication for which you find it useful, in addition to citing a Mandelshtam reference (either the one above or the review article below).
To install, please see our installation instructions. The main requirement is that you first install the BLAS and LAPACK libraries.
Most users will employ Harminv via the standalone harminv program, which is described by the harminv man page. For instructions on calling the libharminv library directly from your own C/C++ code, see the README file included with Harminv.
Essentially, the FDM algorithm works by considering the time-series to be the autocorrelation function of a fictitious dynamical system, such that the problem of finding the frequencies and decay constants is re-expressed as the problem of finding the eigenvalues of the complex-symmetric time-evolution operator of this system. The key point is that, if you are only interested in frequencies within a known band-limited region, the matrix elements of this operator can be expressed purely in terms of Fourier transforms (or, really, z transforms) of your time-series. Then, one can simply diagonalize a small matrix (size proportional to the bandwidth and the number of frequencies) to find the desired result.
H. Hu, Q. N. Van, V. A. Mandelshtam, and A. J. Shaka, "Reference deconvolution, phase correction, and line listing of NMR spectra by the 1D filter diagonalization method," J. Magnetic Resonance 134, 76-87 (1998).
Rongqing Chen and Hua Guo, "Efficient calculation of matrix elements in low storage filter diagonalization," J. Chem. Phys. 111 (2), 464-471 (1999).
J. Chen and V. A. Mandelshtam, "Multiscale filter diagonalization method for spectral analysis of noisy data with nonlocalized features," J. Chem. Phys. 112 (10), 4429-4437 (2000).
V. A. Mandelshtam, "On harmonic inversion of cross-correlation functions by the filter diagonalization method," J. Theoretical and Computational Chemistry 2 (4), 497-505 (2003).
...as well as other variants, such as multidimensional FDM, described in the 2001 review article above.

References: V. 
 V. 

V. 
 V. 
 V. 

V.