Stability of Single-Channel Rolling Aerospace Vehicle with Semi-Automatic Command to Line of Sight

1. Introduction

In modern anti-armor warfare, guided missiles are a primary tool, often employing "three-point guidance" methods that rely on the alignment of the carrier, missile, and target. These systems can be categorized by their guidance mechanism, with wire-guided missiles typically using angular error control, while laser or radio-beam-guided missiles utilize linear error control. This paper focuses on the former, specifically on wire-guided missiles that require a direct hit on heavily armored targets. This is achieved through a Semi-Automatic Command to Line of Sight (SACLOS) guidance method, where an operator tracks the target, and a guidance signal is automatically generated and transmitted to the missile via a wire.

For subsonic flight regimes, where aerodynamic control surfaces have limited effectiveness, these missiles often rely on gasodynamic controls, such as thrust vector control, to perform necessary maneuvers. The use of slow-rolling missiles in engaging heavily armored targets has become increasingly prevalent in modern weaponry. This design offers several advantages, including a reduction in the missile's size and mass. The rolling motion also serves to cancel out disturbing forces and torques that arise from aerodynamic or gasodynamic asymmetries. Furthermore, the roll allows for control to be executed on a single channel with a single pair of mobile nozzles, simplifying the control system and reducing energy requirements, which is conducive to the development of portable missile systems.

Although these systems are known and technically developed [1,2,3], they represent a special aircraft class, for which, due to the specifics of the movement, theoretical models necessary to clarify some aspects regarding the improvement of their functionality are still in development.

The guidance and control of rolling missiles have been extensively studied. Research has addressed two-channel rolling missiles [4,5], including dual-spin projectiles with high-velocity rotating bodies for gyroscopic stability and lower-velocity rotating tips with canards for trajectory control [6]. Other studies have focused on single-channel missiles with time-modulated commands [7,8,9,10,11]. Recent work on two-channel missiles has explored autopilot design with double [12] or triple [13] feedback loops, the stability of precession motion (conical) [4,5], and the influence of signal errors from the target tracking system [12,14]. Another area of research has been the development of sensors to determine the missile's attitude, particularly the roll angle, which is crucial for command signal generation [15,16]. For instance, some have proposed a new combination of silicon micro-gyroscopes to demodulate yaw, pitch, and spin angular velocities. Another strategy suggested in the literature is the application of a two-step adaptive augmented unscented Kalman filter (TS-AAUKF) [16], which serves to calibrate the magnetometer for accurate roll angle estimation.

To address what was done in the specific case of the single-channel missile, the following section analyses the results from previous research [7–11, 17] in detail.

Reference [7] proposes a digital controller as an alternative to the conventional analogue time-modulation controller. The authors demonstrate that by substituting the analogue phase control with a digital amplitude control system, it is possible to achieve decoupling of the yaw and pitch channels. This method effectively mitigates errors arising from actuator delay and sensor noise. However, the primary limitation of this solution (at the time of proposal) was the increased complexity of the missile's control system, which consequently led to a higher overall cost.

The stability of a rolling rocket with single-channel control is examined in [8] through the development of closed-form analytical solutions from linearized dynamics. This analysis uses a mixed reference frame: the control command is processed in the body frame, while the system dynamics is modelled in the non-rolling Resal frame. Furthermore, the command signal is generated via a series expansion, facilitating a comparative analysis between an idealised rectangular command and a practical implementation that accounts for actuator delays.

In reference [9], a linear model for the commanded motion of a slow-rotating rocket is developed from the equations of motion in the Resal frame. This model uses complex coefficients to account for the coupling between channels caused by gyroscopic and Magnus terms, which allows for the characterization of both nutation and precession motion and the definition of flight quality parameters. This linear analysis concludes in a structural scheme of the commanded object that is specific to slow-rolling missiles with complex coefficients. Unlike the serial development approach in [8], the proposed solution in [9] generates the switch command by considering the sign of a trigonometric expression that depends on the roll angle, as detailed in [18]. The paper analyses two types of single-channel, slow-rolling guided missiles: a portable anti-aircraft homing missile and a portable anti-tank wire-guided missile. In contrast to [8], which formulates the command in the body frame, [9] derives an equivalent command in the Resal frame from the time-modulated command. This equivalent command is then used to construct structural diagrams for both the homing and wire-guided missile types. Due to the complexity introduced by the structural schemes' complex coefficients, the scope was limited to determining the characteristic polynomial's roots and the phase-frequency/amplitude-frequency diagrams for the wire-guided rocket, emphasizing the differences from a non-rolling rocket.

The feasibility of controlling an autonomous rolling rocket via a single-channel reactive command is explored in [10]. The work contrasts a simplified two-nozzle Reaction Control System (RCS) with a conventional six-nozzle configuration. While the traditional system ensures complete three-axis attitude control, the proposed two-nozzle alternative provides actuation for only the pitch and yaw channels. This latter approach leaves the roll angle uncontrolled, thus allowing for free roll rotation of the vehicle.

In [11], starting from the results in paper [9], the case of the anti-aircraft rolling homing missile was analysed, performing a stability analysis with the Frank-Wall criterion.

Reference [17] addresses the optimal control problem for the minimum-time attitude manoeuvring of a rolling missile. In this formulation, a single reaction jet is used to generate the necessary transverse control moments. The problem considers a rocket with an arbitrary initial transverse angular velocity. The objective is to steer the missile to a specified final attitude in the shortest possible time, while simultaneously nullifying the final transverse angular velocity.

The current paper builds upon previous work on single-channel missiles, particularly the findings in [9], by analyzing the stability of a single-channel, slow-rolling SACLOS missile using the Frank-Wall stability criterion [19]. A key distinction between two-channel and single-channel missiles lies in the command formation. While two-channel systems use proportional deflection of both canards, single-channel missiles employ a time-modulated maximum deflection of one canard, managed by a switching function. We will also consider the influence of the guidance gain ($k$) and time guidance (${\tau }_c$) of the guided device system, leading to a fifth-order characteristic polynomial, which introduces additional stability conditions that will be analyzed in detail. The effect of phase shift, introduced by actuator delay and inaccuracies in roll angle determination, will also be examined.

Since this paper represents an alternative development of the paper [11], for the case of single-channel slow-rolling missiles with thrust vector control and SACLOS guidance, we will present only the specific elements that appear in this new case, for the common aspects referring to the paper [11].

Contributions of This Paper

This paper presents several key contributions to the study of single-channel, slow-rolling missiles. The primary novelty lies in the comprehensive stability analysis of a SACLOS missile using the comparative Routh-Hurwitz (R-H) and Frank-Wall (F-W) stability criteria and nonlinear analysis. This analysis examines the influence of key guidance system parameters, leading to the definition of specific restrictions and stability domains. To support this analysis, the paper establishes a method for determining the F-W stability parameters for the resulting 5th-order characteristic polynomial, with detailed calculation procedures provided in Appendix E. A detailed structural scheme for a single-channel, slow-rolling SACLOS missile is also developed and presented. Furthermore, this work includes a thorough analysis of the stability and command derivative magnitudes derived from the linearized model, as well as the flight quality parameters. Finally, the stability domains identified through the linear analysis using the R-H and F-W criteria are validated against a nonlinear 6-DOF model, confirming the accuracy and applicability of the results.

The novelty of the paper consists of:

• Comparative analysis between the exact solution and the approximate solution for the equilibrium flight parameters (Section 2);
• The detailed structural scheme for the SACLOS slow-rolling missile with a single-channel (Figure 12);
• Synthesis of the guidance system using the allocation of poly-zeros according to the methodology “Standard coefficient method” (Section 5.6);
• The analysis of the stability of SACLOS missiles using the Routh-Hurwitz criterion for the simplified case without roll. This analysis follows the influence of the main parameters in the missile's guidance system, and stability domains will be defined (Section 7.1).
• The stability analysis for slowly rolling, single-channel SACLOS missiles, employing the Frank-Wall stability criterion. The investigation examines the influence of key guidance system parameters on stability, utilizing a comprehensive model that incorporates higher-order dynamic terms. Based on this analysis, operational constraints are identified, and stability domains are defined for several of these critical parameters. (Section 7.2).
• The evaluation of the Frank-Wall (F-W) stability parameters for 5th order polynomials (the methodology is detailed in Annex E of the paper).
• The validation of linear analysis with R-H and F-W stability criteria using a nonlinear model (Section 7.3)

2. Basic Movement

Since the general equations of motion in the body frame and the Resal frame were presented in the paper [11], to linearize the equations of motion, we will now analyze the basic motion in the Resal frame, because in this case, unlike in [11], we have gasodynamic control.

We will consider a symmetrical evolution in the vertical plane as basic motion. The lateral movement elements are set to zero: $\begin{array}{ccc}{r}^{\mathrm{*}}=0;& \beta =0;& {\delta }_n^{\mathrm{*}}=0\end{array}$. From Figure 1 we can write the equilibrium equations in the vertical plane:

$$\begin{gathered} {ma}_z^{\mathrm{*}}=F_0{C}_z^A+T_0{C}_z^T+G\mathit{cos}\theta ; m{a}_y^{\mathrm{*}}=F_0{C}_y^A+T_0{C}_y^T; \\ {H}_0{C}_m^A+U_0{C}_m^T+r^{\mathrm{*}}pA=0;H_0{C}_n^A+U_0{C}_n^T-q^{\mathrm{*}}pA=0. \end{gathered}$$

(1)

If we detail the lateral accelerations: $a_z^{*}=-q^{*}{u}^{*}$; $a_y^{*}=r^{*}{u}^{*}$, the relations become:

$$\begin{gathered} C_z^A+{{T_{0}F}_0^{-1}C}_z^T+{mq}^{\mathrm{*}}{u}^{\mathrm{*}}{F}_0^{-1}+G{F}_0^{-1}\mathit{cos}\theta =0; C_y^A+{T_{0}F}_0^{-1}{C}_y^T-m{r}^{\mathrm{*}}{u}^{\mathrm{*}}{F}_0^{-1}=0 \\ {C}_m^A+{T_{0}F}_0^{-1}{C}_m^T+r^{\mathrm{*}}pA{H}_0^{-1}=0;C_n^A+{T_{0}F}_0^{-1}{C}_n^T-q^{\mathrm{*}}pA{H}_0^{-1}=0 \end{gathered}$$

(2)

If the terms of gasodynamic damping [20] are neglected and it is noted

${\widehat{x}}_T=x_T/l$; $C_T=T{T}_0^{-1},$ the thrust components become:

$$\begin{gathered} C_z^T=C_{T}sin{\delta }_m^{*}; C_y^T=-C_{T}cos{\delta }_m^{*}sin{\delta }_n^{*}\cong -C_{T}sin{\delta }_n^{*} \\ {C}_m^T={-\widehat{x}}_T{C}_z^T={-\widehat{x}}_T{C}_T{sin\delta }_m^{*};C_n^T={\widehat{x}}_T{C}_y^T={-\widehat{x}}_T{C}_T{sin\delta }_n^{*}; \end{gathered}$$

(3)

where, as shown in [9,21], the equivalent gasodynamic deflection to the pitch axis of the Resal frame can be taken for twin-channel missiles approximately equal to the deflection in a control plane, and for single-channel missiles it is obtained by integrating the projection of the deflection from the control plane over a period of rotation along the direction of the pitch axis of the quasi-linked frame.

Note: Since in the body frame the position of application of the gasodynamic command is behind the center of mass, it follows $x_T<0$.

The aerodynamic coefficients are developed in this way:

$$\begin{gathered} C_z^A=b_1\alpha +b_2{\alpha }^3+b_4{\widehat{q}}^{\mathrm{*}}+b_{13}\widehat{p}\beta ; C_y^A={-b}_1\beta -b_2{\beta }^3-b_4{\widehat{r}}^{\mathrm{*}}+b_{13}\widehat{p}\alpha \\ {C}_m^A=d_1\alpha +d_4{\widehat{q}}^{\mathrm{*}}+d_{13}\widehat{p}\beta \end{gathered}$$

(4)

$$; C_n^A=d_1\mathsf{\beta }+d_4{\widehat{r}}^{\mathrm{*}}-d_{13}\widehat{p}\alpha$$

(4)

For analyzing the magnitude terms in these expressions, and those that will be elaborated later, we will use, according to [22,23], the following parameters:

Small values:

$t^{*}=l/V$- reference time;

$\widehat{g}=gl/V^2$- reduced gravitational acceleration

Large values:

$\tilde{m}=2m/\left(\rho Sl\right)$- the reduced mass of the aircraft

$i_B=2B/\left(\rho S{l}^3\right)$-the reduced moment of inertia in pitch

Introducing these parameters into relations (2), we obtain:

$$\begin{gathered} b_1\alpha +b_2{\alpha }^3+b_4{\widehat{q}}^{\mathrm{*}}+b_{13}\widehat{p}\beta +{{\tilde{T}}_{0}C}_{T}sin{\delta }_m^{\mathrm{*}}+\tilde{m}{\widehat{q}}^{\mathrm{*}}\mathit{cos}{\gamma }^{\mathrm{*}}+ \\ \tilde{m}\widehat{g}\mathit{cos}(\gamma +\alpha )=0 \\ {b}_1\beta +b_2{\beta }^3+b_4{\widehat{r}}^{\mathrm{*}}-b_{13}\widehat{p}\alpha +{\tilde{T}}_0{C}_{T}sin{\delta }_n^{\mathrm{*}}+\tilde{m}{\widehat{r}}^{\mathrm{*}}\mathit{cos}{\gamma }^{\mathrm{*}}=0; \\ {d}_1\alpha {+d}_2{\alpha }^3+d_4{\widehat{q}}^{\mathrm{*}}+d_{13}\widehat{p}\beta -{\tilde{T}}_0{\widehat{x}}_T{C}_T{sin\delta }_m^{\mathrm{*}}+{2i}_B{\widehat{\omega }}_p{\widehat{r}}^{\mathrm{*}}=0; \\ {d}_1\beta +d_2{\beta }^3+d_4{\widehat{r}}^{\mathrm{*}}-d_{13}\widehat{p}\alpha -{\tilde{T}}_0{\widehat{x}}_T{C}_T{sin\delta }_n^{\mathrm{*}}-{2i}_B{\widehat{\omega }}_p{\widehat{q}}^{\mathrm{*}}=0 \end{gathered}$$

(5)

where we noted ${\tilde{T}}_0={T_{0}F}_0^{-1}$; ${\mathit{tan}}^2{\gamma }^{*}={\mathit{tan}}^2\alpha +{\mathit{tan}}^2\beta$ while ${\widehat{\omega }}_p=A\widehat{p}/\left(2B\right)$ is the dimensionless angular rate of precession, and $\gamma =\theta -\alpha$ is the flight-path inclination angle.

If we impose angular velocities ${\widehat{q}}^{*}$, ${\widehat{r}}^{\mathrm{*}}$ and the angle of inclination of the trajectory $\gamma$, the system (5) can be solved numerically with a Newton algorithm [24] in relation to the sizes $\alpha$, $\beta$, ${\delta }_m^{*}$, ${\delta }_n^{*}$ which represent the flight parameters for the basic motion under consideration.

If we impose a translational motion with zero angular velocity of pitch and yaw (${\widehat{q}}^{\mathrm{*}}=0; {\widehat{r}}^{\mathrm{*}}=0)$, and neglecting the aerodynamic terms of higher order ($b_2=0; ;b_{13}=0;d_2=0; d_6=0; d_{13}=0)$, the system decouples and becomes:

$$\begin{gathered} b_1\alpha +{{\tilde{T}}_{0}C}_{T}sin{\delta }_m^{*}+\tilde{m}\widehat{g}\mathit{cos}(\gamma +\alpha )=0; d_1\alpha -{\tilde{T}}_0{\widehat{x}}_T{C}_T{sin\delta }_m^{*}=0; \\ {b}_1\beta +{\tilde{T}}_0{C}_{T}sin{\delta }_n^{*}=0;d_1\beta -{\tilde{T}}_0{\widehat{x}}_T{C}_T{sin\delta }_n^{*}=0 \end{gathered}$$

(6)

One can notice that for lateral motion the last two homogeneous equations lead to the trivial solution: $\beta =0$; ${\delta }_n^{\mathrm{*}}=0,$while the first two equations, corresponding to longitudinal motion, for the given $\gamma$ lead to a transcendent equation in the angle $\alpha$ and a calculus relation for the equivalent thrust deflection in pitch ${\delta }_m^{\mathrm{*}}$.

Thus, eliminating the incidence between the equations of longitudinal motion, we obtain:

$${sin\delta }_{me}^{\mathrm{*}}=-{\displaystyle \frac{\zeta }{{\widehat{x}}_T+\zeta }}{\displaystyle \frac{G}{T}}\mathit{cos}(\gamma +\alpha )$$

(7)

where the static stability reserve was noted.$\zeta =d_1/b_1$.

If the terms of the thrust are eliminated, the following is obtained:

$${\alpha }_e=-{\displaystyle \frac{G\mathit{cos}(\gamma +{\alpha }_e)}{F_0{b}_1}}{\displaystyle \frac{{\widehat{x}}_T}{\left({\widehat{x}}_T+\zeta \right)}}$$

(8)

Imposing the value of the climb angle $\gamma$, implicitly for the particular case of horizontal flight $\gamma =0$, the system is reduced to the relation (8) which becomes a transcendent equation that is solved iteratively in relation to the incidence to which is added the explicit relation (7) which allows the determination of the equivalent thrust angular deflection.

Figure 2. Comparison between the exact solution (5) and the approximate solution (7) (8) for the equilibrium command - initial mass values.

Figure 2. Comparison between the exact solution (5) and the approximate solution (7) (8) for the equilibrium command - initial mass values.

Figure 3. Comparison between the exact solution (5) and the approximate solution (7) (8) for the equilibrium command - final mass values.

Figure 3. Comparison between the exact solution (5) and the approximate solution (7) (8) for the equilibrium command - final mass values.

These parameters (${\alpha }_{ e}{, \delta }_{me}^{\mathrm{*}},{\beta }_{ e}{, \delta }_{ne}^{\mathrm{*}})$represent the basic movement of the rolling missile, presented in Appendix B.

From Figure 4 it can be seen that the equilibrium incidence in the lateral plane (β) is negative, which will lead to a deviation of the trajectory to the left of the firing plane, a phenomenon specific to aerodynamically stabilized rolling missiles.

3. The Linear Equations

To obtain the linear equations, we will consider the hypotheses from the work [11].

Also, for the symmetry, for the gasodynamic commands, we will consider the small angular deflections so that:

$${sin\delta }_m^{\mathrm{*}}\cong {\delta }_m^{\mathrm{*}};{cos\delta }_m^{\mathrm{*}}\cong 1; {sin\delta }_n^{\mathrm{*}}\cong {\delta }_n^{\mathrm{*}}$$

(9)

$$;{cos\delta }_n^{\mathrm{*}}\cong 1$$

(9)

Linearizing translational equations in the Resal frame [11], we obtain:

$$\begin{gathered} ∆\dot{\alpha }={\displaystyle \frac{F_o}{m{u}^{\mathrm{*}}}}{b}_1∆\alpha +\left(1+{\displaystyle \frac{F_{o}l}{mV{u}^{\mathrm{*}}}}{b}_4\right)∆q^{\mathrm{*}}+{\displaystyle \frac{F_{o}l}{mV{u}^{\mathrm{*}}}}{b}_9∆\dot{\alpha }+{\displaystyle \frac{F_{o}l}{mV{u}^{\mathrm{*}}}}{b}_{13}p∆\beta +{\displaystyle \frac{F_o{l}^2}{m{u}^{\mathrm{*}}{V}^2}}{b}_{14}p∆r^{\mathrm{*}} \\ +{\displaystyle \frac{T_o}{m{u}^{\mathrm{*}}}}{C}_T∆{\delta }_m^{\mathrm{*}}-{\displaystyle \frac{F_o}{m{u}^{\mathrm{*}}}}{b}_1∆{\alpha }_W+{\displaystyle \frac{1}{m{u}^{\mathrm{*}}}}∆Z_p^{\mathrm{*}} \\ ∆\dot{\beta }={\displaystyle \frac{F_o}{m{u}^{\mathrm{*}}}}{b}_1∆\beta +\left(1+{\displaystyle \frac{F_{o}l}{mV{u}^{\mathrm{*}}}}{b}_4\right)∆r^{\mathrm{*}}+{\displaystyle \frac{F_{o}l}{mV{u}^{\mathrm{*}}}}{b}_9∆\dot{\beta }-{\displaystyle \frac{F_{o}l}{mV{u}^{\mathrm{*}}}}{b}_{13}p∆\propto \\ -{\displaystyle \frac{F_o{l}^2}{m{u^{\mathrm{*}}V}^2}}{b}_{14}p∆q^{\mathrm{*}}-{\displaystyle \frac{T_o}{m{u}^{\mathrm{*}}}}{C}_T∆{\delta }_n^{\mathrm{*}}-{\displaystyle \frac{F_o}{m{u}^{\mathrm{*}}}}{b}_1∆{\beta }_W-{\displaystyle \frac{1}{m{u}^{\mathrm{*}}}}∆Y_p^{\mathrm{*}} \end{gathered}$$

(10)

On the other hand, linearizing rotational equations in the Resal frame [11], we obtain:

$$\begin{gathered} ∆{\dot{q}}^{*}={\displaystyle \frac{H_0}{B}}{d}_1∆\alpha +{\displaystyle \frac{H_{0}l}{BV}}{d}_4∆q^{*}+{\displaystyle \frac{H_{0}l}{BV}}{d}_9∆\dot{\alpha }+{\displaystyle \frac{H_{0}l}{BV}}{d}_{13}p∆\beta +\left({\displaystyle \frac{H_0{l}^2}{B{V}^2}}{d}_{14}-{\displaystyle \frac{A}{B}}\right)p∆r^{*}-{\displaystyle \frac{U_0}{B}}{{\widehat{x}}_{T}C}_T∆{\delta }_m^{*}-{\displaystyle \frac{H_0}{B}}{d}_1∆{\alpha }_W; \\ ∆{\dot{r}}^{*}={\displaystyle \frac{H_0}{C}}{d}_1∆\beta +{\displaystyle \frac{H_{0}l}{CV}}{d}_4∆r^{*}+{\displaystyle \frac{H_{0}l}{CV}}{d}_9∆\dot{\beta }-{\displaystyle \frac{H_{0}l}{CV}}{d}_{13}p∆\alpha -\left({\displaystyle \frac{H_0{l}^2}{C{V}^2}}{d}_{14}-{\displaystyle \frac{A}{C}}\right)p∆q^{*}-{\displaystyle \frac{U_0}{C}}{\widehat{x}}_T{C}_T∆{\delta }_n^{*}-{\displaystyle \frac{H_0}{C}}{d}_1∆{\beta }_W. \end{gathered}$$

(11)

For coupling the equations, similar to the paper [11], the following dual notations are made:

$$\begin{gathered} a_{\gamma }^{\alpha }={\displaystyle \frac{F_o}{m{u}^{\mathrm{*}}}}{b}_1; a_{\gamma }^q=1+{\displaystyle \frac{F_{o}l}{mV{u}^{\mathrm{*}}}}{b}_4; a_{\gamma }^{\dot{\alpha }}={\displaystyle \frac{F_{o}l}{mV{u}^{\mathrm{*}}}}{b}_9; a_{\gamma }^{\beta }={\displaystyle \frac{F_{o}lp}{mV{u}^{\mathrm{*}}}}{b}_{13}; a_{\gamma }^r={\displaystyle \frac{F_o{l}^{2}p}{m{V}^2{u}^{\mathrm{*}}}}{b}_{14} \\ ; b_{\gamma }^{\delta }={\displaystyle \frac{T_o}{m{u}^{\mathrm{*}}}}{C}_T={\displaystyle \frac{T_o}{m{u}^{\mathrm{*}}}}{C}_z^T=-{\displaystyle \frac{T_o}{m{u}^{\mathrm{*}}}}{C}_y^T; a_{\varpi }^{\alpha }={\displaystyle \frac{H_o}{B}}{d}_1; a_{\varpi }^q={\displaystyle \frac{H_{o}l}{BV}}{d}_4;a_{\varpi }^{\dot{\alpha }}={\displaystyle \frac{H_{o}l}{BV}}{d}_9; a_{\varpi }^{\beta }={\displaystyle \frac{H_{o}lp}{BV}}{d}_{13} \\ ;a_{\varpi }^r=\left({\displaystyle \frac{H_o{l}^2}{B{V}^2}}{d}_{14}-{\displaystyle \frac{A}{B}}\right)p;b_{\varpi }^{\delta }=-{\displaystyle \frac{U_o{\widehat{x}}_T}{B}}{C}_T={\displaystyle \frac{U_o}{B}}{C}_m^T={\displaystyle \frac{U_o}{B}}{C}_n^T, \end{gathered}$$

(12)

These parameters represent the stability and command derivatives that are presented in Appendix B.

Next, we will rewrite the coefficients of the equations of motion, in which we will include the previously defined parameters for the basic movement, obtaining:

$$\begin{gathered} a_{\gamma }^{\alpha }={\displaystyle \frac{b_1}{t^{*}\tilde{m}}}; a_{\gamma }^q=1+{\displaystyle \frac{b_4}{\tilde{m}}}\cong 1; a_{\gamma }^{\dot{\alpha }}={\displaystyle \frac{b_9}{\tilde{m}}}\cong 0; a_{\gamma }^{\beta }=-{\displaystyle \frac{\widehat{p}}{\tilde{m}{t}^{*}}}{b}_{13}\cong 0; a_{\gamma }^r={\displaystyle \frac{\widehat{p}}{\tilde{m}}}{b}_{14}\cong 0; b_{\gamma }^{\delta }={\displaystyle \frac{{\tilde{T}}_0}{\tilde{m}{t}^{*}}}{C}_T; \\ {a}_{\varpi }^{\alpha }={\displaystyle \frac{d_1}{t^{*2}{i}_B}}; a_{\varpi }^q={\displaystyle \frac{d_4}{t^{*}{i}_B}}; a_{\varpi }^{\dot{\alpha }}={\displaystyle \frac{d_9}{t^{*}{i}_B}}; a_{\varpi }^{\beta }={\displaystyle \frac{\widehat{p}}{t^{*2}{i}_B}}{d}_{13}; a_{\varpi }^r=\left({\displaystyle \frac{d_{14}}{t^{*}{i}_B}}-{\displaystyle \frac{A}{B{t}^{*}}}\right)\widehat{p}\cong -{\displaystyle \frac{A}{B{t}^{*}}}\widehat{p}=-{\displaystyle \frac{2{\widehat{\omega }}_p}{t^{*}}}; b_{\varpi }^{\delta }=-{\displaystyle \frac{{\tilde{T}}_0{\widehat{x}}_T}{t^{*2}{i}_B}}{C}_T. \end{gathered}$$

(13)

Also, we define the complex quantities:

$$\begin{gathered} \begin{array}{ccccc}{\gamma }^{\mathrm{*}}=\alpha -i\beta ;& {\varpi }^{\mathrm{*}}=q^{\mathrm{*}}-i{r}^{\mathrm{*}};& {\delta }^{\mathrm{*}}={\delta }_m^{\mathrm{*}}-i{\delta }_n^{\mathrm{*}};& {\gamma }_W^{\mathrm{*}}={\alpha }_W-{i\beta }_W;& F_p^{\mathrm{*}}=Z_p^{\mathrm{*}}+i{Y}_p^{\mathrm{*}};\end{array} \\ {a}_{\gamma }^{\gamma }=a_{\gamma }^{\alpha }+i{a}_{\gamma }^{\beta }; a_{\gamma }^{\varpi }=a_{\gamma }^q+i{a}_{\gamma }^r; a_{\varpi }^{\gamma }=a_{\varpi }^{\alpha }+i{a}_{\varpi }^{\beta }; a_{\varpi }^{\varpi }=a_{\varpi }^q+i{a}_{\varpi }^r; \end{gathered}$$

(14)

with: $i=\sqrt{-1}$.

In this case the equations (10) and (11) are:

$$\begin{gathered} \Delta {\dot{\gamma }}^{*}=a_{\gamma }^{\gamma }\Delta {\gamma }^{*}+a_{\gamma }^{\varpi }\Delta {\varpi }^{*}+a_{\gamma }^{\dot{\alpha }}\Delta {\dot{\gamma }}^{*}+b_{\gamma }^{\delta }\Delta {\delta }^{*}-a_{\gamma }^{\alpha }\Delta {\gamma }_W^{*}+{\displaystyle \frac{\Delta F_p^{*}}{m{u}^{*}}} \\ \Delta {\dot{\varpi }}^{*}=a_{\varpi }^{\gamma }\Delta {\gamma }^{*}+a_{\varpi }^{\varpi }\Delta {\varpi }^{*}+a_{\varpi }^{\dot{\alpha }}\Delta {\dot{\gamma }}^{*}+b_{\varpi }^{\delta }\Delta {\delta }^{*}-a_{\varpi }^{\alpha }\Delta {\gamma }_W^{*} \end{gathered}$$

(15)

Applying the Laplace transform, the homogeneous form is:

$$\left[\begin{array}{cc}(1-a_{\gamma }^{\dot{\alpha }})s-a_{\gamma }^{\gamma }& -a_{\gamma }^{\varpi }\\ -(a_{\varpi }^{\dot{\alpha }}s+a_{\varpi }^{\gamma })& s-a_{\varpi }^{\varpi }\end{array}\right]\left[\begin{array}{c}\Delta {\gamma }^{\mathrm{*}}\\ \Delta {\varpi }^{\mathrm{*}}\end{array}\right]=\left[\begin{array}{c}0\\ 0\end{array}\right]$$

(16)

with the characteristic polynomial:

$$P\left(s\right)=\left(1-a_{\gamma }^{\dot{\alpha }}\right)s^2-\left(a_{\gamma }^{\gamma }+a_{\varpi }^{\varpi }+a_{\varpi }^{\dot{\alpha }}{a}_{\gamma }^{\varpi }-a_{\gamma }^{\dot{\alpha }}{a}_{\varpi }^{\varpi }\right)s+a_{\gamma }^{\gamma }{a}_{\varpi }^{\varpi }-a_{\gamma }^{\varpi }{a}_{\varpi }^{\gamma }$$

(17)

Separating complex terms, we will obtain:

$$\begin{gathered} P\left(s\right)=\left(1-a_{\gamma }^{\dot{\alpha }}\right)s^2-\left[a_{\gamma }^{\alpha }+a_{\varpi }^q+a_{\varpi }^{\dot{\alpha }}{a}_{\gamma }^q-a_{\gamma }^{\dot{\alpha }}{a}_{\varpi }^q+i\left(a_{\gamma }^{\beta }+a_{\varpi }^r+a_{\varpi }^{\dot{\alpha }}{a}_{\gamma }^r-a_{\gamma }^{\dot{\alpha }}{a}_{\varpi }^r\right)\right]s+a_{\gamma }^{\alpha }{a}_{\varpi }^q-a_{\gamma }^q{a}_{\varpi }^{\alpha } \\ +a_{\gamma }^r{a}_{\varpi }^{\beta }-a_{\gamma }^{\beta }{a}_{\varpi }^r+i\left(a_{\gamma }^{\alpha }{a}_{\varpi }^r-a_{\gamma }^q{a}_{\varpi }^{\beta }+a_{\varpi }^q{a}_{\gamma }^{\beta }-a_{\varpi }^{\alpha }{a}_{\gamma }^r\right). \end{gathered}$$

(18)

With the approximations (13), the characteristic polynomial becomes:

$$\begin{gathered} P\left(s\right)\cong s^2+{\displaystyle \frac{1}{t^{\mathrm{*}}}}\left\{{\displaystyle \frac{1}{i_B}}\left[\eta \left(-d_1\right)+\left({-d}_4\right)+\left({-d}_9\right)\right]+i{\displaystyle \frac{A}{B}}\widehat{p}\right\}s+{\displaystyle \frac{1}{t^{\mathrm{*}2}}}\left[{\displaystyle \frac{\left(-d_1\right)}{i_B}}+i{\displaystyle \frac{\left(-b_1\right)}{\tilde{m}}}{\displaystyle \frac{A}{B}}\widehat{p}-{\displaystyle \frac{\widehat{p}}{i_B}}{d}_{13}\right] \\ =s^2+{\displaystyle \frac{1}{i_B{t}^{\mathrm{*}}}}\left[\eta \left(-d_1\right)+\left({-d}_4\right)+\left({-d}_9\right)+i2i_B{\widehat{\omega }}_p\right]s+ \\ {\displaystyle \frac{1}{i_B{t}^{\mathrm{*}2}}}\left[\left(-d_1\right)+i\left(2\eta \left(-d_1\right){\widehat{\omega }}_p-\widehat{p}{d}_{13}\right)\right]. \end{gathered}$$

(19)

where $\eta$ will be explained later.

To verify the correctness of relations (13), in Figure 5 and Figure 6 the roots of the characteristic polynomial obtained with the approximate relation (19) and the exact relation (18) are presented comparatively. The diagrams are obtained for a range of Mach numbers: 0.3...0.45.

Based on these results, the approximate relations (13) will be used to obtain the simplified structural scheme for the commanded object.

4. The Flight Quality Parameters

For the considered configuration, similar to the work [11], the following notations are defined:

$${\omega }_p={\displaystyle \frac{A}{2B}}p; \varphi ={\displaystyle \frac{{\omega }_p}{\Omega }}; \eta ={\displaystyle \frac{i_B}{\zeta \tilde{m}}}; \chi =2\eta {\widehat{\omega }}_p-\widehat{p}{\displaystyle \frac{d_{13}}{\left(-d_1\right)}}; \zeta ={\displaystyle \frac{d_1}{b_1}},$$

(20)

where ${\omega }_p$is the angular rate of precession, with dimensionless form ${\widehat{\omega }}_p={\omega }_p{t}^{*}$ and $\zeta$ represents static stability (Figure B17).

In this case, the characteristic polynomial (19) has the form:

$$P\left(s\right)=s^2+2\Omega (\xi +i\varphi )s+{\Omega }^2(1+i\chi ),$$

(21)

with natural pulsation $\Omega$ and damping factor $\xi$ having the expressions:

$$\Omega =\sqrt{{\displaystyle \frac{a_{\gamma }^{\alpha }{a}_{\varpi }^q-a_{\gamma }^q{a}_{\varpi }^{\alpha }}{1-a_{\gamma }^{\dot{\alpha }}}}}\cong {\displaystyle \frac{1}{t^{\mathrm{*}}}}\sqrt{{\displaystyle \frac{\left(-d_1\right)}{i_B}}};$$

(22)

$$\xi ={\displaystyle \frac{\Omega \left(a_{\gamma }^{\alpha }+a_{\varpi }^q+a_{\varpi }^{\dot{\alpha }}{a}_{\gamma }^q-a_{\gamma }^{\dot{\alpha }}{a}_{\varpi }^q\right)}{2\left(a_{\gamma }^{\alpha }{a}_{\varpi }^q-a_{\gamma }^q{a}_{\varpi }^{\alpha }\right)}}\cong {\displaystyle \frac{1}{2\sqrt{i_B\left(-d_1\right)}}}\left[\eta \left(-d_1\right)+\left(-d_4\right)+\left(-d_9\right)\right]$$

(22)

Sometimes, the inverse of natural pulsation is used:

$T=1/\Omega$ - time constant of nutation movement.

In Figure B18, Figure B19 and Figure B21 from Appendix B the natural pulsation, damping factor and time constant according to the Mach number are presented.

The roots of the characteristic polynomial (21) are:

$$s_{\mathrm{1,2}}=-\xi \Omega -i{\omega }_p\pm i\Omega \sqrt{{1-\xi }^2+{\varphi }^2+i\left(\chi -2\xi \varphi \right)},$$

(23)

Developing we obtain:

$$s_{\mathrm{1,2}}=-\Omega \left(\xi +i\varphi \right)\mp \Omega \sigma \sqrt{1+K^2}\left(sin\varrho -icos\varrho \right)$$

(24)

where: $K={\displaystyle \frac{\left(\chi -2\xi \varphi \right)}{{1-\xi }^2+{\varphi }^2}}$; ${\sigma }^2=1-{\xi }^2+{\varphi }^2$; $\varrho ={\displaystyle \frac{{\mathrm{t}\mathrm{a}\mathrm{n}}^{-1}\left(K\right)}{2}}$To obtain the command parameter, we put the relation (15) in matrix form:

$$\left[\begin{array}{c}\Delta {\gamma }^{\mathrm{*}}\\ \Delta {\varpi }^{\mathrm{*}}\end{array}\right]={\mathit{A}}^{-1}\left(s\right)\left[\begin{array}{ccc}{b}_{\gamma }^{\delta }& -a_{\gamma }^{\alpha }& 1/\left(m{u}^{\mathrm{*}}\right)\\ b_{\varpi }^{\delta }& -a_{\varpi }^{\alpha }& 0\end{array}\right]{\left[\begin{array}{ccc}\Delta {\delta }^{\mathrm{*}}& \Delta {\gamma }_W^{\mathrm{*}}& \Delta F_p^{\mathrm{*}}\end{array}\right]}^T,$$

(25)

with the inverse matrix:

$${\mathit{A}}^{-1}\left(s\right)={\displaystyle \frac{1}{P\left(s\right)}}\left[\begin{array}{cc}s-a_{\varpi }^{\varpi }& a_{\gamma }^{\varpi }\\ a_{\varpi }^{\dot{\alpha }}s+a_{\varpi }^{\gamma }& \left(1-a_{\gamma }^{\dot{\gamma }}\right)s-a_{\gamma }^{\gamma }\end{array}\right],$$

(26)

where $P\left(s\right)$is the characteristic polynomial.

Based on this relation, as shown in [11], after processing, the following transfer function is obtained:

$$H_{\omega }^{\delta }\left(s\right)={\displaystyle \frac{-b_{\gamma }^{\delta }{s}^2+\left[\left(1-a_{\gamma }^{\varpi }-a_{\gamma }^{\dot{\alpha }}\right)b_{\varpi }^{\delta }+\left(a_{\varpi }^{\varpi }+a_{\varpi }^{\dot{\alpha }}\right)b_{\gamma }^{\delta }\right]s+a_{\varpi }^{\gamma }{b}_{\gamma }^{\delta }-a_{\gamma }^{\gamma }{b}_{\varpi }^{\delta }}{P\left(s\right)}}.$$

(27)

where the input ${\delta }^{\mathrm{*}}={\delta }_m^{\mathrm{*}}-i{\delta }_n^{\mathrm{*}}$ is the complex form of gasodynamic deflection in the Resal frame and the output ${\omega }^{\mathrm{*}}={\omega }_m-i{\omega }_n$ represents the complex form of the angular velocity of the velocity vector.

With the notations:

$$k_{\omega }^{\delta }={\displaystyle \frac{a_{\varpi }^{\alpha }{b}_{\gamma }^{\delta }-a_{\gamma }^{\alpha }{b}_{\varpi }^{\delta }}{a_{\gamma }^{\alpha }{a}_{\varpi }^q-a_{\gamma }^q{a}_{\varpi }^{\alpha }}}; 2\lambda T_{\omega }={\displaystyle \frac{a_{\varpi }^r{b}_{\gamma }^{\delta }-a_{\gamma }^r{b}_{\varpi }^{\delta }}{a_{\varpi }^{\alpha }{b}_{\gamma }^{\delta }-a_{\gamma }^{\alpha }{b}_{\varpi }^{\delta }}}; \nu ={\displaystyle \frac{a_{\varpi }^{\beta }{b}_{\gamma }^{\delta }-a_{\gamma }^{\beta }{b}_{\varpi }^{\delta }}{a_{\varpi }^{\alpha }{b}_{\gamma }^{\delta }-a_{\gamma }^{\alpha }{b}_{\varpi }^{\delta }}},$$

(28)

the transfer function, due to the large value of the parameter $\tilde{m}$, becomes:

$$H_{\omega }^{\delta }\left(s\right)\cong {\displaystyle \frac{k_{\omega }^{\delta }(1+i\nu +2i\lambda T_{\omega }s)}{T^2{s}^2+2T(\xi +i\varphi )s+1+i\chi }},$$

(29)

where the following notations were made:

$k_{\omega }^{\delta }$ the command factor and the parameter

$T_{\omega }$ the advance time on command or the aircraft's time constant [25]

and $\nu$ is a Magnus term.

Dropping the secondary terms for the gasodynamic commanded missile, the parameters from the relations (28) become:

$$\begin{gathered} k_{\omega }^{\delta }\cong -{\displaystyle \frac{{\widehat{x}}_T+\zeta }{t^{*}\tilde{m}\zeta }}{C}_T{\tilde{T}}_0=-{\displaystyle \frac{{\widehat{x}}_T+\zeta }{t^{*}\tilde{m}\zeta }}{\displaystyle \frac{T}{F_0}}=-{\displaystyle \frac{T\left(\zeta +{\widehat{x}}_T\right)}{mV\zeta }}\left[{\displaystyle \frac{1}{s}}\right]; \\ 2\lambda T_{\omega }\cong -{\displaystyle \frac{C_T{\tilde{T}}_0}{b_1}}{\displaystyle \frac{2i_B{t}^{*}{\widehat{\omega }}_p}{C_T{\tilde{T}}_0{\widehat{x}}_T+\zeta C_T{\tilde{T}}_0}}={\displaystyle \frac{2}{\left({-b}_1\right)}}{\displaystyle \frac{i_B{t}^{*}{\widehat{\omega }}_p}{\left(\zeta +{\widehat{x}}_T\right)}}={\displaystyle \frac{1}{\left(-b_1\right)}}{\displaystyle \frac{pA{H}_0^{-1}}{\left(\zeta +{\widehat{x}}_T\right)}}\left[s\right]; \\ \nu \cong {\displaystyle \frac{d_{13}{C}_T{\tilde{T}}_0}{b_1}}{\displaystyle \frac{\widehat{p}}{C_T{\tilde{T}}_0{\widehat{x}}_T+\zeta C_T{\tilde{T}}_0}}=-{\displaystyle \frac{d_{13}}{\left({-b}_1\right)}}{\displaystyle \frac{\widehat{p}}{\left(\zeta +{\widehat{x}}_T\right)}}. \end{gathered}$$

(30)

where we recall that it was noted: ${\widehat{x}}_T=x_T/l$ where $x_T$represents the distance between the point of application of the gasodynamic control and the center of mass, in the body frame, having a negative value.

We can also evaluate the complex acceleration in the Resal frame:

$$\Delta a^{*}\cong V\Delta {\omega }^{*},$$

(31)

with:

$$a^{*}=-\left(a_z^{*}+i{a}_y^{*}\right).$$

(32)

Based on these relations, we can build the structural scheme of the commanded movement for slow-rolling missiles, with gasodynamic command, shown in Figure 7.

5. The Guided Flight Model

5.1. The Guidance Relations

As shown in [21,26,27], the guidance kinematic relations will be constructed into a frame linked to the carrier - target line of sight called the guiding frame (Figure 8). The guiding frame $A{X}_T{Y}_T{Z}_T$ originates at point ($A)$, which represents the position of the carrier, the axis $X_T$ pointed target $\left(T\right)$, the axis $Z_T$ is in a vertical plane facing downwards, and the axis $Y_T$ horizontally, it follows from the condition of the right frame. The guiding frame is non-inertial.

In the case of angular deviation control, the guidance parameters are: $\begin{array}{cc}{\epsilon }_y;& {\epsilon }_z\end{array}$, for determining which relations are used:

$$\begin{array}{cc}{\epsilon }_y={\tan }^{-1}{\displaystyle \frac{z_M}{x_M}};& {\epsilon }_z=-{\tan }^{-1}{\displaystyle \frac{y_M}{x_M}};\end{array}$$

(33)

and:

$$\begin{array}{cc}{P}_L=\sqrt{x_M^2+z_M^2}& P_N=\sqrt{x_M^2+y_M^2}\end{array}$$

(34)

From Figure 9 and Figure 10 it is observed that the angular error in the two guidance planes can be expressed as follows:

$${\epsilon }_y={\sigma }_y-{\phi }_y;{\epsilon }_z={\sigma }_z-{\phi }_z.$$

(35)

Assuming that the target, the carrier and the mobile evolve independently in the two planes, the geometric connection between the angles can be derived, obtaining:

$${\upsilon }_y={\omega }_y-{\omega }_L;$$

(36)

$${\upsilon }_z={\omega }_z-{\omega }_N,$$

(36)

where the angular velocities were noted:

$$\begin{array}{ccc}{\upsilon }_y={\dot{\epsilon }}_y;& {\upsilon }_z={\dot{\epsilon }}_z;{\omega }_y={\dot{\sigma }}_y;{\omega }_z={\dot{\sigma }}_z;& {\omega }_L={\dot{\phi }}_y\end{array};{\omega }_N={\dot{\phi }}_z.$$

(37)

5.2. Linear Form of Kinematic Guidance Relations

In the case of semiautomatic command to line of sight (SACLOS), which is based on angular error control, the collinearity of the three points (A-M-T) is achieved by cancelling the angular deviation in the two guidance planes:

$$\begin{array}{cc}{\epsilon }_y=0;& {\epsilon }_z=0\end{array}.$$

(38)

Considering the geometric connections between angles in the guiding planes (35), one can write:

$$\begin{array}{cc}{\phi }_y={\sigma }_y;& {\phi }_z={\sigma }_z\end{array},$$

(39)

which represents the collinearity of the position vectors for target $\mathit{D}$ and missile $\mathit{P}$Figure 8.

If the condition of cancelling the rate of the deviations is also imposed:

$$\begin{array}{cc}{\upsilon }_y=0;& {\upsilon }_z=0\end{array},$$

(40)

relation (36) show:

$$\begin{array}{cc}{\omega }_y={\omega }_L;& {\omega }_z={\omega }_N\end{array},$$

(41)

which is the condition that the vectors $\mathit{D}$and $\mathit{P}$ rotate with the same angular velocity.

To construct a model of guidance with angular deviation control, we start from the observation that the linear shape of the guide controls is decoupled for the case of a general movement.

Thus, if the definition relation of the angular errors is linearized:

$${\epsilon }_y={\tan }^{-1}(z_M/x_M),{\epsilon }_z=-{\tan }^{-1}(y_M/x_M),$$

(42)

under the conditions of the basic movement, the following is obtained:

$$\Delta {\epsilon }_y=\Delta h_z/P,\Delta {\epsilon }_z=-\Delta h_y/P.$$

(43)

By entering the expression of the linear deviation previously established, we obtain:

$$\Delta {\epsilon }_y={\displaystyle \frac{\Delta h_{zM}}{P}}-{\displaystyle \frac{\Delta h_{zT}}{D}};\Delta {\epsilon }_z=-{\displaystyle \frac{\Delta h_{yM}}{P}}+{\displaystyle \frac{\Delta h_{yT}}{D}}.$$

(44)

On the other hand:

$$\Delta {\epsilon }_y=-\Delta {\phi }_y+\Delta {\sigma }_y; \Delta {\epsilon }_z=-\Delta {\phi }_z+\Delta {\sigma }_z$$

(45)

In this case, the deviation of the angular errors becomes:

$$\Delta {\epsilon }_y={\displaystyle \frac{\Delta h_{zM}}{P}}+\Delta {\sigma }_y={\displaystyle \frac{\Delta z_{Mo}}{P}}-{\displaystyle \frac{\Delta z_{Ao}}{P}}+\Delta {\sigma }_y;\Delta {\epsilon }_z=-{\displaystyle \frac{\Delta h_{yM}}{P}}+\Delta {\sigma }_y=-{\displaystyle \frac{\Delta y_{Mo}}{P}}+\Delta y_{Ao}/P+\Delta {\sigma }_z,$$

(46)

From whence:

$$\Delta {\epsilon }_y=-{\displaystyle \frac{u_M}{P{s}^2}}\Delta {\omega }_m-{\displaystyle \frac{\Delta z_{Ao}}{P}}+\Delta {\sigma }_y; \Delta {\epsilon }_z=-{\displaystyle \frac{u_M}{P{s}^2}}\Delta {\omega }_n+{\displaystyle \frac{\Delta y_{Ao}}{P}}+\Delta {\sigma }_z.$$

(47)

Considering the particularity of the basic movement adopted, in which the velocity vector of the mobile is collinear with the target's line of sight, namely:

$$V_M=u_M$$

(48)

and noting the input functions of order "0" due to the carrier's position deviations:

$$\Delta f_{oz}=-{\displaystyle \frac{\Delta z_{Ao}}{P}};\Delta f_{oy}=-{\displaystyle \frac{\Delta y_{Ao}}{P}},$$

(49)

the deviation of the angular error can be put in the form:

$$\Delta {\epsilon }_y=\Delta {\sigma }_y-\left({\displaystyle \frac{V_M}{P{s}^2}}\Delta {\omega }_m-\Delta f_{oz}\right);\Delta {\epsilon }_z=\Delta {\sigma }_z-\left({\displaystyle \frac{V_M}{P{s}^2}}\Delta {\omega }_n-\Delta f_{oy}\right).$$

(50)

thus, obtaining the cinematic block.

From this relation it can be seen that for the case of the fixed carrier the deviation of the absolute angle of the line of sight of the missile is given by:

$$\Delta {\phi }_y={\displaystyle \frac{V_M}{P{s}^2}}\Delta {\omega }_m-\Delta f_{oz};\Delta {\phi }_z={\displaystyle \frac{V_M}{P{s}^2}}\Delta {\omega }_n+\Delta f_{oy}.$$

(51)

On the other hand, the guidance commands in the guidance planes are:

$$\Delta u_y=(k_u^{\epsilon }\Delta {\epsilon }_y+k_u^{\upsilon }\Delta {\upsilon }_y)P, \Delta u_z=(k_u^{\epsilon }\Delta {\epsilon }_z+k_u^{\upsilon }\Delta {\upsilon }_z)P,$$

(52)

where:

$$\Delta {\dot{\epsilon }}_y=\Delta {\upsilon }_y; \Delta {\dot{\epsilon }}_z=\Delta {\upsilon }_z$$

(53)

By applying the Laplace transform, the guidance command can be put in the form:

$$\Delta u_y=P\left(k_u^{\upsilon }s+k_u^{\epsilon }\right)\Delta {\epsilon }_y; \Delta u_z=P\left(k_u^{\upsilon }s+k_u^{\epsilon }\right)\Delta {\epsilon }_z$$

(54)

For the linearization of the guidance command, a particular basic motion is considered, where the target, the mobile and the carrier are collinear and the velocity vector of the missile is oriented along this line.

In this case, the angle of lateral inclination of the missile must be zero, the head angle of the missile is equal to the angle of rotation in the horizontal plane of the guidance frame $\psi ={\psi }_A$, which overlaps the pitch plane of the missile with the vertical plane of guidance.

The axis of the rocket is above the velocity vector, with which it forms a positive incidence that balances the force of weight.

For the simplification of the structural scheme, in the case of the fixed carrier, the zero-order input functions $f_{0z};f_{0y}$ can be considered cumulated with the inputs due to the target through the angles ${\sigma }_y, {\sigma }_z$, obtaining in this case a simplified transfer function, with a single input.

5.3. The Guidance Command

As shown in the paper [11], neglecting the angle of incidence in the vertical plane, and assuming that the pitch angle of the rocket $\theta$ overlaps ${\theta }_A$ angle of the guiding frame, the command can be put in the following form:

$${\left[\begin{array}{cc}{u}_m& u_n\end{array}\right]}^T={\mathit{K}}_{\varphi }{\left[\begin{array}{cc}{u}_y& u_z\end{array}\right]}^T.$$

(55)

where the following notation was made:

$${\mathit{K}}_{\varphi }=\left[\begin{array}{ccc}0& 1& 0\\ 0& 0& 1\end{array}\right]{\mathit{A}}_{\varphi }\left[\begin{array}{cc}0& 0\\ 1& 0\\ 0& 1\end{array}\right]=\left[\begin{array}{cc}\mathit{cos}\varphi & \mathit{sin}\varphi \\ -\mathit{sin}\varphi & \mathit{cos}\varphi \end{array}\right]$$

(56)

${\mathit{A}}_{\varphi }$ being the direct rotation matrix with roll angle.

Next, we will seek to bring matrix expressions to complex form using a polar coordinate system that highlights the module and phase of the guidance command. Thus, if the commands are put in the complex form:

$$u=u_m-i{u}_n;u^{\mathrm{*}}=u_y-i{u}_z$$

(57)

where $u$ is the command in body frame and $u^{\mathrm{*}}$ is the command in the guiding frame / Resal frame.

In this case, the relation (55) becomes:

$$u=u^{\mathrm{*}}{e}^{i\varphi }$$

(58)

Base on Figure 11, the guidance command can be put in polar coordinates:

$$u^{\mathrm{*}}=\rho e^{i\phi }$$

(59)

where:

$\rho =\sqrt{u_y^2+u_z^2}$ – is the command module; $\phi =\mathrm{arctg}{\displaystyle \frac{{-u}_z}{u_y}}$ – is the command phase.

In this case, the command for the two channels (pitch and yaw) in the body frame becomes:

$$u=\rho e^{i(\phi +\varphi +\varsigma )},$$

(60)

where $\varsigma$ is the phase shift, the significance of this angular error being detailed in work [11].

As we showed in papers [9] and [11], for the single-channel case we can give up one component of the command in the body frame, obtaining in the guiding frame a command oriented in both directions due to the roll rotation.

To obtain the time modulation of the command for a single-channel missile, it is necessary to use a switching function of the form indicated in work [18]:

$$f(\tilde{\phi },x)=\mathit{sin}\mathsf{\Phi }(x+\mathit{cos}\mathsf{\Phi })$$

(61)

where $x$, called the fill factor, is given by the relation:

$$x=\rho /u_{max}$$

(62)

in which $\rho$ is the command module, $u_{max}$ is the maximum guidance command, and the argument $\mathsf{\Phi }$, called the relative roll angle, is given by:

$$\mathsf{\Phi }=\phi +\varphi +\pi /2+\varsigma ,$$

(63)

where $\phi$ represents the command phase (Figure 11), $\varphi$ is the roll angle and $\varsigma$ is the phase shift.

The switching function (61) is applied instantaneously through a nonlinear element, relay-type, the pitch canard deflection angle being:

$${\delta }_m(\mathsf{\Phi },x)={\delta }_{max}\mathrm{s}\mathrm{i}\mathrm{g}\mathrm{n}\left(f(\mathsf{\Phi },x)\right)$$

(64)

As shown in the [18] for the expression of the equivalent command, the weight compensation command is also entered. This is achieved by the equilibrium fill factor for a horizontal base motion, which is a function of the ratio of the weight compensation control component to the maximum order value:

$$k_e={\displaystyle \frac{2}{\pi }}{\sin }^{-1}{\displaystyle \frac{{\delta }_e^{\mathrm{*}}}{{\delta }_{max}^{\mathrm{*}}}}$$

(65)

Unlike the roll-stabilized missile, where the balance command can be achieved directly by changing the zero position of the command surfaces (canard, elevator), in the case of rolling missiles, the equivalent balance command is introduced by the guidance term; in this case, the balance control is non-zero. Resuming the linear relation established in the works [11] and [9] can be written as:

$$\Delta {\delta }^{*}=k_{\delta }^u\Delta u^{*},$$

(66)

which expresses the connection in deviations between the command and the thrust equivalent deflection.

Similarly, for the absolute angles of the target's line of sight, the quantity is defined:

$${\sigma }^{*}={\sigma }_y-i{\sigma }_z$$

(68)

and from the angular deviations in the guiding planes, the following are obtained:

$${\epsilon }^{*}={\epsilon }_y-i{\epsilon }_z$$

(69)

where:

$${\epsilon }_y={\sigma }_y-{\phi }_y; {\epsilon }_z={\sigma }_z-{\phi }_z$$

(70)

Between the three complex quantities defined above there is the relationship:

$${\epsilon }^{*}={\sigma }^{*}-{\phi }^{*}$$

(71)

Finally, it is also possible to define the complex form of the zero-order entry function in the system due exclusively to the operator or the sighting device:

$$f_0^{*}=f_{0z}+i{f}_{0y}$$

(72)

function that will be cumulated with the absolute angle of the target's line of sight (${\sigma }^{*}$).

Processing of the error signal and formation of the guidance control signal is done using a P.D. (derivative proportional) type element that has the following known structure:

$$H_u^{\epsilon }\left(s\right)={\displaystyle \frac{(k_u^{\epsilon }+s{k}_u^{\upsilon })}{\tau s+1}}P$$

(73)

where $k_u^{\epsilon } \left[m^{-1}\right]$and $k_u^{\upsilon }\left[m^{-1}s\right]$ are guidance gains, and $\tau$is a time constant that eliminates high frequencies.

NOTE - In case of Manual Command to Line of Sight (MCLOS), the control is done manually by a human operator, as indicated in the paper [18], and the signal can be shaped by a transfer function:

$$H_u^{\epsilon }\left(s\right)=e^{-{\tau }_{0}s}{\displaystyle \frac{(k_u^{\epsilon }+s{k}_u^{\upsilon })}{\tau s+1}}P$$

(74)

an expression which, as indicated in the same paper, can be approximated with:

$$H_u^{\epsilon }\left(s\right)={\left(1-s{\displaystyle \frac{{\tau }_0}{3}}\right)}^3{\displaystyle \frac{(k_u^{\epsilon }+s{k}_u^{\upsilon })}{\tau s+1}}P$$

(75)

5.4. The Structural Scheme

Next, in order to build the structural scheme of the SACLOS, it is necessary to define the complex shape of some pairs of kinematic parameters in the guiding planes.

Thus, for the absolute angles of the rocket's line of sight, the complex size can be defined:

$${\phi }^{*}={\phi }_y-i{\phi }_z$$

(67)

5.5. Determination of the Characteristic Polynomial

For the analysis of the stability, the transfer function of the open system is first determined:

$$H\left(s\right)={\displaystyle \frac{k\left(1+i\nu +2i\lambda T_{\omega }s\right)\left(1+{\tau }_{c}s\right)e^{i\varsigma }}{s^2\left(T^2{s}^2+2T(\xi +i\varphi )s+1+i\chi \right)\left)\right(1+\tau s)}},$$

(76)

where we considered null the phase shift introduced by the actuator $\varsigma =0$, and we noted the guidance gain:

$$k=k_u^{\epsilon }{k}_{\delta }^u{k}_{\omega }^{\delta }{V}_M\left[s^{-2}\right]$$

(77)

and the guidance time:

$${\tau }_c=k_u^{\nu }/k_u^{\epsilon }\left[s\right]$$

(78)

In this case, the transfer function of the closed system becomes:

$$H_0\left(s\right)={\displaystyle \frac{k{\Omega }^2(1+i\nu +2i\lambda T_{\omega }s)(1+{\tau }_{c}s)e^{i\varsigma }}{s^2\left(s^2+2\Omega (\xi +i\varphi )s+{\Omega }^2\left(1+i\chi \right)\right)(1+\tau s)+k{\Omega }^2(1+i\nu +2i\lambda T_{\omega }s)(1+{\tau }_{c}s)e^{i\varsigma }}},$$

(79)

where $\mathsf{\Omega }=1/\mathrm{T}$ represents one's pulsation.

For the stability check we develop the characteristic polynomial from the denominator of the transfer function.

The characteristic polynomial, after ordering, becomes:

$$\begin{gathered} P\left(s\right)=s^5+\left[\left({T\tau }^{-1}+2\xi \right)+2i\varphi \right]\Omega s^4+\left[\left(1+2\xi T{\tau }^{-1}\right)+i\left(\chi +2\varphi T{\tau }^{-1}\right)\right]{\Omega }^2{s}^3 \\ +\left[\left(1-2k\lambda T_{\omega }{\tau }_c\sin \left(\mathsf{\varsigma }\right)\right)+i\left(\chi +2k\lambda T_{\omega }{\tau }_c\cos \left(\mathsf{\varsigma }\right)\right)\right]{\Omega }^2{\tau }^{-1}{s}^2 \\ +k\left[\left(\mathrm{c}\mathrm{o}\mathrm{s}\right(\mathsf{\varsigma })-\mathrm{s}\mathrm{i}\mathrm{n}(\mathsf{\varsigma }\left)\right(\nu +2\lambda T_{\omega }{\tau }_c^{-1}\left)\right)+i\left(\mathrm{s}\mathrm{i}\mathrm{n}\right(\mathsf{\varsigma })+\mathrm{c}\mathrm{o}\mathrm{s} \\ (\mathsf{\varsigma }\left)\right(\nu +2\lambda T_{\omega }{\tau }_c^{-1}\left)\right)\right]{\tau }_c{\tau }^{-1}{\Omega }^{2}s \\ +k\left[\left(\mathrm{c}\mathrm{o}\mathrm{s}\right(\mathsf{\varsigma })-\nu \mathrm{s}\mathrm{i}\mathrm{n}(\mathsf{\varsigma }\left)\right)+i\left(\mathrm{s}\mathrm{i}\mathrm{n}\right(\mathsf{\varsigma })+\nu \mathrm{c}\mathrm{o}\mathrm{s}(\mathsf{\varsigma }\left)\right)\right]{\Omega }^2{\tau }^{-1} \end{gathered}$$

(80)

For the stability analysis, we put the polynomial in the form:

$$P\left(s\right)=s^5+\left(p_1+i{q}_1\right)s^4+\left(p_2+i{q}_2\right)s^3+\left(p_3+i{q}_3\right)s^2+\left(p_4+i{q}_4\right)s+\left(p_5+i{q}_5\right)$$

(81)

where:

$$\begin{gathered} p_1=\left(T{\tau }^{-1}+2\xi \right)\Omega ;q_1=2\varphi \Omega \\ {p}_2=\left(1+2\xi {T\tau }^{-1}\right){\Omega }^2;q_2=\left(\chi +2\varphi T{\tau }^{-1}\right){\Omega }^2 \\ {p}_3=\left[1-2k\lambda T_{\omega }{\tau }_c\mathrm{s}\mathrm{i}\mathrm{n}\left(\mathsf{\varsigma }\right)\right]{\Omega }^2{\tau }^{-1};q_3=\left[\chi +2k\lambda T_{\omega }{\tau }_c\mathrm{c}\mathrm{o}\mathrm{s}\left(\mathsf{\varsigma }\right)\right]{\Omega }^2{\tau }^{-1}; \\ {p}_4=k\left[\cos \left(\mathsf{\varsigma }\right)-\sin \left(\mathsf{\varsigma }\right)\left(\nu +2\lambda T_{\omega }{\tau }_c^{-1}\right)\right]{\Omega }^2{\tau }_c{\tau }^{-1}; \\ {q}_4=k\left[\mathrm{s}\mathrm{i}\mathrm{n}\left(\mathsf{\varsigma }\right)+\mathrm{c}\mathrm{o}\mathrm{s}\left(\mathsf{\varsigma }\right)\left(\nu +2\lambda T_{\omega }{\tau }_c^{-1}\right)\right]{\Omega }^2{\tau }_c{\tau }^{-1}; \\ {p}_5=k\left[\left(\mathrm{c}\mathrm{o}\mathrm{s}\right(\mathsf{\varsigma })-\nu \mathrm{s}\mathrm{i}\mathrm{n}(\mathsf{\varsigma })\right]{\Omega }^2{\tau }^{-1};q_5=k\left[\mathrm{s}\mathrm{i}\mathrm{n}\left(\mathsf{\varsigma }\right)+\nu \mathrm{c}\mathrm{o}\mathrm{s}\left(\mathsf{\varsigma }\right)\right]{\Omega }^2{\tau }^{-1}. \end{gathered}$$

(82)

5.6. The Guidance System Synthesis

At this point we aim to determine the three parameters of the guidance system:$k,{\tau }_c, \tau$.

To begin with, we will neglect the complex coefficients, considering a regular missile without roll rotation. In this case, the transfer function of the closed system (79) becomes:

$$H_0\left(s\right)={\displaystyle \frac{k{\tau }_c{{\tau }^{-1}\Omega }^{2}s+k{\tau }^{-1}{\Omega }^2}{s^5+\left({\tau }^{-1}+2\Omega \xi \right)s^4+\left({\Omega }^2+2\Omega \xi {\tau }^{-1}\right)s^3+{{\Omega }^2{\tau }^{-1}s}^2+k{\tau }_c{\tau }^{-1}{\Omega }^{2}s+{k\tau }^{-1}{\Omega }^2}}$$

(83)

Next, we will use an allocation of poly-zeros according to the methodology “Standard coefficient method” from the paper [28]. This analysis can give us indicative information on the main control parameters of the system.

According to the methodology, for the transfer function obtained (84), with a characteristic polynomial of order 5, the recommended transfer function is of the form:

$$H\left(s\right)={\displaystyle \frac{18{\Omega }_0^{4}s+{\Omega }_0^5}{s^5+9{\Omega }_0{s}^4+29{\Omega }_0^2{s}^3+38{\Omega }_0^3{s}^2+18{\Omega }_0^{4}s+{\Omega }_0^5}},$$

(84)

where ${\Omega }_0={\displaystyle \frac{9}{t_r}}$ with $t_r$ the desired system response time.

By identifying the coefficients of the characteristic polynomial of the two transfer functions, the following relations are obtained:

$${k\tau }^{-1}{\Omega }^2={\Omega }_0^5; k{\tau }_c{\tau }^{-1}{\Omega }^2=18{\Omega }_0^4;{\Omega }^2{\tau }^{-1}=38{\Omega }_0^3; {\Omega }^2+2\Omega \xi {\tau }^{-1}=29{\Omega }_0^2; {\tau }^{-1}+2\Omega \xi =9{\Omega }_0.$$

(85)

The main quantities that intervene in these relations and that are related to the dynamics of the missile are:

$\Omega$ – pulsation; $\xi$ – damping factor and $k_{\omega }^{\delta }$ – amplification factor that intervenes in the gain $k$.

The values of these quantities are presented graphically according to the Mach number in the initial phase (with fuel) and in the final phase (without fuel) in Annex B of the work

From the diagrams (Annex B), for the analyzed velocity range $V_M\cong 100\dots .150$, the parameters have the values: $\Omega =10\dots 16.$; $\xi =0.16\dots 0.19$; $k_{\omega }^{\delta }\cong 0.5\dots 0.45$;

From the last two relations (85) is obtained:

$${\Omega }^2+2\Omega \xi {\tau }^{-1}=29{\Omega }_0^2$$

(86)

$$; {\tau }^{-1}=\left(9{\Omega }_0-2\Omega \xi \right)$$

(86)

Eliminating ${\tau }^{-1}$, the equation follows:

$$29{\Omega }_0^2-18\Omega \xi {\Omega }_0-{\Omega }^2\left(1-4{\xi }^2\right)=0$$

(87)

with the positive solution:

$${\Omega }_0={\displaystyle \frac{\Omega }{29}}\left(9\xi +\sqrt{29-35{\xi }^2}\right)$$

(88)

For the range of values of pulsation $\Omega$and damping factor $\xi$, the solution has the following values: ${\Omega }_0=2.32$...3.85

From ${\Omega }_0$, the response time can be estimated:

$$t_r={\displaystyle \frac{9}{{\Omega }_0}}={\displaystyle \frac{9}{2.32\dots 3.85}}=3.87\dots .2.34\left[\mathrm{s}\right]$$

Since the response time of the system must be low enough to combat close targets, we choose $t_r=3$, from which it follows ${\Omega }_0=3$.

From the third relation (85) we get:

$$\tau ={\displaystyle \frac{{\Omega }^2}{38{\Omega }_0^3}}$$

(89)

which for the pulsation values range $\Omega =10\dots 16.$, has values in the range of:

$$\tau =0.1\dots 0.25 \left[\mathrm{s}\right]$$

Combining the second relation (85) with the relation (89) yields:

$${\tau }_c={\displaystyle \frac{18{\Omega }_0}{38k}}$$

(90)

for which an estimate of the gain $k$ is needed.

Considering approximate relations (30) we obtain:

$$k=k_u^{\epsilon }{k}_{\delta }^u{k}_{\omega }^{\delta }{V}_M\cong -k_u^{\epsilon }{k}_{\delta }^u{\displaystyle \frac{T\left(\zeta +{\widehat{x}}_T\right)}{m{V}_M\zeta }}{V}_M=-k_u^{\epsilon }{k}_{\delta }^{u}T{\displaystyle \frac{\left(\zeta +{\widehat{x}}_T\right)}{m\zeta }}$$

(91)

The parameters that intervene in the expression of gain are $k_u^{\epsilon }=0.035\left[{\mathrm{m}}^{-1}\right]$; $k_{\delta }^u={\displaystyle \frac{2}{\pi }}\cong 0.64$; $T\cong 90N$; $\zeta \cong 0.0314\dots 0.0385$; ${\widehat{x}}_T=-0.295\dots -0.313;m=10.63\dots 8.63kg$from which it follows: $k\cong -2{\displaystyle \frac{\left(\zeta +{\widehat{x}}_T\right)}{m\zeta }}\cong 1.59\dots 1.67\left[s^{-2}\right]$;

For this range of variation of the gain $k,$ from (90) we obtain:

$$k\cong 1.63 \left[s^{-2}\right]; \tau =0.2 \left[s\right]$$

$$; {\tau }_c=0.87 \left[s\right]$$

Based on the results obtained, the values adopted are:

Based on the relations thus obtained, the structural scheme for the SACLOS single-channel rolling missile controlling the angular deviation in the case of the fixed carrier are presented in Figure 12

Figure 12. Structural scheme for rolling single-channel SACLOS missile.

Figure 12. Structural scheme for rolling single-channel SACLOS missile.

6. Calculus Model

To obtain numerical results, we will further specify the data of the calculus model.

6.1. Aerodynamic Characteristics

The aerodynamic characteristics were determined using the configuration from Figure 13.

Considering the specifics of the configuration and the flight regimes, starting from the regular aerodynamic coefficients indicated in [29] and [30], according [31] we introduce additional terms that contain influence of rotational and non-stationary motion and combinations with incidences.

Thus, the development of aerodynamic coefficients in the body frame for the slow-rolling airframe from Figure 13 is given by:

$$\begin{gathered} C_x=a_1+a_2\left({\alpha }^2+{\beta }^2\right)+a_3\left({\alpha }^4+{\beta }^4\right)+a_4{\alpha }^2{\beta }^2 \\ {C}_y={-b}_1\beta -b_2{\beta }^3-b_4\widehat{r}-b_9\widehat{\beta }+b_{13}\widehat{p}\alpha +b_{14}\widehat{p}\widehat{q} \\ {C}_z=b_1\alpha +b_2{\alpha }^3+b_4\widehat{q}{+b}_9\widehat{\alpha }+b_{13}\widehat{p}\beta +b_{14}\widehat{p}\widehat{r} \\ {C}_l=c_1+c_3\widehat{p} \\ {C}_m=d_1\alpha +d_2{\alpha }^3+d_4\widehat{q}{+d}_9\widehat{\alpha }+d_{13}\widehat{p}\beta +d_{14}\widehat{p}\widehat{r} \\ {C}_n=d_1\beta +d_2{\beta }^3+d_4\widehat{r}+d_9\widehat{\beta }-d_{13}\widehat{p}\alpha -d_{14}\widehat{p}\widehat{q} \end{gathered}$$

(92)

where, the rotational velocities and non-stationary translational variables are considered dimensionless, according to the standards [22,31] are:

$$\widehat{p}=p{t}^{\mathrm{*}};\widehat{q}=q{t}^{\mathrm{*}};\widehat{r}=r{t}^{\mathrm{*}};\widehat{\alpha }=\dot{\alpha }{t}^{\mathrm{*}};\widehat{\beta }=\dot{\beta }{t}^{\mathrm{*}}.$$

(93)

The terms of development of aerodynamic coefficients are presented against the Mach number in Appendix A.

6.2. Mechanical and Reference Characteristics

The mechanical characteristics of the model are indicated in Table 1

Table 1. Mechanical characteristics.

Phase	Mass m[kg]	Center of Mass ${\mathrm{x}}_{\mathrm{cm}}$[m]	Roll Inertial Moment A [$\mathrm{kg}{\mathrm{m}}^2$]	Yaw/Pitch Inertial Moment B $\left[\mathrm{kg}{\mathrm{m}}^2\right]$	Time [s]
Start	11.4	0.49	0.021	0.373	0
Initial (1	10.6	0.50	0.018	0.364	0.65
Final	8.6	0.49	0.016	0.341	27.

¹After finishing the booster.

Table 1. Mechanical characteristics.

Phase	Mass m[kg]	Center of Mass ${\mathrm{x}}_{\mathrm{cm}}$[m]	Roll Inertial Moment A [$\mathrm{kg}{\mathrm{m}}^2$]	Yaw/Pitch Inertial Moment B $\left[\mathrm{kg}{\mathrm{m}}^2\right]$	Time [s]
Start	11.4	0.49	0.021	0.373	0
Initial (1	10.6	0.50	0.018	0.364	0.65
Final	8.6	0.49	0.016	0.341	27.

¹After finishing the booster.

For the considered configuration, the reference length is $l=0.830 m, \mathrm{the}$ reference surface is $S=0.012 m^2$, and the maximum angular thrust deflection is ${\delta }_{max}=17°$.

6.3. Thrust Characteristic

The thrust characteristic of the model is shown in Figure 14.

Given the thrust diagram and the table of mechanical characteristics, the missile will have two flight phases: a first phase of acceleration (boosting) and a second phase of marching flight, with approximately constant velocity, as shown in the velocity diagram in Figure 33 and Figure 34.

6.4. Time Constants and Controller Gains

As basic data for the linear model, the following parameters were considered:

• Guidance gain $k\cong 1.63$ indicated as the optimal value in the previous analysis
• The guidance time constant for the derivative term of the guidance law ${\tau }_c=0.87\left[s\right]$
• The time constant for guidance control $\tau =0.2 \left[s\right]$
• Rolling rotational velocity close to that indicated in Figure 35 and Figure 36 $p=\dot{\varphi }\cong 50 \left[{\displaystyle \frac{rad}{s}}\right]\cong 8 \left[{\displaystyle \frac{rot}{s}}\right]$
• Calculation altitude $H=0 \left[km\right]$
• Phase shift $\varrho =20°$

7. Stability Analysis of the Single-Channel SACLOS Slow-Rolling Missile

7.1. Stability Analysis with the Routh–Hurwitz Criterion

Further, neglecting complex terms, we will do a stability analysis with the Routh-Hurwitz criterion to highlight restrictions on the guidance and control parameters.

The Routh–Hurwitz conditions for the 5^th order polynomial with real coefficients are obtained from the matrix of coefficients constructed as follows:

$$\left|\begin{array}{ccccc}{p}_1& p_3& p_5& 0& 0\\ 1& p_2& p_4& 0& 0\\ 0& p_1& p_3& p_5& 0\\ 0& 1& p_2& p_4& 0\\ 0& 0& p_1& p_3& p_5\end{array}\right|$$

(94)

From condition: ${∆}_1=p_1>0,$ results successively:

$$p_1=\Omega \left(T{\tau }^{-1}+2\xi \right)>0=>-2\xi \Omega {<\tau }^{-1}=>$$

(95)

$$-{\left(2\xi \Omega \right)}^{-1}<\tau$$

(95)

From condition: ${∆}_2={p_{1}p}_2-p_3>0$ is obtained successively:

$$\Omega \left(T{\tau }^{-1}+2\xi \right)\left(1+2T\xi {\tau }^{-1}\right)-{\tau }^{-1}>0$$

(96)

$${\tau }^2+2T\xi \tau +T^2>0$$

(97)

For stabile missile, with $0<\xi <1;$and $T>0$ relation (97) is always satisfied, being quadratic equations in$\tau$with complex conjugate roots.

From the condition ${∆}_3>0$, results successively:

$$\begin{gathered} {∆}_3={p_{1}p}_2{p}_3-p_3^2+{p_{1}p}_5{{-p}_{4}p}_1^2 \\ {∆}_3=\left[-T\left(T+2\xi \tau \right)\left({T\tau }_c-T\tau +2\xi \tau {\tau }_c\right)k+2\xi \left(T^2+2\xi T\tau +{\tau }^2\right)\right]T^{-5}{\tau }^{-3}>0 \end{gathered}$$

(98)

Since ${∆}_3$is linear in$k$, condition (98) becomes:

$$k<k_3={\displaystyle \frac{2\xi \left(T^2+2\xi T\tau +{\tau }^2\right)}{T\left(T+2\xi \tau \right)\left({T\tau }_c-T\tau +2\xi \tau {\tau }_c\right)}}$$

(99)

relation that allows obtaining the stability domain related to ${∆}_3$for $\tau$fixed.

From condition ${∆}_4>0$ results successively:

$$\begin{gathered} {∆}_4={p_{1}p}_2{p}_3{p}_4-p_3^2{p}_4+{p_{1}p}_4{p}_5-p_1^2{p}_4^2-p_1{p}_2^2{p}_5+{p_{1}p}_4{p}_5+{p_{2}p}_3{p}_5-p_5^2>0 \\ {∆}_4=\left[-T{\left({T\tau }_c-T\tau +2\xi \tau {\tau }_c\right)}^2{k}^2+2\xi \left(T^2+2\xi T\tau +{\tau }^2\right)\left({\tau }_c-\tau -2\xi T\right)k\right]T^{-7}{\tau }^{-4}>0 \end{gathered}$$

(100)

Since ${∆}_4$is quadratic in $k$, the second condition (100) becomes:

$$0<k<k_4={\displaystyle \frac{2\xi \left(T^2+2\xi T\tau +{\tau }^2\right)\left({\tau }_c-\tau -2\xi T\right)}{T{\left({T\tau }_c-T\tau +2\xi \tau {\tau }_c\right)}^2}}$$

(101)

The relations (99) and (101) allow obtaining the stability domains related to ${∆}_3$and ${∆}_4$.

Figure 15 and Figure 16 show the corresponding stability domains $k=k\left({\tau }_c\right)$ for ${∆}_3$and ${∆}_4$, obtained with relations (99) and (101), for $\tau$fixed at operating point ($\tau =0.2),$ in the initial and final of the second phase.

It can be seen from the two diagrams that the stability domain for ${∆}_4$is more restrictive than that for ${∆}_3$, and hence only the domain corresponding to ${∆}_4$ will be taken into consideration.

In addition, the domain ${∆}_4$ highlights the inner limit of ${\tau }_c$ :

$${\tau +2\xi T<\tau }_c,$$

(102)

where the gain $k$ becomes zero.

Applying relations (99) and (101) over the entire Mach number range, the diagrams from Figure 17 and Figure 18 are obtained. It can be seen from these diagrams that the $k4$ limit is always more restrictive than $k3$ and that the current value of the gains $k$ falls within the stability domain over the entire Mach number range.

To summarize, for a regular system with a subunitary damping factor (0$<\xi <1$) and the time constant of the nutation positive $(T>0)$, the only constraints that remain active are those given by ${∆}_4>0$ which limits the guidance gain $\left(0<k<k_4\right)$ and lower limits the guidance time ${(\tau +2\xi T<\tau }_c),$ and $p_5$which lower limits the time constant ($0<\tau$).

The result obtained is equivalent to the criterion Liénard–Chipart, which for the 5^th order polynomial with real coefficients imposes the following conditions [32]:

$$p_1>0 ;p_3>0;p_5>0;{∆}_2>0;$$

(103)

$${∆}_4>0$$

(103)

As an additional check, in Figure D1 and Figure D2, the roots of the characteristic polynomial with real coefficients at the operating point, analyzed by the conditions R-H, are presented. It is observed that the polynomial has four conjugate complex roots with negative real parts and a negative real root, which confirms the stability of the analyzed system.

Next, we will analyze the characteristic polynomial of the system, considering the complex part of its terms using the stability criterion Frank-Wall.

7.2. Stability Analysis with the Frank-Wall Criterion

Based on the relations (E 5), (E 6), (E 7) and the characteristic polynomial coefficients, the Frank-Wall parameters ($A_1,A_2,A_3,A_4,A_5)$ were determined and graphically shown in Appendix C. For the basic input data considered, the stability criterion F-W is met: ${\mathrm{A}}_1>0;{\mathrm{A}}_2>0;{\mathrm{A}}_3>0;{\mathrm{A}}_4>0;{\mathrm{A}}_5>0$, as we can see from Figure C9, Figure C10, Figure C11 and Figure C12. It is observed that they are positive over the entire domain, which means that the roots of the characteristic polynomial will have real negative parts.

Next, we will analyze the root locus of the characteristic polynomial (Appendix D) To begin, we shall consider the basic test case set out in section 7.3, Table 1, for the second phase of flight at two points, at the beginning ($t=0.65 \left[s\right]$) and the end of the phase ($t=27 \left[s\right]$). The analysis is done on the range of Mach numbers corresponding to phase 2 ($Mach=[0.3,\dots , 0.45]$ ) at ground level. A first observation of a general nature is that the roots do not have a symmetrical distribution with respect to the real axis, as in the case of stabilized roll rockets, where the characteristic polynomial has real coefficients. There are five roots of the polynomial, thus:

- Roots 2 and 4 are close to those of the characteristic polynomial of the commanded object (23) or (24)(24) which has a negative real part and large complex parts that are almost conjugate.

- Root 3 is large in module, with a negative real part and a small complex part, due to guidance control.

- Roots 1 and 5 are small in module, with negative real parts and small complex parts due to the guidance loop

Figure D3 and Figure D4 present the roots of the characteristic polynomial with complex coefficients for the initial point of the second phase without phase shift $(\varsigma =0°$) and with phase shift ($\varsigma =20°)$.

Figure D5 and Figure D6 present the roots of the characteristic polynomial with complex coefficients for the final point of the second phase, without phase shift ($\varsigma =0°$) and with a phase shift $\left(\varsigma =20°\right).$In this case, unlike the case with real coefficients previously analyzed (Figure D1 and Figure D2 ), we have five complex roots arranged non-symmetrically related to the real axis, that is, we no longer have complex conjugate roots. However, it is observed that the real part of the roots is negative in all cases, which ensures the stability of the analyzed system.

It is of interest to verify whether, in the case of the missile with rolling rotation (polynomial with complex coefficients), we obtain the same stability domains as in the case of the R-H analysis for a missile without rolling (polynomial with real coefficients).

Unfortunately, because the F-W calculation relations are complicated, we cannot find an analytic solution like in the R-H case. Therefore, we will summarize some numerical results for verification.

Thus, in Figure 19, Figure 20, Figure 21 and Figure 22 the variation of the F-W parameters is presented depending on the guidance gain $\left(k\right)$. It is observed that the cancellation of these parameters ($A3,A4,A5$), that is, the limit of the stability domain, is done by a jump from the positive to the negative value. The jump is explained by the fact that these points are close to the points of cancellation of the R-H determinants and by the relation (E 8), between the determinants R-H and the parameters F-W, in which the determinant R-H is in the denominator of the parameter F-W. Due to the jump at the point of cancellation, determining an exact solution numerically for these points is difficult.

From the analysis of the diagrams, it is observed that for $A_3$ the condition is similar to ${∆}_3$ is found, i.e. $k < k_{3W}$, and for $A_4$ the stability conditions similar to ${∆}_4$ are found, i.e. we have a domain for $k$ between $0$ and $k_{4w}$. The condition $k_{4w}< k_{3W}$ is also found, i.e. the stability condition for $A_4$ is more restrictive than $A_3$. As for $A_5$, this introduces a condition similar to $A_3$. Regarding $A_1$ and $A_2$, these parameters were not analyzed because they do not depend on the gain $k$.

It can be seen from the figures that in the initial phase (Figure 19,Figure 20) the cancellation values $k_{4W}$ are smaller and in the final phase (Figure 21, Figure 22), which means that the domain of stability is more restricted in the initial phase. Also, in the initial phase, a third cancellation point $k_{4W}$ is observed, which is not found in the final phase nor in the R-H analysis. Regarding the shift phase, it is observed that this leads to a decrease in the upper cancellation point $k_{4W}$, therefore there is a slight decrease in the stability domain.

It can be seen from the two diagrams that the stability domain for ${\mathrm{A}}_4$is more restrictive than that for ${\mathrm{A}}_3$, and hence only the domain corresponding to ${\mathrm{A}}_4$ will be taken into consideration.

In addition, the domain ${\mathrm{A}}_4$ highlights the inner limit of ${\tau }_c$ from relation (102).

The stability domains are presented for the initial moment (Mach=0.3) and the final moment (Mach=0.45). Similar to the R-H analysis, it is observed that in the case of the F-W analysis, the stability domain is smaller in the initial phase than in the final phase of the guided flight.

Figure 23. Stability F-W domains Mach=0.3;$\varsigma =0°;\tau =0.2$- initial.

Figure 23. Stability F-W domains Mach=0.3;$\varsigma =0°;\tau =0.2$- initial.

Figure 24. Stability F-W domains Mach=0.45, $;\varsigma =0°$ ; $\tau =0.2$– final.

Figure 24. Stability F-W domains Mach=0.45, $;\varsigma =0°$ ; $\tau =0.2$– final.

From the previous analysis it follows that the most restrictive stability domain F-W, which must be considered, depends on the first point of cancellation of $A_4$, namely $k_{4w}.$To have a comparative picture, in Figure 25 and Figure 26, the stability range ${∆}_4$ obtained with the R-H criterion and the stability range $A_4$ obtained with the F-W criterion are presented together.

Thus, in Figure 25 and Figure 26 the stability range of the gain $k$ as a function of the guiding time ${\tau }_c$ for a fixed time constant$\tau$. It is also observed from all the diagrams, that for the considered operating range, the F-W stability domain is reduced compared to the R-H domain, which represents an important result in the comparative analysis performed.

To exemplify the result of the previous analysis, we will start from the operating point and increase the gain $k$ until we obtain instability on the F-W criterion. Then we will check whether stability is maintained for the R-H criterion at that new operating point.

Thus, if we increase $k_u^{\epsilon }$ from $k_u^{\epsilon }=0.035$ to $k_u^{\epsilon }=0.105$, and the other parameters remain unchanged (${\tau }_c=0.87$, $\tau =0.2$), the F-W parameter that will become negative will be $A_4$, as can be seen from Figure 27.

At the same time, one root of the characteristic polynomial with complex coefficients has a positive real part, as seen in Figure 28. The diagram is obtained for a range of Mach numbers 0.3...0.45.

On the other hand, from Figure 29 it can be seen that for $\mathrm{M}\mathrm{a}\mathrm{c}\mathrm{h}> 0.35$ the gain k value is maintained in the stability range (R-H) determined on the polynomial with real coefficients but is out of the stability range (F-W) determined on the polynomial with complex coefficients.

Moreover, for this case, the roots of the characteristic polynomial with real coefficients have a negative real part, as can be seen from Figure 30. The diagram is obtained for a range of Mach numbers: 0.3...0.45.

This means that for the case of a rolling missile, the stability analysis on the polynomial with real coefficients, in this case R-H, is not sufficient, requiring an additional analysis on the polynomial with complex coefficients of the F-W type and an evaluation on the nonlinear model, as we will proceed further.

7.3. Stability Analysis Based on the Nonlinear Model

Based on equations of motion in the body frame [11], and kinematic guidance relations (33), (34), as well as guidance control equations (73), switching function (61) with the relative roll angle definition relation (63) and expression of pitch canard deflection angle (64), the nonlinear model of guided flight of the single-channel, slow-rolling, SACLOS missile is obtained.

A nonlinear model in the Resal frame, based on relations from the paper [11] would be much easier to build and use, but we would lose essential information related to command formation such as the fill factor or the command phase, parameters useful in the analysis as will be seen below.

The objective of the nonlinear model analysis is to verify the limits of the guidance gain ($k)$ set out in the previous point, such as the guidance time $\left({\tau }_c\right).$ We will consider parameters similar to linear analysis: $\tau =0.2 \left[s\right]; \varsigma =0° or \varsigma =20°$and a tactical situation with a fixed target placed at a distance of 3 km and a height of 10 m. For the gain we will use the relation (91):

$$k=k_u^{\epsilon }{k_{\delta }^{u}k}_{\omega }^{\delta }{V}_M$$

(104)

From this relation, only the first factor $k_u^{\epsilon }$ can be imposed because it is related to the constructive solution guidance device. At the same time, the second factor is constant $k_{\delta }^u=2/\pi$ and the third factor $k_{\omega }^{\delta }$ depends on the dynamics of the missile, being difficult to control. In this case, we will choose $k_u^{\epsilon }$ so that, when multiplied by the other two gains and velocity, a value close to the desired $k$ is obtained. Given the objectives of this section of the paper, we will seek to obtain three gains, the first less than the operating value $k\cong 0.7$ the second close to the operating value $k\cong 1.6$, considerate optimal, and used in linear analysis, and the third bigger than the operating value, at the limit of the stability F-W range, as can be seen from Figure 29, namely $k\cong 6$. As previously stated, these values will be determined as a function of time, simultaneously with the integration of nonlinear equations of motion that will be done in the body frame. Obviously, due to the variation of the third factor $k_{\omega }^{\delta }$, and the velocity $V_M$ the guidance gain will vary over time, as can be seen in $\mathrm{F}\mathrm{i}\mathrm{g}\mathrm{u}\mathrm{r}\mathrm{e}31$and Figure 32 where the phase shift was $\varsigma =0° \mathrm{o}\mathrm{r} \varsigma =20°$. The analysis is performed considering three identical missiles but with different guidance gains, the first: M1, below the optimal value (the under-gain case), the second: M2 close to the optimal value defined in the linear analysis (the optimal-gain case), and the third: M2 bigger than the optimal (the over-gain case).

Figure 31. Guidance gain; $\varsigma =0°$; $\tau =0.2$; ${\tau }_c=0.87.$

Figure 31. Guidance gain; $\varsigma =0°$; $\tau =0.2$; ${\tau }_c=0.87.$

Figure 32. Guidance gain; $\varsigma =20°$; $\tau =0.2$; ${\tau }_c=0.87.$

Figure 32. Guidance gain; $\varsigma =20°$; $\tau =0.2$; ${\tau }_c=0.87.$

For the three rocket models defined above Figure 33 and Figure 34 show the velocity diagram that was previously used to define the analysis range on the linear model $\mathrm{M}\mathrm{a}\mathrm{c}\mathrm{h}=[0.3,\dots ,0.45]$, which corresponds to the second phase of flight. There are no significant differences between M1 and M2. For M3 (over-gain case), a loss of velocity is observed. This phenomenon occurs regardless of the phase shift, as can be seen from the two diagrams.

Figure 33. Velocity diagram; $\varsigma =0°$; $\tau =0.2$; ${\tau }_c=0.87.$

Figure 33. Velocity diagram; $\varsigma =0°$; $\tau =0.2$; ${\tau }_c=0.87.$

Figure 34. Velocity diagram; $\varsigma =20°$; $\tau =0.2$; ${\tau }_c=0.87.$

Figure 34. Velocity diagram; $\varsigma =20°$; $\tau =0.2$; ${\tau }_c=0.87.$

Figure 35 and Figure 36 show the angular roll velocity diagram based on which the rotational velocity of $p\cong 8\left[{\displaystyle \frac{\mathrm{r}\mathrm{o}\mathrm{t}}{\mathrm{s}}}\right]$, corresponding to the second phase of flight, used in linear analysis, was chosen. It is observed that the roll velocity follows the velocity profile and does not differ substantially between M1 and M2. For M3 (over-gain case) a loss of roll velocity is observed. Similarly, for roll velocity, the phase shift does not have a significant influence.

Figure 35. Roll velocity diagram; $\varsigma =0°$; $\tau =0.2$; ${\tau }_c=0.87.$

Figure 35. Roll velocity diagram; $\varsigma =0°$; $\tau =0.2$; ${\tau }_c=0.87.$

Figure 36. Roll velocity diagram; $\varsigma =20°$; $\tau =0.2$; ${\tau }_c=0.87.$

Figure 36. Roll velocity diagram; $\varsigma =20°$; $\tau =0.2$; ${\tau }_c=0.87.$

Figure 37 and Figure 38 show the vertical and horizontal projections of the trajectories for the three analyzed cases without the shift phase.

Observation: In the trajectory diagrams, the following were noted, according to [26]: A- carrier; T - Target; M – Missile.

Although all missiles hit the target, for M3 (the over-gain case), from vertical projection Figure 37 an overshoot can be seen at the entry into the second phase which caused the loss of velocity and roll velocity observed in the previous diagrams.

From the horizontal projection Figure 38, pronounced oscillations are observed for M3 (the over-gain), which can also lead to loss of velocity. Also, from both diagrams (the vertical and horizontal projection), it is observed that for M1, a greater final distance from the target is obtained but still within the imposed final distance of 3 m.

Figure 39 and Figure 40 show the vertical and horizontal projection of the trajectories for the three analyzed cases with shift phase ($\varsigma =20°$). Unlike the previously analyzed cases, if the shift phase is considered M1, with under-gain, it no longer reaches the target. The other features of the movement are preserved. For M3, an overshoot and oscillations along the trajectory are observed, which leads to a loss of velocity with a decrease in the performance of combating distant targets (over 3 km).

From the vertical projection (Figure 39) it can be seen that M1, after missing the target, continues to fly until it hits the ground.

From Figure 38 and Figure 40 a deviation of the trajectory to the left of the firing plane can be observed, a phenomenon specific to aerodynamically stabilized rolling missiles.

Figure 41 and Figure 42 show the command fill factor without and with the shift phase. It is observed that in both situations, for all missiles, the filling factor stabilizes at the equilibrium value, which achieves weight compensation. For both cases (Figure 41 and Figure 42), for M3, the filling factor saturates over a longer period of time, which leads to a decrease in the missile's maneuverability, therefore a decrease in performance, including the close targets. Also, for M3, strong oscillations are observed during the trajectory, which leads to a loss of velocity and, implicitly, a decrease in performance regarding distant targets. As for M1, it shows pronounced long-period oscillations of the filling factor, especially under shift phase conditions (Figure 42).

Figure 43 and Figure 44 show the command phase without and with the shift phase.

It is observed that in both situations, for all missiles, the phases stabilize at a small negative equilibrium value, which compensates for the Magnus and gyroscopic effects in the lateral plane.

For both cases, the phase of M3 has strong oscillations all the time, which leads to a loss of velocity, therefore a decrease in performance regarding distant targets.

Figure 45, Figure 47 and Figure 48 shows the influence of the guidance time ${(\tau }_c)$. Unlike the cases previously analyzed, in these last situations, the three gain cases will be resumed, but with a reduced value of the guidance time ${(\tau }_c=0.3)$.

Thus, in Figure 45 and Figure 46 the vertical and horizontal projection of the trajectories for the case without the shift phase, and ${\tau }_c=0.3$ is presented. Let's compare it with the corresponding diagrams in Figure 36 and Figure 37 (with tau sub c equals 0.87). We can see that the decrease to ${\tau }_c=0.3$ leads to strong oscillations of the trajectory, which denotes a tendency to destabilase the system. In cases without the shift phase, M2 hit the target, while M1 and M3 missed the target. In contrast, for the cases with the shift phase (Figure 47 and Figure 48) with ${\tau }_c=0.3$ all three missiles miss the target, which shows the importance of a correct choice of guidance parameters.

Figure 45. Target hitting – vertical view; $\varsigma =0°$ ;$\tau =0.2; {\tau }_c=0.3.$

Figure 45. Target hitting – vertical view; $\varsigma =0°$ ;$\tau =0.2; {\tau }_c=0.3.$

Figure 46. Target hitting – horizontal view; $\varsigma =0°$ ;$\tau =0.2; {\tau }_c=0.3.$

Figure 46. Target hitting – horizontal view; $\varsigma =0°$ ;$\tau =0.2; {\tau }_c=0.3.$

Figure 47. Target missing – vertical view; $\varsigma =20°$ ;$\tau =0.2; {\tau }_c=0.3.$

Figure 47. Target missing – vertical view; $\varsigma =20°$ ;$\tau =0.2; {\tau }_c=0.3.$

Figure 48. Target missing – horizontal view; $\varsigma =20°$ ; $\tau =0.2; {\tau }_c=0.3.$

Figure 48. Target missing – horizontal view; $\varsigma =20°$ ; $\tau =0.2; {\tau }_c=0.3.$

8. Organization of Results from the Appendix

Appendix A illustrates the variation of the primary aerodynamic coefficients as a function of the Mach number.

Appendix B tabulates the thrust incidence, equilibrium thrust deflection angle, and key flight quality parameters for the start and end of the flight's second phase (marching mode), based on the mechanical data in Table 1. The analysis domain, selected with reference to the velocity profiles in Figure 33 and Figure 34, covers Mach numbers from 0.34 to 0.45 within the operational flight altitude.

The characteristic polynomial coefficients were determined for the start and end of the second flight phase using the relation (82) and the flight parameters from Appendix B. From these coefficients and relations (E 5), (E 6)(E 6),(E 7), the Frank-Wall (F-W) parameters ($A_1,A_2,A_3,A_4,A_5)$ were calculated as a function of the Mach number. As shown graphically in Appendix C, the results confirm the F-W stability criterion. ${\mathrm{A}}_1>0;{\mathrm{A}}_2>0;{\mathrm{A}}_3>0;{\mathrm{A}}_4>0;{\mathrm{A}}_5>0$, is met, with the positivity of the parameters demonstrated in Figure C9, Figure C10, Figure C11 and Figure C12.

Appendix D contains the root locus of the characteristic polynomial in both cases with real coefficients and complex coefficients.

Appendix E contains the algorithm details for Frank-Wall stability criterion for 5th order polynomials.

9. Conclusions

We present a stability analysis for a slowly rolling, single-channel Semi-Automatic Command to Line-of-Sight (SACLOS) missile. The analysis begins in Section 2 with the definition of the basic motion for a single-channel slow-rolling rocket with gasodynamic command. Next, in section 3, the coupled linear form of the commanded motion of the rolling missile was obtained.

Leveraging the missile's symmetric configuration, the analysis first identifies the symmetrical and cross-coupling terms for the pitch and yaw channels, treating the latter as complex quantities. Building on this, Section 4 derives the flight quality parameters and presents a structural model for the commanded motion of the rolling missile, which is characterized by complex parameters due to the inherent channel coupling.

Subsequently, Section 5 constructs the complete guided flight model. The development begins by defining the guidance kinematics and proceeds to a detailed analysis of the single-channel actuator, including its command generation logic and error compensation methods. The formulation of the final guidance equations integrates these components into the non-linear 6-DOF framework, creating a comprehensive guided flight model.

To complete the linear analysis, the commanded motion model from Section 4 is connected with the transfer functions for the guidance device and actuator within a closed guidance loop, yielding the final closed-loop structural diagram for the SACLOS system. Within this scheme, three parameters are emphasized, specifically the guidance gain. $\left(k\right)$, the guidance time (${\tau }_c$) and the time constant ($\tau$), that will be analyzed in the following sections. Also, in Section 5, we obtained the characteristic polynomial of the transfer function associated with the closed-loop scheme. Next, based on the simplified structural scheme, a synthesis of the parameters of interest was carried out through a poly-zeros allocation method: “Standard coefficient method”.

Because in stability analysis from the following sections, we will use numerical values, in section 6, a numerical model was defined for a similar rocket in this class, for which aerodynamic, mechanical, thrust, time constants and amplification constants were provided, including the values for the guidance gain $\left(k\right)$, the guidance time (${\tau }_c$) and the time constant ($\tau$).

Section 7, dedicated to stability analysis, contains three parts. First, it is dedicated to a simplified analysis using the R-H criterion for the characteristic polynomial with real coefficients. The second is the complex analysis using the F-W stability criterion for the characteristic polynomial with complex coefficients, and the third is for a realistic 6DOF model. The results obtained were compared between the three models, resulting in significant similarity between the linear model with complex coefficients and the nonlinear model. The presence of complex-valued elements in the system's structural scheme results in a characteristic polynomial with complex coefficients. Consequently, stability was assessed using the Frank-Wall (F-W) criterion, a generalization of the Routh-Hurwitz (R-H) criterion applicable to such polynomials. For the resulting 5th-order characteristic polynomial, the F-W criterion requires five distinct parameters to be positive to ensure system stability. Due to the analytical complexity of these F-W parameters, a simplified case was first considered where the missile's roll was neglected. This simplification enabled the use of the standard R-H criterion, yielding determinants with much simpler expressions. The subsequent analysis of these five R-H parameters resulted in the definition of five corresponding system constraints. Thus, from the first condition R-H resulted the first lower limitation for the time constant $(-{\left(2\xi \Omega \right)}^{-1}<\tau )$, the second condition for regular case with damping factor subunitary ($\xi <1)$ and the time constant of nutation positive $(T>0)$ is always fulfilled. The third condition resulted in an upper limitation for the guidance gain $\left(k<k_3\right)$, and the fourth condition imposed a lower and the upper limit of the guidance gain $(0<k<k_4$) more restrictive than the third condition. Also, a fifth condition imposes a lower limit for the guidance time $\left({\tau +2\xi T<\tau }_c\right)$.

The simplified Routh-Hurwitz (R-H) analysis yielded a fifth, more restrictive constraint on the time constant (τ>0). While the R-H parameters were treated analytically, the more complex Frank-Wall (F-W) parameters required numerical evaluation (Section 7). This numerical study included evaluating the F-W stability parameters and the influence of key parameters on the system's characteristic roots. A key finding from comparing the R-H (non-rolling) and F-W (rolling) results is that roll motion contracts the stable range for the guidance gain, particularly at the upper boundary. Finally, the linear F-W stability predictions were verified through non-linear simulations of a tactical scenario. Three cases were tested: under-gained (M1), optimal-gained (M2), and an over-gained case (M3) set at the predicted stability limit.

The results of the nonlinear model confirmed the limitations obtained with the F-W stability criterion. In the first and last cases, the system performance is degraded, for the first case (the under-gain case) obtains a greater final distance from the target, with the possibility of missing the target and for the third case, having large oscillations of the guidance parameters (phase and fill factor) during the flight. Moreover, for the third case, a significant overshoot is observed in the initial phase, which, together with the oscillations along the way, leads to a loss of velocity, with the loss of the ability to combat targets at long distances. Instead, in the second case, with the guidance gain close to the operating value, the missile hit the target, the guidance parameters having an asymptotic behavior during the trajectory. In parallel with the guidance gain, the behavior of the missile in the conditions of decreasing the guidance time from ${\tau }_c=$0.87 [s] to ${\tau }_c=$0.3 [s], for the three cases that were analyzed, finding that ${\tau }_c$ decrease leads to strong oscillations of the trajectory, which denotes a tendency to destabilize the system. In cases without the shift phase, M2 hit the target, while M1 and M3 missed the target. In contrast, for the cases with a shift phase with ${\tau }_c=0.3$ all three missiles miss the target, which shows the importance of a correct choice of guidance parameters.

In summary, this paper develops a comprehensive stability analysis for a slowly rolling, single-channel SACLOS missile, starting from a non-linear 6-DOF model formulated in the Resal frame. This model is first linearized to derive a structural scheme for the commanded motion, which yields flight quality parameters with a complex-valued form due to channel coupling. Subsequently, the guidance loop components—including the kinematic relations, guidance device, and a single-channel actuator with a specific switching function—are defined and linearized. These elements are then integrated to construct a complete, closed-loop linear model of the SACLOS system. From this linear model, the system's characteristic polynomial is determined. The presence of complex coefficients necessitates the use of the Frank-Wall (F-W) stability criterion, which is compared against the simpler Routh-Hurwitz (R-H) criterion for the non-rolling case. This analysis establishes a stability domain for the guidance gain and lower limits for the guidance time constant. A key finding is that the stability domain predicted by the F-W criterion (with roll) is more restrictive than that obtained with the R-H criterion. Finally, the results of the linear stability analysis are verified through simulation with the original non-linear 6-DOF model. The non-linear system's behavior shows strong agreement with the predictions of the more accurate F-W criterion, thus validating the derived stability boundaries for the guidance gain and time constant.

The paper contains five Appendixes, as follows: Appendix A presents graphical the main terms of development of aerodynamics coefficients; Appendix B presents graphical the linear model flight quality parameters used for the characteristic polynomial of the transfer function; Appendix C presents graphical the real and imaginary part of the characteristic polynomial coefficients as well as the five F-W stability parameters; Appendix D contain the root locus of the characteristic polynomial. Appendix E detailing the calculation of F-W stability parameters.