【机器学习】面向自闭症辅助诊断的无监督模糊特征学习新方法

第14卷第5期 智能系统学报 Vol.14 No.5 2019年9月 CAAI Transactions on Intelligent Systems Sep.2019 D0:10.11992/tis.201808005 网络出版地址：http:/kns.cnki.net/kcms/detail/23.1538.tp.20190527.1715.012.html 面向自闭症辅助诊断的无监督模糊特征学习新方法 张英，王骏，鲍国强，张春香，王士同 （江南大学数字媒体学院，江苏无锡214122) 摘要：自闭症患者的行为和认知缺陷与潜在的脑功能异常有关。对于静息态功能磁振图像(functional magnetic resonance imaging,fMRI)高维特征，传统的线性特征提取方法不能充分提取其中的有效信息用于分类。为此， 本文面向MRI数据提出一种新型的无监督模糊特征映射方法，并将其与多视角支持向量机相结合，构建分类 模型应用于自闭症的计算机辅助诊断。该方法首先采用多输出TSK模糊系统的规则前件学习方法，将原始特 征数据映射到线性可分的高维空间：然后引入流形正则化学习框架，提出新型的无监督模糊特征学习方法，从 而得到原输出特征向量的非线性低维嵌入表示：最后使用多视角SVM算法进行分类。实验结果表明：本文方 法能够有效提取静息态MRI数据中的重要特征，在保证模型具有优越且稳定的分类性能的前提下，还可以提 高模型的可解释性。 关键词：自闭症；功能磁共振成像；功能连接：皮尔森相关性；特征选择：无监督模糊特征映射；流形正则化框 架；支持向量机 中图分类号：TP391文献标志码：A文章编号：1673-4785(2019)05-0882-07 中文引用格式：张英，王骏，鲍国强，等.面向自闭症辅助诊断的无监督模糊特征学习新方法智能系统学报，2019,14(5)： 882-888. 英文引用格式：ZHANG Ying,WANGJun,BAO Guoqiang,et al.A novel unsupervised fuz四y feature learning method for com puter-aided diagnosis of autism[J].CAAI transactions on intelligent systems,2019,14(5):882-888 A novel unsupervised fuzzy feature learning method for computer-aided diagnosis of autism ZHANG Ying,WANG Jun,BAO Guoqiang,ZHANG Chunxiang,WANG Shitong (School of Digital Media,Jiangnan University,Wuxi 214122,China) Abstract:Studies have shown that the behavioral and cognitive defect of patients with autism have a close relationship with potential brain dysfunction.For the high-dimensional rs-fMRI features,traditional linear feature extraction method cannot always discriminatively extract the important information for classification.To this end,a novel method for fMRI data based on both unsupervised fuzzy feature mapping and multi-view support vector machine is proposed in this study,which aims to build a classification model for computer aided diagnosis of autism.In this method,the original features are first mapped to a linear separable high-dimensional space using the rule precursor learning method of multi- output Takagi-Sugeno-Kang(TSK)fuzzy system;then the manifold regularization learning framework is introduced.On the basis of this,a novel unsupervised fuzzy feature learning method is used to obtain the nonlinear low-dimensional embedding representation of the original output eigenvector.Finally,a multi-view support vector machine(SVM)al- gorithm is used for classification.The experimental results show that the proposed method can effectively extract im- portant features from the rs-fMRI data and improve the interpretability of the model on the premise of ensuring a superi- or and stable classification performance of the model. Keywords:autism;functional magnetic resonance imaging;functional connectivity;Pearson's correlation; feature selection;unsupervised fuzzy feature mapping;manifold regularization framework;support vector machine 收稿日期：2018-08-08.网络出版日期：2019-05-29 自闭症又叫孤独症，是一种广泛性神经失调 基金项目：江苏省自然科学基金项目(BK20181339). 通信作者：王骏.E-mail:wangjun_sytu@hotmail.com 的发展障碍疾病，其临床表现主要为社会交往障

DOI: 10.11992/tis.201808005 网络出版地址: http://kns.cnki.net/kcms/detail/23.1538.tp.20190527.1715.012.html 面向自闭症辅助诊断的无监督模糊特征学习新方法 张英，王骏，鲍国强，张春香，王士同 （江南大学 数字媒体学院，江苏 无锡 214122） 摘 要：自闭症患者的行为和认知缺陷与潜在的脑功能异常有关。对于静息态功能磁振图像 (functional magnetic resonance imaging, fMRI) 高维特征，传统的线性特征提取方法不能充分提取其中的有效信息用于分类。为此， 本文面向 fMRI 数据提出一种新型的无监督模糊特征映射方法，并将其与多视角支持向量机相结合，构建分类 模型应用于自闭症的计算机辅助诊断。该方法首先采用多输出 TSK 模糊系统的规则前件学习方法，将原始特 征数据映射到线性可分的高维空间；然后引入流形正则化学习框架，提出新型的无监督模糊特征学习方法，从 而得到原输出特征向量的非线性低维嵌入表示；最后使用多视角 SVM 算法进行分类。实验结果表明：本文方 法能够有效提取静息态 fMRI 数据中的重要特征，在保证模型具有优越且稳定的分类性能的前提下，还可以提 高模型的可解释性。 关键词：自闭症；功能磁共振成像；功能连接；皮尔森相关性；特征选择；无监督模糊特征映射；流形正则化框 架；支持向量机 中图分类号：TP391 文献标志码：A 文章编号：1673−4785(2019)05−0882−07 中文引用格式：张英, 王骏, 鲍国强, 等. 面向自闭症辅助诊断的无监督模糊特征学习新方法 [J]. 智能系统学报, 2019, 14(5): 882–888. 英文引用格式：ZHANG Ying, WANG Jun, BAO Guoqiang, et al. A novel unsupervised fuzzy feature learning method for com￾puter-aided diagnosis of autism[J]. CAAI transactions on intelligent systems, 2019, 14(5): 882–888. A novel unsupervised fuzzy feature learning method for computer-aided diagnosis of autism ZHANG Ying，WANG Jun，BAO Guoqiang，ZHANG Chunxiang，WANG Shitong (School of Digital Media, Jiangnan University, Wuxi 214122, China) Abstract: Studies have shown that the behavioral and cognitive defect of patients with autism have a close relationship with potential brain dysfunction. For the high-dimensional rs-fMRI features, traditional linear feature extraction method cannot always discriminatively extract the important information for classification. To this end, a novel method for fMRI data based on both unsupervised fuzzy feature mapping and multi-view support vector machine is proposed in this study, which aims to build a classification model for computer aided diagnosis of autism. In this method, the original features are first mapped to a linear separable high-dimensional space using the rule precursor learning method of multi￾output Takagi-Sugeno-Kang (TSK) fuzzy system; then the manifold regularization learning framework is introduced. On the basis of this, a novel unsupervised fuzzy feature learning method is used to obtain the nonlinear low-dimensional embedding representation of the original output eigenvector. Finally, a multi-view support vector machine (SVM) al￾gorithm is used for classification. The experimental results show that the proposed method can effectively extract im￾portant features from the rs-fMRI data and improve the interpretability of the model on the premise of ensuring a superi￾or and stable classification performance of the model. Keywords: autism; functional magnetic resonance imaging; functional connectivity; Pearson’s correlation; feature selection; unsupervised fuzzy feature mapping; manifold regularization framework; support vector machine 自闭症又叫孤独症，是一种广泛性神经失调 的发展障碍疾病，其临床表现主要为社会交往障 收稿日期：2018−08−08. 网络出版日期：2019−05−29. 基金项目：江苏省自然科学基金项目 (BK20181339). 通信作者：王骏. E-mail：wangjun_sytu@hotmail.com. 第 14 卷第 5 期 智 能 系 统 学 报 Vol.14 No.5 2019 年 9 月 CAAI Transactions on Intelligent Systems Sep. 2019

第5期 张英，等：面向自闭症辅助诊断的无监督模糊特征学习新方法 ·883· 碍、言语交流障碍和动作刻板重复等。自1943 表1NYU实验对象信息 年利奥·坎纳医生首次对自闭症做出经典性描述 Table 1 The information of NYU experimental object 以来，该疾病因其病因不明确，又集多种障碍于 实验对象信息 自闭症患者数 正常对照组 一身，已经得到了广泛的关注。目前，在美国平 人数（男女） 68/11 7726 均每68名儿童中就有1名自闭症儿童；在世界范 年龄 7.13-39.1 6.47-31.78 围内，自闭症患者超过1%。据统计，约2/3孤 表2USM实验对象信息 独症患儿成年后无法独立生活，需要终生照顾和 Table 2 The information of USM experimental object 看护。为此，人们迫切希望寻找科学而高效的手 实验对象信息 自闭症患者数 正常对照组 段对自闭症进行准确诊断以达到早发现早治疗的 人数（男/女） 56/0 43/0 目的。 年龄 11.442.3 8.8-39.4 已有研究表明，自闭症患者行为和认知缺陷 与潜在的脑功能异常有关。静息态MRI技术 每个对象的fMRI数据使用DPARSF(data pro 通过使用血氧水平依赖(blood oxygenation level cessing assistant for resting-state fMRI)lul工具进行 dependent,.BOLD)信号可以有效地反映患者在静 处理。处理步骤主要包括：1)去除功能磁共振图 息状态下脑部血流和代谢活动等功能性变化，因 像序列前10个时间点的数据；2)时间层校正和头 部运动校正；3)数据通过统一使用T1加权图像 此在自闭症的辅助诊断中起着重要的作用。基于 分割并将其归一化到MNIl52(montreal neurologic- 静息态的医学影像诊断大致可分为特征提取和分 al institute152)标准空间中；4)使用AAL(anatom- 类器设计这两个阶段。在特征提取阶段，通过使 ical automatic labeling)模板将大脑划分成Il6个 用相关性特征选择等方法从高维医学影像特征中 脑区，每个区域重新采样3×3×3mm3的体素；5) 得到有鉴别意义的特征；在分类器设计阶段，通 使用半高全宽高斯核进行空间平滑处理；6)应用 过使用线性或非线性分类器进行分类。但是目前 带通滤波(0.01~0.1Hz)去除噪声；7)去线性漂移 多采用线性特征学习方法，虽然具有简单高效等 并进行全局信号校正去干扰变量：8)计算每个脑 优点，但是尚不能有效挖掘出医学影像数据中的 区的平均时间序列。 复杂信息。 2基于fMRI数据的无监督模糊特征 本文基于静息态fMRI提出了一种新型的无 学习方法 监督模糊特征学习方法，并在此基础上构造多视 角分类器从而实现自闭症的辅助诊断。其技术路 2.1方法流程 线可概括为：首先，将静息态fMRI经过预处理 基于静息态fMRI进行自闭症辅助诊断包括 后，基于Pearson相关系数构造低阶功能连接矩 训练和测试两个阶段。本文所采用的训练过程如 阵，并使用文献[7]的方法构造高阶功能连接矩 图1所示，各步骤包括： 阵：然后，提出新型的无监督模糊特征学习方法， 1)根据预处理后的静息态功能磁共振成像 数据，提取出各脑区的平均时间序列信号，计算 将2种功能连接矩阵分别映射到不同的低维嵌入 脑区之间的Pearson系数，得到低阶功能连接矩阵 空间中；最后，以低阶和高阶功能连接在脑疾病 M。 诊断中的互补性为出发点，设计多视角SVMo分 2)将低阶功能连接矩阵的每一行作为各脑区 类器进行分类。 的特征描述，再次计算脑区之间的Pearson系数， 1实验数据集及处理 得到高阶功能连接矩阵M。 3)分别取低阶和高阶功能连接矩阵M和 本文使用的数据集来自ABIDE(autism brain M的上三角阵，按行串联形成新的特征向量，然 imaging datae xchange))数据库(http://preprocessed- 后将所有对象的特征向量进行排列，生成基于低 connectomes-project.org/abide/download.html) 阶和高阶功能连接的所有对象的特征矩阵；对两 NYU子集和USM子集，采集对象的FIQ(full-scale 种不同的特征矩阵计算各特征与类标的相关性， intelligence quotient)范围为65~l50。实验所处理 选择与类标相关性最高的D个特征，构成矩阵 对象的相关信息如表1和表2所示。 X=（)(x);…）月）eRw和X=(x)）:

碍、言语交流障碍和动作刻板重复等[1-3]。自 1943 年利奥·坎纳医生首次对自闭症做出经典性描述 以来，该疾病因其病因不明确，又集多种障碍于 一身，已经得到了广泛的关注。目前，在美国平 均每 68 名儿童中就有 1 名自闭症儿童；在世界范 围内，自闭症患者超过 1%。据统计，约 2/3 孤 独症患儿成年后无法独立生活，需要终生照顾和 看护。为此，人们迫切希望寻找科学而高效的手 段对自闭症进行准确诊断以达到早发现早治疗的 目的。 已有研究表明，自闭症患者行为和认知缺陷 与潜在的脑功能异常有关。静息态 fMRI 技术[4-9] 通过使用血氧水平依赖 (blood oxygenation level dependent, BOLD) 信号可以有效地反映患者在静 息状态下脑部血流和代谢活动等功能性变化，因 此在自闭症的辅助诊断中起着重要的作用。基于 静息态的医学影像诊断大致可分为特征提取和分 类器设计这两个阶段。在特征提取阶段，通过使 用相关性特征选择等方法从高维医学影像特征中 得到有鉴别意义的特征；在分类器设计阶段，通 过使用线性或非线性分类器进行分类。但是目前 多采用线性特征学习方法，虽然具有简单高效等 优点，但是尚不能有效挖掘出医学影像数据中的 复杂信息。 本文基于静息态 fMRI 提出了一种新型的无 监督模糊特征学习方法，并在此基础上构造多视 角分类器从而实现自闭症的辅助诊断。其技术路 线可概括为：首先，将静息态 fMRI 经过预处理 后，基于 Pearson 相关系数构造低阶功能连接矩 阵，并使用文献 [7] 的方法构造高阶功能连接矩 阵；然后，提出新型的无监督模糊特征学习方法， 将 2 种功能连接矩阵分别映射到不同的低维嵌入 空间中；最后，以低阶和高阶功能连接在脑疾病 诊断中的互补性为出发点，设计多视角 SVM[10] 分 类器进行分类。 1 实验数据集及处理 本文使用的数据集来自 ABIDE(autism brain imaging datae xchange) 数据库 (http://preprocessed￾connectomes-project.org/abide/download.html) 的 NYU 子集和 USM 子集，采集对象的 FIQ(full-scale intelligence quotient) 范围为 65～150。实验所处理 对象的相关信息如表 1 和表 2 所示。 表 1 NYU 实验对象信息 Table 1 The information of NYU experimental object 实验对象信息 自闭症患者数 正常对照组 人数(男/女) 68/11 77/26 年龄 7.13~39.1 6.47~31.78 表 2 USM 实验对象信息 Table 2 The information of USM experimental object 实验对象信息 自闭症患者数 正常对照组 人数(男/女) 56/0 43/0 年龄 11.4~42.3 8.8~39.4 每个对象的 fMRI 数据使用 DPARSF(data pro￾cessing assistant for resting-state fMRI)[11] 工具进行 处理。处理步骤主要包括：1）去除功能磁共振图 像序列前 10 个时间点的数据；2）时间层校正和头 部运动校正；3）数据通过统一使用 T1 加权图像 分割并将其归一化到 MNI152(montreal neurologic￾al institute 152) 标准空间中；4）使用 AAL (anatom￾ical automatic labeling) 模板将大脑划分成 116 个 脑区，每个区域重新采样 3×3×3 mm3 的体素；5） 使用半高全宽高斯核进行空间平滑处理；6）应用 带通滤波 (0.01～0.1 Hz) 去除噪声；7）去线性漂移 并进行全局信号校正去干扰变量；8）计算每个脑 区的平均时间序列。 2 基于 fMRI 数据的无监督模糊特征 学习方法 2.1 方法流程 基于静息态 fMRI 进行自闭症辅助诊断包括 训练和测试两个阶段。本文所采用的训练过程如 图 1 所示，各步骤包括： Ml tr 1）根据预处理后的静息态功能磁共振成像 数据，提取出各脑区的平均时间序列信号，计算 脑区之间的 Pearson 系数，得到低阶功能连接矩阵 。 Mh tr 2）将低阶功能连接矩阵的每一行作为各脑区 的特征描述，再次计算脑区之间的 Pearson 系数， 得到高阶功能连接矩阵 。 Ml tr Mh tr X l tr = (( x tr,l 1 )T ; ( x tr,l 2 )T ;··· ; ( x tr,l N )T ) ∈ R N×D X h tr = (( x tr,h 1 )T ; 3）分别取低阶和高阶功能连接矩阵 和 的上三角阵，按行串联形成新的特征向量，然 后将所有对象的特征向量进行排列，生成基于低 阶和高阶功能连接的所有对象的特征矩阵；对两 种不同的特征矩阵计算各特征与类标的相关性， 选择与类标相关性最高的 D 个特征，构成矩阵 和 第 5 期 张英，等：面向自闭症辅助诊断的无监督模糊特征学习新方法 ·883·

·884· 智能系统学报 第14卷 (x)(x)月）eRwo,其中xeRo,x∈R,i= 得到相应的嵌入矩阵E和。 1,2,…,N。 5)基于4)得到的嵌入矩阵E和E础，分别计 算相应的核矩阵。 4)使用无监督模糊特征映射方法对低阶和高 6)将2种不同的核矩阵进行线性组合得到复 阶功能连接的特征矩阵进行特征学习，得到相应 合核矩阵。 的变换矩阵B。,并将：和冲变换到低维空间中 7)构造相应的SVM分类器。 样本1,2，…，N 低阶 ASD ⑥ ⑦ NC 高阶 图1无监督模糊特征学习方法的框架流程图 Fig.1 The framework of an unsupervised fuzzy feature learning method 22低阶和高阶功能连接矩阵的生成 式中：A表示输入特征向量x第d维特征所对应 虽然大脑各脑区在空间上相互隔离，但它们 的第k条模糊规则的模糊子集；A表示模糊合取 之间的神经活动相互影响，通常使用功能连接矩 操作；k=1,2,…,K;d=1,2,…,D,d表示统称； 阵来表示脑区之间的时间相关性。本文使用文 s=1,2,…,S;K表示模糊规则数；S表示模型输出 献[8]中的方法，基于每个脑区fMRI时间序列来 的个数。 计算每个对象各脑区之间的Pearson相关系数，得 本文采用高斯函数表示模糊规则的模糊集，即 到低阶功能连接矩阵。将低阶功能连接矩阵 的每一行作为一个脑区的特征描述，使用同样的 A;(xa)=exp -(xa-c) 26 (1) 方法再次计算脑区之间的Pearson相关系数，从而 得到高阶功能连接矩阵M。 并使用式(2)、式(3)来计算每个样本相对于 2.3基于流形正则化约束的模糊特征学习 每条模糊规则的点火强度： 为了在无监督环境中将高维空间中的样本点 (2) 映射到S维的低维空间，本文引人多输出TSK模 -u 糊建模方法。TSK模糊系统是由形式为“IF- t(x)= 止() THEN”的模糊规则组成的，对于多输出TSK模糊 (3) 系统，给定任意一个D维输入向量x= 2rw [x2…xJ「，第k条模糊规则可以表示为： 参数c告和的值可以使用聚类技术或其他 if x is Af A x2 is AA...AxpA Ab,then 划分方法计算得出。本文使用FCM聚类算 f(x)=+x+…+哈xo 法进行数据集的模糊划分，其中c和可通过 5x)=哈2+2x+…+路xD 式(4)、式（⑤）计算得到： (4) f(x)=哈+x1+…+哈xD

( x tr,h 2 )T ;··· ; ( x tr,h N )T ) ∈ R N×D x tr,l i ∈ R D x tr,h i ∈ R D 1,2,··· ,N ，其中 ， ，i = 。 βg X l tr X h tr 4）使用无监督模糊特征映射方法对低阶和高 阶功能连接的特征矩阵进行特征学习，得到相应 的变换矩阵 ，并将 和 变换到低维空间中 E l tr E h 得到相应的嵌入矩阵 和 tr。 E l tr E h 5）基于 4）得到的嵌入矩阵 和 tr，分别计 算相应的核矩阵。 6）将 2 种不同的核矩阵进行线性组合得到复 合核矩阵。 7）构造相应的 SVM 分类器。 样本1,2,…,N … … 低阶 高阶 … … ASD O x y NC 3 3 1 2 4 4 5 6 + 7 5 图 1 无监督模糊特征学习方法的框架流程图 Fig. 1 The framework of an unsupervised fuzzy feature learning method 2.2 低阶和高阶功能连接矩阵的生成 Ml trMh tr 虽然大脑各脑区在空间上相互隔离，但它们 之间的神经活动相互影响，通常使用功能连接矩 阵来表示脑区之间的时间相关性[12]。本文使用文 献 [8] 中的方法，基于每个脑区 fMRI 时间序列来 计算每个对象各脑区之间的 Pearson 相关系数，得 到低阶功能连接矩阵 。将低阶功能连接矩阵 的每一行作为一个脑区的特征描述，使用同样的 方法再次计算脑区之间的 Pearson 相关系数，从而 得到高阶功能连接矩阵 。 2.3 基于流形正则化约束的模糊特征学习 x = [x1 x2 ··· xD] T 为了在无监督环境中将高维空间中的样本点 映射到 S 维的低维空间，本文引入多输出 TSK 模 糊建模方法。TSK 模糊系统是由形式为“IF￾THEN”的模糊规则组成的，对于多输出 TSK 模糊 系统，给定任意一 个 D 维输入向量 ，第 k 条模糊规则可以表示为： if x1 is A k 1 ∧ x2 is A k 2 ∧ ··· ∧ xD ∧ A k D ,then    f k 1 (x) = β k,1 0 +β k,1 1 x1 +···+β k,1 D xD f k 2 (x) = β k,2 0 +β k,2 1 x1 +···+β k,2 D xD . . . f k S (x) = β k,S 0 +β k,S 1 x1 +···+β k,S D xD A k d x d k ∧ k = 1,2,··· ,K d = 1,2,··· ,D s = 1,2,··· ,S K S 式中： 表示输入特征向量 第 维特征所对应 的第 条模糊规则的模糊子集； 表示模糊合取 操作； ； ， d 表示统称； ； 表示模糊规则数； 表示模型输出 的个数。 本文采用高斯函数表示模糊规则的模糊集，即 µA k d (xd) = exp   − ( xd −c k d )2 2δ k d   (1) 并使用式 (2)、式 (3) 来计算每个样本相对于 每条模糊规则的点火强度： µ k (x) = ∏D d=1 µA k d (xd) (2) µ˜ k (x) = µ k (x) ∑K k=1 µ k (x) (3) c k d δ k d c k d δ k d 参数 和 的值可以使用聚类技术或其他 划分方法[13-14] 计算得出。本文使用 FCM 聚类算 法进行数据集的模糊划分，其中 和 可通过 式 (4)、式 (5) 计算得到： c k d = ∑N i=1 uik xid/ ∑N i=1 uik (4) ·884· 智 能 系 统 学 报 第 14 卷

第5期 张英，等：面向自闭症辅助诊断的无监督模糊特征学习新方法 ·885· =h∑w-2 N 开始 (5) 式中：“：表示每个样本x属于第k类的模糊隶属 输人：X=(：X2;;x)∈R 度；h是隶属度函数的可调尺度参数。多输出 TSK模糊系统的输出可以表示为 引人拉普拉斯矩阵L 模糊映射得到高维数据矩阵X y3= i(x)f (x) (6) 式中s=1,2,…,S。 构造目标函数：minl川A=l+t(X:LXP,) 慧ERo 令X=(x:x:…;xW),1为全1列向量，多输 出TSK模糊系统的建模过程可用式(7)(13)的 计算方法表述： 解此优化问题：（仙x+LXX)u=XX4 X.=(1,X) (7) *=diag(t(x),t(x),…,t(rwX。 (8) 得到最优解：B。=Xg2,,…,as) x=(到)…（月） (9) B=映…P瑞) 计算出嵌入矩阵：E=XB。 (10) B=(82….0B-门 (11) 结束 B=BB…B] (12) 图2无监督模糊特征学习流程图 f(X)=X B: (13) Fig.2 Flowchart of unsupervised fuzzy feature learning 式(7)(13)表述的多输出TSK模糊系统建模 在模型训练过程中，本文分别对低阶功能连 过程可以看作将K(D+1)维空间中的数据通过变 接数据X=(x;x；…;）)eRw和高阶功能连 换矩阵B。降维到S维数据空间的过程，通常 接数据X=(x;x;…;x)eR0采用上述过程 SN。式(14)可以转化为式 =E(E) (20) (15)所描述的广义特征值问题： 使用式(21)对核矩阵进行线性组合，得到最 (IN +LX X)u=YXX:u (15) 终的复合核矩阵K.: 在得到最小S个特征值y1,Y2,…,Ys(y1≤Y2≤…≤ K.=0K+02K (21) ys)及其对应的广义特征向量1,2，…，us(因 式中：0，和2是复合核矩阵的组合系数，满足0+02=1。 1=0,故不用)后，可以求得最优解为： 3实验结果与分析 Bg=Xg(a2,m3…,as+1) (16) 式中：u=u/川XXual表示归一化特征向量，d=2, 3.1性能指标及参数设置 3,…,S+1。 本节通过实验验证本文方法的有效性并与已 无监督模型模糊特征学习的过程如图2所示。 有方法进行比较。本文选取精确度(ACC)、AUC

δ k d = h · ∑N i=1 uik( xid −c k d )2 / ∑N i=1 uik (5) uik xi k h 式中： 表示每个样本 属于第 类的模糊隶属 度 ； 是隶属度函数的可调尺度参数。多输出 TSK 模糊系统的输出可以表示为 ys = ∑K k=1 µ˜ k (x)· f k s (x) (6) 式中 s = 1,2,··· ,S。 X = ( x T 1 ; x T 2 ;··· ; x T N ) 令 ，1 为全 1 列向量，多输 出 TSK 模糊系统的建模过程可用式 (7)~ (13) 的 计算方法表述[15] ： Xe = (1,X) (7) X˜ k = diag( µ˜ k (x1),µ˜ k (x2),··· ,µ˜ k (xN) ) Xe (8) Xg = (( X˜ 1 )T( X˜ 2 )T ···( X˜ K )T )T (9) β s,k = ( β s,k 0 , βs,k 1 ,··· , βs,k D )T (10) β s = (( β s,1 )T , ( β s,2 )T ,··· , ( β s,K )T )T (11) βg = [ β 1 β 2 ··· β S ] (12) f (X) = Xgβg (13) K(D+1) βg S N 本文中所涉及的场景是典型的高维小样本问 题，故满足 。式 (14) 可以转化为式 (15) 所描述的广义特征值问题： ( IN +λLXgX T g ) u = γXgX T g u (15) γ1, γ2,··· , γS (γ1 ⩽ γ2 ⩽ ··· ⩽ γS ) u1,u2,··· ,uS u1 = 0 在得到最小S个特征值 及其对应的广义特征向量 ( 因 ，故不用) 后，可以求得最优解为： βg = X T g (u˜ 2,u˜ 3,··· ,u˜ S+1) (16) u˜ d = ud/||XgX T g ud|| 3,··· ,S +1 式中： 表示归一化特征向量，d = 2, 。 无监督模型模糊特征学习的过程如图 2 所示。 βg∈R(D+1)K×S 构造目标函数：min || βF ||F+λtr (βg TXg T LXg βg ) 解此优化问题：(IN+λLXg Xg T ) u=γXg Xg Tu 开始 结束 引入拉普拉斯矩阵 L 模糊映射得到高维数据矩阵 Xg 计算出嵌入矩阵：E =Xg βg 得到最优解：βg =Xg T ( ~ 2 ,u ~ 3 ,···,u ~ S+1) 输入：X = (x1 ; x2 ; ···; xN )∈RN×D u 图 2 无监督模糊特征学习流程图 Fig. 2 Flowchart of unsupervised fuzzy feature learning X l tr = ( x tr,l 1 ; x tr,l 2 ;··· ; x tr,l N ) ∈ R N×D X h tr = ( x tr,h 1 ; x tr,h 2 ;··· ; x tr,h N ) ∈ R N×D β l g β h g 在模型训练过程中，本文分别对低阶功能连 接数据 和高阶功能连 接数据 采用上述过程 进行处理，分别将它们映射到 S 维的低维空间中 去，记相应的变换矩阵分别为 和 ，则低维空 间中的嵌入矩阵由式 (17) 和式 (18) 计算而得： E l tr= ( X l tr) g β l g (17) E h tr= ( X h tr) g β h g (18) 2.4 基于功能连接矩阵的多视角 SVM E l tr E h tr K l tr K h tr 基于上述模糊特征学习过程得到低阶和高阶 嵌入矩阵 和 ，本文采用文献 [16] 中的方法 构造 SVM 分类器。对于拟构造的 SVM 分类器， 采用线性核计算其核矩阵 和 ： K l tr = E l tr( E l tr)T (19) K h tr = E h tr( E h tr)T (20) Ktr 使用式 (21) 对核矩阵进行线性组合，得到最 终的复合核矩阵 ： Ktr = θ1K l tr +θ2K h tr (21) 式中： θ1和 θ2是复合核矩阵的组合系数，满足 θ1+θ2 = 1。 3 实验结果与分析 3.1 性能指标及参数设置 本节通过实验验证本文方法的有效性并与已 有方法进行比较。本文选取精确度 (ACC)、AUC 第 5 期 张英，等：面向自闭症辅助诊断的无监督模糊特征学习新方法 ·885·

·886· 智能系统学报 第14卷 值、敏感度(SEN)、特异性(SPE)作为实验的评价 进行寻优。相关方法的表述及各参数的详细设置 指标，实验的参数设置使用嵌套五折交叉验证来 等信息如表3所示。 表3各算法说明及参数详细设置 Table 3 Details and parameter settings of algorithms 算法名称 算法说明 特征选择 相关参数设置 仅使用低阶功能连接数据 通过简单特征选择，得到 惩罚项系数c={106,10-3,·,105,10} LOFCC 进行SVM分类 200维特征 高斯核带宽g=(106,105,…,10的，101 仅使用低阶功能连接数据 通过简单特征选择，得到 惩罚项系数c={10-6,10-5,…,10,101 HOFCC 进行SVM分类 200维特征 高斯核带宽g=(10-6,105,…,10,10) 惩罚项系数c=(10-6,10-5,…,10的，10) M2sVC购同时使用低阶功能连接数据和高阶对每一个视角分别采用简单特征选高斯核带宽g=10-,10,10,10 功能连接数据进行SVM分类 择，提取得到200维特征 组合系数01={0,0.1,0.2，…，0.9,1， 62={1,0.9,…,0.2,0.1,01 惩罚项系数c=106,10-5,…,10,101 高斯核带宽g=(10-6,105,…,10的，109 UFFM- 本文方法 对每一个视角分别采用简单特征选尺度参数h=(10-6,10-5,…,105,10) M2SVC 择，提取得到200维特征 组合系数81={0,0.1,0.2，·，0.9,1} 02={1,0.9,…,0.2,0.1,01 模糊指数k=(1.2,1.4,…,2.8,3.01 本文采用式(22)(24)计算性能指标，以此进 M2SVC进行分类，从而取得了比仅使用M2SVC 行算法性能评估： 更好的效果，这表明本文所使用的模糊特征学习 TP+TN ACC= (22) 技术能够得到更有鉴别能力的特征，这对自闭症 TP+FN+TN+FP 的诊断具有更优的判别能力。 TP SEN-TP+FN (23) 表4各算法在自闭症分类中的性能比较NYU TN Table 4 The performance comparison of each algorithm in SPE=TN+FP (24) the classification of autism(NYU) 式中：TP和TN分别代表正确分类为自闭症患者 算法 ACC AUC SEN SPE 的人数和正确分类为正常对照组的人数；FP和 LOFCC 0.67570.7619 0.4375 0.8571 FN分别代表正常个体被分类为自闭症患者的人 HOFCC 0.7027 0.72920.6875 0.7143 数以及自闭症患者被分类为正常对照组的人数； M2SVC 0.7322 0.75770.6632 0.8862 灵敏度SEN作为衡量该方法准确预测出自闭症 UFFM-M2SVC 0.78380.82440.6607 0.8231 患者的能力的指标；特异性SPE是衡量该方法准 确预测出健康个体的能力指标。本文实验使用 表5各算法在自闭症分类中的性能比较(USM LIBSVM8工具包中的高斯核SVM。 Table 5 The performance comparison of each algorithm in 3.2分类结果及分析 the classification of autism(USM) 从表4和表5中可以看到，UFFM-M2SVC使 算法 ACC AUC SEN SPE 用无监督模糊特征映射算法UFFM进行特征学 LOFCC 0.69080.72550.58550.7699 习，在此基础上构造分类器M2SVC,可以取得比 HOFCC 0.71840.73570.5361 0.8353 LOFCC、HOFCC及M2SVC等传统分类器更好的 M2SVC 0.73680.73860.37500.3636 性能。 UFFM-M2SVC 0.78950.81820.36250.9364 M2SVC算法联合使用两个视角的数据进行 分类，其分类结果要明显优于基于单视角的分类 3.3 特征数目对实验结果的影响 器LOFCC和HOFCC;而本文方法采用新型的模 本节研究UFFM所提取的特征数目S对 糊特征学习技术进行特征学习，在此基础上使用 M2SVC性能的影响，图3和图4中给出了算法在

值、敏感度 (SEN)、特异性 (SPE) 作为实验的评价 指标，实验的参数设置使用嵌套五折交叉验证来 进行寻优。相关方法的表述及各参数的详细设置 等信息如表 3 所示。 表 3 各算法说明及参数详细设置 Table 3 Details and parameter settings of algorithms 算法名称 算法说明 特征选择 相关参数设置 LOFCC 仅使用低阶功能连接数据 进行SVM分类 通过简单特征选择，得到 200维特征 c= {10−6 ,10−5 ,··· ,105 ,106 } g={10−6 ,10−5 ,··· ,105 ,106 } 惩罚项系数 高斯核带宽 HOFCC 仅使用低阶功能连接数据 进行SVM分类 通过简单特征选择，得到 200维特征 c={10−6 ,10−5 ,··· ,105 ,106 } g={10−6 ,10−5 ,··· ,105 ,106 } 惩罚项系数 高斯核带宽 M2SVC[16] 同时使用低阶功能连接数据和高阶 功能连接数据进行SVM分类 对每一个视角分别采用简单特征选 择，提取得到200维特征 c={10−6 ,10−5 ,··· ,105 ,106 } g={10−6 ,10−5 ,··· ,105 ,106 } θ1= {0,0.1,0.2,···,0.9,1} θ2= {1,0.9,···,0.2,0.1,0} 惩罚项系数 高斯核带宽 组合系数 ， UFFM￾M2SVC 本文方法 对每一个视角分别采用简单特征选 择，提取得到200维特征 c={10−6 ,10−5 ,··· ,105 ,106 } g={10−6 ,10−5 ,··· ,105 ,106 } h={10−6 ,10−5 ,··· ,105 ,106 } θ1= {0,0.1,0.2,···,0.9,1} θ2= {1,0.9,···,0.2,0.1,0} k = {1.2,1.4,··· ,2.8,3.0} 惩罚项系数 高斯核带宽 尺度参数 组合系数 ， 模糊指数 本文采用式 (22)~(24) 计算性能指标，以此进 行算法性能评估： ACC = TP+TN TP+FN+TN+FP (22) SEN = TP TP+FN (23) SPE = TN TN+FP (24) 式中：TP 和 TN 分别代表正确分类为自闭症患者 的人数和正确分类为正常对照组的人数；FP 和 FN 分别代表正常个体被分类为自闭症患者的人 数以及自闭症患者被分类为正常对照组的人数[17] ； 灵敏度 SEN 作为衡量该方法准确预测出自闭症 患者的能力的指标；特异性 SPE 是衡量该方法准 确预测出健康个体的能力指标。本文实验使用 LIBSVM[18] 工具包中的高斯核 SVM。 3.2 分类结果及分析 从表 4 和表 5 中可以看到，UFFM-M2SVC 使 用无监督模糊特征映射算法 UFFM 进行特征学 习，在此基础上构造分类器 M2SVC，可以取得比 LOFCC、HOFCC 及 M2SVC 等传统分类器更好的 性能。 M2SVC 算法联合使用两个视角的数据进行 分类，其分类结果要明显优于基于单视角的分类 器 LOFCC 和 HOFCC；而本文方法采用新型的模 糊特征学习技术进行特征学习，在此基础上使用 M2SVC 进行分类，从而取得了比仅使用 M2SVC 更好的效果，这表明本文所使用的模糊特征学习 技术能够得到更有鉴别能力的特征，这对自闭症 的诊断具有更优的判别能力。 表 4 各算法在自闭症分类中的性能比较 (NYU) Table 4 The performance comparison of each algorithm in the classification of autism(NYU) 算法 ACC AUC SEN SPE LOFCC 0.675 7 0.761 9 0.437 5 0.857 1 HOFCC 0.702 7 0.729 2 0.687 5 0.714 3 M2SVC 0.732 2 0.757 7 0.663 2 0.886 2 UFFM-M2SVC 0.783 8 0.824 4 0.660 7 0.823 1 表 5 各算法在自闭症分类中的性能比较 (USM) Table 5 The performance comparison of each algorithm in the classification of autism(USM) 算法 ACC AUC SEN SPE LOFCC 0.690 8 0.725 5 0.585 5 0.769 9 HOFCC 0.718 4 0.735 7 0.536 1 0.835 3 M2SVC 0.736 8 0.738 6 0.375 0 0.363 6 UFFM-M2SVC 0.789 5 0.818 2 0.362 5 0.936 4 3.3 特征数目对实验结果的影响 本节研 究 UFFM 所提取的特征数 目 S 对 M2SVC 性能的影响，图 3 和图 4 中给出了算法在 ·886· 智 能 系 统 学 报 第 14 卷

第5期 张英，等：面向自闭症辅助诊断的无监督模糊特征学习新方法 ·887· 数据集NYU和USM上S取不同值时M2SVC的 chopharmacology and biological psychiatry,2017,79: 分类结果。从图3和图4可以看出，UFFM- 434-441 M2SVC算法的分类准确率随着S取值的不同而 [2]DICKIE E W.AMEIS S H.SHAHAB S,et al.Personal- 发生变化，图3中最高达到78.38%，图4中最高 ized intrinsic network topography mapping and functional 达到78.95%，均比M2SVC算法提高了约5%。可 connectivity deficits in autism spectrum disorder[J].Biolo- 知，在特征提取数S合理取值的情况下，本文使 gical psychiatry,2018,84(4):278-286. 用UFFM进行无监督模糊特征学习，不仅可以有 [3]徐云，杨健.自闭症早期发现研究进展).中国临床心理 效保留原始有效的特征信息，而且去除了其中包 学杂志，2014,22(6：1023-1027 含的冗余信息，使得自闭症患者和正常对照组之 XU Yun,YANG Jian.The research progress of the early 间的特征差异更有区分性。 recognization of autism[].Chinese journal of clinical psy- chology,2014,22(6:1023-1027. 100 90 [4]BULLMORE E,SPORNS O.Complex brain networks: 80 graph theoretical analysis of structural and functional sys- 70 一一准确率 60 tems[J].Nature reviews neuroscience,2009,10(3): 写类 -。-AUC值 50 186-198 05060708090100110120130140 特征数日S [5]BULLMORE E,SPORNS O.The economy of brain net- work organization[J].Nature reviews neuroscience,2012. 图3特征数目对分类结果的影响NYU) Fig.3 The influence of the number of feature learning on 13(5):336-349 the classification results (NYU) [6]PLIS S M,SUI Jing,LANE T,et al.High-order interac- 100 tions observed in multi-task intrinsic networks are domin- 一。-准确率 90 --AUC值 ant indicators of aberrant brain function in schizophrenia] 80 70 Neurolmage,.2014,102:35-48 60 [7]CHEN Xiaobo,ZHANG Han,GAO Yue,et al.High-order 50 30354045505560657075 resting-state functional connectivity network for MCI clas- 特征数目S sification[J].Human brain mapping,2016,37(9): 图4特征数目对分类结果的影响(USM) 3282-3296. Fig.4 The influence of the number of feature learning on [8]MEUNIER D.ACHARD S,MORCOM A,et al.Age-re- the classification results(USM) lated changes in modular organization of human brain functional networks[J].Neurolmage,2009,44(3): 4结束语 715-723. 自闭症患者的脑区神经活动异于常人。以此 [9]DSOUZA A M.ABIDIN A Z,CHOCKANATHAN U.et 为出发点，本文基于静息态fMRI提出了一种面 al.Mutual connectivity analysis of resting-state functional 向自闭症辅助诊断的无监督模糊特征学习方法 MRI data with local models[J].Neurolmage,2018,178: UFFM,并与M2SVC相结合，从而实现了自闭症 210-223 的辅助诊断。相对于目前传统的分类算法， [10]CORTES C,VAPNIK V.Support-vector networks[J]. Machine learning,1995,20(3):273-297 UFFM-M2SVC算法训练出的分类模型具有优越 [11]YAN Chaogan,ZANG Yufeng.DPARSF:a MATLAB 且稳定的分类性能。实验结果表明，本文所使用 toolbox for "pipeline"data analysis of resting-state 的无监督模糊特征学习方法能够有效地提取出 fMRI[J].Frontiers in systems neuroscience,2010. MRI数据中的重要特征，从而有助于自闭症的临 413:13 床辅助诊断。但是，本文仅基于静息态fMRI数 [12]LAUFS H.DUNCAN J S.Electroencephalography/func- 据进行实验，针对任务态MRI数据的分类效果 tional MRI in human epilepsy:what it currently can and 还未可知，有待进一步的研究。 cannot do[J].Current opinion in neurology,2007,20(4): 参考文献： 417-423 [13]KESHAVARZI A,SARMADIAN F,SHIRI J,et al.Ap- [1]DUAN Xujun,CHEN Heng,HE Changchun,et al.Rest- plication of ANFIS-based subtractive clustering al- ing-state functional under-connectivity within and between gorithm in soil Cation Exchange Capacity estimation us- large-scale cortical networks across three low-frequency ing soil and remotely sensed data[J].Measurement,2017, bands in adolescents with autism[].Progress in neuro-psy- 95:173-180

数据集 NYU 和 USM 上 S 取不同值时 M2SVC 的 分类结果。从 图 3 和 图 4 可以看出， UFFM￾M2SVC 算法的分类准确率随着 S 取值的不同而 发生变化，图 3 中最高达到 78.38%，图 4 中最高 达到 78.95%，均比 M2SVC 算法提高了约 5%。可 知，在特征提取数 S 合理取值的情况下，本文使 用 UFFM 进行无监督模糊特征学习，不仅可以有 效保留原始有效的特征信息，而且去除了其中包 含的冗余信息，使得自闭症患者和正常对照组之 间的特征差异更有区分性。 特征数目S 准确率 40 50 60 70 80 90 100 110 120 130 140 ,AUC 值/% 50 60 70 80 90 100 准确率 AUC值 图 3 特征数目对分类结果的影响 (NYU) Fig. 3 The influence of the number of feature learning on the classification results (NYU) 特征数目S 准确率 30 35 40 45 50 55 60 65 70 75 ,AUC 值/% 50 60 70 80 90 100 准确率 AUC值 图 4 特征数目对分类结果的影响 (USM) Fig. 4 The influence of the number of feature learning on the classification results (USM) 4 结束语 自闭症患者的脑区神经活动异于常人。以此 为出发点，本文基于静息态 fMRI 提出了一种面 向自闭症辅助诊断的无监督模糊特征学习方法 UFFM，并与 M2SVC 相结合，从而实现了自闭症 的辅助诊断。相对于目前传统的分类算法， UFFM-M2SVC 算法训练出的分类模型具有优越 且稳定的分类性能。实验结果表明，本文所使用 的无监督模糊特征学习方法能够有效地提取出 fMRI 数据中的重要特征，从而有助于自闭症的临 床辅助诊断。但是，本文仅基于静息态 fMRI 数 据进行实验，针对任务态 fMRI 数据的分类效果 还未可知，有待进一步的研究。 参考文献： DUAN Xujun, CHEN Heng, HE Changchun, et al. Rest￾ing-state functional under-connectivity within and between large-scale cortical networks across three low-frequency bands in adolescents with autism[J]. Progress in neuro-psy- [1] chopharmacology and biological psychiatry, 2017, 79: 434–441. DICKIE E W, AMEIS S H, SHAHAB S, et al. Personal￾ized intrinsic network topography mapping and functional connectivity deficits in autism spectrum disorder[J]. Biolo￾gical psychiatry, 2018, 84(4): 278–286. [2] 徐云, 杨健. 自闭症早期发现研究进展 [J]. 中国临床心理 学杂志, 2014, 22(6): 1023–1027. XU Yun, YANG Jian. The research progress of the early recognization of autism[J]. Chinese journal of clinical psy￾chology, 2014, 22(6): 1023–1027. [3] BULLMORE E, SPORNS O. Complex brain networks: graph theoretical analysis of structural and functional sys￾tems[J]. Nature reviews neuroscience, 2009, 10(3): 186–198. [4] BULLMORE E, SPORNS O. The economy of brain net￾work organization[J]. Nature reviews neuroscience, 2012, 13(5): 336–349. [5] PLIS S M, SUI Jing, LANE T, et al. High-order interac￾tions observed in multi-task intrinsic networks are domin￾ant indicators of aberrant brain function in schizophrenia[J]. NeuroImage, 2014, 102: 35–48. [6] CHEN Xiaobo, ZHANG Han, GAO Yue, et al. High-order resting-state functional connectivity network for MCI clas￾sification[J]. Human brain mapping, 2016, 37(9): 3282–3296. [7] MEUNIER D, ACHARD S, MORCOM A, et al. Age-re￾lated changes in modular organization of human brain functional networks[J]. NeuroImage, 2009, 44(3): 715–723. [8] DSOUZA A M, ABIDIN A Z, CHOCKANATHAN U, et al. Mutual connectivity analysis of resting-state functional MRI data with local models[J]. NeuroImage, 2018, 178: 210–223. [9] CORTES C, VAPNIK V. Support-vector networks[J]. Machine learning, 1995, 20(3): 273–297. [10] YAN Chaogan, ZANG Yufeng. DPARSF: a MATLAB toolbox for “pipeline” data analysis of resting-state fMRI[J]. Frontiers in systems neuroscience, 2010, 4(13): 13. [11] LAUFS H, DUNCAN J S. Electroencephalography/func￾tional MRI in human epilepsy: what it currently can and cannot do[J]. Current opinion in neurology, 2007, 20(4): 417–423. [12] KESHAVARZI A, SARMADIAN F, SHIRI J, et al. Ap￾plication of ANFIS-based subtractive clustering al￾gorithm in soil Cation Exchange Capacity estimation us￾ing soil and remotely sensed data[J]. Measurement, 2017, 95: 173–180. [13] 第 5 期 张英，等：面向自闭症辅助诊断的无监督模糊特征学习新方法 ·887·

·888· 智能系统学报 第14卷 [14]DAMAYANTI A.Fuzzy learning vector quantization, 作者简介： neural network and fuzzy systems for classification 张英，女，1992年生，硕士研究 fundus eye images with wavelet transformation[C]//Pro- 生，主要研究方向为模式识别与数据 ceedings of the 2nd International Conferences on Inform- 挖掘。 ation Technology,Information Systems and Electrical En- gineering(ICITISEE).Wuhan,China,2017:331-336. [15]DENG Zhaohong,JIANG Yizhang,CHOI K S,et al. Knowledge-leverage-based TSK Fuzzy System modeling[J]. IEEE transactions on neural networks and learning sys- 王骏.男.1978年生，副教授，博 tems,2013,24(8):1200-1212. 土，主要研究方向为智能计算与数据 [16]ZHANG Daoqiang,WANG Yaping,ZHOU Luping,et al. 挖掘。主持国家自然科学基金项目 1项，江苏省自然科学基金项目1项， Multimodal classification of Alzheimer's disease and mild 申请并获得国家发明专利5项。于 cognitive impairment[J].Neurolmage,2011,55(3): 2016年获江苏省高校科研成果自然科 856-867. 学一等奖。发表学术论文50余篇。 [17]POWERS D M W.Evaluation:from precision,recall and F-factor to ROC,informedness,markedness and correla- 鲍国强，男，1992年生，硕士研究 生，主要研究方向为智能计算与模式 tion[J].Journal of machine learning technologies,2011, 识别。 2(1)37-83. [18]CHANG CC,LIN C J.LIBSVM:a library for support vector machines[J].ACM transactions on intelligent sys- tems and technology,2011,2(3):27

DAMAYANTI A. Fuzzy learning vector quantization, neural network and fuzzy systems for classification fundus eye images with wavelet transformation[C]//Pro￾ceedings of the 2nd International Conferences on Inform￾ation Technology, Information Systems and Electrical En￾gineering (ICITISEE). Wuhan, China, 2017: 331−336. [14] DENG Zhaohong, JIANG Yizhang, CHOI K S, et al. Knowledge-leverage-based TSK Fuzzy System modeling[J]. IEEE transactions on neural networks and learning sys￾tems, 2013, 24(8): 1200–1212. [15] ZHANG Daoqiang, WANG Yaping, ZHOU Luping, et al. Multimodal classification of Alzheimer's disease and mild cognitive impairment[J]. NeuroImage, 2011, 55(3): 856–867. [16] POWERS D M W. Evaluation: from precision, recall and F-factor to ROC, informedness, markedness and correla￾tion[J]. Journal of machine learning technologies, 2011, 2(1): 37–83. [17] CHANG C C, LIN C J. LIBSVM: a library for support vector machines[J]. ACM transactions on intelligent sys￾tems and technology, 2011, 2(3): 27. [18] 作者简介： 张英，女，1992 年生，硕士研究 生，主要研究方向为模式识别与数据 挖掘。 王骏，男，1978 年生，副教授，博 士，主要研究方向为智能计算与数据 挖掘。主持国家自然科学基金项目 1 项，江苏省自然科学基金项目 1 项， 申请并获得国家发明专利 5 项。于 2016 年获江苏省高校科研成果自然科 学一等奖。发表学术论文 50 余篇。 鲍国强，男，1992 年生，硕士研究 生，主要研究方向为智能计算与模式 识别。 ·888· 智 能 系 统 学 报 第 14 卷