LMIs in Control/pages/LMI for Nonrotating Missiles Attitude Control yaw roll channel

LMI for Attitude Control of Nonrotating Missles, Yaw/Roll Channel

Deriving the exact dynamic modeling of a missile is a very complicated procedure. Thus, a simplified model is used to model the missile dynamics. To do so, we consider a simplified attitude system model for the yaw/roll channel of the system. We aim to achieve a non-rotating motion of missiles. Note that the attitude control design for the yaw/roll channel and the pitch channel can be solved exactly in the same way except for different representing matrices of the system.

The System

The state-space representation for the yaw/roll channel can be written as follows:

{\begin{aligned}{\dot {x}}(t)&=A(t)x(t)+B_{1}(t)u(t)+B_{2}(t)d(t)\\y(t)&=C(t)x(t)+D_{1}(t)u(t)+D_{2}(t)d(t)\end{aligned}} where $x=[\beta \quad w_{y}\quad w_{x}\quad \delta _{x}\quad \delta _{y}]^{\text{T}}$ , $u=[\delta _{xc}\quad \delta _{yc}]^{\text{T}}$ , $y=[n_{z}\quad w_{x}]^{\text{T}}$ , and $d=\delta _{z}$ are the state variable, control input, output, and disturbance vectors, respectively. The paprameters $\alpha$ , $w_{z}$ , $\delta _{z}$ , $\delta _{zc}$ , $n_{y}$ , $\beta$ , and $w_{y}$ stand for the attack angle, pitch angular velocity, the elevator deflection, the input actuator deflection, the overload on the side direction, the sideslip angle, and the yaw angular velocity, respectively.

The Data

In the aforementioned yaw/roll channel system, the matrices $A(t),B_{1}(t),B_{2}(t),C(t),D_{1}(t),$ and $D_{2}(t)$ are given as:

{\begin{aligned}A(t)={\begin{bmatrix}A_{11}(t)&A_{12}(t)\\0&A_{22}(t)\end{bmatrix}}\end{aligned}} where

{\begin{aligned}A_{11}(t)={\begin{bmatrix}-b_{4}(t)&1&{\frac {\alpha (t)}{57.3}}\\-{\acute {b}}_{1}(t)b_{4}(t)-b_{2}(t)&-{\acute {b}}_{1}(t)-{\acute {b}}_{1}(t)&{\frac {J_{y}-J_{z}}{57.3J_{x}}}{w}_{z}(t)-{\frac {{\acute {b}}_{1}\alpha (t)}{57.3}}\\c_{2}(t)&{\frac {J_{y}-J_{z}}{57.3J_{x}}}{w}_{z}(t)&-c_{1}(t)\end{bmatrix}}\end{aligned}} {\begin{aligned}A_{12}(t)={\begin{bmatrix}0&-b_{5}(t)\\0&-{\acute {b}}_{1}(t){\acute {b}}_{5}-b_{3}(t)\\-c_{3}(t)&c_{4}(t)\end{bmatrix}}\end{aligned}} {\begin{aligned}A_{22}(t)=-{\frac {1}{\tau _{x}\tau _{y}}}{\begin{bmatrix}\tau _{y}&0\\0&\tau _{x}\end{bmatrix}}\end{aligned}} and

{\begin{aligned}B_{1}(t)={\begin{bmatrix}0&0\\0&0\\0&0\\{\frac {1}{\tau _{x}}}&0\\0&{\frac {1}{\tau _{y}}}\end{bmatrix}},\quad B_{2}(t)={\begin{bmatrix}a_{6}(t)\\-{\acute {b}}_{1}(t)a_{6}(t)\\0\\0\\0\end{bmatrix}}\end{aligned}} {\begin{aligned}C(t)=-{\frac {1}{57.3g}}{\begin{bmatrix}V(t)b_{4}(t)&0&0&0&V(t)b_{5}(t)\\0&0&-57.3g&0&0\end{bmatrix}}\end{aligned}} {\begin{aligned}D_{1}(t)=0,\quad D_{2}(t)=-{\frac {V(t)}{57.3g}}{\begin{bmatrix}b_{6}(t)\\0\end{bmatrix}}\end{aligned}} where $a_{1}(t)\sim a_{6}(t),\quad b_{1}(t)\sim b_{7}(t),{\acute {a}}_{1}(t),{\acute {b}}_{1}(t)$ and $c_{1}(t)\sim c_{4}(t)$ are the system parameters. Moreover, $V$ is the speed of the missle and $J_{x}$ , $J_{y}$ , and $J_{z}$ are the rotary inertia of the missle corresponding to the body coordinates.

The Optimization Problem

The optimization problem is to find a state feedback control law $u=Kx$ such that:

1. The closed-loop system:

{\begin{aligned}{\dot {x}}&=(A+B_{1}K)x+B_{2}d\\z&=(C+D_{1}K)x+D_{2}d\end{aligned}} is stable.

2. The $H_{\infty }$ norm of the transfer function:

$G_{zd}(s)=(C+D_{1}K)(sI-(A+B_{1}K))^{-1}B_{2}+D_{2}$ is less than a positive scalar value, $\gamma$ . Thus:

$||G_{zd}(s)||_{\infty }<\gamma$ The LMI: LMI for non-rotating missle attitude control

Using Theorem 8.1 in , the problem can be equivalently expressed in the following form:

{\begin{aligned}&{\text{min}}\quad \gamma \\&{\text{s.t.}}\quad X>0\\&{\begin{bmatrix}(AX+B_{1}W)^{T}+AX+B_{1}W&B_{2}&(CX+D_{1}W)^{T}\\B_{2}^{T}&-\gamma I&D_{2}^{T}\\CX+D_{1}W&D_{2}&-\gamma I\end{bmatrix}}<0\end{aligned}} Conclusion:

As mentioned, the aim is to attenuate the disturbance on the performance of the missile. The parameter $\gamma$ is the disturbance attenuation level. When the matrices $W$ and $X$ are determined in the optimization problem, the controller gain matrix can be computed by:

$K=WX^{-1}$ Implementation

A link to Matlab codes for this problem in the Github repository: