<<数值分析>> 第三章线性方程组的迭代解法

这篇具有很好参考价值的文章主要介绍了<<数值分析>> 第三章线性方程组的迭代解法。希望对大家有所帮助。如果存在错误或未考虑完全的地方,请大家不吝赐教,您也可以点击"举报违法"按钮提交疑问。

        线性方程组的理论求解公式——,在实际应用中面临着两大问题,1是计算过程复杂,2是无法保证算法的稳定性。同时初始数据存在误差,需要寻求能达到精度要求的、操作和计算过程相对简单的求解方法——迭代法。   

目录

一.迭代法的基本思想

二.基本迭代法

2.1雅可比迭代        

2.1.1三阶线性方程组的雅可比(Jacobi)迭代法

2.1.2雅可比迭代法的矩阵描述

2.2高斯-赛德尔迭代法

2.2.1三阶线性方程组的高斯-赛德尔迭代法

2.2.2高斯-赛德尔迭代法的矩阵描述

三.迭代法的收敛性

3.1迭代法的收敛性判别

3.1.1迭代矩阵的范数​编辑       

3.1.2迭代矩阵的谱半径

3.2收敛性的性质

四.参考资料


一.迭代法的基本思想

        基本思想是在逐步逼近方程组的精确解的迭代过程中获得满足精度要求的近似解。为了避免困难的矩阵求逆运算,将方程组写成等价形式。                                                                                        等价格式如下:

高斯赛德尔迭代,数值分析,算法,线性代数

        迭代格式如下:

高斯赛德尔迭代,数值分析,算法,线性代数

        对于给定的线性方程组,可以写成不同的(无穷多)迭代格式,有意义的(可用的)迭代格式应具有收敛性——生成的解向量序列收敛于方程组的解;而好的迭代法应具有较高的收敛速度。接下来介绍几种常见的迭代法。

二.基本迭代法

2.1雅可比迭代        

2.1.1三阶线性方程组的雅可比(Jacobi)迭代法

例题:

        解下列方程组,精确解x*=。

高斯赛德尔迭代,数值分析,算法,线性代数

(1)将三阶线性方程组写成等价形式:

高斯赛德尔迭代,数值分析,算法,线性代数

(2)构造迭代公式,即雅可比迭代公式:

高斯赛德尔迭代,数值分析,算法,线性代数

(3)取初始向量,即。代入迭代公式,求出:

高斯赛德尔迭代,数值分析,算法,线性代数  

(4)作为初始向量继续代入迭代公式进行迭代,求出。

(5)不断迭代下去,从表中看出迭代解不断接近精确解. 

高斯赛德尔迭代,数值分析,算法,线性代数

2.1.2雅可比迭代法的矩阵描述

        雅可比迭代公式的矩阵形式为高斯赛德尔迭代,数值分析,算法,线性代数 ,其中B是迭代矩阵。                                    将方程组的系数矩阵A,分成A=D-L-U,具体如下:

高斯赛德尔迭代,数值分析,算法,线性代数

        则方程组可以写成这种形式

高斯赛德尔迭代,数值分析,算法,线性代数

2.2高斯-赛德尔迭代法

2.2.1三阶线性方程组的高斯-赛德尔迭代法

例题:

        解下列方程组,精确解x*=。

高斯赛德尔迭代,数值分析,算法,线性代数

(1)将三阶线性方程组写成等价形式:

高斯赛德尔迭代,数值分析,算法,线性代数

(2)构造迭代公式,即高斯-赛德尔迭代公式:

高斯赛德尔迭代,数值分析,算法,线性代数

(3)取初始向量,即代入迭代公式,迭代得到的新,替代初始,代入迭代中,求出:

高斯赛德尔迭代,数值分析,算法,线性代数

(4)作为初始向量继续代入迭代公式进行迭代,迭代得到的新,替代初始,代入迭代中——求出。

(5)不断迭代下去, 从表中看出迭代解不断接近精确解. 

高斯赛德尔迭代,数值分析,算法,线性代数

2.2.2高斯-赛德尔迭代法的矩阵描述

        高斯-赛德尔迭代公式的矩阵形式为 高斯赛德尔迭代,数值分析,算法,线性代数。其中B是迭代矩阵。              同样分解A为A=D-L-U。 则方程组可以写成这种形式:

高斯赛德尔迭代,数值分析,算法,线性代数

        进一步写成迭代形式:

高斯赛德尔迭代,数值分析,算法,线性代数

        展开为:

高斯赛德尔迭代,数值分析,算法,线性代数

三.迭代法的收敛性

3.1迭代法的收敛性判别

3.1.1迭代矩阵的范数       

        迭代矩阵的范数是迭代法收敛的充分条件,迭代矩阵的范数小于1是迭代法收敛的充分条件,

3.1.2迭代矩阵的谱半径

        迭代法收敛的基本定理,迭代公式收敛的充要条件是迭代矩阵B的谱半径小于1。

例题:

        对于下面方程组,证明:雅可比迭代法收敛而高斯-赛德尔迭代法发散。

高斯赛德尔迭代,数值分析,算法,线性代数

证明 1) 对于雅可比迭代矩阵

高斯赛德尔迭代,数值分析,算法,线性代数

        的特征方程为

高斯赛德尔迭代,数值分析,算法,线性代数

         解得,为的特征根,显然,因此雅可比迭代收敛。

证明 2) 对于高斯-赛德尔迭代矩阵

        其迭代矩阵为

高斯赛德尔迭代,数值分析,算法,线性代数

         的特征方程为

高斯赛德尔迭代,数值分析,算法,线性代数

         解得,因此高斯-赛德尔迭代法发散。

3.2收敛性的性质

        对于一个线性方程组Ax=b,不同性质的A导致两种迭代法有不同的收敛。

(1)若A为严格对角占优矩阵,则解Ax = b的雅可比迭代过程和高斯-赛德尔迭代过程均收敛。      (2)若A是对称正定方阵,则解Ax=b的高斯-赛德尔迭代过程收敛。

四.参考资料

1.<<数值分析>>北京理工大学出版,感谢曾繁慧 胡行华两位老师。

2.辽宁工程技术大学,工程数学考试参考。 文章来源地址https://www.toymoban.com/news/detail-823915.html

到了这里,关于<<数值分析>> 第三章线性方程组的迭代解法的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处: 如若内容造成侵权/违法违规/事实不符,请点击违法举报进行投诉反馈,一经查实,立即删除!

领支付宝红包 赞助服务器费用

相关文章

  • <<数值分析>>第二章线性方程组的直接解法

              解线性方程组是工程数学中最常见的模型之一。所说的“最常见”有两方面的含义: 1)一部分工程问题的本身建立的就是线性方程组模型; 2)较多工程问题建立的非线性方程组模型需要转化为线性方程组的求解。           线性方程组为Ax=b,以下介绍

    2024年02月08日
    浏览(41)
  • 【数值分析不挂科】第三章 | 数值积分

    为什么要学习数值积分? 数值积分,把积分求值问题归结于被积函数值的计算,从而避开了 牛顿-莱布尼兹 公式需要寻找原函数的困难。 需要特别注意:① 区别于第二章中 n代表点的个数。**本章中的 n 指的是【区间数】**而不是点的个数!【区间数 = 点的个数 - 1】 ②所有

    2024年02月13日
    浏览(31)
  • 【数值分析实验】(五)线性方程组的迭代解法(含matlab代码)

            迭代法就是用某种极限过程去逐步逼近线性方程精确解的方法。迭代法具有需要计算机的存储单元较少、程序设计简单、原始系数矩阵在计算过程中始终不变等优点,但存在收敛性及收敛速度问题。 3.1.1 算法过程 3.1.2 代码 3.1.3 计算结果 3.2.1 算法过程 3.2.2 代码

    2024年02月03日
    浏览(32)
  • MATLAB数值分析学习笔记:线性代数方程组的求解和高斯消元法

    工程和科学计算的许多基本方程都是建立在守恒定律的基础之上的,比如质量守恒等,在数学上,可以建立起形如 [A]{x}={b} 的平衡方程。其中{x}表示各个分量在平衡时的取值,它们表示系统的 状态 或 响应; 右端向量{b}由无关系统性态的常数组成通常表示为 外部激励。 矩阵

    2023年04月15日
    浏览(53)
  • 数值分析第二次作业-求解系数矩阵为Hilbert 矩阵的线性方程组

    现要求解系数矩阵由16 阶Hilbert 方程组构成的线性方程组,右端项为  即要求解方程组Ax = b,其中 A=A0,b=b0  分别用高斯-赛德尔方法、最速下降法、共轭梯度法求解如下。 2.1. 高斯-赛德尔方法     2.2. 最速下降法    2.3. 共轭梯度法  在最速下降法中,搜索方向p取的是函数减

    2024年02月05日
    浏览(62)
  • MATLAB数值分析学习笔记:线性代数方程组的求解和高斯-赛德尔方法

    迭代法是前面介绍的消元法的有效替代,线性代数方程组常用的迭代法有 高斯-赛德尔方法 和 雅克比迭代法, 下面会讲到二者的不同之处,大家会发现两者的实现原理其实类似,只是方法不同,本篇只重点介绍高斯-赛德尔方法。 看了我之前的笔记的同学应该已经对迭代法不

    2024年02月05日
    浏览(47)
  • 数值分析·学习 | 解线性方程组的直接方法(高斯消去法以及LU求解)matlab实现

    目录 一、前言: 二、算法描述: 三、实现代码: 1、高斯消去法: 2、高斯消去法-列主元消去法: 3、LU分解: 4、求逆矩阵: 四、总结: 个人学习内容分享 1、高斯消去法:         设有线性方程组         或写为矩阵形式

    2024年02月05日
    浏览(46)
  • 【数值计算方法(黄明游)】解线性代数方程组的迭代法(一):向量、矩阵范数与谱半径【理论到程序】

       注意:速读可直接跳转至“4、知识点总结”及“5、计算例题”部分   当涉及到线性代数和矩阵理论时, 向量、矩阵范数以及谱半径 是非常重要的概念,下面将详细介绍这些内容: a. 定义及性质   考虑一个 n n n 维向量 x x x ,定义一个实值函数 N ( x ) N(x) N ( x ) ,

    2024年01月25日
    浏览(34)
  • 线性代数强化第三章

    目录 一,关于A伴随,A逆与初等矩阵  二,分块矩阵 三,矩阵方程 ​    如何证明行列式的值不能是0;    此秩为1.              法一:  法二:      不用看是列变换还是行变换,我们需要从秩不变的性质去解题                         左乘同行右乘同列,再加负号

    2024年02月11日
    浏览(29)
  • 线性代数 第三章 向量

    一、运算 加法、数乘、内积 施密特正交化 二、线性表出 概念:如果 ,则称可由线性表出(k不要求不全为0) 判定: 非齐次线性方程组 有解 无关,相关 如果两个向量组可以互相线性表出,则称这两个 向量组等价 。向量组等价,向量组的秩相等(反过来不成立,秩相等向

    2024年02月07日
    浏览(32)

觉得文章有用就打赏一下文章作者

支付宝扫一扫打赏

博客赞助

微信扫一扫打赏

请作者喝杯咖啡吧~博客赞助

支付宝扫一扫领取红包,优惠每天领

二维码1

领取红包

二维码2

领红包