• <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| 婷婷精品国产亚洲av在线| 最新中文字幕久久久久| 日韩av在线大香蕉| 99热这里只有精品一区| 麻豆av噜噜一区二区三区| 一级毛片久久久久久久久女| 男女那种视频在线观看| 久久人妻av系列| 国产淫片久久久久久久久| 午夜福利在线观看吧| 18禁裸乳无遮挡免费网站照片| 国产片特级美女逼逼视频| 1024手机看黄色片| 夜夜夜夜夜久久久久| 女的被弄到高潮叫床怎么办| 特大巨黑吊av在线直播| 国产91av在线免费观看| 国产又黄又爽又无遮挡在线| 一卡2卡三卡四卡精品乱码亚洲| 欧美成人一区二区免费高清观看| 欧美激情在线99| 搡老妇女老女人老熟妇| 欧美性猛交黑人性爽| 国产精品人妻久久久影院| 搡女人真爽免费视频火全软件 | 久久99热这里只有精品18| 国产不卡一卡二| 欧美不卡视频在线免费观看| 99热这里只有精品一区| 校园春色视频在线观看| 午夜福利成人在线免费观看| 久久鲁丝午夜福利片| 寂寞人妻少妇视频99o| 看十八女毛片水多多多| 人人妻人人看人人澡| 一级毛片电影观看 | 免费看a级黄色片| 欧美日韩综合久久久久久| 成人精品一区二区免费| 91午夜精品亚洲一区二区三区| 久久精品国产亚洲av天美| 免费看日本二区| 99久国产av精品| 可以在线观看的亚洲视频| 欧美色视频一区免费| 一本久久中文字幕| 久久人人爽人人片av| 亚洲成a人片在线一区二区| 一个人观看的视频www高清免费观看| av免费在线看不卡| 亚洲电影在线观看av| 黄色欧美视频在线观看| 日本黄色片子视频| 禁无遮挡网站| 露出奶头的视频| 精华霜和精华液先用哪个| 亚洲人成网站在线播放欧美日韩| 精品久久久噜噜| 色视频www国产| 最近手机中文字幕大全| 狂野欧美激情性xxxx在线观看| 成人国产麻豆网| 日韩欧美 国产精品| 日韩成人伦理影院| 别揉我奶头 嗯啊视频| 国产人妻一区二区三区在| 天美传媒精品一区二区| 18禁黄网站禁片免费观看直播| 欧美日本视频| 午夜福利成人在线免费观看| 美女内射精品一级片tv| 色综合站精品国产| 国产午夜精品久久久久久一区二区三区 | 舔av片在线| 欧美日韩精品成人综合77777| 免费观看在线日韩| 少妇熟女aⅴ在线视频| 别揉我奶头 嗯啊视频| 一级黄片播放器| 国产成人a∨麻豆精品| 高清毛片免费观看视频网站| 大型黄色视频在线免费观看| 又黄又爽又免费观看的视频| 亚洲在线自拍视频| 国产淫片久久久久久久久| 色播亚洲综合网| 精品久久久久久久久av| 91av网一区二区| 12—13女人毛片做爰片一| 又黄又爽又免费观看的视频| av专区在线播放| 日本在线视频免费播放| 在线免费观看不下载黄p国产| 中文字幕熟女人妻在线| 久久久久国产网址| 男人和女人高潮做爰伦理| 日韩成人av中文字幕在线观看 | 国产成人一区二区在线| 一级毛片aaaaaa免费看小| 国产精品永久免费网站| 久久久久国产网址| 波多野结衣高清无吗| 99久久精品热视频| 此物有八面人人有两片| 男女下面进入的视频免费午夜| 91麻豆精品激情在线观看国产| 国产一区二区亚洲精品在线观看| 久久久精品欧美日韩精品| 国产精品免费一区二区三区在线| 草草在线视频免费看| 亚洲最大成人av| 精品久久久久久久末码| av女优亚洲男人天堂| 极品教师在线视频| 日本黄色片子视频| 国产久久久一区二区三区| 联通29元200g的流量卡| 国产色婷婷99| 亚洲av成人精品一区久久| 精品久久久久久久久久免费视频| 国产白丝娇喘喷水9色精品| 欧美又色又爽又黄视频| 菩萨蛮人人尽说江南好唐韦庄 | 久久久国产成人精品二区| 成人性生交大片免费视频hd| 午夜精品在线福利| 99久久成人亚洲精品观看| 小蜜桃在线观看免费完整版高清| 国产视频一区二区在线看| 精品久久久久久久久亚洲| 听说在线观看完整版免费高清| 亚洲最大成人手机在线| 99久国产av精品| 午夜免费男女啪啪视频观看 | 卡戴珊不雅视频在线播放| 亚洲精品影视一区二区三区av| 丰满乱子伦码专区| 91久久精品国产一区二区三区| 国产一区亚洲一区在线观看| 日本爱情动作片www.在线观看 | 国产精品综合久久久久久久免费| 一区福利在线观看| 久久综合国产亚洲精品| 简卡轻食公司| 亚洲图色成人| 99久久九九国产精品国产免费| 日韩欧美精品v在线| 久久久久久久午夜电影| 免费观看的影片在线观看| 蜜桃久久精品国产亚洲av| 小蜜桃在线观看免费完整版高清| av在线蜜桃| 国产一区二区在线观看日韩| 国产一区二区在线观看日韩| 偷拍熟女少妇极品色| 一本久久中文字幕| 国产亚洲精品久久久com| 秋霞在线观看毛片| 色视频www国产| 亚洲av不卡在线观看| 在线播放无遮挡| 又黄又爽又免费观看的视频| 中国美女看黄片| 国产高清激情床上av| 狠狠狠狠99中文字幕| 美女高潮的动态| 国产伦精品一区二区三区四那| 免费在线观看影片大全网站| 亚洲国产欧美人成| 亚洲婷婷狠狠爱综合网| 99久国产av精品| 国产午夜精品论理片| 99热全是精品| 日日摸夜夜添夜夜添小说| 又黄又爽又免费观看的视频| 国产三级中文精品| 别揉我奶头~嗯~啊~动态视频| 99久久中文字幕三级久久日本| 亚洲av熟女| 亚洲第一电影网av| 三级经典国产精品| www.色视频.com| 性色avwww在线观看| 久久久久久久久久久丰满| 国产免费一级a男人的天堂| 久久久精品欧美日韩精品| 国产欧美日韩一区二区精品| 久久久久精品国产欧美久久久| av福利片在线观看| 久久精品久久久久久噜噜老黄 | 欧美丝袜亚洲另类| 国产在线男女| 国产精品一区二区三区四区免费观看 | 蜜桃久久精品国产亚洲av| 免费人成在线观看视频色| 婷婷精品国产亚洲av| 少妇被粗大猛烈的视频| 天堂√8在线中文| 午夜福利视频1000在线观看| 亚洲va在线va天堂va国产| 欧美性猛交╳xxx乱大交人| 国产黄a三级三级三级人| avwww免费| 寂寞人妻少妇视频99o| 一本一本综合久久| 亚洲内射少妇av| 一级毛片我不卡| 国产久久久一区二区三区| 国产精品99久久久久久久久| 网址你懂的国产日韩在线| 99热精品在线国产| 亚洲av不卡在线观看| 精品久久久久久久久久久久久| 免费人成在线观看视频色| 国产一区二区在线观看日韩| 少妇人妻一区二区三区视频| 99国产极品粉嫩在线观看| 日本色播在线视频| 精品午夜福利在线看| 久久久久久久午夜电影| 日韩在线高清观看一区二区三区| 麻豆一二三区av精品| 国产精品人妻久久久影院| 久久久色成人| 久久久久精品国产欧美久久久| 熟女电影av网| 精品国产三级普通话版| a级毛片a级免费在线| 日韩av在线大香蕉| 99久久成人亚洲精品观看| 日韩国内少妇激情av| 色噜噜av男人的天堂激情| 亚洲国产精品成人综合色| 中出人妻视频一区二区| 久久精品国产鲁丝片午夜精品| 亚洲性久久影院| 国产精品乱码一区二三区的特点| 亚洲精品日韩在线中文字幕 | 女的被弄到高潮叫床怎么办| 在线观看午夜福利视频| 国产精品久久久久久久久免| 亚洲三级黄色毛片| 日本五十路高清| 亚州av有码| 国产精品精品国产色婷婷| 麻豆乱淫一区二区| 看黄色毛片网站| 97碰自拍视频| 听说在线观看完整版免费高清| 男人和女人高潮做爰伦理| 乱码一卡2卡4卡精品| 你懂的网址亚洲精品在线观看 | 一卡2卡三卡四卡精品乱码亚洲| 日韩一区二区视频免费看| 不卡一级毛片| 国产私拍福利视频在线观看| 午夜精品在线福利| 国产av一区在线观看免费| 国产亚洲精品av在线| 两个人的视频大全免费| 欧美精品国产亚洲| 午夜激情福利司机影院| 国产精品野战在线观看| 午夜福利在线观看免费完整高清在 | 伦精品一区二区三区| 男人舔女人下体高潮全视频| 亚洲精品日韩在线中文字幕 | 黄色日韩在线| 国产真实伦视频高清在线观看| 五月玫瑰六月丁香| 精品一区二区三区视频在线| 精品久久久噜噜| 久久精品人妻少妇| 少妇裸体淫交视频免费看高清| 日本黄色片子视频| a级一级毛片免费在线观看| 亚洲中文字幕一区二区三区有码在线看| 国产精品国产三级国产av玫瑰| 国产精品嫩草影院av在线观看| 欧美最黄视频在线播放免费| 国产精品一区二区三区四区免费观看 | 久久天躁狠狠躁夜夜2o2o| 欧美成人一区二区免费高清观看| 日本五十路高清| 国产一区二区亚洲精品在线观看| 一个人观看的视频www高清免费观看| 精品99又大又爽又粗少妇毛片| 亚洲经典国产精华液单| 婷婷色综合大香蕉| 日韩av在线大香蕉| 午夜久久久久精精品| av在线蜜桃| 噜噜噜噜噜久久久久久91| 久久草成人影院| 一个人看的www免费观看视频| 99久久中文字幕三级久久日本| 亚洲成人久久爱视频| 国产精品亚洲美女久久久| 自拍偷自拍亚洲精品老妇| 日韩欧美免费精品| 国产成人a∨麻豆精品| 少妇的逼水好多| 久久精品国产自在天天线| 久久韩国三级中文字幕| 尤物成人国产欧美一区二区三区| 国产熟女欧美一区二区| 深夜精品福利| 成人性生交大片免费视频hd| 国产在线男女| 亚洲成人中文字幕在线播放| 精品免费久久久久久久清纯| 女的被弄到高潮叫床怎么办| 国产精品国产高清国产av| 亚洲熟妇中文字幕五十中出| 国产精品一二三区在线看| 波多野结衣巨乳人妻| 国产真实伦视频高清在线观看| 一夜夜www| 免费黄网站久久成人精品| 精品一区二区三区视频在线观看免费| 亚洲18禁久久av| 狂野欧美白嫩少妇大欣赏| 美女被艹到高潮喷水动态| av国产免费在线观看| 99久久无色码亚洲精品果冻| 亚洲不卡免费看| 日韩欧美 国产精品| 精品久久国产蜜桃| 男人狂女人下面高潮的视频| 亚洲一区二区三区色噜噜| 女生性感内裤真人,穿戴方法视频| 日本黄色片子视频| 一级黄色大片毛片| 久久精品影院6| 国产日本99.免费观看| 久久鲁丝午夜福利片| 国产精品免费一区二区三区在线| 麻豆精品久久久久久蜜桃| 婷婷精品国产亚洲av| 香蕉av资源在线| 欧美色视频一区免费| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲精品成人久久久久久| 春色校园在线视频观看| 少妇的逼水好多| 一区福利在线观看| 国产精品一区www在线观看| 欧美激情国产日韩精品一区| 日韩欧美三级三区| 国产精品嫩草影院av在线观看| 日韩在线高清观看一区二区三区| 亚洲国产精品sss在线观看| 欧美xxxx黑人xx丫x性爽| 看片在线看免费视频| 欧美又色又爽又黄视频| 99久久无色码亚洲精品果冻| 国产在线男女| 五月伊人婷婷丁香| av女优亚洲男人天堂| 18+在线观看网站| 免费av毛片视频| 国产日本99.免费观看| 成人亚洲欧美一区二区av| 久久久精品欧美日韩精品| 精品人妻一区二区三区麻豆 | 国产蜜桃级精品一区二区三区| 欧美色视频一区免费| 99久久九九国产精品国产免费| 久久精品影院6| 成人二区视频| 国产精品久久久久久亚洲av鲁大| 国产高清视频在线观看网站| 日本三级黄在线观看| 欧美区成人在线视频| 男女边吃奶边做爰视频| 久久久久九九精品影院| 成人精品一区二区免费| 美女免费视频网站| 69av精品久久久久久| 精品一区二区三区av网在线观看| 国产精品无大码| 99riav亚洲国产免费| 大又大粗又爽又黄少妇毛片口| 一区二区三区高清视频在线| 岛国在线免费视频观看| 久久久久久久亚洲中文字幕| 日产精品乱码卡一卡2卡三| 99riav亚洲国产免费| 久久亚洲国产成人精品v| 男人舔女人下体高潮全视频| 在线免费十八禁| 成人综合一区亚洲| 久久国产乱子免费精品| 午夜福利视频1000在线观看| 午夜免费男女啪啪视频观看 | 免费在线观看成人毛片| 大香蕉久久网| 欧美性猛交黑人性爽| 真实男女啪啪啪动态图| 深夜a级毛片| 免费观看的影片在线观看| 精品少妇黑人巨大在线播放 | 少妇高潮的动态图| 日韩强制内射视频| 中文字幕人妻熟人妻熟丝袜美| av在线观看视频网站免费| 免费观看的影片在线观看| 国产精品不卡视频一区二区| 99久久精品热视频| 欧美色欧美亚洲另类二区| 亚洲电影在线观看av| 最近视频中文字幕2019在线8| 久久久久久久久久久丰满| 午夜福利在线在线| 国产精华一区二区三区| 国产白丝娇喘喷水9色精品| 日韩三级伦理在线观看| 久久久欧美国产精品| 五月伊人婷婷丁香| 亚洲人与动物交配视频| 三级男女做爰猛烈吃奶摸视频| 欧美一区二区亚洲| 亚洲欧美日韩高清在线视频| 亚洲美女视频黄频| 97超级碰碰碰精品色视频在线观看| 91久久精品电影网| 一区二区三区免费毛片| 两性午夜刺激爽爽歪歪视频在线观看| 国产国拍精品亚洲av在线观看| 精品久久久久久成人av| 久久欧美精品欧美久久欧美| 亚洲精品国产av成人精品 | 一个人免费在线观看电影| a级一级毛片免费在线观看| 久久婷婷人人爽人人干人人爱| 久久久久性生活片| 婷婷亚洲欧美| or卡值多少钱| 真实男女啪啪啪动态图| 国产不卡一卡二| 久久午夜亚洲精品久久| 伦精品一区二区三区| 女生性感内裤真人,穿戴方法视频| 国产精品永久免费网站| 搡女人真爽免费视频火全软件 | 热99re8久久精品国产| 身体一侧抽搐| 久久6这里有精品| 欧美日韩国产亚洲二区| 国产一区二区激情短视频| 观看免费一级毛片| 日韩精品青青久久久久久| 日韩高清综合在线| 91在线观看av| 久久精品影院6| 男女那种视频在线观看| 亚洲真实伦在线观看| 国产片特级美女逼逼视频| 人妻制服诱惑在线中文字幕| 精品久久久久久久久av| 午夜久久久久精精品| 精品久久久久久久久亚洲| 夜夜夜夜夜久久久久| 你懂的网址亚洲精品在线观看 | 国产av在哪里看| 我要搜黄色片| 亚洲精华国产精华液的使用体验 | 变态另类丝袜制服| 校园春色视频在线观看| 精品人妻偷拍中文字幕| 久久久久国内视频| 国产成人freesex在线 | 国产精品一区二区免费欧美| 无遮挡黄片免费观看| 最新在线观看一区二区三区| 亚洲国产高清在线一区二区三| 日日干狠狠操夜夜爽| 久久久久免费精品人妻一区二区| 九色成人免费人妻av| 嫩草影院精品99| 天堂av国产一区二区熟女人妻| 亚洲18禁久久av| 十八禁国产超污无遮挡网站| 九九久久精品国产亚洲av麻豆| 久久精品国产亚洲网站| 国产视频一区二区在线看| 久久韩国三级中文字幕| 亚洲,欧美,日韩| 狂野欧美激情性xxxx在线观看| 亚洲av成人精品一区久久| 久久精品久久久久久噜噜老黄 | 国产 一区精品| 久久久色成人| 国产熟女欧美一区二区| 日本爱情动作片www.在线观看 | 久久精品人妻少妇| а√天堂www在线а√下载| 精品久久久久久久人妻蜜臀av| 波多野结衣高清无吗| 精品久久久久久久末码| 国产午夜精品论理片| 亚洲国产欧美人成| 在线播放国产精品三级| 国产男靠女视频免费网站| 国产真实乱freesex| 99视频精品全部免费 在线| 午夜老司机福利剧场| 欧美人与善性xxx| 欧美一区二区国产精品久久精品| 精品久久久久久久久久久久久| 午夜a级毛片| 直男gayav资源| 九九久久精品国产亚洲av麻豆| 欧美成人精品欧美一级黄| 国产精品无大码| 最新在线观看一区二区三区| 亚洲国产精品sss在线观看| 特大巨黑吊av在线直播| 少妇人妻一区二区三区视频| 成人二区视频| 97人妻精品一区二区三区麻豆| 啦啦啦啦在线视频资源| 日本与韩国留学比较| 最近2019中文字幕mv第一页| 日韩人妻高清精品专区| 最近的中文字幕免费完整| 97碰自拍视频| 久久人人爽人人片av| 亚洲经典国产精华液单| 日日摸夜夜添夜夜爱| 看免费成人av毛片| 日本a在线网址| 免费av观看视频| 麻豆国产97在线/欧美| 国产伦精品一区二区三区四那| 欧美一区二区国产精品久久精品| 国产av麻豆久久久久久久| 久久久精品94久久精品| 成人美女网站在线观看视频| 亚洲av免费在线观看| 淫妇啪啪啪对白视频| 久久久久久久久大av| 国产精品一区二区免费欧美| 日韩欧美三级三区| 男女那种视频在线观看| 十八禁网站免费在线| 在线播放国产精品三级| 亚洲综合色惰| 精品久久久噜噜| 男女边吃奶边做爰视频| 国产在线男女| 免费看美女性在线毛片视频| 日本熟妇午夜| 日日干狠狠操夜夜爽| 伦理电影大哥的女人| 国产精品久久久久久av不卡| av天堂中文字幕网| 一边摸一边抽搐一进一小说| 变态另类成人亚洲欧美熟女| 精品欧美国产一区二区三| 国产极品精品免费视频能看的| 午夜福利视频1000在线观看| 日本黄大片高清| 亚洲国产色片| 男人舔奶头视频| 麻豆av噜噜一区二区三区| 国产午夜精品久久久久久一区二区三区 | 欧美xxxx黑人xx丫x性爽| 在线观看免费视频日本深夜| 最后的刺客免费高清国语| 国产视频一区二区在线看| 性色avwww在线观看| 国产极品精品免费视频能看的| 午夜福利视频1000在线观看| 99视频精品全部免费 在线| 午夜福利视频1000在线观看| 国内精品美女久久久久久| 亚洲四区av| 亚洲第一区二区三区不卡| 老司机影院成人| 亚洲第一区二区三区不卡| 老司机影院成人| 亚洲精品亚洲一区二区| aaaaa片日本免费| 国产午夜精品久久久久久一区二区三区 | 亚洲婷婷狠狠爱综合网| 亚洲av二区三区四区| 女人十人毛片免费观看3o分钟| 特级一级黄色大片| av视频在线观看入口| 三级经典国产精品| 日韩强制内射视频| 亚洲,欧美,日韩| 亚洲精品影视一区二区三区av| 国产熟女欧美一区二区| 国产片特级美女逼逼视频| 波多野结衣高清无吗| 成熟少妇高潮喷水视频| www.色视频.com| 校园春色视频在线观看| 久久精品国产清高在天天线| 午夜久久久久精精品| 欧美一区二区国产精品久久精品| 国产探花极品一区二区| 一边摸一边抽搐一进一小说| 色哟哟哟哟哟哟| 1024手机看黄色片| 嫩草影视91久久| 午夜福利成人在线免费观看| 一进一出好大好爽视频| 成人综合一区亚洲| 欧美日韩乱码在线| 免费av不卡在线播放|