%生成正弦信號
N=1000;
t=1:N;
x=sin(0.03*t);
%加噪聲
load noissin;
ns=noissin;
%顯示波形
subplot(3,1,1);
plot(t,x);
xlabel('樣本序號 n');
ylabel('(原始信號)幅值A');
subplot(3,1,2);
plot(ns);
xlabel('樣本序號 n');
ylabel('(含噪信號)幅值A');
%小波消噪
xd=wden(ns,'minimaxi','s','one',5,'db3');
subplot(3,1,3);
plot(xd);
xlabel('樣本序號 n');
ylabel('(消噪信號)幅值A');