Single-epoch GNSS (Global Navigation Satellite System) double-difference parameter conversion method for any main satellite

Abstract

The invention discloses a single-epoch GNSS (Global Navigation Satellite System) double-difference parameter conversion method for any main satellites, which relates to the technical field of satellite navigation and positioning and comprises the following steps of: constructing a main satellite conversion matrix, a double-difference observed quantity sequence conversion matrix and an inverse matrix thereof between any two main satellites based on a single-epoch single-frequency single-system GNSS double-difference mathematical model; the method comprises the following steps: establishing a single-epoch multi-frequency multi-system main satellite transformation matrix with inverse transformation invariance and orthogonality, obtaining the invariance and orthogonality of matrix inverse transformation, constructing a single-epoch multi-frequency multi-system main satellite transformation matrix with inverse transformation invariance and orthogonality, a double-difference observed quantity sequence transformation matrix with inverse transformation invariance and orthogonality and an inverse matrix of the double-difference observed quantity sequence transformation matrix by establishing a single-epoch multi-frequency multi-system double-difference mathematical model and analogy derivation, and establishing a single-epoch multi-frequency multi-system and constructing a single-epoch GNSS double-difference model parameter and estimation parameter conversion and inverse transformation formula between any two or two groups of main satellites, and performing single-epoch double-difference parameter conversion.

Description

Technical Field

[0001] The present invention relates to the field of satellite navigation and positioning technology, and in particular to a single-epoch GNSS double-difference parameter conversion method for any main satellite. Background Art

[0002] Single-epoch differential technology based on the Global Navigation Satellite System (GNSS) can provide all-weather, all-day, fully automatic real-time positioning services. It has broad application prospects in shared vehicle management, unmanned driving, smart city construction, geological disaster monitoring, and other aspects. The correct fixation of the single-epoch ambiguity vector is a prerequisite for using GNSS differential technology to provide precise, real-time positioning services for the above applications.

[0003] Due to the short observation time and lack of verification information of single-epoch GNSS positioning, the accuracy of its ambiguity vector estimation and the reliability of its fixation face severe challenges; although it is possible to use the common double-difference ambiguity vector of the epoch that has been correctly fixed before the current epoch to enhance the estimation accuracy and fixation reliability of the double-difference ambiguity vector corresponding to the current epoch, the problem of primary satellite transformation between observation epochs will lead to changes in the model parameters (such as double-difference ambiguity vector, double-difference pseudorange / carrier observation, etc.) and estimation parameters (such as double-difference floating point / fixed deambiguation vector, ambiguity precision dilution ADOP, etc.) of the single-epoch GNSS double-difference mathematical model. It is difficult to directly use the ambiguity vector correctly fixed in the previous epoch to enhance the estimation accuracy and fixation reliability of the ambiguity vector of the current epoch. It is necessary to convert the double-difference parameters (model parameters and estimation parameters) of two or two groups of arbitrary different primary satellites between the observation epochs. However, the existing methods are only for the conversion of double-difference model parameters and estimated parameters between the first observation satellite as the main satellite and other observation satellites as the main satellite. There is a lack of double-difference mathematical model parameter and estimated parameter conversion methods for any two or two groups of main satellites, making it difficult to realize the conversion of double-difference ambiguity vectors and other estimated parameters between any different main satellites in adjacent epochs. Summary of the Invention

[0004] The embodiment of the present invention provides a single-epoch GNSS double-difference parameter conversion method for any master satellite, which can solve the problem in the prior art that it is difficult to convert double-difference ambiguity vectors and other estimated parameters between any different master satellites in adjacent epochs.

[0005] An embodiment of the present invention provides a single-epoch GNSS double-difference parameter conversion method for any primary satellite, comprising the following steps: Constructed with i The single-epoch single-frequency single-system GNSS double-difference mathematical model of the observation satellite as the main satellite is obtained. iSingle-epoch single-frequency single-system GNSS double-difference estimation parameters of the main satellite; Constructed with i The observation satellite is converted to the main satellite j The master satellite conversion matrix and double-difference observation quantity sequence transformation matrix of the single-epoch single-frequency single-system GNSS double-difference mathematical model when the observation satellite is the master satellite; Construct the master satellite conversion matrix and the inverse matrix of the double-difference observation sequence transformation matrix of the single-epoch single-frequency single-system GNSS double-difference mathematical model, and obtain the invariance of the inverse transformation of the master satellite conversion matrix and the orthogonality of the inverse transformation of the double-difference observation sequence transformation matrix; Constructed with A single-epoch multi-frequency multi-system GNSS double-difference mathematical model with an observation satellite as the main satellite is provided. The single-epoch multi-frequency multi-system GNSS is composed of multiple frequencies of multiple GNSS systems. According to the independence of double-difference observations between multiple frequencies of the single-epoch multi-frequency multi-system GNSS, the invariance of the inverse transformation of the main satellite transformation matrix, and the orthogonality of the inverse transformation of the double-difference observation quantity sequence transformation matrix, the main satellite transformation matrix and the double-difference observation quantity sequence transformation matrix of the double-difference mathematical model of multiple frequencies in the single-epoch multi-frequency multi-system GNSS are respectively arranged diagonally in blocks to obtain the double-difference mathematical model to be formed by the first observation satellite. The observation satellite was converted to The main satellite conversion matrix and double-difference observation quantity sequence transformation matrix of the single-epoch multi-frequency multi-system GNSS double-difference mathematical model with the observation satellite as the main satellite, as well as the invariance of the inverse transformation of the main satellite conversion matrix and the orthogonality of the inverse transformation of the double-difference observation quantity sequence transformation matrix for the single-epoch multi-frequency multi-system GNSS; and Both represent multiple observation satellites; According to the invariance of the main satellite conversion matrix and the double-difference observation sequence transformation matrix of the single-epoch multi-frequency multi-system GNSS double-difference mathematical model, as well as the orthogonality of the inverse transformation of the main satellite conversion matrix and the double-difference observation sequence transformation matrix, a new method based on the first The observation satellite was converted to The conversion formulas and inverse transformation formulas of the model parameters and estimated parameters of the single-epoch multi-frequency multi-system GNSS double-difference mathematical model when the observation satellite is the main satellite are used to perform the conversion between the single-epoch GNSS double-difference model parameters and estimated parameters.

[0006] The embodiment of the present invention provides a single-epoch GNSS double-difference parameter conversion method for any primary satellite. Compared with the prior art, the method has the following advantages: The present invention constructs a master satellite conversion matrix and an observation quantity sequence transformation matrix of a single-epoch single-frequency single-system GNSS double-difference mathematical model between any two master satellites, and constructs the inverse matrix of the master satellite conversion matrix and the observation quantity sequence transformation matrix and obtains the inverse transformation invariance and orthogonality. Then, based on the construction process of the master satellite conversion matrix and the observation quantity sequence transformation matrix of the single-frequency single-epoch single-system GNSS double-difference mathematical model, a single-epoch multi-frequency multi-system GNSS double-difference mathematical model is constructed to obtain a single-epoch multi-frequency multi-system GNSS double-difference with inverse transformation invariance and orthogonality. The main satellite conversion matrix and observation sequence transformation matrix of the mathematical model and their inverse matrix are used to obtain the model parameter and estimation parameter conversion formulas and inverse transformation formulas for the single-epoch single-frequency single-system GNSS double-difference mathematical model and the single-epoch multi-frequency multi-system GNSS double-difference mathematical model, and perform double-difference parameter conversion, thereby realizing the conversion of double-difference ambiguity vectors and other estimated parameters between any two or two groups of different main satellites in adjacent epochs. This has important theoretical value and practical guiding significance for constructing a single-epoch GNSS parameter estimation enhancement model and improving the ambiguity fixation success rate and fixation reliability. BRIEF DESCRIPTION OF THE DRAWINGS

[0007] Figure 1 A schematic diagram of the main flow of a single-epoch GNSS double-difference parameter conversion method for any primary satellite provided by an embodiment of the present invention. DETAILED DESCRIPTION

[0008] To make the above-mentioned objects, features, and advantages of the present invention more clearly understood, specific embodiments of the present invention are described in detail below with reference to the accompanying drawings. The following description sets forth many specific details to facilitate a full understanding of the present invention. However, the present invention can be implemented in many other ways than those described herein, and those skilled in the art can make similar modifications without violating the scope of the present invention. Therefore, the present invention is not limited to the specific embodiments disclosed below.

[0009] See also Figure 1 The embodiment of the present invention provides a single-epoch GNSS double-difference parameter conversion method for any master satellite, comprising: constructing a master satellite conversion matrix of a single-epoch single-frequency single-system GNSS double-difference mathematical model between any two master satellites. Z ; Construct the double-difference observation sequence transformation matrix of the single-epoch single-frequency single-system GNSS double-difference mathematical model L ; Establish the above transformation matrix Z and the transformation matrix L The inverse matrix of the inverse transformation is analyzed and the invariance and orthogonality of its inverse transformation are analyzed; by analogy with the above-mentioned single-epoch single-frequency single-system transformation matrix and observation sequence transformation matrix, the main satellite transformation matrix of the single-epoch multi-frequency multi-system GNSS with inverse transformation invariance and orthogonality is derived. Zand double difference observation order transformation moment L and its inverse matrix; establish a single-epoch GNSS double-difference model parameter and estimated parameter conversion method for any two or two groups of main satellites.

[0010] Specifically: 1. Single epoch GNSS double-difference mathematical model and parameter estimation.

[0011] This part mainly includes the single-epoch double-difference mathematical model (function model and random model) based on the Global Navigation Satellite System (GNSS), double-difference parameter estimation, ambiguity precision index and positioning precision index, and ambiguity fixation reliability theory based on R-ratio test.

[0012] 1.1GNSS double-difference positioning model

[0013] If GNSS is used to observe at the same time m satellites, they can be composed of Single-epoch single-frequency single-system double-difference carrier observations and pseudorange observations , based on the observed and The single-epoch single-frequency single-system GNSS double-difference mathematical model (function model E and random models D ) can be expressed as: (1) (2) in, and Respectively represent The expectation and variance of and Respectively represent The expectation and variance of b and a are the double difference baseline vector and double difference ambiguity vector respectively; B for b The column-full rank double difference coefficient matrix of ; , is the carrier wavelength; is the cofactor matrix of double-difference observations; and are the standard deviations of the undifferenced pseudorange and carrier observations, respectively. and its inverse matrix are all symmetric positive definite matrices. Using the least squares method, according to formula (1) and formula (2), we can get and Floating-point solution of and and its variance-covariance matrix and : (3) (4) in: and They are all symmetric positive definite matrices, and the least squares ambiguity reduction correlation adjustment method can be obtained and Fixed solution and : (5) in: for and The mutual factor matrix of .

[0014] 1.2 Single epoch ADOP and PDOP theory.

[0015] The Dilution of Precision (DOP) is a scalar factor that is easy to calculate. The DOPs involved in this invention include the Ambiguity DOP (ADOP) and the Positioning DOP (PDOP). The ADOP is primarily used to measure the accuracy of ambiguity estimation and the success rate of ambiguity fixation. It is defined as follows: (6) in: for The determinant of m - t is the fuzziness vector dimension, t The number of main satellites. For single-frequency single-system GNSS, t = 1. PDOP is mainly used to measure GNSS positioning accuracy and is defined as follows: (7) Where: trace( )for Generally speaking, the smaller the PDOP is, the higher the positioning accuracy is.

[0016] 1.3 Fuzziness Test Based on R-ratio Ambiguity checking is an important step in high-precision positioning using GNSS, which determines whether to accept the optimal integer solution of the ambiguity vector. is a fixed solution. Fuzzy test methods mainly include the success rate / failure rate index method, the statistical difference test method, and a combination of the two methods. The most commonly used one is the R-ratio test method, which belongs to the statistical difference test method and is defined as follows: (8) in: is the suboptimal integer solution for the ambiguity vector, is the threshold. In general, The values are empirical values such as 1.5, 2.0, 2.5, 3.0, etc.

[0017] 2. Single-epoch single-frequency single-system GNSS double-difference conversion matrix Z and the transformation matrix L .

[0018] According to formula (1) and formula (2), i Single-epoch single-frequency single-system GNSS double-difference function model with two observation satellites as main satellites E and random models D They can be written as: (9) (10) in: and Respectively The expectation and variance of and Respectively The expectation and variance of 、 、 、 、 and Respectively i Parameters of the single-epoch single-frequency single-system GNSS double-difference model with the observation satellite as the main satellite 、 、 B 、 Q 、 a and b ; b and a denote the double difference baseline vector and double difference ambiguity vector respectively; B express b The column-full rank double difference coefficient matrix of ; ; Indicates the carrier wavelength; The cofactor matrix representing the double-difference observations; and Represent the standard deviation of the undifferenced pseudorange and carrier observations respectively; single epoch single frequency single system GNSS double difference model parameters 、 、 、 and They are: (11) in: 、 and ( l =1, 2,…, m ) are respectively l The single-difference pseudorange observations, single-difference carrier observations and single-difference ambiguities between the receivers of the observation satellites; For the first i Single-epoch single-frequency single-system GNSS differential coefficient matrix of the observation satellite as the main satellite , for The identity matrix of order, The elements are all 1 dimensional column vector, For all elements to be 0 rank matrix, , , in particular, when i =1, ,when i = m hour, ; G is the non-difference unit direction vector between the receiver and the satellite. According to equations (3)-(5) and (9)-(11), i The single-epoch single-frequency single-system GNSS double-difference parameter estimation result of the observation satellite as the main satellite can be written as: (12) in: 、 、 、 、 、 、 、 and For the first i Single-epoch single-frequency single-system GNSS double-difference estimation parameters of the main satellite 、 、 、 、 、 、 、 and , for and In addition, according to formula (6)-formula (8), the first i The ADOP, PDOP and R-ratio of the single-epoch single-frequency single-system GNSS with two observation satellites as the main satellite can be written as 、 and : (13)

[0019] 2.1 Single-epoch single-frequency single-system GNSS primary satellite conversion matrix Z .

[0020] When the i The observation satellite is converted to the main satellite j When the single-epoch single-frequency single-system GNSS double-difference mathematical model of the main satellite is used ( ), due to the replacement of the main satellite, all model parameters in Equations (9) to (11) need to be converted to the main satellite, which involves the main satellite conversion matrix of the single epoch single frequency single system GNSS double difference mathematical model , The construction of mainly faces two situations, namely and ; can be expressed as: ① When The main satellite conversion matrix of the single-epoch single-frequency single-system GNSS double-difference mathematical model is Expressed as: (14)

[0021] in: Representation matrix of( x , y )element;" " indicates single epoch single frequency single system GNSS with the i The observation satellite is converted to the main satellite j The symbol of the observation satellite is the main satellite.

[0022] ②When The main satellite conversion matrix of the single-epoch single-frequency single-system GNSS double-difference mathematical model is Expressed as: (15) To facilitate the construction of the matrix , you can make The initial matrix , and then use formula (14) and formula (15) to assign values to the corresponding elements.

[0023] 2.2 Single-epoch single-frequency single-system GNSS double-difference observation sequence transformation matrix L .

[0024] ① When When , the matrix defined by formula (14) The order of double-difference observations in the converted single-epoch single-frequency single-system GNSS double-difference mathematical model will change, and it needs to be transformed. j The order of the double-difference observations is adjusted to i -1, original j +1 to i -1 double-difference observation quantity order minus 1, the corresponding double-difference observation quantity order transformation matrix of the single-epoch single-frequency single-system GNSS double-difference mathematical model It can be written as: (16) in: Representation matrix of( x , y ) elements; in particular, when hour, , that is, no need to The double-difference observations of the converted mathematical model are transformed sequentially.

[0025] ②When When , the matrix defined by formula (15) The order of the converted single-epoch single-frequency single-system GNSS double-difference observations also needs to be changed. j- 1 double difference observation is adjusted to the first i , the original i to j -The order of the two double-difference observations is increased by 1, and the corresponding double-difference observation order transformation matrix of the single-epoch single-frequency single-system GNSS double-difference mathematical model is It can be written as: (17) In particular, when hour, , that is, no need to The double-difference observations of the converted mathematical model are transformed sequentially.

[0026] To facilitate the construction of the matrix , you can make The initial matrix , and then use formula (16) and formula (17) to assign values to the corresponding elements.

[0027] 2.3 Single-epoch single-frequency single-system GNSS conversion matrix Z With the transformation matrix L Properties and its inverse matrix.

[0028] The master satellite conversion matrix of the single-epoch single-frequency single-system GNSS double-difference mathematical model described above and the double-difference observation order transformation matrix It has the following three important properties: ①Matrix Z and L The integer property of and From the definition of , we can see that and The elements in are integers such as 0, 1 or -1, that is, Z conversion and L The transformed ambiguity vector still has integer characteristics.

[0029] ②Matrix Z and L The reversibility and integer characteristics of the inverse matrix elements; According to the properties of elementary row / column transformations of determinants, 、 and Established, that is and are all reversible matrices; according to the adjoint matrix and the matrix determinant Inverse matrix Formula It can be seen that and The elements in the inverse matrix are also integers. This property ensures that the matrix Z and matrix L The inverse transformed ambiguity vector still has integer characteristics.

[0030] ③Matrix Z and L Invariance and orthogonality of inverse matrices; with the help of reversible matrices A The relationship , Established, that is It has inverse transformation invariance; by block reversible matrix Perform elementary row transformations, Established, that is Has orthogonal invariance. The matrix and The above properties ensure that Z and L Transformation or Z and LThe inverse transformed ambiguity vector still maintains the integer property.

[0031] 3. Single epoch multi-frequency multi-system GNSS double-difference conversion matrix Z and the transformation matrix L .

[0032] Assume that the same frequency or different frequencies of different systems select the main satellite in each system's respective frequency to form a double difference mathematical model and the multi-frequency multi-system GNSS includes d frequency, then d The double difference observations of the frequencies are independent or irrelevant to each other and are based on the first Double-difference function model of single-epoch multi-frequency multi-system GNSS with two observation satellites as main satellites E and random models D They can be written as: (18) (19) in: 、 、 、 and Respectively Model parameters of single-epoch multi-frequency multi-system GNSS with observation satellites as main satellites 、 、 、 and , ( k =1, 2,…, d ) is the k Frequency i observation satellites; , and ( k =1,2,…, d ) are the first k The wavelength of the frequency and the dimension of the double-difference ambiguity vector; The single-epoch multi-frequency multi-system GNSS differential coefficient matrix and double-difference cofactor matrix of the observation satellite as the main satellite are and ; 、 、 、 、 and ( k =1, 2,…, d ) are respectively k Frequencyi Double difference model parameters of the observation satellite as the main satellite 、 、 B 、 a 、 and Q , can be constructed with the help of the definition of relevant parameters in formula (11); , According to equations (18) and (19), the least square method and the least square ambiguity reduction correlation adjustment method are used to obtain Double-difference estimation parameters of single-epoch multi-frequency multi-system GNSS with two observation satellites as main satellites 、 、 、 、 、 、 、 、 and , using the model parameters and estimated parameters of the single-epoch multi-frequency multi-system GNSS double-difference mathematical model in equations (18) and (19) to replace the model parameters and estimated parameters of the single-epoch single-frequency single-system GNSS double-difference mathematical model in equation (12), we can obtain The specific analytical expressions of the single-epoch multi-frequency multi-system GNSS double-difference estimation parameters with the observation satellite as the main satellite are given.

[0033] With the help of the single epoch multi-frequency multi-system GNSS double difference mathematical model d The independence of the double-difference observations of the frequencies is obtained by diagonally arranging the main satellite conversion matrix and the double-difference observation sequence transformation matrix of multiple single-frequency single-system double-difference mathematical models in the multi-frequency multi-system GNSS. The observation satellite is converted to the main satellite The double-difference mathematical model of single-epoch multi-frequency multi-system GNSS when the observation satellite is the main satellite and the main satellite conversion matrix and the double-difference observation order transformation matrix , respectively expressed as: (20) (twenty one) in:" " indicates that the single epoch multi-frequency multi-system GNSS is The observation satellite is converted to the main satellite The symbol of the observation satellite as the main satellite; , ( k =1, 2,…, d ) is thek Frequency j observation satellites, ; s , For the k The number of satellites observing at the same frequency, and ( k =1, 2,…, d ) are respectively k The master satellite conversion matrix of the frequency and the double difference observation order transformation matrix are used to convert the first The double difference mathematical model of the observation satellite as the main satellite is converted into the The double difference mathematical model of the observation satellite as the main satellite can be based on and The numerical relationship between them is constructed using formula (14)-formula (17). In particular, when hour, .

[0034] Factors (20) and (21) are both diagonal matrices, and the matrix and Properties ①-③ are also applicable to the master satellite conversion matrix of the single epoch multi-frequency multi-system GNSS double-difference mathematical model and the double-difference observation order transformation matrix ; Therefore, the matrix and The ambiguity vector after transformation or inverse transformation still has integer characteristics.

[0035] 4. Single-epoch GNSS double-difference parameter conversion method for any two or two groups of main satellites.

[0036] According to the master satellite conversion matrix of the single epoch multi-frequency multi-system GNSS double-difference mathematical model in Equations (20) and (21), and the double-difference observation order transformation matrix ,Depend on" The general formula for model parameter conversion of the single-epoch multi-frequency multi-system GNSS double-difference mathematical model can be written as: (twenty two) in: 、 、 and Respectively Parameters of the single-epoch multi-frequency multi-system GNSS double-difference model with an observation satellite as the main satellite 、 、 and , ( k =1, 2,…, d ) is the k Frequency j observation satellites; The single-epoch multi-frequency multi-system GNSS differential coefficient matrix and double-difference cofactor matrix of the observation satellite as the main satellite are and ; 、 、 、 、 and ( k =1, 2,…, d ) are respectively k Frequency j Double difference model parameters of the observation satellite as the main satellite 、 、 B 、 a 、 and Q ; and .

[0037] According to the equations (18) and (19), The single-epoch multi-frequency multi-system GNSS double-difference estimation parameters of the observation satellite as the main satellite are obtained by The general formula for parameter conversion of single-epoch multi-frequency multi-system GNSS double-difference estimation can be written as: (twenty three) (twenty four) in, 、 、 、 、 、 and They are respectively for single epoch multi-frequency multi-system GNSS Observation satellites are the main satellites 、 、 、 、 、 and ; Can be further simplified to At the same time, ADOP definition (6), PDOP definition (7) and R-ratio definition (8) can be obtained by " The general formula for the conversion of single-epoch multi-frequency multi-system GNSS double-difference estimation parameters ADOP, PDOP and R-ratio is: (25) in, 、 and For single epoch multi-frequency multi-system GNSS The ADOP, PDOP and R-ratio values of the observation satellites as the main satellites; In particular, when it is a single-epoch single-frequency single-system GNSS, according to equations (22) to (25), by “ The general formulas for the conversion of single-epoch single-frequency single-system GNSS double-difference model parameters and estimated parameters can be expressed as: (26) (27) Among them: superscript " j " means the j An observation satellite is the main satellite; .

[0038] According to the invariance of the inverse transformation of the main satellite conversion matrix and the orthogonality of the inverse transformation of the double difference observation sequence transformation matrix and equations (22) to (25), we can get “ The general formulas for the inverse transformation of model parameters and estimated parameters of the single-epoch multi-frequency multi-system GNSS double-difference mathematical model are respectively expressed as: (28) (29) in: ;" "express" "Inverse transformation; other single epoch multi-frequency multi-system GNSS double difference estimation parameters 、 、 、 , ADOP, PDOP and R-ratio and Remain unchanged.

[0039] When it is a single-epoch single-frequency single-system GNSS, according to equations (26)-(27), “ The general formulas for the inverse transformation of the parameters and the estimated parameters of the single-epoch single-frequency single-system GNSS double-difference model can be expressed as follows: (30) (31) in:" "express" "Inverse transformation of ;Other single-epoch single-frequency single-system GNSS double-difference estimation parameters 、 、 、 , ADOP, PDOP and R-ratio and remain unchanged; According to equations (22) to (31), the single epoch GNSS double difference estimation parameters are 、 、 、 , ADOP, PDOP and R-ratio are all independent of the double-difference master satellite, that is, their values are independent of the selection of the master satellite; single epoch GNSS double-difference estimation parameters 、 and and double-difference mathematical model parameters (double-difference pseudorange observations , double-difference carrier observation , double difference coefficient matrix , double difference ambiguity vector parameters , difference coefficient matrix and the cofactor matrix of double-difference observations Q ) have double difference elementary transformation consistency, that is, their values vary with the selection of the main satellite, but they can be transformed by the main satellite matrix Z and the double-difference observation order transformation matrix L to convert.

[0040] The above-described embodiments merely illustrate several implementations of the present invention, and while their descriptions are relatively specific and detailed, they should not be construed as limiting the scope of the patent. It should be noted that a person skilled in the art would be able to make numerous variations and improvements without departing from the spirit of the present invention, all of which fall within the scope of protection of the present invention. Therefore, the scope of protection of the patent for this invention shall be determined by the appended claims.

Claims

1. A single-epoch GNSS double-difference parameter conversion method for any primary satellite, characterized in that: The following steps are involved: Constructed with i A single-epoch single-frequency single-system GNSS double-difference mathematical model with the observation satellite as the main satellite is established, and its double-difference estimation parameters are obtained; Build from i The primary satellite switches to the j The main satellite conversion matrix of the single-epoch single-frequency single-system GNSS double-difference mathematical model Z and the double-difference observation order transformation matrix L ; Constructing the Matrix Z and matrix L The inverse matrix of , and obtain the matrix Z Invariance of inverse transformations and matrices L Orthogonality of inverse transforms; Constructed with The mathematical model of single epoch multi-frequency multi-system GNSS double difference with the observation satellite as the main satellite; Based on the independence of double difference observations between multiple frequencies of single epoch multi-frequency multi-system GNSS, and the matrix Z Invariance of inverse transformations and matrices L The orthogonality of the inverse transformation is the matrix of the double-difference mathematical model of multiple frequencies in single-epoch multi-frequency multi-system GNSS. Z and matrix L Arrange the diagonal blocks respectively and obtain The primary satellite of the group switches to the The master satellite conversion matrix of the single-epoch multi-frequency multi-system GNSS double-difference mathematical model of the master satellite group Z' and the double-difference observation order transformation matrix L' , and the matrix Z' Invariance of inverse transformations and matrices L' Orthogonality of inverse transforms; According to the main satellite conversion matrix Z' and the double-difference observation order transformation matrix L' , and the matrix Z' Invariance of inverse transformations and matrices L' The orthogonality of the inverse transform is constructed from the The primary satellite of the group switches to the The conversion formulas and inverse transformation formulas of the model parameters and estimated parameters of the single-epoch multi-frequency multi-system GNSS double-difference mathematical model of the main satellite of the group are used to convert the single-epoch GNSS double-difference model parameters and estimated parameters.

2. The single-epoch GNSS double-difference parameter conversion method for any master satellite according to claim 1, characterized in that: The construction is based on i The single-epoch single-frequency single-system GNSS double-difference mathematical model with the observation satellite as the main satellite is constructed, and its double-difference estimation parameters are obtained, including: When using a single GNSS system to observe at the same time m satellites, they can form The first i Single-epoch single-frequency single-system GNSS double-difference carrier observations of the main satellite and pseudorange observations ,Depend on and Double difference function model E and double difference random model D Composed of i The mathematical model of single-epoch single-frequency single-system GNSS double difference with an observation satellite as the main satellite is: ； ； in: and Respectively The expectation and variance of and Respectively The expectation and variance of 、 、 、 、 and Respectively i Parameters of the single-epoch single-frequency single-system GNSS double-difference model with the observation satellite as the main satellite 、 、 B 、 Q 、 a and b ; b and a denote the double difference baseline vector and double difference ambiguity vector respectively; B express b The column-full rank double difference coefficient matrix of ; ; Indicates the carrier wavelength; The cofactor matrix representing the double-difference observations; and Represent the standard deviation of the undifferenced pseudorange and carrier observations respectively; single epoch single frequency single system GNSS double difference model parameters 、 、 、 and They are: ； in: 、 and Respectively l The single-difference pseudorange observations, single-difference carrier observations and single-difference ambiguities between the receivers of the observation satellites are obtained. l =1, 2,…, m ; For the first i Single-epoch single-frequency single-system GNSS differential coefficient matrix of the observation satellite as the main satellite , for The identity matrix of order, The elements are all 1 dimensional column vector, For all elements to be 0 rank matrix, , , in particular, when i =1, the differential coefficient is ,when i = m hour, ; G is the undifference unit direction vector between receiver and satellite; According to the single epoch single frequency single system GNSS double difference function model E and double difference random model D The formula is obtained by using the least square method. i Single-epoch single-frequency single-system GNSS parameters of the observation satellite as the main satellite and Floating-point solution of and and its variance-covariance matrix and , respectively expressed as: ； ； in: and are all symmetric positive definite matrices, and the least squares ambiguity reduction correlation adjustment method is used to obtain the i Single-epoch single-frequency single-system GNSS parameters of the observation satellite as the main satellite and Fixed solution and and The variance-covariance matrix of , expressed as: ； in, for and The mutual correlation factor matrix is obtained according to the definition of parameters ADOP, PODP and R-ratio. i Single-epoch single-frequency single-system GNSS double-difference estimation parameters of the main satellite 、 and , respectively expressed as: 。 3. The single-epoch GNSS double-difference parameter conversion method for any master satellite according to claim 2, characterized in that: The construction starts from i The primary satellite switches to the j The main satellite conversion matrix of the single-epoch single-frequency single-system GNSS double-difference mathematical model Z and the double-difference observation order transformation matrix L ,include: When the i The observation satellite is converted to the main satellite j When the single epoch single frequency single system GNSS double difference mathematical model is used as the main satellite, all model parameters in the double difference mathematical model need to be converted to the main satellite due to the replacement of the main satellite, involving the main satellite conversion matrix of the single epoch single frequency single system GNSS double difference mathematical model , The construction of mainly includes two cases, namely and ; They are respectively expressed as: when The main satellite conversion matrix of the single-epoch single-frequency single-system GNSS double-difference mathematical model is Expressed as: ； in: Representation matrix of( x , y )element;" " indicates single epoch single frequency single system GNSS with the i The observation satellite is converted to the main satellite j The symbol of the observation satellite as the main satellite; Main satellite conversion matrix The order of double-difference observations in the converted single-epoch single-frequency single-system GNSS double-difference mathematical model will change, and it needs to be transformed. j The order of the double-difference observations is adjusted to i -1, original j +1 to i -1 double-difference observation quantity order minus 1, the corresponding double-difference observation quantity order transformation matrix of the single-epoch single-frequency single-system GNSS double-difference mathematical model Expressed as: ； in: Representation matrix of( x , y ) elements; in particular, when hour, , that is, no need to The double-difference observations of the converted mathematical model are transformed sequentially; when The main satellite conversion matrix of the single-epoch single-frequency single-system GNSS double-difference mathematical model is Expressed as: ； Main satellite conversion matrix The order of the converted single-epoch single-frequency single-system GNSS double-difference observations needs to be changed. j- 1 double difference observation is adjusted to the first i , the original i to j -The order of the two double-difference observations is increased by 1, and the corresponding double-difference observation order transformation matrix of the single-epoch single-frequency single-system GNSS double-difference mathematical model is Expressed as: ； when hour, , no need to The double difference observations of the converted mathematical model are transformed in sequence; in order to facilitate the construction of the main satellite conversion matrix ,make The initial matrix , and then use the above The two formulas assign values to the corresponding elements; to construct the sequential transformation matrix ,make The initial matrix , using the above The two formulas assign values to the corresponding elements.

4. The single-epoch GNSS double-difference parameter conversion method for any master satellite according to claim 3, characterized in that: The obtained matrix Z Invariance of inverse transformations and matrices L Orthogonality of inverse transforms, including: Depend on and The definition of and The elements in are integers of 0, 1 or -1. According to the properties of elementary row / column transformation of determinant, and The determinants are expressed as: ； but and are all reversible matrices; according to the adjoint matrix and the matrix determinant Inverse matrix Formula , and The elements in the inverse matrix are also integers, and The inverse matrix of is expressed as: ； By partitioning the block reversible matrix Perform elementary row transformations, The inverse matrix of is expressed as: ； Then the main satellite conversion matrix of the single epoch single frequency single system GNSS double difference mathematical model is The inverse transformation is invariant and the double difference observation order transformation matrix The inverse transform of has orthogonality.

5. The single-epoch GNSS double-difference parameter conversion method for any master satellite according to claim 1, characterized in that: The obtained The primary satellite of the group switches to the The master satellite conversion matrix of the single-epoch multi-frequency multi-system GNSS double-difference mathematical model of the master satellite group Z' and the double-difference observation order transformation matrix L' , and the matrix Z' Invariance of inverse transformations and matrices L' Orthogonality of inverse transforms, including: The same frequency or different frequencies of different GNSS systems all select the main satellite within their respective frequencies to form a double difference mathematical model and the multi-frequency multi-system GNSS includes d frequency, then d The double difference observations of the frequencies are independent or irrelevant to each other and are based on the first Double-difference function model of single-epoch multi-frequency multi-system GNSS with two observation satellites as main satellites E and random models D Can be expressed as: ； ； in: 、 、 、 and Respectively Model parameters of single-epoch multi-frequency multi-system GNSS with observation satellites as main satellites 、 、 、 and , ( k =1, 2,…, d ) is the k Frequency i observation satellites; , and ( k =1,2,…, d ) are the first k The wavelength of the frequency and the dimension of the double-difference ambiguity vector; The single-epoch multi-frequency multi-system GNSS differential coefficient matrix and double-difference cofactor matrix of the observation satellite as the main satellite are and ; 、 、 、 、 and ( k =1, 2,…, d ) are respectively k Frequency i Double difference model parameters of the observation satellite as the main satellite 、 、 B 、 a 、 and Q ; , According to the above two equations, the least square method and the least square ambiguity reduction correlation adjustment method are used to obtain the Double-difference estimation parameters of single-epoch multi-frequency multi-system GNSS with two observation satellites as main satellites 、 、 、 、 、 、 、 、 and ; According to the single epoch multi-frequency multi-system GNSS double difference mathematical model d The independence of the double-difference observations of the frequencies is obtained by diagonally arranging the main satellite conversion matrix and the double-difference observation sequence transformation matrix of multiple single-frequency single-system double-difference mathematical models in the multi-frequency multi-system GNSS. The observation satellite was converted to The double-difference mathematical model of single-epoch multi-frequency multi-system GNSS when the observation satellite is the main satellite and the main satellite conversion matrix and the double-difference observation order transformation matrix , respectively expressed as: ； ； in:" " indicates that the single epoch multi-frequency multi-system GNSS is The observation satellite was converted to The symbol of the observation satellite as the main satellite; , ( k =1, 2,…, d ) is the k Frequency j observation satellites, ; s , For the k The number of satellites observing at the same frequency, and ( k =1, 2,…, d ) are respectively k The master satellite conversion matrix of the frequency and the double difference observation order transformation matrix are used to convert the first The double difference mathematical model of the observation satellite as the main satellite is converted into the The double difference mathematical model of the observation satellite as the main satellite is based on and The numerical relationship between them is used and Definitional construction and ,when hour, ; According to the main satellite conversion matrix of the single-epoch single-frequency single-system GNSS double-difference mathematical model Invariance of Inverse Transformation and Sequential Transformation Matrix of Double-Difference Observations Orthogonality of inverse transformation, master satellite conversion matrix of single epoch multi-frequency multi-system GNSS double-difference mathematical model The inverse transformation is invariant and the double difference observation order transformation matrix The inverse transform of has orthogonality.

6. The single-epoch GNSS double-difference parameter conversion method for any primary satellite according to claim 5, characterized in that: The construction starts from The primary satellite of the group switches to the The conversion formulas for the model parameters and estimated parameters of the single-epoch multi-frequency multi-system GNSS double-difference mathematical model of the primary satellite of the group include: According to the main satellite conversion matrix and the double-difference observation order transformation matrix , then by " The model parameter conversion formula of the single-epoch multi-frequency multi-system GNSS double-difference mathematical model is expressed as: ； in: 、 、 and Respectively Parameters of the single-epoch multi-frequency multi-system GNSS double-difference model with an observation satellite as the main satellite 、 、 and , ( k =1, 2,…, d ) is the k Frequency j observation satellites; The single-epoch multi-frequency multi-system GNSS differential coefficient matrix and double-difference cofactor matrix of the observation satellite as the main satellite are and ; 、 、 、 、 and ( k =1, 2,…, d ) are respectively k Frequency j Double difference model parameters of the observation satellite as the main satellite 、 、 B 、 a 、 and Q ; and ; Depend on" The conversion formula of the single-epoch multi-frequency multi-system GNSS double-difference estimation parameters is expressed as: ； ； in, 、 、 、 、 、 and They are respectively for single epoch multi-frequency multi-system GNSS Observation satellites are the main satellites 、 、 、 、 、 and ; Simplified to ; According to the definition of ADOP, PDOP and R-ratio, The conversion formula of single-epoch multi-frequency multi-system GNSS double-difference estimation parameters ADOP, PDOP and R-ratio is expressed as follows: ； in: 、 and They are respectively for single epoch multi-frequency multi-system GNSS The ADOP, PDOP and R-ratio values of the observation satellite as the main satellite; when it is a single-epoch single-frequency single-system GNSS, it is represented by" The conversion formulas of the double difference model parameters and estimated parameters are expressed as follows: ； ； in:" j " means the j An observation satellite is the main satellite; .

7. The single-epoch GNSS double-difference parameter conversion method for any master satellite according to claim 6, characterized in that: The construction starts from The primary satellite of the group switches to the The inverse transformation formulas for the model parameters and estimated parameters of the single-epoch multi-frequency multi-system GNSS double-difference mathematical model of the primary satellite of the group include: According to the invariance of the inverse transformation of the main satellite conversion matrix and the orthogonality of the inverse transformation of the double-difference observation sequence transformation matrix, " The formulas for inverse transformation of model parameters and inverse transformation of estimated parameters of the single-epoch multi-frequency multi-system GNSS double-difference mathematical model are respectively expressed as: ； ； in: ; "express" "Inverse transformation; other single epoch multi-frequency multi-system GNSS double difference estimation parameters 、 、 、 , ADOP, PDOP and R-ratio and remain unchanged; When it is a single-epoch single-frequency single-system GNSS, the The formulas for inverse transformation of parameters and estimated parameters of the single-epoch single-frequency single-system GNSS double-difference model are respectively expressed as: ； ； in:" "express" "Inverse transformation of ;Other single-epoch single-frequency single-system GNSS double-difference estimation parameters 、 、 、 , ADOP, PDOP and R-ratio and Remain unchanged.

8. The single-epoch GNSS double-difference parameter conversion method for any primary satellite according to claim 7, characterized in that: When converting the single epoch GNSS double difference model parameters and the estimated parameters, the single epoch GNSS double difference estimated parameters 、 、 、 , ADOP, PDOP and R-ratio all have double-difference main satellite independence, which means that the double-difference parameter values are independent of the selection of the main satellite; Single epoch GNSS double-difference estimation parameters 、 and And double-difference mathematical model parameters: double-difference pseudorange observation , double-difference carrier observation , double difference coefficient matrix , double difference ambiguity vector parameters , difference coefficient matrix and the cofactor matrix of double-difference observations Q They all have double difference elementary transformation consistency, which means that the double difference parameter values vary with the selection of the main satellite, and are converted by the main satellite transformation matrix Z and the double-difference observation order transformation matrix L to convert.

Images