1.一种基于形态学滤波与提升小波变换的心电信号去噪方法,其特征在于按照以下步骤进行:步骤1:将心电信号f进行第一级提升小波分解,得到第一层低频系数CA1和第一层高频系数CD1;
步骤2:对步骤1获得的第一层低频系数CA1进行第二级提升小波分解,得到第二层低频系数CA2和第二层高频系数CD2;
步骤3:对步骤2获得的第二层低频系数CA2进行第三级提升小波分解得到第三层低频系数CA3和第三层高频系数CD3;
步骤4:采用第一提升阈值去噪法对高频系数CD3进行去噪处理,得到去噪后的高频系数CD3’,并将高频系数CD3’与步骤3获得的低频系数CA3进行提升小波重构得到系数CA2’;
步骤5:采用第二提升阈值去噪法对高频系数CD2进行去噪处理,得到去噪后的高频系数CD2’,并将高频系数CD2’与步骤4获得的系数CA2’进行提升小波重构得到系数CA10;
步骤6:采用形态学滤波法对步骤5获得的系数CA10进行处理,去除系数CA10中的高频分量f1得到系数CA1’;
步骤7:采用第三提升阈值去噪法对步骤1获得的高频系数CD1进行去噪处理,得到去噪后的高频系数CD1’,并将高频系数CD1’与步骤6获得的系数CA1’进行第三次提升小波重构,得到去噪后的心电信号f’。
2.根据权利要求1所述的基于形态学滤波与提升小波变换的心电信号去噪方法,其特征在于:所述第一提升阈值去噪法、第二提升阈值去噪法以及第三提升阈值去噪法所采用的阈值去噪函数均为:其中,CD(i)为对应高频系数第i个采样点值,CD'(i)为CD(i)去噪后的值,sign()为符号函数,λ为常数,TL与TH为两个阈值,i=1~N,N为信号的长度。
3.根据权利要求2所述的基于形态学滤波与提升小波变换的心电信号去噪方法,其特征在于:所述常数λ取值为3.5,所述阈值TL与TH的计算公式为:其中, median(CD)为对应高频系数的中值;
当δ≤0.121时,TL=0;当δ>0.121时,
4.根据权利要求1所述的基于形态学滤波与提升小波变换的心电信号去噪方法,其特征在于:步骤6中所述形态学滤波法按照以下步骤进行:步骤6-1:将步骤5获得的系数CA10同时进行一路开—闭运算和一路闭—开运算,并将两路运算结果进行算术平均得到高频分量f1;
步骤6-2:将所述系数CA10与步骤6-1获得的所述高频分量f1进行求差运算,得到系数CA1’。
5.根据权利要求4所述的基于形态学滤波与提升小波变换的心电信号去噪方法,其特征在于:所述的形态学滤波法采用直线形结构元素。