《MATLAB及应用》参考答案

《MATLAB及应用》上机作业

学院名称： 机械工程学院 专业班级： 测控1201 学生姓名： 学生学号：

201 年 4 月

1

《MATLAB及应用》上机作业要求及规范

一、作业提交方式：word文档打印后提交。 二、作业要求：

1．封面：按要求填写学院、班级、姓名、学号，不要改变封面原有字体及大小。 2．内容：只需解答过程（结果为图形输出的可加上图形输出结果），不需原题目；为便于批阅，题与题之间应空出一行；每题答案只需直接将调试正确后的M文件内容复制到word中（不要更改字体及大小），如下所示：

%作业1_1 clc

A=[1 2 3 4;2 3 5 7;1 3 5 7;3 2 3 9;1 8 9 4];

B=[1+4*i 4 3 6 7 7;2 3 3 5 5 4+2*i;2 6+7*i 5 3 4 2;1 8 9 5 4 3]; C=A*B D=C(4:5,4:6)

三、大作业评分标准：

1．提交的打印文档是否符合要求； 2．作业题的解答过程是否完整和正确；

3．答辩过程中阐述是否清楚，问题是否回答正确；

4．作业应独立完成，严禁直接拷贝别人的电子文档，发现雷同者都以无成绩论处。

2

作业1

1、用MATLAB可以识别的格式输入下面两个矩阵

?1??2A??1??3?1?234?4?1?4i?357??23357?，B???26?7i?239???18??894?7??3554?2i? ?5342?9543??367再求出它们的乘积矩阵C，并将C矩阵的右下角2?3子矩阵赋给D矩阵。赋值完成后，调用相应的命令查看MATLAB工作空间的占有情况。 解：

A=[1 2 3 4;2 3 5 7;1 3 5 7 ;3 2 3 9 ;1 8 9 4;]

B=[1+4i 4 3 6 7 7;2 3 3 5 5 4+2i;2 6+7i 5 3 4 2;1 8 9 5 4 3;] B=[1+4i 4 3 6 7 7;2 3 3 5 5 4+2i;2 6+7i 5 3 4 2;1 8 9 5 4 3;] C=A*B D=C(4:5,4:6);

?162313???511108?，求A，A?1，A3，2A?A?1，3A?1?A'，并求矩阵A的2、设矩阵??97612???414152??特征值和特征向量。 解：

A=[16 2 3 13;5 11 10 8;9 7 6 12;4 14 15 2;] det(A) inv(A)

3

A.^3 2*A+inv(A) 3*A-A' [V,D]=eig(A) abs (A)

3、解下列矩阵方程:

?010??100??1?43??????? 100X001?20?1

???????001??010??1?20???????解：A=[0 1 0;1 0 0;0 0 1;]; B=[1 0 0;0 0 1;0 1 0;]; C=[1 -4 3; 2 0 -1;1 -2 0;]; X=inv(A)*C*inv(B)

4、求多项式f(x)?(x?5x?35x?13)(x?45x?3)(31x?3)当x?3时的值，并求

322f(x)的导数。

解：a1=[1 5 35 13]

a2=[0 1 45 3] a3=[0 0 31 3] p1=conv(a1,a2) p1=conv(p1,a3) polyval(p1,3) polyder(p1)

4

5、求多项式f?x?5x?14x?10x?3的根和导数。 解:

432p1=[1 -5 -14 -10 -3] roots(p1) polyder(p1)

10(s4?4s2?5s2?6s?7)6、对于有理多项式 f(s)? (1)计算该多项式相除的结果;(2)

(s?1)(s?2)(s?3)将该多项式展开为部分分式的形式;(3)计算

df。 ds解(1)A=[1 0 4 5 6 7]; B=A.*10; C=[1 1]; D=[1 2]; E=[1 3];

F=conv(conv(C,D),E); G=deconv(B,F); (2)[r,p,k]=residue(B,F) (3)[p,q]=ployder(B,F)

7、在某次传感器输入输出特性实验中测得输入输出的一组数据如下表所示:

x（输入） 1 1。3 2 1。8 3 2。2 4 2。9 5 3。5 y（输出） 5

本文来源：https://www.bwwdw.com/article/pe4v.html