AR(P)模型的P取值試驗

曾山

(重慶育才工程咨詢監理有限公司，重慶 400074)

AR(P)模型的P取值試驗

曾山?

(重慶育才工程咨詢監理有限公司，重慶 400074)

AR(P)模型的p值求解，嚴格而言應通過計算解決，但無論哪一種計算方法，都要求一定的，甚至足夠大的樣本量，而土木工程施工監控的相關規范規定，監控落實的最晚期限，也不能超過3 d～5 d，顯然僅依據3個～5個(組)樣本，算不出正確的p值來，因此在具體工作中，p值往往人為設置。不同p值的AR(P)模型各自具有不同的特點。筆者分別應用p=2、p=3的AR(P)模型，進行隧道邊墻收斂變形分析，認為，若需要盡早落實監控工作，預測精度要求又不是太高，模型的p宜取2；如果項目工期較長，預報精度期望較高的情況下，p取3更合適。

AR(P)模型；預測預報；精度

1 問題的提出

影響AR(P)模型的適應性的因素，包括模型辨識、參數估計、模型定階(p的取值)。求p值，可用偏相關函數法、F準則、FPE準則、AIC準則等等，但無論哪一種方法，都要求一定的，甚至足夠大的樣本量，而土木工程的施工監控，卻是越早落實越好，最長期限也不超過3 d～5 d[1]，樣本貧乏，導致p值無法計算，而只能人為設定。筆者針對不同p值的AR(P)模型，各自具有的特征，采用施工現場數據，作了試驗。

2 AR(P)模型簡介[2～4]

AR(P)模型是時間序列分析模型的一種，其基本思想是，平穩、正態、零均值的時間序列值{xt}，不僅與前n步的各xt-1，xt-2，…xt-n值有關，而且與前m步的干擾at-1，at-2，…at-m有關(m，n=1，2…)，據多元線性回歸的思路，構成一般ARMA模型。

式中，φi(i=1，2…n)為回歸參數，{at}為白噪聲序列。若θj=0，可得AR(P)模型一般式:

式中＾xi為估值(預報值)、△i為觀測誤差。

相應的誤差方程組為:

φ的最小二乘估計:

最后得AR(P)模型的預報方程與殘差方程:

3 試驗分析

水麻公路SM隧道邊墻的相對收斂值實測數據[5]，列于表1(其中t為觀測時間，x為以天為周期的收斂數據)；依次取1 d～4 d的、1 d～5 d的、1 d～6 d的……直至1 d～14 d的收斂數據，分別代入式(8)，一一求出:＾x2(p=2時的收斂預報值)、v2(p=2時的預報殘差值)、＾x3(p=3時預報值)、v3(p=3預報殘差值)，并先后填入表1的相應欄；圖1為根據表1繪制的分析曲線(橫標為時間，縱標為收斂值)，圖中:1線根據x繪出，2線根據＾x2繪出，3線根據＾x3繪出。計算擬合度，知p=2時σ2=2.69， p=3時σ3=1.09，都未超3 mm的擬合度容許誤差，可認為兩個系列的回歸模型均滿足要求的精度，但后者的擬合精度相對更高[6]。

圖1 分析圖

分析表 表1

[1]中，沒有預報殘差的限定，規定隧墻日安全位移限差為3 mm～5 mm，筆者據此取5 mm為預

通過報殘差容許值。表1顯示，無論p取2或3，絕對殘差值均小于5 mm；p=2時，殘差取值范圍[-4.57，0.03]，均值為-2.46，p=3時，殘差取值范圍[-1.83，2.43]，均值為-0.36，可認為p取2和取3，都滿足限定的預報精度，但仍然是p取3的預報精度高。

如圖1所示，1線與2線，前4天高度吻合，但從第9天始，漸行漸遠；1線與3線，前兩天的距離稍大，也是從第9天始，后續點的吻合度明顯提高，并且波動幅度較小。

4 結 語

p取2時，只需4組原始數據即可建模，預測預報第4天即可進行，且前期的預報精度較高，模型整體擬合度尚可，但隨著時間推移，預報精度顯著下降，穩定度也不太好。因此，需要較早落實監控工作時，模型的p宜取2。

p取3時，建模需6組觀測數據，預測預報第7天才可進行，前期“空白”時間較長，但模型的整體擬合度高，預報精度頭兩天低一些，隨數據量的增加，迅速、明顯提升，趨向亦穩定。如果項目工期較長，現場條件也不復雜，預報精度期望高的情況下，p取3的模型更合用。

當然，兩種模型分別先后采用，更能發揮其優長。

參考文獻

[1] JTG F60-2009.公路隧道施工技術規范[S].

[2] 白鳳山，么煥民，李春玲等.數學建模[M].哈爾濱：哈爾濱工業大學出版社，2003.

[3] 吳懷宇.時間序列分析與綜合[M].武漢：武漢大學出版社，2004.

[4] 梅紅，岳東杰.時間序列分析在變形監測數據處理中的應用[J].現代測繪，2005(6).

[5] 徐兮.隧道邊墻水平收斂變形監測及分析研究[J].測繪通報，2007(7)：58～60.

[6] 王成.隧道工程[M].北京：人民交通出版社，2009.

Experimentation that Find the Value of P in AR(P)M odel

Zeng Shan
(Chongqing Yucai Project Supervision Limited Company，Chongqing 400074，China)

The usualmethod to find value of p in the AR(P)model ismathematical calculation，but all calculation methods need enough(even huge)amount of examples.By another side，related criterions for civil engineering project monitoring prescribe that the limited period of projectmonitoring is not over 3 or 5 days.So only by examples of 3-5 days，people can’t calculate the right value of p.In actualworking，the p adopts artificial value.What are characters in different AR(P)models with different values of p?The author adopts p=2、p=3 into AR(P)models to analyzing the tunnel side well convergence distortion.The author considers that if the project needmonitoring starting as early as possible and low forecasting precision requirement，the value of p in model equal to 2 is better；if the project has long construction period and high forecasting precision requirement，the p’s value equal to 3 is better.

AR(P)model；Forecasting；Precision

1672-8262(2013)02-158-02

TU196+.1

B

2012—09—06

曾山(1985—)，男，工程師，主要從事土木工程建設工作。