本发明属于惯性导航系统精度估计技术领域,尤其涉及一种三参数椭圆概率误差描述导航落点精度的估计方法及系统。
背景技术:
惯性导航系统精度的衡量方法通常有三种,包括圆概率误差cep(circularerrorprobability)、不相关两维的落点精度估计方法和相关两维的落点精度估计方法。cep是以理论位置为圆心,围绕圆心包含50%的导弹弹头弹着点的圆的半径,是一种简单且近似的精度衡量方法。而不相关、相关二维的椭圆概率误差落点精度估计方法是基于cep的改进方法,主要解决在横、纵向标准方差在相关和不相关的情况下,对落点精度分布特性进行更加精准地估计。
不相关二维的椭圆概率误差落点精度估计方法具体方法如下所示:
1.确定椭圆在横向位置误差方向的半轴长度a满足方程
2.确定椭圆在纵向误差位置方向的半轴长度b满足方程
其中
求解后得到
其中θ和r分别为极坐标的参数。
相关二维的椭圆概率误差落点精度估计方法具体方法如下所示:
1.椭圆概率误差方程满足
2.由步骤(1)解得半径r为
3.将r表达式带入到
4.利用eep=(a,b)表示相关二维的落点精度
cep、不相关和相关两维的椭圆概率误差落点精度估计是目前评估惯性导航系统精度估计的常用的三种方法,但是存在以下问题:
(1)cep对于标准误差相等时候效果较好,在标准误差不等的时候,落点位置在两个方向分布的距离偏差较大,cep的估计方法不够合理;
(2)不相关二维的椭圆概率误差落点精度估计没有充分考虑相关系数的影响,在横向和纵向有相互影响的时候,不能够很好地解决落点精度估计在相关性方面的问题。
(3)相关两维的椭圆概率误差落点精度虽然考虑了横向和纵向之间的相关系数ρ,但是得到椭圆的面积不是最小的,因此仍不能精确表示落点分布的方向性和相关性特征。
通过以上分析可以看出,目前惯性导航系统精度估计研究中的cep、不相关和相关二维的落点精度估计方法不能给工程人员提供通用、直观以及更加准确的精度评估方法。因此需要对惯性导航系统精度估计继续深入研究,以提供更加完善合理的方法。
技术实现要素:
本发明解决的技术问题是:克服现有技术的不足,提供了一种三参数椭圆概率误差描述惯性导航落点精度的估计方法及系统,既充分考虑了落点分布在不同方向的差异性和方向的相关性,并且与cep和其他椭圆概率误差的精度评定方法相比,三参数的椭圆概率误差模型所得面积最小,同时给出了通用表达式,严格表达在横向和纵向相关的情况下,导航落点精度在解耦坐标轴下横向和纵向的精度估计。
本发明目的通过以下技术方案予以实现:根据本发明的一个方面,提供了一种三参数椭圆概率误差描述惯性导航落点精度的估计方法,所述方法包括如下步骤:
(1)利用惯性导航系统得到运载体n个落点位置,n个落点位置的集合设为落点域;
(2)根据落点域得到落点域中心的坐标o,统计落点位置与落点域中心的横向位置误差和纵向位置误差,并建立直角坐标系oxy,使横向位置误差的方向落在x轴,纵向位置误差的方向落在y轴,并计算得到横向位置误差的标准差σx、纵向位置误差的标准差σy和两者的相关系数ρxy;
(3)通过横向位置误差的标准差σx、纵向位置误差的标准差σy和两者的相关系数ρxy计算得出旋转角度θ;以原点o为旋转中心,旋转θ角得到新坐标系ox'y';计算得到新坐标系下的横向位置标准方差σx′和纵向位置标准方差σy′;
(4)根据步骤(3)得到的横向位置标准方差σx′和纵向位置标准方差σy′,以原点o为中心,以
(5)根据步骤(4)绘制的椭圆得到该惯性导航系统相关两维的三参数描述的落点精度
上述三参数椭圆概率误差描述惯性导航落点精度的估计方法中,在步骤(3)中,旋转角度θ的公式如下:
上述三参数椭圆概率误差描述惯性导航落点精度的估计方法中,在步骤(3)中,新坐标系下的横向位置标准方差σx′满足如下公式:
上述三参数椭圆概率误差描述惯性导航落点精度的估计方法中,在步骤(3)中,新坐标系下的纵向位置标准方差σy′满足如下公式:
上述三参数椭圆概率误差描述惯性导航落点精度的估计方法中,在步骤(3)中,直角坐标系oxy和新坐标系ox'y'满足如下关系:
其中,(x',y')为ox'y'坐标系下的坐标;(x,y)为oxy坐标系下的坐标;θ旋转角度。
上述三参数椭圆概率误差描述惯性导航落点精度的估计方法中,在步骤(4)中,椭圆方程在坐标系ox′y′下表示为:
根据本发明的一个方面,还提供了一种三参数椭圆概率误差描述惯性导航落点精度的估计系统,包括:第一模块,用于利用惯性导航系统得到运载体若干个落点位置,若干个落点位置的集合设为落点域;第二模块,用于根据落点域得到落点域中心的坐标o,统计落点位置与落点域中心的横向位置误差和纵向位置误差,并建立直角坐标系oxy,使横向位置误差的方向落在x轴,纵向位置误差的方向落在y轴,并计算得到横向位置误差的标准差σx、纵向位置误差的标准差σy和两者的相关系数ρxy;第三模块,用于通过横向位置误差的标准差σx、纵向位置误差的标准差σy和两者的相关系数ρxy计算得出旋转角度θ;以原点o为旋转中心,旋转θ角得到新坐标系ox'y';计算得到新坐标系下的横向位置标准方差σx′和纵向位置标准方差σy′;第四模块,用于根据横向位置标准方差σx′和纵向位置标准方差σy′,以原点o为中心,以
上述三参数椭圆概率误差描述惯性导航落点精度的估计系统中,旋转角度θ的公式如下:
上述三参数椭圆概率误差描述惯性导航落点精度的估计系统中,新坐标系下的纵向位置标准方差σy′满足如下公式:
上述三参数椭圆概率误差描述惯性导航落点精度的估计系统中,直角坐标系oxy和新坐标系ox'y'满足如下关系:
其中,(x',y')为ox'y'坐标系下的坐标;(x,y)为oxy坐标系下的坐标;θ旋转角度。
本发明与现有技术相比具有如下有益效果:
(1)本发明充分考虑横向和纵向相关的情况下,利用横向位置误差的标准差、纵向位置误差的标准差和相关系数,对惯性导航系统落点精度的进行估计,更具有一般性和通用性。
(2)本发明给出了三参数椭圆的长轴、短轴和旋转角度描述的相关两维的落点精度的表达式(
(3)本发明的精度估计推理严格,不同条件考虑全面,同时给出明确表达式,简单、易于计算,适合工程特性。
附图说明
通过阅读下文优选实施方式的详细描述,各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的,而并不认为是对本发明的限制。而且在整个附图中,用相同的参考符号表示相同的部件。在附图中:
图1是本发明实施例提供的三参数椭圆概率误差描述惯性导航落点精度的估计方法的流程示意图;
图2是本发明实施例提供的相关条件下落点分布示意图;
图3是本发明实施例提供的不相关二维椭圆概率误差模型的示意图;
图4是本发明实施例提供的ox’y’坐标系下的椭圆概率误差模型的示意图;
图5是本发明实施例提供的oxy坐标下的三参数椭圆概率误差模型的示意图;
图6是本发明实施例提供的相关条件下两参数和三参数椭圆概率误差模型面积对比图。
具体实施方式
下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例性实施例,然而应当理解,可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了能够更透彻地理解本公开,并且能够将本公开的范围完整的传达给本领域的技术人员。需要说明的是,在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
图1是本发明实施例提供的三参数椭圆概率误差描述惯性导航落点精度的估计方法的流程示意图。图1所示,该方法包括如下步骤:
(1)利用惯性导航系统得到运载体n个落点位置,称为落点域;
(2)根据落点域得到落点域中心的坐标o。统计落点位置与落点域中心的横向位置误差和纵向位置误差,并建立直角坐标系oxy,使横向位置误差的方向落在x轴,纵向位置误差的方向落在y轴。同时,计算得到横向位置误差的标准差σx、纵向位置误差的标准差σy和两者的相关系数ρ,落点分布如图2所示;
(3)通过横向位置误差的标准差σx、纵向位置误差的标准差σy和两者的相关系数ρ计算得出旋转角度θ。以o为旋转中心,旋转θ角得到新坐标系ox'y'。计算得到新坐标系下的横向位置标准方差σx′和纵向位置标准方差σy′,两者由相关转变为不相关;
(4)根据步骤(3)得到的横向位置标准方差σx′和纵向位置标准方差σy′,以原点o为中心,以
(5)根据步骤(4)绘制的椭圆得到该惯性导航系统相关两维的三参数描述的落点精度
所述步骤(4)中建立的椭圆满足如下形式:
首先,椭圆相对于相关二维的原始坐标系的旋转矩阵满足如下方程:
其中,x’,y’为ox’y’坐标系下的椭圆坐标;x,y为oxy坐标系下的椭圆坐标;θ旋转角度。
其中,旋转角度进一步表示为:
在坐标系ox′y′下椭圆方程表示为
其中,
椭圆长、短半轴分别为
本实施例采用三参数椭圆概率误差描述惯性导航落点精度的估计方法,进行仿真验证,仿真结果为图2、图3、图4、图5、图6所示。
图2为相关条件下落点分布示意图。由图2知,导弹落点在横向、纵向的分布不均匀,纵向的落点分布比横向落点分布的距离更长,约为两倍关系,相反,密集度就更低。另外,落点在oxy直角坐标系中的分布呈现一定的角度,说明落点在横向和纵向上存在一定的相关性。
图3为不相关二维椭圆概率误差模型。由图3知,模拟打靶1000发,落点分布沿x轴和y轴方向标准差分别为σx=100m、σy=400m,相关系数ρxy=0.5。根据标准差相关时的两参数椭圆概率误差可得,eep=(113.6581,454.6326)。该椭圆包含了50%的落点概率,但是此椭圆的面积并非最小。
图4为ox’y’坐标系下的椭圆概率误差模型。由图4知,经过坐标旋转后,在坐标系ox′y′中落点x轴和y轴方向标准差分别为σx’=403.2628和σy’=85.9018,在旋转后的ox′y′坐标系中,椭圆概率误差eep=(403.2628,85.9018)。由ox′y′坐标系中的椭圆可知,旋转后的椭圆,解决了落点分布相关性的问题,准确的表示了落点误差精度。
图5为oxy坐标下的三参数椭圆概率误差模型。由图5知,沿椭圆的长轴和短轴,概率为50%的落点的位置分别在短轴=85.9018m和长轴403.2628m的椭圆以内;椭圆的旋转方向为82.5343°,即在82.5343°方向上落点分布密集度最低,误差相对较大。经过旋转的椭圆概率误差模型,不仅可以描述落点分布沿两个轴分布的不均匀性,也可以描述落点分布的相关性。
图6为相关条件下两参数和三参数椭圆概率误差模型面积对比图。由图6知,上面的曲线为两参数椭圆概率误差计算的面积s1,下面的曲线为三参数椭圆概率误差计算的面积s2。可以看出,无论ρxy取何值,始终有s1≥s2,也就是说,两个椭圆都包含了50%的落点,三参数椭圆概率误差为最优。
本实施例还提供了一种三参数椭圆概率误差描述惯性导航落点精度的估计系统,包括:第一模块,用于利用惯性导航系统得到运载体若干个落点位置,若干个落点位置的集合设为落点域;第二模块,用于根据落点域得到落点域中心的坐标o,统计落点位置与落点域中心的横向位置误差和纵向位置误差,并建立直角坐标系oxy,使横向位置误差的方向落在x轴,纵向位置误差的方向落在y轴,并计算得到横向位置误差的标准差σx、纵向位置误差的标准差σy和两者的相关系数ρxy;第三模块,用于通过横向位置误差的标准差σx、纵向位置误差的标准差σy和两者的相关系数ρxy计算得出旋转角度θ;以原点o为旋转中心,旋转θ角得到新坐标系ox'y';计算得到新坐标系下的横向位置标准方差σx′和纵向位置标准方差σy′;第四模块,用于根据横向位置标准方差σx′和纵向位置标准方差σy′,以原点o为中心,以
本实施例充分考虑横向和纵向相关的情况下,利用横向位置误差的标准差、纵向位置误差的标准差和相关系数,对惯性导航系统落点精度的进行估计,更具有一般性和通用性;本实施例给出了三参数椭圆的长轴、短轴和旋转角度描述的相关两维的落点精度的表达式(
以上所述的实施例只是本发明较优选的具体实施方式,本领域的技术人员在本发明技术方案范围内进行的通常变化和替换都应包含在本发明的保护范围内。