存储问题的数学模型

存储问题的数学模型第一页，共五十七页，2022年，8月28日一、存贮问题1、工厂订购原料，存入仓库供生产所用；2、车间一次加工生产一批零件，供装配线每天生产之需；3、放在货柜里以备零售；4、水库在雨季蓄水，用于旱季灌溉和发电；5、自动取款机每天早上一次性放入若干货币，以供消费者自由取款；6、计算机的硬盘里的内存的预留问题；等等。

这些问题，都存在一个存储多少才合适的问题，存储过多，存储费用、积压货物本身的资金的流动、资金本身的利息等太多；存储太少，会导致订购货物次数太多而增加采购费或不能及时满足消费需求而增加其它费用。[姜起源.谢金星.数学模型.高等教育出版社.59页]第二页，共五十七页，2022年，8月28日二、与存贮问题有关的几个术语1、需求

所谓需求，就是将存贮在“仓库”的“货物”提供给需要的一方。存贮的目的就是为了满足需求，随着需求的发生，存贮将减少。需求分类：（1）连续性需求和间断性需求

若需求是连续的，则存贮就是连续减少；若需求是间断的，则存贮是跳跃式减少（需求是瞬间完成的）。（2）确定性需求和随机性需求

需求发生的时间和数量是确定的，就是确定性需求；否则就是随机性需求。第三页，共五十七页，2022年，8月28日2、补充

所谓补充就是向“仓库”里添加存贮物，使得存储量得以增加或恢复。补充不及时或者补充不足，就无法满足新的需求。补充的方式：（1）提前定货，货到就即刻补充，并使得存贮量得到规定的量；且一般认为这种补充是瞬间完成的。（2）自己安排生产，一般生产速度（单位时间内生产的货物数量）大于需求速度（单位时间内对货物的需求数量）。3、存贮状态图

在一个存贮系统中，存贮量因需求而减少，因补充而增加。在直角坐标系下，以时间t为横轴，以实际存贮量S为纵轴，将每个时刻t对应的存贮货物量S反应在坐标平面，就得到一条反应实际存贮量S的动态变化规律的曲线，称为存贮状态图。第四页，共五十七页，2022年，8月28日例如，开始有Q这么多的货物量，每单时间的货物需求数量（连续需求速度）为r，若在[0，T]这段时间内货物全部供应给需求者，则仓库的货物数量的状态图为tS0TQS(t)=Q-rt连续需求的存贮状态曲线从图中可以看出，存贮状态随时间变化而变化，且S(0)=Q，S(T)=0。当S=0时，瞬间补充货物量，使得存贮量又恢复到Q（补充量为Q）补充第五页，共五十七页，2022年，8月28日又如，刚开始有500单位货物，需求是离散需求，每10天需求100单位货物，则存贮状态图为St0T补充500离散需求的存贮状态图第六页，共五十七页，2022年，8月28日3、费用

在研究存贮问题中，总是用费用来评价存贮策略的优劣性。为了正确评价一个系统的不同策略的优劣性或评价不同存贮系统的优劣性，必须满足可比性要求：时间一致的可比性和计算口径（项目一致）的可比性。常用一个存贮周期内的总费用或单位时间平均总费用来衡量一个存贮系统的优劣性。在研究存贮问题时，经常要考虑的是如下一些项目的费用：（1）存贮费：存贮货物的资金利息、保险、仓库租借费用、保管费、维修费、物资自然耗损费等，一般与存贮时间和存贮数量成比例；（2）缺货费：存贮不能满足需求而造成的损失，如失去销售机会、停工待料的损失、延期交货的额外支付、对需求方的赔偿。缺货费与缺货数量和时间成比例，当不允许缺货时，缺货费为无穷大。第七页，共五十七页，2022年，8月28日（3）订货费或生产费（3.1）当向外采购货物补充存贮量时，就发生订货费，这包括两部分：一部分是订购费，例如手续费、差旅费等，它与订货次数有关，而与订货数量无关；另一部分就是物资进货成本，如货款、运费等，只与订货数量有关（货款与运费可能会有所波动）。（3.2）如果自行生产货物而补充物资存贮量，则需要生产费，它也包括两部分：一部分是固定费用（装配生产线、生产准备费用、或组织调整生产线的费用等），它与生产次数有关，与生产数量无关；另一部分是与生产量有关的费用，如原材料和零部件的成本、直接加工费用等。第八页，共五十七页，2022年，8月28日4、存贮策略

什么叫存贮策略呢？就是决定以什么方式对存贮进行补充，什么时间补充，补充多少。常见的存贮策略有以下几种：（1）t-循环策略：每隔固定的时间t，补充一个固定的量Q；（2）(t,S)策略：每隔固定时间t补充一次，补充前对仓库盘点，如果物资量为I，仓库的最大库存容量为S，则补充量Q=S-I。（3）（s,S)策略：设仓库的最大容量为S，时刻监控物质的余量I，如果I大于等于s（称为警戒点、安全存储量、保险存储量），则不补充；如果I小于s，则补充货物，补充量为Q=S-I。（4）(t,s,S)策略：（自行补充）。第九页，共五十七页，2022年，8月28日三、存贮模型1不允许缺货模型三、存储模型1

配件厂为装配线生产若干种部件，轮换生产不同的部件时，因更换设备需要付出生产准备费用(与产量无关)，同一部件的产量大于需求时，因资金积压、占用仓库需要存储费用。今已知某一部件日需求量100件，生产准备费用为5000元，存储费用每天每件1元。如果生产能力远远大于需求，并且不允许出现缺货，试安排该产品的生产计划，即多少天生产一次(称为生产周期)，每次产量多少，可使得总费用或平均费用(注意费用口径一致：时间一致，项目一致)最少。模型分析

如果每天安排生产一次，每次生产100件，则一天费用为5000元；但是如果这样连续生产x天，则x天总费用为5000x元。如果每隔x天生产一次，一次生产100x件，则生产准备费用为5000元，存储费用为[姜起源.谢金星.数学模型.高等教育出版社.59页]第十页，共五十七页，2022年，8月28日则每隔x天生产一次，这x天中总费用为元元平均每天的费用为其函数曲线如下(见附件1)第十一页，共五十七页，2022年，8月28日每x天生产一次的平均费用,x*为最优点x*第十二页，共五十七页，2022年，8月28日

由上面的图可以看出，存在x*，每隔x*天安排一次生产，这x*天内每天的费用都低于其它安排。称x*为生产周期。如何求这个x*，其它类似问题是否可以一样求解呢？模型假设1、为了方便，考虑连续情况，即需求是连续，生产周期和产量也是连续；2、产品的需求速度是一个确定常数；3、生产能力无穷大(相对需求)，即可以很短时间内实现任何产量，即永远不会缺货；符号设置R需求速度；Q周期内的产量；T生产周期；C1

生产准备费用；C2单位产品每天的存储费用；S货物存储量；第十三页，共五十七页，2022年，8月28日模型建立t=0时刻开始生产的货物量为Q。存储量S是时间t的函数，记为S(t)，且开始时S(0)=Q（产量），随着需求的进行，S(t)是减函数，且S(T)=0，又开始生产补充。这一过程可以用如下动态过程表示。RT=QT斜率为R0tS(t)存储状态变化图一个周期内的生产准备费用为C1，一个周期内的存储费用为第十四页，共五十七页，2022年，8月28日QTR0tS(t)存储状态变化图存储费，所以周期T内总费用周期T内每天平均费用为《1》第十五页，共五十七页，2022年，8月28日对《1》函数的驻点，得《2》将其带入Q=RT，得到《3》《2》和《3》是经济学上著名的经济订货批量公式(EOQ:)将《2》带入《1》，得到平均最小费用《4》模型求解第十六页，共五十七页，2022年，8月28日结果解释由《2》可以看出1、生产准备费用C1越大，生产周期越长；单位存储费用C2和需求速度R越大，生产周期越短。

由《3》可以看出2、生产准备费用C1越大，需求速度R越大，产量Q越大；单位存储费用C2越小，产量Q越大。

这两个结果符合常识，但是，它们之间关系呈平方根关系，只有靠数学模型得到（用途？）。

将问题已知条件C1=5000，R=100，C2=1代入《2》、《3》和《4》得到第十七页，共五十七页，2022年，8月28日敏感性分析由公式对C1,C2,R,分别求偏导数，得注：函数改变的百分比/自变量改变的百分比=函数的弹性，即函数y=f(x)在x0附近有定义，如果下式成立，《5》《2》第十八页，共五十七页，2022年，8月28日称《5》为函数y=f(x)在点x0处的弹性，反映函数y随自变量x变化的剧烈程度。即自变量变化1%时，函数y变化的百分比，若E(y)为正，表示增加的百分比；如果E(y)为负，表示下降的百分比，总是是反映函数对自变量的灵敏程度。

下面就用弹性这个概念来计算平均费用对各个参数的灵敏程度。这里C1=5000，C2=1，R=100。1、周期T对生产准备费用C1的灵敏度分析结果说明，当生产准备费用变动1%时，周期t用只是变动0.5%，说明周期对准备费用反应不灵敏。第十九页，共五十七页，2022年，8月28日2、周期T对存储费用的灵敏度分析计算结果表明，存储费用每增加1%，周期T则下降0.5%，反应不太灵敏。3、周期T对需求速度R的灵敏度分析计算结果表明，存储费用每增加1%，周期T则下降0.5%，反应不太灵敏。第二十页，共五十七页，2022年，8月28日模型改进1、在假设3中，假设生产能力无限大，即缺货瞬间补充。如果生产能力有限，且生产速度P(大于需求速度)也是一个常数，怎么建立公式《1》~《4》?0tTT1Q斜率R斜率P-R2、在此生产存储模型中，没有提到生产产品的费用？为什么？在什么情况下不考虑生产费用也可以求生产周期？（在生产过程中，生产要素所涉及的价格不改变，且需求速度不改变，即可不考虑）第二十一页，共五十七页，2022年，8月28日四、存贮模型2四、存贮模型2允许缺货的存贮模型

在某些情况下用户允许短时间缺货，虽然会造成一定的损失，但是只要损失费用不超过不允许缺货的准备费用和存储费用的话，允许缺货应该是可以采取的策略。下面讨论允许缺货的模型模型假设1、为了方便，考虑连续情况，即需求是连续，生产周期和产量也是连续；2、产品的需求速度是一个确定常数；3、生产能力无限大（相对于需求），允许缺货，每天每件产品缺货损失费用为C3，但缺货的数量在下次生产(或订货)时补足；[姜起源.谢金星.数学模型.高等教育出版社.61页]第二十二页，共五十七页，2022年，8月28日符号设置R需求速度；A最大存储量；T生产周期；C1

生产准备费用；C2单位产品每天的存储费用；S货物存储量；C3

单位缺货费用；T1

开始缺货点；B最大缺货量；Q货物产量；模型建立St0AT1T-B斜率R第二十三页，共五十七页，2022年，8月28日St0AT1T-B斜率RS(t)=A-Rt

根据问题的描述，[0,T1]上货物的存储量S(t)=A-Rt，S(0)=A，S(T1)=0,S(T)=-B。[0,T1]有货物存储，[T1,T]货物短缺，分别产生存储费和缺货费。且《6》第二十四页，共五十七页，2022年，8月28日生产准备费用C1存储费用下面求[0,T]时间短所发生的费用St0AT1T-B斜率RS(t)=A-Rt缺货费用总费用为《7》第二十五页，共五十七页，2022年，8月28日[0,T]上每天的平均费用为《8》求《8》的驻点，由方程求解得到驻点《9》第二十六页，共五十七页，2022年，8月28日《10》将《9》、《10》带入《6》，得到周期内最大存储量周期内产量《11》《12》令，通过对比《9》~《12》和《2》、《3》，有《13》第二十七页，共五十七页，2022年，8月28日RT=QT*斜率为R0tS(t)不允许存储状态变化图St0AT1T’-B斜率R允许缺货存储状态变化图《13》第二十八页，共五十七页，2022年，8月28日模型1和模型2之间关系的解释1、当T=T1时，模型2就是模型1；2、由于这与C2>0矛盾；3、由允许缺货变成不允许缺货，

是由信誉机制或者惩罚机制自然形成，若让缺货代价趋向无穷大，系统自然朝着不会缺货方向发展。即第二十九页，共五十七页，2022年，8月28日只要缺货费C3越来越大（相对于存储费C2）,则有第三十页，共五十七页，2022年，8月28日五、用规划方法研究存贮系统的方法

配件厂为装配线生产若干种部件，轮换生产不同的部件时，因更换设备需要付出生产准备费用(与产量无关)，同一部件的产量大于需求时，因资金积压、占用仓库需要存储费用。今已知某一部件日需求量100件，生产准备费用为5000元，存储费用每天每件1元。如果生产能力远远大于需求，并且不允许出现缺货，试安排该产品的生产计划，即多少天生产一次(称为生产周期)，每次产量多少，可使得总费用或平均费用(注意费用口径一致：时间一致，项目一致)最少。[姜起源.谢金星.数学模型.高等教育出版社.59页]例题1不允许缺货，补充瞬间完成，需求连续均匀第三十一页，共五十七页，2022年，8月28日模型假设1、为了方便，考虑连续情况，即需求是连续，生产周期和产量也是连续；2、产品的需求速度是一个确定常数；3、生产能力无穷大(相对需求)，即可以很短时间内实现任何产量，即永远不会缺货；符号设置R需求速度；Q周期内的产量；T生产周期；C1

生产准备费用；C2单位产品每天的存储费用；S货物存储量；第三十二页，共五十七页，2022年，8月28日RT=QT斜率为R0tS(t)存储状态变化图存储费，所以周期T内总费用周期T内每天平均费用为《1》第三十三页，共五十七页，2022年，8月28日（1）计算最佳订购量Q*和周期T*的规划模型 第三十四页，共五十七页，2022年，8月28日（1）求解求解程序：min=c1/t+1/2*c2*r*t;q=r*t;@gin(q);c1=5000;c2=1;r=100;计算结果：Globaloptimalsolutionfoundatiteration:218Objectivevalue:1000.000VariableValueReducedCostC15000.0000.000000T10.000000.000000C21.0000000.000000R100.00000.000000Q1000.0000.2489742E-06即最佳周期为T*=10(天），最佳经济订购批量Q*=1000（件）第三十五页，共五十七页，2022年，8月28日sets:qianhou/1..2/:c1,c2,r,t,c,q;endsetsdetc=c(2)-c(1);@for(qianhou:c=c1/t+1/2*c2*r*t);@for(qianhou:q=r*t);@free(detc);data:c1=5000,5000;c2=1,1;r=100,100;t=10,9;enddata（2）灵敏度分析(2.1)其它因素不变化，只要周期缩短一天，平均费用的增量计算程序。计算结果显示VariableValueDETC5.555556根据计算结果，只是周期缩减一天，平均费用增加5.56元。第三十六页，共五十七页，2022年，8月28日sets:qianhou/1..2/:c1,c2,r,t,c,q;endsetsdetc=c(2)-c(1);@for(qianhou:c=c1/t+1/2*c2*r*t);@for(qianhou:q=r*t);@free(detc);data:c1=5000,5000;c2=1,1;r=100,100;t=10,11;enddata周期增加一天的计算结果：VariableValueDETC4.545455根据计算结果，其它因素不变化，只是周期增加一天，平均费用增加4.55元。第三十七页，共五十七页，2022年，8月28日（2.2）其它因素不变化，只是需求速度有所变化，引起的平均费用的改变的计算sets:qianhou/1..2/:c1,c2,r,t,c,q;endsetsdetc=c(2)-c(1);@for(qianhou:c=c1/t+1/2*c2*r*t);@for(qianhou:q=r*t);@free(detc);data:c1=5000,5000;c2=1,1;r=100,101;t=10,10;enddata需求速度增加一个单位引起的平均费用的改变量为：VariableValueDETC5.000000根据计算结果，平均费用增加5.0（元）。第三十八页，共五十七页，2022年，8月28日sets:qianhou/1..2/:c1,c2,r,t,c,q;endsetsdetc=c(2)-c(1);@for(qianhou:c=c1/t+1/2*c2*r*t);@for(qianhou:q=r*t);@free(detc);data:c1=5000,5000;c2=1,1;r=100,99;t=10,10;enddata其它因素不改变，只是需求速度降低一个单位，根据计算：VariableValueDETC-5.000000

此时平均费用减少5元。第三十九页，共五十七页，2022年，8月28日（2.3）其它因素不改变，只是单位存储费c2改变一个单位，引起的平均费用的改变。sets:qianhou/1..2/:c1,c2,r,t,c,q;endsetsdetc=c(2)-c(1);@for(qianhou:c=c1/t+1/2*c2*r*t);@for(qianhou:q=r*t);@free(detc);data:c1=5000,5000;c2=1,1.01;r=100,100;t=10,10;enddata单位存贮费用增加1%，即由1变成1.01，平均费用的改变计算结果如下：VariableValueDETC5.000000即平均费用增加5元。第四十页，共五十七页，2022年，8月28日sets:qianhou/1..2/:c1,c2,r,t,c,q;endsetsdetc=c(2)-c(1);@for(qianhou:c=c1/t+1/2*c2*r*t);@for(qianhou:q=r*t);@free(detc);data:c1=5000,5000;c2=1,0.99;r=100,100;t=10,10;enddata单位存贮费用减少1%，即由1变成0.99，平均费用的改变计算结果如下：VariableValueDETC-5.000000即平均费用减少5元。第四十一页，共五十七页，2022年，8月28日例题2允许缺货，补充时间较短，需求连续均匀

企业生产某种产品，正常生产条件下每天可以生产10件。根据供货合同，需要按照每天7件供货。存储费每天每件0.13元，缺货费每天每件0.5元，每次生产准备费用（装配费）为80元，求最优存储策略。模型假设：（1）需求是连续均匀的，即需求速度是常数；（2）补充需要一定时间，一旦需要，即刻可以开工生产，但需要一定的生产周期，并且设生产速度是均匀的，即生产速度是常数，且生产速度大于需求速度；符号设置r，需求速度；p，生产速度；c1，单位存储费用；c2，单位缺货费用；c3，生产准备费用；S，仓库最大存贮量；s0，最大缺货量；T，周期；C，存储周期内的总费用；Cbar，周期内平均费用；q，一次产量；不考虑货物本身的价值（市场稳定情况）。第四十二页，共五十七页，2022年，8月28日绘制存贮状态图tSS00T1T2T3T区间说明：[0,T1]不缺货，存货均匀连续从S减少到0，只产生存贮费；[T1,T2]开始缺货，需求持续，达到最大缺货量S0，只产生缺货费；斜率为-r斜率为p-r第四十三页，共五十七页，2022年，8月28日[T2,T3]边生产，边满足需求，这期间产生的费用有生产费用，缺货费，没有存储费；[T3，T]边生产，边满足需求，边补充库存量，直到库存量达到最大值S为止，产生存贮费用；[0，T]为一个存贮周期；[T1，T3]为缺货周期；[T2，T]为生产周期。第四十四页，共五十七页，2022年，8月28日先计算各个阶段产生的费用：[0,T1]存储费[T1,T2]缺货费第四十五页，共五十七页，2022年，8月28日[T2,T3]缺货费与生产准备费用[T3,T]存贮费用第四十六页，共五十七页，2022年，8月28日[0，T]周期内的总费用为（用面积法就很好理解）第四十七页，共五十七页，2022年，8月28日建立数学模型s.t.第四十八页，共五十七页，2022年，8月28日计算程序：min=cbar;cbar=c/t;c=1/2*c1*s*t1+1/2*c1*s*(t-t3)+c3+1/2*s0*c2*(t3-t1);s=r*t1;s=(p-r)*(t-t3);s0=r*(t2-t1);s0=(p-r)*(t3-t2);Q=r*t;t>t3;t3>t2;t2>t1;t1>1;c1=0.13;c2=0.5;c3=80;p=10;r=7;计算结果为VariableValueReducedCostCBAR5.8878410.000000C159.99990.000000T27.174640.000000S45.291060.000000T16.4701510.000000T312.077620.000000S011.775680.000000T28.1523910.000000Q190.22240.000000

最佳周期T*=27.17天；最佳产量Q*=190.22件；最大存储量S=45.29件；最大缺货量s0=11.77件；开始在生产周期T-T2=19.02天。第四十九页，共五十七页，2022年，8月28日Globaloptimalsolutionfoundatiteration:10854Objectivevalue:5.887841VariableValueReducedCostCBAR5.8878410.000000C159.99990.000000T27.174640.000000C10.13000000.000000S45.291060.000000T16.4701510.000000T312.077620.000000C380.000000.000000S011.775680.000000C20.50000000.000000R7.0000000.000000P10.000000.000000T28.1523910.000000Q190.22240.000000第五十页，共五十七页，2022年，8月28日第五十一页，共五十七页，2022年，8月28日例题3不允许缺货，补充时间较长

某商店经销某商品，月需求量为30件，需求速度为常数。该商品每件进价300元，月存贮费为进价的2%。向工厂订购该商品时订购费每次20元，订购后5天才开始到货，到货速度为常数，每天2件，求最优订购策略。绘