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

    基于模糊時(shí)間序列模型的股票市場(chǎng)預(yù)測(cè)

    2010-10-21 06:25:26藺玉佩楊一文
    統(tǒng)計(jì)與決策 2010年8期
    關(guān)鍵詞:模糊化論域聚類

    藺玉佩,楊一文

    (西北工業(yè)大學(xué) 管理學(xué)院,西安 710129)

    0 引言

    股票價(jià)格序列預(yù)測(cè)不僅對(duì)投資者而言有巨大的市場(chǎng)利益,同時(shí)對(duì)學(xué)術(shù)界研究市場(chǎng)運(yùn)行規(guī)律無(wú)疑也具有重要的理論價(jià)值。因此研究采用何種方法對(duì)股票價(jià)格進(jìn)行較準(zhǔn)確的預(yù)測(cè)就成為是國(guó)內(nèi)外學(xué)者都非常關(guān)注的課題。時(shí)間序列預(yù)測(cè)方法是常用的股票價(jià)格預(yù)測(cè)方法。一般的時(shí)間序列預(yù)測(cè)方法,如ARIMA等,都是以時(shí)間序列的精確值為基礎(chǔ)并且認(rèn)為時(shí)間序列的未來(lái)值與當(dāng)前值、過(guò)去值以及白噪聲之間存在著確定、明確的函數(shù)關(guān)系。但是,證券市場(chǎng)系統(tǒng)非常復(fù)雜,市場(chǎng)的行為特征是不明確的,精確的記錄往往會(huì)丟失部分有用的信息,許多變量之間的關(guān)系也難以用確定的函數(shù)關(guān)系來(lái)描述。因此,本文將模糊時(shí)間序列模型用于股票指數(shù)價(jià)格的預(yù)測(cè),將觀測(cè)值表示為具有模糊特征的語(yǔ)言變量并且建立了多前件的模糊關(guān)系,最后通過(guò)去模糊化得到精確值。結(jié)果表明此法更能全面反映系統(tǒng)的特征并且有利于復(fù)雜環(huán)境中提高預(yù)測(cè)精度。

    1 模糊時(shí)間序列與模糊邏輯關(guān)系

    本文應(yīng)用模糊時(shí)間序列模型對(duì)上證綜指和深證成指進(jìn)行預(yù)測(cè),首先對(duì)模糊時(shí)間序列的定義作一簡(jiǎn)要敘述。

    令U為給定論域,將論域劃分為n個(gè)子區(qū)間,則U={u1,u2,…,un}。一個(gè)定義在論域U中的模糊集合A表示如下:

    其中,fA(·)是 ui對(duì)模糊集合 A 的隸屬函數(shù),fA(ui)是 ui對(duì)模糊集合 A 的隸屬度,fA(ui)∈[0,1],i=1,2,…,n。

    定義1 令R中一子集Y(t)(t=…,0,1,2,…)為給定論域,fi(t)(i=1,2,…)為定義在其上的模糊集合,由 f1(t),f2(t)(組成的集合F(t)稱作定義在Y(t)(t=…,0,1,2,…)上的模糊時(shí)間序列。

    定義 2 如果 F(t)僅僅由 F(t-1)引起,即 F(t-1)→F(t),這個(gè)關(guān)系可表示為F(t)=F(t-1)°R(t,t-1),稱F(t)作的一階模型,R(t,t-1)為F(t-1)與F(t)之間的模糊關(guān)系。

    定義2中,F(xiàn)(t-1)為模糊關(guān)系的前件,F(xiàn)(t)為后件,此關(guān)系式是具有單個(gè)前件的一階模糊關(guān)系式。

    定義3 F(t)為給定論域上的模糊時(shí)間序列,如果F(t)同時(shí)由F(t-1),F(t-2),…,F(t-n)引起,則n階模糊邏輯關(guān)系可表示為 F(t-n)、…、F(t-2)、F(t-1)→F(t),其中 F(t-n)、…、F(t-2)、F(t-1)為n階模糊邏輯關(guān)系的當(dāng)前狀態(tài)。

    定義4 假設(shè)R中一子集Y(t)(t=…,0,1,2,…)為給定論域,F(xiàn)1(t),F2(t),…,Fn(t)是各自給定論域上的n個(gè)模糊時(shí)間序列,如果 F(t)同時(shí)由 F1(t),F2(t-1),F3(t-1),…,Fn(t-1)引起,具有 n個(gè)前件的一階模糊關(guān)系可表示為Fn(t-1),…,F2(t-1),F1(t-1)→F(t)。

    綜合以上定義,具有n個(gè)前件的n階模糊邏輯關(guān)系式可定義如下:

    定義5 F1(t),F2(t),…,Fn(t)和F(t)為各自給定論域上的n+1個(gè)模糊時(shí)間序列,則具有n個(gè)前件的n階模糊邏輯關(guān)系式可表示為:(F1(t-n),F2(t-n),…,Fn(t-n))、…、(F1(t-2),F2(t-2),…,Fn(t-2))、(F1(t-1),F2(t-1),…,Fn(t-1))→F(t)。

    2 模糊時(shí)間序列預(yù)測(cè)模型

    本文建立的模糊時(shí)間序列模型主要包括以下步驟:模糊化、建立模糊關(guān)系、預(yù)測(cè)并去模糊化。

    步驟1 模糊化

    定義論域 U。 論域 U=[Dmin-σ,Dmax+σ],其中 Dmin和 Dmax分別為訓(xùn)練樣本數(shù)據(jù)中的最小值和最大值,σ為訓(xùn)練樣本標(biāo)準(zhǔn)差。采用模糊C-均值聚類算法將數(shù)據(jù)聚類,并根據(jù)聚類結(jié)果論域劃分為多個(gè)子區(qū)間,首先計(jì)算聚類數(shù)k,則

    其中,n為訓(xùn)練樣本個(gè)數(shù),Dmin和Dmax分別為訓(xùn)練樣本序列中的最大值和最小值,X(t)和X(t-1)分別為t時(shí)刻和t-1時(shí)刻的訓(xùn)練樣本的觀測(cè)值。[]為四舍五入取整運(yùn)算。其次通過(guò)模糊C-均值聚類算法將數(shù)據(jù)分為k類,相應(yīng)地得到k個(gè)聚類中心ji(i=1,2,…,k)。將相鄰兩個(gè)聚類中心的中點(diǎn)做為論域子區(qū)間的邊界點(diǎn),定義邊界點(diǎn)為di(i=1,2,…,k-1)。劃出k個(gè)子區(qū)間:(Dmin,d1)、(d1,d2)、(d2,d3)、…、(dk-1,Dmax),分別用 u1、u2、…、uk表示。可以看出,各子區(qū)間長(zhǎng)度各異,這是因?yàn)闃颖緮?shù)據(jù)在其論域上的分布不是均勻的,而是有其復(fù)雜的內(nèi)部結(jié)構(gòu),用聚類結(jié)果來(lái)指導(dǎo)子區(qū)間的劃分顯然比等分論域能更好地反映數(shù)據(jù)結(jié)構(gòu),有利于提高預(yù)測(cè)精度。

    最后定義子區(qū)間對(duì)應(yīng)的模糊集合,并將數(shù)據(jù)模糊化。定義模糊集合Ai為:

    其中,fij表示uj對(duì)模糊集合Ai的隸屬度(i=1,2,…,k;j=1,2,…,k)。數(shù)據(jù)的模糊化規(guī)則為:如果某樣本數(shù)據(jù)屬于ui,且ui對(duì)Ai的隸屬度在其所有隸屬度中是最大值,則可將該數(shù)據(jù)值模糊化為Ai。將所有樣本數(shù)據(jù)模糊化,得到對(duì)應(yīng)的模糊時(shí)間序列。將股票每日收益率、5日移動(dòng)均線變化值、成交量和DHL四個(gè)時(shí)間序列分別定義為{R(t)}、{M(t)}、{V(t)}和{D(t)}(t=1,2,…,n),其中n為訓(xùn)練樣本數(shù)。根據(jù)式(1)得出聚類數(shù)分別為k1,k2,k3和k4,在各自的論域上將上述4個(gè)時(shí)間序列模糊化后得到的模糊時(shí)間序列分別為 A(t,i)(i=1,2,…,k1),B(t,j)(j=1,2,…,k2),C(t,p)(p=1,2,…,k3)和D(t,q)(q=1,2,…,k4;t=1,2,…,n)。

    步驟2 建立模糊關(guān)系

    選用n天的樣本數(shù)據(jù)作為建模的訓(xùn)練數(shù)據(jù),設(shè)當(dāng)天為t,前三天分別為t-3,t-2和t-1,t=1,2,…,n。利用上述模糊時(shí)間序列建立具有四前件、三階模糊關(guān)系的模糊時(shí)間序列預(yù)測(cè)模型:

    其中 i1,i2,i3,i4=1,2,…,k1;j1,j2,j3=1,2,…,k2;p1,p2,p3=1,2,…,k3;q1,q2,q3=1,2,…,k4。則按時(shí)間順序可以依次建立共計(jì)n-3個(gè)模糊關(guān)系,形成訓(xùn)練數(shù)據(jù)模糊關(guān)系表,如表1所示。

    步驟3 預(yù)測(cè)并去模糊化

    (1)預(yù)測(cè)

    預(yù)測(cè)未來(lái)T時(shí)刻的模糊收益率A(T,i4),其模糊邏輯關(guān)系式為:

    表1 訓(xùn)練數(shù)據(jù)模糊關(guān)系

    將上式左邊的前件依次與表1中列出的n-3個(gè)關(guān)系式的前件逐個(gè)進(jìn)行對(duì)比。具體地,本文提出的對(duì)比原則是:計(jì)算前件相同位置上第二下標(biāo)的差值的絕對(duì)值,例如,模糊集合A(T-3,i3)中 i3=m3(1≤m3≤k1),A(t-3,i3)中 i3=n3(1≤n3≤k1),則差值絕對(duì)值為|m3-n3|,將12對(duì)模糊集合第二下標(biāo)的差值絕對(duì)值求和,如果和小于等于某一設(shè)定值(如20),則認(rèn)為關(guān)系匹配成功。這樣做的優(yōu)點(diǎn)是整體度量?jī)蓚€(gè)模糊關(guān)系前件的差異,設(shè)定的值越小,差異度量越精確,但匹配成功的次數(shù)將隨之降低,反之亦然。如果一次預(yù)測(cè)關(guān)系匹配成功的個(gè)數(shù)為N,則這N個(gè)訓(xùn)練數(shù)據(jù)模糊關(guān)系式的后件構(gòu)成了本次預(yù)測(cè)結(jié)果集合A(·,·),該集合中模糊集合出現(xiàn)的次數(shù)不一定相同,有些可能多次出現(xiàn),有些則可能僅出現(xiàn)一次。設(shè)各可能出現(xiàn)結(jié)果值的頻數(shù)為 fi(i=1,2,…,k1),預(yù)測(cè)結(jié)果集合如表 2所示,集合中的每一個(gè)模糊集合對(duì)應(yīng)一個(gè)聚類中心。

    (2)去模糊化

    采用重心法去模糊化。根據(jù)表2結(jié)果去模糊化得到預(yù)測(cè)值

    3 應(yīng)用實(shí)例

    3.1 數(shù)據(jù)說(shuō)明

    選取上證綜合指數(shù)和深證成份指數(shù)每日收益率、收盤價(jià)、最高價(jià)、最低價(jià)和成交量為樣本,區(qū)間分別為2005-12-14至2009-4-15和2005-12-26至2009-4-15①數(shù)據(jù)來(lái)源于清華金融研究數(shù)據(jù)庫(kù)(http://211.157.28.243/terminal/system/gotoLogin.action)。。收益率采用對(duì)數(shù)收益率,DHL由最高價(jià)、最低價(jià)和收盤價(jià)計(jì)算得到。由于成交量數(shù)值比較大,將其單位改為萬(wàn)手。最后得到上證綜指有效數(shù)據(jù)共800天,深證成指共794天。

    3.2 預(yù)測(cè)步驟

    依照前文的建模步驟,使用上證指數(shù)的數(shù)據(jù)來(lái)說(shuō)明預(yù)測(cè)過(guò)程。

    步驟1 模糊化

    取前750個(gè)交易日數(shù)據(jù)作為訓(xùn)練樣本,后50個(gè)交易日數(shù)據(jù)用于預(yù)測(cè)。以上證指數(shù)收益率的模糊化過(guò)程為例。

    表2 一次模糊預(yù)測(cè)所有可能的結(jié)果

    首先定義論域U。收益率訓(xùn)練樣本中最大值和最小值分別為 0.0903 和-0.0926,σ=0.0258。 得到 U=[-0.1184,0.1161]。由于市場(chǎng)實(shí)施10%的漲跌幅限制,在[-0.1,0.1]范圍以外沒(méi)有值,所以我們定義收益率論域U=[-0.1,0.1]。

    其次劃分子區(qū)間。根據(jù)式(1)可得分類數(shù)k1=10,使用matlab模糊邏輯工具箱函數(shù)中的模糊C-均值聚類函數(shù)fcm將訓(xùn)練數(shù)據(jù)進(jìn)行分類,得到10個(gè)聚類中心分別為:-0.0789,-0.0483,-0.03,-0.0181,-0.0067,0.0015,0.0089,0.0167,0.0267,0.0446,以相鄰兩個(gè)聚類中心中點(diǎn)作為子區(qū)間的分界點(diǎn),將論域劃分成10個(gè)子區(qū)間分別為u1=[-0.1,-0.0636),u2=[-0.0636,-0.0392),u3=[-0.0392,-0.024),u4=[-0.024,-0.0124),u5=[-0.024,-0.0026),u6=[-0.0026,0.0052),u7=[0.0052,0.0128),u8=[0.0128,0.0217),u9=[0.0217,0.0356),u10=[0.0356,0.1]。

    再根據(jù)式(2)定義模糊集合為:

    最后將數(shù)據(jù)模糊化。例如,2006-11-17的收益率為0.0155,它落在區(qū)間u8內(nèi),u8對(duì)于A8的隸屬度在所有區(qū)間中最大,所以將0.0155模糊化為A8。將收益序列中所有數(shù)據(jù)模糊化后可得對(duì)數(shù)收益率的模糊時(shí)間序列。同樣方法將其它三個(gè)序列的數(shù)據(jù)在各自論域上模糊化為模糊時(shí)間序列。

    步驟2 建立模糊關(guān)系

    利用750日的訓(xùn)練數(shù)據(jù)可構(gòu)造747(=750-3)個(gè)模糊關(guān)系如表3所示,然后將預(yù)測(cè)T日的模糊關(guān)系式(5)中的前件與這747個(gè)模糊關(guān)系的前件進(jìn)行比對(duì),選出匹配成功的模糊關(guān)系。需要指出的是,為了書寫簡(jiǎn)潔,又不致產(chǎn)生歧義,后文中省略了模糊時(shí)間序列中模糊集合的第一個(gè)(時(shí)序)下標(biāo)。

    表3 上證指數(shù)訓(xùn)練數(shù)據(jù)模糊關(guān)系

    步驟3 預(yù)測(cè)并去模糊化

    (1)預(yù)測(cè)

    以2009-2-4這一天為例說(shuō)明預(yù)測(cè)過(guò)程。當(dāng)日收益率為Ar。用前3天歷史數(shù)據(jù)并根據(jù)式(5)可建立模糊邏輯關(guān)系式為:

    將其前件中的模糊集合與表3中各關(guān)系式前件中對(duì)應(yīng)的模糊集合根據(jù)對(duì)比原則進(jìn)行一一對(duì)比,最后匹配成功的關(guān)系式個(gè)數(shù)為58個(gè),這58個(gè)模糊關(guān)系式的后件構(gòu)成了預(yù)測(cè)結(jié)果集合如表4所示。

    表4 上證指數(shù)2009-2-4模糊預(yù)測(cè)結(jié)果

    (2)去模糊化

    根據(jù)公式(4)可以得出2009-2-4日的收益率預(yù)測(cè)結(jié)果為0.00142239。當(dāng)天的實(shí)際收益率為0.023。

    4 預(yù)測(cè)結(jié)果分析

    對(duì)上證指數(shù)和深證指數(shù)從2009-2-4至2009-4-15的50天價(jià)格的預(yù)測(cè)結(jié)果如圖1、圖2和表5所示。為客觀評(píng)價(jià)本文提出的模糊時(shí)間序列模型的預(yù)測(cè)效果,我們特意將其與模糊時(shí)間序列中的標(biāo)志模型進(jìn)行對(duì)比,結(jié)果如表5所示。

    由表5(第2、4列)可以看出,對(duì)于同樣兩市指數(shù)的50日預(yù)測(cè),Chen的模型[4]的平均預(yù)測(cè)誤差接近或高于2.5% ,而本文建立的模型的平均預(yù)測(cè)誤差都明顯低于2%,預(yù)測(cè)結(jié)果的改進(jìn)還是很明顯的。同時(shí)也說(shuō)明本文在論域劃分以及模型變量的選擇方面所作的努力還是有一定成效的,同時(shí)也表明,對(duì)于處理復(fù)雜問(wèn)題的模糊時(shí)間序列而言,在論域劃分和模型變量選擇應(yīng)予以高度重視。

    價(jià)格指數(shù)數(shù)值的預(yù)測(cè)固然重要,但是退一步,其未來(lái)的變化趨勢(shì)即漲跌無(wú)疑也是非常重要的信息。應(yīng)用本文建立的模糊時(shí)間序列模型分別對(duì)上證指數(shù)、深證指數(shù)未來(lái)漲跌趨勢(shì)進(jìn)行預(yù)測(cè),預(yù)測(cè)區(qū)間:2009-2-4至2009-2-24。預(yù)測(cè)的過(guò)程同上,只是省了去模糊化這一步驟,直接在論域上找出預(yù)測(cè)得到的模糊集合(后件)所對(duì)應(yīng)的子區(qū)間ui(i=1,2,…,k1),進(jìn)而判斷指數(shù)變動(dòng)趨勢(shì)。同樣地,我們也將其與Chen的模型進(jìn)行了對(duì)比,結(jié)果如表5(第3、5列)所示。Chen的模型對(duì)兩市指數(shù)的趨勢(shì)預(yù)測(cè)正確率均低于60%,而本文模型的預(yù)測(cè)正確率分別接近或高于70%。

    5 結(jié)論

    本文將時(shí)間序列分析與模糊規(guī)則結(jié)合起來(lái)所產(chǎn)生的模糊時(shí)間序列模型具有兩者的優(yōu)點(diǎn),是時(shí)間序列分析方法的拓展,將其應(yīng)用于股票價(jià)格的預(yù)測(cè)是一種有益的嘗試。

    模糊時(shí)間序列模型的預(yù)測(cè)精度主要取決于兩點(diǎn):論域的子區(qū)間劃分和模糊關(guān)系的建立。前者要求合理劃分論域,反映真實(shí)的數(shù)據(jù)分布結(jié)構(gòu),并以此為基礎(chǔ)選擇適當(dāng)?shù)娜ツ:椒?;后者則要求不僅對(duì)數(shù)據(jù)的擬合精度高,而且還要求模型具有較好的推廣能力,為達(dá)到這一目的,必須選擇適當(dāng)?shù)哪P妥兞浚辞凹兞?,同時(shí)還要建立有效的模糊規(guī)則。因此本文選擇模糊聚類方法實(shí)現(xiàn)對(duì)論域的子區(qū)間劃分,力求客觀地反映樣本數(shù)據(jù)的結(jié)構(gòu)分布,并為劃分論域提供一個(gè)比較客觀的方法。相應(yīng)地,采用重心法區(qū)模糊化,實(shí)現(xiàn)對(duì)收益的精確預(yù)測(cè)。影響價(jià)格的因素很多,但模型中又不可能引入過(guò)多的變量,否則不僅導(dǎo)致運(yùn)算量過(guò)大,而且使模糊關(guān)系過(guò)于復(fù)雜,使模型的推廣能力下降?;谏鲜隹紤],本文選擇歷史收益率、5日移動(dòng)平均線、成交量以及多空雙方強(qiáng)弱對(duì)比4個(gè)變量作為模型的變量,建立了多前件的模糊時(shí)間序列模型,并特別設(shè)計(jì)了能反映多空雙方力量對(duì)比的指標(biāo)DHL以及模糊規(guī)則的確定原則。預(yù)測(cè)結(jié)果表明,本文建立的模糊時(shí)間序列模型較典型的模糊時(shí)間序列模型具有較高的預(yù)測(cè)精度。

    表5 上證、深證指數(shù)預(yù)測(cè)誤差

    由于所面臨問(wèn)題的復(fù)雜性,模糊時(shí)間序列模型的運(yùn)用是問(wèn)題導(dǎo)向的。即論域的劃分、模糊關(guān)系的建立必須根據(jù)具體問(wèn)題設(shè)計(jì)合理的方法,這些方面仍有許多工作要做。例如模糊關(guān)系前件比對(duì)過(guò)程中,本文選擇的20這一數(shù)值僅僅是經(jīng)過(guò)累試后得到的,工作量比較大;另外,論域的劃分與模糊規(guī)則的建立目前還是獨(dú)立的兩部分內(nèi)容,如果將二者聯(lián)合起來(lái)綜合考慮,對(duì)提高模型精度無(wú)疑具有重要意義。

    [1]Song,Q.,Chissom,B.S.Fuzzy Time Series and Its Models[J].Fuzzy Sets and Systems,1993,(54).

    [2]Song,Q.,Chissom,B.S.Forecasting Enrollments with Fuzzy Time series-Part I[J].Fuzzy Sets and Systems,1993,(54).

    [3]吳銘鋒,蔣勛.基于模糊時(shí)間序列的預(yù)測(cè)模型-以上證指數(shù)為例[J].價(jià)值工程,2008,(11).

    [4]何云峰,楊燕.基于模糊時(shí)間序列-股票走勢(shì)的建模與應(yīng)用[J].微計(jì)算機(jī)信息,2006,(33).

    [5]張智星,孫春在.神經(jīng)-模糊和軟計(jì)算[M].西安:西安交通大學(xué)出版社,2000.

    猜你喜歡
    模糊化論域聚類
    ([0,1],[0,1])-模糊擬陣的基和秩函數(shù)
    基于變論域模糊控制的Taylor逼近型內(nèi)模PID算法
    餐飲娛樂(lè)空間的“邊界模糊化”態(tài)勢(shì)探討——餐飲娛樂(lè)空間設(shè)計(jì)專輯
    三角模糊數(shù)去模糊化對(duì)VIKOR妥協(xié)解的影響研究
    變論域自適應(yīng)模糊PID控制系統(tǒng)仿真與應(yīng)用
    基于DBSACN聚類算法的XML文檔聚類
    雙論域粗糙集在故障診斷中的應(yīng)用
    微生物燃料電池的變論域自適應(yīng)模糊控制研究
    基于改進(jìn)的遺傳算法的模糊聚類算法
    一種層次初始的聚類個(gè)數(shù)自適應(yīng)的聚類方法研究
    久久精品国产鲁丝片午夜精品| 精品免费久久久久久久清纯| 七月丁香在线播放| 亚洲欧美精品专区久久| 亚洲第一区二区三区不卡| 亚洲色图av天堂| 国产探花极品一区二区| 欧美变态另类bdsm刘玥| 中文字幕人妻熟人妻熟丝袜美| 三级毛片av免费| av女优亚洲男人天堂| 国产亚洲av片在线观看秒播厂 | 日本黄色视频三级网站网址| 熟妇人妻久久中文字幕3abv| 天堂av国产一区二区熟女人妻| 乱码一卡2卡4卡精品| 蜜桃亚洲精品一区二区三区| 国产亚洲av嫩草精品影院| 毛片一级片免费看久久久久| 国产一区二区亚洲精品在线观看| 亚洲精品乱久久久久久| 国产精品,欧美在线| 欧美日韩精品成人综合77777| 免费看a级黄色片| 91久久精品电影网| 两个人的视频大全免费| 夫妻性生交免费视频一级片| 看十八女毛片水多多多| 国产伦理片在线播放av一区| 日韩在线高清观看一区二区三区| 久久久精品大字幕| 亚洲av成人精品一区久久| 亚洲不卡免费看| 成年av动漫网址| 91精品国产九色| 啦啦啦啦在线视频资源| 非洲黑人性xxxx精品又粗又长| 国产麻豆成人av免费视频| 亚洲欧美日韩高清专用| 最近中文字幕高清免费大全6| 69人妻影院| 成人漫画全彩无遮挡| 国产精品精品国产色婷婷| 中文字幕熟女人妻在线| 97在线视频观看| 九色成人免费人妻av| 五月玫瑰六月丁香| 国产黄a三级三级三级人| 五月玫瑰六月丁香| 天堂√8在线中文| 五月玫瑰六月丁香| 国产成年人精品一区二区| 日日撸夜夜添| 日韩制服骚丝袜av| 久久久久久伊人网av| 免费av毛片视频| 美女国产视频在线观看| 亚洲美女视频黄频| 亚洲美女视频黄频| 国产私拍福利视频在线观看| 国产精品一及| 长腿黑丝高跟| 一级毛片电影观看 | 能在线免费看毛片的网站| 毛片女人毛片| 中文乱码字字幕精品一区二区三区 | 国产白丝娇喘喷水9色精品| 日韩国内少妇激情av| 亚洲人成网站在线观看播放| 国产久久久一区二区三区| 亚洲人成网站在线播| 亚洲av不卡在线观看| 免费看日本二区| 女人被狂操c到高潮| 成年免费大片在线观看| 老司机福利观看| 色尼玛亚洲综合影院| 国产精品久久电影中文字幕| 亚洲久久久久久中文字幕| 国产91av在线免费观看| 亚洲美女搞黄在线观看| 欧美成人午夜免费资源| 久久精品国产鲁丝片午夜精品| 国产综合懂色| 男女那种视频在线观看| 夜夜爽夜夜爽视频| 天堂av国产一区二区熟女人妻| 97超碰精品成人国产| 国产精品99久久久久久久久| 亚洲欧美日韩东京热| 黄色一级大片看看| av福利片在线观看| 久久精品91蜜桃| 国产精品一区二区性色av| 中文天堂在线官网| 汤姆久久久久久久影院中文字幕 | 精品国内亚洲2022精品成人| 免费一级毛片在线播放高清视频| 午夜免费男女啪啪视频观看| 日日摸夜夜添夜夜添av毛片| 熟妇人妻久久中文字幕3abv| 国产精品国产高清国产av| 亚洲av二区三区四区| 九九在线视频观看精品| 欧美97在线视频| 男人狂女人下面高潮的视频| 国产91av在线免费观看| 丝袜喷水一区| 成人欧美大片| 亚洲成av人片在线播放无| 不卡视频在线观看欧美| 波野结衣二区三区在线| 亚洲电影在线观看av| 精品人妻偷拍中文字幕| 国产精品国产三级专区第一集| 亚洲自拍偷在线| 久久99热这里只有精品18| 亚洲精品色激情综合| 国产精品日韩av在线免费观看| 嫩草影院精品99| 91久久精品国产一区二区成人| 男人的好看免费观看在线视频| 熟女人妻精品中文字幕| 国产69精品久久久久777片| 国产精品久久电影中文字幕| 老女人水多毛片| 小蜜桃在线观看免费完整版高清| 熟妇人妻久久中文字幕3abv| 插逼视频在线观看| 纵有疾风起免费观看全集完整版 | 国产 一区 欧美 日韩| 联通29元200g的流量卡| 特级一级黄色大片| av福利片在线观看| 中文天堂在线官网| 欧美日韩在线观看h| 亚洲av成人精品一二三区| 美女国产视频在线观看| av在线天堂中文字幕| 色播亚洲综合网| 2021少妇久久久久久久久久久| 国产精品久久久久久精品电影小说 | 国产成人免费观看mmmm| 国产av不卡久久| 我要搜黄色片| 丰满少妇做爰视频| 亚洲国产色片| 免费看av在线观看网站| 日本三级黄在线观看| 天堂网av新在线| 91精品一卡2卡3卡4卡| 国产精品国产高清国产av| 日本黄色视频三级网站网址| 亚洲精品自拍成人| 少妇被粗大猛烈的视频| 亚洲精品乱久久久久久| 99九九线精品视频在线观看视频| 午夜福利在线观看吧| 中文字幕精品亚洲无线码一区| 成人鲁丝片一二三区免费| 亚洲精品乱码久久久v下载方式| 国产成人精品一,二区| 国产又色又爽无遮挡免| 大香蕉久久网| 国产爱豆传媒在线观看| 日本午夜av视频| www.色视频.com| 亚洲国产欧美人成| 欧美zozozo另类| 在线免费十八禁| 18禁动态无遮挡网站| 国产精品电影一区二区三区| 国产成人一区二区在线| 国产极品精品免费视频能看的| 亚洲乱码一区二区免费版| 久久这里只有精品中国| 男女边吃奶边做爰视频| 国产成人a∨麻豆精品| 国产精品美女特级片免费视频播放器| 免费黄色在线免费观看| 国产91av在线免费观看| 可以在线观看毛片的网站| 日韩 亚洲 欧美在线| 欧美xxxx性猛交bbbb| 国产成人a∨麻豆精品| 69av精品久久久久久| 精品欧美国产一区二区三| 国产精品伦人一区二区| 尤物成人国产欧美一区二区三区| 黄色欧美视频在线观看| 久久鲁丝午夜福利片| 国产成人精品久久久久久| 2021天堂中文幕一二区在线观| 精品一区二区免费观看| 国产黄a三级三级三级人| 国产免费视频播放在线视频 | 男人的好看免费观看在线视频| 99热这里只有是精品在线观看| 亚洲精品一区蜜桃| 麻豆成人av视频| 免费观看在线日韩| 亚洲成人中文字幕在线播放| 一边亲一边摸免费视频| 欧美又色又爽又黄视频| 国产真实伦视频高清在线观看| 在线免费观看的www视频| 亚洲av电影在线观看一区二区三区 | 久99久视频精品免费| 高清毛片免费看| 日日摸夜夜添夜夜添av毛片| 高清午夜精品一区二区三区| 成人三级黄色视频| АⅤ资源中文在线天堂| 日韩av在线大香蕉| 一级黄片播放器| 久久久精品大字幕| 黄色一级大片看看| 亚洲精品日韩在线中文字幕| 美女黄网站色视频| 噜噜噜噜噜久久久久久91| 亚洲av免费高清在线观看| 欧美另类亚洲清纯唯美| 亚洲无线观看免费| 永久免费av网站大全| 久久亚洲精品不卡| av在线天堂中文字幕| 亚洲欧美中文字幕日韩二区| 午夜福利成人在线免费观看| 在线播放无遮挡| 国产免费福利视频在线观看| 国产精品av视频在线免费观看| 乱人视频在线观看| 狂野欧美激情性xxxx在线观看| 春色校园在线视频观看| 老司机福利观看| 天堂av国产一区二区熟女人妻| 国产熟女欧美一区二区| 日韩人妻高清精品专区| 欧美一区二区精品小视频在线| 精品久久久久久久久亚洲| 在线观看av片永久免费下载| 久久久久久伊人网av| 综合色丁香网| 少妇人妻精品综合一区二区| 国产av一区在线观看免费| 天天一区二区日本电影三级| 久久精品人妻少妇| 在线免费观看的www视频| 91精品伊人久久大香线蕉| 久久99热6这里只有精品| 成人欧美大片| 中文天堂在线官网| 丰满人妻一区二区三区视频av| 国产精品久久久久久av不卡| 久久99热6这里只有精品| 久久久欧美国产精品| 日韩av不卡免费在线播放| 天堂√8在线中文| 亚洲精品影视一区二区三区av| 晚上一个人看的免费电影| 免费人成在线观看视频色| 欧美日韩一区二区视频在线观看视频在线 | 身体一侧抽搐| 免费不卡的大黄色大毛片视频在线观看 | 国产午夜福利久久久久久| 亚洲av成人精品一二三区| 视频中文字幕在线观看| a级毛色黄片| 免费一级毛片在线播放高清视频| 日本-黄色视频高清免费观看| 日韩精品有码人妻一区| 国产精品精品国产色婷婷| 午夜爱爱视频在线播放| 一级爰片在线观看| 欧美日韩综合久久久久久| 五月伊人婷婷丁香| 最近中文字幕高清免费大全6| 麻豆乱淫一区二区| 可以在线观看毛片的网站| 三级国产精品欧美在线观看| 久久久成人免费电影| 亚洲国产精品成人综合色| 啦啦啦啦在线视频资源| 在线播放无遮挡| 婷婷色麻豆天堂久久 | 久久久a久久爽久久v久久| 三级男女做爰猛烈吃奶摸视频| 久久亚洲精品不卡| 最近中文字幕高清免费大全6| 免费看日本二区| 亚洲无线观看免费| 波多野结衣高清无吗| 久久久久久国产a免费观看| 国产精品一区www在线观看| 国产精品麻豆人妻色哟哟久久 | 三级经典国产精品| 久久这里只有精品中国| 丝袜美腿在线中文| 欧美高清性xxxxhd video| 看十八女毛片水多多多| 在线观看一区二区三区| 国产精品女同一区二区软件| 久久精品久久久久久久性| 亚洲欧美成人精品一区二区| 日韩av在线免费看完整版不卡| 淫秽高清视频在线观看| 日韩 亚洲 欧美在线| 美女被艹到高潮喷水动态| 精品一区二区三区人妻视频| 国产精品一区二区三区四区久久| 一级毛片我不卡| 精品久久久久久久久久久久久| 波野结衣二区三区在线| 日韩制服骚丝袜av| 亚洲综合精品二区| 一级爰片在线观看| 中文字幕av成人在线电影| 永久免费av网站大全| 国内精品美女久久久久久| 亚洲综合色惰| 最后的刺客免费高清国语| 亚洲天堂国产精品一区在线| 一级黄色大片毛片| 日韩欧美三级三区| 在现免费观看毛片| 日韩一区二区视频免费看| 日韩av在线免费看完整版不卡| 国产精品爽爽va在线观看网站| 韩国av在线不卡| 特级一级黄色大片| 亚洲精品成人久久久久久| 成人午夜高清在线视频| 免费黄网站久久成人精品| 国产白丝娇喘喷水9色精品| 亚洲伊人久久精品综合 | 国产老妇伦熟女老妇高清| 成年免费大片在线观看| 特级一级黄色大片| 最新中文字幕久久久久| videossex国产| 久久人妻av系列| 午夜福利在线观看免费完整高清在| 一边亲一边摸免费视频| 久久久久性生活片| 精品少妇黑人巨大在线播放 | 国产黄色视频一区二区在线观看 | 国产高清视频在线观看网站| 日本欧美国产在线视频| 国产探花在线观看一区二区| 日韩一本色道免费dvd| 舔av片在线| 午夜精品国产一区二区电影 | 26uuu在线亚洲综合色| 国产精品久久久久久av不卡| 99久久成人亚洲精品观看| 日本黄色视频三级网站网址| 九色成人免费人妻av| 免费在线观看成人毛片| 日韩国内少妇激情av| 欧美bdsm另类| videossex国产| 国产成人精品一,二区| ponron亚洲| 综合色丁香网| 国产综合懂色| 精品人妻偷拍中文字幕| 欧美一区二区亚洲| 亚洲国产成人一精品久久久| 99久久九九国产精品国产免费| 亚洲av男天堂| 亚洲欧美一区二区三区国产| 亚洲av中文字字幕乱码综合| 亚洲乱码一区二区免费版| 欧美激情在线99| 99久久成人亚洲精品观看| 国产免费福利视频在线观看| 成人特级av手机在线观看| 久久精品国产99精品国产亚洲性色| 大香蕉久久网| 色吧在线观看| 水蜜桃什么品种好| 国产伦精品一区二区三区四那| 亚洲精品,欧美精品| 亚洲aⅴ乱码一区二区在线播放| 在线播放国产精品三级| 午夜福利在线观看免费完整高清在| 国产精品三级大全| 色视频www国产| 国产精品.久久久| 亚洲av.av天堂| 简卡轻食公司| 麻豆成人av视频| 国产精品av视频在线免费观看| 午夜激情福利司机影院| www.av在线官网国产| av免费在线看不卡| 国产美女午夜福利| 欧美人与善性xxx| 搡女人真爽免费视频火全软件| 久久久色成人| 成人午夜精彩视频在线观看| 欧美日韩精品成人综合77777| 久久久色成人| 国产午夜精品论理片| 26uuu在线亚洲综合色| 午夜免费激情av| 久久亚洲国产成人精品v| 亚洲精品日韩av片在线观看| 嫩草影院入口| 国产淫语在线视频| 亚洲成人中文字幕在线播放| 日韩亚洲欧美综合| 久久久欧美国产精品| 男人舔女人下体高潮全视频| 最近手机中文字幕大全| a级一级毛片免费在线观看| 国产 一区精品| 能在线免费看毛片的网站| 成人一区二区视频在线观看| www日本黄色视频网| 美女cb高潮喷水在线观看| 国产午夜精品论理片| 美女黄网站色视频| 少妇熟女欧美另类| 国产黄片美女视频| 国产欧美另类精品又又久久亚洲欧美| 国产乱来视频区| 国产精品久久久久久久电影| 噜噜噜噜噜久久久久久91| 亚洲久久久久久中文字幕| 精品人妻偷拍中文字幕| 男女视频在线观看网站免费| 青春草视频在线免费观看| 亚洲av男天堂| 干丝袜人妻中文字幕| 黄色一级大片看看| АⅤ资源中文在线天堂| 国产乱来视频区| 免费观看在线日韩| 超碰av人人做人人爽久久| 午夜精品在线福利| 欧美日本视频| 国产精品一及| 超碰97精品在线观看| 国产精品久久久久久久电影| 听说在线观看完整版免费高清| 国产久久久一区二区三区| 最近2019中文字幕mv第一页| 欧美高清成人免费视频www| av播播在线观看一区| 午夜福利视频1000在线观看| 亚洲,欧美,日韩| 亚洲av男天堂| 国产精品一二三区在线看| 波野结衣二区三区在线| 久久精品综合一区二区三区| 国产91av在线免费观看| 男女啪啪激烈高潮av片| 视频中文字幕在线观看| 日本一本二区三区精品| 国产真实伦视频高清在线观看| 国产淫片久久久久久久久| 在线播放无遮挡| 国产私拍福利视频在线观看| 建设人人有责人人尽责人人享有的 | 免费观看性生交大片5| 亚洲丝袜综合中文字幕| 日韩制服骚丝袜av| 欧美丝袜亚洲另类| 国产高清不卡午夜福利| 日韩强制内射视频| 久久精品国产鲁丝片午夜精品| 国产三级在线视频| 在现免费观看毛片| 国产精品99久久久久久久久| 国产成人免费观看mmmm| 欧美3d第一页| 国产成人精品久久久久久| 天美传媒精品一区二区| 亚洲三级黄色毛片| 久久这里只有精品中国| a级毛色黄片| 精品不卡国产一区二区三区| 少妇裸体淫交视频免费看高清| 成人美女网站在线观看视频| 看非洲黑人一级黄片| 国产人妻一区二区三区在| 欧美最新免费一区二区三区| 91久久精品国产一区二区成人| 日韩成人伦理影院| 少妇的逼好多水| 级片在线观看| 久久国产乱子免费精品| 亚洲国产色片| 美女脱内裤让男人舔精品视频| 久久精品久久精品一区二区三区| 免费一级毛片在线播放高清视频| 丰满人妻一区二区三区视频av| 亚洲欧美精品综合久久99| 午夜日本视频在线| 亚洲精品国产av成人精品| 国产免费男女视频| 一级毛片电影观看 | 亚洲国产精品合色在线| 有码 亚洲区| 国产 一区 欧美 日韩| 少妇丰满av| 国产视频内射| 日本免费一区二区三区高清不卡| av国产久精品久网站免费入址| 欧美激情国产日韩精品一区| 亚洲欧洲日产国产| 久久久久久久亚洲中文字幕| 日本色播在线视频| 久久久久性生活片| 午夜福利网站1000一区二区三区| 插阴视频在线观看视频| 国产精品一区二区三区四区免费观看| 亚洲av电影在线观看一区二区三区 | 美女脱内裤让男人舔精品视频| 最近中文字幕高清免费大全6| 欧美变态另类bdsm刘玥| 亚洲熟妇中文字幕五十中出| 久久久久久久久久久丰满| 最新中文字幕久久久久| 成年免费大片在线观看| 22中文网久久字幕| 1024手机看黄色片| 欧美三级亚洲精品| 国产精品嫩草影院av在线观看| 人体艺术视频欧美日本| 亚洲婷婷狠狠爱综合网| 91狼人影院| av在线天堂中文字幕| 赤兔流量卡办理| 国产成年人精品一区二区| 国产综合懂色| 国产伦精品一区二区三区视频9| 色视频www国产| 只有这里有精品99| 热99re8久久精品国产| 久久久久久九九精品二区国产| 男插女下体视频免费在线播放| 黄色日韩在线| 少妇裸体淫交视频免费看高清| 中文字幕av在线有码专区| 国产亚洲91精品色在线| 午夜久久久久精精品| 亚洲av.av天堂| 国产久久久一区二区三区| 久久亚洲精品不卡| 欧美极品一区二区三区四区| 亚洲国产精品sss在线观看| 建设人人有责人人尽责人人享有的 | 99久国产av精品国产电影| 小说图片视频综合网站| 六月丁香七月| 夫妻性生交免费视频一级片| 99久久无色码亚洲精品果冻| 伦理电影大哥的女人| 天天躁日日操中文字幕| 白带黄色成豆腐渣| 国产精品久久久久久av不卡| 岛国毛片在线播放| 国产黄片视频在线免费观看| 国产免费男女视频| 观看美女的网站| 熟妇人妻久久中文字幕3abv| 噜噜噜噜噜久久久久久91| 亚洲精华国产精华液的使用体验| 久久久精品94久久精品| 午夜免费男女啪啪视频观看| 国国产精品蜜臀av免费| 赤兔流量卡办理| or卡值多少钱| 97超碰精品成人国产| 久久久久久久午夜电影| 成人午夜精彩视频在线观看| a级毛片免费高清观看在线播放| 中文字幕av在线有码专区| 成人特级av手机在线观看| 成人午夜精彩视频在线观看| 床上黄色一级片| a级一级毛片免费在线观看| 久久久久久大精品| 小蜜桃在线观看免费完整版高清| 国产精品一区二区三区四区久久| 亚洲人成网站高清观看| 丰满人妻一区二区三区视频av| 老女人水多毛片| 亚洲18禁久久av| 免费观看精品视频网站| 国产一区二区在线av高清观看| 秋霞伦理黄片| 国产精品爽爽va在线观看网站| 91狼人影院| 日日干狠狠操夜夜爽| 国产白丝娇喘喷水9色精品| videossex国产| 亚洲电影在线观看av| 欧美色视频一区免费| 人妻夜夜爽99麻豆av| 国产黄色视频一区二区在线观看 | 最近手机中文字幕大全| 国产成人精品久久久久久| 亚洲国产色片| 又粗又爽又猛毛片免费看| 日韩一区二区三区影片| 黄片无遮挡物在线观看| 特大巨黑吊av在线直播| 国内揄拍国产精品人妻在线| 国产av一区在线观看免费| 又粗又硬又长又爽又黄的视频| 秋霞伦理黄片| 激情 狠狠 欧美|