CN114280544B - 一种基于松弛优化的最小过渡带宽度方向图赋形方法 - Google Patents
一种基于松弛优化的最小过渡带宽度方向图赋形方法 Download PDFInfo
- Publication number
- CN114280544B CN114280544B CN202111457818.2A CN202111457818A CN114280544B CN 114280544 B CN114280544 B CN 114280544B CN 202111457818 A CN202111457818 A CN 202111457818A CN 114280544 B CN114280544 B CN 114280544B
- Authority
- CN
- China
- Prior art keywords
- relaxation
- constraint
- theta
- pattern
- shaping
- 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
- 238000000034 method Methods 0.000 title claims abstract description 46
- 230000007704 transition Effects 0.000 title claims abstract description 31
- 238000010586 diagram Methods 0.000 title claims abstract description 17
- 238000005457 optimization Methods 0.000 title claims abstract description 16
- 238000007493 shaping process Methods 0.000 title claims description 25
- 239000013598 vector Substances 0.000 claims description 22
- 238000000354 decomposition reaction Methods 0.000 claims description 16
- 230000003595 spectral effect Effects 0.000 claims description 15
- 238000005311 autocorrelation function Methods 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 5
- 238000001228 spectrum Methods 0.000 claims description 2
- 238000012545 processing Methods 0.000 abstract description 9
- 230000015572 biosynthetic process Effects 0.000 abstract description 5
- 238000003786 synthesis reaction Methods 0.000 abstract description 5
- 238000004088 simulation Methods 0.000 description 18
- 239000011159 matrix material Substances 0.000 description 5
- 230000000694 effects Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000001308 synthesis method Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
该发明公开了一种基于松弛优化的最小过渡带宽度方向图赋形方法,属于阵列雷达信号处理领域。本发明采用其他约束代替原本的雷达阵列无法求解的非凸约束,通过引入松弛变量和相依你的约束条件求解出重构系数,从而实现发明目的。本发明不需要通过经验来设置过渡带的宽度,可以通过迭代来生成最小宽度的过渡带;通过较少的迭代次数就可以生成理想的过渡带,收敛比较快;相对于其他需要迭代的方向图综合算法,运行时间较短。
Description
技术领域
本发明属于阵列雷达信号处理领域,具体地说是一种利用松弛优化的方式,可以不依靠经验地生成带有最小宽度过渡带的赋形方向图方法。
背景技术
天线在日常生活中随处可见,其中,阵列天线凭借其方向性强和高增益的特性在工程应用当中颇受欢迎,且阵列信号处理作为信号处理的一个重要分支,在雷达探测、无线通信、地质勘探等诸多军用及民用方面有着广泛的应用。阵列方向图综合在阵列信号处理中也是一个很重要的研究方向,它是一种使用传感器阵列定向发送和接收信号的信号处理技术,对阵列***的高性能发挥起到重要的作用。近几十年来,国内外对阵列方向图综合的研究一直在进行,并出现了多种有效方法,大体可以分为优化方法和解析方法两大类,而由于凸优化理论的发展,可以针对凸函数最小化的问题求得全局最优解,并且伴随着各种数值求解工具的出现,凸优化方法将更高效且可靠的应用阵列信号处理等各工程领域之中。本发明将松弛思想引入方向图综合中,利用凸优化理论进行最优求解最后得到具有最小过渡带宽度的赋形方向图。
发明内容
本发明的目的在于针对上述已有技术的不足,提出一种基于松弛优化实现最小过渡带宽度的赋形方向图综合方法,来生成理想的方向图。
为实现上述目的,本发明的技术方案为一种基于松弛优化的最小过渡带宽度方向图赋形方法,该方法包括如下步骤:
步骤1:定义N个阵元的均匀线阵雷达,生成的方向图函数为f(θ)=wHa(θ),其中,上标H表示共轭转置运算,w为方向图的权向量,a(θ)为阵列的导向矢量;主瓣区域的方向图赋形约束为l(θ)≤|f(θ)|2≤u(θ),θ∈Ωm,其中l(θ)和u(θ)分别表示主瓣电平的上下界;旁瓣区域的方向图赋形约束为|f(θ)|2≤ρ(θ),θ∈Ωs,其中ρ(θ)表示旁瓣电平的上界,Ωm和Ωs分别代表主瓣区域以及旁瓣区域,详细的方向图赋形区间示意图可见图2;
步骤3:设置两个松弛区间Ωr,1,Ωr,2,对这两个区间内的方向图施加松弛变量γ(θ),得到新的松弛约束上标T表示转置运算,且在Ωr,1区间上的松弛变量保持单调不减,在Ωr,2区间内的松弛变量保持单调不增;
步骤4:采用如下目标函数代替过渡带宽度最小化问题;
步骤5:对步骤4的公式,施加相应的方向图赋形约束,通过CVX工具包进行求解得到理想的rw,根据谱分解的原理,以rw的每个元素为系数构造多项式P(ejθ),并取其单位圆内的根,记为σ1,σ2,…,σN-1,再以这一组数为根构造出新的多项式S(ejθ),结合步骤2中w与rw的对应关系重构出此时S(ejθ)的系数便是权向量w的元素。
进一步的,对所述步骤4中目标函数施加相应的方向图赋形约束后的为:
rw(0)≥0
Dγ≤0
γ(θ)≥0
其中,Dγ≤0表示松弛变量的单调性约束,rw(·)表示rw中的元素。
进一步的,所述步骤5中权向量w的重构方法为:
得到rw后,以rw的元素为系数构造多项式:
并将其单位圆内的根表示为σ1,σ2,…,σN-1,根据谱分解原理得到如下关系式:
其中
S(ejθ)的系数即为权向量w的元素,C为常数。
本发明与现有技术相比具有以下优点:
本发明不需要通过经验来设置过渡带的宽度,可以通过迭代来生成最小宽度的过渡带;通过较少的迭代次数就可以生成理想的过渡带,收敛比较快;相对于其他需要迭代的方向图综合算法,运行时间较短。
附图说明
图1是本发明的流程图;
图2是本发明的方向图赋形区间示意图;
图3是本发明与半正定松弛和谱分解方法在不同条件下生成的方向图;
图4是本发明与半正定松弛和谱分解方法在旁瓣有凹口的条件下生成的方向图;
图5是本发明与半正定松弛和谱分解方法在线性顶赋形的条件下生成的方向图。
具体实施方式
参照图1,本发明的具体实现步骤如下:
步骤1,通过约束权向量的自相关函数rw进行主瓣和旁瓣区间的方向图赋形。
考虑一个N个阵元的均匀线阵,它的方向图函数表示为:
f(θ)=wHa(θ)
其中w=[w0,w1,…,wN-1]T为方向图的权向量,a(θ)=[1,ejΦ,…,ej(N-1)Φ]T为阵列的导向矢量,参照图2,通过以下步骤实现:
赋形方向图的和旁瓣区间的约束表示为:
l(θ)≤|f(θ)|2≤u(θ),θ∈Ωm
|f(θ)|2≤ρ(θ),θ∈Ωs
其中l(θ)和u(θ)分别表示主瓣电平的上下界,ρ(θ)表示旁瓣电平的上界,Ωm和Ωs分别代表主瓣区域以及旁瓣区域;
步骤2:引入一个长度为2N-1的新的导向矢量:
b(θ)=[e-j(N-1)Φ,…,e-jΦ,1,ejΦ,…,ej(N-1)Φ]T
由谱分解原理得到如下关系式:
其中
rw=[rw(-N+1),…,rw(0),…,rw(N-1)]T
表示权向量的自相关函数,将步骤1中对方向图函数的约束转化为对rw的约束:
这样不仅满足了对幅度响应的约束而且还将非凸约束转化为了凸约束。
此外,由于自相关函数的性质,不难发现:
rw(0)≥0
步骤3,设定松弛区间Ωr,1,Ωr,2并对其进行松弛约束。
设定如图2中所示的松弛区间Ωr,1,Ωr,2,并在其中引入松弛变量γ(θ),且同步骤1中一样的对自相关函数进行约束得到对于松弛区间内方向图的松弛约束:
步骤4,对松弛变量进行单调性约束。
参照图2,根据方向图过渡带的形状通过如下步骤实现单调性约束:
3a)不妨先将γ(θ)离散化:
由过渡带的变化趋势可知,Ωr,1中的松弛变量应该是单调不减的,Ωr,2中的松弛变量应该是单调不增的,即满足
其中p=1,2,…,R1-1和q=1,2,…,R2-1。
3b)为了更直观的表达,且便于操作,将松弛变量表示为:
并且构建一个差分结构矩阵:
其中,
矩阵V的下标R对应看矩阵的列数。
通过矩阵的方式将松弛变量的单调性约束表示为:
Dγ≤0
步骤5,将过渡带宽度最小化问题转化为对松弛变量γ的稀疏优化问题。
4a)由于松弛变量的元素γ(θ)是可以取到0的,松弛变量的0值越多,过渡带宽度就会相应的越窄,所以,将其转化为如下的稀疏优化问题:
4b)看出l0范数的最小化问题并不是凸优化问题,用l1范数最小化来替代它,即:
min||γ||1
为了进一步的强化γ的稀疏性,采用l1重加权的方法进行迭代求解,再加上之前设计的方向图赋形以及松弛约束得到一个完整的凸优化问题:
rw(0)≥0
Dγ≤0
γ(θ)≥0
进一步通过CVX工具包求得所需的rw。
利用谱分解原理重构权向量w并画出方向图。
得到rw后,以rw的元素为系数构造多项式:
并将其单位圆内的根表示为σ1,σ2,…,σN-1;根据谱分解原理得到如下关系式:
其中
结合1b)中的谱分解关系式,得到,S(ejθ)的系数即为权向量w的元素,进而得到设计的方向图。
一.仿真条件及仿真数据处理
1.仿真条件
设置仿真参数如表1所示:
表1仿真参数表
2.仿真数据处理
仿真1,将本发明与半正定松弛和谱分解方法在不同条件下生成的方向图对比。
半正定松弛(Semi-Definite Relaxation,SDR)方法也是通过约束主瓣区域的上下界进行赋形,并且通过构造一个对称半正定矩阵来对方向图的响应进行松弛优化;谱分解(Spectral Factorization,SF)方法是利用谱分解的原理构造出权向量直接画出方向图。
将旁瓣上限电平设为-30dB,主瓣赋形区域设置为[-30°,30°]。第一种情况下,将SDR和SF方法中的过渡带区间设置为[-50°,-30°]以及[30°,50°];第二种情况下,将SDR和SF方法中的过渡带区间设置为[-40°,-30°]以及[30°,40°];第三种情况下,将SDR和SF方法中的过渡带区间设置为[-35°,-30°]以及[30°,35°]。对于本发明的仿真则不需要设置过渡带宽度,只需要设定一个较宽的松弛区间,这里设为[-50°,-30°]以及[30°,50°]。仿真结果见图3。在上述的三种不同情况下,三种算法生成的主瓣宽度和运行时间如表2所示。
表2本发明与SDR、SF方法仿真对比
由仿真结果可见,相对于SDR方法和SF方法,本发明生成的主瓣宽度要窄一些,且在同样的我迭代次数下,要比SDR方法运行时间快。
仿真2,主瓣左侧旁瓣上限电平不变,在[-60°,-40°]上加入一个-40dB的凹口,主瓣右侧旁瓣上限电平设为-35dB,主瓣赋形区域为[-20°,30°],对SDR和SF方法的过渡带宽度设为[-30°,-20°]以及[30°,40°],其结果如图4所示。同样的,在这种情况下,三种方法的主瓣宽度及运行时间如表3。
表3本发明与SDR、SF方法仿真对比
由仿真结果可见,本发明在相对复杂的条件下依然可以保持良好的效果,生成最小的过渡带宽度。
仿真3,将主瓣赋形区域设为[-30°,30°],且在这个范围内令其线性的从0dB降为-12dB,旁瓣上限设为-30dB,对于SDR和SF方法的过渡带宽度设为[-40°,-30°]以及[30°,40°],仿真结果见图5。同样地,表4给出不同方法的过渡带宽度和运行时间。
表4本发明与SDR、SF方法仿真对比
由仿真结果可见,本发明在线性主瓣可以保持良好的效果,生成最小的过渡带宽度。
综上,采用本发明生成的赋形方向图不需要通过经验设定过渡带宽度,可以通过少量迭代生成最小宽度的过渡带。
Claims (2)
1.一种基于松弛优化的最小过渡带宽度方向图赋形方法,该方法包括如下步骤:
步骤1:定义N个阵元的均匀线阵雷达,生成的方向图函数为f(θ)=wHa(θ),其中,上标H表示共轭转置运算,w为方向图的权向量,a(θ)为阵列的导向矢量;主瓣区域的方向图赋形约束为l(θ)≤|f(θ)|2≤u(θ),θ∈Ωm,其中l(θ)和u(θ)分别表示主瓣电平的上下界;旁瓣区域的方向图赋形约束为|f(θ)|2≤ρ(θ),θ∈Ωs,其中ρ(θ)表示旁瓣电平的上界,Ωm和Ωs分别代表主瓣区域以及旁瓣区域;
步骤3:设置两个松弛区间Ωr,1,Ωr,2,对这两个区间内的方向图施加松弛变量γ(θ),得到新的松弛约束上标T表示转置运算,且在Ωr,1区间上的松弛变量保持单调不减,在Ωr,2区间内的松弛变量保持单调不增;
步骤4:采用如下目标函数代替过渡带宽度最小化问题;
目标函数施加相应的方向图赋形约束后的为:
rw(0)≥0
γ(θ)≥0
步骤5:对步骤4的公式,施加相应的方向图赋形约束,通过CVX工具包进行求解得到理想的rw,根据谱分解的原理,以rw的每个元素为系数构造多项式P(ejθ),并取其单位圆内的根,记为σ1,σ2,…,σN-1,再以这一组数为根构造出新的多项式S(ejθ),结合步骤2中w与rw的对应关系重构出此时S(ejθ)的系数便是权向量w的元素。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111457818.2A CN114280544B (zh) | 2021-12-02 | 2021-12-02 | 一种基于松弛优化的最小过渡带宽度方向图赋形方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111457818.2A CN114280544B (zh) | 2021-12-02 | 2021-12-02 | 一种基于松弛优化的最小过渡带宽度方向图赋形方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114280544A CN114280544A (zh) | 2022-04-05 |
CN114280544B true CN114280544B (zh) | 2023-06-27 |
Family
ID=80871032
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111457818.2A Active CN114280544B (zh) | 2021-12-02 | 2021-12-02 | 一种基于松弛优化的最小过渡带宽度方向图赋形方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114280544B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107728118A (zh) * | 2017-09-25 | 2018-02-23 | 西北工业大学 | 无需拟合协方差矩阵的低旁瓣发射波束图设计方法 |
CN110032805A (zh) * | 2019-04-16 | 2019-07-19 | 电子科技大学 | 一种地理约束下的稀疏线性阵列优化布阵方法 |
CN113704998A (zh) * | 2021-08-27 | 2021-11-26 | 电子科技大学 | 一种基于交替方向乘子法的波束形成方法 |
Family Cites Families (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB0906269D0 (en) * | 2009-04-09 | 2009-05-20 | Ntnu Technology Transfer As | Optimal modal beamformer for sensor arrays |
US20140180629A1 (en) * | 2012-12-22 | 2014-06-26 | Ecole Polytechnique Federale De Lausanne Epfl | Method and a system for determining the geometry and/or the localization of an object |
US9261592B2 (en) * | 2014-01-13 | 2016-02-16 | Mitsubishi Electric Research Laboratories, Inc. | Method and system for through-the-wall imaging using compressive sensing and MIMO antenna arrays |
CN105655727B (zh) * | 2015-09-28 | 2018-07-24 | 中国电子科技集团公司第二十九研究所 | 一种基于凸优化的增益损失恒定波束的形成方法及装置 |
CN105467365B (zh) * | 2015-12-08 | 2018-07-17 | 中国人民解放军信息工程大学 | 一种改善mimo雷达doa估计性能的低旁瓣发射方向图设计方法 |
US10545232B2 (en) * | 2016-08-24 | 2020-01-28 | Mitsubishi Electric Research Laboratories, Inc. | System and method for radar imaging using distributed moving platforms under unknown position perturbations |
CN106886656B (zh) * | 2017-03-15 | 2020-12-25 | 南京航空航天大学 | 一种三维阵列天线方向图旁瓣抑制方法 |
CN107728112B (zh) * | 2017-09-29 | 2020-08-14 | 西安电子科技大学 | 在目标导向矢量严重失配情况下的稳健波束形成方法 |
CN108337030B (zh) * | 2018-02-07 | 2020-06-30 | 东南大学 | 多天线***中高功效波束成形方法、装置及设备 |
CN108809397B (zh) * | 2018-06-27 | 2020-06-30 | 东南大学 | 多天线***中高功效数模混合波束成形方法、装置及设备 |
CN110133613B (zh) * | 2019-05-17 | 2021-03-16 | 电子科技大学 | 一种基于半定松弛技术的频控阵稳健波束形成方法 |
CN111553095B (zh) * | 2020-06-09 | 2024-03-19 | 南京航空航天大学 | 基于序列二阶锥算法的时间调制阵列边带抑制方法 |
CN111817765B (zh) * | 2020-06-22 | 2021-12-03 | 电子科技大学 | 一种基于频率约束的广义旁瓣对消宽带波束形成方法 |
CN113094904B (zh) * | 2021-04-12 | 2022-09-16 | 电子科技大学 | 一种非凸多约束下的阵列天线波束赋形优化方法 |
CN113673317B (zh) * | 2021-07-12 | 2023-04-07 | 电子科技大学 | 基于原子范数最小化可降维的二维离格doa估计方法 |
CN113655436B (zh) * | 2021-07-23 | 2024-04-16 | 公安部第一研究所 | 一种带通道校准粒子群优化宽带波束形成方法及装置 |
-
2021
- 2021-12-02 CN CN202111457818.2A patent/CN114280544B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107728118A (zh) * | 2017-09-25 | 2018-02-23 | 西北工业大学 | 无需拟合协方差矩阵的低旁瓣发射波束图设计方法 |
CN110032805A (zh) * | 2019-04-16 | 2019-07-19 | 电子科技大学 | 一种地理约束下的稀疏线性阵列优化布阵方法 |
CN113704998A (zh) * | 2021-08-27 | 2021-11-26 | 电子科技大学 | 一种基于交替方向乘子法的波束形成方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114280544A (zh) | 2022-04-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108375752B (zh) | 基于全角度搜索的幅相误差单辐射源测向方法 | |
Fan et al. | MIMO radar waveform design for quasi-equiripple transmit beampattern synthesis via weighted $ l_p $-minimization | |
CN107589399B (zh) | 基于多采样虚拟信号奇异值分解的互质阵列波达方向估计方法 | |
CN107257254B (zh) | 一种安全保密通信***中的恒模波束赋形方法 | |
KR101095996B1 (ko) | 부호교차 위상제어 방식의 위상배열 안테나 빔 널링 방법 | |
CN107395255B (zh) | 一种基于凸优化的稳健混合波束成形方法 | |
CN104539340B (zh) | 一种基于稀疏表示和协方差拟合的稳健波达角估计方法 | |
CN112653496B (zh) | 一种毫米波大规模mimo***的混合预编码方法 | |
CN107231180B (zh) | 一种基于波束成形的动态方向调制方法 | |
Bakhar et al. | Robust blind beam formers for smart antenna system using window techniques | |
CN103715518B (zh) | 泰勒-指数复合不等距模组化阵列天线的设计方法 | |
CN106526565A (zh) | 一种基于支持向量机的单比特空间谱估计方法 | |
CN114280544B (zh) | 一种基于松弛优化的最小过渡带宽度方向图赋形方法 | |
CN113671439B (zh) | 基于非均匀智能超表面阵列的无人机集群测向***及方法 | |
Guo et al. | Atmospheric turbulence compensation for OAM-carrying vortex waves based on convolutional neural network | |
CN106021766B (zh) | 基于机电集成优化的圆极化网状天线高指向精度设计方法 | |
CN110444903A (zh) | 基于多环阵列的低副瓣涡旋波束产生方法 | |
CN113871899A (zh) | 非均匀间隔线阵的宽带非频变扫描波束快速生成方法 | |
CN111162878B (zh) | 一种基于子阵降维带约束的多域联合抗干扰方法 | |
CN110083923B (zh) | 一种基于高阶泰勒展开的低副瓣阵列天线的优化布局方法 | |
Zhang et al. | Channel estimation for mmWave massive MIMO systems with mixed-ADC architecture | |
CN108170888B (zh) | 基于最小化加权向量动态范围的波束图综合设计方法 | |
CN115986404A (zh) | 一种宽带mimo雷达波束扫描与天线位置选择方法 | |
CN105044679A (zh) | 一种基于正交投影的发射波束零陷展宽方法 | |
CN115510610A (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 |