exam9_1.m

《小波分析理论与Matlab7实现例程》一书里面各章节的源码

clc;
clear;
% 计算与bior3.5相关的分解滤波器和重构滤波器
[Rf,Df] = biorwavf('bior3.5');
% 计算需要的4个滤波器
[Lo_D,Hi_D,Lo_R,Hi_R] = biorfilt(Df,Rf); 
subplot(221); stem(Lo_D); 
title('bior3.5分解低通滤波器'); 
subplot(222); stem(Hi_D); 
title('bior3.5分解高通滤波器'); 
subplot(223); stem(Lo_R); 
title('bior3.5重构低通滤波器'); 
subplot(224); stem(Hi_R); 
title('bior3.5重构高通滤波器');
fftld=fft(Lo_D);ffthd=fft(Hi_D);
frq=[1:length(Lo_D)]/length(Lo_D);
subplot(221);plot(frq,abs(fftld),frq,abs(ffthd));
title('分解滤波器变换模型')
fftlr=fft(Lo_R);ffthr=fft(Hi_R);
frq=[1:length(Lo_R)]/length(Lo_R);
subplot(222);plot(frq,abs(fftlr),frq,abs(ffthr));
title('重构滤波器变换模型')
subplot(212);plot(frq,abs(fftld.*fftlr+ffthd.*ffthr));
title('一种双正交性情况')