数学建模论文两种随机存贮管理模型的建立和求解

1、两种随机存贮管理模型的建立和求解 摘 要：本文建立了仓库容量有限条件下单品种、多品种的允许缺货随机存贮模型。采用连续的时间变量更合理地描述了问题，简化了模型的建立。模型的求解是一个以分段的平均损失费用函数作为目标的带约束最优化问题。针对题目中的具体数据对随机量送货滞后时间的密度函数进行了估计，解出了单品种、多品种条件下最优订货点的值和存贮方案。通过分情况讨论把单品种存贮模型推广为多品种（种）存贮模型，论证了目标函数的独立变量为个,使模型更加清晰、求解方便。类比控制论中的相关理论提出了一定条件下多品种存贮的最优性原理，给出了证明，指出该原理简化模型和验证模型求解结果的作用。讨论了销售速率具有随机

2、性时的存贮模型，实际当中调整修正订货点的方法，以及仓库最大存贮量的一种预测办法。最后指出了模型的优缺点。0问题重述工厂生产需定期地定购各种原料，商家销售要成批地购进各种商品。无论是原料或商品，都有一个怎样存贮的问题。存得少了无法满足需求，影响利润；存得太多，存贮费用就高。因此说存贮管理是降低成本、提高经济效益的有效途径和方法。问题1 某商场销售的某种商品。市场上这种商品的销售速率假设是不变的，记为；每次进货的订货费为常数与商品的数量和品种无关；使用自己的仓库存贮商品时，单位商品每天的存贮费用记为，由于自己的仓库容量有限，超出时需要使用租借的仓库存贮商品，单位商品每天的存贮费用记为，且；允许商品

3、缺货，但因缺货而减少销售要造成损失，单位商品的损失记为；每次订货，设货物在天后到达，交货时间是随机的；自己的仓库用于存贮该商品的最大容量为，每次到货后使这种商品的存贮量补充到固定值为止，且；在销售过程中每当存贮量降到时即开始订货。 请你给出求使总损失费用达到最低的订货点（最优订货点）的数学模型。问题2 现给出来自某个大型超市的关于三种商品的真实数据，按你的模型分别计算出这三种商品各自相应的最优订货点。问题3 问题1是只有一种商品需要订货的情形。实际上常遇到在库存容量有限的情况下，有多种商品需要同时订货的情形，这时需考虑充分利用存贮体积的问题。设有种商品需要订货，它们每次一同从一个供应站订货，每

4、次进货的订货费为常数与商品的数量和品种无关；订购的货物同时到达，到货天数如问题1所述是随机的。这种商品的销售速率分别为（袋或盒/天），每袋（或盒）的体积分别为。使用自己的仓库和租借的仓库时单位体积商品每天的存贮费分别记成和，单位体积商品每天的缺货损失记成，自己的仓库用于存贮这种商品的总体积容量为，每次到货后这种商品的存贮量总体积补充到固定体积容量为止，且。每当这种商品的存贮量总体积降到时即开始订货。试通过建立数学模型说明应如何确定最优订货点和自己的仓库用于存贮这种商品的各自体积容量以及在订货到达时使这种商品各自存贮量补充到的固定体积 ，才能使总损失费用达到最低？问题4 如果把问题2中的三种商品

5、按问题3的方法同时订货，其中立方米，立方米，立方米，自己的仓库用于存贮这3种商品的总体积容量立方米，每次到货后这3种商品的存贮量总体积补充到固定体积容量立方米为止，且该供应站从接到订货通知到货物送达商场的天数服从在1天到3天之间的均匀分布。其余数据同问题2中相应的商品中所列出的数据。试按问题3的模型求出这3种商品的最优订货点和自己的仓库用于存贮这3种商品的各自体积容量以及在订货到达时使这3种商品各自存贮量补充到的固定体积 。问题5商品的销售经常是随机的、订货情况在一段时间后是会发生变化的，相应地商家就应该调整订货和存贮策略。你们能否对此建立数学模型加以讨论。1问题假设1不考虑商品销售率的变化；

6、2当订购货物到达时，可无限量瞬时补充，直至，即不考虑供给方的供给能力限制；3时间是连续变化的；4. 不考虑从仓库到超市的时间延迟和运输费用，即认为仓库和超市是一体的；5在多品种存贮问题中不考虑仓库之间的动态调配；2 单存贮随机问题2.1符号约定商品的销售速率：（）；每次进货的订货费：（元）；使用自己仓库存贮时，单位商品每天的存贮费：（）；租借仓库存贮时，单位商品每天的存贮费：（）；缺货时单位商品的损失为：（）每次订货交货时间为：（天）；每次交货时间为的概率密度：；自己仓库最大容量：（盒）；商品存贮量达到的固定值：（盒）；订货点：（盒）2.2单品种存贮问题分析问题1、2属于运筹学中存贮论一支，是

7、一个仓库容量有限、单品种随机存贮的最优化问题。可以通过建立起目标规划模型进行求解。它以总损失费用最低作为目标，订货点为要其中的变量。总损失费用由三个方面组成：订货费用、仓库贮存费、缺货费。订货费用每次都是固定的，即不考虑订货费用与订货数量、品种的关系。由于自己仓库容量有限，所以要租借外面的仓库存贮。存贮费用包括两个方面：使用自己仓库的存贮费和租借仓库的存贮费。而租借仓库存贮费用要高于自己仓库费用，即，所以在销售时应该先售出外仓库内存贮部分，再售出自己仓库存贮部分。给定一个订货点，当存贮量降到订货点时要发出订货单。题中订货不是立即送到，而是需要一定时间的，在存贮论中将这段时间称为滞后期，在后面的

8、论文当中都将采用该名词来描述订货送到时间。滞后期的取值是随机的，因此不能保证在订货到达时商品一定还有剩余，而有可能出现订货未到时商品已全部售出，从而发生商品短缺现象。而商品的短缺直接造成了商家的损失，这一损失有因为减少销售带来的损失、商家信誉受损、客户减少等。题中给出了缺货时单位商品每天的损失费用，说明当天缺少的商品其影响会延续到后面的时间中，直到得到补偿。在此只考虑因为减少销售带来的损失，缺货量增加的速率为销售速率。题中指出，每次订货到达后将商品存贮量补充到固定值为止。所以对于滞后期内的缺货量，在补充货物时不考虑进行补偿。滞后期的不确定，还使得订货到达时商品的剩余量是不确定的，即在订货时商家

9、不能给出具体的订货量，而要在送货到达时保证商品能够补充到固定值，则送货厂商必须要具有一定的供给能力，在题中将不考虑厂商供给能力的限制，而认为它是无限大的，商品必能补充至。因为每次订、送货情况不一样，交货时间是随机的，使得唯一的订货点不能保证每次的损失都是最低的，故采用一个概率平均值来描述损失费用。在概率平均的情况下，以连续两次收到补充订货的时间间隔作为一个时间周期。则应取损失费=，即用单位时间（题中取天）内的平均损失费用来进行评价。利用该损失值最小的约束条件来得出最优订货点。现先讨论一个周期内商品的总损失。商品的总损失由三部分组成：1） 每次进货时的订货费用；2） 商品在未卖出时需要存贮而由此

10、产生的存贮费用；3） 因缺货减少了销售量，由此造成的损失费。则总费用表示为，其中、分别表示上述各项的费用。由题知，每次进货的订货量为常数，即是固定的，它不受送货时间随机性的影响。存货费用由两部分组成：租借仓库存贮费+自己仓库存贮费，表示为。滞后期的改变会使得一个周期的长度改变，即商品存贮时间发生变化，这使的值是不固定的。是滞后期变化的。缺货损失费用是由滞后期的随机性引起的，它也是的函数。故要使一周内的总损失最小，只需考虑存货费用和缺货损失费的影响。下面将对它们进行具体的讨论分析，建立起总损失费用的具体模型。2.3单品种存贮模型的建立在问题2中给出的滞后期是一些离散数据，由之得到的滞后期随机分布

11、函数也是一些离散的值，但是考虑到实际情况缺货半天和缺货一天的损失肯定是不同的。以及我们已经假设仓库与超市是一体的，那么超市的货物是随时可以补充的，所以我们在建立模型的时候，将时间看作是连续的，采用连续的方法来进行分析建模，将库存量、缺货量、库存费用、缺货费用、滞后期分布函数等在时间上进行连续化。这样对模型进行分析和求解也比较简单。由于自己仓库容量的限制，需要租借仓库来进行存贮，这与只存在就增加了问题的经分析，根据订货点的取值将问题分为两种情况：；。1，订货点大于自己仓库的最大容量这又可以分为两种情况，分别如图2-1中左、右所示。图2-1 时库存量及日损失费随时间的变化曲线I：，即不会发生缺货现

12、象，在送货到达时商品还有剩余或恰好售完。这一情况如图2-1中左边部分所示。看左上图，以每个周期开始时刻为0点，此时总的库存量为，它以不变的销售速率均匀减少。到时刻降为，也即租借仓库内的商品销售完全。时刻库存降至，此时发出订货单。在时刻，商品仍有库存或恰好售完，而补充货物送到，将其补充至，故也即下一周期的起点。段长度即为滞后时间。定义一个日损失费，它是指单位时间（天）内除订货费之外其他损失费之和。这一值是随着时间改变的。在一个周期的时间上对它进行积分，可得到一个周期存贮和缺货损失费之和。在时段内，单位时间内自己仓库内的库存不变，其单位时间内的损失费用不变，为；而租借仓库内的商品以速率减少，其库存

13、费用的减少速率为，在左下图中即表现为AB段的斜率为。在A即周期起点处，总库存为，其中部分存在自己仓库内，部分存在租借仓库内，。在时段，库存只剩下自己仓库里的，日库存费用从以速率减少，至时刻库存降为。而点为该周期的终点，此时送货到达，而。很直观地，段的累积损失费为梯形A0B的面积，段累积损失费用为梯形BC的面积，段的累积损失费用为梯形CD的面积。计算如下：则在无缺货情况下， 一个周期内总的损失费用是面积、之和加上订货费，即：（2-1）II：，要发生缺货的情况。、段与I的情况相同，时刻库存降至0，而订货尚未送到，此后时段都处于缺货状态，单位时间的缺货费用（即缺货费用的变化率）是由0开始按的速率增加

14、。在时段，右上图库存量变化折线延伸到了0点以下，但它并不表示库存量为负，而是为了体现出缺货的状态。到时刻订货送达，缺货状态解除，而此时缺货费用增长到最高，。三角形DE的面积即该周期内总的缺货损失费。同上可计算得出：故在有缺货时，一个周期内总的损失费为面积、之和加上订货费。（2-2）观察发现，I、II情况中求得的总损失费用表达式的右端存在着相同的部分，即其是相等的，表示租借仓库存贮量至销售完时所积累存贮费用与自己仓库存贮量降至时积累存贮费之和。不同的只是滞后时间部分，该部分开始时间是存贮量降为的时刻，结束时间为送货到达时刻。对于情况I，不包含缺货费用，而II中多了缺货费用一项。故可以把式子写成：

15、 （2-3）其中，。在上面的式子当中，费用是送货时间的函数，即在每一个周期里，如果确定了，就能确定。而如题中所述，是一个随机变量，对于损失费用的评价最好使用一个概率平均值。在的密度函数已知情况下，可得到损失费用为 （2-4）2，订货点大于自己仓库的最大容量图2-2 时库存量及日损失费随时间的变化曲线如图2-2所示，此时还要分三种情况考虑：I：，送货到达时租借仓库内仍然存有商品，无缺货；II：，送货到达时租借仓库内已无存贮商品，而自己仓库内还存有部分商品，无缺货；III：，送货到达时租借仓库和自己仓库内均空，有缺货现象。下面分别加以推导这三种情况下总损失费用表达式。I：由于分析过程与时相同，故以

16、下只给出计算表达式，而不给出推导过程。；：总的损失费用为： （2-5）II：；： （2-6）III：；： （2-7）观察三种情况下的总费用的表达式，发现它们有共同的一项，其相同。对总费用求其期望值得到如下： （2-8）其中，至此，一个周期内总损失费用的数学模型建立起来。要对总损失进行评价，应该取单位时间内的平均总损失费用，故有，表示单位时间内总损失费用的概率平均值，为上面所求一个周期内的总损失费用概率平均值，而周期取的是概率平均值，此处表示送货时间的数值期望。订货点在之下和订货点在之上有不同的表达式。表达式中除，和外所有其他量均为已知定值。为一积分变量，若、确定，则通过积分可消去。而的概率密度

17、在具体的题目中可由多次订送货相关资料求得，在此视为已知量。这就说明，表达式中只有是变量，它的取值决定了总损失费用的大小。确定了，总损失费用就可以确定。得到一个目标函数：（2-9）在此目标函数中具有积分项，它是一个非线性目标优化模型。通过求解该优化模型可得出的最优解即最低订货点，它使得总损失费用最小。至此，问题1得到解答。2.4单品种存贮问题的求解与结果1概率密度函数的确定：题中指出为随机变量，没有具体给出其变化规律，但在问题2中各个具体商品却提供了一些滞后期的数据。通过这些数据，可以确定出的分布函数，即其概率密度函数。1）商品一：康师傅精装巧碗香菇炖鸡面统计连续的36次订货后到达时间天数纪录，

18、得到表2-1。表2-1 商品一交货时间X在不同时间点的概率分布时间/天01234567出现次数245155311概率5.56%11.11%13.89%41.67%13.89%8.33%2.78%2.78%采用SAS软件对进行正态分布的假设检验，其结果如图2-3和图2-4所示。由QQ图可见各点基本在直线附近，说明分布的正态性较好。利用统计学中的参数估计得知该样本符合均值为2.9722，方差为1.521的正态分布。图2-3 商品一交货时间的概率密度图图2-4 验证商品一交货时间是否正态分布的QQ图2）商品二：心相印手帕纸统计连续的43次订货后到达时间天数纪录，得到表2-2。表2-2商品二交货时间x

19、在不同时间点的概率分布时间天数12345出现次数2231251概率4.65%53.49%27.91%11.63%2.33%由QQ图可知商品二的交货时间服从威布尔分布。其密度分布函数，由最大似然估计确定Welbull分布参数a=0.0380 ，b=3.1434。图2-5 商品二交货时间的概率密度图图2-6 验证商品二交货时间是否威布尔分布的QQ图3）商品三：中汇香米5KG装表2-3商品三交货时间x在不同时间点的概率分布时间天数123456出现次数27208321概率44.26%32.79%27.91%13.11%3.28%1.64%采用与（2）相同的方法得到图2-7和2-8，图2-8QQ图可知商

20、品三的交货时间也服从威布尔分布，其分布参数为a= 0.2312 ，b= 1.8424。图2-7商品三交货时间的概率密度图图2-8 验证商品三交货时间是否正态分布的QQ图2求解方法和结果：该问题是一个比较复杂的单变量求极值的问题，因为含有复杂而繁琐的积分项，常用的通过求导找解析解的方法就难以适用。对于问题2中给出的具体问题，由于数据取值均不大，可以采用遍历搜索方法。从=0开始直至=，依次得出不同值下的总损失费用，比较选出其中最小值即为最低总损失费用，此时对应的即为最低订货点。另外，从优化角度看，本题属于非线性优化问题，也可从利用matlab中的相关函数，将上面列出的目标函数及其约束条件代入求解。

21、表2-4中列出了三种商品在采用遍历搜索和优化方法求解的结果。可以看出两者相差不大，但优化方法能找到更精确的解。在附表中给出了遍历搜索的过程。表2-4 题目2的求解结果遍历搜索方法优化方法订货点周期日平均损失费订货点周期日平均损失费康师傅444.30563.389743.99664.30563.3897心相印443.60164.346944.47473.56994.3459中汇香米352.200811.358235.15352.193111.35803 多品种存贮问题 3.1符号约定：商品种类；：第种商品的销售速率（），；：第种商品单个商品的体积（）；：体积的销售速率（），;：每次进货的定货费（

22、元）；：使用自己仓库存贮时，第种商品每天的存贮费（），；：租借仓库存贮时，第种商品每天的存贮费（），；：缺货时第种商品的损失为（），；：每次订货后交货的滞后时间为（天）；：的概率密度函数；：自己仓库最大容量（）；：自己仓库用于存贮第种商品的最大容量（），；：商品存贮量达到的固定值（）；：第种商品订货到达时存贮量达到的固定值（），；：订货点（）；：达到订货点时每种商品的库存量（），；：总日均损失费用；：每种商品的日均损失费用， ；3.2多品种存贮问题的分析和建模1 问题分析和模型建立同第一题中的分析，这里仍然认为时间和库存量都是连续的变量。总损失费用最低等价于单位时间（即天）的平均损失费用（即日

23、均损失费用）最低。当有种商品同时订货时，单独看每种商品的进货、存贮、销售、订货、缺货流程和第一题中建立的模型是相同的，故当和时，的表达式推导过程与单品种模型中完全相同，这里就不再做讨论了。在第一题中不会出现商品缺货了还没有订货的情况，同理在这个问题中也不会出现种商品都缺货了还没有订货的情况。但是也有两点不同：1）这里的存贮费用和缺货费用都是按体积衡量的，而不是按袋（盒）。2）在这里由于目标函数是多种商品的费用总和的平均，要进行总体的规划使总费用最小，当某种商品的缺货费用较低时也可能出现这种商品已经缺货，但别的商品仍在销售，订货还没有发出。即相对于第一题的模型，不仅可能出现，的情况，还可能出现的

24、情况。类似于第一题，作图进行分析: 图3-1 物品库存量及日损失费随时间的变化曲线上图3-1是一个周期内物品库存量和损失费用变化速率随时间的变化曲线。对应租借仓库的货恰销售完，而自己仓库库存为的时刻；对应物品存贮货物全部售完的时刻；为总库存降至订货点，发出订货单的时刻；为订货到达的时刻。在整个周期内，库存量以速率减少。在段内，物品保持缺货状态，但是订货单未发出。段长度为滞后期。设，即从0时刻开始至发出订货单经历了时间。这一时间可由订货点不为0的物品求出。设物品有，则。A点日损失费为，B点为，C点为0，D点。与单物品存贮模型相同方法分析，算得：则时，在一个周期内的总损失费用为： （3-1）这里没

25、有把订货费用写进，这是因为这种商品只付一次订货费用就可以了。这样对第种商品，在一个周期内的平均每天损失费用为 （3-2）其中：下标区分不同的商品种类，上标区分的不同取值范围,前两种情形下的计算式完全类似于第一题中得到的计算式，且有：。，得到这种商品平均损失费用为 （3-3）对的说明：记维向量，则的表达式与每一个分量的取值范围都有关。其中每一个的取值范围都有三种，根据排列组合理论的取值范围有种，再减掉同时为0，即同时卖完货时才订货的那种可能性，的表达式共有种可能。不同的表达式对应的取值范围不同，这种可能涵盖了所有的取值。在的每一个取值范围内即的表达式是确定的，即是一个分段函数，它的表达式具体写出

26、为： （3-4）其中，代表周期平均值，()的表达式见式（3-2）。得到了目标函数，再加上根据题意得到的各个变量约束条件，则可以得到一个多变量的最优化问题： （3-5）2 对于上述最优化问题自由变量个数的讨论 在上述最有化问题中，目标函数中含有个变量，。它们有隐含的互相决定关系。比如当确定，确定时（个变量），所有的就是唯一确定的。这是因为各种商品的最初存贮体积确定，体积销售速率已知，则每种商品的存贮量变化曲线是确定的，则总存贮量的变化曲线也是确定的（把各种商品的销售曲线按相同时间点相加即可），则唯一地对应于一个时刻，在每一种商品的销售曲线上，对应的存贮量就是。当自己仓库的存贮量（个变量）也确定时

27、，目标函数的值即日均损失费用也就确定了。即自由变量为个。再根据模型中还有两个等式约束，得到所建立的最优化模型的目标函数的自由变量为个。3 模型的求解思路解法一：用现成软件的带约束最优化函数求解上面的最优化模型，因为实际上商品的体积是离散的，不可能有半个商品出现，还需要再将得到的最优解转换为相应的离散最优解。解法二：分别写出每一个取值范围上的的表达式，和相应的约束条件。选取合适的个变量目标函数分别对它们求偏导，令结果为0，则可以解出此范围内使最小的解。将所有种情形都求解之后，进行比较，最小的值对应的就是这个问题的最优解。但是这种方法只适用于取值小而且密度函数的表示比较简单的情形。3.4多品种存贮

28、问题的求解与结果可以利用以上建立的多品种存贮模型来求解具体的问题。下面将对问题4进行求解。1 具体问题的分析及模型简化针对题目四提供的相关数据，分析得到下表 所示的三重商品各类费用比较。可以看出任何一种商品的缺货费用都远远高于其他费用。如果某种商品出现了缺货，缺货费的增加是其他存储费的减少所不能抵消的；另一方面，缺货的出现会使销售周期增大，日平均订货费用会相应减少，但是本题中一次订货、也就是一个周期内的订货费用为10元，故日平均订货费用的减少也远不足以抵消缺货费的增加。也就是说，本题中要实现总损失费用降低，在订货时各种商品都要有一定量的存储，也即订货点都大于0。表3-1三种商品单位时间内各类费

29、用的比较商品（元/日）（元/日）（元/日）康师傅面0.240.1211.4心相印手帕纸0.450.1522.5中汇香米1.61.225另外，如前面分析，3种商品有23-1=5个自由变量。这里取为、，其他的参数可以通过这5个自由变量求取：记订货时，各产品已销售的时间为，易知，。将、代入式（3-2）得到商品一的损失费用，同理可得商品二、三的损失费用、，至此得到数学模型如下：(，) （3-6）2 滞后期概率密度题中指出滞后期服从在1天到3天之间的均匀分布，我们采用连续的方法来处理问题，认为的取值在1天至3天内是连续的。在下图给出了取1天至3天的概率密度，01之间表示1天，12表示2天，23表示3天。

30、图3-2 滞后期的概率密度曲线写成表达式为：3 求解结果Matlab优化工具箱中有函数fmincon，它可以用于求解非线性带约束条件最优化问题。对于问题3中建立起来的模型，将问题4中给出的具体数据带入，调用此函数得到的结果如下。表3-2寻优结果（单位：）商品康师傅面1.45161.37172.0334心相印手帕纸1.42951.36712.0114中汇香米4.01533.26145.9552总计6.8967610从表中看出，各个量的取值均精确到了小数点后四位，而、，、应为的整数倍。所以要把表格中的各值调整到相应的整数倍。于是得到表3-3。表3-3调整结果（单位：）商品康师傅面1.451.402

31、.00心相印手帕纸1.441.402.00中汇香米4.03.26.0总计6.89610表中数据虽然不是数值上的最优解，但是符合实际情况，是最贴近最优数值解的真实解。可见，最优订货点为，中汇香米的存贮体积要比其他商品大，这是因为它的单位商品体积最大，而销售速率最高。另外算得相应的最低日平均损失为4.18元/日。3.5多品种存贮问题的最优性原理的提出和应用1 原理的提出和证明借鉴控制论中的最优性原理的想法，我们提出如下最优性原理：若模型中有一个局部的小模型，它含有个变量，是原模型所含个变量的一部分。且局部模型的变量取值对模型的其他部分无影响，则原模型的最优解必定是局部模型的最优解。证明：反证法记为

32、原模型的目标函数，最优解使取最小值。为局部模型的目标函数，M为模型其他部分，它的值由最优解中另外个变量决定，对它没有影响。若原模型的最优解的一部分，不是小模型的最优解，则求解小模型得到一个最优解满足：则这个结果与是最优解矛盾。故定理成立。根据题意，我们想找到与其他部分的费用无关的局部模型，那么得到的最优解必定也是使这个局部模型费用最低的最优解。考虑当模型符合在货物送达时租借仓库的产品已经全部卖完的特定情形下，这时租借仓库中各种货物所占的体积除了对租借仓库比用自己的仓库多花的那部分贮存费有影响之外，对其他的费用都没有影响。从现实意义上理解这主要是因为多花的那部分费用相当于用自己的仓库免费、租仓库

33、的费率为时，存贮需要花的钱。在我们的前提假设下当商品总量一定时贮存在那个仓库对于销售，缺货等流程完全没有影响。即这部分的费用与其他部分的费用是独立的。应用上面提出的最优性原理，得到（3-2）的最优解应使 ， （3-7）也取最小值。将约束条件代入（3-2）式，则（3-2）式中含有个变量，分别对它们求偏导数可以得到个线性方程，解之就可以求得。2 最优性原理的应用1） 在满足原理前提条件时，和可以相互决定，自由变量减少个。2） 对结果的验证：在第四题中用计算机程序求解得到的值为6.4254，即订货在租借仓库内的商品将要卖完时发出，而到货时间服从1天到3天的均匀分布，我们计算时将这段时间连续化为天。所

34、以基本符合租借仓库的商品卖完之后货才送到的前提。将第四题中数据带入(3-7)式中求导之后解线性方程可得最优解应当基本满足如下一组关系式：，列表比较如下：表3-4：在最优方案中租借仓库存贮量计算值和验证值的比较商品种类康师傅心相印中汇香米计算值0.71090.69892.5901验证值0.77920.62342.5974计算值与验证值之间稍有差别，这是因为这是一个随机模型，计算滞后时间时是从0开始取的，模型有大约3%的概率不符合原理的前提。两组数据基本相符，这从另外一方面说明了本模型的计算结果是比较可信的。4进一步讨论对于问题5将进行讨论。题中对于销售速率、订货费用、存贮费用、缺货费用等假定是固

35、定的，这些都是做了一定的简化。而实际的情况要复杂的多。以销售速率为例，题中假设销售速率固定，即单位时间内销售量不改变。而显见的，在实际当中，销售情况不可能时时一样，而是要变化三的。只能说它是以某个规律、某个分布函数随机变化的。这一分布函数可以由商家在此前统计的资料进行分析确定。假定销售速率取值的概率密度函数为，的出现增加了问题的复杂度。在单品种存贮模型中，建立总损失费用的模型为：取的是对的概率平均值，在考虑销售速率的随机性后，还应该对取上的概率平均值，可得到：（3-8）至此得到考虑了销售速率随机性的总损失费用模型：（3-9）同样，多品种存贮模型里也可如上考虑。要取求得的总损失费用在上的概率平均

36、值，这才是最终的目标函数。求出的总损失费用是单位时间内的平均值，它是长时间内该值的一个趋势。而实际的，随着时间的变化，总损失费用是要发生变化的。商品的滞后期和销售速率都是随机变量，因此有可能随着时间增加，模型变得偏离实际情况。此时就要对模型进行修正。于是，在经过某段时间后，商家应根据此前搜集的各个变量的数据对模型进行调整。按照问题1、3的方法，重新确定例如变量、的概率密度函数。再将其带入目标函数中对模型求解，得出新的订货点。销售速率改变时，商品的最大存贮量也应随之改变。销售速率大时，应多存一些以备需求；销售速率小时，应少存一些以避免不必要的存贮费。即存贮量要根据需求来确定。由商家的资料，可以得

37、到销售速率关于时间的函数，利用可以对接下来一个周期内商品最大存贮量进行预测。考虑单品种模型。设、为接下来一个周期的起点和终点时刻，在上对进行积分，该积分值就是预测的下一周期的最大存贮量，即。需要说明的是，此处取周期是一个概率平均值，随着、的变化，也要不时地进行调整修正。5模型评价文中建立了仓库容量有限条件下单品种随机存贮模型和多品种随机存贮模型。时间变量的连续化更合理地反映了实际情况。单品种问题建模时对总损失费用取滞后时间的概率平均值，比较合理地描述了随机存贮模型。多品种问题建模时采用分情况讨论的方法把单品种模型作了推广，并讨论了目标函数中自由变量的个数，使复杂模型清晰化，建模过程简洁，求解方

38、便。模型的稳定性好，实用性强。本模型的缺点是由于样本点的取值比较少，概率分布密度函数估计的不够精确，这可能会对求解结果造成比较大的影响。参考文献1黄洁纲，存贮论原理及其应用，上海科学技术文献出版社，19842胡运权，运筹学基础及应用，哈尔滨工业大学出版社，20023韩宇鑫 孙晖 杨青双，瞬时进货、允许缺货的租赁库存系统存贮模型，辽宁工学院学报，23:2，20034邢继祥 张春蕊 徐洪泽，最优控制应用基础，科学出版社，2003附表1 对康师傅面进行L的遍历搜索结果12345677.88897.80567.72227.63897.55567.47227.38899.14378.89318.6470

39、8.40548.16867.93667.70958910111213147.30567.22227.13897.05566.97226.88896.80567.48757.27057.05886.85236.65136.45586.2658151617181920216.72226.63896.55566.47226.38896.30566.22226.08165.90305.73035.56355.40265.24775.0989222324252627286.13896.05565.97225.88895.80565.72225.63894.95624.81964.68924.56504.

40、44704.33534.2297293031323334355.55565.47225.38895.30565.22225.13895.05564.13044.03743.95053.86973.79523.72673.6643363738394041424.97224.88894.80564.72224.63894.55564.47223.60793.55743.51293.47423.44143.41093.3914434445464748494.38894.30564.22224.13894.05563.97223.88893.39043.38973.39213.39393.41783.45253.4990505152535455563.80563.72223.63893.55563.47223.38893.30563.55863.63253.72223.82953.95624.10444.2765575859603.22223.13893.05562.97224.47524.70344.96435.2615附表2 对心相印面进行L的遍历搜索结果12345676.46826.40166.33496.26826.20166.13496.068214.053713.630013.213312.803612.401212.00