BiSpectrum Routine

Summary

Compute a part of the bispectrum.

Declaration

Procedure BiSpectrum(Bispectrum, Spectrum: TVec; Index: integer; Zeroed: boolean);

Description

Compute a single horizontal bispectral line (triple product) from Spectrum and place the result in Bispectrum. The index parameter defines for which frequency you want to compute the triple product. The triple product will computed for one frequency against all others. When used together with RunningAverage routine, bicoherence can be estimated. Bicoherence has 1 where a frequency pair is synchronized in phase. Frequency pairs which are unrelated in phase have a zero in the bicoherence matrix.