数学建模作业(1)

习题一在3.1节存储模型中的总费用中增加购买货物本身的费用，重新确定最优订货周期和订货批量。

证明在不允许缺货模型和允许缺货模型中结果都与原来一样。

一、不允许缺货的存储模型
问题分析若生产周期短、产量少，会使存储费用小，准备费用大，货物价格不变；
而周期长、产量多，会使存储费大，准备费小，货物价格不变。

所以必然存在一个最佳周期，使总费用最小。

显然，应建立一个优化模型。

模型假设为了处理的方便，考虑连续模型，即设生产周期T和产量Q为连续量。

根据问题性质作如下假设：
（1）产品每天的需求量为常数r。

（2）每次生产费用为c1，每天每件产品存储费为c2，购买每件货物所需费用为c3.（3）生产能力为无限大（相对于需求量），当存储量降为零时，Q件产品立即生产出来供给需求，即不允许缺货。

模型建立将存储量表示为时间t的函数q（t），t=0生产Q件，存储量q(0)=Q，q(t)以需求速率r递减，直到q(T)=0，如图，显然有：Q=rT
图（1）不允许缺货模型的存储量q（t）
一个周期内的存储费是c2∫q(t)dt,其中积分恰好等于图中三角形面积QT/2，因为一个周期的准备费是c1，购买每件货物的费用为c3，得到一个周期的总费用为：
C=c1+c2QT/2+r Tc3=c1+c2 r T2/2+ r T c3
则每天的平均费用是
C(T)=c1/T+r c3+c2 r T/2
上式为这个优化模型的目标函数。

模型求解求T使上式的C最小。

容易得到
T=√2c1/（c2r）则Q=√2c1r/c2
二、允许缺货的存储模型
(1) 模型假设产品每天的需求量为常数r。

(2) 每次生产费用为c1，每天每件产品存储费为c2，购买每件货物所需费用为c3.
(3) 生产能力为无限大（相对于需求量），允许缺货，每天每件损失费为c4,但缺货数量需在下次生产（或订货）时补足。

，
模型建立因存储量不足造成缺货时，可以认为存储量函数q(t)为负值，如图所示，周期仍记为T,Q是每周期初的存储量，当t=T1时q(t)=0,于是有Q=r T1
t
图（2）允许缺货模型的存储量q(t)
在T1到T这段时间内需求率r不变，q(t)按原斜率继续下降。

由于规定缺货量需补足，所以在t=T时数量为R的产品立即到达，使下周期初的存储量恢复为Q.
所以C=c1+c2QT1/2+ r Tc3+c4r(T-T1)2/2
将模型的目标函数------每天的平均费用------记作T和Q的二元函数
C(T,Q)=c1/T+c2Q2/(2rT)+ +r c3+c3(Rt-Q)2/(2Tr)
模型求解利用微分法求T和Q使C(T,Q)最小，令dC/dT=0 ,dC /dQ =0,可得
T’=√2c1(c2+c4)/(rc2c4) , Q’=√2c1rc4/(c2(c2+c3))
由以上两个模型可以看出在不允许缺货模型和缺货模型中结果都与原来一样
存储模型
问题：建立不允许缺货存储模型。

设生产速率为常数k，销售速率为常数r,k＞r,在每个生产周期T内,开始的一段时间（0＜t＜T0﹚一边生产一边销售,后来的一段时间﹙T0＜t＜T﹚只销售不生产.画出储存量q(t)的图形,设每次生产准备费为c1,单位时间每件产品储存费为c2,以总费用最小为目标确定最优生产周期。

讨论K＞＞r和K≈r的情况。

问题分析：在t< T0时间内k＞r有储存量以k－r速率增加，在T0＜t＜T时间内，储存量以r 速率递减，而一个周期的总费用为C（t）与生产周期，产量与需求量，生产准备费，储存费之间的关系，从而建立数学模型。

可根据数学最值定理求出最优周期。

模型假设：为处理方便考虑连续模型，即生产周期T和生产量Q都是连续量，故作出如下假设
1.生产速率为常数k，销售速率为常数r,k＞r。

2.每次生产准备费为c1,单位时间每件产品储存费为c2。

3.开始的一段时间（0＜t＜T0﹚一边生产一边销售,后来的一段时间﹙T0＜t＜T﹚
只销售不生产.
4.不允许缺货。

模型建立：将存储量表示为时间T的函数q(t)，在0＜t＜T0，q(t)以的速率增加；在T0＜t＜T 时间内，q(t)以r的速率递减，直到q(t) ＝0 ，如图所示，显然有﹙k－r﹚t 0＜t＜T0
q(t)=
﹙k－r﹚T0－r﹙t－T0﹚T0＜t＜T
﹙k－r﹚T0
0 T0T t
则一个周期内的储存费为c2∫0t q(t)dt
一个周期的总费用为C（t）= c1+0.5 c2﹙k－r﹚T0 T
则平均费用为g(t)=C(t)/T= c1 /T+0.5 c2﹙k－r﹚T02/T +0.5 c2﹙k－r﹚T0（T－T0）/T
由于﹙k－r﹚T0=r（T－T0）
故g(t)= c1 /T+c2﹙k－r﹚rT/2k
模型求解：
使平均费用最小的最优周期为T*=√2c1k/c2r（k-r）
模型分析：当k＞＞r时，T*=√2c1/c2r 相当于不考虑生产的情况。

当K≈r时，T* ∞，因为产量被销售量抵消，无法形成储存量。

model:
min=100*x1+40*x5+40*x6;
x1+x6>4;
x2+x4>6;
x3+x5>5;
x1+x5>8;
x1+x5-x4>8;
x2+x3=x1;
end
Global optimal solution found.
Objective value: 640.0000
Infeasibilities: 0.000000
Total solver iterations: 7
Variable Value Reduced Cost X1 4.000000 0.000000 X5 6.000000 0.000000 X6 0.000000 20.00000 X2 4.000000 0.000000 X4 2.000000 0.000000 X3 0.000000 40.00000
Row Slack or Surplus Dual Price
1 640.0000 -1.000000
2 0.000000 -20.00000
3 0.000000 -40.00000
4 1.000000 0.000000
5 2.000000 0.000000
6 0.000000 -40.00000
7 0.000000 40.00000。