李曉東

（沈陽鐵路局 沈陽客運(yùn)段，遼寧 沈陽 110013）

我國經(jīng)濟(jì)的快速發(fā)展、社會的進(jìn)步改變了人們的工作方式和生活習(xí)慣，同時也改變了對運(yùn)輸?shù)男枨?。旅客的?jīng)濟(jì)能力和時間價值觀念逐漸增強(qiáng)，更加關(guān)注運(yùn)輸質(zhì)量，包括出行的便捷性、旅行的舒適性、旅行時間的縮短、安全性的提高、運(yùn)輸服務(wù)水平等[1]。目前，鐵路運(yùn)輸面臨著運(yùn)輸數(shù)量和質(zhì)量的雙重壓力。為解決鐵路運(yùn)輸?shù)摹捌款i”問題，我國進(jìn)行了既有線提速和高速鐵路建設(shè)，預(yù)計 2012年投產(chǎn)新線 6 366 km。“十二五”期間，全國鐵路運(yùn)營里程將達(dá)到12萬 km，其中以高速鐵路為主骨架的快速鐵路網(wǎng)達(dá)4.5萬 km，西部地區(qū)鐵路達(dá) 5萬 km。這將使鐵路既有的優(yōu)勢得以充分發(fā)揮，同時也為鐵路的發(fā)展注入新的活力。

1 客運(yùn)量預(yù)測模型的選擇

鐵路客運(yùn)量既是衡量運(yùn)營效果的一項重要指標(biāo)，運(yùn)輸部門可以根據(jù)客運(yùn)量的變化及時調(diào)整戰(zhàn)略部署，同時也是鐵路項目可行性研究中的重點(diǎn)內(nèi)容，是論證項目建設(shè)必要性和經(jīng)濟(jì)性的基礎(chǔ)。國內(nèi)外采用的鐵路客流量預(yù)測方法較多，其中定量預(yù)測方法可分為因果關(guān)系分析法、產(chǎn)運(yùn)銷平衡法、四階段法等。這些方法在預(yù)測精度上都有待提高，以便為決策者提供更可靠的數(shù)據(jù)。在鐵路既有線提速前，客運(yùn)量預(yù)測只是針對趨勢客運(yùn)量。但是，既有線實(shí)施提速和高速鐵路建成投入運(yùn)營，鐵路客運(yùn)綜合服務(wù)質(zhì)量得到提高，如果照搬某一種方法，預(yù)測結(jié)果將會出現(xiàn)偏差，并且隨著時間的推移，理論預(yù)測值和未來實(shí)際運(yùn)量的誤差也會越來越大[2]。

回歸分析是一種應(yīng)用極為廣泛的數(shù)量分析方法，主要用于確定2種或2種以上變量間相互依賴的定量關(guān)系?；貧w分析的基本思想是：雖然自變量和因變量之間沒有嚴(yán)格的、確定性的函數(shù)關(guān)系，但可以設(shè)法找出最能代表它們之間關(guān)系的數(shù)學(xué)表達(dá)形式。因此，可以通過因變量和自變量來確定變量之間的因果關(guān)系，建立線性回歸模型，并根據(jù)實(shí)測數(shù)據(jù)求解模型的各個參數(shù)，然后評價回歸模型是否能夠很好地擬合實(shí)測數(shù)據(jù)，如果能夠很好地擬合，則可以根據(jù)自變量作進(jìn)一步預(yù)測[3-4]。

線性回歸預(yù)測模型對經(jīng)濟(jì)社會系統(tǒng)預(yù)測具有重要的意義，但由于客運(yùn)量的影響因素復(fù)雜，其變化呈現(xiàn)一定的隨機(jī)波動特點(diǎn)，因而預(yù)測結(jié)果會產(chǎn)生誤差。而馬爾可夫鏈?zhǔn)歉鶕?jù)系統(tǒng)狀態(tài)之間的轉(zhuǎn)移矩陣來描述一個隨機(jī)動態(tài)系統(tǒng)未來的發(fā)展?fàn)顟B(tài)，轉(zhuǎn)移矩陣反映了各個狀態(tài)之間的規(guī)律性。馬爾可夫鏈過程最主要的特征就是無后效性，即t時刻之后的狀態(tài)只與 t時刻的狀態(tài)有關(guān)，而與以前的狀態(tài)無關(guān)。因此，可以通過對線性回歸預(yù)測模型的結(jié)果進(jìn)行馬爾可夫鏈改進(jìn)來提高其預(yù)測的準(zhǔn)確性[5]。

2 鐵路客運(yùn)量預(yù)測

2.1 線性回歸模型確定預(yù)測值

鐵路客運(yùn)量的影響因素包括國民生產(chǎn)總值、人口數(shù)量、居民消費(fèi)水平、運(yùn)價、運(yùn)輸服務(wù)水平、運(yùn)輸方式結(jié)構(gòu)等。目前，國家制定了“十二五”規(guī)劃，但以上因素仍存在一定的不確定性。因此，選取鐵路營業(yè)里程這一相對確定的指標(biāo)進(jìn)行線性回歸分析。2001—2010年我國鐵路客運(yùn)量和鐵路營業(yè)里程如表1所示[6]。根據(jù)數(shù)據(jù)作散點(diǎn)圖，其趨勢大致呈線性關(guān)系，如圖1所示。

表1 2001—2010年全國鐵路客運(yùn)量和鐵路營業(yè)里程

圖1 2001—2010年全國鐵路客運(yùn)量和鐵路營業(yè)里程關(guān)系散點(diǎn)圖

按數(shù)據(jù)趨勢可以建立線性方程，應(yīng)用 MATLAB 軟件，回歸方程為：

式中：X 為營業(yè)里程，萬km；Y 為鐵路客運(yùn)量，萬人。

回歸模型擬合度檢驗(yàn)結(jié)果顯示，相關(guān)系數(shù)為0.950，而決定系數(shù) R2為0.903，回歸方程對樣本數(shù)據(jù)點(diǎn)擬合較好。

根據(jù)公布的數(shù)據(jù)：2011年鐵路新增線路 2 022 km，2012年計劃新增線路 6 366 km，2015年計劃鐵路線路全長達(dá)到 12萬 km。根據(jù)上述公式可以得到鐵路客運(yùn)量的回歸預(yù)測值及相對誤差，如表2所示。

表2 鐵路客運(yùn)量回歸預(yù)測相關(guān)數(shù)據(jù)

2.2 回歸-馬爾可夫預(yù)測模型確定預(yù)測值

馬爾可夫鏈?zhǔn)歉鶕?jù)系統(tǒng)狀態(tài)之間的轉(zhuǎn)移矩陣來描述一個隨機(jī)動態(tài)系統(tǒng)未來的發(fā)展?fàn)顩r，轉(zhuǎn)移概率(即在事件的發(fā)展變化過程中，從某一種狀態(tài)出發(fā)下一時刻轉(zhuǎn)移到其他狀態(tài)的可能性) 則反映了各狀態(tài)之間某種內(nèi)在的規(guī)律性。其中，n 步狀態(tài)轉(zhuǎn)移矩陣的計算為：

P(n) 描述了 m個狀態(tài)相互轉(zhuǎn)移的概率分布。其中，pij為由狀態(tài) i 轉(zhuǎn)移到狀態(tài)j 的概率[7-8]。如 p12表示由狀態(tài)1轉(zhuǎn)移到狀態(tài)2的概率。下一時刻向各狀態(tài)轉(zhuǎn)移的概率 P(E) 為：

其中 E1×z是一個1行 z 列的矩陣，為當(dāng)前時刻的狀態(tài)矩陣，z 為劃分的狀態(tài)區(qū)間。

回歸—馬爾可夫預(yù)測模型的基本思想是先用回歸預(yù)測模型進(jìn)行預(yù)測。然后根據(jù)回歸預(yù)測數(shù)據(jù)與原始數(shù)據(jù)的相對誤差將各原始數(shù)據(jù)年劃分狀態(tài)。再根據(jù)各年所處狀態(tài)求得狀態(tài)轉(zhuǎn)移矩陣，進(jìn)行未來相關(guān)數(shù)據(jù)的預(yù)測。

2.2.1 狀態(tài)劃分

根據(jù)表2中的相對誤差列，將原始數(shù)據(jù)按年劃分為4個狀態(tài)，劃分標(biāo)準(zhǔn)如表3所示。例如，2009年的回歸預(yù)測相對誤差為 0.56%，則該年屬于狀態(tài) 3。

表3 狀態(tài)劃分區(qū)間及各狀態(tài)含義

2.2.2 狀態(tài)轉(zhuǎn)移

鐵路客運(yùn)量在 2001—2010年中各類狀態(tài)的一步轉(zhuǎn)移情況如表4所示。例如，表4中第二行的數(shù)據(jù)含義為：由狀態(tài) 1一步轉(zhuǎn)移到狀態(tài)1和一步轉(zhuǎn)移到狀態(tài)3的情況各出現(xiàn)了1次，合計為2次。表4中的其他數(shù)據(jù)依次類推。

表4 回歸預(yù)測結(jié)果馬爾可夫狀態(tài)轉(zhuǎn)移

2.2.3 數(shù)據(jù)預(yù)測

根據(jù) 2001—2010年鐵路客運(yùn)量數(shù)據(jù)，以2010年為當(dāng)前時間點(diǎn)預(yù)測 2011年鐵路客運(yùn)量，需要計算一步狀態(tài)轉(zhuǎn)移矩陣；預(yù)測 2012 與 2015年的鐵路客運(yùn)量，需要分別計算二步狀態(tài)轉(zhuǎn)移矩陣與五步狀態(tài)轉(zhuǎn)移矩陣。下面以運(yùn)用回歸—馬爾可夫模型預(yù)測2011年鐵路客運(yùn)量為例，說明該模型的具體計算過程。

從狀態(tài)的轉(zhuǎn)移表可確定馬爾可夫一步狀態(tài)轉(zhuǎn)移矩陣P (1)。

由于 2010年鐵路客運(yùn)量所處的狀態(tài)為狀態(tài) 3，對應(yīng) 2010年的狀態(tài)矩陣為 (0，0，1，0)，因而 2011年鐵路客運(yùn)量轉(zhuǎn)向各狀態(tài)的概率為：

即 2011年鐵路客運(yùn)量有 0.25 的概率處于狀態(tài)2和狀態(tài) 4，有 0.5 的概率處于狀態(tài) 3。

由各狀態(tài)區(qū)間的預(yù)測中值，分別乘以預(yù)測年可能處于該狀態(tài)的概率再求和，即：175 238×0.25＋184 225×0.50＋196 806×0.25＝185 124(萬人)，得到 2011年鐵路客運(yùn)量回歸—馬爾可夫預(yù)測值為185 124萬人。同理，計算得到 2012年和 2015年鐵路客運(yùn)量的預(yù)測結(jié)果，如表5所示。

表5 鐵路客運(yùn)量回歸-馬爾可夫預(yù)測值

3 結(jié)束語

基于回歸—馬爾可夫預(yù)測模型的 2011年、2012年鐵路客運(yùn)量預(yù)測值分別為 185 124萬人和 199 427萬人，與 2011年實(shí)際數(shù)據(jù)和 2012年預(yù)計數(shù)據(jù)的誤差為-0.04% 和 -1.72%。因此，在鐵路客運(yùn)量預(yù)測方面，對線性回歸預(yù)測模型結(jié)果進(jìn)行馬爾可夫鏈改進(jìn)的方法能夠提高預(yù)測的準(zhǔn)確性。但還需通過更多的實(shí)際數(shù)據(jù)對預(yù)測模型的可靠性進(jìn)行驗(yàn)證，以不斷完善回歸—馬爾可夫預(yù)測模型。

[1]郭孜政. 鐵路客運(yùn)量影響因素層次結(jié)構(gòu)分析[J]. 交通運(yùn)輸工程與信息學(xué)報，2007(4)：68-71.

[2]王傳峰. 既有線提速后的客運(yùn)量預(yù)測研究[D]. 四川：西南財經(jīng)大學(xué)，2008.

[3]侯麗敏，馬國峰. 基于灰色線性回歸組合模型鐵路客運(yùn)量預(yù)測[J]. 計算機(jī)仿真，2011(7)：1-3，30.

[4]謝孝如，蔣惠園，申耀偉. 基于并聯(lián)灰色—線性回歸組合模型的客運(yùn)量預(yù)測[J]. 鐵道運(yùn)輸與經(jīng)濟(jì)，2008(8)：92-94.

[5]田自力，劉碧發(fā). 灰色馬爾可夫鏈預(yù)測模型[C]. 臺灣第一屆灰色系統(tǒng)理論與應(yīng)用研討，1996：157-161.

[6]國家統(tǒng)計局. 中國統(tǒng)計年鑒：1999—2011 [M]. 北京：中國統(tǒng)計出版社，1999—2011.

[7]關(guān) 靜. 中國民航旅客吞吐量的灰色馬爾可夫模型預(yù)測[J].中國民航大學(xué)學(xué)報，2010(2)：45-47.

[8]張 誠，張廣勝. 基于優(yōu)化灰色—馬爾可夫鏈模型的鐵路貨運(yùn)量預(yù)測[J]. 物流技術(shù)，2011(13)：129-142.