基于小波变换的滤波算法matlab程序

function ynewq = func_wavelet_filter(Signal);

num    = 220;
coef   = 3;
ynewq  = zeros(1,length(Signal));
x      = zeros(1,1024);
xx     = zeros(1,num);
b      = zeros(1,21);
 
for i=1:1:length(Signal)
    if i<=1024
       ynewq(i) = Signal(i);
       x(i)     = Signal(i);
    else
        if i==1025
           for m=1:1:1023-num
               x(m) = x(m+1+num);
           end  
           x(1024-num)=Signal(i);  
           for p = 1025-num:1:1024
               x(p)=x(2*(1024-num)-p+1);
           end
           [ynewq(i),xx] = liftingwave(x,1024-num,1024-2*num,1023-num);
           for m=1:1:num
               ynewq(1025-m)=xx(num-m+1);
           end
        else
           for m=1:1:1024-num
               x(m)=x(m+1);
           end  
           x(1025-num)=Signal(i);
           a=2*ynewq(i-1)-ynewq(i-2);
           for m=1:1:num
               x(1025-num+m)=a+coef*(Signal(i-m)-ynewq(i-m));
           end
          [ynewq(i),xt]=liftingwave(x,1024-num,1023,1023);
        end
    end
end