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

    基于ARIMA模型的余額寶收益率的預(yù)測(cè)

    2017-09-03 10:10:04羅裕凡劉小蓉朱淑雅趙偉帆李敏
    山東科學(xué) 2017年4期
    關(guān)鍵詞:余額差分殘差

    羅裕凡,劉小蓉,朱淑雅,趙偉帆,李敏

    (山東師范大學(xué)數(shù)學(xué)與統(tǒng)計(jì)學(xué)院,山東 濟(jì)南 250014)

    基于ARIMA模型的余額寶收益率的預(yù)測(cè)

    羅裕凡,劉小蓉,朱淑雅,趙偉帆,李敏*

    (山東師范大學(xué)數(shù)學(xué)與統(tǒng)計(jì)學(xué)院,山東 濟(jì)南 250014)

    選取2013年5月30日—2016年7月27日期間的余額寶收益率數(shù)據(jù),進(jìn)行數(shù)據(jù)處理并建立ARIMA模型,從而預(yù)測(cè)出余額寶在2016年7月28日—2016年8月21日期間的收益率。通過實(shí)際結(jié)果與預(yù)測(cè)結(jié)果之間的對(duì)比,確定了在允許誤差范圍內(nèi)該模型的有效性。

    余額寶;收益率;ARIMA模型

    2013年6月,阿里巴巴集團(tuán)和天弘基金合作在支付寶平臺(tái)上正式推出了余額寶,由于其操作的簡(jiǎn)便性和相對(duì)較高的利潤(rùn)而受到很多人的青睞。目前余額寶使用人數(shù)已達(dá)到2.6億,并成為我國(guó)市場(chǎng)上規(guī)模最大的貨幣基金。余額寶是互聯(lián)網(wǎng)的快速發(fā)展帶來的中國(guó)利率市場(chǎng)化的產(chǎn)物,因而受到了很多專家學(xué)者的廣泛關(guān)注。通過分析余額寶收益率的變化情況,并做出準(zhǔn)確合理的預(yù)測(cè),將有利于市場(chǎng)參與方做出自己的決策。王瑩[1]對(duì)余額寶的流動(dòng)性、收益性及風(fēng)險(xiǎn)性進(jìn)行綜合分析,并對(duì)投資者和余額寶自身提出相應(yīng)建議。相對(duì)于銀行的活期存款利息,余額寶的收益率比較高,但是,張大海[2]認(rèn)為將這兩者作比較是不恰當(dāng)?shù)?,并分析說真正和余額寶業(yè)務(wù)相似的是銀行發(fā)行的短期貨幣基金以及開放式無固定期限理財(cái)產(chǎn)品。李恒等[3]從監(jiān)管層面、商業(yè)銀行、產(chǎn)品自身和客戶對(duì)余額寶業(yè)務(wù)的認(rèn)知等多個(gè)方面對(duì)余額寶風(fēng)險(xiǎn)性作了研究。以上文獻(xiàn)都在定性方面對(duì)余額寶作了充分的分析,但是定量方面的研究并未涉及。另一方面,王雪飛等[4-9]基于自四歸積分滑動(dòng)平均模型(auto-regressive integrated moving average model,ARIMA)[10]對(duì)中國(guó)鋼材市場(chǎng)價(jià)格、土建專業(yè)群人才需求、短期風(fēng)速、商品住宅銷售價(jià)格、人口與醫(yī)療需求及中國(guó)人均糧食等方面做了有效的預(yù)測(cè)。本文將通過時(shí)間序列的方法,建構(gòu)ARIMA模型,對(duì)未來短期內(nèi)的余額寶收益率走勢(shì)進(jìn)行預(yù)測(cè)。為了避免因余額寶收益率有時(shí)變化較大而帶來的誤差,本文選取每天的余額寶七日年化收益率作為研究對(duì)象,給出一個(gè)合適的預(yù)測(cè)模型來進(jìn)行預(yù)測(cè)。

    1 ARIMA模型相關(guān)介紹

    ARIMA模型是一種常用的隨機(jī)時(shí)間序列分析模型,這種模型是Box和Jenkins提出的,所以又被稱為B-J模型[10]。這是一種精度相對(duì)較高的短期的時(shí)間序列預(yù)測(cè)方法,其基本思想表現(xiàn)為,即便構(gòu)成時(shí)間序列t的任意一組隨機(jī)變量存在具有不確定性的單個(gè)序列值,但是縱觀整個(gè)時(shí)間序列,其變化卻具有某種規(guī)律。根據(jù)這種思想,我們可以使用相應(yīng)的數(shù)學(xué)模型近似地進(jìn)行描述,進(jìn)而建立數(shù)學(xué)模型,得到最小方差意義下的最優(yōu)預(yù)測(cè)。根據(jù)原序列是否平穩(wěn)以及回歸中所含部分的不同,可以將ARIMA模型分成自回歸模型(auto-regressive model,AR)、移動(dòng)平均模型(moving average model,MA)和自回歸移動(dòng)平均模型(auto-regressive moving average model,ARMA)。MARIMA(p,d,q)稱為差分自回歸移動(dòng)平均模型,其中AR代表自回歸,p代表自回歸項(xiàng);MA代表移動(dòng)平均,q代表移動(dòng)平均項(xiàng)數(shù);d代表時(shí)間序列成為平穩(wěn)時(shí)所做的差分次數(shù)。

    如果時(shí)間序列yt是其前期值和隨機(jī)項(xiàng)組成的線性函數(shù),即可表示為:

    yt=φ1yt-1+φ2yt-2+…+φpyt-p+ut,

    (1)

    令Bk為k步滯后算子,即Bkyt=yt-k,則(1)式可表示成為:

    yt=φ1Byt+φ2B2yt+…+φpBpyt+ut,

    表示滯后多項(xiàng)式,則上述模型可簡(jiǎn)化為:

    如果時(shí)間序列yt是由其當(dāng)期隨機(jī)誤差項(xiàng)和前期隨機(jī)誤差項(xiàng)組成的線性函數(shù),即可表示為:

    yt=ut-ξ1ut-1-ξ2ut-2-…-ξqut-q,

    (2)

    其中,ξ1,ξ2,…,ξq為實(shí)數(shù),稱為移動(dòng)平均系數(shù),則yt是一個(gè)移動(dòng)平均序列,相應(yīng)的(2)式就是一個(gè)q階的移動(dòng)平均模型,將其記作MMA(q)。

    則(2)式可以簡(jiǎn)化為:

    yt=ξ(B)ut。

    為了使AR過程和MA過程能相互轉(zhuǎn)化,即為可逆過程,就必須保證移動(dòng)平均過程無條件平穩(wěn)。因此要求滯后多項(xiàng)式ξ(B)的根都在單位圓外。

    如果時(shí)間序列yt是由其當(dāng)期隨機(jī)誤差項(xiàng),前期隨機(jī)誤差項(xiàng)和前期值組成的線性函數(shù),即可表示為:

    yt=φ1yt-1+φ2yt-2+…+φpyt-p+ut-ξ1ut-1-ξ2ut-2-…-ξqut-q,

    (3)

    則該yt是一個(gè)自回歸序列,相應(yīng)地,(3)式就是一個(gè)(p,q)階的自回歸移動(dòng)平均模型,將其記作MARMA(p,q)。其中,φ1,φ2,…,φp是自回歸系數(shù),ξ1,ξ2,…,ξq是移動(dòng)平均系數(shù)。

    引入滯后算子B,(3)式可以簡(jiǎn)化為:

    2 實(shí)證與預(yù)測(cè)

    2.1 ARIMA模型預(yù)測(cè)的主要步驟

    (1) 對(duì)序列的平穩(wěn)性加以判別。具體是根據(jù)時(shí)間序列的散點(diǎn)圖、自相關(guān)函數(shù)以及偏自相關(guān)函數(shù)圖,對(duì)序列的方差、趨勢(shì)進(jìn)行檢驗(yàn)。如果時(shí)間序列是平穩(wěn)的,直接建立ARIMA模型;否則,必須通過適當(dāng)處理使其滿足平穩(wěn)性要求,具體過程是對(duì)于不平穩(wěn)的數(shù)據(jù)序列,因其具有一定程度的上升或下降的趨勢(shì),那么數(shù)據(jù)需要進(jìn)行差分處理,要保證處理之后的數(shù)據(jù)自相關(guān)函數(shù)值和偏相關(guān)函數(shù)值不顯著為0。

    (2)建立相應(yīng)的模型并進(jìn)行參數(shù)估計(jì),避免高階的MA模型或者包含高階移動(dòng)平均項(xiàng)的ARMA。

    (3) 進(jìn)行模型檢驗(yàn),檢驗(yàn)殘差序列是否為白噪聲。

    (4)進(jìn)行預(yù)測(cè)分析[11]。

    2.2 數(shù)據(jù)

    本文數(shù)據(jù)來源于天天基金網(wǎng)中的天弘余額寶貨幣終端,選取對(duì)象是2013年5月30日—2016年7月27日期間的余額寶七日年化收益率。我們使用ARIMA模型給出2016年7月28日—2016年8月21日這25天的余額寶七日年化收益率。

    2.3 建模過程

    2.3.1 一階差分單位根檢驗(yàn)

    對(duì)序列進(jìn)行一階差分之后,可得到序列圖(圖1),差分后的結(jié)果顯示圖像在零線上下波動(dòng),即模型的擬合值與實(shí)際值的變動(dòng)較小,對(duì)此,我們可假定序列已經(jīng)達(dá)到了平穩(wěn),接下來,為了探究其是否真的達(dá)到平穩(wěn),還需要對(duì)其進(jìn)行單位根檢驗(yàn)。

    圖1 一階差分后序列圖Fig.1 The sequence diagram after the first order difference

    通過單位根檢驗(yàn),得出結(jié)果如表1所示,其中原假設(shè)一階差分有一個(gè)單位根,外生變量是常數(shù),并且滯后長(zhǎng)度是13,該值自動(dòng)基于最大滯后長(zhǎng)度是22的SIC準(zhǔn)則。由表1可得,一階差分之后的序列t檢驗(yàn)的統(tǒng)計(jì)量明顯小于各個(gè)顯著性水平提供的臨界值,所以可以得出,經(jīng)過一階差分之后的時(shí)間序列是平穩(wěn)序列[12]。

    表1 一階差分后序列單位根檢驗(yàn)結(jié)果

    2.3.2 差分序列的相關(guān)性檢驗(yàn)

    因?yàn)榧冸S機(jī)序列不具有建模的必要性,所以為了避免一階差分后的序列呈現(xiàn)純隨機(jī)序列的情況,我們需要對(duì)差分后的序列做相關(guān)性檢驗(yàn),根據(jù)差分后的自相關(guān)系數(shù)來判定其是否為純隨機(jī)序列。進(jìn)行一階差分之后,得出的自相關(guān)函數(shù),偏自相關(guān)函數(shù)以及相關(guān)統(tǒng)計(jì)量的情況如表2所示。

    表2 自相關(guān)函數(shù)和偏自相關(guān)函數(shù)結(jié)果圖

    從表2中我們可以發(fā)現(xiàn),自相關(guān)函數(shù)和偏自相關(guān)函數(shù)都是拖尾的。綜上所述,該時(shí)間序列不是純隨機(jī)序列,因此有建模必要。

    2.3.3 模型最優(yōu)判定

    表3 模型定階原則

    我們通過對(duì)2016年7月27日之前的數(shù)據(jù)進(jìn)行一次差分之后序列的自相關(guān)函數(shù)和偏自相關(guān)函數(shù)的分析,從一階差分的偏自相關(guān)函數(shù)中可以觀察到,在k=3以后的值都在置信區(qū)間內(nèi),偏自相關(guān)函數(shù)很快趨于0,偏自相關(guān)函數(shù)呈拖尾狀態(tài),因此可以確定p=3;從一階差分自相關(guān)函數(shù)可以觀察到,在k=5以后的值都在置信區(qū)間內(nèi),自相關(guān)函數(shù)也很快趨于0,自相關(guān)函數(shù)呈拖尾狀態(tài),因此可以初步確定q取5。在此基礎(chǔ)上我們建立MARIMA(3,1,5)模型,其中關(guān)于MARIMA(3,1,5)模型的相關(guān)參數(shù)和指標(biāo)分別見表4和表5。

    表4 MARIMA(3,1,5)模型的相關(guān)參數(shù)

    表5 MARIMA(3,1,5)模型的指標(biāo)

    2.3.4 模型檢驗(yàn)

    進(jìn)行參數(shù)估計(jì)后,模型的殘差序列還需要進(jìn)行白噪聲檢驗(yàn),主要是檢驗(yàn)其殘差序列是否屬于白噪聲序列[4]。如果殘差序列不屬于白噪聲序列,則還需對(duì)模型進(jìn)行進(jìn)一步的改進(jìn)。我們對(duì)MARIMA(3,1,5)模型殘差的自相關(guān)函數(shù)和偏自相關(guān)函數(shù)進(jìn)行置信區(qū)間檢驗(yàn),LM檢驗(yàn)結(jié)果如圖2所示。

    由圖2可以看出,殘差的自相關(guān)和偏自相關(guān)函數(shù)都處于置信區(qū)間內(nèi),這說明殘差通過了白噪聲檢驗(yàn),可以用于預(yù)測(cè)。

    2.3.5 模型預(yù)測(cè)

    通過上述一系列的分析以及文獻(xiàn)[14]的理論知識(shí),我們可以確定最佳擬合預(yù)測(cè)模型就是MARIMA(3,1,5),圖3是預(yù)測(cè)結(jié)果與實(shí)際結(jié)果,即預(yù)測(cè)和實(shí)際收益率的對(duì)比圖。

    圖2 殘差的LM檢驗(yàn)Fig.2 LM test of the residual errors

    圖3 預(yù)測(cè)結(jié)果與實(shí)際結(jié)果的對(duì)比圖 Fig.3 The comparison of the predicted results with the actual results

    在表6中,我們給出了2016年7月28日—2016年8月21日期間每天的實(shí)際值、預(yù)測(cè)值以及其相對(duì)誤差,最大的相對(duì)誤差為1.146%,因此認(rèn)為該模型是有效的。

    表6 余額寶收益率預(yù)測(cè)結(jié)果與實(shí)際結(jié)果的對(duì)比圖

    續(xù)表6

    日期實(shí)際值預(yù)測(cè)值相對(duì)誤差/%2016?08?182.3672.3820.6342016?08?192.3632.3810.7622016?08?202.3602.3810.8902016?08?212.3592.3800.890

    3 結(jié)論

    本文使用的ARIMA模型,能近似地描述隨時(shí)間變化的序列,基本避免了由于其他經(jīng)濟(jì)原因帶來的誤差,因而本文的預(yù)測(cè)對(duì)象,即余額寶的收益率所組成的序列可以被有效地識(shí)別出來,從而建立相應(yīng)的數(shù)學(xué)模型,通過研究時(shí)間序列的現(xiàn)在值和過去值,對(duì)未來值進(jìn)行預(yù)測(cè)。

    本文首先對(duì)余額寶收益率的原始序列進(jìn)行一階差分,并通過相關(guān)函數(shù)來分析序列的平穩(wěn)性。其次,按照序列的不同性質(zhì)來選用最合適的模型,通過對(duì)模型的運(yùn)用,可以對(duì)序列進(jìn)行短期預(yù)測(cè),并將其與實(shí)際情況進(jìn)行誤差分析,結(jié)果表明,我們使用的模型預(yù)測(cè)效果可信。本文可為更多新型金融衍生品的風(fēng)險(xiǎn)函數(shù)研究提供參考價(jià)值。

    [1]王瑩. 余額寶的流動(dòng)性、收益性及風(fēng)險(xiǎn)分析[J]. 中國(guó)商貿(mào),2013(35):65-66.

    [2]張大海. “余額寶”的機(jī)會(huì)與風(fēng)險(xiǎn)[J]. 創(chuàng)新時(shí)代,2013(8):37-37.

    [3]李恒,郝國(guó)剛. 余額寶的風(fēng)險(xiǎn)及防范研究[J]. 時(shí)代金融旬刊,2014(5):60.

    [4]王雪飛,劉志偉. 基于ARIMA模型的中國(guó)鋼材市場(chǎng)價(jià)格預(yù)測(cè)[J]. 中國(guó)城市經(jīng)濟(jì),2011,(1):20-23.

    [5]肖凱成,蔣春霞,楊波. ARIMA模型在土建專業(yè)群人才需求預(yù)測(cè)中的應(yīng)用[J]. 科技信息,2013(6):42,44.

    [6]陳平,夏敏. 基于ARIMA模型的短期風(fēng)速預(yù)測(cè)方法[J]. 華章,2011(15):266-267.

    [7]谷秀娟,梁潤(rùn)平. 基于ARIMA模型的鄭州市商品住宅銷售價(jià)格預(yù)測(cè)研究[J]. 金融理論與實(shí)踐,2012(1):51-54.

    [8]鄭紅云.ARIMA模型在深圳人口與醫(yī)療需求預(yù)測(cè)的應(yīng)用[J]. 齊魯工業(yè)大學(xué)學(xué)報(bào)(自然科學(xué)版),2012,26(4):44-49.

    [9]吳秀芝,劉穎博. 基于ARIMA模型對(duì)中國(guó)人均糧食的預(yù)測(cè)[J]. 金融經(jīng)濟(jì),2010(24):109-110.

    [10]易丹輝.數(shù)據(jù)分析與Eviews應(yīng)用[M].2版.北京:中國(guó)人民大學(xué)出版社,2014.

    [11]高鐵梅.計(jì)量經(jīng)濟(jì)分析方法和建?!狤views應(yīng)用及實(shí)例[M]. 北京:清華大學(xué)出版社,2009:78-80.

    [12]王周偉,崔百勝,朱敏. 經(jīng)濟(jì)計(jì)量研究指導(dǎo)——實(shí)證分析與軟件實(shí)現(xiàn)[M]. 北京:北京大學(xué)出版社,2015.

    [13]李正宏. 波羅的海運(yùn)價(jià)指數(shù)波動(dòng)規(guī)律與預(yù)測(cè)[J]. 上海海事大學(xué)學(xué)報(bào),2004,25(4):69-72.

    [14]段洪濤,曾繁聲,李景春.一種預(yù)測(cè)頻段占用度的時(shí)間序列分析方法[J].無線電工程,2011,41(7):17-20.

    Prediction of the Yu Ebao yield based on ARIMA model

    LUO Yu-fan, LIU Xiao-rong, ZHU Shu-ya, ZHAO Wei-fan, LI Min*

    (School of Mathematics and Statistics, Shandong Normal University, Jinan 250014, China)

    ∶The research object of this paper is the yield of Yu Ebao. By selecting the yield data during May 30, 2013 to July 27, 2016, the ARIMA model was established for data processing, thus predicating the yield of Yu Ebao between July 28, 2016 and August 21, 2016. Finally, by comparing the actual results with the predicted results, the validity of the model within the allowable error range was determined.

    ∶Yu Ebao; prediction of the Yu Ebao yield ; ARIMA model

    10.3976/j.issn.1002-4026.2017.04.016

    2016-10-31

    教育部留學(xué)回國(guó)人員科研啟動(dòng)基金資助項(xiàng)目(教外司留[2015]1098)

    羅裕凡(1995—),女,研究方向?yàn)樾畔⑴c計(jì)算科學(xué)。

    *通信作者,李敏。E-mail: liminemily@sdnu.edu.cn

    O2.3

    A

    1002-4026(2017)04-0099-07

    猜你喜歡
    余額差分殘差
    2024年兩融余額月度數(shù)據(jù)
    基于雙向GRU與殘差擬合的車輛跟馳建模
    數(shù)列與差分
    2020,余額不足
    基于殘差學(xué)習(xí)的自適應(yīng)無人機(jī)目標(biāo)跟蹤算法
    基于遞歸殘差網(wǎng)絡(luò)的圖像超分辨率重建
    余額寶的感知風(fēng)險(xiǎn)
    滬港通一周成交概況 (2015.5.8—2015.5.14)
    平穩(wěn)自相關(guān)過程的殘差累積和控制圖
    河南科技(2015年8期)2015-03-11 16:23:52
    基于差分隱私的大數(shù)據(jù)隱私保護(hù)
    又粗又爽又猛毛片免费看| 日韩欧美国产一区二区入口| 成人美女网站在线观看视频| 亚洲av成人精品一区久久| 亚洲avbb在线观看| 色av中文字幕| 又爽又黄无遮挡网站| 无人区码免费观看不卡| 国内精品久久久久精免费| 国内精品一区二区在线观看| 亚洲精品成人久久久久久| 久久久久久久午夜电影| 国产精品美女特级片免费视频播放器| 国产精品永久免费网站| 在线观看66精品国产| 久久国产乱子免费精品| 亚洲av第一区精品v没综合| 国产成人aa在线观看| 老女人水多毛片| 国产精品99久久久久久久久| 日韩成人在线观看一区二区三区| 在线观看美女被高潮喷水网站 | 国产精品伦人一区二区| 免费一级毛片在线播放高清视频| 成人av一区二区三区在线看| 亚洲av成人av| 精品一区二区三区视频在线| 脱女人内裤的视频| 日本免费一区二区三区高清不卡| 国产精品三级大全| 狂野欧美白嫩少妇大欣赏| 丰满乱子伦码专区| 欧美性猛交╳xxx乱大交人| 一区二区三区免费毛片| 精品人妻一区二区三区麻豆 | 亚洲人成电影免费在线| 免费一级毛片在线播放高清视频| 此物有八面人人有两片| 淫秽高清视频在线观看| aaaaa片日本免费| 色视频www国产| 精品久久久久久久人妻蜜臀av| 99久久精品一区二区三区| 国产精品电影一区二区三区| 亚洲,欧美,日韩| 好看av亚洲va欧美ⅴa在| 人人妻人人看人人澡| 亚洲av中文字字幕乱码综合| 一本综合久久免费| 黄色日韩在线| 一区二区三区高清视频在线| 最好的美女福利视频网| 欧美bdsm另类| 国产精品一区二区三区四区久久| 亚洲自偷自拍三级| 国产精品久久久久久亚洲av鲁大| 麻豆一二三区av精品| aaaaa片日本免费| av视频在线观看入口| 亚洲综合色惰| 国产高潮美女av| 久久久久免费精品人妻一区二区| 在线国产一区二区在线| 中文字幕久久专区| 欧美高清成人免费视频www| 日韩欧美国产在线观看| 日韩大尺度精品在线看网址| 少妇的逼水好多| 亚洲欧美日韩东京热| 欧美日韩综合久久久久久 | 丁香六月欧美| 天天一区二区日本电影三级| 欧美日韩综合久久久久久 | 好看av亚洲va欧美ⅴa在| 国产免费一级a男人的天堂| 久久精品国产清高在天天线| 成人无遮挡网站| 亚洲美女搞黄在线观看 | 精品久久久久久久人妻蜜臀av| 99久久久亚洲精品蜜臀av| 黄片小视频在线播放| 亚洲aⅴ乱码一区二区在线播放| 国产v大片淫在线免费观看| 国产精品99久久久久久久久| 热99在线观看视频| 男插女下体视频免费在线播放| 高清在线国产一区| 国产一区二区三区视频了| 在线免费观看的www视频| 97热精品久久久久久| 亚洲不卡免费看| 伦理电影大哥的女人| 婷婷色综合大香蕉| 日本 欧美在线| 日本三级黄在线观看| 老熟妇乱子伦视频在线观看| 亚洲av.av天堂| 日韩欧美精品免费久久 | 麻豆国产97在线/欧美| 成年人黄色毛片网站| 又黄又爽又刺激的免费视频.| 精品不卡国产一区二区三区| 国产av不卡久久| 国产真实乱freesex| 99热精品在线国产| 午夜福利高清视频| 色综合站精品国产| 在线免费观看不下载黄p国产 | 日韩人妻高清精品专区| 观看美女的网站| 亚洲精品456在线播放app | 亚洲成人精品中文字幕电影| 亚洲精华国产精华精| 久久精品国产亚洲av天美| 精品福利观看| 国产国拍精品亚洲av在线观看| 亚洲欧美日韩高清专用| 成人一区二区视频在线观看| 一进一出抽搐gif免费好疼| 久久精品综合一区二区三区| 又爽又黄无遮挡网站| 亚洲欧美日韩东京热| 一a级毛片在线观看| 99久久无色码亚洲精品果冻| 亚洲三级黄色毛片| 成人av一区二区三区在线看| 免费观看人在逋| 99久国产av精品| 伦理电影大哥的女人| 久久久成人免费电影| 成人无遮挡网站| 亚洲av成人精品一区久久| 亚洲成人久久性| 床上黄色一级片| 国产成人av教育| 99久久精品一区二区三区| 欧美又色又爽又黄视频| 久久久久久大精品| 色综合欧美亚洲国产小说| 国产精品一区二区免费欧美| 中文字幕高清在线视频| 国产一级毛片七仙女欲春2| 99国产极品粉嫩在线观看| 久久婷婷人人爽人人干人人爱| 精品久久久久久久久av| 69人妻影院| 禁无遮挡网站| 亚洲男人的天堂狠狠| 高潮久久久久久久久久久不卡| 一本综合久久免费| 国产一区二区三区视频了| 亚洲成人中文字幕在线播放| 99国产综合亚洲精品| 午夜免费激情av| 国产精品人妻久久久久久| 欧美日韩黄片免| 丰满乱子伦码专区| 极品教师在线视频| 一级黄色大片毛片| 欧美精品国产亚洲| 成年免费大片在线观看| 中文亚洲av片在线观看爽| 哪里可以看免费的av片| 国内精品久久久久久久电影| 国产精品1区2区在线观看.| 我的女老师完整版在线观看| 欧美色欧美亚洲另类二区| 在线十欧美十亚洲十日本专区| 久久精品国产亚洲av香蕉五月| 亚洲成av人片免费观看| 特级一级黄色大片| 男女床上黄色一级片免费看| 成人鲁丝片一二三区免费| 国产精品99久久久久久久久| 精品人妻1区二区| 久久99热这里只有精品18| 久久久久久久久久黄片| 搡老妇女老女人老熟妇| 亚洲欧美日韩高清专用| 亚洲第一电影网av| 国产伦在线观看视频一区| 12—13女人毛片做爰片一| 91麻豆精品激情在线观看国产| 一夜夜www| 美女 人体艺术 gogo| 天堂网av新在线| 女同久久另类99精品国产91| 91av网一区二区| 亚洲av电影在线进入| 好男人在线观看高清免费视频| 在线播放无遮挡| 欧美bdsm另类| 少妇人妻精品综合一区二区 | 一二三四社区在线视频社区8| 美女 人体艺术 gogo| 黄色配什么色好看| 国产精品日韩av在线免费观看| 成年女人毛片免费观看观看9| 变态另类成人亚洲欧美熟女| 亚洲狠狠婷婷综合久久图片| 午夜久久久久精精品| 一进一出抽搐gif免费好疼| 一个人免费在线观看电影| 在线播放国产精品三级| 最好的美女福利视频网| 在现免费观看毛片| 亚洲经典国产精华液单 | 两性午夜刺激爽爽歪歪视频在线观看| 天堂√8在线中文| 欧美黄色淫秽网站| 一进一出抽搐gif免费好疼| 精品免费久久久久久久清纯| 成人一区二区视频在线观看| 国产精品久久久久久精品电影| 精品一区二区三区视频在线| 免费av不卡在线播放| 国产国拍精品亚洲av在线观看| 亚洲精品456在线播放app | 88av欧美| 欧美高清性xxxxhd video| 每晚都被弄得嗷嗷叫到高潮| 男女视频在线观看网站免费| 真人做人爱边吃奶动态| 国产精品免费一区二区三区在线| 如何舔出高潮| 国产精品久久久久久亚洲av鲁大| 成年人黄色毛片网站| 亚洲,欧美,日韩| 欧美黄色片欧美黄色片| 亚洲avbb在线观看| 自拍偷自拍亚洲精品老妇| 国产精品久久视频播放| 韩国av一区二区三区四区| 久久99热这里只有精品18| 少妇人妻精品综合一区二区 | 99国产综合亚洲精品| 亚洲天堂国产精品一区在线| 午夜福利18| 国产亚洲av嫩草精品影院| 久久久久免费精品人妻一区二区| 久久人人爽人人爽人人片va | 国产精品久久久久久亚洲av鲁大| 久久久成人免费电影| 热99re8久久精品国产| 欧美不卡视频在线免费观看| 久久久精品欧美日韩精品| 一本精品99久久精品77| 久久国产乱子免费精品| 51国产日韩欧美| 成人国产综合亚洲| 国产精品三级大全| 亚洲人成网站高清观看| 亚洲美女视频黄频| 女人被狂操c到高潮| 午夜福利在线观看免费完整高清在 | 69av精品久久久久久| 脱女人内裤的视频| 国产蜜桃级精品一区二区三区| aaaaa片日本免费| 亚洲中文字幕日韩| 久久亚洲真实| 一个人看视频在线观看www免费| 免费高清视频大片| 99久久成人亚洲精品观看| 国产三级在线视频| 精品人妻熟女av久视频| av福利片在线观看| 免费看美女性在线毛片视频| 中文字幕av在线有码专区| 国产免费av片在线观看野外av| x7x7x7水蜜桃| 国内少妇人妻偷人精品xxx网站| 免费搜索国产男女视频| 日韩国内少妇激情av| 精品久久国产蜜桃| 久久国产精品人妻蜜桃| 久久精品久久久久久噜噜老黄 | 成人欧美大片| www.色视频.com| 国产视频内射| 中文字幕av成人在线电影| 18禁裸乳无遮挡免费网站照片| 欧美绝顶高潮抽搐喷水| 久久久久久国产a免费观看| 又爽又黄a免费视频| 2021天堂中文幕一二区在线观| 国产精华一区二区三区| 黄色日韩在线| 狠狠狠狠99中文字幕| 丁香六月欧美| 亚州av有码| 人妻久久中文字幕网| 欧美+亚洲+日韩+国产| 超碰av人人做人人爽久久| 久久久久久久亚洲中文字幕 | 精品久久久久久久久久免费视频| 久久这里只有精品中国| 欧美另类亚洲清纯唯美| 无遮挡黄片免费观看| 精品熟女少妇八av免费久了| 三级国产精品欧美在线观看| 18禁裸乳无遮挡免费网站照片| 免费看光身美女| a级毛片免费高清观看在线播放| 中文在线观看免费www的网站| 丁香欧美五月| 美女高潮的动态| 国产高潮美女av| 十八禁人妻一区二区| 少妇熟女aⅴ在线视频| 久久久久久久午夜电影| 看免费av毛片| 五月玫瑰六月丁香| 级片在线观看| 男女下面进入的视频免费午夜| 一区二区三区免费毛片| 又黄又爽又刺激的免费视频.| 婷婷精品国产亚洲av| 久久精品国产清高在天天线| 国产成人影院久久av| 亚洲在线观看片| 好男人电影高清在线观看| 欧美性感艳星| 国产v大片淫在线免费观看| 亚洲av第一区精品v没综合| 国产一区二区在线av高清观看| 啦啦啦韩国在线观看视频| 亚洲av二区三区四区| 我的老师免费观看完整版| 丰满人妻一区二区三区视频av| 国产淫片久久久久久久久 | 女人被狂操c到高潮| 在线播放无遮挡| 亚洲在线自拍视频| 国产探花在线观看一区二区| 69人妻影院| 中出人妻视频一区二区| 两性午夜刺激爽爽歪歪视频在线观看| 人人妻,人人澡人人爽秒播| 我要搜黄色片| 国产成人啪精品午夜网站| 日本成人三级电影网站| 精品乱码久久久久久99久播| 给我免费播放毛片高清在线观看| 免费人成在线观看视频色| 日本五十路高清| 69av精品久久久久久| 国产aⅴ精品一区二区三区波| 色哟哟哟哟哟哟| 亚洲av电影在线进入| 国产单亲对白刺激| 久久国产精品影院| 哪里可以看免费的av片| 97超级碰碰碰精品色视频在线观看| 欧美色视频一区免费| 国产久久久一区二区三区| 国产真实伦视频高清在线观看 | 蜜桃亚洲精品一区二区三区| 国产真实伦视频高清在线观看 | 国产精品不卡视频一区二区 | 一个人看视频在线观看www免费| 久久久精品欧美日韩精品| 天堂av国产一区二区熟女人妻| 亚洲精品色激情综合| 乱码一卡2卡4卡精品| 热99在线观看视频| 99久久成人亚洲精品观看| 搡女人真爽免费视频火全软件 | 嫩草影院新地址| 久久久久久久久中文| 美女黄网站色视频| 国产成人a区在线观看| 久久国产精品影院| 日韩高清综合在线| 亚洲成人免费电影在线观看| 久久性视频一级片| 亚洲中文字幕一区二区三区有码在线看| 99riav亚洲国产免费| 小说图片视频综合网站| 男插女下体视频免费在线播放| 欧美日韩综合久久久久久 | 男人和女人高潮做爰伦理| 精品国产亚洲在线| 午夜精品在线福利| 一进一出抽搐gif免费好疼| 久9热在线精品视频| 悠悠久久av| 一级作爱视频免费观看| 国产爱豆传媒在线观看| 一级黄片播放器| 此物有八面人人有两片| 色视频www国产| 村上凉子中文字幕在线| 又爽又黄无遮挡网站| 国产午夜精品久久久久久一区二区三区 | 中文字幕av成人在线电影| 国产精品美女特级片免费视频播放器| 国产精品久久久久久亚洲av鲁大| 美女高潮的动态| 嫩草影院精品99| 日日摸夜夜添夜夜添小说| 成年女人永久免费观看视频| 中文字幕高清在线视频| 一级毛片久久久久久久久女| 亚洲国产精品合色在线| 一区二区三区免费毛片| 精品国产亚洲在线| 能在线免费观看的黄片| 国产男靠女视频免费网站| 美女高潮的动态| 此物有八面人人有两片| 欧美中文日本在线观看视频| 久久精品国产自在天天线| 色av中文字幕| 香蕉av资源在线| 国产aⅴ精品一区二区三区波| 97超视频在线观看视频| 成人欧美大片| 极品教师在线免费播放| 又黄又爽又免费观看的视频| 国产精品久久视频播放| 国产一区二区三区视频了| 1024手机看黄色片| 欧美又色又爽又黄视频| 亚洲欧美精品综合久久99| 日韩中字成人| 日韩欧美国产在线观看| 变态另类丝袜制服| 久久久久免费精品人妻一区二区| 亚洲精品粉嫩美女一区| 91麻豆av在线| 国产老妇女一区| 人妻制服诱惑在线中文字幕| 国产午夜精品久久久久久一区二区三区 | 久久精品影院6| 欧美3d第一页| 精品国产亚洲在线| 51午夜福利影视在线观看| 国内精品久久久久精免费| eeuss影院久久| 久久天躁狠狠躁夜夜2o2o| 天堂网av新在线| 国产真实伦视频高清在线观看 | 国产成人aa在线观看| 麻豆成人av在线观看| a级毛片a级免费在线| 久久久久精品国产欧美久久久| 自拍偷自拍亚洲精品老妇| 日韩欧美精品免费久久 | 精品免费久久久久久久清纯| 国产在视频线在精品| 91字幕亚洲| 桃色一区二区三区在线观看| 日本 av在线| 久久国产精品人妻蜜桃| 欧美一区二区国产精品久久精品| 国产探花极品一区二区| 性色avwww在线观看| 色哟哟哟哟哟哟| 非洲黑人性xxxx精品又粗又长| 国产真实伦视频高清在线观看 | 欧美+日韩+精品| 日韩av在线大香蕉| 全区人妻精品视频| 国产精品美女特级片免费视频播放器| 自拍偷自拍亚洲精品老妇| 国产精品伦人一区二区| 国产蜜桃级精品一区二区三区| 男女那种视频在线观看| 麻豆国产97在线/欧美| 色综合婷婷激情| 久久精品综合一区二区三区| 亚洲人成伊人成综合网2020| 啦啦啦韩国在线观看视频| 嫁个100分男人电影在线观看| 狠狠狠狠99中文字幕| 国产精品亚洲av一区麻豆| 免费看光身美女| 又粗又爽又猛毛片免费看| eeuss影院久久| 97超视频在线观看视频| 99在线人妻在线中文字幕| 乱人视频在线观看| 99在线视频只有这里精品首页| 日韩大尺度精品在线看网址| 永久网站在线| 亚洲av免费在线观看| 男插女下体视频免费在线播放| av黄色大香蕉| 99久久久亚洲精品蜜臀av| 一个人免费在线观看的高清视频| 成人鲁丝片一二三区免费| 少妇高潮的动态图| 在线免费观看的www视频| 麻豆成人午夜福利视频| 人人妻人人看人人澡| 18禁黄网站禁片免费观看直播| 毛片一级片免费看久久久久 | 在线看三级毛片| 美女免费视频网站| 国产色爽女视频免费观看| 亚洲精品粉嫩美女一区| 国产探花极品一区二区| 高清毛片免费观看视频网站| 又粗又爽又猛毛片免费看| 国产v大片淫在线免费观看| 国产精品久久久久久精品电影| 看免费av毛片| 日韩欧美免费精品| 精品不卡国产一区二区三区| 国产美女午夜福利| 亚洲国产精品999在线| 欧美高清性xxxxhd video| 国产精品一区二区三区四区久久| 91九色精品人成在线观看| 男女视频在线观看网站免费| 每晚都被弄得嗷嗷叫到高潮| 久久久国产成人精品二区| 国产欧美日韩一区二区三| 日日夜夜操网爽| 性色av乱码一区二区三区2| 又紧又爽又黄一区二区| 狂野欧美白嫩少妇大欣赏| 3wmmmm亚洲av在线观看| 国产高清三级在线| av在线天堂中文字幕| 国产黄片美女视频| 国产淫片久久久久久久久 | 欧美日韩福利视频一区二区| 亚洲国产精品sss在线观看| 在线播放国产精品三级| 国产一区二区亚洲精品在线观看| 国内精品一区二区在线观看| 亚洲av免费高清在线观看| av天堂在线播放| 国产精品永久免费网站| 成年女人永久免费观看视频| 国产精品一区二区三区四区免费观看 | 丁香欧美五月| 亚洲在线自拍视频| 欧美日韩乱码在线| 动漫黄色视频在线观看| 国产亚洲欧美98| 国产精品久久久久久精品电影| 真人做人爱边吃奶动态| 国产精品久久久久久人妻精品电影| 动漫黄色视频在线观看| 夜夜夜夜夜久久久久| 亚洲综合色惰| 亚洲午夜理论影院| 国产精品久久视频播放| 国产三级中文精品| 欧美国产日韩亚洲一区| 综合色av麻豆| 在线观看66精品国产| 美女 人体艺术 gogo| 婷婷六月久久综合丁香| 亚洲中文日韩欧美视频| 一区二区三区免费毛片| 国产精品影院久久| av黄色大香蕉| 欧美性感艳星| 人人妻,人人澡人人爽秒播| 成人亚洲精品av一区二区| 长腿黑丝高跟| 757午夜福利合集在线观看| 真人一进一出gif抽搐免费| 97人妻精品一区二区三区麻豆| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 校园春色视频在线观看| 麻豆国产97在线/欧美| 激情在线观看视频在线高清| 可以在线观看的亚洲视频| eeuss影院久久| av专区在线播放| 欧美zozozo另类| 免费在线观看影片大全网站| 国产精品伦人一区二区| 欧美黄色片欧美黄色片| 两性午夜刺激爽爽歪歪视频在线观看| 男人和女人高潮做爰伦理| 91av网一区二区| 怎么达到女性高潮| 亚洲aⅴ乱码一区二区在线播放| 久久精品91蜜桃| 可以在线观看的亚洲视频| 最好的美女福利视频网| 老司机福利观看| 香蕉av资源在线| 精品国内亚洲2022精品成人| 中文字幕熟女人妻在线| 99国产精品一区二区三区| 亚洲精华国产精华精| 中文亚洲av片在线观看爽| 国产欧美日韩一区二区精品| 亚洲av成人不卡在线观看播放网| 日本撒尿小便嘘嘘汇集6| 18禁黄网站禁片免费观看直播| 亚洲精品乱码久久久v下载方式| 免费电影在线观看免费观看| 亚洲三级黄色毛片| 九九热线精品视视频播放| 亚洲,欧美,日韩| 免费搜索国产男女视频| 亚洲午夜理论影院| 日韩精品青青久久久久久| 亚洲专区国产一区二区| 波多野结衣巨乳人妻| 成年女人看的毛片在线观看| 亚洲乱码一区二区免费版| 国产三级黄色录像| 老司机午夜福利在线观看视频| 亚洲色图av天堂| 琪琪午夜伦伦电影理论片6080| 校园春色视频在线观看|