CN117761745A - information optimization method for few-star positioning based on historical information fusion - Google Patents
information optimization method for few-star positioning based on historical information fusion Download PDFInfo
- Publication number
- CN117761745A CN117761745A CN202410195600.1A CN202410195600A CN117761745A CN 117761745 A CN117761745 A CN 117761745A CN 202410195600 A CN202410195600 A CN 202410195600A CN 117761745 A CN117761745 A CN 117761745A
- Authority
- CN
- China
- Prior art keywords
- satellite
- information
- time
- list
- combination
- 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
- 238000000034 method Methods 0.000 title claims abstract description 40
- 230000004927 fusion Effects 0.000 title claims abstract description 18
- 238000005457 optimization Methods 0.000 title claims abstract description 16
- 238000012216 screening Methods 0.000 claims abstract description 27
- 230000002028 premature Effects 0.000 claims abstract description 6
- 239000011159 matrix material Substances 0.000 claims description 9
- 238000005259 measurement Methods 0.000 claims description 4
- 238000001914 filtration Methods 0.000 claims 1
- 238000005070 sampling Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention provides an information optimization method for satellite-less positioning based on historical information fusion, and belongs to the field of satellite positioning. The method comprises the steps of firstly screening visible satellites at a single moment based on pseudo-range observation values by using a single-point information screening method, and selecting available satellite information at the single moment; then, the available satellite information obtained by single-point information screening is utilized to form a resolving list containing a historical information sequence, and the optimal resolving combination of positioning resolving is screened out from the resolving list; if the information combination newly added into the resolving list exists, redundant information and premature information in the resolving list are removed, and the optimal resolving list is ensured. Under the conditions that the number of available satellites is limited and the satellite signal interference is serious, satellite information with effective satellite signals and relatively good signal quality is optimized, and effective information utilization is realized; meanwhile, effective satellite information is introduced into fusion positioning, so that positioning accuracy and continuity are improved.
Description
Technical Field
The invention belongs to the field of satellite positioning, and particularly relates to an information optimization method for satellite-less positioning based on historical information fusion.
Background
At present, the outdoor positioning is mainly performed by a GNSS system mainly represented by a GPS, but in places such as underground passages, under overpasses, street tree shadows, urban canyons and the like, the GNSS cannot realize positioning due to signal blockage and failure to receive enough satellite signals. Therefore, some students develop positioning methods in the case of few stars, and the main research direction is to assist GNSS by using barometric altimetry, PDR, WIFI, vision, and the like. The learner tries to fuse the historical information of the satellite signals, so that continuous positioning under the condition of few satellites is realized, more information can be utilized to realize positioning with higher precision, and accumulated errors brought by methods such as PDR, vision and the like can be effectively assisted, so that the precision of continuous positioning is ensured. However, how to screen the fusion information in the positioning process under the condition that the number of available satellites is less than 4, so that the positioning with higher precision is realized, and no systematic method exists yet.
In the prior art, partial scholars adopt an optimal geometric error factor method to perform positioning calculation, the method traverses all visible satellites, and performs geometric precision factor solving on each possible combination condition respectively, so as to select satellite combinations with minimum geometric precision factor values. The method has strong reliability and can obtain good positioning precision; the method needs to traverse all possible satellite combinations, has complex operation and long time, and has very poor real-time performance especially when the satellite system is large and the number of visible satellites is large. Meanwhile, the method is effective based on the condition that the number of satellites is large and the satellite signal quality is good, and is not applicable to the condition that the number of satellites is less than 4.
The maximum tetrahedral volume method is adopted by partial scholars, and the method calculates the volume of the polyhedron formed by the visible satellite and the user, so that the combination with the maximum volume is selected, and the geometric precision factor value is the minimum. However, the method is similar to the optimal geometric error factor method, and the method also traverses all polyhedral combinations, so that the problems of complex operation and long time consumption are also caused. For less than 4 satellites, this approach is not applicable.
And the maximum determinant method is adopted by some scholars, and is characterized in that matrix determinant values formed by directional cosine are calculated, satellite combinations corresponding to the maximum determinant values are selected, and the geometric precision factor values are reduced along with the increase of determinant absolute values due to the inverse relation between the geometric precision factor values and determinant absolute values. Compared with the first two methods, the method reduces the calculation task, reduces the time consumption required by calculation, but still can not realize quick star selection. It is considered that effective information selection cannot be performed in the case where the number of satellites is limited and satellite signal interference is serious.
Disclosure of Invention
the invention provides an information optimization method for the less star positioning based on history information fusion, which is based on screening of fusion information of the less star positioning based on history information fusion and designs a corresponding optimization method, thereby realizing effective utilization of information and improving positioning accuracy and continuity.
in order to achieve the above purpose, the invention adopts the following technical scheme:
An information optimization method for a few-star positioning based on historical information fusion, comprising the steps of:
single-point information screening, namely screening visible satellites at a single moment by using a pseudo-range observation value, and selecting available satellite information at the single moment;
based on information optimization of the historical information, the available satellite information obtained by screening single-point information is used for forming a resolving list comprising a historical information sequence, and the optimal resolving combination of positioning resolving is screened out from the resolving list;
And eliminating redundant information, namely eliminating redundant information and premature information in the resolving list when information newly added into the resolving list exists.
the invention has the beneficial effects that:
(1) Under the conditions that the number of available satellites is limited and the satellite signal interference is serious, satellite information with effective satellite signals and relatively good signal quality is optimized, so that effective information utilization is realized;
(2) The effective satellite information is introduced into fusion positioning, so that the positioning accuracy and the continuity are improved.
Drawings
Fig. 1 is a schematic diagram of an information optimization method for a few-star positioning based on historical information fusion according to the present invention.
Detailed Description
The invention relates to an information optimization method for few-star positioning based on historical information fusion, which comprises the steps of selecting single-point available observation satellite information by using a single-point information screening method, and screening by mainly using pseudo-range observation values. And then optimizing the most resolved combined sequence for position resolving based on the information of the time sequence, finally removing redundant information, removing repeated information and premature information, and ensuring the optimization of the historical sequence. Fig. 1 is a schematic diagram of an information optimization method for the star-less positioning based on history information fusion, which specifically includes the following steps:
Step 1: single-point information screening, correcting satellite observation information when new satellite observation information is input at time k+1, and correcting pseudo-range value after correctionWherein/>Representing the true distance between satellite j and the user at time k+1,/>Is the seed difference of satellite receiver,/>Is a measurement error including multipath error and white noise; then calculate the difference between the corrected pseudorangesWherein/>corrected pseudo-range values at k+1st and k respectively,/>, respectivelyRespectively represent the real distance between the satellite j and the user at the k+1 time and the k time,/>Clock difference of satellite receiver at k+1th time and k time,/>, respectivelyrespectively represent the measurement errors of the k+1 time and the k time,/>satellite number for satellite j. Assuming that the user is moving at a low speed, for a user moving at a low speed, the distance moved by two consecutive sampling moments is limited. At the same time, the variation of the receiver seed difference between two consecutive sampling instants is smallThe difference between corrected pseudorange values may be equivalently/>Due to/>is measurement noise including multipath error, white noise, etc., so that the variation in two consecutive sampling moments will not be great, and the value will be at a certain threshold/>, within a certain timeIn range, therefore, using the difference between corrected pseudorange values/>To screen whether the information of satellite j at time k+1 is available, if soWherein/>Representing a difference determination threshold between corrected pseudorange values,/>Representing taking an absolute value, then the information of satellite j at time k+1 is considered to be available; screening all satellites at the k+1 time, and combining the screened information into satellite information combination/>Step 2 is entered; otherwise, repeating the step 1, and waiting for new observation information to be accessed. Wherein/>，/>representing the number of satellites available at time k+1,/>A satellite number combination indicating the available satellites at time k +1,Representing a combination of corrected pseudorange values for all available satellites.
step 2: updating the resolving list, and judging the resolving list formed by combining the satellite information at the kth momentif not, a new resolving list/> isestablished，/>Representing a matrix transpose; otherwise newly built solution list/>step 3 is entered.
step 3: updating satellite number sequence and judging satellite number sequenceIf it is empty, then use satellite number combination/>New satellite numbering sequence/>The method comprises the steps of carrying out a first treatment on the surface of the Otherwise, use satellite numbering combination/>And satellite numbering sequence/>Judging whether new available satellite numbers appear, namely judging satellite number combination/>, at the k+1 timeWhether or not satellite numbering sequence/>If there is no satellite number, the satellite number/>Is the satellite numbering sequence/>If there is no satellite number, the satellite number sequence is updated to/>I.e. satellite numbering sequence/>there are g+1 different satellite numbers, wherein/>And/>Representing the satellite number sequence/>, before time k+1satellite numbering in (3), and satellite numbering sequence/>)there are g satellite numbers and then go to step 4.
Step 4: satellite combination screening according to satellite number sequenceFrom the solution list/>Find satellite number/>satellite in solution list/>All time combinations of occurrence, constituting satellite number/>Time list of (2)Wherein/>Represents the/>Time,/>The time number, k+1 represents the k+1st time, and p represents the satellite number/>Is a sequence number of (2); screening the satellite numbers of the g+1 number to obtain a sequence related to the satellite numberTime list of all satellite numbers/>time list combination for composing all satellite numbersWherein satellite numbering sequence/>in (2), g+1 satellite numbers are given at the k+1 time, and the process goes to step 5.
Step 5: repeated combination rejection, for g+1 time listsScreening, if a time list/>, is presentSo thatDelete time list/>. After screening, a new time list combination/>, is formedAnd/>wherein v represents v combinations and v is not greater than the number of sequences in satellite/>Number g+1 of medium satellite numbers; step 6 is entered.
Step 6: screening optimal solution combinations, utilizing new time list combinationsThe corrected pseudo-range value, ephemeris information and other observation data at corresponding time are respectively combined into a solving combination to calculate a geometric precision factorWhere H is the observation matrix, tr [ ] is the trace operation,/>representing matrix inversion,/>representing the transpose of the observation matrix H; and selecting the smallest GDOP as an optimal solution combination to form a solution sequence for position solution, and entering step 7.
Step 7: redundant information is removed, repeated information is screened, and if satellite information combination exists, the method comprises the steps of selecting the redundant information, removing the redundant information, and screening the redundant informationSatisfy the following requirementsI.e./>Information redundancy is considered to be calculated from the solution list/>Deleting satellite information combinations/>While satellite information combination/>Become new satellite information combination/>. Wherein the method comprises the steps of，/>Representing satellite information combinations,/>indicating the number of satellites available at time i,Satellite number/>, representing the satellite available at time iCombination of/(I)corrected pseudorange values representing all available satellites/>Is a combination of (a) and (b). Step 8 is entered.
Step 8: premature information rejection, for solution listJudgment is made, and time/>, at time 1 is calculatedTime/>, with the current k+1 momentTime difference between/>If/>Wherein/>representing a time threshold, then from the solution list/>Satellite information combination/>, deleting time 1I.e. solution list/>Translating back 1 moment; if/>And returning to the step 1.
The above examples are preferred embodiments of the present invention, and the embodiments of the present invention are not limited to the above examples, and any other changes, modifications, substitutions, combinations, and simplifications that do not depart from the spirit and basic principles of the present invention should be made in the equivalent manner, and the present invention is also within the scope of the present invention.
Claims (4)
1. an information optimization method for a few-star positioning based on historical information fusion, comprising the steps of:
single-point information screening, namely screening visible satellites at a single moment by using a pseudo-range observation value, and selecting available satellite information at the single moment;
based on information optimization of the historical information, the available satellite information obtained by screening single-point information is used for forming a resolving list comprising a historical information sequence, and the optimal resolving combination of positioning resolving is screened out from the resolving list;
And eliminating redundant information, namely eliminating redundant information and premature information in the resolving list when information newly added into the resolving list exists.
2. The information optimization method for the star-less positioning based on the history information fusion according to claim 1, wherein the single-point information filtering specifically comprises:
When new satellite observation information is input at time k+1, correcting the satellite observation information, and correcting pseudo-range valuecalculating differences between corrected pseudorange valuesWherein/>respectively represent the corrected pseudo-range values at the k+1 time,/>, at the k timerespectively represent the k+1 time, the real distance between the satellite j and the user at the k time,/>respectively represent the k+1 time and the clock difference of the satellite receiver at the k time,/>Respectively represent the k+1 time, the measurement error including multipath error and white noise at the k time,/>satellite number for satellite j; with differences between corrected pseudorange values/>To screen whether the information of the satellite j at the k+1 moment is available, if soWherein/>Representing a difference determination threshold between corrected pseudorange values,/>Representing taking an absolute value, then the information of satellite j is considered to be available at time k+1; screening all satellite information at the k+1 time, and forming satellite information combination/> by the screened information，/>representing the number of satellites available at time k+1,/>Satellite number combination indicating available satellites at time k+1,/>Representing a combination of corrected pseudorange values for all available satellites.
3. An information preference method for a few stars positioning based on historical information fusion according to claim 1, wherein the information preference based on historical information specifically comprises:
Updating the resolving list, and judging the resolving list formed by combining the satellite information at the kth momentif not, a new resolving list/> isestablished，/>Representing a matrix transpose; otherwise newly built solution list/>；
Updating satellite number sequence and judging satellite number sequenceIf it is empty, then use satellite number combination/>New satellite numbering sequence/>The method comprises the steps of carrying out a first treatment on the surface of the Otherwise, use satellite numbering combination/>And satellite numbering sequence/>judging whether new available satellite numbers appear, namely judging satellite number combination/>, of available satellites at the k+1 timeWhether or not satellite numbering sequence/>If there is no satellite number, the satellite number/>Is the satellite numbering sequence/>If there is no satellite number, the satellite number sequence is updated to/>I.e. satellite numbering sequence/>there are g+1 different satellite numbers, wherein/>And/>Representing the satellite number sequence/>, before time k+1satellite numbering in (3), and satellite numbering sequence/>)G satellite numbers are included;
satellite combination screening according to satellite number sequenceFrom the solution list/>Find satellite number/>satellite in solution list/>All time combinations of occurrence, constituting satellite number/>Time list of (2)Wherein/>Represents the/>Time,/>The time number, k+1 represents the k+1st time, and p represents the satellite number/>is a sequence number of (2); screening the g+1 different satellite numbers to obtain a sequence/>, related to the satellite numbertime lists of all satellite numbers in the system, and time list combination/>Wherein satellite numbering sequence/>G+1 satellite numbers are given at the k+1 time;
repeated combination rejection, for g+1 time listsScreening, if a time list/>, is presentMake/>Delete time list/>After screening, a new time list combination/> isformedAnd/>wherein v represents v combinations and v is not greater than the number of sequences in satellite/>Number g+1 of medium satellite numbers;
screening optimal solution combinations, utilizing new time list combinationsthe corrected pseudo-range value and ephemeris information observation data at the corresponding time in the model are respectively combined into a solving combination to calculate the geometric precision factor/>Where H is the observation matrix, tr [ ] is the trace operation,/>representing matrix inversion,/>Representing the transpose of the observation matrix H; and selecting the smallest GDOP as the optimal solution combination to form a solution sequence for position solution.
4. The information optimization method for star-less positioning based on history information fusion according to claim 1, wherein the redundant information rejection specifically comprises:
Redundant information is removed, repeated information is screened, and if satellite information combination exists, the method comprises the steps of selecting the redundant information, removing the redundant information, and screening the redundant informationSatisfy the following requirementsI.e./>Information redundancy is considered to be calculated from the solution list/>Deleting satellite information combinations/>Wherein/>，/>representing the number of satellites available at time i,/>Satellite number/>, representing the satellite available at time iCombination of/(I)corrected pseudorange values representing all available satellites/>Is a combination of (a);
premature information rejection, for solution listJudgment is made, and time/>, at time 1 is calculatedTime/>, with the current k+1 momentTime difference between/>If/>Wherein/>representing a time threshold, then from the solution list/>Satellite information combination/>, deleting time 1I.e. solution list/>And translate 1 moment backward.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410195600.1A CN117761745B (en) | 2024-02-22 | 2024-02-22 | Information optimization method for few-star positioning based on historical information fusion |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410195600.1A CN117761745B (en) | 2024-02-22 | 2024-02-22 | Information optimization method for few-star positioning based on historical information fusion |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117761745A true CN117761745A (en) | 2024-03-26 |
| CN117761745B CN117761745B (en) | 2024-04-30 |
Family
ID=90326080
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410195600.1A Active CN117761745B (en) | 2024-02-22 | 2024-02-22 | Information optimization method for few-star positioning based on historical information fusion |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117761745B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102654407A (en) * | 2012-04-17 | 2012-09-05 | 南京航空航天大学 | Multiple-fault detecting device and detecting method for tightly-integrated inertial satellite navigation system |
| CN104181561A (en) * | 2013-05-24 | 2014-12-03 | 凹凸电子（武汉）有限公司 | Receiver and satellite positioning and speed measuring method |
| CN109932735A (en) * | 2019-03-25 | 2019-06-25 | 中国铁路设计集团有限公司 | The localization method of the short baseline single-frequency simple epoch solution of Beidou |
| CN113640838A (en) * | 2021-08-19 | 2021-11-12 | 中国电力科学研究院有限公司 | Multi-system combined satellite navigation positioning time service device and method |
| CN116399351A (en) * | 2023-04-23 | 2023-07-07 | 燕山大学 | Vehicle position estimation method |
| CN116859415A (en) * | 2023-05-23 | 2023-10-10 | 北京时代民芯科技有限公司 | Quick, stable and high-precision multi-fault satellite identification and positioning method |
| CN117040594A (en) * | 2023-08-10 | 2023-11-10 | 武汉大学 | Internet remote sensing satellite real-time service system oriented to mobile terminal user |
-
2024
- 2024-02-22 CN CN202410195600.1A patent/CN117761745B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102654407A (en) * | 2012-04-17 | 2012-09-05 | 南京航空航天大学 | Multiple-fault detecting device and detecting method for tightly-integrated inertial satellite navigation system |
| CN104181561A (en) * | 2013-05-24 | 2014-12-03 | 凹凸电子（武汉）有限公司 | Receiver and satellite positioning and speed measuring method |
| CN109932735A (en) * | 2019-03-25 | 2019-06-25 | 中国铁路设计集团有限公司 | The localization method of the short baseline single-frequency simple epoch solution of Beidou |
| CN113640838A (en) * | 2021-08-19 | 2021-11-12 | 中国电力科学研究院有限公司 | Multi-system combined satellite navigation positioning time service device and method |
| CN116399351A (en) * | 2023-04-23 | 2023-07-07 | 燕山大学 | Vehicle position estimation method |
| CN116859415A (en) * | 2023-05-23 | 2023-10-10 | 北京时代民芯科技有限公司 | Quick, stable and high-precision multi-fault satellite identification and positioning method |
| CN117040594A (en) * | 2023-08-10 | 2023-11-10 | 武汉大学 | Internet remote sensing satellite real-time service system oriented to mobile terminal user |
Non-Patent Citations (1)
| Title |
|---|
| 刘海颖 等: "基于ERAIM的惯性辅助卫星导航系统完好性检测", 中国惯性技术学报, vol. 18, no. 06, 31 December 2010 (2010-12-31), pages 686 - 690 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117761745B (en) | 2024-04-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109471146B (en) | Self-adaptive fault-tolerant GPS/INS integrated navigation method based on LS-SVM | |
| CN109581453B (en) | GNSS segmentation smoothing filtering method based on cycle slip detection and repair | |
| CN103163533B (en) | Seamless fusion expression and correction method of global navigation satellite system (GNSS) global and regional ionospheric delay | |
| CN109212556B (en) | Satellite navigation message parameter calculation method based on database description | |
| CN110187375A (en) | A kind of method and device improving positioning accuracy based on SLAM positioning result | |
| CN109059911B (en) | Data fusion method of GNSS, INS and barometer | |
| CN103064092A (en) | Selection method of navigational satellite | |
| CN103197335A (en) | Method using improved regularization method to restrain difference global positioning system (DGPS) integer ambiguity ill-condition | |
| CN112083457B (en) | Neural network optimized IMM satellite positioning navigation method | |
| CN116626730B (en) | Marine region CORS enhanced PPP method considering NWP | |
| CN109799518A (en) | Being navigated using low rail enhances the broadcast message layout broadcasting method of satellite | |
| CN111156986A (en) | Spectrum red shift autonomous integrated navigation method based on robust adaptive UKF | |
| CN117761745B (en) | Information optimization method for few-star positioning based on historical information fusion | |
| CN114779301A (en) | Satellite navigation real-time precise single-point positioning method based on broadcast ephemeris | |
| CN117272812A (en) | Low latitude small area ionosphere model construction method | |
| CN111694021B (en) | GNSS environment model-based single-station landslide deformation monitoring and early warning method | |
| CN113093240A (en) | Positioning correction method, device and equipment and computer readable storage medium | |
| CN111679297A (en) | Noise point drift removal method for GPS positioning track | |
| CN116147622A (en) | Combined navigation system fusion positioning method based on graph optimization | |
| CN115561793A (en) | Real-time Beidou phase decimal deviation rapid estimation method based on parallel computation | |
| CN113916217B (en) | Starlight positioning method based on regional stratosphere atmospheric refraction model | |
| CN115524731A (en) | High-precision aircraft track playback method based on post-calculation | |
| CN112700007B (en) | Training method, forecasting method and device of ionosphere electronic content forecasting model | |
| CN102129066B (en) | Fast wide satellite SAR geocoding method | |
| CN111427076B (en) | Vehicle performance evaluation method and system based on dual-mode satellite information deep fusion |
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 |