a53.asv

来自「MATLAB中常用的程序」· ASV 代码 · 共 29 行

%给定一个正弦信号并图示之
t=0:1000;
s=sin(0.03*t);
subplot(4,1,1);
plot(s);
axis([0 1000 -1 1]);
title('原始信号');
%=============================
%给该信号加噪声
load noissin;
ns=noissin;
subplot(4,1,2);
plot(ns);
title('染噪信号');
%=============================
%进行消噪处理
xd=wden(ns,'minimaxi','s','one',5,'db3');
subplot(4,1,3);
plot(xd);
title('消噪信号');

xd1=wden(ns,'minimaxi','s','one',5,'db6');
subplot(4,1,4);
plot(xd1);
title('消噪信号1');

p1=1/length(s)*norm(s).^2; 
p2=1/length(s)*norm(ns-xd).^2; 
snr=10*log(p1/p2)