一种自适应滤波及计算脉搏的系统和方法
【专利摘要】本发明的技术方案包括一种自适应滤波及计算脉搏的系统，包括：采样模块，对传感器输出信号进行采集；变换模块，对采集信号小波分解及重构；计算模块，用于计算信号的RMS值；自适应模块，计算重构信号与总信号的RMS比值，通过对比值大小判断确并计算脉搏。本发明的技术方案包括一种自适应滤波及计算脉搏的放，包括：对传感器输出信号进行采集；对采集信号小波分解及重构；用于计算信号的RMS值；计算重构信号与总信号的RMS比值，通过对比值大小判断确并计算脉搏。本发明的有益效果为：能够自适应滤波以及对脉搏进行精准的计算，能够为后续的分析和计算奠定良好的基础。
【专利说明】
一种自适应滤波及计算脉搏的系统和方法
技术领域
[0001] 本发明涉及一种自适应滤波及计算脉搏的系统和方法，属于计算机医用领域。
【背景技术】
[0002] 在医疗电子领域生理信号处理占据非常重要的位置，各种传感器输出的生理信号 非常微弱并且被噪声淹没，所以获得有用的生理信息提取这些微弱信号是关键的一个环 节，典型的生理信号提取如下：
[0003] (1)用于检测人体含氧量的光电容积波信号。
[0004] (2)用于无创测量血压的脉搏波信号。
[0005] 光电容积波的产生是借光电手段在活体组织中检测血液容积变化的一种无创检 测方法，当一定波长的光束照射到指端皮肤表面时，光束将通过透射或反射方式传送到光 电接收器，心脏收缩时外周血容量最多光吸收量也最大检测到的光强度最小；心脏舒张时 正好相反检测到的光强度最大，将此光强度变化信号转换成电信号便可获得光电容积波。 光电容积波被多种噪声干扰，主要的有环境光、暗电流、工频干扰、运动伪差、呼吸基线漂 移、电磁设备干扰、肌电干扰。常用的数字滤波方法有IIR滤波器，FIR滤波器，卡曼滤波器， 多项式拟合滤波等，自适应能力的滤波器尚未应用。脉搏的计算方法主要有微分法，模式匹 配法，滤波器法。微分法提取突变的部分，与阀值比较，过阀值则为一次脉搏，微分法很容易 被干扰加上阀值不好确定因此效果不太理想，模式匹配法选择匹配的模型比较困难，滤波 法也要选择阀值容易干扰，因此脉搏的准确计算存在困难。
[0006] 无创血压的脉搏波信号的产生流程是先将袖带充气以阻断动脉血流，然后在放气 过程中用压力传感器检测袖带内的气体压力会接收到微弱的压力小脉冲，这种小脉冲就是 脉搏波。脉搏波主要的干扰源有呼吸基线漂移、运动干扰、工频干扰、电磁设备干扰等，滤波 方法、脉搏计算方法与容积波有类似的缺陷。
[0007] 小波变换(wavelettransform，WT)是一种新的变换分析方法，它继承和发展了短 时傅立叶变换局部化的思想，同时又克服了窗口大小不随频率变化等缺点，能够提供一个 随频率改变的"时间-频率"窗口，是进行信号时频分析和处理的理想工具。它的主要特点是 通过变换能够充分突出问题某些方面的特征，能对时间（空间）频率的局部化分析，通过伸 缩平移运算对信号（函数)逐步进行多尺度细化，最终达到高频处时间细分，低频处频率细 分，能自动适应时频信号分析的要求，从而可聚焦到信号的任意细节，解决了Fourier变换 的困难问题，成为继Four i er变换以来在科学方法上的重大突破。
【发明内容】
[0008] 针对现有技术的不足，本发明的技术方案提供了一种自适应滤波及计算脉搏的系 统和方法，用于人体生理信号的自适应滤波及脉搏计算。在人体生理信号处理方面脉搏的 计算是后续分析计算的基础，准确的脉搏计算非常重要，本发明利用小波变换准确的实现 了脉搏的计算及自适应滤波为后续的分析计算奠定了良好的基础。
[0009]本发明的技术方案包括一种自适应滤波及计算脉搏的系统，其特征在于，该系统 包括:采样模块，用于对传感器输出信号进行采样，将采样的信号S存入缓存;变换模块，用 于获取缓存中的采样信号S，对采样信号S进行小波分解，获得近似系数cA6及细节系数cDl、 cD2、cD3、cD4、cD5、cD6，进一步对细节系数cD3、cD4、cD5、cD6分别进行重构得到信号d3、d4、 d5、d6 ;计算模块，用于对信号(13、(14、(15、(16分别计算1^3值，获得对应的(131'1118、(141'1118、 d5rms、d6rms，进一步计算重构总信号SI = d3+d4+d5+d6的RMS，获得Slrms;自适应模块，用 于分别计算d3rms、d4rms、d5rms、d6rms与Slrms的比值，获得四个比值，通过对比值大小判 断确定有用信号滤除噪声完成自适应滤波以及确定脉搏所在的信号，进一步在此信号上过 零法计算脉搏。
[0010]根据所述的自适应滤波及计算脉搏的系统，所述的采样模块还包括:其设置的采 样频率为62.5HZ，采样的信号S的缓存大小为512B。
[0011]根据所述的自适应滤波及计算脉搏的系统，所述的变换模块:对采样信号S进行6 尺度小波分解，其中的基小波为db6。
[0012] 根据所述的自适应滤波及计算脉搏的系统，所述的计算模块:其中RMS的运算公式 为
[0013]
[0014] 根据所述的自适应滤波及计算脉搏的系统，所述的自适应模块:用于计算d6rms与 Slrms的比值，若d6rms与Slrms比值大，则脉搏所在的信号在d6上，d6上使用过零法计算脉 搏，得到有用的信号为d3+d4+d5+d6;若d6rms与Slrms比值小，且d5rms与Slrms比值大，则脉 搏所在的信号在d5上，d5上使用过零法计算脉搏，得到有用的信号为d3+d4+d5;若d6rms与 Slrms比值小，且d5rms与Slrms比值小，同时d4rms与Slrms比值大，则脉搏所在的信号在d4 上，d4上过零法计算脉搏，则有用的信号为d3+d4。
[0015] 本发明的技术方案还包括一种自适应滤波及计算脉搏的方法，该方法包括:对传 感器输出信号进行采样，将采样的信号S存入缓存;获取缓存中的采样信号S，对采样信号S 进行小波分解，获得近似系数cA6及细节系数cDl、cD2、cD3、cD4、cD5、cD6，进一步对细节系 数〇03、〇04、〇05、〇06分别进行重构得到信号(13、(14、(15、(16 ;对信号(13、(14、(15、(16分别计算冊3 值，获得对应的d3rms、d4rms、d5rms、d6rms，进一步计算重构总信号
[0016] 51 = (13+(14+(15+(16的1?]\13，获得311'1118;分别计算(131'1118、(141'1118、(151'1118、(161'1118与311'1118 的比值，获得四个比值，通过对比值大小判断确定有用信号滤除噪声完成自适应滤波以及 确定脉搏所在的信号，进一步在此信号上过零法计算脉搏。
[0017] 进一步，该方法还包括：设置的采样频率为62.5HZ，采样的信号S的缓存大小为 512B。
[0018]进一步，该方法还包括:对采样信号S进行6尺度小波分解，其中的基小波为db6。 [0019]进一步，该方法还包括:RMS的运算公式为
[0020]
[0021 ] 进一步，该方法还包括:计算d6rms与Slrms的比值，若d6rms与Slrms比值大，则脉 搏所在的信号在d6上，d6上使用过零法计算脉搏，得到有用的信号为d3+d4+d5+d6;若d6rms 与Slrms比值小，且d5rms与Slrms比值大，则脉搏所在的信号在d5上，d5上使用过零法计算 脉搏，得到有用的信号为d3+d4+d5 ;若d6rms与Slrms比值小，且d5rms与Slrms比值小，同时 d4rms与SI rms比值大，则脉搏所在的信号在d4上，d4上过零法计算脉搏，则有用的信号为d3 +d4〇
[0022] 本发明的有益效果为：能够自适应滤波以及对脉搏进行精准的计算，能够为后续 的分析和计算奠定良好的基础。
【附图说明】
[0023] 图1所示为根据本发明实施方式的自适应滤波及脉搏计算流程图；
[0024]图2所示为根据本发明实施方式的小波阀值去噪流程图；
[0025]图3a_图3i所示为根据本发明实施方式的基于MATLAB软件实现的图。
【具体实施方式】
[0026] 为了使本发明的目的、技术方案和优点更加清楚，下面结合附图和具体实施例对 本发明进行详细描述。本发明的一种自适应滤波及计算脉搏的系统和方法用于人体生理信 号(血氧容积波信号、无创血压脉搏波信号）的自适应滤波及脉搏计算。
[0027] 图1所示为根据本发明实施方式自适应滤波及脉搏计算流程图。人体生理信号(血 氧容积波、血压脉搏波）的频率主要分布在0.5HZ~8HZ之间，根据奈奎斯特采样定理，当采 样频率fs .max大于信号中最高频率fmax的2倍时（fs .max>2fmax)，采样之后的数字信号完 整地保留了原始信号中的信息，一般实际应用中保证采样频率为信号最高频率的5~10倍， 因此，在本发明中采样频率为62.5HZ。一个长度为N的序列，最多能进行log2N次的分解，采 样的信号S的缓存为512最多可以进行9尺度分解，本发明进行6尺度小波分解，获得近似系 数cA6和细节系数〇01、〇02、〇03、〇04、〇05、〇06，对应的频带依次是〇六6(0~0.53!12)、〇06 (0.53-1.06HZ)、cD5(1.06-2.12HZ)、cD4(2.12-4.24HZ)、cD3(4.24-8.48HZ)、cD2(8.48-16.96HZ)、cDl (16.96-33.92HZ)，人体生理信号(血氧容积波、血压脉搏波）的频带在0.5HZ ~8HZ之间，正好落在了细节系数cD3、cD4、cD5、cD6上，分别重构这些系数获得信号d3、d4、 d5、d6，SI = d3+d4+d5+d6获得人体生理信号SI，细节系数cD2、cDl为高频噪声，近似系数cA6 为呼吸基线干扰。脉搏代表了人体生理信号(血氧容积波、血压脉搏波)基波的频率，人的脉 搏一般40-180次/分，即人体生理信号(血氧容积波、血压脉搏波)基波的频率0.66-3HZ，小 波分解重构后 d6(0.53-1.06HZ)、d5( 1.06-2.12HZ)、d4(2.12-4.24HZ), d6、d5、d4 都有可能 是基波，基波是d4,重构信号d5、d6是没有能量的，基波是d5,重构信号d6是没有能量的，基 波是d6,重构信号d4、d5、d6都有能量，通过计算信号d6、d5、d4的能量与信号Sl能量的比值 及大小的判断就可以确定基波所在的信号计算脉搏，完成自适应滤波，具体的逻辑是 d6rms/Slrms比值大，脉搏所在的信号在d6上，d6上过零法计算脉搏，有用的信号为d3+d4+ d5+d6; d6rms/Slrms比值小，d5rms/Slrms比值大，脉搏所在的信号在d5上，d5上过零法计算 脉搏，有用的信号为d3+d4+d5;d6rms/Slrms比值小，d5rms/Slrms比值小，d4rms/Slrms比值 大，脉搏所在的信号在d4上，d4上过零法计算脉搏，有用的信号为d3+d4。
[0028] 图2所示为根据本发明实施方式小波阀值去噪流程图。对小波阀值去噪的思想做 一下延伸，就会构成自适应滤波器，所谓自适应滤波就是根据不同噪声环境选择不同的滤 波参数，由以上可知，小波阀值去噪的关键是阀值的选择以及所滤系数层的选择，选择的依 据就是能量大小，如果一个系数层的能量占总能量的比值小，可以将此层作为噪声滤除，本 发明应用此原理进行自适应滤波及脉搏的计算。
[0029] 图3a_图3i所示为根据本发明实施方式基于MATLAB软件实现的图。图3a_图3i表示 三组实际采集的血压脉搏波数据为例通过MATLAB软件提供的功能函数和命令来分析实现 小波变换的自适应滤波及脉搏计算方法。
[0030] 以下为MATLAB的程序描述：
[0031] %第一组采样的序列s，长度为512,采样频率62.5HZ。
[0032] plot(s)；
[0033] %获得的彳目号s如图3a所不
[0034] %对8进行6尺度分解，母小波db6
[0035] [c，1] =wavedec(s，6，，db6'）；
[0036] %分解信号s获得的图3b所示
[0037] %分别重构〇03、〇04、〇05、〇06系数
[0038] d3=wrcoef(，d，，c, I，，db6,，3);
[0039] d4=wrcoef(，d，，c, I，，db6,，4);
[0040] d5=wrcoef('d'，c，I，，db6'，5);
[0041 ] d6=wrcoef(，d，，c，I，，db6,，6);
[0042] %分别计算d3、d4、d5、d6的能量与总能量的比值
[0043] sum(d6.'2)/sum((d3+d4+d5+d6).'2)
[0044] ans = 0.6250
[0045] sum(d5.'2)/sum((d3+d4+d5+d6).'2)
[0046] ans = 0.3148
[0047] sum(d4.'2)/sum((d3+d4+d5+d6).'2)
[0048] ans = 0.0565
[0049] % d6能量比较大说明d6是基波，在d6上过零法计算脉搏，最终得到的信号d3+d4+ d5+d6
[0050] subplot(2,1,1)
[0051] plot(s)；
[0052] subplot(2,1,2)
[0053] plot(d6)；
[0054] %获得的信号s和信号d6如图3c所示
[0055] %8与(16有相似的频率，在d6上过零点就可以计算脉搏，最终得到的信号d3+d4+d5 +d6
[0056] 第二组采集序列s I，长度为512，采样频率62.5HZ。
[0057] pl〇t(sl)
[0058] %获得的信号si如图3d所示
[0059] %对81进行6尺度分解，母小波db6
[0060] [c，1] = wavedec(sl，6, 'db6'）；
[0061] %分解信号Si获得的信号如图3e所示
[0062] %分别重构 〇03、〇04、〇05、〇〕6系数
[0063] d3=wrcoef ('d'，c，I，'db6'，3);
[0064] d4=wrcoef ('d'，c，I，'db6'，4);
[0065] d5=wrcoef ('d'，c，I，'db6'，5);
[0066] d6=wrcoef ('d'，c，I，'db6'，6);
[0067] %分别计算(13、(14、(15、(16的能量总能量的比值
[0068] sum(d6 厂2)/sum((d3+d4+d5+d6)厂2)
[0069] ans = 0.0084
[0070] sum(d5.，'2)/sum( (d3+d4+d5+d6). "2)
[0071] ans = 0.6056
[0072] sum(d4.，'2)/sum( (d3+d4+d5+d6) / 2)
[0073] ans = 0.3151
[0074] %d6能量小d5能量大说明d5是基波，在d5上过零法计算脉搏，最终得到的信号d3+ d4+d5
[0075] subplot(2,I，I)
[0076] plot(sl);
[0077] subplot(2,I，2)
[0078] pl〇t(d5)；
[0079] %获得的信号si与d5如图3f所示
[0080] %sl与d5有相似的频率，在d5上过零点就可以计算脉搏，最终得到的信号d3+d4+ d5
[0081] 第三组采集序列s2，长度为512，采样频率62.5HZ。
[0082] plot(s2)
[0083] %得到信号s2,如图3g所示
[0084] %对82进行6尺度分解，母小波db6
[0085] [c，1] = wavedec(s2,6, 'db6'）；
[0086] %S2分解获得的信号如图3h所示
[0087] %分别重构 〇03、〇04、〇05、〇06系数
[0088] d3=wrcoef ('d'，c，I，'db6'，3);
[0089] d4=wrcoef ('d'，c，I，'db6'，4);
[0090] d5=wrcoef ('d'，c，I，'db6'，5);
[0091 ] d6=wrcoef ('d'，c，I，'db6'，6);
[0092] %分别计算(13、(14、(15、(16的能量总能量的比值
[0093] sum(d6厂2)/sum((d3+d4+d5+d6)厂2)
[0094] ans = 0.0034
[0095] sum(d5.，'2)/sum( (d3+d4+d5+d6) / 2)
[0096] ans = 0.0218
[0097] sum(d4.'2)/sum((d3+d4+d5+d6).'2)
[0098] ans = 0.6387
[0099] %d6能量小d5能量小d4能量大说明d4是基波，在d4上过零点法计算脉搏，最终得 到的信号d3+d4
[0100] subplot(2,1,1)
[0101] plot(s2)；
[0102] subplot(2,1,2)
[0103] plot(d4)；
[0104] %获得的信号s2与d4如图3i所示
[0105] %s2与d4有相似的频率，在d4上过零点就可以计算脉搏，最终得到的信号d3+d4
[0106] 以上所述，只是本发明的较佳实施例而已，本发明并不局限于上述实施方式，只要 其以相同的手段达到本发明的技术效果，都应属于本发明的保护范围。在本发明的保护范 围内其技术方案和/或实施方式可以有各种不同的修改和变化。
【主权项】
1. 一种自适应滤波及计算脉搏的系统，其特征在于，该系统包括： 采样模块，用于对传感器输出信号进行采样，将采样的信号S存入缓存； 变换模块，用于获取缓存中的采样信号S，对采样信号S进行小波分解，获得近似系数 cA6及细节系数cDl、cD2、cD3、cD4、cD5、cD6，进一步对细节系数cD3、cD4、cD5、cD6分别进行 重构得到信号d3、d4、d5、d6; 计算模块，用于对信号(13、(14、(15、(16分别计算1^8值，获得对应的(131'1113、(141'1113、(151'1113、 dGrms，进一步计算重构总信号S1 =d3+d4+d5+d6的RMS，获得Slrms; 自适应模块，用于分别计算(131'1113、(141'1113、(151'1113、(161'1118与511'1118的比值，获得四个比值， 通过对比值大小判断确定有用信号滤除噪声完成自适应滤波W及确定脉搏所在的信号，进 一步在此信号上过零法计算脉搏。2. 根据权利要求1所述的自适应滤波及计算脉搏的系统，其特征在于，所述的采样模块 还包括： 其设置的采样频率为62.5监，采样的信号S的缓存大小为512B。3. 根据权利要求1所述的自适应滤波及计算脉搏的系统，其特征在于，所述的变换模 块： 对采样信号S进行6尺度小波分解，其中的基小波为化6。4. 根据权利要求1所述的自适应滤波及计算脉搏的系统，其特征在于，所述的计算模块 包括： 其中MS的运算公式为5. 根据权利要求1所述的自适应滤波及计算脉搏的系统，其特征在于，所述的自适应模 块包括： 用于计算derms与Slrms的比值，若derms与Slrms比值大，则脉搏所在的信号在d6上，d6 上使用过零法计算脉搏，得到有用的信号为d3+d4+d5+d6; 若dGrms与Slrms比值小，且dSrms与Slrms比值大，则脉搏所在的信号在d5上，d5上使用 过零法计算脉搏，得到有用的信号为d3+d4+d5; 若dGrms与Slrms比值小，且dSrms与Slrms比值小，同时(Mrms与Slrms比值大，则脉搏所 在的信号在d4上，d4上过零法计算脉搏，则有用的信号为d3+d4。6. -种自适应滤波及计算脉搏的方法，其特征在于，该方法包括： 对传感器输出信号进行采样，将采样的信号S存入缓存； 获取缓存中的采样信号S，对采样信号S进行小波分解，获得近似系数cA6及细节系数 cDl、cD2、cD3、cD4、cD5、cD6，进一步对细节系数cD3、cD4、cD5、cD6分别进行重构得到信号 d3、d4、d5、d6; 对信号(13、(14、(15、(16分别计算1^8值，获得对应的(131'1113、(141'1113、(151'1113、(161'1113，进一步计 算重构总信号S1 =d3+d4+d5+d6的RMS，获得Slrms; 分别计算d化ms、d4;rms、d虹ms、(Wrms与Slrms的比值，获得四个比值，通过对比值大小 判断确定有用信号滤除噪声完成自适应滤波W及确定脉搏所在的信号，进一步在此信号上 过零法计算脉搏。7. 根据权利要求6所述的自适应滤波及计算脉搏的方法，其特征在于，该方法还包括： 设置的采样频率为62.5监，采样的信号S的缓存大小为512B。8. 根据权利要求6所述的自适应滤波及计算脉搏的方法，其特征在于，该方法还包括： 对采样信号S进行6尺度小波分解，其中的基小波为化6。9. 根据权利要求6所述的自适应滤波及计算脉搏的方法，其特征在于，该方法还包括： MS的运算公式为10. 根据权利要求6所述的自适应滤波及计算脉搏的方法，其特征在于，该方法还包括： 计算derms与Slrms的比值，若derms与Slrms比值大，则脉搏所在的信号在d6上，d6上使 用过零法计算脉搏，得到有用的信号为d3+d4+d5+d6; 若dGrms与Slrms比值小，且dSrms与Slrms比值大，则脉搏所在的信号在d5上，d5上使用 过零法计算脉搏，得到有用的信号为d3+d4+d5; 若dGrms与Slrms比值小，且dSrms与Slrms比值小，同时(Mrms与Slrms比值大，则脉搏所 在的信号在d4上，d4上过零法计算脉搏，则有用的信号为d3+d4。
【文档编号】A61B5/021GK106073741SQ201610375711
【公开日】2016年11月9日
【申请日】2016年5月30日
【发明人】金瑞军, 符文剑
【申请人】珠海脉动时代健康科技有限公司