【转】由滤波器系数绘制尺度函数和小波函数图像的Matlab程序

原帖地址：http://blog.csdn.net/LiuHB1127/archive/2007/12/02/1912474.aspx

function ScaleWaveFig(h)
% -- 函数描述 : 由滤波器系数绘制尺度函数和小波函数图像
% M : 标准化常数
% h : (尺度)滤波器系数
% g : 小波滤波器系数
% a : 尺度函数初始化
% w : 小波函数初始化
% -- 时间 : 2007-12-02
% -- 作者 : 刘恒冰(LIUHB) 版权所有(C)

M = 2;
g = fliplr(h);
for i = 1 : length(h)
g(i) = (-1) ^ (i + 1) * g(i);
end
a = h;
w = g;

% 绘制尺度函数图像
b = [ ];
for i = 1 : 7
L = M * length(a);
b(1 : M : L - M + 1) = a;
for j = 2 : M
b(j : M : L - M + j) = zeros(1, L / M);
end
a = b;
a = conv(h, a);
% a = sqrt(M) * a; || a = sqrt(M) * a; ?
n = length(a);
a = a(1, 1 : n - 1);
end
n = length(a);
x = linspace(0, 3, n);
subplot(221);
plot(x, a); grid on;

% 绘制小波函数图像
b = [ ];
for i = 1 : 7
L = M * length(w);
b(1 : M : L - M + 1) = w;
for j = 2 : M
b(j : M : L - M + j) = zeros(1, L / M);
end
w = b;
w = conv(h, w);
% w = sqrt(M) * w; || w = sqrt(M) * w; ?
n = length(w);
w = w(1, 1 : n - 1);
end
n = length(w);
x = linspace(0, 3, n);
subplot(222);
plot(x, w); grid on;