样本熵方法1.txt

function [shang]=ybs(xdate)

% clc,clear,

% xdate = randn(10,1);

m=2;

n=length(xdate);

r=0.2*std(xdate);

cr=[];

gn=1;

gnmax=m;

% while gn<=gnmax

x2m=zeros(n-m+1,m);%存放变换后的向量

d=zeros(n-m+1,n-m);% 存放距离结果的矩阵

cr1=zeros(1,n-m+1);%存放

k=1;

for i=1:n-m+1