CN111551923B - 一种多约束下的均匀线性阵列低旁瓣波束形成优化方法 - Google Patents
一种多约束下的均匀线性阵列低旁瓣波束形成优化方法 Download PDFInfo
- Publication number
- CN111551923B CN111551923B CN202010459360.3A CN202010459360A CN111551923B CN 111551923 B CN111551923 B CN 111551923B CN 202010459360 A CN202010459360 A CN 202010459360A CN 111551923 B CN111551923 B CN 111551923B
- Authority
- CN
- China
- Prior art keywords
- array
- constraint
- representing
- algorithm
- uniform linear
- 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
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 36
- 238000000034 method Methods 0.000 title claims abstract description 22
- 230000014509 gene expression Effects 0.000 claims description 8
- 230000015572 biosynthetic process Effects 0.000 claims description 7
- 230000003190 augmentative effect Effects 0.000 claims description 5
- 239000011159 matrix material Substances 0.000 claims description 4
- 230000000694 effects Effects 0.000 abstract description 8
- 238000012545 processing Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 6
- 238000000342 Monte Carlo simulation Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 4
- 230000009977 dual effect Effects 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 3
- 238000013461 design Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000003491 array Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000002452 interceptive effect Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/40—Means for monitoring or calibrating
- G01S7/4004—Means for monitoring or calibrating of parts of a radar system
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S2013/0236—Special technical features
- G01S2013/0245—Radar with phased array antenna
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明公开一种多约束下的均匀线性阵列低旁瓣波束形成优化方法,应用于阵列信号处理领域,针对现有的ADMM算法求解该类优化问题时惩罚因子初值的选择对算法收敛速度影响较大,易出现算法收敛慢或不收敛的情况;本发明首先建立均匀线阵接收信号模型,以阵列增益作为目标函数,考虑主瓣电平、旁瓣电平和权值幅度约束,然后采用交替方向惩罚算法最大化阵列增益,求解最优阵列天线加权系数;相比于现有的ADMM算法收敛速度更快,收敛效果更好。
Description
技术领域
本发明属于阵列信号处理领域,特别涉及一种多约束下的均匀线性阵列低旁瓣波束形成优化技术。
背景技术
波束形成是一种空域滤波方法,它使用一组传感器来增强所需方向的信号,同时抑制来自空间不同方向的干扰信号,故关于波束形成的优化问题得到了广泛的研究。通过设计阵列天线的加权系数得到的期望的波束方向图通常关注于主瓣电平和旁瓣电平,从而实现对来自空间所需方向外不同方向的干扰信号的抑制;此外,在实际工程应用中,阵列天线的加权系数的优化设计需要加入硬件器件带来的权值幅度的约束,并考虑接收端阵列增益(即输出信噪比与输入信噪比的比值)需求。
为了通过设计阵列天线的加权系数得到低旁瓣的波束方向图,有一种波束形成的优化方法以峰值旁瓣电平作为目标函数,同时考虑主瓣电平约束和权向量恒模约束,采用交替方向乘子(Alternate Direction Method of Multipliers,ADMM)算法最小化代价函数求解得出对应的阵列天线加权系数,由此得到期望的波束方向图。然而,这种方法没有考虑到实际工程中对阵列增益的要求;同时,ADMM算法求解该类优化问题时惩罚因子初值的选择对算法收敛速度影响较大,易出现算法收敛慢或不收敛的情况。
发明内容
为解决上述技术问题,本发明提出一种多约束下的均匀线性阵列低旁瓣波束形成优化方法,以阵列增益作为优化标准建立目标函数,考虑阵列天线硬件实际情况和需求,加入主瓣电平、旁瓣电平和权值幅度的约束,保证了在满足低旁瓣波束形成的同时使得阵列增益最大。
本发明采用的技术方案为:一种多约束下的均匀线性阵列低旁瓣波束形成优化方法,包括:
S1、建立均匀线阵接收信号模型;
S2、以阵列增益作为目标函数;
S3、以主瓣电平、旁瓣电平和权值幅度为约束,建立优化问题;
S4、采用交替方向惩罚(Alternating Direction Penalty Method,ADPM)算法最大化阵列增益,求解优化问题得到最优阵列天线加权系数。
步骤S2所述目标函数表达式为:
其中,w表示权向量,am表示导向矢量,(·)H表示共轭转置,R表示噪声协方差矩阵;
阵列增益的定义为:
其中,E(·)表示求均值,(·)H表示共轭转置,|·|表示复数的模,SNRout和SNRin分别表示输出信噪比和输入信噪比,x(t)和n(t)分别表示输入的有用信号和噪声信号,E(|·|2)表示信号的平均功率。
步骤S3所述主瓣电平约束表达式为:
Pl≤|wHam|2≤Pu
其中,Pl和Pu分别表示主瓣电平约束的下界和上界;
旁瓣电平约束表达式为:
|wHas|2≤η,s=1,…,S
其中,as=a(θs),s=1,…,S,S表示将所观察的空域划分为若干离散的角度后,除了波束主瓣所在方向区域外旁瓣区域所有角度的总数,θs表示旁瓣方向,η表示旁瓣电平约束的上界;
权值幅度约束为:
σ≤|wn|≤1,n=1,…,N
其中,σ=1/2M表示权值幅度取值范围的下界,N表示天线阵元的个数,M表示衰减器的比特位数。
步骤S3所述优化问题为:
步骤S4具体为:引入辅助变量与两个惩罚因子,构建增广拉格朗日函数;将优化问题分解为三个并行的子优化问题,每个子问题中只包含一个约束条件;交替方向惩罚规则用于逐次迭代运算中,直到满足迭代停止条件;得到最优阵列天线加权系数。
引入的两个惩罚因子根据相应的残差进行自适应改变。
所述两个惩罚因子的迭代表达式为:
本发明的有益效果:本发明提供了一种多约束下的均匀线性阵列低旁瓣波束形成优化方法;针对低旁瓣波束形成问题,本发明建立了均匀线阵的接收信号模型,以阵列增益作为优化标准建立目标函数,考虑阵列天线硬件实际情况和需求,加入主瓣电平、旁瓣电平和权值幅度的约束,保证了在满足低旁瓣波束形成的同时使得阵列增益最大;相对于原有的求解低副瓣波束形成的ADMM算法,本发明所采用的交替方向惩罚(ADPM)算法收敛速度快,收敛效果好。
附图说明
图1为本发明实施例的流程图。
图2为图1中步骤S4的流程图。
图3为本发明实施例的阵列方向图对比图。
图4为本发明实施例的惩罚因子变化曲线图。
图5为本发明实施例的阵列增益对比图。
图6为本发明实施例的收敛曲线对比图。
图7为本发明实施例的残差Δry变化曲线对比图。
图8为本发明实施例的残差Δrm变化曲线对比图。
具体实施方式
为便于本领域技术人员理解本发明的技术内容,下面结合附图对本发明内容进一步阐释。
如图1所示,本发明提出的一种多约束下的均匀线性阵列低旁瓣波束形成优化方法,通过以下步骤实现:
S1、均匀线阵的相控阵雷达包含N个天线阵元,建立接收信号模型。
本实施例中,有N个各向同性辐射阵元组成的均匀线阵,其阵元间隔为d,其权向量w表示为:
w=[w1,w2,…,wN]T
其中,(·)T表示转置,wn表示第n个阵元的复加权系数,n=1,…,N。则对于远场情况,对于给定的θ方向的波束函数表示为
y(θ)=wHa(θ)
其中,(·)H表示共轭转置;a(θ)表示θ方向的导向矢量,可以表示为
S2、建立目标函数。
本实施例中,以阵列增益作为波束形成设计准则,阵列增益取决于输出和输入信噪比的比值,阵列增益的定义为:
其中,E(·)表示求均值,(·)H表示共轭转置,|·|表示复数的模,E(|·|2)表示信号的平均功率,SNRout和SNRin分别表示输出信噪比和输入信噪比,x(t)和n(t)分别表示输入的有用信号和噪声信号,主瓣方向导向矢量am=a(θm),θ表示角度,θm表示主瓣波束的方向,R表示噪声协方差矩阵;因此,考虑零均值互不相关高斯白噪声,目标函数可以表示为:
S3、加入主瓣电平、旁瓣电平和权值幅度约束,满足三个约束的同时最大化目标函数,建立优化问题。
本实施例中,步骤S3通过以下子步骤实现:
S31、计算主瓣电平约束;
当给定的目标方向不能准确确定时,有必要使得主瓣电平在一定范围内波动。因此有主瓣电平约束:
Pl≤|wHam|2≤Pu
其中,Pl和Pu分别表示主瓣电平约束的下界和上界,技术人员可以根据应用需要人为设定。在实际设置中考虑以下三点:
1.对于设计不同波束指向的波束形成方向图时,主瓣电平上界和下界设定一致,并且上界和下界之间的波动范围要小,保证方向图中波束指向处的峰值电平在不同指向波动较小,有利于探测目标的方向;
2.通常希望设定的主瓣波动范围接近理想情况下(常规波束形成:相位加权)主瓣电平的最大值;
3.设定主瓣电平的范围,便于衡量副瓣电平的效果(低副瓣的效果是看主瓣电平和副瓣电平之间的差值)。
S32、计算旁瓣电平约束;
低旁瓣波束形成要求得到的天线方向图中旁瓣电平在给定的电平以下,因此旁瓣电平约束为:
|wHas|2≤η,s=1,…,S
其中,η表示旁瓣电平约束的上界,旁瓣方向导向矢量as=a(θs),s=1,…,S,θs表示旁瓣方向,{s=1,2,…,S}为将所观察的空域划分为若干离散的角度后,除了波束主瓣所在方向区域外旁瓣区域所有角度,共S个角度。
S33、计算权值幅度约束;
由于加权系数的幅度受到硬件器件中的衰减器的影响,每个权值的幅度范围受衰减器的比特位数的限制。因此,权值幅度约束为:
σ≤|wn|≤1,n=1,…,N
其中,σ=1/2M表示权值幅度取值范围的下界,M表示衰减器的比特位数,权向量w=[w1,w2,…,wN]T共有N个权值,wn表示第n个复数权值,|wn|取模运算得到幅度。
S34、建立优化问题;
低旁瓣波束形成优化问题为满足所述主瓣电平、旁瓣电平和权值幅度约束的同时,最大化所述目标函数,求解阵列的最优加权系数,优化问题P为:
S4、求解优化问题,得到阵列优化加权系数。
本实施例中,引入辅助变量,采用的ADPM算法架构,构建增广拉格朗日函数;原优化问题分解为三个并行的子优化问题,每个子问题中只包含一个约束条件,各子问题易分别求解;引入的两个惩罚因子根据相应的残差进行自适应改变;所述交替方向惩罚规则用于逐次迭代运算中,直到满足迭代停止条件;基于所述交替方向惩罚算法求解所述优化问题。
如图2所示,ADPM算法根据以下子步骤实现:
S41、引入辅助变量,构建增广拉格朗日函数;参数初始化;
由此构建增广拉格朗日函数:
其中,和分别表示复数的实部和虚部,αr和αi分别为约束的对偶辅助变量的实部和虚部,上标r,i用于区分实部和虚部,和分别为约束ys-wHas=0,s=1,…,S的对偶辅助变量的实部和虚部,ρm>0和ρy>0表示可以自适应的变化的惩罚因子。
S42、原优化问题分解为三个并行的子优化问题,每个子问题中只包含一个约束条件,依次更新各子问题中的变量、惩罚因子以及对偶辅助变量;
本实施例中,预设最多进行Kmax次迭代,Kmax的值可以人为设定,通常将该值设的较大,比如Kmax=104,从而保证找到收敛的解,假设第k次迭代后得到w(k),χ(k),则第(k+1)次迭代结果按照下列公式计算:
其中,δ1<1和δ2>1都是很接近1的数,ν是一个设定的较大的数,比如ν=105,和表示残差。从和的更新公式中可以看出,如果残差和不随着迭代更新逐渐减小,惩罚因子和就会逐渐增大迫使其减小。如果和两个判断条件在迭代过程中总是满足,则惩罚因子保持不变,可以看作是ADMM算法。ADPM算法中惩罚因子可以根据残差的变化自适应的改变,使得本发明的方法比ADMM算法更易找到可行解。
针对更新w(k+1)的子优化问题P1,化简并删去与w(k+1)无关的项,w(k+1)可以按照下式求解得到:
其中,||·||2表示二范数,并且有
其中,0表示1×N的行向量。由此,求解出P1′目标函数的解为w*=(AAH)-1AbH。因此,w(k+1)的解为:
P3″是关于χ的一元问题,因此,将区间[Pl,Pu]分成Kχ个离散的点,易在P3″目标函数中搜索得到χ(k+1)。
S43、判断是否满足预设终止条件,若不满足,回到所述步骤S42;若满足,进入步骤S44;
如果迭代次数超过预设的最大迭代次数Kmax,或者残差同时满足和则迭代停止。其中,εm和εy分别为预设的最大容忍残差,当残差小于预设的最大容忍误差时,认为算法是收敛的,可以得到符合约束条件的解。通常最大容忍误差较小,例如可以设定为εy=εm=10-5。
S44、输出满足迭代停止条件时得到的加权系数w,为该问题阵列最优加权系数。
以下为本发明的两个具体实施例。
实施例一:均匀线阵阵元个数为16,阵元间隔为半波长,波束指向0°,主瓣电平的上下界电平范围为Pl=21.78dB和Pu=21.98dB,旁瓣电平的上界电平为η=0.083dB,权值幅度范围的下界为σ=1/27,最大迭代次数为Kmax=10000,最大容忍残差为εy=εm=10-5,惩罚因子的初始值为ρy=ρm=10-4,对偶辅助变量的取值不超过ν=500000,δ1=0.99,δ2=1.02,其他参数初始化在[0,1]之间随机产生。考虑噪声协方差矩阵为R=I,则在没有其他约束条件下,最大阵列增益为12.0411dB。比较本发明采用的ADPM算法与交替方向乘子(ADMM)算法的优化效果,惩罚因子固定为ρy=ρm=10-4,其他条件保持相同。
图3为采用ADPM和ADMM两种算法求解的最优加权系数对应的方向图;图4为采用ADPM算法迭代求解过程中惩罚因子ρm,ρy的自适应变化结果;图5为进行500次蒙特卡罗仿真实验,采用ADPM和ADMM两种算法求解的最优加权系数对应的阵列增益结果;图6为进行500次蒙特卡罗仿真实验,采用ADPM和ADMM两种算法求解最优加权系数过程进行迭代的次数。图3描述的是两种算法得到的最优加权系数对应的方向图,都满足主瓣电平和旁瓣电平的约束条件。图5描述的是经过500次蒙特卡罗仿真实验两种算法得到的最大阵列增益的结果,说明ADPM算法求解得到最大阵列增益的效果比ADMM算法好。图6描述的是经过500次蒙特卡罗仿真实验两种算法求解最优加权系数过程进行迭代的次数,说明ADPM算法收敛速度更快,收敛效果更好。
实施例二:为分析本发明方法的性能优势,设置不同的惩罚因子初值,基于以上参数,比较两种算法的收敛速度。
图7为惩罚因子ρy分别设置为ρy=0.0001,ρy=0.001和ρy=0.01的情况下,两种算法对应的残差Δry的收敛曲线;图8为惩罚因子ρm分别设置为ρm=0.0001,ρm=0.001和ρm=0.01的情况下,两种算法对应的残差Δrm的收敛曲线。从图7和图8中可以看出,设置不同的惩罚因子初值,对ADPM算法的收敛速度性能影响较小,且设置相同的惩罚因子初值情况下,ADPM算法收敛速度更快。
综上所述,本发明所提的多约束下的均匀线性阵列低旁瓣波束形成优化方法,能满足主瓣电平、旁瓣电平和权值幅度约束的要求,同时优化所得的加权系数与现有技术相比具有更大的阵列增益,同时收敛速度和收敛效果好。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内。
Claims (6)
1.一种多约束下的均匀线性阵列低旁瓣波束形成优化方法,其特征在于,包括:
S1、建立均匀线阵接收信号模型;
S2、以阵列增益作为目标函数;步骤S2所述目标函数表达式为:
其中,w表示权向量,am表示导向矢量,()H表示共轭转置,R表示噪声协方差矩阵;
阵列增益的定义为:
其中,E()表示求均值,()H表示共轭转置,||表示复数的模,SNRout和SNRin分别表示输出信噪比和输入信噪比,x(t)和n(t)分别表示输入的有用信号和噪声信号,E(||2)表示信号的平均功率;
S3、以主瓣电平、旁瓣电平和权值幅度为约束,建立优化问题;
S4、采用交替方向惩罚算法最大化阵列增益,求解优化问题得到最优阵列天线加权系数。
2.根据权利要求1所述的一种多约束下的均匀线性阵列低旁瓣波束形成优化方法,其特征在于,步骤S3所述主瓣电平约束表达式为:
Pl≤|wHam|2≤Pu
其中,Pl和Pu分别表示主瓣电平约束的下界和上界;
旁瓣电平约束表达式为:
|wHas|2≤η,s=1,…,S
其中,as表示旁瓣方向导向矢量,as=a(θs),s=1,…,S,θs表示旁瓣方向,S表示将所观察的空域划分为若干离散的角度后,除了波束主瓣所在方向区域外旁瓣区域所有角度的总数,η表示旁瓣电平约束的上界;
权值幅度约束为:
σ≤|wn|≤1,n=1,…,N
其中,σ=1/2M表示权值幅度取值范围的下界,N表示天线阵元的个数,M表示衰减器的比特位数。
4.根据权利要求3所述的一种多约束下的均匀线性阵列低旁瓣波束形成优化方法,其特征在于,步骤S4具体为:引入辅助变量与两个惩罚因子,构建增广拉格朗日函数;将优化问题分解为三个并行的子优化问题,每个子问题中只包含一个约束条件;交替方向惩罚规则用于逐次迭代运算中,直到满足迭代停止条件;得到最优阵列天线加权系数。
5.根据权利要求4所述的一种多约束下的均匀线性阵列低旁瓣波束形成优化方法,其特征在于,引入的两个惩罚因子根据相应的残差进行自适应改变。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010459360.3A CN111551923B (zh) | 2020-05-27 | 2020-05-27 | 一种多约束下的均匀线性阵列低旁瓣波束形成优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010459360.3A CN111551923B (zh) | 2020-05-27 | 2020-05-27 | 一种多约束下的均匀线性阵列低旁瓣波束形成优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111551923A CN111551923A (zh) | 2020-08-18 |
CN111551923B true CN111551923B (zh) | 2022-11-04 |
Family
ID=71999047
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010459360.3A Active CN111551923B (zh) | 2020-05-27 | 2020-05-27 | 一种多约束下的均匀线性阵列低旁瓣波束形成优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111551923B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112099017B (zh) * | 2020-08-29 | 2022-10-14 | 西北工业大学 | 一种基于过扫描的圆环阵低旁瓣波束优化方法 |
CN112214895A (zh) * | 2020-10-14 | 2021-01-12 | 中国船舶重工集团公司第七二四研究所 | 一种变换点数优化设计的傅里叶多波束形成方法 |
CN113094904B (zh) * | 2021-04-12 | 2022-09-16 | 电子科技大学 | 一种非凸多约束下的阵列天线波束赋形优化方法 |
CN113268853B (zh) * | 2021-04-14 | 2022-12-16 | 中国电子科技集团公司电子科学研究院 | 一种天线方向图优化方法、装置及可读存储介质 |
CN113114322B (zh) * | 2021-04-19 | 2022-07-26 | 华南理工大学 | 无线携能通信系统中基于moead的波束赋形方法 |
CN113326471B (zh) * | 2021-05-27 | 2022-12-09 | 哈尔滨工业大学 | 一种星载相控阵部分阵元失效后的方向图重构方法 |
CN113704998B (zh) * | 2021-08-27 | 2023-03-03 | 电子科技大学 | 一种基于交替方向乘子法的波束形成方法 |
CN114166339B (zh) * | 2021-12-07 | 2022-10-18 | 昆明理工大学 | 一种低频结合高频二次波束形成定位方法 |
CN114580159B (zh) * | 2022-02-23 | 2024-09-17 | 中国电子科技集团公司第五十四研究所 | 一种深空测控天线组阵可扩展规模布局优化设计方法 |
CN114996935B (zh) * | 2022-05-31 | 2023-06-06 | 武汉工程大学 | 基于近似admm的波束形成方法、装置、设备及介质 |
CN115276745B (zh) * | 2022-07-12 | 2024-05-10 | 湖北工业大学 | 一种基于交替方向惩罚的鲁棒波束形成方法及系统 |
CN115276746B (zh) * | 2022-07-12 | 2023-05-30 | 湖北工业大学 | 基于交替方向惩罚的频率一致宽带波束形成方法及系统 |
CN115038012B (zh) * | 2022-08-10 | 2022-11-08 | 湖北工业大学 | 基于admm的麦克风阵列鲁棒频率不变波束形成方法 |
CN115833887B (zh) * | 2022-11-08 | 2024-09-20 | 南京理工大学 | 一种动态超表面天线阵列天线选择与波束赋形方法 |
CN117688789B (zh) * | 2024-02-04 | 2024-04-26 | 南京隼眼电子科技有限公司 | 基于admm的相控阵扫描雷达稀疏阵列综合方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103837861A (zh) * | 2014-03-19 | 2014-06-04 | 北京理工大学 | 基于特征子空间的子阵级线性约束自适应波束形成方法 |
CN108337030A (zh) * | 2018-02-07 | 2018-07-27 | 东南大学 | 多天线系统中高功效波束成形方法、装置及设备 |
CN108446437A (zh) * | 2018-02-08 | 2018-08-24 | 电子科技大学 | 一种阵列天线宽波束功率增益优化方法 |
CN110501675A (zh) * | 2019-07-16 | 2019-11-26 | 北京工业大学 | 一种基于mimo雷达低旁瓣发射方向图设计方法 |
CN111062142A (zh) * | 2019-12-30 | 2020-04-24 | 电子科技大学 | 一种基于线性规划的阵列天线宽波束增益优化方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7280070B2 (en) * | 2004-11-30 | 2007-10-09 | Unnikrishna Sreedharan Pillai | Robust optimal shading scheme for adaptive beamforming with missing sensor elements |
-
2020
- 2020-05-27 CN CN202010459360.3A patent/CN111551923B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103837861A (zh) * | 2014-03-19 | 2014-06-04 | 北京理工大学 | 基于特征子空间的子阵级线性约束自适应波束形成方法 |
CN108337030A (zh) * | 2018-02-07 | 2018-07-27 | 东南大学 | 多天线系统中高功效波束成形方法、装置及设备 |
CN108446437A (zh) * | 2018-02-08 | 2018-08-24 | 电子科技大学 | 一种阵列天线宽波束功率增益优化方法 |
CN110501675A (zh) * | 2019-07-16 | 2019-11-26 | 北京工业大学 | 一种基于mimo雷达低旁瓣发射方向图设计方法 |
CN111062142A (zh) * | 2019-12-30 | 2020-04-24 | 电子科技大学 | 一种基于线性规划的阵列天线宽波束增益优化方法 |
Non-Patent Citations (3)
Title |
---|
Beampattern Synthesis With Sidelobe Control and Applications;Ashenafi Yadessa Gemechu,et al;《IEEE TRANSACTIONS ON ANTENNAS AND PROPAGATION》;20200131;第299页 * |
MIMO Radar Transmit-Receive Design for Moving Target Detection in Signal-Dependent Clutter;Xianxiang Yu,et al;《IEEE》;20191231;第2页左侧栏第2段;第5页右侧栏第1段 * |
Power-Gain Pattern Synthesis of Array Antenna With Dynamic Range Ratio Restriction;Shiwen Lei,et al;《IEEE ANTENNAS AND WIRELESS PROPAGATION LETTERS》;20191231;第2691-2692页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111551923A (zh) | 2020-08-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111551923B (zh) | 一种多约束下的均匀线性阵列低旁瓣波束形成优化方法 | |
CN105137399B (zh) | 基于斜投影滤波的雷达自适应波束形成方法 | |
CN103245956B (zh) | 一种基于稳健波束形成算法的gps抗多径方法 | |
CN109254261B (zh) | 基于均匀圆阵epuma的相干信号零陷加深方法 | |
CN102830387A (zh) | 一种基于数据预处理的协方差矩阵正交化波束形成方法 | |
CN107024681A (zh) | 基于杂波知识未确知条件下mimo雷达收发联合优化方法 | |
CN105182302A (zh) | 一种抗快速运动干扰稳健零陷展宽波束形成方法 | |
CN113325385B (zh) | 一种相控阵-mimo雷达模式收发波束形成抗干扰方法 | |
CN108445486A (zh) | 基于协方差矩阵重建和导向矢量修正的波束形成方法 | |
CN110708103B (zh) | 一种无需预延迟的宽带波束形成方法 | |
CN113075622A (zh) | 一种雷达通信一体化中多用户通信的发射波束赋形方法 | |
CN110138430A (zh) | 期望信号doa误差下基于导向矢量估计的稳健宽带波束形成方法 | |
CN104539331B (zh) | 一种基于改进的混合入侵杂草算法阵列天线波束合成方法 | |
CN106125039B (zh) | 基于局域联合处理的改进空时自适应单脉冲测角方法 | |
CN113376584B (zh) | 基于改进对角加载的稳健自适应波束形成方法 | |
CN110208748A (zh) | 基于阵列共轭对称和双迭代算法的雷达波束形成方法 | |
CN104346532B (zh) | 一种mimo雷达降维自适应波束形成方法 | |
CN111817765B (zh) | 一种基于频率约束的广义旁瓣对消宽带波束形成方法 | |
CN110895327B (zh) | 基于直接凸优化建模的鲁棒性自适应波束形成方法 | |
CN111162878B (zh) | 一种基于子阵降维带约束的多域联合抗干扰方法 | |
CN109633600B (zh) | 一种最小冗余线阵mimo-othr的doa估计方法 | |
CN107167803A (zh) | 基于导向矢量失配估计的鲁棒波束域自适应波束形成方法 | |
CN116148777A (zh) | 一种基于阵列虚拟扩展的阵列雷达主瓣干扰抑制方法 | |
CN112347681B (zh) | 一种基于宏基函数阵列互耦特性预测的稳健波束形成方法 | |
Jagtap et al. | Classification and Trends in Adaptive Beamforming Techniques of MIMO Radar |
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 |