Spectrum continuation analysis
2617956
134864639
2007-05-31T18:01:18Z
Taggard
4523507
[[Wikipedia:WikiProject Stub removal|You can help!]]
'''Spectrum continuation analysis (SCA)''' is a generalization of the concept of [[Fourier series]] to '''non-periodic''' functions of which only a fragment has been sampled in the time domain.
Recall that a Fourier series is only suitable to the analysis of '''periodic''' (or finite-domain) functions ''f''(''x'') with period 2π. It can be expressed as an infinite series of sinusoids:
:<math>f(x) = \sum_{n=-\infty}^{\infty} F_n \,e^{inx}</math>
where <math>F_n</math> is the amplitude of the individual harmonics.
'''In SCA however''', one decomposes the spectrum into optimized discrete frequencies. As a consequence, and as the period of the sampled function is supposed to be infinite or not yet known, each of the discrete periodic functions that compose the sampled function fragment can not be considered to be a multiple of the fundamental frequency:
:<math>f(x) = \sum_{n=-\infty}^{\infty} F_n \,e^{i \omega_n x}.</math>
As such, SCA does not necessarily deliver <math> 2 \pi </math> periodic functions, as would have been the case in Fourier analysis.
For real-valued functions, the SCA series can be written as:
:<math>f(x) = \sum_{n=0}^\infty\left[A_n \cos(\omega_n x)+B_n \sin(\omega_n x)\right] + C(x) </math>
where ''A''<sub>''n''</sub> and ''B''<sub>''n''</sub> are the series amplitudes. The amplitudes can only be solved if the series of values <math>\omega_n</math> is previously optimized for a desired objective function (usually least [[errors and residuals in statistics|residuals]]).
<math>C(x)</math> is not necessarily the average value over the sampled interval: one might prefer to include predominant information on the behavior of the offset value in the time domain.
== Etymology ==
SCA deals with the prediction problem of continuing a frequency spectrum beyond a sampled (usually [[stochastic]]) time series fragment. Unlike ordinary Fourier analysis that infinitely repeats an observed function period or time domain, SCA filters the exact composing frequencies out of the observed spectrum and let them continue (resp. precede) in the time domain.
In the scientific terminology, therefore preference is given to the term ''continuation'' rather than for instance ''extrapolation''.
== Algorithm ==
An algorithm is required to cope with several problems: detrending, decomposition, frequency resolution optimization, superposition, transformation and computational efficiency.
* '''Detrending or trend estimation'''.
* '''Decomposition'''.
Since [[discrete Fourier transform]] is inherently related to Fourier analysis, this type of spectral analysis is by definition not suitable for spectrum decomposition in SCA. DFT (or [[Fast Fourier transform|FFT]]) may provide however an initial approximation, which often speeds up the decomposition.
* '''Improving frequency resolution'''.
After decomposition of a discrete frequency, it should be filtered for optimal resolution (i.e. varying three parameters: frequency value, amplitude and phase).
* '''Transformation'''.
== Spectrum dispersion ==
Compared to [[DFT]] (or [[Fast Fourier transform|FFT]]), which is characterized by perfect spectral resolution, but poor temporal information, SCA favours temporal information, but yields higher spectrum dispersion. This property shows where the analytic strength of SCA is located. For instance, discrete composing frequency resolution is by definition far better in SCA than in DFT.
[[Category:Fourier analysis]]
[[Category:Mathematical series]]
[[Category:Digital signal processing]]
[[Category:Transforms]]
[[Category:Functional analysis]]
[[Category:Econometrics]]
[[nl:Spectrum continuatie analyse]]