導(dǎo)航衛(wèi)星發(fā)射端多徑的高精度測量技術(shù)*

鄭長剛，黃智剛，康成斌(.北京航空航天大學(xué)電子信息工程學(xué)院,北京000;2.北京市第十三中學(xué),北京00009; .中國空間技術(shù)研究院,北京00094)

?

導(dǎo)航衛(wèi)星發(fā)射端多徑的高精度測量技術(shù)*

鄭長剛**1，2，黃智剛1，康成斌3
(1.北京航空航天大學(xué)電子信息工程學(xué)院,北京100011;2.北京市第十三中學(xué),北京100009; 3.中國空間技術(shù)研究院,北京100094)

**通信作者:zhjinjun@ vip. sina. com Corresponding author:zhjinjun@ vip. sina. com

摘 要:針對國內(nèi)外學(xué)者關(guān)于星體是導(dǎo)致“北斗”衛(wèi)星碼載波偏差隨觀測仰角存在相關(guān)性的主要原因的懷疑,提出了一種導(dǎo)航衛(wèi)星發(fā)射端星體多徑的高精度測試方法。利用連續(xù)頻率調(diào)制信號的時域合成,獲得窄時域脈沖,通過時域窗函數(shù)截取多徑信號,從而獲得多徑的高精度延時和幅度。理論分析和試驗(yàn)驗(yàn)證表明:所提方法簡單有效,獲得了較高的時域分辨性能和測量精度,由此得出了導(dǎo)航衛(wèi)星星體多徑不是產(chǎn)生碼載波偏差主要原因的結(jié)論。

關(guān)鍵詞:導(dǎo)航衛(wèi)星；發(fā)射端；多徑測量；碼載波偏差

1 引 言

“斗衛(wèi)”導(dǎo)航系統(tǒng)(BD Satellite System,BDS)包括14顆中高軌道衛(wèi)星(Medium Earth Orbit,MEO)、5顆傾斜軌道同步衛(wèi)星(Inclined Geosynchronous Satellite Orbit,IGSO)和5顆靜止軌道衛(wèi)星(Geosynchronous Earth Orbit,GEO),播發(fā)的B1、B2和B3 3路導(dǎo)航信號被全球各地多模全球?qū)Ш叫l(wèi)星系統(tǒng)(Global Navigation Satellite System,GNSS)接收機(jī)接收和分析研究。Hauschild、Perello Gisbert和Montenbruk等人發(fā)現(xiàn)BDS存在碼-載波偏差(同一顆衛(wèi)星播發(fā)的偽碼相位與載波相位之差)變化,偏差值與觀測仰角存在相關(guān)性,變化范圍約1 m。由于該偏差值基于遍布全球的不同類型接收機(jī)觀測數(shù)據(jù)而獲得的[1-3],因此被認(rèn)為是與衛(wèi)星密切相關(guān)。德國人Wanninger[4]分別給出了BDS中MEO和IGSO衛(wèi)星的碼-載波偏差修正模型,該模型分別適用于所有MEO衛(wèi)星和IGSO衛(wèi)星。類似的碼-載波偏差曾經(jīng)在GPS SVN49衛(wèi)星上出現(xiàn)過,GPS確認(rèn)是由于衛(wèi)星內(nèi)部信號反射引起的碼多徑導(dǎo)致[5]。

本文研究衛(wèi)星星本體多徑的高精度測試方法和技術(shù),利用BDS衛(wèi)星模型進(jìn)行多徑模擬測試試驗(yàn),掌握BDS衛(wèi)星星體多徑數(shù)據(jù),進(jìn)而給出碼-載波偏差與星體多徑的關(guān)聯(lián)性?；舅悸肥窍仁崂矶鄰綔y量方法,識別出適用于衛(wèi)星發(fā)射端的,并針對衛(wèi)星測量條件進(jìn)行改進(jìn),提高測量精度;然后,根據(jù)BDS系統(tǒng)MEO衛(wèi)星星表技術(shù)狀態(tài),建立物理模型和試驗(yàn)環(huán)境,驗(yàn)證衛(wèi)星多徑影響,導(dǎo)出與碼-載波偏差的關(guān)聯(lián)性;最后,針對衛(wèi)星提出消除自身多徑的建議。

2 多徑測量方法

多徑測量一般采用時域脈沖法[6]、多徑估計法[7]和零值法。時域脈沖法是在地面對星表多徑測量較為簡便的方法,但由于超窄帶的時域脈沖在頻域上將呈現(xiàn)超寬帶特性,受信道特性影響較大,在高精度測量中的實(shí)用性受到限制。多徑估計是利用衛(wèi)星發(fā)播的導(dǎo)航信號,地面接收機(jī)接收后通過各種估計算法獲得多徑的估值,一方面估值精度與信號帶寬、估值算法密切相關(guān),另一方面無法排除環(huán)境和地面設(shè)備硬件的影響,不宜采用。零值法是通過檢測接收機(jī)偽距測量扣除真實(shí)距離后的零值和的變化來判斷多徑效應(yīng)的變化,零值的穩(wěn)定性決定了多徑效應(yīng)變化的分辨率,且受時鐘影響較大,也不宜采用。

本文提出一種改進(jìn)的測量方法,采用連續(xù)頻率調(diào)制信號進(jìn)行時域合成,規(guī)避信道特性不理想的影響,從而獲得超窄帶時域脈沖,獲得高精度的多徑測量數(shù)據(jù)。

2. 1 信號模型與時頻處理

假設(shè)衛(wèi)星發(fā)播信號為( ) s t ,那么接收到的多徑信號可表示為

僅考慮星表多徑影響,可忽略衛(wèi)星導(dǎo)航信號星上內(nèi)部生成過程中反射多徑。采用連續(xù)頻率調(diào)制信號作為發(fā)播信號,即

式中:α(t)為瞬時幅度;φ(t)為瞬時相位,其瞬時頻率ω(t)= dφ(t) / dt,掃頻速率v=dω/ dt,起止頻率分別為ωL和ωH。

假定掃頻過程為平穩(wěn)隨機(jī)狀態(tài),將接收到的信號在時間尺度上疊加,則可認(rèn)為發(fā)播信號在頻域上表現(xiàn)為寬帶矩形脈沖信號,疊加的時域信號表現(xiàn)為sinc函數(shù),如圖1所示。

圖1 時間尺度疊加后信號的時頻特性Fig. 1 Time-frequency characteristic of signal overlaying time information

如果在時間尺度上覆蓋了起止頻率,則合成信號的帶寬最寬,此時,獲得的時域脈沖寬度最窄。當(dāng)脈沖寬度足夠小,則可以分辨出多徑信號中的所有時延因子,實(shí)現(xiàn)高精度的多徑測量。

該方法要求掃頻過程平穩(wěn)、信道為寬帶,但由于實(shí)際信道特性一般是帶限且有帶內(nèi)非平坦,導(dǎo)致實(shí)際收到信號是理想信號通過了一個有色信道。設(shè)有色信道傳遞函數(shù)為( ) H t ,則有

此時,需要對信道進(jìn)行白化,發(fā)播信號通過一個預(yù)失真濾波器,預(yù)失真濾波器的傳遞函數(shù)為H-1(t )。此時,可以獲得預(yù)期的接收信號:

多徑信號從主信號中分離出來后,通過時間域門函數(shù)可將多徑信號提取。對多徑信號作傅里葉變換,得到多徑信號的頻域函數(shù)F(ω),即可獲得關(guān)心頻率對應(yīng)的多徑信號強(qiáng)度。

2. 2 測試方法及驗(yàn)證

利用信號源發(fā)射線性調(diào)頻信號為多徑測試信號,經(jīng)發(fā)射天線發(fā)射后,主路徑信號和多徑信號被接收天線接收,利用幅相檢測接收機(jī)檢測寬帶掃頻信號在不同頻點(diǎn)的幅度和相位,對數(shù)據(jù)進(jìn)行反傅里葉變換時域合成,從而分辨出多徑信號時延和相對主信號強(qiáng)度。測試系統(tǒng)示意如圖2所示。

圖2 多徑信號測試系統(tǒng)Fig. 2 Measurement system for multi-path signal

對測試方法和所采用技術(shù)進(jìn)行桌面有線試驗(yàn)驗(yàn)證,構(gòu)建試驗(yàn)框圖如圖3所示。

圖3 桌面有線試驗(yàn)Fig. 3 Experiment of desktop test

試驗(yàn)前采用矢量網(wǎng)絡(luò)分析儀對路徑1與路徑2時延差和損耗差進(jìn)行標(biāo)定,路徑1與路徑2時延差約2.35 ns,路徑1與路徑2損耗差約19.4 dB。分別采用15～20 GHz和8～20 GHz兩種掃頻信號帶寬進(jìn)行試驗(yàn),得到的結(jié)果如圖4所示,具體數(shù)據(jù)見表1。

表1 接收不同掃頻信號的多徑估計值Tab. 1 Value of assessment of multipath to different scanning signals

表1數(shù)據(jù)顯示利用掃頻法進(jìn)行測試的結(jié)果與用矢網(wǎng)標(biāo)定的結(jié)果一致,時延測試誤差在0. 15 ns以內(nèi),強(qiáng)度誤差在1. 5 dB以內(nèi)。

3 測試試驗(yàn)

星體模擬件放置在微波暗室的轉(zhuǎn)臺上,利用平面近場構(gòu)建遠(yuǎn)場條件[8],根據(jù)天線的收發(fā)互易特性,星上導(dǎo)航信號發(fā)射天線接收測試信號,形成多徑測試系統(tǒng),如圖5所示,通過轉(zhuǎn)動轉(zhuǎn)臺模擬不同用戶的觀測方位,獲取各個方向的測量數(shù)據(jù)。

圖4 不同掃描帶寬的試驗(yàn)結(jié)果圖Fig 4 Result of different bandwidth of scanning

由圖4可見,當(dāng)掃頻帶寬足夠?qū)挄r,可在時域上將不同時延的多徑信號分離,帶寬越寬,分離精度越高。

圖5 暗室掃頻測試系統(tǒng)示意圖Fig. 5 Diagram of test system for frequency scanning in unreflected chamber

試驗(yàn)前,在微波暗室利用標(biāo)準(zhǔn)喇叭天線替代模擬件進(jìn)行暗室環(huán)境多徑測試,用于消除環(huán)境對測試結(jié)果的影響。

為測試星體安裝面及安裝其上的其他設(shè)備對B1/ B2/ B3天線的影響,最好是先能夠獨(dú)立對導(dǎo)航信號天線進(jìn)行多徑測試,獲取基礎(chǔ)數(shù)據(jù)。但由于無法保證天線安裝到星體模擬件后與原狀態(tài)位置完全一致,需要采取其他不影響狀態(tài)的辦法達(dá)到相應(yīng)的效果。

采用星體突出物包敷吸波材料方法,獲取削弱星體安裝面影響的多徑測試數(shù)據(jù),將該數(shù)據(jù)作為基礎(chǔ)數(shù)據(jù)與不包敷時測試結(jié)果對消,即可得到星體對B1/ B2/ B3天線多徑環(huán)境的影響。

按照上述步驟,分別測得暗室環(huán)境、包敷吸波材料和未包敷吸波材料3組數(shù)據(jù)。利用第一組數(shù)據(jù)檢查環(huán)境多徑大小,相對強(qiáng)度均在-60 dB以下,可忽略。通過第三組數(shù)據(jù)與第二組數(shù)據(jù)做差,獲得星體引入的多徑數(shù)據(jù),如圖6所示。

圖6 做差處理后的多徑信號與方位角關(guān)系Fig. 6 Correlation between multipath signals after subtraction and azimuth angle

數(shù)據(jù)處理結(jié)果表明,最強(qiáng)的多徑信號出現(xiàn)在延遲3～5 ns的位置,其強(qiáng)度遠(yuǎn)遠(yuǎn)小于主信號,相對強(qiáng)度均在-34 dB以下。多徑與B1/ B2/ B3天線的方位角和離軸角沒有明顯的關(guān)聯(lián)性,僅在高離軸角時存在時延的波動。

在多徑信號強(qiáng)度較弱的時候,由于兩次測量存在各自獨(dú)立的測量誤差,做差處理會引入隨機(jī)誤差的疊加,表現(xiàn)為系統(tǒng)性的多徑誤差淹沒在噪聲中,導(dǎo)致處理后數(shù)據(jù)的不可信。因此,仍需要針對測量原始數(shù)據(jù)進(jìn)行獨(dú)立分析。

同時,采用吸波材料包敷并非理想,材料對該頻段信號的衰減不大會導(dǎo)致系統(tǒng)性的多徑數(shù)據(jù)被抵消,導(dǎo)致做差處理后數(shù)據(jù)結(jié)果偏離真實(shí)情況。因此,需要對兩組原始數(shù)據(jù)進(jìn)行獨(dú)立分析。

提取方位角為0°、5°和離軸角為45°、5°兩種狀態(tài)的原始測量數(shù)據(jù),如圖7所示,均顯示在5 ns左右有較強(qiáng)多徑信號,相對強(qiáng)度最大為-28 dB,與做差處理后數(shù)據(jù)相比更加可靠。

圖7 星體多徑測試結(jié)果時域圖Fig. 7 Time domain result of multipath test for satellite

重新對未包敷吸波材料的數(shù)據(jù)進(jìn)行分析,針對所有方向的原始數(shù)據(jù)進(jìn)行處理,結(jié)果如圖8所示。

圖8 第三組數(shù)據(jù)的多徑信號與方位角關(guān)系Fig. 8 Correlation between multipath signals and azimuth angle of the third set of data

第三組數(shù)據(jù)的處理結(jié)果與做差結(jié)果存在一定差異,如表2所示。由表2可知各方向最強(qiáng)多徑信號的相對強(qiáng)度均值增加了12 dB。

表2 不同測試數(shù)據(jù)處理方法結(jié)果Tab. 2 Result of different data processing methods

盡管星體多徑強(qiáng)度相對主徑的絕對數(shù)值較小,但可能會對用戶產(chǎn)生一定影響,需要進(jìn)行評估。人為產(chǎn)生不同強(qiáng)度的多徑信號,利用基于標(biāo)準(zhǔn)處理算法的地面終端進(jìn)行接收,不做多徑抑制處理,結(jié)果顯示:當(dāng)多徑信號相對強(qiáng)度達(dá)到-20 dB時,僅會引起偽距測量數(shù)據(jù)發(fā)生6 cm的變化,不會導(dǎo)致接近1 m的碼載波偏差變化?？梢?星體多徑不是產(chǎn)生碼載波偏差隨觀測仰角變化的主要原因。

4 結(jié)束語

本文利用連續(xù)頻率調(diào)制信號時域合成所獲得的窄時域脈沖實(shí)現(xiàn)了導(dǎo)航衛(wèi)星發(fā)射端多徑的高精度測量,可以達(dá)到亞納秒級時域分辨性能和2 dB以內(nèi)的測量精度;通過對模擬星體的高精度多徑測量,發(fā)現(xiàn)盡管衛(wèi)星可能是導(dǎo)致碼載波偏差隨觀測仰角存在相關(guān)性的原因,但星體反射多徑不是其主要成因;同時,也注意到星體多徑是存在的,會引起高精度用戶偽距測量精度產(chǎn)生厘米級的系統(tǒng)偏差,該誤差呈現(xiàn)系統(tǒng)特性,是可以消除的。由于試驗(yàn)采用的是模擬件,無法直接對碼載波偏差進(jìn)行測量,如能夠?qū)φ鎸?shí)衛(wèi)星產(chǎn)品在地面進(jìn)行測量,將能夠給出碼載波偏差是否與衛(wèi)星直接關(guān)聯(lián)的結(jié)論,并對其他可能產(chǎn)生碼載波偏差的可能,如星上天線輻射特性隨觀測仰角變化情況等,進(jìn)行測量。建議衛(wèi)星研制方下一步深入開展相關(guān)工作。

參考文獻(xiàn):

[1] HAUSCHILD A,MONTENBRUCK O,SLEEWAEGEN J M,et al. Characterization of compass M-1 signals[J]. GPS Solutions,2012,16(1):117-126.

[2] HAUSCHILD A,MONTENBRUCK O,THOELERT S,et al. A multi - technique approach for characterizing the SVN49 signal anomaly,part I:receiver tracking and IQ constellation[J]. GPS Solutions,2012,16(1):19-28.

[3] GISBERT J V P,BATZILIS N,RISUENO G L,et al. GNSS payload and signal characterization using a 3mdish antenna[C] / / Proceedings of the 25th International Technical Meeting of The Satellite Division of the Institute of Navigation(ION GNSS 2012). Nashville,TN:ION, 2012:347-356.

[4] WANNINGER L, BEER S. BeiDou satellite - induced code pseudo-range variations:diagnosis and therapy[J]. GPS Solutions,2015,19(4):1-10.

[5] SPINGER T,DILSSNER F. SVN49 and other GPS anomalies[J]. Inside GNSS,2009(1):32-36

[6] 陳祝明,丁義元,向敬成.采用Chirp-Z變換提高LFMCW雷達(dá)的測距離精度[J].信號處理,2002(2): 110-112.

CHEN Zhuming,DING Yiyuan,XIANG Jingcheng. Improving range precision of LFMCW radar by Chirp - Z transform[J]. Signal Processing,2002,18(2):110-112. (in Chinese)

[7] 李林,姬紅兵.基于L-Wigner分布的多徑ESM信號估計[J].系統(tǒng)工程與電子技術(shù),2009,31(11):2618 -2621.

LI Lin,JI Hongbing. Estimation of multipath ESM signals based on L-Wigner distribution[J]. Systems Engineering and Electronics,2009,31(11):2618-2621. (in Chinese)

[8] 劉慶云,陸飛飛,朱偉強(qiáng),等.輻射源細(xì)微特征用于個體識別的可行性分析[J].航天電子對抗,2008,24 (2):40-42.

LIU Qingyun,LU Feifei,ZHU Weiqiang,et al. Feasibility analysis of emitter identification based on emitter fine feature[J]. Aerospace Electronic Warfare,2008,24(2):40-42. (in Chinese)

鄭長剛(1998—),男,江蘇人,主要研究方向?yàn)橥ㄐ排c電子技術(shù)和測量技術(shù);

ZHENG Changgang was born in Jiangsu Province,in 1998. His research concerns communication and electronic technology,and measurement technology.

Email:zhjinjun@ vip. sina. com

黃智剛(1962—),男,博士,教授,主要研究方向?yàn)橥ㄐ排c信息系統(tǒng)技術(shù);

HUANG Zhigang was born in 1962. He is now a professor with the Ph. D. degree. His research concerns communication and information system technology.

康成斌(1981—),男,博士,高級工程師,主要研究方向?yàn)橥ㄐ排c電子技術(shù)。

KANG Chengbin was born in 1981. He is now a senior engineer with the Ph. D. degree. His research concerns communication and electronic technology.

簡訊

關(guān)于《電訊技術(shù)》投稿的重要說明

為了規(guī)范投稿,保證“雙向匿名審稿”的順利實(shí)施,特作如下說明:

(1)稿件需自行脫密處理,在初審階段不需要向編輯部提供相關(guān)材料;

(2)投稿流程中必填項(xiàng)必須準(zhǔn)確、完整填寫,“備注”欄可填寫與所投稿件或作者相關(guān)的說明;

(3)作者提供的信息必須真實(shí),不得存在學(xué)術(shù)不端行為,歡迎廣大作者和讀者監(jiān)督,一經(jīng)查實(shí),將嚴(yán)肅處理;

(4)本刊對稿件排版格式?jīng)]有特別要求,但可參照本刊已發(fā)表論文格式或數(shù)字優(yōu)先出版模版(可在本刊網(wǎng)站“相關(guān)下載”欄目下載)編寫;

(5)來稿送審時,編輯部將刪除稿件中所有作者及其單位或機(jī)構(gòu)信息,但投稿時仍需提供完整信息;

(6)如推薦審稿人,請務(wù)必將其信息準(zhǔn)確、完整填寫,無email的審稿人信息無效。

本刊編輯部從未委托任何機(jī)構(gòu)或個人代理征稿事宜,請勿上當(dāng)!

非本刊網(wǎng)站www. teleonline. cn投稿一切與本刊無關(guān)!

本刊編輯部

High Precision Measurement of Multipath for Emission Platform of Navigation Satellite

ZHENG Changgang1,2,HUANG Zhigang1,KANG Chengbin3
(1. School of Electronic and Information Engineering,Beihang University,Beijing 100191,China; 2. Beijing No. 13 High School,Beijing 100009,China; 3. China Academy of Space Technology,Beijing 100094,China)

Abstract:In view of the scholars' review that satellite is the main cause of correlation between divergences of code-carrier of BD satellites(BDS) and observation angle,a high precision method of measuring multipath caused by satellite is presented. This method gets narrow time pulse by using continous frequency modulated signal,and acquires interested multipath signal by time domain window function,in order to obtain the high precision time delay and intensity of multipath signals. By theoretical analysis and test verification,the method is proved to be simple and effective,and high performances of time domain distinguishment and measurement precision are achieved. Based on the method,the conclusion is drawn that multipath of satellite is not the cause of the divergences of code-carrier.

Key words:navigation satellite;emission platform;measurement of multipath;code-carrier divergence

doi:10. 3969/ j. issn. 1001-893x. 2016. 02. 007引用格式:張劍,周興建,盧建川.基于Takagi-Sugeno-Kang模糊集合的噪聲干擾檢測方法[J].電訊技術(shù),2016,56(2):151-155. [ZHANG Jian, ZHOU Xingjian,LU Jianchuan. A Takagi-Sugeno-Kang fuzzy zpproach to noise jamming detection[J]. Telecommunication Engineering,2016, 56(2):151-155. ]

作者簡介:

中圖分類號:TN966

文獻(xiàn)標(biāo)志碼:A

文章編號:1001-893X(2016)02-0145-06

*收稿日期:2015-11-18;修回日期:2016-01-14 Received date:2015-11-18;Revised date:2016-01-14