[发明专利]一种小波去噪方法

权利要求书

1.一种小波去噪方法，其特征在于，包括：

步骤S1：对染噪信号进行小波分解得到各层小波系数和尺度系数；

步骤S2：分别计算各层小波系数和尺度系数中各系数的峭度；

步骤S3：将峭度小于预设阈值的系数识别为噪声对应系数；

步骤S4：利用3σ准则分别计算各层小波系数和尺度系数中噪声对应系数的阈值作为该层小波阈值a_thr＝3σ，其中σ为各层小波系数和尺度系数中噪声对应系数的标准差；

步骤S5：结合小波阈值和卡方能量窗法从各层小波系数和尺度系数中提取脉冲对应系数；

步骤S6：将各层小波系数和尺度系数中除了脉冲对应系数以外的系数置零得到处理后各层小波系数和尺度系数；

步骤S7：利用处理后各层小波系数和尺度系数重构得到去噪信号。

2.根据权利要求1所述的一种小波去噪方法，其特征在于，步骤S2中所述的各系数的峭度的计算公式为：

其中，x_i(i＝1,2,…,N)为以待计算峭度的各系数为中心的系数序列，N为系数序列的长度，N取值应小于脉冲信号半波长度，μ是系数序列的平均值，σ为系数序列的标准差。

3.根据权利要求2所述的一种小波去噪方法，其特征在于，所述小波分解为二进小波分解，系数序列的长度N满足：最高层小波系数所用的N和最高层尺度系数所用的N相同，次高层小波系数所用的N为最高层小波系数所用的N的两倍，以此类推确定各层小波系数所用的N。

4.根据权利要求1所述的一种小波去噪方法，其特征在于，所述预设阈值的取值范围为[2,5]。

5.根据权利要求1所述的一种小波去噪方法，其特征在于，所述卡方能量窗法包括：

步骤S51：根据脉冲宽度下限设定能量窗宽度M；

步骤S52：逐点比较系数的绝对值与小波阈值a_thr，绝对值大于a_thr的系数x(n)为脉冲对应系数；

步骤S53：从识别的x(n)开始，以x(n)为终点的能量窗逐点向前滑动，当能量窗标称累积能量E小于能量阈值e_thr时，能量窗终点对应系数为脉冲起点，其中

步骤S54：从识别的x(n)开始，以x(n)为起点的能量窗逐点向后滑动，当能量窗标称累积能量E小于能量阈值e_thr时，能量窗起点对应系数为脉冲终点，脉冲终点和步骤S53所述脉冲起点之间系数均为脉冲对应系数，其中

步骤S55：从与步骤S54识别的脉冲终点间隔M点的下一个系数开始，重复步骤S52至步骤S55直到系数结尾；

其中，x(i)为能量窗内系数序列，σ为各层小波系数或最高层尺度系数中噪声对应系数的标准差；

所述能量阈值e_thr由卡方能量窗宽度M查询卡方分布表得到。

6.一种非暂态计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现如权利要求1至5任一项所述的小波去噪方法。

7.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现如权利要求1至5任一项所述的小波去噪方法。