供应链|经典论文解读:(s,S) 策略在动态库存下的最优性

在这里插入图片描述
文章考虑了具有订购成本(由单位成本加上重新订购成本组成)的动态库存问题。具体而言,对于每个时期,系统在中期开始是做出一系列采购决策——这些采购有助于库存的积累,并在随后的周期被需求所消耗。每时期系统会产生各种成本,一般包括:采购或订购成本 c ( z ) c(z) c(z),其中 z z z是购买量;持有成本 h ( ⋅ ) h(\cdot) h(),为时期末供应量超过需求量的函数;以及缺货成本 p ( ⋅ ) p(\cdot) p(),为周期末需求量超过供应量的函数。每个时期结束时会计算持有成本或缺货成本,而订购成本则在进行采购时计算。库存问题的目标便是选择采购决策以最小化所有成本的折现值的期望。文章在最初假设采购仅在周期开始时进行,并且交货是即时的;后续将假设拓展考虑延迟交货的情况。

设交货后的库存水平立即为 y y y,则该周期内将收取的预期持有和缺货成本为 L ( y ) L(y) L(y):

其中 φ \varphi φ是需求分布的密度函数。若考虑 n n n个时期、初始库存为 x x x的库存问题,令 C n ( x ) C_n(x) Cn(x)表示在最优配给下 n n n个时期计划的折现成本的期望值(折现因子 α ∈ [ 0 , 1 ] \alpha\in[0,1] α[0,1])。则可得出 C n ( x ) C_n(x) Cn(x)
C n ( x ) = min ⁡ y ≥ x { c ( y − x ) + L ( y ) + α ∫ 0 ∞ C n − 1 ( y − ξ ) φ ( ξ ) d ξ } , C_n(x)=\min_{y\geq x}\{c(y-x)+L(y)+\alpha\int_0^{\infty}C_{n-1}(y-\xi)\varphi(\xi)\mathrm{d}\xi\}, Cn(x)=yxmin{c(yx)+L(y)+α0Cn1(yξ)φ(ξ)dξ},
如果 y n ( x ) y_n(x) yn(x)为最优解,则 y n ( x ) − x y_n(x)-x yn(x)x表示最优的初始购买量。

文章证明,如果持有成本和缺货成本是线性的,则每个时期中的最优策略总是 (s, S) 类型。

首先回顾单时期问题,即 n = 1 n=1 n=1 C 0 ≡ 0 C_0\equiv 0 C00, 而单时期问题本质上是一个微积分问题。

  • 最简单的情况是订购成本是线性的,即 c ( z ) = c ⋅ z c(z)=c\cdot z c(z)=cz. 在这种情况下,单时期模型的最优策略通常由单个临界数 x ˉ \bar{x} xˉ定义如下:如果 x < x ˉ x<\bar{x} x<xˉ,则购买 x ˉ − x \bar{x}-x xˉx;如果 x > x ˉ x>\bar{x} x>xˉ,则不购买。在多时期问题中,类似的结果通常成立,最优策略由一系列临界数 x ˉ 1 , x ˉ 2 , ⋯ \bar{x}_1,\bar{x}_2,\cdots xˉ1,xˉ2,定义。这些结果成立的充分条件是 L ( y ) L(y) L(y)是凸的,而这一条件在持有和短缺成本均为原点处消失的凸递增函数时成立。
  • 当订购成本不再是线性时,文章集中讨论了最简单类型的非线性成本:


其中 K K K通常被描述为固定订购成本。在该订购成本下,单时期模型中的最优策略通常由一对临界数 (s, S) 定义:如果$ x<s $,则订购 ( S − x ) (S-x) (Sx);如果$ x>s ,则不订购。在单时期模型中有一些示例显示这样的策略并非最优。然而,如果持有和短缺成本为线性函数 ,则不订购。在单时期模型中有一些示例显示这样的策略并非最优。然而,如果持有和短缺成本为线性函数 ,则不订购。在单时期模型中有一些示例显示这样的策略并非最优。然而,如果持有和短缺成本为线性函数h(u)=h\cdot u 和 和 p(u)=p\cdot u , 或者更一般地,如果 , 或者更一般地,如果 ,或者更一般地,如果L(y)$是凸的,那么多时期模型的最优策略是 (s, S).

然而,即使假设持有和短缺成本是线性的,在原文章发表的60年代前,关于动态模型最优策略的文献也非常少。Bratten
表示(引用章节9),如果需求的密度是递减的,动态模型的最优策略是由一系列临界数对 ( s 1 , S 1 s_1, S_1 s1,S1), ( s 2 , S 2 s_2, S_2 s2,S2),… 定义。唯一的另一种结果由 Karlin提出,即如果 φ \varphi φ有一个单调似然比,当持有和缺货成本是线性的,且 c + h > α p c+h>\alpha p c+h>αp,那么最优策略也是同样的类型。这两个结果都相当限制性,前者要求一个递减的密度,后者对成本有严格的约束。

而文章证明了,如果持有成本和缺货成本是线性的,或更一般地,当 L ( y ) L(y) L(y)是凸的并且订购成本如上非线性时,则每个时期中的最优策略总是 (s, S) 类型,无需任何附加条件。

上述两个结果的推导基于对函数 G n ( y ) G_n(y) Gn(y)的研究,其中
G n ( y ) = c y + L ( y ) + α ∫ 0 ∞ C n − 1 ( y − ξ ) φ ( ξ ) d ξ , G_n(y)=cy+L(y)+\alpha\int_0^{\infty}C_{n-1}(y-\xi)\varphi(\xi)\mathrm{d}\xi, Gn(y)=cy+L(y)+α0Cn1(yξ)φ(ξ)dξ,
当且仅当 存在 y > x y>x y>x, G n ( x ) > K + G n ( y ) G_n(x)>K+G_n(y) Gn(x)>K+Gn(y), 从 x x x订购是最优的;若要求从 x x x订购,则订购到 使 G n ( y ) G_n(y) Gn(y)最小的 y y y ( y > x ) (y>x) (y>x). 当假设 Bratten 的条件或 Karlin 的条件时, G n ( y ) G_n(y) Gn(y)递减到一个最小值然后再递增。现将最小满足情况的 y y y记作 S n S_n Sn且定义 s n s_n sn:
G n ( s n ) = G n ( S n ) + K , G_n(s_n)=G_n(S_n)+K, Gn(sn)=Gn(Sn)+K,
则由 ( s n , S n s_n,S_n sn,Sn) 定义的策略是最优的。然而,少量的数值计算就足以显示函数 G n G_n Gn并不总是符合次规律——实际上,可能有许多极大值和极小值。文章给出的证明思路是,尽管 G n G_n Gn可能有很多极大值和极小值,但其变化不会致使导致偏离 (s,S) 策略。

如果 L ( y ) L(y) L(y)是凸的,对于任一 a ≥ 0 a\geq 0 a0, 则 K + G n ( a + x ) − G n ( x ) − a G n ′ ( x ) ≥ 0 K+G_n(a+x)-G_n(x)-aG'_n(x)\geq 0 K+Gn(a+x)Gn(x)aGn(x)0.

文章分情况讨论了即时与延时交货的情况来证明上述不等式的成立。假设 G n ( x ) G_n(x) Gn(x)的变化如图1所示:系统会在区间 I 订购到点 S S S,在区间 II 不订购,在区间 III 订购到 S S S,在区间 IV 不订购,在区间 V 订购到 S ′ S' S,在区间 VI 不订购。但若以上不等式成立,则该曲线不能成立:若有让 x + a = S x+a=S x+a=S并且 x x x是在 III 中达到相对最大值的点,而对于这个值的 x x x, G n ′ = 0 G'_n=0 Gn=0并且以上不等式意味着 K + G n ( x ) − G n ( x ) ≥ 0 K+G_n(x)-G_n(x)\geq 0 K+Gn(x)Gn(x)0,这与图1矛盾。同样的论证可以应用于点 S ′ S' S. 接下来针对订购货物到达的时间分情况讨论。


图1. 在假设的订购方案下的期望成本变化

零交货延迟的情况

对于即时交货的情况,文章提出以下定义:

K ≥ 0 K\geq 0 K0并且 f ( x ) f(x) f(x)是一个可微分函数。如果对于所有 a > 0 a>0 a>0和所有 x x x,有
K + f ( a + x ) − f ( x ) − a f ′ ( x ) ≥ 0 , K+f(a+x)-f(x)-af'(x)\geq 0, K+f(a+x)f(x)af(x)0,
则称 f ( x ) f(x) f(x)是 K-凸的。如果不假设可微分性,则 K-凸性的定义为
K + f ( a + x ) − f ( x ) − a [ f ( x + b ) − f ( x ) b ] ≥ 0. K+f(a+x)-f(x)-a[\frac{f(x+b)-f(x)}{b}]\geq 0. K+f(a+x)f(x)a[bf(x+b)f(x)]0.
K-凸函数的简单属性包括:

  1. 0-凸性等同于普通凸性。
  2. 如果 f ( x ) f(x) f(x)是 K-凸的,那么对于所有 h h h, f ( x + h ) f(x+h) f(x+h)也是 K-凸的。
  3. 如果 f f f g g g分别是 K-凸和 M-凸的,且有 α \alpha α β \beta β为正,那么 α f + β g \alpha f+\beta g αf+βg ( α K + β M ) (\alpha K+\beta M) (αK+βM)-凸的。

运用数学归纳法,首先, G 1 ( x ) = c x + L ( x ) G_1(x) = cx+L(x) G1(x)=cx+L(x)显然是 K-凸的(0-凸)。

其次,假设 G 2 , … , G n G_2,\ldots,G_n G2,,Gn也是 K-凸的。根据 G n ( x ) G_n(x) Gn(x)的定义,为证明 G n + 1 ( x ) G_{n+1}(x) Gn+1(x)的 K-凸性,则只需证明 ∫ 0 ∞ C n ( x − ξ ) φ ( ξ ) d ξ \int_{0}^{\infty}C_n(x-\xi)\varphi(\xi)\mathrm{d}\xi 0Cn(xξ)φ(ξ)dξ是 K-凸的,而通过上述属性2和3,足够证明 C n ( x ) C_n(x) Cn(x)是 K-凸的,推导如下。

首先,如果 S n S_n Sn G n ( x ) G_n(x) Gn(x)的绝对最小值,并且已定义 s n s_n sn为满足 K + G n ( S n ) = G n ( s n ) K+G_n(S_n)=G_n(s_n) K+Gn(Sn)=Gn(sn)的$ x<S_n ,那么最优策略是:当 ,那么最优策略是:当 ,那么最优策略是:当 x<s_n 则订购到 则订购到 则订购到S_n$,否则不订购。因此,

  • 当$ x<s_n 时, 时, 时,C_n(x)=K+c(S_n-x)+C_n(S_n)=K-cx+G_n(S_n)$;
  • 当$ x>s_n 时, 时, 时,C_n(x) = -cx+G_n(x)$.

情况 1 x > s n x>s_n x>sn. 在这个区域, C n ( x ) C_n(x) Cn(x)等于线性函数加上K-凸函数,因此是 K-凸的。

情况 2:$ x<s_n<x+a $. 在这种情况下,
K + C n ( x + a ) − C n ( x ) − a C n ′ ( x ) = K + C n ( x + a ) − C n ( x ) + a c K+C_n(x+a)-C_n(x)-aC'_n(x)=K+C_n(x+a)-C_n(x)+ac K+Cn(x+a)Cn(x)aCn(x)=K+Cn(x+a)Cn(x)+ac
为正,因为
C n ( x ) = min ⁡ y > x { K + c ( y − x ) + L ( y ) + α ∫ 0 ∞ C n − 1 ( y − ξ ) φ ( ξ ) d ξ } C_n(x)=\min_{y>x}\{K+c(y-x)+L(y)+\alpha\int_0^{\infty}C_{n-1}(y-\xi)\varphi(\xi)\mathrm{d}\xi\} Cn(x)=y>xmin{K+c(yx)+L(y)+α0Cn1(yξ)φ(ξ)dξ}
≤ K + c a + K ( x + a ) + α ∫ 0 ∞ C n − 1 ( x + a − ξ ) φ ( ξ ) d ξ \leq K+ca +K(x+a)+\alpha\int_0^{\infty}C_{n-1}(x+a-\xi)\varphi(\xi)\mathrm{d}\xi K+ca+K(x+a)+α0Cn1(x+aξ)φ(ξ)dξ
= K + c a + C n ( x + a ) . =K+ca+C_n(x+a). =K+ca+Cn(x+a).
情况 3:$ x + a < s_n . 在这个区域, . 在这个区域, .在这个区域, C_n(x) $是线性的,因此是 K-凸的。归纳完成,并证明了在这一部分考虑的情况下 (s, S) 策略的最优性。

交货延迟的情况

当存在交货延迟时,最优策略的性质很大程度上取决于是超额需求被积压还是被加急处理。如果超额需求被积压,那么已知最优策略是现有库存加上已订购但尚未交付的库存的函数;如果超额需求被加急处理,最优策略则将复杂得多——文章仅考虑积压为满足需求的情况。

令交货延迟记作 λ \lambda λ,即在时期开始时下的订单在 λ \lambda λ个时期后的开始时交付。考虑 n n n个时期,令 x x x表示当前库存, x 1 x_1 x1表示下个时期开始时要交付的库存, x j x_j xj表示 j j j个时期后交付的库存,其中 j = 1 , 2 , … , λ − 1 j=1,2,\ldots,\lambda-1 j=1,2,,λ1. 令 C n ( x , x 1 , … , x λ − 1 ) C_n(x,x_1,\ldots,x_{\lambda-1}) Cn(x,x1,,xλ1)表示最小期望成本,则
C n ( x , x 1 , … , x λ − 1 ) = min ⁡ z ≥ 0 { c ( z ) + L ( x ) + α ∫ − ∞ C n − 1 ( x + x 1 − ξ , x 2 , … , z ) φ ( ξ ) d ξ } . C_n(x,x_1,\ldots,x_{\lambda-1})=\min_{z\geq 0}\{c(z)+L(x)+\alpha\int_-^{\infty}C_{n-1}(x+x_1-\xi,x_2,\ldots,z)\varphi(\xi)\mathrm{d}\xi\}. Cn(x,x1,,xλ1)=z0min{c(z)+L(x)+αCn1(x+x1ξ,x2,,z)φ(ξ)dξ}.
接下来将证明,如 L ( x ) L(x) L(x)是凸的,并且购买成本如前文所示非线性,最优策略可由 S n S_n Sn s n s_n sn描述:如果$ x+x_1+\ldots +x_{\lambda-1}>s_n 则不订购;如果 则不订购;如果 则不订购;如果 x+x_1+\ldots +x_{\lambda-1}<s_n 则订购到 则订购到 则订购到S_n$.

根据 C n ( x , x 1 , … , x λ − 1 ) C_n(x,x_1,\ldots,x_{\lambda-1}) Cn(x,x1,,xλ1)定义,对于 n ≥ λ n\geq \lambda nλ
C n ( x , x 1 , x 2 , … , x λ − 1 ) C_n(x,x_1,x_2,\ldots,x_{\lambda-1}) Cn(x,x1,x2,,xλ1)
= L ( x ) + α ∫ 0 ∞ L ( x + x 1 − ξ ) φ ( ξ ) d ξ + ⋯ =L(x)+\alpha\int_0^{\infty}L(x+x_1-\xi)\varphi(\xi)\mathrm{d}\xi + \cdots =L(x)+α0L(x+x1ξ)φ(ξ)dξ+
+ α λ − 1 ∫ 0 ∞ ⋯ ∫ 0 ∞ L ( x + ⋯ + x λ − 1 − ∑ i = 1 λ − 1 ξ i ) φ ( ξ 1 ) ⋯ φ ( ξ λ − 1 ) d ξ ⋯ d ξ λ − 1 + f n ( x + x 1 + ⋯ + x λ − 1 ) , +\alpha^{\lambda-1}\int_0^{\infty}\cdots\int_0^{\infty}L(x+\cdots +x_{\lambda-1}-\sum_{i=1}^{\lambda-1}\xi_i)\varphi(\xi_1)\cdots\varphi(\xi_{\lambda-1})\mathrm{d}\xi\cdots\mathrm{d}\xi_{\lambda-1} +f_n(x+x_1+\cdots+x_{\lambda-1}), +αλ100L(x++xλ1i=1λ1ξi)φ(ξ1)φ(ξλ1)dξdξλ1+fn(x+x1++xλ1),
其中 f n ( u ) f_n(u) fn(u)满足
f n ( u ) = min ⁡ x ≥ 0 { c ( z ) + α λ ∫ 0 ∞ ⋯ ∫ 0 ∞ L ( u + z − ∑ i = 1 λ ξ i ) φ ( ξ 1 ) ⋯ φ ( ξ λ ) d ξ 1 ⋯ d ξ λ } . f_n(u)=\min_{x\geq0}\{c(z)+\alpha^{\lambda}\int_0^{\infty}\cdots\int_0^{\infty}L(u+z-\sum_{i=1}^{\lambda}\xi_i)\varphi(\xi_1)\cdots\varphi(\xi_{\lambda})\mathrm{d}\xi_1\cdots\mathrm{d}\xi_{\lambda}\}. fn(u)=x0min{c(z)+αλ00L(u+zi=1λξi)φ(ξ1)φ(ξλ)dξ1dξλ}.
因此若令 y = u + z y=u+z y=u+z,即相当于 L ( y ) L(y) L(y)
α λ ∫ 0 ∞ ⋯ ∫ 0 ∞ L ( u + z − ∑ i = 1 λ ξ i ) φ ( ξ 1 ) ⋯ φ ( ξ λ ) d ξ 1 ⋯ d ξ λ } \alpha^{\lambda}\int_0^{\infty}\cdots\int_0^{\infty}L(u+z-\sum_{i=1}^{\lambda}\xi_i)\varphi(\xi_1)\cdots\varphi(\xi_{\lambda})\mathrm{d}\xi_1\cdots\mathrm{d}\xi_{\lambda}\} αλ00L(u+zi=1λξi)φ(ξ1)φ(ξλ)dξ1dξλ}
替换。然而,如果 L ( y ) L(y) L(y)是凸的,那么替换它的上式也是凸的,这就足以完成证明。于是可以得出结论,如果持有成本和缺货成本是线性的,则每个时期中的最优策略总是 (s, S) 类型。

参考文献​


ARROW, K. J., T. HARRIS, and J. MARSCHAK. “Optimal Inventory Policy,” Econometrica, 19(1951), 250-72.​

ARROW, K. J., S. KARLIN, and H. SCARF. Studies in the Mathematical Theory of Inventory and Production, Stanford, Calif.: Stanford University Press, 1958.​

BELLMAN, R., I. GLICKSBERG, and O. GROSS. “On the Optimal Inventory Equation,” Management Science, 2(1955), 83-104.​

DVORETZKY, A., J. KIEFER, and J. WOLFOWITZ. “The Inventory Problem, I. Case of Known Distributions of Demand,” Econometrica, 20(1952), 187-222.​

DVORETZKY, A., J. KIEFER, and J. WOLFOWITZ. “On the Optimal Character of the (S,s) Policy in Inventory Theory,” Econometrica, 21(1953), 586-96.

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.hqwc.cn/news/670529.html

如若内容造成侵权/违法违规/事实不符,请联系编程知识网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

什么是PXE

文章目录 在局域网内搭建PXE服务器PXE 启动组件PXE的优点实验一、搭建PXE服务器&#xff0c;实现远程部署CentOS系统环境准备server关闭防火墙安装组件准备 Linux 内核、初始化镜像文件及PXE引导文件配置启用TFTP 服务配置启动DHCP服务准备CentOS 7 安装源配置启动菜单文件 Cli…

面试笔记——JVM组成

基本介绍 JVM: Java Virtual Machine Java程序的运行环境&#xff08;java二进制字节码的运行环境&#xff09; 使用JVM的好处&#xff1a; 一次编写&#xff0c;到处运行自动内存管理&#xff0c;垃圾回收机制 JVM的组成及运行流程&#xff1a; 程序计数器 程序计数器&a…

手动交互式选点提取三维点云轮廓边界线 附python代码

一种新的三维点云轮廓边界提取方案: 1 手动选择一个边界或者其附近的点 2 自动搜索临近区域,并找到附近的平面和进行平面分割 3 提取平面的交点 4 该交点就是点云的轮廓边界点,把它往两边延展,就是完整的点云轮廓边界 import open3d as o3d import numpy as np import …

[激光原理与应用-92]:振镜的光路图原理

目录 一、振镜的光路 二、振镜的工作原理 2.1 概述 2.2 焊接头 2.3 准直聚焦头-直吹头 2.4 准直聚焦头分类——按应用分 2.4.1 准直聚焦头分类——功能分类 2.4.2 准直聚焦头镜片 2.4.3 振镜焊接头 2.4.4 振镜分类&#xff1a; 2.4.5 动态聚焦系统演示&#xff08;素…

AtCoder Regular Contest 176(ARC176)A、B

题目&#xff1a;AtCoder Regular Contest 176 - tasks 官方题解&#xff1a;AtCoder Regular Contest 176 - editorial 参考&#xff1a;atcoder regular 176 (ARC176) A、B题解 A - 01 Matrix Again 题意 给一个nn的方格&#xff0c;给出m个坐标(x,y)m&#xff0c;在方格中…

一个简单的仓库出入库管理软件的流程是什么样的?有哪些功能?

身为仓库文员&#xff0c;我深知仓库管理对于公司运营的重要性。仓库是公司物资的中转站&#xff0c;其管理的好坏直接关系到公司的运营效率和成本控制。然而&#xff0c;传统的仓库管理方式往往存在着效率低下、易出错等问题&#xff0c;为了解决这些问题&#xff0c;我们需要…

stm32f103zet6_DAC_1_介绍

STM32微控制器系列的DAC&#xff08;数字到模拟转换器&#xff09;功能是其片上外设之一&#xff0c;用于将数字信号转换为模拟信号。DAC在许多应用中都非常有用&#xff0c;例如音频输出、模拟信号生成、闭环控制系统中作为模拟输出等。 STM32微控制器的DAC功能特点包括&…

Sy9-dhcp/dns服务配置

前言、 课程需要&#xff08;进入服务器综合配置使用阶段了&#xff09;&#xff0c;这里还是沿用桌面版的ubuntu来配置dhcp和dns&#xff0c;这里updated了新的版本。2024.5 server端环境&#xff1a; Win11VMS&#xff1a;192.168.105.1192.168.105.128 &#xff08;ubuntu…

搜索算法系列之四(斐波那契)

以下算法被验证过&#xff0c;如有什么问题或有补充的欢迎留言。 前言 斐波那契数列&#xff0c;又称黄金分割数列&#xff0c;是由意大利数学家&#xff08;Leonardo Fibonacci&#xff09;在1202年提出的。这个数列的递推关系是F(0)1&#xff0c;F(1)1&#xff0c;F(n)F(n-…

Vscode 总是发现不到包package

Vscode 右下角&#xff0c;选择3.8.6就不再报错&#xff0c;其他最新版本的3.11.7就不行。

嵌入式学习67-C++(多线程,自定义信号合槽,串口通信)

知识零碎&#xff1a; QmessageBox 报错提示框 GPS传感器获取到的 经纬度信息并不是真实的物理坐标&#xff0c;还需要转换 signals&#xff1a; …

NEO 学习之session7

文章目录 选项 A&#xff1a;它涉及学习标记数据。 选项 B&#xff1a;它需要预定义的输出标签进行训练。 选项 C&#xff1a;它涉及在未标记的数据中寻找模式和关系。 选项 D&#xff1a;它专注于根据输入-输出对进行预测。 答案&#xff1a;选项 C 描述了无监督学习的本质&am…