資源簡介
程序-MATLAB-小波分解(6層),通過編譯的.m文件,直接運行,希望對你有幫助

代碼片段和文件信息
function?[EjgEjfaca6faca6lcd6cd5cd4cd3cd2cd1]=computerh_zixitong(blon)
s=blon;
[cl]=wavedec(s6‘db5‘);
%%%%%%%%%%%高頻系數
ca6=appcoef(cl‘db5‘6);
%%%%%%%%%%%分解6層各層的低頻系數
cd6=detcoef(cl6);
cd5=detcoef(cl5);
cd4=detcoef(cl4);
cd3=detcoef(cl3);
cd2=detcoef(cl2);
cd1=detcoef(cl1);
%%%%%%%%%%直接利用強制削噪的方法重構信號
%cdd6=zeros(1length(cd6));
%cdd5=zeros(1length(cd5));
%cdd4=zeros(1length(cd4));
%cdd3=zeros(1length(cd3));
%cdd2=zeros(1length(cd2));
%cdd1=zeros(1length(cd1));
%c1=[ca6‘cdd6cdd5cdd4cdd3cdd2cdd1];
%s1=waverec(c1l‘db3‘);
%%%%%%%%各層的能量
Ejg=sum(ca6.^2)/length(ca6);
Ej(6)=sum(cd6.^2)/length(cd6);
Ej(5)=sum(cd5.^2)/length(cd5);
Ej(4)=sum(cd4.^2)/length(cd4);
Ej(3)=sum(cd3.^2)/length(cd3);
Ej(2)=sum(cd2.^2)/length(cd2);
Ej(1)=sum(cd1.^2)/length(cd1);
%%%%%%%%%%%%%%計算該子系統的熵
%Pj=Ej/sum(Ej);
%Hwe(i)=-sum(Pj.*log(Pj));
%hs=(blon./Hwe‘)/sum(1./Hwe‘);
%%%%%%%%%%%%計算各層分解后的小波信號的均方差
faca6=std(ca6);
fa(6)=std(cd6);
fa(5)=std(cd5);
fa(4)=std(cd4);
fa(3)=std(cd3);
fa(2)=std(cd2);
fa(1)=std(cd1);
?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件???????1141??2009-04-21?10:15??小波分解\computerh_new.m
?????文件???????2453??2009-04-24?10:57??小波分解\hzuhe_new.m
?????目錄??????????0??2009-07-16?15:38??小波分解
-----------?---------??----------?-----??----
?????????????????3594????????????????????3
評論
共有 條評論