CN107703489B - 一种mimo雷达恒模波形和接收机的联合设计方法 - Google Patents
一种mimo雷达恒模波形和接收机的联合设计方法 Download PDFInfo
- Publication number
- CN107703489B CN107703489B CN201710878926.4A CN201710878926A CN107703489B CN 107703489 B CN107703489 B CN 107703489B CN 201710878926 A CN201710878926 A CN 201710878926A CN 107703489 B CN107703489 B CN 107703489B
- Authority
- CN
- China
- Prior art keywords
- receiver
- waveform
- radar
- signal
- target
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000013528 artificial neural network Methods 0.000 claims abstract description 6
- 239000013598 vector Substances 0.000 claims description 13
- 238000005457 optimization Methods 0.000 claims description 12
- 210000002569 neuron Anatomy 0.000 claims description 11
- 230000003190 augmentative effect Effects 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 4
- 238000002592 echocardiography Methods 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 abstract description 13
- 238000001514 detection method Methods 0.000 description 6
- 239000011159 matrix material Substances 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000001052 transient effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/28—Details of pulse systems
- G01S7/282—Transmitters
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供了一种MIMO雷达恒模波形和接收机的联合设计方法,涉及MIMO雷达技术领域,本发明提供一种基于拉格朗日规划神经网络算法的MIMO雷达恒模波形和接收机的联合设计,因为LPNN算法是对MIMO雷达恒模波形的直接求解,因此避免了SDR方法中因摒弃秩为1的条件使得求解的信号不是最优的情况,从而提高了输出SINR比,改善了目标参数的估计能力,本发明能准确地确定目标位置,同时最大限度地抑制干扰,而且相比于SDR方法,本发明在直接求解非凸非线性的问题时发挥了很大的优势,这是其他方法所不能达到的。
Description
技术领域
本发明涉及MIMO雷达技术领域,尤其是一种联合设计方法。
背景技术
与传统的相控阵雷达不同,多输入多输出(MIMO)雷达通过同时发射多个波形以提高目标检测、参数估计等多方面的能力,基于天线阵列配置的差异,MIMO雷达可分为集中式MIMO雷达和分布式MIMO雷达。MIMO雷达的性能优越性来自于两个方面:一个是天线阵列的配置,一个是波形分集能力,而多波形设计是实现MIMO雷达波形分集的重要途径,有助于进一步提高MIMO雷达在目标检测、参数估计等方面的能力,因此MIMO雷达的多波形设计问题越来越受到人们的关注。
工程实际应用中为了使雷达发射机能够工作在饱和状态以便发挥其最大效能,同时为避免放大器的非线性使得发射波形失真,通常要求雷达发射波形具有恒定包络。对于MIMO雷达,这一实际要求也不例外,而且由于MIMO雷达包含了多个发射机,每个发射机均可以发射相互独立的波形,因此要求每一路发射波形都满足恒定包络特性。
基于输出信干噪比(Signal-to-Interference-plus-Noise Ratio,SINR)研究MIMO雷达的恒模波形设计问题为NP-hard问题,不易得到全局最优解,此外该问题是非凸的、非线性的。目前该问题的求解主要是通过交替优化算法联合设计波形和接收滤波器,其中该算法在对恒模波形设计时,引入了半正定松弛方法(SDR),SDR方法的原理是通过引入一个秩为1的赫尔米特矩阵变量将关于矢量信号的问题转变成一个关于矩阵信号的问题,同时为了构造一个可通过SDR方法求解的凸优化问题需要摒弃秩为1的非凸约束条件,最后,从求解得到的矩阵信号中恢复矢量信号。显而易见,SDR方法求解得到的是原始问题的近似解,除非有先验信息,否则恢复的矢量信号不是最优。SDR方法使得输出的SINR无法达到最大,而且对干扰的抑制能力和目标的侦测能力也会有所下降,因此有必要寻找一种算法可以基于最大化输出SINR直接求解波形矢量信号,进而提高对目标参数的估计能力。
发明内容
为了克服现有技术的不足,本发明提供一种基于拉格朗日规划神经网络(Lagrange Programming Neural Network,LPNN)算法的MIMO雷达恒模波形和接收机的联合设计,因为LPNN算法是对MIMO雷达恒模波形的直接求解,因此避免了SDR方法中因摒弃秩为1的条件使得求解的信号不是最优的情况,从而提高了输出SINR比,改善了目标参数的估计能力。
本发明解决其技术问题所采用的技术方案包括如下步骤:
步骤一、集中式窄带MIMO雷达***由NT个发射天线和NR个接收天线组成,每一部发射机对应一路发射天线,每个发射机独立地通过全方位发射的方式发射波形,将NT个天线发射的波形记为s(n),n=1,2,...,N,N为每一路波形采样后的长度,到达位于方位角θ处的目标的信号为:
at(θ)Ts(n) (1)
其中at(θ)为发射阵列导向矢量;
假设目标的潜在区域中共有K+1个散射点,其中方位角θ0方向的一个散射点为目标,方位角θk方向的散射点为干扰,其中k=1,2,...,K,且θ0≠θk,则雷达接收机中的信号为所有散射点回波之和,即
其中,x(n)是雷达接收机接收的信号,α0和αk分别为方位角θ0方向的目标的幅度和第k个干扰源的幅度,为接收阵列的导向矢量,为发射阵列导向矢量,d为两天线之间的距离,λ是波长,v(n)为高斯噪声,服从均值为0,方差为
令x=[x(1)T,...,x(N)T]T,s=[s(1)T,...,s(N)T]T,v=[v(1)T,...,v(N)T]T,则式(2)写为:
步骤二、采用最大化***输出SINR作为雷达波形设计的准则,将接收机采用的权值记作ω,则经过接收机处理后的雷达信号输出r为:
其中,H为共轭转置;
故***输出信噪比模型为:
由于SINR同时依赖于接收机权值ω以及发射波形s,故为了使SINR最大,联合优化发射波形和接收机权值,故优化问题为:
将式(6)改写为如下形式:
步骤三、将复数形式的优化问题转化为实值形式的优化问题,因为|ωHA(θm)s|2=|sHAH(θm)ω|2,m=0,1,...,K,将公式(7)中的复数实数化等价为:
其中,λi,i=1,...,NTN为拉格朗日乘子,C0是增广拉格朗日常数项,C0大于0;
式(8)的求解过程为使用公式(14)、(15)和(16)循环迭代更新求解,直至公式(14)、(15)和(16)前后两次迭代值的差值小于10-8,则停止循环迭代。
本发明的有益效果在于能准确地确定目标位置,同时最大限度地抑制干扰,而且相比于SDR方法,本发明避免了SDR方法中因摒弃秩为1的条件使得求解的信号不是最优的情况,从而提高了输出SINR比,改善了目标参数的估计能力,在直接求解非凸非线性的问题时发挥了很大的优势,这是其他方法所不能达到的。
附图说明
图1是本发明求解发射波形和接收滤波器权值图。
图2是本发明MIMO雷达***输出端的SINR图。
图3是本发明在不同角度上的方向图响应曲线。
具体实施方式
下面结合附图和实施例对本发明进一步说明。
本发明主要构建拉格朗日函数和定义两种神经元,即,变量神经元和拉格朗日神经元。变量神经元的作用是为了找到拉格朗日函数的最优点,也就是满足KKT条件的平衡点,而拉格朗日神经元的作用是驱动轨迹进入可行域,通过拉格朗日函数以及两种神经元的交替作用从而确定参数,得出理想波形和接收滤波器权值,从而最大化输出端的SINR比,改善目标探测性能。
本发明技术方案的详细步骤如下:
步骤一、集中式窄带MIMO雷达***由NT个发射天线和NR个接收天线组成,每一部发射机对应一路发射天线,每个发射机独立地通过全方位发射的方式发射波形,将NT个天线发射的波形记为s(n),n=1,2,...,N,N为每一路波形采样后的长度,到达位于方位角θ处的目标的信号为:
at(θ)Ts(n) (1)
其中at(θ)为发射阵列导向矢量;
假设目标的潜在区域中共有K+1个散射点,其中方位角θ0方向的一个散射点为目标,方位角θk方向的散射点为干扰,其中k=1,2,...,K,且θ0≠θk,则雷达接收机中的信号为所有散射点回波之和,即
其中,x(n)是雷达接收机接收的信号,α0和αk分别为方位角θ0方向的目标的幅度和第k个干扰源的幅度,为接收阵列的导向矢量,为发射阵列导向矢量,d为两天线之间的距离,λ是波长,v(n)为高斯噪声,服从均值为0,方差为
令x=[x(1)T,...,x(N)T]T,s=[s(1)T,...,s(N)T]T,v=[v(1)T,...,v(N)T]T,则式(2)写为:
步骤二、在雷达目标检测中,SINR的大小对***检测概率有很大的影响,因此采用最大化***输出SINR作为雷达波形设计的准则,将接收机采用的权值记作ω,则经过接收机处理后的雷达信号输出r为:
其中,H为共轭转置;
故***输出信噪比模型为:
由于SINR同时依赖于接收机权值ω以及发射波形s,故为了使SINR最大,联合优化发射波形和接收机权值,故优化问题为:
将式(6)改写为如下形式:
步骤三、因为LPNN算法适用于求解实值优化问题,故需要将复数形式的优化问题转化为实值形式的优化问题,因为|ωHA(θm)s|2=|sHAH(θm)ω|2,m=0,1,...,K,将公式(7)中的复数实数化等价为:
其中,λi为拉格朗日乘子,i=1,...,NTN,C0是增广拉格朗日常数项,C0大于0;
式(8)的求解过程是利用式(14)、(15)和(16)循环迭代更新求解,直至公式(14)、(15)和(16)前后两次迭代值的差值小于10-8,则停止循环迭代。
本发明的流程图如图1所示,算法流程如下:
2.for t=1:T;
3.利用式(11)、(12)、(13)求解暂态行为;
5.end for;
6.输出接收端的beampattern图。
本发明实验由8个发射天线和8个接收天线组成的集中式窄带MIMO雷达***,阵元间隔为半波长,每一路波形采样后的长度为16。位于0°处的目标功率为20dB,噪声功率为0dB,位于-50°、-10°和10°处的三个干扰的功率为30dB。图1描述的是LPNN算法实现的模拟形式,为了观察神经网络的动态行为,实施例以图1描述的LPNN算法的离散方式求解发射波形和接收滤波器的权值及其拉格朗日乘子,即
其中,k表示迭代次数,ρ表示步长。
通过式(17)、(18)和(19)的迭代计算直到SINR前后两次迭代值的差值小于10-8时算法停止。从图2中看到以本发明的方法求得的SINR达到57.05dB,从图3中看到输出端在不同角度上的方向响应曲线图在三个干扰处形成很深的零陷,因此可以充分抑制干扰,并且可以准确定位出目标所在方向。
Claims (1)
1.一种MIMO雷达恒模波形和接收机的联合设计方法,其特征在于包括下述步骤:
步骤一、集中式窄带MIMO雷达***由NT个发射天线和NR个接收天线组成,每一部发射机对应一路发射天线,每个发射机独立地通过全方位发射的方式发射波形,将NT个天线发射的波形记为s(n),n=1,2,...,N,N为每一路波形采样后的长度,到达位于方位角θ处的目标的信号为:
at(θ)Ts(n) (1)
其中at(θ)为发射阵列导向矢量;
假设目标的潜在区域中共有K+1个散射点,其中方位角θ0方向的一个散射点为目标,方位角θk方向的散射点为干扰,其中k=1,2,...,K,且θ0≠θk,则雷达接收机中的信号为所有散射点回波之和,即
其中,x(n)是雷达接收机接收的信号,α0和αk分别为方位角θ0方向的目标的幅度和第k个干扰源的幅度,为接收阵列的导向矢量,为发射阵列导向矢量,d为两天线之间的距离,λ是波长,v(n)为高斯噪声,服从均值为0,方差为
令x=[x(1)T,...,x(N)T]T,s=[s(1)T,...,s(N)T]T,v=[v(1)T,...,v(N)T]T,则式(2)写为:
步骤二、采用最大化***输出SINR作为雷达波形设计的准则,将接收机采用的权值记作ω,则经过接收机处理后的雷达信号输出r为:
其中,H为共轭转置;
故***输出信噪比模型为:
由于SINR同时依赖于接收机权值ω以及发射波形s,故为了使SINR最大,联合优化发射波形和接收机权值,故优化问题为:
将式(6)改写为如下形式:
步骤三、将复数形式的优化问题转化为实值形式的优化问题,因为|ωHA(θm)s|2=|sHAH(θm)ω|2,m=0,1,...,K,将公式(7)中的复数实数化等价为:
其中,λi,i=1,...,NTN为拉格朗日乘子,C0是增广拉格朗日常数项,C0大于0;
式(8)的求解过程为使用公式(14)、(15)和(16)循环迭代更新求解,直至公式(14)、(15)和(16)前后两次迭代值的差值小于10-8,则停止循环迭代。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710878926.4A CN107703489B (zh) | 2017-09-26 | 2017-09-26 | 一种mimo雷达恒模波形和接收机的联合设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710878926.4A CN107703489B (zh) | 2017-09-26 | 2017-09-26 | 一种mimo雷达恒模波形和接收机的联合设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107703489A CN107703489A (zh) | 2018-02-16 |
CN107703489B true CN107703489B (zh) | 2021-05-18 |
Family
ID=61174914
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710878926.4A Expired - Fee Related CN107703489B (zh) | 2017-09-26 | 2017-09-26 | 一种mimo雷达恒模波形和接收机的联合设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107703489B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110208770B (zh) * | 2019-07-09 | 2023-06-23 | 中国人民解放军空军工程大学 | 一种基于深层神经网络的雷达波形设计方法 |
CN111142080B (zh) * | 2019-12-28 | 2023-05-30 | 西安电子科技大学 | 一种针对干扰抑制的快速恒模mimo雷达跟踪波形合成方法 |
CN111060902B (zh) * | 2019-12-30 | 2022-10-14 | 电子科技大学 | 一种基于深度学习的mimo雷达波形设计方法 |
CN111337918A (zh) * | 2020-02-17 | 2020-06-26 | 南京航空航天大学 | 一种基于神经网络的机载雷达射频隐身波形选择方法 |
CN112558021B (zh) * | 2020-11-18 | 2022-03-04 | 中国人民解放军国防科技大学 | 用于脉冲压缩雷达***的波形及滤波器设计方法 |
CN113242069B (zh) * | 2021-05-10 | 2022-06-17 | 东南大学 | 一种基于神经网络的码本设计方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101039290A (zh) * | 2007-04-19 | 2007-09-19 | 西安电子科技大学 | 基于自适应训练序列的mimo相关信道估计方法 |
EP2685219A1 (en) * | 2012-07-09 | 2014-01-15 | ST-Ericsson SA | Radiofrequency signal setting reaching a condition of jamming or clipping |
CN104618291A (zh) * | 2014-12-31 | 2015-05-13 | 重庆邮电大学 | 基于余量自适应准则的资源分配方法及装置 |
CN105319545A (zh) * | 2015-11-09 | 2016-02-10 | 大连大学 | 提高stap检测性能的mimo-ofdm雷达波形设计方法 |
CN106501778A (zh) * | 2016-05-16 | 2017-03-15 | 南京航空航天大学 | 一种基于射频隐身性能优化的dmrs稳健波形设计方法 |
CN106680776A (zh) * | 2016-12-13 | 2017-05-17 | 西北工业大学 | 对多普勒信息不敏感的低旁瓣波形设计方法 |
CN106842113A (zh) * | 2016-12-12 | 2017-06-13 | 西北工业大学 | 高采样1比特量化情况下的信号到达角高精度估计方法 |
-
2017
- 2017-09-26 CN CN201710878926.4A patent/CN107703489B/zh not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101039290A (zh) * | 2007-04-19 | 2007-09-19 | 西安电子科技大学 | 基于自适应训练序列的mimo相关信道估计方法 |
EP2685219A1 (en) * | 2012-07-09 | 2014-01-15 | ST-Ericsson SA | Radiofrequency signal setting reaching a condition of jamming or clipping |
CN104618291A (zh) * | 2014-12-31 | 2015-05-13 | 重庆邮电大学 | 基于余量自适应准则的资源分配方法及装置 |
CN105319545A (zh) * | 2015-11-09 | 2016-02-10 | 大连大学 | 提高stap检测性能的mimo-ofdm雷达波形设计方法 |
CN106501778A (zh) * | 2016-05-16 | 2017-03-15 | 南京航空航天大学 | 一种基于射频隐身性能优化的dmrs稳健波形设计方法 |
CN106842113A (zh) * | 2016-12-12 | 2017-06-13 | 西北工业大学 | 高采样1比特量化情况下的信号到达角高精度估计方法 |
CN106680776A (zh) * | 2016-12-13 | 2017-05-17 | 西北工业大学 | 对多普勒信息不敏感的低旁瓣波形设计方法 |
Non-Patent Citations (4)
Title |
---|
Adaptive Beamforming for Uniform Linear Arrays With Unknown Mutual Coupling;Bin Liao;《IEEE ANTENNAS AND WIRELESS PROPAGATION LETTERS》;20121231;第11卷;464-467 * |
Waveform Design With Unit Modulus and Spectral Shape Constraints via Lagrange Programming Neural Network;Junli Liang;《IEEE JOURNAL OF SELECTED TOPICS IN SIGNAL PROCESSING》;20151231;第9卷(第8期);1377-1386 * |
基于随机LPNN网络的优化设计研究;周盛强;《机械科学与技术》;20071231;第26卷(第3期);309-402 * |
杂波条件下稳健的 MIMO 雷达发射波形和接收滤波器联合优化设计;王玉玺;《电子 与信息学报》;20170831;第39卷(第8期);1879-1886 * |
Also Published As
Publication number | Publication date |
---|---|
CN107703489A (zh) | 2018-02-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107703489B (zh) | 一种mimo雷达恒模波形和接收机的联合设计方法 | |
CN101369014B (zh) | 应用于多输入多输出雷达的双边约束自适应波束形成方法 | |
US9413474B2 (en) | Efficient large-scale multiple input multiple output communications | |
CN102830387B (zh) | 一种基于数据预处理的协方差矩阵正交化波束形成方法 | |
CN106021637B (zh) | 互质阵列中基于迭代稀疏重构的doa估计方法 | |
CN109407055B (zh) | 基于多径利用的波束形成方法 | |
CN105319545B (zh) | 提高stap检测性能的mimo-ofdm雷达波形设计方法 | |
CN103605122A (zh) | 相干mimo雷达的收发稳健降维自适应波束形成方法 | |
Wang et al. | Enhanced automotive sensing assisted by joint communication and cognitive sparse MIMO radar | |
KR101603630B1 (ko) | 송신신호 파형의 직교성을 이용한 재머 신호 억제 및 원신호 도래각 추정 방법 및 장치 | |
CN106646388B (zh) | 基于嵌套阵列的mimo雷达抗干扰方法 | |
Varade et al. | Robust algorithms for DOA estimation and adaptive beamforming for smart antenna application | |
CN103837861A (zh) | 基于特征子空间的子阵级线性约束自适应波束形成方法 | |
CN104833959A (zh) | 基于目标先验信息的mimo雷达波形优化方法 | |
CN105487054B (zh) | 提高基于mimo-ofdm雷达stap最差检测性能的稳健波形设计方法 | |
CN105807275A (zh) | 基于部分杂波先验知识的mimo-ofdm-stap稳健波形设计 | |
CN105049382A (zh) | 一种抗期望信号导向矢量失配的零陷展宽自适应天线波束成形方法 | |
Cheng et al. | Double-phase-shifter based hybrid beamforming for mmWave DFRC in the presence of extended target and clutters | |
CN109861933B (zh) | 一种基于music算法和预编码的毫米波mimo信道估计方法 | |
CN104360316A (zh) | 一种基于协方差矩阵锥化的阵列天线自适应波束形成方法 | |
CN106646420B (zh) | 基于lfm信号的mimo雷达发射方向图设计方法 | |
Sun et al. | A robust secure hybrid analog and digital receive beamforming scheme for efficient interference reduction | |
JP2012112712A (ja) | 送受信装置 | |
CN107290732A (zh) | 一种量子大***的单基地mimo雷达测向方法 | |
CN114726686B (zh) | 一种均匀面阵毫米波大规模mimo信道估计方法 |
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 |
Granted publication date: 20210518 Termination date: 20210926 |
|
CF01 | Termination of patent right due to non-payment of annual fee |