近似熵算法

%求近似熵函数
function  apen=myAPEN3（input_signalinput_minput_r） %input_signal为要求近似熵的序列，input_m为给定的模式的维数，input_r为相似容限
x = input_signal;
r = input_r*std（x）;
[input_rowinput_col] = size（x）; %输入的信号序列一般为1行n列
n = input_col; %信号的点数
for kkk=1:2
    m = input_m+kkk-1;
    rownum = n-m+1; %构造的m维矢量的行数
    X = zeros（rownumm）; %构造信号矩阵
    for i=1:1:rownum
        for j=1:1:m
            X（ij） = x（i+j-1）;
        end
    end
    sumC（kkk） = 0;
    for ii=1:1