南京邮电大学 《数学建模》课程的一次书面大作业:

题目

试建立不允许缺货的生产销售存贮模型。设生产速率为常数 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 rkrk≈rk\approx rkr 的情况。

解答

为了建模方便,假设生产量和销售量的变化都是连续的,生产周期T可以不是整数。

  1. 题目内有以下常数:

    含义 符号
    生产速率 kkk
    销售速率 rrr
    生产准备费 c1c_1c1
    单位时间每件产品贮存费 c2c_2c2
  2. T0T_0T0前,q(t)以k−rk-rkr的速率增长,达到最大值。

  3. T0T_0T0停止生产后,库存数逐渐减少,直至库存数变为0,此时一个周期结束。

image-20240320223801671

所以,这段周期内的总费用为:总存储费用+生产准备费。

有:
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=c20T0q(t)dt+c2T0Tq(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(kr)T02+c22r(TT0)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 (kr)T0=(TT0)rkT0rT0=rTrT0(2)接着化简
kT0=rT即T0=rTk(3) kT_0=rT\\ \tag{3} 即T_0=\frac{rT}{k} kT0=rTT0=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(kr)T02+c22r(TT0)2+c1=2c2[rT0(TT0)+(kr)T0(TT0)]+c1=2c2[kT0(TT0)]+c1=2c2×rT×[TkrT]+c1=2c2×[kkr]×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×[kkr]×rT+Tc1=2kc2r(kr)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(kr)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(kr)2kc1T=c2r(kr)2kc1 ()(6)
式子(6)即为我们要求的最佳生产周期。

  1. k≫rk\gg rkr 时,
    T=2c1c2r T = \sqrt{\frac{2c_1}{c_2r}} T=c2r2c1
    生产周期与k无关。

  2. k≈rk\approx rkr 时,T→+∞T\to + \inftyT+ ,生产周期无限长,几乎无法产生贮存量。

Logo

GitCode 天启AI是一款由 GitCode 团队打造的智能助手,基于先进的LLM(大语言模型)与多智能体 Agent 技术构建,致力于为用户提供高效、智能、多模态的创作与开发支持。它不仅支持自然语言对话,还具备处理文件、生成 PPT、撰写分析报告、开发 Web 应用等多项能力,真正做到“一句话,让 Al帮你完成复杂任务”。

更多推荐