寬帶雙頻雙極化星載賦形反射面天線優(yōu)化設(shè)計

張新剛,吳 剛,鐘 鷹

(中國空間技術(shù)研究院西安分院,西安 710100)

1 引 言

在新一代的衛(wèi)星通信系統(tǒng)中不僅要求天線波束與服務(wù)區(qū)的形狀匹配,還需滿足大容量通信、收發(fā)一體化的要求,以處理同步進(jìn)行接收與發(fā)射的兩個分離頻段的信號,這要求天線能夠同時實現(xiàn)寬帶、雙頻、雙極化工作。賦形反射面天線具有結(jié)構(gòu)簡單、成本低、損耗小、易實現(xiàn)雙頻和雙極化等諸多優(yōu)點,在衛(wèi)星通信系統(tǒng)中應(yīng)用十分廣泛[1]。但要求賦形反射面天線在按照服務(wù)區(qū)的形狀對波束進(jìn)行賦形的同時實現(xiàn)寬帶、雙頻、雙極化特性,具有一定的難度,因為此時對反射面天線的設(shè)計是一個多目標(biāo)優(yōu)化問題,優(yōu)化過程中需要在多個優(yōu)化目標(biāo)之間進(jìn)行折衷。目前,國內(nèi)外現(xiàn)有星載賦形反射面天線的帶寬一般為500MHz或800MHz,限制了它在衛(wèi)星通信領(lǐng)域的應(yīng)用潛力。

目前,對賦形反射面天線進(jìn)行優(yōu)化設(shè)計的算法主要有最速下降法[2]、遺傳算法[3]、Minimax算法等[4],但在實際工程應(yīng)用中要求服務(wù)區(qū)內(nèi)任意一點的天線增益都要滿足設(shè)計指標(biāo),因此Minimax算法的設(shè)計理念最符合實際工程需要。以往在用Minimax算法進(jìn)行優(yōu)化時通常是通過各個評估站點的電場來構(gòu)建目標(biāo)函數(shù),然而設(shè)計要求一般是以增益的形式給出的。直接以各個評估站點的電場來構(gòu)建目標(biāo)函數(shù),必然導(dǎo)致同樣的增益偏差在增益要求較高的站點引起較大的電場偏差,優(yōu)化的重點過多地放在了增益要求較高的站點。為了消除各個評估站點增益對優(yōu)化過程的影響,盡可能拓寬天線帶寬,本文采用歸一化的電場值來構(gòu)建目標(biāo)函數(shù)。最后以一個服務(wù)區(qū)為美國國土、工作在Ku頻段寬帶雙頻雙極化的單偏置拋物反射面天線為例,在一定增益要求下對其寬帶性能進(jìn)行了優(yōu)化設(shè)計,從而充分說明了本文方法的有效性。

2 天線分析方法

由于反射面天線尺寸相對于工作波長(λ)屬于電大尺寸,一般采用PO法進(jìn)行分析計算[5-6]。反射面天線的遠(yuǎn)場表達(dá)式為

式中,Hi是饋源在反射面上的入射磁場,Je是反射面表面的感生電流,n是反射面上法線的單位矢量,μ0是自由空間的磁導(dǎo)率,ω是角頻率,k是波數(shù),k=2π/λ,r是觀察點所在的位置矢量,r= r ,r0=r/r,r′是場源所在的位置矢量,S′是反射面表面。

如果要考慮交叉極化增益或者遠(yuǎn)區(qū)旁瓣,PO法的計算精度就不能滿足要求,必須把反射面邊緣的繞射作用考慮進(jìn)去。本文采用物理繞射理論[7-8]對PO的計算結(jié)果進(jìn)行修正,以提高計算結(jié)果的精度。

式中,IT和MT分別是反射面邊緣等效電流和磁流,s是散射波單位矢量,t是反射面邊緣切向單位矢量,R= r-r′ ,Z0是自由空間特征阻抗,C是反射面邊緣的閉合曲線。

3 目標(biāo)函數(shù)及優(yōu)化算法

投影口徑為圓形的單偏置反射面天線的結(jié)構(gòu)示意圖如圖1所示,其中 a是投影面口徑的半徑,H是偏置距離,f是焦距,θ0是饋源的偏置角。

圖1 單偏置賦形反射面天線的幾何結(jié)構(gòu)Fig.1 Geometry of an offset single-reflector antenna

在如圖1所示的坐標(biāo)系中,反射面可以通過下式表示:

式中,t和φ為在口徑投影坐標(biāo)系下的坐標(biāo)變量,Cnm和Dnm是反射面的展開系數(shù),Fnm(t)是修正的Jacobi多項式。

在對反射面天線進(jìn)行優(yōu)化設(shè)計前,a、H、f、θ0以及饋源的具體形式作為基本設(shè)計參數(shù)已經(jīng)確定。優(yōu)化所要進(jìn)行的工作就是根據(jù)具體的設(shè)計要求,對由Cnm和Dnm構(gòu)成的優(yōu)化變量x進(jìn)行調(diào)節(jié),從而產(chǎn)生一個滿足設(shè)計要求的賦形波束。其中:

為了保證優(yōu)化后反射面的形狀與標(biāo)準(zhǔn)的單偏置相比不發(fā)生太大的形狀變化,C00保持不變。

天線服務(wù)區(qū)域是一個連續(xù)的區(qū)域,優(yōu)化時不可能對服務(wù)區(qū)域內(nèi)所有位置的性能都進(jìn)行評估,必須把這個連續(xù)的區(qū)域離散化,具體的做法是按照一定的規(guī)則在服務(wù)區(qū)域內(nèi)放置足夠多的觀測站點[9-10]。

對天線的帶寬性能有要求時,優(yōu)化目標(biāo)函數(shù)中要考慮頻率的影響。具體優(yōu)化計算時,在頻帶內(nèi)選擇一定數(shù)目的頻點進(jìn)行優(yōu)化,當(dāng)這些頻點的性能滿足要求時,就可以認(rèn)為整個頻帶內(nèi)的性能都滿足設(shè)計要求。當(dāng)以主極化增益和帶寬作為主要評判標(biāo)準(zhǔn),通常采用下式所示的方式來構(gòu)建目標(biāo)函數(shù):

式中,y為在服務(wù)區(qū)內(nèi)放置的J個觀測站點所形成的觀測變量,y=(y1,y2,…,yJ)T;K是優(yōu)化的頻段內(nèi)所取頻點的數(shù)目;wjk和Dk(yj)分別是第k個頻點時第j個觀測站點的權(quán)重系數(shù)和設(shè)計要求的主極化電場值;fk(x,yj)是第 k個頻點設(shè)計參量為x的情況下第j個觀測站點的主極化電場的計算結(jié)果。

Minimax方法通過改變參數(shù)向量x來使得目標(biāo)函數(shù)F(x,y)的最大值最小化,可以寫成如下形式:

目標(biāo)函數(shù)F(x,y)確定之后,就可以通過合適的算法來尋找其最小點從而使得設(shè)計結(jié)果滿足設(shè)計指標(biāo)要求。但是由于天線的方向性以及實際需求,在工作頻帶內(nèi)對各個站點的增益要求是不同的。假定對編號為1和2的兩個站點的設(shè)計增益分別為G1、G2,優(yōu)化過程中兩個站點的計算增益與設(shè)計增益的偏差分別為ΔG1、ΔG2,計算所得的電場值分別為 f1(x)、f2(x),則此時如果采用式(7)所示的方法來建立目標(biāo)函數(shù),此時兩個站點偏差函數(shù)的數(shù)值比為

從上式可以看出,偏差函數(shù)的數(shù)值比不僅與權(quán)重系數(shù)、增益偏差有關(guān),還與各個站點的設(shè)計增益有關(guān);當(dāng)權(quán)重系數(shù)和增益偏差相同時,設(shè)計增益較高的站點偏差函數(shù)的值較大。

Minimax優(yōu)化算法的核心思想是使所有站點偏差函數(shù)的最大值最小化,因此優(yōu)化時把優(yōu)化重點放在偏差函數(shù)值較大的那些站點。采用式(7)所示的方法來建立目標(biāo)函數(shù),站點的設(shè)計增益對其偏差函數(shù)的數(shù)值存在正比的影響,優(yōu)化重點將放在增益較高的站點,必然會影響到優(yōu)化進(jìn)行的方向。Minimax算法是一種局部搜索迭代算法,優(yōu)化過程中不同的進(jìn)行方向可能導(dǎo)致不同的結(jié)果。為了使優(yōu)化結(jié)果盡可能地靠近全局最優(yōu)解,本文采用歸一化的電場來構(gòu)建目標(biāo)函數(shù),具體形式如下式所示:

如果采用式(11)所示的方法來建立目標(biāo)函數(shù),此時1和2兩個站點的偏差函數(shù)值之比為

從上式可以看出,偏差函數(shù)的值僅與權(quán)重系數(shù)、增益偏差有關(guān);在權(quán)重系數(shù)相同的條件下,同樣的目標(biāo)函數(shù)值在不同的站點引起的增益偏差是相同的,消除了各站點不同的設(shè)計增益對目標(biāo)函數(shù)的影響,有利于獲得更好的優(yōu)化結(jié)果。

4 仿真實例及結(jié)果

下面通過對一個工作在Ku頻段寬帶的單偏置拋物反射面進(jìn)行優(yōu)化設(shè)計來說明算法的有效性。天線的服務(wù)區(qū)域為美國國土,具體形狀如圖2所示。圖中小黑點是優(yōu)化過程中部分觀測站點位置,陰影區(qū)域是降雨量較大的區(qū)域,圖中的坐標(biāo)有±0.1°的誤差[11]。

圖2 服務(wù)區(qū)形狀示意圖Fig.2 Sketch map of the coverage area

具體設(shè)計要求是:天線上、下行的帶寬都要達(dá)到1000MHz;工作頻段為 :上行 13.75～ 14.75GHz,下行11.45～12.45GHz;上行為 y極化,下行為 x極化。在整個頻段內(nèi),服務(wù)區(qū)的上行增益為30dBi(允許服務(wù)區(qū)邊界區(qū)域降低0.5dB);下行增益為29dBi(允許服務(wù)區(qū)邊界區(qū)域降低0.5dB);為補償雨衰引起的信號衰減,陰影區(qū)的下行增益需要加強(qiáng)3dB,達(dá)到32dBi。

進(jìn)行優(yōu)化設(shè)計時,天線選用單偏置拋物反射面天線,具體結(jié)構(gòu)形式如圖1所示,其中a為0.762 m,f為1.506 m,H 為1.245 m,θ0為 42.77°。饋源為雙極化的(cosθ)q式饋源,q=14.28。優(yōu)化時 N=9,M=5,優(yōu)化變量的個數(shù)為89(不含C00)。

當(dāng)設(shè)計要求和基本參數(shù)確定后,通過自編的Matlab程序來進(jìn)行優(yōu)化設(shè)計。當(dāng)采用式(7)的方式來構(gòu)建目標(biāo)函數(shù)時,優(yōu)化后服務(wù)區(qū)邊緣的電場值不能在要求頻帶內(nèi)達(dá)到設(shè)計指標(biāo)要求。采用式(11)的方式來構(gòu)建目標(biāo)函數(shù)時,優(yōu)化后得到的遠(yuǎn)場等值線圖如圖3和圖4所示。其中,圖3為在下行頻段的遠(yuǎn)場方向圖,圖4為上行頻段內(nèi)的遠(yuǎn)場方向圖。圖中的“×”點與圖2中黑點是對應(yīng)的,表示觀察站點的位置,外圍的“×”點代表國土地邊界。

從圖中可以看出,在11.45～12.45GHz下行頻帶內(nèi),服務(wù)區(qū)內(nèi)任意一點的增益都在29dBi以上,并且在95%以上的服務(wù)區(qū)域內(nèi)都高于30dBi,陰影區(qū)域的增益都在32dBi以上;在13.75～14.75GHz上行頻帶內(nèi),服務(wù)區(qū)內(nèi)任意一點的增益都在30dBi以上;天線的上、下行帶寬都為1000MHz,具有良好的寬帶性能。

圖3 優(yōu)化后得到的下行頻段遠(yuǎn)場方向圖Fig.3 Contoured patterns in downlink band of the optimized antenna

圖4 優(yōu)化后得到的上行頻段遠(yuǎn)場方向圖Fig.4 Contoured patterns in uplink band of the optimized antenna

優(yōu)化后反射面相對于標(biāo)準(zhǔn)單偏置拋物面形變量的等值線圖如圖5所示。從圖中可以看出優(yōu)化后反射面形變量在(-0.7λ,0.6λ)范圍內(nèi),且沒有發(fā)生突變的區(qū)域,整個表面比較光滑,滿足實際加工要求。

圖5 優(yōu)化后反射面形變量等值線圖(單位:λ)Fig.5 Contoured patterns of the deviation in z direction

5 結(jié) 論

本文在用Minimax算法對賦形反射面天線帶寬性能進(jìn)行優(yōu)化設(shè)計過程中,采用歸一化的電場值來構(gòu)建目標(biāo)函數(shù),消除了設(shè)計增益對優(yōu)化過程的影響,有效地拓寬了天線的帶寬。從對實際工程設(shè)計實例的仿真結(jié)果可以看出,在滿足服務(wù)區(qū)特定增益要求的前提下,賦形反射面天線的上、下行頻帶的帶寬都達(dá)到了1000MHz,比實際應(yīng)用中的同頻段天線帶寬至少提高了200MHz,具有良好的寬帶性能;并且優(yōu)化得到的反射面表面比較光滑,滿足實際加工要求。優(yōu)化結(jié)果對工程應(yīng)用具有重要的參考價值,為將來更寬的Ka頻段天線優(yōu)化設(shè)計奠定了基礎(chǔ)。

[1] Roederer A G,Jensen N E,Crone G A E.Some European satellite-antenna developments and trends[J].IEEEAntennas and Propagation Magazine,1996,38(2):9-21.

[2] Chou H,Chou H T.Fast SDM for shaped reflector antenna synthesis via patch decompositions in PO integrals[J].Progress in Electromagnetic Research,2009,92:361-375.

[3] Acila S L,Carpes W P,Vasconcelos JA.Optimization of an offset Reflector Antenna Using Genetic Algorithms[J].IEEE Transactions on Magnetics,2004,40(2):1256-1259.

[4] Klein C.Design of shaped-beam antennas through minimax gain optimization[J].IEEE Transactions on Antennas and Propagation,1984,32(9):963-968.

[5] 李鵬,張志華,王為,等.基于不同網(wǎng)格形式的反射面天線方向圖分析[J].系統(tǒng)工程與電子技術(shù),2009,31(2):347-351.LI Peng,ZHANG Zhi-hua,WANG Wei,et al.Far field pattern numerical analysis of reflector antenna based on different grid forms[J].Systems Engineering and Electronics,2009,31(2):347-351.(in Chinese)

[6] 王偉,段寶巖,馬伯淵.重力作用下天線反射面形變及其調(diào)整角度的確定[J].電波科學(xué)學(xué)報,2008,23(4):645-650.WANG Wei,DUAN Bao-yan,MA Bo-yuan.Gravity deformation and best rigging angle for surface adjustment of large reflector antennas[J].Chinese Journal of Radio Science,2008,23(4):645-650.(in Chinese)

[7] Johansen P M.Uniform physical theory of diffraction equivalent edge currents for truncated wedge strips[J],IEEE Transactions on Antennas Propagation,1996,44(7):989-995.

[8] 肖疆,徐曉文,董濤.多孔毀傷時賦形反射面天線輻射特性的研究[J].電波科學(xué)學(xué)報,2003,18(6):648-651.XIAO Jiang,XU Xiao-wen,DONG Tao.Study on the radiation characteristics of a shaped beam reflector antenna under the multi-hole damnification[J].Chinese Journal of Radio Science,2003,18(6):648-651.(in Chinese)

[9] Merad,Bendimerad F,Meriah S.Design of linear antenna arrays for side lobe reduction using the tabu search method[J].The International Arab Journal of Information Technology,2008,5(3):219-222.

[10] Phjlara D A.Improving the beam efficiency of an offset parabolic reflector antenna for space-borne radiometric application[J].Progress In Electromagnetic Research C,2009,10:143-150.

[11] Duan D W.A generalized diffraction synthesis technique for high performance reflector antenna[J].IEEE Transactions on Antennas and Propagation,1995,43(1):27-39.