多尺度熵第三段

function [e,A,B]=sampencp(y,M,r);%M为模版匹配数 n=length(y); LL=zeros(1,n);%zeros生成全0数组 L=zeros(1,n); A=zeros(M,1); B=zeros(M,1); p=zeros(M,1); e=zeros(M,1); for i=1:(n-1) nj=n-i; y1=y(i); for jj=1:nj j=jj+i; if abs(y(j)-y1)<r*std(y)%计算两数值之差小于规定误差 L(jj)=LL(jj)+1;

function [eAB]=sampencp（yMr）;%M为模版匹配数
n=length（y）;
LL=zeros（1n）;%zeros生成全0数组
L=zeros（1n）;
A=zeros（M1）;
B=zeros（M1）;
p=zeros（M1）;
e=zeros（M1）;
for i=1:（n-1）
   nj=n-i;
   y1=y（i）;
   for jj=1:nj
      j=jj+i;      
      if abs（y（j）-y1）         L（jj）=LL（jj）+1;
         M1=min（ML（jj））;%计算每一个数值与另一固定值之差小于误差的个数
         for m=1:M1