混沌系统 MATLAB仿真
核心问题：解微分方程组
以经典的Lorenz系统为例子
dxdt= -10*(x-y)
dydt= 30*x-x*z-y
dzdt= -8/3*z+x*y];
首先我们要把微分方程组用MATLAB表示出来
function dydt = lorenz_diff(t,y)
dydt = [ -10*(y(1)-y(2));
30*y(1)-y(1)*y(3)-y(2)
-8/3*y(3)+y(1)*y(2)];
这就是lorzen系统的微分方程组。
用龙格库塔法解上述微分方程组，画出相图
clear
clc
y0 = [0,2,9];
[t,y] = ode45('lorenz_diff',[0,200],y0);
plot(y(:,1),y(:,3),'.');
grid on;
结果展示（一）xz相图
结果展示（一）xy相图
结果展示（一）yz相图
