描述需要求解的微分方程组: 不需附加变量的格式 function xd=funname(t, x) 可以使用附加变量 function xd-funname(t, x, flag, p, p2,.) ot是时间变量或自变量(必须给),x为状态向量 X为状态向量的导数。fag用来控制求解过程,指定初 值,即使初值不用指定,也必须有该变量占位 修改变量: options唯一结构体变量,用 odeset()修改 options=odeset(RelTol,,le-7) options= odeset; options. RelTol-= le-7
描述需要求解的微分方程组： 不需附加变量的格式 function xd=funname(t,x) 可以使用附加变量 function xd=funname(t,x,flag,p1 ,p2 ,…) ％ t是时间变量或自变量（必须给），x为状态向量， xd为状态向量的导数。flag用来控制求解过程,指定初 值，即使初值不用指定，也必须有该变量占位。 修改变量：options唯一结构体变量，用odeset( )修改。 options=odeset(‘RelTol’,1e-7); options= odeset; options. RelTol= 1e-7;
例:求解 Lorenz模型的状态方程,B=8/3,p=10, G=28,初值为x1(0)=x2(0)=0,x3(0)=,∈=10-10 x1(1)=-Bx1(1)+x2(1)x3(t) px2(t)+px3(t) (1)=-x1(1)x2(1)+Ox2(1)-x3( 自变函数 function dot = lorenzeq(t, x) dot:=-8/3*x(1)+x(2)*x(3);-10*x(2)+10*x(3) X(1)*x(2)+28*x(2)-x(3)
• 例： 自变函数 function xdot = lorenzeq(t,x) xdot=[-8/3*x(1)+x(2)*x(3); -10*x(2)+10*x(3);… -x(1)*x(2)+28*x(2)-x(3)];