系统建模与仿真实验
实验一
(1)
num=([1.3 2 2.5]);
den=([1 0.5 1.2 1]);
sys2=tf(num,den,'inputdelay',2,'inputName','flow','OutputName','Temp') (2)
clc;
num={[1 1],[1 0];1,2};
den={[1 2 1],[1 0 2];[2 1],[1 1]};
sys4=tf(num,den,0.2)
get(sys4);
(3)
num={1 1;2 [1 1] };
den={[1 1],[1 2];1,[1 2]};
sys3=tf(num,den)
(4)
num=([1.3 2 2.5]);
den=([1 0.5 1.2 1]);
sys2=tf(num,den,'inputdelay',2)
sys6=zpk(sys2)
z=sys6.z;
Z=z{ : }
p=sys6.p;
P=p{ : }
k=sys6.k
(5)A=[0 1 0;0 0 1;-5 -20 -1];
B=[0;0;1];
C=[1 0 0];
D=0;
[num,den]=ss2tf(A,B,C,D);
printsys(num,den)
[z,p,k]=ss2zp(A,B,C,D);
zpk(z,p,k)
2(有错误)
sys1=tf(10,[1 5]);
A=[-9.201 17.7791;-1.6943 3.2138];B=[-0.5112 0.5362;-0.002 -1.8470];
C=[-3.2897 2.4544;-13.5009 18.0745];D=[-0.5476 -0.1410;-0.6459 0.2958]; Q=[U Y Y;3 1 -4;4 3 0];
sys2=ss(A,B,C,D);
sys3=tf([2 2],[1 2]);
sys=append(sys1,sys2,sys3)
inputs=[1 2];
outputs=[2 3];
sysc=connect(sys,Q,inputs,outputs);
tfsys=(sysc);
clear
char u y;
sys1=tf(10,[1 5]);
A = [ -9.201 17.7791; -1.6943 3.2138 ];
B = [ -0.5112 0.5362; -0.002 -1.8470];
C = [ -3.2897 2.4544; -13.5009 18.0745];
D = [-0.5476 -0.1410;-0.6459 0.2958 ]; sys2=ss(A,B,C,D);
sys3=tf([2 2],[1 2]);
sys=append(sys1,sys2,sys3);
Q=[1 1 1;3 1 -4;4 3 0];
inputs=[1 2];
outputs=[2 3];
sysc= connect(sys,Q,inputs,outputs)
tfsys=tf(sysc)
3
A=[0 1 0;-2 -3 0;-1 1 -3];
B=[1 1 2]';
C=[1 0 1];
[Ab,Bb,Cb,Db]=canon(A,B,C,0,'modal') %产生对角标准型
[Ab,Bb,Cb,Db]=canon(A,B,C,0,'companion') %产生伴随矩阵型
4(1)
n1=30;d1=[0.5 1];
n2=[0.2,0.4];d2=[0.25,1,0];
Gkn=conv(n1,n2);Gkd=conv(d1,d2);
[num,den]=cloop(Gkn,Gkd);
P=roots(den);
disp('系统的闭环极点为'),disp(p)
ss=find(real(P)>0);tt=length(ss);
if(tt>0)
disp('系统不稳定')
disp('位于右半复平面的极点为'),disp(P(ss))
else disp('系统是稳定的')
end
(2)
num=[2,3,-1,0.6,3,2];
den=[6,4,-1,0.6,3,0.8];
[Z,P]=tf2zp(num,den);
ss=find(abs(P)>1);tt=length(ss);
if(tt>0)
disp('系统不稳定')
disp(['位于Z平面单位圆之外的极点有' int2str(tt) '个,分别为']),disp(P(ss))
else disp('系统是稳定的')
end
(3)
A=[1,2,0;-6,-2,3;-3,-4,0];
[m,n]=size(A);
if(n~=m)
disp('输入错误,系统矩阵不是方阵')
else
if(rank(A) disp('系统平衡状态不止一个') else Q=eye(size(A)); P=lyap(A,Q); for ii=1:m detp(ii)=det(P([1:ii],[1:ii])); end ss=find(detp<=0);tt=length(ss); if(tt>0) disp('系统平衡状态是不稳定') else disp('P正定,系统在原点处平衡状态是渐进稳定的') end end end 实验二 (1) syms s y=ilaplace(1*(25/(s^2+2*s+25))) num=[25]; den=[1,2,25]; impulse(num,den) grid (2) numk=25;denk=[1,4,0]; [num,den]=cloop(numk,denk); printsys(num,den) t=0:0.1:5; u1=1+0.2*sin(4*t); u2=0.3*t+0.3*sin(5*t); y1=lsim(num,den,u1,t); y2=lsim(num,den,u2,t); plot(t,[y1,y2]) grid xlabei('t(s)'); ylabei('y(t)'); gtext(0.6,1.1,'u1'); gtext(0.7,1.5,'y1'); gtext(4.5,1.3,'u2'); gtext(4.2,1.5,'y2'); (3) A=[-4,-2.5,-0.5;1.5,0,0;0,2.5,0]; B=[1;2;3];C=[0,1.5,1;-3,0,-6];D=[0;0];x0=[-1,1,-1]; [G,H,Cd,Dd]=c2dm(A,B,C,D,0.5,'zoh'); dinitial(G,H,Cd,Dd,x0) (4) num=[4.8,28.8,24]; den=[1,9,26,24]; step(num,den) grid (5) z=-0.5;p=[0,-2,1];k=6; sys=zpk(z,p,k);sys1=feedback(sys,1); t=0:0.01:10; u=t; lsim(sys1,u,t) kv=limit(s*(6*s+3)/(s^3+s^2-2*s),s,0); ess=1/kv num=1;den=[1,1]; %输入开环传递函数模型[ncloop,dcloop]=cloop(num,den), %求闭环传递函数 w=[2,3];rphase=[0,30];rmag=[2,4]; %输入正弦信号参数 G=freqresp(ncloop,dcloop,sqrt(-1)*w) %频率特性计算 Gmag=abs(G), %闭环系统幅频特性 Gphase=57.3*angle(G) %闭环系统相频特性 Cmag=Gmag'.*rmag %频率响应幅值向量 Cphase=Gphase'+rphase %频率响应初始相位向量 Css(t)=sprintf('=C1s(t)+C2s(t)') Css(t)=sprintf('=%dsin(%dt+%d)+%dcos(%dt+%d)',Cmag(1),w(1),Cphase(1),Cmag( 2),w(2),Cphase(2)) 实验三 1. num=1;den=[1,1]; %输入开环传递函数模型[ncloop,dcloop]=cloop(num,den), %求闭环传递函数 w=[2,3];rphase=[0,30];rmag=[2,4]; %输入正弦信号参数 G=freqresp(ncloop,dcloop,sqrt(-1)*w) %频率特性计算 Gmag=abs(G), %闭环系统幅频特性 Gphase=57.3*angle(G) %闭环系统相频特性 Cmag=Gmag'.*rmag %频率响应幅值向量 Cphase=Gphase'+rphase %频率响应初始相位向量 Css(t)=sprintf('=C1s(t)+C2s(t)') Css(t)=sprintf('=%dsin(%dt+%d)+%dcos(%dt+%d)',Cmag(1),w(1), Cphase(1),Cmag(2),w(2),Cphase(2)) 2. K=1;num=1; den=poly([0,-1,-2]); %由给定的极点求取传递函数的分母多项式 w=logspace(-1,1,100); %在和10频率范围内产生100个点。[m,p]=bode(K*num,den,w); subplot(2,1,1); semilogx(w,20*log10(m)); %绘制对数幅频特性曲线 grid title('对数幅频特性曲线') %加标题 xlabel('\omega(1/s)') %横轴加标题,“\omega”为特殊字符ylabel('L(\omega)(dB)') %纵轴加标题 subplot(2,1,2); semilogx(w,p); %绘制对数相频特性曲线 grid title('对数相频特性曲线') xlabel('\omega(1/s)') ylabel('\phi(\omega)(deg)') %“\phi”为特殊字符 [Kh,r,wg,wc]=margin(m,p,w) [ncloop,dcloop]=cloop(num,den) ; %求取闭环传递函数的分母多项式向量和分子多项式向量 [mc,pc,w]=bode(ncloop,dcloop); %计算闭环幅值和相角向量 subplot(2,1,1) %分割图形窗口为2×1,选中图形区域2 plot(w, mc) %绘制闭环幅频特性曲线 subplot(2,1,2) %分割图形窗口为2×1,选中图形区域2 plot(w,pc) %绘制闭环相频特性曲线 3. num=1; den=conv([1,1,0],[1,2,2]); ngrid('new') %清图后画Nichols图线nichols(num,den) %绘制系统Nichols图grid 4. num=[1,1]; den=conv([1,-1,0],[1,4,16]); sys=tf(num,den) [P, Z]=pzmap(sys) rlocus(sys) K=[0:0.05:80]; rlocus (sys, K); [Kd1,pd1]=rlocfind (sys) [Kd2,pd2]=rlocfind (sys) [Kp1,pp1]=rlocfind (sys) [Kp2,pp2]=rlocfind (sys) 5. num=1; den=[1,12,35,0]; sys=tf(num,den); rlocus(sys) clpoles=rlocus (num,den, 76.3) k=76.3/35; ess=1/k 系统建模与仿真习题二 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); Flexsim应用案例示例 示例一港口集装箱物流系统仿真 (根据:肖锋,基于Flexsim集装箱码头仿真平台关键技术研究,武汉:武汉理工大学硕士学位论文,2006改编) 1、港口集装箱物流系统概述与仿真目的 1.1港口集装箱物流系统概述 1.2港口集装箱物流系统仿真的目的 2、港口集装箱物流系统的作业流程 2.1港口集装箱物流系统描述 2.2港口集装箱物流系统作业流程 2.3港口集装箱物流系统离散模型分析 3、港口集装箱物流系统仿真模型 3.1港口集装箱物流系统布局模型设计 3.2港口集装箱物流系统设备建模 3.3港口集装箱物流系统仿真 4、仿真运行及数据分析 4.1仿真运行及数据处理 4.2仿真数据的结果分析 小结与讨论 示例二物流配送中心仿真 (根据:XXX改编) 1、物流配送中心概述与仿真目的 1.1物流配送中心简介 1.2仿真目的 2、配送中心的作业流程描述 2.1配送中心的功能 2.2配送中心的系统流程 3、配送中心的仿真模型 3.1配送中心的仿真布局模型设计 3.2配送中心的设备建模 3.3配送中心的仿真 4、仿真运行及数据分析 4.1仿真运行及数据处理 4.2仿真数据结果分析 4.3系统优化 小结与讨论 “我也来编书”示例 示例一第X章排队系统建模与仿真学习要点 1、排队系统概述 2、排队系统问题描述 3、排队系统建模 4、排队系统仿真 5、模型运行与结果分析 小结 思考题与习题(3-5题) 参考文献 1、李文锋,袁兵,张煜.2010.物流系统建模与仿真(第6章) 北京:科学出版社 2、王红卫,谢勇,王小平,祁超.2009.物流系统仿真(第6章) 北京:清华大学出版社 3、马向国,刘同娟.2012.现代物流系统建模、仿真及应用案例(第5章) 图书管理系统UML建模: 1.1、确定系统涉及的总体信息 (1)读者: ?借书 ?还书 ?书籍预定 (2)图书馆管理员: ?书籍借出处理 ?书籍归还处理 ?预定信息处理 (3)系统管理员: ?增加书目 ?删除或更新书目 ?增加书籍 ?减少书籍 ?增加读者帐户信息 ?删除或更新读者帐户信息 ?书籍信息查询 ?读者信息查询 1.2.确定系统的参与者 (1)分析系统所涉及的问题领域和系统运行的主要任务:?分析使用该系统主要功能部分的是哪些人 ?谁将需要该系统的支持以完成其工作 ?系统的管理者与维护者 (2)图书馆管理系统的参与者: ?读者(借阅者) ?图书馆管理员 ?图书馆管理系统维护者 1.3.确定系统的用例 1.3.1借阅者请求服务的用例 (1)查询借阅者信息 (2)查询书籍信息 (3)增加书目 (4)删除或更新书目 (5)增加书籍 (6)删除书籍 (7)添加借阅者帐户 (8)删除或更新借阅者帐户 1.3.2 图书馆管理员处理借书、还书等的用例 (1)处理书籍借阅 (2)处理书籍归还 (3)删除预定信息 1.3.3系统管理员进行系统维护的用例 (1)查询借阅者信息 (2)查询书籍信息 (3)增加书目 (4)删除或更新书目 (5)增加书籍 (6)删除书籍 (7)添加借阅者帐户 (8)删除或更新借阅者帐户 1.4.使用Rational Rose绘制用例图的步骤(具体详见教材P83-92) 1.创建用例图 2.用例图工具栏按钮简介 3.工具栏的定制 4.添加参与者与用例 5.添加参与者与用例之间的关系 6.添加用例之间的关系 1.5.图书馆管理系统的用例图 1.5.1借阅者请求服务的用例图 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)(_确定输入数据 第一届(2013)复杂管理系统建模与仿真国际研讨会征文通知 仿真建模、分析与优化,已经成为解决复杂管理系统的重要手段,正日益受到理论界和实践界的广泛关注。本会议的目标,旨在为国内外从事复杂管理系统建模与仿真的研究及实践人员提供一个高水平的、专业的论坛,通过思想碰撞与信息交流,探讨系统仿真领域的最新理论和实践,促进相互合作,进而推动仿真建模、分析与优化技术在复杂管理系统中的研究与应用。本会议邀请了国内外在仿真领域卓有建树的知名学者参加,如美国建模与仿真学会(SCS )主席J. Fowler 教授,美国冬季仿真会议(WSC )主任委员会委员J. Smith 教授,2012年WSC 会务主席美国亚利桑那大学Y .J. Son 教授等;以及中国系统仿真学会、广东省系统工程学会推荐的众多专家学者。会议将请各位专家学者就仿真理论、技术及其应用的最新进展做主题报告,大会还将组织针对性强的分组专题报告与研讨。第一届(2013)复杂管理系统建模与仿真国际研讨会期待您的参与和交流。 主办单位: 中国系统仿真学会离散仿真专业委员会 广东省系统工程学会 深圳大学 承办单位:深圳大学(管理学院) 赞助单位:深圳本斯集团 会议地点:中国·广东·深圳 会议网站:https://www.360docs.net/doc/ad17229131.html,/ismscs13/ 会议时间:2013年6月1日至2日 会议主席:李凤亮教授,深圳大学副校长 大会执行主席: 陈智民教授,深圳大学管理学院院长 周泓教授,北京航空航天大学 张光宇教授,广东省系统工程学会 Dr. J. Fowler ,Arizona StateUniv . (USA) Dr. J. Smith ,Auburn Univ . (USA) Dr. Y .J. Son ,Univ . of Arizona (U SA) 学术委员会主席:周泓教授,北京航空航天大学 学术委员会委员: Dr. J. Fowler ,Arizona StateUniv . (USA) Dr. J. Smith ,Auburn Univ . (USA) Dr. Y .J. Son ,Univ . of Arizona (U SA) Dr. P. Ahrweiler, Univ . of Dublin (Ireland) 范文慧教授,清华大学 何世伟教授,北京交通大学 胡斌教授,华中科技大学 隽志才教授,上海交通大学 任佩瑜教授,四川大学 卫军胡教授,西安交通大学 魏新教授,广东工业大学 徐哲教授,北京航空航天大学 徐宗昌教授,装甲兵工程学院 张光宇教授,广东工业大学 周明教授,深圳大学 朱一凡教授,国防科技大学 戈鹏副教授,四川大学 刘蕾副教授,电子科技大学 龚晓光副教授,华中科技大学 潘燕春副教授,深圳大学 赵晗萍副教授,北京师范大学 会议主题:复杂管理系统建模、仿真与分析---理论研究与应用实践。 议题范围(会议议题包括但不限于以下方面): 1. 复杂管理系统建模:基于系统科学/系统工程的方法; 2. 仿真建模与分析的理论和方法,如离散事件仿真、系 统动力学仿真、多智能体仿真、嵌入式仿真等; 3. 仿真技术与工具; 4. 基于仿真的复杂系统优化; 5. 基于仿真的风险决策与分析; 6. 面向可持续发展的绿色生产与服务:基于仿真的研 究; 7. 系统仿真与信息系统整合:决策支持与智能化管理; 8. 仿真在复杂管理系统中的应用,包括生产制造系统、 供应链与物流系统、交通运输系统、计算机/通讯网络管理系统、医疗服务管理系统、旅游与智慧景区管理系统、作战与综合保障系统,以及循环经济、项目管理、流程再造、工程或技术管理、战略管理、信息管理与电子商务等各大领域。 会议出版物:会议拟将录用的论文以光盘形式出版,并申请国际权威检索机构(EI/ISTP )审查收录。 论文投递:论文请用英文撰写,MS W ord 编辑,采用电子投稿方式,投递至大会邮箱ismscs2013@https://www.360docs.net/doc/ad17229131.html, ,论文格式规范详见会议网站https://www.360docs.net/doc/ad17229131.html,/ismscs13/。 最佳论文奖:大会将评选最佳论文奖(最多3篇),并颁发获奖证书。 重要时间:2013年3月1日,论文扩展摘要或全文(最多6页)投稿截止;2013年4月1日,论文录用与否通知;2013年5月15日,大会注册截止;2013年6月1日-2日,会议召开。 费用:版面费800元(5页内),每超一页加100元;会务费700元每位。版面费和会务费学生减半,详见会议网站。 知识点2 1. 结合具体制造系统或服务系统,分析离散事件动态系统的基本特征。 2. 什么叫“状态空间爆炸”?产生状态空间爆炸的原因是什么?它给系统性能分析带来哪些 挑战? 3. 常用的离散事件系统建模方法有哪些,它们是如何分类的? 4. 什么是马尔可夫特性?它在离散事件系统建模与分析中有什么作用? 5. 根据功能不同,仿真模型(程序)可以分为哪三个层次?分析三个层次之间的关系。 6. 分析事件调度法、活动循环法、进程交互法和消息驱动法等仿真调度方法的特点,在分 析每种调度方法基本原理的基础上,阐述几种仿真调度方法之间的区别与联系,并绘制每种仿真调度方法的流程图。 7. 结合具体的离散事件系统,如银行、理发店、餐厅、超市、医院、作业车间等,采用事 件调度法、活动循环法或进程交互法分析建立此类系统的仿真模型,试分析仿真模型中的建模元素以及仿真调度流程。 8. 从系统描述、建模要点、仿真时钟推进机制等层面,比较事件调度法、活动循环法和进 程交互法的异同之处。 9. 什么叫仿真时钟,它在系统仿真中有什么作用?什么叫仿真时钟推进机制?常用的仿真 时钟推进机制有哪些?它们的主要特点是什么,分别适合于怎样的系统? 10.结合具体的离散事件系统,分析若采用固定步长时间推进机制、下次事件时间推进机制 或混合时间推进机制时,分别具有哪些优点和缺点,以图形或文字等形式分析时钟推进流程。 11.什么叫仿真效率?什么叫仿真精度?分析影响仿真效率和仿真精度的因素? 12.从仿真效率和仿真精度的角度,分析和比较三种仿真时钟推进机制的特点,并分析三种 仿真时钟推进机制分别适合于什么样的系统? 13. 什么是蒲丰投针试验?绘制蒲丰投针试验原理图,通过推导蒲丰投针试验中针与任一直 线相交的概率,分析采用随机投针试验方法来确定圆周率π的原理。 14. 按照蒲丰投针试验的条件和要求,完成投针试验,在统计投针次数、针与直线的相交次 数的基础上,求解π的估计值,并以报表或图形等形式表达试验结果。具体要求如下: ①自行确定针的长度、直线之间的距离。 ②投针10次、20次、30次、40次、50次、…、100次、…、200次、…,分别计算针 与直线相交的概率、π的估计值。 ③以一随机变量描述上述试验结果,并通过编程或采用商品化软件,以图形、报表等形 式表示投针试验结果,分析其中的规律,并给出结论。 ④写出试验报告。 ⑤在熟悉投针试验原理的基础上,编制投针试验仿真程序,动态运行投针试验的过程。15.什么是蒙特卡洛仿真?它有什么特点,蒙特卡洛仿真应用的基本步骤是什么? 16.采用C或C++等语言,分别编写产生均匀分布、正态分布、指数分布以及威布尔分布的伪随机数序列,通过改变每种分布中参数的数值,分析不同参数数值对随机数值的影响;通过对所产生的伪随机数分布区间的统计、分析和绘图,检验伪随机数的特性及其数值特征。 17. 对于制造系统而言,库存有哪些作用和功能? 18. 在制造企业中,库存大致可以分成四种类型。简要论述四种库存的名称和功能。 19. 什么是安全库存、订货提前期?确定安全库存和订货提前期时分别需要考虑哪些因素? 20. 什么叫“订货点法”?要确定订货点,需要哪些条件?订货点法适合于怎样的库存系统? 系统建模与仿真习题三及答案 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: 生产系统建模与及仿真 实验报告 实验一Witness仿真软件认识 一、实验目的 1、学习、掌握Witness仿真软件的主要功能与使用方法; 2、学习生产系统的建模与仿真方法。 二、实验内容 学习、掌握Witness仿真软件的主要功能与使用方法 三、实验报告要求 1、写出实验目的: 2、写出简要实验步骤; 四、主要仪器、设备 1、计算机(满足Witness仿真软件的配置要求) 2、Witness工业物流仿真软件。 五、实验计划与安排 计划学时4学时 六、实验方法及步骤 实验目的: 1、对Witness的简单操作进行了解、熟悉,能够做到基本的操作,并能够进行简单的基础建模。 2、进一步了解Witness的建模与仿真过程。 实验步骤: Witness仿真软件是由英国lanner公司推出的功能强大的仿真软件系统。它可以用于离散事件系统的仿真,同时又可以用于连续流体(如液压、化工、水力)系统的仿真。目前已成功运用于国际数千家知名企业的解决方案项目,有机场设施布局 优化、机场物流规划、电气公司的流程改善、化学公司的供应链物流系统规划、工厂布局优化和分销物流系统规划等。 ◆Witness的安装与启动: ?安装环境:推荐P4 1.5G以上、内存512MB及以上、独立显卡64M以上显存,Windows98、Windows2000、Windows NT以及Windows XP的操作系统支持。 ?安装步骤:⑴将Witness2004系统光盘放入CD-ROM中,启动安装程序; ⑵选择语言(English);⑶选择Manufacturing或Service;⑷选择授权方式(如加密狗方式)。 ?启动:按一般程序启动方式就可启动Witness2004,启动过程中需要输入许可证号。 ◆Witness2004的用户界面: ?系统主界面:正常启动Witness系统后,进入的主界面如下图所示: 主界面中的标题栏、菜单栏、工具栏状态栏等的基本操作与一般可视化界面操作大体上一致。这里重点提示元素选择窗口、用户元素窗口以及系统布局区。 ?元素列表窗口:共有五项内容,分类显示模型中已经建立和可以定义的模型元素。Simulation中显示当前建立的模型中的所有元素列表;Designer中显示当前Designer Elements中的所有元素列表;System中显示系默认的特殊地点;Type中系统建模与仿真习题2
系统仿真示例
图书管理系统UML建模
系统建模与仿真考试题
中文--Call For Papers-ISMSCS-第一届复杂管理系统建模与仿真国际研讨会-
制造系统建模与仿真知识点2
系统建模与仿真习题3及答案
生产系统建模与及仿真实验报告
系统建模与仿真课后作业