matlab求解Lorenz 用matlab解下列微分方程组： dx/dy=10(x-y); dy/dt=28x-y—xz; dz/dt=xy-8/3z; 详细点~
答案:2 悬赏:20 手机版
解决时间 2021-03-02 19:52
- 提问者网友：太高姿态
- 2021-03-02 08:31
非常感谢。用ode45
最佳答案
- 五星知识达人网友：北城痞子
- 2021-03-02 09:37
%%%% 你的lorenz方程貌似不对，以下是按照我掌握的方程
function dx=ahei(t,x) % 在当前文件夹下保存为 ahei.m 文件
dx=[-8/3*x(1)+x(2)*x(3);-10*x(2)+10*x(3);-x(1)*x(2)+28*x(2)-x(3)];
-------------------------------------------------------------
运行以下：
clear;clc
t_final=100; x0=[0;0;1e-10]; % t_final为设定的仿真终止时间
[t,x]=ode45('ahei',[0,t_final],x0);
plot(t,x)
figure
plot3(x(:,1),x(:,2),x(:,3))
function dx=ahei(t,x) % 在当前文件夹下保存为 ahei.m 文件
dx=[-8/3*x(1)+x(2)*x(3);-10*x(2)+10*x(3);-x(1)*x(2)+28*x(2)-x(3)];
-------------------------------------------------------------
运行以下：
clear;clc
t_final=100; x0=[0;0;1e-10]; % t_final为设定的仿真终止时间
[t,x]=ode45('ahei',[0,t_final],x0);
plot(t,x)
figure
plot3(x(:,1),x(:,2),x(:,3))
全部回答
- 1楼网友：孤独的牧羊人
- 2021-03-02 10:40
（28+3）x（x-0.6）=28x 31x(x-0.6)=28x 31x-18.6=28x 31x-28x=18.6 3x=18.6 x=6.2
我要举报
如以上问答信息为低俗、色情、不良、暴力、侵权、涉及违法等信息，可以点下面链接进行举报！
大家都在看
推荐资讯