f10_3.m

来自「小波变换在信号处理中的应用,包括分解」· M 代码 · 共 26 行

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