• <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ì)比及分析
    欧美另类亚洲清纯唯美| 欧美黄色片欧美黄色片| 欧美日韩成人在线一区二区| 国产精品一区二区在线不卡| 热99国产精品久久久久久7| 999久久久国产精品视频| 深夜精品福利| 两性夫妻黄色片| av网站在线播放免费| 日韩免费高清中文字幕av| 自线自在国产av| 深夜精品福利| 久久久国产成人免费| 久久热在线av| 欧美亚洲日本最大视频资源| 久久久久精品人妻al黑| 午夜精品在线福利| 黄色毛片三级朝国网站| 亚洲精品在线观看二区| 久久精品人人爽人人爽视色| 性少妇av在线| 亚洲欧美一区二区三区久久| 亚洲成人免费av在线播放| 久久久国产精品麻豆| 欧美人与性动交α欧美精品济南到| 最近最新中文字幕大全电影3 | 黑丝袜美女国产一区| 欧美人与性动交α欧美精品济南到| 黄片播放在线免费| 亚洲精品国产区一区二| 国产一区二区三区综合在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 一级毛片高清免费大全| 国产欧美日韩一区二区三区在线| 91成年电影在线观看| 成人黄色视频免费在线看| 伦理电影免费视频| 动漫黄色视频在线观看| 在线天堂中文资源库| 亚洲色图 男人天堂 中文字幕| 成人18禁高潮啪啪吃奶动态图| 精品福利观看| 欧美激情 高清一区二区三区| 老鸭窝网址在线观看| 美女扒开内裤让男人捅视频| 麻豆成人av在线观看| 久热这里只有精品99| av超薄肉色丝袜交足视频| 法律面前人人平等表现在哪些方面| 欧美精品啪啪一区二区三区| 欧美午夜高清在线| 久久九九热精品免费| 亚洲在线自拍视频| 成人av一区二区三区在线看| 久久香蕉国产精品| 99热网站在线观看| 色老头精品视频在线观看| 巨乳人妻的诱惑在线观看| 黄频高清免费视频| 成人三级做爰电影| 久久国产精品影院| 天堂√8在线中文| 亚洲熟女精品中文字幕| 亚洲av片天天在线观看| 村上凉子中文字幕在线| 男女下面插进去视频免费观看| 色综合欧美亚洲国产小说| 99国产精品一区二区三区| 在线观看一区二区三区激情| 久久精品91无色码中文字幕| 欧美久久黑人一区二区| 嫩草影视91久久| 午夜福利免费观看在线| 老司机午夜福利在线观看视频| 免费日韩欧美在线观看| 久久精品国产清高在天天线| 在线观看日韩欧美| 国产亚洲精品一区二区www | 国产aⅴ精品一区二区三区波| 精品少妇一区二区三区视频日本电影| 欧美黑人欧美精品刺激| 99热国产这里只有精品6| 99国产精品99久久久久| 丁香欧美五月| 99riav亚洲国产免费| 91大片在线观看| 亚洲精品在线美女| 日本黄色日本黄色录像| 久久久久久亚洲精品国产蜜桃av| 老熟妇乱子伦视频在线观看| 高清视频免费观看一区二区| 精品久久久久久电影网| 国产精品99久久99久久久不卡| 国产欧美日韩一区二区三| 欧美日韩成人在线一区二区| 色播在线永久视频| av欧美777| 午夜福利在线观看吧| av天堂久久9| 亚洲精品中文字幕在线视频| 人人妻人人澡人人爽人人夜夜| 日本vs欧美在线观看视频| 成人永久免费在线观看视频| 国产免费av片在线观看野外av| 18在线观看网站| 久久人人爽av亚洲精品天堂| 婷婷精品国产亚洲av在线 | 国产精品免费视频内射| 欧美激情极品国产一区二区三区| 中出人妻视频一区二区| 成年版毛片免费区| 黄色女人牲交| 无人区码免费观看不卡| 激情在线观看视频在线高清 | 日日夜夜操网爽| 韩国精品一区二区三区| 久久久久久久精品吃奶| 国产aⅴ精品一区二区三区波| 老汉色∧v一级毛片| 国产高清视频在线播放一区| 国产一区二区激情短视频| 精品无人区乱码1区二区| 香蕉久久夜色| 日韩欧美三级三区| 99热国产这里只有精品6| 香蕉久久夜色| 日韩欧美一区视频在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 又黄又爽又免费观看的视频| 91精品三级在线观看| 日本wwww免费看| 国产在视频线精品| 久久精品国产亚洲av香蕉五月 | 村上凉子中文字幕在线| 老鸭窝网址在线观看| 人人澡人人妻人| 男人舔女人的私密视频| 国产成人精品久久二区二区91| 在线播放国产精品三级| 亚洲全国av大片| 中文字幕另类日韩欧美亚洲嫩草| 国产精品永久免费网站| 伊人久久大香线蕉亚洲五| 国产av精品麻豆| 人妻久久中文字幕网| 女同久久另类99精品国产91| 国产又色又爽无遮挡免费看| 亚洲av美国av| 亚洲精品一卡2卡三卡4卡5卡| 黄色毛片三级朝国网站| netflix在线观看网站| 身体一侧抽搐| 人人妻人人爽人人添夜夜欢视频| 高清毛片免费观看视频网站 | 久久香蕉国产精品| 国产黄色免费在线视频| 很黄的视频免费| 久久久国产一区二区| 国产在视频线精品| 91麻豆精品激情在线观看国产 | 久久国产亚洲av麻豆专区| 女人高潮潮喷娇喘18禁视频| 国产精华一区二区三区| 在线十欧美十亚洲十日本专区| 久久精品国产清高在天天线| 可以免费在线观看a视频的电影网站| 黄色成人免费大全| 久久国产精品大桥未久av| 精品午夜福利视频在线观看一区| 超碰97精品在线观看| 91成年电影在线观看| www.精华液| 免费在线观看影片大全网站| 老鸭窝网址在线观看| 日日夜夜操网爽| 国产精品久久久av美女十八| 成人手机av| 王馨瑶露胸无遮挡在线观看| 亚洲精品一二三| 99香蕉大伊视频| 91九色精品人成在线观看| 午夜福利乱码中文字幕| x7x7x7水蜜桃| 国产精品1区2区在线观看. | 欧美大码av| 丰满饥渴人妻一区二区三| 一个人免费在线观看的高清视频| 777米奇影视久久| 丰满饥渴人妻一区二区三| 国产av精品麻豆| 亚洲欧洲精品一区二区精品久久久| 男人舔女人的私密视频| 黄片大片在线免费观看| 好看av亚洲va欧美ⅴa在| 国产av一区二区精品久久| 欧美日韩福利视频一区二区| 一级a爱视频在线免费观看| 黄色视频不卡| 正在播放国产对白刺激| 涩涩av久久男人的天堂| 丝袜美腿诱惑在线| 国产成人欧美在线观看 | 国产精品香港三级国产av潘金莲| 欧美黑人欧美精品刺激| 色94色欧美一区二区| 18禁美女被吸乳视频| 一级,二级,三级黄色视频| 久久午夜亚洲精品久久| √禁漫天堂资源中文www| 亚洲欧美日韩另类电影网站| 国产一卡二卡三卡精品| 国产色视频综合| 搡老乐熟女国产| 精品久久久精品久久久| www.自偷自拍.com| 久久久久久久午夜电影 | 国产精品永久免费网站| 国产亚洲精品久久久久久毛片 | 久久青草综合色| 老熟妇仑乱视频hdxx| 国产成人精品在线电影| e午夜精品久久久久久久| 两性午夜刺激爽爽歪歪视频在线观看 | 国产激情久久老熟女| 亚洲av成人一区二区三| 男人舔女人的私密视频| 1024视频免费在线观看| 在线十欧美十亚洲十日本专区| 日本精品一区二区三区蜜桃| 91麻豆av在线| 午夜精品国产一区二区电影| 在线观看免费日韩欧美大片| 手机成人av网站| 美女高潮到喷水免费观看| 日韩欧美在线二视频 | 一级片'在线观看视频| 女人被躁到高潮嗷嗷叫费观| 啪啪无遮挡十八禁网站| 一本综合久久免费| 91成年电影在线观看| 91大片在线观看| 国产精品久久久av美女十八| 日本撒尿小便嘘嘘汇集6| 日韩中文字幕欧美一区二区| 老鸭窝网址在线观看| 午夜91福利影院| 欧美老熟妇乱子伦牲交| 女人被狂操c到高潮| 999精品在线视频| cao死你这个sao货| 亚洲专区字幕在线| 亚洲av电影在线进入| 欧美大码av| 国产又爽黄色视频| 极品人妻少妇av视频| 国产成人免费无遮挡视频| 国产在线观看jvid| 色婷婷av一区二区三区视频| 三上悠亚av全集在线观看| 麻豆av在线久日| 亚洲自偷自拍图片 自拍| 久久香蕉国产精品| 一级毛片女人18水好多| 久久国产精品男人的天堂亚洲| 99香蕉大伊视频| 一级黄色大片毛片| 夜夜爽天天搞| 亚洲一区中文字幕在线| 叶爱在线成人免费视频播放| 丝袜美足系列| 亚洲熟妇中文字幕五十中出 | 下体分泌物呈黄色| 国产亚洲精品久久久久5区| 啦啦啦在线免费观看视频4| 欧美午夜高清在线| 免费av中文字幕在线| 免费看a级黄色片| 后天国语完整版免费观看| 97人妻天天添夜夜摸| 精品国产乱码久久久久久男人| 一级毛片精品| 啦啦啦 在线观看视频| av天堂久久9| 国产又色又爽无遮挡免费看| 精品久久蜜臀av无| 久久国产乱子伦精品免费另类| av福利片在线| 久久久久久人人人人人| 如日韩欧美国产精品一区二区三区| 老司机在亚洲福利影院| 狠狠狠狠99中文字幕| 黑人操中国人逼视频| 日韩三级视频一区二区三区| 精品国内亚洲2022精品成人 | 老司机福利观看| 女人高潮潮喷娇喘18禁视频| 精品一区二区三区四区五区乱码| 欧美+亚洲+日韩+国产| 香蕉国产在线看| 一级a爱视频在线免费观看| 91av网站免费观看| 天天躁日日躁夜夜躁夜夜| 久久香蕉国产精品| 国产1区2区3区精品| 宅男免费午夜| 欧美精品av麻豆av| 精品一区二区三卡| 一级片'在线观看视频| 色在线成人网| 亚洲成国产人片在线观看| 亚洲全国av大片| 亚洲熟女精品中文字幕| 久久香蕉精品热| 亚洲精品中文字幕一二三四区| 丝袜美足系列| 国产男女超爽视频在线观看| 亚洲av片天天在线观看| 国产一区二区激情短视频| 91老司机精品| 亚洲久久久国产精品| 最新美女视频免费是黄的| 人妻 亚洲 视频| 中出人妻视频一区二区| 99re在线观看精品视频| 欧美成人免费av一区二区三区 | 大香蕉久久成人网| 又紧又爽又黄一区二区| 欧美精品人与动牲交sv欧美| 在线十欧美十亚洲十日本专区| 俄罗斯特黄特色一大片| av电影中文网址| 91成年电影在线观看| 黄色毛片三级朝国网站| 午夜福利免费观看在线| 看片在线看免费视频| 99国产精品免费福利视频| 成人18禁高潮啪啪吃奶动态图| 80岁老熟妇乱子伦牲交| 黑人猛操日本美女一级片| 美女午夜性视频免费| 国产区一区二久久| 色在线成人网| 婷婷丁香在线五月| 亚洲片人在线观看| 男女高潮啪啪啪动态图| 亚洲成人国产一区在线观看| 精品久久久久久久久久免费视频 | 午夜福利免费观看在线| 91在线观看av| 免费女性裸体啪啪无遮挡网站| av电影中文网址| 大陆偷拍与自拍| 91在线观看av| 男女之事视频高清在线观看| 女性被躁到高潮视频| 亚洲七黄色美女视频| 又黄又粗又硬又大视频| 大型av网站在线播放| 亚洲国产精品合色在线| 久久亚洲真实| 午夜两性在线视频| 国产av又大| 热re99久久精品国产66热6| a在线观看视频网站| 麻豆成人av在线观看| 久久精品成人免费网站| 精品国产亚洲在线| 99香蕉大伊视频| 国产精品综合久久久久久久免费 | 久久国产精品大桥未久av| 色综合婷婷激情| 成人亚洲精品一区在线观看| 精品卡一卡二卡四卡免费| 精品久久久久久久毛片微露脸| 看黄色毛片网站| 露出奶头的视频| 王馨瑶露胸无遮挡在线观看| 国产亚洲精品久久久久5区| av不卡在线播放| 狠狠婷婷综合久久久久久88av| 欧美人与性动交α欧美软件| 人成视频在线观看免费观看| 黑人巨大精品欧美一区二区mp4| 最新的欧美精品一区二区| 亚洲三区欧美一区| 久久久国产成人免费| 三上悠亚av全集在线观看| 99久久国产精品久久久| 精品久久久久久,| 欧美一级毛片孕妇| av视频免费观看在线观看| 伦理电影免费视频| 91精品三级在线观看| 国产激情欧美一区二区| 亚洲av电影在线进入| 咕卡用的链子| 99国产精品一区二区蜜桃av | 免费在线观看黄色视频的| 午夜视频精品福利| 国产精品秋霞免费鲁丝片| 国产男靠女视频免费网站| 国产单亲对白刺激| 一区福利在线观看| 成人18禁高潮啪啪吃奶动态图| 人人妻人人爽人人添夜夜欢视频| 国产成人免费无遮挡视频| 露出奶头的视频| 咕卡用的链子| 精品福利永久在线观看| 国产亚洲一区二区精品| 操出白浆在线播放| 无遮挡黄片免费观看| 99久久综合精品五月天人人| 另类亚洲欧美激情| 搡老乐熟女国产| 亚洲熟妇中文字幕五十中出 | 99国产极品粉嫩在线观看| 首页视频小说图片口味搜索| 国产主播在线观看一区二区| 欧美大码av| 老汉色av国产亚洲站长工具| 亚洲av成人一区二区三| 满18在线观看网站| 亚洲中文字幕日韩| 女警被强在线播放| 在线视频色国产色| 亚洲精品国产色婷婷电影| 欧美黑人欧美精品刺激| 亚洲中文日韩欧美视频| 亚洲av第一区精品v没综合| 欧美精品亚洲一区二区| 国产成人影院久久av| 悠悠久久av| 免费在线观看亚洲国产| 欧美精品高潮呻吟av久久| 国产激情久久老熟女| a级片在线免费高清观看视频| 亚洲在线自拍视频| 日韩视频一区二区在线观看| 不卡av一区二区三区| 老鸭窝网址在线观看| 日本黄色视频三级网站网址 | 免费少妇av软件| 午夜成年电影在线免费观看| 亚洲av美国av| 不卡av一区二区三区| 亚洲av欧美aⅴ国产| 黄网站色视频无遮挡免费观看| 国产熟女午夜一区二区三区| 涩涩av久久男人的天堂| 色尼玛亚洲综合影院| 1024香蕉在线观看| 中文字幕av电影在线播放| 国内毛片毛片毛片毛片毛片| 在线观看66精品国产| 亚洲熟女精品中文字幕| 91成人精品电影| 亚洲综合色网址| 国产蜜桃级精品一区二区三区 | 欧美乱码精品一区二区三区| 日韩欧美免费精品| 久久久久国产精品人妻aⅴ院 | 韩国精品一区二区三区| 久久久国产一区二区| 午夜福利,免费看| 亚洲成人手机| 91av网站免费观看| 国产精品亚洲一级av第二区| 国产国语露脸激情在线看| 日韩成人在线观看一区二区三区| www.熟女人妻精品国产| 国产精品.久久久| 一本大道久久a久久精品| 国产精品久久视频播放| 999久久久国产精品视频| 狂野欧美激情性xxxx| 午夜免费观看网址| 亚洲成人手机| 亚洲国产欧美一区二区综合| 纯流量卡能插随身wifi吗| 久久精品国产99精品国产亚洲性色 | 久久久久久久久免费视频了| 国产男女超爽视频在线观看| 一二三四在线观看免费中文在| 午夜精品久久久久久毛片777| 亚洲成人免费av在线播放| 国内久久婷婷六月综合欲色啪| 国产精品二区激情视频| 日韩免费高清中文字幕av| 在线观看午夜福利视频| 国产深夜福利视频在线观看| 夫妻午夜视频| 国产精品.久久久| 女性被躁到高潮视频| 中文欧美无线码| 国产精品影院久久| 国内久久婷婷六月综合欲色啪| 亚洲精品粉嫩美女一区| 亚洲一区中文字幕在线| 一级作爱视频免费观看| 国产精品久久久av美女十八| 久久久久久久久免费视频了| 亚洲九九香蕉| 一区二区三区激情视频| 久久久久久亚洲精品国产蜜桃av| a级片在线免费高清观看视频| 色播在线永久视频| 亚洲精华国产精华精| 一进一出抽搐gif免费好疼 | 亚洲专区国产一区二区| 三级毛片av免费| 国产不卡一卡二| 精品国内亚洲2022精品成人 | 亚洲av成人一区二区三| 亚洲精华国产精华精| 日韩免费av在线播放| 高清av免费在线| 亚洲熟妇熟女久久| 美女国产高潮福利片在线看| 在线观看日韩欧美| av欧美777| 日日爽夜夜爽网站| 人成视频在线观看免费观看| 日韩免费av在线播放| 免费久久久久久久精品成人欧美视频| 欧美激情极品国产一区二区三区| 曰老女人黄片| 精品熟女少妇八av免费久了| 国产又色又爽无遮挡免费看| 黄色毛片三级朝国网站| 一级片免费观看大全| 亚洲专区中文字幕在线| svipshipincom国产片| 国产免费现黄频在线看| 又紧又爽又黄一区二区| 99热网站在线观看| 99热只有精品国产| 成人永久免费在线观看视频| 成年人黄色毛片网站| 天天影视国产精品| xxxhd国产人妻xxx| 亚洲一码二码三码区别大吗| 日韩免费av在线播放| 国产成人啪精品午夜网站| av有码第一页| 一级黄色大片毛片| 国产成人免费观看mmmm| 99精品欧美一区二区三区四区| 欧美人与性动交α欧美软件| 麻豆国产av国片精品| 国产一区在线观看成人免费| 久久 成人 亚洲| 一级片免费观看大全| 黑人巨大精品欧美一区二区蜜桃| 黑人猛操日本美女一级片| 免费高清在线观看日韩| 建设人人有责人人尽责人人享有的| 女人精品久久久久毛片| 亚洲久久久国产精品| 欧美激情久久久久久爽电影 | 国产av又大| 日韩 欧美 亚洲 中文字幕| 国产成人免费无遮挡视频| 亚洲精华国产精华精| 婷婷精品国产亚洲av在线 | 黄色视频,在线免费观看| 国产亚洲一区二区精品| 婷婷丁香在线五月| 日韩熟女老妇一区二区性免费视频| 97人妻天天添夜夜摸| 亚洲欧美激情综合另类| 日韩欧美三级三区| www日本在线高清视频| 国产极品粉嫩免费观看在线| av视频免费观看在线观看| 91九色精品人成在线观看| 久久九九热精品免费| av在线播放免费不卡| 国产成人免费无遮挡视频| 免费看十八禁软件| 免费黄频网站在线观看国产| 国产99久久九九免费精品| 国产精品 国内视频| 男女午夜视频在线观看| 精品少妇久久久久久888优播| 欧美精品人与动牲交sv欧美| 日本五十路高清| 50天的宝宝边吃奶边哭怎么回事| 精品福利永久在线观看| av中文乱码字幕在线| 高潮久久久久久久久久久不卡| 国产成人精品久久二区二区91| 黑人巨大精品欧美一区二区蜜桃| 久久国产精品人妻蜜桃| 国产成人精品久久二区二区91| 日韩制服丝袜自拍偷拍| 丰满迷人的少妇在线观看| 亚洲一码二码三码区别大吗| 乱人伦中国视频| 精品一区二区三区av网在线观看| 欧美精品人与动牲交sv欧美| 精品亚洲成a人片在线观看| 精品国产乱子伦一区二区三区| 一级片免费观看大全| 亚洲美女黄片视频| 中文字幕av电影在线播放| a在线观看视频网站| 99re在线观看精品视频| 王馨瑶露胸无遮挡在线观看| 后天国语完整版免费观看| 午夜福利一区二区在线看| 亚洲自偷自拍图片 自拍| 视频区图区小说| 人人妻人人爽人人添夜夜欢视频| 欧美日韩精品网址|