recon.m

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

function y=recon(w,p,NJ,Jstart)



N=length(w)
Nj=2^Jstart;

L=length(p)
q=fliplr(p);
a=w(1:Nj);
for j=Jstart:(NJ-1)
    b=w(Nj+1:2*Nj);
    m=mod((0:L/2-1),Nj)+1;
    Nj=2*Nj;
    ua(2:2:Nj+L)=[a a(1,m)]; %periodize the data and upsample
    ub(2:2:Nj+L)=[b b(1,m)];
    
    ca=conv(ua,p);
    ca=[ca(Nj:Nj+L-1) ca(L:Nj-1)];
    
    cb=conv(ub,q);
    cb=cb(L:Nj+L-1);
    cb(1:2:Nj)=-cb(1:2:Nj);    %sign change on the odd entries
    
    a=ca+cb;
end

y=a;
yy=[y(N) y];
t=linspace(0,1,N+1);
plot(t,yy)