类Lorenz混沌系统仿真

大佬帮忙写一下仿真程序，与Lorenz类似的方程，有偿谢谢

1.首先把微分方程组用MATLAB表示出来，这是lorzen系统的微分方程组。

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)];  

2.用龙格库塔法解上述微分方程组，画出相图

clear

clc

y0 = [0,2,9];          

[t,y] = ode45('lorenz_diff',[0,200],y0); 

plot(y(:,1),y(:,3),'.');

grid on;

您好，我是有问必答小助手，您的问题已经有小伙伴解答了，您看下是否解决，可以追评进行沟通哦~

如果有您比较满意的答案 / 帮您提供解决思路的答案，可以点击【采纳】按钮，给回答的小伙伴一些鼓励哦～～

ps:问答VIP仅需29元，即可享受5次/月 有问必答服务，了解详情>>>https://vip.csdn.net/askvip?utm_source=1146287632