控制系统数字仿真与CAD习题全答案

《控制系统数字仿真》练习题及答案1. 单选题1. 某系统的函数关系式为y=1/（x3-2x+4），绘制x在0至10之间变化的图形，正确的是（）。

A. fplot('1/(x*3-2*x+4)',[0 10])B. fplot('1/(x.^3-2*x+4)',[0 10])C. plot('1/(x.^3-2*x+4)',[0 10])D. plot('1/(x*3-2*x+4)',[0 10])正确答案：B2. 绘制系统零极点图的命令是（）。

A. stepB. pzmapC. rlocusD. sgrid正确答案：B3. 将系统零极点形式变换为传递函数形式的命令是（）。

A. tf2zpB. ss2tfC. ss2zpD. zp2tf正确答案：D4. AUTOCAD的坐标体系，包括世界坐标和（）坐标系。

A. 绝对坐标B. 平面坐标C. 相对坐标D. 用户坐标正确答案：D5. 在MATLAB工作空间中，表示圆周率的特殊变量是（）。

A. piB. ansC. iD. eps正确答案：A6. 下列哪条指令是求矩阵的行列式的值（）。

A. invB. diagC. detD. eig正确答案：C7. 在CAD网络系统中，以下说法不正确的是（）。

A. 设计资料可以共享B. 硬件可以共享C. 电脑文件可以共享D. 可以方便管理设计进度正确答案：C8. i=2; a=2i;b=2*i;c=2*sqrt(-1);程序执行后；a, b, c的值分别是（）。

A. a=4,b=4,c=2.0000iB. a=4,b=2.0000i, c=2.0000iC. a=2.0000i, b=4,c=2.0000iD. a=2.0000i,b=2.0000i,c=2.0000i正确答案：C9. 在循环结构中跳出循环，执行循环后面代码的命令为（）。

A. returnB. breakC. continueD. keyboard正确答案：B10. figure命令的作用是（）。

球车系统问题的解答1. 建模问题若给定球车系统的参数为小车质量3kg M =；小球质量2kg m =；圆弧轨道半径长0.9m R =；小球半径长0.1m r =；重力加速度取210m/s g =，则可进一步简化/线性化模型为40.28100.2x F F ϕϕϕ=-+⎧⎨=-⎩ (1) 对其进行拉氏变换得到系统的传递函数模型12222()0.2()()10() 5.414()()s G s F s s X s s G s s s Φ-⎧==⎪-⎪⎨-+⎪==⎪Φ⎩(2) 由传递函数模型得到系统动态结构图为图6-8 球车系统动态结构图由此可以得到球车系统状态方程为001000001004000.28010000.2x x F x x ϕϕϕϕ⎡⎤⎡⎤⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥=+⎢⎥⎢⎥⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦⎣⎦⎣⎦ (3) 经过与一阶倒立摆的数学模型比较发现，如果将R r +看做倒立摆的半杆长，将()2[/]J R r r +看作是摆杆的转动惯量，则球车系统的数学模型与一阶倒立摆的数学模型完全相同。

因此，它们的传递函数的形式也完全相同。

这又一次证明了球车系统是一阶倒立摆系统的一种变形，所以对一阶倒立摆的控制方法是完全可以用在球车系统的控制上的。

2. 检测问题传统的位置测量方法，如光电编码器及旋转电位器等是无法应用在对本系统小球在弧形导轨上的位置测量上的。

所以必须考虑采用其它的测量方法。

一种合理的设计如下：小球是导电材料制成的，金属弧形导轨分为上下两层。

两层之间互相隔离，且它们都与小球相接触。

沿弧形导轨上层装有电阻线圈，在线圈的两端分别施加±5V 的电压。

输出端的负极在电阻线圈中央位置，输出端的正极为弧形导轨下层。

这样，小球与弧形导轨上下层共同构成了一个电位器，如下图所示。

当小球处于最左端位置时，输出端输出-5V电压，当小球处于最右端位置时，输出端输出+5V电压，输出电压的值与小球的位置呈线性关系。

3-1．求解下列线性方程，并进行解得验证：(1)7 2 1 -249 15 3 -27-2 -2 11 511 32 130x⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥-⎢⎥⎢⎥⎣⎦⎣⎦,(2)5 76 5 1247 10 8 7 2346 8 10 9 3365 7 9 10 4351 2 3 4 515x⎡⎤⎡⎢⎥⎢⎢⎥⎢⎢⎥⎢=⎢⎥⎢⎢⎥⎢⎢⎥⎢⎣⎦⎣9613614414060⎤⎥⎥⎥⎥⎥⎥⎦由A*X=B得：X=A\B解：>> a=[7 2 1 -2;9 15 3 -2;-2 -2 11 5;1 3 2 13]a =7 2 1 -29 15 3 -2-2 -2 11 51 32 13>> b=[4 7 -1 0]'b =47-1>> x=a\bx =0.49790.14450.0629-0.0813(2)解：>> a=[5 7 6 5 17 10 8 7 26 8 10 9 35 7 9 10 41 2 3 4 5]a =5 76 5 17 10 8 7 26 8 10 9 35 7 9 10 41 2 3 4 5>> b=[24 9634 136 36 144 35 140 15 60] b =24 96 34 136 36 144 35 140 15 60>> x=a\b x =1.0000 4.0000 1.0000 4.0000 1.0000 4.0000 1.0000 4.00001.0000 4.00003-2．进行下列计算，给出不使用for 和while 等循环语句的计算方法。

（1）632ii k ==∑解:根据等比数列求和方法，在利用matlab 中的m 文件，编写程序求解。

M 文件为 n=64;q=2;k=(1-q^n)/(1-q); disp('k 的值为'); disp(k);保存文件q1.m在matlab 命令框中输入 >> q1k 的值为 1.8447e+019(2)求出y=x*sin(x) 在0<x<100条件下的每个峰值解：画出图形 >> x=0:0.01:100;>> y=x.*sin(x); >> plot(x,y); >> grid on>> title('y=x*sin(x)') >> xlabel('x') >>ylabel('y')方法1。

课后答案网 您最真诚的朋友
网团队竭诚为学生服务，免费提供各门课后答案，不用积分，甚至不用注
册，旨在为广大学生提供自主学习的平台！
课后答案网：
视频教程网：
PPT 课件网：
课后答案网：
若侵犯了您的版权利益，敬请来信告知！
课后答案网 w w w .h a c k s h p .c n
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网
课后答案网 。

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

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

1-2在系统分析与设计中仿真法与解析法有何区别？各有什么特点？答：解析法就是运用已掌握的理论知识对控制系统进行理论上的分析，计算。

它是一种纯物理意义上的实验分析方法，在对系统的认识过程中具有普遍意义。

由于受到理论的不完善性以及对事物认识的不全面性等因素的影响，其应用往往有很大局限性。

仿真法基于相似原理，是在模型上所进行的系统性能分析与研究的实验方法。

1-3数字仿真包括那几个要素？其关系如何？答: 通常情况下，数字仿真实验包括三个基本要素，即实际系统，数学模型与计算机。

由图可见，将实际系统抽象为数学模型，称之为一次模型化，它还涉及到系统辨识技术问题，统称为建模问题；将数学模型转化为可在计算机上运行的仿真模型，称之为二次模型化，这涉及到仿真技术问题，统称为仿真实验。

1-4为什么说模拟仿真较数字仿真精度低？其优点如何？。

答：由于受到电路元件精度的制约和容易受到外界的干扰，模拟仿真较数字仿真精度低但模拟仿真具有如下优点：（1）描述连续的物理系统的动态过程比较自然和逼真。

（2）仿真速度极快，失真小，结果可信度高。

（3）能快速求解微分方程。

模拟计算机运行时各运算器是并行工作的，模拟机的解题速度与原系统的复杂程度无关。

（4）可以灵活设置仿真试验的时间标尺，既可以进行实时仿真，也可以进行非实时仿真。

（5）易于和实物相连。

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

控制系统数字仿真离线作业
1什么是仿真，它所遵循的原则是什么？
2在系统分析与设计中仿真法与解析法有何区别？各有什么特点？
3什么是CAD技术？控制系统CAD技术可以解决哪些问题？
4什么是虚拟现实技术？它与仿真技术的关系如何？
5为什么称计算机控制系统为采样控制系统或离散控制系统，其仿真有何特点通常把系统中的离散信号是脉冲序列形成的离散系统,称为采样控制系统或脉冲控制系统;而把数字序列形成的离散系统,称为采样控制系统或计算机控制系统.
采样系统中既有离散信号,又有连续信号. 采样开关接通时刻,系统处于闭环工作状态.而在采样开关断开时刻,系统处于开环工作状态.
特点：高精度,高可靠,有效抑制干扰,良好的通用性.
6用MATLAB语言求下列系统的状态方程、传递函数、零极点增益和部分分式的模型参数，并分别写出其相应的数学模型表达式。

7单位反馈系统的开环传递函数如下
用MATLAB语句、函数求取系统闭环零极点，并求取系统闭环状态方程的可控标准型实现。

8某反馈控制系统的开环传递传递函数为
试绘制其根轨迹。

控制系统数字仿真题库一、填空题1.定义一个系统时，首先要确定系统的边界；边界确定了系统的范围，边界以外对系统的作用称为系统的输入，系统对边界以为环境的作用称为系统的输出。

2.系统的三大要素为：实体、属性和活动。

3.人们描述系统的常见术语为：实体、属性、事件和活动。

4.人们经常把系统分成四类，它们分别为：连续系统、离散系统、采样数据系统和离散-连续系统。

5、根据系统的属性可以将系统分成两大类：工程系统和非工程系统。

6.根据描述方法不同，离散系统可以分为：离散时间系统和离散事件系统。

7.系统是指相互联系又相互作用的实体的有机组合。

8.根据模型的表达形式，模型可以分为物理模型和数学模型二大类，其中数学模型根据数学表达形式的不同可分为二种，分别为：静态模型和动态模型。

9、采用一定比例按照真实系统的样子制作的模型称为物理模型，用数学表达式来描述系统内在规律的模型称为数学模型。

10.静态模型的数学表达形式一般是代数方程和逻辑关系表达式等，而动态模型的数学表达形式一般是微分方程和差分方程。

11.系统模型根据描述变量的函数关系可以分类为线性模型和非线性模型。

12 仿真模型的校核是指检验数字仿真模型和数学模型是否一致。

13．仿真模型的验证是指检验数字仿真模型和实际系统是否一致。

14．计算机仿真的三个要素为：系统、模型与计算机。

15．系统仿真的三个基本活动是系统建模、仿真建模和仿真试验。

16．系统仿真根据模型种类的不同可分为：物理仿真、数学仿真和数学-物理混合仿真。

17．根据仿真应用目的的不同，人们经常把计算机仿真应用分为四类，分别为：系统分析、系统设计、理论验证和人员训练。

18.计算机仿真是指将模型在计算机上进行实验的过程。

19.仿真依据的基本原则是:相似原理。

20.连续系统仿真中常见的一对矛盾为计算速度和计算精度。

21．保持器是一种将离散时间信号恢复成连续信号的装置。

22．零阶保持器能较好地再现阶跃信号。

23.一阶保持器能较好地再现斜坡信号。

5-1 设控制系统的开环传递函数为2(1)()()(1)(416)K s G s H s s s s s +=-++试画出该系统的根轨迹。

解: 在Matlab 窗口中输入下列命令： num=[1 1]; a=[1 0]; b=[1 -1]; c=[1 4 16]; d=conv(a,b); den=conv(d,c); rlocus(num,den) grid on可得到系统的根轨迹如下图所示：5-2 某反馈控制系统的开环传递函数为()()()()24420KG s H s s s s s =+++ 试绘制其根轨迹。

解：在MATLAB 命令窗口中输入下列命令： num=1;den=conv(conv([1,0],[1,4]),[1,4,20]); rlocus(num,den) grid on 运行结果为：5-3.已知某系统传递函数为2180(1)100()11(1)[()20.31]40200200s W s s s s +=++⨯⨯+ 试绘制其伯德图。

解：分子分母同乘100*200得到280200(100)()(2.5100)(20.3200)200s W s s s s ⨯+=++⨯+在Matlab 窗口中输入下列命令：k=80*200; num=[1 100]; a=[2.5 100];b=[(1/200) 2*0.3 200]; 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)); grid;xlabel('Frequency(rad/sec)'); ylabel('Gain(dB)'); subplot(2,1,2); semilogx(w,p); grid;xlabel('Frequency(rad/s)'); ylabel('Phase(deg)');可绘制该系统的伯德图如下所示。

第七节电力电子器件建模一、问题的提出上一节“电力电子系统建模”中所涉及到的电力电子器件(GTO、MOSFET、IGBT)都是理想开关模型（“0”、“1”状态），如表1。

然而，当我们在研究微观时间尺度下的（电压电流）系统响应或者电力电子器件特性的时候，我们就必须对电力电子器件建立更精确的模型。

这里的电力电子器件模型将不再是状态空间表达式或者传递函数的形式，这是因为简单形式的状态空间表达式或者传递函数已经无法精确表达出器件的动、静态过程。

电力电子器件的精确模型主要应用在：器件模型换向过程（微观时间尺度上）、元器件张力、功率消耗、设计器件缓冲电路等情况下。

从某种意义上说电力电子器件建模是电力电子系统建模的补充。

表1 理想开关与实际功率开关对比二、建模机理1.电力电子器件建模需考虑的问题对于功率半导体器件模型的发展，除了考虑半导体器件在建模时所考虑的一般问题和因素之外，在建立比较精确的仿真模型时，以下几个问题必须优先考虑，这些问题在低功率器件中不成问题，但在功率电子器件中这几个问题它们支配了器件的静态和动态特性：(1). 阻系数的调制为了承受较高的电压，功率半导体器件一般都有一个稍微厚度搀杂半导体层，当器件导通时，这个层决定导通压降和功率损失。

这个电阻随电压和电流变化而变化，具有非线性电阻的特性。

单极型器件（MOSFET）中，电阻的变化是由有效电流导通区域变化所引起，另外随着外电场的增加迁移率的降低也会引起导通电阻的变化。

双极器件中，当器件导通时，电子和空穴充满了低搀杂层，此时注入的载流子密度比搀杂浓度还要高，这个区域的电阻明显的降低了。

在区域边界X 1到X r ，面积为A 的区域电阻由下式表示：⎰+=rX X p n p n qA dx R 1)(μμ 这里n 和p 分别是电子和空穴的密度，n μ和p μ是载流子的迁移率，载流子并不是均匀分布的，它们的密度也不是均匀的。

(2). 电荷存储量对于双极型器件而言，当处于导通状态时，载流子电荷被存储在低搀杂区域，这些载流子电荷在器件阻断之前，必须尽快地被移走，这过程是引起开关延时和开关损耗的根本原因。

控制系统仿真与cad第二版课后答案王燕平1、47.下列选项中，不属于显示器主要技术指标的是（）[单选题] *A.分辨率B.重量(正确答案)C.像素的点距D.显示器的尺寸2、15.下列叙述中，错误的是（）。

[单选题] *A.硬磁盘可以与CPU之间直接交换数据(正确答案)B.硬磁盘在主机箱内，可以存放大量文件C.硬磁盘是外存储器之一D.硬磁盘的技术指标之一是每分钟的转速rpm3、光纤通信中使用的复用方式是（）。

易[单选题] *A. 时分多路B. 空分多路C. 频分多路D. 波分多路(正确答案)4、1根据冯.诺依曼原理，准确地说计算机硬件各部件如何动作具体是由（）决定的。

[单选题] *A.存储器中的内容B.操作系统C.用户D.CPU 中所执行的指令(正确答案)5、C：程序设计语言D：符号语言一个完整的计算机系统的组成部分的确切提法应该是______。

[单选题] *A：计算机主机、键盘、显示器和软件B：计算机硬件和应用软件6、在Windows 的"我的电脑" 窗口中，若已选定了文件或文件夹，为了设置其属性，可以打开属性对话框的操作是（）。

[单选题] *A．用鼠标右键单击“文件”菜单中的“属性”命令B．用鼠标右键单击该文件或文件夹名，然后从弹出的快捷菜单中选“属性”项(正确答案) C．用鼠标右键单击“任务拦”中的空白处，然后从弹出的快捷菜单中选择“属性”项D．用鼠标右键单击“查看”菜单中“工具栏”下的“属性”图标7、在WPS文字中,编排完一个文件后,要想知道其打印效果,可以（）。

[单选题] *A.选择“模拟显示”命令B.选择“打印预览”命令(正确答案)C.按F键8、45.下列关于电子邮件的说法，正确的是（）。

[单选题] *A.收件人必须有E-mai1地址，发件人可以没有E-mai1地址B.发件人必须有E-mai1地址，收件人可以没有E-mai1地址C.发件人和收件人都必须有E-mai1地址(正确答案)D.发件人必须知道收件人的邮政编码9、20.计算机网络是一个（）。

习题如图7所示的基于视觉传感器的磁悬浮轴承运动控制系统，它将磁悬浮轴承的控制问题抽象成利用电磁铁将一个金属轴杆悬挂起来的问题。

此系统能磁悬浮轴承实验装置图7 基于视觉传感器的磁悬浮轴承系统结构图够利用摄像头采集的图像计算出轴杆与电磁铁之间的距离x ，然后由计算机计算出为实现平衡所需要的控制电流，最后通过驱动电路对驱动电流实施控制。

假设令0i 表示平衡点的工作电流，则本运动控制系统的微分方程可以在平衡点处线性化为：12mxk x k i =+ 一组合理的参数值为0.02kg m =，120N/m k =，20.4N/A k =。

试解答下述问题：(1) 用极点配置法为此磁悬浮轴承控制系统设计一个控制器，使闭环系统满足如下指标：调整时间0.25s s t ≤，对x 的某一初始偏移的超调量小于20％。

(2) 为你设计的系统绘制关于1k 的根轨迹，并讨论能否用你设计的闭环系统平衡各种质量的轴杆。

(3) 假定一个初始阶跃位移扰动作用于轴杆，且传感器仅能测量距离工作点1/4cm ±范围内的x ，驱动电路仅能提供1A 的控制电流，那么可能控制的最大位移是多少？解答(1)由磁悬浮轴承运动控制系统的微分方程可知，它是一个典型的二阶系统，所以系统的特征方程在复平面上有两个根。

根据自动控制理论，若要使系统具有稳定的暂态响应，需要两个根的实部小于0，也就是都位于左半平面。

同时，根据系统暂态响应特性的不同，系统可以分为四种状态：临界阻尼状态、过阻尼状态、无阻尼状态以及欠阻尼状态。

临界阻尼和过阻尼状态下的二阶系统响应虽无超调，但过程缓慢，不适用于像磁悬浮轴承系统这种快速响应的场合。

无阻尼状态下系统的暂态响应将是恒定振幅的周期函数，不会收敛到稳定值，也不适用于磁悬浮轴承控制系统。

因此，需要设计一个使系统处于欠阻尼状态的控制器。

设控制系统特征方程的两个极点分别为a ，b ，则系统的特征方程可以写成如下形式：2220n n s s ξωω++= (1)其中，()/a b ξ=-+n ω=，是二阶系统的无阻尼自然振荡角频率。

计算机仿真技术与CAD第三版课后练习题含答案前言计算机仿真技术与CAD第三版是近年来广受欢迎的CAD教材。

本书涵盖了CAD 的基础知识、CAD的应用、计算机仿真技术等内容，适合CAD初学者和高级用户使用。

本文主要介绍本书的课后练习题及答案，希望对大家有所帮助。

第一章课后练习题及答案1.1 选择题1.以下哪项不是CAD应用的领域？A. 工程学B. 数学C. 植物学D. 医学答案：C2.CAD是什么的缩写？A. Computer-ded Design（计算机辅助设计）B. Computer-ded Distance（计算机辅助距离）C. Capsule-med Destruction（胶囊定向破坏）D. Capability-Assisted Driving（能力协助驾驶）答案：A1.2 填空题1.CAD软件主要分为三种类型，分别是2D __________ 、3D__________ 、仿真软件。

答案：制图软件、制造软件2.AutoCAD是一款 __________ 系统。

答案：二维制图3.CAD软件通常可以完成______________ 、___________ 、____________ 等操作。

答案：制图、设计、分析1.3 综合应用题设计一件简单的图形，包括设备名称、设备布局图、视图图和详细图。

并在AutoCAD中实现。

答：见图。

简单图形简单图形图中的简单图形包括了设备名称、设备布局图、视图图和详细图。

其中，设备名称为“简单图形”，设备布局图包含了该设备的位置、大小和方向信息，视图图包含了设备的3D视角，详细图则针对设备内部的结构进行了细节描绘。

使用AutoCAD可以方便地实现这件图形的设计和制图。

第二章课后练习题及答案……第三章课后练习题及答案……第四章课后练习题及答案……结语以上是《计算机仿真技术与CAD第三版》的部分课后练习题及答案，希望对大家有所帮助。

本书内容丰富、简明易懂，是一本非常优秀的CAD教材。

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

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

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

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

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

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

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

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

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

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

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

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

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

控制系统计算机仿真及辅助设计课程上机作业1月5日上机作业一：2-1 思考题：（1）数学模型的微分方程，状态方程，传递函数，零极点增益和部分分式五种形式，各有什么特点？（2）数学模型各种形式之间为什么要互相转换？（3）控制系统建模的基本方法有哪些？他们的区别和特点是什么？（4）控制系统计算机仿真中的“实现问题”是什么含意？（5）数值积分法的选用应遵循哪几条原则？（1）微分方程是直接描述系统输入和输出量之间的制约关系，是连续控制系统其他数学模型表达式的基础。

状态方程能够反映系统内部各状态之间的相互关系，适用于多输入多输出系统。

传递函数是零极点形式和部分分式形式的基础。

零极点增益形式可用于分析系统的稳定性和快速性。

利用部分分式形式可直接分析系统的动态过程。

（2）不同的控制系统的分析和设计方法，只适用于特定的数学模型形式。

（3）控制系统的建模方法大体有三种：机理模型法，统计模型法和混合模型法。

机理模型法就是对已知结构，参数的物理系统运用相应的物理定律或定理，经过合理的分析简化建立起来的各物理量间的关系。

该方法需要对系统的内部结构和特性完全的了解，精度高。

统计模型法是采用归纳的方法，根据系统实测的数据，运用统计规律和系统辨识等理论建立的系统模型。

该方法建立的数学模型受数据量不充分，数据精度不一致，数据处理方法的不完善，很难在精度上达到更高的要求。

混合法是上述两种方法的结合。

（4）“实现问题”就是根据建立的数学模型和精度，采用某种数值计算方法，将模型方程转换为适合在计算机上运行的公式和方程，通过计算来使之正确的反映系统各变量动态性能，得到可靠的仿真结果。

（5）数值积分法应该遵循的原则是在满足系统精度的前提下，提高数值运算的速度和并保证计算结果的稳定。

2-2.用matlab语言求下列系统的状态方程、传递函数、零极点增益、和部分分式形式的模型参数，并分别写出其相应的数学模型表达式：(1) G（s）=324327242410355024s s ss s s s+++++++(2) Equation Section (Next).X=2.25 -5 -1.25 -0.542.25 -4.25 -1.25 -0.2520.25 -0.5 -1.25 -121.25 -1.75 -0.25 -0.75 0X⎡⎤⎡⎤⎢⎥⎢⎥⎢⎥⎢⎥+⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦uy=[0 2 0 2] X解：（1）程序代码如下：clear;clc;num=[1 7 24 24];den=[1 10 35 50 24];[Z,P,K]=tf2zp(num,den);[A,B,C,D]=zp2ss(Z,P,K);[R,P1,H]=residue(num,den);G1=tf(num,den)G2=zpk(Z,P,K)G3=ss(A,B,C,D)程序的输出结果如下：Transfer function:s^3 + 7 s^2 + 24 s + 24---------------------------------s^4 + 10 s^3 + 35 s^2 + 50 s + 24Zero/pole/gain:(s+1.539) (s^2 + 5.461s + 15.6)-------------------------------(s+4) (s+3) (s+2) (s+1)a =x1 x2 x3 x4 x1 -3 -1.414 0 0x2 1.414 0 0 0x3 1 1.088 -7 -3.464x4 0 0 3.464 0b =u1x1 1x2 0x3 0x4 0c =x1 x2 x3 x4 y1 1 1.088 -1.539 1.038d =u1y1 0Continuous-time model.（2）程序代码如下：clear;clc;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=[4220];C=[0 2 0 2];D=0;[num,den]=ss2tf(A,B,C,D);[Z,P,K]=ss2zp(A,B,C,D);[R,P1,H]=residue(num,den);G1=tf(num,den)G2=zpk(Z,P,K)G3=ss(A,B,C,D)程序输出结果如下：Transfer function:4 s^3 + 14 s^2 + 22 s + 15--------------------------------------s^4 + 4 s^3 + 6.25 s^2 + 5.25 s + 2.25Zero/pole/gain:4 (s+1.5) (s^2 + 2s + 2.5)--------------------------(s+1.5)^2 (s^2 + s + 1)a =x1 x2 x3 x4x1 2.25 -5 -1.25 -0.5x2 2.25 -4.25 -1.25 -0.25x3 0.25 -0.5 -1.25 -1x4 1.25 -1.75 -0.25 -0.75b =u1x1 4x2 2x3 2x4 0c =x1 x2 x3 x4y1 0 2 0 2d =u1y1 0Continuous-time model.2-3.用欧拉法求下面系统的输出响应y(t)在0≤t≤1上，h=0.1时的数值。

《控制系统仿真与CAD》习题集第1章仿真软件——MATLAB1-1：利用MATLAB的文本文件求函数y=1-2：利用MATLAB的函数文件，求函数：2??y1=3x1+x2+x3 ?2=??y3xxx?123?2x3在x=?4时的值。

在x1=?2,x2=3,x3=1时的值。

1-3：利用图形窗口分割方法将下列极坐标方程：ρ=sin(2θ)cos(2θ)用四种绘图方式画在不同的窗口中。

1-4：绘制函数图形：y(t)=1?2e?tsin(t)（0≤t≤30）且在x轴写上“Time”标号，y轴写上“Amplitude”标号，图形的标题为“Decaying-oscillating Exponential”。

1-5：求S=∑i的值。

i=11001-6：求满足∑i&gt;1000的最小m值。

i=1m1-7：分别利用for和while循环语句求满足S=∑i2&lt;1000的m 的最大值。

i=1m1-8：求函数f(x)=3x?x3的最小值。

1-9：求以下非线性方程组的解。

?sinx+y2+lnz?7=0?3 ?3x+2y?z+1=0?x+y+z?5=0?1-10：求以下线性方程组的解。

?x?y=0 ?xy60+?=???(0)=0下的解。

1-11：求下列微分方程在初始条件x(0)=1,x??????+x=0 x?(1?x2)x?12?1-13：对于矩阵A=??，MATLAB以下四条命令：34??1-12 ：求积分y=∞?x∫e2/2dx的解。

A.^0.5;A^0.5;sqrt(A);sqrtm(A)所得结果相同吗？哪些结果是复数矩阵？1-14：一个3位整数各位数字的立方和等于该数本身则称该数为水仙花数，输出100~999之间的全部水仙花数。

1-15：编写函数文件求半径为r的圆的面积和周长。

1-16：已知y=1111+2+2+????2，当n=100时，求y的值。

2123n1-17：利用函数的递归调用，求n!。