任桃紅+趙聯(lián)文

【摘要】居民消費(fèi)價(jià)格指數(shù)是我國(guó)物價(jià)指數(shù)體系中極其重要的一個(gè)指數(shù)，主要反映消費(fèi)者支付商品和勞務(wù)的價(jià)格變化情況，也是一種度量通貨膨脹水平的工具，以百分比變化為表達(dá)形式。我國(guó)改革開(kāi)放以來(lái)，社會(huì)經(jīng)濟(jì)的各方面發(fā)生巨大的變換，居民消費(fèi)價(jià)格指數(shù)的變動(dòng)也顯示出自身的特點(diǎn)，本文主要是應(yīng)用非平穩(wěn)時(shí)間序列和數(shù)學(xué)統(tǒng)計(jì)軟件Eviews對(duì)歷年的居民消費(fèi)價(jià)格指數(shù)進(jìn)行相關(guān)分析、處理并建立模型，這有利于我們認(rèn)識(shí)它與社會(huì)經(jīng)濟(jì)發(fā)展相聯(lián)系的變動(dòng)規(guī)律，文章最后根據(jù)所建立的模型舉例進(jìn)行了預(yù)測(cè)。

【關(guān)鍵詞】時(shí)間序列 非平穩(wěn)序列 乘積季節(jié)模型 建模

一、相關(guān)知識(shí)

時(shí)間序列分析是一種動(dòng)態(tài)數(shù)據(jù)處理的統(tǒng)計(jì)方法。一個(gè)時(shí)間序列{xt}是長(zhǎng)期趨勢(shì)變動(dòng)Tt，季節(jié)效應(yīng)St，周期變動(dòng)Ct和不規(guī)則變動(dòng)因素It共同作用的結(jié)果。確定性時(shí)間序列分析方法的原理是用確定性函數(shù)對(duì)長(zhǎng)期趨勢(shì)變動(dòng)Tt，季節(jié)效應(yīng)St，周期變動(dòng)Ct進(jìn)行擬合，將非平穩(wěn)時(shí)間序列平穩(wěn)化。本文我們?yōu)榱似椒€(wěn)化引入乘積季節(jié)模型。

乘積季節(jié)模型為：

ψ（B）U（BS）（1-B）d（1-BS）DXt=θ（B）V（BS）εt

E（εt）=0，Var（εt）=σ2ε，E（εSεK）=0，（s≠k）

EXkεt=0，（？坌k

其中：

U（BS）=1-u1BS-u2B2s-…-uPBPS

V（BS）=1-v1BS-v2B2s-…-uQBQS

可以對(duì)不同周期的同一周期點(diǎn)之間的相關(guān)性進(jìn)行擬合：

ψ（B）=1-ψ1B-ψ2B2-…-ψPBP

θ（B）=1-θ1B-θ2B2-…-θqBq

用來(lái)消除同一周期不同周期點(diǎn)之間的相關(guān)性。

季節(jié)模型簡(jiǎn)記為ARIMA（p，d，q）×（P，D，Q）S其中p和q是消除同一周期不同周期點(diǎn)之間相關(guān)性的自回歸階數(shù)和移動(dòng)平均階數(shù)，P和Q是消除不同周期的同一周期點(diǎn)之間的相關(guān)性的自回歸階數(shù)和移動(dòng)平均階數(shù)，s為周期長(zhǎng)度，d為差分的階數(shù)，D為季節(jié)差分的階數(shù)。

二、數(shù)據(jù)預(yù)處理

本文選取的是2000年1月到2011年12月全國(guó)居民消費(fèi)價(jià)格指數(shù)的數(shù)據(jù)，記居民消費(fèi)價(jià)格指數(shù)為具體數(shù)據(jù)如下表1：

通過(guò)Eviews軟件可得到時(shí)間序列y（t）的樣本自相關(guān)系數(shù)圖形可以發(fā)現(xiàn)，自相關(guān)系數(shù)呈線性衰減，衰減速度緩慢，這說(shuō)明該序列y（t）含有一定的趨勢(shì)性。同時(shí)可以看出該序列不僅具有趨勢(shì)性，還具有周期長(zhǎng)度大約為12的季節(jié)效應(yīng)。為消除趨勢(shì)性成分，首先對(duì)序列y1（t）進(jìn)行一階逐期差分，差分序列的線性圖像如圖1：

圖1 圖2

從圖中可以看出，差分后序列y1（t）的均值是穩(wěn)定在零點(diǎn)附近的，即通過(guò)一階差分，原序列的趨勢(shì)性基本消除了。

通過(guò)Eviews軟件得到的一階差分序列的樣本自相關(guān)系數(shù)圖形可知，當(dāng)滯后期k=12時(shí)，該序列的樣本自相關(guān)系數(shù)是-0.497，與0有顯著差異，這表明序列具有周期為12個(gè)月的季節(jié)波動(dòng)。再對(duì)y1（t）進(jìn)行一階季節(jié)差分得到季節(jié)差分序列y2（t），y2（t）的時(shí)序圖為圖2，自相關(guān)圖顯示延遲12步自相關(guān)系數(shù)顯著大于2倍標(biāo)準(zhǔn)差范圍，這說(shuō)明差分后序列中仍蘊(yùn)含著非常顯著的季節(jié)效應(yīng)。延遲1步，2步的自相關(guān)系數(shù)也大于2倍標(biāo)準(zhǔn)差，這說(shuō)明差分后序列還具有短期相關(guān)性，觀察偏相關(guān)圖得到的結(jié)論與此一致。

如果用AR模型和MA模型來(lái)模擬y2（t）的話，我們可以由擬合的殘差圖看出模型殘差非白噪聲，模型信息提取不充分，這里就不贅述了。

考慮到該序列既具有短期相關(guān)性又具有季節(jié)效應(yīng)，短期相關(guān)性和季節(jié)效應(yīng)不能簡(jiǎn)單地，可加性地提取，因而估計(jì)該系列的季節(jié)效應(yīng)和短期相關(guān)性之間具有復(fù)雜的關(guān)聯(lián)性。這是，通常假定短期相關(guān)性和季節(jié)效應(yīng)之間具有乘積關(guān)系，故我們用乘積季節(jié)模型來(lái)擬合。

三、建立模型

在選擇模型時(shí)ARIMA（p，d，q）×（P，D，Q）S，通常需要對(duì)觀測(cè)數(shù)據(jù)進(jìn)行對(duì)數(shù)變換，序列y（t）經(jīng)過(guò)對(duì)數(shù)變換記為序列l(wèi)ogy（t），由的自相關(guān)圖與偏相關(guān)圖（見(jiàn)圖2）可以知道：除了個(gè)別點(diǎn)外，自相關(guān)系數(shù)在1階結(jié)尾，偏相關(guān)系數(shù)在2階結(jié)尾。根據(jù)Box-Jenkins建模思想，我們可以嘗試用，模型1：ARIMA（2，1，0）×（1，1，1）12，模型2：ARIMA（2，1，1）×（1，1，1）12，模型3：ARIMA（1，1，0）×（1，1，1）12，模型4：ARIMA（1，1，1）×（1，1，1）12來(lái)擬合。

相應(yīng)的輸出結(jié)果為分別如下：

從此可以看出：根據(jù)AIC定階原則，選擇AIC值最小的模型。因此，確定ARIMA（1，1，1）（1，1，1）12為該居民消費(fèi)價(jià)格指數(shù)的模型。由圖3可以得到ARIMA（1，1，1）（1，1，1）的預(yù)測(cè)模型是：

（1+0.538445B12）（1-0.874479B）×（1-B）（1-B12）logy（t）=（1-0.922111B12）εt

對(duì)此模擬模型進(jìn)行檢驗(yàn)，檢驗(yàn)結(jié)果顯示該模型順利通過(guò)殘差白噪聲檢驗(yàn)和參數(shù)顯著性檢驗(yàn)。說(shuō)明這個(gè)擬合充分提取了原數(shù)據(jù)的信息。

四、模型的預(yù)測(cè)

我們利用上述模型對(duì)2011年10月～12月的數(shù)據(jù)進(jìn)行預(yù)測(cè)，得到的預(yù)測(cè)結(jié)果為圖4：

圖4 圖3

由圖中可以看出：Theil不等相關(guān)系數(shù)為0.378332，說(shuō)明模型對(duì)數(shù)據(jù)的擬合還是很好的。因此這一模型對(duì)居民價(jià)格指數(shù)的預(yù)測(cè)有明顯的參考價(jià)值。

參考文獻(xiàn)

[1]王沁.時(shí)間序列分析及其應(yīng)用.成都西南交通大學(xué)出版社.

[2]楊位欽，顧蘭著.時(shí)間序列分析與動(dòng)態(tài)數(shù)據(jù)建模.北京：北京理工大學(xué)出版社，1998.

[3]王燕.應(yīng)用時(shí)間序列分析. 中國(guó)人民大學(xué)出版社.

[4]涂雄苓，黃月玲.旅游需求預(yù)測(cè)的ARIMA乘積季節(jié)模型構(gòu)建及實(shí)證分析[J].西財(cái)經(jīng)學(xué)院學(xué)報(bào)2011.

[5]程毛林.動(dòng)態(tài)數(shù)據(jù)擬合的疊合模型及其應(yīng)用[J].運(yùn)籌與管理. 2005.

[6]Mario Milanese，and Carlo Novara.Set Membership Prediction of Nonlinear Time Series. IEEE Transactions on Automatic Control.2005.

[7]Young，P.Data-based mechanistic modeling of Engineering Systems.Journal of Vibration and Control.1998.

【摘要】居民消費(fèi)價(jià)格指數(shù)是我國(guó)物價(jià)指數(shù)體系中極其重要的一個(gè)指數(shù)，主要反映消費(fèi)者支付商品和勞務(wù)的價(jià)格變化情況，也是一種度量通貨膨脹水平的工具，以百分比變化為表達(dá)形式。我國(guó)改革開(kāi)放以來(lái)，社會(huì)經(jīng)濟(jì)的各方面發(fā)生巨大的變換，居民消費(fèi)價(jià)格指數(shù)的變動(dòng)也顯示出自身的特點(diǎn)，本文主要是應(yīng)用非平穩(wěn)時(shí)間序列和數(shù)學(xué)統(tǒng)計(jì)軟件Eviews對(duì)歷年的居民消費(fèi)價(jià)格指數(shù)進(jìn)行相關(guān)分析、處理并建立模型，這有利于我們認(rèn)識(shí)它與社會(huì)經(jīng)濟(jì)發(fā)展相聯(lián)系的變動(dòng)規(guī)律，文章最后根據(jù)所建立的模型舉例進(jìn)行了預(yù)測(cè)。

【關(guān)鍵詞】時(shí)間序列 非平穩(wěn)序列 乘積季節(jié)模型 建模

一、相關(guān)知識(shí)

時(shí)間序列分析是一種動(dòng)態(tài)數(shù)據(jù)處理的統(tǒng)計(jì)方法。一個(gè)時(shí)間序列{xt}是長(zhǎng)期趨勢(shì)變動(dòng)Tt，季節(jié)效應(yīng)St，周期變動(dòng)Ct和不規(guī)則變動(dòng)因素It共同作用的結(jié)果。確定性時(shí)間序列分析方法的原理是用確定性函數(shù)對(duì)長(zhǎng)期趨勢(shì)變動(dòng)Tt，季節(jié)效應(yīng)St，周期變動(dòng)Ct進(jìn)行擬合，將非平穩(wěn)時(shí)間序列平穩(wěn)化。本文我們?yōu)榱似椒€(wěn)化引入乘積季節(jié)模型。

乘積季節(jié)模型為：

ψ（B）U（BS）（1-B）d（1-BS）DXt=θ（B）V（BS）εt

E（εt）=0，Var（εt）=σ2ε，E（εSεK）=0，（s≠k）

EXkεt=0，（？坌k

其中：

U（BS）=1-u1BS-u2B2s-…-uPBPS

V（BS）=1-v1BS-v2B2s-…-uQBQS

可以對(duì)不同周期的同一周期點(diǎn)之間的相關(guān)性進(jìn)行擬合：

ψ（B）=1-ψ1B-ψ2B2-…-ψPBP

θ（B）=1-θ1B-θ2B2-…-θqBq

用來(lái)消除同一周期不同周期點(diǎn)之間的相關(guān)性。

季節(jié)模型簡(jiǎn)記為ARIMA（p，d，q）×（P，D，Q）S其中p和q是消除同一周期不同周期點(diǎn)之間相關(guān)性的自回歸階數(shù)和移動(dòng)平均階數(shù)，P和Q是消除不同周期的同一周期點(diǎn)之間的相關(guān)性的自回歸階數(shù)和移動(dòng)平均階數(shù)，s為周期長(zhǎng)度，d為差分的階數(shù)，D為季節(jié)差分的階數(shù)。

二、數(shù)據(jù)預(yù)處理

本文選取的是2000年1月到2011年12月全國(guó)居民消費(fèi)價(jià)格指數(shù)的數(shù)據(jù)，記居民消費(fèi)價(jià)格指數(shù)為具體數(shù)據(jù)如下表1：

通過(guò)Eviews軟件可得到時(shí)間序列y（t）的樣本自相關(guān)系數(shù)圖形可以發(fā)現(xiàn)，自相關(guān)系數(shù)呈線性衰減，衰減速度緩慢，這說(shuō)明該序列y（t）含有一定的趨勢(shì)性。同時(shí)可以看出該序列不僅具有趨勢(shì)性，還具有周期長(zhǎng)度大約為12的季節(jié)效應(yīng)。為消除趨勢(shì)性成分，首先對(duì)序列y1（t）進(jìn)行一階逐期差分，差分序列的線性圖像如圖1：

圖1 圖2

從圖中可以看出，差分后序列y1（t）的均值是穩(wěn)定在零點(diǎn)附近的，即通過(guò)一階差分，原序列的趨勢(shì)性基本消除了。

通過(guò)Eviews軟件得到的一階差分序列的樣本自相關(guān)系數(shù)圖形可知，當(dāng)滯后期k=12時(shí)，該序列的樣本自相關(guān)系數(shù)是-0.497，與0有顯著差異，這表明序列具有周期為12個(gè)月的季節(jié)波動(dòng)。再對(duì)y1（t）進(jìn)行一階季節(jié)差分得到季節(jié)差分序列y2（t），y2（t）的時(shí)序圖為圖2，自相關(guān)圖顯示延遲12步自相關(guān)系數(shù)顯著大于2倍標(biāo)準(zhǔn)差范圍，這說(shuō)明差分后序列中仍蘊(yùn)含著非常顯著的季節(jié)效應(yīng)。延遲1步，2步的自相關(guān)系數(shù)也大于2倍標(biāo)準(zhǔn)差，這說(shuō)明差分后序列還具有短期相關(guān)性，觀察偏相關(guān)圖得到的結(jié)論與此一致。

如果用AR模型和MA模型來(lái)模擬y2（t）的話，我們可以由擬合的殘差圖看出模型殘差非白噪聲，模型信息提取不充分，這里就不贅述了。

考慮到該序列既具有短期相關(guān)性又具有季節(jié)效應(yīng)，短期相關(guān)性和季節(jié)效應(yīng)不能簡(jiǎn)單地，可加性地提取，因而估計(jì)該系列的季節(jié)效應(yīng)和短期相關(guān)性之間具有復(fù)雜的關(guān)聯(lián)性。這是，通常假定短期相關(guān)性和季節(jié)效應(yīng)之間具有乘積關(guān)系，故我們用乘積季節(jié)模型來(lái)擬合。

三、建立模型

在選擇模型時(shí)ARIMA（p，d，q）×（P，D，Q）S，通常需要對(duì)觀測(cè)數(shù)據(jù)進(jìn)行對(duì)數(shù)變換，序列y（t）經(jīng)過(guò)對(duì)數(shù)變換記為序列l(wèi)ogy（t），由的自相關(guān)圖與偏相關(guān)圖（見(jiàn)圖2）可以知道：除了個(gè)別點(diǎn)外，自相關(guān)系數(shù)在1階結(jié)尾，偏相關(guān)系數(shù)在2階結(jié)尾。根據(jù)Box-Jenkins建模思想，我們可以嘗試用，模型1：ARIMA（2，1，0）×（1，1，1）12，模型2：ARIMA（2，1，1）×（1，1，1）12，模型3：ARIMA（1，1，0）×（1，1，1）12，模型4：ARIMA（1，1，1）×（1，1，1）12來(lái)擬合。

相應(yīng)的輸出結(jié)果為分別如下：

從此可以看出：根據(jù)AIC定階原則，選擇AIC值最小的模型。因此，確定ARIMA（1，1，1）（1，1，1）12為該居民消費(fèi)價(jià)格指數(shù)的模型。由圖3可以得到ARIMA（1，1，1）（1，1，1）的預(yù)測(cè)模型是：

（1+0.538445B12）（1-0.874479B）×（1-B）（1-B12）logy（t）=（1-0.922111B12）εt

對(duì)此模擬模型進(jìn)行檢驗(yàn)，檢驗(yàn)結(jié)果顯示該模型順利通過(guò)殘差白噪聲檢驗(yàn)和參數(shù)顯著性檢驗(yàn)。說(shuō)明這個(gè)擬合充分提取了原數(shù)據(jù)的信息。

四、模型的預(yù)測(cè)

我們利用上述模型對(duì)2011年10月～12月的數(shù)據(jù)進(jìn)行預(yù)測(cè)，得到的預(yù)測(cè)結(jié)果為圖4：

圖4 圖3

由圖中可以看出：Theil不等相關(guān)系數(shù)為0.378332，說(shuō)明模型對(duì)數(shù)據(jù)的擬合還是很好的。因此這一模型對(duì)居民價(jià)格指數(shù)的預(yù)測(cè)有明顯的參考價(jià)值。

參考文獻(xiàn)

[1]王沁.時(shí)間序列分析及其應(yīng)用.成都西南交通大學(xué)出版社.

[2]楊位欽，顧蘭著.時(shí)間序列分析與動(dòng)態(tài)數(shù)據(jù)建模.北京：北京理工大學(xué)出版社，1998.

[3]王燕.應(yīng)用時(shí)間序列分析. 中國(guó)人民大學(xué)出版社.

[4]涂雄苓，黃月玲.旅游需求預(yù)測(cè)的ARIMA乘積季節(jié)模型構(gòu)建及實(shí)證分析[J].西財(cái)經(jīng)學(xué)院學(xué)報(bào)2011.

[5]程毛林.動(dòng)態(tài)數(shù)據(jù)擬合的疊合模型及其應(yīng)用[J].運(yùn)籌與管理. 2005.

[6]Mario Milanese，and Carlo Novara.Set Membership Prediction of Nonlinear Time Series. IEEE Transactions on Automatic Control.2005.

[7]Young，P.Data-based mechanistic modeling of Engineering Systems.Journal of Vibration and Control.1998.

【摘要】居民消費(fèi)價(jià)格指數(shù)是我國(guó)物價(jià)指數(shù)體系中極其重要的一個(gè)指數(shù)，主要反映消費(fèi)者支付商品和勞務(wù)的價(jià)格變化情況，也是一種度量通貨膨脹水平的工具，以百分比變化為表達(dá)形式。我國(guó)改革開(kāi)放以來(lái)，社會(huì)經(jīng)濟(jì)的各方面發(fā)生巨大的變換，居民消費(fèi)價(jià)格指數(shù)的變動(dòng)也顯示出自身的特點(diǎn)，本文主要是應(yīng)用非平穩(wěn)時(shí)間序列和數(shù)學(xué)統(tǒng)計(jì)軟件Eviews對(duì)歷年的居民消費(fèi)價(jià)格指數(shù)進(jìn)行相關(guān)分析、處理并建立模型，這有利于我們認(rèn)識(shí)它與社會(huì)經(jīng)濟(jì)發(fā)展相聯(lián)系的變動(dòng)規(guī)律，文章最后根據(jù)所建立的模型舉例進(jìn)行了預(yù)測(cè)。

【關(guān)鍵詞】時(shí)間序列 非平穩(wěn)序列 乘積季節(jié)模型 建模

一、相關(guān)知識(shí)

時(shí)間序列分析是一種動(dòng)態(tài)數(shù)據(jù)處理的統(tǒng)計(jì)方法。一個(gè)時(shí)間序列{xt}是長(zhǎng)期趨勢(shì)變動(dòng)Tt，季節(jié)效應(yīng)St，周期變動(dòng)Ct和不規(guī)則變動(dòng)因素It共同作用的結(jié)果。確定性時(shí)間序列分析方法的原理是用確定性函數(shù)對(duì)長(zhǎng)期趨勢(shì)變動(dòng)Tt，季節(jié)效應(yīng)St，周期變動(dòng)Ct進(jìn)行擬合，將非平穩(wěn)時(shí)間序列平穩(wěn)化。本文我們?yōu)榱似椒€(wěn)化引入乘積季節(jié)模型。

乘積季節(jié)模型為：

ψ（B）U（BS）（1-B）d（1-BS）DXt=θ（B）V（BS）εt

E（εt）=0，Var（εt）=σ2ε，E（εSεK）=0，（s≠k）

EXkεt=0，（？坌k

其中：

U（BS）=1-u1BS-u2B2s-…-uPBPS

V（BS）=1-v1BS-v2B2s-…-uQBQS

可以對(duì)不同周期的同一周期點(diǎn)之間的相關(guān)性進(jìn)行擬合：

ψ（B）=1-ψ1B-ψ2B2-…-ψPBP

θ（B）=1-θ1B-θ2B2-…-θqBq

用來(lái)消除同一周期不同周期點(diǎn)之間的相關(guān)性。

季節(jié)模型簡(jiǎn)記為ARIMA（p，d，q）×（P，D，Q）S其中p和q是消除同一周期不同周期點(diǎn)之間相關(guān)性的自回歸階數(shù)和移動(dòng)平均階數(shù)，P和Q是消除不同周期的同一周期點(diǎn)之間的相關(guān)性的自回歸階數(shù)和移動(dòng)平均階數(shù)，s為周期長(zhǎng)度，d為差分的階數(shù)，D為季節(jié)差分的階數(shù)。

二、數(shù)據(jù)預(yù)處理

本文選取的是2000年1月到2011年12月全國(guó)居民消費(fèi)價(jià)格指數(shù)的數(shù)據(jù)，記居民消費(fèi)價(jià)格指數(shù)為具體數(shù)據(jù)如下表1：

通過(guò)Eviews軟件可得到時(shí)間序列y（t）的樣本自相關(guān)系數(shù)圖形可以發(fā)現(xiàn)，自相關(guān)系數(shù)呈線性衰減，衰減速度緩慢，這說(shuō)明該序列y（t）含有一定的趨勢(shì)性。同時(shí)可以看出該序列不僅具有趨勢(shì)性，還具有周期長(zhǎng)度大約為12的季節(jié)效應(yīng)。為消除趨勢(shì)性成分，首先對(duì)序列y1（t）進(jìn)行一階逐期差分，差分序列的線性圖像如圖1：

圖1 圖2

從圖中可以看出，差分后序列y1（t）的均值是穩(wěn)定在零點(diǎn)附近的，即通過(guò)一階差分，原序列的趨勢(shì)性基本消除了。

通過(guò)Eviews軟件得到的一階差分序列的樣本自相關(guān)系數(shù)圖形可知，當(dāng)滯后期k=12時(shí)，該序列的樣本自相關(guān)系數(shù)是-0.497，與0有顯著差異，這表明序列具有周期為12個(gè)月的季節(jié)波動(dòng)。再對(duì)y1（t）進(jìn)行一階季節(jié)差分得到季節(jié)差分序列y2（t），y2（t）的時(shí)序圖為圖2，自相關(guān)圖顯示延遲12步自相關(guān)系數(shù)顯著大于2倍標(biāo)準(zhǔn)差范圍，這說(shuō)明差分后序列中仍蘊(yùn)含著非常顯著的季節(jié)效應(yīng)。延遲1步，2步的自相關(guān)系數(shù)也大于2倍標(biāo)準(zhǔn)差，這說(shuō)明差分后序列還具有短期相關(guān)性，觀察偏相關(guān)圖得到的結(jié)論與此一致。

如果用AR模型和MA模型來(lái)模擬y2（t）的話，我們可以由擬合的殘差圖看出模型殘差非白噪聲，模型信息提取不充分，這里就不贅述了。

考慮到該序列既具有短期相關(guān)性又具有季節(jié)效應(yīng)，短期相關(guān)性和季節(jié)效應(yīng)不能簡(jiǎn)單地，可加性地提取，因而估計(jì)該系列的季節(jié)效應(yīng)和短期相關(guān)性之間具有復(fù)雜的關(guān)聯(lián)性。這是，通常假定短期相關(guān)性和季節(jié)效應(yīng)之間具有乘積關(guān)系，故我們用乘積季節(jié)模型來(lái)擬合。

三、建立模型

在選擇模型時(shí)ARIMA（p，d，q）×（P，D，Q）S，通常需要對(duì)觀測(cè)數(shù)據(jù)進(jìn)行對(duì)數(shù)變換，序列y（t）經(jīng)過(guò)對(duì)數(shù)變換記為序列l(wèi)ogy（t），由的自相關(guān)圖與偏相關(guān)圖（見(jiàn)圖2）可以知道：除了個(gè)別點(diǎn)外，自相關(guān)系數(shù)在1階結(jié)尾，偏相關(guān)系數(shù)在2階結(jié)尾。根據(jù)Box-Jenkins建模思想，我們可以嘗試用，模型1：ARIMA（2，1，0）×（1，1，1）12，模型2：ARIMA（2，1，1）×（1，1，1）12，模型3：ARIMA（1，1，0）×（1，1，1）12，模型4：ARIMA（1，1，1）×（1，1，1）12來(lái)擬合。

相應(yīng)的輸出結(jié)果為分別如下：

從此可以看出：根據(jù)AIC定階原則，選擇AIC值最小的模型。因此，確定ARIMA（1，1，1）（1，1，1）12為該居民消費(fèi)價(jià)格指數(shù)的模型。由圖3可以得到ARIMA（1，1，1）（1，1，1）的預(yù)測(cè)模型是：

（1+0.538445B12）（1-0.874479B）×（1-B）（1-B12）logy（t）=（1-0.922111B12）εt

對(duì)此模擬模型進(jìn)行檢驗(yàn)，檢驗(yàn)結(jié)果顯示該模型順利通過(guò)殘差白噪聲檢驗(yàn)和參數(shù)顯著性檢驗(yàn)。說(shuō)明這個(gè)擬合充分提取了原數(shù)據(jù)的信息。

四、模型的預(yù)測(cè)

我們利用上述模型對(duì)2011年10月～12月的數(shù)據(jù)進(jìn)行預(yù)測(cè)，得到的預(yù)測(cè)結(jié)果為圖4：

圖4 圖3

由圖中可以看出：Theil不等相關(guān)系數(shù)為0.378332，說(shuō)明模型對(duì)數(shù)據(jù)的擬合還是很好的。因此這一模型對(duì)居民價(jià)格指數(shù)的預(yù)測(cè)有明顯的參考價(jià)值。

參考文獻(xiàn)

[1]王沁.時(shí)間序列分析及其應(yīng)用.成都西南交通大學(xué)出版社.

[2]楊位欽，顧蘭著.時(shí)間序列分析與動(dòng)態(tài)數(shù)據(jù)建模.北京：北京理工大學(xué)出版社，1998.

[3]王燕.應(yīng)用時(shí)間序列分析. 中國(guó)人民大學(xué)出版社.

[4]涂雄苓，黃月玲.旅游需求預(yù)測(cè)的ARIMA乘積季節(jié)模型構(gòu)建及實(shí)證分析[J].西財(cái)經(jīng)學(xué)院學(xué)報(bào)2011.

[5]程毛林.動(dòng)態(tài)數(shù)據(jù)擬合的疊合模型及其應(yīng)用[J].運(yùn)籌與管理. 2005.

[6]Mario Milanese，and Carlo Novara.Set Membership Prediction of Nonlinear Time Series. IEEE Transactions on Automatic Control.2005.

[7]Young，P.Data-based mechanistic modeling of Engineering Systems.Journal of Vibration and Control.1998.