没有合适的资源?快使用搜索试试~ 我知道了~
首页小行星轨道方程问题——线性方程组直接解法
小行星轨道方程问题——线性方程组直接解法

引入实际问题 介绍线性方程组的直接解法并加入例题,能更好地应用方法去求解线性方程组
资源详情
资源评论
资源推荐

第五章 小行星轨道方程计算问题
——线性方程组的数值解法(一)
小行星轨道方程问题
问题的引入
一天文学家要确定一颗小行星绕太阳运行的轨道,他在轨道平面内建立以
太阳为原点的直角坐标系,其单位为天文测量单位,在5个不同的时间对小行星
作了5次观察,测得轨道上的5个点的坐标数据如下表:
表 轨道上的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
试确立小行星的轨道方程,并画出小行星的运动轨线图形。
模型的分析
由开普勒第一定律知,小行星轨道为一椭圆,设椭圆的一般方程为:
,需要确定系数
利用已知的数据,不妨设 欲确定系数 等价于求解一个
线性方程组:
可写成矩阵的形式:

其中,
, ,
模型的建立
该问题的模型为:
可见,解答上述问题就是对线性方程进行求解。
线性方程组直接解法概述
直接法:利用一系列递推公式计算有限步,能直接得到方程组的精确解。
当然,实际计算结果仍有误差,譬如舍入误差。舍入误差的积累有时甚至会严
重影响解的精度。
求解线性方程组最基本的一种直接法是消去法。这是一个众所周知的古老
方法,但用在现代电子计算机上仍然十分有效。消去法的基本思想是,通过将
一个方程乘或除以某个常数,以及将两个方程相加减这两种手续,逐步减少方
程中的变元的数目,最终使每个方程仅含一个变元,从而得出所求的解。其中
高斯消去法是广泛应用的方法,其求解过程分为消元过程和回代过程两个环节。
消元过程将所给的方程组加工成上三角方程组。所归结的方程组再通过回代过
程得出它的解。高斯消去法由于添加了回代的过程,算法结构稍复杂,但这种
算法的改进明显减少了计算量。
直接法比较适用于中小型方程组。对高阶方程组,即使系数矩阵是稀疏的,
但在运算中很难保持稀疏性,因而有存储量大,程序复杂等不足。
直接解法

高斯消去法
消去法是一个古老的求解线性方程组的方法。由它改进得选主元法
是目前计算机上常用的有效的求解低阶稠密矩阵线性方程组的方法。
例 用 消去法解方程组
解:第 步, 加到 , 加到 ,得等价方
程组:
第 步, 加到 得等价的方程组:
第 步,回代法求解 即可求得该方程组的解为:
用矩阵法描述的约化过程即为:
.
这种求解过程称为具有回代的 消去法。
此例可见 消去法的基本思想是:用矩阵的初等行变换将系数矩阵 化为
具有简单形式的矩阵(如上三角阵,单位矩阵等),而三角形方程组是很容易
回代求解的。
一般的,设有 个未知数的线性方程组为:
设 则 化为: 为

方便, 。则消去法为:
第 步: 计算 ,用 乘 的第一方程加
到第 个方程中去
(即实行行的初等变换) 消去第 个到第 个方程
中的未知数 得与 等价方程组:
记为:
其 中 式 中 元 素 为 进 一 步 需 要 计 算 的 元 素 , 公 式 为 :
第 步, 继续上述过程消元。设第 步到第 步计算已完成,
得到与原方程组等价的方程组:
记为 ,下面进行第 步消元法:
设 ,计算乘数 用 中第 个方
程加到第 个方程 消去 中第 个方程 的
未知数 ,得到与原方程组等价的方程组:

记为 其中 中元素计算公式为:
最后 , 重复上 述过程 , 即 且设 共完 成
步消元计算,得到与 等价的三角形方程组。
再用回代法求解 的解,计算公式为:
元素 称为约化的主元素。将 化为 的过程称为消元过程。
由消元过程和回代过程求解线性方程组的方法称为 消去法。 的
求解过程 称为回代过程。
定 理 ( 消 去 法 ) 设 若 约 化 的 主 元 素
则可通过 消去法(不进行两行的初等变换—两行交换
位置)将方程组化为等价的三角形方程组 。且消元和求解的计算公式
为:
消元计算
剩余25页未读,继续阅读


















安全验证
文档复制为VIP权益,开通VIP直接复制

评论2