实验指导3 矩阵的基本运算

  • 格式:pdf
  • 大小:77.31 KB
  • 文档页数:6

实验指导3 矩阵的基本运算
一、实验目的
熟悉Matlab软件中关于矩阵运算的各种命令
二、实验内容
已知矩阵A、B、b如下:

b=[1 3 5 7 9 11 ]
应用Matlab软件进行矩阵的输入及各种基本操作。

三、实验仪器和设备
1.计算机若干台(装有matlab6.5及以上版本软件)
2.打印机
四、实验要求
1.独立完成各个实验任务;
2.实验的过程保存成 .m 文件,以备检查;
3.实验结果保存成 .mat 文件
五、实验原理
(一)预备知识
z线性代数中的矩阵运算;
z本实验所用Matlab命令提示:
(1)矩阵的输入格式:;b=初始值 :步长 :终值; (2)求A的转置:A ’;
(3)求A加B:A+B;
(4)求A减B:A-B;
(5)求数k乘以A:k*A;
(6)求A乘以B:A*B;
(7)求A的行列式:det(A);
(8)求A的秩:rank(A);
(9)求A的逆:inv(A);
(10)B右乘A的逆: B / A;
(11)B左乘A的逆:A\B;
(12)求A的特征值:eig(A);
(13)求A的特征向量矩阵X及对角矩阵D:[X,D]=eig(A);
(14)方阵A的n次幂: A^ n ;
(15)A与B的对应元素相乘:A .* B;
(16)存储工作空间变量:save’文件名’’变量名’;
(17)列出工作空间所有变量:whos
(二)实验内容及要求
1.输入矩阵A、B、b;
2.作X21=A ’ 、X22=A+B 、X23=A-B 、X24=A*B;
3.作X31=|A| 、X32=|B|;
4.作X41=R(A) 、X42=R(B);
5.作X5=;
6.求满足矩阵方程:XA=C的解矩阵X6。

其中C为A的第i列乘i 所得矩阵;
7.求满足方程组:AX=b ’ 的解向量X7 ;
8.作X6的特征值向量X8 、X6 的特征向量组 X 及对角阵D ;
9.作X9=;
10.存储工作空间变量 A 、B 。

(三)操作提示
1.计算过程(程序名sy1.m)
A=[3 4 -1 1 -9 10;6 5 0 7 4 -16;1 -4 7 -1 6 -8;...
2 -4 5 -6 12 -8; -
3 6 -7 8 -1 1;8 -
4 9 1 3 0]
B=[1 2 4 6 -3 2;7 9 16 -5 8 -7;8 11 20 1 5 5;...
10 15 28 13 -1 9;12 19 36 25 -7 23;2 4 6 -3 0 5]
b=1:2:11
x21=A'
x22=A+B
x23=A-B
x24=A*B
x31=det(A)
x32=det(B)
x41=rank(A)
x42=rank(B)
x5=inv(A)
for i=1:6
c(:,i)=i*A(:,i); end
c
x6=c/A
x7=A\b'
x8=eig(x6)
[x,D]=eig(x6)
x9=B^2*(A^(-1))^2 save sy1 A B
2.计算结果
A =
3 4 -1 1 -9 10
6 5 0
7 4 -16
1 -4 7 -1 6 -8
2 -4 5 -6 12 -8 -
3 6 -7 8 -1 1
8 -4 9 1 3 0
B =
1 2 4 6 -3 2
7 9 16 -5 8 -7
8 11 20 1 5 5
10 15 28 13 -1 9 12 19 36 25 -7 23 2 4 6 -3 0 5
b =
1 3 5 7 9 11
x21 =
3 6 1 2 -3 8
4 5 -4 -4 6 -4
-1 0 7 5 -7 9
1 7 -1 -6 8 1
-9 4 6 12 -1 3
10 -16 -8 -8 1 0
x22 =
4 6 3 7 -12 12
13 14 16 2 12 -23 9 7 27 0 11 -3
12 11 33 7 11 1
9 25 29 33 -8 24
10 0 15 -2 3 5
x23 =
2 2 -5 -5 -6 8
-1 -4 -16 12 -4 -9
-7 -15 -13 -2 1 -13
-8 -19 -23 -19 13 -17
-15 -13 -43 -17 6 -22
6 -8 3 4 3 -5
x24 =
-55 -85 -180 -245 80 -175
127 174 348 250 -13 52
75 110 220 194 -41 154
82 129 260 283 -91 239
53 76 138 21 21 -29
98 151 284 165 -33 167
x31 =
245295
x32 =
x41 =
6
x42 =
4
x5 =
-0.0737 0.0604 -0.2297 0.0067 -0.0804 0.1042 0.3142 0.0036 0.2408 0.1605 0.1259 -0.1436 0.2099 -0.0395 0.3155 0.0364 0.0834 -0.0663 -0.0827 -0.0123 0.0088 -0.0777 0.0779 0.0878 0.0134 -0.0335 -0.0159 0.1129 0.1061 0.0337 0.0377 -0.0525 -0.0110 0.0469 0.0698 0.0411
c =
3 8 -3
4 -4
5 60
6 10 0 28 20 -96
1 -8 21 -4 30 -48
2 -8 15 -24 60 -48
-3 12 -21 32 -5 6
8 -8 27 4 15 0
x6 =
2.9937 -1.3602 0.3790 -1.3862 0.2409 0.6621 -2.9686 4.4209 2.0166 -2.7720 -1.6217 -1.6197
0.7423 0.7632 4.4879 0.9367 0.1844 -1.4517
1.4636 0.3018 1.7126 5.6670 1.2292 -1.6950
-2.9029 0.1517 -2.7542 -1.6300 2.3827 2.2451 2.4342 -1.1655 4.5517 1.1351 2.5048 1.0478
x7 =
-0.5714
2.2062
1.9446
1.0471
1.9493
1.2331
x8 =
1.0000
2.0000
6.0000
3.0000
4.0000
5.0000
x =
-0.2705 0.3578 -0.4541 -0.0698 0.0811 -0.5313 -0.5410 0.4472 0.7265 -0.0000 0.5678 0.2361 -0.0902 -0.3578 0.3633 0.4889 -0.0811 0.3542 -0.1803 -0.3578 0.3633 0.3492 -0.4867 0.7083 0.2705 0.5367 -0.0454 -0.4889 0.6489 -0.0590 -0.7213 -0.3578 0.0000 0.6286 0.0811 0.1771
D =
1.0000 0 0 0 0 0
0 2.0000 0 0 0 0
0 0 6.0000 0 0 0
0 0 0 3.0000 0 0
0 0 0 0 4.0000 0
0 0 0 0 0 5.0000
x9 =
4.8329 1.1095 4.3899 2.1612 3.5801 6.1010
16.7262 1.3541 15.0714 9.5847 13.0914 18.5071
22.0322 3.1761 19.5899 11.9050 17.4229 27.6293
31.6981 5.3951 28.3698 16.2275 24.5831 39.8312
41.8370 8.3265 37.2783 20.7091 32.4948 55.0544
2.1785 1.3564 1.8128 0.6693 2.2719 5.5885
六、实验任务
(1)对本实验中所得的C矩阵求C’ 、|C|、R(C)、、C 的特征值及对应的特征向量。

(2)创建从2开始,公差为4的等差数列的前15项构成的行向量;
(3)将本实验中矩阵A、B的对应元素相乘、对应元素相除并观察分母为0时的结果;
(4)求b每个元素自身次幂所得行向量;
(5)列出本实验中所有变量。