小波去噪.rar_小波 语音去噪_小波语音降噪_小波降噪方法_语音 小波降噪_语音降噪对比

clear all close all clc x = xlsread('data.xls'); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%小波去噪 %%%%%%%%%%%db4小波去噪(分解五层) [c,l]=wavedec(x,5,'db4'); % sigma=wnoisest(c,l,1); % %%%%penalty策略确定阈值 % alpha=2;%选择参数为2 % thr1=wbmpen(c,l,sigma,alpha) % xd1=wdencmp('gbl',c,l,'db4',5,thr1,'s',1); % figure;plot(f); title('原始信号波形（penalty策略）');xlabel('时间');ylabel('幅度'); % figure;plot(x); title('含噪声信号波形（penalty策略）');xlabel('时间');ylabel('幅度'); % figure;plot(f); title('降噪后信号波形（penalty策略）');xlabel('时间');ylabel('幅度'); %%%%Birge-Massart策略确定阈值 alpha2=3;%%%由于正常情况参数取3 [thr2,nkeep]=wdcbm(c,l,alpha2) xd2=wdencmp('lvd',c,l,'db4',5,thr2,'h'); %figure;plot(f); title('原始信号波形（Birge-Massart策略）');xlabel('时间');ylabel('幅度'); figure;plot(x); title('含噪声信号波形（Birge-Massart策略）');xlabel('时间');ylabel('幅度');grid on; figure;plot(xd2); title('降噪后信号波形（Birge-Massart策略）');xlabel('时间');ylabel('幅度');grid on; %%%%%运用全局阈值 [thr3,sorh,keepapp]=ddencmp('den','wv',x); xd3=wdencmp('gbl',x,'db4',5,thr3,sorh,keepapp); %figure;plot(f); title('原始信号波形（全局阈值）');xlabel('时间');ylabel('幅度'); figure;plot(x); title('含噪声信号波形（全局阈值）');xlabel('时间');ylabel('幅度');grid on; figure;plot(xd3); title('降噪后信号波形（全局阈值）');xlabel('时间');ylabel('幅度');grid on; %%%%%抑制细节系数 ca=wrcoef('a',c,l,'db4',5); index=l(2)+1:l(7); cl=c; cl(index)=c(index)/3; xd4=waverec(cl,l,'db4'); %figure;plot(f); title('原始信号波形（抑制细节）');xlabel('时间');ylabel('幅度'); figure;plot(x); title('含噪声信号波形（抑制细节）');xlabel('时间');ylabel('幅度');grid on; figure;plot(xd4); title('降噪后信号波形（抑制细节）');xlabel('时间');ylabel('幅度');grid on; % per1=norm(xd1)/norm(x) % err1=norm(xd1-x) per2=norm(xd2)/norm(x) err2=norm(xd2-x) per3=norm(xd3)/norm(x) err3=norm(xd3-x) per4=norm(xd4)/norm(x) err4=norm(xd4-x) per5=norm(ca)/norm(x) err5=norm(ca-x)