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

    用塔康系統(tǒng)的區(qū)域?qū)Ш蕉ㄎ凰惴?/h1>
    2012-08-27 08:16:06李曉東趙修斌
    電光與控制 2012年12期
    關(guān)鍵詞:航路測(cè)距投影

    李曉東, 趙修斌, 王 嘯

    (空軍工程大學(xué)信息與導(dǎo)航學(xué)院,西安 710077)

    0 引言

    塔康系統(tǒng)(Tactical Air Navigation System,TACAN)是軍用標(biāo)準(zhǔn)近程無(wú)線電導(dǎo)航系統(tǒng),作用距離為400~500 km,采用多瓣技術(shù),極坐標(biāo)體制定位,能在一種設(shè)備、一個(gè)頻道上同時(shí)測(cè)向和測(cè)距。塔康測(cè)向原理與民航伏爾(VOR)信標(biāo)相似,測(cè)距原理與測(cè)距儀(DME)相同。出于優(yōu)化資源利用和軍民融合的考慮,有些地區(qū)將VOR全向信標(biāo)與塔康地面臺(tái)合建于一處,稱為伏塔克(VORTAC),實(shí)現(xiàn)軍民航共用塔康測(cè)距功能[1]。

    為滿足民航空中流量增加、經(jīng)濟(jì)高效和安全飛行的要求,國(guó)際民航組織(ICAO)提出了基于性能導(dǎo)航(PBN)方案,大力發(fā)展區(qū)域?qū)Ш?RNAV)。不同于早期的航路導(dǎo)航圍繞導(dǎo)航臺(tái)對(duì)飛機(jī)進(jìn)行極坐標(biāo)定位,區(qū)域?qū)Ш睫D(zhuǎn)而關(guān)注飛機(jī)與航路點(diǎn)及計(jì)劃航線的相對(duì)位置關(guān)系。另外,PBN方案中的必備導(dǎo)航性能(RNP)規(guī)范要求機(jī)載設(shè)備具備實(shí)際導(dǎo)航性能監(jiān)視能力[2-3]。

    如果軍用飛機(jī)能夠?qū)崿F(xiàn)區(qū)域?qū)Ш?,將大大提高飛行員的態(tài)勢(shì)感知能力和作戰(zhàn)效能,然而軍機(jī)一般沒(méi)有足夠的空間容納已發(fā)展成熟的民航設(shè)備,必須立足軍用裝備進(jìn)行改進(jìn)升級(jí)。出于戰(zhàn)略安全考慮,軍用飛機(jī)不能依賴美國(guó)軍方主導(dǎo)的GPS全球衛(wèi)星導(dǎo)航系統(tǒng),而塔康系統(tǒng)可以為區(qū)域?qū)Ш教峁┯辛Φ难b備支撐,即使在北斗二代全球?qū)Ш较到y(tǒng)建立之后,依然可以作為重要的備份導(dǎo)航源。

    正常模式下,機(jī)載塔康設(shè)備主要測(cè)量飛機(jī)與塔康地面信標(biāo)的斜距和信標(biāo)臺(tái)磁方位角進(jìn)行極坐標(biāo)定位,提供的導(dǎo)航信息相對(duì)較少,并且定位誤差隨著距離增大而顯著增加,因而遠(yuǎn)遠(yuǎn)無(wú)法滿足注重態(tài)勢(shì)感知能力和高精度安全飛行的區(qū)域?qū)Ш叫阅芤?。本文針?duì)此問(wèn)題,提出塔康系統(tǒng)的3個(gè)改進(jìn)點(diǎn)以適應(yīng)區(qū)域?qū)Ш叫枰?1)求解飛機(jī)與航路點(diǎn)的平面位置關(guān)系,并結(jié)合大氣高程數(shù)據(jù)進(jìn)行相應(yīng)的坐標(biāo)變換,得到飛機(jī)在地球坐標(biāo)系下的三維位置估計(jì);2)通過(guò)ρ-ρ測(cè)距平面定位提高定位精度,滿足精確導(dǎo)航要求;3)位置估計(jì)與定位誤差算法一體化,實(shí)現(xiàn)機(jī)載實(shí)時(shí)導(dǎo)航性能監(jiān)測(cè)。最后通過(guò)數(shù)值仿真驗(yàn)證對(duì)其進(jìn)行應(yīng)用分析。

    1 塔康系統(tǒng)在地球坐標(biāo)系中定位

    不同于傳統(tǒng)的塔康極坐標(biāo)定位模式,本文以目標(biāo)航路點(diǎn)地面投影為原點(diǎn),建立東北天坐標(biāo)系及飛機(jī)位置量測(cè)方程,利用牛頓迭代法進(jìn)行平面位置估計(jì),結(jié)合大氣數(shù)據(jù)中的高程信息進(jìn)行坐標(biāo)系變換,得到飛機(jī)在地球坐標(biāo)系的三維位置估計(jì)。

    1.1 以航路點(diǎn)為中心的東北天(ENU)坐標(biāo)系下的位置解

    1.1.1 單臺(tái)定位:ρ-θ方式

    以航路點(diǎn)W在水平面投影W'為坐標(biāo)原點(diǎn),建立東北天坐標(biāo)系,假設(shè)飛機(jī)在過(guò)東向和北向平面上的投影為P',塔康信標(biāo)臺(tái)投影點(diǎn)為T',東北平面坐標(biāo)系內(nèi)三者的位置關(guān)系如圖1所示。

    圖1 東北平面坐標(biāo)系Fig.1 EN reference frame

    將P'與T'間距離真值記作r,用信標(biāo)臺(tái)磁方位角近似真方位角,記作θ。r和θ是飛機(jī)位置坐標(biāo)e、n的函數(shù),設(shè)r和θ的測(cè)量值分別為d、θc,不考慮測(cè)量誤差的情況下,量測(cè)方程為

    其中,(et,nt)表示塔康信標(biāo)臺(tái)在EN坐標(biāo)系中的位置,通過(guò)飛行管理系統(tǒng)(FMS)和導(dǎo)航數(shù)據(jù)庫(kù)可以預(yù)先解算和存儲(chǔ)[4-5]。

    由于牛頓迭代法具有計(jì)算量小、收斂速度快的優(yōu)點(diǎn),故本文采用該方法求解方程,基本步驟[6]如下所述。

    1) 設(shè)飛機(jī)的初始位置為(ek-1,nk-1),其與真值的誤差分別為 Δe、Δn,即 e=ek-1+ Δe,n=nk-1+ Δn。其中,k=1,2,3,…,表示迭代次數(shù)。

    2) 將測(cè)量方程在(ek-1,nk-1)點(diǎn)用泰勒級(jí)數(shù)展開(kāi),忽略二次及以上項(xiàng),則有

    由于式(2)左邊為Δe、Δn的線性方程組,故將其變形可得

    其中飛機(jī)坐標(biāo)e、n的二元函數(shù)矩陣G為

    G滿秩時(shí)則有

    G只與信標(biāo)臺(tái)和飛機(jī)的相對(duì)位置有關(guān),稱為幾何矩陣(也稱為雅可比矩陣)。

    4) 更新估計(jì)值 ek、nk,并重復(fù) 2)、3)步驟;ek、nk的更新值分別為

    1.1.2 雙臺(tái)定位:ρ-ρ方式

    當(dāng)飛機(jī)航行范圍內(nèi)有兩個(gè)或兩個(gè)以上的塔康信標(biāo)可用,理論上可以進(jìn)行ρ-ρ測(cè)距定位(塔康測(cè)距原理與民航大量使用的測(cè)距儀相同,因而可以進(jìn)一步考慮利用豐富的民航資源)。

    同樣在航路點(diǎn)投影為原點(diǎn)的東北天坐標(biāo)系中,兩個(gè)塔康信標(biāo)臺(tái)在東向和北向平面上投影為T'1(et1,nt1)、T'2(et2,nt2),飛機(jī)在平面內(nèi)投影 P'到信標(biāo)臺(tái)投影距離真值分別為r1、r2,設(shè)測(cè)量值分別為 d1、d2,暫時(shí)不考慮誤差項(xiàng),建立量測(cè)方程為

    1.2 飛機(jī)在地心地固坐標(biāo)系ECEF下的位置解

    在通過(guò)塔康系統(tǒng)進(jìn)行平面定位后,結(jié)合高程信息(如大氣數(shù)據(jù)),確定飛機(jī)在東北天坐標(biāo)系中的天向坐標(biāo)u,此時(shí)通過(guò)導(dǎo)航數(shù)據(jù)庫(kù)獲取航路點(diǎn)投影W'在大地坐標(biāo)系下的經(jīng)度λW,緯度φW信息,以及在地心地固參考系下的坐標(biāo)(xW,yW,zW)T,則飛機(jī)在地心地固坐標(biāo)系下的位置為

    其中,坐標(biāo)變換矩陣A為

    代入W'處的經(jīng)緯度λW、φW即可求解。需要指出的是,從飛機(jī)相對(duì)塔康信標(biāo)的斜距和信標(biāo)臺(tái)方位角信息轉(zhuǎn)換到飛機(jī)的地心地固直角坐標(biāo),對(duì)于區(qū)域?qū)Ш絹?lái)說(shuō)具有重要意義:它提供了更多的導(dǎo)航定位信息,并且是計(jì)算相對(duì)于航線和目標(biāo)航路點(diǎn)的偏航誤差、理想航跡角、航程和斜距等區(qū)域?qū)Ш絽?shù)的基礎(chǔ),此外,多導(dǎo)航傳感器信息融合時(shí),這一轉(zhuǎn)換也是前端進(jìn)行初始數(shù)據(jù)配準(zhǔn)的重要步驟。

    另外,還可以根據(jù)需要,通過(guò)坐標(biāo)轉(zhuǎn)換公式得出飛機(jī)直觀的經(jīng)緯高位置信息[7]。

    2 實(shí)際導(dǎo)航性能評(píng)估

    性能基導(dǎo)航PBN的RNP規(guī)范要求飛機(jī)可以實(shí)時(shí)評(píng)估導(dǎo)航位置估計(jì)誤差,也稱為位置不確定度(EPU)估計(jì)。本節(jié)主要對(duì)第1節(jié)所提出算法的定位誤差進(jìn)行分析。

    2.1 投影誤差

    假設(shè)誤差無(wú)偏,可推導(dǎo)得到誤差傳遞公式為

    上式是平面定位誤差分析研究的基礎(chǔ)。不難看出,在保持巡航高度h不變的情況下,飛機(jī)與信標(biāo)臺(tái)距離越遠(yuǎn),塔康測(cè)距誤差與水平投影誤差越接近;反之,飛機(jī)與信標(biāo)臺(tái)距離較近時(shí)誤差變化波動(dòng)較大。

    2.2 平面定位誤差

    實(shí)際導(dǎo)航性能評(píng)估實(shí)質(zhì)上就是機(jī)載設(shè)備實(shí)時(shí)計(jì)算定位誤差的過(guò)程[8-9],在第1節(jié)利用牛頓迭代方法解非線性方程組時(shí),得到幾何矩陣G,位置解算過(guò)程結(jié)束的同時(shí)就可以利用它進(jìn)入定位誤差的計(jì)算。

    2.2.1 ρ-θ方式定位誤差分析

    當(dāng)保留塔康測(cè)向誤差和測(cè)距誤差項(xiàng)時(shí),相應(yīng)地引起塔康與飛機(jī)水平距離r的誤差為εr,塔康信標(biāo)方位角θ的誤差為εθ,則有

    假設(shè)定位誤差均很小,對(duì)方程組線性化的影響可以忽略不計(jì),并且在進(jìn)行第k+1次迭代時(shí)收斂,飛機(jī)的位置坐標(biāo)此時(shí)變?yōu)橐阎浚?e=ek+Δe,n=nk+Δn。

    上式表明了測(cè)量誤差和定位誤差的關(guān)系。假設(shè)各個(gè)測(cè)量誤差εr、εθ均呈相同的正態(tài)分布,均值為0,方差為σ2,且各測(cè)量誤差互不相關(guān),則有

    進(jìn)一步可得標(biāo)準(zhǔn)定位誤差為

    2.2.2 ρ-ρ方式定位誤差分析

    設(shè)測(cè)量誤差 εd1、εd2,與 ρ-θ方式類似,可以得到

    式(12)表明了定位誤差量和測(cè)量誤差量的關(guān)系。不難看出,飛機(jī)與信標(biāo)臺(tái)的相對(duì)位置是影響定位精度的重要因素。θ1、θ2表示兩個(gè)信標(biāo)臺(tái)方位角估值,定位解算前可以用機(jī)載塔康系統(tǒng)的輸出量測(cè)值近似,從而很方便地估計(jì)幾何矩陣G以及標(biāo)準(zhǔn)定位誤差。

    假設(shè)測(cè)距誤差εd1、εd2均呈正態(tài)分布,均值為0,且不相關(guān),計(jì)算定位誤差協(xié)方差陣為

    標(biāo)準(zhǔn)定位誤差為

    進(jìn)一步假設(shè)測(cè)量誤差量εd1、εd2的方差都為σ2時(shí),則有

    D定義為權(quán)系數(shù)矩陣,利用三角函數(shù)關(guān)系合并化簡(jiǎn),得權(quán)系數(shù)矩陣的跡為

    其中,γ表示飛機(jī)與兩個(gè)信標(biāo)臺(tái)連線投影水平面的夾角,權(quán)系數(shù)的跡也可稱為幾何誤差因子GDOP。標(biāo)準(zhǔn)定位誤差可以表示為

    2.3 坐標(biāo)轉(zhuǎn)換誤差分析

    計(jì)算坐標(biāo)轉(zhuǎn)換誤差協(xié)方差矩陣為

    東北天坐標(biāo)系到地心地固直角坐標(biāo)的轉(zhuǎn)換矩陣A為單位正交矩陣,即A-1等于AT,容易證明

    即坐標(biāo)變換前后的誤差協(xié)方差矩陣的跡相等,定位誤差不變。

    3 仿真驗(yàn)證及應(yīng)用分析

    參考某型塔康設(shè)備技術(shù)參數(shù),設(shè)塔康的標(biāo)準(zhǔn)測(cè)角誤差2σθ=2,標(biāo)準(zhǔn)測(cè)距誤差2σs=200 m,飛機(jī)巡航高度10 km,塔康信標(biāo)臺(tái)在航路點(diǎn)東北天坐標(biāo)系中位置確定,利用Matlab進(jìn)行位置解算和誤差仿真。

    圖2表示在保持巡航高度不變的情況下,隨著飛機(jī)與信標(biāo)臺(tái)水平距離增大,水平距離誤差由較大值快速下降,并逐漸逼近斜距測(cè)量誤差;在60 km處,水平距離誤差與斜距誤差相差約為1%;需要說(shuō)明的是,圖中水平距離增大的過(guò)程實(shí)質(zhì)上體現(xiàn)的是飛機(jī)仰角減小的過(guò)程,因而在塔康信標(biāo)臺(tái)信號(hào)覆蓋范圍內(nèi),飛行高度與定位誤差沒(méi)有直接關(guān)系,不直接影響算法的適用性。

    圖2 σr/σs與r的關(guān)系Fig.2 Relation between σr/σsand r

    圖3 和圖4是采用極坐標(biāo)方式測(cè)角測(cè)距單臺(tái)定位誤差仿真分析。圖3表示飛機(jī)與信標(biāo)臺(tái)水平距離為100 km時(shí)(斜距100.5 km),信標(biāo)臺(tái)方位角由0變換至2π,定位誤差恒定為較大值3.496 km,結(jié)合第2節(jié)的理論分析可以推斷,當(dāng)飛機(jī)與信標(biāo)臺(tái)投影距離r一定時(shí),各個(gè)方向不同θ角處的標(biāo)準(zhǔn)定位誤差均相同。

    圖3 標(biāo)準(zhǔn)定位誤差2σ與信標(biāo)臺(tái)方位角θ的關(guān)系Fig.3 Relation between 2σ and θ

    圖4 ρ-θ模式標(biāo)準(zhǔn)定位誤差2σ與r的關(guān)系Fig.4 Relation between 2σ and r of ρ-θ mode

    圖4 顯示的是進(jìn)行單臺(tái)測(cè)角測(cè)距定位時(shí),飛機(jī)保持巡航高度不變的情況下,標(biāo)準(zhǔn)定位誤差2σ與水平距離r的關(guān)系。分析可得,飛機(jī)由塔康信標(biāo)臺(tái)水平距離1 km附近增加至7.6 km處,隨著飛機(jī)與信標(biāo)臺(tái)的仰角逐漸減小,水平投影距離r的誤差2σr迅速減小,導(dǎo)致定位誤差由1.468 km急劇減小至0.5491 km,并取得極小值(此時(shí)投影距離誤差并不是最小值,見(jiàn)圖1);隨著飛機(jī)進(jìn)一步遠(yuǎn)離塔康信標(biāo)臺(tái),水平投影距離測(cè)量誤差受高度影響可以忽略,可看作某一常數(shù)值,定位誤差與水平距離r呈近似線性遞增關(guān)系,變化率約為0.035,即每遠(yuǎn)離信標(biāo)臺(tái)1 km,標(biāo)準(zhǔn)定位誤差增加35 m。

    圖5和圖6是采用雙臺(tái)測(cè)距平面定位的標(biāo)準(zhǔn)誤差分析。圖5表示在飛機(jī)與兩個(gè)塔康信標(biāo)臺(tái)距離都超過(guò)60 km,忽略巡航高度對(duì)水平距離誤差的影響時(shí),標(biāo)準(zhǔn)定位誤差與測(cè)距位置線夾角γ的關(guān)系。容易看出,曲線在90°處取得極小值0.28 km,并以其為中心左右對(duì)稱。民航總局規(guī)定實(shí)施區(qū)域?qū)Ш綍r(shí),兩個(gè)DME臺(tái)站與飛機(jī)連線夾角應(yīng)在30°到150°之間[10],圖中做出了相應(yīng)的標(biāo)注。

    圖5 標(biāo)準(zhǔn)定位誤差2σ與位置線夾角γ的關(guān)系Fig.5 Relation between 2σ and γ

    圖6 表示的是測(cè)距位置線夾角為90°時(shí),標(biāo)準(zhǔn)定位誤差與水平距離的關(guān)系曲線,為控制變量,假設(shè)飛機(jī)與兩個(gè)塔康信標(biāo)臺(tái)距離相等??梢钥闯?,圖6曲線與圖2曲線在形狀上具有一致性,正是由于水平距離誤差在塔康信標(biāo)附近的非線性導(dǎo)致了定位誤差的跳變。

    民航總局規(guī)定使用DME進(jìn)行測(cè)距定位時(shí),從DME臺(tái)位置看,飛機(jī)與信標(biāo)臺(tái)仰角β應(yīng)低于40°[10],作為參考,圖中標(biāo)出了某型塔康裝備在仰角40°處的標(biāo)準(zhǔn)定位誤差值(約為0.3694 km)。

    圖6 ρ-ρ模式下標(biāo)準(zhǔn)定位誤差2σ與水平距離r的關(guān)系Fig.6 Relation between 2σ and r of ρ-ρ mode

    綜合以上仿真圖形,結(jié)合理論分析不難得出,塔康系統(tǒng)進(jìn)行單臺(tái)測(cè)角測(cè)距定位時(shí),定位誤差隨信標(biāo)距離的增大而增大,不適合平行航路的建立,也不能滿足高精度的RNP類型的需要(例如,RNP2要求導(dǎo)航系統(tǒng)誤差不能超過(guò)1 n mile,即1.852 km);雙臺(tái)測(cè)距定位在幾何構(gòu)型合理時(shí)有穩(wěn)定且較小的定位誤差,可以滿足大多數(shù)的航路RNP類型的精度要求,但是相對(duì)空間位置不理想時(shí),定位誤差急劇跳變。因此,考慮以測(cè)距定位為主,當(dāng)幾何誤差因子增大超出門限或者只有唯一的塔康信標(biāo)臺(tái)可用時(shí),可以切換到單臺(tái)測(cè)角測(cè)距定位。定位的同時(shí)在機(jī)載顯示終端實(shí)時(shí)給出實(shí)際導(dǎo)航性能參數(shù),以及相應(yīng)的告警。

    4 小結(jié)

    在區(qū)域?qū)Ш胶铰凤w行中,將機(jī)載塔康設(shè)備看作導(dǎo)航傳感器,輸出斜距和信標(biāo)臺(tái)方位角信息,在東北天和地心地固坐標(biāo)系下通過(guò)合理的重構(gòu)和解析可以提供更加豐富和有效的位置信息和區(qū)域?qū)Ш絽?shù)。理論分析和仿真表明,常規(guī)的測(cè)角測(cè)距定位方式因?yàn)殡S著與信標(biāo)臺(tái)距離的增大,定位誤差不斷增加,有些情況下不能滿足區(qū)域?qū)Ш狡叫泻铰芬约熬_導(dǎo)航的需要;可以考慮合理的技術(shù)改進(jìn),通過(guò)塔康組網(wǎng)實(shí)現(xiàn)雙臺(tái)甚至多臺(tái)測(cè)距平面定位,從而提高定位精度;此外,與 VOR/DME不同的是,由于塔康測(cè)角工作頻率較高,測(cè)角最大作用距離與測(cè)距是一致的,角度信息在機(jī)載實(shí)際導(dǎo)航性能監(jiān)測(cè)中可以得到有效利用。

    [1] 張忠興.無(wú)線電導(dǎo)航理論與系統(tǒng)[M].西安:陜西科學(xué)技術(shù)出版社,1998.

    [2] ICAO.Performance Based Navigation(PBN)manual[S].3rd ed.2008.

    [3] Federal Aviation Administration.Roadmap for Performance-Based Navigation[S].Washington DC,July 2006.

    [4] HERNDON A A.Analysis of advanced Flight Management Systems(FMS),F(xiàn)light Management Computer(FMC)field observations trials:Area Navigation(RNAV)holding patterns[C]//Digital Avionics Systems Conference(DASC),2011 IEEE/AIAA 30th,Oct.2011:1-17.

    [5] ATKINS S,CAPOZZI B.Relative position indicator for merging mixed RNAV and vectored arrival traffic[C]//Digital Avionics Systems Conference(DASC),2011 IEEE/AIAA 30th,Oct.2011:1-13.

    [6] 田孝華,周義建.無(wú)線電定位理論與技術(shù)[M].北京:國(guó)防工業(yè)出版社,2011.

    [7] 馬霞.TACAN數(shù)據(jù)在慣導(dǎo)參數(shù)校正中的坐標(biāo)變換[J].電訊技術(shù),2011,41(6):98-101.

    [8] CHAMLOU R.TIS-B:Calculation of navigation accuracy category for position and velocity parameters[C]//Proceeding of Digital Avionics Systems Conference,Salt Lake City,USA: [s.n.],2004,1:11-13.

    [9] 孫淑光,戴博,張鵬.機(jī)載組合導(dǎo)航系統(tǒng)實(shí)際導(dǎo)航性能計(jì)算方法[J].控制工程,2011(3):262-266.

    [10] 中國(guó)民用航空局飛行標(biāo)準(zhǔn)司.AC-91-FS在航路和終端區(qū)實(shí)施 RNAV1和 RNAV2的運(yùn)行指南[S].北京,2008,6.

    猜你喜歡
    航路測(cè)距投影
    解變分不等式的一種二次投影算法
    基于實(shí)時(shí)航路的PFD和ND的仿真研究
    基于最大相關(guān)熵的簇稀疏仿射投影算法
    類星體的精準(zhǔn)測(cè)距
    科學(xué)(2020年3期)2020-01-06 04:02:51
    找投影
    找投影
    淺談超聲波測(cè)距
    電子制作(2017年7期)2017-06-05 09:36:13
    應(yīng)召反潛時(shí)無(wú)人機(jī)監(jiān)聽(tīng)航路的規(guī)劃
    托勒密世界地圖與新航路的開(kāi)辟
    基于PSOC超聲測(cè)距系統(tǒng)設(shè)計(jì)

    黑人猛操日本美女一级片| 全区人妻精品视频| 日韩伦理黄色片| 激情五月婷婷亚洲| 99热6这里只有精品| 免费黄频网站在线观看国产| 波野结衣二区三区在线| 免费播放大片免费观看视频在线观看| 一二三四中文在线观看免费高清| 极品教师在线视频| 麻豆成人av视频| 久久久亚洲精品成人影院| 国产精品人妻久久久影院| 美女福利国产在线 | 97超碰精品成人国产| 欧美变态另类bdsm刘玥| 国产一区二区三区综合在线观看 | 街头女战士在线观看网站| 国产精品国产三级专区第一集| 国内少妇人妻偷人精品xxx网站| 亚洲欧美精品自产自拍| 亚洲精品视频女| 人妻夜夜爽99麻豆av| 免费观看av网站的网址| 久久人妻熟女aⅴ| 三级经典国产精品| 久久久久久久亚洲中文字幕| 激情 狠狠 欧美| 高清在线视频一区二区三区| 免费观看的影片在线观看| 伦理电影免费视频| 在线观看一区二区三区| 亚洲欧美精品自产自拍| 日本欧美视频一区| 中国国产av一级| 高清黄色对白视频在线免费看 | 麻豆乱淫一区二区| 青春草国产在线视频| 国产精品嫩草影院av在线观看| 插阴视频在线观看视频| 亚洲欧美中文字幕日韩二区| 色婷婷久久久亚洲欧美| 日韩亚洲欧美综合| 国产美女午夜福利| 国产视频内射| 亚洲成人一二三区av| 欧美日韩亚洲高清精品| 内地一区二区视频在线| 免费黄频网站在线观看国产| 国产免费一级a男人的天堂| 亚洲av.av天堂| 男男h啪啪无遮挡| 国产v大片淫在线免费观看| 国产伦在线观看视频一区| 这个男人来自地球电影免费观看 | 3wmmmm亚洲av在线观看| 欧美日韩视频高清一区二区三区二| 99热6这里只有精品| 久久久精品免费免费高清| 国产精品一区二区性色av| 黑人高潮一二区| 婷婷色综合www| 亚洲国产精品国产精品| 插逼视频在线观看| 亚洲av福利一区| 我要看黄色一级片免费的| 伊人久久国产一区二区| 久久人妻熟女aⅴ| 日韩大片免费观看网站| 日韩亚洲欧美综合| 另类亚洲欧美激情| 精品人妻视频免费看| 国产精品一区二区在线不卡| 国产无遮挡羞羞视频在线观看| 久久毛片免费看一区二区三区| 午夜福利网站1000一区二区三区| 成人特级av手机在线观看| 国产 精品1| 青春草视频在线免费观看| 成年美女黄网站色视频大全免费 | 深夜a级毛片| 99九九线精品视频在线观看视频| 一级av片app| 日韩一区二区视频免费看| 国产成人精品福利久久| 大陆偷拍与自拍| 亚洲最大成人中文| 久久久久久久久久久免费av| 精品少妇黑人巨大在线播放| 黄色日韩在线| 永久免费av网站大全| 只有这里有精品99| 一区二区av电影网| 啦啦啦在线观看免费高清www| 五月玫瑰六月丁香| 多毛熟女@视频| 嫩草影院入口| 精品亚洲成a人片在线观看 | 蜜桃亚洲精品一区二区三区| 亚洲成色77777| 国产成人免费观看mmmm| 精品一区二区三区视频在线| 久久精品熟女亚洲av麻豆精品| 制服丝袜香蕉在线| 成年免费大片在线观看| 久久久久久久国产电影| 亚洲人成网站高清观看| 免费看光身美女| 最近中文字幕2019免费版| 久久久久久九九精品二区国产| 中文欧美无线码| 国产精品一及| 久久久精品94久久精品| 美女cb高潮喷水在线观看| 青春草国产在线视频| av在线app专区| 亚洲人成网站在线观看播放| 一级毛片黄色毛片免费观看视频| 日韩不卡一区二区三区视频在线| 国产综合精华液| 狂野欧美白嫩少妇大欣赏| 国产又色又爽无遮挡免| 一区在线观看完整版| 国产成人免费无遮挡视频| 涩涩av久久男人的天堂| 成人毛片60女人毛片免费| 国产亚洲午夜精品一区二区久久| 一级av片app| 一个人看的www免费观看视频| 成人国产av品久久久| 免费久久久久久久精品成人欧美视频 | 51国产日韩欧美| 亚洲精品一二三| 3wmmmm亚洲av在线观看| 一级毛片我不卡| 日韩视频在线欧美| 久久国产乱子免费精品| 国产亚洲精品久久久com| 久久精品熟女亚洲av麻豆精品| 这个男人来自地球电影免费观看 | 在线观看国产h片| 欧美日韩亚洲高清精品| 国产高潮美女av| 91久久精品国产一区二区三区| 干丝袜人妻中文字幕| 大香蕉久久网| 日韩av免费高清视频| 麻豆乱淫一区二区| 高清午夜精品一区二区三区| 久久影院123| 久热久热在线精品观看| 大话2 男鬼变身卡| 国产av国产精品国产| 伊人久久国产一区二区| 欧美亚洲 丝袜 人妻 在线| 国产av国产精品国产| 国产亚洲av片在线观看秒播厂| 久久久久久九九精品二区国产| 日韩成人av中文字幕在线观看| 欧美成人午夜免费资源| 亚洲av欧美aⅴ国产| 亚洲国产精品一区三区| 精品亚洲乱码少妇综合久久| 久久久亚洲精品成人影院| 啦啦啦中文免费视频观看日本| 中文字幕人妻熟人妻熟丝袜美| 三级国产精品欧美在线观看| 女人十人毛片免费观看3o分钟| 精品亚洲成a人片在线观看 | 直男gayav资源| 日本黄大片高清| 91久久精品国产一区二区三区| 欧美极品一区二区三区四区| 亚洲欧美精品自产自拍| 久久6这里有精品| 国产亚洲av片在线观看秒播厂| 亚洲精品视频女| 久久国产亚洲av麻豆专区| 国产色爽女视频免费观看| 22中文网久久字幕| 能在线免费看毛片的网站| h视频一区二区三区| 欧美日韩在线观看h| 好男人视频免费观看在线| 日韩电影二区| 亚洲欧洲国产日韩| 国产精品一二三区在线看| 夫妻性生交免费视频一级片| 国产精品女同一区二区软件| 亚洲第一av免费看| 国产一区二区在线观看日韩| 国产黄色视频一区二区在线观看| .国产精品久久| 爱豆传媒免费全集在线观看| 欧美 日韩 精品 国产| 日本免费在线观看一区| 高清不卡的av网站| 久久久久久久亚洲中文字幕| 欧美日韩在线观看h| 亚洲综合精品二区| 看免费成人av毛片| 国产高清国产精品国产三级 | 一级a做视频免费观看| 久热久热在线精品观看| 亚洲av中文av极速乱| 久久国产亚洲av麻豆专区| 一级毛片aaaaaa免费看小| 欧美日韩在线观看h| 日本一二三区视频观看| 久久鲁丝午夜福利片| 五月玫瑰六月丁香| 亚洲av在线观看美女高潮| 亚洲av欧美aⅴ国产| 日日摸夜夜添夜夜爱| 成人一区二区视频在线观看| 亚洲伊人久久精品综合| 久热久热在线精品观看| 日日啪夜夜撸| 精品一区二区三区视频在线| 最近中文字幕2019免费版| 久久影院123| 丝袜喷水一区| 欧美bdsm另类| 热re99久久精品国产66热6| 成年美女黄网站色视频大全免费 | 免费av中文字幕在线| 中文字幕亚洲精品专区| 亚洲av国产av综合av卡| 中文字幕人妻熟人妻熟丝袜美| 日本免费在线观看一区| 麻豆国产97在线/欧美| 国产精品一二三区在线看| 免费看不卡的av| 国产一区二区在线观看日韩| 日本一二三区视频观看| 午夜免费男女啪啪视频观看| av国产免费在线观看| 久久99热这里只有精品18| 国产成人a∨麻豆精品| 亚洲国产精品国产精品| 三级经典国产精品| 丰满人妻一区二区三区视频av| 夜夜骑夜夜射夜夜干| 国产日韩欧美亚洲二区| 最近中文字幕2019免费版| 日日撸夜夜添| 欧美日韩国产mv在线观看视频 | 国产精品国产三级专区第一集| 亚洲精品乱久久久久久| 日本av手机在线免费观看| 久久久久精品性色| 成人综合一区亚洲| 欧美激情极品国产一区二区三区 | 舔av片在线| 精品人妻熟女av久视频| 亚洲丝袜综合中文字幕| 国内少妇人妻偷人精品xxx网站| 嫩草影院入口| 亚洲最大成人中文| 欧美一区二区亚洲| 亚洲精华国产精华液的使用体验| 精品午夜福利在线看| 99热这里只有精品一区| 国产午夜精品久久久久久一区二区三区| 天堂8中文在线网| 51国产日韩欧美| 国产成人免费观看mmmm| 91精品伊人久久大香线蕉| 免费黄网站久久成人精品| .国产精品久久| 国产男人的电影天堂91| 五月天丁香电影| 国产伦精品一区二区三区四那| 婷婷色麻豆天堂久久| 一级a做视频免费观看| 三级国产精品片| 九草在线视频观看| 久久久久精品久久久久真实原创| 天堂8中文在线网| 男人狂女人下面高潮的视频| 国产在线男女| 熟女电影av网| xxx大片免费视频| 性色av一级| 高清视频免费观看一区二区| 日产精品乱码卡一卡2卡三| 国产亚洲午夜精品一区二区久久| 最近手机中文字幕大全| 成人二区视频| 岛国毛片在线播放| 免费不卡的大黄色大毛片视频在线观看| 免费观看的影片在线观看| 蜜桃久久精品国产亚洲av| 内射极品少妇av片p| 日韩三级伦理在线观看| 最近2019中文字幕mv第一页| 国产成人午夜福利电影在线观看| 日本与韩国留学比较| 亚洲高清免费不卡视频| 国产深夜福利视频在线观看| 人体艺术视频欧美日本| 日韩三级伦理在线观看| 国产成人a∨麻豆精品| 亚洲精品一区蜜桃| 日日摸夜夜添夜夜添av毛片| 午夜福利在线观看免费完整高清在| 欧美日韩精品成人综合77777| 多毛熟女@视频| 亚洲av免费高清在线观看| 亚洲第一av免费看| 亚洲精品aⅴ在线观看| 视频区图区小说| 青春草亚洲视频在线观看| 韩国高清视频一区二区三区| 在线观看免费视频网站a站| 这个男人来自地球电影免费观看 | 免费观看无遮挡的男女| videos熟女内射| 国产乱人视频| 国精品久久久久久国模美| av黄色大香蕉| 一区二区三区精品91| 国模一区二区三区四区视频| 青春草视频在线免费观看| 亚洲欧洲国产日韩| 简卡轻食公司| 亚洲欧美日韩东京热| 日韩大片免费观看网站| 中文乱码字字幕精品一区二区三区| 色5月婷婷丁香| 久久99热这里只有精品18| 乱系列少妇在线播放| 久久鲁丝午夜福利片| 高清av免费在线| 成年av动漫网址| 久久99精品国语久久久| 久久久久网色| 国产男女超爽视频在线观看| 超碰av人人做人人爽久久| 一级爰片在线观看| 亚洲欧美成人综合另类久久久| 99热全是精品| 久久人妻熟女aⅴ| 日日摸夜夜添夜夜爱| 久久毛片免费看一区二区三区| 国产精品国产三级国产av玫瑰| 亚洲av不卡在线观看| 男的添女的下面高潮视频| 性色avwww在线观看| 欧美人与善性xxx| 一级毛片电影观看| 一级爰片在线观看| 国产精品av视频在线免费观看| 建设人人有责人人尽责人人享有的 | 老司机影院毛片| 在线观看美女被高潮喷水网站| 青春草视频在线免费观看| 3wmmmm亚洲av在线观看| a 毛片基地| 国产成人freesex在线| 色吧在线观看| 日韩一区二区三区影片| 99视频精品全部免费 在线| 国产精品爽爽va在线观看网站| 老熟女久久久| freevideosex欧美| 免费观看的影片在线观看| 免费看光身美女| 高清av免费在线| 国产深夜福利视频在线观看| 欧美高清性xxxxhd video| 免费大片黄手机在线观看| 一边亲一边摸免费视频| 国产爽快片一区二区三区| 精品人妻偷拍中文字幕| 欧美性感艳星| 中文字幕人妻熟人妻熟丝袜美| 婷婷色麻豆天堂久久| 国产成人免费无遮挡视频| 日韩欧美一区视频在线观看 | 免费不卡的大黄色大毛片视频在线观看| 99热6这里只有精品| 人体艺术视频欧美日本| 久久人人爽av亚洲精品天堂 | 久久久成人免费电影| 乱系列少妇在线播放| 成人免费观看视频高清| 一本色道久久久久久精品综合| 国产精品一二三区在线看| 你懂的网址亚洲精品在线观看| 舔av片在线| av在线老鸭窝| 国产精品成人在线| 精品久久久久久久末码| 蜜桃在线观看..| 少妇人妻久久综合中文| 一本一本综合久久| 国产高清三级在线| 国产爱豆传媒在线观看| 成人影院久久| 国产乱来视频区| 日日摸夜夜添夜夜添av毛片| 国产精品偷伦视频观看了| 免费观看的影片在线观看| av视频免费观看在线观看| 亚洲欧美日韩无卡精品| 毛片女人毛片| 少妇的逼水好多| 亚洲美女黄色视频免费看| 99久久精品国产国产毛片| 色综合色国产| 精品久久久久久久末码| 亚洲欧美一区二区三区黑人 | 久久婷婷青草| 丰满迷人的少妇在线观看| 国产在线视频一区二区| 国模一区二区三区四区视频| 国产免费一级a男人的天堂| 日韩 亚洲 欧美在线| 这个男人来自地球电影免费观看 | 亚洲国产色片| 亚洲精品自拍成人| 欧美丝袜亚洲另类| 26uuu在线亚洲综合色| 国产在线一区二区三区精| 免费不卡的大黄色大毛片视频在线观看| 欧美97在线视频| 欧美少妇被猛烈插入视频| 中国美白少妇内射xxxbb| 最近中文字幕高清免费大全6| h日本视频在线播放| 18禁裸乳无遮挡免费网站照片| 制服丝袜香蕉在线| 免费在线观看成人毛片| 少妇丰满av| 涩涩av久久男人的天堂| 色哟哟·www| 蜜桃在线观看..| 啦啦啦中文免费视频观看日本| 黄色日韩在线| av女优亚洲男人天堂| 久久久亚洲精品成人影院| 欧美3d第一页| 日本vs欧美在线观看视频 | 三级国产精品欧美在线观看| 人妻系列 视频| 亚洲欧美日韩无卡精品| 又爽又黄a免费视频| 国产一区亚洲一区在线观看| 80岁老熟妇乱子伦牲交| 在线观看美女被高潮喷水网站| 成人免费观看视频高清| 国产精品嫩草影院av在线观看| 亚洲一区二区三区欧美精品| 国产精品久久久久久久久免| 一级毛片久久久久久久久女| 国产久久久一区二区三区| 久久久午夜欧美精品| 午夜老司机福利剧场| 日本黄色日本黄色录像| 免费久久久久久久精品成人欧美视频 | 永久网站在线| 在线观看美女被高潮喷水网站| 日韩一区二区视频免费看| 十分钟在线观看高清视频www | 精品久久久精品久久久| 欧美+日韩+精品| 免费久久久久久久精品成人欧美视频 | 午夜老司机福利剧场| 亚洲经典国产精华液单| 一区二区三区四区激情视频| 国产欧美亚洲国产| 国产黄片视频在线免费观看| 99国产精品免费福利视频| 观看美女的网站| 另类亚洲欧美激情| 激情 狠狠 欧美| 尾随美女入室| 人人妻人人添人人爽欧美一区卜 | 夜夜看夜夜爽夜夜摸| 欧美性感艳星| 人人妻人人看人人澡| 舔av片在线| 另类亚洲欧美激情| 一级二级三级毛片免费看| 人人妻人人爽人人添夜夜欢视频 | 丰满人妻一区二区三区视频av| 老师上课跳d突然被开到最大视频| 我的女老师完整版在线观看| 极品少妇高潮喷水抽搐| 亚洲精品久久久久久婷婷小说| 久久精品国产亚洲网站| 国产乱人视频| 免费在线观看成人毛片| 亚洲一区二区三区欧美精品| 99热这里只有是精品在线观看| 久久99蜜桃精品久久| 久久久久性生活片| 亚洲精品久久久久久婷婷小说| 日本黄色片子视频| 国产精品成人在线| 天堂俺去俺来也www色官网| 国产午夜精品久久久久久一区二区三区| 91精品国产九色| 自拍偷自拍亚洲精品老妇| 午夜激情福利司机影院| 欧美日韩在线观看h| 国产精品久久久久久久久免| 男女下面进入的视频免费午夜| 久久久久久人妻| 最黄视频免费看| 欧美极品一区二区三区四区| 视频区图区小说| 亚洲av在线观看美女高潮| 高清日韩中文字幕在线| 精品亚洲成a人片在线观看 | 丝袜喷水一区| 在线观看三级黄色| 日韩电影二区| 夫妻午夜视频| 三级国产精品片| 一区在线观看完整版| 99热这里只有精品一区| 免费观看a级毛片全部| 九九爱精品视频在线观看| 日日撸夜夜添| 男人添女人高潮全过程视频| 午夜激情久久久久久久| 精品一品国产午夜福利视频| 有码 亚洲区| 尾随美女入室| 久久99热这里只频精品6学生| 国产欧美亚洲国产| 毛片一级片免费看久久久久| 三级经典国产精品| 六月丁香七月| 久久人妻熟女aⅴ| 51国产日韩欧美| 自拍偷自拍亚洲精品老妇| 国产精品久久久久久久久免| 日本wwww免费看| 乱码一卡2卡4卡精品| 日韩av免费高清视频| 久久这里有精品视频免费| 亚洲精品国产色婷婷电影| 国产精品99久久久久久久久| 91精品国产国语对白视频| 国产av码专区亚洲av| 亚洲欧美中文字幕日韩二区| 人妻夜夜爽99麻豆av| 亚洲精品成人av观看孕妇| 国产午夜精品久久久久久一区二区三区| 亚洲精品中文字幕在线视频 | 在线天堂最新版资源| 97热精品久久久久久| 欧美xxⅹ黑人| 中文字幕精品免费在线观看视频 | 一区在线观看完整版| 亚洲人成网站在线播| 亚洲av.av天堂| 免费人成在线观看视频色| 久久影院123| 黄色一级大片看看| 欧美日韩国产mv在线观看视频 | 国产片特级美女逼逼视频| 亚洲欧美一区二区三区黑人 | 精品久久久久久久久av| 国产成人午夜福利电影在线观看| 九色成人免费人妻av| 国产日韩欧美亚洲二区| 我的女老师完整版在线观看| 在线观看一区二区三区激情| 日韩av免费高清视频| 成人综合一区亚洲| 欧美另类一区| 中国美白少妇内射xxxbb| 国产高清有码在线观看视频| 国产精品一区二区在线观看99| 一本一本综合久久| av在线观看视频网站免费| 欧美 日韩 精品 国产| 少妇人妻久久综合中文| 秋霞伦理黄片| 精品少妇久久久久久888优播| 国产白丝娇喘喷水9色精品| 美女cb高潮喷水在线观看| 欧美精品一区二区大全| 婷婷色av中文字幕| 高清日韩中文字幕在线| 亚洲丝袜综合中文字幕| 偷拍熟女少妇极品色| 波野结衣二区三区在线| 美女cb高潮喷水在线观看| 精品人妻偷拍中文字幕| 日本黄色片子视频| 亚洲av成人精品一二三区| 少妇精品久久久久久久| 国产精品三级大全| 麻豆成人av视频| 天天躁日日操中文字幕| 日本-黄色视频高清免费观看| 日韩 亚洲 欧美在线| 成人黄色视频免费在线看| 国产在线一区二区三区精| 久久人人爽人人爽人人片va| 97超碰精品成人国产| 毛片一级片免费看久久久久| 高清午夜精品一区二区三区| 国产 精品1| 哪个播放器可以免费观看大片| 纯流量卡能插随身wifi吗| 自拍欧美九色日韩亚洲蝌蚪91 | 日韩成人伦理影院| 一级毛片我不卡| 这个男人来自地球电影免费观看 |