求排列熵的matlab源程序？

wizer

本人需要求出手上脑电数据的排列熵，哪位大侠能赐教排列熵的matlab源程序？

gghhjj

1. function [pe hist] = pec(y,m,t)
   
2. 
   
3. %  Calculate the permutation entropy
   
4. 
   
5. %  Input:   y: time series;
   
6. %           m: order of permuation entropy
   
7. %           t: delay time of permuation entropy,
   
8. 
   
9. % Output:
   
10. %           pe:    permuation entropy
    
11. %           hist:  the histogram for the order distribution
    
12. 
    
13. %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
    
14. %     X Li, G Ouyang, D Richards, Predictability analysis of absence seizures with permutation entropy, Epilepsy %     %                            Research,  Vol. 77pp. 70-74, 2007
    
15. 
    
16. 
    
17. 
    
18. ly = length(y);
    
19. permlist = perms(1:m);
    
20. c(1:length(permlist))=0;
    
21.    
    
22. for j=1:ly-t*(m-1)
    
23.      [a,iv]=sort(y(j:t:j+t*(m-1)));
    
24.      for jj=1:length(permlist)
    
25.          if (abs(permlist(jj,:)-iv))==0
    
26.              c(jj) = c(jj) + 1 ;
    
27.          end
    
28.      end
    
29. end
    
30. 
    
31. hist = c;
    
32. 
    
33. c=c(find(c~=0));
    
34. p = c/sum(c);
    
35. pe = -sum(p .* log(p));

复制代码

wizer

gghhjj 发表于 2012-12-5 01:10

                               
登录/注册后可看大图

谢谢，我先在数据上运行试试。

aqqwbj

Multiscale Permutation Entropy (MPE)

matlabcentral/fileexchange/authors/256154

jf00332001

这个程序计算时有很多冗余，在算m=7以上的排列熵时，慢得要命

gogspe

jf00332001 发表于 2013-7-23 10:21

                               
登录/注册后可看大图

这个程序计算时有很多冗余，在算m=7以上的排列熵时，慢得要命

不对吧,大于1呢

zhangnan3509

五楼说的对，计算量有点大