CN110907960B - Cycle slip detection method and device based on K-Means dynamic clustering analysis - Google Patents
Cycle slip detection method and device based on K-Means dynamic clustering analysis Download PDFInfo
- Publication number
- CN110907960B CN110907960B CN201811085727.9A CN201811085727A CN110907960B CN 110907960 B CN110907960 B CN 110907960B CN 201811085727 A CN201811085727 A CN 201811085727A CN 110907960 B CN110907960 B CN 110907960B
- Authority
- CN
- China
- Prior art keywords
- chi
- cycle slip
- square
- epoch
- matrix
- 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
- 238000001514 detection method Methods 0.000 title claims abstract description 32
- 238000004458 analytical method Methods 0.000 title claims abstract description 28
- 239000011159 matrix material Substances 0.000 claims abstract description 49
- 238000000546 chi-square test Methods 0.000 claims abstract description 40
- 238000000354 decomposition reaction Methods 0.000 claims abstract description 17
- 238000007689 inspection Methods 0.000 claims abstract description 16
- 238000012545 processing Methods 0.000 claims abstract description 7
- 238000012360 testing method Methods 0.000 claims description 25
- 230000007123 defense Effects 0.000 claims description 7
- 238000007621 cluster analysis Methods 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 4
- 238000003556 assay Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 3
- 238000011524 similarity measure Methods 0.000 claims description 3
- 238000012795 verification Methods 0.000 claims 1
- 230000001934 delay Effects 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 238000000034 method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 2
- 241001061260 Emmelichthys struhsakeri Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003068 static effect Effects 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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/24—Acquisition or tracking or demodulation of signals transmitted by the system
- G01S19/27—Acquisition or tracking or demodulation of signals transmitted by the system creating, predicting or correcting ephemeris or almanac data within the receiver
-
- 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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/20—Integrity monitoring, fault detection or fault isolation of space segment
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Computer Security & Cryptography (AREA)
- Power Engineering (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention provides a cycle slip detection method and a cycle slip detection device based on K-Means dynamic clustering analysis, wherein the cycle slip detection method comprises the following steps: forming an inter-epoch difference observation equation by the current epoch data and the previous epoch data, carrying out QR decomposition on the obtained coefficient matrix, and constructing chi-square inspection volume; checking the chi-square checking quantity, if the chi-square checking quantity passes, all the observed values of the current epoch do not have cycle slip, and then processing the observed value of the next epoch; if the chi-square test does not pass, identifying the observed value of the cycle slip by a K-Means dynamic clustering analysis method, removing the observed value and reconstructing chi-square test quantity; and circularly executing the steps until the chi-square passes the inspection. The invention can effectively solve the problem that the threshold value is difficult to determine in cycle slip detection.
Description
Technical Field
The invention relates to the technical field of cycle slip detection, in particular to a cycle slip detection method and device based on K-Means dynamic clustering analysis.
Background
In the existing single-station single-frequency cycle slip detection and restoration method based on inter-epoch difference, a reference station of a previous epoch is assumed, a current epoch is a rover, observed values of cycle slip are identified according to residual errors after testing of the obtained observed values by combining an RTK positioning principle with robust least square estimation, and when at least 4 cycles of cycle slip does not occur, the cycle slip detection success rate is high.
And performing QR decomposition on the coefficient matrix, wherein the observed quantity with gross errors has strong correlation with the residual error after the experiment, namely the correlation distance is minimum, and the observed value with gross errors can be identified through common cluster analysis.
Cycle slip detection is a difficult problem that must be solved in GNSS (Global Navigation Satellite System) high-precision positioning technology. To acquire high-precision position information, a high-precision carrier phase observation value is required to be used, cycle slip inevitably exists in the carrier phase observation value, particularly, the cycle slip is frequent in complex environments such as urban canyons, and the finding of an effective cycle slip detection method is the key for realizing high-precision positioning. In the current cycle slip detection method, a cycle slip detection quantity is constructed, and a certain function model is selected to process the cycle slip detection quantity, so that a proper threshold value is selected to judge whether cycle slip occurs or not, wherein the threshold value can be constant or time-varying. For a static station with a wide observation environment, a better effect can be achieved, however, urban environments are complex and various, shielding is serious, pseudo-range and multi-path are frequent, and the effect of cycle slip detection by observation values similar to MW combination can be influenced. For the complex environment, the threshold value of cycle slip detection is often difficult to determine, and if the selected threshold value is too small, the 'normal observation value' is easily misjudged as the observation value with cycle slip; if the selected threshold is too large, the cycle slip may be missed for 1-2 weeks. Both of the above situations occur, which affect the performance of the positioning.
Disclosure of Invention
The invention provides a cycle slip detection method and device based on K-Means dynamic clustering analysis, which improve the GNSS positioning performance in a complex environment and solve the technical problems.
The technical scheme adopted by the invention is as follows:
a cycle slip detection method based on K-Means dynamic clustering analysis comprises the following steps:
forming an inter-epoch difference observation equation by the current epoch data and the previous epoch data, carrying out QR decomposition on the obtained coefficient matrix, and constructing chi-square inspection volume;
carrying out chi-square test on chi-square test quantity, if chi-square test is passed, all observed values of the current epoch have no cycle slip, and then processing the observed value of the next epoch; if the chi-square test does not pass, identifying the observed value of the cycle slip by a K-Means dynamic clustering analysis method, removing the observed value and reconstructing chi-square test quantity;
and circularly executing the steps until the chi-square passes the test.
Further, the current epoch data and the previous epoch data constitute an inter-epoch difference observation equation according to the following formula:
wherein f, i, s respectively represent frequency, epoch and satellite number; lambda [ alpha ]fRepresents a carrier wavelength; Represents the carrier phase observed value, with week as unit;representing a defense distance; t is a unit ofi sRepresenting tropospheric delay;indicating ionospheric delay; t is tiAnd ts,iRespectively a receiver clock error and a satellite clock error; biAnd bs,iCarrier phase hardware delays at the receiver end and the satellite end respectively;for blurring the whole weekDegree; epsiloniMeasurement noise that is a carrier phase observation; Δ represents the single difference operator between epochs;is the cycle slip value in weeks.
Further, the differential observation equation is linearized and written in matrix form:
wherein H is a coefficient matrix; δ x is a parameter to be estimated, including a relative position parameter and a receiver clock drift parameter; l is the pre-test residual error; v is the post-test residual.
Further, carrying out QR decomposition on the H coefficient matrix to obtain:
H=QR。
wherein Q is an orthogonal matrix; r is an upper triangular matrix with positive diagonal elements.
Further, the chi-square test quantity is equal to TLTTL, chi square assay compliance2(n-m) distribution, wherein m is the number of parameters to be estimated, n is the number of observed values, and T is taken as a matrix QT(QTThe transpose of the Q matrix), TL ═ T × L.
Further, identifying the observed value of the cycle slip by a K-Means dynamic clustering analysis method specifically comprises the following steps: constructing a data object; data objects are divided into two classes: the observed value of the cycle slip which does not occur and the observed value of the cycle slip which occurs; an observed value for occurrence of cycle slip is identified.
Further, the data object is constructed by the following formula:
where w is n-m, w represents the number of redundant observations, and each column in the formula represents a data object, i.e., d1、d2Each is expressed asThe similarity measure of the same data object is Euclidean distance, and d is selectedn+1And dn+1The most distant Euclidean vector is the initial clustering center and the vector dn+1The same category is the observed value of the occurrence of cycle slip.
The invention also provides a cycle slip detection device based on K-Means dynamic clustering analysis, which comprises:
the chi-square inspection volume construction unit is used for forming an inter-epoch difference observation equation based on the current epoch data and the previous epoch data, carrying out QR decomposition on the obtained coefficient matrix and constructing chi-square inspection volume;
the chi-square testing unit is used for carrying out chi-square testing on chi-square testing quantity, if the chi-square testing quantity passes, all the observed values of the current epoch do not have cycle slip, and then the observed value of the next epoch is processed; and if the chi-square test does not pass, identifying the observed value of the cycle slip by a K-Means dynamic clustering analysis method, removing the observed value and reconstructing chi-square test quantity.
The invention also provides a memory, in which a computer program is stored, the computer program performing the steps of:
Forming an inter-epoch differential observation equation by the current epoch data and the previous epoch data, carrying out QR decomposition on the obtained coefficient matrix, and constructing chi-square inspection variables;
carrying out chi-square test on chi-square test quantity, if chi-square test is passed, all the observed values of the current epoch do not have cycle slip, and then processing the observed value of the next epoch; if the chi-square test does not pass, identifying the observed value of the cycle slip by a K-Means dynamic clustering analysis method, removing the observed value and reconstructing chi-square test quantity;
and circularly and repeatedly executing the steps until the chi-square passes the test.
The invention has the following beneficial effects: the observed value of the cycle slip is identified by applying K-Means dynamic clustering analysis and chi-square test in mathematical statistics, and the problem that the threshold value is difficult to determine in cycle slip detection can be effectively solved.
Drawings
FIG. 1 is a cycle slip detection flow chart according to the present invention.
Detailed Description
The invention solves the technical problem that the threshold value is difficult to determine in cycle slip detection by applying chi-square test and K-Means dynamic clustering analysis. The invention is further illustrated below with reference to the figures and examples.
The first embodiment is as follows:
the invention provides a cycle slip detection method based on K-Means dynamic clustering analysis, a flow chart is shown in figure 1, and the cycle slip detection method comprises the following steps:
Step 1, a difference observation equation between epochs can be formed by current epoch data and previous epoch data according to a formula (2), QR decomposition (orthogonal triangle decomposition, Q is an orthogonal matrix, and R is an upper triangle matrix with positive diagonal elements) is carried out on an obtained coefficient array, and then a chi-square inspection quantity Test is constructed according to a formula (5).
Wherein f, i, s respectively represent frequency, epoch and satellite number; lambda [ alpha ]fRepresents a carrier wavelength;represents a carrier phase observation in units of weeks;representing a defense distance; t isi sRepresenting tropospheric delay;indicating ionospheric delay; t is tiAnd ts,iRespectively a receiver clock error and a satellite clock error; biAnd bs,iCarrier phase hardware delays at the receiver end and the satellite end respectively;is the integer ambiguity; epsiloniMeasurement noise that is a carrier phase observation.
The differential observation equation between epoch i +1 and epoch i can be obtained according to equation (1):
wherein,
Δ represents the single difference operator between epochs;is the cycle slip value in weeks. And equation (2) ignores the amount that can be considered constant in a short time.
Linearize equation (2) and write in matrix form:
wherein H is a coefficient matrix (design matrix); and deltax is a parameter to be estimated, including a relative position parameter and a receiver clock drift parameter, L is a residual before the experiment, and v is a residual after the experiment. The cycle slip parameter is omitted from equation (3) assuming that no cycle slip occurs for all observations.
Performing QR decomposition on the H matrix in equation (3) to obtain:
H=QR (4)
wherein Q is an orthogonal matrix; r is an upper triangular matrix with positive diagonal elements.
Assuming that the number of the parameters to be estimated is m, the number of the observed values is n, QTIs the transpose of the Q matrix. Get matrix QTThe lower half (denoted T) of (n-m) × n is multiplied by the upper matrix to obtain:
TL=T*L (5)
let Test be TLTTL de: test compliance chi2(n-m) distribution, and chi-square Test is performed on Test. If passing the testIf the result is verified, the observation data of the epoch is considered to have no cycle slip; if the test is not passed, the cycle slip of the epoch data is shown, and the cycle slip observed value is identified by the K-Means dynamic clustering analysis method in the step 2.
And 2, if the chi-square passes the checking, all the observed values of the epoch have no cycle slip, and then the observed data of the next epoch is processed. And if the chi-square test is not passed, identifying the observed value of the cycle slip by a K-Means dynamic clustering analysis method, and removing the observed value.
K-Means cluster analysis classifies a plurality of data objects, and the data objects can be classified into K classes. While for cycle slip detection the data can be categorized into two categories: one is an observed value in which cycle slip does not occur; the other is the observed value of cycle slip. Let w be n-m, w denote the number of redundant observations, and now construct the data object as follows:
Each column in equation (6) represents a data object, i.e., d1、d2Etc. respectively representing different data objects, and the similarity measure is Euclidean distance. Now set the value of k to 2, and divide equation (6) into two categories. Selection of dn+1And dn+1The vector with the farthest Euclidean distance is the initial clustering center. Classifying the data object by adopting a K-Means dynamic clustering analysis method and vector dn+1The same category is the observed value of the occurrence of cycle slip.
And 3, circularly and repeatedly executing the step 1 and the step 2 until the chi-square passes the inspection.
Often, only one observed value of cycle slip can be identified for one K-Means cluster analysis, the observed value of cycle slip is often the largest, and in order to solve the problem that multiple observed values simultaneously generate cycle slip, iteration processing needs to be performed on the step 1 and the step 2, that is: and eliminating cycle slip observed values identified by each K-Means dynamic analysis, and then carrying out chi-square test, and repeating the steps until the chi-square test is passed.
The second embodiment:
the invention also provides a cycle slip detection device based on K-Means dynamic clustering analysis, which comprises:
the chi-square inspection volume construction unit is used for forming an inter-epoch differential observation equation based on the current epoch data and the previous epoch data, carrying out QR decomposition on the obtained coefficient matrix and constructing chi-square inspection volume;
The chi-square testing unit is used for carrying out chi-square testing on chi-square testing quantity, if the chi-square testing quantity passes, all the observed values of the current epoch do not have cycle slip, and then the observed value of the next epoch is processed; and if the chi-square test does not pass, identifying the observed value of the cycle slip by a K-Means dynamic clustering analysis method, removing the observed value and reconstructing chi-square test quantity.
Further, the differential observation equation is as follows:
wherein f, i, s respectively represent frequency, epoch and satellite number; lambda [ alpha ]fRepresents a carrier wavelength;represents a carrier phase observation in units of weeks;representing a defense distance; t isi sRepresenting tropospheric delay;indicating ionospheric delay; t is tiAnd ts,iRespectively a receiver clock error and a satellite clock error; biAnd bs,iCarrier phase hardware delays at the receiver end and the satellite end respectively;is the integer ambiguity; epsiloniMeasurement noise that is a carrier phase observation; Δ represents the single difference operator between epochs;is the cycle slip value in weeks.
Example three:
the invention also provides a memory, in which a computer program is stored, the computer program performing the steps of:
forming an inter-epoch difference observation equation by the current epoch data and the previous epoch data, carrying out QR decomposition on the obtained coefficient matrix, and constructing chi-square inspection volume;
Checking the chi-square checking quantity, if the chi-square checking quantity passes, all the observed values of the current epoch do not have cycle slip, and then processing the observed value of the next epoch; if the chi-square test does not pass, identifying the observed value of the cycle slip by a K-Means dynamic clustering analysis method, removing the observed value and reconstructing chi-square test quantity;
and circularly executing the steps until the chi-square passes the test.
Although the present invention has been described with reference to the preferred embodiments, it is not intended to limit the present invention, and those skilled in the art can make variations and modifications of the present invention without departing from the spirit and scope of the present invention by using the methods and technical contents disclosed above.
Claims (6)
1. A cycle slip detection method based on K-Means dynamic clustering analysis is characterized by comprising the following steps:
the current epoch data and the previous epoch data form an inter-epoch differential observation equation, the differential observation equation is linearized and written into a matrix form:
Wherein H is a coefficient matrix; δ x is a parameter to be estimated, including a relative position parameter and a receiver clock drift parameter; l is the pre-test residual error; v is the post-test residual; lambda [ alpha ]fRepresents a carrier wavelength;representing a carrier phase observation;representing a defense distance;indicating ionospheric delay; t is ts,iIs the satellite clock error; f, i, s respectively represent frequency, epoch and satellite number; Δ represents the single difference operator between epochs;
carrying out QR decomposition on the obtained coefficient matrix, and constructing chi-square inspection quantity; and carrying out QR decomposition on the H coefficient matrix to obtain: where Q is an orthogonal matrix and R is an upper triangular matrix with positive diagonal elements;
carrying out chi-square test on chi-square test quantity, if chi-square test is passed, all observed values of the current epoch have no cycle slip, and then processing the observed value of the next epoch; if the chi-square test does not pass, identifying the observed value of the cycle slip by a K-Means dynamic clustering analysis method, removing the observed value and reconstructing chi-square test quantity;
circularly executing the steps until the chi-square passes the inspection;
the chi-square test quantity is equal to TLTTL, chi square assay compliance2(n-m) distribution, wherein m is the number of parameters to be estimated, n is the number of observed values, and T is taken as a matrix Q TLower half of (n-m) × n, QTTL is the transpose of the Q matrix.
2. The cycle slip detection method based on K-Means dynamic cluster analysis of claim 1, wherein the current epoch data and the previous epoch data form an inter-epoch differential observation equation according to the following formula:
wherein f, i, s respectively represent frequency, epoch and satellite number; lambda [ alpha ]fRepresents a carrier wavelength;represents a carrier phase observation in units of weeks;representing a defense distance;representing tropospheric delay;indicating ionospheric delay; t is tiAnd ts,iRespectively a receiver clock error and a satellite clock error; biAnd bs,iAre respectively connected toThe carrier phase hardware delay of the receiver end and the satellite end;is the integer ambiguity; epsiloniMeasurement noise that is a carrier phase observation; Δ represents the single difference operator between epochs;is the cycle slip value in weeks.
3. The cycle slip detection method based on K-Means dynamic cluster analysis as claimed in claim 1, wherein the step of identifying the observed value of the cycle slip by the K-Means dynamic cluster analysis specifically comprises the following steps:
constructing a data object;
data objects are divided into two classes: the observed value of the cycle slip which does not occur and the observed value of the cycle slip which occurs;
An observed value for occurrence of cycle slip is identified.
4. The cycle slip detection method based on K-Means dynamic clustering analysis as claimed in claim 1, wherein the data object is constructed by the following formula:
where w is n-m, w represents the number of redundant observations, and each column in the formula represents a data object, i.e., d1、d2The similarity measure is Euclidean distance, d is selectedn+1And dn+1The most distant Euclidean vector is the initial clustering center and the vector dn+1The same category is the observed value of the occurrence of cycle slip.
5. A cycle slip detection device based on K-Means dynamic clustering analysis is characterized by comprising:
the chi-square inspection volume construction unit is used for forming an inter-epoch difference observation equation based on the current epoch data and the previous epoch data, linearizing the difference observation equation and writing the difference observation equation into a matrix form:
wherein H is a coefficient matrix; δ x is a parameter to be estimated, including a relative position parameter and a receiver clock drift parameter; l is the pre-test residual error; v is the residual error after the test, λfRepresents a carrier wavelength;representing a carrier phase observation;representing a defense distance;indicating ionospheric delay; t is ts,iIs the satellite clock error; f, i, s respectively represent frequency, epoch and satellite number; Δ represents the single difference operator between epochs; carrying out QR decomposition on the obtained coefficient matrix, and constructing chi-square inspection quantity; carrying out QR decomposition on the H coefficient matrix to obtain: where Q is an orthogonal matrix and R is an upper triangular matrix with positive diagonal elements;
The chi-square checking unit is used for checking chi-square checking quantity, if the chi-square checking quantity passes, all the observed values of the current epoch do not have cycle slip, and then the observed value of the next epoch is processed; if the card party verification does not pass,
identifying the observed value of the cycle slip by a K-Means dynamic clustering analysis method, removing the observed value and reconstructing chi-square test quantity;
the chi-square test quantity is equal to TLTTL, chi square test compliance χ2(n-m) distribution, wherein m is the number of parameters to be estimated, n is the number of observed values, and T is a matrix QTLower half of (n-m) × n, QTTL is the transpose of Q matrix.
6. A memory storing a computer program, the computer program performing the steps of:
the current epoch data and the previous epoch data form an inter-epoch differential observation equation, the differential observation equation is linearized and written into a matrix form:
wherein H is a coefficient matrix; δ x is a parameter to be estimated, including a relative position parameter and a receiver clock drift parameter; l is the pre-test residual error; v is the residual error after the test, λfRepresents a carrier wavelength;representing a carrier phase observation;representing a defense distance; Representing ionospheric delay; t is ts,iIs the satellite clock error; f, i, s respectively represent frequency, epoch and satellite number; Δ represents the single difference operator between epochs; carrying out QR decomposition on the obtained coefficient matrix, and constructing chi-square inspection quantity; carrying out QR decomposition on the H coefficient matrix to obtain: where Q is an orthogonal matrix and R is an upper triangular matrix with positive diagonal elements;
carrying out chi-square test on chi-square test quantity, if chi-square test is passed, all the observed values of the current epoch do not have cycle slip, and then processing the observed value of the next epoch; if the chi-square test does not pass, identifying the observed value of the cycle slip by a K-Means dynamic clustering analysis method, removing the observed value and reconstructing chi-square test quantity;
circularly executing the steps until the chi-square is checked to pass;
the chi-square test quantity is equal to TLTTL, chi square assay compliance2(n-m) distribution, wherein m is the number of parameters to be estimated, n is the number of observed values, and T is taken as a matrix QTLower half of (n-m) × n, QTTL is the transpose of the Q matrix.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811085727.9A CN110907960B (en) | 2018-09-17 | 2018-09-17 | Cycle slip detection method and device based on K-Means dynamic clustering analysis |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811085727.9A CN110907960B (en) | 2018-09-17 | 2018-09-17 | Cycle slip detection method and device based on K-Means dynamic clustering analysis |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110907960A CN110907960A (en) | 2020-03-24 |
CN110907960B true CN110907960B (en) | 2022-06-28 |
Family
ID=69813669
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811085727.9A Active CN110907960B (en) | 2018-09-17 | 2018-09-17 | Cycle slip detection method and device based on K-Means dynamic clustering analysis |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110907960B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112305574A (en) * | 2020-05-21 | 2021-02-02 | 国网浙江省电力有限公司 | Beidou GNSS satellite real-time positioning and orientation data preprocessing system and method |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101433908B1 (en) * | 2013-09-24 | 2014-08-27 | 한국과학기술원 | Method and system for data quality check of gnss observation |
CN105068092A (en) * | 2015-07-17 | 2015-11-18 | 北京航空航天大学 | Cycle slip detection and restoration method applied to satellite based augmentation system airborne receiver |
CN107544082A (en) * | 2016-06-29 | 2018-01-05 | 武汉大学 | The step modeling of Big Dipper IGSO/MEO satellite pseudorange codes deviation one |
CN107728181A (en) * | 2017-10-10 | 2018-02-23 | 北京无线电计量测试研究所 | A kind of Detection of Cycle-slip restorative procedure in real time |
CN108363077A (en) * | 2017-12-29 | 2018-08-03 | 中国电子科技集团公司第二十研究所 | A kind of carrier phase cycle slip rehabilitation method in Static Precise Point Positioning device |
CN108519045A (en) * | 2018-05-14 | 2018-09-11 | 桂林电子科技大学 | A kind of Big Dipper precision deformation monitoring and early warning system |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10416306B2 (en) * | 2015-08-17 | 2019-09-17 | Texas Instruments Incorporated | Methods and apparatus to measure and analyze vibration signatures |
-
2018
- 2018-09-17 CN CN201811085727.9A patent/CN110907960B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101433908B1 (en) * | 2013-09-24 | 2014-08-27 | 한국과학기술원 | Method and system for data quality check of gnss observation |
CN105068092A (en) * | 2015-07-17 | 2015-11-18 | 北京航空航天大学 | Cycle slip detection and restoration method applied to satellite based augmentation system airborne receiver |
CN107544082A (en) * | 2016-06-29 | 2018-01-05 | 武汉大学 | The step modeling of Big Dipper IGSO/MEO satellite pseudorange codes deviation one |
CN107728181A (en) * | 2017-10-10 | 2018-02-23 | 北京无线电计量测试研究所 | A kind of Detection of Cycle-slip restorative procedure in real time |
CN108363077A (en) * | 2017-12-29 | 2018-08-03 | 中国电子科技集团公司第二十研究所 | A kind of carrier phase cycle slip rehabilitation method in Static Precise Point Positioning device |
CN108519045A (en) * | 2018-05-14 | 2018-09-11 | 桂林电子科技大学 | A kind of Big Dipper precision deformation monitoring and early warning system |
Non-Patent Citations (4)
Title |
---|
GNSS完备性监测方法、技术与应用;秘金钟;《中国博士学位论文全文数据库 基础科学版》;20101015(第10期);第103-104、110-111页 * |
INS辅助GPS模糊度快速恢复与周跳探测技术研究;谢宏飞;《中国优秀硕士学位论文全文数据库 基础科学辑》;20140215(第2期);第6-9、55-58页、图5.1 * |
K-Means聚类算法的研究;冯超;《中国优秀硕士学位论文全文数据库 信息科技辑》;20080515(第5期);第31-33页 * |
联合惯导的GPS周跳探测与修复算法研究;郑蓉;《中国优秀硕士学位论文全文数据库 基础科学辑》;20170715(第7期);第15-16页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110907960A (en) | 2020-03-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Arora et al. | Redescending estimators for data reconciliation and parameter estimation | |
Borzyszkowski et al. | liger: mock relativistic light cones from Newtonian simulations | |
Castaldo et al. | P‐RANSAC: An Integrity Monitoring Approach for GNSS Signal Degraded Scenario | |
US9213100B1 (en) | Bearing-only tracking for horizontal linear arrays with rapid, accurate initiation and a robust track accuracy threshold | |
Shapiro et al. | Will multiple probes of dark energy find modified gravity? | |
CN111308418B (en) | Steady method for two-dimensional positioning of target with unknown height | |
Carbone et al. | PySE: Software for extracting sources from radio images | |
CN110907960B (en) | Cycle slip detection method and device based on K-Means dynamic clustering analysis | |
Schroth et al. | Failure detection and exclusion via range consensus | |
Costantino et al. | Seismic source characterization from GNSS data using deep learning | |
Angrisano et al. | A resampling strategy based on bootstrap to reduce the effect of large blunders in GPS absolute positioning | |
CN111025354A (en) | Medium-long baseline RTK positioning method based on single-differential ionosphere weighting model | |
García-Alfonso et al. | Comparison of uncertainty analysis of the Montecarlo and Latin Hypercube algorithms in a camera calibration model | |
CN112083463A (en) | Method and device for detecting whether ambiguity is fixed correctly or not and positioning terminal | |
Wyszkowska et al. | Performance of Msplit estimates in the context of vertical displacement analysis | |
US20120050728A1 (en) | Identifying defects | |
CN106873003A (en) | A kind of dual mode receiver RAIM methods based on the accumulation of many epoch parity vectors | |
US20190260572A1 (en) | Efficient computation of bivariate statistical moments for side channel vulnerability evaluation | |
Grigoli et al. | Relative earthquake location procedure for clustered seismicity with a single station | |
Almagbile | Geometric and statistical interpretation of correlation between fault tests in integrated GPS/INS systems | |
US11515995B2 (en) | Efficient computation of univariate statistical moments for side channel vulnerability evaluation | |
Dilmaç et al. | New approaches for outlier detection: The least trimmed squares adjustment | |
Ni et al. | An improved RAIM scheme for processing multiple outliers in GNSS | |
Granat et al. | Statistical approaches to detecting transient signals in GPS: Results from the 2009–2011 transient detection exercise | |
Bekiroglu et al. | Hankel matrix rank as indicator of ghost in bearing-only tracking |
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 | ||
CB02 | Change of applicant information |
Address after: 200438 9 / F, 10 / F, 11 / F, 12 / F, 38 Lane 1688, Guoquan North Road, Yangpu District, Shanghai Applicant after: QIANXUN SPATIAL INTELLIGENCE Inc. Address before: Room j165, 1st floor, building 64, 1436 Jungong Road, Yangpu District, Shanghai, 200433 Applicant before: QIANXUN SPATIAL INTELLIGENCE Inc. |
|
CB02 | Change of applicant information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |