控制系统数字仿真与CAD第一二章习题答案.docx

1、控制系统数字仿真与CAD第一二章习题答案1-1什么是仿真？它所遵循的基本原则是什么？答：仿真是建立在控制理论，相似理论，信息处理技术和计算技术等理论基础之上的，以计算机和 其他专用物理效应设备为工具，利用系统模型对真实或假想的系统进行试验，并借助专家经验 知识，统汁数据和信息资料对试验结果进行分析和研究，进而做出决策的一门综合性的试验性 科学。它所遵循的基本原则是相似原理。1-2在系统分析与设计中仿真法与解析法有何区別？各有什么特点？答：解析法就是运用已掌握的理论知识对控制系统进行理论上的分析，il算。它是一种纯物理意义 上的实验分析方法，在对系统的认识过程中具有普遍意义。由于受到理论的不完善

2、性以及对事 物认识的不全而性等因素的影响，其应用往往有很大局限性。仿真法基于相似原理，是在模型上所进行的系统性能分析与研究的实验方法。1-3数字仿真包括那几个要素？其关系如何？答：通常情况下，数字仿真实验包括三个基本要素，即实际系统，数学模型与让算机。由图可见， 将实际系统抽象为数学模型，称之为一次模型化，它还涉及到系统辨识技术问题，统称为建模 问题：将数学模型转化为可在计算机上运行的仿真模型，称之为二次模型化，这涉及到仿真技 术问题，统称为仿真实验。1-4为什么说模拟仿真较数字仿真精度低？其优点如何？ o答：由于受到电路元件精度的制约和容易受到外界的下扰，模拟仿真较数字仿真精度低但模拟仿真具

3、有如下优点：（1） 描述连续的物理系统的动态过程比较自然和逼真。（2） 仿真速度极快，失真小，结果可信度髙。（3） 能快速求解微分方程。模拟汁算机运行时0运算器是并行工作的，模拟机的解题速度与原系统的复杂程度无关。（4） 可以灵活设置仿真试验的时间标尺，既可以进行实时仿真，也可以进行非实时仿真。（5） 易于和实物相连。1-5什么是CAD技术？控制系统CAD可解决那些问题？答：CAD技术，即计算机辅助设计（Computer Aided Design），是将计算机高速而精确的计算能力, 大容量存储和数据的能力与设讣者的综合分析，逻辑判断以及创造性思维结合起来，用以快速 设计进程，缩短设计周期，提髙

4、设计质量的技术。控制系统CAD可以解决以频域法为主要内容的经典控制理论和以时域法为主要内容的现代控制 理论。此外，自适应控制，自校正控制以及最优控制等现代控制测略都可利用CAD技术实现有 效的分析与设计。1-6什么是虚拟现实技术？它与仿真技术的关系如何？答：虚拟现实技术是一种综合了计算机图形技术，多媒体技术，传感器技术，显示技术以及仿真技 术等多种学科而发展起来的高新技术。1-7什么是离散系统？什么是离散事件系统？如何用数学的方法描述它们？答：本书所讲的“离散系统”指的是离散时间系统，即系统中状态变量的变化仅发生在一组离散时 刻上的系统*它一般采用差分方程.离散状态方程和脉冲传递函数来描述。离

5、散事件系统是系统中状态变量的改变是由离散时刻上所发生的事件所驱动的系统。这种系统 的输入输出是随机发生的，一般采用概率模型来描述。1-8如图1-16所示某卫星姿态控制仿真实验系统，试说明：（1） 若按模型分类，该系统属于那一类仿真系统？（2） 图中“混合汁算机”部分在系统中起什么作用？（3） 与数字仿真相比该系统有什么优缺点？答：（1）按模型分类，该系统属于物理仿真系统“（2）混合计算机集中了模拟仿真和数字仿真的优点，它既可以与实物连接进行实时仿真，计算 一些复杂函数，又可以对控制系统进行反复迭代讣算。其数字部分用来模拟系统中的控制器， 而模拟部分用于模拟控制对象。（4） 与数字仿真相比，物理

6、仿真总是有实物介入，效果逼真，精度高，具有实时性与在线性的特点,但其构成复杂，造价较髙，耗时过长，通用性不强。2-1思考题：（1） 数学模型的微分方程，状态方程，传递函数，零极点增益和部分分式五种形式，各有什么特点？（2） 数学模型各种形式之间为什么要互相转换？（3） 控制系统建模的基本方法有哪些？他们的区别和特点是什么？（4） 控制系统计算机仿真中的“实现问题”是什么含意？（5） 数值积分法的选用应遵循哪几条原则？答:（1）微分方程是直接描述系统输入和输出量之间的制约关系，是连续控制系统其他数学模型表 达式的基础。状态方程能够反映系统内部各状态之间的相互关系，适用于多输入多输出系统匚传递 函

7、数是零极点形式和部分分式形式的基础。零极点增益形式可用于分析系统的稳左性和快速性。利 用部分分式形式可直接分析系统的动态过程。（2）不同的控制系统的分析和设计方法，只适用于特左的数学模型形式。（3控制系统的建模方法大体有三种：机理模型法，统汁模型法和混合模型法。机理模型法就 是对已知结构，参数的物理系统运用相应的物理立律或左理，经过合理的分析简化建立起来的各物 理量间的关系。该方法需要对系统的内部结构和特性完全的了解，精度髙。统讣模型法是采用归纳 的方法，根据系统实测的数据，运用统讣规律和系统辨识等理论建立的系统模型。该方法建立的数 学模型受数据咼不充分，数据精度不一致，数据处理方法的不完善，

8、很难在精度上达到更髙的要求。 混合法是上述两种方法的结合。（4） “实现问题就是根拯建立的数学模型和精度，采用某种数值计算方法，将模型方程转换 为适合在计算机上运行的公式和方程，通过计算来使之正确的反映系统各变量动态性能，得到可靠 的仿真结果。（5） 数值积分法应该遵循的原则是在满足系统精度的前提下，提髙数值运算的速度和并保证计 算结果的稳左。22用matlab语言求下列系统的状态方程.传递函数.零极点增益、和部分分式形式的模型参数, 并分别写出其相应的数学模型表达式：?+7?+245 + 24G(S) Z + 1053 +35?+505 + 242.25 -5 -1.25 -0.542.25

9、 -4.25 -1.25 -0.25x +20.25 -0.5 -1.25 -121.25 -1.75 -0.25 -0.750错误!未找到引用源。X二y=0 2 0 2 X(1)解：（1）状态方程模型参数:编写matlab程序如下num二1 7 24 24; den=l 10 35 50 24;ABC D=tf2ss(num, den)-10得到结果：A二100-10所以模型为：X二num=Ll10u, y=l10-350-350,c=l7 24 24,D 二07 24 24 X（2）零极点增益:编写程序7 24 24; den=l 10 35 50 24; Z P K=tf2zp(num,

10、 den)得到结果 Z二-2.7306 + 2.8531 f -2. 7306 - 2. 8531i ,-1.5388P二-4, -3 ,-2 ,-lK=1(3)部分分式形式:编写程序 num二1 7 24 24; den=l 10 35 50 24; R P H二residue(num, den)得到结果 R二 4. 0000 , -6. 0000, 2. 0000, 1. 0000P二 -4. 0000, -3. 0000 , -2. 0000 ,-1. 0000H A二225 -5 -1.25 -0.52. 25 -4.25 -1.25 -0. 250. 25 -0. 5 -1. 25

11、-11.25-1.75 -0. 25 -0. 75: B二4 2 2 0; C二0 2 0 2; D二0; num den =ss2tf (A, B, C, D)得到结果num = 0 4. 0000 14. 0000 22. 0000 15. 0000den =1. 0000 4. 0000 6. 2500 5. 2500 2. 2500 4 s3 + 14 s2 + 22 s + 15G($) = -7 Z s4 +4s3+6.25 s2 + 5.25 s + 2.25(2)零极点增益模型参数:编写程序 A=2. 25 -5 -1.25 -0.52. 25 -4.25 -1.25 -0.

12、250. 25 -0. 5 -1 25 -11.25-1.75 -0. 25 -0.75： B二4 2 2 0; C二0 2 0 2; D二0; Z, P, K=ss2zp(A, B, C, D)得到结果 Z =-1. 0000 + 1. 2247i -1. 0000 - 1. 2247i -1. 5000P二-0. 5000 + 0. 8660i -0.5000 - 0. 8660i -1. 5000-1. 5000K = 4. 0000表达式 G(沪 空亠2247i)(s+l+1.2247i)(s+O.5-O.866i)(s+O.5+O.866i)(s+1.5)(3)部分分式形式的模型参数

13、：编写程序 A二2. 25 -5 -1.25 -0.52. 25 -4.25 -1.25 -0. 250. 25 -0. 5 -1. 25 -11.25-1.75 -0. 25 -0.75； B二4 2 2 0; C二0 2 0 2; D二0; num den=ss2tf (A, B, C, D) R, P, H =residue (num, den)得到结果 R = 4. 0000 -0. 0000 0. 0000 - 2. 3094i 0. 0000 +2. 3094iP = -1. 5000 -1. 5000 -0. 5000 + 0. 8660i -0. 5000 -0.8660i H

14、 二、 4 2.3094/ 2.3094ZG(s) = + 5 + 1.5 5+ 0.5 - 0.866/ 5 + 0.5 + 0.866/2-3.用欧拉法求下面系统的输出响应y(t)在OWtWl上，h=0.1时的数值。y=-2(0)=1要求保留4位小数，并将结果与真解y(r)= h比较。儿+L儿+*/(，儿)解：欧拉法y=f(tkiyk) (前向欧拉法，可以自启动)其几何意义：把f(t,y)在山，儿yM = o区间内的曲边而积用矩形面穆近似代替。利用matlab提供的m文件编程，得到算法公式。如卞所示(1)m文件程序为h二0.1;dispC函数的数值解为)；显示中间的文字disp( y二);

15、%同上y二 1；for t=O:h：1m=y;disp(y); %显示y的当前值y=m-m*h;end保存文件q2m在matalb命令行中键入 q2得到结果函数的数值解为y= 1 0.9000 0.8100 0.7290 0.6561 0.5905 0.5314 0.4783 0.4305 0.38740. 3487(2)另建一个m文件求解y = k在te0,l的数值 (是)r = y,y(o)= l的真解%)程序为h=0. 1; dispC函数的离散时刻解为)； dispC y二);for t=O:h:1y=exp(-t);disp(y);end 保存文件q3. m在matalb命令行中键入

16、 q3函数的离散时刻解为y= 1 0.9048 0.8187 0.7408 0.6703 0.6065 0.5488 0.4966 0.4493 0.40660. 3679比较欧拉方法求解与真值的差别欧 拉10. 90000.81000. 72900. 65610. 59050. 53140. 47830. 43050. 38740. 3487真值10. 90480. 81870. 74080. 67030. 60650. 54880. 49660. 44930. 40660. 3679误0-0. 004-0. 000一一一一一一-0. 019-0. 019差870.01180. 01420.

17、 01600.01740.01830.018822显然误差与/为同阶无穷小，欧拉法具有一阶计算精度，精度较低，但算法简单。2-4用二阶龙格库塔法求解2-3的数值解,并于欧拉法求得的结果比较。儿+1 = ”+*(+)解：我们经常用到预报二狡疋法的二阶龙格库塔法,儿) 此方法可以自启k2 = f(tk + 九” + hkj动，具冇二阶计算精度，几何意义:把f(t,y)在也，儿区间内的曲边面积用上下底为人和/叶、 髙为h的梯形面积近似代替。利用mat lab提供的m文件编程，得到算法公式。如下所示(1)m文件程序为h=0. 1;dispC函数的数值解为)；disp( y=，);y二 1;for t=

18、0:h:1disp(y);kl=-y;k2=- (y+kl*h); y=y+(kl+k2)*h/2;end保存文件q4mmat lab的命令行中键入 q4 显示结果为函数的数值解为y= 1 0.9050 0.8190 0.7412 0.6708 0.6071 0.5494 0.4972 0.4500 0.40720. 3685(2) 比较欧拉法与二阶龙格-库塔法求解.(误差为绝对值)一真值10. 90480.81870. 74080. 67030. 60650. 54880. 49660. 44930. 40660. 3679龙 库10. 90500.81900. 74120. 67080.

19、60710. 54940. 49720. 45000. 40720. 3685误 差00. 00020. 00030. 00040. 000L0. 00060. 00060. 00060. 00070. 00060. 0006明显误差为/F得同阶无穷小，具有二阶汁算精度，而欧拉法具有以阶计算精度，二阶龙格-库塔法 比欧拉法讣算精度高。2-5.用四阶龙格-库塔法求解题2-3数值解，并与前两题结果相比较。=儿+?&+2人+2人+心) o解：四阶龙格-库塔法表达式其截断误差为斥同阶无穷小,k=Mtk，yJ込=Wk +-o +i)免=+-O- + 牙斤2)k4=f(tk+h,yk+hk3) 当h步距取

20、得较小时，误差是很小的.(1)编辑m文件程序h二0.1;dispC四阶龙格-库塔方法求解函数数值解为);disp( y二);y二 1；for t=0:h:1disp(y);kl=-y;k2=-(y+kl*h,/2);k3=-(y+k2*h/2);k4=-(y+k3*h);y=y+(kl+2*k2+2*k3+k4)*h/6;end 保存文件q5. m在mat lab命令行里键入 q5得到结果四阶龙格库塔方法求解函数数值解为y= 1 0.9048 0.8187 0.7408 0.6703 0.6065 0.5488 0.4966 0.4493 0.40660. 3679(2)比较这几种方法：对于四

21、阶龙格-库塔方法一真值10. 90480.81870. 74080. 67030. 606L00. 54880. 49660. 44930. 40660. 3679龙 库10. 90480.81870. 74080. 67030. 606L0. 54880. 49660. 44930. 40660. 3679误 差00000000000显然四阶龙格-库塔法求解精度很高，基本接近真值。三种方法比较可以得到 精度(四阶)精度(二阶)稱度(欧拉)2-6.已知二阶系统状态方程为“I】a2 iri+blw；召(0)-xio5 ” ILx-4禺(0)写出取计算步长为h时，该系统状态变量X=片,兀的四阶龙格

22、-库塔法递推关系式。儿严儿+纟& + + 2/+心)o解：四阶龙格库塔法表达式、k? = f Ok +了，儿 +-i)k、= f(fk +儿 +_2)人=/(+九儿+g) 所以状态变疑的递推公式可以写作：xix2x知严人+&+2+2+ o& = AXk + Bii =A(X*+Ei*/i/2) + B“ k3 = A(Xk +k2* h! 2) + Bu k4 = A(Xk +k3*h) + Bu可以写成X=AX+Bu,B=.如 12则递推形式彳21 “22.2-7单位反懺系统的开环传递函数已知如下G(s) =5s + 1005(5 + 4.6)(52 +3.4.? +16.35)用matla

23、b语句.函数求取系统闭环零极点.并求取系统闭环状态方程的可控标准型实现。 解：已知开环传递函数，求得闭环传递函数为G(s) = J 二四 5(5 + 4.6)(52 +3.45 + 16.35) + 5$ +100 在mat lab命令行里键入 a=l 0; b二1 4.6; c二1 3.4 16. 35; d=conv(a,b): e=conv(d,c)e = 1.0000 8.0000 31.9900 75.2100 0 f二0 005 100; g=e+fg = 1. 0000 8. 0000 31. 9900 80. 2100 100. 0000%以上是计算闭环传递函数的特征多项式 P

24、=roots(g) %计算特征多项式的根，就是闭环传递函数的极点%p = -0. 9987 + 3. 00911-0.9987 - 3.0091i-3. 0013 + 0. 9697i-3. 0013 - 0. 9697i m=5 100; z=roots(m)z = -20 %计算零点综上；当闭环传函形如G(s)= 戶 _九 时,可控标准型为:s +qs + + 4_$ + %0 1 0 . 0,0 0 100 0 0104 dJ丄C= bnb ;D = 010001000-100-80.21-31.99所以可控标准型是r = -ioo 5 o 0勺 +0uX32-8用matlab语言编制单

25、变量系统三阶龙格库塔法求解程序，程序入口要求能接收状态方程各系 数阵(A,B,C,D),和输入阶跃函数r(t)=R*l(t);WT出口应给出输出量y (t)的动态响应数值解 序列凡,，儿。解:m文件为:function y=hs(A,B, C,D, R, T, h) %T为观测时间,h为计算步长,R为输入信号幅值% dispC数值解为)；y=0；r=R;x=0;0;0;01;N=T/h;for t=l:N;kl=A*x+B*R;k2=A*(x+h*kl/3)+B*R;k3=A*(x+2*h*k2/3)+B*R;x=x+h*(kl+3*k3)/4;y(t)=C*x+D*R;end在命令行里键入A

26、= B二 C= D二R= T= h=y=hs (A, B, C, D, R, T, h) 得到结果。010(TO0010,B=000010-100-80.21-31.99_8.12-9.用题2-8仿真程序求解题2-7系统的闭环输出响应y (t).,C=-100 5 0 0,D=0解：A=1 0 00 1 00 0 1在命令行里键入 A二000-100 -80. 21 -31.99 B二0 0 0 1; C=-100 5 0 0; D二0; T=l; R=l; h=0.01; y=hs (A, B, C, D, R, T, h)数值解为08. 3333e-0075.8659e-0061.8115

27、e-0053. 9384e-0057.0346e-005%仅取一部分%12-10用式（2-34）梯形法求解试验方程y = y ,分析对计算步长h有何限制，说明h对数值解：编写梯形法程序为儿+1 =力+伙1+灯）f 1-9 亠 tb m.l ht Ith的选取不能超出上述范用，否则系统不稳左。得到儿+| =儿（1一+经）稳泄系统最终渐进收敛。 T 2 Z*计算得0 h2tc系统稳泄则1 2 +伫V1 271如图2-27所示斜梁滚球系统，若要研究滚球在梁上的位置可控性，需首先建立其数学模型， 已知力矩电机的输出转矩M与其电流i成正比，横梁为均匀可自平衡梁（即当电机不通电且无滚球 时，横梁可处于0=0的水平状态），是建立系统的数学模型，并给出简化后系统的动态结构图。解：设球的质心到杆的距离为0,该系统为特殊情况下的球棒系统。另令/rm./2分别表示棒的 惯量、球的质量和球的惯量。则球质心的位垃和速度为xe = （xcos 0.xsin0）ve = （ucos 6 一 AWsin 0. vsin 0 + aw cos 0）其中X = v, 6 = （0.因而动能的移动部分为因而动能的移动部分为 Klran, =2昉=*“（+*%）球棒系统的旋转动能为 心=g何+2（卡因而，系统总的动