- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
赛德尔迭代法迭代法的收敛性重点:高斯消去法
计算机数学基础(下)第5编 数值分析 第10章 线性方程组的数值解法 本章主要内容: 高斯消去法 列主元消去法 雅可比迭代法 高斯—赛德尔迭代法 迭代法的收敛性 重点:高斯消去法、雅可比迭代法 难点:迭代法收敛的判定 n元线性方程组的有关概念 由n个未知量,n道所有的未知量都是一次的方程组成的方程组称为n元线性方程组。 习惯上我们用xj表示这些未知量,用aij表示它们的系数,用bj表示方程等号右边的常数。 因此,n元线性方程组的一般形式为: n元线性方程组可以写成矩阵的形式 AX=B 其中: A称为系数矩阵,是n×n矩阵,X称为未知量矩阵,b称为常数矩阵,它们都是n×1的列矩阵或称为n维列向量。 当|A|≠0时,方程组的解唯一存在。 线性方程组的增广矩阵 由系数矩阵和常数矩阵并列构成的矩阵称为线性方程组的增广矩阵,记作[A|b] [A|b]= 10.1 高斯消去法 10.1.1 高斯顺序消去法的基本思想 高斯消去法的基本思想是对线性方程组 的增广矩阵进行行初等变换,使增广矩阵中 的系数矩阵变为上三角矩阵,从而解出最后 一个未知数 xn的值,然后在依次回代,求出 其余的未知数的值。 如果消元是按方程组的自然顺序进行的 就称为高斯顺序消去法。 例1.解线性方程组 解:将方程组写成矩阵的形式 AX=b 其中: [A|b]= 它相当于方程组: 解得x3=5,回代到第2式可解得x2 = 4,再已将 求出的x3,x2的值回代到第1式可解得x1 = -1 方程组的解为(-1,4,5)T 10.1.2 高斯顺序消去法公式 记初始方程组 AX=b为 A(0)X=b(0) 第一次消元后的方程组记为A(1)X=b(1) 消元公式是: 消元后的增广矩阵是: 第二次消元后的方程组记为A(2)X=b(2) 第二次消元公式是: 由此可得, 第k次消元后的方程组记为A(k)X=b(k),公式是: 消元进行到第n-1次时结束, 可求出xn,然后回代。得到公式: 小结: 高斯顺序消去法解线性方程组的步骤: ⑴消元 ⑵对k=1到n-1,若 进行消元 第k次消元公式前面已经给出。 ⑶若 回代,方程组的解为 [定理1] 线性方程组AX=b能用高斯消去法求解的充分必要条件是A的各阶顺序主子式不为0。 例2.用高斯消去法求解线性方程组 解:写出增广矩阵 [A|b]= = [A(0)|b(0)] 第1次消元: 各元素分别加到第2、3行上,得到 [A(1)|b(1)]= 第2次消元: [A(2)|b(2)]= 此时,系数矩阵部分已化为上三角矩阵,因为 于是可以回代,得到: 所以,原方程组的解为 [2001年7月试卷计算题11] 用高斯消去法求解线性方程组 解: [A|b]= 10.1.3 列主元消去法 在高斯消去法中 我们把 分别称为第1步、第2步、 …第k步的主元,如果这些主元中的某一个主元≈0, 当它做除数时就会使舍入误差增大,导致解的严重 失真。 为了提高算法的稳定性,应选取绝对值尽可能大的 元素作为主元,这种消去法称为列主元消去法。 例3. 用列主元消去法解线性方程组 解: [A|b]= 在第1列中选取绝对值最大的元素-18作为主元, 把第1行与第2行互换后再作第1次消元 [A(0)|b(0)]= 消元后再选取第2列主元 [A(1)|b(1)]= 再选取第2列主元,把第2行与第3行互换 [A(1)|b(1)]= 再消元得 [A(2)|b(2)]= 回代得 [2002年1月试卷填空题8] 用列主元消去法解线性方程组 在作第1次消元之前,应选择主元 。 [2002年1月试卷填空题8] 用列主元消去法解线性方程组 作第1次消元之后的第3个方程是 。 解: [A|b]= 如果矩阵满足 即主对角线上每一元素的绝对值均大于同行(列)其它 元素绝对值之和,这样的矩阵称为严格对角占优矩阵。 如果矩阵满足⑴ ,⑵各阶顺 序主子式的值均为正
文档评论(0)