• <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ì)研究
    欧美黑人精品巨大| 成人手机av| 叶爱在线成人免费视频播放| 黄片小视频在线播放| 新久久久久国产一级毛片| 亚洲五月色婷婷综合| 欧美不卡视频在线免费观看 | 他把我摸到了高潮在线观看| 欧美亚洲日本最大视频资源| 亚洲av成人一区二区三| 久久久久精品国产欧美久久久| 国产成人免费观看mmmm| 无人区码免费观看不卡| 欧美成狂野欧美在线观看| 亚洲九九香蕉| 五月开心婷婷网| 19禁男女啪啪无遮挡网站| 91麻豆精品激情在线观看国产 | 日韩免费av在线播放| 在线观看66精品国产| 一二三四在线观看免费中文在| 曰老女人黄片| 中文字幕精品免费在线观看视频| 91九色精品人成在线观看| 亚洲国产欧美网| 丝袜美足系列| 脱女人内裤的视频| 精品第一国产精品| 国产aⅴ精品一区二区三区波| av一本久久久久| 国产成人啪精品午夜网站| 国产精品 欧美亚洲| 免费久久久久久久精品成人欧美视频| 99热网站在线观看| x7x7x7水蜜桃| 国产国语露脸激情在线看| 午夜老司机福利片| 制服诱惑二区| tocl精华| 午夜精品在线福利| 大香蕉久久成人网| 国产熟女午夜一区二区三区| 制服诱惑二区| 女警被强在线播放| 两人在一起打扑克的视频| 精品国产一区二区久久| 久久久国产成人免费| 久久午夜亚洲精品久久| 国产av又大| 久久香蕉激情| 国产精品美女特级片免费视频播放器 | 成人永久免费在线观看视频| 王馨瑶露胸无遮挡在线观看| 久久 成人 亚洲| 国产激情欧美一区二区| 久久青草综合色| 免费观看精品视频网站| av国产精品久久久久影院| 免费观看精品视频网站| 久久国产精品人妻蜜桃| 国产在线精品亚洲第一网站| 国产99久久九九免费精品| 免费观看精品视频网站| 亚洲国产欧美网| 免费日韩欧美在线观看| 国产无遮挡羞羞视频在线观看| 欧美日韩亚洲高清精品| 亚洲精品国产区一区二| 久久婷婷成人综合色麻豆| 午夜精品国产一区二区电影| 亚洲专区中文字幕在线| 国产精品影院久久| 国产精品影院久久| 日韩欧美国产一区二区入口| av天堂在线播放| videos熟女内射| 亚洲欧美色中文字幕在线| 国产真人三级小视频在线观看| 国产又色又爽无遮挡免费看| 国产真人三级小视频在线观看| 久久亚洲精品不卡| 国产高清激情床上av| 宅男免费午夜| 五月开心婷婷网| 午夜福利,免费看| a级片在线免费高清观看视频| 久久中文看片网| 美女高潮喷水抽搐中文字幕| 日本精品一区二区三区蜜桃| 亚洲精品粉嫩美女一区| 国产精品一区二区在线观看99| 久久人妻福利社区极品人妻图片| 亚洲七黄色美女视频| 久久久久久久精品吃奶| 男人的好看免费观看在线视频 | 美女国产高潮福利片在线看| 欧美精品啪啪一区二区三区| 一边摸一边做爽爽视频免费| 国产视频一区二区在线看| 一区二区三区激情视频| 国产欧美日韩一区二区三区在线| 99国产极品粉嫩在线观看| 99久久人妻综合| 婷婷精品国产亚洲av在线 | 女警被强在线播放| e午夜精品久久久久久久| 亚洲第一欧美日韩一区二区三区| 黄色丝袜av网址大全| 久久热在线av| 欧美另类亚洲清纯唯美| 中文字幕最新亚洲高清| 日日摸夜夜添夜夜添小说| 亚洲一区二区三区不卡视频| 正在播放国产对白刺激| 国产精品亚洲av一区麻豆| 多毛熟女@视频| 欧美久久黑人一区二区| 中文字幕高清在线视频| 两性午夜刺激爽爽歪歪视频在线观看 | 18禁观看日本| 国产真人三级小视频在线观看| 人妻 亚洲 视频| 在线天堂中文资源库| 日韩视频一区二区在线观看| 一夜夜www| 欧美精品人与动牲交sv欧美| 在线观看免费日韩欧美大片| 国产欧美日韩综合在线一区二区| 日本撒尿小便嘘嘘汇集6| 午夜精品在线福利| 99精品欧美一区二区三区四区| 99精品在免费线老司机午夜| 国产免费现黄频在线看| 欧洲精品卡2卡3卡4卡5卡区| 新久久久久国产一级毛片| 免费黄频网站在线观看国产| 久久天堂一区二区三区四区| 国产精品.久久久| 免费观看人在逋| 高清毛片免费观看视频网站 | 日韩 欧美 亚洲 中文字幕| 欧美老熟妇乱子伦牲交| 在线天堂中文资源库| 露出奶头的视频| 亚洲欧美一区二区三区黑人| av超薄肉色丝袜交足视频| 久久香蕉激情| 精品无人区乱码1区二区| 成人三级做爰电影| 正在播放国产对白刺激| 国产99久久九九免费精品| 一本综合久久免费| 国产精品99久久99久久久不卡| 亚洲五月色婷婷综合| 中文字幕高清在线视频| 精品一区二区三区av网在线观看| 在线观看舔阴道视频| 少妇 在线观看| 男人舔女人的私密视频| 日本欧美视频一区| 亚洲成av片中文字幕在线观看| 欧美日韩亚洲综合一区二区三区_| 国产黄色免费在线视频| 波多野结衣一区麻豆| 女人被狂操c到高潮| 欧洲精品卡2卡3卡4卡5卡区| 窝窝影院91人妻| 国产精品久久久久成人av| 亚洲熟女毛片儿| 人人妻,人人澡人人爽秒播| 国产主播在线观看一区二区| 成年人免费黄色播放视频| 男女下面插进去视频免费观看| 免费观看a级毛片全部| 久久精品国产99精品国产亚洲性色 | 日本wwww免费看| 久久热在线av| 欧美乱码精品一区二区三区| av免费在线观看网站| 一进一出抽搐gif免费好疼 | 在线观看一区二区三区激情| 中文字幕人妻丝袜一区二区| 日韩视频一区二区在线观看| 视频在线观看一区二区三区| 国产精品秋霞免费鲁丝片| 国产精品 欧美亚洲| 午夜福利影视在线免费观看| 国产亚洲一区二区精品| 新久久久久国产一级毛片| 国产精品免费一区二区三区在线 | 后天国语完整版免费观看| 欧美黄色淫秽网站| 欧美久久黑人一区二区| 又黄又爽又免费观看的视频| 女人被躁到高潮嗷嗷叫费观| 在线观看免费视频网站a站| 免费高清在线观看日韩| 18禁裸乳无遮挡免费网站照片 | 午夜福利在线免费观看网站| 久久精品91无色码中文字幕| 国产成人av教育| 久久午夜综合久久蜜桃| a级毛片黄视频| 亚洲美女黄片视频| 男女床上黄色一级片免费看| 久久久久久久久久久久大奶| 日本wwww免费看| 淫妇啪啪啪对白视频| 色老头精品视频在线观看| 亚洲全国av大片| 久久精品国产清高在天天线| 性色av乱码一区二区三区2| 18禁裸乳无遮挡免费网站照片 | 电影成人av| 日本vs欧美在线观看视频| 精品第一国产精品| a级毛片黄视频| 亚洲精品在线观看二区| 1024香蕉在线观看| 欧美日韩乱码在线| 亚洲一卡2卡3卡4卡5卡精品中文| 国产男女内射视频| 91在线观看av| 国产高清视频在线播放一区| 热99久久久久精品小说推荐| 午夜亚洲福利在线播放| 国产精品美女特级片免费视频播放器 | 国产有黄有色有爽视频| 久久性视频一级片| 久久久久视频综合| 精品国产乱子伦一区二区三区| 亚洲色图av天堂| 新久久久久国产一级毛片| 亚洲精品美女久久av网站| 国产熟女午夜一区二区三区| 国产区一区二久久| 亚洲欧洲精品一区二区精品久久久| 亚洲av成人一区二区三| 丰满饥渴人妻一区二区三| 一区福利在线观看| 亚洲av熟女| 欧美黄色淫秽网站| 在线国产一区二区在线| 激情视频va一区二区三区| 亚洲国产精品sss在线观看 | 精品人妻1区二区| 免费在线观看完整版高清| 亚洲一码二码三码区别大吗| 亚洲人成伊人成综合网2020| 久久人人爽av亚洲精品天堂| 久热爱精品视频在线9| 国产一区有黄有色的免费视频| 亚洲中文av在线| 欧美人与性动交α欧美精品济南到| 国产精品一区二区免费欧美| 黄色a级毛片大全视频| 一边摸一边抽搐一进一出视频| 免费看a级黄色片| 国产免费现黄频在线看| 国产麻豆69| 亚洲国产欧美网| 久久精品国产a三级三级三级| 亚洲熟妇熟女久久| 国产精品偷伦视频观看了| 久久草成人影院| 黄色毛片三级朝国网站| 999久久久精品免费观看国产| 国产蜜桃级精品一区二区三区 | 99精品在免费线老司机午夜| 国产99白浆流出| 黄色女人牲交| 国产精品久久久久久人妻精品电影| 三级毛片av免费| 丁香六月欧美| 悠悠久久av| 人人妻人人爽人人添夜夜欢视频| 精品福利永久在线观看| 久久久久久亚洲精品国产蜜桃av| 王馨瑶露胸无遮挡在线观看| 伊人久久大香线蕉亚洲五| 国产午夜精品久久久久久| 青草久久国产| 最近最新中文字幕大全电影3 | 一级毛片女人18水好多| 热re99久久精品国产66热6| 他把我摸到了高潮在线观看| 久久精品91无色码中文字幕| 色尼玛亚洲综合影院| 国产精品1区2区在线观看. | 国产欧美亚洲国产| 成年人黄色毛片网站| 国产精品一区二区在线不卡| 亚洲熟女毛片儿| 日韩欧美一区二区三区在线观看 | 少妇的丰满在线观看| 国产成人精品久久二区二区免费| 国产深夜福利视频在线观看| 国产亚洲欧美在线一区二区| 久久这里只有精品19| 亚洲欧美一区二区三区久久| 天天躁夜夜躁狠狠躁躁| 91成年电影在线观看| 日本vs欧美在线观看视频| 欧美乱码精品一区二区三区| 波多野结衣av一区二区av| 国内毛片毛片毛片毛片毛片| 亚洲第一欧美日韩一区二区三区| 欧美 亚洲 国产 日韩一| 在线视频色国产色| 久久精品国产亚洲av高清一级| 在线av久久热| 精品高清国产在线一区| 嫩草影视91久久| 美女高潮到喷水免费观看| 女警被强在线播放| 亚洲成人手机| 精品人妻1区二区| 99精国产麻豆久久婷婷| 日韩一卡2卡3卡4卡2021年| 久久精品国产综合久久久| 91成人精品电影| 在线看a的网站| 国产精品久久电影中文字幕 | 国产真人三级小视频在线观看| 久久人妻熟女aⅴ| 国产激情久久老熟女| 精品少妇一区二区三区视频日本电影| 久久天堂一区二区三区四区| 国产av又大| 看黄色毛片网站| 日韩成人在线观看一区二区三区| 高清黄色对白视频在线免费看| 国产精品久久电影中文字幕 | 国产蜜桃级精品一区二区三区 | 黄网站色视频无遮挡免费观看| 国产欧美日韩精品亚洲av| 热99re8久久精品国产| av电影中文网址| 成年动漫av网址| 日韩人妻精品一区2区三区| 丰满迷人的少妇在线观看| 悠悠久久av| 极品教师在线免费播放| 可以免费在线观看a视频的电影网站| 欧美不卡视频在线免费观看 | 色播在线永久视频| 波多野结衣av一区二区av| 交换朋友夫妻互换小说| 久久影院123| 一进一出抽搐gif免费好疼 | 日韩熟女老妇一区二区性免费视频| 少妇的丰满在线观看| 老司机亚洲免费影院| 亚洲视频免费观看视频| 大香蕉久久成人网| 一区二区三区国产精品乱码| 人人妻,人人澡人人爽秒播| 日本a在线网址| 国产激情久久老熟女| 欧美日韩黄片免| 国产欧美亚洲国产| 18禁裸乳无遮挡动漫免费视频| 一边摸一边抽搐一进一出视频| 国产精品免费一区二区三区在线 | xxxhd国产人妻xxx| 国产亚洲精品第一综合不卡| 亚洲国产精品一区二区三区在线| 高清在线国产一区| 激情视频va一区二区三区| 少妇被粗大的猛进出69影院| 国产伦人伦偷精品视频| 亚洲精品国产色婷婷电影| 国产在视频线精品| 亚洲伊人色综图| 国产成人精品久久二区二区91| 亚洲avbb在线观看| 日本a在线网址| 国产单亲对白刺激| 精品福利观看| 天堂√8在线中文| 18禁裸乳无遮挡动漫免费视频| 在线观看66精品国产| 国产一区二区三区综合在线观看| 亚洲免费av在线视频| a级毛片黄视频| 深夜精品福利| 日本vs欧美在线观看视频| 九色亚洲精品在线播放| 天堂动漫精品| 99久久精品国产亚洲精品| 欧美国产精品va在线观看不卡| 一进一出好大好爽视频| 制服诱惑二区| 又黄又粗又硬又大视频| 日韩欧美三级三区| e午夜精品久久久久久久| 一边摸一边做爽爽视频免费| 国产精品.久久久| 精品欧美一区二区三区在线| 黄网站色视频无遮挡免费观看| 90打野战视频偷拍视频| 亚洲第一欧美日韩一区二区三区| 亚洲欧美激情在线| 视频区欧美日本亚洲| 国内久久婷婷六月综合欲色啪| 女人被躁到高潮嗷嗷叫费观| 久久99一区二区三区| 午夜免费鲁丝| 国产成人系列免费观看| 母亲3免费完整高清在线观看| 免费一级毛片在线播放高清视频 | 高清av免费在线| 精品一区二区三区视频在线观看免费 | 久久 成人 亚洲| 欧美在线黄色| 欧美人与性动交α欧美精品济南到| 国产免费男女视频| 亚洲av美国av| 最新在线观看一区二区三区| 在线观看66精品国产| 国产av精品麻豆| 亚洲国产中文字幕在线视频| 黄色 视频免费看| 亚洲av第一区精品v没综合| 淫妇啪啪啪对白视频| 涩涩av久久男人的天堂| 夜夜爽天天搞| 高潮久久久久久久久久久不卡| 日本欧美视频一区| 如日韩欧美国产精品一区二区三区| 高清欧美精品videossex| 亚洲视频免费观看视频| x7x7x7水蜜桃| 变态另类成人亚洲欧美熟女 | 人人妻人人添人人爽欧美一区卜| 日本五十路高清| 亚洲精华国产精华精| 日日爽夜夜爽网站| 亚洲国产毛片av蜜桃av| 久久精品成人免费网站| 亚洲成人手机| 久久精品aⅴ一区二区三区四区| 欧美日韩福利视频一区二区| 精品福利观看| av有码第一页| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品在线观看二区| 高清毛片免费观看视频网站 | 香蕉丝袜av| 天天操日日干夜夜撸| 精品第一国产精品| 欧美精品高潮呻吟av久久| 午夜精品在线福利| 欧美人与性动交α欧美软件| 在线免费观看的www视频| 人妻久久中文字幕网| 久久人人97超碰香蕉20202| 亚洲精品自拍成人| 精品久久久精品久久久| 人妻一区二区av| 99久久综合精品五月天人人| 久久久久久人人人人人| 一级片免费观看大全| 少妇猛男粗大的猛烈进出视频| 国产蜜桃级精品一区二区三区 | 美女视频免费永久观看网站| 男女下面插进去视频免费观看| 亚洲成人手机| 免费在线观看完整版高清| 99久久综合精品五月天人人| 国产深夜福利视频在线观看| 91国产中文字幕| 成年人午夜在线观看视频| 美女高潮到喷水免费观看| 伦理电影免费视频| 国产成人系列免费观看| 国产av又大| 欧美激情久久久久久爽电影 | 一本大道久久a久久精品| www.999成人在线观看| 国产激情久久老熟女| 免费黄频网站在线观看国产| 午夜精品久久久久久毛片777| 国产乱人伦免费视频| 天天添夜夜摸| 91老司机精品| 久久精品aⅴ一区二区三区四区| 十分钟在线观看高清视频www| 免费在线观看影片大全网站| 亚洲精品自拍成人| 美女国产高潮福利片在线看| 亚洲人成77777在线视频| 夫妻午夜视频| 亚洲伊人色综图| tocl精华| 在线观看66精品国产| 午夜亚洲福利在线播放| 国产激情欧美一区二区| 人人妻人人爽人人添夜夜欢视频| 校园春色视频在线观看| av线在线观看网站| 久久精品亚洲av国产电影网| 欧美精品亚洲一区二区| 亚洲欧美一区二区三区黑人| 亚洲美女黄片视频| 国产精华一区二区三区| 老司机在亚洲福利影院| 精品免费久久久久久久清纯 | 91精品三级在线观看| 国产激情久久老熟女| 国产成人系列免费观看| 精品无人区乱码1区二区| 精品国产超薄肉色丝袜足j| av网站免费在线观看视频| 91九色精品人成在线观看| 五月开心婷婷网| 最近最新中文字幕大全电影3 | 亚洲中文字幕日韩| 亚洲av熟女| 757午夜福利合集在线观看| 狠狠婷婷综合久久久久久88av| 精品卡一卡二卡四卡免费| 日韩三级视频一区二区三区| 久久人人爽av亚洲精品天堂| 国产精品国产av在线观看| 成熟少妇高潮喷水视频| av中文乱码字幕在线| 国产精品美女特级片免费视频播放器 | 精品卡一卡二卡四卡免费| 国产一区二区三区综合在线观看| 欧洲精品卡2卡3卡4卡5卡区| 亚洲第一av免费看| 国产成人欧美在线观看 | 别揉我奶头~嗯~啊~动态视频| 妹子高潮喷水视频| 欧美日韩成人在线一区二区| av免费在线观看网站| 男女免费视频国产| 欧美乱码精品一区二区三区| 日本精品一区二区三区蜜桃| 国产亚洲欧美98| 国产精品免费视频内射| 国产午夜精品久久久久久| 国产精品亚洲av一区麻豆| 美国免费a级毛片| 日韩成人在线观看一区二区三区| 一夜夜www| 欧美中文综合在线视频| 国产精品乱码一区二三区的特点 | 国产高清videossex| 午夜福利一区二区在线看| 成人18禁在线播放| 亚洲视频免费观看视频| 制服人妻中文乱码| 热99久久久久精品小说推荐| 久久ye,这里只有精品| 亚洲成人国产一区在线观看| 国产日韩一区二区三区精品不卡| 日韩欧美一区二区三区在线观看 | 亚洲欧美精品综合一区二区三区| 久久久久久久久久久久大奶| 精品福利观看| 丝袜人妻中文字幕| 黑人巨大精品欧美一区二区mp4| 国产成人影院久久av| 亚洲专区国产一区二区| 18在线观看网站| 天堂俺去俺来也www色官网| 黄片大片在线免费观看| 超碰97精品在线观看| 成人黄色视频免费在线看| 日本黄色视频三级网站网址 | 亚洲成人免费电影在线观看| 高清毛片免费观看视频网站 | 久9热在线精品视频| 国产亚洲精品久久久久久毛片 | 男男h啪啪无遮挡| 亚洲中文av在线| 啦啦啦在线免费观看视频4| 成年人免费黄色播放视频| 怎么达到女性高潮| 日韩欧美三级三区| 久久热在线av| 黄片小视频在线播放| 国产精品一区二区在线不卡| 日韩欧美一区二区三区在线观看 | 亚洲在线自拍视频| 精品人妻在线不人妻| 亚洲第一欧美日韩一区二区三区| tube8黄色片| 久久久久久亚洲精品国产蜜桃av| 国产亚洲精品第一综合不卡| 首页视频小说图片口味搜索| 国产日韩欧美亚洲二区| 亚洲人成伊人成综合网2020| 婷婷成人精品国产| 欧美人与性动交α欧美精品济南到| 亚洲精品中文字幕在线视频| 久久精品亚洲av国产电影网| 日本黄色视频三级网站网址 | 一本综合久久免费| 国产av精品麻豆| 欧美乱色亚洲激情| 免费女性裸体啪啪无遮挡网站| 亚洲va日本ⅴa欧美va伊人久久| 免费女性裸体啪啪无遮挡网站| 老司机深夜福利视频在线观看| 免费女性裸体啪啪无遮挡网站| 欧美国产精品va在线观看不卡| 亚洲一码二码三码区别大吗| www日本在线高清视频| 99国产精品一区二区蜜桃av | 捣出白浆h1v1| 亚洲精品av麻豆狂野| 啦啦啦 在线观看视频|