[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN110448273B - 一种基于支持向量机的低功耗癫痫预测电路 - Google Patents

一种基于支持向量机的低功耗癫痫预测电路 Download PDF

Info

Publication number
CN110448273B
CN110448273B CN201910808552.8A CN201910808552A CN110448273B CN 110448273 B CN110448273 B CN 110448273B CN 201910808552 A CN201910808552 A CN 201910808552A CN 110448273 B CN110448273 B CN 110448273B
Authority
CN
China
Prior art keywords
support vector
module
clock
decision
vector machine
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910808552.8A
Other languages
English (en)
Other versions
CN110448273A (zh
Inventor
顾晓峰
田青
虞致国
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jiangnan University
Original Assignee
Jiangnan University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Jiangnan University filed Critical Jiangnan University
Priority to CN201910808552.8A priority Critical patent/CN110448273B/zh
Publication of CN110448273A publication Critical patent/CN110448273A/zh
Application granted granted Critical
Publication of CN110448273B publication Critical patent/CN110448273B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/369Electroencephalography [EEG]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/40Detecting, measuring or recording for evaluating the nervous system
    • A61B5/4076Diagnosing or monitoring particular conditions of the nervous system
    • A61B5/4094Diagnosing or monitoring seizure diseases, e.g. epilepsy
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271Specific aspects of physiological measurement analysis
    • A61B5/7282Event detection, e.g. detecting unique waveforms indicative of a medical condition
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/74Details of notification to user or communication with user or patient ; user input means
    • A61B5/746Alarms related to a physiological condition, e.g. details of setting alarm thresholds or avoiding false alarms

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Molecular Biology (AREA)
  • Animal Behavior & Ethology (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • General Health & Medical Sciences (AREA)
  • Surgery (AREA)
  • Physics & Mathematics (AREA)
  • Neurology (AREA)
  • Physiology (AREA)
  • Neurosurgery (AREA)
  • Psychiatry (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Signal Processing (AREA)
  • Psychology (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Complex Calculations (AREA)

Abstract

本发明公开了一种基于支持向量机的低功耗癫痫预测电路,属于智能医疗应用领域。所述电路在该预测电路中采用线型SVM群模型对脑电特征空间建立超平面并进行分割预测,取代了传统采用高斯核函数SVM,最后通过加权投票机制和阈值判定策略提高线型预测的准确率,发出警报标志,本发明提供的基于支持向量机的癫痫预测电路通过硬件电路实现线性SVM,在保证预测准确率的基础上大大降低了功耗,相对于传统采用高斯核函数实现SVM的方法大大降低了功耗,能够更好的适应智能医疗应用的要求。

Description

一种基于支持向量机的低功耗癫痫预测电路
技术领域
本发明涉及一种基于支持向量机的低功耗癫痫预测电路,属于智能医疗应用领域。
背景技术
针对癫痫发作预测,大量的实验表明,其发作会经历一个过程,它的发作存在很大的预测性。脑电信号作为生理信号的一种,研究者发现可以根据脑电信号的分析实现对癫痫发作的预测,研究者把脑电信号(Electroencephalogram,EEG)分为4个阶段,发作间期、发作前期、发作期和发作后期。癫痫预测的关键在于尽可能早的识别出发作前期的脑电信号。
随着物联网技术和智能医疗技术的发展,目前对生理信号的通用处理方式是:
1)利用传感器将人体生理信号转化为电信号;
2)利用模数转换器将模拟电信号转化为数字信号;
3)利用射频模块将数字信号无线发射至智能终端等监控设备;
4)采用大数据信号处理算法对数字信号进行实时处理。
这种处理方式优点在于,终端可存储着大量生理信号,后期可进行多种算法分析,但是功耗较大。
另一方面,EEG特征和诊断结果之间的关系难以描述,因此早期传统的脑电处理器仅负责采集EEG并上传至云端。近年来,机器学习蓬勃发展,利用机器学习算法模型可以从海量数据中学习出输入到输出的对应关系,从而进行识别和诊断。在关于脑电处理的机器学习算法中,一般使用准确率较高的支持向量机(Support Vetor Machine,SVM)作为算法模型。SVM能够同时最小化经验误差与最大化几何边缘区,其模型依赖于无法观测的隐藏变量。使用SVM算法进行癫病发作预测,前人已做过相关研究。2009年,Netoff等人将6导联颅内EEG的9个不同频段的功率谱作为特征,使用高斯核CSVM(cost-sensitive SVM)识别癫病发作前期和间期的状态,达到了77.8%的敏感度。2014年,Teixeira等人使用多种特征提取如Hjorth统计指数、功率谱、边沿指数等方法提取6导联的头皮EEG特征,将多种特征相结合,用人工神经网络和高斯核SVM进行分类,分类的平均敏感度为73.55%,误警率为0.28/h。2016年,Parvez等人利用6导联EEG的相关相位作为特征向量,使用SVM进行分类,该方法达到了91.95%预测准确度和较低的误警率。2016年,明尼苏达大学张自胜提出一种患者定制的癫痫预测算法模型。模型提取16导联的功率谱特征和交叉相关系数特征,结果显示AdaBoost和SVM的曲线下面积(Area Under Curve,AUC)分别为0.7603和0.8472。
但是上述基于支持向量机的癫痫预测因为采用高斯核函数,所以存在着运算量复杂,功耗大的问题。
发明内容
为了解决目前存在对于癫痫预测存在的运算量复杂、功耗大的问题,本发明提供了一种基于支持向量机的低功耗癫痫预测电路。
一种基于支持向量机的低功耗癫痫预测电路,所述电路包括:时钟生成模块、特征提取模块、支持向量机群模块以及决策模块,所述时钟生成模块分别与特征提取模块、支持向量机群模块及决策模块连接;特征提取模块、支持向量机群模块、决策模块依次连接;
在所述特征提取模块的输入端输入脑电信号,以便所述特征提取模对输入的脑电信号进行特征提取,并将提取到的特征传输至支持向量机群模块;所述支持向量机群模块包括K个线型支持向量机,所述K个线型支持向量机同时根据提取到的特征进行预测,将预测结果传输至决策模块;所述决策模块依次采用加权投票机制和阈值判定策略对K个线型支持向量机的预测结果进行决策得出最终预测结果。
可选的,所述K个线型支持向量机的超平面具有互补性,每个线性支持向量机采用下述公式进行预测:
Figure BDA0002184388060000021
其中,αi
Figure BDA0002184388060000022
和b分别对应支持向量系数、支持向量和超平面平移系数;Nsv为支持向量机训练完成后支持向量的个数;yi表示每组支持向量的标签0或1;
Figure BDA0002184388060000023
表示脑电提取的小波系数能量子特征向量;
Figure BDA0002184388060000024
表示支持向量机的核函数;yu取值为-1或1,分别表示癫痫不发作和发作;sig代表的是符号函数。
可选的,所述支持向量机群模块中的每个线型支持向量机包括支持向量和测试向量存储器、乘加运算单元、加法器、寄存器和控制模块;预测过程中,测试向量存储器中的特征向量
Figure BDA0002184388060000025
通过乘加运算单元与支持向量完成线性内积运算,内积运算结果复用乘加运算单元与支持向量系数αi再次完成内积运算,其结果通过加法器加上超平面平移系数b,然后经过符号判断得出最终预测结果yu放在寄存器中,整个过程中,控制模块控制支持向量和测试向量存储器、乘加运算单元、加法器、寄存器的打开和闭合。
可选的,所述时钟生成模块生成频率为fsamp及2分频时钟f2、4分频时钟f4、8分频时钟f8、16分频时钟f16和32分频时钟f32、预测时钟fp、决策时钟fd;所述时钟生成模块将频率fsamp及2分频时钟f2、4分频时钟f4、8分频时钟f8、16分频时钟f16和32分频时钟f32输入特征提取模块中,将预测时钟fp的时钟接入支持向量机群模块中,将决策时钟fd接入决策模块中。
可选的,所述特征提取模块用于提取脑电信号的频域特征;频域特征为小波分解系数能量子R2、R3、R4;所述特征提取模块提取到的频域特征为由小波分解系数能量子组成的3维特征向量Z=(R2、R3、R4),所述特征提取模块提取到上述3维特征向量Z后将其传输到支持向量机群模块。
可选的,所述支持向量机群模块在接收到所述特征提取模块传输过来的3维特征向量Z后,利用已训练好的预测模型进行预测计算,预测计算过程中采用并行计算方式,在预测时钟fp的时钟下,利用K个线型支持向量机同时预测。
可选的,所述决策模块包括一个存储单元,用于存储K个线型支持向量机的权重βk;在决策时钟fd下,来自支持向量机群模块各个线型支持向量机的预测结果与其对应的权重βk进行乘加运算,从而得出支持向量机群模块的预测结果。
可选的,所述决策模块得出支持向量机群模块的预测结果后,采取阈值判定策略,当连续给出M个预发作状态,则认定脑电信号的癫痫即将发作并发出警报;否则,忽略本次预测结果。
可选的,所述已训练好的预测模型由计算机软件Matlab训练得出,支持向量机群模块实现初步决策;训练过程中导出训练模型,包括支持向量系数、支持向量和超平面平移系数,通过训练模型系数构建预测超平面,实现初步决策功能。
可选的,所述输入的脑电信号为N×1维癫痫患者脑电数据。
本发明有益效果是:
通过提供一种基于支持向量机的低功耗癫痫预测电路,并在该预测电路中采用线型SVM群模型对脑电特征空间建立超平面并进行分割预测,取代了传统采用高斯核函数SVM,最后通过加权投票机制和阈值判定策略提高线型预测的准确率,发出警报标志,本发明提供的基于支持向量机的癫痫预测电路互补线性核SVM的线性组合,在保证预测准确率的基础上大大降低了功耗。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明提供的一种基于SVM的低功耗癫痫预测电路及其方法的模块图。
图2是本发明特征提取模块的电路结构图。
图3是本发明特征提取模块的小波分解运算电路结构图。
图4是本发明中SVM群模块的结构图。
图5是本发明中决策模块的结构图。
图6是本发明中单个线型SVM的结构图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
实施例一:
本实施例提供一种基于SVM的低功耗癫痫预测电路,参见图1,所述预测电路包括:时钟生成模块(1)、特征提取模块(2)、SVM群模块(3)以及决策模块(4)。
所述时钟生成模块(1)分别与特征提取模块(2)、SVM群模块(3)及决策模块(4)连接,特征提取模块(2)、SVM群模块(3)、决策模块(4)依次连接;
使用该基于SVM的低功耗癫痫预测电路进行癫痫发作预测时,将患者脑电信号输入到特征提取模块的输入端,特征提取模块对输入的脑电信号进行特征提取,特征提取模块的输出端与SVM群模块连接,并将处理结果传输到决策模块,决策模块计算最终的预测结果,其中,输入的患者脑电信号为N×1维脑电数据。
预测过程中,时钟生成模块生成频率为fsamp及其2分频时钟f2、4分频时钟f4、8分频时钟f8、16分频时钟f16和32分频时钟f32、预测时钟fp、决策时钟fd
由于特征提取模块内部包含4层小波变换,需要多次降采样获取小波系数,时钟生成模块将频率fsamp的采样时钟及其分频时钟接入特征提取模块中;将预测时钟fp=f2接入SVM群模块中,将决策时钟fd=f2接入决策模块中。
本申请提出的预测电路中的支持向量机群模块包括K个线型支持向量机,K个线型支持向量机同时根据提取到的特征进行预测,故在决策模块中,会出现K个决策结果,其中有K次乘法运算和(K-1)次加法运算。
在每一个线型SVM中,有3×Nsv个内积运算,其中Nsv为支持向量的个数,因此选取频率较高的二分频时钟f2
图2为特征提取模块的结构图,根据图2可知,特征提取模块包括4层小波变换和能量子计算。
特征提取模块对输入的N×1维脑电数据进行特征提取,生成特征向量:Z=(R2、R3、R4),分别表示小波变换的第二、三和四层的小波系数能量子。
常见癫痫脑电数据采样数据为256Hz,对其进行4层小波分解,由小波分解原理可得ca4、cd4、cd3、cd2分别对应0-8Hz、8-16Hz、16-32Hz、32-64Hz,这与脑电信号的节律:delta‘δ’(5-4Hz),theta‘θ’(4-8Hz),alpha‘α’(8-15Hz),Beta‘β’(15-30Hz),gamma‘γ(30Hz以上)基本相吻合。
特征提取过程中,小波分解的过程可以看作是脑电信号和小波滤波器系数的卷积过程。本申请使用的是db2小波,可得出:
低通滤波器系数LD=(-0.1294、0.2241、0.8364、0.4830);
高通滤波器系数HD=(-0.4830、0.8364、-0.2241、-0.1294)。
其中第二层、第三层、第四层只需要对上层分解结果的低频结果作为下层的输入,重复第一层的过程即可。
具体实现过程中,本申请将脑电信号小波分解过程转换为如图3所示的实际运算过程。图3的乘法器_1、乘法器_2、乘法器_3、乘法器_4对应低频滤波系数(-132、229、856、494)以完成低频部分卷积运算,乘法器_5、乘法器_6、乘法器_7、乘法器_8对应高频滤波系数(-494、856、-229、-132)以完成高频部分卷积运算。
(-132、229、856、494)是定点后(210)Db2小波分解低通滤波器系数,(-494、856、-229、-132)是定点后(210)Db2小波分解高通滤波器系数。输入脑电数据分别与小波分解低通和高通滤波器系数进行滑动窗口相乘累加运算,通过分频时钟计算偶数位置的数据,然后对数据进行位数转换,舍弃数据的后10位,得到第一层小波变换的结果。
得到第一层小波变换的结果后,将低频系数再次输入,重复以上过程依次得到第二、第三和第四层结果,此过程复用8个乘法器。
需要注意的是,一方面,EEG信号的指定频率带宽在0.5-100Hz,因此第一层小波变换的细节系数d1(64-128Hz)不作为本申请提取的特征,只需要得到近似系数a1作为下一层的输入即可;另一方面,在第四层小波变换中,只需要提取细节系数d1。因此,可以省去一部分运算,符合低功耗的特点。
在得到小波分解系数d2、d3、d4后,需要计算小波系数能量子;小波能量是将信号进行小波分解后计算出的细节参量的能量值,如果小波基函数是一组正交基函数,那么小波变换具有能量守恒的性质。由于癫痫发作前后的脑电能量变化相应的也会在频谱中有所变换。定义尺度j(即分解层数)下的小波能量为该尺度下细节系数d(k)可由以下公式表示:
Ej=∑k|dj(k)|2 (1)
小波系数总能量Ej为所有尺度的小波系数能量的和:
Et=∑jk|dj(k)|2 (2)
小波系数能量子:
Figure BDA0002184388060000061
为了保证提取足够多的脑电信号特征,且满足低功耗设计的特点,将公式中所有的平方运算转化成取正运算,并不丢失特征信息,即小波系数能量子:
Figure BDA0002184388060000062
因此,特征提取模块可得出一段脑电数据的Z=(R2,R3,R4)小波系数能量子的三维向量。
图4为SVM群模块的结构图,根据图4可知,所述SVM群模块包括K个线型SVM,所述SVM群模块利用已训练好的预测模型将输入的特征向量进行预测计算,采用并行计算方式,在频率f2的时钟下,利用K个线型SVM同时预测。
由于使用的是线型核计算方式,避免了高斯核的e指数运算和乘方运算,极大的降低运算复杂度,因此降低了功耗。
图5为决策模块的结构图,根据图5可知,决策模块包括一个存储模块,用于存储K个线型SVM的权重βi,在f2的时钟频率下,来自SVM群模块(3)的结果与权重βi进行乘加运算,结果存放在寄存器D中。
另外,本申请为消除错误警报采取阈值判定策略,当连续出现M个预测结果,确定患者癫痫将要发作。
SVM群模块前期通过计算机软件仿真,确定K个特异的SVM,所谓模型特异是指SVM的超平面具有互补性。本申请的机器学习算法以低功耗为目标,一方面,线型核运算复杂度最低,因此选取线型SVM。另一方面,机器学习训练部分由计算机软件执行,大幅度地降低了硬件开销。
所述SVM群模块中SVM是最小单元,以一个SVM为例。通过计算机软件导出训练模型参数,包括支持向量系数、支持向量和超平面平移系数。
预测模型公式如下:
Figure BDA0002184388060000063
其中yi、αi
Figure BDA0002184388060000071
和b分别对应支持向量便签、支持向量系数、支持向量和超平面平移系数。sig代表的是符号函数。
因此本申请根据预测的数学模型实现SVM,实现方法如下所述:
①、在硬件实现之前,首先对SVM进行训练并进行数据处理。
支持向量系数、支持向量和超平面平移系数都是浮点数,需要对数据进行放大取整采取定点方式,在取整过程中会丢失有效位数。因此选取合适的放大倍数对SVM的准确率至关重要。通过多次仿真,在保证准确率的同时,可得出最小放大倍数。
②、图6为本发明中单个SVM模块的结构图,根据图6可知,单个SVM模块包括支持向量和测试向量存储器、乘加运算单元、加法器、寄存器和控制模块,其中支持向量和测试向量存储器可采用rom IP核直接例化实现;乘加运算单元、加法器、寄存器、控制模块可采用硬件描述语言verilog综合实现。
单个SVM模块对于从特征提取模块传输来的1×3的特征向量,首先进行的是特征向量和支持向量的内积运算,由于训练生成的支持向量是一个Nsv×3的稀疏矩阵,其中Nsv是训练阶段得到的支持向量的个数。
因此在时钟频率f2下,共进行Nsv×3个乘法运算;接下来对该结果进行Nsv×2个加法运算。采用串行计算的方式,以复用运算单元,从而减小了电路的面积,使得该结构更适用于低功耗设备。对于矩阵乘法运算可以写成如式(5):
Figure BDA0002184388060000072
在时钟频率f2下,当接收到到特征向量后,将计算Z1×SVS(1)1+Z2×SVS(1)2+Z3×SVS(1)3,在第一个时钟沿到来时,首先进行乘法运算Z1×SVS(1)1,运算结果存储到寄存器MEM1,在第二个时钟沿到来时,计算Z2×SVS(1)2+MEM1,以此类推并复用一个乘法单元,在第三个时钟沿计算后,将结果存储到模块4中的MEM1寄存器中,当第六个时钟后计算Z1×SVS(2)1+Z2×SVS(2)2+Z3×SVS(2)3,并将结果分别存储到模块4中的MEM2
以此类推,最终Nsv×3个时钟上升沿之后,MEM1、MEM2…MEMN寄存器中存放计算结果,完成本次所有的特征向量的矩阵运算。另外由于支持向量是一个稀疏矩阵,存在大量0元素,因此也进一步减少计算量。
③、对于上述的结果再次进行矩阵运算:
Figure BDA0002184388060000073
在步骤②完成后,在时钟频率f2下,复用乘法单元。
在第一个时钟沿到来时,计算D1×COE1,运算结果存储到寄存器T中,在第二个时钟沿到来时,计算D2×COE2+T,以此类推并复用乘法单元,在第N个时钟沿计算后,将结果存储到模块4中的寄存器T中。
④、最后实现符号函数sig,判断T+b和0的大小。
如果大于0,标志位Y为-1,反之,标志位Y为1。由于SVM群模块有K个SVM,因此最后传出K个SVM结果。
在使用本申请提供的基于支持向量机的低功耗癫痫预测电路对脑电信号进行分析处理时,包括如下步骤:
步骤一、通过时钟生成模块(1)生成频率为fsamp及其2分频时钟f2、4分频时钟f4、8分频时钟f8、16分频时钟f16和32分频时钟f32
步骤二、在fsamp及其分频时钟频率下,患者脑电信号输入特征提取模块(2)中,特征提取模块(2)生成3维小波系数能量子特征向量Z并传输到SVM机模块(3);
步骤三、通过时钟生成模块,在频率f2的时钟下,所述SVM群模块将输入的特征向量Z利用已建立模型进行预测。所述SVM群模块得出K个结果,再将结果并串转换传输到决策模块中。
步骤四、接收来自SVM群模块(3)的处理结果,采取加权投票机制。由于所训练的SVM是特异的,因此存在互补性。在时钟f2下,读出的每一个SVM的结果都和权重缓存里对应的权重进行相乘,此处权重即是SVM的准确度,最终K个结果相加得出最后结果,如下所示:
Figure BDA0002184388060000081
将最终结果Flag和0比较大小:如果大于0,则单个窗口预测结果将有癫痫发作,反之无癫痫发作。但脑电有时会伴随着噪声,噪声往往只持续几秒钟,因此本申请采取阈值判定策略,当连续给出M个预发作状态,则认定原始脑电信号的癫痫即将发作并发出警报。
为验证本申请提供的基于支持向量机的低功耗癫痫预测电路对脑电信号分析处理的准确性,特实验如下:
使用波士顿儿童医院与麻省理工大学的头皮脑电开源数据库(CHB-MIT ScalpEEG Datebase),采样频率为256Hz,选取1号病人的癫痫数据库,并采用C3-P3单导联,连续截取4s无重叠数据作为一组窗口,即每次输入数据为1024×1维。发作前期定义为距离专家标记发作前35分钟,发作间期定义为距离专家标记发作前至少1小时和发作后至少一小时。预测发作有效范围即为发作前35分钟之内。在训练阶段,发作前期数据来自距离多次专家标记发作前不连续7分钟数据,共取100组。发作间期数据来自距离多次专家标记发作1小时以外的不连续7分钟数据,共取100组。
训练组数据和预测组数据比例为2:1。特征提取模块中fsamp选用20kHz,因为256×4=1024,20k>>1024。在SVM群模块中,本申请经过大量仿真,采用3个线型SVM,即K=3,3个线型SVM分别计为SVM1、SVM2和SVM3,在决策模块中,通过仿真和比对噪声一般不超过15s,因此使用4个窗口可滤除大部分误报并不影响真实情况,即阈值判断策略中M选为4。
为了验证实用性,选取该患者某天从1:44:44到19:24:46连续的脑电信号,该时间段内存在着4次间断的发作,对应数据库中15-21组的数据。
结果显示,电路输入端接入连续的脑电信号后,均能在接下来5个癫痫发作之前发出警报,去噪后平均误警率为0.105/h,预测准确度为100%,详细结果如下表1所示:
表1 患者1:44:44-19:24:46癫痫预测结果
Figure BDA0002184388060000091
为验证本申请提供的基于支持向量机的低功耗癫痫预测电路对脑电信号分析处理时的功耗低于现有采用高斯核函数的方法,使用Synopsis综合软件DC计算本申请所实现的电路功耗,对比结果如下表2所示:
表2 本申请预测电路与传统高斯核支持向量机电路功耗对比结果
Figure BDA0002184388060000101
相比于现有方法,本申请具有如下改进:精简优化癫痫特征提取算法,使用运算度较小的小波系数能量子作为特征。其次,利用多个线型SVM加权机制和阈值判定策略,使得确保本发明所述的电路的预测准确性的前提下,大幅度的减少运算复杂度,符合了低功耗的要求。综上,本发明能更好的适应智能医疗应用的要求。
本发明实施例中的部分步骤,可以利用软件实现,相应的软件程序可以存储在可读取的存储介质中,如光盘或硬盘等。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (9)

1.一种基于支持向量机的低功耗癫痫预测电路,其特征在于,所述电路包括:时钟生成模块、特征提取模块、支持向量机群模块以及决策模块,所述时钟生成模块分别与特征提取模块、支持向量机群模块及决策模块连接;特征提取模块、支持向量机群模块、决策模块依次连接;
在所述特征提取模块的输入端输入脑电信号,以便所述特征提取模块对输入的脑电信号进行特征提取,并将提取到的特征传输至支持向量机群模块;所述支持向量机群模块包括K个线型支持向量机,所述K个线型支持向量机同时根据提取到的特征进行预测,将预测结果传输至决策模块;所述决策模块依次采用加权投票机制和阈值判定策略对K个线型支持向量机的预测结果进行决策得出最终预测结果;
所述K个线型支持向量机的超平面具有互补性,每个线性支持向量机采用下述公式进行预测:
Figure FDA0002581097010000011
其中,αi
Figure FDA0002581097010000012
和b分别对应支持向量系数、支持向量和超平面平移系数;Nsv为支持向量机训练完成后支持向量的个数;yi表示每组支持向量的标签-1或1;
Figure FDA0002581097010000013
表示脑电提取的小波系数能量子特征向量;
Figure FDA0002581097010000014
表示支持向量机的核函数;
yu取值为-1或1,分别表示癫痫不发作和发作。
2.根据权利要求1所述的电路,其特征在于,所述支持向量机群模块中的每个线型支持向量机包括支持向量和测试向量存储器、乘加运算单元、加法器、寄存器和控制模块;
预测过程中,测试向量存储器中的特征向量
Figure FDA0002581097010000015
通过乘加运算单元与支持向量完成内积运算,内积运算结果复用乘加运算单元与支持向量系数αi再次完成内积运算,其结果通过加法器加上超平面平移系数b,然后经过符号函数判断得出最终预测结果yu放在寄存器中,整个过程中,控制模块控制支持向量和测试向量存储器、乘加运算单元、加法器、寄存器的打开和闭合。
3.根据权利要求1所述的电路,其特征在于,所述时钟生成模块生成频率为fsamp及2分频时钟f2、4分频时钟f4、8分频时钟f8、16分频时钟f16和32分频时钟f32、预测时钟fp、决策时钟fd;所述时钟生成模块将频率fsamp及2分频时钟f2、4分频时钟f4、8分频时钟f8、16分频时钟f16和32分频时钟f32输入特征提取模块中,将预测时钟fp的时钟接入支持向量机群模块中,将决策时钟fd接入决策模块中。
4.根据权利要求3所述的电路,其特征在于,所述特征提取模块用于提取脑电信号的频域特征;频域特征为小波分解系数能量子R2、R3、R4;所述特征提取模块提取到的频域特征为由小波分解系数能量子组成的3维特征向量Z=(R2、R3、R4),所述特征提取模块提取到上述3维特征向量Z后将其传输到支持向量机群模块。
5.根据权利要求4所述的电路,其特征在于,所述支持向量机群模块在接收到所述特征提取模块传输过来的3维特征向量Z后,利用已训练好的预测模型进行预测计算,预测计算过程中采用并行计算方式,在预测时钟fp的时钟下,利用K个线型支持向量机同时预测。
6.根据权利要求5所述的电路,其特征在于,所述决策模块包括一个存储单元,用于存储K个线型支持向量机的权重βk;在决策时钟fd下,来自支持向量机群模块各个线型支持向量机的预测结果与其对应的权重βk进行乘加运算,从而得出支持向量机群模块的预测结果。
7.根据权利要求6所述的电路,其特征在于,所述决策模块在得出支持向量机群模块的预测结果后,采取阈值判定策略,当连续给出M个预发作状态,则认定脑电信号的癫痫即将发作并发出警报;否则,忽略本次预测结果。
8.根据权利要求7所述的电路,其特征在于,所述已训练好的预测模型由计算机软件MATLAB训练得出,支持向量机群模块实现初步决策;训练过程中导出训练模型,包括支持向量系数、支持向量和超平面平移系数,通过训练模型系数构建预测超平面,实现初步决策功能。
9.根据权利要求8所述的电路,其特征在于,所述输入的脑电信号为N×1维癫痫患者脑电数据。
CN201910808552.8A 2019-08-29 2019-08-29 一种基于支持向量机的低功耗癫痫预测电路 Active CN110448273B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910808552.8A CN110448273B (zh) 2019-08-29 2019-08-29 一种基于支持向量机的低功耗癫痫预测电路

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910808552.8A CN110448273B (zh) 2019-08-29 2019-08-29 一种基于支持向量机的低功耗癫痫预测电路

Publications (2)

Publication Number Publication Date
CN110448273A CN110448273A (zh) 2019-11-15
CN110448273B true CN110448273B (zh) 2021-03-30

Family

ID=68489975

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910808552.8A Active CN110448273B (zh) 2019-08-29 2019-08-29 一种基于支持向量机的低功耗癫痫预测电路

Country Status (1)

Country Link
CN (1) CN110448273B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113539472B (zh) * 2021-04-01 2024-07-12 西湖大学 生物信号处理器

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101161204A (zh) * 2007-11-20 2008-04-16 深圳市迈科龙影像技术有限公司 超声子宫节育环图像智能识别方法及设备
CN103699697A (zh) * 2014-01-16 2014-04-02 北京航空航天大学 一种支持向量机的fpga并行计算实现装置
CN105046273A (zh) * 2015-07-07 2015-11-11 南京邮电大学 基于多尺度样本熵的癫痫脑皮层电图信号分类方法
CN107095669A (zh) * 2017-05-10 2017-08-29 北京大学 一种癫痫患者脑电信号的处理方法及系统

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2976860C (en) * 2015-02-16 2023-10-17 Nathan Intrator Systems and methods for brain activity interpretation
CN104887222B (zh) * 2015-05-11 2018-02-02 重庆大学 可逆化脑电信号分析方法
CN107049239A (zh) * 2016-12-28 2017-08-18 苏州国科康成医疗科技有限公司 基于可穿戴设备的癫痫脑电特征提取方法
CN109934089B (zh) * 2018-10-31 2020-10-23 北京航空航天大学 基于监督梯度提升器的多级癫痫脑电信号自动识别方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101161204A (zh) * 2007-11-20 2008-04-16 深圳市迈科龙影像技术有限公司 超声子宫节育环图像智能识别方法及设备
CN103699697A (zh) * 2014-01-16 2014-04-02 北京航空航天大学 一种支持向量机的fpga并行计算实现装置
CN105046273A (zh) * 2015-07-07 2015-11-11 南京邮电大学 基于多尺度样本熵的癫痫脑皮层电图信号分类方法
CN107095669A (zh) * 2017-05-10 2017-08-29 北京大学 一种癫痫患者脑电信号的处理方法及系统

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
EEG-based neonatal seizure detection with Support Vector Machines;A. Temko 等;《Clinical Neurophysiology》;20100814;第122卷(第3期);摘要,正文第465页左栏最后一段-第473页左栏第2段 *
Hardware design of multiclass SVM classification for epilepsy and epileptic seizure detection;Yuanfa Wang 等;《IET Circuits, Devices & Systems》;20180118;第12卷(第1期);图2-7 *
基于多个支持向量机组合的入侵检测;花小朋 等;《中原工学院学报》;20071031;第18卷(第5期);全文 *

Also Published As

Publication number Publication date
CN110448273A (zh) 2019-11-15

Similar Documents

Publication Publication Date Title
CN112353402B (zh) 心电信号分类模型的训练方法、心电信号分类方法及装置
CN112244873B (zh) 一种基于混合神经网络的脑电时空特征学习与情感分类方法
CN112001306A (zh) 基于深度卷积对抗生成神经网络的脑电信号解码方法
CN111310570B (zh) 一种基于vmd和wpd的脑电信号情感识别方法及系统
CN104771163A (zh) 基于csp和r-csp算法的脑电信号特征提取方法
CN112800928B (zh) 融合通道和频谱特征的全局自注意力残差网络的癫痫发作预测方法
CN113065526B (zh) 一种基于改进深度残差分组卷积网络的脑电信号分类方法
CN110289081B (zh) 深度网络堆栈模型自适应加权特征融合的癫痫检测方法
CN112200016A (zh) 基于集成学习方法AdaBoost的脑电信号情感识别
CN113011330B (zh) 一种基于多尺度神经网络和空洞卷积的脑电信号分类方法
Li et al. Patient-specific seizure prediction from electroencephalogram signal via multichannel feedback capsule network
Ullah et al. An End‐to‐End Cardiac Arrhythmia Recognition Method with an Effective DenseNet Model on Imbalanced Datasets Using ECG Signal
CN113052099B (zh) 一种基于卷积神经网络的ssvep分类方法
CN113647908B (zh) 波形识别模型的训练、心电波形识别方法、装置及设备
Ziani et al. Extraction of fetal electrocardiogram by combining deep learning and SVD-ICA-NMF methods
Senturk From signal to image: An effective preprocessing to enable deep learning-based classification of ECG
Übeyli Implementing wavelet transform/mixture of experts network for analysis of electrocardiogram beats
CN110448273B (zh) 一种基于支持向量机的低功耗癫痫预测电路
CN114595725B (zh) 一种基于加法网络和监督对比学习的脑电信号分类方法
CN115251909A (zh) 基于时空卷积神经网络的脑电信号评估听力的方法及装置
CN115238796A (zh) 基于并行damscn-lstm的运动想象脑电信号分类方法
Liu et al. Automatic detection of ECG abnormalities by using an ensemble of deep residual networks with attention
CN111603135B (zh) 一种基于主从支持向量机的低功耗癫痫检测电路
CN116898454B (zh) 基于脑电特征融合深度学习模型的癫痫分类方法及系统
CN115017960B (zh) 一种基于时空联合mlp网络的脑电信号分类方法及应用

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant