例程9-5.m

来自「《小波分析理论与Matlab7实现》书中例程」· M 代码 · 共 31 行

% 装载尺度滤波器
load db8; w = db8; 
subplot(421); stem(w); 
title('原始尺度滤波器');
% 计算4个滤波器
[Lo_D,Hi_D,Lo_R,Hi_R] = orthfilt(w); 
subplot(423); stem(Lo_D); 
title('分解低通滤波器'); 
subplot(424); stem(Hi_D); 
title('分解高通滤波器'); 
subplot(425); stem(Lo_R); 
title('重构低通滤波器'); 
subplot(426); stem(Hi_R); 
title('重构高通滤波器');
% 检验正交性
df = [Lo_D;Hi_D];
rf = [Lo_R;Hi_R];
id = df*df'
id = rf*rf'
% 检验二元转换后的正交性
df = [Lo_D 0 0;Hi_D 0 0]; 
dft = [0 0 Lo_D; 0 0 Hi_D]; 
zer = df*dft'
% 高频和低频解释
fftld = fft(Lo_D); ffthd = fft(Hi_D); 
freq = [1:length(Lo_D)]/length(Lo_D); 
subplot(427); plot(freq,abs(fftld)); 
title('转换模量:低通');
subplot(428); plot(freq,abs(ffthd)); 
title('转换模量:高通')