混沌系统 MATLAB仿真

核心问题：解微分方程组 内容来自www.wkfxw.com

以经典的Lorenz系统为例子

dxdt= -10*(x-y)

dydt= 30*x-x*z-y

dzdt= -8/3*z+x*y];

工具/原料

matlab2015B

方法/步骤

首先我们要把微分方程组用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相图

注意事项

其他混沌系统解法类似，可以私聊。