主权项 |
一种非线性数据全局互相关性的并行定量计算方法,其特征在于,包括以下步骤:步骤1:对M通道原始脑电数据shuru进行直接分解,M>2,以两两通道数据最为原始数据进入下层循环;设定初始通道i=1,j=1;步骤2:取第i通道原始数据shurui;步骤3:取第j通道原始数据shuruj;步骤4:设置数据时间窗Epoch和滑动距离Overlap;步骤5:设定初始循环次数为t=1,计算循环次数N:N=length(shuru‑Epoch)/(Epoch‑Overlap)步骤4:按照时间窗从数据输入shurui和shuruj里取数据,计算当前窗口数据偏移量(t‑1)*(Epoch‑Overlap),得到当前数据X<sub>t</sub>(k)和Y<sub>t</sub>(k);步骤5:运用NLI算法计算两个数据的相关性,得到独立性测量S、H、N、M;步骤6:判断t>N是否成立;若是,则顺序执行下述步骤6;若否,则t=t+1,并回转执行所述步骤4;步骤7:判断j>M是否成立;若是,则顺序执行下述步骤8;若否,则j=j+1,并回转执行所述步骤3;步骤8:判断i>M是否成立;若是,则顺序执行下述步骤9;若否,则i=i+1,并回转执行所述步骤2;步骤9:对独立性测量S、H、N、M组成的相关性矩阵进行特征值特征向量分解;对数据矩阵X<sub>M×T</sub>={x<sub>i</sub>(k)},i=1,…,M,k=1,…,T的协方差矩阵C进行特征值分解,其中M表示通道数量,T表示时间窗内的数据点数,获得特征值λ<sub>i</sub>,并对其进行归一化,得到归一化特征值:<maths num="0001"><math><![CDATA[<mrow><msup><msub><mi>λ</mi><mi>i</mi></msub><mo>′</mo></msup><mo>=</mo><mfrac><msub><mi>λ</mi><mi>i</mi></msub><mrow><munderover><mo>Σ</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msub><mi>λ</mi><mi>i</mi></msub></mrow></mfrac><mo>=</mo><mfrac><msub><mi>λ</mi><mi>i</mi></msub><mrow><mi>t</mi><mi>r</mi><mrow><mo>(</mo><mi>C</mi><mo>)</mo></mrow></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000907429150000011.GIF" wi="424" he="211" /></maths>式中tr(C)为协方差矩阵的迹;步骤10:基于特征值特征向量计算S估计器,S估计器的定义为:<maths num="0002"><math><![CDATA[<mrow><mi>S</mi><mo>=</mo><mn>1</mn><mo>+</mo><mfrac><mrow><munderover><mo>Σ</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>M</mi></munderover><msup><msub><mi>λ</mi><mi>i</mi></msub><mo>′</mo></msup><mi>l</mi><mi>o</mi><mi>g</mi><mrow><mo>(</mo><msup><msub><mi>λ</mi><mi>i</mi></msub><mo>′</mo></msup><mo>)</mo></mrow></mrow><mrow><mi>log</mi><mrow><mo>(</mo><mi>M</mi><mo>)</mo></mrow></mrow></mfrac><mo>∈</mo><mo>[</mo><mn>0</mn><mo>,</mo><mn>1</mn><mo>]</mo><mo>.</mo></mrow>]]></math><img file="FDA0000907429150000021.GIF" wi="615" he="212" /></maths> |