预处理和小波降噪

代码说明：

function z2=yuchuli(x)%% 奇异值剔除_拉依达方法%x=dlmread("55.dat");n=length(x);n1=round(n/16384);l=1;L=1;for ii=0:n1-1    for i=0:1        t1=16384*ii+8192*i+1;        t2=16384*ii+8192*i+8192;        if  t2>n     %防止采样值缺失            t1=t1-(t2-n);            t2=n;            x1=x(t1:t2);        else           x1=x(t1:t2);        end        A=mean(x1);        B=std(x1);        for num=1:8192            if  abs(x1(num)-A)>3*B                x1(num)=A;