用于分布式多点定位监视系统的时差定位方法 |
|||||||
申请号 | CN201611204604.3 | 申请日 | 2016-12-23 | 公开(公告)号 | CN106842118A | 公开(公告)日 | 2017-06-13 |
申请人 | 安徽四创电子股份有限公司; | 发明人 | 杨琳; 徐瑾; 刘锐; 马磊; 李朋; | ||||
摘要 | 本 发明 属于分布式多点 定位 监视系统的实时定位技术领域,特别涉及一种用于分布式多点定位监视系统的时差定位方法。本发明根据GDOP时差定位 精度 选择定位方法;依据定位方法确定接收站主站,建立观测模型;利用半正定松弛的时差定位 算法 解算出目标的 位置 信息,并将定位结果通过终端显示 监控系统 呈现给用户,从而提高机场管制员对机场场面飞机和引导车监控管制能 力 。本发明的所有计算分析过程都是自发的,最大限度的减少了人为的因素对飞机和引导车定位结果的影响,而且所有的 数据处理 过程都是并行的,极大地提高了数据处理的效率,使用户执行操作时能最快的获得分析的结果。本发明定位准确,在测量误差较大时仍能对 信号 源的位置进行较好的估计。 | ||||||
权利要求 | 1.用于分布式多点定位监视系统的时差定位方法,其特征在于,包括以下步骤: |
||||||
说明书全文 | 用于分布式多点定位监视系统的时差定位方法技术领域[0001] 本发明属于分布式多点定位监视系统的实时定位技术领域,特别涉及一种用于分布式多点定位监视系统的时差定位方法。 背景技术[0002] 随着我国经济的飞速发展,航班量也在飞速增长,对空管系统管制员的指挥和设施设备的维护提出了更高的要求,机场场面监视雷达是一次雷达,主要用于监视机场场面的飞机及车辆,利用目标对电磁波的自身辐射或反射特性发现目标。设备维护人员在维护机场场面监视雷达的过程中包括如下三方面问题:飞机目标分裂、存在一定的盲区、容易出现假目标。这些问题的存在扰乱了管制员正常的指挥工作,因此,分布式多点定位系统的引入必要而且必须,这项技术是机场场面监视雷达的升级和补充,并且能够通过增加接收站或者改变接收站的布局实现系统的扩展。 [0003] 传统的时差定位方法在进行定位精度解算时,考虑到各个目标接收站误差元素之间是相关的,降低了定位精度解算的效率,提高了定位精度解算复杂程度,而且通常使用泰勒级数算法、Chan算法、MDS算法来确定目标的位置信息,使用泰勒级数算法对迭代初始值的选取有要求,会出现局部收敛或者发散的现象;Chan算法和MDS算法在测量误差较大时导致对信号源的位置判断不准确,从而导致时差定位方法经常出现高度非线性的问题,使得定位不准确。 发明内容[0004] 本发明为了克服上述现有技术的不足,提供了用于分布式多点定位监视系统的时差定位方法,本时差定位方法极大地降低了定位精度解算方法的复杂性,提高了定位精度解算的效率,而且对迭代初始值的选取没有要求,能够收敛到全局极小点,不会出现局部收敛或者发散的现象;在测量误差较大时仍能对信号源的位置进行较好的估计,定位准确。 [0005] 为实现上述目的,本发明采用了以下技术措施: [0006] 用于分布式多点定位监视系统的时差定位方法,包括以下步骤: [0007] S1、分布式多点定位监视系统中目标接收站主站和接收站辅站获取目标应答信号,对所述目标应答信号经过该分布式多点定位监视系统处理得到解析后的测量参数数据; [0009] S3、对配对处理后的测量参数数据进行主站位置选择,得到待定位区域的定位精度; [0010] S4、根据目标的测量参数数据以及待定位区域的定位精度,选择目标的接收站主站,确定时差定位模型; [0011] S5、利用接收站主站和接收站辅站位置信息测量参数数据,通过半正定松弛的时差定位算法解算出目标的位置信息; [0012] S6、将所述目标的位置信息送至终端显示。 [0013] 优选的,所述测量参数数据包括目标应答码、置信度、目标应答方式、目标应答时间、目标的幅度值。 [0014] 优选的,步骤S3的具体步骤包括: [0015] S31、确定目标的接收站主站的位置信息; [0016] 获取配对处理后的目标的接收站的测量数据,选取其中一个接收站作为目标的接收站主站位置,建立观测模型计算该观测模型条件下定位区域的GDOP值,实现对目标的接收站主站的位置信息确定;通过选取不同的接收站作为目标的接收站主站,分别计算不同观测模型条件下定位区域的GDOP;通过选取不同观测模型条件下GDOP最小值即为时差定位精度最小值,实现对分布式多点定位监视系统时差定位模型的构建; [0017] S32、对待定位区域的位置信息求微分; [0018] 通过目标的接收站主站的区域覆盖图、目标的接收站辅站的区域覆盖图,获取可定位区域的位置信息,利用所述可定位区域的位置信息以及到达目标的第i个接收站辅站与目标的接收站主站之间的时间差,通过时差定位方程,对待定位区域的位置信息求微分: [0019] c·Δti=c·(ti-t0)=ri-r0(i=1,2) [0020] 其中,c为光速、Δti为到达目标的第i个接收站辅站与目标的接收站主站之间时间差、ti为到达目标的第i个接收站辅站的时间、t0为到达目标的接收站主站时间、ri为到达目标的第i个接收站辅站的距离、r0为到达目标的接收站主站距离、i为到达目标的接收站辅站的编号; [0021] S33、求出所述待定位区域内目标与目标的接收站主站位置、目标与目标的接收站辅站位置的相关系数矩阵; [0022] 通过所述待定位区域内目标到目标的各接收站辅站与目标到目标的接收站主站的位置差,得出待定位区域内目标与站址位置的相关系数矩阵: [0023] [0024] x,y均为目标的位置信息、x0,y0均为目标的接收站主站位置信息、xi、yi均为目标的第i个接收站辅站位置信息、r0为目标到目标的接收站主站距离、ri为目标到目标的第i个接收站辅站距离, [0025] S34、利用伪逆法计算所述待定位区域的误差估计值; [0026] 定位误差方程式为: [0027] 其中, 为定位误差在x方向上定位方差、 为定位误差在y方向上定位方差、i为目标的第i个接收站辅站、j为目标的第j个接收站辅站, 其中,ηij为Δti与Δtj之间的相关系数, 为第i站的时间测量误差的标准差, 为第j站的时间测量误差的标准差,σs为站址测量误差的标准差; [0028] S35、求出待定位区域的定位精度; [0029] 根据时间测量误差的相关系数ηij,各站址测量误差的标准差σs,利用线性组合求出σij,然后求得x方向上定位方差 y方向上定位方差 通过 公式,计算得到待定位区域的定位精度GDOP。 [0030] 优选的,步骤S5的具体步骤包括: [0031] S51、构建时差定位方程; [0032] 利用接收站主站位置信息测量参数数据建立相应的观测模型,确定的接收站主站的站址信息,利用目标到接收站主站的距离、目标到接收站辅站的距离差构建距离差定位方程: [0033] 其中, 为到达接收站主站和接收站辅站时间差的测量值,c为光速,ri,k-r1,k为信号源发射的应答信号从第k个目标位置uk到接收站si和接收站主站s1之间的真实时间差,ni1,k为测量误差,i为接收站辅站的索引,k为目标索引; [0034] S52、对信号源位置即待定位的目标位置进行最大似然估计; [0035] 利用目标到接收站主站和第i个接收站辅站的距离差测量误差的协方差矩阵表示目标位置的最大似然估计: [0036] [0037] 其中,J(u)为信号源位置u的高度非线性、非凸代价函数,为 所构成的距离差估计组成的向量,f(u)为信号源发射的应答信号从第k个目标位置uk到接收站si和接收站主站s1之间的真实时间差组成的向量,Q为测量误差ni1,k服从均值为零的协方差矩阵,T为矩阵的转置, 为使得代价函数J(u)最小时所需信号源位置u的值; [0038] S53、引入辅助矢量,将距离差定位方程转化为约束最小二乘问题; [0040] [0041] 其中, η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的到达接收站主站和接收站辅站距离差估计组成的向量; [0042] S54、对所述辅助矢量进行加权最小二乘求解; [0043] 对原代价函数整理得到的矩阵进行加权最小二乘求解,得到的辅助矢量估计值如下: [0044] [0045] 其中,加权矩阵W=E[ηηT]-1=(BQBT)-1,B=diag{B1,B2}表示以B1和B2为对角块的矩阵,Q=E[nnT]为噪声协方差矩阵,n为到达接收站主站和接收站辅站距离差测量误差的矩阵形式, Gy-h=η; [0046] S55、利用加权最小二乘解算的辅助矢量初始估计松弛等式约束,构造新的代价函数; [0047] 根据约束条件,构造新的代价函数: [0048] [0049] s.t y(6+k)=||y(3(k-1)+1:3k)||2,k=1,2 [0050] 将(Gy-h)TW(Gy-h)转化为 [0051] 其中,Y=yyT, [0052] 由此可知,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; [0053] S56、利用凸半正定规划优化求解辅助矢量和辅助矢量转置的变量的值,并通过特征值分解得到辅助矢量的值; [0054] 利用已构造的新的代价函数进行优化求解时,把相应的等式约束转化凸半正定规划优化问题, [0055] 待求解的公式为: 其约束条件为: [0056] s.t Y(6+k,6+k)=trace{Y(3k-2:3k,3k-2:3k)} [0057] [0058] S57、根据求得的辅助矢量的值与目标位置之间的关系获取目标的位置信息; [0059] 对Y进行特征值分解可得: [0060] 其中,λi为Y的特征值,i=1,...,r,qi为相应的特征向量,设i=1,...,r已按从大到小排列,即λ1≥λ2≥...≥λr>0,根据Y的定义可得信号源位置 估计为: [0061] 本发明的有益效果在于: [0062] 1)、本发明根据GDOP时差定位精度选择定位方法;依据定位方法确定接收站主站,建立观测模型;利用半正定松弛的时差定位算法解算出目标的位置信息,并将定位结果通过终端显示监控系统呈现给用户,从而提高机场管制员对机场场面飞机和引导车监控管制能力。本发明的所有计算分析过程都是自发的,最大限度的减少了人为的因素对飞机和引导车定位结果的影响,而且所有的数据处理过程都是并行的,极大地提高了数据处理的效率,使用户执行操作时能最快的获得分析的结果。 [0063] 2)、确定待定位区域的定位精度是利用多个目标的接收站的位置信息和内部噪声、环境噪声和时间测量的量化噪声信息对目标的定位估计值进行定位精度解算,通过不同的组合和不同的布站方式,得到测量目标区域覆盖的定位精度GDOP,从而与不同布站方式下的GDOP比对,选择最优布站方式,为后续时差定位解算模型的构建提供依据。 [0064] 3)、本发明将高度非线性的时差定位问题转化为约束加权最小二乘估计问题,然后在此基础上通过半正定松弛转化为凸的半正定规划问题,进行优化求解,避免了传统迭代算法中出现的局部收敛和求解发散的问题,大大地提高了定位的精度。附图说明 [0065] 图1为本发明的流程图; [0066] 图2为本发明的步骤S3的工作流程图; [0067] 图3为本发明的步骤S5的工作流程图; [0068] 图4为本发明的Y形布站仿真图; [0069] 图5为本发明的倒三角形布站仿真图; [0070] 图6为本发明的菱形布站仿真图; [0071] 图7为利用本半正定松弛时差定位方法、Chan算法、泰勒级数法对信号源u1和u2位置估计随 变化时的均方根误差的统计结果示意图; [0072] 图8为利用本半正定松弛时差定位方法、Chan算法、MDS算法、泰勒级数法对信号源u3位置估计随 变化时的均方根误差的统计结果示意图。 具体实施方式[0073] 下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。 [0074] 如图1所示,用于分布式多点定位监视系统的时差定位方法,包括以下步骤: [0075] S1、分布式多点定位监视系统中目标接收站主站和接收站辅站获取目标应答信号,对所述目标应答信号经过该分布式多点定位监视系统处理得到解析后的测量参数数据; [0076] S2、根据所述目标接收站主站和接收站辅站获取所述测量参数数据,对测量参数数据中相同目标应答码和大于50%的置信度进行配对处理,获取目标在同一时刻发出的测量参数数据; [0077] S3、对配对处理后的测量参数数据进行主站位置选择,得到待定位区域的定位精度; [0078] S4、根据目标的测量参数数据以及待定位区域的定位精度,选择目标的接收站主站,确定时差定位模型; [0079] S5、利用接收站主站和接收站辅站位置信息测量参数数据,通过半正定松弛的时差定位算法解算出目标的位置信息; [0080] S6、将所述目标的位置信息送至终端显示。 [0081] 所述目标即为待定位目标。 [0082] 所述测量参数数据包括目标应答码、置信度、目标应答方式、目标应答时间、目标的幅度值。 [0083] 如图2所示,对配对处理后的测量参数数据进行主站位置选择,得到待定位区域的定位精度的具体操作步骤包括: [0084] S31、确定目标的接收站主站的位置信息; [0085] 获取配对处理后的目标的接收站的测量数据,选取其中一个接收站作为目标的接收站主站位置,建立观测模型计算该观测模型条件下定位区域的GDOP值,实现对目标的接收站主站的位置信息确定;通过选取不同的接收站作为目标的接收站主站,分别计算不同观测模型条件下定位区域的GDOP;通过选取不同观测模型条件下GDOP最小值即为时差定位精度最小值,实现对分布式多点定位监视系统时差定位模型的构建; [0086] S32、对待定位区域的位置信息求微分; [0087] 通过目标的接收站主站的区域覆盖图、目标的接收站辅站的区域覆盖图,获取可定位区域的位置信息,利用所述可定位区域的位置信息以及到达目标的第i个接收站辅站与目标的接收站主站之间的时间差,通过时差定位方程,对待定位区域的位置信息求微分: [0088] c·Δti=c·(ti-t0)=ri-r0(i=1,2) [0089] 其中,c为光速、Δti为到达目标的第i个接收站辅站与目标的接收站主站之间时间差、ti为到达目标的第i个接收站辅站的时间、t0为到达目标的接收站主站时间、ri为到达目标的第i个接收站辅站的距离、r0为到达目标的接收站主站距离、i为到达目标的接收站辅站的编号; [0090] S33、求出所述待定位区域内目标与目标的接收站主站位置、目标与目标的接收站辅站位置的相关系数矩阵; [0091] 通过所述待定位区域内目标到目标的各接收站辅站与目标到目标的接收站主站的位置差,得出待定位区域内目标与站址位置的相关系数矩阵: [0092] [0093] x,y均为目标的位置信息、x0,y0均为目标的接收站主站位置信息、xi、yi均为目标的第i个接收站辅站位置信息、r0为目标到目标的接收站主站距离、ri为目标到目标的第i个接收站辅站距离, [0094] S34、利用伪逆法计算所述待定位区域的误差估计值; [0095] 定位误差方程式为: [0096] 其中, 为定位误差在x方向上定位方差、 为定位误差在y方向上定位方差、i为目标的第i个接收站辅站、j为目标的第j个接收站辅站,其中,ηij为Δti与Δtj之间的相关系数, 为第i站的时间测量误差的标准差, 为第j站的时间测量误差的标准差,σs为站址测量误差的标准差; [0097] S35、求出待定位区域的定位精度; [0098] 根据时间测量误差的相关系数ηij,各站址测量误差的标准差σs,利用线性组合求出σij,然后求得x方向上定位方差 y方向上定位方差 通过 公式,计算得到待定位区域的定位精度GDOP。 [0099] 当主站的坐标取(0,0)时,三个辅站的坐标分别取不同的值时,如表1所示,形成的仿真图分别为Y形布站、倒三角布站、菱形布站。 [0100] 表1: [0101] [0102] [0103] 从图4~图6可知,对于同一种布站形式来说,定位精度GDOP随着目标和接收站的分布区域之间的距离的增加而降低,对于不同的布站方式,全方位内定位误差有所不同,用于分布式多点定位监视系统的定位精度方法,可以充分考虑Y形、倒三角形、菱形布站方法对不同区域定位精度的不同,来进行目标的接收站主站和辅站位置的选择。力求目标在不同的区域里,选择最优布站方案,实现目标的高精度定位。 [0104] 如图3所示,利用接收站主站和接收站辅站位置信息测量参数数据,通过半正定松弛的时差定位算法解算出目标的位置信息的具体操作步骤为: [0105] S51、构建时差定位方程; [0106] 利用接收站主站位置信息测量参数数据建立相应的观测模型,确定的接收站主站的站址信息,利用目标到接收站主站的距离、目标到接收站辅站的距离差构建距离差定位方程: [0107] 其中, 为到达接收站主站和接收站辅站时间差的测量值,c为光速,ri,k-r1,k为信号源发射的应答信号从第k个目标位置uk到接收站si和接收站主站s1之间的真实时间差,ni1,k为测量误差,i为接收站辅站的索引,k为目标索引。 [0108] S52、对信号源位置即待定位的目标位置进行最大似然估计; [0109] 利用目标到接收站主站和第i个接收站辅站的距离差测量误差的协方差矩阵表示目标位置的最大似然估计: [0110] [0111] 其中,J(u)为信号源位置u的高度非线性、非凸代价函数,为 所构成的距离差估计组成的向量,f(u)为信号源发射的应答信号从第k个目标位置uk到接收站si和接收站主站s1之间的真实时间差组成的向量,Q为测量误差ni1,k服从均值为零的协方差矩阵,T为矩阵的转置, 为使得代价函数J(u)最小时所需信号源位置u的值; [0112] S53、引入辅助矢量,将距离差定位方程转化为约束最小二乘问题; [0113] 引入辅助矢量即为目标与接收站主站的位置和距离,转化为约束最小二乘问题,并在此基础上转换为凸的半正定规划问题,将原代价函数整理得到的矩阵如下: [0114] [0115] 其中, η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的到达接收站主站和接收站辅站距离差估计组成的向量; [0116] S54、对所述辅助矢量进行加权最小二乘求解; [0117] 对原代价函数整理得到的矩阵进行加权最小二乘求解,得到的辅助矢量估计值如下: [0118] [0119] 其中,加权矩阵W=E[ηηT]-1=(BQBT)-1,B=diag{B1,B2}表示以B1和B2为对角块的矩阵,Q=E[nnT]为噪声协方差矩阵,n为到达接收站主站和接收站辅站距离差测量误差的矩阵形式, Gy-h=η,y为辅助矢量; [0120] S55、利用加权最小二乘解算的辅助矢量初始估计松弛等式约束,构造新的代价函数; [0121] 根据约束条件,构造新的代价函数: [0122] [0123] s.t y(6+k)=||y(3(k-1)+1:3k)||2,k=1,2 [0124] 将(Gy-h)TW(Gy-h)转化为 [0125] 其中,Y=yyT, [0126] 由此可知,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; [0127] S56、利用凸半正定规划优化求解辅助矢量和辅助矢量转置的变量的值,并通过特征值分解得到辅助矢量的值; [0128] 利用已构造的新的代价函数进行优化求解时,把相应的等式约束转化凸半正定规划优化问题, [0129] 待求解的公式为: 其约束条件为: [0130] s.t Y(6+k,6+k)=trace{Y(3k-2:3k,3k-2:3k)} [0131] [0132] S57、根据求得的辅助矢量的值与目标位置之间的关系获取目标的位置信息; [0133] 对Y进行特征值分解可得: [0134] [0135] 其中,λi为Y的特征值,i=1,...,r,qi为相应的特征向量,设i=1,...,r已按从大到小排列,即λ1≥λ2≥...≥λr>0,根据Y的定义可得信号源位置 估计为: [0136] 为了检验本半正定松弛时差定位方法对信号源位置估计的性能,我们将本半正定松弛时差定位方法同传统的Chan算法、泰勒级数法、MDS算法及克拉美罗界的仿真结果进行比较,5个接收站位置坐标分别为:s1=[300,100,150]T,s2=[400,150,100]T,s3=[300,500,200]T,s4=[350,200,100]T,s5=[-100,-100,-100]T,信号源坐标单位均为米。假设各个信号源的TDOA测量值服从均值为零、方差为 的高斯分布,则Qk满足对角线元素为非对角线元素为 为了对比的公平,利用对所述辅助矢量进行加权最小二乘求解的具体步骤中的所得初始定位结果计算利用加权最小二乘解算的辅助矢量初始估计松弛等式约束,构造新的代价函数的具体步骤中的F和Chan算法、MDS算法的权值以及作为泰勒级数法的迭代初始值。三个信号源位置坐标分别为u1=[314,483,209]T,u2=[600,650,550]T和u3=[285,325,275]T。 [0137] 如图7所示,分别给出了本半正定松弛时差定位方法、Chan算法及泰勒级数法对信号源u1和u2位置估计随 变化时的均方根误差的统计结果,较低的一组曲线对应信号源u1,较高的一组对应u2。可以看出,在噪声功率较小时,三种方法对两信号源位置的估计均方根误差均十分接近克拉美罗界。随着测量误差的增大,各算法的均方根误差均有所增加;其中泰勒级数法最先偏离克拉美罗界,这是因为在测量误差较大时,由步骤S4求得的初始值偏离真实值较远,泰勒级数法在进行迭代求解时,很容易陷入局部极小点或者发散;而本半正定松弛时差定位方法对两信号源位置的估计均方根误差均明显小于Chan算法,具有更为准确的定位性能。 [0138] 如图8所示,分别给出了本半正定松弛时差定位方法、Chan算法、MDS算法及泰勒级数法对信号源u3位置估计随 变化时的均方根误差的统计结果。可以看出,在噪声功率较小时,几种方法对信号源位置估计的均方根误差均十分接近克拉美罗界。随着噪声功率的增加,各算法的估计均方根误差都会有所增加。在测量误差较大时,泰勒级数法的估计均方根误差出现急剧增加的现象;MDS算法偏离克拉美罗界要明显晚于Chan算法,但是在测量误差较大时,出现MDS算法的均方根误差大于Chan算法的现象;而本半正定松弛时差定位方法的估计均方根误差一直保持最小。 [0139] 本发明将定位结果通过终端显示监控系统呈现给用户,从而提高机场管制员对机场场面飞机和引导车监控管制能力。本发明的所有计算分析过程都是自发的,最大限度的减少了人为的因素对飞机和引导车定位结果的影响,而且所有的数据处理过程都是并行的,极大地提高了数据处理的效率,使用户执行操作时能最快的获得分析的结果,本发明定位准确,在测量误差较大时仍能对信号源的位置进行较好的估计。 |