|
|
| Line 1: |
Line 1: |
| In [[mathematics]], in the area of [[statistical analysis]], '''bicoherence''' is a squared normalised version of the [[bispectrum]]. The bicoherence takes values bounded between 0 and 1, which make it a convenient measure for quantifying the extent of phase coupling in a signal. It is also known as '''bispectral coherency'''. The prefix ''bi-'' in ''bispectrum'' and ''bicoherence'' refers not to two time series ''x''<sub>''t''</sub>, ''y''<sub>''t''</sub> but rather to two frequencies of a single signal.
| | Greetings. The author's title is Phebe and she feels comfortable when individuals use the [http://my.clevelandclinic.org/disorders/venereal_diseases/hic_genital_warts.aspx complete] [http://computerrepairredlands.com/solid-advice-terms-candida-albicans/ visit the up coming article] name. Hiring is her day job now and she will not alter it whenever quickly. Doing ceramics is what her family members and her enjoy. Puerto Rico is where he's always been residing but she std testing at home std [http://3bbc.com/index.php?do=/profile-548128/info/ over the counter std test] test kit home [http://Www.Feminist.com/resources/ourbodies/std.html requirements] to move because of her family.<br><br>Feel free to surf to my homepage home std test kit [[http://nxnn.info/blog/284637 visit this website link]] |
| | |
| The ''bispectrum'' is a statistic used to search for nonlinear interactions. The [[Fourier transform]] of the second-order [[cumulant]], i.e., the [[autocorrelation]] function, is the traditional [[power spectrum]]. The Fourier transform of C<sub>3</sub>(t<sub>1</sub>,t<sub>2</sub>) (third-order [[cumulant]]) is called bispectrum or '''bispectral density'''. They fall in the category of ''Higher Order Spectra'', or ''Polyspectra'' and provide supplementary information to the power spectrum. The third order polyspectrum (bispectrum) is the easiest to compute, and hence the most popular.
| |
| | |
| The difference with measuring ''coherence'' (coherence analysis is an extensively used method to study the correlations in frequency domain, between two simultaneously measured signals) is the need for both input and output measurements by estimating two auto-spectra and one cross spectrum. On the other hand, bicoherence is an auto-quantity, i.e. it can be computed from a single signal. The coherence function provides a quantification of deviations from linearity in the system which lies between the input and output measurement sensors. The bicoherence measures the proportion of the signal energy at any bifrequency that is quadratically phase coupled. It is usually normalized in the range similar to correlation coefficient and classical (second order) coherence. It was also used for depth of anasthesia assesement and widely in plasma physics (nonlinear energy transfer) and also for detection of gravitation waves.
| |
| | |
| Bispectrum and bicoherence may be applied to the case of non-linear interactions of a continuous spectrum of propagating waves in one dimension [http://www.iop.org/EJ/abstract/0741-3335/30/5/005].
| |
| | |
| Bicoherence measurements have been carried out for [[Electroencephalography|EEG]] [[signals (biology)|signals]] monitoring in [[sleep]], [[Awake|wakefulness]] and [[seizures]] [http://cogprints.org]. Search for papers by Bullock and Achimowicz.
| |
| | |
| == Definition ==
| |
| The bispectrum is defined as the triple product
| |
| :<math>B(f_1,f_2) = F(f_1)F(f_2)F^*(f_1+f_2)</math>
| |
| where <math>B</math> is the bispectrum evaluated at frequencies <math>f_1</math> and <math>f_2</math>, <math>F</math> is the Fourier transform of the signal, and <math>^*</math> denotes the complex conjugate. The Fourier transform is a complex quantity, and so is the bispectrum. From complex multiplication, the magnitude of the bispectrum is equal to the product of the magnitudes of each of the frequency components, and the phase of the bispectrum is the sum of the phases of each of the frequency components.
| |
| | |
| Suppose that the three Fourier components <math>F(f_1)</math>, <math>F(f_2)</math> and <math>F(f_1 + f_2)</math> were perfectly phase locked. Then if the Fourier transform was calculated several times from different parts of the time series, the bispectrum will always have the same value. If we add together all of the bispectra, they will sum without cancelling. On the other hand, suppose that the phases of each of these frequencies was random. Then, the bispectrum will have the same magnitude (assuming that the magnitude of the frequency components is the same) but the phase will be randomly oriented. Adding together all of the bispectra will result in cancellation, because of the random phase orientation, and so the sum of the bispectra will have a small magnitude. Detecting phase coupling requires summation over a number of independent samples- this is the first motivation for defining the bicoherence. Secondly, the bispectrum is not normalized, because it still depends on the magnitudes of each of the frequency components. The bicoherence includes a normalization factor that removes the magnitude dependence.
| |
| | |
| There is some inconsistency with the definition of the bicoherence normalization constant. Some of the definitions that have been used are
| |
| :<math>b(f_1,f_2) = \frac{\left| \sum\limits_{n} F_n(f_1)F_n(f_2)F_n^*(f_1+f_2) \right|}{ \sqrt{\sum\limits_{n} |F_n(f_1)|^2|F_n(f_2)|^2|F_n^*(f_1+f_2)|^2}} </math>
| |
| which was provided in Sigl & Chamoun 1994, but does not appear to be correctly normalized. Alternatively, plasma physics typically uses
| |
| :<math>b^2(f_1,f_2) = \frac{\langle F_n(f_1)F_n(f_2)F_n^*(f_1+f_2) \rangle}{ \langle |F_n(f_1)F_n(f_2)|^2 \rangle\langle |F_n^*(f_1+f_2)|^2\rangle} </math> | |
| where the angle brackets denote averaging. Note that this is the same as using a sum, because <math>n</math> is the same in the numerator and the denominator. This definition is directly from Nagashima 2006, and is also referred to in He 2009 and Maccarone 2005.
| |
| | |
| Finally, one of the most intuitive definitions comes from Hagihira 2001 and Hayashi 2007, which is
| |
| | |
| :<math>b(f_1,f_2) = \frac{\left| \sum\limits_{n} F_n(f_1)F_n(f_2)F_n^*(f_1+f_2) \right|}{ \sum\limits_{n} |F_n(f_1)F_n(f_2)F_n^*(f_1+f_2)|} </math>
| |
| | |
| The numerator contains the magnitude of the bispectrum summed over all of the time series segments. This quantity is large if there is phase coupling, and approaches 0 in the limit of random phases. The denominator, which normalizes the bispectrum, is given by calculating the bispectrum after setting all of the phases to 0. This corresponds to the case where there is perfect phase coupling, because all of the samples have zero phase. Therefore, the bicoherence has a value between 0 (random phases) and 1 (total phase coupling).
| |
| | |
| ==References== | |
| * Hagihira, S., Takashina, M., Mori, T., Mashimo, T., & Yoshiya, I. (2001). Practical Issues in Bispectral Analysis of Electroencephalographic Signals. Anesthesia & Analgesia, 93(4), 966-970. Retrieved from http://www.anesthesia-analgesia.org/content/93/4/966.abstract
| |
| * Hayashi, K., Tsuda, N., Sawa, T., & Hagihira, S. (2007). Ketamine increases the frequency of electroencephalographic bicoherence peak on the alpha spindle area induced with propofol. British journal of anaesthesia, 99(3), 389-95. doi:10.1093/bja/aem175
| |
| * Nagashima, Y., Itoh, K., Itoh, S.-I., Hoshino, K., Fujisawa, A., Ejiri, A., Takase, Y., et al. (2006). Observation of coherent bicoherence and biphase in potential fluctuations around geodesic acoustic mode frequency on JFT-2M. Plasma Physics and Controlled Fusion, 48(5A), A377-A386. doi:10.1088/0741-3335/48/5A/S38
| |
| * He, H. (2009). The Canonical Bicoherence - Part I : Definition, Multitaper Estimation, and Statistics. Signal Processing, IEEE Transactions, 57(4), 1273-1284. Retrieved from http://ieeexplore.ieee.org/xpls/abs_all.jsp?arnumber=4749274
| |
| * Maccarone, T. J., & Schnittman, J. D. (2004). The bicoherence as a diagnostic for models of high-frequency quasi-periodic oscillations. Monthly Notices of the Royal Astronomical Society, 357(1), 12-16. doi:10.1111/j.1365-2966.2004.08615.x
| |
| * Mendel JM. "Tutorial on higher-order statistics (spectra) in signal processing and system theory: theoretical results and some applications." ''Proceedings of the IEEE'', '''79''', 3, 278-305
| |
| * M J Hinich, "Testing for Gaussianity and linearity of a stationary time series", ''Journal of Time Series Analysis'' '''3'''(3), 1982 pp 169–176.
| |
| *[http://www.mathworks.com/matlabcentral/fileexchange/loadFile.do?objectId=3013&objectType=file HOSA - Higher Order Spectral Analysis Toolbox]. ''(shareware for Microsoft Windows-type personal computers.)''
| |
| * T.H. Bullock, J.Z. Achimowicz ''et al.'', "Bicoherence of intracranial EEG in awake, sleep and seizures", Journal of Clinical Neurophysiology and EEG, 1997,vol.231,pp. 130–142.
| |
| | |
| [[Category:Complex analysis]]
| |
| [[Category:Integral transforms]]
| |
| [[Category:Fourier analysis]]
| |
Greetings. The author's title is Phebe and she feels comfortable when individuals use the complete visit the up coming article name. Hiring is her day job now and she will not alter it whenever quickly. Doing ceramics is what her family members and her enjoy. Puerto Rico is where he's always been residing but she std testing at home std over the counter std test test kit home requirements to move because of her family.
Feel free to surf to my homepage home std test kit [visit this website link]