CN111887834A - 基于多示例学习和进化优化的逐搏心率检测方法 - Google Patents

基于多示例学习和进化优化的逐搏心率检测方法 Download PDF

Info

Publication number
CN111887834A
CN111887834A CN202010682668.4A CN202010682668A CN111887834A CN 111887834 A CN111887834 A CN 111887834A CN 202010682668 A CN202010682668 A CN 202010682668A CN 111887834 A CN111887834 A CN 111887834A
Authority
CN
China
Prior art keywords
signal
heartbeat
template signal
training sample
sample set
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
CN202010682668.4A
Other languages
English (en)
Other versions
CN111887834B (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.)
Xidian University
Original Assignee
Xidian 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 Xidian University filed Critical Xidian University
Priority to CN202010682668.4A priority Critical patent/CN111887834B/zh
Publication of CN111887834A publication Critical patent/CN111887834A/zh
Application granted granted Critical
Publication of CN111887834B publication Critical patent/CN111887834B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
    • A61B5/024—Detecting, measuring or recording pulse rate or heart rate
    • A—HUMAN NECESSITIES
    • A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00—Measuring for diagnostic purposes; Identification of persons
    • A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
    • A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
    • A61B5/1102—Ballistocardiography
    • 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/7225—Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
    • 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/7235—Details of waveform analysis
    • 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/7235—Details of waveform analysis
    • A61B5/7264—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • A61B5/7267—Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Physiology (AREA)
  • Public Health (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Artificial Intelligence (AREA)
  • General Health & Medical Sciences (AREA)
  • Animal Behavior & Ethology (AREA)
  • Signal Processing (AREA)
  • Psychiatry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Cardiology (AREA)
  • Evolutionary Computation (AREA)
  • Fuzzy Systems (AREA)
  • Mathematical Physics (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Power Engineering (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

本发明公开了一种基于多示例学习和进化优化的逐搏心率检测方法，主要解决现有技术对人工标记依赖性高和BCG信号心率估计正确率低的问题。其方案为：采集原始心冲击描图信号和手指电信号；提取心冲击图信号的心跳信号特征，将其先划分为正、负包，再划分为训练集和测试集：对训练样本集进行学习得到初始化心跳模板信号，并对其降维得到降维后的心跳模板信号；对降维后的心跳模板信号进行迭代优化，得到最优心跳模板信号；利用最优心跳模板信号对测试样本集进行分类检测，得到最终心率检测结果。本发明提高了对BCG信号心率的估计正确率，且对心跳特征初始化要求低，人工标记成本低，可用于非精确标记的心冲击描记图信号的心跳检测。

Description

基于多示例学习和进化优化的逐搏心率检测方法
技术领域
本发明属于医学信号处理技术领域，更进一步涉及一种心率估计方法，可用于非精确标记的心冲击描记图信号的心跳检测。
背景技术
心率是评价人们，特别是心血管疾病患者身体健康状况的重要生命体征之一。大多数心脏病患者需要终生服药，即使临床症状消失，心脏问题也可能随时复发。实时心率监测是预防心脏病的必要手段。目前用于测量心率的方法有：心电描记法ECG、光学体积描记法PPG、心音描记法PCG及心冲击描记图法BCG。每种方法都通过测量心搏或心跳周期内人体中发生的不同现象来确定心率。
与ECG相比，BCG以非侵入性和舒适的方式记录了心冲击描记图信号，并且比PPG、PCG精度更高，具有无创性、易操作、传感器和嵌入式***布置等优点，更适合在家中进行长期监测。一种使用心脏冲击成像的液压床可用于在睡眠期间收集人的心跳数据，该液压床有四个传感器，所获得的信号是BCG和呼吸信号的叠加。它有四个液压传感器平行放置在床垫下面，为心率和其他生命体征，如血压和呼吸频率，提供灵活、低成本和稳定的长期监测解决方案。
多示例学习是一种解决非精确标记问题的弱监督方法，是以多示例包为训练单元的学习问题。在多示例学习中，训练集由一组具有分类标签的多示例包组成，每个多示例包含有若干个没有分类标签的示例。如果多示例包至少含有一个正示例，则该包被标记为正类多示例包，即正包。如果多示例包的所有示例都是负示例，则该包被标记为负类多示例包即负包。多示例学习的目的是，通过对具有分类标签的多示例包的学习，建立多示例分类器，并将该分类器应用于未知多示例包的预测。在多示例学习中，多示例包的标记是已知的，示例的标记是未知的。
目前，针对BCG信号心率估计的研究方法中，主要有基于频率分析的方法、基于聚类的方法、基于深度学习网络的方法等。例如：
西安电子科技大学在其申请的专利“一种心冲击图信号的深度回归心率估计方法”(专利申请号：CN201910688377.3，公开号：CN110420019A)中提出了一种基于双向循环神经网络和回归网络的心率估计方法。该方法提出的一种心率估计方法是使用双向循环神经网络获取心跳信号的周期性特征和幅值特征,通过回归网络同时利用心冲击图信号的周期性特征和幅值特征估计心率值,简化了心冲击图信号心率估计的步骤,有效地降低了心冲击图信号心率的估计误差。但是，该方法是有监督的方法，获取准确标注的数据集需要耗费大量人工，在实际应用中受限。
桂林电子科技大学在其申请的专利“一种基于高阶谱和卷积神经网络的心冲击信号的分类方法”(专利申请号：CN201910690101.9，公开号：CN110327055A)中提出了一种基于高阶谱和卷积神经网络的心跳信号分类方法。该方法提出的一种心冲击信号分类方法包括：采用切比雪夫和小波变换的滤波方法对获得的数据集训练样本进行去噪预处理,得到纯净的心冲击信号；对心冲击信号进行高阶谱特征分析,得到信号的幅值和相位的特征信息；构建卷积神经网络模型,将步信号的幅值和相位的特征信息,作为卷积神经网络模型的输入,得到分类结果。该方法利用高阶谱的时移不变性、尺度变化性和相位保持性进行特征提取，不仅可以保留更多的信号信息,还能够抑制高斯有色噪声,从而提升心冲击信号分类的性能；同时,该方法具有更好的泛化性能,并有效解决了高阶谱谱应用的二维模板匹配问题。但是该方法的不足之处是：针对不同数据期望得到同样良好的分类精度，需要进行大量特征提取以及特征选择，这显然大大增加了科研人员的工作量，而且，特征选择的不同会对该方法的分类结果造成一定影响。
发明内容
本发明的目的在于针对上述现有技术的不足，提出一种基于多示例学习和进化优化的心率估计方法，以实现对心冲击描图信号心率的准确估计，且充分利用多示例学习弱监督和多目标优化算法的特性，减少对信号不同特征的选择，降低人工成本，提高BCG信号心跳分类的准确性及心率检测的精确性。
为实现上述目的，本发明技术方案包括如下：
(1)采集原始心冲击描图信号和手指电信号，采样频率为100Hz，并对其进行滤波处理，获得滤波后的心冲击描图信号b和手指电信号f；
(2)提取心冲击图信号b的心跳信号特征fb；
(3)将手指电信号f每个波峰的同时刻的心冲击图信号b的波峰及其左右两个波峰对应的心跳信号特征，组成一个多示例正包，将两个正包间心冲击图信号b其余波峰对应的心跳信号特征记为负包，并将这些正负包按照1:1的比例划分为训练样本集x和测试样本y；
(4)对训练样本集进行学习得到初始化心跳模板信号s：
(4a)计算训练样本集x所有负包中示例的均值μb和方差σb，所有训练样本协方差矩阵
Figure BDA0002586406570000031
并通过对该
Figure BDA0002586406570000032
做特征值分解得到特征向量U和特征值D；
(4b)根据(4a)中的结果，对训练样本集x依次进行白化和归一化处理的预处理，得到预处理后的训练样本集
Figure BDA0002586406570000033
(4c)设心跳模板信号为s，根据(4a)中的结果，对s依次进行白化和归一化的预处理，得到预处理后的心跳模板信号
Figure BDA0002586406570000034
(4d)计算预处理后的心跳模板信号
Figure BDA0002586406570000035
与训练样本集
Figure BDA0002586406570000036
中正包的余弦相似度统计量
Figure BDA0002586406570000037
及
Figure BDA0002586406570000038
与训练样本集
Figure BDA0002586406570000039
中负包的余弦相似度统计量
Figure BDA00025864065700000310
以
Figure BDA00025864065700000311
与
Figure BDA00025864065700000312
中正包示例相似度最大，与负包示例相似度最小为目标方程，将初始化心跳模板信号s'的模值为1作为约束条件，对该目标方程进行求解，得到初始化心跳模板信号s'；
(5)对初始化心跳模板信号s'进行降维处理，得到降维后的心跳模板信号s″；
(6)对降维后的心跳模板信号s″用约束进化算法进行迭代优化，得到最优心跳模板信号ubest；
(7)利用最优心跳模板信号ubest对测试样本集y进行分类检测，得到最终心率检测结果。
本发明与现有技术相比，具有以下优点：
1)人工标注成本低
本发明在训练过程中基于包的概念进行学习，是一种对非精确标记数据的学习方法，不需要精确标注BCG数据信号的每个心跳位置，大大减少了人工标记的成本；
2)心跳分类结果更加精确，心率估计结果更加准确
由于本发明采用进化优化对多示例学习得到的初始化心跳模板信号进行优化，利用多目标进化算法搜索空间大、不易陷入局部最优解的优势，更易找到性能优良的全局最优解，能够对心跳位置与非心跳位置实现准确分类，进而得到更加精确的心率检测结果。
附图说明
图1是本发明的实现流程图。
图2是技术方案3中，多示例正包和负包的概念说明。
具体实施方式
下面结合附图对本发明的实施例和效果作进一步详细描述。
参照图1，本发明的实现步骤如下：
步骤1，采集原始心冲击描图信号和手指电信号。
所述原始心冲击描图信号，是使用v个液压传感器以采样频率fc采集受试者长度为T的v个心冲击图信号；
所述手指电信号，是由指夹式脉冲传感器以采样频率fc采集受试者长度为T的心脏脉冲信号；
同时采集受试者的v个原始心冲击描图信号和手指电信号，为了减弱心冲击图信号中的呼吸分量和高频噪声对心率估计性能的影响，采用截断频率下限为f1、上限为f2的二阶巴特沃斯带通滤波器分别对v个心冲击描图信号及手指电信号进行滤波，得到滤波后的v个心冲击图信号b和手指电信号f，其中f1的取值范围的依据为心冲击图信号中呼吸分量的频率上限约为0.3Hz～0.8Hz，f2的取值范围的依据为心冲击图信号中高频噪声分量的频率下限约为8Hz～12Hz，本实例取但不限于f1＝0.4Hz，f2＝10Hz，v＝4，fc＝100Hz,T＝60000。
步骤2提取心冲击图信号的心跳信号特征。
从检测心冲击描图信号b的所有波峰中，以每个波峰的时间位置为中心点，左右各取45个采样点，组成多个长为91的信号段，即为心跳信号特征fb。
步骤3将心跳信号特征划分正包负包，并将正包负包划分为训练集和测试集。
参照图2，用手指电信号f每个波峰在同时刻心冲击图信号b的波峰及其左右两个波峰对应的心跳信号特征，组成一个多示例正包，将两个正包间心冲击图信号b其余波峰对应的心跳信号特征记为负包；
将这些正负包按照1:1的比例划分为训练样本集x和测试样本y。
步骤4对训练样本集进行学习得到初始化心跳模板信号s。
(4.1)计算训练样本集x所有负包中示例的均值μb和方差σb，公式如下：
Figure BDA0002586406570000041
Figure BDA0002586406570000042
其中，xi为训练样本集中所有负包的样本中的第i个样本，1≤i≤n,n为训练样本集中所有负包的样本总和；
(4.2)计算所有训练样本协方差矩阵
Figure BDA0002586406570000051
表示如下：
Figure BDA0002586406570000052
其中，
Figure BDA0002586406570000053
表示训练样本集x中第p个样本和第q个样本间的协方差，即
Figure BDA0002586406570000054
的第p行第q列元素，1≤p≤m，1≤q≤m，m为训练样本集x中样本总数，xpz为第p个样本中的第z个值，1≤z≤k，k为每个样本中值的总个数，
Figure BDA0002586406570000055
为第p个样本的样本均值，xqz为第q个样本中的第z个值，
Figure BDA0002586406570000056
为第q个样本的样本均值；
(4.3)对该训练样本协方差矩阵
Figure BDA0002586406570000057
做特征值分解得到特征向量U和特征值D；
(4.4)根据(4.1)中的结果，对训练样本集x先进行白化处理，得到白化处理后的
Figure BDA0002586406570000058
再对
Figure BDA0002586406570000059
进行归一化处理，得到预处理后的训练样本集
Figure BDA00025864065700000510
公式表示如下：
Figure BDA00025864065700000511
Figure BDA00025864065700000512
其中，
Figure BDA00025864065700000513
是对特征值D的开方，UT是对特征向量U的转置，
Figure BDA00025864065700000514
是
Figure BDA00025864065700000515
的L1范数；
(4.5)设心跳模板信号为s，根据(4.1)中的结果，对s先进行白化处理，得到白化处理后的
Figure BDA00025864065700000516
再对
Figure BDA00025864065700000517
进行归一化处理，得到预处理后的心跳模板信号
Figure BDA00025864065700000518
公式表示如下：
Figure BDA00025864065700000519
Figure BDA00025864065700000520
其中，
Figure BDA00025864065700000521
是对特征值D的开方，UT是对特征向量U的转置，
Figure BDA00025864065700000522
是
Figure BDA00025864065700000523
的L1范数；
(4.6)计算预处理后的心跳模板信号
Figure BDA00025864065700000524
分别与训练样本集
Figure BDA00025864065700000525
中正包的余弦相似度统计量
Figure BDA0002586406570000061
公式如下：
Figure BDA0002586406570000062
其中，
Figure BDA0002586406570000063
为训练样本集x的第h个正包
Figure BDA0002586406570000064
中与预处理后的心跳模板信号
Figure BDA0002586406570000065
最相似的示例，
Figure BDA0002586406570000066
为第h个正包中的第v个示例，1≤h≤N+，N+为训练样本集
Figure BDA0002586406570000067
中正包数量；
(4.7)计算预处理后的心跳模板信号
Figure BDA0002586406570000068
分别与训练样本集
Figure BDA0002586406570000069
中负包的余弦相似度统计量
Figure BDA00025864065700000610
公式如下：
Figure BDA00025864065700000611
其中，
Figure BDA00025864065700000612
为训练样本集x中第k个负包的第g个示例，1≤k≤N-，N-为负包数量，
Figure BDA00025864065700000613
表示第k个负包内示例数量；
(4.8)为得到初始化心跳模板信号s'，以预处理后的心跳模板信号
Figure BDA00025864065700000614
与训练样本集
Figure BDA00025864065700000615
中正包示例相似度最大，与负包示例相似度最小构建目标方程，同时以初始化心跳模板信号s'的模值为1作为约束条件，该约束的目的是防止在目标方程求解过程中
Figure BDA00025864065700000616
趋于无限大，求解目标方程得到s'，表示如下：
Figure BDA00025864065700000617
其中，
Figure BDA00025864065700000618
为预处理后的心跳模板信号，s'T为对初始化心跳模板信号s'的转置，
Figure BDA00025864065700000619
为预处理后的心跳模板信号
Figure BDA00025864065700000620
与训练样本集
Figure BDA00025864065700000621
中正包的余弦相似度统计量，
Figure BDA00025864065700000622
为预处理后的心跳模板信号
Figure BDA00025864065700000623
与训练样本集
Figure BDA00025864065700000624
中正包的余弦相似度统计量
Figure BDA00025864065700000625
步骤5对初始化心跳模板信号进行降维，得到降维后的心跳模板信号。
由于心跳模板信号维度过大，在进化优化过程中难以收敛至最优解，因此需对初始化心跳模板信号s'进行降维处理，得到降维后的心跳模板信号s″，过程如下：
(5.1)对初始化心跳模板信号s'以采样频率fs做N点傅里叶变换，得到变换后所有频域分量的幅值、频率和相位，本实例中N＝91，fs＝128；
(5.2)将变换后所有频域分量按幅值由大至小排序，取排序后的L维频域分量，得到降维后的心跳模板信号s″：
Figure BDA0002586406570000071
其中al，wl，θl分别为第l维分量对应的幅值，频率，相位，1≤l≤L，2≤L≤12。
步骤6对降维后的心跳模板信号s″用约束进化算法进行迭代优化，得到最优心跳模板信号ubest。
现有约束进化算法有罚函数法，可行解偏好法，多目标优化法。本实例使用但不限于多目标优化法的进化算法，其实现如下：
(6.1)设置最大迭代次数tmax，随机生成种群规模为M＝100心跳模板信号的初始种群At；
(6.2)设心跳模板信号种群At中的个体为u，将(4d)中的目标方程改写为个体的适应度函数:G(u)＝[F(u),p(u)]，其中：
Figure BDA0002586406570000072
式中，
Figure BDA0002586406570000073
为训练样本集x的第h个正包
Figure BDA0002586406570000074
中与预处理后的心跳模板信号
Figure BDA0002586406570000075
最相似的示例，1≤h≤N+，N+为训练样本集
Figure BDA0002586406570000076
中正包数量，
Figure BDA0002586406570000077
为
Figure BDA0002586406570000078
的转置，
Figure BDA0002586406570000079
为训练样本集x中第k个负包的第g个示例，1≤k≤N-，N-为负包数量，
Figure BDA00025864065700000710
表示第k个负包内示例数量，
Figure BDA00025864065700000711
为
Figure BDA00025864065700000712
的转置；
(6.3)通过非一致性变异算子变异操作，得到变异后的种群Bt：
由于传统的变异算子的作用和变异代数没有直接关系,当进化到一定代数以后,由于缺乏局部搜索或全程搜索范围相同,很难获得最优解；非一致性变异算子在进化初期,变异的范围相对较大,随着进化到后期,变异的范围越来越小,从而起到一种对***微调的作用，因此本实例采用非一致性变异算子对初始种群At进行变异操作，得到变异后的种群Bt；
(6.4)通过交叉操作，产生子代种群Ct：
由于单纯形交叉算子SPX基于均匀概率分布产生后代，不需要任何适应度信息，计算复杂度低，因此本实例采用单纯形交叉算子SPX对变异后的种群Bt进行交叉操作，产生子代种群Ct；
(6.5)将初始种群At与子代种群Ct合并，得到合集CAt＝Ct∪At；
(6.6)根据(6.2)中的个体适应度函数G(u)计算CAt种群中个体的适应度；
(6.7)对初始种群与子代种群合集CAt利用进化策略NSES根据CAt种群中个体的适应度进行选择，得到t+1代种群At+1；
(6.8)对当前种群At进行判断：
如果At的种群代数t满足(6.1)中的最大迭代次数tmax，则停止迭代，根据帕累托最优法则从种群At中选出最优心跳模板信号ubest；
否则，将t+1代种群At+1转至(6.3)继续迭代。
步骤7利用最优心跳模板信号ubest对测试样本集y进行分类检测，得到最终心率检测结果。
本发明的效果可通过以下仿真进一步说明。
1.仿真条件
仿真在windows10专业版上进行，CPU基本频率为3.3GHZ*2，利用MATLAB R2019a软件进行仿真。
2.仿真内容
仿真一，用现有的基于希尔伯特变换的心冲击图信号心率估计方法HT对6个受试者进行心率估计，分别获取每个受试者测试集样本的估计心率，并分别计算心率估计误差。
仿真二，用现有的基于能量的方法EN对6个受试者进行心率估计，别获取每个受试者测试集样本的估计心率。
仿真三，用本发明方法分别为每个受试者学***均值。
心率估计误差的计算公式为：
Figure BDA0002586406570000081
其中，P为受试者测试集的样本总数，yi和
Figure BDA0002586406570000082
分别为测试集中第i个样本的真实心率和估计心率，·表示绝对值运算。
上述三个仿真实验每个个体的心率估计误差，如表1。
表1本方法和HT算法、EN算法的心率估计误差对比
受试者 现有HT方法 现有EN方法 本发明方法
1 0.85 14.69 0.19±0.04
2 0.79 1.83 0.66±0.16
3 0.48 1.65 0.60±0.38
4 1.26 0.77 0.27±0.10
5 1.07 0.53 0.087±0.09
6 1.05 0.61 0.32±0.05
3.仿真效果分析
从表1中可以看出，本发明方法在六位受试者的测试集样本上的平均心率估计误差为0.35，而HT方法在六位受试者的测试集样本上的平均心率估计误差为0.92，EN方法在六位受试者的测试集样本上的平均心率估计误差为3.35，显然与HT方法、EN方法相比，本方法的估计误差较小。并且从上表中可以看出本发明的心率估计误差方差更小，说明本发明方法的鲁棒性更强。
以上的实验结果与实验分析可以表明，本发明方法可以从非精确标记的数据中学习得判别性概念，由于进化算法可搜索的解空间到更加广阔，因此不易陷入局部最优，可以得到比其他传统迭代优化方法的更优的解，能够BCG信号心跳位置实现有效的分类，证明了该方法在心跳数据检测上的优越性及鲁棒性。

Claims (10)

1.一种基于多示例学习和进化优化的逐搏心率检测方法，其特征在于，包括如下：
(1)采集原始心冲击描图信号和手指电信号，采样频率为100Hz，并对其进行滤波处理，获得滤波后的心冲击描图信号b和手指电信号f；
(2)提取心冲击图信号b的心跳信号特征fb；
(3)将手指电信号f每个波峰的同时刻的心冲击图信号b的波峰及其左右两个波峰对应的心跳信号特征，组成一个多示例正包，将两个正包间心冲击图信号b其余波峰对应的心跳信号特征记为负包，并将这些正负包按照1:1的比例划分为训练样本集x和测试样本y；
(4)对训练样本集进行学习得到初始化心跳模板信号s：
(4a)计算训练样本集x所有负包中示例的均值μb和方差σb，所有训练样本协方差矩阵
Figure FDA0002586406560000011
并通过对该
Figure FDA0002586406560000012
做特征值分解得到特征向量U和特征值D；
(4b)根据(4a)中的结果，对训练样本集x依次进行白化和归一化处理的预处理，得到预处理后的训练样本集
Figure FDA0002586406560000013
(4c)设心跳模板信号为s，根据(4a)中的结果，对s依次进行白化和归一化的预处理，得到预处理后的心跳模板信号
Figure FDA0002586406560000014
(4d)计算预处理后的心跳模板信号
Figure FDA0002586406560000015
与训练样本集
Figure FDA0002586406560000016
中正包的余弦相似度统计量
Figure FDA0002586406560000017
及
Figure FDA0002586406560000018
与训练样本集
Figure FDA0002586406560000019
中负包的余弦相似度统计量
Figure FDA00025864065600000110
以
Figure FDA00025864065600000111
与
Figure FDA00025864065600000112
中正包示例相似度最大，与负包示例相似度最小为目标方程，将初始化心跳模板信号s'的模值为1作为约束条件，对该目标方程进行求解，得到初始化心跳模板信号s'；
(5)对初始化心跳模板信号s'进行降维处理，得到降维后的心跳模板信号s”；
(6)对降维后的心跳模板信号s”用约束进化算法进行迭代优化，得到最优心跳模板信号ubest；
(7)利用最优心跳模板信号ubest对测试样本集y进行分类检测，得到最终心率检测结果。
2.根据权利要求1所述的方法，其特征在于，(2)中提取心冲击图信号b的心跳信号特征fb，是从检测心冲击描图信号b的所有波峰中，以每个波峰的时间位置为中心，左右各取45个采样点，组成多个长为91的信号段，即为心跳信号特征fb。
3.根据权利要求1所述的方法，其特征在于，(4a)中计算训练样本集x所有负包中示例的均值μb和方差σb，公式如下：
Figure FDA0002586406560000021
Figure FDA0002586406560000022
其中，xi为训练样本集中所有负包的样本中的第i个样本，1≤i≤n,n为训练样本集中所有负包的样本总和。
4.根据权利要求1所述的方法，其特征在于，(4a)中所有训练样本协方差矩阵
Figure FDA0002586406560000023
表示如下：
Figure FDA0002586406560000024
其中，σ(xp,xq)表示
Figure FDA0002586406560000025
为训练样本集x中第p个样本和第q个样本间的协方差，即
Figure FDA0002586406560000026
的第p行第q列元素，1≤p≤m，1≤q≤m，m为训练样本集x中样本总数，xpz为第p个样本中的第z个值，1≤z≤k，k为每个样本中值的总个数，
Figure FDA0002586406560000031
为第p个样本的样本均值，xqz为第q个样本中的第z个值，
Figure FDA0002586406560000032
为第q个样本的样本均值。
5.根据权利要求1所述的方法，其特征在于，(4b)中所述对训练样本集x依次依次进行白化和归一化的预处理，通过如下公式进行：
Figure FDA0002586406560000033
Figure FDA0002586406560000034
其中，
Figure FDA0002586406560000035
是对特征值D的开方，UT是对特征向量U的转置，
Figure FDA00025864065600000324
是对训练样本集x白化处理后的结果，
Figure FDA0002586406560000036
是
Figure FDA0002586406560000037
的L1范数，
Figure FDA0002586406560000038
是对
Figure FDA0002586406560000039
归一化处理后的结果。
6.根据权利要求1所述的方法，其特征在于，(4c)中所述对心跳模板信号为s依次进行白化和归一化的预处理，通过如下公式进行：
Figure FDA00025864065600000310
Figure FDA00025864065600000311
其中，
Figure FDA00025864065600000312
是对特征值D的开方，UT是对特征向量U的转置，
Figure FDA00025864065600000313
是对心跳模板信号s白化处理后的结果，
Figure FDA00025864065600000314
是
Figure FDA00025864065600000315
的L1范数，
Figure FDA00025864065600000316
是对
Figure FDA00025864065600000317
归一化处理后的结果。
7.根据权利要求1所述的方法，其特征在于，(4d)中计算所述预处理后的心跳模板信号
Figure FDA00025864065600000318
分别与训练样本集
Figure FDA00025864065600000319
中正包的余弦相似度统计量
Figure FDA00025864065600000320
和负包的余弦相似度统计量
Figure FDA00025864065600000321
公式如下：
Figure FDA00025864065600000322
Figure FDA00025864065600000323
其中，
Figure FDA0002586406560000041
为训练样本集x的第h个正包
Figure FDA0002586406560000042
中与预处理后的心跳模板信号
Figure FDA0002586406560000043
最相似的示例，
Figure FDA0002586406560000044
为第h个正包中的第v个示例，1≤h≤N+，N+为训练样本集
Figure FDA0002586406560000045
中正包数量，
Figure FDA0002586406560000046
为训练样本集x中第k个负包的第g个示例，1≤k≤N-，N-为负包数量，
Figure FDA0002586406560000047
Figure FDA0002586406560000048
表示第k个负包内示例数量。
8.根据权利要求1所述的方法，其特征在于，(4d)中对目标方程进行求解，得到初始化心跳模板信号s'，公式如下：
Figure FDA0002586406560000049
其中，
Figure FDA00025864065600000410
为预处理后的心跳模板信号，s'T为对初始化心跳模板信号s'的转置，
Figure FDA00025864065600000411
为预处理后的心跳模板信号
Figure FDA00025864065600000412
与训练样本集
Figure FDA00025864065600000413
中正包的余弦相似度统计量，
Figure FDA00025864065600000414
为预处理后的心跳模板信号
Figure FDA00025864065600000415
与训练样本集
Figure FDA00025864065600000416
中正包的余弦相似度统计量
Figure FDA00025864065600000417
9.根据权利要求1所述的方法，其特征在于，(5)中对初始化心跳模板信号s'进行降维处理，实现如下：
(5a)对初始化心跳模板信号s'以采样频率fs做N点傅里叶变换，得到变换后所有频域分量的幅值、频率和相位，其中N＝91，fs＝128；
(5b)将变换后所有频域分量按幅值由大至小排序，取排序后的L维频域分量，得到降维后的心跳模板信号：
Figure FDA00025864065600000418
其中al，wl，θl分别为第l维分量对应的幅值，频率，相位，1≤l≤L，2≤L≤12。
10.根据权利要求1所述的方法，其特征在于，(6)中对降维后的心跳模板信号s”用多目标进化算法进行迭代优化，得到降维后的心跳模板信号s”，实现如下：
(6a)设置最大迭代次数tmax，随机生成种群规模为M＝100心跳模板信号的初始种群At；
(6b)设心跳模板信号种群At中的个体为u，将(4d)中的目标方程改写为个体的适应度函数:G(u)＝[F(u),p(u)]，其中：
Figure FDA0002586406560000051
式中，
Figure FDA0002586406560000052
为训练样本集x的第h个正包
Figure FDA0002586406560000053
中与预处理后的心跳模板信号
Figure FDA0002586406560000054
最相似的示例，1≤h≤N+，N+为训练样本集
Figure FDA0002586406560000055
中正包数量，
Figure FDA0002586406560000056
为
Figure FDA0002586406560000057
的转置，
Figure FDA0002586406560000058
为训练样本集x中第k个负包的第g个示例，1≤k≤N-，N-为负包数量，
Figure FDA0002586406560000059
Figure FDA00025864065600000510
表示第k个负包内示例数量，
Figure FDA00025864065600000511
为
Figure FDA00025864065600000512
的转置；
(6c)采用非一致性变异算子对初始种群At进行变异操作，得到变异后的种群Bt；
(6d)采用单纯形交叉算子SPX对变异后的种群Bt进行交叉操作，产生子代种群Ct；
(6e)将初始种群At与子代种群Ct合并，得到合集CAt＝Ct∪At；
(6f)根据(6b)中的个体适应度函数G(u)计算CAt种群中个体的适应度；
(6g)对初始种群与子代种群合集CAt根据进化策略NSES进行选择，得到t+1代种群At+1；
(6f)对当前种群At进行判断：
如果At的种群代数t满足(6a)中的最大迭代次数tmax，则停止迭代，根据帕累托最优法则从种群At中选出最优心跳模板信号ubest；
否则，将t+1代种群At+1转至(6c)继续迭代。
CN202010682668.4A 2020-07-15 2020-07-15 基于多示例学习和进化优化的逐搏心率检测方法 Active CN111887834B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010682668.4A CN111887834B (zh) 2020-07-15 2020-07-15 基于多示例学习和进化优化的逐搏心率检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010682668.4A CN111887834B (zh) 2020-07-15 2020-07-15 基于多示例学习和进化优化的逐搏心率检测方法

Publications (2)

Publication Number Publication Date
CN111887834A true CN111887834A (zh) 2020-11-06
CN111887834B CN111887834B (zh) 2021-11-02

Family

ID=73193070

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010682668.4A Active CN111887834B (zh) 2020-07-15 2020-07-15 基于多示例学习和进化优化的逐搏心率检测方法

Country Status (1)

Country Link
CN (1) CN111887834B (zh)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106202256A (zh) * 2016-06-29 2016-12-07 西安电子科技大学 基于语义传播及混合多示例学习的Web图像检索方法
CN108154927A (zh) * 2017-12-21 2018-06-12 华子昂 一种智能中医机器人软件架构
US20180285774A1 (en) * 2017-03-31 2018-10-04 Yahoo! Inc. Collaborative personalization via simultaneous embedding of users and their preferences
US10217488B1 (en) * 2017-12-15 2019-02-26 Snap Inc. Spherical video editing
CN109977994A (zh) * 2019-02-02 2019-07-05 浙江工业大学 一种基于多示例主动学习的代表性图像选取方法
CN110420019A (zh) * 2019-07-29 2019-11-08 西安电子科技大学 一种心冲击图信号的深度回归心率估计方法
US10561253B2 (en) * 2016-07-29 2020-02-18 Bryte, Inc. Adaptive sleep system using data analytics and learning techniques to improve individual sleep conditions
CN111241965A (zh) * 2020-01-06 2020-06-05 重庆邮电大学 一种基于多示例学习进行遮挡检测的目标跟踪方法
US20200178887A1 (en) * 2016-04-29 2020-06-11 Fitbit, Inc. Sleep monitoring system with optional alarm functionality

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20200178887A1 (en) * 2016-04-29 2020-06-11 Fitbit, Inc. Sleep monitoring system with optional alarm functionality
CN106202256A (zh) * 2016-06-29 2016-12-07 西安电子科技大学 基于语义传播及混合多示例学习的Web图像检索方法
US10561253B2 (en) * 2016-07-29 2020-02-18 Bryte, Inc. Adaptive sleep system using data analytics and learning techniques to improve individual sleep conditions
US20180285774A1 (en) * 2017-03-31 2018-10-04 Yahoo! Inc. Collaborative personalization via simultaneous embedding of users and their preferences
US10217488B1 (en) * 2017-12-15 2019-02-26 Snap Inc. Spherical video editing
CN108154927A (zh) * 2017-12-21 2018-06-12 华子昂 一种智能中医机器人软件架构
CN109977994A (zh) * 2019-02-02 2019-07-05 浙江工业大学 一种基于多示例主动学习的代表性图像选取方法
CN110420019A (zh) * 2019-07-29 2019-11-08 西安电子科技大学 一种心冲击图信号的深度回归心率估计方法
CN111241965A (zh) * 2020-01-06 2020-06-05 重庆邮电大学 一种基于多示例学习进行遮挡检测的目标跟踪方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
CHANGZHE JIAO ET.AL.: "Heart Beat Characterization from Ballistocardiogram Signals using Extended Functions of Multiple Instances", 《2016 38TH ANNUAL INTERNATIONAL CONFERENCE OF THE IEEE ENGINEERING IN MEDICINE AND BIOLOGY SOCIETY (EMBC)》 *
CHANGZHE JIAO,ET.AL: "Multiple Instance Dictionary Learning for Beat-to-Beat Heart Rate Monitoring from Ballistocardiograms", 《IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING》 *
PELLEGRINI T , CANCES L .: "Cosine-similarity penalty to discriminate sound classes in weakly-supervised sound event detection", 《2019 INTERNATIONAL JOINT CONFERENCE ON NEURAL NETWORKS (IJCNN)》 *
安苏阳: "基于多示例学习的计算机辅助肺结节检测研究", 《中国优秀硕士学位论文全文数据库 (医药卫生科技辑)》 *
王敏: "基于智能床垫的心冲击图信号处理及其在心血管健康评估中的应用研究", 《中国优秀硕士学位论文全文数据库》 *

Also Published As

Publication number Publication date
CN111887834B (zh) 2021-11-02

Similar Documents

Publication Publication Date Title
Atal et al. Arrhythmia classification with ECG signals based on the optimization-enabled deep convolutional neural network
Zhang et al. Automated detection of myocardial infarction using a gramian angular field and principal component analysis network
Sangaiah et al. An intelligent learning approach for improving ECG signal classification and arrhythmia analysis
Li et al. Classification of ECG signals based on 1D convolution neural network
CN111449645B (zh) 心电图心拍智能分类识别方法
Ai et al. Classification of parkinsonian and essential tremor using empirical mode decomposition and support vector machine
Chen et al. Atrial fibrillation detection using a feedforward neural network
Belgacem et al. ECG based human authentication using wavelets and random forests
Li et al. Robust ECG biometrics using GNMF and sparse representation
Wu et al. A novel method to detect multiple arrhythmias based on time-frequency analysis and convolutional neural networks
CN113397555A (zh) 生理参数监测用c-lstm的心率失常分类算法
Zhang et al. ECG signal classification with deep learning for heart disease identification
Altan et al. ECG based human identification using second order difference plots
Kayikcioglu et al. Time-frequency approach to ECG classification of myocardial infarction
Pal et al. Increasing the accuracy of ECG based biometric analysis by data modelling
CN113901893A (zh) 基于多重级联深度神经网络的心电信号的识别与分类方法
Prakash et al. A system for automatic cardiac arrhythmia recognition using electrocardiogram signal
Philip et al. Identifying arrhythmias based on ecg classification using enhanced-PCA and enhanced-SVM methods
Singh et al. Short and noisy electrocardiogram classification based on deep learning
Alam et al. Wearable respiration monitoring: interpretable inference with context and sensor biomarkers
Ganguly et al. A non-invasive approach for fetal arrhythmia detection and classification from ecg signals
Zhang et al. Automated localization of myocardial infarction from vectorcardiographic via tensor decomposition
Allam et al. A deformable CNN architecture for predicting clinical acceptability of ECG signal
Wu et al. ECG identification based on neural networks
Dang et al. A deep biometric recognition and diagnosis network with residual learning for arrhythmia screening using electrocardiogram recordings

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