sinsym4.m

该程序基于小波变换

%给定并显示一个正弦信号
t=0:1000;
s=sin(0.05*t);
subplot(4,2,1);
plot(s);
axis([0 1000 -1 1]);
title('原始信号');
%===============================
%给信号加噪并显示
ns=s+randn(1,length(t));
subplot(4,2,2);
plot(ns);
axis([0 1000 -1 1]);
title('加噪信号');
%=================================
%对信号进行消噪并显示
%用小波函数sym4对信号进行5尺度小波分解
[c,l]=wavedec(ns,6,'sym4');
%对尺度1上的低频部分进行重构
a1=wrcoef('a',c,l,'sym4',1);
subplot(4,2,3);
plot(a1);
axis([0 1000 -1 1]);
title('第一次消噪信号');
%对尺度2上的低频部分进行重构
a2=wrcoef('a',c,l,'sym4',2);
subplot(4,2,4);
plot(a2);
axis([0 1000 -1 1]);
title('第二次消噪信号');
%对尺度3上的低频部分进行重构
a3=wrcoef('a',c,l,'sym4',3);
subplot(4,2,5);
plot(a3);
axis([0 1000 -1 1]);
title('第三次消噪信号');
%对尺度4上的低频部分进行重构
a4=wrcoef('a',c,l,'sym4',4);
subplot(4,2,6);
plot(a4);
axis([0 1000 -1 1]);
title('第四次消噪信号');
%对尺度5上的低频部分进行重构
a5=wrcoef('a',c,l,'sym4',5);
subplot(4,2,7);
plot(a5);
axis([0 1000 -1 1]);
title('第五次消噪信号');
%对尺度5上的低频部分进行重构a6=wrcoef('a',c,l,'sym4',6);
subplot(4,2,8);
a6=wrcoef('a',c,l,'sym4',6);
plot(a6);
axis([0 1000 -1 1]);
title('第六次消噪信号');