View source for Frequency spectra and cross-spectra of shear and vibrations

Auto- and cross-spectra should be calculated using the windowed and overlapped periodogram method.
 
The window should be a periodic cosine window. The overlap should be 50%. Although this is not the optimal overlap, it is a convenient one and it is one that provides 90% of the maximum degrees of freedom attainable with the cosine window 
<ref> Nuttall, A., (1971). Spectral estimation by means of overlapped fast fourier transform processing of windowed data. NUSC Tech. Rep. 4169. [Available online at https://apps.dtic.mil/dtic/tr/fulltext/u2/a182402.pdf] </ref>. 
The number of fft-segments within a diss-length segment of data is <math> N_f = 2 \frac{L_D}{L_f}-1 </math> where <math>L_D</math> is the length of a dissipation rate estimate (diss-length) and <math>L_f</math> is the length of an fft-segment (fft-length).

The auto-spectrum is given by  <math> E_i = \frac{2}{MN_ff_s} \sum^{N_f}_{k=1} |\mathscr{F}_{ik}|^2 </math> where <math> \mathscr{F}_{ik} </math> is the fast Fourier transform of the k-th fft-segment, i is the frequency index ranging from <math> i=1</math> to <math> i=1+M/2</math>, <math>M</math> is the number of samples in an fft-segment, <math>N_f</math> is the number of fft-segments within a dissipation segment, and <math>f_s</math> is the sampling rate of the data. 
The frequency of each spectral value is  <math>f_i = \frac{i-1}{M} f_s</math>.
It is assumed that the cosine window applied to each fft-segment has mean-square of 1, i.e., it does not change the variance of the windowed data with respect to the original data. 
A suitable cosine window is <math>W = \sqrt{\frac{2}{3}}\, \left(1 + \cos\left(\frac{\pi}{M}\, [2j-2-M]\right) \right) </math> where <math> j = 1, ..., M</math> is the index to the samples in an fft-segment.

This spectrum, <math>E_i</math> has the property that its integral (obtained by a suitable numerical approximation) equals the variance of the signal. The very first spectral value, <math>E_1</math> represents the mean and should be set to zero. The very last spectral value represents the Nyquist frequency and should be very small if the data have been properly sampled using an anti-aliasing filter.  

Cross-spectra are calculated similarly to the auto-spectra, except that the <math>|\mathscr{F}_{ik}|^2</math> is replaced by <math>\mathscr{F}_{xik} \mathscr{F}^*_{yik} </math> where <math>\mathscr{F}_{xik}</math> and <math>\mathscr{F}^*_{yik} </math> are the fast Fourier transforms of the two signals x and y and the asterisk indicates a complex conjugate. Cross-spectra are complex and are used in the vibration-coherent noise removal algorithm of [https://journals.ametsoc.org/view/journals/atot/23/7/jtech1889_1.xml Goodman et al (2006)]<ref> Goodman, L., Levine, E. R., & Lueck, R. G. (2006). On measuring the terms of the turbulent kinetic energy budget from an AUV. Journal of Atmospheric and Oceanic Technology, 23(7), 977-990. </ref>.


==References==
<references />


----------------------------
return to [[Flow chart for shear probes]]

[[Category:Shear probes]]