function r=hudie3(t,x)
global a;
global b;
global c;
a=8/3;b=10;c=28;
r=zeros(3,1);
r(1)=-a*x(1)+x(2)*x(3);
r(2)=-bx(2)+bx(3);
r(3)=-x(1)x(2)+cx(2)-x(3);
end
function r=hudie3(t,x)
global a;
global b;
global c;
a=8/3;b=10;c=28;
r=zeros(3,1);
r(1)=-a*x(1)+x(2)*x(3);
r(2)=-bx(2)+bx(3);
r(3)=-x(1)x(2)+cx(2)-x(3);
end
和网上一模一样的题,不过你忘记加一些运算符号了
a=8/3;b=10;c=28;
t0=[0,100];f0=[0,0,1e-10];
[t,x]=ode45('hudie3',t0,f0);
plot3(x(:,1),x(:,2),x(:,3));
title('Lorenz模型');xlabel('x');ylabel('y');zlabel('z');
grid on
function r=hudie3(t,x)
global a;
global b;
global c;
a=8/3;b=10;c=28;
r=zeros(3,1);
r(1)=-a*x(1)+x(2)*x(3);
r(2)=-b*x(2)+b*x(3);
r(3)=-x(1)*x(2)+c*x(2)-x(3);
end