CN104622446A - 一种基于khm聚类算法的心率变异信号优化方法 - Google Patents
一种基于khm聚类算法的心率变异信号优化方法 Download PDFInfo
- Publication number
- CN104622446A CN104622446A CN201510056129.9A CN201510056129A CN104622446A CN 104622446 A CN104622446 A CN 104622446A CN 201510056129 A CN201510056129 A CN 201510056129A CN 104622446 A CN104622446 A CN 104622446A
- Authority
- CN
- China
- Prior art keywords
- heart rate
- rate variability
- khm
- variability signals
- clustering algorithm
- 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.)
- Pending
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/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
- A61B5/02405—Determining heart rate variability
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Cardiology (AREA)
- Biomedical Technology (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Engineering & Computer Science (AREA)
- Physiology (AREA)
- Heart & Thoracic Surgery (AREA)
- Physics & Mathematics (AREA)
- Molecular Biology (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Veterinary Medicine (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
本发明提供一种基于KHM聚类算法的心率变异信号优化方法,包括以下步骤:步骤1,对由非接触式生命体征监护***获得的心跳信号进行峰值检测和周期计算得到心率变异信号;步骤2,对步骤1获得的心率变异信号采用KHM聚类算法对心率变异信号进行优化处理,剔除由呼吸谐波、随机运动等造成的无关干扰项,从而实现心率变异性信号的优化;步骤3,对步骤2中优化后的心率变异信号利用三次样条插值原理形成心率变异信号时域波形图。本发明通过KHM聚类分析算法能够有效剔除心跳信号中的无关干扰项,从而实现心率变异性信号的优化。
Description
技术领域
本发明属于雷达领域,特别设计一种基于KHM聚类算法的心率变异信号优化方法。
背景技术
心率变异是指逐波心动周期之间的时间变异数。正常人的心律也并非绝对整齐,其两次心脏跳动之间会有微小的时间差异,通过测量连续时间的正常的心脏跳动期间的变异数,即可以了解心率变异性。心率变异信号作为一种易获得心跳情况的无创手段,再加上其也能适当反映人体心理状态的特点,已经广为研究者接受和使用。
由于***所获取的心跳信号中,有很大可能性存在呼吸信号高次谐波、身体随机运动以及***自抖动等干扰项。所以尝试利用基于统计学中的聚类分析的方法,选择合适的点,同时去除不良峰值点,以此尽可能对此类误差的影响进行消除。
聚类分析方法的主要依据是“物以类聚”的基本原理,根据数据集合中固有的特点计算数据间的相似性进行划分,即相似性满足一定条件的数据划分到一组,不满足条件的划分到不同的组里。聚类分析不同于普通分类之处在于,聚类分析在分类之前不知道每个样本的类别或者类别的特征。实际应用中,聚类分析通常是作为数据处理的预处理过程,为其后续的分析与处理打下良好基础。
KHM算法(The k-harmonic means algorithm,K调和均值算法)是基于KM算法(The k means algorithm,K均值算法)的改进,也是FCM算法(fuzzy c-means,模糊聚类算法)的特殊情况。
调和平均数(Harmonic Average)定义为对各个单位值取导数,然后其倒数相加平均,之后再做倒数。主要用于仅知道端点,而不知道中间连续变量的情况下。
在传统的K均值聚类分析中,初始值的选取是一个非常重要的问题,不同的初始中心可以出现不同的聚类结果。因此,选择对初始聚类中心不敏感的算法逐渐受到重视。由于调和平均数对初始值不敏感,因此,KHM算法在不断地迭代过程中,对聚类中心不断的进行更新,从而相对传统KM算法有所改进。
发明内容
为了克服现有技术存在的问题,本发明通过KHM聚类分析算法能够有效剔除心跳信号中的无关干扰项,从而实现心率变异性信号的优化。
实现本发明目的的技术解决方案为:
步骤1,对由非接触式生命体征监护***获得的心跳信号进行峰值检测和周期计算得到心率变异信号;
步骤2,对步骤1获得的心率变异信号采用KHM聚类算法对心率变异信号进行优化处理,剔除由呼吸谐波、随机运动等造成的无关干扰项,从而实现心率变异性信号的优化;
步骤3,对步骤2中优化后的心率变异信号利用三次样条插值原理形成心率变异信号时域波形图。
(1)步骤1中心率变异信号的获取的具体方法为:
步骤1.1,对原始心跳信号进行峰值检测,获得心跳信号的各个峰值点;
步骤1.2,利用得到的峰值点进行心跳信号的周期计算,获得心跳信号的周期。
(2)步骤2中KHM聚类算法主要包括如下步骤:
步骤2.1,定义心率变异信号的数据点集合X={x1,x2...xn},聚类中心的集合C={c1,c2,...,cn}定义阈值ε;
步骤2.2,任意取聚类中心cj并且随机初始化聚类中心集合C;
步骤2.3,计算目标函数值k表示数据分类个数,p≥2表示模糊指数,l表示此次迭代次序;
步骤2.4,计算第i个数据点xi到第j个聚类中心cj的隶属度∪(cj/xi)
步骤2.5,计算每个数据点的权重w(xi)
步骤2.6,根据数据点的隶属度∪(cj/xi)和权重w(xi)更新聚类中心公式:
步骤2.7,重新计算KHMl+1值;
步骤2.8,如果||KHMl+1-KHMl||<ε,则计算停止,输出优化后的心率变异性信号;否则转回步骤2.3。
本发明与现有技术相比,具有以下优点:本发明先采取传统方法获取原始心率变异信号,然后利用KHM聚类算法剔除心率变异信号中的无关干扰项,从而实现心率变异性信号的优化。本发明有效可靠,优化效果明显。
下面结合说明书附图对本发明做进一步说明。
附图说明
图1为本发明的实现流程框图。
图2为本发明的由采集数据拟合成的原始心跳信号波形图。
图3为本发明的未经KHM聚类算法处理的心率变异信号。
图4为本发明的未经KHM聚类算法处理的心率变异信号时域波形图。
图5为本发明的经过KHM聚类算法处理的心率变异信号。
图6为本发明的经过KHM聚类算法处理的心率变异信号时域波形图。
具体实施方式
结合图1,一种基于KHM聚类算法的心率变异信号优化方法,方法步骤如下:
步骤1,对由非接触式生命体征监护***获得的心跳信号进行峰值检测和周期计算得到心率变异信号;
步骤2,对步骤1获得的心率变异信号采用KHM聚类算法对心率变异信号进行优化处理,剔除由呼吸谐波、随机运动等造成的无关干扰项,从而实现心率变异性信号的优化;
步骤3,对步骤2中优化后的心率变异信号利用三次样条插值原理形成心率变异信号时域波形图;
步骤1中心率变异信号的获取的具体方法为:
步骤1.1,对原始心跳信号进行峰值检测,获得心跳信号的各个峰值点;
步骤1.2,利用得到的峰值点进行心跳信号的周期计算,获得心跳信号的周期。
步骤2中KHM聚类算法主要包括如下步骤:
步骤2.1,定义X={x1,x2...xn}表示被聚类的数据,在这里即表示之前获得的原始心率变异信号的数据集。∪(cj/xi)定义数据中样本xi隶属于聚类中心cj的隶属度,C={c1,c2,...,cn}表示聚类中心的集合,p表示KHM算法的模糊指数,一般情况下模糊指数p≥2,w(xij)表示权重,即下一次迭代过程中,数据样本对聚类中心重新计算的影响程度。
输入端:心率变异信号的数据集以及分类的个数k;
输出端:满足KHM最小隶属度的矩阵和聚类中心矩阵。
步骤2.2,任意取聚类中心cj并且随机初始化聚类中心集合C,设迭代次数l=1;
步骤2.3,计算目标函数值 k表示数据分类个数,p表示模糊指数,l表示此次迭代次序;
步骤2.4,计算第i个数据点xi到第j个聚类中心cj的隶属度∪(cj/xi),公式如下:
步骤2.5,计算每个数据点的权重w(xi),公式如下:
步骤2.6,根据数据点的隶属度∪(cj/xi)和权重w(xi)更新聚类中心公式:
步骤2.7,l=l+1,重新计算KHMl+1值;
步骤2.8,如果||KHMl+1-KHMl||<ε,则计算停止,其中ε是我们所规定的一个值;否则转回步骤2.3。
通过以上步骤最后就能得到获取所需优化后的心率变异信号。
结合图2、图3、图4、图5、图6:
首先通过非接触式生命体征监护***采集心跳信号数据,然后通过曲线拟合方法拟合心跳信号,图2是为由采集数据拟合成的原始心跳信号波形图;接着对拟合的原始心跳信号进行峰值检测和周期计算得到未经KHM聚类算法处理的心率变异信号即如图3所示;然后对图3的心率变异信号利用三次样条插值原理形成心率变异信号时域波形图即图4;由于信号中存在由呼吸谐波、随机运动等引起的干扰点,为了剔除不必要的干扰点,引进KHM聚类算法进行处理。对图3利用KHM聚类算法进行优化处理得到经过KHM聚类算法处理的心率变异信号即图5,再对图5利用三次样条插值法进行处理得到经过KHM聚类算法处理的心率变异信号时域波形图即图6。
Claims (3)
1.一种基于KHM聚类算法的心率变异信号优化方法,其特征在于,方法步骤如下:
步骤1,对由非接触式生命体征监护***获得的心跳信号进行峰值检测和周期计算得到心率变异信号;
步骤2,对步骤1获得的心率变异信号采用KHM聚类算法对心率变异信号进行优化处理,剔除由呼吸谐波、随机运动等造成的无关干扰项,从而实现心率变异性信号的优化;
步骤3,对步骤2中优化后的心率变异信号利用三次样条插值法形成心率变异信号时域波形图。
2.根据权利要求书1所述的基于KHM聚类算法的心率变异信号优化方法,其特征在于:步骤1中心率变异信号的获取的具体方法为:
步骤1.1,对原始心跳信号进行峰值检测,获得心跳信号的各个峰值点;
步骤1.2,利用得到的峰值点进行心跳信号的周期计算,获得心跳信号的周期。
3.根据权利要求书1所述的基于KHM聚类算法的心率变异信号优化方法,其特征在于,步骤2中KHM聚类算法主要包括如下步骤:
步骤2.1,定义心率变异信号的数据点集合X={x1,x2...xn},聚类中心的集合C={c1,c2,...,cn},定义阈值ε;
步骤2.2,任意取聚类中心cj并且随机初始化聚类中心集合C;
步骤2.3,计算目标函数值k表示数据分类个数,p≥2表示模糊指数,l表示此次迭代次序;
步骤2.4,计算第i个数据点xi到第j个聚类中心cj的隶属度∪(cj/xi)
步骤2.5,计算每个数据点的权重w(xi)
步骤2.6,根据数据点的隶属度∪(cj/xi)和权重w(xi)更新聚类中心公式:
步骤2.7,重新计算KHMl+1值;
步骤2.8,如果||KHMl+1-KHMl||<ε,则计算停止,输出优化后的心率变异性信号;否则转回步骤2.3。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510056129.9A CN104622446A (zh) | 2015-02-03 | 2015-02-03 | 一种基于khm聚类算法的心率变异信号优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510056129.9A CN104622446A (zh) | 2015-02-03 | 2015-02-03 | 一种基于khm聚类算法的心率变异信号优化方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104622446A true CN104622446A (zh) | 2015-05-20 |
Family
ID=53202065
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510056129.9A Pending CN104622446A (zh) | 2015-02-03 | 2015-02-03 | 一种基于khm聚类算法的心率变异信号优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104622446A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107115103A (zh) * | 2016-08-31 | 2017-09-01 | 哈尔滨工业大学深圳研究生院 | 一种基于阵列传感器设计的脉象宽度检测及3d脉图构建方法 |
CN107307846A (zh) * | 2016-04-27 | 2017-11-03 | 南京理工大学 | 非接触式睡眠分期方法 |
CN108601546A (zh) * | 2016-02-04 | 2018-09-28 | 日本电信电话株式会社 | 生物信号处理方法和生物信号处理设备 |
CN109567866A (zh) * | 2018-10-15 | 2019-04-05 | 广东宝莱特医用科技股份有限公司 | 一种胎心率周期变异的处理方法 |
CN110141205A (zh) * | 2019-05-27 | 2019-08-20 | 深圳市是源医学科技有限公司 | 基于hrv技术的抗压数据、疲劳数据的测试方法和装置 |
CN117807551A (zh) * | 2024-02-29 | 2024-04-02 | 深圳市魔样科技有限公司 | 基于智能戒指的心率异常捕获方法及*** |
-
2015
- 2015-02-03 CN CN201510056129.9A patent/CN104622446A/zh active Pending
Non-Patent Citations (2)
Title |
---|
刘伟: "K调和均值聚类分析原理及应用", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
卢山: "基于时频分析的心率变异性研究", 《中国优秀硕士学位论文全文数据库 医药卫生科技辑》 * |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108601546A (zh) * | 2016-02-04 | 2018-09-28 | 日本电信电话株式会社 | 生物信号处理方法和生物信号处理设备 |
US10918302B2 (en) | 2016-02-04 | 2021-02-16 | Nippon Telegraph And Telephone Corporation | Biological signal processing method and biological signal processing apparatus |
CN107307846A (zh) * | 2016-04-27 | 2017-11-03 | 南京理工大学 | 非接触式睡眠分期方法 |
CN107115103A (zh) * | 2016-08-31 | 2017-09-01 | 哈尔滨工业大学深圳研究生院 | 一种基于阵列传感器设计的脉象宽度检测及3d脉图构建方法 |
CN109567866A (zh) * | 2018-10-15 | 2019-04-05 | 广东宝莱特医用科技股份有限公司 | 一种胎心率周期变异的处理方法 |
CN110141205A (zh) * | 2019-05-27 | 2019-08-20 | 深圳市是源医学科技有限公司 | 基于hrv技术的抗压数据、疲劳数据的测试方法和装置 |
CN110141205B (zh) * | 2019-05-27 | 2022-03-11 | 深圳市是源医学科技有限公司 | 基于hrv技术的抗压数据、疲劳数据的测试方法和装置 |
CN117807551A (zh) * | 2024-02-29 | 2024-04-02 | 深圳市魔样科技有限公司 | 基于智能戒指的心率异常捕获方法及*** |
CN117807551B (zh) * | 2024-02-29 | 2024-05-24 | 深圳市魔样科技股份有限公司 | 基于智能戒指的心率异常捕获方法及*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104622446A (zh) | 一种基于khm聚类算法的心率变异信号优化方法 | |
CN110443798B (zh) | 一种基于磁共振图像的自闭症检测方法、装置及*** | |
CN109077715A (zh) | 一种基于单导联的心电信号自动分类方法 | |
CN107822622A (zh) | 基于深度卷积神经网络的心电图诊断方法和*** | |
CN102999905A (zh) | 基于自适应pcnn的眼底图像血管自动检测方法 | |
CN108056773A (zh) | 基于改进变分模态分解的心电信号qrs波群检测算法 | |
CN103985099B (zh) | 一种弥散张量磁共振图像张量域非局部均值去噪方法 | |
CN103989485A (zh) | 基于脑电波的人体疲劳度评价方法 | |
CN109745033A (zh) | 基于时频二维图像和机器学习的动态心电质量评估方法 | |
CN102824166B (zh) | 一种用于脉搏波分类处理的装置 | |
CN108056770A (zh) | 一种基于人工智能的心率检测方法 | |
CN104887224A (zh) | 面向癫痫脑电信号的特征提取与自动识别方法 | |
CN104952210A (zh) | 一种基于决策级数据融合的疲劳驾驶状态检测***和方法 | |
CN103714281A (zh) | 一种基于心电信号的身份识别方法 | |
CN111407315B (zh) | 一种针对超声多普勒胎心信号的质量评估方法 | |
CN112861813B (zh) | 基于复值卷积神经网络的墙后人体行为识别方法 | |
CN113273978B (zh) | 一种基于超宽带雷达的人体呼吸和心跳频率的检测方法 | |
Jin et al. | Classification of normal and abnormal ECG records using lead convolutional neural network and rule inference | |
CN110327055A (zh) | 一种基于高阶谱和卷积神经网络的心冲击信号的分类方法 | |
CN106709497A (zh) | 一种基于pcnn的红外运动弱小目标检测方法 | |
CN112711979A (zh) | 基于生物雷达的慢速随机运动下非接触生命体征监测 | |
CN101822534A (zh) | 基于相空间重构的脉搏波观测方法 | |
CN108771534A (zh) | 一种基于多小波变换融合下的脉搏信号特征提取方法 | |
CN109124620A (zh) | 一种房颤检测方法、装置及设备 | |
CN114052744A (zh) | 基于脉冲神经网络的心电信号分类方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20150520 |
|
WD01 | Invention patent application deemed withdrawn after publication |