reconstructed.m

信号压缩与解压程序代码和文件

function Xk=reconstructed(X,Tx,DD,Td,step)
M=length(Tx);H=length(Td);
k0=0;k1=0;Xh=zeros(1,M);D=zeros(1,H);
for m=1:M
    if Tx(m) > 0
        k0=k0+1;
        Xh(m)=X(k0);
    end
end
for h=1:H
    if Td(h) > 0
        k1=k1+1;
        D(h)=DD(k1);
    end
end
[h0,h1,g0,g1]=filtercoefficience;
n=length(Xh);
K=length(h0);
Xn=h0(1:K)-i*h1(1:K);
for t=1:step
    N=length(Xh);
    d=D(((2^(t-1)-1)*n+1):(2^t-1)*n);
    Xt=twoinsert(Xh);Dt=twoinsert(d);
    Yn=Xt(1:2*N)+i*Dt(1:2*N);
    Yk=fft(Yn,2*N);Xk=fft(Xn,2*N);
    XYn=ifft(Yk(1:2*N).*Xk(1:2*N));
    Xh=real(XYn);
end      
Xk=Xh';