Imaginary part of mscohere

Question

Eric 2015-3-8

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/182120-imaginary-part-of-mscohere

评论： David Haslacher 2017-7-28

Hi everybody,

Is there a way to extract the imaginary part of the coherence when using mscohere? I can't find the original code anywhere. It seems to be protected by MATLAB.

I'm very interested in using the Welch method.

Best,

Eric

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Christiaan 2015-3-10

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/182120-imaginary-part-of-mscohere#answer_170714

编辑：Christiaan 2015-3-10

在 MATLAB Online 中打开

Dear Eric,

The formula for the Magniture Squared Coherence can be found here. Pxx and Pyy are the power spectral densities of the signals x and y. (PSD of Autocorrelation) These two quantities are real and positive. The magnitude of the cross power spectral density is (as any magnitude of complex numbers) also real and positive. Therefore you can from definition not extract a 'complex part of the MSC'.

If you still interested in the background of the mscohere function, you can type in the matlab prompt:

edit mscohere

Then you can find that the mscohere function uses the welch function. Since the function is private type in:

which welch -all

Now you have found the folder where the welch function is defined. The last step is to type in edit('.....link of previous output') and look at 'case 'mscohere'.

Good luck! Christiaan

Kind regards, Christiaan

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

David Haslacher 2017-7-28

I would wager that he's talking about coherency, Cab/sqrt(Caa * Cbb), which is complex [1]. Coherence is then the modulus of coherency.

Nolte G., Bai O., Wheaton L., Mari Z., Vorbach S., Hallett M. (2004). Identifying true brain interaction from EEG data using the imaginary part of coherency. Clin. Neurophysiol. 115, 2292–2307. 10.1016/j.clinph.2004.04.029 [PubMed] [Cross Ref]

请先登录，再进行评论。