Fourier Spectral Analysis Object and Template – Periodogram (Spectral Analysis Option)

23.08.2021

The Periodogram spectral procedure computes an averaged frequency spectrum by taking the individual FFTs of multiple (and usually overlapping) segments of the data stream. The segmenting results in a smaller size data record, and consequently in a reduced spectral resolution. However, the averaging reduces the variance that would arise from using only one FFT. Data are assumed stationary. Use the Short-Time Fourier Transform spectral analysis to check stationarity.

The individual FFTs are normalized to match the input data. The overall power of the averaged spectrum will exactly match that of the input data only when each segment contains the same power.

Spectrum Type

The frequency domain information can be output in a variety of formats. In the following table, Re is the real component of a given segment's real (single sided) FFT at a given frequency, Im is the imaginary component, δF is the frequency spacing of the spectrum, n is the data set size, δX is the sampling interval, and σ² is the data set variance.

Spectrum Type

Formula/Description

Amplitude

sqrt(Re² + Im²) / n

RMS Amplitude

sqrt((Re² + Im²) / 2) / n

Amplitude²

(Re² + Im²) / n²

dB

20 * log10(sqrt(Re² + Im²) / n / Aref)

Aref = Reference amplitude, which is assigned 0 dB

dB normalized

20 * log10(sqrt(Re² + Im²) / n) - dBmax

dBmax = dB value of the spectral line with maximum amplitude

PSD - Power Spectral Density

(Re² + Im²) / n² / δF / 2

TISA - Time Integral Amplitude²

δX * (Re² + Im²) / n / 2

MSA - Mean Amplitude²

(Re² + Im²) / n² / 2

SSA - Sum Amplitude²

(Re² + Im²) / n / 2

Variance

(Re² + Im²) / (n * σ²) / 2

Magnitude²

Re² + Im²

Magnitude

sqrt(Re² + Im²)

Third octaves (mean values)

The mean of the amplitudes in a third octave band is formed.

Third octaves (sums)

The amplitudes in a third octave band are summed up.

Third octaves (RMS)

The square mean or RMS for each third octave band is calculated, i.e. the squared amplitudes are averaged and from this the square root is calculated.

Third octaves (quadratic means)

The square amplitudes in a third octave band are summed up.

Octaves (mean values)

The mean of the amplitudes in an octave band is formed.

Octaves (sums)

The amplitudes in an octave band are summed up.

Octaves (RMS)

The square mean or RMS for each octave band is calculated, i.e. the squared amplitudes are averaged and from this the square root is calculated.

Octaves (quadratic mean values)

The square amplitudes in an octave band are summed up.

For the dB normalized type, the averaged spectrum will have a maximum at 0 dB. However, the associated peak will probably be slightly positive as a consequence of the bin interpolation.

Note also that averaging dB values is not an arithmetic average of the individual spectra, but more of a logarithmically weighted average. With a simple average of power, intermittent harmonics or noise bursts in one or more segments can overwhelm an otherwise low power trend. When dB values are averaged, the influence of intermittent elements is significantly lessened. A dB average is thus a robust measure of power, although it may miss a harmonic that appears only briefly in time. An arithmetic average in power, or even amplitude, is more likely to catch an intermittent harmonic. This is one of the reasons it is important to first use the Short-Time Fourier Transform (STFT) or Continuous Wavelet Transform (CWT) analysis object when stationarity is not assured.

For the third octaves and octaves spectral types an amplitude spectrum is computed first, which is then evaluated with the ThirdOctaveAnalysis or OctaveAnalysis FPScript functions respectively.

Windows

FlexPro offers a variety of tapering windows to reduce the spectral leakage. The Window adjustment field is used to set the spectral width, and thus the dynamic range, of adjustable windows. This field will be disabled for fixed windows.

The list box Normalization offers two options to normalize after applying the window. Selecting Amplitude normalizes to the gain of the used window function, i.e. the sum of all values is divided by their number. This compensates the damping of the amplitudes caused by applying the window. This is especially useful to measure peaks within the spectrum. If you select Power the loss of power is compensated. The ratio of the sum of the squared data before and after applying the window is used as a normalization factor. The total power within the spectrum therefore always corresponds to the power of the data before applying the window.

Parameters

The Best Exact N composite algorithm is used for the FFT.

The length of individual data segments, Segment Length and the amount of overlap, Overlap in % can be specified. You should set the segment size based on the resolution needed. You can enter 0 for the segment length to set it to the data length / 4. Values for the overlap that produce the minimum variance are reported to be in the range of 50 to 70%. With the setting Gap in samples the data continue to be ignored. This setting should only be selected for very long time series with slowly changing spectral content.

To accommodate zero padding, the FFT Length can be specified separately. Zero padding occurs when you set the FFT length to a value greater than the segment length. You can enter 0 for the FFT length to set it to the segment length. When a data tapering window is used, then zero-padding causes very little spectral leakage. Zero-padding is especially useful for interpolating peak frequencies with this algorithm, given the loss in resolution incurred by the reduced size of the segments.

Options - Peaks (Analysis Wizard Only)

The spectral peaks are identified by a local maxima detection algorithm. Both the amplitude and the frequency locations of the detected peaks are based upon a cubic spline bin interpolation procedure.

The peaks can be set with a maximum peak count or a dB threshold below the largest peak. Peaks are ranked by interpolated amplitude. Note that a target signal component count may not be realized as fewer peaks than this target may be detected.

You can view the Y  and/or X values of the peaks in the spectrum by pressing Toggle Labels.

Options - Set/Clear Reference (Analysis Wizard Only)

This function lets you compare various spectral procedures and settings. You can view a copy of the currently displayed spectrum in the lower pane by pressing Set Reference. Next, you can adjust additional settings that affect the display in the upper pane. With Clear Reference you can remove the copy and the time signal will appear again.

Harmonic Table (Analysis Wizard Only)

The option Optional numeric results on page three of the Analysis Wizard produces a table of frequencies, amplitudes, amplitudes standard deviations (SD) and PSD for all peaks in the spectrum. Absolute and relative percentages are also given for the component powers. These are often the quantities of interest when comparing strengths of signal components.

FPScript Functions Used

Periodogram

OctaveAnalysis

ThirdOctaveAnalysis

See Also

Analysis Objects

Spectral Analysis Option

Fourier Spectral Analysis Object

FFT Algorithms

Fourier Spectral Analysis

Data Tapering Window

Significance Levels

FFTn Function

Fourier Spectral Analysis Tutorial

Share article or send as email:

You might be interested in these articles