混沌系统MATLAB仿真核心问题:解微分方程组以经典的Lorenz系统为例子dxdt=-10*(x-y)dydt=30*x-x*z-ydzdt=-8/3*z+x*y];
工具/原料
matlab2015B
方法/步骤
1、首先我们要把微抵足谛垴分方程组用MATLAB表示出来functiondydt=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系统的微分方程组。
2、用龙格库塔法解上述微分方程组,画出相图clearclcy0=[0,2,9];咯悝滩镞[t,y]=ode45('lorenz_diff',缪梨痤刻[0,200],y0);plot(y(:,1),y(:,3),'.');gridon;
3、结果展示(一)xz相图
4、结果展示(一)xy相图
5、结果展示(一)yz相图