基于相关系数与自适应滤波的心率估计.pdf
《基于相关系数与自适应滤波的心率估计.pdf》由会员分享,可在线阅读,更多相关《基于相关系数与自适应滤波的心率估计.pdf(12页珍藏版)》请在咨信网上搜索。
1、Advances in Applied Mathematics 应用数学进展应用数学进展,2023,12(8),3518-3529 Published Online August 2023 in Hans.https:/www.hanspub.org/journal/aam https:/doi.org/10.12677/aam.2023.128350 文章引用文章引用:许时佳.基于相关系数与自适应滤波的心率估计J.应用数学进展,2023,12(8):3518-3529.DOI:10.12677/aam.2023.128350 基于相关系数与基于相关系数与自适应滤波自适应滤波的的心率估计心率估
2、计 许时佳许时佳 长安大学理学院,陕西 西安 收稿日期:2023年7月13日;录用日期:2023年8月3日;发布日期:2023年8月14日 摘摘 要要 近年来,基于光电容积脉搏波近年来,基于光电容积脉搏波(PPG)实时监测心率受到各界的广泛关注,但由于在运动状态下光电容积实时监测心率受到各界的广泛关注,但由于在运动状态下光电容积脉搏波信号容易受到运动伪迹脉搏波信号容易受到运动伪迹(MA)的叠加影响,尤其当运动幅度越大时,脉搏波信号产生的叠加干扰越的叠加影响,尤其当运动幅度越大时,脉搏波信号产生的叠加干扰越严重,使得难以利用严重,使得难以利用PPG信号准确计算出计心率值。为了保证能够从脉搏波信号
3、中准确信号准确计算出计心率值。为了保证能够从脉搏波信号中准确地地估计心率值,估计心率值,为此提出一种去除脉搏波信号中运动伪迹的新算法为此提出一种去除脉搏波信号中运动伪迹的新算法CC-RLS,此算法利用复合归一化含噪脉搏波信号和均,此算法利用复合归一化含噪脉搏波信号和均方根形式的加速度信号之间的相关系数法方根形式的加速度信号之间的相关系数法(CC)初步减少运动伪迹对初步减少运动伪迹对PPG信号的叠加;选取合适的信号作信号的叠加;选取合适的信号作为递归最小二乘为递归最小二乘(RLS)自适应滤波的输入信号再次滤除相关运动伪迹,增加了与心率相关的谱峰峰值,自适应滤波的输入信号再次滤除相关运动伪迹,增加
4、了与心率相关的谱峰峰值,最终得到较为纯净的脉搏波信号,并且使用谱峰追踪算法进行实时心率估计。根据实验结果表明,相关最终得到较为纯净的脉搏波信号,并且使用谱峰追踪算法进行实时心率估计。根据实验结果表明,相关系数法能够快速的消除运动伪迹谱峰峰值,并且在不同的运动状态下,该算法系数法能够快速的消除运动伪迹谱峰峰值,并且在不同的运动状态下,该算法在在10个个PPG数据集上的心数据集上的心率估计误差较小、相比其他算法具有更少计算时间、较高的估计精度、算法稳定性高、鲁棒性强,并带率估计误差较小、相比其他算法具有更少计算时间、较高的估计精度、算法稳定性高、鲁棒性强,并带来较强的实用价值。来较强的实用价值。关
5、键词关键词 脉搏波信号,相关系数,自适应滤波,运动伪迹,参考信号脉搏波信号,相关系数,自适应滤波,运动伪迹,参考信号 Heart Rate Estimation Based on Correlation Coefficient and Adaptive Filtering Shijia Xu School of Science,Changan University,Xian Shaanxi Received:Jul.13th,2023;accepted:Aug.3rd,2023;published:Aug.14th,2023 Abstract In recent years,real-time
6、 monitoring of heart rate based on photocapacitance pulse wave(PPG)许时佳 DOI:10.12677/aam.2023.128350 3519 应用数学进展 has received widespread attention from various fields.However,due to the fact that the photoca-pacitance pulse wave signal is easily affected by the superposition of motion artifacts(MA)du
7、ring exercise,especially when the amplitude of exercise is larger,the superposition interference gen-erated by the pulse wave signal becomes more severe,making it difficult to accurately calculate the heart rate value using PPG signals.In order to ensure accurate estimation of heart rate values from
8、 pulse wave signals,a new algorithm CC-RLS is proposed to remove motion artifacts from pulse wave signals.This algorithm utilizes the correlation coefficient method(CC)between com-posite normalized noisy pulse wave signals and root mean square acceleration signals to prelimi-narily reduce the superp
9、osition of motion artifacts on PPG signals;Selecting a suitable signal as the input signal for recursive least squares(RLS)adaptive filtering,the relevant motion artifacts are filtered again,increasing the spectral peak value related to heart rate,and finally obtaining a relatively pure pulse wave s
10、ignal.The spectral peak tracking algorithm is used for real-time heart rate estimation.According to the experimental results,the correlation coefficient method can quickly eliminate the peak values of motion artifact spectra.Moreover,under different motion states,the algorithm has smaller heart rate
11、 estimation errors on 10 PPG datasets,has less compu-tational time,higher estimation accuracy,high algorithm stability,strong robustness,and brings strong practical value compared to other algorithms.Keywords Photoplethysmography(PPG),Correlation Coefficient,Adaptive Filtering,Motion Artifact(MA),Re
12、ference Signals Copyright 2023 by author(s)and Hans Publishers Inc.This work is licensed under the Creative Commons Attribution International License(CC BY 4.0).http:/creativecommons.org/licenses/by/4.0/1.引言引言 心率是衡量人们身体健康的重要生理参数,监测心率变化可以为人体亚健康状态提供科学依据。监测心率的方法目前主要使用的是心电图信号(ECG)1和光电容积脉搏波信号(PPG)2,心电信号在
13、连续实时监测具有局限性,而光电容积脉搏波信号监测心率方便简易的优点备受青睐,检测设备多为智能手环、手表3等。但由于采集设备与受测者皮肤存在间隙,在运动过程中会使脉搏波波形被运动噪声叠加,心率计算发生偏差。因此,如何在不同 MA 的场景下仍能获得纯净的 PPG 信号是心率估计的一个重要研究方向。现阶段消除 MA 的算法中,自适应滤波4 5具有良好的滤波效果,此方法非常依赖参考信号,参考信号往往决定滤波的效果,自适应滤波通常和其他方法结合使用,如小波阈值去噪、经验模态分解(EMD)6等方法。Wan C 7使用小波阈值与递归最小二乘(RLS)自适应滤波结合起来去除干扰信号,该方法在强 MA 条件下仍
14、能保持较高的估计精度。Arunkumar K R 8采用递归最小二乘(RLS)和归一化最小均方(NLMS)自适应滤波器结合对 PPG 信号进行去噪。Islam M S 9提出一种基于约束递推最小二乘(cRLS)的去噪算法来消除 PPG 信号中的 MA。Yang D 10提出自适应频谱噪声消除(ASNC),用于去除 PPG 信号中的 MA。独立分量分析(ICA)11是一种经典的盲源分离算法。赵爱东12采用经验模态分解(EMD)结合 ICA 去除头部运动伪迹。Roy V 13通过 ICA 和规范相关分析(CCA)以及离散小波变换和平稳小波变换方法的组合,每一种组合方法都应用于脑电波运动伪迹抑制。L
15、iu G 14提出一种基于支持向量机(SVR)和集合经验模态分解(EEMD)方法(SVR-EEMD),该方法有效地解决“末端Open AccessOpen Access许时佳 DOI:10.12677/aam.2023.128350 3520 应用数学进展 效应”。Zhang 15提出一种 TROIKA 的方法,先用奇异值分解将 PPG 信号分解为多个分量,使用加速度信号识别与 MA 相关的分量。去除分量后,对 PPG 信号与剩余分量进行稀疏信号重构,获得干净的 PPG 信号。此后 Zhang 引入一种称为 JOSS 的方法,提供了较低的平均误差。文武16运用奇异谱分析将矫正后的三轴加速度信号
16、分组为不同频率成分的信号,作为三级快速横向递归最小二乘(FTRLS)算法的参考信号自适应消除运动伪迹,保留了重博波信息。Nathan V 17通过预测和更新过程使用粒子滤波,利用概率密度估计 HR 值。该算法通过校正异常值提高了估计精度。周光祥18通过小波变换对 PPG 信号进行去噪。王杰华19通过 EMD 和奇异值分解对 PPG 进行去噪。胡芳凝20提出一种基于快速滑动平均滤波的去噪算法。基于上述算法,提出一种新算法 CC-RLS 消除 PPG 信号中的运动伪迹。首先计算含噪 PPG 信号和加速度信号不同延时下的相关系数,选择相关系数最大的含噪 PPG 信号和加速度信号,构造拟合加速度信号,
17、使含噪 PPG 信号与拟合加速度信号的差与加速度信号的相关系数最小,得到较为纯净的信号。其次使用递归最小二乘(RLS)再次去除运动伪迹。结果表明,该算法取得较好的滤波效果,使用去除运动伪迹的脉搏波信号进行心率估计,最终得到较为准确的心率值。2.运动伪迹消除算法运动伪迹消除算法 检测者在不同程度的运动状态下,运动伪迹的叠加使脉搏波波形变得混乱无序,难以利用脉搏波信号计算心率值,为了消除运动伪迹,下文介绍了两种算法。2.1.信号预处理及相关性分析信号预处理及相关性分析 设()()(),xyzatatat为采集器采集的,x y z三个坐标轴的加速度信号,()()12,y tyt为传感器获得的两通道
18、PPG 信号,将()()12,y tyt平均形成复合 PPG 信号()y t,能够减少随机噪声。将三轴加速度信号利用平方和的均方根形式转化为一维加速度信号()a t,一维加速度的具体表达式为()()()()222xyza tatatat=+;可简化滤波器结构,提高滤波器参数收敛性。人在静止或运动状态下 HR 范围为 30200,单位为 BPM/min,对应的频率范围为 0.53.5 Hz,设计带通滤波器来消除()y t和()a t中超出频率范围内的干扰噪声,滤波器通带为 0.53.5 Hz,并对()y t和()a t进行归一化统一量纲。图 1 为某窗口 PPG 信号时频域图,图 1(a)(b)
19、信号变得更加“光滑”,图 1(c)(d)去除了正常 HR 范围以外的频率谱峰;“*”表示真实心率对应的谱峰,显然运动干扰带来的谱峰占据了信号频谱的主要最高部分。(a)复合 PPG 时域信号 许时佳 DOI:10.12677/aam.2023.128350 3521 应用数学进展 (b)归一化+带通滤波后复合 PPG 时域信号 (c)复合 PPG 频域信号 (d)归一化+带通滤波后复合 PPG 频域信号 Figure 1.Time frequency domain diagram of a certain window PPG signal 图图 1.某窗口 PPG 信号时频域图 许时佳 DOI
20、:10.12677/aam.2023.128350 3522 应用数学进展 加速度()a t对心率影响存在延时性,为了确定延时长度,设延时长度为 T,通过计算()y tT+和()a t的相关系数()R T,如图 2 所示,观察出()R T呈现一定的周期性,周期约为 50。“*”表示40T=,此时相关系数较大,相关性较大,下文提出相关系数法消除运动伪迹。Figure 2.Curve plot of correlation coefficient about T 图图 2.相关系数关于 T 的曲线图 2.2.基于相关系数基于相关系数(CC)的噪声分离的噪声分离 加速度信号()a t造成运动伪迹在
21、PPG 信号的叠加,最大延时40T=,通过将延时后的 PPG 信号减去加速度信号的一个倍数,消除这种叠加的影响,其中 w 为倍数,假设消除叠加影响后的 PPG 信号表示为()()()40Y ty twa t=+;()Y t和()a t的相关系数为()()(),wRY ta t,公式如下:()()()()()()()()()()1223212cov,2wY ta tRY ta tD Y tD a tCwCCw CwCC=+(1)其中1C是()40y t+和()a t的协方差、2C是()a t的方差、3C是()40y t+的方差。Figure 3.Graph of the objective fu
22、nction with respect to w 图图 3.目标函数关于 w 的曲线图 许时佳 DOI:10.12677/aam.2023.128350 3523 应用数学进展 确定叠加系数 w,以相关系数()()(),wRY ta t最小化为目标函数;目标函数越小,则去除运动伪迹的效果越好,目标函数公式如下:1223212min2wCwCCw CwCC+(2)目标函数的图像如图3 所示,经计算,0.27w=时相关系数()()()0.0022,wRY ta t=,此时目标函数最小。()y t与()Y t的频谱图如下,“*”代表真实心率对应的谱峰,“o”代表运动干扰对应的谱峰。(a)原始信号的频
23、谱图 (b)消除叠加影响后信号的频谱图 Figure 4.Spectrum diagram of the original PPG signal and the signal after elimination of superposition 图图 4.原始 PPG 信号与消除叠加后信号的频谱图 运动伪迹对 PPG 信号的干扰应为多个时刻加速度信号叠加的影响,人在静止状态下,频域上的最大谱峰即为真实心率所对应的谱峰21。图 4(a)发现受到严重噪声干扰时,运动干扰的谱峰为信号频谱的最高谱峰,图 4(b)发现经相关系数法后消除了对 PPG 信号影响最大的加速度信号的谱峰,消除了部分运动伪迹但没
24、有完全消除;为了保证心率估算的准确性,考虑用自适应滤波对信号进行进一步去噪。2.3.RLS 自适应滤波去噪自适应滤波去噪 加速度信号作为自适应滤波的参考信号来降低 MA 是目前广泛采用的方法,设()()a twa t=,将()a t许时佳 DOI:10.12677/aam.2023.128350 3524 应用数学进展 作为自适应滤波器的参考信号,()Y t作为自适应滤波器的期望信号。RLS 自适应滤波算法是一种基于最小二乘准则的方法,具有速度收敛和稳定滤波的特点。RLS 自适应滤波系统的框图如图 5 所示。Figure 5.RLS Block diagram of adaptive filt
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 相关系数 自适应 滤波 心率 估计
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。