精选优质文档-倾情为你奉上m=load(A.txt);N=length(m);for i=1:N-1 ;q(i,1)=m(i,1); end;d=q;s1=d;change=1000;c,l = wavedec(d,3,db4);%提取小波分解后的低频系数ca3=appcoef(c,l,db4,3);%提取各层小波分解后的高频系数cd3=detcoef(c,l,3);cd2=detcoef(c,l,2);cd1=detcoef(c,l,1);%对信号强制消噪cdd3=zeros(1,length(cd3);%第三层高频系数cd3全置0cdd2=zeros(1,length(cd2);%第二层高频系数cd2全置0cdd1=zeros(1,length(cd1);%第一层高频系数cd1全置0c1=ca3,cdd3,cdd2,cdd1;%构建新的系数矩阵s2=waverec(c1,l,db4)%分解新的结构%thr,sorh,keepapp=ddencmp(den,wv,d)
