a53.m

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

%给定一个正弦信号并图示之
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('消噪信号');

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

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

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