• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    回歸法設(shè)定指數(shù)平滑模型的最優(yōu)預(yù)測(cè)參數(shù)

    2021-09-19 09:10:28徐齊利
    統(tǒng)計(jì)學(xué)報(bào) 2021年4期
    關(guān)鍵詞:枚舉法參數(shù)估計(jì)預(yù)測(cè)

    徐齊利

    (江西財(cái)經(jīng)大學(xué) 經(jīng)濟(jì)學(xué)院,江西 南昌330013)

    一、引言

    指數(shù)平滑預(yù)測(cè)法(Exponential Smoothing Prediction Method)是一種常用的時(shí)間序列預(yù)測(cè)技術(shù),主要由三個(gè)預(yù)測(cè)模型組成,即一次指數(shù)平滑預(yù)測(cè)模型、二次指數(shù)平滑預(yù)測(cè)模型和三次指數(shù)平滑預(yù)測(cè)模型。在預(yù)測(cè)實(shí)踐中,指數(shù)平滑法的應(yīng)用面臨兩個(gè)難題,一是平滑次數(shù)的設(shè)定(吳德會(huì),2007)[1],二是平滑參數(shù)的設(shè)置(吳德會(huì),2008)[2]。

    關(guān)于平滑次數(shù)的設(shè)定,目前主要憑經(jīng)驗(yàn)判斷。在時(shí)間序列中,若波動(dòng)占優(yōu),則主要考慮一次指數(shù)平滑模型(王長(zhǎng)江,2006)[3];若趨勢(shì)占優(yōu),則主要考慮三次指數(shù)平滑模型(馮金巧等,2007)[4];若波動(dòng)和趨勢(shì)平分秋色,則主要考慮二次指數(shù)平滑模型(Khairina et al.,2021)[5]。

    關(guān)于平滑參數(shù)的設(shè)置,現(xiàn)行的方法主要有經(jīng)驗(yàn)法、試算法、枚舉法。

    其一,經(jīng)驗(yàn)法。該方法對(duì)時(shí)間序列的波動(dòng)態(tài)勢(shì)、發(fā)展趨勢(shì)做出經(jīng)驗(yàn)性判斷,據(jù)此給出平滑參數(shù)設(shè)置的參考區(qū)間。(1)當(dāng)時(shí)間序列呈現(xiàn)較穩(wěn)定的水平趨勢(shì)時(shí),平滑參數(shù)可在0.05~0.20之間取值;(2)當(dāng)時(shí)間序列有波動(dòng),但長(zhǎng)期趨勢(shì)變化并不大時(shí),平滑參數(shù)可在0.1~0.4之間取值;(3)當(dāng)時(shí)間序列波動(dòng)很大,長(zhǎng)期趨勢(shì)變化幅度較大,呈現(xiàn)出上升或下降趨勢(shì)時(shí),平滑參數(shù)可在0.6~0.8之間取值;(4)當(dāng)時(shí)間序列呈明顯上升或下降趨勢(shì)時(shí),平滑參數(shù)可在0.6~1之間取值。在預(yù)測(cè)操作上,無(wú)論靜態(tài)預(yù)測(cè)還是動(dòng)態(tài)預(yù)測(cè),通常將所給區(qū)間的中間值選為具體的平滑參數(shù)值(崔世杰等,2016;Joseph,2019;Ani,2019)[6-8]。

    其二,試算法。該方法在經(jīng)驗(yàn)法給出平滑參數(shù)的參考區(qū)間內(nèi),選取少量幾個(gè)(通常為3個(gè))具體數(shù)值作為待選的平滑參數(shù),對(duì)每個(gè)待選的平滑參數(shù)分別試算其預(yù)測(cè)的均方誤差,選擇使預(yù)測(cè)的均方誤差最小的平滑參數(shù)作為實(shí)際應(yīng)用的預(yù)測(cè)參數(shù)。在預(yù)測(cè)操作上,無(wú)論靜態(tài)預(yù)測(cè)還是動(dòng)態(tài)預(yù)測(cè),通常在經(jīng)驗(yàn)法所選平滑參數(shù)值的左右兩邊再各取一個(gè)數(shù)值做第二輪優(yōu)選,在左中右三個(gè)待選數(shù)值中確定一個(gè)使預(yù)測(cè)的均方誤差最小的數(shù)值。然后,再在該最優(yōu)值的左右兩邊進(jìn)一步各選一個(gè)數(shù)值,做第三輪優(yōu)選。以此類推,循環(huán)下去,直至最優(yōu)的平滑參數(shù)收斂或預(yù)測(cè)的均方誤差收斂為止(高春雷等,2014;Hsieh et al.,2020;Shi et al.,2020)[9-11]。

    其三,枚舉法。該方法既放棄經(jīng)驗(yàn)法的先驗(yàn)判斷,又借鑒試算法的基本思路,待選的平滑參數(shù)是一個(gè)區(qū)間更長(zhǎng)的等差數(shù)列,該區(qū)間通常是0~1,對(duì)數(shù)列中的每個(gè)待選參數(shù)分別試算其預(yù)測(cè)的均方誤差,選擇使預(yù)測(cè)的均方誤差最小的平滑參數(shù)為實(shí)際應(yīng)用的預(yù)測(cè)參數(shù)。在預(yù)測(cè)操作上,無(wú)論靜態(tài)預(yù)測(cè)還是動(dòng)態(tài)預(yù)測(cè),根據(jù)業(yè)務(wù)需要,若對(duì)預(yù)測(cè)精度要求并不高,但對(duì)預(yù)測(cè)速度要求很高,則通常枚舉共9個(gè)數(shù)值組成的等差數(shù)列待選;若對(duì)預(yù)測(cè)精度要求較高,且對(duì)預(yù)測(cè)速度要求也較高,則通常枚舉共19個(gè)數(shù)值組成的等差數(shù)列待選;若對(duì)預(yù)測(cè)精度要求很高,但對(duì)預(yù)測(cè)速度要求并不高,則通常枚舉共39個(gè)數(shù)值組成的等差數(shù)列待選(羅辰輝等,2018;Seong,2020;Sulandari,2021)[12-14]。

    經(jīng)驗(yàn)法、試算法、枚舉法各有優(yōu)劣(李守金等,2018;黎鎖平、劉坤會(huì),2004;Gustriansyah,2019)[15-17],但三者有四個(gè)共同的弊端。(1)弊端1:非全域取值。平滑參數(shù)的取值范圍皆被限定在0~1的區(qū)間內(nèi)。學(xué)界和業(yè)界從沒(méi)有解釋過(guò)平滑參數(shù)為什么必須嚴(yán)格限定在0~1的范圍內(nèi),選擇0~1區(qū)間取值僅因當(dāng)初的主觀直覺(jué)和方便操作。根據(jù)指數(shù)平滑模型的結(jié)構(gòu)和性質(zhì),平滑參數(shù)完全可以在1以上取值,也完全可以在0以下取值。(2)弊端2:非連續(xù)取值。平滑參數(shù)本是一個(gè)連續(xù)型變量,三種方法均將其變更為離散型變量來(lái)操作。經(jīng)驗(yàn)法表面上是給出了取值區(qū)間,但在實(shí)踐操作上則是將所給區(qū)間的中間值作為平滑參數(shù)值,試算法和枚舉法則更是直接將連續(xù)變量進(jìn)行了離散化處理。(3)弊端3:非自適應(yīng)取值。平滑參數(shù)皆是人為賦值,預(yù)測(cè)人員自己列示出若干備選答案,然后自己從中挑一個(gè)最合適的答案。經(jīng)驗(yàn)法僅給出四個(gè)備選數(shù)值,試算法也僅是在此基礎(chǔ)上多增加了兩個(gè)備選數(shù)值,枚舉法無(wú)非是再多增加幾個(gè)數(shù)值待選。(4)弊端4:非最優(yōu)取值。平滑參數(shù)的取值皆不是最優(yōu),肯定不是全局最優(yōu),充其量是局部最優(yōu)。經(jīng)驗(yàn)法、試算法、枚舉法均按照殘差平方和最小或均方誤差最小的原則確定最優(yōu)平滑參數(shù),但其可選集皆是可行集的子集的子集,即實(shí)數(shù)集上0~1區(qū)間內(nèi)的若干數(shù)值。

    為克服經(jīng)驗(yàn)法、試算法和枚舉法的弊端,本文提出了回歸法。不同于陳武(2016)[18]對(duì)二次指數(shù)平滑預(yù)測(cè)模型回歸系數(shù)計(jì)算方法的探討,本文將回歸所得的平滑參數(shù)估計(jì)值作為指數(shù)平滑預(yù)測(cè)法的預(yù)測(cè)參數(shù)。其理論過(guò)程為:針對(duì)平滑次數(shù)分別為一次、二次、三次的指數(shù)平滑預(yù)測(cè)模型,對(duì)應(yīng)構(gòu)造平滑次數(shù)分別為一次、二次、三次的指數(shù)平滑隨機(jī)過(guò)程,以此建立平滑次數(shù)分別為一次、二次、三次的指數(shù)平滑回歸模型,隨之進(jìn)行平滑次數(shù)分別為一次、二次、三次的指數(shù)平滑參數(shù)估計(jì)。

    在實(shí)驗(yàn)?zāi)M中,事先設(shè)定平滑參數(shù)的真實(shí)值,由一次、二次、三次指數(shù)平滑隨機(jī)過(guò)程生成對(duì)應(yīng)的指數(shù)平滑時(shí)間序列,采用對(duì)應(yīng)的指數(shù)平滑回歸模型對(duì)相應(yīng)的指數(shù)平滑時(shí)間序列進(jìn)行回歸分析。參數(shù)估計(jì)的結(jié)果說(shuō)明,回歸法設(shè)定指數(shù)平滑模型的最優(yōu)預(yù)測(cè)參數(shù)是可行的。假設(shè)檢驗(yàn)的結(jié)果證明,回歸法設(shè)定指數(shù)平滑模型的最優(yōu)預(yù)測(cè)參數(shù)是可信的。

    在實(shí)證應(yīng)用中,將回歸法用于GDP和稅收的年度時(shí)間序列預(yù)測(cè),并與經(jīng)驗(yàn)法、試算法和枚舉法的預(yù)測(cè)效果進(jìn)行對(duì)比。在預(yù)測(cè)效果上,回歸法的預(yù)測(cè)效果確實(shí)能夠全面優(yōu)于經(jīng)驗(yàn)法、試算法和枚舉法。在預(yù)測(cè)技術(shù)上,回歸法確實(shí)能夠有效規(guī)避經(jīng)驗(yàn)法、試算法和枚舉法在平滑參數(shù)設(shè)置上存在的非全域取值、非連續(xù)取值、非自適應(yīng)取值、非最優(yōu)取值等弊端。

    二、指數(shù)平滑預(yù)測(cè)模型

    (一)一次指數(shù)平滑預(yù)測(cè)模型

    對(duì)時(shí)間序列{Xt}t=1,2,…進(jìn)行一次指數(shù)平滑預(yù)測(cè)的模型為:

    其中,預(yù)測(cè)參數(shù)α需事前設(shè)定。回歸法設(shè)定該一次指數(shù)平滑模型最優(yōu)預(yù)測(cè)參數(shù)的邏輯過(guò)程分為三步:第一步,根據(jù)一次指數(shù)平滑預(yù)測(cè)模型構(gòu)造一次指數(shù)平滑隨機(jī)過(guò)程;第二步,基于一次指數(shù)平滑隨機(jī)過(guò)程建立一次指數(shù)平滑回歸模型;第三步,估計(jì)一次指數(shù)平滑回歸模型的回歸系數(shù),以獲得最優(yōu)的預(yù)測(cè)參數(shù)。

    (二)二次指數(shù)平滑預(yù)測(cè)模型

    對(duì)時(shí)間序列{Xt}t=1,2,…進(jìn)行二次指數(shù)平滑預(yù)測(cè)的模型為:

    其中,預(yù)測(cè)參數(shù)α也需事前設(shè)定?;貧w法設(shè)定該二次指數(shù)平滑模型最優(yōu)預(yù)測(cè)參數(shù)的邏輯過(guò)程分為三步:第一步,根據(jù)二次指數(shù)平滑預(yù)測(cè)模型構(gòu)造二次指數(shù)平滑隨機(jī)過(guò)程;第二步,基于二次指數(shù)平滑隨機(jī)過(guò)程建立二次指數(shù)平滑回歸模型;第三步,估計(jì)二次指數(shù)平滑回歸模型的回歸系數(shù),以獲得最優(yōu)的預(yù)測(cè)參數(shù)。

    (三)三次指數(shù)平滑預(yù)測(cè)模型

    對(duì)時(shí)間序列{Xt}t=1,2,…進(jìn)行三次指數(shù)平滑預(yù)測(cè)的模型為:

    其中,預(yù)測(cè)參數(shù)α仍需事前設(shè)定?;貧w法設(shè)定該三次指數(shù)平滑模型最優(yōu)預(yù)測(cè)參數(shù)的邏輯過(guò)程分為三步:第一步,根據(jù)三次指數(shù)平滑預(yù)測(cè)模型構(gòu)造三次指數(shù)平滑隨機(jī)過(guò)程;第二步,基于三次指數(shù)平滑隨機(jī)過(guò)程建立三次指數(shù)平滑回歸模型;第三步,估計(jì)三次指數(shù)平滑回歸模型的回歸系數(shù),以獲得最優(yōu)的預(yù)測(cè)參數(shù)。

    三、指數(shù)平滑隨機(jī)過(guò)程

    (一)一次指數(shù)平滑隨機(jī)過(guò)程

    回歸法設(shè)定一次指數(shù)平滑模型最優(yōu)預(yù)測(cè)參數(shù)的邏輯過(guò)程的第一步:根據(jù)一次指數(shù)平滑預(yù)測(cè)模型構(gòu)造一次指數(shù)平滑隨機(jī)過(guò)程。為使(1)式一次指數(shù)平滑模型的預(yù)測(cè)值X?t+1與該模型的條件期望E(Xt+1|Xt,…,X1;α)建立起等價(jià)關(guān)系E(Xt+1|Xt,…,X1;α)=X?t+1,特建立如下一次指數(shù)平滑隨機(jī)過(guò)程為:

    如此設(shè)置一次指數(shù)平滑隨機(jī)過(guò)程,則預(yù)測(cè)參數(shù)α可以通過(guò)建立回歸模型進(jìn)行估計(jì)。

    (二)二次指數(shù)平滑隨機(jī)過(guò)程

    回歸法設(shè)定二次指數(shù)平滑模型最優(yōu)預(yù)測(cè)參數(shù)的邏輯過(guò)程的第一步:根據(jù)二次指數(shù)平滑預(yù)測(cè)模型構(gòu)造二次指數(shù)平滑隨機(jī)過(guò)程。為使(2)式二次指數(shù)平滑模型的預(yù)測(cè)值1與該模型的條件期望E(Xt+1|Xt,…,X1;α)建立起等價(jià)關(guān)系E(Xt+1|Xt,…,X1;α)=X?t+1,特建立如下二次指數(shù)平滑隨機(jī)過(guò)程為:

    如此設(shè)置二次指數(shù)平滑隨機(jī)過(guò)程,則預(yù)測(cè)參數(shù)α可以通過(guò)建立回歸模型進(jìn)行估計(jì)。

    (三)三次指數(shù)平滑隨機(jī)過(guò)程

    回歸法設(shè)定三次指數(shù)平滑模型最優(yōu)預(yù)測(cè)參數(shù)的邏輯過(guò)程的第一步:根據(jù)三次指數(shù)平滑預(yù)測(cè)模型構(gòu)造三次指數(shù)平滑隨機(jī)過(guò)程。為使(3)式三次指數(shù)平滑模型的預(yù)測(cè)值X?t+1與該模型的條件期望E(Xt+1|Xt,…,X1;α)建立起等價(jià)關(guān)系E(Xt+1|Xt,…,X1;α)=X?t+1,特建立如下三次指數(shù)平滑隨機(jī)過(guò)程為:

    如此設(shè)置三次指數(shù)平滑隨機(jī)過(guò)程,則預(yù)測(cè)參數(shù)α可以通過(guò)建立回歸模型進(jìn)行估計(jì)。

    四、指數(shù)平滑回歸模型

    (一)一次指數(shù)平滑回歸模型

    回歸法設(shè)定一次指數(shù)平滑模型最優(yōu)預(yù)測(cè)參數(shù)的邏輯過(guò)程的第二步:基于一次指數(shù)平滑隨機(jī)過(guò)程建立一次指數(shù)平滑回歸模型。針對(duì)由(4)式一次指數(shù)平滑隨機(jī)過(guò)程生成的時(shí)間序列{X1,X2,…,Xn},可建立如下回歸模型來(lái)估計(jì)未知參數(shù)α:

    模型中,直接決定可觀測(cè)時(shí)間序列{X1,X2,…,Xn}基本走勢(shì)的平滑因子{S1,S2,…,Sn}是不可直接觀測(cè)的時(shí)間序列。

    (二)二次指數(shù)平滑回歸模型

    回歸法設(shè)定二次指數(shù)平滑模型最優(yōu)預(yù)測(cè)參數(shù)的邏輯過(guò)程的第二步:基于二次指數(shù)平滑隨機(jī)過(guò)程建立二次指數(shù)平滑回歸模型。針對(duì)由(5)式二次指數(shù)平滑隨機(jī)過(guò)程生成的時(shí)間序列{X1,X2,…,Xn},可建立如下回歸模型來(lái)估計(jì)未知參數(shù)α:

    模型中,直接決定可觀測(cè)時(shí)間序列{X1,X2,…,Xn}基本走勢(shì)的平滑因子是不可直接觀測(cè)的時(shí)間序列。

    (三)三次指數(shù)平滑回歸模型

    回歸法設(shè)定三次指數(shù)平滑模型最優(yōu)預(yù)測(cè)參數(shù)的邏輯過(guò)程的第二步:基于三次指數(shù)平滑隨機(jī)過(guò)程建立三次指數(shù)平滑回歸模型。針對(duì)由(6)式三次指數(shù)平滑隨機(jī)過(guò)程生成的時(shí)間序列{X1,X2,…,Xn},可建立如下回歸模型來(lái)估計(jì)未知參數(shù)α:

    模型中,直接決定可觀測(cè)時(shí)間序列{X1,X2,…,Xn}基本走勢(shì)的平滑因子是不可直接觀測(cè)的時(shí)間序列。

    五、指數(shù)平滑參數(shù)估計(jì)

    (一)一次指數(shù)平滑參數(shù)估計(jì)

    對(duì)時(shí)間序列{X1,X2,…,Xn}的一次指數(shù)平滑模型(7)式、二次指數(shù)平滑模型(8)式、三次指數(shù)平滑模型(9)式,皆采取最小二乘的方法進(jìn)行估計(jì),從而得到最優(yōu)預(yù)測(cè)參數(shù)α?為:

    求解非線性最小二乘,通用的迭代算法可參考Seber和Wild(2003)的研究[19]。

    以往n-2期預(yù)測(cè)的均方誤差MSE是隨機(jī)擾動(dòng)項(xiàng)εt+1~N(0,σ2)方差σ2的估計(jì)值,即:

    其中,n-2維雅閣比向量J=(J3,…,Jn)T的分量Jt+1為:

    與通常的非線性回歸模型不同,本文所建的回歸模型其雅閣比向量各分量之間存在遞推關(guān)系。

    回歸法設(shè)定一次指數(shù)平滑模型最優(yōu)預(yù)測(cè)參數(shù)的邏輯過(guò)程的第三步:估計(jì)一次指數(shù)平滑回歸模型的回歸系數(shù),以獲得最優(yōu)的預(yù)測(cè)參數(shù)。對(duì)一次指數(shù)平滑回歸模型(7)式的參數(shù)α進(jìn)行最小二乘估計(jì),算法在迭代過(guò)程中雅閣比向量各分量之間的遞推關(guān)系為:

    (二)二次指數(shù)平滑參數(shù)估計(jì)

    回歸法設(shè)定二次指數(shù)平滑模型最優(yōu)預(yù)測(cè)參數(shù)的邏輯過(guò)程的第三步:估計(jì)二次指數(shù)平滑回歸模型的回歸系數(shù),以獲得最優(yōu)的預(yù)測(cè)參數(shù)。對(duì)二次指數(shù)平滑回歸模型(8)式的參數(shù)α進(jìn)行最小二乘估計(jì),算法在迭代過(guò)程中雅閣比向量各分量之間的遞推關(guān)系為:

    (三)三次指數(shù)平滑參數(shù)估計(jì)

    回歸法設(shè)定三次指數(shù)平滑模型最優(yōu)預(yù)測(cè)參數(shù)的邏輯過(guò)程的第三步:估計(jì)三次指數(shù)平滑回歸模型的回歸系數(shù),以獲得最優(yōu)的預(yù)測(cè)參數(shù)。對(duì)三次指數(shù)平滑回歸模型(9)式的參數(shù)α進(jìn)行最小二乘估計(jì),算法在迭代過(guò)程中雅閣比向量各分量之間的遞推關(guān)系為:

    六、實(shí)驗(yàn)?zāi)M

    (一)一次指數(shù)平滑時(shí)間序列

    為檢驗(yàn)上述回歸法設(shè)定一次指數(shù)平滑模型的最優(yōu)預(yù)測(cè)參數(shù)是否可行、是否可信,需采用蒙特卡洛方法做模擬實(shí)驗(yàn)。針對(duì)(4)式的一次指數(shù)平滑隨機(jī)過(guò)程,令初值X1=1 000.0,X2=1 025.0,參數(shù)α=0.10,生成兩個(gè)樣本容量n=1 000的一次指數(shù)平滑時(shí)間序列,如圖1所示。

    圖1 一次指數(shù)平滑隨機(jī)過(guò)程的實(shí)現(xiàn)

    對(duì)這兩個(gè)一次指數(shù)平滑時(shí)間序列,按照(7)式的一次指數(shù)平滑回歸模型,估計(jì)其參數(shù)α,結(jié)果如表1所示。

    表1 一次指數(shù)平滑時(shí)間序列的估計(jì)結(jié)果

    表1顯示,一次指數(shù)平滑回歸模型的估計(jì)結(jié)果與參數(shù)實(shí)際值非常接近,這從參數(shù)估計(jì)上說(shuō)明,針對(duì)一次指數(shù)平滑時(shí)間序列,通過(guò)回歸的方法設(shè)定其最優(yōu)預(yù)測(cè)參數(shù)是可行的。根據(jù)表1的參數(shù)估計(jì)結(jié)果,構(gòu)造原假設(shè)H0、備擇假設(shè)H1以及檢驗(yàn)統(tǒng)計(jì)量τ,分別為:

    其中,Sα為α的標(biāo)準(zhǔn)誤,以此考察估計(jì)結(jié)果的置信水平。表1顯示,兩序列檢驗(yàn)統(tǒng)計(jì)量的取值皆充分接近0,在自由度n-3=997的t分布下以高置信水平支持原假設(shè)H0:α=0.10。這從假設(shè)檢驗(yàn)上證明,針對(duì)一次指數(shù)平滑時(shí)間序列,通過(guò)回歸的方法設(shè)定其最優(yōu)預(yù)測(cè)參數(shù)是可信的。

    (二)二次指數(shù)平滑時(shí)間序列

    為檢驗(yàn)上述回歸法設(shè)定二次指數(shù)平滑模型的最優(yōu)預(yù)測(cè)參數(shù)是否可行、是否可信,也需采用蒙特卡洛方法做模擬實(shí)驗(yàn)。針對(duì)(5)式的二次指數(shù)平滑隨機(jī)過(guò)程,令初值X1=1 000.0,X2=1 025.0,參數(shù)α=0.10,生成兩個(gè)樣本容量n=1 000的二次指數(shù)平滑時(shí)間序列,如圖2所示。

    圖2 二次指數(shù)平滑隨機(jī)過(guò)程的實(shí)現(xiàn)

    對(duì)這兩個(gè)二次指數(shù)平滑時(shí)間序列,按照(8)式的二次指數(shù)平滑回歸模型,估計(jì)其參數(shù)α,結(jié)果如表2所示。

    表2 二次指數(shù)平滑時(shí)間序列的估計(jì)結(jié)果

    表2顯示,二次指數(shù)平滑回歸模型的估計(jì)結(jié)果與參數(shù)實(shí)際值非常接近,這從參數(shù)估計(jì)上說(shuō)明,針對(duì)二次指數(shù)平滑時(shí)間序列,通過(guò)回歸的方法設(shè)定其最優(yōu)預(yù)測(cè)參數(shù)是可行的。根據(jù)表2的參數(shù)估計(jì)結(jié)果,按照(10)式進(jìn)行假設(shè)檢驗(yàn),以此考察估計(jì)結(jié)果的置信水平。表2進(jìn)一步顯示,兩序列檢驗(yàn)統(tǒng)計(jì)量的取值皆充分接近0,在自由度n-3=997的t分布下以高置信水平支持原假設(shè)H0:α=0.10。這從假設(shè)檢驗(yàn)上證明,針對(duì)二次指數(shù)平滑時(shí)間序列,通過(guò)回歸的方法設(shè)定其最優(yōu)預(yù)測(cè)參數(shù)是可信的。

    (三)三次指數(shù)平滑時(shí)間序列

    為檢驗(yàn)上述回歸法設(shè)定三次指數(shù)平滑模型的最優(yōu)預(yù)測(cè)參數(shù)是否可行、是否可信,仍需采用蒙特卡洛方法做模擬實(shí)驗(yàn)。針對(duì)(6)式的三次指數(shù)平滑隨機(jī)過(guò)程,令初值X1=1 000.0,X2=1 025.0,參數(shù)α=0.10,生成兩個(gè)樣本容量n=1 000的三次指數(shù)平滑時(shí)間序列,如圖3所示。

    對(duì)這兩個(gè)三次指數(shù)平滑時(shí)間序列,按照(9)式的三次指數(shù)平滑回歸模型,估計(jì)其參數(shù)α,結(jié)果如表3所示。

    表3顯示,三次指數(shù)平滑回歸模型的估計(jì)結(jié)果與參數(shù)實(shí)際值都非常接近,這從參數(shù)估計(jì)上說(shuō)明,針對(duì)三次指數(shù)平滑時(shí)間序列,通過(guò)回歸的方法設(shè)定其最優(yōu)預(yù)測(cè)參數(shù)是可行的。根據(jù)表3中的參數(shù)估計(jì)結(jié)果,按照(10)式進(jìn)行假設(shè)檢驗(yàn),以此考察估計(jì)結(jié)果的置信水平。表3進(jìn)一步顯示,兩序列檢驗(yàn)統(tǒng)計(jì)量的取值皆充分接近0,在自由度n-3=997的t分布下以高置信水平支持原假設(shè)H0:α=0.10。這從假設(shè)檢驗(yàn)上證明,針對(duì)三次指數(shù)平滑時(shí)間序列,通過(guò)回歸的方法設(shè)定其最優(yōu)預(yù)測(cè)參數(shù)是可信的。

    圖3 三次指數(shù)平滑隨機(jī)過(guò)程的實(shí)現(xiàn)

    表3 三次指數(shù)平滑時(shí)間序列的估計(jì)結(jié)果

    七、實(shí)證應(yīng)用

    (一)GDP預(yù)測(cè)

    既然實(shí)驗(yàn)?zāi)M顯示回歸法設(shè)定指數(shù)平滑模型的預(yù)測(cè)參數(shù)是可行和可信的,那么接下來(lái)的工作是實(shí)證檢驗(yàn)在預(yù)測(cè)實(shí)踐中回歸法較之已有的經(jīng)驗(yàn)法、試算法、枚舉法的邊際進(jìn)步性。實(shí)證目標(biāo)之一是:在平滑參數(shù)的設(shè)定上,檢驗(yàn)回歸法是否確實(shí)能夠有效規(guī)避經(jīng)驗(yàn)法、試算法、枚舉法的四個(gè)弊端。實(shí)證目標(biāo)之二是:在指數(shù)平滑的預(yù)測(cè)效果上,檢驗(yàn)回歸法是否確實(shí)明顯優(yōu)于經(jīng)驗(yàn)法、試算法、枚舉法。雖然前面的邏輯過(guò)程已經(jīng)在理論上做了肯定性回答,但現(xiàn)實(shí)應(yīng)用效果到底如何、能否達(dá)到理論所述的預(yù)期效果,需要做出驗(yàn)證性回答。俗話說(shuō)“是騾子是馬,拉出來(lái)遛遛”,哲學(xué)所言“實(shí)踐是檢驗(yàn)真理的唯一標(biāo)準(zhǔn)”,講的就是這個(gè)道理。

    以GDP的年度預(yù)測(cè)為例,針對(duì)表4的GDP時(shí)間序列數(shù)據(jù),分別按照(7)式的一次指數(shù)平滑回歸模型、(8)式的二次指數(shù)平滑回歸模型、(9)式的三次指數(shù)平滑回歸模型,估計(jì)各自的參數(shù)α,結(jié)果如表5所示。

    表4 對(duì)GDP進(jìn)行二次指數(shù)平滑和三次指數(shù)平滑的預(yù)測(cè)結(jié)果

    (續(xù)表4)

    表5顯示,三個(gè)模型的回歸系數(shù)皆顯著非零,說(shuō)明一次平滑模型、二次平滑模型、三次平滑模型均可用于預(yù)測(cè)GDP的時(shí)間序列走勢(shì)。進(jìn)一步考察預(yù)測(cè)效果,利用均方誤差MSE構(gòu)造F檢驗(yàn)統(tǒng)計(jì)量MSEi/MSEj~F(n-3,n-3),可從邏輯上做出判斷,利用指數(shù)平滑模型預(yù)測(cè)GDP的時(shí)間序列走勢(shì)時(shí),一次平滑模型的預(yù)測(cè)效果可能較差,二次平滑模型和三次平滑模型的預(yù)測(cè)效果基本相當(dāng)。

    利用表5的參數(shù)取值,分別按照(1)式的一次指數(shù)平滑預(yù)測(cè)模型、(2)式的二次指數(shù)平滑預(yù)測(cè)模型、(3)式的三次指數(shù)平滑預(yù)測(cè)模型,對(duì)GDP的走勢(shì)進(jìn)行向前一步預(yù)測(cè),結(jié)果如表4所示。表4顯示,與回歸分析時(shí)的判斷一致,對(duì)GDP時(shí)間序列進(jìn)行指數(shù)平滑預(yù)測(cè)時(shí),一次平滑模型的預(yù)測(cè)效果較差,二次平滑模型和三次平滑模型的預(yù)測(cè)效果基本相當(dāng)。

    表5 對(duì)GDP進(jìn)行指數(shù)平滑預(yù)測(cè)的模型參數(shù)設(shè)置

    從表5的參數(shù)估計(jì)值可以看出,對(duì)GDP進(jìn)行指數(shù)平滑預(yù)測(cè),在平滑參數(shù)的設(shè)定上,回歸法確實(shí)能夠有效規(guī)避經(jīng)驗(yàn)法、試算法、枚舉法的弊端。圖4的黑點(diǎn)顯示回歸法的預(yù)測(cè)效果,白點(diǎn)顯示經(jīng)驗(yàn)法、試算法與枚舉法的預(yù)測(cè)效果,可以看出,對(duì)GDP進(jìn)行指數(shù)平滑預(yù)測(cè)時(shí),在預(yù)測(cè)效果的對(duì)比上,平滑參數(shù)即便在0~1區(qū)間內(nèi)取值,經(jīng)驗(yàn)法、試算法、枚舉法也難以達(dá)到最優(yōu)預(yù)測(cè),而回歸法幾乎是最優(yōu)預(yù)測(cè)。

    圖4 回歸法與經(jīng)驗(yàn)法、試算法、枚舉法對(duì)GDP進(jìn)行三次指數(shù)平滑預(yù)測(cè)的效果比較

    (二)稅收預(yù)測(cè)

    以稅收的年度預(yù)測(cè)為例,針對(duì)表6的稅收時(shí)間序列數(shù)據(jù),分別按照(7)式的一次指數(shù)平滑回歸模型、(8)式的二次指數(shù)平滑回歸模型、(9)式的三次指數(shù)平滑回歸模型,估計(jì)各自的參數(shù)α,結(jié)果如表7所示。

    表7顯示,三個(gè)模型的回歸系數(shù)皆顯著非零,說(shuō)明一次平滑模型、二次平滑模型、三次平滑模型均可用于預(yù)測(cè)稅收的時(shí)間序列走勢(shì)。進(jìn)一步考察預(yù)測(cè)效果,利用均方誤差MSE,構(gòu)造F檢驗(yàn)統(tǒng)計(jì)量MSEi/MSEj~F(n-3,n-3),可從邏輯上做出判斷,利用指數(shù)平滑模型用于預(yù)測(cè)稅收的時(shí)間序列走勢(shì)時(shí),一次平滑模型的預(yù)測(cè)效果可能較差,二次平滑模型和三次平滑模型的預(yù)測(cè)效果基本相當(dāng)。

    表6 對(duì)稅收進(jìn)行二次指數(shù)平滑和三次指數(shù)平滑的預(yù)測(cè)結(jié)果

    (續(xù)表6)

    表7 對(duì)稅收進(jìn)行指數(shù)平滑預(yù)測(cè)的模型參數(shù)設(shè)置

    利用表7的參數(shù)取值,分別按照(1)式的一次指數(shù)平滑預(yù)測(cè)模型、(2)式的二次指數(shù)平滑預(yù)測(cè)模型、(3)式的三次指數(shù)平滑預(yù)測(cè)模型,對(duì)稅收的走勢(shì)進(jìn)行向前一步預(yù)測(cè),結(jié)果如表6所示。表6顯示,與回歸分析時(shí)的判斷一致,在對(duì)稅收時(shí)間序列進(jìn)行指數(shù)平滑預(yù)測(cè)時(shí),一次平滑模型的預(yù)測(cè)效果較差,二次平滑模型和三次平滑模型的預(yù)測(cè)效果基本相當(dāng)。

    圖5 回歸法與經(jīng)驗(yàn)法、試算法、枚舉法對(duì)稅收進(jìn)行三次指數(shù)平滑預(yù)測(cè)的效果比較

    從表7的參數(shù)估計(jì)值可以看出,對(duì)稅收進(jìn)行指數(shù)平滑預(yù)測(cè),在平滑參數(shù)的設(shè)定上,回歸法確實(shí)能夠有效規(guī)避經(jīng)驗(yàn)法、試算法、枚舉法的弊端。圖5的黑點(diǎn)顯示回歸法的預(yù)測(cè)效果,白點(diǎn)顯示經(jīng)驗(yàn)法、試算法與枚舉法的預(yù)測(cè)效果,可以看出,對(duì)稅收進(jìn)行指數(shù)平滑預(yù)測(cè)時(shí),在預(yù)測(cè)效果的對(duì)比上,平滑參數(shù)即便在0~1區(qū)間內(nèi)取值,經(jīng)驗(yàn)法、試算法、枚舉法也難以達(dá)到最優(yōu)預(yù)測(cè),而回歸法幾乎是最優(yōu)預(yù)測(cè)。

    八、結(jié)論

    對(duì)于指數(shù)平滑預(yù)測(cè)模型的參數(shù)設(shè)置方法,為克服現(xiàn)行的經(jīng)驗(yàn)法、試算法和枚舉法的弊端,本文提出了回歸法。將回歸所得的平滑參數(shù)估計(jì)值作為指數(shù)平滑預(yù)測(cè)法的預(yù)測(cè)參數(shù),其理論過(guò)程為:針對(duì)一次指數(shù)平滑預(yù)測(cè)模型、二次指數(shù)平滑預(yù)測(cè)模型、三次指數(shù)平滑預(yù)測(cè)模型,對(duì)應(yīng)構(gòu)造一次指數(shù)平滑隨機(jī)過(guò)程、二次指數(shù)平滑隨機(jī)過(guò)程、三次指數(shù)平滑隨機(jī)過(guò)程,以此建立一次指數(shù)平滑回歸模型、二次指數(shù)平滑回歸模型、三次指數(shù)平滑回歸模型,隨之進(jìn)行一次指數(shù)平滑參數(shù)估計(jì)、二次指數(shù)平滑參數(shù)估計(jì)、三次指數(shù)平滑參數(shù)估計(jì)。

    實(shí)驗(yàn)?zāi)M中,事先給定平滑參數(shù)真實(shí)值,由一次、二次、三次指數(shù)平滑隨機(jī)過(guò)程的一次、二次、三次指數(shù)平滑時(shí)間序列進(jìn)行回歸分析。參數(shù)估計(jì)結(jié)果皆說(shuō)明,回歸法設(shè)定指數(shù)平滑模型對(duì)應(yīng)生成一次、二次、三次指數(shù)平滑時(shí)間序列,利用一次、二次、三次指數(shù)平滑回歸模型估計(jì)相應(yīng)的最優(yōu)預(yù)測(cè)參數(shù)是可行的。假設(shè)檢驗(yàn)結(jié)果證明,回歸法設(shè)定指數(shù)平滑模型的最優(yōu)預(yù)測(cè)參數(shù)是可信的。

    實(shí)證應(yīng)用中,將回歸法用于GDP和稅收年度時(shí)間序列的一次、二次、三次指數(shù)平滑預(yù)測(cè),并與經(jīng)驗(yàn)法、試算法和枚舉法的預(yù)測(cè)效果進(jìn)行對(duì)比。結(jié)果顯示,在預(yù)測(cè)效果上,回歸法的預(yù)測(cè)效果確實(shí)能夠全面優(yōu)于經(jīng)驗(yàn)法、試算法和枚舉法。在預(yù)測(cè)技術(shù)上,回歸法確實(shí)能夠有效規(guī)避經(jīng)驗(yàn)法、試算法和枚舉法在平滑參數(shù)設(shè)置上存在的非全域取值、非連續(xù)取值、非自適應(yīng)取值、非最優(yōu)取值等弊端。

    猜你喜歡
    枚舉法參數(shù)估計(jì)預(yù)測(cè)
    無(wú)可預(yù)測(cè)
    黃河之聲(2022年10期)2022-09-27 13:59:46
    選修2-2期中考試預(yù)測(cè)卷(A卷)
    選修2-2期中考試預(yù)測(cè)卷(B卷)
    基于新型DFrFT的LFM信號(hào)參數(shù)估計(jì)算法
    枚舉法的程序?qū)崿F(xiàn)及優(yōu)化
    應(yīng)重視用枚舉法解題
    不必預(yù)測(cè)未來(lái),只需把握現(xiàn)在
    Logistic回歸模型的幾乎無(wú)偏兩參數(shù)估計(jì)
    基于向前方程的平穩(wěn)分布參數(shù)估計(jì)
    基于競(jìng)爭(zhēng)失效數(shù)據(jù)的Lindley分布參數(shù)估計(jì)
    波多野结衣av一区二区av| 少妇 在线观看| 最新美女视频免费是黄的| 午夜精品在线福利| 免费人成视频x8x8入口观看| 久久中文字幕一级| 日韩精品免费视频一区二区三区| 欧美日韩亚洲高清精品| 一级a爱视频在线免费观看| 精品久久久久久久久久免费视频 | 精品午夜福利视频在线观看一区| 三级毛片av免费| 亚洲一卡2卡3卡4卡5卡精品中文| 电影成人av| 亚洲狠狠婷婷综合久久图片| 十八禁网站免费在线| 国产精品香港三级国产av潘金莲| 国产亚洲av高清不卡| 午夜福利影视在线免费观看| 最近最新中文字幕大全电影3 | 午夜福利影视在线免费观看| 精品欧美一区二区三区在线| 欧美黑人欧美精品刺激| 国产欧美日韩综合在线一区二区| 亚洲熟妇中文字幕五十中出 | 欧美日韩一级在线毛片| 色综合站精品国产| 日韩欧美一区二区三区在线观看| avwww免费| 新久久久久国产一级毛片| 成人av一区二区三区在线看| 91九色精品人成在线观看| 欧美最黄视频在线播放免费 | 一本大道久久a久久精品| 黑人操中国人逼视频| 国产成人免费无遮挡视频| 精品久久久久久久久久免费视频 | 欧美日韩av久久| 成人永久免费在线观看视频| 好看av亚洲va欧美ⅴa在| 黄色怎么调成土黄色| 午夜成年电影在线免费观看| 黄色片一级片一级黄色片| 久久久久九九精品影院| 男女之事视频高清在线观看| 在线天堂中文资源库| 91精品三级在线观看| 国产欧美日韩一区二区三| 国产91精品成人一区二区三区| 国产一区在线观看成人免费| 久久人人精品亚洲av| 亚洲成人久久性| 国产99白浆流出| 久久国产精品影院| 丰满饥渴人妻一区二区三| 亚洲欧美精品综合一区二区三区| 国产精品亚洲av一区麻豆| 欧美日韩亚洲综合一区二区三区_| 宅男免费午夜| 日韩欧美一区视频在线观看| 亚洲伊人色综图| 国产欧美日韩一区二区三区在线| 黄网站色视频无遮挡免费观看| 一边摸一边抽搐一进一出视频| 国产免费男女视频| 一本大道久久a久久精品| 曰老女人黄片| 亚洲色图av天堂| 精品高清国产在线一区| 高潮久久久久久久久久久不卡| 窝窝影院91人妻| av电影中文网址| a级毛片黄视频| 亚洲精品一卡2卡三卡4卡5卡| 久久婷婷成人综合色麻豆| 韩国精品一区二区三区| 日本欧美视频一区| 国产主播在线观看一区二区| 在线视频色国产色| 美女午夜性视频免费| 欧美黑人精品巨大| 黄片小视频在线播放| 19禁男女啪啪无遮挡网站| 91av网站免费观看| 精品福利永久在线观看| 欧美成狂野欧美在线观看| tocl精华| 91国产中文字幕| 99热国产这里只有精品6| 国产精品秋霞免费鲁丝片| 久久这里只有精品19| 久久久精品国产亚洲av高清涩受| 老司机午夜福利在线观看视频| 日韩欧美在线二视频| 久久影院123| 国产aⅴ精品一区二区三区波| 一区二区三区国产精品乱码| 亚洲av片天天在线观看| 看黄色毛片网站| 人人澡人人妻人| 久久 成人 亚洲| 91精品国产国语对白视频| 精品久久久精品久久久| 99久久99久久久精品蜜桃| av电影中文网址| 在线观看66精品国产| 嫩草影视91久久| 欧美在线黄色| 两个人免费观看高清视频| 免费av毛片视频| 天堂动漫精品| 99在线人妻在线中文字幕| 成人亚洲精品一区在线观看| 人成视频在线观看免费观看| 真人一进一出gif抽搐免费| 嫁个100分男人电影在线观看| 在线av久久热| 色尼玛亚洲综合影院| 91精品国产国语对白视频| 757午夜福利合集在线观看| 老司机靠b影院| 村上凉子中文字幕在线| 人人妻人人爽人人添夜夜欢视频| 亚洲激情在线av| 天堂动漫精品| 亚洲av第一区精品v没综合| 成人手机av| 久久久国产成人精品二区 | 亚洲成人精品中文字幕电影 | 久久精品成人免费网站| 中国美女看黄片| 日韩一卡2卡3卡4卡2021年| 久久久久国产一级毛片高清牌| 国产成人av激情在线播放| 亚洲专区字幕在线| 成人特级黄色片久久久久久久| 人人妻,人人澡人人爽秒播| 在线观看免费视频网站a站| 亚洲av成人av| 亚洲精品中文字幕一二三四区| 如日韩欧美国产精品一区二区三区| 亚洲成国产人片在线观看| 久久影院123| 欧美成人免费av一区二区三区| 国产成人欧美在线观看| 精品欧美一区二区三区在线| 亚洲精品国产一区二区精华液| 精品国产一区二区久久| 国产精品秋霞免费鲁丝片| 国产欧美日韩一区二区精品| 91麻豆av在线| 99精品在免费线老司机午夜| 757午夜福利合集在线观看| 天天躁夜夜躁狠狠躁躁| 男人舔女人下体高潮全视频| 黄色视频,在线免费观看| 亚洲色图 男人天堂 中文字幕| 性欧美人与动物交配| 极品人妻少妇av视频| av片东京热男人的天堂| 国产精品久久电影中文字幕| 男男h啪啪无遮挡| 国产99白浆流出| 美女午夜性视频免费| 日日干狠狠操夜夜爽| 国产欧美日韩一区二区精品| 午夜久久久在线观看| 久久国产亚洲av麻豆专区| 侵犯人妻中文字幕一二三四区| 18禁观看日本| 国产免费av片在线观看野外av| 国产一区二区三区在线臀色熟女 | 久久精品亚洲熟妇少妇任你| 热re99久久精品国产66热6| 天天躁狠狠躁夜夜躁狠狠躁| 19禁男女啪啪无遮挡网站| 视频在线观看一区二区三区| 日韩大码丰满熟妇| 亚洲人成电影免费在线| 高清av免费在线| 黄色丝袜av网址大全| 热99国产精品久久久久久7| 黑人巨大精品欧美一区二区mp4| 亚洲 国产 在线| 婷婷六月久久综合丁香| 90打野战视频偷拍视频| 老鸭窝网址在线观看| 十分钟在线观看高清视频www| 可以在线观看毛片的网站| 757午夜福利合集在线观看| 亚洲av成人不卡在线观看播放网| 久久中文看片网| 高清黄色对白视频在线免费看| 亚洲人成电影免费在线| 高清在线国产一区| 欧美乱码精品一区二区三区| 国产精品亚洲一级av第二区| 美女午夜性视频免费| 51午夜福利影视在线观看| 99热国产这里只有精品6| 黑人猛操日本美女一级片| 午夜日韩欧美国产| 国产亚洲av高清不卡| 男女下面进入的视频免费午夜 | 国产单亲对白刺激| av国产精品久久久久影院| 久久国产精品影院| 在线视频色国产色| 精品久久久久久久毛片微露脸| 久久性视频一级片| 一级片'在线观看视频| 热re99久久国产66热| 一级黄色大片毛片| 巨乳人妻的诱惑在线观看| 国产91精品成人一区二区三区| 久久99一区二区三区| 三级毛片av免费| 淫妇啪啪啪对白视频| 美国免费a级毛片| av中文乱码字幕在线| 国产精品久久久人人做人人爽| 免费观看人在逋| 热99re8久久精品国产| 欧美久久黑人一区二区| 999久久久精品免费观看国产| 亚洲片人在线观看| 日本a在线网址| 人人妻人人澡人人看| 高清在线国产一区| 又黄又爽又免费观看的视频| 在线视频色国产色| 男女高潮啪啪啪动态图| a级毛片在线看网站| 国产av又大| 自拍欧美九色日韩亚洲蝌蚪91| 脱女人内裤的视频| 老司机在亚洲福利影院| 国产精品影院久久| 老汉色av国产亚洲站长工具| 国产精品久久久久久人妻精品电影| 国产成人一区二区三区免费视频网站| 国产av一区在线观看免费| 又黄又爽又免费观看的视频| 亚洲精品在线观看二区| 人人妻人人爽人人添夜夜欢视频| 亚洲黑人精品在线| 亚洲成国产人片在线观看| 91九色精品人成在线观看| 99久久久亚洲精品蜜臀av| 岛国视频午夜一区免费看| 亚洲人成电影观看| 丰满人妻熟妇乱又伦精品不卡| 久久久久国产精品人妻aⅴ院| 一级a爱片免费观看的视频| 国产精品自产拍在线观看55亚洲| 69精品国产乱码久久久| 国产高清国产精品国产三级| 亚洲专区字幕在线| 亚洲 国产 在线| 国产黄a三级三级三级人| 日韩大码丰满熟妇| 欧洲精品卡2卡3卡4卡5卡区| 女性被躁到高潮视频| 村上凉子中文字幕在线| av有码第一页| 欧美日韩视频精品一区| 伦理电影免费视频| 99国产精品一区二区三区| 婷婷精品国产亚洲av在线| 久久精品影院6| 精品熟女少妇八av免费久了| 国产在线精品亚洲第一网站| 波多野结衣av一区二区av| 精品国产美女av久久久久小说| cao死你这个sao货| 美女扒开内裤让男人捅视频| 成人永久免费在线观看视频| 村上凉子中文字幕在线| 国产av精品麻豆| 久久久国产成人精品二区 | 精品一品国产午夜福利视频| 男女之事视频高清在线观看| 在线观看免费午夜福利视频| 久久热在线av| 丰满饥渴人妻一区二区三| 国产欧美日韩一区二区三| 99riav亚洲国产免费| 久久久久国产一级毛片高清牌| 久久精品国产综合久久久| 国产精品综合久久久久久久免费 | 欧美激情久久久久久爽电影 | 波多野结衣av一区二区av| videosex国产| 国产极品粉嫩免费观看在线| av天堂久久9| 黄网站色视频无遮挡免费观看| a级片在线免费高清观看视频| 一级片免费观看大全| 国产激情久久老熟女| 无限看片的www在线观看| 亚洲精品在线美女| 人人妻人人爽人人添夜夜欢视频| 老鸭窝网址在线观看| 欧美另类亚洲清纯唯美| 久久精品国产综合久久久| 黑人巨大精品欧美一区二区蜜桃| 国产精品爽爽va在线观看网站 | 国产1区2区3区精品| 亚洲成人免费电影在线观看| 久久精品aⅴ一区二区三区四区| 51午夜福利影视在线观看| 欧美中文综合在线视频| 日日爽夜夜爽网站| 精品一区二区三区视频在线观看免费 | bbb黄色大片| 黄色a级毛片大全视频| 亚洲欧美激情在线| 国产精品久久久人人做人人爽| 18美女黄网站色大片免费观看| 欧美另类亚洲清纯唯美| 少妇 在线观看| 久久亚洲精品不卡| 国产精品 国内视频| 黄片播放在线免费| 12—13女人毛片做爰片一| 国产精品免费一区二区三区在线| 在线观看午夜福利视频| 高清在线国产一区| 国产精品亚洲一级av第二区| 人人妻人人爽人人添夜夜欢视频| 久9热在线精品视频| 亚洲 欧美一区二区三区| 一进一出抽搐gif免费好疼 | 国产成人系列免费观看| 国产一区二区在线av高清观看| 两性夫妻黄色片| 亚洲中文日韩欧美视频| 国产精品久久久av美女十八| 在线看a的网站| 美女扒开内裤让男人捅视频| 黄色成人免费大全| 丁香欧美五月| 69av精品久久久久久| xxx96com| 亚洲av成人不卡在线观看播放网| 亚洲专区中文字幕在线| 欧美老熟妇乱子伦牲交| 国产黄色免费在线视频| 制服人妻中文乱码| 老司机靠b影院| 国内毛片毛片毛片毛片毛片| 人妻久久中文字幕网| 成人特级黄色片久久久久久久| 国产91精品成人一区二区三区| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲 欧美一区二区三区| 午夜免费成人在线视频| 亚洲一码二码三码区别大吗| 欧洲精品卡2卡3卡4卡5卡区| 少妇 在线观看| 在线免费观看的www视频| a级毛片在线看网站| 丰满迷人的少妇在线观看| 伊人久久大香线蕉亚洲五| 国产av精品麻豆| 免费久久久久久久精品成人欧美视频| 热re99久久精品国产66热6| 国产午夜精品久久久久久| 真人一进一出gif抽搐免费| 99精品在免费线老司机午夜| 日韩欧美国产一区二区入口| 国产成人精品久久二区二区91| 国产片内射在线| 精品电影一区二区在线| 少妇裸体淫交视频免费看高清 | 青草久久国产| 日韩大尺度精品在线看网址 | 国产精品影院久久| 亚洲九九香蕉| 人妻久久中文字幕网| 国产乱人伦免费视频| 亚洲,欧美精品.| 三级毛片av免费| 曰老女人黄片| 精品一区二区三卡| 午夜精品久久久久久毛片777| 热99国产精品久久久久久7| 99热只有精品国产| 操出白浆在线播放| 久久精品亚洲熟妇少妇任你| 亚洲精品一二三| 久久久精品国产亚洲av高清涩受| 丁香六月欧美| 亚洲精品在线观看二区| 久久久久久久久久久久大奶| 久久人人爽av亚洲精品天堂| 一进一出抽搐gif免费好疼 | 欧美乱码精品一区二区三区| 国产三级在线视频| 99riav亚洲国产免费| 亚洲一区二区三区不卡视频| 亚洲精品一二三| 亚洲免费av在线视频| 麻豆成人av在线观看| 日韩视频一区二区在线观看| 美女 人体艺术 gogo| 老司机深夜福利视频在线观看| 国产免费av片在线观看野外av| 后天国语完整版免费观看| 啦啦啦在线免费观看视频4| 亚洲第一av免费看| 中文欧美无线码| 韩国av一区二区三区四区| 日韩欧美三级三区| 免费不卡黄色视频| 变态另类成人亚洲欧美熟女 | 国产欧美日韩精品亚洲av| 亚洲欧美激情综合另类| 久热爱精品视频在线9| 国产精品av久久久久免费| 99久久综合精品五月天人人| 国产一区二区三区视频了| 午夜亚洲福利在线播放| svipshipincom国产片| 大型黄色视频在线免费观看| 黄频高清免费视频| 午夜免费鲁丝| 在线观看一区二区三区激情| 日韩 欧美 亚洲 中文字幕| 99re在线观看精品视频| 久99久视频精品免费| av免费在线观看网站| 日韩欧美一区二区三区在线观看| 久久久精品国产亚洲av高清涩受| 国产精品98久久久久久宅男小说| 午夜免费观看网址| 国产精品秋霞免费鲁丝片| 久久人人爽av亚洲精品天堂| 变态另类成人亚洲欧美熟女 | 亚洲av日韩精品久久久久久密| 悠悠久久av| 很黄的视频免费| 久久精品91无色码中文字幕| 国产精品一区二区免费欧美| 亚洲伊人色综图| 色精品久久人妻99蜜桃| 色综合婷婷激情| 久久精品影院6| svipshipincom国产片| 午夜免费成人在线视频| 久久久久国内视频| 国产99久久九九免费精品| 日韩欧美国产一区二区入口| a在线观看视频网站| 夜夜躁狠狠躁天天躁| 在线观看免费视频网站a站| 99在线视频只有这里精品首页| 在线观看免费视频日本深夜| 一夜夜www| 午夜91福利影院| 国产一区二区激情短视频| 黄片播放在线免费| 久久久久亚洲av毛片大全| 日韩大码丰满熟妇| 麻豆一二三区av精品| 两个人免费观看高清视频| 黄色片一级片一级黄色片| 精品国产亚洲在线| 1024视频免费在线观看| 国产精品野战在线观看 | 久久草成人影院| 波多野结衣高清无吗| 91国产中文字幕| 老熟妇乱子伦视频在线观看| 国产精品自产拍在线观看55亚洲| 黄色a级毛片大全视频| 黄色女人牲交| 精品久久蜜臀av无| 久久人人97超碰香蕉20202| 国产色视频综合| 水蜜桃什么品种好| 少妇裸体淫交视频免费看高清 | 精品午夜福利视频在线观看一区| 欧美日韩av久久| 天天添夜夜摸| 精品久久久久久久久久免费视频 | 中文字幕人妻丝袜制服| 在线国产一区二区在线| 18禁国产床啪视频网站| 欧美日韩亚洲国产一区二区在线观看| 神马国产精品三级电影在线观看 | 久久香蕉激情| av免费在线观看网站| а√天堂www在线а√下载| 日日摸夜夜添夜夜添小说| 在线免费观看的www视频| 多毛熟女@视频| 涩涩av久久男人的天堂| av福利片在线| 国产一区二区三区视频了| 国产色视频综合| 精品少妇一区二区三区视频日本电影| 电影成人av| 欧美不卡视频在线免费观看 | 国产深夜福利视频在线观看| 久久久久久久精品吃奶| 男男h啪啪无遮挡| 亚洲自拍偷在线| 一区二区三区精品91| 成年版毛片免费区| 日日干狠狠操夜夜爽| 好看av亚洲va欧美ⅴa在| videosex国产| xxxhd国产人妻xxx| 天堂动漫精品| 免费av中文字幕在线| 1024香蕉在线观看| 亚洲avbb在线观看| 无遮挡黄片免费观看| 黄色片一级片一级黄色片| 99久久久亚洲精品蜜臀av| 一个人观看的视频www高清免费观看 | 欧美另类亚洲清纯唯美| 99国产极品粉嫩在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜两性在线视频| 亚洲人成电影观看| 一级黄色大片毛片| 久久精品91无色码中文字幕| 国产精品一区二区精品视频观看| 五月开心婷婷网| 丁香六月欧美| 中文字幕人妻丝袜一区二区| 久久热在线av| 日韩欧美三级三区| 亚洲情色 制服丝袜| 亚洲成人国产一区在线观看| 国产片内射在线| 欧美午夜高清在线| 欧洲精品卡2卡3卡4卡5卡区| 精品人妻1区二区| 亚洲人成伊人成综合网2020| av天堂在线播放| 在线av久久热| 搡老岳熟女国产| 大码成人一级视频| 午夜福利欧美成人| 亚洲男人的天堂狠狠| 国产精品香港三级国产av潘金莲| 可以免费在线观看a视频的电影网站| 亚洲成人精品中文字幕电影 | 国产熟女xx| 日韩欧美免费精品| 9191精品国产免费久久| 午夜福利在线免费观看网站| 国产亚洲精品综合一区在线观看 | av网站在线播放免费| 日韩欧美国产一区二区入口| 视频区图区小说| 日韩欧美国产一区二区入口| 亚洲人成伊人成综合网2020| 国产高清视频在线播放一区| 香蕉国产在线看| 12—13女人毛片做爰片一| 午夜精品在线福利| 激情视频va一区二区三区| 午夜免费激情av| 亚洲五月色婷婷综合| 国产精品偷伦视频观看了| 欧美亚洲日本最大视频资源| 在线观看午夜福利视频| 欧美日韩av久久| 亚洲五月婷婷丁香| 日本黄色视频三级网站网址| 国产精品成人在线| 国产无遮挡羞羞视频在线观看| 一区在线观看完整版| 国产欧美日韩一区二区三区在线| 视频区图区小说| 最好的美女福利视频网| 91麻豆av在线| 午夜a级毛片| 久久午夜综合久久蜜桃| 国产蜜桃级精品一区二区三区| 国产aⅴ精品一区二区三区波| 人人妻人人爽人人添夜夜欢视频| 淫妇啪啪啪对白视频| 国产欧美日韩精品亚洲av| 看黄色毛片网站| 日韩视频一区二区在线观看| 午夜福利影视在线免费观看| 久久久水蜜桃国产精品网| 啦啦啦在线免费观看视频4| 久久精品人人爽人人爽视色| 在线观看舔阴道视频| 满18在线观看网站| 一级a爱视频在线免费观看| 午夜福利在线观看吧| 午夜两性在线视频| 国产色视频综合| 女人被躁到高潮嗷嗷叫费观| 国产成人一区二区三区免费视频网站| 久久人人精品亚洲av| 一级片'在线观看视频| 男人舔女人下体高潮全视频| 国内久久婷婷六月综合欲色啪| 男女床上黄色一级片免费看| 看片在线看免费视频| 久久久精品欧美日韩精品| 欧美国产精品va在线观看不卡| 在线看a的网站| 国产精品亚洲一级av第二区| 久9热在线精品视频| av中文乱码字幕在线| 少妇的丰满在线观看|