□文/肖亞明

（安徽工業(yè)大學(xué)研究生學(xué)院 安徽·馬鞍山）

一、模型構(gòu)建

為了便于研究，首先構(gòu)建一個(gè)簡單的二級供應(yīng)鏈系統(tǒng)，我們做如下假設(shè)：

*供應(yīng)鏈上有2 個(gè)企業(yè)，它們分別是1 個(gè)零售商和1 個(gè)制造商；

*零售商直接滿足顧客需求，沒有運(yùn)輸延遲時(shí)間；

*零售商不準(zhǔn)缺貨，制造商的生產(chǎn)能力無限；

*訂單延遲以及運(yùn)輸延遲時(shí)間均為1 周，制造商生產(chǎn)周期為2 周；

*各節(jié)點(diǎn)企業(yè)安全庫存以及初始庫存都為100，單位件；

*各節(jié)點(diǎn)企業(yè)均采用安全庫存策略；

*模擬周期為50 周；

*各個(gè)周期的訂貨量以及發(fā)貨量大于或等于零；

*單位庫存以及訂貨費(fèi)用均為1 元，單位缺貨費(fèi)用為2 元。

相關(guān)說明：設(shè)定第三個(gè)假設(shè)的意義在于保證零售商與制造商的總時(shí)間延遲相同，目的是為了在時(shí)間延遲相同的條件下分析市場需求在向上傳遞過程中的放大現(xiàn)象。零售商的訂單延時(shí)與制造商運(yùn)輸延遲時(shí)間之和為2 周，這就意味著零售商從發(fā)出訂單到最后收到貨物實(shí)際上經(jīng)過了2 周的延遲時(shí)間。制造商的生產(chǎn)周期也為2 周，表示制造商從制定生產(chǎn)計(jì)劃到貨物進(jìn)入倉庫成為庫存的時(shí)間也為2 周。這就剔除了因?yàn)檠舆t時(shí)間不同對需求放大所造成的影響。

該模型中相關(guān)變量以及公式定義如下：

*市場需求：市場需求是用MATLAB 隨機(jī)函數(shù)產(chǎn)生的一組均勻分布隨機(jī)變量。具體公式如下：NORMINV（RAND（）），80，10）

*市場銷售=MIN（上周庫存+本周到貨，市場需求）

*訂貨量=MAX（（安全庫存-庫存）/庫存調(diào)整時(shí)間+本周需求，0），初始值為0

*本周需求=下游企業(yè)訂貨量（延遲時(shí)間以前）

*庫存量=上周庫存量+本周到貨-本周發(fā)貨量，各節(jié)點(diǎn)企業(yè)庫存初始值都為100

*本周到貨=上游發(fā)貨量（運(yùn)輸延遲時(shí)間以前）

*發(fā)貨量=MAX（MIN（上周庫存+本周到貨，本周需求+上周缺貨），0）

*缺貨量=本周需求-本周發(fā)貨量

*生產(chǎn)需求=（安全庫存-庫存）/庫存調(diào)整時(shí)間+本周需求

說明：上述公式中“本周到貨”是由上游供應(yīng)商的“發(fā)貨量”經(jīng)過運(yùn)輸延遲時(shí)間后傳入的，它們在數(shù)值上是相等的?！氨局苄枨蟆笔怯上掠涡枨蠓降摹坝嗀洝苯?jīng)過訂貨延遲傳入的，它們在數(shù)值上也是相等的。

二、模型仿真

首先針對上述模型，運(yùn)用matlab 進(jìn)行仿真，得出牛鞭效應(yīng)圖，如圖1所示。（圖1）

圖1說明：黑線為市場需求；紅線為零售商需求；藍(lán)線為制造商生產(chǎn)需求。

當(dāng)需求出現(xiàn)波動時(shí)，運(yùn)用matlab 仿真，得出如下牛鞭效應(yīng)圖，如圖2所示。（圖2）

圖2說明：黑線為市場需求；紅線為零售商需求；藍(lán)線為制造商生產(chǎn)需求。

比較圖1和圖2可以看出，當(dāng)需求出現(xiàn)波動時(shí)，制造商的生產(chǎn)需求峰值達(dá)到了170，牛鞭效應(yīng)加劇。

三、在市場預(yù)測環(huán)節(jié)加入移動平均法，分析其牛鞭效應(yīng)情況

移動平均法是取預(yù)測對象最近一組歷史數(shù)據(jù)的平均值作為預(yù)測值的方法。這種方法不是僅取最近一期的歷史數(shù)據(jù)作為下一期的預(yù)測值，而是取最近一組歷史數(shù)據(jù)的平均值作為下一期的預(yù)測值，這一方法使近期歷史數(shù)據(jù)參與預(yù)測，使歷史數(shù)據(jù)的隨機(jī)成分有可能互相抵消，平均值所含的隨機(jī)成分就會相應(yīng)減少。

移動平均法的“平均”是指對歷史數(shù)據(jù)的“算術(shù)平均”，而“移動”是指參與平均的歷史數(shù)據(jù)隨預(yù)測值的推進(jìn)而不斷更新。當(dāng)一個(gè)新的歷史數(shù)據(jù)進(jìn)入平均值時(shí)，要剔除原先參與預(yù)測平均的最陳舊的一個(gè)歷史數(shù)據(jù)，并且每一次參與平均的歷史數(shù)據(jù)的個(gè)數(shù)是相同的。

其計(jì)算公式為：Ft+1=

其中，F(xiàn) 表示預(yù)測值，X 表示歷史數(shù)據(jù)，n表示參與移動平均的數(shù)據(jù)個(gè)數(shù)。

為了便于研究對比，在此取前面的研究數(shù)據(jù)運(yùn)用移動平均法進(jìn)行預(yù)測，仿真結(jié)果如下：

運(yùn)用matlab 對幾個(gè)參數(shù)進(jìn)行仿真，得出牛鞭效應(yīng)曲線，如圖3所示。（圖3）

圖3說明：黑線為市場需求；紅線為零售商需求；藍(lán)線為制造商生產(chǎn)需求。

根據(jù)仿真結(jié)果我們可以看到使用移動平均法后，訂貨以及庫存波動的趨勢都有所減弱，對照圖1與圖3我們可看到制造商生產(chǎn)需求的峰值從140 減少到115，振動明顯減弱，有效地減小了系統(tǒng)的牛鞭效應(yīng)。據(jù)此可以看出，在需求預(yù)測環(huán)節(jié)，采用移動平均法能夠減小牛鞭效應(yīng)，有效地增加了供應(yīng)鏈系統(tǒng)的穩(wěn)定性。

當(dāng)需求出現(xiàn)波動時(shí)，使用移動平均法的有效性研究。針對需求突變的情形，現(xiàn)使用移動平均法對輸入變量進(jìn)行仿真，仿真結(jié)果如圖4所示。（圖4）

圖4說明：黑線零售商庫存；紅線為制造商庫存；藍(lán)線為制造商生產(chǎn)需求。

比較圖2及圖4我們可以看出，在需求出現(xiàn)突變也即波動時(shí)，圖2振動加劇也即系統(tǒng)的不穩(wěn)定性明顯加劇，牛鞭效應(yīng)增大，制造商的最大需求量達(dá)到了172。在使用移動平均法對需求進(jìn)行預(yù)測后，觀察圖4，系統(tǒng)振幅減弱，需求曲線振動減小，制造商的最大需求也減小到了150 以下。據(jù)此可以看出，在需求預(yù)測環(huán)節(jié)，即使出現(xiàn)需求波動，采用移動平均法也能夠減小牛鞭效應(yīng)，有效地增加了供應(yīng)鏈系統(tǒng)的穩(wěn)定性?！?/p>