亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        快速均值回歸隨機(jī)波動(dòng)率模型參數(shù)估計(jì)及應(yīng)用

        2020-10-24 02:16:58李蓬實(shí)楊建輝
        運(yùn)籌與管理 2020年2期
        關(guān)鍵詞:標(biāo)的模型

        李蓬實(shí), 楊建輝, 林 焰

        (1.東莞理工學(xué)院 經(jīng)濟(jì)與管理學(xué)院,廣東 東莞 523106; 2.華南理工大學(xué) 工商管理學(xué)院,廣東 廣州 510640; 3.北京大學(xué) 匯豐商學(xué)院,深圳 100871)

        0 引言

        傳統(tǒng)的期權(quán)定價(jià)模型假設(shè)波動(dòng)率是常數(shù), 根據(jù)該假設(shè), 不同執(zhí)行價(jià)格的同種標(biāo)的資產(chǎn)的看漲期權(quán)應(yīng)具有相同的波動(dòng)率。 然而實(shí)證研究表明,從金融市場(chǎng)上的看漲期權(quán)價(jià)格推算出的波動(dòng)率在平值期權(quán)(at the money)附近的波動(dòng)率相對(duì)較低, 但隨著期權(quán)的虛值程度(out of the money)或?qū)嵵党潭?in the money)增加波動(dòng)率也逐漸增大。 這種現(xiàn)象被稱為隱含波動(dòng)率的“微笑曲線”(smile)現(xiàn)象。 隱含波動(dòng)率的“微笑曲線”現(xiàn)象說明了在給定到期期限條件下, 隱含波動(dòng)率與期權(quán)的在值程度或執(zhí)行價(jià)格之間的關(guān)系并非固定不變的, 因此傳統(tǒng)的期權(quán)定價(jià)模型中波動(dòng)率為常數(shù)的假設(shè)并不符合現(xiàn)實(shí)。 除此之外, 隱含波動(dòng)率“微笑曲線”還可能是由于標(biāo)的資產(chǎn)價(jià)格的杠桿效應(yīng)(標(biāo)的資產(chǎn)價(jià)格與波動(dòng)率之間的負(fù)相關(guān)關(guān)系)與標(biāo)的資產(chǎn)價(jià)格的跳躍現(xiàn)象引起的, 然而傳統(tǒng)的期權(quán)定價(jià)模型并沒有考慮杠桿效應(yīng)與資產(chǎn)價(jià)格的跳躍現(xiàn)象的存在。

        為了解決傳統(tǒng)期權(quán)模型中存在的問題, 學(xué)者們提出了各種方法來修正傳統(tǒng)期權(quán)定價(jià)模型的缺陷。 其中兩種主要的模型分別是跳躍擴(kuò)散模型和隨機(jī)波動(dòng)率模型。 跳躍擴(kuò)散模型比較適合解決到期期限較短的期權(quán)定價(jià)問題, 而隨機(jī)波動(dòng)率模型比較適合到期期限較長(zhǎng)的期權(quán)定價(jià)問題。 Bakshi等[1]證明了跳躍擴(kuò)散模型和隨機(jī)波動(dòng)率模型能夠改進(jìn)期權(quán)定價(jià)誤差,同時(shí)也發(fā)現(xiàn)在考慮隨機(jī)波動(dòng)率之后引入隨機(jī)跳躍只能夠有限地減少期權(quán)定價(jià)誤差。 因此, 本文主要研究隨機(jī)波動(dòng)率模型。

        大量實(shí)證研究表明金融資產(chǎn)具有隨機(jī)波動(dòng)率的特點(diǎn)。 Chernov等[2]發(fā)現(xiàn)金融市場(chǎng)中股票價(jià)格的收益率具有明顯的隨機(jī)波動(dòng)率性質(zhì)。 劉志東等[3]研究了上證50股票指數(shù)及其成分股的高頻數(shù)據(jù)后認(rèn)為, 中國(guó)股市中約43%的風(fēng)險(xiǎn)來源于資產(chǎn)收益過程的隨機(jī)波動(dòng)風(fēng)險(xiǎn)。 通過對(duì)中國(guó)股市的研究,吳鑫育等[4]認(rèn)為隨機(jī)波動(dòng)率模型能夠更好地增加模型擬合度。 研究證明標(biāo)的資產(chǎn)的波動(dòng)率不僅具有隨機(jī)的性質(zhì)同時(shí)具有均值回歸的特點(diǎn),同時(shí)吳鑫育等[5]運(yùn)用雙杠桿門限隨機(jī)波動(dòng)率模型對(duì)中國(guó)股市進(jìn)行實(shí)證研究發(fā)現(xiàn),股市具有強(qiáng)波動(dòng)率持續(xù)性以及顯著的杠桿效應(yīng)。 Liu等[6]通過對(duì)臺(tái)灣金融市場(chǎng)數(shù)據(jù)的研究,發(fā)現(xiàn)外匯和股指市場(chǎng)中都存在均值回歸的性質(zhì)。 在實(shí)證研究的支撐下, 一般隨機(jī)波動(dòng)率模型假設(shè)波動(dòng)率是由一個(gè)具有均值回歸性質(zhì)的擴(kuò)散過程驅(qū)動(dòng)。 Hull等[7]首先研究了具有隨機(jī)波動(dòng)率的金融資產(chǎn)的歐式期權(quán)定價(jià)問題。 Heston[8]假設(shè)隨機(jī)方差(stochastic variance)服從CIR過程, 通過特征函數(shù)的方法推導(dǎo)出歐式看漲和看跌期權(quán)的解析解。

        Fouque等[9]利用S&P500的數(shù)據(jù)證實(shí)了股指中存在的波動(dòng)率聚集(volatility clustering)和波動(dòng)率的快速均值(fast mean-reverting)回歸現(xiàn)象, 在此基礎(chǔ)上他們提出了一種基于快速均值回歸的隨機(jī)波動(dòng)率模型, 該模型在衍生品定價(jià)方面得到廣泛的應(yīng)用。 Fouque等[10]在隨機(jī)波動(dòng)率是快速均值回歸的擴(kuò)散過程假設(shè)下, 利用奇異攝動(dòng)分析方法推導(dǎo)出歐式看漲期權(quán)價(jià)格的近似解析解。 Wong等[11]研究了快速均值回歸以及常數(shù)彈性波動(dòng)率假設(shè)下渦輪期權(quán)(Turbo warrant)的定價(jià)問題。 通過引入違約風(fēng)險(xiǎn)中的隨機(jī)波動(dòng)率效應(yīng)。 Fouque等[12]研究了可違約債券(defaultable bond)的定價(jià)問題, 并推導(dǎo)出近似解析解。 Yang等[13]在隨機(jī)波動(dòng)率框架下得到脆弱期權(quán)(vulnerable option)的近似解公式, 并研究了隨機(jī)波動(dòng)率對(duì)脆弱期權(quán)價(jià)格的影響。 Kim等[14]考慮了帶有隨機(jī)利率的隨機(jī)波動(dòng)率模型, 并證明了在較短的到期期限內(nèi)新模型的對(duì)歐式期權(quán)定價(jià)效果要優(yōu)于常數(shù)利率模型。

        在上述研究成果的基礎(chǔ)上, 本文研究基于快速均值回歸隨機(jī)波動(dòng)率的雙限期權(quán)(collar option)定價(jià)問題。 雙限期權(quán)在到期時(shí)間T的收益可以表示為min(max(ST,K1),K2),其中K2>K1>0是期權(quán)的執(zhí)行價(jià)格;ST是標(biāo)的資產(chǎn)在T時(shí)刻的價(jià)格。 雙限期權(quán)可幫助投資者鎖定從投資的標(biāo)的資產(chǎn)中賺得的利潤(rùn), 還可以作為一種對(duì)沖策略使投資者避免標(biāo)的資產(chǎn)價(jià)格劇烈變化所引起的投資風(fēng)險(xiǎn)。 除了可以對(duì)沖金融產(chǎn)品的風(fēng)險(xiǎn), 雙限期權(quán)還可以用于項(xiàng)目中的收益管理。 Yim等[15]提出將雙限期權(quán)用于項(xiàng)目的物料采購(gòu), 用雙限期權(quán)控制建筑項(xiàng)目中原材料的價(jià)格波動(dòng)。

        本文基于快速均值回歸波動(dòng)率的假設(shè), 同時(shí)考慮了標(biāo)的資產(chǎn)價(jià)格過程和驅(qū)動(dòng)波動(dòng)率擴(kuò)散過程之間的相關(guān)性, 研究雙限期權(quán)的定價(jià)問題。 由于考慮了波動(dòng)率的隨機(jī)性質(zhì), 期權(quán)價(jià)格難以直接獲得, 因此本文采用Fouque等[10]的方法研究快速均值回歸波動(dòng)率框架下的期權(quán)定價(jià)。 該方法優(yōu)點(diǎn)在于, 首先可以有效地減少隨機(jī)波動(dòng)率模型中的參數(shù);其次可以通過較為簡(jiǎn)便的方法對(duì)簡(jiǎn)化后的模型參數(shù)進(jìn)行估計(jì);最后可以得到期權(quán)價(jià)格的近似解析解表達(dá)式。 本文的主要理論貢獻(xiàn)在于:首先在快速均值回歸隨機(jī)波動(dòng)率框架下研究雙限期權(quán), 由于雙限期權(quán)在資產(chǎn)風(fēng)險(xiǎn)管理中的作用, 研究雙限期權(quán)的定價(jià)問題具有一定現(xiàn)實(shí)意義; 其次, 對(duì)快速均值回歸隨機(jī)波動(dòng)率模型中的參數(shù)進(jìn)行了估計(jì), 這些參數(shù)可以進(jìn)一步應(yīng)用到其它期權(quán)的定價(jià)研究中。

        本文共分為五個(gè)部分。 第二部分介紹快速均值回歸隨機(jī)波動(dòng)率模型;第三部分給出了基于快速均值回歸隨機(jī)波動(dòng)率模型的參數(shù)估計(jì)公式;第四部分用金融市場(chǎng)數(shù)據(jù)估計(jì)了模型參數(shù), 推導(dǎo)了考慮隨機(jī)波動(dòng)率的雙限期權(quán)定價(jià)公式, 同時(shí)對(duì)雙限期權(quán)價(jià)格進(jìn)行了數(shù)值分析; 第五部分對(duì)本文進(jìn)行總結(jié)。

        1 模型

        記標(biāo)的資產(chǎn)在t時(shí)刻的價(jià)格為St,并假設(shè)標(biāo)的資產(chǎn)滿足以下隨機(jī)微分方程

        dSt=μStdt+f(Yt)StdWt

        (1)

        其中μ是標(biāo)的資產(chǎn)的期望收益率,Wt是標(biāo)準(zhǔn)布朗運(yùn)動(dòng),f(y)是大于零的光滑有界的函數(shù),Yt是一個(gè)擴(kuò)散過程。 與Black-Scholes模型中的常數(shù)波動(dòng)率的假設(shè)不同, 式中的波動(dòng)率是隨機(jī)波動(dòng)率。 假設(shè)驅(qū)動(dòng)隨機(jī)波動(dòng)率的擴(kuò)散過程是一個(gè)快速均值回歸的OU過程,Yt滿足以下隨機(jī)微分方程

        (2)

        記為P*風(fēng)險(xiǎn)中性測(cè)度,γt為光滑的有界函數(shù),表示波動(dòng)率風(fēng)險(xiǎn)的市場(chǎng)價(jià)格。 根據(jù)Girsanov定理, 可得到在風(fēng)險(xiǎn)中性測(cè)度下的隨機(jī)波動(dòng)率模型

        (3)

        V(t,s,y)=E*[e-r(T-t)φ(ST)|St=s,Yt=y]

        (4)

        根據(jù)Feynman-Kac定理可知, 期權(quán)在t時(shí)刻的價(jià)格V(t,s,y)滿足以下偏微分方程和終值條件

        (5)

        V(T,s,y)=φ(s)

        (6)

        將歐式期權(quán)在t時(shí)刻的價(jià)格V(t,s,y)以如下形式展開

        (7)

        利用Fouque等[10]提出的奇異攝動(dòng)方法, 可以由上述偏微分方程和價(jià)格展開式得到關(guān)于歐式期權(quán)價(jià)格近似表達(dá)式。 由該方法得到的關(guān)于歐式期權(quán)價(jià)格的三個(gè)主要結(jié)論:

        結(jié)論1展開式(7)中的第一項(xiàng)與y無關(guān),且V0(t,s)滿足以下偏微分方程和終值條件

        (8)

        V0(T,s)=φ(s)

        (9)

        (10)

        其中τ是期權(quán)合約的到期期限,即τ=T-t;C1和C2是待定系數(shù)。

        (11)

        這種方法的優(yōu)點(diǎn)在于可以有效地減少隨機(jī)波動(dòng)率模型中需要估計(jì)的參數(shù)的數(shù)量。 隨機(jī)波動(dòng)率模型涉及到的模型參數(shù)包括波動(dòng)率擴(kuò)散過程的均值水平、波動(dòng)率的方差、波動(dòng)率的均值回歸速率、模型中兩個(gè)隨機(jī)過程的相關(guān)系數(shù)以及波動(dòng)率風(fēng)險(xiǎn)的市場(chǎng)價(jià)格。 利用Fouque等[10]的方法,可以將上述五個(gè)模型參數(shù)包含到C1和C2兩個(gè)系數(shù)中。 而C1和C2這兩個(gè)系數(shù)可以通過金融市場(chǎng)中的看漲期權(quán)的隱含波動(dòng)率的估計(jì)值得到。 這兩個(gè)參數(shù)的估計(jì)方法在第3部分介紹。

        2 模型參數(shù)估計(jì)

        在給定期權(quán)到期時(shí)間T和執(zhí)行價(jià)格K的條件下,可以從金融市場(chǎng)上觀測(cè)到看漲期權(quán)在當(dāng)前時(shí)刻的市場(chǎng)價(jià)格。 記該看漲期權(quán)的市場(chǎng)價(jià)格為Cmarket(K,T)。 所謂隱含波動(dòng)率(Implied Volatility)指的是使金融市場(chǎng)觀測(cè)到的實(shí)際期權(quán)價(jià)格與模型計(jì)算得到的理論價(jià)格相等的波動(dòng)率數(shù)值。 一般將Black-Scholes期權(quán)定價(jià)公式作為計(jì)算隱含波動(dòng)率的理論值的基準(zhǔn)模型。 記隱含波動(dòng)率為IV,根據(jù)隱含波動(dòng)率的定義以下關(guān)系成立

        Cmarket(K,T)=BSC(t,s,IV,K,T)

        (12)

        =sN(d1)-Ke-r(T-t)N(d2)+

        (13)

        其中,

        =sN(d1)-Ke-r(T-t)N(d2)

        (14)

        (15)

        (16)

        BSC(t,s,IV,K,T)

        (17)

        將式(17)和市場(chǎng)實(shí)際價(jià)格的近似表達(dá)式代入式(12)后, 可知隱含波動(dòng)率滿足以下關(guān)系式

        (18)

        期權(quán)的在值程度(moneyness)指的是期權(quán)的實(shí)值或虛值程度, 定義為期權(quán)執(zhí)行價(jià)格與標(biāo)的資產(chǎn)價(jià)格的比值。 期權(quán)的隱含波動(dòng)率變化可以通過期權(quán)的在值程度和期權(quán)的到期期限反映出來。 由(18)式能夠進(jìn)一步推導(dǎo)出隱含波動(dòng)率IV,期權(quán)的在值程度K/s與期權(quán)的到期期限τ之間的近似線性關(guān)系。

        引理1隱含波動(dòng)率IV,期權(quán)的在值程度K/s與期權(quán)的到期期限τ滿足以下線性關(guān)系。

        (19)

        (20)

        (21)

        根據(jù)式(19),可以利用線性回歸方法估計(jì)出斜率α和β截距的值, 再利用(20)和(21)推算出C1和C2這兩個(gè)參數(shù)。

        3 實(shí)證分析及應(yīng)用

        表1 斜率和截距估計(jì)值系數(shù)

        得到斜率和截距的估計(jì)值后,可以根據(jù)(20)和(21)計(jì)算出快速均值回歸隨機(jī)波動(dòng)率模型中的期權(quán)定價(jià)參數(shù)C1和C2:

        (22)

        (23)

        (24)

        從(24)可知,雙限期權(quán)在到期時(shí)刻的價(jià)格等價(jià)于包含K1單位貨幣,一個(gè)多頭頭寸執(zhí)行價(jià)格為K1的看漲期權(quán)和一個(gè)空頭頭寸執(zhí)行價(jià)格為K2的看漲期權(quán)的資產(chǎn)組合。 因此雙限期權(quán)在t時(shí)刻的價(jià)格為

        (25)

        其中C(Ki)為執(zhí)行價(jià)格為Ki的歐式看漲期權(quán)在t時(shí)刻的價(jià)格。 眾所周知,C(Ki)的表達(dá)式為

        (26)

        (27)

        由于(25)是兩個(gè)歐式看漲期權(quán)的線性組合,易知(25)應(yīng)滿足偏微分方程(8)。因此根據(jù)第2部分結(jié)論1可知,式(25)式快速均值回歸隨機(jī)波動(dòng)率框架下,歐式雙限期權(quán)近似解析解的第一項(xiàng)。在得到近似解析解的第一項(xiàng)之后,可以根據(jù)第2部分結(jié)論2得到近似解析解的第二項(xiàng)表達(dá)式

        (28)

        在快速均值回歸隨機(jī)波動(dòng)率假設(shè)下,雙限期權(quán)的價(jià)格Pt可以通過(25)和(28)式近似表示,即

        (29)

        下面通過控制變量的方法進(jìn)一步分析變量和參數(shù)對(duì)模型的影響。 取當(dāng)期標(biāo)的資產(chǎn)的價(jià)格的變動(dòng)范圍為[200,300],其它參數(shù)保持不變。根據(jù)式(25)和(29)可以分別得到常數(shù)波動(dòng)率雙限期權(quán)價(jià)格、隨機(jī)波動(dòng)率雙限期權(quán)價(jià)格與標(biāo)的資產(chǎn)價(jià)格之間的關(guān)系,如圖1所示。根據(jù)式(28)可以得到雙限期權(quán)近似表達(dá)式的修正項(xiàng)與標(biāo)的資產(chǎn)之間的關(guān)系,如圖2所示。

        由圖1可以看出,在考慮了隨機(jī)波動(dòng)率之后,雙限期權(quán)的價(jià)格比常數(shù)波動(dòng)率模型的價(jià)格較低。尤其在標(biāo)的資產(chǎn)價(jià)格超過雙期權(quán)合約兩個(gè)執(zhí)行價(jià)格的均值之后,兩種不同模型給出的價(jià)格差異會(huì)更加明顯。這個(gè)結(jié)果說明了,若只通過歷史波動(dòng)率和Black-Scholes常數(shù)波動(dòng)率模型定價(jià),很可能會(huì)高估了雙限期權(quán)的價(jià)格,從而造成交易損失。由圖2可知,在標(biāo)的資產(chǎn)變化范圍內(nèi),雙限期權(quán)近似價(jià)格的修正項(xiàng)取值為負(fù)數(shù),并且有隨著標(biāo)的資產(chǎn)價(jià)格增加而減少的趨勢(shì)。因此考慮了快速回歸隨機(jī)波動(dòng)率的雙限期權(quán)價(jià)格小于常數(shù)波動(dòng)率模型的價(jià)格主要是因?yàn)槎▋r(jià)公式中的修正項(xiàng)的影響,反映了隨機(jī)波動(dòng)率對(duì)期權(quán)價(jià)格的影響程度。

        圖1 期權(quán)價(jià)格比較

        圖2 雙限期權(quán)修正項(xiàng)

        圖3 執(zhí)行價(jià)格對(duì)期權(quán)價(jià)格的影響

        取當(dāng)期標(biāo)的資產(chǎn)的價(jià)格的變動(dòng)范圍為[100,300],期權(quán)執(zhí)行價(jià)格K2的變動(dòng)范圍為[200,300],其它參數(shù)保持不變。根據(jù)式(29)可以得到具有不同執(zhí)行價(jià)格的隨機(jī)波動(dòng)率雙限期權(quán)價(jià)格與標(biāo)的資產(chǎn)價(jià)格之間的關(guān)系,如圖3所示。從圖3可知,隨著執(zhí)行價(jià)格K2的增大,隨機(jī)波動(dòng)率雙限期權(quán)的價(jià)格也較高。這個(gè)結(jié)論與常數(shù)波動(dòng)率模型的定價(jià)是一致的,因此在隨機(jī)波動(dòng)率模型中,執(zhí)行價(jià)格的變化并不會(huì)對(duì)雙限期權(quán)的定價(jià)產(chǎn)生影響。

        在上述的實(shí)證分析和應(yīng)用中,我們首先利用金融市場(chǎng)上流動(dòng)性較大的看漲期權(quán)的隱含波動(dòng)率曲線估計(jì)出均值回歸隨機(jī)波動(dòng)率模型中的兩個(gè)參數(shù),進(jìn)而應(yīng)用到雙限期權(quán)的定價(jià)中。該方法與其它模型相比,具有兩個(gè)較為明顯的特點(diǎn)。首先減少了隨機(jī)波動(dòng)率模型(3)中需要估計(jì)的參數(shù)的個(gè)數(shù),提高了效率;其次,可以先從流動(dòng)性較大的期權(quán)(如歐式看漲期權(quán))中估計(jì)出求解偏微分方程的參數(shù),再應(yīng)用到其它流動(dòng)性較小的期權(quán)(雙限期權(quán))定價(jià)中,解決了某些奇異期權(quán)由于交易數(shù)據(jù)稀缺無法進(jìn)行有效參數(shù)估計(jì)的問題。

        4 結(jié)論

        本文基于快速均值回歸隨機(jī)波動(dòng)率模型,研究了雙限期權(quán)的定價(jià)問題。雙限期權(quán)可以對(duì)沖投資者由于標(biāo)的資產(chǎn)價(jià)格波動(dòng)造成的投資風(fēng)險(xiǎn),因此研究雙限期權(quán)定價(jià)問題對(duì)于金融市場(chǎng)中的現(xiàn)貨和期權(quán)交易,以及金融市場(chǎng)的風(fēng)險(xiǎn)管理和控制具有一定的理論和現(xiàn)實(shí)意義。由于傳統(tǒng)的期權(quán)定價(jià)模型中關(guān)于常數(shù)波動(dòng)率的假設(shè)無法解釋隱含波動(dòng)率的“微笑曲線”現(xiàn)象,而快速均值回歸隨機(jī)波動(dòng)率模型不但能夠較好地解釋隱含波動(dòng)率的“微笑曲線”現(xiàn)象,還能夠較好地解釋標(biāo)的資產(chǎn)收益的“尖峰厚尾”分布和波動(dòng)率的“聚集”效應(yīng),因此在快速均值回歸隨機(jī)波動(dòng)率模型框架下研究雙限期權(quán)定價(jià)比較切合實(shí)際。本文通過市場(chǎng)隱含波動(dòng)率的數(shù)據(jù),估計(jì)了快速均值回歸隨機(jī)波動(dòng)率模型中的兩個(gè)參數(shù),并應(yīng)用到雙限期權(quán)的定價(jià)計(jì)算中。在得到考慮了隨機(jī)波動(dòng)率的雙限期權(quán)定價(jià)公式之后,對(duì)期權(quán)價(jià)格與標(biāo)的資產(chǎn)價(jià)格之間的關(guān)系做了數(shù)值分析?;诳焖倬祷貧w隨機(jī)波動(dòng)率模型的期權(quán)定價(jià)方法可以進(jìn)一步應(yīng)用到其它的衍生品定價(jià)和資產(chǎn)組合管理中,這也是下一步研究的重要方向。

        猜你喜歡
        標(biāo)的模型
        一半模型
        重要模型『一線三等角』
        重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
        具有兩個(gè)P’維非線性不可約特征標(biāo)的非可解群
        紅周刊績(jī)優(yōu)指數(shù)100只標(biāo)的股一覽
        紅周刊績(jī)優(yōu)指數(shù)100只標(biāo)的股一覽
        紅周刊績(jī)優(yōu)指數(shù)100只標(biāo)的股一覽
        紅周刊績(jī)優(yōu)指數(shù)100只標(biāo)的股一覽
        紅周刊績(jī)優(yōu)指數(shù)100只標(biāo)的股一覽
        3D打印中的模型分割與打包
        国产成人精品三级麻豆| 天天干天天日夜夜操| 久久久无码精品亚洲日韩蜜臀浪潮| 精品久久久久久777米琪桃花| 国产一级黄色录像| 亚洲一区二区三区av无| 蜜芽亚洲av无码精品色午夜| 中文字幕v亚洲日本| 亚洲无亚洲人成网站77777| 国产人妖视频一区二区| 久久婷婷国产综合精品| 男女动态91白浆视频| 国产三级精品三级| 天天影视色香欲综合久久| 欧美亚洲国产丝袜在线| 久久精品国产亚洲av久五月天| 台湾佬中文娱乐网22| 亚洲色精品aⅴ一区区三区| 国产无线乱码一区二三区| 久久亚洲国产成人精品v| 日本国主产一区二区三区在线观看| 91成人自拍在线观看| 九九热线有精品视频86| 久久亚洲国产中v天仙www| 成人性生交大片免费看7| 色大全全免费网站久久| 黄色a级国产免费大片| 亚洲黄视频| 少妇被搞高潮在线免费观看| 人妻少妇精品视频专区vr| 成年女人vr免费视频| 精品无码国产污污污免费| 久久免费网站91色网站| 亚洲精品一区二区三区新线路| 国产精品亚洲片在线观看不卡| 国产裸体歌舞一区二区| 一本久久精品久久综合桃色| 日韩精品极品免费视频观看| 玩弄放荡人妻少妇系列视频| 久久狠色噜噜狠狠狠狠97| 久久久国产熟女综合一区二区三区|