MATLAB作业3参考答案
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
MATLAB 作业三参考答案
1、 请将下面给出的矩阵A 和B 输入到MA TLAB 环境中,并将它们转换成符号矩阵。若某
一矩阵为数值矩阵,另以矩阵为符号矩阵,两矩阵相乘是符号矩阵还是数值矩阵。
57651653
550123231001432546
2564206441211
346,3
9636623515
212107600774101
20
172440773473
781248672171107681
5A B ⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢
⎥==⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢
⎥---⎢⎥⎢⎥⎢⎥⎢⎥--⎣⎦⎣
⎦
【求解】矩阵的输入与转换是很直接的。
>> 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]
2、 利用MA TLAB 语言提供的现成函数对习题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.0000 0
0.6000 -0.0290 -0.9444 1.0000 0 0 0
0.3000 1.0000 0 0 0 0 0
1.0000 0 0 0 0 0 0
0.7000 -0.4203 1.0000 0 0 0 0
0.4000 0.7536 0.2778 0.6484 1.0000 0 0
U =
10.0000 7.0000 6.0000 0 0 7.0000 7.0000
0 6.9000 4.2000 3.0000 6.0000 3.9000 -0.1000