苗雪平，張林讓，雷 宇

(西安電子科技大學(xué)雷達(dá)信號(hào)處理國(guó)家重點(diǎn)實(shí)驗(yàn)室，陜西西安 710071)

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

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

1 三角調(diào)頻脈沖信號(hào)模型

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

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

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

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

脈沖時(shí)域參數(shù)的測(cè)量以脈沖包絡(luò)為基礎(chǔ)，設(shè)脈沖實(shí)、虛部信號(hào)分別為I(t)和Q(t)，計(jì)算信號(hào)包絡(luò)的表達(dá)式為

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

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

2.1 脈沖個(gè)數(shù)測(cè)量

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

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

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

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

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

(4)重復(fù)步驟(2)至步驟(3)直到標(biāo)志位flag=0，則結(jié)束。

2.2 脈沖寬度測(cè)量

在測(cè)量脈沖個(gè)數(shù)過程中，記錄所有近似上升沿、近似下降沿判決時(shí)刻，對(duì)信號(hào)的某個(gè)脈沖，近似上升沿判決時(shí)刻為t'升，近似下降沿判決時(shí)刻為t'降，則計(jì)算脈沖寬度表達(dá)式為

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

在每個(gè)脈沖的脈沖寬度內(nèi)，運(yùn)用密度分布平均法，計(jì)算每個(gè)脈沖精確的頂值top和底值base。并求得每個(gè)脈沖精確幅度amp為二者之差。

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

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

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

3 脈內(nèi)調(diào)制參數(shù)測(cè)量

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

3.1 短時(shí)傅里葉變換

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

3.2 Wigner－Vile分布

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

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

3.3 相位差分法

假定單一分量連續(xù)時(shí)間信號(hào)s(t)解析表達(dá)式為

對(duì)于采集的離散信號(hào)可通過信號(hào)的相位差分來估計(jì)信號(hào)的瞬時(shí)頻率。采用中心有限差分估計(jì)信號(hào)的瞬時(shí)頻率fc(n)表達(dá)式為

其中，fs為信號(hào)的采樣率;φ(n)為第n點(diǎn)的瞬時(shí)相位。

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

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

4 仿真結(jié)果

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

表1 三角調(diào)頻脈沖信號(hào)時(shí)域參數(shù)表

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

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

圖1 不同方法下三角調(diào)頻脈沖信號(hào)的時(shí)頻分布

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

圖2 SNR=20 dB時(shí)三角調(diào)頻相位差分法曲線

圖3 SNR=7 dB時(shí)三角調(diào)頻相位差分法曲線

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

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

表2 三角調(diào)頻脈沖信號(hào)脈內(nèi)調(diào)制參數(shù)表

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

5 結(jié)束語

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

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

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

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

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

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

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