发明名称 基于深度神经网络的心电图智能诊断方法
摘要 本发明的基于深度神经网络的心电图智能诊断方法,包括:a).信号归一化处理;b).确定训练样本空间;c).确定验证样本空间;d).确定神经网络结构;e).确定激活函数和目标函数;f).训练神经网络;g).心电图信号的自动分析。本发明的基于深度神经网络的心电图智能诊断方法,以MIT-BIH心率失常数据库中的数据为样本,采用Logistic函数作为神经元的激活函数、交叉熵成本函数为目标函数对神经网络进行训练,对心电图信号分析时,利用训练出的神经网络即可得出分析结果,即使诊断医生不具有十分丰富的临床经验,也可获取精确的诊断结果,无需耗费心脏专家大量的精力,减轻医生负担。
申请公布号 CN105105743A 申请公布日期 2015.12.02
申请号 CN201510520104.X 申请日期 2015.08.21
申请人 山东省计算中心(国家超级计算济南中心) 发明人 舒明雷;高岩;王英龙;杨明;王春梅;孔祥龙;王海燕;许继勇;陈长芳;单珂
分类号 A61B5/0402(2006.01)I 主分类号 A61B5/0402(2006.01)I
代理机构 济南泉城专利商标事务所 37218 代理人 褚庆森
主权项 一种基于深度神经网络的心电图智能诊断方法,其特征在于,通过以下步骤来实现:a).信号归一化处理,在MIT‑BIH心率失常数据库中,选取数据库中具有更高信噪比的预处理之后的信号作为训练神经网络使用,利用公式(1)对选取的信号进行归一化处理:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msup><mi>S</mi><mo>&prime;</mo></msup><mo>=</mo><mfrac><mrow><mi>S</mi><mo>-</mo><mi>M</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow></mrow><mrow><mi>M</mi><mi>a</mi><mi>x</mi><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow><mo>-</mo><mi>M</mi><mi>i</mi><mi>n</mi><mrow><mo>(</mo><mi>X</mi><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000786415070000011.GIF" wi="853" he="152" /></maths>其中,S为原始信号,S′为归一化之后的信号,Min(X)为对样本空间进行作用,返回样本空间内数值最小的信号;Max(X)为对样本空间进行作用,返回样本空间内数值最大的信号;b).确定训练样本空间X,心电图信号经过公式(1)的归一化处理后,所有信号数值S∈[0,1],之后将每个信号序列按顺序划分为长度为10000的输入向量x,舍弃长度不足10000的信号组,最后得到满足输入要求的训练样本空间X:X={x<sub>i</sub>|x<sub>i</sub>∈[0,1]<sup>m</sup>,i=1,2,...,m=10000};每个训练样本x<sub>i</sub>对应的心电图异常种类向量y<sub>i</sub>为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>y</mi><mi>i</mi></msub><mo>=</mo><mfenced open = "[" close = "]"><mtable><mtr><mtd><msub><mi>p</mi><mn>1</mn></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>p</mi><mi>j</mi></msub></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><mo>.</mo></mtd></mtr><mtr><mtd><msub><mi>p</mi><mn>10</mn></msub></mtd></mtr></mtable></mfenced><mo>,</mo><mi>j</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>...</mn><mo>,</mo><mn>10</mn><mo>,</mo></mrow>]]></math><img file="FDA0000786415070000012.GIF" wi="517" he="396" /></maths><img file="FDA0000786415070000013.GIF" wi="735" he="160" />其中,10项心电异常种类如下表所示:<img file="FDA0000786415070000014.GIF" wi="1963" he="551" /><img file="FDA0000786415070000021.GIF" wi="1962" he="160" />所有的y<sub>i</sub>就形成了训练样本空间X对应的期望输出空间Y;c).确定验证样本空间,采用与步骤a)和b)相同的方法,确定出检验用的验证样本空间X′和对应的验证样本期望输出空间Y′;d).确定神经网络结构,设L1、Ln分别为神经网络的输入层、输出层,L1至Ln之间的隐藏层分别为L2、L3、…、L(n‑1);L1、L2、…、Ln层中的神经元数目依次递减,并设其分别为n1、n2、…、nn;n1<n2<…<nn;e).确定激活函数和目标函数,选择如公式(3)所示的Logistic函数作为神经元的激活函数;<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>&sigma;</mi><mrow><mo>(</mo><mi>Z</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>+</mo><msup><mi>e</mi><mi>Z</mi></msup></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000786415070000022.GIF" wi="597" he="131" /></maths>选择如公式(4)所示的交叉熵成本函数为目标函数;<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>C</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>w</mi><mo>,</mo><mi>b</mi><mo>)</mo></mrow><mo>=</mo><mo>-</mo><mfrac><mn>1</mn><mi>n</mi></mfrac><munder><mo>&Sigma;</mo><mi>x</mi></munder><mo>&lsqb;</mo><mi>y</mi><mi>ln</mi><mi>a</mi><mo>+</mo><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>y</mi><mo>)</mo></mrow><mi>l</mi><mi>n</mi><mrow><mo>(</mo><mn>1</mn><mo>-</mo><mi>a</mi><mo>)</mo></mrow><mo>&rsqb;</mo><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000786415070000023.GIF" wi="1108" he="141" /></maths>其中,x为输入心电图信号向量,w为神经网络权值,b为神经网络偏置项,n为训练样本的样本容量,y为x对应的心电图异常向量,a为神经网络输出向量;f).训练神经网络,从训练样本空间X中选取n个训练样本,在训练过程中采用随机梯度下降法计算梯度,采用步骤c)中确定出的验证样本空间X′和对应的验证样本期望输出空间Y′进行检验,当正确率超过98%时终止训练,此时得到权值w和偏置项b确定的神经网络;获取如公式(5)所示的心电图异常种类与心电图信号存在的函数关系式:y=Γ(x;w,b)    (5);g).心电图信号的自动分析,假设有一待分析心电图信号,格式化之后的序列为<img file="FDA0000786415070000024.GIF" wi="64" he="59" />将<img file="FDA0000786415070000025.GIF" wi="36" he="51" />带入公式(7)即可得到分析结果<img file="FDA0000786415070000026.GIF" wi="70" he="62" />
地址 250014 山东省济南市历下区科院路19号山东省计算中心