非相干累加中檢波方式的性能分析與參數選取?

羅 奎，程 劍，李仰志，賴廣鋒

（解放軍理工大學通信工程學院，南京210007）

非相干累加中檢波方式的性能分析與參數選取?

羅 奎，程 劍，李仰志，賴廣鋒

（解放軍理工大學通信工程學院，南京210007）

闡述了GNSS信號捕獲過程中非相干能量累加的兩種檢波方式的實現原理，利用雷達檢測理論對兩者的性能進行了分析對比與數值計算，給出了兩種檢波方式的性能曲線。提出了兩種相干累加時間與累加次數的選取方法并進行了仿真，該方法和仿真結果可作為GNSS系統信號捕獲與跟蹤過程方案設計的參考依據。仿真結果還表明，平方律檢波方式相比于包絡檢波方式，在能量累加性能上要略差，差值在0.5 dB以內。

GNSS；信號捕獲；擴頻；非相干累加；包絡檢波；平方律檢波

1 引言

在GNSS（全球導航衛星系統）及其它擴頻系統的信號捕獲過程中，信號能量的累加是所有捕獲算法的基礎。在給定的檢測性能指標（包括檢測概率、虛警概率及輸入信號載噪比等）要求下，只有當信號的能量累加到檢測器所需要的信噪比要求時，才能正確捕獲到信號，實現信號的同步。因此，如何高效地實現能量累加，盡量減小累加時間，對減小接收機復雜度，提高捕獲速度具有十分重要的意義。相干與非相干累加是最常用的能量累加方式，特別是對于弱信號的檢測來說，非相干累加更是必不可少的過程。

包絡檢波與平方律檢波是通常采用的兩種非相干檢波方式，文獻［1-3］對這兩種檢波方式進行了理論分析與仿真。其中，文獻［1］給出了平方損耗的簡單計算方法，文獻［2］給出了包絡檢波與理想相干檢波在具有相同累加時間時的處理損耗。根據雷達檢測理論中的定義，文中所指多次非相干累加輸出信噪比實質為在相同檢測概率與虛警概率下采用理想相干累加時所需信噪比，其推導的處理損耗實際上就是N次累加的包絡檢波損耗。文獻［3］對理想檢測能力因子的理解不準確，以致其對平方律檢波處理損耗的計算有錯誤。本文根據雷達檢測理論給出了包絡檢波時的各種檢測性能的理論推導與仿真，并依據平方律檢波時的概率分布，利用仿真軟件得出了在與包絡檢波有相同的檢測概率、虛警概率以及累加次數條件下所需的檢波器前端信噪比及相干累加時間。

2 能量累加原理

在GNSS系統的信號捕獲過程中，需要將接收信號與本地產生的參考信號進行相關運算，其實質是能量的累加過程。其中，相干累加是最有效的能量累加方式。接收到的信號經過與接收機本地產生的載波與偽碼相乘后，設得到的結果為i（n）和q（n），對兩者分別進行積分，得到：

I（n）、Q（n）即為相干累加結果。其中，d（n）為導航數據比特，τ為本地碼與接收碼之間的相位差，R（τ）為碼自相關函數。Tcoh為相干累加時間，Ncoh表示輸入到I和Q支路的積分器的相關結果個數，φd為載波相位。

相干累加時間越長，相干增益越大，但相干累加長度也不是無限度的，它會受到多普勒頻移、導航數據比特寬度和計算量等的限制［4］。為了解決這個問題，通常采用非相干累加的方式。非相干累加即對Nncoh個相干累加的結果進行平方相加，如式（3）所示：

式中，Nncoh為非相干累加個數。相對于相干累加，這種累加方式存在非相干損耗（即包絡檢波損耗）。但非相干累加不受符號跳變和載波相位的影響，可以在多個導航電文數據位上進行能量累加。但是非相干累加的長度也受碼多普勒的限制，而且隨著非相干累加次數增加，能量損耗也會隨之增大，累加的效果會越來越差。很多時候，為了減小硬件實現的復雜度，會采用以下平方律檢波方式代替包絡檢波，可省去開平方的操作，如式（4）所示：

3 檢測性能分析

根據雷達檢測理論，理想相干累加時，信號加噪聲服從正態分布［5］。由其積分可得理想相干檢測時系統的虛警概率Pfa和檢測概率Pd：

式中，Es為信號幅度，S為信號功率，N為噪聲功率，Etc為檢測門限。由此可得在給定虛警概率與檢測概率的情況下，理想相干檢測的檢測因子（也即所需的信噪比）Dc() 1為

3．1 包絡檢波

實際接收機通常采用非相干包絡檢波或平方律檢波以消除載波相位模糊所帶來的問題。在單次積分累加的情況下，僅有噪聲時的包絡檢測輸出為瑞利（Rayleigh）分布，信號加噪聲的包絡檢測輸出為萊斯（Ricean）分布。設門限電平為Et1，則其虛警概率為

由此得檢測門限為

則檢測概率Pd為

于是，包絡檢波時一次非相干累加的可檢測因子（所需信噪比）D0() 1為

由此，可得包絡檢波相對于理想相干檢波的損耗（檢波器損耗）Cx（1）［6］為

Marcum推導出了n次非相干累加后的精確解［7］，并可由此得出檢測因子（所需信噪比）D0（n）。此時，檢波器損耗Cx（n）定義為

式中，D0（n）表示n次累加中每次累加的可檢測因子（每段相干積分后的信噪比）。于是積累損耗

積累增益

D0（n）可由式（14）近似求得：

檢測概率Pd=0.85，虛警概率Pfa=10-6時，Gi（n）與累加次數n的關系如圖1所示。

3．2 平方律檢波

對于平方律檢波，當信號存在時，Nncoh個I2+ Q2的和的概率密度函數是自由度為2Nncoh的中心χ2分布；當信號不存在時，概率密度函數是自由度為2Nncoh的非中心χ2分布［3］。在H0（僅有噪聲）和H1（信號+噪聲）假設下的統計判決變量z的概率密度函數為

式中，λ=2Nncoh×RSNR為非中心參量，Γ()·為伽馬函數，INncoh-1（·）為第一類Nncoh-1階修正貝塞爾函數。設判決門限為Etn，則虛警概率Pfa為

檢測概率Pd為

同樣，當給定虛警概率與檢測概率后，可由虛警概率公式求出檢測門限Etn，然后由檢測概率公式求出平方律檢波前所需的信噪比。但是，由于上述公式求逆較為困難，因此，我們在實際工程設計中可以采用仿真計算的方法，即利用仿真軟件做出一系列需要的圖表或曲線，再通過查表法得到結果。

以下對這一過程進行仿真，仿真結果如圖2～3所示。

從圖2可以看出，采用平方律檢波，當非相干累加次數為4時，若要使虛警概率達到10-6，則檢測門限應為43。從圖3可以得出，當門限為43、累加次數為4時，要使檢測概率達到85%，則輸入信噪比應至少為8 dB。其它情況可依此類推，可得出給定條件下所需要的輸入信噪比。由于上述平方律檢波的概率密度積分式求逆較為困難，因此，這種方法用于系統設計時計算最小輸入信噪比，確定相干積分時間，還是比較方便的。通過此過程，可以把仿真結果與上述包絡檢波方法所得結果進行比較，比較結果如表1所示。表中，RSNRenv表示包絡檢波時所需的信噪比，RSNRsqrt表示平方律檢波時所需的信噪比，RSNRsqrt-RSNRenv表示兩者的差值。

對比表1兩種不同檢波方式所需的信噪比可得如下結論：在同樣的條件下，采用平方律檢波比采用包絡檢波所要求的輸入信噪比要稍大。但采用包絡檢波的復雜度明顯大于平方律檢波，因為其需要開平方操作，在硬件實現中較為困難。因此，實際工程中更適合采用平方律檢波。

4 相干累加時間與累加次數的選取

在工程中，常常涉及到在給定檢測概率與虛警概率以及需要達到的檢測能力（接收信號載噪比）的條件下，如何設計信號捕獲中能量累加方式的問題。本節將對這一問題進行分析與仿真。

4．1 給定累加次數，選擇相干累加時間

在給定非相干累加次數的條件下，我們可以按照如下步驟來選取相干累加時間：

（1）由給定檢測概率Pd與虛警概率Pfa，根據式（15）求出不同的累加次數下包絡檢波所需的檢波前信噪比RSNR或D0（n）（平方律檢波時用式（17）、式（18）計算）；

（2）根據相干積分的輸出信噪比RSNR與積分時間Tcoh和多普勒頻偏fd的關系求出所需的相干積分時間Tcoh，如式（19）所示：

（3）總的積分時間Ttotle=Tcoh×Nncoh，Nncoh為非相干累加次數。

當多普勒頻偏fd=0時，相干積分時間Tcoh= RSNR／（C／n0）=D0（n）／（C／n0），總積分時間Ttotle= Nncoh×D0（n）／（C／n0）。由于Nncoh越大，D0（n）越小，單次相干累加的時間就越短，但因為每次檢波（包括包絡檢波或平方律檢波）都會有損耗，累加次數越多，損耗越大，故總的累加時間會增加，以下仿真說明了這一點。例如，在fd=0，C／n0=27 dB·Hz的情況下，要達到虛警概率Pfa=10-6，檢測概率Pd=0.9，進行4次、8次和16次非相干累加的時間分別為50 ms、59 ms、73 ms，結果如圖4所示。因此，在可能的條件下，盡量減少累加次數，更有利于能量的累加，由此可以減小累加信號長度，提高捕獲速度。

當fd≠0時，信噪比RSNR不再是積分時間Tcoh的線性函數，其關系如圖5和圖6所示。

由圖5和圖6可知，當C／n0=45 dB·Hz時，要達到8 dB的信噪比，在多普勒為1 000 Hz以內時，0.25 ms到0.5 ms之間的相干累加時間都可以達到要求，但最多不能超過5 ms，否則將會導致累加能量的下降，反而達不到所需的信噪比；若C／n0= 27 dB·Hz，當fd=0時，大約需要12.5 ms相干積分，當fd=18 Hz時，大約需要18 ms的相干積分時間，當fd＞18 Hz時，通過延長相干累加時間已經不能達到8 dB的信噪比了。

4．2 給定相干累加時間，確定累加次數

由于非相干累加次數越多，積累損耗越大，因此，應該盡量減小非相干累加次數，增加每次相干累加時間。在實際工程中，往往更容易首先根據信號特點、計算量大小以及系統可以容忍的最大多普勒頻偏來確定最長的相干累加時間，而不是首先確定非相干累加次數。在這種給定最長相干累加時間的情況下，可以根據以下步驟來選取非相干累加次數：

（1）根據式（12）求出在給定Pfa和Pd條件下的D0（1）；

（2）根據式（19）求出在給定Tcoh、fd和C／n0的情況下能達到的相干累加輸出信噪比RSNR或D0（n）；

（3）對于包絡檢波，根據Gi（n）=D0（1）-D0（n）和圖1得出最小累加次數Nncoh；

（4）對于平方律檢波，根據Gi（n）≈D0（1）-D0（n）+0.5和圖1得出最小累加次數Nncoh；

（5）總的累加時間Ttotle=Tcoh×Nncoh。

例如，設需要捕獲C／n0=27 dB·Hz的信號，要求達到虛警概率Pfa=10-6，檢測概率Pd=0.85。利用以上兩種方法確定相干累加時間與非相干累加次數：相干累加時間Tcoh=20 ms，平方律檢波，非相干累加次數Nncoh=4。捕獲仿真結果如圖7所示。由圖可知，此時可很好地捕獲到信號的相關峰。

5 總結

針對現有文獻對GNSS系統擴頻信號捕獲中的能量累加方式與參數選擇未有系統的分析，且部分文獻對概念理解的錯誤，本文利用雷達檢測原理詳細分析了擴頻信號捕獲中非相干能量累加時的兩種檢波方式，根據這兩種檢波方式輸出信號的分布特性推導了各自的檢測概率與虛警概率，并以此為基礎仿真得出兩種檢波方式下檢測性能與不同的參數之間的關系，給出了在相同條件下兩種檢波方式的檢測性能的差異。針對工程應用實際，給出了能量累加過程中相干累加時間與累加次數的選擇方法，并進行了舉例分析，這些方法對擴頻信號接收機的設計具有重要的指導作用。

［1］James Bao-Yen Tsui.GPS軟件接收機基礎［M］.2版.陳軍，潘高峰，李飛，等，譯.北京：電子工業出版社，2007：198-199.

James Bao-Yen Tsui.Fundamentals of Global Positioning System Receivers A Software Approach［M］.2nd ed.Translated by CHEN Jun，PAN Gao-feng，LI Fei，et al.Beijing：Publishing House of Electronics Industry，2007：198-199.（in Chinese）

［2］王飛雪，郭桂蓉.二相編碼信號分段相關-視頻積累檢測的最優中頻積累時間［J］.國防科技大學學報，1999，21（1）：71-75.

WANG Fei-xue，GUO Gui-rong.The Optimum IF Integra ion Time Length for Binary Phase-coded Signals Detect ion Based on Segment Correlation-Video Integration Method［J］.Journal of National University of Defense Technology，1999，21（1）：71-75.（in Chinese）

［3］朱祥維，王飛雪.基于平方律檢波的二相編碼信號分段相關-視頻積累方法研究［J］.電子學報，2005，33（3）：545-548.

ZHU Xiang-wei，WANG Fei-xue.Study on Segment Correlation-Video Integration Method Using Sqare-law Detector for Binary Phase-coded Signals Detection［J］.Acta Electronica Sinica，2005，33（3）：545-548.（in Chinese）

［4］謝鋼.GPS原理與接收機設計［M］.北京：電子工業出版社，2009：360-362.

XIE Gang.Principles of GPS and Receiver Design［M］.Beijing：Publishing House of Electronics Industry，2009：360-362.（in Chinese）

［5］David K Barton.雷達系統分析與建模［M］.北京：電子工業出版社，2007：29-38.

David K.Barton.Radar System Analysis and Modeling［M］. Beijing：Publishing House of Electronics Industry，2007：29-38.（in Chinese）

［6］ Barton D K.Simple Procedures for Radar Detection Calculation［J］.IEEE Transactions on Aerospace and Electronic Systems，1969，5（5）：837-846.

［7］ Marcum F I.A Statistical Theory of Target Detection by Pulsed Radar［J］.Rand Research Memo，1948，754：1211-1215.

LUO Kui was born in Mianyang，Sichuan Province，in 1984. He received the B.S.degree in 2007.He is now a graduate student.His research interests include satellite communication and satellite navigation.

Email：luo-83＠126.com

程劍（1976—），男，湖南醴陵人，2005年獲博士學位，現為解放軍理工大學通信工程學院副教授，主要研究方向為衛星通信抗干擾、衛星導航；

CHENG Jian was born in Liling，Hunan Province，in 1976.He received the Ph.D.degree in 2005.He is now an associate professor in PLA University of Science and Technology.His research concerns anti -jamming in satellite communication and satellite navigation.

Email：chjatice＠yahoo.com

李仰志（1982—），男，福建龍巖人，2005年獲學士學位，現為解放軍理工大學通信工程學院博士研究生，主要研究方向為衛星通信、衛星導航；

LI Yang-zhi was born in Longyan，Fujian Province，in 1982.He received the B.S.degree in 2005.He is currently working toward the Ph.D.degree in PLA University of Science and Technology.His research interests include satellite communication and satellite navigation.

Email：liyangzhi121＠163.com

賴廣鋒（1984—），男，福建龍巖人，2006年獲學士學位，現為碩士研究生，主要研究方向為衛星通信、衛星導航。

LAI Guang-feng was born in Longyan，Fujian Province，in 1984.He received the B.S.degree in 2006.He is now a graduate student in PLA University of Science and Technology.His research interests include satellite communication and satellite navigation.

Email：

laiguangfeng-88＠163.com

制導與引信信息網11011技術交流會通知

在我國導航、制導與控制技術蓬勃發展之際，為在“十二五”期間進一步推動我國航天技術的進步，促進相關航天應用產業的發展，制導與引信信息網定于2011年6月22日-24日在成都召開“制導與引信11011技術情報交流會”。

本次會議由副網長單位中國航天科工集團二院25所主辦，中國電子科技集團公司第29所、北京高博特廣告有限公司協辦。

有關事項通知祥見廣告（前插3）。

聯系方式：

地 址：北京市海淀區西四環北路15號依斯特大廈709室郵編：100195

電話總機：010-88587506／07 傳真：010-88587416

聯系人：陳運：分機-807 手機：13488876342

孫柏山：分機-813 手機：13641190351

Email：ZDYYX2011＠163.com

Performance Analysis and Parameters Selection of Detection Methods Based on Non-coherent Accumulation

LUO Kui，CHENG Jian，LI Yang-zhi，LAI Guang-feng
（Institute of Communication and Engineering，PLA University Science and Technology，Nanjing 210007，China）

Two detection methods in non-coherent energy accumulation in GNSS（Global Navigation Satellite System）signal acquisition process are discussed，the properties of which are analysed and compared by using radar detection theories，the performance curves of the two methods are plotted by Matlab.Two methods to chose coherent accumulation time and non-coherent accumulation times are proposed and simulated.The methods and simulation results can be used as a reference in GNSS design.The simulation results also show that the envelop detection performs better than square law detection on energy accumulation less than 0.5 dB.

GNSS；signal acquisition；spread spectrum；non-coherent accumulation；envelop detection；square law detection

The National High-tech R＆D Program of China（863 Program）（2008AA12Z307）

TN96

A

10.3969／j.issn.1001-893x.2011.05.012

羅奎（1984—），男，四川綿陽人，2007年獲學士學位，現為碩士研究生，主要研究方向為衛星通信、衛星導航；

1001-893X（2011）05-0056-06

2010-12-06；

2011-03-02

國家高技術研究發展計劃（863計劃）項目（2008AA12Z307）