- 1、本文档共64页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数值分析第2章_解线性方程组的直接法
第3章 解线性方程组的直接法 LU分解求解线性方程组 小 结 3.2.2 Crout(克劳特)分解 假设 其中L为下三角矩阵, U为单位上三角矩阵, Crout Factorization 克劳特分解求方程组步骤:略 3.2.3 Cholesky(楚列斯基)分解 ——平方根法 定义3-2 若矩阵 则称A为对称正定矩阵。 A为对称正定阵,则存在唯一的LU分解 U = uij = u11 uij / uii 1 1 1 u22 unn 记为 A 对称 即 记 D1/2 = 则 仍是下三角阵 定理 设矩阵A对称正定,则存在非奇异下三角阵 使得 。若限定L对角元为正,则分解唯一。 例:对矩阵 进行楚列斯基分解。 解: A为对称矩阵,且 A对称正定。 矩阵LLT分解的实际计算公式: for for for for 解 例 解线性方程组 ? Cholesky分解法求解方程组中需说明的几个问题 ?工作量:约为LU分解的一半; ?不必选主元:A的正定性和算法的稳定性 ?稳定性:是数值稳定的; ?缺陷:存在开平方运算。 ?改进方法: LDLT 分解 改进的平方根法 改进的平方根法 ——避免开方运算 即 l. 对于 2. 引进 按行计算 元素的公式: 1. 2. 3. 对于 Step1 Step2 Step3 Stepn 元素的储存比消去法节省一半,但比平方根法多用一个单元内存, LDLT的计算顺序 LDLT分解公式: 求解方程组 等价方程组 for for for 例 用LDLT解方程组 解 例 用LDLT解方程组 解 3.3 三对角方程组的追赶法 系数矩阵为对角占优的三对角线方程组 简记为 其中,当 LU分解会是什么样子? 其中L为下三角矩阵,U为单位上三角矩阵. 由系数阵A的特点,可以将A分解为两个三角阵的乘积, 即 设 其中 为待定系数. 由矩阵乘法 或分解为:L为单位下三角矩阵,U为上三角矩阵. 计算公式 i = 2, 3, …, n-1 求解 等价于求解两个三角形方程组: 三对角线方程组的追赶法公式: 1. 计算 的递推公式 2. 解 3. 解 赶: 追: 解 例3-7 应用追赶法求解线性方程组 * 数值分析 * 第2章 解线性方程组的直接法 3.1 Gauss消去法 3.2 LU分解 3.3 三对角方程组的追赶法 Direct Methods for Solving Linear Systems 求解 其中 本章讨论n元线性方程组 的解法。 Ax=b 方程组的矩阵形式为 求解 克莱姆(Cramer法则): 所需乘除法的运算量大约为(n+1)!+n n=20时,每秒1亿次运算速度的计算机要算30多万年! 其中,D=det(A), Di 是用向量b 代替A的第i 列后所得矩阵的行列式。 所谓直接解法是指若不考虑计算过程中的舍入误差, 经过有限次算术运算就能求出线性方程组的精确解的方法。 线性方程组Ax=b的一般数值解法 1.直接法: 适用于低阶稠密方程组 Gauss消去法、三角分解法 2.迭代法:通过构造迭代方程组进行迭代。 适用于大型稀疏方程组 非零元素较多,零元素较少 上万阶,零元素很多,非零元素很少 雅克比迭代法 赛德尔迭代法 3.1 Gauss消去法 3.1.1 高斯消去法 高斯消元法是一个古老的直接法, 由它改进得到的选主元法, 是目前计算机上常用于求低阶稠密矩阵方程组的有效方法, 其特点就是通过消元将一般线性方程组的求解问题转化为三角方程组的求解问题。 /* Gaussian Elimination */ 1.消元过程:用初等行变换将原方程组的系数矩阵化为 上三角形矩阵(简称上三角阵)。 2.回代过程:对上三角形方程组的最后一个方程求解,将求得的解逐步往上一个方程代入求解。 思路 首先将A化为上三角阵/* upper-triangular matrix,再回代求解。 = 1. Gauss(顺序)消去法的基本思想 不作行交换! 例3-1 用消去法解方程组 第2步. 回代过程 解 第1步. 消元过程,对增广矩阵进行消元 得等价的三角形方程组 解为 记 Step
您可能关注的文档
- 实验3离散时间系统的频域分析.doc
- 实验报告10_记忆错觉.doc
- 实验报告实验01 C语言概述.doc
- 实验三选择结构程序设计1.doc
- 实验1-1 Windows Server 2008 IIS 70配置php,MySQL,Zend,phpMyadmin教程.doc
- 家居装修必须注意的133个细节.doc
- 宫颈炎6.ppt
- 家庭娱乐及综合信息服务需求与策略研究报告2009版.doc
- 家政管理01绪论.ppt
- 宽带连接错误代码大全600-955.doc
- 江苏省灌云高级中学2023-2024学年高三上学期第一次月考地理试题含答案.docx
- 备战2024年高考物理一轮重难点复习 第十三章 近代物理.docx
- 江苏省高邮市一中2021-2022学年高三上学期第一次月考语文试题(原卷版).docx
- 江苏省华罗庚中学2023-2024学年高三上学期12月阶段检测数学试卷含答案.docx
- 江苏省海门中学2023-2024学年高三上学期第一次调研考试地理试题及答案.docx
- 备战2024年高考物理一轮重难点复习 第九章 恒定电流.docx
- 江苏省高邮市2023-2024学年高三下学期3月学情调研数学试题含答案.docx
- 备战2024年高考物理一轮重难点复习 实验十一:探究影响感应电流方向的因素.docx
- 备战2024年高考物理一轮重难点复习 实验三:验证平行四边形定则.docx
- 备战2024年高考物理一轮重难点复习 第十四章 光.docx
文档评论(0)