天线方向图的分布函数优化方法概要.docx
- 文档编号:7480380
- 上传时间:2023-01-24
- 格式:DOCX
- 页数:22
- 大小:45.17KB
天线方向图的分布函数优化方法概要.docx
《天线方向图的分布函数优化方法概要.docx》由会员分享,可在线阅读,更多相关《天线方向图的分布函数优化方法概要.docx(22页珍藏版)》请在冰豆网上搜索。
天线方向图的分布函数优化方法概要
第50卷第7期
2010年7月电讯技术TelecommunicationV01.50No.7EngineeringJul.2010文章编号:
1001—893X(2010)07—0091—06
圆形阵列天线方向图的分布函数优化方法睾
赵万春1,师君2
(1.中国人民解放军驻157厂军事代表室,四川彭州611930;2.电子科技大学电子工程学院,成都610054)
摘要:
以离散有限长复指数和的计算为基础,讨论了二雏阵列的方向图问题。
通过引入分布函数的概念,可知离散有限长复指数和为其分布函数离散傅里叶变换在1点的值。
基于此结论,分析了二维阵天线阵元分布于天线方向图之间的关系:
二维阵列任意方向的方向图为其在该方向分布函数的离散傅里叶变换,该分布函数可通过统计其投影在该方向的临近区域处阵元数目获得。
然后,通过对混合圆形阵列天线方向图的分析,验证了该关系的正确性,并说明了上述关系在分析阵列方向图的应用。
在此基础上,进一步分析了阵元权重系数对方向的影响,发现该方法与通过改变阵列天线阵元分布的方法等效。
最后,提出了基于分布函数的圆形阵列方向图优化方法,并通过仿真实验验证了该方法的可行性。
该方法可通过改变圆形阵列阵元沿径向的权重系数使得其对应的任意方向的方向图与期望的窗函数近似。
关键词:
圆形阵列;天线方向图;分布函数;权系数优化
中图分类号:
,IN95文献标识码:
Adoi:
10.3969/j.issn.100l一893x.2010.07.019
ResearchofDensityFunctionOptimizationMethod
forCircularArrayPatterns
ZHAOWan.chun1,SHIJun2
(1.PLAMilitaryDelegationOfficefor157Factory,Pengzhou611930,China;2.SchoolofElectronic
Engineering,UniversityofElectronicScienceandTechnologyofChina,Chengdu
Abstract:
Thispaperdiscussestheradiationpatternof2一Dcirculararmy610054,China)thecalculationofexpo—accordingto
nentialsumoffinitediscrete.timefunction.Byintroducingthedensityfunction.itisdemonstratedthatt}levalueofallexponentialsumoffinitediscrete-timefunctioniSequaltothatofthediscreteFouriertransformofitsden—
onsityfunctionatpoint1.Basedthisconclusion.allapproachtoanalysetheradiationpatternof2一DarrayUS-
ingdensityfunctioniSpresented.anditisshownthatthedirectionalradiationpatternofa2一DarrayistheFouriertransformofitsdensityfunctionint}lecorrespondingdirection.Thentheradiationpatternsofhybrid2一Dcirculararraysarediscussedtodemonstratethevalidityandapplicationofthis
onapproach.Then,theeffectofpositiveweightcoefficientthemdiationpatternisdiscussed。
anditisfoundthattheeffectofpositiveweight
coefficientisequivalenceto山edistributionofelements.Finally.aweighoptimizationmethodforcirculararrayispresented,whichensulL,StheoptimizedPSFsimilartotheFourier
approximatelybytransformofthedesiredwindowfunctionadjUStingtheweightcoefficients.
Keywords:
circulararray;antennapattern;densityfunction;weightcoefficientoptimization
天线方向图[1—3]描述的阵列天线的分辨率特征是阵
l引言列天线设计的重要指标之一。
数学上,天线方向图
的分析可看作带参数的振荡积分的求解问题。
但阵列理论广泛应用于雷达、声纳及通信等领域,
鼍收稿日期:
2010—01—26;修回日期:
2010—04—19
万方数据
Mnm.teleonline.cn
电讯技术
2010定
是,由于振荡积分相位函数具有多样性,可能为连续函数、周期函数或随机函数,该振荡积分的求解较为困难。
目前,方向图的优化主要采用基于统计模型的优化方法,如模拟退火算法(SA)t4-8J及遗传算法(GA)【9-13]等。
虽然上述方法可以对天线旁瓣进行优化设计,但无法揭示天线阵元布局与方向图之间的直接联系。
另外,对于高维阵列,当天线阵元数目较多时,上述方法运算量较大。
在实际工程中,阵列天线阵元的数目有限,振荡积分的求解问题可近似为离散有限长复指数和的计算。
通过引入分布函数的概念,本文介绍了离散有限长复指数和的计算方法,并将该方法应用于二维阵列天线方向图的分析,揭示了二维阵列天线阵元布局与其方向图的对应关系。
然后,利用该关系分析了满阵圆形阵列、稀疏圆形阵列及混合圆形阵列的布局与天线方向图的关系,验证了该分析方法的正确性。
在上述研究的基础上,本文进一步分析了加权对天线方向图的影响,并在此基础上提出了一种基于分布函数的圆形阵列方向图优化方法。
该方法通过对不同半径阵元权重的设计使得其对应的天线方向图达到与期望的窗函数相近的效果。
2离散有限长复指数和的计算
阵列天线方向图优化问题数学上可以看作是振荡积分求解问题,由于该问题分析较为复杂,目前主要采用驻定相位原理进行近似求解。
但是对于离散有限长信号,其指数和的计算可以采用本节给出的方法进行计算。
假设已知离散有限长信号.厂(n)=(0,1/4,2/4,0,1/4,1/4),需求解该信号的复指数和,其公式可表示为
∑exp(j・2rt・八,1)):
ej‘h・o+ei‘2Ⅱ‘1/4+ej。
2丌・扒+
ej‘孙‘0+ei’2n.1/4+ej‘2”1/4
(1)
根据加法的交换律和结合律,式(1)可改写为
∑exp(j・27r・八n)):
2・ej‘2Ⅱ‘o+3・ej。
h・1届+
1.ej’2"‘2/4+O・ej’2丌‘3/4(2)
式中,系数2、3、1和0表示各复指数项出现的频率(概率意义上)。
因此,通过引入密度(分布)函数D(i)={2,3,1,0}的概念,可得到:
・92・
万方数据
∑exp(j・2rr・八n)):
∑刀(i)ej咖州(3)
rt
i=0
明显地,公式(3)满足离散傅里叶变换的定义式。
定义D(i)的傅里叶变换为b(k),可得到:
∑exp(j・21r・/(n)):
D(1)
(4)
类似地,带有参数的离散有限长复指数和可表
示为
3
∑exp(j・27r・k・八,1))=∑D(i)eJ’知‘‘。
‘届:
D(后)
n
i=0
(5)
根据上面的例子可看出,带参数的离散有限长复指数和可通过对其分布函数进行傅里叶变换得到,该性质为本文提出的圆形阵列方向图优化方法的技术。
实际上,采用Lebesgue测度等方法,该性质可很容易地扩展到连续函数的振荡积分的计算。
由于离散有限长复指数和的性质已经满足阵列优化问题的需要,连续函数的振荡积分将在其它的研究中详细论述。
3阵列天线方向图
天线方向图是阵列天线的重要指标,反映了天线辐射电磁波在空间中的分布情况。
阵列天线设计的一个关键问题为通过对阵型、阵元布局及权重系数的优化以获得期望的天线方向图。
本节将利用上节介绍的计算方法推导二维阵列天线的方向图,并分析几种典型圆形阵列的方向图。
3.1信号模型
二维阵列天线的几何关系如图l所示。
图1阵列天线几何结构
rig.1
Geometryofarray
f1.1ltenna
不失一般性,可将二维阵列天线描述为一组二
维坐标构成的集合,记作:
而s4={;^:
j;;^=(;&,Yk)},其中,k表示阵元节点的序号,(未^,丸)表示阵元k
相对阵元中心的位置。
则而跗精确地描述了任意二
第7期赵万春等:
圆形阵列天线方向图的分布函数优化方法总第260期
维阵列天线的形状。
假设阵元中心相对目标中心的位置为
P。
。
恤=(xo,Yo,go)
(6)
式中,茗o、Yo和zo表示该位置在筇轴、Y轴和z轴的分量。
则二维阵列相对目标中心的位置可表示为
Os4={m}
(7)
p^全(戈I,Yk,gO)2
Pk+p∞n雠=
(髫o+筑,Yo+丸,Zo)
(8)
式中,Osa描述了二维阵列天线相对观测场景的位置。
假设目标点相对目标中心的相对位置为p缸,且雷达向目标区域发射单频电磁波,则第k个阵元接收到的回波信号为
s(p^,t)=exp(一j・27r・^・(t—Rtff/c))(9)
R铲=0五^一p洫|I
2(10)
式中,如为载波频率,R}为阵元到目标的距离,C表示光速,II・02表示向量2范数。
通过接收系统去载波处理,则阵列天线的接收回波可表示为
s(p^;Pm)=exp(j・2n・RV/a)pt∈ns4
(11)
式中,A为载波波长。
3.2方向图
利用上节建立的信号模型,二维阵列天线的方向图可表示为
z(;。
)全面I∑s(五;石)∥(五;;。
)(12)一‘k
9全∑ls(;&;石)J
2
(13)
%
式中,石表示参考点的位置,上标*表示复共轭算
子,缈表示信号s(;I;一O)的能量,P一。
表示图像空间
中的坐标。
由于z(;。
)为多元函数,为了便于分析,定义}
方向的方向图为
z}(口)全击∑s(jI;0・})_(五;石)(14)
一‘b
}=(玉,玉,o)
(15)口全r/Ro
(16)
式中,}表示图像空间中任意方向的单位矢量,如
图1所示;r表示像素;。
距离图像空间中心位置处
的距离;R。
表示二维阵列天线距离图像空间中心的距离;0为r和Ro的比,当r《Ro时,其近似为像素P。
对应的波束角。
万方数据
将式(11)代入式(14)中,并将△磁‘}用泰勒级
数展开,可得到}方向的方向图为
船(r)=面1・∑exp[j・2n・0・(五・阿州(17)
一
‘b
根据式(8),有:
jt・F=;一・F+血・}T
(18)
则式(17)可写为
z}(口)=召(口)・∑exp[J・27t・0・(j;;^・}T)以】
ns^
(19)
B(口)全面1・exp[j。
2n。
0・(;一。
尹)以】
(20)
将(;I・r)以看作式(5)中的相位函数,利用式(5)
可得到:
z}(口)=B(r)・∑ei∞8‘魄(i)=
B(口)・晚(一0)
(21)
由于IB(口)I=面l,则有:
zt(o)l=面1・lbt(一圳
(22)
式(22)的物理意义如图2所示。
从图中可以看出,对于任意二维阵列,可通过统计投影在i位置临
近区域处阵元数目获得该阵列在}方向的方向图,
而其该方向的方向图可通过对该分布函数进行离散傅里叶变换获得。
参一
令血命DFT上{,
0
0
图2二维阵列与方向图的关系
ris.2
Relatiomhipbetweentheradiationpattern
andthe8hapeofarray
对于圆形阵列,由于其阵列关于圆形对称(或在统计学意义上对称),其任意方向的一维方向图可完全描述该阵列的二维方向图。
3.3混合圆形阵列分析
上一节建立了任意阵列天线阵元布局和其对应的天线方向图的对应关系,该对应关系可用于分析各种平面阵列的方向图。
由于篇幅所限,本节以混合阵列为例,说明该方法的有效性。
混合阵列[14]是一种新型阵列天线结构,本文将考虑一种较为简单的混合圆形阵列,即阵列中部为满阵子阵列,周围为稀疏环形阵列,如图3和图4所示。
・93・
WWW.teleonline.cn
电讯技术
2010正
图3混合圆形阵列
Fig.3hybridcirculararray
图4混合圆形阵列分布函数
Fig.4Densityfunctionofhybridarray
图4为混合圆形阵列天线的方向分布函数。
可以看出,该分布函数可看作某窗函数被噪声调制后的结果,该窗函数可在一定程度上降低方向图的旁瓣,而噪声则在天线方向图中引入了随机性。
图5为混合圆形阵列天线的二维方向图与一维方向图,其峰值旁瓣比和积分旁瓣比分别为一21.77dB(存在随机性)和一15.45dB。
而圆形满阵元阵列的其峰值旁瓣比和积分旁瓣比分别为一15.87dB和
一13.94
dB。
因此,混合圆形阵列可在一定程度上
降低圆形阵列的旁瓣,达到节省天线阵元并优化阵列天线的方向图的目的。
n
l—Ps7
I…FT[DF】J
1
.
黜7㈣l{
酮捌
抵。
}l诣’
激图5混合圆形阵列某方向的方向图
rig.5
1一Dradiationpatternofhybridarray
・94・
万方数据
4基于分布函数的圆形阵列优化
根据上面的分析可以看出阵列分布函数与方向图存在对应关系,因此,可通过对分布函数的设计实现对方向图的优化。
进一步分析表明,通过阵元分布优化对方向图的优化效果与改变阵元权系数对方向图的优化效果等效。
本节将通过对阵列系数的设计实现基于分布函数的圆形阵列优化。
4.1权系数对方向图的影响
假设每个天线阵元上存在一个正权重系数w(pt),则阵列信号回波可写为
s(pk;ptar)=w(pk)・exp(j・27c・尺ⅣA)(23)
相应地,其方向图可表示为
z(j。
)=面1∑埘(五I)・s(;I;石)√(;。
;j。
)‘%
(24)
假设所有权系数均为自然数,则式(24)可表示为
z(五。
)=面1∑∑sG。
;石)∥(五;;。
)(25)‘%”
该公式的物理意义如图6所示,表明在某个阵元上调制权重系数埘相当于在该位置放置W个权重系数为1的阵元。
、、
//
、~J一一
图6权系数对方向图的影响
Fig.6Effectofweightcoefficient
on
radiation
pattern
当权重系数为正实数时,通过将式(25)乘以某个足够大的正实数,并对所有权重系数进行取整近似,可很容易地将上述结论扩展到正实数的情况,即改变权重系数即可改变阵列天线阵元的等效分布。
另外,与改变阵元数目相比,加权的方法可为任意正实数,更便于设计实现。
下面将给出一种基于分布函数的圆形阵列天线优化方法。
4.2基于分布函数的优化方法
由于圆形阵列关于圆点对称,如果其权系数也
第7期
赵万春等:
圆形阵列天线方向图的分布函数优化方法总第260期
关于圆点对称,则其方向图关于圆点对称,此时,只需对其一维方向图优化即可实现对该圆形阵列的优化。
因此,本文假设其权系数关于圆点对称,即其为半径的函数r,而与角度无关日。
假设该权重函数为加(r),则加权累积函数的积分(即加权分布函数的积分,其关系于概率密度函
数和概率累积函数的关系近似)%(Y)可表示为
2j事一1缈矗)
j们(r)d舢=
C谚(y)=
曩
2I
cos一(y/R)w(r)dr,Y≤0(26)
c∥(∞)一C0(一Y),Y>0
(27)
显然,%(Y)对应的加权分布函数为偶函数。
为了实现对圆形阵列的优化,可通过对权重系数加(r)的设计,使得其对应的加权分布函数满足特定的窗函数,即:
劬(y)=I
f
ccJ(u)du(28)
式中,∞(M)为期望的窗函数,如汉明窗、Chebyshev窗等。
由于%(Y)为偶函数,只需考虑Y≤0的情
况,则式(28)可表示为
2
—
COS
,Ly/
R、,加,L
r
、,
dr
=
置r.吖
,出m,,,≤o
(29)
为了求解权系数埘(r),可将式(29)离散化为矩阵形式:
G・W=C
(30)
oos一1(1一管)
O
6=2・
嘶一1(1—2等)嘲-l(1—2等)
10
o嗨一1(o)
啷一(o)
;00s一1(o)00B一1(0)o
(31)
W=【埘(尺)
埘(R—Ar)
…
埘(0)】T
(32)
C=【C,(一尺)C?
(一尺+Y)…cf(0)】T(33)式中,w为离散化的权系数向量,C为离散化的期望窗函数,G为离散化的积分矩阵。
则最优化权系数向量可求解为
W=G一1・C
(34)
万方数据
式中,G-1为G的逆矩阵。
假设期望的窗函数为Chebyshev窗,其旁瓣高度为一70dB旁瓣,可计算得到其对应的权系数。
图7优化加权分布函数
Fig.7
0ptilIlized
density
flln商Oll
图8一维方向图
rig.81一Dradiation
pattern
从图7可以看出,其形状与Chebyshev窗函数较为相似。
从图8可以看出,虽然该函数的峰值旁瓣比为一30.25dB,高于Chebyshev窗(但优于前面分析的各种圆形阵列),但在第二旁瓣处出现了陡降,达到了将近一50dB。
另外,由于存在旁瓣的陡降,其积分旁瓣比为一30.44dB,与峰值旁瓣比接近,要远远优于其它各种圆形阵列(如混合圆形阵列、圆形满阵元阵列等)。
5结论
通过引人分布函数,本文分析了二维阵列天线的方位图,可得到如下结论:
(1)二维阵列任意方向的方向图为其在该方向分布函数的离散傅里叶变换,该分布函数可通过统计其投影在某一方向的临近区域处阵元数目获得;
(2)该方法可较容易地分析各种二维阵列的方向图特征;
・95・
www.teleonline.Cln
电讯技术2010焦
(3)通过改变阵元权重系数即可改变阵列天线对应的分布函数,进而改变该阵列天线的方向图,该方法与通过改变阵列天线阵元分布的方法等效;
(4)利用分布函数可实现对圆形阵列阵元权重系数的优化,使其近似满足特定的窗函数,从而大大改善圆形阵列的波束性能。
参考文献:
I=.
NortonS
J.Syntheticapertureimagingwitharraysofarbitrary
sIl目Ipe:
PartI.Generalcase[J].IEEETramaedons
on
Ultra一
80hi08,FerroeleetricsandFrequencyControl,2002,49(4):
399—403.
匕KozickR
J,Kassam
SA.Coarraysynthesiswithcirelllarand
emptiealboundary
arrays[J].IEEE
Transactions
olrl
Image
Proce龉ing,1992,1(3):
391—405.b1J
Azevedo
JAR,CasimiroAMES.NOII-unifonllsamplingand
polyn,,mialinterpolationforarray
synth商8[J].IET
bli—
crowllves,Antennas&Propagation,2007,1(4):
867—873.
I!
-1J
TruccoA.MurinoV.Stochasticoptimizationoflinear
8pal'鸵
arrays[J].IEEEJoumalof
OC,e痢C
E画neeri.g,1999,24
(3):
291—299.B
31"uccoA.311inningandannealing[J].IEEE噩粼60璐On
weightingoflargeplanartmyS
by
simulated
Ultrasonics,
FerroelectricsandFrequencyControl,1999,46(2):
347—355.
№
Fen-eLm
JA,AresF.Patternsynthesisofeonformalarraysby
thesimulatedannealingtechnique[J].Electronics
letters,
1997,33(14):
1187一1189.一TrueeoA.Synlhesifing
aSylll腑ic
beam
patterns[J].IEEE
JoumalofOceanicE,lgi,leeri,lg,2000,25(3):
347—350..竺l1J
BevelacquaPJ。
BalanisCA.OotimizingAntennaArrayC,e—omelryforInterference
Suppression[J].IEEE
Transactions
On
Antennasand
Propagation,2007,55(3)・637—641.
b1J
PetkoJS,Werner
DH.1heevolution
0f删linear
polyfrac-
talarraysI玛inggenetic
al鲥tll瞄[J].IEEEⅢm践Iclio惜oll
Anternlsand
Pll3滞ltion,2005,55(11)"-3604—3615.
万方数据
叫
BrayM
G,WernerDH,BoeringerDW,eta1.Optimiz日lion
ofthinnedaperiodiclinearplIased
arraysusinggeneticalgo・
fithms
to
reduce
grati,lglobesduringscanning[J].IEEE
Tramaetions
Oil
Antennasand
Propagation,2002,50(12):
1732一1742.
Ⅲ
AUardR
J,WernerDH,WernerPL.Radiationpatternsyn—
thesisforarraysofeonformalantennas
mounted
on
arbitrari・
ly—shaped
three-dimensional
platformsusiI唱geneticalgo—
rithms[JJ.IEEETramactiom
On
AntennasandPropaga・
lion,2003,5l(5):
1054—1062.
例
CHENKe-song,YUNXiao-hua,HEZ,eta1.SynthesisofSparsePlanarArraysUsingModifiedRealGeneticAlgorithm
[J].IEEE
Transactions
013
Antennasand
h删∞,
2007,55(4):
1067—1073.
Ⅲ
Ares-PenaF
J,llodriguez-ComalezJA,Villanueva-I_opez
E,eta1.Geneticalgorithmsinthedesign
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 天线方向图 分布 函数 优化 方法 概要