CN114852375A - Method and device for estimating relative orbit change of formation satellite - Google Patents
Method and device for estimating relative orbit change of formation satellite Download PDFInfo
- Publication number
- CN114852375A CN114852375A CN202210302728.4A CN202210302728A CN114852375A CN 114852375 A CN114852375 A CN 114852375A CN 202210302728 A CN202210302728 A CN 202210302728A CN 114852375 A CN114852375 A CN 114852375A
- Authority
- CN
- China
- Prior art keywords
- relative
- orbit
- target
- eccentricity
- satellite
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 19
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 15
- 238000004364 calculation method Methods 0.000 claims abstract description 14
- 238000001914 filtration Methods 0.000 claims abstract description 13
- 238000002347 injection Methods 0.000 claims description 11
- 239000007924 injection Substances 0.000 claims description 11
- 230000005484 gravity Effects 0.000 claims description 6
- 238000004590 computer program Methods 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 12
- 238000013213 extrapolation Methods 0.000 abstract description 3
- 238000004088 simulation Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 239000000243 solution Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/244—Spacecraft control systems
- B64G1/245—Attitude control algorithms for spacecraft attitude control
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Radar, Positioning & Navigation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Automation & Control Theory (AREA)
- Navigation (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
A method and a device for estimating relative orbit change of formation satellites do not depend on high-precision relative orbit root estimation of inter-satellite measurement links, real-time extrapolation calculation is carried out by injecting target satellite orbit root parameters on the ground, a recursive filtering algorithm based on relative orbit kinematics is proposed to carry out on-line estimation on a long-period drifting item, the problem of high-precision estimation of relative orbit change is solved under the condition that inter-satellite measurement is infeasible, and the high-precision relative orbit root estimation which does not depend on the inter-satellite measurement links is realized.
Description
Technical Field
The invention relates to a method and a device for estimating relative orbit change of formation satellites, and belongs to the technical field of spacecraft control.
Background
And generating a new model of the earth gravitational field by the low-low tracking gravity measurement satellite every 15-30 days in the task period. The mission requires that two satellites are 170-270 km apart and are formed into a formation to fly in tandem. The variation of the earth gravitational field is represented as the difference of the perturbations of the two common rail satellite orbits. The traditional method needs to measure the change of the distance between two satellites by relying on a high-precision high-frequency inter-satellite measurement link, the accurate measurement value of a triaxial accelerometer is used for eliminating the contribution of a non-gravity component to the change of the distance between the satellites after being processed, and a high-precision GNSS is used for accurately determining the orbit of the satellite, so that the spatial distribution of gravity parameters is measured, and an earth gravitational field model is estimated. In some satellite application scenes, a real-time inter-satellite communication means does not exist between two satellites, so that when double-satellite inter-satellite pointing control is carried out, double-satellite relative distance information cannot be directly obtained through measurement, and new requirements for relative distance estimation are provided.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: the method and the device for estimating the relative orbit change of the formation satellites overcome the defects of the prior art, do not depend on the high-precision relative orbit root estimation of the inter-satellite measurement link, perform real-time extrapolation calculation by injecting target satellite orbit root parameters on the ground, and provide a recursive filtering algorithm based on relative orbit kinematics to perform online estimation on a long-period drift term, solve the problem of high-precision estimation of the relative orbit change under the condition that inter-satellite measurement is infeasible, and realize the high-precision relative orbit root estimation which does not depend on the inter-satellite measurement link.
The purpose of the invention is realized by the following technical scheme:
the embodiment of the invention provides a method for estimating relative orbit change of formation satellites, which comprises the following steps:
obtaining effective local satellite orbit parameters and target satellite orbit parameters injected from the ground in the period;
estimating relative flat track angular velocity, relative track eccentricity component and relative track eccentricity change rate by adopting a filtering algorithm;
calculating the relative star spacing, the relative semi-major axis and the relative eccentricity variation;
the relative semi-major axis target amount and the drift velocity target amount are calculated.
Preferably, it is first determined whether the local satellite orbit parameter and the target satellite orbit parameter injected from the ground are received in the present period, and if the local satellite orbit parameter and the target satellite orbit parameter are received, it is further determined whether the local satellite orbit parameter and the target satellite orbit parameter are valid.
Preferably, the effective local star orbit parameters and target star orbit parameters of the surface injection at least include: the reference time of the relative orbit and the variation thereof, the relative star spacing of the two stars and the variation thereof, the relative orbit eccentricity component, the relative orbit semi-major axis of the orbit injection, and the average value of the reference time of the two stars orbit.
Preferably, the filtering algorithm is used to estimate the relative flat track angular velocity, the relative track eccentricity component, and the relative track eccentricity change rate, and the iteration process is as follows:
e xRM =e*cos(ω)-e T *cos(ω T )
e yRM =e*sin(ω)-e T *sin(ω T )
e xdotRM =k fedot (e xRM -e xRM0 )/Dt RM +(1-k fedot )e xdotRM0
e ydotRM =k fedot (e yRM -e yRM0 )/Dt RM +(1-k fedot )e ydotRM0
in the formula, k fedot For filter gain, e is eccentricity, ω is orbital argument of perigee, v RM As drift velocity, v dotRM As rate of change of drift velocity, e xRM And e yRM As a component of relative track eccentricity, e xdotRM And e ydotRM For the amount of change thereof,is a semi-major axis flat root,the orbital angular velocity, Re, the earth radius, μ 398600km 3 /s 2 Subscript 0 represents the state variable of the last beat, and subscript T represents the state variable of the target star, which is the gravity constant; t is t RM0 Reference time, Dt, relative to the track RM0 Is t RM0 Amount of change of (Dr) RM0 Is the relative star spacing of two stars, Dr RM Is Dr RM0 Amount of change of e xRM 、e yRM Are all relative orbital eccentricity components, Δ a kRM Relative orbit semi-major axis, t, for orbit injection midRM Is the average of the two-star orbit reference time.
Preferably, the calculating the relative inter-satellite distance r RM Relative semi-major axis Deltaa RM Relative eccentricity change amount Δ e xRM And Δ e yRM The method comprises the following steps:
Δt RM =t-t RM0
Δe xRM =e xRM +e xdotRM Δt RM
Δe yRM =e yRM +e ydotRM Δt RM 。
preferably, the calculation is relative to the semimajor axis target amount Δ a TRM And a target amount v of drift velocity TRM The method comprises the following steps:
wherein, Δ L TRM Is the target star spacing variation.
The embodiment of the invention provides a device for estimating relative orbit change of formation satellites, which comprises:
the acquisition module is used for receiving effective local satellite orbit parameters and target satellite orbit parameters injected into the ground in the period;
the estimation module estimates the relative flat track angular velocity, the relative track eccentricity component and the relative track eccentricity change rate by adopting a filtering algorithm;
the first calculation module is used for calculating the relative star spacing, the relative semi-major axis and the relative eccentricity variation;
and the second calculation module is used for calculating the relative semi-major axis target quantity and the drift velocity target quantity.
Preferably, the obtaining module first determines whether the local satellite orbit parameter and the target satellite orbit parameter injected from the ground are received in the current period, and if the local satellite orbit parameter and the target satellite orbit parameter are received, further determines whether the local satellite orbit parameter and the target satellite orbit parameter are valid.
Preferably, the effective local star orbit parameters and target star orbit parameters of the surface injection at least include: the reference time of the relative orbit and the variation thereof, the relative star spacing of the two stars and the variation thereof, the relative orbit eccentricity component, the relative orbit semi-major axis of the orbit injection, and the average value of the reference time of the two stars orbit.
Embodiments of the present invention provide a computer-readable storage medium having stored thereon computer program instructions, which, when loaded and executed by a processor, cause the processor to perform the above-mentioned method for estimating relative orbital change of a formation satellite.
Compared with the prior art, the invention has the following beneficial effects:
(1) the invention provides a method for estimating relative orbit change of formation satellites, which is a novel solving method aiming at the relative orbit position change of the formation satellites and breaks through the limitation that the satellites such as high-precision gravity measurement need to rely on inter-satellite links to measure the distance change between two satellites.
(2) The method fully excavates the relative orbit kinematics characteristics of the low-low tracking satellite, skillfully solves the problem of estimating the relative distance of two satellites in formation by combining nineteen relative orbit extrapolation with a novel filtering design, realizes high-precision estimation of the inter-satellite distance without depending on inter-satellite distance measurement, and is simple and feasible.
(3) The invention provides a new solution for the problem of relative distance estimation in the low-low tracking satellite formation process, no additional data input is needed, the calculation is simple, the requirement of double-satellite relative position prediction accuracy can be well met by using the algorithm, and the algorithm has strong engineering practicability.
Drawings
FIG. 1 is a control flow chart of the present invention.
Fig. 2 is a prediction error curve of relative positions of two stars.
FIG. 3 is a plot of star spacing versus relative semi-major axis.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, embodiments of the present invention will be described in detail with reference to the accompanying drawings.
As shown in fig. 1, the method for estimating relative orbital change of a formation satellite of the present invention includes the following steps:
firstly, initializing orbit parameters on the premise that the orbit parameters of a local satellite and the orbit parameters of a target satellite which are injected from the ground are received and confirmed to be effective in the period;
calculating relative flat track angular velocity, relative track eccentricity vector and relative track eccentricity change rate based on filtering estimation;
calculating relative star spacing, relative semi-major axis and relative eccentricity variation;
and step four, calculating the relative semi-major axis target quantity and the drift velocity target quantity.
The specific process of the step one is as follows:
judging whether the local satellite orbit parameter and the target satellite orbit parameter injected from the ground are received in the period, if so, further confirming whether the local satellite orbit parameter and the target satellite orbit parameter are effective, and performing initialization assignment on the orbit parameters under the effective condition, wherein the relevant parameters comprise: reference time t relative to the track RM0 And its variation Dt RM0 Double star relative star spacing Dr RM0 And its variation Dr RM Relative track eccentricity component e xRM 、e yRM Relative orbit semi-major axis Deltaa of orbit injection kRM Mean value t of reference time of two-star orbit midRM 。
In the second step, the angular velocity of the relative flat track, the component of the relative track eccentricity and the change rate of the relative track eccentricity are estimated by using a filtering algorithm, namely
e xRM =e*cos(ω)-e T *cos(ω T )
e yRM =e*sin(ω)-e T *sin(ω T )
e xdotRM =k fedot (e xRM -e xRM0 )/Dt RM +(1-k fedot )e xdotRM
e ydotRM =k fedot (e yRM -e yRM0 )/Dt RM +(1-k fedot )e ydotRM
In the formula, k fedot For filter gain, e is eccentricity, ω is orbital argument of perigee, v RM As drift velocity, v dotRM As rate of change of drift velocity, e xRM And e yRM As a component of relative track eccentricity, e xdotRM And e ydotRM For the variation, a is the half-major axis flat root, n is the orbital angular velocity, Re is the earth radius, μ 398600km 3 /s 2 Is the gravitational constant, () 0 State variable of the last beat, () T Is the state variable of the target star.
In the third step, on the basis of the calculation of the second variable in the step, the relative satellite distance r is calculated RM Relative semi-major axis Deltaa RM And relative eccentricity change amount Δ e xRM 、Δe yRM The calculation is carried out according to the following specific algorithm
Note Δ t RM =t-t RM0
Δe xRM =e xRM +e xdotRM Δt RM
Δe yRM =e yRM +e ydotRM Δt RM
In the fourth step, the target amount Δ a is calculated with respect to the semimajor axis TRM And a target amount v of drift velocity TRM Is calculated, i.e.
Wherein, Δ L TRM A modification may be injected for the target star spacing variation.
The invention is further illustrated by the following examples.
Example 1:
taking a certain satellite as an example, mathematical simulation analysis is carried out on the relative orbit change estimation algorithm provided by the invention, the satellite runs in a sun synchronous circular orbit of 500km, and two satellite orbit parameters are set at the beginning and shown in table 1.
TABLE 1
In the simulation, star A is in front, the star spacing is 269km, the semi-major axis of star A is slightly lower, and two stars continue to be far away. Performing 200-day simulation, wherein simulation results are shown in fig. 2 to fig. 3, fig. 2 is a relative semi-major axis curve, mutation at X165.6 is caused by the star pitch maintaining orbit control, and fig. 3 is a relationship curve between the star pitch and the relative semi-major axis; the change in relative semimajor axis difference is not absolutely uniform, which is related to the change in atmospheric density.
Example 2:
a method for estimating relative orbit change of formation satellites comprises the following steps:
obtaining effective local satellite orbit parameters and target satellite orbit parameters injected into the ground in the period;
estimating relative flat track angular velocity, relative track eccentricity component and relative track eccentricity change rate by adopting a filtering algorithm;
calculating the relative star spacing, the relative semi-major axis and the relative eccentricity variation;
the relative semi-major axis target amount and the drift velocity target amount are calculated.
An apparatus for estimating relative orbital variation of a formation satellite, comprising:
the acquisition module is used for receiving effective local satellite orbit parameters and target satellite orbit parameters injected into the ground in the period;
the estimation module estimates the relative flat track angular velocity, the relative track eccentricity component and the relative track eccentricity change rate by adopting a filtering algorithm;
the first calculation module is used for calculating the relative star spacing, the relative semi-major axis and the relative eccentricity variation;
and the second calculation module is used for calculating the relative semi-major axis target quantity and the drift velocity target quantity.
A computer readable storage medium having stored thereon computer program instructions which, when loaded and executed by a processor, cause the processor to perform the above-described method of estimating relative orbital change of a formation satellite.
Those skilled in the art will appreciate that those matters not described in detail in the present specification are well known in the art.
Although the present invention has been described with reference to the preferred embodiments, it is not intended to limit the present invention, and those skilled in the art can make variations and modifications of the present invention without departing from the spirit and scope of the present invention by using the methods and technical contents disclosed above.
Claims (10)
1. A method for estimating relative orbital change of formation satellites is characterized by comprising the following steps:
obtaining effective local satellite orbit parameters and target satellite orbit parameters injected into the ground in the period;
estimating relative flat track angular velocity, relative track eccentricity component and relative track eccentricity change rate by adopting a filtering algorithm;
calculating the relative star spacing, the relative semi-major axis and the relative eccentricity variation;
the relative semi-major axis target amount and the drift velocity target amount are calculated.
2. The estimation method according to claim 1, wherein it is determined whether the local satellite orbit parameter and the target satellite orbit parameter are received during the local period, and if the local satellite orbit parameter and the target satellite orbit parameter are received, it is further determined whether the local satellite orbit parameter and the target satellite orbit parameter are valid.
3. The estimation method according to claim 1, wherein the effective local star orbit parameters and target star orbit parameters of the surface injection at least comprise: the reference time of the relative orbit and the variation thereof, the relative star spacing of the two stars and the variation thereof, the relative orbit eccentricity component, the relative orbit semi-major axis of the orbit injection, and the average value of the reference time of the two stars orbit.
4. The estimation method according to claim 1, wherein the filtering algorithm is used to estimate the relative flat track angular velocity, the relative track eccentricity component, and the relative track eccentricity change rate, and the iteration process is as follows:
e xRM =e*cos(ω)-e T *cos(ω T )
e yRM =e*sin(ω)-e T *sin(ω T )
e xdotRM =k fedot (e xRM -e xRM0 )/Dt RM +(1-k fedot )e xdotRM0
e ydotRM =k fedot (e yRM -e yRM0 )/Dt RM +(1-k fedot )e ydotRM0
in the formula, k fedot For filter gain, e is eccentricity, ω is orbital argument of perigee, v RM As drift velocity, v dotRM As rate of change of drift velocity, e xRM And e yRM Is a counter railComponent of track eccentricity, e xdotRM And e ydotRM For the amount of change thereof,is a semi-major axis flat root,the orbital angular velocity, Re, the earth radius, μ 398600km 3 /s 2 Subscript 0 represents the state variable of the last beat, and subscript T represents the state variable of the target star, which is the gravity constant; t is t RM0 Reference time, Dt, relative to the track RM0 Is t RM0 Amount of change of (Dr) RM0 Is the relative star spacing of two stars, Dr RM Is Dr RM0 Amount of change of e xRM 、e yRM Are all relative orbital eccentricity components, Δ a kRM Relative orbit semi-major axis, t, for orbit injection midRM Is the average of the two-star orbit reference time.
5. The estimation method according to claim 4, characterized in that said calculation of the relative inter-satellite distance r is carried out RM Relative semi-major axis Deltaa RM Relative eccentricity change amount Δ e xRM And Δ e yRM The method comprises the following steps:
Δt RM =t-t RM0
Δe xRM =e xRM +e xdotRM Δt RM
Δe yRM =e yRM +e ydotRM Δt RM 。
7. An apparatus for estimating relative orbital change of a formation satellite, comprising:
the acquisition module is used for receiving effective local satellite orbit parameters and target satellite orbit parameters injected into the ground in the period;
the estimation module estimates the relative flat track angular velocity, the relative track eccentricity component and the relative track eccentricity change rate by adopting a filtering algorithm;
the first calculation module is used for calculating the relative star spacing, the relative semi-major axis and the relative eccentricity variation;
and the second calculation module is used for calculating the relative semi-major axis target quantity and the drift velocity target quantity.
8. The estimation apparatus as claimed in claim 7, wherein the obtaining module first determines whether the local satellite orbit parameter and the target satellite orbit parameter of the ground injection are received in the current period, and if so, further determines whether the local satellite orbit parameter and the target satellite orbit parameter are valid.
9. The estimation apparatus according to claim 7, wherein the surface injected effective and target star orbit parameters at least comprise: the reference time of the relative orbit and the variation thereof, the relative star spacing of the two stars and the variation thereof, the relative orbit eccentricity component, the relative orbit semi-major axis of the orbit injection, and the average value of the reference time of the two stars orbit.
10. A computer readable storage medium having stored thereon computer program instructions which, when loaded and executed by a processor, cause the processor to perform the method of any of claims 1 to 6.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210302728.4A CN114852375A (en) | 2022-03-24 | 2022-03-24 | Method and device for estimating relative orbit change of formation satellite |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210302728.4A CN114852375A (en) | 2022-03-24 | 2022-03-24 | Method and device for estimating relative orbit change of formation satellite |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114852375A true CN114852375A (en) | 2022-08-05 |
Family
ID=82629139
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210302728.4A Pending CN114852375A (en) | 2022-03-24 | 2022-03-24 | Method and device for estimating relative orbit change of formation satellite |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114852375A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117311394A (en) * | 2023-11-28 | 2023-12-29 | 中国西安卫星测控中心 | Inter-satellite phase overrun control method in non-strict multi-satellite serial formation system |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001153682A (en) * | 1999-11-26 | 2001-06-08 | Fujitsu Ltd | Orbit calculation device and orbit calculating method |
KR101294623B1 (en) * | 2012-12-28 | 2013-08-08 | 한국항공우주연구원 | Enhanced ground based precise attitude determination method of imaging satellite |
CN104142686A (en) * | 2014-07-16 | 2014-11-12 | 北京控制工程研究所 | Autonomous formation flight control method for satellites |
CN105138005A (en) * | 2015-06-16 | 2015-12-09 | 北京航空航天大学 | Method for determining relative orbit elements based on inter-satellite distance |
CN110595486A (en) * | 2019-09-05 | 2019-12-20 | 上海航天控制技术研究所 | High-precision semimajor axis deviation calculation method based on double-star on-orbit telemetry data |
CN110816896A (en) * | 2019-10-28 | 2020-02-21 | 中国空间技术研究院 | Satellite on-satellite simple orbit extrapolation method |
US20200319350A1 (en) * | 2019-04-02 | 2020-10-08 | Technology And Engineering Center For Space Utilization, Chinese Academy Of Sciences | Method for achieving space-based autonomous navigation of global navigation satellite system (gnss) satellites |
CN112208798A (en) * | 2020-10-22 | 2021-01-12 | 上海卫星工程研究所 | Flight-around formation high code rate inter-satellite link switching method and system |
CN113189619A (en) * | 2021-04-02 | 2021-07-30 | 中国空间技术研究院 | Low-orbit constellation phase retention parameter estimation method |
JP2021160635A (en) * | 2020-04-01 | 2021-10-11 | 三菱電機株式会社 | Orbit control device, ground station device and orbit control program |
CN113665849A (en) * | 2021-09-29 | 2021-11-19 | 长光卫星技术有限公司 | Autonomous phase control method combining EKF filtering algorithm and neural network |
US20220002006A1 (en) * | 2020-04-30 | 2022-01-06 | Space Engineering University | Low-orbit satellite deorbit control method and system based on particle swarm algorithm |
US20220017239A1 (en) * | 2020-07-17 | 2022-01-20 | The Aerospace Corporation | Methods and systems for orbit estimation of a satellite |
CN114002709A (en) * | 2021-10-20 | 2022-02-01 | 上海航天空间技术有限公司 | Singularity removing method suitable for recursive calculation of satellite orbit parameters |
-
2022
- 2022-03-24 CN CN202210302728.4A patent/CN114852375A/en active Pending
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2001153682A (en) * | 1999-11-26 | 2001-06-08 | Fujitsu Ltd | Orbit calculation device and orbit calculating method |
KR101294623B1 (en) * | 2012-12-28 | 2013-08-08 | 한국항공우주연구원 | Enhanced ground based precise attitude determination method of imaging satellite |
CN104142686A (en) * | 2014-07-16 | 2014-11-12 | 北京控制工程研究所 | Autonomous formation flight control method for satellites |
CN105138005A (en) * | 2015-06-16 | 2015-12-09 | 北京航空航天大学 | Method for determining relative orbit elements based on inter-satellite distance |
US20200319350A1 (en) * | 2019-04-02 | 2020-10-08 | Technology And Engineering Center For Space Utilization, Chinese Academy Of Sciences | Method for achieving space-based autonomous navigation of global navigation satellite system (gnss) satellites |
CN110595486A (en) * | 2019-09-05 | 2019-12-20 | 上海航天控制技术研究所 | High-precision semimajor axis deviation calculation method based on double-star on-orbit telemetry data |
CN110816896A (en) * | 2019-10-28 | 2020-02-21 | 中国空间技术研究院 | Satellite on-satellite simple orbit extrapolation method |
JP2021160635A (en) * | 2020-04-01 | 2021-10-11 | 三菱電機株式会社 | Orbit control device, ground station device and orbit control program |
US20220002006A1 (en) * | 2020-04-30 | 2022-01-06 | Space Engineering University | Low-orbit satellite deorbit control method and system based on particle swarm algorithm |
US20220017239A1 (en) * | 2020-07-17 | 2022-01-20 | The Aerospace Corporation | Methods and systems for orbit estimation of a satellite |
CN112208798A (en) * | 2020-10-22 | 2021-01-12 | 上海卫星工程研究所 | Flight-around formation high code rate inter-satellite link switching method and system |
CN113189619A (en) * | 2021-04-02 | 2021-07-30 | 中国空间技术研究院 | Low-orbit constellation phase retention parameter estimation method |
CN113665849A (en) * | 2021-09-29 | 2021-11-19 | 长光卫星技术有限公司 | Autonomous phase control method combining EKF filtering algorithm and neural network |
CN114002709A (en) * | 2021-10-20 | 2022-02-01 | 上海航天空间技术有限公司 | Singularity removing method suitable for recursive calculation of satellite orbit parameters |
Non-Patent Citations (3)
Title |
---|
刘其睿: "卫星姿轨控系统设计与分析平台软件方案及实现", 空间电子技术, 10 February 2016 (2016-02-10), pages 74 - 77 * |
杨盛庆;杜耀珂;贾艳胜;王文妍;: "基于约化相对轨道拟平根数的长期稳定高精度卫星编队导航技术", 空间控制技术与应用, no. 01, 15 February 2017 (2017-02-15), pages 30 - 35 * |
雪丹;曹喜滨;: "卫星编队相对轨道的分布式估计方法", 吉林大学学报(工学版), no. 06, 10 November 2006 (2006-11-10), pages 252 - 255 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117311394A (en) * | 2023-11-28 | 2023-12-29 | 中国西安卫星测控中心 | Inter-satellite phase overrun control method in non-strict multi-satellite serial formation system |
CN117311394B (en) * | 2023-11-28 | 2024-03-12 | 中国西安卫星测控中心 | Inter-satellite phase overrun control method in non-strict multi-satellite serial formation system |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105224737B (en) | A kind of first value correction method of extraterrestrial target improvement of orbit | |
CN107797130B (en) | Method for calculating uplink data of multi-point and multi-parameter orbit of low-orbit spacecraft | |
Canuto et al. | All-propulsion design of the drag-free and attitude control of the European satellite GOCE | |
Li et al. | A machine learning-based approach for improved orbit predictions of LEO space debris with sparse tracking data from a single station | |
CN110378012B (en) | Strict regression orbit design method, system and medium considering high-order gravity field | |
CN102495950B (en) | Acquisition method of inclination biased quantity for sun synchronous orbit | |
CN108875244B (en) | Orbit prediction precision improvement method based on random forest | |
CN115143955B (en) | Method for determining initial orbit of geosynchronous orbit with spacecraft based on astronomical angle measurement data | |
Li et al. | Innovative Mars entry integrated navigation using modified multiple model adaptive estimation | |
CN113761809B (en) | Passive detection orbit determination method based on deep neural network | |
CN114852375A (en) | Method and device for estimating relative orbit change of formation satellite | |
CN111731513B (en) | Method for maintaining regression orbit in high-precision gravitational field based on monopulse orbit control | |
CN107246883A (en) | A kind of Rotating Platform for High Precision Star Sensor installs the in-orbit real-time calibration method of matrix | |
CN108490973B (en) | Method and device for determining relative orbit of spacecraft formation | |
CN110203424B (en) | Method and device for estimating spacecraft spin motion by using speed measurement data | |
CN106959457A (en) | A kind of GLONASS almanac parameters methods of estimation for satellite navigation | |
CN108516107A (en) | A kind of online motor power and specific impulse method of estimation and system | |
CN109682383B (en) | Real-time filtering positioning method for measuring distance and data by using deep space three-way | |
CN112797988B (en) | Unknown maneuvering spacecraft orbit determination method based on neural network | |
CN113375697A (en) | X-ray pulsar signal time delay estimation method based on maximum posterior estimation | |
CN116401833A (en) | Satellite orbit calculation method based on SGP4 model | |
CN113128032B (en) | Intersection point time and position solving algorithm based on orbit analysis perturbation solution | |
Sabzevari et al. | Observability analysis and design of two nested filters for the satellite attitude estimation with magnetometer‐only | |
CN104713561B (en) | A kind of lunar orbiter precise orbit determination method | |
CN114002710A (en) | On-satellite orbit position autonomous prediction method for small-eccentricity low-orbit satellite |
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 |