CN104197954B - Method for estimating precision of drop points of inertial navigation system in three-dimensional space - Google Patents
Method for estimating precision of drop points of inertial navigation system in three-dimensional space Download PDFInfo
- Publication number
- CN104197954B CN104197954B CN201410397823.2A CN201410397823A CN104197954B CN 104197954 B CN104197954 B CN 104197954B CN 201410397823 A CN201410397823 A CN 201410397823A CN 104197954 B CN104197954 B CN 104197954B
- Authority
- CN
- China
- Prior art keywords
- sigma
- integral
- ellipsoid
- inertial navigation
- navigation system
- 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.)
- Active
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
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
The invention discloses a method for estimating the precision of drop points of an inertial navigation system in the three-dimensional space. According to the method, after standard deviations of a longitudinal deviation, a lateral deviation and a height deviation of the inertial navigation system are known, a space rectangular coordinate system OXYZ is established, an ellipsoid can be obtained by virtue of the standard deviations, and 50% of drop points of the inertial navigation system in the three-dimensional space are within the range of the ellipsoid; compared with a traditional sphere error probability (SEP) description method, the method has the advantages that an expression formula is accurate, and the characteristics of the drop points are accurately described.
Description
Technical field
The present invention relates to a kind of impact accuracy method of estimation, belong to inertial navigation system Accuracy extimate field.
Background technology
In inertial navigation system, the important indicator evaluating its three-dimensional impact accuracy is based on three-dimensional ball probability by mistake
Difference sep (sphere error probablity).This Accuracy Measure method is to provide one for ballistic missile point of impact positional precision
Plant simple tolerance, sep is defined as comprising the radius of 50% ball of guided missile warhead point of impact around realistic objective.
The theoretical calculation method of sep is as follows:
Define coordinate system oxyz, point centered on wherein o point in three dimensions, be located at point in this coordinate system coordinate (x,
Y, z) obey three-dimensional normal distribution and separate, note
Its average is
Its standard error is
This point is in initial point as the center of circle, the probability in the ball as r for the radius by error component Gauss distribution triple integral
Be given.For x, y and z of Gauss distribution, zero-mean, its probability is
Above formula is the function of variable r, i.e. this point is in the size of the probability in ball, relevant with the size of radius r.If
With polar coordinate, triple integral is just calculated in a suitable form.Polar parameter is made to be θ and r.Then
Using triple integral polar coordinate transform formula
Probability integral becomes
Using above formula, can in known probability p its corresponding r value of approximate calculation, but solution procedure is very loaded down with trivial details.
If σx=10m, σy=30m, σz=40m, then p (x2+y2+z2<r2) as shown in Figure 1 with the Changing Pattern of r.From figure
In can obtain, in r=42.32, p (x2+y2+z2<r2)=0.5.
By above standard deviation simulation practice shooting 1000 times distribution as shown in Figure 2, in figure ball be sep ball.
From accompanying drawing 2 as can be seen that point of impact is substantially both less than the radius of this ball in a direction, but have in another two direction
A large amount of discrete points are all in the outside of this ball.In addition, only the size according to sep value not can determine that the closeness in which direction is bigger.
Therefore, sep such as is suitable at precision analysis during standard error, lacks directivity.Meanwhile, in the product of probability asking formula (1) and formula (3)
Timesharing, does not have accurate expression formula, leads to gained sep value for an approximation, rather than precise results.
For navigation drop point, actual falling point the plane of riee the distance between subpoint and impact point referred to as longitudinally by mistake
Difference;Distance between drop point and the plane of riee is referred to as lateral error;Deviation perpendicular to horizontal plane is referred to as height tolerance.
Content of the invention
The technology solve problem of the present invention is: overcomes the deficiencies in the prior art, provides a kind of inertial navigation system three-dimensional space
Between impact accuracy method of estimation, be both given inertial navigation three-dimensional impact accuracy estimate relatively precisely expression formula, can sentence simultaneously
Break the distribution character point of impact in longitudinal bias direction, lateral deviation direction and height tolerance direction.
The technical solution of the present invention: a kind of inertial navigation system three dimensions impact accuracy method of estimation, step is such as
Under:
(1) the n positional information that after the n navigation of collection inertial navigation system, carrier stops, n position composition position
Group, n is at least 6;
(2) the n positional information being obtained according to step (1), obtains the coordinate at set of locations center, then measures each position
With respect to the longitudinal bias at set of locations center, lateral deviation and height tolerance, it is calculated all longitudinal bias on this basis
Standard deviation sigma1, the standard deviation sigma of lateral deviation2Standard deviation sigma with height tolerance3;
(3) with set of locations center for initial point o, set up rectangular coordinate system in space oxyz, so that the direction of longitudinal bias is fallen in x-axis
Or on one of axle of y-axis, the direction of height tolerance falls in z-axis;
(4) standard deviation sigma of the longitudinal bias being obtained according to step (2)1, the standard deviation sigma of lateral deviation2With height tolerance
Standard deviation sigma3, setting up an ellipsoid centered on initial point o in the coordinate system of step (3), make this ellipsoid cover n/2 carrier
Stop position;
(5) this inertial navigation system three-dimensional impact accuracy esep=(a, b, c) is obtained according to the ellipsoid that step (4) is set up,
Wherein a, b, c are respectively the ellipsoid set up of step (4) the half of longitudinal bias direction, lateral deviation direction and height tolerance direction
Shaft length, thus complete the estimation of inertial navigation system three-dimensional impact accuracy.
The ellipsoid set up in described step (4) meets following form:
Ellipsoid meets equation in the half shaft length a in longitudinal bias direction
Ellipsoid meets equation in the half shaft length b in lateral deviation direction
Ellipsoid meets equation in the half shaft length c in height tolerance direction
WhereinRepresent that the n position that carrier stops falls vertical
Probability in in bias direction [- a ,+a];Represent that carrier stops
The probability that falls on lateral deviation direction [- b ,+b] of n position;
Represent the probability that the n position that carrier stops falls on height tolerance direction [- c ,+c];
Obtain after solution
A=1.539 σ1
B=1.539 σ2
C=1.539 σ3
Wherein θ, β and r are respectively polar parameter.
Present invention advantage compared with the conventional method is:
(1) because calculating standard error σ in three directions in the present invention respectively1、σ2、σ3, and calculate the three of ellipsoid with it
Individual axial length, can directly obtain the three-dimensional impact accuracy of inertial navigation system, and result has precise forms, non-approximated result;
(2) because employing three parameters to describe the distribution character of drop point, that is, using esep=(1.539 σ1,1.539
σ2,1.539σ3) as index, can more accurately understand inertial navigation result distribution character spatially, it is navigation error analysis
More accurately foundation is provided.
(3) this method calculates simplicity, fast it is not necessary to carry out troublesome calculation can obtain inertial navigation system performance.
Brief description
Fig. 1 is ball probability integral schematic diagram;
Fig. 2 is that simulation is practiced shooting the distribution of 1000 times and sep schematic diagram;
Fig. 3 is the inventive method flow chart;
Fig. 4 is that simulation is practiced shooting the distribution of 1000 times and ellipsoid schematic diagram.
Specific embodiment
A kind of inertial navigation system three dimensions impact accuracy method of estimation proposed by the present invention, flow chart as shown in figure 3,
Specifically comprise the following steps that
(1) the n positional information that after using gps to gather n navigation of inertial navigation system, carrier stops, n position letter
Breath constitutes position group information, and n is at least 6;
(2) the n positional information being obtained according to step (1), obtains the coordinate at set of locations center, then measures each position
With respect to the longitudinal bias at set of locations center, lateral deviation and height tolerance, on this basis by statistical computation obtain n vertical
Standard deviation sigma to deviation1, the standard deviation sigma of n lateral deviation2Standard deviation sigma with n height tolerance3.
For example, test through 10 times, collected with longitude, latitude, 10 final positions highly representing using gps
Point, by being averaged to three position coordinateses respectively, it is possible to obtain the longitude at set of locations center, latitude, height coordinate, Ran Houtong
Cross longitudinal error, lateral error and the height error asking for measurement point and set of locations center, obtain n longitudinal error, n horizontal
Error and n height error, can obtain corresponding longitudinal bias standard deviation, lateral deviation standard deviation and height based on this three class value
Degree deviation standard is poor.
(3) with set of locations center for initial point o, set up rectangular coordinate system in space oxyz, so that the direction of longitudinal bias is fallen in x-axis
Or on one of axle of y-axis, the direction of height tolerance falls in z-axis;
(4) standard deviation sigma being obtained according to step (2)1、σ2And σ3, being built centered on initial point o in the coordinate system of step (3)
A vertical ellipsoid, makes this ellipsoid cover n/2 carrier stop position;
Wherein, ellipsoid meets equation in the half shaft length a in longitudinal bias direction
Ellipsoid meets equation in the half shaft length b in lateral deviation direction
Ellipsoid meets equation in the half shaft length c in height tolerance direction
WhereinRepresent that the n position that carrier stops falls in longitudinal direction
Probability in [- a ,+a] in bias direction;Represent that n carrier stops
The probability that falls on lateral deviation direction [- b ,+b] of n position;
Represent the probability that the n position that carrier stops falls on height tolerance direction [- c ,+c].
Obtain after solution
A=1.539 σ1
B=1.539 σ2
C=1.539 σ3
Wherein θ, β and r are respectively polar parameter.
(5) this inertial navigation system three-dimensional impact accuracy esep=(a, b, c) is obtained according to the ellipsoid that step (4) is set up,
Wherein a, b, c be respectively the ellipsoid set up of step (4) longitudinally, laterally, the half shaft length in height tolerance direction.
The inventive method can be applied in the three dimensions impact accuracy of guided missile is estimated, for example, when being existed according to this method
The impact accuracy of certain guided missile gathers the drop point site of 1000 guided missiles in estimating, be then passed through being calculated σ1=30m, σ2=
10m、σ3=40m, sets up coordinate system oxyz according to said method, and the set of locations center making the drop point site of 1000 guided missiles is former
Point o, the direction of longitudinal bias falls in x-axis, and the direction of height tolerance falls in z-axis, sets up ellipsoid, make this in this coordinate system
Ellipsoid covers the drop point site of 500 guided missiles, and three semiaxis of the ellipsoid set up are respectively 46.17m, 15.39m, 61.56m,
I.e. esep=(46.17m, 15.39m, 61.56m), wherein, 46.17m is the half shaft length in longitudinal bias direction for the ellipsoid,
15.39m is the half shaft length in lateral deviation direction for the ellipsoid, and 61.56m is the half shaft length in height tolerance direction for the ellipsoid.Mould
The distribution that plan target practice is 1000 times is as shown in Figure 4.It can be seen that describing distribution with ellipsoid relatively compared with Fig. 2
Accurately, and sep assessment result can not accurate description drop point feature.
The non-detailed description of the present invention is known to the skilled person technology.
Claims (1)
1. a kind of inertial navigation system three dimensions impact accuracy method of estimation is it is characterised in that step is as follows:
(1) the n positional information that after the n navigation of collection inertial navigation system, carrier stops, n position constitutes set of locations, and n is extremely
It is 6 less;
(2) the n positional information being obtained according to step (1), obtains the coordinate at set of locations center, then measures each position relatively
In the longitudinal bias at set of locations center, lateral deviation and height tolerance, it is calculated the mark of all longitudinal bias on this basis
Quasi- difference σ1, the standard deviation sigma of lateral deviation2Standard deviation sigma with height tolerance3;
(3) with set of locations center for initial point o, set up rectangular coordinate system in space oxyz, so that the direction of longitudinal bias is fallen in x-axis or y
On one of axle of axle, the direction of height tolerance falls in z-axis;
(4) standard deviation sigma of the longitudinal bias being obtained according to step (2)1, the standard deviation sigma of lateral deviation2Standard with height tolerance
Difference σ3, setting up an ellipsoid centered on initial point o in the coordinate system of step (3), make this ellipsoid cover n/2 carrier stopping
Position;
The ellipsoid set up meets following form:
Ellipsoid meets equation in the half shaft length a in longitudinal bias direction
Ellipsoid meets equation in the half shaft length b in lateral deviation direction
Ellipsoid meets equation in the half shaft length c in height tolerance direction
WhereinRepresent that the n position that carrier stops falls in longitudinal direction
Probability in [- a ,+a] in bias direction;Represent what carrier stopped
The probability that n position falls on lateral deviation direction [- b ,+b];
Represent the probability that the n position that carrier stops falls on height tolerance direction [- c ,+c];
Obtain after solution
A=1.539 σ1
B=1.539 σ2
C=1.539 σ3
Wherein θ, β and r are respectively polar parameter;
(5) this inertial navigation system three-dimensional impact accuracy esep=(a, b, c) is obtained according to the ellipsoid that step (4) is set up, thus
Complete the estimation of inertial navigation system three-dimensional impact accuracy.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410397823.2A CN104197954B (en) | 2014-08-13 | 2014-08-13 | Method for estimating precision of drop points of inertial navigation system in three-dimensional space |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410397823.2A CN104197954B (en) | 2014-08-13 | 2014-08-13 | Method for estimating precision of drop points of inertial navigation system in three-dimensional space |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104197954A CN104197954A (en) | 2014-12-10 |
CN104197954B true CN104197954B (en) | 2017-01-18 |
Family
ID=52083282
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410397823.2A Active CN104197954B (en) | 2014-08-13 | 2014-08-13 | Method for estimating precision of drop points of inertial navigation system in three-dimensional space |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104197954B (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104848857B (en) * | 2015-04-30 | 2017-11-28 | 北京航天控制仪器研究所 | Ballistic missile inertial measurement system precision index auto-allocation method |
CN107270943B (en) * | 2017-07-28 | 2020-02-14 | 北京航天控制仪器研究所 | Two-dimensional correlation inertial navigation drop point precision evaluation method |
CN108548537B (en) * | 2018-02-11 | 2020-09-18 | 北京航天控制仪器研究所 | Estimation method and system for describing navigation drop point precision by six-parameter ellipsoid probability error |
CN108458727B (en) * | 2018-02-11 | 2020-09-18 | 北京航天控制仪器研究所 | Precision index self-adaptive distribution method and system for ballistic missile inertia measurement system |
CN108692731B (en) * | 2018-05-22 | 2020-11-20 | 北京航天控制仪器研究所 | Non-equal standard error and zero-mean sphere probability error representation method |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1999030102A1 (en) * | 1997-12-04 | 1999-06-17 | Eurocopter Deutschland Gmbh | Device and method for determining the impact point of a ballistic missile |
EP1291600A1 (en) * | 2001-09-07 | 2003-03-12 | Tda Armements S.A.S. | Method for guiding a device, especially a munition |
-
2014
- 2014-08-13 CN CN201410397823.2A patent/CN104197954B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1999030102A1 (en) * | 1997-12-04 | 1999-06-17 | Eurocopter Deutschland Gmbh | Device and method for determining the impact point of a ballistic missile |
EP1291600A1 (en) * | 2001-09-07 | 2003-03-12 | Tda Armements S.A.S. | Method for guiding a device, especially a munition |
Non-Patent Citations (4)
Title |
---|
导弹落点精度的鉴定方法-概率圆方法;程光显等;《国防科技大学学报》;20010531;第23卷(第5期);全文 * |
导弹落点精度评估的仿真分析;张士峰等;《***仿真学报》;20021130;第14卷(第11期);全文 * |
弹道导弹落点精度分析与仿真方法研究;王肖飞等;《海军航空工程学院学报》;20060131;第21卷(第1期);全文 * |
空间预警***落点预报精度分析;童丽等;《导弹与航天运载技术》;20070430(第2期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN104197954A (en) | 2014-12-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104197954B (en) | Method for estimating precision of drop points of inertial navigation system in three-dimensional space | |
CN104197953B (en) | A kind of inertial navigation system bidimensional impact accuracy method of estimation | |
CN103076618B (en) | Method for performing ground check and performance evaluation on satellite-borne measurement-type GNSS (Global Navigation Satellite System) receiver | |
CN103486984B (en) | The detection method of profile right alignment in a kind of wind-tunnel | |
CN102707726B (en) | Target location method of unmanned plane | |
CN103954953B (en) | The blind source error compensation method of a kind of airborne laser radar based on data-driven | |
CN104316079B (en) | Drop point precision estimation method for inertia measurement system based on rocket sled test | |
CN108061889A (en) | AIS and the correlating method of radar angular system deviation | |
CN110221244B (en) | Robust positioning method based on arrival time difference under non-line-of-sight condition | |
CN106446422A (en) | Log likelihood estimation based novel passive locating and tracking method | |
CN103235304B (en) | SAR (synthetic aperture radar) geometric correction method for modifying error equivalent RD (range-Doppler) model | |
CN106934113A (en) | Suitable for the modeling method of the improved polynomial of the vertical total electron content modeling in region ionosphere | |
CN105509765A (en) | Inertial/DVL/USBL installation error calibration method | |
CN104066179A (en) | Improved method for positioning WSN nodes through adaptive iterative UKF | |
CN106199580A (en) | A kind of Singer model refinement algorithm based on fuzzy inference system | |
CN105021199B (en) | Multi-model self-adapting method for estimating state and system based on LS | |
CN105046046A (en) | Ensemble Kalman filter localization method | |
CN107270943B (en) | Two-dimensional correlation inertial navigation drop point precision evaluation method | |
CN108917698B (en) | Azimuth angle calculation method | |
CN103954949A (en) | Multi-model oval Hough transformation accumulation detection method for invisible ski-jump maneuvering target | |
CN107621632A (en) | Adaptive filter method and system for NSHV tracking filters | |
CN106353756A (en) | Descending track spotlight SAR (synthetic aperture radar) positioning method based on image matching | |
CN106018201A (en) | Mixed field particle size measuring method based on average filtering | |
CN105466391A (en) | A tower base section generating method by utilization of a digital elevation model and field-data correction | |
CN110221245B (en) | Robust TDOA (time difference of arrival) positioning method for jointly estimating target position and non-line-of-sight error |
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 |