一种基于先导信息的MIMO雷达波形正交方法

一种基于先导信息的mimo雷达波形正交方法
技术领域
1.本发明涉及雷达信号处理技术领域，尤其是一种基于先导信息的mimo雷达波形正交方法。


背景技术：

2.随着各种隐身技术以及干扰技术的发展，雷达对目标探测和跟踪面临着前所未有的挑战，寻求新的雷达体制和信号处理手段，来提升雷达探测性能和自身抗干扰能力，是雷达工程人员追寻的目标。mimo雷达可以实现多波束，且波束扫描范围较宽，对于存在较大起伏目标的探测效果显著。与具有m个相同阵元的相控阵雷达相比，mimo雷达峰值功率和发射天线增益都仅为相控阵雷达的1/m，因此mimo雷达具有抗截获优势。mimo雷达可利用少数实体阵元形成多个虚拟阵列，进而扩大了阵列孔径，提高了阵列角度分辨率。mimo雷达每个子阵发射的信号分量正交，其灵活性较强，自由度高。因此，mimo雷达能满足强杂波环境下的若目标检测和多目标检测与跟踪。正交波形影响着mimo雷达对弱目标的检测、影响其匹配滤波后信号的恢复效果以及雷达系统的距离分辨率，是有效提升mimo雷达目标参数检测和跟踪性能的急需解决的难题。


技术实现要素：

3.本发明所要解决的技术问题在于，提供一种基于先导信息的mimo雷达波形正交方法，能够实现mimo雷达波形正交，进而提高雷达系统多普勒分辨率、抗截获能力以及对弱小目标的探测与跟踪。
4.为解决上述技术问题，本发明提供一种基于先导信息的mimo雷达波形正交方法，包括如下步骤：
5.步骤1、将窄带信号先导信息先调制在mimo雷达各子阵发射波形上，通过雷达数字发射机发射波形，其他子阵通过数字接收机接收波形并通过中频滤波器将其他子阵的先导信息提取出来；
6.步骤2、将步骤1中提取出的mimo雷达其他各个子阵的单载频信号先导信息经过解调后，对其进行短时傅里叶变换stft，得到mimo雷达各个子阵先导信息对应的时频图；
7.步骤3、通过对步骤2的先导信息的时频图特征进行提取，若发现不同mimo雷达子阵之间存在相同的工作参数，即不同子阵之间波形不正交，通过代价函数对波形的参数进行优化，使得各个子阵的工作参数独立，保持mimo雷达波形正交；
8.步骤4、步骤3涉及的波形优化根据mimo雷达波形设计正交原则构建波形优化代价函数，用模拟退火算法对波形编码参数进行优化，构造mimo雷达发射波形；
9.步骤5、构建mimo雷达发射波形数学模型，提供步骤4所需的波形设计所需优化的频率编码参数脉位编码参数和相位编码参数
10.步骤6、提供步骤4中所需要的代价函数，对步骤5中的波形参数进行优化；
11.步骤7、将步骤6优化后的波形用于mimo雷达各个子阵进行目标探测，探测波束的回波信号进行数字波束合成、脉冲压缩以及fft等信号处理方式后经过恒虚警检测生成相应的目标点迹，在探测过程中，mimo雷达各子阵间周期性广播先导信息，再次执行步骤1至步骤6，进行mimo雷达波形正交设计。
12.优选的，步骤2中，短时傅里叶变换公式如下所示：
[0013][0014]
其中k＝0,1,...,k-1,k∈n，k表示多普勒频率索引，u＝0,1,...,u-1,u∈n,u表示时间索引，w(
·
)是一个时域的窗，m＝0,1,...,m-1表示慢时间的索引。
[0015]
优选的，步骤5中，构建mimo雷达发射波形数学模型，提供步骤4所需的波形设计所需优化的编码参数具体为：设mimo雷达有q个子阵，一个cpi内的脉冲个数为n，每个脉冲内采用码长为m的相位编码，第q个子阵发射信号的数学模型为：
[0016][0017]
其中，t
p
为脉冲宽度，tr为子脉冲的平均脉冲周期，δtr为最小脉位跳变间隔，ts为相位编码码元宽度，ξn为n个子脉冲的脉位编码，为第n个子脉冲的载波，δf为频率最小间隔，μ为调频斜率，为第n个子脉冲的第m个码元，从发射信号的数学模型可知，波形优化参数为用向量表示记为ω＝[ξ1,ξ2,
…
,ξn]、
[0018]
优选的，步骤6中，采用峰值旁瓣电平准则(psl)进行正交波形设计,以发射波形的相位编码矩阵、载波的频率编码序列以及脉位捷变序列作为优化变量，构造的代价函数为：
[0019][0020][0021][0022]
ω∈perms{[ξ1,ξ2,
…
,ξ
p
]}
[0023]
i＝1,2,
…
,m,k＝1,2,
…
,n,j＝1,2,
…
,q
[0024][0025]
其中，s
p
(t)和sq(t)为不同子阵的发射信号，当p＝q时，c
pq
(τ)为发射信号的自相关函数；当p≠q时，c
pq
(τ)为发射信号的互相关函数。
[0026]
本发明的有益效果为：本发明能够解决mimo雷达波形设计过程中，如何通过各mimo雷达子阵间使用引导信息，广播各自的波形工作参数，使得各个子阵避开相同的波形
工作参数，提高子阵之间的正交性，以免削弱mimo雷达系统的主副瓣比、正交性以及主瓣宽度，最终实现mimo雷达波形正交，进而提高雷达系统多普勒分辨率、抗截获能力以及对弱小目标的探测与跟踪。
附图说明
[0027]
图1为本发明的方法流程示意图。
[0028]
图2为本发明优化后的mimo雷达发射波形时域图。
[0029]
图3为本发明优化后的mimo雷达发射波形频域图。
具体实施方式
[0030]
如图1所示，一种基于先导信息的mimo雷达波形正交方法，包括如下步骤：
[0031]
步骤1、将窄带信号先导信息先调制在mimo雷达各子阵发射波形上，通过雷达数字发射机发射波形，其他子阵通过数字接收机接收波形并通过中频滤波器将其他子阵的先导信息提取出来。
[0032]
步骤2、将步骤1中提取出的mimo雷达其他各个子阵的单载频信号先导信息经过解调后，对其进行短时傅里叶变换(stft)，得到mimo雷达各个子阵先导信息对应的时频图。短时傅里叶变换公式如下所示：
[0033][0034]
其中k＝0,1,...,k-1,k∈n，k表示多普勒频率索引，u＝0,1,...,u-1,u∈n,u表示时间索引。w(
·
)是一个时域的窗。m＝0,1,...,m-1表示慢时间的索引。
[0035]
步骤3、通过对步骤2的先导信息的时频图特征进行提取，若发现不同mimo雷达子阵之间存在相同的工作参数，即不同子阵之间波形不正交，通过代价函数对波形的参数进行优化，使得各个子阵的工作参数独立，保持mimo雷达波形正交。
[0036]
步骤4、步骤3涉及的波形优化根据mimo雷达波形设计正交原则构建波形优化代价函数，用模拟退火算法对波形编码参数进行优化，构造mimo雷达发射波形,优化后的波形如图2所示。
[0037]
步骤5、构建mimo雷达发射波形数学模型，提供步骤4所需的波形设计所需优化的频率编码参数脉位编码参数和相位编码参数
[0038]
设mimo雷达有q个子阵，一个cpi内的脉冲个数为n，每个脉冲内采用码长为m的相位编码，第q个子阵发射信号的数学模型为：
[0039][0040]
其中，t
p
为脉冲宽度，tr为子脉冲的平均脉冲周期，δtr为最小脉位跳变间隔，ts为相位编码码元宽度，ξn为n个子脉冲的脉位编码，为第n个子脉冲的载波，δf为频率最小间隔，μ为调频斜率，为第n个子脉冲的第m个码元，从发射信号的数学模
型可知，波形优化参数为用向量表示记为ω＝[ξ1,ξ2,
…
,ξn]、
[0041]
步骤6、提供步骤4中所需要的代价函数，对步骤5中的波形参数进行优化，此处采用峰值旁瓣电平准则(psl)进行正交波形设计,以发射波形的相位编码矩阵、载波的频率编码序列以及脉位捷变序列作为优化变量，构造的代价函数为：
[0042][0043][0044][0045]
ω∈perms{[ξ1,ξ2,
…
,ξ
p
]}
[0046]
i＝1,2,
…
,m,k＝1,2,
…
,n,j＝1,2,
…
,q
[0047][0048]
其中，s
p
(t)和sq(t)为不同子阵的发射信号，当p＝q时，c
pq
(τ)为发射信号的自相关函数；当p≠q时，c
pq
(τ)为发射信号的互相关函数。
[0049]
步骤7、将步骤6优化后的波形用于mimo雷达各个子阵进行目标探测，探测波束的回波信号进行数字波束合成、脉冲压缩以及fft等信号处理方式后经过恒虚警检测生成相应的目标点迹，在探测过程中，mimo雷达各子阵间周期性广播先导信息，再次执行步骤1至步骤6，进行mimo雷达波形正交设计。
[0050]
本发明是为了解决mimo雷达波形设计过程中，如何通过各mimo雷达子阵间使用引导信息，广播各自的波形工作参数，使得各个子阵避开相同的波形工作参数，提高子阵之间的正交性，以免削弱mimo雷达系统的主副瓣比、正交性以及主瓣宽度，最终实现mimo雷达波形正交，进而提高雷达系统多普勒分辨率、抗截获能力以及对弱小目标的探测与跟踪。
[0051]
当mimo雷达子阵使用相同的波形工作参数(如载频、带宽、调频斜率和相位)时，mimo雷达波形不能保持正交，进而影响雷达对目标探测和跟踪的精度。因此mimo雷达子阵之间可以先各自转发含有自己波形工作参数的先导信息，先导信息调制到雷达波形中，通过雷达子阵天线发射出去，各子阵接收机经过中频滤波和解调后得到含有其他子阵先导信息的信号，将该信号通过短时傅里叶变换后得到各自子阵先导信息的时频图，由时频图判断各子阵之间波形是否正交。若波形之间不正交，使用代价函数优化波形中的参数，再将优化后的波形用于mimo雷达，在mimo雷达进行目标探测过程中，周期性地发射先导信息来不断感知mimo雷达各子阵之间是否正交。
[0052]
本发明解决了目标探测过程中mimo雷达各子阵因为工作参数相同而引起波形正交性差的问题，有效提升了mimo雷达目标参数检测和跟踪性能。利用mimo雷达各子阵先发射一段包含各自工作参数(如载频、调频斜率、带宽和相位等)的先导信息，其他mimo雷达子
阵接收到这段先导信息后会避免使用相同的工作参数，若mimo雷达其他子阵也在使用相同的工作参数时，通过波形优化对mimo雷达各子阵的发射波形进行频率、相位和脉位联合编码，编码后的mimo雷达各子阵的波形正交。优化后的波形用于目标探测，在探测过程中，mimo雷达各子阵再周期性发射和接收先导信息来感知相互间波形正交性。使用先导信息极大的提高了mimo雷达波形的正交性，进而提高雷达系统多普勒分辨率、抗截获能力以及对弱小目标的探测与跟踪。

技术特征：
1.一种基于先导信息的mimo雷达波形正交方法，其特征在于，包括如下步骤：步骤1、将窄带信号先导信息先调制在mimo雷达各子阵发射波形上，通过雷达数字发射机发射波形，其他子阵通过数字接收机接收波形并通过中频滤波器将其他子阵的先导信息提取出来；步骤2、将步骤1中提取出的mimo雷达其他各个子阵的单载频信号先导信息经过解调后，对其进行短时傅里叶变换stft，得到mimo雷达各个子阵先导信息对应的时频图；步骤3、通过对步骤2的先导信息的时频图特征进行提取，若发现不同mimo雷达子阵之间存在相同的工作参数，即不同子阵之间波形不正交，通过代价函数对波形的参数进行优化，使得各个子阵的工作参数独立，保持mimo雷达波形正交；步骤4、步骤3涉及的波形优化根据mimo雷达波形设计正交原则构建波形优化代价函数，用模拟退火算法对波形编码参数进行优化，构造mimo雷达发射波形；步骤5、构建mimo雷达发射波形数学模型，提供步骤4所需的波形设计所需优化的频率编码参数脉位编码参数和相位编码参数步骤6、提供步骤4中所需要的代价函数，对步骤5中的波形参数进行优化；步骤7、将步骤6优化后的波形用于mimo雷达各个子阵进行目标探测，探测波束的回波信号进行数字波束合成、脉冲压缩以及fft等信号处理方式后经过恒虚警检测生成相应的目标点迹，在探测过程中，mimo雷达各子阵间周期性广播先导信息，再次执行步骤1至步骤6，进行mimo雷达波形正交设计。2.如权利要求1所述的基于先导信息的mimo雷达波形正交方法，其特征在于，步骤2中，短时傅里叶变换公式如下所示：其中k＝0,1,...,k-1,k∈n，k表示多普勒频率索引，u＝0,1,...,u-1,u∈n,u表示时间索引，w(
·
)是一个时域的窗，m＝0,1,...,m-1表示慢时间的索引。3.如权利要求1所述的基于先导信息的mimo雷达波形正交方法，其特征在于，步骤5中，构建mimo雷达发射波形数学模型，提供步骤4所需的波形设计所需优化的编码参数具体为：设mimo雷达有q个子阵，一个cpi内的脉冲个数为n，每个脉冲内采用码长为m的相位编码，第q个子阵发射信号的数学模型为：其中，t
p
为脉冲宽度，t
r
为子脉冲的平均脉冲周期，δt
r
为最小脉位跳变间隔，t
s
为相位编码码元宽度，ξ
n
为n个子脉冲的脉位编码，为第n个子脉冲的载波，δf为频率最小间隔，μ为调频斜率，为第n个子脉冲的第m个码元，从发射信号的数学模型可知，波形优化参数为用向量表示记为ω＝[ξ1,ξ2,
…
,
ξ
n
]、4.如权利要求1所述的基于先导信息的mimo雷达波形正交方法，其特征在于，步骤6中，采用峰值旁瓣电平准则(psl)进行正交波形设计,以发射波形的相位编码矩阵、载波的频率编码序列以及脉位捷变序列作为优化变量，构造的代价函数为：编码序列以及脉位捷变序列作为优化变量，构造的代价函数为：编码序列以及脉位捷变序列作为优化变量，构造的代价函数为：ω∈perms{[ξ1,ξ2,
…
,ξ
p
]}i＝1,2,
…
,m,k＝1,2,
…
,n,j＝1,2,
…
,q其中，s
p
(t)和s
q
(t)为不同子阵的发射信号，当p＝q时，c
pq
(τ)为发射信号的自相关函数；当p≠q时，c
pq
(τ)为发射信号的互相关函数。

技术总结
本发明公开了一种基于先导信息的MIMO雷达波形正交方法，MIMO雷达子阵之间可以先各自转发含有自己波形工作参数的先导信息，先导信息调制到雷达波形中，通过雷达子阵天线发射出去，各子阵接收机经过中频滤波和解调后得到含有其他子阵先导信息的信号，将该信号通过短时傅里叶变换后得到各自子阵先导信息的时频图，由时频图判断各子阵之间波形是否正交。若波形之间不正交，使用代价函数优化波形中的参数，再将优化后的波形用于MIMO雷达，在MIMO雷达进行目标探测过程中，周期性地发射先导信息来不断感知MIMO雷达各子阵之间是否正交。本发明能够实现MIMO雷达波形正交，进而提高雷达系统多普勒分辨率、抗截获能力以及对弱小目标的探测与跟踪。与跟踪。与跟踪。


技术研发人员：董浩 胡文 姚翼荣 王原正 张逸雯
受保护的技术使用者：南京航空航天大学
技术研发日：2021.11.15
技术公布日：2022/3/8