VIP专享控制系统数字仿真自考题型举例与解答

1、填空题基本概念题型4 35填空题基本概念题型4 35?= ?。阶跃题型举例与总复习一、A类1、系统是指相互联系又相互作用的 实体的有机组合。2、定义一个系统时，首先要确定系统的 边界；边界确定了系统的范围，边界以外对系统的作用称为系统的 输入，系统对边界以为环境的作用称为系统的 输出。3、系统的三大要素为： 实体、属性和活动。4、根据系统的 属性可以将系统分成两大类： 工程系统 和非工程系统 。5、相似原理用于仿真时，对仿真建模方法的三个基本要求是 稳定性 、准确性 和快速性 。6、根据模型种类不同，系统仿真可分为三种： 物理仿真 、数字仿真 和半实物仿真 。7、按照系统模型特征分类，仿真可分

2、为 连续系统 仿真及离散事件系统 仿真两大类。8、采用一定比例按照真实系统的样子制作的模型称为 物理模型 ，用数学表达式来描述系统内在规律的模型称为 数学模型 。9、计算机仿真是指将 模型在计算机上进行试验的过程。10、系统仿真的三个基本活动是 系统建模 、仿真建模 和仿真试验 ，计算机仿真的三个要素为：系统、模型与计算机。11、如果某数值计算方法的计算结果对初值误差和计算误差不敏感，则称该计算方法是 稳定的。12、数值积分法步长的选择应遵循的原则为计算 稳定性 及计算精度。13、采样数值积分方法时有两种计算误差，分别为 截断误差 和舍入误差 。14、三阶隐式啊达姆氏算法的截断误差为 O(?

3、)，二阶龙格 -库塔法的局部截断误差为 O(? )，四阶龙格 -库塔法的局部截断误差为 O(? )。15、在判定数值积分方法的稳定域时，使用的测试方程为16、龙格 -库塔法的基本思想是用几个点上函数值的 线性组合 来避免计算函数的高阶导数，提高数值计算的精度。17、连续系统仿真中常见的一对矛盾为计算 速度和计算 精度。18、离散相似法在采样周期的选择上应当满足 采样定理。19、保持器是一种将离散时间信号恢复成 连续信号 的装置，零阶保持器能较好地再现信号，一阶保持器能较好地再现 斜坡信号。20、实际信号重构器不可能无失真地重构信号，具体表现为信号重构器会对被重构的信号产生相位的 滞后和幅度的

4、衰减。21、一般将采样控制系统的仿真归类为 连续系统仿真。22、在控制理论中，由系统传递函数来建立系统状态方程的问题被称为 “实现问题 ”。23、常用的非线性环节包括 :饱和非线性、 失灵非线性、 迟滞回环 非线性。1简单计算题型D(?)0.02s，则该校正环节的数字仿真模型为：?(?)?(?)?0.5?1+ 0.06? = ?1?(?) ?= 0.5? ?0.06? + ? 0 ? 0.4。由题可知?(?) + ?(?) =?- + ?(?+ 1)?(?)= ?(?(?)?(?) ?+ 1p=-1；对于一阶系统，采用阶跃信号输入时，其稳态输G(z)的极点对应 ?(?)?：?=简单计算题型D(

5、?)0.02s，则该校正环节的数字仿真模型为：?(?)?(?)?0.5?1+ 0.06? = ?1?(?) ?= 0.5? ?0.06? + ? 0 ? 0.4。由题可知?(?) + ?(?) =?- + ?(?+ 1)?(?)= ?(?(?)?(?) ?+ 1p=-1；对于一阶系统，采用阶跃信号输入时，其稳态输G(z)的极点对应 ?(?)?：?= 1-?(?) -?(?)?(?)?(?)?- ?)?u?= =(?2?0.5?+ 0.06 1?0.5?1+ 0.06?21 2=?5?，y(0)?=?(?) ?(-1?)，系统 S域的传递函数为11?+ 1的极点，并将无穷远点作为?1?(1 )?

6、(1- )?-1?(?) ?(?)?= 0.5?1?0.06? + ?12?-0.5?1+ 0.06?2= ?11= 1-，(1- )?u?G(z)的零点，构建 G(z)：)?(1-?2?0.5?+ 0.062 1。，用欧拉法仿真，为保证计算稳定，则对计算步长5， 带入判定方程，可以计算出稳定域为0) = 1?）? ?Z传递函数为：?(?)(? = ?(?) -?，采样周期为)。h?(0,0.4)，用根匹配法求得的离散1? 1?1 ?) ?， (1 )?(?)(1 ?1?- = ?(?) -?)?)(124、已知某采样控制系统的数字校正环节为T=分析：由控制器的 Z 传递函数：D(?)= =

7、=(1 )?(?)经整理可得系统的差分数字模型为：?25、系统微分方程的要求为： ( 。分析：根据一阶显示方程 ?= ?的稳定性判定方程：|1+ ?| 1满足判定方程的解即为稳定域当步长的取值在 0到0.4 之间时，用欧拉方法仿真计算是稳定的。26、一个连续系统的微分方程为化模型为： （?=分析：先对原微分方程取拉氏变换，得?(?) = =由传递函数可知，系统无零点，有一个一级极点出可以由终值定理求得：?() = ?(?)?(?) = ? ?1?= 1? ?作根匹配替换，令G(z) = ? = ? ?再由 Z域的终值定理，系统在同样输入下稳态输出相同，求出?() = ?(1 ?1)?(?)?(

8、?)= ?(1? ? ? = ? =? ?由此求得 ，于是求得离散G(z) = =最终根据离散传递函数求输出序列：G(z) = =可得离散化差分模型为：?= +2?(?) =?=?= ?1?+?=?= ? ?(?)=1 ?(?+ 1) ?+ ?2 1 2(?(?) =?=?= ?1?+?=?= ? ?(?)=1 ?(?+ 1) ?+ ?2 1 2(?1) + ?(?+ 1) (?+ 2)?+ (?2)?(?)= = ?+ 0.25?3?(?)= =，= 2.25= ?0+ ?1= 2+ 0.252.25 +0.25)(2+ 0.25)?(?)?0 = 2 0.25简答题1?+ 1的近似脉冲传递

9、函数为?1?2? 1；?；2? ?+1代入 ?+1? + 1?(?)?(?)?(?)?(?)?1=（?(?)；1 1中，有：0.45?0.45?为小数部分。根据线性插值方法，取以下线性组合作为输出：=Z域模型，则由于步长?+ ?(?+ 2)?+ (?2)T=0.2，延迟时）。分析：常用的替换公式有1、 欧拉替换：2、 双线性替换：3、 根匹配替换：题目要求对连续系统作双线性替换，将连续模型转换为离散事件模型，可将双线性替换公式直接带入连续时间表达式求得。将?(?) = = =? ?+ 128、某纯延迟环节的输入为 u，输出为 y，传递函数 ，若取步长T=0.2，则这个环节的仿真模型为： （0.

10、75?2 ）分析：将 通过替换公式 ?= ?变换为离散时间间?= 0.45?= 0.45T 0.2?(?)= ? = ?(2?(?) = ?(?)?(?)= ?式中， 为整数部分，?= (1 ?0.25)?2+ 0.25?3= 0.75?2+ 0.25?3二、29、试简述为什么需要采用系统仿真方法对系统进行试验？答：这是因为1） 在系统建成前或设计阶段，系统没有建立起来，因此不可能在真实系统上进行试验；3模型运行块显示仿真试验仿真建模输出程序块打印计算机绘图2） 当在真实系统上做试验，会破坏系统的运行时，可采用仿真方法；模型运行块显示仿真试验仿真建模输出程序块打印计算机绘图3） 在真实系统上试

11、验，难以保证试验结果真实性和再现性，可采用仿真方法；4） 如果在真实系统上试验，会导致时间太长、费用太大或者有危险时，可选用仿真方法。30、简单介绍仿真程序的基本功能，并画出仿真程序结构图答：仿真程序包括以下基本功能模块1） 主程序 仿真运行控制，以便修改参数，选择算法等；2） 置初值 设置初始条件，设定系统参数；3） 模型运行 调用相应的仿真算法程序，完成仿真运算功能；4） 输出处理 输出仿真结构及仿真结果的处理。主程序（仿真逻辑控制）置初值（初始条件和参数）仿真算法子程序31、简述计算机仿真三要素及其三个基本活动，并用图表示答：系统仿真的三要素是系统、模型和计算机；对应的三个基本活动是系统

12、建模、仿真建模和仿真试验。其关系图示如下：系统系统建模模型32、试述系统仿真的一般步骤答：系统仿真一般步骤包括4以及输入矩阵 。?1= ?01 ?6?2= ?1?5?6+ ?01?以及输入矩阵 。?1= ?01 ?6?2= ?1?5?6+ ?01?3= ?2+ ?02?4= ?3?4?5= ?6 ?6= ?2+ ?4+ ?02?+ ?0?0W0010000001001000100000-1010-10000-1-10010 ?0=1100000010012） 仿真建模 根据系统及模型的特点选择合适的仿真算法；3） 程序设计 将仿真模型以合适的方式转换为可执行的计算机程序；4） 程序检验 检验仿

13、真算法的合理性和正确性；5） 仿真试验 运行仿真程序，并得到输出数据；6） 结果分析 根据仿真运行的结果，对系统进行分析，并形成报告。33、已知系统结构如下图所示，其中方框内的数字表示环节的编号，试写出系统的连接矩阵?解：由结构图，有：即?= ，其中：?=34、简述仿真建模方法的基本要求答：对仿真建模方法的基本要求为1） 稳定性 若原连续系统是稳定的，则离散化后得到的仿真模型也是稳定的；2） 准确性 计算结果的绝对误差或相对满足一定的误差要求；3） 快速性 每一步计算时间决定了仿真速度，仿真速度应满足实际仿真问题的需求。35、给出采样控制系统的典型结构，给出其仿真程序框图，并说明在采样系统的数

14、字仿真中，应如何处理计算步距和采样周期的关系？答：采样控制系统框图如下：5?+ 3?+ 2?= ?+ 2?+:?(?)?(?)?1?2 0?32?(?) =?+?+?，零初值，步长1?(1?+ 3?+ 2?= ?+ 2?+:?(?)?(?)?1?2 0?32?(?) =?+?+?，零初值，步长1?(1+ ?)= ?1+ ?11 + ?3?，已知系统为零初值，求系统的状态?2+ 2?+ 3?3+ 3?2+ 2?0=011?(1+h=0.1，求系统的状态方程和输出方程。1?110?2?)，若用面向结构图的数字仿真方法，典型环节取110 01 0?3 1?1?2?3?2 ?1?2?3+36、描述系统

15、的微分方程为空间表达式。解：首先将微分方程转换为传递函数。对原方程取拉氏变换，有(?3+ 3?2+ 2?)?(?)= (?2+ 2?+ 3)?(?)并由此求得系统的传递函数为?(?)= =由系统传递函数标准形式可以直接写出系统的能控标准形状态表达式：?(0)= 0输出方程：?= 337、已知系统传递函数? ?(?)=解：将传递函数分解为基本环节：?(?) = ?即两个基本环节级联，如下图所示，取状态变量如图标注：?6?1 = ?1?1? 11?=?+ ? ?(0) = 1 ?= 0.1? ?值，并说明造成差异的原因。? ?= ?(?,?) =?+ ?1=11?1 = ?1?1? 11?=?+

16、? ?(0) = 1 ?= 0.1? ?值，并说明造成差异的原因。? ?= ?(?,?) =?+ ?1=11?1 = ?(?,?)=? ?2? + ?2 ?(?+ ?1?,?+ =?(?+ ?1?) + (?+ )?3 ?+ ?2?,?+ =?(?+ ?2?) + (?+ )?4= ?(?+ ?3?,?+ ?)=?(?+ ?3?) + ?)?= ?0= 1?1 = ?20 + ?0=?12+ 0=?1 ?2= ?(?0+ ?)2+ ?0+ ?(1 ?1 0.1) + 0 + =?0.8525 ?3= ?(?0+ ?)2+ ?0+ ?(1 ?0.8525 0.1) + 0+ =?0.8666?4

17、 ?(?0+ ?3?)2+ + ?)=? (1 ?0.8666 0.1) + (0+ 0.1) =?0.7342 = ?=222(?(?6=2= 0.1 = 0 ?1= ?1 ? 1 0.12?1 2)= 2 21 ? 1 0.12?2 2)= 2 2=?2 ?2 ?2?2，试分别用欧拉法、四阶龙格201 ? 1 ?2 2) 2 2(?+?0( ( )( ( )(?00?1 0u1库塔法(21 ? 1 ?2 2) 2 2时?12220 12的值 ?1?2+ 又由题目条件知，系统时零状态，于是?(0)= 0系统输出方程为：?= 038、已知 ， ，取计算步距求?= 时的解：被求函数 的导函数 ，

18、以下分别用两种方法求解（1） 欧拉法由欧拉法的递推公式?+ = ?+ ?(?,?)?= ?+ ?2? + ?)?得：?1 ?0+ + ?0)?= 1+ ?12+ 0)0.1= 0.9（2） 四阶龙格 库塔法RK-4 的递推公式为：?+ = ?+ (?1+ 2?2+ 2?3+ ?4) ?其中由已知条件， ，? ，由 递推出7116RK-4 方法精度更高。?1?0(?)?1 ?1?1 1X = ?+?)= sin ?+ = ?(?)= sin ?+ = cos?1 22 1X?+ = 0?+(?6?1?2?(?)0?1=? ?0?，?(0)，用增广矩阵法将其转换为齐次方程，并标明初值。1 2=?1

19、210?10?116RK-4 方法精度更高。?1?0(?)?1 ?1?1 1X = ?+?)= sin ?+ = ?(?)= sin ?+ = cos?1 22 1X?+ = 0?+(?6?1?2?(?)0?1=? ?0?，?(0)，用增广矩阵法将其转换为齐次方程，并标明初值。1 2=?1210?10?1 ?+ 3= ?0，? 00 ?00 1?3 ?2 0?3，?= 0 ?= 10113113(?，?为?维状态变量，若输，? 0 ?+ 0(? 1)0 ?3(?)?11，1?3?1)?+2001?(?+ 3)?3(?)0?= ?3?)?1+1?+ 33? 19(1 31 ?= 1+ ( ?1?

20、20.8525 ?20.8666 ?0.7342)0.1 = 0.9138（3）计算结果产生差异是由于两种方法的精度不一样，39、（时域离散相似法）已知线性定常系统的状态方程为： ?= 1求 ， ，及离散化状态方程。解：根据题意，有：?=(?) = =系统的状态转移矩阵：(?)= ?= ?1(?)系统的输入转移矩阵：?(?) = (?系统的离散化状态方程及离散化输出方程为：?+ = (?)?+ ?(?)?= ?40、已知线性时不变系统的状态方程为入?(解：设则有：?+ = ?cos?= ?+?+ ?sin?=?+将其增广到系统状态方程上，有： 8X(0)?+?+?(?) = =?(?) T?=

21、 ?1?1 ?2(?+ 1)22 ?1 2 ?1 4(?X(0)?+?+?(?) = =?(?) T?= ?1?1 ?2(?+ 1)22 ?1 2 ?1 4(?1) + 10?(?1)(?+ 1) + 6(?+ 1)?+ 1) ?+ 1?(?)= (1 )(1?1)(? 3(1? ?3?)1?2?2?)根匹配法?(?) 1 1?(?) (?+ 2)(?+ 3)?1=?2(?K?；对于连续模型，当系统输入为阶跃信号时，应用终值定理1 1? 6?1(0)2(0)?(?) 1?(?)2? +1，代入系统传递函数作替换2 221? 1 1? ?(?+ 2)(?+ 3)12? 1?3(?3?)?2+ 5

22、?+ 6?2，?2?2?)(?3?)+ 5?+ 60=1?2+5?+ 6，分别用下述方法求取与之近似等效?+ 5?+ 61?2? 1?) )(?1?=?3023?3?，无有限零点；根据根匹配法，有系统离散传递函数：(1 41、二阶连续系统的传递函数为的脉冲传递函数 ，计算步长取 。解：（1） 双线性替换法根据双线性替换公式?(?)= =( + 5 + 6（2） 频域离散相似法（加虚拟采样开关及零阶保持器）串入零阶保持器 并求系统的 变换?(?)= ?(?)?(?) = ?1 ? = ?1)? = ? ?2=2(?（3）?(?) = = =系统有两个一阶极点?(?)= K?现根据终值相等，确定增

23、益?() = lim ?(?)?(?)= lim ?(?)1= lim =? ? ?2对于离散模型，同样阶跃输入时，应有相同的稳态输出，应用终值定理916(11 ?2(1? )(1?2? 3?(?)= 2.620.1秒下进行数字仿真，求此时?1?=1?=?= 2.62?0.98?()= ?(?：(12? 3?(1? )?(? ?)(? )?0.98?16(11 ?2(1? )(1?2? 3?(?)= 2.620.1秒下进行数字仿真，求此时?1?=1?=?= 2.62?0.98?()= ?(?：(12? 3?(1? )?(? ?)(? )?0.98?0.64，采样周期为?(?，由于?(?)?，于是?= ? = ?0.98?0.645?0.64 0.1074) (1?1(1 ) (?)(1? )? ?2