[发明专利]人体表面肌电信号单通道盲源分离法

权利要求书

1.一种人体表面肌电信号单通道盲源分离法，其特征是采用的模块包括有：多通道输入单通道输出的肌电信号预处理模块和单通道输入多通道输出的肌电信号单通道盲源分离模块，实施步骤是：

A、将采集到多路肌电信号通过屏蔽导线送到肌电信号预处理模块直接相加，得到预处理单通道肌电信号                                                ；

B、将预处理所得到的单通道肌电信号送到肌电信号盲源分离模块，进行总体经验模态分解（Ensemble Empirical Mode Decomposition），即EEMD分解，主成份分析PCA降维和ICA分析，实现多路表面肌电信号通过一个输入口采集，多路输出口输出：

B.1、总体经验模态分解EEMD处理得到本征模态函数IMFs： 

a、将预处理所得到的单通道肌电信号,多次加入具有均值为零、标准差为常数的白噪声,即为,其中为加入白噪声之后的信号，为第i次加入的白噪声，白噪声n_i(t)的准则为，，，其中，表示加入的高斯白噪声标准差，表示肌电信号中有效高频成分的幅值标准差，表示肌电信号幅值标准差，为比例系数，通常情况下，=能有效避免肌电信号分解中的模式混淆；

b、对所得到的信号进行经验模态分解EMD，得到各自的IMF记为,和余项，其中表示加入白噪声后分解得到的第个IMF；

c、将步骤b所得到的IMF进行总体平均运算，得到IMFs为，其中，为对原信号进行EEMD分解后所得到的第个IMF，由此可以得到IMF分量矩阵,其中，为单通道肌电信号经EEMD分解后的IMF个数，上标T为转置运算；

B.2、 对得到IMF分量进行PCA降维：

，其中A为经过EEMD分解后得到的IMF分量矩阵，R为m个变量IMFs的自相关矩阵，V为R的阶特征向量矩阵，其列向量是R的正交规一化的特征向量，为R的特征对角矩阵, ，i=1,2, …,m为第i个对角线上的元素，构造m个不相关的新变量，Y = {y₁,y₂,…,y_m}^T,对（i=1,2, …,m）按降序排列后，取前面p个较大特征值所对应的特征向量，得到降维后阶向量矩阵B，其中；

B.3、将PCA降维所得到矩阵B进行ICA处理，采用的是效果好且速度快的基本定点迭代算法FastICA进行ICA处理：

a、白化数据，给出观测矩阵X，其中，A为信号的混合矩阵，为源信号；

b、随机选取权向量作为解混矩阵W的权值，W为解混矩阵，即为混合矩阵A的虚拟反矩阵；

c、令，其中非二次函数G的导数，，是常数，常取做1；

d、循环迭代，直到收敛为止，最后得到阶向量矩阵，通过观察选取分离后肌电信号。