例6.11 利用卡尔曼滤波估计一个未知常数 题目： 设已知一个未知常数x 的噪声观测集合，已知噪声v(n)的均值为零， 方差为 ，v(n)与x 不相关，试用卡尔曼滤波估计该常数 题目分析： 回忆Kalman 递推估计公式 由于已知x 为一常数，即不随时间n 变化，因此可以得到： 状态方程： x(n)=x(n-1) 观测方程： y(n)=x(n)+v(n) 得到A(n)=1,C(n)=1, , 将A(n)=1,代入迭代公式 得到：P(n|n-1)=P(n-1|n-1) 用P(n-1)来表示P(n|n-1)和P(n-1|n-1),这是卡尔曼增益表达式变为 从而 2v σ1??(|1)(1)(1|1)(|1)(1)(1|1)(1)()()(|1)()[()(|1)()()]???(|)(|1)()[()()(|1)](|)[()()](|1)H w H H v x n n A n x n n P n n A n P n n A n Q n K n P n n C n C n P n n C n Q n x n n x n n K n y n C n x n n P n n I K n C n P n n --=----=----+=--+=-+--=--2()v v Q n σ=()0w Q n =(|1)(1)(1|1)(1)()H w P n n A n P n n A n Q n -=----+21 ()(|1)[(|1)]v K n P n n P n n σ-=--+22(1)()[1()](1)(1)v v P n P n K n P n P n σσ-=--=-+
一、基本概念 1、数字正弦载波调制 在通信中不少信道不能直接传送基带信号，必须用基带信号对载波波形的某些参量进行控制，使得载波的这些参量随基带信号的变化而变化，即所谓数字正弦载波调制。 2、数字正弦载波调制的分类。 在二进制时, 数字正弦载波调制可以分为振幅键控（ASK）、移频键控（FSK）和移相键控（PSK）三种基本信号形式。如黑板所示。 2、高斯白噪声信道 二、实验原理 1、实验系统组成 2、实验系统结构框图
图 1 2FSK信号在高斯白噪声信道中传输模拟框图 各个模块介绍p12 3、仿真程序 x=0:15;% x表示信噪比 y=x;% y表示信号的误比特率，它的长度与x相同FrequencySeparation=24000;% BFSK调制的频率间隔等于24KHz BitRate=10000;% 信源产生信号的bit率等于10kbit/s SimulationTime=10;% 仿真时间设置为10秒SamplesPerSymbol=2;% BFSK调制信号每个符号的抽样数等于2 for i=1:length(x)% 循环执行仿真程序 SNR=x(i);% 信道的信噪比依次取中的元素 sim('project_1');% 运行仿真程序得到的误比特率保存在工作区变量BitErrorRate中 y(i)=mean(BitErrorRate); end hold off% 准备一个空白的图 semilogy(x,y);%绘制的关系曲线图，纵坐标采用对数坐标 三、实验结论
图 4 2FSK信号误比特率与信噪比的关系曲线图 系统建模与仿真（二） ——BFSK在多径瑞利衰落信道中的传输性能 一、基本概念 多径瑞利衰落信道 二、实验原理 1、实验系统组成
系统建模与仿真习题二 1. 考虑如图所示的典型反馈控制系统框图 (1)假设各个子传递函数模型为 66.031.05 .02)(232++-+=s s s s s G ，s s s G c 610)(+=，2 1)(+=s s H 分别用feedback （）函数以及G*Gc/(1+G*Gc*H)（要最小实现）方法求该系统的传递函数模型。 (2) 假设系统的受控对象模型为s e s s s G 23 )1(12 )(-+=，控制器模型为 s s s G c 32)(+=，并假设系统是单位负反馈，分别用feedback （）函数以及G*Gc/(1+G*Gc*H)（要最小实现）方法能求出该系统的传递函数模型？如果不能，请近似该模型。 2. 假定系统为： )(0001)(111000100001024269)(t u t x t x ????? ???????+????????????----= [])(2110)(t x t y = 请检查该系统是否为最小实现，如果不是最小实现，请从传递函数的角度解释该模型为何不是最小实现，并求其最小实现。 3. 双输入双输出系统的状态方程：
)(20201000)()(20224264)(75.025.075.125 .1125.15.025.025.025.125.425.25.025.1525.2)(t x t y t u t x t x ??????=????? ???????+????????????------------= （1）试将该模型输入到MATLAB 空间，并求出该模型相应的传递函数矩阵。 （2）将该状态空间模型转化为零极点增益模型，确定该系统是否为最小实现模型。如果不是，请将该模型的传递函数实现最小实现。 （3）若选择采样周期为s T 1.0=，求出离散后的状态方程模型和传递函数模型。 （4）对离散的状态空间模型进行连续变化，测试一下能否变回到原来的系统。 4. 假设系统的传递函数模型为： 222 )(2+++=s s s s G 系统状态的初始值为?? ????-21，假设系统的输入为t e t u 2)(-=。 （1）将该传递函数模型转化为状态空间模型。 （2）利用公式 ?--+=t t t A t t A d Bu e t x e t x 0 0)()()()(0)(τττ求解],0[t 的状态以及系统输出的解析解。 （3）根据上述的解析解作出s ]10,0[时间区间的状态以及系统输出曲线。 （4）采用lsim 函数方法直接作出s ]10,0[时间区间的状态以及系统输出曲线，并与（3）的结果作比较。 5. 已知矩阵 ???? ??????----=212332110A （1）取1:1.0:0=t ，利用expm(At)函数绘制求A 的状态转移矩阵，看运行的速度如何？ （2）采用以下程序绘制A 的状态转移矩阵的曲线，看运行的速度如何？ clc;clear; A=[0 1 -1;-2 -3 3;2 1 -2]; t=0:0.1:2; Nt=length(t);
物流系统仿真 期末作业 题目：Manufacturing System Planning and Scheduling 班级：物流工程131 学号：1311393003 1311393008 姓名：黎宇帆张力夫 日期：2015-09-19 成绩：
制造系统规划与调度 翻译 2.1引言 现代生产调度工具是非常强大的，提供了广阔的范围内调整工具的行为的真实过程要求的选项和参数。 然而，更多的选项的存在，它就在实践中找到的工具的最佳配置更加困难。 即专家们经常无法预测的多种可能性的影响。 测试甚至一小部分在现实中可能的配置，对实际生产过程的影响可能需要几个月的时间，可能会严重降低整体性能。 因此，这样的试验在实践中是不可行的。 优化的生产调度仿真模型比使用真正的过程更安全，更便宜,更快，更容易测试。为了在一个中等规模的制造公司充分使用先进的调度工具的优势，找到它的一个最佳的规则和参数的优化配置。 模块化仿真模型的整个业务的制造系统和生产过程中阳极氧化阶段是建立以测试不同的调度配置的影响。调度工具的配置测试和优化进行了离线使用的仿真模型。实际生产过程不受干扰，可以非常快速、低成本的找到最优配置。 2.2问题描述 位于英国的一个中型制造商，生产一系列的不同的小压铝零件和一系列大批量的其他面向消费者的产品。典型的应用包括香水的喷雾组件和哮喘患者的分配器。这是一个高度竞争的行业，成功取决于是否能实现高效率和低成本制造。所以生产调度是非常重要的。 在过去，该公司安装的软件工具可以支持生产过程中的各个区域调度。全面提高公司绩效，增加产量和减少产品的交货时间，他们计划建立自动电抗器的供应链规划服务器–总调度系统协调当地所有的业务和生产区。为了提供最好的解决方案，调度工具供应商，预优国际（https://www.doczj.com/doc/6d7585703.html,）决定使用模拟求解调度工具的优化配置。 问题是建立一个仿真工具，它将接受的到来客户订单和生产订单排序以满足这些需求。一个重要的地方是模型的生产过程本身，以确保它的主要阶段的最佳时刻加载。阳极氧化阶段是整个生产过程中特别重要的，因此，它必须是非常详细的模拟，以测试到整体订单的交货时间可以通过阳极氧化过程阶段优化减少到什么程度。 在这种情况下的研究主要目标是以下几个： （1）为了确定公司模型间的相关业务和生产过程和确定订单和交货时间， （2）在规划部门分析和优化业务流程，为了处理传入的需求和规划生产订单。 （3）测试的整体生产时间，提高灵敏度，特别是确定是否引入特定排序规则的生产订单将减少在阳极氧化处理阶段总的处理时间。
卫星轨迹 一．问题提出 设卫星在空中运行的运动方程为： 其中是k 重力系数(k=401408km3/s)。卫星轨道采用极坐标表示，通过仿真，研究发射速度对卫星轨道的影响。实验将作出卫星在地球表面（r=6400KM ，θ=0）分别以v=8KM/s,v=10KM/s,v=12KM/s 发射时，卫星绕地球运行的轨迹。 二．问题分析 1．卫星运动方程一个二阶微分方程组，应用Matlab 的常微分方程求解命令ode45求解时，首先需要将二阶微分方程组转换成一阶微分方程组。若设，则有： 2．建立极坐标如上图所示，初值分别为：卫星径向初始位置，即地球半径：y(1,1)=6400；卫星初始角度位置：y(2,1)=0；卫星初始径向线速度：y(3,1)=0；卫星初始周向角速度：y(4,1)=v/6400。 3．将上述一阶微分方程及其初值带入常微分方程求解命令ode45求解，可得到一定时间间隔的卫星的径向坐标值y(1)向量；周向角度坐标值y(2)向量；径向线速度y(3)向量；周向角速度y(4)向量。 4．通过以上步骤所求得的是极坐标下的解，若需要在直角坐标系下绘制卫星的运动轨迹，还需要进行坐标变换，将径向坐标值y(1)向量；周向角度坐标值y(2)向量通过以下方程转换为直角坐标下的横纵坐标值X,Y 。 5．卫星发射速度速度的不同将导致卫星的运动轨迹不同，实验将绘制卫星分别以v=8KM/s ，v=10KM/s ，v=12KM/s 的初速度发射的运动轨迹。 三．Matlab 程序及注释 1．主程序 v=input('请输入卫星发射速度单位Km/s ：\nv=');%卫星发射速度输入。 axis([-264007000-1000042400]);%定制图形输出坐标范围。 %为了直观表达卫星轨迹，以下语句将绘制三维地球。 [x1,y1,z1]=sphere(15);%绘制单位球。 x1=x1*6400;y1=y1*6400;???????-=+-=dt d dt dr r dt d dt d r r k dt r d θ θθ2)(2 22222θ==)2(,)1(y r y ?????????????**-=**+*-===)1(/)4()3(2)4()4()4()1()1()1()3()4()2() 3()1(y y y dt dy y y y y y k dt dy y dt dy y dt dy ???*=*=)] 2(sin[)1(Y )]2(cos[)1(X y y y y
1.信息时代认识世界（科学研究）的三种方法是：理论研究、（_实验研究_）、（__ 仿真___）。 2.根据系统状态随时间变化是连续性还是间断性的，可将系统划分为（_连续系统_）、 （__离散系统__）。 3.系统仿真中的三个基本概念是系统、（__模型_）、仿真。 4.拟对某系统进行研究，首先要对系统作出明确的描述，即确定系统各个要素：实体、 属性、活动、（__状态_）、（_事件___）。 ?阶段性知识测试 5.系统仿真有三个基本的活动，即系统建模、仿真建模和（__仿真实验__），联系这 三个活动的是系统仿真的三要素，即系统、模型和计算机（硬件和软件）。 6.系统仿真的一般步骤是：（1）调研系统，明确问题、（2）（___设立目标，收集数据 __）、（3）建立仿真模型、（4）编制程序、（5）运行模型，计算结果、（6）（_统计分析，进行决策__） ?阶段性知识测试 7.仿真软件发展经历了四个阶段（1）高级程序语言阶段；（2）仿真程序包、初级仿 真语言阶段；（3）商业化仿真语言阶段；(4) （_一体化建模与仿真环境_）阶段。 8.常用的仿真软件有Arena、Automod、MATLAB、Promodel、（__WITNESS______）、 （______FLEXSIM___）。 9.求解简单系统问题的“原始”方法是（___解析解决____），借助（___实验__）可大大 提高该方法的效率和精度。 ?阶段性知识测试 10.排队系统可简化表示为A/B/C/D/E。其中A为到达模式；B为（服务模式）、C为服 务台数量、D为系统容量；E为排队规则。 11.常见的排队规则有：先到先服务、后到后服务、优先级服务、最短处理时间优先服 务、随机服务等。请以连线方式将下列排队规则名称的中英文对照起来。 先进先出FIFO 后进先出LIFO 随机服务SIRO 最短处理时间优先SPT 优先级服务PR ?阶段性知识测试 12.模型中，习惯称实体为成分。成分可分为主动成分和被动成分。请问排队系统中的 随机到达的顾客属于（主动）成分（主动/被动）。 13.事件是改变系统状态的瞬间变化的事情。一般指活动的开始和结束。事件可分为必 然事件（主要）、条件事件（次要）、系统事件。其中（______）一般不出现在将来事件表中（FEL）。 14.活动是具有指定长度的持续时间，其开始时间是确定。排队系统主要活动有 （_______）和服务活动。 ?阶段性知识测试 15.仿真时钟表示仿真时间的变量。Witness仿真系统中仿真钟用系统变量（TIME）表 示。 仿真策略，也称仿真算法。离散事件系统适用的仿真策略有（_事件调度法_）、活动扫描法、进程交互法、三阶段法等。 16.建立输入数据模型需要4个步骤：（1）从现实系统收集数据；（2）（_确定输入数据
一曲柄滑块机构运动学仿真 1、设计任务描述 通过分析求解曲柄滑块机构动力学方程，编写matlab程序并建立Simulink 模型，由已知的连杆长度和曲柄输入角速度或角加速度求解滑块位移与时间的关系，滑块速度和时间的关系，连杆转角和时间的关系以及滑块位移和滑块速度与加速度之间的关系，从而实现运动学仿真目的。 2、系统结构简图与矢量模型 下图所示是只有一个自由度的曲柄滑块机构，连杆与长度已知。 图2-1 曲柄滑块机构简图 设每一连杆（包括固定杆件）均由一位移矢量表示，下图给出了该机构各个杆件之间的矢量关系 图2-2 曲柄滑块机构的矢量环
3．匀角速度输入时系统仿真 3.1 系统动力学方程 系统为匀角速度输入的时候，其输入为输出为；。 (1) 曲柄滑块机构闭环位移矢量方程为: (2)曲柄滑块机构的位置方程 (3)曲柄滑块机构的运动学方程 通过对位置方程进行求导，可得 由于系统的输出是与，为了便于建立A*x=B形式的矩阵，使x=[], 将运动学方程两边进行整理，得到 将上述方程的v1与w3提取出来，即可建立运动学方程的矩阵形式 3.2 M函数编写与Simulink仿真模型建立 3.2.1 滑块速度与时间的变化情况以及滑块位移与时间的变化情况 仿真的基本思路：已知输入w2与，由运动学方程求出w3和v1，再通过积分，即可求出与r1。 （1）编写Matlab函数求解运动学方程 将该机构的运动学方程的矩阵形式用M函数compv（u）来表示。 设r2=15mm，r3=55mm，r1（0）=70mm，。 其中各个零时刻的初始值可以在Simulink模型的积分器初始值里设置
M函数如下： function[x]=compv(u) %u(1)=w2 %u(2)=sita2 %u(3)=sita3 r2=15; r3=55; a=[r3*sin(u(3)) 1;-r3*cos(u(3)) 0]; b=[-r2*u(1)*sin(u(2));r2*u(1)*cos(u(2))]; x=inv(a)*b; （2）建立Simulink模型 M函数创建完毕后，根据之前的运动学方程建立Simulink模型，如下图： 图3-1 Simulink模型 同时不要忘记设置r1初始值70，如下图： 图3-2 r1初始值设置
系统建模与仿真习题三及答案 1.已知系统 )24(32)(21+++=s s s s s G 、2 103)(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 s Transfer function: s^4 + s^3 + 10 s^2 + 28 s + 6 ------------------------------ 10 s^4 + 42 s^3 + 28 s^2 + 4 s Transfer function: 20 s^2 + 34 s + 6 -------------------------------- 10 s^4 + 42 s^3 + 30 s^2 + s – 9 2.某双闭环直流电动机控制系统如图所示：
利用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:
《生产系统建模与仿真》教学大纲 （理论课程） 开课系（部）：工程学院课程编号：010396 课程类型：专业课总学时：48 学分：3 适用专业：工业工程开课学期：2014-2015学年第一学期 先修课程：概率论与数理统计、C语言程序设计、系统工程导论 一、课程简述 《生产系统建模与仿真》是面向工程实际的应用型课程，是工业工程系的主导课程之一。学生通过本课程的学习能够初步运用仿真技术来发现生产系统中的关键问题，并通过改进措施的实现，提高生产能力和生产效率。 本课程具有较强的理论性，同时具有较强的实践性和应用性，能够有效增强学生的系统仿真理论基础，提高学生对系统仿真、分析工作的适应性，培养其开发创新能力。 本课程的教学目标是培养学生的设计能力、创新能力和工程意识。课程以制造型生产企业为核心，通过理论教学和实践环节相结合，阐述了离散事件系统建模与仿真技术在生产企业分析中的基本原理和方法。其容涉及计算机仿真技术在生产系统分析中的作用和原理、仿真软件的介绍，重点介绍排队系统、库存系统、加工系统以及输入、输出数据分析。本课程的目的是要求学生通过学习、课堂教育和上机训练，能了解如何运用计算机仿真技术模拟生产系统的布置和调度管理；并熟悉和掌握计算机仿真软件的基本操作和能够实现的功能；使学生了解计算机仿真的基本步骤。 二、课程要求 （一）教学方法 1、启发式课堂讨论 针对关键知识点、典型题和难题，通过教师提问，鼓励学生回答问题或请到讲台前做题，并请其他学生评判或提出不同的答案或不同的解决方法。目的是加强学生自主学习的能力和判断能力，培养主动思考的习惯，启发学生的探索精神。 2、重视在教学中加强知识演进的逻辑规律的讲解 提高学生的逻辑思维能力，培养学生分析问题、解决问题的能力。 3、加强计算机辅助设计、分析 将Flexsim仿真软件引入教学中。应用计算机辅助设计、分析，能方便的改变系统
补充内容：模拟调制系统的MATLAB 仿真 1.抽样定理 为了用实验的手段对连续信号分析，需要先对信号进行抽样（时间上的离散化），把连续数据转变为离散数据分析。抽样（时间离散化）是模拟信号数字化的第一步。 Nyquist 抽样定律：要无失真地恢复出抽样前的信号，要求抽样频率要大于等于两倍基带信号带宽。 抽样定理建立了模拟信号和离散信号之间的关系，在Matlab 中对模拟信号的实验仿真都是通过先抽样，转变成离散信号，然后用该离散信号近似替代原来的模拟信号进行分析的。 【例1】用图形表示DSB 调制波形)4cos()2cos(t t y ππ= 及其包络线。 clf %%计算抽样时间间隔 fh=1;%%调制信号带宽(Hz) fs=100*fh;%%一般选取的抽样频率要远大于基带信号频率，即抽样时间间隔要尽可能短。 ts=1/fs; %%根据抽样时间间隔进行抽样，并计算出信号和包络 t=(0:ts:pi/2)';%抽样时间间隔要足够小，要满足抽样定理。 envelop=cos(2*pi*t);%%DSB 信号包络 y=cos(2*pi*t).*cos(4*pi*t);%已调信号 %画出已调信号包络线 plot(t,envelop,'r:','LineWidth',3); hold on plot(t,-envelop,'r:','LineWidth',3); %画出已调信号波形 plot(t,y,'b','LineWidth',3); axis([0,pi/2,-1,1])% hold off% xlabel('t'); %写出图例 【例2】用图形表示DSB 调制波形)6cos()2cos(t t y ππ= 及其包络线。 clf %%计算抽样时间间隔 fh=1;%%调制信号带宽(Hz) fs=100*fh;%抽样时间间隔要足够小，要满足抽样定理。 ts=1/fs; %%根据抽样时间间隔进行抽样
物流系统建模与仿真考前复习题 1、名词解释（5*4分） （1）系统：系统是由若干可以相互区别、相互联系而又相互作用的要素所组成，在一定的阶层结构形成中分布，在给定的环境约束下，为达到整体的目的而存在的有机集合体。 （2）物流系统模型：物流系统模型是对物流系统特征要素、有关信息和变化规律的一种抽象表达，描述了系统各要素之间的相互关系、系统与环境之间的相互作用，以反映系统的某些本质。 （3）系统仿真：应用数学模型、相应的实用模型的装置、计算机系统、部分实物的仿真系统，对某一给定系统进行数学模拟、半实物模拟、实物模拟，以便分析、设计、研究这种给定系统；或者利用这种仿真训练给定系统的专业人员。 （4）离散事件系统：指系统状态在某些随机时间点上发生离散变化的系统。离散事件动态系统，本质上属于人造系统 （4）实体：实体是描述系统的三个基本要素(实体、属性、活动)之一。在离散事件系统中的实体可分为两大类：临时实体及永久实体。在系统中只存在一段时间的实体叫临时实体。这类实体由系统外部到达系统，通过系统，最终离开系统。临时实体按一定规律不断地到达(产生)，在永久实体作用下通过系统，最后离开系统，整个系统呈现出动态过程。 （5）事件：事件就是引起系统状态发生变化的行为。从某种意义上说，这类系统是由
事件来驱动的。在一个系统中，往往有许多类事件，而事件的发生一般与某一类实体相联系，某一类事件的发生还可能会引起别的事件发生，或者是另一类事件发生的条件等，为了实现对系统中的事件进行管理，仿真模型中必须建立事件表，表中记录每一发生了的或将要发生的事件类型和发生时问，以及与该事件相联的实体的有关属性等。 （6）仿真时钟：仿真钟用于表示仿真时间的变化。离散事件动态系统的状态是在离散时间点上发生变化的，并且由于引起状态变化的事件发生时间的随机性，仿真钟的推进步长是随机的。如果两个相邻发生的事件之间系统状态不发生任何变化，则仿真钟可以跨过这些“不活动”周期。从一个事件发生时刻推进到下一事件发生时刻，仿真钟的推进呈跳跃性，推进速度具有随机性。 （7）事件调度法：仿真模型中的时间控制部件用于控制仿真钟的推进。在事件调度法中，事件表按事件发生时间先后顺序安排事件。时间控制部件始终从事件表中选择具有最早发生时问的事件记录，然后将仿真钟修改到该事件发生时刻。对每一类事件，仿真模型有相应的事件子程序。每一个事件记录包含该事件的若干个属性，其中事件类型是必不可少的，要根据事件类型调用相应的事件子程序。在事件子程序中，处理该事件发生时系统状态的变化，进行用户所需要的统计计算；如果是条件事件，则应首先进行条件测试，以确定该事件是否确能发生。如果条件不满足，则推迟或取消该事件。该事件子程序处理完后返回时问控制部件。 （8）进程交互法：一个进程包含若干个有序事件及有序活动。进程交互法采用进程描述系统，它将模型中的主动成分所发生的事件及活动按时间顺序进行组合，从而形成进程表，一个成分一旦进入进程，它将完成该进程的全部活动。 （9）连接：通过对象之间的连接定义仿真模型的流程，模型中对象之间是通过端口来
上海海洋大学试卷 姓名：学号：专业班名： 一．简述题（共40分） 1.什么是事件？在单通道排队系统中，哪两个典型事件影响系统的状态？这两个典型事件分别发生时，可能会改变系统哪些状态？（5分） 事件是指引起系统状态发生变化的行为或者事情 在单通道派对系统中的典型事件是：顾客到达和服务结束 顾客到达发生，系统可能会由闲开始变为忙，可能引起队长发生变化 服务结束，系统的状态可能有忙变为闲，可能引起队长发生变化 2.分析FMS（柔性制造系统）中的实体、状态、事件和活动。要求每一项写出2个。（8分） 实体：机床、工件 状态：空闲、加工 事件：工件到达、加工结束 活动：工件到达与工件加工开始这之间的一段事件是一个活动
3．在排队模型中，假定用链表来存放排队等待服务的顾客。链表中只有“到达时间”这样的单属性，当前CLOCK ＝10，已用空间表和可用空间表的情形见下图1，并且任何时候队列中的顾客数不会超过4位。若已知排队系统中依次发生的事件如下表1。 请根据表1中列出的事件画出CLOCK ＝15，CLOCK ＝20，CLOCK ＝25时的已用空间表和可用空间表的情形（注意：画出的图形中必须标上行号）。（8分）
4．库存系统仿真中有哪4种类型的事件？当这4种事件同时发生时，系统如何处理4种事件？（4分） 1 货物到达 2 顾客需求 3 仿真结束 4 月初清库 5．请问输入数据分析的基本步骤有哪些，并简述各个步骤的基本内容？（6分） 输入数据收集 分布的识别 参数估计 拟合度检验 6．在稳态仿真中，哪两种方法能够提高仿真结果的精度？（4分） 重复运行次数和增加运行长度
新疆财经大学实验报告 课程名称：物流管理综合实验 实验项目名称：系统建模与仿真 学号： 2013104059 姓名：曾华艳 班级：物流管理11-1 指导教师：林秋平 2014年 6月 2日
新疆财经大学实验报告
《铁路局联通营业厅排队仿真分析实验报告》 一、实验目的 (一)通过对铁路局联通营业厅运作的观察，建立计算机仿真全过程，对营业厅运作进行数据采集、建模和仿真分析，为联通营业厅提出改进和优化方案的建议。 (二)通过这次实验活动，全面了解计算机仿真技术在物流领域、生产制造领域等离散事件系统中的应用，理解仿真技术如何辅助管理人员进行决策。 (三)通过分组合作的形式，提供一种系统仿真工作中常见的团队协作方式的实践体验，培养协调工作、共同完成任务的能力。 二、系统描述 人们进入联通营业厅，首先要通过取票系统拿到自己的号，先在等待区等待叫号系统报自己的号。一共有2个服务台，2个服务台同时工作，哪个服务台叫到几号，拿这个号码的人就去哪个服务台，叫号系统按顺序叫号，2个服务台叫号不会发生重复现象。我们组决定针对铁路局联通营业厅叫号排队办理业务的过程进行研究，因此我们采集了仿真模型相关数据。记录了每位顾客到达时间、等待时间和离开时间。将收集的数据整理，录入excel中，并计算出了顾客的到达时间间隔和被服务时间，再利用flexsim建立仿真模型进行仿真分析与优化。 三、小组分工 （一）本组成员 1.组长：曾华艳 2.组员：晁芙蓉、陈磊、阿尔孜姑丽、宗泽宁、张振恒 （二）小组分工 1.调查收集数据和模型优化：全体成员 2.数据录入：晁芙蓉、张振恒、阿尔孜姑丽 3.数据处理：宗泽宁、阿尔孜姑丽 4.仿真模型建立与分析：陈磊、曾华艳 5.实验报告：曾华艳、晁芙蓉、宗泽宁 6.PPT 制作：张振恒、陈磊
控制系统仿真大作业 1、曲线拟合的Matlab实现和优化度检验 通过一个实际的例子，介绍最小二乘曲线拟合法的基本原理，对最小二乘曲线拟合法的Matlab实现方法进行研究，并给出曲线拟合Matlab实现的源程序。论述了Matlab软件在做曲线拟合时的用法，并进行曲线的拟合和相应的图像。 2、基于Matlab的液位串级控制系统 运用组态王和Matlab混合编程的方法设计了一个双容(两个水箱串联)液位串级在线控制系统，由组态王编制人机交互界面，用Matlab完成控制算法，二者通过DDE进行实时数据交换；采用串级控制策略,减小二次干扰的影响，验证其方法的有效性。 3、基于Matlab的变压器差动保护闭环仿真研究 应用Matlab建立了微机保护仿真系统,并对不同原理的变压器差动保护进行了仿真和比较.仿真系统采用积木式结构,根据微机保护的实现原理构建模块,实现保护的闭环仿真,对保护的动作过程进行分析. 4、基于MATLAB/SIMULINK的交流电机调速系统建模与仿真 根据直接转矩控制原理，利用MATLAB/SIMULINK软件构造了一个交流电机调速系统，该系统能够很好地模拟真实系统，实现高效的调速系统设计。仿真结果验证该方法的有效性。 5、基于MCGS和MATLAB的薄膜厚度控制系统仿真 以MCGS组态软件和MATLAB为平台,设计和仿真了一个薄膜厚度控制系统.MCGS完成硬件接口的设置、数据的实时采集、人机对话、以动画的方式显示控制系统的运行情况,MATLAB完成PID参数的自动整定,并利用动态数据交换(DDE)技术建立两者间的通讯.并分析其仿真结果。 6、Matlab在动态电路分析中的应用 用Matlab计算动态电路,可得到解析解和波形图.一阶电路先计算3要素,后合成解
第一章习题 1-1什么是仿真？它所遵循的基本原则是什么？ 答：仿真是建立在控制理论，相似理论，信息处理技术和计算技术等理论基础之上的，以计算机和其他专用物理效应设备为工具，利用系统模型对真实或假想的系统进行试验，并借助专家经验知识，统计数据和信息资料对试验结果进行分析和研究，进而做出决策的一门综合性的试验性科学。 它所遵循的基本原则是相似原理。 1-2在系统分析与设计中仿真法与解析法有何区别？各有什么特点？ 答：解析法就是运用已掌握的理论知识对控制系统进行理论上的分析，计算。它是一种纯物理意义上的实验分析方法，在对系统的认识过程中具有普遍意义。由于受到理论的不完善性以及对事物认识的不全面性等因素的影响，其应用往往有很大局限性。 仿真法基于相似原理，是在模型上所进行的系统性能分析与研究的实验方法。 1-3数字仿真包括那几个要素？其关系如何？ 答: 通常情况下，数字仿真实验包括三个基本要素，即实际系统，数学模型与计算机。由图可见，将实际系统抽象为数学模型，称之为一次模型化，它还涉及到系统辨识技术问题，统称为建模问题；将数学模型转化为可在计算机上运行的仿真模型，称之为二次模型化，这涉及到仿真技术问题，统称为仿真实验。 1-4为什么说模拟仿真较数字仿真精度低？其优点如何？。 答：由于受到电路元件精度的制约和容易受到外界的干扰，模拟仿真较数字仿真精度低 但模拟仿真具有如下优点： （1）描述连续的物理系统的动态过程比较自然和逼真。 （2）仿真速度极快，失真小，结果可信度高。 （3）能快速求解微分方程。模拟计算机运行时各运算器是并行工作的，模拟机的解题速度与原系统的复杂程度无关。 （4）可以灵活设置仿真试验的时间标尺，既可以进行实时仿真，也可以进
第一章系统建模与仿真概述 系统：系统是由两个以上相互区别或相互作用的单元有机的结合在起来，完成某一功能的综合体。 系统的特征：1.系统的整体性 2.系统的层次性 3.系统的相关系 4.系统的目的性 5.系统对环境的适应性系统： 模型：模型是对系统的特征要素，有关信息和变化规律的一种抽象表述、它反映 了系统某些本质属性，描述了系统各要素间的相互关系，系统与环境之间的相互 作用。 模型的意义：1.客观实体系统很难做试验，或者根本不能做实验。 2.对象问题虽然可以做试验，但是利用模型更便于理解。 3.模型易于操作，利用模型的参数变化来了解现实问题的本质和规 律更加经济方便。 系统模型的种类：抽象模型和形象模型 抽象模型：数学模型图形模型计算机模型概念模型 形象模型：模拟模型实体模型 建立模型的步骤： 1.弄清问题，掌握实际情况 2.搜集资料 3.确定因素之间的关系 4.构造建模 5.求解模型 6.检验模型的正确性 系统建模预防针的一般方法和步骤（P17） 仿真的发展趋势：建模方法面对对象仿真分布交互仿真人工智能与 计算机仿真虚拟现实仿真 Internet网上仿真 第二章商贸物流系统建模与仿真 商贸流通在社会经济中的地位与作用：1，商贸流通是连接生产和消费的纽带； 2，商贸流通对生产具有反作用； 3，商贸流通是国民经济现代化的支柱。 商贸活动的内容： 1，商流，对象物所有权转移的活动称为商流。 2，物流，是指事物从供给方向需求方的转移。
3，资金流，主要是指资金流的转移过程，包括付款，转账等过程，是 整个商贸活动的目的。 4，信息流，指商品信息的提供，商品促销信息，技术支持，售后服务 等内容，也包括诸如询单价，报单价，付款通知单，转账通知单等商业贸易单证以及交易 方的支付能力和支付信誉。 预测：所谓预测就是人们对某一不确定的或未知事件的表述。 预测的作用：从变化的事物中找出使事物发生变化的变化的固有规律，寻找和研究各种变化现象的背景及其演变的逻辑关系，从而去揭示事物未来的面貌。 判断预测方法：一，部门负责人评判预测法；二，销售人员估计法；三，德尔菲法；四， 历时类比法。 德尔菲法：依靠技术专家小组背靠背景来判断，来代替面对面的会议，是不同专家将分歧的幅度和理由都能够表达出来，经过客观分析以求达到客观规律的一致意见。 时间序列预测技术：一，移动平均预测法（计算题p30例2）； 二，指数平均预测法。 DRP：是分销需求计划的简称，它是MRP原理和技术在流通领域中的应用。该技术主要解决分销物资的应用和调度问题，其基本目标是合理进行分销物资和资源配置，以达到既有效 地满足市场需求优势的配置费用最省的目的。 *DRP的基本概念 1.库存：指仓库或物流中心实际存在的物资数量。 2.安全库存：为便于生产经营活动正常进行，防止因需求货供应的波动 引起缺货或停工待料，经常在仓库各项目保持一定数量的计划库存量， 成为安全库存。 3.期初和期末库存：指在论述的时间段开始和结束时本单位的实际库存。 4.进货提前期：指从发出订货到所定货物运回并入库所需要的时间长度。 5.送货提前期：指从接收订单到货物送到用户手中并接收入库的时间长度。 6.在途物资：指供应商已经接受订单备货，但尚未来到本单位入库的物资。 7.订货批量：指一次订货所订的物资数量。 8.时间周期：就是根据实际需要划分的时间段信息，如一日，周，月划分。 9.计划期：是指DRP尽心运算的整个时间段，可能是一个月，一个季度 或一年，他可划分为几个计划周期。 10.物流中心：从事物流活动的具有完善的信息网络的场所或组织。 BOD简介：B OD是MRP中物料清单BOM的概念和结构在分销领域的运用，它同BOM在产品结构树中连接各零件和成品一样，在供应方和各个需求方之间架起了一座沟通的桥梁。 DRP在分销网络中的运作原理（p43DRP原理图）
《系统建模与仿真》复习题样例：考试内容主要但不 限于如下内容 一、单项选择题(每题.5分,共32题) 1、下列哪个图标表示输送链Conveyor元素（C）。----序号17 A、 B、 C、 D、 2、某条生产线生产产品A，生产速率为1件/3分钟，生产的产品将送入仓库Buf 存储起来，假设生产线产出的第一件A在仿真时刻3，则运行至仿真时刻60，统计进入Buf的零件A的数量可以使用下面的函数(B)。----序号507 A、NPARTS(Buf) B、NPARTS(A) C、NPARTS2(Buf,A,1) D、APARTS(Buf) 3、一次能处理多个部件，即n个部件输入n个部件输出的是：（B ）。----序号218 A、单处理机 B、批处理机 C、装配机 D、生产机 4、在模型中有一属性元素process_time，表示不同的零件在某一机器上所需要的加工时间，那么，在机器详细设计中，对机器的加工时间cycle time栏中应输入（）----序号144 A、process_time B、process_time（） C、match D、cycle time 5、对缓冲器（buffer）中几个缓冲区用矩形框框起来的可视化设计，其所需要使用的可视化属性是（B ）。----序号134 A、name B、rectangle C、patch D、part queue 6、零件（part）到达系统的时间间隔规律在零件详细设计对话框的（ C）中进行设置。----序号148 A、type B、first arrival C、inter arrival D、to 7、有3个零件nut一次性进入系统缓冲区buf1中等待机器加工，机器加工该零部件的时间为3分钟，则计算Bmaxtime（buf1）的结果是（）。----序号261 A、3 B、6 C、8 D、9 8、可以用于机器（machine）输入（from）规则的是（ A）。----序号80 A、pull B、push C、send D、take 9、下列哪个图标表示时间序列曲线Timeseries元素（A ）。----序号21 A、 B、 C、 D、
西安理工大学 研究生课程论文/研究报告 课程名称：电力电子系统建模与分析 任课教师： 完成日期：2016 年7 月 5 日 专业：电力电子与电力传动 学号： 姓名： 同组成员： 成绩：
题目要求 某用户需要一直流电源，要求：直流输出24V/200W，输出电压波动及纹波均<1%。用户有220V交流电网（±10%波动变化）可供使用： (1) 设计电源主电路及其参数； (2) 建立电路数学模型，获得开关变换器传函模型； (3) 设计控制器参数，给出控制补偿器前和补偿后开环传递函数波特图，分 析系统的动态和稳态性能； (4) 根据设计的控制补偿器参数进行电路仿真，实现电源要求； (5) 讨论建模中忽略或近似因素对数学模型的影响，得出适应性结论(量化 性结论：如具体开关频率、具体允许扰动幅值及频率等）。 主要工作 本次设计主要负责电源主电路及其参数的的设计，以及建立电路数学模型并获得开关变换器传函模型这两部分内容，具体如下： (1) 本次设计电源主电路及其参数，采用从后向前的逆向设计思想。首先根据系统输出要求，设计了后级DC/DC型Buck电路的参数。接着设计了前级不控整流电路以及工频变压器的参数。考虑到主电路启动运行时的安全性，在主电路中加入了软启动电路； (2) 本次DC/DC变换器的建模并没有采用传统的状态空间平均方法，而是采用更为简单、直观的平均开关建模方法，建立了Buck变换器小信号交流模型。最后，推到出了开关变换器的传递函数模型，并给出了Buck电路闭环控制框图。
1 设计主电路及其参数 1.1主电路设计 根据题目要求，系统为单相交流220V/50Hz 输入，直流24V/200W 输出。对于小功率单相交流输入的场合，由于二极管不控整流电路简单，可靠性高，产生的高次谐波较少，广泛应用于不间断电源(UPS)、开关电源等场合。所以初步确定本系统主电路拓扑为：前级AC-DC 电路为电源经变压器降压后的二极管不控整流，后级DC-DC 电路为Buck 斩波电路，其中Buck 电路工作在电感电流连续模式（CCM ），前后级之间通过直流母线和直流电容连接在一起。系统主电路结构如图1-1所示。 AC 220V/50Hz L C 1 C 2R D S 图1-1 系统主电路结构图 1.2主电路参数设计 本次设计电源主电路参数，采用从后向前的逆向设计思想。先对后级DC/DC 型Buck 电路的参数进行设计，接着对前级不控整流电路以及工频变压器的参数进行设计。下面分别对后级的Buck 电路和前级经变压器降压后的不控整流电路各参数进行分析设计。 1.2.1 输出电阻计算 根据系统电路参数：220,50;24;200i o U V Hz U V P W ===，可计算： 输出电流： /200/248.33O O I P U W V A ==≈ （1-1） 负载等值电阻： /24/8.33 2.88O O R U I V A ==≈Ω （1-2）
根据质量守恒定律，血液中药物变化量等于该时刻药物进入血液速率与从血液排泄出去的速率之差，得： 由于静脉推注时输入f10=D δ(t) 得： 求解此微分方程，得： 那么，药物血药浓度为： 三、计算题 6.在标准状况下，常人进行一次有效呼吸约吸入500ml 空气，其中氧含量约为21%，二氧化碳含量为0.03%，经过一次气体交换呼出气体中氧含量变为15%，二氧化碳量占20%。 试求：呼出气体容量E V 、耗氧量2Q V 及二氧化碳产生量2CO V 解：呼出气体容量 E V =+-2O I V V 2CO V 其中耗氧量 2Q V =I ICO E ECO V F V F ..22-（其中F.为气体含量百分比） 其中吸入气体中二氧化碳量很少，在计算中可忽略不计，所以可得二氧化碳产生量为 2CO V =E ECO V V .2 由已知数据代入以上三式得： ? ?? ??=-?=+-=E CO E O CO O E V V V V V V V 2.015.050021.0500222 2 10 1011 )() (f t x k dt t dx +-=??? ??=-=+D x t x k dt t dx )0()() (11011t k De t x 01)(1-=t k e V D t C 011)(-=
可解得： ??? ??===ml V ml v ml V co o E 5.12195.136072 2 （2）（心电正问题）是研究心脏电兴奋在不同的心脏状态下是如何传播及形成体表电位的；（心电逆问题）是指从体表电位分布推断心脏内的电活动进程即求取心电源的分布。 计算题 主动脉模型中，有3个胸主动脉段内含有气囊，故在这三段的建模中，其容积下限设定为该段内气囊的瞬时体积。由于气囊的介入，在这三个胸主动脉段内产生血流等效粘滞阻力和惯性项。那么血流等效粘滞阻力和惯性项的计算公式是什么？ L n =L 0/(+) R n =R 0/[1.333r b +0.667 式中和分别为第n 段主动脉和其内气囊的半径。 L 0和R 0由下式给出 L 0=L* R 0=R* 1、建立模型一般过程为（实验设计）、（模型结构的确定）、（参数估计）、（模型验证）。 2、体温控制系统热交换系统需从热量在体内的(产生)、（传导）、（散出）过程中分析规律。 选择 1、古典生物膜理论建立基础是（A 、D ） A 、扩散 B 、布朗运动 C 、定向运动 D 、漂移 2、LFX 仿真能得到（A 、B 、D ）信息。 1、体表点位分布图 B 、12导联心电图C 、心磁D 、心脏兴奋时序图 3、半知模型称为（C ） A 、白箱 B 、黑箱 C 、灰箱 D 、透箱 以下几种算法哪一种训练神经网络收敛速度最快（ b ） A. 模拟退火算法 B.带有免疫算子的遗传算法 C.蒙特卡洛算法 D.遗传算法 2．以下那个选项不属于呼吸过程（ d ） A.外呼吸 B.气体在血液中的运输 C.内呼吸 D.琥珀酸循环 3．生物系统建模时常用四种模型是物理近似模型、物理模拟模型、图解文字或符号式模型、数学模型。 4.常用来解决非线性模式识别问题的生物系统模型是神经网络。 5.在标准状况下测得某人吸入空气后，呼出氧气16％ 氮气78％ 二氧化碳4％ 稀有气体1％ 较多的水汽。测得当时空气各成分含量氧气21％ 氮气78％ 二氧化碳0.03％ 稀有气体0.94％ 较少的水汽。假设吸入呼出气体体积不变。请计算呼吸气体交换比。 气体交换比R=[4%*(1-21%)]/(21%-21%*4%-4%)=1.96 选择： 1， 建立模型的步骤有：○ 1试验设计○2模型结构的确定○3参数估计