CN106125100A - A kind of spacing wave seriality appraisal procedure - Google Patents

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

A kind of spacing wave seriality appraisal procedure

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 reference_BDT), S_AV (t_BDT)={ (t_BDT, ' X ') }, wherein, t_BDTFor north During bucket；X is ' 0 ' or ' 1 ', and ' 0 ' represents available, and ' 1 ' represents unavailable；

And according to interruption source, by S_AV (t_BDT) it is divided into plan interrupt sequence S_TS (t_BDT) and unplanned interrupt sequence S_ TU(t_BDT)；

Step 2, calculates successional Theoretical Prediction value C caused by unplanned interruption_cal-1:

C_cal-1=exp (-T/MTBO)

Wherein, T is the seriality calculation window set；MTBO_iFor 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 benchmark_BDT), S_H (t_BDT)={ (t_BDT, ' 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 interruption_mon:

$C_{mon}=\underset{t_{BDT}=t_{s\mathrm{ta}rt},inc=t_p}{\overset{t_{end}-T}{\Σ}}\left\{\underset{t=t_{BDT},inc=t_p}{\overset{t_{BDT}+T}{\Σ}}bool(S\_H_{TU}(t\left)\right)\right\}/\underset{t_{BDT}=t_{s\mathrm{ta}rt},inc=t_p}{\overset{t_{end}-T}{\Σ}}bool(S\_H_{TU}(t_{BDT}\left)\right)$

Wherein, [t_start,t_end] it is testing time section；t_pFor 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_H_TU(t_BDT) for getting rid of in the works Disconnected healthy word state sequence, S_H_TU(t_BDT)=S_H (t_BDT)-S_TS(t_BDT)；

Step 5, according to the unplanned interrupt sequence S_TU (t in the signal upstate sequence in step 1_BDT), the company of calculating Theoretical calculation value C of continuous property_cal-2:

$C_{cal-2}=\underset{t_{BDT}=t_{s\mathrm{ta}rt},inc=t_p}{\overset{t_{end}-T}{\Σ}}\left\{\underset{t=t_{BDT},inc=t_p}{\overset{t_{BDT}+T}{\Σ}}bool(S\_TU(t\left)\right)\right\}/\underset{t_{BDT}=t_{s\mathrm{ta}rt},inc=t_p}{\overset{t_{end}-T}{\Σ}}bool(S\_TU(t_{BDT}\left)\right)$

Step 6, to Theoretical Prediction value C_cal-1, theoretical calculation value C_cal-2With measured value C_monContrast, 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(t_BDT), S_AV (t_BDT)={ (t_BDT, ' 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 (t_BDT) and unplanned interruption Sequence S_TU (t_BDT)；

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 Dipper_BDTIt is expressed as " wwww.ssssss ", signal upstate sequence S_AV (t_BDT)= { (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 interrupted_i, 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 interruptions_i, 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 C_cal-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(t_BDT), S_H (t_BDT)={ (t_BDT, ' 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 (t_BDT)={ (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 contrast_BDT) and the plan interrupt sequence S_ of spacing wave upstate TS(t_BDT), obtain measured data health word state sequence S_H that eliminating plan is interrupted_TU(t_BDT), wherein, S_H_TU(t_BDT)=S_ H(t_BDT)-S_TS(t_BDT), S_H_TU(t_BDT)={ (t_BDT, ' 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 [t_start,t_end], the subscriber computer sampling interval is t_p, use formula

$C_{mon}=\underset{t_{BDT}=t_{s\mathrm{ta}rt},inc=t_p}{\overset{t_{end}-T}{\Σ}}\left\{\underset{t=t_{BDT},inc=t_p}{\overset{t_{BDT}+T}{\Σ}}bool(S\_H_{TU}(t\left)\right)\right\}/\underset{t_{BDT}=t_{s\mathrm{ta}rt},inc=t_p}{\overset{t_{end}-T}{\Σ}}bool(S\_H_{TU}(t_{BDT}\left)\right)$

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 1_BDT) carry out Calculate.Assume that whole timing statistics section is for [t_start,t_end], the signal condition sampling interval is t_p, use formula

$C_{cal-2}=\underset{t_{BDT}=t_{s\mathrm{ta}rt},inc=t_p}{\overset{t_{end}-T}{\Σ}}\left\{\underset{t=t_{BDT},inc=t_p}{\overset{t_{BDT}+T}{\Σ}}bool(S\_TU(t\left)\right)\right\}/\underset{t_{BDT}=t_{s\mathrm{ta}rt},inc=t_p}{\overset{t_{end}-T}{\Σ}}bool(S\_TU(t_{BDT}\left)\right)$

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 C_cal-1、C_cal-2With C_monAcquired results, wherein C_cal-1With C_cal-2For under same data source The comparing result of two kinds of methods, C_cal-2With C_monFor 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 benchmark_BDT), S_AV (t_BDT)={ (t_BDT, ' X ') }, wherein, t_BDTFor 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 (t_BDT) it is divided into plan interrupt sequence S_TS (t_BDT) and unplanned interrupt sequence S_TU (t_BDT)；

Step 2, calculates successional Theoretical Prediction value C caused by unplanned interruption_cal-1:

C_cal-1=exp (-T/MTBO)

Wherein, T is the seriality calculation window set；MTBO_iAverage 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 Zhun_BDT), S_H (t_BDT)={ (t_BDT, ' 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 interruption_mon:

$C_{mon}=\underset{t_{BDT}=t_{start},inc=t_p}{\overset{t_{end}-T}{\mathrm{\Σ}}}\left\{\underset{t=t_{BDT},inc=t_p}{\overset{t_{BDT}+T}{\mathrm{\Π}}}bool(S\_H_{TU}(t\left)\right)\right\}/\underset{t_{BDT}=t_{start},inc=t_p}{\overset{t_{end}-T}{\mathrm{\Σ}}}bool(S\_H_{TU}(t_{BDT}\left)\right)$

Wherein, [t_start,t_end] it is testing time section；t_pFor 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_H_TU(t_BDT) for getting rid of what plan was interrupted Healthy word state sequence, S_H_TU(t_BDT)=S_H (t_BDT)-S_TS(t_BDT)；

Step 5, according to the unplanned interrupt sequence S_TU (t in upstate sequence in step 1_BDT), calculate successional theory Measuring and calculating value C_cal-2:

$C_{cal-2}=\underset{t_{BDT}=t_{start},inc=t_p}{\overset{t_{end}-T}{\mathrm{\Σ}}}\left\{\underset{t=t_{BDT},inc=t_p}{\overset{t_{BDT}+T}{\mathrm{\Π}}}bool(S\_TU(t\left)\right)\right\}/\underset{t_{BDT}=t_{start},inc=t_p}{\overset{t_{end}-T}{\mathrm{\Σ}}}bool(S\_TU(t_{BDT}\left)\right)$

Step 6, to Theoretical Prediction value C_cal-1, theoretical calculation value C_cal-2With measured value C_monContrast, it is achieved to spacing wave Seriality is assessed.