原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
第4章 《Matlab软件基础及应用》.doc
数值计算
矩阵的LU分解、行列式、逆，线性方程组的解
4．1．1 常用算法
B=A‘ % 求矩阵的转置矩阵
det(A) %求矩阵A对应的行列式值
inv(A) %求矩阵A的逆矩阵
x=A\b %求线性方程Ax=b的解
x=b/A %求线性方程xA=b的解
[L，U，P]=lu(A) %矩阵A的上三角分解，U是上三角矩阵，L*U=P*X
[说明]
第四条指令中得斜杠”\”是”左除”符号, 对Ax=b, A在变量x得左边,所以指令中得A必须在”\”得左边,不可放错. 如要解xC=d形式得方程,则要用”右除”指令, x=d/C
【例1】“求逆”法和“左除”法解线性方程的性能对比
（1）
2368x1+0.2471x2+0.2568x3+1.2671x4=1.8471
1968x1+0.2071x2+1.2168x3+0.2271x4=1.7471
1582x1+1.1675x2+0.1768x3+0.1871x4=1.6471
1161x1+0.1254x2+0.1397x3+0.1490x4=1.5471
A=[0.2368 0.2471 0.2568 1.2671
0.1968 0.2071 1.2168 0.2271
0.1582 1.1675 0.1768 0.1871
1.1161 0.1254 0.1397 0.1490]
b=[1.8471;1.7471;1.6471;1.5471]
x=A\b
x =
1.04058380083522
0.98695649396012
0.93505250521627
88129691655365
（2）C=inv(A);
x=C*b
1.04058380083522
0.98695649396012
0.93505250521627
88129691655365
[L U P]=lu(A)
L =
1.00000000000000 0 0 0
0.14174357136457 1.00000000000000 0 0
0.17632828599588 0.16089793266405 1.00000000000000 0
0.21216736851537 0.19177989839796 0.16886625005882 1.00000000000000
U =
1.11610000000000 0.12540000000000 0.13970000000000 0.14900000000000
0 1.14972535615088 0.15699842308037 0.16598020786668
0 0 1.16690621674123 0.17412121307772
0 0 0 1.17425219838234
P =
0 0 0 1
0 0 1 0
0 1 0 0
1 0 0 0
【例2】
（1）
randn(state,0);
A=gallery(randsvd,100,2e13,2); %产生条件数为2×1013的100阶随机矩阵
x=ones(100,1); %指定真解
b=A*x; %为使Ax=b方程一致，用A和x生成b向量
cond(A) %计算条件数
ans =
1.9990e+013
（2）“求逆”法解方程的误差、残差、运算次数和所用的时间
tic %启动计时器Stopwatch Timer
xi=inv(A)*b; %xi是用“求逆”法获得的解
ti=toc %关闭计时器，并获得解方程所用的时间
eri
文档评论（0）