一种频谱匹配独立成分分析方法及系统与流程

1.本发明涉及脑磁图医疗设备数据处理技术领域，尤其涉及一种频谱匹配独立成分分析方法及系统。


背景技术：

2.由于脑磁图具有高时间分辨率与空间分辨率特性，其在脑科学研究和大脑疾病诊断等方面具有广泛应用前景。然而，现有的基于超导量子干涉装置(squid)的脑磁图系统需要在低温下运行，因此必须固定在低温杜瓦瓶内，使得系统变得庞大、笨重且它们无法适应不同的头部形状或尺寸，为脑磁图系统的信号测量精度和使用场景带来了巨大的挑战。近年来量子技术的进步促进了新型磁场传感器的发展，基于量子传感技术的光泵磁强计(opm)可提供亚ft级极弱磁磁场测量，其灵敏度与传统meg中使用的squid传感器相似，但它们不需要冷却，可贴近头皮进行测量，而且可以集成芯片化,其体积轻便，因此基于opm（optical pump magnetometer：光泵磁强计）的脑磁图系统被广泛的研究。
3.目前基于opm的可穿戴meg（magnetoencephalogram：脑磁图）系统已经被研制出来，其具备能够提供较高的数据质量、适应任何受试者、允许在移动过程中捕获数据等特点，已经被许多研究机构引入到脑科学研究领域进行信号测量。但是，应用基于opm的可穿戴meg系统进行信号测量分析，关键的步骤是进行源分离，因为测到的信号中不仅包含了大脑活动信号源，还会包含眨眼、心跳、线噪声、肌肉活动等干扰源，这些干扰源的存在会对源定位和分析产生较大影响，导致分析结果不准确甚至出错。
4.在脑磁图信号分析领域进行源分离较为流行的方法是独立成分分析，但现存的独立成分分析方法大多采用非高斯独立源的无噪声混合模型，其存在强制估计与传感器一样多的源，依赖于损失信号的pca降维操作等限制条件，导致其对噪声变化较大的移动情况下捕捉到的数据进行源分离准确性较差，使得基于opm的脑磁图不能较好的应用于人体移动情况下的数据测量与分析。


技术实现要素：

5.（一）要解决的技术问题鉴于现有技术的上述缺点、不足，本发明提供一种频谱匹配独立成分分析方法及系统，其解决了现有的独立成分分析方法大多采用非高斯独立源的无噪声混合模型，其存在强制估计与传感器一样多的源，依赖于损失信号的pca降维操作等限制条件，导致其对噪声变化较大的移动情况下捕捉到的数据进行源分离准确性较差技术问题。
6.（二）技术方案为了达到上述目的，本发明采用的主要技术方案包括：第一方面，本发明实施例提供一种频谱匹配独立成分分析方法，包括：s1、获取meg系统所采集到的被试者数据中单个频带内的噪声功率ba参数的具体值、单个频带内的源功率aa参数的具体值、meg系统所采集到的被试者数据频域内的噪声项
的频谱协方差矩阵和meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵m的具体值；所述噪声功率ba参数为meg系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵在预先划分的频带中的第a个频带内的对应平均值；所述源功率aa参数为meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵在预先划分的频带中的第a个频带内的对应平均值；s2、基于所述meg系统所采集到的被试者数据中单个频带内的噪声功率ba参数的具体值、单个频带内的源功率aa参数的具体值、meg系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵和meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵m的具体值，获取该meg系统所采集到的被试者数据中预先设定的频带内的进行滤波的信号矩阵xa中的源信号。
7.优选的，所述s2具体包括：基于所述meg系统所采集到的被试者数据中单个频带内的噪声功率ba参数的具体值、单个频带内的源功率aa参数的具体值、meg系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵和meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵m的具体值，根据公式（1）计算得到该meg系统所采集到的被试者数据中第a频带范围上的进行滤波的信号矩阵xa中的源信号；所述公式（1）为：其中，为该meg系统所采集到的被试者数据中预先划分的频带中的第a个频带内的进行滤波的信号矩阵xa对应恢复的源信号矩阵；为meg系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵和meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵m的转置；为所述meg系统所采集到的被试者数据中单个频带内的噪声功率ba参数的逆；为单个频带内的源功率aa参数的逆。
8.优选的，在s1之前还包括：s0、基于预先获取的meg系统所采集到的被试者数据、预先获取的背景噪音信号，建立频谱模型，所述频谱模型中包括噪声功率ba参数、单个频带内的源功率aa参数、meg系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵和meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵m参数；相应的，所述s1具体包括：基于预先获取的meg系统所采集到的被试者数据、预先获取的meg系统采集的空房间中的背景噪声信号和所述频谱模型，获取meg系统所采集到的被试者数据中单个频带内的噪声功率ba参数的具体值、单个频带内的源功率aa参数的具体值、meg系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵和meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵m的具体值。
9.优选的，
所述频谱模型为：其中，ca为c（f）在所划分的单个频带内的对应平均值；c（f）为x（t）的自协方差矩阵的傅里叶变换；其中，所述x（t）为p维向量，表示从具有p个探头传感器的meg系统所采集到的被试者数据；其中，其中，n（t）为meg系统所采集到的被试者数据时域内的背景噪声项且为p维向量；表示meg系统所采集到的被试者数据时域内的q个独立源且为q维向量；m是一个的混合矩阵，表示假设通过该矩阵源信号与噪声信号实现加性混合。
10.优选的，所述s1具体包括：s11、基于预先获取的meg系统所采集到的被试者数据，获取ca参数的具体值；s12、基于预先获取的背景噪声信号，获取ba参数的具体值；s13、基于所述ca参数的具体值和所述ba参数的具体值，获取单个频带内的源功率aa参数的具体值、meg系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵和meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵m的具体值。
11.优选的，所述s11具体包括：基于预先获取的meg系统所采集到的被试者数据，采用公式（2）获取预先获取的meg系统所采集到的被试者数据x（t）的离散样本x（r）的傅里叶系数；所述公式（2）为：；为meg系统所采集到的被试者数据x（t）的离散样本x（r）的傅里叶系数；r表示预先获取的meg系统所采集到的被试者数据时域采样点数；r表示离散样本x（r）中的第r个样本；基于所述x（r）的傅里叶系数，采用公式（3）获取的具体值，并将所述的具体值作为ca的具体值；所述公式（3）为：；为ca的估计值；h表示对矩阵做共轭转置；为的共轭转置；其中，na为预先划分meg系统所采集到的被试者数据的频带中的第a个频带内可用的傅里叶系数数目，并且，；
其中，#表示计数；ia为预先划分的频带中的第a个频带；k表示总的傅里叶系数数目。优选的，所述s12具体包括：基于预先获取的meg系统所采集到空房间中的背景噪声信号，采用公式（4）获取预先获取的meg系统所采集到的背景噪声信号n（t）的离散样本n（w）的傅里叶系数；所述公式（4）为：；为meg系统所采集到的背景噪声信号n（t）的离散样本n（w）的傅里叶系数；w表示meg系统所采集到的背景噪声信号时域采样点数；w表示离散样本n（w）中的第w个样本；基于所述n（w）的傅里叶系数，采用公式（5）获取的具体值，并将所述的具体值作为ba的具体值；所述公式（5）为：；为ba的估计值；为的共轭转置；其中，za为预先划分meg系统所采集到的背景噪声信号的频带中的第a个频带内可用的傅里叶系数数目，并且，。
12.优选的，所述s13具体包括：基于kl散度构造频谱匹配标准，所述频谱匹配标准的形式为：其中，m为频带数目；表示和的kl散度；基于所述频谱匹配标准、所述ca参数的具体值和所述ba参数的具体值，采用期望最大化算法获取单个频带内的源功率aa参数的具体值和meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵m的具体值。
13.优选的，所述meg系统所采集到的被试者数据为用户处于背景剩磁消减后的磁屏蔽房内并在相同声音信号的刺激下静止坐下、静止站立以及身体和头部移动时所采集的meg被试者
数据；所述背景剩磁消减后的磁屏蔽房内的背景磁场低于50nt；预先获取的背景噪音信号为背景剩磁消减后的磁屏蔽房内的预先设定时间段内的空的磁屏蔽房内的磁信号。
14.第二方面，本发明实施例提供一种频谱匹配独立成分分析系统，所述系统包括：至少一个处理器；以及与所述处理器通信连接的至少一个存储器，其中，所述存储器存储有可被所述处理器执行的程序指令，所述处理器调用所述程序指令能够执行上述任一的频谱匹配独立成分分析方法。
15.（三）有益效果本发明的有益效果是：本发明的一种频谱匹配独立成分分析方法及系统，由于通过获取的meg系统所采集到的被试者数据中单个频带内的噪声功率ba参数的具体值、单个频带内的源功率参数aa的具体值、meg系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵和meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵m的具体值，得到该meg系统所采集到的被试者数据中预先设定的频带内的进行滤波的信号矩阵xa中的源信号，使得源信号恢复更为准确。
附图说明
16.图1为本发明的一种频谱匹配独立成分分析方法流程图；图2为本发明实施例二中的meg系统中的磁屏蔽系统结构示意图；图3为本发明实施例二中的meg系统中的opm-meg头盔阵列数据获取系统结构示意图；图4为本发明实施例二中的meg系统中的听力刺激系统结构示意图；图5为本发明实施例二中的meg系统在采集被试者数据时的示意图；图6为本发明实施例二中的频谱匹配独立成分分析方法流程图。
17.【附图标记说明】101：磁屏蔽房；102：磁场消零线圈；103：参考阵列；104：线圈驱动模块；105：ni daq；1051：dac；1052：adc；106：传感器采集控制电路；107：传感器、线圈控制设备；201：头盔阵列；203：数据获取设备；301：无磁性塑料细管耳机；302：听力刺激设备。
具体实施方式
18.为了更好的解释本发明，以便于理解，下面结合附图，通过具体实施方式，对本发明作详细描述。
19.为了更好的理解上述技术方案，下面将参照附图更详细地描述本发明的示例性实施例。虽然附图中显示了本发明的示例性实施例，然而应当理解，可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反，提供这些实施例是为了能够更清楚、透彻地理解本发明，并且能够将本发明的范围完整的传达给本领域的技术人员。
20.实施例一参见图1，本实施例提供一种频谱匹配独立成分分析方法，包括：s1、获取meg系统所采集到的被试者数据中单个频带内的噪声功率ba参数的具体值、单个频带内的源功率aa参数的具体值、meg系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵和meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵m的具体值。
21.所述噪声功率ba参数为meg系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵在预先划分的频带中的第a个频带内的对应平均值。
22.所述源功率aa参数为meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵在预先划分的频带中的第a个频带内的对应平均值。
23.具体的，本实施例中meg系统所采集到的被试者数据的频率范围为0-70hz，然后被预先划分为40个频带，每个频带的频率间隔范围都相同。
24.s2、基于所述meg系统所采集到的被试者数据中单个频带内的噪声功率ba参数的具体值、单个频带内的源功率aa参数的具体值、meg系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵和meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵m的具体值，获取该meg系统所采集到的被试者数据中预先设定的频带内的进行滤波的信号矩阵xa中的源信号。
25.在本实施例的实际应用中，所述s2具体包括：基于所述meg系统所采集到的被试者数据中单个频带内的噪声功率ba参数的具体值、单个频带内的源功率aa参数的具体值、meg系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵和meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵m的具体值，根据公式（1）计算得到该meg系统所采集到的被试者数据中第a频带范围上的进行滤波的信号矩阵xa中的源信号。
26.所述公式（1）为：其中，为该meg系统所采集到的被试者数据中预先划分的频带中的第a个频带内的进行滤波的信号矩阵xa对应恢复的源信号矩阵。
27.为meg系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵和meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵m的转置。
28.为所述meg系统所采集到的被试者数据中单个频带内的噪声功率ba参数的逆。
29.为单个频带内的源功率aa参数的逆。
30.在本实施例的实际应用中，在s1之前还包括：s0、基于预先获取的meg系统所采集到的被试者数据、预先获取的背景噪音信号，建立频谱模型，所述频谱模型中包括噪声功率ba参数、单个频带内的源功率aa参数、meg系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵和meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵m参数。
31.相应的，所述s1具体包括：基于预先获取的meg系统所采集到的被试者数据、预先获取的meg系统采集的空房间中的背景噪声信号和所述频谱模型，获取meg系统所采集到的被试者数据中单个频带内的噪声功率ba参数的具体值、单个频带内的源功率aa参数的具体值、meg系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵和meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵m的具体值。
32.在本实施例的实际应用中，所述频谱模型为：其中，ca为c（f）在所划分的单个频带内的对应平均值；c（f）为x（t）的自协方差矩阵的傅里叶变换；其中，所述x（t）为p维向量，表示从具有p个探头传感器的meg系统所采集到的被试者数据。
33.其中，其中，n（t）为meg系统所采集到的被试者数据时域内的背景噪声项且为p维向量；表示meg系统所采集到的被试者数据时域内的q个独立源且为q维向量；m是一个的混合矩阵，表示假设通过该矩阵源信号与噪声信号实现加性混合。
34.在本实施例的实际应用中，所述s1具体包括：s11、基于预先获取的meg系统所采集到的被试者数据，获取ca参数的具体值。
35.s12、基于预先获取的背景噪声信号，获取ba参数的具体值。
36.s13、基于所述ca参数的具体值和所述ba参数的具体值，获取单个频带内的源功率aa参数的具体值、meg系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵和meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵m的具体值。
37.在本实施例的实际应用中，所述s11具体包括：基于预先获取的meg系统所采集到的被试者数据，采用公式（2）获取预先获取的meg系统所采集到的被试者数据x（t）的离散样本x（r）的傅里叶系数。
38.所述公式（2）为：；为meg系统所采集到的被试者数据x（t）的离散样本x（r）的傅里叶系数。r表示预先获取的meg系统所采集到的被试者数据时域采样点数；r表示离散样本x（r）中的第r个样本。
39.基于所述x（r）的傅里叶系数，采用公式（3）获取的具体值，并将所述的具体值作为ca的具体值。
40.所述公式（3）为：
。
41.为ca的估计值。
42.h表示对矩阵做共轭转置；为的共轭转置。
43.其中，na为预先划分meg系统所采集到的被试者数据的频带中的第a个频带内可用的傅里叶系数数目，并且，；其中，#表示计数。
44.ia为预先划分的频带中的第a个频带。
45.k表示总的傅里叶系数数目。
46.在本实施例的实际应用中，所述s12具体包括：基于预先获取的meg系统所采集到空房间中的背景噪声信号，采用公式（4）获取预先获取的meg系统所采集到的背景噪声信号n（t）的离散样本n（w）的傅里叶系数。
47.所述公式（4）为：；为meg系统所采集到的背景噪声信号n（t）的离散样本n（w）的傅里叶系数。
48.w表示meg系统所采集到的背景噪声信号时域采样点数。
49.w表示离散样本n（w）中的第w个样本。
50.基于所述n（w）的傅里叶系数，采用公式（5）获取的具体值，并将所述的具体值作为ba的具体值。
51.所述公式（5）为：。
52.为ba的估计值。
53.为的共轭转置。
54.其中，za为预先划分meg系统所采集到的背景噪声信号的频带中的第a个频带内可用的傅里叶系数数目，并且，。
55.在本实施例的实际应用中，所述s13具体包括：基于kl散度构造频谱匹配标准，所述频谱匹配标准的形式为：
其中，m为频带数目。
56.表示和的kl散度。
57.基于所述频谱匹配标准、所述ca参数的具体值和所述ba参数的具体值，采用期望最大化算法获取单个频带内的源功率aa参数的具体值和meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵m的具体值。
58.在本实施例的实际应用中，所述meg系统所采集到的被试者数据为用户处于背景剩磁消减后的磁屏蔽房内并在相同声音信号的刺激下静止坐下、静止站立以及身体和头部移动时所采集的meg被试者数据。
59.所述背景剩磁消减后的磁屏蔽房内的背景磁场低于50nt。
60.预先获取的背景噪音信号为背景剩磁消减后的磁屏蔽房内的预先设定时间段内的空的磁屏蔽房内的磁信号。
61.另一方面，本实施例还提供一种频谱匹配独立成分分析系统，所述系统包括：至少一个处理器；以及与所述处理器通信连接的至少一个存储器，其中，所述存储器存储有可被所述处理器执行的程序指令，所述处理器调用所述程序指令能够执行上述频谱匹配独立成分分析方法。
62.由于本发明上述实施例所描述的系统，为实施本发明上述实施例的方法所采用的系统，故而基于本发明上述实施例所描述的方法，本领域所属技术人员能够了解该系统的具体结构及变形，因而在此不再赘述。凡是本发明上述实施例的方法所采用的系统都属于本发明所欲保护的范围。
63.本实施例中的一种频谱匹配独立成分分析方法及系统，由于通过获取的meg系统所采集到的被试者数据中单个频带内的噪声功率ba参数的具体值、单个频带内的源功率参数aa的具体值、meg系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵和meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵m的具体值，得到该meg系统所采集到的被试者数据中预先设定的频带内的进行滤波的信号矩阵xa中的源信号，使得源信号恢复更为准确。
64.实施例二参见图2-5，本实施例中的meg系统包括：磁屏蔽系统、opm-meg头盔阵列数据获取系统、听力刺激系统。
65.其中，本实施例中磁屏蔽系统包括磁屏蔽房101、磁场消零线圈102、由xyz方向构成的参考阵列103、ni daq（ni数据采集系统）105、传感器采集控制电路106、传感器和线圈控制设备107、线圈驱动模块104等构件。
66.其中，ni daq105包括：dac（数模转换器）1051，adc（模数转换器）1052。
67.opm-meg头盔阵列数据获取系统包括头盔阵列201（由3d打印头盔和一套32通道opm探头组成）、ni daq105、数据获取设备203等构件。
68.听力刺激系统包括听力刺激设备302、无磁性塑料细管耳机301等构件。
69.meg系统采集信号前，会启动磁屏蔽房进行主动消磁，在完成主动消磁后打开opm
传感器进行参考信号测量，利用测量的参考信号和剩磁消零算法进行背景剩磁消减，以实现最大程度的背景磁场抑制。
70.开始进行空房间的磁信号测量，记录时长与后面被试者进入房间到完成实验的计划时长一致（也就是预先设定时间段），将其视为背景噪声信号。
71.然后 ，受试者进入房间，背上提前准备好的双肩背包，戴上3d打印的无磁性头盔，并将头盔线固定在受试者所背的双肩背包上用夹子夹上。根据实验需要将32个探头分别插入对应插槽，上下两层排列阵列完全相同，最后将塑料细管耳机戴在耳朵上。
72.最后关上磁屏蔽房门，再重复上述剩磁补偿步骤，至此所述背景剩磁消减后的磁屏蔽房内的背景磁场低于50nt，完成数据采集前所有准备。
73.发出声音信号，指示受试者按声音信号要求做出相应动作，并播放声音信号，声音信号由高保真声卡通过无磁性塑料细管耳机呈现给受试者。同样的声音信号的刺激作用下，受试者身体状况在静止坐下、静止站立、身体和头部移动情况下进行数据采集，得到meg系统所采集到的被试者数据。
74.本实施例中通过测量空房间下的信号作为背景噪声信号，使得模型解算更为接近实际情况，源分离准确度提高。
75.本实施例在声音刺激下采用了静态坐立、静态站立、身体和大脑移动等三种状态进行数据采集，增强了数据的可参照性。
76.本实施例中meg系统的数据采集能够实现较高质量自然状态的数据采集，使得测量数据具有可对比性和参照性。
77.参见图6，本实施例中的频谱匹配独立成分分析方法包括：q0、基于预先获取的meg系统所采集到的被试者数据、预先获取的背景噪音信号，建立频谱模型，所述频谱模型中包括噪声功率ba参数、单个频带内的源功率aa参数、meg系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵和meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵m参数。
78.所述频谱模型为：其中，ca为c（f）在所划分的单个频带内的对应平均值。
79.c（f）为x（t）的自协方差矩阵的傅里叶变换。
80.其中，所述x（t）为p维向量，表示从具有p个探头传感器的meg系统所采集到的被试者数据。
81.其中，其中，n（t）为meg系统所采集到的被试者数据时域内的背景噪声项且为p维向量。
82.表示meg系统所采集到的被试者数据时域内的q个独立源且为q维向量。
83.m是一个的混合矩阵，表示假设通过该矩阵源信号与噪声信号实现加性混合。
84.本实施例中的频谱模型能够将预先获取的meg系统所采集到的被试者数据降维嵌套入模型，克服了传统独立成分分析需要假设源与噪声相同可能导致的模型退化问题。频谱模型构造假设源和噪声信号为平稳的高斯时间序列，考虑到了源和噪声信号的时间相关性，使得源信号恢复更准确。
85.q1、基于预先获取的meg系统所采集到的被试者数据、预先获取的meg系统采集的
空房间中的背景噪声信号和所述频谱模型，获取meg系统所采集到的被试者数据中单个频带内的噪声功率ba参数的具体值、单个频带内的源功率aa参数的具体值、meg系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵和meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵m的具体值。
86.所述噪声功率ba参数为meg系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵在预先划分的频带中的第a个频带内的对应平均值。
87.所述源功率aa参数为meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵在预先划分的频带中的第a个频带内的对应平均值。
88.q2、基于所述meg系统所采集到的被试者数据中单个频带内的噪声功率ba参数的具体值、单个频带内的源功率aa参数的具体值、meg系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵和meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵m的具体值，根据公式（1）计算得到该meg系统所采集到的被试者数据中第a频带范围上的进行滤波的信号矩阵xa中的源信号。
89.所述公式（1）为：其中，为该meg系统所采集到的被试者数据中预先划分的频带中的第a个频带内的进行滤波的信号矩阵xa对应恢复的源信号矩阵。
90.为meg系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵和meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵m的转置。
91.为所述meg系统所采集到的被试者数据中单个频带内的噪声功率ba参数的逆。
92.为单个频带内的源功率aa参数的逆。
93.具体的，在本实施例中，所述q1具体包括：q11、基于预先获取的meg系统所采集到的被试者数据，获取ca参数的具体值。
94.q12、基于预先获取的背景噪声信号，获取ba参数的具体值。
95.q13、基于所述ca参数的具体值和所述ba参数的具体值，获取单个频带内的源功率aa参数的具体值、meg系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵和meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵m的具体值。
96.具体的，在本实施例中，所述q11具体包括：基于预先获取的meg系统所采集到的被试者数据，采用公式（2）获取预先获取的meg系统所采集到的被试者数据x（t）的离散样本x（r）的傅里叶系数。
97.所述公式（2）为：；为meg系统所采集到的被试者数据x（t）的离散样本x（r）的傅里叶系数。
98.r表示预先获取的meg系统所采集到的被试者数据时域采样点数；r表示离散样本x（r）中的第r个样本。
99.基于所述x（r）的傅里叶系数，采用公式（3）获取的具体值，并将所述的具体值作为ca的具体值。
100.所述公式（3）为：；为ca的估计值。
101.h表示对矩阵做共轭转置；为的共轭转置。
102.其中，na为预先划分meg系统所采集到的被试者数据的频带中的第a个频带内可用的傅里叶系数数目，并且，；其中，#表示计数。
103.ia为预先划分的频带中的第a个频带。
104.k表示总的傅里叶系数数目。具体的，在本实施例中，所述q12具体包括：基于预先获取的meg系统所采集到空房间中的背景噪声信号，采用公式（4）获取预先获取的meg系统所采集到的背景噪声信号n（t）的离散样本n（w）的傅里叶系数。
105.所述公式（4）为：；为meg系统所采集到的背景噪声信号n（t）的离散样本n（w）的傅里叶系数。
106.w表示meg系统所采集到的背景噪声信号时域采样点数。
107.w表示离散样本n（w）中的第w个样本。
108.基于所述n（w）的傅里叶系数，采用公式（5）获取的具体值，并将所述的具体值作为ba的具体值。
109.所述公式（5）为：；为ba的估计值。
110.为的共轭转置。
111.其中，za为预先划分meg系统所采集到的背景噪声信号的频带中的第a个频带内可用的傅里叶系数数目，并且，
。
112.具体的，在本实施例中，所述q13具体包括：基于kl散度构造频谱匹配标准，所述频谱匹配标准的形式为：其中，m为频带数目。
113.表示和的kl散度。
114.基于所述频谱匹配标准、所述ca参数的具体值和所述ba参数的具体值，采用期望最大化算法获取单个频带内的源功率aa参数的具体值和meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵m的具体值。
115.本实施例中的一种频谱匹配独立成分分析方法，由于通过获取的meg系统所采集到的被试者数据中单个频带内的噪声功率ba参数的具体值、单个频带内的源功率aa参数的具体值、meg系统所采集到的被试者数据频域内的噪声项的频谱协方差矩阵和meg系统所采集到的被试者数据频域内q个独立源频谱协方差矩阵的加性混合矩阵m的具体值，得到该meg系统所采集到的被试者数据中预先设定的频带内的进行滤波的信号矩阵xa中的源信号，使得源信号恢复更为准确。
116.本领域内的技术人员应明白，本发明的实施例可提供为方法、系统或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例，或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质（包括但不限于磁盘存储器、cd-rom、光学存储器等）上实施的计算机程序产品的形式。
117.本发明是参照根据本发明实施例的方法、设备（系统）和计算机程序产品的流程图和／或方框图来描述的。应理解可由计算机程序指令实现流程图和／或方框图中的每一流程和／或方框、以及流程图和／或方框图中的流程和／或方框的结合。
118.应当注意的是，在权利要求中，不应将位于括号之间的任何附图标记理解成对权利要求的限制。词语“包含”不排除存在未列在权利要求中的部件或步骤。位于部件之前的词语“一”或“一个”不排除存在多个这样的部件。本发明可以借助于包括有若干不同部件的硬件以及借助于适当编程的计算机来实现。在列举了若干装置的权利要求中，这些装置中的若干个可以是通过同一个硬件来具体体现。词语第一、第二、第三等的使用，仅是为了表述方便，而不表示任何顺序。可将这些词语理解为部件名称的一部分。
119.此外，需要说明的是，在本说明书的描述中，术语“一个实施例”、“一些实施例”、“实施例”、“示例”、“具体示例”或“一些示例”等的描述，是指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
120.尽管已描述了本发明的优选实施例，但本领域的技术人员在得知了基本创造性概念后，则可对这些实施例作出另外的变更和修改。所以，权利要求应该解释为包括优选实施例以及落入本发明范围的所有变更和修改。
121.显然，本领域的技术人员可以对本发明进行各种修改和变型而不脱离本发明的精神和范围。这样，倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也应该包含这些修改和变型在内。