第五章科学计算与数学建模.pptVIP

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