控制系统数字仿真题型举例与总复习一、填空题A类基本概念题型1、系统是指相互联系又相互作用的实体的有机组合。
2、定义一个系统时，首先要确定系统的边界；边界确定了系统的范围，边界以外对系统的作用称为系统的输入，系统对边界以为环境的作用称为系统的输出。
3、系统的三大要素为：实体、属性和活动。
4、根据系统的属性可以将系统分成两大类：工程系统和非工程系统。
5、相似原理用于仿真时，对仿真建模方法的三个基本要求是稳定性、准确性和快速性。
6、根据模型种类不同，系统仿真可分为三种：物理仿真、数字仿真和半实物仿真。
7、按照系统模型特征分类，仿真可分为连续系统仿真及离散事件系统仿真两大类。
8、采用一定比例按照真实系统的样子制作的模型称为物理模型，用数学表达式来描述系统内在规律的模型称为数学模型。
9、计算机仿真是指将模型在计算机上进行试验的过程。
10、系统仿真的三个基本活动是系统建模、仿真建模和仿真试验，计算机仿真的三个要素为：系统、模型与计算机。
11、如果某数值计算方法的计算结果对初值误差和计算误差不敏感，则称该计算方法是稳定的。
12、数值积分法步长的选择应遵循的原则为计算稳定性及计算精度。
13、采样数值积分方法时有两种计算误差，分别为截断误差和舍入误差。
14、三阶隐式啊达姆氏算法的截断误差为O(ℎ4)，二阶龙格-库塔法的局部截断误差为O(ℎ3)，四阶龙格-库塔法的局部截断误差为O(ℎ5)。
15、在判定数值积分方法的稳定域时，使用的测试方程为ẏ=μy。
16、龙格-库塔法的基本思想是用几个点上函数值的线性组合来避免计算函数的高阶导数，提高数值计算的精度。
17、连续系统仿真中常见的一对矛盾为计算速度和计算精度。
18、离散相似法在采样周期的选择上应当满足采样定理。
19、保持器是一种将离散时间信号恢复成连续信号的装置，零阶保持器能较好地再现阶跃信号，一阶保持器能较好地再现斜坡信号。
20、实际信号重构器不可能无失真地重构信号，具体表现为信号重构器会对被重构的信号产生相位的滞后和幅度的衰减。
21、一般将采样控制系统的仿真归类为连续系统仿真。
22、在控制理论中，由系统传递函数来建立系统状态方程的问题被称为“实现问题”。
23、常用的非线性环节包括:饱和非线性、失灵非线性、迟滞回环非线性。
B类简单计算题型24、已知某采样控制系统的数字校正环节为D(z)=Y(z)U(z)=zz2−0.5z+0.06，采样周期为T=0.02s，则该校正环节的数字仿真模型为：( y k=0.5y k−1−0.06y k−2+u k−1)。
分析：由控制器的 Z传递函数：D(z)=Y(z)U(z)=zz2−0.5z+0.06=z−11−0.5z−1+0.06z−2(1−0.5z−1+0.06z−2)Y(z)=z−1U(z)y k−0.5y k−1+0.06y k−2=u k−1经整理可得系统的差分数字模型为：y k=0.5y k−1−0.06y k−2+u k−1。
25、系统微分方程dydt=−5y，y(0)=1，用欧拉法仿真，为保证计算稳定，则对计算步长h 的要求为：( 0<h<0.4 )。
分析：根据一阶显示方程ẏ=μy的稳定性判定方程：|1+μℎ|<1满足判定方程的解即为稳定域。
由题可知μ=−5，带入判定方程，可以计算出稳定域为h∈(0,0.4)，当步长的取值在0到0.4之间时，用欧拉方法仿真计算是稳定的。
26、一个连续系统的微分方程为y′(t)+y(t)=u(t)，y(0)=1，用根匹配法求得的离散化模型为：（y k=e−T∙y k−1+(1−e−T)∙u k）分析：先对原微分方程取拉氏变换，得(s+1)Y(s)=U(s)，系统S域的传递函数为G(s)=Y(s)U(s)=1s+1由传递函数可知，系统无零点，有一个一级极点p=-1；对于一阶系统，采用阶跃信号输入时，其稳态输出可以由终值定理求得：y(∞)=lims→0[sG(s)U(s)]=lims→0[s∙1s+1∙1s]=1作根匹配替换，令G(z)的极点对应G(s)的极点，并将无穷远点作为G(z)的零点，构建G(z)：G(z)=K z∙zz−e pT=K z∙zz−e−T再由Z域的终值定理，系统在同样输入下稳态输出相同，求出K z：y(∞)=limz→1[(1−z−1)G(s)U(s)]=limz→1[(1−z−1)∙K z∙zz−e−T∙11−z−1]=limz→1[K z∙zz−e−T]=1由此求得K z=1−e−T，于是求得离散Z传递函数为：G(z)=Y(z)U(z)=(1−e−T)zz−e−T最终根据离散传递函数求输出序列：G(z)=Y(z)U(z)=(1−e−T)zz−e−TY(z)(z−e−T)=U(z)(1−e−T)z Y(z)(1−z−1∙e−T)=U(z)(1−e−T)可得离散化差分模型为：y k =e −T ∙y k−1+(1−e −T )∙u k27、用双线性替换法求得的系统G (s )=1s+1的近似脉冲传递函数为 （G (z )=Tz+T(T+2)z+(T−2)）。
分析：常用的替换公式有 1、 欧拉替换：s =z−1T ；2、 双线性替换： s =2T∙z−1z+1；3、 根匹配替换： z =e sT ；题目要求对连续系统作双线性替换，将连续模型转换为离散事件模型，可将双线性替换公式直接带入连续时间表达式求得。
将 s =2T ∙z−1z+1代入G (s )=1s+1中，有：G (z )=12T ∙z −1z +1+1=T (z +1)2(z −1)+T(z +1)=Tz +T(T +2)z +(T −2)28、某纯延迟环节的输入为u ，输出为y ，传递函数G (s )=Y (s )U (s )=e −0.45s ，若取步长T=0.2，则这个环节的仿真模型为：（0.75u k−2+0.25u k−3）分析：将G (s )=Y (s )U (s )=e −0.45s 通过替换公式z =e sT 变换为离散时间Z 域模型，则由于步长T=0.2，延迟时间τ=0.45，τT =0.450.2=2.25=C 0+C 1=2+0.25 G (z )=z −2.25=z −(2+0.25) Y (z )=G (z )U (z )=z −(2+0.25)U(z)式中，C 0=2为整数部分，C 1=0.25为小数部分。
根据线性插值方法，取以下线性组合作为输出：y k =(1−0.25)u k−2+0.25u k−3=0.75u k−2+0.25u k−3二、简答题29、试简述为什么需要采用系统仿真方法对系统进行试验？答：这是因为1）在系统建成前或设计阶段，系统没有建立起来，因此不可能在真实系统上进行试验；2）当在真实系统上做试验，会破坏系统的运行时，可采用仿真方法；3）在真实系统上试验，难以保证试验结果真实性和再现性，可采用仿真方法；4）如果在真实系统上试验，会导致时间太长、费用太大或者有危险时，可选用仿真方法。
30、简单介绍仿真程序的基本功能，并画出仿真程序结构图答：仿真程序包括以下基本功能模块1）主程序——仿真运行控制，以便修改参数，选择算法等；2）置初值——设置初始条件，设定系统参数；3）模型运行——调用相应的仿真算法程序，完成仿真运算功能；4）输出处理——输出仿真结构及仿真结果的处理。
31、简述计算机仿真三要素及其三个基本活动，并用图表示答：系统仿真的三要素是系统、模型和计算机；对应的三个基本活动是系统建模、仿真建模和仿真试验。
其关系图示如下：32、试述系统仿真的一般步骤答：系统仿真一般步骤包括1）系统建模——根据研究目的，建立实际系统的模型；2）仿真建模——根据系统及模型的特点选择合适的仿真算法；3）程序设计——将仿真模型以合适的方式转换为可执行的计算机程序；4）程序检验——检验仿真算法的合理性和正确性；5）仿真试验——运行仿真程序，并得到输出数据；6）结果分析——根据仿真运行的结果，对系统进行分析，并形成报告。
33、已知系统结构如下图所示，其中方框内的数字表示环节的编号，试写出系统的连接矩阵W以及输入矩阵W0。
解：由结构图，有：{u1=y01−y6u2=y1−y5−y6+y01u3=y2+y02u4=y3−y4u5=y6u6=y2+y4+y02即U=WY+W0Y0，其中：W=[ 0 0 01 0 00 1 00 0−10−1−10 0 00 0 10 0 00 1 0−1 0 000 110 0 ]W0=[101001000001]34、简述仿真建模方法的基本要求答：对仿真建模方法的基本要求为1）稳定性——若原连续系统是稳定的，则离散化后得到的仿真模型也是稳定的；2）准确性——计算结果的绝对误差或相对满足一定的误差要求；3）快速性——每一步计算时间决定了仿真速度，仿真速度应满足实际仿真问题的需求。
35、给出采样控制系统的典型结构，给出其仿真程序框图，并说明在采样系统的数字仿真中，应如何处理计算步距和采样周期的关系？答：采样控制系统框图如下：三、计算题36、描述系统的微分方程为y⃛+3ÿ+2ẏ=ü+2u̇+3u，已知系统为零初值，求系统的状态空间表达式。
解：首先将微分方程转换为传递函数。
对原方程取拉氏变换，有(s3+3s2+2s)Y(s)=(s2+2s+3)U(s)并由此求得系统的传递函数为:G(s)=Y(s)()=s2+2s+332由系统传递函数标准形式可以直接写出系统的能控标准形状态表达式：[x1x2 x3]=[0100010−2−3][x1x2x3]+[1]X(0)=0输出方程：y=[321][x1x2x3]37、已知系统传递函数G(s)=1s(1+s)，若用面向结构图的数字仿真方法，典型环节取G i(s)=C i+D i sA i+B i s，零初值，步长h=0.1，求系统的状态方程和输出方程。
解：将传递函数分解为基本环节：G (s )=1s(1+s)=1s ∙11+s即两个基本环节级联，如下图所示，取状态变量如图标注：ux 1x 2y由图可知，x 1=u ，x 2=x 1−x 2，y =x 2，写成矩阵形式，有系统状态方程：[x 1x 2]=[001−1][x 1x 2]+[10]u 又由题目条件知，系统时零状态，于是X (0)=0系统输出方程为：y =[01][x 1x 2]38、已知ẏ=−y 2+t ，y (0)=1，取计算步距ℎ=0.1，试分别用欧拉法、四阶龙格—库塔法求t =ℎ时的y 值，并说明造成差异的原因。
解：被求函数y 的导函数ẏ=f (y,t )=−y 2+t ，以下分别用两种方法求解（1） 欧拉法由欧拉法的递推公式y n+1=y n +f (y n ,t n )∙ℎ=y n +(−y n 2+t n )∙ℎ得：y 1=y 0+(−y 02+t 0)∙ℎ=1+(−12+0)×0.1=0.9（2） 四阶龙格—库塔法 RK-4的递推公式为：y n+1=y n +16(K 1+2K 2+2K 3+K 4)∙ℎ其中{K 1=f (y n ,t n )=−y n 2+t n K 2=f (y n +12K 1∙ℎ,t n +ℎ2)=−(y n +12K 1∙ℎ)2+(t n +ℎ2)K 3=f (y n +12K 2∙ℎ,t n +ℎ2)=−(y n +12K 2∙ℎ)2+(t n +ℎ2)K 4=f (y n +K 3∙ℎ,t n +ℎ)=−(y n +K 3∙ℎ)2+(t n +ℎ)1s11+s由已知条件，y n=y0=1，ℎ=0.1，由t0=0递推出t1=ℎ时y1的值{K1=−y02+t0=−12+0=−1 K2=−(y0+12K1∙ℎ)2+(t0+ℎ2)=−(1−1×12×0.1)2+(0+0.12)=−0.8525K3=−(y0+12K2∙ℎ)2+(t0+ℎ2)=−(1−0.8525×12×0.1)2+(0+0.12)=−0.8666K4=−(y0+K3∙ℎ)2+(t0+ℎ)=−(1−0.8666×0.1)2+(0+0.1)=−0.7342y1=y0+16(K1+2K2+2K3+K4)∙ℎ=1+16(−1−2×0.8525−2×0.8666−0.7342)×0.1=0.9138（3）计算结果产生差异是由于两种方法的精度不一样，RK-4方法精度更高。
