- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 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)