% 裝載尺度濾波器
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('轉換模量:高通')