CN109490916A - A kind of GNSS receiver autonomous integrity monitoring method - Google Patents

A kind of GNSS receiver autonomous integrity monitoring method Download PDF

Info

Publication number
CN109490916A
CN109490916A CN201910052794.9A CN201910052794A CN109490916A CN 109490916 A CN109490916 A CN 109490916A CN 201910052794 A CN201910052794 A CN 201910052794A CN 109490916 A CN109490916 A CN 109490916A
Authority
CN
China
Prior art keywords
epoch
moment
variance
extrapolation
gnss receiver
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910052794.9A
Other languages
Chinese (zh)
Other versions
CN109490916B (en
Inventor
孙蕊
张文宇
何秀凤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Aeronautics and Astronautics
Hohai University HHU
Original Assignee
Nanjing University of Aeronautics and Astronautics
Hohai University HHU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing University of Aeronautics and Astronautics, Hohai University HHU filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201910052794.9A priority Critical patent/CN109490916B/en
Publication of CN109490916A publication Critical patent/CN109490916A/en
Application granted granted Critical
Publication of CN109490916B publication Critical patent/CN109490916B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
    • G01S19/01Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
    • G01S19/13Receivers
    • G01S19/23Testing, monitoring, correcting or calibrating of receiver elements

Abstract

The invention discloses a kind of GNSS receiver autonomous integrity monitoring methods, are related to receiver autonomous integrity monitoring technology, belong to the technical field of measurement test.This method uses the monitoring that GNSS receiver autonomous integrity is realized based on CKF and new breath extrapolation, the test statistics at moment of extrapolating each time is constructed by newly ceasing the new breath generated in extrapolation process and its variance-covariance battle array, small gradual pseudorange biases can be examined, dynamic updates the Mishap Database for generating machine learning rule while differentiated and classified using test statistics of the machine learning algorithm to each extrapolation moment, the judging rules of the unmodeled fault type of on-line study, it is less with calculation amount, precision is suitable, the advantages of suitable for nonlinear system and capable of monitoring small gradual pseudorange biases to high confidence level.

Description

A kind of GNSS receiver autonomous integrity monitoring method
Technical field
The invention discloses a kind of GNSS receiver autonomous integrity monitoring methods, are related to receiver autonomous integrity monitoring Technology belongs to the technical field of measurement test.
Background technique
To RAIM (Receiver Autonomous Integrity Monitoring, receiver autonomous integrity monitoring) Research start from last century the mid-80.In navigation association, the U.S. in 1986, Y.C.Lee proposes Agonists by Distance Comparison Method, That is, when making integrity decision every time only with the observed quantity of current epoch, according to the weighted sum of five satellite distance errors Failure is judged whether there is, this method, which is only used for detection failure, not can determine that fault satellites but.Next year, R.M.Kalafus is just Minimum satellite vehicle number needed for formula proposes the concept of RAIM and demonstrates fault-finding and Fault Isolation.1988, Bradford W.Parkinson proposes least square method by carrying out processing to residual error, and this method handles pseudorange by least-squares estimation After obtaining residual error, the probability density characteristics according to residual error carry out decision;The same year, Mark A.Sturzar propose parity vector Noise error is placed in the progress fault detection of odd even space by dexterously carrying out QR decomposition to observation battle array by method, verified, minimum Square law and parity vector method are equivalent in result.In view of the suggestiveness of parity vector method computationally, this method becomes The recommendation RAIM algorithm of intact sex service in GPS application.
The design pattern of Agonists by Distance Comparison Method, least square method and odd_even adjudgement rule can be included in instant method, also be " snapshot The scope of method " carries out the influence that not will receive the conception of history measured data when decision every time.Work as since the judgement to integrity is based only on The observation at preceding moment, the obvious fault occurred suddenly can be detected in traditional RAIM algorithm, in stationary random process Failure is then not sensitive enough.
In order to improve the sensitivity of RAIM algorithm, scholar is made that more effort.Wang Lingxuan et al. proposes general using priori Rate model is to fault modeling, by minimum risk cost arbitrary target Functional Analysis residual error statistic, to optimize threshold value to mention It is highly sensitive.Li Chao et al. then devises RAIM algorithm according to least-square residuals method, and using Beidou satellite navigation system into Verifying is gone.Ran Jianhua etc. establishes the RAIM algorithm based on polynary accumulation sum, which has carried out place to parity vector residual error Reason not only can detect biggish saltant type failure with a step, can also detect the accumulation type failure of moderate, but for larger Odd even residual error average drifting then lose the ability further detected.Wang Shitai proposes the RAIM algorithm based on M estimation, should The method that the robust characteristic that algorithm is estimated using M devises fault-finding and Fault Isolation, and pass through simulating, verifying algorithm Performance is substantially better than least square method.It is devised for the larger problem of Multipath Errors, tension et al. under city and canyon environment Adaptive RAIM algorithm tests to positioning accuracy before and after troubleshooting using precision filtered device is excluded, so that enhancing is fixed The accuracy of position.
With BDS (BeiDou Navigation Satellite System, Chinese Beidou satellite navigation system), The global position systems such as GLONASS (Global Navigation Satellite System, Global Navigation Satellite System) It emerges and opens, more constellation RAIM algorithms also result in the attention of some scientific research personnel.Chen Ting is in paper by single constellation With the performance evaluation of RAIM algorithm under more constellation combinations, it is indicated that the RAIM performance under more constellation combinations is significantly stronger than single constellation feelings Condition.Wang Haijun proposes more epoch accumulation double star failure methods, which increases and improve with epoch number, to micro- Poor, middle difference and rough error can be detected.Lu Deqian et al. proposes the RAIM algorithm dynamically distributed based on integrity index, leads to Crossing dynamic allocation omission factor reduces vertical protection class, to improve the availability of RAIM.Wang Ershen et al., which gives, to be based on BDS/GPS positioning weighted average solution algorithm, the algorithm to optimal weighting be averaged resolving Algorithm weighted factor optimize from And test statistics is designed by positioning solution, it is finally debugged using weighted least-squares method, which can detecte same Double star failure of the Shi Fasheng in various constellations.
Traditional RAIM algorithm is carried out after by system linearization, with the proposition and application of various algorithm, people Be also to develop the RAIM algorithm of suitable nonlinear system to make efforts.Song Jiancai et al. is proposed to be filtered based on spreading kalman The RAIM algorithm of wave, the algorithm carry out recursion using instant data and historical data simultaneously, can take into account time-domain information, in addition, The step of handling observed quantity using proper orthogonal decomposition in paper, eliminating finding the inverse matrix effectively reduces calculation amount, and by imitative Really demonstrate the ability of its fault detection.The RAIM that Wang Ershen et al. proposes the particle filter based on chaotic particle swarm optimization is calculated Method, the algorithm can detect and debug under conditions of coloured noise system, they also proposed based on hereditary grain The RAIM algorithm of son filtering improves particle filter accuracy come the problem coped with sample degeneracy with sample exhaustion to enhance RAIM calculation The detection performance of method.
In addition, various enhancing RAIM algorithms more attract attention.Shi Yibing gives the RAIM algorithm of clock auxiliary, the party Method predicts receiver clock deviation using discrete Grey System Model, and the receiver clock-offsets sequence of prediction is equivalent to One satellite, to promote the performance of RAIM in the case of few star.Sun Minghan et al. is also to three kinds of judgements of RAIM algorithm availability Method (δ Hmax, ARP and HPL method) carries out analysis comparison, and demonstrates the equivalence of δ Hmax and HPL method.
On the whole, according to the analysis of above-mentioned domestic and foreign literature it is found that existing RAIM technology has the disadvantages that (1) passes RAIM algorithm of uniting is established on the basis of linearizing observational equation, and precision is still improved space because rounding error;(2) traditional RAIM algorithm is concerned only with current time information, is not concerned with time-domain information, therefore can ignore the control information in time domain.Traditional RAIM algorithm calculates be independently of each historical data every time, therefore then reacts not sensitive enough to the failure in stationary random process, that is, Traditional RAIM algorithm cannot find the failure that error is gradually increased at any time in time, to unmodeled in fault satellites exclusion process Fault type is difficult to be detected.
Summary of the invention
Goal of the invention of the invention is the deficiency for above-mentioned background technique, and it is independently intact to provide a kind of GNSS receiver Property monitoring method, by combine volume Kalman filtering algorithm and information extrapolation on-line study fault type, in terms of less Calculation amount and suitable precision preferably monitor small gradual pseudorange biases, realize the detection to slow unmodeled failure, solve Integrity monitoring is because being difficult to extract the technical issues of feature of every kind of satellite failure is without adapting to fault detection demand.
The present invention adopts the following technical scheme that for achieving the above object
A kind of GNSS receiver autonomous integrity monitoring method filters the shape of algorithm prediction navigation system using volume Kalman State estimated value simultaneously updates measurement information, according to before volume Kalman's Navigation Filter several epoch it is new breath and its variance-association Variance matrix carries out the test statistics at new breath extrapolation building each extrapolation moment, collect the test statistics at multiple extrapolation moment with Integrated database is constructed, the test statistics at current extrapolation moment is differentiated using machine learning method and then determines that fault satellites are compiled Number, when machine learning method differentiates current extrapolation moment failure, the Mishap Database in integrated database is updated, from updated The rule for being used for machine learning next time is extracted in integrated database.
Scheme is advanced optimized as GNSS receiver autonomous integrity monitoring method, is calculated using volume Kalman filtering The state estimation of method prediction navigation system simultaneously updates measurement information comprising time update and measures update:
Time updates:
To the variance-covariance matrix P of k epoch state errorkCholesky is carried out to decompose to obtain k epoch state error The lower triangular matrix S of variance-covariance matrixk, according to expression formula:Choose the of k epoch state estimation I volume point For the state estimation of k epoch, ξiFor the weight parameter of i-th of volume point,2n is volume point sum, k epoch state estimation I-th of volume pointI-th of volume point of the k+1 epoch state estimation transmitted through state transition functionForF () is the state transition function of navigation system, wk、ΓkFor the system noise and its drive of k epoch Dynamic matrix, according to the state estimation of the k+1 epoch of k epoch status predicationForK+1 is gone through The variance-covariance battle array P of first state errork+1/kFor QkFor the system noise variance matrix of k epoch;
It measures and updates:
To the variance-covariance matrix P of k+1 epoch state errork+1/kCholesky is carried out to decompose to obtain k+1 epoch state The lower triangular matrix S of the variance-covariance matrix of errork+1/k, i-th of volume point of k+1 epoch state estimationForI-th of volume point of k+1 epoch state estimationFunction passes are measured to obtain K+1 epoch measure i-th of volume point of estimated valueForH () is the amount of navigation system Survey function, vkFor the measurement noise of k epoch, the measurement estimated value of k+1 epochForK+1 epoch State estimationForFor k+1 epoch measuring value and Measure estimated value, Kk+1For the filtering gain of k+1 epoch, For k+1 epoch prediction it is mutual Variance-covariance battle array is closed,
As the further prioritization scheme of GNSS receiver autonomous integrity monitoring method, navigated according to volume Kalman New breath and its variance-covariance battle array before filter several epoch carry out the inspection system at new breath extrapolation building each extrapolation moment The expression formula of metering are as follows: Tavg,tIt is outer Push away the test statistics of moment t, vk-m、Ak-mIt is newly ceased for k-m epoch and its variance-covariance battle array, Aavg -1The sum average value of variance-covariance inverse matrix is newly ceased for k-m epoch,Rk-mFor the measuring noise square difference battle array of k-m epoch.
Scheme is advanced optimized as GNSS receiver autonomous integrity monitoring method, is worked as using the differentiation of machine learning method The method of the test statistics at preceding extrapolation moment and then determining fault satellites number is to be determined using machine learning rule current outer When pushing away moment failure, newly ceasing the corresponding satellite number of greatest member on variance matrix-covariance matrix leading diagonal is fault satellites Number.
As the scheme that advanced optimizes of GNSS receiver autonomous integrity monitoring method, the event in integrated database is updated The method of barrier database is when the fault data at current extrapolation moment Mishap Database least significant end moment is replaced with current extrapolation The test statistics at quarter.
As the scheme that advanced optimizes of GNSS receiver autonomous integrity monitoring method, measurement information is the pseudorange of GPS Observation.
A kind of computer readable storage medium, is stored thereon with computer program, realization when which is executed by processor Above-mentioned GNSS receiver autonomous integrity monitoring method.
A kind of GNSS receiver, comprising: the meter that memory, processor and storage are run on a memory and on a processor Calculation machine program, processor realize above-mentioned GNSS receiver autonomous integrity monitoring method when executing computer program.
The present invention by adopting the above technical scheme, has the advantages that
(1) present invention is using the monitoring for realizing GNSS receiver autonomous integrity based on CKF and new breath extrapolation, by newly ceasing The new variance-covariance battle array for ceasing and newly ceasing generated in extrapolation process constructs the test statistics at each extrapolation moment, can examine Small gradual pseudorange biases, dynamic is updated for generating machine while differentiating each epoch test statistics using machine learning method The Mishap Database of learning rule, the judging rules of the unmodeled fault type of on-line study, the faulty model without establishing, Without extracting, faulty feature, that is, be suitable for the requirement of complete detection failure.
(2) it can adapt to the non-linear behavior of navigation system using volume Kalman filtering algorithm estimating system state, no It only avoids in EKF algorithm and the first approximation of nonlinear system is handled, and estimated accuracy can be made to be promoted to three ranks, be suitable for The UAV system of high dynamic.
Detailed description of the invention
Fig. 1 is the flow chart of GNSS receiver autonomous integrity monitoring method disclosed in the present application.
Specific embodiment
It is described in detail below with reference to technical solution of the Fig. 1 to invention.
Core content of the present invention is the RAIM algorithm based on CKF and new breath extrapolation, and is produced according in new breath extrapolation process New breath and its variance-covariance battle array before raw several epoch establish the test statistics at extrapolation moment, with test statistics For the database of fault verification, real-time update dicision rules while whether determining failure, with calculation amount, less, precision is fitted Preferably, the advantages of monitoring small gradual pseudorange biases suitable for nonlinear system and energy high confidence level.
CKF (Cubature Kalman Filter, volume Kalman filtering) algorithm relies on gaussian filtering frame, according to three Rank spherical surface-phase diameter volume rule, Posterior Mean and variance-association side from the angle approximate non-linear function passes of numerical integration Difference.By the priori mean value and variance-covariance of state, volume point is chosen by volume rule, volume point is through nonlinear function Transmitting, then do weighting processing and come approximation state Posterior Mean and variance-covariance.
Core of the invention content is filter building, considers following non-linear Gaussian system:
xk=f (xk-1)+Γk-1wk-1(1),
Zk=h (xk)+vk(2),
Wherein, xkFor system mode vector;ZkTo measure vector;F () and h () is respectively the state of nonlinear system Transfer function and measurement function;Γk-1For the driving matrix of system noise;wk-1And vkRespectively system noise and measurement noise, It is white Gaussian noise and uncorrelated.
(1) CKF algorithm
Time updates:
Known k moment system mode vector xkStatistical property beTo this, the state vector at moment is missed The variance-covariance battle array P of differencekCarrying out Cholesky decomposition has: Pk=SkSk T, enable ξ1, ξ2..., ξ2nFor the point of 2n equal weight,
Choose volume point are as follows:
The volume point transmitted by system are as follows:
Obtain the status predication value at k+1 moment:
QkFor the system noise variance matrix at k moment, the State error variance-covariance matrix at k+1 moment are as follows:
It measures and updates:
To Pk+1/kCarrying out Cholesky decomposition has: Pk+1/k=Sk+1/kSk+1/k T, calculating volume point has:
Have by measuring function passes volume point:
Obtain the observation predicted value at k+1 moment:
Obtain the error in measurement variance-covariance battle array at k+1 moment:
K+1 moment one-step prediction cross-correlation variance-covariance battle array are as follows:
The filtering gain at k+1 moment are as follows:
The state estimation at k+1 moment are as follows:
The State error variance at k+1 moment-covariance estimated value are as follows:
(2) system modelling
1. determining state equation:
Building system mode vector is xk=[r τ d τ]T, wherein r is no-manned plane three-dimensional spatial position vector, and τ is to receive Machine clock deviation, d τ are receiver clock frequency drift value.
At this point, first-order system state-transition matrix are as follows:
2. determining measurement equation:
zk=h (xk)+vk(15),
Measuring value is the Pseudo-range Observations of GPS, it may be assumed that
Wherein, ρiFor Pseudo-range Observations, (Xi, Yi, Zi) it is position of i-th moonscope moment in space;(X, Y, Z) Moment is observed in the position in space for receiver;VionAnd VtropRespectively Ionospheric delay correcting number and tropospheric delay correction Number;C is the light velocity;WithRespectively receiver clock-offsets and satellite clock correction;εiTo measure noise.
3. determining test statistics based on new breath extrapolation:
For small gradual pseudorange biases, the present invention improves traditional new breath detection method, using fusion engineering The new breath extrapolation practised carries out the monitoring of receiver consistency.In newly breath extrapolation, using two volume Kalman filter, One is Navigation Filter: under system non-failure conditions, whole measurement informations are updated for state, and filter result is only used for leading Boat;Secondly being Recursive Filtering device: the state recursion before collecting T epoch of Navigation Filter obtains the system shape of current time k State information saves the new breath generated in recursive process.
By the system state estimation before T epochWith State error variance-covariance matrix Pk-TWhen being extrapolated to current K is carved, the state estimation and State error variance-covariance matrix generated in extrapolation process is denoted asWithWherein, 1≤m≤ T。
When system fault-free, Normal Distribution, v are newly ceasedk~N (0, Vk), VkFor the variance matrix of multiple normal distribution.Newly BreathVariance-covariance matrix be Ak-m, convolution (10) has,
Provide the test statistics of extrapolation moment t
Wherein:
4. the realization of the fault detection based on machine learning method:
Present invention determine that the method for failure is different from carrying out the conventional method of fault verification using threshold value comparison, using machine Device learning method determines whether there is failure, and the fault type new based on big data on-line study, has lower false detection rate With higher confidence level.
1) Primary Stage Data prepares: the unmanned plane during flying of (such as urban canyons) is corresponded under scene according to application scenarios pre-acquired Data newly cease v according to above-mentioned volume kalman filter method building systemk, collect the test statistics T at each extrapolation momentAvg, t, Establish data warehouse under initial line huge enough;
2) data markers and training: according to the satellite health status parameter decision satellite in satellite ephemeris with the presence or absence of event Barrier, label failure extrapolate test statistics corresponding to the moment as " 1 ", and test statistics corresponding to the fault-free extrapolation moment is " 0 " is marked to construct integrated database D0 database under initial line, wherein including initial Mishap Database D1 (each data markers of the inside are " 1 ") and normal data library D2 (each data markers of the inside are " 0 ");
3) it fault verification: completes primary new breath extrapolation every and is currently extrapolated the test statistics T of moment tAvg, tAfterwards, it adopts With the subordinating degree function based on fuzzy logic to the test statistics T of current extrapolation moment tAvg, tIt is marked, works as inspection statistics Measure TAvg, tWhen being marked as " 0 ", then it is judged to fault-free, when being marked as " 1 ", is then judged to failure, is marked as the inspection system of " 1 " Metering will be used for the real-time update of integrated database, be extracted from updated integrated database using neural metwork training device Out 4) fuzzy rule for next extrapolation moment fault verification, the update of integrated database are detailed in.It is new to cease when system fault-free The equal Normal Distribution of each element in sequence.When a fault has occurred, the variance of corresponding residual error error increases.Therefore, work as system In there are when failure, it is new that cease the corresponding satellite number of greatest member on variance-covariance battle array leading diagonal be that fault satellites are compiled Number.Failure determine after, system by so that complete paired fault alarm.
4) real-time update is used for the Mishap Database of neural metwork training: will currently extrapolate moment Mishap Database D1tMost The fault data at end moment replaces with the data that " 1 " is labeled as by the fuzzy rule that currently extrapolation moment neural network extracts, Obtain updated Mishap Database D1t+1, then merge the normal data library D2 remained unchanged and rebuild for neural network instruction Practice the integrated database D0 that device is further trainedt+1, to realize the online updating of failure tranining database.

Claims (8)

1. a kind of GNSS receiver autonomous integrity monitoring method, which is characterized in that filter algorithm prediction using volume Kalman and lead The state estimation of boat system simultaneously updates measurement information, according to before volume Kalman's Navigation Filter several epoch it is new breath and Its variance-covariance battle array carries out the test statistics at new breath extrapolation building each extrapolation moment, collects the inspection at multiple extrapolation moment The test statistics and then determining event that statistic is tested to construct integrated database, using the machine learning method differentiation current extrapolation moment Hinder satellite number, when machine learning method differentiates current extrapolation moment failure, updates the Mishap Database in integrated database, from The rule for being used for machine learning next time is extracted in updated integrated database.
2. a kind of GNSS receiver autonomous integrity monitoring method according to claim 1, which is characterized in that use volume card The state estimation of Kalman Filtering algorithm prediction navigation system simultaneously updates measurement information comprising time update and measures update:
Time updates:
To the variance-covariance matrix P of k epoch state errorkCholesky is carried out to decompose to obtain the variance-of k epoch state error The lower triangular matrix S of covariance matrixk, according to expression formula:Choose i-th of appearance of k epoch state estimation Plot pointFor the state estimation of k epoch, ξiFor the weight parameter of i-th of volume point,2n is volume point sum, k epoch state estimation I-th of volume pointI-th of volume point of the k+1 epoch state estimation transmitted through state transition functionForF () is the state transition function of navigation system, wk、ΓkFor k epoch system noise and It drives matrix, according to the state estimation of the k+1 epoch of k epoch status predicationFork The variance-covariance battle array P of+1 epoch state errork+1/kForQkFor the system noise variance matrix of k epoch;
It measures and updates:
To the variance-covariance matrix P of k+1 epoch state errork+1/kCholesky is carried out to decompose to obtain k+1 epoch state error Variance-covariance matrix lower triangular matrix Sk+1/k, i-th of volume point of k+1 epoch state estimationForI-th of volume point of k+1 epoch state estimationFunction passes are measured to obtain K+1 epoch measure i-th of volume point of estimated valueForH () is the amount of navigation system Survey function, vkFor the measurement noise of k epoch, the measurement estimated value of k+1 epochForK+1 epoch State estimationFor For the measuring value and amount of k+1 epoch Survey estimated value, Kk+1For the filtering gain of k+1 epoch, For k+1 epoch prediction it is mutual Variance-covariance battle array is closed,
3. a kind of GNSS receiver autonomous integrity monitoring method according to claim 2, which is characterized in that according to volume card New breath and its variance-covariance battle array before Germania Navigation Filter several epoch carry out new breath extrapolation and construct each extrapolation moment Test statistics expression formula are as follows: TAvg, tFor the test statistics for the moment t that extrapolates, vk-m、Ak-mIt is newly ceased for k-m epoch and its variance-covariance battle array,Aavg -1The sum average value of variance-covariance inverse matrix is newly ceased for k-m epoch,Rk-mFor the measuring noise square difference battle array of k-m epoch.
4. a kind of GNSS receiver autonomous integrity monitoring method according to claim 1, which is characterized in that use engineering Habit method differentiates the test statistics at current extrapolation moment and then the method for determining fault satellites number is, using machine learning rule When determining current extrapolation moment failure, newly ceasing the corresponding satellite number of greatest member on variance matrix-covariance matrix leading diagonal is For fault satellites number.
5. a kind of GNSS receiver autonomous integrity monitoring method according to claim 1, which is characterized in that update comprehensive number Method according to the Mishap Database in library is to replace with the fault data at current extrapolation moment Mishap Database least significant end moment The test statistics at current extrapolation moment.
6. a kind of GNSS receiver autonomous integrity monitoring method according to claim 1, which is characterized in that the measurement letter Breath is the Pseudo-range Observations of GPS.
7. a kind of computer readable storage medium, is stored thereon with computer program, which is characterized in that the program is held by processor GNSS receiver autonomous integrity monitoring method described in claim 1 is realized when row.
8. a kind of GNSS receiver, comprising: the calculating that memory, processor and storage are run on a memory and on a processor Machine program, which is characterized in that the processor realizes GNSS receiver described in claim 1 when executing the computer program Autonomous integrity monitoring method.
CN201910052794.9A 2019-01-21 2019-01-21 GNSS receiver autonomous integrity monitoring method Active CN109490916B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910052794.9A CN109490916B (en) 2019-01-21 2019-01-21 GNSS receiver autonomous integrity monitoring method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910052794.9A CN109490916B (en) 2019-01-21 2019-01-21 GNSS receiver autonomous integrity monitoring method

Publications (2)

Publication Number Publication Date
CN109490916A true CN109490916A (en) 2019-03-19
CN109490916B CN109490916B (en) 2020-01-17

Family

ID=65714825

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910052794.9A Active CN109490916B (en) 2019-01-21 2019-01-21 GNSS receiver autonomous integrity monitoring method

Country Status (1)

Country Link
CN (1) CN109490916B (en)

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110068840A (en) * 2019-05-15 2019-07-30 北京航空航天大学 A kind of ARAIM fault detection method based on pseudo range measurement characteristics extraction
CN110109162A (en) * 2019-03-26 2019-08-09 西安开阳微电子有限公司 A kind of Kalman filtering positioning calculation method that GNSS receiver is adaptive
CN110175682A (en) * 2019-04-22 2019-08-27 广东技术师范大学 A kind of optimization core pivot element analysis fault monitoring method based on Chaos-Particle Swarm Optimization
CN110321366A (en) * 2019-06-27 2019-10-11 中国科学院自动化研究所 Method and system are determined based on the statistic of on-line study
CN110398764A (en) * 2019-08-01 2019-11-01 上海交通大学 A kind of Beidou Navigation System two-band fusion and positioning method and device
CN110596736A (en) * 2019-10-15 2019-12-20 中国电子科技集团公司第五十四研究所 GNSS observation abnormal value detection and isolation method
CN111060133A (en) * 2019-12-04 2020-04-24 南京航空航天大学 Integrated navigation integrity monitoring method for urban complex environment
CN111323793A (en) * 2020-03-30 2020-06-23 中国矿业大学 GNSS pseudo-range single-point positioning state domain integrity monitoring method
CN112230247A (en) * 2020-09-25 2021-01-15 南京航空航天大学 GNSS integrity monitoring method used in urban complex environment
CN112526549A (en) * 2020-12-01 2021-03-19 北京航空航天大学 Method and system for identifying integrity fault of foundation enhancement system
CN112558112A (en) * 2020-11-17 2021-03-26 中国矿业大学 GNSS state domain slow-change slope fault integrity monitoring method
CN112733315A (en) * 2019-10-29 2021-04-30 千寻位置网络有限公司 Integrity risk verification method and system
CN113283861A (en) * 2021-05-18 2021-08-20 上海示右智能科技有限公司 Method for constructing intelligent enterprise compliance
CN113406670A (en) * 2021-06-17 2021-09-17 哈尔滨工程大学 Beidou non-motorized broadcast ephemeris fault real-time monitoring method
CN113670337A (en) * 2021-09-03 2021-11-19 东南大学 Method for detecting slow-changing fault of GNSS/INS combined navigation satellite
CN115291253A (en) * 2022-08-02 2022-11-04 东北大学 Vehicle positioning integrity monitoring method and system based on residual error detection
CN116881088A (en) * 2023-09-06 2023-10-13 长沙金维信息技术有限公司 System monitoring method and device, storage medium and electronic device

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2648018A1 (en) * 2012-04-02 2013-10-09 Astrium GmbH An improved RAIM algorithm
CN104344970A (en) * 2014-10-13 2015-02-11 中国船舶重工集团公司第七二六研究所 Equipment system fault integral diagnosis method and device
CN107193016A (en) * 2017-04-06 2017-09-22 广州中硕云空间信息技术有限公司 A kind of method and system of city GNSS navigation quality evaluations and prediction
CN107479069A (en) * 2017-08-01 2017-12-15 天津博创金成技术开发有限公司 A kind of slow change slope failure completeness monitoring method

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2648018A1 (en) * 2012-04-02 2013-10-09 Astrium GmbH An improved RAIM algorithm
CN104344970A (en) * 2014-10-13 2015-02-11 中国船舶重工集团公司第七二六研究所 Equipment system fault integral diagnosis method and device
CN107193016A (en) * 2017-04-06 2017-09-22 广州中硕云空间信息技术有限公司 A kind of method and system of city GNSS navigation quality evaluations and prediction
CN107479069A (en) * 2017-08-01 2017-12-15 天津博创金成技术开发有限公司 A kind of slow change slope failure completeness monitoring method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
贾昌磊: "GPS/GLONASS/BDS多星座导航接收机自主完好性监视", 《中国优秀硕士学位论文全文数据库 信息科技辑》 *
钟丽娜: "惯性/卫星组合导航自适应容错与精度性能增强关键技术研究", 《中国博士学位论文全文数据库 信息科技辑》 *

Cited By (25)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110109162A (en) * 2019-03-26 2019-08-09 西安开阳微电子有限公司 A kind of Kalman filtering positioning calculation method that GNSS receiver is adaptive
CN110175682A (en) * 2019-04-22 2019-08-27 广东技术师范大学 A kind of optimization core pivot element analysis fault monitoring method based on Chaos-Particle Swarm Optimization
CN110068840B (en) * 2019-05-15 2020-12-29 北京航空航天大学 ARAIM fault detection method based on pseudo-range measurement characteristic value extraction
CN110068840A (en) * 2019-05-15 2019-07-30 北京航空航天大学 A kind of ARAIM fault detection method based on pseudo range measurement characteristics extraction
CN110321366B (en) * 2019-06-27 2021-06-25 中国科学院自动化研究所 Statistical quantity determining method and system based on online learning
CN110321366A (en) * 2019-06-27 2019-10-11 中国科学院自动化研究所 Method and system are determined based on the statistic of on-line study
CN110398764A (en) * 2019-08-01 2019-11-01 上海交通大学 A kind of Beidou Navigation System two-band fusion and positioning method and device
CN110596736A (en) * 2019-10-15 2019-12-20 中国电子科技集团公司第五十四研究所 GNSS observation abnormal value detection and isolation method
CN110596736B (en) * 2019-10-15 2021-04-02 中国电子科技集团公司第五十四研究所 GNSS observation abnormal value detection and isolation method
CN112733315A (en) * 2019-10-29 2021-04-30 千寻位置网络有限公司 Integrity risk verification method and system
CN111060133A (en) * 2019-12-04 2020-04-24 南京航空航天大学 Integrated navigation integrity monitoring method for urban complex environment
CN111323793B (en) * 2020-03-30 2021-02-05 中国矿业大学 GNSS pseudo-range single-point positioning state domain integrity monitoring method
CN111323793A (en) * 2020-03-30 2020-06-23 中国矿业大学 GNSS pseudo-range single-point positioning state domain integrity monitoring method
CN112230247A (en) * 2020-09-25 2021-01-15 南京航空航天大学 GNSS integrity monitoring method used in urban complex environment
CN112558112A (en) * 2020-11-17 2021-03-26 中国矿业大学 GNSS state domain slow-change slope fault integrity monitoring method
CN112558112B (en) * 2020-11-17 2022-09-16 中国矿业大学 GNSS state domain slow-change slope fault integrity monitoring method
CN112526549A (en) * 2020-12-01 2021-03-19 北京航空航天大学 Method and system for identifying integrity fault of foundation enhancement system
CN113283861A (en) * 2021-05-18 2021-08-20 上海示右智能科技有限公司 Method for constructing intelligent enterprise compliance
CN113283861B (en) * 2021-05-18 2024-04-16 上海示右智能科技有限公司 Method for constructing compliance of intelligent enterprise
CN113406670A (en) * 2021-06-17 2021-09-17 哈尔滨工程大学 Beidou non-motorized broadcast ephemeris fault real-time monitoring method
CN113670337A (en) * 2021-09-03 2021-11-19 东南大学 Method for detecting slow-changing fault of GNSS/INS combined navigation satellite
CN115291253A (en) * 2022-08-02 2022-11-04 东北大学 Vehicle positioning integrity monitoring method and system based on residual error detection
CN115291253B (en) * 2022-08-02 2023-12-05 东北大学 Vehicle positioning integrity monitoring method and system based on residual error detection
CN116881088A (en) * 2023-09-06 2023-10-13 长沙金维信息技术有限公司 System monitoring method and device, storage medium and electronic device
CN116881088B (en) * 2023-09-06 2023-11-28 长沙金维信息技术有限公司 System monitoring method and device, storage medium and electronic device

Also Published As

Publication number Publication date
CN109490916B (en) 2020-01-17

Similar Documents

Publication Publication Date Title
CN109490916A (en) A kind of GNSS receiver autonomous integrity monitoring method
CN109521443B (en) Method for detecting ephemeris abnormity
CN107807371B (en) Data dependence processing method based on big-dipper satellite observation data
CN103760571B (en) For the fragility monitoring method based on influence factor feature of GPS
CN111060133B (en) Integrated navigation integrity monitoring method for urban complex environment
US20090273511A1 (en) Method of operating a satellite navigation receiver
CN109308518B (en) Monitoring system based on probabilistic neural network and smoothing parameter optimization method thereof
Pesonen A framework for Bayesian receiver autonomous integrity monitoring in urban navigation
Zhu et al. Extended Kalman filter (EKF) innovation-based integrity monitoring scheme with C/N 0 weighting
CN114545454A (en) Fusion navigation system integrity monitoring method for automatic driving
KR102203960B1 (en) Gnss integrity monitoring method based on machine learning for satellite navigation augmentation system
Du et al. A comparative study on active SLAM and autonomous exploration with particle filters
de Oliveira et al. Recent advances in sensor integrity monitoring methods—A review
CN115420284A (en) Fault detection and identification method for combined navigation system
CN103926596A (en) Steady GNSS anti-spoofing method based on particle filter
CN115561782B (en) Satellite fault detection method in integrated navigation based on odd-even vector projection
Ahn et al. GPS integrity monitoring method using auxiliary nonlinear filters with log likelihood ratio test approach
Favenza et al. A machine learning approach to GNSS scintillation detection: Automatic soft inspection of the events
Chen et al. A SVM based GNSS performance assessment with reliable vulnerability degree model
Harbaoui et al. Navigation context adaptive fault detection and exclusion strategy based on deep learning & information theory: Application to a gnss/imu integration
Gao et al. A method for estimating the positioning accuracy of vehicles in urban environments
CN111999750A (en) Real-time single-station cycle slip detection improvement method aiming at inaccurate lever arm
Gao et al. An improved detection method of GNSS faults with fractional information divergence
No et al. Machine learning based overbound modeling of multipath error for safety critical urban environment
Tanil et al. Kalman filter partial innovation sequence monitor

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