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

    灰線性馬爾科夫模型在建筑物變形監(jiān)測(cè)中的應(yīng)用

    2016-08-02 03:38:34李克昭李志偉趙磊杰
    測(cè)繪工程 2016年10期
    關(guān)鍵詞:馬爾科夫預(yù)測(cè)值殘差

    李克昭,李志偉,趙磊杰

    (1.河南理工大學(xué) 測(cè)繪與國(guó)土信息工程學(xué)院,河南 焦作 454000;2.北斗導(dǎo)航應(yīng)用技術(shù)協(xié)同創(chuàng)新中心,河南 鄭州 450052)

    ?

    灰線性馬爾科夫模型在建筑物變形監(jiān)測(cè)中的應(yīng)用

    李克昭1,2,李志偉1,趙磊杰1

    (1.河南理工大學(xué) 測(cè)繪與國(guó)土信息工程學(xué)院,河南 焦作 454000;2.北斗導(dǎo)航應(yīng)用技術(shù)協(xié)同創(chuàng)新中心,河南 鄭州 450052)

    摘要:針對(duì)傳統(tǒng)灰色GM(1,1)預(yù)測(cè)模型在建筑物變形監(jiān)測(cè)預(yù)報(bào)中的擬合精度較差、預(yù)測(cè)精度較低和預(yù)測(cè)時(shí)間較短的問(wèn)題,文中以傳統(tǒng)GM(1,1)、線性回歸和馬爾科夫模型為理論基礎(chǔ),構(gòu)建了灰線性馬爾科夫預(yù)測(cè)模型,并結(jié)合某建筑物變形監(jiān)測(cè)的觀測(cè)數(shù)據(jù),運(yùn)用新陳代謝的計(jì)算模式進(jìn)行預(yù)測(cè)。結(jié)果表明,灰線性馬爾科夫預(yù)測(cè)模型的擬合精度和預(yù)測(cè)精度優(yōu)于單一的灰色GM(1,1)預(yù)測(cè)模型和線性回歸預(yù)測(cè)模型,灰線性馬爾科夫預(yù)測(cè)模型具有預(yù)測(cè)精度高、預(yù)測(cè)時(shí)間長(zhǎng)和穩(wěn)定性高的優(yōu)勢(shì)。

    關(guān)鍵詞:GM(1,1)模型;線性回歸模型;馬爾科夫模型;新陳代謝

    對(duì)建筑物的變形監(jiān)測(cè)分析與預(yù)報(bào)是非常必要的。變形監(jiān)測(cè)的預(yù)報(bào)可用于指導(dǎo)建筑施工、確保施工質(zhì)量和得到建筑物變形的先驗(yàn)信息。灰色系統(tǒng)理論是研究數(shù)據(jù)貧乏、資料較少、不確定性問(wèn)題的理論,非常適合于建筑物的變形監(jiān)測(cè)預(yù)報(bào)工作。

    GM(1,1)模型是灰色理論的最簡(jiǎn)單的預(yù)測(cè)模型之一,它打破了回歸分析和概率統(tǒng)計(jì)的局限性,以灰色生成函數(shù)為基礎(chǔ),以微分?jǐn)M合為核心的一種建模方法。但是GM(1,1)模型在構(gòu)建模型過(guò)程中受到隨機(jī)數(shù)據(jù)擾動(dòng)影響較大,構(gòu)建的模型穩(wěn)定性較差。近年來(lái),很多學(xué)者就針對(duì)GM(1,1)模型的初值確定[1]、GM(1,1)背景值的構(gòu)建[2]和GM(1,1)模型殘差的修正[3-4]方面做了很多方法的嘗試。但是對(duì)單一模型進(jìn)行改進(jìn),對(duì)于提高預(yù)測(cè)精度相對(duì)來(lái)說(shuō)是緩慢的,同時(shí)組合幾個(gè)互補(bǔ)的預(yù)測(cè)模型可使提高預(yù)測(cè)精度達(dá)到立竿見(jiàn)影的效果。

    線性回歸預(yù)測(cè)模型[5-7]根據(jù)事物成長(zhǎng)的規(guī)律性、事物發(fā)展的連續(xù)性以及事物因果的相關(guān)性對(duì)短時(shí)期的預(yù)測(cè)能夠取得非常好的效果,但對(duì)長(zhǎng)期預(yù)測(cè)的效果并不明顯。馬爾科夫預(yù)測(cè)模型的轉(zhuǎn)移概率矩陣可以有效反映出隨機(jī)數(shù)據(jù)的波動(dòng)程度,很大程度上彌補(bǔ)了GM(1,1)和線性回歸預(yù)測(cè)模型的局限性[8]。3種預(yù)測(cè)模型組合成灰線性馬爾科夫預(yù)測(cè)模型不僅提高了預(yù)測(cè)模型的精度和適用范圍,同時(shí)還能保證預(yù)測(cè)模型的穩(wěn)定性。本文以傳統(tǒng)GM(1,1)、線性回歸和馬爾科夫預(yù)測(cè)模型為理論基礎(chǔ),將提出灰線性馬爾科夫預(yù)測(cè)模型,并結(jié)合某建筑物的變形監(jiān)測(cè)觀測(cè)資料,對(duì)組合預(yù)測(cè)模型進(jìn)行驗(yàn)證和分析。

    1灰線性馬爾科夫組合預(yù)測(cè)模型的建立

    1.1灰線性GM(1,1)模型

    假設(shè)一組原始數(shù)據(jù)序列,記為X(0),

    (1)

    對(duì)數(shù)據(jù)X(0)序列進(jìn)行一次累加,生成的新數(shù)據(jù)序列記為X(1),

    (2)

    (3)

    式中:a為發(fā)展系數(shù),b為灰作用量,其形式可以記為

    (4)

    (5)

    將式(5)進(jìn)行累減計(jì)算得到灰線性GM(1,1)模型的擬合值和預(yù)測(cè)值。從式(5)中不難看出:當(dāng)m1=0時(shí),預(yù)測(cè)方程為線性回歸方程;當(dāng)m2=0時(shí),預(yù)測(cè)方程為傳統(tǒng)GM(1,1)預(yù)測(cè)方程。因此,灰線性GM(1,1)預(yù)測(cè)模型繼承了傳統(tǒng)GM(1,1)和線性回歸的優(yōu)點(diǎn)。

    1.2灰線性馬爾科夫預(yù)測(cè)模型

    1)計(jì)算波動(dòng)指數(shù)序列

    (6)

    2)劃分馬爾科夫狀態(tài)

    (7)

    式(7)表示為第k對(duì)象的波動(dòng)指數(shù)處于第i種狀態(tài)Ei;a1i,a2i分別表示為狀態(tài)Ei的上下界。因此,總的狀態(tài)集合表示為

    3)構(gòu)造狀態(tài)轉(zhuǎn)移概率矩陣

    (8)

    式中:mij(n)為狀態(tài)Ei經(jīng)過(guò)n步轉(zhuǎn)移到達(dá)狀態(tài)Ej的次數(shù),Mi為Ei出現(xiàn)的次數(shù)。由n步轉(zhuǎn)移概率元素Pij(n)構(gòu)成n步轉(zhuǎn)移概率矩陣P(n)。

    4)構(gòu)造未來(lái)狀態(tài)概率矩陣

    (9)

    式中選取離預(yù)測(cè)目標(biāo)最近的s個(gè)原始對(duì)象,按照從近到遠(yuǎn)的順序,所需的轉(zhuǎn)移步數(shù)分別為1,2,…,s,分別以各個(gè)對(duì)象所對(duì)應(yīng)的狀態(tài)為初始狀態(tài),在n步轉(zhuǎn)移概率矩陣P(n)中取各自所對(duì)應(yīng)的行向量Pi(n)=[Pi1(n),Pi2(n),…,Pis(n)],i∈s,從而構(gòu)成新的概率矩陣。

    5)得出預(yù)測(cè)方程。

    (10)

    2灰線性馬爾科夫預(yù)測(cè)模型精度評(píng)定

    表1 預(yù)測(cè)模型精度評(píng)價(jià)[8]

    3實(shí)例計(jì)算與結(jié)果分析

    本文采用文獻(xiàn)[7]中的居民樓變形監(jiān)測(cè)工程數(shù)據(jù),該工程為廣州市某地產(chǎn)開(kāi)發(fā)公司開(kāi)發(fā)的新建居民樓。按照二等水準(zhǔn)的測(cè)量規(guī)范,獲得1#建筑物CJ1號(hào)點(diǎn)13期沉降累計(jì)觀測(cè)數(shù)據(jù)。本文利用前8期沉降累計(jì)數(shù)據(jù)進(jìn)行建模,分別用GM(1,1)模型、灰色線性回歸模型、灰線性回歸預(yù)測(cè)第9~13期沉降累計(jì)數(shù)據(jù)。CJ1號(hào)點(diǎn)的觀測(cè)數(shù)據(jù)如表2所示。

    表2 CJ1點(diǎn)的實(shí)測(cè)累計(jì)沉降量

    3.1灰線性馬爾科夫模型數(shù)據(jù)計(jì)算

    文獻(xiàn)[5-7]中表明,當(dāng)建模過(guò)程中的建模期數(shù)n=8,參數(shù)m=1,2,3,4時(shí),模型擬合精度最優(yōu)。參考文獻(xiàn)[12]的新陳代謝計(jì)算模式,用表2的前8期數(shù)據(jù)構(gòu)建灰線性馬爾科夫模型,預(yù)測(cè)第9期數(shù)據(jù)為例,計(jì)算過(guò)程如下:

    1)利用前8期數(shù)據(jù)構(gòu)建灰線性GM(1,1)模型

    X(0)=[2.82,3.51,4.59,5.87,6.89,

    8.37,10.99,12.59].

    2)灰線性GM(1,1)模型建模生成模型擬合序列、殘差序列和波動(dòng)指數(shù)序列,同時(shí)計(jì)算可得第9期的預(yù)測(cè)值,如表3所示。

    表3 灰線性GM(1,1)的計(jì)算結(jié)果

    3)劃分馬爾科夫模型狀態(tài)。根據(jù)表3中的波動(dòng)指數(shù)劃分為4種狀態(tài):E1:[-4,-2],E2:[-2,0],E3:[0,2],E4:[2,6]。

    模型值所處的狀態(tài),如表4所示。

    表4 模型值狀態(tài)的劃分

    4)計(jì)算轉(zhuǎn)移概率矩陣

    5)計(jì)算第9期數(shù)據(jù)所處的狀態(tài)。第9期狀態(tài)的預(yù)測(cè)計(jì)算表如表5所示。

    表5 第9期狀態(tài)預(yù)測(cè)表

    從表5中可得,第9期的預(yù)測(cè)狀態(tài)為E4。

    6)計(jì)算第9期的預(yù)測(cè)值。第9期的預(yù)測(cè)值為

    3.2灰線性馬爾科夫模型數(shù)據(jù)結(jié)果分析

    利用Matlab7.0軟件為平臺(tái),按照上述計(jì)算步驟編寫程序,計(jì)算第9~13期的預(yù)測(cè)值,即:利用前8期觀測(cè)數(shù)據(jù)建模,得到第9期的預(yù)測(cè)值;然后去掉建模數(shù)據(jù)中第1期數(shù)據(jù),加入第9期的預(yù)測(cè)值重新建模,計(jì)算第10期的預(yù)測(cè)值;依次類推,計(jì)算出第9~13期的預(yù)測(cè)值;最后,求取新陳代謝過(guò)程中所產(chǎn)生的擬合值和預(yù)測(cè)值的平均值,作為灰線性馬爾科夫模型的擬合值和預(yù)測(cè)值。其計(jì)算結(jié)果和分析如下:

    1)GM(1,1)、灰線性GM(1,1)、灰線性馬爾科夫3種模型的擬合值和預(yù)測(cè)值見(jiàn)表6、表7。

    表6 3種預(yù)測(cè)模型擬合值結(jié)果 mm

    從表6得出,灰線性馬爾科夫的模型擬合殘差中誤差0.030 2mm,灰線性GM(1,1)模型的擬合殘差中誤差0.233 5mm,傳統(tǒng)的GM(1,1)模型的擬合殘差中誤差0.569 1mm?;揖€性馬爾科夫組合模型的擬合精度優(yōu)于前兩種預(yù)測(cè)模型。

    表7 3種預(yù)測(cè)模型預(yù)測(cè)值結(jié)果 mm

    從表7中得到,灰線性馬爾科夫預(yù)測(cè)模型的預(yù)測(cè)值殘差中誤差1.181 5mm,灰線性GM(1,1)預(yù)測(cè)模型的預(yù)測(cè)值殘差中誤差1.676 9mm,傳統(tǒng)GM(1,1)預(yù)測(cè)模型的預(yù)測(cè)值殘差中誤差3.354 3mm?;揖€性馬爾科夫組合模型的預(yù)測(cè)精度優(yōu)于前兩種預(yù)測(cè)模型。

    2)GM(1,1)、灰線性GM(1,1)、灰線性馬爾科夫3種模型的相對(duì)誤差見(jiàn)表8。

    相對(duì)誤差更能直觀反映出模型的相對(duì)精度,區(qū)分出模型精度的優(yōu)劣。

    從表8和圖1中可以看出,無(wú)論是擬合值還是預(yù)測(cè)值,灰線性馬爾科夫預(yù)測(cè)模型均優(yōu)于灰線性GM(1,1)預(yù)測(cè)模型和傳統(tǒng)GM(1,1)預(yù)測(cè)模型。因此,灰線性馬爾科夫預(yù)測(cè)模型更加準(zhǔn)確,外推的預(yù)測(cè)值更多,外推的預(yù)測(cè)值更可靠,模型更加穩(wěn)定。

    3)GM(1,1)、灰線性GM(1,1)、灰線性馬爾科夫3種模型的精度評(píng)定見(jiàn)表9。

    表8 3種預(yù)測(cè)模型的相對(duì)誤差結(jié)果

    圖1 3種預(yù)測(cè)模型相對(duì)誤差比較

    傳統(tǒng)GM(1,1)灰線性GM(1,1)灰線性馬爾科夫P111C0.07450.07070.0149

    經(jīng)計(jì)算,傳統(tǒng)GM(1,1)預(yù)測(cè)模型的小殘差概率P=1,方差比C=0.074 5;灰線性GM(1,1)預(yù)測(cè)模型的小殘差概率P=1,方差比C=0.070 7;灰線性馬爾科夫預(yù)測(cè)模型的小殘差概率P=1,方差比C=0.014 9。由于所選取的數(shù)據(jù)變化比較平滑,所以3種預(yù)測(cè)模型精度結(jié)果皆為優(yōu)。但是從方差比可以看出,灰線性馬爾科夫預(yù)測(cè)模型精度優(yōu)于前兩種預(yù)測(cè)模型。

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

    本文綜合傳統(tǒng)GM(1,1)模型、線性回歸模型和馬爾科夫模型的優(yōu)點(diǎn),構(gòu)建了灰線性馬爾科夫預(yù)測(cè)模型,運(yùn)用新陳代謝的計(jì)算模式對(duì)建筑物變形監(jiān)測(cè)進(jìn)行預(yù)測(cè)。得出的結(jié)論為:灰線性馬爾科夫預(yù)測(cè)模型無(wú)論是擬合精度還是預(yù)測(cè)精度都要優(yōu)于單一的灰色GM(1,1)預(yù)測(cè)模型和線性回歸預(yù)測(cè)模型,灰線性馬爾科夫預(yù)測(cè)模型具有預(yù)測(cè)精度高、預(yù)測(cè)時(shí)間長(zhǎng)和模型穩(wěn)定性高的優(yōu)勢(shì)。

    對(duì)于變形監(jiān)測(cè)工作中實(shí)測(cè)的觀測(cè)數(shù)據(jù)建立合適的預(yù)測(cè)模型。在滿足工程應(yīng)用的精度要求下,預(yù)測(cè)模型需要的觀測(cè)資料較少,減輕了大量的外業(yè)實(shí)測(cè)工作,提高了工作效率,同時(shí)提供了可靠的參考資料。本文建立的灰線性馬爾科夫預(yù)測(cè)模型,為了保證預(yù)測(cè)模型的可靠性,建議建模數(shù)據(jù)不要低于6期實(shí)測(cè)數(shù)據(jù),同時(shí)建議連續(xù)外推的預(yù)測(cè)值不要超過(guò)5期為宜。

    參考文獻(xiàn):

    [1]張彬,西桂權(quán).基于背景值和邊值修正的GM(1,1)模型優(yōu)化[J].系統(tǒng)工程理論與實(shí)踐,2013,33(3):682-688.

    [2]袁豹,岳東杰,李成任.基于總體最小二乘的改進(jìn)GM(1,1)模型及其在建筑物沉降預(yù)測(cè)中應(yīng)用[J].測(cè)繪工程,2013,22(3):52-55.

    [3]李秀珍,孔紀(jì)名,王成華.灰色GM(1,1)殘差修正模型在滑坡預(yù)測(cè)中的對(duì)比應(yīng)用[J].山地學(xué)報(bào),2007,25(6):741-746.

    [4]郭蘭蘭,鄒志紅,安巖.基于殘差修正的GM(1,1)模型在水質(zhì)預(yù)測(cè)中的應(yīng)用[J].數(shù)學(xué)的實(shí)踐與認(rèn)識(shí),2014,44(19):176-181.

    [5]高寧,崔希民,高彩云.高層建筑物沉降變形的灰線性預(yù)測(cè)[J].測(cè)繪科學(xué),2012,37(3):96-98.

    [6]陳懷爭(zhēng),王勇.灰色線性回歸組合模型在變形監(jiān)測(cè)數(shù)據(jù)處理和預(yù)測(cè)中的應(yīng)用[J].資源與產(chǎn)業(yè),2011,13(增1):165-167.

    [7]張俊中,雷偉偉,王睿,等.灰線性回歸模型在建筑物沉降分析中的應(yīng)用[J].河南科學(xué),2015,33(3):416-420.

    [8]王磊,武術(shù)靜,李長(zhǎng)青.灰色馬爾科夫模型對(duì)煤自燃發(fā)火預(yù)測(cè)的研究[J].河南理工大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,34(1):35-39.

    [9]楊燦,楊敏華,汪奇松.基于灰色馬爾科夫鏈預(yù)測(cè)系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)[J].測(cè)繪通報(bào),2009(6):182-183.

    [10] 沈哲輝,黃騰,唐佑輝.灰色-馬爾科夫模型在大壩內(nèi)部變形預(yù)測(cè)中的應(yīng)用[J].測(cè)繪工程,2015,24(2):69-74.

    [11] 楊錦偉,孫寶磊.基于灰色馬爾科夫模型在平頂山市空氣污染物濃度預(yù)測(cè)[J].數(shù)學(xué)的實(shí)踐與認(rèn)識(shí),2014,44(2):64-70.

    [12] 岳仁賓,騰德貴,胡波,等.灰色模型在深基坑變形監(jiān)測(cè)中的應(yīng)用研究[J].測(cè)繪通報(bào),2014(S1):85-87.

    [責(zé)任編輯:劉文霞]

    DOI:10.19349/j.cnki.issn1006-7949.2016.10.002

    收稿日期:2015-06-26

    基金項(xiàng)目:國(guó)家自然科學(xué)基金資助項(xiàng)目(41202245;41272373);河南理工大學(xué)骨干教師資助項(xiàng)目(72105/090)

    作者簡(jiǎn)介:李克昭(1977-),男,副教授,博士.

    中圖分類號(hào):TV698.1

    文獻(xiàn)標(biāo)識(shí)碼:A

    文章編號(hào):1006-7949(2016)10-0005-05

    Application of grey line Markov model to the deformation monitoring of buildings

    LI Kezhao1,2,LI Zhiwei1,ZHAO Leijie1

    (1.School of Surveying and Landing Information Engineering,Henan Polytechnic University,Jiaozuo 454000,China;2.Collaborative Innovation Center of BDS Research Application,Zhengzhou 450052,China)

    Abstract:Traditional gray GM (1,1) predicting model has the problems of poor fitting accuracy,lower prediction accuracy and shorter prediction time in deformation monitoring and forecasting of buildings.In this paper,a combination of the traditional GM (1,1) model,linear regression and Markov model is constructed of grey linear Markov model.Combined with observations data of the deformation monitoring of buildings,the metabolism computing model is used to predict.The results show that:the fitting accuracy and model prediction accuracy of the linear Markov model gray model are better than a single gray GM (1,1) forecast model and linear regression forecasting model.Gray linear Markov model has the advantages of high accuracy,long time and high stability of prediction.

    Key words:GM(1,1) model;linear regression models;Markov model;metabolism

    猜你喜歡
    馬爾科夫預(yù)測(cè)值殘差
    IMF上調(diào)今年全球經(jīng)濟(jì)增長(zhǎng)預(yù)期
    企業(yè)界(2024年8期)2024-07-05 10:59:04
    基于雙向GRU與殘差擬合的車輛跟馳建模
    加拿大農(nóng)業(yè)部下調(diào)2021/22年度油菜籽和小麥產(chǎn)量預(yù)測(cè)值
    ±800kV直流輸電工程合成電場(chǎng)夏季實(shí)測(cè)值與預(yù)測(cè)值比對(duì)分析
    基于疊加馬爾科夫鏈的邊坡位移預(yù)測(cè)研究
    基于殘差學(xué)習(xí)的自適應(yīng)無(wú)人機(jī)目標(biāo)跟蹤算法
    基于改進(jìn)的灰色-馬爾科夫模型在風(fēng)機(jī)沉降中的應(yīng)用
    法電再次修訂2020年核發(fā)電量預(yù)測(cè)值
    基于遞歸殘差網(wǎng)絡(luò)的圖像超分辨率重建
    馬爾科夫鏈在教學(xué)評(píng)價(jià)中的應(yīng)用
    久久久久视频综合| 中国三级夫妇交换| 欧美xxⅹ黑人| 99精国产麻豆久久婷婷| 亚洲精华国产精华液的使用体验| 免费黄色在线免费观看| av在线播放精品| 免费久久久久久久精品成人欧美视频| 亚洲av中文av极速乱| 精品少妇黑人巨大在线播放| 一区福利在线观看| 久久久久久人妻| 曰老女人黄片| 亚洲婷婷狠狠爱综合网| 天堂8中文在线网| 成人国语在线视频| 日韩中字成人| 91aial.com中文字幕在线观看| 亚洲第一区二区三区不卡| 日韩视频在线欧美| 如日韩欧美国产精品一区二区三区| 一本色道久久久久久精品综合| 欧美av亚洲av综合av国产av | 91成人精品电影| 欧美日韩av久久| 欧美日韩精品成人综合77777| 日韩av不卡免费在线播放| 久久久国产欧美日韩av| 精品久久久精品久久久| 色哟哟·www| 69精品国产乱码久久久| 久久精品国产亚洲av涩爱| 自线自在国产av| 一区二区三区四区激情视频| 这个男人来自地球电影免费观看 | 国产精品嫩草影院av在线观看| 日韩伦理黄色片| 女人久久www免费人成看片| 精品国产一区二区三区四区第35| 18在线观看网站| 人人妻人人爽人人添夜夜欢视频| 欧美变态另类bdsm刘玥| 亚洲精品一二三| 亚洲精品美女久久久久99蜜臀 | 亚洲国产欧美网| 制服人妻中文乱码| 亚洲精品乱久久久久久| 亚洲综合色惰| 香蕉国产在线看| 天堂中文最新版在线下载| 久久精品国产亚洲av高清一级| 国产在线免费精品| 五月开心婷婷网| 国产精品二区激情视频| 丝袜喷水一区| 男人舔女人的私密视频| 免费黄色在线免费观看| 欧美 亚洲 国产 日韩一| 最近手机中文字幕大全| 亚洲成人手机| 午夜免费鲁丝| 国产综合精华液| 制服人妻中文乱码| 欧美国产精品一级二级三级| 丰满乱子伦码专区| 丰满乱子伦码专区| av片东京热男人的天堂| 欧美精品国产亚洲| 午夜久久久在线观看| 中文精品一卡2卡3卡4更新| 色婷婷av一区二区三区视频| 午夜av观看不卡| a级片在线免费高清观看视频| 男女午夜视频在线观看| tube8黄色片| 久久精品夜色国产| 一级毛片我不卡| 欧美精品亚洲一区二区| 一边摸一边做爽爽视频免费| 99久久精品国产国产毛片| 国产精品久久久久久精品电影小说| 一边摸一边做爽爽视频免费| 成年av动漫网址| 天天躁狠狠躁夜夜躁狠狠躁| 精品视频人人做人人爽| 国产探花极品一区二区| videosex国产| 男女免费视频国产| 亚洲一区二区三区欧美精品| 日韩中字成人| 大码成人一级视频| 成人漫画全彩无遮挡| 久久人人爽av亚洲精品天堂| 一区二区三区乱码不卡18| 欧美日韩精品网址| 国产无遮挡羞羞视频在线观看| 国产无遮挡羞羞视频在线观看| 欧美 日韩 精品 国产| 天天躁夜夜躁狠狠躁躁| 国产淫语在线视频| 深夜精品福利| 国产有黄有色有爽视频| 永久免费av网站大全| 汤姆久久久久久久影院中文字幕| 欧美人与性动交α欧美精品济南到 | 丝袜美足系列| 国产黄色视频一区二区在线观看| 飞空精品影院首页| 巨乳人妻的诱惑在线观看| 亚洲第一青青草原| 日本黄色日本黄色录像| 在线天堂最新版资源| 亚洲一码二码三码区别大吗| 这个男人来自地球电影免费观看 | 久久亚洲国产成人精品v| 日韩电影二区| 捣出白浆h1v1| 国产免费福利视频在线观看| 欧美人与性动交α欧美精品济南到 | 免费观看a级毛片全部| 国产精品秋霞免费鲁丝片| 亚洲国产欧美在线一区| 综合色丁香网| 一本大道久久a久久精品| 激情视频va一区二区三区| 婷婷色av中文字幕| 老熟女久久久| 久久av网站| 90打野战视频偷拍视频| 香蕉丝袜av| 免费观看a级毛片全部| 少妇人妻久久综合中文| 激情五月婷婷亚洲| 好男人视频免费观看在线| 高清在线视频一区二区三区| 黑人巨大精品欧美一区二区蜜桃| 熟女av电影| 天堂中文最新版在线下载| 国产精品无大码| 国产成人精品婷婷| 国产av国产精品国产| 国产日韩欧美视频二区| 国产亚洲精品第一综合不卡| av视频免费观看在线观看| 深夜精品福利| 大话2 男鬼变身卡| 午夜av观看不卡| 国产激情久久老熟女| 91精品国产国语对白视频| 99九九在线精品视频| 侵犯人妻中文字幕一二三四区| av线在线观看网站| 男人添女人高潮全过程视频| 看十八女毛片水多多多| 男女边摸边吃奶| 精品久久久久久电影网| 有码 亚洲区| 色播在线永久视频| 久久久久久免费高清国产稀缺| 免费少妇av软件| www.av在线官网国产| 伊人亚洲综合成人网| 日韩一卡2卡3卡4卡2021年| 国产人伦9x9x在线观看 | 亚洲成av片中文字幕在线观看 | 美女大奶头黄色视频| 观看av在线不卡| 亚洲av男天堂| 丰满迷人的少妇在线观看| 中文字幕精品免费在线观看视频| 亚洲av综合色区一区| 国产免费又黄又爽又色| 在线精品无人区一区二区三| 国产福利在线免费观看视频| 一区二区三区精品91| 日本猛色少妇xxxxx猛交久久| 久久青草综合色| 曰老女人黄片| 久久精品国产鲁丝片午夜精品| 国产精品二区激情视频| 免费观看av网站的网址| 免费高清在线观看视频在线观看| 成年女人毛片免费观看观看9 | 成人国语在线视频| av在线播放精品| 老司机亚洲免费影院| 人妻系列 视频| 久久国产精品男人的天堂亚洲| 久久精品国产亚洲av涩爱| 男的添女的下面高潮视频| 青草久久国产| 欧美国产精品一级二级三级| 夜夜骑夜夜射夜夜干| 国产成人欧美| kizo精华| 一级毛片黄色毛片免费观看视频| 久久久亚洲精品成人影院| 国产麻豆69| 我的亚洲天堂| 国产老妇伦熟女老妇高清| av线在线观看网站| 国产成人一区二区在线| 你懂的网址亚洲精品在线观看| 亚洲天堂av无毛| 亚洲av福利一区| 少妇猛男粗大的猛烈进出视频| 久久精品夜色国产| 九九爱精品视频在线观看| 精品一区二区三卡| videosex国产| 日韩av在线免费看完整版不卡| 女人精品久久久久毛片| 久久99热这里只频精品6学生| 黄色视频在线播放观看不卡| 久久久久久久国产电影| 亚洲视频免费观看视频| 肉色欧美久久久久久久蜜桃| 26uuu在线亚洲综合色| 国产高清国产精品国产三级| www日本在线高清视频| 成年av动漫网址| 亚洲中文av在线| 深夜精品福利| 岛国毛片在线播放| 男女啪啪激烈高潮av片| 中文字幕人妻丝袜制服| 香蕉国产在线看| 国产男人的电影天堂91| 久久国产亚洲av麻豆专区| 亚洲欧美日韩另类电影网站| 免费黄频网站在线观看国产| 久久av网站| 久久久久久久久久久久大奶| 亚洲国产av影院在线观看| 亚洲av电影在线观看一区二区三区| 美女国产高潮福利片在线看| 丝袜喷水一区| 国产免费现黄频在线看| 日本色播在线视频| 成人漫画全彩无遮挡| 亚洲国产精品成人久久小说| 波多野结衣av一区二区av| 欧美激情 高清一区二区三区| 国产黄色免费在线视频| 中国国产av一级| 国产黄色免费在线视频| 亚洲综合色惰| 国产成人a∨麻豆精品| 久久人人97超碰香蕉20202| 久久久久久久久久久免费av| 一区二区av电影网| 人妻人人澡人人爽人人| 久久久a久久爽久久v久久| 欧美激情极品国产一区二区三区| 久久99精品国语久久久| 国产日韩一区二区三区精品不卡| 97精品久久久久久久久久精品| 欧美变态另类bdsm刘玥| 一级片免费观看大全| √禁漫天堂资源中文www| 丝袜脚勾引网站| 欧美变态另类bdsm刘玥| 国产精品无大码| av天堂久久9| 99热全是精品| 欧美黄色片欧美黄色片| 亚洲av国产av综合av卡| 久久综合国产亚洲精品| 亚洲精品在线美女| 亚洲精品美女久久久久99蜜臀 | 亚洲精品在线美女| kizo精华| 国产黄频视频在线观看| 91在线精品国自产拍蜜月| 晚上一个人看的免费电影| 精品人妻熟女毛片av久久网站| 91国产中文字幕| 成人漫画全彩无遮挡| 下体分泌物呈黄色| 亚洲图色成人| 99久久精品国产国产毛片| 国产深夜福利视频在线观看| 性高湖久久久久久久久免费观看| 中文天堂在线官网| 激情五月婷婷亚洲| 久久国内精品自在自线图片| 欧美 日韩 精品 国产| 国产精品二区激情视频| 老司机影院毛片| 91国产中文字幕| 一区二区三区精品91| 欧美中文综合在线视频| av女优亚洲男人天堂| 欧美精品高潮呻吟av久久| 伦理电影大哥的女人| a级毛片在线看网站| 亚洲国产av新网站| 久久热在线av| 久久久久精品性色| 99热全是精品| 国产精品麻豆人妻色哟哟久久| 欧美在线黄色| 色94色欧美一区二区| 亚洲一码二码三码区别大吗| 人人妻人人澡人人看| av天堂久久9| 精品少妇内射三级| 一区二区三区激情视频| 女人高潮潮喷娇喘18禁视频| 一个人免费看片子| 捣出白浆h1v1| 亚洲一码二码三码区别大吗| 爱豆传媒免费全集在线观看| 我的亚洲天堂| 女性被躁到高潮视频| 久久人人爽人人片av| 国产深夜福利视频在线观看| 亚洲第一av免费看| 在线天堂最新版资源| 少妇人妻 视频| 国产亚洲最大av| 欧美日韩精品网址| 亚洲精品一区蜜桃| 免费大片黄手机在线观看| 亚洲av国产av综合av卡| 亚洲精品国产一区二区精华液| 在现免费观看毛片| 丁香六月天网| 亚洲经典国产精华液单| 午夜日本视频在线| freevideosex欧美| 国产一区亚洲一区在线观看| 观看美女的网站| 亚洲精品第二区| 两个人免费观看高清视频| 美女脱内裤让男人舔精品视频| 亚洲国产av影院在线观看| 久久 成人 亚洲| 国产精品久久久久久精品古装| 免费看不卡的av| 国产女主播在线喷水免费视频网站| av免费在线看不卡| 一级毛片我不卡| 日韩中文字幕欧美一区二区 | 不卡视频在线观看欧美| 亚洲中文av在线| 亚洲av欧美aⅴ国产| av又黄又爽大尺度在线免费看| 少妇猛男粗大的猛烈进出视频| 精品少妇一区二区三区视频日本电影 | 亚洲国产毛片av蜜桃av| 女人久久www免费人成看片| 啦啦啦在线观看免费高清www| 9热在线视频观看99| 亚洲欧美成人精品一区二区| freevideosex欧美| 精品一区在线观看国产| 蜜桃国产av成人99| 在线观看免费视频网站a站| 亚洲美女黄色视频免费看| 国产成人av激情在线播放| 韩国av在线不卡| 国产精品麻豆人妻色哟哟久久| 日韩中文字幕视频在线看片| 青青草视频在线视频观看| 伊人久久大香线蕉亚洲五| 在线天堂最新版资源| 黄片无遮挡物在线观看| 欧美 亚洲 国产 日韩一| 久久久久国产网址| 免费高清在线观看视频在线观看| 久久精品国产综合久久久| videos熟女内射| 制服丝袜香蕉在线| 色视频在线一区二区三区| 午夜免费鲁丝| 免费看av在线观看网站| 在线天堂中文资源库| 哪个播放器可以免费观看大片| 亚洲国产精品一区三区| 亚洲伊人色综图| 久久精品国产亚洲av高清一级| 欧美中文综合在线视频| 三上悠亚av全集在线观看| 国产在线免费精品| 69精品国产乱码久久久| 国产精品嫩草影院av在线观看| 久久精品熟女亚洲av麻豆精品| 国产成人精品久久久久久| 久久久精品国产亚洲av高清涩受| 老司机影院成人| 1024香蕉在线观看| 亚洲av国产av综合av卡| 欧美国产精品一级二级三级| 亚洲国产精品成人久久小说| 久久狼人影院| 欧美日韩一区二区视频在线观看视频在线| 一本色道久久久久久精品综合| 丝袜脚勾引网站| 高清视频免费观看一区二区| 一区在线观看完整版| 一边亲一边摸免费视频| www日本在线高清视频| 国产av一区二区精品久久| 汤姆久久久久久久影院中文字幕| 免费观看性生交大片5| www.自偷自拍.com| 国产男女超爽视频在线观看| 久久久久久伊人网av| 91久久精品国产一区二区三区| 国产精品久久久久久精品古装| 国产成人精品久久二区二区91 | 韩国精品一区二区三区| 亚洲精品美女久久久久99蜜臀 | 国产免费一区二区三区四区乱码| 女人高潮潮喷娇喘18禁视频| 国产精品av久久久久免费| 日韩av免费高清视频| 秋霞在线观看毛片| 国产精品二区激情视频| 午夜日本视频在线| 多毛熟女@视频| 国产成人av激情在线播放| 欧美成人精品欧美一级黄| 男人舔女人的私密视频| 精品国产一区二区三区久久久樱花| 亚洲成国产人片在线观看| 一本久久精品| 中文字幕制服av| 久久久国产精品麻豆| 久久久久久伊人网av| 美国免费a级毛片| 女人高潮潮喷娇喘18禁视频| 熟女av电影| 汤姆久久久久久久影院中文字幕| 亚洲综合色惰| 久久精品国产综合久久久| 精品国产一区二区三区四区第35| 一级片免费观看大全| 国产精品欧美亚洲77777| 一边摸一边做爽爽视频免费| 18+在线观看网站| 视频在线观看一区二区三区| 国产人伦9x9x在线观看 | h视频一区二区三区| 电影成人av| 国产亚洲欧美精品永久| 中文乱码字字幕精品一区二区三区| 纯流量卡能插随身wifi吗| 在线天堂最新版资源| 日韩大片免费观看网站| 女的被弄到高潮叫床怎么办| 看十八女毛片水多多多| 亚洲天堂av无毛| 日韩欧美精品免费久久| 日本欧美视频一区| 国产黄色免费在线视频| 久久国产亚洲av麻豆专区| av在线老鸭窝| 最近最新中文字幕大全免费视频 | 亚洲人成电影观看| 卡戴珊不雅视频在线播放| 少妇精品久久久久久久| av网站免费在线观看视频| 亚洲情色 制服丝袜| 久久久久精品性色| a级片在线免费高清观看视频| 欧美日韩精品成人综合77777| 久久久久人妻精品一区果冻| 黄片播放在线免费| 999精品在线视频| 性色av一级| 欧美激情极品国产一区二区三区| 婷婷色综合大香蕉| av国产久精品久网站免费入址| 男女无遮挡免费网站观看| 汤姆久久久久久久影院中文字幕| 在线观看美女被高潮喷水网站| 国产精品久久久久久精品电影小说| 精品人妻熟女毛片av久久网站| 国产一区二区在线观看av| 久久国内精品自在自线图片| 欧美黄色片欧美黄色片| 国产亚洲欧美精品永久| 韩国av在线不卡| 狠狠婷婷综合久久久久久88av| 99九九在线精品视频| 亚洲第一青青草原| 色婷婷av一区二区三区视频| 狠狠精品人妻久久久久久综合| 国产女主播在线喷水免费视频网站| 高清欧美精品videossex| 日韩精品免费视频一区二区三区| 亚洲伊人久久精品综合| 九色亚洲精品在线播放| 久久人妻熟女aⅴ| 观看美女的网站| 亚洲在久久综合| 国产一区亚洲一区在线观看| 乱人伦中国视频| 大话2 男鬼变身卡| 成人18禁高潮啪啪吃奶动态图| 日韩精品免费视频一区二区三区| 一级片免费观看大全| 日韩中字成人| 九色亚洲精品在线播放| 日本黄色日本黄色录像| 天美传媒精品一区二区| 在线观看一区二区三区激情| 午夜日本视频在线| 韩国av在线不卡| 免费女性裸体啪啪无遮挡网站| 热99国产精品久久久久久7| 国产精品二区激情视频| 9色porny在线观看| 女人高潮潮喷娇喘18禁视频| 国产黄频视频在线观看| 久久97久久精品| 成人午夜精彩视频在线观看| 精品一区二区三区四区五区乱码 | 欧美日韩av久久| 午夜福利乱码中文字幕| 成人黄色视频免费在线看| 热99国产精品久久久久久7| 亚洲国产精品一区二区三区在线| 嫩草影院入口| 久久久久久久久久久免费av| 伦理电影免费视频| 丝瓜视频免费看黄片| 日韩一卡2卡3卡4卡2021年| 国产精品国产三级国产专区5o| 夫妻性生交免费视频一级片| 一本—道久久a久久精品蜜桃钙片| 欧美精品亚洲一区二区| 亚洲精品久久久久久婷婷小说| 久久久精品区二区三区| 欧美精品一区二区大全| 我要看黄色一级片免费的| 永久免费av网站大全| 婷婷色综合www| 免费少妇av软件| 免费在线观看黄色视频的| 狠狠精品人妻久久久久久综合| 汤姆久久久久久久影院中文字幕| 精品卡一卡二卡四卡免费| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 在线看a的网站| 亚洲情色 制服丝袜| 人妻 亚洲 视频| 国产极品粉嫩免费观看在线| 一级毛片我不卡| 久久久久国产网址| 18禁裸乳无遮挡动漫免费视频| 亚洲国产欧美日韩在线播放| 国产熟女午夜一区二区三区| a级毛片黄视频| 亚洲成色77777| av网站在线播放免费| 欧美成人精品欧美一级黄| 高清黄色对白视频在线免费看| 欧美精品国产亚洲| 伊人久久大香线蕉亚洲五| 超色免费av| 色播在线永久视频| 黄片小视频在线播放| 高清av免费在线| 性色av一级| 亚洲欧美一区二区三区黑人 | 韩国精品一区二区三区| 婷婷色麻豆天堂久久| 久久热在线av| 亚洲国产看品久久| 美女主播在线视频| 国产精品免费视频内射| 寂寞人妻少妇视频99o| 高清不卡的av网站| 国产亚洲av片在线观看秒播厂| 91aial.com中文字幕在线观看| 最近2019中文字幕mv第一页| 色网站视频免费| 天美传媒精品一区二区| 亚洲人成网站在线观看播放| 日韩制服骚丝袜av| 老汉色∧v一级毛片| 亚洲精品日韩在线中文字幕| 午夜日韩欧美国产| 欧美黄色片欧美黄色片| 日日爽夜夜爽网站| 黄片小视频在线播放| 中文精品一卡2卡3卡4更新| 亚洲av综合色区一区| 亚洲精品视频女| 免费日韩欧美在线观看| 女性生殖器流出的白浆| 色哟哟·www| 午夜福利视频在线观看免费| 欧美精品一区二区免费开放| 九色亚洲精品在线播放| 丝袜美腿诱惑在线| 国产免费一区二区三区四区乱码| 亚洲中文av在线| 777米奇影视久久| 丰满饥渴人妻一区二区三| 日韩人妻精品一区2区三区| 老司机影院成人| 亚洲成人av在线免费| 日本av手机在线免费观看| 丝袜喷水一区| 久久精品国产自在天天线| 日韩中文字幕欧美一区二区 | 亚洲精品,欧美精品| 咕卡用的链子| 久久久久久久久久久免费av| 男男h啪啪无遮挡| 中文字幕人妻丝袜制服|