多尺度熵 matlab程序

%%  求解多尺度熵 Multiscale  Entropy
clc
clear all
x=random（1100）;  %假设随机信号
k=0;
t=2;% 尺度因子
N=floor（length（x）/t）;
p= 0.20;   % p  可以选择  在 0.1~0.25
r=p*std（x）
for i=1:N
    for j=1+（i-1）*t:i*t
        k=k+x（j）;
    end
    y（i）=k/t;
    k=0;
end

%%  求解m=2
m=2;
ph2=0;
num=zeros（N-m+1）;
c=zeros（N-m+1）;
for  i=1:N-m+1
    for  j=1:N-m+1
        if  j~=i
        d（ij）=max（[abs（y（i+0）-y（j+0）） abs（y（i+1）-y（j+1））]）;
        end            
    end
end
for  i=1:N-m+1
    for  j=1:N-m+1
    if d（ij）