% 裝載正交小波相關的尺度濾波器
load db10; 
subplot(321); stem(db10); title('db10 低通濾波器');
% 計算二次鏡像濾波器
qmfdb10 = qmf(db10); 
subplot(322); stem(qmfdb10); title('QMF db10 濾波器');
% 檢查頻率條件（對于正交很有必要），對于每個頻率點：
% abs(fft(filter))^2 + abs(fft(qmf(filter))^2 = 1
m = fft(db10); 
mt = fft(qmfdb10); 
freq = [1:length(db10)]/length(db10); 
subplot(323); plot(freq,abs(m)); 
title('db10轉換模量')
subplot(324); plot(freq,abs(mt)); 
title('QMF db10的轉換模量')
subplot(325); plot(freq,abs(m).^2 + abs(mt).^2); 
title('檢查db10 and QMF db10的QMF條件') 
xlabel(' abs(fft(db10))^2 + abs(fft(qmf(db10))^2 = 1')