近似熵的计算

我们计算两个函数，一个是henon，另外一个是TriMap：

代码：

构造henon函数：

function [x,y]=Henon(x0,y0,a,b,M)
M=M+10000;
x = zeros(1,M+1);
y = zeros(1,M+1);
x(1)=x0;
y(1)=y0;
for i = 1:M
    x(i+1) = 1+y(i)-a*x(i)^2;
    y(i+1) = b*x(i);
end
 
x=x(10001:M);
y=y(10001:M);
end

构造TriMap：

% 2-D TRIGONOMETRIC MAP
% A New Chaotic Map With Dynamic Analysis and Encryption Application in Internet of Health Things
 
function [x,y]=triMAP(x0,y0,a,b,M)
M=M+10000;
x = zeros(1,M+1);
y = zeros(1,M+1);
x(1)=x0;
y(1)=y0;
for i = 1:M
    x(i+1) = sin(a*x(i))-b*sin(a*y(i));
    y(i+1) = cos(a*x(i));
end
 
x=x(10001:M);
y=y(10001:M);
en