用于分布式多点定位监视系统的半正定松弛时差定位方法技术方案

本发明专利技术属于分布式多点定位监视系统的实时定位技术领域，特别涉及一种用于分布式多点定位监视系统的半正定松弛时差定位方法。本发明专利技术首先构建时差定位方程，再对信号源位置即待定位的目标位置进行最大似然估计，引入辅助矢量，将距离差定位方程转化为约束最小二乘问题，再对所述辅助矢量进行加权最小二乘求解，再利用加权最小二乘解算的辅助矢量初始估计松弛等式约束，构造新的代价函数，利用凸半正定规划优化求解辅助矢量和辅助矢量转置的变量的值，并通过特征值分解得到辅助矢量的值，最后根据求得的辅助矢量的值与信号源位置之间的关系获取信号源的位置信息。本发明专利技术避免了传统迭代算法中出现的局部收敛和求解发散的问题，提高了定位的精度。

【技术实现步骤摘要】
用于分布式多点定位监视系统的半正定松弛时差定位方法
本专利技术属于分布式多点定位监视系统的实时定位
，特别涉及一种用于分布式多点定位监视系统的半正定松弛时差定位方法。
技术介绍
随着我国经济的飞速发展，航班量也在飞速增长，对空管系统管制员的指挥和设施设备的维护提出了更高的要求，机场场面监视雷达是一次雷达，主要用于监视机场场面的飞机及车辆，利用目标对电磁波的自身辐射或反射特性发现目标。设备维护人员在维护机场场面监视雷达的过程中包括如下三方面问题：飞机目标分裂、存在一定的盲区、容易出现假目标。这些问题的存在扰乱了管制员正常的指挥工作，因此，分布式多点定位系统的引入必要而且必须，这项技术是机场场面监视雷达的升级和补充，并且能够通过增加接收站或者改变接收站的布局实现系统的扩展。半正定松弛时差定位方法是时差定位方法非常重要的一部分，利用半正定松弛时差定位方法来确定目标的位置信息，并将定位结果通过终端显示监控系统呈现给用户，从而提高机场管制员对机场场面飞机和引导车监控管制能力。传统的半正定松弛时差定位方法通常使用泰勒级数算法、Chan算法、MDS算法来确定目标的位置信息，使用泰勒级数算法对迭代初始值的选取有要求，会出现局部收敛或者发散的现象；Chan算法和MDS算法在测量误差较大时导致对信号源的位置判断不准确，从而导致时差定位方法经常出现高度非线性的问题。
技术实现思路
本专利技术为了克服上述现有技术的不足，提供了用于分布式多点定位监视系统的半正定松弛时差定位方法，本方法对迭代初始值的选取没有要求，能够收敛到全局极小点，不会出现局部收敛或者发散的现象；而且定位准确，在测量误差较大时仍能对信号源的位置进行较好的估计。为实现上述目的，本专利技术采用了以下技术措施：用于分布式多点定位监视系统的半正定松弛时差定位方法，包括以下步骤：S1、构建时差定位方程；S2、对信号源位置即待定位的目标位置进行最大似然估计；S3、引入辅助矢量，将距离差定位方程转化为约束最小二乘问题；S4、对所述辅助矢量进行加权最小二乘求解；S5、利用加权最小二乘解算的辅助矢量初始估计松弛等式约束，构造新的代价函数；S6、利用凸半正定规划优化求解辅助矢量和辅助矢量转置的变量的值，并通过特征值分解得到辅助矢量的值；S7、根据求得的辅助矢量的值与信号源位置之间的关系获取信号源位置信息。优选的，步骤S1的具体步骤包括：利用接收站主站位置信息测量参数数据建立相应的观测模型，确定的接收站主站的站址信息，利用目标到接收站主站的距离、目标到接收站辅站的距离差构建距离差定位方程：其中，为到达接收站主站和接收站辅站时间差的测量值，c为光速，ri,k-r1,k为信号源发射的应答信号从第k个目标位置uk到接收站si和接收站主站s1之间的真实时间差，ni1,k为测量误差，i为接收站辅站的索引，k为目标索引。优选的，步骤S2的具体步骤包括：利用目标到接收站主站和第i个接收站辅站的距离差测量误差的协方差矩阵表示目标位置的最大似然估计：其中，J(u)为信号源位置u的高度非线性、非凸代价函数，为所构成的距离差估计组成的向量，f(u)为信号源发射的应答信号从第k个目标位置uk到接收站si和接收站主站s1之间的真实时间差组成的向量，Q为测量误差ni1,k服从均值为零的协方差矩阵，T为矩阵的转置，为使得代价函数J(u)最小时所需信号源位置u的值。优选的，步骤S3的具体步骤包括：引入辅助矢量即为目标与接收站主站的位置和距离，转化为约束最小二乘问题，并在此基础上转换为凸的半正定规划问题，将原代价函数整理得到的矩阵如下：其中，ηk＝Bknk，Bk＝diag{r2,k,...,rN,k}，k＝1,2，R＝[s2-s1,...,sN-s1]T，O(N-1)×3为(N-1)×3的0矩阵，0N-1为N-1的0行向量，为目标1的到达接收站主站和接收站辅站距离差估计组成的向量，为目标2的到达接收站主站和接收站辅站距离差估计组成的向量。进一步的，步骤S4的具体步骤包括：对原代价函数整理得到的矩阵进行加权最小二乘求解，得到的辅助矢量估计值如下：其中，加权矩阵W＝E[ηηT]-1＝(BQBT)-1，B＝diag{B1,B2}表示以B1和B2为对角块的矩阵，Q＝E[nnT]为噪声协方差矩阵，n为到达接收站主站和接收站辅站距离差测量误差的矩阵形式，Gy-h＝η，y为辅助矢量。进一步的，步骤S5的具体步骤包括：根据约束条件，构造新的代价函数：s.ty(6+k)＝||y(3(k-1)+1:3k)||2,k＝1,2将(Gy-h)TW(Gy-h)转化为其中，Y＝yyT，由此可知，Y(6+k,6+k)＝trace{Y(3k-2:3k,3k-2:3k)}，其中，k＝1,2；由Y＝yyT，其中y＝[(u1-s1)T,(u2-s1)T,r1,1,r1,2]T，可知Y(7,8)＝r1,1r1,2。进一步的，步骤S6的具体步骤包括：利用已构造的新的代价函数进行优化求解时，把相应的等式约束转化凸半正定规划优化问题，待求解的公式为：其约束条件为：s.tY(6+k,6+k)＝trace{Y(3k-2:3k,3k-2:3k)}更进一步的，步骤S7的具体步骤包括：对Y进行特征值分解可得：其中，λi为Y的特征值，i＝1,...,r，qi为相应的特征向量，设i＝1,...,r已按从大到小排列，即λ1≥λ2≥...≥λr＞0，根据Y的定义可得信号源位置估计为：本专利技术的有益效果在于：1)、本专利技术首先构建时差定位方程，再对信号源位置即待定位的目标位置进行最大似然估计，引入辅助矢量，将距离差定位方程转化为约束最小二乘问题，再对所述辅助矢量进行加权最小二乘求解，再利用加权最小二乘解算的辅助矢量初始估计松弛等式约束，构造新的代价函数，利用凸半正定规划优化求解辅助矢量和辅助矢量转置的变量的值，并通过特征值分解得到辅助矢量的值，最后根据求得的辅助矢量的值与信号源位置之间的关系获取信号源的位置信息。本专利技术针对时差定位方法中出现的非凸优化问题所提出的，将高度非线性的时差定位问题转化为约束加权最小二乘估计问题，然后在此基础上通过半正定松弛转化为凸的半正定规划问题，进行优化求解，避免了传统迭代算法中出现的局部收敛和求解发散的问题，提高了定位的精度。2)、本方法在测量误差较大时仍能对信号源位置进行较好的估计，具有较强的稳健性。附图说明图1为本专利技术的流程图；图2为利用本半正定松弛时差定位方法、Chan算法、泰勒级数法对信号源u1和u2位置估计随变化时的均方根误差的统计结果示意图；图3为利用本半正定松弛时差定位方法、Chan算法、MDS算法、泰勒级数法对信号源u3位置估计随变化时的均方根误差的统计结果示意图。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。用于分布式多点定位监视系统的半正定松弛时差定位方法是时差定位方法中非常重要的一个步骤，时差定位方法的具体步骤包括：首先分布式多点定位监视系统中目标接收站主站和接收站辅站获取目标应答信号，对所述目标应答信号经本文档来自技高网...

【技术保护点】
用于分布式多点定位监视系统的半正定松弛时差定位方法，其特征在于，包括以下步骤：S1、构建时差定位方程；S2、对信号源位置即待定位的目标位置进行最大似然估计；S3、引入辅助矢量，将距离差定位方程转化为约束最小二乘问题；S4、对所述辅助矢量进行加权最小二乘求解；S5、利用加权最小二乘解算的辅助矢量初始估计松弛等式约束，构造新的代价函数；S6、利用凸半正定规划优化求解辅助矢量和辅助矢量转置的变量的值，并通过特征值分解得到辅助矢量的值；S7、根据求得的辅助矢量的值与信号源位置之间的关系获取信号源的位置信息。

【技术特征摘要】
1.用于分布式多点定位监视系统的半正定松弛时差定位方法，其特征在于，包括以下步骤：S1、构建时差定位方程；S2、对信号源位置即待定位的目标位置进行最大似然估计；S3、引入辅助矢量，将距离差定位方程转化为约束最小二乘问题；S4、对所述辅助矢量进行加权最小二乘求解；S5、利用加权最小二乘解算的辅助矢量初始估计松弛等式约束，构造新的代价函数；S6、利用凸半正定规划优化求解辅助矢量和辅助矢量转置的变量的值，并通过特征值分解得到辅助矢量的值；S7、根据求得的辅助矢量的值与信号源位置之间的关系获取信号源的位置信息。2.如权利要求1所述的用于分布式多点定位监视系统的半正定松弛时差定位方法，其特征在于，步骤S1的具体步骤包括：利用接收站主站位置信息测量参数数据建立相应的观测模型，确定的接收站主站的站址信息，利用目标到接收站主站的距离、目标到接收站辅站的距离差构建距离差定位方程：其中，为到达接收站主站和接收站辅站时间差的测量值，c为光速，ri,k-r1,k为信号源发射的应答信号从第k个目标位置uk到接收站si和接收站主站s1之间的真实时间差，ni1,k为测量误差，i为接收站辅站的索引，k为目标索引。3.如权利要求2所述的用于分布式多点定位监视系统的半正定松弛时差定位方法，其特征在于，步骤S2的具体步骤包括：利用目标到接收站主站和第i个接收站辅站的距离差测量误差的协方差矩阵表示目标位置的最大似然估计：其中，J(u)为信号源位置u的高度非线性、非凸代价函数，为所构成的距离差估计组成的向量，f(u)为信号源发射的应答信号从第k个目标位置uk到接收站si和接收站主站s1之间的真实时间差组成的向量，Q为测量误差ni1,k服从均值为零的协方差矩阵，T为矩阵的转置，为使得代价函数J(u)最小时所需信号源位置u的值。4.如权利要求3所述的用于分布式多点定位监视系统的半正定松弛时差定位方法，其特征在于，步骤S3的具体步骤包括：引入辅助矢量即为目标与接收站主站的位置和距离，转化为约束最小二乘问题，并在此基础上转换为凸的半正定规划问题，将原代价函数整理得到的矩阵如下：其中，ηk＝Bknk，Bk＝diag{r2,k,…,rN,k}，k＝1,2，R＝[s2-s1,...,sN...

【专利技术属性】
技术研发人员：杨琳，徐瑾，刘锐，马磊，崔扬，
申请(专利权)人：安徽四创电子股份有限公司，
类型：发明
国别省市：安徽,34