KR20160070573A - Real-time prediction method of impact point of guided missile - Google Patents
Real-time prediction method of impact point of guided missile Download PDFInfo
- 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
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
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F41—WEAPONS
- F41G—WEAPON SIGHTS; AIMING
- F41G9/00—Systems for controlling missiles or projectiles, not provided for elsewhere
-
- F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
- F41—WEAPONS
- F41H—ARMOUR; ARMOURED TURRETS; ARMOURED OR ARMED VEHICLES; MEANS OF ATTACK OR DEFENCE, e.g. CAMOUFLAGE, IN GENERAL
- F41H13/00—Means 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
Description
본 발명은 유도탄의 요격지점(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]
여기서, 는 유도탄 기준 요격지점(IP)의 수평축 위치이고, d는 요격지점 (IP)과 경로점(WP)의 상대거리이다. here, 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).
이다. to be.
여기서, 는 현 시점의 유도탄의 비행 경로각, 은 유도탄의 경로점 통과각, 는 현 시점의 표적 기준 유도탄의 수직축 위치이다. here, Is the flight path angle of the present guided missile, Is the path passing angle of the missile, 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) "
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) "
여기서 는 부스팅 단계 종료시점, 는 부스팅 단계 종료시점의 유도탄 속력으로 이며, 는 부스팅 단계의 평균 가속도이고, 는 글라이딩 단계의 평균 가속도이다. . here At the end of the boosting phase, Is the speed of the missile at the end of the boosting phase Lt; Is the average acceleration of the boosting phase, 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) "
여기서, 는 현시점의 유도탄 속력이고, here, Is the current speed of the touring car,
이다.
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) "
4. 요격지점에 관한 관계식4. Relation on Intercept Point
수학식 2의 유도탄 비행궤적을 고려한 비행거리(SR)와 수학식 4의 유도탄 비행속도를 고려한 비행거리(SL)는 동일해야 하며, 수학식 4의 잔여비행시간(tgo)은 수학식 5와 같으므로, 이를 모두 만족시키는 식은 다음의 수학식 6과 같이 무차원 요격지점()에 대한 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) ) Is expressed by a cubic equation.
[수학식 6]&Quot; (6) "
여기서, here,
이고, ego,
이며, Lt;
이다. to be.
5. 요격지점에 관한 해5. A description of the intercept point
5-1. 유도탄의 속도가 일정한 경우5-1. If the speed of the missile is constant
유도탄 속도가 일정한 경우 이므로 수학식 6은 수학식 7과 같이 나타낼 수 있다. Constant speed of missile Equation (6) can be expressed by Equation (7).
[수학식 7]&Quot; (7) "
여기서, here,
이다. to be.
따라서, 2차 방정식의 근의 공식으로부터 다음과 같이 해를 구할 수 있다.Therefore, from the root equation of the quadratic equation, the solution can be obtained as follows.
① 인 경우① If
이 경우 두 실근이 존재하는데 인 경우 이어야 하므로 그의 해인 무차원 요격지점()은 다음의 수학식 8과 같다.In this case, there are two exoskeletons If So it's his year, the dimensionless intercept point ( ) Is expressed by the following equation (8).
[수학식 8]&Quot; (8) "
② 인 경우② If
이 경우 해인 무차원 요격지점()은 다음의 수학식 9와 같다.In this case, a non-dimensional intercept point ( ) ≪ / RTI >
[수학식 9]&Quot; (9) "
③ 그 외의 경우③ 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) "
5-2. 유도탄의 속도가 변하는 경우5-2. When the speed of the missile is changed
유도탄의 속도가 변하는 경우 이므로 수학식 6을 다음의 수학식 11과 같이 다시 쓸 수 있다.When the speed of the missile is changed The equation (6) can be rewritten as the following equation (11).
[수학식 11]&Quot; (11) "
여기서,here,
이다. 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) "
여기서, here,
이다. to be.
따라서, 해는 다음과 같다. Therefore, the solution is as follows.
① D > 0인 경우① When D> 0
[수학식 13]&Quot; (13) "
여기서,이다.here, 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) "
③ 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) "
여기서, here,
이다. 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,
② 유도탄 비행속도를 고려한 비행거리(SL)는 0 이상이다. 즉,② Flight distance ( SL ) considering missile flying speed is more than zero. In other words,
여기서, 이다.here, 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,
여기서, 이고, 는 양의 작은 값이다. here, ego, 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) "
여기서, 는 유도탄이 등속비행()인 경우의 해로서, 수학식 8~10에 의해 결정된 값이다.here, (A) ), 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)
무차원 요격지점()을 이용하여 다음의 수학식 17과 같이 요격예상지점을 산출한다.Non-dimensional intercept point ( ) To calculate an intercept prediction point as shown in the following Equation (17).
[수학식 17]&Quot; (17) "
여기서, 은 유도탄의 수평축 절대좌표, 는 표적의 수평 및 수직 절대좌표이다. 따라서, 경로점(WP)의 위치는 다음의 수학식 18과 같다.here, The horizontal axis absolute coordinates of the missile, 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) "
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) "
여기서,이다.here, 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.
무차원 요격지점에 대한 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.
유도탄의 속도가 등속인지 판단하는 단계; 및
유도탄의 속도가 등속인 경우 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차 방정식의 요격지점 관계식의 판별식을 산출하는 단계; 및
상기 산출된 판별식의 크기에 따라 하나의 해 또는 다중해를 산출하는 단계;를 포함하는 것을 특징으로 하는 유도탄의 요격지점 실시간 예측 방법. 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.
상기 선택된 해중에서 소정 평가지수를 만족하는 해를 요격지점에 대한 해로 최종 결정하는 단계;를 더 포함하는 것을 특징으로 하는 유도탄의 요격지점 실시간 예측 방법. 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.
요격지점이 유도탄과 표적 사이에 위치하고,
유도탄 비행속도를 고려한 비행거리는 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.
다음 식에 의해 계산하는 것을 특징으로 하는 유도탄의 요격지점 실시간 예측 방법.
여기서, 는 구해진 요격지점에 대한 해석해, 은 유도탄의 수평축 절대 좌표, 는 표적의 절대 좌표이다.The method according to claim 1, , )silver
And calculating by using the following equation.
here, Interpreted the obtained intercept point, The horizontal axis absolute coordinates of the missile, Is the absolute coordinate of the target.
다음 식에 의해 계산하는 것을 특징으로 하는 유도탄의 요격지점 실시간 예측 방법.
여기서, 이고, , 은 현 시점의 유도탄의 속력, 는 글라이딩 단계의 유도탄의 평균 가속도, 는 표적의 수평축 위치이다. The method of claim 1,
And calculating by using the following equation.
here, ego, , The speed of the current missile, The average acceleration of the missile in the gliding phase, Is the horizontal axis position of the target.
여기서, , 는 요격지점의 절대 위치, d는 요격지점과 경로점의 상대거리, 은 유도탄의 경로점 통과각이다.2. The method according to claim 1, further comprising: calculating a route point of the guided car using the obtained solution, , ) Is calculated by the following equation.
here, , Is the absolute position of the intercept point, d is the relative distance between the intercept point and the path point, Is the path passing angle of the missile.
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)
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)
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 |
-
2014
- 2014-12-10 KR KR1020140177705A patent/KR101658464B1/en active IP Right Grant
Patent Citations (2)
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)
Title |
---|
논문자료(2010.05.27) * |
Cited By (8)
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 |