章小行星轨道方程计算问题——线性方程组求解的直接法.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 线性方程组直接解法概述 直接法：利用一系列递推公式计算有限步，能直接得到方程组的精确解。当然，实际计算结果仍有误差，譬如舍入误差，舍入误差的积累有时甚至会严重影响解的精度。 求解线性方程组最基本的一