第二章常微分方程的初值问题解析.pptVIP

练习：二阶Runge-Kutta法处理二阶常微分方程 2.4 算法的稳定性 以欧拉法的一个简单扩展为例 对微分方程进行积分时，一个首要的考虑是所用算法的数值稳定性，也就是说，舍入误差或数值计算中的其它误差能被放大的程度。 为了启动上面的递推关系，还需要 y1 的值，这可以由泰勒展开获得 当把这一方法应用到下述问题 这个问题的解析解是 y=e-x 设解为指数形式 y=Arn，代入递推关系得 正根略小于1，它对应于我们要找的按指数减小的解，但是负根略小于-1，因而它对应于一个虚假的解 上述方程的解为 其大小随n增大并且在格点上逐点发生振荡。 线性差分方程的通解正是这两个指数解的一个线性组合。虽然可以精心安排初值 y0 和 y1， 使得当 x 值小时只呈现指数衰减的解，但是在递推过程中的舍入误差将会掺入一个小的“坏”解，它最后将增长到在解中占压倒地位。 一个好的经验规则是，每当积分一个随着迭代过程急剧衰减的解时，就应当小心不稳定性和舍入误差。 例 1. 强迫钟摆 一根长度为 l 钟摆被限制在一个垂直的平面内，在强迫外力 fd 和阻力 fr 的作用下振荡运动。钟摆的运动可以通过 Newton 方程来描述 其中 fg = -mg sinθ 是重力在运动方向的分力，a=l d2 θ/dt2 是沿切线方向的加速度， θ 是杆和垂线的夹角。 2.5 动力学中的有序和混沌 那么 Newton