文档介绍:北京工业大学
Matlab实验报告
姓名:
学号:
指导教师:吕策
实验二、Mat lab的基本计算
(一) 实验目的
掌握建立矩阵的方法。
掌握Matlab各种表达式的书写规则以及常用函数的使用。
能用Mat
(1) KU = A+6*B和K12 = A_3 + / (其中[为单位矩阵)
» k 12=a-b+eye (3, 3)
kl2 =
12
31
-3
32
8
84
0
0
1
» kll=a+6*b
kll =
18
52
-10
46
7
105
21
455
49
(2) K21 = A^B和K22 = A*8
>> k21=a*b
k21 =
68
-224
62
309
5757
596
154
464
241
» k22=a. *b
k22 =
12
102
4
68
0
261
9
4225
49
(3)K31 = A"3 和 K32 = A"3
» k31=aA3
k31 =
37226
233824
48604
247370
149188
600766
78688
454142
118820
» k32=a. “3
k32 =
1728
39304
-64
39304
343
658503
27
274625
343
⑷ K41 = A/8和K42 = 8\A
>> k41=a/b
k41 =
0. 1118
-
0
1. 0000
» k42=b\a
k42 =
19. 5265
-0. 7794
-
-
⑸ K51 = [A B]和 K52 = [A([l,3],:);矿 2]
k51 =
12
34
-4
1
3
-1
34
7
87
2
0
3
3
65
7
3
65
7
» k52=[a([l,3], :) ;b*2]
k52 =
12
34
-4
3
65
7
4
-62
1
11
201
19
154
464
241
:
1/2
1/3
1/4'
xl
「「
1/3
1/4
1/5
x2
=
_1/4
1/5
1/6
x3
(1)求方程的解(矩阵除法和求逆法)
» z=[0. 95 0. 67 ]
0. 9500
» x=[l/2 1/3
1/4;1/3
1/4 1/5;1/4 1/5 1/6]
x =
0. 3333
0. 2500
0. 2000
0. 2000
0. 1667
» y=z/x
y =
1. 2000
0. 6000
0. 6000
yi =
⑵,再求解,并比较饥的变化和解的相对变化。
» z=[0. 95 0. 67 ]
0. 9500 0. 5300
» x=[l/2 1/3 1/4;1/3
1/4 1/5;1/4 1/5 1/6]
0. 2000
» y=z/x
3. 0000
-6. 6000
6. 6000
» yl=z*inv(x)
yi =
3. 0000 -6. 6000 6. 6000
,方差为1的6X6正态分布随机矩阵C,然后统计C中大于 -,。
0. 0493
0. 6343
0. 7546
0.