r11.m

该程序基于小波变换

snr=1;
init=2055615866;
[s,ns]=wnoise(2,11,snr,init);
s=s(1:2000);
ns=ns(1:2000);
subplot(3,1,1);
plot(s);
title('原始信号');
%axis([0,2000,1,4]);
subplot(3,1,2);
plot(ns);
title('加噪信号');
%axis([0,2000,-1,4]);
%用小波对加噪信号进行lev层分解。
wname='sym4';lev=5;
[c,l]=wavedec(ns,lev,wname);
%对信号进行消噪
alpha=2;
[thr,nkeep]=wdcbm(c,l,alpha);
ds=wdencmp('lvd',c,l,wname,lev,thr,'s');
subplot(3,1,3);;
plot(ds);
title('消噪信号');
%axis([0,2000,1,4]);
per=norm(ds)/norm(s)
err=norm(ds-s)