数值分析试题及解答2009.docVIP

北航数值分析2009年试题计算题解答 四、计算 1、用二级二阶R-K方法解初值问题 时，由绝对稳定性对步长有何限制？ 解：微分方程组的Jacobi矩阵为 由解得的特征值为 当时，都是实数。对于，由绝对稳定区间，得到的取值范围 对于，又得到 易知，当时，对的限制为。 当时，，由 得 如果直接求解此不等式，将非常复杂。因此，我们采取绘图的方法。当在的范围内变化时，特征根的轨迹如图1中的红线所示。 图1　绘图分析 当时，绝对稳定条件为 （1） 式（1）表示的区域是图1中绿色曲线内部的区域。由图1可以看出，红线轨迹被完全包含在绿线围成的区域之内。因此，当时，用步长求解是绝对稳定的。当时，绝对稳定区域比图1中绿线围成的区域大，且都包含图1中绿线围成的区域（如图2所示），因此，当时，用步长求解是绝对稳定的。 图2　绘图分析（2） 综上所述，如果在整个求解区间内用不变的步长，则应使得满足 第　3　页