原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
MATLAB作业3参考解答
MATLAB作业三参考答案
请将下面给出的矩阵A与B输入到MATLAB环境中，并将它们转换成符号矩阵。若某一矩阵为数值矩阵，另以矩阵为符号矩阵，两矩阵相乘是符号矩阵还是数值矩阵。
【求解】矩阵的输入与转换是很直接的。
A=[5,7,6,5,1,6,5; 2,3,1,0,0,1,4; 6,4,2,0,6,4,4; 3,9,6,3,6,6,2;
10,7,6,0,0,7,7; 7,2,4,4,0,7,7; 4,8,6,7,2,1,7]; A=sym(A)
A =
[ 5, 7, 6, 5, 1, 6, 5]
[ 2, 3, 1, 0, 0, 1, 4]
[ 6, 4, 2, 0, 6, 4, 4]
[ 3, 9, 6, 3, 6, 6, 2]
[ 10, 7, 6, 0, 0, 7, 7]
[ 7, 2, 4, 4, 0, 7, 7]
[ 4, 8, 6, 7, 2, 1, 7]
B=[3,5,5,0,1,2,3; 3,2,5,4,6,2,5; 1,2,1,1,3,4,6; 3,5,1,5,2,1,2;
4,1,0,1,2,0,1; -3,-4,-7,3,7,8,12; 1,-10,7,-6,8,1,5]; B=sym(B)
B =
[ 3, 5, 5, 0, 1, 2, 3]
[ 3, 2, 5, 4, 6, 2, 5]
[ 1, 2, 1, 1, 3, 4, 6]
[ 3, 5, 1, 5, 2, 1, 2]
[ 4, 1, 0, 1, 2, 0, 1]
[ -3, -4, -7, 3, 7, 8, 12]
[ 1, -10, 7, -6, 8, 1, 5]
利用MATLAB语言提供的现成函数对习题1中给出的两个矩阵进行分析，判定它们是否为奇异矩阵，得出矩阵的秩、行列式、迹与逆矩阵，检验得出的逆矩阵是否正确。
【求解】以A 矩阵为例，可以对其进行如下分析。
A=[5,7,6,5,1,6,5; 2,3,1,0,0,1,4; 6,4,2,0,6,4,4; 3,9,6,3,6,6,2;
10,7,6,0,0,7,7; 7,2,4,4,0,7,7; 4,8,6,7,2,1,7]; A=sym(A);
rank(A)
ans =
7
det(A)
ans =
-35432
trace(A)
ans =
27
B=inv(A);
A*B
ans =
[ 1, 0, 0, 0, 0, 0, 0]
[ 0, 1, 0, 0, 0, 0, 0]
[ 0, 0, 1, 0, 0, 0, 0]
[ 0, 0, 0, 1, 0, 0, 0]
[ 0, 0, 0, 0, 1, 0, 0]
[ 0, 0, 0, 0, 0, 1, 0]
[ 0, 0, 0, 0, 0, 0, 1]
3、试求出习题1中给出的A与B矩阵的特征多项式、特征值与特征向量，并对它们进行LU分解。
【求解】仍以A 矩阵为例。
A=[5,7,6,5,1,6,5; 2,3,1,0,0,1,4; 6,4,2,0,6,4,4; 3,9,6,3,6,6,2;
10,7,6,0,0,7,7; 7,2,4,4,0,7,7; 4,8,6,7,2,1,7]; A=sym(A);
eig(A)
ans =
5.0093966800793665262158730069552
28.679593193974410579078264020229
.27480714110743938760483528351799e-1+1.1755376247101009492093136044131*i
-1.6336795424500642956747726147329+6.9740721596526560301948635104611*i
-3.4765922173751363914655588544224
-1.6336795424500642956747726147329-6.9740721596526560301948635104611*i
.27480714110743938760483528351799e-1-1.1755376247101009492093136044131*i
p=poly(A)
p =
x^7-27*x^6-18*x^5-1000*x^4+3018*x^3+24129*x^2+2731*x+35432
p=sym2poly(p)
p =
1 -27 -18 -1000 3018 24129 2731 35432
[L U]=lu(A)
L =
0.5000 0.5072 0.5556 0.1099 0.5000 0.5376 1.0000
0.2000 0.2319 -0.7500 0.6429 0.5702 1.
文档评论（0）