Loading [MathJax]/jax/output/SVG/jax.js

基于交替方向惩罚法的低精度量化MIMO雷达恒模波形设计方法

万环 余显祥 全智 廖斌

万环, 余显祥, 全智, 等. 基于交替方向惩罚法的低精度量化MIMO雷达恒模波形设计方法[J]. 雷达学报, 2022, 11(4): 557–569. doi: 10.12000/JR22072
引用本文: 万环, 余显祥, 全智, 等. 基于交替方向惩罚法的低精度量化MIMO雷达恒模波形设计方法[J]. 雷达学报, 2022, 11(4): 557–569. doi: 10.12000/JR22072
WAN Huan, YU Xianxiang, QUAN Zhi, et al. Constant modulus waveform design for low-resolution quantization MIMO radar based on an alternating direction penalty method[J]. Journal of Radars, 2022, 11(4): 557–569. doi: 10.12000/JR22072
Citation: WAN Huan, YU Xianxiang, QUAN Zhi, et al. Constant modulus waveform design for low-resolution quantization MIMO radar based on an alternating direction penalty method[J]. Journal of Radars, 2022, 11(4): 557–569. doi: 10.12000/JR22072

基于交替方向惩罚法的低精度量化MIMO雷达恒模波形设计方法

DOI: 10.12000/JR22072
基金项目: 国家自然科学基金(62171292),广东省自然科学基金(2020A1515010410, 2022A1515010188)
详细信息
    作者简介:

    万 环(1992-),女,江西南昌人,博士。主要研究方向为阵列信号处理、雷达波形设计以及最优化理论算法等

    余显祥(1991-),男,四川人,博士。主要研究方向为雷达波形设计与处理、最优化理论算法以及阵列信号处理等。目前发表IEEE Trans期刊论文10余篇

    全 智(1978-),男,广西柳州人,深圳大学特聘教授,博士生导师,国家优秀青年科学基金获得者。主要研究方向为远距离宽带无线通信系统、射频系统校准与测试、数据驱动的信号处理等。担任IEEE Transactions on Signal Processing等期刊编委

    廖 斌(1983-),男,江西萍乡人,深圳大学特聘研究员,博士生导师。主要研究方向为阵列信号处理、自适应信号处理、雷达信号处理等。担任IEEE Transactions on Aerospace and Electronic Systems等期刊编委

    通讯作者:

    廖斌 binliao@szu.edu.cn

  • 责任主编:梁军利 Corresponding Editor: LIANG Junli
  • 中图分类号: TN958

Constant Modulus Waveform Design for Low-resolution Quantization MIMO Radar Based on an Alternating Direction Penalty Method

Funds: The National Natural Science Foundation of China (62171292), Guangdong Basic and Applied Basic Research Foundation (2020A1515010410, 2022A1515010188)
More Information
  • 摘要: 在MIMO雷达中配备大量有源天线单元可以获得优异的波束形成性能,但会导致系统能耗大、电路复杂及成本高等问题。采用低精度的DAC组件可有效克服上述问题,但现有基于无限精度DAC条件所设计的MIMO雷达波形往往难以直接适用于低精度DAC系统。为此,该文提出了一种离散相位约束下基于最小化积分副主瓣比的低精度量化MIMO雷达恒模波形设计方法。该方法首先采用丁克尔巴赫(Dinkelbach)算法将目标函数二次分数形式转换成减法形式,再利用交替方向惩罚法求解非凸恒模离散相位约束问题。最后通过数值仿真与其他方法进行对比,分析了所提方法的发射方向图与积分副主瓣比性能,验证了该方法的有效性。

     

  • 多输入多输出(Multiple-Input Multiple-Output, MIMO)雷达采用多信号同时发射与接收模式[1-3],可按需调整发射波形,实现高增益窄波束或全区域覆盖的宽波束。相比传统相控阵雷达,MIMO雷达基于波形分集,可为系统发射与接收端提供更多的自由度,提高波达角(Directions of Arrival, DOA)估计精度、目标检测能力、干扰抑制能力[4-6]。MIMO雷达波形设计的一个主要目的是使发射能量集中在目标空间区域(主瓣)范围内,同时降低非目标区域(副瓣)能量辐射,实现波束能量的最佳匹配发射,从而提高雷达目标探测、参数估计和检测性能[7-9]。为了提升性能,MIMO雷达系统通常会配置大量有源天线单元,但采用高精度(超过10位)数模转换器(Digital-to-Analog Converters, DAC)组件会大幅增加系统电路的复杂度、能耗及成本。相反地,低精度DAC组件可以显著降低系统功耗与成本[10-12],已广泛应用于各种场景,例如低功耗超宽带通信系统、大规模或超大规模MIMO系统[13,14]等。

    现有的MIMO雷达波形设计方法主要面向采用高精度DAC组件的系统[15-20],针对采用低精度DAC组件的MIMO雷达波形设计方法相对较少[21-26]。若直接对现有的基于无限精度DAC算法所设计的波形进行量化来适配低精度DAC组件,系统性能将会严重下降。低精度量化波形设计的核心难点在于离散相位或离散幅度的非凸约束求解。最直接的方法是采用穷搜法,但随着信号维度和快拍数的增加,计算复杂度与时间也呈指数量级增加。为了克服计算效率问题,目前主要有两类解决方法,第1类为改进式穷搜法,通过改变搜索策略与范围来减少搜索次数从而降低计算复杂度与运算时间。例如,文献[21]提出了基于块坐标下降(Block Coordinate Descent, BCD)搜索法,使搜索次数从指数级LNtN(L为低精度离散相位或幅度的个数,Nt为发射阵列天线数,N为快拍数)降低至TLNtN(T为该算法收敛时最大迭代次数)。基于此思路,文献[22]提出了两种广义似然上升搜索算法,第t+1次迭代结果可通过似然判断从第t次迭代结果翻转符号得到,进一步降低BCD算法的收敛迭代次数,但该方案仅对极低精度(1比特)波形有较好的求解效果。第2类将离散约束近似为连续函数进行求解,再将无限量化精度解重新量化或映射为低精度。文献[23]采用最小化积分副主瓣比(Integrated Sidelobe-to-Mainlobe Ratio, ISMR)准则,结合特殊的矩阵块结构,利用半正定松弛(Semidefinite Relaxation, SDR)技术,获得高精度发射信号,再通过量化器得到1比特信号。文献[24]采用连续可微函数来逼近1比特信号,并利用交替方向乘子法(Alternation Direction Method of Multipliers, ADMM)求得高精度解,再映射为1比特信号。基于文献[24],文献[25,26]利用辅助变量与发射信号向量相互约束关系,将1比特信号的离散约束转换成区间连续约束,再通过ADMM算法使得信号不断逼近1比特后再映射为1比特,从而提高信号对1比特DAC的适配能力。

    文献[23-26]考虑的是极低精度1比特波形,波形序列元素的实部与虚部的取值均为±E/2NtN(E为发射波形总功率),对应在极坐标轴上体现为有限个相位点(离散相位),即{π/4,3π/4,5π/4,7π/4}。现有的低精度量化波形设计方法主要针对1比特,面向低精度多比特(2~5比特)的波形设计方法相对匮乏。另外,在雷达应用中,为了最大化发射机效率,防止发射机功放非线性失真,通常要求发射波形具有恒模特性,即每个码元的模值是恒定的[27-30]。因此,本文研究基于ISMR最小化准则的低精度(包括极低精度1比特)有限相位恒模波形设计方法。所建模的波形优化问题包含二次分式目标函数和非凸离散约束,难以求解。为此,本文提出了一种基于丁克尔巴赫交替方向惩罚法(Dinkelbach Alternating Direction Penalty Method, DADPM)的优化求解方法。该方法首先利用Dinkelbach算法[31]将目标函数二次分数形式转换成减法形式,再通过ADPM框架求解非凸离散相位约束问题,实现了低精度量化恒模发射波形设计问题的高效求解。

    图1所示,考虑一个共置窄带MIMO雷达,包含Nt个发射天线,按照半波长间距均匀排布。每个通道配置两个B比特DAC组件,分别对基带信号sm(n), m=1,2,,Ntn=1,2,,N的实部与虚部进行量化,第m个天线n时刻的发射形式表示为sB,m(n)=QB(sm(n)),其中QB()表示B比特量化器。第n个快拍的发射波形表示为˜sB,n=[sB,1(n),sB,2(n),,sB,Nt(n)]T。为方便设计与分析,将N个采样快拍下的发射波形序列堆叠为一个向量,即sB=[˜sTB,1,˜sTB,2,,˜sTB,N]T。则B比特恒模发射波形序列第i个元素可以表示为

    图  1  配置低精度DAC组件的MIMO雷达发射端系统结构图
    Figure  1.  System structure diagram of MIMO radar transmitter with low-resolution DACs
    sB(i)=qejφi,i=1,2,,NtN (1)

    式中,q=E/(NtN)为发射波形元素的模值,其中E为发射波形总功率,φi表示s中第i个元素的相位。

    B比特恒模波形元素的相位个数共有L=2B+1个,相位集合XB表示为

    XB{πL+2πLl|l=0,1,,L1} (2)

    B比特恒模发射波形有L个相位均匀分布在以半径为q的极坐标圆周上。图2展示了1比特与2比特波形元素可行域。

    图  2  1比特与2比特波形元素可行域(红点)
    Figure  2.  Feasible areas of 1 bit and 2 bit waveform entries (red dots)

    基于上述系统模型,在远场方向θ处,N个采样快拍下的信号可表示为[16]

    xt(θ)=(INaTt(θ))sB (3)

    式中,代表Kroneker积,()T代表转置。at(θ)Nt维发射阵列导向矢量,表示为

    at(θ)=[1,ej2πdsinθ/λ,,ej2πd(Nt1)sinθ/λ]T (4)

    其中,λ为波长,d为天线阵元间距,且d=λ/2。根据式(3),雷达远场空间功率谱可表示为[7]

    Pt(θ)=xHt(θ)xt(θ)=sHB(INat(θ)aTt(θ))sB =sHBR(θ)sB (5)

    式中,()H代表共轭转置,()*代表共轭,矩阵R(θ)

    R(θ)=INat(θ)aTt(θ) (6)

    当MIMO雷达发射相干波形,式(5)表示相控阵方向图。如果发射相互正交的波形,表明各个方向辐射的功率相等,实现了空域均匀覆盖。若发射相关波形,发射波束方向图取决于波形具体形式。

    定义主瓣区域表示为Θm,副瓣区域为Θs,则B比特发射信号的ISMR表示为[16,20]

    ISMR(sB)=ΘssHBR(θ)sBdθΘmsHBR(θ)sBdθ=sHBΩssBsHBΩmsB (7)

    式中,Ωs=ΘsR(θ)dθ, Ωm=ΘmR(θ)dθ

    本文采用ISMR最小化准则来设计B比特恒模发射信号,因此该优化问题模型描述为

    minsB ISMR(sB)  s.t. angle(sB(i))XB,|sB(i)|=q, i=1,2,,NtN (8)

    式(8)中, angle()表示输入变量的相位。上述问题的目标函数为二次分式,约束条件包含非凸离散相位约束,该问题为非确定性多项式-难(Nondeterministic Polynomial-hard, NP-hard),难以求解。

    本节将提出一种基于ADPM的优化算法对问题(8)进行求解。该方法首先通过Dinkelbach算法将目标函数二次分数形式转换成减法形式,再基于ADPM框架,引入辅助变量,将离散相位约束转换为NtN个独立并行的三角函数问题,通过迭代逐步逼近最优解。

    基于Dinkelbach算法原理,可以将目标函数二次分数形式转换成减法形式,即

    f(sB,ξ)=sHBΞsB (9)

    式中,Ξ表示为

    Ξ=ΩsξΩm (10)

    式中,参数ξ0,在Dinkelbach方法中通过式(11)不断更新:

    ξ(k+1)=ISMR(s(k)B) (11)

    式中,k为迭代次数。这里需注意的是,Ξ可能不是正定矩阵。当Ξ为非正定矩阵时,对该矩阵进行对角加载,使其满足正定,即

    Ξ={Ξ,ΞΞ+μINtN, (12)

    式中,μ>λmin(Ξ)λmin()表示矩阵最小特征值。当Ξ非正定时,目标函数sHB(Ξ+μINtN)sB=sHBΞsB+μNtN,其中μNtN为常数,因此,对Ξ矩阵的对角进行加载操作,不会影响问题(8)的解。在Dinkelbach方法[31]中已证明,问题(8)的目标函数的解等价于f(s(k+1)B,ξ(k+1))=0的解。因此,基于每次迭代更新的ξ(k),可以通过以下问题对sB进行优化:

    minsB f(sB,ξ(k))  s.t. angle(sB(i))XB,|sB(i)|=q, i=1,2,,NtN (13)

    问题(13)中约束条件包括离散相位约束以及恒模约束,可以利用ADPM算法[32]进行求解。

    ADPM算法与利用固定惩罚因子的传统ADMM算法不同,ADPM算法采用动态更新惩罚因子的方式,使惩罚项(˜sBsB22)趋近0,在确保算法的收敛性的同时,可以使得算法能够找到相对较优的可行解。基于该算法框架,引入一个辅助变量˜sB,问题(13)等效表示为

    minsB,˜sB f(sB,ξ(k)) s.tsB=˜sB,angle(˜sB(i))XB,|˜sB(i)|=q, i=1,2,,NtN (14)

    根据ADPM算法原理,式(14)的增广拉格朗日函数表达式为

    L(˜sB,sB,ϱ,p)=f(sB,ξ(k))+{pH(˜sBsB)}+ϱ/2˜sBsB22 (15)

    式中,{}代表实部,pCNtNϱ>0分别为拉格朗日乘子向量与惩罚因子。

    基于ADPM算法框架,在离散相位约束与恒模约束条件下,通过最小化L(˜sB,sB,ϱ,p)准则交替优化更新变量˜sBsB,以及更新变量ϱp。在本文中,用˜s(t)B,s(t)B,ϱ(t),p(t)表示第t次迭代时˜sB,sB,ϱ,p的优化结果。通过ADPM算法求解问题(10)的逼近解,主要有以下优化步骤:

    步骤1 优化变量˜sB

    固定s(t)B,ϱ(t),p(t)的值,最小化增广拉格朗日函数L(˜sB,s(t)B,ϱ(t),p(t)),针对˜sB的优化问题可表示为

    ˜s(t+1)B min˜sBL(˜sB,s(t)B,p(t),ϱ(t))s.t. angle(˜sB(i))XB,|˜sB(i)|=q, i=1,2,,NtN (16)

    忽略目标函数L(˜sB,s(t)B,ϱ(t),p(t))˜sB无关的项,可进行以下变换与化简:

    argmin˜sBL(˜sB,s(t)B,p(t),ϱ(t))=argmin˜sB{p(t)H(˜sBs(t)B)}+ϱ(t)2˜sBs(t)B22=argmin˜sB{(p(t))H˜sB}ϱ(t){s(t)HB˜sB}=argmin˜sB{(p(t)ϱ(t)s(t)B)H˜sB} (17)

    因此,问题(16)可等价表示为

    min˜sB {(p(t)ϱ(t)s(t)B)H˜sB}  s.t. angle(˜sB(i))XB,|˜sB(i)|=q, i=1,2,,NtN (18)

    由问题(18)不难发现,目标函数和约束条件对于{˜sB(i)}是可分离的,即问题(18)中的NtN个变量可采用并行优化。针对˜sB中的第i个元素,可以得到:

    max˜sB(i){(ψ(t)(i))H˜sB(i)}  s.t. angle(˜sB(i))XB,|˜sB(i)|=q,i=1,2,,NtN (19)

    式中,ψ(t)=ϱ(t)s(t)Bp(t)。问题(19)可进一步化简为

    maxωicos(ωiαi)  s.tωiXB (20)

    式中,ωiαi分别为˜sB(i)ψ(t)(i)的相位。问题(20)离散相位闭式解ωi

    ωi=πL+2πL˜l (21)

    式中,

    ˜l=argmaxl=0,1,,L1{τ0,τ1,,τL1} (22)

    其中,τl=cos(π/L+(2π/L)lαi), l=0,1,,L1。通过式(21)得到相位ωi后,代入式(1),可得到˜s(t+1)B 

    步骤2 优化变量sB

    固定˜s(t+1)B,ϱ(t),p(t)的值,最小化增广拉格朗日函数L(˜s(t+1)B,sB,p(t),ϱ(t)),针对sB的更新优化问题可表示为

    s(t+1)BminsBL(˜s(t+1)B,sB,p(t),ϱ(t)) (23)

    与式(17)类似,将问题(23)目标函数L(˜s(t+1)B,sB,p(t),ϱ(t))中与sB无关项去除,并进行向量与矩阵等式变换与化简,可得

    minsBsHBΞ(t+1)sB{(˜ψ(t+1))HsB} (24)

    式中,˜ψ(t+1)=ϱ(t)˜s(t+1)B +p(t), Ξ(t+1)=Ξ+ϱ(t)INtN/2。令式(24)目标函数1阶导数为0可求得最优解为

    s(t+1)B=(Ξ(t+1))1˜ψ(t+1)/2 (25)

    为避免在每次迭代中直接对Ξ(t+1)进行求逆,从而减少计算量,设Ξ的特征分解如下:

    Ξ=UΛUH (26)

    式中,ΛΞ的特征值对角矩阵,U为特征向量酉矩阵,可知Ξ(t+1)的逆可以表示为

    (Ξ(t+1))1=U(Ξ+ϱ(t)INtN/2)1UH (27)

    由于ΛU在更新过程中保持不变,因此可在ADPM算法开始前预先获得,后续更新过程中只需进行对角矩阵Ξ+ϱ(t)INtN/2求逆运算和矩阵相乘运算。因此,设迭代次数为tmax,计算量可由O(tmaxN3tN3)降为O(tmaxN2tN2+N3tN3)

    步骤3 更新惩罚因子ϱ与拉格朗日乘子向量p

    本文采用的ADPM算法基于原始残差值Δˉs(t)=˜s(t+1)Bs(t)B2对惩罚因子ϱ进行更新。若Δˉs(t)未随着迭代次数的增加而减小,则加大ϱ(t),迫使惩罚项趋近0,从而寻找到可行解。如果Δˉs(t)随着迭代次数的增加而减小,则ϱ(t)可以保持不变。即ϱ(t+1)表示为[32]

    ϱ(t+1)={ϱ(t),Δˉs(t+1)Δˉs(t)ϱ(t)δ, (28)

    式中,δ为正实数,且满足δ>1

    固定˜s(t+1)B,s(t + 1)B,ϱ(t+1),拉格朗日乘子向量p的更新值表达式为[32]

    p(t+1)={˜p(t+1),p(t+1)maxν˜p(t+1)/p(t+1)max, (29)

    式中,ν为足够大的正数,˜p(t+1)p(t+1)max表达式为

    ˜p(t+1)=p(t)+ϱ(t+1)(˜s(t+1)B s(t+1)B) (30)
    p(t+1)max=max{|˜p(t+1)(1)|,|˜p(t+1)(2)|,,|˜p(t+1)(NtN)|} (31)

    由式(30)可发现,˜p(t+1)的模值会随着ϱ(t+1)的增大而增大,但是当˜p(t+1)模值过大时,容易导致ADPM算法不收敛。因此,式(29)在更新p(t+1)时,首先判断p(t+1)中元素绝对值的最大值是否超过设定的门限值,若大于该门限值,则进行归一化处理。

    如果ϱ(t+1)p(t+1)的更新值分别为ϱ(t+1)=ϱ(t), p(t+1)=p(t)+ϱ(t+1)(˜s(t+1)B s(t+1)B),则该ADPM算法与传统ADMM算法等价。相比固定惩罚因子的ADMM算法,ADPM算法的惩罚因子动态更新策略可以使得算法在确保收敛的同时有望找到相对更优的可行解。

    根据文献[33],可将迭代终止条件设置为

    {(a): ˜s(t+1)Bs(t+1)B2δ(t+1)pri (b): ϱ(0)˜s(t)B˜s(t+1)B2δ(t+1)dual (32)

    式中,δ(t+1)pri>0δ(t+1)dual>0分别为第t+1次迭代的原始残差与对偶残差的可行性容忍度,这两个值可以根据绝对和相对标准进行选择,即

    δ(t+1)pri=2NtNδabs+δrelmax{s(t+1)B2,˜s(t+1)B2} (33)
    δ(t+1)dual=2NtNδabs+δrel˜s(t+1)B 2 (34)

    其中,δabs>0是绝对误差,δrel>0是相对误差。当满足式(32)中任一个终止条件,ADPM算法迭代停止,得到ˆsB=s(t+1)B。进入到外循环,令s(k+1)B=ˆsB并根据式(11)获得ξ(k+1)。当f(s(k+1)B,ξ(k+1))=0时,输出问题的解sB

    另外,初始值ϱ(0)会影响算法的收敛速度。文献[34]利用正则最小化和二次规划约束找到ADMM迭代收敛因子最小的最优初始值参数。本文采用的ADPM算法为ADMM算法的改进算法,对ADMM初始值的选取方式在ADPM算法中同样适用。因此,本文所根据文献[34]所提方法思想,将ADPM算法的惩罚因子初始值ϱ(0)设置为

    ϱ(0)=λmin(Ξ)λmax(Ξ) (35)

    式中,λmax()表示矩阵最大特征值。

    初始发射波形信号采用正交线性调频信号S(0)S(0)矩阵的第(m,n)个元素表示为

    S(0)(m,n)=ej2πm(n1)/Nejπ(n1)2/NNtN (36)

    式中,m=1,2,,Ntn=1,2,,NNtN×1维初始发射信号向量形式s(0)B可通过堆叠S(0)的列来获得。

    本文所提算法外层循环采用Dinkelbach法,内层循环采用ADPM算法,算法伪代码如表1所示。为分析算法的收敛性,首先证明序列{ξ(k)}是单调减小的。

    表  1  丁克尔巴赫交替方向惩罚法的低精度量化MIMO雷达恒模波形设计算法
    Table  1.  MIMO radar constant modulus waveform design algorithm with low-precision quantized based on DADPM
     输入:s(0)B, ξ(0), B, δ, ν, ϵ
     输出:sB
     步骤1:设置 k=0
     步骤2:初始化:ϱ(0), p(0)
     步骤3:计算Ξ(k)=Ωsξ(k)Ωm
     步骤4:设置 t=0
     步骤5:更新˜s(t+1)Bs(t+1)B,分别通过解问题(16)与问题(23);
     步骤6:更新ϱ(t+1)p(t+1),通过式(28)与式(29);
     步骤7:更新内循环迭代次数,令t=t+1
     步骤8:重复步骤5—步骤7,直到满足式(32)中任一停止条件,存
         储s(t+1)B
     步骤9:令s(k+1)B=s(t+1)B,计算ξ(k+1)=ISMR(s(k+1)B)
     步骤10:更新外循环迭代次数,令k=k+1
     步骤11:重复步骤2—步骤10,直到f(s(k+1)B,ξ(k+1))ϵ
     步骤12:返回 :问题(8)的解sB=s(k+1)B
    下载: 导出CSV 
    | 显示表格

    定义:h(ξ)=minsBf(sB,ξ)=sHBΩssBξsHBΩmsB。由ξ(k+1)=ISMR(s(k)B) = s(k)HBΩss(k)Bs(k)HBΩms(k)B,可推导得

    h(ξ(k))=s(k)HBΩss(k)Bξ(k)(s(k)B)HΩms(k)B=s(k)HBΩms(k)B(ξ(k+1)ξ(k)) (37)

    又因s(k)B=argminsBsHBΩssBξ(k)sHBΩmsB,可得

    h(ξ(k))=s(k)HBΩss(k)Bξ(k)s(k)HBΩms(k)Bs(k1)HBΩss(k1)Bξ(k)s(k1)HBΩms(k1)B=0 (38)

    结合式(37)与式(38),可得s(k)HBΩms(k)B(ξ(k+1)ξ(k))0,即

    ξ(k+1)ξ(k) (39)

    由式(39)得出,Dinkelbach法具有严格的单调性。内层循环采用ADPM算法,其中惩罚因子基于原始残差值动态更新,避免了传统ADMM算法在处理NP-hard问题时依赖惩罚因子初始值选取而存在不收敛问题,保证任意初始值情况下的收敛性[32]。综上,本文提出的DADPM算法中外循环Dinkelbach法具有严格的单调性,内循环ADPM算法具有强收敛性,可得出DADPM算法具有良好的收敛性。

    本文所提出的低精度恒模发射波形设计方法计算复杂度主要与迭代次数、离散相位个数、发射阵列天线数和信号快拍数有关。B比特信号离散相位符号表个数为L=2B+1个,Nt个发射天线,N个快拍下,内循环ADPM算法迭代的计算复杂度为O(tmaxN2tN2+N3tN3)。假设外循环Dinkelbach算法迭代收敛时,最大迭代次数为kmax,则整个算法的计算复杂度为O(tmaxkmaxN2tN2+kmaxN3tN3)

    在本文实验中,测试了不同参数下所提方法的性能。发射阵列天线数设置为Nt=15,样本数为N=100。发射信号总功率固定为E=1。以1°为采样间隔,在整个空域(θ[90,90])均匀采样。主瓣对称情况下,单主瓣区域设置为Ωm=[10,10],副瓣为Ωs=[90,9][11,90]。双主瓣区域设置为Ωm=[40,30][30,40],副瓣Ωs=[90,41][29,29][41,90]。主瓣非对称情况下,双主瓣区域设置为Ωm=[40,30][15,30],副瓣为Ωs=[90,41][29,14][31,90]

    为了方便算法性能分析,针对极低精度(1比特)量化的波形分析,本文提出的基于ADPM算法的1比特量化的波形(DADPM-1bit)对比了基于无穷比特(无量化/无相位约束)的ADMM优化算法设计的恒模发射波形[20](ADMM-bit)与该无穷比特波形直接运用符号函数量化后得到的1比特量化的波形(QADMM-1bit)。还对比了5种针对1比特量化DAC的设计方法,分别为:基于BCD算法的设计方法[21](BCD-1bit)、两种基于广义似然上升搜索算法的设计方法[22](GLAS1-1bit, GLAS2-1bit)、基于SDR算法的设计方法[23](SDR-1bit)、基于ADMM算法的1比特波形设计方法[25](ADMM1-1bit)。

    针对低精度(2~5比特)量化波形,测试了主瓣对称与非对称情况下本文提出的基于DADPM算法的2~5比特量化的波形(DADPM-Bbit,B=2,3,,5)性能,同时与基于ADMM优化算法的无穷比特恒模发射波形直接量化为2~5比特的波形(QADMM-Bbit, B=2,3,,5)进行分析比较。

    本节测试分析极低精度1比特量化波形的性能。对于1比特量化的波形,其B=1,相位符号数为L=2B+1=4个,根据式(2)计算,相位符号表示为{π/4,3π/4,5π/4,7π/4}图3为极低精度量化的对称单主瓣波形序列相位分布图,展示了不同方法的相位分布情况。从图3可见,ADMM-bit算法波形序列元素的相位个数远远超过4个相位,其他本文测试的所有1比特量化的波形方案相位均属于{π/4,3π/4,5π/4,7π/4}中。

    图  3  极低精度1比特量化的对称单主瓣波形序列相位分布图
    Figure  3.  1-bit quantized waveform for single symmetrical mainlobe element phase diagram

    极低精度1比特量化的波形对称单主瓣与双主瓣波形性能分别如图4图5所示。图4(a)图5(a)分别为1比特量化的对称单主瓣与双主瓣波形方向图,横坐标均为空间角度,纵坐标为发射信号在该方向上的平均功率,可通过式(5)计算得到。图4(b)图5(b)分别为1比特量化波形的对称单主瓣与双主瓣ISMR与迭代次数k关系图,横坐标为迭代次数,纵坐标为系统ISMR。

    图  4  极低精度1比特对称单主瓣波形方向图和ISMR与迭代次数关系图
    Figure  4.  1-bit quantized waveform for single symmetrical mainlobe beampattern and the relationship between ISMR versus iteration number
    图  5  极低精度1比特量化的对称双主瓣波形方向图和ISMR与迭代次数关系图
    Figure  5.  1-bit quantized waveform for two symmetrical mainlobe beampattern and the relationship between ISMR versus iteration number

    图4(a)可明显观察到,无穷比特的恒模方案(ADMM-bit)具有最低的副瓣,且主瓣区域较宽、有良好的增益,但该无穷比特波形直接量化为1比特波形后(QADMM-1bit),副瓣明显提高,高于本文测试的所有基于1比特DAC设计的波形,且主瓣中心出现轻微下陷现象。相比另外5种基于1比特DAC设计的波形,本文所提方法具有最低的副瓣,在主瓣区间也有良好的增益。这种波形现象对应在波形ISMR表现为ADMM-bit具有最低ISMR,当该无穷量化精度波形直接量化应用在1比特DAC组件时(QADMM-1bit),ISMR提高了将近11 dB,皆高于其他1比特波形的ISMR。本文提出的DADPM-1bit波形,相较于无穷比特波形,ISMR提高了大约6 dB,相较于其他1比特波形,ISMR值最低。

    图5(a)可观察到,ADMM-bit波形具有最低副瓣,其他1比特量化算法波形方向图大致重合。因此,图5(b)中,ADMM-bit波形ISMR明显最低,其他1比特波形ISMR相差不大,但仍可以看出,无穷比特波形直接量化的QADMM-1bit波形的ISMR最高,本文提出的DADPM-1bit波形的ISMR最低。

    表2为极低精度1比特量化的主瓣对称情况下不同算法1000次蒙特卡罗实验性能统计表。从表格运算时间可发现,相同条件下,SDR-1bit算法运算时间最长,GLAS1-1bit算法运算时间最短。本文所提DADPM-1bit算法性能明显优于其他极低算法性能,但运算时间相对较长。

    表  2  主瓣对称下极低精度量化波形算法性能统计表
    Table  2.  Performance statistics table of the extreme low precision quantized waveform algorithm for symmetrical mainlobe
    主瓣对称情况下方法最小ISMR (dB)最大ISMR (dB)平均ISMR (dB)运算时间(s)
    单主瓣双主瓣单主瓣双主瓣单主瓣双主瓣单主瓣双主瓣
    ADMM-bit–15.7192–7.682511.17829.4329–15.1842–7.30312.35272.0368
    QADMM-1bit–4.1057–3.995310.1032–1.3280–4.2891–3.98722.38652.1003
    ADMM-1bit–7.1965–6.72442.3044–4.6130–5.3700–5.13293.14572.6269
    GLAS1-1bit–7.6930–7.7134–6.9940–6.6514–7.3283–7.14670.01060.0104
    GLAS2-1bit–7.6170–8.3327–6.5995–7.6447–7.0952–7.99700.07360.0748
    SDR-1bit–5.9636–7.0137–5.5988–6.7182–5.7154–6.825637.825038.0630
    BCD-1bit–6.8706–3.0849–2.0389–0.2157–6.8327–3.05870.19720.1964
    DADPM-1bit–8.8275–3.9973–2.5693–1.7831–8.6527–3.98128.51478.3249
    下载: 导出CSV 
    | 显示表格

    本节测试分析了本文所提方法和基于ADMM优化算法的无穷比特(无相位约束)恒模发射波形与其直接量化为不同量化精度的低精度波形性能。测试的量化精度为2比特、3比特、4比特、5比特和无穷比特的ADMM算法。图6(a)图7(a)分别为低精度2~5比特量化的对称单主瓣与双主瓣波形方向图。图6(b)图7(b)分别为低精度2~5比特量化波形的对称单主瓣与双主瓣ISMR与迭代次数k关系图。

    图  6  低精度(2~5比特)量化的对称单主瓣波形方向图和ISMR与迭代次数关系图
    Figure  6.  Low precision quantized waveform for symmetrical single mainlobe beampattern and the relationship between ISMR versus iteration number
    图  7  低精度(2~5比特)量化的对称双主瓣波形方向图和ISMR与迭代次数关系图
    Figure  7.  Low precision quantized waveform for two symmetrical mainlobe beampattern and the relationship between ISMR versus iteration number

    图6(a)可明显观察到,ADMM方案经过低精度量化后,主瓣会出现中心下凹的现象,且本文所提方案明显比低精度量化后的ADMM方案具有更低的副瓣。图6(b)为对称单主瓣ISMR与迭代次数关系图,从图6可发现,相较于ADMM-bit波形,本文所提DADPM-5bit波形ISMR差距约5 dB。与相同精度量化的ADMM波形相比,本文所提低精度量化的DADPM算法设计的波形具有更低的ISMR。

    图7(a)可观察到,本文所提低精度波形与量化后的ADMM波形较为接近,主瓣大体重合,ADMM-bit波形具有较低的副瓣。从图7(b)可观察到,对称双主瓣时,相较于ADMM-bit波形,本文所提DADPM-5bit波形ISMR差距约1 dB。相同量化精度下,本文所提低精度量化的DADPM算法设计的波形比直接量化后的ADMM算法波形具有更低的ISMR。通过图6(b)图7(b)可观察到,随着DAC量化精度提高,波束的ISMR减小,但ISMR下降的幅度值越来越小。

    图8为低精度(2~5比特)非称双主瓣波形方向图和ISMR与迭代次数关系图。从图8(a)可观察到,当两个主瓣不对称时,其中一个主瓣的峰值会下降,甚至主瓣区域内产生零陷。主要由于优化准则为最小化积分(离散累加和)副瓣与主瓣的比值,且矩阵R(θ)具有共轭性,使得Pt(θ)+Pt(θ)=sHB(R(θ)+R(θ))sB=2sHB{R(θ)}sB,因此,非对称情况下无法保证每个主瓣都有一个较高的峰值,而主瓣对称情况下,可使得两个主瓣都具有较好的增益,如图5(a)图7(a)。从图8(b)可观察到,主瓣非对称时,ADMM-bit波形直接量化后,ISMR值提高超过6 dB以上,而本文所提2~5比特低精度DADPM波形与ADMM-bit波形的ISMR相差约1 dB。

    图  8  低精度(2~5比特)量化的非称双主瓣波形方向图和ISMR与迭代次数关系图
    Figure  8.  Low precision quantized waveform for two asymmetrical mainlobe beampattern and the relationship between ISMR versus iteration number

    表3为低精度量化的对称主瓣波形算法在1000次蒙特卡罗实验下的性能统计表。从表3运算时间可发现,本文所提DADPM算法,不同精度对算法的运算时间没有太大影响,精度越高波形ISMR越小。DADPM算法的ISMR明显低于同精度量化的ADMM算法,但DADPM算法运算时间明显长于同精度量化的ADMM算法。主要是因为内循环采用的ADPM算法在迭代的同时动态更新惩罚因子,从而保证算法的收敛性。而ADMM算法惩罚因子直接根据设计人员经验或者实验总结给定,省略了算法寻找合适惩罚因子的过程,算法运算时间更短,但算法性能表现过于依赖惩罚因子。

    表  3  低精度量化的对称主瓣波形算法性能统计表
    Table  3.  Performance statistics table of the low precision algorithm for symmetrical mainlobe
    主瓣对称情况下算法最小ISMR (dB)最大ISMR (dB)平均ISMR (dB)运算时间(s)
    单主瓣双主瓣单主瓣双主瓣单主瓣双主瓣单主瓣双主瓣
    QADMM-2bit–5.5427–4.324511.02377.7345–5.4178–4.07842.38142.1377
    QADMM-3bit–6.3020–4.956011.05327.7081–6.0587–4.65642.34682.1597
    QADMM-4bit–6.5840–5.109111.08637.6597–6.1687–4.97902.15272.3519
    QADMM-5bit–6.6815–5.149011.10277.6038–6.2214–5.00742.25272.4368
    ADMM-bit–15.7192–7.682511.17829.4329–15.1842–7.30312.35272.0368
    DADPM-2bit–9.0532–5.0738–4.8751–4.9024–8.7875–4.70158.77568.2487
    DADPM-3bit–9.5309–5.8123–9.4311–5.7812–9.1178–5.47258.78748.1834
    DADPM-4bit–9.6694–6.2968–9.6103–6.1025–9.2789–6.10348.87438.2981
    DADPM-5bit–10.1160–6.5541–9.9715–6.3251–9.8321–6.27538.89758.3546
    下载: 导出CSV 
    | 显示表格

    本文提出了一种基于低精度量化的MIMO雷达发射波形设计方法。通过设计B比特恒模发射波形序列,使发射波形对低精度量化DAC组件有更好的适配性,实现任意精度波形的最佳匹配发射。为解决所建模的恒模离散相位约束非凸优化问题,首先通过Dinkelbach算法将二次分式转换成减法形式,再运用ADPM算法框架,将离散相位约束转换为并行的三角函数问题,通过交替迭代逐步逼近最优解。低精度DAC组件可降低MIMO雷达电路结构复杂度与能耗,但这些低精度组件也会导致一定程度的性能下降。本文所提出的低精度恒模发射模型设计方法可适用于任意量化精度的发射波形设计,相对其他低精度设计方法,取得了更低的ISMR性能表现,可为实际工程应用中对波形性能要求与DAC量化精度的选择提供理论依据与参考价值。

  • 图  1  配置低精度DAC组件的MIMO雷达发射端系统结构图

    Figure  1.  System structure diagram of MIMO radar transmitter with low-resolution DACs

    图  2  1比特与2比特波形元素可行域(红点)

    Figure  2.  Feasible areas of 1 bit and 2 bit waveform entries (red dots)

    图  3  极低精度1比特量化的对称单主瓣波形序列相位分布图

    Figure  3.  1-bit quantized waveform for single symmetrical mainlobe element phase diagram

    图  4  极低精度1比特对称单主瓣波形方向图和ISMR与迭代次数关系图

    Figure  4.  1-bit quantized waveform for single symmetrical mainlobe beampattern and the relationship between ISMR versus iteration number

    图  5  极低精度1比特量化的对称双主瓣波形方向图和ISMR与迭代次数关系图

    Figure  5.  1-bit quantized waveform for two symmetrical mainlobe beampattern and the relationship between ISMR versus iteration number

    图  6  低精度(2~5比特)量化的对称单主瓣波形方向图和ISMR与迭代次数关系图

    Figure  6.  Low precision quantized waveform for symmetrical single mainlobe beampattern and the relationship between ISMR versus iteration number

    图  7  低精度(2~5比特)量化的对称双主瓣波形方向图和ISMR与迭代次数关系图

    Figure  7.  Low precision quantized waveform for two symmetrical mainlobe beampattern and the relationship between ISMR versus iteration number

    图  8  低精度(2~5比特)量化的非称双主瓣波形方向图和ISMR与迭代次数关系图

    Figure  8.  Low precision quantized waveform for two asymmetrical mainlobe beampattern and the relationship between ISMR versus iteration number

    表  1  丁克尔巴赫交替方向惩罚法的低精度量化MIMO雷达恒模波形设计算法

    Table  1.   MIMO radar constant modulus waveform design algorithm with low-precision quantized based on DADPM

     输入:s(0)B, ξ(0), B, δ, ν, ϵ
     输出:sB
     步骤1:设置 k=0
     步骤2:初始化:ϱ(0), p(0)
     步骤3:计算Ξ(k)=Ωsξ(k)Ωm
     步骤4:设置 t=0
     步骤5:更新˜s(t+1)Bs(t+1)B,分别通过解问题(16)与问题(23);
     步骤6:更新ϱ(t+1)p(t+1),通过式(28)与式(29);
     步骤7:更新内循环迭代次数,令t=t+1
     步骤8:重复步骤5—步骤7,直到满足式(32)中任一停止条件,存
         储s(t+1)B
     步骤9:令s(k+1)B=s(t+1)B,计算ξ(k+1)=ISMR(s(k+1)B)
     步骤10:更新外循环迭代次数,令k=k+1
     步骤11:重复步骤2—步骤10,直到f(s(k+1)B,ξ(k+1))ϵ
     步骤12:返回 :问题(8)的解sB=s(k+1)B
    下载: 导出CSV

    表  2  主瓣对称下极低精度量化波形算法性能统计表

    Table  2.   Performance statistics table of the extreme low precision quantized waveform algorithm for symmetrical mainlobe

    主瓣对称情况下方法最小ISMR (dB)最大ISMR (dB)平均ISMR (dB)运算时间(s)
    单主瓣双主瓣单主瓣双主瓣单主瓣双主瓣单主瓣双主瓣
    ADMM-bit–15.7192–7.682511.17829.4329–15.1842–7.30312.35272.0368
    QADMM-1bit–4.1057–3.995310.1032–1.3280–4.2891–3.98722.38652.1003
    ADMM-1bit–7.1965–6.72442.3044–4.6130–5.3700–5.13293.14572.6269
    GLAS1-1bit–7.6930–7.7134–6.9940–6.6514–7.3283–7.14670.01060.0104
    GLAS2-1bit–7.6170–8.3327–6.5995–7.6447–7.0952–7.99700.07360.0748
    SDR-1bit–5.9636–7.0137–5.5988–6.7182–5.7154–6.825637.825038.0630
    BCD-1bit–6.8706–3.0849–2.0389–0.2157–6.8327–3.05870.19720.1964
    DADPM-1bit–8.8275–3.9973–2.5693–1.7831–8.6527–3.98128.51478.3249
    下载: 导出CSV

    表  3  低精度量化的对称主瓣波形算法性能统计表

    Table  3.   Performance statistics table of the low precision algorithm for symmetrical mainlobe

    主瓣对称情况下算法最小ISMR (dB)最大ISMR (dB)平均ISMR (dB)运算时间(s)
    单主瓣双主瓣单主瓣双主瓣单主瓣双主瓣单主瓣双主瓣
    QADMM-2bit–5.5427–4.324511.02377.7345–5.4178–4.07842.38142.1377
    QADMM-3bit–6.3020–4.956011.05327.7081–6.0587–4.65642.34682.1597
    QADMM-4bit–6.5840–5.109111.08637.6597–6.1687–4.97902.15272.3519
    QADMM-5bit–6.6815–5.149011.10277.6038–6.2214–5.00742.25272.4368
    ADMM-bit–15.7192–7.682511.17829.4329–15.1842–7.30312.35272.0368
    DADPM-2bit–9.0532–5.0738–4.8751–4.9024–8.7875–4.70158.77568.2487
    DADPM-3bit–9.5309–5.8123–9.4311–5.7812–9.1178–5.47258.78748.1834
    DADPM-4bit–9.6694–6.2968–9.6103–6.1025–9.2789–6.10348.87438.2981
    DADPM-5bit–10.1160–6.5541–9.9715–6.3251–9.8321–6.27538.89758.3546
    下载: 导出CSV
  • [1] FORSYTHE K W, BLISS D W, and FAWCETT G S. Multiple-input multiple-output (MIMO) radar: Performance issues[C]. The Thirty-Eighth Asilomar Conference on Signals, Systems and Computers, Pacific Grove, USA, 2004: 310–315.
    [2] XU Haisheng, BLUM R S, WANG Jian, et al. Colocated MIMO radar waveform design for transmit beampattern formation[J]. IEEE Transactions on Aerospace and Electronic Systems, 2015, 51(2): 1558–1568. doi: 10.1109/TAES.2014.140249
    [3] 朱圣棋, 余昆, 许京伟, 等. 波形分集阵列新体制雷达研究进展与展望[J]. 雷达学报, 2021, 10(6): 795–810. doi: 10.12000/JR21188

    ZHU Shengqi, YU Kun, XU Jingwei, et al. Research progress and prospect for the noval waveform diverse array radar[J]. Journal of Radars, 2021, 10(6): 795–810. doi: 10.12000/JR21188
    [4] GODRICH H, HAIMOVICH A M, and BLUM R S. Target localization accuracy gain in MIMO radar-based systems[J]. IEEE Transactions on Information Theory, 2010, 56(6): 2783–2803. doi: 10.1109/TIT.2010.2046246
    [5] STOICA P, HE Hao, and LI Jian. Optimization of the receive filter and transmit sequence for active sensing[J]. IEEE Transactions on Signal Processing, 2012, 60(4): 1730–1740. doi: 10.1109/TSP.2011.2179652
    [6] CAO Siyang and ZHENG Yuanfang. Recent developments in radar waveforms[J]. Journal of Radars, 2014, 3(5): 603–621. doi: 10.3724/SP.J.1300.2014.14044
    [7] GAO Xiang, EDFORS O, RUSEK F, et al. Massive MIMO performance evaluation based on measured propagation data[J]. IEEE Transactions on Wireless Communications, 2015, 14(7): 3899–3911. doi: 10.1109/TWC.2015.2414413
    [8] AHMED S and ALOUINI M S. MIMO radar transmit beampattern design without synthesising the covariance matrix[J]. IEEE Transactions on Signal Processing, 2014, 62(9): 2278–2289. doi: 10.1109/TSP.2014.2310435
    [9] FAN Wen, LIANG Junli, YU Guoyang, et al. MIMO radar waveform design for quasi-equiripple transmit beampattern synthesis via weighted lp -minimization[J]. IEEE Transactions on Signal Processing, 2019, 67(13): 3397–3411. doi: 10.1109/TSP.2019.2917871
    [10] SINGH J, DABEER O, and MADHOW U. On the limits of communication with low-precision analog-to-digital conversion at the receiver[J]. IEEE Transactions on Communications, 2009, 57(12): 3629–3639. doi: 10.1109/TCOMM.2009.12.080559
    [11] LANDAU L and FETTWEIS G. On reconstructable ASK-sequences for receivers employing 1-Bit quantization and oversampling[C]. 2014 IEEE International Conference on Ultra-WideBand (ICUWB), Paris, France, 2014: 180–184.
    [12] RISI C, PERSSON D, and LARSSON E G. Massive MIMO with 1-bit ADC[EB/OL]. http://arxiv.org/abs/1404.7736v1, 2014.
    [13] HOYOS S, SADLER B M, and ARCE G R. Monobit digital receivers for ultrawideband communications[J]. IEEE Transactions on Wireless Communications, 2005, 4(4): 1337–1344. doi: 10.1109/TWC.2005.850270
    [14] MO Jianhua and HEATH R W. Capacity analysis of one-bit quantized MIMO systems with transmitter channel state information[J]. IEEE Transactions on Signal Processing, 2015, 63(20): 5498–5512. doi: 10.1109/TSP.2015.2455527
    [15] STOICA P, LI Jian, and XIE Yao. On probing signal design for MIMO radar[J]. IEEE Transactions on Signal Processing, 2007, 55(8): 4151–4161. doi: 10.1109/TSP.2007.894398
    [16] AUBRY A, DE MAIO A, and HUANG Yongwei. MIMO radar beampattern design via PSL/ISL optimization[J]. IEEE Transactions on Signal Processing, 2016, 64(15): 3955–3967. doi: 10.1109/TSP.2016.2543207
    [17] ZHANG Xiaojun, HE Zishu, RAYMAN-BACCHUS L, et al. MIMO radar transmit beampattern matching design[J]. IEEE Transactions on Signal Processing, 2015, 63(8): 2049–2056. doi: 10.1109/TSP.2015.2398841
    [18] ZHANG Weijian, HU Jinfeng, ZHOU Qihang, et al. Constant modulus waveform design for colocated MIMO radar: A convex relaxation approach[J]. Digital Signal Processing, 2021, 117: 103141. doi: 10.1016/j.dsp.2021.103141
    [19] YU Xianxiang, QIU Hui, YANG Jing, et al. Multispectrally constrained MIMO radar beampattern design via sequential convex approximation[J]. IEEE Transactions on Aerospace and Electronic Systems, 2022, 58(4): 2935–2949. doi: 10.1109/TAES.2022.3150619
    [20] CHENG Ziyang, HAN Chunlin, LIAO Bin, et al. Communication-aware waveform design for MIMO radar with good transmit beampattern[J]. IEEE Transactions on Signal Processing, 2018, 66(21): 5549–5562. doi: 10.1109/TSP.2018.2868042
    [21] HONG Mingyi, RAZAVIYAYN M, LUO Zhiquan, et al. A unified algorithmic framework for block-structured optimization involving big data: With applications in machine learning and signal processing[J]. IEEE Signal Processing Magazine, 2016, 33(1): 57–77. doi: 10.1109/MSP.2015.2481563
    [22] DENG Minglong, CHENG Ziyang, LU Xiaoying, et al. Binary waveform design for MIMO radar with good transmit beampattern performance[J]. Electronics Letters, 2019, 55(19): 1061–1063. doi: 10.1049/el.2019.1602
    [23] WEI Tong, CHU Ping, CHENG Ziyang, et al. Transmit beampattern design for MIMO radar with one-bit DACs via block-sparse SDR[C]. 2020 IEEE 11th Sensor Array and Multichannel Signal Processing Workshop (SAM), Hangzhou, China, 2020: 1–5.
    [24] CHENG Ziyang, LIAO Bin, HE Zishu, et al. Transmit signal design for large-scale MIMO system with 1-bit DACs[J]. IEEE Transactions on Wireless Communications, 2019, 18(9): 4466–4478. doi: 10.1109/TWC.2019.2925343
    [25] WEI Tong, CHENG Ziyang, and LIAO Bin. Transmit beampattern synthesis for MIMO radar with one-bit digital-to-analog converters[J]. Signal Processing, 2021, 188: 108228. doi: 10.1016/j.sigpro.2021.108228
    [26] CHENG Ziyang, SHI Shengnan, HE Zishu, et al. Transmit sequence design for dual-function radar-communication system with one-bit DACs[J]. IEEE Transactions on Wireless Communications, 2021, 20(9): 5846–5860. doi: 10.1109/TWC.2021.3070586
    [27] FAN Wen, LIANG Junli, and LI Jian. Constant modulus MIMO radar waveform design with minimum peak sidelobe transmit beampattern[J]. IEEE Transactions on Signal Processing, 2018, 66(16): 4207–4222. doi: 10.1109/TSP.2018.2847636
    [28] ALDAYEL O, MONGA V, and RANGASWAMY M. Tractable transmit MIMO beampattern design under a constant modulus constraint[J]. IEEE Transactions on Signal Processing, 2017, 65(10): 2588–2599. doi: 10.1109/TSP.2017.2664040
    [29] CHENG Ziyang, HE Zishu, ZHANG Shengmiao, et al. Constant modulus waveform design for MIMO radar transmit beampattern[J]. IEEE Transactions on Signal Processing, 2017, 65(18): 4912–4923. doi: 10.1109/TSP.2017.2718976
    [30] HU Jinfeng, ZHANG Weijian, ZHU Haoming, et al. Constant modulus waveform design for MIMO radar via manifold optimization[J]. Signal Processing, 2022, 190: 108322. doi: 10.1016/j.sigpro.2021.108322
    [31] DINKELBACH W. On nonlinear fractional programming[J]. Management Science, 1967, 13(7): 492–498. doi: 10.1287/mnsc.13.7.492
    [32] YU Xianxiang, CUI Guolong, YANG Jing, et al. Quadratic optimization for unimodular sequence design via an ADPM framework[J]. IEEE Transactions on Signal Processing, 2020, 68: 3619–3634. doi: 10.1109/TSP.2020.2998637
    [33] BOYD S, PARIKH N, CHU E, et al. Distributed optimization and statistical learning via the alternating direction method of multipliers[J]. Foundations and Trends ® in Machine Learning, 2011, 3(1): 1–122. doi: 10.1561/2200000016
    [34] GHADIMI E, TEIXEIRA A, SHAMES I, et al. Optimal parameter selection for the alternating direction method of multipliers (ADMM): Quadratic problems[J]. IEEE Transactions on Automatic Control, 2015, 60(3): 644–658. doi: 10.1109/TAC.2014.2354892
  • 期刊类型引用(2)

    1. 杨诗兴,张国鑫,梁雲飞,易伟,孔令讲. 动平台分布式雷达系统动目标低比特数据检测算法. 雷达学报. 2024(03): 584-600 . 本站查看
    2. 刘德顺,夏德平,陈露,马艳峰. 基于神经网络的FDA-MIMO低截获发射波形和接收波束形成联合设计. 雷达学报. 2024(06): 1239-1251 . 本站查看

    其他类型引用(2)

  • 加载中
图(8) / 表(3)
计量
  • 文章访问数: 1284
  • HTML全文浏览量: 383
  • PDF下载量: 189
  • 被引次数: 4
出版历程
  • 收稿日期:  2022-04-24
  • 修回日期:  2022-08-15
  • 网络出版日期:  2022-08-24
  • 刊出日期:  2022-08-28

目录

/

返回文章
返回