decomposition.m

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

% 小波分解与程序,Xk0是要分解的原始信号,step是表示要分解的层数
function [Xh,D]=decomposition(Xk0,step)
[h0,h1,g0,g1]=filtercoefficience;
Xh=Xk0';D=0;
for lstep=1:step
    N=length(Xh);
    Yk=fft(Xh,N);
    H0k=fft(h0,N);
    H1k=fft(h1,N);
    Xh=drawwing(ifft(Yk(1:N).*conj(H0k(1:N))));
    Xh=real(Xh);D=real(D);
    d=drawwing(ifft(Yk(1:N).*conj(H1k(1:N))));
    thrd= ddencmp('den','wv',d);%用默阈值进行消噪处理
    Td=abs(d)>thrd;
    d=Td.*d;
    D=[d,D];
end