CN105005651A - Gradient separate zone optimization design method for spacecraft pulse rendezvous trajectory - Google Patents
Gradient separate zone optimization design method for spacecraft pulse rendezvous trajectory Download PDFInfo
- Publication number
- CN105005651A CN105005651A CN201510400314.5A CN201510400314A CN105005651A CN 105005651 A CN105005651 A CN 105005651A CN 201510400314 A CN201510400314 A CN 201510400314A CN 105005651 A CN105005651 A CN 105005651A
- Authority
- CN
- China
- Prior art keywords
- interval
- queue
- decision variable
- sub
- lnew
- 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
- 238000000034 method Methods 0.000 title claims abstract description 86
- 238000005457 optimization Methods 0.000 title claims abstract description 62
- 238000013461 design Methods 0.000 title claims abstract description 60
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 34
- 238000012545 processing Methods 0.000 claims abstract description 8
- 238000000926 separation method Methods 0.000 claims abstract description 5
- 230000000694 effects Effects 0.000 claims description 18
- 230000009471 action Effects 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 14
- 238000004458 analytical method Methods 0.000 claims description 13
- 238000012360 testing method Methods 0.000 claims description 9
- 238000012217 deletion Methods 0.000 claims description 7
- 230000037430 deletion Effects 0.000 claims description 7
- 239000000446 fuel Substances 0.000 claims description 6
- 238000007689 inspection Methods 0.000 claims description 6
- 230000001419 dependent effect Effects 0.000 claims description 3
- 230000000541 pulsatile effect Effects 0.000 claims description 3
- 238000004513 sizing Methods 0.000 claims description 3
- 230000011218 segmentation Effects 0.000 abstract description 5
- 238000005516 engineering process Methods 0.000 abstract description 4
- 238000004364 calculation method Methods 0.000 abstract description 2
- 238000013178 mathematical model Methods 0.000 abstract 1
- 230000006870 function Effects 0.000 description 68
- 238000012938 design process Methods 0.000 description 6
- 230000002068 genetic effect Effects 0.000 description 6
- 230000009286 beneficial effect Effects 0.000 description 5
- 238000011160 research Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 239000002245 particle Substances 0.000 description 3
- 230000002159 abnormal effect Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 239000011664 nicotinic acid Substances 0.000 description 2
- 230000009467 reduction Effects 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000005574 cross-species transmission Effects 0.000 description 1
- 238000010141 design making Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000003032 molecular docking Methods 0.000 description 1
- 239000003380 propellant Substances 0.000 description 1
- 238000002922 simulated annealing Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Navigation (AREA)
Abstract
The invention provides a gradient separate zone optimization design method for a spacecraft pulse rendezvous trajectory, belongs to orbit design and optimization technologies in the field of spacecraft orbit dynamics and control, and in particular belongs to a rendezvous trajectory design technology in which a spacecraft rendezvous trajectory global optimization solution is obtained through calculation based on an interval arithmetic mathematical model. For a global issue of the method, a gradient optimization result is only used to determine an interval separation point and update a value of an upper bound of a target function, thereby preserving globalness of an interval optimization algorithm. For a large computation amount in the interval optimization algorithm, a gradient optimization method is introduced, so that the value of the upper bound of the target function can be quickly updated, and a gradient segmentation method can be used; in combination with a symbol segmentation method and an interval optimization strategy, intervals that do not comprise an optimal solution are effectively separated and removed, thereby improving operational efficiency. To solve the problem that the interval optimization algorithm has a high storage requirement, only a limited number of intervals are selected for processing in each iteration, so as to ensure that a storage space overflowing phenomenon does not occur in an operation.
Description
Technical field
The invention belongs to the Track desigh in spacecraft orbit dynamic and control field and optimisation technique, relating to a spacecraft utilizes the thrust of the jet thrust device generations such as the rocket engine carried by it and predetermined natural celestial body or another spacecraft to arrive the designing technique of the Orbit of Rendezvous of same position at one time with identical speed, particularly based on the Orbit of Rendezvous designing technique of the Spacecraft Rendezvous track global optimization solution utilizing Interval Computation mathematic(al) mode to calculate.
Background technology
In the task analysis of spacecraft, Track desigh occupies critical positions.In track design objective, a spacecraft applies acting force, the design making itself and predetermined nature celestial body or another spacecraft arrive the Orbit of Rendezvous of same position at one time with identical speed is the condition precedent that multiple space mission realizes.These tasks comprise: the maintenance in-orbit etc. of the detection of the spacecrafts rendezvous of space platform (as space station) in orbit and other spacecraft (as manned shuttle vehicle, space transoportation freighter), natural celestial body and landing, in-orbit spacecraft.The quantity realizing the possible track of intersection is infinite many in theory, so Orbit of Rendezvous design objective just becomes a track optimizing task, and the result that the track namely finding out certain Least-cost from above-mentioned infinite multiple track designs as Orbit of Rendezvous.So-called Least-cost can be that the flight time is the shortest, fuel consumption is minimum.The benefit that this optimization brings comprises the reduction of space mission cost, the raising etc. of space mission efficiency.Therefore, the design optimizing of Spacecraft Rendezvous track has very important position in aerospace engineering.
The realization of Spacecraft Rendezvous track generally can be divided into two kinds of approach.Use various rocket engines for traditional, ripe use chemical propellant, by instruction rocket engine predetermined instant carry out the short time, the igniting that can be considered pulse action, produce the acting force that spacecraft is flown by predetermined Orbit of Rendezvous.Another kind of use advanced, not yet full-fledged, that low thrust can be exported continuously various thrusters at present, in flight course all or in part, carry out long thrust output by instruction thruster, produce the acting force realizing intersection.
The Optimization Design of Spacecraft Rendezvous track can be divided into indirect method, direct method and mixing method three major types.Indirect method utilizes variational principle and Pang Te lia king maximal principle that Spacecraft Rendezvous Orbit Optimized problem is converted into the two-point boundary value problem meeting necessary condition for optimality and solves.Under given conditions, use indirect method can obtain the analytic solution of Spacecraft Rendezvous Orbit Optimized, but in most cases, become abnormal difficult by the susceptibility of the non-linear integral operation caused to association's state initial value of problem by making solving of problem.Direct method utilizes and the nonlinear programming problem that Spacecraft Rendezvous Orbit Optimized problem is converted into multidimensional is solved the discretize of quantity of state and control function and parametrization.Research about direct method uses based on the Sequential Quadratic Programming method of gradient or based on bionic intelligent optimization algorithm mostly, as genetic algorithm, simulated annealing, Artificial Immune Algorithm, particle cluster algorithm, ant group algorithm.The shortcoming of Sequential Quadratic Programming method is that it can only ensure to obtain the characteristic of locally optimal solution to the susceptibility of the initial guess of parameter and its.The shortcoming of intelligent optimization algorithm is that it cannot ensure its Global Optimality separated.Direct method and indirect method are combined use by mixing method, introduce association's state on the one hand to determine necessary condition for optimality, introduce parametric method on the other hand and solve optimal control problem.Research about mixing method often uses the method for solving identical with direct method, therefore also has the shortcoming similar to direct method research.
Above-mentionedly all belong to randomness optimization method based on bionic intelligent optimization algorithm.The ultimate principle of randomness optimization method is: from selected initial guess, by generating new conjecture value to the random perturbation of initial guess, and repeat this process until reach termination condition given in advance.When long enough operation time, randomness optimization method can obtain the solution of approaching globally optimal solution with high probability in theory.But in actual use, the Global Optimality of the operation result of randomness optimization method is undemonstrable.Relative to randomness optimization method, what deterministic optimization method was determined by structure puts by limited point or to unlimited that delimits interior the sequence formed, and by iterative process, this sequence converges is included in enough little compass in globally optimal solution or by globally optimal solution.
Interval Optimization Method is exactly a kind of determinacy global optimization method received much concern in recent years.Interval Optimization Method is based on the Study on Interval Analysis Theory proposed clearly and systematically in the works " Interval Analysis " of publication in 1966 by American scholar Ramon Edgar Moore, by the segmentation to interval, progressively get rid of the interval that can not comprise globally optimal solution, and the final list obtaining interval that contain globally optimal solution, that meet pre-provisioning request.Within 2010, Erik-Jan van Kampen is at AIAAGuidance, Navigation, the paper " Optimization of spacecraftrendezvous and docking using interval analysis " that and Control Conference delivers gives comparing of an interval Optimization Method and gradient optimizing method, result shows that gradient optimizing method converges on non-globally optimal solution when processing the non-convex optimization problem closer to reality, and interval Optimization Method still can obtain globally optimal solution.But range optimization algorithm is also because storage space is excessive needed for computing, calculated amount is large, the shortcomings such as speed of convergence is slow, pind down in practicality, Chen Cheng in 2014 give also above-mentioned viewpoint at the paper of Shanghai University's application master's degree " realizing based on the research of the parallel global optimizing algorithm of intervl mathematics and system " by name.Much research is attempted by being combined to solve the problem with other optimization methods by interval Optimization Method.The paper of keeping gentle room few pure " a kind of novel interval-particle swarm optimization algorithm " by name delivered on " Northeastern University's journal (natural science edition) " 2012 end of the year discloses a kind of method be combined with interval Optimization Method by particle cluster algorithm, but the result of calculation provided in literary composition shows, objective function codomain interval is very close but do not comprise actual optimum value.The people such as Tong Chen in 2013 are at " Journal of Guidance, Control andDynamics " on the paper of by name " Optimization of time-open constrained Lambert rendezvous usinginterval analysis " delivered indicate that interval Optimization Method storage demand is excessive, in computing, internal memory such as easily to overflow at the deficiency, and disclose a kind of for spacecraft pulse intersection Optimal design of trajectory, the method that interval Optimization Method and gradient optimizing method combine, the method first uses interval Optimization Method to try to achieve the still satisfied predetermined interval comprising global optimum required of some interval width, then use gradient optimizing method instead and carry out local optimum calculating, improve optimization efficiency, but its cost is the Global Optimality that cannot ensure optimum results.
In summary, although the Optimization Design of spacecraft pulse intersection track has a lot of existing method available, each have their own weakness, and the Global Optimality that generally cannot ensure optimum results.And can ensure that the weakness of the existing interval Optimization Method of trying to achieve globally optimal solution is mainly that operand when being applied to practical problems is large, storage demand this 2 point high in theory, the method for existing this two problems of solution is often to destroy the guarantee of Global Optimality for cost.
Summary of the invention
For the deficiency that above-mentioned prior art exists, the present invention proposes a kind of for spacecraft pulse intersection Optimal design of trajectory, based on the method for range optimization theory.
For solving the problems of the technologies described above, technical scheme of the present invention is as follows:
According to the spacecraft pulse intersection track optimizing model that adopts or set up, construct the Global Optimal Problem based on Study on Interval Analysis Theory and according to the following step, set up Global Optimal Problem solved:
Step 1, intervalization process is carried out to spacecraft pulse intersection track optimizing model, comprises following sub-step:
Sub-step 1.1, the impulse force effect times N that should use in given spacecraft pulse intersection process (N be more than or equal to 2 natural number), get last impulse force t action time
f, last impulse force t action time
fand t interval time between each subpulse time of thrust application except last impulse force
i → fcomponent (the v of the speed increment that (i=1,2...N-1), each subpulse thrust are brought to pursuit spacecraft
ix, v
iy, v
iz) (i=1,2 ..., N) as decision variable, and form decision variable interval [u]=[[v by the possible span of each decision variable
1x], [v
1y], [v
1z] ..., [v
(N-1) x], [v
(N-1) y], [v
(N-1) z], [t
1 → f], [t
2 → f] ..., [t
n-1 → f], [t
f]].
Sub-step 1.2, by the relative distance constraint S of X, Y, the Z tri-under relative coordinate system during Spacecraft Rendezvous on direction of principal axis
iu (i=x, y, z) is converted into inequality constrain in ()=0 | [S
i] ([u]) | (ε is little arithmetic number to < ε; I=x, y, z) and equation characteristic formp 0 ∈ [S
i] ([u]) (i=x, y, z).[S
i] ([u]) (i=x, y, z) represent the relative distance function in i direction;
Sub-step 1.3, the interval extension getting fuel optimal conditions is objective function
and the upper bound J of given objective function value
mininitial value;
Sub-step 1.4, decision variable interval number M, the predetermined positive number δ of each inspection of setting
uwith predetermined positive number δ
j;
Step 2, according to spacecraft pulse intersection track optimizing model, interal separation is carried out to decision-making range of variables [u], be called " symbol segmentation method ", comprise following sub-step:
Sub-step 2.1, by the component ([v of the speed increment that each subpulse thrust is brought to pursuit spacecraft
ix], [v
iy], [v
iz]) (i=1,2 ..., N) be that boundary is respectively divided into only containing on the occasion of interval with the component only containing two speed increments of negative value with 0;
Sub-step 2.2, by last impulse force [t action time
f], last [t interval time between impulse force action time and each subpulse time of thrust application except last impulse force
i → f] (i=1,2...N-1), with the quadrant of used spacecraft pulse intersection track optimizing model intermediate cam function for boundary is respectively divided into multiple time interval;
Sub-step 2.3, by the component of the speed increment of sub-step 2.1 and sub-step 2.2 gained, to obtain Q decision variable interval for interval and time interval combination in any, and Q described decision variable is interval forms interval group's queue L;
Step 3, interval selection strategy, according to the interval number M of the decision variable of each inspection of setting, if Q≤M, then make M=Q, select M decision variable interval composition sub-range queue L1 at the most from interval group's queue L assigned address, from interval group's queue L, delete M at the most the decision variable selected interval simultaneously;
Step 4, uses on each the decision variable interval in the sub-range queue L1 of the optimized algorithm based on gradient respectively selected by step 3 and solves spacecraft pulse intersection track optimizing problem, one of proceed as follows respectively by solving result:
Step 4 result 1: any one the decision variable interval [u in sub-range queue L1
□] on, the optimized algorithm based on gradient has solution, carries out the operation of following sub-step 4.1 and sub-step 4.2:
Sub-step 4.1, gradient split plot design, namely centered by the decision variable value of this solution correspondence, decision variable value expanded to [I between the expansion area of the solution of the optimized algorithm included based on gradient for radius to the numerical value of sizing
Δ], with [I between expansion area
Δ] the boundary value of each interval variable split described decision variable interval [u respectively
□] each corresponding interval variable, by decision variable interval [u
□] each interval variable be divided into 3 intervals or 2 intervals, and combined as interval queue Lnew;
Sub-step 4.2, the upper bound based on the objective function value of gradient upgrades, [I between the expansion area namely decision variable value being expanded to the solution of the optimized algorithm included based on gradient centered by the decision variable value of this solution correspondence, with very little numerical value for radius
Δ 1], utilize Study on Interval Analysis Theory [I between expansion area
Δ 1] upper calculating, check whether the interval upper bound of the objective function of its constraint condition whether meeting spacecraft pulse intersection track optimizing model, its correspondence is less than the upper bound J of objective function value
min, when above-mentioned two results checked are all "Yes", with [I between expansion area
Δ 1] upper bound in corresponding objective function interval upgrades the upper bound J of objective function value
min, otherwise keep the upper bound J of existing objective function value
minconstant;
Step 4 result 2: any one the decision variable interval [u in sub-range queue L1
□] on, the optimized algorithm based on gradient is not separated, with decision variable interval [u
□] on the interval midpoint of each interval variable split each interval variable respectively, by decision variable interval [u
□] each interval variable be divided into 2 intervals, and combined as interval queue Lnew;
Step 5, interval deletion strategy, namely carries out interval analysis respectively, and carries out the operation of following sub-step 5.1 ~ sub-step 5.4 on each decision variable interval of the interval queue Lnew of step 4 gained:
Sub-step 5.1, if relative distance constraint S
iu in ()=0 (i=x, y, z), each impulse speed is linear correlation, then utilize the decision variable except the component of the speed increment brought to pursuit spacecraft except jth subpulse thrust in decision variable interval [u] to pass through
[S
i] to solve the component of the speed increment that jth subpulse thrust is brought to spacecraft interval in ([u])=0 (i=x, y, z)
[[v
jx_f], [v
jy_f], [v
jz_f]] (j=1,2...N-1), and calculate itself and [[v in [u]
jx], [v
jy], [v
jz]] common factor
[[v
jx_new], [v
jy_new], [v
jz_new]], if this common factor is empty set, then this decision variable interval [u] is deleted from interval queue Lnew, otherwise the component interval upgrading the speed increment that the jth subpulse thrust in decision variable interval [u] is brought to pursuit spacecraft is [[v
jx_new], [v
jy_new], [v
jz_new]], and obtain the interval queue Lnew of renewal simultaneously;
Sub-step 5.2, checks whether each the decision variable interval in the interval queue Lnew after step 5.1 upgrades meets relative distance constraint condition, and is deleted from interval queue Lnew in the decision variable interval not meeting relative distance constraint condition;
Sub-step 5.3, appoints and gets 3 speed increment variablees corresponding to pulsatile once effect in decision variable, according to relative distance constraint S
i(u)=0 (i=x, y, z), 3 speed increment variablees of specifying can be expressed as the function of other decision variables, be the spacecraft pulse intersection track optimizing model that is decision variable with its dependent variable except 3 named variables by spacecraft pulse intersection track optimizing model conversation with this, check and whether each the decision variable interval in the interval queue Lnew of objective function J after processing through step 5.2 comprises 0 to the first-order partial derivative interval of new decision variable, if check result is no, then corresponding decision variable interval is deleted from interval queue Lnew;
Sub-step 5.4, to the interval upper bound J checking objective function interval and objective function value respectively of each decision variable in the interval queue Lnew after sub-step 5.3 process
minbetween relation, and carry out the operation of following grandson step 5.4.1 and grandson step 5.4.2:
Grandson step 5.4.1, if the lower bound in objective function interval is greater than the upper bound J of objective function value
min, then corresponding decision variable interval is deleted from interval queue Lnew;
Grandson step 5.4.2, if the upper bound in objective function interval is less than the upper bound J of objective function value
min, then the upper bound J of objective function value is upgraded
minfor the upper bound in objective function interval;
Step 6, interval deflation strategy, namely checks relative distance function [S
i] in ([u]) (i=x, y, z) interval queue Lnew after processing through step 5 on each decision variable interval [u] to [u] in the monotonicity of each interval variable, and carry out the operation of following sub-step 6.1:
Sub-step 6.1, if interval to any one decision variable in interval queue Lnew, relative distance function [S
i] ([u]) to the interval variable [u of the kth in decision-making range of variables [u]
k] be dull, and [u
k] lower bound place [S
i] ([u
k] .inf) or [u
k] upper bound place
[S
i] ([u
k] .sup) and symbol be just or be negative, then at [u
k] in find be greater than [u
k] lower bound and the [S of correspondence
i] and [u
k] lower bound place
[S
i] ([u
k] .inf) u of jack per line
k1if do not have, then make u
k1=[u
k] .inf; Searching is less than [u
k] upper bound and the [S of correspondence
i] and [u
k] on
Place of boundary [S
i] ([u
k] .sup) u of jack per line
k2if do not have, then make u
k2=[u
k] .sup; With interval [u
k1, u
k2] upgrade the kth interval variable [u of the corresponding decision range of variables [u] in this interval queue Lnew
k], wherein, [u
k] .inf represents [u
k] lower bound, [u
k] .sup represents [u
k] the upper bound;
Step 7, the interval end condition of each decision variable between test zone in queue Lnew, the width B in the width A in each the decision variable interval namely between test zone in queue Lnew and the objective function interval corresponding to each the decision variable interval in interval queue Lnew, wherein the width in decision variable interval is the maximal value of the width of all interval variables in decision variable interval, performs the operation of following sub-step 7.1 and sub-step 7.2:
Sub-step 7.1, when above-mentioned width A is less than predetermined positive number δ
u, or width B is less than predetermined positive number δ
jtime, corresponding decision variable interval is deleted from interval queue Lnew, and described corresponding decision variable interval is inserted in the interval queue R of design result;
Sub-step 7.2, after the operation to the interval completing steps 7.1 of each decision variable in interval queue Lnew, upgrade interval group's queue L, the deletion point that interval group's queue L deletes M at the most interval of selection is in step 3 inserted into, using the interval group's queue L after queue Lnew between insert district as new interval group's queue L by interval queue Lnew;
Step 8, the number in decision variable interval in group's queue L between test zone, and carry out one of following operation respectively according to result:
Step 8 result 1, in interval group's queue L, the number in decision variable interval is not 0, then proceed to step 3, continues design operation;
Step 8 result 2, in interval group's queue L, the number in decision variable interval is 0, then design end;
Step 9, appoint from the interval queue R of design result and get a design result interval, in the interval of each interval variable in selected design result interval, any value combines, and the speed increment that the last impulse force effect accordingly that calculates brings to pursuit spacecraft, just obtain an optimal design solution of aforementioned spacecraft pulse intersection track.
Preferably, assigned address described in the step 3 of said method is the tail of the queue of interval group's queue L, in the sub-step 7.2 corresponded, interval queue Lnew is placed in interval group's queue L tail of the queue to form new interval group's queue L.
Preferably, assigned address described in the step 3 of said method is the head of the queue of interval group's queue L, in the sub-step 7.2 corresponded, interval group's queue L is placed in interval queue Lnew tail of the queue to form new interval group's queue L.
Two kinds of preferred versions of above-mentioned steps 3 make the deletion of interval group's queue L and renewal rewards theory more convenient.
Preferably, at [u in the sub-step 6.1 of said method
k] in find be greater than [u
k] lower bound and the [S of correspondence
i] and [u
k] lower bound place [S
i] ([u
k] .inf) the maximal value u of jack per line
k3if do not have, then make u
k3=[u
k] .inf, find and be less than [u
k] upper bound and the [S of correspondence
i] and [u
k] upper bound place [S
i] ([u
k] .sup) jack per line minimum value u
k4if do not have, then make u
k4=[u
k] .sup, with interval [u
k3, u
k4] upgrade the kth interval variable [u of the corresponding decision range of variables [u] in this interval queue Lnew
k].
The preferred version of above-mentioned sub-step 6.1 will obtain better interval deflation result.
Preferably, in the step 9 of said method, the upper bound of getting a corresponding objective function interval from design result interval queue R is the upper bound J of objective function value
mindesign result interval, in selected design result interval, in the interval of each decision variable, any value combines, and the speed increment that the last impulse force effect accordingly that calculates brings to pursuit spacecraft, just obtain an optimal design solution of aforementioned spacecraft pulse intersection track.
The preferred version of above-mentioned steps 9 by the target function value that ensures corresponding to design result not higher than J
min.
Compared with prior art, beneficial effect of the present invention is the Global Optimality that ensure that gained design result.In the present invention, can not ensure that because depending on Initial value choice the optimization method based on gradient obtaining globally optimal solution is used in step 4.But based on the operation result of the optimization method of gradient only for determining interal separation point and upgrading the upper bound of objective function value, the present invention is a kind of interval Optimization Method completely, and therefore inherits the Global Optimality of interval Optimization Method.
Compared with prior art, beneficial effect of the present invention is also that to the maximal value of the demand in Computer Storage space be pre-determined, and this will cause the abnormal situation of interrupting of design process to occur to thus avoid memory space inadequate (spilling) that range optimization algorithm may occur.This beneficial effect provides primarily of the step 3 in technique scheme.According to step 3, each limited individual (maximum M) interval of only selecting from pending interval group processes, and M value can be arranged according to the storage space of used computing machine, and therefore can ensure not occur storage space spillover in design process.
Compared with prior art, beneficial effect of the present invention is also the raising of design efficiency, namely obtains the minimizing of design result required time.This beneficial effect provides primarily of the step 2 in technique scheme and step 4-step 6.According to step 2, the whole region of search is divided into multiple interval, avoid because interval comprises 0 and the too extensive interval extension caused, be conducive to carrying out in early days and the quick minimizing of consequent pending interval quantity of interval deletion action, and thereby reduce design process required time.According to step 4, based on the use of the operation result of the optimization method of gradient, be conducive to utilizing the information of locally optimal solution or globally optimal solution to upgrade the upper bound of objective function value fast.According to step 4, the optimization method computing acquired results based on gradient is utilized to carry out interal separation, be conducive to obtaining fast the decision variable interval having and comprise local optimum and there is less interval width, be also conducive to other and do not comprise giving up fast of optimum solution interval.These two advantages provided by step 4 are obviously conducive to the minimizing of design process required time.According to the operation of step 5, be conducive to reducing the number in pending interval and upgrading the upper bound of objective function value, and be therefore conducive to the minimizing of design process required time.According to step 6, local monotonicity in constraint condition is used to the width reducing pending interval further, be conducive to the reduction of interval extension in interval analysis computing, and be therefore conducive to carrying out in early days of interval deletion action, thereby produce the effect reducing design process required time.
Accompanying drawing explanation
Fig. 1 be spacecraft pulse intersection track provided by the invention gradient cut section between the schematic flow sheet of Optimization Design;
Fig. 2 is the schematic diagram that the step 3-step 7 of method provided by the invention performs interval group's queue change in the process of twice continuously;
Fig. 3 is the renewal process schematic diagram based on the upper bound of gradient objective function value in method provided by the invention;
Fig. 4 be gained according to one embodiment of present invention two spacecrafts between the schematic diagram of relative movement orbit;
Fig. 5 is the schematic diagram of the interval queue of Optimum Design Results of gained according to one embodiment of present invention.
Embodiment
Below in conjunction with accompanying drawing and embodiment, the present invention is described in further detail.
For the dipulse spacecrafts rendezvous optimization problem of on-fixed time, relative coordinate system X-axis in orbit plane along passive space vehicle heading, Z axis points to the earth's core from passive space vehicle barycenter, Y-axis meets right hand orthogonal coordinate system, perpendicular to orbit plane, and following CW equation approximate description is adopted to set up spacecraft pulse intersection track optimizing model:
In formula, ω is passive space vehicle orbit angular velocity, and x, y, z is respectively pursuit spacecraft relative position in the X, Y, Z direction,
be respectively pursuit spacecraft relative velocity in the X, Y, Z direction,
be respectively pursuit spacecraft relative acceleration in the X, Y, Z direction.
Through type (1), can obtain t pursuit spacecraft relative to the relative position of passive space vehicle and speed as follows:
In formula, (x
0, y
0, z
0) and
be respectively initial position and the initial velocity in pursuit spacecraft 0 moment.
If the moment of pursuit spacecraft and passive space vehicle intersection is t
f.From 0 moment, at t
1in the moment, pursuit spacecraft applies first time pulse, and impulse force effect is v to the speed increment that pursuit spacecraft brings
1(v
1x, v
1y, v
1z).Through t
1 → f=t
f-t
1time, at t
fin the moment, pursuit spacecraft applies second time pulse, and impulse force effect is v to the speed increment that pursuit spacecraft brings
2(v
2x, v
2y, v
2z), two Spacecraft Rendezvous.
Through type (2) can obtain intersection moment t
frelative distance S and relative velocity V, and make S=0 and V=0.
S(t
1,t
f,v
1)=0
(3)
V(t
1,t
f,v
1,v
2)=0
Then (3) formula is with t
1, t
f, v
1, v
2for the constraint condition of the spacecraft pulse intersection track optimizing model of decision variable, be respectively relative distance constraint and relative velocity constraint.
If only consider the requirement to fuel consumption, then the objective function of spacecraft pulse intersection track optimizing model is:
J gets minimum value, then dipulse intersection fuel consumption is minimum.
Adopt Optimization Design between the gradient cut section of spacecraft pulse intersection track to solve above Optimized model, between the gradient cut section of spacecraft pulse intersection track, the schematic flow sheet of Optimization Design, as shown in Figure 1, specifically comprises the steps:
Step 1, intervalization process is carried out to spacecraft pulse intersection track optimizing model, comprises following sub-step:
Sub-step 1.1, assuming that two spacecrafts move in the orbit plane of passive space vehicle, then has decision variable interval
[u]=[[v
1x],[v
1z],[t
1→f],[t
f]]。
Sub-step 1.2, the relative velocity constraint V=0 speed increment (v that can be brought to pursuit spacecraft by the effect of second time impulse force at Spacecraft Rendezvous place
2x, v
2y, v
2z) ensured.And for relative distance constraint S=0, because two spacecrafts move in the orbit plane of passive space vehicle, so S=0 can be decomposed into S
x=0, S
z=0.By S
x=0, S
z=0 is converted into inequality constrain form | [S
x] ([u]) | < ε, | [S
z] ([u]) | < ε (ε is little arithmetic number) and equation characteristic formp 0 ∈ [S
x] ([u]), 0 ∈ [S
z] ([u]).[S
x] ([u]) represent the relative distance function in x direction, [S
z] ([u]) represent the relative distance function in z direction.
Sub-step 1.3, the interval extension getting fuel optimal conditions is objective function
And the upper bound J of given objective function value
mininitial value; N=2.
Sub-step 1.4, decision variable interval number M, the predetermined positive number δ of each inspection of setting
uwith predetermined positive number δ
j;
Step 2, uses symbol segmentation method, comprises following sub-step:
Sub-step 2.1, by the component [v of the speed increment that first time impulse force effect brings to pursuit spacecraft
1x] and [v
1z] be that boundary is respectively divided into only containing on the occasion of interval with the component only containing two speed increments of negative value with 0;
Sub-step 2.2, will second time impulse force [t action time
f], [t interval time between second time impulse force action time and impulse force action time first time
1 → f], with the quadrant of used spacecraft pulse intersection track optimizing model intermediate cam function for boundary is respectively divided into multiple time interval;
Sub-step 2.3, by the component of the speed increment of sub-step 2.1 and sub-step 2.2 gained, to obtain Q decision variable interval for interval and time interval combination in any, and Q described decision variable is interval forms interval group's queue L;
Step 3, according to the interval number M of the decision variable of each inspection of setting, if Q≤M, then make M=Q, M decision variable interval composition sub-range queue L1 is at the most selected from interval group's queue L tail of the queue (or optional position), from interval group's queue L, delete M at the most the decision variable selected interval, step 3 ~ step 7 performs the signal of interval group's queue change in the process of twice continuously, as shown in Figure 2 simultaneously;
Step 4, uses on each the decision variable interval in the sub-range queue L1 of the optimized algorithm based on gradient respectively selected by step 3 and solves spacecraft pulse intersection track optimizing problem, one of proceed as follows respectively by solving result:
Step 4 result 1: any one the decision variable interval [u in sub-range queue L1
□] on, the optimized algorithm based on gradient has solution, carries out the operation of following sub-step 4.1 and sub-step 4.2:
Sub-step 4.1, centered by the decision variable value of this solution correspondence, decision variable value expanded to [I between the expansion area of the solution of the optimized algorithm included based on gradient for radius to the numerical value of sizing
Δ], with [I between expansion area
Δ] the boundary value of each interval variable split described decision variable interval [u respectively
□] each corresponding interval variable, by decision variable interval [u
□] each interval variable be divided into 3 intervals or 2 intervals, and combined as interval queue Lnew;
Sub-step 4.2, based on the renewal in the objective function value upper bound of gradient, [I between the expansion area namely decision variable value being expanded to the solution of the optimized algorithm included based on gradient centered by the decision variable value of this solution correspondence, with very little numerical value for radius
Δ 1], utilize Study on Interval Analysis Theory [I between expansion area
Δ 1] on calculate, check whether the interval upper bound of the objective function of its constraint condition whether meeting spacecraft pulse intersection track optimizing model, its correspondence is less than the upper bound J of objective function value
min, when above-mentioned two results checked are all "Yes", with [I between expansion area
Δ 1] upper bound in corresponding objective function interval upgrades the upper bound J of objective function value
min, otherwise keep the upper bound J of existing objective function value
minconstant, process schematic as shown in Figure 3;
Step 4 result 2: the decision variable interval [u in sub-range queue L1
□] on, the optimized algorithm based on gradient is not separated, with decision variable interval [u
□] on each interval variable mid point split each interval variable respectively, by decision variable interval [u
□] each interval variable be divided into 2 intervals, and combined as interval queue Lnew;
Step 5, interval deletion strategy, namely carries out interval analysis respectively, and carries out the operation of following sub-step 5.1 ~ sub-step 5.4 on each decision variable interval of the interval queue Lnew of step 4 gained:
Sub-step 5.1, if relative distance constraint S
iu in ()=0 (i=x, y, z), each impulse speed is linear correlation, then utilize interval variable [t in decision variable interval [u]
1 → f] and [t
f], by [S
i] ([u])=0 (i=x, z) solve component the interval [[v of the speed increment that first time impulse force effect brings to pursuit spacecraft
1x_f], [v
1z_f]], and calculate itself and [[v in [u]
1x], [v
1z]] common factor
[[v
1x_new], [v
1z_new]], if this common factor is empty set, then this decision variable interval [u] is deleted from interval queue Lnew, otherwise the component interval upgrading the speed increment that the first time impulse force effect in decision variable interval [u] brings to pursuit spacecraft is [[v
1x_new], [v
1z_new]], and obtain the interval queue Lnew of renewal simultaneously;
Sub-step 5.2, checks whether each the decision variable interval in the interval queue Lnew after step 5.1 upgrades meets relative distance constraint condition, and is deleted from interval queue Lnew in the decision variable interval not meeting relative distance constraint condition;
Sub-step 5.3, appoints and gets 3 speed increment variablees corresponding to pulsatile once effect in decision variable, according to relative distance constraint S
i(u)=0 (i=x, y, z), 3 speed increment variablees of specifying can be expressed as the function of other decision variables, be the spacecraft pulse intersection track optimizing model that is decision variable with its dependent variable except 3 named variables by spacecraft pulse intersection track optimizing model conversation with this, check and whether each the decision variable interval in the interval queue Lnew of objective function J after processing through sub-step 5.2 comprises 0 to the first-order partial derivative interval of new decision variable, if check result is no, then corresponding decision variable interval is deleted from interval queue Lnew,
Sub-step 5.4, to the interval upper bound J checking objective function interval and objective function value respectively of each decision variable in the interval queue Lnew after sub-step 5.3 process
minbetween relation, and carry out the operation of following grandson step 5.4.1 and grandson step 5.4.2:
Grandson step 5.4.1, if the lower bound in objective function interval is greater than the upper bound J of objective function value
min, then corresponding decision variable interval is deleted from interval queue Lnew;
Grandson step 5.4.2, if the upper bound in objective function interval is less than the upper bound J of objective function value
min, then the upper bound J of objective function value is upgraded
minfor the upper bound in objective function interval;
Step 6, interval deflation strategy, namely checks relative distance function [S
i] on each decision variable interval [u] in ([u]) (i=x, z) interval queue Lnew after processing through step 5 to [u] in each interval variable monotonicity, and carry out the operation of following sub-step 6.1:
Sub-step 6.1, if interval to any one decision variable in interval queue Lnew, relative distance function [S
i] ([u]) to the interval variable [u of the kth in decision-making range of variables [u]
k] be dull, and [u
k] lower bound place [S
i] ([u
k] .inf) or [u
k] upper bound place
[S
i] ([u
k] .sup) and symbol be just or be negative, then at [u
k] in find be greater than [u
k] lower bound and the [S of correspondence
i] and [u
k] lower bound place
[S
i] ([u
k] .inf) the maximal value u of jack per line
k1if do not have, then make u
k1=[u
k] .inf; Searching is less than [u
k] upper bound and the [S of correspondence
i] with
[u
k] upper bound place [S
i] ([u
k] .sup) the minimum value u of jack per line
k2if do not have, then make u
k2=[u
k] .sup; With interval [u
k1, u
k2] upgrade the kth interval variable [u of the corresponding decision range of variables [u] in this interval queue Lnew
k], wherein, [u
k] .inf represents [u
k] lower bound, [u
k] .sup represents [u
k] the upper bound;
Step 7, the interval end condition of each decision variable between test zone in queue Lnew, the width B in the width A in each the decision variable interval namely between test zone in queue Lnew and the objective function interval corresponding to each the decision variable interval in interval queue Lnew, wherein the width in decision variable interval is the maximal value of the width of all interval variables in decision variable interval, performs the operation of following sub-step 7.1 and sub-step 7.2:
Sub-step 7.1, when above-mentioned width A is less than predetermined positive number δ
u, or width B is less than predetermined positive number δ
jtime, corresponding decision variable interval is deleted from interval queue Lnew, and described corresponding decision variable interval is inserted in the interval queue R of design result;
Sub-step 7.2, after the operation to the interval completing steps 7.1 of each decision variable in interval queue Lnew, upgrade interval group's queue L, namely interval queue Lnew is inserted into interval group's queue L tail of the queue;
Step 8, the number in decision variable interval in group's queue L between test zone, and carry out one of following operation respectively according to result:
Step 8 result 1, in interval group's queue L, the number in decision variable interval is not 0, then proceed to step 3, continues design operation;
Step 8 result 2, in interval group's queue L, the number in decision variable interval is 0, then design end.
Step 9, chooses Optimum Design Results, gets the upper bound J that the interval upper bound of the interval corresponding objective function of a design result is objective function value from the interval queue R of design result
mindesign result interval, in selected design result interval, in the interval of each decision variable, any value combines, and the speed increment that the last impulse force effect accordingly that calculates brings to spacecraft, just obtain an optimal design solution of aforementioned spacecraft pulse intersection track.
The design effect of method provided by the invention is shown further below in conjunction with a concrete example.
If passive space vehicle is on the circular orbit that 400km is high, relative distance precision ε is 0.015m, δ
uand δ
jbe 0.01, if initial decision range of variables [u]=[[-30,30], [-30,30], [0,8000], [0,8000]], objective function estimated value initial value is 100ms
-1, pursuit spacecraft initial position is (-10000,0,8000) m, and initial velocity is (10,0 ,-15) ms
-1.
Simulation calculating show that the upper bound of objective function value and the solution interval of correspondence are as shown in table 1 under current accuracy
The upper bound J of table 1 objective function value
mincorresponding solution is interval
In this group decision variable interval, any value, its target function value all in [J] interval, the wherein upper bound J of objective function value
minthe higher limit being [J] is 8.750384691131890ms
-1.
As shown in Figure 4, As time goes on relative distance, be progressively reduced to 0, and error is within the scope of the 0.015m specified.From path, demonstrate the correctness that the present invention solves intuitively.
As shown in Figure 5, within the scope of designated precision, the feasible solution interval that may comprise optimization solution satisfied condition has 5520, each [t
f] corresponding to [J] all contain the upper bound of objective function value.
Adopt this problem of genetic algorithm for solving, decision variable is t
1 → fand t
f, hunting zone is [0,8000] s, if initial population is 100, in maximum genetic algebra 100 generation, solve 10 times, result is as table 2.
Table 2 genetic algorithm result table
As shown in Table 2, genetic algorithm, due to algorithmic characteristic, depends on Initial value choice, is easily absorbed in local optimum, is unable to undergo computing repeatedly.In table 2,10 genetic algorithm objective function optimal values are the objective function optimization value 8.750384691363504 of the 6th computing, the upper dividing value 8.750384691131890ms of the objective function value of the global optimization approach proposed with the present invention
-1very nearly the same, prove the correctness that the present invention solves, also provide support for advantage of overall importance of the present invention.
Claims (5)
1. an Optimization Design between the gradient cut section of spacecraft pulse intersection track, is characterized in that, said method comprising the steps of:
Step 1, intervalization process is carried out to spacecraft pulse intersection track optimizing problem, comprises following sub-step:
Sub-step 1.1, the impulse force effect times N that should use in given spacecraft pulse intersection process, gets last impulse force t action time
f, last impulse force t action time
fand t interval time between each subpulse time of thrust application except last impulse force
i → f, the component (v of speed increment that brings to pursuit spacecraft of each subpulse thrust
ix, v
iy, v
iz) as decision variable, and it is interval to form decision variable by the possible span of each decision variable
[u]=[[v
1x], [v
1y], [v
1z] ..., [v
(N-1) x], [v
(N-1) y], [v
(N-1) z], [t
1 → f], [t
2 → f] ..., [t
n-1 → f], [t
f]]; N be more than or equal to 2 natural number; Interval time t
i → fmiddle i=1,2...N-1, the component (v of speed increment
ix, v
iy, v
iz) middle i=1,2 ..., N;
Sub-step 1.2, by the relative distance constraint S of X, Y, the Z tri-under relative coordinate system during Spacecraft Rendezvous on direction of principal axis
iu ()=0 is converted into inequality constrain | [S
i] ([u]) | < ε and equation characteristic formp 0 ∈ [S
i] ([u]); ε is little arithmetic number, i=x, y, z;
Sub-step 1.3, the interval extension getting fuel optimal conditions is objective function
and the upper bound J of given objective function value
mininitial value;
Sub-step 1.4, decision variable interval number M, the predetermined positive number δ of each inspection of setting
uwith predetermined positive number δ
j;
Step 2, according to spacecraft pulse intersection track optimizing model, interal separation is carried out to decision-making range of variables [u], comprise following sub-step:
Sub-step 2.1, by the component ([v of the speed increment that each subpulse thrust is brought to pursuit spacecraft
ix], [v
iy], [v
iz]) be that boundary is respectively divided into only containing on the occasion of interval with the component only containing two speed increments of negative value with 0; Wherein i=1,2 ..., N;
Sub-step 2.2, by last impulse force [t action time
f], last [t interval time between impulse force action time and each subpulse time of thrust application except last impulse force
i → f], i=1,2...N-1, with the quadrant of used spacecraft pulse intersection track optimizing model intermediate cam function for boundary is respectively divided into multiple time interval;
Sub-step 2.3, by the component of the speed increment of sub-step 2.1 and sub-step 2.2 gained, to obtain Q decision variable interval for interval and time interval combination in any, and Q described decision variable is interval forms interval group's queue L;
Step 3, according to the interval number M of the decision variable of each inspection of setting, if Q≤M, then make M=Q, select M interval composition sub-range queue L1 at the most from interval group's queue L assigned address, from interval group's queue L, delete M at the most the decision variable selected interval simultaneously;
Step 4, uses on each the decision variable interval in the sub-range queue L1 of the optimized algorithm based on gradient respectively selected by step 3 and solves spacecraft pulse intersection track optimizing problem, one of proceed as follows respectively by solving result:
Step 4 result 1: any one the decision variable interval [u in sub-range queue L1
□] on, the optimized algorithm based on gradient has solution, carries out the operation of following sub-step 4.1 and sub-step 4.2:
Sub-step 4.1, gradient split plot design, namely centered by the decision variable value of this solution correspondence, decision variable value expanded to [I between the expansion area of the solution of the optimized algorithm included based on gradient for radius to the numerical value of sizing
Δ], with [I between expansion area
Δ] the boundary value of each interval variable split described decision variable interval [u respectively
□] each corresponding interval variable, by decision variable interval [u
□] each interval variable be divided into 3 intervals or 2 intervals, and combined as interval queue Lnew;
Sub-step 4.2, [I between expansion area decision variable value being expanded to the solution of the optimized algorithm included based on gradient centered by the decision variable value of this solution correspondence, with very little numerical value for radius
Δ 1], utilize Study on Interval Analysis Theory [I between expansion area
Δ 1] upper calculating, check whether the interval upper bound of the objective function of its constraint condition whether meeting spacecraft pulse intersection track optimizing model, its correspondence is less than the upper bound J of objective function value
min, when above-mentioned two results checked are all "Yes", with [I between expansion area
Δ 1] the interval upper bound of corresponding objective function upgrades the upper bound J of objective function value
min, otherwise keep the upper bound J of existing objective function value
minconstant;
Step 4 result 2: any one the decision variable interval [u in sub-range queue L1
□] on, the optimized algorithm based on gradient is not separated, with decision variable interval [u
□] on the interval midpoint of each interval variable split each interval variable respectively, by decision variable interval [u
□] each interval variable be divided into 2 intervals, and combined as interval queue Lnew;
Step 5, interval analysis is carried out respectively in each decision variable interval of the interval queue Lnew of step 4 gained, and carries out the operation of following sub-step 5.1 ~ sub-step 5.4:
Sub-step 5.1, if relative distance constraint S
iu in ()=0, each impulse speed is linear correlation, i=x, y, z, then the decision variable except the component utilizing the speed increment brought to pursuit spacecraft except jth subpulse thrust in decision variable interval [u] is by [S
i] ([u])=0 solves component the interval [[v of the speed increment that jth subpulse thrust is brought to spacecraft
jx_f], [v
jy_f], [v
jz_f]], i=x, y, z, j=1,2...N-1, and calculate itself and [[v in [u]
jx], [v
jy], [v
jz]] common factor
[[v
jx_new], [v
jy_new], [v
jz_new]], if this common factor is empty set, then this decision variable interval [u] is deleted from interval queue Lnew, otherwise the component interval upgrading the speed increment that the jth subpulse thrust in decision variable interval [u] is brought to pursuit spacecraft is [[v
jx_new], [v
jy_new], [v
jz_new]], and obtain the interval queue Lnew of renewal simultaneously;
Sub-step 5.2, checks whether each the decision variable interval in the interval queue Lnew after step 5.1 upgrades meets relative distance constraint condition, and is deleted from interval queue Lnew in the decision variable interval not meeting relative distance constraint condition;
Sub-step 5.3, appoints and gets 3 speed increment variablees corresponding to pulsatile once effect in decision variable, according to relative distance constraint S
i(u)=0, i=x, y, z, 3 speed increment variablees of specifying are expressed as the function of other decision variables, be the spacecraft pulse intersection track optimizing model that is decision variable with its dependent variable except 3 named variables by spacecraft pulse intersection track optimizing model conversation with this, check and whether each the decision variable interval in the interval queue Lnew of objective function J after processing through step 5.2 comprises 0 to the first-order partial derivative interval of new decision variable, if check result is no, then corresponding decision variable interval is deleted from interval queue Lnew;
Sub-step 5.4, to the interval upper bound J checking objective function interval and objective function value respectively of each decision variable in the interval queue Lnew after sub-step 5.3 process
minbetween relation, and carry out the operation of following grandson step 5.4.1 and grandson step 5.4.2:
Grandson step 5.4.1, if the lower bound in objective function interval is greater than the upper bound J of objective function value
min, then corresponding decision variable interval is deleted from interval queue Lnew;
Grandson step 5.4.2, if the upper bound in objective function interval is less than the upper bound J of objective function value
min, then the upper bound J of objective function value is upgraded
minfor the upper bound in objective function interval;
Step 6, interval deflation strategy, namely checks relative distance function [S
i] in ([u]) interval queue Lnew after processing through step 5 on each decision variable interval [u] to [u] in the monotonicity of each interval variable, i=x, y, z, and the operation carrying out following sub-step 6.1:
Sub-step 6.1, if interval to any one decision variable in interval queue Lnew, relative distance function [S
i] ([u]) to the interval variable [u of the kth in decision-making range of variables [u]
k] be dull, and [u
k] lower bound place [S
i] ([u
k] .inf) or [u
k] upper bound place
[S
i] ([u
k] .sup) and symbol be just or be negative, then at [u
k] in find be greater than [u
k] lower bound and the [S of correspondence
i] and [u
k] lower bound place
[S
i] ([u
k] .inf) u of jack per line
k1if do not have, then make u
k1=[u
k] .inf, find and be less than [u
k] upper bound and the [S of correspondence
i] and [u
k] on
Place of boundary [S
i] ([u
k] .sup) u of jack per line
k2if do not have, then make u
k2=[u
k] .sup, with interval [u
k1, u
k2] upgrade the kth interval variable [u of the corresponding decision range of variables [u] in this interval queue Lnew
k], wherein, [u
k] .inf represents [u
k] lower bound, [u
k] .sup represents [u
k] the upper bound;
Step 7, the interval end condition of each decision variable between test zone in queue Lnew, the width B in the width A in each the decision variable interval namely between test zone in queue Lnew and the objective function interval corresponding to each the decision variable interval in interval queue Lnew, wherein the width in decision variable interval is the maximal value of the width of all interval variables in decision variable interval, performs the operation of following sub-step 7.1 and sub-step 7.2:
Sub-step 7.1, when above-mentioned width A is less than predetermined positive number δ
u, or width B is less than predetermined positive number δ
jtime, corresponding decision variable interval is deleted from interval queue Lnew, and described corresponding decision variable interval is inserted in the interval queue R of design result;
Sub-step 7.2, after the operation to the interval completing steps 7.1 of each decision variable in interval queue Lnew, upgrade interval group's queue L, the deletion point that interval group's queue L deletes M at the most interval of selection is in step 3 inserted into, using the interval group's queue L after queue Lnew between insert district as new interval group's queue L by interval queue Lnew;
Step 8, the number in decision variable interval in group's queue L between test zone, and carry out one of following operation respectively according to result:
Step 8 result 1, in interval group's queue L, the number in decision variable interval is not 0, then proceed to step 3, continues design operation;
Step 8 result 2, in interval group's queue L, the number in decision variable interval is 0, then design end;
Step 9, appoint from the interval queue R of design result and get a design result interval, in the interval of each interval variable in selected design result interval, any value combines, and the speed increment that the last impulse force effect accordingly that calculates brings to pursuit spacecraft, just obtain an optimal design solution of aforementioned spacecraft pulse intersection track.
2. Optimization Design between the gradient cut section of a kind of spacecraft pulse intersection track according to claim 1, it is characterized in that, assigned address described in step 3 is the tail of the queue of interval group's queue L, accordingly, in sub-step 7.2, interval queue Lnew is placed in interval group's queue L tail of the queue to form new interval group's queue L.
3. Optimization Design between the gradient cut section of a kind of spacecraft pulse intersection track according to claim 1, it is characterized in that, assigned address described in step 3 is the head of the queue of interval group's queue L, accordingly, in sub-step 7.2, interval group's queue L is placed in interval queue Lnew tail of the queue to form new interval group's queue L.
4. according to the Optimization Design of the spacecraft pulse intersection track in claims 1 to 3 described in any one claim, it is characterized in that, in step 6.1, at a kth interval variable [u
k] in find be greater than [u
k] lower bound and the [S of correspondence
i] and [u
k] lower bound place [S
i] ([u
k] .inf) the maximal value u of jack per line
k3if do not have, then make u
k3=[u
k] .inf, find and be less than [u
k] upper bound and the [S of correspondence
i] and [u
k] upper bound place [S
i] ([u
k] .sup) jack per line minimum value u
k4if do not have, then make u
k4=[u
k] .sup, with interval [u
k3, u
k4] upgrade the kth interval variable [u of the corresponding decision range of variables [u] in this interval queue Lnew
k].
5., according to the Optimization Design of the spacecraft pulse intersection track in claims 1 to 3 described in any one claim, it is characterized in that, step 9 is as follows:
The upper bound of getting a corresponding objective function interval from design result interval queue R is the upper bound J of objective function value
mindesign result interval, in the interval of each interval variable in selected design result interval, any value combines, and the speed increment that the last impulse force effect accordingly that calculates brings to spacecraft, just obtain an optimal design solution of spacecraft pulse intersection track.
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510400314.5A CN105005651B (en) | 2015-07-09 | 2015-07-09 | Optimization Design between the gradient cut section of spacecraft pulse intersection track |
PCT/CN2016/082052 WO2017005052A1 (en) | 2015-07-09 | 2016-05-13 | Optimization and design method for gradient segmentation of intervals of spacecraft pulse rendezvous trajectory |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510400314.5A CN105005651B (en) | 2015-07-09 | 2015-07-09 | Optimization Design between the gradient cut section of spacecraft pulse intersection track |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105005651A true CN105005651A (en) | 2015-10-28 |
CN105005651B CN105005651B (en) | 2016-04-27 |
Family
ID=54378327
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510400314.5A Expired - Fee Related CN105005651B (en) | 2015-07-09 | 2015-07-09 | Optimization Design between the gradient cut section of spacecraft pulse intersection track |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN105005651B (en) |
WO (1) | WO2017005052A1 (en) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017005052A1 (en) * | 2015-07-09 | 2017-01-12 | 北京航空航天大学 | Optimization and design method for gradient segmentation of intervals of spacecraft pulse rendezvous trajectory |
CN106407619A (en) * | 2016-11-16 | 2017-02-15 | 哈尔滨工业大学 | Linear-feedback global stabilization method for controlling limited spacecraft rendezvous control system |
CN107526368A (en) * | 2017-09-12 | 2017-12-29 | 北京理工大学 | A kind of multiple-pulse ring moon satellites formation initial method for considering error |
CN107992936A (en) * | 2017-12-26 | 2018-05-04 | 河北科技大学 | The hedging method of hedging method and flying object based on population |
CN115719035A (en) * | 2022-11-03 | 2023-02-28 | 哈尔滨工业大学 | Single-pair multi-satellite continuous fly-by observation trajectory optimization method and system considering sunlight constraint |
CN117262252A (en) * | 2023-09-21 | 2023-12-22 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | Spacecraft autonomous intersection and docking control method capable of realizing fuel optimization |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110032768B (en) * | 2019-03-15 | 2022-10-04 | 中国西安卫星测控中心 | Four-pulse orbit intersection optimization method using accurate dynamic model |
CN110986974B (en) * | 2019-10-28 | 2023-03-17 | 南京航空航天大学 | Complex dynamics environment-oriented multi-spacecraft task intelligent planning and control method |
CN111045447B (en) * | 2019-11-21 | 2023-08-29 | 浙江大学 | High-precision hypersonic aircraft track optimization multi-scale optimal control system |
CN113485396B (en) * | 2021-06-29 | 2023-11-10 | 北京控制工程研究所 | Spacecraft intersection butt joint final approximation segment relative orbit and attitude tracking control method |
CN113722821B (en) * | 2021-08-27 | 2024-01-23 | 北京理工大学 | Projection method for spacecraft intersection docking trajectory planning event constraint |
CN113886947B (en) * | 2021-09-13 | 2023-04-14 | 北京航空航天大学 | Aircraft static aeroelastic system output state quantity interval determination method based on iteration strategy |
CN114722711B (en) * | 2022-04-11 | 2022-12-06 | 中国科学院空天信息创新研究院 | Method and system for determining Lambert transfer orbit under photographic condition |
CN114861320B (en) * | 2022-05-19 | 2023-02-10 | 北京航天飞行控制中心 | Spacecraft attitude control thrust modeling and orbit determination resolving method |
CN114955009B (en) * | 2022-05-20 | 2023-03-28 | 北京航天飞行控制中心 | Selection method, system, medium and equipment of earth orbit intersection butt joint |
CN116449714B (en) * | 2023-04-20 | 2024-01-23 | 四川大学 | Multi-spacecraft pursuit game track control method |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050080521A1 (en) * | 2002-08-08 | 2005-04-14 | Horst-Dieter Fischer | Method for the fuel-optimized selection of a thruster configuration |
CN103455707A (en) * | 2013-07-22 | 2013-12-18 | 西北工业大学 | Finite-thrust driven spacecraft automatic-rendezvous trajectory planning method on basis of convex optimization technology |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103488830B (en) * | 2013-09-16 | 2016-08-31 | 北京航空航天大学 | The task simulation system that a kind of ground based on Cycler track moon comes and goes |
CN104229161B (en) * | 2014-08-22 | 2016-08-24 | 北京控制工程研究所 | A kind of spacecrafts rendezvous Trajectory Safety band based on control pulse determines method |
CN105005651B (en) * | 2015-07-09 | 2016-04-27 | 北京航空航天大学 | Optimization Design between the gradient cut section of spacecraft pulse intersection track |
-
2015
- 2015-07-09 CN CN201510400314.5A patent/CN105005651B/en not_active Expired - Fee Related
-
2016
- 2016-05-13 WO PCT/CN2016/082052 patent/WO2017005052A1/en active Application Filing
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050080521A1 (en) * | 2002-08-08 | 2005-04-14 | Horst-Dieter Fischer | Method for the fuel-optimized selection of a thruster configuration |
CN103455707A (en) * | 2013-07-22 | 2013-12-18 | 西北工业大学 | Finite-thrust driven spacecraft automatic-rendezvous trajectory planning method on basis of convex optimization technology |
Non-Patent Citations (2)
Title |
---|
伍升钢 等: "交会对接远程引导轨道控制状态规划与优化算法", 《飞行器测控学报》, vol. 31, no. 4, 31 August 2012 (2012-08-31), pages 90 - 94 * |
闫循良 等: "基于任务规划的有限推力燃料最优交会", 《固体火箭技术》, vol. 33, no. 6, 31 December 2010 (2010-12-31) * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017005052A1 (en) * | 2015-07-09 | 2017-01-12 | 北京航空航天大学 | Optimization and design method for gradient segmentation of intervals of spacecraft pulse rendezvous trajectory |
CN106407619A (en) * | 2016-11-16 | 2017-02-15 | 哈尔滨工业大学 | Linear-feedback global stabilization method for controlling limited spacecraft rendezvous control system |
CN106407619B (en) * | 2016-11-16 | 2019-08-02 | 哈尔滨工业大学 | A kind of linear feedback Stabilization method of Control constraints Spacecraft Rendezvous control system |
CN107526368A (en) * | 2017-09-12 | 2017-12-29 | 北京理工大学 | A kind of multiple-pulse ring moon satellites formation initial method for considering error |
CN107526368B (en) * | 2017-09-12 | 2020-02-11 | 北京理工大学 | Error-considered multi-pulse ring-moon satellite formation initialization method |
CN107992936A (en) * | 2017-12-26 | 2018-05-04 | 河北科技大学 | The hedging method of hedging method and flying object based on population |
CN115719035A (en) * | 2022-11-03 | 2023-02-28 | 哈尔滨工业大学 | Single-pair multi-satellite continuous fly-by observation trajectory optimization method and system considering sunlight constraint |
CN115719035B (en) * | 2022-11-03 | 2023-06-30 | 哈尔滨工业大学 | Single-to-multi-star continuous glancing observation track optimization method and system considering sunlight constraint |
CN117262252A (en) * | 2023-09-21 | 2023-12-22 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | Spacecraft autonomous intersection and docking control method capable of realizing fuel optimization |
CN117262252B (en) * | 2023-09-21 | 2024-06-11 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | Spacecraft autonomous intersection and docking control method capable of realizing fuel optimization |
Also Published As
Publication number | Publication date |
---|---|
CN105005651B (en) | 2016-04-27 |
WO2017005052A1 (en) | 2017-01-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105005651B (en) | Optimization Design between the gradient cut section of spacecraft pulse intersection track | |
CN108563243B (en) | Unmanned aerial vehicle track planning method based on improved RRT algorithm | |
CN104635741B (en) | Control System for Reusable Launch Vehicle reenters attitude control method | |
CN103995540A (en) | Method for rapidly generating finite time track of hypersonic aircraft | |
CN103869701A (en) | Attitude sequence resolving-based air vehicle novel real-time guide method | |
JP5822675B2 (en) | Multistage rocket guidance device, multistage rocket guidance program, multistage rocket guidance method, and multistage rocket guidance system | |
Feng et al. | Optimal control for a cooperative rendezvous between two spacecraft from determined orbits | |
CN103853047A (en) | Low thrust tracking guidance method based on state quantity feedback | |
Ma et al. | A unified trajectory optimization framework for lunar ascent | |
Fidan et al. | Longitudinal motion control of air-breathing hypersonic vehicles based on time-varying models | |
Ma et al. | Three-dimensional trajectory optimization for lunar ascent using Gauss pseudospectral method | |
Feng et al. | Optimization control for the far-distance rapid cooperative rendezvous of spacecraft with different masses | |
Yu et al. | Approximate analytical solutions for launch-vehicle ascent trajectory | |
CN106295218B (en) | A kind of quick determining energetic optimum intercepts the numerical optimization of predicted set-forward position | |
Cummings | Summary of Progress for the DoD HPCMP Hypersonic Vehicle Simulation Institute | |
Pontani et al. | Ascent Trajectories of Multistage Launch Vehicles: Numerical Optimization with Second‐Order Conditions Verification | |
Zhang et al. | Optimal satellite formation reconfiguration using pseudospectral method and co-evolutionary particle swarm optimization | |
Lu et al. | Closed-loop atmospheric ascent guidance based on finite element method | |
Wang et al. | A Fuel Optimization Method for the Pursuer in the Spacecraft Pursuit-Evasion Game | |
Saranathan | Algorithmic Advances to Increase the Fidelity of Conceptual Hypersonic Mission Design | |
Huang et al. | Cooperative Entry Trajectory Planning for Multiple Hypersonic Gliding Vehicles Using Convex Programing Method | |
Jiang et al. | Adaptive genetic algorithm in rendezvous orbit design | |
Sun et al. | Optimization method of multi-stages launch vehicle trajectory based on segmental Chebyshev pseudo spectral method | |
Gunnell et al. | Powered Descent and Landing of an Orbital-Class Rocket | |
Lingpu et al. | The optimal guidance law for rendezvous and docking of spacecraft |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20160427 |
|
CF01 | Termination of patent right due to non-payment of annual fee |