雷达相对系统误差估值有效性评价方法与流程

本发明属于多雷达数据融合
技术领域：
，具体涉及一种雷达相对系统误差估值有效性评价方法。
背景技术：
：雷达相对系统误差是指雷达在对目标进行测量时存在的相对于指定参照系的固定误差，包括雷达动态测量中的测向误差Δθ、测距误差Δρ、雷达站址标定误差(Δx,Δy)和雷达组网应用中存在的坐标转换误差等。在雷达联网观测中，动态目标的绝对位置多数情况下是不可知的，比较容易得到的是在同一时间段内不同雷达对同一飞行目标的离散观测值。而相对于网内其它雷达来说，总会有一部雷达的观测是精确的，此时可以以该雷达(命名为主站)的测量值作为对目标位置的真实描述，其它雷达(命名为次站)以此为参照，从而求得次站雷达相对于主站雷达的相对系统误差。多雷达组网应用中，不同雷达之间存在的相对系统误差会造成同一目标观测结果的空间分裂，严重时就会妨碍来自同一目标的航迹关联，而又造成对应不同目标的航迹的错误关联。因此研究系统误差的规律性，尽可能准确地定位、估计系统误差，对于提高雷达网目标状态估计和多雷达数据融合的准确性极为重要。估计、消除雷达测量中存在的系统误差有两种途径：设备校准和数据校准。通常，雷达完成阵地部署后进行的入网前阵地检飞校准多采用设备校准手段。设备校准一般依靠固定地物回波或飞机检飞方式进行，在工程实施上有可能引入新的误差且效费比低，对多雷达组网后的系统误差校准，情况就更为复杂。因此，必须寻求科学合理、易于实现的方法评价一部新入网雷达经过阵地检飞校准后系统误差估值的有效性。为了克服设备校准的局限性，人们提出了具有普遍意义的数据校准方法，如实时精度控制法、最小平方法、最大似然法、广义最小平方法和最小二乘相对系统误差修正法等等。它从数据处理的角度出发，在选定的参照系中对雷达测量数据进行分析，对其中存在的系统误差进行估计，并以此估计结果为依据，反过来对测量数据进行修正，以此达到同一坐标系内观测结果的一致性。这是一种直接、有效的校准方法，尤其在雷达组网后的误差修正中，显示出不可替代的优越性。但是如何评价某种数据校准方法的有效性，也成为一个亟待解决的问题。与雷达测量中存在的随机误差相比，系统误差不随时间改变，是相对固定的误差。对于一部已经完成阵地部署的雷达而言，在一段时期内，相对于指定参照系的系统误差不会发生较大改变。但是，雷达在不同气候环境(春、夏、秋、冬、云、雨等)、不同探测区域(按照方位、距离分格)的系统误差可能存在差异，在一个系统误差估值不能“包打天下”的情况下，就要具体问题具体分析，“因地制宜”地建立雷达相对系统误差估值表。在这一过程中，研究雷达相对系统误差估值有效性评价方法就显得十分必要了。传统的系统误差估值有效性评价多采用“取点评价法”，就是将同一时刻，雷达测量位置点与目标在参照系中的位置点(这两个点又称为时间配准点)之间的直线距离作为基本评价依据。为了提高评价的可信度，往往需要取多个时间配准点，将多个时间配准点修正前后距离的平均值之比作为最终的系统误差估值有效性评价指标。这种方法原理上比较科学，但在工程实现上，由于雷达异步工作，“多个时间配准点”在原始测量中几乎不存在，虽然采用卡尔曼滤波等方法可以外推得到心怡时刻的目标位置点，但外推过程又不可避免的引入了新的计算误差，因此取点评价法的可操作性和准确性不高。技术实现要素：(一)要解决的技术问题本发明要解决的技术问题是：如何提出一种雷达相对系统误差估值有效性评价方法，以提高雷达阵地检飞、数据校准方法选择和雷达相对系统误差估值表计算过程中对系统误差估值评价的科学性、准确性和可操作性，有助于提升雷达网目标状态估计和多雷达数据融合的质量。(二)技术方案为解决上述技术问题，本发明提供一种雷达相对系统误差估值有效性评价方法，所述评价方法包括如下步骤：步骤1：选取主站雷达、次站雷达对同一空中目标的一段直线航迹线观测数据；步骤2：对主站雷达观测数据进行坐标变换，得到一组对应的统一直角坐标{(Xtzi,Ytzi)，i＝1,2,…n}，n表示主站雷达观测数据的总点数；步骤3：分别对次站雷达观测数据、次站雷达系统误差修正数据进行坐标变换，得到两组统一直角坐标{(Xtcj,Ytcj)，j＝1,2,…m}和{(Xtxj,Ytxj)，j＝1,2,…m}，m表示次站雷达观测数据的总点数；步骤4：在统一直角坐标系中，使用加权直线航迹线参数估计模型，分别对主站雷达观测数据、次站雷达观测数据和次站雷达系统误差修正数据进行直线参数迭代估计，得到主站观测航迹线参数(kz,dz)、次站观测航迹线参数(kc,dc)和次站修正航迹线参数(kx,dx)；所述kz、kc、kx分别表示直线的斜率，dz、dc、dx分别表示直线的截距；步骤5：分别计算次站观测航迹线、次站修正航迹线与主站观测航迹线的相对熵JSc和JSx；相对熵越小，表明两条航迹线越相似；步骤6：次站修正航迹线系统误差估值有效性指标V表示为：(三)有益效果与现有技术相比较，本发明通过选取主、次站雷达对某个典型航路目标的一组观测数据，经中心统一直角坐标转换后，先后采用不加权直线航迹线参数估计模型和加权直线航迹线参数估计模型对两雷达观测航迹线、次站雷达系统误差修正航迹线进行直线参数迭代估计，通过计算次站观测航迹线系统误差修正前后与主站观测航迹线的相对熵，构建次站修正航迹线系统误差估值有效性指标，用以评价系统误差估值的有效性以及修正效果。本发明原理方法科学、实施步骤合理，与传统系统误差有效性评价方法相比，在提高准确性和可操作性的同时，简化了复杂度，方便了工程上的实现。本发明有助于提高雷达阵地检飞、数据校准方法选择和雷达相对系统误差估值表计算过程中对系统误差估值评价的科学性、准确性和可操作性，从而提升雷达网目标状态估计和多雷达数据融合的质量。本发明所提供的方法的时间复杂度和空间复杂度都很低，可操作性和实用性很强。附图说明图1为本发明技术方案中评价方法的主体流程示意图。图2为本发明实施例中主、次站雷达观测数据和两组次站雷达系统误差修正数据在统一直角坐标系中的显示图。具体实施方式为使本发明的目的、内容、和优点更加清楚，下面结合附图和实施例，对本发明的具体实施方式作进一步详细描述。本发明针对传统方法从局部(单点)入手带来的问题与不足，引入信息论中相对熵的概念，通过选取典型航路(目标保持一定高度沿直线飞行)上的一组主、从雷达测量数据，整体比较次站观测航迹线系统误差修正前后与主站观测航迹线的相似性来评价系统误差估值的有效性以及修正效果。这样在提高准确性和可操作性的同时，简化了评价方法的复杂程度，方便了工程上的实现。具体而言，为解决现有技术问题，本发明提供一种雷达相对系统误差估值有效性评价方法，所述方法应用于雷达阵地检飞、数据校准方法选择和雷达相对系统误差估值表计算过程中，属于雷达网目标状态估计和多雷达数据融合系统的前期数据准备工作范畴；所述评价方法包括如下步骤：步骤1：选取主站雷达、次站雷达对同一空中目标的一段直线航迹线观测数据；步骤2：对主站雷达观测数据进行坐标变换，得到一组对应的统一直角坐标{(Xtzi,Ytzi)，i＝1,2,…n}，n表示主站雷达观测数据的总点数；步骤3：分别对次站雷达(原始)观测数据、次站雷达系统误差修正数据进行坐标变换，得到两组统一直角坐标{(Xtcj,Ytcj)，j＝1,2,…m}和{(Xtxj,Ytxj)，j＝1,2,…m}，m表示次站雷达观测数据的总点数；步骤4：在统一直角坐标系中，使用加权直线航迹线参数估计模型，分别对主站雷达观测数据、次站雷达观测数据和次站雷达系统误差修正数据进行直线参数迭代估计，得到主站观测航迹线参数(kz,dz)、次站观测航迹线参数(kc,dc)和次站修正航迹线参数(kx,dx)；所述kz、kc、kx分别表示直线的斜率，dz、dc、dx分别表示直线的截距；步骤5：分别计算次站观测航迹线、次站修正航迹线与主站观测航迹线的相对熵JSc和JSx；相对熵越小，表明两条航迹线越相似；步骤6：次站修正航迹线系统误差估值有效性指标V表示为：V∈[0,1]时，V越大，表示系统误差估值越有效，航迹修正效果越好；V<0时，表示系统误差估值不仅完全无效，反而加剧了主、次站雷达观测航迹的分裂程度。其中，所述步骤1包括如下步骤：步骤1.1：选取空中目标处于直线航迹线时，主站雷达、次站雷达同时段上报的一段目标观测数据；步骤1.2：选取的主站雷达观测数据为：{(ρzi,θzi,tzi)，i＝1,2,…n}表示tzi时刻主站雷达测得的目标距离ρzi和方位θzi，其中θzi∈[0,360),单位为度，雷达正北方向为0度，正东为90度，顺时针为正；n表示主站雷达观测数据的总点数；步骤1.3：选取的次站雷达观测数据为：{(ρcj,θcj,tcj)，j＝1,2,…m}，表示tcj时刻次站雷达测得的目标距离ρcj和方位θcj，且|tc1-tz1|≤T，|tcm-tzn|≤T，m表示次站雷达观测数据的总点数。其中，所述步骤1.1中，每部雷达观测数据点数量一般不少于10点。其中，所述步骤1.1中，所述同时段是指：主站雷达观测数据、次站雷达观测数据的首点和末点时间差均不大于1个雷达探测周期T。其中，所述雷达探测周期T一般为10或20秒。其中，所述步骤2包括如下步骤：步骤2.1：将主站雷达观测数据中的(ρzi,θzi)，i＝1,2,…n，转换为以本站为中心的二维直角坐标(Xzi,Yzi)：Xzi＝ρzisinθziYzi＝ρzicosθzi步骤2.2：将(Xzi,Yzi)，i＝1,2,…n转换为中心统一直角坐标(Xtzi,Ytzi)：Xtzi＝Xzicosδxz-Yzisinδxz+XzxYtzi＝Xzisinδxz+Yzicosδxz+Yzx其中：(Xzx,Yzx)为主站雷达在中心统一直角坐标系中的坐标；δxz为主站雷达站址与直角坐标系中心点的经度差，单位为弧度。其中，所述步骤3包括如下步骤：步骤3.1：将次站雷达观测数据中的(ρcj,θcj)，j＝1,2,…m转换为以本站为中心的二维直角坐标(Xcj,Ycj)：Xcj＝ρcjsinθcjYcj＝ρcjcosθcj步骤3.2：将(Xcj,Ycj)，j＝1,2,…m转换为中心统一直角坐标(Xtcj,Ytcj)：Xtcj＝Xcjcosδxc-Ycjsinδxc+XcxYtcj＝Xcjsinδxc+Ycjcosδxc+Ycx其中：(Xcx,Ycx)为次站雷达在中心统一直角坐标系中的坐标；δxc为次站雷达站址与直角坐标系中心点的经度差，单位为弧度；步骤3.3：将次站雷达系统误差修正数据(ρxj,θxj)，j＝1,2,…m转换为以本站为中心的二维直角坐标(Xxj,Yxj)：Xxj＝ρxjsinθxjYxj＝ρxjcosθxj此处若假设次站系统误差估值为(Δρ,Δθ,Δx,Δy)，表示次站雷达在测距、测向和站址定位上存在的相对于主站的系统误差，则有如下关系存在：ρxj＝ρcj-Δρθxj＝θcj-Δθ步骤3.4：将(Xxj,Yxj)，j＝1,2,…m转换为中心统一直角坐标(Xtxj,Ytxj)：Xtxj＝Xxjcosδxc-Yxjsinδxc+Xcx-ΔxYtxj＝Xxjsinδxc+Yxjcosδxc+Ycx-Δy其中：(Xcx,Ycx)为次站雷达在中心统一直角坐标系中的坐标；δxc为次站雷达站址与直角坐标系中心点的经度差，单位为弧度；(Δx,Δy)表示次站相对于主站的站址定位系统误差。其中，所述步骤4包括如下步骤：步骤4.1：采用不加权直线航迹线参数估计模型粗略估计主站雷达观测到的直线航迹线y-k1x-d1＝0，其中k1为直线的斜率，d1为直线在Y轴上的截距，步骤4.1包括如下步骤：步骤4.1.1：用主站雷达所有观测点{(Xtzi,Ytzi),i＝1,2,...n}到该直线的距离li的平方和最小作为条件构造直线，计算在此条件下的这条直线的最佳参数(k′,d′)，其中k′为待求直线的斜率，d′为待求直线的截距；此处，主站雷达所有观测点{(Xtzi,Ytzi),i＝1,2,...n}简记为：{(xi,yi)，i＝1,2,…n}；对于(1)式，应有f′(k′,d′)分别对k′和d′求偏导数，并等于零，即有下式成立：令：a′＝c′0-a′1b′1，b′＝a′2-b′2-a1′2+b1′2，c′＝a′1b′1-c′0，则(2)式的解为：步骤4.1.2：按照距离最小原则确定方程的合理解；(k′1,d′1)和(k′2,d′2)都是方程(2)的实根，且k′1×k′2＝-1，即解得的两条直线相互垂直；按照观测点{(xi,yi)，i＝1,2,…n}到所求直线的距离的平方和最小原则，确定合理的直线参数值；该问题可简化为：计算测量点(x1,y1)分别到直线y＝k′1×x+d′1和直线y＝k′2×x+d′2的距离l′1,l′2：若|l′1|<|l′2|，则取(k′1,d′1)，否则取(k′2,d′2)作为所求直线的合理参数，记为(k1,d1)；步骤4.2：采用加权直线航迹线参数估计模型迭代估计主站雷达观测到的直线航迹线y-kzx-dz＝0，其中kz为直线的斜率，dz为直线的截距，步骤4.2包括如下步骤：步骤4.2.1：计算主站雷达各观测点(xi,yi)到直线y-k1x-d1＝0的距离之和；n为观测点数；步骤4.2.2：求主站雷达各观测点(xi,yi)到直线y-k(s)x-d(s)＝0的距离li；式中s表示迭代次数，n表示观测点数；s初始值为1，即：k(1)＝k1，d(1)＝d1；步骤4.2.3：求|li|的倒数；步骤4.2.4：求各点的权值vi；步骤4.2.5：求解加权直线航迹线参数估计模型；用主站雷达所有观测点{(xi,yi)，i＝1,2,…n}到该直线的加权距离(vi×li)的平方和最小作为条件构造直线，计算在此条件下的这条直线的最佳参数(k,d)；对于(3)式，应有f(k,d)分别对k和d求偏导数，并等于零，即有下式成立：令：c″＝c0+a1b1，则(4)式的解为：步骤4.2.6：按照距离最小原则确定方程的合理解；按照点{(xi,yi)，i＝1,2,…n}到所求直线的距离的平方和最小原则，确定合理的直线参数值；该问题可简化为：计算点(x1,y1)分别到直线y＝k1×x+d1和直线y＝k2×x+d2的距离l1,l2：若|l1|+|l2|>Lmin，Lmin初值为106，则输出(k(s-1),d(s-1))作为所求直线的合理参数，并简记为(kz,dz)，迭代过程结束；否则Lmin＝|l1|+|l2|；s值加1，若|l1|<|l2|，则取(k1,d1)，否则取(k2,d2)作为所求直线的合理参数，记为(k(s),d(s))，s表示迭代次数；步骤4.2.7：计算所有观测点到新直线y-k(s)x-d(s)＝0的加权距离之和f(s)(k(s),d(s))；式中s表示迭代次数，n表示观测点数；步骤4.2.8：判别是否为最佳解；若f(s)(k(s),d(s))≥f(s-1)(k(s-1),d(s-1))，则输出解(k(s-1),d(s-1))，并简记为(kz,dz)；否则重复步骤4.2.2至步骤4.2.7；式中s表示迭代次数；步骤4.3：参照步骤4.1的过程，采用不加权直线航迹线参数估计模型粗略估计次站雷达观测到的直线航迹线；步骤4.4：参照步骤4.2的过程，采用加权直线航迹线参数估计模型迭代估计次站雷达观测到的直线航迹线；步骤4.5：参照步骤4.1的过程，采用不加权直线航迹线参数估计模型粗略估计次站雷达数据进行系统误差修正后的直线航迹线y-kx1x-dx1＝0；步骤4.6：参照步骤4.2的过程，采用加权直线航迹线参数估计模型迭代估计次站雷达数据进行系统误差修正后的直线航迹线y-kxx-dx＝0。其中，所述步骤5包括如下步骤：步骤5.1：计算次站观测航迹线与主站观测航迹线的相对熵JSc；求解步骤如下:①计算f(xi′)和g(xi′)。f(xi′)＝kzxi′+dz，g(xi′)＝kcxi′+dc，i′＝1,2,...n+m。其中：xi′∈{(Xtzi|i＝1,2,...n)∪(Xtcj|j＝1,2,...m)}，表示主站观测点与次站观测点横坐标的并集，n和m分别表示主站雷达、次站雷达的观测点数；②计算KL(f(x)||g(x))和KL(g(x)||f(x))；③计算JSc；步骤5.2：计算次站修正航迹线与主站观测航迹线的相对熵JSx；求解步骤如下：①计算f′(x′i′)和g′(x′i′)。f′(x′i′)＝kzx′i′+dz，g′(x′i′)＝kxx′i′+dx，i′＝1,2,...n+m。其中：x′i′∈{(Xtzi|i＝1,2,...n)∪(Xtxj|j＝1,2,...m)}，表示主站观测点与次站修正点横坐标的并集，n和m分别表示主、次站雷达的观测点数；②计算KL′(f′(x)||g′(x))和KL′(g′(x)||f′(x))；③计算JSx；下面结合具体实施例来详细描述本发明。实施例1本实施例具体描述一种雷达相对系统误差估值有效性评价方法，并进行了模拟对比试验过程。所述评价方法应用于雷达阵地检飞、数据校准方法选择和雷达相对系统误差估值表计算过程中，属于雷达网目标状态估计和多雷达数据融合系统的前期数据准备工作范畴。所述评价方法包括如下步骤：步骤1：按照表1中的模拟雷达与目标参数，仿真产生主、次站雷达对同一空中目标的一段直线航迹线观测结果，见表2，表3。表1：基本参数表2：主站雷达观测数据序号时间t距离ρ(km)方位θ(度)10249.92179.96210247.25179.53320244.62178.92430241.87178.28540239.06177.85650236.60177.31760234.07176.59870231.59175.90980229.01175.541090226.55174.7811100223.96174.1012110221.38173.6213120218.97172.7614130216.63172.2415140214.20171.43表3：次站雷达观测数据步骤2：对表2的主站雷达观测数据进行坐标变换，得到统一直角坐标{(Xtzi,Ytzi)，i＝1,2,…15}。步骤2.1：将主站雷达观测数据点(ρzi,θzi)，i＝1,2,…15转换为以本站为中心的二维直角坐标(Xzi,Yzi)，计算公式为：Xzi＝ρzisinθziYzi＝ρzicosθzi计算结果见表4。表4：以主站为中心的二维直角坐标序号Xz(km)Yz(km)10.1803-249.915922.0088-247.240634.6122-244.581047.2612-241.762558.9662-238.8894611.0985-236.3348713.9182-233.6584816.5745-230.9913917.8050-228.31771020.5925-225.61651123.0172-222.77621224.5882-220.01041327.5795-217.22431429.2639-214.64781531.9351-211.8052步骤2.2：将(Xzi,Yzi)，i＝1,2,…11转换为中心统一直角坐标(Xxzi,Yxzi)，计算公式为：Xxzi＝Xzicosδxz-Yzisinδxz+XzxYxzi＝Xzisinδxz+Yzicosδxz+Yzx计算结果见表5。表5：主站测量点的中心统一直角坐标序号Xtz(km)Ytz(km)11500.18032300.084121502.00882302.759431504.61222305.419041507.26122308.237551508.96622311.110661511.09852313.665271513.91822316.341681516.57452319.008791517.80502321.6823101520.59252324.3835111523.01722327.2238121524.58822329.9896131527.57952332.7757141529.26392335.3522151531.93512338.1948步骤3：分别对次站雷达原始观测数据、次站系统误差修正数据进行坐标变换，得到两组统一直角坐标{(Xtcj,Ytcj)，j＝1,2,…15}和{(Xtxj,Ytxj)，j＝1,2,…15}。步骤3.1：将次站雷达观测数据点(ρcj,θcj)，j＝1,2,…15转换为以本站为中心的二维直角坐标(Xcj,Ycj)，计算公式为：Xcj＝ρcjsinθcjYcj＝ρcjcosθcj计算结果见表6。表6：以次站为中心的二维直角坐标步骤3.2：将(Xcj,Ycj)，j＝1,2,…15转换为中心统一直角坐标(Xtcj,Ytcj)，计算公式为：Xtcj＝Xcjcosδxc-Ycjsinδxc+XcxYtcj＝Xcjsinδxc+Ycjcosδxc+Ycx计算结果见表7。表7：次站测量点的中心统一直角坐标序号Xtc(km)Ytc(km)11492.93612292.391921497.89422294.575431499.40642297.932641500.76562298.880551504.39882304.213161508.11472306.197871508.00982308.603981512.57142311.267691514.03662313.9060101517.10702316.1215111520.08672317.5894121521.99672319.2396131525.33492324.1517141525.67042324.5845151530.83722328.8352步骤3.3：为了便于比较本发明的应用效果，我们使用表1中通过不同途径获得的需要评价的两组次站雷达系统误差估值1和估值2对表3中的次站雷达原始测量数据进行系统误差修正，得到次站雷达经测距、测向系统误差(Δρ,Δθ)修正后的数据集1和数据集2。见表8。其中，次站雷达系统误差修正数据1(ρ1xj,θ1xj)和数据2(ρ2xj,θ2xj)与次站雷达原始观测数据(ρcj,θcj)之间有如下关系：ρ1xj＝ρcj-Δρ1，θ1xj＝θcj-Δθ1ρ2xj＝ρcj-Δρ2，θ2xj＝θcj-Δθ2表8：次站雷达系统误差修正数据集1和数据集2然后，将表8中的次站雷达系统误差修正数据1和数据2转换为以本站为中心的二维直角坐标(X1xj,Y1xj)和(X2xj,Y2xj)：X1xj＝ρ1xjsinθ1xj，Y1xj＝ρ1xjcosθ1xjX2xj＝ρ2xjsinθ2xj，Y2xj＝ρ2xjcosθ2xj计算结果见表9。表9：次站雷达系统误差修正数据的站直角坐标步骤3.4：将(X1xj,Y1xj)，(X2xj,Y2xj)，j＝1,2,…15转换为中心统一直角坐标(X1txj,Y1txj)和(X2txj,Y2txj)：X1txj＝X1xjcosδxc-Y1xjsinδxc+Xcx-Δx1Y1txj＝X1xjsinδxc+Y1xjcosδxc+Ycx-Δy1X21txj＝X2xjcosδxc-Y2xjsinδxc+Xcx-Δx2Y2txj＝X2xjsinδxc+Y2xjcosδxc+Ycx-Δy2计算结果见表10。表10：次站雷达系统误差修正数据的统一直角坐标主、次站雷达原始测量点和两组次站雷达系统误差修正点在统一直角坐标系中的显示如图2所示。步骤4：在统一直角坐标系中使用加权直线航迹线参数估计模型分别对主、次站雷达观测数据和次站雷达系统误差修正数据进行直线参数迭代估计，得到主站观测航迹线参数(kz,dz)、次站观测航迹线参数(kc,dc)和两组次站修正航迹线参数(k1x,d1x)、(k2x,d2x)。步骤4.1：采用不加权直线航迹线参数估计模型粗略估计主站雷达观测到的直线航迹线y-k1x-d1＝0(k1为直线的斜率，d1为直线在Y轴上的截距)，包括如下步骤：步骤4.1.1：用主站雷达所有观测点{(Xtzi,Ytzi),i＝1,2,...15}(简记为：{(xi,yi)，i＝1,2,…15})到该直线的距离的平方和最小作为条件构造直线，计算在此条件下的这条直线的最佳参数(k,d)。对(1)式的求解步骤如下。①计算a′1,a′2,b′1,b′2,c′0(n＝15)。②计算a′,b′,c′。a′＝c′0-a′1b′1＝115.0896，b′＝a′2-b′2-a1′2+b1′2＝41.8806，c′＝a′1b′1-c′0＝-115.0896。③解方程，计算所有解。d′1＝b′1-a′1k′1＝502.4071，d′2＝b′1-a′1k′2＝3584.1047。步骤4.1.2：按照距离最小原则确定方程的合理解。计算测量点(x1,y1)分别到直线y＝k′1×x+d′1和直线y＝k′2×x+d′2的距离l′1,l′2：因|l′1|<|l′2|，则取(k′1,d′1)＝(1.1984，502.4071)作为所求直线的合理参数，记为(k1,d1)。步骤4.2：采用加权直线航迹线参数估计模型迭代估计主站雷达观测到的直线航迹线y-kzx-dz＝0，包括如下步骤：步骤4.2.1：计算主站各观测点(xi,yi)到直线y-k1x-d1＝0的距离之和(n＝15)。步骤4.2.2：求各观测点(xi,yi)到直线y-k(s)x-d(s)＝0的距离li。计算结果如表11所示。式中s表示迭代次数，初始值为1，即：k(1)＝k1，d(1)＝d1。步骤4.2.3：求|li|的倒数。计算结果如表11所示。步骤4.2.4：求各点的权值vi。计算结果如表11所示。表11距离和权值计算结果序号liPivi1-0.05630817.7594870.12277520.2538683.9390550.0272323-0.04100924.3847930.1685774-0.2690433.7168760.02569650.2626663.8071170.02631960.2622703.8128700.0263597-0.1879805.3197260.0367768-0.5185771.9283530.01333190.2495624.0070150.02770110-0.1599226.2530350.04322911-0.2018664.9537790.034247120.3639622.7475400.01899413-0.1476626.7722200.046818140.2098484.7653470.03294415-0.01980950.4833500.349002步骤4.2.5：求解加权直线航迹线参数估计模型。用主站雷达所有观测点{(xi,yi)，i＝1,2,…15}到该直线的加权距离(vi×li)的平方和最小作为条件构造直线，计算在此条件下的这条直线的最佳参数(k,d)。对(3)式的求解步骤如下。①计算a0,a1,a2,b1,b2,c0(n＝15)。②计算a,b,c。a＝-c0-a1b1＝189.5367，c″＝c0+a1b1＝-189.5367。③解方程，计算所有解。d1＝b1-a1k1＝500.5643，d2＝b1-a1k2＝3599.1113。步骤4.2.6：按照距离最小原则确定方程的合理解。计算测量点(x1,y1)分别到直线y＝k1×x+d1和直线y＝k2×x+d2的距离l1,l2：因|l1|+|l2|＝37.192271<Lmin(初值为106)，则Lmin＝37.192271。s加1后值为2，且|l1|<|l2|，则取(k1,d1)＝(1.1995，500.5643)作为所求直线的合理参数，记为(k(2),d(2))。步骤4.2.7：计算所有观测点到新直线y-k(2)x-d(2)＝0的加权距离之和f(2)(k(2),d(2))。步骤4.2.8：判别是否为最佳解。因f(1)(k(1),d(1))＝3.204352，f(2)(k(2),d(2))＝0.083881,f(2)(k(2),d(2))<f(1)(k(1),d(1))，则重复步骤4.2.2至步骤4.2.7。经计算，s＝3、执行步骤4.2.6时，|l1|+|l2|＝739.048264，大于上一次迭代时的Lmin＝37.192271，因此输出解(k(2),d(2))，并简记为(kz,dz)＝(1.1995，500.5643)。步骤4.3：采用不加权直线航迹线参数估计模型粗略估计次站雷达观测到的直线航迹线y-kc1x-dc1＝0。具体过程与步骤4.1类同。步骤4.4：采用加权直线航迹线参数估计模型迭代估计次站雷达观测到的直线航迹线y-kcx-dc＝0。具体过程与步骤4.2类同。计算结果为：(kc,dc)＝(1.004519，791.838008)。步骤4.5：采用不加权直线航迹线参数估计模型粗略估计次站雷达数据采用系统误差估值1修正后的直线航迹线y-k1x1x-d1x1＝0。具体过程与步骤4.1类同。步骤4.6：采用加权直线航迹线参数估计模型迭代估计次站雷达数据采用系统误差估值1修正后的直线航迹线y-k1xx-d1x＝0。具体过程与步骤4.2类同。计算结果为：(k1x,d1x)＝(1.153464，570.324616)。步骤4.7：采用不加权直线航迹线参数估计模型粗略估计次站雷达数据采用系统误差估值2修正后的直线航迹线y-k2x1x-d2x1＝0。具体过程与步骤4.1类同。步骤4.8：采用加权直线航迹线参数估计模型迭代估计次站雷达数据采用系统误差估值2修正后的直线航迹线y-k2xx-d2x＝0。具体过程与步骤4.2类同。计算结果为：(k2x,d2x)＝(1.095545，655.384635)。步骤5：分别计算次站观测航迹线、两组次站修正航迹线与主站观测航迹线的相对熵JSc、JS1x和JS2x。相对熵越小，表明两条航迹线越相似。所述步骤5包括如下步骤：步骤5.1：计算次站观测航迹线与主站观测航迹线的相对熵JSc。①计算f(xi′)和g(xi′)。f(xi′)＝kzxi′+dz，g(xi′)＝kcxi′+dc，i′＝1,2,...30；xi′∈{(Xtzi|i＝1,2,...15)∪(Xtcj|j＝1,2,...15)}计算结果如表12所示。表12f(xi')与g(xi')计算结果②计算KL(f(x)||g(x))和KL(g(x)||f(x))；③计算JSc。步骤5.2：计算第一条次站修正航迹线与主站观测航迹线的相对熵JS1x。①计算f′(x′i′)和g′(x′i′)。f′(x′i′)＝kzx′i′+dz，g′(x′i′)＝kxx′i′+dx，i′＝1,2,...30。其中：x′i′∈{(Xtzi|i＝1,2,...15)∪(X1txj|j＝1,2,...15)}，表示主站观测点与次站第一组修正点横坐标的并集。计算结果如表13所示。表13f′(x′i′)与g′(x′i′)计算结果②计算KL′(f′(x)||g′(x))和KL′(g′(x)||f′(x))。③计算JS1x。步骤5.3：计算第二条次站修正航迹线与主站观测航迹线的相对熵JS2x。具体过程与步骤5.2类同。计算结果为：JS2x＝0.084058。由相对熵JSc>JS2x>JS1x的结果可以看出，JS1x＝0.002087对应的次站修正航迹线与主站观测航迹线更加相似(接近)，JS2x＝0.084058对应的次站修正航迹线次之；相比之下，JSc＝0.188534对应的次站观测航迹线与主站观测航迹线相似度最差，与图2的观察结果一致。步骤6：则次站两组修正航迹线系统误差估值有效性指标V1和V2计算结果分别为：由V1>V2>0可以得出以下结论：第一，V1和V2对应的系统误差估值均能有效减小主、次站雷达观测航迹的分裂程度。第二，在减小主、次站雷达观测航迹的分裂程度上，V1对应的系统误差估值(Δρ1,Δθ1,Δx1,Δy1)＝(1.4,4.0,0.4,0.1)要优于V2对应的系统误差估值(Δρ2,Δθ2,Δx2,Δy2)＝(0.5,2.5,-0.3,0.1)，取得了较好的航迹修正效果,与图2的观察结果一致。本发明通过选取主、次站雷达对某个典型航路目标的一组观测数据，经中心统一直角坐标转换后，先后采用不加权直线航迹线参数估计模型和加权直线航迹线参数估计模型对两雷达观测航迹线、次站雷达系统误差修正航迹线进行直线参数迭代估计，通过计算次站观测航迹线系统误差修正前后与主站观测航迹线的相对熵，构建次站修正航迹线系统误差估值有效性指标，用以评价系统误差估值的有效性以及修正效果。本发明原理方法科学、实施步骤合理，与传统系统误差有效性评价方法相比，在提高准确性和可操作性的同时，简化了复杂度，方便了工程上的实现。本发明有助于提高雷达阵地检飞、数据校准方法选择和雷达相对系统误差估值表计算过程中对系统误差估值评价的科学性、准确性和可操作性，从而提升雷达网目标状态估计和多雷达数据融合的质量。本发明所提供的方法的时间复杂度和空间复杂度都很低，可操作性和实用性很强。以上所述仅是本发明的优选实施方式，应当指出，对于本
技术领域：
的普通技术人员来说，在不脱离本发明技术原理的前提下，还可以做出若干改进和变形，这些改进和变形也应视为本发明的保护范围。例如但不限于以下几点：(1)步骤3.3中引用的次站系统误差估值(Δρ,Δθ,Δx,Δy)仅仅是系统误差的某一种表述方式。系统误差在不同表述模型中具有多种样式，如增加时间误差、高度误差、测距系数误差等，在采用GPS等精确测量仪器后可以忽略定位误差(Δx,Δy)等。无论采取何种系统误差模型，都是用以将次站测量数据转换为修正数据，最终转换成统一直角坐标，进而转入步骤4。所以，系统误差模型的多样性不影响本发明的有效性，而应看成是本发明的某种变形应用。(2)一定要确保步骤2和步骤3中将极坐标转换为中心统一直角坐标时，横、纵坐标值均大于零，必要时需要对中心统一直角坐标系原点进行平移。(3)步骤5中，当航迹线与X轴夹角大于75度时，对于f(xi)＝kzxi+dz、g(xj)＝kcxj+dc和g(xj)＝kxxj+dx，要分别替换为：和(4)本发明所述实施步骤2和步骤3是针对两坐标雷达的表述，而对于三坐标雷达测量数据本发明同样适用。只要将相应的极坐标到雷达站二维直角坐标转换公式在考虑目标高度测量值的情况下稍作调整即可。(5)本发明所述步骤1选取主、次站雷达极坐标值(ρzi,θzi,tzi)作为观测数据适用于普遍情况。如果在雷达检飞或雷达站具备ADS-B接收设备的情况下，可以获得来自GPS或ADS-B设备更加准确、数据率更高的目标三维位置信息(通常以地理坐标形式给出)。此时只要将步骤2和步骤3中的极坐标到统一直角坐标转换公式替换为地理坐标到统一直角坐标的转换公式即可。当前第1页1 2 3