Bispectrum Indirect FFT method
7 次查看(过去 30 天)
显示 更早的评论
1. Mostly, the bispectrum were computed using the indirect FFT method (I hope I am right) ([bspec, waxis] = bispeci(y, nlag, samp_seg, overlap, flag, wind)). Here, I do not understand what “nlag” specifies. In the toolbox, it is mentioned as “number of cumulant lags to be computed”. Also I do not understand the term “flag” (biased or unbiased). Can you please explain with an example? How to set these parameters? Apart from default window, is there any way to do our analysis with other windows (e.g., hanning)?
0 个评论
回答(1 个)
Wayne King
2013-9-29
编辑:Wayne King
2013-9-29
This function estimates the bispectrum by first estimating the third-order cumulants of the random process, x(t), which is formally
E\{x*(t)x(t+k)x(t+l)\}
The asterisk denotes complex conjugation. k and l are lags.
'biased' (the default) means that the average in the expectation is obtained by dividing by the number of samples 'unbiased' means the average in the expectation is obtained by dividing by the number of samples-1
Look at biased vs. unbiased estimates in statistics and you'll see where that comes from. I would recommend using 'biased' in this case for the same reasons you should use biased in the autocorrelation sequence estimates. This is well documented in the literature.
Without the application of a smoothing window, the bispectrum estimate is an inconsistent estimator of the true bispectrum. The nlags parameter determines how many samples are used in the smoothing window. I would start with approximately 1/4 of your input time series length as a first pass.
另请参阅
类别
在 Help Center 和 File Exchange 中查找有关 Dynamic System Models 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!