機電掃描天線陣列接收信號調(diào)制效應(yīng)分析

謝洪森,周 鵬,遲志艷,朱春平

(1.海軍航空工程學(xué)院青島校區(qū),山東青島266041;2.海軍航空兵學(xué)院,遼寧葫蘆島125001;3.中國人民解放軍91395部隊,北京102443)

0 引言

利用天線方向圖的方向性在空域進行周期性的掃描,是雷達或目標獲取角度參數(shù)的主要方式。就天線而言,要在空域形成動態(tài)的波束掃描,通常有兩種方式:一是通過伺服系統(tǒng)和傳動機構(gòu)使天線反射面周期性機械轉(zhuǎn)動,從而在空域形成動態(tài)的掃描波束,這種方式稱之為機械掃描,其主要特點是系統(tǒng)設(shè)計簡單,技術(shù)成熟,但機械運動慣性大,掃描速度低,如環(huán)視雷達、跟蹤雷達等;二是電掃描,采用這種方式時天線反射體和饋源不進行機械運動,而是通過周期性地改變各個陣元的相位、頻率或幅度等參量,實現(xiàn)天線的波束在空域中的周期掃描[1]。其主要特點是掃描速度高,波束控制靈活精確,但系統(tǒng)設(shè)計復(fù)雜,設(shè)備相對昂貴。

作為雷達系統(tǒng)的一個特殊種類——塔康(Tactical Air Navigation,TACAN)系統(tǒng),即戰(zhàn)術(shù)空中導(dǎo)航系統(tǒng),它屬于相位時間復(fù)合的航空近程無線電導(dǎo)航系統(tǒng)。系統(tǒng)利用天線方向圖的周期性旋轉(zhuǎn)和發(fā)射周期性基準脈沖,為飛機提供角度信息,利用發(fā)射詢問脈沖為飛機提供距離信息。該系統(tǒng)可安裝在機場、航路點、運動載體上,可為以導(dǎo)航臺為中心,半徑幾百公里范圍內(nèi)的飛機或運載體提供導(dǎo)航服務(wù),機動性和靈活性較高,特別適合軍事應(yīng)用[2]。

本文針對TACAN系統(tǒng)的一種特殊應(yīng)用場合,其天線既存在機械掃描又實施電掃描方式情況下,首先對TACAN系統(tǒng)測角原理進行了簡要分析,建立了TACAN系統(tǒng)接收信號模型,并對其頻譜特性進行了分析,計算機仿真結(jié)果與理論分析相吻合,最后對TACAN系統(tǒng)性能進行了分析,給出了論文的結(jié)論。分析研究系統(tǒng)天線波束空域掃描的調(diào)制效應(yīng)具有重要的現(xiàn)實意義。

1 TACAN測角原理及其接收信號模型

TACAN電掃天線采用若干個陣元均勻排列構(gòu)成一個圓形陣列結(jié)構(gòu),通過控制各個陣元的饋電相位形成空域的周期性掃描波束。為了進一步改善俯仰方向上的測角性能,TACAN天線是由若干個層疊的圓周陣列構(gòu)成。而對于存在水平方向的機械運動,層疊陣列的分析方法與單層陳列的分析結(jié)果是一致的,因此本文假定TACAN天線簡化為單層圓周陣列構(gòu)成。

目標要獲得相對于TACAN的距離坐標,是通過向TACAN發(fā)射經(jīng)過編碼的詢問脈沖,TACAN經(jīng)過適當延遲向目標發(fā)射應(yīng)答脈沖,目標通過分析應(yīng)答脈沖的延遲判斷相對TACAN的距離。而目標要獲得相對于TACAN的角度坐標,則是通過測量TACAN發(fā)射的周期性脈沖和TACAN周期性旋轉(zhuǎn)的波束進行位置的判定[2]。

TACAN天線水平方向圖如圖1(a)所示,其數(shù)學(xué)表達式可寫為

式中,A和B與陣元特性、天線結(jié)構(gòu)以及饋電信號有關(guān),K為一確定整數(shù)。該方向圖在θ=0時存在一個最大值,在處存在極大值,一共有K個極大值。對TACAN而言,參數(shù)K=9。更為形象的解釋是TACAN天線方向圖是一個心臟形狀(圖1(a)虛線所示)上疊加9個瓣而已。如果該天線方向圖進行水平周期性旋轉(zhuǎn),那么接收方將接收到一個類似AM調(diào)制的信號,假定天線方向圖的旋轉(zhuǎn)頻率為那么接收信號中將存在2個非零點頻,分別為Ω和KΩ。如果天線方向圖在進行水平旋轉(zhuǎn)的同時,還以頻率為Ω發(fā)射定向基準脈沖,該脈沖的起始相位與天線方向圖指向某個特定方向同步。這樣接收機將接收信號的最大值與基準脈沖的相位差就可以判斷出目標的相對位置,為了提高測角精度,還可以利用KΩ頻率信號,得到更精確的測角精度[2]。對TACAN,Ω=30π。

要實現(xiàn)天線方向圖的旋轉(zhuǎn),現(xiàn)以電掃TACAN為主要研究對象,其幾何關(guān)系如圖1(b)所示,接收機位于P點,某個陣元位于A點,Q為P點到XOY面的垂點,B和C分別是A點到線段OQ和OP的垂點。電掃天線是由Ne個陣元均勻排列構(gòu)成的圓形陣列,要實現(xiàn)天線方向圖以Ω速率旋轉(zhuǎn),則第n個陣元的饋電信號為

式中,A1和A2為調(diào)幅度,各陣元饋電信號初始相位TACAN接收機接收各個陣元發(fā)射的信號,其幾何關(guān)系如圖1(b)所示。接收信號經(jīng)過包絡(luò)檢波后為

2 電掃描與機械掃描的調(diào)制效應(yīng)分析

現(xiàn)在考慮一種特殊情況,就是圖1(b)所示的圓形陣列以速率Φ勻速水平旋轉(zhuǎn),假定旋轉(zhuǎn)方向滿足笛卡兒坐標右手系方向。那么簡單分析可知,接收信號將不僅受到饋電信號的影響,還可能受到陣列結(jié)構(gòu)和旋轉(zhuǎn)速率等多種因素的影響,而這些因素均有一個共同特征,即周期性,從而導(dǎo)致接收信號仍然是周期信號,其頻譜由若干線譜構(gòu)成。而TACAN測角主要依賴接收信號中的兩個線譜分量的時域測量,為此必須充分研究接收信號的線譜擴展情況,才能對TACAN測角性能進行分析。

圖1 TACAN測角原理

先討論陣列結(jié)構(gòu)對接收信號的影響問題。如圖2所示,假定觀測方向箭頭所指方向,那么觀測到的陣列結(jié)構(gòu)存在兩種狀態(tài),狀態(tài)1如圖2(a)所示,此狀態(tài)表示某個陣元位于觀測方向上,接收機可以同時接收到N+1個陣元發(fā)射的信號,考慮到對稱性,N取偶數(shù)。顯然該狀態(tài)表現(xiàn)為離散形式,在陣列旋轉(zhuǎn)一周的過程中,均勻出現(xiàn)Ne次。因此該狀態(tài)的出現(xiàn)頻率為旋轉(zhuǎn)速率Φ的Ne倍。兩個相鄰的狀態(tài)1之間的過程稱之為狀態(tài)2,該狀態(tài)是一個連續(xù)的過程,在這個過程內(nèi)接收機可以同時接收到N個陣元發(fā)射的信號。在陣列旋轉(zhuǎn)一周的過程中,均勻出現(xiàn)Ne次,因此該狀態(tài)出現(xiàn)的頻率與狀態(tài)1一樣。有了這種概念,并結(jié)合式(3)可以初步分析陣列結(jié)構(gòu)的周期旋轉(zhuǎn)對信號周期性的影響[5]。假定起始時刻為圖2(a)所示,假定1號陣元的饋電相位為0,經(jīng)過時間T之后,m號陣元位于圖2(a)所示的1號陣元位置,其饋電相位也為0,即

考慮到陣列的周期性旋轉(zhuǎn),時間T還滿足:

綜合式(4)和式(5),可以求得公共周期:

這說明接收信號的基準頻率是饋電信號頻率和陣列旋轉(zhuǎn)頻率的和。而對于TACAN系統(tǒng)而言,接收機接收的信號頻率發(fā)生了改變將導(dǎo)致嚴重的性能惡化,解決的辦法是改變饋電信號頻率。僅僅分析其基本周期是不夠的,這是因為接收信號是若干個基準周期整數(shù)倍的正弦信號疊加,接收信號的頻率發(fā)生了周期延拓,頻率分量可能發(fā)生交疊。為此必須進一步分析接收信號的數(shù)學(xué)模型。

圖2 觀測方向上的陣列結(jié)構(gòu)

為了簡化分析,假定陣列饋電信號為相位調(diào)制的正弦信號,目標的方位角度為0。狀態(tài)1是離散狀態(tài),因此在時域上表現(xiàn)為多個狄拉克函數(shù)的求和形式,陣列結(jié)構(gòu)在狀態(tài)1中保持不變,即求和公式中n的取值范圍保持不變,且指數(shù)部分也保持不變,而饋電信號的相位卻發(fā)生周期性變化[6],因此容易得到狀態(tài)1時刻的接收信號可以表示為

式中,mod表示取余算子。為了適合計算機仿真需要,需進行接收信號的數(shù)字化采樣,而對于線譜信號,其采樣率必須選擇合適,否則容易出現(xiàn)頻譜泄露問題。為此,假定狀態(tài)2采樣M-1個點,狀態(tài)1對應(yīng)一個采樣點,那么采樣率[7]可表示為

故狀態(tài)1的接收信號數(shù)學(xué)模型可表示為

其中饋電信號相位中的取余算子并不需要[8],因此可進一步表示為

再分析狀態(tài)2時的接收信號模型,如圖2所示的狀態(tài)2中,陣元個數(shù)始終為N個,因此求和范圍與狀態(tài)1不同,其次狀態(tài)2為一個連續(xù)狀態(tài),因此其饋電信號的相位也是連續(xù)變化的,饋電相位的初始相位受到由于陣列旋轉(zhuǎn)引入的陣列編號的變化。而指數(shù)部分在每個狀態(tài)2中均是相同的:

利用式(8)可得到數(shù)字信號模型:

從式(12)可以看出,信號的基準周期與式(6)分析結(jié)論一致,其次考慮到旋轉(zhuǎn)效應(yīng)和陣列結(jié)構(gòu)的周期性,狀態(tài)2的周期性表現(xiàn)為頻率發(fā)生了周期延拓,頻率周期為

接收信號的數(shù)學(xué)模型可用式(10)和(12)進行完整表述。利用該模型可以得出這樣的結(jié)論:當饋電信號頻率為Ω,陣列旋轉(zhuǎn)速率為Φ,那么接收信號的線譜位置為

由此可見,接收信號的線譜結(jié)構(gòu)發(fā)生了兩個方面的畸變,一是基準頻率不再是饋電信號的頻率,而是加上了旋轉(zhuǎn)速率,二是由于陣列的周期結(jié)構(gòu),導(dǎo)致了頻率發(fā)生了周期延拓[8]。

除了考慮譜線出現(xiàn)的位置,還必須對譜線的強度進行分析。信號的周期性結(jié)構(gòu)可以表示為

該信號是某個非周期信號r(t)的矩形窗截斷后的周期擴展,接收信號的頻譜就是r(t)的頻譜與sinc函數(shù)的卷積后再進行頻域抽樣后的線譜構(gòu)成,獲得一個譜線幅度的閉合表達式比較困難,不過r(t)的頻譜結(jié)構(gòu)是一個低通函數(shù),而sinc函數(shù)也是低通的,因此接收信號的線譜結(jié)構(gòu)表現(xiàn)為隨著階數(shù)的增加而逐漸衰減的過程。也就是說,譜線擴展的階數(shù)越高,其衰減越大。表1給出了3個點頻的各階擴展譜線的衰減程度。頻率范圍選擇為0～180 Hz。從表1的數(shù)據(jù)可以看出,135 Hz譜線的各階擴展分量衰減較小,其3階和4階擴展譜線出現(xiàn)在15 Hz附近,可能會對15 Hz信號的處理產(chǎn)生影響。

表1 各譜線衰減程度

3 調(diào)制效應(yīng)的計算機仿真分析

依據(jù)前面討論的信號模型,針對TACAN系統(tǒng)分析其接收機測角性能。由式(2)可知,發(fā)射信號中包含的調(diào)制信號頻率分別為0,Ω,9Ω,這三種信號滿足線性疊加原理,可以分別進行仿真和分析。假定天線陣列半徑為1 m,陣列由36個陣元構(gòu)成,觀測陣元數(shù)目為6～7個。在進行數(shù)字仿真過程中,假定狀態(tài)1與狀態(tài)2共采樣100個點,轉(zhuǎn)速為1 Hz,饋電信號頻率為14 Hz。

圖3給出了各個陣元等幅同相饋電時,接收信號的波形及其頻譜。接收信號的幅度表現(xiàn)為狀態(tài)1和狀態(tài)2交替,狀態(tài)2持續(xù)的時間為轉(zhuǎn)速的36倍。接收信號的頻譜由若干線譜構(gòu)成,分別為0,36,72,144等。

圖4給出了當饋電信號為14 Hz,各個陣元相位相差10°,接收信號的時域和頻域波形。從時域中可以看出,存在離散的脈沖信號,它是由于狀態(tài)1引入的。頻譜分析如式(13)分析結(jié)論一致。其中15 Hz信號是饋電信號頻率14 Hz與轉(zhuǎn)速1 Hz的和。21 Hz和51 Hz分別是36 Hz±15 Hz產(chǎn)生的。

圖5給出了饋電信號如式(2)所示時,接收信號的時域和頻域波形?？梢钥闯?時域和頻域波形均受到嚴重的破壞,時域波形不再是如圖1(a)所示的規(guī)則信號,從信號的頻域分析來看,其頻譜譜線增加了,這是由于饋電的3個頻率發(fā)生的周期延拓和交疊。一般而言,我們感興趣的是接收信號15 Hz和135 Hz附近的線譜結(jié)構(gòu),而這些線譜結(jié)構(gòu)顯然最終影響了TACAN系統(tǒng)的測角性能。從圖5可以看出,對于15 Hz的信號,其周圍存在9 Hz和21 Hz的譜線,這兩個譜線分別是36×4-135和36-15產(chǎn)生的,顯然9 Hz的譜線是135 Hz譜線4階擴展,因此其幅度低于有15 Hz的1階擴展譜線21 Hz分量。再觀察135 Hz附近的譜線,分別是129 Hz和144 Hz,它們分別是4×36-15和4×36產(chǎn)生的,顯然129 Hz譜線是15 Hz譜線的4階擴展,而144 Hz譜線是0 Hz的4階擴展,因此這兩個譜線的幅度相對較低。

圖3 等幅同相饋電時接收信號的時域頻域波形

圖4 等幅不同相位饋電時接收信號的時域頻域波形

圖5 TACAN接收信號的時域頻域波形

4 調(diào)制效應(yīng)對TACAN系統(tǒng)性能的影響分析

TACAN天線饋電信號包含3種線譜,即0頻率,15 Hz和135 Hz。3種頻率分量是線性疊加過程,因此針對每一個頻譜成分均可以利用前述的數(shù)學(xué)模型進行分析。而考慮到TACAN系統(tǒng)的測角原理是對15 Hz信號和135 Hz信號的處理,為此利用前面的信號模型,重點分析15 Hz和135 Hz附近的譜線分量。

從前面分析可知,3種頻率會產(chǎn)生周期擴展,而對于TACAN系統(tǒng)所關(guān)心的是15 Hz和135 Hz附近的譜線干擾。TACAN接收機采用的是15 Hz和135 Hz的窄帶模擬濾波器,其帶寬有一定的限制,當這兩個頻率分量附近存在其他譜線時可能會對其產(chǎn)生影響。

對于0頻率,其譜線擴展的位置出現(xiàn)在:

假定轉(zhuǎn)速為1 Hz,陣元個數(shù)為36,則在15 Hz最近的線譜分量為0階擴展0 Hz和1階擴展36 Hz,對135 Hz最近的譜線為3階擴展分量108 Hz和4階擴展分量144 Hz。對于15 Hz分量,其主要影響到135 Hz的譜線為3階擴展分量123 Hz,4階擴展分量129 Hz和159 Hz。135 Hz分量對15 Hz譜線的影響為3階擴展分量27 Hz,4階擴展分量9 Hz。

以上分析僅對某種一定情況而言,如果陣元個數(shù)發(fā)生改變,或者陣列旋轉(zhuǎn)速度發(fā)生改變,前面的分析結(jié)果要發(fā)生相應(yīng)的改變。一般而言,轉(zhuǎn)速越大,低階譜線擴展分量的影響就越大,例如若轉(zhuǎn)速提高到4倍,則135 Hz的1階擴展分量就會出現(xiàn)在9 Hz上,從而會增加對15 Hz信號的干擾。其次是當某些特定的轉(zhuǎn)速時,可能會發(fā)生某階擴展譜線分量出現(xiàn)的15 Hz或135 Hz分量上,如當其他條件不變,轉(zhuǎn)速為10/3 Hz,容易分析出,135 Hz的1階擴展分量出現(xiàn)在15 Hz上,由于相位的差異,會導(dǎo)致15 Hz信號分量發(fā)生畸變。

考慮到譜線的位置是已知的,可以設(shè)計合理的濾波器實現(xiàn)15 Hz譜線和135 Hz譜線的提取。由于接收信號的頻譜寬度受到接收機帶寬的限制,故15 Hz通道的輸出頻譜只含有15 Hz及其附近的譜線,135 Hz通道的輸出只含有135 Hz及其附近的譜線。設(shè)計一個含有對應(yīng)零點的FIR濾波器就可將這些位置已知的譜線加以濾除。濾波器的幅頻響應(yīng)如圖6所示,其設(shè)計方法是將15 Hz或135 Hz附近的已知譜線位置作為濾波器幅頻響應(yīng)的零點,從而能夠有效抑制干擾信號的影響。

圖6 濾波器幅頻響應(yīng)

5 結(jié)束語

天線波束的空域掃描方式包括電掃描和機械掃描。本文以TACAN系統(tǒng)為研究對象,在天線陣列實施電掃描同時存在陣列的機械旋轉(zhuǎn)運動時,給出了接收信號的數(shù)學(xué)模型。經(jīng)理論推導(dǎo)和計算機仿真結(jié)果表明,接收信號的頻率是饋電信號頻率與轉(zhuǎn)速頻率的和,其譜線發(fā)生了周期性擴展。針對這種情況,通過設(shè)計合理的數(shù)字濾波器可以抵消譜線擴展引入的接收信號失真問題。研究還表明,當轉(zhuǎn)速處于某些特定頻率時,擴展的譜線之間會發(fā)生交疊現(xiàn)象,可能會對TACAN系統(tǒng)測角產(chǎn)生不良影響,其影響的量化分析需進一步研究。其次,陣列的機械旋轉(zhuǎn)速率可能存在不穩(wěn)定的情況,表現(xiàn)在接收信號的譜線發(fā)生抖動或展寬,因此必須采用自適應(yīng)濾波器替代原先的FIR濾波器。

[1]吳德偉.航空無線電導(dǎo)航系統(tǒng)[M].北京:電子工業(yè)出版社,2010:61-67.

[2]謝洪森.艦載機著艦無線電引導(dǎo)關(guān)鍵技術(shù)研究[D].煙臺:海軍航空工程學(xué)院,2012:8-11.

[3]彭政諭.陣列天線波束賦形技術(shù)研究與應(yīng)用[D].杭州:浙江大學(xué),2014:35-38.

[4]戴幻堯,李永禎,陳志杰,等.電掃偶極子相控陣天線的空域極化特性分析[J].國防科技大學(xué)學(xué)報,2010,32(1):84-89.

[5]趙萬春,師君.圓形陣列天線方向圖的分布函數(shù)優(yōu)化方法[J].電訊技術(shù),2010,50(7):91-96.

[6]葛佩,李明,劉翔.一種低旁瓣子陣線自適應(yīng)波束形成方法[J].雷達科學(xué)與技術(shù),2012,10(5):519-523.GE Pei,LI Ming,LIU Xiang.An Adaptive Beamforming Method at Subarray Level for Planar Array with Low Sidelobes[J].Radar Science and Technology,2012,10(5):519-523.(in Chinese)

[7]CHEN K,YUN X,HE Z,et al.Synthesis of Sparse Planar Arrays Using Modified Real Genetic Algorithm[J].IEEE Trans on Antennas and Propagation,2007,55(4):1067-1073.

[8]WILSON M J,MCHUGH R.Sparse-Periodic Hybrid Array Beamformer[J].IET Radar,Sonar and Navigation,2007,1(2):116-123.