CN112257343B - High-precision ground track repetitive track optimization method and system - Google Patents
High-precision ground track repetitive track optimization method and system Download PDFInfo
- Publication number
- CN112257343B CN112257343B CN202011138645.3A CN202011138645A CN112257343B CN 112257343 B CN112257343 B CN 112257343B CN 202011138645 A CN202011138645 A CN 202011138645A CN 112257343 B CN112257343 B CN 112257343B
- Authority
- CN
- China
- Prior art keywords
- track
- satellite
- orbit
- earth
- correction
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 58
- 238000005457 optimization Methods 0.000 title claims abstract description 37
- 230000003252 repetitive effect Effects 0.000 title claims abstract description 20
- 238000012937 correction Methods 0.000 claims abstract description 67
- 230000005484 gravity Effects 0.000 claims abstract description 67
- 238000004458 analytical method Methods 0.000 claims abstract description 13
- 238000004364 calculation method Methods 0.000 claims abstract description 12
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 8
- 238000013461 design Methods 0.000 claims description 14
- 230000001174 ascending effect Effects 0.000 claims description 9
- 230000000630 rising effect Effects 0.000 claims description 8
- 230000010354 integration Effects 0.000 claims description 7
- 238000004088 simulation Methods 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 5
- 238000006073 displacement reaction Methods 0.000 claims description 5
- 230000014759 maintenance of location Effects 0.000 claims description 5
- 238000012360 testing method Methods 0.000 claims description 4
- 229910002056 binary alloy Inorganic materials 0.000 claims description 3
- 230000009191 jumping Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 4
- 230000003044 adaptive effect Effects 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 230000008014 freezing Effects 0.000 description 2
- 238000007710 freezing Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/27—Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/004—Artificial life, i.e. computing arrangements simulating life
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Software Systems (AREA)
- Artificial Intelligence (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computational Linguistics (AREA)
- Biophysics (AREA)
- Biomedical Technology (AREA)
- Data Mining & Analysis (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Medical Informatics (AREA)
- Computer Hardware Design (AREA)
- Geometry (AREA)
- Navigation (AREA)
- Position Fixing By Use Of Radio Waves (AREA)
Abstract
The invention provides a high-precision ground track repetitive orbit optimization method and system, wherein a corresponding satellite orbit dynamics recursion model is established according to the requirement analysis of an earth gravity field non-spherical perturbation order; the method comprises the steps of providing heavy rail interference track parameters under a J2 low-order gravity field according to an analytical formula, and obtaining correction parameters under J4 perturbation through iterative calculation to serve as initial values of a high-precision ground track repeated track parameter self-adaptive setting algorithm; and according to the requirement of the repeated track precision, carrying out self-adaptive optimization setting on the repeated track parameters of the high-precision ground track. The method solves the problems that the traditional low-order gravity field orbit model is low in ground track repetition precision, the high-order gravity field orbit model is high in nonlinearity, many in iterative correction parameters, incapable of resolving and the like, and has more general engineering practicability.
Description
Technical Field
The invention relates to astronavigation aircraft orbital dynamics, in particular to a high-precision ground track repeated orbit optimization method and system.
Background
The revisiting characteristic of the track of the sub-satellite points is an important index for the design of the earth observation satellite track, the requirement of the traditional satellite ground track on the return precision of the track is not high, the track deviation can be from several kilometers to dozens of kilometers, and the flying task of the satellite is not influenced. In recent years, with the requirement of imaging precision greatly improved, the limitation of the existing regression orbit design method is increasingly obvious, and particularly for a surveying and mapping task, in order to realize accurate surveying and mapping and efficient task planning, high-precision regression of a satellite reference orbit ground track needs to be ensured.
At present, a plurality of regression orbit design schemes are applied at home and abroad, but most of the regression orbit design schemes are based on an analysis method of a low-order gravity field model, and the regression precision deviation is as high as 3km to 10km. In the Chinese patent "a method for determining a strict regression orbit of a near earth satellite" (CN 106092105A), the popchengqing, ducuke, wang etime and the like perform iterative correction by establishing a relationship between a semi-major axis of the orbit, an inclination angle of the orbit and a longitude and latitude of a sub-satellite by taking regression accuracy as an index; and aiming at the characteristic of the eccentricity ratio vector limit cycle, repeatedly carrying out iterative correction on the eccentricity ratio and the perigee argument by adopting an averaging method. In the method, the perturbation of the earth 2-order gravity field is only considered in the functional relation between the satellite points for iterative correction and the orbit semi-major axis and inclination angle, so that the regression precision is influenced; each iteration correction of eccentricity and argument of near place needs to acquire 4-month period data for averaging, the calculated amount is large, and the iteration convergence speed is influenced; the regression precision of the obtained orbit can only reach the meter level, and engineering application is influenced.
In summary, a high-precision ground track repetitive track design under the influence of high-order gravity field perturbation needs to be developed for centimeter-level regression precision requirements.
Disclosure of Invention
Aiming at the defects in the prior art, the invention aims to provide a high-precision ground track repetitive track optimization method and system.
The invention provides a high-precision ground track repetitive track optimization method, which comprises the following steps:
step A: according to the requirement analysis of the earth gravity field non-spherical perturbation order, establishing a corresponding satellite orbit dynamics recursion model;
and B: the method comprises the steps of providing heavy rail interference track parameters under a J2 low-order gravity field according to an analytical formula, and obtaining correction parameters under J4 perturbation through iterative calculation to serve as initial values of a high-precision ground track repeated track parameter self-adaptive setting algorithm;
step C: and according to the requirement of the repeated track precision, carrying out self-adaptive optimization setting on the repeated track parameters of the high-precision ground track.
Preferably, the precision of the heavy rail track is less than 0.01m.
Preferably, the step a includes:
taking a numerical simulation test result as a basis, comprehensively considering ground track repetition precision and track recursion calculation amount, establishing a satellite track dynamic model corresponding to the following steps by adopting an expression method of an earth gravity field potential function:
wherein,is the earth-centered position vector of the satellite;is the function of the field potential of the earth gravity;
the gravity field potential function comprises two parts of earth central gravity and earth non-spherical gravity, and if the earth is considered as a rigid body and an equatorial plane is superposed with a basic plane of an epoch inertial system, the gravity field potential function is expanded into a series form in the earth central inertial system as follows:
wherein mu is the gravitational constant of the earth, r is the earth center distance of the satellite, C n0 And C nm And S nm All are spherical harmonic coefficients, R e Is the equatorial radius,The geocentric latitude and the lambda are geocentric longitude which are obtained by calculating a position vector in a satellite geostationary system;andthe terms are Legendre and associated Legendre polynomials respectively, are correction parts of a real earth gravitational potential pair uniform sphere and comprise a band harmonic term and a field harmonic term;
and obtaining coefficient values of each order of a gravitational field potential function through an earth gravitational field table, calculating to obtain a perturbation force expression under a satellite geostationary system, and further integrating to obtain the position and the speed of the satellite.
Preferably, the step B includes:
step S4.1: according to the overall design constraint, giving a semi-major axis analytic solution a of the heavy rail interference track under the J2 low-order gravity field J2 Orbit dip angle analytic solution i J2 A heavy rail period T and a rail turn number Q;
according to the load working power, the satellite weight, the carrying and launching capacity and other overall constraints, the semi-major axis a of the heavy-rail interference orbit can be determined J2 ;
Orbit dip analytic solution i J2 :
Track number of turns Q, heavy rail period T:
T=T n ×Q
wherein, T n Is a point of intersectionPeriod, w e Is the rotational angular velocity of the earth;
step S4.2: according to the J2 perturbation analysis orbit, a J4 perturbation correction orbit based on a Newton iteration method is provided, and the method comprises the following steps:
step S4.2.1: selecting t according to the ground track arrangement requirement of the satellite task 0 Latitude and longitude lambda of time sub-satellite point 0 、Establishing a functional relation of longitude and latitude about a semi-long axis and an inclination angle after the following track is repeated;
wherein,
M=n(t-t 0 )
G 0 is t 0 At time Greenwich mean sidereal time omega 0 Is t 0 The right ascension at the ascending crossing point of the moment,is J 4 Perturbation of the rising point right ascension drift rate;
step S4.2.2: obtaining the following heavy rail interference track eccentricity e by utilizing the constraint relation of the frozen track J3 And estimating the initial value of the argument w of the near place:
step S4.2.3: first with the satellite t 0 Time ground track initial longitude and latitude lambda 0 、Taking the starting point as the point of the second track, performing the integration of the second track period T to obtain T 0 Satellite ground track longitude and latitude lambda at + T moment n 、Then calculating to obtain the longitude and latitude deviation delta lambda = lambda of the initial end and the final end of the heavy rail n -λ 0 ,Judging whether the deviation meets the heavy rail interference precision index or not; if the index is satisfied, outputting a semimajor axis and a tilt angle correction value a J4 、i J4 If the iteration parameter is out of tolerance, the following parameter correction is carried out, and the integral judgment process is converted back again until the iteration parameter meets the precision index;
a k =a k-1 +Δa k-1
i k =i k-1 +Δi k-1
wherein,
a k for the semi-major axis of the orbit at time k, a k-1 Is the half-major axis of the track at time k-1, i k For the track inclination at time k, i k-1 For track inclination at time k-1;
Step S4.2.4: calculating to obtain t by using the functional relationship between the semi-major axis and the inclination angle of the track in the step S4.2.1 and the longitude and latitude of the subsatellite point 0 Time a J4 、i J4 Corresponding rising point right ascension omega J4 Angle M close to the mean J4 。
Preferably, the step C includes:
step S5.1: the method comprises the steps of converting a track repeated nonlinear parameter solving problem under a high-order gravity field into a multivariable and multi-target optimization problem, guiding parameter self-adaptive setting by combining target track characteristic information, and describing an optimization model as follows:
optimizing variables: [ Δ a, Δ e, Δ i, Δ Ω, Δ w, Δ M ]
Initial conditions were as follows: [ a ] A J4 ,e J3 ,i J4 ,Ω J4 ,w=90°,M J4 ]
Constraint conditions are as follows:
delta a is correction of track semimajor axis, delta e is correction of track eccentricity, delta i is correction of track inclination angle, delta omega is correction of right ascension at track ascending intersection point, delta w is correction of amplitude at track perigee, delta M is correction of angle at track mean perigee,is a vector of the acceleration of the satellite,in order to be a satellite position deviation,for the satellite position vector at time T,is t 0 A time satellite position vector;
step S5.2: selecting individuals through a binary system tournament method, and performing crossing and variation to generate a new population;
step S5.3: calculating and updating a new population objective function value, namely the position and speed deviation of the satellite earth-fixed system;
step S5.4: generating a new combined population by a merging method, and carrying out non-dominant sorting;
step S5.5: selecting individuals to form a new generation of population through a displacement and elite retention strategy;
step S5.6: and skipping to the step S5.2, and circularly updating until the termination condition is met.
The invention provides a high-precision ground track repetitive track optimization system, which comprises:
a module A: according to the demand analysis of the earth gravity field non-spherical perturbation order, establishing a corresponding satellite orbit dynamics recursion model;
and a module B: according to an analytical formula, heavy-rail interference track parameters under a J2 low-order gravity field are given, correction parameters under J4 perturbation are obtained through iterative calculation and serve as initial values of a high-precision ground track repeated track parameter self-adaptive setting algorithm;
and a module C: and performing self-adaptive optimization setting of the high-precision ground track repeated track parameters according to the repeated track precision requirement.
Preferably, the precision of the heavy rail track is less than 0.01m.
Preferably, the module a comprises:
taking the numerical simulation test result as a basis, comprehensively considering the ground track repetition precision and the orbit recursion computation amount, establishing a satellite orbit dynamics model corresponding to the following steps, and adopting an expression method of an earth gravity field potential function:
wherein,is the earth-centered position vector of the satellite;is a function of the field potential of the earth gravity;
the gravity field potential function comprises two parts of earth central gravity and earth non-spherical gravity, and if the earth is considered as a rigid body and an equatorial plane is superposed with a basic plane of an epoch inertial system, the gravity field potential function is expanded into a series form in the earth central inertial system as follows:
wherein mu is the gravitational constant of the earth, r is the earth center distance of the satellite, C n0 And C nm And S nm All are spherical harmonic coefficients, R e Is the equatorial radius,The geocentric latitude and the lambda are geocentric longitude which are obtained by calculating a position vector in a satellite geostationary system;andthe terms are Legendre and associated Legendre polynomials respectively, are correction parts of a real earth gravitational potential pair uniform sphere and comprise a band harmonic term and a field harmonic term;
and obtaining coefficient values of each order of a gravitational field potential function through an earth gravitational field table, calculating to obtain a perturbation force expression under a satellite geostationary system, and further integrating to obtain the position and the speed of the satellite.
Preferably, the module B comprises:
module S4.1: according to the overall design constraint, giving a semi-major axis analytic solution a of the heavy rail interference track under the J2 low-order gravity field J2 Orbit dip angle analytic solution i J2 And heavy track period T, trackThe number of turns Q;
according to the load working power, the satellite weight, the carrying and launching capacity and other overall constraints, the semi-major axis a of the heavy-rail interference orbit can be determined J2 ;
Orbit dip analytic solution i J2 :
Track turn number Q, heavy track period T:
T=T n ×Q
wherein, T n Is a period of intersection, w e Is the rotational angular velocity of the earth;
module S4.2: according to the J2 perturbation analysis orbit, a J4 perturbation correction orbit based on a Newton iteration method is provided, and the method comprises the following steps:
module s4.2.1: selecting t according to the ground track arrangement requirement of the satellite task 0 Latitude and longitude lambda of time point under star 0 、Establishing a functional relation of longitude and latitude about a semi-major axis and an inclination angle after the following tracks are repeated;
wherein,
M=n(t-t 0 )
G 0 is t 0 At time Greenwich mean sidereal time omega 0 Is t 0 The right ascension at the ascending crossing point of the moment,is J 4 Perturbation of the rising point right ascension drift rate;
module s4.2.2: obtaining the following heavy rail interference track eccentricity e by utilizing the constraint relation of the frozen track J3 And estimating the initial value of the argument w of the near place:
module s4.2.3: firstly, with the satellite t 0 Time ground track initial longitude and latitude lambda 0 、Taking the starting point as the point of integration of the heavy track period T to obtain T 0 Satellite ground track longitude and latitude lambda at + T moment n 、Then calculating to obtain the longitude and latitude deviation delta lambda = lambda of the initial end of the heavy rail n -λ 0 ,Judging whether the deviation meets the heavy rail interference precision index or not; if the index is satisfied, outputting a semimajor axis and inclination angle correction value a J4 、i J4 If the iteration parameter is out of tolerance, the following parameter correction is carried out, and the integral judgment process is converted back again until the iteration parameter meets the precision index;
a k =a k-1 +Δa k-1
i k =i k-1 +Δi k-1
wherein,
a k for the semi-major axis of the orbit at time k, a k-1 Is the half-major axis of the track at time k-1, i k For the track inclination at time k, i k-1 Is the inclination angle of the orbit at the moment of k-1;
module s4.2.4: t is obtained by calculation according to the function relation between the semi-major axis and the inclination angle of the orbit in the module S4.2.1 and the longitude and latitude of the satellite point 0 Time a J4 、i J4 Corresponding rising point right ascension omega J4 Angle M close to the mean J4 。
Preferably, the module C comprises:
module S5.1: the method comprises the following steps of converting a track repetition nonlinear parameter solving problem under a high-order gravity field into a multivariable and multi-target optimization problem, guiding parameter self-adaptive setting by combining target track characteristic information, and describing an optimization model as follows:
optimizing variables: [ Δ a, Δ e, Δ i, Δ Ω, Δ w, Δ M ]
Initial conditions were as follows: [ a ] A J4 ,e J3 ,i J4 ,Ω J4 ,w=90°,M J4 ]
Constraint conditions are as follows:
delta a is correction of track semimajor axis, delta e is correction of track eccentricity, delta i is correction of track inclination angle, delta omega is correction of right ascension at track ascending intersection point, delta w is correction of amplitude at track perigee, delta M is correction of angle at track mean perigee,as a vector of the acceleration of the satellite,in order to be a satellite position deviation,for the satellite position vector at time T,is t 0 A time satellite position vector;
module S5.2: selecting individuals through a binary tournament method, and performing crossing and variation to generate a new population;
module S5.3: calculating and updating a new population objective function value, namely the position and speed deviation of the satellite earth-fixed system;
module S5.4: generating a new combined population by a merging method, and carrying out non-dominant sorting;
module S5.5: selecting individuals to form a new generation of population through a displacement and elite retention strategy;
module S5.6: and jumping to a module S5.2, and circularly updating until a termination condition is met.
Compared with the prior art, the invention has the following beneficial effects:
the method solves the problems that the traditional low-order gravity field orbit model is low in ground track repetition precision, the high-order gravity field orbit model is high in nonlinearity, many in iterative correction parameters, incapable of resolving and the like, and has more general engineering practicability.
Drawings
Other features, objects and advantages of the invention will become more apparent upon reading of the detailed description of non-limiting embodiments with reference to the following drawings:
FIG. 1 is a schematic block diagram of a high-precision ground track repetitive orbit optimization method;
FIG. 2 is a schematic block diagram of iterative correction of an initial value of a high-precision ground track repetitive orbit J4;
FIG. 3 is a schematic diagram of high-precision ground track repetitive orbit parameter adaptive setting optimization;
FIG. 4 is a high-precision ground track repetitive orbit high-order gravity field model order determination simulation diagram.
Detailed Description
The present invention will be described in detail with reference to specific examples. The following examples will aid those skilled in the art in further understanding the present invention, but are not intended to limit the invention in any manner. It should be noted that variations and modifications can be made by persons skilled in the art without departing from the concept of the invention. All falling within the scope of the present invention.
As shown in the attached figure 1, the invention provides a high-precision ground track repetitive track optimization method with self-adaptive parameter setting under the influence of high-order gravity field perturbation. The method specifically comprises the following steps:
step A: according to the requirement analysis of the earth gravity field non-spherical perturbation order, establishing a corresponding satellite orbit dynamics recursion model;
and B, step B: giving out heavy-rail interference track parameters under a J2 low-order gravity field according to an analytical formula, and iteratively calculating on the basis to obtain correction parameters under J4 perturbation as initial values of a high-precision ground track repeated track parameter self-adaptive setting algorithm;
and C: and according to the requirement of the repeated track precision, carrying out self-adaptive optimization setting on the repeated track parameters of the high-precision ground track.
In this embodiment, the regression accuracy is better than 0.01m.
In the step A: as shown in fig. 4, the orbit repetition position precision of the regression orbit under each other low-order perturbation model is calculated by a numerical simulation test method with 120 × 120-order earth gravity field as a standard, and then the orbit recursion computation amount is considered comprehensively, so that a suggestion that the order of the gravity field model is selected to be 90 × 90 is given, and the following corresponding satellite orbit dynamics model is established.
In order to facilitate the description of the numerical integration of the satellite orbit in the gravity field, a representation method of the potential function of the earth gravity field is adopted, namely:
wherein,is the earth-centered position vector of the satellite;is a function of the field potential of the earth's gravity.
The gravity field potential function includes two parts of earth central gravity and earth non-spherical gravity, and if the earth is considered as a rigid body and an equatorial plane is superposed with a basic plane of an epoch inertial system, the gravity field potential function can be expanded into a series form in the earth central inertial system as follows:
wherein mu is the gravitational constant of the earth, r is the earth center distance of the satellite, C n0 And C nm And S nm All are spherical harmonic coefficients, R e Is the equatorial radius,Latitude of geocentric, λ of geocentricLongitude, both calculated from the position vector in the satellite earth's system;andthe terms are Legendre and associated Legendre polynomials respectively, are correction parts of a real earth gravitational potential pair uniform sphere, comprise a band harmonic term and a field harmonic term, and reflect the unevenness of the earth.
The coefficient values of each order of the gravitational field potential function are obtained through the earth gravitational field table, the perturbation force expression under the satellite earth-fixed system can be obtained through calculation, and then the satellite position and the satellite speed are obtained through integration.
Preferably, the step B comprises the following steps:
step S4.1: according to the overall design constraint, giving a semi-major axis analytic solution a of the heavy rail interference track under the J2 low-order gravity field J2 Orbit dip angle analytic solution i J2 And a heavy track period T, a track turn number Q.
According to the load working power, the satellite weight, the carrying and launching capacity and other overall constraints, the semi-major axis a of the heavy-rail interference orbit can be determined J2 ;
Orbit dip analytic solution i J2 :
Track turn number Q, heavy track period T:
T=T n ×Q
wherein, T n Is a period of intersection, w e Is the rotational angular velocity of the earth.
Step S4.2: and (4) according to the J2 perturbation analysis orbit, providing a J4 perturbation correction orbit based on a Newton iteration method. The research idea of the method is as follows:
step S4.2.1: selecting t according to the requirement of satellite task on ground track arrangement 0 Latitude and longitude lambda of time sub-satellite point 0 、And establishing the following functional relationship of longitude and latitude of the repeated track with respect to the semi-major axis and the inclination angle.
Wherein,
M=n(t-t 0 )
G 0 is t 0 At time Greenwich mean sidereal time omega 0 Is t 0 The right ascension at the ascending crossing point of the moment,is J 4 Ascension point right ascension drift rate in perturbation.
Step S4.2.2: and obtaining the following heavy rail interference track eccentricity and near-place argument initial value estimation by utilizing the constraint relation of the frozen track.
Step S4.2.3: as shown in the calculation flow of FIG. 2, firstly, the satellite t is used 0 Time ground track initial longitude and latitude lambda 0 、Taking the starting point as the point of integration of the heavy track period T to obtain T 0 Satellite ground track longitude and latitude lambda at + T moment n 、Then calculating to obtain the longitude and latitude deviation delta lambda = lambda of the initial end and the final end of the heavy rail n -λ 0 ,Judging whether the deviation meets the heavy rail interference precision index or not; if the index is satisfied, outputting a semimajor axis and a tilt angle correction value a J4 、i J4 And if the iteration parameter is out of tolerance, the following parameter correction is carried out, and the integral judgment process is converted back until the iteration parameter meets the precision index.
a k =a k-1 +Δa k-1
i k =i k-1 +Δi k-1
Wherein,
a k for the semi-major axis of the orbit at time k, a k-1 Is the semimajor axis of the track at time k-1, i k For the track inclination at time k, i k-1 Is the inclination angle of the orbit at the moment of k-1;
step S4.2.4: by utilizing the functional relationship between the semi-major axis and the inclination angle of the orbit and the longitude and latitude of the satellite point in the step S4.2.1, t can be obtained by calculation 0 Time a J4 、i J4 Corresponding rising point right ascension omega J4 Angle M close to the mean J4 。
As shown in fig. 3, the step C includes the following steps:
step S5.1: the problem of solving the track repeated nonlinear parameters under the high-order gravity field is converted into a multivariable and multi-target optimization problem, and the parameters are guided to be self-adaptively adjusted by combining target track characteristic information (freezing performance and regression period). The optimization model is described below.
optimizing variables: [ Δ a, Δ e, Δ i, Δ Ω, Δ w, Δ M ]
Initial conditions: [ a ] A J4 ,e J3 ,i J4 ,Ω J4 ,w=90°,M J4 ]
Constraint conditions are as follows:
delta a is correction of track semimajor axis, delta e is correction of track eccentricity, delta i is correction of track inclination angle, delta omega is correction of right ascension at track ascending intersection point, delta w is correction of amplitude at track perigee, delta M is correction of angle at track mean perigee,is a vector of the acceleration of the satellite,in order to be a satellite position deviation,for the satellite position vector at time T,is t 0 A time satellite position vector;
step S5.2: selecting individuals through a binary tournament method, and performing crossing and variation to generate a new population;
step S5.3: calculating and updating a new population target function value, namely satellite geostationary system position and speed deviation;
step S5.4: generating a new combined population by a merging method, and carrying out non-dominated sorting;
step S5.5: selecting individuals to form a new generation of population through a displacement and elite retention strategy;
step S5.6: and skipping to the step S5.2, and circularly updating until the termination condition is met.
In this embodiment, the design input is the orbit height 6989.90km, and the satellite initial starting point longitude and latitude are 0.135 ° S and 90.019 ° W. The 90 × 90-order earth gravity field model of the EGM2008 model is selected for orbit recursion, and the initial epoch is 0 minute 0 second (UTCG) at 6 months, 1 day, 12 hours, 2023 years.
As shown in table 1, firstly, according to the design constraint, obtaining an analytic solution and regression characteristics under the perturbation of the regression trajectory J2 according to step a; then obtaining a correction value under J4 perturbation on the basis according to the step B, and taking the correction value as an initial value of high-precision ground track repetitive orbit parameter adaptive optimization design; and finally, according to the step C, combining the track freezing characteristic and the regression characteristic determined in the initial value solving process to perform parameter self-adaptive optimization of the high-precision ground track repeated track until the regression precision meets centimeter-level design requirements (shown in a table 2), and outputting a group of high-precision ground track repeated track parameters.
TABLE 1 high-precision ground track repetitive orbit initial value and high-order optimization solution
TABLE 2 ground track position repeat accuracy
The invention provides a high-precision ground track repetitive track optimization system, which comprises:
a module A: and establishing a corresponding satellite orbit dynamics recursion model according to the requirement analysis of the earth gravity field non-spherical perturbation order.
And a module B: and (3) giving out heavy rail interference track parameters under the J2 low-order gravity field according to an analytical formula, and performing iterative calculation to obtain correction parameters under J4 perturbation, wherein the correction parameters are used as initial values of a high-precision ground track repeated track parameter adaptive setting algorithm.
And a module C: and according to the requirement of the repeated track precision, carrying out self-adaptive optimization setting on the target track parameter of the high-precision ground track repeated track parameter self-adaptive setting algorithm.
Those skilled in the art will appreciate that, in addition to implementing the system and its various devices, modules, units provided by the present invention as pure computer readable program code, the system and its various devices, modules, units provided by the present invention can be fully implemented by logically programming method steps in the form of logic gates, switches, application specific integrated circuits, programmable logic controllers, embedded microcontrollers and the like. Therefore, the system and various devices, modules and units thereof provided by the invention can be regarded as a hardware component, and the devices, modules and units included in the system for realizing various functions can also be regarded as structures in the hardware component; means, modules, units for realizing various functions can also be regarded as structures in both software modules and hardware components for realizing the methods.
The foregoing description of specific embodiments of the present invention has been presented. It is to be understood that the present invention is not limited to the specific embodiments described above, and that various changes or modifications may be made by one skilled in the art within the scope of the appended claims without departing from the spirit of the invention. The embodiments and features of the embodiments of the present application may be combined with each other arbitrarily without conflict.
Claims (6)
1. A high-precision ground track repetitive track optimization method is characterized by comprising the following steps:
step A: according to the demand analysis of the earth gravity field non-spherical perturbation order, establishing a corresponding satellite orbit dynamics recursion model;
and B, step B: the method comprises the steps of providing heavy rail interference track parameters under a J2 low-order gravity field according to an analytical formula, and obtaining correction parameters under J4 perturbation through iterative calculation to serve as initial values of a high-precision ground track repeated track parameter self-adaptive setting algorithm;
and C: according to the requirement of the repeated track precision, carrying out self-adaptive optimization setting on the repeated track parameters of the high-precision ground track;
the step A comprises the following steps:
taking the numerical simulation test result as a basis, comprehensively considering the ground track repetition precision and the orbit recursion computation amount, establishing a satellite orbit dynamics model corresponding to the following steps, and adopting an expression method of an earth gravity field potential function:
wherein,is the earth-centered position vector of the satellite;is the function of the field potential of the earth gravity;
the gravity field potential function comprises two parts of earth central gravity and earth non-spherical gravity, and if the earth is considered as a rigid body and an equatorial plane is superposed with a basic plane of an epoch inertial system, the gravity field potential function is expanded into a series form in the earth central inertial system as follows:
wherein mu is the gravitational constant of the earth, r is the earth center distance of the satellite, C n0 And C nm And S nm All are spherical harmonic coefficients, R e Is the equatorial radius,The geocentric latitude and the lambda are geocentric longitude which are obtained by calculating a position vector in a satellite geostationary system;andthe terms are Legendre and associated Legendre polynomials respectively, are correction parts of a real earth gravitational potential pair uniform sphere and comprise a band harmonic term and a field harmonic term;
acquiring each order coefficient value of a gravitational field potential function through an earth gravitational field table, calculating to obtain a perturbation force expression under a satellite geostationary system, and further integrating to obtain the position and the speed of the satellite;
the step B comprises the following steps:
step S4.1: according to the overall design constraint, giving a semi-major axis analytic solution a of the heavy rail interference track under the J2 low-order gravity field J2 And resolving the track inclination angle i J2 A heavy rail period T and a rail turn number Q;
according to load working power, satellite weight and carrying and launching capacity total constraints, a heavy-orbit interference orbit semi-major axis analytical solution a can be determined J2 ;
Orbit dip analytic solution i J2 :
Track number of turns Q, heavy rail period T:
T=T n ×Q
wherein, T n Is a period of intersection, w e Is the angular velocity of rotation of the earth i J2 Resolving the orbit inclination angle;
step S4.2: according to the J2 perturbation analysis orbit, a J4 perturbation correction orbit based on a Newton iteration method is provided, and the method comprises the following steps:
step S4.2.1: selecting t according to the ground track arrangement requirement of the satellite task 0 Latitude and longitude lambda of time point under star 0 、Establishing a functional relation of longitude and latitude about a semi-long axis and an inclination angle after the following track is repeated;
wherein,
M=n(t-t 0 )
G 0 is t 0 At time Greenwich mean sidereal time omega 0 Is t 0 The right ascension point of the moment,is the rising point right ascension drift rate under J4 perturbation;
step S4.2.2: obtaining the following heavy rail interference track eccentricity e by utilizing the constraint relation of the frozen track J3 And estimating the initial value of the argument w of the near place:
step S4.2.3: first with the satellite t 0 Time ground track initial longitude and latitude lambda 0 、Taking the starting point as the point of integration of the heavy track period T to obtain T 0 Satellite ground track longitude and latitude lambda at + T moment n 、Then calculating to obtain the longitude and latitude deviation delta lambda = lambda of the initial end of the heavy rail n -λ 0 ,Judging whether the deviation meets the heavy rail interference precision index or not; if the index is satisfied, outputting a semimajor axis and a tilt angle correction value a J4 、i J4 If the iteration parameter is out of tolerance, the following parameter correction is carried out, and the integral judgment process is converted back again until the iteration parameter meets the precision index;
a k =a k-1 +Δa k-1
i k =i k-1 +Δi k-1
wherein,
a k for the semi-major axis of the orbit at time k, a k-1 Is the semimajor axis of the track at time k-1, i k For the track inclination at time k, i k-1 Is the inclination angle of the orbit at the moment of k-1;
step S4.2.4: using the function relationship between the semi-major axis and the inclination angle of the orbit and the longitude and latitude of the satellite point in the step S4.2.1 to calculate and obtain t 0 Time a J4 、i J4 Corresponding rising point right ascension omega J4 Angle M close to the mean J4 。
2. A high accuracy ground track repeat trajectory optimization method as claimed in claim 1, wherein said repeat trajectory accuracy is less than 0.01m.
3. The high accuracy ground track repetitive orbit optimization method of claim 1, wherein the step C comprises:
step S5.1: the method comprises the steps of converting a track repeated nonlinear parameter solving problem under a high-order gravity field into a multivariable and multi-target optimization problem, guiding parameter self-adaptive setting by combining target track characteristic information, and describing an optimization model as follows:
optimizing variables: [ Δ a, Δ e, Δ i, Δ Ω, Δ w, Δ M ]
Initial conditions: [ a ] A J4 ,e J3 ,i J4 ,Ω J4 ,w=90°,M J4 ]
Constraint conditions are as follows:
delta a is correction of track semimajor axis, delta e is correction of track eccentricity, delta i is correction of track inclination angle, delta omega is correction of right ascension at track ascending intersection point, delta w is correction of amplitude at track perigee, delta M is correction of angle at track mean perigee,as a vector of the acceleration of the satellite,in order to be a satellite position deviation,for the satellite position vector at time T,is t 0 A time satellite position vector;
step S5.2: selecting individuals through a binary system tournament method, and performing crossing and variation to generate a new population;
step S5.3: calculating and updating a new population objective function value, namely the position and speed deviation of the satellite earth-fixed system;
step S5.4: generating a new combined population by a merging method, and carrying out non-dominated sorting;
step S5.5: selecting individuals to form a new generation of population through a displacement and elite retention strategy;
step S5.6: and jumping to the step S5.2, and circularly updating until the termination condition is met.
4. A high-precision ground track repeat trajectory optimization system, comprising:
a module A: according to the requirement analysis of the earth gravity field non-spherical perturbation order, establishing a corresponding satellite orbit dynamics recursion model;
and a module B: according to an analytical formula, heavy-rail interference track parameters under a J2 low-order gravity field are given, correction parameters under J4 perturbation are obtained through iterative calculation and serve as initial values of a high-precision ground track repeated track parameter self-adaptive setting algorithm;
and a module C: according to the requirement of the repeated track precision, carrying out self-adaptive optimization setting on the repeated track parameters of the high-precision ground track;
the module A comprises:
taking the numerical simulation test result as a basis, comprehensively considering the ground track repetition precision and the orbit recursion computation amount, establishing a satellite orbit dynamics model corresponding to the following steps, and adopting an expression method of an earth gravity field potential function:
wherein,is the earth-centered position vector of the satellite;is the function of the field potential of the earth gravity;
the gravity field potential function comprises two parts of earth central gravity and earth non-spherical gravity, and if the earth is considered as a rigid body and an equatorial plane is superposed with a basic plane of an epoch inertial system, the gravity field potential function is expanded into a series form in the earth central inertial system as follows:
wherein mu is the gravitational constant of the earth, r is the earth center distance of the satellite, C n0 And C nm And S nm All are spherical harmonic coefficients, R e Is the equatorial radius,The geocentric latitude and the lambda are geocentric longitude which are obtained by calculating a position vector in a satellite geostationary system;andlegendre and associated Legendre polynomials are respectively a correction part of a real earth gravitational potential to a uniform sphere, and comprise a band harmonic term and a field harmonic term;
obtaining each order coefficient value of a gravitational field potential function through an earth gravitational field table, calculating to obtain a perturbation force expression under a satellite geostationary system, and further integrating to obtain the position and the speed of the satellite;
the module B comprises:
module S4.1: according to the overall design constraint, giving a semi-major axis analytic solution a of the heavy rail interference track under the J2 low-order gravity field J2 And resolving the track inclination angle i J2 A heavy rail period T and a rail turn number Q;
according to load working power, satellite weight and total constraint of carrying and launching capacity, a heavy-orbit interference orbit semi-major axis analytic solution a can be determined J2 ;
Orbit dip analytic solution i J2 :
Track number of turns Q, heavy rail period T:
T=T n ×Q
wherein, T n Is a period of intersection, w e Is the angular velocity of rotation of the earth i J2 Resolving the orbit inclination angle;
module S4.2: according to the J2 perturbation analysis orbit, a J4 perturbation correction orbit based on a Newton iteration method is provided, and the method comprises the following steps:
module s4.2.1: selecting t according to the requirement of satellite task on ground track arrangement 0 Latitude and longitude lambda of time sub-satellite point 0 、Establishing a functional relation of longitude and latitude about a semi-major axis and an inclination angle after the following tracks are repeated;
wherein,
M=n(t-t 0 )
G 0 is t 0 At time Greenwich mean sidereal time omega 0 Is t 0 The right ascension at the ascending crossing point of the moment,is the ascent point right ascension drift rate under J4 perturbation;
module s4.2.2: obtaining the following heavy rail interference track eccentricity e by utilizing the constraint relation of the frozen track J3 And estimating the initial value of the argument w of the near place:
module s4.2.3: first with the satellite t 0 Time ground track initial longitude and latitude lambda 0 、Taking the starting point as the point of integration of the heavy track period T to obtain T 0 Satellite ground track longitude and latitude lambda at + T moment n 、Then calculateObtaining the longitude and latitude deviation delta lambda = lambda of the initial end and the final end of the heavy rail n -λ 0 ,Judging whether the deviation meets the heavy rail interference precision index or not; if the index is satisfied, outputting a semimajor axis and a tilt angle correction value a J4 、i J4 If the iteration parameter is out of tolerance, the following parameter correction is carried out, and the integral judgment process is converted back again until the iteration parameter meets the precision index;
a k =a k-1 +Δa k-1
i k =i k-1 +Δi k-1
wherein,
a k for the semi-major axis of the orbit at time k, a k-1 Is the semimajor axis of the track at time k-1, i k For the track inclination at time k, i k-1 Is the inclination angle of the orbit at the moment of k-1;
module s4.2.4: t is obtained by calculation according to the function relation between the semi-major axis and the inclination angle of the orbit in the module S4.2.1 and the longitude and latitude of the satellite point 0 Time a J4 、i J4 Corresponding rising point right ascension omega J4 Angle M close to the mean J4 。
5. A high accuracy ground track repetitive track optimization system as set forth in claim 4, wherein said repetitive track accuracy is less than 0.01m.
6. The high accuracy ground track repetitive track optimization system of claim 4, wherein the module C comprises:
module S5.1: the method comprises the steps of converting a track repeated nonlinear parameter solving problem under a high-order gravity field into a multivariable and multi-target optimization problem, guiding parameter self-adaptive setting by combining target track characteristic information, and describing an optimization model as follows:
optimizing variables: [ Δ a, Δ e, Δ i, Δ Ω, Δ w, Δ M ]
Initial conditions: [ a ] A J4 ,e J3 ,i J4 ,Ω J4 ,w=90°,M J4 ]
Constraint conditions are as follows:
delta a is correction of track semimajor axis, delta e is correction of track eccentricity, delta i is correction of track inclination angle, delta omega is correction of right ascension at track ascending intersection point, delta w is correction of amplitude at track perigee, delta M is correction of angle at track mean perigee,as a vector of the acceleration of the satellite,in order to be a satellite position deviation,for the satellite position vector at time T,is t 0 A time satellite position vector;
module S5.2: selecting individuals through a binary system tournament method, and performing crossing and variation to generate a new population;
module S5.3: calculating and updating a new population objective function value, namely the position and speed deviation of the satellite earth-fixed system;
module S5.4: generating a new combined population by a merging method, and carrying out non-dominant sorting;
module S5.5: selecting individuals to form a new generation of population through a displacement and elite retention strategy;
module S5.6: and jumping to a module S5.2, and circularly updating until a termination condition is met.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011138645.3A CN112257343B (en) | 2020-10-22 | 2020-10-22 | High-precision ground track repetitive track optimization method and system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011138645.3A CN112257343B (en) | 2020-10-22 | 2020-10-22 | High-precision ground track repetitive track optimization method and system |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112257343A CN112257343A (en) | 2021-01-22 |
CN112257343B true CN112257343B (en) | 2023-03-17 |
Family
ID=74263648
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011138645.3A Active CN112257343B (en) | 2020-10-22 | 2020-10-22 | High-precision ground track repetitive track optimization method and system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112257343B (en) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113190911B (en) * | 2021-03-11 | 2023-05-05 | 上海卫星工程研究所 | Regional multi-target satellite detection simulation method and system |
CN113378290B (en) * | 2021-05-12 | 2022-11-11 | 北京航空航天大学 | Orbit maintaining method for ultra-low orbit satellite |
CN113779765B (en) * | 2021-08-12 | 2023-07-07 | 深圳市魔方卫星科技有限公司 | Heavy orbit satellite orbit optimization method, system, computer equipment and storage medium |
CN113627029B (en) * | 2021-08-24 | 2022-07-12 | 南京大学 | Method for designing track near earth-moon triangle translational point |
CN113734468B (en) * | 2021-08-30 | 2023-02-03 | 北京宇航系统工程研究所 | Orbital plane accurate control method based on iterative guidance |
CN113779788B (en) * | 2021-09-02 | 2023-12-12 | 上海卫星工程研究所 | Method and system for determining combined body separation track under large track-in deviation condition |
CN114383619B (en) * | 2021-12-07 | 2023-09-05 | 上海航天控制技术研究所 | High-precision track calculation method |
CN114440886B (en) * | 2021-12-30 | 2023-09-05 | 上海航天控制技术研究所 | High-accuracy track calculation method for large-eccentricity track |
CN114547549B (en) * | 2022-01-18 | 2024-06-21 | 上海卫星工程研究所 | Atmospheric density autonomous estimation method and system based on low-orbit satellite orbit attenuation |
CN115096319B (en) * | 2022-08-24 | 2022-11-18 | 航天宏图信息技术股份有限公司 | Method and device for determining initial orbit of satellite in star chain based on optical angle measurement data |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110378012A (en) * | 2019-07-16 | 2019-10-25 | 上海交通大学 | A kind of stringent regression orbit design method considering high-order gravitational field |
CN110689505A (en) * | 2019-12-11 | 2020-01-14 | 长沙天仪空间科技研究院有限公司 | Scene-based satellite-borne remote sensing instrument self-adaptive correction method and system |
CN110733671A (en) * | 2019-11-22 | 2020-01-31 | 北京理工大学 | small celestial body spin angular velocity dynamics correction method |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105867119B (en) * | 2016-01-15 | 2018-08-28 | 南京航空航天大学 | A kind of big envelope curve method for handover control of re-entry space vehicle using protection mapping theory |
CN106092105A (en) * | 2016-06-03 | 2016-11-09 | 上海航天控制技术研究所 | A kind of determination method of the strict regression orbit of near-earth satellite |
CN107797130B (en) * | 2017-10-16 | 2021-01-05 | 中国西安卫星测控中心 | Method for calculating uplink data of multi-point and multi-parameter orbit of low-orbit spacecraft |
CN109032176B (en) * | 2018-07-25 | 2021-06-22 | 西北工业大学 | Geosynchronous orbit determination and parameter determination method based on differential algebra |
CN110686684B (en) * | 2019-11-22 | 2021-09-24 | 北京理工大学 | Optical collaborative orbit determination method for small celestial body surrounding detector |
-
2020
- 2020-10-22 CN CN202011138645.3A patent/CN112257343B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110378012A (en) * | 2019-07-16 | 2019-10-25 | 上海交通大学 | A kind of stringent regression orbit design method considering high-order gravitational field |
CN110733671A (en) * | 2019-11-22 | 2020-01-31 | 北京理工大学 | small celestial body spin angular velocity dynamics correction method |
CN110689505A (en) * | 2019-12-11 | 2020-01-14 | 长沙天仪空间科技研究院有限公司 | Scene-based satellite-borne remote sensing instrument self-adaptive correction method and system |
Also Published As
Publication number | Publication date |
---|---|
CN112257343A (en) | 2021-01-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112257343B (en) | High-precision ground track repetitive track optimization method and system | |
CN110378012B (en) | Strict regression orbit design method, system and medium considering high-order gravity field | |
CN106697333B (en) | A kind of robust analysis method of spacecraft orbit control strategy | |
CN109255096B (en) | Geosynchronous satellite orbit uncertain evolution method based on differential algebra | |
CN107797130B (en) | Method for calculating uplink data of multi-point and multi-parameter orbit of low-orbit spacecraft | |
CN109032176B (en) | Geosynchronous orbit determination and parameter determination method based on differential algebra | |
CN101059349A (en) | Minitype combined navigation system and self-adaptive filtering method | |
CN110053788B (en) | Constellation long-term retention control frequency estimation method considering complex perturbation | |
CN103853887B (en) | The satellite orbit of the eccentricity of a kind of Frozen Orbit determines method | |
CN110816896B (en) | Satellite on-satellite simple orbit extrapolation method | |
CN105136166A (en) | Strapdown inertial navigation error model simulation method with specified inertial navigation position precision | |
CN111731513B (en) | Method for maintaining regression orbit in high-precision gravitational field based on monopulse orbit control | |
CN108021138A (en) | A kind of Geomagnetic Field Model simplifies design method | |
CN110059285B (en) | Consider J2Item-influenced missile free-section trajectory deviation analysis and prediction method | |
CN113609708A (en) | Mars global remote sensing track design method and system based on near-fire drift | |
CN111814313B (en) | Regression orbit design method in high-precision gravitational field | |
CN111428912B (en) | Mars detector orbit prediction method and system based on support vector machine | |
CN110231619B (en) | Radar handover time forecasting method and device based on Enk method | |
CN115392540A (en) | Rapid forecasting method for lunar orbit rendezvous guidance | |
CN115265540A (en) | Method and device for acquiring strict regression orbit parameters | |
Allasio et al. | GOCE mission: design phases and in-flight experiences | |
CN112327333A (en) | Satellite position calculation method and device | |
CN113282097A (en) | Control method for non-spherical perturbation error of relative position of GEO game spacecraft | |
Feng et al. | A framework of gravity field online modeling and trajectory optimization in asteroid soft-landing mission scenarios | |
CN111547274A (en) | Spacecraft high-precision autonomous target forecasting method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |