系统建模与仿真课后作业
系统建模与仿真习题2及答案
系统建模与仿真习题二及答案1. 考虑如图所示的典型反馈控制系统框图(1)假设各个子传递函数模型为66.031.05.02)(232++-+=s s s s s G ，s s s G c 610)(+=，21)(+=s s H 分别用feedback （）函数以及G*Gc/(1+G*Gc*H)（要最小实现）方法求该系统的传递函数模型。
(2) 假设系统的受控对象模型为s e s s s G 23)1(12)(-+=，控制器模型为 ss s G c 32)(+=，并假设系统是单位负反馈，分别用feedback （）函数以及G*Gc/(1+G*Gc*H)（要最小实现）方法能求出该系统的传递函数模型？如果不能，请近似该模型。
解：（1）clc;clear;G=tf([2 0 0.5],[1 -0.1 3 0.66]);Gc=tf([10 6],[1 0]);H=tf(1,[1 2]);G1=feedback(G*Gc,H)G2=G*Gc/(1+G*Gc*H)Gmin=minreal(G2)结果：Transfer function:20 s^4 + 52 s^3 + 29 s^2 + 13 s + 6s^5 + 1.9 s^4 + 22.8 s^3 + 18.66 s^2 + 6.32 s + 3Transfer function:20 s^8 + 50 s^7 + 83.8 s^6 + 179.3 s^5 + 126 s^4 + 57.54 s^3 + 26.58 s^2 + 3.96 ss^9 + 1.8 s^8 + 25.61 s^7 + 22.74 s^6 + 74.11 s^5 + 73.4 s^4 + 30.98 s^3+ 13.17 s^2 + 1.98 s Transfer function:20 s^4 + 52 s^3 + 29 s^2 + 13 s + 6s^5 + 1.9 s^4 + 22.8 s^3 + 18.66 s^2 + 6.32 s + 3(2)由于s c e s s s s G s G 232)1(3624)(*)(-++= 方法1：将s e 2-转换为近似多项式。
系统建模控制与仿真 习题及解答
1、工业控制系统可分为几种大类型，各有什么特点?适合的应用领域。
答：工业控制系统可分为分布式控制系统（DCS ）和可编程逻辑控制器（PLC ）两大类型。
分布式控制系统（DCS ））是以微处理机为基础，以危险分散控制，操作和管理集中为特性的新型控制系统，它具有高可靠性、开放性、灵活性、协调性、易于维护、控制功能齐全等特点，属于过程控制系统，主要控制手段是PID ，适用于流程工业； 可编程逻辑控制器（PLC ）是一种专门为在工业环境下应用而设计的数字运算操作的电子装置，它的特点有可靠性高、抗干扰能力强、硬件配套齐全、功能完善、适用性强、易学易用、容易改造、体积小、重量轻、能耗低等，属于离散控制系统，主要控制手段是顺序与逻辑控制，适用于制造业，目前，已广泛应用于钢铁、石油、化工、电力、建材、机械制造、汽车、轻纺、交通运输、环保及文化娱乐等各个行业。
2、根据自己的理解简述现代控制理论的发展历史，分析为什么现代控制理论在过程控制系统中难以应用？答：现代控制理论是为了分析多输入多输出系统、非线性系统和时变系统而出现的，先是贝尔曼等人提出状态分析法，接着卡尔曼等人提出状态空间法，后来，罗森布洛克等人将经典控制理论传递函数的概念推广到多变量系统，并探讨了传递函数矩阵与状态方程之间的等价转换关系，与此同时，系统辨识、最优控制、离散时间系统和自适应控制的发展大大丰富了现代控制理论的内容。
现代控制理论之所以难以在过程控制系统中应用，有以下几个主要原因：1) 现代控制理论获得较好效果的前提是系统内部结构参数完全已知，并且很精确，而过程控制系统中系统参数一般都是经常变化的，在每次重新开机之后系统参数都会发生变化2) 现代控制理论对系统状态变量采用的是微分运算，对各种干扰非常敏感，而过程控制系统中的干扰非常多，因此在过程控制系统中运用现代控制理论经常得不到好的结果，甚至出现错误的结果，现代控制理论适用于航天、实验室等干扰很少的场合3) x Ax Bu ∙=+；y Cx =中，输出并不在闭环内部，因此现代控制理论不能保证输出具有稳定性4) 现代控制理论在解决现实问题中为了保证系统的稳定性，一般会采用模糊化，这样做的代价是大大降低系统的精确性、快速性5) 在过程控制系统中应用现代控制理论后得到的结果往往不如直接运用PID 后得到的结果好3、什么是串级控制系统？其两个回路各有什么特点？答：串级控制系统是两只调节器串联起来工作，其中一个调节器的输出作为另一个调节器的给定值的系统。
《系统建模与仿真》作业题
第1题——最小二乘法的具体应用例子[题目]已知某一单输入单输出线性系统的差分方程形式为101()(1)()(1)()y k a y k b u k b u k k ξ=--++-+但其参数1a ，0b ，1b 为未知数，且()k ξ为不相关的随机序列。
经过辨识试验，测得5组输入输出数据为(1) 1.5u =，(2) 3.8u =-，(3) 3.1u =，(4) 2.99u =，(5) 5.12u =和(1)0.8y =，(2)0.1y =-，(3)0.45y =，(4)0.34y =，(5)0.12y =。
试求出其最优参数估计。
解：编写MATLAB 程序如下：u(1)=1.5;u(2)=-3.8;u(3)=3.1;u(4)=2.99;u(5)=5.12;y(1)=0.8;y(2)=-0.1;y(3)=0.45;y(4)=0.34;y(5)=0.12;c=10000;P=diag([c,c,c]);f=[-y(1) u(2) u(1)];K=P*f'*inv((1+f*P*f'));x=K*y(3);for n=1:1:2f=[-y(n+1) u(n+2) u(n+1)];%f2P=P-P*f'*inv((1+f*P*f'))*f*P;%P1K=P*f'*inv((1+f*P*f'));%K2x=x+K*(y(n+3)-f*x);%x2end其中x 中存储最终的结果：a1=-0.023;b0=-0.018;b1=0;第2题——以M 文件的形式编写一个MATLAB 仿真程序[题目]给定被控系统模型(1)()0.8()()y k y k u k f k +=++，其中干扰()0.20.1sin 0.01f k k =+，初值(0)0.2y =，(0)0u =。
试编写一个M 文件，对该系统在PD 控制律(1)()()u k u k u k +=+∆12()[()()]{[(1)()][(1)()]}r r r u k B y k y k B y k y k y k y k ∆=-++--+-作用下的输出()y k 进行仿真，将系统的实际输出()y k 与参考输出()r y k 画在同一张坐标图上。
(完整)系统建模与仿真习题答案(forstudents)
第一章习题1-1什么是仿真？它所遵循的基本原则是什么？答：仿真是建立在控制理论，相似理论，信息处理技术和计算技术等理论基础之上的,以计算机和其他专用物理效应设备为工具，利用系统模型对真实或假想的系统进行试验，并借助专家经验知识，统计数据和信息资料对试验结果进行分析和研究，进而做出决策的一门综合性的试验性科学。
它所遵循的基本原则是相似原理。
1-2在系统分析与设计中仿真法与解析法有何区别?各有什么特点?答：解析法就是运用已掌握的理论知识对控制系统进行理论上的分析，计算。
它是一种纯物理意义上的实验分析方法,在对系统的认识过程中具有普遍意义。
由于受到理论的不完善性以及对事物认识的不全面性等因素的影响，其应用往往有很大局限性.仿真法基于相似原理，是在模型上所进行的系统性能分析与研究的实验方法.1-3数字仿真包括那几个要素？其关系如何？答: 通常情况下，数字仿真实验包括三个基本要素，即实际系统，数学模型与计算机。
由图可见，将实际系统抽象为数学模型，称之为一次模型化，它还涉及到系统辨识技术问题,统称为建模问题；将数学模型转化为可在计算机上运行的仿真模型，称之为二次模型化，这涉及到仿真技术问题，统称为仿真实验.1—4为什么说模拟仿真较数字仿真精度低？其优点如何?.答：由于受到电路元件精度的制约和容易受到外界的干扰,模拟仿真较数字仿真精度低但模拟仿真具有如下优点：（1）描述连续的物理系统的动态过程比较自然和逼真。
（2）仿真速度极快，失真小，结果可信度高。
（3）能快速求解微分方程.模拟计算机运行时各运算器是并行工作的，模拟机的解题速度与原系统的复杂程度无关.（4）可以灵活设置仿真试验的时间标尺，既可以进行实时仿真，也可以进行非实时仿真.（5）易于和实物相连。
1-5什么是CAD技术？控制系统CAD可解决那些问题？答：CAD技术，即计算机辅助设计（Computer Aided Design)，是将计算机高速而精确的计算能力,大容量存储和处理数据的能力与设计者的综合分析，逻辑判断以及创造性思维结合起来，用以加快设计进程，缩短设计周期，提高设计质量的技术.控制系统CAD可以解决以频域法为主要内容的经典控制理论和以时域法为主要内容的现代控制理论。
《系统建模与仿真》作业题1
《系统建模与仿真》作业第一次（共两次）布置作业时间：“经典建模法”结束后 要求交作业时间：从布置日开始不超过1周 作业量：共3道题第1题——体现电气系统经典建模[题目] 在如图所示的电路中，R 表示一个电阻，L 表示一个电感，C 表示一个电容，i 表示电流强度，u 表示输入电压，c u 表示电容器的输出电压。
试写出一个状态空间数学模型。
u图1 典型的RLC 电路第2题——体现机械系统经典建模[题目]如图2是一个文字处理器打印轮轴控制系统的简化图。
打印轮轴由一个直流电动机通过皮带和滑轮进行控制。
假设皮带是刚性的，电动机与皮带轮之间的连接也是刚性的，并定义如下的参数和变量：m ()T t 是电动机的力矩；m ()t 是电动机的角位移；()y t 是打印轮轴的线性位移；m J 式电动机的惯量；m B 是电动机的粘性摩擦系数；K 是扭转轴的刚性系数；r 是滑轮的半径；M 是打印轮轴的质量。
（1）写出系统的微分方程；（2）写出系统的传递函数模型m ()()Y s T s 。
T打印轮轴图2 打印轮轴控制系统简化图第3题——体现热工过程经典建模汽轮机高压加热器疏水系统的原理框图如图3所示，其中各段抽汽的压力大小关系为321p p p >>，抽汽温度大小关系为321T T T >>。
给水流量w 和给水温度T 一般来说为两个随机变量。
三个疏水管道阀门的开度为归一化量，即]1,0[,,321∈u u u 。
三个高压加热器的疏水水位分别为1y ，2y ，3y 。
它们的关系可描述为),,,(111111T w T p f u k y+-= ),,,(22223122T w T p f u k u k y+-= ),,,(T w T p f u k u k y+-=图3 汽机高加疏水系统原理框图式中的),,,(111T w T p f ，),,,(222T w T p f ，),,,(333T w T p f 表示系统的不确定干扰，1k ，2k ，3k ，4k ，5k 表示适当的正常数。
《控制系统建模与仿真》课后习题-2021版
《控制系统建模与仿真》课程习题（1）一、“投针实验”的历史价值在人类数学文化史中，对圆周率 精确值的追求吸引了许多学者的研究兴趣。
在众多的圆周率计算方法中，最为奇妙的是法国物理学家布丰（Boffon）在1777年提出的“投针实验”。
试回答下列问题：1、试对“投针实验”的机理给出一种直观形象的物理解释？2、有人说“布丰/ Boffon（投针实验）是仿真技术的奠基者”，为什么？3、试用MATLAB语言编制“投针实验”的仿真程序，仿真证明之。
二、自平衡式两轮电动车的安全问题近年来，自平衡式两轮电动车产品成为“抢眼”的代步工具，但也出现很多问题（如上图所示）；试根据你所了解的情况就“平衡车产品是否可以合法上路？”问题，给出你的意见与建议。
提示：可从“技术、安全、法律、可持续”等方面，有理有据地展开讨论。
参考书：张晓华《控制系统数字仿真与CAD》 (第4版) 机械工业出版社 2020张晓华《系统建模与仿真》（第2版）清华大学出版社 2016《控制系统建模与仿真》课程习题（2）一、一阶直线倒立摆系统的建模问题对于教材中图2-7所示的一阶直线倒立摆系统，基于牛顿定律所建立的数学模型（如教材的图2-8所示），试问：这个数学模型是否正确，给出你的分析与证明。
提示：（1）基于MATLAB仿真进行模型验证（参见教材第四章第三节）；（2）应用“拉格朗日方程”方法建模，进行结果对比。
二、一阶直线双倒立摆系统的可实现问题如下图所示的一阶直线双倒立摆系统，试问：能否通过控制力F实现“在保持两杆不倒的条件下，使小车在直线X方向的位置任意移动”？提示：（1）建立系统数学模型；（2）应用现代控制理论的“能控性定理”进行分析。
参考书：张晓华《控制系统数字仿真与CAD》 (第4版) 机械工业出版社 2020张晓华《系统建模与仿真》（第2版）清华大学出版社 2016《控制系统建模与仿真》课程习题（3）一、水箱液位控制系统设计问题如下图所示的“水箱液位系统”，试回答下列问题：1、试给出含有（控制器+传感器）的“水箱液位控制系统”方案；2、试依据“流体力学”的基本概念，建立系统的数学模型；3、若使系统液位控制实现稳态无静差，试给出PID控制器设计方案；二、水箱液位控制的拓展问题试回答下述问题：1、某人在上述“水箱液位控制系统”中，采用单片机作控制器，程序设计为“增量式PI控制算法”，如果控制系统在“阶跃给定”下存在稳态误差，试问这种情况是否合理？为什么？2、对于上图所示的“水箱液位系统”，在下排水出口处流体呈“紊流”状态,试证明：其流量与液位高度的关系为Q=K∙√H。
数学建模与系统仿真网课答案.docx
数学建模与系统仿真网课答案问：计算机的运算速度只与机器的主频相关。
答：错问：玩物丧志这句话对于创新人才来说是不正确的。
（）答：√问：在决定是否相信一个断言时不必考察它和我们的背景信息是否冲突。
（）答：错问：存储器的容量应该包括主存容量和辅存容量。
答：对问：代表宋代写实绘画标本的画家是：（）答：郭熙问：就规模和影响力而言，成长最快的信息来源是（）答：谷歌问：当流言这样一种信息如果与群体无意识结合的时候，就容易形成一种（）社会心理答：集群的冲动问：中医“四大经典”包括（）。
答：《黄帝内经》《黄帝八十一难经》《神农本草经》《伤寒杂病论》问：中国梦的实现路径不包括( )。
答：实现中国梦必须加快经济体制改革问：《人生》中暗示了一种青年文学主题的转折,曾经作为改天换地的主力军青年一代,从外部的世界中回到自己的个人世界,他们将依照个人的利益行使自己的主动权。
()答：正确问：“海权论”理论的提出者是（）。
答：A问：控制器用来完成算术运算和逻辑运算。
答：错问：考古中，区分文化类型的区系类型表示的意义是:（）答：区代表空间，系代表时间问：中华民族近代以来最伟大的梦想,就是（）。
答：B问：沙皇即俄罗斯最高统治者。
答：错问：中唐时期诗坛上出现了比盛唐时期更多的风格流派，具体来说主要有（）答：现实主义诗派浪漫诗派田园诗派边塞诗派问：1978年成龙主演的哪两部喜剧片标志着功夫喜剧片的开端？答：《蛇形刁手》《醉拳》问：输入设备将机器运算结果转换成人们熟悉的信息形式。
答：错问：购物车易用的实现方式答：1. 1.“放入购物车”的按钮必须十分明显。
2.购物车按钮文字。
3.随时放入购物车。
4.编辑购物车。
5.非注册用户。
6.随时随地查看购物车内容。
7.尽可能提供各种可能的付款及配送方式。
8.不要打扰用户付款。
9.在购物过程中尽早显示产品价格。
10.订单信息完整。
11.为所有购买用户建立账号 12.付款进程提示。
问：人们与众不同的生活方式，本质是来自于他们独特的（）答：心智模式。
系统建模与仿真 作业
>> [m,n]=size(A) m = 8 n = 8 >> B=A(2:2:floor(m/2)*2,:)
B = 9 40 41 8 2-3 建立一个m文件，如下： function y=fenduan(x,D,h) if(x>D) y=h; elseif(x<=D&x>=-D); y=h/(D*x) else y=(-h); end 在命令窗口作如下命令输入： Command window： >> x=2 x = 2 >> fenduan(x,D,h) ans = 0.1000 >> x=-6 x = -6 55 26 23 58 54 27 22 59 12 37 44 5 13 36 45 4 51 30 19 62 50 31 18 63 16 33 48 1
3 4 2
2 1 1
2-31: clear all clc x=1:24; y=[14 14 14 14 15 16 18 20 22 23 25 28 31 32 31 29 27 25 24 22 20 18 17 16]; plot(x,y,'o');hold on p1=polyfit(x,y,3);%三次多项式拟合 y1=polyval(p1,x); plot(x,y1,'r') p2=polyfit(x,y,4);%四次多项式拟合 y2=polyval(p2,x); plot(x,y2,'b') p3=polyfit(x,y,5);%五次多项式拟合 y3=polyval(p3,x); plot(x,y3,'m') p4=polyfit(x,y,6);%6次多项式拟合 y4=polyval(p4,x); plot(x,y4,'g')
浙江工商大学系统建模与仿真课后习题答案
参考教材：离散事件系统建模仿真及GPSSWorld 教程，译作者：谢毅缪亚萍，出版社：清华大学出版社，年代：2011部分习题答案：（具体题目信息省略）第二章：2.16（双理发师手工仿真）理发师1利用率：49/58理发师2利用率：28/58时钟系统状态时间表统计计数器StQ （t ）StS1（t ）StS2（t ）C B1B2W 0000A （1,5）00005010D （1,23），A （2,11）00011012D （1,23），D （2,31），A （3,19）060019112D （1,23），D （2,31），A （4,38）0148023032D （2,31），D （3,38），A （4,38）11812431030D （3,38），A （4,38）2262043800A （4,38）333204040D （4,54），A （5,50）33320450045D （4,54），D （5,66），A （6,58）34520454005D （5,66），A （6,58）4492445865D （6,74），D （5,66），A （7,…）449284平均等待时间：4/6（min）第三章：3.3（货物出库）GENERATE10,5 TERMINATE10 GENERATE15 TERMINATE20 GENERATE30,10 TERMINATE50START20003.4（零件加工）GENERATE20,5ADD1QUEUE Q_ASEIZE ADEPART Q_AADVANCE16,5RELEASE AQUEUE Q_BSEIZE BDEPART Q_BADVANCE15,10RELEASE BQUEUE Q_CSEIZE CDEPART Q_CADVANCE10,2RELEASE C TRANSFER0.05,ADD2,ADD1ADD2TERMINATEGENERATE60TERMINATE1START40（仿真4个小时）若为生产200个合格零件，则ADD2改为：ADD2TERMINATE1START200（仿真4个小时）3.6（加工中心）S_A STORAGE2S_B STORAGE1S_C STORAGE3S_D STORAGE5GENERATE20,10QUEUE Q_AENTER S_ADEPART Q_AADVANCE30,15LEAVE S_ATRANSFER0.6,ADDC,ADDBADDC ENTER S_CADVANCE70,20LEAVE S_CTRANSFER,ADDDADDB ENTER S_BADVANCE20,10LEAVE S_BTRANSFER0.2,ADDD,ADDC ADDD ENTER S_DADVANCE90,30LEAVE S_DTERMINATEGENERATE60TERMINATE1START1003.7（汽车清洗店，顾客有偏好）注：可用BOTH或是ALL模式，但是相对比较复杂，可用SELECT，实现更简单。
系统建模与仿真习题1及答案
系统建模与仿真习题一及答案1. 有源网络如图所示（1） 列些输出0u 与输入1u 之间的微分方程。
（2） Ω=101R 、Ω=52R 、Ω=23R 、Ω=34R 、F C 2=，在零初始条件下，将（1）中的微分方程表示为传递函数、状态空间形式、零极点增益形式。
（3）求（2）中方程在输入1u 为单位阶跃响应下的输出曲线。
解：（1） 由运算放大器的基本特点以及电压定理)4()3()(1)2()()1(2132021421320111R i R i u dt i i Cu R i i u R i u R i u c c -=+=+++==⎰（3）式代入（2）式得：42121320)()(1R i i dt i i C R i u ++++=⎰ （5）消去中间变量21,i i 有13142430114131230111120)(1u R R R R R R u u R R dt u R R R R u R u C u R R u ++++++=⎰ 两边求导整理后得（2）代入数据可以得到微分方程为：11007.02.610u u u u--=+ 程序如下：clc;clear;num=[-6.2 -0.7]; den=[10 1]; Gtf=tf(num,den) Gss=ss(Gtf) Gzpk=zpk(Gtf)结果：Transfer function: -6.2 s - 0.7 ------------ 10 s + 1状态空间形式： a =x1 x1 -0.1 b =u1 x1 0.125 c =x1 y1 -0.064 d =u1 y1 -0.62Continuous-time model.Zero/pole/gain: -0.62 (s+0.1129) ---------------- (s+0.1)（3）由（2）知系统的传递函数为-6.2 s - 0.7 ------------ 10 s + 1系统的输入信号为单位阶跃函数，则其Laplace 变换为1/s ，这样系统的输出信号的Laplace 变换为Y(s)=-6.2 s - 0.7 ------------ 10 s^2 + s编写程序，将其表示为（R,P,Q ）形式 clc;clear; s=tf('s')Gtf=(-6.2*s-0.7)/(10*s^2+s) [num,den]=tfdata(Gtf,'v') [R,P,Q]=residue(num,den) R =0.0800 -0.7000 P =-0.1000 0 Q = []于是得到：7.008.0)(1.0-=-t e t y 绘制曲线程序： clc;clear; t=0:0.1:100;y=0.08*exp(-0.1*t)-0.7; plot(t,y)2.已知系统的框图如下：其中：G1=1/(s+1)，G2=s/(s^2+2)，G3=1/s^2，G4=(4*s+2)/(s+1)^2，G5=(s^2+2)/(s^3+14)。
系统建模与仿真习题3及答案
系统建模与仿真习题三及答案1.已知系统)24(32)(21+++=s s s s s G 、2103)(2+-=s s s G 求G 1(s)和G 2(s)分别进行串联、并联和反馈连接后的系统模型。
解：clc;clear;num1=[2 3];den1=[1 4 2 0];num2=[1 -3];den2=[10 2];G1=tf(num1,den1);G2=tf(num2,den2);Gs1=series(G1,G2)Gp1=parallel(G1,G2)Gf=feedback(G1,G2)结果：Transfer function:2 s^2 -3 s - 9------------------------------10 s^4 + 42 s^3 + 28 s^2 + 4 sTransfer function:s^4 + s^3 + 10 s^2 + 28 s + 6------------------------------10 s^4 + 42 s^3 + 28 s^2 + 4 sTransfer function:20 s^2 + 34 s + 6--------------------------------10 s^4 + 42 s^3 + 30 s^2 + s – 92.某双闭环直流电动机控制系统如图所示：利用feedback( )函数求系统的总模型。
解：模型等价为：编写程序：clc;clear;s=tf('s');G1=1/(0.01*s+1);G2=(0.17*s+1)/(0.085*s);G3=G1;G4=(0.15*s+1)/(0.051*s);G5=70/(0.0067*s+1);G6=0.21/(0.15*s+1);G7=(s+2)/s;G8=0.1*G1;G9=0.0044/(0.01*s+1);sys1=feedback(G6*G7,0.212);sys2=feedback(sys1*G4*G5,G8*inv(G7)); sys=G1*feedback(sys2*G2*G3,G9)结果：Transfer function:3.749e-005 s^6 + 0.008117 s^5 + 0.5024 s^4 + 6.911 s^3 + 36.57 s^2 + 78.79 s + 58.8 -------------------------------------------------------------------------------------------------------4.357e-014 s^10 + 2.432e-011 s^9 +5.43e-009 s^8 +6.303e-007 s^7 + 4.145e-005 s^6 + 0.001578 s^5 + 0.03217 s^4 + 0.2098 s^3 + 0.4116 s^2 + 0.3467 s + 0.2587根据需要可忽略高阶项。
西工大2022年4月机考《系统建模与仿真》作业参考答案
西工大2022年4月机考《系统建模与仿真》作业参考答案试卷总分:100 得分:100本科目3次作答机会，每次试题内容相同，只是题目和选项顺序是随机调整的，大家可放心下载使用一、单选题(共20 道试题,共40 分)1.数学模型根据模型的状态变量可以分为（）。
A.连续变化模型和离散变化模型B.连续时间模型和离散时间模型C.确定性模型和随机性模型D.同构模型和同态模型正确答案:A2.在仿真模型一样，所要仿真的时间长度也一样的情况下，采用（）可获得最高的效率。
A.固定步长时间推进机制B.下次事件时间推进机制C.混合时间推进机制D.随机步长时间推进机制正确答案:B3.忽略具体事物的特殊性，着眼于整体和一般规律，这种研究方法是（）。
A.抽象B.归纳C.演绎D.推导正确答案:A4.（）是把过程调用和响应调用执行码结合在一起的过程A.汇编B.联编C.调试D.执行正确答案:B5.在系统与模型之间，如果在行为一级等价，则称之为（）。
A.同构模型B.同态模型C.数学模型D.本构模型正确答案:B6.一种产品进入市场之后，一般会经过销售速度先不断增加然后又逐渐下降的过程，这称为产品的（）。
A.生命周期B.保质期C.生产周期D.销售周期正确答案:A7.由于大多数微分方程是求不出其解析解的，因此研究其（）和数值解法是十分重要的手段。
A.离散性B.连续性C.非稳定性D.稳定性正确答案:D8.根据事件调度法建立的仿真模型称为（）仿真模型。
A.面向事件的B.面向对象的C.面向用户的D.面向系统的正确答案:A9.能够预定事件发生时间的策略方法是（）。
A.事件调度法B.活动扫描法C.进程交互法D.结果预测法正确答案:A10.系统在有确定输入时，得到的输出却不确定，这种事物发展变化没有确定因果关系的模型是（）。
A.连续变化模型B.离散变化模型C.随机性模型D.因果模型正确答案:C11.系统数学模型的建立需要按照模型论对输入、输出状态变量及其间的函数关系进行抽象，这种抽象理论称为（）。
系统建模与仿真习题5及答案
系统建模与仿真习题五及答案1.已知系统的开环传递函数为)62)(5)(33(12)(222++++++=s s s s s s s s G ，根据相角裕度，幅值裕度判断单位负反馈下闭环系统的稳定性，并用时域响应验证结论。
解：clc;clear;num=[2 1];den=conv([1 0 0],conv([1 3 3],conv([1 5],[1 2 6])));G1=tf(num,den); margin(G1)由图知20lgh=24.8dB ，h>1；r=2.67deg 。
因此对应的闭环系统稳定。
下面由负反馈的阶跃响应验证闭环系统的稳定性clc;clear;num=[2 1];den=conv([1 0 0],conv([1 3 3],conv([1 5],[1 2 6])));G1=tf(num,den);step(feedback(G1,1))2. 已知某系统的开环传递函数为)32)(5()5(2)(2++++=s s s s s G （1）绘制系统的奈奎斯特曲线，判断闭环系统的稳定性。
（2）求出系统的单位阶跃响应，证明（1）中稳定性的判断。
解：clc;clear;num=[2 10];den=conv([1 2],[1 2 3]);G=tf(num,den);subplot(2,1,1)Nyquist(G)subplot(2,1,2)step(feedback(G,1))结论：开环传递函数位于s右半平面的极点数p=0。
由Nyquist图知：Nyquist曲线逆时针包围临界点(-1,j0)的圈数R=0。
因此，闭环正实部特征根个数Z=0，说明系统是稳定的。
仿真曲线也表明闭环系统是稳定的。
3. 已知一个离散系统的输入、输出数据如下：u=[0.9103;0.7622;0.2625;0.0475;0.7361;0.3282;0.6326;0.7564;0.9910;0.3653;0.2470 ;0.9826;0.7227;0.7534;0.6515;0.0727;0.6316;0.8847;0.2727;0.4364;0.7665;0.4777;0. 2378;0.2749]y=[0;18.4984;31.4285;32.3228;28.5690;39.1704;39.8825;46.4963;54.5252;65.9972;6 2.9181;57.5592;67.6080;70.7397;73.7718;74.0165;62.1589;63.0000;68.6356;60.8267 ;57.1745;60.5321;57.3803;49.6011]请用最小二乘法辨识出系统的脉冲传递函数模型，要求该模型的分子、分母的阶次分别为2、3次。
系统建模与仿真作业
病菌传染人数动态变化模型的仿真专业：机械电子工程姓名：王勇（10S030039）日期：2010年11月8日摘要本文利用已知的模型，运用MATLAB中Simulink工具箱对模型进行的准确的描述，然后进行仿真分析。
Simulink的每个子模型库包含有相应的功能模块，用户也可以定制和创建用户自已的模块,模型化图形输入是Simulink提供了一些按功能分类的基本的系统模块，通过对这些基本模块的调用，再将它们连接起来就可以构成所需要的系统模型，进而进行仿真与分析。
通过分析对传染情况有了准确的了解，利于对传染情况的控制。
关键字:建模,MATLAB/Simulink,分析AbstractThis paper using the known model, using MATLAB Simulink toolbox of model of accurate description, then the simulation analysis.Each submode Simulink this repository contains a corresponding function module, users can also customize and create user own module, modeling graph input is Simulink provides some according to the basic function classification system module, through to these basic modules calls, and then connect them up can form required system model, and then, a simulation and analysis.Through the analysis of infectious diseases have accurate understanding, benefit of infection status of control.Keywords：Modeling，MATLAB/Simulink，Analysis引言传染病是致病性(微)生物在人与人、动物与人及动物与动物之间相互传播的疾病,其流行既有隐蔽性又有突发性.不论急性还是慢性传染病都给人类健康带来极大灾难、给社会经济发展造成很大的损失。
系统建模与仿真课后作业
1.347
1.565
0.220
0.017
1.731
-2.831
0.285
0.343
-0.873
1.048
0.554
0.357
-0.678
1.290
0.372
0.356
-0.865
1.081
0.910
0.466
-0.424
-0.573
0.426
0.304
-0.433
-0.635
0.976
0.807
0.759766
19
2093
45
0.087891
44
1948
412
0.804688
20
228
228
0.445313
45
2063
15
0.029297
21
1143
119
0.232422
46
78
78
0.152344
22
598
86
0.167969
47
393
393
0.767578
23
433
433
0.845703
166
1318
294
0.574219
191
1433
409
0.798828
167
1473
449
0.876953
192
2048
0
0
168
2248
200
0.390625
193
3
3
0.005859
169
1003
491
0.958984
系统建模与仿真课后作业
所以，蒙特卡罗方法的基本思想就是：当试验次数充分多时，某一事件出现的频率近似等于该事件发生的概率。
一般步骤：（1）构造或描述概率过程
（2）以已知概率分布进行抽样
（3）建立各种估计量
2.8、简述离散事件系统仿真的一般步骤。
（1）阐明问题与设定目标
（2）仿真建模
29
1023
511
0.998047
5
103
103
0.201172
30
2558
510
0.996094
6
518
6
0.011719
31
2553
505
0.986328
7
33
33
0.064453
32
2528
480
0.9375
8
168
168
0.328125
33
2403
355
0.693359
9
843
331
0.646484
34
1778
242
0.472656
10
1658
122
0.238281
35
1213
189
0.369141
11
613
101
0.197266
36
948
436
0.851563
12
508
508
0.992188
37
2183
135
0.263672
13
2543
495
0.966797
38
678
166
0.324219
系统建模与仿真课后作业
、系统、模型和仿真三者之间具有怎样的相互关系答：系统是研究的对象，模型是系统的抽象，仿真通过对模型的实验以达到研究系统的目的。
、通过因特网查阅有关蒲丰投针实验的文献资料，理解蒙特卡罗方法的基本思想及其应用的一般步骤。
答：蒲丰投针实验内容是这样的：在平面上画有一组间距为a的平行线，将一根长度为L（L<a）的针任意掷在这个平面上，求此针与平行线中任一条相交的概率。
”布丰本人证明了，这个概率是：p=2L/（πa) (π为圆周率)利用这个公式可以用概率的方法得到圆周率的近似值。
所以，蒙特卡罗方法的基本思想就是：当试验次数充分多时，某一事件出现的频率近似等于该事件发生的概率。
一般步骤：（1）构造或描述概率过程（2）以已知概率分布进行抽样（3）建立各种估计量、简述离散事件系统仿真的一般步骤。
（1）阐明问题与设定目标（2）仿真建模（3）数据采集（4）仿真模型的验证（5）仿真程序的编制与校核（6）仿真模型的运行（7）仿真输出结果的统计分析、以第二章图2-5所示的并行加工中心系统为对象，试分别画出相应的实体流图和活动循环图，并比较它们两者有何区别和练习。
（1）实体流图(2)活动循环图、以第二章中图2-5所示的并行加工中心系统为对象，建立Petri 网模型。
3214Petri 网模型的运行过程，并将分析结果同例3-5相比较。
、任取一整数作为种子值，采用第三题中得到的随机数发生器生成随机数序列的前200项数据，并对其统计性能进行检验。
解：由第3题可得到一个随机数发生器： a=5 b=9 c=3 m=512取种子值，生成的随机数序列前200项数据如下： nn1500000332326458458t 4t 3 P 1 t 1P 2P 6 P 3 P 5 t 2 P 4(2)t 3发生后 t 4t 3 P 1 t 1P 2P 6 P 3 P 5 t 2P 4(3)t 2发生后 (4)t 1不能发生t 4t 3 P 1t 1 P 2 P 6 P 3 P 5 t 2 P 4 (5)t 4发生后2161882272293245 3413413281228204 4206820291023511 5103103302558510 65186312553505 73333322528480 8168168332403355 9843331341778242 101658122351213189 1161310136948436 12508508372183135 13254349538678166 14247843039833321 15215310540160872 165281641363363 178383421818282 18418418431413389 19209345441948412 2022822845206315 211143119467878 225988647393393 23433433481968432 242168120492163115 25603915057866n n 5133333376828316 52166813277158347 5366315178238238 54758246791193169 55123320980848336 56104824811683147 5712312382738226 58618106831133109 59533218454836 6010810885183183 615433186918406 62158158872033497 63793281882488440 641408384892203155 65192338790778266 661938402911333309 67201347792154812682388340936363 69170316794318318 7083832695159357 7116339796288288 72488488971443419 73244339598209850 74197844299253253 7522131651001268244n n 1011223199126478478 1029984861272393345 10324333851281728192 1041928392129963451 10519634271302258210 106213890131105329 107453453132148148 1082268220133743231 1091103791341158134 110398398135673161 1111993457136808296 11222882401371483459 11312031791382298250 1148983861391253229 11519333971401148124 1161988452141623111 117226321514255846 118107854143233233 1192732731441168144 1201368344145723211 1211723187146105834 122938426147173173 123213385148868356 1244284281491783247 1252143951501238214n n 1511073491764848 152248248177243243 15312432191781218194 154109874179973461 1553733731802308260 1561868332181130327915716631271821398374 1586381261831873337 1596331211841688152 16060896185763251 1614834831861258234 16224183701871173149 1631853317188748236 1641588521891183159 165263263190798286 16613182941911433409 1671473449192204800 168224820019333 16910034911941818 17024584101959393 17120535196468468 17228281972343295 1731431431981478454 1747182061992273225 175103392001128104对上述数据进行参数检验如下：经计算可知，===因此可知统计量=()==()=假定显著性水平，则查表可知故可以认为：在显著性水平时，该随机数序列总体的均值和方差与均匀分布U（0,1）的均值和方差没有显著性的差异。
离散事件系统建模与仿真第二版课后习题答案
a 8L 3
x
取为奇数，则可以达到最长的周期
p2
k 2
（1）能达到； （2）不能达到； （3）当 x 取奇数时，能取得最长周期
0
（4）当
x0
取奇数时，能取得最长周期；
2.9[解答] 利用反变化法求解，求出
1
F ( x) y 的反函数得：
1 0 y 2 1 y 1 2
6y 3 x F ( y) 4 2 y 1
1 , 0 x 30 30
0 ,其 它
x〈0
x/30 ， 0≤x≤30 1 ， x〉30 乘客等待时间小于等于5min的概率为 p( x 5) p( x 0) p(15 x 20 )
则乘客等待时间大于5min的概率为
1 4 p( 5) 1 p( 5) 1 5 5
2
2 1 1 S ( n n ) 1 2
1 2
1 1 26.553 ( ) 8 8
1 2
1.261
由于 0.05，查表有 t0.05 (14) 2.145 。由于 t 1.261 2.145 所以不能拒绝均值相等的假设。
建模与仿真
作业
1、某公共汽车站按规定从上午6:40至上午8:40内每20min有一班公 共汽车到站，某个乘客不了解其调度规律，而是每天早上7:00到7:30 均匀地随机到达车站，问旅客等待公共汽车时间多于5min的概率是 多少?（10分） 解：乘客到达概率密度函数为
p( x)
0， 则分布函数F（x）=
F (0) F (20) F (15) 1 5 1 30 30 5
系统建模与仿真习题4及答案
系统建模与仿真习题四及答案1. 假设函数，通过以下两种方法求该函数的laplace变换。
（1）F= int(f(t)*exp(-s*t),0,inf)计算（2）laplace(f)解：clc;clear;syms a t sf=(2*t^2+a*sin(t));F1=int(f*exp(-s*t),'t',0,inf)F2=laplace(2*t^2+a*sin(t))结果：F1 =limit(-(2*exp(-s*t)*s^4*t^2+2*exp(-s*t)*s^2*t^2+4*exp(-s*t)*s^3*t+4*s*t*exp(-s*t)+4*exp(-s*t)*s^2+4*exp(-s*t)+a*exp(-s*t)*s^3*cos(t)+a*exp(-s*t)*s^4*sin(t)-4-4*s^2-a*s^3)/s^3/(s^2+1),t = Inf) F2 =4/s^3+1/(s^2+1)*a2. 已知系统的传递函数为：求其输入为响应输出的解析解。
解：clc;clear;syms sG=(2*s^2+1)/((s+1)*(s+2));R=1/(s^2+1);y=ilaplace(G*R)结果y = -9/5*exp(-2*t)+3/10*cos(t)-1/10*sin(t)+3/2*exp(-t)3.考虑下面给出的多变量系统：（1）试求出系统的零点与极点；（2）绘制零极点图，判断系统的稳定性及系统是否为最小相位系统。
解：clc;clear;A=[-3 1 2 1;0 -4 -2 -1;1 2 -1 1;-1 -1 1 2]; B=[1 0;0 2;0 3;1 1];C=[1 2 2 -1;2 1 -1 2];D=[0 0;0 0];G=ss(A,B,C,D)[p,z] = pzmap(G)pzmap(G)结果：p =2.2361-2.2361-3.0000 + 0.0000i-3.0000 - 0.0000iz =-3.19830.6427系统不稳定，为非最小相位系统。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
、系统、模型和仿真三者之间具有怎样的相互关系答：系统是研究的对象，模型是系统的抽象，仿真通过对模型的实验以达到研究系统的目的。
、通过因特网查阅有关蒲丰投针实验的文献资料，理解蒙特卡罗方法的基本思想及其应用的一般步骤。
答：蒲丰投针实验内容是这样的：在平面上画有一组间距为a的平行线，将一根长度为L（L<a）的针任意掷在这个平面上，求此针与平行线中任一条相交的概率。
”布丰本人证明了，这个概率是：p=2L/（πa) (π为圆周率)利用这个公式可以用概率的方法得到圆周率的近似值。
所以，蒙特卡罗方法的基本思想就是：当试验次数充分多时，某一事件出现的频率近似等于该事件发生的概率。
一般步骤：（1）构造或描述概率过程（2）以已知概率分布进行抽样（3）建立各种估计量、简述离散事件系统仿真的一般步骤。
（1）阐明问题与设定目标（2）仿真建模（3）数据采集（4）仿真模型的验证（5）仿真程序的编制与校核（6）仿真模型的运行（7）仿真输出结果的统计分析、以第二章图2-5所示的并行加工中心系统为对象，试分别画出相应的实体流图和活动循环图，并比较它们两者有何区别和练习。
（1）实体流图(2)活动循环图、以第二章中图2-5所示的并行加工中心系统为对象，建立Petri 网模型。
3214所示Petri 网模型的运行过程，并将分析结果同例3-5相比较。
、任取一整数作为种子值，采用第三题中得到的随机数发生器生成随机数序列的前200项数据，并对其统计性能进行检验。
解：由第3题可得到一个随机数发生器：a=5 b=9 c=3 m=512取种子值，生成的随机数序列前200项数据如下：n n5000032458 4ttPtPPPPtP(2)t3发生后ttPtPPPPtP(3)t2发生后(4)t1不能发生ttPtPPPPtP(5)t4发生后2 3 4334334819684322 4 21681204921631152 5 60391557866n n5 1 333333768283165 2 1668132771583475 3 663151782382385 4 7582467911931695 5 123320988483365 6 1048248116831475 7 123123827382265 8 6181068311331095 9 5332184548366 0 108108851831836 1 54331869184066 2 1581588720334976 3 7932818824884406 4 14083848922031556 5 192338797782666 6 19384029113333096 7 20134779215481262388 3963 68 40 3 36 9 1703167943183187 0 838326951593577 1 163397962882887 2 4884889714434197 3 244339598209857 4 1978442992532537 5 22131651001268244n n1 01 12231991264784781 02 99848612723933451 03 243338512817281921 04 19283921299634511 05 196342713022582101 06 213891311053291 07 4534531321481481 08 22682201337432311 09 11037913411581341 10 3983981356731611 11 19934571368082961 12 228824013714834591 13 120317913822982501 14 89838613912532291 15 193339714011481241 16 19884521416231111 17 2263215142558461 18 1078541432332331 19 27327314411681441 20 13683441457232111 21 17231871461058341 22 9384261471731731 23 2133851488683561 24 42842814917832471 25 2143951501238214n n1 51 10734917648481 52 2482481772432431 53 124321917812181941 54 1098741799734611 55 37337318023082601 56 186833218113032791 57 166312718213983741 58 63812618318733371633 111688 159 21 84 521 60 608961857632511 61 48348318612582341 62 241837018711731491 63 18533171887482361 64 15885218911831591 65 2632631907982861 66 131829419114334091 67 14734491922048 0 01 68 22482001933 31 69 100349119418181 70 245841019593931 71 2053 51964684681 72 282819723432951 73 14314319814784541 74 71820619922732251 75 1033 92001128104对上述数据进行参数检验如下：经计算可知，===因此可知统计量=()==()=假定显著性水平，则查表可知故可以认为：在显著性水平时，该随机数序列总体的均值和方差与均匀分布U（0,1）的均值和方差没有显著性的差异。
、三角分布的概率密度函数为试写出其相应的分布函数，并采用反变换法给出生成该三角分布随机变量的算法步骤。
解：根据密度函数f(x)可计算得到x的分布函数如下：采用反变换法生成该三角分布随机变量的算法步骤如下：计算其反函数令u=F(x),则其反函数x=则数列{}即为所求的指数分布的随机变量、根据第4章复习思考题第6题中得到的结论，生成标准正态分布N（0,1）的前200项数据，并根据这些数据分别绘制相应的相关图、散点图和直方图，以检验样本数据的独立性及其分布形式是否为正态分布。
解：由已知条件可生成如下的随机数所以，由上述数据可生成如下图形：的由图形可知，样本数据大体符合正态分布的形式，虽然从图形上来看，X1 分布与标准正态分布的形式有一定的差距，但其偏差应该在误差范围内，所以可以认为样本数据是独立的、正态分布。
、分析终态仿真与稳态仿真这两种仿真方式的异同。
答：相同点：都是对系统进行仿真及输出分析的方式不同点：终态仿真结果与系统初始条件有关，而稳态仿真的最终结果是不受初始条件影响的；终态仿真主要研究的是在规定时间内的系统行为，稳态仿真更侧重于对系统长期运行的稳态行为的关注。
、对图2-1所示的简单加工系统，进行独立的重复仿真10次，每次仿真运行的长度为200，初始条件为初始队长q(0)=0,且钻床设备处于空闲状态。
仿真运行的结果如下：平均等待时间平均队长试计算求解该简单加工系统平均等待时间平均队长这两个性能指标的置信度为的置信区间。
解：1-=（1）求的置信区间由已知数据可得，=-= +=+*=平均等待时间的置信度为的置信区间为（,）（2）求的置信区间由已知数据可得，=-= +=+*=平均队长的置信度为的置信区间为（,）、对图2-1所示的简单加工系统，进行独立的重复仿真运行10次，按批统计分别得到平均队长10个批次的批平均值如下表所示。
运行次数批次1 2 3 4 5 6 7 8 911 234567891试运用批平均值法计算其平均队长Q的90%的置信区间。
解：根据表中数据可得每批十个观测值的批平均值为：总的样本均值为：平均值批次1234567891而=（++…+）=置信区间为（ +）=（,）。