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

原帖地址：http://blog.youkuaiyun.com/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;