function?msen=MsEn（mrXnt）

%本函數計算時間序列的多尺度熵
%input
%??m：維數
%??r：閾值
%??Xn：時間序列
%??t：尺度因子
%output
%??msen：計算所得的樣本熵值

%2013.6.9??趙冰心
M=fix（length（Xn）/t）;
for?i=1:M
???????Yn（i）=?sum?（?Xn（（（（i-1）*t）+1）?:（i*t）?）?）/?t?;?
end
for?j=1:（length（Yn）-m+1）
????Ym（:j）=Yn（（（j-1）+1）:（j-1+m））;
end
for?k=1:（length（Yn）-m）
????for?g=（k+1）:（length（Yn）-m+1）
????????mm=abs（Ym（:k）-Ym（:g））;
????????d（gk）=max（mm）;
????end
end
Bm=0;

for?k=1:（length（Yn）-m）
????i=0;
????for?g=（k+1）:（length（Yn）-m+1）
????????if?d（gk）???????????i=i+1;?
???????????Bm（k）=i;
????????end
????end
end
Cmr=Bm/（M-m）;
meanCmr=mean（Cmr）;
for?j=1:（length（Yn）-m）
????Ym1（:j）=Yn（（（j-1）+1）:（j+m））;
end
for?k=1:（length（Yn）-m-1）
????for?g=（k+1）:（length（Yn）-m）
????????d1（gk）=max（abs（Ym1（:k）-Ym1（:g））

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1100??2020-07-27?23:41??MsEn.m