CN114357767B - 一种宽带宽角波束扫描的和差波束稀疏阵列综合方法 - Google Patents

一种宽带宽角波束扫描的和差波束稀疏阵列综合方法 Download PDF

Info

Publication number
CN114357767B
CN114357767B CN202210004076.6A CN202210004076A CN114357767B CN 114357767 B CN114357767 B CN 114357767B CN 202210004076 A CN202210004076 A CN 202210004076A CN 114357767 B CN114357767 B CN 114357767B
Authority
CN
China
Prior art keywords
array
sum
difference
constraint
optimization
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202210004076.6A
Other languages
English (en)
Other versions
CN114357767A (zh
Inventor
刘颜回
刘泱
陈礼阳
杨仕文
胡俊
聂在平
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN202210004076.6A priority Critical patent/CN114357767B/zh
Publication of CN114357767A publication Critical patent/CN114357767A/zh
Application granted granted Critical
Publication of CN114357767B publication Critical patent/CN114357767B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • HELECTRICITY
    • H01ELECTRIC ELEMENTS
    • H01QANTENNAS, i.e. RADIO AERIALS
    • H01Q21/00Antenna arrays or systems
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/04Constraint-based CAD
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

本发明提供了一种宽带宽角波束扫描的和差波束稀疏阵列综合方法,涉及阵列天线方向图综合以及天线阵元位置优化。该方法包括如下步骤:步骤1.引入“设计频率”概念,将宽带宽角波束扫描的和差波束阵列综合问题转换为基于单一频点的和差波束阵列综合问题;步骤2.采用升幂级数(RPS)解析法对阵元进行初始布局;步骤3.在“设计频率”采用迭代凸优化方法进行多步阵元位置微扰优化,并在优化过程中设置多个约束条件,例如和波束及差波束的副瓣电平约束、目标方向上的差波束斜率约束、阵元位置微扰量幅值约束和阵元最小间距约束等。本发明可通过优化阵元布局同时生成具有宽带宽角波束扫描性能的和波束及差波束,并设置多个约束条件以满足实际的应用需求。

Description

一种宽带宽角波束扫描的和差波束稀疏阵列综合方法
技术领域
本发明涉及天线,尤其是涉及一种宽带宽角波束扫描的和差波束稀疏阵列综合方法。
背景技术
和差波束阵列天线在雷达和通信系统中对目标的精确跟踪和探测具有重要的作用,因而引起了广泛的关注。关于和差波束的阵列综合问题已经存在许多研究,大多数方法是通过优化激励分布或者阵元位置实现和波束与差波束的低副瓣辐射性能或者简化波束形成网络的设计复杂度。文献1和文献2分别采用Dolph Chebyshev和Bayliss权向量法通过两个不相关的激励分别综合和差波束,但是该方法的馈电网络设计复杂。文献3提出了两段波束成形网络方法,通过优化右半部分的激励分布,再对阵列的左右两部分以同相或反相馈入,由两种辐射模式产生和波束和差波束,但设计自由度受到较大的限制。文献4采用了子阵划分与进化差分算法结合的混合方法优化各个子阵的激励分布,从而简化馈电网络设计的复杂度和减少数字通道数量,但是只用于设计工作于单一频点的和差波束。
但是,上述这些方法无法抑制当和差波束阵列天线处于较宽的工作频段以及波束宽角扫描时所引入的栅瓣电平。众所周知,宽带阵列天线由于其较宽的工作频段,会导致栅瓣的出现以及副瓣尖峰电平的升高;特别是,对于生成和差波束的宽带阵列天线,需要同时对和波束及差波束的副瓣尖峰电平和栅瓣进行抑制。为了解决上述问题,常用的解决方法分别是紧耦合技术以及稀疏阵列技术,但是采用紧耦合技术将导致系统成本高、阵列通道数量多和存在大有源驻波比等缺点。相反稀疏阵列技术通过优化阵元布局可实现栅副瓣电平的抑制,具有散热性能较好、有效驻波比较小且减少系统成本等优点。稀疏阵列技术是指阵元之间不再以工作频点的半波长为单元间距,而是根据某种规则进行布阵,这使得稀疏阵列天线具有波束变窄、空间角分辨率高和阵元之间的互耦减弱等特点,这些特点促使阵列天线大量用于雷达、声纳等系统。目前,在和差波束阵列综合中,采用稀疏阵列技术以简化波束成形网络,而将稀疏阵列技术实现具有宽带及宽角波束扫描性能的和差波束研究却很少。
中国专利CN201611199282.8公开了一种稀布阵天线及其设计方法,该专利首先将天线阵面分为四个象限,每个象限划分为不同的子阵,每个象限的每个子阵与其它三个象限的对应子阵利用和差器得到子阵级和差波束。这些子阵级和差波束经过分别加权后再合成,就得到整个阵列的和差波束。但是该方法所综合的和差波束仍然不具有宽带及宽角波束扫描性能。
参考文献:
[1]Dolph C L.A current distribution for broadside arrays whichoptimizes the relationship between beam width and side-lobe level[J].Proceedings of the IRE,1946,34(6):335-348.
[2]Bayliss E T.Design of monopulse antenna difference patterns withlow sidelobes[J].Bell System Technical Journal,1968,47(5):623-650.
[3]Kwak S,Chun J,Park D,et al.Asymmetric sum and difference beampattern synthesis with a common weight vector[J].IEEE Antennas and WirelessPropagation.
[4]Letters,2016,15:1622-1625.Cui C Y,Jiao Y C,Li Z,et al.Synthesis ofsubarrayed monopluse arrays with contiguous elements using a DE Algorithm[J].IEEE Transactions on Antennas and Propagation,2017,PP(8):1-1.
发明内容
针对上述现有技术的不足,本发明提供一种宽带宽角波束扫描的和差波束稀疏阵列综合方法,通过优化后的阵元布局可同时生成具有宽带宽角波束扫描性能的和波束及差波束,并在优化过程中同时考虑和波束及差波束的副瓣电平约束,目标方向上的差波束斜率约束、阵元位置微扰量幅值约束、阵元最小单元间距约束和阵列孔径约束,从而满足实际的应用需求。
实现本发明目的的技术思路是:首先,引入“设计频率”概念,将宽带宽角波束扫描的和差波束阵列综合问题转换为基于单一频点的和差波束阵列综合问题。其次,采用升幂级数 (RPS)解析法对阵元进行初始布局,并且阵元位置关于阵列相位中心对称,从而仅以右半部分的阵元布局即可表示整个阵列的和波束及差波束阵因子。最后,为同时抑制和波束与差波束的最大副瓣电平及高频栅瓣电平且满足实际应用需求,以右半部分的阵元位置为优化变量,在“设计频率”采用迭代凸优化方法进行多步位置微扰优化,并在优化过程中嵌入多个约束条件,例如和波束及差波束的副瓣电平约束,目标方向上的差波束斜率约束、阵元位置微扰量幅值约束、阵元最小单元间距约束和阵列孔径约束。
为达上述技术目的,本发明采用如下技术方案予以实现。
一种宽带宽角波束扫描的和差波束稀疏阵列综合方法,包括以下步骤:
步骤1,为了减少时间和计算成本,将宽带扫描方向图转换为设计频率fD=RfRufL下的窄带侧射方向图从而进行综合;其中,Rf=fH/fL宽带比用来衡量工作频段的最高频率fH与最低频率fL之比;由波束指向角θ0∈[π/2-θmax,π/2+θmax]可以得到
Figure RE-RE-GDA0003547400530000021
采用Ru=max{u}/max{u0}=(1+|sinθmax|)波束扫描范围比衡量波束扫描后进入可见空间的增量,θmax是波束扫描范围的最大角度,u0∈[-1,1]是侧射波束可视空间的范围;设计频率下的侧射方向图包含了宽带扫描方向图的全部信息,即整个工作频率范围内所有波束扫描角的全部波束特征。
步骤2,以工作频段最低频的半波长为最小单元间距,采用基于升幂级数(RPS)解析法生成关于坐标原点对称且沿Z轴分布的2N个阵元初始位置,再通过选择合适的调节参数提升阵列的宽带性能。由此所生成的右半部分阵列单元和左半部分阵列单元位置表示为Rn与Ln。基于“设计频率”,以右半部分阵列单元位置Rn对右半部分阵列和左半部分阵列的阵因子分别作和运算以及差运算,则可得到和波束与差波束,其阵因子表达式分别如式(1)和(2) 所示:
Figure RE-RE-GDA0003547400530000022
Figure RE-RE-GDA0003547400530000023
其中,βD=2πfD/c表示基于“设计频率”的自由空间波数,u=cosθ-cosθ0,θ表示从x轴观测的角度,θ0为和波束及差波束的指向角,Rn表示右半部分阵列单元的位置。
步骤3,为了进一步降低栅瓣/副瓣区域的峰值电平,采用多步数值微扰优化法,所述多步数值微扰优化法的具体步骤为:将式(1)和(2)中阵元位置Rn(n=0,1,2,...,N)替换为Rnn,δn代表每一个阵元的位置微扰量。采用一阶泰勒近似sin(u0βDδn)≈u0βDδn,cos(u0βDδn)≈1对式(1)和(2)展开可得式(3)和(4),其中|βDδn|≤1。则式(1)和(2)可改写为式(3)和(4):
Figure RE-RE-GDA0003547400530000031
Figure RE-RE-GDA0003547400530000032
步骤4,为满足实际的需求,在优化过程中需要同时考虑和波束及差波束的多个约束条件,具体如下:
(1)和波束的副瓣电平约束:为了保证在优化过程中和波束的副瓣电平不超过指定的范围,引入辅助变量ε1约束和波束在副瓣区域
Figure RE-RE-GDA0003547400530000033
的取值范围为[-ε1,ε1]。
(2)差波束的副瓣电平约束:为了保证在优化过程中差波束的副瓣电平不超过指定的范围,引入辅助变量ε2约束差波束在副瓣区域
Figure RE-RE-GDA0003547400530000034
的取值范围为[-ε2,ε2]。
(3)最小间距约束:在实际天线阵布局过程中,阵元间距需不小于所设计的单元天线尺寸,因此在微扰过程中,需要满足最小间距约束dn-dn-1≥dmin
(4)位置微扰量幅值约束:为了保证位置微扰的近似精度,要求δn(n=0,1,2,…,N)均满足 -μ≤βDδn≤μ,此处的μ决定了位置微扰的幅度;
(5)阵列孔径约束:在某些可用空间受限的应用中,阵列空间约束的存在十分有必要。该约束可以表示为2(RNN)≤A,其中A为可接受的最大阵列孔径;
(6)目标方向上的差波束斜率约束:对于目标检测和追踪的应用中,需要维持一个较窄和深的零点,差波束斜率被要求尽可能大以及零点深度尽可能小,该约束表示为
Figure RE-RE-GDA0003547400530000035
步骤5,以右半部分的阵元位置为优化变量,在“设计频率”采用迭代凸优化方法进行多步位置微扰优化,所设置的约束条件如步骤4所述,优化目标为最小化和波束及差波束在副瓣区域的尖峰电平之和,表达式为:
min(ε12)
Figure RE-RE-GDA0003547400530000041
在迭代凸优化过程中,ε1表示和波束在副瓣区域的峰值电平,ε2表示差波束在副瓣区域的峰值电平,
Figure RE-RE-GDA0003547400530000042
表示和波束的副瓣区域,
Figure RE-RE-GDA0003547400530000043
表示差波束的副瓣区域,η表示差波束斜率的最小值,优化目标是最小化和波束与差波束副瓣区域的峰值电平之和(ε12)。阵元布局的多步微扰后,当栅瓣/副瓣区域的峰值电平之和收敛到最小值时,且阵列布局也保持不变时,得到最优化阵列布局结果。
与现有技术相比,本发明的有益效果在于:通过引入“设计频率”概念,采用解析法与迭代凸优化法相结合的阵列布局优化方法,并在优化过程中同时考虑和波束及差波束的副瓣电平、差波束在目标方向上的斜率及最小单元间距等约束条件。可实现和差波束稀疏阵列在较宽的工作频段以及波束扫描范围较大时的栅瓣/副瓣电平抑制以及降低了计算复杂度,并满足了实际的应用需求。在对阵元布局的优化完成后,即可同时生成具有宽带宽角波束扫描性能的和波束及差波束。
附图说明
图1为本发明的技术方案流程图。
图2为不同r值下的和波束与差波束的栅瓣/副瓣峰值电平之和。
图3为阵列布局优化过程中和波束与差波束的栅瓣/副瓣峰值电平之和。
图4为最优化稀疏阵列正半轴的阵元位置图。
图5为基于RPS方法布局的和波束在3GHz时波束指向135度的扫描方向图。
图6为基于RPS方法布局的差波束在3GHz时波束指向135度的扫描方向图。
图7为最优化稀疏阵列的和波束在3GHz时波束指向135度的扫描方向图。
图8为最优化稀疏阵列的差波束在3GHz时波束指向135度的扫描方向图。
图9为最优化稀疏阵列的和波束在2GHz时波束指向135度的扫描方向图。
图10为最优化稀疏阵列的差波束在2GHz时波束指向135度的扫描方向图。
图11为最优化稀疏阵列的和波束在1GHz时波束指向135度的扫描方向图。
图12为最优化稀疏阵列的差波束在1GHz时波束指向135度的扫描方向图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施案例,并参照附图,对本发明进一步详细说明。
下面结合一个实例对本发明做进一步解释。设计目标:频率范围在1GHz~3GHz之间,波束扫描范围为[45°,135°],栅瓣/副瓣区域的峰值电平保持不变的和差波束无栅瓣稀疏线阵。
按下面步骤实施:
步骤1:由于fL=1GHz和fH=3GHz,采用宽带比Rf=fH/fL=3衡量工作频段的最高频率与最低频率之比;由于波束指向角θ0∈[π/2-θmax,π/2+θmax],θmax=45°且 u∈[-1-|sinθmax|,1+|sinθmax|],采用波束扫描距离比 Ru=max{u}/max{u0}=(1+|sinθmax|)=1.707衡量波束扫描后进入可见空间的增量,其中 u0=cosθ,u0∈[-1,1]是侧射波束可视空间u的范围。为了减少时间和计算成本,将对整个工作频带以及波束扫描范围内的方向图进行迭代优化转换为对设计频率 fD=RfRufL=5.1213GHz下的阵列方向图进行综合。由于设计频率下的侧射方向图包含了阵列综合所需要的全部信息,即整个工作频率范围内所有波束扫描角的全部波束特征。
步骤2:在笛卡尔坐标系中,以工作频段最低频的半波长(0.15m)为最小单元间距dmin=λL/2,采用基于升幂级数(RPS)解析法生成关于坐标原点对称且沿Z轴分布的2N个阵元初始位置。同时,为了避免过度稀疏化,调节参数r值不可选取的过大或过小,因此 0.75≤r≤1.25。为了选取合适的r值,由图2得知,当r取1.25时,在整个工作频率范围[1GHz,3GHz]及波束扫描范围[45°,135°]内,和波束与差波束副瓣峰值电平之和最小 -16.62dB。因此,r的取值为1.25,基于设计频率下的和波束与差波束的阵因子表达式为:
Figure RE-RE-GDA0003547400530000051
Figure RE-RE-GDA0003547400530000052
其中,从X轴观测的角度为θ,采样密度△θ是0.02°,βD=2πfD/c表示“设计频率”的自由空间波数。
步骤3,为了进一步降低栅瓣/副瓣区域的峰值电平,采用多步数值微扰优化法对阵元位置进一步优化,所述多步数值微扰优化法的具体步骤为:将式(1)和(2)中阵元位置Rn(n=0,1,2,...,N)替换为Rnn,δn代表每一个阵元的位置微扰量。采用一阶泰勒近似sin(u0βDδn)≈u0βDδn,cos(u0βDδn)≈1对式(1)和(2)展开可得式(3)和(4),其中|βDδn|≤1。则式(1)和(2)可改写为式(3)和(4):
Figure RE-RE-GDA0003547400530000053
Figure RE-RE-GDA0003547400530000061
步骤4,为满足实际的需求,在优化过程中需要同时考虑和波束及差波束的多个约束条件,具体如下:
(1)和波束的副瓣电平约束:为了保证在优化过程中和波束的副瓣电平不超过指定的范围,引入辅助变量ε1约束和波束在副瓣区域
Figure RE-RE-GDA0003547400530000062
的取值范围为[-ε1,ε1]。
(2)差波束的副瓣电平约束:为了保证在优化过程中差波束的副瓣电平不超过指定的范围,引入辅助变量ε2约束差波束在副瓣区域
Figure RE-RE-GDA0003547400530000063
的取值范围为[-τ2,τ2]。
(3)最小间距约束:在实际天线阵布局过程中,阵元间距需不小于所设计的单元天线尺寸,因此在微扰过程中,需要满足最小间距约束dn-dn-1≥dmin;本例中设置的最小单元间距约束为工作最低频点的半波长λL/2,即dmin=λL/2=0.15m。
(4)位置微扰量幅值约束:为了保证位置微扰的近似精度,要求即对于δn(n=0,1,2,…,N)均满足-μ≤βDδn≤μ,此处的μ决定了位置微扰的幅度;本例中相对于工作波长位置微扰为 1/20,即μ=π/20。
(5)阵列孔径约束:在某些可用空间受限的应用中,阵列空间约束的存在十分有必要。该约束可以表示为2(RNN)≤A,其中A为可接受的最大阵列孔径;本例中最大阵列孔径为 A=64.68λL
(6)目标方向上的差波束斜率约束:对于目标检测和追踪的应用中,需要维持一个较窄和深的零点,差波束斜率被要求尽可能大以及零点深度尽可能小,该约束表示为式(5);本例中η为52.59dB。
Figure RE-RE-GDA0003547400530000064
步骤5,遵循以上约束条件,以右半部分的阵元位置为优化变量,在“设计频率”采用迭代凸优化法对和波束阵列及差波束阵列的布局进行多步位置微扰优化,所设置的约束条件如步骤4所述,优化目标为最小化和波束及差波束在副瓣区域的尖峰电平之和,表达式为(6):
min (ε12)
Figure RE-RE-GDA0003547400530000071
在迭代凸优化过程中,ε1表示和波束在副瓣区域的峰值电平;ε2表示差波束在副瓣区域的峰值电平;
Figure RE-RE-GDA0003547400530000072
表示和波束的副瓣区域;
Figure RE-RE-GDA0003547400530000073
表示差波束的副瓣区域;η表示差波束斜率的最小值;μ表示位置微扰量幅值约束;表示在每一步迭代凸优化过程中,优化目标是最小化和波束与差波束最小化栅瓣/副瓣区域的峰值电平之和;大约40次迭代优化后,当栅瓣/副瓣区域的峰值电平之和收敛到最小值-22.04dB时,且阵列布局也保持不变时,得到最优化阵列布局结果。
图2给出了对于RPS方法最优的调节参数r的取值,r为1.25时和波束和差波束的最大副瓣电平之最小,为-16.62dB。图3给出了每步微扰结果的栅瓣/副瓣区域的峰值电平值,大约40次迭代后,栅瓣/副瓣区域的峰值电平降到大约-22.04dB,随后基本保持不变,通过迭代凸优化之后的和波束与差波束栅瓣/副瓣区域的峰值电平总共降低了5.42dB。图4给出了最优化的稀疏阵列布局图,由于阵元关于原点对称分布,这里仅画出正半轴的阵元。图5和图6是基于RPS方法在3GHz工作频点进行布局后的和波束与差波束阵因子图,和波束与差波束的最大副瓣电平为-8.23dB和-8.38dB。图7和图8分别是和波束与差波束的最优化稀疏阵列在3GHz时的阵因子图,此时的和波束与差波束的栅瓣/副瓣区域峰值电平之和为-22.04dB。图7~11给出了3个不同频率(1GHz、2GHz和3GHz)下的和波束与差波束阵因子图,和波束在 3个不同频点的最大副瓣电平为-10.96dB、-10.97dB和-10.96dB,差波束在3个不同频点的最大副瓣电平为-11.04dB、-11.02dB和-11.03dB。
以上所述实施例仅表达了本发明的具体实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形和改进,这些都属于本发明的保护范围。

Claims (6)

1.一种宽带宽角波束扫描的和差波束稀疏阵列综合方法,其特征在于,包括以下步骤:
步骤1.根据目标工作频段以及最大波束扫描范围,引入“设计频率”概念到和差波束的综合过程,将宽带波束可扫描的和差波束阵列综合问题转换为基于单一频点的和差波束阵列综合问题;
步骤2.确定最小单元间距,即工作频段最低频点的半波长,以及阵元数目2N,采用基于升幂级数解析法生成关于阵列相位中心呈对称分布的阵元初始位置,再通过选择合适的调节参数提升阵列的宽带性能,并以右半部分的阵元布局表示整个阵列的和波束及差波束阵因子;从而,仅以右半部分的阵元布局作为优化变量,即可实现和波束及差波束的宽带及宽角扫描性能的同时优化;
步骤3.采用微扰法对阵元位置进行优化,并通过一阶泰勒近似保证可凸优化求解;
步骤4.同时考虑多个和波束及差波束的实际设计要求,在优化过程中设置多个约束条件,包括和波束及差波束的副瓣电平约束,目标方向上的差波束斜率约束、阵元位置微扰量幅值约束、阵元最小单元间距约束和阵列孔径约束;
步骤5.为同时抑制和波束与差波束的最大副瓣电平及高频栅瓣电平,并满足设置的约束条件,以右半部分的阵元位置为优化变量,在“设计频率”采用迭代凸优化方法进行多步阵元位置微扰优化,多次迭代之后当和波束及差波束副瓣/栅瓣区域的峰值电平之和收敛到最小值时,得到最优化阵列布局结果。
2.根据权利要求1所述的一种宽带宽角波束扫描的和差波束稀疏阵列综合方法,其特征在于:所述步骤1具体如下:
为了减少时间和计算成本,将宽带扫描方向图转换为设计频率fD=RfRufL下的窄带侧射方向图从而进行综合;其中,Rf=fHfL宽带比用来衡量工作频段的最高频率fH与最低频率fL之比;由波束指向角θ0∈[π/2-θmax,π/2+θmax]可以得到
Figure FDA0004053784060000011
采用Ru=max{u}max{u0}=(1+sinθmax)波束扫描范围比衡量波束扫描后进入可见空间的增量,θmax是波束扫描范围的最大角度,u0∈[-1,1]是侧射波束可视空间的范围;设计频率下的侧射方向图包含了宽带扫描方向图的全部信息,即整个工作频率范围内所有波束扫描角的全部波束特征。
3.根据权利要求1所述的一种宽带宽角波束扫描的和差波束稀疏阵列综合方法,其特征在于:所述步骤2具体如下:
以工作频段最低频的半波长为最小单元间距,采用基于升幂级数解析法生成关于坐标原点对称且沿Z轴分布的2N个阵元初始位置,再通过选择合适的调节参数提升阵列的宽带性能,由此所生成的右半部分阵列单元和左半部分阵列单元位置表示为Rn与Ln;基于“设计频率”,以右半部分阵列单元位置Rn对右半部分阵列的阵因子AFR(u,βD)和左半部分阵列的阵因子AFL(u,βD)分别作和运算以及差运算,则可得到和波束与差波束,其阵因子表达式分别如式(1)和(2)所示:
Figure FDA0004053784060000021
Figure FDA0004053784060000022
其中,βD=2πfDc表示基于“设计频率”的自由空间波数,u=cosθ-cosθ0,θ表示从x轴观测的角度,θ0为和波束及差波束的指向角,Rn表示右半部分阵列单元的位置。
4.根据权利要求3所述的一种宽带宽角波束扫描的和差波束稀疏阵列综合方法,其特征在于:所述步骤3具体如下:
为了进一步降低栅瓣/副瓣区域的峰值电平,采用多步数值微扰优化法,所述多步数值微扰优化法的具体步骤为:将式(1)和(2)中阵元位置Rn(n=0,1,2,...,N)替换为Rnn,δn代表每一个阵元的位置微扰量;采用一阶泰勒近似sin(u0βDδn)≈u0βDδn,cos(u0βDδn)≈1对式(1)和(2)展开可得式(3)和(4),其中βDδn≤1;则式(1)和(2)可改写为式(3)和(4):
Figure FDA0004053784060000023
Figure FDA0004053784060000024
5.根据权利要求1所述的一种宽带宽角波束扫描的和差波束稀疏阵列综合方法,其特征在于:所述步骤4具体如下:
为满足实际的需求,在优化过程中需要同时考虑和波束及差波束的多个约束条件,具体如下:
(1)和波束的副瓣电平约束:为了保证在优化过程中和波束的副瓣电平不超过指定的范围,引入辅助变量ε1约束和波束在副瓣区域
Figure FDA0004053784060000025
的取值范围为[-ε1,ε1];
(2)差波束的副瓣电平约束:为了保证在优化过程中差波束的副瓣电平不超过指定的范围,引入辅助变量ε2约束差波束在副瓣区域
Figure FDA0004053784060000026
的取值范围为[-ε2,ε2];
(3)最小间距约束:在实际天线阵布局过程中,阵元间距需不小于所设计的单元天线尺寸,因此在微扰过程中,需要满足最小间距约束dn-dn-1≥dmin
(4)位置微扰量幅值约束:为了保证位置微扰的近似精度,要求δn(n=0,1,2,…,N)均满足-μ≤βDδn≤μ,此处的μ决定了位置微扰的幅度;
(5)阵列孔径约束:在某些可用空间受限的应用中,阵列空间约束的存在十分有必要;该约束可以表示为2(RNN)≤A,其中A为可接受的最大阵列孔径;
(6)目标方向上的差波束斜率约束:对于目标检测和追踪的应用中,该约束表示为
Figure FDA0004053784060000031
6.根据权利要求5所述的一种宽带宽角波束扫描的和差波束稀疏阵列综合方法,其特征在于:所述步骤5具体如下:
以右半部分的阵元位置为优化变量,在“设计频率”采用迭代凸优化方法进行多步位置微扰优化,所设置的约束条件如步骤4所述,优化目标为最小化和波束及差波束在副瓣区域的尖峰电平之和,表达式为:
min(ε12)
Figure FDA0004053784060000032
在迭代凸优化过程中,ε1表示和波束在副瓣区域的峰值电平,ε2表示差波束在副瓣区域的峰值电平,
Figure FDA0004053784060000033
表示和波束的副瓣区域,
Figure FDA0004053784060000034
表示差波束的副瓣区域,η表示差波束斜率的最小值,优化目标是最小化和波束与差波束副瓣区域的峰值电平之和(ε12);阵元布局的多步微扰后,当栅瓣/副瓣区域的峰值电平之和收敛到最小值时,且阵列布局也保持不变时,得到最优化阵列布局结果。
CN202210004076.6A 2022-01-04 2022-01-04 一种宽带宽角波束扫描的和差波束稀疏阵列综合方法 Active CN114357767B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210004076.6A CN114357767B (zh) 2022-01-04 2022-01-04 一种宽带宽角波束扫描的和差波束稀疏阵列综合方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210004076.6A CN114357767B (zh) 2022-01-04 2022-01-04 一种宽带宽角波束扫描的和差波束稀疏阵列综合方法

Publications (2)

Publication Number Publication Date
CN114357767A CN114357767A (zh) 2022-04-15
CN114357767B true CN114357767B (zh) 2023-05-12

Family

ID=81107757

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210004076.6A Active CN114357767B (zh) 2022-01-04 2022-01-04 一种宽带宽角波束扫描的和差波束稀疏阵列综合方法

Country Status (1)

Country Link
CN (1) CN114357767B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115084874B (zh) * 2022-07-19 2023-12-26 电子科技大学 基于异构子阵非均匀布局的波束扫描阵列优化设计方法
CN116315735B (zh) * 2023-03-07 2023-12-05 电子科技大学 一种任意极化可编程共形阵列天线构建及综合方法
CN117634115B (zh) * 2024-01-26 2024-04-16 安徽大学 一种小型化天线阵列方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4874917B2 (ja) * 2007-09-28 2012-02-15 株式会社東芝 アレーアンテナ、アンテナ素子配置方法及びビーム形成方法
CN106650260A (zh) * 2016-12-22 2017-05-10 厦门大学 最小间距可控的超宽带无栅瓣稀疏线阵设计方法
CN113032989B (zh) * 2021-03-19 2021-10-22 电子科技大学 一种约束方向性系数的最小化阵元数目的稀布阵优化方法

Also Published As

Publication number Publication date
CN114357767A (zh) 2022-04-15

Similar Documents

Publication Publication Date Title
CN114357767B (zh) 一种宽带宽角波束扫描的和差波束稀疏阵列综合方法
Wu et al. Proactive conformal antenna array for near-field beam focusing and steering based on curved substrate integrated waveguide
CN105789877B (zh) 基于超表面的四波束微带透射阵天线及其设计方法
Mailloux Electronically scanned arrays
CN107422569B (zh) 一种二维光学相控阵
US10566698B2 (en) Multifocal phased array fed reflector antenna
Poveda-García et al. Ambiguity resolution in amplitude-monopulse systems using broad-beam patterns
CN113252998B (zh) 相控阵天线和、差波束信号电平的平坦度优化方法
CN109346843B (zh) 一种空间一维扫描透镜天线的设计方法及波束扫描方法
CN116244940A (zh) 一种双波束超宽带阵列天线优化布局方法
CN110854547B (zh) 一种阵列馈电式大范围波束扫描反射面天线
CN110176665B (zh) 一种相控阵天线及相控阵雷达
Giusto et al. Phase-only optimization for the generation of wide deterministic nulls in the radiation pattern of phased arrays
CN108242600B (zh) 一种线极化单脉冲平板缝隙天线
CN113489523A (zh) 基于分离校准迭代fft的唯相位多波束方向图综合方法
CN112949193A (zh) 一种子阵级稀疏阵列天线方向图数值方法及系统
Mohammed Simplified rectangular planar array with circular boundary for side lobe suppression
CN112558065B (zh) 一种基于可重构电磁表面阵列的三维成像方法
CN115133291A (zh) 非规则天线子阵、相控阵天线及相控阵天线的设计方法
CN113113784A (zh) 超大间距阵列无栅瓣大角度扫描布阵方法
CN111291493B (zh) 一种用于机载预警共形阵俯仰探测波束赋形的设计方法
CN113656747B (zh) 基于分支定界的多期望信号下的阵列自适应波束形成方法
Hassan et al. Comparative study of different power distribution methods for array antenna beamforming for soil moisture radiometer
CN117895245B (zh) 基于ga-sqp算法的反射面天线设计方法
Lee A simple method of dual‐reflector geometrical optics synthesis

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant