田志勇，周 麗，袁瑞萍

（北京物資學(xué)院 信息學(xué)院，北京 101149）

考慮收益波動性的風(fēng)險厭惡型報童問題

田志勇，周 麗，袁瑞萍

（北京物資學(xué)院 信息學(xué)院，北京 101149）

將報童收益波動性融入目標(biāo)函數(shù)，以價格和庫存量作為決策變量，針對單期、加法形式的隨機(jī)需求函數(shù)和風(fēng)險厭惡型報童建立模型。分別利用順序優(yōu)化和同時優(yōu)化的方法對模型進(jìn)行研究，分析了模型最優(yōu)解與風(fēng)險系數(shù)和剩余庫存處理單價等參數(shù)的關(guān)系，利用彈性概念描述了模型最優(yōu)解存在的條件及實際意義。研究表明兩種方法的最優(yōu)解具有相同的性質(zhì)—價格和庫存量都是風(fēng)險系數(shù)的減函數(shù)、是剩余庫存處理單價的增函數(shù)。關(guān)于最優(yōu)解存在條件，同時優(yōu)化的局部最優(yōu)解要求最為寬松，全局最優(yōu)解要求最嚴(yán)格，而順序優(yōu)化最優(yōu)解的要求則介于兩者之間。這些研究結(jié)論為企業(yè)辨析決策環(huán)境、進(jìn)行有效決策提供了有益的參考。

報童模型；收益波動性；風(fēng)險厭惡型；彈性

1 引言

報童問題的經(jīng)典模型是以報童期望收益為目標(biāo)函數(shù)、庫存水平和價格為決策變量，研究其最優(yōu)化問題[1-2]。一般地，單期報童收益包括三部分：報童銷售收入、剩余存貨價值和潛在銷售損失。對于需求確定性問題來說，無論目標(biāo)函數(shù)包含哪部分，問題的抽象模型不會出現(xiàn)本質(zhì)差異。而對于包含隨機(jī)性問題來說，報童收益三部分的數(shù)學(xué)表現(xiàn)和性質(zhì)有較大差異，會對模型分析的復(fù)雜程度、最優(yōu)解特征甚至結(jié)論產(chǎn)生重要影響。多數(shù)的研究，如xu等[3-4]、Dogan[5]、Jadidi等[6]、Khouja[7]，目標(biāo)函數(shù)包含了三部分。少數(shù)研究則只包含部分收益，如Gerchak等[8]包含第二、三部分，Chua等[9]和丁小東等[10]只考慮了第一、二部分，而Javier等[11]則僅考慮了第一部分。采用收益的哪部分作為目標(biāo)函數(shù)，需要根據(jù)研究的目的、進(jìn)展和應(yīng)用場景而定。但一般來說，采用部分報童收益的研究多是一種新的研究應(yīng)用的開始，隨著研究的深入，會逐步將報童收益的三個部分包含進(jìn)來。

但僅以期望收益作為模型優(yōu)化的目標(biāo)函數(shù)，與報童決策的真實情景仍有較大差距。因為收益的波動顯然也是報童決策時要考量的因素，在某種情況下，甚至可能會成為影響報童決策的關(guān)鍵因素。比如，著名的6σ管理，其核心就是“不是關(guān)于平均數(shù)的問題，而是關(guān)于方差或波動的問題”[12]。對風(fēng)險的度量主要有三種[13]：均值方差類應(yīng)用、在險價值VaR和條件風(fēng)險價值CVaR。目前，VaR和CVaR應(yīng)用于報童風(fēng)險監(jiān)控的研究文獻(xiàn)比較多，如文獻(xiàn)[14]、[15]、[16]和[17]等。但VaR和CVaR均是測度和控制報童收益處于風(fēng)險（達(dá)到或達(dá)不到預(yù)定點）的概率，并未直接考量收益的波動性。而均值方差的研究文獻(xiàn)則集中在需求或供應(yīng)層面，直接應(yīng)用在成本或收益層面的文獻(xiàn)很少。需求波動雖在一定程度上可以看作利用需求概率函數(shù)間接地考量了收益波動（方差），但由于需求的波動性需要通過一定的函數(shù)關(guān)系傳遞到收益層面上，這種函數(shù)關(guān)系往往比較復(fù)雜，必然會導(dǎo)致直接在收益層面研究波動性在技術(shù)上更為復(fù)雜。但需求-收益函數(shù)關(guān)系使兩個層面的波動性呈現(xiàn)出什么樣的不同特征？對報童決策又將產(chǎn)生什么樣的影響？這些都需要針對收益的均值方差專門研究。據(jù)作者所知，文獻(xiàn)[11]是這一領(lǐng)域第一篇直接在收益層面考慮波動性的研究，但該文獻(xiàn)考慮的報童收益還僅僅是報童銷售收入，也即前面提及的第一部分。本文在以下方面進(jìn)行了擴(kuò)展：（1）目標(biāo)函數(shù)中增加剩余庫存價值，即包含了報童收益的第一和第二部分；（2）對順序優(yōu)化的兩種情況進(jìn)行了較系統(tǒng)地分析；（3）從本質(zhì)上分析順序優(yōu)化和同時優(yōu)化之間的區(qū)別與聯(lián)系。

2 模型構(gòu)建

2.1 需求函數(shù)

本文采用需求函數(shù)的加法形式

D表示需求,是價格p和隨機(jī)因素θ的函數(shù)。y(p)=a-bp為確定性需求函數(shù)，p為價格,a＞0為市場規(guī)模，b＞0為需求對價格反饋系數(shù)。θ為隨機(jī)變量，取值范圍為[A,B]，且 A≤0,B≥0，期望 E(θ)=μ ，方差V(θ)=σ2,概率密度和概率分布函數(shù)分布為f(.)和F(.)，設(shè)概率分布函數(shù)的補(bǔ)函數(shù)為

2.2 目標(biāo)函數(shù)

設(shè)q為庫存水平或訂貨批量，s為剩余庫存處理單價。一般來說，有c＞s。在單期情況下，考慮剩余存貨價值的報童利潤函數(shù)為：

函數(shù)Θ(z)是期望存貨剩余量函數(shù)，E(z)是期望銷售量函數(shù)。V(z)是利潤函數(shù)方差中與z相關(guān)的部分，或者由于z與p的可分離性，在利潤函數(shù)的方差中剔除p后的變動部分。

對于綜合報童期望利潤及利潤方差的目標(biāo)函數(shù)，我們采用文獻(xiàn)[11]的形式，即兩者的線性和，如下：

式中λ是風(fēng)險系數(shù)，為利潤方差賦予的權(quán)重。由于本文僅考慮風(fēng)險厭惡型報童，所以λ＞0。這意味著λ越大，報童對收益的波動性越敏感，越傾向于減少波動。

將確定性需求函數(shù)y()p=a-bp代入式（1），將目標(biāo)函數(shù)進(jìn)一步明確為：

3 順序優(yōu)化

順序優(yōu)化適于特定的情景，如文獻(xiàn)[18]所描述的兩種情景：公司內(nèi)的銷售營銷部門在給定庫存水平z的情況下進(jìn)行價格p(z)決策，然后生產(chǎn)運營部門在已知p(z)的情況下進(jìn)行庫存水平的決策；或者與此過程相反，先由生產(chǎn)運營部門在給定價格p情況下進(jìn)行庫存水平z(p)決策，然后由銷售營銷部門進(jìn)行定價決策。

在進(jìn)行分析前，需先做如下假設(shè)：

(A1)意味著，對于風(fēng)險厭惡型報童來說，其對風(fēng)險（方差）的加權(quán)系數(shù)有一個限度，不能無限大。因為那樣地話，將導(dǎo)致方差權(quán)重過大，在決策時過度考慮風(fēng)險因素，而對期望利潤的考慮過輕。根據(jù)2.1部分的分析，(A2)意味著嚴(yán)格限制意義下的價格最大值與期望概念下價格最大值的距離要小于與成本的距離。

3.1 順序優(yōu)化I

順序優(yōu)化也可以分為兩種情景：情景I先由生產(chǎn)運營部門確定庫存水平z再由銷售營銷部門確定價格p；情景II先由銷售營銷部門確定價格p再由生產(chǎn)運營部門確定庫存水平z。

對于情景I，銷售營銷部門在z確定的情況下進(jìn)行價格p決策。其一階條件為

其二階條件可以保證目標(biāo)函數(shù)Σ是關(guān)于z的凹函數(shù)。所以可得最優(yōu)解p：

引理1：p*(z)是z的增函數(shù)；對于 ?z∈[A,B]，有p*(z)∈[c,pmax]。

對于第二部分，

第一個不等式是分子、分母分別取其可能的最大、最小值得到，第二個不等式是由假設(shè)(A2)得到。

可見，報童的風(fēng)險系數(shù)λ越小，剩余庫存對于報童的決策影響越小，報童越傾向于在最優(yōu)解處采取高價格。剩余庫存單價越高，報童越能容忍更多的剩余庫存量，就越傾向于在最優(yōu)解處采取更高的價格。

接著研究z對目標(biāo)函數(shù)的影響，得目標(biāo)函數(shù)對z的一階條件：

還需檢驗關(guān)于z的二階條件，如文獻(xiàn)[11]和[18]，先定義一個概念—銷售損失率對價格彈性。

定義：根據(jù)經(jīng)濟(jì)學(xué)彈性概念的一般定義，對于給定庫存水平q的情況下，銷售損失率對價格彈性為：

將上式代入銷售損失率對價格的彈性表達(dá)式，可得：

在本文中，由于考慮到處理剩余庫存的單價，我們將銷售損失率對價格的彈性定義為對“絕對價格”即=p-s的彈性，則上式變?yōu)椋?/p>

將在最優(yōu)點p*()z處的彈性表示為：

然后，有以下定理：

證明：

對z*關(guān)于兩個重要參數(shù)λ和s的關(guān)系，有如下引理：

引理3：z*關(guān)于λ單調(diào)減；z*關(guān)于s單調(diào)增。

根據(jù)式（5）應(yīng)用隱函數(shù)定理，可以分別求解z*關(guān)于λ和s的導(dǎo)數(shù)，易證。

3.2 順序優(yōu)化II

再分析順序優(yōu)化情景II：先確定價格p再確定庫存水平z。其決策過程仍是逆向歸納求解：生產(chǎn)經(jīng)營部門先對庫存水平z進(jìn)行優(yōu)化，然后銷售和營銷部門在已知z為p函數(shù)的情況下進(jìn)行定價p的優(yōu)化。

生產(chǎn)經(jīng)營部門的庫存水平優(yōu)化，得到目標(biāo)函數(shù)關(guān)于z的一階條件：

式（6）與式（5）形式相同，只是p的取值有差異。利用假設(shè)(A1)，容易驗證二階條件成立，確保了一階條件可以獲得最優(yōu)解z*()p。與定理1相比，本例最優(yōu)解的存在條件相對來說要寬松一些。

根據(jù)一階條件式（5）和隱函數(shù)定理分析z與p、λ及s的關(guān)系。

引理4：最優(yōu)解z*關(guān)于p單調(diào)增；關(guān)于λ單調(diào)減；關(guān)于s單調(diào)增。

證明過程同引理3。

雖然我們不能得到z*關(guān)于p的解析解，但仍可以將其視為p的函數(shù)，利用萊布尼茨鏈?zhǔn)角髮?dǎo)法則求解目標(biāo)函數(shù)對p的一階和二階導(dǎo)數(shù)。目標(biāo)函數(shù)關(guān)于p的一階條件為

證明：p的二階條件為：

關(guān)于p*與參數(shù)λ和s的關(guān)系，有以下引理：

引理5：最優(yōu)解p*關(guān)于λ單調(diào)減；關(guān)于s單調(diào)增。

證明過程同引理3、4。

比較順序優(yōu)化II的引理4和順序優(yōu)化I的引理1、3，兩種情況下的最優(yōu)值z*與p、λ和s的相關(guān)關(guān)系完全一致。再比較順序優(yōu)化II的引理5與順序優(yōu)化I的引理2，兩種情況下的最優(yōu)值p*與λ和s的相關(guān)關(guān)系完全一致。如前面所述，這絕非偶然。因為兩種情況的最優(yōu)解是相同的，均是目標(biāo)函數(shù)相對于p和z的一階條件組成的方程組的解。

4 同時優(yōu)化

根據(jù)多變量最優(yōu)化問題理論，同時優(yōu)化的最優(yōu)解也是由各變量的一階條件組成的方程組的解。不過，最優(yōu)解的存在需要檢驗海塞矩陣，根據(jù)海塞矩陣的正定性有定理3和4。

證明：需證明在點(p*,z*)處，海塞矩陣為半負(fù)定。即

同時關(guān)于z的二階偏導(dǎo)：

上式不等式成立是由于假設(shè)(A1)。

然后海塞行列式的值為：

證明：需證明海塞矩陣負(fù)定，其過程與定理3完全相同，區(qū)別僅為相應(yīng)的二階偏導(dǎo)數(shù)及海塞矩陣行列式值的取值點，本例中為任意點(p,z)而非駐點(p*,z*)。證畢。

為了便于比較幾種情況下的優(yōu)化，將定理1—4的條件和結(jié)論予以總結(jié)，見表1。

由表1，定理1和2的充分條件介于定理3和4之間，強(qiáng)于定理3但弱于定理4。從這個角度，可以把順序優(yōu)化的最優(yōu)解稱為“偏全局最優(yōu)解”，順序優(yōu)化I為“偏z全局最優(yōu)解”，順序優(yōu)化II為“偏p全局最優(yōu)解”。本質(zhì)上，四個定理得到的最優(yōu)解是相同的。只是定理1和定理2給出了該方程組的求解方法—代入消元法，定理1是先用函數(shù)p*(z)消去p得到z的值，然后代入p*(z)得到p；定理2是先用函數(shù)z*(p)消去z求解p，然后代入z*(p)得到z。定理3和4只是列出方程組，指出方程組的解即為最優(yōu)解，并未指出求解方法。

根據(jù)此處對三種方法、四個定理最優(yōu)點的比較分析，可知前面引理1和2的結(jié)論具有通用性。所以，同時優(yōu)化的最優(yōu)解(p*,z*)對相關(guān)參數(shù)λ和s的關(guān)系同引理1和2，此處不再贅述。

表1 順序優(yōu)化和同時優(yōu)化條件及結(jié)論

5 總結(jié)

本文研究了以包含收益波動的收益函數(shù)為目標(biāo)函數(shù)、以價格和庫存為決策變量的報童問題，其中包含隨機(jī)因素的需求函數(shù)是加法形式的顯函數(shù)。論文先對這種研究思路與以往以期望收益（成本）為目標(biāo)函數(shù)的區(qū)別與聯(lián)系進(jìn)行了分析，然后對文獻(xiàn)[11]進(jìn)行擴(kuò)展，在報童收益中增加了剩余庫存價值，針對風(fēng)險厭惡型報童分析了順序優(yōu)化和同時優(yōu)化的最優(yōu)解、最優(yōu)解相對于主要參數(shù)的比較靜態(tài)分析以及最優(yōu)解存在的充分條件，并對其實際意義進(jìn)行了解釋。由于需求函數(shù)的顯式性，得到一些相對較明確的結(jié)論，但由于隨機(jī)因素的概率密度函數(shù)和累計概率函數(shù)以隱函數(shù)形式給出，所以最優(yōu)解的明確性并不徹底。幸運的是，最優(yōu)解對主要參數(shù)的比較靜態(tài)分析的方向性關(guān)系還是明確的。本文在研究上進(jìn)行了拓展和深入，在對順序優(yōu)化和同時優(yōu)化進(jìn)行分析時，特別是兩種情景的順序優(yōu)化，都應(yīng)用了最優(yōu)解求解—最優(yōu)解對主要參數(shù)λ、s（另有p對z，或z對p）的比較靜態(tài)分析—二階條件判斷的分析框架一以貫之，用統(tǒng)一的框架較清晰、系統(tǒng)地開展了相應(yīng)研究。最后，通過對順序優(yōu)化和同時優(yōu)化兩種方法的比較分析，得到兩種方法間的本質(zhì)區(qū)別與聯(lián)系：順序優(yōu)化是同時優(yōu)化的局部最優(yōu)解，兩者的最優(yōu)解及最優(yōu)解的比較靜態(tài)分析結(jié)論是一致的，但二階判斷條件不同，順序優(yōu)化存在性條件的復(fù)雜性介于同時優(yōu)化的局部最優(yōu)和全局最優(yōu)之間，相當(dāng)于“偏全局最優(yōu)”。此外，對于風(fēng)險偏好型報童，最優(yōu)解及比較靜態(tài)分析與風(fēng)險厭惡相同或相似，但最優(yōu)解成立的充分條件—二階條件相對更為復(fù)雜，也可以得到與風(fēng)險厭惡型報童相對應(yīng)的二階條件。但由于該條件無論從形式上還是實際意義上，都不如風(fēng)險厭惡型報童情景中表現(xiàn)得簡潔直觀，所以本文不再對風(fēng)險偏好報童的情形進(jìn)行討論。此外，在目標(biāo)函數(shù)中考慮收益波動性的報童問題這一方向還有不少問題有待研究，如在需求函數(shù)方面，對于常見的乘法形式和指數(shù)形式等函數(shù)形式對問題的影響；在波動性方面，將測度波動性的常用工具變異系數(shù)融入報童模型，研究衡量波動性對決策的影響等，都是在這個方向上值得深入思考和進(jìn)一步研究的問題。

[1]Petruzzi N,Dada M.Pricing and the News Vendor Problem:a Review with Extensions[J].OperationsResearch,1999,47(2):183-194.

[2]Qin Y,Wang R,Vakharia A,et al.The newsvendor problem:Review and directions for future research[J].European Journal of Operational Research,2011,213(2):361-374.

[3]Xu M,Chen Y,Xu X.The effect of demand uncertainty in a price-setting newsvendor model[J].European Journal of Operational Research,2011,207(2):946-957.

[4]Xu M,Lu Y.The effect of supply uncertainty in price-setting newsvendor models[J].European Journal of Operational Research,2013,227(3):423-433.

[5]Dogan A.Production and pricing policies in dual sourcing supply chains[J].Transportation Research Part E Logistics&Transportation Review,2015,76:1-12.

[6]Jadidi O,Taghipour S,Zolfaghari S.A two-price policy for a newsvendor product supply chain with time and price sensitive demand[J].European Journal of Operational Research,2016,253(1):132-143.

[7]Khouja M.Forward-buying and the naive newsvendor[J].Operations Research Letters,2016,44(2):238-242.

[8]Gerchak Y,Mossman D.On the Effect of Demand Randomness on Inventories and Costs[J].Operations Research,1992,40(4):804-807.

[9]Chua G,Liu Y.On the Effect of Demand Randomness on Inventory,Price and Profit[J].Operations Research Letters,2015,43(5):514-518.

[10]丁小東,徐菱,蔣葛夫,等.考慮浪費規(guī)避偏差和缺貨規(guī)避偏差的零售商訂購行為[J].計算機(jī)集成制造系統(tǒng),2014,20(10):2 582-2 598.

[11]Javier R,Melike B,Anna J.A price-setting newsvendor problem under mean-variance criteria[J].European Journal of Operational Research,2015,247(2):575-587.

[12]杰克·韋爾奇,蘇茜·韋爾奇.贏[M].北京:中信出版社,2013.

[13]葉飛,林強(qiáng),李怡娜.基于CVaR的“公司＋農(nóng)戶”型訂單供應(yīng)鏈協(xié)調(diào)契約機(jī)制[J].系統(tǒng)工程理論與實踐,2011,31(3):450-460.

[14]朱傳波,季建華,包興.供應(yīng)風(fēng)險規(guī)避下基于VaR的零售商訂貨策略[J].系統(tǒng)管理學(xué)報,2014,23(6):861-866.

[15]程露,萬仲平,侯闊林,蔣威.CVaR準(zhǔn)則下的雙層報童問題模型研究[J].運籌學(xué)學(xué)報,2008,12(4):83-93.

[16]許明輝,于剛,張漢勤.帶有缺貨懲罰的報童模型中的CVaR研究[J].系統(tǒng)工程理論與實踐,2006,(10):1-8.

[17]馮艷剛,吳軍.突發(fā)事件環(huán)境下風(fēng)險厭惡型報童博弈模型[J].系統(tǒng)工程理論與實踐,2015,35(3):598-607.

[18]Kocabiyikoglu A,Popescu I.An Elasticity Approach to the Newsvendor with Price-Sensitive Demand[J].Operations Research,2011,59(2):301-312.

Study on Risk Aversive Newsboy Problem Considering Earnings Fluctuation

Tian Zhiyong,Zhou Li,Yuan Ruiping
(School of Informaiton,Beijing Wuzi University,Beijing 101149,China)

In this paper,targeting at the fluctuation of the newsboy earnings,we built the risk aversive newsboy model with singlephased additive stochastic demand function and price and inventory as the decision-making variables.We analyzed the model respectively through sequential optimization and simultaneous optimization,investigated the relationship between the optimal solution of the model and such parameters as risk coefficient and disposal unit price of the residual inventory,etc.,and used the concept of elasticity to describe the condition of existence and practical significance of the optimal solution of the model.Through the analysis,we found that the optimal solutions yielded by both methods were of the same nature,in that both the price and the inventory quantity were decreasing function to the risk coefficient and increasing function to the disposal unit price of the residual inventory.As for the condition of existence of the optimal solution,the local optimal solution through simultaneous optimization was the least stringent in prerequisite,followed by the optimal solution through sequential optimization,and further followed by the global optimal solution through simultaneous optimization.

newsboy model;earnings fluctuation;risk aversive;elasticity

F224.0

A

1005-152X(2017)07-0066-06

10.3969/j.issn.1005-152X.2017.07.015

2017-05-06

國家自然科學(xué)基金資助項目(71501015)；北京智能物流系統(tǒng)協(xié)同創(chuàng)新中心；智能物流系統(tǒng)北京市重點實驗室(BZ0211)

田志勇，男，北京物資學(xué)院信息學(xué)院副教授，博士，研究方向：低碳物流、電子商務(wù)。