控制系统数字仿真 要点

1、中国地质大学（北京）继续教育学院2017年03课程考试第 页（共5页）第 页（共5页）词汇表解析法：就是运用已经掌握的理论知识对控制系统进行理论上的分析、计算。它是一种纯理论上的试验分析方法，在对系统的认识过程中具有普遍意义。实验法：对于已经建立的实际系统，利用各种仪器仪表及装置，对系统施加一定类型的信号，通过测取系统的响应来确定系统性能的方法。仿真分析法：就是在模型的基础上所进行的系统性能分析与研究的实验方法，它所遵循的基本原则是相似原理。模拟仿真：采用数学模型在计算机上进行的试验研究称之为模拟仿真。数字仿真：采用数学模型，在数字计算机上借助于数值计算的方法所进行的仿真试验称之为数字仿真。混

2、合仿真：将模拟仿真和数字仿真结合起来的仿真方法。数值计算：有效使用数字计算机求数学问题近似解的方法与过程。数值计算主要研究如何利用计算机更好的解决各种数学问题，包括连续系统离散化和离散形方程的求解，并考虑误差、收敛性和稳定性等问题。病态问题：闭环极点差异非常大的控制系统叫做病态系统，解决这类系统的问题就叫病态问题。显式算法：在多步法中，若计算第k+1次的值时，需要的各项数据均是已知的，那么这种算法就叫做显式算法。隐式算法：在多步法中，若计算第k+1次的值时，又需要用到第k+1次的值，即算式本身隐含着当前正要计算的量，那么这种算法就叫做隐式算法。数值稳定性：数值积分法求解微分方程，实质上是通过差

3、分方程作为递推公式进行的。在将微分方程离散为差分方程的过程中，有可能将原本稳定的系统变为不稳定系统。如果某个数值计算方法的累积误差不随着计算时间无限增大，则这种数值方法是稳定的，反之是不稳定的。实体：就是存在于系统中的具有实际意义的物体。属性：就是实体所具有的任何有效特征。活动：系统内部发生的任何变化过程称之为内部活动；系统外部发生的对系统产生影响的任何变化过程称之为外部活动。描述模型：是一种抽象的、无实体的，不能或者很难用数学方法精确表示的，只能用语言描述的系统模型。16.CAD技术：ComputerAidedDesign的英文缩写，CAD技术就是将计算机高速而精确的计算能力、大容量存储和处

4、理数据的能力与设计者的综合分析、逻辑判断及创造性思维结合起来，用以加快设计进程、缩短设计周期、提高设计质量的技术。机理模型法：采用由一般到特殊的推理演绎方法，对已知结构、参数的物理系统运用相应的物理定律，经过合理的分析简化而建立起来的描述系统各物理量动、静态变化性能的数学模型。统计模型法：采用由特殊到一般的逻辑、归纳方法，根据一定数量的在系统运行中实测的物理数据，运用统计规律、系统辨识等理论合理估计出反映系统各物理量相互制约的关系的数学模型。系统零极点：使系统传递函数分母为零的根叫做系统的极点；使系统传递函数分子为零的根叫做系统的零点。微分方程：表示未知函数的导数以及自变量之间关系的方程，就叫

5、做微分方程。紊流：当流体的雷诺系数大于2000时，流体的流态叫紊流。层流：当流体的雷诺系数小于2000时，流体的流态叫层流。实现问题：根据控制系统的传递函数描述求取其相应的状态空间描述。最小二乘曲线拟合：拟合出的曲线为数据点最小误差的二次和最小。插值：根据已知数据点，估计中间数据点的方法。线性系统：能用线性微分方程描述其输入和输出关系的称为线性系统。非线性系统：系统中包含非线性特性的元件，即为非线性系统。连续系统：系统的各部分的输入和输出信号都是连续函数的模拟量。离散系统：系统中一处或多处的信号以脉冲列或者数码的形式传递的系统。传递函数：初始条件为零的线性定常系统书出的拉普拉斯变换与输入的拉普

6、拉斯变换之比。鲁棒性：是指控制系统在一定（结构，大小）的参数摄动下，维持某些性能的特性。根据对性能的不同定义，可分为稳定鲁棒性和性能鲁棒性。DCS控制系统：即分布式控制系统，又称集散控制系统。是相对于集中控制系统而言的一种控制系统，它依靠网络将控制系统的各个部分连接起来。非最小相位系统：控制系统传递函数的零点或极点具有正实部，或有延迟环节。FFT：即快速傅氏变换，是离散傅氏变换的快速算法，它是根据离散傅氏变换的奇、偶、虚、实等特性，对离散傅立叶变换的算法进行改进获得的。Z变换：离散控制系统的拉普拉斯变换。与连续系统的拉氏变换相对应。广义Z变换：一种能够将两个采样时刻间的信息保存的Z变换方法，是

7、扩展的Z变换。时域分析：以时间为横坐标，利用关于时间的函数曲线分析系统特性的方法叫时域分析。频域分析：以频率为横坐标，利用关于频率的函数曲线分析系统特性的方法叫频域分析。根轨迹：当控制系统的某一参数发生变化时，特征方程的根在s平面上运动的轨迹叫根轨迹。超调：在系统的阶跃响应过程中，实际输出超过给定值的现象。超前/滞后校正：如果校正部分的传递函数的极点小于零点，则该校正部分在高于某一频率范围时会给系统一个正相移，称作相位超前校正；如果校正部分的传递函数的极点大于零点，则该校正部分会给系统一个负相移，称作相位滞后校正。目标函数：目标函数是衡量控制系统控制性能的评判工具，它的大小就说明了控制性能的优

8、劣。单纯形：所谓单纯形是指变量空间内最简单的形体。如在二维平面内，正三角形就是单纯形。ITAE准则：这是一种目标函数，它是指系统实际输出与系统给定的误差绝对值与时间的乘积的积分。可控性：如果在有限时间区间内存在容许控制向量u（t），能使此系统从初始状态转移到原点，则称该初始状态是能控的；若状态空间中所有的点都是能控的，那么就称该系统是能控的。可观性：对于连续线性系统，若存在一个时刻t1，使得系统以某非零状态为初始状态的输出恒为零，则称该非零状态不能观测。若系统的所有非零状态都不为不能观测，则系统是完全能观的。全状态反馈控制：是以系统的所有状态变量为反馈变量的控制形式。输出反馈控制：是以系统的输

9、出作为反馈变量的控制形式。时间最优控制：最优控制的一种，以时间为目标函数，实现调节时间最短的控制方法。相平面：以状态变量及其导数为坐标轴构成的平面叫相平面。答疑库问题1：什么是仿真？它所遵循的基本原则是什么？解答：仿真是建立在控制理论，相似理论，信息处理技术和计算技术等理论基础之上的，以计算机和其他专用物理效应设备为工具，利用系统模型对真实或假想的系统进行试验，并借助专家经验知识，统计数据和信息资料对试验结果进行分析和研究，进而做出决策的一门综合性的试验性科学。它所遵循的基本原则是相似原理。问题2：在系统分析与设计中仿真法与解析法有何区别？各有什么特点？解答：解析法就是运用已掌握的理论知识对控

10、制系统进行理论上的分析，计算。它是一种纯物理意义上的实验分析方法，在对系统的认识过程中具有普遍意义。由于受到理论的不完善性以及对事物认识的不全面性等因素的影响，其应用往往有很大局限性。仿真法基于相似原理，是在模型上所进行的系统性能分析与研究的实验方法。问题3：数字仿真包括那几个要素？其关系如何？解答：通常情况下，数字仿真实验包括三个基本要素，即实际系统，数学模型与计算机。由图可见，将实际系统抽象为数学模型，称之为一次模型化，它还涉及到系统辨识技术问题，统称为建模问题；将数学模型转化为可在计算机上运行的仿真模型，称之为二次模型化，这涉及到仿真技术问题，统称为仿真实验。问题4：为什么说模拟仿真较数

11、字仿真精度低？其优点如何？解答：由于受到电路元件精度的制约和容易受到外界的干扰，模拟仿真较数字仿真精度低，但模拟仿真具有如下优点：（1）描述连续的物理系统的动态过程比较自然和逼真。（2）仿真速度极快，失真小，结果可信度高。（3）能快速求解微分方程。模拟计算机运行时各运算器是并行工作的，模拟机的解题速度与原系统的复杂程度无关。（4）可以灵活设置仿真试验的时间标尺，既可以进行实时仿真，也可以进行非实时仿真。易于和实物相连。问题5：什么是CAD技术？控制系统CAD可解决那些问题？解答：CAD技术，即计算机辅助设计（ComputerAidedDesign）,是将计算机高速而精确的计算能力，大容量存储和

12、数据的能力与设计者的综合分析，逻辑判断以及创造性思维结合起来，用以快速设计进程，缩短设计周期，提高设计质量的技术。控制系统CAD可以解决以频域法为主要内容的经典控制理论和以时域法为主要内容的现代控制理论。此外，自适应控制，自校正控制以及最优控制等现代控制测略都可利用CAD技术实现有效的分析与设计。问题6：什么是虚拟现实技术？它与仿真技术的关系如何？解答：虚拟现实技术是一种综合了计算机图形技术，多媒体技术，传感器技术，显示技术以及仿真技术等多种学科而发展起来的高新技术。问题7：什么是离散系统？什么是离散事件系统？如何用数学的方法描述它们？解答：本书所讲的“离散系统”指的是离散时间系统，即系统中状

13、态变量的变化仅发生在一组离散时刻上的系统。它一般采用差分方程，离散状态方程和脉冲传递函数来描述。离散事件系统是系统中状态变量的改变是由离散时刻上所发生的事件所驱动的系统。这种系统的输入输出是随机发生的，一般采用概率模型来描述。若按模型分类，该系统属于那一类仿真系统？图中“混合计算机”部分在系统中起什么作用？与数字仿真相比该系统有什么优缺点？问题8：如图所示某卫星姿态控制仿真实验系统，试说明2）（1）3）指令与控制台PtH射頻模担S聽1-8卫1妻裁制访真试聲磁解答：（1）按模型分类，该系统属于物理仿真系统。（2）混合计算机集中了模拟仿真和数字仿真的优点，它既可以与实物连接进行实时仿真，计算一些复

14、杂函数，又可以对控制系统进行反复迭代计算。其数字部分用来模拟系统中的控制器，而模拟部分用于模拟控制对象。（3）与数字仿真相比，物理仿真总是有实物介入，效果逼真，精度高，具有实时性与在线性的特点，但其构成复杂，造价较高，耗时过长，通用性不强。问题9：数学模型的微分方程，状态方程，传递函数，零极点增益和部分分式五种形式，各有什么特点？解答：微分方程是直接描述系统输入和输出量之间的制约关系，是连续控制系统其他数学模型表达式的基础。状态方程能够反映系统内部各状态之间的相互关系，适用于多输入多输出系统。传递函数是零极点形式和部分分式形式的基础。零极点增益形式可用于分析系统的稳定性和快速性。利用部分分式形

15、式可直接分析系统的动态过程。问题10：数学模型各种形式之间为什么要互相转换？解答：不同的控制系统的分析和设计方法，只适用于特定的数学模型形式。问题11：控制系统建模的基本方法有哪些？他们的区别和特点是什么？解答：控制系统的建模方法大体有三种：机理模型法，统计模型法和混合模型法。机理模型法就是对已知结构，参数的物理系统运用相应的物理定律或定理，经过合理的分析简化建立起来的各物理量间的关系。该方法需要对系统的内部结构和特性完全的了解，精度高。统计模型法是采用归纳的方法，根据系统实测的数据，运用统计规律和系统辨识等理论建立的系统模型。该方法建立的数学模型受数据量不充分，数据精度不一致，数据处理方法的

16、不完善，很难在精度上达到更高的要求。混合法是上述两种方法的结合。问题12：控制系统计算机仿真中的“实现问题”是什么含意？解答：“实现问题”就是根据建立的数学模型和精度，采用某种数值计算方法，将模型方程转换为适合在计算机上运行的公式和方程，通过计算来使之正确的反映系统各变量动态性能，得到可靠的仿真结果。问题13：数值积分法的选用应遵循哪几条原则？解答：数值积分法应该遵循的原则是在满足系统精度的前提下，提高数值运算的速度和并保证计算结果的稳定。问题14：用matlab语言求下列系统的状态方程、传递函数、零极点增益、和部分分式形式的模型参数，并分别写出其相应的数学模型表达式c_s3+7s2+24+2

17、4S_s4+10s3+35s2+50s+24解答：（1）状态方程模型参数:编写matlab程序如下num=172424;den=l10355024;得到结果：A=-50001,C=172424;D=0-1O-35-50-241000藍+00100000100所以模型宵：産=U,y=172424X(2)零极点增益:编写程序num=l72424;den=l10355024;ZPK=tf2zp(num,den)得到结果Z=-2.7306+2.8531?-2.7306-2.85311,-1.5388P二-4?-3,-2?-lK=1(3)部分分式形式:编写程序num=l72424;den=l103550

18、24;RPH=residue(num,den)得到结果R=4.0000?-6.0000?2.0000,1.0000P=-4.0000,-3.0000,-2.0000,-1.0000H=z.4-621G(s)=+s+4s+3s+2s+1问题15：用matlab语言求下列系统的状态方程、传递函数、零极点增益、和部分分式形式的模型参数，并分别写出其相应的数学模型表达式2.25-5-1.25-0.5_2.25-4.25-1.254).2520.25-0.5-1.25-121.25-1.754).254).750解答：（1）传递函数模型参数:编写程序_A.ECD=tf2ss(num,den)A=2.25

19、-5-1.25-0.5-4.25-1.25-0.250.25-0.5-1.25-1-1.75-0.25-0.75；E=42201;C=0202;D=0;numden=ss2tf(A,B?C?D)得到结果num=0den=1.00004.000014.000022.000015.00004.00006.25005.25002.25004s3+14s2+22s+15s4+4s3+6.25s2+5.25s+2.25(2)零极点增益模型参数:编写程序A=2.25-5-1.25-0.5-4.25-1.25-0.250.25-0.5-1.25-1-1.75-0.25-0.75;B=422O1;C=0202

20、;D=0;Z?RK=SS2zp(A?B?C?D)得到结果Z=-1.0000+1.22471-1.0000-1.22471-1.5000P=-0.5000+0.86601-0.5000-0.86601-1.5000-1.5000K=4.0000表达式4(s+l-1.2247i)(s+1+1.22471)(s+0.5-0.866i)(s+0.5+0.866i)(s+1.5)(3)部分分式形式的模型参数：编写程序A=2.25-5-1.25-0.5-4.25-1.25-0.250.25-0.5-1.25-1-1.75-0.25-0.75;B=422O1;C=0202;D=OJ;numden=ss2tf

21、(A,B,C,D)R?RH=residue(num,den)中国地质大学（北京）继续教育学院2017年03课程考试中国地质大学（北京）继续教育学院2017年03课程考试中国地质大学(北京)继续教育学院2017年03课程考试第 页（共5页）第 页（共5页）第 页(共5页)得到结果尺二4.0000-0.00000.0000-2.309410.0000+2.30941P=-1.5000-1.5000-0.5000+0.86601-0.5000-0.86601H=42.309412.30941O(s)=+s+1.5s+0.5-0.866is+0.5+0.866/问题16:用欧拉法求下面系统的输出响应y

22、(t)在0t止，h=01时的数值。”=一”风)=1。要求保留4位小数，并将结果与真解比较。A+i=A+解答：欧拉法b曲)=耳(前向欧拉法，可以自启动)其几何意义：把f(t,y)在如九区间内的曲边面积用矩形面积近似代替。利用matlab提供的m文件编程，得到算法公式(1)m文件程序为h=O.l;蚯就函数的数值解为，)；显示中间的文字dispCy1);%同上尸1；fort=O:h:lm=y；disp(y);%显示y的当前值y=m-m*h;end保存文件q2.m在matalb命令行中键入q2得到结果函数的数值解为y=lO.9OOOO.8lOOO.729OO.656lO.59O5O.53l4O.478

23、3O.43O5O.3874O.3487(2)另建一个m文件求解y=Q在诈成1的数值(=Q是yr=-yy(0)=1的真解程序为h=D.l;蚯就函数的离散时亥懈那；disp(y-);fort=0:h:ldip(y);end保存文件q3.m在matalb命令行中键入q3函数的离散时刻解为y=10.90480.81870.74080.67030.60650.54880.49660.44930.40660.3679比较欧拉方法求堺与真值的差别欧拉10.90000.S1000.72900.65610.59050.53140.47S30.43050.38740.3487真值10.90480.81870.74

24、080.67030.60650.54880.49660.44930.40660.3679误差0-0.0048-0.0007-0.0118-0.0142-0.0160-0.0174-0.0183-0.0188-0.0192-0.0192显然误差与h2为同阶无穷小，欧拉法具有一阶计算精度，精度较低，但算法简单。问题17：用二阶龙格库塔法求解题16的数值解，并于欧拉法求得的结果比较。凡骚二片+临+禺）$=了（.4$焉=了（氐+鼠几+血対）解答：经常用到预报-校正法的二阶龙-格库塔法，也刃此方法可以自启动，具有二阶计算精度，几何意义：把f（t,y）在仏几区间内的曲边面积用上下底为比和九1、高为h的梯形

25、面积近似代替。利用matlab提供的m文件编程，得到算法公式(1)m文件程序为h=0.1?蚯就函数的数值解划)；disp(y=);fort=0:h:ldisp(y);kl=-y;k2=-(y+kl*h);y=y+(kl+k2)*h/2;end保存文件q4.m在matlab的命令行中键入q4显示结果为函数的数值解为尸10.90500.81900.4120.67DS0.60710.54940.49720.45000407203685比较欧拉法与二阶龙格-库塔法求解.（误差为绝对值）真值10.90430.81870.74030.67030.60650.54380.49660.44930.40660.

26、3679龙库10.90500.81900.74120.67080.60710.54940.49720.45000.40720.3685误差00.00020.00030.00040.00050.00060.00060.00060.00070.00060.0006明显误差为F得同阶无穷小，具有二阶计算精度，而欧拉法具有以阶计算精度，二阶龙格-库塔法比欧拉法计算精度高。问题18：用四阶龙格-库塔法求解题16数值解，并与前两题结果相比较。几+1二片+（样+2闰+2池+七J俎=了（氐必）禺=于应+于九+阳）解答：四阶龙格-库塔法表达式由=了（4+也丹+眺），其截断误差为同阶无穷小，当h步距取得较小时，误

27、差是很小的.（1）编辑m文件程序h=0.1;蚯就四阶龙格-库塔方法求解函数数值解为disp(y=);尸1；fort=0:h:ldisp(y);kl=-y;k2=-(y+kl*h/2);k3=-(y+k2*h/2);k4=-(y+k3*h);y=7+(kl+2*k2+2*k3+k4)*h/6;end保存文件q5.m在matlab命令行里键入q5得到结果四阶龙格-库塔方法求解函数数值解宵y=10.90480.81870.74080.67030.60650.54880.49660.44930.40660.36792）比较这几种方法：对于四阶龙格-库塔方法真值10.90480.81870.74080.

28、67030.60650.54880.49660.44930.40660.3679龙库10.90480.81870.74080.67030.60650.54880.49660.44930.40660.3679误差00000000000显然四阶龙格-库塔法求解精度很高，基本接近真值。三种方法比较可以得到精度（四阶）精度（二阶）精度（欧拉）问题19：已知二阶系统状态方程为写出取计算步长为h时，该系统状态变量X=%毛的四阶龙格-库塔法递推关系式。解答：四阶龙格-库塔法表达式A+1二丹+31+2怠+妬）耳-了（氐必）&=于&+品）忍=/+，片+妬）俎=+鼠儿+盹）所以状态变量的递推公式可以写作：中国地质

29、大学（北京）继续教育学院2017年03课程考试中国地质大学(北京)继续教育学院2017年03课程考试第 页（共5页）第 页(共5页)A=说_2112?B=英、可以写成X=AX+弧兀+1二兀+2妬+2&+札)麻=AXk+Bn则递推形式d=conv(a,b)；e=conv(d?c)e=1.00008.000031.990075.21000f=0005100;g=e+fg=1.00008.000031.990080.2100100.0000悯以上是计算闭环传递函数的特征塞项式悯p=roots(g)%计算特征务项式的根，就是闭环传递函数的极点P=-0.9987+3.0091i-0.99S7-3J009

30、1i-3.0013+0.9697i-3.0013-05697im=5100;z=roots(m)z=-20%计算零点时，可控标准型为：综上：当闭环传函形如俎护+.+跳_谆+如sK+.+住皆押+監(10.o001.00A=00101%bJ所以可控标准型是叮0=0问题21：用matlab语言编制单变量系统三阶龙格-库塔法求解程序，_010o-_0_0010%0+00010-100-80.21-31.99-8_1_程序入口要求能接收状态方程各系数阵(A,B,C,D),和输入阶跃函数r(t)=R*1(t);程序出口应给出输出量y(t)的动态响应数值解序列解答：m文件为:functiony=hs(A,B

31、,C,D,R,T,h)%T为观测时间,h为计算步长R为输入信号幅值%蚯就数值解为戈r=R;0；0;0;0;N=T/h;fort=l:N;kl=A*x+k2=A*(x+h*kl/3)+B*R;k3=A*(x+2*h*k2/3)+B*R;x+h*(kl+3*k3)/4;y(t)=C*x+D*R;end在命令行里键入A=E=C=D=R=T=h=y=hs(AB?C?D?艮T,h)得到结果o问题22：如课本上图2-27所示斜梁滚球系统，若要研究滚球在梁上的位置可控性，需首先建立其数学模型，已知力矩电机的输出转矩M与其电流i成正比，横梁为均匀可自平衡梁（即当电机不通电且无滚球时，横梁可处于日=0的水平状态

32、），是建立系统的数学模型，并给出简化后系统的动态结构图。解答：设球的质心到杆的距离为0，该系统为特殊情况下的球棒系统另令兀您&分别表示棒的惯量、球的质量和球的惯量。则球质心的位置和速度为兀=(xcossin0)叫=(vcos-xiysin?vsini9+xcos0其中x=v,8=5因而动能的移动部分为因而动能的移动部分为瓦m1肋；=；總（，+沧）球棒系统的旋转动能为因而，系统总的动能其中岛汨为常数。此系统的拉格朗日方程组为中国地质大学（北京）继续教育学院2017年03课程考试中国地质大学（北京）继续教育学院2017年03课程考试第 页（共5页）第 页（共5页）d（dTdT.4dt-T（）-=-

33、mSsin3,dt张唸di色（竺）竺=4胞呂co昉dtdede综合以上公式的系统的方程组为mAx-mx2+mgsin（0=0（右+mx2）3+2mxx3+mgxcos（0=id设系统在平衡点附近办0，cosl,sin&刘乩则系统方程可化为mAx+mg0=02（h+mx）O+mgx=kt对上式进行拉普拉斯变换并化简后可得到需。问题23：求解下列线性方程，并进行解的验证：解答：x=abx=0.49790.14450.0629-0.0813问题24：求解下列线性方程，并进行解的验证_57651-2496_71087234136681093X=3614457910435140_12345_1560_中

34、国地质大学(北京)继续教育学院2017年03课程考试中国地质大学(北京)继续教育学院2017年03课程考试第 页(共5页)第 页(共5页)解答：a=5765171087268109357910412345765710878109579101234b=2496136361441401560b=2496136361441401560 x=abx=1.00001.00001.00001.00001.00004.00004.00004.00004.00004.0000123斗5问题25：进行下列计算，给出不使用for和while等循环语句的计算方法。I解答：根据等比数列求和方法，在利用matlab中的

35、m文件，编写程序求解。M文件为n=64;q=2;k=(l-qM)/(l-q);disp(k的值为);disp(k);中国地质大学（北京）继续教育学院2017年03课程考试中国地质大学（北京）继续教育学院2017年03课程考试第 页（共5页）第 页（共5页）保存文件ql.m。在matlab命令框中输入k的值为1.8447e+019问题26:求出y=x*sin(x)在盐100条件下的每个峰值a=size(y)a=l1001b=(y(2:1000)y(l:999)&(y(2:1000)y(3:1001);aWmd(b=l);解答：disp(y(at)问题27:绘制下面的图形：解答：y=sm(l./t

36、);plot(t?y)gridonxlabelCt1)ylabelCy1)title(ysmCl/t)1)%注盍是环是/%问题28:绘制下面的图形：1-cos3(7f)-ltgridonxlabelft1)ylabelCy1)解答：title(7=1yi=polyval(p?xi);plot(x?y?xi;yi)gridon红色：采样曲线绿色：拟合曲线问题30：使用simulink求解具有如下闭环传递函数的系统的阶跃响应欧町=/+8s3+36s2+40s+10解答：在simulink模型窗口中建立如下模型，键入该题的传递函数。start后，观察scope中的仿真波形如下：占eq必口與e6s3+

37、26s2+6s+20，试分析问题31:已知系统的闭环传递函数?+3?+4s2+2s+2该系统的稳定性解答：由稳定性判据：当闭环传递函数所有极点都位于虚轴左半平面时，该系统稳定。传递函数的特征方程为：+菇+加+春+2,解此方程，得到特征根，即闭环极点p=l3422;在matlab命令行里键入得到r=-1.4734+1.02561-1.4734-1.02561-0.0266+0.78731-0.0266-0.78731闭环极点的实部都小于零，即都位于虚轴左半平面，所以系统稳定问题32：某小功率随动系统动态结构如图所示，已知：=0.01,爲=0.0疋=1,疋严300，瓦=1,瓦=0.02.若系统输入

38、分别为务=临），试用simulink分析系统的输出监打？解答：问题33:设典型闭环结构控制系统如图所示，当阶跃输入幅值应=20时，用sp4_l.m求取输出讥的响应。解答：用曲4.m求解过程如下：在MATLAB语言环境下，输入以下命令语句a=0.0160.8643.273.421;b=3025;X0=0000;V=2,n=4;T0=0;Tf=10;%系统状态向量初值为零%反愦系数v=2h=0.01;R=20;%仿真步长h=0.01,阶跃输入幅值R=20sp4_lplot(t,y)%调用聊4_l.ni函数运行结果为:附：sp4_1.m函数为b=b/a(1);a=a/a(1);A=a(2:n+1);

39、A=rot90(rot90(eye(n-1,n);-fliplr(A);B=zeros(1,n-1),1;m1=length(b);C=fliplr(b),zeros(1,n-m1);Ab=A-B*C*V;X=X0;y=0;t=T0;N=round(Tf-T0)/h);fori=1:NK1=Ab*X+B*R;K2=Ab*(X+h*K1/2)+B*R;K3=Ab*(X+h*K2/2)+B*R;K4=Ab*(X+h*K3)+B*R;X=X+h*(K1+2*K2+2*K3+K4)/6;y=y,C*X;t=t,t（i）+h;end碍，兀拓扑连结关系，可写出每个环节输入碍受哪些环问题34:系统结构图如图

40、，写出该系统的联结矩阵莊和吒，并写出联结矩阵非零元素阵解答：根据图中节输出的影响，现列如入下:=凡&4=”_必吗幵遐io=y?把环节之间的关系和环节与参考输入的关系分别用矩阵表示出来，U=WY+Wpl10000*+0000如_0_即莊=L00000000100001000010000100000000000000000000000000-10000000000-10000000000000-110000010000001000001000012234681001-12118-11-161617171问题35：若系统为图4-5b双输入-双输出结构，试写出该系统的联接矩阵莊，吒，说明应注意什么？解

41、答：根据图4-5b中碍，兀拓扑连结关系，可列写如下关系式:转换成矩阵形式为所以联接矩阵莊=此时应注意输入联接矩阵变为6滋型。问题36：求图4-49非线性系统的输出响应y（t）,并与无非线性环节情况进行比较。解答：（1）不考虑非线性环节影响时，求解过程如下：1）先将环节编号标入图中。2）在MATLAB命令窗口下，按编号依次将环节参数输入P阵；P=0.110.51;01200;2110;10110;3）按各环节相对位置和联接关系，有联接矩阵如下001001000-1000010所以非零元素矩阵014-12131WIJ=101;14-1;211;321;431;4）由于不考虑非线性影响，则非线性标志

42、向量和参数向量均应赋零中国地质大学(北京)继续教育学院2017年03课程考试中国地质大学(北京)继续教育学院2017年03课程考试第 页(共5页)第 #页(共5页)值；Z=0000;S=0000；5)输入运行参数：开环截至频率厶零约为1，故计算步长h取经验公Z2h=0.01;L1=25;n=4;T0=0Tf=20;nout=4;Y0=10;sp4_4;plot(t,y,r)holdon运行结果如图中红色实线所示。(2)考虑非线性环节N影响时，只需将非线性标志向量Z和参数向量S的相应分量正确输入即可。在MATLAB命令窗口中输入下列语句：Z=4000;S=5000;%第一个线性环节后有饱和非线性

43、，参数值为5。sp4_4;plot(t,y,-)运行结果如图中蓝色虚线所示。从图中可以清楚的地看出，饱和非线性环节对线性系统输出响应的影响。中国地质大学（北京）继续教育学院2017年03课程考试中国地质大学（北京）继续教育学院2017年03课程考试第 页（共5页）第 页（共5页）附：sp4_4函数为:A=P(:,1);B=P(:,2);C=P(:,3);D=P(:,4);m=length(WIJ(:,1);W0=zeros(n,1);W=zeros(n,n);fork=1:mif(WIJ(k,2)=0);W0(WIJ(k,1)=WIJ(k,3);elseW(WIJ(k,1),WIJ(k,2)=

44、WIJ(k,3);end;end;fori=1:nif(A(i)=0);FI(i)=1;FIM(i)=h*C(i)/B(i);FIJ(i)=h*h*(C(i)/B(i)/2;FIC(i)=1;FID(i)=0;if(D(i)=0);FID(i)=D(i)/B(i);elseendelseFI(i)=exp(-h*A(i)/B(i);FIM(i)=(1-FI(i)*C(i)/A(i);FIJ(i)=h*C(i)/A(i)-FIM(i)*B(i)/A(i);FIC(i)=1;FID(i)=0;if(D(i)=0);FIC(i)=C(i)/D(i)-A(i)/B(i);FID(i)=D(i)/B(i

45、);elseendendendY=zeros(n,1);X=Y;y=0;Uk=zeros(n,1);Ubb=Uk;t=T0:h*L1:Tf;N=length(t);fork=1:N-1fori=1:L1Ub=Uk;Uk=W*Y+W0*Y0;fori=1:nif(Z(i)=0)if(Z(i)=1)Uk(i)=satu(Uk(i),S(i);endif(Z(i)=2)Uk(i)=dead(Uk(i),S(i);endif(Z(i)=3)Uk(i),Ubb(i)=backlash(Ubb(i),Uk(i),Ub(i),S(i);endendendUdot=(Uk-Ub)/h;Uf=2*Uk-Ub;X

46、=FI.*X+FIM.*Uk+FIJ.*Udot;Yb=Y;Y=FIC.*X+FID.*Uf;fori=1:nif(Z(i)=0)if(Z(i)=4)Y(i)=satu(Y(i),S(i);endif(Z(i)=5)Y(i)=dead(Y(i),S(i);endif(Z(i)=6)Y(i),Ubb(i)=backlash(Ubb(i),Y(i),Yb(i),S(i);endendendendy=y,Y(nout);end附：饱和非线性函数satu.m为functionUc=satu(Ur,S1)if(abs(Ur)=S1)if(Ur0)Uc=S1;elseUc=-S1;endelseUc=Ur

47、;end问题37：采样控制系统如图所示，编写程序实现对该系统的仿真分析（提示：连续部分按环节离散化考虑）图中，训应+勻J为典型数字PID控制器；宀=065为比例系数；爲=0.7为积分时间常数；兀=0.2为微分时间常(卑+1)(那+1)为具有纯滞后特性的典型二阶控制对象；数；解答：在控制对象前引入零阶保持器，将连续环节部分按环节离散化:1厂冷严设忑,为简化运算及编程，取耳为T的整数倍。2口）=对上式进行Z逆变换，得到T3T3=(1一-aT)U(_k-y-T)+严+必严)2)+2s(k_亍_1)_由此可编写仿真程序。在MATLAB命令窗口中输入下列语句KP=0.65;TI=0.7;TD=0.2;T

48、1=0.3;a=1/T1;T3=0.4;T=0.1;h=0.001;Tf=10;hh编写M脚本文件，存为hh.m。%离散化后各参数为：A=1-a*h*exp(-a*h)-exp(-a*h);B=exp(-2*a*h)-exp(-a*h)+a*h*exp(-a*h);C=2*exp(-a*h);D=exp(-2*a*h);P=KP*(1+T/TI+TD/T);H=KP*(1+2*TD/T);M=KP*TD/T;%系统初始值为：E=zeros(1,3);U=zeros(1,2+T3/T+1);Y=zeros(1,2+T3/h+1);R=1;yk=0;yt=0;t=0;%仿真迭代运算：forK1=1

49、:Tf/Tek=R-Y(1);E=ek,E(1:2);uk=P*E(1)-H*E(2)+M*E(3)+U(1);U=uk,U(1:(2+T3/T);forK2=1:T/h-D*Y(T3/h+2);yk=A*U(T3/T+1+1)+B*U(T3/T+2+1)+C*Y(T3/h+1)Y=yk,Y(1:(2+T3/h);endyt=yt,yk;t=t,K1*T;end%输出波形：plot（t,yt）运行结果为:此题可以用SIMULINK仿真进行验证建立SIMULINK仿真模型:运行结果为:中国地质大学(北京)继续教育学院2017年03课程考试中国地质大学(北京)继续教育学院2017年03课程考试第

50、页(共5页)第 #页(共5页)G(s)H(s)=習巴问题38:设控制系统的开环传递函数为试画出该系统的根轨迹。解答：在Matlab窗口中输入下列命令：num=11;a=10;b=1-1;c=1416;d=conv(a,b);den=conv(d,c);rlocus(num,den)gridon可得到系统的根轨迹如下图所示：中国地质大学（北京）继续教育学院2017年03课程考试中国地质大学（北京）继续教育学院2017年03课程考试第 页（共5页）第 页（共5页）d总aieii_beui问题39:某反馈控制系统的开环传递函数为K龙+呗/+牝+20）试绘制其根轨迹。解答：在MATLAB命令窗口中输入

51、下列命令num=1;den=conv(conv(1,0,1,4),1,4,20);rlocus(num,den)gridon运行结果为：问题40：已知某系统传递函数为试绘制其伯德图。Eg200佔+100)(2.5s+100)(+2x0.3s+200)Ws)解答：分子分母同乘100*200得到在Matlab窗口中输入下列命令k=80*200;num=1100;a=2.5100;b=(1/200)2*0.3200;den=conv(a,b);w=logspace(-1,1,100);m,p=bode(k*num,den,w);subplot(2,1,1);semilogx(w,20*log10(m

52、);grid;xlabel(Frequency(rad/sec);ylabel(Gain(dB);subplot(2,1,2);semilogx(w,p);grid;xlabel(Frequency(rad/s);ylabel(Phase(deg);可绘制该系统的伯德图如下所示。15858-3-93s-151010Frequerciiiraiis)1037.8510Freqpnui.itrad点er?i问题41：设控制系统具有如下的开环传递函数恥吩心+恥+5)试求取当K=10时的相角裕度和幅值裕度，并画出其伯德图。解答：在MATLAB命令窗口中输入下列命令：k=10;num=1;den=pol

53、y(0,-1,-5);m,p,w=bode(num,den);subplot(2,1,1);semilogx(w,20*log10(m);gridon;ylabel(Gain(dB);subplot(2,1,2);semilogx(w,p);gridon;xlabel(Frequency(rad/sec);ylabel(Phase(deg);gm,pm,wcg,wcp=margin(m,p,w)这里gm,wcg为幅值裕度值与相应的频率pm,wcp为相角裕度值与相应的频率，运行结果为：gm=30.0000,pm=76.5911,wcg=2.2361,wcp=0.1960。因此，系统的幅值裕度和相

54、角裕度分别为30dB和系统的伯德图如下所示。问题42：已知某单位反馈系统开环传递函数为_s(0.15+1)(0.02s+1)(0.01s+1)(0.005s+1)若性能指标要求如下：丫二仃尼=2g%=13d,试确定校正装置。解答：由速度误差常数的要求可求得K=200，则可绘制阿的波德图。在Matlab窗口中输入下列命令：k=200*10*50*100*200;num=1;den=poly(0-10-50-100-200);w=logspace(-1，2，200);m，p=bode(k*num，den，w);subplot(2，1，1);semilogx(w，20*log10(m);grid;y

55、label(Gain(dB);subplot(2，1，2);semilogx(w，p);grid;xlabel(Frequency(rad/s);ylabel(Phase(deg);总uraCD从图中可以看出，相角裕度为-53度，剪切频率为37.7rad/s，这显然EapvEEqd不能满足性能指标要求。由于补偿角为=5?+45=98，单独使用超前校正装置或是滞后校正装置都不能很好的达到设计要求，因此这里考虑使用串联相位滞后超前校正。设校正装置的传递函数为二（砂+1）（加+1）“（0五&+1）（殆/0+1）1=Ur1乌=1用2，O按照给定的要求选定系统剪切频率为=15rad/s，为使系统稳定，应

56、保证开环对数幅频特性以-20dB/dec的斜率穿越频率轴，即中频段的斜率为-20dB/dec。为使校正后系统的开环增益不低于200，校正前后系统的低频段特性应保持一致。中频段与低频段之间用斜率为-40dB/dec的直线连接，连接线与中频段相交的交接频率叭与气不宜离的太近，取=/10=1.5rad/s；同时取口航，则巧=1/1.5=0.6667=1/10=0.15。最后来确定Q的取值。由于校正后的剪切频率为=15,则对应I=15rad/S处的增益为201g_20;而未校正系统在I=15rad/S时的增益201g=42.5为L5。两者相减就得到串联校正装置在1.5rad/s111111111111

57、i1ii11ii1111III111111II111111111111111111i1ii11ii11I111|111111II1111111111Lu丄叫-.丄Lll-L丄Llli-!-LUL-1-_Jill11111iniii1i1ii1I*Ii|i|11ii11ii1i11111ii111ini1iiiiiii1ii|i|11ii11ii1ii111111iliiiiiii1illiilii11iiiil111iiiiiiTi-i-tri+in-nii-u士巴_I-+TriTlli1iiliiMr+l-iii11iiii111-卜-*mnn-t-i-mtit+8D-U-46运行程序后得到

58、校正后系统的剪切频率畋i4Q49,丫二沁嘶,这些都满足期望的性能指标，同时系统伯德图如下所示。Frequercy（rad/sei：）校正后闭环系统的阶跃响应如下图所示，可见系统具有较好的动态性能与稳态性能。问题43：某过程控制系统如图5-71所示，试设计PID调节器参数使该系统动态性能达到最佳。解答：本题选用ITSE准则为目标函数。首先，编写M函数文件，存为optm.m。functionss=optm(x)globalkp;globalki;globalkd;kp=x(1);ki=x(2);kd=x(3);tt,xx,yy=sim(opt,40,);ss=yy(end);其次，输入恰当的模块参

59、数，使用SIMULINK建立仿真模型文件，存为opt.mdl。在MATLAB命令窗口中输入下列指令:globalkp;globalki;globalkd;result=fminsearch（optm,2,1,2）运行结果为：result=2.69140.46102.1267,即kp=2.6914,ki=0.4610；kd=2.1267。用SIMULINK仿真验证控制器的控制效果，得到阶跃响应曲线如下图所示（蓝色实线所示），与初始值kp=2,ki=1,kd=2时的阶跃响应比较（图中红色虚线所示），动态性能得到了很好的改善。仿真结果如下图所示。问题44：试采用smith预估控制方法对题5-6所述系

60、统进行重新设计并用仿真的方法分析滞后参数变化对系统动态性能的影响。解答：设为过程控制通道传递函数，其中：为过程中不包含纯滞后部分的传递函数，为控制器的传递函数。在本题中,莎a二幻-技+?丁-士-行，P。不使用Smith预估控制方法时，系统闭环传递由于在特征方程中引入了它：“项，使得闭环系统的品质大大恶化。Smith预估补偿控制方法的实质是通过预估补偿装置的引入，将实际工业过程中不可分割的与y在传函形式上分开，并以为过程控制通道的传递函数，以其输出信号作为反馈信号，以此来改善控制质量。预估补偿装置的传递函数为匸匸一匕，使用后控制系统的框图如下，系统闭环传递函数为:1、-八丄-丿可见，经过Smit