CN112693631B - Initial trajectory generation method and system in online sequential convex optimization of aircraft - Google Patents

Initial trajectory generation method and system in online sequential convex optimization of aircraft Download PDF

Info

Publication number
CN112693631B
CN112693631B CN202011352629.4A CN202011352629A CN112693631B CN 112693631 B CN112693631 B CN 112693631B CN 202011352629 A CN202011352629 A CN 202011352629A CN 112693631 B CN112693631 B CN 112693631B
Authority
CN
China
Prior art keywords
planning
initial
convex optimization
track
online
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202011352629.4A
Other languages
Chinese (zh)
Other versions
CN112693631A (en
Inventor
周祥
何睿智
谢磊
张洪波
汤国建
郑伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
National University of Defense Technology
Original Assignee
National University of Defense Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by National University of Defense Technology filed Critical National University of Defense Technology
Priority to CN202011352629.4A priority Critical patent/CN112693631B/en
Publication of CN112693631A publication Critical patent/CN112693631A/en
Application granted granted Critical
Publication of CN112693631B publication Critical patent/CN112693631B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B64AIRCRAFT; AVIATION; COSMONAUTICS
    • B64GCOSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
    • B64G1/00Cosmonautic vehicles
    • B64G1/22Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
    • B64G1/24Guiding or controlling apparatus, e.g. for attitude control
    • B64G1/242Orbits 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)
  • Traffic Control Systems (AREA)
  • Feedback Control In General (AREA)

Abstract

The invention discloses an initial trajectory generation method and system in aircraft online sequence convex optimization, aiming at the problem of guidance of a first-level recovery landing segment of a carrier rocket, firstly, generating an initial trajectory of a first planning period based on the number of fixed discrete points; secondly, determining a time range of online sequence convex optimization of the current planning cycle, including initial time and terminal time, and discretizing the time range; if the previous planning period can obtain an effective reference track, performing linear interpolation on the reference track; otherwise, carrying out linear interpolation on the reference track of the last but one planning cycle so as to obtain the state quantity of each discrete point in the time range; and finally, forming an initial track for online sequential convex optimization in the current planning period based on the state quantities at the discrete points. The method can effectively solve the problem of low convergence speed of the online sequence convex optimization algorithm caused by large initial track error.

Description

Initial trajectory generation method and system in online sequential convex optimization of aircraft
Technical Field
The invention relates to the technical field of aircraft guidance control, in particular to an initial trajectory generation method and system in online sequence convex optimization.
Background
The carrier rocket sublevel recovery technology is one of key technologies for realizing the reusability of an engine and reducing the space entering cost. The vertical landing section is the last stage of the recovery process, and the vertical landing section is controlled by a variable thrust engine to realize rocket fixed-point landing, which puts higher requirements on the robustness and reliability of the guidance method. The current successful guidance scheme is a guidance method based on reference track on-line generation, the method carries out on-line generation of the reference track according to the current state in each planning period, the adaptability to interference factors is wider, and the robustness of the guidance method is improved. The core of the guidance method is the online generation of the reference track, the current track planning method with online application potential mainly is a sequence convex method, the method makes full use of the rapidity of solving a convex optimization problem, and the reference track is planned online based on the current state under the conditions of a given initial track, the number of discrete points and the like. The selection of the initial track directly influences the convergence speed and precision of sequence convexity, thereby influencing the updating frequency and precision of the reference track, and finally reflecting the influence on the performance of the guidance method. The traditional initial trajectory is obtained on the basis of linear hypothesis, namely, all motion state quantities are considered to be linearly changed, and the initial trajectory obtained by the method does not meet the dynamic constraint and has larger deviation with a real trajectory, so that the convergence speed of the sequence convexity is lower, even is not converged.
Disclosure of Invention
The invention aims to solve the technical problem that aiming at the defects of the prior art, the invention provides the initial track generation method and the system in the online sequence convex optimization, solves the problem that the initial track generation method based on the linear hypothesis easily causes the slow convergence speed of the sequence convex and even does not converge, and realizes the rapid convergence of the online sequence convex.
In order to solve the technical problems, the technical scheme adopted by the invention is as follows: an initial trajectory generation method in online sequential convex optimization of an aircraft comprises the following steps:
in a first planning period, generating an initial track by adopting a linear hypothesis method;
starting from the second planning period, determining the initial track used for online sequence convex optimization in the current period based on the reference track successfully planned last time.
The invention provides an initial track generation method in online sequence convex optimization. The method effectively improves the precision of generating the initial track in the subsequent planning period, reduces the number of sequence convex iteration times, improves the convergence speed, and simultaneously ensures the guidance precision.
In a first planning period, a specific implementation process for generating an initial trajectory by using a linear hypothesis method includes:
set from the current state
Figure BDA0002801761990000021
To terminal state S f Is linearly changed with time from the current state
Figure BDA0002801761990000022
To terminal state S f Time range of
Figure BDA0002801761990000023
Each uniformly discrete time is
Figure BDA0002801761990000024
The superscript 1 represents the 1 st planning cycle, the subscript represents the discrete time sequence number, N is the number of discrete points, then the initial trajectory of the current planning cycle is
Figure BDA0002801761990000025
Let the sequence number k of the planning period be 1, wherein
Figure BDA0002801761990000026
The values of the x-direction component of the position at different discrete times,
Figure BDA0002801761990000027
for the values of the y-direction component of the position at different discrete times,
Figure BDA0002801761990000028
the values of the z-direction component of the position at different discrete times,
Figure BDA0002801761990000029
the values of the velocity x-direction component at different discrete times,
Figure BDA00028017619900000210
for the values of the velocity y-direction component at different discrete times,
Figure BDA00028017619900000211
the subscripts of the position component and the velocity component represent the serial numbers of the discrete points,
Figure BDA00028017619900000212
current state
Figure BDA00028017619900000213
Terminal state S f =[0 0 0 0 0 0] Τ . The initial trajectory generation method adopting the linear hypothesis has the advantages of simple engineering realization, small calculation amount and the like.
Figure BDA00028017619900000214
The value is 0, and the dynamic characteristic analysis of the rocket sublevel recovery vertical landing section is carried out to obtain the target
Figure BDA00028017619900000215
The value range is set to be 20-60 s, and the largest part of the flight time of the landing segment can be contained, so that the feasible reference track can be successfully planned through an online sequence convex optimization algorithm.
The implementation process of determining the initial track used for online sequence convex optimization in the current period based on the reference track successfully planned last time from the second planning period comprises the following steps: determining initial time and terminal time for performing online sequence convex optimization in the current planning period, and discretizing the time range; if the previous planning period obtains an effective reference track, performing linear interpolation on the reference track; otherwise, carrying out linear interpolation by using the reference track of the last but one planning cycle so as to obtain the state quantity of each discrete point in the time range; and the state quantities at the discrete points form an initial track for performing online sequential convex optimization in the current planning period. The problem of low convergence speed of online sequence convex optimization caused by large initial track error is further effectively solved.
The specific implementation process of acquiring the state quantity at the discrete point comprises the following steps:
1) based on initialized trajectories
Figure BDA0002801761990000031
Executing an online sequence convex optimization algorithm, and outputting a reference track R obtained in the kth planning cycle after the algorithm is converged (k) Laying a foundation for the initial track generation of the next planning period;
2) the current motion state is obtained by measurement
Figure BDA0002801761990000032
Will be provided with
Figure BDA0002801761990000033
As an initial time of the (k + 1) th planning cycle, and
Figure BDA0002801761990000034
T planning indicating planning cycle, motion state
Figure BDA0002801761990000035
As the initial state of the (k + 1) th planning period, providing an initial value for the (k + 1) th planning period to execute an online sequence convex optimization algorithm;
3) when in use
Figure BDA0002801761990000036
When so, ending; otherwise, turning to the step 4); wherein epsilon H Judging the tolerance for the landing condition; the step is used for judging whether the on-line track planning is terminated;
4) to be provided with
Figure BDA0002801761990000037
The time range for carrying out online sequence convex processing on the k +1 th planning period is equally divided into N-1 intervals to obtain corresponding moments at each discrete point in the time range, and the corresponding moments are recorded as
Figure BDA0002801761990000038
And is provided with
Figure BDA0002801761990000039
The step realizes the discretization of the continuous problem and is the basis for executing the online sequence convex optimization algorithm;
5) if the reference track R is obtained in the k-th track planning period (k) If the path is not empty, namely the path planning is successful, executing step 6); otherwise, turning to step 7); the step considers different initial track generation strategies under two conditions, and different steps are subsequently executed for different strategies;
6) reference trajectory obtained based on kth planning cycle
Figure BDA00028017619900000310
Linear interpolation is performed to obtain a reference trajectory R (k) At each discrete time
Figure BDA00028017619900000311
Is taken as
Figure BDA0002801761990000041
Will be provided with
Figure BDA0002801761990000042
As the (k + 1) th planPeriodically carrying out an initialization track of online sequence convex, adding 1 to the value of k, and turning to the step 1); the step obtains a higher-precision initialization track by performing linear interpolation on the reference track obtained in the last planning cycle, and improves the convergence speed of the online sequence convex optimization algorithm;
7) reference track obtained based on k-1 planning cycle
Figure BDA0002801761990000043
Performing linear interpolation to obtain R (k-1) At each discrete time
Figure BDA0002801761990000044
Is taken as
Figure BDA0002801761990000045
Will be provided with
Figure BDA0002801761990000046
And (4) taking the initial trajectory as the k +1 th planning period for online sequence convex, adding 1 to the value of k, and turning to the step 1).
ε H The allowable value of the landing height is expressed in a physical sense, strictly speaking, 0 should be taken, but considering that certain deviation necessarily exists in algorithm execution in actual engineering, 0.5-1 is recommended in the invention, and the unit is m.
The invention also provides an initial trajectory generation system in the online sequence convex optimization of the aircraft, which comprises computer equipment; the computer device is configured or programmed for performing the steps of the above-described method.
To facilitate the acquisition of the relevant data, the computer device communicates with a navigation device of the aircraft.
Compared with the prior art, the invention has the following beneficial effects: the method can effectively solve the problem of low convergence speed of the on-line sequence convex optimization algorithm caused by large initial track error, effectively improves the precision of generating the initial track in a subsequent planning period (not including a first planning period), and reduces the iteration times of the sequence convex optimization algorithm, thereby improving the convergence speed, and ensuring the landing guidance precision.
Detailed Description
Aiming at the problem of guidance of a first-stage recovery landing segment of a carrier rocket, the method firstly generates an initial track of a first planning period based on the number of fixed discrete points; secondly, determining the initial time and the terminal time of the online sequential convex optimization of the current planning cycle, and discretizing the time range; if the previous planning period can obtain an effective reference track, performing linear interpolation on the reference track; otherwise, carrying out linear interpolation by using the reference track of the last but one planning cycle so as to obtain the state quantity of each discrete point in the time range; and finally, forming an initial track for online sequential convex optimization in the current planning period based on the state quantities at the discrete points.
In the embodiment of the invention, the initial state of the landing segment is x by taking one-stage recovery of a certain type of carrier rocket as an example 0 =3528.5m,y 0 =-715.7m,z 0 =0,V x0 =-207.0m/s,V y0 =87.7m/s,V z0 The mass was 23000kg at-0.1 m/s. The planning period is T planning =2s,ε H 1 m. The sequence convex optimization algorithm can refer to: m Szmuk, B Acikmese, A W Berning. Successive conformation for fuel-optimal powered mapping with aerogenic drag and non-conditional constraints [ C].AIAA Guidance,Navigation,and Control Conference.San Diego,California USA,2016。
1) The number of discrete points adopted by the online sequence convex optimization method is set to be N, N is a positive integer, generally 10-100 can be taken, and the flight time range of the 1 st track planning is
Figure BDA0002801761990000051
Figure BDA0002801761990000052
The value of 0 is generally taken to be 0,
Figure BDA0002801761990000053
any real number between 20-60s may be taken. All position speed state quantities areDescribing in a land coordinate system, the terminal state of the landing segment is set to S f =[0 0 0 0 0 0] Τ The first three components represent the position in three directions and the last three components represent the velocity in three directions. After entering the landing segment, the initial state corresponding to the 1 st planning period measured by the navigation equipment is
Figure BDA0002801761990000054
Assuming a current state
Figure BDA0002801761990000055
To terminal state S f The position velocity of (a) varies linearly with time, each uniformly discrete time in the time range being
Figure BDA0002801761990000056
The superscript 1 represents the 1 st planning cycle, the subscript represents the discrete time sequence number, and the initial trajectory of the 1 st planning cycle is
Figure BDA0002801761990000057
Corresponding to a time variation range of
Figure BDA0002801761990000058
Let the planning period number k equal to 1, wherein
Figure BDA0002801761990000059
The values of the x-direction component of the position at different discrete times,
Figure BDA00028017619900000510
the values of the y-direction component of the position at different discrete instants,
Figure BDA00028017619900000511
the values of the z-direction component of the position at different discrete times,
Figure BDA00028017619900000512
the values of the velocity x-direction component at different discrete times,
Figure BDA00028017619900000513
for the values of the velocity y-direction component at different discrete times,
Figure BDA00028017619900000514
the value of the component in the z direction of the speed at different discrete moments is obtained, and the corresponding moment at each discrete point is
Figure BDA00028017619900000515
Wherein
Figure BDA00028017619900000516
2) Based on an initialization trajectory
Figure BDA0002801761990000061
Executing an online sequence convex optimization algorithm, and outputting a reference track R obtained in the kth planning cycle after the algorithm is converged (k)
3) The navigation equipment measures to obtain the current motion state as
Figure BDA0002801761990000062
Representing the position and velocity components in three directions, respectively. Will be provided with
Figure BDA0002801761990000063
As an initial time of the (k + 1) th planning cycle, and
Figure BDA0002801761990000064
state of motion
Figure BDA0002801761990000065
As the initial state of the (k + 1) th planning cycle;
4) when in use
Figure BDA0002801761990000066
When the rocket is in the ground state, the rocket is explained to be in the ground state, and the program is ended; otherwise, turning to the step 5;
5) to be provided with
Figure BDA0002801761990000067
The time range of the online sequence projection is carried out for the (k + 1) th planning period, the interval is equally divided into N-1 intervals, the corresponding time of each discrete point in the time range is obtained and recorded as
Figure BDA0002801761990000068
And is
Figure BDA0002801761990000069
6) If the reference track R is obtained in the k-th track planning period (k) If the path is not empty, namely the path planning is successful, executing the step 7, otherwise, executing the step 8;
7) reference trajectory obtained based on kth planning cycle
Figure BDA00028017619900000610
Linear interpolation is carried out, including linear interpolation of six state quantities of position and speed, so as to obtain six state quantities at each discrete moment
Figure BDA00028017619900000611
The value of (A) is recorded as
Figure BDA00028017619900000612
And taking the initial trajectory as an initialization trajectory for performing online sequence convex processing on the k +1 th planning period. Changing k to k +1, and turning to the step 2;
8) reference track obtained based on k-1 planning cycle
Figure BDA00028017619900000613
Linear interpolation is carried out, including linear interpolation of six state quantities of position and speed, so as to obtain six state quantities at each discrete moment
Figure BDA00028017619900000614
The value of (A) is recorded as
Figure BDA00028017619900000615
And taking the initial trajectory as an initial trajectory for performing online sequence convex processing on the k +1 th planning cycle. Changing k to k +1, and turning to the step 2;
the invention also provides an initial trajectory generation system in the online sequential convex optimization of the aircraft, which comprises computer equipment; the computer device is configured or programmed for performing the steps of the above-described embodiment method. To facilitate the acquisition of the relevant data, the computer device communicates with a navigation device of the aircraft. The navigation device is used for measuring the current motion state of the aircraft.
An experiment is given below to verify the beneficial effects of the present invention.
In order to verify that the method provided by the invention can improve the convergence rate of the sequence convex optimization algorithm, the following experiment compares the method of the invention with an initial trajectory generation method adopting linear hypothesis. Table 1 compares the results of the number of iterations required for the two methods applied in the sequence convex optimization algorithm. It can be seen that the method for generating the initial trajectory is beneficial to reducing the iteration times of the sequence convex optimization algorithm and improving the convergence rate. The landing guidance accuracy comparison under the two methods is shown in table 2, and it can be seen that the two methods can both meet the requirement of the landing guidance accuracy. The method can enable the sequence convex optimization algorithm to have higher convergence speed under the condition that the landing guidance accuracy is basically consistent.
TABLE 1 comparison of sequence convex iterations under different methods
Figure BDA0002801761990000071
Figure BDA0002801761990000081
Note that "X + 1" in table 1 above indicates that the linear hypothesis generation method is performed 1 time, and the sequence convex is iterated X times; "X" indicates that the method provided by the invention is adopted to generate an initial track, and the sequence is subjected to the convex iteration for X times.
TABLE 2 comparison of landing guidance accuracy under different methods
Linear hypothesis generation method Method of the invention
x 0.0340 0.0323
y -0.0625 0.0720
z 0.0118 0.0135
Vx 0.0058 0.0073
Vy -0.2070 -0.2392
Vz 0.0393 0.0450

Claims (6)

1. An initial trajectory generation method in online sequential convex optimization of an aircraft is characterized by comprising the following steps:
in a first planning period, generating an initial track by adopting a linear hypothesis method;
starting from the second planning period, determining an initial track used for online sequence convex optimization in the current period based on the reference track successfully planned last time; the implementation process comprises the following steps: determining a time range for performing online sequence convex optimization in the current planning period, wherein the time range comprises initial time and terminal time, and discretizing the time range; if the previous planning period obtains an effective reference track, performing linear interpolation on the reference track; otherwise, carrying out linear interpolation by using the reference track of the last but one planning cycle so as to obtain the state quantity of each discrete point in the time range; the state quantities at the discrete points form an initial track for performing online sequential convex optimization in the current planning period;
the specific implementation process of acquiring the state quantity at the discrete point comprises the following steps:
1) based on an initialization trajectory
Figure FDA0003584983570000011
Executing an online sequence convex optimization algorithm, and outputting a reference track R obtained in the kth planning cycle after the algorithm is converged (k)
2) The current motion state is obtained by measurement
Figure FDA0003584983570000012
Will be provided with
Figure FDA0003584983570000013
As an initial time of the (k + 1) th planning cycle, and
Figure FDA0003584983570000014
T planning indicating planning cycle, motion state
Figure FDA0003584983570000015
As the initial state of the (k + 1) th planning cycle;
3) when in use
Figure FDA0003584983570000016
When so, ending; otherwise, turning to the step 4); wherein epsilon H Judging the tolerance for the landing condition;
4) to be provided with
Figure FDA0003584983570000017
The time range of online sequence convex optimization is carried out for the (k + 1) th planning cycle, the time range is equally divided into N-1 intervals, and the corresponding time of each discrete point in the time range is obtained and recorded as
Figure FDA0003584983570000018
And is
Figure FDA0003584983570000019
5) If the reference track R is obtained in the k-th track planning period (k) If the path is not empty, namely the path planning is successful, executing step 6); otherwise, turning to step 7);
6) reference trajectory obtained based on kth planning cycle
Figure FDA00035849835700000110
Linear interpolation is performed to obtain a reference trajectory R (k) At each discrete time
Figure FDA0003584983570000021
Is taken as
Figure FDA0003584983570000022
Will be provided with
Figure FDA0003584983570000023
As the (k + 1) th planning cycleThe line online sequence convex optimization initialization track is carried out, the value of k is added with 1, and the step 1) is carried out;
7) reference track obtained based on k-1 planning cycle
Figure FDA0003584983570000024
Performing linear interpolation to obtain R (k-1) At each discrete time
Figure FDA0003584983570000025
Is taken as
Figure FDA0003584983570000026
Will be provided with
Figure FDA0003584983570000027
And (4) taking the initial trajectory as the k +1 th planning period for online sequence convex, adding 1 to the value of k, and turning to the step 1).
2. The method for generating the initial trajectory in the aircraft online sequential convex optimization according to claim 1, wherein in the 1 st planning period, the specific implementation process for generating the initial trajectory by using the linear hypothesis method comprises: set from the current state
Figure FDA0003584983570000028
To terminal state S f Is linearly changed with time from the current state
Figure FDA0003584983570000029
To terminal state S f Time range of
Figure FDA00035849835700000210
Each uniformly discrete time is
Figure FDA00035849835700000211
The 1 st programming cycle is indicated by the superscript 1, N is the number of discrete pointsTo that end, the initial trajectory of the current planning cycle is
Figure FDA00035849835700000212
Let the sequence number k of the planning period be 1, wherein
Figure FDA00035849835700000213
The values of the x-direction position component at different discrete times,
Figure FDA00035849835700000214
for the values of the y-direction position component at different discrete times,
Figure FDA00035849835700000215
for the values of the z-direction position component at different discrete times,
Figure FDA00035849835700000216
the values of the x-direction velocity component at different discrete times,
Figure FDA00035849835700000217
for the values of the y-direction velocity component at different discrete times,
Figure FDA00035849835700000218
the subscripts of the position component and the velocity component represent the serial numbers of discrete points,
Figure FDA00035849835700000219
current state
Figure FDA00035849835700000220
Terminal state
S f =[0 0 0 0 0 0] T
3. The aircraft online sequence of claim 2An initial trajectory generation method in convex optimization, characterized in that,
Figure FDA0003584983570000031
The value of the oxygen is 0, and the oxygen concentration is less than or equal to zero,
Figure FDA0003584983570000032
the value range is 20-60 s.
4. The method for generating initial trajectory in aircraft online sequential convex optimization according to claim 1, wherein ε H Taking 0.5-1, and taking the unit as m.
5. An initial trajectory generation system in an aircraft online sequential convex optimization is characterized by comprising computer equipment; the computer device is configured or programmed for carrying out the steps of the method according to one of claims 1 to 4.
6. The system of claim 5, wherein the computer device is in communication with a navigation device of the aircraft.
CN202011352629.4A 2020-11-27 2020-11-27 Initial trajectory generation method and system in online sequential convex optimization of aircraft Active CN112693631B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011352629.4A CN112693631B (en) 2020-11-27 2020-11-27 Initial trajectory generation method and system in online sequential convex optimization of aircraft

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011352629.4A CN112693631B (en) 2020-11-27 2020-11-27 Initial trajectory generation method and system in online sequential convex optimization of aircraft

Publications (2)

Publication Number Publication Date
CN112693631A CN112693631A (en) 2021-04-23
CN112693631B true CN112693631B (en) 2022-07-29

Family

ID=75506542

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011352629.4A Active CN112693631B (en) 2020-11-27 2020-11-27 Initial trajectory generation method and system in online sequential convex optimization of aircraft

Country Status (1)

Country Link
CN (1) CN112693631B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102890506A (en) * 2011-07-19 2013-01-23 北京理工大学 Small body approaching section guidance control method based on constraint programming
FR2991276A1 (en) * 2012-06-04 2013-12-06 Renault Sa DEVICE FOR CONTROLLING THE TRACK OF A VEHICLE
CN109470252A (en) * 2018-10-23 2019-03-15 哈尔滨工业大学 A kind of quick track optimizing method of VTOL Reusable Launch Vehicles based on convex optimization
CN110736470A (en) * 2019-11-06 2020-01-31 北京理工大学 method for hybrid search of continuous thrust tracks near small irregular celestial bodies
CN111580535A (en) * 2020-05-18 2020-08-25 中国人民解放军国防科技大学 Reentry trajectory three-dimensional profile planning method and system based on convex optimization
CN111931131A (en) * 2020-08-12 2020-11-13 北京航空航天大学 Online trajectory planning method and system for power landing segment of planetary probe

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6666410B2 (en) * 2001-10-05 2003-12-23 The Charles Stark Draper Laboratory, Inc. Load relief system for a launch vehicle

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102890506A (en) * 2011-07-19 2013-01-23 北京理工大学 Small body approaching section guidance control method based on constraint programming
FR2991276A1 (en) * 2012-06-04 2013-12-06 Renault Sa DEVICE FOR CONTROLLING THE TRACK OF A VEHICLE
CN109470252A (en) * 2018-10-23 2019-03-15 哈尔滨工业大学 A kind of quick track optimizing method of VTOL Reusable Launch Vehicles based on convex optimization
CN110736470A (en) * 2019-11-06 2020-01-31 北京理工大学 method for hybrid search of continuous thrust tracks near small irregular celestial bodies
CN111580535A (en) * 2020-05-18 2020-08-25 中国人民解放军国防科技大学 Reentry trajectory three-dimensional profile planning method and system based on convex optimization
CN111931131A (en) * 2020-08-12 2020-11-13 北京航空航天大学 Online trajectory planning method and system for power landing segment of planetary probe

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Rapid trajectory optimization for hypersonic entry using convex optimization and pseudospectral method;Wang, JB; Cui, NG and Wei, CZ;《AIRCRAFT ENGINEERING AND AEROSPACE TECHNOLOGY》;20190614;第91卷(第4期);669-679 *
变信赖域序列凸规划RLV再入轨迹在线重构;宗群,李智禹,叶林奇,田栢苓;《哈尔滨工业大学学报》;20200331;第52卷(第3期);147-155 *
基于凸优化的再入轨迹三维剖面规划方法;周祥,张洪波,何睿智,汤国建,包为民;《航空学报》;20201125;第41卷(第11期);623842-1至623842-17 *

Also Published As

Publication number Publication date
CN112693631A (en) 2021-04-23

Similar Documents

Publication Publication Date Title
CN109472418B (en) Maneuvering target state prediction optimization method based on Kalman filtering
CN109159151A (en) A kind of mechanical arm space tracking tracking dynamic compensation method and system
CN107024863A (en) A kind of UUV Trajectory Tracking Control methods for avoiding differential from exploding
CN107102617A (en) A kind of high-precision spatial elliptic curve Real-time Interpolation
CN112550770B (en) Rocket soft landing trajectory planning method based on convex optimization
CN104573190B (en) A kind of method for tracking target based on interactive multi-model
CN109375515B (en) Online identification method for dynamic characteristics of online trajectory planning of vertical take-off and landing rocket
CN110309909A (en) A kind of intelligent real-time predicting method of a wide range of maneuvering target track of high speed
CN106950999A (en) A kind of fitup Trajectory Tracking Control method of use Auto Disturbances Rejection Control Technique
CN103399488B (en) Multiple Model Control Method based on self study
CN114370793A (en) Rocket sublevel return and vertical landing guidance method
CN116184813A (en) Method, device, equipment and storage medium for controlling posture of boosting gliding rocket
CN112693631B (en) Initial trajectory generation method and system in online sequential convex optimization of aircraft
CN114660940A (en) Grinding constant force control method based on deep reinforcement learning PPO algorithm
CN112560343B (en) J2 perturbation Lambert problem solving method based on deep neural network and targeting algorithm
CN117289709A (en) High-ultrasonic-speed appearance-changing aircraft attitude control method based on deep reinforcement learning
CN114967472A (en) Unmanned aerial vehicle trajectory tracking state compensation depth certainty strategy gradient control method
CN116360258A (en) Hypersonic deformed aircraft anti-interference control method based on fixed time convergence
CN114800523B (en) Mechanical arm track correction method, system, computer and readable storage medium
CN113110559B (en) Optimal control method for small celestial body surface bouncing movement
CN114815878B (en) Hypersonic aircraft collaborative guidance method based on real-time optimization and deep learning
CN113670315B (en) Variation iterative Kalman filtering-based method for estimating attitude of Liqun heavy tail interference noise dynamic aircraft
Yuqi et al. Time-varying parameters estimation with adaptive neural network EKF for missile-dual control system
CN115576317A (en) Multi-preview-point path tracking control method and system based on neural network
CN112596537B (en) Model error compensation method, system and storage medium for online trajectory planning

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