CN104197953A - 一种惯性导航***两维落点精度估计方法 - Google Patents
一种惯性导航***两维落点精度估计方法 Download PDFInfo
- Publication number
- CN104197953A CN104197953A CN201410397787.XA CN201410397787A CN104197953A CN 104197953 A CN104197953 A CN 104197953A CN 201410397787 A CN201410397787 A CN 201410397787A CN 104197953 A CN104197953 A CN 104197953A
- Authority
- CN
- China
- Prior art keywords
- sigma
- inertial navigation
- navigation system
- along track
- ellipse
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C25/00—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass
- G01C25/005—Manufacturing, calibrating, cleaning, or repairing instruments or devices referred to in the other groups of this subclass initial alignment, calibration or starting-up of inertial devices
Landscapes
- Engineering & Computer Science (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Navigation (AREA)
Abstract
本发明公开了一种惯性导航***两维落点精度估计方法,通过采集惯性导航***n次导航后运载体停止的位置信息得到惯性导航纵向偏差和横向偏差的标准偏差,建立平面直角坐标系OXY,利用标准偏差建立一个椭圆,在该椭圆范围内包含了惯性导航在水平面内50%的落点。本发明提供了一种对水平面内的惯性导航落点精度的描述方法,相对传统的圆概率误差(CEP)的描述方法,具有可得到精确表达式、准确描述落点特征的优点。
Description
技术领域
本发明涉及一种两维落点精度估计方法,属于惯性导航***精度估计领域。
背景技术
惯性导航***的精度可用与标准误差的一致程度来表示,最常用的惯性导航***的精度度量方法是基于二维平面的圆概率误差CEP(Circular ErrorProbablity)。这种精度度量法为弹道导弹弹着点位置精度提供了一种简单的度量,它定义为围绕实际目标包含50%的导弹弹头弹着点的圆的半径,CEP圆的圆心是实际水平面的位置。
CEP的理论计算方法如下:
在水平面内定义坐标系OXY,其中O点为中心点,设在该坐标系内点的坐标(x,y)服从二维正态分布且相互独立,记
其均值为
其标准误差为
该点处于以原点为圆心、半径为R的圆内的概率由误差分量高斯分布的二重积分给出。对于高斯分布、零均值的x和y,其概率为
上式是变量R的函数,即,该点处于圆内的概率的大小,与半径R的大小有关。如果用极坐标,式(1)的二重积分就能够以适当的形式计算出来。令极坐标的参数为θ和r。则
利用二重积分极坐标变换公式
式(1)变成
利用上式,可以在已知概率P时近似计算其对应的R值,但求解过程十分繁琐。
设σx=100m、σy=400m,则P(x2+y2<R2)随着R的变化过程如图1所示。当CEP=R=290.2m时,有P(x2+y2<R2)=0.5。此时,可保证50%的弹着点在此圆内。
设σx=100m、σy=400m,则模拟打靶1000次的落点分布如图2所示。
从图2可以看出,弹着点在X方向基本都小于此圆的半径,但在Y方向有大量离散点都在此圆的外部。另外,只根据CEP值的大小不能确定X或Y轴哪个方向的密集度更大。因此,CEP适合于等标准误差时的精度分析,缺少方向性。同时,在求式(1)和式(3)的概率积分时,没有精确的表达式,导致所得CEP值为一个近似值,而非精确结果。
对于导航落点,实际落点在射击平面的投影点与目标点之间的距离称为纵向偏差;落点与射击平面间的距离称为横向偏差。
发明内容
本发明的技术解决问题是:克服现有技术的不足,提供一种惯性导航***两维落点精度估计方法,既给出惯性导航两维落点精度估计的相对精确地表达式,同时能判断出弹着点在纵向偏差方向和横向偏差方向的分布特性。
本发明的技术解决方案:一种惯性导航***两维落点精度估计方法,步骤如下:
(1)采集惯性导航***n次导航后运载体停止的n个位置信息,n个位置构成位置组,n至少为6;
(2)根据步骤(1)得到的n个位置信息得到位置组中心的坐标,然后测量每个位置相对于位置组中心的纵向偏差和横向偏差,在此基础上计算得到所有纵向偏差的标准差σ1和横向偏差的标准差σ2;
(3)以位置组中心为原点O,建立两维直角坐标系OXY,使纵向偏差的方向落在X轴或Y轴的其中一个轴上;
(4)根据步骤(2)得到的纵向偏差标准差σ1和横向偏差标准差σ2,在步骤(3)的坐标系内以原点O为中心建立一个椭圆,使该椭圆覆盖n/2个运载体停止位置;
(5)根据步骤(4)建立的椭圆得到该惯性导航***两维落点精度ECEP=(a,b),其中a和b分别为步骤(4)建立的椭圆在纵向偏差方向和横向偏差方向的半轴长度,从而完成落点精度估计。
所述步骤(4)中建立的椭圆满足如下形式:
椭圆在纵向偏差方向的半轴长度a满足方程
椭圆在横向偏差方向的半轴长度b满足方程
其中表示运载体停止的n个位置落在纵向偏差方向上[-a,+a]内的概率;表示运载体停止的n个位置落在横向偏差方向上[-b,+b]内的概率;
求解后得到
其中θ和r分别为极坐标的参数。
本发明与现有方法相比的优点在于:
(1)因为在本发明中分别计算两个方向的标准误差σ1、σ2,并以其计算椭圆长轴和短轴值,所以相对于CEP描述的惯性导航***的落点精度更精确;
(2)因为采用了两个参量来描述落点的分布特性,即使用作为指标,可更准确了解惯性导航结果在空间上的分布特性;
(3)本方法计算简便、快捷,不需要进行繁琐计算即可获得惯性导航***性能。
附图说明
图1为圆概率积分示意图;
图2为模拟打靶1000次的落点分布及CEP示意图;
图3为本发明方法流程图;
图4为模拟打靶1000次的落点分布及椭圆示意图。
具体实施方式
本发明涉及一种惯性导航***两维落点精度估计方法,其流程如图3所示,具体步骤如下:
(1)使用GPS采集惯性导航***n次导航后运载体停止的n个位置信息,n个位置信息构成位置组信息,n至少为6;
(2)根据步骤(1)得到的位置信息,得到位置组中心的坐标,然后测量每个位置相对于位置组中心的纵向偏差和横向偏差,在此基础上通过统计计算得到n个纵向偏差的标准差σ1和n个横向偏差的标准差σ2。
例如,经过10次试验,得到使用GPS采集得到的用经纬度表示的10个最终位置点,通过分别对经度和纬度取平均,可以获得位置组中心的经纬度坐标,然后求取测量点与位置组中心的纵向误差和横向误差,得到n个纵向偏差和n个横向偏差,基于这两组值可以得到对应的纵向偏差标准差和横向偏差标准差。
(3)以位置组中心为原点O,建立两维直角坐标系OXY,使纵向偏差的方向落在X轴或Y轴的其中一个轴上;
(4)根据步骤(2)得到的n个运载体停止位置的纵向偏差标准差σ1和横向偏差标准差σ2,在步骤(3)的坐标系内以原点O为中心建立一个椭圆,使该椭圆覆盖n/2个运载体停止位置。
其中,椭圆在纵向偏差方向的半轴长度a满足方程
椭圆在横向偏差方向的半轴长度b满足方程
上式中,表示n个运载体停止位置中,在纵向偏差方向上落在[-a,+a]内的概率;表示n个运载体停止位置中,在横向偏差方向上落在[-b,+b]内的概率;
求解后得到
θ和r分别为极坐标的参数。
(5)根据步骤(4)建立的椭圆得到该惯性导航***两维落点精度ECEP=(a,b),其中a和b分别为步骤(4)建立的椭圆在纵向偏差方向和横向偏差方向的半轴长度。
本发明方法可以应用在导弹的落点精度估计中,例如,当根据本方法在某次导弹的落点精度估计中采集1000次导弹的落点位置,然后经过计算得到σ1=100m、σ2=400m,根据上述方法建立坐标系OXY,使1000次导弹的落点位置的位置组中心为原点O,纵向偏差的方向落在Y轴上,在该坐标系内建立椭圆,使该椭圆覆盖500个导弹的落点位置,所建立的椭圆的长半轴为470.964m,短半轴为117.741m,即ECEP=(470.964m,117.741m),即,470.964m为椭圆在纵向偏差方向上的半轴长度,117.741m为椭圆在横向偏差方向上的半轴长度。模拟打靶1000次的落点分布如图4所示,图4中同时加入了利用CEP进行评估的结果和利用本发明方法(椭圆)进行评估的结果。从图4中可以看出,用椭圆来描述落点分布相对精确,而CEP评估结果并不能准确描述落点特征。
本发明未详细描述内容为本领域技术人员公知技术。
Claims (2)
1.一种惯性导航***两维落点精度估计方法,其特征在于步骤如下:
(1)采集惯性导航***n次导航后运载体停止的n个位置信息,n个位置构成位置组,n至少为6;
(2)根据步骤(1)得到的n个位置信息得到位置组中心的坐标,然后测量每个位置相对于位置组中心的纵向偏差和横向偏差,在此基础上计算得到所有纵向偏差的标准差σ1和横向偏差的标准差σ2;
(3)以位置组中心为原点O,建立两维直角坐标系OXY,使纵向偏差的方向落在X轴或Y轴的其中一个轴上;
(4)根据步骤(2)得到的纵向偏差标准差σ1和横向偏差标准差σ2,在步骤(3)的坐标系内以原点O为中心建立一个椭圆,使该椭圆覆盖n/2个运载体停止位置;
(5)根据步骤(4)建立的椭圆得到该惯性导航***两维落点精度ECEP=(a,b),其中a和b分别为步骤(4)建立的椭圆在纵向偏差方向和横向偏差方向的半轴长度,从而完成落点精度估计。
2.根据权利要求1所述的一种惯性导航***两维落点精度估计方法,其特征在于:所述步骤(4)中建立的椭圆满足如下形式:
椭圆在纵向偏差方向的半轴长度a满足方程
椭圆在横向偏差方向的半轴长度b满足方程
其中表示运载体停止的n个位置落在纵向偏差方向上[-a,+a]内的概率;表示运载体停止的n个位置落在横向偏差方向上[-b,+b]内的概率;
求解后得到
其中θ和r分别为极坐标的参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410397787.XA CN104197953B (zh) | 2014-08-13 | 2014-08-13 | 一种惯性导航***两维落点精度估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410397787.XA CN104197953B (zh) | 2014-08-13 | 2014-08-13 | 一种惯性导航***两维落点精度估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104197953A true CN104197953A (zh) | 2014-12-10 |
CN104197953B CN104197953B (zh) | 2016-08-24 |
Family
ID=52083281
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410397787.XA Active CN104197953B (zh) | 2014-08-13 | 2014-08-13 | 一种惯性导航***两维落点精度估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104197953B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105629253A (zh) * | 2015-12-25 | 2016-06-01 | 北京航天控制仪器研究所 | 一种靶场末区多方式融合落点测量*** |
CN106225790A (zh) * | 2016-07-13 | 2016-12-14 | 百度在线网络技术(北京)有限公司 | 一种无人车定位精度的确定方法及装置 |
CN107270943A (zh) * | 2017-07-28 | 2017-10-20 | 北京航天控制仪器研究所 | 一种相关两维的惯性导航落点精度评估方法 |
CN108458727A (zh) * | 2018-02-11 | 2018-08-28 | 北京航天控制仪器研究所 | 弹道导弹惯性测量***精度指标自适应分配方法及*** |
CN108519104A (zh) * | 2018-02-11 | 2018-09-11 | 北京航天控制仪器研究所 | 三参数椭圆概率误差描述导航落点精度的估计方法及*** |
CN108548537A (zh) * | 2018-02-11 | 2018-09-18 | 北京航天控制仪器研究所 | 六参数椭球概率误差描述导航落点精度的估计方法及*** |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1796230A (zh) * | 2004-12-21 | 2006-07-05 | 西昌卫星发射中心 | 液体火箭***碎片散布范围的确定方法 |
WO2007016905A1 (de) * | 2005-08-09 | 2007-02-15 | Eads Deutschland Gmbh | Verfahren zur flugführung mehrerer im verband fliegender flugzeuge |
CN101614507A (zh) * | 2008-06-29 | 2009-12-30 | 范海涛 | 一种导弹反拦截方法 |
CN101620672A (zh) * | 2009-08-14 | 2010-01-06 | 华中科技大学 | 一种利用立体地标定位识别地面立体建筑物的方法 |
CN101629827A (zh) * | 2009-08-14 | 2010-01-20 | 华中科技大学 | 一种飞行器前视末制导导航定位方法 |
CN102865880A (zh) * | 2011-07-04 | 2013-01-09 | 中国科学院沈阳自动化研究所 | 一种预计水下机器人导航***圆概率偏差的方法 |
CN103955227A (zh) * | 2014-04-29 | 2014-07-30 | 上海理工大学 | 无人机精准降落的控制方法 |
-
2014
- 2014-08-13 CN CN201410397787.XA patent/CN104197953B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1796230A (zh) * | 2004-12-21 | 2006-07-05 | 西昌卫星发射中心 | 液体火箭***碎片散布范围的确定方法 |
WO2007016905A1 (de) * | 2005-08-09 | 2007-02-15 | Eads Deutschland Gmbh | Verfahren zur flugführung mehrerer im verband fliegender flugzeuge |
CN101614507A (zh) * | 2008-06-29 | 2009-12-30 | 范海涛 | 一种导弹反拦截方法 |
CN101620672A (zh) * | 2009-08-14 | 2010-01-06 | 华中科技大学 | 一种利用立体地标定位识别地面立体建筑物的方法 |
CN101629827A (zh) * | 2009-08-14 | 2010-01-20 | 华中科技大学 | 一种飞行器前视末制导导航定位方法 |
CN102865880A (zh) * | 2011-07-04 | 2013-01-09 | 中国科学院沈阳自动化研究所 | 一种预计水下机器人导航***圆概率偏差的方法 |
CN103955227A (zh) * | 2014-04-29 | 2014-07-30 | 上海理工大学 | 无人机精准降落的控制方法 |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105629253A (zh) * | 2015-12-25 | 2016-06-01 | 北京航天控制仪器研究所 | 一种靶场末区多方式融合落点测量*** |
CN106225790A (zh) * | 2016-07-13 | 2016-12-14 | 百度在线网络技术(北京)有限公司 | 一种无人车定位精度的确定方法及装置 |
CN106225790B (zh) * | 2016-07-13 | 2018-11-02 | 百度在线网络技术(北京)有限公司 | 一种无人车定位精度的确定方法及装置 |
CN107270943A (zh) * | 2017-07-28 | 2017-10-20 | 北京航天控制仪器研究所 | 一种相关两维的惯性导航落点精度评估方法 |
CN107270943B (zh) * | 2017-07-28 | 2020-02-14 | 北京航天控制仪器研究所 | 一种相关两维的惯性导航落点精度评估方法 |
CN108458727A (zh) * | 2018-02-11 | 2018-08-28 | 北京航天控制仪器研究所 | 弹道导弹惯性测量***精度指标自适应分配方法及*** |
CN108519104A (zh) * | 2018-02-11 | 2018-09-11 | 北京航天控制仪器研究所 | 三参数椭圆概率误差描述导航落点精度的估计方法及*** |
CN108548537A (zh) * | 2018-02-11 | 2018-09-18 | 北京航天控制仪器研究所 | 六参数椭球概率误差描述导航落点精度的估计方法及*** |
CN108548537B (zh) * | 2018-02-11 | 2020-09-18 | 北京航天控制仪器研究所 | 六参数椭球概率误差描述导航落点精度的估计方法及*** |
CN108519104B (zh) * | 2018-02-11 | 2020-12-18 | 北京航天控制仪器研究所 | 三参数椭圆概率误差描述导航落点精度的估计方法及*** |
Also Published As
Publication number | Publication date |
---|---|
CN104197953B (zh) | 2016-08-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104197953B (zh) | 一种惯性导航***两维落点精度估计方法 | |
CN104197954B (zh) | 一种惯性导航***三维空间落点精度估计方法 | |
CN108061889A (zh) | Ais与雷达角度***偏差的关联方法 | |
CN104390646A (zh) | 水下潜器地形辅助惯性导航***的位置匹配方法 | |
CN104316079B (zh) | 一种基于火箭橇试验的惯性测量***落点精度估计方法 | |
CN108089148A (zh) | 一种基于时差信息的无源多站多目标测向交叉定位方法 | |
CN102426018A (zh) | 基于混合地形轮廓匹配tercom算法和粒子滤波的地形辅助导航方法 | |
CN106446422A (zh) | 一种基于对数似然估计的无源定位跟踪新方法 | |
CN104504728B (zh) | 多机动目标跟踪方法、***及其广义联合概率数据关联器 | |
CN105242252B (zh) | 基于图像匹配的下降轨聚束sar雷达定位方法 | |
CN103558602B (zh) | 一种用于多基地声纳配置方式的模拟退火定位方法 | |
CN113342059B (zh) | 基于位置和速度误差的多无人机跟踪移动辐射源方法 | |
CN107270943B (zh) | 一种相关两维的惯性导航落点精度评估方法 | |
CN103869279A (zh) | 一种多传感器平台的多目标定位跟踪方法 | |
CN103983969B (zh) | 基于二次拟合距离方程的地面加速运动目标成像方法 | |
CN106597428B (zh) | 一种海面目标航向航速估算方法 | |
CN104913743A (zh) | 基于惯性测量的电力铁塔变形监测方法 | |
CN108594193A (zh) | 一种基于固定目标与非合作目标的雷达***偏差估计方法 | |
CN103616024B (zh) | 一种行星探测进入段自主导航***可观测度确定方法 | |
CN106353756A (zh) | 基于图像匹配的下降轨聚束合成孔径雷达定位方法 | |
Gemerek et al. | Fugitive gas emission rate estimation using multiple heterogeneous mobile sensors | |
CN108917698B (zh) | 一种方位角计算方法 | |
CN104237862A (zh) | 基于ads-b的概率假设密度滤波雷达***误差融合估计方法 | |
CN104777465A (zh) | 基于b样条函数任意扩展目标形状及状态估计方法 | |
CN103491628A (zh) | 一种基于tpoa的nlos传输环境无线定位方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |