版块导航: 正在加载中...

登录注册

应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜

返回列表

madfrog李

铁虫 (初入文坛)

应助: 0 (幼儿园)
金币: 41
帖子: 18
在线: 10.7小时
虫号: 4241499
注册: 2015-11-23
性别: GG
专业: 机械动力学

[求助] 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
这是什么问题呢？

回复此楼

» 猜你喜欢

1楼 2016-03-19 17:00:56

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

霜小妹二

木虫 (正式写手)

哈哈

应助: 154 (高中生)
金币: 1806.7
红花: 25
帖子: 531
在线: 66小时
虫号: 3237838
注册: 2014-05-27
专业: 机械动力学

【答案】应助回帖

首先，这个笑脸真不知道是什么鬼；其次，根据判断，你写的该是一个参数减去iv，可是，iv是向量吧，哪怕其中只有一个值！

赞一下

回复此楼

没事儿就进行交流~

2楼2016-03-28 21:03:03