KR20160070573A - Real-time prediction method of impact point of guided missile - Google Patents

Real-time prediction method of impact point of guided missile Download PDF

Info

Publication number
KR20160070573A
KR20160070573A KR1020140177705A KR20140177705A KR20160070573A KR 20160070573 A KR20160070573 A KR 20160070573A KR 1020140177705 A KR1020140177705 A KR 1020140177705A KR 20140177705 A KR20140177705 A KR 20140177705A KR 20160070573 A KR20160070573 A KR 20160070573A
Authority
KR
South Korea
Prior art keywords
missile
flight
solution
intercept point
speed
Prior art date
Application number
KR1020140177705A
Other languages
Korean (ko)
Other versions
KR101658464B1 (en
Inventor
김성준
이용인
문관영
Original Assignee
국방과학연구소
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 국방과학연구소 filed Critical 국방과학연구소
Priority to KR1020140177705A priority Critical patent/KR101658464B1/en
Publication of KR20160070573A publication Critical patent/KR20160070573A/en
Application granted granted Critical
Publication of KR101658464B1 publication Critical patent/KR101658464B1/en

Links

Images

Classifications

    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F42AMMUNITION; BLASTING
    • F42BEXPLOSIVE CHARGES, e.g. FOR BLASTING, FIREWORKS, AMMUNITION
    • F42B15/00Self-propelled projectiles or missiles, e.g. rockets; Guided missiles
    • F42B15/01Arrangements thereon for guidance or control
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F41WEAPONS
    • F41GWEAPON SIGHTS; AIMING
    • F41G9/00Systems for controlling missiles or projectiles, not provided for elsewhere
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F41WEAPONS
    • F41HARMOUR; ARMOURED TURRETS; ARMOURED OR ARMED VEHICLES; MEANS OF ATTACK OR DEFENCE, e.g. CAMOUFLAGE, IN GENERAL
    • F41H13/00Means of attack or defence not otherwise provided for

Landscapes

  • Engineering & Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Chemical & Material Sciences (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Combustion & Propulsion (AREA)
  • Aiming, Guidance, Guns With A Light Source, Armor, Camouflage, And Targets (AREA)

Abstract

The present invention relates to a method to predict an impact point of a guided missile in real time. The method includes: a step of deriving a flight distance in consideration of the trace of a guided missile when inputting a flight condition; a step of analytically deriving a flight distance in consideration of a flight speed of the guided missile by linearizing the flight speed of the guided missile by section and integrating the speed for a residual flight time; a step of deriving an impact point interaction formula by combining the flight distance, in consideration of the flight speed, and the flight distance in consideration of the trace of the guided missile; a step of calculating a value of the derived impact point interaction formula; and a step of predicting an expected impact point and a residual flight time by using the calculated value of the impact point interaction formula. Therefore, the present invention is capable of enabling effective midflight.

Description

유도탄의 요격지점 실시간 예측 방법{REAL-TIME PREDICTION METHOD OF IMPACT POINT OF GUIDED MISSILE}REAL-TIME PREDICTION METHOD OF IMPACTION POINT OF GUIDED MISSILE

본 발명은 유도탄의 요격지점(Impact Point)을 실시간으로 예측할 수 있는 유도탄의 요격지점 실시간 예측 방법에 관한 것이다. The present invention relates to a method for real-time prediction of an intercept point of a missile capable of predicting an impact point of a missile in real time.

일반적으로 유도탄은 유도장치에 의해 비행하여 목표지점까지 유도되는 무기를 말하며, 요격 성능을 향상시키기 위해서는 요격지점 (Impact Point, IP)을 정확히 예측하는 것이 필요하다. In general, the missile is a weapon that is guided by the guiding device to the target point. In order to improve the interception performance, it is necessary to accurately predict the impact point (IP).

종래에는 유도탄의 표적 요격지점(IP)을 예측하는 기술로서 반복적인 계산 방법이 사용되었다. 이 방법은 반복적 계산에 기반을 두고 있어 정확한 해를 구하는데 많은 시간이 소요되기 때문에 비행 중 실시간 계산 시에는 반복 계산 횟수에 제한을 두고 있다.Conventionally, repetitive calculation methods have been used as a technique for predicting the target intercept point (IP) of a guided missile. Since this method is based on iterative calculation, it takes a lot of time to find the exact solution.

또한, 종래의 방법은 유한 반복횟수 이내 해의 수렴성을 보장할 수 없어 예측오차가 큰 단점이 있다. 따라서, 실시간성과 정확성이 보장되는 해석적 방법 기반의 실시간 요격지점 예측 방법의 필요성이 대두되었다.In addition, the conventional method can not guarantee convergence within a finite number of iterations, resulting in a large prediction error. Therefore, the necessity of real - time and accurate real - time interceptor location prediction method based on analytical method has emerged.

본 반명의 목적은 실시간성과 정확성이 보장되는 해석적 방법 기반의 요격지점 예측 방법을 제공하는데 있다. The purpose of this half is to provide an analytical method based interceptor location prediction method that guarantees real - time and accuracy.

본 발명의 다른 목적은 원활한 종말 호밍을 위해 초/중기 유도단계에서 요격예상지점을 산출하여 비행하는데 있다. Another object of the present invention is to calculate and fly an intercept prediction point in the early / mid-stage induction stage for smooth end homing.

상기와 같은 목적을 달성하기 위하여 본 발명의 실시예에 따른 유도탄의 요격지점 실시간 예측 방법은, 비행조건 입력시 유도탄의 비행궤적을 고려한 비행거리를 도출하는 단계; 유도탄의 비행속도를 구간 선형화하고 잔여비행시간 동안 적분하여 유도탄의 비행속도를 고려한 비행거리를 해석적으로 도출하는 단계; 상기 도출된 유도탄의 비행궤적을 고려한 비행거리와 유도탄의 비행속도를 고려한 비행거리를 연립하여 요격지점 관계식을 도출하는 단계; 상기 도출된 요격지점 관계식의 해를 산출하는 단계; 및 상기 산출된 요격지점 관계식의 해를 이용하여 요격예상지점 및 잔여비행시간을 실시간 예측하는 단계;를 포함한다. According to another aspect of the present invention, there is provided a method for real-time prediction of an intercept point of a guided missile in accordance with an embodiment of the present invention, comprising the steps of: The step of linearizing the flight speed of the missile and integrating it during the remaining flight time to analytically derive the flight distance considering the flight speed of the missile; Deriving an intercept point relationship based on the derived distance in consideration of the trajectory of the guided missile and the flight distance in consideration of the speed of the missile; Calculating a solution of the derived intercept point relationship; And estimating an intercept prediction point and a residual flight time in real time using the calculated solution of the intercept point relationship.

상기 요격지점 관계식은 무차원 요격지점에 대한 3차 방정식이다.The interceptor point relation is a cubic equation for the non-dimensional intercept point.

상기 요격지점 관계식의 해를 산출하는 단계는, 유도탄의 속도가 등속인지 판단하는 단계; 유도탄의 속도가 등속인 경우 3차 방정식의 요격지점 관계식을 2차 방정식으로 변환한 후 근의 공식을 이용하여 요격지점의 관계식의 해를 구하는 단계;를 포함할 수 있다.The step of calculating the solution of the intercept point relationship may include the steps of: determining whether the speed of the guided vehicle is constant; And if the velocity of the missile is constant, converting the intercept point relationship of the cubic equation into a quadratic equation and finding the solution of the relational expression of the intercept point using the root formula.

상기 요격지점 관계식의 해를 산출하는 단계는 유도탄의 속도가 등속이 아닌 경우 3차 방정식의 요격지점 관계식의 판별식을 산출하는 단계; 및 상기 산출된 판별식의 크기에 따라 하나의 해 또는 다중해를 산출하는 단계; 상기 산출된 해에서 다중해 선택논리를 만족하는 해를 선택하는 단계; 및 상기 선택된 해 중에서 소정 평가지수를 만족하는 해를 요격지점에 대한 해로 최종 결정하는 단계;를 포함할 수 있다.Calculating the solution of the interceptor point relational equation includes calculating a discriminant of the intercept point relation of the cubic equation when the velocity of the guided vehicle is not constant; And calculating one solution or multiple solutions according to the size of the calculated discrimination equation; Selecting a solution satisfying the multi-solution selection logic in the solution solution; And finally determining a solution satisfying a predetermined evaluation index from the selected solution as a solution to the intercept point.

상기 다중해 선택논리는 요격지점이 유도탄과 표적 사이에 위치하고, 유도탄 비행속도를 고려한 비행거리는 0 이상이며, 유도탄의 비행속도를 고려한 비행거리와 비행궤적을 고려한 비행거리는 유사하다는 조건을 포함할 수 있다.The multi-solution selection logic may include a condition that the intercept point is located between the guided missile and the target, the flying distance considering the missile flying velocity is 0 or more, and the flying distance considering the flying speed of the missile is similar to the flying distance considering the flying trajectory .

상기 구해진 해를 이용하여 유도탄의 경로점을 산출하는 단계;를 더 포함할 수 있다. And calculating a route point of the guided vehicle using the obtained solution.

본 발명은 비행조건이 주어지면 시뮬레이션을 수행할 필요 없이 즉시 요격지점을 예측할 수 있는 수식을 도출한 후 이를 이용하여 초/중기 유도단계에서 미리 요격예상지점을 실시간으로 산출함으로써 종말호밍 단계에서 표적 요격성능을 향상시키기 위한 효과적인 중기비행이 가능한 효과가 있다. In the present invention, when a flight condition is given, a formula for predicting an intercept point can be derived immediately without having to perform a simulation, and then an intercept prediction point is calculated in real time at a mid / mid- There is an effective mid-flight effect to improve performance.

또한 본 발명은 실시간 구현이 매우 용이한 요격지점 예측 방법을 유도탄 탑재컴퓨터에 구현함으로써 종말호밍 단계에서 표적 요격성능을 향상시키기 위한 효과적인 중기비행이 가능한 효과가 있다. In addition, the present invention realizes an effective mid-flight flight for improving the target interception performance in the end-homing step by implementing an intercept point prediction method which is very easy to real-time realization in a missile loading computer.

도 1은 도 1은 유도탄의 유도기하를 나타낸 그래프.
도 2는 유도탄 속도선도를 나타낸 그래프.
도 3은 본 발명에 따른 요격지점 실시간 예측 방법을 나타낸 순서도.
FIG. 1 is a graph showing the induction geometry of a guided vehicle; FIG.
Fig. 2 is a graph showing a guide shaft speed graph. Fig.
3 is a flowchart showing an intercept point real-time prediction method according to the present invention.

본 발명은 유도탄의 표적 요격지점(Impact Point, IP)을 예측하는 기술로서 종래의 반복적 계산에 기반한 방법의 어려움을 극복하기 위하여 해석적 방법 기반의 실시간 요격지점 예측 방법을 제공한다. 즉, 본 발명은 비행조건이 주어지면 시뮬레이션을 수행할 필요 없이 즉시 요격지점을 예측할 수 있는 수식을 도출한 후 이를 탑재 컴퓨터에 구현하여 초/중기 유도단계에서 미리 요격예상지점을 실시간 산출할 수 있는 알고리즘을 제공한다.The present invention provides a method for predicting a target intercept point (IP) of a missile, and a real-time intercept point prediction method based on an analytical method in order to overcome the difficulty of the method based on the conventional iterative calculation. That is, according to the present invention, when a flight condition is given, a formula for predicting an intercept point can be immediately derived without performing a simulation, and then the calculated formula can be implemented in a mounted computer, and an intercept prediction point can be calculated in real- Algorithm.

일반적으로 원활한 종말 호밍을 위해 초/중기 유도단계에서 요격예상지점을 산출한 후 해당 지점을 향해 유도탄이 비행하도록 하는 것이 중요하다. 이를 위해서는 발사 초기부터 요격지점을 정확히 예측하는 것이 필요하다. Generally, it is important to calculate the predicted point of interception during the early / mid-stage induction phase so that guided missiles fly toward the point. To do this, it is necessary to accurately predict the intercept point from the beginning of the launch.

요격지점에 대한 해석해를 이끌어내기 위하여 몇 가지 가정이 필요하다. Several assumptions are needed to elucidate the interpretation of the intercept point.

첫째, 표적은 유도탄과 동일 평면에서 기동하며 등고도 및 등속으로 유도탄 발사점을 향해 비행한다(도 1). First, the target is maneuvered in the same plane as the guided car, flying at the uniform height and constant velocity toward the guided firing point (Fig. 1).

둘째, 유도탄의 속도 선도(profile)는 부스팅(boosting) 단계와 글라이딩 (gliding) 단계로 나누고, 각 단계에서 시간에 대해 선형으로 변화하는 것으로 간주한다(도 2). Second, the velocity profile of the missile is divided into a boosting stage and a gliding stage, and is considered to change linearly with respect to time at each stage (FIG. 2).

셋째, 유도탄은 중기유도단계에서 거리에 대한 3차 다항식 형태의 궤적으로 비행하고 종말 호밍단계에서는 직선 비행한다. Third, the guided missiles fly in the third - order polynomial trajectory in the mid - stage induction phase and fly straight at the end homing stage.

이와 같은 가정 하에 본 발명은, 유도탄의 비행궤적을 고려한 비행거리의 해석해와 유도탄의 속도 선도를 적분하여 구한 비행거리의 해석해를 연립하고, 유도탄의 비행시간과 표적의 비행시간이 동일하다는 전제 하에 요격지점에 대한 3차 방정식을 도출한다. 이후 도출된 3차 방정식을 해를 구하여 요격지점에 대한 해석해를 구하고 이를 탑재 컴퓨터에 구현하여 요격예상지점 및 잔여비행시간을 실시간 산출한다.Under this assumption, the present invention is based on the assumption that the interpretation solution of the flight distance in consideration of the flight path of the missile is integrated with the velocity diagram of the missile, Derive a cubic equation for the point. Then, the solution of the cubic equation is derived, and the interpretation solution for the intercept point is obtained. Then, it is implemented in the onboard computer to calculate the intercept expectation point and the remaining flight time in real time.

이후 본 발명의 실시예에 따른 유도탄의 요격지점 실시간 예측 방법을 보다 상세히 설명하면 다음과 같다. Hereinafter, a method for real-time prediction of an intercept point of a missile according to an embodiment of the present invention will be described in detail.

1. 유도탄 비행궤적을 고려한 비행거리 1. Flight distance considering missile flight trajectory

도 1은 유도탄의 유도기하를 나타낸 그래프이다. 1 is a graph showing the induction geometry of a guided vehicle.

도 1에서 경로점(WP)까지 비행하는 중기유도구간에서 비행궤적은 거리의 3차 다항식 형태라고 가정하면, 중기유도단계의 곡선거리(S1)는 다음 수학식 1과 같이 나타낼 수 있다. In FIG. 1, assuming that the flight trajectory is a third-order polynomial shape of the distance in the midway induction section flying to the point WP, the curve distance S 1 in the middle stage induction phase can be expressed by the following Equation 1:

[수학식1] [Equation 1]

Figure pat00001
Figure pat00001

여기서,

Figure pat00002
는 유도탄 기준 요격지점(IP)의 수평축 위치이고, d는 요격지점 (IP)과 경로점(WP)의 상대거리이다. here,
Figure pat00002
Is the horizontal axis position of the guided missile reference intercept point (IP), and d is the relative distance between the intercept point (IP) and the route point (WP).

Figure pat00003
이다.
Figure pat00003
to be.

여기서,

Figure pat00004
는 현 시점의 유도탄의 비행 경로각,
Figure pat00005
은 유도탄의 경로점 통과각,
Figure pat00006
는 현 시점의 표적 기준 유도탄의 수직축 위치이다. here,
Figure pat00004
Is the flight path angle of the present guided missile,
Figure pat00005
Is the path passing angle of the missile,
Figure pat00006
Is the vertical axis position of the target guided missile at the present time.

따라서, 요격지점(IP)까지의 비행거리(SR)는 다음 수학식 2와 같이 나타낼 수 있다.Therefore, the flying distance S R to the intercepting point IP can be expressed by the following equation (2).

[수학식 2]&Quot; (2) "

Figure pat00007
Figure pat00007

2. 유도탄 비행속도를 고려한 비행거리2. Flight distance considering missile flight speed

도 2는 유도탄의 속도선도를 도시한 그래프이다. 2 is a graph showing a speed diagram of a guided vehicle.

도 2를 참조하면, 유도탄의 속도선도는 부스팅 단계와 글라이딩 단계로 나누어 지며, 이를 수식으로 표현하면 수학식 3과 같다.Referring to FIG. 2, the speed diagram of the missile is divided into a boosting stage and a gliding stage.

[수학식 3]&Quot; (3) "

Figure pat00008
Figure pat00008

여기서

Figure pat00009
는 부스팅 단계 종료시점,
Figure pat00010
는 부스팅 단계 종료시점의 유도탄 속력으로
Figure pat00011
이며,
Figure pat00012
는 부스팅 단계의 평균 가속도이고,
Figure pat00013
는 글라이딩 단계의 평균 가속도이다. . here
Figure pat00009
At the end of the boosting phase,
Figure pat00010
Is the speed of the missile at the end of the boosting phase
Figure pat00011
Lt;
Figure pat00012
Is the average acceleration of the boosting phase,
Figure pat00013
Is the average acceleration of the gliding phase. .

따라서, 현 시점부터 요격시점까지 잔여비행시간을 tgo라고 할 때, 유도탄 속도를 적분한 비행거리 즉, 유도탄 비행속도를 고려한 비행거리(SL)는 다음 수학식 4와 같다. Therefore, when the remaining flight time from the current point to the intercept point is t go , the flying distance (S L ), which takes the missile speed into account, ie, the missile flying speed, is expressed by the following equation (4).

[수학식 4]&Quot; (4) "

Figure pat00014
Figure pat00014

여기서,

Figure pat00015
는 현시점의 유도탄 속력이고, here,
Figure pat00015
Is the current speed of the touring car,

Figure pat00016
이다.
Figure pat00016
to be.

3. 잔여비행시간 계산3. Calculate remaining flight time

도 1에서 유도탄이 표적을 요격하는데 걸리는 시간은 표적이 요격지점까지 이동하는 시간과 동일하므로 잔여비행시간(tgo)은 수학식 5와 같이 나타낼 수 있다.In FIG. 1, the time taken for the missile to intercept the target is equal to the time required for the target to travel to the intercept point, so that the remaining flight time t go can be expressed by Equation (5).

[수학식 5] &Quot; (5) "

Figure pat00017

Figure pat00017

4. 요격지점에 관한 관계식4. Relation on Intercept Point

수학식 2의 유도탄 비행궤적을 고려한 비행거리(SR)와 수학식 4의 유도탄 비행속도를 고려한 비행거리(SL)는 동일해야 하며, 수학식 4의 잔여비행시간(tgo)은 수학식 5와 같으므로, 이를 모두 만족시키는 식은 다음의 수학식 6과 같이 무차원 요격지점(

Figure pat00018
)에 대한 3차 방정식으로 표현된다.And will be the same mathematical fly Considering missile flight trajectory of the formula (2) Flight distance (S L) taking into account the missile flight speed of (S R) and the equation (4), and the remaining flight time (t go) of the equation (4) is Equation 5, the equation that satisfies all of them is expressed by the following equation (6)
Figure pat00018
) Is expressed by a cubic equation.

[수학식 6]&Quot; (6) "

Figure pat00019
Figure pat00019

여기서, here,

Figure pat00020
이고,
Figure pat00020
ego,

Figure pat00021
Figure pat00021

Figure pat00022
이며,
Figure pat00022
Lt;

Figure pat00023
이다.
Figure pat00023
to be.

5. 요격지점에 관한 해5. A description of the intercept point

5-1. 유도탄의 속도가 일정한 경우5-1. If the speed of the missile is constant

유도탄 속도가 일정한 경우

Figure pat00024
이므로 수학식 6은 수학식 7과 같이 나타낼 수 있다. Constant speed of missile
Figure pat00024
Equation (6) can be expressed by Equation (7).

[수학식 7]&Quot; (7) "

Figure pat00025
Figure pat00025

여기서, here,

Figure pat00026
이다.
Figure pat00026
to be.

따라서, 2차 방정식의 근의 공식으로부터 다음과 같이 해를 구할 수 있다.Therefore, from the root equation of the quadratic equation, the solution can be obtained as follows.

Figure pat00027
인 경우①
Figure pat00027
If

이 경우 두 실근이 존재하는데

Figure pat00028
인 경우
Figure pat00029
이어야 하므로 그의 해인 무차원 요격지점(
Figure pat00030
)은 다음의 수학식 8과 같다.In this case, there are two exoskeletons
Figure pat00028
If
Figure pat00029
So it's his year, the dimensionless intercept point (
Figure pat00030
) Is expressed by the following equation (8).

[수학식 8]&Quot; (8) "

Figure pat00031
Figure pat00031

Figure pat00032
인 경우②
Figure pat00032
If

이 경우 해인 무차원 요격지점(

Figure pat00033
)은 다음의 수학식 9와 같다.In this case, a non-dimensional intercept point (
Figure pat00033
) ≪ / RTI >

[수학식 9]&Quot; (9) "

Figure pat00034
Figure pat00034

③ 그 외의 경우③ Other cases

이 경우 해가 존재하지 않으므로 실제 구현시에는 아래 수학식 10과 같이 유도탄이 요격지점(IP)을 향해 직선비행한다고 가정한 결과로 대체한다. In this case, since there is no solution, the actual implementation is replaced with a result that assumes that the missile will fly straight to the intercept point (IP) as shown in Equation 10 below.

[수학식 10]&Quot; (10) "

Figure pat00035
Figure pat00035

5-2. 유도탄의 속도가 변하는 경우5-2. When the speed of the missile is changed

유도탄의 속도가 변하는 경우

Figure pat00036
이므로 수학식 6을 다음의 수학식 11과 같이 다시 쓸 수 있다.When the speed of the missile is changed
Figure pat00036
The equation (6) can be rewritten as the following equation (11).

[수학식 11]&Quot; (11) "

Figure pat00037
Figure pat00037

여기서,here,

Figure pat00038
이다.
Figure pat00038
to be.

3차 방정식의 근은 판별식에 따라 근의 형태가 다르게 나타나므로 먼저 판별식을 다음 수학식 12와 같이 산출한다.Since the roots of the cubic equation appear differently in the shape of the roots according to the discriminant, first, the discriminant is calculated according to the following equation (12).

[수학식 12]&Quot; (12) "

Figure pat00039
Figure pat00039

여기서, here,

Figure pat00040
이다.
Figure pat00040
to be.

따라서, 해는 다음과 같다. Therefore, the solution is as follows.

① D > 0인 경우① When D> 0

[수학식 13]&Quot; (13) "

Figure pat00041
Figure pat00041

여기서,

Figure pat00042
이다.here,
Figure pat00042
to be.

② D = 0인 경우② When D = 0

이 경우에는 중근을 포함한 세 실근이 존재하므로 아래의 수학식 14에 도시된 두 근 중 다중해 선택논리에 따라 적절한 값을 선택한다.In this case, since there are three real roots including the root, appropriate values are selected according to the multi-solution selection logic of the two roots shown in Equation (14) below.

[수학식 14] &Quot; (14) "

Figure pat00043
Figure pat00043

③ D < 0인 경우③ When D <0

이 경우 서로 다른 세 실근이 존재하므로 아래의 수학식 15에 도시된 세 근 중에서 다중해 선택논리에 따라 적절한 값을 선택한다. In this case, since there are three different real roots, an appropriate value is selected according to the multi-solution selection logic among the three roots shown in Equation (15) below.

[수학식 15]&Quot; (15) &quot;

Figure pat00044
Figure pat00044

Figure pat00045
Figure pat00045

여기서, here,

Figure pat00046
이다.
Figure pat00046
to be.

5-3. 다중해 선택 논리5-3. Multipath selection logic

3차 방정식인 수학식 11에서 판별식이 0 이하인 경우 산출되는 다중해로부터 아래의 물리적 관계식을 만족하는 해를 골라낸다.In Equation (11), which is a cubic equation, when the discriminant equation is 0 or less, a solution satisfying the following physical relation is selected from the multiplexed solution.

① 요격지점(IP)은 유도탄과 표적 사이에 위치한다. 즉,① The intercept point (IP) is located between the missile and the target. In other words,

Figure pat00047
Figure pat00047

② 유도탄 비행속도를 고려한 비행거리(SL)는 0 이상이다. 즉,② Flight distance ( SL ) considering missile flying speed is more than zero. In other words,

Figure pat00048
Figure pat00048

여기서,

Figure pat00049
이다.here,
Figure pat00049
to be.

③ 유도탄의 비행속도를 고려한 비행거리(SL)와 비행궤적을 고려한 비행거리(SR)는 유사하다. 즉, (3) The flying distance (S L ) considering the flight speed of the missile is similar to the flying distance (S R ) considering the flight trajectory. In other words,

Figure pat00050
Figure pat00050

여기서,

Figure pat00051
이고,
Figure pat00052
는 양의 작은 값이다. here,
Figure pat00051
ego,
Figure pat00052
Is a small positive value.

상기 ①~③ 조건을 모두 만족시키는 값 중에서 아래의 수학식 16과 같은 평가지수를 만족하는 해를 선택한다.Among the values satisfying all of the conditions (1) to (3), a solution satisfying the evaluation index as shown in the following equation (16) is selected.

[수학식 16]&Quot; (16) &quot;

Figure pat00053
Figure pat00053

여기서,

Figure pat00054
는 유도탄이 등속비행(
Figure pat00055
)인 경우의 해로서, 수학식 8~10에 의해 결정된 값이다.here,
Figure pat00054
(A)
Figure pat00055
), Which is a value determined by Equations (8) to (10).

따라서, 본 발명은 상기에서 구해진 해석해를 탑재 컴퓨터에 적용하여 하기와 같이 요격예상지점과 잔여비행시간을 실시간으로 산출한다. Accordingly, the present invention applies the interpretation solution obtained above to the on-board computer to calculate the intercept prediction point and the remaining flight time in real time as follows.

6. 요격지점(IP) 및 경로점(WP) 예측6. Prediction of Intercept Point (IP) and Path Point (WP)

무차원 요격지점(

Figure pat00056
)을 이용하여 다음의 수학식 17과 같이 요격예상지점을 산출한다.Non-dimensional intercept point (
Figure pat00056
) To calculate an intercept prediction point as shown in the following Equation (17).

[수학식 17]&Quot; (17) &quot;

Figure pat00057
Figure pat00057

여기서,

Figure pat00058
은 유도탄의 수평축 절대좌표,
Figure pat00059
는 표적의 수평 및 수직 절대좌표이다. 따라서, 경로점(WP)의 위치는 다음의 수학식 18과 같다.here,
Figure pat00058
The horizontal axis absolute coordinates of the missile,
Figure pat00059
Is the horizontal and vertical absolute coordinates of the target. Therefore, the position of the path point WP is expressed by the following equation (18).

[수학식 18]&Quot; (18) &quot;

Figure pat00060
Figure pat00060

7. 잔여비행시간 예측7. Predict remaining flight time

수학식 5를 이용하면 쉽게 잔여비행시간(tgo)을 산출할 수 있으나, 표적 속도가 0인 경우에는 계산이 불가능하므로 수학식 4를 이용하여 잔여비행시간(tgo)을 다음의 수학식 19와 같이 산출한다. Using the equation (5) easier, but can be used to calculate the remaining time of flight (t go), the target speed is 0, the calculation is only possible by using the equation (4) the remaining flight time (t go), and then equation (19) of the case .

[수학식 19]&Quot; (19) &quot;

Figure pat00061
Figure pat00061

여기서,

Figure pat00062
이다.here,
Figure pat00062
to be.

도 3은 본 발명의 실시예에 따른 유도탄의 요격지점 실시간 예측 방법을 나타낸 순서도이다. 3 is a flowchart illustrating a method for real-time prediction of an intercept point of a missile according to an embodiment of the present invention.

도 3에 도시된 바와 같이, 비행조건(유도탄 및 표적정보)이 입력되면(S100), 초/중기 유도단계에서 유도탄의 비행궤적을 고려한 비행거리를 해석적으로 도출하고(S110)(수학식 2), 유도탄의 비행속도를 구간 선형화하고 잔여비행시간 동안 적분하여 유도탄의 비행속도를 고려한 비행거리를 해석적으로 도출한다(수학식 4)(S120). 3, when the flight condition (missile and target information) is inputted (S100), the flight distance considering the flight path of the missile is calculated (S110) ), The flight speed of the missile is linearized and integrated during the remaining flight time to calculate the flight distance considering the flight speed of the missile (S120).

이어서, 상기 도출된 유도탄의 비행궤적을 고려한 비행거리와 유도탄의 비행속도를 고려한 비행거리를 연립하여 수학식 6과 같이 요격지점의 관계식을 도출한다(S130). Then, a relational expression of an intercept point is derived as shown in Equation (6) by taking the flight distance in consideration of the derived trajectory of the guided missile and the flight distance of the missile in consideration of the flight speed of the missile (S130).

그리고, 상기 도출된 요격지점의 관계식의 해를 구하여 요격지점에 대한 해석해를 구하는데, 상기 요격지점의 관계식의 해는 유도탄의 속도가 일정한 경우(등속)와 변하는 경우에 따라 달라진다. Then, the solution of the relational expression of the derived intercept point is obtained to obtain the interpretation solution for the intercept point. The solution of the relational expression of the intercept point varies depending on the case where the speed of the guided car is constant (constant speed).

따라서, 유도탄의 속도가 등속인지 체크하여(S140), 등속이 아닌 경우에는 3차방정식인 요격지점 관계식의 다중해를 구한 후 다중해 선택논리에 근거하여 하나의 최종해를 결정한다(S150, S160). 반면에 등속인 경우에는 2차방정식으로 변환된 요격지점 관계식의 해를 산출한다(S170). Accordingly, it is checked whether the velocity of the missile is constant (S140). If the velocity is not constant, a final solution is determined on the basis of the multiple solution selection logic after obtaining the multiple solution of the interceptor point relationship which is the cubic equation (S150, S160 ). On the other hand, if the velocity is constant, a solution of the intercept point relationship transformed into the quadratic equation is calculated (S170).

일단 요격지점 관계식의 해가 구해지면 해당 해를 이용하여 요격예상지점, 잔여비행시간 및 경로점을 실시간으로 산출한다(S180)(수학식 17~19). Once the solution of the intercept point relation is obtained, the intercept prediction point, the remaining flight time and the route point are calculated in real time using the corresponding solution (S180) (Equations 17 to 19).

상술한 바와 같이, 본 발명은 비행조건이 주어지면 시뮬레이션을 수행할 필요 없이 즉시 요격지점을 예측할 수 있는 수식을 도출한 후 이를 이용하여 초/중기 유도단계에서 미리 요격예상지점을 실시간으로 산출함으로써 종말호밍단계에서 표적 요격성능을 향상시키기 위한 효과적인 중기비행이 가능하다. As described above, according to the present invention, when a flight condition is given, a formula for predicting an intercept point can be immediately derived without performing a simulation, and then an intercept prediction point is calculated in real time at a mid / mid- Effective mid-flight flight is possible to improve target interception performance at the homing stage.

상기와 같이 설명된 유도탄의 요격지점 실시간 예측 방법은 상기 설명된 실시예들의 구성과 방법이 한정되게 적용될 수 있는 것이 아니라, 상기 실시예들은 다양한 변형이 이루어질 수 있도록 각 실시예들의 전부 또는 일부가 선택적으로 조합되어 구성될 수도 있다.The method of real-time prediction of the intercepting point of the guided missiles described above can be applied to the construction and method of the embodiments described above in a limited manner, but the embodiments can be applied to all or some of the embodiments As shown in FIG.

Claims (9)

비행조건 입력시 유도탄의 비행궤적을 고려한 비행거리를 도출하는 단계;
유도탄의 비행속도를 구간 선형화하고 잔여비행시간 동안 적분하여 유도탄의 비행속도를 고려한 비행거리를 해석적으로 도출하는 단계;
상기 도출된 유도탄의 비행궤적을 고려한 비행거리와 유도탄의 비행속도를 고려한 비행거리를 연립하여 요격지점 관계식을 도출하는 단계;
상기 도출된 요격지점 관계식의 해를 산출하는 단계; 및
상기 산출된 요격지점 관계식의 해를 이용하여 요격예상지점 및 잔여비행시간을 실시간 예측하는 단계;를 포함하는 것을 특징으로 하는 유도탄의 요격지점 실시간 예측 방법.
Deriving a flight distance in consideration of the flight path of the missile when entering the flight condition;
The step of linearizing the flight speed of the missile and integrating it during the remaining flight time to analytically derive the flight distance considering the flight speed of the missile;
Deriving an intercept point relationship based on the derived distance in consideration of the trajectory of the guided missile and the flight distance in consideration of the speed of the missile;
Calculating a solution of the derived intercept point relationship; And
And estimating an intercept prediction point and a residual flight time in real time using the solution of the calculated intercept point relationship.
제1항에 있어서, 상기 요격지점의 관계식은
무차원 요격지점에 대한 3차 방정식으로,
유도탄의 비행궤적을 고려한 비행거리와 유도탄의 비행속도를 고려한 비행거리는 동일하고, 유도탄이 표적을 요격하는데 걸리는 시간과 표적이 요격지점까지 이동하는 시간이 동일하다는 것을 모두 만족시키는 것을 특징으로 하는 유도탄의 요격지점 실시간 예측 방법.
The method according to claim 1, wherein the relational expression of the intercept point
As a cubic equation for a non-dimensional intercept point,
Wherein the flying distance in consideration of the flight path in consideration of the flight path of the missile and the flight speed of the missile are the same and the time required for the missile to intercept the target and the time for the target to travel to the interception point are the same Real time prediction method of intercept point.
제1항에 있어서, 상기 요격지점 관계식의 해를 산출하는 단계는
유도탄의 속도가 등속인지 판단하는 단계; 및
유도탄의 속도가 등속인 경우 3차 방정식의 요격지점 관계식을 2차 방정식으로 변환한 후 근의 공식을 이용하여 요격지점의 관계식의 해를 구하는 단계;를 포함하는 것을 특징으로 하는 유도탄의 요격지점 실시간 예측 방법.
2. The method of claim 1, wherein the step of calculating the solution of the intercept point relationship
Determining whether the speed of the guided vehicle is constant; And
And a step of finding a solution of a relational expression of an intercept point by using a formula of a root after transforming an intercept point relationship of a cubic equation into a quadratic equation when the speed of the guided vehicle is constant, Prediction method.
제3항에 있어서, 상기 요격지점 관계식의 해를 구하는 단계는
유도탄의 속도가 등속이 아닌 경우 3차 방정식의 요격지점 관계식의 판별식을 산출하는 단계; 및
상기 산출된 판별식의 크기에 따라 하나의 해 또는 다중해를 산출하는 단계;를 포함하는 것을 특징으로 하는 유도탄의 요격지점 실시간 예측 방법.
4. The method of claim 3, wherein the step of obtaining the solution of the intercept point relationship
Calculating a discriminant of the intercept point relation of the cubic equation when the speed of the guided vehicle is not constant; And
And calculating one solution or multiple solution according to the calculated size of the discrimination equation.
제4항에 있어서, 상기 산출된 해에서 다중해 선택논리를 만족하는 해를 선택하는 단계; 및
상기 선택된 해중에서 소정 평가지수를 만족하는 해를 요격지점에 대한 해로 최종 결정하는 단계;를 더 포함하는 것을 특징으로 하는 유도탄의 요격지점 실시간 예측 방법.
5. The method of claim 4, further comprising: selecting a solution satisfying the multi-solution selection logic in the solution solution; And
And finally determining a solution satisfying a predetermined evaluation index in the selected sea as a solution for an intercept point.
제5항에 있어서, 상기 다중해 선택논리는
요격지점이 유도탄과 표적 사이에 위치하고,
유도탄 비행속도를 고려한 비행거리는 0 이상이며,
유도탄의 비행속도를 고려한 비행거리와 비행궤적을 고려한 비행거리는 유사하다는 조건을 포함하는 것을 특징으로 하는 유도탄의 요격지점 실시간 예측 방법.
6. The method of claim 5,
An intercept point is located between the missile and the target,
The flying distance considering the missile flying speed is more than 0,
Wherein the estimated distance includes a condition that the distance of the missile in consideration of the speed of the missile is similar to the distance of the missile in consideration of the trajectory of the missile.
제1항에 있어서, 상기 요격예상지점(
Figure pat00063
,
Figure pat00064
)은
다음 식에 의해 계산하는 것을 특징으로 하는 유도탄의 요격지점 실시간 예측 방법.
Figure pat00065

여기서,
Figure pat00066
는 구해진 요격지점에 대한 해석해,
Figure pat00067
은 유도탄의 수평축 절대 좌표,
Figure pat00068
는 표적의 절대 좌표이다.
The method according to claim 1,
Figure pat00063
,
Figure pat00064
)silver
And calculating by using the following equation.
Figure pat00065

here,
Figure pat00066
Interpreted the obtained intercept point,
Figure pat00067
The horizontal axis absolute coordinates of the missile,
Figure pat00068
Is the absolute coordinate of the target.
제1항에 있어서, 상기 잔여비행시간은
다음 식에 의해 계산하는 것을 특징으로 하는 유도탄의 요격지점 실시간 예측 방법.
Figure pat00069

여기서,
Figure pat00070
이고,
Figure pat00071
Figure pat00072
Figure pat00073
,
Figure pat00074
은 현 시점의 유도탄의 속력,
Figure pat00075
는 글라이딩 단계의 유도탄의 평균 가속도,
Figure pat00076
는 표적의 수평축 위치이다.
The method of claim 1,
And calculating by using the following equation.
Figure pat00069

here,
Figure pat00070
ego,
Figure pat00071
Figure pat00072
Figure pat00073
,
Figure pat00074
The speed of the current missile,
Figure pat00075
The average acceleration of the missile in the gliding phase,
Figure pat00076
Is the horizontal axis position of the target.
제1항에 있어서, 상기 구해진 해를 이용하여 유도탄의 경로점을 산출하는 단계;를 더 포함하며, 상기 경로점(
Figure pat00077
,
Figure pat00078
)은 다음 식에 의해 산출하는 것을 특징으로 하는 유도탄의 요격지점 실시간 예측 방법.
Figure pat00079

여기서,
Figure pat00080
,
Figure pat00081
는 요격지점의 절대 위치, d는 요격지점과 경로점의 상대거리,
Figure pat00082
은 유도탄의 경로점 통과각이다.
2. The method according to claim 1, further comprising: calculating a route point of the guided car using the obtained solution,
Figure pat00077
,
Figure pat00078
) Is calculated by the following equation.
Figure pat00079

here,
Figure pat00080
,
Figure pat00081
Is the absolute position of the intercept point, d is the relative distance between the intercept point and the path point,
Figure pat00082
Is the path passing angle of the missile.
KR1020140177705A 2014-12-10 2014-12-10 Real-time prediction method of impact point of guided missile KR101658464B1 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020140177705A KR101658464B1 (en) 2014-12-10 2014-12-10 Real-time prediction method of impact point of guided missile

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020140177705A KR101658464B1 (en) 2014-12-10 2014-12-10 Real-time prediction method of impact point of guided missile

Publications (2)

Publication Number Publication Date
KR20160070573A true KR20160070573A (en) 2016-06-20
KR101658464B1 KR101658464B1 (en) 2016-09-22

Family

ID=56354390

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020140177705A KR101658464B1 (en) 2014-12-10 2014-12-10 Real-time prediction method of impact point of guided missile

Country Status (1)

Country Link
KR (1) KR101658464B1 (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20180026269A (en) * 2016-09-02 2018-03-12 국방과학연구소 Apparatus and Method for extracting available reaction time based on applying responsive domain having irregular shape to origin shift
CN107942673A (en) * 2017-12-11 2018-04-20 北京理工大学 A kind of martian atmosphere approach section parsing method of guidance of parachute-opening point height tracing
EP3415860A1 (en) * 2017-06-16 2018-12-19 Thales Method for predicting the trajectory of a hostile aircraft, particularly in the context of anti-air defence
KR102031929B1 (en) * 2019-04-02 2019-10-14 엘아이지넥스원 주식회사 Apparatus and method for terminal lead angle control with Time Varying Continuous Biased PNG
CN114754628A (en) * 2022-03-31 2022-07-15 南京理工大学 Flight body trajectory control method based on drop point prediction and virtual tracking

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007010273A (en) * 2005-07-01 2007-01-18 Toshiba Corp Meeting time calculation device, blasting unit of warhead and missile
JP2007178049A (en) * 2005-12-27 2007-07-12 Toshiba Corp Flying object guidance system and target meeting time estimating method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007010273A (en) * 2005-07-01 2007-01-18 Toshiba Corp Meeting time calculation device, blasting unit of warhead and missile
JP2007178049A (en) * 2005-12-27 2007-07-12 Toshiba Corp Flying object guidance system and target meeting time estimating method

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
논문자료(2010.05.27) *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20180026269A (en) * 2016-09-02 2018-03-12 국방과학연구소 Apparatus and Method for extracting available reaction time based on applying responsive domain having irregular shape to origin shift
EP3415860A1 (en) * 2017-06-16 2018-12-19 Thales Method for predicting the trajectory of a hostile aircraft, particularly in the context of anti-air defence
FR3067840A1 (en) * 2017-06-16 2018-12-21 Thales METHOD OF PREDICTING THE TRACK OF A HOSTILE AIRCRAFT IN PARTICULAR IN THE CONTEXT OF AN ANTIAERIAN DEFENSE
EP3415860B1 (en) 2017-06-16 2022-11-16 Thales Method for predicting the trajectory of a hostile aircraft, particularly in the context of anti-air defence
CN107942673A (en) * 2017-12-11 2018-04-20 北京理工大学 A kind of martian atmosphere approach section parsing method of guidance of parachute-opening point height tracing
KR102031929B1 (en) * 2019-04-02 2019-10-14 엘아이지넥스원 주식회사 Apparatus and method for terminal lead angle control with Time Varying Continuous Biased PNG
CN114754628A (en) * 2022-03-31 2022-07-15 南京理工大学 Flight body trajectory control method based on drop point prediction and virtual tracking
CN114754628B (en) * 2022-03-31 2023-08-04 南京理工大学 Flying body trajectory control method based on drop point prediction and virtual tracking

Also Published As

Publication number Publication date
KR101658464B1 (en) 2016-09-22

Similar Documents

Publication Publication Date Title
KR101658464B1 (en) Real-time prediction method of impact point of guided missile
CN106020231B (en) Hypersonic aircraft reentry trajectory optimization method based on reentry point parameter
CN108153323B (en) A kind of high-altitude unmanned vehicle high-precision reentry guidance method
CN103486905B (en) Determining method for terminal guidance shift-exchange conditions of reenter vehicle
EP3004786B1 (en) Method of fire control for gun-based anti-aircraft defence
CN111027206B (en) Self-adaptive sliding mode control method for intercepting maneuvering target with specified performance
CN103592946B (en) A kind of powered phase program self-adaptation longitudinal guidance method measured based on apparent acceleration
CN104077469A (en) Speed prediction based segmentation iteration remaining time estimation method
CN109033493B (en) Unscented Kalman filtering-based method for identifying pneumatic parameters of high-speed rotation bullet
CN109063242B (en) Guidance tool error identification method based on particle swarm optimization
CN108073742B (en) Method for estimating flight state of intercepted missile tail section based on improved particle filter algorithm
CN107367941B (en) Method for observing attack angle of hypersonic aircraft
CN104503471A (en) Terminal guidance method for maneuvering aircraft multi-terminal constraint backstepping sliding mode
CN109827478B (en) guidance method with falling angle constraint and overload constraint
Xu et al. Application of situation function in air combat differential games
CN111832121A (en) Multi-aircraft cooperative detection and guidance integrated method and system
EP3669248B1 (en) Gbias for rate based autopilot
KR101492830B1 (en) Apparatus and method for analyzing weapon delivery of aircraft
KR101874257B1 (en) In-flight Impact Angle Determination Apparatus for Biased Proportional Navigation Guidance, Method thereof, and Computer readable storage having the same method
Chen et al. Impact time and angle constrained guidance via range‐based line‐of‐sight shaping
CN104627388B (en) A kind of self adaptation ballistic prediction method of reentry vehicle
CN110728010A (en) Method for calculating firing data of antiaircraft gun fire control system
CN109977543A (en) Three-dimensional profile coverage area calculation method, system and medium based on lateral priority
KR102501275B1 (en) Precise trajectory prediction method and apparatus of a highly maneuvering target through flight dynamics estimation and multipath suppression
Huang et al. Radar tracking for hypersonic glide vehicle based on aerodynamic model

Legal Events

Date Code Title Description
E902 Notification of reason for refusal
E701 Decision to grant or registration of patent right