稀疏波束图综合设计方法技术

技术编号:15391812 阅读:115 留言:0更新日期:2017-05-19 04:57
本发明专利技术公开的稀疏波束图综合设计方法,具体为:步骤1、在实际阵列信号系统中,依据设计要求,建立以p范数(0<p<2)为稀疏约束的权值优化目标模型,并依据要求设置主瓣宽度,纹波宽度以及旁瓣限制的约束条件;步骤2、待步骤1完成后,针对该非凸优化模型,采用交替方法乘子法来迭代求出稀疏权值系数w;步骤3、利用步骤2得到的稀疏权值系数w并结合计算机仿真软件按照设计要求作出波束图本发明专利技术稀疏波束图综合设计方法,解决了现有技术中存在的局限性和不精确性,并且降低了实际阵列系统的复杂度以及实际应用的成本。

Sparse beam pattern synthesis design method

The invention discloses a sparse beam pattern synthesis method, in particular: Step 1, in the actual array signal system, according to the design requirements, based on the p norm (0< p< 2) for weight optimization target model sparse constraint, and according to the requirements set the width of the main lobe width and the side lobe, ripple limit constraints step 2, to step 1; after the completion of the non convex optimization model to iterative sparse weight coefficient W method using alternating multiplier method; step 3 step, using sparse weight coefficients of W 2 are combined with computer simulation software to make beam whose invention beam pattern synthesis design method of sparse solution according to design requirements. The limitations of the existing technology and accuracy, and reduce the actual array system complexity and cost of practical application.

【技术实现步骤摘要】
稀疏波束图综合设计方法
本专利技术属于阵列信号处理方法
,具体涉及一种稀疏波束图综合设计方法。
技术介绍
近年来,阵列信号处理技术在无线通信领域应用广泛,其中一个重要部分就是波束形成技术,其实质是通过对各阵元加权进行空域滤波,来达到增强期望信号及抑制干扰的目的。波束图综合设计技术是阵列信号处理领域的关键技术,已广泛应用于雷达、声呐、电子监视及地震探测等领域。当期望来波方向估计不精确或不明确时,就可以通过控制辐射波束图的主瓣波束宽度和响应纹波来提高估计的鲁棒性。传统的波束图合成方法并不能完全任意的设计旁瓣下界或者仅仅能应用在均匀间隔的阵列中,抑或不能精确的控制波束合成要求。另外,对于波束合成,要尽可能使用少量天线,来减轻重量、节省成本以及减少阵列系统的复杂度。因此,稀疏波束图综合设计技术在实际应用中有很多优势,克服了传统合成方法的缺点,可以精确设计任意旁瓣要求并且减少了实际应用的成本。
技术实现思路
本专利技术的目的在于提供一种稀疏波束图综合设计方法,解决了现有技术中存在的局限性和不精确性,并且降低了实际阵列系统的复杂度及实际应用的成本。本专利技术所采用的技术方案是,稀疏波束图综合设计方法,具体按照以下步骤实施:步骤1、在实际阵列信号系统中,依据设计要求,建立以p范数0<p<2为稀疏约束的权值优化目标模型,并依据要求设置主瓣宽度,纹波宽度以及旁瓣限制的约束条件;步骤2、待步骤1完成后,针对得到的非凸优化模型,采用交替方法乘子法来迭代求出稀疏权值系数w;步骤3、利用步骤2得到的稀疏权值系数w并结合计算机内的仿真软件,按照设计要求作出波束图。本专利技术的特点还在于:步骤1具体按照以下方法实施:对于一维阵列,有N个传感器,则在角度θ方向上的波束模式具体如下:p(θ)=wHa(θ)(1);在式(1)中:w=[w1w2,…wN]T为权值向量,a(θ)=[a1(θ),…,aN(θ)]T为导向矢量;设定将主瓣区域均匀划分为M个角度,其宽度为:{θ1,…,θm},m=1,…,M,则主瓣纹波宽度约束具体为:1-ε≤|wHa(θm)|2≤1+ε(2);在式(2)中:ε为纹波系数;同理,设定将旁瓣区域被均匀划分为S个角度,其宽度为s=1,…,S,则旁瓣上界约束具体如下:式(3)中,η为上界系数;建立如下稀疏约束的权重优化目标模型,具体如下:步骤2具体按照以下步骤实施:步骤2.1、根据ADMM算法,引入辅助变量um,vs,则式(4)变形为如下形式:构造如下的拉格朗日函数:其中ρ是迭代步长,是拉格朗日乘子,分别为取复数实部和虚部;步骤2.2、进行优化处理,具体按照以下步骤实施:步骤①、固定w(t),由如下算法求解出vs,um,具体算法如下:将式(6)简化为如下形式:在式(7)中:将式(7)用复数形式表示为如下形式:在式(9)中:显然,式(9)能分解为M+S个子问题来求解,具体如下:和对式(12)进行分解,设定式(12)的最小值点在|um|2=μ处,μ是区间[1-ε,1+ε]处的常数,则目标函数表示为如下形式:在式(14)中:(·)*是满足KKT条件的平衡点;式(14)表明最优点um与xm同相位且幅值为即:将式(15)代入式(14)中,则将式(12)优化为如下形式:式(16)的解具体如下:将式(17)代入式(15),则得到如下形式:同理,得到式(13)的解具体如下:步骤②、固定um(t+1),vs(t+1),由如下算法求解出w:进一步的得到:式(20)的矩阵表示为如下形式:在式(22)中:引入辅助变量z,以求解式(22),令z=w,则得到如下形式:根据式(24),构造拉格朗日函数:应用ADMM算法求解式(25),具体按照以下步骤实施:步骤I、固定w(t),λr(t),λi(t),利用如下算法求解出z:式(26)中:进一步,式(26)等价于如下的最小二乘问题:在式(28)中:I为单位阵;则得到:步骤II、待步骤I完成后,固定z(t+1),λr(t),λi(t),利用如下算法求解出w:在式(30)中:式(30)能分离为N个子问题:在式(31)中:w=[w1,w2…,wN]T,在式(31)最小时,wn与有相同的相位,即βn,则wn能表示为如下形式:式(32)中,hn是wn的幅值,根据wn的幅值大小hn,即能得到式(31)的最优点,则能将式(31)简化为如下形式:式(33)中,是的幅值;由于S(hn)在上分别单调递增和递减,则S(hn)的全局最小值出现在上,即式(33)等价于如下形式:判断S(hn)在上的凹凸性,并分别对S(hn)求1阶导、2阶导,3阶导:对p分三种情况讨论式(34)的最小值点:第一种情况:1<p≤2;时,故S(hn)是凸的,并且所以S(hn)在区间有唯一解,使用二分法对式(35)求零点,即为S(hn)的最小值点;第二种情况:p=1;当p=1时,式(34)等价于如下算法:则S(hn)的最小值点为:第三种情况:0<p<1;0<p<1时故是单增的,S(hn)的凹凸性取决于的符号,记的解如下:当S(hn)是凹的,故S(hn)的最小值点在0或之间;当既可能是正的也可能是负的,所以要分和两个区间讨论,具体分别如下:在上,S(hn)是凹的,故S(hn)的最小值点在0或之间;在上,S(hn)是凸的,且若S(hn)是单增的,故上的局部最小值点若使用二分法求取上的局部最小值点;显然,S(hn)的全局最小值点通过比较和上的局部最小值而获得;步骤III、固定z,w,利用如下算法求解出λ:λ(t+1)=λ(t)+ρ(w-z)(41);循环迭代步骤I、步骤II及步骤III,直到z,w,λ收敛。步骤③、固定um(t+1),vs(t+1),w,利用如下算法求解出λm,κs:λm(t+1)=λm+ρ(um-wHa(θm))(42);步骤2.3、经步骤2.2中的步骤①、步骤②及步骤③循环迭代一定次数,直到um,vs,w,λm,κs收敛,即得到稀疏权值系数w。本专利技术的有益效果在于:(1)本专利技术稀疏波束图综合设计方法有效克服了现有技术存在的局限性和不精确性,并且降低了实际阵列系统的复杂度。(2)本专利技术稀疏波束图综合设计方法,能精确设计任意旁瓣要求。(3)本专利技术稀疏波束图综合设计方法减少了实际应用的成本。(4)本专利技术稀疏波束图综合设计方法,非常适合在阵列信号处理中推广使用。附图说明图1是在41个均匀间隔分布非对称的阵列系统中稀疏16个传感器的仿真波束图;图2是图1波束图主瓣区间的放大图;图3是在41个非均匀间隔对称阵列系统中稀疏14个传感器的仿真波束图;图4是图3波束图主瓣区间的放大图。具体实施方式下面结合附图和具体实施方式对本专利技术进行详细说明。本专利技术稀疏波束图综合设计方法,具体按照以下步骤实施:步骤1、在实际阵列信号系统中,依据设计要求,建立以p范数(0<p<2)为稀疏约束的权值优化目标模型,并依据要求设置主瓣宽度,纹波宽度以及旁瓣限制的约束条件,具体按照以下方法实施:对于一维阵列,有N个传感器,则在角度θ方向上的波束模式具体如下:p(θ)=wHa(θ)(1);在式(1)中:w=[w1w2,…wN]T为权值向量,a(θ)=[a1(θ),…,aN(θ)]T为导向矢量;设定将主瓣区域均匀划分为M个角度,其宽度为本文档来自技高网...
稀疏波束图综合设计方法

【技术保护点】
稀疏波束图综合设计方法,其特征在于,具体按照以下方法实施:步骤1、在实际阵列信号系统中,依据设计要求,建立以p范数0<p<2为稀疏约束的权值优化目标模型,并依据要求设置主瓣宽度,纹波宽度以及旁瓣限制的约束条件;步骤2、待步骤1完成后,针对得到的非凸优化模型,采用交替方法乘子法来迭代求出稀疏权值系数w;步骤3、利用步骤2得到的稀疏权值系数w并结合计算机内的仿真软件,按照设计要求作出波束图。

【技术特征摘要】
1.稀疏波束图综合设计方法,其特征在于,具体按照以下方法实施:步骤1、在实际阵列信号系统中,依据设计要求,建立以p范数0<p<2为稀疏约束的权值优化目标模型,并依据要求设置主瓣宽度,纹波宽度以及旁瓣限制的约束条件;步骤2、待步骤1完成后,针对得到的非凸优化模型,采用交替方法乘子法来迭代求出稀疏权值系数w;步骤3、利用步骤2得到的稀疏权值系数w并结合计算机内的仿真软件,按照设计要求作出波束图。2.根据权利要求1所述的稀疏波束图综合设计方法,其特征在于,所述步骤1具体按照以下方法实施:对于一维阵列,有N个传感器,则在角度θ方向上的波束模式具体如下:p(θ)=wHa(θ)(1);在式(1)中:w=[w1w2,…wN]T为权值向量,a(θ)=[a1(θ),…,aN(θ)]T为导向矢量;设定将主瓣区域均匀划分为M个角度,其宽度为:{θ1,…,θm},m=1,…,M,则主瓣纹波宽度约束具体为:1-ε≤|wHa(θm)|2≤1+ε(2);在式(2)中:ε为纹波系数;同理,设定将旁瓣区域被均匀划分为S个角度,其宽度为则旁瓣上界约束具体如下:式(3)中,η为上界系数;建立如下稀疏约束的权重优化目标模型,具体如下:3.根据权利要求1所述的稀疏波束图综合设计方法,其特征在于,所述步骤2具体按照以下步骤实施:步骤2.1、根据ADMM算法,引入辅助变量um,vs,则式(4)变形为如下形式:构造如下的拉格朗日函数:其中ρ是迭代步长,是拉格朗日乘子,分别为取复数实部和虚部;步骤2.2、进行优化处理,具体按照以下步骤实施:步骤①、固定由如下算法求解出vs,um,具体算法如下:将式(6)简化为如下形式:在式(7)中:将式(7)用复数形式表示为如下形式:在式(9)中:显然,式(9)能分解为M+S个子问题来求解,具体如下:对式(12)进行分解,设定式(12)的最小值点在|um|2=μ处,μ是区间[1-ε,1+ε]处的常数,则目标函数表示为如下形式:在式(14)中:(·)*是满足KKT条件的平衡点;式(14)表明最优点um与xm同相位且幅值为即:将式(15)代入式(14)中,则将式(12)优化为如下形式:式(16)的解具体如下:将式(17)代入式(15),则得到如下形式:同理,得到式(13)的解具体如下:步骤②、固定um(t+1),vs(t+1),由如下算法求解出w:

【专利技术属性】
技术研发人员:梁军利张旋
申请(专利权)人:西北工业大学
类型:发明
国别省市:陕西,61

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1