多尺度熵第三段

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)