阅读说明：本技术 基于多导联脑电最优导联选择的抑郁检测系统 (Depression detection system based on optimal lead selection of multi-lead electroencephalogram ) 是由 胡斌 张晓炜 沈健 李娜 潘静 李俊磊 吴漫兮 于 2019-11-22 设计创作，主要内容包括：本发明提供的一种基于多导联脑电最优导联选择的抑郁检测系统,有效去除多导联脑电信号引起的信息冗余,降低计算复杂度,提高抑郁检测的识别率。包括：多导联脑电信号获取模块、脑电信号预处理模块、特征提取模块、导联选择模块、抑郁检测模块；多导联脑电信号获取模块用于获取静息状态采集的多导联脑电信号,脑电信号预处理模块用于对获取的多导联脑电信号进行预处理,特征提取模块用于提取进行最优导联选择的有效特征,所述导联选择模块用于采用修改版本的核排列作为目标函数,利用新的二元粒子群优化法来对目标函数进行导联选择,求取最优化导联子集,抑郁检测模块用于对选择出来的最优化导联子集进行识别,构建最优化抑郁检测分类模型,来进行抑郁检测。(The depression detection system based on the optimal lead selection of the multi-lead electroencephalogram effectively removes information redundancy caused by the multi-lead electroencephalogram signals, reduces the calculation complexity and improves the identification rate of depression detection. The method comprises the following steps: the device comprises a multi-lead electroencephalogram signal acquisition module, an electroencephalogram signal preprocessing module, a feature extraction module, a lead selection module and a depression detection module; the device comprises a multi-lead electroencephalogram signal acquisition module, an electroencephalogram signal preprocessing module, a characteristic extraction module, a lead selection module, a new binary particle swarm optimization method and a depression detection module, wherein the multi-lead electroencephalogram signal acquisition module is used for acquiring multi-lead electroencephalogram signals acquired in a resting state, the electroencephalogram signal preprocessing module is used for preprocessing the acquired multi-lead electroencephalogram signals, the characteristic extraction module is used for extracting effective characteristics for optimal lead selection, the lead selection module is used for adopting modified version nuclear arrangement as a target function, conducting lead selection on the target function by utilizing the new binary particle swarm optimization method, obtaining an optimized lead subset, and the depression detection module is used for identifying the selected optimized lead subset, constructing an optimized depression.)

基于多导联脑电最优导联选择的抑郁检测系统

技术领域

本发明涉及基于多导联脑电最优导联选择的抑郁检测系统，属于抑郁症辅助医疗技术领域。

背景技术

抑郁症是一种可以持续数周、数月甚至数年的主要的精神疾病，它会对个体的生活能力造成损害。根据世界卫生组织预测，抑郁症将在2020年成为威胁人类生命的第二大疾病。抑郁症可以显著地影响人类的思维、行为、感觉和幸福感。此外，根据世界卫生组织预测，全世界有3.5亿人受抑郁症影响，并且每年有近100万抑郁患者***。抑郁症已经严重威胁到了人类的生命。虽然抑郁症可以通过一些有效的临床手段进行治疗，但是由于忽视、诊断不及时以及误诊，全世界有大量的抑郁症患者正在忍受抑郁症带来的痛苦。因此，尽早地进行抑郁症检测对于治愈抑郁症从而改善抑郁症患者的生活质量、保障其身心健康是至关重要的，同时也可以直接降低抑郁症带来的社会和经济压力。

目前，诊断和检测抑郁症的主要临床手段是患者自评量表和临床医生访谈，但是患者自评量表通常带有主观性，临床医生访谈往往受限于他们的临床技能和诊断方法。此外，抑郁诊断和检测的临床实践中没有客观的评估标准。随着抑郁症患者数量的增长，通过早期的检测和评估来对抑郁症患者进行有效的跟踪治疗已经刻不容缓。为了获得更为客观的检测结果，国内外研究者们通过表情、语音以及姿势等行为信号来进行抑郁检测，运用多种机器学习、深度学习方法已经取得了一些成果，但是由于行为信号是间接信号，具有易控性和可伪装性，信号本身可能会受到被试的主观因素的干扰，有时可能无法反映被试潜在的、真实的精神状态。脑电信号可以直接反映大脑内部的工作状态，它是自发的不受人为控制的大脑内部皮层放电活动，具有很高的时间分辨率。而且，脑电具有采集安全、无创、成本低廉以及采集简便等特性。因此，使用脑电信号来提供客观的抑郁检测是可行的并具有前景的，基于脑电的抑郁检测研究为临床诊断和治疗抑郁症提供了更加便捷的新途径。近年来，基于脑电的抑郁检测研究也受到了专家和学者们的广泛关注。

当前基于脑电信号的抑郁检测研究，识别率尚待提高。并且，由于临床采集的脑电信号往往具有多导联、高采样率以及高密度等特性，临床脑电信号的信息冗余和计算复杂度大大提高，给机器学习和深度学习技术带来了巨大的挑战。大量的导联可能会造成机器学习和深度学习模型过拟合，甚至，不相关的导联可能会造成抑郁症检测识别率的损失。此外，在临床采集脑电信号时，设置大量的导联会耗费大量的时间，同时也会引起被试不适，从而导致采集到的被试脑电信号不能准确反映被试的精神状态。

发明内容

本发明针对现有技术存在的问题，提供一种基于多导联脑电最优导联选择的抑郁检测系统，针对脑电信号在抑郁症临床检测中存在的识别率不高以及多导联脑电信号信息冗余和计算复杂度高的问题，有效去除多导联脑电信号引起的信息冗余，降低计算复杂度，提高抑郁检测的识别率。

本发明的技术方案是：

1.一种基于多导联脑电最优导联选择的抑郁检测系统，其特征在于，包括多导联脑电信号获取模块、脑电信号预处理模块、特征提取模块、导联选择模块、抑郁检测模块；所述多导联脑电信号获取模块用于获取静息状态采集的多导联脑电信号，所述脑电信号预处理模块用于对获取的多导联脑电信号进行预处理，所述特征提取模块用于提取进行最优导联选择的有效特征，所述导联选择模块用于采用修改版本的核排列mKTA作为目标函数，利用新的二元粒子群优化法NBPSO来对目标函数进行导联选择，求取最优化导联子集，所述抑郁检测模块用于对选择出来的最优化导联子集进行识别，构建最优化抑郁检测分类模型，来进行抑郁检测。

2.所述多导联脑电信号获取模块获取的64导脑电信号是在Brain Products平台上采集完成的，采集设备使用的是64导10-20国际标准电极帽，传感器的头皮阻抗低于20千欧，采样率设置为1000Hz，参考电极和接地电极分别设置为FCz和AFz。

3.所述多导联脑电信号获取模块获取的是健康被试和抑郁症患者在相同条件下的安静闭眼静息状态下5分钟的静息状态64导脑电信号。

4.脑电信号预处理模块对采集到的64导脑电信号，进行数据预处理，包括改变参考电极、降采样、带通滤波、去除伪迹四个步骤：

1)改变参考电极，是将原参考电极位FCz重设为两侧乳突部位的导联电位；

2)降采样，是将脑电信号的采样率由1000Hz降低至256Hz；

3)带通滤波是进行1Hz-40Hz的带通滤波以除去直流干扰以及高频噪声信号；

4)去除伪迹是对原始信号进行目视检查，挑选出含伪迹信号和不良电极等最小的70秒静息状态脑电信号，再对这些挑选出来的70秒静息状态脑电信号使用独立成分分析方法进一步去除伪迹信号。

5.所述特征提取模块包括线性特征提取的步骤：线性特征提取是对预处理后的脑电信号提取三个线性特征，包括功率谱最大功率、功率谱平均功率、功率谱中心频率；首先对每段脑电信号使用Welch方法求出功率谱，再使用对应段功率谱信号，求取其最大值，得到功率谱最大功率；求取其平均值，得到功率谱平均功率；求取中心频率，得到功率谱中心频率。

6.所述特征提取模块包括非线性特征提取的步骤：非线性特征提取是对预处理后的脑电信号提取两个非线性特征，包括柯尔莫哥洛夫熵、LZ复杂度；

柯尔莫哥洛夫熵通过若干增长嵌入维度的关联积分来表征系统的混沌程度：

其中m是嵌入维度，C_m(r)为该嵌入维度的关联积分；

LZ复杂度是以测量值的粗粒度化为基础求取信号复杂度的，经证明，它的上界为：

并且c(n)可以通过b(n)被归一化为：

其中n为信号长度，c(n)为信号的复杂度，b(n)为信号二值化后的复杂度。

7.所述导联选择模块包括最优化导联权重选择的步骤：把包含所有脑电的特征集划分为训练集和测试集，对训练集实施如下的导联权重选择步骤：首先将训练集中的特征数据按照导联组织特征矩阵f，由类标签y^T计算目标矩阵L；再由初始化的P个w，使用NBPSO方法生成P组导联选择权重w，计算P组导联子集S_i，并计算P组导联子集的导联选择核矩阵K_S，再使用mKTA作为目标函数来度量通过选择得到的导联子集的优劣，计算P组mKTA核排列值；接着，根据上述P组mKTA值，更新表示P组个体最优mKTA值对应的所有导联的权重p_best和表示全局最优mKTA值对应的所有导联的权重g_best，并使用NBPSO来更新上述P组导联权重并更新导联子集；再使用更新后的P组导联权重w更新P组导联子集S_i，并计算导联选择核矩阵K_S，如此进行迭代，直至输出最优化导联权重w。

8.所述抑郁检测模块包括构建最优化抑郁检测分类模型的步骤：首先，根据求出的最优化导联权重分别计算训练集和测试集相应的最优化导联子集训练集S_tr和最优化导联子集测试集S_te；使用十折交叉验证法进行抑郁检测参数寻优来构建最优化的抑郁检测分类模型：将求得的最优化导联子集训练集S_tr和训练类标签

随机划分为10等份，轮流使用其中的9份作为训练集训练模型，剩下的1份作为验证集验证模型，寻找最优化的抑郁检测模型参数；然后，使用寻找到的最优化参数在求得的最优化导联子集训练集上构建最优化的抑郁检测分类模型f(S_tr)，f(S_tr)＝W_modelk(S_tr,S_tr)+b，其中W_model是权重系数，k为核函数，b为偏置项。

9.所述抑郁检测模块包括对最优化抑郁检测分类模型测试的步骤：将求得的最优化导联子集测试集S_te输入构建的最优化的抑郁检测分类模型中：f(S_te)＝W_modelk(S_te,S_tr)+b，通过计算得到抑郁测试结果，比较测试类标签和最优化的抑郁检测模型输出的抑郁测试结果，并计算相应的准确率、敏感度、特异性等模型性能评价指标。

10.所述抑郁检测模块包括抑郁检测辅助诊断的步骤：首先对采集到的脑电信号数据进行预处理，并提取有效特征，得到特征矩阵f_new，使用最优化导联权重，计算辅助诊断最优化导联特征子集S_new，然后将求得的辅助诊断最优化导联特征子集输入最优化的抑郁检测分类模型，计算抑郁检测结果：f(S_new)＝W_modelk(S_new,S_tr)+b。

本发明的技术效果：

本发明提供的一种基于多导联脑电最优导联选择的抑郁检测系统，针对脑电信号在抑郁症临床检测中存在的识别率不高以及多导联脑电信号信息冗余和计算复杂度高的问题，有效去除多导联脑电信号引起的信息冗余，降低计算复杂度，提高抑郁检测的识别率。本发明也解决了目前存在的导联选择方法存在的解释性差、识别效果不佳、对导联组合考虑不够、计算损耗大以及易过拟合等问题。

本发明的基于多导联脑电最优导联选择的抑郁检测方法，使用修改版本的核排列(modified Kernel Target Alignment,mKTA)作为导联选择的目标函数，并通过新的二元粒子群优化(New Binary Particle Swarm Optimization,NBPSO)来优化选择的导联子集，最终解得最优化的导联子集。由于在求取最优化导联子集的过程中，会最小化mKTA值，换言之，就是最大化抑郁症患者和正常被试脑电在特征空间中的类间距，因此，使用经过本发明的导联选择方法选出的最优化的导联子集进行抑郁诊断会获得较为理想的抑郁诊断识别率，可以准确地区分开抑郁症患者和正常被试的脑电信号。同时，本发明提供的目标函数mKTA是用以度量核函数和目标函数的相似程度的，具有良好的解释性，同样，本发明提供的NBPSO也是简单有效的函数优化方法，也具有良好的解释性，而且，本发明提供的mKTA也可以有效地避免过拟合。本发明提供的NBPSO是从整体上来选择导联子集的，会充分考查不同导联间的组合关系来选择最优化的导联子集。众所周知，粒子群优化是一种进化方法，且可以并行地计算多组所有的导联权重，因此本发明提供的新的二元粒子群优化也可以解决传统导联选择算法计算损耗大的问题。

附图说明

图1是本发明的基于多导联脑电最优导联选择的抑郁检测系统的工作流程图。

图2是64导10-20国际标准电极位置图。

图3是导联选择模块进行导联权重选择的流程图。

图4是抑郁检测模块进行抑郁检测的流程图。

具体实施方式

以下结合附图和实施例对本发明作进一步详细说明。

一种基于多导联脑电最优导联选择的抑郁检测系统，包括：多导联脑电信号获取模块、脑电信号预处理模块、特征提取模块、导联选择模块、抑郁检测模块；多导联脑电信号获取模块用于获取静息状态采集的多导联脑电信号，脑电信号预处理模块用于对获取的多导联脑电信号进行预处理，特征提取模块用于提取进行最优导联选择的有效特征，导联选择模块用于采用修改版本的核排列作为目标函数，利用新的二元粒子群优化法来对目标函数进行导联选择，求取最优化导联子集，抑郁检测模块用于对选择出来的最优化导联子集进行识别，构建最优化抑郁检测分类模型，来进行抑郁检测。

图1是本发明的基于多导联脑电最优导联选择的抑郁检测系统的工作流程图。

本发明的基于多导联脑电最优导联选择的抑郁检测系统的工作流程包括如下5个阶段：

(1)多导联脑电信号获取模块的数据采集阶段，采集年龄、性别、教育程度匹配的健康被试和抑郁症患者的在相同条件下的安静闭眼静息状态64导脑电信号。

(2)脑电信号预处理模块的数据预处理阶段，对采集到的脑电信号进行改变参考电极、降采样、带通滤波、去除伪迹四个步骤的处理。

(3)特征提取模块的特征提取阶段，对预处理后的脑电信号提取三个线性特征，包括功率谱最大功率、功率谱平均功率、功率谱中心频率，以及两个非线性特征包括柯尔莫哥洛夫熵、LZ复杂度。

(4)导联选择模块的导联选择阶段，使用修改版本的核排列方法作为目标函数，利用新的二元粒子群优化算法来对该目标函数进行导联选择，求出最优化导联子集。

(5)抑郁检测模块的抑郁检测阶段，完成导联选择后，使用支持向量机等分类器对选择出来的最优化导联子集进行识别，构建最优化分类模型，来进行抑郁检测。

下面分别说明本发明基于多导联脑电最优导联选择的抑郁检测的5个阶段：

(1)数据采集阶段

数据采集阶段，我们通过心理量表以及医生诊断选择了15名抑郁症患者以及20名年龄、性别、教育程度匹配的健康被试。所有的被试都是右手利，视力正常或者矫正正常，听力、智力正常，无其他神经疾病史。在数据采集之前，所有的被试均未服用任何神经药物或精神药物，数据采集时，被试精神状态良好。所有被试的数据采集均在相同的环境条件下进行。

数据采集是在Brain Products平台上完成的，采集被试5分钟在安静闭眼静息状态下的静息状态脑电信号，采集设备使用的是64导10-20国际标准电极帽，电极位置如图2所示。传感器的头皮阻抗低于20千欧。采样率设置为1000Hz。参考电极和接地电极分别设置为FCz和AFz。

(2)数据预处理阶段

对采集到64导脑电信号，进行数据预处理。数据预处理主要包括改变参考电极、降采样、带通滤波、去除伪迹四个步骤。

由于采集脑电信号时参考电极设置为FCz，这导致前额叶、顶叶区域导联的信号幅值很低。因此，首先将原参考电极位重设为两侧乳突部位的导联电位，增大前额叶、顶叶区域导联的信号幅值，便于后续处理。因为脑电信号的原始采样频率为1000Hz，采样率过高会造成数据量过大而影响后续的处理和分析。因此，要对采集到的脑电信号进行降采样，将脑电信号的采样率由1000Hz降低至256Hz。本发明进行1Hz-40Hz的带通滤波以除去直流感染以及高频噪声信号。虽然本发明是在安静闭眼静息状态下采集脑电信号的，但是，依然很难避免眼电、眼动以及肌电信号带来的影响。因此需要对采集到的脑电信号进行去除伪迹。去除伪迹是对原始信号进行目视检查，挑选出含伪迹信号如：身体运动，眨眼信号，眼球运动信号，肌电信号和不良电极等最小的70秒静息状态脑电信号，再对这些挑选出来的70秒静息状态脑电信号使用独立成分分析方法进一步去除伪迹信号。

(3)特征提取阶段

本发明使用在抑郁诊断中常用并且被证明为有效的三个线性特征以及两个非线性特征对每一段经预处理后的脑电信号提取特征。

1)三个线性特征：包括功率谱最大功率、功率谱平均功率、功率谱中心频率；

对每段脑电信号使用Welch方法求出功率谱，再使用对应段功率谱信号，分别求取其最大值、平均值、中心频率，得到功率谱最大功率、功率谱平均功率、功率谱中心频率。

2)两个非线性特征：包括柯尔莫哥洛夫熵、LZ复杂度；

柯尔莫哥洛夫熵可以通过若干增长嵌入维度的关联积分来表征系统的混沌程度：

其中m是嵌入维度，C_m(r)为该嵌入维度的关联积分；

LZ复杂度是以测量值的粗粒度化为基础求取信号复杂度的，经证明，它的上界为：

并且c(n)可以通过b(n)被归一化为：

其中n为信号长度，c(n)为信号的复杂度，b(n)为信号二值化后的复杂度。

(4)导联选择阶段

首先把包含所有脑电的特征集划分为训练集和测试集。对训练集实施如下的导联权重选择步骤，如图3所示，是导联选择模块进行导联权重选择的流程图。首先将训练集中的特征数据按照导联进行组织特征矩阵f_i＝[c_i1,...,c_iE]，f＝[f₁,...,f_n]^T，其中c_ie表示第e个导联的特征数据，E表示脑电信号的导联数量，本实施例E为64，n为训练集中的脑电信号数，n表示脑电信号的段数；由类标签y^T得到目标矩阵L，L＝yy^T,y＝[y₁,...,y_n]^T,y_i∈{+1,-1}，y_i表示类标签；再由随机初始化的P个w，使用NBPSO方法生成P组导联选择权重w，计算P组导联子集S_i，S_i＝wof_i(o表示哈达玛积)，并计算P组导联子集的导联选择核矩阵K_S，K_Sij＝k(S_i,S_j)＝k(wof_i,wof_j)，k表示核函数。

使用修改版本的核排列(modified Kernel Target Alignment,mKTA)作为目标函数来度量通过选择得到的导联子集的优劣，计算P组mKTA值：

其中，L是目标矩阵，L＝yy^T,y＝[y₁,...,y_n]^T,y_i∈{+1,-1}，y_i表示类标签，n表示脑电信号的段数，f＝[f₁,...,f_n]^T,f_i＝[c_i1,...,c_iE],w＝[w₁,...,w_E],w_e∈{0,1},S_i＝wof_i，f_i是第i段脑电信号所有导联的特征，c_ie表示第i段脑电信号第e个导联的特征，w_e表示第e个导联的权重，当w_e＝1时，表明第e个导联被选择，而当w_e＝0时，表明第e个导联未被选择，w表示所有导联的权重(因此，w可以被认为是选择出的导联子集)，S_i表示第i段脑电信号经过导联选择后的特征子集，其中未被选中的导联的特征值被置为0而被选中的导联的特征值不变。K_S是用于导联选择的核矩阵，K_Sij＝k(S_i,S_j)＝k(wof_i,wof_j)。使用mKTA作为目标函数，可以通过优化被选择的导联子集，最小化目标矩阵L与导联选择核矩阵K_S之间的差异。根据mKTA定义，mKTA值越小，表明目标矩阵L与导联选择核矩阵K_S之间的相似度越高，也就意味着经选择的导联子集在抑郁诊断中能够获得良好的识别率。因此，最小化mKTA值也就等价于在特征空间中最大化类间距离。

接着，根据上述P组mKTA值，使用新的二元粒子群优化(New Binary ParticleSwarm Optimization,NBPSO)更新表示P组个体最优mKTA值对应的所有导联的权重p_best和表示全局最优mKTA值对应的所有导联的权重g_best，并使用NBPSO来更新上述P组导联权重并更新导联子集：

S(v_t+1)＝2×|Sig(v_t+1)-0.5|

其中，w_t表示第t轮迭代时的所有导联的权重(导联子集)，v_t表示第t轮迭代时的所有导联的权重w_t的速度，v_t+1表示经过更新后的第t+1轮迭代时的所有导联的权重w_t的速度，

是惯性系数常量，c₁、c₂是正常数，rand()和Rand()是以[0,1]为区间的随机函数，p_best表示个体最优mKTA值对应的所有导联的权重，g_best表示全局最优mKTA值对应的所有导联的权重，Sig(v_t+1)表示速度v_t+1的Sigmod函数，S(v_t+1)是对速度v_t+1进行转换的转换函数，exchange(w_t)表示对第t轮所有导联的权重w_t值进行转换，将1转换为0，反之亦然，w_t+1是经过更新后的第t+1轮迭代时的所有导联的权重。使用NBPSO选择导联子集，求出该子集对应的mKTA的值来更新个体最优mKTA值以及全局最优mKTA值，并根据个体最优mKTA值以及全局最优mKTA值所对应的所有导联的权重来更新当前选择的导联子集的下一轮更新的速度v_t+1，并对该速度进行转换，再通过该速度更新下一轮所有导联的权重w_t+1直到求出最优化的导联权重。本发明提供的NBPSO是一种并行化的优化方法，可以并行地计算和更新多组导联权重，如在初始化时构建P组导联权重，就可以在更新迭代中并行地计算这P组导联权重，直至求出最优化导联权重。

使用更新后的P组导联权重更新P组导联子集S_i，并计算导联选择核函数K_S，如此进行迭代，直至求出最优化导联子集S，本发明实施例求出的最优化导联子集为Fp1，C4，P3，O2，F7，CP1，FC5，TP9，F1，F2，C1，FC3，F5，C5，AF7，FT8，TP7，FT9，FT10，Fpz。这些导联所在区域与已知的抑郁相关的脑部区域一致。

(5)抑郁检测阶段

图4是抑郁检测模块进行抑郁检测的流程图。

根据阶段(4)求出的最优化导联权重得出最优化导联子集，分别计算训练集和测试集相应的最优化导联子集S_tr＝wof_tr和S_te＝wof_te，其中w是求得的最优化导联权重，f_tr是训练特征集，S_tr是最优化导联子集训练集，f_te是测试特征集，S_te是最优化导联子集测试集。在辅助诊断中，使用最优化导联权重对辅助诊断脑电特征计算辅助诊断最优化导联子集S_new＝wof_new，其中f_new是辅助诊断特征集，S_new是辅助诊断最优化导联子集特征集。类标签y^T划分为测试类标签

和训练类标签

首先包括构建最优化抑郁检测分类模型的步骤：使用十折交叉验证法进行抑郁检测参数寻优来构建最优化的抑郁检测分类模型。将求得的最优化导联子集训练集S_tr和训练类标签随机划分为10等份，轮流使用其中的9份作为训练集训练模型，剩下的1份作为验证集验证模型，寻找最优化的抑郁检测模型参数，换言之，每份数据都会被当作验证集来验证使用其他数据生成的训练模型的性能，通常，选择平均验证精度最高的模型参数作为最优化抑郁检测模型参数；然后，使用寻找到的最优化参数在求得的最优化导联子集训练集上构建最优化的抑郁检测分类模型f(S_tr)，本发明使用的是支持向量机分类模型，通过该分类器构建的分类模型为：

f(S_tr)＝W_modelk(S_tr,S_tr)+b，其中W_model是权重系数，k为核函数，b为偏置项。

其次，包括对最优化抑郁检测分类模型测试的步骤：将求得的最优化导联子集测试集S_te输入构建的最优化的抑郁检测分类模型中，f(S_te)＝W_modelk(S_te,S)+b，通过计算得到抑郁测试结果，比较测试类标签

和最优化的抑郁检测模型输出的抑郁测试结果，并计算相应的准确率、敏感度、特异性等模型性能评价指标。

然后，包括抑郁检测辅助诊断的步骤：首先对采集到的脑电信号数据进行预处理，并提取有效特征，得到特征矩阵f_new，使用最优化导联权重，计算辅助诊断最优化导联特征子集S_new，然后将求得的辅助诊断最优化导联特征子集输入最优化的抑郁检测分类模型，通过计算可以得到抑郁检测结果：f(S_new)＝W_modelk(S_new,S)+b。

在导联选择阶段和抑郁检测阶段，本发明使用留一被试交叉验证方法。本发明的抑郁检测识别率为80％，敏感度为80％，特异性为80％，而不使用本发明进行导联选择的抑郁识别率仅为65.71％。

应当指出，以上所述具体实施方式可以使本领域的技术人员更全面地理解本发明创造，但不以任何方式限制本发明创造。一切不脱离本发明创造的精神和范围的技术方案及其改进，其均涵盖在本发明创造的保护范围当中。