专利名称：：基于tis-b的估计位置误差的计算方法和系统的制作方法
技术领域：
：本发明涉及航空监视技术，尤其涉及一种基于TIS-B的估计位置误差的计算方法和系统。
背景技术：
：随着人们生活水平的不断提高，越来越多的人选择飞机作为出行的主要交通工具，而安全性是乘客首要考虑的问题。因此，和地面交通管制类似，空中交通管制的目的是使航线上的飞机安全、高效、有计划地在空域中飞行，通过管制员对所管制空域内飞机的飞行动态进行实时监一见，以准确掌握空中的交通活动信息。现有技术中存在的监视技术包括雷达监视手段、自动相关监视(AutomaticDependentSurveillance;以下简称ADS)技术、广播式自动相关监视(AutomaticDependentSurvei1lance-Broadcast.，以下简称ADS-B)技术、广播式交通信息月良务(TrafficInformationService-Broadcast;以下简称TIS-B)技术等。其中，雷达监视手段采用询问/应答的方式对目标进行探测，但雷达系统自身存在诸多局限性，如雷达波束为直线传播，造成大量盲区，某些特殊地貌地区无法进行探测，雷达旋转周期在一定程度上限制了监视精度的提高，而且无法获得飞机的计划航路、速度等态势数据，对跟踪精度的提高和短期冲突告警(ShortTermConflictAlert;以下简称STCA)能力。ADS技术通过卫星数据链或甚高频空-地数据链，将机载导航系统获得的导航信息自动实时地发送到地面接收和处理系统，通过显示设备提供伪雷达画面，以供地面监视飞机的运行状态。ADS-B技术作为ADS技术中的一种，其以导航设备及其他机载设备产生的信息为数据源，采用先进的地5空/空空数据链作为通信手段，通过对外自动广播自身的状态参数，实现地面对飞机的实时监视；同时它也接收其他飞机的广播信息，达到飞机间的相互感知，实现对周边空域交通状况全面、详细的了解。但ADS-B技术不能作为唯一的监视手段，由于一旦导航系统出现问题，监视功能将丧失，造成无法估计的损失。因此，二次监3见雷达(SecondarySurveillanceRadar;以下简称SSR)技术和ADS-B技术的结合成为必然，以实现资源的最大化利用，TIS-B技术则用于解决在监视系统更新换代的过程中ADS-B与雷达监视手段的兼容性问题。而ADS-B与TIS-B的数据源的定位原理存在本质的区别，ADS-B的监视信息来源于卫星定位导4元系统(GlobalPositioningSystem;以下简称GPS),在其监视信息中本身带有相应的表征精度的数据项，TIS-B的监视信息则来源于SSR，且SSR数据中没有相应的表征精度的数据项，已有的TIS-B系统利用标准差比的三次方的反比例函数与单一的数字补充项来计算数据精度，但该计算方法计算不高，使得得到的飞机位置存在较大误差。
发明内容本发明的目的在于提供一种基于TIS-B的估计位置误差的计算方法和系统，提供适用于TIS-B系统的估计位置误差计算方法，实现ADS-B系统与TIS-B系统的兼容性，提高数据精度。为了实现上述目的，本发明提供了一种基于TIS-B的估计位置误差的计算方法，包括根据目标位置测量点和系统误差建立估计位置误差圆，以所述目标位置测量点为中心，以随机误差分布所在的高斯椭圆的长轴、短轴为坐标轴建立直角坐标系；根据系统误差矢量和随机误差的高斯分布的协方差矩阵计算标准差值当所述系统误差矢量与所述高斯椭圆的长轴在同一条直线上时，根据所述高斯椭圆被包含在所述估计位置误差圆内的积分概率计算变量因子值&95;对目标位置的实测值和目标位置的计算值的差值进行拟合，根据阻尼最小二乘估计理论估计得到^f务正因子值M;根据所述标准差值cr^、所述变量因子值&M和所述修正因子值A"十算估计位置误差值EPU。本实施还提供了一种基于TIS-B的估计位置误差的计算系统，包括模型建立模块，用于根据目标位置测量点和系统误差建立估计位置误差圆，以所述目标位置测量点为中心，以随才几误差分布所在的高斯椭圆的长轴、短轴为坐标轴建立直角坐标系；标准差值计算模块，用于根据系统误差矢量和随机误差的高斯分布的协方差矩阵计算标准差值cr^;变量因子值计算模块，用于当所述系统误差矢量与所述高斯椭圆的长轴在同一条直线上时，根据所述高斯椭圆被包含在所述估计位置误差圓内的积分概率计算变量因子值&95;修正因子值计算4莫块，用于对目标位置的实测值和目标位置的计算值的差值进行拟合，根据阻尼最小二乘估计理论估计得到修正因子值M;估计位置误差值计算模块，用于根据所述标准差值计算模块计算的所述标准差值cr^、所述变量因子值计算模块计算的所述变量因子值^95和所述修正因子值计算模块计算的所述修正因子值M计算估计位置误差值EPU。本发明提供的一种基于TIS-B的估计位置误差的计算方法和系统，提供了适用于TIS-B系统的估计位置误差计算方法，实现了ADS-B系统与TIS-B系统的兼容性，提高了数据精度。图1为本发明基于TIS-B的估计位置误差的计算方法实施例的流程图；图2为本发明基于TIS-B的估计位置误差的计算方法实施例中的NACp扩展模型示意图3为本发明基于TIS-B的估计位置误差的计算系统实施例的结构图。具体实施例方式下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。图1为本发明基于TIS-B的估计位置误差的计算方法实施例的流程图，如图1所示，本实施例提供了一种基于TIS-B的估计位置误差的计算方法，本实施例提供的基于TIS-B的估计位置误差的计算方法为基于TIS-B系统的飞机位置误差的计算方法，TIS-B系统用来解决ADS-B系统与雷达监视手段的兼容性问题，TIS-B系统对雷达监视信息进行处理后广播发送给其服务范围内的ADS-B用户，从而实现ADS-B用户对周围非ADS-B用户的监控。尽管TIS-B系统可以实现雷达监视手段与ADS-B系统的结合，而基于TIS-B系统本身的特殊性，现有的基于GPS的ADS-B航迹估计位置误差(EstimatePositionUncertainty;以下简称EPU)或/ZFOM鄉的量化算法不适用于TIS-B系统，因此，本实施例提供的方法即为适用于TIS-B系统的估计位置误差的计算方法。本实施例的基于TIS-B的估计位置误差的计算方法具体包括如下步骤步骤101,根据目标位置测量点和系统误差建立估计位置误差圆，以目标位置测量点为中心，以随机误差分布所在的高斯椭圆的长轴、短轴为坐标轴建立直角坐标系。在实施本发明的基于TIS-B的估计位置误差的计算方法时，先建立基于TIS-B系统的导航位置精度类另)J(NavigationAccuracyCategory-Position;以下简称NACp)扩展模型，该NACp扩展模型在传统的NACp模型的基础上，除考虑随机误差外，将系统误差也作为测量误差考虑在内。在以雷达为数据源的TIS-B的实际系统中，系统误差主要包括雷达传感器轴系误差、数据传递误差、动态滞后误差、数学变换的近似误差以及大气传输误差等。由于上述系统误差的存在，使得计算得到的EPU的值扩大，最终导致NACp级别的降低，因此，系统误差是不可忽略的。图2为本发明基于TIS-B的估计位置误差的计算方法实施例中的NACp扩展;f莫型示意图，如图2所示，布ii殳点A所在位置为目标位置测量点，根据目标位置测量点和系统误差建立估计位置误差EPU圆，该EPU圆即图2中的圆形虚线，该虚线圆的半径即为EPU。传统的NACp模型不考虑系统误差，则EPU圆的圆心即为目标位置测量点，而本NACp扩展模型则将系统误差考虑在内，因此，目标位置测量点则不再是EPU圆的圓心。而目标位置真实点在该NACp扩展模型中圆形虚线所包含的范围内，根据中心极限定理，TIS-B系统的随^/L误差分布为一个以目标位置测量点为中心的平面二维高斯椭圓。以目标位置测量点为中心，以随才几误差分布所在的高斯椭圆的长轴、短轴为坐标轴建立直角坐标系A-xy,将直角坐标系建立在目标位置测量点，而不是EPU圆的圆心，这样可以方便后续的分析，减少解决问题时的计算量。步骤102,根据系统误差矢量和随机误差的高斯分布的协方差矩阵计算标准差值cr^。完成上述NACp扩展模型的建立之后，在此NACp扩展才莫型下进行估计位置误差的计算，先根据系统误差矢量和随机误差的高斯分布的协方差矩阵V计算标准差值cj^。其中，系统误差矢量A为对目标位置测量点进行实时测量时存在的系统误差矢量，随才几误差的高斯分布的协方差矩阵V为在平面坐标系(x，y)中随机误差的高斯分布的平面坐标系的协方差矩阵。具体地，本步骤可以包括如下步骤首先，计算目标位置测量点的系统误差矢量A,可以采用下述公式(1)进行计算Z)=0/,d)r=(jc-m，_y—w)r(1)9其中，(x，y)为目标位置测量点的坐标值，(、，)为目标位置真实点的坐标值。其次，根据公式(1)计算的系统误差矢量久计算随机误差的高斯分布的协方差矩阵V,可以采用下述公式(2)进行计算<formula>formulaseeoriginaldocumentpage10</formula>(2)再次，根据上述公式(2)计算的协方差矩阵V获取标准差值"即矩阵V中第二行第二列的元素即为cx^的平方值。步骤103，当系统误差矢量与高斯椭圆的长轴在同一条直线上时，根据高斯椭圆被包含在估计位置误差圆内的积分概率计算变量因子值&95。在建立如图2所示的NACp扩展模型之后，假设EPU圆的半径为、o;,则包含有系统误差的NACp扩展模型可表示为如下形式<formula>formulaseeoriginaldocumentpage10</formula>(3)其中，A为高斯椭圆在EPU圆包含范围内的积分概率，、为一个变量因子，《为平面坐标系"力下x的方差，《为平面坐标系"，力下y的方差，积分区域S为EPU圓，且S由如下公式确定<formula>formulaseeoriginaldocumentpage10</formula>(4)町乂由公式(3)可知，EPU的计算转化成了当积分值&为0.95时EPU圆的半径的计算，即求A为0.95时&的值，因此将其设为&95。本步骤103可以具体包括如下步骤首先，当高斯椭圆被包含在估计位置误差圆内的积分概率为0.95时，获取变量因子的多项式表达式。由于等式(3)和(4)在积分区域S内进行积分没有解析解，因此通过本步骤103计算其数值积分，即假设系统误差矢量Dm与高斯椭圆的长轴在同一条直线上时，估算^95的最大化值。当系统误差矢量与高斯椭圓的长轴在同一条直线上时，根据公式(3)和(4)获取变量因子、的多项式表达式。其次，将变量因子的多项式表达式在点(1,0,0)附近进行泰勒级数展开，获取变量因子的泰勒级数展开式，如下式所示<formula>formulaseeoriginaldocumentpage11</formula>其中，W=/o^，"=仄为拉格朗日余项。具体地，在进行实验时，设初始值h(l，0卜2.4477，由非线性最优化理论的利弗博格-马夸特(Levenberg-Marquardt;以下简称L-M)算法可以得出，在最高阶r^4时即可满足精度要求，则泰勒级数展开式如下所示"0.95—<formula>formulaseeoriginaldocumentpage11</formula>公式(6)中各项对应的系数如下表1所示表1泰勒级数展开式中各项对应的系数<table>tableseeoriginaldocumentpage11</column></row><table>再次，根据变量因子的泰勒级数展开式计算变量因子值、95,根据上式(6)和表1中的各项对应的系数值便可计算得到变量因子值、95。步骤104,对目标位置的实测值和目标位置的计算值的差值进行拟合，根据阻尼最小二乘估计理论估计得到修正因子值M。由于实际中系统误差矢量A与高斯椭圆的长轴不一定在同一直线上，因此上述的估计方法存在一定程度的误差，、95的计算是基于总体测量数据拟合而成的，因此在某些区域内会产生突然上升或下降的现象，故引入修正值因子AA:,以得到更加精确的结果。本步骤104可以包括如下步骤首先，根据目标位置的实测值和目标位置的计算值计算差值，并对差值进行拟合，得到与差值相关的拟合公式。对获得的目标位置的多个实测值和目标位置的多个计算值进行差值计算，对差值进行拟合，得到与差值相关的拟合公式，如下式所示牟l=1,A:2)=胁"1+卢2'胁"2+/3sin(。&2+"3)+/4(7)其次，根据泰勒公式和阻尼最小二乘估计理论确定拟合公式(7)中的各项系数，假设&1=0_/0^=1，结果如下表2所示表2拟合公式中各项对应的系数<table>tableseeoriginaldocumentpage12</column></row><table>再次，根据高斯椭圆被包含在EPU圆内的积分概率的下降趋势设定如下修正函数(8):5(H):=l/exp{A:l-1}(8)。由于EPU圆所包含的范围内的积分概率会在某些区域出现突然上升或下降的现象，修正函数(8)的设定正是为了抑制这个现象。上述函数万(W)与W的指数函数成反比例函数的关系，用来描述估计位置误差EPU的下降趋势。最后，根据拟合公式(7)与修正函数(8)的乘积计算修正因子值Al将函数J(H=1，A:2)与B(A:1)相乘，确定4,正因子值M。步骤105，根据标准差值"变量因子值&M和修正因子值M计算估计位置误差值EPU。在通过上述步骤计算获得标准差值、变量因子值*。95和修正因子值M之后，将、；t,和A/t代入下述公式(9)中，便可根据计算得到估计位置误差值EPU:<formula>formulaseeoriginaldocumentpage13</formula>进一步地，在步骤105之后，还可以包括如下步骤对估计位置误差EPU进行误差级别量化。通常情况下，ADS-B系统通过NACp来表征系统获取的定位信息的4青度，也可以通过导4元位置误差类另1(NavigationUncertaintyCategory-Position;以下简称NUCp)来综合表征定位信息的精度和完好性。由于EPU表示的是目标位置的真实值在NACp扩展^f莫型中圆形虚线所包含的EPU圆的半径大小，在实际应用中，为了更好地方^f更用户，更直"l妄、迅速地得到结果，需要将EPU的误差级别进行量化。即将上述步骤得到EPU转换成NUCp或NACp的相应级别，通过筒单的数字级别来直》见地代表EPU的大小。具体地，使用者可以通过查询NUCp或NACp级别与所支持的飞行器间隔保证(AircraftSeperationAssurance;以下简称ASA)应用类型的对照表(如下表3和表4所示)，即可获取EPU相对应的NUCp或NACp级别。通常，NUCp的级别要求至少大于4，即NACp大于等于5时，再确定其是否满足所支持的特定ASA监视的应用需要。表3NACp级别与ASA应用类型对照表<table>tableseeoriginaldocumentpage13</column></row><table>表4NUCp级别与ASA应用类型对照表<table>tableseeoriginaldocumentpage14</column></row><table>本实施例提供了一种基于TIS-B的估计位置误差的计算方法，通过提供一种适用于TIS-B系统的估计位置误差计算方法，在符合关于ADS-B系统与TIS-B系统的数据精度参数定义的前提下，分析了针对SSR系统的EPU估算方法，提出了同时考虑系统误差和随机误差的条件下的NACp扩展模型，从而弥补了传统方法中未考虑系统误差存在的不足。同时，针对建立的NACp扩展模型，创造性地使用多元泰勒函数近似法得到EPU的估计算式，在算法的研究中使用了L-M最优化估计理论，在一定程度上提高了数据精度，确保了结果的准确度。本实施例实现了ADS-B系统与TIS-B系统的兼容性，大大提高了数据精度。图3为本发明基于TIS-B的估计位置误差的计算系统实施例的结构图，如图3所示，本实施例提供了一种基于TIS-B的估计位置误差的计算系统，基于TIS-B的估计位置误差的计算系统具体包括模型建立模块1、标准差值计算模块2、变量因子值计算模块3、修正因子值计算模块4和估计位置误差值计算模块5。其中，模型建立模块1用于根据目标位置测量点和系统误差建立估计位置误差圓，以目标位置测量点为中心，以随^/L误差分布所在的高斯椭圆的长轴、短轴为坐标轴建立直角坐标系。标准差值计算模块2用于根据系统误差矢量和随^/L误差的高斯分布的协方差矩阵计算标准差值CT^。变量因子值计算模块3用于当系统误差矢量与所述高斯椭圆的长轴在同一条直线上时，根据高斯椭圓被包含在所述估计位置误差圆内的积分概率计算变量因子值*。95。修正因子值计算模块4用于对目标位置的实测值和目标位置的计算值的差值进行拟合，根据阻尼最小二乘估计理论估计得到修正因子值AL估计位置误差值计算模块5用于根据标准差值计算模块2计算的标准差值cr^、变量因子值计算模块3计算的变量因子值&95和修正因子值计算模块4计算的修正因子值Ayt计算估计位置误差值EPU。具体地，标准差值计算模块2还可以包括系统误差矢量计算单元21、协方差矩阵计算单元22和标准差获取单元23。其中，系统误差矢量计算单元21用于计算目标位置测量点的系统误差矢量。协方差矩阵计算单元22用于根据系统误差矢量计算单元21计算的系统误差矢量计算随机误差的高斯分布的协方差矩阵。标准差获取单元23用于才艮据协方差矩阵计算单元22计算的协方差矩阵获取标准差值cx^。变量因子值计算模块3可以具体包括多项式获取单元31、展开式获取单元32和变量因子值计算单元33。其中，多项式获取单元31用于当高斯椭圓被包含在估计位置误差圆内的积分概率为0.95时，获取变量因子的多项式表达式。展开式获取单元32用于将多项式获取单元31获取的变量因子的多项式表达式进行泰勒级数展开，获取变量因子的泰勒级数展开式。变量因子值计算单元33用于根据展开式获取单元32获取的变量因子的泰勒级数展开式计算变量因子值&95。修正因子值计算模块4可以具体包括拟合单元41、系数确定单元42、设定单元43和修正因子值计算单元44。其中，拟合单元41用于根据目标位置的实测值和目标位置的计算值计算差值，并对差值进行拟合，得到与差值相关的拟合公式。系数确定单元42用于根据泰勒公式和阻尼最小二乘估计理论15确定拟合单元41得到的拟合公式中的各项系数。设定单元43用于根据高斯椭圆被包含在位置误差圆内的积分概率的下降趋势确定修正函数。修正因子值计算单元44用于根据拟合单元41得到的拟合公式与设定单元43得到的修正函数的乘积计算修正因子值。进一步地，本实施例提供的一种基于TIS-B的估计位置误差的计算系统还可以包括量化模块6，用于对估计位置误差值计算模块5计算得到的估计位置误差值EPU进行误差级别量化。本实施例提供了一种基于TIS-B的估计位置误差的计算系统，通过提供一种适用于TIS-B系统的估计位置误差计算系统，通过设置模型建立模块、标准差值计算模块、变量因子值计算模块、修正因子值计算模块和估计位置误差值计算模块，实现了ADS-B系统与TIS-B系统的兼容性，大大提高了数据精度。最后应说明的是以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解其依然可以对前述实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围。权利要求1、一种基于TIS-B的估计位置误差的计算方法，其特征在于，包括根据目标位置测量点和系统误差建立估计位置误差圆，以所述目标位置测量点为中心，以随机误差分布所在的高斯椭圆的长轴、短轴为坐标轴建立直角坐标系；根据系统误差矢量和随机误差的高斯分布的协方差矩阵计算标准差值σmy；当所述系统误差矢量与所述高斯椭圆的长轴在同一条直线上时，根据所述高斯椭圆被包含在所述估计位置误差圆内的积分概率计算变量因子值k0.95；对目标位置的实测值和目标位置的计算值的差值进行拟合，根据阻尼最小二乘估计理论估计得到修正因子值Δk；根据所述标准差值σmy、所述变量因子值k0.95和所述修正因子值Δk计算估计位置误差值EPU。2、根据权利要求1所述的方法，其特征在于，所述根据系统误差矢量和随机误差的高斯分布的协方差矩阵计算标准差值cr^包括计算所述目标位置测量点的系统误差矢量；根据所述系统误差矢量计算随机误差的高斯分布的协方差矩阵；根据所述协方差矩阵获取标准差值(；。3、根据权利要求l所述的方法，其特征在于，所述根据所述高斯椭圆被包含在所述估计位置误差圆内的积分概率计算变量因子值*。95包括当所述高斯椭圓被包含在所述估计位置误差圆内的积分概率为0.95时，获取变量因子的多项式表达式；将所述变量因子的多项式表达式进行泰勒级凄緣开，获取所述变量因子的泰勒级l緣开式；根据所述变量因子的泰勒级数展开式计算变量因子值yt。95。4、根据权利要求1所述的方法，其特征在于，所述对目标位置的实测值和目标位置的计算值的差值进行拟合，根据阻尼最小二乘估计理论估计得到修正因子值aa包括根据目标位置的实测值和目标位置的计算值计算差值，并对所述差值进行拟合，得到与所述差值相关的拟合公式；根据泰勒公式和阻尼最小二乘估计理论确定所述拟合公式中的各项系数；根据所述高斯椭圆被包含在所述位置误差圆内的积分概率的下降趋势设定修正函数；根据所述拟合公式与所述修正函数的乘积计算修正因子值m。5、根据权利要求1-4中任一项所述的方法，其特征在于，在所述计算估计位置误差值EPU之后，还包括对所述估计位置误差值EPU进行误差级别量化。6、一种基于TIS-B的估计位置误差的计算系统，其特征在于，包括模型建立模块，用于根据目标位置测量点和系统误差建立估计位置误差圆，以所述目标位置测量点为中心，以随机误差分布所在的高斯椭圆的长轴、短轴为坐标轴建立直角坐标系；标准差值计算模块，用于根据系统误差矢量和随机误差的高斯分布的协方差矩阵计算标准差值(；；变量因子值计算模块，用于当所述系统误差矢量与所述高斯椭圆的长轴在同一条直线上时，根据所述高斯椭圆被包含在所述估计位置误差圓内的积分概率计算变量因子值、95;修正因子值计算模块，用于对目标位置的实测值和目标位置的计算值的差值进行拟合，根据阻尼最小二乘估计理论估计得到修正因子值;估计位置误差值计算模块，用于根据所述标准差值计算模块计算的所述标准差值"所述变量因子值计算模块计算的所述变量因子值^95和所述修正因子值计算模块计算的所述修正因子值Ait计算估计位置误差值EPU。7、根据权利要求6所述的系统，其特征在于，所述标准差值计算模块包括系统误差矢量计算单元，用于计算所述目标位置测量点的系统误差矢量；协方差矩阵计算单元，用于根据所述系统误差矢量计算单元计算的所述系统误差矢量计算P逸机误差的高斯分布的协方差矩阵；标准差获取单元，用于根据所述协方差矩阵计算单元计算的所述协方差矩阵获取标准差值cr^。8、根据权利要求6所述的系统，其特征在于，所述变量因子值计算模块包括多项式获取单元，用于当所述高斯椭圆被包含在所述估计位置误差圆内的积分概率为0.95时，获取变量因子的多项式表达式；展开式获取单元，用于将所述多项式获取单元获取的所述变量因子的多项式表达式进行泰勒级M开，获取所述变量因子的泰勒级数展开式；变量因子值计算单元，用于根据所述展开式获取单元获取的所述变量因子的泰勒级数展开式计算变量因子值&95。9、根据权利要求6所述的系统，其特征在于，所述修正因子值计算模块包括拟合单元，用于根据目标位置的实测值和目标位置的计算值计算差值，并对所述差值进行拟合，得到与所述差值相关的拟合公式；系数确定单元，用于根据泰勒公式和阻尼最小二乘估计理论确定所述拟合单元得到的所述拟合7>式中的各项系数；设定单元，用于根据所述高斯椭圆被包含在所述位置误差圆内的积分概率的下降趋势设定修正函数；修正因子值计算单元，用于根据所述拟合单元得到的所述拟合公式与所述设定单元所述修正函数的乘积计算修正因子值M。10、根据权利要求6-9中任一项所述的系统，其特征在于，还包括量化模块，用于对所述估计位置误差值计算模块计算得到的所述估计位置误差值EPU进行误差级别量化。全文摘要本发明公开了一种基于TIS-B的估计位置误差的计算方法和系统，基于TIS-B的估计位置误差的计算方法，包括根据目标位置测量点和系统误差建立估计位置误差圆，以目标位置测量点为中心，以随机误差分布所在的高斯椭圆的长轴、短轴为坐标轴建立直角坐标系；根据系统误差矢量和随机误差的高斯分布的协方差矩阵计算标准差值σ<sub>my</sub>；当系统误差矢量与所述高斯椭圆的长轴在同一条直线上时，根据高斯椭圆被包含在估计位置误差圆内的积分概率计算变量因子值k<sub>0.95</sub>；对目标位置的实测值和目标位置的计算值的差值进行拟合，根据阻尼最小二乘估计理论估计得到修正因子值Δk；根据标准差值σ<sub>my</sub>、变量因子值k<sub>0.95</sub>和修正因子值Δk计算估计位置误差值EPU。本发明提高了数据精度。文档编号G01S5/02GK101533094SQ20091008255公开日2009年9月16日申请日期2009年4月23日优先权日2009年4月23日发明者伟刘,军张,张青竹,朱衍波,熙林,嘉高申请人:民航数据通信有限责任公司;北京航空航天大学