Next Article in Journal
Performance Analysis of Global Navigation Satellite System Signal Acquisition Aided by Different Grade Inertial Navigation System under Highly Dynamic Conditions
Previous Article in Journal
Chemiresistor Devices for Chemical Warfare Agent Detection Based on Polymer Wrapped Single-Walled Carbon Nanotubes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

An Adaptive B-Spline Neural Network and Its Application in Terminal Sliding Mode Control for a Mobile Satcom Antenna Inertially Stabilized Platform

1
School of Instrumentation Science and Opto-Electronics Engineering, Beihang University, Beijing 100191, China
2
Beijing Institute of Control & Electronic Technology, Beijing 100038, China
*
Author to whom correspondence should be addressed.
Sensors 2017, 17(5), 978; https://doi.org/10.3390/s17050978
Submission received: 10 January 2017 / Revised: 5 April 2017 / Accepted: 25 April 2017 / Published: 28 April 2017
(This article belongs to the Section Physical Sensors)

Abstract

:
The mobile satcom antenna (MSA) enables a moving vehicle to communicate with a geostationary Earth orbit satellite. To realize continuous communication, the MSA should be aligned with the satellite in both sight and polarization all the time. Because of coupling effects, unknown disturbances, sensor noises and unmodeled dynamics existing in the system, the control system should have a strong adaptability. The significant features of terminal sliding mode control method are robustness and finite time convergence, but the robustness is related to the large switching control gain which is determined by uncertain issues and can lead to chattering phenomena. Neural networks can reduce the chattering and approximate nonlinear issues. In this work, a novel B-spline curve-based B-spline neural network (BSNN) is developed. The improved BSNN has the capability of shape changing and self-adaption. In addition, the output of the proposed BSNN is applied to approximate the nonlinear function in the system. The results of simulations and experiments are also compared with those of PID method, non-singularity fast terminal sliding mode (NFTSM) control and radial basis function (RBF) neural network-based NFTSM. It is shown that the proposed method has the best performance, with reliable control precision.

1. Introduction

A mobile satcom antenna (MSA) is a kind of satellite communication antenna which can maintain continuous and reliable communications between satellites and users in movement. MSAs are widely applied in vehicle carriers, flying platforms and vessels. The Jet Propulsion Laboratory (JPL) has developed several mobile vehicular antenna systems for the satellite-based applications [1]. Due to the uncertain working environment, such as bumpy roads, sea waves and airflows, there might be strong disturbances in actual systems [2]. MSAs seeking to achieve high data-rate communications require the inertial antenna to be capable of pointing to within fractions of a degree [3]. Stabilization is used in the antenna pointing control system to maintain the line of sight by isolating the disturbances caused by vehicle motion [4]. Generally speaking, the attitude sensor can be an attitude and heading reference system (AHRS) or an inertial measurement unit (IMU) [5]. The outputs of these inertial sensors are widely used as feedbacks for stabilization in antenna control systems [6]. The accuracy of control is influenced by many factors, mainly including the existence of unmodeled nonlinear uncertainty, random bias of states and sensor noises [7] and friction restriction, which decreases the accuracy most. The disturbances and uncertainties can lead to the degradation of tracking and stabilization accuracy of inertial stabilized platform in engineering applications [8].
Various decoupling control and disturbance rejecting methods have been developed to enable MSA with high performance. Though the PID control method is widely applied in engineering practice due to its simplicity and reliability, it is difficult to reduce the disturbances and uncertainties rapidly [8]. Therefore, this method offers low control accuracy and robustness. A feedforward scheme for a two axis inertial stabilized platform (ISP) was adopted to reject the periodic disturbing torque acting on the payload due to the static mass unbalance [9,10], but the accelerometer signal may be easily distorted by the centripetal acceleration [11]. The H control method presents strong robustness, but it is highly conservative and the control accuracy is usually sacrificed [12]. Internal mode control can improve the disturbance rejection capability and robustness against model uncertainties, but the exact structural property of the disturbance needs to be known in advance [13]. The sliding mode control method is studied in depth for its unique robust control performance for nonlinear and uncertain systems. In the linear sliding mode, the deviation between the system and the desired trajectory is asymptotically convergent in exponential forms [14]. TSM control is improved and compared with the linear hyperplane based sliding modes, showing that the TSM offers some superior properties such as fast response, finite time convergence [15,16].
However, the conventional TSM has many defects such as singular problem and chattering phenomena. An initial discussion about avoiding the singularity in TSM control was presented in [17]. In order to avoid the singular problem, many non-singular TSMs were proposed [14,15,18,19,20]. The good robustness of TSM is based on the large switching control gain, which is determined by the uncertain disturbances [20]. If the switching control gain is small, the uncertainty issues cannot be eliminated, whereas if the gain is large enough, the chattering phenomena would be enhanced. In many practical control systems, including DC motors and aircraft control system, it is important to avoid control chattering by providing continuous or smooth control signals [21]. According to the high order sliding mode algorithm proposed in [22], “chattering removal” can be achieved by combining the arbitrary-order sliding mode controller with the dynamic sliding mode. Another well-known way to reduce chattering problem in TSM is to use tanh function instead of the sign function [23]. Neural networks can learn and approximate any arbitrary nonlinear functions [24]. Therefore, neural networks are implemented to approximate uncertain disturbances in systems as well as to reduce the chattering [25].
A radial basis function (RBF) neural network has several important features such as simple structure, fast learning and better approximation capabilities. It has been adopted in TSM control systems [24,26]. Wavelet neural networks, which are applied to terminal sliding mode control [27,28], are a new class of neural networks that have been developed using a combined method of multi-layer artificial neural networks and wavelet analysis. The B-spline function is a piecewise polynomial function widely used in computer-aided design and computer graphics. For its excellent local features, the B-spline function can be used as an activation function of neural network. The B-spline neural networks (BSNNs) can be divided into two groups: one group is based on the B-spline basis function [29,30,31,32] and the other group is based on the B-spline curve [33,34,35,36,37,38].
The B-spline basis function is only related to the internal knots, hence the definite internal knots can derive the definite B-spline basis curve. For those kinds of BSNNs, only the weights are trained in some studies [29,31]. The internal knots are trained by a heuristic algorithm in some other studies [30,32], making the shapes of B-spline basis function changeable. In this group, the input of the neural network can be considered as the value of internal knots, then according to the b-Spline basis function, the output of the neural network can be obtained.
In the group where the B-spline curve is adopted, it is known that the B-spline curve is a linear combination of control points and B-spline basis functions. If the positions of the control points are changed, the shapes of B-spline curves are also changed. In some studies, the weights are considered as control points [38]. Through the training of weights, the BSNN can then become adaptive. In other references, control points are trained by rules. In [33,37], the position of control points can be changed to a known nearby value, but only one control point can be changed once by this method. Additionally, the input and output of neural network are not relevant directly in this group.
In this paper, a novel B-spline neural network with the capability of shape adjustment is proposed and applied to approximate the external disturbance and the unmolded dynamics of the system. The equation of the B-spline curve is reconstructed in numeric form, therefore the new B-spline curve possesses the properties of the radial basis function, making it different from the previous approaches. In this research, the activation function of proposed BSNN is a pp-form spline.
The paper is organized as follows: in Section 2, the dynamic model of MSA and the control block diagram are established respectively. Then a new terminal sliding mode control strategy is proposed and proved. In Section 3, a novel BSNN is proposed and the parameter updating rules are built. In Section 4, experiments and simulations are carried out. The performances of the proposed methods are also compared with other methods. Finally, conclusions and contributions are summarized in Section 5.

2. Modeling of a Mobile Satcom Antenna

The paper focuses on a double gimbal MSA. The innermost gimbal is an elevation gimbal on which the antenna array is mounted, while the outmost gimbal is an azimuth gimbal which is mounted on the base. According to the coordinates system definition and the conversion relationships between different coordinates, a dynamic model of the mobile satcom antenna is derived.

2.1. Gimbal Coordinates Definition

An orthogonal coordinate system is defined: Base coordinate frame ( o x b y b z b ) , azimuth gimbal frame ( o x a y a z a ) and elevation gimbal frame ( o x e y e z e ) . The rotating transformations are shown in Figure 1. As we can see from the Figure 1, θ a and θ ˙ a are the relative angular and relative angular rate between the azimuth gimbal and the base carrier, respectively. θ e and θ ˙ e are the relative angular and relative angular rate between the elevation gimbal and the azimuth gimbal, respectively:
Rotation matrices can be written as C b a = [ cos θ a sin θ a 0 sin θ a cos θ a 0 0 1 1 ] and C a e = [ 1 0 0 0 cos θ e sin θ e 0 sin θ e cos θ e ] , which represent the transformation from the base coordinate to the azimuth coordinate and transformation from azimuth coordinate to the elevation coordinate, respectively.
According to the rotation matrix C b a and C a e , the angular rate relations of the two gimbals can be obtained as follows:
ω i a a = C b a ω i b b + ω b a a = [ p cos θ a + q sin θ a p sin θ a + q cos θ a r + θ ˙ a ]
ω i e e = C a e ω i a a + ω a e e = [ p cos θ a + q sin θ a + θ ˙ e p cos θ e sin θ a + q cos θ a cos θ e + r sin θ e + θ ˙ a sin θ e p sin θ a sin θ e q cos θ a sin θ e + r cos θ e + θ ˙ a cos θ e ]
in which, ω i b b is the angular rate of the base carrier; ω i a a is the angular rate of the azimuth gimbal; ω i e e is the angular rate of the elevation gimbal; ω b a a = [ 0 0 θ ˙ a ] T represents the angular rate of the azimuth gimbal with respect to the base; ω a e e = [ θ ˙ e 0 0 ] T represents the angular rate of the elevation gimbal with respect to the azimuth gimbal; ω i b b = [ p q r ] T is the angular rate of the mobile carrier.

2.2. Dynamic of the Azimuth Gimbal

The gimbals in this system are all rigid, therefore the basic Newton-Euler rotation equation can be written as follows:
M = H ˙ + ω × H H = J ω
in which M is the resultant moment of the force added to the rigid body; H is the inertial angular momentum; J is the inertial moment; ω is the absolute angular rate.
Since the structures of the two gimbals are similar (the difference is that the azimuth gimbal is effected by coupling of elevation gimbal), only the azimuth gimbal is considered in this paper.
According to the Equation (3), the dynamic model of azimuth gimbal can be obtained as follows:
( J a z + J e z ) ( r ˙ + θ ¨ a ) = M a m M a d M a e z ( J a y + J e y J a x J e x ) ω i a x a ω i a y a
where, M a d is the disturbance torque added to elevation gimbal; M a m is the driving torque of elevation motor; M a e z is the coupling torque that effected by elevation gimbal. They are given by Equations (5) and (6):
M a e z = J e z ω ˙ i e z e sin θ e + ( J e x J e y ) ω i e x e ω i e y e sin θ e + J e z ω ˙ i e z e cos θ e + ( J e x J e y ) ω i e x e ω i e y e cos θ e
M a m = k a t R a a ( N a u a a N a 2 k a b θ ˙ a ) ( N a 2 J a d + J a f ) θ ¨ a ( N a 2 f a d + f a f ) θ ˙ a
J a x , J a y , J a z are the moments of inertia of the azimuth gimbal related to the X , Y , Z axes in azimuth coordinates, respectively; J e x , J e y , J e z are the moments of inertia of the elevation gimbal related to the X , Y , Z axes in the elevation gimbal, respectively; k a t is the azimuth motor torque constant; k a b is the back-EMF coefficient of azimuth motor; N a is the azimuth gear ratio; R a a is the motor resistance; J a d is the moment of inertial of driver gear; f a d is the viscous friction coefficient of driver gear; J a f is the moment of inertial of passive gear; f a f is the viscous friction coefficient of passive gear; u a a is the input driven voltage of the motor.
Based on the preceding analysis, system (4) can be simplified into a standard second order control system:
{ x ˙ a 1 ( t ) = x a 2 ( t ) x ˙ a 2 ( t ) = A a x a 2 ( t ) + B a u a ( t ) + C a d a ( t ) + D a w a ( t )
in which, x a 1 = θ a , x a 2 = θ ˙ a , d a ( t ) is external disturbance and w a ( t ) is the coupling effect torque. To further simplify the analysis, we define that | C a d a ( t ) + D a w a ( t ) | l d . Then the proposed control block diagram can be obtained as follows:
A a = [ 0 1 0 k a t k a b N a 2 + R a a N a 2 f a d + R a a f a f ( J a z + J e z + N a 2 J a d + J a f ) R a a ] , B a = [ 0 k a t N a ( J a z + J e z + N a 2 J a d + J a f ) R a a ]
C a = [ 0 1 J a z + J e z + N a 2 J a d + J a f ] , D a = [ 0 1 J a z + J e z + N a 2 J a d + J a f ]
and d a ( t ) = M a d + ( J a z + J e z ) r ˙ + ( J a y + J e y J a x J e x ) ω i a x a ω i a y a .
Then the control block diagram shown in Figure 2 can be obtained.
The NFTSM control applied in this paper can guarantee that the system state arrives at the equilibrium point in a finite time with fast response and higher precision. In the following section, a new NFTSM controller is proposed and proved.

2.3. Non-Singularity Fast Terminal Sliding Mode

The nonlinear system is defined by Equation (7). In order to improve the control performance of terminal sliding mode, the paper proposed an improved NFTSM as shown in Equation (8), which is similar to reference [14].
s = e ^ + 1 α e ^ r + 1 s i g n ( e ^ ) + 1 β e ^ ˙ p + 1 s i g n ( e ^ ˙ ) = 0
in which, α , β > 0 , 0 < p < 1 and r > p , e ^ = θ a θ a c . The first derivative can be expressed as
s ˙ = e ^ ˙ + r + 1 α e ^ r e ^ ˙ + p + 1 β e ^ ˙ p e ^ ¨
Theorem 1.
For System (7) with the adopted NFTSM, if the control law is designed as:
u n = B a 1 [ β p + 1 e ^ ˙ 1 p s i g n ( e ^ ˙ ) ( 1 + r + 1 α e ^ r ) + A a e ^ ˙ + ( l d + δ ) s i g n ( s ) + ( A a θ ˙ a c θ ¨ a c ) ]
the NFTSM manifold will be reached in a finite time. In addition, the tracking error on the sliding mode surface will also converge to zero in a finite time.
Proof. 
Consider the Lyapunov candidate function as:
V = 1 2 s 2
The derivative of V is given by the following function,
V ˙ = s s ˙ = s ( e ^ ˙ + r + 1 α e ^ r e ^ ˙ + p + 1 β e ^ ˙ p A a e ^ ˙ + p + 1 β e ^ ˙ p B a u n ( t ) + p + 1 β e ^ ˙ p C a d n ( t ) + p + 1 β e ^ ˙ p D a w n ( t ) + p + 1 β e ^ ˙ p ( A a θ ˙ a c θ ¨ a c ) ) = s ( p + 1 β e ^ ˙ p ( ( l d + δ ) s i g n ( s ) ) + p + 1 β e ^ ˙ p C a d n ( t ) + p + 1 β e ^ ˙ p D a w n ( t ) ) s ( δ s i g n ( s ) ) = | s | δ
It is obvious that V ˙ < 0 if σ ˙ = 0 . Thus, the condition for Lyapunov stability is satisfied. Then the system state will converge to zero within finite time. As shown in Equation (10), the value of 1 p < 0 is always greater than zero, therefore, the controller is non-singular. ☐
We can also obtain the reaching time t r satisfies (see Appendix A)
t r p + 1 p β 1 p + 1 | e ( t 0 ) | p p + 1
That completes the proof. Therefore, the system can reach the terminal sliding mode surface within a finite time.
For terminal sliding mode control methods, the upper disturbance bound is generally required, but in practical applications, a large upper bound can be used to eliminate the external disturbance item and guarantee the robustness. However, the chattering phenomena is enhanced correspondingly. The paper attempts to approximate nonlinear disturbance by using BSNN. A novel BSNN is proposed in this paper. As its name implies, the B-spline curve is adopted as activation function.

3. B-Spline Neural Network

Artificial neural networks have high nonlinear approximation abilities. A spline function is a piecewise polynomial function of degree k . The spline neural network can play a role in local features. In this section, a new kind of B-spline function is introduced herin, then the proposed BSNN is deduced.

3.1. B-Spline Basis Function Definition

A B-spline curve differs from a Hermite or Bérzier curve, because a B-spline curve usually consists of more than one curve segment. The B-spline curve is widely used in the field of computer graphics for its excellent performance. The B-spline basis function is the foundation of BSNN.
Definition 1
[39]. Given a knot vector T = { t 0 t 1 t m } , the i th B-spline basis function of degree k can be written as:
B i , 0 ( t ) = { 1 ,   i f   t i < = t < t i + 1 0 , o t h e r w i s e
B i , k ( t ) = t t i t i + k 1 t i B i , k 1 ( t ) + t i + k t t i + k t i + 1 B i + 1 , k 1 ( t )
Remark 1.
If t i + 1 t i ( 0 i m 1 ) is a constant, the B-spline curve is called uniform; otherwise, it is called non-uniform.
Definition 1 is usually referred to the Cox-de Boor recursion formula.
The j th B-spline basis curve with degree k can be written as B j , k i ( t ) on an arbitrary vector [ t i , t i + 1 ) :
B j , k i ( t ) = a j , 0 i + a j , 1 i t + a j , 2 i t 2 + + a j , k i t k = [ t k t k 1 1 ] [ a j , k i a j , 1 i a j , 0 i ]
An order k + 1 B-spline is formed by joining several pieces of polynomials of degree k with at most C k 1 continuity at the breakpoints [40], hence the following formula can be obtained:
[ t i k t i k 1 t i 1 t i + 1 k t i + 1 k 1 t i 1 k t i + 1 k 1 ( k 1 ) t i + 1 k 2 1 0 k ! t i + 1 ( k 1 ) ! 0 0 ] [ a j , k i a j , k 1 i a j , 1 i a j , 0 i ] = [ B j , k i ( t i ) B j , k i ( t i + 1 ) B j , k i ( 1 ) ( t i + 1 ) B j , k i ( k 1 ) ( t i + 1 ) ]
The derivative of Equation (15) can be written as follows:
d d u B i , k ( t ) = B i , k ( t ) = k t i + k t i B i , k 1 ( t ) k t i + k + 1 t i + 1 B i + 1 , k 1 ( t )
An arbitrary order derivative formula can be obtained as Equation (19) shows:
B i , k ( d ) ( t ) = k ( B i , k 1 ( d 1 ) ( t ) t i + k t i B i + 1 , k 1 ( d 1 ) ( t ) t i + k + 1 t i + 1 )
According to Equations (14), (15), (18) and (19), the coefficients a j , k i , a j , k 1 i , , a j , 1 i , a j , 0 i in Equation (16) can be obtained. Then the j th B-Spline basis curve B j , k i ( t ) with degree k on vector [ t i , t i + 1 ) can be calculated.
The B-Spline basis function B j , k ( t ) with degree k can also be represented in matrix form in the interval [ t 0 , t m ] ,
B j , k ( t ) = [ B j , k 0 ( t ) B j , k 1 ( t ) B j , k m 1 ( t ) ] = [ t k , t k 1 , , 1 ] [ a j , k 0 a j , k 1 a j , k m 1 a j , 1 0 a j , 1 1 a j , 1 m 1 a j , 0 0 a j , 0 1 a j , 0 m 1 ] = v t M e
The expression of different B-spline basis functions in an interval [ t i , t i + 1 ) can be written as follows:
B k i ( t ) = [ B 0 , k i ( t ) B 1 , k i ( t ) B n , k i ( t ) ] = [ t k   t k 1     1 ] [ a 0 , k i a 1 , k i a n , k i a 0 , 1 i a 1 , 1 i a n , 1 i a 0 , 0 i a 1 , 0 1 a n , 0 i ] = v t M e
The coefficient matrixes M e in Equation (20) and M e in Equation (21) can be obtained by Equations (14), (15), (18) and (19).
The B-spline curve is a linear combination of control points P i and a B-spline basis function B i , k ( t ) . The definition is as follows:

3.2. B-Spline Curve Definitions

Definition 2
[39]. Given n + 1 control points P i , a knot vector T = { t 0 t 1 t m } and B-spline basis functions B i , k ( t ) , a B-spline curve is given by:
C ( t ) = i = 0 n B i , k ( t ) P i
Remark 2.
Indexes m , n , k must satisfy m = n + k + 1 .
Remark 3.
If the first knot and the last knot have multiplicity with value k + 1, the B-spline curve is called closed clamped; otherwise, it is called an open clamped B-spline curve. More specific illustrations can be seen in Figure 3.
From Equation (22) on the i th vector the B-spline function can be written as:
C i ( t ) = [ B 0 , k i   B 1 , k i     B n , k i ] [ P 0   P 1     P n ] T = v t M e [ P 0   P 1     P n ] T = v t M e [ ( P x 0 , P y 0 )   ( P x 1 , P y 1 )     ( P x n , P y n ) ] T
If the horizontal value x i n of an arbitrary point ( x i n , y i n ) on the B-spline curve is known, the following formula can be obtained by:
x i n = v t M e [ P x 0   P x 1     P x n ] T
Therefore, the unknown internal vector v t ^ can be calculated. Furthermore, the vertical value y i n can be obtained:
y i n = v t ^ M e [ P y 0   P y 1     P y n ] T
In this paper, order 3 B-spline curve is adopted. Then Equation (23) can be simplified to the following form:
C j ( u ) = [ t 2   t 1   1 ] [ a 2 1 a 2 2 a 2 3 a 1 1 a 1 2 a 1 3 a 0 1 a 0 2 a 0 3 ] [ P j   P j + 1   P j + 2 ] T = [ t 2   t 1   1 ] M e [ P j   P j + 1   P j + 2 ] T
x i n can be represented as x i n = [ t 2 t 1 ] A β , in which M e = ( α 1 , α 2 , α 3 ) T , β = [ P x 0 P x 1 P x 2 ] T . According to the known horizontal coordinate value x i n , the corresponding internal knot t ^ can be obtained:
t ^ = ( α 2 · β ) ± ( α 2 · β ) 2 4 ( α 1 · β ) [ ( α 3 · β ) x i n ] 2 ( α 1 · β ) , ( α 1 · β ) 0
Furthermore:
y i n = [ t ^ 2 t ^ 1 ] [ ( α 1 γ ) ( α 2 γ ) ( α 3 γ ) ] T = ( α 1 γ ) t ^ 2 + ( α 2 γ ) t ^ + ( α 3 γ )
in which γ = [ P y 0 P y 1 P y 2 ] T .
According to the Definition 2, the B-spline function is improved in this paper by proposing two deformation factors called translation factor and scaling factor.
Definition 3.
Given n + 1 control points P i and a knot vector T = { t 0 t 1 t m } , a new B-spline function of degree k can be written as:
C ( t ) = i = 0 n P i B i , k ( t ) = i = 0 n ( P x i + λ t i + κ w i , P y i ) B i , k ( t ) = i = 0 n [ P i + ( λ b i , 0 ) + ( κ w i , 0 ) B i , k ( t ) ]
Remark 4.
At the beginning, the value of the two deformation factors λ and κ are zero. thus, the formulas in Definition 2 and 3 are equivalent. λ is called translation factor, and κ is called scaling factor.
Assumption 1.
The two factors specify how the shape changes, assuming that b = [ 1 1 1 ] 1 × n , w = [ 1 1 1 ρ 0 τ 1 1 1 ρ ] 1 × n and 2 × ρ + τ = n ( τ can be 0 or 1).

3.3. The B-Spline Neural Network

Compared with some previous studies, the proposed B-spline adopts an order 3 B-spline function as activation function rather than a B-spline basis function as in some references [29,31]. Generally speaking, neural networks always have more than one hidden layer. The structure of our BSNN is shown in Figure 4.
As shown in Figure 4, the neural network system has multiple inputs and one output. Each input corresponds to several activation functions. The output of BSNN can be mapped from inputs x i to an output L d by using B-splines as activation functions. The output of the hidden layer can be written as follows:
h l = i = 1 a η i , j ( x i ) , j = 1 , 2 , , ρ
In which, η i , j ( x i ) means the i th input corresponding to the j th activation function; ρ means the number of inputs; a means the number of activation functions. Therefore, the output of BSNN is as follows:
L d = l = 1 L ω l h l = l = 1 L ω l ( i = 1 ρ η i , j ( x i ) ) , L = a ρ
The training method is to minimize the error function which is defined as follows:
E ( t ) = V
The chain rules of parameter update can be written as follows:
E s s l d = p + 1 β e ˙ p s
Δ ω i = E ω i = E s s l d l d ω i = p + 1 β e ˙ p s h l
ω i ( t ) = ω i ( t 1 ) + η Δ ω i + α [ ω i ( t 1 ) ω i ( t 2 ) ]
Δ λ i , j = E λ i , j = E s s l d l d h l h l η i , j η i , j λ i , j
λ i , j ( t ) = λ i , j ( t 1 ) + η Δ λ i , j + α [ λ i , j ( t 1 ) λ i , j ( t 2 ) ]
Δ κ i , j = E κ i , j = E s s l d l d h l h l η i , j η i , j κ i , j
κ i , j ( t ) = κ i , j ( t 1 ) + η Δ κ i , j + α [ κ i , j ( t 1 ) κ i , j ( t 2 ) ]
l d h l h l η i , j = w l s = 1 s i ρ 1 η s , j ( x s ) , j = 1 , 2 , a
In these equations, η is the learning rate; α is the inertial coefficient.
η i , j ( x j ) λ i , j = [ s = 0 n ( P s i , j + λ i , j t s + κ i , j b s ) B s , k ( t ^ ) ] λ i , j = [ s = 0 n P s i , j B s , k ( t ^ ) + s = 1 n λ i , j t s B s , k ( t ^ ) + s = 1 n κ i , j b s B s , k ( t ^ ) ] λ i , j = s = 0 n t s B s , k ( t ^ )
The B-spline basis function has an important property which is called “Partition of Unity”, indicating that the sum of all non-zero degree k basis functions on span [ t i , t i + 1 ) is 1, i.e., s = 0 n B s , k ( t ) = 1 .
According to Assumption 1, the following Equations (42) and (44) can be obtained:
η i , j ( x j ) λ i , j = s = 0 n B s , k ( t ^ ) = 1
η i , j ( x j ) κ i , j = [ s = 0 n ( P s i , j + λ i , j t s + κ i , j w s ) B s , k ( t ^ ) ] κ i , j = [ s = 0 n P i , j B s , k ( t ^ ) + λ i , j t s B s , k ( t ^ ) + κ i , j w s B s , k ( t ^ ) ] κ i , j = s = 0 n w s B s , k ( t ^ )
η i , j ( x j ) κ i , j = { s = 0 ( n 1 ) 2 B s , k ( t ^ ) s = ( n + 1 ) / 2 n B s , k ( t ^ ) ,   n   is   odd s = 0 ( n 2 ) 2 B s , k ( t ^ ) s = ( n + 2 ) / 2 n B s , k ( t ^ ) , n   is   even
Therefore, an adaptive BSNN is proposed. Experiments and simulations carried out to verify the performance of the BSNN are discussed in the following section.

4. Results of Experiment and Simulation

Experiments and simulations were carried out to validate the performance of the proposed BSNN-based NFTSM control method in this section. The experimental system is shown in Figure 5. Comparisons are made between the results of the conventional PID method, the TSM in reference [19], the NTSM in reference [15], the NFTSM in this paper, the RBF neural network-based NFTSM. The parameters of the MSA control system are listed in Table 1.
In this system, the BSNN has two inputs and one output. Each of the inputs corresponds to three B-spline activation functions which have n + 1 = 7 control points separately. The internal knot vector of B-Spline curve is [ 1 1 1 2 3 4 5 6 6 6 ] and the degree of B-spline curve is 2. Moreover, the learning rate η = 0.6 , the inertial coefficient α = 0.05 and c = 5 . l d = 7 is adopted for NFTSM control method.

4.1. Description of the Coupling Effect

In this section, the amplitude of the elevation gimbal is 0 degrees in the first 47.9 s (Stage A); the gimbal is assigned to perform a sinusoidal motion 0.5 × sin ( 0.334 × π × t ) from 47.9 s to 89.8 s (Stage B); the gimbal is assigned to perform a sinusoidal motion 1.0 × sin ( 0.1336 × π × t ) from 89.8 s to 200 s (Stage C), as shown in Figure 6. In order to describe the coupling effect caused by the elevation and vehicle movement, there is no control and no input applied to the azimuth gimbal. The vehicle is assigned to move with speed of 0.5 ° / s . Simulation and experiment results are demonstrated in Figure 7.
As shown in Figure 7, the azimuth gimbal has a large angular output in both the simulation and experiment. The azimuth gradually increases due to the motion of the elevation gimbal. Apparently, the results are not what the controller engineer expected. Therefore, a decoupling controller is required.

4.2. Description of De-Coupling Effect

In order to validate the decoupling ability of the proposed method, the azimuth gimbal is assigned to perform a sinusoidal motion with amplitude 0 in the first 48.8 s (Stage A); a sinusoidal motion 1.0 × sin ( 0.164 × π × t ) is applied in the following period (Stage B), as shown in Figure 8. The elevation gimbal is still assigned to perform a sinusoidal motion as in Section 4.1.

4.2.1. Simulation Results

The simulation results are shown in Figure 9.
From the simulation results, the following conclusions can be obtained: the conventional PID method has the worst performance, the peak to peak (pp) value of tracking error is 0.131 ° , the tracking error at time 48.8 s is 0.058 ° ; while the tracking error Δ θ a of the conventional TSM method is 1.75 × 10 3 ( ° ) and the tracking error at time 48.8 s is 26.1 × 10 3 ( ° ) ; the tracking error Δ θ a of the NTSM method in reference is 6.9 × 10 4 ( ° ) , the tracking error is 20.0 × 10 3 ( ° ) ; the tracking error Δ θ a of NFTSM is 6.7 × 10 5 ( ° ) , and the tracking error at time 48.8 s is 18.0 × 10 3 ( ° ) ; for RBFNN based NFTSM the tracking error Δ θ a is 7.3 × 10 5 ( ° ) and the tracking error at time 48.8 s is 23.0 × 10 3 ( ° ) ; the best performance corresponds to the BSNN-based NFTSM with a tracking error of 1.1 × 10 5 ( ° ) and the tracking error at time 48.8 s of 7.8 × 10 3 ( ° ) . At time 48.8 s, the azimuth gimbal begins to perform the assigned sinusoidal motion from a still state. The simulation results show that the proposed BSNN has a better adaptive capability and the proposed controller shows better decoupling effect than other controllers.
The forms of the B-spline curves are shown in Figure 10. It can be seen from the figure that both the positions and the shapes of the B-spline curves are changed.

4.2.2. Experiment Results

Decoupling effect experiments are carried out in this section. The results are shown in Figure 11. As shown in Figure 11, the proposed BSNN-based NFTSM performs better than the other three methods. The pp value of Δ θ a is 0.162 ° when the PID controller is applied; the pp value of Δ θ a is 0.12 ° when the NFTSM controller is applied; the pp value of Δ θ a is 0.077 ° when the RBFNN- based NFTSM is applied; the pp value of Δ θ a is 0.051 ° when the BSNN-based NFTSM is applied.

4.3. Disturbance Rejecting Ability Simulations and Experimental Results

In order to validate the disturbance rejecting ability of the proposed method, pulse disturbances of 50 N·m and 100 N·m are added at time 70 s and 130 s, respectively, as shown in Figure 12.

4.3.1. Simulation Results

The disturbance rejecting ability simulation results are shown in Figure 13. From the simulation results in Figure 13, the pp values of Δ θ a are 0.5 ° and 1.0 ° when the disturbances are added to the conventional PID controller at time 70 s and 130 s, respectively, while the tracking errors of the conventional TSM method are 2.0 × 10 3 ( ° ) and 0.48 ° at time 70 s and 130 s, respectively; the tracking errors Δ θ a of the NTSM method in reference are 9.9 × 10 4 ( ° ) and 0.44 ° ; the pp values of Δ θ a are 8.6 × 10 5 ( ° ) and 0.42 ° at time 70 s and 130 s when NFTSM is applied, respectively.
However, when the RBFNN-based NFTSM is applied, the pp values of Δ θ a are 7.7 × 10 5 ( ° ) and 17 × 10 3 ( ° ) at time 70 s and 130 s, respectively; the pp values of Δ θ a are 4.0 × 10 5 ( ° ) and 6.2 × 10 5 ( ° ) at time 70 s and 130 when the BSNN-based NFTSM is applied, respectively. Therefore, it is apparent that the neural network is good at nonlinear approximation, and the proposed BSNN performs better. The forms of B-spline curves are shown in Figure 14.

4.3.2. Experimental Results

Four different methods are applied to the control system of MSA. From the experimental results in Figure 15, the following conclusions can be obtained: the pp values of Δ θ a are 0.52 ° and 1.01 ° when the disturbances are added to PID controller at time 70 s and 130 s, respectively; whereas the tracking values for NFTSM are 0.072 ° and 0.14 ° at time 70 s and 130 s, respectively; the tracking values for RBFNN based NFTSM are 0.03 ° and 0.082 ° at time 70 s and 130 s, respectively; when BSNN based NFTSM is applied, the pp values of are 0.013 ° at time 70 s and 0.059 ° at time 130 s.
According to the simulations and experiments, it can be seen that the proposed BSNN based NFTSM is good at nonlinear approximate and has strong self-adaptability.

5. Conclusions

The paper focuses on the inertial sensor-based two gimbal mobile satcom antenna. In order to obtain high quality communications, the MSA should point to the specific satellite when the carrier is moving. The dynamic model of MSA is established based on traditional Newton-Euler method and the corresponding control block diagram is built. In this paper, the non-singular fast terminal sliding mode control is adopted and developed to increase the line of sight stabilization accuracy. Meanwhile, the features of existence and convergence in finite time are proved. Then a neural network is employed to approximate the nonlinear item in the system. In addition, a novel BSNN is proposed and used in this paper. A brief study of B-spline basis and B-spline function is also carried out, then the computational function used to obtain the arbitrary point on the curve is derived. The B-spline function is reformed to enhance its adaptive capacity. To validate the effectiveness of the proposed NFTSM and BSNN, simulations and experiments are conducted. Results of different methods, including PID, NFTSM, NFTSM-RBF, NFTSM-BSNN, are compared in this paper. It is shown that the proposed method has better decoupling effects and disturbance rejecting ability than the others. Because the B-spline curve has an excellent ability called “local control”, it can be used to approximate arbitrarily shaped curves. According to the analysis in this paper, we can conclude that the proposed BSNN is good at nonlinear approximation owing to its local features. The robustness of the system can also be improved by applying this method.

Acknowledgments

This work was supported in part by the Natural Science Foundation of China under Grant 61233005 and National 973 Project under Grant 2014CB744200. The author would like to thank the reviewers whose critical analysis and suggestions brought insight into the final editing of this paper.

Author Contributions

Xiaolei Zhang, Yan Zhao, Kai Guo, Gaoliang Li and Nianmao Deng conceived the study idea and designed the experiments. Xiaolei Zhang and Yan Zhao developed the simulation platform. Xiaolei Zhang and Kai Guo performed the experiments. Xiaolei Zhang drafted the manuscript. All authors read and approved the manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The Terminal sliding mode surface (8) can be written as follows:
s = e ^ s i g n ( e ^ ) + 1 α e ^ r + 1 s i g n ( e ^ ) + 1 β e ^ ˙ p + 1 s i g n ( e ^ ˙ ) = 0
The following formula can be obtained:
s i g n ( e ^ ) = s i g n ( e ^ ˙ )
Then Equation (A1) can be further simplified as:
e ^ + 1 α e ^ r + 1 = 1 β e ^ ˙ p + 1
It is apparent that:
( β e ^ + β α e ^ r + 1 ) 1 p + 1 = e ^ ˙ = e ^ ˙ s i g n ( e ^ ˙ )
Then the following formula can be obtained:
( β e ^ + β α e ^ r + 1 ) 1 p + 1 s i g n ( e ^ ) = e ^ ˙
If e ^ 0 ,
( β e ^ + β α e ^ r + 1 ) 1 p + 1 = e ^ ˙
Else if e ^ < 0 ,
( β ( e ^ ) + β α ( e ^ ) r + 1 ) 1 p + 1 = e ^ ˙
Equation (A5) can be integrated and written as follows:
0 e ^ ( t 0 ) d ( e ^ ) ( β ( e ^ ) + β α ( e ^ ) r + 1 ) 1 p + 1 + 0 e ^ ( t 0 ) d e ^ ( β e ^ + β α e ^ r + 1 ) 1 p + 1 = 0 t d t
0 | e ^ ( t 0 ) | d e ^ ( β e ^ ) 1 p + 1 0 t d t
Then the convergence time can be obtained:
t p + 1 p β 1 p + 1 | e ^ ( t 0 ) | p p + 1
This ends the proof.

References

  1. Densmore, A.C.; Jamnejad, V. A satellite-tracking k- and ka-band mobile vehicle antenna system. IEEE Trans. Veh. Technol. 1993, 42, 502–513. [Google Scholar] [CrossRef]
  2. Cheng, W.; Minzhe, T.; Weizhou, S. An h2/h control design for mobile satcom antenna servo systems. In Proceedings of the 2016 35th Chinese Control Conference (CCC), Chengdu, China, 27–29 July 2016; Volume 2016, pp. 3023–3028. [Google Scholar]
  3. Marsh, E.A. Inertially Stabilized Platforms for Satcom on-the-Move Applications: A Hybrid Open/Closed-Loop Antenna Pointing Strategy. Master’s degree, Massachusetts Institute of Technology, Cambridge, MA, USA, 2008. [Google Scholar]
  4. Debruin, J. Control systems for mobile satcom antennas. IEEE Control Syst. 2008, 28, 86–101. [Google Scholar] [CrossRef]
  5. Tian, F.; Wang, R.; Gao, F.; Yao, M. A beam stabilization algorithm based on nonlinear observer for low cost satcom-on-the-move. China Commun. 2014, 11, 16–22. [Google Scholar] [CrossRef]
  6. Rue, A.K. Stabilization of precision electrooptical pointing and tracking systems. IEEE Trans. Aerosp. Electr. Syst. 1969, AES-5, 805–819. [Google Scholar]
  7. Shuang, C.; Ke, D.; Shang, W. Modeling analysis on the gyro stabilized. Sci. Technol. Rev. 2011, 241, 1606–1615. [Google Scholar]
  8. Fang, J.; Yin, R.; Lei, X. An adaptive decoupling control for three-axis gyro stabilized platform based on neural networks. Mechatronics 2015, 27, 38–46. [Google Scholar] [CrossRef]
  9. Ito, K.; Maebashi, W.; Ikeda, J.; Iwasaki, M. Fast and precise positioning of rotary table systems by feedforward disturbance compensation considering interference force. In Proceedings of the IECON 2011 37th Annual Conference on IEEE Industrial Electronics Society, Melbourne, Australia, 7–10 November 2011; pp. 3382–3387. [Google Scholar]
  10. Řezáč, M.; Hurák, Z. Vibration rejection for inertially stabilized double gimbal platform using acceleration feedforward. In Proceedings of the IEEE International Conference on Control Applications (CCA), Denver, CO, USA, 28–30 September 2011; pp. 363–368. [Google Scholar]
  11. Mu, Q.; Liu, G.; Lei, X. A RBFNN-based adaptive disturbance compensation approach applied to magnetic suspension inertially stabilized platform. Math. Probl. Eng. 2014, 2014, 1–9. [Google Scholar] [CrossRef]
  12. Moorty, J.A.R.K.; Marathe, R.; Sule, V.R. H control law for line-of-sight stabilization for mobile land vehicles. Opt. Eng. 2002, 41, 41–45. [Google Scholar]
  13. Li, S.; Zhong, M.; Zhao, Y. Estimation and compensation of unknown disturbance in three-axis gyro-stabilized camera mount. Trans. Inst. Meas. Control 2014. [Google Scholar] [CrossRef]
  14. Ma, D.; Lin, H. Chattering-free nonsingular fast terminal sliding-mode control for permanent magnet synchronous motor servo system. In Proceedings of the 35th Chinese Control Conference (CCC), Chengdu, China, 27–29 July 2016. [Google Scholar]
  15. Feng, Y.; Yu, X.; Man, Z. Non-singular terminal sliding mode control of rigid manipulators. Automatica 2002, 38, 2159–2167. [Google Scholar] [CrossRef]
  16. Zak, M. Terminal attractors for addressable memory in neural networks. Phys. Lett. A 1988, 133, 18–22. [Google Scholar] [CrossRef]
  17. Wu, Y.; Yu, X.; Man, Z. Terminal sliding mode control design for uncertain dynamic systems. Syst. Control Lett. 1998, 34, 281–287. [Google Scholar] [CrossRef]
  18. Madani, T.; Daachi, B.; Djouani, K. Non-singular terminal sliding mode controller: Application to an actuated exoskeleton. Mechatronics 2016, 33, 136–145. [Google Scholar] [CrossRef]
  19. Yu, S.; Yu, X.; Stonier, R. Continuous finite-time control for robotic manipulators with terminal sliding modes. In Proceedings of the 6th International Conference of Information Fusion, Budapest, Hungary, 15–17 October 2003; pp. 1957–1964. [Google Scholar]
  20. Ding, F.; Cui, Y.; Wang, C.; Zhang, X. Course control of air cushion vessel based on terminal sliding mode control with rbf neural network. In Proceedings of the 35th Chinese Control Conference (CCC), Chengdu, China, 27–29 July 2016; pp. 3496–3502. [Google Scholar]
  21. Shtessel, Y.; Edwards, C.; Fridman, L.; Levant, A. Sliding Mode Control and Observation; Springer: New York, NY, USA, 2015; pp. 213–249. [Google Scholar]
  22. Levant, A. Universal single-input-single-output (siso) sliding-mode controllers with finite-time convergence. IEEE Trans. Autom. Control 2001, 46, 1447–1451. [Google Scholar] [CrossRef]
  23. Aghababa, M.P. Design of a chatter-free terminal sliding mode controller for nonlinear fractional-order dynamical systems. Int. J. Control 2013, 86, 1744–1756. [Google Scholar] [CrossRef]
  24. Tran, M.D.; Kang, H.J. Adaptive terminal sliding mode control of uncertain robotic manipulators based on local approximation of a dynamic system. Neurocomputing 2017, 228, 231–240. [Google Scholar] [CrossRef]
  25. Jia, T.; Kang, G. An rbf neural network-based nonsingular terminal sliding mode controller for robot manipulators. In Proceedings of the 2012 Third International Conference on Intelligent Control and Information Processing (ICICIP), Dalian, China, 15–17 July 2012; pp. 72–76. [Google Scholar]
  26. Wang, H.; Yang, Z.; Zhou, Z. Rbf-based terminal sliding mode control for a class of underactuated mechanical system. In Proceedings of the 2016 Chinese Control and Decision Conference (CCDC), Yinchuan, China, 28–30 May 2016; pp. 664–667. [Google Scholar]
  27. Fang, Y.; Liu, L.; Li, J.; Xu, Y. Decoupling control based on terminal sliding mode and wavelet network for the speed and tension system of reversible cold strip rolling mill. Int. J. Control 2015, 88, 1–17. [Google Scholar] [CrossRef]
  28. Lin, C.K. Nonsingular terminal sliding mode control of robot manipulators using fuzzy wavelet networks. IEEE Trans. Fuzzy Syst. 2006, 14, 849–859. [Google Scholar] [CrossRef]
  29. Ferch, M.; Zhang, J.; Knoll, A. Robot skill transfer based on B-Spline fuzzy controllers for force-control tasks. In Proceedings of the 1999 IEEE International Conference on Robotics and Automation, Detroit, MI, USA, 10–15 May 1999; Volume 1172, pp. 1170–1175. [Google Scholar]
  30. Yiu, K.F.; Wang, S.; Teo, K.L.; Tsoi, A.C. Nonlinear system modeling via knot-optimizing B-Spline networks. IEEE Trans. Neural Netw. 2001, 12, 1013–1022. [Google Scholar] [PubMed]
  31. Cabrita, C.; Ruano, A.E.; Fonseca, C.M. Single and multi-objective genetic programming design for B-Spline neural networks and Neuro-Fuzzy systems. In Proceedings of the Ifac Workshop on Advanced Fuzzy-Neural Control, Valencia, Spain, 15–16 October 2001; pp. 75–80. [Google Scholar]
  32. Leandro, D.S.C.; Guerra, F.A. B-spline neural network design using improved differential evolution for identification of an experimental nonlinear process. Appl. Soft Comput. 2008, 8, 1513–1522. [Google Scholar]
  33. Qamar, S.; Khan, L.; Qamar, Z. Online adaptive full car active suspension control using b-spline fuzzy-neural network. In Proceedings of the 2013 11th International Conference on Frontiers of Information Technology (FIT), Islamabad, Pakistan, 16–18 December 2013; pp. 205–210. [Google Scholar]
  34. Qamar, S.; Khan, L.; Ali, S. Adaptive B-Spline based Neuro-Fuzzy control for full car active suspension system. Middle East J. Sci. Res. 2013, 16, 1348–1360. [Google Scholar]
  35. Wang, C.H.; Wang, W.Y.; Lee, T.T.; Tseng, P.S. Fuzzy B-Spline membership function and its applications in Fuzzy-Neural control. IEEE Trans. Syst. Man Cybern. 1995, 25, 841–851. [Google Scholar] [CrossRef]
  36. Mar, J.; Lin, F.J.; Lin, H.T.; Hsu, L.C. Evolutionary learning of BMF Fuzzy-Neural networks using a reduced-form genetic algorithm. IEEE Trans. Syst. Man Cybern. Part B Cybern. 2003, 33, 966–976. [Google Scholar]
  37. Mao, J.; Wang, Y.; Sun, W. Remote sensing images classification using fuzzy B-Spline function neural network. J. Electr. Meas. Instrom. 2002, 3, 2159–2163. [Google Scholar]
  38. Coelho, L.D.S.; Pessôa, M.W. Nonlinear identification using a B-Spline neural network and chaotic immune approaches. Mech. Syst. Signal Proc. 2009, 23, 2418–2434. [Google Scholar] [CrossRef]
  39. Boor, C.D. A practical Guide to Splines/Arl De Boor; Wiley: Hoboken, NJ, USA, 1985. [Google Scholar]
  40. B-Splines. Available online: http://web.mit.edu/hyperbook/Patrikalakis-Maekawa-Cho/node16.html (accessed on 1 January 2017).
Figure 1. Coordinate definition and transformation.
Figure 1. Coordinate definition and transformation.
Sensors 17 00978 g001
Figure 2. Control block diagram of azimuth gimbal.
Figure 2. Control block diagram of azimuth gimbal.
Sensors 17 00978 g002
Figure 3. (a) Order 3 B-spline basis function in which the internal knot vector is [ 1 2 3 4 5 6 ] ; (b) Order 3 B-spline basis function that has multiple knots in which the internal knot vector is [ 1 1 1 2 3 4 5 6 6 6 ] ; (c) Order 3 B-spline functions with multiple knots; the internal knot vector is [ 1 1 1 2 3 4 5 6 6 6 ] ; (d) Amplitude standardization of c; the internal knot vector is [ 1 1 1 2 3 4 5 6 6 6 ] .
Figure 3. (a) Order 3 B-spline basis function in which the internal knot vector is [ 1 2 3 4 5 6 ] ; (b) Order 3 B-spline basis function that has multiple knots in which the internal knot vector is [ 1 1 1 2 3 4 5 6 6 6 ] ; (c) Order 3 B-spline functions with multiple knots; the internal knot vector is [ 1 1 1 2 3 4 5 6 6 6 ] ; (d) Amplitude standardization of c; the internal knot vector is [ 1 1 1 2 3 4 5 6 6 6 ] .
Sensors 17 00978 g003
Figure 4. The general structure of BSNN.
Figure 4. The general structure of BSNN.
Sensors 17 00978 g004
Figure 5. Experimental two axis MSA system.
Figure 5. Experimental two axis MSA system.
Sensors 17 00978 g005
Figure 6. The sinusoidal motion of the elevation gimbal.
Figure 6. The sinusoidal motion of the elevation gimbal.
Sensors 17 00978 g006
Figure 7. The coupling effect of the azimuth gimbal: (a) simulated coupling effect result; (b) experimental coupling effect result.
Figure 7. The coupling effect of the azimuth gimbal: (a) simulated coupling effect result; (b) experimental coupling effect result.
Sensors 17 00978 g007
Figure 8. The sinusoidal motion of the azimuth gimbal.
Figure 8. The sinusoidal motion of the azimuth gimbal.
Sensors 17 00978 g008
Figure 9. Comparative azimuth gimbal absolute angle simulation results with four different methods: (a) PID method; (b) TSM method in reference [19]; (c) NTSM method in reference [15]; (d) the proposed NFTSM method; (e) NFTSM-RBFNN method; (f) the proposed NFTSM-BSNN method.
Figure 9. Comparative azimuth gimbal absolute angle simulation results with four different methods: (a) PID method; (b) TSM method in reference [19]; (c) NTSM method in reference [15]; (d) the proposed NFTSM method; (e) NFTSM-RBFNN method; (f) the proposed NFTSM-BSNN method.
Sensors 17 00978 g009
Figure 10. Form change of the B-spline activation function: (a) initial forms of the activation curves; (b) final forms of the activation curves after training.
Figure 10. Form change of the B-spline activation function: (a) initial forms of the activation curves; (b) final forms of the activation curves after training.
Sensors 17 00978 g010
Figure 11. Comparative azimuth gimbal absolute angle experimental results with four different methods: (a) PID method; (b) NFTSM method; (c) NFTSM-RBFNN method; (d) the proposed method.
Figure 11. Comparative azimuth gimbal absolute angle experimental results with four different methods: (a) PID method; (b) NFTSM method; (c) NFTSM-RBFNN method; (d) the proposed method.
Sensors 17 00978 g011
Figure 12. Pulse disturbances added to the system.
Figure 12. Pulse disturbances added to the system.
Sensors 17 00978 g012
Figure 13. Comparative azimuth gimbal absolute angle simulation results with four different methods: (a) PID method; (b) TSM method in reference [19]; (c) NTSM method in reference [15]; (d) the proposed NFTSM method; (e) NFTSM-RBFNN method; (f) the proposed NFTSM-BSNN method.
Figure 13. Comparative azimuth gimbal absolute angle simulation results with four different methods: (a) PID method; (b) TSM method in reference [19]; (c) NTSM method in reference [15]; (d) the proposed NFTSM method; (e) NFTSM-RBFNN method; (f) the proposed NFTSM-BSNN method.
Sensors 17 00978 g013
Figure 14. Form change of the B-spline activation function: (a) initial forms of the curves; (b) final forms of the curves after training.
Figure 14. Form change of the B-spline activation function: (a) initial forms of the curves; (b) final forms of the curves after training.
Sensors 17 00978 g014
Figure 15. Comparative azimuth gimbal absolute angle experimental results with four different methods: (a) PID method; (b) NFTSM method; (c) NFTSM-RBFNN method; (d) the proposed method.
Figure 15. Comparative azimuth gimbal absolute angle experimental results with four different methods: (a) PID method; (b) NFTSM method; (c) NFTSM-RBFNN method; (d) the proposed method.
Sensors 17 00978 g015
Table 1. Parameters of MSA system.
Table 1. Parameters of MSA system.
ParametersValuesUnitParametersValuesUnitParametersValuesUnit
J ex 0.51 kg · m 2 J ey 0.54 kg · m 2 J ez 0.59 kg · m 2
J ax 1.49 kg · m 2 J ay 1.56 kg · m 2 J az 1.42 kg · m 2
J ad 0.96 kg · m 2 J af 1.37 kg · m 2 N a 2.8 N · m / A
k at 0.42 k ab 0.42 R aa 1.23 Ω
f ad 0.3 f af 0.5

Share and Cite

MDPI and ACS Style

Zhang, X.; Zhao, Y.; Guo, K.; Li, G.; Deng, N. An Adaptive B-Spline Neural Network and Its Application in Terminal Sliding Mode Control for a Mobile Satcom Antenna Inertially Stabilized Platform. Sensors 2017, 17, 978. https://doi.org/10.3390/s17050978

AMA Style

Zhang X, Zhao Y, Guo K, Li G, Deng N. An Adaptive B-Spline Neural Network and Its Application in Terminal Sliding Mode Control for a Mobile Satcom Antenna Inertially Stabilized Platform. Sensors. 2017; 17(5):978. https://doi.org/10.3390/s17050978

Chicago/Turabian Style

Zhang, Xiaolei, Yan Zhao, Kai Guo, Gaoliang Li, and Nianmao Deng. 2017. "An Adaptive B-Spline Neural Network and Its Application in Terminal Sliding Mode Control for a Mobile Satcom Antenna Inertially Stabilized Platform" Sensors 17, no. 5: 978. https://doi.org/10.3390/s17050978

APA Style

Zhang, X., Zhao, Y., Guo, K., Li, G., & Deng, N. (2017). An Adaptive B-Spline Neural Network and Its Application in Terminal Sliding Mode Control for a Mobile Satcom Antenna Inertially Stabilized Platform. Sensors, 17(5), 978. https://doi.org/10.3390/s17050978

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop