《2022年《计算机仿真技术与CAD》习题答案 .pdf》由会员分享，可在线阅读，更多相关《2022年《计算机仿真技术与CAD》习题答案 .pdf（22页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第 0 章绪论0-1 什么是仿真？它所遵循的基本原则是什么？答：仿真是建立在控制理论、相似理论、信息处理技术和计算机技术等理论基础之上的，以计算机和其他专用物理效应设备为工具，利用系统模型对真实或假想的系统进行试验，并借助专家经验知识、统计数据和信息资料对试验结果进行分析和研究，进而做出决策的一门综合性的试验性科学。它所遵循的基本原则是相似原理。0-2 仿真的分类有几种？为什么？答：依据相似原理来分：物理仿真、数学仿真和混合仿真。物理仿真：就是应用几何相似原理，制作一个与实际系统相似但几何尺寸较小或较大的物理模型（例如飞机模型放在气流场相似的风洞中）进行实验研究。数学仿真：就是应用数学相似原理
2、，构成数学模型在计算机上进行研究。它由软硬件仿真环境、动画、图形显示、输出打印设备等组成。混合仿真又称数学物理仿真，它是为了提高仿真的可信度或者针对一些难以建模的实体，在系统研究中往往把数学仿真、物理仿真和实体结合起来组成一个复杂的仿真系统， 这种在仿真环节中有部分实物介入的混合仿真也称为半实物仿真或者半物理仿真。0-3 比较物理仿真和数学仿真的优缺点。答：在仿真研究中，数学仿真只要有一台数学仿真设备（如计算机等），就可以对不同的控制系统进行仿真实验和研究，而且，进行一次仿真实验研究的准备工作也比较简单，主要是受控系统的建模、控制方式的确立和计算机编程。数学仿真实验所需的时间比物理仿真大大缩短
3、，实验数据的处理也比物理仿真简单的多。与数学仿真相比，物理仿真总是有实物介入，效果直观逼真，精度高， 可信度高，具有实时性与在线性的特点；但其需要进行大量的设备制造、安装、 接线及调试工作，结构复杂，造价较高，耗时过长，灵活性差，改变参数困难，模型难以重用，通用性不强。0-4 简述计算机仿真的过程。答：第一步：根据仿真目的确定仿真方案根据仿真目的确定相应的仿真结构和方法，规定仿真的边界条件与约束条件。第二步：建立系统的数学模型对于简单的系统，可以通过某些基本定律来建立数学模型。而对于复杂的系统，则必须利用实验方法通过系统辩识技术来建立数学模型。数学模型是系统仿真的依据，所以，数学模型的准确性是
4、十分重要。精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 1 页，共 22 页 - - - - - - - - - - 第三步：建立仿真模型即通过一定算法对原系统的数学模型进行离散化处理，就连续系统言，就是建立相应的差分方程。第四步：编制仿真程序对于非实时仿真，可用一般高级语言或仿真语言。对于快速的实时仿真，往往需要用汇编语言。第五步：进行仿真实验并输出仿真结果通过实验对仿真系统模型及程序进行校验和修改，然后按系统仿真的要求输出仿真结果。0-5 什么是CAD 技术？控制系统CAD 可解决哪些问题？答：CA
5、D 技术，即计算机辅助设计（Computer Aided Design） ，是将计算机高速而精确的计算能力、大容量存储和数据处理能力与设计者的综合分析、逻辑判断以及创造性思维结合起来，以加快设计进程、缩短设计周期、提高设计质量的技术。控制系统CAD 可以解决以频域法为主要内容的经典控制理论和以时域法为主要内容的现代控制理论。此外，自适应控制、自校正控制以及最优控制等现代控制策略都可利用CAD 技术实现有效的分析和设计。精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 2 页，共 22 页 - - - - -
6、 - - - - - 第 1 章仿真软件 MATLAB 1-1 对于矩阵A=1 2;3 4， MATLAB以下四条命令：A.(0.5) ； A(0.5) ； sqrt(A) ； sqrtm(A) 所得结果相同吗？它们中哪个结果是复数矩阵，为什么？答：A.(0.5)=1.0000 1.4142;1.7321 2.0000 ；A(0.5)= 0.5537 + 0.4644i 0.8070 - 0.2124i;1.2104 - 0.3186i 1.7641 + 0.1458i；sqrt(A)= 1.0000 1.4142;1.7321 2.0000 ；sqrtm(A) = 0.5537 + 0.46
7、44i 0.8070 - 0.2124i;1.2104 - 0.3186i 1.7641 + 0.1458i；其中，“ A.(0.5) ”表示向量的乘方，“ A(0.5) ”表示矩阵的乘方，“ sqrt(A) ”只定义在矩阵的单个元素上，即分别对矩阵的每个元素进行运算，“ sqrtm(A) ”表示对矩阵（方阵）的超越函数进行运算。1-4 求二元函数方程组：sin(x-y)=0,cos(x+y)=0 的解。答：x,y=solve(sin(x-y)=0,cos(x+y)=0,x,y) x = -1/4*pi 1/4*pi y = -1/4*pi 1/4*pi 1-5 求函数y(t)=exp(-t)
8、*|sincost|的最大值（0=tf=(-1)*exp(-(abs(x)*abs(sin(cos(abs(x); x=fminsearch(f,0),ymax=exp(-(abs(x)*abs(sin(cos(abs(x) x = 0 ymax = 0.8415 1-6 设 D2y-3Dy+2y=x,y(0)=1,Dy(0)=0，求 y(0.5) 的值。答： f=D2y-3*Dy+2*y=x;g=dsolve(f,y(0)=1,Dy(0)=0,x);x=0.5;y=eval(g) y = 0.6100 1-7 求方程cos(t)2*exp(-0.1t)=0.5t的解。答：t1=solve(c
9、os(t)2*exp(-0.1*t)=0.5*t,t);t=eval(t1) 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 3 页，共 22 页 - - - - - - - - - - t = 0.8329 1-8 求方程组：x2+y2=1,xy=2 的解。答：x,y=solve(x2+y2=1,x*y=2,x,y) x = -1/2*(1/2*5(1/2)+1/2*i*3(1/2)3+1/4*5(1/2)+1/4*i*3(1/2) -1/2*(1/2*5(1/2)-1/2*i*3(1/2)3+1/4*
10、5(1/2)-1/4*i*3(1/2) -1/2*(-1/2*5(1/2)+1/2*i*3(1/2)3-1/4*5(1/2)+1/4*i*3(1/2) -1/2*(-1/2*5(1/2)-1/2*i*3(1/2)3-1/4*5(1/2)-1/4*i*3(1/2) y = 1/2*5(1/2)+1/2*i*3(1/2) 1/2*5(1/2)-1/2*i*3(1/2) -1/2*5(1/2)+1/2*i*3(1/2) -1/2*5(1/2)-1/2*i*3(1/2) 1-9 求 f(kT)=kexp(-akT)的 Z 变换表达式。答：syms k t z;f=k*exp(-a*t);F=ztran
11、s(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)/a 1-12 求以下方程组边值问题的解。Df=3f+4g, Dg=-4f+3g, f(0)=0, g(0)=1 答：f=Dx1=3*x1+4*x2,Dx2=-4*x1+3*x2;x1,x2=dsolve(f,x1(0)=0,x2(0)=1,t) x1 = exp(3*t)*sin(4*t) x2 = exp(3*t
12、)*cos(4*t) 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 4 页，共 22 页 - - - - - - - - - - 第 2 章 控制系统的数学模型及其转换2-1 已知系统的传递函数为23216116ssG ssss试用 MATLAB建立其状态空间表达式。答：num=1 1 1;den=1 6 11 6;A,B,C,D=tf2ss(num,den) A = -6 -11 -6 1 0 0 0 1 0 B = 1 0 0 C = 1 1 1 D = 0 2-2 已知系统的状态空间表达式为111
13、22212 0 11 0-2 -31 11 01 1xtxtutxtutxtxty txt试用 MATLAB求其传递函数阵。答： A=0 1;-2 -3;B=1 0;1 1;C=1 0;1 1;D=zeros(2,2); num1,den1=ss2tf(A,B,C,D,1),num2,den2=ss2tf(A,B,C,D,2) num1 = 0 1.0000 4.0000 0 2.0000 2.0000 den1 = 1 3 2 num2 = 0 0.0000 1.0000 0 1.0000 1.0000 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名
14、师归纳 - - - - - - - - - -第 5 页，共 22 页 - - - - - - - - - - den2 = 1 3 2 2-3 已知两子系统的传递函数分别为1112Gsss，213Gss s试利用MATLAB求两子系统串联和并联时系统的传递函数。答： num1=1;den1=1 3 2;num2=1;den2=1 3 0; num,den=series(num1,den1,num2,den2) num = 0 0 0 0 1 den = 1 6 11 6 0 num1=1;den1=1 3 2;num2=1;den2=1 3 0; num,den=parallel(num1,
15、den1,num2,den2) num = 0 0 2 6 2 den = 1 6 11 6 0 2-4 设系统的状态空间表达式为112212 0 11-2 -323 0 xtxtu txtxtxty txt若取线性变换阵1 11 -1P设新的状态变量为1xPx，则利用 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 0 3 -1 B1 = 3 -1 C1 = 1.5000 1.5000 精品资料 - - - 欢迎下载 - - - -
16、 - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 6 页，共 22 页 - - - - - - - - - - D1 = 0 2-5 已知离散系统状态空间表达式11221210 101 3111 1xkxku txkxkxky kxk试用 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.9598 7.3357 7.3357 24.9669 B1 = 1.9598 7.3357 C1 = 1 1 D1 = 0 精品资料 - - - 欢
17、迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 7 页，共 22 页 - - - - - - - - - - 第 3 章 连续系统的数字仿真3-1 已知线性定常系统的状态空间表达式为112212 0 12-5 -601 2xtxtu txtxtxty txt且初始状态为零，试利用四阶-龙格库塔法求系统的单位阶跃响应。答：%ex3_1.m r=1; A=0 1;-5 -6; B=2;0; C=1 2; d=0; Tf=5; h=0.1; x=zeros(length(A),1); y=0; t=0; for i=1:Tf/h K1
18、=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; end plot(t,y) 3-2 设单位反馈系统的开环传递函数42G ss s试利用二阶-龙格库塔法求系统的单位阶跃响应。精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 8 页，共 22 页 - - - - - - - - - - 答：%ex3_2.m r=1; nu
19、mo=4; deno=1,2,0; num,den=cloop(numo,deno); A,b,C,d=tf2ss(num,den);Tf=5; h=0.1; x=zeros(length(A),1); y=0; t=0; for i=1:Tf/h K1=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 利用 input( ) 函数修改例3-1 所给程序ex3_1.m ， 将其中给定的参数r， numo ， deno，numh 和 denh 利用键盘输入，使其变为连续控制系统
20、面向传递函数的通用数字仿真程序。答：精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 9 页，共 22 页 - - - - - - - - - - 3-5 利用 input( ) 函数修改例3-2 所给程序ex3_2.m ，将其中给定的参数r，P，W，W0和 Wc 利用键盘输入，使其变为连续控制系统面向结构图的通用数字仿真程序。答：精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 10 页，共 22 页 - - - -
21、- - - - - - 第 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;T=0.02;A=P(:,1);B=P(:,2);C=P(:,3);D=P(:,4);FZ=P(:,5)
22、;S=P(:,6);n=length(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);else E(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)=1;L(i)=D(i)/B(i);Q(i)=0;endelsei
23、f (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;else disp(A(i)=B(i)=0);endendendx=zeros(length(A),1);x0=x;z=x;u=zeros(length(A),1);u0=u;y=zeros(length(Wc(:,1),1);t=0;for j=1:Tf/T u1=u; u=W*x+W0*R;for i=1:nif (FZ(i)=0)if (FZ(i)=1) u(i)=saturation(u(i),S(i);endif (FZ(i)=
24、2) u(i)=deadzone(u(i),S(i);end精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 11 页，共 22 页 - - - - - - - - - - if (FZ(i)=3) u(i),u0(i)=backlash(u0(i),u(i),u1(i),S(i); endif (FZ(i)=4) u(i)=sign1(u(i),S(i);endendendx1=x;for i=1:n z(i)=E(i)*z(i)+F(i)*u(i)+G(i)*u1(i);x(i)=H(i)*z(i)+L
25、(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);endif (FZ(i)=6) x(i)=deadzone(x(i),S(i);endif (FZ(i)=7) x(i),x0(i)=backlash(x0(i),x(i),x1(i),S(i); endif (FZ(i)=8) x(i)=sign1(x(i),S(i);end endendy=y,Wc*x;t=t,t(j)+T;endplot(t,y)%saturation.m %子程序function x=saturatio
26、n(u,s)if (abs(u)=s)if (u0) x= s;else x=-s;endelsex= 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_1 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 12 页，共 22 页 - - - - - - - - - - 4-2 针对例3-2 所给线性定常系统，试利用第4 章所给程序，求系统的单位阶跃响应，并对其结果进行比较。答：ex3_2 ex4_1 精品资料 -
27、 - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 13 页，共 22 页 - - - - - - - - - - 4-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_1 4-4 利用 input( )函数修改例4-1 所给程序ex4_1.m ，将其中给定的参数R， P，W，W0 和 Wc 利用键盘输入，使其变为连续控制系
28、统按环节离散化的通用数字仿真程序。答：略精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 14 页，共 22 页 - - - - - - - - - - 第 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=
29、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);n=length(A);n1=length(Fr);m1=length(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);else E(i)=exp(-A(i)*T/ B(i); F(i)=(D(i)/B(i)- C(i)/ A(i)*(1- E(i)* B(i)/( A(i)
30、*T)-1); G(i)=(D(i)/B(i)- C(i)/ A(i)*(1+( E(i)-1)*(1+ B(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;else disp(A(i)= B(i)=0);endendendx=zeros(length(A),1); x0=x;z=x;u=zeros(length(A),1); u0=u;y=zeros(length(Wc(:,1
31、),1);t=0;Ur=zeros(n1,1); Er=zeros(m1,1);for ij=0:Tf/Tm;e=R-x(n);Er=e;Er(1:m1-1);ur=-Fr*Ur+ Gr*Er;Ur= ur;Ur(1:n1-1);for j=1:Tm/T u1= u; u = W*x+W0*ur;精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 15 页，共 22 页 - - - - - - - - - - for i=1:nif (FZ(i)=0)if (FZ(i)=1) u(i)=saturation(
32、u(i), S(i);endif (FZ(i)=2) u(i)=deadzone(u(i), S(i);endif(FZ(i)=3) u(i), u0(i)=backlash(u0(i), u(i), u1(i), S(i);endif (FZ(i)=4) u(i)=sign1(u(i), S(i);endendend x1= x;for i=1:n z(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)=satura
33、tion(x(i),S(i);endif(FZ(i)=6) x(i)=deadzone(x(i),S(i);endif (FZ(i)=7) x(i),x0(i)=backlash(x0(i),x(i),x1(i),S(i);endif(FZ(i)=8) x(i)=sign1(x(i),S(i);endendendy=y,Wc*x; t= t,t(length(t)+T;endendplot(t,y)ex5_1 5-2 针对例3-2 和例 4-1 所给连续系统，试利用第5 章所给程序，求系统的单位阶跃响应，并对其结果进行比较分析。答：ex3_2 精品资料 - - - 欢迎下载 - - - - -
34、 - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 16 页，共 22 页 - - - - - - - - - - ex4_1 ex5_2 5-4 略精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 17 页，共 22 页 - - - - - - - - - - 第 6 章 动态仿真集成环境Simulink 6-1 已知单变量系统如图题6-1 所示（图略） ，试利用Simulink求输出量y 的动态响应。答：6-2 假设某一系由图题6-2 所示的四个典型环节组成（图略），
35、试利用Simulink求输出量 y 的动态响应。答：6-3 已知非线性系统如图题6-3 所示，试利用Simulink求输出量y 的动态响应。精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 18 页，共 22 页 - - - - - - - - - - 答：6-4 已知采样系统结构如图题6-4 所示，试利用Simulink求输出量y 的动态响应。答：6-5 已知非线性系统如图6-5 所示，试利用Simulink分析非线性环节的c 值与输入幅值对系统输出性能的影响。答：精品资料 - - - 欢迎下载 - -
36、- - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 19 页，共 22 页 - - - - - - - - - - （1）r=1 c=0 （2）r=0.5 c=0 （3）r=1 c=1 （4）r=1 c=2 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 20 页，共 22 页 - - - - - - - - - - 6-6 已知线性定常系统的状态方程为112212 0 10-2 -310 1=-10 xtxtu txtxtxx试利用Simulink求 u(t
37、) 为单位阶跃函数时系统状态方程的解。答：t,x,y=sim(ex6_6,10); plot(t,y(:,1),:b,t,y(:,2),-r);legend(y1,y2) 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 21 页，共 22 页 - - - - - - - - - - 文档编码：KDHSIBDSUFVBSUDHSIDHSIBF-SDSD587FCDCVDCJUH 欢迎下载 精美文档欢迎下载 精品资料 - - - 欢迎下载 - - - - - - - - - - - 欢迎下载 名师归纳 - - - - - - - - - -第 22 页，共 22 页 - - - - - - - - - -