具体实施方式
参照图1,为本发明的一种基于二阶锥规划的宽带发射数字波束形成干扰设计方法流程图;其中所述基于二阶锥规划的宽带发射数字波束形成干扰设计方法,包括以下步骤:
步骤1,确定均匀线阵和宽带数字信号,所述均匀线阵包括M个阵元,每个阵元分别对应一个阵元通道,对所述宽带数字信号进行分路,得到M路宽带数字信号。
然后基于二阶锥规划的算法计算得到每个阵元通道对应的FIR滤波器冲激响应,记为M个FIR滤波器冲激响应;M为大于1的正整数;M路宽带数字信号与M个FIR滤波器冲激响应一一对应。
通过计算得到每个阵元通道对应的FIR滤波器的冲激响应hm,其中m=1,2,…,M;M>1为均匀线阵阵元的个数。
1.1确定宽带数字信号S,所述宽带数字信号S为均匀线阵发射的数字信号中数字信号带宽与经过上变频处理后的数字信号中心频率之间的比值在7%-15%之间的数字信号;将所述宽带数字信号S作为待发射干扰信号,其中宽带数字信号S的长度为Γ,Γ≥1,宽带数字信号S的带宽范围为[f
l,f
u],f
l表示宽带数字信号S的最低频率,f
u表示宽带数字信号S的最高频率;则宽带数字信号S的角频率范围为[ω
l,ω
u],ω
l表示宽带数字信号的角频率最小值,ω
u表示宽带数字信号的角频率最大值;将宽带数字信号的角频率范围离散化为N
f个频率栅格点,得到N
f个频率栅格点的角频率
其满足:
其中,i=1,2,…,Nf,ωl=2πfl,ωu=2πfu。
所述宽带数字信号发射到空中的表现形式为波束,将波束主瓣指向与M个阵元连线之间的夹角,记为期望形成波束主瓣指向φd,0°<φd<180°;如图2所示,以M个阵元连线为0°方向、均匀线阵的法线方向为90°方向建立坐标系来确定期望形成波束主瓣指向φd的大小。
1.2根据图2,宽带数字信号S经过分路器(即根据功率对宽带数字信号S进行分路),得到M路宽带数字信号,M路宽带数字信号功率之和与宽带数字信号S功率相等,且每路宽带数字信号的长度均为Γ;现对如何求解每个阵元通道对应的FIR滤波器冲激响应进行推导。
1.2.1将期望形成波束主瓣指向φd上的宽带数字信号表示为:
其中,ω=2πf',f'为频率变量,f'∈[f
l,f
u];τ
m(φ
d)表示第m个阵元的宽带数字信号延迟,第m个阵元的宽带数字信号延迟τ
m(φ
d)由均匀线阵布局决定,在均匀线阵中,
M表示均匀线阵包括的阵元总个数;c表示光速,
表示阵元间距,h
m[k]表示第m个阵元通道对应的FIR滤波器冲激响应中第k个FIR滤波器系数,k=0,1,…,L-1,L表示每个阵元通道对应的FIR滤波器冲激响应包括的FIR滤波器系数总个数,L≥1;s[n]表示将宽带数字信号分路为M路后每路宽带数字信号中的第n个元素值,s[n-k]表示对s[n]延迟长度k后的元素值,n=1,2,…,Γ,Γ表示将宽带数字信号分路为M路后每路宽带数字信号包括的元素值总个数,上标*表示共轭。
1.2.2将期望形成波束主瓣指向φd上的宽带数字信号频谱表示为Y(ejω):
其中,n-k=0,1,2,…,Γ-1,k=0,1,…,L-1,m=1,2,…,M,
令hm=[hm[0],…,hm[L-1]],m=1,…,M,则得到M个FIR滤波器冲激响应h,M个FIR滤波器冲激响应h是ML×1维向量,h=[h1,…,hM]T,上标T表示转置。
然后得到M×1维向量x(ejω),x(ejω)=[X1(ejω),...,Xm(ejω),...,XM(ejω)]T,那么波束形成器(即图2所示结构)的输出信号频谱Y(ejω)为:
其中,
I表示M×M维单位矩阵,
代表张量积运算,L×1维向量e
T(e
jω)的定义如下:
e(ejω)=[1,e-jω,e-j2ω,…,e-j(L-1)ω]T
1.2.2波束形成器的输出信号功率谱密度Sy(ω)如下式所示:
其中,S
xx(ω)=E{x(e
jω)x
H(e
jω)}表示待发射干扰信号的自相关矩阵,维度为M×M;上标H表示共轭转置,上标T表示转置,上标*表示共轭,
表示张量积。
然后计算得到波束形成器最终的功率谱密度
表达式如下:
那么,在频带ω∈[ωl,ωu]内,宽带波束形成器的输出信号功率Py为:
1.2.3对波束形成器最常见的优化方法是在满足期望方向上信号增益的同时使得其输出总功率达到最小;此外根据不同情况,还会附加一些约束条件,例如旁瓣衰减等。
为了解决这些优化问题,须使用待发射干扰信号的自相关矩阵,待发射干扰信号的自相关矩阵Sxx(ω)可以进一步通过下式求得:
从上式可以看出,待发射干扰信号的自相关矩阵S
xx(ω)由期望形成波束主瓣指向φ
d发射的宽带数字信号频谱自相关矩阵S
d(ω),以及与干扰信号不相关的传感器噪声的自相关矩阵
组成;其中,
σ
n表示噪声的标准差,取0.01即可;I表示M×M维单位矩阵,S
d(ω)表示宽带数字信号的频谱幅度,即
令S
d(ω)≡1,得到经过波束形成器后的期望宽带数字信号频谱
其中:
为了将有限的干扰功率尽可能地集中在主瓣方向,从而保证干扰的总体效益,可以采用基于线性约束最小功率准则(LCMP)对波束形成器进行优化,即在满足期望形成干扰波束的方位角上为单位增益的同时,对特定方向附加线性约束,宽带波束形成器的输出信号功率Py的值达到最小;进而得到基于线性约束最小功率准则(LCMP)的波束形成优化问题为:
其中,subject to表示约束条件,C为约束矩阵,f为约束值矢量,约束矩阵和约束值矢量在此看作未知变量;约束条件CTh=f保证了在期望方向上的无失真响应;当宽带数字信号通过波束形成器到达M个阵元并在期望形成波束主瓣指向φd形成波束时,假设其没有产生畸变,则得到的期望宽带数字信号频谱为Yd(ejω):
将上式拆分成实部和虚部的形式,那么约束条件CTh=f可以表示为
1.3利用拉格朗日乘子算法求解基于线性约束最小功率准则(LCMP)的波束形成优化问题,得到:
h=Rxx -1C(CTRxx -1C)-1f
其中,CT=Ad,f=bd,令ωi表示第i个频率栅格点的角频率,i=1,2,…,Nf,Nf表示将宽带数字信号的角频率范围离散化后的频率栅格点总个数,L表示每个阵元通道对应的FIR滤波器冲激响应包括的FIR滤波器系数总个数,L≥1;Re表示取实部操作,Im表示取虚部操作;
I表示M×M维单位矩阵,
R
xx表示ML×ML维矩阵,
ω
l表示宽带数字信号的角频率最小值,ω
u表示宽带数字信号的角频率最大值,e(e
jω)=[1,e
-jω,e
-j2ω,…,e
-j(L-1)ω]
T,S
xx(ω)=E{x(e
jω)x
H(e
jω)},上标T表示转置,x(e
jω)=[X
1(e
jω),...,X
m(e
jω),...,X
M(e
jω)]
T,
ω表示角频率变量,ω∈[ω
l,ω
u],E表示求期望,
n-k=0,1,2,…,Γ-1,s[n]表示将宽带数字信号分路为M路后每路宽带数字信号中的第n个元素值,s[n-k]表示对s[n]延迟长度k后的元素值,n=0,1,…,Γ-1,Γ表示将宽带数字信号分路为M路后每路宽带数字信号包括的元素值总个数,k=0,1,…,L-1,
φ
q表示第q个期望形成零点方向,q=1,…,Q,Q表示期望形成零点方向总个数,Q≥1,m=1,...,M,c表示光速,
表示阵元间距,上标H表示共轭转置,e(e
jω)=[1,e
-jω,e
-j2ω,…,e
-j(L-1)ω]
T,上标*表示共轭,
表示张量积。
此外,通过增加约束条件,可以在特定方向φ
q,q=1,…,Q上形成零点,Q≥1;其中φ
q≠φ
d且0°<φ
q<180°,令
a
q(e
jω)表示第q个期望形成零点方向φ
q对应的导向矢量,f
add=[0,…,0]
T,其中Q为期望形成零点方向总个数,则有:
其中,
或者是在计算待发射干扰信号的自相关矩阵Sxx(ω)时,在第q个期望形成零点方向φq上包含一带宽与期望发射宽带数字信号带宽相同的信号,则在使得hTRxxh最小化的过程中,同样会在期望形成零点的方向形成零点,即
其中,σ
q为针对第q个期望形成零点方向φ
q对应的控制参数,取值范围0<σ
q<1,
例如在均匀线阵中,
φ
q表示第q个期望形成零点方向,q=1,…,Q,Q表示期望形成零点方向总个数,Q≥1,m=1,...,M,c表示光速,
表示阵元间距。
1.4对基于线性约束最小功率准则(LCMP)的波束形成器进行进一步优化,可以从获得更窄的主瓣宽度和更低的旁瓣电平进行考虑,从而进一步提高干扰功率的有效利用率,最直接的方法就是增加旁瓣约束;那么根据所述基于线性约束最小功率准则(LCMP)的波束形成优化问题得到基于线性约束最小功率准则(LCMP)的波束形成优化问题转化形式:
subject to Adh=bd
其中,
其中,subject to表示约束条件,
表示第ξ个方位角栅格、第i个频率栅格点对应的导向矢量,上标*表示共轭,
表示张量积,上标H表示共轭转置;ε
i表示第i个频率栅格点的角频率ω
i对应的波束旁瓣的约束值,为未知量;
表示h
TR
xxh取最小值时对应的h。
表示波束旁瓣区域内包括的方位角栅格总个数,对于均匀线阵,波束旁瓣区域内包括的方位角栅格总个数
一般取值为200即可;SLR表示波束旁瓣区域,
φ
ξ表示第ξ个方位角栅格方向,
表示波束旁瓣区域内包括的方位角栅格总个数;ω
i∈[ω
l,ω
u]。
需要说明的是,波束是信号经过阵元发射到空间中的表示形式,分为主瓣和副瓣两部分,主瓣有一定宽度,其宽度范围为主瓣区域,除去主瓣区域剩余的区域为副瓣区域,例如波束主瓣指向为60°,波束主瓣宽度为20°,那么在60°±20°的角度范围为主瓣区域,(0°,50°)∪(70°,180°)的范围为波束旁瓣区域。
像干扰机阵列这样的宽带波束形成器,通常工作在很大的频带范围内,为获得一组对整个频带范围都适用的滤波器系数,通常需要增大将宽带数字信号的角频率范围离散化后的频率栅格点总个数N
f的值,这会导致旁瓣约束式
规模很大,以至于很难获得一个可行的优化解。
为了解决这个问题,可以通过消除
的频率依赖性来优化宽带波束方向图,即将ε
i(i=1,2,…,N
f)共N
f个变量转化为一个变量ε,因此将其改写为:
其中,‖·‖表示欧几里得范数,
α表示波束旁瓣区域功率约束值的真值大小,通常取α<0.22,N
f×ML维矩阵
定义如下:
对ML×ML维矩阵Rxx进行Cholesky分解,可得:
Rxx=VTV
其中,V表示对ML×ML维矩阵Rxx进行Cholesky分解后的上三角矩阵,因此有:
hTRxxh=(Vh)T(Vh)=‖Vh‖2
1.5引入辅助变量η,其中η为未知变量;然后利用二阶锥规划的算法将上述基于线性约束最小功率准则(LCMP)的波束形成优化问题转化形式转化为二阶锥形式问题:
其中,subject to表示约束条件;求解上述二阶锥形式问题,可以获得M个FIR滤波器冲激响应h,M个FIR滤波器冲激响应h是ML×1维向量,h=[h1,...,hm,...,hM]T,hm表示第m个阵元通道对应的FIR滤波器冲激响应,m=1,2,…,M,M表示均匀线阵包括的阵元总个数。
步骤2,将每路宽带数字信号分别经过对应的FIR滤波器冲激响应进行FIR滤波后,得到M路经过FIR滤波后的宽带数字信号。
步骤3,对每路经过FIR滤波后的宽带数字信号分别进行数/模转换、上变频处理、射频滤波处理和放大处理后,得到M路经过幅相处理后的宽带信号;M路经过幅相处理后的宽带信号与M个阵元一一对应。
步骤4,将每路经过幅相处理后的宽带信号分别通过对应阵元发射到空中,然后利用波束形成技术得到宽带干扰波束,所述宽带干扰波束为基于二阶锥规划的宽带发射数字波束形成干扰设计结果。
通过以下仿真实验对本发明效果进行进一步验证说明。
1.仿真条件
本发明的仿真是在主频2.3GHz的Intel Core i5-5500HQ、内存8.0GB的硬件环境和MATLAB R2016以及工具箱CVX1.22的软件环境下进行的。
2.仿真内容
仿真1:
仿真要求:假设期望在70°方向形成波束,发射信号频率范围为[1.9GHz,2.1GHz],即干扰波束的频域带宽为200MHz,主瓣零点波束宽度设定为20°,波束旁瓣对应的归一化输出功率低于-20dB。
主要仿真参数:仿真中采用均匀线阵,阵元个数为16,阵元间距为
λ
min/2=0.0714m(λ
min为最高频率对应的波长),FIR滤波器长度L=100(滤波器长度的设置需保证其通带范围大于或等于干扰信号频带范围),N
f=40(即将干扰波束的频域带宽范围离散化频率栅格点数为40),因要求波束旁瓣对应的归一化输出功率低于-20dB(20lgα=-20dB,则α=0.1),所以旁瓣功率约束参数
***噪声的标准差σ
n=0.01。
仿真效果分析:
由图3(a)至3(c)可知,采用二阶锥规划(SOCP)对附加旁瓣约束的线性约束最小功率准则(LCMP)进行优化,其主瓣零点波束宽度和波束旁瓣功率两个指标的设计结果基本满足了期望实现的干扰波束的设计要求,形成的波束图与传统的线性约束最小功率准则(LCMP)形成的波束图相比,可以明显看出,发射功率集中到了主瓣区域。
虽然设计结果对整个频域范围内基本达到了总的归一化旁瓣功率小于-20dB的要求,但从图3(b)空域剖面图可知,仍然存在频率分量的旁瓣功率值在一定方位角上略大于预设值;由于在波束优化设计问题中,旁瓣功率和主瓣宽度是一对矛盾的设计指标,因此在实际工程应用中,若需要所有的频率分量的旁瓣功率值均低于-20dB或者是其它阈值时,可以通过调整主瓣宽度和旁瓣功率约束参数ε的大小来进一步降低整个频带范围内的旁瓣功率。
仿真2:
仿真要求:在50°、90°和110°方向形成波束零点,其余仿真要求与仿真1中相同。
取零点控制参数σq=0.5,其余主要仿真参数与仿真1中相同
仿真效果分析:
由图4(a)可知,本发明的方法能够在期望形成零点的方向形成较深的零陷,其零陷值均小于或等于-60dB,效果比较理想;结合图4(b)可知,增加波束零点,会使发射波束方向图的副瓣性能发生一定畸变,并且导致主瓣有一定程度的展宽;仿真结果表明,该方法可以较好地解决相应约束条件下雷达宽带干扰波束形成优化设计问题,实现了雷达宽带干扰数字波束形成,验证了该方法的有效性。
综上所述,仿真实验验证了本发明的正确性,有效性和可靠性,本发明的设计方法均能满足预先要求的设计指标,达到了优化设计的预期要求。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围;这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。。