由濾波器系數(shù)繪制尺度函數(shù)和小波函數(shù)圖像的Matlab程序 收藏
function ScaleWaveFig(h)
% -- 函數(shù)描述 : 由濾波器系數(shù)繪制尺度函數(shù)和小波函數(shù)圖像
%    M : 標(biāo)準(zhǔn)化常數(shù)
%    h : (尺度)濾波器系數(shù)
%    g : 小波濾波器系數(shù)
%    a : 尺度函數(shù)初始化
%    w : 小波函數(shù)初始化
% -- 時(shí)間 : 2007-12-02
% -- 作者 : 劉恒冰(LIUHB)  版權(quán)所有(C)

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

% 繪制尺度函數(shù)圖像
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;

% 繪制小波函數(shù)圖像
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;

本文來(lái)自CSDN博客，轉(zhuǎn)載請(qǐng)標(biāo)明出處：http://blog.csdn.net/LiuHB1127/archive/2007/12/02/1912474.aspx