Lab07.库存系统的仿真.doc

Lab07库存系统的仿真实验目的： 1了解库存系统问题； 2掌握随机事件的计算机仿真方法。实验内容：1某仓库的最大库存水平为11(单位)，订货周期= 5(天)，每天需求的单位数是随机变量。数据如下表：需求概率累积概率随机数区间00.100.10011010.250.35113520.350.70367030.210.91719140.091.009200另外仓库的订货一般必须提前时间(即从订货到货物到达)，这个提前时间也是随机变量。数据如下表：提前时间概率累积概率随机数区间10.60.61620.30.97930.11.00设开始库存量为3个单位，并订货8个单位，安排在2天内到达。每个周期的第5天订一次货使得库存量达到11个单位。1) 编写对这个库存系统前5个周期进行仿真的M文件；2) 运用编写的M文件对该库存系统进行仿真，求出每天库存量的平均水平及缺货天数所占的比例；3) 检查仿真所得结果的可靠性；4) 若订货单价为c=1万元/单位，需求单价为p=2万元/单位，库存费每天w=0.1万元/单位，缺货罚款为每天0.2万元/单位。进行库存系统仿真，计算5个周期后的总利润。2. 某自行车商店的仓库保管人员采用一种简单的库存策略。当库存量降低到P辆时，就向生产厂订货，订货量为Q。如果某一天需求量超过库存量，将对商店造成损失(利润损失、信誉损失等)；库存量太多，又造成资金积压，保管费增多。因此希望选择一种合适的库存策略，使总管理费用最少。现给出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种策略，找出最优的一种。实验仪器与软件： 1CPU主频在2GHz以上，内存在512Mb以上的PC；2Matlab 2010a及以上版本。实验讲评：实验成绩：评阅教师：20 年 月 日班级 ：数应（2）班 姓名： 江美珍 学号：2010051303Lab07库存系统的仿真一、第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:Z for j=1:D m=(i-1)*Z+j; T(1)=j;T(3)=Rout(i,j);T(4)=Out(Rout(i,j); T(7)=-1; T(8)=-1; if m1 if TO(9)=0 T(9)=TO(9)-1; else T(9)=-1; end; if TO(9)=0 kucun=kucun+InAmount; T(2)=kucun; else T(2)=TO(5); end; end; kucun=kucun-T(4); if kucun=1 & r=11 & r=36 & r=71 & r=1&r=7&r m12(5,5,11,3)ans(:,:,1) = 1 3 75 3 0 0 -1 -1 1 2 0 85 3 0 3 -1 -1 0 3 5 2 0 5 0 -1 -1 -1 4 5 71 3 2 0 -1 -1 -1 5 2 68 2 0 0 11 9 1ans(:,:,2) = 1 0 45 2 0 2 -1 -1 0 2 9 53 2 7 0 -1 -1 -1 3 7 68 2 5 0 -1 -1 -1 4 5 43 2 3 0 -1 -1 -1 5 3 30 1 2 0 9 9 1ans(:,:,3) = 1 2 93 4 0 2 -1 -1 0 2 7 20 1 6 0 -1 -1 -1 3 6 38 2 4 0 -1 -1 -1 4 4 30 1 3 0 -1 -1 -1 5 3 54 2 1 0 10 6 1ans(:,:,4) = 1 1 47 2 0 1 -1 -1 0 2 9 67 2 7 0 -1 -1 -1 3 7 83 3 4 0 -1 -1 -1 4 4 19 1 3 0 -1 -1 -1 5 3 15 1 2 0 9 5 1ans(:,:,5) = 1 2 42 2 0 0 -1 -1 0 2 9 84 3 6 0 -1 -1 -1 3 6 50 2 4 0 -1 -1 -1 4 4 19 1 3 0 -1 -1 -1 5 3 70 2 1 0 10 9 13) 仿真结果分析 答：一般而言，仿真的次数越多，则结论可能更加可靠和合理一些。每天库存量的平均水平为：“剩余库存量”列上的数字之和除以总的天数。即：68/25=2.72 而25天中有8天可能缺货，也就是说：缺货天数所占的比例为：8/25*100%=32%。4) 调整参数的仿真计算 m12(5,5,11,3)ans(:,:,1) = 1 3 82 3 0 0 -1 -1 1 2 0 29 1 0 1 -1 -1 0 3 7 84 3 4 0 -1 -1 -1 4 4 44 2 2 0 -1 -1 -1 5 2 52 2 0 0 11 9 1ans(:,:,2) = 1 0 64 2 0 2 -1 -1 0 2 9 34 1 8 0 -1 -1 -1 3 8 57 2 6 0 -1 -1 -1 4 6 69 2 4 0 -1 -1 -1 5 4 88 3 1 0 10 7 1ans(:,:,3) = 1 1 82 3 0 2 -1 -1 0 2 8 53 2 6 0 -1 -1 -1 3 6 37 2 4 0 -1 -1 -1 4 4 62 2 2 0 -1 -1 -1 5 2 17 1 1 0 10 1 1ans(:,:,4) = 1 1 66 2 0 1 -1 -1 0 2 9 73 3 6 0 -1 -1 -1 3 6 70 2 4 0 -1 -1 -1 4 4 79 3 1 0 -1 -1 -1 5 1 98 4 0 3 14 0 3ans(:,:,5) = 1 0 34 1 0 4 -1 -1 2 2 0 31 1 0 5 -1 -1 1 3 0 55 2 0 7 -1 -1 0 4 7 96 4 3 0 -1 -1 -1 5 3 27 1 2 0 9 9 1总利润为：总需求量的费用-总订货费-总库存费-总缺货罚款费而总需求费为86万，总订货费为64万，总库存费为6.5万，总缺货罚款为4万，所以总利润为11.5万。二、第2题的仿真仿真程序：clear days=30;P=125,150;Q=150,250;cost=0,0;arrivalinterval=2;storagefee=0.5;lossfee=1.6;bookfee=75;storage0=110;booknumber=0;arrivedate=0;nr=rand(days,1);for i=1:2 storage(1)=storage0; n=round(99*nr(1); sale=n; remain=storage(1)-n; if remain=P(i) booknumber=Q(i) arrivedate=4; orderfee=bookfee; else oederfee=0; end storage(1)=remain; cost(i)=cost(i)+remain*storagefee+orderfee; for j=2:days dh=j; if abs(dh-arrivedate)=n sale=n; remain=storage(j)-n; shortagenumber=0; else sale=storage(j); remain=0; shortagenumber=n-storage(j); end storage(j)=remain; if remain+booknumber=P(i); booknumber=Q(i); arrivedate=dh+arrivalinterval; orderfee=bookfee; else orderfee=0; end cost(i)=cost(i)+remain*storagefee+shortagenumber*lossfee+orderfee; end; mincost=min(cost);endcostmincost运行结果为：booknumber = 150booknumber = 250cost = 1.0e+003 * 2.4115 3.3115min