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

    附加原子鐘物理模型的PPP時(shí)間傳遞算法

    2016-12-07 03:18:18于合理郝金明劉偉平田英國(guó)
    測(cè)繪學(xué)報(bào) 2016年11期
    關(guān)鍵詞:原子鐘鐘差協(xié)方差

    于合理,郝金明,劉偉平,田英國(guó),鄧 科

    1. 信息工程大學(xué)導(dǎo)航與空天目標(biāo)工程學(xué)院,河南 鄭州 450001; 2. 北斗導(dǎo)航應(yīng)用技術(shù)河南省協(xié)同創(chuàng)新中心,河南 鄭州 450001

    ?

    附加原子鐘物理模型的PPP時(shí)間傳遞算法

    于合理1,2,郝金明1,2,劉偉平1,2,田英國(guó)1,2,鄧 科1

    1. 信息工程大學(xué)導(dǎo)航與空天目標(biāo)工程學(xué)院,河南 鄭州 450001; 2. 北斗導(dǎo)航應(yīng)用技術(shù)河南省協(xié)同創(chuàng)新中心,河南 鄭州 450001

    傳統(tǒng)精密單點(diǎn)定位(PPP)時(shí)間傳遞算法通常把接收機(jī)鐘差當(dāng)作相互獨(dú)立的白噪聲逐歷元進(jìn)行估計(jì),而忽略了鐘差參數(shù)歷元間的相關(guān)性。針對(duì)這一問(wèn)題,本文提出了一種附加原子鐘物理模型的PPP時(shí)間傳遞算法。該算法通過(guò)利用Kalman濾波對(duì)高穩(wěn)定度的原子鐘鐘差進(jìn)行建模,拓展傳統(tǒng)PPP時(shí)間傳遞模型中的接收機(jī)鐘差參數(shù),并給出了Kalman濾波過(guò)程噪聲協(xié)方差和初始狀態(tài)向量的確定方法。試驗(yàn)結(jié)果表明:該算法可以有效避免傳統(tǒng)算法時(shí)間傳遞結(jié)果需要一定收斂時(shí)間的問(wèn)題,使解算結(jié)果更加符合原子鐘的物理特性,能夠顯著提高時(shí)間傳遞結(jié)果的精度和穩(wěn)定性,可將單站時(shí)間傳遞精度平均提高58%,站間時(shí)間傳遞精度平均提高51%。

    原子鐘;物理模型;精密單點(diǎn)定位;時(shí)間傳遞

    在高精度遠(yuǎn)程時(shí)間傳遞中衛(wèi)星雙向時(shí)間頻率傳遞(two way satellite time and frequency transfer,TWSTFT)和精密單點(diǎn)定位(precise point positioning,PPP)時(shí)間傳遞技術(shù)是精度最高的兩種方法[1-3]。TWSTFT需要專門的中介衛(wèi)星建立通信鏈路,且參與雙方都需要專門的發(fā)射接收設(shè)備,設(shè)備價(jià)格昂貴,目前世界上只有少數(shù)較大的實(shí)驗(yàn)室應(yīng)用。而PPP時(shí)間傳遞算法是單站作業(yè),操作靈活方便,時(shí)間傳遞的精度可達(dá)亞納秒量級(jí),其短期穩(wěn)定度甚至優(yōu)于TWSTFT技術(shù)[1-3],已經(jīng)被廣泛應(yīng)用于時(shí)間傳遞領(lǐng)域。隨著國(guó)際GNSS服務(wù)組織(International GNSS Services,IGS)和國(guó)際時(shí)頻機(jī)構(gòu)的深入合作,許多實(shí)驗(yàn)室已經(jīng)將PPP方法用于時(shí)間傳遞和原子時(shí)的計(jì)算[3],國(guó)際計(jì)量局還將其用于維護(hù)協(xié)調(diào)世界時(shí)的運(yùn)作[1]。文獻(xiàn)[1—5]對(duì)PPP時(shí)間傳遞技術(shù)進(jìn)行了研究,并取得了良好效果。文獻(xiàn)[6—7]還對(duì)非組合PPP時(shí)間傳遞方法進(jìn)行了研究,該方法能有效避免傳統(tǒng)PPP算法由于觀測(cè)值組合過(guò)程中觀測(cè)噪聲及多路徑效應(yīng)誤差被放大而造成的對(duì)接收機(jī)鐘差解算的不利影響,解算結(jié)果優(yōu)于傳統(tǒng)PPP算法。傳統(tǒng)的PPP時(shí)間傳遞算法通常是將各歷元接收機(jī)鐘差看作相互獨(dú)立的白噪聲逐歷元進(jìn)行估計(jì),未充分發(fā)掘接收機(jī)鐘差參數(shù)在相鄰歷元間變化可能較為平穩(wěn)的約束信息。同時(shí)由于接收機(jī)鐘差、測(cè)站高程和對(duì)流層參數(shù)的高度相關(guān),解算得到的鐘差序列往往會(huì)偏離接收機(jī)鐘差序列多項(xiàng)式的特性[8-9],時(shí)間傳遞結(jié)果還會(huì)存在一定時(shí)間的收斂過(guò)程。為此,本文利用原子鐘的物理特性對(duì)接收機(jī)鐘差建模來(lái)進(jìn)一步提高PPP時(shí)間傳遞的精度和穩(wěn)定性。

    利用接收機(jī)鐘物理特性對(duì)鐘差進(jìn)行建模的思想并不是最新提出的,美國(guó)噴氣推進(jìn)實(shí)驗(yàn)室(Jet Propulsion Laboratory,JPL)在20世紀(jì)80年代為了提高GNSS軌道等相關(guān)參數(shù)的精度就采用不同的隨機(jī)過(guò)程模型來(lái)對(duì)接收機(jī)鐘差和對(duì)流層進(jìn)行建模[10]。但由于早期的接收機(jī)一般采用穩(wěn)定度較差的石英鐘,即使采用復(fù)雜的數(shù)學(xué)模型也無(wú)法準(zhǔn)確描述石英鐘的運(yùn)行特征,鐘差建模對(duì)解算結(jié)果的影響微乎其微。而目前許多時(shí)間實(shí)驗(yàn)室都已經(jīng)配備了高穩(wěn)定度的原子鐘,這使通過(guò)鐘差物理建模提高時(shí)間傳遞精度變得可能。近年來(lái)眾多學(xué)者對(duì)鐘差建模進(jìn)行了大量的研究,取得了豐碩的成果。在衛(wèi)星原子鐘建模方面,文獻(xiàn)[11]在顧及原子鐘物理模型的基礎(chǔ)上提出了適用于原子鐘的Kalman濾波器;文獻(xiàn)[12—15]進(jìn)一步研究了原子鐘濾波過(guò)程中狀態(tài)矩陣和觀測(cè)噪聲矩陣的確定方法;文獻(xiàn)[16]利用原子鐘物理特性對(duì)鐘差建模實(shí)現(xiàn)了對(duì)銣原子鐘的狀態(tài)監(jiān)測(cè)。在接收機(jī)鐘差建模方面,文獻(xiàn)[17—19]對(duì)低軌衛(wèi)星星載接收機(jī)端的鐘差建模方法進(jìn)行了研究;文獻(xiàn)[9]對(duì)高穩(wěn)定性地面接收機(jī)鐘的隨機(jī)模型進(jìn)行了研究;文獻(xiàn)[8,20—21]對(duì)通過(guò)接收機(jī)鐘物理建模來(lái)提高定位精度和收斂速度進(jìn)行了研究,并取得了一定的效果。目前針對(duì)接收機(jī)鐘差建模提高時(shí)間傳遞性能方面的研究較少,且處于起步階段。文獻(xiàn)[22]通過(guò)接收機(jī)鐘差建模對(duì)接收機(jī)鐘差參數(shù)估計(jì)進(jìn)行了一些研究,但主要是針對(duì)石英晶振接收機(jī),并不能有效提高時(shí)間傳遞精度。文獻(xiàn)[1]也對(duì)基于原子鐘鐘差模型的時(shí)間傳遞算法進(jìn)行了初步的研究。

    因此,本文利用原子鐘的物理特性,采用Kalman濾波方法對(duì)接收機(jī)鐘差進(jìn)行建模,給出了Kalman濾波過(guò)程噪聲協(xié)方差矩陣、初始狀態(tài)向量及其協(xié)方差矩陣的確定方法,并重點(diǎn)分析了附加原子鐘物理模型對(duì)單站時(shí)間傳遞及站間時(shí)間同步結(jié)果的影響。

    1 附加原子鐘物理模型的PPP時(shí)間傳遞算法

    1.1 基本觀測(cè)模型

    對(duì)于雙頻接收機(jī),PPP時(shí)間傳遞算法一般采用消電離層組合觀測(cè)量[1],觀測(cè)方程如下

    (1)

    (2)

    1.2 附加原子鐘物理模型

    原子鐘具有較高的頻率穩(wěn)定度,其鐘差通常采用簡(jiǎn)單的一次多項(xiàng)式或二次多項(xiàng)式進(jìn)行描述[23]。限于篇幅,本文僅以二次多項(xiàng)式為例進(jìn)行分析

    x(t)=a0+a1(t-t0)+a2(t-t0)2+εx(t)

    (3)

    式中,x(t)為t時(shí)刻的接收機(jī)鐘差;a0、a1、a2為確定性分量,分別表示參考時(shí)刻t0的初始鐘差、鐘速、鐘漂,當(dāng)a2=0時(shí),式(3)可簡(jiǎn)化為線性模型;εx(t)為隨機(jī)變化的不確定性分量。

    PPP時(shí)間傳遞算法通常采用Kalman濾波進(jìn)行參數(shù)估計(jì),一般是將接收機(jī)鐘差當(dāng)作獨(dú)立的白噪聲,逐歷元估計(jì)鐘差參數(shù),未充分發(fā)掘接收機(jī)鐘差參數(shù)在相鄰歷元間變化可能較為平穩(wěn)的約束信息。若考慮接收機(jī)鐘差歷元間的關(guān)系,可附加原子鐘的物理模型進(jìn)行約束,拓展傳統(tǒng)PPP模型中的接收機(jī)鐘差參數(shù),使解算的鐘差序列仍然體現(xiàn)原子鐘的物理特性。具體做法是在PPP時(shí)間傳遞基本觀測(cè)方程式(1)與式(2)估計(jì)接收機(jī)鐘差參數(shù)的基礎(chǔ)上,同時(shí)估計(jì)鐘速和鐘漂參數(shù)。此時(shí)觀測(cè)方程中未知參數(shù)為測(cè)站三維坐標(biāo)參數(shù)、鐘差、鐘速、鐘漂、對(duì)流層和模糊度參數(shù)。為了便于推導(dǎo)和直觀地顯示附加鐘差模型的方法,公式推導(dǎo)過(guò)程僅考慮鐘差模型參數(shù),觀測(cè)方程簡(jiǎn)化為

    Lk=AkXk+Δk

    (4)

    式中,Lk為觀測(cè)量;Ak=[100]為設(shè)計(jì)矩陣;Xk=[x(t+τ)y(t+τ)z(t+τ)]T為tk時(shí)刻的三維狀態(tài)向量,x(t)、y(t)、z(t)為接收機(jī)鐘差(時(shí)差)、鐘速(頻率)和鐘漂(頻漂),τ為tk時(shí)刻與tk-1時(shí)刻的時(shí)間間隔;Δk為觀測(cè)噪聲。

    Kalman濾波三維狀態(tài)空間方程為

    (5)

    (6)

    可得狀態(tài)空間方程的簡(jiǎn)化形式

    Xk=Φk,k-1Xk-1+Wk

    (7)

    式中,Φk,k-1為狀態(tài)轉(zhuǎn)移矩陣;Wk為三維過(guò)程噪聲向量,其協(xié)方差矩陣Q可以表示成過(guò)程噪聲參數(shù)的函數(shù)[13]

    (8)

    式中,q1、q2和q3分別可用調(diào)頻白噪聲、調(diào)頻隨機(jī)游走噪聲和調(diào)頻隨機(jī)奔跑噪聲來(lái)描述。

    若原子鐘線性頻漂不明顯,可采用二維狀態(tài)模型來(lái)描述接收機(jī)的動(dòng)態(tài)過(guò)程,此時(shí)狀態(tài)轉(zhuǎn)移矩陣和協(xié)方差陣需要根據(jù)三維狀態(tài)模型進(jìn)行相應(yīng)的調(diào)整。

    2 過(guò)程噪聲協(xié)方差矩陣的確定

    系統(tǒng)過(guò)程噪聲協(xié)方差矩陣的確定是Kalman濾波的關(guān)鍵,可由原子鐘的濾波過(guò)程噪聲參數(shù)來(lái)表示。過(guò)程噪聲參數(shù)可以通過(guò)哈達(dá)瑪總方差或阿倫方差進(jìn)行估計(jì)。

    2.1 基于哈達(dá)瑪總方差的過(guò)程噪聲估計(jì)

    原子鐘頻漂較明顯時(shí),可采用哈達(dá)瑪總方差來(lái)估計(jì)過(guò)程噪聲參數(shù)。基于時(shí)差數(shù)據(jù)的哈達(dá)瑪總方差計(jì)算公式為[5,24]

    (9)

    當(dāng)同時(shí)考慮調(diào)相白噪聲和狀態(tài)方程中存在的3種調(diào)頻噪聲,且假設(shè)調(diào)相白噪聲與其他3種噪聲不相關(guān)時(shí),可推導(dǎo)出公式[5]

    (10)

    式中,q0為對(duì)應(yīng)于調(diào)相白噪聲的測(cè)量噪聲參數(shù)。

    2.2 基于阿倫方差的過(guò)程噪聲估計(jì)

    原子鐘頻漂不明顯時(shí),可以采用阿倫方差來(lái)估計(jì)過(guò)程噪聲參數(shù)。基于時(shí)差數(shù)據(jù)阿倫方差計(jì)算公式為[5,24-25]

    (11)

    當(dāng)同時(shí)考慮調(diào)相白噪聲和狀態(tài)方程中存在的3種調(diào)頻噪聲時(shí),可推導(dǎo)出公式[5]

    (12)

    設(shè)定不同的采樣間隔,就可以得到不同的哈達(dá)瑪總方差(或阿倫方差)公式,聯(lián)立各式利用最小二乘算法即可解得過(guò)程噪聲參數(shù)q0、q1、q2和q3,進(jìn)而確定系統(tǒng)過(guò)程噪聲協(xié)方差矩陣。

    3 初始狀態(tài)向量及協(xié)方差確定

    Kalman濾波算法需要給定初始狀態(tài)向量及其協(xié)方差矩陣。在無(wú)先驗(yàn)信息時(shí),通常令初始狀態(tài)向量為零向量,協(xié)方差矩陣為aI,其中a為很大的正數(shù),I為單位矩陣,以確保濾波結(jié)果的無(wú)偏性。但在沒(méi)有先驗(yàn)信息時(shí),濾波器初始化會(huì)比較困難,且會(huì)影響濾波器的穩(wěn)定性。在進(jìn)行連續(xù)解算或已知先驗(yàn)信息時(shí),若已知3個(gè)時(shí)差測(cè)量值l(0)、l(1)、l(2),τ1、τ2為這3個(gè)觀測(cè)量的采樣間隔,則根據(jù)式(13)—式(21)即可確定狀態(tài)向量初值X(2)=[x(2)y(2)z(2)]以及協(xié)方差矩陣P(2),其中協(xié)方差矩陣為對(duì)角陣[15-16]

    x(2)=l(2)

    (13)

    (14)

    (15)

    P11(2)=q0

    (16)

    (17)

    (18)

    (19)

    (20)

    (21)

    4 試驗(yàn)分析

    由JPL發(fā)布的“l(fā)oglist.txt”文件可知,截至2016年6月2日,IGS已經(jīng)有超過(guò)130個(gè)跟蹤站配備了高頻穩(wěn)度的原子鐘,其中71個(gè)配備了氫原子鐘[26],這為本文算法性能的驗(yàn)證提供了很大的便利。為了驗(yàn)證附加原子鐘物理模型的PPP時(shí)間傳遞算法的應(yīng)用效果,選取2015年1月5日配備了氫原子鐘的WAB2、CRO1、PTBB、NRC1和MDVJ 5個(gè)IGS跟蹤站的觀測(cè)數(shù)據(jù)進(jìn)行試驗(yàn)分析。觀測(cè)數(shù)據(jù)采樣間隔為30 s,采用偽距和載波相位消電離層組合觀測(cè)量,偽距與載波相位的權(quán)比為1/10 000,高度截止角設(shè)為10°。試驗(yàn)采用靜態(tài)模式進(jìn)行處理,利用卡爾曼濾波估計(jì)三維坐標(biāo)、接收機(jī)鐘差、鐘速、對(duì)流層延遲和模糊度參數(shù)。衛(wèi)星軌道和鐘差采用IGS提供的事后精密星歷和鐘差產(chǎn)品,為了避免衛(wèi)星和接收機(jī)硬件延遲偏差對(duì)接收機(jī)鐘差造成的系統(tǒng)性偏差影響,利用CODE發(fā)布的DCB產(chǎn)品對(duì)其進(jìn)行改正[27]。分別采用以下兩種方案進(jìn)行PPP解算。

    方案1:將接收機(jī)鐘差當(dāng)作白噪聲,采用傳統(tǒng)的逐歷元估計(jì)一維鐘差參數(shù)的方法。

    方案2:顧及接收機(jī)鐘差歷元間的相關(guān)性,采用附加原子鐘物理模型的方法進(jìn)行解算,過(guò)程噪聲由阿倫方差根據(jù)IGS事后精密產(chǎn)品解算得到。

    目前IGS提供的事后接收機(jī)鐘差產(chǎn)品精度為75 ps,可將其作為真值進(jìn)行比對(duì)。由于IGS事后精密鐘差產(chǎn)品的采樣間隔為5 min,因此本文以5 min為采樣間隔對(duì)方案1和方案2的時(shí)間傳遞精度進(jìn)行分析。限于篇幅,僅給出了WAB2和CRO1兩跟蹤站采用兩種方案得到的結(jié)果與IGS時(shí)間序列對(duì)比,分別如圖1與圖2所示。同時(shí)將兩種方案的鐘差結(jié)果分別與IGS最終結(jié)果作差可以得到對(duì)應(yīng)的差值序列,圖3與圖4分別給出了WAB2和CRO1兩跟蹤站單站時(shí)間傳遞結(jié)果的差值序列。

    從圖1—圖4可知:

    (1) 方案2的收斂速度明顯快于方案1,單歷元就可以達(dá)到很高的精度。而方案1由于未使用原子鐘物理模型作為先驗(yàn)信息進(jìn)行約束,需要一定時(shí)間的收斂過(guò)程才能達(dá)到較高的精度,WAB2和CRO1兩測(cè)站分別需要22 min、35 min才能收斂。另外,對(duì)方案2得到的坐標(biāo),對(duì)流層和模糊度結(jié)果進(jìn)行分析,其收斂速度均有一定程度的提高。

    圖1 WAB2-GPST兩種方案結(jié)果與IGS時(shí)間序列對(duì)比Fig.1 WAB2-GPST results of the two schemes compared with IGS time series

    圖2 CRO1-GPST兩種方案結(jié)果與IGS時(shí)間序列對(duì)比Fig.2 CRO1-GPST results of the two schemes compared with IGS time series

    圖3 WAB2-GPST兩種方案時(shí)間傳遞結(jié)果與IGS差值序列Fig.3 Difference between WAB2-GPST results and IGS from two schemes

    (2) 方案2由于附加了原子鐘物理模型,得到的接收機(jī)鐘差時(shí)間序列的連續(xù)性和穩(wěn)定性明顯優(yōu)于方案1,與IGS發(fā)布的事后精密鐘差結(jié)果的一致性較高,更能準(zhǔn)確地描述接收機(jī)鐘差序列的變化規(guī)律。

    圖4 CRO1-GPST兩種方案時(shí)間傳遞結(jié)果與IGS差值序列Fig.4 Difference between CRO1-GPST results and IGS from two schemes

    (3) 從圖1明顯可以看出,方案1解算的鐘差序列略微偏離了接收機(jī)鐘本身的多項(xiàng)式特性,這是由于方案1解算的鐘差序列之間是相互獨(dú)立的,解算的過(guò)程中吸收了許多不屬于鐘差的其他誤差源而造成。方案2的濾波結(jié)果更能準(zhǔn)確可靠地反映原子鐘的物理特性。

    表1中給出了WAB2、CRO1、PTBB、NRC1和MDVJ5個(gè)跟蹤站相應(yīng)的誤差統(tǒng)計(jì)結(jié)果。表中SD為標(biāo)準(zhǔn)差(standarddeviation),用來(lái)衡量鐘差誤差序列自身的離散程度,RMS為均方根誤差用來(lái)衡量鐘差估計(jì)值同真值之間的誤差。

    表1 兩種方案結(jié)果時(shí)間傳遞精度

    Tab.1 Accuracy of the time transferring results from different schemes ns

    由表1可知,各測(cè)站方案2時(shí)間傳遞結(jié)果的SD和RMS均優(yōu)于方案1,方案1的鐘差解相對(duì)較為離散,方案2的鐘差解穩(wěn)定度較高,與IGS結(jié)果的符合程度更高。方案2附加原子鐘物理模型后,能夠顯著提高PPP時(shí)間傳遞的精度,可將WAB2、CRO1、PTBB、NRC1和MDVJ 5個(gè)跟蹤站的單站時(shí)間傳遞精度(RMS)分別提高71%、62%、73%、53%和33%,平均提高58%。該算法對(duì)于促進(jìn)PPP時(shí)間傳遞技術(shù)在國(guó)際時(shí)間實(shí)驗(yàn)室等高精度時(shí)間同步領(lǐng)域的應(yīng)用具有一定意義。

    為了進(jìn)一步分析本文算法應(yīng)用于站間時(shí)間同步中的精度,將各站利用PPP獲得的原子鐘的鐘差結(jié)果相減即可得到站間時(shí)間同步結(jié)果,表2給出了兩種方案得到的不同時(shí)間鏈路的時(shí)間同步序列的精度統(tǒng)計(jì)結(jié)果,為更好地觀察時(shí)間同步序列變化情況,圖5以CRO1-WAB2時(shí)間鏈路為例給出了站間時(shí)間同步結(jié)果同IGS結(jié)果的差值序列圖。

    表2 兩種方案時(shí)間同步結(jié)果精度

    Tab.2 Accuracy of the time transferring results from different schemes ns

    圖5 CRO1-WAB2鏈路兩方案時(shí)間同步結(jié)果與IGS結(jié)果差值序列Fig.5 Difference between CRO1-WAB2 results and IGS from two schemes

    根據(jù)表2和圖5可知,方案2得到的站間時(shí)間同步結(jié)果的精度明顯優(yōu)于方案1,與IGS時(shí)間同步結(jié)果一致性較高,能夠可靠地描述測(cè)站間時(shí)間同步序列的變化規(guī)律。由圖5可知,方案1采用傳統(tǒng)算法CRO1-WAB2鏈路需要大約35 min才能收斂,而方案2附加原子鐘物理模型后幾乎不需要收斂時(shí)間即可達(dá)到很高的精度。方案2相比方案1,能夠?qū)RO1-WAB2、CRO1-PTBB、CRO1-NRC1和CRO1-MDVJ鏈路的站間時(shí)間同步精度分別提高51%、49%、29%和73%,平均提高51%。

    5 結(jié) 論

    鑒于當(dāng)前許多國(guó)際時(shí)頻實(shí)驗(yàn)室和IGS跟蹤站均配置有高性能原子鐘,本文通過(guò)附加原子鐘物理模型對(duì)現(xiàn)有PPP時(shí)間傳遞模型進(jìn)行改進(jìn),使用Kalman濾波對(duì)接收機(jī)原子鐘進(jìn)行鐘差建模,研究了Kalman濾波中系統(tǒng)過(guò)程噪聲協(xié)方差矩陣、初始狀態(tài)向量及協(xié)方差的確定方法,并重點(diǎn)分析了附加原子鐘物理模型對(duì)單站時(shí)間傳遞及站間時(shí)間同步結(jié)果的影響。試驗(yàn)結(jié)果表明:附加原子鐘物理模型后,解決了傳統(tǒng)算法中收斂過(guò)程的耗時(shí)問(wèn)題,顯著提高單站時(shí)間傳遞和站間時(shí)間同步的精度,可將單站時(shí)間傳遞精度平均提高58%,尤其是WAB2和PTBB站精度提高水平高達(dá)71%和73%,同時(shí)可將站間時(shí)間同步精度平均提高51%。研究附加原子鐘物理模型的PPP時(shí)間傳遞算法,對(duì)于國(guó)際時(shí)間實(shí)驗(yàn)室實(shí)現(xiàn)高精度的時(shí)間傳遞具有一定的參考意義,對(duì)于拓展GNSS在時(shí)頻傳遞領(lǐng)域的應(yīng)用具有重要意義。本文算法并非適用于所有的GNSS接收機(jī),其前提是要判定接收機(jī)鐘本身的物理特性。它要求接收機(jī)配備有高精度和高穩(wěn)定的原子鐘,否則加入鐘差物理模型后不但不能提高時(shí)間傳遞精度,反而會(huì)引起比較大的誤差。

    [1] 李紅濤. 基于GPS和GLONASS的單站授時(shí)和時(shí)差監(jiān)測(cè)研究[D]. 西安: 長(zhǎng)安大學(xué), 2012.

    LI Hongtao. Research on the Single Station Time Service and Time Offset Based on the GPS and GLONASS Data[D]. Xi’an: Chang’an University, 2012.

    [2] DEFRAIGNE P, AERTS W, HARMEGNIES A, et al. Advances in Multi-GNSS Time Transfer[C]∥2013 Joint European Frequency and Time Forum & International Frequency Control Symposium (EFTF/IFC). Prague, Czech Republic: IEEE, 2013: 508-512.

    [3] DEFRAIGNE P, AERTS W, POTTIAUX E. Monitoring of UTC(k)’s Using PPP and IGS Real-time Products[J]. GPS Solutions, 2015, 19(1): 165-172.

    [4] 張小紅, 蔡詩(shī)響, 李星星, 等. 利用GPS精密單點(diǎn)定位進(jìn)行時(shí)間傳遞精度分析[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2010, 35(3): 274-278.

    ZHANG Xiaohong, CAI Shixiang, LI Xingxing, et al. Accuracy Analysis of Time and Frequency Transfer Based on Precise Point Positioning[J]. Geomatics and Information Science of Wuhan University, 2010, 35(3): 274-278.

    [5] 黃觀文. GNSS星載原子鐘質(zhì)量評(píng)價(jià)及精密鐘差算法研究[D]. 西安: 長(zhǎng)安大學(xué), 2012.

    HUANG Guanwen. Research on Algorithms of Precise Clock Offset and Quality Evaluation of GNSS Satellite Clock[D]. Xi’an: Chang’an University, 2012.

    [6] 閆偉, 袁運(yùn)斌, 歐吉坤, 等. 非組合精密單點(diǎn)定位算法精密授時(shí)的可行性研究[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2011, 36(6): 648-651.

    YAN Wei, YUAN Yunbin, OU Jikun, et al. Feasibility of Precise Timing with Uncombined PPP[J]. Geomatics and Information Science of Wuhan University, 2011, 36(6): 648-651.

    [7] 張寶成, 歐吉坤, 袁運(yùn)斌, 等. 基于GPS雙頻原始觀測(cè)值的精密單點(diǎn)定位算法及應(yīng)用[J]. 測(cè)繪學(xué)報(bào), 2010, 39(5): 478-483.

    ZHANG Baocheng, OU Jikun, YUAN Yunbin, et al. Precise Point Positioning Algorithm Based on Original Dual-frequency GPS Code and Carrier-phase Observations and Its Application[J]. Acta Geodaetica et Cartographica Sinica, 2010, 39(5): 478-483.

    [8] 張小紅, 陳興漢, 郭斐. 高性能原子鐘鐘差建模及其在精密單點(diǎn)定位中的應(yīng)用[J]. 測(cè)繪學(xué)報(bào), 2015, 44(4): 392-398. DOI: 10.11947/j.AGCS.2015.20140287.ZHANG Xiaohong, CHEN Xinghan, GUO Fei. High-performance Atomic Clock Modeling and Its Application in Precise Point Positioning[J]. Acta Geodaetica et Cartographica Sinica, 2015, 44(4): 392-398. DOI: 10.11947/j.AGCS.2015.20140287.

    [9] WANG Kan, ROTHACHER M. Stochastic Modeling of High-stability Ground Clocks in GPS Analysis[J]. Journal of Geodesy, 2013, 87(5): 427-437.

    [10] LICHTEN S M, BORDER J S. Strategies for High-precision Global Positioning System Orbit Determination[J]. Journal of Geophysical Research: Solid Earth, 1987, 92(B12): 12751-12762.

    [11] JONES R H, TRYON P V. Continuous Time Series Models for Unequally Spaced Data Applied to Modeling Atomic Clocks[J]. SIAM Journal on Scientific and Statistical Computing, 1987, 8(1): 71-81.

    [12] 林旭, 羅志才. 一種新的衛(wèi)星鐘差Kalman濾波噪聲協(xié)方差估計(jì)方法[J]. 物理學(xué)報(bào), 2015, 64(8): 080201.

    LIN Xu, LUO Zhicai. A New Noise Covariance Matrix Estimation Method of Kalman Filter for Satellite Clock Errors[J]. Acta Physica Sinica, 2015, 64(8): 080201.

    [13] 郭海榮, 楊元喜, 何海波, 等. 導(dǎo)航衛(wèi)星原子鐘Kalman濾波中噪聲方差-協(xié)方差的確定[J]. 測(cè)繪學(xué)報(bào), 2010, 39(2): 146-150.

    GUO Hairong, YANG Yuanxi, HE Haibo, et al. Determination of Covariance Matrix of Kalman Filter Used for Time Prediction of Atomic Clocks of Navigation Satellites[J]. Acta Geodaetica et Cartographica Sinica, 2010, 39(2): 146-150.

    [14] SU W, FILLER R L. Application of Kalman Filtering Techniques to the Precision Clock with Non-constant Aging[C]∥Proceedings of the 46th Frequency Control Symposium. Hershey, PA: IEEE, 1992: 231-237.

    [15] STEIN S R, FILLER R L. Kalman Filter Analysis for Real Time Applications of Clocks and Oscillators[C]∥Proceedings of the 42nd Annual Frequency Control Symposium. Baltimore, MD: IEEE, 1988: 447-452.

    [16] 張清華, 隋立芬, 賈小林. 應(yīng)用Jones-Tryon Kalman濾波器對(duì)在軌GPS Rb鐘進(jìn)行狀態(tài)監(jiān)測(cè)[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2012, 37(4): 436-440.

    ZHANG Qinghua, SUI Lifen, JIA Xiaolin. Monitor State of GPS Rb Clock Using Jones-Tryon Kalman Filter[J]. Geomatics and Information Science of Wuhan University, 2012, 37(4): 436-440.

    [17] YANG Yang, YUE Xiaokui, YUAN Jianping, et al. Enhancing the Kinematic Precise Orbit Determination of Low Earth Orbiters Using GPS Receiver Clock Modelling[J]. Advances in Space Research, 2014, 54(9): 1901-1912.

    [18] WEINBACH U, SCH?N S. Improved GRACE Kinematic Orbit Determination Using GPS Receiver Clock Modeling[J]. GPS Solutions, 2013, 17(4): 511-520.

    [19] WEINBACH U, SCHON S. Improved GPS Receiver Clock Modeling for Kinematic Orbit Determination of the GRACE Satellites[C]∥Proceedings of European Frequency and Time Forum. Gothenburg, Sweden: IEEE, 2012: 157-160.

    [20] 肖國(guó)銳, 隋立芬, 陳泉余, 等. 利用接收機(jī)鐘差建模提升PPP收斂速度及精度[J]. 測(cè)繪科學(xué)技術(shù)學(xué)報(bào), 2015, 32(6): 555-558, 564.XIAO Guorui, SUI Lifen, CHEN Quanyu, et al. Improving PPP Convergence and Accuracy Using Receiver Clock Modeling[J]. Journal of Geomatics Science and Technology, 2015, 32(6): 555-558, 564

    [21] WEINBACH U, SCH?N S. GNSS Receiver Clock Modeling When Using High-precision Oscillators and Its Impact on PPP[J]. Advances in Space Research, 2011, 47(2): 229-238.

    [22] FILHO E A M, KUGA H K, LOPES R V F. Real Time Estimation of GPS Receiver Clock Offset by the Kalman Filter[J]. Personal Communication, 2003.

    [23] 于合理, 郝金明, 劉偉平, 等. 一種衛(wèi)星鐘差異常實(shí)時(shí)監(jiān)測(cè)算法[J]. 武漢大學(xué)學(xué)報(bào)(信息科學(xué)版), 2016, 41(1): 106-110.

    YU Heli, HAO Jinming, LIU Weiping, et al. A Real-time Anomaly Monitoring Algorithm for Satellite Clock[J]. Geomatics and Information Science of Wuhan University, 2016, 41(1): 106-110.

    [24] RILEY W J. Handbook of Frequency Stability Analysis[R]. NIST Special Publication 1065, 2007: 1-123.

    [25] MANNING D M. AF/NGA GPS Monitor Station High-performance Cesium Frequency Standard Stability 2005/2006: from NGA Kalman Filter Clock Estimates[C]∥Proceedings of the 38th Annual Precise Time and Time Interval Systems and Applications Meeting. Reston, Virginia: [s.n.], 2006: 137-152.

    [26] CLKLOG. A Summary File of the Deployment History for GPS Receiver, Antenna, Frequency Standards, and Other Equipment at IGS Stations[EB/OL]. (2010-01-20) [2016-06-02]. ftp:∥igscb.jpl.nasa.gov/igscb/station/general/loghist.txt.

    [27] 于合理, 郝金明, 謝建濤, 等. 硬件延遲偏差對(duì)衛(wèi)星鐘差解算的影響[J]. 導(dǎo)航定位學(xué)報(bào), 2015, 3(1): 71-73.

    YU Heli, HAO Jinming, XIE Jiantao, et al. Impact of DCB on the Estimate of Satellite Clock Bias[J]. Journal of Navigation and Positioning, 2015, 3(1): 71-73.

    (責(zé)任編輯:叢樹平)

    A Time Transfer Algorithm of Precise Point Positioning with Additional Atomic Clock Physical Model

    YU Heli1,2,HAO Jinming1,2,LIU Weiping1,2,TIAN Yingguo1,2,DENG Ke1

    1. School of Navigation and Aerospace Engineering, Information Engineering University, Zhengzhou 450001, China; 2. BeiDou Navigation Technology Collaborative Innovation Center of Henan, Zhengzhou 450001, China

    In the traditional time transfer algorithms of precise point positioning, the receiver clock offset is estimated epoch by epoch as independent white noise and the short-term correlation between the epoch of the atomic clock offset is ignored. In order to solve this problem, a time transfer algorithm of precise point positioning with additional atomic clock physical model is proposed. The receiver clock parameters are expanded in the traditional model of precise point positioning by using Kalman filter to model the clock errors of high stability atomic clock, and the method of determining the covariance matrices of process noise and the initial state vector of the Kalman filter is presented. Experimental results show that the proposed algorithm can effectively avoid the convergence process of the traditional precise point positioning algorithm, and the result of the algorithm is more consistent with the physical properties of the atomic clock, this algorithm can significantly improve the accuracy and stability of the time transfer results, the accuracy of time transfer results of single station can be increased on average by 58%, and the accuracy of time synchronization between stations can be improved on average by 51%.

    atomic clock; physical model; precise point positioning; time transfer

    Specific Project of Chinese Second-generation Satellite Navigation System (No. GFZX0301040308); The Open Research Foundation of State Key Laboratory of Geo-information Engineering (No. SKLGIE2015-M-1-6)

    YU Heli(1989—),male,PhD candidate,majors in GNSS time and frequency transfer.

    于合理,郝金明,劉偉平,等.附加原子鐘物理模型的PPP時(shí)間傳遞算法[J].測(cè)繪學(xué)報(bào),2016,45(11):1285-1292.

    10.11947/j.AGCS.2016.20160217.

    YU Heli,HAO Jinming,LIU Weiping,et al.A Time Transfer Algorithm of Precise Point Positioning with Additional Atomic Clock Physical Model[J]. Acta Geodaetica et Cartographica Sinica,2016,45(11):1285-1292. DOI:10.11947/j.AGCS.2016.20160217.

    P228

    A

    1001-1595(2016)11-1285-08

    中國(guó)第二代衛(wèi)星導(dǎo)航系統(tǒng)重大專項(xiàng)(GFZX0301040308);地理信息工程國(guó)家重點(diǎn)實(shí)驗(yàn)室開放研究基金(SKLGIE2015-M-1-6)

    2016-05-05

    修回日期: 2016-09-05

    于合理(1989—),男,博士生,研究方向?yàn)镚NSS時(shí)頻傳遞方法。

    E-mail: yuhl1989@163.com

    猜你喜歡
    原子鐘鐘差協(xié)方差
    超高精度計(jì)時(shí)器——原子鐘
    IGS快速/超快速衛(wèi)星鐘差精度評(píng)定與分析
    用于小型銣如原子鐘中介質(zhì)諧振腔激勵(lì)分析
    實(shí)時(shí)干涉測(cè)量中對(duì)流層延遲與鐘差精修正建模
    載人航天(2016年4期)2016-12-01 06:56:24
    基于拉格朗日的IGS精密星歷和鐘差插值分析
    原子鐘頻跳快速探測(cè)方法
    不確定系統(tǒng)改進(jìn)的魯棒協(xié)方差交叉融合穩(wěn)態(tài)Kalman預(yù)報(bào)器
    一種基于廣義協(xié)方差矩陣的欠定盲辨識(shí)方法
    北斗星載原子鐘異常數(shù)據(jù)處理的研究
    基于區(qū)域CORS的實(shí)時(shí)精密衛(wèi)星鐘差估計(jì)研究
    欧美日韩亚洲国产一区二区在线观看| 久久这里只有精品19| 国产亚洲欧美精品永久| 欧美精品啪啪一区二区三区| 丁香欧美五月| 性色av乱码一区二区三区2| av福利片在线| 啦啦啦观看免费观看视频高清| 亚洲色图av天堂| 无遮挡黄片免费观看| 90打野战视频偷拍视频| 黄色片一级片一级黄色片| 国产aⅴ精品一区二区三区波| 亚洲最大成人中文| 一a级毛片在线观看| 亚洲av电影在线进入| 别揉我奶头~嗯~啊~动态视频| 久久午夜综合久久蜜桃| 婷婷精品国产亚洲av| 国产精品久久久av美女十八| 国产高清有码在线观看视频 | 精品久久久久久久末码| 少妇的丰满在线观看| 国产亚洲精品综合一区在线观看 | 变态另类丝袜制服| 制服诱惑二区| 精品一区二区三区av网在线观看| 亚洲 欧美 日韩 在线 免费| 天堂动漫精品| 精品熟女少妇八av免费久了| 久久久久久久久免费视频了| 丝袜在线中文字幕| 亚洲av日韩精品久久久久久密| 俺也久久电影网| 久久久国产成人精品二区| 人人妻人人看人人澡| 久久久久久九九精品二区国产 | 狠狠狠狠99中文字幕| 夜夜夜夜夜久久久久| 99国产精品99久久久久| 国产三级在线视频| 亚洲精品中文字幕一二三四区| 国产日本99.免费观看| 丁香六月欧美| 少妇粗大呻吟视频| 美女大奶头视频| 美女免费视频网站| 亚洲男人天堂网一区| 长腿黑丝高跟| а√天堂www在线а√下载| 一区二区三区精品91| 亚洲精品国产一区二区精华液| 99久久99久久久精品蜜桃| 亚洲中文字幕日韩| 久久狼人影院| or卡值多少钱| 久久久久九九精品影院| 久久久精品国产亚洲av高清涩受| 午夜成年电影在线免费观看| 香蕉久久夜色| 一进一出抽搐gif免费好疼| 亚洲精品久久成人aⅴ小说| 变态另类丝袜制服| 性色av乱码一区二区三区2| 好男人在线观看高清免费视频 | 欧美一区二区精品小视频在线| 欧美一区二区精品小视频在线| 国产私拍福利视频在线观看| 午夜免费激情av| 高清在线国产一区| 人人妻人人看人人澡| 日本免费一区二区三区高清不卡| 色综合站精品国产| 久久久精品国产亚洲av高清涩受| 又黄又爽又免费观看的视频| 久久99热这里只有精品18| 亚洲精华国产精华精| 国产精品永久免费网站| 一区二区三区激情视频| 国产欧美日韩精品亚洲av| 久久香蕉激情| 中文字幕精品免费在线观看视频| 男女午夜视频在线观看| 亚洲黑人精品在线| 黑人欧美特级aaaaaa片| 91字幕亚洲| 成人亚洲精品一区在线观看| 欧美黄色淫秽网站| 777久久人妻少妇嫩草av网站| 久热这里只有精品99| 亚洲精华国产精华精| 国产精品永久免费网站| 久久精品91蜜桃| 一二三四社区在线视频社区8| 欧美成狂野欧美在线观看| 脱女人内裤的视频| 一级a爱片免费观看的视频| 国产高清有码在线观看视频 | 成人欧美大片| 亚洲国产精品999在线| 国产高清有码在线观看视频 | 色综合欧美亚洲国产小说| 超碰成人久久| 亚洲男人天堂网一区| 久久香蕉激情| 女人爽到高潮嗷嗷叫在线视频| 动漫黄色视频在线观看| 男女下面进入的视频免费午夜 | 欧美 亚洲 国产 日韩一| 国产在线观看jvid| 在线观看免费日韩欧美大片| 亚洲av成人一区二区三| 国产主播在线观看一区二区| 亚洲人成电影免费在线| 日本熟妇午夜| 国产精品二区激情视频| 亚洲免费av在线视频| 欧美精品啪啪一区二区三区| 亚洲熟妇中文字幕五十中出| 在线av久久热| 在线观看www视频免费| 国产黄a三级三级三级人| x7x7x7水蜜桃| 欧美日韩精品网址| av有码第一页| 午夜福利成人在线免费观看| 国产亚洲av高清不卡| 长腿黑丝高跟| 亚洲av电影不卡..在线观看| 国产亚洲av高清不卡| 国产黄a三级三级三级人| xxxwww97欧美| 国产精品亚洲一级av第二区| 精品久久久久久久毛片微露脸| 亚洲成国产人片在线观看| 老熟妇乱子伦视频在线观看| 99在线人妻在线中文字幕| 久久久精品国产亚洲av高清涩受| 又紧又爽又黄一区二区| 久久中文字幕人妻熟女| 久久人妻av系列| 自线自在国产av| 久久人妻福利社区极品人妻图片| 最新美女视频免费是黄的| 成人av一区二区三区在线看| 天堂√8在线中文| 在线观看日韩欧美| 757午夜福利合集在线观看| 国产精品久久久久久精品电影 | 无遮挡黄片免费观看| 一级片免费观看大全| 精品第一国产精品| 丁香欧美五月| 在线观看免费午夜福利视频| 久久伊人香网站| 女同久久另类99精品国产91| 国产一区二区在线av高清观看| 在线看三级毛片| 国产1区2区3区精品| www日本在线高清视频| 久久午夜综合久久蜜桃| 女性被躁到高潮视频| 自线自在国产av| 一进一出好大好爽视频| 丝袜人妻中文字幕| 级片在线观看| 老司机在亚洲福利影院| 老熟妇仑乱视频hdxx| 在线免费观看的www视频| 母亲3免费完整高清在线观看| 国产97色在线日韩免费| 欧美乱色亚洲激情| 人人妻,人人澡人人爽秒播| 亚洲九九香蕉| 亚洲免费av在线视频| www日本在线高清视频| 日韩欧美一区二区三区在线观看| 嫩草影院精品99| 首页视频小说图片口味搜索| 香蕉久久夜色| 搡老岳熟女国产| 两个人视频免费观看高清| 亚洲欧美精品综合久久99| 女性被躁到高潮视频| av欧美777| 麻豆av在线久日| 成人18禁高潮啪啪吃奶动态图| 一级a爱视频在线免费观看| 日韩欧美一区二区三区在线观看| 午夜福利在线观看吧| www国产在线视频色| 欧美最黄视频在线播放免费| 日韩国内少妇激情av| 麻豆久久精品国产亚洲av| 亚洲国产毛片av蜜桃av| 久久久国产精品麻豆| 亚洲一区高清亚洲精品| 亚洲国产欧美网| 色哟哟哟哟哟哟| 18禁国产床啪视频网站| 青草久久国产| 男男h啪啪无遮挡| 大型黄色视频在线免费观看| 亚洲中文av在线| e午夜精品久久久久久久| 国产成+人综合+亚洲专区| 国产成人系列免费观看| 欧美性猛交╳xxx乱大交人| 真人一进一出gif抽搐免费| 精品日产1卡2卡| 亚洲午夜精品一区,二区,三区| 午夜影院日韩av| 国产精品一区二区三区四区久久 | 亚洲五月婷婷丁香| 亚洲在线自拍视频| 欧美精品亚洲一区二区| 亚洲激情在线av| 日本 欧美在线| 观看免费一级毛片| 日本黄色视频三级网站网址| 国产精品九九99| 亚洲成a人片在线一区二区| 成人欧美大片| 在线观看66精品国产| 美女高潮喷水抽搐中文字幕| 桃红色精品国产亚洲av| 看免费av毛片| 国产成人欧美| 国产精品日韩av在线免费观看| 热99re8久久精品国产| 窝窝影院91人妻| 国产成人系列免费观看| 久久天躁狠狠躁夜夜2o2o| 亚洲在线自拍视频| 自线自在国产av| 精品国产乱码久久久久久男人| 91麻豆av在线| 又黄又爽又免费观看的视频| 亚洲国产看品久久| 两个人视频免费观看高清| 午夜福利成人在线免费观看| 亚洲av五月六月丁香网| 国内揄拍国产精品人妻在线 | 99热这里只有精品一区 | av视频在线观看入口| 老司机福利观看| 亚洲专区国产一区二区| 亚洲激情在线av| 久久九九热精品免费| 国产精品野战在线观看| 国产亚洲精品第一综合不卡| 色播亚洲综合网| www.熟女人妻精品国产| 在线观看www视频免费| 哪里可以看免费的av片| 在线播放国产精品三级| www日本黄色视频网| 黄色片一级片一级黄色片| 黑人巨大精品欧美一区二区mp4| 中文字幕另类日韩欧美亚洲嫩草| 亚洲在线自拍视频| 国产蜜桃级精品一区二区三区| 日本撒尿小便嘘嘘汇集6| 精品国产亚洲在线| 日韩高清综合在线| 白带黄色成豆腐渣| 国产精品精品国产色婷婷| 看免费av毛片| 99精品欧美一区二区三区四区| 老司机福利观看| 亚洲精品在线观看二区| 在线观看免费午夜福利视频| 色哟哟哟哟哟哟| 国产精品久久视频播放| 国产乱人伦免费视频| 国产精品自产拍在线观看55亚洲| 中文字幕高清在线视频| 老汉色av国产亚洲站长工具| 国产精品永久免费网站| 一区二区日韩欧美中文字幕| www国产在线视频色| 嫩草影视91久久| 在线观看免费日韩欧美大片| 国产视频一区二区在线看| 最近最新中文字幕大全免费视频| 国产伦一二天堂av在线观看| 亚洲av电影在线进入| 看免费av毛片| 长腿黑丝高跟| 午夜免费鲁丝| 亚洲久久久国产精品| 成人免费观看视频高清| 精品一区二区三区av网在线观看| 欧美激情 高清一区二区三区| 97超级碰碰碰精品色视频在线观看| 婷婷精品国产亚洲av| 亚洲性夜色夜夜综合| 夜夜躁狠狠躁天天躁| www日本黄色视频网| 日韩一卡2卡3卡4卡2021年| 国产99白浆流出| 免费人成视频x8x8入口观看| 淫妇啪啪啪对白视频| 日本免费一区二区三区高清不卡| or卡值多少钱| aaaaa片日本免费| 亚洲专区字幕在线| 欧美黑人精品巨大| 成人亚洲精品av一区二区| 午夜影院日韩av| 美女 人体艺术 gogo| 亚洲一区二区三区色噜噜| 麻豆av在线久日| 曰老女人黄片| 久久亚洲真实| 午夜激情av网站| 天天添夜夜摸| 精品卡一卡二卡四卡免费| 国产黄a三级三级三级人| 国产在线精品亚洲第一网站| 宅男免费午夜| 精品第一国产精品| 99精品欧美一区二区三区四区| 黄色丝袜av网址大全| 国产又色又爽无遮挡免费看| 亚洲男人的天堂狠狠| 一级a爱片免费观看的视频| 非洲黑人性xxxx精品又粗又长| 身体一侧抽搐| 在线观看午夜福利视频| 午夜日韩欧美国产| 亚洲熟妇熟女久久| 亚洲 欧美一区二区三区| 亚洲欧美日韩无卡精品| 欧美色视频一区免费| 亚洲 欧美 日韩 在线 免费| 一本一本综合久久| 国产高清视频在线播放一区| 1024手机看黄色片| 国内少妇人妻偷人精品xxx网站 | 禁无遮挡网站| 色哟哟哟哟哟哟| 亚洲成av人片免费观看| 亚洲精品国产区一区二| 精品欧美国产一区二区三| 亚洲人成网站在线播放欧美日韩| 国产精华一区二区三区| 妹子高潮喷水视频| tocl精华| 每晚都被弄得嗷嗷叫到高潮| av福利片在线| 日韩精品青青久久久久久| 18禁黄网站禁片午夜丰满| 日本免费一区二区三区高清不卡| 99热6这里只有精品| 亚洲精品一卡2卡三卡4卡5卡| 日本 av在线| 天天躁狠狠躁夜夜躁狠狠躁| 丝袜美腿诱惑在线| 色尼玛亚洲综合影院| 悠悠久久av| 最近最新中文字幕大全免费视频| 日韩精品中文字幕看吧| 国产成人精品久久二区二区91| 国产精品乱码一区二三区的特点| 日韩av在线大香蕉| 日本成人三级电影网站| 久久性视频一级片| 老司机靠b影院| 制服诱惑二区| 天天躁夜夜躁狠狠躁躁| 午夜激情福利司机影院| 国产男靠女视频免费网站| 国产视频内射| 亚洲国产毛片av蜜桃av| 国产精品久久电影中文字幕| 国产精品久久久久久精品电影 | 午夜a级毛片| 国产亚洲欧美精品永久| 国产精品久久久久久精品电影 | or卡值多少钱| 免费在线观看视频国产中文字幕亚洲| 狠狠狠狠99中文字幕| 美女国产高潮福利片在线看| 女性生殖器流出的白浆| 国产亚洲欧美精品永久| 国产免费av片在线观看野外av| 伦理电影免费视频| 亚洲熟妇中文字幕五十中出| 色综合婷婷激情| 黑丝袜美女国产一区| 国产亚洲精品av在线| 国产精品精品国产色婷婷| 欧美乱色亚洲激情| 夜夜爽天天搞| 国产午夜福利久久久久久| 99在线视频只有这里精品首页| 老熟妇仑乱视频hdxx| 午夜老司机福利片| 国产精品久久久av美女十八| 国产精品 欧美亚洲| 免费av毛片视频| 最近最新中文字幕大全电影3 | 成熟少妇高潮喷水视频| 久久精品国产清高在天天线| 50天的宝宝边吃奶边哭怎么回事| 久久香蕉国产精品| 亚洲一区二区三区不卡视频| 色播亚洲综合网| 少妇 在线观看| 禁无遮挡网站| 嫩草影视91久久| 夜夜夜夜夜久久久久| 亚洲自拍偷在线| 在线av久久热| 午夜福利一区二区在线看| 九色国产91popny在线| АⅤ资源中文在线天堂| 91九色精品人成在线观看| 久久中文字幕人妻熟女| 国产黄片美女视频| 国产av又大| 欧美三级亚洲精品| 中文字幕精品亚洲无线码一区 | 亚洲人成电影免费在线| 亚洲美女黄片视频| 色尼玛亚洲综合影院| 岛国视频午夜一区免费看| 窝窝影院91人妻| 午夜福利免费观看在线| 精品高清国产在线一区| 叶爱在线成人免费视频播放| 欧美日本视频| 人成视频在线观看免费观看| 精品一区二区三区视频在线观看免费| 一级作爱视频免费观看| 久久国产乱子伦精品免费另类| 99热这里只有精品一区 | 久久婷婷成人综合色麻豆| 国产成人av教育| 亚洲欧美日韩无卡精品| 亚洲欧洲精品一区二区精品久久久| 国产精品一区二区三区四区久久 | 亚洲久久久国产精品| 真人做人爱边吃奶动态| 亚洲欧美一区二区三区黑人| 亚洲熟女毛片儿| 一级作爱视频免费观看| 亚洲国产毛片av蜜桃av| 制服人妻中文乱码| 国产精品九九99| 国产伦一二天堂av在线观看| 欧美成人一区二区免费高清观看 | 午夜激情av网站| 给我免费播放毛片高清在线观看| 亚洲全国av大片| 免费在线观看成人毛片| 国产一区在线观看成人免费| 亚洲天堂国产精品一区在线| 90打野战视频偷拍视频| 白带黄色成豆腐渣| 久久久国产成人免费| 午夜影院日韩av| 色尼玛亚洲综合影院| 老汉色av国产亚洲站长工具| 国产高清激情床上av| 国产成人一区二区三区免费视频网站| www日本黄色视频网| 国产aⅴ精品一区二区三区波| 欧美激情极品国产一区二区三区| 欧美绝顶高潮抽搐喷水| 中国美女看黄片| 亚洲国产欧美日韩在线播放| 国产97色在线日韩免费| 99re在线观看精品视频| 精品久久蜜臀av无| 亚洲第一av免费看| 亚洲国产欧美一区二区综合| 亚洲免费av在线视频| 每晚都被弄得嗷嗷叫到高潮| 在线观看午夜福利视频| 国产精品一区二区精品视频观看| 满18在线观看网站| 搡老熟女国产l中国老女人| 亚洲av成人av| 亚洲在线自拍视频| 99热只有精品国产| 级片在线观看| 久久天堂一区二区三区四区| 欧美国产日韩亚洲一区| 美女午夜性视频免费| www国产在线视频色| 日韩成人在线观看一区二区三区| 两人在一起打扑克的视频| 亚洲真实伦在线观看| 欧美色欧美亚洲另类二区| 欧美+亚洲+日韩+国产| 国产国语露脸激情在线看| 久久精品人妻少妇| 99久久国产精品久久久| 大型黄色视频在线免费观看| 亚洲人成电影免费在线| 自线自在国产av| 国内精品久久久久久久电影| 久久99热这里只有精品18| 丁香六月欧美| 人成视频在线观看免费观看| 国产高清视频在线播放一区| 色av中文字幕| 黄色视频不卡| 老熟妇乱子伦视频在线观看| 免费在线观看黄色视频的| 国产又色又爽无遮挡免费看| 国产激情欧美一区二区| 69av精品久久久久久| 亚洲国产欧美日韩在线播放| 精品午夜福利视频在线观看一区| 亚洲精品美女久久久久99蜜臀| 中文字幕精品免费在线观看视频| 久久九九热精品免费| 色在线成人网| 黄色a级毛片大全视频| 99国产综合亚洲精品| 熟女电影av网| 日本撒尿小便嘘嘘汇集6| 丰满的人妻完整版| 欧美日本视频| 大香蕉久久成人网| 色在线成人网| 夜夜躁狠狠躁天天躁| 黄频高清免费视频| 欧美亚洲日本最大视频资源| 特大巨黑吊av在线直播 | 观看免费一级毛片| 精品国产乱码久久久久久男人| 成人特级黄色片久久久久久久| 成年免费大片在线观看| 1024香蕉在线观看| 国产亚洲精品av在线| 18禁裸乳无遮挡免费网站照片 | 久9热在线精品视频| 亚洲国产精品999在线| 我的亚洲天堂| 大型av网站在线播放| 18禁观看日本| 99热这里只有精品一区 | 久久久国产欧美日韩av| 国产成人av激情在线播放| 18禁美女被吸乳视频| 悠悠久久av| 午夜久久久在线观看| 成人亚洲精品一区在线观看| 国产视频一区二区在线看| 国产区一区二久久| 亚洲精品美女久久久久99蜜臀| 欧美激情久久久久久爽电影| 国产国语露脸激情在线看| 免费av毛片视频| 国产熟女xx| 黄片大片在线免费观看| 国产精品久久久av美女十八| 成人欧美大片| 免费看十八禁软件| 亚洲自拍偷在线| 中文字幕久久专区| 一级片免费观看大全| 亚洲精品在线观看二区| 黑人欧美特级aaaaaa片| 最近最新中文字幕大全免费视频| 亚洲全国av大片| 中文字幕最新亚洲高清| 叶爱在线成人免费视频播放| 欧美亚洲日本最大视频资源| 国产成人啪精品午夜网站| 亚洲狠狠婷婷综合久久图片| 国产午夜福利久久久久久| 精品日产1卡2卡| 久久久国产精品麻豆| 亚洲人成网站在线播放欧美日韩| 少妇的丰满在线观看| 亚洲精品粉嫩美女一区| 一级a爱视频在线免费观看| 国产区一区二久久| 久久精品影院6| av免费在线观看网站| 欧美色欧美亚洲另类二区| 国产精品一区二区三区四区久久 | 少妇裸体淫交视频免费看高清 | 国产不卡一卡二| 99国产精品99久久久久| 男人的好看免费观看在线视频 | 老鸭窝网址在线观看| 女人高潮潮喷娇喘18禁视频| 黄片播放在线免费| cao死你这个sao货| av天堂在线播放| 人人澡人人妻人| 免费搜索国产男女视频| 伊人久久大香线蕉亚洲五| 久久人妻av系列| 久久久久久久精品吃奶| 婷婷精品国产亚洲av在线| 午夜免费鲁丝| 久久天躁狠狠躁夜夜2o2o| 日韩成人在线观看一区二区三区| 成人欧美大片| 黄片大片在线免费观看| av欧美777| 亚洲一区二区三区不卡视频| 欧美日韩一级在线毛片| 国产一卡二卡三卡精品| 欧美 亚洲 国产 日韩一| 人人妻人人澡人人看| 欧美黄色片欧美黄色片| 国语自产精品视频在线第100页|