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

    基于最大熵卡爾曼濾波算法的液壓支架調(diào)直方法

    2022-12-07 17:24:02宋單陽盧春貴陶心雅楊金衡王培恩鄭文強宋建成
    工礦自動化 2022年11期
    關(guān)鍵詞:直線度高斯軌跡

    宋單陽,盧春貴,陶心雅,楊金衡,王培恩,鄭文強,宋建成

    (1. 太原理工大學(xué) 礦用智能電器技術(shù)國家地方聯(lián)合工程實驗室,山西 太原 030024;2. 太原理工大學(xué) 煤礦電氣設(shè)備與智能控制山西省重點實驗室,山西 太原 030024;3. 太原科技大學(xué) 電子信息工程學(xué)院,山西 太原 030024)

    0 引言

    液壓支架作為綜采工作面的重要生產(chǎn)設(shè)備之一,不僅起著支撐頂板的作用,還與采煤機、刮板輸送機等綜采設(shè)備互相配合、協(xié)同作業(yè),將工作面不斷向前推進。工作面煤壁、刮板輸送機和液壓支架都保持直線狀態(tài)是綜采工作面綠色高效自動化生產(chǎn)的技術(shù)保障。然而在實際生產(chǎn)過程中,由于測量誤差、推移誤差、累計誤差等原因,液壓支架排列并非在一條直線上,這樣會直接影響刮板輸送機直線度,從而影響綜采工作面生產(chǎn)效率、增大能源消耗、加速設(shè)備損壞,另外還存在極大的安全隱患[1-2]。因此,如何實現(xiàn)液壓支架直線度精確檢測、快速調(diào)直,是提高綜采工作面三機設(shè)備運行效率,加速綜采工作面自動化建設(shè)亟需解決的關(guān)鍵問題[3]。

    近年來,煤礦領(lǐng)域許多學(xué)者對液壓支架的調(diào)直方法進行了大量研究。文獻(xiàn)[4]通過在架間行走機器人上布置傳感器來測量液壓支架底座多維位置偏移數(shù)據(jù),從而獲得液壓支架直線度信息,但該方法需要改裝液壓支架的人行板和布置輔助測量板,不適宜在井下實際使用。文獻(xiàn)[5]基于支持向量機和遺傳算法構(gòu)建了液壓支架調(diào)直系統(tǒng)模型,但該模型容易受到頂板壓力和乳化液泵站壓力的影響,對壓力傳感器的精度要求比較高。文獻(xiàn)[6-9]主要是通過在液壓支架之間安裝大量傳感器來判斷液壓支架的相對位置,通過電液控制系統(tǒng)控制液壓支架進行“拉架”和“推溜”,以實現(xiàn)液壓支架之間的位置調(diào)整,從而控制整個工作面液壓支架群的直線度;這些方法的缺點是傳感器數(shù)量多,極易產(chǎn)生累計誤差,且在液壓支架調(diào)直過程中沒有考慮推移系統(tǒng)誤差的影響,無法達(dá)到理想的調(diào)直效果。文獻(xiàn)[10]通過對液壓支架超調(diào)量進行不完全補償來實現(xiàn)液壓支架的調(diào)直,但當(dāng)液壓支架推移系統(tǒng)的系統(tǒng)噪聲為非高斯噪聲時,該方法是否具有可行性有待驗證。文獻(xiàn)[11]提出了一種基于卡爾曼濾波(Kalman Filter,KF)的刮板輸送機調(diào)直方法,通過實時監(jiān)測刮板輸送機中部槽的位置完成對刮板輸送機的調(diào)直,最后通過電液控制系統(tǒng)控制液壓支架進行“拉架”動作,實現(xiàn)液壓支架直線度的調(diào)整,但該方法易受系統(tǒng)非高斯量測噪聲的影響。

    綜上可知,現(xiàn)有液壓支架調(diào)直方法受到傳感器測量誤差和液壓支架推移誤差的影響,使得調(diào)直誤差較大;在非高斯量測噪聲環(huán)境下,傳統(tǒng)基于KF 算法的調(diào)直方法對液壓支架軌跡的預(yù)測準(zhǔn)確度低,無法達(dá)到理想的調(diào)直效果。針對上述問題,本文在KF迭代過程中引入最大熵準(zhǔn)則(Maximum Correntropy Criterion,MCC)處理非高斯量測噪聲,提出了一種基于最大熵卡爾曼濾波(MCC Kalman Filter,MCKF)算法的液壓支架調(diào)直方法,以提高調(diào)直精度。

    1 液壓支架調(diào)直方法

    基于MCKF 算法的液壓支架調(diào)直方法的基本原理:首先根據(jù)液壓支架的位置坐標(biāo)和工作面推進方向確定調(diào)直參考直線;然后構(gòu)建液壓支架線性推移系統(tǒng)的狀態(tài)方程和觀測方程,經(jīng)MCKF 算法處理后得到液壓支架推移后的預(yù)測軌跡;最后根據(jù)液壓支架預(yù)測軌跡與調(diào)直參考直線解算出每架液壓支架的推移距離補償量,從而達(dá)到調(diào)直目的。

    液壓支架調(diào)直過程如圖1 所示。采煤機完成第k-1 次截割后,液壓支架的實際軌跡為TD(k-1),取實際軌跡TD(k-1)上最滯后點向工作面推進方向作垂線得到參考直線LC(k-1) ,將參考直線LC(k-1)向工作面推進方向平移距離h即可得到采煤機第k次截割時的液壓支架理想軌跡TG(k)。取理想軌跡TG(k)與實際軌跡TD(k-1)上各移架點的差值作為液壓支架第k次推移距離,在檢測誤差和推移誤差的影響下推移液壓支架,形成下一次液壓支架的軌跡TD(k)。反復(fù)循環(huán)以上過程,不斷對液壓支架的直線度進行修正,使其最終穩(wěn)定在一定范圍內(nèi)。

    圖 1 液壓支架調(diào)直過程Fig. 1 Hydraulic support straightening process

    2 基于MCKF 算法的液壓支架軌跡預(yù)測

    綜采工作面液壓支架推移系統(tǒng)屬于線性系統(tǒng),其狀態(tài)方程和觀測方程可以反映液壓支架的推移軌跡狀態(tài),本文利用MCKF 算法對液壓支架推移系統(tǒng)的狀態(tài)變量進行更新,以實現(xiàn)液壓支架軌跡預(yù)測,通過引入MCC 獲得濾波誤差的高階統(tǒng)計信息,使系統(tǒng)跟蹤性能得到極大改善,提高噪聲濾除效果和調(diào)直精度。

    2.1 MCC

    對于給定的2 個隨機變量U和V,其相關(guān)熵定義如下:

    式中:φ(u,v)為核函數(shù),u、v為隨機變量U、V的具體值;ρ(u,v)為 隨機變量U和V的聯(lián)合概率密度函數(shù)。

    φ(u,v)核函數(shù)具有線性、多項式、高斯和指數(shù)核函數(shù)等多種形式[12-13],本文選擇常用的高斯核函數(shù)。

    式中:Gσ()為帶寬為σ的核函數(shù);t=u-v。

    MCC 不受高斯分布假設(shè)的影響[14-15],本文在KF 迭代過程中引入MCC,在非高斯噪聲環(huán)境下,使得液壓支架推移系統(tǒng)的跟蹤性能得到極大改善。

    2.2 MCKF 算法

    液壓支架推移系統(tǒng)的狀態(tài)方程和觀測方程可以描述為

    式中:A為轉(zhuǎn)移矩陣;Y(k)為采煤機第k-1 次割煤后液壓支架在推進方向上的坐標(biāo);B為控制矩陣;D(k)為采煤機第k次割煤后液壓支架的推進距離;w(k)為過程噪聲;Z(k)為采煤機第k-1 次割煤后液壓支架在推進方向的測量坐標(biāo);H為測量矩陣;r(k)為量測噪聲。

    KF 算法的遞推表達(dá)式是通過最小化系統(tǒng)代價函數(shù)實現(xiàn)的。KF 算法的代價函數(shù)為式中:Y(k|k-1) 為狀態(tài)變量預(yù)測值;P(k|k-1)為預(yù)測誤差協(xié)方差矩陣;R(k)為量測噪聲的協(xié)方差矩陣。

    為了加強KF 算法在非高斯量測噪聲下的跟蹤效果,將相關(guān)熵的概念引入代價函數(shù)的量測部分。定義ek=(R(k))-1/2(HY(k|k-1)-Z(k)),根據(jù)式(4),可得MCKF 算法的代價函數(shù):

    式中ek,i為ek的第i(i=1,2,···,m,m為分量個數(shù))個分量。

    對式(5)在Y(k)處求導(dǎo),并令其為零,可得

    式中ψ(k) 為權(quán)矩陣,ψ(k)=diag(Gσ(ek,i))。

    式(6)可以視為JMCKF在Y(k)處的導(dǎo)數(shù):

    對比式(4)和式(7)可看出,MCKF 算法在KF 算法的基礎(chǔ)上對量測噪聲的協(xié)方差進行了修改,從而改善系統(tǒng)的跟蹤性能。

    基于線性系統(tǒng)方程,MCKF 算法運算步驟如下:

    (1) 初始化:根據(jù)MCKF 算法迭代要求,選擇合適的核函數(shù)帶寬 σ,設(shè)定狀態(tài)變量初始值為Y(0)和對應(yīng)的誤差協(xié)方差矩陣為P(0)。

    (2) 狀態(tài)變量預(yù)測計算:

    (3) 預(yù)測誤差協(xié)方差矩陣計算:

    式中Q為過程噪聲的方差。

    (4) 修改量測噪聲協(xié)方差矩陣:MCKF 算法首先利用MCC 處理非高斯量測噪聲,根據(jù)狀態(tài)變量預(yù)測值的信息計算ek,然后基于ek計算權(quán)矩陣 ψ(k),最后用權(quán)矩陣修改量測噪聲協(xié)方差矩陣:

    (5) 濾波增益矩陣計算:

    (6) 狀態(tài)更新:

    (7) 誤差協(xié)方差矩陣更新:

    式中I為單位矩陣。

    與KF 算法相比,MCKF 算法基于MCC 修改量測噪聲協(xié)方差矩陣,使其在非高斯噪聲條件下改善濾波效果,從而解決了非高斯噪聲下液壓支架調(diào)直精度低的問題。

    3 仿真及結(jié)果分析

    為了檢驗基于MCKF 算法的液壓支架調(diào)直方法的有效性,基于Matlab 對液壓支架調(diào)直過程進行了數(shù)值仿真。為了簡化計算,假設(shè)在東北天坐標(biāo)系下,液壓支架沿東方向排列,推進方向向北。以100 架液壓支架作為仿真對象,選擇液壓支架中心距為1.5 m,標(biāo)準(zhǔn)推移距離為0.8 m。由于井下液壓支架過程噪聲和量測噪聲的方差變化范圍較大,為0~1,本文選用該變化范圍內(nèi)的噪聲方差數(shù)據(jù),在量測噪聲分別為高斯噪聲和非高斯噪聲環(huán)境下進行仿真實驗,并與基于KF 算法的調(diào)直方法進行對比。

    3.1 高斯噪聲環(huán)境下的仿真

    假設(shè)液壓支架過程噪聲 δa和量測噪聲 δb為隨機噪聲,相互獨立并服從高斯分布(δa~N(0,0.006),δb~N(0,0.006))。液壓支架的初始軌跡如圖2 所示,將液壓支架推移100 次,得到液壓支架的推移軌跡,如圖3 所示。選擇預(yù)測軌跡數(shù)據(jù)與真實軌跡數(shù)據(jù)之間的均方誤差(Mean Square Error,MSE)作為衡量液壓支架軌跡預(yù)測效果的評判依據(jù),MSE 越小,說明預(yù)測效果越好。由圖3 可看出,在高斯噪聲條件下,KF 算法的預(yù)測軌跡與真實軌跡之間的MSE 為1.78 mm,MCKF 算法的預(yù)測軌跡與真實軌跡之間的MSE 為1.26 mm,可看出基于KF 算法和MCKF 算法的調(diào)直方法對于液壓支架直線度檢測均有很好的效果。

    圖 2 高斯噪聲條件下液壓支架初始軌跡Fig. 2 Initial trajectory of hydraulic support under condition of Gaussian noise

    圖 3 高斯噪聲條件下液壓支架推移軌跡Fig. 3 Moving trajectory of hydraulic support under condition of Gaussian noise

    在液壓支架推移100 次過程中,可以得到KF 算法預(yù)測軌跡與真實軌跡之間的MSE 及MCKF 算法預(yù)測軌跡與真實軌跡之間的MSE 隨推移次數(shù)變化的曲線,如圖4 所示。KF 算法的MSE 平均值為2.36 mm,MCKF 算法的MSE 平均值為1.70 mm。2 種算法的預(yù)測精度均較高,但MCKF 算法的預(yù)測結(jié)果更加逼近實際值。

    圖 4 高斯噪聲條件下KF 算法與MCKF 算法預(yù)測軌跡的MSEFig. 4 Mean square errors of prediction trajectory of KF algorithm and MCKF algorithm under condition of Gaussian noise

    選取液壓支架在北方向坐標(biāo)最大值與最小值之差作為液壓支架最大直線度誤差。在高斯噪聲條件下,利用KF 算法和MCKF 算法對液壓支架直線度進行軌跡預(yù)測后,推移液壓支架可以得到其最大直線度誤差隨推移次數(shù)變化的曲線,如圖5 所示??煽闯鯧F 算法的預(yù)測軌跡直線度誤差和MCKF 算法的預(yù)測軌跡直線度誤差均為38 mm 左右,相比于接近100 mm 的液壓支架初始軌跡直線度誤差,2 種算法的調(diào)直效果均很好,說明MCKF 算法在高斯噪聲環(huán)境下的預(yù)測效果也較好。

    圖 5 高斯噪聲條件下液壓支架直線度誤差隨推移次數(shù)變化曲線Fig. 5 Variation curves of the straightness error of hydraulic support with moving number under condition of Gaussian noise

    3.2 非高斯噪聲環(huán)境下的仿真

    假設(shè)液壓支架過程噪聲δa服從高斯分布,量測噪聲δb為 非高斯噪聲,即δb~0.8N(0,0.006)+0.2N(0,0.03)。液壓支架的初始軌跡如圖6 所示,將液壓支架推移100 次,得到液壓支架的推移軌跡如圖7 所示。在非高斯噪聲條件下,KF 算法的預(yù)測軌跡與真實軌跡之間的MSE 為16.82 mm,MCKF 算法的預(yù)測軌跡與真實軌跡之間的MSE 為4.35 mm,可以看出,相比KF 算法,MCKF 算法對液壓支架軌跡的估計一致性相對較好。

    圖 6 非高斯噪聲條件下液壓支架初始軌跡Fig. 6 Initial trajectory of hydraulic support under condition of non-Gaussian noise

    圖 7 非高斯噪聲條件下液壓支架推移軌跡Fig. 7 Moving trajectory of hydraulic support undercondition of non-Gaussian noise

    液壓支架推移100 次,可以得到KF 算法預(yù)測軌跡與真實軌跡之間的MSE 及MCKF 算法預(yù)測軌跡與真實軌跡之間的MSE 隨推移次數(shù)變化的曲線,如圖8 所示。KF 算法的MSE 平均值為16.74 mm,MCKF算法的MSE 平均值為4.76 mm。MCKF 算法的MSE遠(yuǎn)遠(yuǎn)小于KF 算法的MSE,能更加真實地反映液壓支架的真實狀態(tài)。

    圖 8 非高斯噪聲條件下KF 算法預(yù)測軌跡與MCKF 算法預(yù)測軌跡的MSEFig. 8 Mean square errors of prediction trajectory of KF algorithm and MCKF algorithm under condition of non-Gaussian noise

    在非高斯噪聲條件下,利用KF 算法和MCKF 算法對液壓支架直線度進行軌跡預(yù)測后,推移液壓支架可以得到其最大直線度誤差隨推移次數(shù)變化的曲線,如圖9 所示??煽闯鯧F 算法的預(yù)測軌跡直線度誤差為72 mm 左右,MCKF 算法的預(yù)測軌跡直線度誤差為46 mm 左右。相比于接近100 mm的液壓支架初始軌跡直線度誤差,經(jīng)MCKF 算法濾波后的軌跡比經(jīng)KF 算法濾波后的軌跡直線度誤差減小了36%,說明基于MCKF 算法的調(diào)直方法的效果比基于KF 算法的調(diào)直方法更理想。此外,還可看出,液壓支架調(diào)直誤差不受推移次數(shù)的影響,即直線度誤差只與本次調(diào)直過程有關(guān),與之前的調(diào)直過程無關(guān),有效避免了累計誤差。

    圖 9 非高斯噪聲條件下液壓支架直線度誤差隨推移次數(shù)變化曲線Fig. 9 Variation curves of the straightness error of hydraulic support with moving number under condition of non-Gaussian noise

    4 結(jié)論

    (1) 提出了一種基于MCKF 算法的液壓支架調(diào)直方法,通過MCKF 算法精確預(yù)測液壓支架推移后的軌跡,以計算液壓支架下一次的推移距離,根據(jù)兩者間的差值進行調(diào)直,解決了非高斯量測噪聲干擾下傳統(tǒng)KF 算法對液壓支架直線度檢測精度低,導(dǎo)致調(diào)直效果差的問題。

    (2) 仿真結(jié)果表明:與基于KF 算法的調(diào)直方法相比,基于MCKF 算法的液壓支架調(diào)直方法能夠有效降低量測噪聲和過程噪聲對液壓支架直線度的影響,特別當(dāng)量測噪聲服從非高斯分布時,該方法的MSE 平均值僅為4.76 mm,遠(yuǎn)小于基于KF 算法的調(diào)直方法的MSE,能夠更加準(zhǔn)確地預(yù)測液壓支架直線度。MCKF 算法的預(yù)測軌跡直線度誤差比KF 算法的預(yù)測軌跡直線度誤差減小了36%,調(diào)直效果更佳,且液壓支架的直線度誤差只與本次調(diào)直過程有關(guān),而與之前的調(diào)直過程無關(guān),有效避免了累計誤差。

    猜你喜歡
    直線度高斯軌跡
    小高斯的大發(fā)現(xiàn)
    長導(dǎo)軌直線度分段測量拼接方法研究
    軌跡
    軌跡
    天才數(shù)學(xué)家——高斯
    軌跡
    基于西門子840D數(shù)控系統(tǒng)的直線度自動計算程序研究
    進化的軌跡(一)——進化,無盡的適應(yīng)
    中國三峽(2017年2期)2017-06-09 08:15:29
    直線度誤差曲線形成機理與形位特性研究
    重型機械(2016年1期)2016-03-01 03:42:06
    評定導(dǎo)軌直線度誤差的方法
    亚洲三级黄色毛片| 丰满乱子伦码专区| av.在线天堂| 女的被弄到高潮叫床怎么办| 日韩亚洲欧美综合| 国产不卡一卡二| 久久精品夜夜夜夜夜久久蜜豆| 日韩精品青青久久久久久| 婷婷色av中文字幕| 最近最新中文字幕免费大全7| 三级国产精品片| 纵有疾风起免费观看全集完整版 | 久久久久久久久久久免费av| 高清毛片免费看| 亚洲欧美中文字幕日韩二区| 亚洲欧洲日产国产| 国内精品美女久久久久久| 日韩,欧美,国产一区二区三区| 国产大屁股一区二区在线视频| 亚洲国产精品国产精品| 看黄色毛片网站| 一个人看视频在线观看www免费| 免费观看精品视频网站| 秋霞在线观看毛片| 久久97久久精品| 一边亲一边摸免费视频| 97超视频在线观看视频| 少妇人妻一区二区三区视频| 天堂中文最新版在线下载 | 国产精品国产三级专区第一集| 亚洲18禁久久av| 国产成年人精品一区二区| 久久国产乱子免费精品| 日本黄大片高清| 国产免费福利视频在线观看| 在线播放无遮挡| 欧美97在线视频| 亚洲精品国产av蜜桃| 美女高潮的动态| 大香蕉97超碰在线| 五月天丁香电影| 别揉我奶头 嗯啊视频| 又爽又黄a免费视频| av黄色大香蕉| 午夜免费激情av| 99九九线精品视频在线观看视频| 国产不卡一卡二| 91av网一区二区| 汤姆久久久久久久影院中文字幕 | 在线观看美女被高潮喷水网站| 乱人视频在线观看| 成人国产麻豆网| 精品人妻偷拍中文字幕| 成年版毛片免费区| 黄色配什么色好看| 看黄色毛片网站| 亚洲,欧美,日韩| 1000部很黄的大片| 日本与韩国留学比较| 街头女战士在线观看网站| 欧美不卡视频在线免费观看| 男女那种视频在线观看| 午夜久久久久精精品| 日韩电影二区| 成人美女网站在线观看视频| 国产精品麻豆人妻色哟哟久久 | 美女内射精品一级片tv| 熟妇人妻不卡中文字幕| 欧美+日韩+精品| 秋霞伦理黄片| 国产91av在线免费观看| 午夜精品在线福利| 青春草亚洲视频在线观看| 精品一区在线观看国产| 一级片'在线观看视频| 国产黄a三级三级三级人| av在线播放精品| 亚洲精品影视一区二区三区av| 91久久精品国产一区二区成人| 亚洲国产欧美人成| 自拍偷自拍亚洲精品老妇| 免费无遮挡裸体视频| 成年版毛片免费区| 色综合色国产| 人人妻人人看人人澡| 久久精品国产亚洲av涩爱| 久久99精品国语久久久| 精品久久国产蜜桃| 2021少妇久久久久久久久久久| 最近2019中文字幕mv第一页| 国产亚洲av嫩草精品影院| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲欧美精品专区久久| 国产精品一区二区在线观看99 | av免费在线看不卡| 五月玫瑰六月丁香| 免费播放大片免费观看视频在线观看| 国产熟女欧美一区二区| 中文字幕久久专区| 免费人成在线观看视频色| 精品久久久久久久人妻蜜臀av| 国产午夜精品一二区理论片| 国产精品女同一区二区软件| 免费av毛片视频| 少妇被粗大猛烈的视频| 在线天堂最新版资源| 黄色一级大片看看| 人妻制服诱惑在线中文字幕| 欧美日韩一区二区视频在线观看视频在线 | 爱豆传媒免费全集在线观看| 免费黄色在线免费观看| 人人妻人人澡人人爽人人夜夜 | 卡戴珊不雅视频在线播放| 成人无遮挡网站| 搡老妇女老女人老熟妇| 在线观看一区二区三区| 国产 一区精品| 天天躁夜夜躁狠狠久久av| 亚洲精品aⅴ在线观看| 97超视频在线观看视频| 男女视频在线观看网站免费| 婷婷色综合www| 精品一区二区三区视频在线| 亚洲欧美日韩卡通动漫| 黑人高潮一二区| 丰满人妻一区二区三区视频av| 噜噜噜噜噜久久久久久91| 在线观看免费高清a一片| 国产综合精华液| 国产成人a∨麻豆精品| 欧美丝袜亚洲另类| 极品少妇高潮喷水抽搐| 精品国内亚洲2022精品成人| 老师上课跳d突然被开到最大视频| 成年免费大片在线观看| 国产精品久久久久久久久免| 少妇的逼水好多| 夜夜爽夜夜爽视频| 日本三级黄在线观看| 乱码一卡2卡4卡精品| 在线观看人妻少妇| 国产精品人妻久久久久久| 国产精品久久久久久久电影| 26uuu在线亚洲综合色| 日日撸夜夜添| 亚洲av成人av| 高清日韩中文字幕在线| 搡女人真爽免费视频火全软件| 亚洲自偷自拍三级| 青春草亚洲视频在线观看| 国产精品精品国产色婷婷| 日日摸夜夜添夜夜爱| 成年免费大片在线观看| 国产午夜精品一二区理论片| 国产毛片a区久久久久| 天天一区二区日本电影三级| 中文字幕制服av| 久久午夜福利片| 激情五月婷婷亚洲| 少妇裸体淫交视频免费看高清| 能在线免费看毛片的网站| 成人亚洲精品av一区二区| 国产精品久久久久久久电影| 最近中文字幕高清免费大全6| 国产综合懂色| 精品不卡国产一区二区三区| 亚洲怡红院男人天堂| 免费不卡的大黄色大毛片视频在线观看 | 亚洲av国产av综合av卡| 日韩亚洲欧美综合| 五月天丁香电影| 观看免费一级毛片| 国产成人freesex在线| 久久久久国产网址| 高清午夜精品一区二区三区| 久久久久国产网址| 国产精品99久久久久久久久| 内射极品少妇av片p| 国内精品一区二区在线观看| 三级经典国产精品| 国产极品天堂在线| 菩萨蛮人人尽说江南好唐韦庄| 国产成年人精品一区二区| 美女黄网站色视频| 国产成年人精品一区二区| 美女xxoo啪啪120秒动态图| 嫩草影院精品99| 日日干狠狠操夜夜爽| 一本一本综合久久| 成人欧美大片| 日韩欧美 国产精品| 亚洲乱码一区二区免费版| 街头女战士在线观看网站| 亚洲精品自拍成人| 国产真实伦视频高清在线观看| 国产一区亚洲一区在线观看| 97超视频在线观看视频| 日韩强制内射视频| 午夜激情欧美在线| www.av在线官网国产| 男女啪啪激烈高潮av片| 草草在线视频免费看| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久久精品欧美日韩精品| 男插女下体视频免费在线播放| 亚洲精品456在线播放app| 成人综合一区亚洲| 成年版毛片免费区| 美女脱内裤让男人舔精品视频| av播播在线观看一区| 一二三四中文在线观看免费高清| 国产乱人偷精品视频| 亚洲精品日韩在线中文字幕| 禁无遮挡网站| 国产精品1区2区在线观看.| 男女国产视频网站| 乱码一卡2卡4卡精品| 自拍偷自拍亚洲精品老妇| 最近手机中文字幕大全| 精品酒店卫生间| 2018国产大陆天天弄谢| 中文字幕免费在线视频6| av.在线天堂| 国产精品人妻久久久影院| 男女边吃奶边做爰视频| 国产精品一及| 天堂√8在线中文| 国产亚洲午夜精品一区二区久久 | 成人国产麻豆网| 一级二级三级毛片免费看| 精品99又大又爽又粗少妇毛片| 国产精品综合久久久久久久免费| 成人无遮挡网站| av专区在线播放| 老司机影院毛片| 日本wwww免费看| 免费大片18禁| 能在线免费观看的黄片| av女优亚洲男人天堂| 亚洲精品第二区| 亚洲精品成人av观看孕妇| 国产综合懂色| 精品不卡国产一区二区三区| 国产 一区精品| 成人特级av手机在线观看| 一本一本综合久久| 亚洲色图av天堂| 久久久精品94久久精品| av在线老鸭窝| 草草在线视频免费看| 国产精品国产三级国产专区5o| 精品国产露脸久久av麻豆 | 中文精品一卡2卡3卡4更新| 在线免费观看的www视频| 国产一区二区亚洲精品在线观看| av国产久精品久网站免费入址| 日本黄色片子视频| 亚洲三级黄色毛片| .国产精品久久| av国产免费在线观看| 久久久色成人| 国产老妇女一区| 18禁在线播放成人免费| 狠狠精品人妻久久久久久综合| 18+在线观看网站| 亚洲成色77777| 男女啪啪激烈高潮av片| 99久久中文字幕三级久久日本| 嫩草影院精品99| 少妇人妻精品综合一区二区| 精品久久久噜噜| 一个人看的www免费观看视频| 五月玫瑰六月丁香| 久久久a久久爽久久v久久| 久久久久久九九精品二区国产| 777米奇影视久久| 在线观看免费高清a一片| 久99久视频精品免费| 国产综合懂色| 国产探花极品一区二区| eeuss影院久久| 91精品伊人久久大香线蕉| 国产视频内射| 80岁老熟妇乱子伦牲交| 久久精品夜夜夜夜夜久久蜜豆| 久久久久精品性色| 成人午夜精彩视频在线观看| 纵有疾风起免费观看全集完整版 | 国产精品久久久久久精品电影| 亚洲国产欧美人成| 五月天丁香电影| 联通29元200g的流量卡| 精品国产一区二区三区久久久樱花 | 最后的刺客免费高清国语| 久久久久久久久中文| 亚洲国产欧美在线一区| 亚洲在线自拍视频| 丰满人妻一区二区三区视频av| 欧美极品一区二区三区四区| 国产一区二区三区av在线| 午夜激情欧美在线| 亚洲精品视频女| av女优亚洲男人天堂| 看免费成人av毛片| 99热网站在线观看| 青青草视频在线视频观看| 大片免费播放器 马上看| 日韩 亚洲 欧美在线| 精品久久久噜噜| 欧美xxⅹ黑人| 国产精品爽爽va在线观看网站| 国内精品宾馆在线| 亚洲天堂国产精品一区在线| 精品不卡国产一区二区三区| 国产av国产精品国产| 国产又色又爽无遮挡免| 2018国产大陆天天弄谢| 亚洲在线自拍视频| 中文精品一卡2卡3卡4更新| 91精品国产九色| ponron亚洲| 亚洲最大成人av| 亚洲精品第二区| 尾随美女入室| 22中文网久久字幕| 国产高清国产精品国产三级 | 国内精品美女久久久久久| 精品午夜福利在线看| 禁无遮挡网站| 九草在线视频观看| freevideosex欧美| 最后的刺客免费高清国语| 亚洲经典国产精华液单| 久久午夜福利片| 国产白丝娇喘喷水9色精品| 日日干狠狠操夜夜爽| 国产精品99久久久久久久久| 亚洲精品久久午夜乱码| 欧美激情久久久久久爽电影| 91在线精品国自产拍蜜月| 特级一级黄色大片| 日韩av免费高清视频| 中文字幕免费在线视频6| 国产高潮美女av| 久久久久久伊人网av| 日韩成人av中文字幕在线观看| 国产精品一二三区在线看| 国产欧美日韩精品一区二区| 免费大片黄手机在线观看| 亚洲美女视频黄频| 好男人视频免费观看在线| a级一级毛片免费在线观看| 边亲边吃奶的免费视频| 我的女老师完整版在线观看| 男人和女人高潮做爰伦理| 日本熟妇午夜| 激情五月婷婷亚洲| 亚洲国产日韩欧美精品在线观看| 男女国产视频网站| 91狼人影院| 国产精品精品国产色婷婷| 在线播放无遮挡| 成人漫画全彩无遮挡| 国产成人午夜福利电影在线观看| 偷拍熟女少妇极品色| 18+在线观看网站| 久久久精品欧美日韩精品| 日韩在线高清观看一区二区三区| 欧美日韩视频高清一区二区三区二| 赤兔流量卡办理| 一级a做视频免费观看| 极品教师在线视频| 白带黄色成豆腐渣| av播播在线观看一区| av一本久久久久| 最近2019中文字幕mv第一页| 国产精品一区二区三区四区免费观看| 汤姆久久久久久久影院中文字幕 | 日日干狠狠操夜夜爽| 国产精品一二三区在线看| 网址你懂的国产日韩在线| 亚洲av免费在线观看| 国产精品一二三区在线看| 黄色日韩在线| 日韩欧美三级三区| 国产麻豆成人av免费视频| 国产精品1区2区在线观看.| 卡戴珊不雅视频在线播放| 成人亚洲精品一区在线观看 | 男女下面进入的视频免费午夜| 三级毛片av免费| 色网站视频免费| 观看美女的网站| 精品久久久久久久久久久久久| 国产女主播在线喷水免费视频网站 | 51国产日韩欧美| 国产探花极品一区二区| 久久97久久精品| 三级经典国产精品| 午夜福利视频1000在线观看| 国产一区二区三区综合在线观看 | 在线观看一区二区三区| 色综合站精品国产| 中文字幕亚洲精品专区| 亚洲成色77777| 欧美一区二区亚洲| 亚洲不卡免费看| 午夜日本视频在线| 国产在视频线精品| 国产午夜精品论理片| 天美传媒精品一区二区| 国产伦理片在线播放av一区| 一二三四中文在线观看免费高清| 国产熟女欧美一区二区| 国产精品美女特级片免费视频播放器| 日韩av在线免费看完整版不卡| 青春草国产在线视频| 亚洲精品日韩av片在线观看| 亚洲精品456在线播放app| 成人国产麻豆网| 国产亚洲av片在线观看秒播厂 | 大香蕉97超碰在线| 国产精品.久久久| 免费大片18禁| 欧美高清性xxxxhd video| 波野结衣二区三区在线| 一边亲一边摸免费视频| 中文字幕免费在线视频6| av.在线天堂| av在线天堂中文字幕| 亚洲精品一区蜜桃| 久久热精品热| 在现免费观看毛片| 精品酒店卫生间| 丰满少妇做爰视频| 国产爱豆传媒在线观看| 亚洲精品中文字幕在线视频 | 97精品久久久久久久久久精品| 亚洲精品亚洲一区二区| 一本一本综合久久| 天天一区二区日本电影三级| 狂野欧美激情性xxxx在线观看| 久久草成人影院| 亚洲欧美一区二区三区黑人 | 久久久久久久久久久免费av| 午夜亚洲福利在线播放| 三级男女做爰猛烈吃奶摸视频| 日本欧美国产在线视频| 欧美xxⅹ黑人| 三级国产精品片| 久久精品国产自在天天线| 蜜桃久久精品国产亚洲av| 久久国产乱子免费精品| 久久精品夜色国产| 精品不卡国产一区二区三区| 日本免费a在线| 欧美潮喷喷水| 日韩不卡一区二区三区视频在线| 人妻系列 视频| 国产熟女欧美一区二区| 舔av片在线| a级一级毛片免费在线观看| 中文精品一卡2卡3卡4更新| 欧美极品一区二区三区四区| 国产永久视频网站| 久久久色成人| 中文字幕av成人在线电影| 免费看日本二区| 美女国产视频在线观看| 欧美变态另类bdsm刘玥| 一二三四中文在线观看免费高清| 大话2 男鬼变身卡| 美女高潮的动态| 国产三级在线视频| 久久99热6这里只有精品| 国产色爽女视频免费观看| 内射极品少妇av片p| 卡戴珊不雅视频在线播放| 欧美xxⅹ黑人| 一本一本综合久久| 精品一区二区三区人妻视频| 国产在视频线精品| 人体艺术视频欧美日本| 亚洲aⅴ乱码一区二区在线播放| 高清视频免费观看一区二区 | 中文字幕亚洲精品专区| 亚洲欧美精品专区久久| 一级黄片播放器| 免费看光身美女| 成人特级av手机在线观看| 成人综合一区亚洲| 亚洲自偷自拍三级| 日本一本二区三区精品| 99热网站在线观看| 看黄色毛片网站| 色播亚洲综合网| 两个人的视频大全免费| 亚洲精品日本国产第一区| 人妻夜夜爽99麻豆av| 国产亚洲av嫩草精品影院| 一级毛片黄色毛片免费观看视频| 卡戴珊不雅视频在线播放| 亚洲无线观看免费| 久久久久久久午夜电影| 免费观看精品视频网站| 少妇的逼好多水| 丝瓜视频免费看黄片| 国产真实伦视频高清在线观看| 99热全是精品| 寂寞人妻少妇视频99o| 日本爱情动作片www.在线观看| 高清视频免费观看一区二区 | 亚洲图色成人| 欧美bdsm另类| 中文字幕免费在线视频6| 少妇丰满av| 又爽又黄a免费视频| av一本久久久久| 久久久久久久久久久免费av| 亚洲欧美一区二区三区国产| 久久精品夜夜夜夜夜久久蜜豆| 免费人成在线观看视频色| 两个人视频免费观看高清| 人妻夜夜爽99麻豆av| 非洲黑人性xxxx精品又粗又长| 2022亚洲国产成人精品| 久久人人爽人人片av| 婷婷色麻豆天堂久久| 三级毛片av免费| 中文字幕av成人在线电影| 男人和女人高潮做爰伦理| 色网站视频免费| 久久韩国三级中文字幕| 人妻制服诱惑在线中文字幕| 一个人免费在线观看电影| 国产探花极品一区二区| 一夜夜www| 久久鲁丝午夜福利片| 国产av国产精品国产| 黄色配什么色好看| 欧美激情在线99| 成年版毛片免费区| 国产麻豆成人av免费视频| 国产精品av视频在线免费观看| 亚洲在线观看片| 免费大片18禁| 大话2 男鬼变身卡| 国产成人aa在线观看| 日产精品乱码卡一卡2卡三| 欧美丝袜亚洲另类| 亚洲欧美一区二区三区黑人 | 美女大奶头视频| 久久精品久久久久久久性| 国产高清不卡午夜福利| 久久精品久久精品一区二区三区| 日韩,欧美,国产一区二区三区| 性插视频无遮挡在线免费观看| 欧美最新免费一区二区三区| 精品久久久久久久末码| 日本黄大片高清| 亚洲成人一二三区av| 中文精品一卡2卡3卡4更新| 国产精品人妻久久久久久| 国产成人精品婷婷| 男女那种视频在线观看| 精品久久久久久久久av| 成人美女网站在线观看视频| 97超视频在线观看视频| 亚洲在线观看片| 国产男人的电影天堂91| 尤物成人国产欧美一区二区三区| 国产伦精品一区二区三区视频9| 丝瓜视频免费看黄片| 又大又黄又爽视频免费| 高清在线视频一区二区三区| av又黄又爽大尺度在线免费看| 黄色欧美视频在线观看| 久久久国产一区二区| 亚洲精华国产精华液的使用体验| 欧美性猛交╳xxx乱大交人| 日韩国内少妇激情av| 97热精品久久久久久| 免费高清在线观看视频在线观看| 十八禁国产超污无遮挡网站| 免费观看性生交大片5| 国产乱人视频| 国产亚洲av片在线观看秒播厂 | 狠狠精品人妻久久久久久综合| 色综合亚洲欧美另类图片| 成人美女网站在线观看视频| 色哟哟·www| 久久这里有精品视频免费| 99九九线精品视频在线观看视频| 日韩一区二区视频免费看| 麻豆国产97在线/欧美| a级一级毛片免费在线观看| 亚洲内射少妇av| 欧美高清成人免费视频www| 青青草视频在线视频观看| www.av在线官网国产| 最近手机中文字幕大全| 边亲边吃奶的免费视频| 国产乱来视频区| 最新中文字幕久久久久| 天堂影院成人在线观看| 少妇的逼好多水| 免费看光身美女| 最近的中文字幕免费完整| 美女国产视频在线观看| 伊人久久精品亚洲午夜| 午夜爱爱视频在线播放| 国产 亚洲一区二区三区 | 成人一区二区视频在线观看| 亚洲国产最新在线播放| 一区二区三区高清视频在线|