【matlab混沌理论】1.2.洛伦兹吸引子

        Lorenz洛伦兹吸引子。定义洛伦兹函数组后，通过ode45函数求解此微分方程方程。

input:

% Lorenz函数的洛伦兹吸引子% 2.定义模型参数
sigma = 10; 
beta = 8/3;
rho = 28; 
% 定义一组初始条件和一个时间间隔
x0 = [1, 0, 0.5];
tspan = 0:.01:100;% 3.解常微分方程,定义一个选项结构来设置ode45的参数
options = odeset('RelTol', 1e-4, 'AbsTol', [1e-4 1e-4 1e-4]);
% ode45函数求解洛伦茨方程
[t, x] = ode45(@(t,x) LorenzFunc(t, x, sigma, rho, beta), tspan, x0, options);% 4.绘制洛伦兹吸引子
plot3(x(:,1), x(:,2), x(:,3), 'linewidth', 0.5); 
xlabel("x"); ylabel("y"); zlabel("z");
title('洛伦兹吸引子图(Lorenz Attractor)');
grid on% 1.先定义洛伦兹方程，使用ode45函数来求解常微分方程组。
function dxdt = LorenzFunc(t, x, sigma, rho, beta)dxdt = [sigma*(x(2)-x(1));      % dx/dt = sigma*(y-x)x(1)*(rho-x(3))-x(2);           % dy/dt = x*(rho-z)-yx(1)*x(2) - beta*x(3)];         % dz/dt = xy - betaz
end

output:

本文来自互联网用户投稿，文章观点仅代表作者本人，不代表本站立场，不承担相关法律责任。如若转载，请注明出处。 如若内容造成侵权/违法违规/事实不符，请点击【内容举报】进行投诉反馈！