function?ApEn_value?=?ApEn（signalmr_factor）
%
%?函數名稱：ApEn（ApproximateEntropy）
%?函數功能：求信號的近似熵?
%?函數輸入：輸入時間序列???signal
%??????????模式維數???????m?:?m=1?or?m=?2?
%?函數輸出：近似熵????????ApEn_value
%?編寫作者：?Rong
%?編寫時間：?2012.10.18
%?其他說明：N?（signal?length）?between?75?and?5000;
%??????????r?between?0.1*STD?and?0.25*STD?where?STD?is?the?signal?standard?deviation
%
signal?=?signal（:）‘;
N?=?length（signal）;?
%?C?computation?for?the?“m“?pattern.
[C_m]?=?C_m_computation（signalmr_factor）;
%?C?computation?for?the?“m+1“?pattern.
[C_m_1]?=?C_m_computation（signalm+1r_factor）;
%?Phi’s?computation.
phi_m?=?mean（log（C_m））;
phi_m_1?=?mean（log（C_m_1））;
%?Final?ApEn?computation.
ApEn_value?=?[phi_m-phi_m_1];

function?[C_im]?=?C_m_computation（signalmr_factor）
X?=?[];
C_im?=?[];
n_im?=?[];
max_dif?=?[];
N?=?length（signal）;?
%?Construction?of?the?X’s?vectors.
for?j?=?1:N-m+1
????X（j:）?=?signal（j:j+m-1）;
end
%?C?computation.
for?j?=?1:N-m+1
????aux1?=?repmat（X（j:）N-m+11）;
????dif_aux?=?abs（X-aux1）;
????n_im?=?0;
????for?k?=?1:N-m+1
????????if?max（abs（dif_aux（k:）））?????????????n_im?=?n_im+1;
????????end;
????end
????C_im?=?[C_im;?n_im/（N-m+1）];
end

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

?????文件???????1268??2012-10-18?13:15??approximateEntropyAlgorithm\ApEn.m

?????文件????????501??2012-10-18?13:28??approximateEntropyAlgorithm\mainApEn.m

?????文件?????286208??2012-10-18?10:43??approximateEntropyAlgorithm\近似熵值的計算.doc

?????文件????1224704??2012-10-18?12:12??approximateEntropyAlgorithm\近似熵應用.doc

?????目錄??????????0??2012-10-18?13:30??approximateEntropyAlgorithm

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

??????????????1512681????????????????????5