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!