计算机仿真作业.pptVIP

例: 机器修理系统仿真 有m台机器，由c个工人共同负责修理，并设： ①各台机器质量相同，机器的连续运转时间相互独立且服从同一负指数分布，平均寿命为1/v（v0）。 ②每个工人技术相同,且修理时间相互独立并服从同一负指数分布,平均修理时间为1/u(u0). ③修复后的机器其寿命分布不变。 ④机器停止运转每单位时间的损失费为c1元，工人单位时间的产值为c2元。 若机器的等待时间为E，工人总的空闲时间为F，则系统总的损失费为 S = c1E + c2F 试求当机器数m固定时，为使系统的总损失费最小，应配备多少工人为最优？ 羌神俘茹峦沉纬名钮榷裕凹墙纺糯匙惭瑚怂剁酞故冈跳简逐冕赵宦后系豆计算机仿真作业计算机仿真作业 解:采用最短时间事件步长法进行仿真，首先要确定一个初始状态。假定开始时所有的机器都正常运转，工人处于空闲状态。 设T为仿真终止时间。依次仿真计算修理工人数C分别取1、2、3、…时的系统损失费。最后根据系统损失费S的极小值来确定最优工人数。 根据以上思想并采用最短时间事件步长法仿真。 啃桓兆邯玉澎现菠歉册悍疵倾衫坊滞政便禾租邱舟董武坏辉段庄厉史狰再计算机仿真作业计算机仿真作业 ㈡ 输入原始数据c1,c2,U,V,m,T 给出m台机器的指数寿命 B(H)=1? 结束 F=F+x-K ② P= T3/T 输出C,S C=1 E=0 F=0 初始状态:B(i)=1,i=1,2,…,m;D(j)=0,j=1,2,..,c 选取寿命最短的机器i:H=i,K=A(i) 选取最早释放时刻的工人j:y=j,x=D(j) xK? ③ ① F=F+x-K 改变第y个人工的释 放时刻:D(y)=A(H) 机器由故障释放 改为正常:B(H)=1 机器由正常改为 故障:B(H)=0 C=C+1 K≥T? 选取寿命最短的 机器i:H=i,K=A(i) 给出机器修复时间 仿真时钟前进一个步长 A(H)=K-1/vLOG(RND(1)) 给出机器正常运行时间 仿真时钟前进一个步长 A(H)=K-1/vLOG(RND(1)) ㈠ ㈠ ㈡ ① ② ③ 是 否 否 是 是 否 豫译琶檄嚎城煌杆睁颤捌妓最祭乡限章态狐拂怜伙控峪痰护御刃腿限彩叫计算机仿真作业计算机仿真作业 框图中各标识符号含义如下： T:仿真终止时间 m:机器数 c:工人数 B(i):第i台机器所处状态 A(i):第i台机器现在状态的当前时刻 D(j):第j个工人修复机器的释放时刻 c1:一个工人单位时间的产值 c2:一台机器单位故障时间的损失费 { B(i)= 1 机器正常运转 0 机器发生故障 攀馒唯符坍硼次财茁茸络稀痛筐酚摇嘘践蝶盾笨兼黎糠狡梳需标仕氛肆喜计算机仿真作业计算机仿真作业 U:一个工人单位时间平均能修复的机器数 V:一台机器在单位运转时间内发生故障的平均次数 S:总损失费用 E:机器故障时间累计 F:工人空闲时间累计 K:最短时间机器的当前时刻 H:最短时间机器对应的机器号 L:工人修完机器的释放时刻 x:最短释放时间工人的当前时刻 y:最短释放时间对应的工人序号 灰崩番海庚怒沥渠携藐冶乏鸵会馅蔗怒卖挺捕弟灭个桔也陀选吉近寝炎官计算机仿真作业计算机仿真作业 下表列出了当 m= 86，1/V= 500小时，1/U= 34小时，c1=3.46元/小时，c2＝3.2元/小时时的仿真结果，其中每次预定仿真时间为一万个小时，连续仿真五次取其平均值作为仿真结果。由仿真结果可知当工人数c≤10时，最优工人数为7． 工人数C 每小时平均损失E 2 183.0 3 128.9 4 76.0 5 37.6 6 24.5 7 9.0 8 9.4 9 10.2 10 15.2 瘪墙卿鸿澈识死杨愤沿仪延新喷缕吭矿新激碑渣磷京愧甸帽娃蛙宝枚沸尖计算机仿真作业计算机仿真作业 应用举例-可靠性问题 例：一设备上三个相同的轴承，每个轴承正常工作寿命为随机变量，其概率分布如表所示． 寿命h 概率 1000 0.10 1100 0.13 1200 0.25 1300 0.13 1400 0.09 1500 0.12 1600 0.02 1700 0.06 1800 0.05 1900 0.05 在任何一个轴承损坏都可以使设备停止工作，从有轴承损坏，设备停止工作，到检修工到达开始更换部件为止，称为一个延迟时间．延迟时间也是随机变量，其概率分布如下表所示． 延迟时间min 概率 5 0.6 10 0.3 15 0.1 营秧畔漠柔属指蜂找柔缘绣侗礁彭秩绥谈贺幂豪痞纂书雨晦盂芜猾事陡庆计算机仿真作业计算机仿真作业 应用举例-可靠性问题 设备停工时每分钟损失5元,检修工每小时工时费12元，轴承每个成本 16元．更换一个轴承需要 20 min，同时更换两个轴承需要30mi