24小时热门版块排行榜    

查看: 878  |  回复: 1

xiaochao2138

新虫 (初入文坛)

[交流] 请问有没有虫友在搞多尺度熵(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)

%% %% %% %%
回复此楼

» 猜你喜欢

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
2楼2016-02-04 15:13:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xiaochao2138 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 346求调剂[0856] +4 WayneLim327 2026-03-16 7/350 2026-03-21 04:02 by JourneyLucky
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +3 晨昏线与星海 2026-03-18 3/150 2026-03-21 00:46 by JourneyLucky
[考研] 308求调剂 +3 阿姐阿姐家啊 2026-03-18 3/150 2026-03-20 23:24 by JourneyLucky
[考研] 一志愿苏州大学材料求调剂,总分315(英一) +5 sbdksD 2026-03-19 5/250 2026-03-20 22:10 by luoyongfeng
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 求调剂 +5 Mqqqqqq 2026-03-19 5/250 2026-03-20 20:46 by zhukairuo
[考研] 295材料求调剂,一志愿武汉理工085601专硕 +5 Charlieyq 2026-03-19 5/250 2026-03-20 20:35 by JourneyLucky
[考研] 0856调剂,是学校就去 +8 sllhht 2026-03-19 9/450 2026-03-20 14:25 by 无懈可击111
[考研] 能源材料化学课题组招收硕士研究生8-10名 +5 脱颖而出 2026-03-16 14/700 2026-03-20 09:30 by kkcoco25
[考研] 0703化学调剂 +10 妮妮ninicgb 2026-03-15 14/700 2026-03-19 22:59 by 学员8dgXkO
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 304求调剂 +12 小熊joy 2026-03-14 13/650 2026-03-18 12:34 by Linda Hu
[考研] 301求调剂 +4 A_JiXing 2026-03-16 4/200 2026-03-17 17:32 by ruiyingmiao
[考研] 考研化学学硕调剂,一志愿985 +4 张vvvv 2026-03-15 6/300 2026-03-17 17:15 by ruiyingmiao
[考博] 26申博 +4 八6八68 2026-03-16 4/200 2026-03-17 13:00 by 轻松不少随
[考研] 304求调剂 +4 ahbd 2026-03-14 4/200 2026-03-16 16:48 by 我的船我的海
[考研] 277材料科学与工程080500求调剂 +3 自由煎饼果子 2026-03-16 3/150 2026-03-16 14:10 by 运气yunqi
[考研] 327求调剂 +6 拾光任染 2026-03-15 11/550 2026-03-15 22:47 by 拾光任染
[考研] 070305求调剂 +3 mlpqaz03 2026-03-14 4/200 2026-03-15 11:04 by peike
信息提示
请填处理意见