%給定一個(gè)正弦信號(hào)并圖示之
t=0:1000;
s=sin(0.03*t);
subplot(3,1,1);
plot(s);
axis([0 1000 -1 1]);
title('原始信號(hào)');
%=============================
%給該信號(hào)加噪聲
load noissin;
ns=noissin;
subplot(3,1,2);
plot(ns);
title('染噪信號(hào)');
%=============================
%進(jìn)行消噪處理
xd=wden(ns,'minimaxi','s','one',5,'db3');
subplot(3,1,3);
plot(xd);
title('消噪信號(hào)');