实验07库存系统的仿真

1、班级： 姓名：学号：实验07 库存系统的仿真实验目的：1. 了解库存系统问题；2 掌握随机事件的计算机仿真方法。实验内容：11（单位），订货周期=5（天），每天需求的单位1 某仓库的最大库存水平为=数是随机变量。数据如下表：需求概率累积概率随机数区间00.100.1001 1010.250.3511 3520.350.7036 7030.210.9171 9140.091.0092 00另外仓库的订货一般必须提前时间（即从订货到货物到达），这个提前时间也是 随机变量。数据如下表:提前时间概率累积概率随机数区间10.60.61620.30.97930.11.00设开始库存量为3个单位，并订货8个

2、单位，安排在2天内到达。每个周期的第5 天订一次货使得库存量达到11个单位。1）编写对这个库存系统前5个周期进行仿真的M文件；2）运用编写的M文件对该库存系统进行仿真，求出每天库存量的平均水平及缺 货天数所占的比例；3）检查仿真所得结果的可靠性；4）若订货单价为c=1万元/单位，需求单价为p=2万元/单位，库存费每天w=0.1 万元/单位，缺货罚款为每天0.2万元/单位。进行库存系统仿真，计算5个周期后的 总利润。2. 某自行车商店的仓库保管人员采用一种简单的库存策略。当库存量降低到 P 辆时，就向生产厂订货，订货量为 Q如果某一天需求量超过库存量，将对商店造 成损失 (利润损失、信誉损失等

3、) ；库存量太多，又造成资金积压，保管费增多。因 此希望选择一种合适的库存策略，使总管理费用最少。现给出5种策略( 以下数据单位为辆 )：(1) P=125 ， Q=150；(2) P=125 ， Q=250；(3) P=150 ， Q=250；(4) P=175 ， Q=250；(5) P=175 ， Q=300。已知初始库存S=115(辆)，从每次订货至到货需3天(延迟时间)，保管费为0.75 元/辆，缺货损失为1.8元/辆，每次需订货费为75元，需求量为099之间的均匀分 布。假定第一天没有订货。考察每天仓库的货物 ( 自行车) 变化情况，仿真 150天，比较这 5种策略，找出最优的 一

4、种。实验仪 器与软件：1. CPU主频在2GHz以上，内存在 512Mb以上的PC;2Matlab 2010a 及以上版本。实验讲评：实验成绩：评阅教师：20 年 月 日实验 07库存系统的仿真一、第 1 题的仿真1) 仿真程序function y=Fangzhen(Z,D,M,K)Rout=round(rand(Z,D)*100);Rin=round(rand(Z+1,1)*10);kucun=K;InAmount=M-K;T(9)=In(Rin(1);T(2)=kucun;for i=1:Zfor j=1:Dm=(i-1)*Z+j;T(1)=j; T(3)=Rout(i,j); T(4)=

5、Out(Rout(i,j);T(7)=-1;T(8)=-1;if m1if T0(9)=0T(9)=T0(9)-1;elseT(9)=-1;end;if T0(9)=0 kucun=kucun+InAmount; T(2)=kucun;elseT(2)=T0(5);end;end;kucun=kucun-T(4);if kucun=1 & r=11 & r=36 & r=71 & r=1 & r=7 & r=9y=2; elsey=3;end2) 仿真实验 Fangzhen(5,5,11,3) 运行结果如下：班级： 姓名：学号：3）仿真结果分析根据这10个星期的仿真，得到结论：平均每天库存为“

6、剩余库存量”列上的 数字之和除以天数，即67 / 50=1.34 （单位）而50天中有16天可能缺货。也就是说：可能缺货的天数占考察总天数的比例 =3 / 50=6%4）检查仿真所得结果的可靠性：a、当初始库存量为4时，运行结果如下：Fan gzhe n（5,5,11,4）r MATLAB 7.10.0 (R2010a) I = | 回根据这10个星期的仿真，得到结论：平均每天库存为“剩余库存量”列上的 数字之和除以天数，即69 / 50=1.38 (单位)而50天中有16天可能缺货。也就是说：可能缺货的天数占考察总天数的比例 =5 / 50=10%b、当初始库存量为7时，运行结果如下：Fan

7、 gzhe n(5,5,11,7)File Edit Debug Parallel Desktop Window Help： &扮ED:MATLABQ袒* MATLAB 7.10.0 (R2010a)Shortcuts 创 How to Add 主 Whats NewI MATLAB 7.10.0 (R2010a)File Edit Debug Parallell Desktop Window Help:Q晶掷直韵0D:MATLAB*口由Shortcuts Z1 How to Add X Whats New根据这10个星期的仿真，得到结论：平均每天库存为“剩余库存量”列上的数字之和除以天数，即

8、67 / 50=1.34 （单位）而50天中有16天可能缺货。也就是说：可能缺货的天数占考察总天数的比例 =3 / 50=6% 5）调整参数的仿真计算：订货单价c=1万元/单位，需求单价p=2万元/单位，库存费w=0.1万元/单 位，缺货罚款为每天0.2万元/单位。总利润=销售所得-成本所支成本部分： 订购费用C = cQ = 1 54 = 5(力丿元) 库存费用：Cq = wQs其中Q = 60 127.5 80 65 20 = 352.5(单位)因此，库存保管费为Cq =wQs =0.1 352.5 = 35.25(万兀) 缺货罚款费用：q = 0.2 10 = 2(万兀)而销售所得为R=

9、 pQ = 2疋 352.5 = 705(万元)故5个周期的总利润为P = R-(C0 +CQ +q) =705-(54+35.25+2)=613.75万元)二、第2题的仿真程序如下：days=30;P=125,150;Q=150,250;cost=0,0;arrivali nterval=2;storagefee=0.5;lossfee=1.6;bookfee=75;storage0=110;book nu mber=0;arrivedate=0;nr=ran d(days,1);for i=1:2storage(1)=storage0;n=rou nd(99* nr(1);sale=n;r

10、emai n=storage(1)-n;ifremai n=P(i)book nu mber=Q(i)arrivedate=4;orderfee=bookfee;elseoederfee=0;endstorage(1)=remai n;cost(i)=cost(i)+rema in *storagefee+orderfee;for j=2:daysdh=j;if abs(dh-arrivedate)=nsale=n;remain=storage(j)-n; shortagenumber=0; elsesale=storage(j);remain=0;shortagenumber=n-stora

11、ge(j);endstorage(j)=remain;if remain+booknumber=P(i);booknumber=Q(i);arrivedate=dh+arrivalinterval;orderfee=bookfee;elseorderfee=0;endcost(i)=cost(i)+remain*storagefee+shortagenumber*lossfee+orderfee; end ;mincost=min(cost);endcostmincost运行结果如下：当库存量降低到150辆时，就向生产厂订货，订货量为250辆。因此选择P=150, Q=250勺库存策略，使总管理费用最少，最少费用为2.1306 1 03元三、