《计算机仿真技术与CAD》习题答案

第0章绪论0-1 什么是仿真？它所遵循的基本原则是什么？答：仿真是建立在控制理论、相似理论、信息处理技术和计算机技术等理论基础之上的，以计算机和其他专用物理效应设备为工具，利用系统模型对真实或假想的系统进行试验，并借助专家经验知识、统计数据和信息资料对试验结果进行分析和研究，进而做出决策的一门综合性的试验性科学。

它所遵循的基本原则是相似原理。

0-2 仿真的分类有几种？为什么？答：依据相似原理来分：物理仿真、数学仿真和混合仿真。

物理仿真：就是应用几何相似原理，制作一个与实际系统相似但几何尺寸较小或较大的物理模型（例如飞机模型放在气流场相似的风洞中）进行实验研究。

数学仿真：就是应用数学相似原理，构成数学模型在计算机上进行研究。

它由软硬件仿真环境、动画、图形显示、输出打印设备等组成。

混合仿真又称数学物理仿真，它是为了提高仿真的可信度或者针对一些难以建模的实体，在系统研究中往往把数学仿真、物理仿真和实体结合起来组成一个复杂的仿真系统，这种在仿真环节中有部分实物介入的混合仿真也称为半实物仿真或者半物理仿真。

0-3 比较物理仿真和数学仿真的优缺点。

答：在仿真研究中，数学仿真只要有一台数学仿真设备（如计算机等），就可以对不同的控制系统进行仿真实验和研究，而且，进行一次仿真实验研究的准备工作也比较简单，主要是受控系统的建模、控制方式的确立和计算机编程。

数学仿真实验所需的时间比物理仿真大大缩短，实验数据的处理也比物理仿真简单的多。

与数学仿真相比，物理仿真总是有实物介入，效果直观逼真，精度高，可信度高，具有实时性与在线性的特点；但其需要进行大量的设备制造、安装、接线及调试工作，结构复杂，造价较高，耗时过长，灵活性差，改变参数困难，模型难以重用，通用性不强。

0-4 简述计算机仿真的过程。

答：第一步：根据仿真目的确定仿真方案根据仿真目的确定相应的仿真结构和方法，规定仿真的边界条件与约束条件。

第二步：建立系统的数学模型对于简单的系统，可以通过某些基本定律来建立数学模型。

而对于复杂的系统，则必须利用实验方法通过系统辩识技术来建立数学模型。

数学模型是系统仿真的依据，所以，数学模型的准确性是十分重要。

第三步：建立仿真模型即通过一定算法对原系统的数学模型进行离散化处理，就连续系统言，就是建立相应的差分方程。

第四步：编制仿真程序对于非实时仿真，可用一般高级语言或仿真语言。

对于快速的实时仿真，往往需要用汇编语言。

第五步：进行仿真实验并输出仿真结果通过实验对仿真系统模型及程序进行校验和修改，然后按系统仿真的要求输出仿真结果。

0-5 什么是CAD技术？控制系统CAD可解决哪些问题？答：CAD技术，即计算机辅助设计（Computer Aided Design），是将计算机高速而精确的计算能力、大容量存储和数据处理能力与设计者的综合分析、逻辑判断以及创造性思维结合起来，以加快设计进程、缩短设计周期、提高设计质量的技术。

控制系统CAD可以解决以频域法为主要内容的经典控制理论和以时域法为主要内容的现代控制理论。

此外，自适应控制、自校正控制以及最优控制等现代控制策略都可利用CAD技术实现有效的分析和设计。

第1章仿真软件——MATLAB1-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.4644i 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*pi1/4*piy =-1/4*pi1/4*pi1-5 求函数y(t)=exp(-t)*|sin[cost]|的最大值（0<=t<inf）。

答：>>f='(-1)*exp(-(abs(x)))*abs(sin(cos(abs(x))))';>>x=fminsearch(f,0),ymax=exp(-(abs(x)))*abs(sin(cos(abs(x))))x =ymax =0.84151-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.61001-7 求方程cos(t)^2*exp(-0.1t)=0.5t的解。

答：>>t1=solve('cos(t)^2*exp(-0.1*t)=0.5*t','t');t=eval(t1)t =0.83291-8 求方程组：x^2+y^2=1,xy=2 的解。

答：>>[x,y]=solve('x^2+y^2=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*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=ztrans(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='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)*cos(4*t)第2章 控制系统的数学模型及其转换2-1 已知系统的传递函数为()23216116s s G s s s s ++=+++ 试用MATLAB 建立其状态空间表达式。

答：>>num=[1 1 1];den=[1 6 11 6];[A,B,C,D]=tf2ss(num,den)A =-6 -11 -61 0 00 1 0 B =10 C = 1 1 1D =2-2 已知系统的状态空间表达式为()()()()()()()()()11122212 0 1 1 0-2 -3 1 11 01 1x t x t u t x t u t x t x t y t x t ⎧⎡⎤⎡⎤⎡⎤⎡⎤⎡⎤⎪⎢⎥=+⎢⎥⎢⎥⎢⎥⎢⎥⎪⎢⎥⎣⎦⎣⎦⎢⎥⎢⎥⎪⎣⎦⎣⎦⎢⎥⎣⎦⎨⎪⎡⎤⎡⎤⎪=⎢⎥⎢⎥⎪⎣⎦⎢⎥⎣⎦⎩g g 试用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.00000 2.0000 2.0000den1 =1 3 2num2 =0 0.0000 1.00000 1.0000 1.00001 3 22-3 已知两子系统的传递函数分别为()()()1112G s s s =++，()()213G s s 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 1den =1 6 11 6 0>> num1=1;den1=[1 3 2];num2=1;den2=[1 3 0];>>[num,den]=parallel(num1,den1,num2,den2)num =0 0 2 6 2den =1 6 11 6 02-4 设系统的状态空间表达式为()()()()()()[]()()112212 0 11-2 -323 0x t x t u t x t x t x t y t x t ⎧⎡⎤⎡⎤⎡⎤⎡⎤⎪⎢⎥=+⎡⎤⎢⎥⎢⎥⎢⎥⎣⎦⎪⎢⎥⎣⎦⎣⎦⎢⎥⎪⎣⎦⎢⎥⎣⎦⎨⎪⎡⎤⎪=⎢⎥⎪⎢⎥⎣⎦⎩g g 若取线性变换阵 1 11 -1P ⎡⎤=⎢⎥⎣⎦ 设新的状态变量为1x P x ∧-=，则利用MA TLAB 求在新状态变量下，系统状态空间表达式。