CN107346986A - 一种基于稀疏频控传感器天线阵列的多波束形成方法 - Google Patents
一种基于稀疏频控传感器天线阵列的多波束形成方法 Download PDFInfo
- Publication number
- CN107346986A CN107346986A CN201710566383.2A CN201710566383A CN107346986A CN 107346986 A CN107346986 A CN 107346986A CN 201710566383 A CN201710566383 A CN 201710566383A CN 107346986 A CN107346986 A CN 107346986A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- array
- frequency control
- 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.)
- Granted
Links
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/06—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station
- H04B7/0613—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission
- H04B7/0615—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal
- H04B7/0617—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal for beam forming
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Variable-Direction Aerials And Aerial Arrays (AREA)
Abstract
本发明公开了一种基于稀疏频控传感器天线阵列的多波束形成方法,首先确定初始化频控阵、频率增量及其针对各波束相应的距离维加权向量和角度维加权向量,再根据由合成距离维加权向量和合成角度维加权向量之差构成的联合角度‑距离向量确定稀疏频控阵的阵元位置参数和各阵元频率增量组,进而确定目的频控阵中各阵元位置及相应阵元各频率增量组,最后通过凸优化方法确定稀疏频控阵各放大器的增益参数。本发明由于采用稀疏重建理论的稀疏优化布阵方法,在获得尽可能好的频控阵***性能的条件下节约了阵元数,因而具有可利用更少的天线阵元达到同样低的阵列最高旁瓣高度,以及提高信号处理的速度、降低阵列***的成本等优点。
Description
技术领域
本发明属于频控传感器天线阵列设计技术领域,具体涉及一种基于稀疏频控传感器天线阵列的多波束形成方法的设计。
背景技术
有关无线传播信号的天线阵列信号处理技术已广泛应用于雷达、天文、地震、声纳、医学信号检测、地质探测、无线通信等领域。然而阵列***中阵元数的多少严重地影响到阵列***的成本、控制器和信号处理装置的复杂性以及信号处理的数据量和处理速度。因此,在给定条件下如何设置阵列***中各传感器(阵元天线)的空间位置、频率增量大小及其放大器增益参数以实现最佳的空域滤波和空间谱估计具有重要的意义。优化布阵问题属于经典的高维非线性优化难题,目前针对相位阵普遍采用的布阵方法包括:逐对消去法、半定规划法、模拟退火法、遗传设计法、粒子群优化法等。其中半定规划法、模拟退火法、遗传设计法及粒子群优化设计方法在布阵时均是首先制定一个使旁瓣高度最小的目标函数,通过对目标函数进行处理找到最优的阵列加权或最优的阵元位置,或者同时搜索最优的阵列加权及阵元位置,但由于目标函数不是凸函数,所以涉及多维非线性优化问题的求解或多维参数搜索需进行大量的计算,并且优化的结果往往是需要较多的天线阵元,针对减少所需阵元数的处理复杂。而逐对消去法是对加权和位置分步优化,即在给定孔径条件下,首先设置一个阵元间隔为半波长且权系数相等的阵列,然后去除使旁瓣高度最小的那对阵元,重复该步骤直至剩余的阵元数达到给定值为止,最后利用线性规划方法求出剩余阵元数的最优加权,该方法虽然改善了传统均匀布阵的性能,但它仍不能获得满意的效果。而目前针对频控阵的布阵方法非常少,现有的离散椭球序列方法仍存在利用较多天线阵元的问题。
发明内容
本发明的目的是针对背景技术存在的缺陷,提出了一种基于稀疏频控传感器天线阵列的多波束形成方法,以达到在给定阵列孔径值及最高旁瓣高度不高于给定值的约束下,利用更少的天线阵元达到较低的阵列最高旁瓣高度,同时提高信号处理的速度,降低阵列***的成本。
本发明的技术方案为:一种基于稀疏频控传感器天线阵列的多波束形成方法,包括以下步骤:
S1、确定初始化频控阵及其方位导向矩阵:根据给定的频控阵列孔径,按相邻阵元间隔1/4波长、相应各阵元第一频率增量依次递增7.5KHz,设置一个含有密集阵元的初始化频控阵,并根据设置的阵元数、频率增量组及将观测范围等栅格划分的观测区间数确定该初始化频控阵的方位导向矩阵。
S2、确定频控阵加权矩阵:在给定频控阵列孔径范围内,且所有波束相应的主瓣高度归一化、距离维最高旁瓣高度以及角度维最高旁瓣高度均不大于给定值的约束条件下,通过1-范数合成距离维加权向量、1-范数合成角度维加权向量以及1-范数联合角度-距离向量确定频控阵加权矩阵。其中联合角度-距离向量为合成距离维加权向量和合成角度维加权向量之差。
S3、确定稀疏频控阵阵列:根据联合角度-距离向量中各个大于最大元素0.001倍的元素所对应的位置分别放置天线阵元,同时确定各天线阵元对应的一组频率增量,以此得到一个稀疏频控阵阵列。
S4、设定稀疏频控阵阵列中各放大器的增益参数:采用凸优化方法确定所述稀疏频控阵阵列的最优距离维加权矩阵和最优角度维加权矩阵,并依次由所述最优距离维加权矩阵和最优角度维加权矩阵中各相应列向量的克罗内克积得到各波束对应的最优角度-距离加权向量,同时根据各波束的期望方位信息,分别将各最优角度-距离加权向量中的元素值作为各波束对应天线阵元中所有载频信号所对应放大器的增益参数,得到稀疏频控阵阵列对应的多波束综合图。
本发明的有益效果是:本发明采用基于稀疏重建理论的稀疏优化布阵方法,在获得尽可能好的阵列***性能的条件下节约了阵元数,因而具有可利用更少的天线阵元达到同样低的阵列最高旁瓣高度,以及提高信号处理的速度、降低阵列***的成本等优点。
附图说明
图1所示为本发明实施例提供的采用27个阵元数时的阵元位置示意图。
图2所示为本发明实施例提供的一种基于稀疏频控传感器天线阵列的多波束形成方法流程图。
图3所示为本发明实施例提供的采用27个阵元数时的角度-距离依赖性多波束综合图。
图4所示为本发明实施例提供的离散椭球序列方法采用32个阵元数时的角度-距离依赖性多波束综合图。
具体实施方式
现在将参考附图来详细描述本发明的示例性实施方式。应当理解,附图中示出和描述的实施方式仅仅是示例性的,意在阐释本发明的原理和精神,而并非限制本发明的范围。
本发明实施例以天线阵列的工作波长15cm、在15.5倍波长的孔径范围内设置27个阵元的非均匀线性频控阵为例(阵元位置如图1所示),提供了一种基于稀疏频控传感器天线阵列的多波束形成方法,如图2所示,包括以下步骤S1-S4:
S1、确定初始化频控阵及其方位导向矩阵:根据给定的频控阵列孔径(15cm×15.5=232.5cm),按相邻阵元间隔1/4λ(波长)、阵元数M=125、相应各阵元第一频率增量依次递增7.5KHz,设置一个全向天线均位于一条水平直线上的初始化频控阵,将平面波的来波中心方向θ0分别设定为50°、20°、-40°方向,距离中心r0相应分别定在12km、10km、8km,并设定该水平直线为x轴,各阵元离原点的距离从近到远依次为d1,d2,…,d125,(dm=λ·(m-1)/8,m=1,...,125),在观测范围Ω={6km≤r≤14km,-90°≤θ≤90°}内等栅格间隔地选择14661(81×181)个观测区间,设x轴上最左侧的阵元为参考阵元,且该参考阵元的第一个频率增量为7.5KHz,中心载频为2GHz,则该初始化频控阵的方位导向矩阵为:
A=[a(r1,θ1) a(r2,θ2) … a(rL,θH)] (1)
式中L表示观测范围距离维等间隔划分的观测个数,H表示观测范围角度维等间隔划分的观测个数,(r1,θ1),(r1,θ2),...,(rL,θH)表示由距离r和方向θ形成的各观测方位,a(ri,θj),i=1,...,L;j=1,...,H为方位导向向量。本发明实施例中,将L=81,H=181带入公式(1)得到:
A=[a(r1,θ1) a(r2,θ2) … a(r81,θ181)]
为阵列的方向向量,其中:
式中设阵列左侧的第一个阵元为参考阵元,Δf为参考阵元的第一频率增量,f0为频控阵的中心载频,c为光速,M为阵元数,d为相邻阵元间隔,表示克罗内克积,表示向量aR(ri)中的每一个元素均与向量aΘ(θj)相乘,[·]T表示向量的转置。本发明实施例中,将Δf=7.5KHz,f0=2GHz,d=7.5cm,c=3×108m/s,M=125带入上式,得:
S2、确定频控阵加权矩阵:在给定频控阵列孔径范围内,且所有波束相应的主瓣高度归一化、距离维最高旁瓣高度以及角度维最高旁瓣高度均不大于给定值的约束条件下,通过1-范数合成距离维加权向量、1-范数合成角度维加权向量以及1-范数联合角度-距离向量确定频控阵加权矩阵。其中联合角度-距离向量为合成距离维加权向量和合成角度维加权向量之差。频控阵加权矩阵为下述凸优化过程的结果:
式中N0为给定阵列孔径内虚拟频控阵的天线个数,K为波束个数,和分别为距离维加权矩阵和角度维加权矩阵的第i行,为合成距离维加权向量,为合成角度维加权向量,为联合角度-距离加权向量;和分别为第k个波束方位相应的距离维旁瓣区域和角度维旁瓣区域对应的方位导向矩阵,aR(rk)和aΘ(θk)分别为第k个波束方位相应的距离维导向向量和角度维导向向量,ε和η分别为综合波束图和相应期望波束图{Gd(rk),Gd(θk)}之间的拟合误差,Gd(rk)和Gd(θk)分别为第k个期望波束的距离维波束图和角度维波束图,Wr(:,k)和Wθ(:,k)分别为第k个期望波束相应的距离维加权向量和角度维加权向量;[·]H表示向量的共轭转置;向量的2-范数||·||2等于向量中各元素的平方和再求均方根,表示各向量的2-范数再求和,等价于求向量的1-范数,向量的∞-范数||·||∞等于向量元素的绝对值的最大值。本发明实施例中,取N0=65,ε=η=0.27,K=3,带入公式(2)后得:
Wr(:,k)HaR(rk)=1
Wθ(:,k)HaΘ(θk)=1,k=1,2,...,3
Wr≠0
Wθ≠0
根据以上线性约束优化方法的结果,确定具有稀疏特性的距离维加权矩阵Wr和角度维加权矩阵Wθ,{aR(12000),aΘ(50°)},{aR(10000),aΘ(20°)},{aR(8000),aΘ(-40°)}分别为第1个波束、第2个波束和第3个波束相应的距离维导向向量和角度维导向向量。
和分别为第1个波束的距离维旁瓣区域和角度维旁瓣区域对应的方位导向矩阵。
和分别为第2个波束的距离维旁瓣区域和角度维旁瓣区域对应的方位导向矩阵。
和分别为第3个波束的距离维旁瓣区域和角度维旁瓣区域对应的方位导向矩阵。
S3、确定稀疏频控阵阵列:根据步骤S2中得到的联合角度-距离向量中各个大于最大元素0.001倍的27个元素的空间索引,依次确定频控阵各阵元空间位置,同时确定天线频控阵中各阵元相应的第一频率增量(本发明实施例中即0、37.5、67.5、105、135、172.5、210、240、277.5、307.5、345、382.5、412.5、450、480、517.5、555、585、622.5、652.5、690、757.5、795、825、892.5,单位为KHz),以此得到一个稀疏频控阵阵列。
S4、设定稀疏频控阵阵列中各放大器的增益参数:采用凸优化方法确定所述稀疏频控阵阵列的最优距离维加权矩阵和最优角度维加权矩阵,并依次由所述最优距离维加权矩阵和最优角度维加权矩阵中各相应列向量的克罗内克积得到各波束对应的最优角度-距离加权向量,同时根据各波束的期望方位信息,分别将各最优角度-距离加权向量中的元素值作为各波束对应天线阵元中所有载频信号所对应放大器的增益参数,得到稀疏频控阵阵列对应的多波束综合图。最优距离维加权矩阵和最优角度维加权矩阵为下述凸优化过程的结果:
式中和υ分别为距离维最大旁瓣高度和角度维最大旁瓣高度所能容忍的值,Asparse_r0(:,k)和分别为针对稀疏频控阵阵列第k个波束方位对应的距离维导向向量和角度维导向向量,Asparse_RSL(:,k)和Asparse_ΘSL(:,k)分别为针对稀疏频控阵阵列第k个波束的距离维旁瓣区域和角度维旁瓣区域对应的方位导向矩阵,Wsparse_r(:,k)和Wsparse_θ(:,k)分别为针对稀疏频控阵阵列第k个波束对应的距离维和角度维加权向量。
本发明实施例中,采用凸优化方法所确定的第1个波束指向时的各放大器相应的距离维增益参数分别为:0.0074、0.0170、0.0168、0.0238、0.0442、0.0426、0.0637、0.0854、0.0267、0.0697、0.0927、0.1263、0.1531、0.0786、0.1681、0.2099、0.1308、0.1289、0.1642、0.1074、0.0481、0.0633、0.0432、0.0128、0.0131、0.0106、0.0048;角度维增益参数分别为:0.0041、0.0150、0.0109、0.0165、0.0228、0.0329、0.0401、0.0432、0.0426、0.0445、0.0574、0.0618、0.0457、0.0442、0.0621、0.0546、0.0635、0.0445、0.0570、0.0519、0.0446、0.0389、0.0282、0.0289、0.0148、0.0171、0.0132。
采用凸优化方法所确定的第2个波束指向时的各放大器相应的距离维增益参数分别为:0.0125、0.0034、0.0278、0.0484、0.0017、0.0708、0.0891、0.0044、0.1030、0.0061、0.1151、0.1278、0.0155、0.1378、0.0343、0.1383、0.0458、0.1228、0.0435、0.0916、0.0301、0.0542、0.0106、0.0236、0.0070、0.0033、0.0010;角度维增益参数分别为:0.0039、0.0000、0.0031、0.0019、0.0033、0.0151、0.0266、0.0310、0.0385、0.0415、0.0554、0.0651、0.0591、0.0643、0.0669、0.0659、0.0746、0.0605、0.0670、0.0626、0.0505、0.0461、0.0404、0.0348、0.0156、0.0119、0.0084。
采用凸优化方法所确定的第3个波束指向时的各放大器相应的距离维增益参数分别为:0.0074、0.0170、0.0168、0.0238、0.0442、0.0426、0.0637、0.0854、0.0267、0.0697、0.0927、0.1263、0.1531、0.0786、0.1681、0.2099、0.1308、0.1289、0.1642、0.1074、0.0481、0.0633、0.0432、0.0128、0.0131、0.0106、0.0048;角度维增益参数分别为:0.0139、0.0233、0.0161、0.0256、0.0320、0.0390、0.0477、0.0454、0.0406、0.0441、0.0600、0.0538、0.0373、0.0328、0.0635、0.0449、0.0599、0.0377、0.0514、0.0484、0.0411、0.0411、0.0233、0.0257、0.0164、0.0220、0.0186。
根据各波束的期望方位信息,分别将各最优角度-距离加权向量中的元素值作为各波束对应天线阵元中所有载频信号所对应放大器的增益参数,最后根据各波束相应的距离-角度维加权向量确定角度-距离依赖性多波束综合图。
采用本发明实施例提供的一种基于稀疏频控传感器天线阵列的多波束形成方法,最终得到的角度-距离依赖性多波束综合图如图3所示,而作为对比,根据现有的离散椭球序列方法,采用由32个阵元构成的均匀频控阵所对应的3波束角度-距离依赖性多波束综合图如图4所示,其中各天线阵元所对应的第一个频率增量依次增加30KHz。通过对比图3和图4可知,采用本发明方法所布置的非均匀线性频控阵,与现有的离散椭球序列方法相比,在15.5个波长的孔径范围内,获得的距离维最大旁瓣高度和角度维最大旁瓣高度分别为-38.81dB和-20.25dB,在阵列最高旁瓣高度相近及孔径范围相同的条件下本发明仅需27个阵元,而离散椭球序列方法需32个阵元,获得的距离维最大旁瓣高度和角度维最大旁瓣高度分别为-38.26dB和-20.43dB,本发明所需阵元数减少了18.52%,且综合性能相当。因而本发明与现有技术相比,具有可利用更少的天线阵元达到同样低的阵列最高旁瓣高度,以及提高信号处理的速度、降低阵列***的成本等优点。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (9)
1.一种基于稀疏频控传感器天线阵列的多波束形成方法,其特征在于,包括以下步骤:
S1、根据给定的频控阵列孔径设置初始化频控阵,并确定该初始化频控阵的方位导向矩阵;
S2、根据初始化频控阵及其方位导向矩阵,在约束条件下确定频控阵加权矩阵;
S3、根据频控阵加权矩阵得到稀疏频控阵阵列;
S4、设定稀疏频控阵阵列中各放大器的增益参数,得到稀疏频控阵阵列对应的多波束综合图。
2.根据权利要求1所述的多波束形成方法,其特征在于,所述步骤S1具体为:
根据给定的频控阵列孔径,按相邻阵元间隔1/4波长、相应各阵元第一频率增量依次递增7.5KHz,设置一个含有密集阵元的初始化频控阵,并根据设置的阵元数、频率增量组及将观测范围等栅格划分的观测区间数确定该初始化频控阵的方位导向矩阵。
3.根据权利要求2所述的多波束形成方法,其特征在于,所述方位导向矩阵为:
A=[a(r1,θ1) a(r2,θ2) … a(rL,θH)] (1)
式中L表示观测范围距离维等间隔划分的观测个数,H表示观测范围角度维等间隔划分的观测个数,(r1,θ1),(r1,θ2),...,(rL,θH)表示由距离r和方向θ形成的各观测方位,a(ri,θj),i=1,...,L;j=1,...,H为方位导向向量;其中:
<mrow>
<msub>
<mi>a</mi>
<mi>R</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mi>i</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mn>2</mn>
<msub>
<mi>&pi;r</mi>
<mi>i</mi>
</msub>
<mi>&Delta;</mi>
<mi>f</mi>
<mo>/</mo>
<mi>c</mi>
</mrow>
</msup>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mi>j</mi>
<mn>2</mn>
<msub>
<mi>&pi;r</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mi>&Delta;</mi>
<mi>f</mi>
<mo>/</mo>
<mi>c</mi>
</mrow>
</msup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
</mrow>
<mrow>
<msub>
<mi>a</mi>
<mi>&Theta;</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>j</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msup>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mi>d</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>f</mi>
<mi>d</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mi>c</mi>
</mrow>
</msup>
</mtd>
<mtd>
<mn>...</mn>
</mtd>
<mtd>
<msup>
<mi>e</mi>
<mrow>
<mo>-</mo>
<mi>j</mi>
<mn>2</mn>
<mi>&pi;</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>f</mi>
<mn>0</mn>
</msub>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<mi>M</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mi>&Delta;</mi>
<mi>f</mi>
<mi>d</mi>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mi>&theta;</mi>
<mo>)</mo>
<mo>/</mo>
<mi>c</mi>
</mrow>
</msup>
</mtd>
</mtr>
</mtable>
</mfenced>
<mi>T</mi>
</msup>
</mrow>
式中设阵列左侧的第一个阵元为参考阵元,Δf为参考阵元的第一频率增量,f0为频控阵的中心载频,c为光速,M为阵元数,d为相邻阵元间隔, 表示克罗内克积,表示向量aR(ri)中的每一个元素均与向量aΘ(θj)相乘,[·]T表示向量的转置。
4.根据权利要求1所述的多波束形成方法,其特征在于,所述步骤S2中的约束条件为:
在给定频控阵列孔径范围内,所有波束相应的主瓣高度归一化、距离维最高旁瓣高度以及角度维最高旁瓣高度均不大于给定值。
5.根据权利要求1所述的多波束形成方法,其特征在于,所述步骤S2中的频控阵加权矩阵包括1-范数合成距离维加权向量、1-范数合成角度维加权向量以及1-范数联合角度-距离向量;所述联合角度-距离向量为合成距离维加权向量和合成角度维加权向量之差。
6.根据权利要求5所述的多波束形成方法,其特征在于,所述频控阵加权矩阵为下述凸优化过程的结果:
<mrow>
<mtable>
<mtr>
<mtd>
<mi>min</mi>
</mtd>
<mtd>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
</munderover>
<mo>|</mo>
<mo>|</mo>
<msubsup>
<mi>W</mi>
<mi>r</mi>
<mi>i</mi>
</msubsup>
<mo>-</mo>
<msubsup>
<mi>W</mi>
<mi>&theta;</mi>
<mi>i</mi>
</msubsup>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>2</mn>
</msub>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
</munderover>
<mo>|</mo>
<mo>|</mo>
<msubsup>
<mi>W</mi>
<mi>r</mi>
<mi>i</mi>
</msubsup>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>2</mn>
</msub>
<mo>+</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<msub>
<mi>N</mi>
<mn>0</mn>
</msub>
</munderover>
<mo>|</mo>
<mo>|</mo>
<msubsup>
<mi>W</mi>
<mi>&theta;</mi>
<mi>i</mi>
</msubsup>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>2</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>s</mi>
<mo>.</mo>
<mi>t</mi>
<mo>.</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>W</mi>
<mi>r</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mo>:</mo>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mi>H</mi>
</msup>
<msubsup>
<mi>A</mi>
<mrow>
<mi>R</mi>
<mi>S</mi>
<mi>L</mi>
</mrow>
<mi>k</mi>
</msubsup>
<mo>|</mo>
<msub>
<mo>|</mo>
<mi>&infin;</mi>
</msub>
<mo>&le;</mo>
<mi>&epsiv;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>W</mi>
<mi>&theta;</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mo>:</mo>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mi>H</mi>
</msup>
<msubsup>
<mi>A</mi>
<mrow>
<mi>&Theta;</mi>
<mi>S</mi>
<mi>L</mi>
</mrow>
<mi>k</mi>
</msubsup>
<mo>|</mo>
<msub>
<mo>|</mo>
<mi>&infin;</mi>
</msub>
<mo>&le;</mo>
<mi>&eta;</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>W</mi>
<mi>r</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mo>:</mo>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mi>H</mi>
</msup>
<msub>
<mi>a</mi>
<mi>R</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>W</mi>
<mi>&theta;</mi>
</msub>
<msup>
<mrow>
<mo>(</mo>
<mo>:</mo>
<mo>,</mo>
<mi>k</mi>
<mo>)</mo>
</mrow>
<mi>H</mi>
</msup>
<msub>
<mi>a</mi>
<mi>&Theta;</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>&theta;</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mi>k</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>...</mn>
<mo>,</mo>
<mi>K</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>W</mi>
<mi>r</mi>
</msub>
<mo>&NotEqual;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow></mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>W</mi>
<mi>&theta;</mi>
</msub>
<mo>&NotEqual;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式中N0为给定阵列孔径内虚拟频控阵的天线个数,K为波束个数,和分别为距离维加权矩阵和角度维加权矩阵的第i行,为合成距离维加权向量,为合成角度维加权向量,为联合角度-距离加权向量;和分别为第k个波束方位相应的距离维旁瓣区域和角度维旁瓣区域对应的方位导向矩阵,aR(rk)和aΘ(θk)分别为第k个波束方位相应的距离维导向向量和角度维导向向量,ε和η分别为综合波束图和相应期望波束图{Gd(rk),Gd(θk)}之间的拟合误差,Gd(rk)和Gd(θk)分别为第k个期望波束的距离维波束图和角度维波束图,Wr(:,k)和Wθ(:,k)分别为第k个期望波束相应的距离维加权向量和角度维加权向量;[·]H表示向量的共轭转置;向量的2-范数||·||2等于向量中各元素的平方和再求均方根,表示各向量的2-范数再求和,等价于求向量的1-范数,向量的∞-范数||·||∞等于向量元素的绝对值的最大值。
7.根据权利要求5所述的多波束形成方法,其特征在于,所述步骤S3具体为:
根据所述联合角度-距离向量中各个大于最大元素0.001倍的元素所对应的位置分别放置天线阵元,同时确定各天线阵元对应的一组频率增量,以此得到一个稀疏频控阵阵列。
8.根据权利要求1所述的多波束形成方法,其特征在于,所述步骤S4具体为:
采用凸优化方法确定所述稀疏频控阵阵列的最优距离维加权矩阵和最优角度维加权矩阵,并依次由所述最优距离维加权矩阵和最优角度维加权矩阵中各相应列向量的克罗内克积得到各波束对应的最优角度-距离加权向量,同时根据各波束的期望方位信息,分别将各最优角度-距离加权向量中的元素值作为各波束对应天线阵元中所有载频信号所对应放大器的增益参数,得到稀疏频控阵阵列对应的多波束综合图。
9.根据权利要求8所述的多波束形成方法,其特征在于,所述最优距离维加权矩阵和最优角度维加权矩阵为下述凸优化过程的结果:
式中和υ分别为距离维最大旁瓣高度和角度维最大旁瓣高度所能容忍的值,Asparse_r0(:,k)和分别为针对稀疏频控阵阵列第k个波束方位对应的距离维导向向量和角度维导向向量,Asparse_RSL(:,k)和Asparse_ΘSL(:,k)分别为针对稀疏频控阵阵列第k个波束的距离维旁瓣区域和角度维旁瓣区域对应的方位导向矩阵,Wsparse_r(:,k)和Wsparse_θ(:,k)分别为针对稀疏频控阵阵列第k个波束对应的距离维和角度维加权向量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710566383.2A CN107346986B (zh) | 2017-07-12 | 2017-07-12 | 一种基于稀疏频控传感器天线阵列的多波束形成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710566383.2A CN107346986B (zh) | 2017-07-12 | 2017-07-12 | 一种基于稀疏频控传感器天线阵列的多波束形成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107346986A true CN107346986A (zh) | 2017-11-14 |
CN107346986B CN107346986B (zh) | 2020-11-06 |
Family
ID=60257656
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710566383.2A Expired - Fee Related CN107346986B (zh) | 2017-07-12 | 2017-07-12 | 一种基于稀疏频控传感器天线阵列的多波束形成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107346986B (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108446437A (zh) * | 2018-02-08 | 2018-08-24 | 电子科技大学 | 一种阵列天线宽波束功率增益优化方法 |
CN108761414A (zh) * | 2018-05-25 | 2018-11-06 | 电子科技大学 | 一种基于频控阵的s形干扰波束测试方法 |
CN108896983A (zh) * | 2018-05-10 | 2018-11-27 | 电子科技大学 | 一种基于频控阵的时不变空间聚焦波束形成方法 |
CN109298395A (zh) * | 2018-09-28 | 2019-02-01 | 西安建筑科技大学 | 一种基于最大信干噪比的稀疏阵列波束形成方法 |
CN110069896A (zh) * | 2019-05-29 | 2019-07-30 | 重庆邮电大学 | 基于稀疏2d线性阵列的涡旋电磁波产生与优化方法 |
CN110176953A (zh) * | 2019-05-28 | 2019-08-27 | 电子科技大学 | 一种基于广义特征值分解的平面频控阵波束形成方法 |
CN111175753A (zh) * | 2020-02-24 | 2020-05-19 | 四川九洲电器集团有限责任公司 | 一种车载防撞雷达分波目标探测方法 |
CN111355517A (zh) * | 2020-03-17 | 2020-06-30 | 电子科技大学 | 一种面向高速移动用户的频控阵基站协作传输方法 |
CN113466847A (zh) * | 2020-03-31 | 2021-10-01 | 南京理工大学 | 基于lfmcw的非均匀频偏频控阵的实现方法和装置 |
CN113821846A (zh) * | 2021-08-02 | 2021-12-21 | 电子科技大学 | 基于流形几何的频控阵阵元位置设计方法 |
CN115130331A (zh) * | 2022-08-30 | 2022-09-30 | 湖北工业大学 | 基于稀疏阵列的鲁棒频率不变波束形成方法 |
CN113466847B (zh) * | 2020-03-31 | 2024-08-02 | 南京理工大学 | 基于lfmcw的非均匀频偏频控阵的实现方法和装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105842702A (zh) * | 2016-05-13 | 2016-08-10 | 南京航空航天大学 | 多波束成像声呐的旁瓣抑制方法、阵列稀疏方法 |
-
2017
- 2017-07-12 CN CN201710566383.2A patent/CN107346986B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105842702A (zh) * | 2016-05-13 | 2016-08-10 | 南京航空航天大学 | 多波束成像声呐的旁瓣抑制方法、阵列稀疏方法 |
Non-Patent Citations (2)
Title |
---|
HUI CHEN 等: "Sparse reconstruction based frequency diverse array transmit beampattern synthesis", 《2015 3RD INTERNATIONAL WORKSHOP ON COMPRESSED SENSING THEORY AND ITS APPLICATIONS TO RADAR, SONAR AND REMOTE SENSING (COSERA)》 * |
陈慧: "稀疏欠采样阵列信号处理方法研究", 《中国博士学位论文全文数据库信息科技辑(月刊)》 * |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108446437A (zh) * | 2018-02-08 | 2018-08-24 | 电子科技大学 | 一种阵列天线宽波束功率增益优化方法 |
CN108446437B (zh) * | 2018-02-08 | 2021-04-30 | 电子科技大学 | 一种阵列天线宽波束功率增益优化方法 |
CN108896983A (zh) * | 2018-05-10 | 2018-11-27 | 电子科技大学 | 一种基于频控阵的时不变空间聚焦波束形成方法 |
CN108761414A (zh) * | 2018-05-25 | 2018-11-06 | 电子科技大学 | 一种基于频控阵的s形干扰波束测试方法 |
CN108761414B (zh) * | 2018-05-25 | 2022-05-20 | 电子科技大学 | 一种基于频控阵的s形干扰波束测试方法 |
CN109298395B (zh) * | 2018-09-28 | 2020-12-22 | 西安建筑科技大学 | 一种基于最大信干噪比的稀疏阵列波束形成方法 |
CN109298395A (zh) * | 2018-09-28 | 2019-02-01 | 西安建筑科技大学 | 一种基于最大信干噪比的稀疏阵列波束形成方法 |
CN110176953A (zh) * | 2019-05-28 | 2019-08-27 | 电子科技大学 | 一种基于广义特征值分解的平面频控阵波束形成方法 |
CN110176953B (zh) * | 2019-05-28 | 2020-07-07 | 电子科技大学 | 一种基于广义特征值分解的平面频控阵波束形成方法 |
CN110069896A (zh) * | 2019-05-29 | 2019-07-30 | 重庆邮电大学 | 基于稀疏2d线性阵列的涡旋电磁波产生与优化方法 |
CN111175753B (zh) * | 2020-02-24 | 2021-08-17 | 四川九洲电器集团有限责任公司 | 一种车载防撞雷达分波目标探测方法 |
CN111175753A (zh) * | 2020-02-24 | 2020-05-19 | 四川九洲电器集团有限责任公司 | 一种车载防撞雷达分波目标探测方法 |
CN111355517B (zh) * | 2020-03-17 | 2021-04-06 | 电子科技大学 | 一种面向高速移动用户的频控阵基站协作传输方法 |
CN111355517A (zh) * | 2020-03-17 | 2020-06-30 | 电子科技大学 | 一种面向高速移动用户的频控阵基站协作传输方法 |
CN113466847A (zh) * | 2020-03-31 | 2021-10-01 | 南京理工大学 | 基于lfmcw的非均匀频偏频控阵的实现方法和装置 |
CN113466847B (zh) * | 2020-03-31 | 2024-08-02 | 南京理工大学 | 基于lfmcw的非均匀频偏频控阵的实现方法和装置 |
CN113821846A (zh) * | 2021-08-02 | 2021-12-21 | 电子科技大学 | 基于流形几何的频控阵阵元位置设计方法 |
CN113821846B (zh) * | 2021-08-02 | 2023-04-21 | 电子科技大学 | 基于流形几何的频控阵阵元位置设计方法 |
CN115130331A (zh) * | 2022-08-30 | 2022-09-30 | 湖北工业大学 | 基于稀疏阵列的鲁棒频率不变波束形成方法 |
Also Published As
Publication number | Publication date |
---|---|
CN107346986B (zh) | 2020-11-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107346986A (zh) | 一种基于稀疏频控传感器天线阵列的多波束形成方法 | |
CN103020363B (zh) | 一种通过改善阵列波束方向图旁瓣特性设计天线的方法 | |
CN107026686B (zh) | 一种零陷跟踪干扰源的任意形状波束快速赋形方法 | |
CN104615854A (zh) | 一种基于稀疏约束的波束展宽和旁瓣抑制方法 | |
CN108919199A (zh) | 多波束成像声纳稀疏阵列的旁瓣抑制方法、阵列稀疏方法 | |
CN111062142B (zh) | 一种基于线性规划的阵列天线宽波束增益优化方法 | |
CN101420068A (zh) | 一种传感器天线阵列的布阵方法 | |
CN106650104A (zh) | 考虑互耦效应的宽带非频变稀疏阵列综合方法 | |
CN111400919A (zh) | 阵列天线中的低旁瓣波束设计方法 | |
CN109214088B (zh) | 一种最小间距可控的大规模超稀疏平面阵列快速布局方法 | |
Sartori et al. | Hybrid design of non-regular linear arrays with accurate control of the pattern sidelobes | |
US8742983B2 (en) | Radar reception signal processing apparatus and method thereof | |
CN104539331B (zh) | 一种基于改进的混合入侵杂草算法阵列天线波束合成方法 | |
CN108776347B (zh) | 一种双极化天线阵列基于零陷展宽技术的高动态gnss干扰抑制方法 | |
CN104346532B (zh) | 一种mimo雷达降维自适应波束形成方法 | |
CN113671439A (zh) | 基于非均匀智能超表面阵列的无人机集群测向***及方法 | |
CN111262612B (zh) | 一种基于笔形波束的阵列赋形波束解析合成方法 | |
CN112858994A (zh) | 一种基于均匀圆阵的比幅测向方法 | |
CN107783081A (zh) | 一种针对宽带雷达的sstmv最小方差波束形成方法 | |
CN109490846B (zh) | 基于空时联合优化的多输入多输出雷达波形设计方法 | |
CN116718999A (zh) | 一种数字子阵划分方法、***、设备及介质 | |
CN111988078B (zh) | 一种基于三维台阶阵列的方向-距离自适应波束形成方法 | |
CN110208830B (zh) | 一种基于空时二维稀疏阵列的导航抗干扰方法 | |
CN111211425B (zh) | 一种超大扫描角的不规则子阵排布优化方法 | |
CN111293438B (zh) | 一种基于均匀线阵宽带波束形成的栅瓣抑制方法及*** |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20201106 Termination date: 20210712 |