CN116466299A - 二维子阵级稀疏阵列fpmimo雷达收发波束合成方法 - Google Patents
二维子阵级稀疏阵列fpmimo雷达收发波束合成方法 Download PDFInfo
- Publication number
- CN116466299A CN116466299A CN202310729386.9A CN202310729386A CN116466299A CN 116466299 A CN116466299 A CN 116466299A CN 202310729386 A CN202310729386 A CN 202310729386A CN 116466299 A CN116466299 A CN 116466299A
- Authority
- CN
- China
- Prior art keywords
- array
- signal
- subarray
- vector
- fpmimo
- 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.)
- Granted
Links
- 238000001308 synthesis method Methods 0.000 title claims abstract description 9
- 239000013598 vector Substances 0.000 claims abstract description 78
- 239000011159 matrix material Substances 0.000 claims abstract description 38
- 238000005070 sampling Methods 0.000 claims abstract description 17
- 238000001228 spectrum Methods 0.000 claims abstract description 12
- 238000000034 method Methods 0.000 claims description 22
- 238000012545 processing Methods 0.000 claims description 21
- 230000005540 biological transmission Effects 0.000 claims description 8
- 230000001427 coherent effect Effects 0.000 claims description 7
- 238000000354 decomposition reaction Methods 0.000 claims description 6
- 230000008054 signal transmission Effects 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000002592 echocardiography Methods 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- 230000002452 interceptive effect Effects 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims 2
- 238000003786 synthesis reaction Methods 0.000 claims 2
- 238000004422 calculation algorithm Methods 0.000 description 6
- 238000004088 simulation Methods 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000001914 filtration Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000003491 array Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002401 inhibitory effect Effects 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Abstract
本发明公开了一种二维子阵级稀疏阵列FPMIMO雷达收发波束合成方法:步骤1:子阵以相控阵方式发射LFM信号;步骤2:获取虚拟扩展阵列样本数据;步骤3:计算协方差矩阵;步骤4,计算虚拟扩展阵列导向矢量;步骤5:获取采样信号空间谱估计;步骤6:计算干扰信号几何中心角度;步骤7:子阵以几何中心角度发射信号,构建虚拟扩展阵列的INCM;步骤8:重构虚拟扩展阵列的信号协方差矩阵;步骤9:估计信号导向矢量;步骤10:计算雷达收发波束形成权矢量;步骤11:计算雷达收发波束形成输出数据。本发明解决了在增强期望目标回波信号的同时无法有效抑制干扰造成的接收信号SINR较低的问题。
Description
技术领域
本发明属于雷达天线阵列数字波束形成设计技术领域,具体涉及一种二维子阵级稀疏阵列FPMIMO雷达收发波束合成方法。
背景技术
目前,基于子阵的阵列天线波束形成算法研究以子阵划分方式研究为主,由于目前的HPMIMO(Hybrid Phased-Multiple-Input-Multiple-Output,混合相控阵-多输入多输出)雷达仅考虑了发射天线的相控阵处理,而没有考虑接收信号时的相干处理增益,其输出信号的SINR(信干噪比)低于FPMIMO雷达,更没有结合二维子阵级稀疏阵列进行研究。现阶段对二维子阵级稀疏阵列的FPMIMO雷达收发波束形成算法的研究成果极少,并且在算法稳健性和稀疏阵列FPMIMO雷达稳健收发波束形成算法的研究方面还存在明显不足,显示出在增强期望目标回波信号的同时无法对干扰信号进行有效抑制造成接收信号的SINR较低的问题,无法满足已经产生的相应的工程应用背景需求。
发明内容
本发明的目的在于,提供一种二维子阵级稀疏阵列FPMIMO雷达收发波束合成方法,以解决现有算法稳健性不足,在增强期望目标回波信号的同时无法有效抑制干扰造成的接收信号SINR较低的问题。
为了实现上述目的,本发明采用如下技术方案予以解决:
一种二维子阵级稀疏阵列FPMIMO雷达收发波束合成方法,包括如下步骤:
步骤1:二维子阵级稀疏阵列FPMIMO雷达中的所有子阵在发射信号权矢量的控制下,以相控阵方式向目标方向发射相互正交的LFM信号;二维子阵级稀疏阵列FPMIMO雷达中的子阵个数为K;
步骤2:每个子阵分别以相控阵方式接收回波得到一路数字信号,共得到K路数字信号;再将每一路数字信号分别进行K路MF处理,共得到K²路数据,作为虚拟扩展阵列样本数据;
步骤3:采用下式计算虚拟扩展阵列样本数据的协方差矩阵:
式中:
—虚拟扩展阵列样本数据的协方差矩阵;
L—回波的采样快拍数;
—虚拟扩展阵列样本数据;
—共轭转置矩阵;
步骤4,计算虚拟扩展阵列的导向矢量:
式中:
—虚拟扩展阵列的导向矢量;
—信号发射的相干处理增益矢量;
—波形分集矢量;
—接收信号的相干处理增益矢量;
φ—方位角;
θ—俯仰角;
w1、w2、…、wk—第1~K个子阵的均匀加权的权矢量,w1= w2=…=wk=;
—共轭转置矩阵;
—转置矩阵;
—第1~K个子阵的发射信号导向矢量,
;
vec— vec函数,用于将一个矩阵向量化处理;
j—虚数单位;
—指数函数;
f c—发射信号载波频率;
d—子阵内相邻阵元间距;
c—光速;
x 1、x 2、…、x k—第1~K个子阵的位置x坐标;
y 1、y 2、…、y k—第1~K个子阵的位置y坐标;子阵的位置定义为该子阵左下角阵元所在的位置;第1个子阵的位置(x 1,y 1) =(0,0);
u1、u2、…、uK—第1~K个子阵的子阵内部的信号接收权矢量,u1=u2=…=uK=;
n 1、n 2—子阵的行数和列数;
步骤5:通过下式获取采样信号的空间谱估计:
式中:
—采样信号的空间谱估计;
φ—方位角;
θ—俯仰角;
—虚拟扩展阵列的导向矢量;
—共轭转置矩阵;
—逆矩阵;
步骤6:对步骤5获得的采样信号的空间谱估计进行峰值提取,得到S+1个回波信号方向:/>,其中,为目标方向,/>为S个干扰信号的来波方向,由下式计算干扰信号的几何中心角度:
式中:
—干扰信号的几何中心角度;
i—第i个干扰信号的来波方向;
—第i个干扰信号的来波的方位角;
—第i个干扰信号的来波的俯仰角;
步骤7:所有子阵均以干扰信号的几何中心角度为波束指向发射信号,构建虚拟扩展阵列的INCM/>;
步骤8:采用下式重构虚拟扩展阵列的信号协方差矩阵:
式中:
—虚拟扩展阵列的信号协方差矩阵;
—虚拟扩展阵列的导向矢量;
—方位向的角度偏差;
—俯仰向的角度偏差;
—目标的方位角;
—目标的俯仰角;
步骤9:采用下式对虚拟扩展阵列的信号协方差矩阵进行特征值分解,并估计信号导向矢量vd;
vd=vmax
式中:
V—由特征向量构成的矩阵;
D—由特征值构成的对角矩阵;
—特征值分解操作;
vd—信号导向矢量;
vmax—最大特征值对应的特征向量;
步骤10:采用下式计算FPMIMO雷达收发波束形成权矢量wtr:
式中:
wtr—雷达收发波束形成权矢量;
—虚拟扩展阵列采样数据的INCM;
vd—信号导向矢量;
—共轭转置矩阵;
步骤11:采用下式计算FPMIMO雷达收发波束形成输出数据yout:
式中:
yout—雷达收发波束形成输出数据;
—虚拟扩展阵列样本数据。
进一步的,步骤1中,K为6~30。
进一步的,方位角φ的取值范围为,俯仰角θ的取值范围为。
进一步的,步骤7包括如下子步骤:
步骤71,调整所有子阵的发射信号波束指向为进行信号发射,采用下式计算子阵的发射波束形成权矢量wt,每个子阵在发射波束形成权矢量wt的控制下以相控阵方式向/>方向发射相互正交的LFM信号;
式中:
wt—子阵的发射波束形成权矢量;
C—约束向量构成的矩阵;
—子阵k在方向/>上的导向矢量;k∈[1,K];
—目标方向;
—目标的方位角;
—目标的俯仰角;
—子阵k在目标方向/>上的导向矢量;
—共轭转置矩阵;
步骤72,每个子阵分别以相控阵方式接收回波,获得干扰加噪声的采样信号;将干扰加噪声的采样信号转换为K路数字信号,然后将每路数字信号分别通过K路MF处理得到估计的虚拟扩展阵列干扰加噪声信号;
步骤73,通过下式计算得到虚拟扩展阵列采样数据的INCM;
式中:
—虚拟扩展阵列采样数据的INCM;
L—回波的采样快拍数;
—虚拟扩展阵列干扰加噪声信号;
—共轭转置矩阵。
进一步的,步骤8中,方位向的角度偏差取3~5°,俯仰向的角度偏差/>取3~5°。
相较于现有技术,本发明具有如下技术效果:
1、通过在子阵内部采用相控阵收发信号模式,增强了在目标方向信号的发送和接收能力,同时通过匹配滤波(MF处理)扩展了信号自由度,提高了干扰抑制性能。
2、通过估计虚拟扩展阵列的导向矢量,提高了收发波束形成的指向精度,增强了对期望信号的接收能力,并通过发射指向干扰信号方向的探测信号,获取虚拟扩展阵列的INCM,提高了对干扰信号的抑制能力,从而有效解决了算法稳健性不足的问题。
综上,本发明的方法中,子阵内部以相控阵方式发射和接收信号,子阵之间工作在MIMO模式,利用MF处理技术对子阵采样数据进行虚拟扩展,通过虚拟阵列导向矢量估计和虚拟扩展阵列INCM的获取,得到FPMIMO雷达子阵的收发波束形成权矢量,从而利用虚拟扩展数据进行空域滤波处理,在增强期望目标回波信号的同时对干扰信号进行有效抑制,提高了接收信号的SINR。
附图说明
图1是二维子阵级稀疏阵列FPMIMO雷达信号发射处理框图;
图2是二维子阵级稀疏阵列FPMIMO雷达信号接收处理框图;
图3是虚拟扩展阵列的INCM获取时的发射波束示意图;
图4是虚拟扩展阵列的空间谱估计;
图5是子阵发射信号波束图;
图6是子阵发射信号波束指向剖面图;
图7是子阵发射信号目标方向剖面图;
图8是二维子阵级稀疏阵列FPMIMO雷达收发波束图;
图9是期望信号方位向波束剖面图;
图10是期望信号俯仰向波束剖面图。
以下结合附图和具体实施方式对本发明进一步解释说明。
具体实施方式
本发明的方法涉及的技术术语如下:
(1)二维子阵级稀疏阵列:在二维平面上稀疏布置的多个相同大小的子阵形成的阵列,子阵的行数和列数记为n1、n2,各子阵内部相邻阵元间距d=λ/2,λ为发射信号的波长;稀疏布置是指二维子阵级稀疏阵列中的任意两个子阵之间不重叠。
(2)FPMIMO雷达:Full Phased-Multiple-Input-Multiple-Output Radar,即全相控阵-多输入多输出雷达。
(3)LFM信号:Linear Frequency Modulation,即线性调频信号。
(4)INCM:Interference-plus-Noise Covariance Matrix,即干扰加噪声协方差矩阵。
(5)MF处理:Matched Filtering,即匹配滤波。
本发明涉及的二维子阵级稀疏阵列FPMIMO雷达的信号发射和接收的原理如下:
二维子阵级稀疏阵列FPMIMO雷达采用收/发共置阵列,即雷达信号的发射和接收使用同一阵列,阵列天线工作在FPMIMO模式,各个子阵内部工作在PAR模式。
如图1所示,二维子阵级稀疏阵列FPMIMO雷达信号发射处理过程:首先由发射信号处理机生成K个相互正交的LFM信号和相应的子阵发射信号权矢量,LFM信号分别经对应的发射通道进行D/A转换、上变频和功率放大处理后,由各个子阵(子阵1~子阵K)在发射信号权矢量的控制下以相控阵方式向目标方向发射出去,在收发波束形成中,所有子阵的发射信号波束指向同一方向。
如图2所示,二维子阵级稀疏阵列FPMIMO雷达信号接收处理过程:当阵列天线接收回波信号时,K个子阵分别以相控阵模式形成接收波束,将接收信号送入对应接收通道进行低噪声放大、带通滤波、混频器、下变频处理、A/D转换后变为数字信号,然后将每路数字信号送入K路MF处理,共得到K2路数字信号。接收信号处理机提供子阵内部以相控阵模式接收的权矢量。
基于上述原理,本发明给出的二维子阵级稀疏阵列FPMIMO雷达收发波束合成方法,包括如下步骤:
步骤1:二维子阵级稀疏阵列FPMIMO雷达中的所有子阵在发射信号权矢量的控制下,以相控阵方式向目标方向发射相互正交的LFM信号;二维子阵级稀疏阵列FPMIMO雷达中的子阵个数为K,K为6~30。
步骤2:每个子阵分别以相控阵方式接收回波得到一路数字信号,共得到K路数字信号;再将每一路数字信号分别进行K路MF处理,共得到K²路数据,作为虚拟扩展阵列样本数据。该步骤实现了将K路数字信号扩展为K²维信号。
步骤3:采用下式计算虚拟扩展阵列样本数据的协方差矩阵:
,
式中:
—虚拟扩展阵列样本数据的协方差矩阵;
L—回波的采样快拍数;
—虚拟扩展阵列样本数据,由步骤2得到;
—共轭转置矩阵。
步骤4,计算虚拟扩展阵列的导向矢量:
式中:
—虚拟扩展阵列的导向矢量;
—信号发射的相干处理增益矢量;
—波形分集矢量;
—接收信号的相干处理增益矢量;
φ—方位角,取值范围为;
θ—俯仰角,取值范围为;
w1、w2、…、wk—第1~K个子阵的均匀加权的权矢量,w1= w2=…=wk=;
—共轭转置矩阵;
—转置矩阵;
—第1~K个子阵的发射信号导向矢量,
;
vec— vec函数,用于将一个矩阵向量化处理;
j—虚数单位;
—指数函数;
f c—发射信号载波频率;
d—子阵内相邻阵元间距;
c—光速;
x 1、x 2、…、x k—第1~K个子阵的位置x坐标;
y 1、y 2、…、y k—第1~K个子阵的位置y坐标;子阵的位置定义为该子阵左下角的阵元所在的位置;第1个子阵的位置(x 1,y 1) =(0,0);
u1、u2、…、uK—第1~K个子阵的子阵内部的信号接收权矢量,u1=u2=…=uK=;
n 1、n 2—子阵的行数和列数;
步骤5:通过下式获取采样信号的空间谱估计:
式中:
—采样信号的空间谱估计;
φ—方位角,取值范围为;
θ—俯仰角,取值范围为;
—虚拟扩展阵列的导向矢量,由步骤4得到;
—共轭转置矩阵;
—逆矩阵。
步骤6:对步骤5获得的采样信号的空间谱估计进行峰值提取,得到S+1个回波信号方向:/>,其中,为目标方向,/>为S个干扰信号的来波方向,由下式计算干扰信号的几何中心角度;
式中:
—干扰信号的几何中心角度;
i—第i个干扰信号的来波方向;
—第i个干扰信号的来波的方位角;
—第i个干扰信号的来波的俯仰角。
步骤7:所有子阵均以为波束指向发射信号,构建虚拟扩展阵列的INCM/>;具体操作如下:
步骤71,调整所有子阵的发射信号波束指向为进行信号发射,采用下式计算子阵的发射波束形成权矢量wt,每个子阵在发射波束形成权矢量wt的控制下以相控阵方式向/>方向发射相互正交的LFM信号;
式中:
wt—子阵的发射波束形成权矢量;
C—约束向量构成的矩阵;
—子阵/>在方向/>上的导向矢量;k∈[1,K];
—目标方向,由步骤6得到;
—目标的方位角;
—目标的俯仰角;
—子阵k在目标方向/>上的导向矢量;
—共轭转置矩阵;
步骤72,每个子阵分别以相控阵方式接收回波,获得干扰加噪声的采样信号;将干扰加噪声的采样信号转换为K路数字信号,然后将每路数字信号分别通过K路MF处理得到估计的虚拟扩展阵列干扰加噪声信号。/>为K2维信号。
步骤73,通过下式计算得到虚拟扩展阵列采样数据的INCM;
式中:
—虚拟扩展阵列采样数据的INCM;
L—回波的采样快拍数;
—虚拟扩展阵列干扰加噪声信号,由步骤72得到;
—共轭转置矩阵。
步骤8:采用下式重构虚拟扩展阵列的信号协方差矩阵:
式中:
—虚拟扩展阵列的信号协方差矩阵;
—虚拟扩展阵列的导向矢量,由步骤4得到;
—方位向的角度偏差,一般取3~5°;
—俯仰向的角度偏差,一般取3~5°;
—目标的方位角;
—目标的俯仰角;
步骤9:采用下式对虚拟扩展阵列的信号协方差矩阵进行特征值分解,并估计信号导向矢量vd;
vd=vmax
式中:
V—由特征向量构成的矩阵;
D—由特征值构成的对角矩阵;
—特征值分解操作;
vd—信号导向矢量;
vmax—最大特征值对应的特征向量。
步骤10:采用下式计算FPMIMO雷达收发波束形成权矢量wtr:
式中:
wtr—雷达收发波束形成权矢量;
—虚拟扩展阵列采样数据的INCM;
vd—信号导向矢量;
—共轭转置矩阵。
步骤11:采用下式计算FPMIMO雷达收发波束形成输出数据yout:
式中:
yout—雷达收发波束形成输出数据;
—虚拟扩展阵列样本数据,由步骤2得到。
为了验证本发明的可行和有效性,本发明进行了基于二维子阵级稀疏阵列FPMIMO雷达的稳健收发波束合成方法的性能仿真验证。
由发射信号处理机生成子阵之间发射的相互正交的LFM信号,脉宽为Tp=10μs,频率间隔为fm=3MHz,带宽为3MHz,复采样频率为f s=40MHz,每条回波的采样波门宽度为512。本发明在重构虚拟阵列的信号加噪声协方差矩阵 (Signal-plus-Noise Covariance Matrix,SNCM)时,取方位角区间为,俯仰角区间为/>,获取虚拟扩展阵列的INCM。
注意,根据实际应用背景需求,在仿真过程中设置方位角的变化范围为,俯仰角的变化范围为/>。假设已知目标真实的方向为,期望目标回波的输入SNR为-10dB,在探测范围内存在2种4个干扰,其中第1种是2个干扰目标(如转发式干扰、无源干扰),其角度分别为/>和/>,干扰信号的输入INR均为60dB,第2种是2个主动式有源干扰,其空间位置分别是和/>方向,主动干扰信号为大功率噪声干扰,其输入INR均为80dB。
仿真实验:波束图仿真
采样信号通过MF处理进行虚拟扩展,可得到如图4所示的空间谱,由此可得期望目标回波信号位于角度,4个干扰信号来波方向分别为/>,,/>和/>,可确定在获取干扰加噪声采样信号时,发射信号方向为/>,发射波束零陷方向为/>,子阵的发射信号方向图如图5所示。为便于观察,截取发射波束图的剖面图如图6、图7所示,其中图6中表示波束指向为/>,图7表示目标所在方向/>是发射信号的零陷位置,可避免期望目标回波信号混入采样信号中,从而高精度地获取干扰信号和噪声。
波束图仿真结果如图8所示,由仿真结果可知,稀疏阵列的栅瓣得到有效抑制,旁瓣增益降低,波束指向精确地位于期望信号方向,证明了本发明的方法的有效性。其中,在期望信号方向的波束剖面图如图9、图10所示,波束指向精确地位于角度处,图9显示方位向的PSLL为-18.12dB,图10显示在俯仰向的PSLL为-16.13dB。/>
Claims (5)
1.一种二维子阵级稀疏阵列FPMIMO雷达收发波束合成方法,其特征在于,包括如下步骤:
步骤1:二维子阵级稀疏阵列FPMIMO雷达中的所有子阵在发射信号权矢量的控制下,以相控阵方式向目标方向发射相互正交的LFM信号;二维子阵级稀疏阵列FPMIMO雷达中的子阵个数为K;
步骤2:每个子阵分别以相控阵方式接收回波得到一路数字信号,共得到K路数字信号;再将每一路数字信号分别进行K路MF处理,共得到K²路数据,作为虚拟扩展阵列样本数据;
步骤3:采用下式计算虚拟扩展阵列样本数据的协方差矩阵:
式中:
—虚拟扩展阵列样本数据的协方差矩阵;
L—回波的采样快拍数;
—虚拟扩展阵列样本数据;
—共轭转置矩阵;
步骤4,计算虚拟扩展阵列的导向矢量:
式中:
—虚拟扩展阵列的导向矢量;
—信号发射的相干处理增益矢量;
—波形分集矢量;
—接收信号的相干处理增益矢量;
φ—方位角;
θ—俯仰角;
w1、w2、…、wk —第1~K个子阵的均匀加权的权矢量,w1= w2=…=wk=;
—共轭转置矩阵;
—转置矩阵;
—第1~K个子阵的发射信号导向矢量,
;
vec— vec函数,用于将一个矩阵向量化处理;
j—虚数单位;
—指数函数;
f c—发射信号载波频率;
d—子阵内相邻阵元间距;
c—光速;
x 1、x 2、…、x k—第1~K个子阵的位置x坐标;
y 1、y 2、…、y k—第1~K个子阵的位置y坐标;子阵的位置定义为该子阵左下角阵元所在的位置;第1个子阵的位置(x 1, y 1) =(0,0);
u1、u2、…、uK—第1~K个子阵的子阵内部的信号接收权矢量,u1=u2=…=uK =;
n 1、n 2—子阵的行数和列数;
步骤5:通过下式获取采样信号的空间谱估计:
式中:
—采样信号的空间谱估计;
φ—方位角;
θ—俯仰角;
—虚拟扩展阵列的导向矢量;
—共轭转置矩阵;
—逆矩阵;
步骤6:对步骤5获得的采样信号的空间谱估计进行峰值提取,得到S+1个回波信号方向:/>,其中,为目标方向,/>为S个干扰信号的来波方向,由下式计算干扰信号的几何中心角度:
式中:
—干扰信号的几何中心角度;
i—第i个干扰信号的来波方向;
—第i个干扰信号的来波的方位角;
—第i个干扰信号的来波的俯仰角;
步骤7:所有子阵均以干扰信号的几何中心角度为波束指向发射信号,构建虚拟扩展阵列的INCM />;
步骤8:采用下式重构虚拟扩展阵列的信号协方差矩阵:
式中:
—虚拟扩展阵列的信号协方差矩阵;
—虚拟扩展阵列的导向矢量;
—方位向的角度偏差;
—俯仰向的角度偏差;
—目标的方位角;
—目标的俯仰角;
步骤9:采用下式对虚拟扩展阵列的信号协方差矩阵进行特征值分解,并估计信号导向矢量vd;
vd =vmax
式中:
V—由特征向量构成的矩阵;
D—由特征值构成的对角矩阵;
—特征值分解操作;
vd—信号导向矢量;
vmax—最大特征值对应的特征向量;
步骤10:采用下式计算FPMIMO雷达收发波束形成权矢量wtr:
式中:
wtr —雷达收发波束形成权矢量;
—虚拟扩展阵列采样数据的INCM;
vd—信号导向矢量;
—共轭转置矩阵;
步骤11:采用下式计算FPMIMO雷达收发波束形成输出数据yout:
式中:
yout—雷达收发波束形成输出数据;
—虚拟扩展阵列样本数据。
2.如权利要求1所述的二维子阵级稀疏阵列FPMIMO雷达收发波束合成方法,其特征在于,步骤1中,K为6~30。
3.如权利要求1所述的二维子阵级稀疏阵列FPMIMO雷达收发波束合成方法,其特征在于,方位角φ的取值范围为,俯仰角θ的取值范围为/>。
4.如权利要求1所述的二维子阵级稀疏阵列FPMIMO雷达收发波束合成方法,其特征在于,步骤7包括如下子步骤:
步骤71,调整所有子阵的发射信号波束指向为进行信号发射,采用下式计算子阵的发射波束形成权矢量wt,每个子阵在发射波束形成权矢量wt的控制下以相控阵方式向/>方向发射相互正交的LFM信号;
式中:
wt —子阵的发射波束形成权矢量;
C—约束向量构成的矩阵;
—子阵k在方向/>上的导向矢量;k∈[1,K];
—目标方向;
—目标的方位角;
—目标的俯仰角;
—子阵k在目标方向/>上的导向矢量;
—共轭转置矩阵;
步骤72,每个子阵分别以相控阵方式接收回波,获得干扰加噪声的采样信号;将干扰加噪声的采样信号转换为K路数字信号,然后将每路数字信号分别通过K路MF处理得到估计的虚拟扩展阵列干扰加噪声信号;
步骤73,通过下式计算得到虚拟扩展阵列采样数据的INCM;
式中:
—虚拟扩展阵列采样数据的INCM;
L—回波的采样快拍数;
—虚拟扩展阵列干扰加噪声信号;
—共轭转置矩阵。
5.如权利要求1所述的二维子阵级稀疏阵列FPMIMO雷达收发波束合成方法,其特征在于,步骤8中,方位向的角度偏差取3~5°,俯仰向的角度偏差/>取3~5°。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310729386.9A CN116466299B (zh) | 2023-06-20 | 2023-06-20 | 二维子阵级稀疏阵列fpmimo雷达收发波束合成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310729386.9A CN116466299B (zh) | 2023-06-20 | 2023-06-20 | 二维子阵级稀疏阵列fpmimo雷达收发波束合成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116466299A true CN116466299A (zh) | 2023-07-21 |
CN116466299B CN116466299B (zh) | 2023-08-18 |
Family
ID=87179287
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310729386.9A Active CN116466299B (zh) | 2023-06-20 | 2023-06-20 | 二维子阵级稀疏阵列fpmimo雷达收发波束合成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116466299B (zh) |
Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07115317A (ja) * | 1993-10-18 | 1995-05-02 | Mitsubishi Heavy Ind Ltd | 電子走査アンテナ |
JP2000180538A (ja) * | 1998-12-14 | 2000-06-30 | Denso Corp | 走査型レーダ装置 |
US20160091598A1 (en) * | 2014-09-26 | 2016-03-31 | The Govemment of the United States of America, as represented by the Secretary of the Navy | Sparse Space-Time Adaptive Array Architecture |
CN105699944A (zh) * | 2016-01-23 | 2016-06-22 | 西安电子科技大学 | 基于lfm基波束的mimo雷达部分相关波形设计方法 |
US20170371030A1 (en) * | 2016-06-23 | 2017-12-28 | GM Global Technology Operations LLC | Low latency decoding in multi-input multi-output radar |
US20190212411A1 (en) * | 2018-01-08 | 2019-07-11 | Hyundai Mobis Co., Ltd. | Method and apparatus for estimating direction of arrival using generation of virtual received signals |
CN111948643A (zh) * | 2020-08-13 | 2020-11-17 | 中国科学院空天信息创新研究院 | 基于多波束非线性鉴角的目标跟踪方法 |
WO2021131601A1 (ja) * | 2019-12-25 | 2021-07-01 | 日立Astemo株式会社 | レーダ装置 |
CN113325385A (zh) * | 2021-07-15 | 2021-08-31 | 中国人民解放军火箭军工程大学 | 一种相控阵-mimo雷达模式收发波束形成抗干扰方法 |
CN114218814A (zh) * | 2022-02-23 | 2022-03-22 | 中国人民解放军火箭军工程大学 | 一种降低距离维波束形成副瓣的稀疏阵列优化配置方法 |
US20220120847A1 (en) * | 2020-10-20 | 2022-04-21 | Aptiv Technologies Limited | Radar System with Modified Orthogonal Linear Antenna Subarrays |
CN114399044A (zh) * | 2021-12-14 | 2022-04-26 | 西安电子科技大学 | 一种子阵级稀疏阵列发射波束旁瓣电平优化方法 |
CN114648041A (zh) * | 2022-03-15 | 2022-06-21 | 中国人民解放军火箭军工程大学 | 一种基于平行稀疏阵列的二维欠定doa估计算法 |
CN115085784A (zh) * | 2022-05-27 | 2022-09-20 | 中国人民解放军火箭军工程大学 | 一种分级加权稀疏子阵天线阵列稳健波束合成方法 |
CN116148777A (zh) * | 2022-12-27 | 2023-05-23 | 北京理工大学 | 一种基于阵列虚拟扩展的阵列雷达主瓣干扰抑制方法 |
-
2023
- 2023-06-20 CN CN202310729386.9A patent/CN116466299B/zh active Active
Patent Citations (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07115317A (ja) * | 1993-10-18 | 1995-05-02 | Mitsubishi Heavy Ind Ltd | 電子走査アンテナ |
JP2000180538A (ja) * | 1998-12-14 | 2000-06-30 | Denso Corp | 走査型レーダ装置 |
US20160091598A1 (en) * | 2014-09-26 | 2016-03-31 | The Govemment of the United States of America, as represented by the Secretary of the Navy | Sparse Space-Time Adaptive Array Architecture |
CN105699944A (zh) * | 2016-01-23 | 2016-06-22 | 西安电子科技大学 | 基于lfm基波束的mimo雷达部分相关波形设计方法 |
US20170371030A1 (en) * | 2016-06-23 | 2017-12-28 | GM Global Technology Operations LLC | Low latency decoding in multi-input multi-output radar |
US20190212411A1 (en) * | 2018-01-08 | 2019-07-11 | Hyundai Mobis Co., Ltd. | Method and apparatus for estimating direction of arrival using generation of virtual received signals |
WO2021131601A1 (ja) * | 2019-12-25 | 2021-07-01 | 日立Astemo株式会社 | レーダ装置 |
CN111948643A (zh) * | 2020-08-13 | 2020-11-17 | 中国科学院空天信息创新研究院 | 基于多波束非线性鉴角的目标跟踪方法 |
US20220120847A1 (en) * | 2020-10-20 | 2022-04-21 | Aptiv Technologies Limited | Radar System with Modified Orthogonal Linear Antenna Subarrays |
CN113325385A (zh) * | 2021-07-15 | 2021-08-31 | 中国人民解放军火箭军工程大学 | 一种相控阵-mimo雷达模式收发波束形成抗干扰方法 |
CN114399044A (zh) * | 2021-12-14 | 2022-04-26 | 西安电子科技大学 | 一种子阵级稀疏阵列发射波束旁瓣电平优化方法 |
CN114218814A (zh) * | 2022-02-23 | 2022-03-22 | 中国人民解放军火箭军工程大学 | 一种降低距离维波束形成副瓣的稀疏阵列优化配置方法 |
CN114648041A (zh) * | 2022-03-15 | 2022-06-21 | 中国人民解放军火箭军工程大学 | 一种基于平行稀疏阵列的二维欠定doa估计算法 |
CN115085784A (zh) * | 2022-05-27 | 2022-09-20 | 中国人民解放军火箭军工程大学 | 一种分级加权稀疏子阵天线阵列稳健波束合成方法 |
CN116148777A (zh) * | 2022-12-27 | 2023-05-23 | 北京理工大学 | 一种基于阵列虚拟扩展的阵列雷达主瓣干扰抑制方法 |
Non-Patent Citations (2)
Title |
---|
JIANFENG LI等: "DOA estimation for sparse nested MIMO radar with velocity receive sensor array", MULTIDIM SYST SIGN PROCESS, pages 1397 - 1410 * |
张逸楠等: "基于无人机集群的近场线性稀疏阵列波束形成研究", 电子与信息学报, vol. 45, no. 1, pages 181 - 190 * |
Also Published As
Publication number | Publication date |
---|---|
CN116466299B (zh) | 2023-08-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106788653B (zh) | 一种基于协方差矩阵重构的自适应波束形成方法 | |
CN108508423B (zh) | 基于异型阵的子阵数字和差单脉冲测角方法 | |
CN101369014B (zh) | 应用于多输入多输出雷达的双边约束自适应波束形成方法 | |
CN103728597B (zh) | 基于辅助阵列的网络化雷达抑制压制式主瓣干扰方法 | |
CN109765562B (zh) | 一种三维前视声像声纳系统和方法 | |
CN109375213B (zh) | 基于子阵划分的频率分集阵列信号处理方法 | |
CN113325385B (zh) | 一种相控阵-mimo雷达模式收发波束形成抗干扰方法 | |
CN108710112B (zh) | 空时编码阵列空间平滑解相干谱估计方法 | |
CN113660017A (zh) | 一种irs辅助的双功能雷达通信系统的sinr最大化方法 | |
Qu et al. | Performance analysis of beamforming for MIMO radar | |
CN117250596A (zh) | 一种频率分集-相控阵mimo雷达的阵列优化方法 | |
CN110196417B (zh) | 基于发射能量集中的双基地mimo雷达角度估计方法 | |
CN105974377B (zh) | 一种对数字阵列雷达自适应调零技术的干扰方法 | |
CN114114188A (zh) | 一种低旁瓣的fda雷达通信一体化波形设计方法 | |
CN113156380B (zh) | 基于收发失配处理的主瓣噪声压制干扰自适应抑制方法 | |
CN108880647B (zh) | 一种基于频率分集阵列天线的波束控制方法 | |
CN113075649B (zh) | 一种适用于分布式网络化雷达的信号级直接定位方法 | |
CN110850383B (zh) | 一种基于共形阵的mimo雷达信号处理方法 | |
CN112130139B (zh) | 一种分布式全相参稀疏直线阵列雷达系统优化布阵方法 | |
CN116466299B (zh) | 二维子阵级稀疏阵列fpmimo雷达收发波束合成方法 | |
CN113820665B (zh) | 一种基于时间分集阵子阵划分的空域覆盖范围增强方法 | |
CN114814739A (zh) | 基于分布式fda-mimo雷达的主瓣干扰抑制方法 | |
CN111142080B (zh) | 一种针对干扰抑制的快速恒模mimo雷达跟踪波形合成方法 | |
CN110954887B (zh) | 基于球不变约束和反对称的相控阵mimo波束形成方法 | |
He et al. | Multibeam amplitude comparison problems for MIMO radar's angle measurement |
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 |