[发明专利]一种脑电信号特征提取方法

权利要求书

1.一种脑电信号特征提取方法，该方法包括有以下步骤：

①小波重构

通过脑电图机采集的脑电信号，对脑电信号进行四级小波重构，第一级分解将采集的脑电信号从频率中点分开为低频、高频两部分，第二级分解将第一级分解得到的低频成分从其频率中点分开为低频、高频两部分，第三级分解将第二级分解得到的低频成分从其频率中点分开为低频、高频两部分，第四级分解将第三级分解得到的低频成分从其频率中点分开为低频、高频两部分，再通过第四级小波重构即可得到脑电信号的五种生理节律；

②加窗水平可视图复杂网络转换

脑电信号的采样频率为f，令时间窗时间长度为t₁,时间窗滑动距离为t₂，则时间窗内采样点个数为N_w=t₁×f，时间窗滑动的采样点个数为N_s=t₂×f，令时间窗内的脑电信号时间序列为{x_i}i=1,2,...,N_w，令i和j是复杂网络中的两个节点，如果对于所有的x_nn∈(i,j)，x_i,x_j>x_n，则认为i和j有链接，否则没有连接，从而建立N_w×N_w维复杂网络邻接矩阵M，其中，当i和j有链接时，M(i,j)=1；否则M(i,j)=0；

对步骤①中的脑电信号数据和所有小波重构获得的各节律脑电数据，或者根据需要选取通过小波重构获得的各节律脑电信号中的任意节律的脑电信号，按照上述方法尽心网络转换，得到脑电数据和各节律脑电数据对应的复杂网络邻接矩阵M；

③复杂网络参数特征提取；

对步骤②得到的复杂网络邻接矩阵M提取复杂网络的网络特性路径长度L和聚类系数C；

网络的平均路径长度定义为：

L=2N(N-1)Σi=2NΣj=1i-1dij---(1)]]>即网络中任意两节点的距离的平均值，其中N为网络节点个数，其中，d_ij定义为连接两个节点i和j的最短路径上的边数，网络的聚类系数定义为所有节点聚类系数的平均值：

C=1NΣi=1NCi---(2)]]>

其中N为网络节点个数，其中，C_i即为节点i的聚类系数，是当节点i的度为k_i时，k_i个节点间实际存在的边数R_i与最多存在的边数E_i=k_i(k_i-1)/2的比值；

将步骤②得到的脑电信号和各节律脑电信号对应的网络矩阵按照公式（1）和公式（2）计算，得到相应的网络特征参数：平均路径长度L和平均聚类系数C；

④通过步骤③，实时获得脑电信号及各节律脑电信号时间窗内数据的网络特征参数平均路径长度L和平均聚类系数C，平均路径长度L和平均聚类系数C作为脑电信号的特征值。