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

    三軸磁傳感器系統(tǒng)的在線標(biāo)定*

    2015-12-26 05:58:25蔚,蔣
    彈箭與制導(dǎo)學(xué)報 2015年3期

    盛 蔚,蔣 勇

    (北京航空航天大學(xué)儀器科學(xué)與光電工程學(xué)院,北京 100191)

    三軸磁傳感器系統(tǒng)的在線標(biāo)定*

    盛蔚,蔣勇

    (北京航空航天大學(xué)儀器科學(xué)與光電工程學(xué)院,北京100191)

    摘要:現(xiàn)有的三軸磁傳感器的標(biāo)定方法對于設(shè)備的要求過高,系統(tǒng)不能夠獨立完成標(biāo)定工作。針對此種情況,提出三軸磁傳感器系統(tǒng)的在線標(biāo)定算法,建立三軸磁傳感器系統(tǒng)誤差模型,采用橢球曲面擬合三軸磁場數(shù)據(jù),利用遞推法計算標(biāo)定參數(shù)。達成了傳感器系統(tǒng)本身在線獨立完成標(biāo)定的目的,并對在線標(biāo)定參數(shù)的可信度提出了評價方法。通過實驗驗證了在線標(biāo)定方法的有效性和在線標(biāo)定參數(shù)可信度評價方法的有效性。

    關(guān)鍵詞:磁傳感器;在線標(biāo)定;橢球擬合;遞推最小二乘法

    0引言

    三軸磁傳感器系統(tǒng)是用來測量地球磁場在載體坐標(biāo)系下三軸磁場分量大小的傳感器系統(tǒng)。既可以與加速度計信息組合起來計算載體姿態(tài)[1],又可以當(dāng)作測量磁場異常的傳感器系統(tǒng)用于磁異常定位[2]。根據(jù)三軸磁傳感器的輸出模型,可將外界干擾等效成三軸磁傳感器系統(tǒng)三軸磁場數(shù)據(jù)的零偏誤差和標(biāo)度因數(shù)。目前對于三軸磁傳感器系統(tǒng)的標(biāo)定,主要有陳偉等人提出的利用GPS信息作為參考進行標(biāo)定[3],郭鵬飛等人提出的十二位置不對北標(biāo)定法[4],以及朱建良等人提出的橢球曲面擬合法[5]。這些方法都有各自的缺點,極大的限制了三軸磁傳感器系統(tǒng)在野外的使用,同時也增加了三軸磁傳感器系統(tǒng)的使用成本。文中提出的三軸磁傳感器系統(tǒng)在線標(biāo)定算法,繼承了橢球曲面擬合標(biāo)定法的靈活性以及無需轉(zhuǎn)臺就能進行等優(yōu)點,并針對橢球曲面擬合算法無法獨立完成計算的缺點,提出了使用遞推算法計算標(biāo)定過程中需要用到的橢球系數(shù)的方法。使三軸磁傳感器系統(tǒng)能夠獨立完成標(biāo)定工作,降低了使用成本,也使其使用范圍更加廣泛。

    1誤差模型

    地磁場的干擾主要來源于周圍環(huán)境中的磁性材料產(chǎn)生的磁場。在地磁場和磁性材料產(chǎn)生的干擾磁場共同作用下,三軸磁傳感器在載體坐標(biāo)系下的輸出模型[6]可表示為:

    (1)

    由地磁場的特性,可以認(rèn)為在一個小區(qū)域內(nèi),地磁場的強度值為固定值。在理想條件下,三軸磁傳感器系統(tǒng)可以測得地磁場在載體坐標(biāo)系下的三軸磁場分量的值。以三軸磁場強度分量為坐標(biāo)的測量點將構(gòu)成一個以地磁場強度為半徑、球心在原點的球面。在數(shù)學(xué)上的表示為:

    (2)

    上述情況是在地磁場未受到外界環(huán)境中磁場干擾的理想條件下才能實現(xiàn)的。事實上,由于存在硬磁干擾和軟磁干擾,再加上三軸磁傳感器系統(tǒng)軸向間存在的非正交性和磁阻傳感器本身存在的零偏誤差,上述球面將畸變成一個橢球面。通過旋轉(zhuǎn)三軸磁傳感器系統(tǒng),獲取足夠多的橢球面上的點,就可以通過擬合算法計算出橢球面方程的系數(shù)參數(shù)。利用這些系數(shù)參數(shù)計算出k和u0,對磁場測量值進行補償,消除三軸磁傳感器系統(tǒng)受到的干擾。

    2標(biāo)定參數(shù)的計算

    將式(2)展開,可得到如下二次標(biāo)準(zhǔn)型方程:

    (3)

    橢球曲面方程的一般表達式可用下式表示:

    (4)

    將橢球曲面方程整理為矢量型式:

    (5)

    對比式(3)和式(5)可得:

    (6)

    由上式可知,只要知道了橢球曲面方程中的9個系數(shù)參數(shù)就能夠計算出標(biāo)定參數(shù)的值,然后通過標(biāo)定參數(shù)對磁場測量值進行補償。橢球曲面的擬合方法有牛頓迭代法等,這些方法的共同缺點就是需要存儲大量的數(shù)據(jù)進行計算,需要其他設(shè)備協(xié)助完成標(biāo)定參數(shù)的計算,限制了其在嵌入式系統(tǒng)中的應(yīng)用。為此文中提出了一種在線標(biāo)定三軸磁傳感器系統(tǒng)的方法——利用遞推最小二乘法對橢球曲面方程進行擬合計算。解決了三軸磁傳感器系統(tǒng)標(biāo)定必須依賴其他設(shè)備的缺點

    3擬合橢球曲面方程參數(shù)的計算

    在采用遞推最小二乘法計算擬合橢球曲面參數(shù)的過程中,有可能會出現(xiàn)系數(shù)全為0的奇點[7],而且受到干擾的磁場測量數(shù)據(jù)還有可能造成擬合橢球曲面發(fā)散為雙曲面。因此需要對遞推過程施加約束條件。不等式形式的約束條件會使參數(shù)的計算變得極為困難。為此,Fitzgibbon提出在二維條件下約束不等式[8]可寫為4ac-b2=1。Nikos Grammalidis等人將其原理推廣到三維條件下[9],即橢球曲面擬合的約束條件為:

    在不考慮約束條件的情況下,將前文的橢球曲面方程看作是參數(shù)待定的二次曲面方程,若測得N組測量數(shù)據(jù),則可寫出N個方程,寫成矩陣形式為:

    (7)

    將上式簡寫為:tN=φNθ。

    根據(jù)遞推最小二乘估計理論的一般公式,θ的遞推最小二乘估計公式如下:

    (8)

    式(8)為不考慮約束條件下的遞推公式。對于帶線性約束條件的遞推最小二乘法,只需要在不帶約束的最小二乘法的基礎(chǔ)上加上一個修正項即可。

    (9)

    式(8)在加上約束條件下的補償公式為式(9),即為約束條件下橢球曲面參數(shù)計算的遞推公式。將采集到的每組數(shù)據(jù)代入到遞推計算公式中進行計算,直到最后一組數(shù)據(jù)遞推計算完畢,此時θ的最小二乘估計值即為求得的橢球曲面方程的系數(shù)參數(shù)值。由橢球曲面參數(shù)系數(shù)即可求得三軸磁傳感器系統(tǒng)的標(biāo)定參數(shù)。

    4標(biāo)定質(zhì)量的評價

    對三軸磁傳感器系統(tǒng)的標(biāo)定過程實際上是為了確定上述橢球方程的10個系數(shù)參數(shù)。橢球參數(shù)的計算是基于地磁場強度為固定值的條件,而在標(biāo)定過程中也無法分辨磁場的來源是地磁場還是其他物體產(chǎn)生的磁場。因此,在三軸磁傳感器系統(tǒng)標(biāo)定時需要一個相對穩(wěn)定的磁場環(huán)境。若外界存在變化的磁場,標(biāo)定的精度將受到嚴(yán)重的影響,標(biāo)定計算得到的標(biāo)定參數(shù)也會失效。為了判斷標(biāo)定結(jié)果的可信程度,選取遞推計算過程中磁場總強度的計算值的平方mTm作為判斷標(biāo)定過程中磁場環(huán)境是否穩(wěn)定參考變量。在穩(wěn)定的磁場環(huán)境中,mTm是一個逐漸收斂的值。遞推得到的參數(shù)越接近最終結(jié)果,mTm的值也會越穩(wěn)定,最終趨向于一個穩(wěn)定的值。如果在標(biāo)定過程中受到變化的磁場的干擾,mTm的值將不能收斂,這就在一定程度上判斷出標(biāo)定結(jié)果的可信程度。

    在三軸磁傳感器系統(tǒng)的標(biāo)定過程中,每次遞推計算后都能得到一組擬合橢球曲面的參數(shù),繼而可以計算出標(biāo)定參數(shù)k和u0。由式(2)可得此時磁場強度的計算值,寫成三軸分量形式為:

    (10)

    由磁場的三軸分量計算出磁場強度的平方mTm。在標(biāo)定結(jié)束時,取mTm的最后100個計算值求其方差V。

    5實驗分析

    用于實驗的三軸磁傳感器系統(tǒng)采用Honeywell磁阻傳感器,具有較高的磁場分辨率和穩(wěn)定的磁場輸出。將三軸磁傳感器系統(tǒng)安裝到六面體工裝上,任意旋轉(zhuǎn)變換三軸磁傳感器系統(tǒng)姿態(tài)。保證采集到的磁場數(shù)據(jù)在擬合橢球曲面的每個方向上都有分布。旋轉(zhuǎn)工裝的速度要緩慢,讓每個位置都能采集到足夠的數(shù)據(jù),如果某個方向的數(shù)據(jù)過少,很可能會影響到標(biāo)定的精度。記錄采集到3 000組數(shù)據(jù)并繪圖。

    在未受到干擾的條件下,采集到的磁場三軸分量測量值數(shù)據(jù)繪圖如圖1所示。

    圖1 無干擾條件下三軸磁場數(shù)據(jù)情況

    在標(biāo)定過程中采集到的每組數(shù)據(jù)作為新息進行遞推計算。在未受到干擾的條件下,遞推計算過程中的磁場總強度平方的計算值繪圖如圖2所示。

    圖2 無干擾條件下磁場總強度變化曲線平方

    在三軸磁傳感器系統(tǒng)的標(biāo)定過程中,在一旁放置一塊磁鐵并不斷改變磁鐵的位置,用于施加一個變化的干擾磁場。在這種條件下,標(biāo)定過程中采集到的數(shù)據(jù)繪圖如圖3所示。

    圖3 干擾條件下三軸磁場數(shù)據(jù)情況

    在外界磁場不穩(wěn)定的條件下,標(biāo)定過程中磁場總強度平方的計算值繪圖如圖4所示。

    圖4 干擾條件下磁場總強度平方變化曲線

    為了驗證三軸磁傳感器系統(tǒng)的標(biāo)定效果,利用其測量得到的三軸磁場值計算出載體的航向。將三軸磁傳感器系統(tǒng)安裝在無磁轉(zhuǎn)臺上,以其初始位置為起始位置,每旋轉(zhuǎn)30°記錄一次角度誤差數(shù)據(jù),直到旋轉(zhuǎn)一周。記錄數(shù)據(jù)繪圖如圖5所示。

    圖5 三種條件下角度誤差曲線

    上圖中數(shù)據(jù)顯示三軸磁傳感器系統(tǒng)測量得到的航向輸出誤差由18.7°減小到0.7°,表明三軸磁傳感器系統(tǒng)在經(jīng)過在線標(biāo)定算法標(biāo)定后,航向輸出誤差顯著減小。在線標(biāo)定算法成功的補償了外界干擾對磁場測量的影響。

    在未受到干擾的條件下,遞推計算過程中的磁場總強度平方的計算值如圖2所示,此時V的值為1.71×10-5。在外界磁場不穩(wěn)定的條件下,標(biāo)定過程中磁場總強度平方的計算值如圖4所示,此時V的值為0.001 5。

    經(jīng)過大量的實驗數(shù)據(jù)分析,V的值越小,標(biāo)定參數(shù)的可靠性越高。在V的值小于1.0×10-4的情況下,標(biāo)定參數(shù)在誤差允許的范圍內(nèi)都是可靠的。V的值在一定程度上表明了標(biāo)定結(jié)果的可信程度。

    6結(jié)論

    利用地磁場在小范圍內(nèi)短時間不變的特性,建立起橢球曲面模型對三軸磁傳感器系統(tǒng)進行補償,是一種十分有效的標(biāo)定方法。但是由于其計算時需要大量的數(shù)據(jù),限制了其應(yīng)用范圍。文中提出的以帶約束的遞推最小二乘法計算擬合橢球曲面參數(shù)的在線標(biāo)定方法,克服了這種缺點,同時保證了在計算過程中參數(shù)的穩(wěn)定性。相當(dāng)于將整體計算擬合橢球曲面參數(shù)的過程拆分成一步一步進行,避免了計算參數(shù)時對其他設(shè)備的依賴,在避免標(biāo)定精度降低的條件下實現(xiàn)了三軸磁傳感器系統(tǒng)的在線標(biāo)定,擴大了其使用范圍,具有很強的實用性。同時文中提出在線標(biāo)定參數(shù)質(zhì)量的評價方法給出了一種標(biāo)定參數(shù)有效性的初步評價標(biāo)準(zhǔn),能夠及時的發(fā)現(xiàn)標(biāo)定過程中是否受到干擾。實驗結(jié)果證明了在線標(biāo)定算法的有效性以及標(biāo)定參數(shù)質(zhì)量評價的有效性。

    參考文獻:

    [1]Qi Zhang, Qixiao Gao, Yongbing Chen. A novel magnetic compass calibration method based on improved ellipse model [C]∥International Conference on Intelligent Control and Information Processing, 2010: 11-15.

    [2]Chao Hu, Tongxing Ma. Sensor arrangement optimization of magnetic localization and orientation system [C]∥ Proceedings of the 2007 IEEE International Conference on Integration Technology, 2007: 311-315.

    [3]Wei Chen, Ruizhi Chen, Yuwei Chen. An adaptive calibration approach for a 2-axis digital compass in a low-cost pedestrian navigation system [C]∥Proceedings of the 2010 IEEE International Conference on Instrumentation and Measurement Technology, 2010: 1392-1397.

    [4]郭鵬飛, 任章, 邱海韜. 一種十二位置不對北的磁羅盤標(biāo)定方法 [J]. 中國慣性技術(shù)學(xué)報, 2007, 15(5): 598-601.

    [5]朱建良, 王興全, 吳盤龍. 基于橢球曲面擬合的三位磁羅盤誤差補償算法 [J]. 中國慣性技術(shù)學(xué)報, 2012, 20(5): 562-566.

    [6]范玉寶, 李杰, 張曉明. 基于橢球擬合的微慣性測量組合現(xiàn)場快速標(biāo)定方法 [J]. 傳感技術(shù)學(xué)報, 2011, 24(11): 1547-1550.

    [7]孫宏偉, 房建成, 李艷. 橢圓擬合方法在磁羅盤羅差校準(zhǔn)中的應(yīng)用 [J]. 光學(xué)精密工程, 2009, 17(12): 3034-3038.

    [8]Andrewfitzgibbon, Mauriziopilu, Robertb. DirectLeast Square Fitting of Ellipses [J]. Pattern Analysis And Machine Intelligence, 1999, 21(5): 476-480.

    [9]Nikosgrammalidis, Michaelg Strintzis. Direct least square fitting of ellipsoids [C]∥International Conference on Pattern Recognition, 2012: 3228-3231. Zhang, Qixiao Gao, Yongbing Chen. A novel magnetic compass calibration method based on improved ellipse model [C]∥International Conference on Intelligent Control and Information Processing, 2010: 11-15.

    [2]Chao Hu, Tongxing Ma. Sensor arrangement optimization of magnetic localization and orientation system [C]∥ Proceedings of the 2007 IEEE International Conference on Integration Technology, 2007: 311-315.

    [3]Wei Chen, Ruizhi Chen, Yuwei Chen. An adaptive calibration approach for a 2-axis digital compass in a low-cost pedestrian navigation system [C]∥Proceedings of the 2010 IEEE International Conference on Instrumentation and Measurement Technology, 2010: 1392-1397.

    [4]郭鵬飛, 任章, 邱海韜. 一種十二位置不對北的磁羅盤標(biāo)定方法 [J]. 中國慣性技術(shù)學(xué)報, 2007, 15(5): 598-601.

    [5]朱建良, 王興全, 吳盤龍. 基于橢球曲面擬合的三位磁羅盤誤差補償算法 [J]. 中國慣性技術(shù)學(xué)報, 2012, 20(5): 562-566.

    [6]范玉寶, 李杰, 張曉明. 基于橢球擬合的微慣性測量組合現(xiàn)場快速標(biāo)定方法 [J]. 傳感技術(shù)學(xué)報, 2011, 24(11): 1547-1550.

    [7]孫宏偉, 房建成, 李艷. 橢圓擬合方法在磁羅盤羅差校準(zhǔn)中的應(yīng)用 [J]. 光學(xué)精密工程, 2009, 17(12): 3034-3038.

    [8]Andrewfitzgibbon, Mauriziopilu, Robertb. DirectLeast Square Fitting of Ellipses [J]. Pattern Analysis And Machine Intelligence, 1999, 21(5): 476-480.

    [9]Nikosgrammalidis, Michaelg Strintzis. Direct least square fitting of ellipsoids [C]∥International Conference on Pattern Recognition, 2012: 3228-3231.

    收稿日期:2014-06-12

    基金項目:國家自然科學(xué)基金(61274117);國防技術(shù)基礎(chǔ)科研(JSJC2013601C005)資助

    作者簡介:盛蔚(1972-),女,陜西西安人,教授,研究方向:硅MEMS慣性傳感器與微小型導(dǎo)航、制導(dǎo)與控制系統(tǒng)。

    中圖分類號:V241.61

    文獻標(biāo)志碼:A

    An On-line Calibration Method of Three-axis Magnetic Sensor System

    SHENG Wei,JIANG Yong

    (School of Instrument Science and Opto-electronics Engineering, BUAA, Beijing 100191, China)

    Abstract:Available calibration methods for three-axis magnetic sensor systems overclaim for equipment, systems cannot complete calibration independently. Regarding the situation, an online calibration method was proposed. With this method, three-axis magnetic sensor system could complete calibration independently. A method for evaluating credibility of the online calibration parameters was also put forward in the paper. The validity of the online calibration method and online calibration parameter evaluation method was proved to be effective by a series of experiments.

    Keywords:magnetic sensor; online calibration; ellipsoid fitting; recursive least-squares method

    久久精品人妻少妇| 90打野战视频偷拍视频| 成人无遮挡网站| 女人被狂操c到高潮| 精品国内亚洲2022精品成人| 久久久久久人人人人人| 亚洲av免费在线观看| 男女视频在线观看网站免费| tocl精华| 成人18禁在线播放| 麻豆成人av在线观看| 久久人人精品亚洲av| 国产av不卡久久| 成人午夜高清在线视频| 可以在线观看的亚洲视频| e午夜精品久久久久久久| 99久国产av精品| 日本撒尿小便嘘嘘汇集6| 精品久久久久久成人av| 国产精品99久久99久久久不卡| 少妇人妻一区二区三区视频| 香蕉丝袜av| 精品无人区乱码1区二区| 啦啦啦观看免费观看视频高清| 国产精品亚洲一级av第二区| 国产高清激情床上av| 黄频高清免费视频| 午夜福利成人在线免费观看| 国产精品一区二区三区四区免费观看 | 免费在线观看成人毛片| 成人亚洲精品av一区二区| 少妇的丰满在线观看| 老熟妇乱子伦视频在线观看| 国产极品精品免费视频能看的| 国产爱豆传媒在线观看| 久久人人精品亚洲av| 亚洲熟妇熟女久久| 高清毛片免费观看视频网站| 最新在线观看一区二区三区| 亚洲精品456在线播放app | 午夜免费激情av| 精品欧美国产一区二区三| 在线看三级毛片| 精品久久久久久成人av| 精品久久久久久久久久免费视频| 日日夜夜操网爽| 少妇的逼水好多| 国产伦一二天堂av在线观看| 熟女电影av网| 午夜免费成人在线视频| 国内揄拍国产精品人妻在线| 亚洲av美国av| 91在线精品国自产拍蜜月 | 美女cb高潮喷水在线观看 | 国产一区二区三区在线臀色熟女| 又爽又黄无遮挡网站| 久久久久免费精品人妻一区二区| 一级a爱片免费观看的视频| 日日干狠狠操夜夜爽| 欧美中文日本在线观看视频| 亚洲最大成人中文| 久久精品国产清高在天天线| 久9热在线精品视频| 久久久精品大字幕| 此物有八面人人有两片| 一夜夜www| 免费看日本二区| 国产午夜精品论理片| 久久国产乱子伦精品免费另类| 听说在线观看完整版免费高清| 免费在线观看亚洲国产| 99国产精品一区二区三区| 久久久国产欧美日韩av| 免费人成视频x8x8入口观看| 在线播放国产精品三级| 悠悠久久av| 久久久久久国产a免费观看| 日韩有码中文字幕| 国产精品亚洲美女久久久| 成人三级黄色视频| 这个男人来自地球电影免费观看| 99热精品在线国产| 亚洲真实伦在线观看| 日韩欧美免费精品| 一本精品99久久精品77| 999久久久国产精品视频| 久久久国产精品麻豆| 精品无人区乱码1区二区| 亚洲无线观看免费| 精品久久久久久,| 免费观看人在逋| 成人三级做爰电影| 亚洲 欧美 日韩 在线 免费| АⅤ资源中文在线天堂| 色视频www国产| 熟妇人妻久久中文字幕3abv| 他把我摸到了高潮在线观看| 久久性视频一级片| 亚洲av成人一区二区三| 国产精品av视频在线免费观看| 最近最新免费中文字幕在线| 日日夜夜操网爽| 人人妻,人人澡人人爽秒播| 精品久久久久久久毛片微露脸| 操出白浆在线播放| 亚洲成av人片在线播放无| 精品久久久久久,| 亚洲黑人精品在线| 国内精品久久久久精免费| 欧美一区二区精品小视频在线| 女生性感内裤真人,穿戴方法视频| 国产欧美日韩一区二区精品| 亚洲精品乱码久久久v下载方式 | 亚洲第一电影网av| www.熟女人妻精品国产| www国产在线视频色| 国产av一区在线观看免费| 精品久久久久久成人av| 免费在线观看视频国产中文字幕亚洲| 老汉色∧v一级毛片| 波多野结衣高清作品| 成年人黄色毛片网站| 精品99又大又爽又粗少妇毛片 | 黄色成人免费大全| 一本综合久久免费| av天堂在线播放| 亚洲成人中文字幕在线播放| 欧美日韩瑟瑟在线播放| 亚洲 欧美 日韩 在线 免费| 中文字幕人成人乱码亚洲影| 99国产综合亚洲精品| x7x7x7水蜜桃| 午夜日韩欧美国产| av在线蜜桃| 久久精品人妻少妇| 欧美一级毛片孕妇| 99精品久久久久人妻精品| 一个人免费在线观看电影 | 一级毛片女人18水好多| 黄色成人免费大全| 日韩欧美精品v在线| 床上黄色一级片| 两性夫妻黄色片| 欧美日韩亚洲国产一区二区在线观看| 国产伦一二天堂av在线观看| 两个人的视频大全免费| 久久久成人免费电影| 最近最新中文字幕大全免费视频| 午夜免费激情av| 成人三级做爰电影| 久久久久久久久中文| 国产视频一区二区在线看| 免费看a级黄色片| 亚洲av美国av| 色精品久久人妻99蜜桃| 欧美黄色淫秽网站| 午夜精品一区二区三区免费看| 亚洲国产精品合色在线| 国产av在哪里看| 一进一出好大好爽视频| 叶爱在线成人免费视频播放| 黄片大片在线免费观看| 91久久精品国产一区二区成人 | 国产精品精品国产色婷婷| 18禁裸乳无遮挡免费网站照片| 久99久视频精品免费| 久久午夜综合久久蜜桃| 色老头精品视频在线观看| 国产99白浆流出| 日韩欧美一区二区三区在线观看| 少妇裸体淫交视频免费看高清| 国产成人aa在线观看| 亚洲在线自拍视频| 欧美不卡视频在线免费观看| 亚洲无线在线观看| 久久久久精品国产欧美久久久| 一区二区三区高清视频在线| 91在线精品国自产拍蜜月 | 不卡一级毛片| 成在线人永久免费视频| 久久人妻av系列| 欧美+亚洲+日韩+国产| 性色av乱码一区二区三区2| 在线播放国产精品三级| 在线视频色国产色| 亚洲成人久久爱视频| 熟妇人妻久久中文字幕3abv| 国产一区二区在线av高清观看| 91麻豆av在线| 日韩免费av在线播放| 91在线精品国自产拍蜜月 | 午夜视频精品福利| 亚洲国产精品成人综合色| 黑人欧美特级aaaaaa片| 美女午夜性视频免费| 亚洲成av人片免费观看| 757午夜福利合集在线观看| 亚洲,欧美精品.| 免费在线观看成人毛片| 国内精品美女久久久久久| 久久中文字幕一级| 国产精品电影一区二区三区| 丰满的人妻完整版| 女人被狂操c到高潮| 九九在线视频观看精品| 国产黄色小视频在线观看| 亚洲人成网站高清观看| 日本 av在线| 麻豆国产av国片精品| 亚洲一区二区三区不卡视频| 免费高清视频大片| 亚洲欧洲精品一区二区精品久久久| 精品电影一区二区在线| a在线观看视频网站| 国产熟女xx| 他把我摸到了高潮在线观看| 亚洲欧美精品综合久久99| 国产高清视频在线播放一区| 又大又爽又粗| 亚洲最大成人中文| 天天添夜夜摸| 一级毛片女人18水好多| 色综合亚洲欧美另类图片| 精品一区二区三区四区五区乱码| 精品国产超薄肉色丝袜足j| 国产黄色小视频在线观看| 天堂网av新在线| 少妇人妻一区二区三区视频| 精品无人区乱码1区二区| 精品不卡国产一区二区三区| 国产成人精品久久二区二区91| 久久久精品大字幕| 久久99热这里只有精品18| 黄色片一级片一级黄色片| 欧美日韩精品网址| 香蕉丝袜av| 夜夜躁狠狠躁天天躁| 91av网一区二区| 别揉我奶头~嗯~啊~动态视频| 又黄又爽又免费观看的视频| 国产精品自产拍在线观看55亚洲| 久久久精品大字幕| 久久亚洲精品不卡| 又爽又黄无遮挡网站| 12—13女人毛片做爰片一| 一个人观看的视频www高清免费观看 | 在线十欧美十亚洲十日本专区| АⅤ资源中文在线天堂| 男女那种视频在线观看| 丰满的人妻完整版| 亚洲国产中文字幕在线视频| 2021天堂中文幕一二区在线观| 一本一本综合久久| 99视频精品全部免费 在线 | 狂野欧美白嫩少妇大欣赏| 特级一级黄色大片| 听说在线观看完整版免费高清| 亚洲人成网站高清观看| 国产欧美日韩一区二区三| 久久精品国产99精品国产亚洲性色| 麻豆久久精品国产亚洲av| 国产成人精品久久二区二区91| 亚洲中文日韩欧美视频| 国产精品一区二区三区四区久久| 国产av一区在线观看免费| 国产男靠女视频免费网站| 午夜福利视频1000在线观看| 精品不卡国产一区二区三区| 国产精品 国内视频| 欧美日韩国产亚洲二区| 亚洲av日韩精品久久久久久密| 色尼玛亚洲综合影院| 国产男靠女视频免费网站| 精品欧美国产一区二区三| 国产精品一区二区精品视频观看| 精品免费久久久久久久清纯| 久久精品亚洲精品国产色婷小说| 国产精品自产拍在线观看55亚洲| 亚洲中文日韩欧美视频| 亚洲国产欧美人成| 九色国产91popny在线| 99re在线观看精品视频| av片东京热男人的天堂| 好男人电影高清在线观看| 国产亚洲精品久久久com| 亚洲最大成人中文| 黄频高清免费视频| 欧美乱色亚洲激情| 中文字幕熟女人妻在线| 99久久精品一区二区三区| 91麻豆精品激情在线观看国产| 国产成人欧美在线观看| 午夜影院日韩av| 亚洲成人久久爱视频| 久久香蕉精品热| 久久久久久久精品吃奶| 不卡av一区二区三区| 琪琪午夜伦伦电影理论片6080| 国产又黄又爽又无遮挡在线| 99久久精品热视频| 日本在线视频免费播放| 最新中文字幕久久久久 | 免费在线观看影片大全网站| 18禁裸乳无遮挡免费网站照片| 国产精品国产高清国产av| 久久久精品欧美日韩精品| 怎么达到女性高潮| 天堂√8在线中文| 国产又黄又爽又无遮挡在线| 免费电影在线观看免费观看| 久久午夜综合久久蜜桃| 免费电影在线观看免费观看| 最新美女视频免费是黄的| 欧美在线黄色| 午夜福利成人在线免费观看| 在线观看一区二区三区| 美女扒开内裤让男人捅视频| 午夜福利欧美成人| 老司机福利观看| av天堂在线播放| bbb黄色大片| 日本黄大片高清| 亚洲五月婷婷丁香| 中文亚洲av片在线观看爽| 亚洲aⅴ乱码一区二区在线播放| 无限看片的www在线观看| or卡值多少钱| 国产成人精品久久二区二区免费| 精品久久久久久久毛片微露脸| 国产av麻豆久久久久久久| 五月玫瑰六月丁香| 亚洲在线自拍视频| 天堂√8在线中文| 淫妇啪啪啪对白视频| 国产91精品成人一区二区三区| 欧美另类亚洲清纯唯美| 国产一区在线观看成人免费| 中出人妻视频一区二区| 国产男靠女视频免费网站| 精品无人区乱码1区二区| 国产乱人视频| 国产成人av教育| 一区二区三区高清视频在线| bbb黄色大片| 午夜福利高清视频| 国产麻豆成人av免费视频| 久久久久久九九精品二区国产| 亚洲色图 男人天堂 中文字幕| 亚洲熟女毛片儿| 又大又爽又粗| 亚洲成人免费电影在线观看| 久久午夜亚洲精品久久| 999久久久精品免费观看国产| 在线观看午夜福利视频| 国产精品亚洲一级av第二区| 两性夫妻黄色片| 成年女人看的毛片在线观看| 好男人在线观看高清免费视频| 成人特级黄色片久久久久久久| 日韩高清综合在线| 亚洲人成网站在线播放欧美日韩| 国产精华一区二区三区| 两个人的视频大全免费| 狂野欧美白嫩少妇大欣赏| 伊人久久大香线蕉亚洲五| 狂野欧美白嫩少妇大欣赏| 大型黄色视频在线免费观看| 欧美色欧美亚洲另类二区| 日本五十路高清| 两人在一起打扑克的视频| 亚洲精品一区av在线观看| 欧美色欧美亚洲另类二区| 人人妻人人澡欧美一区二区| 欧美黄色淫秽网站| 少妇丰满av| 色综合欧美亚洲国产小说| 亚洲欧美日韩无卡精品| 亚洲男人的天堂狠狠| 久久亚洲精品不卡| 国产aⅴ精品一区二区三区波| 日韩欧美三级三区| 窝窝影院91人妻| 色综合站精品国产| 香蕉久久夜色| 99国产精品一区二区蜜桃av| 一级黄色大片毛片| 真人一进一出gif抽搐免费| 99久久久亚洲精品蜜臀av| 久99久视频精品免费| 97碰自拍视频| 日本黄色片子视频| 在线观看一区二区三区| 欧美大码av| 在线免费观看的www视频| 久久国产精品人妻蜜桃| 亚洲国产欧洲综合997久久,| 免费一级毛片在线播放高清视频| 国产一区二区三区在线臀色熟女| 久久精品亚洲精品国产色婷小说| cao死你这个sao货| 国产精品久久久久久久电影 | 淫妇啪啪啪对白视频| АⅤ资源中文在线天堂| 9191精品国产免费久久| 欧美日韩乱码在线| 国产精品1区2区在线观看.| 日本 欧美在线| 精品不卡国产一区二区三区| 一本一本综合久久| 国产精品久久久久久亚洲av鲁大| 性欧美人与动物交配| 91av网一区二区| 亚洲男人的天堂狠狠| 一进一出好大好爽视频| 免费在线观看日本一区| 90打野战视频偷拍视频| 国产精品乱码一区二三区的特点| 99久久无色码亚洲精品果冻| 老司机深夜福利视频在线观看| 国产三级黄色录像| 久久人妻av系列| 国产蜜桃级精品一区二区三区| 精品日产1卡2卡| 精品人妻1区二区| 久久久色成人| 国产毛片a区久久久久| 成人特级av手机在线观看| 国产成人啪精品午夜网站| 免费在线观看影片大全网站| 亚洲国产精品成人综合色| 国产精品99久久99久久久不卡| 久久精品91无色码中文字幕| 999精品在线视频| 亚洲国产色片| 色在线成人网| 欧美中文综合在线视频| 黄片小视频在线播放| 亚洲在线观看片| 国产精品日韩av在线免费观看| 国产亚洲精品综合一区在线观看| 亚洲第一电影网av| 草草在线视频免费看| 国产精品一区二区免费欧美| 久久久久久大精品| 老司机深夜福利视频在线观看| 夜夜爽天天搞| 免费在线观看影片大全网站| 国产97色在线日韩免费| avwww免费| av黄色大香蕉| 哪里可以看免费的av片| 国产午夜精品论理片| 亚洲黑人精品在线| 三级国产精品欧美在线观看 | 一区福利在线观看| 欧美一级a爱片免费观看看| 综合色av麻豆| 欧美精品啪啪一区二区三区| 美女cb高潮喷水在线观看 | 亚洲男人的天堂狠狠| 美女黄网站色视频| 成人欧美大片| 精品国产乱码久久久久久男人| 精品福利观看| tocl精华| 亚洲自拍偷在线| 国产精品永久免费网站| 国产午夜精品久久久久久| 国产精品亚洲美女久久久| 最近最新免费中文字幕在线| ponron亚洲| www国产在线视频色| 国产欧美日韩一区二区三| 丁香六月欧美| 99国产精品一区二区蜜桃av| 身体一侧抽搐| 欧美色欧美亚洲另类二区| 最新中文字幕久久久久 | 国产蜜桃级精品一区二区三区| 看免费av毛片| 制服丝袜大香蕉在线| 久久国产精品影院| 国产成人av教育| 中文字幕人成人乱码亚洲影| 天堂网av新在线| 亚洲九九香蕉| 中文亚洲av片在线观看爽| 国产成人一区二区三区免费视频网站| 少妇丰满av| 高潮久久久久久久久久久不卡| 淫秽高清视频在线观看| 国产毛片a区久久久久| 久久午夜亚洲精品久久| 亚洲欧美日韩无卡精品| 真实男女啪啪啪动态图| 亚洲无线在线观看| h日本视频在线播放| 午夜精品久久久久久毛片777| 欧美日本亚洲视频在线播放| 一本精品99久久精品77| 舔av片在线| 一二三四在线观看免费中文在| 国产美女午夜福利| 亚洲片人在线观看| 天堂av国产一区二区熟女人妻| 亚洲精品中文字幕一二三四区| 国产视频一区二区在线看| 在线免费观看不下载黄p国产 | 99精品在免费线老司机午夜| 国产高清有码在线观看视频| 欧美三级亚洲精品| 女人被狂操c到高潮| 一进一出好大好爽视频| 99久久99久久久精品蜜桃| 亚洲精华国产精华精| 日韩人妻高清精品专区| 日本三级黄在线观看| 床上黄色一级片| 欧美精品啪啪一区二区三区| 99在线视频只有这里精品首页| 亚洲国产精品sss在线观看| 两个人看的免费小视频| 制服丝袜大香蕉在线| 精品国产亚洲在线| 免费看日本二区| 天堂动漫精品| 亚洲自偷自拍图片 自拍| 亚洲最大成人中文| 最近在线观看免费完整版| 97人妻精品一区二区三区麻豆| 精品一区二区三区av网在线观看| 亚洲九九香蕉| 精品一区二区三区视频在线观看免费| 亚洲成a人片在线一区二区| 婷婷六月久久综合丁香| 精品久久久久久久人妻蜜臀av| 国产成人精品久久二区二区91| 每晚都被弄得嗷嗷叫到高潮| 亚洲欧美精品综合久久99| 99在线人妻在线中文字幕| 国产精品国产高清国产av| 在线观看免费午夜福利视频| 国产精品av视频在线免费观看| 成人一区二区视频在线观看| 啪啪无遮挡十八禁网站| 欧美一区二区精品小视频在线| 亚洲avbb在线观看| 久久精品综合一区二区三区| 国语自产精品视频在线第100页| 成人性生交大片免费视频hd| www日本黄色视频网| 午夜精品一区二区三区免费看| 欧美成人性av电影在线观看| 久久久久久九九精品二区国产| 日韩欧美在线二视频| 成年女人永久免费观看视频| 国产黄片美女视频| 美女 人体艺术 gogo| 男女做爰动态图高潮gif福利片| 国产1区2区3区精品| 麻豆av在线久日| 2021天堂中文幕一二区在线观| 国产在线精品亚洲第一网站| 狠狠狠狠99中文字幕| 免费在线观看影片大全网站| 国产成人啪精品午夜网站| 亚洲av电影在线进入| 一区二区三区激情视频| 99久国产av精品| 国产一区二区激情短视频| 亚洲专区国产一区二区| 日本黄色片子视频| 天堂√8在线中文| 亚洲色图av天堂| 热99re8久久精品国产| 宅男免费午夜| 国产精品乱码一区二三区的特点| 日韩高清综合在线| 国内少妇人妻偷人精品xxx网站 | 日本五十路高清| 亚洲男人的天堂狠狠| 亚洲av成人av| 国产午夜精品论理片| 精品久久久久久久末码| 麻豆久久精品国产亚洲av| 日本撒尿小便嘘嘘汇集6| 国产不卡一卡二| 高潮久久久久久久久久久不卡| 色播亚洲综合网| 精品人妻1区二区| 欧美一级a爱片免费观看看| 99久久无色码亚洲精品果冻| 精品人妻1区二区| 午夜精品一区二区三区免费看| 99久久无色码亚洲精品果冻| 免费看a级黄色片| 看免费av毛片| 欧美成人一区二区免费高清观看 | 国产 一区 欧美 日韩| 欧美成人性av电影在线观看| 97碰自拍视频| 99热这里只有是精品50| 亚洲黑人精品在线| 2021天堂中文幕一二区在线观| 精品国内亚洲2022精品成人| 久久亚洲真实| 午夜福利18| 天堂网av新在线| 精品久久久久久久末码| 亚洲国产精品999在线| 免费一级毛片在线播放高清视频| 十八禁人妻一区二区| 97人妻精品一区二区三区麻豆| 中文字幕高清在线视频| 一本一本综合久久|