- 1、本文档共34页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
matlab解线性方程组 线性方程组及MATLAB应用
matlab解线性方程组 线性方程组及MATLAB应用
数值实验 线性方程组与MATLAB应用
王
1.实验目的:理解矩阵的范数与条件数。 实验内容:已知矩阵
11?1?
??11?1
A??
?1?11?
?1?1?1?1??1?
求A1,A2,A?和cond2(A)。 ?1?1??
解:编写了一个M文件来求矩阵A的范数与条件数:
test3_1.m如下:
A=[1 1 1 1;-1 1 -1 1;-1 -1 1 1;1 -1 -1 1]; norm(A,1) norm(A,2) norm(A,inf) cond(A,2)
计算结果依次是: 4 2 4 1.0000
2.实验目的:研究高斯消去法的数值稳定性(出现小主元)。
实验内容:设方程组Ax?b,其中两个矩阵如下,分别对以上两个方程组
?0.3?10?15?
5.291
(1)A1??
?11.2?
1??1??59.17?
??46.78?
?6.130?12?? ,b1??
?1?952?
???
2211????59.14
3
?7?10
??32.09999999999999
(2)A2??
?5?1?
1?001?8??
?5.9000000000?000162?? ?,b2????55?1?
???
102???
(1)计算矩阵的条件数,判断系数矩阵是良态的还是病态的?
解: 本题编写了一个test3_21的M文件如下:
A1=[0.3*1e-15 59.14 3 1;5.291 -6.130 -1 2;11.2 9 5 2;1 2 1 1]; A2=[10 -7 0 1;-3 2.099999999999999 6 2;5 -1 5 -1;0 1 0 2]; cond(A1) cond(A2)
求得两个矩阵的条件数分别为68.4296和8.9939,易知这矩阵A1的条件数远远大于1,而矩阵A2的条件数刚大于1,故这,矩阵A1为病态矩阵,矩阵A2为良态矩阵。 (2)用列主元消去法求得L和U及解向量x1,x2?R;
解:本题利用Matlab的列主元三角分解函数lu();具体求解如下: A1=[0.3*1e-15 59.14 3 1;5.291 -6.130 -1 2;11.2 9 5 2;1 2 1 1]; A2=[10 -7 0 1;-3 2.099999999999999 6 2;5 -1 5 -1;0 1 0 2];
4
b1=[59.17;46.78;1;2];
b2=[8;5.0000000000001;5;1]; [L1,U1]=lu(A1)
L1 = 0.0000 1.0000 0 0 0.4724 -0.1755 1.0000 0 1.0000 0 0 0 0.0893 0.0202 -0.1738 1.0000 U1 = 11.2000 9.0000 5.0000 2.0000 0 59.1400 3.0000 1.0000 0 0 -2.8354 1.2307 0 0 0 1.0151 [L2,U2]=lu(A2)
L2 =1.0000 0 0 0 -0.3000 -0.0000 1.0000 0 0.5000 1.0000 0 0 0 0.4000 -0.3333 1.0000 U2 =10.0000 -7.0000 0 1.0000 0 2.5000 5.0000 -1.5000 0 0 6.0000 2.3000 0 0 0 3.3667 y1=L1\b1; x1=U1\y1 x1 =3.8457 1.6095 -15.4761 10.4113 y2=L2\b2; x2=U2\y2 x2 =0.1337 -0.8218 0.8842 0.9109
?1,x?2?R4x用不选主元的高斯消去法求得L和U及解向量;
解:编写一个LU_Fact的M文件储存不选主元的LU分解法然后调用求解。具
您可能关注的文档
- dnf深渊派对挑战书有什么用.doc
- dnf璀璨的挑战书礼盒.doc
- dnf璀璨的挑战书礼箱.doc
- dnf炒挑战书.doc
- dnf站街送挑战书.doc
- dnf邀请函和挑战书.doc
- dnf邀请函挑战书.doc
- doorman说明书.doc
- dreamweaver毕业设计.doc
- dreamweaver毕业设计_0.doc
- 上课《归园田居(其一)》课件28张.pptx
- 色彩的魅力课件55.ppt
- 山水田园诗导学练习鉴赏.ppt
- 部编版《白鹭》课件.pptx
- 《Excel在财务中的应用》项目七 全面预算编制方法与应用 教学课件.pptx
- 《Excel在财务中的应用》项目六 利润管理决策模型与应用 教学课件.pptx
- 《Excel在财务中的应用》项目五 运营资金的管理 教学课件.pptx
- 《Excel在财务中的应用》项目四 项目投资管理决策与应用 教学课件.pptx
- 《Excel在财务中的应用》项目三 筹资管理决策模型与应用 教学课件.pptx
- 《Excel在财务中的应用》项目二 Excel 常用函数基础 教学课件.pptx
文档评论(0)