lorenzdifeqn1.m

来自「用小数据量方法求lyapunov指数」· M 代码 · 共 10 行

function dx = LorenzDifEqn1(t,x); 
% Lorenz系统微分方程
% 方程如下: 
% dx = -σ*(x-y)
% dy = -x*z+r*x-y
% dz = x*y-b*z
dx_1 = -(8/3)*x(1)+x(2)*x(3);
dx_2 = -10*x(2)+10*x(3);
dx_3 = -x(2)*x(1)+28*x(2)-x(3);
dx = [dx_1;dx_2;dx_3];