CN107917646B - Infrared air-to-air missile anti-interference guidance method based on target terminal reachable area prediction - Google Patents
Infrared air-to-air missile anti-interference guidance method based on target terminal reachable area prediction Download PDFInfo
- Publication number
- CN107917646B CN107917646B CN201710018102.XA CN201710018102A CN107917646B CN 107917646 B CN107917646 B CN 107917646B CN 201710018102 A CN201710018102 A CN 201710018102A CN 107917646 B CN107917646 B CN 107917646B
- Authority
- CN
- China
- Prior art keywords
- missile
- target
- air
- calculating
- next moment
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F42—AMMUNITION; BLASTING
- F42B—EXPLOSIVE CHARGES, e.g. FOR BLASTING, FIREWORKS, AMMUNITION
- F42B15/00—Self-propelled projectiles or missiles, e.g. rockets; Guided missiles
- F42B15/01—Arrangements thereon for guidance or control
Abstract
The anti-interference guidance method of the infrared air-to-air missile based on the target future position prediction comprises the following steps: (1) the infrared air-to-air missile seeker obtains measurement information of a target, filter estimation is carried out on target state information, and the residual flight time of the missile is estimated; (2) assuming that the target acceleration follows the probability distribution of gamma distribution, and calculating the position distribution of the target at the next moment and the meeting moment; (3) calculating the position of the missile which possibly appears at the next moment, and calculating the region which can be reached at the meeting moment by the maximum flight capacity at the position of the next moment; (4) calculating the position of the missile at the next moment corresponding to the area covered by the missile with the maximum probability by comparing the coverage area of the missile with the reachable area of the target; (5) calculating a control instruction of the missile according to the position of the next moment; (6) and (3) exerting control on the missile, and returning to the step (1) until the missile meets the target.
Description
Technical Field
The invention relates to an anti-interference guidance algorithm of an infrared air-to-air missile under the condition of target launching interference, in particular to an anti-interference guidance algorithm under the condition that true and false targets cannot be distinguished.
Background
After a target aircraft is threatened by a missile, an interference missile similar to the spectral characteristic of the aircraft is thrown to improve the survival rate of the aircraft, and the infrared air-to-air missile is used for accurately hitting the target, needing to be anti-interfered by hardware and an algorithm and finally hitting the target accurately.
An infrared imaging guidance anti-interference analysis by Jiakurui et al and an anti-interference technology by an infrared imaging guidance head by Li Juan et al all introduce in detail how an infrared guidance head resists interference when being thrown into a target airplane, the target and the interference are identified by analyzing differences in characteristics such as frequency spectrum, motion, gray scale and spatial distribution, an infrared air-to-air Missile guidance method under the condition of interference is researched in an Arthur Vermeulen paper Missile avaidance manoeuvres with Simultaneous export discovery deployment, and the miss distance relationship of the guidance method for tracking an energy center under different operational environments is analyzed. Dionne studied the Guidance method based on Predictive control in case of a single chaff in the paper "Predictive Guidance for Pursuit-evolution engage investments Decoys". However, with the development of the infrared interference missile, the spectral characteristic and the motion characteristic of the infrared interference missile are more and more similar to those of an airplane, the interference is greatly improved from the appearance and the motion state, particularly when the accuracy of an infrared imaging seeker at a long distance is not high, a certain difficulty and error rate exist in identifying a target airplane and the interference missile through the infrared air-to-air missile seeker, and in order to improve the anti-interference capability of the infrared air-to-air missile and prevent the missile from being trapped by the interference missile, the guidance algorithm of the infrared air-to-air missile is required to be improved.
In order to improve the success rate of intercepting the target by the missile under the condition that the target cannot be accurately identified, the invention provides a guidance method based on prediction guidance. The guidance method is based on a predicted target and a future reachable area of an interfering missile, combines the reachable range of the missile, and calculates a guidance instruction by taking the reached area of a true target and a false target which cannot be distinguished covered by the maximum probability as an index, so that the target can be accurately hit under the condition that the target cannot be distinguished.
Disclosure of Invention
The technical problem to be solved by the invention is as follows: the anti-interference guidance method for the infrared air-to-air missile is provided, and the hit rate of the infrared air-to-air missile after the infrared air-to-air missile is interfered by a target.
The technical solution of the invention is as follows: an anti-interference guidance method of an infrared air-to-air missile based on prediction guidance comprises the following steps:
(1) the infrared air-to-air missile seeker obtains measurement information of a target, filter estimation is carried out on target state information, and the residual flight time of the missile is estimated;
(2) assuming that the target acceleration follows the probability distribution of gamma distribution, and calculating the position distribution region of the target at the next moment and the meeting moment;
(3) calculating the reachable area of the missile at the next moment, and calculating the maximum flight capacity of the missile at each possible position at the next moment at the final moment tfAn area that can be reached;
(4) providing a missile control instruction which is calculated by taking a maximum probability hit target of the missile as a performance index and maximizing the performance;
(6) and (3) applying control to the missile, and returning to the step (1) for circulation until the missile meets the target.
Compared with the prior art, the invention has the advantages that:
(1) compared with the traditional anti-interference method, the method can be used in the interference identification end section, and can supplement the guidance problem that the infrared air-to-air missile cannot identify the true and false targets according to other characteristics;
(2) the method can make up the deficiency of the hardware of the seeker to a certain extent, for example, the target and the interference which are caused by low resolution cannot be identified;
(3) compared with a common anti-interference guidance method such as tracking a target and interfering a mass center, the method has higher accuracy and can improve the hit rate of the missile.
Drawings
FIG. 1 is a guidance flow chart of anti-interference guidance of an infrared air-to-air missile;
FIG. 2 is a two-dimensional planar model of missile versus target engagement;
FIG. 3 is a schematic illustration of an acceleration probability distribution of an object;
FIG. 4 is a schematic diagram of a missile accessible zone and a target likely-to-occur zone;
FIG. 5 shows the time t when the target is thrown into the interfering bombgoBullet trajectory diagram of 4.5s
FIG. 6 is a comparison of the miss distance of the missile when the target throws the missile at different moments and the miss distance of the missile under different guidance laws.
FIG. 7 shows the target for the jamming of a bomb tgo2s bullet trajectory diagram
FIG. 8 shows the different remaining flight times t after the target has thrown the chaffgoAnd comparing the miss distance with the miss distance of the missile under different guidance laws.
Detailed Description
The invention is further described below with reference to the accompanying drawings.
The invention has the following implementation steps:
1. first, the engagement model in FIG. 2 is established
Taking a two-dimensional plane pursuit model as a research object, and selecting x (t) ═ r as a state quantityx(t) ry(t) vx(t) vy(t) atx(t) aty(t)]TWherein r (t) represents the relative distance of the bullet eyes, v (t) represents the approaching speed of the bullet eyesDegree of at(t) represents a target acceleration. The motion model of the target adopts a 'current' statistical model to establish a dynamic model of the pursuit system:
wherein
α represents the inverse of the maneuver time constant, i.e., the maneuver frequency.
u(t)=[amx(t) amy(t)]T
u (t) is the acceleration of the missile;
C=[0 0 0 0 α α]T
w(t)=[0 0 0 0 wtx wty]T
where w (t) is a mean of 0 and a variance ofWhite noise of (2). WhereinσaRespectively, the mean and variance of the target maneuvering acceleration in the "current" statistical model.
The measurement information only comprises the line-of-sight angle information, and the observation equation is as follows:
z(t)=h(x)+v(t)
obtained after discretization
Equation of observation
z(k)=h[x(k)]+v(k)
v (k) is measurement noise, which is a gaussian white random vector sequence.
(2) And performing state estimation on the target by using a particle filter estimation algorithm, wherein the number of particles N is 1000, and obtaining the state of the target, including information such as position, speed, acceleration and the like.
(3) Only the distribution of the targets in the y direction is considered at present, because the targets can be hit only if the missile and the targets are guaranteed to be at the same height. Let it be assumed here that the current acceleration of the target in the y-direction is aE(tk) Next moment target acceleration aE(tk+1) Mainly distributed in the field, and the distribution follows gamma distribution
Wherein a (b-1) ═ aE(tk)+aEmaxAnd the current acceleration satisfies aE(tk)∈[-aEmax,0]For aE(tk)∈[0,aEmax]In the case of (2), the probability density function is symmetrical to the above-mentioned interval along the ordinate. Wherein the acceleration corresponding to the maximum probability is the current acceleration aE(tk) And has aE(tk)=ab+b-8。
The probability density distribution of the predicted position of the target at the next time instant is then described as follows
Pr[yE(tk+1|tk)]=f(yE(tk),Pr[aE(tk+1)],tk)
Fig. 3 shows the probability distribution for different accelerations and different values of a, b.
(4) Missile tk+1Time of day prediction statetfThe missile at all possible states is integratedAll sets reachable in the y-direction are D (t)f,xP(tk+1|tk) The maximum reachable area of the missile is the distribution of a flight envelope formed by the missile flying at the maximum acceleration capacity of 50g in the y direction. FIG. 4 is a schematic view of the flight reach of a missile with interference.
(5) Proposing a performance function according to the interception requirement
The meaning of the performance function is the probability integral of the true and false targets contained in the reachable region of the missile, and the maximum probability is close to the target when the performance function is guaranteed to be maximum.
The optimal position of the missile at the next moment is
(6) The missile control command is
u(tk)=uy(tk)/cosθP
θP=arctan(vPy/vPx)
Dy=[0 1 0 0 0 0]
Dvy=[0 0 0 1 0 0]
The initial simulation conditions are as follows: the x and y directions of the initial position of the target are x respectivelyE(t0)=8000m,yE(t0) 6000m, initial velocity x, y directions are vEx(t0)=400m/s,vEy(t0) 0. The initial position x and y directions of the missile are xP=0,yP6000m, initial speed x, y direction xPx=1000m/s,v Py0. Wherein the burning time of the interference bomb is 3s, and figure 5 shows the target interference bomb throwing time tgoThe projectile trajectory diagram of 4.5s, and fig. 6 is the miss distance comparison between the target throwing the interfering projectile at different time and the miss distance of the missile under different guidance laws. It can be seen that the missile miss distance based on the method has better performance relative to the guidance method of tracking energy center, and fig. 7 shows that t is t when a target throws a plurality of continuous interference missilesgoFig. 8 shows different remaining flight times t after the target throws a plurality of consecutive interfering bombsgoAnd comparing the miss distance with the miss distance of the missile under different guidance laws.
Claims (1)
1. An anti-interference guidance method of an infrared air-to-air missile based on target terminal reachable region prediction is characterized by comprising the following steps:
(1) the infrared air-to-air missile seeker obtains measurement information of a target, filter estimation is carried out on target state information, and the residual flight time of the missile is estimated;
(2) assuming that the target acceleration follows the probability distribution of gamma distribution, and calculating the position distribution region of the target at the next moment and the meeting moment;
(3) calculating the reachable area of the missile at the next moment, and calculating the final moment of the missile at each possible position at the next moment with the maximum flight capacityAn area that can be reached;
(4) providing a missile control instruction which is calculated by taking a maximum probability hit target of the missile as a performance index and maximizing the performance;
(5) And (3) applying control to the missile, and returning to the step (1) for circulation until the missile meets the target.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710018102.XA CN107917646B (en) | 2017-01-10 | 2017-01-10 | Infrared air-to-air missile anti-interference guidance method based on target terminal reachable area prediction |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710018102.XA CN107917646B (en) | 2017-01-10 | 2017-01-10 | Infrared air-to-air missile anti-interference guidance method based on target terminal reachable area prediction |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107917646A CN107917646A (en) | 2018-04-17 |
CN107917646B true CN107917646B (en) | 2020-11-24 |
Family
ID=61898551
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710018102.XA Active CN107917646B (en) | 2017-01-10 | 2017-01-10 | Infrared air-to-air missile anti-interference guidance method based on target terminal reachable area prediction |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107917646B (en) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109255393A (en) * | 2018-09-30 | 2019-01-22 | 上海机电工程研究所 | Infrared Imaging Seeker anti-jamming performance evaluation method based on random forest |
CN109359395B (en) * | 2018-10-23 | 2021-08-27 | 北京理工大学 | Infrared guidance missile efficiency evaluation method and system based on Petri network |
CN112256055B (en) * | 2020-10-12 | 2021-11-05 | 清华大学 | Three-body confrontation defense prediction guidance method adopting fixed point optimization |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101975575A (en) * | 2010-10-15 | 2011-02-16 | 西安电子科技大学 | Multi-target tracking method for passive sensor based on particle filtering |
CN102393913A (en) * | 2011-10-31 | 2012-03-28 | 北京航空航天大学 | Accurate dim and small target tracking method based on spectral fingerprint characteristics |
CN102842047A (en) * | 2012-09-10 | 2012-12-26 | 重庆大学 | Infrared small and weak target detection method based on multi-scale sparse dictionary |
CN103884237A (en) * | 2014-04-08 | 2014-06-25 | 哈尔滨工业大学 | Several-for-one collaborative guidance method based on target probability distribution information |
CN104200495A (en) * | 2014-09-25 | 2014-12-10 | 重庆信科设计有限公司 | Multi-target tracking method in video surveillance |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP3074713B1 (en) * | 2013-11-25 | 2020-01-08 | BAE Systems PLC | System integration |
-
2017
- 2017-01-10 CN CN201710018102.XA patent/CN107917646B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101975575A (en) * | 2010-10-15 | 2011-02-16 | 西安电子科技大学 | Multi-target tracking method for passive sensor based on particle filtering |
CN102393913A (en) * | 2011-10-31 | 2012-03-28 | 北京航空航天大学 | Accurate dim and small target tracking method based on spectral fingerprint characteristics |
CN102842047A (en) * | 2012-09-10 | 2012-12-26 | 重庆大学 | Infrared small and weak target detection method based on multi-scale sparse dictionary |
CN103884237A (en) * | 2014-04-08 | 2014-06-25 | 哈尔滨工业大学 | Several-for-one collaborative guidance method based on target probability distribution information |
CN104200495A (en) * | 2014-09-25 | 2014-12-10 | 重庆信科设计有限公司 | Multi-target tracking method in video surveillance |
Also Published As
Publication number | Publication date |
---|---|
CN107917646A (en) | 2018-04-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR101411699B1 (en) | Systems and methods for disambiguating shooter locations with shockwave-only location | |
CN107917646B (en) | Infrared air-to-air missile anti-interference guidance method based on target terminal reachable area prediction | |
CN111025249B (en) | Method for judging interference effect of foil strip interference in simulation environment | |
US9689965B2 (en) | Hostile fire detection for an airborne platform | |
US20200284566A1 (en) | Search and protect device for airborne targets | |
CN111598232B (en) | Method for estimating complex jog space cone target parameters by using deep learning convolutional neural network | |
Huang et al. | Defense strategy of aircraft confronted with IR guided missile | |
CN114020021A (en) | Design method and system of multi-missile distributed cooperative guidance law | |
CN111880172A (en) | Method for measuring descending section trajectory by staring radar | |
CN114442666B (en) | Heterogeneous fine guidance projectile group collaborative damage assessment method | |
Fu et al. | Air combat assignment problem based on bayesian optimization algorithm | |
CN116294811A (en) | Armed helicopter fire flight cooperative mechanism analysis method based on multi-target wolf swarm algorithm | |
Cheng et al. | New-generation hybrid guidance system based on infrared and millimeter waves | |
CN114819055A (en) | Missile interception point prediction method based on LSTM network | |
You et al. | Air combat command and guidance situation assessment based on attack area | |
Jia et al. | Mode design and control structure of manned/unmanned aerial vehicles cooperative engagement | |
Liu et al. | Research on penetration technology of intelligent cluster missile system | |
Weilin et al. | Decision-making of one-on-one beyond-visual-range air combat based on improved q-network | |
Bao et al. | The influence of dual stealth aircraft formation cooperating with aerial bombing on the performance of monopulse radar | |
RU2771865C1 (en) | Method and device for multifactor protection of objects from miniature unmanned aerial vehicles | |
Zhang et al. | A Firepower Assignment Method for Barragging the JDAM | |
Tao et al. | Modeling study on the operational process of photoelectric-guided munitions | |
KR102252186B1 (en) | Apparatus for target selection of guided air vehicle | |
Tupitsyn | Estimation of the required dimension of net to capture drone | |
RU42302U1 (en) | MISSILE DEFENSE SYSTEM |
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 |