一种基于MDS的被动MIMO时差定位方法技术

本发明专利技术属于电子对抗技术领域，具体涉及一种基于MDS的被动MIMO时差定位方法。本发明专利技术通过在被动MIMO雷达定位系统下，将时差定位中的距离和作为一个纯虚数维度引入到三维位置空间，来构造一个含有纯虚维度的特殊四维空间，从而实现时差定位问题的多维标度分析，得到目标位置估计的加权最小二乘解，最后再通过解相关技术，进一步优化对目标的位置估计。本发明专利技术的方法实现了MIMO雷达系统下时差定位问题的多维标度分析，对于近场/远场辐射源目标均可实现精确定位,定位解算方法计算量小，具有良好的稳健性和准确性。

A passive MIMO time difference location method based on MDS

The invention belongs to the field of electronic countermeasures technology, and specifically relates to a passive MIMO TDOA location method based on MDS. Under the passive MIMO radar positioning system, the distance of the time difference location and a pure virtual dimension are introduced into the three-dimensional position space to construct a special four-dimensional space with pure virtual dimension, thus the multidimensional scaling analysis of the time difference localization problem is realized, and the weighted least square solution of the target location estimation is obtained. Finally, the position estimation of the target is further optimized by solving correlation techniques. The method of the invention realizes the multidimensional scaling analysis of the time difference location problem under the MIMO radar system. It can accurately locate the target of the near field / far field radiation source, and the calculation method of the location calculation method is small, and has good robustness and accuracy.

【技术实现步骤摘要】
一种基于MDS的被动MIMO时差定位方法
本专利技术属于电子对抗
，具体涉及一种基于MDS的被动MIMO时差定位方法。
技术介绍
MIMO(Multiple-InputMultiple-Output)雷达是多输入多输出雷达的简称。由于MIMO技术在移动通信技术的研究取得实际进展，在本世纪初，人们开始将MIMO技术拓展应用于雷达探测
，用于目标检测，定位跟踪等。MIMO雷达系统较传统雷达在信号检测能力、参数估计精度、目标分辨率等方面有明显优势，尤其在抗雷达信号截获、强杂波条件下的弱目标检测及慢速目标检测跟踪等方面的性能较传统雷达有显著改善。被动MIMO雷达系统是指利用已广泛存在的多基站外辐射源如:调频广播信号,电视信号,手机基站信号等作为照射源,在接收端接收经目标反射的信号,从而基于时差定位原理,获得距离和测量值，实现目标检测与定位跟踪。近年来，随着电子对抗技术的迅猛发展，传统雷达的生存能力和战场性能受到极大威胁，被动MIMO雷达系统因其隐蔽性好，定位构型灵活而受到各国军方的广泛关注。多维标度(MDS：multidimensionalscaling)是现代统计学中的重要分析方法，它是用来表示在低维空间中由任意两点的空间距离来刻画对象之间的相似度量。多维标度分析方法在现代信息技术中的无线传感器网络(WSN:wirelesssensornetwork)节点定位、蜂窝无线通信网络移动台定位及位置校正等领域展现出良好的定位稳健性。因此，将多维标度分析方法用于被动MIMO雷达系统中进行目标定位，能够有效提升雷达定位系统的抗电磁打击能力及定位稳健性能，具有很强的工程实践价值。
技术实现思路
本专利技术的目的，就是针对上述问题，提出了一种基于MDS的被动MIMO时差定位方法，该方法通过在被动MIMO雷达定位系统下，将时差定位中的距离和作为一个纯虚数维度引入到三维位置空间，来构造一个含有纯虚维度的特殊四维空间，从而实现时差定位问题的多维标度分析，得到目标位置估计的加权最小二乘解，最后再通过解相关技术，进一步优化对目标的位置估计。本专利技术所采用的技术方案为：本专利技术的方法在MIMO雷达定位系统下，假定MIMO雷达系统有M个发射站和N个接收站，其站址坐标分别表示为k＝1,2,...,M和l＝1,2,...,N，所要定位的目标辐射源的位置为u＝[x,y,z]T，定位模型如图1所示。在本系统中，真实值用等表示，测量值用τkl，dkl等表示，测量噪声用等表示。则时差测量值可表示为：将其转化为距离和表示为:则,记获得的距离和测量值向量为：d＝[d11d12...d1Nd21...dM1]T故，系统测量方程可表示为:记系统测量误差向量为：本专利技术主要包括以下步骤：a、构造一个以距离和为纯虚维度的特殊四维空间，建立时差定位问题的多维标度分析模型；b、通过多维标度分析原理，建立线性方程，得到目标向量的最小二乘估计。c、进一步，得到存在时差测量误差条件下的目标向量的加权最小二乘估计。d、通过对目标向量进行解相关分析计算，获得目标辐射源的位置估计。算法流程图如图2所示。具体的，所述步骤a中，本专利技术基于以下原理：根据上述距离和测量方程，可得：则定义待求的目标向量为：其中j为虚数单位。进行多维标度分析，构造一个含有纯虚维度的特殊四维空间，定义中心化矩阵Z为：则，建立多维标度分析中的内积矩阵B为：B＝ZZT将B分块表示为：其中具体表达式为：上述Amn,Cmn,Dmn分别为矩阵A,C,D中位置下标为(m,n)的元素,下同。进一步的，在所述步骤b中：根据多维标度分析理论，可建立如下线性方程：式中，F＝Ψ(ΨTΨ)-1，其中：定义：则上述线性方程可重写为：Gx＝g进而，得到目标向量的最小二乘解：由于在时差定位系统中，存在时差测量误差，因此上述定位解算过程中均需用测量值代替真实值，即将上述B、Ψ、F矩阵中真实值用测量值代替可得进而得到：则最终目标向量的最小二乘解为：如上所述，在时差定位系统中，存在时差测量误差，则需进一步实现存在时差测量误差条件下的多维标度分析，得到目标向量的加权最小二乘解。更进一步的，在所述步骤c中：测量误差向量ε可表示为：式中，的具体表达式如下：其中：其中：式中，见步骤b中，Φ具体表达式为：～表示不关心的矩阵单元，下同。如上所述，构成中的q(k,m)，γ(j)，h(k)，t(m,j)分别为对应矩阵和向量中位置下标为(k,j)的元素。同时，基于步骤b，上式中取：则：则：E(ε)＝0式中，Qn为测量噪声n的协方差矩阵，具体表达式为：则，最终可得目标向量的加权最小二乘解为：可得到目标位置的加权最小二乘解为：最后，在所述步骤d中：由于在上述解算过程中，我们假定目标向量:中x,y,z与是相互独立的，但式：表明它们是相关的，因此进一步的对目标向量进行解相关计算，来改进步骤c所得的定位估计由上述可知，系统测量方程可表示为:距离和测量值向量为：d＝[d11d12...d1Nd21...dM1]T则，在步骤c所得的基础上,在极大似然意义下,构建目标方程为：令式中：记fkl＝f(u,tk,rl)，则ξ的具体表达式为：其中：k＝1,...,M,l＝1,...,N则通过求解目标函数极小值，可得目标位置解相关后的加权最小二乘解：本专利技术的有益效果为，本专利技术的方法实现了MIMO雷达系统下时差定位问题的多维标度分析，对于近场/远场辐射源目标均可实现精确定位,定位解算方法计算量小，具有良好的稳健性和准确性。附图说明图1为MIMO雷达系统下的时差定位模型图；图2为定位算法流程图；图3为近场辐射源目标下定位算法的均方误差随观测误差的变化曲线；图4为远场辐射源目标下定位算法的均方误差随观测误差的变化曲线；具体实施方式下面结合附图对上述基于MDS的被动MIMO时差定位方法进行仿真验证说明，首先对系统模型作如下合理假定：1.为了便于仿真，假定MIMO雷达系统为3发3收模型，即由3个发射站和3个接收站构成MIMO雷达系统。2.假定测量误差服从均值为零的高斯分布，且误差之间相互独立。如上所述，假定3个发射站的位置坐标分别为t1＝[4300250050]T,t2＝[-43002500150]T，t3＝[0-5000100]T，单位均为m，下同；3个接收站的位置坐标分别为r1＝[05000200]T，r2＝[-4300-2500100]T，r3＝[4300-2500150]T。(1)近场辐射源目标时定位效果：如图2所示，假定近场辐射源目标位置为u＝[3008002000]T，采用上述的定位解算方法在时差测量误差变化下进行MonteCarlo仿真，并将解相关前后的定位均方误差进行了对比。从图2中可以看出，本文提出的方法对于近场辐射源目标可以进行有效的定位解算，随着测量误差的增大，定位的均方误差平稳增加，表现出良好的稳健性和准确性。同时由图中对比可以看出，基于MDS的时差定位估计经过解相关之后定位效果又进一步得到提升。(2)远场辐射源目标时定位效果：如图3所示，假定近场辐射源目标位置为u＝[400060002000]T，采用上述的定位解算方法在时差测量误差变化下进行MonteCarlo仿真，并将解相关前后的定位均方误差进行了对比。从图3中可以看出，本文提出的方法对于远场辐射源目标可本文档来自技高网...

【技术保护点】
1.一种基于MDS的时差定位方法，该方法用于被动MIMO雷达系统，所述被动MIMO雷达系统有M个发射站和N个接收站，其站址坐标分别表示为：发射站为

【技术特征摘要】
1.一种基于MDS的时差定位方法，该方法用于被动MIMO雷达系统，所述被动MIMO雷达系统有M个发射站和N个接收站，其站址坐标分别表示为：发射站为k＝1,2,...,M和接收站为l＝1,2,...,N，并设定所要定位的目标辐射源的位置为u＝[x,y,z]T；将时差测量值表示为：其中，tkl为测量值，为真实值，为测量噪声；将其转化为距离和表示为：其中，dkl为测量值，为真实值，为测量噪声；记获得的距离和测量值向量为：d＝[d11d12...d1Nd21...dM1]T则系统测量方程为：系统测量误差向量为：其特征在于，所述时差定位方法包括以下步骤：a、构造一个以距离和为纯虚维度的四维空间，建立时差定位问题的多维标度分析模型；b、通过多维标度分析原理，建立线性方程，得到目标向量的最小二乘估计；c、根据步骤b的结果，获取存在时差测量误差条件下的目标向量的加权最小二乘估计；d、通过对目标向量进行解相关分析计算，获得目标辐射源的位置估计。2.根据权利要求1所述的一种基于MDS的时差定位方法，其特征在于，所述步骤a的具体方法为：根据距离和测量方程，可得：则定义待求的目标向量为：其中j为虚数单位；进行多维标度分析，构造一个含有纯虚维度的四维空间，定义中心化矩阵Z为：则，建立多维标度分析中的内积矩阵B为：B＝ZZT将B分块表示为：其中：上述Amn,Cmn,Dmn分别为矩阵A,C,D中位置下标为(m...

【专利技术属性】
技术研发人员：李万春，陈锐滨，扶彩霞，郭昱宁，
申请(专利权)人：电子科技大学，
类型：发明
国别省市：四川,51