f10_3.m

matlab时频分析配套源码第11章的程序

%生成正弦信号
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');