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