# Iterative Sliding Mode and Increment Feedback Attitude Control for On-Orbit Capturing Process of Spacecraft

Authors: ChunWu Yin; YouJun Ding; LiJun Song; ZhaoQiang Li; Xing He
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

According to the characteristics of spacecraft capturing noncooperative targets in orbit, an increment feedback controller based on nonlinear iterative sliding mode is presented. Firstly, the attitude tracking error equation is established, and then, an increment feedback control law based on bounded iterative sliding modes is proposed, which does not need to estimate the uncertain moment of inertia and external disturbances. For comparing, an adaptive sliding mode controller has been designed in the paper. Some numerical simulations have been given in the presence of spacecraft on-orbit capturing noncooperative target, and the simulation results show that the increment feedback controller has strong robustness to the unknown parametric variations and external disturbances and has a smaller control input torque in control process.


---

## Body

## 1. Introduction

With the increasing demand of communication and remote inspection, countless of satellites have been launched to space, and a large number of spacecraft debris have been generated simultaneously. In order to promote the sustainable development of the space and to ensure the safety of spacecraft on-orbit, the technology of spacecraft on-orbit acquisition has been proposed and studied by several researchers. On-orbit capturing technology has wide range of applications, such as cleaning space debris and repairing satellite, so the attitude stability controller of spacecraft capturing noncooperation target has an important practical application value. In the process of spacecraft on-orbit capturing noncooperation target, the moment of inertia of spacecraft will change due to the extension and recovery of manipulator. Thus, the essence of spacecraft stability control is attitude stability control under unknown inertial disturbance in the noncooperative process of spacecraft acquisition on-orbit.The problem of spacecraft attitude control has been researched during the past decades, and most of control methods both linear and nonlinear have been provided.The linear system linearizes the dynamics and kinematics equations of spacecraft at the equilibrium point, and the attitude stability controller for linear system is designed by classical linear control methods such as LQR control, H2 control, H∞ control, and H2 /H∞ control based on representation form of quaternion; the spacecraft attitude motion equations are linearized by Silani and Lovera [1] and Yang [2, 3]. Based on representation form of Euler angle, the spacecraft attitude motion equations are linearized by Yang and Sun [4]. The robustness and efficiency of all the classical linear control methods is compared through numerical simulations and practical applications by Silani and Lovera [1], Won [5], and Yang and Sun [4]. Nonlinear control such as PID-like control [6], backstepping control [7], sliding mode control [8–10], adaptive control [11, 12], inverse optimal control [13], fuzzy control [14], neural-network control [15], and other control methods [16] have been employed to design attitude stability controller, and Lyapunov theory has been utilized to analyze the stability of the whole system.Some controllers only consider external disturbance, while others only compensate parametric disturbances. However, the focus of current research is on the design of attitude stabilization controllers with both parametric uncertainties and bounded external disturbances. Generally, Kalman filter [17–19], adaptive controller, or fuzzy/neural-network approximator are used to estimate the disturbance of moment of inertia, and the sliding mode controller is applied to counteract external disturbance.Considering the structure of adaptive controllers or sliding mode controllers, the numerical simulation results based on these controllers are analyzed, it is found that these controllers are complex, and the control input torque is large. The control input torque of the sliding mode adaptive controller designed by Yeh [10] is more than 2000 N·m. The control torque of the controller designed by Silani and Lovera [1] is between −50 N·m and 50 N·m, and these control inputs do not satisfy the physical constraint of spacecraft.Due to the physical and energy constraints, the actual system has state constraints and input saturation. Ignoring these restrictions, some undesired dynamic behaviors would occur in the practical control process. Therefore, the robust controller design with state and input constraint is still a challenge in the field of nonlinear system stability control [20, 21]. In this paper, the design of robust attitude controller is discussed under the conditions of inertial matrix uncertainty, unknown disturbance, and control input constraints.In this paper, an attitude stability controller based on iterative sliding mode and increment feedback is designed in the paper. In Section2, the attitude tracking error equation is derived from attitude dynamics and kinematics equations of rigid spacecraft; in Section 3, the iterative sliding mode and increment feedback controller for attitude tracking error system are proposed. Based on Lyapunov stability theory, the closed-loop system is proved to be exponentially convergent. In Section 4, the adaptive sliding mode attitude stability controller is designed for comparison. In Section 5, the control performance is verified by numerical simulations. Finally, the proposed method is demonstrated by numerical simulation in Section 6.


## 2. Attitude Dynamics of a Rigid Spacecraft

### 2.1. Dynamics Equation

The attitude dynamic motion of rigid spacecraft with the external disturbance can be described by the following set of differential equation:(1)Jω˙=−ω×Jω+τu+τd,where J∈R3×3 denotes the inertia matrix of the body and satisfies J=JT>0, ω∈R3 is the angular velocity vector of the body in the body-fixed frame, τu∈R3 is the control torque vector, and τd∈R3 is the bounded external disturbance vector such as gravity and radiation. Moreover, the matrix ω×∈R3×3 denotes a skew symmetric matrix given by(2)ω×=0ω3t−ω2t−ω3t0ω1tω2t−ω1t0.


### 2.2. Kinematics Equation

The orientation of rigid spacecraft can be described by several parameterizations with a special orthogonal group. In this paper, the orientation of the rigid body is parameterized by the yaw-pitch-rollθ3−θ2−θ1Euler angle. The attitude kinematics equation of rigid spacecraft is described as follows:(3)ω=Rθθ˙−ωcθ,where θ=θ1,θ2,θ3Τ is the attitude Euler angles and ω0 is the orbital rate. Also,(4)Rθ=10−sinθ20cosθ1sinθ1cosθ20−sinθ1cosθ1cosθ2,ωcθ=ω0cosθ2sinθ3cosθ1sinθ3+sinθ1sinθ2sinθ3−sinθ1cosθ3+cosθ1sinθ2sinθ3.Taking the time derivative of (3), the following is obtained:(5)ω˙=Rθθ¨+ddtRθθ˙−ddtωcθ.Substituting (3) and (5) into (1) and premultiplying (5) by the matrix RTθ, the following is obtained:(6)Mθθ¨+Cθ,θ˙θ˙+Gθ,θ˙=u+d,where(7)Mθ=RTθJRθ,Cθ,θ˙=RTθJddtRθ+RTθω×JRθ,Gθ,θ˙=−RTθJddtωcθ−RTθω×Jωcθ,u=RTθτa,d=RTθτd.DenoteJ0 as nominal moment of inertia and ΔJas inertia uncertainties, then the real moment of inertia is J=J0+ΔJ.Accordingly, there are(8)Mθ=RTθJ0Rθ+RTθΔJRθ=M0θ+ΔMθ,(9)Cθ,θ˙=C0θ,θ˙+ΔCθ,θ˙.Substituting (8) and (9) into (6) yields(10)M0θθ¨+C0θ,θ˙θ˙=u+d−ΔMθθ¨−ΔCθ,θ˙θ˙−Gθ,θ˙.Let(11)Y=−ΔMθθ¨−ΔCθ,θ˙θ˙−Gθ,θ˙.Then, equation (10) becomes(12)θ¨=M0−1θu+d+Y−C0θ,θ˙θ˙.


### 2.3. Attitude Tracking Error Equation

Denote the desired attitude angle asθr∈C2; correspondingly, its desired velocity vector is θ˙r and acceleration vector is θ¨r. Denote attitude tracking error as(13)e1=θ−θr.Lete2=e˙1=θ˙−θ˙r, taking the time derivative of e2 with (12):(14)e˙2=θ¨−θ¨r=M0−1u+d+Y−C0θ˙−θ¨r=M0−1u+d+Y−C0e2−C0θ˙r−θ¨r.Then, it will obtain the attitude tracking error equations as(15)e˙1=e2,e˙2=−M0−1C0e2−M0−1C0θ˙r−θ¨r+M0−1u+d+Y.The attitude control objective of spacecraft on-orbit capturing noncooperative target is to design the control law of system (15) to ensure global asymptotic stability of the closed-loop system in the presence of uncertain moment inertia matrix ΔJ, unknown disturbance d, and state and input constraints.


## 2.1. Dynamics Equation

The attitude dynamic motion of rigid spacecraft with the external disturbance can be described by the following set of differential equation:(1)Jω˙=−ω×Jω+τu+τd,where J∈R3×3 denotes the inertia matrix of the body and satisfies J=JT>0, ω∈R3 is the angular velocity vector of the body in the body-fixed frame, τu∈R3 is the control torque vector, and τd∈R3 is the bounded external disturbance vector such as gravity and radiation. Moreover, the matrix ω×∈R3×3 denotes a skew symmetric matrix given by(2)ω×=0ω3t−ω2t−ω3t0ω1tω2t−ω1t0.


## 2.2. Kinematics Equation

The orientation of rigid spacecraft can be described by several parameterizations with a special orthogonal group. In this paper, the orientation of the rigid body is parameterized by the yaw-pitch-rollθ3−θ2−θ1Euler angle. The attitude kinematics equation of rigid spacecraft is described as follows:(3)ω=Rθθ˙−ωcθ,where θ=θ1,θ2,θ3Τ is the attitude Euler angles and ω0 is the orbital rate. Also,(4)Rθ=10−sinθ20cosθ1sinθ1cosθ20−sinθ1cosθ1cosθ2,ωcθ=ω0cosθ2sinθ3cosθ1sinθ3+sinθ1sinθ2sinθ3−sinθ1cosθ3+cosθ1sinθ2sinθ3.Taking the time derivative of (3), the following is obtained:(5)ω˙=Rθθ¨+ddtRθθ˙−ddtωcθ.Substituting (3) and (5) into (1) and premultiplying (5) by the matrix RTθ, the following is obtained:(6)Mθθ¨+Cθ,θ˙θ˙+Gθ,θ˙=u+d,where(7)Mθ=RTθJRθ,Cθ,θ˙=RTθJddtRθ+RTθω×JRθ,Gθ,θ˙=−RTθJddtωcθ−RTθω×Jωcθ,u=RTθτa,d=RTθτd.DenoteJ0 as nominal moment of inertia and ΔJas inertia uncertainties, then the real moment of inertia is J=J0+ΔJ.Accordingly, there are(8)Mθ=RTθJ0Rθ+RTθΔJRθ=M0θ+ΔMθ,(9)Cθ,θ˙=C0θ,θ˙+ΔCθ,θ˙.Substituting (8) and (9) into (6) yields(10)M0θθ¨+C0θ,θ˙θ˙=u+d−ΔMθθ¨−ΔCθ,θ˙θ˙−Gθ,θ˙.Let(11)Y=−ΔMθθ¨−ΔCθ,θ˙θ˙−Gθ,θ˙.Then, equation (10) becomes(12)θ¨=M0−1θu+d+Y−C0θ,θ˙θ˙.


## 2.3. Attitude Tracking Error Equation

Denote the desired attitude angle asθr∈C2; correspondingly, its desired velocity vector is θ˙r and acceleration vector is θ¨r. Denote attitude tracking error as(13)e1=θ−θr.Lete2=e˙1=θ˙−θ˙r, taking the time derivative of e2 with (12):(14)e˙2=θ¨−θ¨r=M0−1u+d+Y−C0θ˙−θ¨r=M0−1u+d+Y−C0e2−C0θ˙r−θ¨r.Then, it will obtain the attitude tracking error equations as(15)e˙1=e2,e˙2=−M0−1C0e2−M0−1C0θ˙r−θ¨r+M0−1u+d+Y.The attitude control objective of spacecraft on-orbit capturing noncooperative target is to design the control law of system (15) to ensure global asymptotic stability of the closed-loop system in the presence of uncertain moment inertia matrix ΔJ, unknown disturbance d, and state and input constraints.


## 3. Adaptive Sliding Mode Attitude Stable Controller Design

### 3.1. Controller Design

When the parameter perturbation satisfies the matching conditions, the sliding mode control has good robustness due to the invariance of the system. Increment feedback control has simple structure and small control torque because the parameter perturbance and external disturbance are directly contained into the controller. Therefore, a nonlinear iterative sliding mode attitude stability control law is designed based on increment feedback control in this section. Design PID attitude controller (16) for system (15):(16)u=−kde˙1−kpe1−ki∫e1dt,where const kp∈R+,kd∈R+,andki∈R+.Modify the controller (16)(17)u=−kde˙1−kpe1−ki∫e1dt=−kde˙1+kp2kd+kp24kd2−kikde1+kp2kd−kp24kd2−kikde1+kp2kd+kp24kd2−kikdkp2kd−kp24kd2−kikd∫e1dt.Letk1=kp/2kd−kp2/4kd2−ki/kd,k2=kp/2kd+kp2/4kd2−ki/kd and obtain(18)u=−kde˙1+k1e1+k2e1+k1k2∫e1dt=−kde˙1+k1e1+k2e1+k1∫e11dt=−kdddte1+k1∫e1dt+k2e1+k1∫e1dt.Definee0=∫e1dt, then e˙0=e1; equation (15) can be modified as(19)e˙0=e1,e˙1=e2,e˙2=−M0−1C0e2−M0−1C0θ˙r−θ¨r+M0−1u+d+Y.Substitutinge0=∫e1dt in equation (18), the following is obtained:(20)u=−kdddte1+k1e0+k2e1+k1e0.Define variates as follows:(21)σ1=k1e0+e˙0,σ2=k2σ1+σ˙1,where σ1 is a sliding mode on the phase plane e0e˙0, e˙0=−k1e0 as σ1⟶0, e0 will converge exponentially, and σ2 is a sliding mode on the phase planeσ1σ˙1. When σ2⟶0, it can obtain σ˙1=−k2σ1, and σ2 also will converge exponentially, so the stability of system (15) is equal to the stabilization control of σ2.From the above analysis, the iterative design method of spacecraft attitude control can be obtained. Firstly, the sliding modeσ1 based on tracking errors e1 and e˙1 can be designed for system (15) to obtain a new system. Then, another sliding mode σ2 can be designed for the new system based on σ1 and σ˙1.In order to satisfy the constraints of bounded state and control input and to eliminate the uncertainties caused by parameter perturbations and external disturbances, the nonlinear interactive sliding mode is designed as follows:(22)σ1=k1tanhk2e1+e˙1,σ2=k3tanhk4σ1+σ˙1,where const k1,k2,k3,k4>0. tanh⋅ is the hyperbolic tangent function, and it can limit the range of the variable. σ1 is a sliding mode on the phase plane e1e˙1, e˙1=−k1tanhk2e1 as σ1⟶0; when the deviation of e1 is big, e1 will convergence with a fixed rate(i.e., k1), and when the deviation of e1 is small, e1 will converge exponentially; k2 is a compression or tension coefficient, which determines the sliding slope in the origin of phase plane together with k1. σ2 is a sliding mode on the phase planeσ1σ˙1; when σ2⟶0, it can obtain σ˙1=−k3tanhk4σ1; when the deviation of σ1 is big, σ1 will converge with a fixed rate approximately(i.e., k3); if the deviation of σ1 is small, σ1 will converge exponentially; from the definition of σ1 and σ2, we can obtain that the stability control of system (15) is equal to the stability control of σ2.Taking the time derivative ofσ1, give(23)σ˙1=k14ek2e1+e−k2e12Τk2e˙1+e¨1=k14ek2e1+e−k2e12Τk2e˙1−M0−1C0e2−M0−1C0θ˙r−θ¨r+M0−1u+d+Y.Substituting (23) into (22) yields(24)σ2=k3tanhk4k1tanhk2e1+k4e2+k14ek2e1+e−k2e12Τk2e˙1−M0−1C0e2−M0−1C0θ˙r−θ¨r+M0−1u+d+Y.In order to avoid the estimation of the unknown moment of inertia and unknown external disturbance and simplify the implementation of attitude stability controller, the following incremental feedback controller is designed to ensure the stability of zero-order scalar system (24):(25)u˙=−kpσ2−εsatσ2,satx=xmin,ifx<xmin,x,ifxmin≤x≤xmax,xmax,ifx>xmax,where constkp,ε∈R+. xmin and xmax are the minimum value and maximum value of x, respectively. The incremental feedback controller does not need to estimate or cancel unknown moment of inertia and unknown external disturbance but only needs to search the solution of σ2=0 along with the opposite direction ofσ2, which avoids the waste of control input energy and reduces the magnitude of the control input torque effectively.


### 3.2. Analysis of the Stability

Barb˘alat Lemma: if a function f(t) is continuous and differentiable and satisfies the following:(1)
limt⟶∞ft=C C is bounded constant.(2)
f˙t is uniformly continuous (meansf¨t is bounded), then limt⟶∞f˙t=0.Theorem 1.
Under the iterative sliding mode (22) and incremental feedback controller (25), the closed-loop system (15) is globally asymptotic stable.Proof.
According to the above analysis, it is known that the stability control of system (15) is equal to the stability control of σ2.
Considering Lyapunov function asV=1/2σ2Τσ2 and taking the time derivative of V yields(26)V˙=σ2Τσ˙2=σ2Τ∂σ2∂uu˙.
Ignoring the variables of independent ofu in (24) and taking the time derivative of (24) about u, we obtain(27)∂σ2∂u=M0−1.
BecauseM0 is a positive matrix, M0−1 is also a positive matrix. Substituting (25) and (27) into (26) yields(28)V˙=σ2ΤM0−1−kpσ2−εsatσ2=−kpσ2ΤM0−1σ2−εσ2ΤM0−1satσ2≤−kpσ2ΤM0−1σ2≤0.
Then, we can obtainV˙=σ2Τσ˙2≤0, which means V˙≤0 when t≥0. According to the Barbalat lemma, we can obtain σ2⟶0 as t⟶∞, and from the definition of σ2and σ1, we can obtain σ1⟶0 and e1⟶0, respectively. Therefore, system (15) is globally asymptotic stable with the incremental feedback controller (25) and interactive sliding mode (22).


## 3.1. Controller Design

When the parameter perturbation satisfies the matching conditions, the sliding mode control has good robustness due to the invariance of the system. Increment feedback control has simple structure and small control torque because the parameter perturbance and external disturbance are directly contained into the controller. Therefore, a nonlinear iterative sliding mode attitude stability control law is designed based on increment feedback control in this section. Design PID attitude controller (16) for system (15):(16)u=−kde˙1−kpe1−ki∫e1dt,where const kp∈R+,kd∈R+,andki∈R+.Modify the controller (16)(17)u=−kde˙1−kpe1−ki∫e1dt=−kde˙1+kp2kd+kp24kd2−kikde1+kp2kd−kp24kd2−kikde1+kp2kd+kp24kd2−kikdkp2kd−kp24kd2−kikd∫e1dt.Letk1=kp/2kd−kp2/4kd2−ki/kd,k2=kp/2kd+kp2/4kd2−ki/kd and obtain(18)u=−kde˙1+k1e1+k2e1+k1k2∫e1dt=−kde˙1+k1e1+k2e1+k1∫e11dt=−kdddte1+k1∫e1dt+k2e1+k1∫e1dt.Definee0=∫e1dt, then e˙0=e1; equation (15) can be modified as(19)e˙0=e1,e˙1=e2,e˙2=−M0−1C0e2−M0−1C0θ˙r−θ¨r+M0−1u+d+Y.Substitutinge0=∫e1dt in equation (18), the following is obtained:(20)u=−kdddte1+k1e0+k2e1+k1e0.Define variates as follows:(21)σ1=k1e0+e˙0,σ2=k2σ1+σ˙1,where σ1 is a sliding mode on the phase plane e0e˙0, e˙0=−k1e0 as σ1⟶0, e0 will converge exponentially, and σ2 is a sliding mode on the phase planeσ1σ˙1. When σ2⟶0, it can obtain σ˙1=−k2σ1, and σ2 also will converge exponentially, so the stability of system (15) is equal to the stabilization control of σ2.From the above analysis, the iterative design method of spacecraft attitude control can be obtained. Firstly, the sliding modeσ1 based on tracking errors e1 and e˙1 can be designed for system (15) to obtain a new system. Then, another sliding mode σ2 can be designed for the new system based on σ1 and σ˙1.In order to satisfy the constraints of bounded state and control input and to eliminate the uncertainties caused by parameter perturbations and external disturbances, the nonlinear interactive sliding mode is designed as follows:(22)σ1=k1tanhk2e1+e˙1,σ2=k3tanhk4σ1+σ˙1,where const k1,k2,k3,k4>0. tanh⋅ is the hyperbolic tangent function, and it can limit the range of the variable. σ1 is a sliding mode on the phase plane e1e˙1, e˙1=−k1tanhk2e1 as σ1⟶0; when the deviation of e1 is big, e1 will convergence with a fixed rate(i.e., k1), and when the deviation of e1 is small, e1 will converge exponentially; k2 is a compression or tension coefficient, which determines the sliding slope in the origin of phase plane together with k1. σ2 is a sliding mode on the phase planeσ1σ˙1; when σ2⟶0, it can obtain σ˙1=−k3tanhk4σ1; when the deviation of σ1 is big, σ1 will converge with a fixed rate approximately(i.e., k3); if the deviation of σ1 is small, σ1 will converge exponentially; from the definition of σ1 and σ2, we can obtain that the stability control of system (15) is equal to the stability control of σ2.Taking the time derivative ofσ1, give(23)σ˙1=k14ek2e1+e−k2e12Τk2e˙1+e¨1=k14ek2e1+e−k2e12Τk2e˙1−M0−1C0e2−M0−1C0θ˙r−θ¨r+M0−1u+d+Y.Substituting (23) into (22) yields(24)σ2=k3tanhk4k1tanhk2e1+k4e2+k14ek2e1+e−k2e12Τk2e˙1−M0−1C0e2−M0−1C0θ˙r−θ¨r+M0−1u+d+Y.In order to avoid the estimation of the unknown moment of inertia and unknown external disturbance and simplify the implementation of attitude stability controller, the following incremental feedback controller is designed to ensure the stability of zero-order scalar system (24):(25)u˙=−kpσ2−εsatσ2,satx=xmin,ifx<xmin,x,ifxmin≤x≤xmax,xmax,ifx>xmax,where constkp,ε∈R+. xmin and xmax are the minimum value and maximum value of x, respectively. The incremental feedback controller does not need to estimate or cancel unknown moment of inertia and unknown external disturbance but only needs to search the solution of σ2=0 along with the opposite direction ofσ2, which avoids the waste of control input energy and reduces the magnitude of the control input torque effectively.


## 3.2. Analysis of the Stability

Barb˘alat Lemma: if a function f(t) is continuous and differentiable and satisfies the following:(1)
limt⟶∞ft=C C is bounded constant.(2)
f˙t is uniformly continuous (meansf¨t is bounded), then limt⟶∞f˙t=0.Theorem 1.
Under the iterative sliding mode (22) and incremental feedback controller (25), the closed-loop system (15) is globally asymptotic stable.Proof.
According to the above analysis, it is known that the stability control of system (15) is equal to the stability control of σ2.
Considering Lyapunov function asV=1/2σ2Τσ2 and taking the time derivative of V yields(26)V˙=σ2Τσ˙2=σ2Τ∂σ2∂uu˙.
Ignoring the variables of independent ofu in (24) and taking the time derivative of (24) about u, we obtain(27)∂σ2∂u=M0−1.
BecauseM0 is a positive matrix, M0−1 is also a positive matrix. Substituting (25) and (27) into (26) yields(28)V˙=σ2ΤM0−1−kpσ2−εsatσ2=−kpσ2ΤM0−1σ2−εσ2ΤM0−1satσ2≤−kpσ2ΤM0−1σ2≤0.
Then, we can obtainV˙=σ2Τσ˙2≤0, which means V˙≤0 when t≥0. According to the Barbalat lemma, we can obtain σ2⟶0 as t⟶∞, and from the definition of σ2and σ1, we can obtain σ1⟶0 and e1⟶0, respectively. Therefore, system (15) is globally asymptotic stable with the incremental feedback controller (25) and interactive sliding mode (22).


## 4. Adaptive Sliding Mode Attitude Stable Controller Design

For comparison, an adaptive sliding mode attitude stable controller is designed in this section. The adaptive law is used to estimate unknown moment of inertia online, and the sliding mode controller is used to eliminate external disturbance. Denoteξ=ΔJ11,ΔJ22,ΔJ33,ΔJ12,ΔJ13,ΔJ23T, and ξ^=ΔJ^11,ΔJ^22,ΔJ^33,ΔJ^12,ΔJ^13,ΔJ^23T is the estimation of ξ; define error of estimation as ξ˜=ξ^−ξ. Because −ΔMθθ¨−ΔCθ,θ˙θ˙−ΔGθ,θ˙ is linear with ξ, denote(29)Y=−ΔMθθ¨−ΔCθ,θ˙θ˙−G0θ,θ˙−ΔGθ,θ˙=−G0θ,θ˙+Ηξ,where His a matrix. Let sliding mode surface as(30)σ=kSL1e1+e2,where parameter kSL1 is a constant, then(31)σ˙=kSL1e˙1+e˙2=kSL1e2+−M0−1C0e2−M0−1C0θ˙r−θ¨r+M0−1u+d−G0θ,θ˙+Ηξ.Design the adaptive sliding mode controller (ASMC) as(32)u=−kSLpM0σ−Ηξ^−M0kSL1e2+C0e2+C0θ˙r+M0θ¨r+G0θ,θ˙−Dsatσ,and adaptive law as(33)ξ^˙=ΓΗΤσ,where kSL1 is a constant and Γ is a constant matrix.To verify the global asymptotic stability of closed-loop system, consider Lyapunov function(34)V=12σΤM0σ+12ξ˜ΤΓ−1ξ˜.Take the time derivative ofV and substituting (32) and (33) into derivative, we get(35)V˙=σΤM0σ˙+ξ˜ΤΓ−1ξ^˙=σΤM0kSL1e2+−M0−1C0e2−M0−1C0θ˙r−θ¨r+M0−1u+d−G0θ,θ˙+Ηξ+ξ˜ΤΓ−1ξ^˙=σΤM0kSL1e2+−M0−1C0e2−M0−1C0θ˙r−θ¨r+M0−1−kSLpM0σ−Ηξ^−M0kSL1e2+C0e2+C0θ˙r+M0θ¨r+G0−Dsatσ+d−G0+Ηξ+ξ˜ΤΓ−1ξ^˙=σΤM0M0−1−kSLpM0σ−Ηξ^−Dsatσ+d+Ηξ+ξ˜ΤΓ−1ξ^˙=−kSLpσΤM0σ−σΤDsatσ−d+ξ˜ΤΓ−1ξ^˙−ΗΤσ=−kSLpsΤM0σ−σΤDsatσ−d≤−kSLpσΤM0σ.Then,V˙≤−kSLpσΤM0σ≤0; therefore, the adaptive sliding mode controller (32) and (33) ensures the global asymptotic stability of the closed-loop system (15).


## 5. Numerical Simulation

### 5.1. Performance Comparison between the above Two Controllers

In this section, the iterative sliding mode and increment feedback controller(ISMIFC) is compared with the adaptive sliding mode controller (ASMC) in the environment of spacecraft on-orbit capturing noncooperation target. The initial attitude angle is θ0=10∘,10∘,10∘Τ, and correspondingly, angular velocity is θ˙0=0,0,0Τ, orbital rate is ω0=0.01rad/s, and nominal inertial matrix is J0=2000120901201700140901401500.The desired attitude angle isθr=1∘,1∘,1∘Τ, desired angle velocity is θ˙r0=0,0,0Τ, and acceleration is θ¨r0=0,0,0Τ. The external disturbance is d=sint,sint,sintΤ, and the perturbation matrix of moment of inertia is ΔJ=0.001×J0×t.When the control parameters are set ask1=1.2,k2=1,k3=0.008,k4=1,kp=1.5,ε=0.1, and the time of simulation is set as 800 s; the control based on incremental feedback controller is simulated. Under the iterative sliding mode and increment feedback controller, the attitude angle tracking error histories and control input torque are shown in Figures 1 and 2, respectively.Figure 1
Attitude angle tracking error with ISMIFCΔJ=0.001×J0×t.Figure 2
Control input torque with ISMIFCΔJ=0.001×J0×t.The adaptive sliding mode controller is simulated, control parameters are set askSLp=1.5,kSL1=1.2,D=1,Γ=diag3,3,3,3,3,3, and ξ^0=0,0,0,0,0,0T, and the time of simulation is set as 10 s. The attitude angle tracking error trajectory and control input torque of the adaptive sliding mode controller are shown in Figures 3 and 4, respectively.Figure 3
Attitude angle tracking error with ASMCΔJ=0.001×J0×t.Figure 4
Control input torque with ASMCΔJ=0.001×J0×t.Comparing with Figures1 and 3, it can see that the attitude angle of the two controllers decreases rapidly, but the convergence time of adaptive sliding mode controller is shorter than that of the iterative sliding and increment feedback controller, the convergence time of adaptive sliding mode controller is being about 7 s, and the iterative sliding controller is being about 700 s. In the process of attitude stability control, the control time is not a strict restrain, and the convergence time of 700 second is not long. Figures 2 and 4 show the control torque of two controllers, respectively. In Figure 2, the control torque is between −0.6 N·m and 0.6 N·m at the early stage and stabilizes at zero. In Figure 4, the control torque ranges from −550 N·m to 200 N·m and gradually declines to zero. Compared with the control torques of two controllers, the control input of iterative sliding controller consumes less than adaptive sliding controller, so the iterative sliding controller is more suitable in practical spacecraft attitude controller.


### 5.2. Robust Analysis of Iterative Sliding and Increment Feedback Controller

From Figures5 to 8, the trajectory of attitude, angle tracking error, and control input torque are given.Figure 5
Attitude angle tracking error with ISMIFCΔJ=0.01×J0×t.Figure 6
Control input torque with ISMIFCΔJ=0.01×J0×t.Figure 7
Attitude angle tracking error with ISMIFCΔJ=0.1×J0×t.Figure 8
Control input torque with ISMIFCΔJ=0.1×J0×t.In Figures5 and 6, the attitude angle error and control input torque are given by the perturbation matrix ΔJ=0.01×J0×t, respectively. When the perturbation matrix isΔJ=0.1×J0×t, the attitude angle error and control input torque are shown in Figures 7 and 8, respectively.Compared with Figures1, 5, and 7, it can be seen that the convergence time and the convergence curve are not the same. It can be seen from Figures 2, 6, and 8, the control input torque curve is similar.In conclusion, the convergence speed, convergence time, and control input torque of system are not affected by inertia disturbance. It is shown that the incremental feedback controller is robust to parameter perturbation and is suitable for spacecraft attitude stability control in the process of on-orbit capturing noncooperative target.


## 5.1. Performance Comparison between the above Two Controllers

In this section, the iterative sliding mode and increment feedback controller(ISMIFC) is compared with the adaptive sliding mode controller (ASMC) in the environment of spacecraft on-orbit capturing noncooperation target. The initial attitude angle is θ0=10∘,10∘,10∘Τ, and correspondingly, angular velocity is θ˙0=0,0,0Τ, orbital rate is ω0=0.01rad/s, and nominal inertial matrix is J0=2000120901201700140901401500.The desired attitude angle isθr=1∘,1∘,1∘Τ, desired angle velocity is θ˙r0=0,0,0Τ, and acceleration is θ¨r0=0,0,0Τ. The external disturbance is d=sint,sint,sintΤ, and the perturbation matrix of moment of inertia is ΔJ=0.001×J0×t.When the control parameters are set ask1=1.2,k2=1,k3=0.008,k4=1,kp=1.5,ε=0.1, and the time of simulation is set as 800 s; the control based on incremental feedback controller is simulated. Under the iterative sliding mode and increment feedback controller, the attitude angle tracking error histories and control input torque are shown in Figures 1 and 2, respectively.Figure 1
Attitude angle tracking error with ISMIFCΔJ=0.001×J0×t.Figure 2
Control input torque with ISMIFCΔJ=0.001×J0×t.The adaptive sliding mode controller is simulated, control parameters are set askSLp=1.5,kSL1=1.2,D=1,Γ=diag3,3,3,3,3,3, and ξ^0=0,0,0,0,0,0T, and the time of simulation is set as 10 s. The attitude angle tracking error trajectory and control input torque of the adaptive sliding mode controller are shown in Figures 3 and 4, respectively.Figure 3
Attitude angle tracking error with ASMCΔJ=0.001×J0×t.Figure 4
Control input torque with ASMCΔJ=0.001×J0×t.Comparing with Figures1 and 3, it can see that the attitude angle of the two controllers decreases rapidly, but the convergence time of adaptive sliding mode controller is shorter than that of the iterative sliding and increment feedback controller, the convergence time of adaptive sliding mode controller is being about 7 s, and the iterative sliding controller is being about 700 s. In the process of attitude stability control, the control time is not a strict restrain, and the convergence time of 700 second is not long. Figures 2 and 4 show the control torque of two controllers, respectively. In Figure 2, the control torque is between −0.6 N·m and 0.6 N·m at the early stage and stabilizes at zero. In Figure 4, the control torque ranges from −550 N·m to 200 N·m and gradually declines to zero. Compared with the control torques of two controllers, the control input of iterative sliding controller consumes less than adaptive sliding controller, so the iterative sliding controller is more suitable in practical spacecraft attitude controller.


## 5.2. Robust Analysis of Iterative Sliding and Increment Feedback Controller

From Figures5 to 8, the trajectory of attitude, angle tracking error, and control input torque are given.Figure 5
Attitude angle tracking error with ISMIFCΔJ=0.01×J0×t.Figure 6
Control input torque with ISMIFCΔJ=0.01×J0×t.Figure 7
Attitude angle tracking error with ISMIFCΔJ=0.1×J0×t.Figure 8
Control input torque with ISMIFCΔJ=0.1×J0×t.In Figures5 and 6, the attitude angle error and control input torque are given by the perturbation matrix ΔJ=0.01×J0×t, respectively. When the perturbation matrix isΔJ=0.1×J0×t, the attitude angle error and control input torque are shown in Figures 7 and 8, respectively.Compared with Figures1, 5, and 7, it can be seen that the convergence time and the convergence curve are not the same. It can be seen from Figures 2, 6, and 8, the control input torque curve is similar.In conclusion, the convergence speed, convergence time, and control input torque of system are not affected by inertia disturbance. It is shown that the incremental feedback controller is robust to parameter perturbation and is suitable for spacecraft attitude stability control in the process of on-orbit capturing noncooperative target.


## 6. Conclusion

The attitude stability control problem of spacecraft on-orbit capturing noncooperative target is analyzed, and a globally asymptotic stable nonlinear controller based on nonlinear iterative sliding mode and increment feedback control is proposed. The controller does not need to estimate the moment of inertia and external disturbances, and the control input torque is small due to the use of utilization of the bounded sliding mode surface. The simulation results show that the controller is robust for parameter variations with fast convergence speed and less control input torque.



---

*Source: 6622685-2021-01-25.xml*