24小时热门版块排行榜

返回列表

xiaochao2138

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 52
帖子: 6
在线: 4.2小时
虫号: 3017302
注册: 2014-03-04
专业: 光子与光电子器件

[交流] 请问有没有虫友在搞多尺度熵（MSE），有没有相关的matlab程序？

我最近在做多尺度熵，网上找到一段程序，不知道可不可用，贴一下，请各位指教：
（如果有人有不同的多尺度熵matlab程序，能不能也贴一下

谢了）
%%  求解多尺度熵 Multiscale  Entropy
clc
clear  all
load X11b.mat s
x=s(1:128,1);
k=0;
t=2;% 尺度因子,可以取1~20
N=floor(length(x)/t);
p= 0.20; % p  可以选择  在 0.1~0.25
r=p*std(x)
for i=1:N
for j=1+(i-1)*t:i*t
      k=k+x(j);
end
y(i)=k/t;
k=0;
end

%%  求解m=2
m=2;
ph2=0;
num=zeros(N-m+1);
c=zeros(N-m+1);
for  i=1:N-m+1
for  j=1:N-m+1
      if  j~=i
      d(i,j)=max([abs(y(i+0)-y(j+0)) abs(y(i+1)-y(j+1))]);
      end
end
end
for  i=1:N-m+1
for  j=1:N-m+1
if d(i,j)<r
      num(i)=num(i)+1;
end
end
c(i)=num(i)/(N-m);
ph2=ph2+c(i)
end
ph2=ph2/(N-m+1);

%%  求解m=3
m=3;
ph3=0;
num=zeros(N-m+1);
c=zeros(N-m+1);
for  i=1:N-m+1
for  j=1:N-m+1
      if  j~=i
      d(i,j)=max([abs(y(i+0)-y(j+0)) abs(y(i+1)-y(j+1)) abs(y(i+2)-y(j+2))]);
      end
end
end
for  i=1:N-m+1
for  j=1:N-m+1
if d(i,j)<r
      num(i)=num(i)+1;
end
end
c(i)=num(i)/(N-m);
ph3=ph3+c(i);
end
ph3=ph3/(N-m+1);
%%  求得多尺度熵
MSEn=log(ph2/ph3)

%% %% %% %%

回复此楼

» 猜你喜欢

假如你的研究生提出不合理要求已经有4人回复
论文终于录用啦！满足毕业条件了已经有27人回复
所感已经有3人回复
要不要辞职读博？已经有7人回复
不自信的我已经有11人回复
北核录用已经有3人回复
实验室接单子已经有3人回复
磺酰氟产物，毕不了业了！已经有8人回复
求助:我三月中下旬出站，青基依托单位怎么办？已经有10人回复
26申博（荧光探针方向，有机合成）已经有4人回复

1楼 2016-02-04 14:43:44

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

8wangle

2楼2016-02-04 15:13:45

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 xiaochao2138 的主题更新

返回列表