如何比较假想的相干值?

信息处理 Python 脑电图 连贯性 阵列信号处理
2022-02-09 23:41:27

嗨,我目前正在写关于催眠易感性的硕士论文。我正在比较低易感个体和高易感个体的 EGG 信号。

想象的连贯性由下式给出:

ImCoh=Im{E{Sxy}}E{Sxx}E{Syy}

我发现我可以scipy.signal.coherence使用以下代码计算类似函数的连贯性:

sfarrA, Pxx = scipy.signal.welch(ex_sig1[0], fs=128)
sfarrB, Pyy = scipy.signal.welch(ex_sig1[9], fs=128)
sfarrAB, Pxy = scipy.signal.csd(ex_sig1[0],ex_sig1[9],fs=128) ```

coherence = abs(Pxy)**2/(Pxx*Pyy)

接下来,我使用以下代码计算虚构的相干性:

Icoh = np.imag(Pxy)/np.sqrt(Pxx*Pyy)

结果,我得到了一组负值和非负值。

array([ 0.        ,  0.15618352,  0.09958203,  0.00483643,  0.02767622,
        0.07590722,  0.06308123,  0.05571067,  0.0312004 , -0.06364436,
       -0.15178242, -0.0313257 ,  0.01023862, -0.02718701, -0.03659143,
        0.04386298, -0.05511231, -0.05728281, -0.03789811, -0.07106073,
       -0.34574961, -0.27614587, -0.27608529, -0.37483751, -0.34951999,

如何以有意义的方式比较这些值?

任何帮助将不胜感激。

非常感谢

0个回答
没有发现任何回复~