MATLAB矩阵实验报告

MATLAB程序设计实验

班级：电信1104班 姓名：龙刚 学号：1404110427

实验内容：了解MATLAB基本使用方法和矩阵的操作

一．实验目的

1.了解MATLAB的基本使用方法。

2.掌握MATLAB数据对象的特点和运算规则。

3.掌握MATLAB中建立矩阵的方法和矩阵的处理方法。

二．实验内容

浏览MATLAB的start菜单，了解所安装的模块和功能。

建立自己的工作目录，使用MATLAB将其设置为当前工作目录。使用path命令和工作区浏览两种方法。

使用Help帮助功能，查询inv、plot、max、round等函数的用法和功能。使用help命令和help菜单。

建立一组变量，如x=0:pi/10:2*pi，y=sin(x)，在命令窗口显示这些变量；在变量窗口打开这些变量，观察其值并使用绘图菜单绘制y。

分多行输入一个MATLAB命令。

求表达式的值

，a=3.5，b=5，c=-9.8

，

7.已知

，

求

A+6B，A2-B+I

A*B，A.*B，B*A

A/B，B/A

[A,B]，[A([1,3], :); B^2]

8.已知

输出A在[10,25]范围内的全部元素

取出A的前三行构成矩阵B，前两列构成矩阵C，右下角3x2子矩阵构成矩阵D，B与C的乘积构成矩阵E

分别求表达式E

9.已知

求A的特征值和特征向量，分析其数学意义。

代码与实现

代码：

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

运行结果：

代码：

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

运行结果为：

代码：

disp('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))

输出结果为：

代码为：

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)

输出结果为：