%求小波包的信息熵
%ECG為信號源，n為分解的層數，wpname為小波名字，我常用的是db3haardb4.?
function?y=waveletentropy（ECGnwpname）?
ECG=imread（‘Empreinte.bmp‘）;
n=3;
wpname=‘bior3.5‘;
wpt1=wpdec2（ECGnwpname）;?%對數據進行小波包分解?
for?i=1:2*n?%wpcoef（wpt1[ni-1]）是求第n層第i個節點的系數?
disp（‘每個節點的能量E（i）‘）;?
E（i）=norm（wpcoef（wpt1[ni-1]）2）*norm（wpcoef（wpt1[ni-1]）2）?
%求第i個節點的范數平方，其實也就是平方和?
end?
disp（‘小波包分解總能量E_total‘）;?
E_total=sum（E）?%求總能量?
disp（‘以下是每個節點的概率P‘）;?
for?i=1:2*n?
p（i）=?E（i）/E_total?%求每個節點的概率?
end?
%以下計算小波熵，即－sum（pj*lnpj），?
disp（‘小波熵的值S_wt‘）;?
for?i=1:2*n?
m（i）=p（i）*log（p（i））;?
end?
S_wt=sum（m）*（-1）?
disp（‘小波包分解圖形‘）;?
subplot（131）;?
text（1010‘原始圖象‘）?;?
plot（ECG）;?
subplot（132）;?
text（10200‘概率分布圖象‘）;?
plot（p）;?
subplot（133）;?
plot（E）;?

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????843??2009-03-16?21:05??wavelet\小波.txt

?????文件?????120054??2007-07-24?09:47??wavelet\Empreinte.bmp

?????文件????????899??2011-04-15?14:55??wavelet\waveletentropy.m

????..A.SH.??????8704??2011-04-15?19:02??wavelet\Thumbs.db

?????目錄??????????0??2011-04-15?10:45??wavelet

-----------?---------??----------?-----??----

???????????????130500????????????????????5