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

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

Info

Publication number
CN110448273A
CN110448273A CN201910808552.8A CN201910808552A CN110448273A CN 110448273 A CN110448273 A CN 110448273A CN 201910808552 A CN201910808552 A CN 201910808552A CN 110448273 A CN110448273 A CN 110448273A
Authority
CN
China
Prior art keywords
module
clock
prediction
frequency
vector
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.)
Granted
Application number
CN201910808552.8A
Other languages
English (en)
Other versions
CN110448273B (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

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316—Modalities, i.e. specific diagnostic methods
    • A61B5/369—Electroencephalography [EEG]
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/40—Detecting, measuring or recording for evaluating the nervous system
    • A61B5/4076—Diagnosing or monitoring particular conditions of the nervous system
    • A61B5/4094—Diagnosing or monitoring seizure diseases, e.g. epilepsy
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7271—Specific aspects of physiological measurement analysis
    • A61B5/7282—Event detection, e.g. detecting unique waveforms indicative of a medical condition
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/74—Details of notification to user or communication with user or patient ; user input means
    • A61B5/746—Alarms 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)
  • Public Health (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • General Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Surgery (AREA)
  • Physiology (AREA)
  • Neurology (AREA)
  • Psychiatry (AREA)
  • Neurosurgery (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Signal Processing (AREA)
  • Artificial Intelligence (AREA)
  • Psychology (AREA)
  • Complex Calculations (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (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个线型支持向量机的超平面具有互补性，每个线性支持向量机采用下述公式进行预测：
其中，αi、和b分别对应支持向量系数、支持向量和超平面平移系数；Nsv为支持向量机训练完成后支持向量的个数；yi表示每组支持向量的标签0或1；表示脑电提取的小波系数能量子特征向量；表示支持向量机的核函数；yu取值为-1或1，分别表示癫痫不发作和发作；sig代表的是符号函数。
可选的，所述支持向量机群模块中的每个线型支持向量机包括支持向量和测试向量存储器、乘加运算单元、加法器、寄存器和控制模块；预测过程中，测试向量存储器中的特征向量通过乘加运算单元与支持向量完成线性内积运算，内积运算结果复用乘加运算单元与支持向量系数α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＝∑j∑k|dj(k)|2 (2)
小波系数能量子：
为了保证提取足够多的脑电信号特征，且满足低功耗设计的特点，将公式中所有的平方运算转化成取正运算，并不丢失特征信息，即小波系数能量子：
因此，特征提取模块可得出一段脑电数据的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为例。通过计算机软件导出训练模型参数，包括支持向量系数、支持向量和超平面平移系数。
预测模型公式如下：
其中yi、αi、和b分别对应支持向量便签、支持向量系数、支持向量和超平面平移系数。sig代表的是符号函数。
因此本申请根据预测的数学模型实现SVM，实现方法如下所述：
①、在硬件实现之前，首先对SVM进行训练并进行数据处理。
支持向量系数、支持向量和超平面平移系数都是浮点数，需要对数据进行放大取整采取定点方式，在取整过程中会丢失有效位数。因此选取合适的放大倍数对SVM的准确率至关重要。通过多次仿真，在保证准确率的同时，可得出最小放大倍数。
②、图6为本发明中单个SVM模块的结构图，根据图6可知，单个SVM模块包括支持向量和测试向量存储器、乘加运算单元、加法器、寄存器和控制模块，其中支持向量和测试向量存储器可采用rom IP核直接例化实现；乘加运算单元、加法器、寄存器、控制模块可采用硬件描述语言verilog综合实现。
单个SVM模块对于从特征提取模块传输来的1×3的特征向量，首先进行的是特征向量和支持向量的内积运算，由于训练生成的支持向量是一个Nsv×3的稀疏矩阵，其中Nsv是训练阶段得到的支持向量的个数。
因此在时钟频率f2下，共进行Nsv×3个乘法运算；接下来对该结果进行Nsv×2个加法运算。采用串行计算的方式，以复用运算单元，从而减小了电路的面积，使得该结构更适用于低功耗设备。对于矩阵乘法运算可以写成如式(5)：
在时钟频率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元素，因此也进一步减少计算量。
③、对于上述的结果再次进行矩阵运算：
在步骤②完成后，在时钟频率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个结果相加得出最后结果，如下所示：
将最终结果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癫痫预测结果
为验证本申请提供的基于支持向量机的低功耗癫痫预测电路对脑电信号分析处理时的功耗低于现有采用高斯核函数的方法，使用Synopsis综合软件DC计算本申请所实现的电路功耗，对比结果如下表2所示：
表2 本申请预测电路与传统高斯核支持向量机电路功耗对比结果
相比于现有方法，本申请具有如下改进：精简优化癫痫特征提取算法，使用运算度较小的小波系数能量子作为特征。其次，利用多个线型SVM加权机制和阈值判定策略，使得确保本发明所述的电路的预测准确性的前提下，大幅度的减少运算复杂度，符合了低功耗的要求。综上，本发明能更好的适应智能医疗应用的要求。
本发明实施例中的部分步骤，可以利用软件实现，相应的软件程序可以存储在可读取的存储介质中，如光盘或硬盘等。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

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

Cited By (1)

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

Citations (8)

* 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并行计算实现装置
CN104887222A (zh) * 2015-05-11 2015-09-09 重庆大学 可逆化脑电信号分析方法
CN105046273A (zh) * 2015-07-07 2015-11-11 南京邮电大学 基于多尺度样本熵的癫痫脑皮层电图信号分类方法
CN107049239A (zh) * 2016-12-28 2017-08-18 苏州国科康成医疗科技有限公司 基于可穿戴设备的癫痫脑电特征提取方法
CN107095669A (zh) * 2017-05-10 2017-08-29 北京大学 一种癫痫患者脑电信号的处理方法及系统
CN107530012A (zh) * 2015-02-16 2018-01-02 内森·英特拉托 用于脑活动解析的系统与方法
CN109934089A (zh) * 2018-10-31 2019-06-25 北京航空航天大学 基于监督梯度提升器的多级癫痫脑电信号自动识别方法

Patent Citations (8)

* 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并行计算实现装置
CN107530012A (zh) * 2015-02-16 2018-01-02 内森·英特拉托 用于脑活动解析的系统与方法
CN104887222A (zh) * 2015-05-11 2015-09-09 重庆大学 可逆化脑电信号分析方法
CN105046273A (zh) * 2015-07-07 2015-11-11 南京邮电大学 基于多尺度样本熵的癫痫脑皮层电图信号分类方法
CN107049239A (zh) * 2016-12-28 2017-08-18 苏州国科康成医疗科技有限公司 基于可穿戴设备的癫痫脑电特征提取方法
CN107095669A (zh) * 2017-05-10 2017-08-29 北京大学 一种癫痫患者脑电信号的处理方法及系统
CN109934089A (zh) * 2018-10-31 2019-06-25 北京航空航天大学 基于监督梯度提升器的多级癫痫脑电信号自动识别方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
A. TEMKO 等: "EEG-based neonatal seizure detection with Support Vector Machines", 《CLINICAL NEUROPHYSIOLOGY》 *
YUANFA WANG 等: "Hardware design of multiclass SVM classification for epilepsy and epileptic seizure detection", 《IET CIRCUITS, DEVICES & SYSTEMS》 *
花小朋 等: "基于多个支持向量机组合的入侵检测", 《中原工学院学报》 *

Cited By (1)

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

Also Published As

Publication number Publication date
CN110448273B (zh) 2021-03-30

Similar Documents

Publication Publication Date Title
Wang et al. Linking attention-based multiscale CNN with dynamical GCN for driving fatigue detection
Zhao et al. A 13.34 μW event-driven patient-specific ANN cardiac arrhythmia classifier for wearable ECG sensors
Sadiq et al. Motor imagery BCI classification based on multivariate variational mode decomposition
CN112353402A (zh) 心电信号分类模型的训练方法、心电信号分类方法及装置
CN107180155A (zh) 一种基于异构集成模型的疾病预测方法及系统
CN114415842B (zh) 一种基于位点等效增强的脑机接口解码方法及装置
CN108937912A (zh) 一种基于深度神经网络的自动心律失常分析方法
CN106951063A (zh) 数据管理方法和使用缓存的设备
Rongjun et al. Collaborative extreme learning machine with a confidence interval for P2P learning in healthcare
CN117009780A (zh) 一种基于对比学习的时空频域有效通道注意力运动想象脑电解码方法
Nakra et al. Deep neural network with harmony search based optimal feature selection of EEG signals for motor imagery classification
CN108805206A (zh) 一种用于模拟电路故障分类的改进型lssvm建立方法
Weerakody et al. Cyclic gate recurrent neural networks for time series data with missing values
Übeyli Implementing wavelet transform/mixture of experts network for analysis of electrocardiogram beats
CN110477865A (zh) 一种癫痫发作检测装置、终端设备及存储介质
Sanjeet et al. Low-energy real FFT architectures and their applications to seizure prediction from EEG
Tefai et al. ASIC implementation of a pre-trained neural network for ECG feature extraction
Mohammad et al. Tri-model classifiers for EEG based mental task classification: hybrid optimization assisted framework
CN110448273A (zh) 一种基于支持向量机的低功耗癫痫预测电路
Lu et al. Identification of arrhythmia by using a decision tree and gated network fusion model
Azam et al. Classification of COVID-19 symptoms using multilayer perceptron
Huynh A Survey of Machine Learning algorithms in EEG
Malgieri Ontologies, Machine Learning and Deep Learning in Obstetrics
Chaithra et al. A Review of Machine Learning Techniques Used in the Prediction of Heart Disease.
CN107887029A (zh) 疾病预测方法及装置

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