王星，曹晉龍，＊＊，趙玉，王士巖，2，李承志

(1.空軍工程大學(xué) 航空航天工程學(xué)院，陜西 西安 710038;2.解放軍93286部隊32分隊，沈陽 110141;3.解放軍93363部隊，沈陽 110141)

1 引言

相位編碼是雷達(dá)脈內(nèi)調(diào)制技術(shù)中一種典型的調(diào)制技術(shù)，其中主要的調(diào)制形式是二相編碼(BPSK)和四相編碼(QPSK)。要正確地獲得相位編碼信號的編碼規(guī)律，必須首先知道信號的載頻、碼速率、初始相位、碼同步等各項參數(shù)，其中載頻估計是對其他參數(shù)進(jìn)行估計的基礎(chǔ)。在電子偵察領(lǐng)域，對接收的相位編碼信號的參數(shù)沒有任何先驗知識，對這類信號的截獲、識別和參數(shù)估計存在一定的困難。針對這一問題，文獻(xiàn)［1］提出了一種基于相位展開和最小二乘多項式擬合的算法，此算法在較低信噪比下載頻估計的均方根誤差依然接近相位編碼信號載頻估計的克拉美－羅限(Cramer－Rao)，但是該算法的復(fù)雜度較高，不宜實時處理。文獻(xiàn)［2］給出了一種簡單快速實用的非線性相位編碼信號載頻盲估計算法，但在電子偵察領(lǐng)域，對接收到的信號無法進(jìn)行匹配接收，同時接收機(jī)帶寬又遠(yuǎn)大于信號帶寬，導(dǎo)致整體估計性能下降，甚至使帶外噪聲完全淹沒信號，最終無法對載頻進(jìn)行估計。文獻(xiàn)［3］對MAT算法進(jìn)行了改進(jìn)，提出了M－MAT算法，可以在較低信噪比條件下對QPSK信號的載頻精確估計，但該算法仍然較為復(fù)雜，不易在RWR/ESM系統(tǒng)上實現(xiàn)。

瞬時測頻(IFM)接收機(jī)測頻分辨力高，瞬時帶寬寬，實時性強(qiáng)，被廣泛應(yīng)用于RWR/ESM測頻系統(tǒng)中。但由于傳統(tǒng)IFM實現(xiàn)機(jī)理以及RWR/ESM測頻資源及結(jié)構(gòu)的限制，使得當(dāng)前RWR/ESM系統(tǒng)無法通過IFM對相位編碼信號的檢測進(jìn)行威脅告警，只是將此類信號歸類為未知雷達(dá)信號。IFM若能實時檢測相位編碼信號，將使RWR/ESM系統(tǒng)具備對相位編碼信號進(jìn)行實時告警的能力，因此對此進(jìn)行研究意義重大。為此本文根據(jù)IFM基本原理，應(yīng)用仿真的方法分析BPSK和QPSK信號經(jīng)過IFM系統(tǒng)的輸出，以此為基礎(chǔ)對傳統(tǒng)的IFM系統(tǒng)進(jìn)行了改進(jìn)，并提出估計BPSK和QPSK信號參數(shù)的方法。

2 RWR/ESM測頻系統(tǒng)概述

當(dāng)前RWR/ESM測頻系統(tǒng)多采用引導(dǎo)式測頻體制，以滿足空域頻域?qū)掗_、實時性強(qiáng)和頻率分辨力高的要求，測頻過程可分為粗分頻段、粗測頻、精測頻。

IFM接收機(jī)對雷達(dá)信號的測頻能力直接決定RWR/ESM測頻系統(tǒng)的測頻能力。而由于IFM接收機(jī)在對雷達(dá)信號測頻時，一個脈寬內(nèi)只進(jìn)行一次采樣，因此對脈內(nèi)具有相位編碼特征的雷達(dá)信號是無法測頻的。而當(dāng)前能夠檢測相位編碼信號的相關(guān)算法和數(shù)字接收機(jī)技術(shù)不適用于RWR/ESM系統(tǒng)，主要受以下兩個因素的限制:

(1)當(dāng)前相位編碼信號檢測算法雖然靈敏度較高，且具有良好的復(fù)雜信號適應(yīng)能力，但是算法數(shù)學(xué)運算復(fù)雜，運算量大，導(dǎo)致其信號測量時間較長，實時性差，不能滿足 RWR/ESM系統(tǒng)實時告警的需求;

(2)盡管當(dāng)前較為先進(jìn)的數(shù)字接收機(jī)的ADC已可達(dá)數(shù) Gsample/s采樣率，瞬時帶寬已經(jīng)大于1 GHz(極限值為5 GHz)，但相對于 RWR/ESM系統(tǒng)的瞬時覆蓋帶寬需求(通常數(shù)十GHz)而言，它仍是一個窄帶接收機(jī)，不能滿足RWR/ESM系統(tǒng)寬瞬時帶寬的要求。目前的解決方法是采用多個窄帶數(shù)字接收機(jī)拼接來滿足大的瞬時帶寬，拼接方式主要有時間交替多通道并行采樣與信道化［4］，但按此方法體積、重量、功耗、成本代價都相當(dāng)高［5］。

綜上所述，目前多數(shù)RWR/ESM系統(tǒng)并不具備對相位編碼信號的實時測頻能力。那么，IFM能否在沒有先驗知識的情況下，實時估計相位編碼信號參數(shù)呢?本文通過對IFM技術(shù)的仿真分析，發(fā)現(xiàn)其具備這一潛力。

3 相位編碼信號的IFM仿真分析

3.1 IFM接收機(jī)原理簡介

典型IFM單元為微波鑒相器(Microwave Phase Discriminator，MPD)，如圖1 所示，它由功分器、延遲線、90°電橋、平方率檢波器和差分放大器組成［6］。

圖1 實用微波鑒相器結(jié)構(gòu)示意圖Fig.1 Block diagram of practical MPD

其中，相角如公式(2)所示:

由式(1)、(2)得信號頻率如公式(3)所示:

式(1)～(3)中，K為檢波器系數(shù)，A為信號幅度，T為延遲線的時間延遲，f為輸入信號的載頻。為了解決頻段覆蓋和分辨率的矛盾，工程上采用圖2所示IFM并聯(lián)結(jié)構(gòu)［7］，其多個延遲線長度符合一定關(guān)系，短延遲線支路保證測頻范圍，長延遲線保證精度;并根據(jù)長短延遲線之間的關(guān)系解IFM頻率測量值的模糊。

圖2 傳統(tǒng)IFM系統(tǒng)Fig.2 Traditional IFM system

若設(shè)圖1端口1進(jìn)入一個BPSK或QPSK信號，則UI、UQ的數(shù)學(xué)表達(dá)式較為復(fù)雜，下面應(yīng)用MALAB仿真的方式對其進(jìn)行分析。

3.2 BPSK信號的IFM仿真

BPSK信號一般模型［8］表示為

式中，q(t)為雙極性偽隨機(jī)序列，A為信號幅度，Tc為chip時寬，f0為載頻，φn為第n個碼元對應(yīng)的相位，等概率時取0和π。

以13位巴克碼序列［1111100110101］構(gòu)造一個BPSK 信號 s(t)，A=1 V，Tc=5 μs，τ=65 μs。為便于分析，載頻(設(shè)為 MHz級別)f0=2 MHz。應(yīng)用MATALB進(jìn)行仿真，建立三路延遲時間分別為 0.1 μs、0.4 μs、1.6 μs的并行 MPD(Maintenance Planning Document)模型。BPSK信號s(t)一個脈寬內(nèi)結(jié)果如圖3所示。其進(jìn)入MPD模型中，輸出的UI、UQ值如圖4所示。

圖3 BPSK信號Fig.3 BPSK signal

圖4 輸入為BPSK信號時MPD的UI、UQFig.4 UI、UQof MPD with BPSK signal

3.3 QPSK信號的IFM仿真

QPSK信號可以視為兩路正交的BPSK信號的疊加，其一般模型［8］可表示為

其中，q1(t)、q2(t)為兩路的雙極性偽隨機(jī)序列，φ1_n、φ2_n為q1(t)、q2(t)的第n個碼元所對應(yīng)的相位。

設(shè)q1(t)、q2(t)兩個序列分別為［1100001100］、［1111001111］，A=1 V，Tc=5 μs，f0=2 MHz，脈沖寬度τ=50μs。一個脈寬內(nèi)的仿真結(jié)果如圖5所示。此QPSK信號s(t)同樣進(jìn)入上述三路并行MPD模型中，輸出的UI、UQ值如圖6所示。

圖5 QPSK信號Fig.5 QPSK signal

3.4 仿真結(jié)果分析

由以上仿真可知，BPSK和QPSK信號經(jīng)過MPD后的UI、UQ為恒定電壓值，而在s(t)的相位變化處存在跳變，跳變的持續(xù)時間為MPD的延遲時間。BPSK信號的UI、UQ跳變幅度恒定，而QPSK信號的UI、UQ跳變幅度變化。由此可見，這些相位跳變是存在一定規(guī)律的，若掌握了其中的規(guī)律，即可對BPSK和QPSK信號的相關(guān)參數(shù)進(jìn)行估計。但傳統(tǒng)的IFM系統(tǒng)對一個脈寬的信號只進(jìn)行一次采樣，無法得到跳變信息，因此需對傳統(tǒng)的IFM系統(tǒng)進(jìn)行改進(jìn)。

圖6 輸入為QPSK信號時MPD的UI、UQFig.6 UI、UQof MPD with QPSK signal

4 IFM系統(tǒng)改進(jìn)

4.1 改進(jìn)的IFM系統(tǒng)概述

由以上分析可知，可用ADC代替?zhèn)鹘y(tǒng)IFM系統(tǒng)中的極性量化器，對MPD輸出的UI、UQ進(jìn)行連續(xù)時間采樣，檢測到相位變化處的跳變信息。根據(jù)跳變信息對BPSK和QPSK信號進(jìn)行區(qū)分識別，采用相關(guān)算法對離散的UI、UQ進(jìn)行數(shù)字解算，進(jìn)而估計出BPSK和QPSK信號的相關(guān)參數(shù)。

改進(jìn)后的IFM如圖7所示，具有如下兩個典型特征:

(1)采用高精度ADC代替量化器結(jié)構(gòu)，通過數(shù)字解算不僅可解模糊，且可以獲得更多信號特征;

(2)信號檢波的包絡(luò)經(jīng)過整形后，作為信號有效信號。在其有效持續(xù)時間內(nèi)，用ADC對UI、UQ進(jìn)行連續(xù)時間采樣，對采樣結(jié)果進(jìn)行結(jié)算。

圖7 改進(jìn)后的IFM系統(tǒng)Fig.7 The modified IFM system

4.2 測頻分辨力分析

由于改進(jìn)IFM系統(tǒng)是對UI、UQ單獨采樣量化，而不是采用傳統(tǒng)IFM系統(tǒng)中的極性量化，因此可對式(3)求全微分，得改進(jìn)后一路鑒相器的頻率分辨力如公式(6)所示:

式中，ΔUI、ΔUQ是 ADC 對 UI、UQ的量化間隔。圖7所示四路鑒相器并行運用的改進(jìn)后IFM系統(tǒng)，其頻率分辨力如公式(7)所示:

5 BPSK、QPSK信號參數(shù)估計

5.1 BPSK、QPSK 信號檢測識別

估計BPSK和QPSK信號參數(shù)的前提是必須區(qū)分這兩種信號。根據(jù)文獻(xiàn)［9］可知，利用不同信號的相位差分變化規(guī)律可以有效地識別BPSK和QPSK信號。本文MPD輸出的相位是指由延遲線引起的相位差，這與文獻(xiàn)［9］中的差分相位本質(zhì)是一樣的，不同的是文獻(xiàn)［9］中的差分相位是應(yīng)用微分的方法求得，是理論方法，而本文中相位差是應(yīng)用IFM鑒相的方法求得，是一種硬件實現(xiàn)方法，更具有實際意義。下面具體分析BPSK和QPSK的識別方法。

BPSK和QPSK信號經(jīng)過MPD后的輸出相位結(jié)果在幅度上明顯不同，前者的輸出結(jié)果幅度只跳變一次，后者輸出結(jié)果幅度跳變?nèi)巍＿@是因為BPSK信號碼元對應(yīng)相位φn為兩個值，一般為0和π，s(t)的相位跳變一次，跳變值為π;而QPSK信號的碼元對應(yīng)相位φn為4個值，假設(shè)φn=［0π/2 π 3π/2］，s(t)的相位跳變?nèi)?，跳變值為?2、π和3π/4。根據(jù)上述規(guī)律，可設(shè)定門限值對BPSK和QPSK信號進(jìn)行區(qū)分。若MPD的輸出結(jié)果等幅跳變，則可判定輸入為BPSK信號;若MPD輸出結(jié)果幅度跳變?nèi)?，則可判定輸入為QPSK信號。區(qū)分BPSK和QPSK信號后，可對信號的載頻、碼元及碼周期進(jìn)行估計。

5.2 載頻估計

為模擬實際情況，在仿真中對BPSK和QPSK信號加入信噪比為10dB噪聲。用采樣頻率為10MHz的 ADC 對 MPD1、MPD2、MPD3輸出的UI、UQ值進(jìn)行采樣量化，由公式(3)可得離散頻率值，應(yīng)用解模糊算法［10］，消除長延遲線支路的模糊，并剔除相位跳變點得結(jié)果如圖8所示。對第三路輸出的精度較高的離散頻率值取平均，BPSK信號載頻估計=1.9979 MHz，QPSK信號載頻估計=2.0012 MHz。

圖8 頻率輸出Fig.8 QPSK signal

改變信噪比，仿真500次的均方根誤差結(jié)果如圖9所示。由圖可知，隨信噪比的增加，誤差逐漸減小。BPSK信號的誤差曲線比QPSK信號的誤差曲線要平滑，這是因為BPSK信號在測頻時所產(chǎn)生的跳變點恒定，而QPSK信號的跳變點不恒定。當(dāng)信噪比大于5 dB時，測頻的均方根誤差收斂于 0.05 MHz。

圖9 測頻誤差與信噪比關(guān)系Fig.9 Relationship between SNR and error frequency measurement

本文只對三路MPD進(jìn)行了仿真，這是由于相位跳變的持續(xù)時間與延遲線產(chǎn)生的延時時間相等，按各路MPD延時間的倍數(shù)關(guān)系，第四路MPD的延時為6.4 μs，此時相位跳變的持續(xù)時間大于碼元寬度，UI、UQ會發(fā)生嚴(yán)重失真，使載頻無法估計。而實際雷達(dá)信號的載頻為幾GHz至幾十GHz之間，因此MPD的延遲時間為ns數(shù)量級別。典型的四路MPD的延遲時間為 0.15 ns、0.6 ns、2.4 ns、9.6 ns，一般遠(yuǎn)小于碼元寬度，因此不會造成UI、UQ的失真。此時第四路會輸出更高精度的測量結(jié)果，一般比第三路高一個數(shù)量級。因此，信噪比大于5 dB時可到達(dá)當(dāng)前IFM的精度要求。

5.3 碼元和碼周期估計

BPSK信號的碼元只有0和1，因此可根據(jù)UI、UQ的幅度跳變信息直接估計得到，即相位跳變判別法。每當(dāng)相位幅度跳變一次，碼元相應(yīng)變化一次，跳變的最小時間間隔即為碼周期。

QPSK信號的碼元估計可采兩種方法。

一種是如上述的相位跳變判別法。QPSK信號的UI、UQ的幅度跳變信息有正值和負(fù)值，不易找出其變化規(guī)律。因此對UI、UQ取反正切，得相位值如圖10所示，此處相位值是由延時線引起的相位差。由于在解模糊的過程中，相位跳變點會被剔除，因此圖10是未解模糊的結(jié)果。

圖10 QPSK的相位差Fig.10 Phase difference of QPSK

下面就其相位跳變規(guī)律進(jìn)行分析。QPSK信號的碼元有4個:00，01，11和10。假設(shè)其相位分別為φn=［0，π/2，π，3π/2］，那么碼元之間的變化規(guī)律與相位跳變之間的變化規(guī)律的關(guān)系如圖11所示，碼元從00到01，相位跳變π/2;碼元從00到11，相位跳變π;碼元從00到10，相位跳變3π/2。因此可得出如下規(guī)律:相鄰碼元之間的相位跳變值最小，間隔兩個碼元之間的相位跳變值最大，間隔一個碼元之間的相位跳變值在上述兩個之間。按此規(guī)律分析圖10中MPD3輸出相位，假設(shè)初始碼元為11，在10μs處相位跳變最小，碼元由11變?yōu)?0;在15 μs處相位跳變最大，碼元由10變?yōu)?1;在20μs處相位跳變在最大值與最小值之間，碼元由11變?yōu)?0。同理，在30μs和40μs處相位也有跳變，碼元依次為11和10。根據(jù)相位跳變的最小時間間隔可得碼周期為5 μs，由此可估計出碼元為［11111011000011111010］。其中碼元高位對應(yīng)q2(t)的估計值，為［1111001111］;低位對應(yīng) q1(t)估計值，為［1100001100］。

圖11 相位跳變與碼元關(guān)系圖Fig.11 Relationship between symbol and phase jump

另一種方法是應(yīng)用相干解調(diào)法，此方法的前提是必須已知信號載頻，而在RWR/ESM系統(tǒng)中，信號載頻屬未知參數(shù)，因此可用載頻估計值^fQPSK代替原始信號載頻，實際應(yīng)用中可用窄帶接收機(jī)輸出的精測頻數(shù)據(jù)對本振進(jìn)行引導(dǎo)。碼元估計原理如圖12所示。

圖12 相干解調(diào)法原理圖Fig.12 Schematic diagram of coherent demodulation

根據(jù)原理框圖進(jìn)行建模仿真，結(jié)果如圖13所示，圖中碼元估計的延遲是由于濾波延遲所造成的。

圖13 相干解調(diào)法估計的QPSK信號碼元Fig.13 QPSK signal symbol of coherent demodulation

綜合上述分析，兩種方法各有其優(yōu)缺點:相位跳變判別法較為簡便，實時性好，適用于RWR/ESM系統(tǒng)，但碼元初始值未知則無法完成估計，因此需要一定的先驗知識才能估計出碼元;相干解調(diào)法不需要任何先驗知識，但對信號載頻的估計精度要求較高，且實時性不如相位跳變判別法，因此這種方法適用于電子情報偵察，得到碼元估計值可作為建立碼元數(shù)據(jù)庫的基礎(chǔ)。

6 結(jié)束語

本文應(yīng)用仿真的方法分析了IFM系統(tǒng)對相位編碼信號(二相編碼和四相編碼信號)處理情況，提出了基于IFM的相位編碼信號的載頻和碼元估計算法，其原理簡單，計算量小，實時性好，適用于RWR/ESM系統(tǒng)中。本文的研究為今后RWR/ESM系統(tǒng)測頻接收機(jī)的研制和改進(jìn)提供了一種新的思路，有較好的實際應(yīng)用價值。

［1］鄧振淼，劉渝.MPSK信號載頻盲估計［J］.通信學(xué)報，2007，28(2):94－100.DENG Zhen－miao，LIU Yu.Blind estimation of MPSK carrier frequency［J］.Journal on Communications，2007，28(2):94－100.(in Chinese)

［2］Ghogho M，Swami A，Durrani T.Blind estimation of frequency offset in the presence of unknown multipath［C］//Proceedings of 2000IEEE International Conference on Personal Wireless Communications. Piscataway， NJ:IEEE，2000:104－108.

［3］朱霞，劉渝，狄慧.一種低信噪比條件下QPSK信號盲處理方法［J］.數(shù)據(jù)采集與處理，2011，26(5):555－558.ZHU Xia，LIU Yu，DI Hui.Blind Processing Met－h(huán)od for QPSK Signal under Low SNR［J］.Journal of Data Acquisition ＆ Processing，2011，26(5):555－558.(in Chinese)

［4］王洪.寬帶數(shù)字接收機(jī)關(guān)鍵技術(shù)研究及系統(tǒng)實現(xiàn)［D］.成都:電子科技大學(xué)，2007.WANG Hong.Research and Implementation of system on Wide band Digital Receiver and Key Technologies［D］.Chengdu:University of Electronic Science and Technology of China，2007.(in Chinese)

［5］王坤達(dá).基于實時數(shù)字信號處理的寬帶單比特瞬時測頻接收技術(shù)［J］.艦船電子對抗，2012，35(3):31－36.WANG Kun－da.Wide band single－bit IFM receiving technology based on real－time digital signal processing［J］.Shipboard Electronic Countermeasure，2012，35(3):31－36.(in Chinese)

［6］王星.航空電子對抗原理［M］.北京:國防工業(yè)出版社，2008.WANG Xing.Principles of aircraft electronic countermeasure［M］.Beijing:National Defense Industry Press，2008.(in Chinese)

［7］張學(xué)成.寬帶大動態(tài)瞬時測頻接收機(jī)的設(shè)計與實現(xiàn)［J］.艦船電子對抗，2011，34(4):29－31.ZHANG Xue－cheng.Design and realization of wide band instantaneous frequency measurement receiver with large dynamic range［J］.Shipboard Electronic Countermeasure，2011，34(4):29－31.(in Chinese)

［8］Skolnik M L.雷達(dá)手冊［M］.南京電子技術(shù)研究所，譯.北京:電子工業(yè)出版社，2010.Skolnik M L.Radar Handbook［M］.Translated by Nanjing Institute of Electronic Technology.Beijing:Publishing House of Electronics Industry，2010.(in Chinese)

［9］黃知濤，周一宇.一種有效的BPSK/QPSK信號調(diào)制識別方法［J］.電子對抗技術(shù)，2005，20(2):10－13.HUANG Zhi－tao，ZHOU Yi－yu.An effective approach to the recognition of BPSK/QPSK signals［J］.Electronic Warfare Technology，2005，20(2):10－13.(in Chinese)

［10］唐永年.雷達(dá)對抗工程［M］.北京:北京航空航天大學(xué)出版社，2012:195－198.TANG Yong－nian.Engineering of radar counter measure［M］.Beijing:Publishing House of Beijing University of Aeronautics and Astronautics，2012:195－198.(in Chinese)