dec.m

文件包含三个子文件：dec.m,test_dec.m,recon.m dec.m文件使用MATLAB语言自己动手编写小波分解算法； test_dec.m使用MATLAB语言测试dec.m文件中自己

function w=dec(f,p,NJ,Jstop)
%f--待分析信号，长度是2^NJ
%p---尺度系数
%Jstop---停止尺度级别
%w-----
L=length(p)
pf=fliplr(p);
q=p;
q(2:2:L)=-q(2:2:L);

a=f;
t=[];
for j=NJ:-1:Jstop+1
    n=length(a)
    mod((-L+1:-1),n)+1
    a(mod((-L+1:-1),n)+1)
    a=[a(mod((-L+1:-1),n)+1) a]; %make the data periodic
    b=conv(a,q);
    b=b(L+1:2:L+n-1)/2;  %down sample
    
    a=conv(a,pf);
    a=a(L:L+n-1)/2;
    ab=a(1:L);
    a=[a(L+1:n) ab];%periodize
    a=a(2:2:n)   %down sample
    
    t=[b,t];
end

w=[a,t];  
JJ=2^Jstop;
ww=[w(JJ) w(1:JJ)];
tt=linspace(0,1,JJ+1);
plot(tt,ww);