# On-Line Compensation for the Disturbance Rejection Rate of a Platform Seeker Based on a High-Gain Extended State Observer

Authors: Li Wei; Liu Shixiang; Zhang Wenjie; Xia Qunli
Journal: International Journal of Aerospace Engineering
 (2019)

---

## Abstract

This paper focuses on the on-line compensation of the disturbance rejection rate (DRR) for a platform seeker. The mathematical model of the typical platform seeker based on the inertial space is established, and the line-of-sight (LOS) rate from different signal extraction points is strictly derived. Considering the spring torque disturbance and damping torque disturbance caused by the missile attitude motion, the seeker DRR transfer functions are deduced and the amplitude and phase characteristics at different frequencies are also analyzed. In order to close the engineering practice, the DRR parasitic loop (DRRPL) model of the seeker is rationally simplified and the stable region of the parasitic loop from different extraction points is also obtained. However, to increase the stability and guidance accuracy of the missile terminal flight, the compensation scheme based on the high-gain extended state observer (ESO) is used to estimate the disturbance torques and eliminate the seeker DRR effect. Numerical simulations are conducted to verify the effectiveness of the proposed scheme. The simulation results show that the seeker DRR effect mainly exists in the middle and low frequencies and the stable region of the parasitic loop at different signal extraction points is different. The proposed compensation scheme can effectively eliminate the parasitic loop effect of the seeker and increase the flight stability of the missile. It can reduce the terminal miss distance of the missile and improve the strike accuracy.


---

## Body

## 1. Introduction

Seekers are the tracking and detecting system of missiles, which are well used to acknowledge and track a target [1, 2]. The accurate line-of-sight information which is extracted and then provided to the guidance and control system of the missile is aimed at making the missile fly to the target according to the predetermined guidance law. Therefore, the stability and validity of the seeker detector directly determines the terminal strike accuracy of a missile [2, 3]. Usually, a seeker equipped with an inertial stable platform (ISP) is available to isolate the coupling effect from the attitude motion of the projectile and other external disturbances such as wind and air-stream-induced torque [4]. Actually, the seeker cannot completely eliminate the disturbance effect of body motion due to the existence of various disturbances such as wire pulling, friction, and mass imbalance, which produce various disturbance torques and reduce the tracking stability of the seeker [5, 6]. Instead, the body motion is coupled into the output LOS rate signals, which then induces the disturbance rejection rate characteristic of the platform seeker [6].The DRR represents the ability of the seeker to isolate the disturbance of the missile body. It not only affects the output accuracy of the LOS rate but also forms an extra parasitic loop in the guidance system, which affects the stability and guidance accuracy [7, 8]. At present, scholars at home and abroad pay less attention to the effect of the DRR parasitic loop on the guidance system. Cui et al. pointed out that the disturbance torque is the major factor of the DRR effect by establishing the platform seeker mathematical model [9]. Li et al. [10] analyzed the DRR model of a strapdown seeker, but the model is relatively simple, which is quite different from the DRR model of the platform seeker. Garnell [11] and Nesline and Zarchan [2] have studied the DRR induced by the radome boresight error slope and analyzed the influence on the guidance loop. Du et al. [12], Li et al. [13], and Song et al. [14] analyzed the influence of the positive feedback of DRR in a strapdown phased array seeker (SPAS) on the guidance loop stability and terminal guidance precision, in which the results of all these studies indicate that the positive feedback of DRR is a more severe parasitic loop.Obviously, the DRRPL effect of the platform seeker cannot be ignored for the missile’s terminal flight. Improving seeker isolation and reducing the DRR effect have also become an important research field. Two effective methods are proposed by researchers, in which one is to improve the machining accuracy and the other is to design improvements to the seeker control system [15, 16]. However, from the perspective of the control system, estimating disturbance torque and then compensating the output information of the seeker will be a feasible method. The decoupling of the LOS rate is realized using the feedforward estimation strategy in reference [17], and this method effectively suppressed the DRR effect. Xu et al. [18] and Sadhu and Ghoshal [3] used a low-pass filter to estimate the disturbance moment and then extracted the additional LOS rate caused by the disturbance. Liu et al. [16] estimated the disturbance torque by the Kalman filtering technique to compensate for the LOS rate in the seeker output; however, the derivation process has some mistakes and the obtained results do not look very intuitive either.As is well known, disturbance observer (DO) method has been introduced to compensate for strong nonlinearity, unmodeled dynamics, and external disturbances of uncertain systems [19–22]. Among these disturbance observers, the extended state observer, as the key module of active disturbance rejection control [23], is widely employed in various engineering control systems, such as the ESO-based finite-time tracking control for the servo system [24] and the ESO-based dynamic surface trajectory control for the quadrotor UAV [25]. Since the high-gain ESO can perform a simultaneous estimation of the state variables and the unknown disturbance [26], its interconnection in the seeker control system allows using the estimates instead of their true values to eliminate the DRR effect of the platform seeker and achieve the accurate LOS rate output.Furthermore, Wangze [27] pointed out that the stable loop of the platform seeker has a wide bandwidth and the LOS rate output can be extracted either from the stable loop command or from the angular rate gyro. Meanwhile, his research shows that the DRR of the seeker is closed related not only to the characteristics of the control system but also to the extraction point of guidance signal and the type of disturbance torque. Based on the above foundations, a novel on-line compensation method for the DRR effect of the platform seeker is proposed in this paper. The structure of this paper is as follows. In Section 2, the mathematical model of the platform seeker considering a different guidance signal extraction point is proposed and verified. The DRR transfer functions that come up with different disturbance torque are rigorously deduced, and the corresponding amplitude and phase characteristics are analyzed in Section 2.2.1. In addition, the dynamic model of DRR is simplified equivalently, and the stable region of DRRPL is also analyzed in Section 2.2.2. In Section 3, the unknown disturbance torque is estimated, while the LOS rate output is compensated based on the high-gain ESO technique. In Section 4, numerical simulations are conducted to verify the effectiveness of the proposed method. In Section 5, relevant conclusions to reduce the effect of DRRPL on the missile terminal guidance property are summarized.


## 2. Mathematical Model of the Platform Seeker and DRR

### 2.1. Extraction of the LOS Rate

The schematic diagram of the platform seeker is as depicted in Figure1, where oxbybzb is the body coordinate and xb stands for axial pointing of projectile body, xs is the pointing seeker detector, xt is the orientation of the target, ϑ is the attitude angle at the pitching direction of the projectile body, qs is the deflection angle between the detector and the inertial space, qt refers to the angle of sight of the missile target, and ε is the measuring error angle [15]. According to Figure 1(b), the deflection angle of stabilized platform can be expressed as
(1)φt=qst−ϑt.Figure 1
The structure of the platform seeker.
(a)
Wire disturbance
(b)
The angle relationship of the pitch planeThe schematic diagram of the wire disturbance is illustrated in Figure1(a), where the stabilized platform has relative motion to the missile body in the process of real-time target detection; thus, the disturbance of the spring torque caused by the pulling of wires will hinder the movement of the frame [28]. The disturbance of the spring torque can be expressed as
(2)Tst=Ksφt.Another major disturbance named as the damping torque is related to the rotation angular rate of the seeker platform, which can be expressed as(3)Tvt=Kwq̇st−ϑ̇t=Kwφ̇t.In Equations (2) and (3), φ̇ is the deflection angular rate of the stabilized platform and Ks and Kw are the spring torque coefficient and damping torque coefficient, respectively. Then, the motor dynamic model with the disturbance torque can be obtained as
(4)Jq̇s=TM−Kwφ̇−Ksφ,where J is the rotational inertia of the platform and TM is the output torque of servo motor. However, to simplify the seeker model, the disturbance of the spring torque and the damping torque can be uniformly represented by GDs, which expressed as
(5)GDs=Kw+Kss.To form a closed loop for the stable platform, the LOS rate in the inertia space feedback is acquired via the rate gyro. The schematic diagram of a typical platform seeker tracking control system is shown in Figure2 [29, 30], where G1s and G2s are the transfer functions of the correction network for the tracking loop and stable loop, respectively. k1 and k2 are the forward gains of the operational amplifier. In addition, L and R, respectively, indicate the inductance and capacitance of the servo machine winding, and KT is the moment coefficient of the servo motor mounted on a platform frame, and KE is the back electromotive force coefficient. GDs is the transfer function of the disturbance torque, while Hs is the transfer function of the angular rate gyro.Figure 2
Schematic diagram of the platform seeker tracking control system.As shown in Figure2, the tracking control system is used to achieve accurate tracking to the theoretical LOS angle, which mainly includes a stable loop, tracking loop, back-EMF loop, and disturbance torque loop. The disturbance torque loop couples the base motion to the seeker and reduces the precision. The back-EMF loop is generally small, and thus, it minimally affects seeker control precision. The LOS rate signals transferred to the missile guidance system can be obtained from the C point and the S point. Thus, the transfer functions of the platform seeker by different signal points can be expressed as
(6)Gq̇Cs=q̇sCsq̇ts=G1sk1Ls+RJs+GDsLs+R+G2sHsk2KTsJsLs+R+GDsLs+R+G2sHsk2KT+G1sk1G2sk2KT,(7)Gq̇Ss=q̇sSsq̇ts=G1sG2sHsk1k2KTsJsLs+R+GDsLs+R+G2sHsk2KT+G1sk1G2sk2KT.High-frequency dynamics, correction networks, and small values can be neglected in engineering [14–16]. Assuming G1s=1,G2s=1,R≈1, and L≈0, the equivalent stable loop gain is K2=k2KTG2s/J and the equivalent tracking loop gain is K1=k1G1s. Then Equations (6) and (7) can be simplified as
(8)Gq̇Cs=q̇sCsq̇ts=K1Js+GDs+K2Js2+GDs+K2s+K1K2,(9)Gq̇Ss=q̇sSsq̇ts=K1K2Js2+GDs+K2s+K1K2.In fact, for a perfect seeker, the influence of disturbanceGDs on the tracking characteristic is usually neglected when it is relatively small [15], and then, Equations (8) and (9) can both be transformed into
(10)Gq̇Cs≈Gq̇Ss≈q̇sq̇t=1s/K1+1,It is obvious from Equation (10) that the seeker dynamics of the LOS rate signal responses from the C point and the S point are consistent when neglecting the effects of zero deviations and detector noises on the control system of the seeker.The typical parameters of a platform seeker are shown as in Table1, and the step responses of the seeker tracking control system without considering any disturbances are shown in Figure 3(a). But for an imperfect seeker, the related zero deviations and detector noises are introduced, and for this purpose, it is assumed that the variances of the seeker detector noise (ΔDs) and rate gyro noise (ΔGs) are 0.0001°/s2 and 0.001°/s2, respectively. The zero deviations of the seeker detector (ΔD0) and servo system (ΔS0) are both chosen as 0.1°/s2. Then, the simulation results under various disturbances are shown in Figure 3(b)–3(f).Table 1
Typical parameters of a platform seeker.
Parameter
Value
Tracking loop gain
k1=12
Stable loop gain
k2=12
Inductance
L=0.0035H
Resistance
R=8Ω
Moment coefficient of servo motor
KT=0.2334N⋅m/A
Rotational inertia
J=0.0020995kg⋅m2
Back-EMF coefficient
KE=0.2344V/rad/AFigure 3
LOS rate extraction from the C point and S point.
(a)
Without any disturbance
(b)
Zero deviation of detector
(c)
Zero deviation of rate gyro
(d)
Zero deviation of servo system
(e)
Detector noise
(f)
Rate gyro noiseIt can be observed from Figure3(a) that the output of the LOS rate from the S point lags slightly behind the C point, while remaining consistent basically in Figure 3(b)–3(d). However, in Figure 3(b)–3(d), the steady-state gains of the output responses exist at a certain deviation compared with the expected LOS rate. In Figure 3(e), the influence of the seeker detector noise on the signal responses from the C point is larger than the S point because of the filtering effect of the seeker control system, while the opposite conclusion can be obtained in Figure 3(f). Finally, compared with Figures 3(e) and 3(f), it can be seen that the detector noise has a greater impact on the seeker control system than on the rate gyro noise.


### 2.2. Analysis of DRRPL

#### 2.2.1. The DRR Model of the Platform Seeker

The DRR represents the capability of the isolation base motion of a missile, which can guarantee a high-precision LOS rate. The smaller the DRR is, the stronger the ability of the isolation base disturbance the seeker has. Thus, the DRR transfer function in Figure2 can be expressed as
(11)Rs=Δq̇sϑ̇s,where ϑ̇s is the attitude angular rate of the missile and Δq̇s is the coupled LOS rate of the seeker induced by the missile motion. According to the different signal extraction methods shown in Figure 2, ignoring the influence of the back-EMF loop and angular rate gyro dynamic Hs, the DRR transfer functions of the seeker are obtained as
(12)RsC=Δq̇Csϑ̇s=−G1sGDsk1Ls+RsJsLs+R+GDsLs+R+G2sk2KT+G1sG2sk1k2KT,(13)RsS=Δq̇Ssϑ̇s=sGDsLs+RsJsLs+R+GDsLs+R+G2sk2KT+G1sG2sk1k2KT.The DRR amplitude ratio of the C point to the S point is−k1/s, which can be obtained by dividing Equation (12) by Equation (13) regardless of the correction link G1s. Thus, no matter what the disturbance torque is, the relation between DRR functions from the C point and the S point is inversely proportional to the operating frequency of the seeker and proportional to the forward gain k1 of the tracking loop, as shown in Figure 4.Figure 4
The DRR amplitude ratio of the C point to the S point.
(a)
Proportional coordinate
(b)
Logarithmic coordinateBy substituting Equations (2) and (3) into Equation (12), the DRR transfer functions from the C point can be obtained, as shown in
(14)RsvC=Δq̇vCsϑ̇s=−Kwk1Ls+RG1ssJsLs+R+KwLs+R+k2KTG2s+k1k2KTG1sG2s,(15)RssC=Δq̇sCsϑ̇s=−Ksk1G1sLs+Rs2JsLs+R+k2KTG2s+KsLs+R+k1k2KTG1sG2ss.Similarly, the DRR transfer functions from the S point can also be obtained as follows:(16)RsvS=Δq̇vSsϑ̇s=sKwLs+RsJsLs+R+KwLs+R+G2sk2KT+G1sG2sk1k2KT,(17)RssS=Δq̇sSsϑ̇s=KsLs+RsJsLs+R+k2KTG2s+k1k2KTG1sG2s+KsLs+R.From Equations (14)–(17), it is obvious that the DRR transfer functions of the seeker are different for different disturbance torques and guidance signal extraction points. Therefore, the frequency characteristic of the seeker DRR is closely related to the frequency characteristic of the seeker control system and cannot be simply regarded as a constant value of a certain frequency point.Based on Equations (14)–(17), and the typical parameters of a platform seeker which are given in Table 1, the DRR transfer function Bode diagrams of different disturbance torques and guidance signal extraction points are shown in Figure 5.Figure 5
Bode diagrams of the DRR characteristic.
(a)
S point under spring torque
(b)
S point under damping torque
(c)
C point under spring torque
(d)
C point under damping torqueAccording to Figure5, the amplitude frequency characteristic of the seeker DRR attenuates seriously in the high-frequency band, so it can be inferred that the influence of DRR on the performance of the seeker control system mainly exists in the middle and low frequency band. From the phase frequency characteristic, the initial phase of the DRR transfer function from the S point under the spring torque and damping torque is 0° and 90°, whereas the initial phase from the C point is −270° and −180°.


#### 2.2.2. Stable Region Analysis of DRRPL

In engineering, the seeker DRR model established in Section2.2.1 can be simplified under the following assumptions [15, 16, 31]:
(a)
High-frequency dynamics, correction links, and small values can be neglected in the seeker DRR analysis, i.e.,G1s=1,G2s=1,R≈1, and L≈0
(b)
The back-EMF loop is generally small, and thus, it minimally affects seeker control precision
(c)
Ignoring the rate gyro dynamics and the feedback gain of the stable loop is equal to 1 (i.e.,kg=1)
(d)
The forward gain of the stable loop is equivalent toK2=k2KTG2s/J, and the forward gain of the tracking loop is equivalent to K1=k1G1sBased on the above assumptions, the schematic diagram of platform seeker shown in Figure2 can be transformed into Figure 6, where q¨D is the disturbance angular acceleration induced by disturbance torques.Figure 6
Equivalent block diagram of the DRR.DefiningKs′=Ks/J and Kw′=Kw/J, Equations (14) and (15) can then be rewritten as Equations (18) and (19), respectively:
(18)RsvC=Δq̇vCsϑ̇s=−Kw′/K21/K1K2s2+Kw′+K2/K1K2s+1,(19)RssC=Δq̇sCsϑ̇s=−Ks′K1/Ks′+K1K2s1/Ks′+K1K2s2+K2/Ks′+K1K2s+1.Similarly, Equations (16) and (17) can be rewritten as Equations (20) and (21), respectively.
(20)RsvS=Δq̇vSsϑ̇s=Kw′s/K1K21/K1K2s2+Kw′+K2/K1K2s+1,(21)RssS=Δq̇sSsϑ̇s=Ks′/Ks′+K1K21/Ks′+K1K2s2+K2/Ks′+K1K2s+1.According to Table1, there exists K1K2≫1 and K1K2≫Ks′, and the DRR transfer functions shown in Equations (18)–(21) can be simplified to a first-order dynamic model, as shown in Table 2.Table 2
First-order DRR transfer function of the seeker.
C point
S point
Equivalent coefficient
Damping torque
RsvC=−Rw/Tds+1
RsvS=Rws/K1/Tds+1
Rw=Kw′/K2, Td=Kw′+K2/K1K2
Spring torque
RssC=−Rs/Tss+1s
RssS=Rs/K1/Tss+1
Rs≈Ks′/K2, Ts=K2/K1K2+Ks′In the terminal guidance phase, the seeker provides the guidance signal for the missile. However, if the seeker cannot isolate the body disturbance perfectly, the guidance signal not only includes the real LOS rate but also includes the additional signalΔq̇ which is caused by missile attitude motion. Then, the parasitic loop is stimulated, which is called the DRRPL. In order to analyze the characteristic of DRRPL intuitively, the simple structure of the DRR parasitic loop is shown in Figure 7. The guidance system is regarded as a five-order system with guidance time constant Tg [30], in which the guidance filter can be simplified as GFs=1/Tgs/5+1, whereas the autopilot is the third-order model Gas=1/Tgs/5+13. The seeker is regarded as the first-order model with Tseeker=Tg/5. Tα is the time constant of the angle of attack, which represents the rapidity of missile attitude response.
(22)Gs=amq̇t=NVcTgs/5+15−NVc/VmTαs+1Tgs/5+1Rs.Figure 7
The simple structure of the DRR parasitic loop.According to Figure7, the closed-loop transfer function of the DRR parasitic loop can be obtained, as shown in Equation (22), where Rs is the transfer function of the DRR model. Furthermore, the stable region of the DRRPL can be obtained, as shown in Figure 8.Figure 8
Stable region of the DRRPL.
(a)
Damping torque
(b)
Spring torqueIt can be observed from Figure8 that the stable region of the DRRPL decreases with the increase of the disturbance torque coefficients Rw and Rs. The stable region under the damping torque from the C point is wider than that from the S point, whereas the stable region under the spring torque from the C point is less than that from the S point. Thus, the stable region of DRRPL should be analyzed considering the main disturbance torques of the seeker caused by the trajectory characteristics, so as to select the appropriate extraction point of the LOS rate during the terminal-guided flight.


## 2.1. Extraction of the LOS Rate

The schematic diagram of the platform seeker is as depicted in Figure1, where oxbybzb is the body coordinate and xb stands for axial pointing of projectile body, xs is the pointing seeker detector, xt is the orientation of the target, ϑ is the attitude angle at the pitching direction of the projectile body, qs is the deflection angle between the detector and the inertial space, qt refers to the angle of sight of the missile target, and ε is the measuring error angle [15]. According to Figure 1(b), the deflection angle of stabilized platform can be expressed as
(1)φt=qst−ϑt.Figure 1
The structure of the platform seeker.
(a)
Wire disturbance
(b)
The angle relationship of the pitch planeThe schematic diagram of the wire disturbance is illustrated in Figure1(a), where the stabilized platform has relative motion to the missile body in the process of real-time target detection; thus, the disturbance of the spring torque caused by the pulling of wires will hinder the movement of the frame [28]. The disturbance of the spring torque can be expressed as
(2)Tst=Ksφt.Another major disturbance named as the damping torque is related to the rotation angular rate of the seeker platform, which can be expressed as(3)Tvt=Kwq̇st−ϑ̇t=Kwφ̇t.In Equations (2) and (3), φ̇ is the deflection angular rate of the stabilized platform and Ks and Kw are the spring torque coefficient and damping torque coefficient, respectively. Then, the motor dynamic model with the disturbance torque can be obtained as
(4)Jq̇s=TM−Kwφ̇−Ksφ,where J is the rotational inertia of the platform and TM is the output torque of servo motor. However, to simplify the seeker model, the disturbance of the spring torque and the damping torque can be uniformly represented by GDs, which expressed as
(5)GDs=Kw+Kss.To form a closed loop for the stable platform, the LOS rate in the inertia space feedback is acquired via the rate gyro. The schematic diagram of a typical platform seeker tracking control system is shown in Figure2 [29, 30], where G1s and G2s are the transfer functions of the correction network for the tracking loop and stable loop, respectively. k1 and k2 are the forward gains of the operational amplifier. In addition, L and R, respectively, indicate the inductance and capacitance of the servo machine winding, and KT is the moment coefficient of the servo motor mounted on a platform frame, and KE is the back electromotive force coefficient. GDs is the transfer function of the disturbance torque, while Hs is the transfer function of the angular rate gyro.Figure 2
Schematic diagram of the platform seeker tracking control system.As shown in Figure2, the tracking control system is used to achieve accurate tracking to the theoretical LOS angle, which mainly includes a stable loop, tracking loop, back-EMF loop, and disturbance torque loop. The disturbance torque loop couples the base motion to the seeker and reduces the precision. The back-EMF loop is generally small, and thus, it minimally affects seeker control precision. The LOS rate signals transferred to the missile guidance system can be obtained from the C point and the S point. Thus, the transfer functions of the platform seeker by different signal points can be expressed as
(6)Gq̇Cs=q̇sCsq̇ts=G1sk1Ls+RJs+GDsLs+R+G2sHsk2KTsJsLs+R+GDsLs+R+G2sHsk2KT+G1sk1G2sk2KT,(7)Gq̇Ss=q̇sSsq̇ts=G1sG2sHsk1k2KTsJsLs+R+GDsLs+R+G2sHsk2KT+G1sk1G2sk2KT.High-frequency dynamics, correction networks, and small values can be neglected in engineering [14–16]. Assuming G1s=1,G2s=1,R≈1, and L≈0, the equivalent stable loop gain is K2=k2KTG2s/J and the equivalent tracking loop gain is K1=k1G1s. Then Equations (6) and (7) can be simplified as
(8)Gq̇Cs=q̇sCsq̇ts=K1Js+GDs+K2Js2+GDs+K2s+K1K2,(9)Gq̇Ss=q̇sSsq̇ts=K1K2Js2+GDs+K2s+K1K2.In fact, for a perfect seeker, the influence of disturbanceGDs on the tracking characteristic is usually neglected when it is relatively small [15], and then, Equations (8) and (9) can both be transformed into
(10)Gq̇Cs≈Gq̇Ss≈q̇sq̇t=1s/K1+1,It is obvious from Equation (10) that the seeker dynamics of the LOS rate signal responses from the C point and the S point are consistent when neglecting the effects of zero deviations and detector noises on the control system of the seeker.The typical parameters of a platform seeker are shown as in Table1, and the step responses of the seeker tracking control system without considering any disturbances are shown in Figure 3(a). But for an imperfect seeker, the related zero deviations and detector noises are introduced, and for this purpose, it is assumed that the variances of the seeker detector noise (ΔDs) and rate gyro noise (ΔGs) are 0.0001°/s2 and 0.001°/s2, respectively. The zero deviations of the seeker detector (ΔD0) and servo system (ΔS0) are both chosen as 0.1°/s2. Then, the simulation results under various disturbances are shown in Figure 3(b)–3(f).Table 1
Typical parameters of a platform seeker.
Parameter
Value
Tracking loop gain
k1=12
Stable loop gain
k2=12
Inductance
L=0.0035H
Resistance
R=8Ω
Moment coefficient of servo motor
KT=0.2334N⋅m/A
Rotational inertia
J=0.0020995kg⋅m2
Back-EMF coefficient
KE=0.2344V/rad/AFigure 3
LOS rate extraction from the C point and S point.
(a)
Without any disturbance
(b)
Zero deviation of detector
(c)
Zero deviation of rate gyro
(d)
Zero deviation of servo system
(e)
Detector noise
(f)
Rate gyro noiseIt can be observed from Figure3(a) that the output of the LOS rate from the S point lags slightly behind the C point, while remaining consistent basically in Figure 3(b)–3(d). However, in Figure 3(b)–3(d), the steady-state gains of the output responses exist at a certain deviation compared with the expected LOS rate. In Figure 3(e), the influence of the seeker detector noise on the signal responses from the C point is larger than the S point because of the filtering effect of the seeker control system, while the opposite conclusion can be obtained in Figure 3(f). Finally, compared with Figures 3(e) and 3(f), it can be seen that the detector noise has a greater impact on the seeker control system than on the rate gyro noise.


## 2.2. Analysis of DRRPL

### 2.2.1. The DRR Model of the Platform Seeker

The DRR represents the capability of the isolation base motion of a missile, which can guarantee a high-precision LOS rate. The smaller the DRR is, the stronger the ability of the isolation base disturbance the seeker has. Thus, the DRR transfer function in Figure2 can be expressed as
(11)Rs=Δq̇sϑ̇s,where ϑ̇s is the attitude angular rate of the missile and Δq̇s is the coupled LOS rate of the seeker induced by the missile motion. According to the different signal extraction methods shown in Figure 2, ignoring the influence of the back-EMF loop and angular rate gyro dynamic Hs, the DRR transfer functions of the seeker are obtained as
(12)RsC=Δq̇Csϑ̇s=−G1sGDsk1Ls+RsJsLs+R+GDsLs+R+G2sk2KT+G1sG2sk1k2KT,(13)RsS=Δq̇Ssϑ̇s=sGDsLs+RsJsLs+R+GDsLs+R+G2sk2KT+G1sG2sk1k2KT.The DRR amplitude ratio of the C point to the S point is−k1/s, which can be obtained by dividing Equation (12) by Equation (13) regardless of the correction link G1s. Thus, no matter what the disturbance torque is, the relation between DRR functions from the C point and the S point is inversely proportional to the operating frequency of the seeker and proportional to the forward gain k1 of the tracking loop, as shown in Figure 4.Figure 4
The DRR amplitude ratio of the C point to the S point.
(a)
Proportional coordinate
(b)
Logarithmic coordinateBy substituting Equations (2) and (3) into Equation (12), the DRR transfer functions from the C point can be obtained, as shown in
(14)RsvC=Δq̇vCsϑ̇s=−Kwk1Ls+RG1ssJsLs+R+KwLs+R+k2KTG2s+k1k2KTG1sG2s,(15)RssC=Δq̇sCsϑ̇s=−Ksk1G1sLs+Rs2JsLs+R+k2KTG2s+KsLs+R+k1k2KTG1sG2ss.Similarly, the DRR transfer functions from the S point can also be obtained as follows:(16)RsvS=Δq̇vSsϑ̇s=sKwLs+RsJsLs+R+KwLs+R+G2sk2KT+G1sG2sk1k2KT,(17)RssS=Δq̇sSsϑ̇s=KsLs+RsJsLs+R+k2KTG2s+k1k2KTG1sG2s+KsLs+R.From Equations (14)–(17), it is obvious that the DRR transfer functions of the seeker are different for different disturbance torques and guidance signal extraction points. Therefore, the frequency characteristic of the seeker DRR is closely related to the frequency characteristic of the seeker control system and cannot be simply regarded as a constant value of a certain frequency point.Based on Equations (14)–(17), and the typical parameters of a platform seeker which are given in Table 1, the DRR transfer function Bode diagrams of different disturbance torques and guidance signal extraction points are shown in Figure 5.Figure 5
Bode diagrams of the DRR characteristic.
(a)
S point under spring torque
(b)
S point under damping torque
(c)
C point under spring torque
(d)
C point under damping torqueAccording to Figure5, the amplitude frequency characteristic of the seeker DRR attenuates seriously in the high-frequency band, so it can be inferred that the influence of DRR on the performance of the seeker control system mainly exists in the middle and low frequency band. From the phase frequency characteristic, the initial phase of the DRR transfer function from the S point under the spring torque and damping torque is 0° and 90°, whereas the initial phase from the C point is −270° and −180°.


### 2.2.2. Stable Region Analysis of DRRPL

In engineering, the seeker DRR model established in Section2.2.1 can be simplified under the following assumptions [15, 16, 31]:
(a)
High-frequency dynamics, correction links, and small values can be neglected in the seeker DRR analysis, i.e.,G1s=1,G2s=1,R≈1, and L≈0
(b)
The back-EMF loop is generally small, and thus, it minimally affects seeker control precision
(c)
Ignoring the rate gyro dynamics and the feedback gain of the stable loop is equal to 1 (i.e.,kg=1)
(d)
The forward gain of the stable loop is equivalent toK2=k2KTG2s/J, and the forward gain of the tracking loop is equivalent to K1=k1G1sBased on the above assumptions, the schematic diagram of platform seeker shown in Figure2 can be transformed into Figure 6, where q¨D is the disturbance angular acceleration induced by disturbance torques.Figure 6
Equivalent block diagram of the DRR.DefiningKs′=Ks/J and Kw′=Kw/J, Equations (14) and (15) can then be rewritten as Equations (18) and (19), respectively:
(18)RsvC=Δq̇vCsϑ̇s=−Kw′/K21/K1K2s2+Kw′+K2/K1K2s+1,(19)RssC=Δq̇sCsϑ̇s=−Ks′K1/Ks′+K1K2s1/Ks′+K1K2s2+K2/Ks′+K1K2s+1.Similarly, Equations (16) and (17) can be rewritten as Equations (20) and (21), respectively.
(20)RsvS=Δq̇vSsϑ̇s=Kw′s/K1K21/K1K2s2+Kw′+K2/K1K2s+1,(21)RssS=Δq̇sSsϑ̇s=Ks′/Ks′+K1K21/Ks′+K1K2s2+K2/Ks′+K1K2s+1.According to Table1, there exists K1K2≫1 and K1K2≫Ks′, and the DRR transfer functions shown in Equations (18)–(21) can be simplified to a first-order dynamic model, as shown in Table 2.Table 2
First-order DRR transfer function of the seeker.
C point
S point
Equivalent coefficient
Damping torque
RsvC=−Rw/Tds+1
RsvS=Rws/K1/Tds+1
Rw=Kw′/K2, Td=Kw′+K2/K1K2
Spring torque
RssC=−Rs/Tss+1s
RssS=Rs/K1/Tss+1
Rs≈Ks′/K2, Ts=K2/K1K2+Ks′In the terminal guidance phase, the seeker provides the guidance signal for the missile. However, if the seeker cannot isolate the body disturbance perfectly, the guidance signal not only includes the real LOS rate but also includes the additional signalΔq̇ which is caused by missile attitude motion. Then, the parasitic loop is stimulated, which is called the DRRPL. In order to analyze the characteristic of DRRPL intuitively, the simple structure of the DRR parasitic loop is shown in Figure 7. The guidance system is regarded as a five-order system with guidance time constant Tg [30], in which the guidance filter can be simplified as GFs=1/Tgs/5+1, whereas the autopilot is the third-order model Gas=1/Tgs/5+13. The seeker is regarded as the first-order model with Tseeker=Tg/5. Tα is the time constant of the angle of attack, which represents the rapidity of missile attitude response.
(22)Gs=amq̇t=NVcTgs/5+15−NVc/VmTαs+1Tgs/5+1Rs.Figure 7
The simple structure of the DRR parasitic loop.According to Figure7, the closed-loop transfer function of the DRR parasitic loop can be obtained, as shown in Equation (22), where Rs is the transfer function of the DRR model. Furthermore, the stable region of the DRRPL can be obtained, as shown in Figure 8.Figure 8
Stable region of the DRRPL.
(a)
Damping torque
(b)
Spring torqueIt can be observed from Figure8 that the stable region of the DRRPL decreases with the increase of the disturbance torque coefficients Rw and Rs. The stable region under the damping torque from the C point is wider than that from the S point, whereas the stable region under the spring torque from the C point is less than that from the S point. Thus, the stable region of DRRPL should be analyzed considering the main disturbance torques of the seeker caused by the trajectory characteristics, so as to select the appropriate extraction point of the LOS rate during the terminal-guided flight.


## 2.2.1. The DRR Model of the Platform Seeker

The DRR represents the capability of the isolation base motion of a missile, which can guarantee a high-precision LOS rate. The smaller the DRR is, the stronger the ability of the isolation base disturbance the seeker has. Thus, the DRR transfer function in Figure2 can be expressed as
(11)Rs=Δq̇sϑ̇s,where ϑ̇s is the attitude angular rate of the missile and Δq̇s is the coupled LOS rate of the seeker induced by the missile motion. According to the different signal extraction methods shown in Figure 2, ignoring the influence of the back-EMF loop and angular rate gyro dynamic Hs, the DRR transfer functions of the seeker are obtained as
(12)RsC=Δq̇Csϑ̇s=−G1sGDsk1Ls+RsJsLs+R+GDsLs+R+G2sk2KT+G1sG2sk1k2KT,(13)RsS=Δq̇Ssϑ̇s=sGDsLs+RsJsLs+R+GDsLs+R+G2sk2KT+G1sG2sk1k2KT.The DRR amplitude ratio of the C point to the S point is−k1/s, which can be obtained by dividing Equation (12) by Equation (13) regardless of the correction link G1s. Thus, no matter what the disturbance torque is, the relation between DRR functions from the C point and the S point is inversely proportional to the operating frequency of the seeker and proportional to the forward gain k1 of the tracking loop, as shown in Figure 4.Figure 4
The DRR amplitude ratio of the C point to the S point.
(a)
Proportional coordinate
(b)
Logarithmic coordinateBy substituting Equations (2) and (3) into Equation (12), the DRR transfer functions from the C point can be obtained, as shown in
(14)RsvC=Δq̇vCsϑ̇s=−Kwk1Ls+RG1ssJsLs+R+KwLs+R+k2KTG2s+k1k2KTG1sG2s,(15)RssC=Δq̇sCsϑ̇s=−Ksk1G1sLs+Rs2JsLs+R+k2KTG2s+KsLs+R+k1k2KTG1sG2ss.Similarly, the DRR transfer functions from the S point can also be obtained as follows:(16)RsvS=Δq̇vSsϑ̇s=sKwLs+RsJsLs+R+KwLs+R+G2sk2KT+G1sG2sk1k2KT,(17)RssS=Δq̇sSsϑ̇s=KsLs+RsJsLs+R+k2KTG2s+k1k2KTG1sG2s+KsLs+R.From Equations (14)–(17), it is obvious that the DRR transfer functions of the seeker are different for different disturbance torques and guidance signal extraction points. Therefore, the frequency characteristic of the seeker DRR is closely related to the frequency characteristic of the seeker control system and cannot be simply regarded as a constant value of a certain frequency point.Based on Equations (14)–(17), and the typical parameters of a platform seeker which are given in Table 1, the DRR transfer function Bode diagrams of different disturbance torques and guidance signal extraction points are shown in Figure 5.Figure 5
Bode diagrams of the DRR characteristic.
(a)
S point under spring torque
(b)
S point under damping torque
(c)
C point under spring torque
(d)
C point under damping torqueAccording to Figure5, the amplitude frequency characteristic of the seeker DRR attenuates seriously in the high-frequency band, so it can be inferred that the influence of DRR on the performance of the seeker control system mainly exists in the middle and low frequency band. From the phase frequency characteristic, the initial phase of the DRR transfer function from the S point under the spring torque and damping torque is 0° and 90°, whereas the initial phase from the C point is −270° and −180°.


## 2.2.2. Stable Region Analysis of DRRPL

In engineering, the seeker DRR model established in Section2.2.1 can be simplified under the following assumptions [15, 16, 31]:
(a)
High-frequency dynamics, correction links, and small values can be neglected in the seeker DRR analysis, i.e.,G1s=1,G2s=1,R≈1, and L≈0
(b)
The back-EMF loop is generally small, and thus, it minimally affects seeker control precision
(c)
Ignoring the rate gyro dynamics and the feedback gain of the stable loop is equal to 1 (i.e.,kg=1)
(d)
The forward gain of the stable loop is equivalent toK2=k2KTG2s/J, and the forward gain of the tracking loop is equivalent to K1=k1G1sBased on the above assumptions, the schematic diagram of platform seeker shown in Figure2 can be transformed into Figure 6, where q¨D is the disturbance angular acceleration induced by disturbance torques.Figure 6
Equivalent block diagram of the DRR.DefiningKs′=Ks/J and Kw′=Kw/J, Equations (14) and (15) can then be rewritten as Equations (18) and (19), respectively:
(18)RsvC=Δq̇vCsϑ̇s=−Kw′/K21/K1K2s2+Kw′+K2/K1K2s+1,(19)RssC=Δq̇sCsϑ̇s=−Ks′K1/Ks′+K1K2s1/Ks′+K1K2s2+K2/Ks′+K1K2s+1.Similarly, Equations (16) and (17) can be rewritten as Equations (20) and (21), respectively.
(20)RsvS=Δq̇vSsϑ̇s=Kw′s/K1K21/K1K2s2+Kw′+K2/K1K2s+1,(21)RssS=Δq̇sSsϑ̇s=Ks′/Ks′+K1K21/Ks′+K1K2s2+K2/Ks′+K1K2s+1.According to Table1, there exists K1K2≫1 and K1K2≫Ks′, and the DRR transfer functions shown in Equations (18)–(21) can be simplified to a first-order dynamic model, as shown in Table 2.Table 2
First-order DRR transfer function of the seeker.
C point
S point
Equivalent coefficient
Damping torque
RsvC=−Rw/Tds+1
RsvS=Rws/K1/Tds+1
Rw=Kw′/K2, Td=Kw′+K2/K1K2
Spring torque
RssC=−Rs/Tss+1s
RssS=Rs/K1/Tss+1
Rs≈Ks′/K2, Ts=K2/K1K2+Ks′In the terminal guidance phase, the seeker provides the guidance signal for the missile. However, if the seeker cannot isolate the body disturbance perfectly, the guidance signal not only includes the real LOS rate but also includes the additional signalΔq̇ which is caused by missile attitude motion. Then, the parasitic loop is stimulated, which is called the DRRPL. In order to analyze the characteristic of DRRPL intuitively, the simple structure of the DRR parasitic loop is shown in Figure 7. The guidance system is regarded as a five-order system with guidance time constant Tg [30], in which the guidance filter can be simplified as GFs=1/Tgs/5+1, whereas the autopilot is the third-order model Gas=1/Tgs/5+13. The seeker is regarded as the first-order model with Tseeker=Tg/5. Tα is the time constant of the angle of attack, which represents the rapidity of missile attitude response.
(22)Gs=amq̇t=NVcTgs/5+15−NVc/VmTαs+1Tgs/5+1Rs.Figure 7
The simple structure of the DRR parasitic loop.According to Figure7, the closed-loop transfer function of the DRR parasitic loop can be obtained, as shown in Equation (22), where Rs is the transfer function of the DRR model. Furthermore, the stable region of the DRRPL can be obtained, as shown in Figure 8.Figure 8
Stable region of the DRRPL.
(a)
Damping torque
(b)
Spring torqueIt can be observed from Figure8 that the stable region of the DRRPL decreases with the increase of the disturbance torque coefficients Rw and Rs. The stable region under the damping torque from the C point is wider than that from the S point, whereas the stable region under the spring torque from the C point is less than that from the S point. Thus, the stable region of DRRPL should be analyzed considering the main disturbance torques of the seeker caused by the trajectory characteristics, so as to select the appropriate extraction point of the LOS rate during the terminal-guided flight.


## 3. DRR Estimation and Compensation

Based on the analysis in Section2, it is known that the existence of the seeker DRR will inevitably affect the terminal strike accuracy of the missile, and the DRR effect is time-varying because the disturbance induced by the body attitude motion is uncertain. Secondly, the stable region of the DRRPL will decrease with the increase of the disturbance torque. The above conclusions indicate that it is necessary for the seeker to compensate for the DRR effect in real time. Therefore, in this section, the additional LOS rate caused by the disturbance torque is quickly estimated based on the ESO to compensate for the output of the seeker.


### 3.1. ESO-Based Estimation of Disturbance Torque

The specific model of the seeker disturbance torque is ignored, and the influence of the disturbance torque on the tracking control system of the seeker is replaced by the interference angular accelerationq¨D, as shown in Figure 9. Consistent with Section 2.2.2, the seeker stable loop gain and tracking loop gain are represented by K2 and K1, respectively. The high-gain ESO is used to estimate the angular acceleration, and then the additional LOS rate signals Δq̇c and Δq̇s are calculated according to the DRR model of the platform seeker. Finally, the compensated guidance signals q̇^c and q̇^s can be transmitted to the autopilot to achieve the flight control of the missile.Figure 9
Compensation scheme of DRR.Based on Figure9, the system state equations of the platform seeker can be established as
(23)ẋ1=x2ẋ2=K2u−q¨Dy=x1,where x1 and x2 are the state variables, which represent qs and q̇s, respectively. u is the control input of the servo motor. Defining an extended state, x3≜−q¨D for Equation (23), and the new extended order system can be written as
(24)ẋ1=x2ẋ2=x3+K2uẋ3=wty=x1,where wt is the rate of change of the disturbance angular acceleration, q¨D, and is assumed to be an unknown but bounded quantity. By making q¨D a state, however, it is now possible to estimate it by using the state estimator. To this end, consider a nonlinear observer of the form [32–34]
(25)e1=z1−yż1=z2−β01εe1α1signe1ż2=z3−β02ε2e1α2signe1+K2uż3=−β03ε3e1α3signe1,where β01>0,β03>0, and β02>β03/β01 and α3=α∈0,1,α2=2α3+1/3,α1=α3+2/3, and ε∈0,1. z1,z2 and z3 stand for the estimate values of qs,q̇s and −q¨D, respectively.Theorem 1.
Supposinge1=z1−y,e2=z2−ẏ, and e3=z3+q¨D stand for the estimation errors. Considering the extended system (24) and observer system (25), and that L>0,σ=1−γ/γ exists, γ∈0,1/2 makes the estimation errors stabilized in the following internals with finite time, as shown:
(26)z1−y≤Lε3σ,z2−ẏ≤Lε3σ−1,z3−x3≤Lε3σ−2Proof.
The error equations between the extended system (see Equation (24)) and the observer system (see Equation (25)) can be obtained as
(27)ė1=e2−β01e1α1signe1εė2=e3−β02e1α2ε2signe1ė3=−β03e1α3ε3signe1−wt.
Then, considering the following coordinate transformation(28)τ=tε,ηt=ε3wt,wiτ=εi−1eit,i=1,2,3.
Thus, Equation (27) can be rewritten as
(29)dw1dτ=w2−β01w1α1signw1dw2dτ=w3−β02w1α2signw1dw3dτ=−β03w1α3signw1−η.
Neglectingη, Equation (29) can be replaced by
(30)dw1dτ=w2−β01w1α1signw1dw2dτ=w3−β02w1α2signw1dw3dτ=−β03w1α3signw1.Definition 1 (see [35]).
Letfx=f1x,⋯,fnxT:Rn→Rn as a vector field. If ∀ε>0, then r1,⋯,rn∈Rn exists, where ri>0,i=1,2,⋯,n, and makes
(31)fiεr1x1,⋯,εrnxn=εk+rifix,i=1,2,⋯,n,where k>−minri,i=1,2,⋯,n, and then the vector field fx is homogeneous of degree k with dilation r=r1,⋯,rn.
According to Definition1, let fw stand for the vector function of system (29). When r1=1,r2=α+2/3,, and r3=2α+1/3, there exists
(32)fwcr1w1,cr2w2,εr3x3=cr2w2−β01cr1w1α1signw1cr3w3−β02cr1w1α2signw1−β03cr1w1α3signw1=cr1+m000cr2+m000cr3+mw2−β01w1α1signw1w3−β02w1α2signw1−β03w1α3signw1,where m=α−1/3<0, which is also homogeneous with vector function fw with dilation (r1,r2,r3).Lemma 1 (see [36]).
Suppose vector fieldfx is homogeneous with degree m with respect to ν. Then, the origin is a finite-time-stable equilibrium under fx if and only if the origin is an asymptotically stable equilibrium under fx and m<0.Lemma 2 (see [36]).
For any continuous functiongt and gt<δ0, and considering the following disturbed system
(33)ẋ=fx,t+gt.
If there is aC1 smooth positive definite function Vx, and for any β1>0 and β2∈0,0.5,which satisfies
(34)V̇x+β1Vβ2x≤0.
So there existsL0>0, γ0∈0,0.5, and a finite adjustment time Treach, and then the solution of system (33) satisfies
(35)x≤L0δ0σ0,where σ0=1−γ0/γ0.In addition, according to the Routh stability criterion,s3+k3s2+k2s+k1 is the Hurwitz polynomial only when k3>0,k1>0, and k2>k1/k3. The origin of system (30) is the equilibrium point of asymptotic stability because of the linear correlation between the vector field fw and polynomial s3+k3s2+k2s+k1. According to Lemma 1, the origin of system (30) is finite-time stable. Assuming that the upper limit of wt is w0, then ητ≤w0ε3. Next, by using Lemma 2, in which there is a positive number μ, this makes the solution of system (29) satisfy wτ≤μw0ε3σ in finite time. Finally, according to Equation (28), it can be observed that the estimation errors are able to meet ei≤wτε−i+1=Lε3σ−i+1, where L=μw0σ, thus completing the stability proof.Furthermore, to weaken the high-frequency chattering phenomenon in the simulation, the switching functionse1αsigne1 in Equation (25) can be replaced by
(36)fale1,αi,δ=e1δαi−1,e1≤δ,e1αisigne1,e1>δ,where δ is a small positive constant determining the interval length of the linear segment.


### 3.2. On-Line Compensation Scheme

As shown in Figure9, the angular acceleration induced by the different disturbance torque is estimated. Then, the additional LOS rate caused by the disturbance angular acceleration is also obtained on the basis of the DRR model of the seeker. Therefore, the LOS rate outputs of the seeker based on the C point and S point are compensated for by the additional LOS rate. The compensation scheme can be expressed as
(37)Gq¨⌢DCs=K1s2+K2kg/Js+K1K2/J,Gq¨⌢DSs=−kgss2+K2kg/Js+K1K2/J.Analysis of the above compensation process shows that the scheme does not need to pay more attention to the specific form of the disturbance torque model but depends on the accurate dynamic model of the platform seeker. The above compensation scheme is different from that proposed in reference [16], where the linear filtering algorithm is adopted and errors exist in its deduction process.


## 3.1. ESO-Based Estimation of Disturbance Torque

The specific model of the seeker disturbance torque is ignored, and the influence of the disturbance torque on the tracking control system of the seeker is replaced by the interference angular accelerationq¨D, as shown in Figure 9. Consistent with Section 2.2.2, the seeker stable loop gain and tracking loop gain are represented by K2 and K1, respectively. The high-gain ESO is used to estimate the angular acceleration, and then the additional LOS rate signals Δq̇c and Δq̇s are calculated according to the DRR model of the platform seeker. Finally, the compensated guidance signals q̇^c and q̇^s can be transmitted to the autopilot to achieve the flight control of the missile.Figure 9
Compensation scheme of DRR.Based on Figure9, the system state equations of the platform seeker can be established as
(23)ẋ1=x2ẋ2=K2u−q¨Dy=x1,where x1 and x2 are the state variables, which represent qs and q̇s, respectively. u is the control input of the servo motor. Defining an extended state, x3≜−q¨D for Equation (23), and the new extended order system can be written as
(24)ẋ1=x2ẋ2=x3+K2uẋ3=wty=x1,where wt is the rate of change of the disturbance angular acceleration, q¨D, and is assumed to be an unknown but bounded quantity. By making q¨D a state, however, it is now possible to estimate it by using the state estimator. To this end, consider a nonlinear observer of the form [32–34]
(25)e1=z1−yż1=z2−β01εe1α1signe1ż2=z3−β02ε2e1α2signe1+K2uż3=−β03ε3e1α3signe1,where β01>0,β03>0, and β02>β03/β01 and α3=α∈0,1,α2=2α3+1/3,α1=α3+2/3, and ε∈0,1. z1,z2 and z3 stand for the estimate values of qs,q̇s and −q¨D, respectively.Theorem 1.
Supposinge1=z1−y,e2=z2−ẏ, and e3=z3+q¨D stand for the estimation errors. Considering the extended system (24) and observer system (25), and that L>0,σ=1−γ/γ exists, γ∈0,1/2 makes the estimation errors stabilized in the following internals with finite time, as shown:
(26)z1−y≤Lε3σ,z2−ẏ≤Lε3σ−1,z3−x3≤Lε3σ−2Proof.
The error equations between the extended system (see Equation (24)) and the observer system (see Equation (25)) can be obtained as
(27)ė1=e2−β01e1α1signe1εė2=e3−β02e1α2ε2signe1ė3=−β03e1α3ε3signe1−wt.
Then, considering the following coordinate transformation(28)τ=tε,ηt=ε3wt,wiτ=εi−1eit,i=1,2,3.
Thus, Equation (27) can be rewritten as
(29)dw1dτ=w2−β01w1α1signw1dw2dτ=w3−β02w1α2signw1dw3dτ=−β03w1α3signw1−η.
Neglectingη, Equation (29) can be replaced by
(30)dw1dτ=w2−β01w1α1signw1dw2dτ=w3−β02w1α2signw1dw3dτ=−β03w1α3signw1.Definition 1 (see [35]).
Letfx=f1x,⋯,fnxT:Rn→Rn as a vector field. If ∀ε>0, then r1,⋯,rn∈Rn exists, where ri>0,i=1,2,⋯,n, and makes
(31)fiεr1x1,⋯,εrnxn=εk+rifix,i=1,2,⋯,n,where k>−minri,i=1,2,⋯,n, and then the vector field fx is homogeneous of degree k with dilation r=r1,⋯,rn.
According to Definition1, let fw stand for the vector function of system (29). When r1=1,r2=α+2/3,, and r3=2α+1/3, there exists
(32)fwcr1w1,cr2w2,εr3x3=cr2w2−β01cr1w1α1signw1cr3w3−β02cr1w1α2signw1−β03cr1w1α3signw1=cr1+m000cr2+m000cr3+mw2−β01w1α1signw1w3−β02w1α2signw1−β03w1α3signw1,where m=α−1/3<0, which is also homogeneous with vector function fw with dilation (r1,r2,r3).Lemma 1 (see [36]).
Suppose vector fieldfx is homogeneous with degree m with respect to ν. Then, the origin is a finite-time-stable equilibrium under fx if and only if the origin is an asymptotically stable equilibrium under fx and m<0.Lemma 2 (see [36]).
For any continuous functiongt and gt<δ0, and considering the following disturbed system
(33)ẋ=fx,t+gt.
If there is aC1 smooth positive definite function Vx, and for any β1>0 and β2∈0,0.5,which satisfies
(34)V̇x+β1Vβ2x≤0.
So there existsL0>0, γ0∈0,0.5, and a finite adjustment time Treach, and then the solution of system (33) satisfies
(35)x≤L0δ0σ0,where σ0=1−γ0/γ0.In addition, according to the Routh stability criterion,s3+k3s2+k2s+k1 is the Hurwitz polynomial only when k3>0,k1>0, and k2>k1/k3. The origin of system (30) is the equilibrium point of asymptotic stability because of the linear correlation between the vector field fw and polynomial s3+k3s2+k2s+k1. According to Lemma 1, the origin of system (30) is finite-time stable. Assuming that the upper limit of wt is w0, then ητ≤w0ε3. Next, by using Lemma 2, in which there is a positive number μ, this makes the solution of system (29) satisfy wτ≤μw0ε3σ in finite time. Finally, according to Equation (28), it can be observed that the estimation errors are able to meet ei≤wτε−i+1=Lε3σ−i+1, where L=μw0σ, thus completing the stability proof.Furthermore, to weaken the high-frequency chattering phenomenon in the simulation, the switching functionse1αsigne1 in Equation (25) can be replaced by
(36)fale1,αi,δ=e1δαi−1,e1≤δ,e1αisigne1,e1>δ,where δ is a small positive constant determining the interval length of the linear segment.


## 3.2. On-Line Compensation Scheme

As shown in Figure9, the angular acceleration induced by the different disturbance torque is estimated. Then, the additional LOS rate caused by the disturbance angular acceleration is also obtained on the basis of the DRR model of the seeker. Therefore, the LOS rate outputs of the seeker based on the C point and S point are compensated for by the additional LOS rate. The compensation scheme can be expressed as
(37)Gq¨⌢DCs=K1s2+K2kg/Js+K1K2/J,Gq¨⌢DSs=−kgss2+K2kg/Js+K1K2/J.Analysis of the above compensation process shows that the scheme does not need to pay more attention to the specific form of the disturbance torque model but depends on the accurate dynamic model of the platform seeker. The above compensation scheme is different from that proposed in reference [16], where the linear filtering algorithm is adopted and errors exist in its deduction process.


## 4. Simulations and Results

In this section, numerical simulations are conducted to validate the effectiveness of the proposed compensation scheme. The DRR compensation effect is analyzed under different disturbance models. Secondly, the proposed DRR compensation scheme is applied to the missile guidance flight process, and its impact on miss distance is also analyzed.


### 4.1. Analysis of DRR Compensation Effect

The desired LOS rate which isq̇t=0.5∘/s and the disturbance angular acceleration are selected as the sinusoidal signal and colored noise signal, respectively, as shown in Figure 10. The proposed compensation scheme is executed to eliminate the effect of disturbance, and the observer gains are selected as β01=10,β02=500,β03=4680,α=0.25, and ε=0.8. The LOS rate outputs from the C point and the S point while considering the DRR compensation are shown in Figures 11 and 12.Figure 10
Disturbance angular acceleration.
(a)
Sinusoidal signal
(b)
Color noise signalFigure 11
LOS rate outputs under sinusoidal disturbance.
(a)
C point output with sinusoidal disturbance
(b)
S point output with sinusoidal disturbanceFigure 12
LOS rate outputs under colored noise.
(a)
C point output with color noise
(b)
S point output with color noiseIn Figures11(a) and 11(b), the estimated effects of the C point and S point are compared under sinusoidal disturbance angular acceleration. Findings show that the LOS rate output from the C point is very smooth with less noise than the S point due to the effect of angular rate gyro noise. Moreover, the seeker DRR shows a sinusoidal trend and can be effectively suppressed by adopting the compensation scheme. However, the basic characteristics are consistent with that of the sinusoidal disturbance when the colored noise signal is used as the disturbance angular acceleration, as shown in Figure 12. Therefore, whatever angular acceleration is used as the disturbance signal, the proposed DRR compensation scheme can be used to estimate the disturbance moment signal and enhance the extraction accuracy of the LOS rate.


### 4.2. Influence of DRR Compensation on Miss Distance

In order to analyze the influence of the seeker DRR compensation scheme on the terminal miss distance of a missile, the block diagram of the missile guidance and control system is established, as shown in Figure13. The decoupled LOS rate is transmitted to the proportional navigation law after first-order filtering and then generating the required overload command signal. The guidance command is used as the input of a third-order autopilot model to generate the available overload for achieving the flight control of the missile. In addition, in Figure 13, errors mainly include target maneuvering and initial pointing error, and the initial pointing error is selected as the source of the whole simulation error.Figure 13
Block diagram of the missile guidance and control system.The typical flight parameters of a missile are listed in Table3, and the initial velocity pointing error angle and trajectory simulation time are chosen as ε=10∘ and tf=5s, respectively. Considering the influence of the seeker DRRPL caused by the damping moment and spring moment, the related simulation results are obtained and compared based on the proposed DRR compensation scheme, as shown in Figures 14–16.Table 3
Typical flight parameters of missile.
Parameter
Value
Parameter
Value
Navigation ratio
N=4
Missile velocity
Vm=450m/s
Closing velocity
Vc=300m/s
Rate gyro gain
kg=1.0
Autopilot time constant
Ta=0.1s
Filtering time constant
TF=0.1s
Guidance time constant
Tg=0.5s
Angle-of-attack time constant
Tα=1.0sFigure 14
Damping torque.
(a)
Disturbance torque
(b)
LOS rateThe LOS rate outputs from the C point and S point are given in Figure14 under damping torque. In Figure 14(a), the disturbance torque curves become smoother after being compensated, which indicates that the effect on the guidance and control system is less. However, in Figure 14(b), it can be observed that the LOS rate outputs using the compensation scheme keep consistent with no DRR situation. In all, the simulation results show that the proposed DRR compensation scheme can weaken the effect of the damping torque and improve the stability of the missile guidance and control system significantly. The same conclusions can be obtained when the spring torque is the disturbance source, as shown in Figure 15.Figure 15
Spring torque.
(a)
Disturbance torque
(b)
LOS rateThe missile flight trajectory curves while considering the influences of the initial velocity pointing error and seeker DRR parasitical loop are shown in Figure16. When the seeker disturbance torque is the damping torque, the results show that the trajectory curve from the S point for guidance signal extraction appears with a certain diverge, and the terminal miss distance is relatively large. However for the C point, the trajectory curve appears to have a relatively large deviation and the miss distance is large. After adopting the DRR compensation scheme, the trajectory curves of the guidance signal extracted from the C point and the S point are both smoother and the miss distances become smaller, which are similar to the situation without considering the DRR effect. In Figure 16(b), the same conclusions can be obtained when considering the spring torque disturbance.Figure 16
Trajectory curves.
(a)
Damping torque
(b)
Spring torqueIn conclusion, all the simulation results above confirm that the seeker DRR compensation scheme based on the high-gain ESO can effectively suppress the negative impact of the DRRPL on the missile terminal flight. The scheme can be compared with reference [16], which proposed a linear filtering algorithm to estimate the damping torque and spring torque and then compensate the LOS rate outputs. Errors about the deduction process exist in reference [16], and related simulation results are given (see the Appendix). Comparing the simulation results of these two papers, the effect of the seeker DRR on both the missile guidance and the control system can be suppressed effectively, but our method has better robustness and is simpler to implement.


## 4.1. Analysis of DRR Compensation Effect

The desired LOS rate which isq̇t=0.5∘/s and the disturbance angular acceleration are selected as the sinusoidal signal and colored noise signal, respectively, as shown in Figure 10. The proposed compensation scheme is executed to eliminate the effect of disturbance, and the observer gains are selected as β01=10,β02=500,β03=4680,α=0.25, and ε=0.8. The LOS rate outputs from the C point and the S point while considering the DRR compensation are shown in Figures 11 and 12.Figure 10
Disturbance angular acceleration.
(a)
Sinusoidal signal
(b)
Color noise signalFigure 11
LOS rate outputs under sinusoidal disturbance.
(a)
C point output with sinusoidal disturbance
(b)
S point output with sinusoidal disturbanceFigure 12
LOS rate outputs under colored noise.
(a)
C point output with color noise
(b)
S point output with color noiseIn Figures11(a) and 11(b), the estimated effects of the C point and S point are compared under sinusoidal disturbance angular acceleration. Findings show that the LOS rate output from the C point is very smooth with less noise than the S point due to the effect of angular rate gyro noise. Moreover, the seeker DRR shows a sinusoidal trend and can be effectively suppressed by adopting the compensation scheme. However, the basic characteristics are consistent with that of the sinusoidal disturbance when the colored noise signal is used as the disturbance angular acceleration, as shown in Figure 12. Therefore, whatever angular acceleration is used as the disturbance signal, the proposed DRR compensation scheme can be used to estimate the disturbance moment signal and enhance the extraction accuracy of the LOS rate.


## 4.2. Influence of DRR Compensation on Miss Distance

In order to analyze the influence of the seeker DRR compensation scheme on the terminal miss distance of a missile, the block diagram of the missile guidance and control system is established, as shown in Figure13. The decoupled LOS rate is transmitted to the proportional navigation law after first-order filtering and then generating the required overload command signal. The guidance command is used as the input of a third-order autopilot model to generate the available overload for achieving the flight control of the missile. In addition, in Figure 13, errors mainly include target maneuvering and initial pointing error, and the initial pointing error is selected as the source of the whole simulation error.Figure 13
Block diagram of the missile guidance and control system.The typical flight parameters of a missile are listed in Table3, and the initial velocity pointing error angle and trajectory simulation time are chosen as ε=10∘ and tf=5s, respectively. Considering the influence of the seeker DRRPL caused by the damping moment and spring moment, the related simulation results are obtained and compared based on the proposed DRR compensation scheme, as shown in Figures 14–16.Table 3
Typical flight parameters of missile.
Parameter
Value
Parameter
Value
Navigation ratio
N=4
Missile velocity
Vm=450m/s
Closing velocity
Vc=300m/s
Rate gyro gain
kg=1.0
Autopilot time constant
Ta=0.1s
Filtering time constant
TF=0.1s
Guidance time constant
Tg=0.5s
Angle-of-attack time constant
Tα=1.0sFigure 14
Damping torque.
(a)
Disturbance torque
(b)
LOS rateThe LOS rate outputs from the C point and S point are given in Figure14 under damping torque. In Figure 14(a), the disturbance torque curves become smoother after being compensated, which indicates that the effect on the guidance and control system is less. However, in Figure 14(b), it can be observed that the LOS rate outputs using the compensation scheme keep consistent with no DRR situation. In all, the simulation results show that the proposed DRR compensation scheme can weaken the effect of the damping torque and improve the stability of the missile guidance and control system significantly. The same conclusions can be obtained when the spring torque is the disturbance source, as shown in Figure 15.Figure 15
Spring torque.
(a)
Disturbance torque
(b)
LOS rateThe missile flight trajectory curves while considering the influences of the initial velocity pointing error and seeker DRR parasitical loop are shown in Figure16. When the seeker disturbance torque is the damping torque, the results show that the trajectory curve from the S point for guidance signal extraction appears with a certain diverge, and the terminal miss distance is relatively large. However for the C point, the trajectory curve appears to have a relatively large deviation and the miss distance is large. After adopting the DRR compensation scheme, the trajectory curves of the guidance signal extracted from the C point and the S point are both smoother and the miss distances become smaller, which are similar to the situation without considering the DRR effect. In Figure 16(b), the same conclusions can be obtained when considering the spring torque disturbance.Figure 16
Trajectory curves.
(a)
Damping torque
(b)
Spring torqueIn conclusion, all the simulation results above confirm that the seeker DRR compensation scheme based on the high-gain ESO can effectively suppress the negative impact of the DRRPL on the missile terminal flight. The scheme can be compared with reference [16], which proposed a linear filtering algorithm to estimate the damping torque and spring torque and then compensate the LOS rate outputs. Errors about the deduction process exist in reference [16], and related simulation results are given (see the Appendix). Comparing the simulation results of these two papers, the effect of the seeker DRR on both the missile guidance and the control system can be suppressed effectively, but our method has better robustness and is simpler to implement.


## 5. Conclusions

In this paper, the mathematical model of the platform seeker is built and the feasibility of the seeker structure concept is verified, which have laid the theoretical basis for the following research. The DRR model of the seeker is obtained considering different disturbance torques, and the amplitude and phase characteristics at different frequencies are also analyzed. Based on the reasonable assumptions, the DRRPL of the seeker is simplified and the stable region from different signal extraction points is obtained. The guidance and control system model of the missile is established considering the parasitic loop coupling, and a novel compensation scheme is used to eliminate the DRR effect. Through numerical simulation validation, the following conclusions can be drawn:(1)
The output of the LOS rate from the S point lags slightly behind the C point. The outputs of the seeker are at a certain offset relative to the desired LOS rate under the influence of the detector zero deviation, while they chatter when the disturbance is detector noise.
(2)
The stable region under the damping torque from the C point is wider than that from the S point, whereas the stable region under the spring torque from the C point is less than that from the S point
(3)
The DRR effect of the seeker mainly manifests in the low and middle frequencies and can be well suppressed by the proposed compensation scheme. The compensation scheme can effectively improve missile guidance accuracy and reduce terminal miss distanceOur future work will focus on a more realistic seeker model and disturbance torque model and then verify the effectiveness of the proposed scheme through mathematical simulation and hardware-in-the-loop simulation.



---

*Source: 3106732-2019-09-08.xml*