- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
Cholesky分解公式 存储??? 因对称性无需存储 Step1 Step2 Step3 Stepn 的计算过程: 逐 列 计 算 元素 仍然存放在矩阵 的相应位置上 矩阵 分解的实际计算公式(算法3.3.1): for for for for 例5:用Cholesky分 解法求解下列方程组 解: 系数矩阵为 Step1 Step2 Step3 求解方程组 求解方程组 ? Cholesky分解法求解方程组中需说明的几个问题 ?工作量:约为 分解的一半; ?不必选主元: 的正定性和算法的稳定性 ?稳定性:是数值稳定的; ?缺陷:存在开平方运算。 ?改进方法: 分解 改进的平方根法 ?改进的Cholesky法(改进的平方根法) 求解方程组 等价方程组 三对角线性方程组 其中 三. 三对角线性方程组的追赶法 给定矩阵 ,如果满足: 且 时, 则称 为上半带宽为 ,下半带宽为 的带状矩阵, 称为带状方程组; 如果 ,则称 为 的半带宽, 并称之为等带宽方程组; 为 的总带宽。 例如 上半带宽为2,下半带宽为1 总带宽为3 ?一些基本概念和理论: (保带状结构定理) 设 为上半带宽为 ,下半带宽为 的带状矩阵, 且其顺序主子式 ,则 有唯一的三角分解 , 其中 是下半带宽 为 的单位下三角阵, 是上半带宽为 的上三角阵。 证明可根据前面讲过的三角分解公式 保带状结构定理说明:矩阵的三角分解中, 和 带外元素为零,因此不必计算,且不必参加求和运算 根据保带状结构定理,系数矩阵可作如下三角分解: ?三对角矩阵 分解的计算公式: ?方程组求解的计算公式: ?解方程组 ?解方程组 “追”的过程 “赶”的过程 例2:用追赶法求解三对角方程组 , 其中: 求解方程组 求解方程组 * 因为 * 因为 * 第三章 线性方程组的直接解法 ?线性方程组: ?线性方程组的求解: Cramer法则: 所需乘除法的运算量大约为(n+1)!+n n=20时,每秒1亿次运算速度的计算机要算30多万年! 在没有舍入误差的情况下,经过有限次 运算可以得到方程组的精确解的方法。 直接法 选取初始值,逐步逼近的方法。迭代法 具有占存储单元少,程序设计简单, 原始系数矩阵在迭代过程中不变等优点, 但存在收敛性及收敛速度等问题 迭代法 基于分裂的经典迭代法; Krylov子空间迭代法。 三角形方程组 三角分解法 高斯消去法 向量和矩阵范数 线性方程组的扰动分析 本章主要内容: §3.1 三角形方程组 一、 三角形方程组的解法 ?考虑下三角形方程组 的计算公式为: 算法3.1.1 下三角形方程组的前代法: for end ?考虑上三角形方程组 的计算公式为: §3.2 高斯消去法 一、 高斯消去法 对方程组,作如下的变换,解不变 ①交换两个方程的次序 ②一个方程的两边同时乘以一个非0的数 ③一个方程的两边同时乘以一个非0数,加到另一个方程 ? 高斯消去法思想: 思路 首先将增广矩阵化为上三角阵,再回代求解 。 = 因此,对应的对增广矩阵(A,b),作如下的变换,解不变 ①交换矩阵的两行 ②某一行乘以一个非0的数 ③某一个乘以一个非0数,加到另一行 高斯消去法就是对增广矩阵作上述行的变换,变为 上三角 方程组,而后求解 设 Ax=b. 记[A b]:=[A(1) b(1)] 1、第一次消元。设 ? 高斯消去法过程: 设第k-1次消元得A(k)x=b(k) 其中 最后 Gauss消元法的可行条件为: 1. 有些方程组,不能用Gauss消元求解 2. 如果某个 很小的话,会引入大的误差 ? 高斯消元法可行条件: 当且仅当A的所有顺序主子式不为0,即 二、 列选主元高斯消去法 在Gauss消元第k步之前,做如下的事情: 若 则交换k行和j行 注: 行的交换有效地克服了Gauss消去地缺陷 例:求解方程组 请用4位浮点数计算方程组的解。 记 为下三角矩阵, 为上三角矩阵. §3.3 三角分解法 ?基本思想: ?在矩阵计算中,矩阵还有很多其它的分解,比如QR分解,奇异值分解,满秩分解等。 ?分解算法推导: 比较第1行: 比较第1列: L为单位下三角,U为上三角 一. Doolittle分解法 ?分解形式: 存储??? 比较第2行: 比较第2列: 存储??? 比
您可能关注的文档
最近下载
- 数学中考总复习.doc VIP
- 限制型心肌病超声诊断与评估.pptx
- _【课件】第四章 三角形 +问题解决策略:特殊化课件北师大版七年级数学下册.pptx VIP
- 中考数学总复习第一轮考点复习(重庆专版).pptx VIP
- 2025在线网课《信息检索与科技写作( 理大)》单元测试考核答案.pdf VIP
- 高频精选:京东快递员ai面试题及答案.doc VIP
- 新北师大版初中七年级数学下册《第四章三角形问题解决策略:特殊化(1)》教学课件.pptx VIP
- 河北 2023年农信社储蓄知识考试真题模拟汇编(共213题).doc VIP
- 河北 2023年农信社基础知识考试真题模拟汇编(共672题).doc VIP
- 四川省宜宾市2023-2024学年高一下学期期末学业质量监测语文试卷(原卷版+解析版).docx VIP
文档评论(0)