3、使用泰勒公式 以此方法为基础,有龙格-库塔法、线性多步法等方 法。 4、数值公式的精度 当一个数值公式的截断误差可表示为O(h+1)时 (k为正整数,h为步长),称它是一个k阶公式 k越大,则数值公式的精度越高。 欧拉法是一阶公式,改进的欧拉法是二阶公式 龙格-库塔法有二阶公式和四阶公式。 线性多步法有四阶阿达姆斯外插公式和内插公式。 返回
3、使用泰勒公式 以此方法为基础，有龙格-库塔法、线性多步法等方 法。 4、数值公式的精度 当一个数值公式的截断误差可表示为O（h k+1）时 （k为正整数，h为步长），称它是一个k阶公式。 k越大，则数值公式的精度越高。 •欧拉法是一阶公式，改进的欧拉法是二阶公式。 •龙格-库塔法有二阶公式和四阶公式。 •线性多步法有四阶阿达姆斯外插公式和内插公式。 返 回
(三)用 Matlab软件求常微分方程的数值解 t,x]=solver (f,, ts, Xo, options) ode45 ode23 自变函数 odell3 由待解ts=[o,t,函数的 量值值 delos 方程写tt为自初值 ode23s成的m-变量的初 文件名值和终值 ode23:组合的2/3阶龙格-库塔-芬尔格算法 ode45:运用组合的45阶龙格-库塔-芬尔格算法 用于设定误差限(缺省时设定相对误差103,绝对误差106 命令为: options= odeset(’ reltol' rt.'abstol,at), rt,at:分别为设定的相对误差和绝对误差
（三）用Matlab软件求常微分方程的数值解 [t，x]=solver（’f’ ,ts,x0 ,options） ode45 ode23 ode113 ode15s ode23s 由待解 方程写 成的m- 文件名 ts=[t0，tf]， t0、tf为自 变量的初 值和终值 函数的 初值 ode23：组合的2/3阶龙格-库塔-芬尔格算法 ode45：运用组合的4/5阶龙格-库塔-芬尔格算法 自变 量值 函数 值 用于设定误差限(缺省时设定相对误差10-3 , 绝对误差10-6), 命令为：options=odeset（’reltol’ ,rt, ’abstol’ ,at）, rt，at：分别为设定的相对误差和绝对误差