《自动控制系统计算机仿真》习题参考答案

《自动控制系统计算机仿真》习题参考答案

1-1 什么是仿真？ 它的主要优点是什么？它所遵循的基本原则是什么？

答：所谓仿真，就是使用其它相似的系统来模仿真实的需要研究的系统。计算机仿真是指以数字计算机为主要工具，编写并且运行反映真实系统运行状况的程序。对计算机输出的信息进行分析和研究，从而对实际系统运行状态和演化规律进行综合评估与预测。它是非常重要的设计自动控制系统或者评价系统性能和功能的一种技术手段。

仿真的主要优点是：方便快捷、成本低廉、工作效率和计算精度都很高。它所遵循的基本原则是相似性原理。

1-2 你认为计算机仿真的发展方向是什么？

答：向模型更加准确的方向发展，向虚拟现实技术，以及高技术智能化、一体化方向发展。向更加广阔的时空发展。

1-3 计算机数字仿真包括哪些要素？它们的关系如何？

答：计算机仿真的三要素是：系统——研究的对象、模型——系统的抽象、计算机——仿真的工具和手段。它们的关系是相互依存。

2-1 控制算法的步长应该如何选择？

答：控制算法步长的选择应该恰当。如果步长太小，就会增加迭代次数，增加计算量；如果步长太大，计算误差将显著增加，甚至造成计算结果失真。 2-2 通常控制系统的建模有哪几种方法？

答：1）机理建模法；2）实验建模法；3）综合建模法。

2-3 用欧拉法求以下系统的输出响应()y t 在0≤t ≤1上，0.1h =时的数值解。

0y

y += ， (0)0.8y = 解：输入以下语句 绘制的曲线图

2-4 用二阶龙格-库塔法对2-3题求数值解，并且比较两种方法的结果。解：输入以下语句绘制的曲线图

经过比较两种方法的结果，发现它们几乎没有什么差别。

3-1 编写两个m文件，分别使用for和while循环语句计算200

3

1

k

k

=∑。

解：第1个m文件，第2个m文件运行结果都是

3-2 求解以下线性代数方程：

1

2

3

1022

1131

3121

x

x

x

⎡⎤⎡⎤⎡⎤

⎢⎥⎢⎥⎢⎥

=

⎢⎥⎢⎥⎢⎥

⎢⎥⎢⎥⎢⎥

⎣⎦⎣⎦⎣⎦

解：输入语句计算结果

3-3 已知矩阵

013

=121

542

⎡⎤

⎢⎥

⎢⎥

⎢⎥

⎣⎦

A，

218

=414

332

⎡⎤

⎢⎥

⎢⎥

⎢⎥

⎣⎦

B

试分别求出A阵和B阵的秩、转置、行列式、逆矩阵以及特征值。

解：求矩阵的秩、行列式、逆矩阵、特征值的函数分别为：rank( ), det( ), inv( ), eig( )。求矩阵转置的命令为“.’ ”。

求出A 阵的秩为3、转置为015=124312⎡⎤

⎢⎥'⎢⎥⎢⎥⎣⎦A 、行列式为15-、逆矩阵10-0.66670.3333=-0.21-0.20.4-0.33330.0667-⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦A 以及特征值 6.1926()30.8074eig A ⎡⎤

⎢⎥=-⎢⎥

⎢⎥⎣⎦

；

求出B 阵的秩为3、转置为243B =113842⎡⎤

⎢⎥'⎢⎥⎢⎥⎣⎦、行列式为56、逆矩阵1-0.17860.3929-0.0714=0.0714-0.35710.42860.1607-0.0536-0.0357B -⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦以及特征值9.1606() 2.0803 1.33632.0803 1.3363eig B i i ⎡⎤

⎢⎥=-+⎢⎥

⎢⎥--⎣⎦

。

3-4 对于3-3题中的A 阵和B 阵，在Command 窗口中分别求出2

C =A 、

D 矩阵为A 中每个元素平方组成的矩阵、

E 矩阵为A 阵乘以B 阵、

F 矩阵为A 阵和B 阵数组乘积（即：对应元素分别相乘的积构成的矩阵）。 解：输入命令

，

，

，

，计算结果为

，

，

，

3-5 已知某系统的闭环传递函数()s Φ如下，试用roots （）命令来判断系统的稳定性。

25432

325

()24576

s s s s s s s s ++Φ=+++++ 解：输入命令

计算结果如下，系统不稳定

3-6 求复数矩阵135********i i i i i i +-+⎡⎤

=⎢⎥++-⎣⎦

C 的转置C1与共轭转置C2。

解：输入语句

计算结果

4-1 某系统的传递函数为

2321.323

()0.5 1.21

s s G s s s s ++=+++

使用MATLAB 求出状态空间表达式和零极点模型。 解：输入语句

计算机返回

表明该系统的状态空间表达式为

0.5 1.211X 100X 00100u ---⎡⎤⎡⎤

⎢⎥⎢⎥=+⎢⎥⎢⎥

⎢⎥⎢⎥⎣⎦⎣⎦

， []1.32

3X

y = 该系统的零极点模型为 1.3(0.7692 1.31)(0.7692 1.31)

()(0.1153 1.1642)(0.1153 1.1642)(0.7307)

s i s i G s s i s i s +-++=

---++