emd分解 边际谱 hilbert包络解调 小波包频带能量等程序

包括emd分解程序，求边际谱程序，hilbert包络解调程序，求取小波包频带能量程序，K均值聚类程序。

function [e] =energy（x0x1x2x3）

temp = wptd（x0）;
a（:1:200） = sample（temp）;
temp = wptd（x1）;
a（:201:400） = sample（temp）;
temp = wptd（x2）;
a（:401:600） = sample（temp）;
temp = wptd（x3）;
a（:601:800） = sample（temp）;
%求频带能量
e1=zeros（8800）;
for i=1:800
wpt1=wpdec（a（:i）3‘db2‘‘shannon‘）;
%d=zeros（82000）;
for j=1:2^3
    d（j:）=wprcoef（wpt1[3j-1]）;
end
for j=1:2^3
      e1（ji）=sum（abs（d（j:））.^2）;
%E1（i）=norm（D_11（i:））;
end
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%频带能量归一化
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
e=zeros（8800）;
for i=1:800
   m=sum（e1（:i））; 
    e（:i）=e1（:i）/m;
end
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1918  2012-03-30 22:42  EMD bianjipu.txt

     文件       1234  2012-07-23 09:29  EMD program.txt

     文件        699  2014-03-25 16:02  energy.m

     文件        943  2012-03-07 20:06  hilbert.txt

     文件        799  2012-04-18 09:01  Kmeans.txt

----------- ---------  ---------- -----  ----

                 5593                    5