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

    基于實測軌跡的高精度捷聯(lián)慣導(dǎo)模擬器

    2015-07-07 00:53:26嚴(yán)恭敏WANGJinling周馨怡
    導(dǎo)航定位學(xué)報 2015年4期
    關(guān)鍵詞:演算法慣導(dǎo)模擬器

    嚴(yán)恭敏,WANG Jinling,周馨怡

    (1.西北工業(yè)大學(xué)自動化學(xué)院,西安 710072; 2.School of Civil and Environment Engineering,University of New South Wales,Sydney,NSW 2052,Australia)

    基于實測軌跡的高精度捷聯(lián)慣導(dǎo)模擬器

    嚴(yán)恭敏1,WANG Jinling2,周馨怡1

    (1.西北工業(yè)大學(xué)自動化學(xué)院,西安 710072; 2.School of Civil and Environment Engineering,University of New South Wales,Sydney,NSW 2052,Australia)

    針對傳統(tǒng)的純數(shù)學(xué)模型SINS慣性器件模擬器仿真效果不夠逼真的問題,本文提出了一種基于已有的真實導(dǎo)航系統(tǒng)飛行記錄數(shù)據(jù)的慣性器件模擬器方法。該方法利用高精度載波相位差分GNSS與中等精度IMU融合的組合導(dǎo)航算法,結(jié)合Kalman濾波的部分反饋修正和三次樣條擬合方法,生成平滑的軌跡參數(shù),再通過推導(dǎo)SINS反演算法實現(xiàn)了慣性器件的模擬采樣仿真,反演算法中充分考慮了姿態(tài)圓錐誤差和速度劃船誤差的補償。通過機載飛行實驗數(shù)據(jù)驗證,結(jié)果表明所提慣性器件模擬器具有極高的精度性能和良好的頻率特性,滿足高精度SINS對數(shù)據(jù)源仿真精度和頻率復(fù)雜性的要求。

    GNSS;SINS模擬器;組合導(dǎo)航;實際飛行數(shù)據(jù)

    0 引言

    在慣性導(dǎo)航系統(tǒng)及其組合導(dǎo)航的仿真研究中,運載體的運動軌跡參數(shù)和慣性器件(陀螺和加速度計)的數(shù)據(jù)源仿真是項基礎(chǔ)性工作,特別是研究高精度的捷聯(lián)慣導(dǎo)系統(tǒng)(strapdown inertial navigation system,SINS)算法,更是需要高精度的數(shù)據(jù)模擬器作為支撐[1]。以機載SINS仿真為例,傳統(tǒng)的基于純數(shù)學(xué)模型的軌跡模擬器,只能給出由靜止、加速、爬升、巡航和盤旋轉(zhuǎn)彎等簡單飛行動作組合的載機運動軌跡,運動形式過于簡單,效果不夠逼真,若要進(jìn)一步考慮載機飛控特性和氣動環(huán)境的影響,物理和數(shù)學(xué)模型往往又過于復(fù)雜,并且理想模型與實際飛行情況之間或多或少總會存在一些差距[2]。

    本文基于已有的真實導(dǎo)航傳感器飛行記錄數(shù)據(jù),利用高精度載波相位差分全球衛(wèi)星導(dǎo)航系統(tǒng)(global navigation satellite system,GNSS)與中等精度慣性測量單元(inertial measurement unit, IMU)融合的組合導(dǎo)航算法,結(jié)合卡爾曼(Kalman)濾波的部分反饋修正和三次樣條函數(shù)擬合和插值方法,生成平滑的軌跡參數(shù),即姿態(tài)、速度和位置信息。借助傳統(tǒng)的高精度SINS數(shù)值更新算法[3-5],推導(dǎo)了由姿態(tài)、速度和位置信息反向求解慣性器件增量采樣輸出的方法,稱其為SINS反演算法,在反演算法中考慮了姿態(tài)圓錐誤差和速度劃船誤差效應(yīng)的影響,進(jìn)行了相應(yīng)誤差的補償。利用SINS反演算法實現(xiàn)了慣性器件的模擬采樣輸出。通過飛行實驗數(shù)據(jù)驗證和對比分析,慣性器件模擬采樣輸出具有很高的導(dǎo)航精度,其精度僅受限于數(shù)值算法的計算誤差;慣性器件的模擬采樣輸出與原始真實采樣輸出在低頻段的功率譜密度基本一致,說明了模擬采樣輸出能夠較好地反映出載機的飛行環(huán)境特性。

    1 傳統(tǒng)的捷聯(lián)慣導(dǎo)算法

    1.1 捷聯(lián)慣導(dǎo)系統(tǒng)微分方程

    首先,簡要說明文中涉及到的幾個坐標(biāo)系:慣性坐標(biāo)系表示為i系,地球坐標(biāo)系表示為e系,“東-北-天”導(dǎo)航坐標(biāo)系為n系,“右-前-上”載體坐標(biāo)系為b系。

    SINS導(dǎo)航算法包含一組微分方程,即姿態(tài)微分方程、速度微分方程和位置微分方程,可分別表示為[6]

    1.2 捷聯(lián)慣導(dǎo)系統(tǒng)數(shù)值更新算法

    SINS數(shù)值更新算法是一種航位推算算法,它根據(jù)上一時刻(tm-1時刻)的已知姿態(tài)、速度和位置導(dǎo)航信息,結(jié)合慣性器件在更新周期內(nèi)的采樣輸出,遞推計算當(dāng)前時刻(tm時刻)的導(dǎo)航信息。記Tm=tm-tm-1為SINS更新周期。文中假設(shè)陀螺采樣輸出為角增量并且加速度計輸出為速度增量,高精度的慣導(dǎo)系統(tǒng)往往都采取這種采樣方式。

    考慮到數(shù)值更新算法的精度和為了后續(xù)推導(dǎo)反演算法的方便,下面以“單子樣+前一周期”誤差補償方法給出SINS數(shù)值更新算法的主體框架[3-6]。

    1.2.1 姿態(tài)更新算法

    下面給出由等效旋轉(zhuǎn)矢量Φ計算變換矩陣C的公式

    1.2.2 速度更新算法

    1.2.3 位置更新算法

    2 GNSS/SINS組合導(dǎo)航

    2.1 GNSS/SINS組合Kalman濾波

    采用高精度載波相位GNSS與SINS構(gòu)成組合導(dǎo)航系統(tǒng),GNSS可達(dá)厘米級的定位精度并且載機的機動大、飛行速度高,因此必需考慮到GNSS 和SINS之間的安裝桿臂誤差δl以及兩者之間的時間不同步誤差δt影響。這里建立19維的GNSS/ SINS組合Kalman濾波器,狀態(tài)選取及濾波器構(gòu)造為

    式(18)~式(20)中各符號的含義以及詳細(xì)系統(tǒng)參數(shù)可參見文獻(xiàn)[6-7]。

    2.2 Kalman濾波的部分反饋修正原理

    為了保持慣導(dǎo)系統(tǒng)的導(dǎo)航誤差始終為小量,即保證慣導(dǎo)誤差方程的線性,提高濾波器模型精度,在組合導(dǎo)航過程中需將Kalman濾波的狀態(tài)估計值不斷反饋給慣導(dǎo)系統(tǒng),進(jìn)行慣導(dǎo)修正。傳統(tǒng)的反饋方法是:選定某些狀態(tài)參與反饋,將Kalman濾波器中的這些狀態(tài)的估計值一次性全部反饋給慣導(dǎo)系統(tǒng),修正慣導(dǎo)參數(shù),再將濾波器中所有參與反饋的狀態(tài)清零。但是,上述全反饋方法容易造成慣導(dǎo)導(dǎo)航輸出劇烈的鋸齒狀波動,從而影響后續(xù)SINS反演算法的效果,在慣性器件模擬器上產(chǎn)生很大的間歇性脈沖失真。為了避免全反饋的不良影響,下面提出部分反饋修正原理,以經(jīng)度誤差的狀態(tài)反饋修正為例,表示為

    雖然GNSS/SINS組合導(dǎo)航的量測更新周期一般遠(yuǎn)大于SINS更新周期,但如果選擇反饋修正周期與SINS更新周期Tm一致,并利用上述部分反饋修正方法,則可得到比較平滑的導(dǎo)航參數(shù)輸出,為慣性模擬器的反演算法提供更加逼真的參考軌跡。當(dāng)然,如欲進(jìn)一步提高參考軌跡參數(shù)的求解精度,還可采用雙向濾波或平滑等算法[8]。

    3 捷聯(lián)慣導(dǎo)反演算法

    反演算法可以看作是常規(guī)SINS數(shù)值更新算法的逆過程。在反演算法中,假設(shè)已知姿態(tài)角和位置序列,前者為角運動信息而后者為線運動信息,再通過改造常規(guī)SINS更新算法,反向求解陀螺角增量和加速度計速度增量,實現(xiàn)慣性器件的高精度仿真和模擬。SINS反演算法的關(guān)鍵技術(shù)敘述如下。

    3.1 姿態(tài)和位置的三次樣條函數(shù)擬合

    以位置信息中的經(jīng)度為例,對組合導(dǎo)航輸出的經(jīng)度序列…,λm-1,λm,λm+1,…進(jìn)行分段三次樣條函數(shù)擬合,假設(shè)在時間區(qū)間t∈[tm-1,tm]上的擬合三次多項式函數(shù)為

    式(24)中,am0,am1,am2,am3為樣條擬合系數(shù)。

    經(jīng)度與慣導(dǎo)東向速度之間的關(guān)系為

    一般區(qū)間[tm-1,tm]很小,將式(25)中的L(t) 和RNh(t)取成tm-1/2時刻的值。由于位置函數(shù)已知,這里L(fēng)m-1/2可通過樣條函數(shù)插值精確獲得,而不象常規(guī)SINS更新算法中的那樣需通過線性外推進(jìn)行估計。

    由式(25)移項,再將式(24)代入,即可求得東向速度函數(shù)

    同理,對緯度和高度序列進(jìn)行處理,可求得北向速度和天向速度,從而獲得任意時刻的位置函數(shù)p(t)和速度函數(shù)vn(t);若對歐拉角序列進(jìn)行同樣處理,可獲得歐拉角分段擬合函數(shù),不妨記為A(t)。值得注意的是,如果歐拉角表示姿態(tài)存在奇異點,則需做特殊處理,這里不再詳述。

    隨著我國教育事業(yè)的不斷發(fā)展,教學(xué)模式也越來越多樣化,在對人才的教育教學(xué)過程中教師應(yīng)該不僅僅是進(jìn)行課堂教學(xué),通過利用數(shù)字媒體技術(shù)來提高學(xué)生的學(xué)習(xí)經(jīng)驗和個人素質(zhì)素養(yǎng)。對于每個人一生的教育學(xué)習(xí)來說,專業(yè)核心能力的教育培養(yǎng)非常重要。加強人才培養(yǎng)過程中的數(shù)字媒體技術(shù)的運用,堅持社會主義教學(xué)事業(yè)的快速發(fā)展,推動科學(xué)的教育事業(yè)發(fā)展,做好對孩子德智體美的全面發(fā)展工作,做好人才的前期的培養(yǎng)工作,為中國未來的發(fā)展奠定有利的教育基礎(chǔ)。

    3.2 陀螺的角增量模擬

    假設(shè)SINS反演算法生成慣性器件的解算周期為Tk=tk-tk-1,它可以不同于Tm,一般情況下小于Tm,Tk越小則解算精度越高。

    3.3 加速度計的速度增量模擬

    這里同樣可假設(shè)初值Δv0=0。至此,根據(jù)式(29)和式(31)即可求解得加速度計的速度增量,完成慣性器件的模擬采樣仿真。

    4 實驗數(shù)據(jù)處理與分析

    現(xiàn)有一組機載導(dǎo)航實驗數(shù)據(jù),其SINS中陀螺的隨機常值漂移約為0.03°/h,加速度計的逐次啟動隨機常值偏值約為200μg(特別說明:g= 9.8 m/s2,下同),原始數(shù)據(jù)輸出頻率200 Hz;載波相位差分GPS的定位精度為5 cm,數(shù)據(jù)輸出頻率1 Hz。載機飛行軌跡包括靜止、滑跑起飛、轉(zhuǎn)彎、巡航和在作業(yè)區(qū)域?qū)嵤┳鳂I(yè)等階段,圖1給出了1 h飛行軌跡的水平投影圖,其中左下角“☆”為軌跡起始點。

    圖1 飛行軌跡

    對實驗數(shù)據(jù)進(jìn)行處理,大致過程如下:

    (1)采用二子樣姿態(tài)圓錐誤差和速度劃船誤差補償算法進(jìn)行SINS導(dǎo)航更新,同時進(jìn)行GPS/ SINS組合導(dǎo)航Kalman濾波,使用部分反饋修正方法修正慣導(dǎo)誤差,這些算法均基于PSINS工具箱實現(xiàn)[9];

    (2)對組合導(dǎo)航輸出的100 Hz姿態(tài)和位置數(shù)據(jù)作三次樣條函數(shù)擬合,獲得分段連續(xù)的軌跡函數(shù)A(t)、vn(t)和p(t);

    (3)對軌跡函數(shù)A(t)、vn(t)和p(t)作200 Hz等間隔插值,得到序列Ak和pk;

    (4)根據(jù)SINS反演算法生成慣性器件模擬采樣值Δθk和Δvk,比較IMU原始采樣數(shù)據(jù)與模擬采樣數(shù)據(jù)的功率譜密度(power spectrum density, PSD),參見圖2,圖中顯示的是x軸陀螺和x軸加速度計在1 000~1 100 s的結(jié)果,其它軸和時間段的效果類似,不再給出;

    (5)選用二子樣誤差補償算法,由模擬采樣Δθk和Δvk進(jìn)行純慣導(dǎo)SINS解算,但對慣導(dǎo)的高度通道進(jìn)行了阻尼,純慣導(dǎo)解算結(jié)果與(1)中的組合導(dǎo)航結(jié)果比較,誤差參見圖3。

    圖2 原始采樣數(shù)據(jù)與模擬采樣數(shù)據(jù)的功率譜密度對比

    圖3 基于模擬采樣數(shù)據(jù)的純慣導(dǎo)誤差

    從圖2可看出,原始采樣數(shù)據(jù)與模擬采樣數(shù)據(jù)的功率譜密度在低頻段吻合得很好,說明了慣性器件模擬器能夠模擬出載機的低頻運動特性。圖3表明,模擬采樣數(shù)據(jù)具有很高的導(dǎo)航精度,在1 h內(nèi)純慣導(dǎo)的姿態(tài)誤差為0.001″量級、位置誤差為0.1 m量級。

    5 結(jié)束語

    根據(jù)傳統(tǒng)的高精度SINS數(shù)值更新算法,它由慣性器件輸出求解姿態(tài)、速度和位置等導(dǎo)航信息,提出了SINS反演算法,即如果已知姿態(tài)和位置軌跡信息,反向模擬求解慣性器件采樣值,實現(xiàn)慣性器件模擬器的仿真。針對實際載機的飛行數(shù)據(jù),需要結(jié)合GNSS進(jìn)行組合導(dǎo)航和使用部分反饋修正技術(shù)處理,提高導(dǎo)航精度和獲得平滑的軌跡參數(shù),才能作為慣性器件模擬器的理想輸入。實驗數(shù)據(jù)處理結(jié)果表明,所提慣性模擬器具有在1 h內(nèi)純慣導(dǎo)位置誤差0.1 m量級的精度,完全滿足高精度慣導(dǎo)系統(tǒng)對數(shù)據(jù)源仿真精度的要求,如果增加模擬器求解頻率,還有望進(jìn)一步減少仿真計算誤差;慣性器件的模擬采樣數(shù)據(jù)頻譜與實際采樣數(shù)據(jù)頻譜在低頻段基本相同,如果有更高的高頻一致性要求,可通過增加實際數(shù)據(jù)的采集頻率或者采用單子樣算法進(jìn)行SINS解算來實現(xiàn)。

    當(dāng)然,在高精度的模擬器數(shù)據(jù)源之后,還可針對特定的任務(wù)需求,再加入隨機漂移誤差和不正交安裝誤差等建模誤差,進(jìn)行更深層次的慣性導(dǎo)航或組合導(dǎo)航仿真。

    [1] HUDDLE J,GRETHEL R.Radio astro inertial Doppler error simulator[EB/OL].[2014-06-12].http://www.industrycortex.com.

    [2] GPSoft.Sat Nav toolbox 3.0 for Matlab[EB/OL].[2014-06-12].http://gpsoftnav.com/products/satellite-navigation-satnav-toolbox-3-0/.

    [3] PAUL G S.Strapdown inertial navigation integration algorithm design part 1:attitude algorithms[J].Journal of Guidance,Control and Dynamics,1998,21(1):19-28.

    [4] PAUL G S.Strapdown inertial navigation integration algorithm design part 2:velocity and position algorithms[J].Journal of Guidance,Control and Dynamics,1998,21(2):208-221.

    [5] WU Yuanxin,WANG Jinling,HU Dewen.A new technique for INS/GNSSattitude and parameter estimation using online optimization[J].IEEE Trans.on Signal Processing,2014,62(10):2642-2655.

    [6] 嚴(yán)恭敏.車載定位定向系統(tǒng)研究[D].西安:西北工業(yè)大學(xué),2006:16-27,82-84.

    [7] 趙春蓮.航空攝影測量中高精度POS導(dǎo)航算法研究[D].西安:西北工業(yè)大學(xué),2014:60-64.

    [8] 郝萬亮,孫付平,崔建勇,等.雙向濾波平滑在GNSS/INS組合導(dǎo)航中的應(yīng)用[C]//中國衛(wèi)星導(dǎo)航系統(tǒng)管理辦公室學(xué)術(shù)交流中心.第四屆中國衛(wèi)星導(dǎo)航學(xué)術(shù)年會論文集.北京:中國衛(wèi)星導(dǎo)航系統(tǒng)管理辦公室,2013:1-4.

    [9] 嚴(yán)恭敏.高精度捷聯(lián)慣性導(dǎo)航系統(tǒng)Matlab工具箱[EB/OL].(2013-09-16)[2014-06-12].http://blog.sina.com.cn/s/ blog_40edfdc90101heg0.html.

    High-Precision Simulator for Strapdown Inertial Navigation Systems Based on Real Dynamics

    YAN Gongmin1,WANG Jinling2,ZHOU Xinyi1
    (1.School of Automation,Northwestern Polytechnical University,Xi’an 710072,China; 2.School of Civil and Environment Engineering,University of New South Wales,Sydney 2052,Australia)

    Traditional strapdown inertial navigation systems(SINS)simulator based on pure mathematical model cannot reflect the reality of complex dynamics of a carrier.In this paper,based on the actual flight data of an airborne navigation system,navigation algorithms of high-precision carrier phase differences global navigation satellite system(GNSS)and mediumprecision inertial measurement unit(IMU)are used,and combined with the partial feedback principle of Kalman filter and cubic spline interpolation method,to generate smooth trajectory parameters.Then,the simulation of inertial sensors is achieved with the inversed SINSalgorithm,and the new algorithm considers the impacts of attitude coning error and velocity sculling error compensation.Numerical verifications of airborne flight test have demonstrated that the inertial sensor simulator has high accuracy and good frequency characteristic,and can meet the high-accuracy SINS requirements for simulated inertial sensor data source and frequency complexity.

    GNSS;SINS simulator;Integrated navigation;Actual flight data

    V249.3

    A

    2095-4999(2015)-04-0027-05

    2014-10-16

    嚴(yán)恭敏(1977—),福建建甌人,男,博士,副教授,主要從事慣性導(dǎo)航與信息融合理論研究。

    注:本文于2015年5月在第六屆中國衛(wèi)星導(dǎo)航學(xué)術(shù)年會(CSNC2015)宣讀過。

    嚴(yán)恭敏,WANG Jinling,周馨怡.基于實測軌跡的高精度捷聯(lián)慣導(dǎo)模擬器[J].導(dǎo)航定位學(xué)報,2015,3(4):27-31,37.YAN Gongmin, WANG Jinling,ZHOU Xinyi.High-Precision Simulator for Strapdown Inertial Navigation Systems Based on Real Dynamics[J].Journal of Navigation and Positioning,2015,3(4):27-31,37.

    10.16547/j.cnki.10-1096.20150406

    猜你喜歡
    演算法慣導(dǎo)模擬器
    《四庫全書總目》子部天文演算法、術(shù)數(shù)類提要獻(xiàn)疑
    國學(xué)(2021年0期)2022-01-18 05:59:08
    了不起的安檢模擬器
    盲盒模擬器
    劃船模擬器
    單多普勒天氣雷達(dá)非對稱VAP風(fēng)場反演算法
    自適應(yīng)模糊多環(huán)控制在慣導(dǎo)平臺穩(wěn)定回路中的應(yīng)用
    無人機室內(nèi)視覺/慣導(dǎo)組合導(dǎo)航方法
    基于Bagging模型的慣導(dǎo)系統(tǒng)誤差抑制方法
    運動平臺下X波段雷達(dá)海面風(fēng)向反演算法
    基于多線程的慣導(dǎo)邏輯仿真器設(shè)計
    計算機工程(2015年4期)2015-07-05 08:28:57
    国产亚洲精品av在线| 国产精品久久电影中文字幕| 可以在线观看毛片的网站| 日韩精品中文字幕看吧| 午夜福利18| 美女黄网站色视频| 日本 欧美在线| 色综合欧美亚洲国产小说| 国产精品一及| 欧美一区二区精品小视频在线| 欧美日本视频| 不卡av一区二区三区| 婷婷精品国产亚洲av| 男人舔奶头视频| 国产午夜精品久久久久久| 国产亚洲av嫩草精品影院| 色av中文字幕| 成人国产综合亚洲| 成人午夜高清在线视频| 99re在线观看精品视频| 九九在线视频观看精品| av天堂在线播放| 日本黄大片高清| 亚洲乱码一区二区免费版| 国产又色又爽无遮挡免费看| 亚洲九九香蕉| 亚洲色图 男人天堂 中文字幕| 99久国产av精品| xxx96com| 免费av不卡在线播放| 伊人久久大香线蕉亚洲五| 欧美午夜高清在线| а√天堂www在线а√下载| 国产真实乱freesex| 特大巨黑吊av在线直播| 国产成+人综合+亚洲专区| 国产一级毛片七仙女欲春2| 欧美av亚洲av综合av国产av| 九九热线精品视视频播放| 最近最新免费中文字幕在线| 校园春色视频在线观看| 无限看片的www在线观看| 欧美三级亚洲精品| 国产免费av片在线观看野外av| 亚洲精品国产精品久久久不卡| 日本a在线网址| 一级作爱视频免费观看| 日韩大尺度精品在线看网址| www日本在线高清视频| 亚洲一区二区三区色噜噜| 国产真人三级小视频在线观看| 久久久久久久午夜电影| 热99re8久久精品国产| 中文字幕最新亚洲高清| 变态另类成人亚洲欧美熟女| 在线视频色国产色| 视频区欧美日本亚洲| 午夜精品在线福利| 男女那种视频在线观看| 色综合站精品国产| 91av网站免费观看| 欧美成人免费av一区二区三区| 欧美高清成人免费视频www| 两个人的视频大全免费| 久久久国产精品麻豆| 精品熟女少妇八av免费久了| netflix在线观看网站| 亚洲专区中文字幕在线| 淫秽高清视频在线观看| 国产亚洲欧美98| 欧美一级a爱片免费观看看| 亚洲人成网站在线播放欧美日韩| 国产精品免费一区二区三区在线| 好男人在线观看高清免费视频| 九色国产91popny在线| 九九在线视频观看精品| 国产欧美日韩精品亚洲av| 欧美xxxx黑人xx丫x性爽| 欧洲精品卡2卡3卡4卡5卡区| 国产亚洲av高清不卡| 精品国产三级普通话版| 69av精品久久久久久| 操出白浆在线播放| 国产亚洲精品av在线| 日韩免费av在线播放| 亚洲18禁久久av| 欧美乱妇无乱码| 九色国产91popny在线| 亚洲 欧美 日韩 在线 免费| 国产亚洲欧美在线一区二区| av片东京热男人的天堂| 国产精品国产高清国产av| 欧美日本亚洲视频在线播放| 欧美在线一区亚洲| xxx96com| 日韩人妻高清精品专区| 欧美性猛交╳xxx乱大交人| 久久久久性生活片| 又黄又爽又免费观看的视频| 香蕉av资源在线| 99riav亚洲国产免费| 制服人妻中文乱码| 久久久久久人人人人人| av在线天堂中文字幕| e午夜精品久久久久久久| 黑人操中国人逼视频| 亚洲真实伦在线观看| 欧美性猛交╳xxx乱大交人| 国产97色在线日韩免费| 一个人看的www免费观看视频| 久久这里只有精品中国| 亚洲成av人片免费观看| 亚洲av片天天在线观看| 国产又色又爽无遮挡免费看| 成年人黄色毛片网站| 人妻丰满熟妇av一区二区三区| 不卡一级毛片| 久久久国产精品麻豆| 亚洲成a人片在线一区二区| 欧美日韩国产亚洲二区| 国产免费av片在线观看野外av| 久久久久国内视频| 国产1区2区3区精品| www.自偷自拍.com| 成人午夜高清在线视频| 夜夜躁狠狠躁天天躁| av在线蜜桃| 国产一区二区在线av高清观看| 国产麻豆成人av免费视频| 亚洲午夜理论影院| 亚洲美女视频黄频| 一级作爱视频免费观看| 亚洲最大成人中文| 美女黄网站色视频| 免费在线观看成人毛片| 亚洲色图av天堂| 国产精品,欧美在线| 久久午夜亚洲精品久久| 国产亚洲精品一区二区www| 亚洲国产精品999在线| 99久久无色码亚洲精品果冻| 欧美大码av| 999精品在线视频| 亚洲中文av在线| 日本撒尿小便嘘嘘汇集6| 特级一级黄色大片| 丰满人妻一区二区三区视频av | 欧美黄色淫秽网站| 小说图片视频综合网站| 欧美一区二区国产精品久久精品| 成人三级做爰电影| 亚洲成人久久性| 亚洲av电影在线进入| aaaaa片日本免费| 两性午夜刺激爽爽歪歪视频在线观看| 久久久久免费精品人妻一区二区| 一级毛片女人18水好多| 手机成人av网站| 少妇人妻一区二区三区视频| 嫁个100分男人电影在线观看| 免费观看精品视频网站| 巨乳人妻的诱惑在线观看| 亚洲精华国产精华精| 淫妇啪啪啪对白视频| 亚洲 国产 在线| 一个人观看的视频www高清免费观看 | 美女高潮的动态| 欧美日韩精品网址| 国产精品一区二区三区四区久久| 老司机福利观看| 熟女电影av网| 久久久久久国产a免费观看| 欧美丝袜亚洲另类 | 禁无遮挡网站| 日韩国内少妇激情av| 免费av毛片视频| 一级作爱视频免费观看| 岛国在线观看网站| 国产一区二区激情短视频| 成年女人毛片免费观看观看9| 欧美日韩综合久久久久久 | 最近最新中文字幕大全电影3| 久久久久久九九精品二区国产| 国产伦人伦偷精品视频| 久久国产精品人妻蜜桃| 美女大奶头视频| 美女cb高潮喷水在线观看 | 亚洲国产精品成人综合色| 看免费av毛片| 18禁黄网站禁片免费观看直播| 久久久久九九精品影院| 两性午夜刺激爽爽歪歪视频在线观看| 久久久水蜜桃国产精品网| 欧美3d第一页| 成人av一区二区三区在线看| 免费搜索国产男女视频| 国产激情欧美一区二区| 国产真人三级小视频在线观看| 国产精华一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 在线免费观看的www视频| 99热6这里只有精品| 黄色日韩在线| 国产精品日韩av在线免费观看| 成人鲁丝片一二三区免费| 女同久久另类99精品国产91| 搡老妇女老女人老熟妇| 午夜精品一区二区三区免费看| 国产av在哪里看| 一边摸一边抽搐一进一小说| 搡老熟女国产l中国老女人| 国产欧美日韩一区二区三| 国产精品久久电影中文字幕| 一级黄色大片毛片| 搡老妇女老女人老熟妇| 一区二区三区国产精品乱码| 99re在线观看精品视频| 久久这里只有精品19| 亚洲色图av天堂| 久久精品综合一区二区三区| 久久人妻av系列| 午夜福利18| 18美女黄网站色大片免费观看| 中文字幕av在线有码专区| 亚洲美女视频黄频| 国产单亲对白刺激| 国产三级中文精品| 日本一二三区视频观看| 曰老女人黄片| 国产精品自产拍在线观看55亚洲| 久久久精品大字幕| 成人三级黄色视频| 国产成人一区二区三区免费视频网站| 国产精品亚洲一级av第二区| 一夜夜www| xxx96com| 久久精品夜夜夜夜夜久久蜜豆| 99精品久久久久人妻精品| 白带黄色成豆腐渣| 制服丝袜大香蕉在线| 久久久久久久久免费视频了| 天天躁狠狠躁夜夜躁狠狠躁| 国产成人一区二区三区免费视频网站| 麻豆久久精品国产亚洲av| 欧美激情在线99| 国产激情偷乱视频一区二区| 久久天堂一区二区三区四区| 精品久久久久久成人av| 日韩欧美精品v在线| 亚洲色图 男人天堂 中文字幕| 成人国产综合亚洲| 91av网站免费观看| 亚洲av片天天在线观看| 老汉色av国产亚洲站长工具| 最新美女视频免费是黄的| 2021天堂中文幕一二区在线观| 给我免费播放毛片高清在线观看| www日本在线高清视频| 女生性感内裤真人,穿戴方法视频| 久久伊人香网站| 国产男靠女视频免费网站| 男插女下体视频免费在线播放| 欧美乱色亚洲激情| 色综合站精品国产| 国产探花在线观看一区二区| 国产91精品成人一区二区三区| 亚洲第一欧美日韩一区二区三区| 又粗又爽又猛毛片免费看| www.www免费av| 在线十欧美十亚洲十日本专区| 国产一区二区在线av高清观看| 国产高清videossex| 国产精品久久久av美女十八| av片东京热男人的天堂| 日韩欧美国产一区二区入口| 一进一出抽搐动态| 成人av一区二区三区在线看| 亚洲欧美日韩高清专用| 动漫黄色视频在线观看| 久久久久国内视频| 久久国产精品影院| 久久天躁狠狠躁夜夜2o2o| 成人永久免费在线观看视频| 成在线人永久免费视频| 亚洲 欧美一区二区三区| 国产免费av片在线观看野外av| 免费看日本二区| 好男人电影高清在线观看| 麻豆av在线久日| 国产免费男女视频| 免费观看人在逋| 亚洲中文日韩欧美视频| 亚洲国产欧洲综合997久久,| 两个人看的免费小视频| 国产免费av片在线观看野外av| 两人在一起打扑克的视频| 老熟妇仑乱视频hdxx| 久久久久九九精品影院| 亚洲熟女毛片儿| 少妇裸体淫交视频免费看高清| 免费一级毛片在线播放高清视频| 天堂动漫精品| 国产高清videossex| 最新中文字幕久久久久 | 99视频精品全部免费 在线 | 人人妻,人人澡人人爽秒播| 这个男人来自地球电影免费观看| 欧美日韩综合久久久久久 | 国产精品99久久久久久久久| 香蕉丝袜av| 小蜜桃在线观看免费完整版高清| 欧美日韩中文字幕国产精品一区二区三区| 可以在线观看毛片的网站| 美女黄网站色视频| xxx96com| 一进一出好大好爽视频| 99久久成人亚洲精品观看| www日本黄色视频网| 97超级碰碰碰精品色视频在线观看| 日韩免费av在线播放| 美女cb高潮喷水在线观看 | 在线免费观看不下载黄p国产 | 亚洲第一欧美日韩一区二区三区| 丁香欧美五月| av片东京热男人的天堂| 午夜福利18| 久久久久性生活片| 国产野战对白在线观看| 国产单亲对白刺激| www日本黄色视频网| 久久久国产成人精品二区| 精品人妻1区二区| 国产真实乱freesex| 日本黄大片高清| 久久性视频一级片| 动漫黄色视频在线观看| 好男人在线观看高清免费视频| 日韩中文字幕欧美一区二区| 两性夫妻黄色片| 一区二区三区激情视频| 后天国语完整版免费观看| 精品国产亚洲在线| 欧美色视频一区免费| 999久久久精品免费观看国产| 久久久久九九精品影院| 欧美日韩黄片免| 两性夫妻黄色片| 亚洲电影在线观看av| 亚洲精品久久国产高清桃花| 九九在线视频观看精品| 亚洲成av人片免费观看| 成人av一区二区三区在线看| 狂野欧美白嫩少妇大欣赏| 国产爱豆传媒在线观看| 99热6这里只有精品| 男人舔奶头视频| 国产麻豆成人av免费视频| 大型黄色视频在线免费观看| 在线看三级毛片| 人人妻,人人澡人人爽秒播| 观看免费一级毛片| 日本 av在线| 亚洲精品456在线播放app | 午夜久久久久精精品| 黄色女人牲交| 欧美国产日韩亚洲一区| 国产精品,欧美在线| 制服人妻中文乱码| 亚洲精品美女久久久久99蜜臀| 真实男女啪啪啪动态图| 99精品欧美一区二区三区四区| 九色国产91popny在线| 亚洲,欧美精品.| 国产精品一区二区免费欧美| 十八禁网站免费在线| 亚洲人成网站高清观看| 亚洲精品国产精品久久久不卡| 国产av麻豆久久久久久久| 国产一区二区在线av高清观看| 亚洲欧美日韩高清在线视频| 男人舔女人下体高潮全视频| 美女 人体艺术 gogo| 婷婷精品国产亚洲av| 在线观看美女被高潮喷水网站 | 黄片小视频在线播放| 国产综合懂色| 国产私拍福利视频在线观看| 久久热在线av| 波多野结衣高清作品| 成人一区二区视频在线观看| 国产一区二区三区在线臀色熟女| 搞女人的毛片| 亚洲av熟女| 国产精品一区二区三区四区免费观看 | 国产亚洲精品久久久com| 精品不卡国产一区二区三区| 99riav亚洲国产免费| 国产精品亚洲一级av第二区| 麻豆国产av国片精品| 老熟妇仑乱视频hdxx| 99久久国产精品久久久| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲国产精品久久男人天堂| 亚洲,欧美精品.| 亚洲在线自拍视频| 久久久国产精品麻豆| 91九色精品人成在线观看| 国产黄色小视频在线观看| 亚洲最大成人中文| 亚洲精品久久国产高清桃花| 12—13女人毛片做爰片一| 1024手机看黄色片| av天堂在线播放| 国产单亲对白刺激| 久久婷婷人人爽人人干人人爱| 久久午夜综合久久蜜桃| 变态另类丝袜制服| netflix在线观看网站| 国产黄色小视频在线观看| 婷婷亚洲欧美| 最近最新免费中文字幕在线| 三级国产精品欧美在线观看 | 五月伊人婷婷丁香| www.熟女人妻精品国产| 亚洲精品美女久久av网站| 丁香六月欧美| 亚洲狠狠婷婷综合久久图片| 欧美日韩瑟瑟在线播放| 久久久久九九精品影院| 大型黄色视频在线免费观看| 久久草成人影院| 婷婷六月久久综合丁香| 亚洲天堂国产精品一区在线| 国产成人福利小说| 无限看片的www在线观看| 91久久精品国产一区二区成人 | 99国产极品粉嫩在线观看| 首页视频小说图片口味搜索| 99久久精品热视频| 精品国内亚洲2022精品成人| 俄罗斯特黄特色一大片| 亚洲欧美精品综合一区二区三区| 亚洲色图 男人天堂 中文字幕| 欧美+亚洲+日韩+国产| 日本黄色片子视频| 久久99热这里只有精品18| 久久久精品大字幕| 国产精品亚洲美女久久久| 麻豆国产97在线/欧美| 国产亚洲精品av在线| av视频在线观看入口| 国产毛片a区久久久久| 欧美成人一区二区免费高清观看 | 久久久久性生活片| 少妇丰满av| 欧美大码av| 国产三级黄色录像| 一本综合久久免费| 欧美成人一区二区免费高清观看 | 成年女人毛片免费观看观看9| 久久这里只有精品中国| 免费观看精品视频网站| 制服丝袜大香蕉在线| 黄色日韩在线| 天堂网av新在线| 琪琪午夜伦伦电影理论片6080| 91九色精品人成在线观看| 99国产精品一区二区蜜桃av| 国产99白浆流出| 免费在线观看成人毛片| 男女下面进入的视频免费午夜| 一本精品99久久精品77| 午夜福利免费观看在线| 白带黄色成豆腐渣| 久久久国产成人免费| 亚洲乱码一区二区免费版| 国产精品 欧美亚洲| 老汉色av国产亚洲站长工具| 欧美绝顶高潮抽搐喷水| 久久久久国内视频| 此物有八面人人有两片| 精品99又大又爽又粗少妇毛片 | 国产亚洲精品av在线| 亚洲黑人精品在线| 日韩欧美 国产精品| 欧美中文日本在线观看视频| 国产私拍福利视频在线观看| 国产亚洲精品综合一区在线观看| 一个人观看的视频www高清免费观看 | 亚洲国产欧洲综合997久久,| 亚洲专区字幕在线| 午夜免费激情av| 久久久久亚洲av毛片大全| 成人鲁丝片一二三区免费| 黄色女人牲交| 中文字幕高清在线视频| 亚洲aⅴ乱码一区二区在线播放| 国产亚洲精品久久久久久毛片| 久99久视频精品免费| 91av网一区二区| 两个人看的免费小视频| 偷拍熟女少妇极品色| 国产精品1区2区在线观看.| 久9热在线精品视频| 亚洲国产精品sss在线观看| 国产aⅴ精品一区二区三区波| 精品久久久久久久毛片微露脸| 欧美不卡视频在线免费观看| 1024手机看黄色片| 巨乳人妻的诱惑在线观看| 91老司机精品| 黄频高清免费视频| 国产一级毛片七仙女欲春2| 最新在线观看一区二区三区| 精品乱码久久久久久99久播| 看片在线看免费视频| 两个人视频免费观看高清| 老汉色∧v一级毛片| 久久草成人影院| 深夜精品福利| 99国产精品一区二区三区| 国产亚洲精品久久久com| 亚洲美女黄片视频| 黑人欧美特级aaaaaa片| 观看美女的网站| 亚洲国产高清在线一区二区三| 97超级碰碰碰精品色视频在线观看| 色综合亚洲欧美另类图片| 色综合婷婷激情| 亚洲精华国产精华精| 欧美+亚洲+日韩+国产| 日本三级黄在线观看| 亚洲av日韩精品久久久久久密| 国产黄片美女视频| 国产美女午夜福利| 久久99热这里只有精品18| 色综合婷婷激情| 人人妻人人澡欧美一区二区| 久久久久国产一级毛片高清牌| 97超视频在线观看视频| 国产成+人综合+亚洲专区| 国产精华一区二区三区| 亚洲国产精品999在线| 国内精品久久久久精免费| 狠狠狠狠99中文字幕| 国产高潮美女av| 特大巨黑吊av在线直播| 亚洲精品粉嫩美女一区| 亚洲av五月六月丁香网| 美女 人体艺术 gogo| 丰满的人妻完整版| netflix在线观看网站| 18美女黄网站色大片免费观看| 国产精品久久久久久亚洲av鲁大| 中文在线观看免费www的网站| 亚洲性夜色夜夜综合| 久久国产精品影院| 国产亚洲欧美98| 男女床上黄色一级片免费看| 欧美精品啪啪一区二区三区| 久99久视频精品免费| 久久婷婷人人爽人人干人人爱| 亚洲av成人一区二区三| 男插女下体视频免费在线播放| 少妇人妻一区二区三区视频| 国产探花在线观看一区二区| 欧美黑人巨大hd| 99国产极品粉嫩在线观看| 国产 一区 欧美 日韩| 偷拍熟女少妇极品色| 欧美日韩黄片免| 人妻久久中文字幕网| 久久久久久久午夜电影| 国产亚洲av高清不卡| 首页视频小说图片口味搜索| 久久精品影院6| 亚洲av五月六月丁香网| 老司机午夜十八禁免费视频| 宅男免费午夜| 两性午夜刺激爽爽歪歪视频在线观看| 麻豆成人午夜福利视频| 亚洲国产精品sss在线观看| 中文字幕久久专区| 免费搜索国产男女视频| 女人高潮潮喷娇喘18禁视频| 夜夜躁狠狠躁天天躁| 免费搜索国产男女视频| 精品人妻1区二区| 天堂网av新在线| 夜夜爽天天搞| 国产野战对白在线观看| 脱女人内裤的视频| 国产精品 国内视频| 人人妻人人澡欧美一区二区| 日韩大尺度精品在线看网址| 男女视频在线观看网站免费| 欧美日本亚洲视频在线播放| 欧美最黄视频在线播放免费| 欧美xxxx黑人xx丫x性爽| 天天躁狠狠躁夜夜躁狠狠躁| 脱女人内裤的视频| 国模一区二区三区四区视频 | 最近最新中文字幕大全电影3| 51午夜福利影视在线观看| 国产爱豆传媒在线观看| 精品久久久久久久久久久久久| 国内毛片毛片毛片毛片毛片| 一区二区三区高清视频在线| 国产人伦9x9x在线观看| 日日干狠狠操夜夜爽| www.自偷自拍.com| 亚洲人成网站高清观看| 免费搜索国产男女视频| 村上凉子中文字幕在线| 脱女人内裤的视频|