程序-MATLAB-小波分解6层，通过编译，直接运行

程序-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