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

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