Simulation arraying method of radar active induced polarization system facing anti-radiation missile

Abstract

The invention discloses a simulation arraying method of a radar active induced polarization system facing an anti-radiation missile, and relates to the technical field of radar active induced polarization arraying. According to the method, the attack process simulation of an anti-radiation missile is realized based on a quaternion method, an anti-radiation missile motion model is constructed, and the four-point-source radar active decoy simulation system is obtained by drawing up different array distances under different array forms of the four-point-source radar active decoy system. Through obtaining an impact point of an anti-radiation missile motion model in the four-point-source radar active decoy simulation system and matching with judgment of an effective killing range, calculation of a survival probability and a maximum induced deviation distance of the four-point-source radar active decoy simulation system is realized according to a probability density function. And determining an optimal arrangement mode and an optimal arrangement distance between the three baits and the radar so as to achieve the purpose of ensuring the maximum survival probability of the radar and the baits while performing maximum coverage on the airspace.

Description

Technical Field

[0001] The invention relates to the technical field of radar active decoy array deployment, in particular to a simulation array deployment method of a radar active decoy system for an anti-radiation missile. Background Art

[0002] Anti-radiation missiles track and attack radiation sources through passive seekers. With the rapid development of anti-radiation missile technology, the safety of active equipment and the lives of radar operators on the battlefield are under great threat, especially the safety of tracking and guidance radars. Anti-radiation missiles are the effective and main means of attacking radar positions, and active decoy systems to counter them have also developed rapidly.

[0003] It is generally believed that after the power gain product of the active decoy covers the side lobes of the radar, it can achieve the effect of decoying anti-radiation missiles by radiating pulse signals similar to radar signals. However, in actual combat use, it is necessary to consider the deployment form and deployment distance of the active decoy system and the radar to achieve maximum coverage of the airspace and maximize the survival probability of the radar and the decoy. For this reason, a simulation deployment method of the radar active decoy system for anti-radiation missiles is proposed, and the deployment form and deployment distance are converted into data, which provides reliable support for achieving maximum coverage of the airspace and maximizing the survival probability of the radar and the decoy. Summary of the invention

[0004] 1. Technical issues to be solved

[0005] In view of the deficiencies in the prior art, the present invention provides a simulation array method for a radar active decoy system for an anti-radiation missile, which solves the above-mentioned problems.

[0006] (II) Technical solution

[0007] To achieve the above object, the present invention is implemented by the following technical scheme: a simulation array method for a radar active decoy system for an anti-radiation missile, specifically comprising the following steps:

[0008] Step 1: Establish a four-point source radar active decoy system including three decoy points and one radar point, and obtain an anti-radiation missile attack process based on the four-point source radar active decoy system;

[0009] Step 2, different array distances under different array forms are formulated for the four-point source radar active deception system to obtain a four-point source radar active deception simulation system;

[0010] Step 3: Establish an anti-radiation missile motion model according to the anti-radiation missile attack process;

[0011] Step 4: Calculate the survival probability of radar points and decoy points using probability density function;

[0012] Step 5: Select the formation with the highest survival probability, and analyze the impact of different formation distances on the survival probability and maximum decoy distance under this formation to obtain the optimal formation distance under high survival probability.

[0013] The present invention is further configured as follows: the anti-radiation missile attack process based on the four-point source radar active deception system in step 1 includes:

[0014] The height of the anti-radiation missile is H. At this time, the axis of the anti-radiation missile points to the energy center O of the radiation source composed of three decoy points and one radar point. The attack angle is α. The distance between the radar point and the energy center O is L. The final impact point of the anti-radiation missile is O1. When the passive radar seeker on the anti-radiation missile is in a random equilibrium position, the passive radar seeker identifies the target position. When the passive radar seeker identifies the radar point position, it is at the maximum overload α. M Turn toward the radar point, and the maximum horizontal distance turned is the horizontal distance between the energy center point O and the final impact point point O1.

[0015] The present invention is further configured as follows: the array forms in step 2 include star array, triangle array and diamond array.

[0016] The present invention is further configured as follows: the establishment of the anti-radiation missile motion model in step 3 includes:

[0017] The position coordinates of the anti-radiation missile are A k (x k ,y k ,z k ), at this time, the energy center measured by the passive radar seeker and the aiming direction formed by the passive radar seeker are The current speed direction is After time Δt, the anti-radiation missile flies to A k+1 (x k+1 ,y k+1 ,z k+1 );

[0018] The speed direction With aiming direction The angle is α k :

[0019]

[0020] and The angle α is the angle that the anti-radiation missile can actually turn z :

[0021] α z =min{α k,α max};

[0022] In the formula, α max The maximum angle that can be rotated theoretically is:

[0023]

[0024] According to the angle α that the anti-radiation missile can actually turn z , use the quaternion method to rotate the anti-radiation missile and update the position coordinate A of the anti-radiation missile k+1 (x k+1 ,y k+1 ,z k+1 ).

[0025] The present invention is further configured as follows: the anti-radiation missile position coordinate A k+1 (x k+1 ,y k+1 ,z k+1 ) include:

[0026] A1. Calculate the quaternion rotation axis

[0027]

[0028] A2. Calculate the rotation vector of the quaternion

[0029]

[0030] A3. Calculate the speed direction:

[0031]

[0032] A4. Calculate length

[0033]

[0034] A5. Update the position coordinates of the anti-radiation missile A k+1 (x k+1 ,y k+1 ,z k+1 ):

[0035]

[0036] A6. Determine whether the anti-radiation missile has landed. After landing, stop updating the position coordinates. Otherwise, if it has not landed, repeat A1-A5 to update the position coordinates.

[0037] The present invention is further configured as follows: the method of calculating the survival probabilities of radar points and decoy points by using the probability density function in step 4 includes:

[0038] The position coordinates of the anti-radiation missile after landing are used as the impact point, and the effective killing radius of the anti-radiation missile is set. When any of the three decoy points and one radar point appears within the effective killing radius, it is determined that the corresponding decoy point or radar point cannot survive. Conversely, when none of the three decoy points and one radar point appears within the effective killing radius, it is determined that the decoy point and the radar point survive.

[0039] The coordinates of the impact point are A end = [XYZ], the estimated mean of the impact point in the X direction of the coordinate axis is μ x :

[0040]

[0041] Calculate the standard deviation σ x :

[0042]

[0043] The estimated mean value of the impact point in the Y direction of the coordinate axis is μ y :

[0044]

[0045] Calculate the standard deviation σ y ;

[0046]

[0047] In the formula, x i ,y i is the coordinate of the landing point in the X and Y directions in the i-th trial, and n is the number of trials;

[0048] Calculate the correlation coefficient r:

[0049]

[0050] Get the probability density function:

[0051]

[0052] According to the probability density function, the radar point and the three decoy points are integrated in the circular area of ​​the effective killing radius of the anti-radiation missile to obtain the survival probability of the radar point and the three decoy points:

[0053] P = 1-∫∫f(x,y)dxdy.

[0054] The present invention is further configured as follows: the maximum decoy distance in step five is the maximum distance between the impact point of the anti-radiation missile and the radar point, and the maximum distance between the impact point of the anti-radiation missile and the decoy point.

[0055] (III) Beneficial effects

[0056] The present invention provides a simulation array method for a radar active decoy system for an anti-radiation missile, which has the following beneficial effects:

[0057] The present invention realizes the attack process simulation of the anti-radiation missile based on the quaternion method, constructs the anti-radiation missile motion model, obtains the four-point source radar active decoy simulation system by formulating different array distances under different array forms for the four-point source radar active decoy simulation system, obtains the impact point of the anti-radiation missile motion model in the four-point source radar active decoy simulation system, cooperates with the judgment of the effective killing range, realizes the calculation of the survival probability and the maximum decoy distance of the four-point source radar active decoy simulation system according to the probability density function, and determines the optimal array mode and array distance between three decoys and the radar, so as to achieve the purpose of maximizing the survival probability of the radar and the decoy while maximizing the coverage of the airspace. BRIEF DESCRIPTION OF THE DRAWINGS

[0058] Figure 1 It is a schematic diagram of the process of the present invention;

[0059] Figure 2 A schematic diagram of an anti-radiation missile attack radar position in an embodiment of the present invention;

[0060] Figure 3 This is a schematic diagram of the anti-radiation missile attacking a radiation source target in an embodiment of the present invention;

[0061] Figure 4 It is a schematic diagram of flight adjustment of an anti-radiation missile in an embodiment of the present invention;

[0062] Figure 5 It is a line graph of array distance and survival probability in an embodiment of the present invention. DETAILED DESCRIPTION

[0063] The technical solutions in the embodiments of the present invention will be described clearly and completely below in conjunction with the accompanying drawings in the embodiments of the present invention.

[0064] See also Figure 1-5 , the embodiment of the present invention provides the following technical solutions:

[0065] Embodiment 1

[0066] A simulation array method for a radar active decoy system for an anti-radiation missile specifically comprises the following steps:

[0067] Step 1: Establish a four-point source radar active decoy system including three decoy points and one radar point, and obtain an anti-radiation missile attack process based on the four-point source radar active decoy system;

[0068] Specifically, the anti-radiation missile uses a passive radar seeker to receive and track the radiation source, where the radiation source is the three decoy points and the radar point. Due to its technical characteristics, the passive radar seeker will track the energy center of the radiation source, which is the midpoint of the radiation area of ​​the three decoy points and the radar point. When the anti-radiation missile is far away from the radiation source, the passive radar seeker resolution angle Δθ R All radiation sources can be covered. At this time, the passive radar seeker tracks the energy center of the radiation source. During the flight, the angle between one or more radiation sources and the aiming direction of the passive radar seeker will increase, and the angle is equal to Δθ R / 2, the seeker reaches the random balance point, and the next time the radiation source cannot deflect the passive radar seeker. Therefore, during the anti-radiation weapon attack, the radiation sources in the field of view change from multiple to one, and finally a target is selected for attack. Because of its fast flight speed, when it can distinguish the target and attack a certain radiation source, the distance is often short. In addition, the overload capacity of the anti-radiation missile is limited by the structure and material, and it is difficult to correct the initial value of the tracking error. It will eventually deviate from the target and land. The attack principle is as shown in the attached figure. Figure 2 As shown;

[0069] To further explain, the anti-radiation missile attack process based on the four-point source radar active decoy system includes:

[0070] The height of the anti-radiation missile is H. At this time, the axis of the anti-radiation missile points to the energy center O of the radiation source composed of three decoy points and one radar point. The attack angle is α. The distance between the radar point and the energy center O is L. The final impact point of the anti-radiation missile is O1. When the passive radar seeker on the anti-radiation missile is in a random equilibrium position, the passive radar seeker identifies the target position. When the passive radar seeker identifies the radar point position, it is at the maximum overload α. M Turn towards the radar point. The maximum horizontal distance covered is the horizontal distance between the energy center point O and the final impact point point O1.

[0071] Step 2, different array distances under different array forms are formulated for the four-point source radar active deception system to obtain a four-point source radar active deception simulation system;

[0072] The array forms include star array, triangle array and diamond array;

[0073] The star array is as follows: three decoy points form a triangular area, and the radar point is located inside the triangular area;

[0074] The triangular array is as follows: three decoy points form a triangular area, and the radar point is on the edge of the triangular area;

[0075] The diamond array is as follows: three bait points and radar points form a flat diamond area.

[0076] Step 3: Establish an anti-radiation missile motion model based on the anti-radiation missile attack process:

[0077] The position coordinates of the anti-radiation missile are A k (x k ,y k ,z k ), at this time, the energy center measured by the passive radar seeker and the aiming direction formed by the passive radar seeker are The current speed direction is After time Δt, the anti-radiation missile flies to A k+1 (x k+1 ,y k+1 ,z k+1 );

[0078] The speed direction With aiming direction The angle is α k :

[0079]

[0080] and The angle α is the angle that the anti-radiation missile can actually turn z :

[0081] α z =min{α k ,α max};

[0082] In the formula, α max The maximum angle that can be rotated theoretically is:

[0083]

[0084] According to the angle α that the anti-radiation missile can actually turn z , use the quaternion method to rotate the anti-radiation missile and update the position coordinate A of the anti-radiation missile k+1 (x k+1 ,y k+1 ,z k+1 );

[0085] Anti-radiation missile position coordinates A k+1 (x k+1 ,y k+1 ,z k+1 ) include:

[0086] A1. Calculate the quaternion rotation axis

[0087]

[0088] A2. Calculate the rotation vector of the quaternion

[0089]

[0090] A3. Calculate the speed direction:

[0091]

[0092] A4. Calculate length

[0093]

[0094] A5. Update the position coordinates of the anti-radiation missile A k+1 (x k+1 ,y k+1 ,z k+1 ):

[0095]

[0096] A6, determine whether the anti-radiation missile has landed. After landing, stop updating the position coordinates. Otherwise, if it has not landed, repeat A1-A5 to update the position coordinates.

[0097] Step 4: Calculate the survival probability of radar points and decoy points using probability density function: Take the position coordinates of the anti-radiation missile after landing as the impact point, set the effective killing radius of the anti-radiation missile, and when any of the three decoy points and one radar point appears within the effective killing radius, determine that the corresponding decoy point or radar point cannot survive; conversely, when none of the three decoy points and one radar point appears within the effective killing radius, determine that the decoy point and the radar point survive;

[0098] The coordinates of the impact point are A end = [XYZ], the estimated mean of the impact point in the X direction of the coordinate axis is μ x :

[0099]

[0100] Calculate the standard deviation σ x :

[0101]

[0102] The estimated mean value of the impact point in the Y direction of the coordinate axis is μ y :

[0103]

[0104] Calculate the standard deviation σ y ;

[0105]

[0106] In the formula, x i ,y i is the coordinate of the landing point in the X and Y directions in the i-th trial, and n is the number of trials;

[0107] Calculate the correlation coefficient r:

[0108]

[0109] Get the probability density function:

[0110]

[0111] According to the probability density function, the radar point and the three decoy points are integrated in the circular area of ​​the effective killing radius of the anti-radiation missile to obtain the survival probability of the radar point and the three decoy points:

[0112] P = 1-∫∫f(x,y)dxdy;

[0113] Step 5: Filter the formation with the highest survival probability, and analyze the influence of different formation distances on the survival probability and the maximum decoy distance under this formation, and obtain the optimal formation distance under high survival probability, where the maximum decoy distance is the maximum distance between the impact point of the anti-radiation missile and the radar point, and the maximum distance between the impact point of the anti-radiation missile and the decoy point.

[0114] Embodiment 2

[0115] According to the method of embodiment 1, a triangular array is used for array formation, and the anti-radiation missile flight speed is set to 1000m / s, the maximum longitudinal overload is 3g, the anti-radiation missile field of view is 60 degrees, and the missile killing radius is 30m. The anti-radiation missile attack direction is in pitch (20°-90°) and azimuth (0°-360°). In the simulation, 10 points are taken within the pitch range, and 100 points are evenly taken within the azimuth unit, and 1000 attack tests are conducted.

[0116] Among them, the radar coordinates are (0, 0, 0) m, the decoy 1 coordinates are (120, 0, 0) m, the decoy 2 coordinates are (0, 207, 0) m, and the decoy 3 coordinates are (-120, 0, 0) m. In the form of a triangular array, the maximum decoy distances between the impact point and the radar, decoy 1, decoy 2, and decoy 3 in 1000 attack tests are 162.36 m, 256.25 m, 310.70 m, and 261.57 m, respectively.

[0117] The effective killing results of radar, decoy 1, decoy 2 and decoy 3 in 1000 attack tests were normalized. The results showed that the survival probabilities of radar, decoy 1, decoy 2 and decoy 3 were 52.04%, 99.96%, 100% and 99.96% respectively.

[0118] Embodiment 3

[0119] According to the method of embodiment 1, the star array is used for array formation, the anti-radiation missile flight speed is set to 1000m / s, the maximum longitudinal overload is 3g, the anti-radiation missile field of view is 60 degrees, and the missile killing radius is 30m. The anti-radiation missile attack direction is in pitch (20°-90°) and azimuth (0°-360°). In the simulation, 10 points are taken in the pitch range, and 100 points are evenly taken in the azimuth unit, and 1000 attack tests are carried out.

[0120] Among them, the radar coordinates are (0, 51, 0) m, the decoy 1 coordinates are (120, 0, 0) m, the decoy 2 coordinates are (0, 207, 0) m, and the decoy 3 coordinates are (-120, 0, 0) m. In the star array form, the maximum decoy distances between the impact point and the radar, decoy 1, decoy 2, and decoy 3 in 1000 attack tests are 130.51 m, 250.76 m, 267.99 m, and 255.12 m respectively;

[0121] The effective killing results of radar, decoy 1, decoy 2 and decoy 3 in 1000 attack tests were normalized. The results showed that the survival probabilities of radar, decoy 1, decoy 2 and decoy 3 were 43.44%, 100%, 99.96% and 100% respectively.

[0122] Embodiment 3

[0123] According to the method of embodiment 1, a diamond array is used for array formation, and the anti-radiation missile flight speed is set to 1000m / s, the maximum longitudinal overload is 3g, the anti-radiation missile field of view is 60 degrees, and the missile killing radius is 30m. The anti-radiation missile attack direction is in pitch (20°-90°) and azimuth (0°-360°). In the simulation, 10 points are taken within the pitch range, and 100 points are evenly taken within the azimuth unit, and 1000 attack tests are conducted.

[0124] Among them, the radar coordinates are (120, 0, 0) m, the decoy 1 coordinates are (-120, 0, 0) m, the decoy 2 coordinates are (0, 207, 0) m, and the decoy 3 coordinates are (0, -207, 0) m. In the form of diamond array, the maximum decoy distances between the impact point and the radar, decoy 1, decoy 2, and decoy 3 in 1000 attack tests are 400.12 m, 295.31 m, 414.29 m, and 305.04 m respectively;

[0125] The effective killing results of radar, decoy 1, decoy 2 and decoy 3 in 1000 attack tests were normalized. The results showed that the survival probabilities of radar, decoy 1, decoy 2 and decoy 3 were 99.8%, 99.78%, 99.98% and 99.98% respectively.

[0126] Judging from the results of the survival probability in the above-mentioned Examples 2, 3 and 4, it can be seen that the survival probability and the maximum decoy distance of the radar under the star array are the worst, with the survival probability being only 43.44% and the maximum decoy distance being only 130.51m. The radar survival probability and the maximum decoy distance under the diamond array are optimal, with the survival probability reaching 99.66% and the maximum decoy distance reaching 400.12m. In the triangular array and star array, although the survival probability can reach more than 99%, these two maximum decoy distances are both smaller than the maximum decoy distance under the diamond array. Therefore, from the perspective of the survival probability and the maximum decoy distance of the radar and the decoy, the diamond array is the optimal array form.

[0127] Embodiment 4:

[0128] Based on Example 3, a diamond array is adopted. The difference is that the distances between the radar, bait 1, bait 2 and bait 3 are set to 100m, 160m, 200m, 240m, 280m, 320m, 360m, 400m, 440m, 480m and 520m respectively, and tests are carried out at different array distances.

[0129] The relationship between the formation distance and the survival probability is shown in the figure below. Figure 5 As shown below, it can be seen that the survival probability of radar and decoy has a trend of first increasing and then decreasing with increasing distance. In the initial stage, although the anti-radiation missile has limited maneuverability, due to the close array distance, the probability of killing the radiation source after the anti-radiation missile maneuvers is relatively high; as the distance increases, the anti-radiation missile cannot kill the radiation source after being deflected. In the distance segment of 240m-320m, the survival probability of radar and decoy is in a higher range; as the distance further increases, the decoy's coverage of the airspace decreases, and the passive radar seeker can more easily identify the radiation source, so the survival probability decreases instead.

Claims

1. A simulation array method for a radar active decoy system for an anti-radiation missile, characterized in that: The specific steps include: Step 1: Establish a four-point source radar active decoy system including three decoy points and one radar point, and obtain an anti-radiation missile attack process based on the four-point source radar active decoy system; Step 2, different array distances under different array forms are formulated for the four-point source radar active deception system to obtain a four-point source radar active deception simulation system; Step 3: Establish an anti-radiation missile motion model according to the anti-radiation missile attack process; Step 4: Calculate the survival probability of radar points and decoy points using probability density function; Step 5: Select the formation with the highest survival probability, and analyze the impact of different formation distances on the survival probability and maximum decoy distance under this formation to obtain the optimal formation distance under high survival probability.

2. The simulation array method of a radar active decoy system for anti-radiation missiles according to claim 1 is characterized in that: The anti-radiation missile attack process based on the four-point source radar active deception system in step 1 includes: The height of the anti-radiation missile is H. At this time, the axis of the anti-radiation missile points to the energy center O of the radiation source composed of three decoy points and one radar point. The attack angle is α. The distance between the radar point and the energy center O is L. The final impact point of the anti-radiation missile is O1. When the passive radar seeker on the anti-radiation missile is in a random equilibrium position, the passive radar seeker identifies the target position. When the passive radar seeker identifies the radar point position, it is at the maximum overload α. M Turn toward the radar point, and the maximum horizontal distance turned is the horizontal distance between the energy center point O and the final impact point point O1.

3. The simulation array method of a radar active decoy system for anti-radiation missiles according to claim 1 is characterized in that: The array forms in step 2 include star array, triangle array and diamond array.

4. The simulation array method of a radar active decoy system for anti-radiation missiles according to claim 1 is characterized in that: The establishment of the anti-radiation missile motion model in step 3 includes: The position coordinates of the anti-radiation missile are A k (x k ,y k ,z k ), at this time, the energy center measured by the passive radar seeker and the aiming direction formed by the passive radar seeker are The current speed direction is After time Δt, the anti-radiation missile flies to A k+1 (x k+1 ,y k+1 ,z k+1 ); The speed direction With aiming direction The angle is α k : and The angle α is the angle that the anti-radiation missile can actually turn z : a z =min{a k ,a max }; In the formula, α max The maximum angle that can be rotated theoretically is: According to the angle α that the anti-radiation missile can actually turn z , use the quaternion method to rotate the anti-radiation missile and update the position coordinate A of the anti-radiation missile k+1 (x k+1 ,y k+1 ,z k+1 ).

5. The simulation array method of the radar active decoy system for anti-radiation missiles according to claim 4 is characterized in that: The anti-radiation missile position coordinate A k+1 (x k+1 ,y k+1 ,z k+1 ) include: A1. Calculate the quaternion rotation axis A2. Calculate the rotation vector of the quaternion A3. Calculate the speed direction: A4. Calculate length A5. Update the position coordinates of the anti-radiation missile A k+1 (x k+1 ,y k+1 ,z k+1 ): A6. Determine whether the anti-radiation missile has landed. After landing, stop updating the position coordinates. Otherwise, if it has not landed, repeat A1-A5 to update the position coordinates.

6. The simulation array method of the radar active decoy system for anti-radiation missiles according to claim 5 is characterized in that: The method of calculating the survival probability of the radar point and the bait point by using the probability density function in step 4 includes: The position coordinates of the anti-radiation missile after landing are used as the impact point, and the effective killing radius of the anti-radiation missile is set. When any of the three decoy points and one radar point appears within the effective killing radius, it is determined that the corresponding decoy point or radar point cannot survive. Conversely, when none of the three decoy points and one radar point appears within the effective killing radius, it is determined that the decoy point and the radar point survive. The coordinates of the impact point are A end = [XYZ], the estimated mean of the impact point in the X direction of the coordinate axis is μ x : Calculate the standard deviation σ x : The estimated mean value of the impact point in the Y direction of the coordinate axis is μ y : Calculate the standard deviation σ y ; In the formula, x i ,y i is the coordinate of the landing point in the X and Y directions in the i-th trial, and n is the number of trials; Calculate the correlation coefficient r: Get the probability density function: According to the probability density function, the radar point and the three decoy points are integrated in the circular area of ​​the effective killing radius of the anti-radiation missile to obtain the survival probability of the radar point and the three decoy points: P = 1-∫∫f(x,y)dxdy.

7. The simulation array method of a radar active decoy system for anti-radiation missiles according to claim 1 is characterized in that: The maximum decoy distance in step five is the maximum distance between the impact point of the anti-radiation missile and the radar point, and the maximum distance between the impact point of the anti-radiation missile and the decoy point.

Images