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

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

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

一.高斯列主元消去法

1.1介绍

1.2例题

1.3特点

 二.LU分解求解方程组

2.1公式介绍

2.2求解思路

2.3例题

三.特殊的LU分解

3.1平方根法

3.2Cholesky分解

3.2.1方法介绍

3.2.2例题

3.3改进的平方根法

3.3.1方法介绍

3.3.2分解过程

3.3.3例题

四.向量和矩阵的范数

4.1向量的范数

 4.2矩阵的范数

4.2.1 例题:

 五.矩阵的条件数与误差分析

 5.1误差原因

5.1.1病态矩阵介绍

5.2线性方程组的误差分析

5.2.1 b有误差而A无误差的情形

5.2.2条件数的性质

5.2.3例题:

六.参考资料

一.高斯列主元消去法

1.1介绍

        方程需要有唯一解,并且不接近不接近奇异矩阵。高斯列主元素消去法就是在消元的每一步选取(列)主元素——列中绝对值最大的元(数)取作主元素。在每轮消元前,根据需要消去的行,确定消元因子Lij(小于1的数)。以下是运算的步骤:

(1)选列主元素:从最左侧列开始选,同一列中最大值。                                                                (2)交换两行:是列主元所在的行为第一行,若不是,则交换两行。                                                (3)消元运算:用列主元将每一列的其余项消为0。                                                                        (4)重复以上步骤                                                                                                                            (5)代入求解:得到一个x的解,代入其他行,求解得到其他的x解。

1.2例题

用高斯列主元消去法解方程组:

lu分解例题及解析,数值分析,算法,线性代数,矩阵

解:

lu分解例题及解析,数值分析,算法,线性代数,矩阵

         易得方程组的解为 

1.3特点

        高斯列主元消去法是数值稳定的方法。

 二.LU分解求解方程组

2.1公式介绍

        高斯消去法的过程,可以看作下三角矩阵左乘系数矩阵A,乘积为可逆的上三角矩阵U。系数方程组为Ax=b,将Ux看作y,先解y,再解x。则有以下的公式:

lu分解例题及解析,数值分析,算法,线性代数,矩阵

        其中A为系数矩阵,L为单位的下三角矩阵,U为可逆的上三角矩阵

lu分解例题及解析,数值分析,算法,线性代数,矩阵

2.2求解思路

(1)首先确定A分成的L和U。由 A 的第一行、第一列元素确定 U 的第一行、L 的第一列元素。再确定U的下一行,L的下一列。可以确定全部的L,U值。                                                                    (2)解方程组Ly=b,求解y的值,再解方程组Ux=y,求解x的值。

2.3例题

用LU分解解方程组

lu分解例题及解析,数值分析,算法,线性代数,矩阵

(1)对A进行LU分解:A=LU

lu分解例题及解析,数值分析,算法,线性代数,矩阵

        易得y的解为,x的解为

三.特殊的LU分解

3.1平方根法

        平方根法是解对称正定方程组的有效方法,系数矩阵A分解为L和U,再将U分解成D和Uo。可以得到A=L·D·Uo 。                                                                                                                                   具体分解如下:

lu分解例题及解析,数值分析,算法,线性代数,矩阵

        对称正定矩阵A有唯一的分解,D是正定对角阵——对角矩阵只在对角线上含有非0元素,其它位置都为0。D的对角元素Uij,记作,其中为下图矩阵

lu分解例题及解析,数值分析,算法,线性代数,矩阵

        则

3.2Cholesky分解

3.2.1方法介绍

        将上述中记作L,则被称为Cholesky分解。利用 Cholesky 直接分
解公式,推导出的解方程组方法,称为 Cholesky 方法或平方根法。

求解思路:

(1)首先将对称正定矩阵A分解成。                                                                                      (2)根据矩阵乘法,逐步确定L的第i行元素,得到和。                                                          (3)方程组转换为。解,得y的解,解,得x的解。

3.2.2例题

用Cholesky 方法解方程组

lu分解例题及解析,数值分析,算法,线性代数,矩阵

 解:对系数矩阵A进行分解得到

lu分解例题及解析,数值分析,算法,线性代数,矩阵

         解,得,解,得。

3.3改进的平方根法

3.3.1方法介绍

        用平方根法解对称正定方程组时,求解L时用到开方的运算。为避免开方运算,改进平方根法,得到改进的平方根法——其适用于A是对称正定,或A对称且顺序主子式全不为0的情况。

3.3.2分解过程

        将A分解为,L是单位下三角矩阵,D是对角矩阵,分解的L、D如下

lu分解例题及解析,数值分析,算法,线性代数,矩阵

求解思路

(1)首先将对称正定矩阵A分解为。                                                                              (2)有Ax=b,x=b。首先解Ly=b,解得y,再解x=y,解得x。

3.3.3例题

用改进的平方根法解方程组

lu分解例题及解析,数值分析,算法,线性代数,矩阵

四.向量和矩阵的范数

        为了对线性方程组的直接法做出误差分析和讨论方程组迭代法的收敛性。引入向量空间中向量及向量序列极限的大小。范数是用来度量向量的大小的。

4.1向量的范数

        对于一个向量x,其满足条件,

lu分解例题及解析,数值分析,算法,线性代数,矩阵

(1)1—范数:值为所有向量的绝对值之和。

lu分解例题及解析,数值分析,算法,线性代数,矩阵

(2)2—范数:值为向量和他的转置向量乘积,再开根号。

lu分解例题及解析,数值分析,算法,线性代数,矩阵

(3)—范数:值为向量中的最大值

lu分解例题及解析,数值分析,算法,线性代数,矩阵

 4.2矩阵的范数

(1)1—范数:值为所有的列向量组中绝对值和最大值。

lu分解例题及解析,数值分析,算法,线性代数,矩阵

(2)2—范数:值为A的转置和A的乘积的矩阵的最大的特征值,再开根号。

lu分解例题及解析,数值分析,算法,线性代数,矩阵

(3) —范数:值为所有的行向量组中绝对值和的最大值。

lu分解例题及解析,数值分析,算法,线性代数,矩阵

(4)F—范数:值为矩阵中所有值的平方和,再开根号。

lu分解例题及解析,数值分析,算法,线性代数,矩阵

4.2.1 例题:

 求矩阵的范数:

lu分解例题及解析,数值分析,算法,线性代数,矩阵

解得: 

lu分解例题及解析,数值分析,算法,线性代数,矩阵

 五.矩阵的条件数与误差分析

 5.1误差原因

直接法解线性方程组产生误差的主要原因:

(1)算法及舍入误差的影响。

(2)方程组本身存在病态,原始数据扰动对解影响。

5.1.1病态矩阵介绍

        如果矩阵A或者常数项b的微小变化,引起了方程组Ax=b的解——x的巨大变化,称此方程组为病态方程组,矩阵A相对于方程组为病态矩阵,反之称为良态方程组、良态矩阵。

5.2线性方程组的误差分析

5.2.1 b有误差而A无误差的情形

        将带有误差的右端项b和带误差的解向量x*代入方程组Ax=b,则有等式:

lu分解例题及解析,数值分析,算法,线性代数,矩阵

        因为 x*为精确值,有A x*=b,则有 。不等式1,取范数为,同时Ax=b,取范数lu分解例题及解析,数值分析,算法,线性代数,矩阵,即不等式2lu分解例题及解析,数值分析,算法,线性代数,矩阵 。不等式1比上不等式2可得

误差估计式:

lu分解例题及解析,数值分析,算法,线性代数,矩阵

        其中方阵A的条件数。

5.2.2条件数的性质

        条件数很大的矩阵称为“病态”矩阵;病态矩阵对应的方程组称为病态方程组。反之,则称矩阵为良态矩阵,对应的方程组为良态方程组。

(1)解的相对误差是右端项b的相对误差的 cond(A)倍;                                                                  (2)如果条件数越大,则解的相对误差就可能越大;
(3)条件数成了刻划矩阵的病态程度和方程组解对A或b扰动的敏感程度。

5.2.3例题:

        已知方程组Ax=b,b有扰动,估计解的相对误差。

lu分解例题及解析,数值分析,算法,线性代数,矩阵

根据误差估计式:

lu分解例题及解析,数值分析,算法,线性代数,矩阵

        分别求解右侧的值cond(A),,。并求解其无穷范数,代入不等式得:

lu分解例题及解析,数值分析,算法,线性代数,矩阵

         解的相对误差相比右端项b的相对误差相差了2015倍。文章来源地址https://www.toymoban.com/news/detail-720336.html

六.参考资料

1.<<数值分析>>北京理工大学出版,感谢曾繁慧 胡行华两位老师
2.辽宁工程技术大学,工程数学考试重点。
3.感谢 鞠采洋对我的支持。

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

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

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

相关文章

  • MATLAB数值分析学习笔记:线性代数方程组的求解和高斯消元法

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

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

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

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

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

    2024年02月05日
    浏览(79)
  • 从零开始学数据分析之——《线性代数》第二章 矩阵

    元素全为实数的矩阵称为实矩阵  元素全为负数的矩阵称为复矩阵 只有一行(列)的矩阵称为行(列)矩阵 元素全为零的矩阵称为零矩阵 行数和列数都等于n的矩阵称为n阶矩阵或n阶方阵 主对角线元素全为1,其余元素全为0的矩阵称为单位矩阵,记作E或I 两个矩阵行数和列数

    2023年04月23日
    浏览(49)
  • 线性代数 第二章 矩阵

    一、概念 个数排成的m行n列的表格 二、运算法则 三、初等变换 (1)用非零常数k乘矩阵的某一行(列); (2)互换矩阵某两行(列)的位置; (3)把某行(列)的k倍加至另一行(列)。 称为矩阵的 初等行(列)变换 ,统称 初等变换 。矩阵经初等行变换后秩不变。 初等

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

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

    2024年01月25日
    浏览(47)
  • 高等数学:线性代数-第二章

    n bm{n} n 元线性方程组 设有 n 个未知数 m 个方程的线性方程组 { a 11 x 1 + a 12 x 2 + ⋯ + a 1 n x n = b 1 a 21 x 1 + a 22 x 2 + ⋯ + a 2 n x n = b 2 ⋯ ⋯ ⋯ ⋯ a m 1 x 1 + a m 2 x 2 + ⋯ + a m n x n = b m begin{cases} a_{11}x_{1} + a_{12}x_{2} + cdots + a_{1n}x_{n} = b_{1} \\\\ a_{21}x_{1} + a_{22}x_{2} + cdots + a_{2n}x_{n} = b

    2024年02月11日
    浏览(42)
  • 【数据结构】第二章——线性表(4)

    大家好,很高兴又和大家见面啦!!! 在前面的内容中我们介绍了线性表的第一种存储方式——顺序存储,相信大家经过前面的学习应该已经掌握了对顺序表的一些基本操作了。今天,我们将开始介绍线性表的第二种存储方式——链式存储。 线性表中的数据元素在存储时,

    2024年02月04日
    浏览(49)
  • 【数据结构】第二章——线性表(3)

    大家好,很高兴又和大家见面了!!! 在上一篇中,咱们介绍了顺序表的基本概念,以及通过C语言实现顺序表的创建和对表长的修改。今天咱们将详细介绍一下使用C语言实现顺序表的增删改查。接下来,跟我一起来看看今天的内容吧!!! 我们先来回顾一下上一篇的内容,

    2024年02月04日
    浏览(54)
  • 【数据结构】第二章——线性表(2)

    大家好,很高兴又和各位见面啦!!!在上一个篇章中,我们简单了解了一下线性表的基础知识以及一下重要的术语。在今天的篇章中我们将来开始正式介绍线性表的顺序存储——又称顺序表。我们将会在本章介绍什么是顺序表,对于顺序表的操作我们又应该如何实现。接下

    2024年02月03日
    浏览(50)

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

支付宝扫一扫打赏

博客赞助

微信扫一扫打赏

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

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

二维码1

领取红包

二维码2

领红包