CN106610495B - A kind of satellite navigation system availability continuity integrity Allocation method - Google Patents
A kind of satellite navigation system availability continuity integrity Allocation method Download PDFInfo
- Publication number
- CN106610495B CN106610495B CN201611027438.4A CN201611027438A CN106610495B CN 106610495 B CN106610495 B CN 106610495B CN 201611027438 A CN201611027438 A CN 201611027438A CN 106610495 B CN106610495 B CN 106610495B
- Authority
- CN
- China
- Prior art keywords
- satellite
- mtbo
- service
- mean time
- availability
- 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.)
- Expired - Fee Related
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
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (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 satellite navigation system availability continuity integrity Allocation methods, include the following steps: the constellation value and service availability according to satellite navigation system, using service availability calculation formula, obtain constellation state probability；Using constellation state probability calculation formula, single star availability is obtained；According to single star availability and service continuity, the conditional probability of visible satellite interruption is obtained, single star signal continuity is further obtained, obtains the Mean Time Between Replacement of short-term unplanned interruption；According to spacing wave fault-free false dismissal probability, the overproof false dismissal probability of spacing wave and service integrity risk, single star signal integrity risk is obtained, mean time between interruptions caused by Mean Time Between Interrupts caused by soft fault and hard fault is further obtained；Further obtain short-term unplanned interruption Mean Time To Recovery；According to short-term unplanned Mean Time Between Interrupts and Mean Time To Recovery, design of interruption is carried out to satellite navigation system.
Description
Technical field
The invention belongs to satellite navigation system reliability fields, and in particular to a kind of satellite navigation system availability continuity
Integrity Allocation method.
Background technique
Satellite navigation system is important space fundamental facilities, directly provides round-the-clock, round-the-clock, high-precision for user
Service.Highly reliable satellite navigation system is built, core is to realize the system services such as good availability, continuity and integrity
Performance.Interruption is the important indicator of direct reflection satellite-signal health status, can build for system development and provide reliability foundation
And guidance.
China's Beidou satellite navigation system has opened operation, and discloses system to the whole world and disclose service performance index,
Space segment is made of three kinds of not co-orbital satellites, belongs to isomery COMPLEX MIXED constellation.GPS of America satellite navigation system uses horse
Er Kefu chain carries out approve- useful index distribution, and this method is suitable for the constellation that homogeneous tracks are constituted, and it is mixed not to be able to satisfy isomery complexity
The Distribution Indexes demand of constellation is closed, and cannot achieve the distribution to continuity, integrity index.How by Beidou satellite navigation system
Service availability, continuity, integrity index decomposition unite as all kinds of interruption indexs, constraint and guiding space section and ground segment are respectively
Reliability design of uniting works, and needs a kind of Allocation method that can be operated, easily realize.
Summary of the invention
In view of this, the present invention provides a kind of satellite navigation system availability continuity integrity Allocation method,
Can be realized system disclosing the availability, continuity, integrity index decomposition of service is all kinds of interruption indexs, thus in satellite
In the design of interruption of navigation system, the design objective of short-term unplanned interruption is provided, which can satisfy system and disclose service
Availability, continuity, integrity index request, further increase the reliability of satellite navigation system.
In order to achieve the above object, the technical solution of the present invention is as follows: a kind of satellite navigation system availability continuity is intact
Property Allocation method, includes the following steps:
Step 1, according to the constellation value α of satellite navigation systemk,nWith service availability Aservice, utilize service availability meter
Formula is calculated, constellation state probability is obtained；
Step 2 obtains single star availability p using constellation state probability calculation formula；
Step 3, according to single star availability p and service continuity Cservice, using service continuity calculation formula,
Obtain the conditional probability of visible satellite interruption；
Step 4, the conditional probability calculation formula interrupted using visible satellite, obtains single star signal continuity c；
Step 5 utilizes being averaged for single star signal continuity c and short-term unplanned interruption according to single star signal continuity c
Interval time MTBOstuRelational expression, MTBO is calculatedstu, realize service continuity CserviceIt distributes to short-term unplanned
The Mean Time Between Replacement MTBO of interruptionstu；
Step 6, according to spacing wave fault-free false dismissal probability Pffmd, the overproof false dismissal probability P of spacing wavemdWith serviced
Good property risk Iservice, using service integrity Risk Calculation formula, obtain spacing wave integrity risk Psis, further according to can
Depending on satellite number N, spacing wave integrity risk P is utilizedsisWith single star signal integrity risk PsatBetween relational expression, calculate
To Psat；
Step 7, according to single star signal false dismissal probability PmdsatWith single star signal integrity risk Psat, single star letter is calculated
Number soft fault probability Psoft, further obtain Mean Time Between Interrupts MTBO caused by soft faultsoft, realization will service intact
Property risk IserviceIt distributes to mean time between interruptions MTBO caused by soft faultsoft；Further according to the short-term non-meter in step 5
Draw Mean Time Between Interrupts MTBOstu, further obtain mean time between interruptions MTBO caused by hard faulthard, realizing will
Service integrity risk IserviceIt distributes to mean time between interruptions MTBO caused by hard faulthard；
Step 8, according to single star availability p in step 2, in prolonged stoppage Mean Time Between Replacement MTBOltWith it is extensive
Multiple time MTTRlt, short-term plan Mean Time Between Interrupts MTBOstsWith recovery time MTTRsts, it is short-term it is unplanned interrupt it is average
Interval time MTBOstuOn the basis of, obtain short-term unplanned interruption Mean Time To Recovery MTTRstu, realize service availability
AserviceIt distributes to short-term unplanned interruption Mean Time To Recovery MTTRstu；
Step 9, according to short-term unplanned Mean Time Between Interrupts MTBOstuWith Mean Time To Recovery MTTRstu, to defending
Star navigation system carries out design of interruption, so that it is MTBO that the short-term unplanned interruption of satellite navigation system, which meets Mean Time Between Replacement,stu
It is MTTR with Mean Time To Recoverystu。
Further, in step 1, according to the constellation value α of satellite navigation systemk,nWith service availability Aservice, wherein M
For the population of satellite in constellation,To have all number of combinations of k satellite failure, α in N satellitek,nFor n combined situation
Under constellation value, utilize service availability calculation formulaObtain the constellation under n combined situation
State probability Pk,n。
Further, in step 2, constellation state probability calculation formula isIts
Middle pn,mThe availability of lower the m operational satellites is combined for n, m takes from 1 to (N-k), (1-pn,i) it is that n combines
The probability of malfunction of specific i-th fault satellites down, the value of i are that N-k+1 takes to N, total k kind situation.
Further, in step 3, according to single star availability p and service continuity Cservice, continuous using service
Property calculation formulaObtain the conditional probability Q of visible satellite interruptioni’,x, it is initial
Stage have total i ' satellite it is visible on the basis of, within the preset stipulated time again have x satellite interrupt conditional probability, wherein
N is the visual satellite number of the user of satellite navigation system, and m is the minimum satellite number for meeting service continuity, Pi’For satellite navigation
The state probability of visual i ' the satellite of the user of system,To advise on the basis of total i ' satellite is visible
The mean availability for thering is x satellite to interrupt again in fixing time,For the number of combinations for thering is x satellite to interrupt in i ' satellite；An
(l, t) is the service availability at observation moment t moment, the position observation position l, and wherein bool () is Boolean function.
Further, in step 4, it is seen that satellite interrupt conditional probability calculation formula be
Wherein c is single star signal continuity.
Further, in step 5, according to single star signal continuity c, formula c=exp (- T/MTBO is utilizedstu), it calculates
Obtain the Mean Time Between Replacement MTBO of short-term unplanned interruptionstu, wherein T is single star signal successional unit time.
Further, in step 6, according to spacing wave fault-free false dismissal probability Pffmd, the overproof false dismissal probability of spacing wave
PmdWith service integrity risk Iservice, utilize service integrity Risk Calculation formula Iservice=Pffmd(1-Psis)+PmdPsis,
Obtain spacing wave integrity risk Psis, further according to satellites in view number N, utilize formulaIt is complete to obtain single star signal
Good property risk Psat；
Wherein, PffmdPosition error occurs when for spacing wave fault-free to transfinite but the probability of receiver false dismissal, PmdFor sky
Between signal fault when position error occur transfinite and the probability of receiver false dismissal, single star signal integrity risk PsatIt describes and defends
The trusting degree of star signal offer information correctness.
Further, step 7 comprises the following specific steps that:
Step 7.1, according to single star signal false dismissal probability PmdsatWith single star signal integrity risk Psat, utilize formula Psat=
PmdsatPsoft, obtain single star signal soft fault probability Psoft, wherein PmdsatExceed integrity alarm limit value for satellite-signal precision
And the probability of alarm is issued not in time；
Step 7.2, time T and single star signal soft fault probability P according to the rulessoft, utilize
Obtain Mean Time Between Interrupts MTBO caused by soft faultsoft, wherein T is the unit time of single star signal integrity risk,
Soft fault refers to that satellite-signal broadcasts normal but signal accuracy and is unsatisfactory for requirement；
Step 7.3, according to the short-term unplanned Mean Time Between Interrupts MTBO in step 5stu, cause in conjunction with soft fault
Mean Time Between Interrupts MTBOsoft, utilize formula 1/MTBOstu=1/MTBOsoft+1/MTBOhard, obtain short-term hard fault
Caused mean time between interruptions MTBOhard, wherein hard fault refers to that satellite-signal broadcasts failure.
Further, step 8 comprises the following specific steps that:
Step 8.1, according to single star availability formula p=MTBO/ (MTBO+MTTR), wherein MTBO and MTTR is to defend
The Mean Time Between Replacement and recovery time that star interrupts, and satellite interruption is divided into prolonged stoppage, short-term plan is interrupted, short-term non-meter
It draws and interrupts three classes, according to following calculation formula:
Wherein, MTBOltAnd MTTRltMean Time Between Replacement and recovery time for prolonged stoppage, MTBOstsAnd MTTRstsFor
The Mean Time Between Replacement and recovery time that short-term plan is interrupted, MTTRstuFor the Mean Time To Recovery of short-term unplanned interruption.
Step 8.2, the Mean Time Between Replacement MTBO of prolonged stoppage is obtained according to satellite design lifetimelt, according to satellite launch,
Orbital detection and joint debugging time are restored time MTTRlt；Short-term plan interruption is obtained according to the frequency of satellite attended operation
Mean Time Between Replacement MTBOsts, it is restored time MTTR according to each time for executing satellite attended operationsts；
Step 8.3, according to single star availability p in step 2, in prolonged stoppage Mean Time Between Replacement MTBOltWith it is extensive
Multiple time MTTRlt, short-term plan Mean Time Between Interrupts MTBOstsWith recovery time MTTRsts, it is short-term it is unplanned interrupt it is average
Interval time MTBOstuOn the basis of, obtain short-term unplanned interruption Mean Time To Recovery MTTRstu, realize service availability
AserviceIt distributes to short-term unplanned Mean Time Between Interrupts MTBOstuWith Mean Time To Recovery MTTRstu。
The utility model has the advantages that
(1) present invention solves the top level services performance indicator such as satellite navigation system service availability, continuity, integrity
The problem of distribution, refers to service availability, continuity, integrity index decomposition for single star availability, continuity, integrity
Mark, and then it is decomposed into prolonged stoppage, short-term plan interruption, short-term unplanned interruption three classes interruption index, and according to short-term non-meter
Draw Mean Time Between Interrupts MTBOstuWith Mean Time To Recovery MTTRstu, design of interruption is carried out to satellite navigation system, that is, is defended
It is MTBO that the short-term unplanned interruption of star navigation system, which meets Mean Time Between Replacement,stuIt is MTTR with Mean Time To Recoverystu, it is this in
Disconnected design method can satisfy the index request that system discloses the availability of service, continuity, integrity, further increase satellite
The reliability of navigation system.
(2) this method is utilized, the availability, continuity, integrity Distribution Indexes of satellite navigation system not only may be implemented,
And for other stars integration star base system (such as communication satellite constellation system, remote sensing constellation systems, early warning constellation systems)
Availability, continuity, integrity Distribution Indexes also have important application value.
Detailed description of the invention
Fig. 1 is flow chart of the present invention.
Specific embodiment
The present invention will now be described in detail with reference to the accompanying drawings and examples.
Embodiment 1, the present invention provides a kind of satellite navigation system availability continuity integrity Allocation methods, will
Service availability, continuity, integrity index decomposition are single star availability, continuity, integrity index, and then are decomposed into
Prolonged stoppage, short-term plan are interrupted, short-term unplanned interruption three classes interrupt index, can provide one kind for service performance Distribution Indexes
The effective ways can operate, easily realized.Detailed process is as shown in Figure 1, include the following steps:
Step 1, constellation state probability calculates.
According to the constellation value α of satellite navigation systemk,nWith service availability Aservice, wherein M is that the satellite in constellation is total
Number,To have all number of combinations of k satellite failure, α in N satellitek,nFor the constellation value under n combined situation, clothes are utilized
Business availability calculations formulaObtain the constellation state probability P under n combined situationk,n。
Wherein, constellation value αk,nThe inherent characteristic of constellation, reflect constellation geometrical property and continuous visibility, refer to
Geometric dilution of precision (DOP) is less than threshold value DOP in the area of coveragemaxRegion account for entire service area area percentage full-time
Average value in section, calculation formula are as follows:
Wherein, Δ T is total simulation time, t0For initial time, L is grid points total number, DOPt,iFor i-th of grid points
In the DOP value of t moment, Area is the coverage gross area, areaiFor the area of i-th of grid points.
Step 2, single star availability calculates.
Constellation state probability calculation formula isWherein pn,mFor under n combination
The availability of the m operational satellites, m take from 1 to (N-k), (1-pn,i) it is that specific i-th failure is defended under n combination
The probability of malfunction of star, the value of i are that N-k+1 takes to N, total k kind situation.
Wherein, constellation state probability, which refers in constellation, there is the not available probability of certain amount satellite-signal, is believed by single star
Number availability impact, and single star availability refers to that the satellite in constellation on regulation orbital position provides healthy space signal
Probability is interrupted to prolonged stoppage, short-term plan, unplanned interruption is related in short term.
Step 3, it is seen that the conditional probability that satellite interrupts calculates.
According to single star availability p and service continuity Cservice, utilize service continuity calculation formulaObtain the conditional probability Q of visible satellite interruptioni’,x, have total i ' for the initial stage
On the basis of satellite is visible, the conditional probability for having x satellite to interrupt again within the preset stipulated time, wherein N leads for satellite
The visual satellite number of the user of boat system, m are the minimum satellite number for meeting service continuity, Pi’For the user of satellite navigation system
The state probability of visual i ' satellite,For on the basis of total i ' satellite is visible, at the appointed time again
The mean availability for thering is x satellite to interrupt,For the number of combinations for thering is x satellite to interrupt in i ' satellite；An(l, t) is to see
Moment t moment is surveyed, the service availability at the position observation position l, wherein bool () is Boolean function.
Wherein, state probability PiIt is calculated by single star availability, user's satellites in view number and height angle of cut-off phase
It closes, is obtained by constellation visibility simulation analysis, minimum satellite number is required to determine by service continuity.
Step 4, single star signal continuity calculates.
The conditional probability calculation formula interrupted using visible satelliteIt is continuous to obtain single star signal
Property c, refer to a healthy open service clearance signal can the stipulated time section in not occur it is unplanned interruption and continue working
Probability.Single star signal continuity and short-term unplanned interruption are closely related.
Step 5, the Mean Time Between Replacement of short-term unplanned interruption calculates.
According to single star signal continuity c, formula c=exp (- T/MTBO is utilizedstu), short-term unplanned interruption is calculated
Mean Time Between Replacement MTBOstu, middle or short term it is unplanned interruption refer to do not notify in advance in the case where cause signal interruption
Event, T is single star signal successional unit time, is generally taken 1 hour.
Step 6, single star signal integrity Risk Calculation.
According to spacing wave fault-free false dismissal probability Pffmd, the overproof false dismissal probability P of spacing wavemdWith service integrity risk
Iservice, utilize service integrity Risk Calculation formula Iservice=Pffmd(1-Psis)+PmdPsis, obtain spacing wave integrity
Risk Psis, further according to satellites in view number N, utilize formulaObtain single star signal integrity risk Psat。
Wherein, PffmdPosition error occurs when for spacing wave fault-free to transfinite but the probability of receiver false dismissal, PmdFor sky
Between signal fault when position error occur transfinite and the probability of receiver false dismissal, single star signal integrity risk PsatIt describes and defends
The trusting degree of star signal offer information correctness.
Step 7, mean time between interruptions caused by soft or hard failure calculates.
Step 7.1, according to single star signal false dismissal probability PmdsatWith single star signal integrity risk Psat, utilize formula Psat=
PmdsatPsoft, obtain single star signal soft fault probability Psoft, wherein PmdsatExceed integrity alarm limit value for satellite-signal precision
And the probability of alarm is issued not in time.
Step 7.2, time T (generally taking 1 hour) and single star signal soft fault probability P according to the rulessoft, utilize Psoft=
exp(-T/MTBOsoft), obtain Mean Time Between Interrupts MTBO caused by soft faultsoft, wherein T is single star signal integrity
The unit time of risk, soft fault refer to that satellite-signal broadcasts normal but signal accuracy and is unsatisfactory for requirement.
Step 7.3, according to the short-term unplanned Mean Time Between Interrupts MTBO in step 5stu, cause in conjunction with soft fault
Mean Time Between Interrupts MTBOsoft, utilize formula 1/MTBOstu=1/MTBOsoft+1/MTBOhard, obtain short-term hard fault
Caused mean time between interruptions MTBOhard, wherein hard fault refers to that satellite-signal broadcasts failure.
Step 8, short-term unplanned interruption Mean Time To Recovery calculates.
Step 8.1, according to single star availability formula p=MTBO/ (MTBO+MTTR), wherein MTBO and MTTR is to defend
The Mean Time Between Replacement and recovery time that star interrupts, and satellite interruption is divided into prolonged stoppage, short-term plan is interrupted, short-term non-meter
It draws and interrupts three classes, according to following calculation formula:
Wherein, MTBOltAnd MTTRltMean Time Between Replacement and recovery time for prolonged stoppage, MTBOstsAnd MTTRstsFor
The Mean Time Between Replacement and recovery time that short-term plan is interrupted, MTTRstuFor the Mean Time To Recovery of short-term unplanned interruption.
By in single star availability formula MTBO and MTTR be replaced, establish in single star availability and three classes
Disconnected relationship.
Step 8.2, the Mean Time Between Replacement MTBO of prolonged stoppage is obtained according to satellite design lifetimelt, according to satellite launch,
Orbital detection and joint debugging time are restored time MTTRlt；Short-term plan interruption is obtained according to the frequency of satellite attended operation
Mean Time Between Replacement MTBOsts, it is restored time MTTR according to each time for executing satellite attended operationsts。
Step 8.3, according to single star availability p in step 2, in prolonged stoppage Mean Time Between Replacement MTBOltWith it is extensive
Multiple time MTTRlt, short-term plan Mean Time Between Interrupts MTBOstsWith recovery time MTTRsts, it is short-term it is unplanned interrupt it is average
Interval time MTBOstuOn the basis of, obtain short-term unplanned interruption Mean Time To Recovery MTTRstu, realize service availability
AserviceIt distributes to short-term unplanned Mean Time Between Interrupts MTBOstuWith Mean Time To Recovery MTTRstu。
Step 9, according to short-term unplanned Mean Time Between Interrupts MTBOstuWith Mean Time To Recovery MTTRstu, to defending
Star navigation system carries out design of interruption, i.e., it is MTBO that the short-term unplanned interruption of satellite navigation system, which meets Mean Time Between Replacement,stuWith
Mean Time To Recovery is MTTRstu。
Embodiment is named to further illustrate method provided by the present invention.
Such as: to Mr. Yu's satellite navigation system, which is made of 24 satellites, and constellation is configured as Walker constellation, clothes
Availability of being engaged in is 0.9, service continuity 0.99, and service integrity risk is 0.999, carries out above-mentioned steps 1~4, obtains list
Star availability is 0.957, and single star signal continuity is 0.9998/ hour；Above-mentioned steps 5 are carried out, are obtained short-term unplanned
The Mean Time Between Replacement of interruption is 4380 hours；Step 6 is carried out, obtaining single star signal integrity risk is 1E-5/ hours；It carries out
Above-mentioned steps 7, obtaining Mean Time Between Interrupts caused by soft fault is 1E5 hours, interrupts equispaced caused by hard fault
Time is 4580 hours；Above-mentioned steps 8 are carried out, the Mean Time Between Replacement for obtaining prolonged stoppage is 52560 hours, and recovery time is
1750 hours, the Mean Time Between Replacement that short-term plan is interrupted was 4380 hours, and recovery time is 12 hours, short-term unplanned interruption
Mean Time To Recovery be 36 hours.
To sum up, the above is merely preferred embodiments of the present invention, it is not intended to limit the scope of the present invention.It is all
Within the spirit and principles in the present invention, any modification, equivalent replacement, improvement and so on should be included in protection of the invention
Within the scope of.
Claims (9)
1. a kind of satellite navigation system availability continuity integrity Allocation method, which comprises the steps of:
Step 1, according to the constellation value α of satellite navigation systemk,nWith service availability Aservice, calculated using service availability public
Formula obtains constellation state probability；
Step 2 obtains single star availability p using constellation state probability calculation formula；
Step 3, according to single star availability p and service continuity Cservice, using service continuity calculation formula, obtaining can
See the conditional probability that satellite interrupts；
Step 4, the conditional probability calculation formula interrupted using visible satellite, obtains single star signal continuity c；
Step 5 utilizes the equispaced of single star signal continuity c and short-term unplanned interruption according to single star signal continuity c
Time MTBOstuRelational expression, MTBO is calculatedstu, realize service continuity CserviceIt distributes to short-term unplanned interruption
Mean Time Between Replacement MTBOstu；
Step 6, according to spacing wave fault-free false dismissal probability Pffmd, the overproof false dismissal probability P of spacing wavemdWith service integrity
Risk Iservice, using service integrity Risk Calculation formula, obtain spacing wave integrity risk Psis, further according to visually defending
Star number N utilizes spacing wave integrity risk PsisWith single star signal integrity risk PsatBetween relational expression, be calculated
Psat；
Step 7, according to single star signal false dismissal probability PmdsatWith single star signal integrity risk Psat, it is soft that single star signal is calculated
Probability of malfunction Psoft, further obtain Mean Time Between Interrupts MTBO caused by soft faultsoft, integrity wind will be serviced by realizing
Dangerous IserviceIt distributes to mean time between interruptions MTBO caused by soft faultsoft；Further according in step 5 it is short-term it is unplanned in
Disconnected Mean Time Between Replacement MTBOstu, further obtain mean time between interruptions MTBO caused by hard faulthard, realizing will service
Integrity risk IserviceIt distributes to mean time between interruptions MTBO caused by hard faulthard；
Step 8, according to single star availability p in step 2, in prolonged stoppage Mean Time Between Replacement MTBOltWhen with restoring
Between MTTRlt, short-term plan Mean Time Between Interrupts MTBOstsWith recovery time MTTRsts, short-term unplanned interruption equispaced
Time MTBOstuOn the basis of, obtain short-term unplanned interruption Mean Time To Recovery MTTRstu, realize service availability Aservice
It distributes to short-term unplanned interruption Mean Time To Recovery MTTRstu；
Step 9, according to short-term unplanned Mean Time Between Interrupts MTBOstuWith Mean Time To Recovery MTTRstu, satellite is led
Boat system carries out design of interruption, so that it is MTBO that the short-term unplanned interruption of satellite navigation system, which meets Mean Time Between Replacement,stuPeace
Equal recovery time is MTTRstu。
2. a kind of satellite navigation system availability continuity integrity Allocation method as described in claim 1, feature
It is, in the step 1, according to the constellation value α of satellite navigation systemk,nWith service availability Aservice, wherein M is in constellation
The population of satellite,To have all number of combinations of k satellite failure, α in N satellitek,nFor the constellation under n combined situation
Value, utilizes service availability calculation formulaObtain the constellation state probability under n combined situation
Pk,n。
3. a kind of satellite navigation system availability continuity integrity Allocation method as claimed in claim 2, feature
It is, in the step 2, constellation state probability calculation formula isWherein pn,mIt is
N kind combines the availability of lower the m operational satellites, and m takes from 1 to (N-k), (1-pn,i) it is specific the under n combination
The probability of malfunction of i fault satellites, the value of i are that N-k+1 takes to N, total k kind situation.
4. a kind of satellite navigation system availability continuity integrity Allocation method as claimed in claim 3, feature
It is, in the step 3, according to single star availability p and service continuity Cservice, calculated using service continuity public
FormulaObtain the conditional probability Q of visible satellite interruptioni’,x, have altogether for the initial stage
On the basis of i ' satellite is visible, the conditional probability for having x satellite to interrupt again within the preset stipulated time, wherein N is satellite
The visual satellite number of the user of navigation system, m are the minimum satellite number for meeting service continuity, Pi’For the use of satellite navigation system
The state probability of visual i ' the satellite in family,For on the basis of total i ' satellite is visible, at the appointed time
The mean availability for thering is x satellite to interrupt again,For the number of combinations for thering is x satellite to interrupt in i ' satellite；An(l, t) be
Moment t moment is observed, the service availability at the position observation position l, wherein bool () is Boolean function.
5. a kind of satellite navigation system availability continuity integrity Allocation method as claimed in claim 4, feature
Be, in the step 4, it is seen that satellite interrupt conditional probability calculation formula beWherein c is single
Star signal continuity.
6. a kind of satellite navigation system availability continuity integrity Allocation method as claimed in claim 5, feature
It is, in the step 5, according to single star signal continuity c,
Utilize formula c=exp (- T/MTBOstu), the Mean Time Between Replacement MTBO of short-term unplanned interruption is calculatedstu, wherein
T is single star signal successional unit time.
7. a kind of satellite navigation system availability continuity integrity Allocation method as claimed in claim 6, feature
It is, in the step 6, according to spacing wave fault-free false dismissal probability Pffmd, the overproof false dismissal probability P of spacing wavemdAnd service
Integrity risk Iservice, utilize service integrity Risk Calculation formula Iservice=Pffmd(1-Psis)+PmdPsis, obtain space
Signal integrity risk Psis, further according to satellites in view number N, utilize formulaObtain single star signal integrity risk
Psat；
Wherein, PffmdPosition error occurs when for spacing wave fault-free to transfinite but the probability of receiver false dismissal, PmdFor space letter
Position error occurs when number failure to transfinite and the probability of receiver false dismissal, single star signal integrity risk PsatDescribe satellite letter
Number provide information correctness trusting degree.
8. a kind of satellite navigation system availability continuity integrity Allocation method as claimed in claim 7, feature
It is, the step 7 comprises the following specific steps that:
Step 7.1, according to single star signal false dismissal probability PmdsatWith single star signal integrity risk Psat, utilize formula Psat=
PmdsatPsoft, obtain single star signal soft fault probability Psoft, wherein PmdsatExceed integrity alarm limit value for satellite-signal precision
And the probability of alarm is issued not in time；
Step 7.2, time T and single star signal soft fault probability P according to the rulessoft,
Utilize Psoft=exp (- T/MTBOsoft), obtain Mean Time Between Interrupts MTBO caused by soft faultsoft, wherein T be
The unit time of single star signal integrity risk, soft fault refer to that satellite-signal is broadcast normal but signal accuracy and is unsatisfactory for using wanting
It asks；
Step 7.3, according to the short-term unplanned Mean Time Between Interrupts MTBO in step 5stu, in conjunction with caused by soft fault in
Disconnected Mean Time Between Replacement MTBOsoft, utilize formula 1/MTBOstu=1/MTBOsoft+1/MTBOhard, obtain short-term hard fault and cause
Mean time between interruptions MTBOhard, wherein hard fault refers to that satellite-signal broadcasts failure.
9. a kind of satellite navigation system availability continuity integrity Allocation method as claimed in claim 8, feature
It is, the step 8 comprises the following specific steps that:
Step 8.1, according to single star availability formula p=MTBO/ (MTBO+MTTR), wherein MTBO and MTTR is in satellite
Disconnected Mean Time Between Replacement and recovery time, and satellite interruption be divided into prolonged stoppage, short-term plan interrupt, in short term it is unplanned in
Disconnected three classes, according to following calculation formula:
Wherein, MTBOltAnd MTTRltMean Time Between Replacement and recovery time for prolonged stoppage, MTBOstsAnd MTTRstsIt is short-term
Plan the Mean Time Between Replacement interrupted and recovery time, MTTRstuFor the Mean Time To Recovery of short-term unplanned interruption；
Step 8.2, the Mean Time Between Replacement MTBO of prolonged stoppage is obtained according to satellite design lifetimelt, according to satellite launch, in-orbit
Test and joint debugging time are restored time MTTRlt；Being averaged for short-term plan interruption is obtained according to the frequency of satellite attended operation
Interval time MTBOsts, it is restored time MTTR according to each time for executing satellite attended operationsts；
Step 8.3, according to single star availability p in step 2, in prolonged stoppage Mean Time Between Replacement MTBOltWhen with restoring
Between MTTRlt, short-term plan Mean Time Between Interrupts MTBOstsWith recovery time MTTRsts, short-term unplanned interruption equispaced
Time MTBOstuOn the basis of, obtain short-term unplanned interruption Mean Time To Recovery MTTRstu,
It realizes service availability AserviceIt distributes to short-term unplanned Mean Time Between Interrupts MTBOstuAnd Mean Time To Recovery
MTTRstu。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611027438.4A CN106610495B (en) | 2016-11-18 | 2016-11-18 | A kind of satellite navigation system availability continuity integrity Allocation method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611027438.4A CN106610495B (en) | 2016-11-18 | 2016-11-18 | A kind of satellite navigation system availability continuity integrity Allocation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106610495A CN106610495A (en) | 2017-05-03 |
| CN106610495B true CN106610495B (en) | 2019-01-08 |
Family
ID=58636174
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611027438.4A Expired - Fee Related CN106610495B (en) | 2016-11-18 | 2016-11-18 | A kind of satellite navigation system availability continuity integrity Allocation method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106610495B (en) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108334416B (en) * | 2018-01-08 | 2021-05-04 | 千寻位置网络有限公司 | Sampling acquisition automatic verification method for broadcast data |
| CN108897021B (en) * | 2018-04-20 | 2022-10-14 | 中国航天标准化研究所 | Method for determining availability of single orbit position of navigation satellite |
| CN111060933B (en) * | 2019-12-11 | 2022-03-04 | 中国航天标准化研究所 | Global continuous monitoring ground tracking station availability distribution method |
| CN111310337B (en) * | 2020-02-18 | 2023-08-29 | 北京空间飞行器总体设计部 | Analysis method for on-orbit interruption frequency of navigation satellite signal |
| CN111323791B (en) * | 2020-04-10 | 2022-04-05 | 航天恒星科技有限公司 | Quality evaluation method and device of navigation enhancement product |
| CN111967210B (en) * | 2020-07-03 | 2023-06-06 | 北京空间飞行器总体设计部 | PHM design method for navigation downlink signal function interruption |
| CN112269194B (en) * | 2020-09-17 | 2023-06-30 | 中国航天标准化研究所 | Navigation satellite system continuity risk assessment method |
| CN112417589B (en) * | 2020-11-05 | 2023-11-10 | 航天东方红卫星有限公司 | Satellite steady-state availability computing method |
| CN113627717B (en) * | 2021-06-30 | 2024-03-29 | 中国空间技术研究院 | Index distribution method for communication satellite service interruption caused by space single event |
| CN114609650B (en) * | 2022-03-09 | 2022-12-16 | 中国人民解放军92728部队 | Integrity testing method based on Beidou full-chain fault excitation |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030184275A1 (en) * | 2002-03-26 | 2003-10-02 | Public Service Company Of New Mexico | System, method, field unit, reference unit and computer program product for phase tracking of electrical conductors |
| CN104215974A (en) * | 2014-09-15 | 2014-12-17 | 中国航天标准化研究所 | Integrity detection availability determining method of satellite navigation system |
| CN105738922A (en) * | 2015-09-21 | 2016-07-06 | 北京圣涛平试验工程技术研究院有限责任公司 | Service reliability analysis method and system of navigation satellite constellation system |
| CN105929414A (en) * | 2016-06-14 | 2016-09-07 | 中国航天标准化研究所 | Navigation signal integrity assessment method |
| CN106125100A (en) * | 2016-06-14 | 2016-11-16 | 中国航天标准化研究所 | A kind of spacing wave seriality appraisal procedure |
-
2016
- 2016-11-18 CN CN201611027438.4A patent/CN106610495B/en not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20030184275A1 (en) * | 2002-03-26 | 2003-10-02 | Public Service Company Of New Mexico | System, method, field unit, reference unit and computer program product for phase tracking of electrical conductors |
| CN104215974A (en) * | 2014-09-15 | 2014-12-17 | 中国航天标准化研究所 | Integrity detection availability determining method of satellite navigation system |
| CN105738922A (en) * | 2015-09-21 | 2016-07-06 | 北京圣涛平试验工程技术研究院有限责任公司 | Service reliability analysis method and system of navigation satellite constellation system |
| CN105929414A (en) * | 2016-06-14 | 2016-09-07 | 中国航天标准化研究所 | Navigation signal integrity assessment method |
| CN106125100A (en) * | 2016-06-14 | 2016-11-16 | 中国航天标准化研究所 | A kind of spacing wave seriality appraisal procedure |
Non-Patent Citations (5)
| Title |
|---|
| Feasibility of Achieving GBAS Category II/III;Rick Cassell et.al;《ION 63rd Annual Meeting》;20071231;全文 |
| The Galileo Integrity Concept;Veit Oehler et.al;《ION GNSS 17th International Technical Meeting》;20071231;全文 |
| 利用监测网中断分析导航***完好性监测;杨卓鹏等;《导航定位学报》;20150630;第3卷(第2期);第38-48页 |
| 卫星导航***性能监测及评估方法研究;李作虎;《中国博士学位论文全文数据库基础科学辑》;20130615;A008-12 |
| 卫星导航***星座可用性分析;郑恒等;《航天控制》;20110630;第29卷(第3期);第87-94页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106610495A (en) | 2017-05-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106610495B (en) | A kind of satellite navigation system availability continuity integrity Allocation method | |
| CN101776762B (en) | Completeness monitoring method, device and system based on multi-foundation enhancement system | |
| US20020120391A1 (en) | Method and system for displaying target vehicle position information | |
| CN101893713B (en) | Device and method for monitoring release information performance of local area augmentation system | |
| CN104215974A (en) | Integrity detection availability determining method of satellite navigation system | |
| Team | Global positioning system (gps) standard positioning service (sps) performance analysis report | |
| CN104331593B (en) | For ground prediction aircraft along the device and method of the feature of the positioning in path | |
| CN105866802A (en) | Method and device for monitoring delay integrity of troposphere | |
| CN107591032A (en) | For calculating required navigational foreseeable apparatus and method | |
| Lee et al. | GPS and Galileo with RAIM or WAAS for Vertically Guided Approaches | |
| Lee | RAIM Availability for GPS Augmented with Barometric Altimeter Aiding and Clock Coasting 1 | |
| Zalewski | GNSS integrity concepts for maritime users | |
| Krzykowska et al. | Signal monitoring as a part of maintenance of navigation support system in civil aviation | |
| Shi et al. | Structural resilience in sewer reconstruction: from theory to practice | |
| Filip et al. | Safety concept of railway signalling based on Galileo Safety-of-Life Service | |
| Savitz et al. | Enhancing US Coast Guard Metrics | |
| Martineau | Performance of receiver autonomous integrity monitoring (RAIM) for vertically guided approaches | |
| Shi et al. | Reviews of resilience theories and mathematical generalization | |
| Konno | Design of an aircraft landing system using dual-frequency GNSS | |
| Lazaro | EUMETSAT LEO Conjunctions Events Operations | |
| Hamzah et al. | Prototyping of Mobile Application and Web Monitoring for Small Vessel Tracking System | |
| Song et al. | GBAS Integrity Performance Evaluation Based on the Mixed g-and-h Distribution | |
| Mink | Performance of receiver autonomous integrity monitoring (RAIM) for maritime operations | |
| Klepsvik et al. | A Critical Look at the IMO Requirements for GNSS | |
| García González | Development of a 3-dimensional mathematical collision risk model based on recorded aircraft trajectories to estimate the safety level in high density en-route airspaces |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190108 Termination date: 20191118 |