- 1、本文档共132页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
利用矩阵乘法,有, 即, 例3-9:用增广矩阵杜利特尔分解法求解方程组。 解: 增广矩阵的杜里特尔分解为, 即: 用LU分解求解方程组时,由于把系数矩阵的计算和对右端项的计算分开了,这使得计算方程组系非常方便。 例3-10:用杜利特尔分解法求下列矩阵的逆矩阵。 解: 求逆矩阵时,增广矩阵的杜里特尔分解为, 利用矩阵乘法,可得, 分别回代求解,得, 即, 第3章 线性方程组的数值解法 * 闫再友 第3章 线性方程组的数值解法 3.1 高斯消去法 3.2 矩阵三角分解法 3.3 平方根法 3.4 向量和矩阵的范数 3.5 方程组的形态和误差分析 3.6 迭代法 3.7 迭代法的收敛性(选) 3.2 矩阵三角分解法 3.2.1 高斯消去法的矩阵描述 矩阵三角分解法是高斯消去法的变形方法。 高斯消去法解线性方程组 先消元,然后 再回代。这些过程可以通过矩阵分解来实现。 高斯消去法用矩阵分解来描述时,是对系数矩阵分解成为一个上三角矩阵和一个下三角矩阵的乘积,即LU分解。 设线性方程组 Ax=b,系数矩阵A的各阶顺序主子式不为0,则可用顺序高斯消去法求解。 令 A(1) = A, b(1) = b 则对系数矩阵A(1)进行除不交换两行位置的初等行变换等价于用初等矩阵M1左乘A(1) ,在对方程组第1次消元后,有。 则对系数矩阵A(1)进行除不交换两行位置的初等行变换等价于用初等矩阵M1左乘A(1) ,在对方程组第1次消元后,有。 其中,M1为下三角矩阵 在对方程组第k次消元后,有 记, ,L是下三角矩阵,U是上三角矩阵。 消元过程是对k=1到n-1进行的,因此有 这样, 定义3-2:将矩阵A分解成一个下三角阵L和一个上三角阵U的乘积,即 A=LU 称为A的三角分解或LU分解。 3.2.2 矩阵的直接三角分解 将矩阵A的LU分解不唯一。 LU分解不需要消元,是直接分解 把矩阵A LU分解为一个单位下三角矩阵L或一个单位上三角矩阵U时, 这两种LU分解是唯一的。 这里A的LU分解有2种,不唯一。 例如: 一般地,任取与A同阶非奇异对角矩阵D,则对LU分解有, 杜里特尔分解 (Doolittle) 常用的两种三角分解 克洛特分解 (Crout) 定理3-5:矩阵A各阶顺序主子式不为零,则可唯一地分解成一个单位下三角阵L和一个非奇异的上三角阵U的乘积。 即: 证明:顺序高斯消去法已表明 A=LU的存在性。 再证明LU分解中,L为单位下三角阵的唯一性。 假设LU分解中,L为单位下三角阵不唯一,则有 其中,L1和L2均为单位下三角阵。 由于 |A|≠0,所以L1、L2、U1、U2均为非奇异三角阵。因此根据上式,有 该等式的左端是一个单位下三角矩阵,而右端是一个上三角矩阵,所以左右两端均为单位矩阵。唯一性得证. 例题:设 讨论a 取何值时,矩阵A可作LU分解。 解: 当A的顺序主子式不为零时,矩阵A有LU分解 所以,有 非奇异矩阵不一定存在LU分解。 例如: A是非奇异矩阵,假设存在LU分解 比较等式两端第一列,有 两式显然矛盾,所以该非奇异矩阵A不存在这样的LU分解。 设 A=LU 杜里特尔(Doolittle)分解步骤 由矩阵乘法,可得 进一步,由矩阵乘法,可得 手算杜里特尔分解的步骤 利用杜里特尔分解,并将A、L、U写在一起,有 例3-7:对下列矩阵A进行杜里特尔分解,并计算其行列式值。 解:可以直接利用式(3-27)、(3-28)、(3-29)计算 对于低价矩阵,也可以利用定义和矩阵乘法计算 设矩阵A的杜里特尔分解为, 利用矩阵乘法,有, 即, 利用手算(紧凑格式) 例题:对下列矩阵A进行杜里特尔分解 解:可以直接利用式(3-27)、(3-28)计算 例题:对下列矩阵A进行杜里特尔分解 解:可以直接利用式(3-27)、(3-28)计算 例题:对下列矩阵A进行杜里特尔分解 解:可以直接利用式(3-27)、(3-28)计算 。 设A非奇异,并有三角分解A=LU,则方程组 Ax=b 就化为 LUx=b 只须求解两个简单的三角形方程组: (1)解Ly=b 顺代求出 y (2)解Ux=y ,回代求出x. 3.2.3 用三角分解求解线性方程组 求得L、U后再求解方程组 Ly = b, Ux = y LU求解方程组的计算量和高斯消去法基本相同, 约为 例3-8:用杜利特尔分解法求解方程组(系数矩阵就是例题3-7)。 解:可以直接利
您可能关注的文档
最近下载
- 2025年长沙民政职业技术学院高职单招职业技能测试近5年常考版参考题库含答案解析.docx
- 数字医疗项目可行性报告.docx
- 110kV变电站预试定检综合项目施工专项方案.doc VIP
- 2025年21年一消防工程师继续教育题 .pdf VIP
- 2024年南昌工学院单招职业技能测试题库word版.docx VIP
- 非煤矿山露天采石场主要风险分级表.pdf VIP
- Unit 2 Making a Difference Understanding ideas The Well that changed the world 课件-2023-2024学年高中英语外研版(2019)必修第三册.pptx
- 防治责任范围矢量化操作流程.docx
- 2025学年湖南省怀化市重点中学高三5月模拟(一模)考试数学试题 .pdf VIP
- 湘少版-英语-四下-Unit1_单元测试卷.pdf
文档评论(0)