MATLAB矩阵实验报告

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

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

disp('E

disp(E

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))

输出结果为：