關(guān)于小波分析的matlab程序-資料下載頁
2025-06-16 03:46本頁面
【正文】 第六層低頻部分進(jìn)行重構(gòu)a6=wrcoef(39。a39。,c,l,39。db539。,6)。subplot(8,1,2)。plot(a6)。Ylabel(39。a639。)。%對分解結(jié)構(gòu)【c，l】中各層高頻部分進(jìn)行重構(gòu)for i=1:6 decmp=wrcoef(39。d39。,c,l,39。db539。,7i)。 subplot(8,1,i+2)。 plot(decmp)。 Ylabel([39。d39。,num2str(7i)])。end%rec=zeros(1,300)。rect=zeros(1,300)。ke=1。u=0。d1=wrcoef(39。d39。,c,l,39。db539。,1)。 figure(2)。plot(d1)。si=0。N1=0。N0=0。sce=0。for n=20:N30 rect(ke)=s(n)。 ke=ke+1。 if(ke=301) if(si==2) rec=rect。 u=2。 end。 si=0。 ke=1。 end。if(d1(n)) % the condition of abnormal append. N1=n。 if(N0==0) N0=n。 si=si+1。 end。 if(N1N0+30) Nlen=N1N0。 Tab=Nlen/10000。 end。 end。if(si==1) for k=N0:N0+99 %testing of 1/4 period signals to sce=sce+s(k)*s(k)/10000。 end。 re=sqrt(sce*200) %re indicate the pike value of . sce=0。 si=si+1。 end。end。 Nlen N0 n=1:300。 figure(3) plot(n,rec)。 基于小波變換的圖象去噪 Normalshrink算法function [T_img,Sub_T]=threshold_2_N(img,levels)% reference :image denoising using wavelet thresholding[xx,yy]=size(img)。HH=img((xx/2+1):xx,(yy/2+1):yy)。delt_2=(std(HH(:)))^2。%(median(abs(HH(:)))/)^2。%T_img=img。for i=1:levels temp_x=xx/2^i。 temp_y=yy/2^i。% belt=*(log(temp_x/(2*levels)))^。 belt=*(log(temp_x/(2*levels)))^。 % %HL HL=img(1:temp_x,(temp_y+1):2*temp_y)。 delt_y=std(HL(:))。 T_1=belt*delt_2/delt_y。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% T_HL=sign(HL).*max(0,abs(HL)T_1)。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% T_img(1:temp_x,(temp_y+1):2*temp_y)=T_HL。 Sub_T(3*(i1)+1)=T_1。 %LH LH=img((temp_x+1):2*temp_x,1:temp_y)。 delt_y=std(LH(:))。 T_2=belt*delt_2/delt_y。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% T_LH=sign(LH).*max(0,abs(LH)T_2)。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% T_img((temp_x+1):2*temp_x,1:temp_y)=T_LH。 Sub_T(3*(i1)+2)=T_2。 %HH HH=img((temp_x+1):2*temp_x,(temp_y+1):2*temp_y)。 delt_y=std(HH(:))。 T_3=belt*delt_2/delt_y。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% T_HH=sign(HH).*max(0,abs(HH)T_3)。 %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% T_img((temp_x+1):2*temp_x,(temp_y+1):2*temp_y)=T_HH。 Sub_T(3*(i1)+3)=T_3。end
語文相關(guān)推薦
鄂ICP備17016276號-1