基于Verhulst與線性回歸組合模型的港口吞吐量預(yù)測研究

代雯強(qiáng)　楊　珩

(中交第二航務(wù)工程勘察設(shè)計院有限公司　武漢　430060)

基于Verhulst與線性回歸組合模型的港口吞吐量預(yù)測研究

代雯強(qiáng)楊珩

(中交第二航務(wù)工程勘察設(shè)計院有限公司武漢430060)

摘要由于港口吞吐量預(yù)測的復(fù)雜性，在許多情況下，單純利用一種特定的預(yù)測方法進(jìn)行預(yù)測往往具有片面性，為盡可能使預(yù)測結(jié)果具有較高可信度，基于“誤差平方和最小”為最優(yōu)準(zhǔn)則，建立灰色Verhulst時序殘差修正模型與一元線性回歸模型的組合模型，將其用于某港口貨物吞吐量預(yù)測，經(jīng)驗證組合模型具有較好的預(yù)測效果，并對該港口未來3年的貨物吞吐量進(jìn)行了預(yù)測。

關(guān)鍵詞吞吐量預(yù)測時序殘差修正Verhulst模型線性回歸模型組合模型

目前，港口吞吐量的預(yù)測方法有很多，定量預(yù)測法主要分為2大類：時間序列法和因果關(guān)系法。前一類預(yù)測法包括指數(shù)平滑法、自回歸法、移動平均法、灰色預(yù)測法等；后一類預(yù)測法包括回歸分析法、系統(tǒng)動力學(xué)法、投入產(chǎn)出法等[1]。

考慮到在中長期港口吞吐量預(yù)測中基于時序殘差修正的Verhulst模型可以很好地解決吞吐量增長到比較飽和階段的預(yù)測問題[2]，以某港為例，基于時序殘差修正Verhulst模型和一元線性回歸模型，以“誤差平方和最小”為準(zhǔn)則建立組合預(yù)測模型，并通過對模型精度進(jìn)行分析，選擇最優(yōu)模型，對港口吞吐量做出預(yù)測。

1數(shù)學(xué)模型

1.1　時序殘差修正Verhulst模型

灰色預(yù)測最大的優(yōu)點是以微分方程形式描述模型，能夠反映出事物發(fā)展的連續(xù)性，較適用于中長期的預(yù)測規(guī)劃，而Verhulst模型可描述未來年具有飽和狀態(tài)的預(yù)測[2]。這些均符合港口吞吐量預(yù)測的特點，因此，選用灰色Verhulst模型對吞吐量進(jìn)行預(yù)測是合適的。

1.1.1灰色X(0)-Verhulst 模型

X(0)表示n年的吞吐量數(shù)據(jù)序列：X(0)= (x(0)(1)，x(0)(1)，…，x(0)(n))

對X(0)做一次累加生成(即1-AGO生成)得到X(1)= (x(1)(1)，x(1)(1)，…，x(1)(n))。

設(shè)非負(fù)原始數(shù)據(jù)序列，其中

(1)

取x(1)(1)=x(0)(1)，則灰色Verhulst模型的時間響應(yīng)函數(shù)為

(2)

利用一次累減還原計算，得到原始數(shù)據(jù)序列擬合值：

(3)

式(2)，(3)即為x(0)-Verhulst模型，x(0)為歷年吞吐量數(shù)據(jù)序列，依據(jù)式(2)，(3)和x(0)可確定參數(shù)a，b，進(jìn)而得到港口吞吐量預(yù)測值。

1.1.2Verhulst 模型的時序殘差修正

利用上述x(0)-Verhulst模型計算并預(yù)測吞吐量存在一定誤差，且時序k按等時序處理。為提高模型精度，可對時序k進(jìn)行殘差修正。假設(shè)k存在殘差ek，令

(4)

(5)

建立時序殘差序列ek(0)的GM(1,1)模型，即

(6)

(7)

由式(2)，(3)可對吞吐量進(jìn)行初步預(yù)測，再由式(6)，(7)對時序k進(jìn)行殘差估計，進(jìn)而對吞吐量預(yù)測值進(jìn)行修正，最終得到未來年的吞吐量預(yù)測值。

1.2　線性回歸模型

回歸分析是通過研究與預(yù)測目標(biāo)相關(guān)的影響因素，并通過因素統(tǒng)計資料運(yùn)用最小二乘法擬合[3]。采用一元線性回歸模型預(yù)測港口吞吐量，考慮港口吞吐量與國內(nèi)生產(chǎn)總值(GDP)往往具有較好的相關(guān)性，選擇GDP作為相關(guān)關(guān)系分析中的自變量，以港口吞吐量作為因變量建立相應(yīng)的線性回歸模型如下：

X(t)=aY(t)+b

式中：t為年份；X(t)為港口吞吐量；Y(t)為GDP。

1.3　“誤差平方和最小”的組合預(yù)測模型

(1) 組合預(yù)測模型建立。設(shè)某一預(yù)測問題在某一時段的實際值為Y(t)，基于上述2種單向預(yù)測，分別設(shè)2種單向預(yù)測法預(yù)測吞吐量值分別為V(0)(t)，X(0)(t)(t=1,2，…，n)，根據(jù)參考文獻(xiàn)[4-5]，可建立如下方程組：

(8)

求解該方程，可以得到w1，w2，則“誤差平方和最小”的組合預(yù)測模型目標(biāo)函數(shù)為：

(9)

(2) 模型精度。采用殘差檢驗，計算其模型精度η。

2某港吞吐量預(yù)測

以某港為例，選取某港2008～2013年港口貨物吞吐量數(shù)值及國內(nèi)生產(chǎn)總值GDP數(shù)值作為建模數(shù)據(jù)(見表1)，先對模型進(jìn)行驗證，包括時序殘差修正Verhulst模型、一元線性回歸模型以及基于“誤差平方和最小”的組合預(yù)測模型。驗證完畢后，對3個模型進(jìn)行比較和分析。

表1　某港2008-2013年港口貨物吞吐量及GDP統(tǒng)計表

2.1　時序殘差修正Verhulst模型

根據(jù)1.1計算得，時序殘差修正Verhulst模型中a=-0.139 250，b=4 549.889，預(yù)測結(jié)果見表2。

經(jīng)模型精度檢驗，模型精度η=0.972>0.9，殘差檢驗結(jié)果為很好。另外，采用后差檢驗法，經(jīng)計算c=0.000 012 4<0.35，p=1，可知模型精度1級(好)，可用所建的模型進(jìn)行吞吐量預(yù)測。

2.2　一元線性回歸模型

考慮港口吞吐量與國內(nèi)生產(chǎn)總值(GDP)往往具有較好的相關(guān)性，選擇GDP作為相關(guān)關(guān)系分析中的自變量，以港口吞吐量作為因變量建立相應(yīng)的線性回歸模型如下：

X(t)=0.975Y(t)+1933.603

(10)

式中:t為年份；X(t)為港口吞吐量；Y(t)為GDP。

經(jīng)相關(guān)系數(shù)檢驗，該回歸模型顯著性較好。預(yù)測結(jié)果見表3。

表3　一元線性回歸模型預(yù)測值(萬t)及精度

2.3　“誤差平方和最小”的組合預(yù)測模型

根據(jù)上述2種模型的預(yù)測結(jié)果，參考公式(10)建立方程組如下：

求解該方程，可以得到w1=0.557 595 101，w2=0.442 404 899。經(jīng)計算，時序殘差修正Verhulst模型預(yù)測的權(quán)系數(shù)為0.557 595 101，一元線性回歸模型預(yù)測的權(quán)系數(shù)為0.442 404 899。將其代入式(9)，得到組合預(yù)測模型：

0.442 404 899X(0)(t)t=1,2,…,n

分別將t=1, 2,…,n代入即可得到組合模型預(yù)測的吞吐量，見表4。

表4　組合模型預(yù)測值(萬t)及精度

采用殘差檢驗，計算其模型精度η=0.973>0.9，殘差檢驗結(jié)果為很好，組合模型可用于吞吐量預(yù)測。

2.4　模型評價

根據(jù)上述結(jié)果，從單一模型來看，時序殘差修正Verhulst模型的預(yù)測效果好于線性回歸模型，表明其更能反映出該港目前的吞吐量發(fā)展趨勢。

組合預(yù)測模型的均方百分比誤差最小，預(yù)測的精度較單一的預(yù)測模型有所提高，其應(yīng)作為該港吞吐量預(yù)測的最優(yōu)模型。

2.5　運(yùn)用組合模型預(yù)測該港港口吞吐量

最后運(yùn)用3種模型預(yù)測2018~2020年共3年某港港口貨物吞吐量。預(yù)測結(jié)果見表5。

表5　某港2018~2020年港口貨物吞吐量預(yù)測值

3結(jié)論

(1) 時序殘差修正Verhulst模型可較好地描述未來年具有飽和狀態(tài)的港口吞吐量預(yù)測，其預(yù)測精度相較于線性回歸模型有所提高，表明其更能反映出該港目前的吞吐量發(fā)展趨勢。

(2) 基于時序殘差修正Verhulst模型和一元線性回歸模型，以“誤差平方和最小”為準(zhǔn)則建立組合預(yù)測模型新建的組合模型，經(jīng)驗證組合模型具有較好的預(yù)測效果。

(3) 經(jīng)比較，建立的組合預(yù)測模型均方百分比誤差最小，預(yù)測的精度較單一的預(yù)測模型有所提高，其應(yīng)作為該港吞吐量預(yù)測的最優(yōu)模型。

(4) 通過對3種模型精度進(jìn)行分析，選擇最優(yōu)模型，嘗試對港口吞吐量做出預(yù)測，由預(yù)測值可以看出，未來某港的貨物吞吐量仍將呈上升趨勢，2018相比2013年遞增幅度達(dá)到15%，2019年、2020年增幅為14%。

參考文獻(xiàn)

[1]翟希東.港口集裝箱吞吐量預(yù)測模型研究[D].大連：大連理工大學(xué)，2006.

[2]葉劍，宋向群，郭子堅.港口吞吐量預(yù)測的時序殘差修正Verhulst模型[J].水運(yùn)工程，2004(12)：18-22.

[3]劉明維，王多銀.港口貨物吞吐量預(yù)測方法探討[J].水運(yùn)工程，2005(3)：53-55.

[4]朱超.組合預(yù)測在港口吞吐量預(yù)測中的應(yīng)用研究[J].港工技術(shù)，2006(9)：10-12.

[5]張浩.基于最優(yōu)線性組合的港口集裝箱吞吐量預(yù)測方法[J].武漢理工大學(xué)學(xué)報:交通科學(xué)與工程版,2007(2):373-376.

Port Throughput Forecasting Method Based on Time Series Error

Corrected Verhulst Model and Linear Regression Combined Model

DaiWenqiang,YangHeng

(CCCC Second Harbor Consultants Co., Ltd., Wuhan 430063, China)

Abstract：Owing to the complexity of port throughput forecast, using one prediction technique always is partial in many cases. In order to increase the reliability of result, aiming the square sum of error, the combined model of gray Verhulst series error correction model and simple linear regression model is constructed to forecast the port throughput capacity. It is proved that combination model's prediction is effective and the port throughput capacity in the next 3 years is predicted.

Key words:throughput forecasting; verhulst model on time series error correction; linear regression; combined model

收稿日期：2015-04-15

DOI 10.3963/j.issn.1671-7570.2015.03.059