阅读说明：本技术 一种同步卫星外测数据的质量监控方法、装置和电子设备 (Method and device for monitoring quality of geostationary satellite external measurement data and electronic equipment ) 是由 张俊丽 董开封 冯卫东 姚凡凡 邢东旭 蒋振伟 张文雅 李峰 吴琛 李昂 崔晓阳 于 2021-03-05 设计创作，主要内容包括：本发明提供了一种同步卫星外测数据的质量监控方法、装置和电子设备,涉及通信的技术领域,在获取到实时测距数据后,首先利用实时测距数据和日常轨道精密定轨与预报输出的卫星星历文件确定出实时测距数据对应的理论测距数据,然后再结合当前气象数据和测距系统误差确定出实时测距数据的随机误差。本发明将轨道动力学模型和观测模型等先验信息用于同步卫星外测数据的实时处理,提高了外测数据质量异常监测的敏感度,基于实时测距数据的随机误差变化可实现对同步卫星外测数据的实时监测评估,为及时掌握外测数据质量,调整定轨策略提供了数据支撑。(The invention provides a quality monitoring method and device for synchronous satellite external measurement data and electronic equipment, and relates to the technical field of communication. The invention uses the prior information such as the orbit dynamics model, the observation model and the like for the real-time processing of the synchronous satellite external data, improves the sensitivity of the abnormal monitoring of the external data quality, can realize the real-time monitoring and evaluation of the synchronous satellite external data based on the random error change of the real-time ranging data, and provides data support for timely mastering the external data quality and adjusting the orbit determination strategy.)

1. A quality monitoring method for synchronous satellite external measurement data is characterized by comprising the following steps:

receiving current meteorological data and real-time ranging data sent by a ground ranging master station;

acquiring satellite ephemeris files and ranging system errors which are output by daily orbit precise orbit determination and prediction;

determining theoretical ranging data corresponding to the real-time ranging data based on the real-time ranging data and the satellite ephemeris file;

determining a random error of the real-time ranging data based on the current meteorological data, the ranging system error, the real-time ranging data, and the theoretical ranging data; and the random error is used for representing the data quality of the real-time ranging data.

2. The method of claim 1, wherein determining theoretical ranging data corresponding to real-time ranging data based on the real-time ranging data and the satellite ephemeris file comprises:

acquiring a ranging time scale of the real-time ranging data;

determining a theoretical ranging reference value corresponding to the real-time ranging data based on the ranging time scale and the satellite ephemeris file;

correcting the ranging time scale based on the theoretical ranging reference value and a preset fuzzy period to obtain the actual transmitting time of the ranging pseudo code signal;

and determining theoretical ranging data corresponding to the real-time ranging data based on the actual transmitting time and the satellite ephemeris file.

3. The method of claim 2, wherein determining theoretical ranging reference values corresponding to the real-time ranging data based on the ranging time stamp and the satellite ephemeris file comprises:

determining a geocentric position vector of a satellite at a corresponding moment of the ranging time scale based on the ranging time scale and the satellite ephemeris file;

if the real-time ranging data are single-pass ranging data or double-pass ranging data, determining a theoretical ranging reference value corresponding to the single-pass ranging data or a theoretical ranging reference value corresponding to the double-pass ranging data based on a geocentric position vector of the satellite at the corresponding moment of the ranging time scale and a geocentric position vector of the ground ranging master station;

and if the real-time ranging data are four-range ranging data, determining a theoretical ranging reference value corresponding to the four-range ranging data based on the geocentric position vector of the satellite at the corresponding moment of the ranging time scale, the geocentric position vector of the ground ranging main station and the geocentric position vector of the ground ranging auxiliary station.

4. The method of claim 2, wherein determining theoretical ranging data corresponding to the real-time ranging data based on the actual transmission time and the satellite ephemeris file comprises:

determining a geocentric position vector of a satellite at the actual transmission time based on the actual transmission time and the satellite ephemeris file;

if the real-time ranging data are single-pass ranging data or double-pass ranging data, determining ranging time delay of the single-pass ranging data or ranging time delay of the double-pass ranging data based on a geocentric position vector of the satellite at the actual transmitting moment and a geocentric position vector of a ground ranging master station;

determining theoretical ranging data corresponding to the one-way ranging data based on the ranging time delay and the light speed of the one-way ranging data; or determining theoretical ranging data corresponding to the two-way ranging data based on the ranging time delay and the light speed of the two-way ranging data;

if the real-time ranging data are four-range ranging data, determining ranging time delay of the four-range ranging data based on the geocentric position vector of the satellite at the actual transmitting moment, the geocentric position vector of the ground ranging main station and the geocentric position vector of the ground ranging auxiliary station;

and determining theoretical ranging data corresponding to the four-range ranging data based on the ranging time delay and the light speed of the four-range ranging data.

5. The method of claim 2, wherein determining the random error of the real-time ranging data based on the current meteorological data, the ranging system error, the real-time ranging data, and the theoretical ranging data comprises:

performing outlier identification and elimination on the current meteorological data and the real-time ranging data to obtain target meteorological data and target ranging data;

determining an atmospheric refraction correction corresponding to the ranging data based on the target meteorological data;

correcting the target ranging data based on the theoretical ranging reference value and a preset fuzzy period to obtain corrected ranging data;

and determining the random error of the real-time ranging data based on the atmospheric refraction correction amount, the ranging system error, the ranging data after correction and the theoretical ranging data.

6. The method of claim 2, wherein the step of modifying the ranging time scale based on the theoretical ranging reference value and a preset ambiguity period to obtain an actual transmission time of the ranging pseudo code signal comprises:

correcting the ranging time scale based on the formula T ═ T-N · Δ T to obtain the actual transmission time of the ranging pseudo code signal; wherein t 'represents the actual transmission time of the ranging pseudo code signal, t represents the corresponding time of the ranging time scale, N represents the ambiguity number, and N ═ ρ'_c/Δρ|，ρ′_cRepresents the theoretical ranging reference value, Δ ρ ═ Δ T × c, Δ T represents the preset blur period, and c represents the speed of light.

7. The method of claim 5, wherein determining an atmospheric refraction correction for ranging data based on the target meteorological data comprises:

equation of utilizationDetermining the atmospheric refraction correction corresponding to the distance measurement data, wherein delta R represents the atmospheric refraction correction, C represents the atmospheric refractive index, E_NWhich represents the elevation angle of a theoretical survey station,which represents the refractive index of the ground,t represents absolute temperature, P represents ground atmospheric pressure, P_eRepresenting the ground water pressure.

8. A quality monitoring device for geostationary satellite telemetry data, comprising:

the receiving module is used for receiving the current meteorological data and the real-time distance measurement data sent by the ground distance measurement master station;

the acquisition module is used for acquiring satellite ephemeris files and ranging system errors which are output by the precision orbit determination and prediction of the daily orbit;

the first determining module is used for determining theoretical ranging data corresponding to the real-time ranging data based on the real-time ranging data and the satellite ephemeris file;

a second determining module for determining a random error of the real-time ranging data based on the current meteorological data, the ranging system error, the real-time ranging data, and the theoretical ranging data; and the random error is used for representing the data quality of the real-time ranging data.

9. An electronic device comprising a memory, a processor, and a computer program stored on the memory and executable on the processor, wherein the processor implements the steps of the method of any of claims 1 to 7 when executing the computer program.

10. A computer-readable medium having non-volatile program code executable by a processor, the program code causing the processor to perform the method of any of claims 1 to 7.

Technical Field

The invention relates to the technical field of communication, in particular to a method and a device for monitoring quality of geostationary satellite external measurement data and electronic equipment.

Background

The orbit measurement and determination are one of the core capabilities of the in-orbit management system of the synchronous satellite and are important bases for ensuring the normal operation and the full and effective exertion of the application efficiency of the synchronous satellite. The orbit measurement is the basis for realizing high-precision orbit determination and prediction, and the precision of the external measurement data processing directly influences the precision of the geostationary satellite orbit determination. However, in the prior art, the overall quality condition of the external measurement data is not monitored by an effective means, the external measurement data processing needs to be manually rechecked by a professional technician, the operation state of the whole process of the geostationary satellite orbit determination processing cannot be rapidly known, the working intensity of the personnel is high, and the personnel is not suitable for the current situation of the follow-up satellite platform in-orbit management.

Disclosure of Invention

The invention aims to provide a quality monitoring method and device for synchronous satellite external measurement data and electronic equipment, which realize real-time monitoring and evaluation of the satellite external measurement data by random error change based on real-time ranging data and provide data support for timely mastering the quality of the external measurement data and adjusting an orbit determination strategy.

In a first aspect, the present invention provides a method for monitoring quality of geostationary satellite external data, including: receiving current meteorological data and real-time ranging data sent by a ground ranging master station; acquiring satellite ephemeris files and ranging system errors which are output by daily orbit precise orbit determination and prediction; determining theoretical ranging data corresponding to the real-time ranging data based on the real-time ranging data and the satellite ephemeris file; determining a random error of the real-time ranging data based on the current meteorological data, the ranging system error, the real-time ranging data, and the theoretical ranging data; and the random error is used for representing the data quality of the real-time ranging data.

In an optional embodiment, determining theoretical ranging data corresponding to real-time ranging data based on the real-time ranging data and the satellite ephemeris file includes: acquiring a ranging time scale of the real-time ranging data; determining a theoretical ranging reference value corresponding to the real-time ranging data based on the ranging time scale and the satellite ephemeris file; correcting the ranging time scale based on the theoretical ranging reference value and a preset fuzzy period to obtain the actual transmitting time of the ranging pseudo code signal; and determining theoretical ranging data corresponding to the real-time ranging data based on the actual transmitting time and the satellite ephemeris file.

In an optional embodiment, determining a theoretical ranging reference value corresponding to the real-time ranging data based on the ranging timestamp and the satellite ephemeris file includes: determining a geocentric position vector of a satellite at a corresponding moment of the ranging time scale based on the ranging time scale and the satellite ephemeris file; if the real-time ranging data are single-pass ranging data or double-pass ranging data, determining a theoretical ranging reference value corresponding to the single-pass ranging data or a theoretical ranging reference value corresponding to the double-pass ranging data based on a geocentric position vector of the satellite at the corresponding moment of the ranging time scale and a geocentric position vector of the ground ranging master station; and if the real-time ranging data are four-range ranging data, determining a theoretical ranging reference value corresponding to the four-range ranging data based on the geocentric position vector of the satellite at the corresponding moment of the ranging time scale, the geocentric position vector of the ground ranging main station and the geocentric position vector of the ground ranging auxiliary station.

In an optional embodiment, determining theoretical ranging data corresponding to the real-time ranging data based on the actual transmission time and the satellite ephemeris file includes: determining a geocentric position vector of a satellite at the actual transmission time based on the actual transmission time and the satellite ephemeris file; if the real-time ranging data are single-pass ranging data or double-pass ranging data, determining ranging time delay of the single-pass ranging data or ranging time delay of the double-pass ranging data based on a geocentric position vector of the satellite at the actual transmitting moment and a geocentric position vector of a ground ranging master station; determining theoretical ranging data corresponding to the one-way ranging data based on the ranging time delay and the light speed of the one-way ranging data; or determining theoretical ranging data corresponding to the two-way ranging data based on the ranging time delay and the light speed of the two-way ranging data; if the real-time ranging data are four-range ranging data, determining ranging time delay of the four-range ranging data based on the geocentric position vector of the satellite at the actual transmitting moment, the geocentric position vector of the ground ranging main station and the geocentric position vector of the ground ranging auxiliary station; and determining theoretical ranging data corresponding to the four-range ranging data based on the ranging time delay and the light speed of the four-range ranging data.

In an alternative embodiment, determining the random error of the real-time ranging data based on the current meteorological data, the ranging system error, the real-time ranging data, and the theoretical ranging data comprises: performing outlier identification and elimination on the current meteorological data and the real-time ranging data to obtain target meteorological data and target ranging data; determining an atmospheric refraction correction corresponding to the ranging data based on the target meteorological data; correcting the target ranging data based on the theoretical ranging reference value and a preset fuzzy period to obtain corrected ranging data; and determining the random error of the real-time ranging data based on the atmospheric refraction correction amount, the ranging system error, the ranging data after correction and the theoretical ranging data.

In an optional embodiment, the modifying the ranging time scale based on the theoretical ranging reference value and a preset fuzzy period to obtain an actual transmission time of the ranging pseudo code signal includes: correcting the ranging time scale based on the formula T ═ T-N · Δ T to obtain the actual transmission time of the ranging pseudo code signal; wherein t 'represents the actual transmission time of the ranging pseudo code signal, t represents the corresponding time of the ranging time scale, N represents the ambiguity number, and N ═ ρ'_c/Δρ|，ρ′_cRepresents the theoretical ranging reference value, Δ ρ ═ Δ T × c, Δ T represents the preset blur period, and c represents the speed of light.

In an optional embodiment, determining an atmospheric refraction correction amount corresponding to the ranging data based on the target meteorological data includes: equation of utilizationDetermining the atmospheric refraction correction corresponding to the distance measurement data, wherein delta R represents the atmospheric refraction correction, C represents the atmospheric refractive index, E_NWhich represents the elevation angle of a theoretical survey station,which represents the refractive index of the ground,t represents absolute temperature, P represents ground atmospheric pressure, P_eRepresenting the ground water pressure.

In a second aspect, the present invention provides a quality monitoring apparatus for geostationary satellite external data, including: the receiving module is used for receiving the current meteorological data and the real-time distance measurement data sent by the ground distance measurement master station; the acquisition module is used for acquiring satellite ephemeris files and ranging system errors which are output by the precision orbit determination and prediction of the daily orbit; the first determining module is used for determining theoretical ranging data corresponding to the real-time ranging data based on the real-time ranging data and the satellite ephemeris file; a second determining module for determining a random error of the real-time ranging data based on the current meteorological data, the ranging system error, the real-time ranging data, and the theoretical ranging data; and the random error is used for representing the data quality of the real-time ranging data.

In a third aspect, the present invention provides an electronic device, comprising a memory and a processor, wherein the memory stores a computer program operable on the processor, and the processor executes the computer program to implement the steps of the method according to any of the foregoing embodiments.

In a fourth aspect, the invention provides a computer readable medium having non-volatile program code executable by a processor, the program code causing the processor to perform the method of any of the preceding embodiments.

The invention provides a quality monitoring method for synchronous satellite external measurement data, which comprises the following steps: receiving current meteorological data and real-time ranging data sent by a ground ranging master station; acquiring satellite ephemeris files and ranging system errors which are output by daily orbit precise orbit determination and prediction; determining theoretical ranging data corresponding to the real-time ranging data based on the real-time ranging data and the satellite ephemeris file; determining a random error of the real-time ranging data based on the current meteorological data, the ranging system error, the real-time ranging data and the theoretical ranging data; the random error is used for representing the data quality of the real-time ranging data.

According to the quality monitoring method for the synchronous satellite external measurement data, after the real-time distance measurement data are obtained, theoretical distance measurement data corresponding to the real-time distance measurement data are determined by using the real-time distance measurement data and a satellite ephemeris file output by daily orbit precise orbit determination and prediction, and then random errors of the real-time distance measurement data are determined by combining the current meteorological data and the distance measurement system errors. According to the method, the prior information such as the orbit dynamics model and the observation model is used for real-time processing of the synchronous satellite external data, the sensitivity of abnormal monitoring of the external data quality is improved, real-time monitoring and evaluation of the satellite external data can be realized based on random error change of real-time ranging data, and data support is provided for timely mastering the external data quality and adjusting the orbit determination strategy.

Drawings

In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.

Fig. 1 is a flowchart of a method for monitoring quality of geostationary satellite external data according to an embodiment of the present invention;

fig. 2 is a flowchart illustrating a process of determining theoretical ranging data corresponding to real-time ranging data based on the real-time ranging data and a satellite ephemeris file according to an embodiment of the present invention;

fig. 3 is a schematic diagram illustrating a pseudo code ranging measurement principle according to an embodiment of the present invention;

fig. 4 is a schematic diagram of two-way distance measurement and four-way distance measurement according to an embodiment of the present invention;

fig. 5 is a quality monitoring interface of multi-station two-way and four-way ranging data according to an embodiment of the present invention;

FIG. 6 is a functional block diagram of an apparatus for monitoring quality of geostationary satellite telemetry data according to an embodiment of the present invention;

fig. 7 is a schematic diagram of an electronic device according to an embodiment of the present invention.

Detailed Description

In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. The components of embodiments of the present invention generally described and illustrated in the figures herein may be arranged and designed in a wide variety of different configurations.

Thus, the following detailed description of the embodiments of the present invention, presented in the figures, is not intended to limit the scope of the invention, as claimed, but is merely representative of selected embodiments of the invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.

Some embodiments of the invention are described in detail below with reference to the accompanying drawings. The embodiments described below and the features of the embodiments can be combined with each other without conflict.

The orbit measurement data comprises one-way distance measurement of a ground distance measurement main station, two-way distance measurement data, four-way distance measurement data of an auxiliary station and meteorological data of each station, and the ground operation and control center uses the orbit observation data to complete the precise orbit determination of the satellite every day. The modification and maintenance of ground station equipment, equipment switching, rainfall and the like can all affect the external measurement data of the synchronous satellite, the orbit determination precision of the satellite can be possibly affected, and the orbit determination strategy is adjusted according to the change of the quality of the external measurement data when a post operator carries out a synchronous satellite orbit determination task every day. At present, the overall quality condition of the external measurement data is not monitored by an effective means, and basically needs to be rechecked by professional technicians, so that the precise orbit determination of the satellite needs excessive manual intervention, the operation state of the whole process of the orbit determination processing of the synchronous satellite cannot be quickly obtained, the working strength of the personnel is high, and the personnel is not suitable for the current situation of the on-orbit management of a subsequent satellite platform. In view of the above, embodiments of the present invention provide a method for quality monitoring of geostationary satellite telemetry data, so as to alleviate the above-mentioned technical problems.

Example one

Fig. 1 is a flowchart of a method for monitoring quality of geostationary satellite external measurement data according to an embodiment of the present invention, and as shown in fig. 1, the method specifically includes the following steps:

and S102, receiving the current meteorological data and the real-time distance measurement data sent by the ground distance measurement master station.

Specifically, the method for monitoring quality of the geostationary satellite external measurement data provided by the embodiment of the present invention is mainly used for monitoring the overall quality of the real-time ranging data, and the real-time ranging data includes: the method provided by the embodiment of the invention is suitable for performing quality evaluation on any one of the ranging data. In order to reduce the influence of climate on the quality of the real-time ranging data, the meteorological data needs to be taken into consideration when the quality of the real-time ranging data is evaluated, so that the current meteorological data needs to be acquired besides the real-time ranging data sent by the ground ranging master station.

And step S104, acquiring satellite ephemeris files and ranging system errors output by the precision orbit determination and prediction of the daily orbit.

And S106, determining theoretical ranging data corresponding to the real-time ranging data based on the real-time ranging data and the satellite ephemeris file.

The error of the ranging system is a parameter output by the precision orbit determination of the daily orbit and is used for correcting real-time ranging data. In order to calculate theoretical ranging data corresponding to the real-time ranging data, when the method is operated, a satellite ephemeris file output by daily orbit prediction needs to be acquired, a geocentric position vector of a satellite at a certain time can be determined based on the satellite ephemeris file, and each piece of ranging data is provided with a ranging time mark in the received real-time ranging data sent by the ground ranging master station, so that the theoretical ranging data can be calculated based on the real-time ranging data and the satellite ephemeris file. The calculation method of the theoretical ranging data will be described in detail below.

And S108, determining the random error of the real-time ranging data based on the current meteorological data, the ranging system error, the real-time ranging data and the theoretical ranging data.

After theoretical ranging data are obtained through calculation, the atmospheric refraction correction amount can be calculated through the current meteorological data, the ranging system error can also correct the real-time ranging data by a certain amount, and the real-time ranging data are subtracted from the theoretical ranging data after the two corrections are performed on the real-time ranging data, so that the random error of the real-time ranging data can be determined; the size of the random error is used for representing the data quality of the real-time ranging data, obviously, the larger the random error is, the larger the difference between the real-time ranging data and the theoretical ranging data is, and the data quality of the real-time ranging data is poorer or the satellite orbit changes; on the contrary, if the random error is smaller, the data quality of the real-time ranging error is better. The embodiment of the invention does not specifically limit the random error, and a user can set a corresponding value range according to actual requirements, so that whether the quality of the real-time ranging data is qualified or not can be quickly judged.

According to the quality monitoring method for the synchronous satellite external measurement data, after the real-time distance measurement data are obtained, theoretical distance measurement data corresponding to the real-time distance measurement data are determined by using the real-time distance measurement data and a satellite ephemeris file output by daily orbit precise orbit determination and forecast, and then random errors of the real-time distance measurement data are determined by combining the current meteorological data and ranging system errors output by daily precise orbit determination. According to the method, the prior information such as the orbit dynamics model and the observation model is used for real-time processing of the satellite external data, the sensitivity of abnormal monitoring of the external data quality is improved, real-time monitoring and evaluation of the satellite external data can be realized based on random error change of real-time ranging data, and data support is provided for timely mastering the external data quality and adjusting the orbit determination strategy.

The method for monitoring the quality of the geostationary satellite external measurement data provided by the embodiment of the present invention is briefly described above, and the related method steps involved therein are specifically described below.

In an optional embodiment, as shown in fig. 2, in the step S106, determining theoretical ranging data corresponding to the real-time ranging data based on the real-time ranging data and the satellite ephemeris file specifically includes the following steps:

step S1061, obtaining a ranging time scale of the real-time ranging data.

Step S1062, determining a theoretical ranging reference value corresponding to the real-time ranging data based on the ranging time stamp and the satellite ephemeris file.

Specifically, as can be seen from the above description, the real-time ranging data all carry the ranging time scale, but the ranging time scale carried by the real-time ranging data is not the actual transmission time of the ranging pseudo code signal, in the embodiment of the present invention, the transmission time delay of the ranging pseudo code signal in the space is not considered when calculating the theoretical ranging reference value, that is, the transmission time of the ranging pseudo code signal is considered to be the same as the time of receiving the ranging pseudo code signal, and for the four-range ranging data, the receiving time of the ground ranging secondary station is also the same as the transmission time, so that by using the above features, the theoretical ranging reference value corresponding to each kind of real-time ranging data (one-way, two-way, four-range) can be calculated.

And step S1063, correcting the ranging time scale based on the theoretical ranging reference value and the preset fuzzy period to obtain the actual transmitting time of the ranging pseudo code signal.

Fig. 3 shows a schematic diagram of a pseudo code ranging measurement principle, as shown in fig. 3, measured ranging data R is obtained by calculating a time difference between a latest receiving time "1" and a latest sending time "1", and a ranging time mark is marked on the latest sending time mark "1", if the ranging data spans 6 ambiguities, when a receiving end receives the 1 st pulse, the sending end has reached the 7 th pulse, so the ranging time mark is marked on the 7 th pulse, the measured R value is a ranging value without ambiguity resolution, and the ranging time mark t is also an uncorrected time mark.

In the embodiment of the invention, a theoretical ranging reference value and a preset fuzzy period are used when the actual transmitting moment of the ranging pseudo code signal is calculated, and the preset fuzzy periods of the one-way ranging data, the two-way ranging data and the four-way ranging data are obtained by a device development party. When calculating the actual transmission time of the ranging pseudo code signals of the one-way, two-way and four-way ranging data, the corresponding theoretical ranging reference value and the preset fuzzy period are needed to be used. The specific calculation method will be described in detail below.

Step S1064, determining theoretical ranging data corresponding to the real-time ranging data based on the actual transmitting time and the satellite ephemeris file.

Specifically, after the actual transmission time of the ranging pseudo code signal is determined, the earth center position vector of the satellite at the actual transmission time of the ranging pseudo code signal can be determined by using the satellite ephemeris file, and for the one-way ranging data or the two-way ranging data, the theoretical ranging data corresponding to the one-way ranging data and the two-way ranging data can be determined only by combining the one-way ranging data or the two-way ranging data with the earth center position vector of the ground ranging master station; for the four-range ranging data, because the ranging process also involves the ground ranging secondary station, after the geocentric position vector of the satellite at the actual transmitting moment of the ranging pseudo code signal is determined, the theoretical ranging data corresponding to the four-range ranging data can be determined by combining the geocentric position vector of the ground ranging primary station and the geocentric position vector of the ground ranging secondary station.

In the above, a brief description is made on how to determine theoretical ranging data corresponding to real-time ranging data based on the real-time ranging data and the satellite ephemeris file, and a method for determining a theoretical ranging reference value corresponding to the real-time ranging data is described in detail below.

In an optional embodiment, in step S1062, the determining a theoretical ranging reference value corresponding to the real-time ranging data based on the ranging time stamp and the satellite ephemeris file specifically includes the following steps:

in step S10621, a geocentric position vector of the satellite at the time corresponding to the ranging time stamp is determined based on the ranging time stamp and the satellite ephemeris file.

Specifically, the satellite ephemeris file is a file for recording satellite operation orbit data, and the earth center position vector of the satellite at the moment corresponding to the ranging time mark can be determined by using the satellite ephemeris file and giving the ranging time mark.

Specifically, in the embodiment of the invention, the ephemeris file is used for carrying out interpolation calculation to solve the geocentric position vector of the satellite, the interpolation method adopts Lagrange polynomial interpolation, and the function is set at x₁,x₂,…,x_n(corresponds to the embodiment of the inventionTime corresponding to the earth center position vector of the satellite) are respectively y₁,y₂,…,y_n(corresponding to the geocentric position vector of the satellite in the embodiment of the present invention), the classical formula of lagrange polynomial interpolation is as follows:and determining the geocentric position vector of the satellite at the corresponding moment of the ranging time scale by using the interpolation formula. Generally, the satellite ephemeris interpolation selects 7-8 order polynomial for interpolation, and optionally, n is 8.

If the real-time ranging data is one-way ranging data or two-way ranging data, executing step S10622; if the real-time ranging data is the four-range ranging data, step S10623 is performed.

Step S10622, determining a theoretical ranging reference value corresponding to the one-way ranging data or a theoretical ranging reference value corresponding to the two-way ranging data based on the geocentric position vector of the satellite at the time corresponding to the ranging time scale and the geocentric position vector of the ground ranging master station.

Because the objects related to the one-way ranging data and the two-way ranging data are both a ground ranging master station and a satellite, on the premise of neglecting the transmission delay of the ranging pseudo code signal in the space, the theoretical ranging reference value corresponding to the two-way ranging data is considered to be 2 times of the theoretical reference value corresponding to the one-way ranging data, namely the theoretical ranging reference value is considered to be 2 times of the theoretical reference value corresponding to the one-way ranging dataWherein, ρ'₂Represents a theoretical ranging reference value, rho ', corresponding to the two-range ranging data'₁Representing theoretical reference values corresponding to single-pass ranging data,indicating the satellite at the corresponding time t on the ranging time scale₀The location vector of the earth's center of the earth,indicating the corresponding time t of the ground ranging master station at the ranging time scale₀The geocentric location vector.

Step S10623, determining theoretical ranging reference values corresponding to the four-range ranging data based on the geocentric position vector of the satellite at the corresponding time of the ranging time scale, the geocentric position vector of the ground ranging primary station and the geocentric position vector of the ground ranging secondary station.

Because the objects related to the four-range ranging data are a ground ranging main station, a satellite and a ground ranging secondary station, on the premise of neglecting the transmission delay of the ranging pseudo code signal in the space, the theoretical ranging reference value corresponding to the four-range ranging data can utilize a formulaTo perform a calculation in which, among others,indicating the corresponding time t of the ground ranging secondary station at the ranging time scale₀Is calculated (typically without change over time).

The process of calculating the theoretical ranging reference value corresponding to the real-time ranging data is described in detail above, and how to determine the theoretical ranging data corresponding to the real-time ranging data is described in detail below.

In an optional embodiment, in step S1064, the determining theoretical ranging data corresponding to the real-time ranging data based on the actual transmission time and the satellite ephemeris file specifically includes the following steps:

step S10641, determining the geocentric position vector of the satellite at the actual transmitting time based on the actual transmitting time and the satellite ephemeris file.

Specifically, in the same manner as in the step S10621, the geocentric position vector of the satellite at the actual transmission time of the ranging pseudo code signal can be obtained by using the lagrange polynomial interpolation formula.

If the real-time ranging data is one-way ranging data or two-way ranging data, executing step S10642-step S10643; if the real-time ranging data is the four-range ranging data, step S10644-step S10645 are performed.

Step S10642, determining the ranging time delay of the one-way ranging data or the ranging time delay of the two-way ranging data based on the geocentric position vector of the satellite at the actual transmitting time and the geocentric position vector of the ground ranging master station.

Step S10643, determining theoretical ranging data corresponding to the one-way ranging data based on the ranging time delay and the light speed of the one-way ranging data; or determining theoretical ranging data corresponding to the two-way ranging data based on the ranging time delay and the light speed of the two-way ranging data.

Specifically, when calculating theoretical ranging data corresponding to any one of the ranging data (one-way, two-way, or four-way), first, the ranging delay of the ranging data is calculated, and then the theoretical ranging data corresponding to the ranging data can be obtained by multiplying the ranging delay of the ranging data by the speed of light. For example, when calculating the theoretical ranging data corresponding to the one-way ranging data, the ranging delay of the one-way ranging data is calculated first, and then the theoretical ranging data corresponding to the one-way ranging data can be obtained by multiplying the ranging delay of the one-way ranging data by the speed of light.

The ranging time delay of the one-way ranging data is the time delay of the ranging pseudo code signal transmitted to the satellite by the ground ranging master station, and the ranging time delay of the one-way ranging data can use a formulaIs shown, wherein, τ₁Representing the ranging delay, t, of one-way ranging data₂Indicating the moment at which the satellite receives the ranging pseudocode signal, t₁Indicating the time at which the ground ranging master station transmits the ranging pseudo code signal,represents t₂The vector of the geocentric position of the satellite at the moment,represents t₁And (5) the earth center position vector of the earth distance measuring master station at the moment. At the moment of transmission t when the ranging pseudo code signal has been determined₁And the earth center position vector of the satellite at the actual transmission time of the ranging pseudo code signalThe precise coordinates of the ground station can be known, and the geocentric position vector of the corresponding moment can be obtained through coordinate conversion, so that the moment t when the satellite receives the ranging pseudo code signal₂And t₂Earth center position vector of time satelliteIs an unknown number.

The embodiment of the invention uses the following method to calculate the ranging time delay of any section of forwarding interval of which the receiving end position of the ranging signal changes along with the time, and assumes that a certain section of forwarding interval of the ranging pseudo code signal is sent from a target A to a target B, and the time when the ranging signal starts from the target A is t_sThen the transmission delay of the ranging signal forwarded from target a to target BWherein the target B is at t_sGeocentric location vector at + τ timeFor unknowns, target A is at t_sGeocentric position vector of timeFor a known number, c denotes the speed of light, τ can be solved using the following algorithm:

1) computingAnd determines whether or not Δ τ ═ τ ″ - τ' is less than 10^-12(ii) a Wherein, when calculating for the first time, let τ' be 0;

2) if not, taking the currently calculated tau ' as the tau ' in the next iterative calculation, and iteratively calculating the tau ' of the next generation by using the formula in the step 1) until the convergence of the delta tau (delta tau) is achieved<10^-12) And taking the tau' meeting the convergence condition as the ranging signal transmission time delay tau of the target A to the target B.

By the above methodDistance measuring time delay tau capable of obtaining one-way distance measuring data₁. Then, will tau₁And multiplying the light speed to obtain theoretical ranging data corresponding to the one-way ranging data. According toThe time t of the satellite receiving the ranging pseudo code signal can be calculated₂And further combining the satellite ephemeris file to calculate t₂Earth center position vector of time satellite

Fig. 4 is a schematic diagram of two-way ranging and four-way ranging according to an embodiment of the present invention, in fig. 4, a solid line represents a flow direction of a ranging pseudo code signal during two-way ranging, and a dotted line represents a flow direction of the ranging pseudo code signal during four-way ranging. The ranging delay of the two-way ranging data should be denoted as tau₁+τ₂Whereinτ₁Indicating the time delay, tau, of the ranging pseudocode signal transmitted by the ground ranging master station to the satellite₂Representing the time delay of the satellite returning the ranging pseudocode signal to the ground ranging master station, tau having been calculated above₁，t₂Andthen, using the above-described method for determining the ranging delay of the forwarding interval in which the receiving end position of any ranging signal changes with time, τ can be determined₂And time t when the ranging pseudo code signal returns to the ground ranging master station₃And then determining the ranging time delay tau of the two-way ranging data₁+τ₂. Then, will tau₁+τ₂And multiplying the light velocity to obtain theoretical ranging data corresponding to the two-way ranging data.

Step S10644, determining the ranging time delay of the four-range ranging data based on the geocentric position vector of the satellite at the actual transmitting time, the geocentric position vector of the primary ground ranging station and the geocentric position vector of the secondary ground ranging station.

Step S10645, determining theoretical ranging data corresponding to the four-range ranging data based on the ranging delay and the light speed of the four-range ranging data.

The ranging delay of the four-range ranging data should be denoted as tau₁+τ₃+τ₄+τ₅Whereinτ₁Indicating the time delay, tau, of the ranging pseudocode signal transmitted by the ground ranging master station to the satellite₃Representing the time delay, τ, of the satellite relaying the ranging pseudocode signal to the ground ranging secondary station₄Representing the time delay, τ, of the ground ranging secondary station relaying the ranging pseudocode signal to the satellite₅Representing the time delay for the satellite to forward the ranging pseudocode signal to the ground ranging master station.

τ has been calculated above₁，t₂Andτ₃in the formula, the unknown parameter is the time t when the ranging pseudo code signal is forwarded to the ground ranging secondary station₄And t and₄earth center position vector of time ground ranging secondary stationTau can be obtained by using the method for obtaining the ranging time delay of the forwarding interval of the receiving end position of any section of ranging signal changing with time described above₃And t₄。

Further, τ₄In the formula (III), the unknown parameter is the time t when the ranging pseudo code signal is transmitted to the satellite₅And t and₅earth center position vector of time satelliteTau can be obtained by using the method for obtaining the ranging time delay of the forwarding interval of the receiving end position of any section of ranging signal changing with time described above₄，t₅Further combining with satellite ephemeris fileCalculate out

Further, τ₅In the formula, the unknown parameter is the time t when the ranging pseudo code signal is transmitted to the ground ranging master station₆And t and₆earth center position vector of time ground ranging master stationTau can be obtained by using the method for obtaining the ranging time delay of the forwarding interval of the receiving end position of any section of ranging signal changing with time described above₅And t₆And further determining the ranging time delay tau of the four-range ranging data₁+τ₃+τ₄+τ₅Then, t is added₁+τ₃+τ₄+τ₅And multiplying the four-range distance measurement data by the light speed to obtain theoretical distance measurement data corresponding to the four-range distance measurement data.

In the above, a detailed description is given to a method for determining theoretical ranging data corresponding to any real-time ranging data, and a detailed description is given to a method for calculating a random error of the real-time ranging data.

In an optional embodiment, in the step S108, determining the random error of the real-time ranging data based on the current meteorological data, the ranging system error, the real-time ranging data, and the theoretical ranging data specifically includes the following steps:

and step S1081, outlier identification and elimination are carried out on the current meteorological data and the real-time ranging data, and target meteorological data and target ranging data are obtained.

The satellite external data is observed by a ground distance measuring master station in real time, the data is processed quickly, outliers contained in the data are identified and rejected, and an external data outlier rejection algorithm which not only saves the memory, but also ensures the requirements on precision and real-time performance must be established. In the embodiment of the invention, a sliding data fitting method is utilized, fitting coefficients are continuously improved, fitting precision is improved, and the reasonability of observed data is checked by using a function through prediction.

In particularLet y_j(j is 1 to n) represents the observed amount of the external measurement data at a certain time, and y is considered to be general in order to avoid loss of generality_jIs a second order function of time t:since the measurements are equally spaced, i.e. t_jJ is a sampling interval (constant value), and j is a sampling point sequence. Suitable fitting coefficients may be determined according to the least squares principle so thatThe minimum variance E is called the optimal binary fit approximation.

Observe data y for the dot sequence j_jRespectively taking the first n data: y is_j-n,…y_j-1Equation of arithmeticFitting coefficient a of polynomial on both sides₀,a₁,a₂The derivative is taken and made zero and the fitting value of the available coefficients is: wherein α is 2ⁿ⁺¹，

From this, a fitting function of the observed data can be determined

Thus, the fitting value can be usedChecking the original observed value y_jWhether or not to be distorted, if so, using the fitting valueThe specific idea of repairing data by fitting method is described above by taking the jth point data as an example.

Because the data observation is continuous, namely because the target is in continuous motion, the fitting function needs to be updated according to the motion change, and the efficiency of outlier identification and elimination can be improved. Therefore, using the newly observed data, the following sliding method is aided, and the coefficient fitting values are continuously improvedThe fitting precision is improved, and the reasonability of the observed data is checked by forecasting through the function.

As can be seen from the above description, the procedure of identifying and rejecting outliers of the external data is as follows:

(1) constructing n external data buffer areas;

(2) n data y accumulated by buffer₁,…y_nConstructing a fitting function

(3) For the new observation data y obtained by the step (2)_n+1Predicted value of

(4) Judgment ofIs true, where σ₂Representing a field value judgment threshold value, and generally taking 17 times of equipment observation precision index values;

(5) if the inequality is true, then y is indicated_n+1Is a normal data point, and the buffer removes the first point y using the first-in-first-out principle₁Supplement y_n+1Go to the buffer to repeat the above process;

(6) if the inequality is not true, useIn place of y_n+1And (4) repeating the steps (1) to (4).

In the judgment process, if a plurality of (for example, 5) new observation data can not make the inequality in the step (4) hold, discarding the data in the buffer, and reselecting and constructing a new fitting function.

Therefore, before calculating the random error, the above method is firstly used to identify and reject the outliers of the current meteorological data and the real-time ranging data to obtain the target meteorological data and the target ranging data.

And step S1082, determining an atmospheric refraction correction amount corresponding to the ranging data based on the target meteorological data.

The meteorological data is a basic basis for correcting the refraction error of the observation data, so that the ranging data can be corrected in real time according to the meteorological data sent by the ground ranging master station.

And step S1083, correcting the target ranging data based on the theoretical ranging reference value and the preset fuzzy period to obtain corrected ranging data.

In the above, the process of correcting the ranging time scale is described, and similar to the principle of time scale correction, the target ranging data also needs to be corrected, and specifically, the embodiment of the present invention uses the formulaTo calculate corrected ranging data, where ρ ″ "_cRepresenting corrected range data, p_cRepresents real-time ranging data ρ'_cAnd c represents the light speed, and delta T represents the preset fuzzy period corresponding to the real-time ranging data. When the single-pass, double-pass and four-pass ranging data are corrected, the corresponding theoretical ranging reference values and the corresponding preset fuzzy periods need to be used.

And step S1084, determining a random error of the real-time ranging data based on the atmospheric refraction correction amount, the ranging system error, the ranging data after correction and the theoretical ranging data.

After calculating the atmospheric refraction correction amount, the ranging system error, the ranging data after correction and the theoretical ranging data, the embodiment of the invention adopts the formulaTo calculate the random error of the real-time ranging data, wherein Δ σ represents the random error of the real-time ranging data, ρ ″_cRepresents the range-finding data after the correction,representing theoretical distance measurement data, Δ ρ, corresponding to real-time distance measurement data_tropIndicating the amount of atmospheric refraction correction, Δ ρ_sIndicating a ranging system error.

For data in a tracking arc segment, random error mean square error can be calculated to further evaluate the quality of the external data, and the estimated value of the random error mean square error isWhere M represents the total number of measurements within the tracking arc segment.

In an optional implementation manner, in step S1063, the step of correcting the ranging time scale based on the theoretical ranging reference value and the preset fuzzy period to obtain the actual transmission time of the ranging pseudo code signal specifically includes the following steps:

correcting the ranging time scale based on the formula T ═ T-N · delta T to obtain the actual transmitting time of the ranging pseudo code signal; where t 'denotes the actual transmission time of the ranging pseudo code signal, t denotes the time corresponding to the ranging time scale, N denotes the ambiguity number, and N ═ ρ'_c/Δρ|，ρ′_cDenotes a theoretical range reference, Δ ρ ═ Δ T × c, Δ T denotes a preset period of blurring, and c denotes the speed of light.

Specifically, the embodiment of the present invention utilizes equationsAnd calculating the actual transmitting time of the ranging pseudo code signal, wherein the delta rho is delta T, and can be understood as fuzzy distance, and the preset fuzzy periods corresponding to the one-way ranging data, the two-way ranging data and the four-way ranging data are related to equipment. When calculating the actual transmission time of the ranging pseudo code signals of the one-way, two-way and four-way ranging data, the corresponding theoretical ranging reference value and the corresponding preset fuzzy period are used.

In an optional embodiment, in step S1082, the atmospheric refraction correction amount corresponding to the ranging data is determined based on the target meteorological data, which specifically includes the following contents:

equation of utilizationDetermining the atmospheric refraction correction corresponding to the distance measurement data, wherein DeltaR represents the atmospheric refraction correction, C represents the atmospheric refractive index, E_NWhich represents the elevation angle of a theoretical survey station,which represents the refractive index of the ground,t represents absolute temperature, P represents ground atmospheric pressure, P_eRepresenting the ground water pressure.

For simplifying meteorological observation, the atmospheric refractive index C can be taken as a statistical average value 1.4142E-4 (1/m), and the absolute temperature T is T ═ T₀+ T', where T₀273.15 deg.C, T' represents the ground temperature and the ground water pressureU represents ground relative humidity.

In summary, the embodiment of the present invention provides a high-precision method for monitoring quality of geostationary satellite external data, in which prior information such as an orbit dynamics model (using satellite ephemeris) and an observation model (combining with ranging system errors) is used for real-time processing of the geostationary satellite external data, so as to improve sensitivity of abnormal monitoring of the external data quality, realize real-time monitoring and evaluation of the satellite external data based on random error change of the real-time ranging data, and provide data support for timely mastering the external data quality and adjusting an orbit determination strategy.

Fig. 5 shows a quality monitoring interface of multi-station two-way and four-way ranging data, a left view in fig. 5 shows original observation data which is sent by a ground ranging master station and is subjected to fuzzy distance and time scale correction and outlier identification and rejection, and displayed numbers are measured values at the current observation time. The right side view in fig. 5 is the random error, and the number shown is the mean of the random errors for the arc statistics. As can be seen from fig. 5, the secondary station 2 is out of lock, and the measured data is abnormal.

Example two

The embodiment of the invention also provides a quality monitoring device for the geostationary satellite external data, which is mainly used for executing the quality monitoring method for the satellite external data provided by the embodiment of the invention, and the quality monitoring device for the geostationary satellite external data provided by the embodiment of the invention is specifically introduced below.

Fig. 6 is a functional block diagram of an apparatus for monitoring quality of geostationary satellite external data according to an embodiment of the present invention, as shown in fig. 6, the apparatus mainly includes: the device comprises a receiving module 10, an obtaining module 20, a first determining module 30 and a second determining module 40, wherein:

and the receiving module 10 is used for receiving the current meteorological data and the real-time ranging data sent by the ground ranging master station.

And the obtaining module 20 is configured to obtain a satellite ephemeris file and a ranging system error output by the precision orbit determination and prediction of the daily orbit.

The first determining module 30 is configured to determine theoretical ranging data corresponding to the real-time ranging data based on the real-time ranging data and the satellite ephemeris file.

A second determining module 40, configured to determine a random error of the real-time ranging data based on the current meteorological data, the ranging system error, the real-time ranging data, and the theoretical ranging data; the random error is used for representing the data quality of the real-time ranging data.

The invention provides a quality monitoring device for synchronous satellite external measurement data, which comprises: and the receiving module 10 is used for receiving the current meteorological data and the real-time ranging data sent by the ground ranging master station. And the obtaining module 20 is configured to obtain a satellite ephemeris file and a ranging system error output by the precision orbit determination and prediction of the daily orbit. The first determining module 30 is configured to determine theoretical ranging data corresponding to the real-time ranging data based on the real-time ranging data and the satellite ephemeris file. A second determining module 40, configured to determine a random error of the real-time ranging data based on the current meteorological data, the ranging system error, the real-time ranging data, and the theoretical ranging data; the random error is used for representing the data quality of the real-time ranging data. The device uses the prior information such as the orbit dynamics model, the observation model and the like for the real-time processing of the synchronous satellite external data, improves the sensitivity of the abnormal monitoring of the external data quality, can realize the real-time monitoring and evaluation of the synchronous satellite external data based on the random error change of the real-time ranging data, and provides data support for timely mastering the external data quality and adjusting the orbit determination strategy.

Optionally, the first determining module 30 includes:

and the acquisition unit is used for acquiring the ranging time scale of the real-time ranging data.

And the first determining unit is used for determining a theoretical ranging reference value corresponding to the real-time ranging data based on the ranging time mark and the satellite ephemeris file.

And the first correcting unit is used for correcting the ranging time scale based on the theoretical ranging reference value and the preset fuzzy period to obtain the actual transmitting time of the ranging pseudo code signal.

And the second determining unit is used for determining theoretical ranging data corresponding to the real-time ranging data based on the actual transmitting time and the satellite ephemeris file.

Optionally, the first determining unit is specifically configured to:

and determining the earth center position vector of the satellite at the corresponding moment of the ranging time mark based on the ranging time mark and the satellite ephemeris file.

And if the real-time ranging data is single-pass ranging data or double-pass ranging data, determining a theoretical ranging reference value corresponding to the single-pass ranging data or a theoretical ranging reference value corresponding to the double-pass ranging data based on the geocentric position vector of the satellite at the corresponding moment of the ranging time scale and the geocentric position vector of the ground ranging master station.

And if the real-time ranging data are four-range ranging data, determining theoretical ranging reference values corresponding to the four-range ranging data based on the geocentric position vector of the satellite at the corresponding moment of the ranging time scale, the geocentric position vector of the ground ranging main station and the geocentric position vector of the ground ranging auxiliary station.

Optionally, the second determining unit is specifically configured to:

and determining the geocentric position vector of the satellite at the actual transmission moment based on the actual transmission moment and the satellite ephemeris file.

And if the real-time ranging data is single-pass ranging data or double-pass ranging data, determining the ranging time delay of the single-pass ranging data or the ranging time delay of the double-pass ranging data based on the geocentric position vector of the satellite at the actual transmitting moment and the geocentric position vector of the ground ranging master station.

Determining theoretical ranging data corresponding to the one-way ranging data based on the ranging time delay and the light speed of the one-way ranging data; or determining theoretical ranging data corresponding to the two-way ranging data based on the ranging time delay and the light speed of the two-way ranging data.

And if the real-time ranging data is the four-range ranging data, determining the ranging time delay of the four-range ranging data based on the geocentric position vector of the satellite at the actual transmitting moment, the geocentric position vector of the ground ranging main station and the geocentric position vector of the ground ranging auxiliary station.

And determining theoretical ranging data corresponding to the four-range ranging data based on the ranging time delay and the light speed of the four-range ranging data.

Optionally, the second determining module 40 includes:

and the outlier identification and rejection unit is used for identifying and rejecting outliers of the current meteorological data and the real-time ranging data to obtain target meteorological data and target ranging data.

And the third determining unit is used for determining the atmospheric refraction correction amount corresponding to the ranging data based on the target meteorological data.

And the second correction unit is used for correcting the target ranging data based on the theoretical ranging reference value and the preset fuzzy period to obtain the corrected ranging data.

And the third determining unit is used for determining the random error of the real-time ranging data based on the atmospheric refraction correction amount, the ranging system error, the ranging data after correction and the theoretical ranging data.

Optionally, the first correction unit is specifically configured to:

correcting the ranging time scale based on the formula T ═ T-N · delta T to obtain the actual transmitting time of the ranging pseudo code signal; where t 'denotes the actual transmission time of the ranging pseudo code signal, t denotes the time corresponding to the ranging time scale, N denotes the ambiguity number, and N ═ ρ'_c/Δρ|，ρ′_cDenotes a theoretical range reference, Δ ρ ═ Δ T × c, Δ T denotes a preset period of blurring, and c denotes the speed of light.

Optionally, the third determining unit is specifically configured to:

equation of utilizationDetermining the atmospheric refraction correction corresponding to the distance measurement data, wherein DeltaR represents the atmospheric refraction correction, C represents the atmospheric refractive index, E_NWhich represents the elevation angle of a theoretical survey station,which represents the refractive index of the ground,t represents absolute temperature, P represents ground atmospheric pressure, P_eRepresenting the ground water pressure.

EXAMPLE III

Referring to fig. 7, an embodiment of the present invention provides an electronic device, including: a processor 60, a memory 61, a bus 62 and a communication interface 63, wherein the processor 60, the communication interface 63 and the memory 61 are connected through the bus 62; the processor 60 is arranged to execute executable modules, such as computer programs, stored in the memory 61.

The memory 61 may include a high-speed Random Access Memory (RAM) and may also include a non-volatile memory (non-volatile memory), such as at least one disk memory. The communication connection between the network element of the system and at least one other network element is realized through at least one communication interface 63 (which may be wired or wireless), and the internet, a wide area network, a local network, a metropolitan area network, and the like can be used.

The bus 62 may be an ISA bus, PCI bus, EISA bus, or the like. The bus may be divided into an address bus, a data bus, a control bus, etc. For ease of illustration, only one double-headed arrow is shown in FIG. 7, but this does not indicate only one bus or one type of bus.

The memory 61 is used for storing a program, the processor 60 executes the program after receiving an execution instruction, and the method executed by the apparatus defined by the flow process disclosed in any of the foregoing embodiments of the present invention may be applied to the processor 60, or implemented by the processor 60.

The processor 60 may be an integrated circuit chip having signal processing capabilities. In implementation, the steps of the above method may be performed by integrated logic circuits of hardware or instructions in the form of software in the processor 60. The Processor 60 may be a general-purpose Processor, and includes a Central Processing Unit (CPU), a Network Processor (NP), and the like; the device can also be a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field-Programmable Gate Array (FPGA), or other Programmable logic devices, discrete Gate or transistor logic devices, discrete hardware components. The various methods, steps and logic blocks disclosed in the embodiments of the present invention may be implemented or performed. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like. The steps of the method disclosed in connection with the embodiments of the present invention may be directly implemented by a hardware decoding processor, or implemented by a combination of hardware and software modules in the decoding processor. The software module may be located in ram, flash memory, rom, prom, or eprom, registers, etc. storage media as is well known in the art. The storage medium is located in a memory 61, and the processor 60 reads the information in the memory 61 and, in combination with its hardware, performs the steps of the above method.

The method, the apparatus, and the computer program product for monitoring quality of geostationary satellite external data provided in the embodiments of the present invention include a computer-readable storage medium storing a non-volatile program code executable by a processor, where instructions included in the program code may be used to execute the method described in the foregoing method embodiments, and specific implementation may refer to the method embodiments, and will not be described herein again.

In addition, functional units in the embodiments of the present invention may be integrated into one processing unit, or each unit may exist alone physically, or two or more units are integrated into one unit.

The functions, if implemented in the form of software functional units and sold or used as a stand-alone product, may be stored in a non-volatile computer-readable storage medium executable by a processor. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.

It should be noted that: like reference numbers and letters refer to like items in the following figures, and thus, once an item is defined in one figure, it need not be further defined and explained in subsequent figures.

In the description of the present invention, it should be noted that the terms "center", "upper", "lower", "left", "right", "vertical", "horizontal", "inner", "outer", etc. indicate orientations or positional relationships based on the orientations or positional relationships shown in the drawings or the orientations or positional relationships that the products of the present invention are conventionally placed in use, and are only used for convenience in describing the present invention and simplifying the description, but do not indicate or imply that the devices or elements referred to must have a specific orientation, be constructed and operated in a specific orientation, and thus, should not be construed as limiting the present invention. Furthermore, the terms "first," "second," "third," and the like are used solely to distinguish one from another and are not to be construed as indicating or implying relative importance.

Furthermore, the terms "horizontal", "vertical", "overhang" and the like do not imply that the components are required to be absolutely horizontal or overhang, but may be slightly inclined. For example, "horizontal" merely means that the direction is more horizontal than "vertical" and does not mean that the structure must be perfectly horizontal, but may be slightly inclined.

In the description of the present invention, it should also be noted that, unless otherwise explicitly specified or limited, the terms "disposed," "mounted," "connected," and "connected" are to be construed broadly and may, for example, be fixedly connected, detachably connected, or integrally connected; can be mechanically or electrically connected; they may be connected directly or indirectly through intervening media, or they may be interconnected between two elements. The specific meanings of the above terms in the present invention can be understood in specific cases to those skilled in the art.

Finally, it should be noted that: the above embodiments are only used to illustrate the technical solution of the present invention, and not to limit the same; while the invention has been described in detail and with reference to the foregoing embodiments, it will be understood by those skilled in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; and the modifications or the substitutions do not make the essence of the corresponding technical solutions depart from the scope of the technical solutions of the embodiments of the present invention.