CN116294837B - Control method for sub-missile head drop point based on perturbation guidance - Google Patents

Control method for sub-missile head drop point based on perturbation guidance Download PDF

Info

Publication number
CN116294837B
CN116294837B CN202310319729.4A CN202310319729A CN116294837B CN 116294837 B CN116294837 B CN 116294837B CN 202310319729 A CN202310319729 A CN 202310319729A CN 116294837 B CN116294837 B CN 116294837B
Authority
CN
China
Prior art keywords
sub
missile
point
head
missile head
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
CN202310319729.4A
Other languages
Chinese (zh)
Other versions
CN116294837A (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.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of 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 Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN202310319729.4A priority Critical patent/CN116294837B/en
Publication of CN116294837A publication Critical patent/CN116294837A/en
Application granted granted Critical
Publication of CN116294837B publication Critical patent/CN116294837B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/16Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F41WEAPONS
    • F41GWEAPON SIGHTS; AIMING
    • F41G3/00Aiming or laying means
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F41WEAPONS
    • F41GWEAPON SIGHTS; AIMING
    • F41G7/00Direction control systems for self-propelled missiles
    • F41G7/34Direction control systems for self-propelled missiles based on predetermined target position data
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F42AMMUNITION; BLASTING
    • F42BEXPLOSIVE CHARGES, e.g. FOR BLASTING, FIREWORKS, AMMUNITION
    • F42B15/00Self-propelled projectiles or missiles, e.g. rockets; Guided missiles
    • F42B15/01Arrangements thereon for guidance or control
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Abstract

A control method of a sub-missile head drop point based on perturbation guidance relates to a control method of an independent re-returning sub-missile head drop point. The invention solves the problem that the existing split missile head drop point control method ignores the influence of the earth rotation and the atmospheric resistance due to the adoption of two-body assumption in orbit mechanics, so that the result model error of solving the speed increment is larger. The method comprises the following steps: s1: binding longitude and latitude of an emission point and a target point; s2: adopting Newton iteration to calculate an emission azimuth angle and a first-stage maximum negative attack angle; s3: starting perturbation guidance to calculate the speed increment required by each sub missile head; s4: and storing and writing the obtained speed increment of each sub-missile head into a local file of a missile-borne computer, respectively applying the corresponding speed increment to the sub-missile head at the designated sub-guide time of each sub-missile head by a sub-guide cabin, and performing full-trajectory flight to finish the drop point control of the sub-missile head. The method has the characteristics of less iteration variable, high calculation speed and high result precision.

Description

Control method for sub-missile head drop point based on perturbation guidance
Technical Field
The invention relates to a control method for an independent re-dividing missile head drop point.
Background
In order to improve the defending capability and the number of hit targets of an intercontinental ballistic missile, a plurality of warheads are arranged in a main cabin so that one missile can hit a plurality of different targets. The launching units bound by the missile can only guarantee aiming at one target, so that the technical requirement for solving the speed increment required by maneuvering and rail transfer of the missile head is provided for realizing high-precision and large-scale striking depth of the missile head. The existing method adopts two body assumptions in orbit mechanics, ignores the influence of earth rotation and atmospheric resistance, solves the result model error obtained by speed increment, and is relatively large and difficult to be practically applied.
Disclosure of Invention
In order to solve the problem that the existing split missile head drop point control method ignores the influence of earth rotation and atmospheric resistance due to the adoption of two-body assumption in orbit mechanics, so that the result model error of solving the speed increment is larger; the provided control method for the missile head drop point based on perturbation guidance adopts a high-precision dynamics model under a launching system.
The control method of the sub-missile head drop point based on perturbation guidance comprises the following steps:
S1: binding longitude and latitude of an emission point and a target point;
s2: adopting Newton iteration to calculate an emission azimuth angle and a first-stage maximum negative attack angle;
s3: starting perturbation guidance to calculate the speed increment required by each sub missile head;
S4: storing and writing the obtained speed increment of each sub-missile head into a local file of a missile-borne computer, respectively applying the corresponding speed increment to the sub-missile head at the designated sub-guide time of each sub-missile head by a sub-guide cabin, and performing full-trajectory flight to finish the drop point control of the sub-missile head;
wherein, the longitudinal range L and the transverse range H of the missile in S3 are written into functional of shutdown time, shutdown point position and speed parameters (t k,vxk,vyk,vzk,xk,yk,zk)
Performing Taylor series expansion on the longitudinal range and the transverse range of the missile near the standard ballistic shutdown point parameters to obtain a first-order linear expression of the range with respect to the speed increment
Taking the deviation (delta x k,Δyk,Δzk) of the spatial position of the warhead relative to the nominal position at the diversion moment (t k) appointed by each sub-missile head as a disturbance term and taking the speed increment (delta v xk,Δvyk,Δvzk) as a control quantity; the longitudinal range is controlled by adopting Deltav xk, the transverse range is controlled by adopting Deltav zk, and the formula (2) is simplified into a form facing perturbation guidance iteration
Inverting the matrix equation to obtain the speed increment to be solved by one iteration of each sub-missile head
And merging the Deltav xk and Deltav zk of each bullet dividing head into the speed of the nominal shutdown point, and continuing to perform iterative calculation until the speed increment meeting the precision requirement of the drop point is obtained.
Wherein t k is the assigned guiding time of the missile head; v xk is the x-axis direction speed of the shutdown point; v yk is the speed of the shutdown point in the y-axis direction; v zk is the speed of the shutdown point in the z-axis direction; x k,yk,zk is the shutdown point location parameter.
Further, S1: binding the longitude and latitude of a launching point of the intercontinental ballistic missile, and binding the longitude and latitude of a target point for each sub-missile head.
Further, the closest target point is selected as the calculation basis of the firing data of the intercontinental ballistic missile. Therefore, the 1 st sub-guide bullet is separated from the sub-guide cabin without the need of providing a speed increment by the sub-guide cabin; the sub-pod provides speed increments from the 2 nd sub-missile head. Wherein, according to the separation sequence of the sub-guiding warhead and the sub-guiding cabin, the 1 st sub-guiding warhead, the 2 nd sub-guiding warhead, the 3 rd sub-guiding warhead are set, and the like.
Further, S2: given the primary maximum negative angle of attack and the initial value of the azimuth of emission, perturbation is applied to the angle of attack and the azimuth of emission, respectively.
Further, in S2, the launching azimuth angle and the first-stage maximum negative attack angle are calculated through the Jacobian matrix and the longitude and latitude difference of the drop point in an iterative mode until the drop point precision meets the requirement of the first warhead target point.
The invention adopts a high-precision dynamics model under a transmitting system, expands the first-order Talter of the longitudinal range and the transverse range based on a perturbation guidance theory, and converts the first-order Talter into a perturbation expression which only needs to iterate 2 speed increments. The method has the characteristics of less iteration variable, high calculation speed and high result precision.
The invention adopts Deltav xk to control the longitudinal range and adopts Deltav zk to control the transverse range; the iteration variable is few, and the speed increment required by the independent re-striking of different target points of the sub-missile head under the maneuvering and launching conditions of the intercontinental ballistic missile can be rapidly calculated only by iteratively calculating Deltav xk and Deltav zk; the number of iteration variables is effectively reduced, and the iterative calculation speed is improved. The two-body assumption is not introduced, the calculation of the velocity increment directly depends on a high-precision ballistic dynamics model of the launching system, so that the calculated velocity increment result has higher precision and is close to an actual value.
Drawings
FIG. 1 is a flow chart of an intercontinental guided missile flight control in accordance with example 1;
FIG. 2 is a ballistic diagram of three split warheads in example 1;
FIG. 3 is a graph showing the speed of three split warheads in example 1;
FIG. 4 is a graph of the ballistic tilt of three split heads of example 1;
FIG. 5 is a graph of the ballistic deflection angle of three split heads of example 1;
FIG. 6 is a graph of the x-axis velocity of the firing train of three split warheads in example 1;
FIG. 7 is a graph of the y-axis velocity of the firing train of three split warheads in example 1;
FIG. 8 is a graph of the z-axis velocity of the firing train of three split warheads in example 1.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
It should be noted that, without conflict, the embodiments of the present invention and features of the embodiments may be combined with each other.
The first embodiment is as follows: the control method of the sub-missile head drop point based on perturbation guidance in the embodiment comprises the following steps:
S1: binding longitude and latitude of an emission point and a target point;
s2: adopting Newton iteration to calculate an emission azimuth angle and a first-stage maximum negative attack angle;
s3: starting perturbation guidance to calculate the speed increment required by each sub missile head;
S4: storing and writing the obtained speed increment of each sub-missile head into a local file of a missile-borne computer, respectively applying the corresponding speed increment to the sub-missile head at the designated sub-guide time of each sub-missile head by a sub-guide cabin, and performing full-trajectory flight to finish the drop point control of the sub-missile head;
wherein, the longitudinal range L and the transverse range H of the missile in S3 are written into functional of shutdown time, shutdown point position and speed parameters (t k,vxk,vyk,vzk,xk,yk,zk)
Performing Taylor series expansion on the longitudinal range and the transverse range of the missile near the standard ballistic shutdown point parameters to obtain a first-order linear expression of the range with respect to the speed increment
Taking the deviation (delta x k,Δyk,Δzk) of the spatial position of the warhead relative to the nominal position at the diversion moment (t k) appointed by each sub-missile head as a disturbance term and taking the speed increment (delta v xk,Δvyk,Δvzk) as a control quantity; the longitudinal range is controlled by adopting Deltav xk, the transverse range is controlled by adopting Deltav zk, and the formula (2) is simplified into a form facing perturbation guidance iteration
Inverting the matrix equation to obtain the speed increment to be solved by one iteration of each sub-missile head
And then merging the Deltav xk and Deltav zk of each bullet dividing head into the speed of the nominal shutdown point for iterative calculation until the speed increment meeting the precision requirement of the drop point is obtained.
And guiding the sub-guide cabin to a space position and a motion state required by the release of the sub-missile head, immediately releasing the sub-guide missile head once the separation condition is met, and transferring the sub-guide cabin to a sub-guide program of the next sub-guide missile head. The sub-guide cabin generally needs to change the flight track once when one sub-guide warhead is released until the last warhead is released and the sub-guide process is finished.
The second embodiment is as follows: the present embodiment differs from the first embodiment in that: s1: binding the longitude and latitude of a launching point of the intercontinental ballistic missile, and binding the longitude and latitude of a target point for each sub-missile head. The other is the same as in the first embodiment.
And a third specific embodiment: the present embodiment differs from the first or second embodiment in that: according to the position characteristics of the target points, the target point closest to the target point is selected as the calculation basis of the firing data of the intercontinental ballistic missile. The other is the same as the first or second embodiment.
The specific embodiment IV is as follows: the present embodiment differs from one or more of the embodiments in that: s2: given the primary maximum negative angle of attack and the initial value of the azimuth of emission, perturbation is applied to the angle of attack and the azimuth of emission, respectively. The other is the same as in one of the first to third embodiments.
Fifth embodiment: the present embodiment differs from the first to fourth embodiments in that: and S2, iteratively calculating an emission azimuth angle and a first-stage maximum negative attack angle through the jacobian matrix and the longitude and latitude difference of the drop point until the drop point precision meets the requirement. The others are the same as in one to one fourth embodiments.
Example 1
1-8, The control method of the head drop point of the split missile based on perturbation guidance simulates launching an intercontinental ballistic missile with 3 split missiles:
step S1: binding longitude and latitude of an emission point and a target point;
binding the longitude and latitude of a launching point of the intercontinental ballistic missile, and binding the longitude and latitude of a target point for each sub-missile head.
Judging whether the target point is within the effective range according to the known performance of the ballistic missile. If the target point is within the range, according to the position characteristics of the target point, the target point closest to the target point is selected as the calculation basis of the firing data of the intercontinental ballistic missile. Therefore, the 1 st sub-guide bullet is separated from the sub-guide cabin without the need of providing a speed increment by the sub-guide cabin; in practice, the sub-pod provides speed increment from the 2 nd sub-missile head.
Step S2: adopting Newton iteration to calculate an emission azimuth angle and a first-stage maximum negative attack angle;
Taking the rotation of the earth and the perturbation influence of the non-spherical J2 term of the earth into consideration, constructing a barycenter dynamics equation of the intercontinental ballistic missile in a launching coordinate system as a ballistic calculation model; the first stage adopts a program negative attack angle for turning, the second stage adopts a given constant negative attack angle for turning, low trajectory and gravity turning during the boosting of the third stage engine. The first warhead is separated from the main cabin after the three-stage engine is shut down.
Two variable controls of the first-order program turning maximum negative attack angle alpha m and the emission azimuth angle A 0 are selected as emission data elements needing iteration. And the initial value of the attack angle is obtained according to an off-line obtained table and the range interpolation of the target point.
The initial value of the emission azimuth A 0 is solved by the following formula
Where O EON represents a vector with the centroid pointing to the north pole, O E F represents a vector with the centroid pointing to the emission point, and O E T represents a vector with the centroid pointing to the target point.
Step S3: starting perturbation guidance to calculate the speed increment required by each sub missile head;
the firing units enable the 1 st sub-missile head to be separated from the sub-missile cabin without providing speed increment by the sub-missile cabin; the split guide is actually started from the split guide 2 bullet. And guiding the sub-guide cabin to a space position and a motion state required by the release of the sub-missile head, immediately releasing the sub-guide missile head once the separation condition is met, and transferring the sub-guide cabin to a sub-guide program of the next sub-guide missile head. The sub-guide cabin generally needs to change the flight track once when one sub-guide warhead is released until the last warhead is released and the sub-guide process is finished.
Writing the longitudinal range L and the transverse range H of the missile in the S3 into functional functions of shutdown time, shutdown point position and speed parameters (t k,vxk,vyk,vzk,xk,yk,zk)
Performing Taylor series expansion on the longitudinal range and the transverse range of the missile near the standard ballistic shutdown point parameters to obtain a first-order linear expression of the range with respect to the speed increment
Taking the deviation (delta x k,Δyk,Δzk) of the spatial position of the warhead relative to the nominal position at the diversion moment (t k) appointed by each sub-missile head as a disturbance term and taking the speed increment (delta v xk,Δvyk,Δvzk) as a control quantity; the longitudinal range is controlled by adopting the Deltav xk, the transverse range is controlled by adopting the Deltav zk, and the formula 2 is simplified into a form facing perturbation guidance iteration
Inverting the matrix equation to obtain the speed increment to be solved by one iteration of each sub-missile head
And then merging the Deltav xk and Deltav zk of each bullet dividing head into the speed of the nominal shutdown point for iterative calculation until the speed increment meeting the precision requirement of the drop point is obtained.
S4: and storing and writing the obtained speed increment of each sub-missile head into a local file of a missile-borne computer, respectively applying the corresponding speed increment to the sub-missile head at the designated sub-guide time of each sub-missile head by a sub-guide cabin, and performing full-trajectory flight to finish the drop point control of the sub-missile head.

Claims (5)

1. A control method of a sub-missile head drop point based on perturbation guidance is characterized by comprising the following steps:
S1: binding longitude and latitude of an emission point and a target point;
s2: adopting Newton iteration to calculate an emission azimuth angle and a first-stage maximum negative attack angle;
s3: starting perturbation guidance to calculate the speed increment required by each sub missile head;
S4: storing and writing the obtained speed increment of each sub-missile head into a local file of a missile-borne computer, respectively applying the corresponding speed increment to the sub-missile head at the designated sub-guide time of each sub-missile head by a sub-guide cabin, and performing full-trajectory flight to finish the drop point control of the sub-missile head;
wherein, the longitudinal range L and the transverse range H of the missile in S3 are written into functional of shutdown time, shutdown point position and speed parameters (t k,vxk,vyk,vzk,xk,yk,zk)
Performing Taylor series expansion on the longitudinal range and the transverse range of the missile near the standard ballistic shutdown point parameters to obtain a first-order linear expression of the range with respect to the speed increment
Taking the deviation (delta x k,Δyk,Δzk) of the space position of the warhead relative to the standard trajectory shutdown point at the diversion time (t k) appointed by each sub-missile head as a disturbance item and taking the speed increment (delta v xk,Δvyk,Δvzk) as a control quantity; the longitudinal range is controlled by adopting Deltav xk, the transverse range is controlled by adopting Deltav zk, and the formula (2) is simplified into a form facing perturbation guidance iteration
Inverting the matrix equation to obtain the speed increment to be solved by one iteration of each sub-missile head
And then merging the Deltav xk and Deltav zk of each bullet dividing head into the speed of the standard ballistic shutdown point for iterative calculation until the speed increment meeting the precision requirement of the drop point is obtained.
2. The method for controlling a missile head drop point based on perturbation guidance according to claim 1, wherein the following steps are S1: binding the longitude and latitude of a launching point of the intercontinental ballistic missile, and binding the longitude and latitude of a target point for each sub-missile head.
3. The method for controlling the head drop point of the split missile based on perturbation guidance according to claim 2, wherein the closest target point is selected as the calculation basis of the firing data of the intercontinental ballistic missile according to the position characteristics of the target point.
4. A method of controlling a missile head landing point based on perturbation guidance according to claim 1 or 3, characterized by S2: given the primary maximum negative angle of attack and the initial value of the azimuth of emission, perturbation is applied to the maximum negative angle of attack and the azimuth of emission, respectively.
5. A method for controlling a missile head drop point based on perturbation guidance according to claim 1 or 3, wherein in S2, the emission azimuth and the first-order maximum negative attack angle are calculated iteratively by jacobian matrix and the drop point longitude and latitude difference until the drop point precision meets the requirement.
CN202310319729.4A 2023-03-28 2023-03-28 Control method for sub-missile head drop point based on perturbation guidance Active CN116294837B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310319729.4A CN116294837B (en) 2023-03-28 2023-03-28 Control method for sub-missile head drop point based on perturbation guidance

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310319729.4A CN116294837B (en) 2023-03-28 2023-03-28 Control method for sub-missile head drop point based on perturbation guidance

Publications (2)

Publication Number Publication Date
CN116294837A CN116294837A (en) 2023-06-23
CN116294837B true CN116294837B (en) 2024-05-03

Family

ID=86822196

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310319729.4A Active CN116294837B (en) 2023-03-28 2023-03-28 Control method for sub-missile head drop point based on perturbation guidance

Country Status (1)

Country Link
CN (1) CN116294837B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109798902A (en) * 2019-03-11 2019-05-24 北京星际荣耀空间科技有限公司 One kind being suitable for carrier rocket and enters the orbit modified interative guidance method
CN110044210A (en) * 2019-04-22 2019-07-23 中国人民解放军国防科技大学 Closed-circuit guidance on-line compensation method considering arbitrary-order earth non-spherical gravitational perturbation

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7566026B2 (en) * 2006-03-29 2009-07-28 Raytheon Company Onboard guidance method for ballistic missiles

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109798902A (en) * 2019-03-11 2019-05-24 北京星际荣耀空间科技有限公司 One kind being suitable for carrier rocket and enters the orbit modified interative guidance method
CN110044210A (en) * 2019-04-22 2019-07-23 中国人民解放军国防科技大学 Closed-circuit guidance on-line compensation method considering arbitrary-order earth non-spherical gravitational perturbation

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
一种空中发射弹道导弹摄动制导方法;郑晓龙等;火力与指挥控制;20130331;第38卷(第03期);第89-91页 *
多头分导的显式制导方程;范强;;航天控制;19861231(第04期);第1-6页 *
弹道导弹零射程线及其应用;郑伟等;飞行力学;20080831;第26卷(第04期);第73-76页 *
摄动制导诸元计算的一种模糊控制方法;马瑞萍等;航天控制;20110430;第29卷(第02期);第23-26、30页 *
摄动多头分导与单弹头制导的统一模式;常文森;;国防科技大学学报;19811231(第03期);第113-126页 *
椭圆轨道编队构型保持控制方法研究;宋申民等;黑龙江大学工程学报;20110531;第02卷(第02期);第109-118页 *

Also Published As

Publication number Publication date
CN116294837A (en) 2023-06-23

Similar Documents

Publication Publication Date Title
CN111473696B (en) Guidance rocket vertical attack guidance method based on drop point estimation
WO2008115216A2 (en) Apparatus, method and computer program product for weapon flyout modeling and target damage assesment
de Celis et al. Guidance and control for high dynamic rotating artillery rockets
CN108279005B (en) Guidance information reconstruction method under seeker data failure mode
CN109737812B (en) Air-to-ground guided weapon side attack method and device
CN104567545B (en) The method of guidance of RLV endoatmosphere powered phase
CN111591470A (en) Aircraft accurate soft landing closed-loop guidance method adapting to thrust adjustable mode
Ohlmeyer et al. Nonlinear integrated guidance-control laws for homing missiles
CN115406312B (en) Missile guidance control integrated method considering view angle and steering engine delay constraint
Park et al. Time-delay control for integrated missile guidance and control
Lukacs et al. Trajectory-shape-varying missile guidance for interception of ballistic missiles during the boost phase
CN116294837B (en) Control method for sub-missile head drop point based on perturbation guidance
Kamel et al. Simulation and modelling of flight missile dynamics and autopilot analysis
Proff et al. Study of impact point prediction methods for zero-effort-miss guidance: Application to a 155 mm spin-stabilized guided projectile
CN117471952A (en) Integrated control method for backstepping supercoiled sliding mode guidance of aircraft
CN116858037A (en) Method for correcting guidance trajectory in atmospheric extra-layer interception based on perturbation guidance
Wan et al. Suboptimal midcourse guidance with terminal-angle constraint for hypersonic target interception
Hahn et al. Predictive guidance of a projectile for hit-to-kill interception
Braun et al. Advances in inertial guidance technology for aerospace systems
Elbaioumy et al. Modelling and Simulation of Surface to Surface Missile General Platform
CN114153143A (en) Design method of guided missile nonsingular fixed time sliding mode guidance law
CN114610057A (en) Design method for maneuvering penetration strategy of high-Mach aircraft
Fairfax et al. Distributed and cooperative control for swarming munitions
Hou et al. Impact angle control guidance to intercept moving targets by virtual target technique
Slegers Predictive control of a munition using low-speed linear theory

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