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

    GNSS測(cè)站速度估計(jì)模型重建

    2017-02-09 03:08:14徐克科伍吉倉(cāng)
    測(cè)繪工程 2017年4期
    關(guān)鍵詞:測(cè)站卡爾曼濾波時(shí)序

    徐克科,伍吉倉(cāng)

    (1.河南理工大學(xué) 測(cè)繪與國(guó)土信息工程學(xué)院,河南 焦作 454000;2.同濟(jì)大學(xué) 測(cè)繪與地理信息學(xué)院,上海 200092)

    GNSS測(cè)站速度估計(jì)模型重建

    徐克科1,伍吉倉(cāng)2

    (1.河南理工大學(xué) 測(cè)繪與國(guó)土信息工程學(xué)院,河南 焦作 454000;2.同濟(jì)大學(xué) 測(cè)繪與地理信息學(xué)院,上海 200092)

    獲取準(zhǔn)確的GNSS測(cè)站速度對(duì)于研究全球板塊運(yùn)動(dòng)、地殼形變、地震活動(dòng)及其地球動(dòng)力學(xué)過程至關(guān)重要。為此,以GNSS基線解算后的基線向量作為觀測(cè)值,重建了最小二乘綜合速度解算模型和卡爾曼濾波速度估計(jì)模型。模型中,考慮了測(cè)站坐標(biāo)、運(yùn)動(dòng)速度、年、半年周期項(xiàng)一同作為參數(shù),在基線解網(wǎng)平差的同時(shí),一并求解獲取速度估值。同時(shí),利用坐標(biāo)時(shí)序分析的方法,顧及白噪聲和冪律噪聲的影響,對(duì)單日解坐標(biāo)時(shí)間序列重建了時(shí)序速度擬合估計(jì)模型,以獲取長(zhǎng)期趨勢(shì)項(xiàng)作為速度值。基于3種模型,對(duì)川滇地區(qū)中國(guó)地殼運(yùn)動(dòng)觀測(cè)網(wǎng)絡(luò)2010—2014年21個(gè)GPS基準(zhǔn)站的觀測(cè)數(shù)據(jù)進(jìn)行了速度估計(jì)與對(duì)比分析。結(jié)果表明:3種模型所估計(jì)測(cè)站速度非常接近,差異基本處于0~1 mm/a范圍之內(nèi);速度估值中誤差均在亞毫米水平。由此得出,3種速度估計(jì)模型具有本質(zhì)的一致性,均可正確估計(jì)測(cè)站運(yùn)動(dòng)速度,能夠滿足高精度地殼形變研究的需要。

    GNSS;速度;基線向量;卡爾曼濾波;時(shí)序分析

    要正確理解全球板塊運(yùn)動(dòng)、地殼形變構(gòu)造運(yùn)動(dòng)過程及其動(dòng)力學(xué)機(jī)制,首先必須先獲取高精度、可靠的測(cè)站運(yùn)動(dòng)速度[1-6]。如果只觀測(cè)兩期,通過GNSS高精度處理軟件經(jīng)赫爾模特坐標(biāo)變換可得到這兩期的ITRF框架下的坐標(biāo),兩期的坐標(biāo)相減除以兩期間隔時(shí)間便可得到測(cè)站實(shí)際運(yùn)動(dòng)速度。然而,地殼形變的研究通常需要經(jīng)過長(zhǎng)時(shí)間連續(xù)站連續(xù)觀測(cè)或區(qū)域站多期的持續(xù)觀測(cè)。隨著GNSS觀測(cè)資料的持續(xù)積累,通過GNSS長(zhǎng)期的觀測(cè)數(shù)據(jù)解算的單日解坐標(biāo)時(shí)間序列可以估計(jì)出更高精度的測(cè)站運(yùn)動(dòng)速度。目前,GNSS測(cè)站速度估計(jì)的方法主要可歸結(jié)于兩大類。一是將測(cè)站坐標(biāo)和速度作為參數(shù),建立關(guān)于基線向量的觀測(cè)方程,在對(duì)基線解網(wǎng)平差的同時(shí)一起求解[7-8]。二是對(duì)網(wǎng)平差后單日解測(cè)站坐標(biāo)時(shí)間序列進(jìn)行時(shí)序分析,擬合出線性趨勢(shì)項(xiàng)作為速度[9-12]。有關(guān)研究表明,當(dāng)GNSS坐標(biāo)時(shí)間序列跨度大于4.5 a時(shí),擬合中是否有季節(jié)項(xiàng)的加入對(duì)速度估計(jì)誤差影響可以忽略不計(jì);低于2.5 a,對(duì)速度估計(jì)影響較大,不可接受;當(dāng)時(shí)間跨度是整數(shù)年加上半年時(shí),季節(jié)項(xiàng)對(duì)速度場(chǎng)的影響最小[13-14]。GAMIT軟件的輸出結(jié)果文件H文件提供了單日解基線向量參數(shù)估計(jì)結(jié)果及其協(xié)方差陣[9,15],這為利用合并單天解估計(jì)速度提供了有利條件。對(duì)此,以基線向量作為觀測(cè)值,引入測(cè)站坐標(biāo)、運(yùn)動(dòng)速度、年、半年周期項(xiàng)一同作為參數(shù),重建了基線向量最小二乘綜合解算模型、基線向量卡爾曼濾波解算模型;顧及白噪聲和冪律噪聲的影響構(gòu)建了坐標(biāo)時(shí)序分析速度擬合估計(jì)模型。利用川滇地區(qū)中國(guó)地殼運(yùn)動(dòng)觀測(cè)網(wǎng)絡(luò)2010—2014年GPS基準(zhǔn)站觀測(cè)數(shù)據(jù),分別基于這3種模型進(jìn)行了GNSS測(cè)站運(yùn)動(dòng)速度的估計(jì)和對(duì)比分析。

    1 基線向量最小二乘綜合解算模型

    (1)

    (2)

    (3)

    上述的觀測(cè)方程中,待求參數(shù)為測(cè)站i和j在參考時(shí)刻m的坐標(biāo)、測(cè)站i,j的運(yùn)動(dòng)速度和年、半年周期項(xiàng)系數(shù),觀測(cè)值為基線三維向量,權(quán)陣為GNSS基線解算的協(xié)方差陣的逆陣。這只是基于其中一條基線所列出的誤差方程。同理,可列出所有時(shí)段整個(gè)GNSS網(wǎng)所有基線的觀測(cè)方程并聯(lián)立求解。因基線解算結(jié)果單天解是松弛解,聯(lián)合解的站坐標(biāo)組成的網(wǎng)沒有固定的參考基準(zhǔn),因此需要確定參考框架。采用若干穩(wěn)定的IGS站的站坐標(biāo)和速度進(jìn)行擬穩(wěn)平差或重心基準(zhǔn)平差。

    2 基線向量卡爾曼濾波估計(jì)模型

    利用GAMIT軟件的輸出結(jié)果的H文件提供的單日解基線向量參數(shù)估計(jì)結(jié)果及其協(xié)方差陣,基于動(dòng)態(tài)卡爾曼濾波理論,將測(cè)站坐標(biāo)、測(cè)站速度、年、半年周期項(xiàng)作為具有方差信息的準(zhǔn)觀測(cè)值,建立了動(dòng)態(tài)卡爾曼濾波漂移速度估計(jì)的狀態(tài)模型。對(duì)一系列的單日解進(jìn)行卡爾曼濾波合并。并通過對(duì)全球IGS 核心站的約束來進(jìn)行參考框架的定義和七參數(shù)轉(zhuǎn)換,從而獲取了ITRF參考框架下的測(cè)站的坐標(biāo)和運(yùn)動(dòng)速度[9,15]。因測(cè)站坐標(biāo)參數(shù)是隨機(jī)變化的,采用馬爾科夫過程來描述這些參數(shù)是隨機(jī)漫步過程。將測(cè)站運(yùn)動(dòng)速度、年、半年周期項(xiàng)系數(shù)作為恒定參數(shù),測(cè)站坐標(biāo)作為時(shí)變參數(shù),列出關(guān)于基線向量的卡爾曼濾波模型的觀測(cè)方程和狀態(tài)方程。觀測(cè)方程同第1節(jié)基線向量最小二乘綜合解算模型中的觀測(cè)方程,見式(3)。建立以測(cè)站坐標(biāo)、速度、年、半年周期項(xiàng)系數(shù)作為狀態(tài)向量的狀態(tài)轉(zhuǎn)移方程,見式(4)。

    (4)

    Dt=

    (5)

    從k到k+1時(shí)刻的預(yù)報(bào)公式為

    (6)

    (7)

    從k到k+1時(shí)刻的修正公式為

    (8)

    (9)

    (10)

    式中,K為卡爾曼濾波增益矩陣。

    上述卡爾曼濾波估計(jì)的時(shí)間順序是逐漸增加的,直到最后參考?xì)v元時(shí)刻的參數(shù)最佳估計(jì)值和協(xié)方差矩陣,這是前向卡爾曼濾波。還需要再運(yùn)行后向卡爾曼濾波估計(jì)進(jìn)行回代。后向卡爾曼濾波估計(jì)與前向卡爾曼濾波估計(jì)公式相同,只是時(shí)間順序相反。

    3 坐標(biāo)時(shí)序擬合估計(jì)模型

    對(duì)于長(zhǎng)時(shí)間連續(xù)GNSS坐標(biāo)時(shí)間序列,為獲取GNSS測(cè)站精確的速度場(chǎng)信息,可以對(duì)此采用時(shí)序分析的方法進(jìn)行速度擬合。GNSS單站、單分量坐標(biāo)時(shí)間序列速度擬合函數(shù)模型可表示為

    y(ti)=a+bti+csin(2πti)+dcos(2πti)+

    esin(4πti)+fcos(4πti)+

    (11)

    式中:ti為以年為單位的GPS站點(diǎn)單日歷元,H表示Heaviside階梯函數(shù),參數(shù)a表示站點(diǎn)的起始位置,b表示線速度,c和d組合表示年周期項(xiàng)系數(shù),e和f表示半年周期項(xiàng)系數(shù),第七項(xiàng)表示發(fā)生在歷元Tgj處大小為gj的ng個(gè)偏移常量,即站點(diǎn)由于地震或者更換接收機(jī)天線等原因產(chǎn)生的后續(xù)歷元位置整體偏移。

    在速度估計(jì)時(shí),考慮到有色噪聲對(duì)速度估計(jì)的影響,通過構(gòu)造最大似然函數(shù)估計(jì)譜指數(shù)。對(duì)于冪律譜噪聲,協(xié)因數(shù)陣Jpl可表示成轉(zhuǎn)換矩陣與其轉(zhuǎn)置的乘積[11,16],即

    (12)

    轉(zhuǎn)換矩陣T(a)可表示為

    (13)

    (14)

    式中:σwh,σpl為白噪聲和有色噪聲分量;Rwh,Jpl為白噪聲和有色噪聲協(xié)方差矩陣;Rwh為單位陣。構(gòu)造噪聲分量和譜指數(shù)的極大似然函數(shù)為

    (15)

    式中,n表示時(shí)間序列的總長(zhǎng)度,其對(duì)數(shù)形式為

    (16)

    采用單純形法求解上式中的最大值,得到相應(yīng)的似然值及白噪聲與有色噪聲分量σwh,σpl。研究結(jié)果表明,原始時(shí)間序列的最優(yōu)噪聲模型與PL+VW噪聲模型的極大似然值相差甚微,冪律噪聲的水平可以反映時(shí)間序列的整體噪聲水平[17]。因?yàn)殚W爍噪聲和隨機(jī)漫步噪聲均屬于冪律噪聲類型,在此,采用PL+VW噪聲模型進(jìn)行了速度擬合。

    4 實(shí)例分析

    利用川滇地區(qū)中國(guó)地殼運(yùn)動(dòng)觀測(cè)網(wǎng)絡(luò)21個(gè)基準(zhǔn)站(包括框架站)自2010—2014年的觀測(cè)數(shù)據(jù),經(jīng)GAMIT10.4軟件處理后分別得到了單日基線解和單日坐標(biāo)解,部分站的結(jié)果見圖1和圖2。因不同軟件,不同方法的解算結(jié)果之間的比較,更能客觀地評(píng)價(jià)數(shù)據(jù)處理結(jié)果的外部一致性與可靠性。為此,同時(shí)利用BERNESE5.2軟件也進(jìn)行了解算,作了比較。由圖1基線結(jié)果看,兩種軟件解算的DX、基線長(zhǎng)度DL離散度均在5 mm之內(nèi),DY,DZ離散度在10 mm之內(nèi);兩者基線時(shí)間序列整體吻合度較高,變化趨勢(shì)較一致。由圖2坐標(biāo)解算結(jié)果看,兩種軟件網(wǎng)平差后單日解坐標(biāo)時(shí)間序列離散度接近,基本都在10 mm之內(nèi),中誤差均優(yōu)于3 mm;兩者坐標(biāo)時(shí)序整體趨勢(shì)變化一致??梢?,數(shù)據(jù)處理過程中所解算基線向量和坐標(biāo)結(jié)果可靠,為后續(xù)測(cè)站運(yùn)動(dòng)速度估計(jì)的準(zhǔn)確性提供了保障。

    圖1 基線SCLH-SCXJ、SCDF-SCYX三分量與基線長(zhǎng)度單日解時(shí)間序列

    圖2 測(cè)站SCLH、SCXJ、SCDF和SCYX單日解坐標(biāo)時(shí)間序列及中誤差

    對(duì)于基線解,分別利用第1節(jié)基線向量最小二乘方法和第2節(jié)基線向量卡爾曼濾波方法進(jìn)行了速度估計(jì)。對(duì)于單日解坐標(biāo)時(shí)間序列,利用第3節(jié)時(shí)序分析的方法擬合了測(cè)站運(yùn)動(dòng)速度。3種方法解算的21個(gè)臺(tái)站的速度結(jié)果見表1。由表1所示,3種模型所解算測(cè)站速度值非常接近,之差大都小于1 mm/a;速度估值中誤差均在亞毫米水平。這表明,3種速度估計(jì)模型均可正確解算測(cè)站運(yùn)動(dòng)速度,能夠滿足高精度地殼形變研究的需要。最小二乘和卡爾曼濾波估計(jì)模型結(jié)果十分接近,兩種模型具有本質(zhì)的一致性。時(shí)序擬合模型結(jié)果表明,所有測(cè)站殘差的有色噪聲功率譜指數(shù)基本處于(-1,0)區(qū)間內(nèi),表現(xiàn)為閃爍噪聲。不考慮時(shí)間序列中的有色噪聲會(huì)導(dǎo)致擬合速度精度的嚴(yán)重高估。顧及有色噪聲的影響,GPS測(cè)站擬合速度值變化不大,但是擬合速度中誤差將擴(kuò)大約2~8倍。同時(shí),對(duì)于5 a的GNSS坐標(biāo)時(shí)間序列,是否有季節(jié)項(xiàng)的加入對(duì)速度估計(jì)影響微乎其微,可以忽略。

    表1 3種速度估計(jì)模型的估計(jì)結(jié)果 mm/a

    5 結(jié) 論

    以GAMIT解算的基線向量作為觀測(cè)值,引入測(cè)站坐標(biāo)、運(yùn)動(dòng)速度、年、半年周期項(xiàng)一同作為參數(shù),重建了基線向量最小二乘綜合解算模型、基線向量卡爾曼濾波解算模型。并顧及白噪聲和冪律噪聲的影響構(gòu)建了坐標(biāo)時(shí)序分析速度擬合估計(jì)模型。利用川滇地區(qū)中國(guó)地殼運(yùn)動(dòng)觀測(cè)網(wǎng)絡(luò)2010—2014年GPS基準(zhǔn)站觀測(cè)數(shù)據(jù),分別基于3種模型進(jìn)行了GNSS測(cè)站運(yùn)動(dòng)速度估計(jì)和對(duì)比分析。

    1)基線向量最小二乘綜合解算模型是利用三維自由網(wǎng)解得到的所有基線向量納入統(tǒng)一綜合模型作為觀測(cè)值,基線向量的驗(yàn)后方差-協(xié)方差陣則被用來確定觀測(cè)值得權(quán)陣,將參考時(shí)刻的測(cè)站坐標(biāo)、速度、周年、半周年系數(shù)作為未知參數(shù),把所有時(shí)段所有基線聯(lián)合起來利用最小二乘一同求解。同時(shí)需附加一定的約束條件,如固定其中的某些點(diǎn)的坐標(biāo)速度,或進(jìn)行擬穩(wěn)平差和重心基準(zhǔn)平差等,求出最終解,這通常是法方程迭加過程。隨著觀測(cè)期數(shù)的增加,A陣變得很大,對(duì)計(jì)算機(jī)存儲(chǔ)空間要求較高,數(shù)據(jù)處理變得復(fù)雜。而基線向量卡爾曼濾波估計(jì)勿需保留用過的觀測(cè)值序列,按照一套遞推算法,把參數(shù)估計(jì)和預(yù)報(bào)有機(jī)地結(jié)合起來,計(jì)算速度較快,節(jié)約內(nèi)存。坐標(biāo)時(shí)序擬合模型是對(duì)GNSS單站N,E,U三分量的坐標(biāo)時(shí)間序列進(jìn)行時(shí)序分析,考慮長(zhǎng)期線性運(yùn)動(dòng)+年/半年周期運(yùn)動(dòng)+階躍,基于譜指數(shù)計(jì)算或最小范數(shù)二次無偏估計(jì)(MINQUE)或極大似然估計(jì),可以顧及白噪聲、閃爍噪聲、隨機(jī)游動(dòng)、冪指數(shù)噪聲和帶通濾波噪聲等不同噪聲組合,采用加權(quán)最小二乘法求取測(cè)站長(zhǎng)時(shí)間的運(yùn)動(dòng)速度。同時(shí)還可以對(duì)不同噪聲特性及產(chǎn)生原因進(jìn)行深入分析。

    2)3種模型速度估計(jì)結(jié)果表明,三者速度值十分接近,差異大都位于0~1 mm/a范圍之內(nèi),速度值中誤差均在亞毫米水平。由此驗(yàn)證了3種速度估計(jì)模型的一致性和可靠性,這說明3種模型所估計(jì)的測(cè)站運(yùn)動(dòng)速度均可以滿足高精度地殼形變研究的需要。

    3)時(shí)序分析結(jié)果表明,所有測(cè)站殘差的有色噪聲功率譜指數(shù)基本處于(-1,0)區(qū)間內(nèi)??梢姡? a的數(shù)據(jù)主要表現(xiàn)為閃爍噪聲。不考慮時(shí)間序列中的有色噪聲會(huì)導(dǎo)致擬合速度精度的嚴(yán)重高估。顧及有色噪聲的影響,相比只考慮白噪聲,GPS測(cè)站擬合速度值變化不大,但是擬合速度中誤差將擴(kuò)大約2~8倍。由于隨機(jī)漫步噪聲的長(zhǎng)延時(shí)相關(guān)性,而川滇地區(qū)連續(xù)站坐標(biāo)時(shí)間序列的時(shí)間過短,僅5 a時(shí)間,很可能還不足以準(zhǔn)確量化隨機(jī)漫步噪聲的特征。同時(shí),對(duì)于5 a的GNSS坐標(biāo)時(shí)間序列,是否有季節(jié)項(xiàng)的加入對(duì)速度估計(jì)影響微乎其微,可以忽略不計(jì)。

    [1] 趙國(guó)強(qiáng),蘇小寧.基于GPS獲得的中國(guó)大陸現(xiàn)今地殼運(yùn)動(dòng)速度場(chǎng)[J]. 地震,2014,34(1):97-103.

    [2] 鄒鎮(zhèn)宇,江在森,武艷強(qiáng),等. 基于GPS速度場(chǎng)變化結(jié)果研究汶川地震前后南北地震帶地殼運(yùn)動(dòng)動(dòng)態(tài)特征[J]. 地球物理學(xué)報(bào),2015,58(5):1597-1609.

    [3] 程鵬飛,文漢江,孫羅慶,等. 中國(guó)大陸GPS速度場(chǎng)的球面小波模型及多尺度特征分析[J]. 測(cè)繪學(xué)報(bào),2015,40(10):1063-1070.

    [4] 魏子卿,劉光明,吳富梅. 2000中國(guó)大地坐標(biāo)系:中國(guó)大陸速度場(chǎng)[J]. 測(cè)繪學(xué)報(bào),2011,40(4):403-410.

    [5] 張風(fēng)霜,占偉.利用GNSS連續(xù)觀測(cè)資料獲取高精度動(dòng)態(tài)速度場(chǎng)的研究[J].地震研究,2015,38(1):75-83.

    [6] 王琪,張培震,馬宗晉.中國(guó)大陸現(xiàn)今構(gòu)造變形GPS觀測(cè)數(shù)據(jù)與速度場(chǎng)[J].地學(xué)前緣,2002,9(2):415-429.

    [7] 王解先.由GPS基線向量解算地面形變[J].同濟(jì)大學(xué)學(xué)報(bào)(自然科學(xué)版),2005,33(7):967-970.

    [8] 王解先.GPS精密定軌定位[M].上海:同濟(jì)大學(xué)出版社,1997.

    [9] HERRING T A,KING R W,MCCLUSKY S C. GLOBK Reference Manual:Global Kalman filter VLBI and GPS analysis program[R].Release 10.4.Department of Earth,Atmospheric,and Planetary Sciences,Massachusset Institute of Technology. 2010.

    [10] WILLIAMS S D P.CATS:GPS coordinate time series analysis software[J].GPS solutions. 2008,12(2):147-153.

    [11] LANGBEIN J.Estimating rate uncertainty with maximum likelihood:differences between power-law and flicker-random-walk models[J].Journal of Geodesy,.2012,86(9):775-783.

    [12] WILLIAMS S D P.The effect of coloured noise on the uncertainties of rates estimated from geodetic time series[J].Journal of Geodesy,2003,76(9-10):483-494.

    [13] BLEWITT G,LAVALLéE D.Effect of annual signals on geodetic velocity[J].Journal of Geophysical Research:Solid Earth (1978-2012),2002,107(B7):ETG 9-1-ETG 9-11.

    [14] BOS M S,BASTOS L,FERNANDES R M S.The influence of seasonal signals on the estimation of the tectonic motion in short continuous GPS time-series[J].Journal of Geodynamics, 2010,49(3):205-209.

    [15] HERRING T A,KING R W,MCCLUSKY S C. GAMIT Reference Manual:GPS Analysis at MIT[R].Release 10.4.Department of Earth,Atmospheric,and Planetary Sciences,Massachusset Institute of Technology,2010.

    [16] 姜衛(wèi)平,周曉慧. 澳大利亞GPS坐標(biāo)時(shí)間序列跨度對(duì)噪聲模型建立的影響分析[J]. 中國(guó)科學(xué)(D輯:地球科學(xué)),2014,44(11):2461-2478.

    [17] 吳偉偉. 華北地區(qū)GPS連續(xù)站坐標(biāo)序列特征研究[D].北京:中國(guó)地震局地震預(yù)測(cè)研究所,2014.

    [責(zé)任編輯:劉文霞]

    Velocity estimation model reconstruction for GNSS observation station

    XU Keke1,WU Jicang2

    (1.School of Surveying and Land Information Engineering,Henan Polytechnic University,Jiaozuo 454000,China;2.College of Surveying and Geo-Information,Tongji University, Shanghai 200092,China)

    To acquire reliable velocity of GNSS observation station is very important for the study of global plate movement,crustal deformation,seismic activity and geodynamic process. Therefore, considering the baseline vector after baselines calculation as observation value, introducing the coordinates,speed,annual and semi-annual seasonal items of GNSS stations as parameters together to be estimated, this paper reconstructs the models of least squares estimation and Kalman filter estimation, by which the velocity is determined along with GNSS network adjustment. Meanwhile, considering the effect of the white noise and power law noise, it is same to the velocity fitting estimation model by the coordinate time series using the method of timing series analysis, by which the long-term trend item is estimated and considered as velocity value. Based on three speed estimation models,the GPS data from the Crustal Movement Observation Network of China from 2010 to 2014 in Sichuan-Yunnan area is processed and analyzed.The result shows that the velocity difference between different models is less than 1mm/a and the root mean square error is at a submillimeter level.The three speed estimation models are veritied to be consistent and reliable,and the estimated velocity can satisfy the need the study of high precision crustal deformation..

    GNSS; velocity;baseline vector;Kalman filter;timing series analysis

    引用著錄:徐克科,伍吉倉(cāng).GNSS測(cè)站速度估計(jì)模型重建[J].測(cè)繪工程,2017,26(4):6-11.

    10.19349/j.cnki.issn1006-7949.2017.04.002

    2016-02-17

    國(guó)家重點(diǎn)973基礎(chǔ)研究發(fā)展計(jì)劃資助項(xiàng)目(2013CB733304);國(guó)家自然科學(xué)基金資助項(xiàng)目(41404023)

    徐克科(1979-),男,副教授,博士.

    P228.4

    A

    1006-7949(2017)04-0006-06

    猜你喜歡
    測(cè)站卡爾曼濾波時(shí)序
    時(shí)序坐標(biāo)
    GNSS鐘差估計(jì)中的兩種測(cè)站選取策略分析
    基于Sentinel-2時(shí)序NDVI的麥冬識(shí)別研究
    全球GPS測(cè)站垂向周年變化統(tǒng)計(jì)改正模型的建立
    測(cè)站分布對(duì)GPS解算ERP的影響分析
    基于遞推更新卡爾曼濾波的磁偶極子目標(biāo)跟蹤
    一種毫米波放大器時(shí)序直流電源的設(shè)計(jì)
    電子制作(2016年15期)2017-01-15 13:39:08
    基于模糊卡爾曼濾波算法的動(dòng)力電池SOC估計(jì)
    基于擴(kuò)展卡爾曼濾波的PMSM無位置傳感器控制
    DPBUS時(shí)序及其設(shè)定方法
    河南科技(2014年15期)2014-02-27 14:12:36
    久久午夜福利片| 精品久久国产蜜桃| 亚洲国产欧洲综合997久久,| 国产私拍福利视频在线观看| 男女做爰动态图高潮gif福利片| 久久久久久九九精品二区国产| 欧美不卡视频在线免费观看| h日本视频在线播放| 国产精品人妻久久久久久| 国产在线精品亚洲第一网站| 国产v大片淫在线免费观看| 亚洲不卡免费看| 97超碰精品成人国产| 日本黄色视频三级网站网址| 国产私拍福利视频在线观看| 能在线免费看毛片的网站| 久久精品国产亚洲网站| 久久这里有精品视频免费| 国产老妇女一区| 国产真实伦视频高清在线观看| 亚洲国产欧美在线一区| 亚洲精品影视一区二区三区av| 中文字幕精品亚洲无线码一区| 国产亚洲av片在线观看秒播厂 | 国产片特级美女逼逼视频| 精品久久久久久成人av| 亚洲国产日韩欧美精品在线观看| 男女啪啪激烈高潮av片| 国产精品日韩av在线免费观看| 好男人视频免费观看在线| 女人被狂操c到高潮| 成人欧美大片| 日本黄大片高清| 亚洲美女视频黄频| 亚洲高清免费不卡视频| 亚洲国产精品久久男人天堂| 国产黄片美女视频| 亚洲欧美成人综合另类久久久 | 麻豆av噜噜一区二区三区| 啦啦啦观看免费观看视频高清| 少妇高潮的动态图| 精品国产三级普通话版| 日日干狠狠操夜夜爽| 久久久国产成人免费| 国产在视频线在精品| 日韩,欧美,国产一区二区三区 | 91久久精品国产一区二区三区| 久久九九热精品免费| 亚洲久久久久久中文字幕| 日日摸夜夜添夜夜添av毛片| 日产精品乱码卡一卡2卡三| 99在线人妻在线中文字幕| 亚洲精品久久国产高清桃花| 在线播放国产精品三级| 亚洲经典国产精华液单| 久久99热这里只有精品18| 狠狠狠狠99中文字幕| 日本撒尿小便嘘嘘汇集6| 青春草视频在线免费观看| 欧美日韩综合久久久久久| 欧美日韩精品成人综合77777| 最近视频中文字幕2019在线8| 好男人在线观看高清免费视频| 国产乱人偷精品视频| 国产毛片a区久久久久| 黄片无遮挡物在线观看| 国产高清有码在线观看视频| 校园春色视频在线观看| 国产高清不卡午夜福利| 色吧在线观看| 99久国产av精品国产电影| 级片在线观看| 天堂中文最新版在线下载 | 啦啦啦啦在线视频资源| 久久久欧美国产精品| 国产免费一级a男人的天堂| 美女黄网站色视频| 国内精品宾馆在线| 啦啦啦观看免费观看视频高清| 亚洲国产高清在线一区二区三| 草草在线视频免费看| 两个人视频免费观看高清| 色视频www国产| 小蜜桃在线观看免费完整版高清| 国产高潮美女av| 亚州av有码| 99riav亚洲国产免费| 欧美高清性xxxxhd video| 黄色配什么色好看| 婷婷六月久久综合丁香| 99久久精品国产国产毛片| 色尼玛亚洲综合影院| 岛国毛片在线播放| 美女国产视频在线观看| 久久国内精品自在自线图片| 青春草视频在线免费观看| 精品一区二区三区人妻视频| 丰满乱子伦码专区| 久99久视频精品免费| 午夜精品一区二区三区免费看| 亚洲中文字幕日韩| 长腿黑丝高跟| 精品人妻熟女av久视频| 国产精品不卡视频一区二区| 97人妻精品一区二区三区麻豆| 免费看美女性在线毛片视频| 亚洲av不卡在线观看| 国产成人一区二区在线| avwww免费| 丝袜美腿在线中文| 久久人妻av系列| 亚洲人成网站在线播放欧美日韩| 久久中文看片网| 午夜免费激情av| 国产不卡一卡二| 观看美女的网站| 婷婷色av中文字幕| 亚洲av二区三区四区| 久久99热这里只有精品18| 久久午夜福利片| 国产在视频线在精品| 国产老妇女一区| 国产私拍福利视频在线观看| 亚洲欧美中文字幕日韩二区| 亚洲无线观看免费| 最新中文字幕久久久久| 毛片女人毛片| 午夜免费男女啪啪视频观看| 又粗又爽又猛毛片免费看| 久久久精品94久久精品| 久久精品国产99精品国产亚洲性色| 中文字幕av在线有码专区| 久久久久久九九精品二区国产| 久久久久久久午夜电影| 久99久视频精品免费| 赤兔流量卡办理| 欧美又色又爽又黄视频| 舔av片在线| 好男人视频免费观看在线| 全区人妻精品视频| 国产日韩欧美在线精品| 亚洲欧美成人精品一区二区| 大香蕉久久网| 精品少妇黑人巨大在线播放 | 国产亚洲精品av在线| 久久久久久久久久久免费av| 欧美+亚洲+日韩+国产| 哪个播放器可以免费观看大片| 久久精品国产清高在天天线| 变态另类丝袜制服| or卡值多少钱| 亚洲婷婷狠狠爱综合网| 成人av在线播放网站| 免费看日本二区| 亚洲在线观看片| 国产精品久久久久久久电影| 国产一区二区亚洲精品在线观看| 中出人妻视频一区二区| 国产精品久久久久久久电影| 亚洲av中文av极速乱| 国产伦在线观看视频一区| 精品一区二区三区人妻视频| 国产久久久一区二区三区| 日本一二三区视频观看| 成人特级av手机在线观看| 成人三级黄色视频| 亚洲,欧美,日韩| 色视频www国产| 插阴视频在线观看视频| 欧美日韩乱码在线| 国产一区二区在线av高清观看| 国产精品久久久久久精品电影| 2021天堂中文幕一二区在线观| 午夜福利在线观看免费完整高清在 | 少妇熟女欧美另类| 干丝袜人妻中文字幕| 亚洲人成网站高清观看| 国产精品1区2区在线观看.| 国产精品一及| 亚洲人与动物交配视频| 看免费成人av毛片| 亚洲av成人av| 中出人妻视频一区二区| av在线蜜桃| 激情 狠狠 欧美| 免费看日本二区| av天堂中文字幕网| 国产一区二区三区在线臀色熟女| 国产精品永久免费网站| 日本黄大片高清| 乱人视频在线观看| 精品久久久久久久久久久久久| 色哟哟哟哟哟哟| 噜噜噜噜噜久久久久久91| 晚上一个人看的免费电影| 女同久久另类99精品国产91| 久久久精品欧美日韩精品| 老司机福利观看| 国产老妇女一区| 国产一区二区激情短视频| 亚洲一区高清亚洲精品| av福利片在线观看| 乱人视频在线观看| 亚洲在线自拍视频| 3wmmmm亚洲av在线观看| 午夜精品在线福利| 久久久久国产网址| 有码 亚洲区| 国产精品伦人一区二区| 高清午夜精品一区二区三区 | 国产精品久久电影中文字幕| 天天一区二区日本电影三级| 神马国产精品三级电影在线观看| 国产一区亚洲一区在线观看| 亚洲精品久久国产高清桃花| 在线国产一区二区在线| 五月伊人婷婷丁香| 99国产精品一区二区蜜桃av| 中文字幕精品亚洲无线码一区| 亚洲人成网站高清观看| 一级黄色大片毛片| 久久精品国产自在天天线| 九九热线精品视视频播放| 给我免费播放毛片高清在线观看| 亚洲av不卡在线观看| 日韩成人伦理影院| 又粗又爽又猛毛片免费看| 亚洲av免费在线观看| av视频在线观看入口| 观看美女的网站| 日日摸夜夜添夜夜爱| 夜夜爽天天搞| 99久久精品国产国产毛片| 亚洲av成人av| 天天躁日日操中文字幕| 精品日产1卡2卡| 小蜜桃在线观看免费完整版高清| 久久精品人妻少妇| 丰满乱子伦码专区| 国产精品久久久久久精品电影小说 | 色综合亚洲欧美另类图片| 国产精品女同一区二区软件| 久久久久性生活片| 亚洲欧美成人精品一区二区| 亚洲真实伦在线观看| 国产黄a三级三级三级人| 18禁裸乳无遮挡免费网站照片| 嘟嘟电影网在线观看| 少妇的逼水好多| 亚洲一区二区三区色噜噜| 丰满的人妻完整版| 夜夜爽天天搞| 久久久久久久久久久丰满| 国产免费男女视频| 丰满人妻一区二区三区视频av| 亚洲欧美日韩高清专用| 三级经典国产精品| 深爱激情五月婷婷| 亚洲丝袜综合中文字幕| 精品一区二区免费观看| 国产单亲对白刺激| 日韩av不卡免费在线播放| 久久韩国三级中文字幕| 91精品一卡2卡3卡4卡| 国产精品久久久久久精品电影| 午夜激情福利司机影院| av.在线天堂| 亚洲av一区综合| 亚洲在线自拍视频| 国产伦精品一区二区三区视频9| 在线免费十八禁| 欧美激情久久久久久爽电影| 噜噜噜噜噜久久久久久91| 高清在线视频一区二区三区 | 草草在线视频免费看| 欧美xxxx性猛交bbbb| 老熟妇乱子伦视频在线观看| 美女xxoo啪啪120秒动态图| 十八禁国产超污无遮挡网站| 久久精品国产自在天天线| 成人av在线播放网站| 欧美不卡视频在线免费观看| av在线天堂中文字幕| 久久久a久久爽久久v久久| 国产精品福利在线免费观看| 亚洲国产色片| 高清午夜精品一区二区三区 | 国产熟女欧美一区二区| 嫩草影院新地址| 亚洲欧美日韩高清在线视频| www日本黄色视频网| 能在线免费观看的黄片| 久久久久久久久久久免费av| 美女黄网站色视频| 国内精品宾馆在线| 日本一本二区三区精品| 国产成人精品久久久久久| 久久精品国产自在天天线| 国产精品女同一区二区软件| 熟女人妻精品中文字幕| 成人午夜精彩视频在线观看| 国产一级毛片七仙女欲春2| 99久久中文字幕三级久久日本| 亚洲婷婷狠狠爱综合网| 国产久久久一区二区三区| 亚洲成人av在线免费| 免费不卡的大黄色大毛片视频在线观看 | 国产日韩欧美在线精品| 12—13女人毛片做爰片一| 亚洲欧洲日产国产| 99热全是精品| 波多野结衣巨乳人妻| av天堂在线播放| 91在线精品国自产拍蜜月| 成人二区视频| 日日啪夜夜撸| 亚洲国产色片| 亚洲中文字幕一区二区三区有码在线看| 一区二区三区四区激情视频 | 岛国在线免费视频观看| 中国美白少妇内射xxxbb| 亚洲国产欧美在线一区| 老师上课跳d突然被开到最大视频| 亚洲欧美日韩无卡精品| 国产中年淑女户外野战色| 日本-黄色视频高清免费观看| 观看美女的网站| 国产精品野战在线观看| 国产精品一区二区三区四区久久| 午夜久久久久精精品| av国产免费在线观看| 中文字幕精品亚洲无线码一区| 日本黄色视频三级网站网址| 在线观看午夜福利视频| 少妇高潮的动态图| av在线播放精品| 国产成人午夜福利电影在线观看| 久久精品影院6| 久久精品国产鲁丝片午夜精品| 村上凉子中文字幕在线| 人妻系列 视频| 99久久久亚洲精品蜜臀av| 免费人成视频x8x8入口观看| 特级一级黄色大片| 中文字幕熟女人妻在线| a级一级毛片免费在线观看| 在线免费观看的www视频| 亚洲av不卡在线观看| 久久久午夜欧美精品| 99久久精品国产国产毛片| 国产白丝娇喘喷水9色精品| 中国国产av一级| 啦啦啦韩国在线观看视频| 亚洲av男天堂| 在线播放国产精品三级| 搡女人真爽免费视频火全软件| 日本黄色片子视频| 国产av在哪里看| 国产精品久久久久久精品电影小说 | 亚洲av中文字字幕乱码综合| 九色成人免费人妻av| 男女那种视频在线观看| 在线免费观看的www视频| 直男gayav资源| 亚洲精品日韩av片在线观看| 色视频www国产| 亚洲精品日韩在线中文字幕 | 国产午夜精品久久久久久一区二区三区| 在线免费十八禁| 人体艺术视频欧美日本| 能在线免费看毛片的网站| 久久久久久久久久久免费av| 夜夜爽天天搞| h日本视频在线播放| 岛国在线免费视频观看| 黄色一级大片看看| 亚洲av熟女| 国产免费男女视频| 搡老妇女老女人老熟妇| 少妇熟女aⅴ在线视频| 伊人久久精品亚洲午夜| 国产黄片视频在线免费观看| 少妇熟女aⅴ在线视频| 欧美日韩乱码在线| 伊人久久精品亚洲午夜| 舔av片在线| 久久久久久久久中文| 国产一区二区激情短视频| 免费观看精品视频网站| 又粗又硬又长又爽又黄的视频 | 男人舔女人下体高潮全视频| 精品日产1卡2卡| 亚洲精品日韩av片在线观看| 国产69精品久久久久777片| 此物有八面人人有两片| 国产成人aa在线观看| 九九在线视频观看精品| 免费观看的影片在线观看| 女同久久另类99精品国产91| 日本一二三区视频观看| 国产熟女欧美一区二区| 禁无遮挡网站| 日本黄大片高清| 亚洲av免费高清在线观看| 日本黄色视频三级网站网址| 好男人在线观看高清免费视频| 日本黄大片高清| 久久久色成人| 精品人妻偷拍中文字幕| 51国产日韩欧美| 日韩人妻高清精品专区| 午夜激情欧美在线| 免费看a级黄色片| 国产日本99.免费观看| 国产极品天堂在线| 干丝袜人妻中文字幕| 91aial.com中文字幕在线观看| 成人鲁丝片一二三区免费| 国产真实伦视频高清在线观看| 又爽又黄无遮挡网站| 欧美日本亚洲视频在线播放| 成人二区视频| 一区福利在线观看| 99久久成人亚洲精品观看| 99久久精品热视频| 久久久a久久爽久久v久久| 在线免费观看的www视频| 男女那种视频在线观看| 少妇的逼好多水| 久久这里有精品视频免费| 国产精品日韩av在线免费观看| 免费观看精品视频网站| 亚洲无线在线观看| 国产伦在线观看视频一区| 日韩成人av中文字幕在线观看| 精品少妇黑人巨大在线播放 | 国产精品久久电影中文字幕| 色综合亚洲欧美另类图片| 中文字幕人妻熟人妻熟丝袜美| 丰满人妻一区二区三区视频av| 又黄又爽又刺激的免费视频.| 爱豆传媒免费全集在线观看| 亚洲丝袜综合中文字幕| 最近视频中文字幕2019在线8| 国产伦一二天堂av在线观看| 一区二区三区四区激情视频 | 国产精品一区www在线观看| 大型黄色视频在线免费观看| 精品无人区乱码1区二区| 别揉我奶头 嗯啊视频| 免费观看a级毛片全部| 亚洲av电影不卡..在线观看| 欧美bdsm另类| 日韩欧美国产在线观看| 国产男人的电影天堂91| 少妇的逼水好多| 深爱激情五月婷婷| 色综合色国产| 高清毛片免费观看视频网站| 悠悠久久av| 午夜激情福利司机影院| 亚洲国产精品sss在线观看| 成熟少妇高潮喷水视频| 亚洲av.av天堂| 国产成人a区在线观看| 国产亚洲91精品色在线| 国产精品人妻久久久久久| 一级毛片久久久久久久久女| 伦理电影大哥的女人| 中国美女看黄片| 精品熟女少妇av免费看| 三级毛片av免费| 中文欧美无线码| 一卡2卡三卡四卡精品乱码亚洲| 搡女人真爽免费视频火全软件| 乱码一卡2卡4卡精品| 免费电影在线观看免费观看| 99热这里只有精品一区| 日本与韩国留学比较| 国产精华一区二区三区| 国产白丝娇喘喷水9色精品| 听说在线观看完整版免费高清| 精品久久久久久久久久免费视频| 免费观看a级毛片全部| 少妇熟女欧美另类| 夫妻性生交免费视频一级片| 国产精品日韩av在线免费观看| 亚洲美女搞黄在线观看| 日韩精品青青久久久久久| 能在线免费看毛片的网站| 亚洲最大成人中文| a级毛片a级免费在线| 午夜精品一区二区三区免费看| 久久99热这里只有精品18| 中出人妻视频一区二区| 国产三级中文精品| 一本久久中文字幕| 亚洲经典国产精华液单| 国产高清有码在线观看视频| 亚洲精华国产精华液的使用体验 | 99热这里只有是精品50| 国产成人91sexporn| 人妻久久中文字幕网| 日本与韩国留学比较| 最近的中文字幕免费完整| 国产免费一级a男人的天堂| 最近手机中文字幕大全| 夫妻性生交免费视频一级片| 少妇猛男粗大的猛烈进出视频 | 精品无人区乱码1区二区| 亚洲精品日韩av片在线观看| 中文在线观看免费www的网站| 精品久久久久久久久av| 一个人看视频在线观看www免费| 美女内射精品一级片tv| 成人毛片a级毛片在线播放| 色视频www国产| 精品国产三级普通话版| 在线播放国产精品三级| 婷婷六月久久综合丁香| 成年av动漫网址| av免费观看日本| 99久国产av精品| 我的女老师完整版在线观看| 赤兔流量卡办理| av视频在线观看入口| 欧美日韩一区二区视频在线观看视频在线 | 欧美一级a爱片免费观看看| 久久久精品欧美日韩精品| 欧美不卡视频在线免费观看| 亚洲欧美日韩高清专用| 欧美日韩在线观看h| 又爽又黄a免费视频| 国国产精品蜜臀av免费| 成人高潮视频无遮挡免费网站| 男插女下体视频免费在线播放| 男人舔奶头视频| 亚洲人成网站在线播| 国产高清视频在线观看网站| 亚洲性久久影院| 好男人视频免费观看在线| 欧美一区二区精品小视频在线| 久久精品国产99精品国产亚洲性色| 日日撸夜夜添| 18+在线观看网站| 欧美成人a在线观看| 婷婷精品国产亚洲av| videossex国产| 国产熟女欧美一区二区| 在线观看午夜福利视频| 九草在线视频观看| 一个人观看的视频www高清免费观看| 少妇熟女aⅴ在线视频| 日日干狠狠操夜夜爽| 精品少妇黑人巨大在线播放 | 熟女人妻精品中文字幕| 精品久久久久久久久亚洲| 日本一本二区三区精品| 国产精品女同一区二区软件| 插逼视频在线观看| 麻豆乱淫一区二区| 国产成人91sexporn| 中文字幕av在线有码专区| 免费搜索国产男女视频| 91久久精品电影网| 午夜精品国产一区二区电影 | 男的添女的下面高潮视频| 成人鲁丝片一二三区免费| 久久精品国产清高在天天线| 日韩在线高清观看一区二区三区| 久久午夜福利片| 欧美+亚洲+日韩+国产| 成年av动漫网址| 精品久久国产蜜桃| 白带黄色成豆腐渣| 国产淫片久久久久久久久| 亚洲成人中文字幕在线播放| 亚洲欧美日韩高清在线视频| 亚洲av第一区精品v没综合| 欧美日韩一区二区视频在线观看视频在线 | 欧美高清成人免费视频www| 美女内射精品一级片tv| 亚洲精品乱码久久久久久按摩| 蜜桃久久精品国产亚洲av| 国产一区二区亚洲精品在线观看| 三级经典国产精品| 51国产日韩欧美| 日韩制服骚丝袜av| 黄色欧美视频在线观看| 日日干狠狠操夜夜爽| 我的老师免费观看完整版| 国产精品久久久久久久久免| 天堂√8在线中文| 国产毛片a区久久久久| 国产人妻一区二区三区在| 黄片wwwwww| av在线老鸭窝| 国产一区二区在线av高清观看| 少妇的逼水好多| 在线观看美女被高潮喷水网站| 国产毛片a区久久久久| 97热精品久久久久久| 国产极品精品免费视频能看的| 亚洲成人久久爱视频| 自拍偷自拍亚洲精品老妇| 亚洲av电影不卡..在线观看| 久久久久久久久大av| 日韩欧美国产在线观看| 久久久久久久久久黄片| 国产色婷婷99| 久久久久久大精品| 国产伦在线观看视频一区| av在线播放精品| 在线观看66精品国产| 简卡轻食公司|