《计算机的仿真技术与CAD》习地训练题目答案详解

1、第0章绪论0-1什么是仿真？它所遵循的基本原则是什么？答：仿真是建立在控制理论、相似理论、信息处理技术和计算机技术等理论基础之上的，以计算机和其他专用物理效应设备为工具，利用系统模型对真实或假想的系统进 行试验，并借助专家经验知识、统计数据和信息资料对试验结果进行分析和研究，进 而做出决策的一门综合性的试验性科学。它所遵循的基本原则是相似原理。0-2仿真的分类有几种？为什么？答：依据相似原理来分：物理仿真、数学仿真和混合仿真。物理仿真：就是应用几何相似原理，制作一个与实际系统相似但几何尺寸较小或较大的物理模型（例如飞机模型放在气流场相似的风洞中）进行实验研究。数学仿真：就是应用数学相似原理，构

2、成数学模型在计算机上进行研究。它由软 硬件仿真环境、动画、图形显示、输出打印设备等组成。混合仿真又称数学物理仿真，它是为了提高仿真的可信度或者针对一些难以建模的实体，在系统研究中往往把数学仿真、物理仿真和实体结合起来组成一个复杂的仿 真系统，这种在仿真环节中有部分实物介入的混合仿真也称为半实物仿真或者半物理 仿真。0-3比较物理仿真和数学仿真的优缺点。答：在仿真研究中，数学仿真只要有一台数学仿真设备（如计算机等），就可以对不同的控制系统进行仿真实验和研究，而且，进行一次仿真实验研究的准备工作也比较 简单，主要是受控系统的建模、控制方式的确立和计算机编程。数学仿真实验所需的 时间比物理仿真大大缩

3、短，实验数据的处理也比物理仿真简单的多。与数学仿真相比， 物理仿真总是有实物介入，效果直观逼真， 精度高，可信度高,具有实时性与在线性的特点；但其需要进行大量的设备制造、安装、接线及调试工作，结构复杂，造价较高，耗时过长，灵活性差，改变参数困难，模型难以重用，通用性 不强。0-4简述计算机仿真的过程。答：第一步：根据仿真目的确定仿真方案根据仿真目的确定相应的仿真结构和方法，规定仿真的边界条件与约束条件。第二步：建立系统的数学模型对于简单的系统，可以通过某些基本定律来建立数学模型。而对于复杂的系统， 则必须利用实验方法通过系统辩识技术来建立数学模型。数学模型是系统仿真的依 据，所以，数学模型的准

4、确性是十分重要。第三步：建立仿真模型即通过一定算法对原系统的数学模型进行离散化处理，就连续系统言，就是建立相应的差分方程。第四步：编制仿真程序对于非实时仿真，可用一般高级语言或仿真语言。对于快速的实时仿真，往往需要用汇编语言。第五步：进行仿真实验并输出仿真结果通过实验对仿真系统模型及程序进行校验和修改，然后按系统仿真的要求输出仿真结果。0-5什么是CAD技术？控制系统 CAD可解决哪些问题？答：CAD技术，即计算机辅助设计(Computer Aided Design ),是将计算机高速而精确的计算能力、大容量存储和数据处理能力与设计者的综合分析、逻辑判断以及创造 性思维结合起来，以加快设计进程

5、、缩短设计周期、提高设计质量的技术。控制系统CAD可以解决以频域法为主要内容的经典控制理论和以时域法为主要 内容的现代控制理论。此外，自适应控制、自校正控制以及最优控制等现代控制策略 都可利用 CAD技术实现有效的分析和设计。第1章仿真软件一一MATLAB1-1对于矩阵 A=1 2;3 4 , MATLAB以下四条命令：A.a(0.5) ； AA(0.5) ； sqrt(A) ； sqrtm(A)所得结果相同吗？它们中哪个结果是复数矩阵，为什么？答：A.a(0.5)=1.00001.4142;1.73212.0000；Aa(0.5)= 0.5537 + 0.4644i 0.8070 - 0.2

6、124i;1.2104 - 0.3186i1.7641 + 0.1458i；sqrt(A)= 1.00001.4142;1.73212.0000；sqrtm(A) = 0.5537 + 0.4644i 0.8070 - 0.2124i;1.2104 - 0.3186i1.7641 + 0.1458i；其中，“ A.a(0.5) ”表示向量的乘方， “ AA(0.5) ”表示矩阵的乘方， “ sqrt(A) ”只定 义在矩阵的单个元素上，即分别对矩阵的每个元素进行运算，“ sqrtm(A) ”表示对矩阵(方阵)的超越函数进行运算。1-4 求二元函数方程组： sin(x-y)=0,cos(x+y)

7、=0 的解。答：x,y=solve(si n(x-y)=0,cos(x+y)=0,x,y)x =-1/4*pi1/4*piy =-1/4*pi1/4*pi1-5 求函数 y(t)=exp(-t)*|sincost|的最大值(0=tf=(-1)*exp(-(abs(x)*abs(si n(cos(abs(x);x=fm in search(f,0),ymax=exp(-(abs(x)*abs(s in( cos(abs(x)x =0ymax =0.84151-6 设 D2y-3Dy+2y=x,y(0)=1,Dy(0)=0 ，求 y(0.5)的值。答： f=D2y-3*Dy+2*y=x;g=dso

8、lve(f,y(0)=1,Dy(0)=0,x);x=0.5;y=eval(g)y =0.61001-7 求方程 cos(t)A2*exp(-0.1t)=0.5t 的解。答：t1=solve(cos(t)A2*exp(-0.1*t)=0.5*t,t);t=eval(t1)t =0.83291-8 求方程组：xA2+yA2=1,xy=2的解。答：x,y=solve(xA2+yA2=1,x*y=2,x,y)x =-1/2*(1/2*5A(1/2)+1/2*i*3A(1/2)A3+1/4*5A(1/2)+1/4*i*3A(1/2) -1/2*(1/2*5人(1/2)-1/2*产3人(1/2)人3+1/

9、4*5人(1/2)-1/4*产3人(1/2)-1/2*(-1/2*5A(1/2)+1/2*i*3A(1/2)A3-1/4*5A(1/2)+1/4*i*3A(1/2)-1/2*(-1/2*5A(1/2)-1/2*i*3A(1/2)A3-1/4*5A(1/2)-1/4*i*3A(1/2)y =1/2*5A(1/2)+1/2*i*3A(1/2)1/2*5A(1/2)-1/2*i*3A(1/2)-1/2*5A(1/2)+1/2*i*3A(1/2)-1/2*5A(1/2)-1/2*i*3A(1/2)1-9 求f(kT)=kexp(-akT) 的Z变换表达式。答：syms k t z;f=k*exp(-a

10、*t);F=ztra ns(f,t,z)f =k*z/exp(-a)/(z/exp(-a)-1)1-10 求一阶微分方程Dx=ax+by(t),x(0)=x0 的解。答：f=Dx=a*x+b*y;x=dsolve(f,x(0)=x0,t)x =-b*y/a+exp(a*t)*(b*y+x0*a)/a1- 12求以下方程组边值问题的解。Df=3f+4g, Dg=-4f+3g, f(0)=0, g(0)=1答：f=Dx 仁3*x1+4*x2,Dx2=-4*x1+3*x2;x1,x2=dsolve(f,x1(0)=0,x2(0)=1,t)x1 =exp(3*t)*si n(4*t)x2 =exp(3

11、*t)*cos(4*t)第2章控制系统的数学模型及其转换2- 1已知系统的传递函数为g /s2 +s+1G(s 尸 S3 +6s2 +11s 衣试用MATLAB 建立其状态空间表达式。答：n um=1 1 1;de n=1 6 11 6;A,B,C,D=tf2ss( num,de n)A =-6-11-6100010B =100C =1 1 1D =02- 2已知系统的状态空间表达式为X1X2一 o 1X1勺1 3t t1 2 :一 o 1x1 (t JL F 吐2(t)yt =1试用MATLAB求其传递函数阵。答： A=0 1;-2 -3;B=1 0;1 1;C=1 0;1 1;D=zero

12、s(2,2); nu m1,de n1=ss2tf(A,B,C,D,1), nu m2,de n2=ss2tf(A,B,C,D,2)num1 =01.00004.000002.00002.0000den1 =132num2 =00.00001.000001.00001.0000den2 =2-3已知两子系统的传递函数分别为1 ,r1)_(s +1 (s +2)，s _s s 3试利用MATLAB 求两子系统串联和并联时系统的传递函数。答: num1=1;de n1=1 3 2; num2=1;de n2=1 3 0;nu m,de n=series( nu m1,de n1,nu m2,de

13、n2)num =00001den =161160 num1=1;de n1=1 3 2; num2=1;de n2=1 3 0;nu m,de n=parallel( nu m1,de n1,nu m2,de n2)num =den =112-4设系统的状态空间表达式为rii粗）y t - 3P=1=0若取线性变换阵-J设新的状态变量为xpx，则利用MATLAB求在新状态变量下，系统状态空间表达式。答:A=0 1;-2 -3;B=1;2;C=3 0;D=0;P=1 1;1 -1;A1,B1,C1,D1=ss2ss(A,B,C,D,P)A1 =-2 03-1B1 =3-1C1 =1.50001.

14、5000D1 =02-5已知离散系统状态空间表达式( k X1X2- yfl-LJrX1X2- .1 3 k k1 1 k k试用MATLAB求其系统的脉冲传递函数。答：A=0 1;1 3;B=0;1;C=1 1;D=0;T=1;A1,B1,C1,D1=c2dm(A,B,C,D,T)A1 =2.95987.33577.335724.9669B1 =1.95987.3357C1 =1 1D1 =0第3章连续系统的数字仿真3- 1已知线性定常系统的状态空间表达式为且初始状态为零，试利用四阶-龙格库塔法求系统的单位阶跃响应。答：%ex3_1.mr=1; A=0 1;-5 -6; B=2;0; C=1

15、 2; d=0;Tf=5; h=0.1;x=zeros(le ngth(A),1); y=0; t=0;for i=1:Tf/hK1=A* x+B*r;K2=A*(x+h*K1/2)+B*r;K3=A*(x+h*K2/2)+B*r;K4=A*(x+h*K3)+B*r;x=x+h*(K1+2*K2+2*K3+K4)/6;y=y;C*x; t=t;t(i)+h;endplot(t,y)0.2 -0.44)6OB-1 -00.511.522.533.5 J J 53- 2设单位反馈系统的开环传递函数试利用二阶-龙格库塔法求系统的单位阶跃响应。答:%ex3_2.mr=1; numo=4; deno=1

16、,2,0; nu m,de n=cloop (num o,de no);A,b,C,d=tf2ss( nu m,den);Tf=5; h=0.1;x=zeros(le ngth(A),1); y=0; t=0;for i=1:Tf/hK仁A* x+b*r;K2=A*(x+h*K1)+b*r; x=x+h*(K1+K2)/2; y=y;C*x; t=t;t(i)+h;end plot(t,y)3-4利用in put()函数修改例 3-1所给程序 ex3_1.m，将其中给定的参数 r, n umo , de no ,numh和denh利用键盘输入，使其变为连续控制系统面向传递函数的通用数字仿真程序

17、。答:Sexr=input 输入给定信号匸二);numinput输入开环传递函数分手多项式系数向量mm2 ); deno= input (辅入幵环传递函数分母多项式系数向sdeno=);Entnii, den=cla(jp (numo, deno) ; A$d =tf 2ss (num, d&n);Tf 二 input仿算时阖);h= input (计算步长h=、k= seros(length(A)3 1) ; y=0k t=0;for i=l:Tf/hK1=A* x+b*r ;K2=A*(s+h#Kl/2)+b*r :K3=A* (x+h*K2/2)+b*r .K4=A+(K+h*K3)+b

18、*r ；s=x+h* (K1+2+K2+2+K3+K4)/6;y= y;C#K ;t=t ;t(iHh; endplot (t,y)3-5利用input()函数修改例 3-2所给程序ex3_2.m，将其中给定的参数r, P, W, WO和Wc利用键盘输入，使其变为连续控制系统面向结构图的通用数字仿真程序。答：%es3_5, hlr= input (输入给定信号;P= input辅入典型环节系数矩阵2 ).Winput谕入各典型环节谕入输出间的连接关系矩阵归J ;W0=input (输入系统的参考输入与各环节输入闾的连接关系矩阵喲三);Wcinput输入系统的输岀与各环节输出间的连接关系矩阵WW

19、 );Tf-input C 仿JI时间Tf=) : h=input C 计算步长h二* );Al=diag(P(:f D) ; Bl=diag (P 2) : Cl=dia(P(:f 3); Dl-diag:H=B1-DUV; Q=C 1 *W-A 1; A=inv(H) *Q ; B=mv *C1*WO.x= zeros (length(A), 1) ; y= zeros (lengthCWc(:? l)j 1);t=0;for i= 1.: Tf/hKl-A*x+B*r: K2=A*(x+h*Kl/2)+B*r,K3A+(i：+h*K2/2)+B+r : K4=A# (z+h*K3)+0*

20、r;x=x+h* (K1+2*E2+2*K3+K4)/S.y= y, Wc*k ; t= t, t (i)-Hi:endplot Ct,y)第4章连续系统按环节离散化的数字仿真4- 1已知非线性习题如图题4-1所示，试利用连续系统按环节离散化的数字仿真方法,求输出量y的动态响应，并与无非线性环节进行比较。(图略)答：%ex4_1.m% 主程序R=10;P=0.1 1 0.5 1 5 5; 0 1 1 0 0 0; 2 1 2 0 0 0;10 1 10 0 0 0;W=0 0 0 -1; 1 0 0 0; 0 1 0 0; 0 0 1 0;W0=1;0;0;0;Wc=0 0 0 1;Tf=25

21、;T=0.02;A=P(:,1);B=P(:,2);C=P(:,3);D=P(:,4);FZ=P(:,5);S=P(:,6);n=len gth(A);for i=1: nif (A(i)=0)if (B(i)=0)E(i)=0;F(i)=0;G(i)=0;H(i)=0;L(i)=(C(i)+D(i)/T)/A(i);Q(i)=-D(i)/(A(i)*T);elseE(i)=exp(-A(i)*T/B(i);F(i)=(D(i)/B(i)-C(i)/A(i)*(1-E(i)*B(i)/(A(i)*T)-1);G(i)=(D(i)/B(i)-C(i)/A(i)*(1+(E(i)-1)*(1+B(

22、i)/(A(i)*T);H(i)=1;L(i)=D(i)/B(i);Q(i)=0;endelseif (B(i)=0)E(i)=1;F(i)=0.5*C(i)*T/B(i);G(i)=F(i);H(i)=1;L(i)=D(i)/B(i);Q(i)=0;elsedisp( A(i)=B(i)=0);endendendx=zeros(le ngth(A),1);x0=x;z=x;u=zeros(le ngth(A),1);u0=u;y=zeros(le ngth(Wc(:,1),1);t=0;for j=1:Tf/Tu1=u; u=W*x+W0*R;for i=1: nif (FZ(i)=0)if

23、 (FZ(i)=1) u(i)=saturatio n(u(i),S(i);endif (FZ(i)=2) u(i)=deadz on e(u(i),S(i);endendif (FZ(i)=3) u(i),uO(i)=backlash(uO(i),u(i),u1(i),S(i);if (FZ(i)=4) u(i)=sig n1(u(i),S(i);endendendx1=x;for i=1: nz(i)=E(i)*z(i)+F(i)*u(i)+G(i)*u1(i);x(i)=H(i)*z(i)+L(i)*u(i)+Q(i)*u1(i);endfor i=1: nif (FZ(i)=0)if

24、(FZ(i)=5) x(i)=saturatio n(x(i),S(i);endendif (FZ(i)=6) x(i)=deadzo ne(x(i),S(i);endif (FZ(i)=7) x(i),xO(i)=backlash(xO(i),x(i),x1(i),S(i);if (FZ(i)=8) x(i)=sig n1(x(i),S(i);endendendy=y,Wc*x;t=t,t(j)+T;endplot(t,y)%saturatio n.m%子程序fun cti onx=saturati on( u,s)if (abs(u)=s)if (u0) x= s;else x=-s;en

25、delsex= u;end修改“ P=0.1 1 0.5 1 0 0; 0 1 1 0 0 0; 2 1 2 0 0 0;10 1 10 0 0 0;ex4_14- 2针对例3-2所给线性定常系统，试利用第4章所给程序，求系统的单位阶跃响应,并对其结果进行比较。答：ex3_2ex4_14-3针对例4-1所给系统，去掉饱和非线性环节后求系统的单位阶跃响应，并与例4-1所得结果进行比较。答：ex4_1修改“ P=0.1 1 0.5 1 0 0; 0 1 1 0 0 0; 2 1 2 0 0 0;10 1 10 0 0 0;ex4_14-4利用input()函数修改例 4-1所给程序 ex4_1.m

26、，将其中给定的参数R，P，W，W0和Wc利用键盘输入，使其变为连续控制系统按环节离散化的通用数字仿真程序。答：略第5章采样控制系统的数字仿真5- 1已知采样控制系统的结构图如图题5-1所示(图略)。试利用采样控制系统的数字仿真方法，求当采样周期T=0.1s，且初始状态为零时，离散系统的单位阶跃响应。答：%ex5_1.mR=1;Gr=1;Fr=0;P=1 1 1 0 0 0;1 2 1 0 0 0;W=0 0;1 0;W0=1;0; Wc=0 1;Tf=25;Tm=0.1;T=0.01;A=P(:,1);B=P(:,2);C=P(:,3);D=P(:,4);FZ=P(:,5);S=P(:,6);

27、n=le ngth(A); n1=le ngth(Fr);m1=le ngth(Gr);for i=1: nif (A(i)=0)if (B(i)=0)E(i)=0;F(i)=0;G(i)=0;H(i)=0;L(i)=(C(i)+D(i)/T)/ A(i); Q(i)=- D(i)/( A(i)*T);elseE(i)=exp(-A(i)*T/ B(i);F(i)=(D(i)/B(i)- C(i)/ A(i)*(1- E(i)* B(i)/( A(i)*T)-1);G(i)=(D(i)/B(i)- C(i)/ A(i)*(1+( E(i)-1)*(1+ B(i)/( A(i)*T);H(i)=

28、1; L(i)=D(i)/ B(i); Q(i)=0;endelseif (B(i)=0)E(i)=1;F(i)=0.5*C(i)*T/B(i);G(i)=F(i);H(i)=1;L(i)=D(i)/B(i);Q(i)=0;elsedisp( A(i)= B(i)=0);endendendx=zeros(le ngth(A),1); x0=x;z=x;u=zeros(length(A),1); u0=u;y=zeros(le ngth(Wc(:,1),1);t=0;Ur=zeros(n1,1); Er=zeros(m1,1);for ij=0:Tf/Tm;e=R-x( n);Er=e;Er(1

29、:m1-1);ur=-Fr* Ur+ Gr*Er;Ur= ur;Ur(1: n1-1);for j=1:Tm/Tu1= u; u = W*x+W0*ur;for i=1: nif (FZ(i)=0)if (FZ(i)=1) u(i)=saturatio n( u(i), S(i);endif (FZ(i)=2) u(i)=deadzo ne(u(i), S(i);endif (FZ(i)=3)u(i), u0(i)=backlash(u0(i),u(i),u1(i),S(i);if (FZ(i)=4) u(i)=sig n1(u(i), S(i);endendendx1 = x;for i=1: nz(i)=E(i)*z(i)+F(i)* u(i)+G(i)*u1(i);x(i)=H(i)*z(i)+L(i)* u(i)+Q(i)*u1(i);endfor i=1: nif (FZ(i)=0)if (FZ(i)=5) x(i)=saturation(x(i),S(i);end