matlab-小波分解與重構(gòu)-濾波

1、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);%提取小波分解后的低頻系數(shù) ca3=appcoef(c,l,db4,3);%提取各層小波分解后的高頻系數(shù)cd3=detcoef(c,l,3);cd2=detcoef(c,l,2);cd1=detcoef(c,l,1);%對信號強制消噪cdd3=zeros(1,length(cd3);% 第三層高頻系數(shù) cd3 全置 0 cdd2=zeros(1,length(cd2);% 第二層高頻系數(shù) cd2

2、 全置 0 cdd1=zeros(1,length(cd1);% 第一層高頻系數(shù) cd1 全置 0c1=ca3,cdd3,cdd2,cdd1;% 構(gòu)建新的系數(shù)矩陣 s2=waverec(c1,l,db4)%分解新的結(jié)構(gòu)%thr,sorh,keepapp=ddencmp(den,wv,d);%s2=wdencmp(gbl,c,l,db4,4,thr,sorh,keepapp);%subplot(413)%plot(1:change,s2(1:change);%title( 默認軟閾值消噪后信號)figure(1)subplot(9,2,1)plot(1:change,s1(1:change)ti

3、tle(原始信號)ylabel(s1)subplot(9,2,2)plot(1:change,s2(1:change)title(強制消噪后信號)ylabel(s2)wpt=wpdec(s1,3,db1,shannon);%plot(wpt);%重構(gòu)第三層 8 個節(jié)點小波系數(shù) s130=wprcoef(wpt,3,0);s131=wprcoef(wpt,3,1);s132=wprcoef(wpt,3,2);s133=wprcoef(wpt,3,3);s134=wprcoef(wpt,3,4);s135=wprcoef(wpt,3,5);s136=wprcoef(wpt,3,6);s137=wp

4、rcoef(wpt,3,7);%計算第三層 8 個節(jié)點小波能量 s10=norm(s130);s11=norm(s131);s12=norm(s132);s13=norm(s133);s14=norm(s134);s15=norm(s135);s16=norm(s136);s17=norm(s137);%計算方差st10=std(s130);st11=std(s131);st12=std(s132);st13=std(s133);st14=std(s134);st15=std(s135);st16=std(s136);st17=std(s137);disp(正常信號的特征向量); snorm1=s10,s11,s12,s1