CN111268176A - Perturbation track four-pulse intersection rapid optimization method - Google Patents
Perturbation track four-pulse intersection rapid optimization method Download PDFInfo
- Publication number
- CN111268176A CN111268176A CN202010052133.9A CN202010052133A CN111268176A CN 111268176 A CN111268176 A CN 111268176A CN 202010052133 A CN202010052133 A CN 202010052133A CN 111268176 A CN111268176 A CN 111268176A
- Authority
- CN
- China
- Prior art keywords
- pulse
- spacecraft
- orbit
- time
- intersection
- 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.)
- Granted
Links
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
- B64G1/242—Orbits and trajectories
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Radar, Positioning & Navigation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Navigation (AREA)
Abstract
The invention discloses a perturbation track four-pulse intersection rapid optimization method, which comprises the following steps: according to known conditions of the spacecraft and the target, obtaining control quantity and minimum speed increment through a control model taking elimination of orbital root difference at the intersection time as an optimization result; establishing a four-pulse rendezvous optimization model according to the control quantity and the speed increment, taking values of the time and the size of the first two pulses in the neighborhood according to the estimated value of the step, and recalculating the last two pulses according to the actual orbit recursion result; and replacing the analytic perturbation kinetic model with a full perturbation accurate kinetic model, taking the optimal solution of the analytic kinetic model as an initial value, and solving again to obtain the optimal solution of the four-pulse intersection under the accurate kinetic model. The problem of low computing speed and the like in the prior art is solved, and computing efficiency is improved.
Description
Technical Field
The invention relates to the technical field of aerospace navigation control, in particular to a perturbation orbit four-pulse intersection rapid optimization method.
Background
On the near-earth orbit, the earth non-spherical gravitational perturbation, the atmospheric resistance and the like are the most main perturbation items. When pulse-cross orbit optimization is carried out, the influence of a perturbation term must be considered, but the corresponding orbit integral calculation amount is obviously increased compared with a two-body analysis dynamic model. In addition, the existing method for solving the optimal speed increment of the rail intersection lacks reasonable initial value setting, the solution space range needing to be searched when an evolutionary algorithm is used is large, the calculation time is long, and the convergence is slow particularly when the intersection transfer time is long. Most of the existing algorithms use analytic perturbation models, but the optimal calculation speed is still not ideal and has a certain difference with an accurate dynamic model.
Disclosure of Invention
The invention provides a perturbation orbit four-pulse intersection rapid optimization method which is used for overcoming the defects of long calculation time, slow convergence and the like in the prior art, realizing the accuracy of an iterative model and improving the calculation efficiency.
In order to achieve the purpose, the invention provides a perturbation track four-pulse intersection rapid optimization method, which comprises the following steps:
step 1, according to known conditions of a spacecraft and a target, obtaining control quantity and minimum speed increment through a control model taking elimination of orbital root difference at a meeting time as an optimization result;
changing the number of orbits when the spacecraft is controlled to run in the first circle of orbits according to the pulse which converts a part of control quantity and speed increment into two fixed periods at intervals; converting the residual control quantity and the speed increment into two other pulse control spacecrafts with fixed intervals to finish intersection when the last circle of the orbit runs; estimating the size of four pulses and components in three directions of radial direction, tangential direction and normal direction;
step 2, establishing a four-pulse rendezvous optimization model, taking values of the time and the size of the first two pulses in the neighborhood according to the estimated value in the step 1, and recalculating the last two pulses according to the actual orbit recursion result;
and 3, replacing the analytic perturbation kinetic model with a full perturbation accurate kinetic model, and solving the step 2 again by taking the optimal solution of the analytic kinetic model as an initial value to obtain the optimal solution of the four-pulse intersection under the accurate kinetic model.
The perturbation orbit four-pulse intersection rapid optimization method provided by the invention is suitable for the optimization problem of the orbit intersection with large ascension difference at the ascending intersection point and more similar numbers of other orbits. The method comprises the steps of firstly, providing a reasonable initial value for the four-pulse rendezvous optimization by adopting an estimation method according to the difference value of the number of orbits of a spacecraft and a rendezvous target and fixed transfer time, then designing an iterative model for quickly obtaining an accurate solution, obtaining the optimal solution of the four-pulse rendezvous through local search, and obtaining the accurate optimal solution of the four-pulse rendezvous by using less optimization variables; the method solves the problems of large kinetic recursion calculated amount and slow convergence when the transfer time is long in the traditional multi-pulse rendezvous optimization algorithm.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the structures shown in the drawings without creative efforts.
FIG. 1 is a flow chart of a perturbation orbit four-pulse intersection fast optimization method proposed by the invention.
The implementation, functional features and advantages of the objects of the present invention will be further explained with reference to the accompanying drawings.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
It should be noted that all the directional indicators (such as up, down, left, right, front, and rear … …) in the embodiment of the present invention are only used to explain the relative position relationship between the components, the movement situation, etc. in a specific posture (as shown in the drawing), and if the specific posture is changed, the directional indicator is changed accordingly.
In addition, the descriptions related to "first", "second", etc. in the present invention are only for descriptive purposes and are not to be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include at least one such feature. In the description of the present invention, "a plurality" means at least two, e.g., two, three, etc., unless specifically limited otherwise.
In the present invention, unless otherwise expressly stated or limited, the terms "connected," "secured," and the like are to be construed broadly, and for example, "secured" may be a fixed connection, a removable connection, or an integral part; the connection can be mechanical connection, electrical connection, physical connection or wireless communication connection; they may be directly connected or indirectly connected through intervening media, or they may be connected internally or in any other suitable relationship, unless expressly stated otherwise. The specific meanings of the above terms in the present invention can be understood by those skilled in the art according to specific situations.
In addition, the technical solutions in the embodiments of the present invention may be combined with each other, but it must be based on the realization of those skilled in the art, and when the technical solutions are contradictory or cannot be realized, such a combination of technical solutions should not be considered to exist, and is not within the protection scope of the present invention.
Example one
As shown in fig. 1, an embodiment of the present invention provides a method for rapidly optimizing perturbation orbit four-pulse intersection, which is characterized by comprising the following steps:
step 1, according to known conditions of a spacecraft and a target, obtaining control quantity and minimum speed increment through a control model taking elimination of orbital root difference at a meeting time as an optimization result;
changing the number of orbits when the spacecraft is controlled to run in the first circle of orbits according to the pulse which converts a part of control quantity and speed increment into two fixed periods at intervals; converting the residual control quantity and the speed increment into two other pulse control spacecrafts with fixed intervals to finish intersection when the last circle of the orbit runs; estimating the size of four pulses and components in three directions of radial direction, tangential direction and normal direction;
given that the spacecraft and the target are at the transfer start time t0In the case of the number of tracks of (1), the calculation is at tfTime-of-day intersection, transferring an optimal four-pulse (including the time and magnitude of application of each pulse) of duration Δ t, said step 1 comprising:
step 11, respectively calculating the orbit root difference of the spacecraft and the target which need to be eliminated by the pulse at the intersection moment according to the known conditions and the J2 analytic dynamic model;
step 12, designing an optimal control model for eliminating the track root number difference, and obtaining control quantity and minimum speed increment required by intersection after pulse elimination;
step 13, taking the solution obtained in the step 12 as an initial value, and calculating the minimum speed increment required by intersection under a high-order approximation model;
step 14, correcting the minimum speed increment obtained in the step 13 according to the phase difference to obtain the minimum speed increment required by intersection after the phase difference is corrected;
and step 15, correcting the minimum speed increment obtained in the step 14 according to the eccentricity vector difference, and obtaining the estimated size of each pulse and components in three directions, namely radial direction, tangential direction and normal direction.
Step 2, establishing a four-pulse rendezvous optimization model, taking values of the time and the size of the first two pulses in the neighborhood according to the estimated value in the step 1, and recalculating the last two pulses according to the actual orbit recursion result;
and 3, replacing the analytic perturbation kinetic model with a full perturbation accurate kinetic model, and solving the step 2 again by taking the optimal solution of the analytic kinetic model as an initial value to obtain the optimal solution of the four-pulse intersection under the accurate kinetic model.
The known spacecraft and the target are at a starting time t0The number of tracks. Let tfFor the meeting time, the transition duration from the transition start time to the meeting time is Δ t. The number of the tracks is respectively a semi-major axis a of the track, an eccentricity e of the track, an inclination angle i of the track, a right ascension channel omega of a rising intersection point, an angle omega of a near center point and an angle M of a horizontal near point.
Respectively calculating the intersection time t according to the J2 analytic dynamic modelfThe number of spacecraft orbits and the number of target orbits. Wherein, the J2 analytic kinetic model is as follows:
wherein:the derivative of the number of tracks with respect to time, J2Being a second order term of Earth's gravity, ReIs the radius of the earth, p ═ a (1-e)2) The semi-diameter of the track is the track semi-diameter,mu is the earth's gravitational constant, orbital angular velocity.
The orbit root difference [ Delta a ] of the spacecraft and the target at the moment of intersection, which needs to be eliminated by the pulse0,Δe0,Δi0,ΔΩ0,Δω0,ΔM0]Is tfSubtracting the orbit number of the spacecraft from the target orbit number at the moment.
Elimination of Delta omega in step 12 of the invention0The optimal control model is as follows: applying control to the spacecraft at the transfer starting moment to enable the number of the orbits of the spacecraft to generate variation delta a, delta i and delta omega ', wherein the delta a is the variation of the semi-major axis of the orbits of the spacecraft at the transfer starting moment, the delta i is the variation of the inclination angle of the orbits of the spacecraft at the transfer starting moment, and the delta omega' is the variation of the right ascension channel of the ascending nodes of the spacecraft at the transfer starting moment; applying control quantities with the same reverse direction to the spacecraft at the meeting moment; Δ a, Δ i, Δ Ω' are required to satisfy the condition that the right ascension of the spacecraft at the intersection relative to the target at the intersection is zeroThe beam, under this constraint, solves for the optimal Δ a, Δ i, Δ Ω' distribution in the control quantities, making the speed increment take a minimum value. The specific implementation method comprises the following steps:
for the semi-major axis a and the ascent point right ascension Ω of the spacecraft orbit, according to the J2 analytic dynamic model, that is, equation (1), when the orbit inclination angle of the spacecraft and the rendezvous target is close and the orbit eccentricity is close to 0, it is considered that:
the control quantity delta i and delta a of the orbit inclination angle and the orbit semi-major axis of the spacecraft can cause the right ascension at the lifting point of the spacecraftChange of (2):
the pulse application mode is as follows: at the transfer start time tfApplication of a single pulse Δ v at0Changes the semi-major axis and the orbit inclination angle of the orbit of the spacecraft, and further changes the ascension point right ascension channel of the spacecraftThen at the meeting time tfApplying a single pulse DeltavfRestoring the semi-major axis and the orbit inclination angle of the spacecraft orbit to the initial values when the control quantity is not applied, and then tfThe right ascension difference of the intersection target relative to the spacecraft at the moment is as follows:
the speed increment required for eliminating delta omega' is divided into two and synthesized into double pulses, and then in order to eliminate delta omega0The required speed increments for the spacecraft are:
Δ v in formula (5)a,Δvi,ΔvΩRespectively the speed increment corresponding to the variation delta a of the semi-major axis of the orbit of the spacecraft, the variation delta i of the inclination angle of the orbit and the variation delta omega' of the right ascension of the ascending intersection point,the orbit average velocity is the initial moment of the spacecraft and sini ≈ sin (i + Δ i) needs to be assumed.
The estimation problem of the optimal speed increment is converted into the solution of Δ vsumI.e. Δ v0The extreme value problem of (2). The estimation problem for the optimal speed increment can be described as:
WhereinThe method for calculating the red channel drift rate of the lifting intersection point at the initial time of the spacecraft is obtained by only retaining a first order term according to a J2 analytic dynamic model, namely formula (1), in linear expansion of the number of orbits of the spacecraft at the initial time.
Considering the normalized variables x Δ a/a, y Δ i, z Δ Ω', f and g are equivalent to:
defining L ═ f + λ g, where λ is a constraint multiplier, then the extremum condition is:
equation (9) is a fourth order linear equation that can be directly solved:
only Δ Ω is eliminated0The minimum speed increments required are:
further, the implementation method of step 13 of the present invention is as follows:
to at tfTime of day elimination of Δ a0And Δ i0Δ v in formula (5)fThe rewrite is:
when Δ a and Δ i are large, the calculation is performed by equation (7)Is large, the calculation of formula (3) is adoptedAnd the semi-major axis of the spacecraft orbit can not be less than a in the process of meeting maneuverminConstraint (a)minTaken as 6600km), a more accurate solution model is obtained as:
wherein a isminIs the minimum spacecraft orbit semi-major axis of constraint.
Likewise, let x be Δ a/acY ═ Δ i, z ═ Δ Ω', and L ═ f + λ g1+χg2Wherein λ and χ are constraint multipliers, respectively, and the corresponding extreme conditions are:
whereinCalculated when a is replaced by a + ax and i is replaced by i + y in equation (2)Namely, it is
At the moment, the equation set (15) is nonlinear, x, y, z and lambda solved in the formula (11) in the step 12 are used as initial values of x, y, z and lambda in the equation set (15), the initial value of chi in the equation set (15) is set to be 0, multi-step iterative correction is carried out, new x, y, z, lambda and chi are obtained, and then new x, y, z, lambda and chi are obtained
Further, the implementation method of step 14 of the present invention is as follows:
let Δ λ0=Δω0+ΔM0The spacecraft and the target are at t without controlfA phase difference of time, applying Δ v0After, tfThe time phase difference is changed to:
normalizing Δ λ to [ - π, π]To satisfy the convergence constraint, Δ a and Δ i need to be modified to produce a relative drift rate of the phase
Considering sin2(i + Δ i) is insensitive to variations in i, changes in Δ a are considered dominant; the correction amount is:
superimpose it on Δ v0I.e. fixing x in formula (5)fix=(Δa+Δamod) A, the system of equations translates to:
using the new y, z, λ obtained by iterative correction in step 13 as the initial values of y, z, λ in equation (20), re-solving Δ i and Δ Ω' to obtain the minimum speed increment required for intersection after correcting the phase difference
Further, the implementation method of step 15 of the present invention is as follows:
applying Δ v to a spacecraft at an initial time0Later, the drift rate of the spacecraft ω changes, corresponding to the crossing time tfThe new ω is:
tfthe difference of the orbit eccentricity of the spacecraft and the target at the moment is
According to Δ v0、ΔvfThe normal component of the two pulses to the change quantity of the orbit inclination angle i and the ascension point right ascension omega of the spacecraft, and the control phase of the two pulses is known to be
Wherein at u0To apply Δ v0The variation of the eccentricity of the spacecraft orbit rotates along with the drift of omega, and t is considered to befThe phase for this change in time is:
if Δ v is to be reduced0,ΔvfSplit into two pulses (defined as av) separated by half a period1、Δv2Equivalent to Δ v0,Δv3、Δv4Equivalent to Δ vf) The distribution proportion is different according to the pulse size, and the change quantity delta e of the eccentricity ratio of the spacecraft orbit isv0And Δ evfComprises the following steps:
-Δa/a<|Δev0|<Δa/a
-(Δa-Δa0)/a<|Δevf|<(Δa-Δa0)/a (25)
let k0、kf∈[0,1]Are respectively Δ v0,ΔvfThe delta v is obtained after the superposition of the variation coefficient of the eccentricity of the spacecraft orbit0,ΔvfChange quantity delta e to orbit eccentricity of spacecraftx,Δey]Comprises the following steps:
the existing control quantity is used for eliminating the eccentricity difference relative to the target, and the method is equivalent to solving the extreme value problem:
after the solution is obtained in min (f),if Δ eremA value of 0, indicating that no additional control is required; if Δ eremValues greater than zero require an additional amount of radial control: Δ vr=ΔeremV, Δ V in step (4) according to equation (28)fIn the middle of superposition, the estimated optimal meeting speed increment is obtained
Wherein c is1For the radial pulse component Δ vrAt Δ v0、ΔvfAccording to the distribution ratio of Δ v obtained in step (4)0、ΔvfValue calculationΔv1,Δv2,Δv3,Δv4The sizes are respectively as follows:
wherein k is0,k0See equation (27), where the semimajor axis component corresponds to the tangential component of the pulse, the dip and lift-off right ascension components correspond to the normal component, and the remainder corresponds to the radial component. By the formula (23), Δ v1,Δv2,Δv3,Δv4The phase angles of (a) are respectively:
wherein u is1And u2First turn after task start, u3And u4At the last turn.
Preferably, the step 2 optimizes the four-pulse intersection optimization model by using a differential evolution algorithm.
Preferably, the step 2 includes:
step 21, directly applying the four-pulse estimation value obtained in the step 1 to a spacecraft, recurrently delivering the orbit to a known rendezvous moment, and defining a local search range coefficient of a speed increment and a local search range coefficient of a control quantity according to the deviation between the recurrently delivering orbit and a target orbit; the step 21 includes:
in step 211, the four-pulse estimation values obtained in step 1 are respectively: Δ v1,Δv2,Δv3,Δv4The phase angles are respectively: u. of1,u2,u3,u4Wherein u is1,u2For the first turn after the start of the rendezvous task, u3,u4The last circle after the rendezvous task starts;
at step 212, four pulses are applied directly to the spacecraft and the orbit is recurred until time tfDeviation from the target trajectory, defining 4 coefficients [ x ]1x2x3x4]In a local search, where x1,x2,x3Local search range coefficients, x, representing the tangential, normal and radial components of the velocity increment, respectively4A local search range coefficient representing the corresponding eccentricity control quantity in the step 1; equations (29) and (30), i.e. the estimated values of the four pulses, are applied directly to the spacecraft and the orbit is recurred to time tfThere will be some small deviation from the target trajectory. The invention considers that the difference between the estimated value and the actual optimal solution is small, and only local search is needed.
Step 213, convert Δ v1、Δv2The changes are as follows:
wherein, delta a is the variation of the semi-major axis of the orbit of the spacecraft at the moment of the transfer initiation, delta i is the variation of the inclination angle of the orbit of the spacecraft at the moment of the transfer initiation, delta omega is the variation of the right ascension channel of the ascending intersection point of the spacecraft at the moment of the transfer initiation, and delta vrTo require an additional amount of radial control, a1For a predetermined upper limit of the local change of the semi-major axis, k0For controlling space flightThe coefficient of the variation of the eccentricity of the orbit, i is the inclination angle of the orbit,the average orbit velocity of the spacecraft.
Step 22, determining the moments of four pulses through an optimization model established by the time of pulse intersection, the phase of the spacecraft and the orbital angular velocity when the spacecraft intersects; preferably, said step 22 comprises:
step 221, defining an optimization model, first, the number of orbits of the spacecraft is deduced to the starting moment t of the handover controlf- Δ t, set to [ as,es,is,Ωs,ωs,Ms]If the phase at that time is u ═ ωs+MsDue to the pulse Δ v1Is u1Then the pulse Deltav is known1Time deltat relative to the start of a session1Comprises the following steps:
Δt1=mod(u1-u,2π)/n(32)
wherein mod represents the sum of1-u is normalized to between 0 and 2 π,for the orbital angular velocity at the beginning of a spacecraft rendezvous, namely, the spacecraft rendezvous needs to wait for delta t1After that, the first pulse is applied, and Δ v is known from the formula (31)1The three-dimensional components in the spacecraft orbit coordinate system are:
the number of orbits of the spacecraft is pushed to tf-Δt+Δt1At the moment, the number of tracks is converted into a position velocity vector, and delta v is superposed on the velocity1Obtaining a spacecraft position and speed vector after the first pulse, and then converting to obtain the orbit number;
step 222, due to Δ v1And Δ v2Phase difference of u2-u1Pi, so that the second pulse is separated from the first by a time at2Comprises the following steps:
Δt2=π/n1(34)
whereinTo add Δ v1New orbital angular velocity of the spacecraft after, i.e. passing at2Applying a second pulse of formula (31), Δ v2The three-dimensional components in the spacecraft orbit coordinate system are:
then the orbit number of the spacecraft after the second pulse is pushed to tfTime is set as [ a ]f,ef,if,Ωf,ωf,Mf]The number of tracks [ a ] at the same time as the targett(tf),et(tf),it(tf),Ωt(tf),ωt(tf),Mt(tf)]The difference shows that the right ascension angle Δ Ω at this timefAnd the difference of inclination angle Δ if;
Step 223, according to the principle of equation (30), the phases u of the third and the fourth pulses are determined3,u4The change is:
u3=atan((Δif)/(ΔΩfsini/2))
u4=atan((Δif)/(ΔΩfsini/2))+π (36)
the time at which the fourth pulse is applied and tfThe moments not necessarily coinciding, with a time difference Δ t4Comprises the following steps:
Δt4=-(mod(u4-(ωf+Mf),2π)-2π)/nf(37)
whereinTrack angular velocity for the rendezvous target, (mod (u)4-(ωf+Mf),2π)-2π) Indicating handle u4-(ωf+Mf) Normalized to between-2 pi and 0, i.e. the fourth pulse has a time tf-Δt4The time advance delta t of the third pulse relative to the fourth pulse3Comprises the following steps:
Δt3=π/n3(38)
wherein n is3For the angular velocity of the spacecraft orbit after the third pulse is applied, since n3Is difficult to be directly calculated, and introduces a correction coefficient x to be solved for compensating eccentricity5Order:
Δt3=x5π/nf(39)
the timings of the third and fourth pulses are determined by equations (37) and (39). But its size has not yet been determined. Calculating Δ ν using Lambert's algorithm taking perturbation into account3And Δ v4The size of (2). The Lambert problem can be described as a given start-stop time instant (t)inAnd tout) And a corresponding start-stop position vector (r)inAnd rout) Calculating to obtain a start-stop velocity vector (v)inAnd vout) The position and speed state of the spacecraft is represented as rinvin]Time passes tout-tinThe long time orbit recursion can reach the state routvout]. The method has an analytic calculation method for a two-body dynamics model, has more descriptions in open documents, and is regarded as a function in the invention, which is shown in a formula (40).
The magnitude of the last two pulses is calculated by the Lambert algorithm taking perturbation into account, step 23. Preferably, the step 23 includes:
step 231, for a given start-stop time tin、toutAnd corresponding start-stop position vector rin、routCalculating a start-stop velocity vector v for a two-body dynamics modelin、vout:
[vinvout]=Lambert(rin,rout,tin,tout) (40)
Applying the Lambert algorithm, i.e. setting the starting time tinIs set to tf-Δt4-Δt3Terminate atTime toutIs set to tf-Δt4At a starting position rinFor the number of orbits of the spacecraft after the second pulse to be pushed to tf-Δt4-Δt3Position vector obtained at time, end position routStepping to t for target trajectoryf-Δt4A position vector obtained at a time;
step 232, firstly, obtaining the v corresponding to the two-body dynamic model by using an analytic calculation methodinAnd voutAnd obtaining v corresponding to the J2 analytic perturbation dynamics model by differential correctioninAnd vout;
Firstly, the v corresponding to the two-body dynamic model is obtained by using a public analytical calculation methodinAnd voutThen, the differential correction is used to obtain v corresponding to the perturbation dynamics model formula (1)inAnd vout;
Step 233, set v3 -For the number of orbits of the spacecraft after the second pulse to be pushed to tf-Δt4-Δt3The velocity vector taken at the moment, the third pulse needs to change the spacecraft velocity to vinThe meeting location constraint can be satisfied; v. the4 +Stepping to t for target trajectoryf-Δt4The velocity vector obtained at the moment, the fourth pulse needs to drive the spacecraft velocity from voutChange to v4 +The meeting speed constraint can be met; thus Δ v3And Δ v4The sizes are respectively as follows:
the total velocity increment for the four pulses is then:
Δvsum=Δv1+Δv2+Δv3+Δv4(42)
to sum up, the local optimization model includes 5 variables: x is the number of1,x2,x3,x4,x5The optimization index is Δ v calculated by the following equations (31) to (42)sumMinimum; x is the number of1,x2,x3,x4,x5The value ranges are as follows:
solving by using a differential evolution algorithm to obtain optimal x1,x2,x3,x4,x5And accurate values of the four pulse moments and the four pulse magnitudes under the analytic dynamics model can be obtained.
The differential evolution algorithm can be easily solved, and is a general algorithm for solving the optimization problem, and is not described herein again. Obtaining the optimal x1,x2,x3,x4,x5And then, accurate values of the time and the size of the four pulses under the analytic kinetic model can be obtained.
Preferably, the step 3 comprises:
step 31, replacing the J2 analytic perturbation kinetic model of formula (1) with the following accurate kinetic equation, and applying to all processes involving orbit recursion calculation in step 2:
wherein r is the three-dimensional position vector of the spacecraft, v is the three-dimensional velocity vector, aPerturbationThe acceleration caused by the gravity forces such as the non-spherical gravity of the earth, the atmospheric resistance, the sun-moon gravity, the sunlight pressure and the like is considered; the calculation formula can be obtained in public teaching materials, and is not described in detail herein. It is obvious that the optimal solution of step (6) cannot be directly used as the corresponding optimal solution of the exact kinetic equation (44).
Step 32, using x solved in step 22,x3,x4Re-optimizing x as an optimal solution for accurate kinetic models1,x5;
And step 33, recalculating the optimization models corresponding to the expressions (31) to (42) by using a differential evolution algorithm, so as to obtain the four pulse moments and the four pulse magnitudes with the optimal precise dynamics.
The main difference between the analytic kinetic equation and the accurate kinetic equation is that accumulated phase difference is generated after long-time recursive calculation and can be eliminated through semi-long axis correction. Therefore, in order to further improve the calculation efficiency, x obtained by the solution of the step (6) is used2,x3,x4Re-optimizing only x as an optimal solution for accurate kinetic models1,x5And (4) finishing. The number of the optimized variables is reduced to 2, and the solving speed is higher. Recalculating the optimization models (fixed x) corresponding to equations (31) to (42) by applying a differential evolution algorithm2,x3,x4) Then, four pulse time moments and sizes optimal for precise dynamics can be obtained.
Next, the method of the present invention is simulated to obtain the optimal speed increment through a specific embodiment.
In this embodiment: spacecraft and target initial time t0The number of orbitals is 0, as shown in table 1:
TABLE 1 spacecraft and target initial time orbital radical
A(km) | e | i (degree) | Omega (degree) | w (degree) | M (degree) | |
Spacecraft | 6638.14 | 0.009039 | 42.05 | 172.6 | 120 | 1 |
Target | 6720.14 | 0.00001 | 42 | 169.2 | 100 | 145 |
To calculate the optimal four-pulse control required for a 14 th day rendezvous starting at the initial time compared to the track transfer. I.e. Δ t equals 14 days, tfDay 14.
Firstly, calculating the orbit root difference of the spacecraft and the target at the meeting time.
First extrapolate the number of tracks to tfThe track root difference [ Delta a ] is obtained at the moment0,Δi0,ΔΩ0]Is [82000 m0.0008730.00585 ]]. Δ i here0,ΔΩ0All are radians.
And secondly, calculating the control quantity required by the intersection according to a first-order approximate model.
According to Δ Ω0Solving the sum formula (11) to obtain [ x y z]=[0.00103 6.61e-5 1.04e-4]. Then the optimal control quantity [ delta a, delta i, delta omega ] at the initial moment required for eliminating the ascension crossing point declination under the first-order approximation model']Are respectively [6814.99m 6.61 e-51.04 e-4 ]]。
And thirdly, calculating the control quantity required by intersection under the high-order approximation model by taking the solution of the second step as an initial value.
Using it as an initial value, substituting it into formula (15), wherein aminThe radius of the earth plus 200km is set as the solution to [ x y z]=[0.00104 -7.32e-5 0]. Substitution of formula (14) to eliminate [ Delta a ]0,Δi0,ΔΩ0]Required optimum control amount Δ v0And Δ vfRespectively 4.06m/s and 44.27 m/s.
And step four, correcting the control quantity obtained in the step three according to the phase difference.
According to the number of orbits of the spacecraft and delta v0Can be recalculated at tfThe phase difference of the time target relative to the spacecraft is-3.31 rad, and the semi-long axis correction quantity required by corresponding delta t is 10381.17 m. Substitute it into xfix=(Δa+Δamod) (ii) when/a is 0.0026, the formula (20) is solved to obtain [ y z%]=[-0.0014 -0.0067]. Re-substitution (14), i.e. correcting phase difference to obtain optimum control quantity delta v0And Δ vf22.75m/s and 41.79m/s, respectively.
And fifthly, correcting the control quantity obtained in the fourth step according to the eccentricity ratio vector difference.
According to the number of orbits of the spacecraft and delta v0Calculated as tfThe eccentricity difference of the target at the moment relative to the spacecraft is as follows:u is represented by the formula (23)0And ufRespectively-0.26 and-1.36, and solving the formula (27) to obtain Delta erem0.0052, i.e. Δ vr=ΔeremV was 40.38 m/s. k is a radical of0,kfAre respectively 1 and-0.64, c1Is 0.35. The four pulse sizes thus estimated are: 26.70m/s, 0.13m/s, 8.87m/s, 40.42 m/s.
And sixthly, establishing a four-pulse intersection optimization model and applying a differential evolution algorithm for optimization.
According to the estimated pulse and the local optimization model in the fifth step, solving X ═ X1x2x3x4x5]Is [ -0.21,0.8,1.0, -0.28,0.56 [)]The optimal four pulse sizes when the analytic kinetic equation is used are respectively as follows: 23.33m/s, 9.02m/s, 15.72m/s, 26.64m/s, the sum being 65.68 m/s. The time is respectively as follows: 1760.71 seconds, 4461.18 seconds, 1206291.32 seconds, 1207800.56 seconds. Local optimization takes only 1 second or less.
And seventhly, replacing the analytic perturbation dynamic model with a full perturbation accurate dynamic model (10-order gravitational field, atmospheric resistance (coefficient 2.2, face value ratio 0.003) and daily and monthly gravitational perturbation), and solving the sixth step again by taking the optimal solution of the analytic dynamic model as an initial value.
The optimal four pulse sizes of the accurate dynamic model are respectively as follows: 20.89m/s, 8.07m/s, 9.98m/s, 38.70m/s, the sum being 69.57 m/s. The time is respectively as follows: 1760.71 seconds, 4461.18 seconds, 1205335.53 seconds, 1207465.19 seconds. Local optimization takes less than 2 minutes.
The method has the advantages of small error of the estimated speed increment and high solving speed. The method solves the problem of fast solving of the four-pulse optimal intersection.
The above description is only a preferred embodiment of the present invention, and is not intended to limit the scope of the present invention, and all modifications and equivalents of the present invention, which are made by the contents of the present specification and the accompanying drawings, or directly/indirectly applied to other related technical fields, are included in the scope of the present invention.
Claims (9)
1. A perturbation orbit four-pulse intersection rapid optimization method is characterized by comprising the following steps:
step 1, according to known conditions of a spacecraft and a target, obtaining control quantity and minimum speed increment through a control model taking elimination of orbital root difference at a meeting time as an optimization result;
changing the number of orbits when the spacecraft is controlled to run in the first circle of orbits according to the pulse control method of converting a part of control quantity and speed increment into two fixed intervals; converting the residual control quantity and the speed increment into two other pulse control spacecrafts with fixed intervals to finish intersection when the last circle of the orbit runs; estimating the size of four pulses and components in three directions of radial direction, tangential direction and normal direction;
step 2, establishing a four-pulse rendezvous optimization model, taking values of the time and the size of the first two pulses in the neighborhood according to the estimated value in the step 1, and recalculating the last two pulses according to the actual orbit recursion result;
and 3, replacing the analytic perturbation kinetic model with a full perturbation accurate kinetic model, and solving the step 2 again by taking the optimal solution of the analytic kinetic model as an initial value to obtain the optimal solution of the four-pulse intersection under the accurate kinetic model.
2. The perturbation orbit four-pulse intersection rapid optimization method according to claim 1, wherein the step 2 adopts a differential evolution algorithm to optimize a four-pulse intersection optimization model.
3. A perturbation orbit four-pulse intersection fast optimization method as claimed in claim 2, wherein said step 2 comprises:
step 21, directly applying the four-pulse estimation value obtained in the step 1 to a spacecraft, recurrently delivering the orbit to a known rendezvous moment, and defining a local search range coefficient of a speed increment and a local search range coefficient of a control quantity according to the deviation between the recurrently delivering orbit and a target orbit;
step 22, determining the moments of four pulses through an optimization model established by the time of pulse intersection, the phase of the spacecraft and the orbital angular velocity when the spacecraft intersects;
the magnitude of the last two pulses is calculated by the Lambert algorithm taking perturbation into account, step 23.
4. A perturbation orbit four-pulse intersection fast optimization method as claimed in claim 3, wherein said step 21 comprises:
in step 211, the four-pulse estimation values obtained in step 1 are respectively: Δ v1,Δv2,Δv3,Δv4The phase angles are respectively: u. of1,u2,u3,u4Wherein u is1,u2For the first turn after the start of the rendezvous task, u3,u4The last circle after the rendezvous task starts;
at step 212, four pulses are applied directly to the spacecraft and the orbit is recurred until time tfDeviation from the target trajectory, defining four coefficients [ x ]1x2x3x4]Searching locallyWherein x is1,x2,x3Local search range coefficients, x, representing the tangential, normal and radial components of the velocity increment, respectively4A local search range coefficient representing the corresponding eccentricity control quantity in the step 1;
step 213, convert Δ v1、Δv2The changes are as follows:
wherein, delta a is the variation of the semi-major axis of the orbit of the spacecraft at the moment of the transfer initiation, delta i is the variation of the inclination angle of the orbit of the spacecraft at the moment of the transfer initiation, delta omega is the variation of the right ascension channel of the ascending intersection point of the spacecraft at the moment of the transfer initiation, and delta vrTo require an additional amount of radial control, a1For a predetermined upper limit of the local change of the semi-major axis, k0I is the orbit inclination angle for the coefficient of the variation of the control quantity to the orbit eccentricity of the spacecraft,the average orbit velocity of the spacecraft.
5. A perturbation orbit four-pulse intersection fast optimization method as claimed in claim 4, wherein said step 22 comprises:
step 221, the number of orbits of the spacecraft is pushed to the starting moment t of the handover controlf- Δ t, set to [ as,es,is,Ωs,ωs,Ms]If the phase at that time is u ═ ωs+MsDue to the pulse Δ v1Is u1Then the pulse Deltav is known1Time deltat relative to the start of a session1Comprises the following steps:
Δt1=mod(u1-u,2π)/n (32)
wherein mod represents the sum of1-u is normalized to between 0 and 2 π,for the orbital angular velocity at the beginning of a spacecraft rendezvous, namely, the spacecraft rendezvous needs to wait for delta t1After that, the first pulse is applied, and Δ v is known from the formula (31)1The three-dimensional components in the spacecraft orbit coordinate system are:
the number of orbits of the spacecraft is pushed to tf-Δt+Δt1At the moment, the number of tracks is converted into a position velocity vector, and delta v is superposed on the velocity1Obtaining a spacecraft position and speed vector after the first pulse, and then converting to obtain the orbit number;
step 222, due to Δ v1And Δ v2Phase difference of u2-u1Pi, so that the second pulse is separated from the first by a time at2Comprises the following steps:
Δt2=π/n1(34)
whereinTo add Δ v1New orbital angular velocity of the spacecraft after, i.e. passing at2Applying a second pulse of formula (31), Δ v2The three-dimensional components in the spacecraft orbit coordinate system are:
then the orbit number of the spacecraft after the second pulse is pushed to tfTime is set as [ a ]f,ef,if,Ωf,ωf,Mf]The number of tracks [ a ] at the same time as the targett(tf),et(tf),it(tf),Ωt(tf),ωt(tf),Mt(tf)]The difference shows that the right ascension angle Δ Ω at this timefAnd the difference of inclination angle Δ if;
Step 223, the phases u of the third and fourth pulses3,u4The change is:
u3=atan((Δif)/(ΔΩfsini/2))
u4=atan((Δif)/(ΔΩfsini/2))+π (36)
the time at which the fourth pulse is applied and tfThe moments not necessarily coinciding, with a time difference Δ t4Comprises the following steps:
Δt4=-(mod(u4-(ωf+Mf),2π)-2π)/nf(37)
whereinTrack angular velocity for the rendezvous target, (mod (u)4-(ωf+Mf) 2 π) -2 π) represents u4-(ωf+Mf) Normalized to between-2 pi and 0, i.e. the fourth pulse has a time tf-Δt4The time advance delta t of the third pulse relative to the fourth pulse3Comprises the following steps:
Δt3=π/n3(38)
wherein n is3For the angular velocity of the spacecraft orbit after the third pulse is applied, since n3Is difficult to be directly calculated, and introduces a correction coefficient x to be solved for compensating eccentricity5Let us order
Δt3=x5π/nf(39)
Thereby determining the time instants of the third and fourth pulses.
6. A perturbation orbit four-pulse intersection fast optimization method as claimed in claim 5, wherein said step 23 comprises:
step 231, for a given start-stop time tin、toutAnd corresponding start-stop position vector rin、routCalculating a start-stop velocity vector v for a two-body dynamics modelin、vout:
[vinvout]=Lambert(rin,rout,tin,tout) (40)
Applying the Lambert algorithm, i.e. setting the starting time tinIs set to tf-Δt4-Δt3End time toutIs set to tf-Δt4At a starting position rinFor the number of orbits of the spacecraft after the second pulse to be pushed to tf-Δt4-Δt3Position vector obtained at time, end position routStepping to t for target trajectoryf-Δt4A position vector obtained at a time;
step 232, firstly, obtaining the v corresponding to the two-body dynamic model by using an analytic calculation methodinAnd voutAnd obtaining v corresponding to the J2 analytic perturbation dynamics model by differential correctioninAnd vout;
Step 233, set v3 -For the number of orbits of the spacecraft after the second pulse to be pushed to tf-Δt4-Δt3The velocity vector taken at the moment, the third pulse needs to change the spacecraft velocity to vinThe meeting location constraint can be satisfied; v. the4 +Stepping to t for target trajectoryf-Δt4The velocity vector obtained at the moment, the fourth pulse needs to drive the spacecraft velocity from voutChange to v4 +The meeting speed constraint can be met; thus Δ v3And Δ v4The sizes are respectively as follows:
the total velocity increment for the four pulses is then:
Δvsum=Δv1+Δv2+Δv3+Δv4(42)
to sum up, the local optimization model includes 5 variables: x is the number of1,x2,x3,x4,x5The optimization index is Δ v calculated by the following equations (31) to (42)sumMinimum; x is the number of1,x2,x3,x4,x5The value ranges are as follows:
solving by using a differential evolution algorithm to obtain optimal x1,x2,x3,x4,x5And accurate values of the four pulse moments and the four pulse magnitudes under the analytic dynamics model can be obtained.
7. A perturbation orbit four-pulse intersection fast optimization method as claimed in claim 6, wherein said step 3 comprises:
step 31, replacing the J2 analytic perturbation dynamics model with an accurate dynamics equation of the following formula, and applying the equation to all the processes related to the orbit recursion calculation in step 2:
wherein r is the three-dimensional position vector of the spacecraft, v is the three-dimensional velocity vector, aPerturbationThe acceleration caused by the gravity forces such as the non-spherical gravity of the earth, the atmospheric resistance, the sun-moon gravity, the sunlight pressure and the like is considered;
step 32, using x solved in step 22,x3,x4Re-optimizing x as an optimal solution for accurate kinetic models1,x5;
And step 33, recalculating the optimization models corresponding to the expressions (31) to (42) by using a differential evolution algorithm, so as to obtain the four pulse moments and the four pulse magnitudes with the optimal precise dynamics.
8. A perturbation orbit four-pulse intersection fast optimization method as claimed in any one of claims 1 to 7, wherein in step 1, the known conditions include: initial time of the spacecraft and the target, the number of orbits of the spacecraft and the target at the initial time and rendezvous time;
the predetermined period of the interval between the two pulses is half a period.
9. A perturbation orbit four-pulse intersection fast optimization method as claimed in claim 8, wherein said step 1 comprises:
step 11, respectively calculating the orbit root difference of the spacecraft and the target which need to be eliminated by the pulse at the intersection moment according to the known conditions and the J2 analytic dynamic model;
step 12, designing an optimal control model for eliminating the track root number difference, and obtaining control quantity and minimum speed increment required by intersection after pulse elimination;
step 13, taking the solution obtained in the step 12 as an initial value, and calculating the minimum speed increment required by intersection under a high-order approximation model;
step 14, correcting the minimum speed increment obtained in the step 13 according to the phase difference to obtain the minimum speed increment required by intersection after the phase difference is corrected;
and step 15, correcting the minimum speed increment obtained in the step 14 according to the eccentricity vector difference, and obtaining the estimated size of each pulse and components in three directions, namely radial direction, tangential direction and normal direction.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010052133.9A CN111268176B (en) | 2020-01-17 | 2020-01-17 | Perturbation track four-pulse intersection rapid optimization method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010052133.9A CN111268176B (en) | 2020-01-17 | 2020-01-17 | Perturbation track four-pulse intersection rapid optimization method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111268176A true CN111268176A (en) | 2020-06-12 |
CN111268176B CN111268176B (en) | 2021-06-11 |
Family
ID=70997354
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010052133.9A Active CN111268176B (en) | 2020-01-17 | 2020-01-17 | Perturbation track four-pulse intersection rapid optimization method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111268176B (en) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111731513A (en) * | 2020-06-15 | 2020-10-02 | 航天东方红卫星有限公司 | Method for maintaining regression orbit in high-precision gravitational field based on monopulse orbit control |
CN112084581A (en) * | 2020-09-24 | 2020-12-15 | 中国人民解放军国防科技大学 | Spacecraft low-thrust perturbation intersection trajectory optimization method and system |
CN112464429A (en) * | 2020-10-09 | 2021-03-09 | 清华大学 | Maximum value estimation method for long-term evolution of orbit number of low-thrust spacecraft |
CN113060306A (en) * | 2021-03-31 | 2021-07-02 | 中国空气动力研究与发展中心设备设计与测试技术研究所 | Multi-pulse intersection iterative guidance method and device for limited thrust and electronic equipment |
CN113148232A (en) * | 2021-02-24 | 2021-07-23 | 中国人民解放军国防科技大学 | Maneuvering fixed monopulse rail aiming method and device |
CN113343442A (en) * | 2021-05-21 | 2021-09-03 | 中国科学院国家空间科学中心 | Method and system for solving fixed-time finite fuel multi-pulse transfer orbit |
CN113602535A (en) * | 2021-08-02 | 2021-11-05 | 北京航空航天大学 | Method for controlling in-orbit autonomous intersection of micro/nano satellite and computer equipment |
CN114063645A (en) * | 2021-11-12 | 2022-02-18 | 中国西安卫星测控中心 | Eccentricity dip angle vector-based method for evaluating tilt fly-around holding control effect |
CN116701823A (en) * | 2023-08-07 | 2023-09-05 | 长沙翔宇信息科技有限公司 | Intersection space range estimation method, device, terminal equipment and storage medium |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107526368A (en) * | 2017-09-12 | 2017-12-29 | 北京理工大学 | A kind of multiple-pulse ring moon satellites formation initial method for considering error |
CN109491406A (en) * | 2019-01-22 | 2019-03-19 | 中国人民解放军战略支援部队航天工程大学 | A kind of spacecraft evasion manoeuvre method based on energy consumption |
CN110032768A (en) * | 2019-03-15 | 2019-07-19 | 中国西安卫星测控中心 | A kind of four pulse-orbits intersection optimization method using precise kinetic model |
-
2020
- 2020-01-17 CN CN202010052133.9A patent/CN111268176B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107526368A (en) * | 2017-09-12 | 2017-12-29 | 北京理工大学 | A kind of multiple-pulse ring moon satellites formation initial method for considering error |
CN109491406A (en) * | 2019-01-22 | 2019-03-19 | 中国人民解放军战略支援部队航天工程大学 | A kind of spacecraft evasion manoeuvre method based on energy consumption |
CN110032768A (en) * | 2019-03-15 | 2019-07-19 | 中国西安卫星测控中心 | A kind of four pulse-orbits intersection optimization method using precise kinetic model |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111731513A (en) * | 2020-06-15 | 2020-10-02 | 航天东方红卫星有限公司 | Method for maintaining regression orbit in high-precision gravitational field based on monopulse orbit control |
CN112084581A (en) * | 2020-09-24 | 2020-12-15 | 中国人民解放军国防科技大学 | Spacecraft low-thrust perturbation intersection trajectory optimization method and system |
CN112464429A (en) * | 2020-10-09 | 2021-03-09 | 清华大学 | Maximum value estimation method for long-term evolution of orbit number of low-thrust spacecraft |
CN113148232A (en) * | 2021-02-24 | 2021-07-23 | 中国人民解放军国防科技大学 | Maneuvering fixed monopulse rail aiming method and device |
CN113060306A (en) * | 2021-03-31 | 2021-07-02 | 中国空气动力研究与发展中心设备设计与测试技术研究所 | Multi-pulse intersection iterative guidance method and device for limited thrust and electronic equipment |
CN113060306B (en) * | 2021-03-31 | 2022-02-08 | 中国空气动力研究与发展中心设备设计与测试技术研究所 | Multi-pulse intersection iterative guidance method and device for limited thrust and electronic equipment |
CN113343442A (en) * | 2021-05-21 | 2021-09-03 | 中国科学院国家空间科学中心 | Method and system for solving fixed-time finite fuel multi-pulse transfer orbit |
CN113602535A (en) * | 2021-08-02 | 2021-11-05 | 北京航空航天大学 | Method for controlling in-orbit autonomous intersection of micro/nano satellite and computer equipment |
CN113602535B (en) * | 2021-08-02 | 2024-03-12 | 北京航空航天大学 | Method for controlling micro-nano satellite in-orbit autonomous intersection and computer equipment |
CN114063645A (en) * | 2021-11-12 | 2022-02-18 | 中国西安卫星测控中心 | Eccentricity dip angle vector-based method for evaluating tilt fly-around holding control effect |
CN116701823A (en) * | 2023-08-07 | 2023-09-05 | 长沙翔宇信息科技有限公司 | Intersection space range estimation method, device, terminal equipment and storage medium |
CN116701823B (en) * | 2023-08-07 | 2023-10-27 | 长沙翔宇信息科技有限公司 | Intersection space range estimation method, device, terminal equipment and storage medium |
Also Published As
Publication number | Publication date |
---|---|
CN111268176B (en) | 2021-06-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111268176B (en) | Perturbation track four-pulse intersection rapid optimization method | |
CN110789739B (en) | Method for quickly estimating optimal speed increment of long-time rail crossing under J2 perturbation | |
CN108710303B (en) | Spacecraft relative attitude control method containing multi-source disturbance and actuator saturation | |
CN110032768B (en) | Four-pulse orbit intersection optimization method using accurate dynamic model | |
CN115649491B (en) | Low orbit optical remote sensing satellite staring imaging control method suitable for multi-source interference | |
CN109080854B (en) | Orbit-changing planning method for large elliptical orbit of spacecraft returning to preset drop point | |
CN109649692B (en) | Method and system for on-orbit correction of orbit control strategy based on thruster deviation estimation | |
CN112395689B (en) | Rocket fault post-online reconstruction method based on convex optimization | |
CN106354147A (en) | Control system and control method for posture rail coupling under shifting state of big barycenter based on engine pre-swinging | |
CN109606739B (en) | Method and device for correcting earth-moon transfer orbit of detector | |
CN113900448A (en) | Aircraft prediction correction composite guidance method based on sliding mode disturbance observer | |
CN113343561B (en) | Method and system for solving optimal moon fly-by transfer orbit | |
CN108516107A (en) | A kind of online motor power and specific impulse method of estimation and system | |
CN114889849A (en) | Estimation method for multi-constraint flying-over optimal speed of spacecraft | |
CN113581494A (en) | GEO (geostationary orbit) satellite electric hybrid propulsion orbital transfer method | |
CN111605736A (en) | Earth-moon L2 point transfer orbit optimal error correction point selection method | |
CN110209190B (en) | Satellite nominal orbit unbiased flight control method | |
CN110955255A (en) | High-precision orbit control attitude maintaining method, system and medium based on CMG | |
CN114967729B (en) | Multi-rotor unmanned aerial vehicle height control method and system | |
CN112506209B (en) | Reentry vehicle prediction control method based on self-adaptive prediction period | |
CN112455725A (en) | Method for transferring and converting pulse orbit transfer direction to limited thrust orbit | |
CN109094818B (en) | Remote rendezvous guidance method for spacecraft | |
CN113701755A (en) | Optical remote sensing satellite attitude determination method without high-precision gyroscope | |
CN116817932B (en) | Integrated method for track maintenance, track determination and real-time mapping of gravitational field | |
CN114996841B (en) | Optimal coplanarity transfer search method based on tangent initial value |
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 |