数值计算课后答案3.docVIP

  1. 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
  2. 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  3. 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  4. 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  5. 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  6. 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  7. 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
PAGE PAGE 40 习 题 三 解 答 1、用高斯消元法解下列方程组。 (1) 解:,消去第二、三个方程的,得: 再由消去此方程组的第三个方程的,得到三角方程组: 回代,得: ,, 所以方程组的解为 注意: ①算法要求,不能化简。化简则不是严格意义上的消元法,在算法设计上就多出了步骤。实际上,由于数值计算时用小数进行的,化简既是不必要的也是不能实现的。无论是顺序消元法还是选主元素消元法都是这样。 ②消元法要求采用一般形式,或者说是分量形式,不能用矩阵,以展示消元过程。 要通过练习熟悉消元的过程而不是矩阵变换的技术。 矩阵形式错一点就是全错,也不利于检查。 一般形式或分量形式: 矩阵形式 向量形式 ③必须是方程组到方程组的变形。三元方程组的消元过程要有三个方程组,不能变形出单一的方程。 ④消元顺序,不能颠倒。按为支援在方程组中的排列顺序消元也是存储算法的要求。实际上,不按顺序消元是不规范的选主元素。 ⑤不能化简方程,否则系数矩阵会变化,也不利于算法设计。 (2) 解:,消去第二、三个方程的,得: 再由消去此方程组的第三个方程的,得到三角方程组: 回代,得: , 所以方程组的解为 2、将矩阵 作LU分解。 解:设 根据矩阵乘法,先求U的第一行,由,得 。 再求L的第一列,由矩阵乘法,因为,所以,而,所以,所以。 再求U的第二行,得 ,则 , ,则 , ,则 , 再求L的第二列,得 ,则 ,则 再求U的第三行,得 ,则 ,则 再求L的第三列,得 ,则 再求U的第四行,得 ,则 所以,矩阵A的LU分解为: 指出: 用分数而表示元素,不能化成近似小数也不化成小数表示。 3、用LU分解紧凑格式分解法解方程组。 解一,用一般格式求解: 将系数矩阵作LU分解得: Ly=b方程组为 解之得 同样地,解方程组Ux=y得 。 解二,用LU紧凑格式分解法求解: 对增广矩阵三角分解: 原方程组化成同解的上三角方程组为: 回代得 。 指出: 紧凑格式是直接应用公式进行计算,计算结果保存在A的相应元素位置。从算法的角度,紧凑格式实际体现在数据的存储方法上。 由于紧凑格式计算时不再需要A的前面的元素,因此可以进行。 4、 用列主元的三角分解法解线性方程组。 解一,列选主元素消元法: 先选第一列主元为,将第一个方程与第二个方程交换,消去得: 再选第二列主元为,交换第二、三两个方程,消去得三角形方程组: 回代求得方程组的解,, 所以方程组的解为 。 解二,列主元素三角分解法: 同解的三角形方程组为 回代求得方程组的解,, 所以方程组的解为 。 说明: 用矩阵讨论中,矩阵元素进行了化简。 5.用追赶法解方程组 。 分析: 三对角矩阵 可以分解如下形式的两个矩阵: 。 即 由矩阵乘法规则,有 , 这样可以求出矩阵L和U的所有元素。 设有系数矩阵为A的方程组: , 这样的方程组称为三对角方程组。 三对角方程组经LU分解分解为 , 求解之 , 这就是所谓追赶法。 解:由公式 由此得下三角方程组 和上三角方程组 解上三角方程组 代入并解上三角方程组 6、用改进的Cholesky分解法解方程组 解:设此方程组的系数矩阵为A,右端向量为b,则 矩阵A是对称正定矩阵,可以进行乔累斯基分解。 设 由矩阵乘法得 由 得 , 再由 得 。 7、用改进的Cholesky分解法解方程组 解: 解下三角方程组 得 解上三角方程组 得 指出: 6、7两题应用一般的乔累斯基分解而没有采用书上的方法。 用MATLAB求解为: format rat a=[4,1,-1,0;1,3,-1,0;-1,-1,5,2;0,0,2,4] a = 4 1 -1 0 1 3 -1 0 -1 -1 5 2 0 0 2 4 [P,q]=chol(a) P = 2 1/2 -1/2 0 0 1985/1197 -379/838 0 0 0 1179/553 1106/1179 0 0 0 561

文档评论(0)

ma982890 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档