题目

用matlab解Lorenz微分方程组的问题
dx/dt=a(y-x)
dy/dt=rx-y-xz
dz/dt=-bx+xy
解微分方程的时候出现问题
S=dsolve('Dx=10*(y-x)','Dy=28*x-y-x*z','Dz=-2.6*x+x*y','x(0)=1','y(0)=2','z(0)=3');
disp([S.x,S.y,S.z])
Warning:Explicit solution could not be found.
> In dsolve at 333
Access to an object's fields is only permitted within its methods.

提问时间：2021-03-13

答案

没有解析解,就计算数值解吧.function myode45[t,xyz]=ode45(@fun,0:1,[1 2 3]) % 'x(0)=1','y(0)=2','z(0)=3'); plot(t,xyz)function f=fun(t,X)x=X(1);y=X(2);z=X(3);dx=10*(y-x);dy=28*x-y-x*z;dz=-2.6*x+x*y;f=[d...