第五章 小行星轨道方程计算问题——线性方程组求解的直接法.docVIP

第五章 小行星轨道方程计算问题 ——线性方程组求解的直接法 5.1 小行星轨道方程问题 5.1.1 问题的引入 一天文学家要确定一颗小行星绕太阳运行的轨道，他在轨道平面内建立以太阳为原点的直角坐标系，其单位为天文测量单位，在5个不同的时间对小行星作了5次观察，测得轨道上的5个点的坐标数据如下表： 表5.1.1 轨道上的5个点的坐标数据 1 2 3 4 5 5.764 6.286 6.759 7.168 7.408 0.648 1.202 1.823 2.526 3.360 试确立小行星的轨道方程，并画出小行星的运动轨线图形。 5.1.2 模型的分析 由开普勒第一定律知，小行星轨道为一椭圆，设椭圆的一般方程为： ， 需要确定系数。 利用已知的数据，不妨设，欲确定系数等价于求解一个线性方程组： 可写成矩阵的形式： 其中， ，，。 5.1.3 模型的假设 假设： (1)小行星轨道方程满足开普勒第一定律； (2)以上所测得数据真实有效。 5.1.4 模型的建立 该问题的模型为： ， 可见，解答上述问题就是对线性方程进行求解。 5.2 线性方程组直接解法概述 直接法：利用一系列递推公式计算有限步，能直接得到方程组的精确解。当然，实际计算结果仍有误差，譬如舍入误差，舍入误差的积累有时甚至会严重影响解的精度。 求解线性方程组最基本的一种直接法是消去法。这是一个众所周知的古老方法，但用在现代电子计算机上仍然十分有效。消去法的基本思想是，通过将一个方程乘以或除以某个常数，以及将两个方程相加减这两种手续，逐步减少方程中的变元的数目，最终使每个方程仅含一个变元，从而得出所求的解。其中高斯消去法是广泛应用的方法，其求解过程分为消元过程和回代过程两个环节。消元过程将所给的方程组加工成上三角方程组。所归结的方程组再通过回代过程得出它的解。高斯消去法由于添加了回代的过程，算法结构稍复杂，但这种算法的改进明显减少了计算量。 直接法比较适用于中小型方程组。对高阶方程组，即使系数矩阵是稀疏的，但在运算中很难保持稀疏性，因而有存储量大，程序复杂等不足。 5.3 直接解法 5.3.1 高斯消去法 消去法是一个古老的求解线性方程组的方法。由它改进的选主元法是目前计算机上常用的有效的求解低阶稠密矩阵线性方程组的方法。 例5.3.1 用消去法解方程组 解 第1步，加到(5.3.2)，加到(5.3.3)，得等价方程组： 第2步，加到(5.3.5)得等价的方程组： 第3步，回代法求解(5.3.6)即可求得该方程组的解为： 用矩阵法描述的约化过程即为： . 这种求解过程称为具有回代的消去法。 此例可见消去法的基本思想是：用矩阵的初等行变换将系数矩阵化为具有简单形式的矩阵（如上三角阵，单位矩阵等），而三角形方程组是很容易回代求解的。 一般的，设有个未知数的线性方程组为： （5.3.7） 设则(5.3.7)化为为方便，。 则消去法为： 第步：计算，用乘(5.3.7)的第一方程加到第个方程中 (即实行行的初等变换) 。消去第2个到第个方程中的未知数得与(5.3.7)等价方程组： (5.3.8) 记为： 其中(5.3.8)式中元素为进一步需要计算的元素，公式为： 。 第步，继续上述过程消元。设第1步到第步计算已完成，得到与原方程组等价的方程组： (5.3.9) 记为，下面进行第步消元法： 设，计算乘数用乘 (5.3.9)中第个方程加到第个方程消去(5.3.9)中第个方程的未知数，得到与原方程组等价的方程组： (5.3.10)记为。其中中元素计算公式为： (5.3.11) 最后，重复上述过程，即且设共完成步消元计算，得到与(5.3.7)等价的三角形方程组。 (5.3.12) 再用回代法求解(5.3.12)的解，计算公式为： (5.3.13) 元素称为约化的主元素。将(5.3.7)化为(5.3.12)的过程称为消元过程。 由消元过程和回代过程求解线性方程组的方法称为消去法。(5.3.12)的求解过程(5.3.13)称为回代过程。 定理5.3.1 （Gauss消去法）设若约化的主元素则可通过Gauss消去法（不进行两行的初等变换—两行交换位置）将方程组化为等价的三角形方程组(5.3.12)。消元和求解的计算公式为： 消元计算 回代计算 5.3.2 矩阵的三角分解 下面用矩阵理论进一步来分析消去法，设约化主元素由于对实行的初等变换相当于用初等矩阵左乘，于是，消去法第1步：,有： 其中： (为