CN110989665A - Remote guidance rocket projectile drop point prediction method based on experimental design and Kriging model - Google Patents
Remote guidance rocket projectile drop point prediction method based on experimental design and Kriging model Download PDFInfo
- Publication number
- CN110989665A CN110989665A CN201911240169.3A CN201911240169A CN110989665A CN 110989665 A CN110989665 A CN 110989665A CN 201911240169 A CN201911240169 A CN 201911240169A CN 110989665 A CN110989665 A CN 110989665A
- Authority
- CN
- China
- Prior art keywords
- point
- rocket projectile
- function
- remotely guided
- guided rocket
- 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.)
- Withdrawn
Links
- 238000000034 method Methods 0.000 title claims abstract description 63
- 238000013401 experimental design Methods 0.000 title claims description 10
- 238000012360 testing method Methods 0.000 claims abstract description 43
- 238000005314 correlation function Methods 0.000 claims abstract description 34
- 238000013461 design Methods 0.000 claims abstract description 18
- 238000012888 cubic function Methods 0.000 claims abstract description 11
- 238000012549 training Methods 0.000 claims description 25
- 230000008859 change Effects 0.000 claims description 19
- 230000000630 rising effect Effects 0.000 claims description 19
- 238000013178 mathematical model Methods 0.000 claims description 14
- 230000005484 gravity Effects 0.000 claims description 12
- 230000001133 acceleration Effects 0.000 claims description 10
- 230000005540 biological transmission Effects 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 6
- 238000004088 simulation Methods 0.000 claims description 6
- 238000013016 damping Methods 0.000 claims description 4
- 230000001419 dependent effect Effects 0.000 claims description 4
- 239000003814 drug Substances 0.000 claims description 4
- 230000003068 static effect Effects 0.000 claims description 4
- 238000013519 translation Methods 0.000 claims description 4
- 238000013400 design of experiment Methods 0.000 claims 1
- 238000012937 correction Methods 0.000 abstract description 4
- 238000001914 filtration Methods 0.000 description 4
- 230000010354 integration Effects 0.000 description 4
- 238000013213 extrapolation Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000012886 linear function Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000001174 ascending effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/10—Simultaneous control of position or course in three dimensions
- G05D1/107—Simultaneous control of position or course in three dimensions specially adapted for missiles
Landscapes
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention belongs to the technical field of guided rocket drop point prediction, and particularly relates to a remote guided rocket drop point prediction method based on test design and a Kriging model. For the drop point prediction corresponding to the arc-rising section, the selected correlation functions include a Spline function, a Matern function and a Cubic function; for the drop point prediction corresponding to the arc-dropping segment, the selected correlation functions include a Spline function, a Matern function, a Gauss function and a Cubic function. The invention provides an effective method for the trajectory correction control system of the remotely guided rocket projectile to predict the drop point in real time.
Description
Technical Field
The invention belongs to the technical field of guided rocket drop point prediction, and particularly relates to a remote guided rocket drop point prediction method based on test design and a Kriging model.
Background
The core technology for realizing the accurate striking of the remotely guided rocket projectile is trajectory correction, and the drop point prediction guidance is one of the main methods of trajectory correction. The method for rapidly and accurately predicting the projectile drop point is one of key technologies of drop point prediction guidance, and the accuracy and the real-time performance of the method can directly influence the ballistic trajectory correction effect. Common methods for predicting the drop point include numerical integration, linearization and filter extrapolation.
The numerical integration method is a method for obtaining a drop point through an iterative ballistic program on the basis of an established ballistic model, and the most typical ballistic model is a six-degree-of-freedom model. Theoretically, the drop point obtained by calculation through a six-degree-of-freedom ballistic model is the most accurate, but in the calculation process, a large amount of time is consumed for tedious iteration, and the requirement on hardware is high, so that a numerical integration method generally adopts ballistic models with other degrees of freedom, such as a four-degree-of-freedom model, a three-degree-of-freedom model and a two-degree-of-freedom model.
The linearization method is a method for obtaining a linear ballistic equation system by approximately linearizing a nonlinear outer ballistic model, solving an analytical solution of the linear ballistic equation system, and predicting a projectile drop point by using the analytical solution. The linearization method can quickly predict the projectile landing point, but is insufficient in accuracy.
The filtering extrapolation method is a method for extrapolating a rocket drop point through an established filtering trajectory model. Because the theoretical basis of the filtering extrapolation method is established on a linear system and a Gaussian noise environment, certain errors are inevitably generated when the method is used for a nonlinear system of a flight trajectory, and even filtering divergence is caused when the noise is non-Gaussian noise.
Compared with the short-range guided rocket projectile, the long-range guided rocket projectile has the advantages of longer range, larger ballistic height and longer flight time. Under the condition, if the numerical integration method is still adopted to predict the drop point of the remotely guided rocket projectile, the real-time requirement cannot be met at all, and the linear method cannot meet the precision requirement. Meanwhile, compared with an outer ballistic model of the short-range guided rocket projectile, the relatively accurate outer ballistic model of the long-range guided rocket projectile needs to comprehensively consider the influences of gravity eccentricity, surface curvature, the change of gravity acceleration along with height and latitude, Coriolis inertia force and other factors, and the influencing factors are determined by the latitude, the elevation and the shooting direction of a gun position, so that the influence of launching conditions needs to be considered when the drop point of the long-range guided rocket projectile is predicted. The method for predicting the drop point of the remotely guided rocket projectile under any launching condition has practical significance.
Disclosure of Invention
Technical problem to be solved
The technical problem to be solved by the invention is as follows: how to utilize the optimized Latin hypercube test design and the Kriging model to carry out high-precision rapid prediction on the drop point of the remotely guided rocket projectile.
(II) technical scheme
In order to solve the technical problem, the invention provides a remote guidance rocket projectile drop point prediction method based on experimental design and a Kriging model, which comprises the following steps:
step 1: establishing a remote rocket projectile motion equation set;
step 2: establishing a nonlinear mathematical model for predicting the landing point of the remotely guided rocket projectile;
and step 3: acquiring a training sample and a test sample;
and 4, step 4: selecting a proper correlation function of the Kriging model to train samples corresponding to the rising arc section and the falling arc section;
and 5: and loading the drop point prediction model meeting the requirements of precision and real-time performance into the missile-borne computer.
In the step 1, it is assumed that the remotely guided rocket projectile is predicted according to a landing point of the uncontrolled ballistic flight, and an uncontrolled ballistic mathematical model in a longitudinal plane of the remotely guided rocket projectile is used as a motion equation set of the remotely guided rocket projectile, and the motion equation set comprises the following steps:
in the formula: m is the flight quality of the remotely guided rocket projectile; t is time; vxAnd VyFor two speeds under the ground coordinate systemA degree component; omegaix1、ωiy1And ωiz1The three components of the rotating angular velocity vector of the projectile body relative to a translation coordinate system are in the projectile body coordinate system; j. the design is a squarex1、Jy1And Jz1Respectively the polar moment of inertia, the equatorial moment of inertia and the equatorial moment of inertia of the remotely guided rocket projectile; omegaeThe angular speed of the rotation of the earth axis; omegaex、ωeyAnd ωezIs omegaeThree components in a transmission coordinate system;
and theta is the pitch angle and ballistic inclination angle respectively, α is the angle of attack;
p is engine thrust; q is dynamic pressure; srefAnd LrefRespectively a reference area and a reference length; cxAnd CyRespectively a drag coefficient and a lift coefficient; x and y are two position components under a ground coordinate system; r0xAnd R0yTwo components of the emission point geocentric radial under an emission coordinate system are taken as the emission point; r is the model of the geocentric radial of any point on the trajectory; g'rThe component of the acceleration of the gravity along the direction of the earth center radial;
is the component of the acceleration of the gravity along the earth axis direction; l is0Is the geocentric latitude of the transmitting point; a. the0Is the transmit azimuth;
is the rate of change of the static moment coefficient with α;
for damping moment coefficient
Of wherein
mcIs the rate of change of the mass of the engine over time.
In the step 2, a nonlinear mathematical model for predicting the drop point of the remotely guided rocket projectile under the standard meteorological condition is established, and the dependent variable of the nonlinear mathematical model is the drop point YRAnd the independent variables are parameters of launching condition and flight state, including the latitude B of the gun position0Altitude H of gun location0Shoot to ATElevation H of target pointT、Vx、Vy、x、y、
The basic form is:
in the step 3, the training samples and the test samples are obtained by the method based on the optimized Latin hypercube test design and the remotely guided rocket projectile motion equation set.
Wherein, in the step 3,
as the acquisition method of the training sample and the test sample is the same, only the acquisition method of the training sample is given, and the method comprises the following steps:
s301: using optimized latin hypercube pairs including the latitude of the gun position B0Height H of gun position0Shoot to ATDistance XGElevation H of target pointTMedicine temperature TSCarrying out numerical test design on all the factors, wherein the value ranges of 6 factors are given according to the launching conditions and the range capability of the remotely guided rocket projectile;
s302: calculating a specific transmission condition B for each numerical test point in S3010、H0、AT、HTAnd TSAt a given range XGEach test point corresponds to a trajectory at a corresponding shooting angle under the standard meteorological condition;
s303: for each test point in S302, every other test point takes the engine shutdown time as a starting pointOutputting flight data of a long-distance rocket projectile at a certain time, including Vx、Vy、x、y、
ωz1An internal ballistic parameter;
s304: randomly generating a corresponding attack angle aiming at the flight data output each time in the S303, giving the change range of the attack angle according to the change range of the attack angle of the remotely guided rocket projectile at different stages in the actual flight process, and utilizing a formula
Recalculating the pitch angle in each flight data and comparing the pitch angle in each flight data in S303
Replacing with the latest pitch angle;
s305: and for the flight data and the corresponding launching conditions in the S304, recalculating the drop points corresponding to the flight data under the standard meteorological conditions by using a ballistic simulation mode in a longitudinal plane, and dividing the drop points into two groups according to the modes of the rising arc section samples and the falling arc section samples.
Wherein, in the step 4, the method is suitable for the correlation function of the Kriging model for predicting the drop points corresponding to the ballistic arc-rising section and the ballistic arc-falling section,
for the drop point prediction corresponding to the arc-rising section, the selected correlation function comprises: a Spline function, a Matern function and a Cubic function;
for the drop point prediction corresponding to the arc-dropping segment, the selected correlation function comprises: a Spline function, a Matern function, a Gauss function, and a Cubic function.
The value range of the non-negative integer q of the Matern function is 3-5.
The method can accurately and quickly predict the falling point of the remote guidance rocket projectile under any launching condition.
(III) advantageous effects
In order to solve the urgent need of fast and high-precision prediction of the drop point of the remotely guided rocket, the invention provides a method for predicting the drop point of the remotely guided rocket based on an optimized Latin hypercube test design and a Kriging model.
Drawings
FIG. 1 is a flow chart of remotely guided rocket projectile drop point prediction under standard meteorological conditions.
FIG. 2 is a graph of the prediction accuracy of a Kriging model with a correlation function of a Matern function to a drop point as a function of a non-negative integer q.
Detailed Description
In order to make the objects, contents, and advantages of the present invention clearer, the following detailed description of the embodiments of the present invention will be made in conjunction with the accompanying drawings and examples.
In order to solve the problems in the prior art, the invention provides a method for predicting the drop point of a remotely guided rocket projectile based on a test design and a Kriging model, which can accurately and quickly predict the drop point of the remotely guided rocket projectile under any launching condition, and as shown in figure 1, the method comprises the following steps:
step 1: establishing a remote rocket projectile motion equation set;
step 2: establishing a nonlinear mathematical model for predicting the landing point of the remotely guided rocket projectile;
and step 3: acquiring a training sample and a test sample;
and 4, step 4: selecting a proper correlation function of the Kriging model to train samples corresponding to the rising arc section and the falling arc section;
and 5: and loading the drop point prediction model meeting the requirements of precision and real-time performance into the missile-borne computer.
In the step 1, it is assumed that the remotely guided rocket projectile is predicted according to a landing point of the uncontrolled ballistic flight, and an uncontrolled ballistic mathematical model in a longitudinal plane of the remotely guided rocket projectile is used as a motion equation set of the remotely guided rocket projectile, and the motion equation set comprises the following steps:
in the formula: m is the flight quality of the remotely guided rocket projectile; t is time; vxAnd VyTwo velocity components under a ground coordinate system; omegaix1、ωiy1And ωiz1The three components of the rotating angular velocity vector of the projectile body relative to a translation coordinate system are in the projectile body coordinate system; j. the design is a squarex1、Jy1And Jz1Respectively the polar moment of inertia, the equatorial moment of inertia and the equatorial moment of inertia of the remotely guided rocket projectile; omegaeThe angular speed of the rotation of the earth axis; omegaex、ωeyAnd ωezIs omegaeThree components in a transmission coordinate system;
and theta is the pitch angle and ballistic inclination angle respectively, α is the angle of attack;
p is engine thrust; q is dynamic pressure; srefAnd LrefRespectively a reference area and a reference length; cxAnd CyRespectively a drag coefficient and a lift coefficient; x and y are two position components under a ground coordinate system; r0xAnd R0yTwo components of the emission point geocentric radial under an emission coordinate system are taken as the emission point; r is the model of the geocentric radial of any point on the trajectory; g'rThe component of the acceleration of the gravity along the direction of the earth center radial;
is the component of the acceleration of the gravity along the earth axis direction; l is0Is the geocentric latitude of the transmitting point; a. the0Is the transmit azimuth;
is the rate of change of the static moment coefficient with α;
for damping moment coefficient
Of wherein
mcIs the rate of change of the mass of the engine over time.
In the step 2, the nonlinear mathematical model for the remote guidance rocket projectile landing point prediction comprehensively considers launching condition factors and flight state parameters influencing the remote guidance rocket projectile landing point prediction precision under standard meteorological conditions, wherein the launching condition factors comprise a shot space latitude, a shot space elevation, a shot direction and a target point elevation, and the flight state parameters comprise the speed and the position in the x direction and the y direction and a pitch angle under a ground coordinate system;
a nonlinear mathematical model for predicting the drop point of the remotely guided rocket projectile under the standard meteorological condition is established, and the dependent variable is the drop point YRAnd the independent variables are parameters of launching condition and flight state, including the latitude B of the gun position0Altitude H of gun location0Shoot to ATElevation H of target pointT、Vx、Vy、x、y、
The basic form is:
in the step 3, the training samples and the test samples are obtained by the method based on the optimized Latin hypercube test design and the remotely guided rocket projectile motion equation set.
Wherein, in the step 3,
as the acquisition method of the training sample and the test sample is the same, only the acquisition method of the training sample is given, and the method comprises the following steps:
s301: using optimized latin hypercube pairs including the latitude of the gun position B0Height H of gun position0Shoot to ATRange of shootingXGElevation H of target pointTMedicine temperature TSCarrying out numerical test design on all the factors, wherein the value ranges of 6 factors are given according to the launching conditions and the range capability of the remotely guided rocket projectile;
s302: calculating a specific transmission condition B for each numerical test point in S3010、H0、AT、HTAnd TSAt a given range XGEach test point corresponds to a trajectory at a corresponding shooting angle under the standard meteorological condition;
s303: for each test point in S302, the flight data of the remote rocket projectile is output at regular intervals by taking the engine shutdown time as a starting point, wherein the flight data comprises Vx、Vy、x、y、
ωz1An internal ballistic parameter;
s304: randomly generating a corresponding attack angle aiming at the flight data output each time in the S303, giving the change range of the attack angle according to the change range of the attack angle of the remotely guided rocket projectile at different stages in the actual flight process, and utilizing a formula
Recalculating the pitch angle in each flight data and comparing the pitch angle in each flight data in S303
Replacing with the latest pitch angle;
s305: and for the flight data and the corresponding launching conditions in the S304, recalculating the drop points corresponding to the flight data under the standard meteorological conditions by using a ballistic simulation mode in a longitudinal plane, and dividing the drop points into two groups according to the modes of the rising arc section samples and the falling arc section samples.
Wherein, in the step 4, the method is suitable for the correlation function of the Kriging model for predicting the drop points corresponding to the ballistic arc-rising section and the ballistic arc-falling section,
for the drop point prediction corresponding to the arc-rising section, the selected correlation function comprises: a Spline function, a Matern function and a Cubic function;
for the drop point prediction corresponding to the arc-dropping segment, the selected correlation function comprises: a Spline function, a Matern function, a Gauss function, and a Cubic function.
The value range of the non-negative integer q of the Matern function is 3-5.
Example 1
In the embodiment, a high-precision rapid prediction method for a remotely guided rocket projectile landing point based on an optimized latin hypercube test design and a Kriging model is provided, as shown in fig. 1, the method comprises the following steps:
s1: and establishing a motion equation set of the remotely guided rocket projectile. The drop point of a remotely guided rocket projectile is primarily determined by motion in the longitudinal plane. On the basis of carrying out stress analysis on the remotely guided rocket projectile, a relatively accurate motion model in a longitudinal plane of the remotely guided rocket projectile is established by combining the conventional ballistic model and based on the principle of being as accurate as possible, and the model comprehensively considers the influences of gravity eccentricity, surface curvature, gravity acceleration along with the change of height and latitude, Coriolis inertia force and the like. Meanwhile, the invention mainly researches the remotely guided rocket projectile according to a drop point prediction method of the uncontrolled ballistic flight, so that only an uncontrolled ballistic mathematical model in a longitudinal plane is given:
in the formula: m is the flight quality of the remotely guided rocket projectile; t is time; vxAnd VyTwo velocity components under a ground coordinate system; omegaix1、ωiy1And ωiz1The three components of the rotating angular velocity vector of the projectile body relative to a translation coordinate system are in the projectile body coordinate system; j. the design is a squarex1、Jy1And Jz1Respectively the polar moment of inertia, the equatorial moment of inertia and the equatorial moment of inertia of the remotely guided rocket projectile; omegaeThe angular speed of the rotation of the earth axis; omegaex、ωeyAnd ωezIs omegaeThree components in a transmission coordinate system;
and theta is the pitch angle and ballistic inclination angle respectively, α is the angle of attack;
S2: and establishing a nonlinear mathematical model for predicting the landing point of the remotely guided rocket projectile. Analyzing the equation set of motion of the remotely guided rocket projectile established in the step S1, wherein under the standard meteorological condition, when the engine of the remotely guided rocket projectile is shut down, the drop point is mainly determined by the current flight state parameters and stress conditions, wherein the flight state is determined by the current flight state parameters and stress conditionsThe state parameters include velocity V in a ground coordinate systemxAnd VyPosition x and y under ground coordinate system, pitch angle
Pitch angle rate omegaz1. For the remote guidance rocket projectile, the pitch angle speed in the longitudinal plane is small, and the influence on the drop point is small, so that the influence of the pitch angle speed can be ignored when the drop point of the remote guidance rocket projectile is predicted. For the stress condition of the current flight state, besides the influence of the current flight state parameters, the stress condition is mainly influenced by the latitude of a gun position, the altitude of the gun position and the direction. Of course, the landing point of the remote rocket projectile is also related to the elevation of the target point. The analysis result is integrated to establish a nonlinear mathematical model for predicting the drop point of the remotely guided rocket projectile under the standard meteorological condition, and the dependent variable is the drop point YRAnd the independent variables are parameters of launching condition and flight state, including the latitude B of the gun position0Altitude H of gun location0Shoot to ATElevation H of target pointT、Vx、Vy、x、y、
The basic form is:
s3: training samples and test samples are obtained. Before establishing the non-linear functional relationship between the remotely guided rocket projectile landing point and each influencing factor in the step S2, a training sample and a test sample need to be obtained. Because the acquisition methods of the training sample and the test sample are the same, only the acquisition method of the training sample is given, and the method mainly comprises the following steps:
s301: using optimized Latin hypercube method to measure various factors (including the latitude B of the gun position)0Height H of gun position0Shoot to ATDistance XGElevation H of target pointTMedicine temperature TS) Carrying out numerical test design, wherein the value range of 6 factors is determined according to the remote guidance rocket projectileThe emission condition and the range capability are given;
s302: for each numerical test point in S301, a specific transmission condition (B) is calculated0、H0、AT、HTAnd TS) At a given range (X)G) Each test point corresponds to a trajectory at a corresponding shooting angle under the standard meteorological condition;
s303: for each test point in S302, the flight data (including V) of the remote rocket projectile is output at regular intervals by taking the engine shutdown time as a starting pointx、Vy、x、y、
ωz1Equal ballistic parameters);
s304: randomly generating a corresponding attack angle aiming at the flight data output each time in the S303, giving the change range of the attack angle according to the change range of the attack angle of the remotely guided rocket projectile at different stages in the actual flight process, and utilizing a formula
Recalculating the pitch angle in each flight data and comparing the pitch angle in each flight data in S303
Replacing with the latest pitch angle;
s305: and for the flight data and the corresponding launching conditions in the S304, recalculating the drop points corresponding to the flight data under the standard meteorological conditions by using a ballistic simulation program in a longitudinal plane, and dividing the drop points into two groups according to the mode of the rising arc section samples and the falling arc section samples.
S4: and selecting a correlation function to construct a Kriging model. And selecting a proper Kriging model correlation function for training the training sample in the step S3. Common correlation functions include Cubic function (Cubic), exponential function (Exp), gaussian function (Gauss), linear function (Lin), Spherical function (sphere), Spline function (Spline), and the like, as shown in table 1. In addition to the usual correlation functions described above, the Matern correlation function may be more efficient. When the positive parameter v takes a half integer, i.e. v q +1/2, where q is a non-negative integer, the expression of the Matern correlation function becomes very simple. In this case, the Matern correlation function becomes the product of an exponential function and a polynomial of order q:
in the formula: r (theta, x)(i),x(j)) Representing the training sample point x for the correlation function with the parameter theta(i)And x(j)Spatial correlation between them;
NDvis x(i)Dimension (d) of (a).
For the drop point prediction corresponding to the rising arc section, the selected correlation functions include a Spline function, a Matern function (the value of the non-negative integer q is 3-5) and a Cubic function, and for the drop point prediction corresponding to the falling arc section, the selected correlation functions include the Spline function, the Matern function (the value of the non-negative integer q is 3-5), a Gauss function and the Cubic function. When the training error meets the requirement or the iteration reaches a certain number of times, extracting a drop point prediction model, testing the prediction precision and the real-time property by using a test sample, outputting the drop point prediction model if the drop point prediction model meets the precision and the real-time property requirement, modifying the type of a correlation function or properly increasing the number of training samples if the drop point prediction model cannot meet the precision or the real-time property requirement, training a new sample, and repeating the process until the drop point prediction model meets the precision and the real-time property requirement at the same time.
S5: the non-linear functional relationship between the drop point and each influencing factor in step S4 is loaded into the missile-borne computer.
S6: during actual flight of the remotely guided rocket projectile, for any given B0、H0、AT、HT、Vx、Vy、x、y、
The corresponding drop point can be obtained through the non-linear function relationship in step S5.
Example 2
In this embodiment, a long-distance guided rocket projectile is taken as an example, and for an ascending arc section after the shutdown of an engine and a descending arc section before the final guidance, Kriging models with relevant functions of Cubic, Exp, Gauss, Lin, Matern, sphere and Spline are respectively used to verify the method. Aiming at all factors (gun position latitude, gun position elevation, shooting, range, target point elevation and chemical temperature), a test point is selected by adopting an optimized Latin hypercube test design, flight simulation data are output at regular intervals by taking the engine shutdown time as an initial point, and a training sample and a test sample are respectively flight simulation data generated on the basis of 1000 trajectories and 2000 trajectories. The training samples and the test samples are divided into two groups according to the mode of the arc rising section samples and the arc falling section samples.
For the rising arc section after the shutdown of the engine and the falling arc section before the end guidance, the change relation of the prediction precision of the Kriging model with the correlation function of the Matern function to the falling point along with the non-negative integer q is shown in FIG. 2, wherein the regression model part is a 4-order polynomial. When q belongs to [1,4], the drop point prediction precision corresponding to the rising arc section and the falling arc section is improved along with the increase of q, and when q belongs to [1,3], the falling point prediction precision corresponding to the rising arc section and the falling arc section is obviously improved along with the increase of q; when q belongs to [4,7], the drop point prediction precision corresponding to the arc rising section and the arc falling section is slightly reduced along with the increase of q. Taking the drop point prediction corresponding to the arc-rising section as an example, when q is 1, the Maximum Absolute Error (MAE) and the Root Mean Square Error (RMSE) of the drop point are 356.49m and 40.83m respectively; when q is 3, both fall to 140.15m and 13.44m, respectively; when q is 4, both further fall to 102.11m and 11.76 m; and when q is 7, they are 178.45m and 18.51m, respectively. Compared with the drop point prediction precision corresponding to the arc rising section, the Kriging model with the correlation function being the Matern function has higher prediction precision for the drop point corresponding to the arc falling section.
When the regression model is a 4 th order polynomial, the influence of the 7 correlation functions on the accuracy of the drop point prediction is shown in table 2, where Matern represents the optimal Matern correlation function (i.e., q is 4). When the number of training samples changes, the optimal nonnegative integer q of the Matern correlation function also changes, and a large amount of numerical simulation shows that the value of q is more suitable for 3-5. For the rising arc section and the falling arc section, the related functions have larger influence on the prediction precision of the falling points, the related functions Cubic, Matern and Spline have higher prediction precision on the falling points corresponding to the rising arc section and the falling arc section, a good interpolation result can be obtained mainly because the 3 related functions have good smoothness, the related functions Exp, Gauss, Lin and Spherical only have higher prediction precision on the falling points corresponding to the falling arc section, and the prediction precision of the related functions Exp, Lin and Spherical on the falling points corresponding to the falling arc section is obviously lower than that of the related functions Cubic, Matern and Spline on the falling points corresponding to the falling arc section. For the rising arc section and the falling arc section, the prediction accuracy of the 3 correlation function functions to the falling point is respectively Spline, Matern and Cubic from high to low, wherein the falling point MAE and RMSE corresponding to the correlation function Spline are respectively 83.17m and 10.46m, and 37.76m and 3.13 m. In summary, when the Kriging model is used to predict the drop points corresponding to the rising arc segment and the falling arc segment, the selected correlation functions include Spline, Matern, and Cubic, and if only the drop points corresponding to the falling arc segment are predicted, the Gauss function may also be selected.
TABLE 1
In the table:
TABLE 2
The above description is only a preferred embodiment of the present invention, and it should be noted that, for those skilled in the art, several modifications and variations can be made without departing from the technical principle of the present invention, and these modifications and variations should also be regarded as the protection scope of the present invention.
Claims (8)
1. A remote guidance rocket projectile drop point prediction method based on experimental design and a Kriging model is characterized by comprising the following steps:
step 1: establishing a remote rocket projectile motion equation set;
step 2: establishing a nonlinear mathematical model for predicting the landing point of the remotely guided rocket projectile;
and step 3: acquiring a training sample and a test sample;
and 4, step 4: selecting a proper correlation function of the Kriging model to train samples corresponding to the rising arc section and the falling arc section;
and 5: and loading the drop point prediction model meeting the requirements of precision and real-time performance into the missile-borne computer.
2. The method for predicting the landing point of the remotely guided rocket projectile based on the experimental design and the Kriging model as claimed in claim 1, wherein in the step 1, the remotely guided rocket projectile is assumed to be predicted according to the landing point of the uncontrolled ballistic flight, and the uncontrolled ballistic mathematical model in the longitudinal plane is used as a motion equation system of the remotely guided rocket projectile, as follows:
in the formula: m is the flight quality of the remotely guided rocket projectile; t is time; vxAnd VyTwo velocity components under a ground coordinate system; omegaix1、ωiy1And ωiz1The three components of the rotating angular velocity vector of the projectile body relative to a translation coordinate system are in the projectile body coordinate system; j. the design is a squarex1、Jy1And Jz1Respectively the polar moment of inertia, the equatorial moment of inertia and the equatorial moment of inertia of the remotely guided rocket projectile; omegaeThe angular speed of the rotation of the earth axis; omegaex、ωeyAnd ωezIs omegaeThree components in a transmission coordinate system;
and theta is the pitch angle and ballistic inclination angle respectively, α is the angle of attack;
p is engine thrust; q is dynamic pressure; srefAnd LrefRespectively a reference area and a reference length; cxAnd CyRespectively a drag coefficient and a lift coefficient; x and y are two position components under a ground coordinate system; r0xAnd R0yTwo components of the emission point geocentric radial under an emission coordinate system are taken as the emission point; r is the model of the geocentric radial of any point on the trajectory; gr' is the component of the acceleration of the earth's gravity along the radial direction of the earth's center;
is the component of the acceleration of the gravity along the earth axis direction; l is0Is the geocentric latitude of the transmitting point; a. the0Is the transmit azimuth;
is the rate of change of the static moment coefficient with α;
for damping moment coefficient
Of wherein
mcIs the rate of change of the mass of the engine over time.
3. The design of experiment and Kriging model based remotely guided rocket launch point prediction of claim 2The method is characterized in that in the step 2, a nonlinear mathematical model for predicting the drop point of the remotely guided rocket projectile under the standard meteorological condition is established, and the dependent variable is the drop point YRAnd the independent variables are parameters of launching condition and flight state, including the latitude B of the gun position0Altitude H of gun location0Shoot to ATElevation H of target pointT、Vx、Vy、x、y、
The basic form is:
4. the method for predicting the landing point of the remotely guided rocket projectile based on the experimental design and the Kriging model as claimed in claim 3, wherein in the step 3, the training samples and the test samples are obtained based on the optimized Latin hypercube experimental design and the remotely guided rocket projectile motion equation set.
5. The remotely guided rocket projectile landing point prediction method based on experimental design and Kriging model as claimed in claim 4, wherein in said step 3,
as the acquisition method of the training sample and the test sample is the same, only the acquisition method of the training sample is given, and the method comprises the following steps:
s301: using optimized latin hypercube pairs including the latitude of the gun position B0Height H of gun position0Shoot to ATDistance XGElevation H of target pointTMedicine temperature TSCarrying out numerical test design on all the factors, wherein the value ranges of 6 factors are given according to the launching conditions and the range capability of the remotely guided rocket projectile;
s302: calculating a specific transmission condition B for each numerical test point in S3010、H0、AT、HTAnd TSAt a given range XGEach test point corresponds to a trajectory at a corresponding shooting angle under the standard meteorological condition;
s303: for each test point in S302, the flight data of the remote rocket projectile is output at regular intervals by taking the engine shutdown time as a starting point, wherein the flight data comprises Vx、Vy、x、y、
ωz1An internal ballistic parameter;
s304: randomly generating a corresponding attack angle aiming at the flight data output each time in the S303, giving the change range of the attack angle according to the change range of the attack angle of the remotely guided rocket projectile at different stages in the actual flight process, and utilizing a formula
Recalculating the pitch angle in each flight data and comparing the pitch angle in each flight data in S303
Replacing with the latest pitch angle;
s305: and for the flight data and the corresponding launching conditions in the S304, recalculating the drop points corresponding to the flight data under the standard meteorological conditions by using a ballistic simulation mode in a longitudinal plane, and dividing the drop points into two groups according to the modes of the rising arc section samples and the falling arc section samples.
6. The method for predicting the drop point of the remotely guided rocket projectile based on the experimental design and the Kriging model as claimed in claim 5, wherein in the step 4, the correlation function of the Kriging model is applied to the drop point prediction corresponding to the ballistic rising arc section and the ballistic falling arc section,
for the drop point prediction corresponding to the arc-rising section, the selected correlation function comprises: a Spline function, a Matern function and a Cubic function;
for the drop point prediction corresponding to the arc-dropping segment, the selected correlation function comprises: a Spline function, a Matern function, a Gauss function, and a Cubic function.
7. The method for predicting the drop point of the remotely guided rocket based on the experimental design and the Kriging model as claimed in claim 6, wherein the non-negative integer q of the Matern function has a value range of 3-5.
8. The trial design and Kriging model-based remotely guided rocket projectile landing point prediction method as claimed in claim 1, wherein the method can accurately and rapidly predict the landing point of the remotely guided rocket projectile under any launching condition.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911240169.3A CN110989665A (en) | 2019-12-06 | 2019-12-06 | Remote guidance rocket projectile drop point prediction method based on experimental design and Kriging model |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911240169.3A CN110989665A (en) | 2019-12-06 | 2019-12-06 | Remote guidance rocket projectile drop point prediction method based on experimental design and Kriging model |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN110989665A true CN110989665A (en) | 2020-04-10 |
Family
ID=70090582
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911240169.3A Withdrawn CN110989665A (en) | 2019-12-06 | 2019-12-06 | Remote guidance rocket projectile drop point prediction method based on experimental design and Kriging model |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110989665A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110991051A (en) * | 2019-12-06 | 2020-04-10 | 北京京航计算通讯研究所 | Remote guidance rocket projectile drop point prediction system based on experimental design and Kriging model |
| CN111737815A (en) * | 2020-05-31 | 2020-10-02 | 南京航空航天大学 | Flight system drop point error correction method based on dimensionless form |
Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5419116A (en) * | 1993-07-02 | 1995-05-30 | The United States Of America As Represented By The Secretary Of The Navy | Miniscale ballistic motor testing method for rocket propellants |
| JP2008275447A (en) * | 2007-04-27 | 2008-11-13 | Mitsubishi Electric Corp | Rocket tracking radar system |
| JP2012106515A (en) * | 2010-11-15 | 2012-06-07 | Mitsubishi Space Software Kk | Rocket guidance computing device, rocket guidance system, rocket, rocket guidance computing program, and rocket guidance computing method of rocket guiding computing device |
| CN103307938A (en) * | 2013-04-23 | 2013-09-18 | 北京电子工程总体研究所 | Method for acquiring aerodynamic parameters of spinning projectile |
| CN107423556A (en) * | 2017-06-09 | 2017-12-01 | 南京理工大学 | A kind of computational methods of the Long Range Rocket Gun launch data based on RBF agent model |
| US20170350914A1 (en) * | 2016-06-02 | 2017-12-07 | Kenneth L. Oehler, III | System and method for long-range ballistic calculation |
| CN107679655A (en) * | 2017-09-15 | 2018-02-09 | 中国人民解放军63816部队 | A kind of space launch rocket drop point forecasting system |
| CN109101717A (en) * | 2018-08-07 | 2018-12-28 | 重庆大学 | Solid propellant rocket Reliability Prediction Method based on reality with the study of fuzzy data depth integration |
| CN110108162A (en) * | 2019-06-18 | 2019-08-09 | 北京电子工程总体研究所 | A kind of drop point that motion platform long distance is thrown automatically amendment Guidance and control method |
| CN110991051A (en) * | 2019-12-06 | 2020-04-10 | 北京京航计算通讯研究所 | Remote guidance rocket projectile drop point prediction system based on experimental design and Kriging model |
-
2019
- 2019-12-06 CN CN201911240169.3A patent/CN110989665A/en not_active Withdrawn
Patent Citations (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5419116A (en) * | 1993-07-02 | 1995-05-30 | The United States Of America As Represented By The Secretary Of The Navy | Miniscale ballistic motor testing method for rocket propellants |
| JP2008275447A (en) * | 2007-04-27 | 2008-11-13 | Mitsubishi Electric Corp | Rocket tracking radar system |
| JP2012106515A (en) * | 2010-11-15 | 2012-06-07 | Mitsubishi Space Software Kk | Rocket guidance computing device, rocket guidance system, rocket, rocket guidance computing program, and rocket guidance computing method of rocket guiding computing device |
| CN103307938A (en) * | 2013-04-23 | 2013-09-18 | 北京电子工程总体研究所 | Method for acquiring aerodynamic parameters of spinning projectile |
| US20170350914A1 (en) * | 2016-06-02 | 2017-12-07 | Kenneth L. Oehler, III | System and method for long-range ballistic calculation |
| CN107423556A (en) * | 2017-06-09 | 2017-12-01 | 南京理工大学 | A kind of computational methods of the Long Range Rocket Gun launch data based on RBF agent model |
| CN107679655A (en) * | 2017-09-15 | 2018-02-09 | 中国人民解放军63816部队 | A kind of space launch rocket drop point forecasting system |
| CN109101717A (en) * | 2018-08-07 | 2018-12-28 | 重庆大学 | Solid propellant rocket Reliability Prediction Method based on reality with the study of fuzzy data depth integration |
| CN110108162A (en) * | 2019-06-18 | 2019-08-09 | 北京电子工程总体研究所 | A kind of drop point that motion platform long distance is thrown automatically amendment Guidance and control method |
| CN110991051A (en) * | 2019-12-06 | 2020-04-10 | 北京京航计算通讯研究所 | Remote guidance rocket projectile drop point prediction system based on experimental design and Kriging model |
Non-Patent Citations (5)
| Title |
|---|
| XIN ZHANG等: "Impact point prediction guidance based on iterative process for dual-spin projectile with fixed canards", 《CHINESE JOURNAL OF AERONAUTICS》 * |
| 吴汉洲等: "基于摄动理论的弹丸落点预测精度仿真研究", 《计算机仿真》 * |
| 薄学纲等: "基于落点预测的火箭弹变系数末制导律", 《探测与控制学报》 * |
| 赵强 等: "基于代理模型的制导火箭炮发射诸元计算方法", 《北京航空航天大学学报》 * |
| 赵捍东等: "利用落点预测制导律提高火箭弹射击精度的方法研究", 《弹箭与制导学报》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110991051A (en) * | 2019-12-06 | 2020-04-10 | 北京京航计算通讯研究所 | Remote guidance rocket projectile drop point prediction system based on experimental design and Kriging model |
| CN110991051B (en) * | 2019-12-06 | 2023-03-10 | 北京京航计算通讯研究所 | Remote guidance rocket projectile drop point prediction system based on experimental design and Kriging model |
| CN111737815A (en) * | 2020-05-31 | 2020-10-02 | 南京航空航天大学 | Flight system drop point error correction method based on dimensionless form |
| CN111737815B (en) * | 2020-05-31 | 2024-02-20 | 南京航空航天大学 | Flight system drop point error correction method based on dimensionless form |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110991051B (en) | Remote guidance rocket projectile drop point prediction system based on experimental design and Kriging model | |
| CN107423556B (en) | Remote rocket gun emission data calculation method based on radial basis function proxy model | |
| US20080206718A1 (en) | Apparatus, method and computer program product for weapon flyout modeling and target damage assessment | |
| CN105608251B (en) | The BNSobol methods of helicopter fire control system precision sensitivity analysis | |
| CN110989665A (en) | Remote guidance rocket projectile drop point prediction method based on experimental design and Kriging model | |
| CN109033493B (en) | Unscented Kalman filtering-based method for identifying pneumatic parameters of high-speed rotation bullet | |
| CN110874055B (en) | Prediction and control method for hypersonic aircraft separation process under action of two-phase flow field | |
| CN115081251B (en) | Virtual battlefield-based artillery striking effect evaluation method | |
| CN115248038B (en) | SINS/BDS combined navigation engineering algorithm under emission system | |
| CN111611648B (en) | Method for evaluating maximum flight distance capability of boosting-gliding aircraft | |
| CN114153143A (en) | Design method of guided missile nonsingular fixed time sliding mode guidance law | |
| RU2265233C1 (en) | Device for determination of coordinates | |
| Fresconi et al. | Aerodynamic characterizations of asymmetric and maneuvering 105mm, 120mm, and 155mm fin-stabilized projectiles derived from telemetry experiments | |
| Baranowski et al. | Comparison of explicit and implicit forms of the modified point mass trajectory model | |
| CN116611160A (en) | Online real-time characteristic parameter identification and trajectory prediction method for uncontrolled aircraft based on measured trajectory parameters | |
| CN110765669A (en) | Method for identifying zero lift resistance coefficient of active section of axisymmetric wingless and rudder-free missile | |
| Keen et al. | Trajectory simulations should match flight tests and other lessons learned in 30 years of store-separation analysis | |
| Zygmunt et al. | Theoretical and experimental research of supersonic missile ballistics | |
| CN114636982A (en) | Equivalent coupling scattering center reconstruction method for analyzing close proximity dynamic warhead | |
| Charubhun et al. | Development of automatic firing angle calculation for ground to ground MLRS | |
| Milinovic et al. | Experimental and simulation testing of flight spin stability for small caliber cannon projectile | |
| Aldoegre | Comparison between trajectory models for firing table application | |
| Trzun et al. | Statistical Analysis of Missile's Trajectory Deviation Due to Production Errors | |
| Wang et al. | Optimization of tactical aircraft weapon delivery using Tactics Templates | |
| Demir et al. | Prediction of projectile impact point using approximate statistical moments |
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 | ||
| WW01 | Invention patent application withdrawn after publication |
Application publication date: 20200410 |
|
| WW01 | Invention patent application withdrawn after publication |