系统建模与仿真综合题目2018自动化、测控专业.doc

1、系统仿真与matlab综合试题题 目： 编 号： 姓 名 班 级 学 号 联系方式 成 绩 15系统仿真与matlab综合试题要求：1 所有程序均要求用matlab完成； 2 程序要有相应的注释，条理性好； 3 调试通过，没有错误； 4 程序每个人独立完成； 5 要有研究报告，研究报告要求有目录，封面用提供的统一封面； 6 题目从所给的题目列表中挑选，也可以自己选择，如果自己选择，应对所给题目进行详细描述，题目形式可以多样，但应体现系统仿真的基本概念及仿真策略和方法，要求用matlab实现该系统；原则上是从所给题目列表中选取； 7 特别提醒：程序每个人自己独立完成，如果发现有两个人的程序或研究

2、报告完全相同，将视同考试作弊，与考场上的考试作弊同等对待；一旦发现有从网上下载的程序，也将视同考试作弊，与考场上的考试作弊同等对待。评分标准：1 新颖性； 2 界面友好性； 3 程序可读性； 4 程序代码量； 5平时作业； 6 实验成绩。提交材料：1 研究报告纸质版（包括：封面、目录、试题建模过程、试题实现中的关键难点、程序运行指南、程序运行实例分析）；2 研究报告单独成一个word文档，命名为：姓名班级学号； 3 实验报告纸质版；4 编写的源程序电子版，于2018年12月29日前发送到； 5 全体同学提交的纸质材料，由一位（或多位）同学在2018年12月29日

3、之前交到南一楼中527室。 题目：(1) 报童问题模型报童每天清晨从报社购进报纸零售，晚上将没有卖掉的报纸退回。设报纸每份的购进价为，零售价为，退回价为，应该自然地假设为。这就是说，报童售出一份报纸赚，退回一份赔。报童每天如果购进的报纸太少，不够卖的，会少赚钱；如果购进太多，卖不完，将要赔钱。请你为报童筹划一下，他应如何确定每天购进报纸的数量，以获得最大的收入。模型假设1. 众所周知，应该根据需求量确定购进量。需求量是随机的，假定报童已经通过自己的经验或其它的渠道掌握了需求量的随机规律，即在他的销售范围内每天报纸的需求量为份的概率是。有了和，b，c，就可以建立关于购进量的优化模型了。2. 假设

4、每天购进量为份，因为需求量是随机的，可以小于，等于或大于，致使报童每天的收入也是随机的，所以作为优化模型的目标函数，不能是报童每天的收入，而应该是他长期(几个月，一年)卖报的日平均收入。从概率论大数定律的观点看，这相当于报童每天收入的期望值，即平均收入。仿真要求 根据上述假设进行系统建模与仿真，系统的输入为报纸每份的购进价，零售价，退回价，报纸需求量概率模型（报童销售范围内每天报纸的需求量为份的概率是）。系统的输出报童的平均收益。要求有输入、输出界面及仿真过程。(2) H1N1流感问题SIR模型大多数传染病如天花、流感、肝炎、麻疹等治愈后均有很强的免疫力，所以病愈的人既非健康者（易感染者），也

5、非病人（已感染者），他们已经退出传染系统。模型假设1. H1N1流感传播期内，总人数为不变，既不考虑生死，也不考虑迁移，人群分为易感染者，发病人群和退出人群 （包括死亡者和治愈者）三类，时刻内这三类人在总人数中所占比例分别为、。2. 每个病人每天有效接触的平均人数是常数，称日接触率。当病人与健康者有效接触时，使健康者受感染变为病人。根据假设，每个病人每天可使个健康者变为病人，因为病人数为所以每天共有个健康者被感染。3. 病人每天被治愈的占病人总数的比例为，称为日治愈率，治愈的病人具有了免疫力，即治愈后不再会成为二次患者。4. 、之和是一个常数1。仿真要求 根据上述假设进行系统建模与仿真，系统输

6、入为和，总人数N，日接触率，日治愈率。系统输出为时刻的健康者和病人人数N和N。要求有输入、输出界面及仿真过程。(3) 药物在体内的分布和排除药物进入机体后，在随血液输送到各器官和组织的过程中，不断地被吸收、分布、代谢，最终排出体外。药物在血液中的浓度（）称血药浓度。血液浓度的大小直接影响到药物的疗效，浓度太低不能达到预期的治疗效果，浓度太高又可能导致中毒、副作用太强或造成浪费。因此研究药物在体内吸收、分布和排除的动态过程，对于新药研制时剂量的确定、给药方案设计等药理学和临床医学的发展具有重要的指导意义和实用价值。为了研究目的，将一个机体划分成若干个房室，每个房室是机体的一部分，比如中心室和周边

7、室。在一个房室内药物呈均匀分布，而在不同的房室之间按一定规律进行转移。如果要求的精度不是太高的情况下，可以只考虑一室模型。模型假设 1. 药物进入机体后，全部进入中心室（血液较丰富的心、肺、肾等器官和组织），中心室的容积在给药过程中保持不变；2. 药物从中心室排出体外，与排除的数量相比，药物的吸收可以忽略；3. 药物排除的速率与中心室的血药浓度成正比，比例系数为；4. 给药速率为，中心室药量为。仿真要求 根据上述假设进行系统建模与仿真，系统输入为和，比例系数，中心室容积。系统输出为时刻的中心室血药浓度。要求有输入、输出界面及仿真过程。(4) 冰山运输问题在以盛产石油著称的波斯湾地区，浩瀚的沙漠

8、覆盖着大地，水资源十分贫乏，不得不采用淡化海水的办法为国民提供用水，成本大约是每立方米淡水0.1英镑。有些专家提出从相距9600km之遥的南极用拖船运送冰山到波斯弯，以取代淡水的办法。模型假设 1. 拖船航行过程中船速不变，航行不考虑天气等任何因素的影响。总航行距离；2. 冰山形状为球形，球面各点的融化速率相同，这是相当无奈的假设，在冰山上各点融化速率相同的条件下，只有球形的形状不变，体积的变化才能简单地计算。3. 记冰山球面半径融化速率为，船速为，拖船与南极的距离为。4. 假设冰山融化速率与各因素有如下关系其中，为系数，根据不同情况，可以改变。仿真要求 根据上述假设进行系统建模与仿真，系统输

9、入为，。系统输出为时刻的冰山半径。要求有输入、输出界面及仿真过程。(5) 存储问题的数学模型在我们的周边有一家配件厂，据我们得知，该厂为装配线生产若干种部件时因更换要付生产准备费（与生产数量无关），同一部件的产量大于需求时因积压资金、占用仓库要付贮存费。模型假设为了处理的方便，考虑连续模型，即设生产周期T和产量Q均为连续量。根据问题性质，我们作如下假设：1. 产品每天的需求量为常数r；2. 每次生产准备费为c1，每天每件产品贮存费为c2；3. 生产能力为无限大（相对于需求量），当贮存量降到零时，Q件产品立即生产出来供给需求；4. 假设允许缺货，每天每件产品缺货损失费为c3，但缺货数量需在下次生

10、产（或订货）时补足。仿真要求 根据上述假设进行系统建模与仿真，系统输入为需求量常数r，每次生产准备费用c1，每天每件产品储存费c2，每天每件产品缺货损失费为c3，生产周期T和产量Q。系统输出为每天平均C。要求有输入、输出界面及仿真过程。(6) 现代社会里的 “义利”问题在一个有互助和合作精神的社会中生活是很愉快和令人向往的。不过在一个追求经济利益的现代社会里，帮助别人的人会吃亏吗，而拒绝提供帮助的人反而会在社会中获益吗？不过在生活中，乐于助人的人往往具有很好的社会形象。其乐于助人的良好形象会在其所在社会传播，这样在他或她需要帮助时，周围会有很多人乐于对其提供帮助从而使其获得利益。另一方面，不喜

11、欢帮助别人的人，其自身的利益不会由于帮助别人而减少，不过这样的人由于其不良形象会使其获得别人帮助的机会大大减少。这比较符合中国人文精神中“义利统一”的价值观。现假设有一虚拟的演化的社会群体：其中的每一个个体都有一表示其富裕程度的财富值；同时每一个个体有一个据其乐于助人的程度而定的形象值，并且每一个人都有固定的根据他人的形象而决定是否对其提供帮助的策略。社会中每一个人的财富会随着帮助别人或接受别人帮助而变化。在该社会的演化过程中，每一个体根据其财富产生下一代，财富越多，产生的下一代的数目也越多；并且新生代将继承其前代的策略和财富。讨论该社会演化成互助合作社会的可能性。（注：所谓互助合作社会是指社

12、会中的大多数成员的策略是乐于助人的。）模型假设1 这一虚拟的社会群体中的所有人按其乐于助人的程度分为四类：极乐于助人型、乐于助人型、不乐于助人型和极不乐于助人型。为简便起见，记为A类、B类、C类和D类；2 施助次数的多少取决于施助方的策略、施助方的人数以及受助方的形象值；3 各类人都是以单代传递，即第（n+1）代产生的同时，第n代人消失；4 各类人的形象值与策略在社会演化过程中是不变的；5 受助与施助时，财富各按固定比例系数增加和减少，且这两个比例系数间存在一定比例关系；6. 假设每类人的财富与人数成正比例；7. 初始情况记为第0代。8. NA(n)、NB(n)、NC(n)、ND(n)代表第n

13、代人时A、B、C、D各类人的总人数；PA(n)、PB(n)、PC(n)、PD(n)代表第n代人时A、B、C、D各类人的财富值；IA、IB、IC、ID代表A、B、C、D各类人的形象值（即：被帮助的概率）；SA、SB、SC、SD代表A、B、C、D各类人的策略值（即：愿意帮助别人的概率）；a代表受助时财富增加系数(常数)；b代表施助时财富减少系数(常数)。仿真要求 根据上述假设进行系统建模与仿真，系统输入为IA、IB、IC、ID，SA、SB、SC、SD，NA(0)、NB(0)、NC(0)、ND(0) ，PA(0)、PB(0)、PC(0)、PD(0)。系统输出为NA(n)、NB(n)、NC(n)、ND

14、(n)，PA(n)、PB(n)、PC(n)、PD(n)。要求有输入、输出界面及仿真过程。(7) 飞机票的预定策略问题航空公司为某次航班发售机票，发售机票数不能太多也不能太少，若太多，乘客不能按时登机，公司不仅要付给乘客一定的赔偿费，而且乘客还将怨声载道；若太少，公司也将受到异地功能的损失。模型假设1. 乘客之间彼此独立；2. N代表飞机容量；3. g代表机票价格；4. f代表飞行费用（与乘客多少无关），f=0.6Ng；5. b代表乘客准时到达机场而未乘上飞机的赔偿费（0.1g）；6. p代表每位乘客迟到的概率p3%；7. m代表发售机票数；8. ES代表公司的平均利润；9. ES/f代表公司满

15、意度；10. P(j)代表超过j个乘客不能按时登记的概率（声誉指标），P(5)=5%；仿真要求 根据上述假设进行系统建模与仿真，系统输入为飞机容量N，机票价格g，乘客准时到达机场而未乘上飞机的赔偿费b，每位乘客迟到的概率p，发售机票数m，乘客不能按时登记人数j。系统输出为公司满意度ES/f，超过j个乘客不能按时登记的概率P(j)。要求有输入、输出界面及仿真过程。(8) 减肥模型随着人们生活水平的提高，肥胖的人越来越多，然而研究表明，体重指数（BMI：体重(kg)/身高(m)的平方）增高，一些疾病的发病率会随之上升。针对东方人的特点，联合国世界卫生组织颁布的体重指数，当18.5 BMI 24为正

16、常，24BMI 29为超重，29BMI为肥胖。为了保持身体健康，建议BMI值偏大的人合理减肥。模型假设1. 假设该人身体状况正常，且肥胖不是遗传性的；2. 体重增加正比于吸收的热量，平均每9000kcal增加体重1kg（kcal为非国际单位制单位1kcal=4.2kj）；3. 正常代谢引起的体重减少正比于体重，每公斤体重每周消耗热量一般在200kcal至320kcal之间，且因人而异；4. 运动引起的体重减少正比于体重，且与运动形式和运动时间有关，假设该人在减肥过程 中采取跑步的运动方式，并且每天运动t小时，即每周运动T小时；5. 为了安全与健康，每周体重减少不宜超过1.5kg，每周吸收热量不

17、要小于10000kcal。6. 假设每周吸收热量逐渐减少（减少量为定值），直至接近安全下限。仿真要求 根据上述假设进行系统建模与仿真，系统输入为初始体重w(0)，第0周吸收热量c(0)，热量转换系数a，正常代谢消耗系数b1，每周运动时间T，运动消耗代谢系数b2，每周吸收热量减少J。系统输出为体重w(k)。要求有输入、输出界面及仿真过程。(9) 传染病问题人群中有病人（带菌者)和健康人(易感染者)，任何两人之间的接触是随机的，当健康人与病人接触时健康人是否被感染也是随机的。模型假设1人群分为易感染者和已感染者，以下简称健康者和病人。时刻这两类人在总人数中所占的比例分别记作和。2每个病人每天有效接

18、触的平均人数是常数，称日接触率。当病人与健康者有效接触时，使健康者受感染变为病人。根据假设，每个病人每天可使个健康者变为病人，因为病人数为所以每天共有个健康者被感染。3病人每天被治愈的占病人总数的比例为，称为日治愈率。病人治愈后成为仍可被感染的健康者，显然是这种传染病的平均传染期。仿真要求 根据上述假设进行系统建模与仿真，系统输入为和，总人数N，日接触率，日治愈率。系统输出为时刻的健康者和病人人数N和N。要求有输入、输出界面及仿真过程。(10) 森林救火问题森林失火了，消防站接到报警后派多少消防队员前去救火呢?派的队员越多，森林的损失越小，但是救援的开支会越大，所以需要综合考虑森林损失费和救援

19、费与消防队员人数之间的关系，以总费用最小来决定派出队员的数目。问题分析 损失费通常正比于森林烧毁的面积，而烧毁面积与失火、灭火(指火被扑灭)的时间有关，灭火时间又取决于消防队员数目，队员越多灭火越快救援费除与消防队员人数有关外，也与灭火时间长短有关。记失火时刻为，开始救火时刻为，灭火时刻为。设在时刻森林烧毁面积为，则造成损失的森林烧毁面积为建模要对函数的形式作出合理的简单假设 研究比更为直接和方便。是单位时间烧毁面积，表示火势蔓延的程度在消防队员到达之前，即火势越来越大，即随的增加而增加；开始救火以后，即如果消防队员救火能力足够强，火势会越来越小，即应减小，并且当时 救援费可分为两部分；一部分

20、是灭火器材的消耗及消防队员的薪金等，与队员人数及灭火所用的时间均有关，另一部分是运送队员和器材等一次性支出，只与队员人数有关模型假设需要对烧毁森林的损失费、救援费及火势蔓延程度的形式作出假设。1损失费与森林烧毁面积成正比，比例系数，即烧毁单位面积的损失费2从失火到开始救火这段时间()内，火势蔓延程度与时间成正比，比例系数称火势蔓延速度。3派出消防队员x名，开始救火以后()火势蔓延速度降为，其中可视为每个队员的平均灭火速度显然应有4每个消防队员单位时间的费用为，于是每个队员的救火费用是；每个队员的一次性支出是第2条假设可作如下解释：火势以失火点为中心，以均匀速度向四周呈圆形蔓延所以蔓延的半径与时

21、间成正比，又因为烧毁面积与成正比，故与成正比，从而与成正比。仿真要求 系统输入为派出消防队员x名，每个队员的平均灭火速度，火势蔓延速度，开始救火时刻，烧毁单位面积的损失费，每个消防队员单位时间的费用，每个队员的一次性支出。系统输出为损失费和救援费以及总费用，灭火时刻，森林烧毁面积。要求有输入、输出界面及仿真过程。(11) 香烟过滤嘴问题香烟制造商既要满足瘾君子的需要，又要顺应减少吸烟危害的潮流还要获取丰厚的利润于是普遍地在香烟上安装了过滤嘴过滤嘴的作用倒底有多大，与使用的材料和过滤嘴的长度有什么关系，要从定量的角度回答这些问题就要建立一个描述吸烟过程的数学模型吸烟时毒物吸入人体的过程大致是这样

22、的：毒物基本上均匀地分布在烟草中吸烟时点燃处的烟草大部分化为烟雾，毒物由烟雾携带着一部分直接进入空中，另一部分沿香烟穿行在穿行过程中又部分地被未点燃的烟草和过滤嘴吸收而沉积下来，剩下的进入人体被烟草吸收而沉积下来的那一部分毒物，当香烟燃烧到那里的时候又通过烟雾部分进入空气，部分沿香烟穿行，这个过程一直继续到香烟燃烧至过滤嘴处为止于是我们看到，原来分布在烟草中的毒物除了进入空气和被过滤嘴吸收的一部分外，剩下的全都被人体吸入。实际的吸烟过程非常复杂并且因人而异为了能建立一个初步的模型，可以认为毒物随烟雾进入空气和沿香烟穿行的数量比例、烟雾穿行的速度、过滤嘴和烟草对毒物的吸收率等在吸烟过程中都是常数

23、模型假设1烟草和过滤嘴的长度分别是和香烟总长毒物M(毫克)均匀分布在烟草中密度为2点燃处毒物随烟雾进入空气和沿香烟穿行的数量比例是，3未点燃的烟草和过滤嘴对随烟雾穿行的毒物的吸附率（单位时间内毒物被吸收的比例)分别是常数和4烟雾沿香烟穿行的速度是常数v，香烟燃烧速度是常数u，且5将一支烟吸完后毒物进入人体的总量(不考虑从空气的烟雾中吸入的)记作Q。仿真要求 系统输入为烟草和过滤嘴的长度和，毒物总质量M，点燃处毒物随烟雾进入空气和沿香烟穿行的数量比例，点燃的烟草和过滤嘴对随烟雾穿行的毒物的吸附率和，烟雾沿香烟穿行的速度v，香烟燃烧速度u系统输出为毒物进入人体的量Q。要求有输入、输出界面及仿真过程

24、。(12) 弱肉强食问题Volterra模型处于同一自然环境中两个种群之间的关系除了相互竞争和相互依存之外，还有一种更为有趣的生存方式：种群甲靠丰富的天然资源生长，而种群乙则靠掠食甲为生地中海里的食用鱼与鲨鱼，加拿大森林中的美洲兔与山猫，阿尔卑斯山中的落叶松与芽虫等都是这种生存方式的典型生态学上种群甲称为食饵 (Prey)，种群乙称为捕食者(Predator)，二者共处组成食饵一捕食者系统(简称PP系统)模型假设食饵和捕食者在时刻的数量分别记作和，因为大海中资源丰富可以假设如果食饵独立生存则将以增长率按指数规律增长，即有捕食者的存在使食饵的增长率降低，设降低的程度与捕食者数量成正比于是满足方程

25、比例系数，反映捕食者掠取食饵的能力。 捕食者离开食饵无法生存，若设它独自存在时死亡率为，即，而食饵为它提供食物的作用相当于使死亡率降低、或使之增长设这个作用与食饵数量成正比，于是满足比例系数，反映食饵对捕食者的供养能力上述方程是在没有人工捕获情况下自然环境中食饵与捕食者之间的制约关系，是Volterra提出的最简单的模型可以看出这个模型没有考虑自身的阻滞作用仿真要求 系统输入为仿真时间T，食饵和捕食者在初始时刻的数量和，食饵独立生存增长率，捕食者掠取食饵的能力，捕食者独自存在时死亡率，食饵对捕食者的供养能力，系统输出为食饵和捕食者在时刻的数量和。要求有输入、输出界面及仿真过程。(13) 独家销

26、售的广告问题在当今信息社会中，广告在商品推销中起着极其重要的作用。当生产者生产出一批产品后，下一步便去思考如何更快更多地售出产品。由于广告的大众性和快捷性，其在促销活动中大受经营者的青睐。当然，经营者在利用广告这一手段时自然要关心：广告与促销到底有何关系，广告在不同时期的效果如何？模型假设1、商品的销售速度会因做广告而增加，但当商品在市场上趋于饱和时，销售速度将趋于极限值，这时销售速度将开始下降；2、自然衰减是销售速度的一种特性，商品销售速度的变化率随商品的销售率的增加而减少；3、设为时刻商品的销售速度，M表示销售速度的上限，为衰减因子常数，即广告作用随时间的增加而自然衰减的速度，为时刻的广告

27、水平（以费用表示）（提示：根据上面的假设，我们建立模型，其中P为影响系数，即对的影响力，P为常数。又由假设1，当销售进行到某个时刻时，无论怎样做广告，都无法阻止销售速度的下降，故选择如下广告策略：，其中A为常数）仿真要求 系统输入初始时刻销售速度，影响系数P，销售速度的上限M，广告策略，衰减因子。系统输出为时刻商品的销售速度要求有输入、输出界面，仿真过程。(14) M/M/1 排队系统（单服务员排队系统）的仿真 对于排队服务系统，顾客往往注重排队顾客是否太多、等待时间是否太长，而服务员则关心她的空闲时间。因此队长、等待时间以及服务利用率等指标可以衡量系统性能。已知顾客到达时间和服务时间的统计规

28、律（往往来自实际数据或一定的概率分布）的情况下，如何仿真排队系统。首先，进行系统的假设：（1）顾客源是无穷的；（2）排队长度没有限制；（3）到达系统的顾客按先后顺序进入服务。按照顾客到达的时间概率分布为泊松分布，顾客服务时间的长短服从负指数分布，试完成M/M/1排队系统的仿真。系统输入为：泊松分布和负指数分布中的参数，系统输出是：平均等待时间、平均队长、服务利用率。要求有输入、输出界面、顾客到达和离开的仿真过程表示。(15) M/M/N 排队系统（多服务员排队系统）的仿真多服务员排队系统在仿真上较单服务员排队系统要复杂的多，在此先对该排队系统作一些必要的假设：（1）顾客源是无穷的；（2）排队长

29、度没有限制；（3）到达系统的顾客按先到先服务原则依次进入服务；（4）服务员在仿真过程中没有休假；（5）顾客到达时排成一队，当有服务台空闲时进入服务状态；按照顾客到达的时间概率分布为泊松分布，顾客服务时间的长短服从负指数分布，试完成M/M/1排队系统的仿真。系统输入为：泊松分布和负指数分布中的参数，服务台个数，系统输出是：平均等待时间、平均队长、服务利用率。要求有输入、输出界面、顾客到达和离开的仿真过程表示。(16) 电梯群控系统的仿真模型假设：假设顾客的到达服从泊松流分布，顾客乘坐电梯是上行还是下行是随机的，顾客到达的层数也是随机的，分布可以参照相关内容假设。试设计一个仿真系统，模拟电梯的运行

30、情况，输入参数有：泊松分布的参数、楼房的层数、楼房每层的高度、电梯的最高速度和加速度、电梯台数、电梯的容量、开门和关门时间长度、平均每位乘客上梯时间、平均每位乘客下梯时间、仿真时间长度；输出参数：乘客平均等待时间、最大乘客侯梯时间、长时间侯梯率（侯梯时间超过60S的人数占总人数的比例）。要求有电梯运行的动画界面，乘客的出现可以自动产生，也可以由键盘产生，乘客到达的层数可以系统随机分配，也可以从键盘输入产生。(17) 电影院售票与咨询系统仿真模型假设：电影院售票窗口假设只有一名服务员，负责现场售票与咨询、以及电话咨询，其中现场售票与咨询的优先级高于电话咨询的优先级。假设现场售票人员的到达、现场咨询人员的到达、电话咨询的到达都服从泊松分布，每种服务时间的长短服从负指数分布。模型输入：现场售票人员的到达率、现场咨询人员的到达率、电话咨询的到达率，每种服务时间的长短的参数。模型输出：每类服务的平均等待时间、服务员在工作时间内的占用率。要求有输入输出界面、人员的动态显示过程。每种服务的产生可以由系统自动产生，也可以定义键盘上的三个按键来分别表示三种服务的产生，也可以定义鼠标的三个按键来分别表示三种服务的产生。(18) 在可视窗口下，模拟停车场车位信息的显示模型假设：每一个车位被占用状态是随机的，服从泊松分布，车的停泊时间也是随机的，服从指数分布。假定停车场有8个