• <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)用
    天堂√8在线中文| 国产97色在线日韩免费| 国产99白浆流出| 十分钟在线观看高清视频www| 免费看日本二区| 日本三级黄在线观看| xxx96com| 国产成人精品久久二区二区免费| 成人永久免费在线观看视频| 国产成+人综合+亚洲专区| 美女国产高潮福利片在线看| 亚洲成人免费电影在线观看| 久久欧美精品欧美久久欧美| 欧美又色又爽又黄视频| 精品一区二区三区av网在线观看| 国产精品一区二区免费欧美| 免费在线观看日本一区| 日韩欧美在线二视频| 天天躁狠狠躁夜夜躁狠狠躁| 黄网站色视频无遮挡免费观看| netflix在线观看网站| 免费在线观看日本一区| 国产亚洲精品久久久久5区| 欧美 亚洲 国产 日韩一| 日日夜夜操网爽| 法律面前人人平等表现在哪些方面| 热99re8久久精品国产| 麻豆av在线久日| 久久天堂一区二区三区四区| 在线永久观看黄色视频| 国产精品久久久人人做人人爽| 久久婷婷人人爽人人干人人爱| 欧美在线黄色| 天堂动漫精品| 欧美激情极品国产一区二区三区| www.999成人在线观看| 一区福利在线观看| 国产一级毛片七仙女欲春2 | 日韩成人在线观看一区二区三区| 日韩大尺度精品在线看网址| 亚洲成人免费电影在线观看| 国产精品国产高清国产av| 久久婷婷成人综合色麻豆| 亚洲五月色婷婷综合| 国产精品一区二区精品视频观看| 亚洲五月天丁香| 国产精品 欧美亚洲| 美国免费a级毛片| 在线av久久热| 久9热在线精品视频| 夜夜爽天天搞| 精品国产乱码久久久久久男人| 91大片在线观看| 久久久久国内视频| 国产精品,欧美在线| 啦啦啦免费观看视频1| 亚洲成人精品中文字幕电影| 成人av一区二区三区在线看| 国产99白浆流出| 国产野战对白在线观看| 国产亚洲精品第一综合不卡| 大型黄色视频在线免费观看| 999精品在线视频| 麻豆久久精品国产亚洲av| 色哟哟哟哟哟哟| 精品国产一区二区三区四区第35| 亚洲一码二码三码区别大吗| 亚洲成a人片在线一区二区| 亚洲在线自拍视频| 久久久久国产精品人妻aⅴ院| 日韩免费av在线播放| 欧美在线黄色| 国产成人精品无人区| 国产v大片淫在线免费观看| 窝窝影院91人妻| 丝袜在线中文字幕| 一二三四社区在线视频社区8| av视频在线观看入口| 大型av网站在线播放| 美女扒开内裤让男人捅视频| www日本在线高清视频| 国产99白浆流出| 中文字幕人妻丝袜一区二区| 欧美成人性av电影在线观看| 麻豆国产av国片精品| 一区二区三区精品91| 最近最新中文字幕大全免费视频| 制服人妻中文乱码| 日韩一卡2卡3卡4卡2021年| 日韩中文字幕欧美一区二区| 首页视频小说图片口味搜索| 久久久久国内视频| 这个男人来自地球电影免费观看| av在线天堂中文字幕| 99久久综合精品五月天人人| 精品欧美国产一区二区三| 搡老妇女老女人老熟妇| 俄罗斯特黄特色一大片| 国产精品日韩av在线免费观看| 久久精品国产亚洲av高清一级| 狠狠狠狠99中文字幕| 国产精品一区二区三区四区久久 | 黄色视频不卡| 久久久久国产一级毛片高清牌| 97碰自拍视频| 欧美 亚洲 国产 日韩一| 正在播放国产对白刺激| 欧美久久黑人一区二区| 在线国产一区二区在线| 午夜日韩欧美国产| 国产免费男女视频| 看黄色毛片网站| 久久伊人香网站| 亚洲午夜理论影院| 亚洲国产日韩欧美精品在线观看 | 老汉色av国产亚洲站长工具| 国产一区二区激情短视频| 国产精品久久久久久人妻精品电影| 亚洲av熟女| 欧美成人一区二区免费高清观看 | 欧美乱妇无乱码| 欧美黑人巨大hd| 一级片免费观看大全| 1024香蕉在线观看| 国产一级毛片七仙女欲春2 | 久久久久久久午夜电影| 亚洲性夜色夜夜综合| 亚洲专区中文字幕在线| 妹子高潮喷水视频| 18禁观看日本| 欧美三级亚洲精品| 精品国产超薄肉色丝袜足j| 麻豆久久精品国产亚洲av| 亚洲午夜理论影院| 天天躁夜夜躁狠狠躁躁| 欧美+亚洲+日韩+国产| 欧美黑人精品巨大| 欧美黄色片欧美黄色片| 国产一区二区激情短视频| 久久午夜亚洲精品久久| 桃红色精品国产亚洲av| 国产区一区二久久| 午夜福利成人在线免费观看| 91国产中文字幕| 国产单亲对白刺激| 可以在线观看的亚洲视频| 欧美日韩福利视频一区二区| 波多野结衣av一区二区av| 亚洲真实伦在线观看| 国产一区二区在线av高清观看| 又大又爽又粗| 91麻豆av在线| 久久久久久久久免费视频了| 久久精品国产综合久久久| 韩国av一区二区三区四区| 精品国产国语对白av| 777久久人妻少妇嫩草av网站| 国产精品国产高清国产av| 女生性感内裤真人,穿戴方法视频| 一级毛片高清免费大全| 欧美一级a爱片免费观看看 | 国产精品,欧美在线| 变态另类丝袜制服| 亚洲欧美一区二区三区黑人| a级毛片在线看网站| 又大又爽又粗| 在线视频色国产色| 宅男免费午夜| 久久久久久久久久黄片| 国产欧美日韩一区二区三| 国产熟女xx| 真人做人爱边吃奶动态| 午夜免费激情av| 国产真人三级小视频在线观看| 亚洲va日本ⅴa欧美va伊人久久| 欧洲精品卡2卡3卡4卡5卡区| 97超级碰碰碰精品色视频在线观看| 日本撒尿小便嘘嘘汇集6| 老司机午夜福利在线观看视频| 国产av一区在线观看免费| 久久午夜综合久久蜜桃| 精品国内亚洲2022精品成人| 亚洲av片天天在线观看| 高清在线国产一区| 国产一区二区在线av高清观看| 久久精品影院6| 亚洲片人在线观看| 国产精品久久久av美女十八| 久久精品人妻少妇| 久久久久亚洲av毛片大全| 18禁黄网站禁片免费观看直播| 美女免费视频网站| 99在线人妻在线中文字幕| 国产精品一区二区三区四区久久 | 韩国精品一区二区三区| 欧美又色又爽又黄视频| 变态另类丝袜制服| 亚洲精品中文字幕一二三四区| 在线观看午夜福利视频| 1024手机看黄色片| 免费在线观看亚洲国产| 男男h啪啪无遮挡| 亚洲欧美精品综合一区二区三区| 999久久久国产精品视频| 美女午夜性视频免费| av中文乱码字幕在线| a级毛片a级免费在线| 又大又爽又粗| 国产97色在线日韩免费| 国产又黄又爽又无遮挡在线| 久久精品亚洲精品国产色婷小说| 亚洲久久久国产精品| 夜夜看夜夜爽夜夜摸| 精品欧美一区二区三区在线| 女警被强在线播放| 国产v大片淫在线免费观看| 国产精品 欧美亚洲| 黑人操中国人逼视频| 色老头精品视频在线观看| 国产精品综合久久久久久久免费| 精华霜和精华液先用哪个| 啦啦啦 在线观看视频| 久久久久久久精品吃奶| 一边摸一边做爽爽视频免费| 久热爱精品视频在线9| 欧美国产精品va在线观看不卡| 身体一侧抽搐| 在线永久观看黄色视频| 俄罗斯特黄特色一大片| 丝袜人妻中文字幕| 18美女黄网站色大片免费观看| 欧美丝袜亚洲另类 | 超碰成人久久| 欧美亚洲日本最大视频资源| 色综合站精品国产| 日韩大码丰满熟妇| 十分钟在线观看高清视频www| 一级a爱视频在线免费观看| 一进一出好大好爽视频| 亚洲av日韩精品久久久久久密| 一级片免费观看大全| 久99久视频精品免费| www国产在线视频色| 美女扒开内裤让男人捅视频| 美女午夜性视频免费| 国产不卡一卡二| 亚洲精品av麻豆狂野| 国产1区2区3区精品| 亚洲国产日韩欧美精品在线观看 | 久久精品成人免费网站| 夜夜夜夜夜久久久久| 国产久久久一区二区三区| 精品久久久久久久人妻蜜臀av| 亚洲最大成人中文| 亚洲电影在线观看av| 成人三级黄色视频| 欧美激情 高清一区二区三区| www.999成人在线观看| 久久久精品欧美日韩精品| 久久久久久国产a免费观看| 嫁个100分男人电影在线观看| www.精华液| 男男h啪啪无遮挡| 99国产综合亚洲精品| 日本一区二区免费在线视频| 人人澡人人妻人| 国产精品久久久久久亚洲av鲁大| 露出奶头的视频| 午夜老司机福利片| 国产在线精品亚洲第一网站| 熟女电影av网| 丁香欧美五月| 禁无遮挡网站| 高清在线国产一区| 免费无遮挡裸体视频| 白带黄色成豆腐渣| 国产精品影院久久| 亚洲精品粉嫩美女一区| 国产精品久久视频播放| 人妻丰满熟妇av一区二区三区| 又紧又爽又黄一区二区| 日韩免费av在线播放| 国产私拍福利视频在线观看| 免费搜索国产男女视频| 18禁黄网站禁片午夜丰满| 少妇熟女aⅴ在线视频| 天天躁夜夜躁狠狠躁躁| 久久中文看片网| 啪啪无遮挡十八禁网站| 99热6这里只有精品| 久久国产精品人妻蜜桃| 成年版毛片免费区| 最好的美女福利视频网| 男人舔女人下体高潮全视频| 国产午夜精品久久久久久| 国产av一区二区精品久久| 午夜福利成人在线免费观看| 国产成人精品久久二区二区免费| 成人特级黄色片久久久久久久| 国产不卡一卡二| 一夜夜www| aaaaa片日本免费| 我的亚洲天堂| 白带黄色成豆腐渣| 18禁黄网站禁片免费观看直播| 成在线人永久免费视频| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美久久黑人一区二区| 在线观看免费视频日本深夜| 国产精品久久久av美女十八| 亚洲av电影在线进入| 高清毛片免费观看视频网站| 十八禁人妻一区二区| 免费在线观看黄色视频的| 精品久久久久久久久久免费视频| 又紧又爽又黄一区二区| 久久精品国产亚洲av香蕉五月| 99国产精品一区二区蜜桃av| 免费看日本二区| 老司机深夜福利视频在线观看| 韩国精品一区二区三区| 一个人免费在线观看的高清视频| 欧美成人一区二区免费高清观看 | 黄色视频,在线免费观看| 日韩欧美免费精品| 久久精品国产99精品国产亚洲性色| 叶爱在线成人免费视频播放| 国产麻豆成人av免费视频| 制服人妻中文乱码| 国产视频一区二区在线看| 欧美黑人精品巨大| 国产精品自产拍在线观看55亚洲| 嫁个100分男人电影在线观看| 一本精品99久久精品77| 老熟妇乱子伦视频在线观看| 国产精品综合久久久久久久免费| 深夜精品福利| 久久午夜亚洲精品久久| 精品一区二区三区av网在线观看| 婷婷亚洲欧美| 久久欧美精品欧美久久欧美| 免费高清在线观看日韩| 免费看十八禁软件| 中出人妻视频一区二区| 老司机深夜福利视频在线观看| 男人舔女人的私密视频| 国产黄a三级三级三级人| 亚洲人成电影免费在线| 国产精品电影一区二区三区| 亚洲五月色婷婷综合| 婷婷六月久久综合丁香| 日韩av在线大香蕉| 狂野欧美激情性xxxx| 日韩国内少妇激情av| 国产亚洲精品第一综合不卡| 叶爱在线成人免费视频播放| 夜夜夜夜夜久久久久| 9191精品国产免费久久| 久久久久亚洲av毛片大全| 美女高潮到喷水免费观看| 欧美乱妇无乱码| 美女扒开内裤让男人捅视频| 国产爱豆传媒在线观看 | av免费在线观看网站| 亚洲片人在线观看| 午夜福利在线观看吧| 成人18禁在线播放| 一边摸一边做爽爽视频免费| 欧美午夜高清在线| 狠狠狠狠99中文字幕| 亚洲第一欧美日韩一区二区三区| 久久久久久久久中文| 久久99热这里只有精品18| 国产亚洲欧美在线一区二区| 99热6这里只有精品| 无限看片的www在线观看| 人人妻人人澡欧美一区二区| 亚洲五月婷婷丁香| 国内久久婷婷六月综合欲色啪| 亚洲专区国产一区二区| 午夜免费成人在线视频| avwww免费| 丰满的人妻完整版| 久久九九热精品免费| 啦啦啦观看免费观看视频高清| 在线国产一区二区在线| 欧美又色又爽又黄视频| 欧美在线一区亚洲| 在线播放国产精品三级| 久热爱精品视频在线9| 一级黄色大片毛片| 免费看a级黄色片| 亚洲一码二码三码区别大吗| 国产区一区二久久| 熟女少妇亚洲综合色aaa.| 中文字幕久久专区| 国产午夜福利久久久久久| 欧美av亚洲av综合av国产av| 丰满人妻熟妇乱又伦精品不卡| 久久久久久国产a免费观看| 18禁裸乳无遮挡免费网站照片 | 欧美激情久久久久久爽电影| 久久国产精品影院| 神马国产精品三级电影在线观看 | 日韩欧美国产在线观看| 亚洲色图 男人天堂 中文字幕| 级片在线观看| 国产伦一二天堂av在线观看| 国产精品综合久久久久久久免费| 日韩欧美一区二区三区在线观看| 亚洲无线在线观看| 久久久久久久精品吃奶| 热99re8久久精品国产| 淫妇啪啪啪对白视频| 国产亚洲精品久久久久久毛片| 美女免费视频网站| 日韩 欧美 亚洲 中文字幕| 变态另类丝袜制服| 久久天躁狠狠躁夜夜2o2o| 午夜老司机福利片| 国产精品电影一区二区三区| 老司机在亚洲福利影院| 日韩欧美国产一区二区入口| 国产精品久久久久久精品电影 | 国产成人影院久久av| 亚洲欧美精品综合久久99| 男女那种视频在线观看| 婷婷亚洲欧美| 国产精品综合久久久久久久免费| 99国产综合亚洲精品| 99久久综合精品五月天人人| 国产精品亚洲一级av第二区| av有码第一页| 亚洲男人天堂网一区| 国产欧美日韩一区二区精品| 女人爽到高潮嗷嗷叫在线视频| 久久精品夜夜夜夜夜久久蜜豆 | 国产黄a三级三级三级人| 18禁黄网站禁片午夜丰满| 午夜日韩欧美国产| 女性生殖器流出的白浆| 18禁黄网站禁片午夜丰满| 亚洲av片天天在线观看| 国产一区二区在线av高清观看| 丰满人妻熟妇乱又伦精品不卡| 一区福利在线观看| 久久精品国产99精品国产亚洲性色| 国产av在哪里看| 1024视频免费在线观看| 午夜免费鲁丝| 黄片小视频在线播放| 久久伊人香网站| 日韩精品免费视频一区二区三区| 中出人妻视频一区二区| 久久久国产欧美日韩av| 黄色片一级片一级黄色片| 不卡一级毛片| 午夜福利视频1000在线观看| 一二三四在线观看免费中文在| 一级片免费观看大全| 亚洲最大成人中文| 男女之事视频高清在线观看| 女人爽到高潮嗷嗷叫在线视频| or卡值多少钱| 久久久久久久久久黄片| 亚洲黑人精品在线| 亚洲成a人片在线一区二区| 午夜影院日韩av| 少妇的丰满在线观看| 欧美av亚洲av综合av国产av| 人妻久久中文字幕网| 欧美又色又爽又黄视频| 午夜激情福利司机影院| 好男人电影高清在线观看| 国内精品久久久久久久电影| 婷婷精品国产亚洲av在线| 老司机午夜十八禁免费视频| 午夜福利在线在线| 亚洲一区高清亚洲精品| 美女大奶头视频| 亚洲第一欧美日韩一区二区三区| 自线自在国产av| 国产成人精品久久二区二区91| 欧美大码av| 欧美另类亚洲清纯唯美| 成人国产综合亚洲| 免费av毛片视频| 国产色视频综合| 看黄色毛片网站| 亚洲人成网站在线播放欧美日韩| a级毛片a级免费在线| 欧美成人免费av一区二区三区| 好看av亚洲va欧美ⅴa在| 天堂影院成人在线观看| 免费看美女性在线毛片视频| 丝袜人妻中文字幕| 国产精品,欧美在线| 欧美在线黄色| 久久精品成人免费网站| 国产亚洲精品综合一区在线观看 | 88av欧美| 女生性感内裤真人,穿戴方法视频| 亚洲av熟女| 日日摸夜夜添夜夜添小说| 久久国产乱子伦精品免费另类| 青草久久国产| 99在线人妻在线中文字幕| 欧美日韩一级在线毛片| 真人做人爱边吃奶动态| 亚洲国产精品sss在线观看| 一个人观看的视频www高清免费观看 | 久久久久九九精品影院| 国产片内射在线| 久久国产精品影院| 亚洲精品美女久久久久99蜜臀| 久久九九热精品免费| 午夜免费鲁丝| 国产午夜福利久久久久久| 中文字幕人成人乱码亚洲影| av超薄肉色丝袜交足视频| 国产在线精品亚洲第一网站| 每晚都被弄得嗷嗷叫到高潮| av有码第一页| 成人国语在线视频| 亚洲国产高清在线一区二区三 | 国内少妇人妻偷人精品xxx网站 | 午夜福利欧美成人| 一本久久中文字幕| 国产精品野战在线观看| 成人三级黄色视频| a级毛片在线看网站| 亚洲欧洲精品一区二区精品久久久| 少妇粗大呻吟视频| 男女之事视频高清在线观看| 一级毛片女人18水好多| 波多野结衣av一区二区av| АⅤ资源中文在线天堂| 午夜亚洲福利在线播放| 国产成人欧美在线观看| 欧美日韩精品网址| 视频区欧美日本亚洲| xxxwww97欧美| 精品国产超薄肉色丝袜足j| 日韩av在线大香蕉| 国产人伦9x9x在线观看| 亚洲成人久久性| 日日干狠狠操夜夜爽| 一区二区三区国产精品乱码| 巨乳人妻的诱惑在线观看| 精品久久久久久久末码| 9191精品国产免费久久| 美女免费视频网站| 一区二区三区精品91| 亚洲国产精品合色在线| 午夜福利一区二区在线看| 亚洲av第一区精品v没综合| 91老司机精品| 欧美+亚洲+日韩+国产| 51午夜福利影视在线观看| 国产激情欧美一区二区| 久久精品国产清高在天天线| 老司机午夜福利在线观看视频| 每晚都被弄得嗷嗷叫到高潮| 国产成人精品无人区| 国产精品免费视频内射| 日本免费a在线| 日日爽夜夜爽网站| 黄片小视频在线播放| 亚洲精品国产一区二区精华液| 久久久久亚洲av毛片大全| 亚洲五月婷婷丁香| 曰老女人黄片| 黄色视频,在线免费观看| 最新美女视频免费是黄的| 国产欧美日韩一区二区三| 日本撒尿小便嘘嘘汇集6| 亚洲av成人一区二区三| 亚洲久久久国产精品| 久久草成人影院| 日本五十路高清| 后天国语完整版免费观看| 亚洲五月色婷婷综合| 国产成人精品久久二区二区91| 国产在线精品亚洲第一网站| 国产又色又爽无遮挡免费看| 好男人在线观看高清免费视频 | 天堂√8在线中文| 欧美绝顶高潮抽搐喷水| 老鸭窝网址在线观看| 国产99白浆流出| 国产熟女午夜一区二区三区| 午夜a级毛片| 怎么达到女性高潮| 婷婷六月久久综合丁香| 一本综合久久免费| a级毛片在线看网站| 亚洲成人免费电影在线观看| 成年女人毛片免费观看观看9| 一进一出好大好爽视频| 亚洲 欧美一区二区三区| 亚洲五月天丁香| 老汉色∧v一级毛片| 黄色片一级片一级黄色片| 久久中文看片网| 久久 成人 亚洲| 欧美日韩乱码在线| 日韩欧美 国产精品| 国产精品亚洲av一区麻豆| 久久99热这里只有精品18| 亚洲专区字幕在线| 亚洲精品久久国产高清桃花| 亚洲国产欧美日韩在线播放| 女警被强在线播放| 久久人人精品亚洲av| 夜夜夜夜夜久久久久|