CN106125100A - A kind of spacing wave seriality appraisal procedure - Google Patents
A kind of spacing wave seriality appraisal procedure Download PDFInfo
- Publication number
- CN106125100A CN106125100A CN201610421087.9A CN201610421087A CN106125100A CN 106125100 A CN106125100 A CN 106125100A CN 201610421087 A CN201610421087 A CN 201610421087A CN 106125100 A CN106125100 A CN 106125100A
- Authority
- CN
- China
- Prior art keywords
- bdt
- seriality
- sequence
- interruption
- spacing wave
- 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.)
- Pending
Links
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/20—Integrity monitoring, fault detection or fault isolation of space segment
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Security & Cryptography (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention discloses a kind of spacing wave seriality appraisal procedure.Use the present invention can utilize navigational satellite system regular maintenance data, downstream signal Monitoring Data and reliability rationale, provide a kind of operable, effective ways of easy realization for the assessment of spacing wave continuity parameter.First the present invention obtains system day-to-day operation and interrupts the most healthy word information of data and receiver actual measurement, situation is interrupted according to the plan analyzed in day-to-day operation, obtain the running status sequence under the influence of the most unplanned interruption and the healthy word sequence of actual measurement, the seriality risk of two kinds of data sources is obtained respectively according to Theoretical Prediction method and actual measurement assessment method, and utilize the running status sequence under the influence of unplanned interruption to use actual measurement statistic law to obtain seriality risk, finally three kinds of Risk Calculation results of contrast, are estimated the seriality of spacing wave.
Description
Technical field
The present invention relates to navigation system reliability security technique field, be specifically related to a kind of spacing wave seriality assessment
Method.
Background technology
After satellite navigation system is built up, systems face provides open service performance to promise to undertake for users, maintains indispensability
The tasks such as energy index.In addition to precision index, the user such as International Civil Aviation, maritime affairs require indispensable performance indications also include integrity,
Seriality and availability.
At present, focus primarily upon the assessment to service continuity for the successional assessment of satellite navigation system, 2008
The U.S. takes the lead in being classified as spacing wave seriality the specification index of assessment GPS spacing wave performance, and navigation and relevant should
Approve with field, but there is no clear and definite, unified spacing wave seriality appraisal procedure.At present, US Federal Aviation Administration (FAA) exists
To the spacing wave continuity parameter be given in GPS performance assessment report quarterly, the method for employing is equivalent to availability and comments
Estimate method, quite different with successional definition and intension.European Union Galileo is in system design document, it is proposed that with reliably
The method that degree formula calculates continuity parameter.Chinese satellite positioning system has opened operation, runs inside and outside prison in system
Survey two aspects and all can obtain signal interruption data, how these data are effectively utilized, accurate evaluation and prediction seriality
Index, needs a kind of ground spacing wave continuity parameter appraisal procedure that realizes operable, easy badly and instructs the prediction to this index and comment
Estimate.
Summary of the invention
In view of this, the invention provides a kind of spacing wave seriality appraisal procedure, it is possible to for spacing wave seriality
Index evaluation provides a kind of operable, effective ways of easy realization.
The spacing wave seriality appraisal procedure of the present invention, comprises the steps:
Step 1, the interruption data operated according to the ground system day-to-day operation of triones navigation system, obtain at dipper system
Spacing wave upstate sequence S_AV (t under time referenceBDT), S_AV (tBDT)={ (tBDT, ' X ') }, wherein, tBDTFor north
During bucket;X is ' 0 ' or ' 1 ', and ' 0 ' represents available, and ' 1 ' represents unavailable;
And according to interruption source, by S_AV (tBDT) it is divided into plan interrupt sequence S_TS (tBDT) and unplanned interrupt sequence S_
TU(tBDT);
Step 2, calculates successional Theoretical Prediction value C caused by unplanned interruptioncal-1:
Ccal-1=exp (-T/MTBO)
Wherein, T is the seriality calculation window set;MTBOiFor under each interruption source
Mean time between interruptions;N is total number of the interruption source of unplanned interruption;
Step 3, gathers aeronautical satellite and can follow the tracks of the healthy word information in segmental arc down navigation signal, it is thus achieved that during dipper system
Between healthy word state sequence S_H (t under benchmarkBDT), S_H (tBDT)={ (tBDT, ' Y ') }, wherein, Y is ' 0 ' or ' 1 ', ' 0 ' table
Showing health, ' 1 ' represents unhealthy;
Step 4, according to healthy word state sequence, calculates the measured value C of the seriality risk caused by unplanned interruptionmon:
Wherein, [tstart,tend] it is testing time section;tpFor the subscriber computer sampling interval;In bool () function, if t,
Healthy word state is healthy, then bool () functional value is 1;Otherwise bool () functional value is 0;S_HTU(tBDT) for getting rid of in the works
Disconnected healthy word state sequence, S_HTU(tBDT)=S_H (tBDT)-S_TS(tBDT);
Step 5, according to the unplanned interrupt sequence S_TU (t in the signal upstate sequence in step 1BDT), the company of calculating
Theoretical calculation value C of continuous propertycal-2:
Step 6, to Theoretical Prediction value Ccal-1, theoretical calculation value Ccal-2With measured value CmonContrast, it is achieved to space
The seriality assessment of signal.
Beneficial effect:
(1) present invention solves satellite navigation system spacing wave continuity parameter and comments to test from concept, theoretical method
The problem estimating realization.
(2) according to the data source difference of different user, propose respectively and be applicable to day-to-day operation data and reliability design
The theoretical calculation method of data, and it is applicable to the practical test statistics method of domestic consumer's receiver measured data, comment conveniently and efficiently
Estimate, add up and determine spacing wave continuity parameter.
(3) utilize this method and appraisal procedure, be possible not only to determine spacing wave continuity parameter, moreover it is possible to expanded application in
Spacing wave difference seriality, and service continuity index determines, supplements and improves Chinese satellite positioning system service reliability
Field Monitoring Indexes evaluation system, improves, for system state monitoring and design, the professional technique support providing important.
Accompanying drawing explanation
Fig. 1 is flow chart of the present invention.
Fig. 2 is two kinds of method seriality statistical result comparison diagrams.
Fig. 3 is two kinds of data source seriality statistical result comparison diagrams.
Detailed description of the invention
Develop simultaneously embodiment below in conjunction with the accompanying drawings, describes the present invention.
The invention provides a kind of spacing wave seriality appraisal procedure, utilize navigational satellite system regular maintenance data,
Downstream signal Monitoring Data and reliability rationale, provide a kind of operable, Yi Shi for the assessment of spacing wave continuity parameter
Existing effective ways.
First the present invention obtains system day-to-day operation and interrupts the most healthy word information of data and receiver actual measurement, according to dividing
Situation is interrupted in plan in analysis day-to-day operation, obtains the running status sequence under the influence of the most unplanned interruption and the healthy word sequence of actual measurement
Row, obtain the seriality risk of two kinds of data sources respectively, and utilize unplanned interruption according to Theoretical Prediction method and actual measurement assessment method
Under the influence of running status sequence use actual measurement statistic law to obtain seriality risk, finally three kinds of Risk Calculation results of contrast.Tool
Body flow process is as it is shown in figure 1, comprise the steps:
Step 1, acquisition system run signal is interrupted data, is obtained signal interruption sequence.
In this step, collect the interruption data of triones navigation system ground system day-to-day operation operation, interrupt data and include
Interrupt starting to describe with end time, interruption source;Under dipper system time reference, obtain spacing wave upstate sequence
S_AV(tBDT), S_AV (tBDT)={ (tBDT, ' X ') }, wherein, X is ' 0 ' or ' 1 ', and ' 0 ' represents available, and ' 1 ' represents unavailable.
Meanwhile, according to interruption source, spacing wave upstate sequence is divided into plan interrupt sequence S_TS (tBDT) and unplanned interruption
Sequence S_TU (tBDT);
Wherein, the dipper system time uses " second in week+week " (" WN+SOW ") mode to represent, WN represents with " wwww ", SOW
Represent with " ssssss ", i.e. t during the Big DipperBDTIt is expressed as " wwww.ssssss ", signal upstate sequence S_AV (tBDT)=
{ (wwww.ssssss, ' X ') }, X ∈ { 0,1}
Step 2, the unplanned interrupt sequence obtained according to step 1, use formula of reliability to calculate by unplanned
The disconnected seriality Risk Theory discreet value caused:
Step 2.1, classifies to interruption according to the interruption source of unplanned interruption, and calculates each respectively by interrupt type
The mean time between interruptions MTBO that type is interruptedi, subscript i represents i-th kind of interruption source;Time between interruptions be of the same race in
Under disconnected type, the interval time between adjacent twice interruption.
Step 2.2, according to the mean time between interruptions MTBO of all types of unplanned interruptionsi, calculate under combined influence
MTBO,N represents the total unplanned interrupt type of n kind.
MTBO under combined influence: spacing wave interrupts to be caused by many reasons, every kind of time between interruptions rule
There is difference, according to Reliability Theories, product fault rate under n kind failure cause combined influenceAnd
Step 2.3, according to formula Ccal-1=exp (-T/MTBO) calculates seriality, and wherein, T represents the seriality meter of setting
Calculate window.
Step 3, gathers the spacing wave health word sequence of user side actual measurement.
In this step, utilize receiver round-the-clock to gather aeronautical satellite and can follow the tracks of segmental arc down navigation signal, it is thus achieved that can be with
In track segmental arc signal, real-time, continuous print health word information, obtains under dipper system time reference, healthy word state sequence, is denoted as
S_H(tBDT), S_H (tBDT)={ (tBDT, ' Y ') }, wherein, Y is ' 0 ' or ' 1 ', and ' 0 ' represents health, and ' 1 ' represents unhealthy.
Down navigation signal: comprise pseudo range observed quantity and navigation message data;
Navigation message data: comprise signal launch time the Big Dipper time, satellite health word, satellite ephemeris, clock correction, ionosphere mould
Shape parameter etc.;
Healthy word: satellite-signal health word only has the binary digit of 1 bit, span be 0,1}, ' 0 ' represents and can use,
' 1 ' represents unavailable;
Healthy word sequence state: S_H (tBDT)={ (wwww.ssssss, ' Y ') }, Y ∈ { 0,1}.
Step 4, according to the receiver measured data in step 3, uses slip recurrence formula to calculate and is drawn by unplanned interruption
The measured value of the seriality risk risen:
Step 4.1, the healthy word state sequence S_H (t of contrastBDT) and the plan interrupt sequence S_ of spacing wave upstate
TS(tBDT), obtain measured data health word state sequence S_H that eliminating plan is interruptedTU(tBDT), wherein, S_HTU(tBDT)=S_
H(tBDT)-S_TS(tBDT), S_HTU(tBDT)={ (tBDT, ' Z ') }, wherein, Z is ' 0 ' or ' 1 ', and ' 0 ' represents health, and ' 1 ' represents
Unhealthy;
Step 4.2, it is assumed that testing time section is [tstart,tend], the subscriber computer sampling interval is tp, use formula
Calculate the measured value of the seriality risk caused by unplanned interruption.Wherein, when T is the window that seriality calculates
Between.In bool () function, if k moment health word state is healthy, then take 1;Otherwise take 0.
Step 5, two kinds of methods of relative analysis, the result of calculation of two kinds of data sources:
Step 5.1, according to the unplanned interrupt sequence S_TU (t in the signal upstate sequence in step 1BDT) carry out
Calculate.Assume that whole timing statistics section is for [tstart,tend], the signal condition sampling interval is tp, use formula
Calculating seriality, T represents the widow time that seriality calculates, if k time-ofday signals state is available, then and bool () letter
Number takes 1, otherwise takes 0.
Step 5.2, is respectively compared Ccal-1、Ccal-2With CmonAcquired results, wherein Ccal-1With Ccal-2For under same data source
The comparing result of two kinds of methods, Ccal-2With CmonFor the comparing result of under same procedure two kinds of data sources, continuous to spacing wave
Property is estimated.
Theoretical calculation method is applicable to grasp system reliability design data and system day-to-day operation data, and lacks reception
The situation of machine measured data, is mainly used according to historical empirical data or designs data to existing satellite in orbit or follow-up new design
Satellite spatial signal seriality is estimated;Practical test statistics method is applicable to carry out long-time Monitoring Data continuously, it is possible to
The domestic consumer of in time notice is interrupted in the plan that obtains, be mainly used according to actual measurement user data to satellite in orbit a period of time
Interior spacing wave seriality is evaluated.The present invention is directed to the reliability field new ideas occurred in satellite navigation system,
According to relating to the definition of parameter, in conjunction with reliability rationale, according to the data source difference of different user, propose applicable respectively
In the theoretical calculation method of day-to-day operation data Yu reliability design data, and it is applicable to domestic consumer's receiver measured data
Actual measurement method, solves the problem determining spacing wave continuity parameter, can be Chinese satellite positioning system service performance
Index Design prediction and test assessment provide important professional technique support.
In sum, these are only presently preferred embodiments of the present invention, be not intended to limit protection scope of the present invention.
All within the spirit and principles in the present invention, any modification, equivalent substitution and improvement etc. made, should be included in the present invention's
Within protection domain.
Claims (1)
1. a spacing wave seriality appraisal procedure, it is characterised in that comprise the steps:
Step 1, the interruption data operated according to the ground system day-to-day operation of triones navigation system, obtain in the dipper system time
Spacing wave upstate sequence S_AV (t under benchmarkBDT), S_AV (tBDT)={ (tBDT, ' X ') }, wherein, tBDTFor the Big Dipper
Time;X is ' 0 ' or ' 1 ', and ' 0 ' represents available, and ' 1 ' represents unavailable;
And according to interruption source, by S_AV (tBDT) it is divided into plan interrupt sequence S_TS (tBDT) and unplanned interrupt sequence S_TU
(tBDT);
Step 2, calculates successional Theoretical Prediction value C caused by unplanned interruptioncal-1:
Ccal-1=exp (-T/MTBO)
Wherein, T is the seriality calculation window set;MTBOiAverage under each interruption source
Time between interruptions;N is total number of the interruption source of unplanned interruption;
Step 3, gathers aeronautical satellite and can follow the tracks of the healthy word information in segmental arc down navigation signal, it is thus achieved that dipper system time base
Healthy word state sequence S_H (t under ZhunBDT), S_H (tBDT)={ (tBDT, ' Y ') }, wherein, Y is ' 0 ' or ' 1 ', and ' 0 ' represents strong
Health, ' 1 ' represents unhealthy;
Step 4, according to healthy word state sequence, calculates the measured value C of the seriality risk caused by unplanned interruptionmon:
Wherein, [tstart,tend] it is testing time section;tpFor the subscriber computer sampling interval;In bool () function, if t, healthy
Word state is healthy, then bool () functional value is 1;Otherwise bool () functional value is 0;S_HTU(tBDT) for getting rid of what plan was interrupted
Healthy word state sequence, S_HTU(tBDT)=S_H (tBDT)-S_TS(tBDT);
Step 5, according to the unplanned interrupt sequence S_TU (t in upstate sequence in step 1BDT), calculate successional theory
Measuring and calculating value Ccal-2:
Step 6, to Theoretical Prediction value Ccal-1, theoretical calculation value Ccal-2With measured value CmonContrast, it is achieved to spacing wave
Seriality is assessed.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610421087.9A CN106125100A (en) | 2016-06-14 | 2016-06-14 | A kind of spacing wave seriality appraisal procedure |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610421087.9A CN106125100A (en) | 2016-06-14 | 2016-06-14 | A kind of spacing wave seriality appraisal procedure |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106125100A true CN106125100A (en) | 2016-11-16 |
Family
ID=57270075
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610421087.9A Pending CN106125100A (en) | 2016-06-14 | 2016-06-14 | A kind of spacing wave seriality appraisal procedure |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106125100A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106610495A (en) * | 2016-11-18 | 2017-05-03 | 中国航天标准化研究所 | Availability, continuity and completeness index distribution method of satellite navigation system |
CN112269194A (en) * | 2020-09-17 | 2021-01-26 | 中国航天标准化研究所 | Navigation satellite system continuity risk assessment method |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103731221A (en) * | 2014-01-09 | 2014-04-16 | 中国航天标准化研究所 | Space and ground integrated network system availability determining method |
CN104215974A (en) * | 2014-09-15 | 2014-12-17 | 中国航天标准化研究所 | Integrity detection availability determining method of satellite navigation system |
-
2016
- 2016-06-14 CN CN201610421087.9A patent/CN106125100A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103731221A (en) * | 2014-01-09 | 2014-04-16 | 中国航天标准化研究所 | Space and ground integrated network system availability determining method |
CN104215974A (en) * | 2014-09-15 | 2014-12-17 | 中国航天标准化研究所 | Integrity detection availability determining method of satellite navigation system |
Non-Patent Citations (3)
Title |
---|
刘帅 等: "GNSS空间信号连续性评估模型与算法", 《导航定位学报》 * |
杨卓鹏 等: "利用检测网中断分析导航***完好性监测", 《导航定位学报》 * |
龚佩佩 等: "空间信号性能理论计算与测试评估方法研究", 《第六届中国卫星导航学术年会论文集—S06北斗/GNSS测试评估技术》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106610495A (en) * | 2016-11-18 | 2017-05-03 | 中国航天标准化研究所 | Availability, continuity and completeness index distribution method of satellite navigation system |
CN106610495B (en) * | 2016-11-18 | 2019-01-08 | 中国航天标准化研究所 | A kind of satellite navigation system availability continuity integrity Allocation method |
CN112269194A (en) * | 2020-09-17 | 2021-01-26 | 中国航天标准化研究所 | Navigation satellite system continuity risk assessment method |
CN112269194B (en) * | 2020-09-17 | 2023-06-30 | 中国航天标准化研究所 | Navigation satellite system continuity risk assessment method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Jiang et al. | Evaluation of Sentinel-3 SRAL SAR altimetry over Chinese rivers | |
Prieto et al. | Flow prediction in ungauged catchments using probabilistic random forests regionalization and new statistical adequacy tests | |
CN103760571B (en) | For the fragility monitoring method based on influence factor feature of GPS | |
CN104238352A (en) | National benchmark controlled local time standard generating system and method | |
Heng et al. | Statistical characterization of GLONASS broadcast ephemeris errors | |
Zhu et al. | Extended Kalman filter (EKF) innovation-based integrity monitoring scheme with C/N 0 weighting | |
CN104035113A (en) | Pseudo-range-based reliable locating method of multimode GNSS receiver | |
CN103235883A (en) | Non-contact river surface velocity measurement based regression calculation method for instantaneous flow rate | |
Gunning et al. | Multi-GNSS constellation anomaly detection and performance monitoring | |
Petrie et al. | Data evaluation for acoustic Doppler current profiler measurements obtained at fixed locations in a natural river | |
CN106125100A (en) | A kind of spacing wave seriality appraisal procedure | |
CN104504247B (en) | A kind of RAIM methods of the double satellite failures of GPS | |
CN113031036B (en) | Ionosphere phase flicker factor construction method based on GNSS 30s sampling frequency data | |
Wang et al. | Enhanced fault detection for GNSS/INS integration using maximum correntropy filter and local outlier factor | |
CN105136126B (en) | The method that tsunami ripple detecting is carried out using deep seafloor pressure data | |
Wang et al. | Convective weather forecast accuracy analysis at center and sector levels | |
Wang et al. | Real-time assessment of GNSS observation noise with single receivers | |
CN102749634A (en) | Pseudo-range acceleration failure detecting method in satellite navigation region reinforcement system | |
Weng et al. | Assessment and mitigation of ionospheric disturbance effects on GPS accuracy and integrity | |
Li et al. | A multi-scale three-dimensional variational data assimilation scheme and its application to coastal oceans | |
Veerman et al. | A tool for GNSS integrity verification based on statistical extreme value theory | |
Roohi | Capability of pulse-limited satellite radar altimetry to monitor inland water bodies | |
CN104330807A (en) | Ionospheric storm monitoring algorithm based on code and carrier separation | |
Over et al. | Uncertainty analysis of index-velocity meters and discharge computations at the Chicago Sanitary and Ship Canal near Lemont, Illinois, water years 2006–16 | |
Perea Diaz | Design of an integrity support message for offline advanced RAIM |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20161116 |
|
WD01 | Invention patent application deemed withdrawn after publication |