苗雪平，張林讓，雷 宇

(西安電子科技大學雷達信號處理國家重點實驗室，陜西西安 710071)

三角調頻脈沖信號是現(xiàn)代雷達電子戰(zhàn)中常采用的一種發(fā)射信號，其具有對多普勒頻移不敏感性，解決信號能量和距離分辨率矛盾的優(yōu)點。測量三角調頻脈沖信號參數(shù)是雷達信號處理工程領域的重點之一。傳統(tǒng)的測量方法僅針對時域或調制域，未將時域和調制域相結合進行參數(shù)測量。且測量三角調頻信號調制域參數(shù)需進行時頻分析，常用的時頻分析方法，如短時傅里葉變換(STFT)和Wigner－Vile分布，其在對三角調頻信號時頻分析時存在著諸多缺陷［1］。STFT存在著窗長和窗函數(shù)選擇困難的問題，由于使用了傅里葉變換，運算量較大，且受測不準準則的限制，不能解決時間分辨率和頻率分辨率的矛盾。Wigner－Vile分布雖不使用窗函數(shù)，不存在STFT的問題，但對多分量信號進行時頻分布時存在著“交叉項”干擾，時頻分布將變得模糊，對頻率估計不準，且其采用了積分運算，運算量較大。

將三角調頻脈沖信號的時域和調制域參數(shù)測量相結合，時域分析時采用統(tǒng)計脈沖信號包絡數(shù)值的方法估計判決門限，進行脈沖個數(shù)、脈沖寬度等時域參數(shù)的測量。在時域測量結果的基礎上，對三角調頻信號的脈內調制參數(shù)進行測量。脈內調制參數(shù)測量采用相位差分法來得到信號的瞬時頻率曲線，并進行線性擬合來獲得線性度等調制域參數(shù)。其中，相位差分法首先通過反正切函數(shù)計算信號的瞬時相位，再對瞬時相位進行差分運算，直接得到信號的瞬時頻率，相對于STFT，運算量大幅減小，且無窗函數(shù)，也不受測不準準則的限制。此外，相位差分法不存在Wigner－Vile分布的“交叉項”干擾問題。最后進行了Matlab仿真，仿真表明，高信噪比下該方法相對于其他方法，測量參數(shù)運算量少、誤差小，適合于工程實踐。

1 三角調頻脈沖信號模型

三角調頻脈沖信號是一種脈內頻率調制信號，其指的是在脈沖寬度內頻率連續(xù)線性變化，信號頻率先向上后向下或先向下后向上掃描脈沖寬度。其信號表達式為

式中，rect(t)為矩形脈沖;A為三角調頻信號的幅度;f0為中心頻率;μ為調頻系數(shù);f0+μt或 f0－μt為瞬時頻率。

2 脈沖時域參數(shù)測量

脈沖時域參數(shù)是描述雷達信號特征的重要參數(shù)，主要包括脈沖寬度、脈沖周期、上升時間、下降時間、關閉時間、占空比、脈沖幅度、幅度90%、幅度10%、頂值、底值、過沖、下沖等，各參數(shù)定義［2］。

脈沖時域參數(shù)的測量以脈沖包絡為基礎，設脈沖實、虛部信號分別為I(t)和Q(t)，計算信號包絡的表達式為

以信號包絡y(t)來估計脈沖的近似幅度，關鍵要估計脈沖波形的近似頂值和近似底值，以獲得判決門限。估計算法主要有兩種:(1)密度分布平均法［3］;(2)密度分布眾數(shù)法［3］。

文獻［2］指出，當測量采樣點數(shù)越多，在脈沖包絡最大值和最小值之間均分的區(qū)間個數(shù)越多，所求得的近似底值和近似頂值越精確;通常密度分布平均法能比密度分布眾數(shù)法得到更好的精度，但其運算量也偏大。

2.1 脈沖個數(shù)測量

估計出脈沖的近似頂值 topvalue'和近似底值basevalue'后，近似幅度ampvalue'表達式為二者之差。

以10%ampvalue'作為上升沿下降沿檢測門限，測量脈沖個數(shù)的步驟為:

(1)脈沖個數(shù)初始化為N=0，從信號的起始位置開始查找第一個近似上升沿，查找成功標志位為flag=1，否則標志位為flag=0。

(2)查找到近似上升沿后，繼續(xù)往后查找近似下降沿，查找成功標志位為flag=1，否則標志位為flag=0。

(3)查找到近似下降沿后，繼續(xù)往后查找近似上升沿，查找成功后脈沖個數(shù)N=N+1，并將標志位賦值flag=1，否則標志位flag=0。

(4)重復步驟(2)至步驟(3)直到標志位flag=0，則結束。

2.2 脈沖寬度測量

在測量脈沖個數(shù)過程中，記錄所有近似上升沿、近似下降沿判決時刻，對信號的某個脈沖，近似上升沿判決時刻為t'升，近似下降沿判決時刻為t'降，則計算脈沖寬度表達式為

2.3 脈沖其他參數(shù)測量

在每個脈沖的脈沖寬度內，運用密度分布平均法，計算每個脈沖精確的頂值top和底值base。并求得每個脈沖精確幅度amp為二者之差。

求出每個脈沖包絡的最大值max和最小值min，則過沖overshoot為最大值max和頂值top之差，下沖undershoot為底值base和最小值min之差，幅度的90%(記為amp90%)和幅度的10%(記為amp10%)表達式為

然后通過幅度的90%和幅度的10%這兩條脈沖幅度參考線，在脈沖寬度τ內查找這兩條參考值所對應的時刻值，即可計算出脈沖的上升時間tr、下降時間tf，脈沖周期T，關閉時間toff及占空比dt等參數(shù)。其計算表達式為

其中，tr90是上升沿幅度90%所對應的時刻;tr10是上升沿幅度10%所對應的時刻;t't10是相鄰下一次上升沿幅度10%所對應的時刻;tf10是下降沿幅度10%所對應的時刻;tf90是下降沿幅度90%所對應的時刻。由脈沖寬度τ和脈沖周期T，可計算關閉時間toff和占空比dt。

3 脈內調制參數(shù)測量

基于已估計出三角調頻脈沖信號的脈沖寬度及脈沖周期，可測量其脈內調制參數(shù)，主要包括分析調頻方向和測量最小頻率、最大頻率、中心頻率、帶寬、左調頻系數(shù)、右調頻系數(shù)、左線性度和右線性度。首先，應提取三角調頻信號每個脈沖的頻率，即要對信號進行時頻分析。常用的時頻分析方法有短時傅里葉變換方法(STFT)、Wigner－Vile分布方法和相位差分法等，本文采用相位差分法。

3.1 短時傅里葉變換

STFT利用一個時間有限的窗函數(shù)來截取信號，然后對截取信號采用傅里葉變換進行分析，以確定在這一窗時間段內存在的頻率，再沿著時間軸移動窗函數(shù)，最終得到信號的瞬時頻率隨時間的變化關系圖。研究表明STFT算法簡單，且為線性時頻分布，不會在頻率分量之間產生交叉項干擾，但采用傅里葉變換，需多次進行積分運算，運算量較大，且由于窗函數(shù)的影響，窗長和窗函數(shù)選擇困難，也無法克服時間分辨率與頻率分辨率之間的矛盾［4］。因此，STFT不適合工程上三角調頻信號時頻分析。

3.2 Wigner－Vile分布

Wigner－Vile分布是實函數(shù)，設某實信號的解析形式為x(t)，則其Wigner－Vile分布為

從上式可看出，Wigner－Vile分布中x(t)出現(xiàn)了兩次，且不含任何窗函數(shù)，是雙線性變換，不存在時間分辨率與頻率分辨率的矛盾，故具有較好的時頻聚集性。但其不滿足可加性，且為了得到時頻曲線，需進行多次積分運算，運算量較大。經研究表明，其還會在時頻平面上虛假分量的位置產生分量之間的交叉項干擾。交叉項會產生有限的能量，其反映了兩個相關項之間的相關性并高度震蕩，難以克服［5－6］。因此，該方法不適與多分量信號三角調頻信號的時頻分析。

3.3 相位差分法

假定單一分量連續(xù)時間信號s(t)解析表達式為

對于采集的離散信號可通過信號的相位差分來估計信號的瞬時頻率。采用中心有限差分估計信號的瞬時頻率fc(n)表達式為

其中，fs為信號的采樣率;φ(n)為第n點的瞬時相位。

根據中心有限差分估計的特點，可看出其對信號的頻率是無偏估計，且群時延為零，對應時頻分布的一階矩。相對于STFT，算法簡單，無多次復雜的積分運算，僅有一次反正切運算和差分運算，運算量大幅減小，且無窗函數(shù)，也不受測不準準則的限制。此外，相位差分法不存在Wigner－Vile分布的“交叉項”干擾問題，后續(xù)的仿真結果將對此驗證。故利用中心有限差分估計三角調頻信號頻率是有效的。

由瞬時頻率曲線可得出調頻方向。根據調頻方向，可確定頻率曲線中心為頻率最大值或是最小值。以頻率曲線為中心將曲線分為左、右兩半部分，分別進行直線擬合，求得左、右線性度和左、右調頻系數(shù)。求得脈沖寬度τ，瞬時頻率最大值freqMax，頻率最小值freqMin。帶寬B為二者差值，中心頻率f0為二者之和的

4 仿真結果

為驗證時域測量方法的有效性，以三角調頻脈沖信號進行Matlab仿真，參數(shù)設置為調頻方式為向上調頻，脈沖寬度為6μs，脈沖重復周期為16μs，完整脈沖個數(shù)為4，采樣率為50 MHz，帶寬10 MHz，信噪比為20 dB。根據上述方法，測得的時域參數(shù)如表1所示。

表1 三角調頻脈沖信號時域參數(shù)表

從表1中可看出，測量的脈沖個數(shù)為4，脈沖寬度約為6.1μs，重復周期為16μs，上升時間為0.04μs，下降時間為0.02μs，與實際設定的參數(shù)基本上一致，測量誤差小。用STFT、Wigner－Vile方法分析三角調頻脈沖信號，時頻分布如圖1所示。

圖1(a)STFT采用窗長為30的矩形窗，每次滑窗均在時頻圖上有明顯的痕跡，STFT方法受窗長影響較大。如圖1(b)所示，采用Wigner－Vile分布方法，三角調頻的時頻分布中間出現(xiàn)了“虛假”，有用的時變譜圖變得模糊，因此STFT和Wigner－Vile分布均不適合分析三角調頻信號。

圖1 不同方法下三角調頻脈沖信號的時頻分布

與STFT方法仿真相同信號，信噪比為20 dB和7 dB時，相位差分法擬合曲線和理想曲線圖對比如圖2和圖3所示。

圖2 SNR=20 dB時三角調頻相位差分法曲線

圖3 SNR=7 dB時三角調頻相位差分法曲線

由仿真結果可看出，直接根據瞬時頻率的定義展開的相位差分法，具有算法簡單、運算量小等特點。不存在“交叉項”干擾，尤其適合工程實踐上的大數(shù)據量測量，但抗噪性能并不理想，當信噪比較小時，測量性能較差，因此更適用于高信噪比的信號。

將算法用Matlab仿真200次，測量所得各脈內調制參數(shù)如表2所示。

表2 三角調頻脈沖信號脈內調制參數(shù)表

其中，B表示帶寬;u1、u2分別表示左、右調頻系數(shù);freqMax、freqMin分別表示脈內最大頻率和最小頻率。通過測量數(shù)據可以看出，調頻方向為上三角調頻擬合后參數(shù)值與理想值相比，測量誤差＜1%，測量效果良好。實際仿真時，Matlab仿真200次，用相位差分法，與使用STFT、Wigner－Vile分布相比，得到時頻曲線的時間大幅減小。因此，相位差分法具有更小的運算量。

5 結束語

針對雷達中常用的三角調頻脈沖信號、聯(lián)合時域、調制域兩方面測量其參數(shù)，調制域測量基于時域測量結果。在脈內調制參數(shù)測量上，舍棄受選擇窗長和窗函數(shù)影響較大的STFT方法與時頻分布有交叉項的Wigner－Vile方法，而采用相位差分法。仿真表明，其對于高信噪比信號，具有算法簡單、運算量少、測試誤差小的優(yōu)點，利于工程實踐的測試測量使用。

［1］劉東霞，趙國慶.脈內調制信號參數(shù)測量與分析［J］.現(xiàn)代雷達，2003，25(11):17 －20.

［2］郭允晟，蘇秉煒，方偉橋.脈沖參數(shù)與時域測量技術［M］.北京:中國計量出版社，1987.

［3］劉明亮，陸福敏，朱江淼.現(xiàn)代脈沖測量［M］.北京:科學出版社，2010.

［4］Czerwinski R N.Adaptive short－ time Fourier analysis［J］.IEEE Signal Processing Letters，1997，4(2):42 －45.

［5］羅利春.用維格納－威利時頻分布進行信號時頻分析［J］.航天電子對抗，2003(4):18 －21.

［6］韓中合，朱霄珣，李文華，等.基于EMD消除Wigner－Vile分布交叉項的研究［J］.汽輪機技術，2010，52(3):211 －214.