[发明专利]一种频谱匹配独立成分分析方法及系统

权利要求书

1.一种频谱匹配独立成分分析方法，其特征在于，包括：

S1、获取MEG系统所采集到的被试者数据中单个频带内的噪声功率B_a参数的具体值、单个频带内的源功率A_a参数的具体值、MEG系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵和MEG系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵M的具体值；

所述噪声功率B_a参数为MEG系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵在预先划分的频带中的第a个频带内的对应平均值；

所述源功率A_a参数为MEG系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵在预先划分的频带中的第a个频带内的对应平均值；

S2、基于所述MEG系统所采集到的被试者数据中单个频带内的噪声功率B_a参数的具体值、单个频带内的源功率A_a参数的具体值、MEG系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵和MEG系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵M的具体值，获取该MEG系统所采集到的被试者数据中预先设定的频带内的进行滤波的信号矩阵X_a中的源信号；

所述S2具体包括：

基于所述MEG系统所采集到的被试者数据中单个频带内的噪声功率B_a参数的具体值、单个频带内的源功率A_a参数的具体值、MEG系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵和MEG系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵M的具体值，根据公式（1）计算得到该MEG系统所采集到的被试者数据中第a频带范围上的进行滤波的信号矩阵X_a中的源信号；

所述公式（1）为：

其中，为该MEG系统所采集到的被试者数据中预先划分的频带中的第a个频带内的进行滤波的信号矩阵X_a对应恢复的源信号矩阵；

为MEG系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵和MEG系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵M的转置；

为所述MEG系统所采集到的被试者数据中单个频带内的噪声功率B_a参数的逆；

为单个频带内的源功率A_a参数的逆；

在S1之前还包括：

S0、基于预先获取的MEG系统所采集到的被试者数据、预先获取的背景噪音信号，建立频谱模型，所述频谱模型中包括噪声功率B_a参数、单个频带内的源功率A_a参数、MEG系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵和MEG系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵M参数；

相应的，所述S1具体包括：

基于预先获取的MEG系统所采集到的被试者数据、预先获取的MEG系统采集的空房间中的背景噪声信号和所述频谱模型，获取MEG系统所采集到的被试者数据中单个频带内的噪声功率B_a参数的具体值、单个频带内的源功率A_a参数的具体值、MEG系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵和MEG系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵M的具体值；

所述频谱模型为：

其中，C_a为C（f）在所划分的单个频带内的对应平均值；

C（f）为X（t）的自协方差矩阵的傅里叶变换；

其中，所述X（t）为p维向量，表示从具有p个探头传感器的MEG系统所采集到的被试者数据；

其中，

其中，N（t）为MEG系统所采集到的被试者数据时域内的背景噪声项且为p维向量；

表示MEG系统所采集到的被试者数据时域内的q个独立源且为q维向量；

M是一个的混合矩阵，表示假设通过该矩阵源信号与噪声信号实现加性混合。