| 2 | 1/1 | 返回列表 |
| 查看: 871 | 回復(fù): 1 | |||
xiaochao2138新蟲 (初入文壇)
|
[交流]
請問有沒有蟲友在搞多尺度熵(MSE),有沒有相關(guān)的matlab程序?
|
|
我最近在做多尺度熵,網(wǎng)上找到一段程序,不知道可不可用,貼一下,請各位指教: (如果有人有不同的多尺度熵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) %% %% %% %% |
|
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 一志愿中海洋材料工程專碩330分求調(diào)劑 +5 | 小材化本科 2026-03-18 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +6 | d如愿上岸 2026-03-12 | 9/450 |
|
|
[考研] 【同濟軟件】軟件(085405)考研求調(diào)劑 +3 | 2026eternal 2026-03-18 | 3/150 |
|
|
[考研] 085601材料工程專碩求調(diào)劑 +6 | 慕寒mio 2026-03-16 | 6/300 |
|
|
[考研] 268求調(diào)劑 +6 | 簡單點0 2026-03-17 | 6/300 |
|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[基金申請] 國自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進步了0608 2026-03-16 | 6/300 |
|
|
[考研] 085600調(diào)劑 +5 | 漾漾123sun 2026-03-12 | 6/300 |
|
|
[考研] 285求調(diào)劑 +6 | ytter 2026-03-12 | 6/300 |
|
|
[考研] 085601材料工程315分求調(diào)劑 +3 | yang_0104 2026-03-15 | 3/150 |
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8 | dream…… 2026-03-12 | 8/400 |
|
|
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機器學(xué)習(xí) 一志愿哈工大 +4 | 林yaxin 2026-03-12 | 4/200 |
|
|
[考研] 考研調(diào)劑 +4 | 芬達46 2026-03-12 | 4/200 |
|
|
[考研] 290求調(diào)劑 +7 | ADT 2026-03-12 | 7/350 |
|
|
[考研] 308求調(diào)劑 +3 | 是Lupa啊 2026-03-12 | 3/150 |
|
|
[考研] 一志愿山大07化學(xué) 332分 四六級已過 本科山東雙非 求調(diào)劑! +3 | 不想理你 2026-03-12 | 3/150 |
|
|
[考研] 289求調(diào)劑 +3 | 李政瑩 2026-03-12 | 3/150 |
|