自己编写的小波去噪函数.docVIP

利用小波去噪的M文件 局部放电试验所采集的信号中往往混有白噪声、周期干扰信号去除。此处采用常用db系列小波中的db6小波进行9尺度的多分辨分解后，根据白噪声能量特性，估算各尺度的阈值大小，采用硬值进行处理，后进行重构。Matlab程序如下： function sd xiao_denoise mix_signal %此函数用于去除白躁信号＆周期性干扰信号 %输入参数mix_signal为采集到的信号波形 p 0.6745; w_dept 9; w_name db6; coef cell 1,w_dept ; thr zeros 1,w_dept+1 ; [c,l] wavedec mix_signal,w_dept,w_name ; %对混合信号S进行db6的9尺度一维分解 coef 1 appcoef c,l,w_name,w_dept ;%计算尺度为9的一维分解低频系数 cs [cs,coef_soft j ]; thr 1 median abs coef 1 /p*sqrt 2*log length coef 1 ;%计算1尺度上的阈值 coef_soft 1 wthresh coef 1 ,h,thr 1 ;%对小波系数进行阈值为thr 1 的硬阈值处理 cs [coef_soft 1 ]; for j 2:w_dept+1 coef j detcoef c,l,w_dept-j+2 ;%计算尺度为9到2的各尺度高频小波系数 coef1 j detcoef c,l,w_dept-j+2 ; thr j median abs coef j /p*sqrt 2*log length coef j ;%计算9到2各尺度上的阈值 coef_soft j wthresh coef j ,h,thr j ;%对小波系数进行阈值为thr j 的硬阈值处理 cs [cs,coef_soft j ]; end sd waverec cs,l,w_name ; %根据小波系数[cs,l]对信号进行重构 实例验证 为了验证去噪的有效性，先仿真产生一个局放脉冲然后叠加0.1倍白噪声和周期干扰，利用前面的程序去造，结果如图1，从图上可以看到去噪后信号与原始信号幅值、相位都基本没有变化程序如下： fc 40e4; %振荡频率 t4 0.8e-3; %脉冲起始时间 tn 1e-3; %总时间 x 0:step:tn; x4 t4:step:tn; %s4 exp t4-x4 *13/t -exp t4-x4 *22/t .*sin 2*pi*fc*x4 ; s4 exp t4-x4 /tr -exp t4-x4 /td .*sin 2*pi*fc*x4 ; s4 [zeros 1,t4/step ,s4]; p tn/step; n 0.1*randn 1,p ; %产生白噪信号 n [n,0]; s5 0.1*sin 2*pix ; %产生周期性干扰信号 s6 s4+n+s5; sd liu_denoise s6 ; subplot 311 ;plot x,s4 ;title 单个局放脉冲仿真波形 ; subplot 312 ;plot x,mix_signal ;title 染噪后波形 ; subplot 313 ;plot x,sd ;title 小波去噪后波形 ; 图1 小波去噪对比图