延期支付与订购量相关的易腐产品库存决策模型
2018-04-02韦才敏李忠萍林小苹
韦才敏,李忠萍,林小苹
(1.汕头大学理学院,广东汕头 515063)
(2.汕头大学数字信号与图像处理技术重点实验室,广东汕头 515063)
1 引言
当供应链成员企业面临资金约束,而银行信贷紧缩、资本市场不发达,企业如何从成员企业间获得内部融资来实现最优的运营决策是供应链融资管理研究的热点问题[1].内部融资包括提前付款融资模式与延期支付融资模式.针对延期支付而言,其形式有两种:在供应商为强势企业的情况下,供应商为了激励零售商的购买行为,往往会给予零售商一些优惠策略[2];在零售商为强势企业的情况下,其通过占用上游企业供应商的资金,扩大市场规模来获得收益[1].
已有很多国内外学者对有关延期支付的问题进行研究,如Goyal[3]最早构建了允许延期支付的EOQ模型,随后又有不少科研工作者对该模型进行扩展.在允许延期支付条件下,Chung[4]发现了一个确定EOQ下的定理.在两阶段的信贷策略下,基于EPQ模型框架,Huang[5]获得了零售商的最优补货决策.在供应商提供的延期支付与零售商订购量相关时,Chen[6]等人建立了零售商的EOQ模型,获得了最优订购量及订购周期.基于EPQ模型框架,针对变质产品,在零售商提前还款及供应商给予其现金折扣优惠策略下,李明芳[7]等人建立了零售商的库存决策模型,获得了限制条件下零售商的最优决策,然而,其是对变质很低的情形下的结果,在实际应用中有一定的限制.
然而,上述的研究都假定需求率是确定的常数,而实际中需求往往和零售价格相关,此种情况下的延期支付问题,也有学者对其进行研究.在数量折扣与延期支付固定的条件下,Sana[8]等人考虑了需求依赖于库存水平、时间及价格等背景下,其获得了零售商的最优决策解.针对随机需求情形,在两水平的信贷策略下,Dye[9]等人建立了确定的EOQ模型,通过粒子群优化算法获得了最优零售价格、补货数及补货策略.针对需求依赖于价格的情形,李明芳[10]等人在供应商允许延期支付和现金折扣的优惠策略下,建立了联合库决策模型,获得了零售商的最优订购周期、最优零售价格及最优付款时间.梁培培[11]等人建立了需求基于瞬时库存的延期支付库存决策模型,并验证了协调的必然性.基于延期支付,Shukla[12]等人建立了需求依赖于时间的EOQ模型,获得了零售商最优订购策略.Shukla[13]等人建立了需求依赖于库存水平的EOQ模型,获得了最优订购周期、订购量,实现了相关成本的最小化.
在以上文献中,虽然考虑了需求率和零售价格、时间、库存水平等相关下的延期支付库存决策问题,但大多在不允许缺货假设下进行的,且很少考虑延期支付与订购量相关情形下的库存决策模型.为了使模型更接近实际情况,在以上的文献基础之上,针对单周期内市场的需求依赖于零售价格、允许缺货发生的情形,建立了延期支付与订购量相关的易腐产品的库存决策模型.讨论了如何制定最优缺货时间、订购周期及零售价格,来实现零售商年利润最大化,并给出了求解该最优决策问题的内点算法,为零售商在库存管理方面提供了理论支持.并通过数值实验分析参数对最优决策的影响.
2 符号与假设
2.1 符号
模型中所使用的符号如表1所示.
表1:符号说明
2.2 假设
在模型的构建中所做的假设如下:
(1)时间期限无限;
(2)交货时间忽略不计;
(3)允许缺货发生,且缺货量完全延期供给;
(4)仅考虑单一易腐产品,如食品、水果、蔬菜等;
(5)单位时间市场的需求量D(p)=a−bp,其中a为市场的潜在需求量,b>0为客户对价格的敏感度;
(6)如果零售商的订购量Q+S≥Qd(Td≤t1),贷款全部获得延期支付,即零售商允许在延期期限M 偿还进货费用c(Q+S);如果零售商的订购量Q+S<Qd(Td>t1),贷款部分延期支付,可以延期到期限M 的费用为αc(Q+S),其余费用在订购产品时立即付清;
(7)在延期期限M 之前零售商将销售收入存入银行获得利息收益,在延期期限M 之后如果零售商不偿还贷款,供应商将向其收取利息,且Ik>Ie.
3 数学模型
本节首先获得t时刻零售商库存水平的表达式I(t),然后分情况讨论零售商单位时间的平均总利润表达式.类似[14],I(t)满足以下微分方程
满足条件I(t1)=0,可得微分方程的解为
且满足
由(3.2)与(3.3)式,可得订购周期内零售商库存水平随时间变化的大致图形如下.
图1:零售商库存水平变化图
零售商单位时间的平均总利润为ATP=ATP(t1,T,p)=(销售收入−总费用)/订购周期.总费用=订购费用+库存费用+腐败损失费用+采购费用+短缺费用+利息收入−利息支出.在一个周期内总销售收入为pD(p)T,订购费用为A,库存费用为腐败损失费用为c[Q−D(p)t1],采购费用为c(Q+S),短缺费用为其中在贷款延期到达之前零售商将销售产品所获得的收益存入银行获得利息收入,而在延期期限到达之后如果零售商不偿还贷款,则其需要支付高额利息费用,而对于订购周期内利息收入与利息支出,根据贷款全部或部分延期支付分为以下两种情形.
情形1Td≤t1.
情形2 Td>t1.
图2:利息收入及支付示意图
综合情形1与情形2,可得零售商单位时间的平均总利润为
其中
4 最优决策分析
本节针对延期支付与订购量相关的易腐产品库存决策模型,且允许缺货发生,缺货完全延期供给.以零售商单位时间内的平均总利润最大化为目标,获得最优零售价格、最优订购周期及最优缺货时间即最优化问题为以下形式
由于(4.1)式是含有约束条件的非线性优化问题,由于计算的复杂性,因此,很难求得该问题的解析解.接下来,在限制一些条件情况下,获得该问题的解析解;然后,又提出一种算法来获取该问题最优数值解.
4.1 易腐率很小时的最优决策
首先,通过一个定理来说明优化问题(4.1)的最优解的存在性.然后,当产品易腐率很低时,运用泰勒展开式定理,在固定T时获得最优零售价格及最优缺货时间表达式,实现了零售商单位时间的平均总利润最大化.
接下来,通过一个定理来说明优化问题(4.1)的最优解的存在性.
定理4.1存在最优零售价格、最优订购周期及最优缺货时间使得ATPi(t1,T,p)达到最大化.
证由于决策目标函数ATPi(t1,T,p)对应的海塞矩阵为可以判定其对角线元素H11,H22,H33均小于零,即证.其中
然后,在产品易腐率很低及固定T的情况下,通过一个定理给出订购周期内零售商的最优零售价格及最优缺货时间的解析解,并实现了零售商单位时间的平均总利润的最大化.
定理4.2当T固定,且B2−4AC≥0时,存在最优零售价格及最优缺货时间使得零售商单位时间的平均总利润ATPi(t1,T,p)达到最大化,且
其中
证在固定T的情况下,决策函数TPi(t1,T,p)关于t1,p求一阶偏导数,可得
令(4.4)与(4.5)式等于零,可得
由(4.8)与(4.6)式,可得(4.3)式,对应的海塞矩阵为
4.2 易腐率比较大时的最优决策
令x=(p,T,t1),则方程(3.4)中的目标函数可转化为
其中 gik≥ 0(i=1,···,4;k=1,···,4)为不等式约束条件;fi(x)= −ATPi(t1,T,p),则公式(3.4)可变为
通过minf(x)构造罚函数,即增广目标函数的无约束子问题为
为了获得优化问题(4.9)的数值解,本节采用内点算法[15−17]进行计算,该算法的具体迭代流程图如下.
图3:算法流程图
注4.1以上算法中,罚因子选择准则为选取数列{τk},满足以下条件
(1)正数列;
(2)单调递减趋于零.
注4.2罚因子的选择对该算法有一定的影响,一般来说不能太大也不能太小,当罚因子过大,可能会使问题(4.10)远离问题(4.9);当罚因子过小,即τk→0,惩罚项→0·∞形式,则会给增广目标函数的求解带来计算上的困难.
假设某物品的周需求量为D=100000−5000p千克,每千克该物品的购买成本为c=3元,每次固定订购费用为A=1000元,每周期每千克该物品的存储费用为h=0.3元,每周期每千克该物品的存储费用为s=2元,该物品的变质率为θ=0.01,延期期限为M=0.5周,年利息收益率为Ie=0.09,年利息支付率为Ie=0.13,贷款部分延期支付的比率α=0.75(每周期按365天计算).运用内点算法,并借助Matlab2014a软件进行求解,可得最优零售价格、最优订购周期及最优缺货时间的数值解,及零售商单位时间的平均总利润最大值如下表2.
表2:各种情形下对应的零售商的最优决策
可得,不同情形下,零售商的最优订购量为
分析可知在贷款全部获得延期支付的最小订购量取对应值下,零售商的最优缺货时间、最优订购周期、最优零售价格及订购周期内年最大利润,如下表3.
表3:贷款全部获得延期支付最小订购量Qd对缺货时间、订购周期、零售价格及零售商利润影响
5 数值分析
本节通过数值实验进一步分析采购费用、库存费用、缺货费用及腐蚀率对最优缺货时间、最优订购周期、最优零售价格及零售商年最大利润的影响,如下表4–7.假设所对应的参数值和4.2节中的相同,不失一般性,仅M ≤t1<T(Q+S≥Qd)情形下进行讨论.
表4:采购费用对缺货时间、订购周期、零售价格及零售商利润的影响
表5:库存成本对缺货时间、订购周期、零售价格及零售商利润的影响
表6:腐蚀率对缺货时间、订购周期、零售价格及零售商利润的影响
由表4、5、6可知,最优缺货时刻随着单周期内固定采购费用、库存成本及腐蚀率的增加而不变;最优订购周期及最优零售价格随着单周期内固定采购费用、库存成本及腐蚀率的增加而增加;零售商年利润随着单周期内固定采购费用、库存成本及腐蚀率的增加而减少.
表7:缺货成本对缺货时间、订购周期、零售价格及零售商利润的影响
由表7可知,最优缺货时刻随着缺货成本的增加而不变;最优零售价格随着单周期内缺货成本的增加而增加;最优订购周期及零售商年利润随着缺货成本的增加而减少.
根据以上的数值分析,可以得到以下管理见解:
第一 单周期固定采购费用、库存成本及腐蚀率越少,缺货成本越多,则缺货时间越短.
第二 单周期固定采购费用、库存成本及腐蚀率越少,缺货成本越多,则订购周期越短.
第三 单周期固定采购费用、库存成本、腐蚀率及缺货成本越多,则零售价格越大.
第四 单周期固定采购费用、库存成本、腐蚀率及缺货成本越多,则零售商年利润越小.
6 结论
当零售商面临资金约束,零售商很难从银行等外部渠道获得融资时,如何实现内部融资是供应链金融研究的热点.因此,供应商对零售商提供延期支付优惠策略来刺激其订购行为,且延期支付的期限与订购量相关.首先,针对易腐产品,且允许缺货发生、零售价格影响市场需求量情形下,建立了相应的库存决策模型.其次,在易腐率很小的情况下,对目标函数进行二阶泰勒展开,并运用微积分及矩阵论知识,在订购周期固定的条件下,获得了的零售商最优缺货时间、最优零售价格及年平均最大利润的解析解.然后,在易腐率比较大的情况下,借助内点算法,获得了零售商最优决策的数值解.最后,通过数值实验分析单周期内固定采购费用、库存成本、腐蚀率及缺货成本对零售商最优决策的影响,可得单周期内固定采购费用、库存成本、腐蚀率及缺货成本越少,对零售商和客户越有利.
该模型在实际生活中有如下几方面的应用价值:①供应商可以根据日常中的具体情况,通过改变贷款全部获得延期支付的最小订购量来调控零售商的订购量.②当产品发生缺货时,零售商可以确定适当的缺货时间、订购周期及零售价格来实现其年利润的最大化,为零售商在库存管理方面提供了理论支持.③对易腐产品库存管理方面的问题,借助文中的算法可以很容易获取零售商的最优决策.同时,本模型可以进一步扩展,如解决多阶段多产品库存问题.
[1]王文利.需求不确定下的供应链融资研究―基于报童模型视角[M].北京:中国金融出版社,2014.
[2]李明芳,王道平.延期支付期限与订货量相关情形下的EOQ模型[J].系统管理学报,2011,20(4):389–397.
[3]Goyal S K.Economic order quantity under conditions of permissible delay in payments[J].J.Oper.Res.Soc.,1985,36(4):335–338.
[4]Chung K J.A theorem on the determination of economic order quantity under conditions of permissible delay in payments[J].Comp.Oper.Res.,1998,25(1):49–52.
[5]Huang Y.Optimal retailer’s replenishment decisions in the EPQ model under two levels of trade credit policy[J].Euro.J.Oper.Res.,2007,176(3):1577–1591.
[6]Chen S,C´ardenas-Barr´on L E.Teng J.Retailer’s economic order quantity when the supplier offers conditionally permissible delay in payments link to order quantity[J].Intern.J.Prod.Econ.,2014,155:284–291.
[7]李明芳,王道平,李锋.基于现金折扣和延期支付条件下变质产品的补货策略[J].管理评论,2011,23(4):122–128.
[8]Sana S S,Chaudhuri K S.A deterministic EOQ model with delays in payments and price-discount offers[J].Euro.J.Oper.Res.,2008,184(2):509–533.
[9]Dye C,Ouyang L.A particle swarm optimization for solving joint pricing and lot-sizing problem with fluctuating demand and trade credit financing[J].Comput.Indus.Engin.,2011,60(1):127–137.
[10]李明芳,王道平.延期支付条件下零售商的最优售价和订货策略[J].运筹与管理,2011,20(06):59–65.
[11]梁培培,孙延明.需求基于瞬时库存的延期支付策略设计[J].统计与决策,2015,30(06):41–43.
[12]Shukla H S,Shukla V and Yadav S K.Optimal ordering policies in the EOQ(economic order quantity)model with time-dependent demand rate under permissible delay in payments[J].Intern.J.Modern Engin.Sci.,2015,4(1):1–3.
[13]Shukla H S,Tripathi R P,Siddiqui A.EOQ model with inventory level dependent demand rate under permissible delay in payments with cash discount[J].Indian J.Sci.Tech.,2015,8(28).
[14]周永务,王圣东.库存控制理论与方法[M].北京:科学出版社,2009.
[15]马昌凤.最优化方法及其Matlab程序设计[M].北京:科学出版社,2009.
[16]雍龙泉,邓方安,陈涛.单调线性互补问题的一种内点算法[J].数学杂志,2009,29(5):681–686.
[17]龚小玉,张明望.非单调线性互补问题的高阶宽领域内点算法[J].数学杂志,2009,29(2):217–223.