微分方程求解器solver说明.docVIP

微分方程求解器solver说明 表一：求解器Solver 函数指令含义函数含义solveode23普通2~3阶法解ODEodefile包含ODE的文件ode23s低阶法解刚性ODE选项odeset创建、更改Solver选项ode23t解适度刚性ODEodeget读取Solver的设置值ode23tb低阶法解刚性ODE输出odeplotODE的时间序列图ode45普通4~5阶法解ODEodephas2ODE的二维相平面图ode15s变阶法解刚性ODEodephas3ODE的三维相平面图ode113普通变阶法解ODEodeprint在命令窗口输出结果表二：不同求解器的特点 SolveODE类型特点说明ode45非刚性一步算法；4、5阶Ruhge-Kutta方程；累计截断误差达(△x)^3大部分场合的首选算法ode23非刚性一步算法；2、3阶Ruhge-Kutta方程；累计截断误差达(△x)^3适用于精度较低的情形ode113非刚性多步法；Adams算法；高低精度均可达到10^(-3)~10^(-6)计算时间比ode45短ode23t适度刚性采用梯形算法适度刚性情形ode15s刚性多步法；Gear’s反向数值微分；精度中等若ode45失败，可尝试ode23s刚性一步法；2阶Rosebrock算法；低精度精度较低时比ode15s计算时间短ode23tb刚性梯形算法；低精度同上表三：具体参数设置 属性名取值含义AbsTol有效值：正实数或向量；默认值：1e-6绝对误差对应于解向量中的所有元素；向量则分别对应于解向量中的每一分量RelTol有效值：正实数；默认值：1e-3相对误差对应于解向量中的所有元素；在每步计算过程中，误差估计为：e(k)=max(RelTol*abs(y(k)), AbsTol(k))NormControl有效值：on，off默认值：off为on时，控制解向量范数的相对误差，使每步计算中有：norm(e)=max(RelTol*norm(y), AbsTol)Events有效值：on，off返回相应的时间记录OutputFcn有效值：odeplot，odephas2，odephas3，odeprint默认值：odeplot若无输出参数，则solver执行下面操作之一：画出解向量中各元素随时间的变化；画出解向量中前两个分量构成的相平面图； 画出解向量中的前三个分量构成的三维相平面图； 随计算过程，显示解向量OutputSel有效值：正整数向量；默认值：[ ]若不使用默认设置，则属性所表现的是那些正整数指定的解向量中的分量的曲线或数据。若为默认值时，则默认地按上面情形进行操作属性名取值含义Refine有效值：正整数k1默认值：k=1若k1，则增加每个积分步中的数据点距离，使解曲线更加光滑Jacobian有效值：on，off默认值：off若为on时，返回相应的ode函数的Jacobi矩阵Jpattren有效值：on，off默认值：off返回相应的ode函数的稀疏Jacobi矩阵Mass有效值：none、M、M(t)、M(t, y)默认值：noneM：不随时间变化的常数矩阵M(t)：随时间变化的矩阵M(t,y)：随时间、地点变化的矩阵MaxStep有效值：正实数默认值：tspans/10最大积分步长