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

    基于Matlab的股票市場(chǎng)收益率波動(dòng)分析實(shí)驗(yàn)

    2014-06-05 09:51:47徐德玲龐世達(dá)孫凱旋
    關(guān)鍵詞:方差收益率波動(dòng)

    叢 超,徐德玲,龐世達(dá),孫凱旋

    基于Matlab的股票市場(chǎng)收益率波動(dòng)分析實(shí)驗(yàn)

    叢 超,徐德玲,龐世達(dá),孫凱旋

    (重慶理工大學(xué) 電子信息與自動(dòng)化學(xué)院,重慶 400054)

    針對(duì)金融風(fēng)暴背景下的股票市場(chǎng)價(jià)格的波動(dòng)特性,應(yīng)用數(shù)學(xué)分析、經(jīng)濟(jì)統(tǒng)計(jì)與計(jì)量知識(shí),對(duì)中國(guó)上海、深圳股票綜合指數(shù)2007~2009年的數(shù)據(jù)進(jìn)行實(shí)驗(yàn)分析,并利用Matlab金融分析工具箱以及廣義自回歸異方差模型編程建模,實(shí)現(xiàn)對(duì)股票市場(chǎng)收益率的分析和預(yù)測(cè)。結(jié)果表明,股票市場(chǎng)收益率序列的波動(dòng)有顯著的異方差性。

    股票市場(chǎng);時(shí)間序列分析;廣義自回歸異方差模型;Matlab編程

    金融時(shí)間序列收益率的波動(dòng)是動(dòng)態(tài)變化的,不同金融市場(chǎng)的波動(dòng)還存在波動(dòng)溢出。股票作為一種重要的金融產(chǎn)品,其價(jià)格行為理論是整個(gè)市場(chǎng)金融理論的基礎(chǔ),股市價(jià)格行為(behavior of stockmarket prices)一詞最早由Fama提出[1],其核心含義是指股價(jià)的行為方式,即變動(dòng)規(guī)律。人們對(duì)股票預(yù)測(cè)也提出了各種不同的方法[2],從最初的圖表分析、技術(shù)指標(biāo)(成交量曲線圖、K線圖、移動(dòng)平均線等),到建立數(shù)學(xué)模型的方法等。隨著經(jīng)濟(jì)學(xué)、數(shù)學(xué)以及計(jì)算機(jī)科學(xué)的發(fā)展,針對(duì)金融市場(chǎng)時(shí)間序列進(jìn)行統(tǒng)計(jì)性的建模與分析已經(jīng)成為一項(xiàng)跨學(xué)科的分析課題,成為國(guó)內(nèi)外研究的重要方向。

    本文通過(guò)Matlab編程建立模型,針對(duì)2007年1月1日到2009年12月31日中國(guó)上海、深圳股票綜合指數(shù)數(shù)據(jù)進(jìn)行實(shí)驗(yàn)分析。首先,利用Matlab金融工具箱對(duì)股票市場(chǎng)的收益率曲線進(jìn)行計(jì)算,并檢驗(yàn)金融分指數(shù)序列的平穩(wěn)性與波動(dòng)性;其次,通過(guò)使用自相關(guān)(auto correlation function,ACF)和偏相關(guān)(partial auto correlation function,PACF)分析的方法檢驗(yàn)序列的自相關(guān)性;隨后,采用Q檢驗(yàn)和自回歸異方差檢驗(yàn)(auto regressive conditional heteroskedasticity,ARCH)進(jìn)行平穩(wěn)性的驗(yàn)證,并使用GARCH(1,1)模型對(duì)收益率曲線進(jìn)行建模分析;最后,本文對(duì)模型的優(yōu)缺點(diǎn)進(jìn)行了評(píng)價(jià),并給出了推廣與改進(jìn)的建議。

    1 基本概念

    1.1 一般時(shí)間序列的統(tǒng)計(jì)特性

    對(duì)于一般的時(shí)間序列,主要的統(tǒng)計(jì)特征量有以下三種[3]。

    1)均值函數(shù)

    2)自協(xié)方差函數(shù)

    式中:當(dāng)t=s時(shí),γt,t為時(shí)間序列的方差,如下式所示。

    3)自相關(guān)函數(shù)

    可用上述三種特征來(lái)分析一個(gè)時(shí)間序列的基本特性,如平穩(wěn)性檢驗(yàn)[4]、正態(tài)性、獨(dú)立性等。

    除此之外,在實(shí)驗(yàn)中通常還會(huì)用到以下幾種類(lèi)型的特性,用于支持或細(xì)化分析時(shí)間序列建模。

    1)樣本峰度,用于度量時(shí)間序列偏離某分布的情況[5]。例如,正態(tài)分布的峰度為3。

    2)樣本偏差,用于衡量時(shí)間序列的對(duì)稱性。若偏差為負(fù),則均值左側(cè)的離散度比右側(cè)高,反之亦然。

    1.2 股票指數(shù)時(shí)間序列的特征及建模方法

    股票指數(shù)作為時(shí)間序列的一種,除了具有時(shí)間序列的特性之外,本身有著金融產(chǎn)品的特殊性。

    1)非平穩(wěn)性。通常來(lái)說(shuō),股票指數(shù)的時(shí)間序列是非平穩(wěn)的,通過(guò)對(duì)股票指數(shù)的三大時(shí)間序列特征的分析,可以很容易得到股票指數(shù)序列具有非平穩(wěn)的特性。

    2)收益率序列具有波動(dòng)聚集性。早在20世紀(jì)60年代,F(xiàn)ama針對(duì)股市價(jià)格行為分析得出結(jié)論,在某一段時(shí)間內(nèi),收益率序列存在持續(xù)偏高或偏低的情況,這就是波動(dòng)聚集性。

    3)收益率序列具有異方差性。異方差性指的是針對(duì)時(shí)間序列進(jìn)行經(jīng)典回歸分析時(shí),隨機(jī)誤差項(xiàng)μt由于金融產(chǎn)品容易受到隨機(jī)因素的影響,比如政策變動(dòng)、自然災(zāi)害或金融危機(jī)等,同時(shí),人的不確定的社會(huì)行為也會(huì)對(duì)股票指數(shù)序列造成影響。因此,針對(duì)金融產(chǎn)品的分析通常不能忽視異方差性。

    在早些時(shí)候,學(xué)者對(duì)股票價(jià)格時(shí)間序列的建模大多為幾何布朗運(yùn)動(dòng)模型[6]。這種建模方法的主要缺陷在于,難以解釋股指時(shí)間序列具有方差時(shí)變、波動(dòng)集群的特征。而股票指數(shù)序列本身具有非平穩(wěn)性,許多傳統(tǒng)的時(shí)間序列建模方法就很難應(yīng)用在股票指數(shù)上。

    1.3 GARCH模型的基本概念

    GARCH模型是Bollerslev等人于1986年提出的廣義自回歸條件異方差模型(generalized auto regressive conditional heteroskedasticity,GARCH),主要針對(duì)金融時(shí)間序列收益率所量體訂做的回歸模型。除去和普通回歸模型的相同之處,GARCH將方差和條件方差區(qū)別開(kāi)來(lái),對(duì)條件方差進(jìn)行自回歸建模,適用于波動(dòng)性的分析和預(yù)測(cè)。

    GARCH模型可表示為GARCH(p,q),模型假定當(dāng)前時(shí)間的條件方差依賴于其滯后值ht-i和殘差的滯后值當(dāng)期的方差決定于三個(gè)因素:常數(shù)項(xiàng)a0、前一期的殘差和前一期的預(yù)測(cè)方差。

    均值方程:

    條件方差方程:

    式(7)給出的均值方程是一個(gè)以Xt以前信息為基礎(chǔ)的向前預(yù)測(cè)方差,稱為條件均值方程。在許多文獻(xiàn)中以自回歸(AR)方程的形式給出。

    式(8)給出的方程中:a0為常數(shù)項(xiàng),其余的部分以自回歸滑動(dòng)平均(ARMA)的形式給出。此方程又稱條件方差方程,說(shuō)明時(shí)間序列條件方差的變化特征。

    式(9)中,N(0,1)為正態(tài)分布。

    1.4 Matlab在金融時(shí)間序列分析實(shí)驗(yàn)中的應(yīng)用

    Matlab是科學(xué)計(jì)算與實(shí)驗(yàn)驗(yàn)證的重要工具,在針對(duì)時(shí)間序列統(tǒng)計(jì)特性的分析以及時(shí)間序列建模分析的過(guò)程中,Matlab提供了大量統(tǒng)計(jì)分析函數(shù),如表1所示。

    表1 時(shí)間序列統(tǒng)計(jì)特征計(jì)算函數(shù)

    除統(tǒng)計(jì)函數(shù),Matlab也提供了金融分析工具箱,能夠進(jìn)行序列的平穩(wěn)性檢驗(yàn)、GARCH檢驗(yàn)、GARCH模型擬合等。

    2 具體實(shí)驗(yàn)分析與Matlab編碼實(shí)現(xiàn)

    2.1 樣本數(shù)據(jù)的選取

    本文選取的股市指數(shù)樣本數(shù)據(jù)來(lái)源于中國(guó)知網(wǎng)國(guó)際金融統(tǒng)計(jì)數(shù)據(jù)。數(shù)據(jù)庫(kù)從2007年1月1日到2009年12月31日,共3 000個(gè)觀測(cè)值。并且運(yùn)用Matlab R2007軟件進(jìn)行分析驗(yàn)證。序列的基本信息如圖1所示。

    圖1 深圳股票指數(shù)

    圖1 中,橫坐標(biāo)表示時(shí)間,縱坐標(biāo)從上至下依次為隨時(shí)間變化的收盤(pán)價(jià)及成交量(萬(wàn)股)。

    首先,針對(duì)收盤(pán)價(jià)進(jìn)行直接分析可以看出,2007年初的收盤(pán)價(jià)較為平穩(wěn),整體呈現(xiàn)上升趨勢(shì),并在2007年10月左右到達(dá)峰頂。但在中間一段,也就是2008年金融危機(jī)爆發(fā)后,波動(dòng)性明顯,所有指數(shù)均有很大的下降趨勢(shì),而到了2009年,數(shù)據(jù)有所回復(fù),說(shuō)明在2008年的金融危機(jī)情況下,股市呈下滑趨勢(shì)。

    2.2 數(shù)據(jù)的統(tǒng)計(jì)特征

    為了對(duì)股票序列進(jìn)一步分析,以成交量為研究對(duì)象,就日收益、日對(duì)數(shù)收益率、峰度、偏度、方差等進(jìn)行了計(jì)算。

    圖2和圖3分別為深成指數(shù)日對(duì)數(shù)收益率及收益率結(jié)果統(tǒng)計(jì)圖。

    圖2 深成指數(shù)日對(duì)數(shù)收益率

    從圖2和圖3可以看出:日收益率以及對(duì)數(shù)收益率都具有波動(dòng)聚集性,深成指數(shù)日對(duì)數(shù)收益率集中在(-0.05,0.05)之間,在波動(dòng)較大的幾個(gè)時(shí)段上收益率絕對(duì)值達(dá)到了0.10,在到達(dá)頂峰值后又迅速回落。在樣本期間,平均收益率(mediam),即均值為0.002 8,偏度(skewness)為-0.402 9,峰度(kurtosis)為4.202 8,數(shù)據(jù)分布向左偏移,有一個(gè)沉重的尾巴。這說(shuō)明從整體上看,收益率低于其均值的時(shí)候較多,波動(dòng)性強(qiáng),頭部和中間包含了大量統(tǒng)計(jì)信息。

    隨后對(duì)收益率曲線進(jìn)行自相關(guān)與偏相關(guān)分析,結(jié)果分別如圖4和圖5所示。

    從日收益率的ACF(自相關(guān)性)和PACF(偏相關(guān)性)看不出日收益率具有明顯自相關(guān)性。我們可以用Q檢驗(yàn)和ARCH檢驗(yàn)來(lái)檢驗(yàn)日收益率序列的自相關(guān)性。

    2.3 相關(guān)性檢驗(yàn)

    相關(guān)性檢驗(yàn)的Matlab實(shí)現(xiàn)如下所示。

    1)Q檢驗(yàn)

    2)ARCH檢驗(yàn)

    garchplot(innovations,sigmas,shouyi)

    >>[H,pValue,Stat,CriticalValue]=…

    archtest(seq-mean(seq),[10 1520]′,0.05);

    圖3 日收益率統(tǒng)計(jì)圖及計(jì)算結(jié)果圖

    圖4 深圳股票自相關(guān)指數(shù)

    圖5 深成指數(shù)偏相關(guān)序列

    相關(guān)性檢驗(yàn)的結(jié)果分別如表2和表3所示。

    從以上兩種檢驗(yàn)結(jié)果可以發(fā)現(xiàn):在Q檢驗(yàn)中,H=0表明,日收益率殘差和收益波動(dòng)不存在自相關(guān)性;ARCH檢驗(yàn)中,H=1表明,收益率序列具有ARCH效應(yīng)。均值方程的ht項(xiàng)的系數(shù)是0.0012,在5的顯著性水平下顯著大于0,這表明日收益率與市場(chǎng)風(fēng)險(xiǎn)水平呈弱正相關(guān),驗(yàn)證了高風(fēng)險(xiǎn)對(duì)應(yīng)于高收益的投資組合理論。

    表2 Q檢驗(yàn)結(jié)果

    表3 ARCH檢驗(yàn)

    1.0 0.000 3 41.176 0 24.995 8 1.0 0.003 4 41.331 8 31.410 4

    2.4 參數(shù)估計(jì)

    我們利用Matlab金融分析工具箱進(jìn)行GARCH建模,Matlab實(shí)現(xiàn)如下。

    1)GARCH擬合

    [coeff,errors,LLF,innovations,sigmas,summary]=garchfit(seq);

    2)GARCH參數(shù)

    garchdisp(coeff,errors)

    擬合的參數(shù)如表4所示。

    條件異方差方程的參數(shù):常數(shù)項(xiàng)a0,K值,GARCH項(xiàng)(前一期的方差的預(yù)測(cè)值ht-1),ARCH項(xiàng)。表4中,第一列為所得值,第二列是標(biāo)準(zhǔn)差,第三列是T檢驗(yàn)的數(shù)值。

    由計(jì)算結(jié)果用GARCH(1,1)擬合觀察值,得出均值方程和條件方差方程。

    均值方程:

    條件方差方程:

    2.5 估計(jì)效果分析

    GARCH(1,1)模型的β1的系數(shù)都比較大且通過(guò)了顯著性檢驗(yàn),說(shuō)明股價(jià)波動(dòng)具有“長(zhǎng)期記憶性”,即過(guò)去價(jià)格的波動(dòng)與其無(wú)限長(zhǎng)期價(jià)格波動(dòng)的大小都有關(guān)系。條件方差方程中,系數(shù)a1和β1都顯著為正,說(shuō)明過(guò)去的波動(dòng)對(duì)市場(chǎng)未來(lái)波動(dòng)有著正向而減緩的影響,從而使股市波動(dòng)出現(xiàn)群聚性現(xiàn)象。a1+β1都接近于1,這說(shuō)明股市波動(dòng)對(duì)外部沖擊的反應(yīng)函數(shù)以一個(gè)相對(duì)較慢的速度遞減,股市一旦出現(xiàn)大的波動(dòng)在短期內(nèi)很難消除。另外,由于GARCH(1,1)模型中a1+β1小于1,說(shuō)明收益率條件方差序列是平穩(wěn)的,模型具有可預(yù)測(cè)性。

    除了參數(shù)估計(jì)值和標(biāo)準(zhǔn)差,garchfit同時(shí)計(jì)算了殘差(innovations)、條件標(biāo)準(zhǔn)差(conditional standard deviations)。從圖6中可以看出,殘差出現(xiàn)了波動(dòng)聚集現(xiàn)象,顯示了很強(qiáng)的非平穩(wěn)性。經(jīng)過(guò)GARCH模型擬合之后,條件標(biāo)準(zhǔn)殘差的波動(dòng)聚集現(xiàn)象已經(jīng)不明顯了。

    圖6 殘差(a)與條件標(biāo)準(zhǔn)差(b)

    3 結(jié)束語(yǔ)

    3.1 本文的主要貢獻(xiàn)

    本文在對(duì)股票指數(shù)進(jìn)行分析時(shí)選用了GARCH模型,準(zhǔn)確地模擬時(shí)間序列變量的波動(dòng)性的變化。本文所選用的模型在金融工程學(xué)的實(shí)證研究中應(yīng)用廣泛,尤其是應(yīng)用在風(fēng)險(xiǎn)價(jià)值(value at risk)理論中,使人們能更加準(zhǔn)確地把握風(fēng)險(xiǎn)。

    GARCH用方差以預(yù)測(cè)時(shí)間序列方差:(1)估計(jì)方差,衡量風(fēng)險(xiǎn);(2)可以計(jì)算均值方差中變量的置信區(qū)間;(3)對(duì)條件異方差正確估計(jì)可以使估計(jì)參數(shù)更準(zhǔn)確,特別適用于波動(dòng)性的分析和預(yù)測(cè),這樣的分析對(duì)投資者的決策能起到非常重要的指導(dǎo)性作用,其意義很多時(shí)候超過(guò)了對(duì)數(shù)值本身的分析和預(yù)測(cè)。

    可以預(yù)見(jiàn),未來(lái)的研究將會(huì)在方法論和工具論兩個(gè)方向進(jìn)一步展開(kāi),特別是其應(yīng)用研究還在不斷拓展,伴隨著市場(chǎng)微觀結(jié)構(gòu)理論的成熟,采用ARCH模型來(lái)模擬波動(dòng)性,將會(huì)對(duì)期貨交易制度設(shè)計(jì)、風(fēng)險(xiǎn)控制制度設(shè)計(jì)和投資組合風(fēng)險(xiǎn)管理策略研究,提供一個(gè)更為廣闊的研究空間。

    3.2 本文的缺陷及改進(jìn)的方法

    GARCH模型適合在計(jì)算量不大時(shí),方便地描述了高階的ARCH過(guò)程,因而具有更大的適用性。但GARCH(1,1)模型在應(yīng)用于資產(chǎn)定價(jià)方面存在以下的不足:

    1)GARCH模型不能解釋股票收益和收益變化波動(dòng)之間出現(xiàn)的負(fù)相關(guān)現(xiàn)象。GARCH(p,q)模型假定條件方差是滯后殘差平方的函數(shù),因此,殘差的符號(hào)不影響波動(dòng),即條件方差對(duì)正的價(jià)格變化和負(fù)的價(jià)格變化的反應(yīng)是對(duì)稱的。然而,在經(jīng)驗(yàn)研究中發(fā)現(xiàn),當(dāng)利空消息出現(xiàn)時(shí),即預(yù)期股票收益會(huì)下降時(shí),波動(dòng)趨向于增大;當(dāng)利好消息出現(xiàn)時(shí),即預(yù)期股票收益會(huì)上升時(shí),波動(dòng)趨向于減小。GARCH模型不能解釋這種非對(duì)稱現(xiàn)象。

    2)GARCH(1,1)模型為了保證非負(fù),假定式(2)中所有系數(shù)均大于零。這些約束隱含著的任何滯后項(xiàng)增大都會(huì)增加,因而排除了隨機(jī)波動(dòng)行為,使得在估計(jì)GARCH模型時(shí)可能出現(xiàn)震蕩現(xiàn)象。

    3)此外,針對(duì)波動(dòng)溢出分析,原始GARCH模型僅能支持單一金融市場(chǎng),難以對(duì)不同市場(chǎng)之間的相互影響進(jìn)行分析。

    3.3 改進(jìn)的方法

    如果要對(duì)波動(dòng)溢出做進(jìn)一步分析,可以考慮如下方法:

    1)主成分分析與獨(dú)立成分分析。通過(guò)多個(gè)金融市場(chǎng)的收益率波動(dòng)數(shù)據(jù)構(gòu)成協(xié)方差矩陣并計(jì)算主成分和獨(dú)立成分,并與GARCH模型相結(jié)合的方法來(lái)進(jìn)行分析。

    2)建立多個(gè)金融市場(chǎng)的時(shí)間序列矩陣,引入多變量GARCH模型或矢量GARCH模型(bekkmgarch),才能達(dá)到理想的分析結(jié)果。

    [1]Fama E F.The behavior of stock-market prices[J].The Journal of Business,1965,38(1):34-105.

    [2]文雪梅,柏滿迎.股票收益率預(yù)測(cè)模型的比較[J].現(xiàn)代商業(yè),2009(11):26.

    [3]趙耀軍.時(shí)間序列分析[J].山西冶金,2012(6):19.

    [4]張樹(shù)京,齊立心.時(shí)間序列分析簡(jiǎn)明教程[M].北京:清華大學(xué)出版社,2003.

    [5]張善文,雷英杰,馮有前.Matlab在時(shí)間序列分析中的應(yīng)用[M].西安:西安電子科技大學(xué)出版社,2007.

    [6]林建華,王福昌,馮敬海.股價(jià)波動(dòng)的指數(shù)O2U過(guò)程模型[J].經(jīng)濟(jì)數(shù)學(xué),2000,17(4):29-32.

    Volatility M odelling of Chinese Stock M arket Based on M atlab

    CONG Chao,XU Deling,PANG Shida,SUN Kaixuan
    (School of Electronic Information and Automation,Chongqing University of Technology,Chongqing 400054,China)

    On the background of a global financial turmoil,mathematical,statistical and economic analysiswas applied in the experiment of analyzing the volatility feature of Shanghaiand Shenzhen stock indexes from 2007 to2009.Then,a Matlab based financial toolbox and GARCH(generalized auto regressive conditional heteroskedasticity)modelingmethodswere used in the experiment.The result showed an obvious heteroscedasticity feature in return series and a very good fitting with GARCH model.

    stock market;time series analysis;generalized auto regressive conditional heteroskedasticity;Matlab program

    TP311;F832.5

    A

    10.3969/j.issn.1672-4550.2014.05.022

    2013-06-06;修改日期:2014-07-12

    叢 超(1984-),男,碩士,助理實(shí)驗(yàn)師,研究方向:圖像處理、軟件開(kāi)發(fā)、高性能并行計(jì)算。

    猜你喜歡
    方差收益率波動(dòng)
    1年期國(guó)債收益率跌至1%
    銀行家(2025年1期)2025-02-08 00:00:00
    方差怎么算
    概率與統(tǒng)計(jì)(2)——離散型隨機(jī)變量的期望與方差
    羊肉價(jià)回穩(wěn) 后期不會(huì)大幅波動(dòng)
    計(jì)算方差用哪個(gè)公式
    微風(fēng)里優(yōu)美地波動(dòng)
    2019年國(guó)內(nèi)外油價(jià)或?qū)⒉▌?dòng)加劇
    方差生活秀
    干濕法SO2排放波動(dòng)對(duì)比及分析
    国产亚洲精品综合一区在线观看 | 9色porny在线观看| 亚洲精品久久午夜乱码| 亚洲中文日韩欧美视频| 久久精品91蜜桃| 久久天堂一区二区三区四区| 久久中文看片网| 国产一区二区三区在线臀色熟女 | 亚洲精品国产一区二区精华液| 国产在线精品亚洲第一网站| 日韩精品中文字幕看吧| 国产在线精品亚洲第一网站| www日本在线高清视频| 国产一区二区三区视频了| 老司机靠b影院| 午夜免费成人在线视频| 久久人人97超碰香蕉20202| av在线播放免费不卡| 在线天堂中文资源库| cao死你这个sao货| 老司机午夜十八禁免费视频| 国产单亲对白刺激| 脱女人内裤的视频| 亚洲va日本ⅴa欧美va伊人久久| 搡老熟女国产l中国老女人| 女性被躁到高潮视频| 久久久久久人人人人人| 少妇被粗大的猛进出69影院| 欧美性长视频在线观看| 老汉色∧v一级毛片| 动漫黄色视频在线观看| 免费少妇av软件| 亚洲av成人av| 成人免费观看视频高清| 国产成人影院久久av| 女警被强在线播放| 91老司机精品| 在线观看66精品国产| 午夜免费成人在线视频| 日本欧美视频一区| 日日干狠狠操夜夜爽| 男人操女人黄网站| 在线播放国产精品三级| 777久久人妻少妇嫩草av网站| 一进一出抽搐gif免费好疼 | 69精品国产乱码久久久| 精品欧美一区二区三区在线| 亚洲国产欧美网| 一级a爱视频在线免费观看| 亚洲 国产 在线| 精品一品国产午夜福利视频| 在线观看免费视频网站a站| 国产精品九九99| 精品福利永久在线观看| 色婷婷av一区二区三区视频| 亚洲成av片中文字幕在线观看| www国产在线视频色| 亚洲在线自拍视频| 久久国产亚洲av麻豆专区| av国产精品久久久久影院| 一区二区三区国产精品乱码| 亚洲国产精品合色在线| 亚洲成人免费av在线播放| 深夜精品福利| 在线看a的网站| 久99久视频精品免费| 国产人伦9x9x在线观看| 亚洲精品一卡2卡三卡4卡5卡| 不卡一级毛片| 亚洲国产精品999在线| 精品国产一区二区三区四区第35| 日本 av在线| 亚洲国产毛片av蜜桃av| 免费观看精品视频网站| 久久人妻av系列| 高潮久久久久久久久久久不卡| 日日爽夜夜爽网站| 搡老熟女国产l中国老女人| 久久精品亚洲熟妇少妇任你| 黑人巨大精品欧美一区二区mp4| 97人妻天天添夜夜摸| 激情在线观看视频在线高清| 亚洲欧美日韩另类电影网站| 亚洲午夜精品一区,二区,三区| 亚洲午夜精品一区,二区,三区| 50天的宝宝边吃奶边哭怎么回事| 可以在线观看毛片的网站| 精品国产国语对白av| 久久国产亚洲av麻豆专区| av免费在线观看网站| 久久久久久久久免费视频了| 国产高清videossex| 制服人妻中文乱码| 丰满迷人的少妇在线观看| 黄片播放在线免费| 久久影院123| 最近最新中文字幕大全电影3 | 精品福利永久在线观看| 国产免费现黄频在线看| 亚洲一区二区三区色噜噜 | 人妻丰满熟妇av一区二区三区| 欧美成人午夜精品| 国产午夜精品久久久久久| 国产极品粉嫩免费观看在线| 亚洲专区字幕在线| 午夜a级毛片| 欧美成人性av电影在线观看| 校园春色视频在线观看| 啦啦啦在线免费观看视频4| 欧美日韩黄片免| 久久精品国产清高在天天线| 国产精品 国内视频| 亚洲va日本ⅴa欧美va伊人久久| 手机成人av网站| 天天躁狠狠躁夜夜躁狠狠躁| 久久精品国产亚洲av高清一级| 国产精品久久电影中文字幕| 欧美日韩黄片免| 88av欧美| 免费不卡黄色视频| 9191精品国产免费久久| 日韩欧美一区视频在线观看| 亚洲人成77777在线视频| 国产成人精品久久二区二区91| 国产精品亚洲av一区麻豆| 在线国产一区二区在线| 国产野战对白在线观看| 性色av乱码一区二区三区2| 精品一区二区三区四区五区乱码| 搡老熟女国产l中国老女人| 日韩精品青青久久久久久| 久久中文看片网| 男人的好看免费观看在线视频 | 国产三级在线视频| 久久久精品欧美日韩精品| bbb黄色大片| 精品国内亚洲2022精品成人| 久久久久久人人人人人| 国产午夜精品久久久久久| 最新在线观看一区二区三区| 黑人操中国人逼视频| 久久人人97超碰香蕉20202| 女生性感内裤真人,穿戴方法视频| 如日韩欧美国产精品一区二区三区| 级片在线观看| 亚洲一区二区三区色噜噜 | 欧美一级毛片孕妇| 97人妻天天添夜夜摸| 在线观看免费视频网站a站| 国产色视频综合| 91精品国产国语对白视频| 午夜影院日韩av| 国产精品久久久久久人妻精品电影| 国产精品一区二区在线不卡| ponron亚洲| 香蕉丝袜av| 日韩欧美在线二视频| 亚洲三区欧美一区| 91av网站免费观看| 正在播放国产对白刺激| 婷婷精品国产亚洲av在线| 免费人成视频x8x8入口观看| 日日夜夜操网爽| 成人三级黄色视频| 国产熟女午夜一区二区三区| a级毛片黄视频| 日本精品一区二区三区蜜桃| 久久久久久亚洲精品国产蜜桃av| 欧美成人午夜精品| 日本免费a在线| av天堂久久9| 亚洲 欧美 日韩 在线 免费| 亚洲中文av在线| 精品熟女少妇八av免费久了| 国产精品野战在线观看 | 午夜精品久久久久久毛片777| 黑人欧美特级aaaaaa片| 三上悠亚av全集在线观看| 日韩三级视频一区二区三区| 一进一出抽搐动态| 岛国在线观看网站| 亚洲第一av免费看| 亚洲五月婷婷丁香| 国产精品野战在线观看 | 无人区码免费观看不卡| 99久久99久久久精品蜜桃| 精品国产一区二区三区四区第35| 免费搜索国产男女视频| 男女高潮啪啪啪动态图| 夜夜看夜夜爽夜夜摸 | 亚洲五月婷婷丁香| 丰满迷人的少妇在线观看| 欧美日韩黄片免| 久久香蕉国产精品| 99久久精品国产亚洲精品| av免费在线观看网站| 1024香蕉在线观看| 99精品久久久久人妻精品| 国产成人免费无遮挡视频| 伊人久久大香线蕉亚洲五| 久久久精品欧美日韩精品| 女人被狂操c到高潮| 在线观看舔阴道视频| 亚洲av成人不卡在线观看播放网| 成人免费观看视频高清| 久久人人爽av亚洲精品天堂| 老司机靠b影院| videosex国产| 精品久久久久久成人av| 亚洲精品国产区一区二| 99国产综合亚洲精品| 黑人猛操日本美女一级片| 妹子高潮喷水视频| 露出奶头的视频| 最新在线观看一区二区三区| 变态另类成人亚洲欧美熟女 | 日日干狠狠操夜夜爽| 欧洲精品卡2卡3卡4卡5卡区| 国产视频一区二区在线看| 亚洲久久久国产精品| 别揉我奶头~嗯~啊~动态视频| 日韩大尺度精品在线看网址 | 变态另类成人亚洲欧美熟女 | 免费搜索国产男女视频| 搡老乐熟女国产| 每晚都被弄得嗷嗷叫到高潮| 国产黄a三级三级三级人| 日本一区二区免费在线视频| a级毛片黄视频| 国产欧美日韩一区二区三| 成年人黄色毛片网站| 91精品国产国语对白视频| 国产野战对白在线观看| 黄片大片在线免费观看| 免费在线观看日本一区| 美女大奶头视频| 交换朋友夫妻互换小说| 亚洲成国产人片在线观看| 精品福利观看| av在线播放免费不卡| 最近最新免费中文字幕在线| 18禁黄网站禁片午夜丰满| 黄色a级毛片大全视频| 国产精品久久电影中文字幕| 97碰自拍视频| av在线天堂中文字幕 | 久久久久久人人人人人| 91成年电影在线观看| 欧美中文日本在线观看视频| 亚洲欧美激情在线| 免费高清在线观看日韩| svipshipincom国产片| 一边摸一边做爽爽视频免费| 久久久国产精品麻豆| 久久久久久亚洲精品国产蜜桃av| 少妇 在线观看| 亚洲av成人av| 成人国语在线视频| 成人精品一区二区免费| 老汉色∧v一级毛片| 久久人人97超碰香蕉20202| 动漫黄色视频在线观看| 精品国产美女av久久久久小说| 亚洲七黄色美女视频| 日本一区二区免费在线视频| 50天的宝宝边吃奶边哭怎么回事| 老司机在亚洲福利影院| 伊人久久大香线蕉亚洲五| 国产黄a三级三级三级人| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美黄色淫秽网站| 9191精品国产免费久久| 男女午夜视频在线观看| 日韩欧美免费精品| 欧美一级毛片孕妇| 女性被躁到高潮视频| 国产精品国产av在线观看| 久久精品91无色码中文字幕| 久久午夜综合久久蜜桃| 不卡一级毛片| 欧美激情极品国产一区二区三区| 热99re8久久精品国产| 色综合婷婷激情| 国产欧美日韩一区二区三区在线| 国产一区二区在线av高清观看| av福利片在线| 在线观看免费高清a一片| ponron亚洲| 日韩三级视频一区二区三区| 久99久视频精品免费| 可以在线观看毛片的网站| 色哟哟哟哟哟哟| 亚洲男人天堂网一区| 精品免费久久久久久久清纯| 精品无人区乱码1区二区| 精品国内亚洲2022精品成人| 侵犯人妻中文字幕一二三四区| 亚洲精品一区av在线观看| 亚洲,欧美精品.| 亚洲一码二码三码区别大吗| 99久久精品国产亚洲精品| 免费观看精品视频网站| 亚洲精品成人av观看孕妇| 免费在线观看影片大全网站| 国产欧美日韩精品亚洲av| 麻豆av在线久日| 老司机午夜福利在线观看视频| 精品高清国产在线一区| 黄网站色视频无遮挡免费观看| 亚洲国产看品久久| 黄频高清免费视频| av网站在线播放免费| 不卡一级毛片| 最近最新免费中文字幕在线| 精品高清国产在线一区| 他把我摸到了高潮在线观看| 日本撒尿小便嘘嘘汇集6| 国产亚洲精品久久久久久毛片| 午夜福利在线免费观看网站| 丰满人妻熟妇乱又伦精品不卡| 亚洲全国av大片| 麻豆av在线久日| 757午夜福利合集在线观看| 精品国产乱码久久久久久男人| 成人精品一区二区免费| 精品国产国语对白av| 午夜免费观看网址| 最新在线观看一区二区三区| ponron亚洲| 欧美日韩亚洲综合一区二区三区_| 国产av一区在线观看免费| 纯流量卡能插随身wifi吗| 久久人妻熟女aⅴ| 人人澡人人妻人| 丁香欧美五月| av片东京热男人的天堂| 中文亚洲av片在线观看爽| 国产亚洲精品久久久久5区| 亚洲精品久久午夜乱码| av在线播放免费不卡| 国产国语露脸激情在线看| 九色亚洲精品在线播放| 丝袜美腿诱惑在线| 人人妻人人添人人爽欧美一区卜| 国产精品香港三级国产av潘金莲| x7x7x7水蜜桃| 丝袜在线中文字幕| 18禁裸乳无遮挡免费网站照片 | 午夜影院日韩av| 亚洲精华国产精华精| 看片在线看免费视频| 国产高清激情床上av| 亚洲人成伊人成综合网2020| 不卡一级毛片| 69av精品久久久久久| 精品高清国产在线一区| 亚洲成人国产一区在线观看| 久久久国产一区二区| 亚洲全国av大片| 国产成人影院久久av| 老司机深夜福利视频在线观看| 十八禁网站免费在线| 精品欧美一区二区三区在线| 嫩草影视91久久| 久久精品亚洲av国产电影网| 色老头精品视频在线观看| 一个人免费在线观看的高清视频| 中文字幕av电影在线播放| 午夜免费观看网址| 久久伊人香网站| 欧美日本亚洲视频在线播放| 国产单亲对白刺激| 日韩精品免费视频一区二区三区| av网站在线播放免费| 亚洲av日韩精品久久久久久密| 国产黄色免费在线视频| 亚洲 欧美 日韩 在线 免费| 国产熟女午夜一区二区三区| 在线观看www视频免费| 精品福利永久在线观看| 在线观看免费视频日本深夜| 亚洲人成电影观看| 欧美日韩精品网址| 18美女黄网站色大片免费观看| 97人妻天天添夜夜摸| 亚洲精华国产精华精| 侵犯人妻中文字幕一二三四区| 一边摸一边抽搐一进一小说| 丰满饥渴人妻一区二区三| 国产av又大| 身体一侧抽搐| 99热国产这里只有精品6| 欧美一区二区精品小视频在线| 在线观看免费午夜福利视频| 亚洲精品一区av在线观看| 两人在一起打扑克的视频| 村上凉子中文字幕在线| 欧美日韩亚洲国产一区二区在线观看| 日韩欧美一区视频在线观看| 精品第一国产精品| 成年人免费黄色播放视频| av国产精品久久久久影院| 美女高潮到喷水免费观看| 50天的宝宝边吃奶边哭怎么回事| 在线天堂中文资源库| 一区二区三区国产精品乱码| 中文字幕人妻丝袜制服| 欧美日韩av久久| 18美女黄网站色大片免费观看| 亚洲成人免费电影在线观看| 在线观看日韩欧美| av中文乱码字幕在线| 两人在一起打扑克的视频| 美女大奶头视频| 精品久久久久久电影网| 亚洲精品在线美女| 国产欧美日韩一区二区三| 91精品三级在线观看| 精品国产一区二区久久| 侵犯人妻中文字幕一二三四区| 黄色视频,在线免费观看| 国产精品久久电影中文字幕| 99久久久亚洲精品蜜臀av| 国产亚洲精品久久久久5区| 丁香六月欧美| 亚洲成人国产一区在线观看| 欧美黄色淫秽网站| 麻豆成人av在线观看| av在线天堂中文字幕 | 欧美国产精品va在线观看不卡| 在线天堂中文资源库| av天堂久久9| 99re在线观看精品视频| 国产精品美女特级片免费视频播放器 | 黄色视频,在线免费观看| 黄色毛片三级朝国网站| 亚洲,欧美精品.| 国产人伦9x9x在线观看| av欧美777| 久久草成人影院| 人人妻,人人澡人人爽秒播| 别揉我奶头~嗯~啊~动态视频| 热re99久久精品国产66热6| 黄片大片在线免费观看| 一边摸一边抽搐一进一出视频| а√天堂www在线а√下载| 欧美国产精品va在线观看不卡| 国产精品自产拍在线观看55亚洲| 国产激情欧美一区二区| 在线天堂中文资源库| 在线国产一区二区在线| 狠狠狠狠99中文字幕| 女性被躁到高潮视频| 久久久久国内视频| 天天躁夜夜躁狠狠躁躁| 中出人妻视频一区二区| 精品乱码久久久久久99久播| 看片在线看免费视频| 日本免费a在线| 国产精品免费视频内射| 欧美日韩福利视频一区二区| 国产在线观看jvid| 国内久久婷婷六月综合欲色啪| 免费在线观看完整版高清| 可以在线观看毛片的网站| 国产有黄有色有爽视频| 精品少妇一区二区三区视频日本电影| 99在线视频只有这里精品首页| av天堂在线播放| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品久久久人人做人人爽| 淫妇啪啪啪对白视频| 久久久久久大精品| 久久久久国产精品人妻aⅴ院| 这个男人来自地球电影免费观看| 一级a爱视频在线免费观看| 久9热在线精品视频| 91成人精品电影| 久久伊人香网站| 日韩三级视频一区二区三区| 久久精品亚洲熟妇少妇任你| 成人手机av| 成人国产一区最新在线观看| 成人三级做爰电影| 91精品三级在线观看| av电影中文网址| 免费看十八禁软件| 日本免费一区二区三区高清不卡 | www日本在线高清视频| 桃红色精品国产亚洲av| 国产精品影院久久| 亚洲国产中文字幕在线视频| 日本精品一区二区三区蜜桃| 无限看片的www在线观看| 又大又爽又粗| 精品人妻1区二区| 国产激情久久老熟女| 十八禁网站免费在线| 亚洲一区二区三区不卡视频| 午夜福利在线观看吧| 欧美日韩亚洲高清精品| 成人三级做爰电影| 欧美日本中文国产一区发布| 在线十欧美十亚洲十日本专区| 国产精品久久视频播放| 亚洲av成人不卡在线观看播放网| √禁漫天堂资源中文www| 亚洲熟妇中文字幕五十中出 | 国产成人影院久久av| 久久欧美精品欧美久久欧美| 多毛熟女@视频| 老汉色av国产亚洲站长工具| 久久精品亚洲熟妇少妇任你| 超色免费av| 亚洲av成人不卡在线观看播放网| 黑人猛操日本美女一级片| 国产亚洲欧美在线一区二区| 精品福利观看| 久9热在线精品视频| 亚洲自拍偷在线| 亚洲欧洲精品一区二区精品久久久| 琪琪午夜伦伦电影理论片6080| 又黄又爽又免费观看的视频| 日韩av在线大香蕉| 日本免费一区二区三区高清不卡 | 欧美精品亚洲一区二区| 三级毛片av免费| 久热这里只有精品99| 亚洲精品av麻豆狂野| cao死你这个sao货| 99久久久亚洲精品蜜臀av| 精品日产1卡2卡| 中文字幕人妻丝袜一区二区| 欧美国产精品va在线观看不卡| 国产精品偷伦视频观看了| 村上凉子中文字幕在线| 精品久久久久久久久久免费视频 | 午夜日韩欧美国产| 亚洲精华国产精华精| 日韩欧美国产一区二区入口| 伦理电影免费视频| videosex国产| 色老头精品视频在线观看| 两性夫妻黄色片| 亚洲第一av免费看| 亚洲精品久久午夜乱码| 欧美在线一区亚洲| 男人舔女人的私密视频| 黄色丝袜av网址大全| 成在线人永久免费视频| 激情视频va一区二区三区| 女人精品久久久久毛片| 天堂动漫精品| 欧美日韩瑟瑟在线播放| 国产极品粉嫩免费观看在线| 午夜福利欧美成人| 日日爽夜夜爽网站| 国产伦人伦偷精品视频| 丰满饥渴人妻一区二区三| 久久草成人影院| 成年人免费黄色播放视频| 少妇的丰满在线观看| 欧美黄色片欧美黄色片| 免费日韩欧美在线观看| 99热国产这里只有精品6| 国产蜜桃级精品一区二区三区| 99久久99久久久精品蜜桃| 超色免费av| 国产精品久久久av美女十八| 免费观看精品视频网站| 母亲3免费完整高清在线观看| 久久欧美精品欧美久久欧美| 久久影院123| 日本a在线网址| 每晚都被弄得嗷嗷叫到高潮| 亚洲国产中文字幕在线视频| 免费高清在线观看日韩| 亚洲精品在线美女| 欧美日韩亚洲国产一区二区在线观看| videosex国产| 精品国产亚洲在线| 久久草成人影院| 亚洲中文av在线| 久久香蕉精品热| 中出人妻视频一区二区| 亚洲中文av在线| 午夜两性在线视频| 成年女人毛片免费观看观看9| 日本黄色日本黄色录像| 国产乱人伦免费视频| 99热国产这里只有精品6| 欧美黑人精品巨大| 午夜91福利影院| 青草久久国产| 男女下面进入的视频免费午夜 | 成人国语在线视频| 日韩三级视频一区二区三区| 变态另类成人亚洲欧美熟女 | 国产精品国产高清国产av| 国产麻豆69| 亚洲,欧美精品.| 一边摸一边抽搐一进一出视频| 精品久久久久久电影网| 国产精品一区二区三区四区久久 | 天堂影院成人在线观看| 久久久精品欧美日韩精品| 精品国产超薄肉色丝袜足j| 自拍欧美九色日韩亚洲蝌蚪91| 露出奶头的视频| 身体一侧抽搐| 成人18禁在线播放| 天堂中文最新版在线下载| 亚洲专区中文字幕在线| 久久人妻av系列| 一个人观看的视频www高清免费观看 |