pls 可以用于计算排列熵

可以用于计算排列熵，很实用

function [pe hist] = pec（ymt）

%  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）;
m=5
t=1
permlist = perms（1:m）;
c（1:length（permlist））=0;
    
 for j=1:ly-t*（m-1）
     [aiv]=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=c（find（c~=0））;
p = c/sum（c）;
pe = -sum（p .* log（p））;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1077  2014-04-13 09:10  pls.m

----------- ---------  ---------- -----  ----

                 1077                    1