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