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

    基于自適應(yīng)EKF的BDS/GPS精密單點(diǎn)定位方法

    2016-09-07 03:18:48張勝宗
    關(guān)鍵詞:分量濾波觀測(cè)

    趙 琳, 張勝宗, 李 亮, 王 雪

    (哈爾濱工程大學(xué)自動(dòng)化學(xué)院, 黑龍江 哈爾濱 150001)

    ?

    基于自適應(yīng)EKF的BDS/GPS精密單點(diǎn)定位方法

    趙琳, 張勝宗, 李亮, 王雪

    (哈爾濱工程大學(xué)自動(dòng)化學(xué)院, 黑龍江 哈爾濱 150001)

    為使精密單點(diǎn)定位(precise point positioning,PPP)獲得更短的收斂時(shí)間和更高的定位精度,多個(gè)導(dǎo)航系統(tǒng)的集成(例如北斗衛(wèi)星導(dǎo)航系統(tǒng)(Beidou navigation satellite system,BDS)與全球定位系統(tǒng)(global positioning system,GPS)的組合)和更優(yōu)的定位方法是兩種可行選擇。針對(duì)傳統(tǒng)最小二乘(least square,LS)法解算孤立各歷元觀測(cè)量之間的關(guān)系以及擴(kuò)展卡爾曼濾波(extended Kalman filter,EKF)解算先驗(yàn)信息不準(zhǔn)的問題,在PPP中,運(yùn)用自適應(yīng)擴(kuò)展卡爾曼濾波(adaptive extended Kalman filter,AEKF)對(duì)過程噪聲進(jìn)行調(diào)整,以達(dá)到對(duì)系統(tǒng)狀態(tài)的最優(yōu)估計(jì)。文章通過實(shí)測(cè)數(shù)據(jù)對(duì)算法進(jìn)行了分析和驗(yàn)證,測(cè)試結(jié)果表明,與傳統(tǒng)的EKF算法相比,基于AEKF算法的PPP收斂速度可提高9 min,定位精度可提高33.7%。

    精密單點(diǎn)定位; 收斂時(shí)間; 定位精度; 自適應(yīng)擴(kuò)展卡爾曼濾波

    0 引 言

    精密單點(diǎn)定位(precise point positioning,PPP)技術(shù)指針對(duì)單臺(tái)接收機(jī)的載波相位觀測(cè)值,分別利用國(guó)際GNSS服務(wù)(international GNSS service,IGS)所發(fā)布的精密衛(wèi)星軌道與精密衛(wèi)星鐘差等產(chǎn)品,以及其他誤差修正模型對(duì)多種主要誤差源進(jìn)行改正來實(shí)現(xiàn)絕對(duì)定位的方法[1]。

    目前PPP技術(shù)應(yīng)用較為廣泛,基于偽距與載波相位觀測(cè)數(shù)據(jù),其靜態(tài)數(shù)據(jù)最終定位結(jié)果達(dá)到厘米級(jí)[2],相關(guān)技術(shù)已經(jīng)比較成熟。文獻(xiàn)[3]利用9個(gè)參考站模擬基于北斗的動(dòng)態(tài)定位,其結(jié)果顯示僅利用北斗系統(tǒng)的PPP,其定位精度可以達(dá)到分米級(jí)。但是,目前針對(duì)北斗衛(wèi)星導(dǎo)航系統(tǒng)(Beidou navigation satellite system,BDS)與全球定位系統(tǒng)(global positioning system,GPS)組合的PPP研究相對(duì)較少,組合BDS/GPS系統(tǒng)進(jìn)行PPP,不僅能夠增加觀測(cè)衛(wèi)星數(shù),提高量測(cè)冗余度,而且,較之單一系統(tǒng),組合系統(tǒng)定位的可靠性、連續(xù)性、精度都有了顯著的改善。隨著全球?qū)Ш叫l(wèi)星系統(tǒng)(global navigation saiellite system,GNSS)的現(xiàn)代化進(jìn)程,雙系統(tǒng)甚至多系統(tǒng)聯(lián)合定位已經(jīng)成為一種趨勢(shì)。通過組合可以充分發(fā)揮各系統(tǒng)的優(yōu)勢(shì),解決單系統(tǒng)在觀測(cè)條件不佳時(shí)可用衛(wèi)星數(shù)不足的問題,以提高定位的連續(xù)性[4]。

    最小二乘(least square,LS)和擴(kuò)展卡爾曼濾波(extended Kalman filter,EKF)是實(shí)現(xiàn)PPP兩種傳統(tǒng)的方法,其中傳統(tǒng)的LS解算方法,算法簡(jiǎn)單有效,但由于忽略了各歷元之間的關(guān)系,尤其是在高采樣率的條件下,短時(shí)間內(nèi)空間信號(hào)誤差存在較強(qiáng)的相關(guān)性[5]。由于忽略了歷元之間的相關(guān)信息,傳統(tǒng)LS方法需要更長(zhǎng)的時(shí)間完成定位收斂的過程??柭鼮V波是在線性無偏最小方差估計(jì)原理下推出的一種遞推濾波方法[6]。然而,傳統(tǒng)基于EKF的PPP方法大都采用經(jīng)驗(yàn)?zāi)P驮O(shè)置濾波先驗(yàn)狀態(tài)誤差協(xié)方差陣[7-8],實(shí)際在動(dòng)態(tài)定位的情況下,先驗(yàn)信息很難與真實(shí)的導(dǎo)航狀態(tài)相匹配,導(dǎo)致定位需要很長(zhǎng)的時(shí)間才能收斂。而在自適應(yīng)擴(kuò)展卡爾曼濾波(adaptive extended Kalman filter,AEKF)算法中,對(duì)過程噪聲和量測(cè)噪聲的估計(jì)可以優(yōu)化濾波算法的性能。系統(tǒng)過程噪聲和量測(cè)噪聲作為系統(tǒng)的主要先驗(yàn)信息,對(duì)卡爾曼濾波的估計(jì)性能和穩(wěn)定性有著重要的意義,對(duì)于兩者估計(jì)不準(zhǔn)確,可能會(huì)導(dǎo)致有偏估計(jì),甚至濾波發(fā)散[6]。來源于真實(shí)量測(cè)信息的新息數(shù)據(jù),對(duì)于卡爾曼濾波性能的評(píng)估而言具有重要的意義[9],通過新息方差對(duì)過程噪聲和量測(cè)噪聲實(shí)時(shí)估計(jì)和調(diào)整,從而更好適應(yīng)外界環(huán)境變化所引起的噪聲統(tǒng)計(jì)特性的不一致性,使得濾波器具有更佳的收斂性能[9]。使用AEKF能夠避免先驗(yàn)信息估計(jì)不準(zhǔn)對(duì)卡爾曼濾波性能的影響,可以進(jìn)一步提高PPP定位的收斂時(shí)間和精度。

    由于動(dòng)態(tài)導(dǎo)航環(huán)境對(duì)傳統(tǒng)卡爾曼模型的擾動(dòng),以及過程噪聲和量測(cè)噪聲的變化等因素的影響,制約了傳統(tǒng)卡爾曼濾波應(yīng)用的普適性[10]。通過采用極大似然準(zhǔn)則下的AEKF算法,對(duì)系統(tǒng)的過程噪聲進(jìn)行實(shí)時(shí)的估計(jì)和調(diào)整,使得濾波器可以對(duì)于系統(tǒng)噪聲統(tǒng)計(jì)特性的變化做出最優(yōu)化的適應(yīng)。針對(duì)BDS/GPS組合系統(tǒng),對(duì)這種自適應(yīng)濾波算法進(jìn)行了分析。實(shí)測(cè)計(jì)算結(jié)果證明該算法的有效性和可行性。

    1 組合精密單點(diǎn)原理與技術(shù)分析

    1.1觀測(cè)方程及誤差處理

    衛(wèi)星導(dǎo)航定位中基本的觀測(cè)方程為

    (1)

    (2)

    式中,φ、P分別為載波相位以及測(cè)碼偽距觀測(cè)量;ρ為用戶與衛(wèi)星之間真實(shí)的幾何距離;c為電磁波在真空中的傳播速度;dtr、dts分別為接收機(jī)鐘差以及衛(wèi)星鐘差;I為電離層延遲誤差;T為對(duì)流層延遲;dorb為衛(wèi)星軌道誤差;dtid為潮汐誤差;drel為相對(duì)論效應(yīng)誤差;N為整周模糊度;dant為天線相位誤差;dphw為相位纏繞誤差;εφ為載波相位觀測(cè)噪聲和多路徑誤差;εP為偽距觀測(cè)噪聲和多路徑誤差。

    實(shí)際應(yīng)用中,根據(jù)電離層介質(zhì)的彌散性質(zhì),采用雙頻觀測(cè)值消除電離層折射影響[11]。考慮組成L1、L2(B1、B2)頻率上的線性組合的觀測(cè)值,其組合系數(shù)α1、α2滿足:

    (3)

    式中,f1、f2分別為載波L1、L2(B1、B2)的頻率。原則上,α1、α2可以是任意一組滿足式(3)的值。但在實(shí)際應(yīng)用中,經(jīng)常選擇[12]

    (4)

    由此,可以利用雙頻觀測(cè)值消除電離層影響。通常使用無電離層組合(ionosphere-free,IF)觀測(cè)值,組合觀測(cè)量如下:

    (5)

    (6)

    式中,φIF、PIF分別為無電離層組合載波相位和偽距觀測(cè)值;NIF為組合整周模糊度;νφ、νP分別為無電離層相位和偽距組合觀測(cè)噪聲,此外

    (7)

    (8)

    式(7)和式(8)中所涉及誤差處理方法如表1所示。

    表1 PPP誤差處理方法

    在對(duì)各部分誤差進(jìn)行處理后,可以得到修正后的IF組合載波相位觀測(cè)值和偽距觀測(cè)值:

    (9)

    (10)

    1.2參數(shù)估計(jì)方法

    在上述對(duì)PPP各部分誤差進(jìn)行相應(yīng)的處理后,PPP的估計(jì)參數(shù)有接收機(jī)三維位置、三維速度(動(dòng)態(tài))、接收機(jī)鐘差、接收機(jī)頻漂、對(duì)流層濕分量延遲、組合模糊度。濾波處理中選取的狀態(tài)向量為

    (11)

    表2 狀態(tài)處理方法

    選取的量測(cè)量為每顆可視衛(wèi)星的IF組合載波相位與偽距觀測(cè)值:

    (12)

    然后對(duì)系統(tǒng)進(jìn)行建模,考慮到系統(tǒng)的非線性,需要先將濾波方程線性化,其中,量測(cè)方程非線性,系統(tǒng)模型表示為

    (13)

    (14)

    (15)

    其中

    (16)

    其中

    (17)

    (18)

    (19)

    則將非線性量測(cè)函數(shù)式(14)一階線性化后得到量測(cè)方程為

    (20)

    至此,可以得到線性化的系統(tǒng)模型式(13)和式(20),其中狀態(tài)轉(zhuǎn)移陣Φ和量測(cè)陣H可表示為

    (21)

    式中,I表示單位陣。靜態(tài)模式時(shí)A=03×3,B=03×3,C=[1,0;0,0];動(dòng)態(tài)模式時(shí)A=TsI3×3,B=I3×3,C=[1,Ts;0,1],Ts是狀態(tài)方程的差分步長(zhǎng),在這里也等于測(cè)量過程的采樣時(shí)間間隔。

    (22)

    式中,G3×k=(lk,mk,nk)(k=1,2,…,m)表示接收機(jī)到衛(wèi)星的觀測(cè)向量的方向余弦矩陣;I表示單位陣;E=[I1×m,01×m,I1×m,01×m];M1×m=(M1,M2,…,Mm)表示對(duì)流層濕分量映射函數(shù)矩陣。

    在得到由精密星歷文件計(jì)算的衛(wèi)星位置誤差方差,由精密鐘差文件計(jì)算的衛(wèi)星鐘差誤差方差,IF組合觀測(cè)噪聲方差以及對(duì)流層誤差方差之后,可以較為準(zhǔn)確地得到量測(cè)噪聲方差陣。此外,同時(shí)對(duì)過程噪聲和量測(cè)噪聲進(jìn)行調(diào)整的AEKF算法,容易導(dǎo)致濾波發(fā)散及濾波過程不穩(wěn)定[6]。因此,下面所推導(dǎo)的自適應(yīng)過程僅對(duì)過程噪聲進(jìn)行調(diào)整估計(jì)。

    本文運(yùn)用基于極大似然準(zhǔn)則的AEKF對(duì)系統(tǒng)狀態(tài)進(jìn)行估計(jì),極大似然準(zhǔn)則綜合新息以及新息協(xié)方差的變化,使其能夠適應(yīng)模型擾動(dòng)以及噪聲統(tǒng)計(jì)特性變化的影響[13]。已知系統(tǒng)的線性模型如式(13)和式(20)所示,基于參考文獻(xiàn)[6]可以計(jì)算自適應(yīng)過程下的過程噪聲為

    (23)

    式中,N為平滑窗口寬度;Δx為狀態(tài)更正量;P為估計(jì)誤差方差陣。

    2 數(shù)據(jù)處理與分析

    為充分驗(yàn)證所提出算法的性能,通過靜態(tài)和動(dòng)態(tài)PPP實(shí)驗(yàn)進(jìn)行了測(cè)試。需要指出的是,測(cè)試過程中,BDS和GPS精密星歷與精密鐘差數(shù)據(jù)分別由中科院測(cè)地所和IGS中心提供,更新時(shí)間分別為15 min和5 min。

    在對(duì)數(shù)據(jù)進(jìn)行處理時(shí),還應(yīng)考慮各個(gè)觀測(cè)值的權(quán)重大小,權(quán)重表示了各個(gè)觀測(cè)值對(duì)整個(gè)解的貢獻(xiàn)大小,正確匹配權(quán)值才能獲得最優(yōu)的結(jié)果,這其中涉及的以衛(wèi)星高度角定權(quán)的一種隨機(jī)模型形式為

    (24)

    (25)

    在濾波過程中,由式(23)可知在估計(jì)過程噪聲Q時(shí)涉及到Δx項(xiàng),從而直接反映系統(tǒng)的狀態(tài)變化,由此可見,準(zhǔn)確地估計(jì)過程噪聲將會(huì)使得濾波器更好地跟蹤動(dòng)態(tài)條件系統(tǒng)狀態(tài)的變化情況。所以在數(shù)據(jù)處理中,必須恰當(dāng)?shù)剡x擇平滑窗口的寬度N,為此,測(cè)試時(shí)對(duì)N進(jìn)行了不同的取值,取N從50開始,按步長(zhǎng)50增加,通過比較分析可知,當(dāng)窗口寬度N取200時(shí),狀態(tài)估計(jì)誤差精度達(dá)到最高。

    2.1靜態(tài)模式

    為測(cè)試靜態(tài)條件下BDS/GPS組合系統(tǒng)及單系統(tǒng)PPP的效果,采用澳大利亞CUT0站2014年6月17日的數(shù)據(jù)進(jìn)行定位解算,其接收機(jī)為Trimble Net R9,該觀測(cè)數(shù)據(jù)采樣率為30 s,截止高度角設(shè)為15°。在AEKF處理過程中,BDS和GPS載波相位觀測(cè)量的觀測(cè)噪聲標(biāo)準(zhǔn)差均設(shè)為0.003 m,偽距觀測(cè)量的觀測(cè)噪聲標(biāo)準(zhǔn)差分別設(shè)為0.6 m和0.3 m。

    本文分別對(duì)BDS單系統(tǒng),GPS單系統(tǒng),BDS/GPS組合系統(tǒng)3種情況進(jìn)行AEKF數(shù)據(jù)處理。以XYZ三方向的位置誤差(分別用errorX、errorY、errorZ表示)均達(dá)到10 cm以內(nèi)定義濾波收斂。對(duì)BDS/GPS組合系統(tǒng)分別應(yīng)用LS、EKF、AEKF 3種方法進(jìn)行解算,得到實(shí)驗(yàn)結(jié)果如圖1~圖4以及表3所示。

    圖1 3種解算方法X方向誤差圖Fig.1 X-component error of three methods

    圖2 3種解算方法Y方向誤差圖Fig.2 Y-component error of three methods

    圖3 3種解算方法Z方向誤差圖Fig.3 Z-component error of three methods

    對(duì)比圖1~圖3中組合系統(tǒng)3種解算方法的位置誤差曲線,可以看出,AEKF算法較之LS、EKF在定位精度及收斂速度方面有很大提升。LS方法由于沒有將不同時(shí)刻的定位值聯(lián)系起來相互制約,因而定位結(jié)果通常顯得相當(dāng)粗糙、雜亂,無法實(shí)現(xiàn)收斂;EKF方法由于初值及先驗(yàn)信息估計(jì)不夠準(zhǔn)確,導(dǎo)致收斂時(shí)間增加,收斂后的定位精度不夠高;而AEKF方法由于通過LS方法得到了較為準(zhǔn)確的濾波初值,濾波時(shí)間較短,定位結(jié)果較為理想。

    圖4 3種定位方式衛(wèi)星數(shù)及PDOP值Fig.4 Numbers of satellites and PDOP value of different systems

    表3 BDS、GPS、BDS/GPS組合定位性能比較

    對(duì)比圖4及表3可以發(fā)現(xiàn),雖然在大部分時(shí)間里BDS系統(tǒng)的衛(wèi)星數(shù)要多于GPS系統(tǒng),但是定位精度因子(positiondilutionofprecision,PDOP)有時(shí)沒有GPS衛(wèi)星的小,這主要是由于目前BDS系統(tǒng)的中地球軌道衛(wèi)星較少,衛(wèi)星空間結(jié)構(gòu)不好造成的。目前BDS系統(tǒng)正處于全球組網(wǎng)階段,其跟蹤站數(shù)量有限,所提供的BDS衛(wèi)星的精密軌道數(shù)據(jù)精度還不夠高,而PPP解的優(yōu)劣與衛(wèi)星軌道和衛(wèi)星鐘差修正的精度有很大關(guān)系,并且,目前BDS系統(tǒng)只有14顆可用衛(wèi)星,而GPS系統(tǒng)有32顆,其幾何強(qiáng)度要遠(yuǎn)比BDS好。因此當(dāng)前條件下BDSPPP的收斂時(shí)間明顯長(zhǎng)于GPSPPP[7]。而BDS/GPS組合系統(tǒng)相較單系統(tǒng)來說,極大提高了可用衛(wèi)星的數(shù)目,從而改善了空間衛(wèi)星幾何分布。另外由圖4可以看出,組合系統(tǒng)對(duì)PDOP值有平滑作用,使得數(shù)據(jù)處理收斂效果較好。這也是組合系統(tǒng)能夠加快收斂時(shí)間的主要原因。定量分析發(fā)現(xiàn),靜態(tài)條件下,BDS系統(tǒng)PPP的X分量大部分優(yōu)于5cm,Y分量大部分優(yōu)于6cm,Z分量大部分優(yōu)于5cm。GPS系統(tǒng)PPP的X分量大部分優(yōu)于3cm,Y分量大部分優(yōu)于3cm,Z分量大部分優(yōu)于3cm。BDS/GPS組合系統(tǒng)PPP的X分量大部分優(yōu)于1cm,Y分量大部分優(yōu)于2cm,Z分量大部分優(yōu)于2cm。

    GPS系統(tǒng)X、Y、Z 3個(gè)方向上的均方根誤差(rootmeansquare,RMS)統(tǒng)計(jì)均小于BDS系統(tǒng),這是因?yàn)槟壳癇DS衛(wèi)星數(shù)較少,導(dǎo)致衛(wèi)星幾何構(gòu)型和可靠性較差[14]以及精密軌道和精密鐘差的精度相對(duì)較低。而BDS/GPS組合系統(tǒng)的衛(wèi)星數(shù)目增加了約一半,衛(wèi)星數(shù)目的增加使得相應(yīng)的PDOP值有較大改善,定位精度可以得到明顯提高。

    2.2動(dòng)態(tài)模式

    為進(jìn)一步驗(yàn)證算法在動(dòng)態(tài)模式下的性能,在校園環(huán)境內(nèi)開展動(dòng)態(tài)試驗(yàn)?;鶞?zhǔn)站和移動(dòng)站均采用NovAtelPP6接收機(jī),基準(zhǔn)站與移動(dòng)站之間的距離小于1km,整個(gè)動(dòng)態(tài)實(shí)驗(yàn)時(shí)間約為22min,接收機(jī)運(yùn)動(dòng)速度較低(小于10m/s),數(shù)據(jù)采樣率為1s,截止高度角為10°,圖5為移動(dòng)站運(yùn)動(dòng)軌跡圖。在AEKF處理過程中,BDS和GPS載波相位觀測(cè)量的觀測(cè)噪聲標(biāo)準(zhǔn)差均設(shè)為0.003m,偽距觀測(cè)量的觀測(cè)噪聲標(biāo)準(zhǔn)差分別設(shè)為0.6m和0.3m。

    圖5 移動(dòng)站運(yùn)動(dòng)軌跡圖Fig.5 Reference trajectory of rover user

    為評(píng)估動(dòng)態(tài)PPP精度,定位參考坐標(biāo)來自于開源軟件包RTKlib中的雙頻載波相位差分定位技術(shù),定位精度可達(dá)到10cm級(jí)別。

    運(yùn)用LS、EKF、AEKF對(duì)BDS/GPS組合系統(tǒng)進(jìn)行解算,以及運(yùn)用AEKF對(duì)BDS單系統(tǒng)、GPS單系統(tǒng)、BDS/GPS組合系統(tǒng)進(jìn)行分析,得解算結(jié)果如圖6~圖9及表4所示。

    圖6 3種解算方法X方向動(dòng)態(tài)誤差圖Fig.6 X-component kinematic error of three methods

    圖7 3種解算方法Y方向動(dòng)態(tài)誤差圖Fig.7 Y-component kinematic error of three methods

    圖8 3種解算方法Z方向動(dòng)態(tài)誤差圖Fig.8 Z-component kinematic error of three methods

    對(duì)比圖6~圖8中組合系統(tǒng)3種解算方法的位置誤差曲線,可以看出,AEKF算法較之LS、EKF在定位誤差及收斂速度方面有很大提升。與靜態(tài)的情況相比,動(dòng)態(tài)模式較為復(fù)雜,實(shí)際系統(tǒng)中,接收機(jī)的運(yùn)動(dòng)狀態(tài)往往是未知的,同時(shí)運(yùn)動(dòng)的背景環(huán)境也可能是比較復(fù)雜的。因此在動(dòng)態(tài)高采樣率的情況下,LS方法由于忽略了各歷元之間的關(guān)系,各歷元定位結(jié)果之間存在突變情況,無法實(shí)現(xiàn)收斂,另外由圖9中PDOP值的突變,解釋了LS方法中定位誤差的突變情況;同樣地,在動(dòng)態(tài)模式下,EKF對(duì)初值及先驗(yàn)信息的估計(jì)準(zhǔn)確度更為敏感,由于先驗(yàn)信息估計(jì)不夠準(zhǔn)確,使得EKF估計(jì)方法收斂時(shí)間增加,收斂后的定位結(jié)果還不夠好;AEKF的濾波初值及先驗(yàn)信息由最小二乘法得到,與通過經(jīng)驗(yàn)值進(jìn)行濾波初值設(shè)置的方法相比,本文所提AEKF算法的濾波性能相對(duì)更好,這在一定程度上彌補(bǔ)了運(yùn)動(dòng)的復(fù)雜性及高動(dòng)態(tài)性,使得濾波時(shí)間相對(duì)較短,定位結(jié)果較為理想。

    圖9 3種動(dòng)態(tài)定位方式衛(wèi)星數(shù)及PDOP值Fig.9 Satellite number and PDOP in the kinematic mode

    對(duì)比圖9及表4可以發(fā)現(xiàn),BDS系統(tǒng)的衛(wèi)星個(gè)數(shù)要少于GPS系統(tǒng),BDS對(duì)應(yīng)PDOP值要明顯大于GPS,然而隨著BDS組網(wǎng)建設(shè),將大大增加中軌道衛(wèi)星的數(shù)目,因此BDS衛(wèi)星的空間幾何結(jié)構(gòu)將會(huì)得到有效改善。此外,BDS測(cè)距精度和信號(hào)穩(wěn)定性較GPS略差[15],相應(yīng)地,在定位精度上BDS系統(tǒng)比GPS系統(tǒng)要差,而BDS/GPS組合系統(tǒng)由于衛(wèi)星數(shù)量上的提高,使得系統(tǒng)的幾何結(jié)構(gòu)及定位可靠性增強(qiáng),PDOP值相對(duì)單系統(tǒng)來說顯著減小,使得定位精度有較大提高。

    表4 BDS、GPS、BDS/GPS組合動(dòng)態(tài)定位性能比較

    從數(shù)值上看,組合系統(tǒng)對(duì)BDS系統(tǒng)PDOP值的改善效果要比GPS大得多,這主要是因?yàn)锽DS尚不完善,衛(wèi)星幾何構(gòu)型還不夠好,因此導(dǎo)致其PDOP值較大。雖然目前BDS系統(tǒng)并不完善,在軌衛(wèi)星較少,但是已經(jīng)可以為GPS觀測(cè)提供冗余觀測(cè)量,對(duì)GPS幾何構(gòu)型的改善效果是顯著的。單一的GPS或者BDS系統(tǒng)PDOP值在一些時(shí)刻變化較大,說明雖然該觀測(cè)時(shí)段內(nèi)系統(tǒng)的可用衛(wèi)星數(shù)平均較多,但是由于衛(wèi)星分布不均勻及信號(hào)失鎖等原因,會(huì)造成單一系統(tǒng)定位結(jié)果不夠準(zhǔn)確,系統(tǒng)組合很好地解決了這個(gè)問題。為分析濾波收斂后的定位精度和穩(wěn)定性,對(duì)結(jié)果進(jìn)行分析,可以發(fā)現(xiàn),動(dòng)態(tài)條件下,收斂之后BDS系統(tǒng)PPP的X分量?jī)?yōu)于11cm,Y分量?jī)?yōu)于12cm,Z分量?jī)?yōu)于6cm。GPS系統(tǒng)PPP的X分量?jī)?yōu)于9cm,Y分量?jī)?yōu)于10cm,Z分量?jī)?yōu)于5cm。BDS/GPS組合系統(tǒng)精密單點(diǎn)定位的X分量?jī)?yōu)于4cm,Y分量?jī)?yōu)于6cm,Z分量?jī)?yōu)于4cm。

    與靜態(tài)模式相比,動(dòng)態(tài)模式下,由于接收機(jī)的運(yùn)動(dòng)以及較高的采樣率,使得對(duì)系統(tǒng)狀態(tài)的估計(jì)更為困難,采用多系統(tǒng)組合相比單系統(tǒng)可以成倍地增加可視衛(wèi)星數(shù),衛(wèi)星數(shù)目的增加使得相應(yīng)的PDOP值有較大改善。另外,衛(wèi)星數(shù)目的增加使得量測(cè)量增加,提高了量測(cè)冗余度。對(duì)比靜態(tài)模式,可以看出,動(dòng)態(tài)模式下,組合系統(tǒng)在性能上較之單系統(tǒng)幅度提高更大,效果更為明顯。系統(tǒng)定位可靠性、連續(xù)性、精度都有了顯著的改善。因此,多系統(tǒng)聯(lián)合精密定位是可行的。

    3 結(jié) 論

    采用AEKF的方法,本文針對(duì)BDS/GPS組合PPP進(jìn)行了改進(jìn)研究與實(shí)現(xiàn)。一方面,與單系統(tǒng)PPP結(jié)果的對(duì)比分析發(fā)現(xiàn),組合系統(tǒng)較單一系統(tǒng)在定位精度與收斂速度上都有顯著提高,而且組合系統(tǒng)具有更多的可用衛(wèi)星,能夠保證定位的可靠性、可用性和連續(xù)性。另一方面,將AEKF應(yīng)用于組合系統(tǒng),相比傳統(tǒng)方法,組合系統(tǒng)性能達(dá)到最優(yōu)化,對(duì)動(dòng)態(tài)的適應(yīng)性有很大提高。

    從實(shí)驗(yàn)結(jié)果來看,雖然單獨(dú)BDS系統(tǒng)可以實(shí)現(xiàn)高精度定位,但是多個(gè)衛(wèi)星導(dǎo)航系統(tǒng)聯(lián)合定位,可以進(jìn)一步提升精密定位的性能,本文的研究可為多個(gè)GNSS組合情況下的PPP提供一定的參考。

    [1]ZumbergeJF.AutomatedGPSdataanalysisservice[J].GPS Solutions, 1999, 2(3): 76-78.

    [2]ZhangXH,LiuJN,ForsbergR.Applicationofprecisepointpositioninginairbornesurvey[J].Geomatics and Information Science of Wuhan University, 2006, 31(1): 19-22.

    [3]LiX,GeM,ZhangH,etal.TheGFZreal-timeGNSSprecisepositioningservicesystemanditsadaptionforCOMPASS[J].Advances in Space Research, 2013, 51(6): 1008-1018.

    [4]HeF,GaoCF,PanSG.ApplicationofrobustHelmertvariancecomponentestimationinGPS/BDS/GLONASSintegratedpositioning[J].Journal of Navigation and Positioning, 2013, 1(3): 56-61.

    [5]ZhaoL,LiL,SunM,etal.NoveladaptivehatchfiltertomitigatetheeffectsofionosphereandmultipathonLAAS[J].Journal of Systems Engineering and Electronics, 2010, 21(6): 1046-1053.

    [6]ZhaoL,WangXX,LiL,etal. Nonlinear system filtering theory[M].Beijing:NationalDefenseIndustrialPress,2012.(趙琳,王小旭,李亮,等.非線性系統(tǒng)濾波理論[M].北京:國(guó)防工業(yè)出版社,2012.)

    [7]ZhangXH,ZuoX,LiP.ConvergencetimeandpositioningaccuracycomparisonbetweenBDSandGPSprecisepointpositioning[J].Acta Geodaetica et Cartographica, 2015, 44(3): 250-256.

    [8]MengXG,GuoJM.ResearchonGPS/GLONASSintegratedprecisepointpositioning[J].Geomatics and Information Science of Wuhan University, 2010, 35(12): 1409-1413.

    [9]LiL.Groundbasedaugmentationsystempositioningandintegritymonitoring[D].Harbin:HarbinEngineeringUniversity,2012. (李亮. 陸基增強(qiáng)系統(tǒng)定位與完好性監(jiān)測(cè)技術(shù)研究[D]. 哈爾濱: 哈爾濱工程大學(xué), 2012.)

    [10]SouzaMFD,LuoZQ.Finite-horizonrobustKalmanfilterdesign[J].IEEE Trans. on Signal Processing, 2001, 49(9): 2103-2112.

    [11]ZhangSX. Technology and application of GPS[M].Beijing:NationalDefenseIndustrialPress, 2004. (張守信.GPS技術(shù)與應(yīng)用[M]. 北京: 國(guó)防工業(yè)出版社, 2004.)

    [12]YeSR.TheoryanditsGPSprecisepointpositioningusingun-differencedphaseobservation[D].Wuhan:WuhanUniversity, 2002.(葉世榕.GPS非差精密單點(diǎn)定位理論與實(shí)現(xiàn)[D]. 武漢: 武漢大學(xué), 2002.)

    [13]YueXK,YuanJP.AnadaptiveKalmanfilteringalgorithmbasedonmaximum-likelihoodcriterion[J].Journal of Northwestem Polytechnical University, 2005, 23(4): 469-474.

    [14]JinB,YangSW,LiuWK.AnalysisofGPS/BDSintegratedsinglepointpositioning[J].Hydrographic Surveying and Charting, 2013, 33(4): 39-41, 65.

    [15]TangWM,XuK,JinL,etal.TestandanalysisofBeidou/GPScombinedpseudorangepointpositioningperformance[J].Geomatics and Information Science of Wuhan University, 2015, 40(4): 529-533. (唐衛(wèi)明, 徐坤, 金蕾, 等. 北斗/GPS組合偽距單點(diǎn)定位性能測(cè)試和分析[J].武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2015, 40(4): 529-533.)

    李亮(1984-),通信作者,男,講師,博士,主要研究方向?yàn)樾l(wèi)星導(dǎo)航定位與完好性監(jiān)測(cè)。

    E-mail:liliang@hrbeu.edu.cn

    王雪(1992-),女,碩士研究生,主要研究方向?yàn)榫軉吸c(diǎn)定位。

    E-mail:wangxue_407@hrbeu.edu.cn

    BDS/GPS integrated precise point positioning based on adaptive extended Kalman filter

    ZHAO Lin, ZHANG Sheng-zong, LI Liang, WANG Xue

    (CollegeofAutomation,HarbinEngineeringUniversity,Harbin150001,China)

    In order to achieve the shorter convergence period and better positioning accuracy for precise point positioning (PPP), the combination of multiple navigation satellite, such as the Beidou navigation satellite system (BDS) and the global positioning system (GPS), and the better positioning method are two ways to choose. However, since the traditional least square (LS) estimator ignores the temporal correlation of observations, and the extended Kalman filter (EKF) is limited by inaccurate prior information, a novel PPP method based on adaptive extended Kalman filter (AEKF) is used to adjust the process noise with the measurement consistency test. The experiment results show that, compared with traditional EKF solution, the convergence period of PPP based on AEKF can be shorten by 9 minutes, and the accuracy of positioning can be increased by 33.7%.

    precise point positioning (PPP); convergence period; positioning accuracy; adaptive extended Kalman filter (AEKF)

    2015-10-14;

    2015-12-21;網(wǎng)絡(luò)優(yōu)先出版日期:2016-02-18。

    國(guó)家自然科學(xué)基金(61304235,61273081, 61304234);中央高?;究蒲袠I(yè)務(wù)費(fèi)(HEUCF160416);黑龍江省博士后科研啟動(dòng)基金(LBH-Q15033, LBH-Q14054) 資助課題

    U 666.1

    A

    10.3969/j.issn.1001-506X.2016.09.25

    趙琳(1968-),男,教授,博士,主要研究方向?yàn)樾l(wèi)星導(dǎo)航、慣性導(dǎo)航。

    E-mail:zhaolin@hrbeu.edu.cn

    張勝宗(1988-),男,碩士,主要研究方向?yàn)樾l(wèi)星導(dǎo)航、精密單點(diǎn)定位。

    E-mail:zhangshengzong@hrbeu.edu.cn

    網(wǎng)絡(luò)優(yōu)先出版地址:http://www.cnki.net/kcms/detail/11.2422.TN.20160218.1209.004.html

    猜你喜歡
    分量濾波觀測(cè)
    觀測(cè)到恒星死亡瞬間
    軍事文摘(2023年18期)2023-11-03 09:45:42
    帽子的分量
    一物千斤
    智族GQ(2019年9期)2019-10-28 08:16:21
    論《哈姆雷特》中良心的分量
    分量
    天測(cè)與測(cè)地VLBI 測(cè)地站周圍地形觀測(cè)遮掩的討論
    可觀測(cè)宇宙
    太空探索(2016年7期)2016-07-10 12:10:15
    高分辨率對(duì)地觀測(cè)系統(tǒng)
    太空探索(2015年8期)2015-07-18 11:04:44
    RTS平滑濾波在事后姿態(tài)確定中的應(yīng)用
    基于線性正則變換的 LMS 自適應(yīng)濾波
    国产99白浆流出| 亚洲 欧美 日韩 在线 免费| 精品高清国产在线一区| 日日爽夜夜爽网站| 亚洲欧美日韩高清在线视频| 久久久久精品国产欧美久久久| 亚洲avbb在线观看| 搞女人的毛片| 亚洲国产毛片av蜜桃av| 久久精品国产亚洲av高清一级| bbb黄色大片| 黄色a级毛片大全视频| 国产精品98久久久久久宅男小说| 亚洲精品久久成人aⅴ小说| 性色av乱码一区二区三区2| 日韩一卡2卡3卡4卡2021年| 日日夜夜操网爽| 亚洲欧美精品综合一区二区三区| 国产蜜桃级精品一区二区三区| 免费看美女性在线毛片视频| 亚洲成a人片在线一区二区| 一级毛片精品| 97碰自拍视频| 亚洲九九香蕉| 成人三级做爰电影| 亚洲专区字幕在线| 啦啦啦免费观看视频1| 免费看美女性在线毛片视频| 在线av久久热| 欧美乱色亚洲激情| 国产xxxxx性猛交| 少妇裸体淫交视频免费看高清 | 欧美黑人欧美精品刺激| 如日韩欧美国产精品一区二区三区| 国产av一区在线观看免费| 亚洲最大成人中文| 国产蜜桃级精品一区二区三区| 亚洲国产精品成人综合色| 中文字幕色久视频| 曰老女人黄片| 1024香蕉在线观看| 99在线视频只有这里精品首页| 欧美大码av| 一二三四在线观看免费中文在| 日韩 欧美 亚洲 中文字幕| 精品国产乱码久久久久久男人| 女警被强在线播放| 91九色精品人成在线观看| 韩国精品一区二区三区| 国产亚洲精品久久久久5区| 97超级碰碰碰精品色视频在线观看| 日日爽夜夜爽网站| 亚洲美女黄片视频| 成人亚洲精品一区在线观看| 十八禁网站免费在线| 大陆偷拍与自拍| 国产成+人综合+亚洲专区| 黄色a级毛片大全视频| 99久久99久久久精品蜜桃| 曰老女人黄片| 成人av一区二区三区在线看| 男人舔女人的私密视频| 女人高潮潮喷娇喘18禁视频| 亚洲美女黄片视频| 色在线成人网| 宅男免费午夜| 久久精品人人爽人人爽视色| 亚洲精品中文字幕一二三四区| 丝袜人妻中文字幕| 国产极品粉嫩免费观看在线| 亚洲成a人片在线一区二区| 成人18禁在线播放| 欧美 亚洲 国产 日韩一| 老熟妇仑乱视频hdxx| 亚洲成人国产一区在线观看| 99精品久久久久人妻精品| 免费在线观看日本一区| 纯流量卡能插随身wifi吗| 亚洲熟妇熟女久久| 亚洲一区高清亚洲精品| 亚洲av五月六月丁香网| 亚洲成人国产一区在线观看| 国产熟女午夜一区二区三区| 热re99久久国产66热| 久久久久亚洲av毛片大全| 亚洲性夜色夜夜综合| 两人在一起打扑克的视频| 国产熟女xx| 午夜福利18| 性色av乱码一区二区三区2| 欧美成人免费av一区二区三区| 亚洲午夜精品一区,二区,三区| 男女下面进入的视频免费午夜 | 亚洲黑人精品在线| 天堂影院成人在线观看| 色哟哟哟哟哟哟| 国产精品二区激情视频| 可以在线观看的亚洲视频| 国产成人免费无遮挡视频| 日韩精品青青久久久久久| 亚洲性夜色夜夜综合| 老司机午夜福利在线观看视频| 色婷婷久久久亚洲欧美| 日韩欧美三级三区| 国产精品秋霞免费鲁丝片| 欧美乱码精品一区二区三区| 啦啦啦免费观看视频1| av超薄肉色丝袜交足视频| 国产av一区在线观看免费| 午夜两性在线视频| 久久久国产精品麻豆| 免费不卡黄色视频| 国语自产精品视频在线第100页| 免费av毛片视频| 久久久久九九精品影院| 免费在线观看日本一区| 巨乳人妻的诱惑在线观看| 亚洲国产欧美一区二区综合| 精品国内亚洲2022精品成人| 亚洲欧美日韩无卡精品| 午夜福利视频1000在线观看 | 叶爱在线成人免费视频播放| √禁漫天堂资源中文www| 国产视频一区二区在线看| 人人妻人人爽人人添夜夜欢视频| 国产成人一区二区三区免费视频网站| 美女国产高潮福利片在线看| 日韩欧美三级三区| 国产区一区二久久| 欧美精品亚洲一区二区| 亚洲熟妇熟女久久| 人妻久久中文字幕网| 在线观看66精品国产| 在线免费观看的www视频| 国产熟女午夜一区二区三区| 超碰成人久久| 在线永久观看黄色视频| 国产麻豆成人av免费视频| 亚洲午夜理论影院| 久久精品aⅴ一区二区三区四区| 丝袜人妻中文字幕| 咕卡用的链子| 国产主播在线观看一区二区| 中文字幕人妻熟女乱码| 国产精品亚洲美女久久久| 欧美在线黄色| 久99久视频精品免费| 精品久久久精品久久久| 一区二区日韩欧美中文字幕| 国产av又大| 丝袜美足系列| 少妇被粗大的猛进出69影院| 很黄的视频免费| 欧美绝顶高潮抽搐喷水| 国产精品二区激情视频| 精品久久蜜臀av无| 亚洲一区二区三区色噜噜| 视频区欧美日本亚洲| 可以免费在线观看a视频的电影网站| 国产精品一区二区精品视频观看| 亚洲av电影在线进入| 久久草成人影院| 精品福利观看| 久久久久久大精品| www.熟女人妻精品国产| 亚洲精品在线观看二区| 18禁裸乳无遮挡免费网站照片 | 欧美日韩瑟瑟在线播放| 国产欧美日韩精品亚洲av| 日本vs欧美在线观看视频| 欧美不卡视频在线免费观看 | 午夜福利欧美成人| 国产精华一区二区三区| 亚洲一区二区三区不卡视频| 日韩欧美三级三区| 欧美丝袜亚洲另类 | 欧美日韩一级在线毛片| 激情在线观看视频在线高清| 成人免费观看视频高清| 丰满人妻熟妇乱又伦精品不卡| 一进一出抽搐gif免费好疼| 日本免费a在线| 黄色丝袜av网址大全| 国产精品综合久久久久久久免费 | 成人精品一区二区免费| 美女午夜性视频免费| 国产色视频综合| 在线观看舔阴道视频| 久99久视频精品免费| 国产免费男女视频| 亚洲av成人av| 国产精品久久电影中文字幕| 999精品在线视频| 美女高潮到喷水免费观看| 欧美成人免费av一区二区三区| 激情视频va一区二区三区| 嫩草影院精品99| 免费久久久久久久精品成人欧美视频| 99国产极品粉嫩在线观看| av超薄肉色丝袜交足视频| 中文字幕av电影在线播放| 国产成人系列免费观看| 欧美日韩一级在线毛片| 欧美成人午夜精品| 一夜夜www| 国产日韩一区二区三区精品不卡| 啦啦啦免费观看视频1| 欧美色视频一区免费| 久久精品亚洲精品国产色婷小说| 亚洲视频免费观看视频| 欧美中文日本在线观看视频| 久久久久久人人人人人| 久久久久九九精品影院| 亚洲一区中文字幕在线| 成人亚洲精品一区在线观看| 91老司机精品| 在线观看一区二区三区| 非洲黑人性xxxx精品又粗又长| 午夜激情av网站| 女生性感内裤真人,穿戴方法视频| 精品欧美一区二区三区在线| 18禁观看日本| 久久天躁狠狠躁夜夜2o2o| 一级毛片精品| 丰满的人妻完整版| 精品欧美一区二区三区在线| 久久欧美精品欧美久久欧美| 国产精华一区二区三区| 91精品国产国语对白视频| 动漫黄色视频在线观看| 黄色片一级片一级黄色片| 国产97色在线日韩免费| 精品国产乱子伦一区二区三区| 久久精品91蜜桃| 日韩欧美一区视频在线观看| 久久久久久久久久久久大奶| 欧美成人免费av一区二区三区| 少妇熟女aⅴ在线视频| 亚洲片人在线观看| 亚洲狠狠婷婷综合久久图片| 美女扒开内裤让男人捅视频| 美女国产高潮福利片在线看| 日本免费一区二区三区高清不卡 | 国产成人影院久久av| 亚洲精品av麻豆狂野| 亚洲美女黄片视频| 1024视频免费在线观看| 欧美一级a爱片免费观看看 | 亚洲国产中文字幕在线视频| 老熟妇乱子伦视频在线观看| 国产av一区二区精品久久| 精品免费久久久久久久清纯| 日本vs欧美在线观看视频| 黄色毛片三级朝国网站| 亚洲一区二区三区色噜噜| 国产成年人精品一区二区| 亚洲精品久久国产高清桃花| 免费久久久久久久精品成人欧美视频| 久久精品亚洲熟妇少妇任你| 搡老岳熟女国产| 少妇裸体淫交视频免费看高清 | 久久久久久人人人人人| 女生性感内裤真人,穿戴方法视频| 久久久国产成人免费| 日日夜夜操网爽| 亚洲专区中文字幕在线| 国产黄a三级三级三级人| 美女 人体艺术 gogo| 一本大道久久a久久精品| 久热这里只有精品99| 中国美女看黄片| 精品高清国产在线一区| 黄色视频不卡| 淫妇啪啪啪对白视频| 亚洲精品久久成人aⅴ小说| 免费不卡黄色视频| 天天躁夜夜躁狠狠躁躁| a在线观看视频网站| 日日摸夜夜添夜夜添小说| 国产精品免费视频内射| 欧美另类亚洲清纯唯美| 亚洲一区二区三区不卡视频| 国产高清视频在线播放一区| 国产黄a三级三级三级人| 精品无人区乱码1区二区| 免费一级毛片在线播放高清视频 | 国产精品自产拍在线观看55亚洲| 国产麻豆69| 一本综合久久免费| 搡老岳熟女国产| 老鸭窝网址在线观看| 正在播放国产对白刺激| 日韩免费av在线播放| 一个人免费在线观看的高清视频| 国产97色在线日韩免费| 国产精品国产高清国产av| 国产免费av片在线观看野外av| 给我免费播放毛片高清在线观看| 在线国产一区二区在线| 国产又爽黄色视频| 好男人在线观看高清免费视频 | 99国产极品粉嫩在线观看| 亚洲色图 男人天堂 中文字幕| 亚洲欧美精品综合一区二区三区| 国产精品 国内视频| 少妇熟女aⅴ在线视频| 黑丝袜美女国产一区| 99精品在免费线老司机午夜| 国产一区二区三区综合在线观看| 国内精品久久久久精免费| 久久天堂一区二区三区四区| 婷婷丁香在线五月| 亚洲国产看品久久| 99久久综合精品五月天人人| 亚洲欧美日韩无卡精品| 成年女人毛片免费观看观看9| 好男人电影高清在线观看| 午夜福利视频1000在线观看 | 99久久国产精品久久久| 久久久久久人人人人人| 母亲3免费完整高清在线观看| 黄色a级毛片大全视频| 免费av毛片视频| 无人区码免费观看不卡| 侵犯人妻中文字幕一二三四区| 51午夜福利影视在线观看| 99久久99久久久精品蜜桃| 久久香蕉国产精品| 精品久久久久久久人妻蜜臀av | 最新美女视频免费是黄的| 国产蜜桃级精品一区二区三区| 人人妻人人澡人人看| 少妇粗大呻吟视频| 日韩有码中文字幕| 此物有八面人人有两片| 亚洲人成伊人成综合网2020| 人人妻,人人澡人人爽秒播| 亚洲人成电影免费在线| 欧美+亚洲+日韩+国产| 日韩有码中文字幕| 黑人巨大精品欧美一区二区蜜桃| av福利片在线| 在线观看www视频免费| 18美女黄网站色大片免费观看| 欧美绝顶高潮抽搐喷水| 高清黄色对白视频在线免费看| 好看av亚洲va欧美ⅴa在| 免费高清在线观看日韩| 99国产精品99久久久久| 久久久久久大精品| 咕卡用的链子| 精品免费久久久久久久清纯| 日韩精品中文字幕看吧| 69av精品久久久久久| 亚洲自偷自拍图片 自拍| 欧美激情高清一区二区三区| 久久精品亚洲精品国产色婷小说| 亚洲七黄色美女视频| 欧美绝顶高潮抽搐喷水| 国产成人精品无人区| 色老头精品视频在线观看| 免费观看人在逋| 亚洲专区字幕在线| 国产成人精品久久二区二区91| 欧美性长视频在线观看| 99国产精品一区二区蜜桃av| 九色亚洲精品在线播放| 欧美成人一区二区免费高清观看 | 美女 人体艺术 gogo| www.自偷自拍.com| 亚洲成a人片在线一区二区| 国产男靠女视频免费网站| 国产麻豆69| 亚洲精品在线观看二区| 久久精品国产99精品国产亚洲性色 | 999精品在线视频| 国产av又大| 日韩中文字幕欧美一区二区| 黄频高清免费视频| 男男h啪啪无遮挡| 欧美乱色亚洲激情| 日韩大尺度精品在线看网址 | 在线av久久热| 老熟妇仑乱视频hdxx| 欧美成人一区二区免费高清观看 | 久久人人精品亚洲av| 一级作爱视频免费观看| 精品少妇一区二区三区视频日本电影| 露出奶头的视频| 日本在线视频免费播放| 精品久久久久久,| 亚洲成av人片免费观看| 两性午夜刺激爽爽歪歪视频在线观看 | 黄色片一级片一级黄色片| 亚洲色图av天堂| 欧美激情 高清一区二区三区| 国产亚洲欧美98| 欧美日韩黄片免| 大型av网站在线播放| 黄色毛片三级朝国网站| 一进一出好大好爽视频| 首页视频小说图片口味搜索| 亚洲精品国产区一区二| 久久久久久久久久久久大奶| 韩国av一区二区三区四区| 一二三四社区在线视频社区8| 欧美日韩亚洲国产一区二区在线观看| 999久久久国产精品视频| 日韩高清综合在线| 亚洲最大成人中文| 亚洲av电影在线进入| 淫妇啪啪啪对白视频| 日韩欧美一区视频在线观看| 欧美丝袜亚洲另类 | 最近最新免费中文字幕在线| 欧美成人午夜精品| 亚洲人成电影免费在线| 91成人精品电影| 亚洲欧美精品综合久久99| 中文字幕人成人乱码亚洲影| 国产片内射在线| 亚洲成人国产一区在线观看| 丁香六月欧美| 免费女性裸体啪啪无遮挡网站| 曰老女人黄片| 欧美乱妇无乱码| 日韩免费av在线播放| 免费不卡黄色视频| 亚洲精品久久成人aⅴ小说| 国产高清videossex| 在线免费观看的www视频| 欧美最黄视频在线播放免费| 99国产综合亚洲精品| www日本在线高清视频| 夜夜爽天天搞| 亚洲一区二区三区色噜噜| 亚洲情色 制服丝袜| 亚洲,欧美精品.| 看免费av毛片| 欧美日韩亚洲综合一区二区三区_| 久久精品亚洲熟妇少妇任你| www.www免费av| 亚洲成人久久性| 午夜老司机福利片| 丁香欧美五月| 国产高清videossex| 天堂动漫精品| 精品久久久精品久久久| 日本 欧美在线| 国产一级毛片七仙女欲春2 | 国产伦一二天堂av在线观看| 亚洲国产欧美一区二区综合| 悠悠久久av| 桃色一区二区三区在线观看| 亚洲国产看品久久| 91老司机精品| 亚洲九九香蕉| 侵犯人妻中文字幕一二三四区| 后天国语完整版免费观看| 波多野结衣高清无吗| 少妇的丰满在线观看| 欧美黑人精品巨大| 午夜免费鲁丝| 日韩欧美免费精品| 国产视频一区二区在线看| 精品熟女少妇八av免费久了| 搡老岳熟女国产| 国产精品一区二区在线不卡| 视频区欧美日本亚洲| 少妇熟女aⅴ在线视频| 日韩大码丰满熟妇| 久久人人97超碰香蕉20202| 性欧美人与动物交配| 日韩欧美一区二区三区在线观看| 欧美最黄视频在线播放免费| 男女之事视频高清在线观看| 国产视频一区二区在线看| 国产极品粉嫩免费观看在线| 99re在线观看精品视频| 50天的宝宝边吃奶边哭怎么回事| 久久中文字幕一级| 亚洲人成网站在线播放欧美日韩| 美女 人体艺术 gogo| 午夜福利欧美成人| 人妻丰满熟妇av一区二区三区| 淫秽高清视频在线观看| 国产精品一区二区精品视频观看| 欧美日韩福利视频一区二区| 日本欧美视频一区| 欧美久久黑人一区二区| 大陆偷拍与自拍| 黄色成人免费大全| 国产麻豆成人av免费视频| 精品国产一区二区三区四区第35| 我的亚洲天堂| 性色av乱码一区二区三区2| 亚洲熟女毛片儿| 99久久国产精品久久久| 一区二区日韩欧美中文字幕| 国产99久久九九免费精品| 午夜亚洲福利在线播放| 亚洲精品粉嫩美女一区| 亚洲专区字幕在线| 99riav亚洲国产免费| 此物有八面人人有两片| 欧美国产日韩亚洲一区| 久久久国产精品麻豆| 亚洲色图av天堂| 99香蕉大伊视频| 一级a爱片免费观看的视频| 国产精品一区二区免费欧美| 精品国产乱码久久久久久男人| 国产亚洲av高清不卡| 很黄的视频免费| av视频免费观看在线观看| 看片在线看免费视频| 亚洲黑人精品在线| 天天躁夜夜躁狠狠躁躁| 午夜久久久在线观看| 长腿黑丝高跟| 日韩精品中文字幕看吧| 精品熟女少妇八av免费久了| 午夜a级毛片| 可以在线观看毛片的网站| 性色av乱码一区二区三区2| 久久久久久免费高清国产稀缺| 久久久久久亚洲精品国产蜜桃av| 免费在线观看完整版高清| 在线观看舔阴道视频| 这个男人来自地球电影免费观看| 日本五十路高清| www.熟女人妻精品国产| 国产精品自产拍在线观看55亚洲| 一边摸一边做爽爽视频免费| 欧美色欧美亚洲另类二区 | 中文字幕高清在线视频| 久久精品91无色码中文字幕| 人人妻,人人澡人人爽秒播| 亚洲精品国产色婷婷电影| 国产精品自产拍在线观看55亚洲| 亚洲熟女毛片儿| 十分钟在线观看高清视频www| 18禁国产床啪视频网站| 亚洲专区字幕在线| 精品欧美国产一区二区三| 国产真人三级小视频在线观看| 欧美最黄视频在线播放免费| 人成视频在线观看免费观看| 伊人久久大香线蕉亚洲五| 法律面前人人平等表现在哪些方面| 琪琪午夜伦伦电影理论片6080| 精品福利观看| 免费在线观看亚洲国产| 欧美激情 高清一区二区三区| 亚洲专区字幕在线| 久9热在线精品视频| 啦啦啦免费观看视频1| 欧美最黄视频在线播放免费| av网站免费在线观看视频| 大型黄色视频在线免费观看| 国产精品 欧美亚洲| 一个人免费在线观看的高清视频| 国产视频一区二区在线看| 首页视频小说图片口味搜索| АⅤ资源中文在线天堂| 亚洲av五月六月丁香网| 很黄的视频免费| 免费在线观看完整版高清| 日韩国内少妇激情av| 一区在线观看完整版| 日韩 欧美 亚洲 中文字幕| 波多野结衣巨乳人妻| 99riav亚洲国产免费| 亚洲国产毛片av蜜桃av| 窝窝影院91人妻| 久久天躁狠狠躁夜夜2o2o| 欧美日本亚洲视频在线播放| 欧美一级a爱片免费观看看 | 男女之事视频高清在线观看| 午夜视频精品福利| www.熟女人妻精品国产| 免费观看人在逋| 国产高清视频在线播放一区| 国产精品永久免费网站| 丝袜美腿诱惑在线| 欧美亚洲日本最大视频资源| 人妻久久中文字幕网| 淫秽高清视频在线观看| 日本黄色视频三级网站网址| 人妻久久中文字幕网| 90打野战视频偷拍视频| 动漫黄色视频在线观看| 欧美另类亚洲清纯唯美| 黑丝袜美女国产一区| 亚洲一区二区三区色噜噜| 国产成人啪精品午夜网站| 极品教师在线免费播放| 亚洲成人免费电影在线观看| 精品午夜福利视频在线观看一区| 国产亚洲欧美在线一区二区| 国产视频一区二区在线看| 中出人妻视频一区二区| 国产精品九九99| 中文亚洲av片在线观看爽| 国产单亲对白刺激| 国产精品日韩av在线免费观看 | x7x7x7水蜜桃| 人人妻人人爽人人添夜夜欢视频| 老熟妇乱子伦视频在线观看| 亚洲成a人片在线一区二区| 国产一区二区三区综合在线观看| 亚洲专区国产一区二区| 欧美丝袜亚洲另类 | 亚洲色图av天堂| 国产97色在线日韩免费|