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

    基于ARMA與GARCH-M模型對我國豆粕期貨價(jià)格波動(dòng)的分析預(yù)測

    2021-04-27 12:34:52原云霄于惠蘭崔靜
    飼料博覽 2021年2期
    關(guān)鍵詞:豆粕對數(shù)方差

    原云霄,于惠蘭,崔靜

    (1.青島農(nóng)業(yè)大學(xué)海都學(xué)院經(jīng)濟(jì)與管理系,山東煙臺(tái)265200;2.城陽區(qū)人民醫(yī)院,山東青島266109)

    豆粕作為富含高蛋白質(zhì)的油粕品種,其總量的85%用作家禽牲畜飼料。中國是全球進(jìn)口豆類農(nóng)產(chǎn)品最大的國家,2017年中國進(jìn)口大豆9 700萬t,約占全球進(jìn)口總量的64%,來自美國、巴西以及阿根廷的大豆占比分別為34%、53%、6.9%,若由我國自己生產(chǎn)則至少需要4億hm2土地,可知中國豆類農(nóng)產(chǎn)品需求嚴(yán)重依賴進(jìn)口。

    2020年伴隨著中美貿(mào)易戰(zhàn)局勢愈演愈烈,根據(jù)大豆收割周期,二、三季度通過轉(zhuǎn)變豆類進(jìn)口結(jié)構(gòu)勉強(qiáng)由巴西供應(yīng),但到了一、四季度即便降低甚至免去豆類進(jìn)口關(guān)稅也很難挖掘其他國家豆類供應(yīng)潛力,因此豆粕期貨交易成為復(fù)雜形勢中緩解我國豆粕需求的必備工具。

    豆粕期貨在我國農(nóng)產(chǎn)品期貨交易規(guī)模、市場定價(jià)中均占有重要地位,受國際資本投機(jī)、國際形勢變化、自然災(zāi)害等影響,我國豆粕期貨價(jià)格波動(dòng)呈現(xiàn)短期聚集性以及杠桿效應(yīng),因此深入研究豆粕期貨價(jià)格波動(dòng)規(guī)律、預(yù)測其未來價(jià)格走勢不論對國內(nèi)豆粕期貨市場穩(wěn)定還是緩沖國內(nèi)豆粕需求缺口均具有重要意義。

    1 文獻(xiàn)綜述

    梳理以往文獻(xiàn),發(fā)現(xiàn)國內(nèi)外學(xué)者對期貨價(jià)格波動(dòng)做了大量研究且大多采用高頻數(shù)據(jù),A.J.Foster[1]研究國際原油期貨價(jià)格得出,期貨價(jià)格波動(dòng)與交易量呈現(xiàn)明顯的正相關(guān);J.Fleming等[2]研究發(fā)現(xiàn),當(dāng)期恐慌指數(shù)波動(dòng)和前期股指收益存在微弱正相關(guān);P.Giot[3]研究納指波動(dòng)指數(shù)、恐慌指數(shù)與未來市場收益間的聯(lián)系發(fā)現(xiàn),高波動(dòng)率指數(shù)可預(yù)示市場超賣;王秀東等[4]研究發(fā)現(xiàn),大豆期貨存在二階自回歸條件異方差,且金融危機(jī)擴(kuò)大了期貨價(jià)格波動(dòng)率;肖忠意等[5]選用周度數(shù)據(jù)研究國際飼料糧期貨價(jià)對國內(nèi)肉產(chǎn)品沖擊,發(fā)現(xiàn)豬牛羊肉價(jià)格均有異方差性且均對國際大豆、豆粕期貨價(jià)格反應(yīng)顯著;王宏磊等[6]研究中美大豆期貨市場價(jià)格,發(fā)現(xiàn)美國大豆期價(jià)引導(dǎo)中國大豆期價(jià),而全球經(jīng)濟(jì)危機(jī)、石油危機(jī)等將沖擊中美大豆期價(jià)關(guān)系。

    本研究結(jié)合以往學(xué)者經(jīng)驗(yàn),選取高頻日度數(shù)據(jù)研究國內(nèi)豆粕期貨價(jià)格波動(dòng)率異方差規(guī)律并做出短期期價(jià)走勢分析,構(gòu)建ARMA模型預(yù)測豆粕期價(jià)長期走勢,以期為有關(guān)部門制定政策提供數(shù)據(jù)支持。

    2 模型構(gòu)建方法與指標(biāo)選取

    2.1 高頻日度數(shù)據(jù)帶均值廣義自回歸條件異方差GARCH-M模型

    金融資產(chǎn)價(jià)格收益率呈現(xiàn)波動(dòng)聚集性,即大波動(dòng)接著大波動(dòng)、小波動(dòng)接著小波動(dòng)的特征,這代表收益率序列自相關(guān)即存在異方差性,R.F.Engle 1982年提出自回歸條件異方差A(yù)RCH模型,可描述金融資產(chǎn)收益率波動(dòng)性,但鑒于ARCH模型違反參數(shù)非負(fù)數(shù)以及滯后階難以確定的弊端,T.Boller?slev 1986年提出GARCH模型,使條件方差包含殘差前期平方和隨機(jī)誤差項(xiàng)條件方差前期值。本研究構(gòu)建GARCH-M模型,將風(fēng)險(xiǎn)因素(特指條件方差的自然對數(shù))加入收益率均值方程中,首先確定線性回歸方程是否存在異方差,其模型基本形式如下:

    式中:γt為金融資產(chǎn)收益,αk為回歸系數(shù),xkt為影響因素,ut為隨機(jī)誤差項(xiàng)。

    隨機(jī)誤差項(xiàng)條件方差模型如下:

    式中檢驗(yàn)條件方差可通過殘差平方圖或者自回歸條件異方差殘差序列LM檢驗(yàn)來確定,其檢驗(yàn)輔助回歸方程如下:

    本研究將豆粕期貨價(jià)格指數(shù)序列dpjg對數(shù)降低估計(jì)時(shí)的誤差,建立隨機(jī)游走模型:

    在確定存在ARCH效應(yīng)的前提下,構(gòu)建MGARCH模型:

    選取2016年3月16日—2018年6月22日豆粕期貨收盤價(jià)日度高頻數(shù)據(jù),數(shù)據(jù)來源WIND咨詢平臺(tái)。

    2.2 低頻月度數(shù)據(jù)ARMA模型

    對于平穩(wěn)的時(shí)間序列(收益率波動(dòng)),ARMA模型能夠分析相關(guān)經(jīng)濟(jì)變量的變動(dòng)規(guī)律,自回歸移動(dòng)平均模型包括一個(gè)自回歸AR(p)與一個(gè)移動(dòng)平均MA(q),表達(dá)式如下:

    式中:p、q為滯后階數(shù),μt為白噪聲序列。首先檢驗(yàn)序列穩(wěn)定性,其次根據(jù)相關(guān)圖以及不同滯后階數(shù)代數(shù)嘗試的方式獲得ARMA模型的階數(shù)(p,q),根據(jù)識(shí)別后的模型選取適當(dāng)范圍的樣本(留出一定樣本對比)進(jìn)行模型估計(jì)并對模型進(jìn)行診斷,可通過模型擬合序列與原序列與殘差序列白噪聲檢驗(yàn),通過白噪聲檢驗(yàn)最終對預(yù)留樣本進(jìn)行預(yù)測,并可通過動(dòng)態(tài)預(yù)測獲得樣本范圍外的預(yù)測值。

    選取2016年3月16日—2018年6月22日豆粕期貨收盤價(jià)日度高頻數(shù)據(jù),經(jīng)過均值處理轉(zhuǎn)化為月度低頻數(shù)據(jù),數(shù)據(jù)來源WIND咨詢平臺(tái)。

    3 實(shí)證分析

    3.1 GARCH-M模型建立及分析

    3.1.1 ARCH效應(yīng)檢驗(yàn)

    首先對豆粕期貨價(jià)格取對數(shù)后建立隨機(jī)游走模型,見式(4),將數(shù)據(jù)帶入后得到回歸估計(jì)結(jié)果,見式(8)。

    對數(shù)似然值L=1 721.846,AIC=-5.81,SC=-5.79

    對隨機(jī)游走模型(8)殘差進(jìn)行ARCH檢驗(yàn),生成殘差序列折線圖,見圖1。

    圖1 隨機(jī)游走模型殘差序列折線圖

    由圖1可知,2016年第一、二、三季度,大的波動(dòng)后緊接著大的波動(dòng),同樣的情況,2016年第四季度至2018年第一季度,小波動(dòng)跟著小波動(dòng),回歸方程方差表現(xiàn)出波動(dòng)的聚集性,極有可能存在ARCH效應(yīng),因此繪制殘差平方相關(guān)圖做進(jìn)一步判斷,見圖2。

    由圖2可知,自相關(guān)函數(shù)Autocorreltion(AC)1,2階超過置信區(qū)間,Q統(tǒng)計(jì)量(Q-Statistic)極其顯著且概率值非常??;拒絕原序列不存在自相關(guān)的原假設(shè),可認(rèn)為殘差序列存在ARCH效應(yīng)。通過殘差序列拉格朗日殘差檢驗(yàn)LM檢驗(yàn)得到同樣結(jié)果,見圖3。

    圖2 隨機(jī)游走模型殘差平方相關(guān)圖

    圖3 隨機(jī)游走模型殘差序列LM檢驗(yàn)

    3.1.2 GARCH-M模型識(shí)別與估計(jì)

    鑒于回歸方程(8)殘差存在ARCH效應(yīng),通常GARCH對金融資產(chǎn)收益率建模,因此先對豆粕期貨指數(shù)收盤價(jià)格進(jìn)行一階自然對數(shù)差分(通過當(dāng)期/前期后再取自然對數(shù)獲得,ldpsylt=ln(dpjgt/dpjgt-1)),對對數(shù)收益率擬合自回歸AR(p),根據(jù)ldpsylt的相關(guān)圖(見圖4)初步判斷p=1。

    圖4 豆粕期貨指數(shù)收盤價(jià)對數(shù)收益率相關(guān)圖

    由圖4可知,豆粕期貨收益率自相關(guān)與偏相關(guān)函數(shù)除了2階外均在置信區(qū)間范圍內(nèi)(判斷q為低階1),Q統(tǒng)計(jì)量概率值均>10%檢驗(yàn)水平,判斷收益率序列不存在自相關(guān),其均值方程形式:ldpsylt=c+ut,隨后估計(jì)GARCH模型,估計(jì)結(jié)果見圖5。

    通過豆粕期貨收益率均值方程(9)可知豆粕期價(jià)日對數(shù)收益率方差通常>0,但很小,因此整體均值基本在零值,條件方差估計(jì)參數(shù)均>0,從而保證方差為正,P值極小說明均顯著,首先意味著當(dāng)期豆粕期價(jià)波動(dòng)率不僅受前期誤差信息(如外部信息差異或難以計(jì)量忽略)影響,還受自身前期波動(dòng)率影響,其次表明豆粕期價(jià)收益率呈現(xiàn)明顯的聚集性即大的波動(dòng)后緊跟著大的波動(dòng),小的波動(dòng)后持續(xù)小波動(dòng)。政府部門應(yīng)在豆粕出現(xiàn)波動(dòng)苗頭初期就干預(yù)指導(dǎo),特別是在中美貿(mào)易戰(zhàn)背景下(中國對美豆征收農(nóng)產(chǎn)品關(guān)稅),密切關(guān)注巴西、南美等主要大豆出口國供給力出現(xiàn)疲軟時(shí)及早制定應(yīng)對措施;ARCH與GARCH項(xiàng)的估計(jì)系數(shù)和α1+β1=0.88<1,符合GARCH模型參數(shù)約束條件,能夠收斂到無條件方差σ2=α0/(1-α1-β1)=0.000 16,系數(shù)和接近1代表前期沖擊持久影響后期,方差對數(shù)波動(dòng)率每增加1百分點(diǎn)豆粕期貨對數(shù)收益率下降0.002 4百分點(diǎn),即昨天豆粕期價(jià)波動(dòng)率會(huì)導(dǎo)致今天收益率的下跌。這與傳統(tǒng)認(rèn)為風(fēng)險(xiǎn)越大金融資產(chǎn)收益越高相違背。鑒于期貨資產(chǎn)屬于以小博大高風(fēng)險(xiǎn)類金融產(chǎn)品,絕大多數(shù)投資者都難以克服人性弱點(diǎn)頂住高風(fēng)險(xiǎn)壓力持續(xù)買開倉。

    3.1.3 模型分析

    對GARCH-M模型進(jìn)行分析,生成實(shí)際值、擬合值以及殘差值,見圖6。

    由圖6可知,收益率擬合值基本處于零值,這與此前收益率回歸方程判斷結(jié)果一致,盡管殘差序列與實(shí)際值非常相似,仍通過ARCH-LM檢驗(yàn)判斷GARCH-M模型殘差序列是否自相關(guān)即未提取的信息,obs*R-squared=9.071 477且P=0.525 3,接受不存在ARCH效應(yīng)的原假設(shè),ARCH-LM檢驗(yàn)見圖7。

    根據(jù)構(gòu)建的模型對樣本外區(qū)間2018.6.22—2018.8.1做動(dòng)態(tài)預(yù)測,結(jié)果見圖8。

    圖5 GARCH-M估計(jì)結(jié)果

    圖6 GARCH-M模型實(shí)際值、擬合值以及殘差值

    圖8 GARCH-M模型動(dòng)態(tài)預(yù)測結(jié)果

    由圖8可知,豆粕期價(jià)指數(shù)對數(shù)收益率最終圍繞在-0.000 3,呈現(xiàn)平直結(jié)構(gòu)。這是因?yàn)閘dpsylt=-0.021-0.002 4 lg(GARCH)中l(wèi)g(GARCH)為負(fù)數(shù)且較小,與系數(shù)相乘加上常數(shù)項(xiàng)接近于零所致,且均在置信區(qū)間內(nèi),方差項(xiàng)呈現(xiàn)下凹(曲線一階導(dǎo)數(shù)為負(fù)且趨向于零),最終收斂于0.000 17~0.000 18之間的無條件方差,表明隨著時(shí)間推移收益率逐漸回升并最終收斂于-0.000 3。

    3.1.4 模型預(yù)測

    動(dòng)態(tài)預(yù)測結(jié)果見下圖8與表1。

    表1 GARCH-M模型豆粕期價(jià)指數(shù)收盤對數(shù)收益率樣本外2018.6.25—2018.8.1動(dòng)態(tài)預(yù)測結(jié)果

    3.2 ARMA模型建立及分析

    3.2.1 平穩(wěn)性檢驗(yàn)

    將2016年3月—2018年4月豆粕期貨指數(shù)收盤價(jià)對數(shù)收益率月度數(shù)據(jù)記為lsyl序列。為了反映收益率指數(shù)的穩(wěn)定性情況可利用統(tǒng)計(jì)分析軟件EVIEWS8.0提供的自相關(guān)分析圖,若序列的自相關(guān)系數(shù)緩慢趨向于0為不平穩(wěn)的。由圖9可知,自相關(guān)函數(shù)AC均處于置信區(qū)間內(nèi),說明lsyl序列平穩(wěn)。

    對lsyl進(jìn)行ADF單位根檢驗(yàn),見圖10。ADF單位根為-3.498,小于5%顯著性水平下的臨界值-2.976,故該序列是平穩(wěn)的。下一步可以建立ARMA模型。

    圖9 序列對數(shù)收益率lsyl自相關(guān)(AC)、偏相關(guān)(PAC)圖

    圖10 序列ADF單位根檢驗(yàn)結(jié)果

    3.2.2 模型的識(shí)別和定階

    穩(wěn)定的lsyl序列應(yīng)建立ARMA模型,模型的階數(shù)可以借助序列自相關(guān)、偏相關(guān)圖來初步判斷,通常ARMA(p,q)模型自相關(guān)函數(shù)與偏相關(guān)函數(shù)均呈現(xiàn)指數(shù)衰減模式,而lsyl序列符合該模式,余下根據(jù)試錯(cuò)法以系數(shù)顯著性為參照標(biāo)準(zhǔn)確定p=10、q=2。

    3.2.3 參數(shù)估計(jì)和檢驗(yàn)

    選用精確度較高的非線性最小二乘法估計(jì)參數(shù),從被估參數(shù)初始值逐步迭代,使參數(shù)殘差平方和逐漸減小直至給定精度下殘差平方和下降幅度很小為止,ARMA(10,2)模型所有解釋變量的系數(shù)在0.05的顯著性水平下均顯著,這說明豆粕期價(jià)收益率會(huì)對未來第10期的收益率產(chǎn)生影響,且負(fù)向波動(dòng)率每增加1百分點(diǎn),10期后的波動(dòng)率增加0.433百分點(diǎn),同樣前兩期的誤差項(xiàng)對當(dāng)期豆粕期指收益率沖擊為負(fù)向誤差每提升1百分點(diǎn),豆粕期指收益率增加0.888百分點(diǎn),見圖11。

    圖11 ARMA(10,2)模型估計(jì)結(jié)果

    3.2.4 白噪聲檢驗(yàn)

    對模型殘差序列做白噪聲檢驗(yàn),若存在自相關(guān)性,表明仍有信息未提取,采用拉格朗日乘數(shù)檢驗(yàn)法(LM檢驗(yàn)),LM檢驗(yàn)結(jié)果見圖12,統(tǒng)計(jì)量值的P值均>0.1檢驗(yàn)通過,模型殘差序列為白噪聲序列。

    圖12 模型殘差序列LM檢驗(yàn)結(jié)果

    3.2.5 預(yù)留樣本動(dòng)態(tài)預(yù)測

    利用所估計(jì)的ARMA模型對2018年5月、6月的豆粕期指對數(shù)收益率進(jìn)行動(dòng)態(tài)預(yù)測,結(jié)果見表2。由表2可知,預(yù)測誤差非常小,模型預(yù)測效果好。

    表2 2018年5—6月豆粕期指對數(shù)收益率lsyl預(yù)測值、實(shí)際值

    3.2.6 靜態(tài)預(yù)測

    對2016.3—2018.4豆粕期指對數(shù)收益率做樣本內(nèi)靜態(tài)預(yù)測,結(jié)果見圖13。由圖13可知,方差比例為0.012,協(xié)方差為0.985,表明模型的預(yù)測能力好,方差比例較小,說明預(yù)測方差與實(shí)際方差偏離較小,預(yù)測效果較理想??蛇\(yùn)用建立的模型ARMA(10,2)對2018年7月—2018年12月的豆粕期指對數(shù)收益率進(jìn)行預(yù)測。

    3.2.7 樣本外動(dòng)態(tài)預(yù)測

    鑒于模型內(nèi)部預(yù)留樣本的動(dòng)態(tài)預(yù)測與靜態(tài)預(yù)測較好的情形,對2018年7月—2018年12月的豆粕期指對數(shù)收益率進(jìn)行預(yù)測,結(jié)果見表3。預(yù)測數(shù)據(jù)發(fā)現(xiàn),2017年11月—2018年6月,2018年7月—2018年12月的豆粕期指對數(shù)收益率基本處于下降區(qū)間,其中7月收益率下降趨勢與短期日度數(shù)據(jù)預(yù)測結(jié)果一致。

    圖13 ARMA(10,2)模型對總體樣本靜態(tài)預(yù)測圖

    表3 2018年7月—2018年12月豆粕期指對數(shù)收益率樣本外動(dòng)態(tài)預(yù)測結(jié)果

    4 結(jié)論與建議

    第一,國內(nèi)豆粕期貨價(jià)格存在一階ARCH效應(yīng),豆粕期貨價(jià)格指數(shù)波動(dòng)率會(huì)受到前期波動(dòng)率的影響以及前期誤差項(xiàng)的影響,符合金融資產(chǎn)高風(fēng)險(xiǎn)高收益的特性,同時(shí)豆粕期價(jià)收益率波動(dòng)性聚集性爆發(fā)并持續(xù)一段時(shí)間才趨于平穩(wěn)。

    第二,豆粕期指收益率不論短期還是中期均呈現(xiàn)下降趨勢,這說明目前豆粕期貨投資市場以風(fēng)險(xiǎn)回避型投資者為主導(dǎo),對豆粕操作需求處于觀望階段。中美貿(mào)易戰(zhàn)以來豆粕期價(jià)指數(shù)收益率波動(dòng)增強(qiáng),未來豆粕期貨市場不確定性隨著貿(mào)易戰(zhàn)局勢持續(xù)發(fā)酵而不斷加大,這個(gè)時(shí)候貿(mào)然進(jìn)入必然不是明智之舉。

    鑒于以上結(jié)論提出如下建議:

    第一,做好豆粕期價(jià)預(yù)警機(jī)制。豆粕期價(jià)收益率存在GARCH效應(yīng),與其他金融產(chǎn)品一樣呈現(xiàn)波動(dòng)聚集性,應(yīng)及時(shí)關(guān)注豆粕期價(jià)變動(dòng),發(fā)現(xiàn)波動(dòng)苗頭階段及時(shí)采取有效措施,如政府提前儲(chǔ)備一定額度豆粕于其期貨價(jià)格上漲時(shí)及時(shí)投放現(xiàn)貨市場以穩(wěn)住現(xiàn)貨價(jià)格,從而間接制約期貨市場價(jià)格走勢,于豆粕期價(jià)下跌時(shí)及時(shí)購買一定豆粕現(xiàn)貨抑制期價(jià)下行。

    第二,充分發(fā)揮套期保值功能。當(dāng)豆粕供應(yīng)商對未來豆粕期價(jià)指數(shù)走勢感到擔(dān)憂時(shí),應(yīng)在期貨市場賣出一定量豆粕后在期貨市場買進(jìn)同等數(shù)量豆粕期貨,屆時(shí)一旦豆粕現(xiàn)貨升值賠錢,在期貨市場以同樣數(shù)額利潤買入賺回,對于未來豆粕需求商而言,為防止豆粕價(jià)格上漲可先買入、賣出定量豆粕期貨,屆時(shí)即便現(xiàn)貨價(jià)格賠錢,期貨市場可高價(jià)賣出賺回等額利潤,

    第三,調(diào)整我國豆粕進(jìn)口結(jié)構(gòu)與提升豆粕提取深加工水平并舉。美國作為中國最大的豆粕出口方隨著貿(mào)易戰(zhàn)的開打,中國必將轉(zhuǎn)向其他國家收購豆粕,其中南美、巴西將逐步取代美國成為中國進(jìn)口豆粕主要來源國,鑒于中國對豆粕需求量較大,提升豆粕提取工藝特別是一浸豆粕生產(chǎn)工藝水平,同時(shí)提高豆粕轉(zhuǎn)換為家禽牲畜飼料的利用率,從源頭與出口兩方共同努力保證我國豆粕的穩(wěn)定供應(yīng)。

    猜你喜歡
    豆粕對數(shù)方差
    方差怎么算
    含有對數(shù)非線性項(xiàng)Kirchhoff方程多解的存在性
    指數(shù)與對數(shù)
    概率與統(tǒng)計(jì)(2)——離散型隨機(jī)變量的期望與方差
    指數(shù)與對數(shù)
    豆粕:養(yǎng)殖飼料需求回升 國內(nèi)豆粕價(jià)格上漲
    豆粕:貿(mào)易談判持續(xù)進(jìn)行 國內(nèi)豆粕價(jià)格振蕩
    豆粕:貿(mào)易談判再生變數(shù) 國內(nèi)豆粕價(jià)格上漲
    計(jì)算方差用哪個(gè)公式
    對數(shù)簡史
    捣出白浆h1v1| 黑人巨大精品欧美一区二区蜜桃| 欧美大码av| 国产精品久久久久久人妻精品电影 | 亚洲五月色婷婷综合| 黄色a级毛片大全视频| 波多野结衣一区麻豆| 久久久国产精品麻豆| 亚洲欧美日韩另类电影网站| 国产av国产精品国产| 久久精品国产99精品国产亚洲性色 | 咕卡用的链子| 国产欧美日韩综合在线一区二区| 久久久久精品人妻al黑| 一区福利在线观看| 日韩免费av在线播放| 国产激情久久老熟女| 91国产中文字幕| 亚洲九九香蕉| a级片在线免费高清观看视频| 国产真人三级小视频在线观看| 欧美国产精品va在线观看不卡| 正在播放国产对白刺激| 在线观看免费午夜福利视频| 91av网站免费观看| 日韩三级视频一区二区三区| 午夜免费成人在线视频| 免费在线观看视频国产中文字幕亚洲| 一级毛片精品| 青青草视频在线视频观看| 午夜成年电影在线免费观看| 欧美精品一区二区大全| av线在线观看网站| 女同久久另类99精品国产91| 一边摸一边抽搐一进一出视频| 久久ye,这里只有精品| 久久久精品国产亚洲av高清涩受| 精品国产乱码久久久久久男人| 国产精品久久久人人做人人爽| av天堂在线播放| 少妇粗大呻吟视频| 99九九在线精品视频| av网站在线播放免费| 十八禁人妻一区二区| 露出奶头的视频| 亚洲欧美一区二区三区黑人| 欧美成人午夜精品| 午夜福利视频精品| 一级片'在线观看视频| 欧美乱妇无乱码| 国产伦人伦偷精品视频| 欧美激情 高清一区二区三区| 十八禁高潮呻吟视频| 亚洲色图av天堂| 亚洲色图综合在线观看| 欧美另类亚洲清纯唯美| 91老司机精品| 亚洲avbb在线观看| 超色免费av| 亚洲九九香蕉| 中文字幕另类日韩欧美亚洲嫩草| 久久毛片免费看一区二区三区| 午夜精品国产一区二区电影| 久久婷婷成人综合色麻豆| 动漫黄色视频在线观看| 少妇 在线观看| 国产精品亚洲一级av第二区| 国产亚洲精品一区二区www | 男女边摸边吃奶| av网站免费在线观看视频| 一二三四社区在线视频社区8| 国产精品国产av在线观看| 女人高潮潮喷娇喘18禁视频| 免费在线观看影片大全网站| 黄色视频在线播放观看不卡| 热99国产精品久久久久久7| 国产男女内射视频| 国产熟女午夜一区二区三区| 天天添夜夜摸| 男女床上黄色一级片免费看| 男女无遮挡免费网站观看| 中文字幕制服av| 麻豆成人av在线观看| 精品卡一卡二卡四卡免费| 国产精品一区二区在线不卡| 日韩中文字幕视频在线看片| 欧美日韩福利视频一区二区| 午夜福利在线观看吧| 精品国产乱码久久久久久男人| 99在线人妻在线中文字幕 | 久久久久久人人人人人| 亚洲欧美日韩另类电影网站| 免费在线观看视频国产中文字幕亚洲| 国产99久久九九免费精品| 午夜久久久在线观看| 老司机午夜十八禁免费视频| 日本av免费视频播放| 亚洲精华国产精华精| 亚洲九九香蕉| 五月开心婷婷网| av视频免费观看在线观看| 91成人精品电影| 国产亚洲精品第一综合不卡| 国产精品熟女久久久久浪| 成人精品一区二区免费| 亚洲国产欧美网| 夜夜骑夜夜射夜夜干| 日韩免费高清中文字幕av| 免费日韩欧美在线观看| 色播在线永久视频| av天堂久久9| e午夜精品久久久久久久| 桃花免费在线播放| 国产97色在线日韩免费| 午夜精品久久久久久毛片777| 国产成人影院久久av| 中文字幕色久视频| 国产片内射在线| 91麻豆av在线| 成人国产av品久久久| 国产一区二区在线观看av| 后天国语完整版免费观看| 精品国产一区二区久久| cao死你这个sao货| 国产亚洲欧美在线一区二区| 91精品三级在线观看| 国产不卡av网站在线观看| 日韩视频在线欧美| 精品亚洲成国产av| 亚洲精品在线观看二区| 国产高清视频在线播放一区| 精品一区二区三区视频在线观看免费 | 一级毛片精品| 国产99久久九九免费精品| 一边摸一边抽搐一进一小说 | av视频免费观看在线观看| 精品国产一区二区三区四区第35| 亚洲少妇的诱惑av| 日韩欧美国产一区二区入口| 每晚都被弄得嗷嗷叫到高潮| 国产一区二区三区视频了| 国产亚洲av高清不卡| 精品少妇一区二区三区视频日本电影| 热re99久久精品国产66热6| 亚洲熟女毛片儿| 99久久人妻综合| 天天躁狠狠躁夜夜躁狠狠躁| videos熟女内射| 国产主播在线观看一区二区| 高清视频免费观看一区二区| 视频区欧美日本亚洲| 18禁美女被吸乳视频| 最近最新中文字幕大全免费视频| 大陆偷拍与自拍| 国产欧美亚洲国产| 国产精品麻豆人妻色哟哟久久| 精品国产一区二区三区四区第35| 欧美激情久久久久久爽电影 | 男女高潮啪啪啪动态图| 欧美精品av麻豆av| 91老司机精品| 高清欧美精品videossex| av电影中文网址| 一本久久精品| 天天躁日日躁夜夜躁夜夜| 黑人操中国人逼视频| 午夜91福利影院| 色播在线永久视频| 欧美乱码精品一区二区三区| 精品卡一卡二卡四卡免费| 9191精品国产免费久久| 色播在线永久视频| 亚洲欧美一区二区三区久久| 欧美日韩国产mv在线观看视频| 亚洲av成人一区二区三| 天天躁日日躁夜夜躁夜夜| 咕卡用的链子| 久久精品亚洲熟妇少妇任你| 老司机亚洲免费影院| 女性被躁到高潮视频| 丰满饥渴人妻一区二区三| 中文亚洲av片在线观看爽 | 9热在线视频观看99| 国产一区二区三区视频了| 国产日韩欧美亚洲二区| 国产免费av片在线观看野外av| 久久久国产欧美日韩av| 亚洲国产毛片av蜜桃av| 午夜福利乱码中文字幕| 亚洲精品国产精品久久久不卡| 久久人妻熟女aⅴ| 久久99一区二区三区| 1024视频免费在线观看| 最黄视频免费看| 无人区码免费观看不卡 | 亚洲精品中文字幕在线视频| 亚洲美女黄片视频| av视频免费观看在线观看| 欧美老熟妇乱子伦牲交| 99re6热这里在线精品视频| 五月开心婷婷网| 中文字幕另类日韩欧美亚洲嫩草| 狠狠狠狠99中文字幕| 欧美在线一区亚洲| 色精品久久人妻99蜜桃| 成人手机av| av福利片在线| 最新在线观看一区二区三区| 国产亚洲欧美精品永久| 一级毛片电影观看| 精品国产乱子伦一区二区三区| 99国产精品99久久久久| 曰老女人黄片| 久热爱精品视频在线9| 亚洲成人免费电影在线观看| 老熟女久久久| 在线观看免费视频网站a站| 亚洲精品国产精品久久久不卡| 久久亚洲真实| 亚洲国产av新网站| 亚洲成国产人片在线观看| 黑人猛操日本美女一级片| 久久热在线av| 女性被躁到高潮视频| 久久人人97超碰香蕉20202| av天堂在线播放| 精品乱码久久久久久99久播| 中文字幕人妻丝袜制服| 亚洲国产欧美一区二区综合| bbb黄色大片| www.精华液| 精品国产国语对白av| 麻豆成人av在线观看| 悠悠久久av| 麻豆av在线久日| 嫁个100分男人电影在线观看| 精品国产乱码久久久久久小说| 桃花免费在线播放| 久久久久精品国产欧美久久久| 成人18禁高潮啪啪吃奶动态图| 亚洲精品av麻豆狂野| 免费日韩欧美在线观看| 成人特级黄色片久久久久久久 | 丁香六月天网| 国产精品久久久人人做人人爽| 欧美精品一区二区大全| 老司机午夜十八禁免费视频| 久久久国产一区二区| 精品国产乱子伦一区二区三区| 精品久久久精品久久久| 国产精品99久久99久久久不卡| 12—13女人毛片做爰片一| 午夜福利视频在线观看免费| videosex国产| 久久天躁狠狠躁夜夜2o2o| 99精品久久久久人妻精品| 亚洲综合色网址| 一本久久精品| 欧美+亚洲+日韩+国产| www.999成人在线观看| 一边摸一边抽搐一进一小说 | 满18在线观看网站| 国产av又大| 成人影院久久| 国产日韩欧美视频二区| 欧美成人午夜精品| 亚洲男人天堂网一区| 国产成人一区二区三区免费视频网站| 国产亚洲精品一区二区www | 久久国产精品男人的天堂亚洲| 99国产精品一区二区蜜桃av | 亚洲一卡2卡3卡4卡5卡精品中文| 法律面前人人平等表现在哪些方面| 黑人操中国人逼视频| 免费看a级黄色片| 欧美老熟妇乱子伦牲交| 中亚洲国语对白在线视频| 中文字幕色久视频| 搡老熟女国产l中国老女人| 一级a爱视频在线免费观看| 91成人精品电影| 国产男女超爽视频在线观看| 99香蕉大伊视频| 淫妇啪啪啪对白视频| 国产精品国产高清国产av | 777久久人妻少妇嫩草av网站| 亚洲 国产 在线| 久久久欧美国产精品| 久久中文字幕人妻熟女| 成年人黄色毛片网站| 少妇的丰满在线观看| 天堂8中文在线网| 在线十欧美十亚洲十日本专区| 悠悠久久av| 午夜福利一区二区在线看| 大陆偷拍与自拍| 后天国语完整版免费观看| 免费久久久久久久精品成人欧美视频| 美女福利国产在线| 国产亚洲av高清不卡| 国产成人系列免费观看| 成年人免费黄色播放视频| 亚洲国产欧美一区二区综合| av片东京热男人的天堂| 黄色a级毛片大全视频| 黄色视频在线播放观看不卡| 亚洲少妇的诱惑av| 亚洲精品美女久久久久99蜜臀| 亚洲黑人精品在线| 日本vs欧美在线观看视频| 青青草视频在线视频观看| 下体分泌物呈黄色| 精品亚洲成a人片在线观看| 十八禁网站网址无遮挡| 一级,二级,三级黄色视频| 首页视频小说图片口味搜索| 国精品久久久久久国模美| 麻豆乱淫一区二区| 久久天堂一区二区三区四区| 亚洲黑人精品在线| 人妻久久中文字幕网| 在线观看舔阴道视频| 亚洲av成人不卡在线观看播放网| 99在线人妻在线中文字幕 | 视频在线观看一区二区三区| 欧美中文综合在线视频| 成年版毛片免费区| 女人被躁到高潮嗷嗷叫费观| 国产免费av片在线观看野外av| 亚洲欧美一区二区三区黑人| 亚洲国产欧美在线一区| 国产精品熟女久久久久浪| 巨乳人妻的诱惑在线观看| 丝袜美腿诱惑在线| 成人国语在线视频| 午夜91福利影院| 波多野结衣av一区二区av| 老司机影院毛片| 亚洲成国产人片在线观看| 一本大道久久a久久精品| 国产精品麻豆人妻色哟哟久久| 黄频高清免费视频| 色婷婷av一区二区三区视频| 高潮久久久久久久久久久不卡| 侵犯人妻中文字幕一二三四区| 国产在线观看jvid| 精品国产乱码久久久久久男人| 欧美精品一区二区免费开放| 色婷婷av一区二区三区视频| 国产aⅴ精品一区二区三区波| 下体分泌物呈黄色| 久久国产精品男人的天堂亚洲| 精品国产一区二区三区久久久樱花| 一进一出抽搐动态| av国产精品久久久久影院| 动漫黄色视频在线观看| 亚洲视频免费观看视频| 久久免费观看电影| 欧美另类亚洲清纯唯美| 欧美国产精品va在线观看不卡| 91成人精品电影| 亚洲九九香蕉| 久久久久精品国产欧美久久久| 嫩草影视91久久| 日本a在线网址| 桃花免费在线播放| 亚洲国产欧美网| 日韩欧美国产一区二区入口| 国产精品九九99| 国产欧美日韩综合在线一区二区| 热re99久久国产66热| 亚洲欧美色中文字幕在线| 极品教师在线免费播放| 极品人妻少妇av视频| 高清欧美精品videossex| 精品一品国产午夜福利视频| 这个男人来自地球电影免费观看| 亚洲情色 制服丝袜| 久久久久久久国产电影| 国产av精品麻豆| 老司机午夜福利在线观看视频 | 久久久久网色| 国产野战对白在线观看| 大型av网站在线播放| 国产高清国产精品国产三级| 欧美 日韩 精品 国产| 国产精品九九99| 亚洲免费av在线视频| 久久久久网色| 黄片小视频在线播放| tocl精华| 国产亚洲精品一区二区www | 999久久久精品免费观看国产| 久久久国产精品麻豆| 一本色道久久久久久精品综合| 岛国在线观看网站| 成年动漫av网址| 一边摸一边抽搐一进一小说 | 人妻久久中文字幕网| 成年版毛片免费区| 欧美午夜高清在线| 香蕉国产在线看| 中文字幕av电影在线播放| 亚洲精品在线观看二区| 国产有黄有色有爽视频| 91麻豆av在线| 一级a爱视频在线免费观看| 亚洲成人国产一区在线观看| 日本av免费视频播放| 国产av国产精品国产| 亚洲精品在线观看二区| 中文字幕人妻丝袜制服| 亚洲精品国产一区二区精华液| 国产精品98久久久久久宅男小说| 国产一卡二卡三卡精品| 欧美精品一区二区免费开放| 啦啦啦视频在线资源免费观看| 99国产极品粉嫩在线观看| 18禁美女被吸乳视频| 亚洲欧洲日产国产| 一边摸一边做爽爽视频免费| 国产一区二区三区在线臀色熟女 | 精品国产超薄肉色丝袜足j| 精品少妇黑人巨大在线播放| 性少妇av在线| 麻豆成人av在线观看| 成人亚洲精品一区在线观看| 大片电影免费在线观看免费| 中文字幕制服av| cao死你这个sao货| 婷婷丁香在线五月| 精品一区二区三区视频在线观看免费 | 啦啦啦免费观看视频1| 免费看十八禁软件| 午夜福利视频精品| 国产不卡一卡二| 美女主播在线视频| 欧美日韩国产mv在线观看视频| 久久精品熟女亚洲av麻豆精品| 免费观看a级毛片全部| 免费日韩欧美在线观看| 男女床上黄色一级片免费看| 在线av久久热| 激情视频va一区二区三区| 色94色欧美一区二区| 欧美成人免费av一区二区三区 | 亚洲综合色网址| 亚洲欧美激情在线| 日韩熟女老妇一区二区性免费视频| e午夜精品久久久久久久| 成年版毛片免费区| 亚洲精品美女久久久久99蜜臀| 免费不卡黄色视频| 国产精品一区二区精品视频观看| videosex国产| 亚洲成国产人片在线观看| 黑人巨大精品欧美一区二区mp4| 捣出白浆h1v1| 精品久久蜜臀av无| 丝袜在线中文字幕| 交换朋友夫妻互换小说| 一级毛片女人18水好多| 亚洲成人免费av在线播放| 香蕉国产在线看| 热99国产精品久久久久久7| 老汉色av国产亚洲站长工具| av网站在线播放免费| 18禁观看日本| 免费看十八禁软件| 欧美激情 高清一区二区三区| av又黄又爽大尺度在线免费看| 18禁观看日本| a级毛片在线看网站| 精品高清国产在线一区| 亚洲第一青青草原| 婷婷丁香在线五月| 在线亚洲精品国产二区图片欧美| 国产精品 国内视频| 精品久久久精品久久久| 人人妻人人爽人人添夜夜欢视频| 日韩大片免费观看网站| 不卡一级毛片| 国产成人精品久久二区二区91| 99国产精品99久久久久| 两个人免费观看高清视频| 久久久精品国产亚洲av高清涩受| 亚洲中文字幕日韩| 亚洲欧美一区二区三区黑人| 窝窝影院91人妻| 亚洲男人天堂网一区| svipshipincom国产片| 精品国产一区二区三区四区第35| 久久久久久亚洲精品国产蜜桃av| 精品国产亚洲在线| 交换朋友夫妻互换小说| 亚洲国产欧美网| 精品乱码久久久久久99久播| 欧美黄色淫秽网站| 热re99久久精品国产66热6| 性色av乱码一区二区三区2| 久久青草综合色| 久久亚洲精品不卡| 亚洲人成伊人成综合网2020| 国产日韩欧美视频二区| 免费观看人在逋| 一级毛片女人18水好多| 日韩精品免费视频一区二区三区| 国产国语露脸激情在线看| 又紧又爽又黄一区二区| 欧美激情高清一区二区三区| 亚洲精品久久午夜乱码| 在线观看免费视频网站a站| 国产成人一区二区三区免费视频网站| 一级黄色大片毛片| 免费女性裸体啪啪无遮挡网站| 亚洲欧美日韩高清在线视频 | 午夜精品久久久久久毛片777| 国产精品久久久久久精品电影小说| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲精品在线观看二区| 国产精品成人在线| 美女国产高潮福利片在线看| 日韩熟女老妇一区二区性免费视频| 国产高清国产精品国产三级| 亚洲欧美色中文字幕在线| 日韩有码中文字幕| 亚洲中文字幕日韩| 天天操日日干夜夜撸| 欧美 亚洲 国产 日韩一| 黑丝袜美女国产一区| 人人妻人人添人人爽欧美一区卜| 男女无遮挡免费网站观看| 亚洲国产av新网站| 久久久久久免费高清国产稀缺| 久久天躁狠狠躁夜夜2o2o| 亚洲一区中文字幕在线| 国产有黄有色有爽视频| 精品国产国语对白av| 大香蕉久久网| 黄色怎么调成土黄色| 久久天堂一区二区三区四区| svipshipincom国产片| 欧美在线黄色| 亚洲成a人片在线一区二区| 三上悠亚av全集在线观看| 18在线观看网站| 午夜激情av网站| cao死你这个sao货| 亚洲va日本ⅴa欧美va伊人久久| 亚洲人成77777在线视频| 精品一区二区三区av网在线观看 | 免费在线观看完整版高清| 伦理电影免费视频| h视频一区二区三区| 亚洲精品乱久久久久久| 欧美av亚洲av综合av国产av| 国产成人av教育| 亚洲全国av大片| 别揉我奶头~嗯~啊~动态视频| 精品国产超薄肉色丝袜足j| 嫁个100分男人电影在线观看| 亚洲精品一卡2卡三卡4卡5卡| 亚洲七黄色美女视频| 午夜日韩欧美国产| 女人久久www免费人成看片| 亚洲精品久久午夜乱码| a级片在线免费高清观看视频| 亚洲avbb在线观看| 免费少妇av软件| 久久久久久亚洲精品国产蜜桃av| 老熟女久久久| 99久久国产精品久久久| 亚洲国产成人一精品久久久| 国产黄频视频在线观看| 女同久久另类99精品国产91| 国产精品 国内视频| a级毛片黄视频| 亚洲一码二码三码区别大吗| 免费在线观看视频国产中文字幕亚洲| 黑丝袜美女国产一区| 王馨瑶露胸无遮挡在线观看| 五月开心婷婷网| 真人做人爱边吃奶动态| 19禁男女啪啪无遮挡网站| 亚洲精品久久午夜乱码| 精品少妇黑人巨大在线播放| 国产精品一区二区免费欧美| 99国产精品免费福利视频| xxxhd国产人妻xxx| 亚洲精品美女久久av网站| 黄色怎么调成土黄色| 一二三四社区在线视频社区8| 色视频在线一区二区三区| 老汉色av国产亚洲站长工具| av国产精品久久久久影院| 欧美激情极品国产一区二区三区| 在线观看免费高清a一片| av有码第一页| 99久久99久久久精品蜜桃| 俄罗斯特黄特色一大片| 两性午夜刺激爽爽歪歪视频在线观看 | 999久久久精品免费观看国产| 欧美乱妇无乱码| 亚洲精品中文字幕在线视频| 日本av免费视频播放| 免费在线观看影片大全网站| 波多野结衣一区麻豆| 在线亚洲精品国产二区图片欧美| 99热国产这里只有精品6| www.自偷自拍.com| 久久这里只有精品19| 精品亚洲成国产av| 少妇的丰满在线观看| av在线播放免费不卡| www.熟女人妻精品国产|