小波源程序.m

这种小波构造方式适合于各种DB系列的小波

mm=16;
%% CONSTRUCT transform matrix ww
h=wfilters('db1','l');
rank=log2(mm);
ww=1;
for jj=1:4%如果是二阶分解的话就是jj=3:4
    nn=2^jj;
    clear w1 p1 p2 p1_0 p2_0;
    for ii=1:length(h)
        p1_0(ii)=h(length(h)-ii+1);
        p2_0(ii)=(-1)^ii*h(ii);
    end
    p1_0=sparse([p1_0,zeros(1,nn-length(h))]);
    p2_0=sparse([p2_0,zeros(1,nn-length(h))]);

    for ii=1:nn/2
        p1(ii,:)=circshift(p1_0',2*(ii-1))';
        p2(ii,:)=circshift(p2_0',2*(ii-1))';
    end
    w1=[p2;p1];
    clear p1;
    clear p2;
    L=2^rank-length(w1);
    w=sparse([eye(L,L),zeros(L,length(w1));zeros(length(w1),L),w1;]);
    ww=ww*w;
    clear w;
end
%这种小波构造方式适合于各种DB系列的小波，但面对db6等比较长的小波，则均衡器的权长要足够长。