The paper considers a method to build optimal conduction law solve the problem of local optimization of the omniscient form based on Letov-Kalman approach. The obtained optimum lead law is similar to the traditional proportional lead law, but the ratio coefficient varies with missile-target distance. The law of conductivity can be realized for the active radio self-guided missile class.
In the field of rocket guidance, number of studies suggests different laws of conduction. The law of proportional approach with the fixed coefficients (traditional proportional conduction law) has been widely applied since its high realization ability. However, the paper [1,2] has shown that this lead law has many limitations in the case of self-guided missiles with maneuverable targets. Therefore, the law of conduction has been continuously studied to find a new law of conduction that can improve the quality of missile guidance in the case of mobile targets [3].
Consider a linear control system with a defined structure as follows:
x˙y(t) = Fy(t)xy(t) + By(t)u(t) + ξy(t) (1)
Where: xy
- the output state vector of the system
ξy
- the systemic noise vector of the Gaussian white noise with the mathematical expectation of 0;
u(t)- control signal vector
Fy, By -Control efficiency matrix.
The index function of the local quality:
I = [xT(t) - xy(t)]TQ[xT(t) - xy(t)] - ∫0tUT(t)Ku(t)dt (2)
Which: XT - the request state vector
Q- matrix of penalty coefficients according to accuracy at time t
K- matrix of the penalty coefficients according to the magnitude of the control signal
Solving the problem of finding optimal control signals for system (1) and slab (2) by Bellman dynamic planning method.
u = K−1BTyQ[xT − xy] (3)
The spatial self-conducting geometry is described by the vector equation
λ→¨ = [− 2R˙R]λ→˙ + [− R¨R]λ→ − 1Ra→ + 1Ramt−→− (4)
Which: amt−→−, a→
, is the total acceleration vector of the target and missile λ→ is the unit vector of missile-target distance vector with R→ = R⋅λ→
To describe the coordinates of the vectors λ→, λ→˙, λ→¨, a→, amt−→−
in the fixed ground coordinate system (inertial coordinate system) we use vectors. The corresponding state silk is as follows:
⎧⎩⎨⎪⎪λq = [λxqλyqλzq]T; λ˙q = [λ˙xqλ˙yqλ˙zq]T;λ¨q = [λ¨xqλ¨yqλ¨zq]Taq = [axqayqazq]T; amtq = [amtxqamtyqamtzq]T;⎫⎭⎬⎪⎪ (5)
The state vector describes the instantaneous state of the line of sight:
Xy = [λTqλ′Tq]T (6)
Put:
Fy = ⎡⎣O3− R¨RI3I3− 2R˙RI3⎤⎦; By = − 1R[O3I3]; u = aq; Bmt = 1R[O3I3]; ξmt = amtq (7)
In which: O3- is a square matrix of order 3, including only zeros; I3 is a 3 × 3 unit matrix.
From equations (5) to (7) we have self-conducting geometry described in terms of equations of state
x˙y = Fyxy + Byu + Bmtξmt (8)
Select the optimal slab as follows:
I = qh⋅ω2λ(t) + ∫01Ku⋅a2(t)dt (9)
In which: qh > 0 is the penalty coefficient according to the instantaneous slip; ku > 0 is the penalty factor for the total acceleration (or overload) required up to time t.
Vector depicts the desired state of sight
xT = [OTqT]T (10)
Put: Qh = qh[O3O3O3I3]; Ku = ku.I3 (11)
Then the (9) equivalent of the slab
I = [xT(t) - xy(t)]TQh[xT(t) - xy(t)] - ∫01uT(t)Kuu(t)dt (12)
When synthesizing the conduction law in the absence of information about the target acceleration, the target acceleration can be considered as an unknown random effect. The state equation without the impact of the target acceleration has the form
x˙y = Fyxy + Byu (13)
Applying optimal control algorithm for system (13) and local criteria blade (12) we have
u = aq = (1R⋅qhku)λ˙q (14)
From (5) and (14) we have optimal lead law corresponding to the criterion (9) or (12)
a→ = (1R⋅qhku)λ˙ (15)
The symbol ωλ−→
is the vector of the angle of rotation of the line of sight, we have ωλ−→ ⊥ λ→
and
λ→ = (ω→λxλ→) = (ωλ.ω0λ−→xλ→) = ωλ.λ0−→ (16)
With λ0−→, ω0λ−→
are unit vectors of λ→ and ω→λ
respectively
From (16) we have (λ→,λ→˙,ω→λ)
perpendicular to each other and set up the perpendicular coordinate system in space, so we can describe the direction of the required acceleration vector of the name. The rocket corresponds to the law of missile guidance (15) as shown in the Figure 1.
According to [4] the line of sight kinetics is described by the equation
R.ω˙λ +˙ 2R˙.ωλ = (a⊥mt − a⊥) (17)
In which: a⊥mt, a⊥
are the projections of the target acceleration vector, the missile is perpendicular to the line of sight and lies in the self-guided plane.
From the picture and according to (16) we have
a⊥ = a→.λ0−→ = (1R⋅qhku)ωλ (18)
Replace (18) into (17) and transform we have
R(1R⋅qhku) − 2Vc⋅ωλ + ωλ = a⊥mt(1R⋅qhku) − 2Vc (19)
In which: V = − R˙
is the speed at which the missile approaches the target
According to [5], the instantaneous slip h of the self-conductive process is determined by
h = R2ΔV⋅ωλ (20)
Where ΔV is the magnitude of the vector ΔV−→− = (V→ − Vmt−→−)
The instantaneous slip (h) of the self-conduction process is proportional to the rotation speed of the line of sight so allowing the self-guided control ring to ensure that the instantaneous slip h → 0 will correspond to the control of the rocket maneuver so that ωλ → 0. Equation (19) shows that when the missile is self-conducting according to the law (15), the self-conducting control ring corresponds to the first order of inertia stage with the time constant:
τωk = R1R⋅qhku − 2Vc (21)
When a missile attacks a mobile target, the total active self-driving time is usually very short. Therefore, to ensure the required slip at the meeting point, we set a requirement to limit the transient time of the self-conducting control ring right from the start of self-conduction:
τωk ≤ Tcp (22)
Tcp- allowable limit of self-conductive control loop time constant.
From (21) and (22) we have the boundary condition of (22) being equivalent
qhku = 2RVc + R2Tcp (23)
For (qh, ku) choose according to (23) then the law of leading (15) becomes
a→ = (2Vc + R2Tcp)λ→˙ (24)
The law of conductivity (24) ensures that the transient time constant of the self-conductive control loop is equal to the limit of the permissible value. Unlike the traditional proportional conduction law, this law of conduction has a proportional coefficient that decreases with distance and depends on the time constant limit [6].
Simulations are performed with traditional proportional conduction law (coefficient K = 3) and optimal conduction law (24) under the same conditions: Rp0 = 6 [km]; target overload nmt = 9, flight time tq0 = 0.3 [s] Figure 2, Figure 3 and Figure 4.
• When the target is strongly mobile, the slip at the meeting point of the orbit corresponding to the law of conduction (24) is significantly smaller than the slip at the meeting point of the conduction trajectory corresponding to the traditional proportional conduction law.
• With traditional proportional conduction law, at the beginning of self-conductivity, the required missile overload is of small value so that the initial slip is not quickly eliminated from the start thanks to the maximum effective overload of the name fire.
• With the law of conductivity (24), because slip is always required to reduce quickly with limited transition time, the day from the start of self-conductivity, according to the requirements of the law of conductivity, the missile using the overload is maximized to quickly reduce the initial slip. Therefore the initial slip is reduced faster and the slip at meeting point is guaranteed to be small enough as required.
The survey simulation results show that the obtained optimal conduction law meets the quality requirements of the self-conductive control loop and is more effective than the traditional ratio approach law when the target is maneuverable. This law of conductivity can be applied to the active radio self-guided missile class, which helps to improve the conductivity of these missiles in aerial combat with highly maneuverable flying targets.