CN113156368B - Error parameter identification co-location method based on factor graph - Google Patents
Error parameter identification co-location method based on factor graph Download PDFInfo
- Publication number
- CN113156368B CN113156368B CN202110306773.2A CN202110306773A CN113156368B CN 113156368 B CN113156368 B CN 113156368B CN 202110306773 A CN202110306773 A CN 202110306773A CN 113156368 B CN113156368 B CN 113156368B
- Authority
- CN
- China
- Prior art keywords
- error
- parameter identification
- speed
- factor graph
- location
- 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
- 238000000034 method Methods 0.000 title claims abstract description 61
- 238000001914 filtration Methods 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims description 8
- 238000013461 design Methods 0.000 claims description 4
- 238000009795 derivation Methods 0.000 claims 1
- 238000005259 measurement Methods 0.000 abstract description 3
- 230000004927 fusion Effects 0.000 abstract 1
- 238000004088 simulation Methods 0.000 description 11
- 238000004422 calculation algorithm Methods 0.000 description 5
- 238000012546 transfer Methods 0.000 description 5
- 238000004891 communication Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000002159 abnormal effect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000012937 correction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000009189 diving Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 239000000178 monomer Substances 0.000 description 1
- 230000010363 phase shift Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/18—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves
Abstract
The invention provides a factor graph-based error parameter identification co-location method. Firstly, an error parameter identification factor graph model is established, then error parameter identification is carried out on process errors comprising speed errors and course errors, information is transmitted between function nodes and variable nodes in the factor graph by utilizing the maximum relevant entropy criterion, error compensation on the speed and the course of the slave boat is achieved, and finally filtering fusion estimation on position state information of the slave boat is achieved. Under the condition that the measurement accuracy of an inertial device in the system is not changed, the process error is identified and compensated, the co-location error is reduced, and the location capability of the co-location system is improved.
Description
Technical Field
The invention relates to an autonomous underwater vehicle (Autonomous Underwater Vehicle, AUV) co-location technology, in particular to a method for transmitting information between function nodes and variable nodes of a factor graph by adopting a maximum relevant entropy criterion, identifying a process error containing speed error and course error of a slave boat in the operation process of a co-location system, restraining the increase of the location error by compensating the process error, and improving the autonomous location capability of the slave boat.
Background
AUV has advantages such as the movable range is big, mobility is strong, intelligent degree is high, has great value in fields such as marine investigation, marine resource exploration, search and rescue under water, diving support, military investigation. For the ocean field, AUV related technology is a mainstream trend and development direction, and AUV cooperative operation can not only bear complex tasks which are difficult to compete with monomers, such as information collection, mine detection, coastal anti-diving, relay communication and the like, but also has the advantages of high efficiency, high reliability, high quality and the like. Therefore, AUV cooperative operation has wide application prospect.
The AUV cooperative positioning technology is a positioning mode based on a network, utilizes the underwater acoustic communication technology, fuses the relative position relation among AUVs by using the measurement information of the shared sensor to improve the positioning precision, and can effectively inhibit the influence of autonomous navigation positioning errors of dead reckoning on the cooperative positioning precision, so that the overall positioning error of a cooperative positioning system is bounded.
In an AUV co-location system, a common mode is a master-slave co-location structure comprising a master boat and a slave boat. The main boat is usually provided with high-precision navigation positioning equipment with high price, the high-precision position information of the main boat can be obtained, the precision of the navigation positioning equipment carried by the slave boat is low, the position estimation is firstly carried out by using a dead reckoning algorithm of the main boat, and then the correction is carried out by using the distance information between the main boat and the slave boat.
In an AUV co-location system, since ranging and communication are achieved using underwater acoustic devices, they are inevitably affected by characteristics of underwater acoustic signals. In an actual underwater environment, the phenomena of multipath propagation and abnormal measured value usually occur in underwater acoustic ranging and communication, and the influence of the multipath and abnormal value can cause the characteristic of non-Gaussian distribution to be presented in the underwater acoustic ranging, so that heavy tail noise with the tail part of probability distribution heavier than the Gaussian distribution occurs. However, the estimation performance of the traditional nonlinear co-location positioning methods, such as an extended kalman filtering method and an unscented kalman filtering method, when processing heavy-tail non-gaussian noise can be greatly reduced, and even the filter can be unstable.
In addition, for dead reckoning algorithms from boats, the key factors affecting autonomous positioning accuracy are speed error and heading error. Under the condition that the lateral movement of the AUV is not considered, when the position is calculated by using the measured speed and heading according to a mathematical model of the system, the speed error and the heading error have larger influence on the position estimation of the AUV, so that the error of self-positioning is increased.
Aiming at the problems, the invention designs an error parameter identification co-location method based on a factor graph, and the maximum correlation entropy criterion is adopted to transfer information between the function node and the variable node of the factor graph, so that the situation of heavy-tail non-Gaussian noise can be better processed. Meanwhile, an error parameter identification method is designed, the process error of the slave boat including the speed error and the course error in the operation process of the co-location system is identified, the increase of the location error is restrained through the compensation of the process error, the autonomous location capability of the slave boat is improved, and the purpose of improving the location precision of the co-location system is achieved.
Disclosure of Invention
The invention aims to design an error parameter identification co-location method based on a factor graph, and on the premise of not changing the precision of an inertial device, the information is transferred between a function node and a variable node of the factor graph by adopting a maximum relevant entropy criterion, so that the process error of a slave boat including a speed error and a course error in the operation process of a co-location system is identified, the increase of the location error is restrained by compensating the process error, and the autonomous location capability of the slave boat is improved.
The object of the invention can be achieved by the following steps:
step 1: establishing an error parameter identification co-location method factor graph model;
step 2: identifying the speed error and the course error of the system by utilizing the maximum correlation entropy criterion;
step 3: and compensating the speed error and the heading error, and further filtering and updating the state information of the system.
In step 1, since depth can be accurately measured by using the depth gauge mounted on the AUV, and accuracy improvement of depth information by co-positioning is limited, position information and distance observation information of each AUV can be projected into a horizontal plane, and the problem of co-positioning of the AUV in a two-dimensional plane can be discussed, namely, only horizontal coordinates (x, y) of the AUV need to be filtered. The factor graph contains two nodes: variable nodes and function nodes. As shown in figure 1, in order to consider an ideal factor graph model of speed error and heading error, a circle in the graph represents a variable node, a rectangle shape represents a function node, and each edge in the graph is connected with a variable node and a function node.
In an ideal case, the self-positioning according to FIG. 1 requires the use of measurementThe obtained speedAnd headingSubtracting the identified speed error +.>And heading error->Obtaining a speed and a course which are closer to the reality and participating in calculation; when the observation information is received from the boat to perform co-location to obtain posterior position estimation, the information is reversely transmitted to obtain the speed +.>And heading->And further estimate the corresponding error. But when the time k (k>l+1) receiving again the observation information and no observation information during the period l to k, the recognized speed error +_ is caused by the difficulty in obtaining a posterior position estimate during the period l to k>And heading error->And become complex. Typically the l-to-k time is short, so the following assumptions are made about speed error and heading error:
the positioning error at time l can be written as:
since no observation information is received from time l to time k, the position update adopts a priori estimation value, and the a priori estimation of x is as follows:
a priori estimates of y as
And (3) finishing to obtain:
then the positioning error at time k can be written as:
The factor graph model for identifying the speed error and the heading error can be redesigned according to the two formulas, and the factor graph model is shown in the figure 2. Meanwhile, the two formulas can be rewritten as follows:
further written as:
when (when)When the coefficient matrix is not full of rank, it is impossible to determine +.> and />But the manoeuvre of the AUV is small and the direction of travel does not change much during the co-location period, so the matrix is typically full rank.
In step 2, the invention designs a factor graph-based error parameter identification co-location method, and adopts the related entropy of variables to transfer information between variable nodes and function nodes of the factor graph. The solid line box in fig. 2 is a time update stage, that is, when the observed quantity sent by the master boat is not received from the slave boat; the part of the dashed box is the observation update phase, when the observations sent by the master are received from the boats and the posterior estimate of the position (x k ,y k ) And carrying out parameter identification on the speed error and the heading error of the slave boat as observation information. The specific update procedure is as follows.
(1) Time update
(1) Time update of speed error and heading error
(2) Time update of position estimate
In the calculation process, if no observed quantity is received at the time n, the time update of the position is utilized as the position estimation result at the current time, as follows:
x n =x′ n
y n =y′ n
(3) temporal updating of coefficients
(2) Observation update
In order to reduce the influence of heavy tail noise in observed quantity on speed and course error parameter identification, the method selects relevant entropy as a cost function in the error parameter identification co-location method based on the factor graph. In the observation update phase, the variable node is firstly used forThrough the function node R k To variable node->The transfer information may be expressed specifically as a process of maximizing the following expression.
Slave variable nodeThrough the function node S k To variable node->The transfer of information can be expressed in particular as a process that maximizes the following formula:
then there is at this time:
when the above is established, it is obvious that
Slave function node I k To variable node delta x k Can be expressed as:
slave function node L k To variable node delta y k Can be expressed as:
variable node δx k Can be expressed as:
when the correlation entropy in the above equation takes the maximum value, there are:
when the above formula is satisfied, it is apparent that:
variable node delta y k Can be expressed as:
when the correlation entropy in the above equation takes the maximum value, there are:
when the above formula is satisfied, it is apparent that:
when the correlation entropy in the above equation takes the maximum value, there are:
the method can be obtained after the above steps are unfolded:
the method can be obtained after finishing the above steps and phase shifting:
due to the right side of the equationAt this point, no update is made, and the calculation is replaced with a priori values, so there are:
thus, there may be:
when the correlation entropy in the above equation takes the maximum value, there are:
the method can be obtained after the above steps are unfolded:
the phase shift and arrangement can be achieved:
due to the right side of the equationAt this point, there is no update, and the calculation is replaced with a priori value, so there is:
thus, there may be:
when the error parameter identification co-location method based on the factor graph provided by the invention is used for error identification, two free parameters in the algorithm need to be determined, namely calculationValue and->Core width in the equation of value. In the method of the invention, since the speed error and heading error change in a short time are small, and the influence of the position error on the method of the invention is large, the large kernel width will cause the result of the algorithm to have large fluctuation, so the small kernel width should be selected, and the calculation is not needed>The kernel width of the value may be designed as:
drawings
FIG. 1 is a schematic diagram of an idealized factor graph model of an established error parameter identification co-location method.
FIG. 2 is a factor graph model of an established error parameter identification co-location method.
Fig. 3 is a diagram of the motion trajectories of the master and slave boats in a simulation experiment.
Fig. 4 is a range error in a simulation experiment.
FIG. 5 shows the speed and heading error recognition results in simulation experiments.
FIG. 6 shows the positioning errors of different algorithms in the simulation experiment.
Detailed Description
The present invention will be described in detail with reference to specific examples.
The invention provides an error parameter identification co-location method based on factor graphs, which adopts the maximum relevant entropy criterion to transfer information between function nodes and variable nodes of the factor graphs, identifies the process error containing speed error and course error of a slave boat in the operation process of a co-location system, inhibits the increase of the location error through the compensation of the process error, and improves the autonomous location capability of the slave boat. The invention aims at realizing the following steps:
1. establishing an error parameter identification co-location method factor graph model;
2. identifying the speed error and the course error of the system by utilizing the maximum correlation entropy criterion;
3. and compensating the speed error and the heading error, and further filtering and updating the state information of the system.
In order to verify the effectiveness of the invention, software is utilized to simulate the error parameter identification co-location method based on the factor graph.
The simulation conditions are as follows: in the simulation, errors of 0.1m/s and 3 degrees are respectively added to an electromagnetic log and a magnetic compass, gaussian white noise interference with standard deviation of 2m is added to underwater sound equipment, and 20m of error simulation heavy tail noise is added to 10% probability. The simulation time is set to 2000s, the cooperative positioning is carried out in the first 1000s in the experiment, and the positioning is carried out in the last 1000s from the boat only by means of self navigation equipment.
In the simulation, two main boats are shared, one secondary boat, the starting points of the two main boats are respectively 450m and 0m and 450m, and the starting points of the secondary boats are 0m and 450m, and the heading of each boat is 90 degrees and the speed is 2m/s. The actual running track of each boat is shown in figure 3. The range error in the simulation is shown in fig. 4, the identification result of the error parameter identification factor graph co-locating method in the simulation on the speed error and the course error is shown in fig. 5, and it can be seen from fig. 5 that the speed error and the course error are effectively identified in the system co-locating period by the method provided by the invention.
Fig. 6 shows positioning errors estimated by using different co-positioning methods, wherein a first method is a factor graph co-positioning method without error parameter identification, and a second method is a factor graph-based error parameter identification co-positioning method according to the present invention. As can be seen from the figure, the positioning error estimated by the method does not diverge with time before 1000s, but when the slave boat cannot receive the observed quantity sent by the master boat after 1000s, the slave boat performs autonomous positioning by means of the navigation equipment of the slave boat, and the positioning error diverges rapidly with the increase of time; and the second method recognizes and compensates the speed error and the course error, so that the positioning error of the ship after 1000s is far smaller than the positioning error estimated by the first method.
The effectiveness of the error parameter identification co-location method based on the factor graph provided by the invention is verified through the experiment, and the purpose of improving the overall location accuracy of the system can be achieved by identifying and compensating the speed error and the course error in the running process of the system on the premise of not improving the measurement accuracy of the inertial components in the co-location system.
Claims (2)
1. The error parameter identification co-location method based on the factor graph is characterized by comprising the following steps of:
step 1: establishing an error parameter identification co-location method factor graph model;
step 2: identifying the speed error and the course error of the system by utilizing the maximum correlation entropy criterion;
step 3: compensating the speed error and the course error, further filtering and updating the state information of the system,
in the step 2, in the process of identifying the speed error and heading error of the system by using the maximum correlation entropy criterion, information is transmitted in a factor graph by using the correlation entropy of the variable to obtain the speed errorAs shown in the formula:
wherein ,representing a priori estimates of the velocity error, the kernel width σ is chosen as follows:
wherein ,representing a priori estimates of heading error, the kernel width σ is selected as follows:
wherein ,δxk -representing the positioning error in the x-direction at time k;
δy k -representing the positioning error in the y-direction at time k;
2. The error parameter identification co-location method as defined in claim 1, wherein in the error parameter identification co-location method established in step 1, the measured velocity is used in an ideal factor graph modelAnd heading->Subtracting the identified speed error +.>And heading error->Obtaining the speed and heading which are closer to the true and participating in calculation, further filtering to obtain the position state information of the carrier,
in the factor model used in the design, the speed and the course after compensation obtained by theoretical derivation are shown in the following formula:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110306773.2A CN113156368B (en) | 2021-03-23 | 2021-03-23 | Error parameter identification co-location method based on factor graph |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110306773.2A CN113156368B (en) | 2021-03-23 | 2021-03-23 | Error parameter identification co-location method based on factor graph |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113156368A CN113156368A (en) | 2021-07-23 |
CN113156368B true CN113156368B (en) | 2023-05-30 |
Family
ID=76887996
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110306773.2A Active CN113156368B (en) | 2021-03-23 | 2021-03-23 | Error parameter identification co-location method based on factor graph |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113156368B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116608864B (en) * | 2023-07-19 | 2023-10-27 | 青岛哈尔滨工程大学创新发展中心 | AUV cooperative positioning method based on factor graph under influence of communication time delay |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1994028435A1 (en) * | 1993-05-28 | 1994-12-08 | Trimble Navigation Limited | Navigation system for adaptive weighting of gps and dead-reckoning information |
CN108489498A (en) * | 2018-06-15 | 2018-09-04 | 哈尔滨工程大学 | A kind of AUV collaborative navigation methods without mark particle filter based on maximum cross-correlation entropy |
CN110837854A (en) * | 2019-10-30 | 2020-02-25 | 东南大学 | AUV multi-source information fusion method and device based on factor graph |
CN112325885A (en) * | 2020-10-30 | 2021-02-05 | 哈尔滨工业大学 | Factor graph co-location algorithm based on mathematical statistical characteristics |
-
2021
- 2021-03-23 CN CN202110306773.2A patent/CN113156368B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1994028435A1 (en) * | 1993-05-28 | 1994-12-08 | Trimble Navigation Limited | Navigation system for adaptive weighting of gps and dead-reckoning information |
CN108489498A (en) * | 2018-06-15 | 2018-09-04 | 哈尔滨工程大学 | A kind of AUV collaborative navigation methods without mark particle filter based on maximum cross-correlation entropy |
CN110837854A (en) * | 2019-10-30 | 2020-02-25 | 东南大学 | AUV multi-source information fusion method and device based on factor graph |
CN112325885A (en) * | 2020-10-30 | 2021-02-05 | 哈尔滨工业大学 | Factor graph co-location algorithm based on mathematical statistical characteristics |
Non-Patent Citations (2)
Title |
---|
Fan Shiwei 等.Cooperative Positioning for Multi-AUVs Based on Factor Graph and Maximum Correntropy.IEEE Access.2019,第7卷153327-153337. * |
范世伟 等.基于因子图的协同定位与误差估计算法.系统工程与电子技术.2021,第43卷(第43期),499-507. * |
Also Published As
Publication number | Publication date |
---|---|
CN113156368A (en) | 2021-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109459040B (en) | Multi-AUV (autonomous Underwater vehicle) cooperative positioning method based on RBF (radial basis function) neural network assisted volume Kalman filtering | |
Morgado et al. | USBL/INS tightly-coupled integration technique for underwater vehicles | |
WO2022088797A1 (en) | Measurement abnormality-considered cooperative localization method for cluster type multi-deep-sea underwater vehicle | |
Bo et al. | Cooperative localisation of AUVs based on huber-based robust algorithm and adaptive noise estimation | |
CN106323291A (en) | Multi-robot cooperative location algorithm based on square root cubature Kalman filtering | |
CN110186461A (en) | A kind of collaborative navigation method based on gravity gradient information ranging | |
CN113156368B (en) | Error parameter identification co-location method based on factor graph | |
CN109765914A (en) | A kind of unmanned surface vehicle collision prevention method based on sliding window population | |
CN109460539A (en) | A kind of object localization method based on simplified volume particle filter | |
Fossen et al. | Exogenous kalman filter (xkf) for visualization and motion prediction of ships using live automatic identification systems (ais) data | |
Zhao et al. | Review of slam techniques for autonomous underwater vehicles | |
Havyarimana et al. | Achieving reliable intervehicle positioning based on redheffer weighted least squares model under multi-GNSS outages | |
Wang et al. | A novel adaptive sliding observation-based cooperative positioning algorithm under factor graph framework for multiple UUVs | |
CN112241589B (en) | Multi-underwater unmanned ship collaborative motion formation design method based on system observability degree | |
CN116026324B (en) | Cross-domain navigation system and method for water-air cross-medium craft | |
CN116182853A (en) | Ship self-adaptive combined navigation equipment and method | |
CN113008235B (en) | Multi-source navigation information fusion method based on matrix K-L divergence | |
CN113218421B (en) | Robust self-adaptive dynamic alignment method for strapdown inertial navigation system under Beidou rejection condition | |
CN115373383A (en) | Autonomous obstacle avoidance method and device for garbage recovery unmanned boat and related equipment | |
CN113687349A (en) | Unmanned ship sea surface target tracking method and device based on multi-sensor fusion | |
Zhang et al. | An outlier-robust Rao–Blackwellized particle filter for underwater terrain-aided navigation | |
CN112749481A (en) | Factor graph co-location algorithm based on maximum correlation entropy | |
Luo et al. | An optimal node depth adjustment method with computation-efficiency for target tracking in UWSNs | |
Chen et al. | Target Tracking of Navigation Radar for Unmanned Surface Vehicle Based on an Improved Adaptive Kalman Filtering | |
Moratuwage et al. | Collaborative multi-vehicle localization and mapping in high clutter environments |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |