基于小波去噪的激光多普勒振动信号处理

谈 渊，甘学辉*，张东剑，刘香玉，廖 壑

(1.上海市高性能纤维复合材料省部共建协同创新中心，上海 201620；2.东华大学 机械工程学院，上海 201620)

引 言

激光多普勒测振技术是基于多普勒效应通过解析被测物体的振动所产生的多普勒频移量，来得到物体的振动状态[1-2]。由于其响应速度快、检测精度高，在运动状态检测、模态分析、结构损伤识别等领域应用越来越广泛[3-6]。然而由于外界环境的干扰，采集的激光干涉信号不仅携带被测物体的振动信号，而且携带更多的是外界环境噪声；甚至在大多数情况下,振动信号会被噪声完全掩盖。因此,如何实现强噪声环境下提取目标振动信号成为首要任务和难点[7]。

小波变换能同时对信号时域和频域进行多尺度细化分析，较好地区分信号中的突变部分和噪声，进而对信号进行去噪[8]。其中小波阈值去噪法被广泛应用于激光多普勒振动信号处理[9-11],但是由于其去噪会产生固定偏差,许多学者也对小波阈值去噪算法做了改进。XIANG等人[12]将样本熵表征噪声强度引入阈值函数，使阈值函数更适应小波系数的分布，但是用表征信号强度的熵表征噪声强度会造成较大偏差。LIU等人[13]构造了随尺度变化的局部阈值，同时引入正弦函数改进阈值函数，解决了软硬阈值函数存在的不足，但是仅将信噪比作为指标会导致无法找到最优分解层数。HAO等人[14]将均方根误差和信噪比作为评价指标选择最优分解层数，解决了单一评价指标存在的不足，但是单纯用现有指标线性叠加仍会产生偏差。上述学者在阈值函数和分解层数评价指标方面对算法做了改进，去噪效果有一定提升，但仍存在问题：难以寻找最优分解层数、忽略有用信号小波系数随尺度的分布。

本文中将尺度引入阈值函数模型、选用局部阈值、将均方根误差变化量、平滑度变化量处理后线性叠加作为新的分解层数评价指标来改进小波阈值去噪算法，将改进后算法用于激光多普勒测振信号处理，并且通过仿真和音叉去噪实验验证了改进算法的有效性。

1 小波阈值去噪原理及模型改进

1.1 小波阈值去噪原理

小波阈值去噪原理是通过估计噪声强度生成阈值，利用阈值函数处理小波系数去除部分噪声相关的分量，从而提高有用信号的比重。主要分为以下3个步骤：(1)通过小波变换将含噪信号分解成若干层的小波系数;(2)利用阈值函数处理各层小波系数去除与噪声有关的分量;(3)通过小波逆变换重构信号[15]。

1.2 改进的去噪模型

阈值是影响去噪效果的一个重要因素，低尺度的小波系数是大于高尺度的，包含的噪声系数也越大，对不同尺度利用相同阈值处理，会影响去噪效果。因此本文中使用局部阈值。局部阈值λj的公式为[16]：

(1)

式中，σj是第j层小波系数的噪声标准差，wj是第j层小波系数数列，λ是根据阈值规则计算得出的阈值。(1)式可以满足小尺度下取较高阈值，避免全局阈值产生恒定偏差。median为中值函数，表示求数列中值。

小波阈值去噪关键在于阈值函数的选取[17]，针对软阈值和硬阈值的缺点，提出了基于局部阈值的改进阈值函数。改进后阈值函数g(x) 形式为:

(2)

式中，x是小波系数数列中的数。

(3)

改进后的阈值函数在x=λj处连续，克服了硬阈值的缺点。

6.牛流行热。潜伏期为3～7 d， 体温升高到40℃以上，稽留2～3 d后恢复正常。结膜充血，水肿，呼吸促迫困难，发出呻吟声，呼吸次数每分钟可达80次以上。妊娠母牛可发生流产、死胎，乳量下降或停止。

当x>0时：

(4)

f=x-g(x)>0

(5)

(6)

当x<0时：

(7)

f=g(x)-x>0

(8)

(9)

由(4)式～(9)式可知,g(x)=x是阈值函数的一条渐近线，随着x变大处理后的小波系数向原小波系数不断逼近，同时随着尺度增加，逼近程度不断增加，符合小波系数的组成随尺度的变化情况，解决了软阈值函数产生较大偏差的问题。

采集的激光多普勒振动信号中噪声并不是稳定的，某一时刻外界噪声可能急剧变大，每次采样都需要寻找最优分解层数，一般来说分解层数为1～8[18]。参考现有融合评价指标，作者提出了自适应分解层数[19-20]，选取均方误差和平滑度融合得到评价指标T，来选择最优分解层数。具体计算方法如下。

Δe(m)=|e(m+1)-e(m)|

(10)

Δr(m)=|r(m+1)-r(m)|

(11)

式中，e(m)表示均方根误差；r(m)表示平滑度，Δe(m)表示m尺度下均方根误差变化量；Δr(m)表示m尺度下平滑度变化量。

将上述两个指标归一化处理：

Δe′(m)=[Δe(m)-Δemin]/(Δemax-Δemin)

(12)

Δr′(m)=[Δr(m)-Δrmin]/(Δrmax-Δrmin)

(13)

计算两个指标的权重系数：

(14)

(15)

最后通过线性组合叠加得到评价指标T：

T=We×Δe′+Wr×Δr′

(16)

2 去噪仿真实验

为验证改进小波阈值去噪算法的效果，通过MATLAB创建信噪比为5dB，有用信号为S=3sin(t+π/4)+4cos(t+π/3)的含噪信号,其中S表示生成的模拟信号，t是时间。利用自己改进的算法和现有软阈值、硬阈值、参考文献[14]中改进算法对生成的模拟信号去噪。实验中采用db4基小波，无偏阈值规则，先利用本文中提出的改进算法去噪可以得到一个信噪比，在保持阈值规则、小波基不变情况下采用软阈值、硬阈值和参考文献[14]中去噪算法，比较去噪后信噪比的大小。图1为含噪信号，图2为参考文献[14]中算法去噪后信号，图3为软阈值法去噪后信号，图4为硬阈值法去噪后信号，图5为本文中提出的改进算法去噪后信号。

从上述图中可以看出,由于模拟信号复杂程度较低，所以采用上述4种方法去噪效果都较好，不存在突出尖峰；但是软硬阈值法去噪后峰值相差较大，并且曲线不光滑；参考文献[14]中算法和软硬阈值法相比,去噪效果有了一定改进，但是仍然存在局部的不光滑曲线;本文中提出的改进去噪算法峰值平齐，还接近实际峰值，同时去噪后曲线平滑程度高。为进一步突出表现去噪效果，分别选择db4,sym6,coif3小波作为小波基比较去噪后信噪比(signal-to-noise ratio,SNR)。表1中为各种算法去噪后信噪比，表2中为改进算法处理仿真信号信噪比比原有算法处理效果提升百分比。

Fig.1 Noisy signal

Fig.2 Signal processed by validate algorithm

Fig.3 Signal processed by soft threshold algorithm

Fig.4 Signal processed by hard threshold algorithm

Fig.5 Signal processed by the improved algorithm

Table 1 SNR of each algorithm after denoising

Table 2 Improvement percentage of SNR

从表1和表2中的数据可以看出：采用改进的算法去噪效果都比原有方法好，用不同小波基时改进效果也会发生改变，改进的算法比原有算法信噪比至少提升19.4%，比参考文献[14]中算法信噪比仍提升约7.7%，表明改进的算法具有更好的去噪效果。

3 信号去噪实验

3.1 实验测试及信号处理

实际采集的激光多普勒振动信号较复杂，仿真实验并不能体现对实际信号的去噪效果，因此需要对实际的采集信号处理验证改进小波去噪算法的效果。本文中基于多普勒效应搭建如图6所示的激光多普勒测振实验平台。光源选择波长为632.8nm、功率为5mW的He-Ne激光器，光电探测器选择卓立汉光公司的DSi300硅光电探测器，声光频移器选择中国电子科技SGT40-633-2PA一体化声光调制器，载波频率为40MHz。以频率为512Hz的音叉为振动物体，调整光路得到干涉牛顿环后敲击音叉，通过光电探测器和采集卡传输振动信号,基于LabVIEW和MATLAB联合编写振动信号处理程序。

Fig.6 Laser Doppler experimental device

通过搭建的实验平台和软件系统，设置采样频率为5000Hz，采样数为2000，分别利用本文中提出的改进去噪算法、参考文献[14]中算法以及现有的软硬阈值函数去噪算法对信号处理，观察时域和频域信号比较去噪效果。图7是测量得到的原始信号(包括时域和频域信号)，图8是软阈值函数算法处理后频域信号，图9是硬阈值函数算法处理后频域信号，图10是参考文献[14]中算法处理后频域信号，图11是改进去噪算法处理后频域信号。

3.2 实验结果分析

如图7所示,采集到的音叉振动信号包含许多噪声干扰难以看出音叉的振动频率，同时具有很强的低频外部噪声，所以在利用去噪算法处理前采用一个低通滤波器滤除低频噪声，分别采用不同的去噪处理算法处理采集的振动信号。考虑到在音叉上粘贴了一个增强反光效率的小镜片，所以采集到的频率会与音叉固有频率(512Hz)有一定偏差。由于采样点数过多，曲线较密，时域信号较难看出去噪效果，所以在频域图里观察频率组成来检验去噪效果，由图8和图9可以看出,现有的软硬阈值函数去噪算法由于固定的分解层数限制，保留的噪声仍然过多，无法实现振动频率提取。由图10可以看出,参考文献[14]中去噪算法相比现有算法有了改进，可以提取到振动频率，但是低频和高频噪声仍然有较高幅值。观察图11可以发现,本文中提出的改进去噪算法低频和高频段振动幅值较小，因此该改进算法去噪效果更优，对高低频噪声都有较好抑制作用。

Fig.7 Raw signal

Fig.8 Frequency-domain signal processed by soft threshold algorithm

Fig.9 Frequency-domain signal processed by hard threshold algorithm

Fig.10 Frequency-domain signal processed by validate algorithm

Fig.11 Frequency-domain signal processed by improved algorithm

4 结 论

利用改进去噪算法处理信噪比为5dB的仿真信号，信噪比变为25dB，改进的算法比原有算法信噪比大约提升19.4%，提升了信号的去噪效果。随后利用改进算法处理频率为512Hz的音叉振动信号，可以得到音叉振动频率为515Hz ，由于在音叉上粘贴一个增强反光效率的小镜片会造成固有频率发生细小改变，所以,该算法在保留有用信号基础上对低频和高频噪声都有较好滤除效果。结合仿真和实验结果分析，去噪后信号信噪比提高，噪声得到去除，有用信号得以保留，验证了利用改进算法处理激光多普勒振动信号，获取振动状态是可行的。