MATLAB程序设计实验
班级：电信1104班姓名：龙刚学号：1404110427
实验内容：了解MA TLAB基本使用方法和矩阵的操作
一．实验目的
1.了解MA TLAB的基本使用方法。
2.掌握MA TLAB数据对象的特点和运算规则。
3.掌握MA TLAB中建立矩阵的方法和矩阵的处理方法。
二．实验内容
1.浏览MATLAB的start菜单，了解所安装的模块和功能。
2.建立自己的工作目录，使用MA TLAB将其设置为当前工作目录。
使用path命令和工作
区浏览两种方法。
3.使用Help帮助功能，查询inv、plot、max、round等函数的用法和功能。
使用help命
令和help菜单。
4.建立一组变量，如x=0:pi/10:2*pi，y=sin(x)，在命令窗口显示这些变量；在变量窗口打
开这些变量，观察其值并使用绘图菜单绘制y。
5.分多行输入一个MA TLAB命令。
6.求表达式的值
)6
10.3424510
w-
=+⨯
()
2
2
tan
b c
a e
abc
x
b c a
π
π
+
+-
+
=
++
，a=3.5，b=5，
c=-9.8
(20.5ln t
z e t
=，
213
50.65
i
t
-
⎡⎤
=⎢⎥
-
⎣⎦
7.已知
154
078
3617
A
--
⎡⎤
⎢⎥
=⎢⎥
⎢⎥
⎣⎦
，
831
253
320
B
-
⎡⎤
⎢⎥
=⎢⎥
⎢⎥
-⎣⎦求
A+6B，A2-B+I
A*B，A.*B，B*A
A/B，B/A
[A,B]，[A([1,3], :); B^2]
8.已知
23100.7780414565532503269.54
54 3.14A -⎡⎤⎢⎥-⎢⎥=⎢⎥⎢⎥-⎣⎦ 输出A 在[10,25]范围内的全部元素
取出A 的前三行构成矩阵B ，前两列构成矩阵C ，右下角3x2子矩阵构成矩阵D ，B 与C 的乘积构成矩阵E
分别求表达式E<D ，E&D ，E|D ，(~E) | (~D)
9.已知
2961820512885A -⎡⎤⎢⎥=⎢⎥⎢⎥-⎣⎦
求A 的特征值和特征向量，分析其数学意义。
三、代码与实现
6.代码：
disp('1.求表达式的值')
w=sqrt(2)*(1+0.34245*10^(-6))
a=3.5;b=5;c=-9.8;
x=(2*pi*a+(b+c)/(tan(b+c)+a))
t=[2 1-3i;5 -0.65];
z=0.5*exp(2*t)*log(t+sqrt(1+t^2))
运行结果：
7.代码：
disp('2.已知A=[-1 5 -4;0 7 8;3 61 7];B=[8 3 -1;2 5 3;-3 2 0];求下列值')
A=[-1 5 -4;0 7 8;3 61 7];
B=[8 3 -1;2 5 3;-3 2 0];
I=[1 0 0;0 1 0;0 0 1];
disp('A+6*B=')
disp(A+6*B)
disp('A.^2-B+I=')
disp(A.^2-B+I)
disp('A*B=')
disp(A*B)
disp('A.*B=')
disp(A.*B)
disp('B*A=')
disp(B*A)
disp('A/B=')
disp(A/B)
disp('A\B=')
disp(A\B)
disp('[A,B]=')
disp([A,B])
disp('[A([1,3],:);B^2]=')
disp([A([1,3],:);B^2])
运行结果为：
8.代码：
d isp('3.已知A=[23 10 -0.778 0;41 -45 65 5;32 5 0 32;6 -9.54 54 3.14];') disp('(1)输出A在[10,25]范围内的全部元素;')
A=[23 10 -0.778 0;41 -45 65 5;32 5 0 32;6 -9.54 54 3.14];
m=find(A>=10&A<=25);
disp(A(m))
disp('(2)取出A的前3行构成矩阵B，前两列构成矩阵C，右下角3X2子矩阵构成矩阵D,B 与C的乘积构成矩阵E;')
B=A([1:3],:)
C=A(:,[1,2])
D=A([2:4],[3,4])
E=B*C
disp('(3)分别求表达式E<D，E&D，E|D，(~E) | (~D)。
')
disp('E<D=')
disp(E<D)
disp('E&D=')
disp(E&D)
disp('E|D=')
disp(E|D)
disp('~E|~D=')
disp((~E)|(~D))
disp('(4)生成一个Hilbert矩阵和Parscal矩阵，求其行列式的值')
H=hilb(5)
P=pascal(5)
disp('det(H)=')
disp(det(H))
disp('det(P)=')
disp(det(P))
输出结果为：
9.代码为：
disp('4.已知A=[-29 6 18;20 5 12;-8 8 5];求A的特征值和特征向量，并分析其数学意义')
A=[-29 6 18;20 5 12;-8 8 5];
disp('特征值为：')
disp(eig(A))
disp('特征向量为：')
[V,D]=eig(A);
disp('A的特征向量分别为：')
a1=V(:,1)
a2=V(:,2)
a3=V(:,3)
输出结果为：。