program_10_16.m

这是基于MATLAB系统上开发的,应用于小波分析方面,很有用的程序模块,书名为"MATLAB在小波分析中的应用"

load noisbloc
s= noisbloc;
ls=length(s);
plot(s);
[swa,swd] = swt(s,1,'db4');
subplot(2,1,1);plot(swa);title('低频')
subplot(2,1,2);plot(swd);title('高频')
A0 = iswt(swa,swd,'db4');
err=norm(s-A0)
figure(2)
subplot(3,1,1);plot(s);title('原始信号')
subplot(3,1,2);plot(A0);title('重构信号')
subplot(3,1,3);plot(s-A0);title('误差信号')
nulcfs = zeros(size(swa));
A1 = iswt(swa,nulcfs,'db4');
D1 = iswt(nulcfs,swd,'db4');
figure(3)
subplot(2,1,1), plot(A1); 
title('重构低频A1');
subplot(2,1,2), plot(D1); 
title('重构高频D1');
[swa,swd] = swt(s,3,'db4');
kp = 0;
for i = 1:3
subplot(3,2,kp+1), plot(swa(i,:));
title(['低频系数：level ',num2str(i)])
subplot(3,2,kp+2), plot(swd(i,:));
title(['高频系数：level ',num2str(i)])
kp = kp + 2;
end
mzero = zeros(size(swd));
A = mzero;
A(3,:) = iswt(swa(3,:),mzero,'db4');  %  重构低频分量
D = iswt(mzero,swd,'db4');  %  重构所有高频分量
s_r=A(3,:)+D;  %  重构源信号
norm(s_r-s)    %  误差
figure(4)
subplot(2,2,1);
plot(s);        %  源信号
title('源信号');
subplot(2,2,2)
plot(A(3,:));   %  信号的低频部分
title('信号的低频部分');
subplot(2,2,3)
plot(D);        %  信号的高频部分
title('信号的高频部分');
subplot(2,2,4)
plot(s_r);      %  重构信号
title('重构信号');