reconstructed.m

小波阈值去噪 VC编译通过,效果很好,有很大参考价值!

function Xk=reconstructed(Xh,D,step)%信号重组函数
[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=insert(Xh);Dt=insert(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';