CN111191359A - Method for constructing air-defense missile space killer area model and calculating shooting data based on small amount of actual data - Google Patents
Method for constructing air-defense missile space killer area model and calculating shooting data based on small amount of actual data Download PDFInfo
- Publication number
- CN111191359A CN111191359A CN201911351195.3A CN201911351195A CN111191359A CN 111191359 A CN111191359 A CN 111191359A CN 201911351195 A CN201911351195 A CN 201911351195A CN 111191359 A CN111191359 A CN 111191359A
- Authority
- CN
- China
- Prior art keywords
- area
- killing
- killing area
- interface
- space
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a method for constructing an air-defense missile space killer zone model and calculating shooting data elements based on a small amount of actual data, which mainly comprises the following steps of: 1) acquiring air defense missile data; 2) establishing a space killer area model based on the air defense missile data; 3) and resolving the space killing area model to obtain shooting data. According to the invention, a killing area model is constructed by adopting polynomial fitting and linear fitting methods and mapping matrix algorithms of curved edge polygons, planar polygons and spatial polygons according to an analytic geometric mathematical model of a killing area and a small amount of public actual data, and the defects of blindness and low efficiency in the construction of the traditional spatial killing area are overcome.
Description
Technical Field
The invention relates to the field of missile weapon system application, in particular to an air-defense missile space killer area model construction and shooting data resolving method based on a small amount of actual data.
Background
The killing zone (killing zone) refers to a space area (a space range for killing a typical target with a given probability) for killing a specified aerial target by the air defense missile weapon system, and is an important operational comprehensive performance index of the air defense missile weapon system. The space shape of the air defense missile weapon system killing area is complex and comprises indexes such as a far boundary, a near boundary, a high boundary, a low boundary, a maximum airway shortcut, a maximum airway angle, a maximum altitude angle and the like. The general outline of the killing area is represented by a data table and a graph, most of the general outlines do not provide a mathematical model of the boundary of the killing area, and the space killing area is usually represented by two plane graphs of a vertical plane killing area and a horizontal plane killing area in engineering. The vertical plane killing area is a graph obtained after the vertical plane cuts the space killing area, the horizontal plane killing area is obtained after the horizontal plane cuts the space killing area, obviously, the two plane killing areas are infinite, and in practical application, only typical situations such as a zero-route shortcut vertical plane killing area and a horizontal plane killing area with a typical height are generally analyzed. The key technology of modeling the killing area is generally divided into four basic method types of a mathematical modeling method, a simulation method, a fitting method and a table look-up method, or a method combining the basic types. The existing mathematical model is not derived from theory, but seven space killing area expressions, five vertical killing area expressions with zero route shortcut and four horizontal killing area expressions with a certain height are formed by adopting a numerical fitting method according to actual shooting and simulation experiment data, so that the change rule of the boundary of the killing area can be theoretically explained, the internal relation and essence between discrete experiment data can be well reflected, and the method is very valuable for researching the air-defense missile shooting theory. However, due to the simple and clear characteristic of a single-part expression of the existing mathematical model of the space killing area, the research on the complex relationship correspondingly generated by the coupling interface between the boundary surfaces is less, and a high-precision and totally-enclosed space killing area engineering model constructed based on the mathematical model and actual data is rare. The simulation calculation method needs to construct and solve a missile simulation differential trajectory equation set consisting of a mass equation, a kinetic equation, a kinematic equation, a guidance equation and a control equation, carry out an interception simulation experiment to obtain interception encounter point data of the air-defense missile, obtain a miss distance range or an equal overload envelope, and then correct to obtain a calculated theoretical kill area. The numerical fitting method has the problems that the function is difficult to determine and the segmentation range of the function is difficult to master; the form of the fitting function is various. For the one-dimensional fitting problem of data, nodes can be smoothly connected through segmentation processing, and fitting errors can be made to be small at will. But for multi-dimensional problems, it is very difficult to process if segmentation is performed, because a curved surface or a curve, or even a higher dimensional situation, will appear at the multi-dimensional problem segmentation. For the occurrence of high dimensional situations, fitting the segmentation is not well handled. The table lookup interpolation method has higher requirements on calculation and storage and is less used in engineering application. A large amount of data of missile launching area envelope can be effectively utilized through BP artificial neural network solution, and integral data fitting is realized; for the case of large data volume or complex input and output relation, a deep neural network architecture and a deep learning method under a large data environment are adopted, and a deep belief network is introduced, so that the fitting performance can be improved, and higher precision can be achieved. However, the uncertainty of local convergence of the method is still little applied in engineering, and fitting data mainly comes from analog simulation data and is not suitable for fitting of a small amount of actual data. A method for generating a space killing area in a plane killing area needs a spline function and a piecewise interpolation method to process data to obtain the boundary of a vertical killing area, the vertical killing area is rotated to be combined with a horizontal killing area which is translated, then the space killing area is obtained by correcting a near interface, a side interface, a low interface and the like, or a discrete cluster of vertical killing areas distributed on an equidistant airway shortcut is used for representing. The method is complicated, the change of the characteristic points is complicated, the difficulty of generating a smooth space killing area by fitting a cluster of vertical killing areas is high, the occupied data space is large, the real-time performance is weak, and the calculation of the firing data of any application route shortcut and the formation of the horizontal killing area are inconvenient. The traditional analysis method of the killing area needs a large amount of actual data, but in most cases, enough actual data cannot be obtained, particularly for some introduced or newer models. The method is based on a small amount of actual data, and adopts a method of combining a mathematical model and numerical fitting to construct an air-defense missile killing area Matalab model, so that the efficiency is improved, and the engineering operability is enhanced. Therefore, the research on modeling of the air defense missile killer area has important significance for air defense combat simulation.
Disclosure of Invention
The present invention is directed to solving the problems of the prior art.
The technical scheme adopted for achieving the purpose of the invention is that the method for constructing the air-defense missile space-killing area model and calculating shooting data base based on a small amount of actual data mainly comprises the following steps:
1) and acquiring air defense missile data.
Further, the air defense missile data comprise a high bound H of a killing areamaxAnd a low boundary H of the killing areaminAnd killing zone far boundary SsyNear boundary S of killing zonesjMaximum height angle epsilon of killing areamaxMaximum route angle q of killer areamaxMaximum airway shortcut PmaxTime t of encounter of air defense missile and targetD practiceAnd encounter point slope distance Ds actual。
2) And establishing a space killing area model based on the air defense missile data.
The method mainly comprises the following steps of:
2.1) establishing a killer area analytic geometric mathematical model based on the killer area boundary parameters.
2.2) inputting the air defense missile data into the analytic geometric mathematical model of the killing area to obtain the engineering model of the closed killing area.
The closed killing area engineering model comprises a killing area low-near interface, a killing area high-near interface, a killing area far interface, a killing area high interface, a killing area low interface, a killing area positive side interface, a killing area negative side interface, a killing area upper front positive interface, a killing area upper front negative interface, a killing area lower front positive interface and a killing area lower front negative interface.
a) The method mainly comprises the following steps of determining the DGG 'D' of the low and near interface of the killing area:
1) establishing a low-near-field spherical equation, namely:
wherein P is the air route shortcut. H is the target height. And q is the target route angle. HminIs the low boundary of the killing area. DsjminThe minimum near-bound slant distance is a killing area. SscjIs the near boundary of the killing area. EpsilonmaxThe maximum elevation angle of the killing area. q. q.smaxThe maximum route angle of the killing area. And S is a killing area interface. HJIs the near-boundary height of the killing area.
2) Correcting the initial low and near sphere by using the find function to obtain a matrix set (S)scj,P,H)DGG′D′。
3) Using surf (S)scj,P,H)DGG′D′The function forms the killing zone low-near interface DGG 'D'.
b) The main steps for determining the high-near interface BDD 'B' of the killing area are as follows:
1) establishing a high-near cone equation, namely:
in the formula, HmaxIs the high boundary of the killing area. SsjIs a near boundary of the killing area.
2) Correcting the initial high and near boundary conical surface by using the find function to obtain a matrix group (S)sj,P,H)BDD′B′。
3) Invoking surf (S)sj,P,H)BDD′B′The function forms the killing zone high-near interface BDD 'B'.
c) The main steps for determining the far boundary EHH 'E' of the kill zone are as follows:
1) establishing a far-bound surface equation, namely:
in the formula, a, b and c are curved surface equation coefficients. SsyIndicating a far boundary of the kill zone.
2) The data (H) of m groups of air defense missilesj,Ssyj) Input into p ═ polyfit (H, S)syR) function, and fitting out far boundary of killing area by polynomial fitting methodThe non-linear function expression of (2). r is the fitting order. r is less than m. p is a vector formed by polynomial coefficients of the function F. j is 1, …, m. (H)j,Ssyj) And representing the jth air defense missile data.
3) In the interpolation interval [ Hmin,Hmax]Uniformly taking n points and substituting into the far boundary of the killing areaIn the function, a bus (H) of the far boundary of the killing area is linearly fittedi,Ssyi). Bus bar (H)i,Ssyi) Around the OH axis of the coordinate system within + -qmaxRotating within the range to obtain an initial far-bound curved surface. i is 1, …, n. (H)i,Ssyi) And representing the air defense missile data corresponding to the ith point.
4) Correcting the initial far-bound surface by using a find function to obtain a matrix group(s)sy,P,H)EHH′E′。
5) Using surf (S)sy,P,H)EHH′E′The function forms a far interface EHH 'E' of the kill zone.
d) The main steps for determining the high interface HABB ' A ' H ' of the killing area are as follows:
1) establishing a high-boundary plane equation, namely:
2) correcting the initial high-boundary plane equation by using a find function to obtain a matrix set (S, P, H)HABB′A′H′。
3) Using surf (S, P, H)HABB′A′H′The function forms the killer high interface HABB ' A ' H '.
e) The main steps for determining the killing area low interface GFF 'G' are as follows:
1) establishing a low-boundary plane equation, namely:
2) correcting the initial high-boundary plane equation by using a find function to obtain a matrix set (S, P, H)GFF′G′。
3) Using surf (S, P, H)GFF′G′The function forms the killer region low interface GFF 'G'.
f) The main steps for determining the interface on the front side of the killing area are as follows:
1) curve of distance bound
in the formula, H (,) represents the coordinates of H (,) in the formula.
Determining a characteristic point E of a curved-edge polygon HAE, wherein the method comprises the following steps: in the interpolation interval [ Hmin,Hmax]Uniformly taking n points in the inner part of the sample,and substituted into the far boundary of the killing areaIn the function, a function is obtainedSubstituting n points into function in sequenceIn (1) obtainingThe k point corresponding to the time so as to obtain the coordinate of the E pointHkThe height corresponding to the k-th point.
2) Respectively mapping the curved-edge polygons HAE into (n-k +1) -dimensional matrixes SHAEA (n-k +1) -dimensional matrix PHAEAnd (n-k +1) -dimensional matrix HHAEObtaining a matrix set (S)HAE,PHAE,HHAE)。
3) Invoking surf (S)HAE,PHAE,HHAE) The function establishes the killing zone positive side interface.
g) The main steps for determining the negative side interface of the killing area are as follows:
1) curve of distance boundStraight lineAnd a straight lineAnd (3) stretching into a curved polygon H ' A ' E ', namely the equation:
2) polygonal with curved edgeH ' A ' E ' are respectively mapped into (n-k +1) -dimensional matrixes SH′A′E′A (n-k +1) -dimensional matrix PH′A′E′And (n-k +1) -dimensional matrix HH′A′E′Obtaining a matrix set (S)H′A′E′,PH′A′E′,HH′A′E′)。
3) Invoking surf (S)H′A′E′,PH′A′E′,HH′A′E′) The function establishes the kill zone negative side interface.
h) The main steps for determining the front positive interface on the killing area are as follows:
1) plane of maximum route angleThe high interface of the killing area, the positive side interface of the killing area and the high-near interface of the killing area are stretched into a plane polygon ABCD, namely:
2) mapping a planar polygonal ABCD surface to an n-dimensional matrix SABCDN-dimensional matrix PABCDAnd n-dimensional matrix HABCDObtaining a matrix set (S)ABCD,PABCD,HABCD)。
3) Invoking surf (S)ABCD,PABCD,HABCD) The function yields the front positive interface on the kill zone.
i) The main steps for determining the front negative interface on the killing area are as follows:
1) plane of maximum route angleThe high boundary of the killing area, the negative side boundary of the killing area and the high-near boundary of the killing area are stretched into a plane polygon A 'B' C 'D', namely:
2) mapping a planar polygon A 'B' C 'D' surface to an n-dimensional matrix SA′B′C′D′、nDimension matrix PA′B′C′D′And n-dimensional matrix HA′B′C′D′Obtaining a matrix set (S)A′B′C′D′,PA′B′C′D,HA′B′C′D′)。
3) Invoking surf (S)A′B′C′D′,PA′B′C′D′,HA′B′C′D′) The function yields the front negative interface on the kill zone.
j) The main steps for determining the front positive interface under the killing area are as follows:
1) plane of maximum route angleThe low boundary of the killing area, the positive boundary of the killing area and the low-near boundary of the killing area are expanded into a spatial polygon CEFGD, namely:
in the formula, the point F is the starting point of the far field.
2) Mapping spatial polygons CEFGD into a n x 3 n-dimensional matrix SCEFGDN x 3n dimensional matrix PCEFGDAnd n x 3n dimensional matrix HCEFGDObtaining a matrix set (S)CEFGD,PCEFGD,HCEFGD)。
3) Invoking surf (S)CEFGD,PCEFGD,HCEFGD) The function yields the front positive interface under the kill zone.
k) The main steps for determining the front negative interface under the killing area are as follows:
1) plane of maximum route angleThe low interface, the front side interface and the low-near interface are flared into a spatial polygon C ' E ' F ' G ' D ', namely:
in the formula, the F' point is the starting point of the far boundary.
2) Will spaceThe polygon C ' E ' F ' G ' D ' is mapped to a matrix S of dimensions n x 3nCE′F′G′D′N x 3n dimensional matrix PC′E′F′G′D′And n x 3n dimensional matrix HC′E′F′G′D′Obtaining a matrix set (S)C′E′F′G′D′,PC′E′F′G′D′,HC′E′F′G′D′)。
3) Invoking surf (S)C′E′F′G′D′,PC′E′F′G′D′,HC′E′F′G′D′) The function yields the front negative interface under the kill zone.
3) And resolving the space killing area model to obtain shooting data.
The main steps of solving the space killing area model are as follows:
3.1) establishing a space transmitting area model based on the space killing area model. And inputting the air defense missile data into a space launching area model to obtain space launching area parameters.
The main steps for obtaining the parameters of the space transmitting area are as follows:
3.1.1) establishing a space transmitting area model, namely:
Sf=Ss+Vt·tD(12)
in the formula, VtIs the target speed. t is tDThe time of the air defense missile encountering the target. DfTo encounter the point skew. SfIs an emitter interface.
Time t of air defense missile encountering with targetDAs follows:
3.1.2) time t of encountering air defense missile with targetD practiceAnd encounter point slope distance Ds actualInput into the formula (14) to obtain the time calculation coefficient b1Time calculation coefficient b2Time calculation coefficient b3And calculating the time t of the air-defense missile encountering the targetD。
3.1.3) time t of encountering air defense missile with targetDAnd inputting the parameters into a space transmitting area model to obtain space transmitting area parameters.
And 3.2) resolving the space killing area model to obtain a horizontal killing area and a transmitting area under the target height.
The steps of resolving the space killing area model to obtain a horizontal killing area and a transmitting area under the target height are as follows:
3.2.1) determining the height interval Δ H of the cutting horizontal killing zone of the space killing zone.
3.2.2) to interpolate the interval [ H ]min,Hmax]Within the range of any height H as a reference and in the height interval
3.2.3) the encounter point and emission point Data of the horizontal killer region, horizontal emitter region at height H are extracted using Malab's Data Cursor function.
And 3.3) resolving the space killing area model to obtain a vertical killing area and a transmitting area under the short-cut of the route where the target is located.
The steps of resolving the space killing area model to obtain a vertical killing area and a transmitting area under the route shortcut where the target is located are as follows:
3.3.1) determining the air route shortcut interval delta P of the space killing area launching area cutting vertical killing area launching area.
3.3.2) short-cut range of air route [ -P ]max,Pmax]Taking any route shortcut P in the range as a reference, and taking the route shortcut P as a referenceAnd internally intercepting a vertical killing area and a transmitting area on the short-cut of the navigation path.
3.3.3) extracting the encounter point and the emission point Data of the vertical killing area and the vertical emission area on the airway shortcut P by utilizing the Data Cursor function of Malab.
3.4) extracting the encountering points and the transmitting points Data in the horizontal killing area, the horizontal transmitting area, the vertical killing area and the vertical transmitting area by using a Data Cursor function to obtain the transmitting time tWhen the emission is leftAnd time t of target flying out of emitting areaTime of flightAnd completing shooting data element calculation. The firing data includes encounter point slope, encounter point height, shot time, time-to-flight, shot distance and shot distance.
The emission remains as follows:
where s denotes the s-coordinate of the target in the firing coordinate system O-SPH. DFYIs the far bound slant distance of the emitting area. DFJThe near-bound slant range of the emitting area. SfyFar bound the emission area. SfjThe emission region is near-bound. t is t0The time the target is at the current location.
In a shooting coordinate system O-SPH, the air defense missile system guidance radar or the launching device is taken as the origin of coordinates. The OS axis is in the horizontal plane passing through the O point and is directed to be parallel to and opposite to the horizontal projection line of the target route. The OH axis is vertical to the horizontal plane and points to the sky. The OP axis is determined according to the right hand rule.
The target departure launch time is as follows:
ttime of flight=tFJ-tFY(16)
Wherein the time t when the target reaches the near boundary of the emission regionFJAs follows:
in the formula, v is a target speed.
It is worth explaining that the construction method of the killing area model has strong data requirement, large calculation amount and low efficiency, so that the construction method of the killing area model overcomes the defects of blindness and low efficiency of the construction method of the killing area model based on a high-precision and totally-enclosed space killing area engineering model constructed on the basis of a mathematical model and actual data aiming at the complex relation correspondingly generated by the coupling interfaces between the boundary surfaces.
The technical effect of the present invention is undoubted. Compared with the prior art, the invention has the advantages that: 1) according to the analytic geometric mathematical model of the killing area and a small amount of public actual data, a killing area model is constructed by adopting a polynomial fitting method, a linear fitting method and a mapping matrix algorithm of a curved polygon, a planar polygon and a space polygon, so that the defects of blindness and low efficiency in the construction of the traditional space killing area are overcome. 2) According to the time and the actual data of the slant range, the space killing area is used for solving the space transmitting area, a small amount of actual data information and characteristics are fully utilized, and the efficiency and the precision of data fitting of the space transmitting area are improved. 3) The constructed space killing area contains the information of all plane killing areas, the horizontal killing area and the launching area can be easily solved, and compared with a simulation calculation method, a BP neural network and a deep learning method, the method has the advantages of higher efficiency and stronger engineering operability. 4) The Data Cursor technology is used for extracting encounter point and launch point Data to complete shooting Data solution, and compared with a traditional method for iteratively solving the midpoint by using a nonlinear equation, the method is high in searching efficiency, high in speed and high in solving precision.
Drawings
FIG. 1 shows a Matlab model construction step of an air defense missile space killer area;
FIG. 2 is a schematic view of an air defense missile system space killing zone;
FIG. 3 is a low-near interface, a high-near interface, a far interface and a high interface of a space killing zone;
FIG. 4 is a front side interface of the kill zone;
FIG. 5 is the upper front face of the kill zone;
FIG. 6 is the lower front face of the kill zone;
FIG. 7 is a process of shot data calculation;
fig. 8 is a horizontal transmitting region model of H ═ 14 km;
fig. 9 is a horizontal kill zone model with H ═ 14 km;
fig. 10 is a vertical kill zone model with P-0;
fig. 11 is a vertical emitter model with P-0;
FIG. 12 is time-slope data for a type of air defense missile;
FIG. 13 is a model of the firing zone of an anti-subsonic missile killer zone;
fig. 14 shows a horizontal emitter model with H ═ 15 km.
Detailed Description
The present invention is further illustrated by the following examples, but it should not be construed that the scope of the above-described subject matter is limited to the following examples. Various substitutions and alterations can be made without departing from the technical idea of the invention and the scope of the invention is covered by the present invention according to the common technical knowledge and the conventional means in the field.
Example 1:
referring to fig. 1 to 7, a method for constructing an air-defense missile space-kill zone model and calculating shooting data base based on a small amount of actual data mainly comprises the following steps:
1) and acquiring air defense missile data.
The air defense missile data comprises a high bound H of a killing areamaxAnd a low boundary H of the killing areaminAnd killing zone far boundary SsyNear boundary S of killing zonesjMaximum height angle epsilon of killing areamaxMaximum route angle q of killer areamaxMaximum airway shortcut PmaxTime t of encounter of air defense missile and targetD practiceAnd encounter point slope distance Ds actual。
2) And establishing a space killing area model based on the air defense missile data.
The method mainly comprises the following steps of:
2.1) establishing a killer area analytic geometric mathematical model based on the killer area boundary parameters.
2.2) inputting the air defense missile data into the analytic geometric mathematical model of the killing area to obtain the engineering model of the closed killing area.
The closed killing area engineering model comprises a killing area low-near interface, a killing area high-near interface, a killing area far interface, a killing area high interface, a killing area low interface, a killing area positive side interface, a killing area negative side interface, a killing area upper front positive interface, a killing area upper front negative interface, a killing area lower front positive interface and a killing area lower front negative interface.
Establishing a shooting parameter rectangular coordinate system O _ SPH (namely O _ XYZ) by taking the air-defense missile system guided radar or the transmitting device as a coordinate origin, wherein an OS axis is parallel to and opposite to a horizontal projection line of a target route in a horizontal plane passing through an O point; the OH axis is vertical to the horizontal plane and points to the sky; the OP axis is determined in a right-hand system. As shown in fig. 2, assuming that the target flies linearly, after obtaining the target height H, the airway shortcut P and the target airway angle (course angle) q, the HABB ' a ' H ' is a high interface of the killing area, typically a part of a horizontal sector; GFF 'G' is a killer zone low interface, typically a portion of a horizontal sector; EHH 'E' is the far interface of the killing area, which is generally a curved surface; BDD 'B' is a high-near interface of a killing area, and is typically a part of a conical surface, and the top of the conical surface is at a point 0; DGG 'D' is a low-near interface of a killing area, and is typically a part of a spherical surface, and the center of the spherical surface is at a point O; HAE and H ' A ' E ' are two side near interfaces of a killing area, and the front side surface respectively has an included angle of +/-q between two vertical planes passing through the O point and having zero short-cut with an airwaymaxIn the vertical plane, the shortcuts of the positive side surface and the negative side surface of the air route are plus or minus P respectivelymaxThe vertical plane of (a).
a) The method mainly comprises the following steps of determining the DGG 'D' of the low and near interface of the killing area:
1) establishing a low-near-field spherical equation, namely:
wherein P is the air route shortcut. H is the target height. And q is the target route angle. HminIs the low boundary of the killing area. DsjminThe minimum near-bound slant distance is a killing area. SscjIs the near boundary of the killing area. EpsilonmaxThe maximum elevation angle of the killing area. q. q.smaxThe maximum route angle of the killing area. And S is a killing area interface. DmaxThe far-range slant distance is a killing area.
2) Correcting the initial low and near sphere by using the find function to obtain a matrix set (S)scj,P,H)DGG′D′。
3) Using surf (S)scj,P,H)DGG′D′Function forming killer regionLow proximity interface DGG′D′。
b) The main steps for determining the high-near interface BDD 'B' of the killing area are as follows:
1) establishing a high-near cone equation, namely:
in the formula, HmaxIs the high boundary of the killing area. SsjIs a near boundary of the killing area. The initial interface equation is in parentheses, and the reserved part of the interface equation to be corrected is in parentheses (the interface equations are the same way below).
2) Correcting the initial high and near boundary conical surface by using the find function to obtain a matrix group (S)sj,P,H)BDD′B′。
3) Invoking surf (S)sj,P,H)BDD′B′The function forms the killing zone high-near interface BDD 'B'.
c) The main steps for determining the far boundary EHH 'E' of the kill zone are as follows:
1) establishing a far-bound surface equation, namely:
in the formula, a, b and c are curved surface equation coefficients. SsyIndicating a far boundary of the kill zone.
2) The data (H) of m groups of air defense missilesj,Ssyj) Input into p ═ polyfit (H, S)syR) function, and fitting out far boundary of killing area by polynomial fitting methodThe non-linear function expression of (2). r is the fitting order. r is less than m. p is a vector formed by polynomial coefficients of the function F. j is 1, …, m. r is 2.
3) In the interpolation interval [ Hmin,Hmax]Uniformly taking n points and substituting into the far boundary of the killing areaIn the function, a bus (H) of the far boundary of the killing area is linearly fittedi,Ssyi). Bus bar (H)i,Ssyi) Around the OH axis of the coordinate system within + -qmaxRotating within the range to obtain an initial far-bound curved surface. i is 1, …, n.
4) Correcting the initial far-bound surface by using a find function to obtain a matrix group (S)sy,P,H)EHH′E′。
5) Using surf (S)sy,P,H)EHH′E′The function forms a far interface EHH 'E' of the kill zone.
d) The main steps for determining the high interface HABB ' A ' H ' of the killing area are as follows:
1) establishing a high-boundary plane equation, namely:
2) correcting the initial high-boundary plane equation by using a find function to obtain a matrix set (S, P, H)HABB′A′H′
3) Using surf (S, P, H)HABB′A′H′The function forms the killer high interface HABB ' A ' H '.
e) The main steps for determining the killing area low interface GFF 'G' are as follows:
1) establishing a low-boundary plane equation, namely:
the upper and lower boundaries of the killing area are Z ═ Hmax、Z=HminIs limited by the maximum way shortcut, the maximum way angle, the maximum altitude angle, the far bound and the near bound.
2) Correcting the initial high-boundary plane equation by using a find function to obtain a matrix set (S, P, H)GFF′G′
3) Using surf (S, P, H)GFF′G′The function forms the killer region low interface GFF 'G'.
f) The main steps for determining the interface on the front side of the killing area are as follows:
1) curve of distance bound
determining a characteristic point E of a curved-edge polygon HAE, wherein the method comprises the following steps: the boundary curve function of the killing area is processed by using a linear interpolation method, an interpolation interval is divided into a plurality of small intervals, and a linear interpolation function is constructed in each small interval so as to achieve the purposes of properly reducing the length of the interpolation interval and improving the interpolation precision. In geometric sense, piecewise linear interpolation is to use polygonal line approximation to replace curveIn the interpolation interval [ Hmin,Hmax]Uniformly taking n points and substituting into the far boundary of the killing areaIn the function, a function is obtainedSubstituting n points into function in sequenceIn (1) obtainingThe k point corresponding to the time so as to obtain the coordinate of the E pointHkThe height corresponding to the k-th point.
2) Respectively mapping the curved-edge polygons HAE into (n-k +1) -dimensional matrixes SHAEA (n-k +1) -dimensional matrix PHAEAnd (n-k +1) -dimensional matrix HHAEObtaining a matrix set (S)HAE,PHAE,HHAE)。
3) Invoking surf (S)HAE,PHAE,HHAE) The function establishes the killing zone positive side interface.
The algorithm is as follows:
I) using round () function pair Ssyi(i ═ k, k +1, …, n) for quantization
qi=round((Ssyi-min(Ssyi))/((max(Ssyi)-min(Ssyi))/(k-1))+1;
II) Linear interpolation function line () pair SHAEAssigning values to the 1-qi columns of the ith row;
SHAE(i,1:qi)=linspace(min(x6),Ssyi,qi));
III) reacting SHAEThe remaining (n-k +1-qi) portion of the ith row of (a) is set to a non-integer Nan;
V)Pmaxbecomes a matrix HHAE。
The key to applying piecewise interpolation is the problem of proper selection of interpolation points. The principle of selecting interpolation nodes in theory is to select interpolation nodes near the interpolation points as much as possible, and the more interpolation nodes, the higher the interpolation precision. However, for the air-defense missile killing area, the selection of interpolation nodes and the determination of the number of the nodes need to consider the operational background of the air-defense missile killing area, and certain mathematical calculation processing is performed on the graph of the killing area by combining the shooting command requirements. After a plane killing area is approximately processed by using a piecewise linear interpolation method, the original irregular killing area is replaced by a polygon formed by linear piecewise functions with a certain form, so that a part of a killing airspace is lost.
In order to minimize the damage airspace loss of the damage with high probability, the approximate polygon should cover the original damage area as much as possible. Obviously, the more interpolation nodes, the more coverage of the approximation polygon and the original killing area, and the higher the approximation precision. However, for the killing area used by the air defense missile shooting command, increasing the interpolation node to improve the approximation precision sometimes does not have great significance. Therefore, determining the number and distribution of interpolation nodes is a core problem of the planar killer area processing. On the premise of taking mathematical calculation analysis as a tool, the actual needs of tactical application of the method are fully considered. Under the precision of meeting certain combat requirements, the workload is simplified as much as possible, and no meaningless work is done.
After the above treatment, the killer area patterns with different shapes are all approximated to regular closed polygons. Thus, the killing area of any one boundary plane is represented as a set of coordinates of points, i.e., a row of points consisting of the vertices of the polygon of the killing area.
g) The main steps for determining the negative side interface of the killing area are as follows:
1) curve of distance boundStraight lineAnd a straight lineAnd (3) stretching into a curved polygon H ' A ' E ', namely the equation:
2) respectively mapping the curved edge polygons H ' A ' E ' into (n-k +1) dimensional matrixes SH′A′E′A (n-k +1) -dimensional matrix PH′A′E′And (n-k +1) -dimensional matrix HH′A′E′Obtaining a matrix set (S)H′A′E′,PH′A′E′,HH′A′E′)。
3) Invoking surf (S)H′A′E′,PH′A′E′,HH′A′E′) The function establishes the kill zone negative side interface.
h) The main steps for determining the front positive interface on the killing area are as follows:
1) plane of maximum route angleThe high interface of the killing area, the positive side interface of the killing area and the high-near interface of the killing area are stretched into a plane polygon ABCD, namely:
2) mapping a planar polygonal ABCD surface to an n-dimensional matrix SABCDN-dimensional matrix PABCDAnd n-dimensional matrix HABCDObtaining a matrix set (S)ABCD,PABCD,HABCD)。
3) Invoking surf (S)ABCD,PABCD,HABCD) The function yields the front positive interface on the kill zone.
The algorithm is as follows:
I)SABCD1 st line A to B, SABCDN-dimensional linear interpolation of S coordinates of the n-th line C to D;
II)SABCDi-2 to (n-1) actions
SABCD(:,i)=linspace(SABCD(1,i),SABCD(1000,i),n);
III)PABCD1 st line A to B, PABCDThe nth behavior of C to D is an n-dimensional linear interpolation of the P coordinates;
IV)PABCDi-2 to (n-1) actions
PABCD(:,i)=linspace(PABCD(1,i),PABCD(1000,i),n);
V)HABCD1 st line A to B, HABCDThe nth behavior of C to D is an n-dimensional linear interpolation of the P coordinates;
VI)HABCDi-2 to (n-1) actions
HABCD(:,i)=linspace(HABCD(1,i),HABCD(1000,i),n)。
i) The main steps for determining the front negative interface on the killing area are as follows:
1) plane of maximum route angleThe high boundary of the killing area, the negative side boundary of the killing area and the high-near boundary of the killing area are stretched into a plane polygon A 'B' C 'D', namely:
2) mapping a planar polygon A 'B' C 'D' surface to an n-dimensional matrix SA′B′C′D′N-dimensional matrix PA′B′C′D′And n-dimensional matrix HA′B′C′D′Obtaining a matrix set (S)A′B′C′D′,PA′B′C′D′,HA′B′C′D′)。
3) Invoking surf (S)A′B′C′D′,PA′B′C′D′,HA′B′C′D′) The function yields the front negative interface on the kill zone.
j) The main steps for determining the front positive interface under the killing area are as follows:
1) plane of maximum route angleThe low boundary of the killing area, the positive boundary of the killing area and the low-near boundary of the killing area are expanded into a spatial polygon CEFGD, namely:
in the formula, the point F is the starting point of the far field.
2) Mapping spatial polygons CEFGD into a n x 3 n-dimensional matrix SCEFGDN x 3n dimensional matrix PCEFGDAnd n x 3n dimensional matrix HCEFGDObtaining a matrix set (S)CEFGD,PCEFGD,HCEFGD)。
3) Invoking surf (S)CEFGD,PCEFGD,HCEFGD) The function yields the front positive interface under the kill zone.
The find function and the surf function are stored in Matlab.
The algorithm is as follows:
I) obtaining matrix groups (S) according to the algorithm of h) respectivelyCEF,PCEF,HCEF)、(SCDF,PCDF,HCDF) And (S)DFG,PDFG,HDFG);
II) constructing a new set S of matricesCEFGD=[SCEF,SCDF,SDFG],
PCEFGD=[PCEF,PCDF,PDFG]And HCEFGD=[HCEF,HCDF,HDFG]。
k) The main steps for determining the front negative interface under the killing area are as follows:
1) plane of maximum route angleThe low interface, the front side interface and the low-near interface are flared into a spatial polygon C ' E ' F ' G ' D ', namely:
in the formula, the F' point is the starting point of the far boundary.
2) Mapping the spatial polygon C ' E ' F ' G ' D ' into a n-x 3 n-dimensional matrix SC′E′F′G′D′N x 3n dimensional matrix PC′E′F′G′D′And n x 3n dimensional matrix HC′E′F′G′D′Obtaining a matrix set (S)C′E′F′G′D′, PC′E′F′G′D′,HC′E′F′G′D′)。
3) Invoking surf (S)C′E′F′G′D′,PC′E′F′G′D′,HC′E′F′G′D′) Function to obtain the lower front negative bound of the killing areaAnd (5) kneading.
3) And resolving the space killing area model to obtain shooting data.
The main steps of solving the space killing area model are as follows:
3.1) establishing a space transmitting area model based on the space killing area model. And inputting the air defense missile data into a space launching area model to obtain space launching area parameters.
As shown in fig. 8, according to the target height, the horizontal killing area and the emitting area at the height are solved by the space killing area and the emitting area. And shearing the air defense missile space killing area with the closed surface and the launching area to obtain the horizontal killing area. The horizontal kill zone is the cross-section resulting from cutting the kill zone at a given height level. The horizontal killing area and the launching area can more visually observe the position relation among the target, the killing area and the launching area, and the shooting data can be more easily solved.
To make the kill zone model take effect better in the computer, and provide real-time support for shoot command decisions, [ H ] can be usedmin,Hmax]A horizontal kill zone of any height within the range. Because the horizontal killing area and the far boundary and the near boundary of the launching area of most air-defense missile weapons are circular arc-shaped, the horizontal killing area and the launching area models with different heights can be conveniently generated by the space killing area model and the launching area model.
The main steps for obtaining the parameters of the space transmitting area are as follows:
3.1.1) establishing a space transmitting area model, namely:
in the formula, VtIs the target speed. t is tDThe time of the air defense missile encountering the target. DfTo encounter the point skew. SfIs an emitter interface.
Time t of air defense missile encountering with targetDAs follows:
3.1.2) time t of encountering air defense missile with targetD practiceAnd encounter point slope distance Ds actualInput into the formula (14) to obtain the time calculation coefficient b1Time calculation coefficient b2Time calculation coefficient b3And calculating the time t of the air-defense missile encountering the targetD。
3.1.3) time t of encountering air defense missile with targetDAnd inputting the parameters into a space transmitting area model to obtain space transmitting area parameters.
And 3.2) resolving the space killing area model to obtain a horizontal killing area and a transmitting area under the target height.
As shown in fig. 10 and 11, according to the target route shortcut, the vertical killing area and the launching area under the route shortcut are solved by the space killing area and the launching area. And shearing the air defense missile space killing area and the launching area with the closed surfaces to obtain a vertical killing area and a launching area. The vertical kill zone and launch zone are cross sections obtained by cutting the space kill zone and launch zone from a given airway shortcut vertical plane. The vertical killing area and the launching area can more intuitively observe the position relation between the maneuvering target in the vertical plane and the killing area and the launching area, and the shooting data can be solved more easily.
For the killing area and transmitting area model to take effect in computer and providing real-time support for shooting command decision, the [ -P ] can be usedmax,Pmax]Vertical kill areas and launch areas of any route shortcut within the range. As the vertical killing area and the far boundary and the near boundary of the launching area of most air-defense missile weapons are all polygonal closed by the multi-section broken lines, the vertical killing area and the launching area models of different air routes and shortcuts can be conveniently generated by the space killing area and the launching area models.
The steps of resolving the space killing area model to obtain a horizontal killing area and a transmitting area under the target height are as follows:
3.2.1) determining the height interval Δ H of the cutting horizontal killing zone of the space killing zone.
3.2.2) to interpolate the interval [ H ]min,Hmax]Within the range of any height H as a reference and in the height interval
3.2.3) the encounter point and emission point Data of the horizontal killer region, horizontal emitter region at height H are extracted using Malab's Data Cursor function.
And 3.3) resolving the space killing area model to obtain a vertical killing area and a transmitting area under the short-cut of the route where the target is located.
The steps of resolving the space killing area model to obtain a vertical killing area and a transmitting area under the route shortcut where the target is located are as follows:
3.3.1) determining the air route shortcut interval delta P of the space killing area launching area cutting vertical killing area launching area.
3.3.2) short-cut range of air route [ -P ]max,Pmax]Taking any route shortcut P in the range as a reference, and taking the route shortcut P as a referenceAnd internally intercepting a vertical killing area and a transmitting area on the short-cut of the navigation path.
3.3.3) extracting the encounter point and the emission point Data of the vertical killing area and the vertical emission area on the airway shortcut P by utilizing the Data Cursor function of Malab.
3.4) extracting the encountering points and the transmitting points Data in the horizontal killing area, the horizontal transmitting area, the vertical killing area and the vertical transmitting area by using the Data Cursor function of Malab to obtain the transmitting residual time tWhen the emission is leftAnd time t of target flying out of emitting areaTime of flightAnd completing shooting data element calculation. The firing data includes encounter point slope, encounter point height, shot time, time-to-flight, shot distance and shot distance.
The emission remains as follows:
where s denotes the s-coordinate of the target in the firing coordinate system O-SPH. DFYIs the far bound slant distance of the emitting area. DFJThe near-bound slant range of the emitting area. SfyFar bound the emission area. SfjThe emission region is near-bound. t is t0The time the target is at the current location.
In a shooting coordinate system O-SPH, the air defense missile system guidance radar or the launching device is taken as the origin of coordinates. The OS axis is in the horizontal plane passing through the O point and is directed to be parallel to and opposite to the horizontal projection line of the target route. The OH axis is vertical to the horizontal plane and points to the sky. The OP axis is determined according to the right hand rule.
The target departure launch time is as follows:
ttime of flight=tFJ-tFY(16)
Wherein the time t when the target reaches the near boundary of the emission regionFJAs follows:
in the formula, v is a target speed.
Example 2:
a method for constructing an air-defense missile space killer zone model and calculating shooting data elements based on a small amount of actual data mainly comprises the following steps:
1) constructing a Matlab model of a space killing area: according to the analytic geometrical mathematical model of the killing area, inputting a small amount of public actual data, and constructing a closed killing area engineering model by adopting a polynomial fitting method, a linear fitting method and a mapping matrix algorithm of a curved polygon, a planar polygon and a spatial polygon according to the sequence from a low-near interface, a high-near interface, a far interface, a high-low interface, a positive-negative side interface from the top to the bottom.
2) Solving shooting data elements by using a killing area model: solving shooting data elements by using a killing area model, and solving a space launching area by using a space killing area according to time and slope distance actual data; according to the target height, solving a horizontal killing area and a transmitting area under the height by the space killing area; according to the route shortcut where the target is located, solving a vertical killing area and a launching area under the route shortcut by a space killing area; and (4) extracting encounter point and emission point Data in the horizontal killing area and the vertical killing area by using a Data Cursor technology to complete shooting Data resolution.
Example 3:
an experiment for verifying a method for constructing an air-defense missile space killer zone model and calculating shooting data based on a small amount of actual data mainly comprises the following steps:
1) acquiring actual data: the actual data of the parameters of the air defense missile killer zone of a certain type are shown in the table 1, and the actual data of time and slope distance are shown in the figure 12.
TABLE 1 practical data table of parameters of air defense missile killer zone
2) Killer and emitter model construction for typical targets
A certain missile weapon system intercepts a horizontal uniform-speed flying anti-ship missile target with V less than or equal to 300m/s, and the calculation results of a killing area and a launching area are shown in figure 13.
2.1) killing zone
a) Low proximity interface
Killing zone low near-bound DGG 'D' equation:
b) high and near interface
Killer zone high-near boundary BDD 'B' equation:
c) remote interface
Killing zone remote boundary EHH 'E' equation:
d) high low interface
① killer zone high interface HABB ' a ' H ' equation:
② killer zone low interface GFF 'G' equation:
e) positive and negative side interface
① front side interface
Curved-edge polygon HAE equation:
② negative side interface
Curved polygon H ' A ' E ' equation:
f) front positive and negative surfaces
① front upper surface
Trapezoidal ABCD equation:
② negative front
Trapezoidal A 'B' C 'D' equation:
g) lower front positive and negative surface
① lower front face
Spatial polygon CEFGD equation:
② lower front negative
The spatial polygon C ' E ' F ' G ' D ' equation:
2.2) emission area
21 sets of data (t) were extracted from table 1 available ExcelD practice,DActual fact of killing) Fitting t by a polyfit () functionDHeel DKilling and killingBy a polynomial of the functional relationship (b), b can be obtained1,b2,b3;
Substituting into the desired space killing area (S)Killing and killing,PKilling and killing,HKilling and killing) Obtaining parameters of the killing areaCorresponding tDFrom (12) the spatial emission region (S) can be derivedLaunching,PLaunching,HLaunching) The model is shown in fig. 13.
3) Data resolution for typical target shooting
3.1) shooting data element calculation of horizontal equal-height straight line flight target
Assuming that the position where the target M is precisely steered is (120, 20, 14), the target flies at 300M/s along the OP axis at the same level. At this time, the horizontal transmitting area H-14 km shown in FIG. 8 and the horizontal killing area shown in FIG. 9 are intercepted. The intersection points of the target and the far and near boundaries of the emitting area are A, B respectively, and the intersection points of the target and the far and near boundaries of the killing area are C, D respectively. Matlab can solve for the A, B, C, D coordinates and solve for the shot data.
a) When emission is left:
b) flying-out time:
ttime of flight=tFJ-tFY=266.50s
Table 2 shows the data of the shots at different heights when the shortcut of the flight path is 20 km.
TABLE 2 data of shots at different heights
H(km) | Sfy(km) | Sfj(km) | TWhen the emission is left(s) | TTime of flight(s) |
10 | 93.57 | 19.8 | 88.1 | 245.9 |
11 | 95.71 | 20.2 | 80.97 | 251.70 |
12 | 98.57 | 20.5 | 71.43 | 260.24 |
13 | 100.0 | 21.0 | 66.67 | 263.33 |
14 | 100.4 | 20.45 | 65.33 | 266.50 |
15 | 100.2 | 21.25 | 66 | 263.17 |
16 | 99.7 | 21.6 | 67.67 | 260.33 |
17 | 97.86 | 22.0 | 73.8 | 252.87 |
18 | 95.71 | 22.3 | 80.97 | 244.70 |
4) Data resolution of horizontal serpentine maneuvering flight target
Assuming that the position where the target M is precisely steered is (120, 0, 15), the target level maneuvers. The truncated H-15 km horizontal transmission area is shown in fig. 14. The intersection points of the target and the far and near boundaries of the emitting area are A, B respectively, and coordinates A (102.8, 3.311, 14.86) and B (13.86, -3.094,14.99) can be obtained by adopting Matlab's DataCursor technology. Although the target adopts horizontal maneuvering flight, the flight speed component of the target along the S axis is 300m/S, so the shooting data of the horizontal maneuvering target is solved in a transmitting area corresponding to a killing area of a high-straight-line flight target with the height of H being 15km horizontally. Then Sfy=102.8km,Sfj=13.86km, TWhen the emission is left=57.33s,TTime of flight=296.47s。
5) Serpentine maneuvering flight target data element calculation in vertical plane
Assuming that the flying velocity component of the target along the S axis is constant, the position of the target which is accurately guided is (120, 0, 10), and the target flies in a snake-shaped maneuvering plane in the vertical plane. The vertical emission region with cut P-0 is shown in fig. 11.
The intersection points of the target and the far and near boundaries of the emitting region are A, B respectively, and coordinates A (99.88, 0, 11.36) and B (7.391, 0, 14.44) can be obtained by adopting the Data Cursor technology of Matlab. Although the target adopts vertical maneuvering flight, the flight speed component of the target along the S axis is 300m/S, so the shooting data of the vertical maneuvering target is solved in a transmitting area corresponding to a killing area of a high-straight-line flight target with the H-0 km horizontal height. Then Sfy=99.88km, Sfj=7.39km,TWhen the emission is left=67.07s,TTime of flight=308.30s。
Claims (10)
1. A method for constructing an air-defense missile space-killing area model and calculating shooting data elements based on a small amount of actual data is characterized by mainly comprising the following steps of:
1) and acquiring the air defense missile data.
2) Establishing a space killer area model based on the air defense missile data;
3) and resolving the space killing area model to obtain shooting data.
2. The method for constructing model of air-defense missile space-kill zone and solving shooting data base on the basis of a small amount of actual data according to claim 1, wherein the air-defense missile data comprises a high bound H of the kill zonemaxAnd a low boundary H of the killing areaminAnd killing zone far boundary SsyNear boundary S of killing zonesjMaximum height angle epsilon of killing areamaxMaximum route angle q of killer areamaxMaximum airway shortcut PmaxTime t of encounter of air defense missile and targetD practiceAnd encounter point slope distance Ds actual。
3. The method for constructing the air-defense missile space-killing area model and shooting data element calculation based on the small amount of actual data according to claim 1, wherein the method for establishing the space-killing area model mainly comprises the following steps:
1) establishing a killer area analytic geometric mathematical model based on the killer area boundary parameters;
2) and inputting the air defense missile data into the analytic geometric mathematical model of the killing area to obtain the engineering model of the closed killing area.
4. The method for constructing the air-defense missile space-killer zone model and calculating shooting data elements based on the small amount of actual data according to claim 3, wherein the method comprises the following steps: the closed killing area engineering model comprises a killing area low-near interface, a killing area high-near interface, a killing area far interface, a killing area high interface, a killing area low interface, a killing area positive side interface, a killing area negative side interface, a killing area upper front positive interface, a killing area upper front negative interface, a killing area lower front positive interface and a killing area lower front negative interface;
a) the method mainly comprises the following steps of determining the DGG 'D' of the low and near interface of the killing area:
1) establishing a low-near-field spherical equation, namely:
in the formula, P is an air route shortcut; h is the target height; q is a target route angle; hminIs the low boundary of the killing area; dsjminMinimum near-bound skew distance for a killing area; sscjIs a near boundary of a killing area; epsilonmaxThe maximum height angle of the killing area; q. q.smaxThe maximum route angle is a killer area; s is a killing area interface;
2) correcting the initial low and near sphere by using the find function to obtain a matrix set (S)scj,P,H)DGG′D′;
3) Using surf (S)scj,P,H)DGG′D′Forming a killing area low and near interface DGG 'D' by a function;
b) the main steps for determining the high-near interface BDD 'B' of the killing area are as follows:
1) establishing a high-near cone equation, namely:
in the formula, HmaxIs the high boundary of the killing area; ssjIs a near boundary of a killing area;
2) correcting the initial high and near boundary conical surface by using the find function to obtain a matrix group (S)sj,P,H)BDD′B′;
3) Invoking surf (S)sj,P,H)BDD′B′A function forming killing area high-near interface BDD 'B';
c) the main steps for determining the far boundary EHH 'E' of the kill zone are as follows:
1) establishing a far-bound surface equation, namely:
in the formula, a, b and c are curved surface equation coefficients;Ssyrepresenting the far boundary of a killing area;
2) the data (H) of m groups of air defense missilesj,Ssyj) Input into p ═ polyfit (H, S)syR) function, and fitting out far boundary of killing area by polynomial fitting methodA non-linear function expression of (a); r is the fitting order; r is less than m; p is a vector formed by polynomial coefficients of the function F; j is 1, …, m;
3) in the interpolation interval [ Hmin,Hmax]Uniformly taking n points and substituting into the far boundary of the killing areaIn the function, a bus (H) of the far boundary of the killing area is linearly fittedi,Ssyi) (ii) a Bus bar (H)i,Ssyi) Around the OH axis of the coordinate system within + -qmaxRotating within the range to obtain an initial far-bound curved surface; 1, …, n;
4) correcting the initial far-bound surface by using a find function to obtain a matrix group (S)sy,P,H)EHH′E′;
5) Using surf (S)sy,P,H)EHH′E′Function form kill zone remote interface EHH 'E';
d) the main steps for determining the high interface HABB ' A ' H ' of the killing area are as follows:
1) establishing a high-boundary plane equation, namely:
2) correcting the initial high-boundary plane equation by using a find function to obtain a matrix set (S, P, H)HABB′A′H′;
3) Using surf (S, P, H)HABB′A′H′Function-forming killer region high interface HABB ' a ' H ':
e) the main steps for determining the killing area low interface GFF 'G' are as follows:
1) establishing a low-boundary plane equation, namely:
2) correcting the initial high-boundary plane equation by using a find function to obtain a matrix set (S, P, H)GFF′G′;
3) Using surf (S, P, H)GFF′G′Forming a killing area low interface GFF 'G' by a function;
f) the main steps for determining the interface on the front side of the killing area are as follows:
1) curve of distance bound
determining a characteristic point E of a curved-edge polygon HAE, wherein the method comprises the following steps: in the interpolation interval [ Hmin,Hmax]Uniformly taking n points and substituting into the far boundary of the killing areaIn the function, a function is obtainedSubstituting n points into function in sequenceIn (1) obtainingThe k point corresponding to the time so as to obtain the coordinate of the E pointHkThe height corresponding to the k point;
2) respectively mapping the curved-edge polygons HAE into (n-k +1) -dimensional matrixes SHAEA (n-k +1) -dimensional matrix PHAEAnd (n-k +1) -dimensional matrix HHAEObtaining a matrix set (S)HAE,PHAE,HHAE);
3) Invoking surf (S)HAE,PHAE,HHAE) Establishing a front side interface of a killing area by a function;
g) the main steps for determining the negative side interface of the killing area are as follows:
1) curve of distance boundStraight lineAnd a straight lineAnd (3) stretching into a curved polygon H ' A ' E ', namely the equation:
2) respectively mapping the curved edge polygons H ' A ' E ' into (n-k +1) dimensional matrixes SH′A′E′A (n-k +1) -dimensional matrix PH′A′E′And (n-k +1) -dimensional matrix HH′A′E′Obtaining a matrix set (S)H′A′E′,PH′A′E′,HH′A′E′);
3) Invoking surf (S)H′A′E′,PH′A′E′,HH′A′E′) Establishing a negative side interface of a killing area by a function;
h) the main steps for determining the front positive interface on the killing area are as follows:
1) plane of maximum route angleThe high interface of the killing area, the positive side interface of the killing area and the high-near interface of the killing area are stretched into a plane polygon ABCD, namely:
2) mapping a planar polygonal ABCD surface to an n-dimensional matrix SABCDN-dimensional matrix PABCDAnd n-dimensional matrix HABCDObtaining a matrix set (S)ABCD,PABCD,HABCD);
3) Invoking surf (S)ABCD,PABCD,HABCD) Obtaining a front positive interface on the killing area by the function;
i) the main steps for determining the front negative interface on the killing area are as follows:
1) plane of maximum route angleThe high boundary of the killing area, the negative side boundary of the killing area and the high-near boundary of the killing area are stretched into a plane polygon A 'B' C 'D', namely:
2) mapping a planar polygon A 'B' C 'D' surface to an n-dimensional matrix SA′B′C′D′N-dimensional matrix PA′B′C′D′And n-dimensional matrix HA′B′C′D′Obtaining a matrix set (S)A′B′C′D′,PA′B′C′D′,HA′B′C′D′);
3) Invoking surf (S)A′B′C′D′,PA′B′C′D′,HA′B′C′D′) Obtaining a front negative interface on the killing area by the function;
j) the main steps for determining the front positive interface under the killing area are as follows:
1) plane of maximum route angleThe low boundary of the killing area, the positive boundary of the killing area and the low-near boundary of the killing area are expanded into a spatial polygon CEFGD, namely:
wherein point F is the starting point of the far field;
2) mapping spatial polygons CEFGD into a n x 3 n-dimensional matrix SCEFGDN x 3n dimensional matrix PCEFGDAnd n x 3n dimensional matrix HCEFGDObtaining a matrix set (S)CEFGD,PCEFGD,HCEFGD);
3) Invoking surf (S)CEFGD,PCEFGD,HCEFGD) Obtaining a front interface under the killing area by the function;
k) the main steps for determining the front negative interface under the killing area are as follows:
1) plane of maximum route angleThe low interface, the front side interface and the low-near interface are flared into a spatial polygon C ' E ' F ' G ' D ', namely:
wherein, the F' point is the starting point of the far boundary;
2) mapping the spatial polygon C ' E ' F ' G ' D ' into a n-x 3 n-dimensional matrix SC′E′F′G′D′N x 3n dimensional matrix PC′E′F′G′D′And n x 3n dimensional matrix HC′E′F′G′D′Obtaining a matrix set (S)C′E′F′G′D′,PC′E′F′G′D′,HC′E′F′G′D′);
3) Invoking surf (S)C′E′F′G′D′,PC′E′F′G′D′,HC′E′F′G′D′) The function yields the front negative interface under the kill zone.
5. The method for constructing the air-defense missile space-killing area model and calculating shooting data base based on the small amount of actual data according to claim 1, wherein the method for calculating the space-killing area model mainly comprises the following steps:
1) establishing a space transmitting area model based on the space killing area model; inputting the air defense missile data into a space launching area model to obtain space launching area parameters;
2) resolving the space killing area model to obtain a horizontal killing area and a transmitting area under the target height;
3) resolving the space killing area model to obtain a vertical killing area and a transmitting area under the short-cut of the route where the target is located;
4) data Cursor function is used for extracting encounter point and emission point Data in a horizontal killing area, a horizontal emission area, a vertical killing area and a vertical emission area to obtain emission residual time tWhen the emission is leftAnd time t of target flying out of emitting areaTime of flightAnd completing shooting data element calculation.
6. The method for constructing the air-defense missile space-killer zone model and shooting data solution based on the small amount of actual data according to claim 1 or 5, wherein the method comprises the following steps: the firing data includes encounter point slope, encounter point height, shot time, time-to-flight, shot distance and shot distance.
7. The method for constructing the air-defense missile space-kill-zone model and calculating shooting data elements based on the small amount of actual data according to claim 1, wherein the method for obtaining the parameters of the space launching zone comprises the following main steps:
1) establishing a space transmitting area model, namely:
Sf=Ss+Vt·tD(12)
in the formula, VtIs the target speed; t is tDThe time of encountering the air defense missile with the target; dfIs the encounter point slope distance; sfIs an emitter region interface;
time t of air defense missile encountering with targetDAs follows:
2) time t for encountering air defense missile with targetD practiceAnd encounter point slope distance Ds actualInput into equation (14) to obtain time calculation coefficient b1Time calculation coefficient b2Time calculation coefficient b3And calculating the time t of the air-defense missile encountering the targetD;
3) Time t for encountering air defense missile with targetDAnd inputting the parameters into a space transmitting area model to obtain space transmitting area parameters.
8. The method for constructing the air-defense missile space-killing area model and shooting data element calculation based on the small amount of actual data according to claim 1, wherein the step of calculating the space-killing area model to obtain the horizontal killing area and the launching area under the target height comprises the following steps:
1) determining the height interval delta H of a cutting horizontal killing area of the space killing area;
2) by interpolation interval [ Hmin,Hmax]Within the range of any height H as a reference and in the height intervalInternally intercepting a horizontal killing area and a transmitting area;
3) and extracting encounter point and emission point Data of the horizontal killing area and the horizontal emission area on the height H by using a Data Cursor function of Malab.
9. The method for constructing the air-defense missile space-killing-area model and calculating shooting data base based on the small amount of actual data according to claim 1, wherein the step of calculating the space-killing-area model to obtain the vertical killing area and the launching area under the route shortcut where the target is located is as follows:
1) determining an airway shortcut interval delta P for cutting the emitting area of the vertical killing area by the emitting area of the space killing area;
2) taking a short-cut range of the air route [ -P [ ]max,Pmax]Taking any route shortcut P in the range as a reference, and taking the route shortcut P as a referenceInternally intercepting a vertical killing area and a transmitting area on the route shortcut;
3) and extracting encounter points and emission point Data of the vertical killing area and the vertical emission area on the airway shortcut P by using a Data Cursor function of Malab.
10. The method for constructing the air-defense missile space-killer zone model and calculating shooting data elements based on the small amount of actual data according to claim 1, wherein the firing residue is as follows:
wherein s represents the s coordinate of the target in the shooting coordinate system O-SPH; dFYA far-bound slant range of the emitting area; dFJThe near-boundary slant distance of the emission region; sfyA far boundary of the emission region; sfjIs the near boundary of the emission region; t is t0Time of the target at the current position;
the target departure launch time is as follows:
ttime of flight=tFJ-tFY(16)
Wherein the time t when the target reaches the near boundary of the emission regionFJAs follows:
in the formula, v is a target speed.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911351195.3A CN111191359B (en) | 2019-12-24 | 2019-12-24 | Method for constructing space killer zone model of air-defense missile and calculating shooting data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911351195.3A CN111191359B (en) | 2019-12-24 | 2019-12-24 | Method for constructing space killer zone model of air-defense missile and calculating shooting data |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111191359A true CN111191359A (en) | 2020-05-22 |
CN111191359B CN111191359B (en) | 2021-03-23 |
Family
ID=70709410
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911351195.3A Active CN111191359B (en) | 2019-12-24 | 2019-12-24 | Method for constructing space killer zone model of air-defense missile and calculating shooting data |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111191359B (en) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112417653A (en) * | 2020-11-03 | 2021-02-26 | 中国人民解放军海军工程大学 | Method for establishing ship-air missile interception model |
CN112504016A (en) * | 2020-09-21 | 2021-03-16 | 上海航天控制技术研究所 | Target non-escape area reliable prediction method adaptive to collaborative task planning |
CN112896558A (en) * | 2021-03-19 | 2021-06-04 | 湖南工学院 | Manual control rendezvous and docking human error identification method based on multi-dimensional data |
CN113361196A (en) * | 2021-06-04 | 2021-09-07 | 上海机电工程研究所 | Missile killing probability evaluation method, system, equipment and readable medium |
CN113483739A (en) * | 2021-07-14 | 2021-10-08 | 中国科学院长春光学精密机械与物理研究所 | Offshore target position measuring method |
CN113706697A (en) * | 2021-09-06 | 2021-11-26 | 北京航空航天大学 | Method and system for determining shooting precision of spherical target object |
CN117371256A (en) * | 2023-12-07 | 2024-01-09 | 中国人民解放军海军航空大学 | Multi-type multi-number platform deployment setting scheme planning method |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100297589A1 (en) * | 2006-08-18 | 2010-11-25 | Saab Ab | Device arranged for illuminate an area |
CN105224778A (en) * | 2014-06-09 | 2016-01-06 | 上海机电工程研究所 | General killing area computing method and general launch site computing method |
CN109357578A (en) * | 2018-10-24 | 2019-02-19 | 刘勇 | School shooting method is aimed at based on the calibration shell of Beidou navigation technology and artilleryman |
CN110597056A (en) * | 2019-08-16 | 2019-12-20 | 南京理工大学 | Large closed-loop calibration control method for antiaircraft gun fire control system |
-
2019
- 2019-12-24 CN CN201911351195.3A patent/CN111191359B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100297589A1 (en) * | 2006-08-18 | 2010-11-25 | Saab Ab | Device arranged for illuminate an area |
CN105224778A (en) * | 2014-06-09 | 2016-01-06 | 上海机电工程研究所 | General killing area computing method and general launch site computing method |
CN109357578A (en) * | 2018-10-24 | 2019-02-19 | 刘勇 | School shooting method is aimed at based on the calibration shell of Beidou navigation technology and artilleryman |
CN110597056A (en) * | 2019-08-16 | 2019-12-20 | 南京理工大学 | Large closed-loop calibration control method for antiaircraft gun fire control system |
Non-Patent Citations (4)
Title |
---|
HUI YAOLUO等: "Dynamic attack zone of air-to-air missile after being launched in random wind field", 《CHINESE JOURNAL OF AERONAUTICS》 * |
张月星等: ""弹炮结合武器系统总体设计研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
王君等: "地空导弹杀伤区数值化模型研究", 《系统工程理论与实践》 * |
耿杰恒: "地空导弹仿真建模中射击诸元计算模型研究", 《舰船电子工程》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112504016A (en) * | 2020-09-21 | 2021-03-16 | 上海航天控制技术研究所 | Target non-escape area reliable prediction method adaptive to collaborative task planning |
CN112417653A (en) * | 2020-11-03 | 2021-02-26 | 中国人民解放军海军工程大学 | Method for establishing ship-air missile interception model |
CN112896558A (en) * | 2021-03-19 | 2021-06-04 | 湖南工学院 | Manual control rendezvous and docking human error identification method based on multi-dimensional data |
CN113361196A (en) * | 2021-06-04 | 2021-09-07 | 上海机电工程研究所 | Missile killing probability evaluation method, system, equipment and readable medium |
CN113483739A (en) * | 2021-07-14 | 2021-10-08 | 中国科学院长春光学精密机械与物理研究所 | Offshore target position measuring method |
CN113483739B (en) * | 2021-07-14 | 2022-05-13 | 中国科学院长春光学精密机械与物理研究所 | Offshore target position measuring method |
CN113706697A (en) * | 2021-09-06 | 2021-11-26 | 北京航空航天大学 | Method and system for determining shooting precision of spherical target object |
CN113706697B (en) * | 2021-09-06 | 2023-05-12 | 北京航空航天大学 | Method and system for determining shooting precision of spherical target |
CN117371256A (en) * | 2023-12-07 | 2024-01-09 | 中国人民解放军海军航空大学 | Multi-type multi-number platform deployment setting scheme planning method |
CN117371256B (en) * | 2023-12-07 | 2024-03-29 | 中国人民解放军海军航空大学 | Multi-type multi-number platform deployment setting scheme planning method |
Also Published As
Publication number | Publication date |
---|---|
CN111191359B (en) | 2021-03-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111191359B (en) | Method for constructing space killer zone model of air-defense missile and calculating shooting data | |
CN110608743B (en) | Multi-unmanned aerial vehicle collaborative route planning method based on multi-population chaotic grayling algorithm | |
CN108549233B (en) | Unmanned aerial vehicle air combat maneuver game method with intuitive fuzzy information | |
CN103557867B (en) | The collaborative path planning method of a kind of many UAV of three-dimensional based on sparse A* search | |
CN108153328B (en) | Multi-missile collaborative track planning method based on segmented Bezier curve | |
CN107392388A (en) | A kind of method for planning no-manned plane three-dimensional flight path using artificial fish-swarm algorithm is improved | |
CN107063255A (en) | A kind of three-dimensional Route planner based on improvement drosophila optimized algorithm | |
CN103942401A (en) | Tool kit and method for optimizing high-precision self-adaptation and modular spacecraft trajectory multi-constrained track | |
CN112666981B (en) | Unmanned aerial vehicle cluster dynamic route planning method based on dynamic group learning of original pigeon group | |
CN112861257B (en) | Aircraft fire control system precision sensitivity analysis method based on neural network | |
CN105608251A (en) | BNSobol method for sensitivity analysis on precision of Helicopter fire control system | |
CN116225066A (en) | Unmanned aerial vehicle path optimization method based on chaotic mapping pely optimization algorithm | |
CN109597839A (en) | A kind of data digging method based on the avionics posture of operation | |
CN111121784A (en) | Unmanned reconnaissance aircraft route planning method | |
CN109558934A (en) | A kind of unmanned plane interpolation moth flutters flame low-level penetration method | |
CN116680509A (en) | Dynamic matching method for multi-spacecraft escape-tracking game task | |
Li et al. | Multi-uav air combat weapon-target assignment based on genetic algorithm and deep learning | |
CN118094779A (en) | Guided rocket large airspace glide increase Cheng Dandao optimization method based on ant colony algorithm | |
CN110728010A (en) | Method for calculating firing data of antiaircraft gun fire control system | |
CN113177583A (en) | Aerial target clustering method | |
CN115357051B (en) | Deformation and maneuvering integrated avoidance and defense method | |
CN107329935B (en) | Solving method for air situation pareto attack and defense strategy | |
Duan et al. | Route planning method design for UAV under radar ECM scenario | |
CN115186378A (en) | Real-time solution method for tactical control distance in air combat simulation environment | |
CN105096379A (en) | Triangular mesh surface reconstruction method based on k-neighborhood |
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 |