Backstepping Sliding Mode Control for Radar Seeker Servo System Considering Guidance and Control System

This paper investigates the design of a missile seeker servo system combined with a guidance and control system. Firstly, a complete model containing a missile seeker servo system, missile guidance system, and missile control system (SGCS) was creatively proposed. Secondly, a designed high-order tracking differentiator (HTD) was used to estimate states of systems in real time, which guarantees the feasibility of the designed algorithm. To guarantee tracking precision and robustness, backstepping sliding-mode control was adopted. Aiming at the main problem of projectile motion disturbance, an adaptive radial basis function neural network (RBFNN) was proposed to compensate for disturbance. Adaptive RBFNN especially achieves online adjustment of residual error, which promotes estimation precision and eliminates the “chattering phenomenon”. The boundedness of all signals, including estimation error of high-order tracking differentiator, was especially proved via the Lyapunov stability theory, which is more rigorous. Finally, in considered scenarios, line of sight angle (LOSA)-tracking simulations were carried out to verify the tracking performance, and a Monte Carlo miss-distance simulation is presented to validate the effectiveness of the proposed method.


Introduction
A radar seeker is the "eye" of a missile, and it's one of the most important parts of and the biggest sensor in the missile. The seeker guides the missile to attack a target because of its ability of detecting and tracking. A radar seeker servo system with platform (RSSSP) is a kind of high-precision servo tracking system, which is mounted on the front of a missile to achieve stable tracking of moving targets, while the object of the controlling system is an inertially stabilized platform (ISP). An ISP is widely used in varieties of seekers, including radar seekers [1] and infrared seekers; the application of ISP also includes aerial shooting, airborne remote sensing systems [2], robotics [3], deep-space exploration, etc. The main factors that affect the tracking performance are projectile motion interference, friction torque between the shafts, and the uncertainty of modelling [4,5]. Thus, ISP is a nonlinear, time-varying system with complex disturbance and parameter perturbation.
To realize a better dynamic response performance and stronger robustness of RSSSP, researchers have tried a variety of approaches. For simple control structure, the proportional-integral-derivative control has been used in many systems [6]; however, the contradiction between robustness and rapidity can't be solved. Similarly, although H ∝ control is introduced to eliminate the sensitivity for disturbance, the robustness conflicts with its performance more sharply [7,8]. Besides that, it has strong reservation towards control performance. Furthermore, sliding mode control (SMC) is also used (1) Differently from existing studies, this paper combined the RSSSP with missile guidance and control systems to design a control algorithm, and a Monte Carlo simulation was carried out to verify the improvement of guidance precision, which is more realistic than analyzing servo systems by themselves; (2) differently from traditional research in which the reference signal is given as a specific function, this paper applies HTD to estimate system states in real time, and all signals involved were generated in real time; (3) different from traditional RBFNN, this paper proposed an adaptive RBFNN that adjusts the residual error in time, which enhances the estimation precision. No training is needed.
The paper is organized as follows: In Section 2, the model containing guidance, control, and RSSSP was built and its working preliminaries are presented. In Section 3, sliding backstepping controller was designed, along with the stability analysis. In Section 4, simulations were carried out to demonstrate effectiveness. Figure 1 shows the schematic diagram of two-axis RSSSP. We can see that stabilized platform consists of two gimbals, which are pitch gimbal and yaw gimbal, respectively. The system is driven by two servo motors; the detective sensor was placed in the inner frame. The radar seeker antenna array is a high-precision compositive sensor that detects targets by emitting radio signals and receiving signals. The length and width of the pitch gimbal were 32.8 cm and 32.8 cm, respectively. While the length and width of yaw gimbal were 28.6 cm and 28.6 cm, respectively. The diameter of the antenna array was 25.0 cm.

Constitution and Operating Principle of Two-Axis RSSSP
From Figure 1, we can see the relationships between two gimbals: gyroscopes measuring the angular rate of the pitch and yaw gimbal, angle sensors measuring the angle of the pitch and yaw gimbal, and moment sensors measuring the moment of the pitch and yaw motor. Those sensors are the foundation of control system, which offer crucial feedback state information of the RSSSP.
ISP is fixed at the projectile body to achieve target angle alignment.  Figure 1 shows the schematic diagram of two-axis RSSSP. We can see that stabilized platform consists of two gimbals, which are pitch gimbal and yaw gimbal, respectively. The system is driven by two servo motors; the detective sensor was placed in the inner frame. The radar seeker antenna array is a high-precision compositive sensor that detects targets by emitting radio signals and receiving signals. The length and width of the pitch gimbal were 32.8 cm and 32.8 cm, respectively. While the length and width of yaw gimbal were 28.6 cm and 28.6 cm, respectively. The diameter of the antenna array was 25.0 cm.

Constitution and Operating Principle of Two-Axis RSSSP
From Figure 1, we can see the relationships between two gimbals: gyroscopes measuring the angular rate of the pitch and yaw gimbal, angle sensors measuring the angle of the pitch and yaw gimbal, and moment sensors measuring the moment of the pitch and yaw motor. Those sensors are the foundation of control system, which offer crucial feedback state information of the RSSSP.
ISP is fixed at the projectile body to achieve target angle alignment.  Because of the low coupling and similar characteristic of the pitch and yaw channel [5], we chose the pitch channel to analyze. Furthermore, the guidance system was analyzed in the longitudinal plane. Figure 2 shows the angle relationship of radar seeker tracking. q is the line of sight angle (LOSA), i.e., the angle between horizontal plane and connection of missile and target. The antenna array was trying to align at the target, while, because of the relative motion between missile and target, there unavoidably existed misalignment angle q Δ , which should be eliminated. a θ is the angle between the missile lengthwise axis and horizontal plane, while g θ is the rotation angle of the antenna array surface in pitch channel.
To eliminate q Δ quickly and accurately, a designed algorithm gave demand to the servo motor to drive the antenna array. In practice, g θ plus a θ is considered as the LOSA, and the changing rate of it is the LOSA rate, which is necessary for the guidance system. Because of the low coupling and similar characteristic of the pitch and yaw channel [5], we chose the pitch channel to analyze. Furthermore, the guidance system was analyzed in the longitudinal plane. Figure 2 shows the angle relationship of radar seeker tracking. q is the line of sight angle (LOSA), i.e., the angle between horizontal plane and connection of missile and target. The antenna array was trying to align at the target, while, because of the relative motion between missile and target, there unavoidably existed misalignment angle ∆q, which should be eliminated. θ a is the angle between the missile lengthwise axis and horizontal plane, while θ g is the rotation angle of the antenna array surface in pitch channel.
To eliminate ∆q quickly and accurately, a designed algorithm gave demand to the servo motor to drive the antenna array. In practice, θ g plus θ a is considered as the LOSA, and the changing rate of it is the LOSA rate, which is necessary for the guidance system.      Figure 3 shows the pitch channel block diagram of the RSSSP, where the block within the red imaginary line stands for the servo motor, while the block within the blue imaginary line stands for the friction disturbance; θ d represents the angle conference signal of the system; k PW M is the power-amplifier coefficient;

Dynamic Model of the RSSSP
. θ is the angular rate of stabilized platform in inertial space; . ϑ is the projectile pitch angle rate; T turb is the disturbance moment; k g is the simplified transfer function of rate gyroscope; T c is the moment output of servo motor; i is the electric current of the servo motor; J L is the rotational inertia of motor load; L a is the inductance of inductance of armature winding; R a is the resistance of armature winding, C m is the moment coefficient of motor; C e is the coefficient of counter electromotive force.    Combined with the dynamic equation of the stable platform and the dynamic equation of the motor, the mathematical model of the RSSSP is acquired as follows [5]: Combined with the dynamic equation of the stable platform and the dynamic equation of the motor, the mathematical model of the RSSSP is acquired as follows [5]:

Dynamic Model of the RSSSP
where x 1 , x 2 , x 3 are state variables, which represent θ, . θ and T c , respectively u is the output of controller. The disturbance torque mainly consists of the spring torque and damping torque. The spring torque mainly results from the drag of wire when the platform is rotating towards the projectile, while the damping torque mainly results from the interaction between platform and base. The torque model is shown as Equation (2), where K w and K N are the proportionality of the spring torque and damp torque, respectively [21]: 2.3. Dynamic Model of Guidance and Control Systems Figure 4 demonstrates the relationship of missile-target motion, where V M , V T represent the velocity of missile and target, respectively, θ m , θ T are the trajectory inclination angular of missile and target, respectively, and R is the distance between missile and target. The relative motion equation between missile and target is given as follows: , , x x x are state variables, which represent , θ θ  and c T , respectively u is the output of controller. The disturbance torque mainly consists of the spring torque and damping torque. The spring torque mainly results from the drag of wire when the platform is rotating towards the projectile, while the damping torque mainly results from the interaction between platform and base. The torque model is shown as Equation (2), where w K and N K are the proportionality of the spring torque and damp torque, respectively [21]: cos cos .

Dynamic Model of Guidance and Control Systems
. .
where θ M = q − η M , θ T = q − η T , and A M , A T represent the command acceleration of missile and target, respectively [26]. The simplified missile dynamics and autopilot model are expressed as Above all, a complete model is established in Figure 5.

Remark 1.
Through the complete model in Figure 5, we can connect the servo system to the guidance system; furthermore, the motion between missile and target will be acquired (as seen in Figures 8b and 9b). Different from traditional research on RSSSP, the reference tracking signal comes from an angle signal detected by the detective sensor (seen in Figures 8a, 9a and 12) in real time. The projectile attitude motion generated by the missile control system is also regarded as disturbance . ϑ (as seen in Figures 10 and 13), which can test the disturbance-isolation ability of RSSSP. The working situation of RSSSP is close to real combat situations. Moreover, a Monte Carlo target miss-distance simulation (seen in Table 2) is accessible via the established model.

Remark 2.
The model describes the complete system relative to seeker servo systems except for signal-processing systems. Considering that the signal-processing system is independent of controller design, the system is ignored in modeling.
detective sensor (seen in Figures 8a, 9a and 12) in real time. The projectile attitude motion generated by the missile control system is also regarded as disturbance ϑ  (as seen in Figures 10 and 13), which can test the disturbance-isolation ability of RSSSP. The working situation of RSSSP is close to real combat situations. Moreover, a Monte Carlo target miss-distance simulation (seen in Table 2) is accessible via the established model.

Remark 2.
The model describes the complete system relative to seeker servo systems except for signalprocessing systems. Considering that the signal-processing system is independent of controller design, the system is ignored in modeling.

Control Problems for RSSSP
There are some troublesome characteristics in the RSSSP: (1) When there exist torque disturbances and angular-rate disturbances generated by projectile motion, high-performance angle tracking is hard to guarantee; (2) due to the change of the flight environment and the limited accuracy of mathematical modeling, coefficient uncertainty, and perturbation, tracking performance may not be guaranteed; (3) to enhance the performance of RSSSP, the states have to be estimated precisely in real time.

Control Problems for RSSSP
There are some troublesome characteristics in the RSSSP: (1) When there exist torque disturbances and angular-rate disturbances generated by projectile motion, high-performance angle tracking is hard to guarantee; (2) due to the change of the flight environment and the limited accuracy of mathematical modeling, coefficient uncertainty, and perturbation, tracking performance may not be guaranteed; (3) to enhance the performance of RSSSP, the states have to be estimated precisely in real time.

Controller Design
The control block diagram of the proposed method is shown in Figure 6.

Controller Design
The control block diagram of the proposed method is shown in Figure 6. To enhance tracking precision and robustness, a sliding mode backstepping control based on an adaptive RBFNN is proposed. The proposed method was designed under the structure of backstepping control. The system can be considered as a three-loop control system, where 1d x is the reference input of system, 2 d x can be considered as the first virtual control term for angular velocity To enhance tracking precision and robustness, a sliding mode backstepping control based on an adaptive RBFNN is proposed. The proposed method was designed under the structure of backstepping control. The system can be considered as a three-loop control system, where x 1d is the reference input of system, x 2d can be considered as the first virtual control term for angular velocity control loop, x 3d is the virtual control term for torque loop, u is the actual control term. To ensure the feasibility of the designed algorithm, newly defined states χ 2 , χ 3 , χ 4 are essential, which will be estimated by designed HTD.
. ϑ is the motion disturbance generated by the missile control system. Aiming at the problem of compound disturbance, an adaptive RBFNN was applied.D is the transferred disturbance estimated online, which would be used to offset disturbance.
Compared with traditional designs of sliding mode control, this paper applies sliding mode control only in an angular velocity loop to enhance the robustness of angular velocity tracking. To avoid a severe "chattering phenomenon", sliding mode control is not used in angular and torque tracking loops. The upper bound of residual error is also estimated online to reduce the "chattering phenomenon".
In this section, high-order tracking differentiator and adaptive RBFNN were designed at first, and were then used in controller design in Section 3.3. Sliding mode backstepping control based on adaptive RBFNN was designed for better performance of RSSSP. Furthermore, stability analysis was carried out in Section 3.4, and the semiglobal uniform ultimate stability of the system was proved.

High-Order Tracking Differentiator
In the subsequent developments, an HTD was designed to estimate newly-defined states. The HTD was formulated as follows: .
The application of HTD can be found in Figure 6, and the validation of HTD is verified in Figure 7.
From the above deduction, the estimation errors i can converge to zero by choosing an adequately large R . If we choose an infinitely large but bounded R for HTD, there exist positive constants ( 1,2,3,4) i e i = , such that ( 1,2,3,4) i i e e i ≥ = [27,28].
The application of HTD can be found in Figure 6, and the validation of HTD is verified in Figure  7.

Adaptive Neural Network
To guarantee the controller's robustness, an adaptive RBFNN is introduced to approximate the compound disturbance owing to its excellent performance and global approximation [29]. The adaptive RBFNN is defined as the mapping relationship between input vector

Adaptive Neural Network
To guarantee the controller's robustness, an adaptive RBFNN is introduced to approximate the compound disturbance owing to its excellent performance and global approximation [29]. The adaptive RBFNN is defined as the mapping relationship between input vector X = [x 1 , x 2 , . . . , x n ] T ∈ R n and the output y ∈ R [30].
where W = [w 1 , w 2 , . . . w m ] T ∈ R m donates weight vector; m and n represent the node number and input number, respectively; and h(X where c = [c 1 , c 2 , . . . , c n ] T ∈ R m and σ = σ j1 , σ j2 , . . . , σ jn T ∈ R n mean a center and a width vector of h j (X), respectively [31].
For an arbitrary continuous unknown function F(X), it has to be proven that there exists an ideal weight vector W = [w 1 , w 2 , . . . w m ] T ∈ R m , such that where ε is approximate error. It should be noted that W * and ε are unknown; their elements w 1 * , w 2 * , . . . w m * , and ε are required to be adjusted adaptively.
Define the error between the ideal weight W * and the estimated weightŴ as ε is the estimated value of ε. The adaptation laws ofŴ andε are designed in next Section. The application of adaptive RBFNN can be found in Figure 6.

Controller Design for RSSSP
Assumption 1. We assume that the LOSA q d is detected accurately.

Assumption 2.
The reference signal x 1d , its derivative .
x 1d are limited.
The algorithm design is taken as follows: Define the state error as The time derivative of z 1 is obtained by Define the virtual control law as Use χ 2 to replace x 1d , Equation (15) becomes where x 2d is available virtual control law owing to tracking differentiator, while x 3d in Equation (28) has the same meaning. Taking the derivative of x 2d , we can get ..
where k 1 is a positive constant value. Considering that T turb is a positive constant value. Considering that can't be acquired directly, the item is ignored, and will be estimated later. Substituting χ 3 , χ 4 for ..
x 2d can't be directly acquired, intermediate variables .
x 2d are defined, and will be used for controller design later.

The first Lyapunov function is chosen as
Differentiating V 1 with respect to time, and we can get To enable that the LOSA rate tracking possesses strong robustness, sliding mode control is adopted to eliminate the effect of uncertain parameters perturbance.
Define a traditional sliding mode and a complementary sliding mode as follows: where k 2 > 0 is the parameter to be designed. And the second Lyapunov function is chosen as The relationship between s and s c can be expressed as Then the virtual law is designed as Furthermore, we can get the derivative that will be used later .
Considering that T turb is an unknown function, x 3d can't be used directly; therefore, virtual control law is chosen as Besides, the derivation of x 3d is .
x 3d is acquired To guarantee the convergence of z 3 , the third Lyapunov is chosen as Taking the derivation of z 3 , we get the following equation: .
Design the actual control law as follows where the whole disturbance D is presented as Considering that D is a complex function with z 1 , z 2 , z 3 as arguments, an adaptive RBFNN is introduced to approximate D on line.  ϑ is seduced by guidance command, which can be written as . ϑ(a n ), while guidance command is determined by the LOSA rate, which can be seen as a n (z 1 , z 2 , z 3 ). Above all, . ϑ can be seen as the function with z 1 , z 2 , z 3 as arguments, which can be rewritten as . ϑ(a n (z 1 , z 2 , z 3 )). Therefore, . ϑ can be estimated by RBFNN.
Therefore, D can be expressed as The control law is rewritten as To develop the adaption laws of W andε, define Lyapunov function Combining with Equations (20), (28), (30) and (31), the derivation of V 4 is acquired. .
Since W * is a constant, it should be noted that W.
According to Equation ( (41), it can be observed that the global boundedness of disturbance is guaranteed. Besides, the chattering phenomenon caused by sliding mode control is eliminated.
and define the following compact sets: It can be seen that . V 4 will be negative if z 1 / ∈ Ω z 1 or z 2 / ∈ Ω z 2 or z 3 / ∈ Ω z 3 . Hence, z 1 , z 2 , z 3 are semiglobally uniformly ultimately bounded. This is the end of proof.

Simulation and Analysis
Simulations were carried out in this section to verify the effectiveness of the proposed control scheme. We consider two different combat scenarios to verify the algorithm. The designed algorithm was applied to the missile RSSSP. In the first scenario, the missile attacked the target from upwards. In the second scenario, the missile attacked the target from downwards and the target was carrying out sine maneuvering.
The parameters of the RSSSP are shown in Table 1. Case 1. In this case, we assume that a missile intercepts a low-altitude target. The positions of missile and target are (0 m, 2000 m) and (5000 m, 100 m), respectively; their speeds are 1000 m/s and (400 + 30 × t)m/s, respectively. Thus, the LOSA is 0.36 rad, and the initial servo system rotation angle was 0.36 rad, while the initial trajectory inclination angles were 30 • and 0 • , respectively. Because of the multipath effect of a radar seeker while interpreting low-altitude targets, Brewster angle restraint has to be guaranteed, and the double sliding mode guidance law [32] is designed as where λ, k, ε are positive design parameters, S 1 = x 1 − x 1d , S 2 = x 2 , S = S 2 + λV M S 1 /R. The coefficients of the missile control system are T α = 1.5, T Auto = 0.06, respectively. The order of HTD is chosen as 2, and the design parameters of HTD are chosen as follows: R 1 = 200, R 2 = 20, R 3 = 5, a 11 = 10, a 21 = 10, a 31 = 10, a 12 = 10, a 22 = 10, a 31 = 10; the control design parameters are k 1 = 10, k 2 = 3, k 3 = 0.0125.
Parameters of RBFNN are chosen as γ 1 = 1000, γ 2 = 3, the number of node is 100, and the center of RBFNN is evenly spaced in c ∈ [−50, 49], the width is chosen as b = 15. In particular, the mentioned parameters are all designed, whileŴ andε are adaptive vector and adaptive parameter that adjust online.
The effectiveness of HTD is shown as follows; it can be observed from Figure 16 that the estimation of LOSA is perfect.
To show the priority, the proposed method is compared with a prescribed performance controller in Reference [20].
From Figures 8a and 9a, we can see that the convergence speed and tracking precision of the proposed method are better than the method in Reference [33]. Figure 10 shows the projectile disturbance generated in real time, and the figure shows that the disturbance is sharp. While Figure 11 shows the disturbance estimated by adaptive RBFNN. The designed adaptive RBFNN achieves precise compensation of disturbance, which will enhance the tracking ability. From Figures 8b and 9b, we can see that a missile with the two mentioned methods hits the target in high precision. We can conclude that in this scenario, although the tracking performance of method in Reference [20] is worse than proposed method, its guidance precision is satisfying. The controller output contrast in Figure 12 shows that the controller output is smooth except for the initial big overshoot, while a sharp "chattering phenomenon" occurs in the sliding mode controller. We can conclude that the proposed method eliminates chattering of the controller output.

Remark 6.
It should be noted that Figure 10 shows the disturbance generated by projectile motion . ϑ, while Figure 11 shows total disturbance D = L a C m k PWM ( z 2 z 3 (s + s c + 2k 1 ) + (T turb + (2 − 2k 2 J L ) .
T turb + C e C m . ϑ L a )/z 3 ). The difference should be clarified.           Guidance law is chosen as the proportional navigation law, and the proportionality coefficient Case 2. In this case, we assume that an air-defense missile hits a plane target. The positions of missile and target are (0 m, 10, 000 m) and (8000 m, 13, 000 m), respectively; their initial speeds are 1000 m/s and 600 m/s, respectively, and the speed of target in y-axis direction is V ty = 100 × cos(0.2πt)(m/s). Thus, the LOSA is −0.358 rad, the initial servo system rotation angle is −0.34 rad, while the initial trajectory inclination angles are −30 • and 0 • , respectively.
Guidance law is chosen as the proportional navigation law, and the proportionality coefficient is K.

Remark 7.
The rules of control parameters are illustrated as follows. For a closed-loop control system, the primary purpose is to guarantee the stability based on the Lyapunov theorem of stability, while the stability interval of system falls in compact sets (Equation (46)). By choosing appropriate control parameters k 1 , k 2 , k 3 within the interval (Equation (45)), the system will be stable. Parameter R determines the response speed and estimation precision of HTD, its positive correlation, while a i determines the convergence characteristic of the ith term, its positive correlation, too. Parameters γ 1 , γ 2 , b of RBFNN determines the convergence speed, besides, γ 1 , γ 2 determines the estimation precision directly. The center c is adjusted by the specific range of state variety, while the number of nodes depends on the complexity of estimation. The parameters of the control system have to be adjusted repeatedly.
It can be seen from Figure 13 that the proposed method converges quicker; besides, when there exists a 30% model parameters perturbation, the proposed method precisely tracks the reference, while the precision and robustness of the method in Reference [20] are worse. Figure 14 shows the projectile generated in real-time, while Figure 15 shows total disturbance estimated by RBFNN. It can be observed from Figure 16 that a higher tracking performance of servo system advances hit time, and the trajectory is gentler.
interval of system falls in compact sets (Equation (46)). By choosing appropriate control parameters 1 2 3 within the interval (Equation (45)), the system will be stable. Parameter R determines the response speed and estimation precision of HTD, its positive correlation, while i a determines the convergence characteristic of the ith term, its positive correlation, too. Parameters 1 2 , ,b γ γ of RBFNN determines the convergence speed, besides, 1 2 , γ γ determines the estimation precision directly. The center c is adjusted by the specific range of state variety, while the number of nodes depends on the complexity of estimation. The parameters of the control system have to be adjusted repeatedly.
It can be seen from Figure 13 that the proposed method converges quicker; besides, when there exists a 30% model parameters perturbation, the proposed method precisely tracks the reference, while the precision and robustness of the method in Reference [20] are worse. Figure 14 shows the projectile generated in real-time, while Figure 15 shows total disturbance estimated by RBFNN. It can be observed from Figure 16 that a higher tracking performance of servo system advances hit time, and the trajectory is gentler.
(a) LOSA tracking of method in [20] (b) LOSA tracking of proposed method  interval of system falls in compact sets (Equation (46)). By choosing appropriate control parameters 1 2 3 , , k k k within the interval (Equation (45)), the system will be stable. Parameter R determines the response speed and estimation precision of HTD, its positive correlation, while i a determines the convergence characteristic of the ith term, its positive correlation, too. Parameters 1 2 , ,b γ γ of RBFNN determines the convergence speed, besides, 1 2 , γ γ determines the estimation precision directly. The center c is adjusted by the specific range of state variety, while the number of nodes depends on the complexity of estimation. The parameters of the control system have to be adjusted repeatedly.
It can be seen from Figure 13 that the proposed method converges quicker; besides, when there exists a 30% model parameters perturbation, the proposed method precisely tracks the reference, while the precision and robustness of the method in Reference [20] are worse. Figure 14 shows the projectile generated in real-time, while Figure 15 shows total disturbance estimated by RBFNN. It can be observed from Figure 16 that a higher tracking performance of servo system advances hit time, and the trajectory is gentler.
(a) LOSA tracking of method in [20] (b) LOSA tracking of proposed method      To verify and demonstrate the advanced designed servo algorithm, a Monte Carlo simulation was carried out [34,35]. Through the Monte Carlo simulation, the effect of the improved servo algorithm is clearly shown through target-miss distance.
It can be seen from Table 2 that the designed method eliminates the miss distance in the considered scenarios.

Conclusions
In this paper, a backstepping sliding mode control based on an adaptive RBFNN is proposed to keep precise and stable tracking of seeker servo systems. A complete model of SGCS was built, and the algorithm design was carried out based on HTD. It is proved that all closed-loop system signals are semiglobally uniformly ultimately bounded. The designed method enhances the tracking robustness and precision in large scale, and solves problems of uncertain disturbance and parameters perturbation. A method in Reference [20] was compared with the proposed method, and the better performance of proposed method was proved through LOSA tracking contrasts. Besides that, Monte Carlo simulations showed that the proposed method enhances guidance precision.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, and in the decision to publish the results.

Data Availability Statement:
The data in this paper include two parts: one is model-parameter data and the other is simulation results generated by digital simulation. The model data are mostly derived from references, which have been remarked. The simulation results are derived from our original simulation program. All necessary data are presented in this paper.