24小时热门版块排行榜    

Znn3bq.jpeg
查看: 1020  |  回复: 1

madfrog李

铁虫 (初入文坛)

[求助] matlab中排列熵程序出错 已有1人参与

请教,程序如下,请大神帮帮忙,
function [pe hist c] = pec(y,m,t)

%  Calculate the permutation entropy

%  Input:   y: time series;
%           m: order of permuation entropy
%           t: delay time of permuation entropy,

% Output:
%           pe:    permuation entropy
%           hist:  the histogram for the order distribution

%Ref: G Ouyang, J Li, X Liu, X Li, Dynamic Characteristics of Absence EEG Recordings with Multiscale Permutation %     %                             Entropy Analysis, Epilepsy Research, doi: 10.1016/j.eplepsyres.2012.11.003
%     X Li, G Ouyang, D Richards, Predictability analysis of absence seizures with permutation entropy, Epilepsy %     %                            Research,  Vol. 77pp. 70-74, 2007



ly = length(y);
permlist = perms(1:m);
c(1:length(permlist))=0;

for j=1:ly-t*(m-1)
     [a,iv]=sort(y(j:t:j+t*(m-1)));
     for jj=1:length(permlist)
         if (abs(permlist(jj,-iv))==0
             c(jj) = c(jj) + 1 ;
         end
     end
end

hist = c;

c=hist(find(hist~=0));
p = c/sum(c);
pe = -sum(p .* log(p));
% normalizedpe=pe/log(factorial(m));


程序如上所示,请问我在运行时出现错误
>> pec(a,2,6)
错误使用  -
矩阵维度必须一致。

出错 pec (line 25)
         if (abs(permlist(jj,-iv))==0
这是什么问题呢?
回复此楼

» 猜你喜欢

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

霜小妹二

木虫 (正式写手)

哈哈

【答案】应助回帖

首先,这个笑脸真不知道是什么鬼;其次,根据判断,你写的该是一个参数减去iv,可是,iv是向量吧,哪怕其中只有一个值!
没事儿就进行交流~
2楼2016-03-28 21:03:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 madfrog李 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿华中农业071010,320求调剂 +14 困困困困坤坤 2026-04-14 16/800 2026-04-16 22:02 by 猪会飞
[考研] 22408 312求调剂 +23 门路摸摸 2026-04-14 25/1250 2026-04-16 21:21 by Art1977
[考研] 290求调剂 +23 luoziheng 2026-04-10 25/1250 2026-04-16 17:25 by 冰雪遇见猪
[考研] 求调剂推荐 +8 小聂爱学习 2026-04-14 8/400 2026-04-16 07:22 by 学员JpLReM
[考研] 求调剂学校 +14 不会吃肉 2026-04-13 16/800 2026-04-15 21:59 by noqvsozv
[考研] 297,工科调剂? +10 河南农业大学-能 2026-04-14 10/500 2026-04-15 21:50 by noqvsozv
[考研] 一志愿A区211,22408 321求调剂 +6 随心所欲☆ 2026-04-15 7/350 2026-04-15 21:45 by lbsjt
[考研] 0854调剂 +13 长弓傲 2026-04-12 16/800 2026-04-15 13:45 by fenglj492
[考研] 复试调剂 +21 积极向上; 2026-04-10 23/1150 2026-04-15 12:50 by 西北望—风沙
[考研] 化工学硕294分,求导师收留 +32 yzyzx 2026-04-12 36/1800 2026-04-14 17:45 by lhj2009
[考研] 调剂 +12 月@163.com 2026-04-11 12/600 2026-04-14 15:37 by zs92450
[教师之家] 转长聘了 +7 简单化xn 2026-04-13 7/350 2026-04-14 14:50 by xindong
[考研] 考研求调剂 +12 子木呐 2026-04-12 13/650 2026-04-14 01:19 by 王珺璞
[考研] 326求调剂 +6 Shansyn 2026-04-10 6/300 2026-04-12 09:46 by hammer3
[考研] 一志愿西北工业大学289 085602 +33 yang婷 2026-04-10 34/1700 2026-04-12 08:11 by Art1977
[考研] 331求调剂 +5 王国帅 2026-04-11 5/250 2026-04-11 22:56 by 溪涧流水
[考研] 359求调剂 +5 胃痉挛累了 2026-04-11 5/250 2026-04-11 19:55 by lbsjt
[考研] 调剂 +5 文道星台 2026-04-11 5/250 2026-04-11 15:01 by 凯凯要变帅
[考研] 广东省 085601 329分求调剂 +14 Eddieddd 2026-04-10 14/700 2026-04-11 09:58 by bljnqdcc
[考研] 309求调剂 +14 wdhw 2026-04-10 15/750 2026-04-10 21:06 by zhouxiaoyu
信息提示
请填处理意见