CN105716605A - Matching method of gravity-aided inertial navigation system - Google Patents

Matching method of gravity-aided inertial navigation system Download PDF

Info

Publication number
CN105716605A
CN105716605A CN201610192227.XA CN201610192227A CN105716605A CN 105716605 A CN105716605 A CN 105716605A CN 201610192227 A CN201610192227 A CN 201610192227A CN 105716605 A CN105716605 A CN 105716605A
Authority
CN
China
Prior art keywords
point
candidate match
inertial navigation
gravity
track
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
Application number
CN201610192227.XA
Other languages
Chinese (zh)
Other versions
CN105716605B (en
Inventor
邓志红
卢文典
付梦印
王博
肖烜
刘彤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Technology BIT
Original Assignee
Beijing Institute of Technology BIT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201610192227.XA priority Critical patent/CN105716605B/en
Publication of CN105716605A publication Critical patent/CN105716605A/en
Application granted granted Critical
Publication of CN105716605B publication Critical patent/CN105716605B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/10Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
    • G01C21/12Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
    • G01C21/16Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Automation & Control Theory (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Navigation (AREA)

Abstract

The invention discloses a matching method of a gravity-aided inertial navigation system. The manner that search regions can be changed is adopted, the search regions change in real time along with changes of inertial navigation precision, the searching range is effectively reduced, the number of alternative matched track points is reduced, and the searching efficiency is improved; the alternative matched track points are further screened by looking for points approaching a gravity abnormal value measured by a gravity meter in each search region; the number of the alternative matched tracks is effectively reduced under the constraint condition that the distance between each matched tracking point and the corresponding inertial navigation output track point deviates increasingly and changes in a certain range, and the matching efficiency is improved.

Description

A kind of Method in Gravity Aided INS system matches method
Technical field
The invention belongs to integrated navigation technology field, relate to a kind of Method in Gravity Aided INS system matches method.
Background technology
As the big country that land-sea has both, "Oceanic" strategy status in China's overall development strategy is more and more important, China also increasingly payes attention to the development of marine technology, and underwater hiding-machine technology is indispensable in the development of marine technology, has highly important strategic importance.Inertial navigation system relies solely on the feature that self can realize navigation by feat of it, usually used as the prime navaid system of underwater hiding-machine.Without from outer received signal during inertial navigation system work, without to external radiated signal, there is disguised high, the passive advantage such as independently.But inertial navigation is as a kind of integral form navigation mode, and error can constantly be accumulated in time, and its positioning precision can not meet the requirement of latent device during long boat.Therefore rely solely on inertial navigation system and cannot realize precisely navigating when growing boat of underwater hiding-machine, adopt other navigation mode aided inertial navigation systems to promote the trend that navigation accuracy is a certainty.
Geophysical field (landform field, earth's magnetic field gravitational field etc.) information is the intrinsic information of the earth, is often referred to the information of the various physical phenomenons measured by various geophysical methods on earth's surface or adjacent ground surface.Utilize the technology that Geophysical Field Information realizes navigation to cause the great attention of people gradually by feat of its passive autonomous feature, and be able to develop rapidly at aided inertial navigation system regions.The conventional mode utilizing Geophysical Field Information aided inertial navigation has terrain aided inertial navigation system, earth's magnetic field aided inertial navigation system and Method in Gravity Aided INS system.For underwater navigation, terrain data measurement difficulty of getting up is very big, earth's magnetic field is easily not high by the interference such as ferromagnetic material and stability, and gravity field data is prone to measure and sufficiently stable, therefore Method in Gravity Aided INS system can greatly meet the requirement of underwater hiding-machine " autonomy, in high precision, disguised ".
Gravity matching algorithm is one of core technology of Method in Gravity Aided INS system, ultimate principle is that the gravimetric data recorded by real-time gravimeter is compared by certain data processing method with the gravimetric data being previously stored in gravimetric database, and the similar of two groups of data or degree of correlation is determined by decision criteria, finally determine final Optimum Matching sequence (point).Gravity matching algorithm is broadly divided into two big classes according to sample mode, and single-point mates and sequences match.Single-point matching algorithm adopts the mode that each sampled point carries out Gravity Matching, suppresses increasing of inertial navigation error with this, and conventional is Sang Diya (SITAN) matching algorithm proposed by U.S.'s Sang Diya laboratory.This algorithm utilizes EKF technology to realize, matching process is realized by gravitational field being carried out local linearization process, but the shortcoming of this algorithm is: non-linear due to gravity field feature, in the obvious region of gravity field feature, EKF linearized stability is bigger, filtering divergence is caused so that coupling loses meaning time serious.Sequence correlation matching algorithm adopts the mode just carrying out a Gravity Matching when sample sequence is accumulated to certain length, and inertial navigation error is modified.Conventional sequence correlation matching algorithm includes nearest contour iterative algorithm (ICCP), relevant extreme value matching algorithm etc..ICCP algorithm comes from image registration algorithm, and relevant extreme value matching algorithm is by the terrain contour matching TERCOM algorithm development in terrain match.Sequences match arithmetic accuracy is high, but owing to coupling each time will carry out after sampling enough counting, so real-time is poor.
Summary of the invention
In view of this, it is an object of the invention to provide a kind of Method in Gravity Aided INS system matches method, with the problem solving conventional sequence matching algorithm poor real.When carrying out sequences match every time, adopt the mode in variable search region, make the size precision real-time change along with inertial navigation of region of search, effectively reduce hunting zone, reduce the quantity of candidate match tracing point;In region of search, find the approximately equalised point of gravimetric data measured with gravimeter successively, complete the screening of candidate match tracing point further;Utilize the range deviation between coupling tracing point and corresponding inertial navigation output trajectory point increasing and change this constraints within the specific limits, effectively reduce the quantity of candidate match track;Finally consider that the error of inertial navigation system is constantly accumulated in time, therefore introduce weighting performance indications, give bigger weights to the data that accuracy and confidence is higher in the matching process.
A kind of Method in Gravity Aided INS system matches method, comprises the steps:
Under step one, off-line state, the gravity datum figure prestored finds the inertial navigation system output trajectory point of correspondence;By the inertial navigation system output trajectory point of acquisition with N number of for one group, it is divided into and organizes matching sequence more;The span of described N is 5 to 20;
Step 2, from navigation start time, select first group of matching sequence, respectively centered by each point in this matching sequence, one by one delimit square aearch region;In square aearch region, each point except inertial navigation system output trajectory point is called candidate match tracing point;
Step 3, for each inertial navigation system output trajectory point, the GRAVITY ANOMALIES measured by the gravimeter of its correspondence asks poor with the GRAVITY ANOMALIES of each candidate match tracing point in the square aearch region at this inertial navigation system output trajectory point place, if difference is less than predetermined threshold value δ, then it is assumed that candidate match tracing point and corresponding inertial navigation system output trajectory point approximate gravity equivalent point each other;Being weeded out by the candidate match tracing point being unsatisfactory for approximate gravity equivalent point condition, tracing point is alternately mated in the continuation satisfied condition;
Step 4: after the candidate match tracing point in each square aearch region has screened, arbitrarily selects a candidate match tracing point in each square aearch region, sequentially in time one candidate match track of composition;As stated above, travel through each candidate match tracing point, obtain all of candidate match track;
For any one in all candidate match tracks, calculate the distance between candidate match tracing point and corresponding inertial navigation output trajectory point on this candidate match track sequentially in time one by one;
Then the candidate match track meeting following matching constraint condition is screened, described matching constraint condition is: elapse in time, in same candidate match track, distance between each candidate match tracing point and corresponding inertial navigation system output trajectory point is increasing, meanwhile, the distance that two adjacent candidate match tracing points are corresponding is less than or equal to default threshold value;
Step 5, for each the candidate match track screened, calculate the GRAVITY ANOMALIES g that each candidate match tracing point is correspondingi(xk,yk) the GRAVITY ANOMALIES g that records with real-time gravimeteri(xi,yi) difference square, then be weighted summation, the weighted value of i-th candidate match tracing point is μi, obtaining error accumulation Δ A, computing formula is:
Δ A = Σ i = 1 N μ i [ g i ( x i , y i ) - g i ( x k , y k ) ] 2
Candidate match track when finding error accumulation Δ A to take minima, is best coupling track.
Step 6, method according to step 2 to step 5, each group of matching sequence that step one is divided processes, and obtains the best coupling track each group corresponding, ultimately forms complete coupling track output.
It is also preferred that the left the length of side in described square aearch region accumulates continuous increase in time.
It is also preferred that the left accumulate in time, described weight is more and more less.
It is also preferred that the left half length of side in each square aearch region is equal to the drift distance of corresponding inertial navigation output trajectory point.
There is advantages that
(1) algorithm of the present invention utilizes the mode in variable search region, makes region of search with inertial navigation precision real-time change, effectively reduces hunting zone, decrease the quantity of candidate match tracing point, and search efficiency is promoted.
(2) algorithm of the present invention by finding the approximately equalised point of GRAVITY ANOMALIES measured with gravimeter in each region of search, completes the screening of candidate match tracing point further;Utilize the range deviation between coupling tracing point and corresponding inertial navigation output trajectory point increasing and change this constraints within the specific limits simultaneously, effectively reduce the quantity of candidate match track, improve matching efficiency.
Accompanying drawing explanation
Fig. 1 is the matching process flow chart of the present invention;
Fig. 2 is the gravitational field Background (containing planning flight path and inertial navigation track) in the embodiment of the present invention;
Fig. 3 is the region of search adopting traditional method to obtain;
Fig. 4 is the region of search adopting the method for the present invention to obtain;
Fig. 5 is the approximately equivalent point adopting the method for the present invention to find;
Fig. 6 is the candidate match track adopting the method for the present invention to find;
Fig. 7 is matching result in the embodiment of the present invention.
Detailed description of the invention
Develop simultaneously embodiment below in conjunction with accompanying drawing, describe the present invention.
A kind of Method in Gravity Aided INS system matches method, comprises the steps:
Under step one, off-line state, the gravity datum figure prestored finds the inertial navigation system output trajectory point of correspondence;By inertial navigation system output trajectory point with N number of for one group, it is divided into and organizes matching sequence more;The span of described N is 5 to 20;
Step 2, from navigation start time, select first group of matching sequence, centered by each point in this matching sequence, with inertial navigation system offset distance for according to determining square aearch region one by one.In square aearch region, each point except inertial navigation system output trajectory point is called candidate match tracing point.When underwater hiding-machine navigated by water in the starting stage, inertial navigation precision is higher, and what now can the length of side in square aearch region be arranged is less.Along with inertial navigation error is constantly accumulated in time, cannot accurately find coupling tracing point in the initial region of search set, now the length of side in square aearch region can be become larger.Wherein, for ensureing to search coupling tracing point, wherein, half length of side in each square aearch region is equal to the drift distance of corresponding inertial navigation output trajectory point.The mode in variable search region effectively reduces hunting zone, and the candidate match point quantity in region of search also greatly reduces.
Step 3, in theory, after overmatching, a coupling tracing point can be found in each moving-square search region, the GRAVITY ANOMALIES that the GRAVITY ANOMALIES that this coupling tracing point reads on gravity datum figure records with corresponding real-time gravimeter should be equal, but owing to there is real-time gravimeter measurement error equal error item, the two should be approximately equivalent relation.
Utilize this feature, for inertial navigation system output trajectory point, the GRAVITY ANOMALIES g measured by gravimeteri(xi,yi) the GRAVITY ANOMALIES g of candidate match point each with in corresponding square aearch regioni(xk,yk) ask poor, if | gi(xi,yi)-gi(xk,yk) |=Δ gi≤ δ (δ is the positive number meeting navigation request, rule of thumb sets), then it is assumed that candidate match point k and corresponding inertial navigation system output trajectory point i is approximate gravity equivalent point each other.Being weeded out by the candidate match tracing point being unsatisfactory for approximate gravity equivalent point and asking for condition, tracing point is alternately mated in the continuation satisfied condition.
Step 4: after the candidate match point in each square aearch region has screened, arbitrarily selects a candidate match point in each square aearch region, sequentially in time, and candidate match point one the candidate match track of composition that will be singled out;As stated above, travel through each candidate match point, obtain all of candidate match track.
Owing to inertial navigation error is constantly accumulated in time, therefore the range deviation between coupling tracing point and corresponding inertial navigation output trajectory point should be increasing.Simultaneously taking account of inertial navigation system and drift about not too large at short notice, therefore change can in certain scope.
From all candidate match tracks obtained above, arbitrarily select one, calculate the distance between candidate match tracing point and corresponding inertial navigation output trajectory point on this candidate match track sequentially in time one by one;
Assume tiMoment is first candidate match point of this candidate match track, tiThe longitude of moment inertial navigation output trajectory point i and latitude are respectivelyWithtiThe longitude of the candidate match tracing point A that moment correspondence inertial navigation output trajectory point i extracts and latitude are respectivelyWithti+1The longitude of moment inertial navigation output trajectory point i+1 and latitude are respectivelyWithti+1The longitude of the candidate match tracing point B that moment correspondence inertial navigation output trajectory point i+1 extracts and latitude are respectivelyWith(k=1,2 ..., Ni+1);The distance of inertial navigation output trajectory point i and candidate match tracing point A is δi;The distance of inertial navigation output trajectory point i+1 and candidate match tracing point B is δi+1;Earth radius is R;ε is a small distance meeting navigation request.Introduce matching constraint condition, in same candidate match track, whether the distance between each candidate match point and corresponding inertial navigation output trajectory point is increasing, simultaneously, distance corresponding to 2 adjacent candidate match points, whether less than or equal to default threshold epsilon, is namely expressed as following formula:
0≤δi+1i≤ε
Wherein
Special instruction, if the warp of 2, latitude are respectively on the earthWithThen between 2, the shortest circular arc distance is:
The candidate match track meeting above-mentioned matching constraint condition is screened.
Step 5, for each the candidate match track screened, calculate the GRAVITY ANOMALIES g that each candidate match point is correspondingi(xk,yk) the GRAVITY ANOMALIES g that records with real-time gravimeteri(xi,yi) difference square, then be weighted summation, weighted value is μi, obtaining error accumulation Δ A, computing formula is:
Δ A = Σ i = 1 N μ i [ g i ( x i , y i ) - g i ( x k , y k ) ] 2
Owing to inertial navigation error is accumulated in time, the data that export than later moment in time of data of previous moment output are compared, and are generally of higher accuracy and confidence, therefore in the matching process the weighted value of previous moment more than later moment in time weighted value, namely accumulating in time, weight is more and more less.Optimum Matching design criteria is to make error accumulation Δ A take minima, obtains best coupling track.
Step 6, method according to step 2 to step 5, each group of matching sequence that step one is divided processes, and obtains the best coupling track each group corresponding, ultimately forms complete coupling track output.
Embodiment:
In the present embodiment, simulated conditions is: the longitude interval of matching area is 138.3-139.3 degree, and latitude interval is 25.1-25.4 degree;Gravity map resolution 0.5' × 0.5';One group of sequence takes 10 points.Gravitational field Background (containing planning flight path and inertial navigation track) is as shown in Figure 2.
The region of search such as Fig. 3 utilizing traditional method (fixing search region) to determine, utilizes the method (variable search region) of the present invention such as Fig. 4, it is seen that hunting zone effectively reduces.Fig. 4 compares Fig. 3, and the total number of candidate match tracing point in region of search is by 104Individual effective minimizing is to 103Individual.
Behind limit search region, utilize | gi(xi,yi)-gi(xk,yk) |=Δ gi≤ δ selects the approximate gravity equivalent point of correspondence, as shown in Figure 5.The process of choosing of approximate gravity station further reduces the number of candidate match tracing point, has compressed it 102The order of magnitude.
After selecting approximate gravity equivalent point, just complete the screening of candidate match point.Candidate match track can be obtained by extracting candidate match point.Utilize coupling tracing point apart from deviation this constraints increasing of corresponding inertial navigation output point, i.e. δi+1i>=0, decrease the quantity of candidate match track.After seeking out approximately equivalent point, candidate match track has 1010Bar, and after utilizing the constraint of this constraints, candidate match track is effectively reduced to 105Bar, as shown in Figure 6.
Finally utilizing the weighting performance indications that the inventive method proposes, select candidate match track when Δ A takes minima as final matching results, matching result is as shown in Figure 7.
Simulation result: in this matching sequence, coupling tracing point latitude error is 0.42' to the maximum, and longitude error is 0.46' to the maximum, is respectively less than a grid, it is seen that matching result precision is significantly high.In coupling real-time and high efficiency, by introducing the mode in variable search region, by the candidate match total number of point in region of search by 104Individual minimizing is to 103Individual;Choose process by what introduce approximate gravity equivalent point, further the number of candidate match point is compressed to 102The order of magnitude;After seeking out approximate gravity equivalent point, alternate trajectory has 1010Bar, by introducing matching constraint condition, alternate trajectory effectively reduces to 105Bar.
As can be seen here, compared with the inventive method extreme value matching algorithm relevant with tradition, it is ensured that be effectively improved matching efficiency while matching precision is significantly high.
In sum, these are only presently preferred embodiments of the present invention, be not intended to limit protection scope of the present invention.All within the spirit and principles in the present invention, any amendment of making, equivalent replacement, improvement etc., should be included within protection scope of the present invention.

Claims (4)

1. a Method in Gravity Aided INS system matches method, it is characterised in that comprise the steps:
Under step one, off-line state, the gravity datum figure prestored finds the inertial navigation system output trajectory point of correspondence;By the inertial navigation system output trajectory point of acquisition with N number of for one group, it is divided into and organizes matching sequence more;The span of described N is 5 to 20;
Step 2, from navigation start time, select first group of matching sequence, respectively centered by each point in this matching sequence, one by one delimit square aearch region;In square aearch region, each point except inertial navigation system output trajectory point is called candidate match tracing point;
Step 3, for each inertial navigation system output trajectory point, the GRAVITY ANOMALIES measured by the gravimeter of its correspondence asks poor with the GRAVITY ANOMALIES of each candidate match tracing point in the square aearch region at this inertial navigation system output trajectory point place, if difference is less than predetermined threshold value δ, then it is assumed that candidate match tracing point and corresponding inertial navigation system output trajectory point approximate gravity equivalent point each other;Being weeded out by the candidate match tracing point being unsatisfactory for approximate gravity equivalent point condition, tracing point is alternately mated in the continuation satisfied condition;
Step 4: after the candidate match tracing point in each square aearch region has screened, arbitrarily selects a candidate match tracing point in each square aearch region, sequentially in time one candidate match track of composition;As stated above, travel through each candidate match tracing point, obtain all of candidate match track;
For any one in all candidate match tracks, calculate the distance between candidate match tracing point and corresponding inertial navigation output trajectory point on this candidate match track sequentially in time one by one;
Then the candidate match track meeting following matching constraint condition is screened, described matching constraint condition is: elapse in time, in same candidate match track, distance between each candidate match tracing point and corresponding inertial navigation system output trajectory point is increasing, meanwhile, the distance that two adjacent candidate match tracing points are corresponding is less than or equal to default threshold value;
Step 5, for each the candidate match track screened, calculate the GRAVITY ANOMALIES g that each candidate match tracing point is correspondingi(xk,yk) the GRAVITY ANOMALIES g that records with real-time gravimeteri(xi,yi) difference square, then be weighted summation, the weighted value of i-th candidate match tracing point is μi, obtaining error accumulation Δ A, computing formula is:
Δ A = Σ i = 1 N μ i [ g i ( x i , y i ) - g i ( x k , y k ) ] 2
Candidate match track when finding error accumulation Δ A to take minima, is best coupling track.
Step 6, method according to step 2 to step 5, each group of matching sequence that step one is divided processes, and obtains the best coupling track each group corresponding, ultimately forms complete coupling track output.
2. a kind of Method in Gravity Aided INS system matches method as claimed in claim 1, it is characterised in that the length of side in described square aearch region accumulates continuous increase in time.
3. a kind of Method in Gravity Aided INS system matches method as claimed in claim 1, it is characterised in that accumulating in time, described weight is more and more less.
4. a kind of Method in Gravity Aided INS system matches method as claimed in claim 1, it is characterised in that half length of side in each square aearch region is equal to the drift distance of corresponding inertial navigation output trajectory point.
CN201610192227.XA 2016-03-30 2016-03-30 A kind of Method in Gravity Aided INS system matches method Active CN105716605B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610192227.XA CN105716605B (en) 2016-03-30 2016-03-30 A kind of Method in Gravity Aided INS system matches method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610192227.XA CN105716605B (en) 2016-03-30 2016-03-30 A kind of Method in Gravity Aided INS system matches method

Publications (2)

Publication Number Publication Date
CN105716605A true CN105716605A (en) 2016-06-29
CN105716605B CN105716605B (en) 2018-07-10

Family

ID=56158484

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610192227.XA Active CN105716605B (en) 2016-03-30 2016-03-30 A kind of Method in Gravity Aided INS system matches method

Country Status (1)

Country Link
CN (1) CN105716605B (en)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106441297A (en) * 2016-08-30 2017-02-22 北京航空航天大学 Gravity error vector acquiring method and device of inertial navigation system
CN106723848A (en) * 2017-01-06 2017-05-31 安徽鼎晖新能源科技有限公司 A kind of anti-knapsack of wandering away of solar energy
CN108362281A (en) * 2018-02-24 2018-08-03 中国人民解放军61540部队 A kind of Long baselines underwater submarine matching navigation method and system
CN108896040A (en) * 2018-03-30 2018-11-27 中国空间技术研究院 Sky sea integrated water diving device inertia/gravity Combinated navigation method and system
CN109141426A (en) * 2018-08-10 2019-01-04 中国空间技术研究院 A kind of preferred method in subaqueous gravity matching navigation adaptation area
CN109557594A (en) * 2018-12-11 2019-04-02 中国人民解放军火箭军工程大学 Gravity datum figure time-varying modification method and system based on gravity anomaly time-varying
CN110906930A (en) * 2019-12-18 2020-03-24 中国人民解放军61540部队 AUV (autonomous underwater vehicle) combined underwater gravity lighthouse submarine navigation method and system
CN113503871A (en) * 2021-05-19 2021-10-15 北京理工大学 Gravity matching method based on correlation filtering
CN113624227A (en) * 2021-07-23 2021-11-09 北京理工大学 Contour line iteration gravity matching algorithm based on Mahalanobis distance
US11268813B2 (en) 2020-01-13 2022-03-08 Honeywell International Inc. Integrated inertial gravitational anomaly navigation system
CN114894164A (en) * 2022-04-08 2022-08-12 广州南方卫星导航仪器有限公司 Inclined image matching screening method and system

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102128625A (en) * 2010-12-08 2011-07-20 北京航空航天大学 Initial matching method for use in gravimetric map matching in gravity-aided inertial navigation system
CN102778230A (en) * 2012-06-14 2012-11-14 辽宁工程技术大学 Gravity gradient auxiliary positioning method of artificial physical optimization particle filtering
CN102788578A (en) * 2012-07-25 2012-11-21 中国人民解放军海军工程大学 Matching navigation method based on local gravity field approximation
CN103822634A (en) * 2014-02-24 2014-05-28 北京理工大学 ICCP algorithm-based gravity matching aided inertial navigation method
CN105043388A (en) * 2015-06-29 2015-11-11 中国船舶重工集团公司第七0七研究所 Vector search iterative matching method based on inertia/gravity matching integrated navigation
CN105180938A (en) * 2015-06-03 2015-12-23 北京理工大学 Particle filter-based gravity sampling vector matching positioning method

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102128625A (en) * 2010-12-08 2011-07-20 北京航空航天大学 Initial matching method for use in gravimetric map matching in gravity-aided inertial navigation system
CN102778230A (en) * 2012-06-14 2012-11-14 辽宁工程技术大学 Gravity gradient auxiliary positioning method of artificial physical optimization particle filtering
CN102788578A (en) * 2012-07-25 2012-11-21 中国人民解放军海军工程大学 Matching navigation method based on local gravity field approximation
CN103822634A (en) * 2014-02-24 2014-05-28 北京理工大学 ICCP algorithm-based gravity matching aided inertial navigation method
CN105180938A (en) * 2015-06-03 2015-12-23 北京理工大学 Particle filter-based gravity sampling vector matching positioning method
CN105043388A (en) * 2015-06-29 2015-11-11 中国船舶重工集团公司第七0七研究所 Vector search iterative matching method based on inertia/gravity matching integrated navigation

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
BO WANG ET AL.: "A Particle Filter-Based Matching Algorithm With Gravity Sample Vector for Underwater Gravity Aided Navigation", 《IEEE/ASM TRANSACTIONS ON MECHATRONICS》 *

Cited By (19)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106441297A (en) * 2016-08-30 2017-02-22 北京航空航天大学 Gravity error vector acquiring method and device of inertial navigation system
CN106441297B (en) * 2016-08-30 2019-11-01 北京航空航天大学 The gravity error vector acquisition methods and device of inertial navigation system
CN106723848A (en) * 2017-01-06 2017-05-31 安徽鼎晖新能源科技有限公司 A kind of anti-knapsack of wandering away of solar energy
CN108362281A (en) * 2018-02-24 2018-08-03 中国人民解放军61540部队 A kind of Long baselines underwater submarine matching navigation method and system
CN108362281B (en) * 2018-02-24 2020-11-24 中国人民解放军61540部队 Long-baseline underwater submarine matching navigation method and system
CN108896040A (en) * 2018-03-30 2018-11-27 中国空间技术研究院 Sky sea integrated water diving device inertia/gravity Combinated navigation method and system
CN108896040B (en) * 2018-03-30 2021-03-26 中国空间技术研究院 Inertia/gravity combined navigation method and system for sky-sea integrated underwater vehicle
CN109141426B (en) * 2018-08-10 2020-11-10 中国空间技术研究院 Method for matching navigation adaptation area by underwater gravity
CN109141426A (en) * 2018-08-10 2019-01-04 中国空间技术研究院 A kind of preferred method in subaqueous gravity matching navigation adaptation area
CN109557594A (en) * 2018-12-11 2019-04-02 中国人民解放军火箭军工程大学 Gravity datum figure time-varying modification method and system based on gravity anomaly time-varying
CN110906930A (en) * 2019-12-18 2020-03-24 中国人民解放军61540部队 AUV (autonomous underwater vehicle) combined underwater gravity lighthouse submarine navigation method and system
CN110906930B (en) * 2019-12-18 2021-08-27 中国人民解放军61540部队 AUV (autonomous underwater vehicle) combined underwater gravity lighthouse submarine navigation method and system
US11268813B2 (en) 2020-01-13 2022-03-08 Honeywell International Inc. Integrated inertial gravitational anomaly navigation system
CN113503871A (en) * 2021-05-19 2021-10-15 北京理工大学 Gravity matching method based on correlation filtering
CN113503871B (en) * 2021-05-19 2024-04-16 北京理工大学 Gravity matching method based on correlation filtering
CN113624227A (en) * 2021-07-23 2021-11-09 北京理工大学 Contour line iteration gravity matching algorithm based on Mahalanobis distance
CN113624227B (en) * 2021-07-23 2023-10-03 北京理工大学 Contour line iteration gravity matching algorithm based on mahalanobis distance
CN114894164A (en) * 2022-04-08 2022-08-12 广州南方卫星导航仪器有限公司 Inclined image matching screening method and system
CN114894164B (en) * 2022-04-08 2023-08-29 广州南方卫星导航仪器有限公司 Oblique image matching screening method and system

Also Published As

Publication number Publication date
CN105716605B (en) 2018-07-10

Similar Documents

Publication Publication Date Title
CN105716605B (en) A kind of Method in Gravity Aided INS system matches method
CN105760811B (en) Global map closed loop matching process and device
Wang et al. A particle filter-based matching algorithm with gravity sample vector for underwater gravity aided navigation
Han et al. An improved TERCOM-based algorithm for gravity-aided navigation
US8718927B2 (en) GPS data repair
US7890262B2 (en) Position estimation for navigation devices
CN102128625B (en) Initial matching method for use in gravimetric map matching in gravity-aided inertial navigation system
Han et al. A combined matching algorithm for underwater gravity-aided navigation
CN105953799B (en) Underwater vehicle path planning method in gravitational field adaptation area based on Information Entropy
WO2018008082A1 (en) Travel lane estimation system
CN104061932B (en) Method for navigation positioning by using gravitation vector and gradient tensor
CN106979778B (en) Positioning method, positioning device and mobile terminal
CN108896040B (en) Inertia/gravity combined navigation method and system for sky-sea integrated underwater vehicle
CN105043380A (en) Indoor navigation method based on a micro electro mechanical system, WiFi (Wireless Fidelity) positioning and magnetic field matching
CN102788578A (en) Matching navigation method based on local gravity field approximation
CN110118988A (en) A kind of combined positioning method and its device
CN103379619A (en) Method and system for positioning
CN105758408A (en) Method and device for building local maps
Chen et al. A new geomagnetic matching navigation method based on multidimensional vector elements of earth’s magnetic field
CN111307143B (en) Bionic navigation algorithm for multi-target evolution search based on geomagnetic gradient assistance
CN107894235A (en) Model error compensation method for autonomous navigation system of ultra-high-speed aircraft
CN110487276B (en) Sampling vector matching positioning method based on correlation analysis
CN109724592A (en) A kind of AUV earth magnetism bionic navigation method based on evolutionary gradient search
Han et al. A matching algorithm based on the nonlinear filter and similarity transformation for gravity-aided underwater navigation
Wang et al. A filtered-marine map-based matching method for gravity-aided navigation of underwater vehicles

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant