图 1-7 万有引力场中质点运动轨迹
由上面例子,我们初步了解了 Matlab 解常微分方程的一般过程,首先是
建立微分方程函数文件,文件的格式如下:
fuction ydot=filename(t,y,p1,p2) %t,y 是积分区间和解矩阵 p1,p2 是参数
ydot=[关于 t,y 的表达式]; %ydot 表示 dy/dt
下面介绍 ode45 命令的用法,ode45 的一般调用格式为:
[T,Y]=ode45('fun',tspan,y0,options,p1,p2,…)
其中含义如下表:
表 1-6 ode45 命令含义
单调递增(减)的积分区间[t0:tstep:tfinal]
用 odeset 建立的优化选项,一般用默认值,
为空矢量“[ ]”
T 是输出的时 间列矢量,矩阵 Y 的每一个
列矢量是解的一个分量
各个项在命令中的位置和顺序不能颠倒,否则程序就会出错。