南京邮电大学 《数学建模》课程的一次书面大作业:
题目
试建立不允许缺货的生产销售存贮模型。设生产速率为常数 k,销售速率为常数 r,k>r,在每个生产周期T内,开始的一段时间( 0 < t < T 0 0<t<T_0 0<t<T0)一边生产一边销售后来的一段时间( T 0 < t < T T_0<t<T T0<t<T)只销售不生产,画出存贮量 q(t) 的图形。设每次生产准备费为 c 1 c_1 c1,单位时间每件产品贮存费为 c 2 c_2 c2,以平均每天总费用最小为目标确定最优生产周期。讨论 k ≫ r k\gg r k≫r 和 k ≈ r k\approx r k≈r 的情况。
解答
为了建模方便,假设生产量和销售量的变化都是连续的,生产周期T可以不是整数。
-
题目内有以下常数:
含义 符号 生产速率 k k k 销售速率 r r r 生产准备费 c 1 c_1 c1 单位时间每件产品贮存费 c 2 c_2 c2 -
T 0 T_0 T0前,q(t)以 k − r k-r k−r的速率增长,达到最大值。
-
T 0 T_0 T0停止生产后,库存数逐渐减少,直至库存数变为0,此时一个周期结束。
![image-20240320223801671](https://img-blog.csdnimg.cn/img_convert/f2b4b90e927cb1f7ce06d5f970f5f4e8.png)
所以,这段周期内的总费用为:总存储费用+生产准备费。
有:
c
=
c
2
∫
0
T
0
q
(
t
)
d
t
+
c
2
∫
T
0
T
q
(
t
)
d
t
+
c
1
c=c_{2}\int_{0}^{T_0}q(t)dt+c_{2}\int_{T_{0}}^{T}q(t)dt+c_{1}
c=c2∫0T0q(t)dt+c2∫T0Tq(t)dt+c1
化简后有:
c
=
c
2
(
k
−
r
)
T
0
2
2
+
c
2
r
(
T
−
T
0
)
2
2
+
c
1
\begin{equation} c=c_{2}\frac{(k-r)T_{0}^{2}}{2}+c_{2}\frac{r\left(T-T_{0}\right)^{2}}{2}+c_{1} \tag{1} \end{equation}
c=c22(k−r)T02+c22r(T−T0)2+c1(1)
根据库存的数量关系,有以下等式:
(
k
−
r
)
T
0
=
(
T
−
T
0
)
r
k
T
0
−
r
T
0
=
r
T
−
r
T
0
(2)
\tag{2} (k-r)T_0=(T-T_0)r\\ kT_0-rT_0=rT-rT_0
(k−r)T0=(T−T0)rkT0−rT0=rT−rT0(2)接着化简
k
T
0
=
r
T
即
T
0
=
r
T
k
(3)
kT_0=rT\\ \tag{3} 即T_0=\frac{rT}{k}
kT0=rT即T0=krT(3)
将(2)和(3)代入到关于式子(1)中有
c
=
c
2
(
k
−
r
)
T
0
2
2
+
c
2
r
(
T
−
T
0
)
2
2
+
c
1
=
c
2
2
[
r
T
0
(
T
−
T
0
)
+
(
k
−
r
)
T
0
(
T
−
T
0
)
]
+
c
1
=
c
2
2
[
k
T
0
(
T
−
T
0
)
]
+
c
1
=
c
2
2
×
r
T
×
[
T
−
r
T
k
]
+
c
1
=
c
2
2
×
[
k
−
r
k
]
×
r
T
2
+
c
1
\begin{align} \tag{} c&=c_{2}\frac{(k-r)T_{0}^{2}}{2}+c_{2}\frac{r\left(T-T_{0}\right)^{2}}{2}+c_{1}\\\tag{} &=\frac{c_2}{2} \left [ rT_0(T-T_0)+(k-r)T_0(T-T_0) \right ] +c_1\\ \tag{} &=\frac{c_2}{2} \left [ kT_0(T-T_0) \right ] +c_1\\ \tag{} &=\frac{c_2}{2} \times rT \times\left [ T-\frac{rT}{k} \right ] +c_1\\ \tag{4} &=\frac{c_2}{2} \times\left [ \frac{k-r}{k} \right ] \times rT^2 +c_1\\ \end{align}
c=c22(k−r)T02+c22r(T−T0)2+c1=2c2[rT0(T−T0)+(k−r)T0(T−T0)]+c1=2c2[kT0(T−T0)]+c1=2c2×rT×[T−krT]+c1=2c2×[kk−r]×rT2+c1()()()()(4)
对(4)式除以周期T,求平均每天的费用
c
ˉ
\bar{c}
cˉ :
c
ˉ
=
c
T
=
c
2
2
×
[
k
−
r
k
]
×
r
T
+
c
1
T
=
c
2
r
(
k
−
r
)
T
2
k
+
c
1
T
\begin{align}\tag{} \bar{c} &= \frac{c}{T} \\ \tag{} &=\frac{c_2}{2} \times\left [ \frac{k-r}{k} \right ] \times rT +\frac{c_1}{T} \\ \tag{5} &=\frac{c_2r(k-r)T}{2k} +\frac{c_1}{T} \\ \end{align}
cˉ=Tc=2c2×[kk−r]×rT+Tc1=2kc2r(k−r)T+Tc1()()(5)
式(5)即为目标函数,为使
c
ˉ
\bar{c}
cˉ 最小,有:
d
c
ˉ
d
T
=
c
2
r
(
k
−
r
)
2
k
−
c
1
T
2
=
0
\frac{\mathrm{d} \bar{c} }{\mathrm{d} T} = \frac{c_2r(k-r)}{2k} -\frac{c_1}{T^2} = 0
dTdcˉ=2kc2r(k−r)−T2c1=0
化简可得:
T
2
=
2
k
c
1
c
2
r
(
k
−
r
)
T
=
2
k
c
1
c
2
r
(
k
−
r
)
\begin{align}\tag{} T^2 = \frac{2kc_1}{c_2r(k-r)} \\ \tag{6} T = \sqrt{\frac{2kc_1}{c_2r(k-r)}} \\ \end{align}
T2=c2r(k−r)2kc1T=c2r(k−r)2kc1()(6)
式子(6)即为我们要求的最佳生产周期。
-
当 k ≫ r k\gg r k≫r 时,
T = 2 c 1 c 2 r T = \sqrt{\frac{2c_1}{c_2r}} T=c2r2c1
生产周期与k无关。 -
当 k ≈ r k\approx r k≈r 时, T → + ∞ T\to + \infty T→+∞ ,生产周期无限长,几乎无法产生贮存量。