[发明专利]脑电比功谱和比相干谱的测量分析方法及设备

权利要求书

1．一种脑电比功谱和比相干谱的测量分析方法，其特征是电极数量为特定的9个并分别安放在左脑F3、C3、P3左侧1cm处的F3＇、C3＇、P3＇，右脑F4、C4、P4右侧1cm处的F4＇、C4＇、P4＇和中线上的Fz、Cz、Pz位置：A／D转换器采用固定采样率为256／s；采样长度为64s或64s的2，3，4倍，采样程序分别为ad，ad2，ad3和ad4；分析程序为SP1，SP2，SP4和Scoh；脑电比功谱的计算方法采用“整数频率间隔内的功率积分法”，步骤为：

a．对于一个数字化的脑电信号序列，X(i)，i=1，2，…，N，先求出其快速哈特莱变换(FHT)系数，即DX(k)，k=1，…，N

b．求出1～70Hz内每个变换点上的脑电功率

PW(k)=DX(k)²+DX(N-k)²

在每个整数频率之间有等于所用脑电信号的样本长SPL(1，2或4秒)的变换点数故k=1，…，70×SPL

c．若样本长SPL大于1秒，则将每个整数频率间隔内的SPL个功率值加和，即 SPW ( k ) = Σ j = 1 SPL PW ( K + j ) , ]]>其中k=1，2，…，70(Hz)K=k+(k-1)SPL

d．将求得的积分功率值取对数变为以dB表示，即

P(k)=10log10(SPW(k))，k=1，2，…，70(Hz)

e．校正因SPL不同而引起的功率值偏移，即

Pc(k)=P(k)+CP(SPL)

其中，CP(1)=2．29，CP(2)=0．96，CP(4)=0．39

f．根据本发明装置的硬件特性将所得比功谱值转化为头皮记录处脑电的真实比功率值：设Vs为采样后数字化的脑电信号值，脑电放大器2的增益为G，考虑到模／数转换器3的精度为12位，即其转换增益为2．5×2×1000／2¹²，则折合到头皮记录部位的脑电信号幅值Vr为：

Vr=Vs×5000／4096×10³／G，以μν表示；

脑电比相干谱的计算方法采用“自互比功率相干法”，步骤如下：

a．对在不同脑区记录的二个数字化的脑电信号序列，X(i)和Y(i)，i=1，2，…N，先求出相应的快速哈特莱变换系数(即DX(k)和DY(k)，k=1，…，N)

b．根据“整数频率间隔内的功率积分法”的步骤a～c，计算出各自的比功率序列

SPWx(k)和SPWy(k)，k=1，2，…，70(Hz)

c．求出二者的比互功率

1)求出1～70Hz内各变换节点上的如下交叉系数：

TXYr(k)=DX(k)×DY(k)+DX(N-k)×DY(N-k)

TXYi(k)=DX(k)×DY(N-k)+DX(N-k)×DY(k)

其中k=1，2，…，70×SPL，SPL为脑电样本长(秒)

2)求出各整数频率间隔内上述TXYr(k)和TXYi(k)的SPL个值的加和，即 TXYR ( i ) = Σ j = 1 8 [ TXYr ( K = j ) ] ]]> TXYI ( i ) = Σ j = 1 8 [ TXYi ( K = j ) ] ]]>

其中，i=1，2，…，70(Hz)，K=i+(i-1)SPL

d．计算比相干幅谱

每个整数频率上的相干幅值为

CohA(i)=[TXYR(i)²+TXYI(i)²]／[SPWx(i)×SPWy(i)]

其中i=1，2，…，70(Hz)

e．计算比相干相谱

每个整数频率上的相干相位值为

CohP(i)=arctg[TXYI(i)／TXYR(i)]，i=1，2，…，70(Hz)

将所得的相干相位值的区间变换为-180°～+180°，即

若TXYR(i)＜0且TXYI(i)＞0，则Cohp(i)=CohP(i)+180

若TXYR(i)＜0且TXYI(i)＜0，则Cohp(i)=CohP(i)-180

若TXYR(i)＜0且TXYI(i)=0，则Cohp(i)=CohP(i)-180。