数学建模 书面大作业 存贮模型
设生产速率为常数 k,销售速率为常数 r,k>r,在每个生产周期T内,开始的一段时间(为了建模方便,假设生产量和销售量的变化都是连续的,生产周期T可以不是整数。)只销售不生产,画出存贮量 q(t) 的图形。停止生产后,库存数逐渐减少,直至库存数变为0,此时一个周期结束。,以平均每天总费用最小为目标确定最优生产周期。所以,这段周期内的总费用为:总存储费用+生产准备费。对(4)式除以周期T,求平均每天
南京邮电大学 《数学建模》课程的一次书面大作业:
题目
试建立不允许缺货的生产销售存贮模型。设生产速率为常数 k,销售速率为常数 r,k>r,在每个生产周期T内,开始的一段时间(0<t<T00<t<T_00<t<T0)一边生产一边销售后来的一段时间(T0<t<TT_0<t<TT0<t<T)只销售不生产,画出存贮量 q(t) 的图形。设每次生产准备费为 c1c_1c1,单位时间每件产品贮存费为 c2c_2c2,以平均每天总费用最小为目标确定最优生产周期。讨论 k≫rk\gg rk≫r 和 k≈rk\approx rk≈r 的情况。
解答
为了建模方便,假设生产量和销售量的变化都是连续的,生产周期T可以不是整数。
-
题目内有以下常数:
含义 符号 生产速率 kkk 销售速率 rrr 生产准备费 c1c_1c1 单位时间每件产品贮存费 c2c_2c2 -
T0T_0T0前,q(t)以k−rk-rk−r的速率增长,达到最大值。
-
T0T_0T0停止生产后,库存数逐渐减少,直至库存数变为0,此时一个周期结束。

所以,这段周期内的总费用为:总存储费用+生产准备费。
有:
c=c2∫0T0q(t)dt+c2∫T0Tq(t)dt+c1 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=c2(k−r)T022+c2r(T−T0)22+c1 \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)T0=(T−T0)rkT0−rT0=rT−rT0(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)接着化简
kT0=rT即T0=rTk(3) kT_0=rT\\ \tag{3} 即T_0=\frac{rT}{k} kT0=rT即T0=krT(3)
将(2)和(3)代入到关于式子(1)中有
c=c2(k−r)T022+c2r(T−T0)22+c1=c22[rT0(T−T0)+(k−r)T0(T−T0)]+c1=c22[kT0(T−T0)]+c1=c22×rT×[T−rTk]+c1=c22×[k−rk]×rT2+c1 \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ˉ=cT=c22×[k−rk]×rT+c1T=c2r(k−r)T2k+c1T \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ˉ 最小,有:
dcˉdT=c2r(k−r)2k−c1T2=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
化简可得:
T2=2kc1c2r(k−r)T=2kc1c2r(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≫rk\gg rk≫r 时,
T=2c1c2r T = \sqrt{\frac{2c_1}{c_2r}} T=c2r2c1
生产周期与k无关。 -
当 k≈rk\approx rk≈r 时,T→+∞T\to + \inftyT→+∞ ,生产周期无限长,几乎无法产生贮存量。
更多推荐
所有评论(0)