4.3 一阶微分方程求解.pptVIP

4.3 一阶微分方程的求解 当两相邻离散点之间的间隔较小时，用一阶差商 取代一阶导数 前向欧拉法的几何意义： 例1. 应用前向欧拉法解初值问题 微分方程 是一阶线性微分方程，可求出其通解： 计算结果列表（ 为前向欧拉法计算近似值， 为精确值） 分析： 当步长不是很小时，前向欧拉法的精度不 是很高。步长取定后，步数越多，误差越 大。 由于前向欧拉法舍弃一阶导数以后诸项， 造成的截断误差是 的数量级，故称为 二阶精度。 二、后向欧拉法 后向欧拉法的几何意义： 例2. 应用后向欧拉法解初值问题 计算结果列表（ 为后向欧拉法计算近似值， 为精确值） 三. 梯形法及其预估-校正法 四、龙格－库塔法（R－K法） 前向欧拉法为显式的一步法，使用方便，但精度较低 四阶R-K法 * * 一阶微分方程的求解可归结为在给定初始条件下， 求微分方程的初值问题 数值解法的基本思想： 在初值问题存在唯一解的时间区间内，在若干个时间离散点上，用差分方程代替微分方程，然后逐点求解差分方程，得到各时间离散点 、 … 处的函数 近似值 、 … 一.前向欧拉法 令步长 ，则 其近似值为： 在任一步长内，用一段直线代替函数 的曲线，此直线段的斜率等于该函数在该步长起点的斜率。 取步长h=0.1，并把计算结果与精确解比较 解：据前向欧拉法 又 有： 则方程的解为： 从而有： 带入初值 可得 1.497477100 7.963873478 6.466396378 1.7 7 1.100143680 5.720961526 4.620817846 1.6 6 0.780221172 3.967666294 3.187445122 1.5 5 0.526811863 2.620359551 2.093547688 1.4 4 0.330236736 1.607215078 1.276978344 1.3 3 0.181886958 0.866642536 0.684755578 1.2 2 0.0740191694 0.345919877 0.271828183 1.1 1 0 0 0 1.0 0 n 用一阶差商近似代替 在一个步长终点的一阶导数，则原微分方程化为： 对于给定初始条件 的微分方程 其近似值： 在任一步长内，用一段直线 代替函数 的曲线，此直 线段的斜率等于该函数在该 步长终点的斜率。 取步长h=0.1，并把计算结果与精确解比较 解：据后向欧拉法 又 -1.922824061 7.963873478 9.886697539 1.7 7 -1.420624330 5.720961526 7.141585856 1.6 6 -1.013245029 3.967666294 4.980911323 1.5 5 -0.688050222 2.620359551 3.308409773 1.4 4 -0.433745534 1.607215078 2.040960612 1.3 3 -0.240212999 0.866642536 1.106855535 1.2 2 -0.098362898 0.345919877 0.444282775 1.1 1 0 0 0 1.0 0 n 用一阶差商近似地代替函数在一个步长起点和终点的 一阶导数的平均值 梯形公式 （欧拉中点公式） 近似值： 显然，梯形公式是隐式法，一般求 需要解方程，常采用迭代法，初值由显式的欧拉公式给出： 然后将 替代梯形公式等式右边出现的 当步长h足够小，且由前向欧拉法计算的已是较好 的近似，则迭代一、二次即可 例3. 应用梯形预估-校正法解初值问题 取步长h=0.1，并把计算结果与精确解比较 解：据前向欧拉法 梯形预估-校正 计算结果列表（ 为梯形预估-校正法计算近似值， 为精确值） 0.054664262 7.963873478 7.909209216 1.7 7 0.042054423 5.720961526 5.678907103 1.6 6 0.031222180 3.967666294 3.936444114 1.5 5 0.022061312 2.620359551 2.598298239 1.4 4 0.014465435 1.607215078 1.592749643 1.3 3 0.008327998 0.866642536 0.858314537 1.2 2 0.003542088 0.345919