一种基于模糊函数的频控阵MIMO雷达目标定位方法与流程

本发明属于频控阵雷达技术领域，尤其涉及一种基于模糊函数的频控阵mimo雷达目标定位方法。

背景技术：

频控阵雷达因其特有的距离-角度依赖性而在雷达系统、无线通信、雷达成像、目标估计与追踪、干扰及抗干扰等领域有着非常广泛的应用。相对于相控阵而言，均匀线形频控阵能够形成具有角度-距离依赖性的发射波束，对远场目标能够进行角度维和距离维的估计。传统的目标估计方法music算法在较低信噪比情况下，对运动目标角度和距离的估计误差较大，且在多目标角度和距离接近的情况下，无法进行有效估计。而目前针对此问题的解决主要通过配置线性频控阵的参数如阵元数目、阵元间隔等，能够提高对目标信息的估计，其实现较为复杂，为此，本次发明提出一种基于模糊函数和music的mimo频控阵目标定位方法，即利用模糊矩阵替代相关矩阵进行特征值分解求解噪声子空间和信号子空间，然后在时频谱进行二维搜索求得目标的角度和距离维信息。

技术实现要素：

针对现有技术中的上述不足，本发明提供的一种基于模糊函数的频控阵mimo雷达目标定位方法解决了频控阵mimo在较低信噪比情况下下无法对目标信息准确估计的问题。

为了达到以上目的，本发明采用的技术方案为：

本方案提供一种基于模糊函数的频控阵mimo雷达目标定位方法，包括以下步骤：

s1、根据天线的频控发射阵列、相控阵接收阵列、频控阵列发射阵元的等间距、相控阵接收阵元的等间距以及远场处运动目标的坐标(θk,rk)构建频控阵mimo，其中，θ表示运动目标的方位角，取值范围为[-90°,90°]，r表示观测点到参考点的距离，取值范围为[rmin,∞)，rmin表示使观测点满足远场假设的最小距离，k表示远场处的第k个运动目标；

s2、根据所述频控阵mimo构建频控阵mimo接收阵列接收的信号矩阵；

s3、根据所述信号矩阵分别利用模糊函数和多信号分类算法music，求解雷达目标的角度维和距离维信息，从而完成对雷达目标的定位。

进一步地，所述步骤s1中频控发射阵列按频控发射阵列的等间距dt分布m个发射阵元，且各发射阵元的频偏值δfx依次按线性递增。

再进一步地，所述步骤s1中相控阵接收阵列按相控阵接收阵列的等间距dr分布n个相控阵接收阵元。

再进一步地，所述步骤s2中mimo频控阵接收阵列接收的信号矩阵x的表达式如下：

x＝[x(t1),x(t2),...x(tn)]

其中，tn表示采样快拍时刻，x(tn)表示目标在采样快拍时接收的信号值。

再进一步地，所述步骤s3包括如下步骤：

s301、分别在所述信号矩阵的第i个虚拟接收阵元和第l个虚拟接收阵元接收到的信号中取互模糊函数，并根据所述互模糊函数取期望值；

s302、根据所述互模糊函数和期望值组成模糊矩阵raf，并对模糊矩阵raf进行特征值分解，得到噪声子空间和信号子空间；

s303、利用多信号分类算法music对所述噪声子空间和信号子空间进行二维谱峰搜索，得到雷达目标的角度和距离，从而完成对雷达目标的定位。

再进一步地，所述步骤s302中对模糊矩阵raf进行特征值分解的表达式如下：

其中，表示由模糊函数矩阵raf的最大的k个特征值对应的特征向量所张成的信号子空间，表示由模糊函数矩阵raf的最小特征值对应的特征向量张成的噪声子空间，∧表示模糊函数矩阵raf的特征值矩阵,m表示频控阵发射阵列阵元个数，n表示相控阵接收阵列阵元个数，n(t)表示加性高斯白噪声，k表示远场处运动目标的总数。

再进一步地，所述步骤s303中雷达目标的角度和距离paf(θ,r)的表达式如下：

其中，θ表示运动目标的方位角，取值范围为[-90°,90°]，r表示观测点到参考点的距离，取值范围为[rmin,∞)，rmin表示使观测点满足远场假设的最小距离，afda表示频控阵mimo雷达的导向矢量，表示频控阵mimo雷达导向矢量的共轭转置，表示由模糊函数矩阵的最小特征值对应的特征向量张成的噪声子空间，h表示向量的共轭转置。

本发明的有益效果：

本发明根据天线的频控发射阵列、相控阵接收阵列、频控阵列发射阵元的等间距、相控阵接收阵元的等间距以及远场处运动目标的坐标构建频控阵mimo，并根据所述频控阵mimo构建频控阵mimo接收阵列接收的信号矩阵，再利用模糊函数和多信号分类算法music，求解雷达目标的角度维和距离维信息，从而完成对雷达目标的定位，即利用模糊矩阵替代相关矩阵进行特征值分解求解噪声子空间和信号子空间，然后在时频谱进行二维搜索求得目标的角度和距离维信息，采用该目标定位方法可对远场目标定位性能有较大改善，解决了mimo频控阵在较低信噪比情况下下无法对目标信息准确估计的问题。

附图说明

图1为本发明的方法流程图。

图2为本发明中频控阵mimo的模型示意图；

图3为利用传统music算法对两个运动目标得到的二维谱图；

图4为采用本发明方法后对两个运动目标得到的二维谱图；

图5为传统music算法和采用本发明方法后对两个运动目标的信噪比snr曲线示意图。

具体实施方式

下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。

实施例

本发明利用模糊矩阵替代相关矩阵进行特征值分解求解噪声子空间和信号子空间，然后在时频谱进行二维搜索求得目标的角度和距离维信息，采用该目标定位方法可对远场目标定位性能有较大改善，解决了mimo频控阵在较低信噪比情况下下无法对目标信息准确估计的问题。

如图1所示，本发明公开了一种基于模糊函数的频控阵mimo雷达目标定位方法，其特征在于，包括以下步骤：

s1、根据天线的频控发射阵列、相控阵接收阵列、频控阵列发射阵元的等间距、相控阵接收阵元的等间距以及远场处运动目标的坐标(θk,rk)构建频控阵mimo，其中，θ表示运动目标的方位角，取值范围为[-90°,90°]，r表示观测点到参考点的距离，取值范围为[rmin,∞)，rmin表示使观测点满足远场假设的最小距离，k表示远场处的第k个运动目标，所述频控发射阵列按频控发射阵列的等间距dt分布m个发射阵元，且各发射阵元的频偏值δfx依次按线性递增，所述相控阵接收阵列按相控阵接收阵列的等间距dr分布n个相控阵接收阵元，在具体实施例中，如图2-图3所示，设置天线的阵元个数，工作频率，相邻天线阵元间的频偏值以及等间距，假设频控阵中的天线均采用全向天线，在不考虑能量随距离衰减情况下，设置频控阵发射阵列阵元个数m＝4，相控阵接收阵列阵元个数n＝4，第一个全向天线(参考阵元)的载波频率为f0＝1ghz，发射端阵元之间的频偏值为δf＝3khz，阵元间距为λ为波长；

s2、根据所述频控阵mimo构建频控阵mimo接收阵列接收的信号矩阵，在具体实施例中，假定距离非常近的频控阵发射阵列和相控阵接收阵列分别放置m和n个阵元，则在远场目标pk处可以认为两个阵列是在同一空间区域位置，对于位于(θk,rk)处的远场目标，两个相邻发射阵元发射的信号到达该目标时期相位差为：

其中，d表示频控阵发射阵元的间距，δf频控阵发射阵元之间的线性频偏，c表示光速。

因此，发射阵列的导向矢量at(θk,rk)为：

其中，以发射阵列中的第一个阵元为参考阵元，1表示参考阵元的导向矢量，表示第2个发射阵元的导向矢量，第m个发射阵元的导向矢量，t表示向量的转置。

发射阵列发射的信号在目标处反射并由相控阵接收阵列接收，接收阵列的导向矢量br(θk)为：

其中，以接收阵列中的第一个阵元为参考阵元，1表示其导向矢量，表示第2个接收阵元导向矢量，表示第n个接收阵元的导向矢量，t表示向量的转置。

此时，两个相邻发射阵元之间的相位差为：

其中，d表示频控阵发射阵元间距，δf表示频控阵发射阵元之间的线性频偏，c表示光速。

发射导向矢量atr(θk,rk)在接收阵列处可以写为：

其中，以发射阵列中的第一个阵元为参考阵元，1表示参考阵元的导向矢量，表示经过目标反射后到达接收阵列时第2个发射阵元的导向矢量，表示经过目标反射后到达接收阵列时第m个发射阵元的导向矢量，t表示向量的转置。

假设目标为运动目标，其散射信号sk(t)可以表示为：

其中，βk表示第k个运动目标的雷达横截面积的复振幅，fk和uk分别表示远场处第k个运动目标的多普勒频移和调频斜率，j表示复向量，t表示时间，即假设目标信号的反射信号属于非平稳信号,则接收阵列接收到的数据矩阵y可以写成：

其中，b表示发射的基带正交波形，n(t)表示加性高斯白噪声，k表示远场处的第k个运动目标，k表示远场处运动目标的总数，br表示相控阵接收阵列的导向矢量，θk表示远场处第k个运动目标的角度，rk表示远场处第k个运动目标的距离，sk表示远场运动目标散射信号，表示频控阵发射阵列在接收阵列处的导向矢量的转置，在接收天线阵列处经过匹配滤波后可以得到输出信号x(t)：

其矩阵形式为：

x(t)＝(br·at)s(t)+n(t)(9)

其中，br＝[br(θ1,r1),br(θ2,r2),...br(θk,rk)]表示远场处第k个运动目标相对于接收阵列的导向矢量，at＝[atr(θ1,r1),atr(θ2,r2),...atr(θk,rk)]表示远场处第k个运动目标相对于发射阵列的导向矢量，取j个快拍组成接收信号矩阵：

x＝[x(t1),x(t2),...x(tn)](10)；

其中，x(tn)表示目标在tn采样快拍时接收信号值，tn表示采样快拍时刻。

s3、根据所述信号矩阵分别利用模糊函数和多信号分类算法music，求解雷达目标的角度维和距离维信息，从而完成对雷达目标的定位，其实施方法如下：

s301、分别在所述信号矩阵的第i个虚拟接收阵元和第l个虚拟接收阵元接收到的信号中取互模糊函数和取期望值；

s302、根据所述互模糊函数和期望值组成模糊矩阵raf，并对模糊矩阵raf进行特征值分解，得到噪声子空间和信号子空间；

s303、利用多信号分类算法music对所述噪声子空间和信号子空间进行二维谱峰搜索，得到雷达目标的角度和距离，从而完成对雷达目标的定位。

在具体实施例中，对于线性均匀排列的频控阵发射，相控阵接收阵列而言，其接收信号x(t)为：

x(t)＝y(t)+n(t)＝as(t)+n(t)(11)

其中，y(t)为波达信号，br＝[br(θ1,r1),br(θ2,r2),...br(θk,rk)]表示远场处的中k个运动目标相对于接收阵列的导向矢量，且k＝1,2,3..k，k表示远场运动目标的总数at＝[atr(θ1,r1),atr(θ2,r2),...atr(θk,rk)]表示远场处的k个运动目标相对于发射阵列的导向矢量，n(t)表示加性高斯白噪声。

对于频控阵mimo雷达而言，分别对信号矩阵中的第i个虚拟接收阵元和第l个虚拟接收阵元接收到的信号取互模糊函数afi,l(td,fd)表示为：

其中，td,fd分别表示时间延迟和频率偏移，i,l＝1,...mn，m表示频控阵发射列阵元个数，n表示相控阵接收阵列阵元个数，表示在时间域的积分，xi(t)表示第i个虚拟接收阵元接收到的信号，表示第l个虚拟接收阵元接收到的信号的共轭，t表示时间，将公式(11)代入公式(12)可以得到：

其中，n(t)表示加性高斯白噪声，yl表示第l个接收阵元的波达信号，yi表示第i个接收阵元的波达信号，a＝br·at表示频控阵mimo雷达的导向矢量，n^*(t)表示加性高斯白噪声的共轭。

对公式(13)求期望，可以得到：

其中，e(afi,l(td,fd))表示对模糊矩阵afi,l(td,fd)求期望，表示第l个接收阵元的波达信号的共轭，e(·)表示求期望。

由于目标散射信号和高斯白噪声不相关，所以有

公式(15)可以为：

对于常见的信号如调频信号等，其模糊函数在时频域是一条过原点的直线，所以对其在时频域采取加窗处理，降低其他信号干扰。然后以公式(16)得到的模糊函数期望组成模糊矩阵raf，对其进行特征值分解，得到噪声子空间和信号子空间：

其中，表示由模糊函数矩阵raf的最大的k个特征值对应的特征向量所张成的信号子空间，表示由模糊函数矩阵raf的最小特征值对应的特征向量张成的噪声子空间，∧表示模糊函数矩阵raf的特征值矩阵,m表示频控阵发射阵列个数，n表示相控阵接收阵列个数，n(t)表示加性高斯白噪声，k表示远场运动目标的总数，和多信号分类算法music的处理方法类似，对空间谱进行二维谱峰搜索，即为目标的角度和距离：

其中，θ表示方位角，取值范围为[-90°,90°]，r表示观测点到参考点的距离，取值范围为[rmin,∞)，rmin表示使观测点满足远场假设的最小距离，afda＝br·at表示频控阵mimo雷达的导向矢量，表示由模糊函数矩阵raf的最小特征值对应的特征向量张成的噪声子空间，h表示向量的共轭转置。

本实施例中，对于多个运动目标而言，在目标之间角度和距离较近的情况下，传统的算法已经不能对目标的角度和距离信息进行有效的估计，假设两个运动目标分别位于(10km,30°)和(10km，35°)时，目标1的散射信号的起始频率为0.4，截止频率为0，目标2的散射信号的起始频率为0.4，截止频率为0。在信噪比为0db的情况下，分别使用传统算法和本文所提算法对目标进行定位，分别画出其二维时频谱图如图3和如图4。如图5所示，最后在同样的仿真条件下通过1000次蒙特卡洛实验得到本文所提算法在角度和距离维关于snr的仿真实验结果。从图3可以看出，信噪比为0db时，对于两个角度相近，距离相同的运动目标而言，传统的music算法的二维空间谱图只存在一个尖峰，说明传统算法已经不能对目标进行有效估计。从图4可以看出，本文算法的二维空间谱图存在两个尖峰，证明本发明中的算法可以对目标进行有效估计。从图5可看出，本发明随着信噪比snr的变换曲线，即随着信噪比的增加，从而提高目标精度。

本发明通过以上设计解决了传统music算法在较低信噪比不能够对两个相近的目标进行准确定位的缺点，实现了可以准确得到目标位置信息，从而对目标进行定位。