陳 嘯，李廣俠，李志強(qiáng)，朱文明（解放軍理工大學(xué)通信工程學(xué)院，江蘇南京210007）

硬件與結(jié)構(gòu)

基于FPGA的低截獲概率衛(wèi)星測控信號捕獲系統(tǒng)*

陳 嘯，李廣俠，李志強(qiáng)，朱文明
（解放軍理工大學(xué)通信工程學(xué)院，江蘇南京210007）

提出了一種基于FPGA的低截獲概率衛(wèi)星測控信號捕獲解決方案。基于AD+FPGA的硬件架構(gòu)完成了帶寬為200 MHz的跳頻和直接序列擴(kuò)頻混合擴(kuò)頻信號的解跳、解擴(kuò)，實現(xiàn)了基于Xi1inx的快速傅里葉變換IP核的擴(kuò)頻碼的快速滑動碼相關(guān)，以及多跳積分結(jié)果的相干累加，并采用相對門限法進(jìn)行了捕獲判決。板上試驗結(jié)果證明了該方案的捕獲性能，最大捕獲時間小于4 s，捕獲概率達(dá)到99%，虛警概率小于0.001。該系統(tǒng)已應(yīng)用于北斗二代導(dǎo)航MEO衛(wèi)星測控鏈路。

FPGA；低截獲概率信號；衛(wèi)星測控；信號捕獲

O 引言

當(dāng)前我國衛(wèi)星測控網(wǎng)主要采用微波直接序列擴(kuò)頻測控體制，具有抗干擾性強(qiáng)、測距精度較高、一站對多星測控等優(yōu)點［1］。一種新的低截獲率混合擴(kuò)頻測控信號綜合了跳頻擴(kuò)頻（FHSS）和直接序列擴(kuò)頻（DSSS）的優(yōu)點，測控信號抗截獲、抗干擾能力得到有效提升，測控精度也隨信號帶寬增加而提高［2］。

低截獲概率信號的合作性接收機(jī)設(shè)計中，信號捕獲設(shè)計最關(guān)鍵，是信號跟蹤的前提，運(yùn)算量通常占同步運(yùn)算的80%以上。其結(jié)合了DSSS和FHSS的低截獲率擴(kuò)頻信號，由于信號頻譜比直擴(kuò)信號進(jìn)一步展寬，因此其自相關(guān)函數(shù)峰更加尖銳，導(dǎo)致捕獲時的時延搜索精度要求更高，另一方面，由于采用跳頻體制，多普勒捕獲問題中又加入了跳頻頻率跳變，需要設(shè)計高效的捕獲算法［3］。

FPGA具有設(shè)計開發(fā)周期短、可在線更新、處理性能優(yōu)異等特點，適用于解決測控通信中的捕獲問題［4-6］。針對低截獲率混合擴(kuò)頻信號的捕獲問題，本文設(shè)計了一種基于FPGA芯片的跳頻直擴(kuò)混合擴(kuò)頻信號捕獲系統(tǒng)，該系統(tǒng)在完成跳頻解跳、直接序列擴(kuò)頻序列剝離的同時，實現(xiàn)快速的跳頻多普勒搜索和時延搜索。對A/D采樣信號進(jìn)行FPGA芯片內(nèi)數(shù)字正交下變頻解跳，使用Xi1inx的IP核進(jìn)行基于FFT-IFFT的并行直擴(kuò)碼解擴(kuò)，用分段積分和相位補(bǔ)償?shù)姆椒▽崿F(xiàn)了大多普勒范圍搜索，實現(xiàn)了較大時間范圍內(nèi)的信號滑動相關(guān)，采用相對門限法進(jìn)行信號捕獲的判決，獲得了優(yōu)異的檢測性能。該捕獲系統(tǒng)芯片資源占用率合理，功耗較低，采用多重高可靠性設(shè)計，適用于衛(wèi)星等空間飛行器的測控系統(tǒng)。

1 捕獲模型

本文研究的低截獲率擴(kuò)頻測控信號是相位連續(xù)的快跳頻加直接序列擴(kuò)頻信號，也稱為DS/FH混合擴(kuò)頻信號，采用二進(jìn)制相位調(diào)制方式（BPSK）。跳頻方式為塊跳頻（b1ock hoPPing），即直擴(kuò)碼相位周期、跳頻相位周期均等于一個調(diào)制符號的時間長度，跳頻圖案和偽碼相位同周期變化，周期均為N跳，信號自相關(guān)峰的周期也為N跳，每個跳頻時段內(nèi)包含一個整周期的直接序列擴(kuò)頻碼。由于衛(wèi)星測控通信屬于高可靠性的低速率通信，一個跳擴(kuò)頻周期內(nèi)的調(diào)制信息是恒定的，那么信號可表示為：

其中，S為發(fā)送信號的平均功率，fi表示跳頻頻點，fRF表示跳頻頻點的射頻基準(zhǔn)頻率，那么實際的射頻跳頻頻率為fi+fRF。Th為一個跳頻頻率的持續(xù)時間，N個Th組成一個跳頻周期。Re［·］為取實部運(yùn)算。sl，i（t）表示第i個跳頻頻率內(nèi)的擴(kuò)頻碼，并且定義：

其中，PN（t）表示偽碼速率為fc的直接序列擴(kuò)頻碼。為了確保各個跳頻時段Th內(nèi)的射頻載波相位起始值相同，跳頻頻率設(shè)計為：

其中，k表示跳頻圖案，fi為1 /Th的k次諧波，射頻跳頻帶寬為BH=K/Th，K為跳頻頻率池的大小，N為所選用的跳頻點的個數(shù)，N?K。擴(kuò)頻碼采用周期長度為M的偽隨機(jī)碼。為了滿足塊跳頻要求，M滿足

TC為一個碼片時間長度。這種擴(kuò)頻信號的處理增益高，能夠?qū)箚我舾蓴_、寬帶干擾等多種樣式的單一和組合干擾。跳頻圖案采用Costas跳頻序列，周期為200跳，PN碼采用周期為1 023的Go1d碼，信號參數(shù)如表1所示。

表1 低截獲概率信號參數(shù)

DS/FH混合信號作為一種擴(kuò)頻信號，信號頻帶寬，且能量均勻分布在整個頻帶內(nèi)，導(dǎo)致信號被噪聲所淹沒，具有很強(qiáng)的抗截獲能力，要實現(xiàn)對此混合信號的捕獲，必須對接收信號進(jìn)行滑動相關(guān)運(yùn)算。如圖1所示，接收機(jī)中頻接收到跳擴(kuò)信號后，經(jīng)過滑動解擴(kuò)和解跳，將信號的跳變頻率和直擴(kuò)偽碼從信號中剝離［8］。當(dāng)滑動相關(guān)時接收信號和本地信號相位重合，即信號恢復(fù)成窄帶信號后，再對信號進(jìn)行累加，濾除帶外噪聲，得到高信噪比的檢測量或檢測信號；最后將其送入檢測器，完成捕獲判決。

圖1 相干累加原理

對于式（1）中的相干跳擴(kuò)頻信號，在背景噪聲為白噪聲、信道中無干擾的條件下，本地跳頻擴(kuò)頻信號以半個碼片（或者更細(xì)?。┑乃阉鏖g隔進(jìn)行滑動，完成載波/偽碼剝離后經(jīng)長度為多跳的相干積分后檢測量的獲得。無信號時檢測量服從自由度為1的中心χ2分布，有信號時檢測量服從自由度為1的非中心χ2分布。

假設(shè)有信號條件下，單跳積分檢測量服從A/σ2=10，σ2=1的高斯分布，并取N=10，那么根據(jù)參考文獻(xiàn)［7］，3種檢測量的包絡(luò)在信號有/無條件下的概率密度函數(shù)如圖2所示。信號捕獲判決準(zhǔn)則采用Neyman-Pearson準(zhǔn)則，即在給定信號虛警概率PF前提下，設(shè)定判斷區(qū)域V，使信號檢測概率PD最大。從圖2中可看出，選多跳相干累加檢測量作為捕獲判決量，在相同的虛警概率條件下能獲得較高的檢測概率。

圖2 檢測量概率分布比較

2 系統(tǒng)設(shè)計方案

2.1 硬件環(huán)境

根據(jù)捕獲運(yùn)算要求，本系統(tǒng)用基于FPGA的硬件設(shè)計方案實現(xiàn)高速并行捕獲邏輯運(yùn)算。FPGA處理的硬件環(huán)境如圖3所示。

圖3 系統(tǒng)硬件環(huán)境

待捕獲的跳擴(kuò)頻信號通過兩路ADC，將模擬信號轉(zhuǎn)換為數(shù)字信號，進(jìn)入FPGA，供FPGA進(jìn)行處理。程序存儲PROM中是本設(shè)計產(chǎn)生的FPGA配置文件，加電后，F(xiàn)PGA讀取配置文件，在時鐘驅(qū)動下，實現(xiàn)開發(fā)人員設(shè)計的邏輯功能。

2.2 分模塊設(shè)計

為實現(xiàn)FPGA內(nèi)部軟件的跳擴(kuò)頻信號捕獲功能，本系統(tǒng)設(shè)計的處理流程如圖4所示。

軟件設(shè)計分為4個功能模塊，其名稱和功能如下：

（1）捕獲控制模塊：AcqCtr1

負(fù)責(zé)捕獲處理流程控制和捕獲結(jié)果的判斷。

（2）正交解跳模塊：AcqFront

在捕獲控制模塊控制下，完成跳擴(kuò)頻信號的解跳。

（3）并行解擴(kuò)模塊：CorrFFT

用FFT-IFFT并行相關(guān)算法完成直接序列擴(kuò)頻信號的解擴(kuò)，負(fù)責(zé)時間長度為1個跳頻時段的信號積分。此模塊負(fù)責(zé)完成分辨率為半個直接序列擴(kuò)頻碼碼片的時延搜索。

（4）多跳累積模塊：PhasCmP

負(fù)責(zé)將N個跳頻時段內(nèi)的積分結(jié)果累加起來，得到半碼片以內(nèi)的時域搜索的捕獲結(jié)果，同時實現(xiàn)（-50 kHz，50 kHz）范圍的多普勒搜索，將累加結(jié)果送給捕獲控制模塊進(jìn)行判斷。綜合以上4個模塊，形成了本軟件的模塊架構(gòu)，如圖5所示。虛線表示控制信號，實線表示信號流。

圖4 軟件捕獲功能流程圖

圖5 FPGA程序結(jié)構(gòu)設(shè)計

2.2.1 正交解跳模塊

為節(jié)省硬件資源，本文算法對應(yīng)的捕獲結(jié)構(gòu)采用了雙下變頻器交替變頻解跳和基于FFT的快速碼相關(guān)解擴(kuò)的方案。圖6表示解跳結(jié)構(gòu)，跳頻帶寬為200 MHz，采樣頻率為280 MHz正交I、Q采樣，用不同圖案表示不同跳頻頻率，τ表示收發(fā)時間差，PN1表示接收第1跳內(nèi)的PN碼，本地下變頻器相位從0開始，到4 kπ結(jié)束，跳頻頻率重合時，解擴(kuò)A、B輸出正確的PN碼。2個下變頻器A、B交替解跳的結(jié)構(gòu)加長了單個跳頻諧波解跳的時長，獲得了2Th時長的解跳結(jié)果，交替解跳結(jié)構(gòu)能夠保留2Th時長的解擴(kuò)輸出，從而保證碼相關(guān)積分時長為Th時，碼滑動相關(guān)的搜索范圍能夠達(dá)到τmax=Th。

圖6 AcqFront模塊工作示意

2.2.2 基于FFT的并行解擴(kuò)模塊

解跳完成后，在芯片內(nèi)對信號進(jìn)行下抽處理，降低后端信號處理速率需求。將280 MHz采樣信號下抽40倍到7 MHz后進(jìn)行處理。圖7表示解擴(kuò)結(jié)構(gòu)，解跳結(jié)果以流水線順序送入FFT核，每一幀點數(shù)為4 096點，同時本地擴(kuò)頻碼也進(jìn)行FFT。經(jīng)處理后，每一跳時間內(nèi)對應(yīng)輸出M= 2 046個解擴(kuò)結(jié)果，存入RAM中。

圖7 基于FFT并行碼相關(guān)結(jié)構(gòu)

2.2.3 多跳累積模塊

將解擴(kuò)模塊輸出的N跳內(nèi)M個碼相位上的共M×N個解擴(kuò)結(jié)果從RAM中讀出之后，按0.1碼片的步進(jìn)進(jìn)行時間搜索和相干累積。圖8所示為相干累積結(jié)構(gòu)。

圖8 N跳相干累積結(jié)構(gòu)

將2NM個數(shù)據(jù)點Pi（m）按跳頻點i分為2M組，逐組對i=1，2，…，N個數(shù)據(jù)進(jìn)行相干累加，相干累加結(jié)果取平方后經(jīng)過一個比較器，保留最大相干累積能量及其對應(yīng)的時間-頻率搜索值。當(dāng)?shù)玫阶畲笙喔衫鄯e能量Pmax及其對應(yīng)的時間-頻率信息后，再取該時間點附近的L個點計算出L個積分結(jié)果，表示碼相位范圍為（-L/2，L/2）個碼片。若Pmax>2·E［P（l）］，則判斷捕獲成功，反之則判定為捕獲失敗。

2.3 實現(xiàn)情況

為實現(xiàn)低截獲概率信號捕獲系統(tǒng)，選用了Xi1inx公司的Kintex7-FFG900芯片，A/D芯片采用2片14位輸入的AD9726芯片。FPGA外部參考時鐘為100 MHz，內(nèi)部工作時鐘為280 MHz。芯片內(nèi)部S1ice使用率為56%，滿足航天設(shè)備芯片降額使用的標(biāo)準(zhǔn)，芯片功耗估計值為4.3 W，包括A/D轉(zhuǎn)換等器件的基帶電路總功耗為7 W。實現(xiàn)了第2節(jié)中低截獲概率跳擴(kuò)頻信號的捕獲，時間捕獲范圍為（-50 ms，50 ms），多普勒頻率捕獲范圍為（-50 kHz，50 kHz），多普勒變化率范圍為（-12 kHz/s，12 kHz/s）。最大捕獲時間不超過4 s，捕獲概率大于99%，虛警概率小于0.001。

3 試驗結(jié)果

對本文的捕獲系統(tǒng)進(jìn)行了基帶電路板仿真驗證，統(tǒng)計了不同的時間搜索精度下，信號源發(fā)射信噪比與檢測概率的關(guān)系如圖9所示。試驗環(huán)境包括安捷倫E4438C信號源2臺、XX型衛(wèi)星應(yīng)答機(jī)射頻通道、基于Xi1inx Kintex-7的基帶板一副，以及電源、連接線若干。信號源負(fù)責(zé)信號上變頻和射頻加噪，基帶板完成發(fā)送和捕獲中頻跳擴(kuò)測控信號的功能。

圖9 相干累積性能比較

電路板驗證時可通過FPGA預(yù)留調(diào)試口調(diào)整時間搜索精度，可選精度有0.01 /0.1 /0.25 /0.5 chiP。由于碼相位搜索間隔以及FFT運(yùn)算的有限字長效應(yīng)［9］，導(dǎo)致板上仿真檢測性能略有下降，距離理論值約有3 dB距離；當(dāng)選取Δτ=0.1 chiP，本文算法在單符號信噪比為17 dB時可達(dá)到99%的捕獲概率；當(dāng)碼相位搜索間隔從0.25 chiP擴(kuò)大到0.5 chiP時，檢測性能急劇下降。

4 結(jié)論

為解決低截獲概率混合擴(kuò)頻的快速捕獲問題，選擇了以FPGA作為處理核心的基帶解決方案［10］，設(shè)計實現(xiàn)了基于FFT分段碼相關(guān)和多跳相干累積的捕獲系統(tǒng)并進(jìn)行了板上試驗。本算法捕獲精度達(dá)到0.01 chiP，捕獲到的DS/FH信號可直接送入載波環(huán)路進(jìn)行載波跟蹤，能減輕跟蹤環(huán)路負(fù)擔(dān)。實際應(yīng)用中，可根據(jù)捕獲精度需求選擇時間捕獲分辨率，大小從0.01 chiP到0.5 chiP。經(jīng)在軌測試后，本系統(tǒng)已成功運(yùn)用于北斗二代導(dǎo)航衛(wèi)星的某新型衛(wèi)星測控系統(tǒng)中。

［1］沈榮駿，趙軍.我國航天測控技術(shù)的發(fā)展趨勢與策略［J］.宇航學(xué)報，2001，22（3）：1-5

［2］楊文革，王金寶，孟生云，等.DS/FH混合擴(kuò)頻信號捕獲方法綜述［J］.測控技術(shù)，2009，28（z1）：16-20.

［3］KAPLAN E D，HEGARTY C J，Understanding GPS，Princi-P1es and aPP1ications（Second Edition）［M］.Norwood，MA：Artech House，Inc.ch.5.4，173-200.

［4］王水魚，馮曉靖.一種新型de1ta-sigma小數(shù)分頻器的FPGA實現(xiàn)［J］.微型機(jī)與應(yīng)用，2013，32（5）：25-27，31.

［5］彭奇，習(xí)友寶.基于FPGA和改進(jìn)CORDIC算法的NCO設(shè)計與實現(xiàn)［J］.微型機(jī)與應(yīng)用，2013，32（5）：60-62.

［6］宋萬均，張安堂.雙基地雷達(dá)目標(biāo)速度計算的FPGA實現(xiàn)［J］.電子技術(shù)應(yīng)用，2014，40（1）：47-49，52.

［7］MARCUM F I.A statistica1 theory of target detection by Pu1sed radar，rand research memo［R］.RM-754，December 1947，with APPendix，RM-753，Ju1y 1948.

［8］潘點飛，程乃平，郝建華，等.DS/FH測控系統(tǒng)抗干擾性能分析與測試［J］.空間科學(xué)學(xué)報，2013，33（5）：540-547.

［9］黃磊，張其善，寇艷紅，等.GPS信號FFT捕獲算法的有限字長效應(yīng)分析［J］.宇航學(xué)報，2006，27（5）：1029-1033.

［10］陳巍，王國富，張法全，等.基于FPGA同步時鐘測量系統(tǒng)的研究及實現(xiàn)［J］.電子技術(shù)應(yīng)用，2014，40（12）：86-88，91.

陳嘯（1986 -），通信作者，男，博士生，主要研究方向：衛(wèi)星導(dǎo)航，衛(wèi)星測控。E-mai1：ty_chenxiao@163.com。

李廣俠（1964 -），男，碩士，教授，博士生導(dǎo)師，主要研究方向：衛(wèi)星通信，衛(wèi)星導(dǎo)航，衛(wèi)星測控。

李志強(qiáng)（1974 -），男，博士，教授，主要研究方向：衛(wèi)星導(dǎo)航，衛(wèi)星測控。

FPGA-based acquisition system of a sate11ite LPITT&C signa1

Chen Xiao，Li Guangxia，Li Zhiqiang，Zhu Wenming
（Co11ege of Communications Engineering，PLA University of Science and Techno1ogy，Nanjing 210007，China）

An FPGA based acquisition so1ution for 1ow-Probabi1ity-of-intercePtion signa1 is ProPosed.W ith AD+ FPGA architecture，the demodu1ation of a 200 MHz-wide frequency hoPPing and direct sequence（FHDS）sPread sPectrum signa1 is achieved，the fast corre1ation of PN code based on the FFT IP core of Xi1inx CorP.is achieved and the accumu1ation of the integra1 in mu1tiP1e hoPs is imP1emented.A re1ative thresho1d method is aPP1ied in signa1acquisition.The test-board exPeriment resu1ts va1idate the acquisition Performance of this so1ution.The maximum acquisition time is 1ess than 4 s，the Probabi1ity of detection is over 99%and the fa1se a1arm Probabi1ity is 1ess than 0.001.This design has been aPP1ied to the TT&C sub-system of MEO sate11ites in Beidou navigation conste11ation.

FPGA；1ow-intercePtion-Probabi1ity signa1；sate11ite tracking；signa1 acquisition

TN927

A

10.19358 /j.issn.1674-7720.2016.09.010

陳嘯，李廣俠，李志強(qiáng)，等.基于FPGA的低截獲概率衛(wèi)星測控信號捕獲系統(tǒng)［J］.微型機(jī)與應(yīng)用，2016，35（9）：30-33.

國家自然科學(xué)基金（61571464）

2016-01-05）