program_09_06.m

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

% 装载正交小波相关的尺度滤波器
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')