a.m

在Matlab环境下

function a=a(m,p)
r=0.2*variance(p);
size1=size(p);
length=size1(2);
x=zeros(length-m+1,m);
for i=1:length-m+1
    t=[];
    for j=1:m
        t=[t p(i+j-1)];
    end
    x(i,:)=t;
end
i=1;j=1;
for k=1:length-m+1
    s=0;
    for i=1:length-m+1
        if d(x(i,:),x(k,:))<r
            s=s+1;
        end
    end
    if s>1
        c(k)=(s-1)/(length-m+1);
    else
        c(k)=1/(length-m+1);
    end
end
i=1;
a=0;
for i=1:(length-m+1)
    a=a+log(c(i));
end
a=a/(length-m+1);