[发明专利]一种基于AR模型谱估计的脉搏信号随机噪声去噪方法

说明书

本发明涉及一种基于AR模型谱估计的脉搏信号随机噪声去噪方法，属信息科学与医学融合领域。本发明首先将采集到的脉搏信号分为长度相同的两段，计算信号的功率信噪比，建立随机噪声信号的模型，设计随机噪声信号；分别给两段信号补零使其长度为N，对补零后信号进行傅里叶变换保留其相位谱；对补零后信号，建立AR模型，估计其功率谱；对噪声信号补零，使其长度亦为N，再估计其功率谱；用混合信号的功率谱减去噪声信号的功率谱，得到有效信号的功率谱，结合去噪前混合信号的相位谱，通过变换得到时域脉搏有效信号。本发明在明显去除了随机噪声的基础上又没有降低信号分辨率及保真度，去噪效果很好。

技术领域

本发明涉及一种基于AR模型谱估计的脉搏信号随机噪声去噪方法，属于信息科学与医学融合技术领域。

背景技术

人体内部各个生理系统之间是相互耦合的。反映人身体健康状态相对最重要、最全面的是心脏血液循环系统，因此通过采集脉搏波进而分析心脏循环系统功能，能较全面地反映人体的健康情况。但从人体采集到的脉搏波信号，由于脉搏波信号信噪比比较低，给后续参数的准确测量带来了困难，所以对于噪声干扰的去除是非常重要而必须的。

常用的脉搏信号去噪方法主要包括仅在时域内进行简单降噪处理的方法或传统的建立在傅立叶变换基础上的去噪方法。仅在时域内进行简单降噪处理，即采用低通、带通、高通等常规的滤波方法对脉搏波信号进行滤波去噪，然而常用的滤波器截止频率固定，在噪声频率超过其截止频率时，无法消除噪声，截止频率过高会滤除有用信息。当噪声与信号的频谱相近或重叠时，传统的单纯时域滤波或频域滤波往往无法达到很好的效果。传统的建立在傅立叶变换基础上的去噪方法，在提高信噪比和提高空间分辨率两项指标上存在矛盾。低通滤波能通过平滑抑制噪声，但同时会使信号的边沿变模糊。高通滤波可以使信号边沿更加的陡峭，但背景噪声也同时被加强了。小波变换作为一种新型的时域分析法，由于其具有良好的时频局部性，并且有快速算法(Mallat算法)加以实现，因而在去噪领域受到了越来越多的这类方法处理比较简单，运算量较小，基于小波变换的信号去噪方法，一般有阈值法、空域相关滤波法和模极大值法。但采用阈值法对脉搏波信号进行去噪，当信号具有突变的不连续点或信噪比较低时，阈值法去噪会出现伪吉布斯现象，即在不连续点附近的信号会在原有信号电平上上下跳变。空域相关滤波它能够保留信号突变点位置携带的重要信息，从而对噪声有很好的滤波效果，但是传统的空域相关滤波算法很容易引入随机噪声。基于模极大值的小波去噪算法，小波分解尺度的选择非常重要，尺度过小，小波系数受噪声影响非常大，产生许多伪极值点，尺度过大又会使信号丢失某些重要的局部奇异性，因此需要选择合适的尺度。

发明内容

本发明提供了一种基于AR模型谱估计的脉搏信号随机噪声去噪方法，通过估计脉搏信号功率信噪比，进而建立随机噪声信号模型，提出了利用AR模型谱估计，在频域内去除随机噪声的方法，用于解决有用信号和噪声信号在频域宽度内重合时，用传统时域滤波和频域滤波方法无法较好分离有用信号和噪声信号的问题；也用于解决基于小波变换的去噪方法会产生伪吉布斯现象或引入随机噪声的问题；且本发明中，确定AR参数模型的阶数为未补零前的信号长度，参数的选择简单且固定，但去噪效果依然很好，解决了参数选择困难及因参数选择不当而影响到去噪效果的问题。

本发明的技术方案是：一种基于AR模型谱估计的脉搏信号随机噪声去噪方法，所述基于AR模型谱估计的脉搏信号随机噪声去噪方法的具体步骤如下：

Step1、将采集到的脉搏信号分为长度相同的两段来计算信号的功率信噪比，通过计算得到的信号功率信噪比估计噪声方差，进而建立随机噪声信号的模型，设计随机噪声信号；

Step2、分别给两段信号补零使两段信号的大小均为最接近原信号长度的2的整数次方，对补零后的混合信号进行傅里叶变换保留其相位谱；

Step3、对补零后的两段混合信号建立AR模型，确定AR模型及模型阶数，根据模型求出模型的参数，代入功率谱密度估计公式估计出补零后的混合信号的功率谱；其中模型阶数选择未补零前的每段信号的长度；