xiaozao.m

小波信号消噪

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');