%%??求解樣本熵??Sample??Entropy
function?[SampEn]=yangbenshang（xdate）
N=length（xdate）;
n=N;
x=xdate;
p=?0.15;???%?p??可以選擇??在?0.1~0.25
r=p*std（xdate）;
%%??求解m=2
m=2;
num=zeros（1N-m+1）;
x2m=zeros（n-m+1m）;%存放變換后的向量
d=zeros（n-m+1n-m）;%?存放距離結果的矩陣
??for?i=1:n-m+1??
????for?j=1:m???
????????x2m（ij）=xdate（i+j-1）;?
????end
??end
??k=1;?
??for?i=1:n-m+1?
????for?j=1:n-m+1
????????if?i~=j
????????????d（ik）=max（abs（x2m（i:）-x2m（j:）））;%計算各個元素和響應元素的距離?
????????????k=k+1;
????????end
????end
????k=1;
??end
for??i=1:N-m+1
????for??j=1:N-m
????if?d（ij）????????num（i）=num（i）+1;
????en