Abstract

This paper presents a composite control algorithm for the radio telescope servo control system to address the target tracking and matched/mismatched disturbance suppression problems. The algorithm consists of the equivalent-input-disturbance (EID) approach and the optimal control method. An EID estimator is developed using the difference between the estimated output of the state observer and the measurement output and then the estimate of the EID is fed forward into the control input to reject the disturbance. A cost function with clear physical meaning is selected and the weighting parameters are adjusted for the optimal controller to improve tracking performance. Considering the nonminimum phase characteristics of the radio telescope system, the state observer gain is computed using the linear matrix inequality (LMI) method. The system stability is analyzed using the small gain theorem, and the linear quadratic regulator (LQR) control method is utilized to determine the state feedback gain. Finally, the composite controller is designed for an identified telescope model. Simulation results show that for the tracking performance, the settling time of proposed method is 1.13 s and reduces by about 0.32 s and 0.87 s than that of the ADRC controller and PI controller, respectively. For the antidisturbance performance, the RMS value and the maximum error of the proposed method are 0.0039 and 0.0128, which are 42.86% and 40.38% of the ADRC controller and 30.71% and 27.77% of the PI controller, respectively, which indicates that the proposed method has better control performance. In addition, the proposed controller has certain robustness to systematic parameter perturbations.

1. Introduction

Radio telescopes play a crucial role in collecting radio waves emitted by celestial objects [1]. Sensitivity and spatial angular resolution are important indicators for measuring the performance of a radio telescope. Increasing the aperture of the reflecting surface is the main way to improve those indicators of the telescope [2], such as the 110 m Qitai radio telescope (QTT) aims to allow high-sensitivity observations from 150 MHz up to 115 GHz [3]. However, this leads to a decrease in the mechanical resonant frequency of the antenna [4]. In addition, the transfer function from the motor’s output torque to the displacement at the end of the flexible structure always displays the characteristics of the nonminimum phase [5], which makes the design of the antenna’s servo controllers challenging.

The primary function of the antenna is to locate and continuously track targets accurately. As a result, the servo system of the antenna should possess high tracking accuracy, excellent transient performance, and strong resistance to disturbances. However, radio telescopes are typically located in outdoors and are susceptible to external disturbances such as wind gusts and temperature fluctuations [6]. These factors can increase the difficulty of designing antenna servo system controllers.

Currently, the main control algorithm used in the antenna servo control system is the proportional-integral-derivative (PID) controller. While the PID controller is easy to adjust and has a simple structure, it is insufficient for meeting the high-frequency observation requirements of large aperture telescopes. With the development of the control theory, linear quadratic Gaussian (LQG) control [7, 8] and control [9] have been successfully applied to improve the tracking and antidisturbance performance of the antenna. However, LQG control requires an accurate mathematical model and lacks robustness [10], while control sacrifices control performance for increased robustness [11]. In addition, these control algorithms are single degree of freedom and there is a tradeoff between tracking performance and disturbance suppression [12].

The active disturbance compensation method has two degrees of freedom, one degree is used to estimate and compensate for the disturbances and the other is designed to achieve the expected control performance. The sliding mode controller and its improvements, combined with some observers, have been applied to flexible joint robots to improve their tracking performance and antidisturbance performance [1316]. However, the design of the controllers is very complex. Active disturbance rejection control (ADRC) has been successfully applied to improve the control performance of the telescope servo system [17, 18]. However, standard ADRC is limited to integral chain systems [19], and the extended state observer (ESO) can only tackle slowly varying disturbances [20]. A disturbance observer-based composite position controller (DOB-CPC) is designed to enhance antidisturbance capability of Leighton Chajnantor telescope [21]. A cascade acceleration feedback control (AFC) enhanced by disturbance observation and compensation is proposed to improve the tracking precision of telescope systems [22]. However, the DOB mainly estimates input disturbances and requires the inverse model of the system, which may destroy the system’s stability when the plant is a nonminimum phase system [23].

Another active disturbance rejection method is the equivalent input disturbance (EID) approach [24]. The fundamental concept is to define an EID in the control input channel which has the same impact on the system output as external disturbances, estimate the EID using an EID estimator, and then introduce the estimated EID into the input channel to compensate for the external disturbance. The EID approach does not require the inverse model of the system and can compensate for many kinds of unmeasurable disturbances, no matter if it is matched or mismatched. It has been successfully applied to time-delay systems [25], nonlinear systems [26], and linear uncertain systems [27]. What is more, some improved methods were developed to enhance the antidisturbance performance of the EID approach. For example, dynamic compensator and high-order low-pass filter were introduced to improve the disturbance rejection performance [28, 29]. Unlike the disturbance observer and the feedback controller are designed simultaneously to guarantee the closed system robust stability [30, 31], the EID-based control method has an important feature; the state-feedback control law can be designed independently of the observer gain and the disturbance estimator when stability is the only concern [24]. Also, the EID approach can compensate for the uncertainties and disturbances, so the designed controller has certain robustness [32].

Applying the EID method to the antenna servo control system can improve the antenna’s antidisturbance ability under unknown disturbances. However, the antenna servo system is a nonminimum phase system, and the perfect regulation method used in [28, 29] cannot be applied to obtain the observer gain. In addition, the method presented in [33] is too complex and needs solving three inequalities.

Based on these considerations, this paper proposes a new parameter design method for the antenna servo controller based on the EID approach. To resolve the nonminimum phase characteristic of the antenna system model, the LMI method based on the bounded real lemma is proposed to obtain the state observer gain. To simplify the problem of difficulty in selecting weight parameters in optimal control, the weighting parameters of a cost function with clear physical meaning are adjusted to get an optimal controller using the LQR control method. The EID estimator is also designed for disturbance rejection, and the system’s stability is analyzed using the small gain theorem.

The primary contributions of this paper are as follows:(1)Utilizing a composite control algorithm with the EID approach and the optimal control method can enhance the tracking performance, suppressing matched and mismatched disturbances of the antenna servo system.(2)Unlike the methods used in [25, 28, 29], this paper proposes a cost function with clear physical meanings and achieves satisfactory tracking performance using the LQR control method by adjusting the weighting parameters.(3)The proposed method of using the LMI technique to get observer gain is very simple and can apply to minimum-phase and nonminimum-phase systems.

The remaining sections of the paper are organized as follows. Section 2 outlines the antenna model and controller structure. In Section 3, the stability of the closed system, based on the separation theorem, is analyzed, and a method for designing the controller parameters is proposed. Section 4 presents simulation results on an antenna servo system to demonstrate the effectiveness of the proposed method. Finally, Section 5 provides conclusions.

2. Antenna Model and Controller Structure

2.1. Antenna Model Analysis

The servo control system of a radio telescope is a typical electromechanical system that consists of driving motors, a reducer, a main reflector, and a feedarm supporting subreflector. In the service mode, the reflector is driven by the azimuth axis and elevation axis motors to ensure the accuracy of pointing and tracking. Figure 1 shows the structure of a fully steerable radio telescope.

As the aperture increases, the structure of the telescope servo system becomes more flexible and complex, making it difficult to model using the traditional mechanism modeling method. We can obtain a close approximation model of the actual antenna by employing system identification methods. Since the azimuth and elevation axes of the antenna are orthogonal, which results in minimal coupling between the two systems [17, 34], this allows us to treat them as independent subsystems. In the paper, we focus on identifying the azimuth axis position open-loop system and designing a position loop controller.

Figure 2 shows the structure of the azimuth axis system used for identification experiments. The antenna model includes the drive model and the antenna structure. The speed loop corrected by the speed controller has good dynamic characteristics, which can be described using a first-order transfer function, namely,where is the time constant of the speed loop. The antenna structure can be described by a transfer function where and , , , are to be determined parameters. Then, the whole antenna model can be described as

Equation (3) can be described by , where , and are the state matrix, input matrix, and output matrix. Through the system identification experiment, the parameters , and can be obtained. Then, we can get the following state space equation to describe the dynamic model of the radio telescope.where , , , , , and . , , and are, respectively, the state vectors, the control input, and the measurement output.

2.2. Structure of the EID-Based Controller

Since the antenna system is affected by different types of disturbances and most of the disturbances do not enter the system through the input channel, the system model and the disturbances can be described by`where , . If , the disturbance is a matched one; otherwise, it is mismatched.

Assumption 1. The state of system (5) is fully controllable and observable, and the controlled object has no zeros at the origin.

Remark 2. The controllability ensures that the designed system is stable, observability is necessary for designing the system state observer, and having no zeros at the origin ensures that the system can accurately track a step signal without any steady-state error.
The system block diagram of the radio telescope servo system, based on the EID approach, is presented in Figure 3. The designed controller consists of internal model control, state feedback gain, Luenberger state observer, and an EID estimator.
As seen in the literature [24], there is an EID, , in the input channel which generates the same output as the unknown mismatched disturbance for the system. Consequently, system (5) can be converted intoA Luenberger state observer is constructed to estimate the ,where is the estimate of the system state, is the estimated output, and is the state observer gain. Defining the state observation error asand incorporating it in (5) yieldsAssuming there exists , letand by combining equations (7)–(10), it can be concluded thatwhere is the Moore–Penrose generalized inverse of , namely,
Since the EID estimate may be influenced by measured noise and there is a causal relationship between and the control input , a first-order low-pass filter,is designed to solve these problems. And, should satisfywhere is the maximum angular frequency of the external disturbance. The cutoff frequency of the low-pass filter is generally chosen as . And then,where are the Laplace transforms of .
For the tracking performance of the radio telescope, based on the reference signal of the system, the internal model is obtained aswhere and are the reference input and the measurement output, are the constant matrices, and is the state. The internal model is used to guarantee zero-error tracking performance in the steady state.
By designing the state feedback of the system, the control law of the whole system is finally obtained aswherewhere is the state feedback gain of the system.

3. Stability Analysis and Controller Parameters Design

3.1. Stability Analysis of Closed-Loop System

Stability is the inherent property of linear systems and is independent of external signals. Therefore, let and , then the system (5) is transformed to (4). Based on the above equation and equations (7), (8), (11), and (16), we get

Then, the transfer function from to is

The above equation suggests that the state feedback gain is independent of the observer gain. According to the separation theorem, the state feedback gain and the observer gain can be independently designed. As shown in Figure 4 which is equivalently converted from Figure 3, the control system is divided into two subsystems; one with the state feedback gain and internal model and the other with an EID estimator and the .

For subsystem 1, according to the small gain theorem, if the system is stable, it should satisfy

Based on the above analysis, we can get the following theorem:

Theorem 3. System (5) satisfies Assumption 1, and the designed system satisfies the following conditions [35]:(a) is a Hurwitz matrix;(b) is a Hurwitz matrix;(c);then, the whole control system designed by the EID approach and the optimal control method is asymptotically stable.

Proof. Since the EID-based system can be decomposed into two subsystems (subsystem 1 and subsystem 2), the stability of the whole system is equivalent to that of these two subsystems. For the subsystem 1, it is composed of and . First, and should be stable, and the condition (2) must be held. According to the small gain theorem, the subsystem 1 is stable if condition (3) holds. For the subsystem 2, it is a state feedback control system, when the state feedback gain has been designed, and the subsystem 2 is stable if condition (1) holds.

3.2. Controller Parameters Design

The out-loop controller of the system is composed of the internal model and optimal state feedback, which is designed through the LQR control method. The generalized controlled plant is described aswhere

Design the control law, , to find the control input that minimizes a cost function. Select a target output aswhere is the position output of the system; is the speed of the output that means how fast the system can reach the inference input; is the integration of the residual between the reference input signal and the output signal, which reflects the tracking performance of the system; and and are the weight coefficient. The cost function of the system is designed as

Adjust , , and , and use the LQR control method to obtain state feedback gains to meet the tracking performance.

The servo control system of the radio telescope has nonminimum phase characteristics. However, the perfect regulation method [36] is only applicable to the minimum phase system. Therefore, in this paper, we utilize the LMI method based on the bounded real lemma to calculate the observer gain. This method is suitable for both minimum-phase and nonminimum-phase systems.

Lemma 4 (Bounded real lemma [37]). For a continuous system, a minimum realization iswhere , , , and are the system matrices. If there is a positive definite symmetric matrix that satisfiesthen the system is stable and . Note: .

For the dual system of (5),where is a virtual external signal of the dual system. By designing a control law and integrating it into equation (25), the transfer function from the input to the output can be obtained as

By comparing equations (19) and (28), we find that the amplitude change of has impact on the amplitude of , thereby allowing us to obtain the state observer gain that satisfies the stability conditions in Theorem 3. The following theorem outlines how to derive the gain of the state observer using the bounded real lemma.

Theorem 5. For a continuous system, the minimum realization iswhere , , , , , and . If there exists a symmetric positive definite matrix and a general matrix such thatholds, then there exists , making that the system is stable and , where is the state feedback gain and is the transfer function from to .

Proof. Based on the bounded real lemma, there exists a symmetric positive definite matrix such thatholds, then the system is stable and . According to the performance criterion for linear matrix inequalities, multiplying both sides of an inequality by a positive definite matrix does not change the negative definite of the inequality. Multiplying the inequality (31) by on the left and on the right, respectively, yields the equivalent matrix inequality of equation (31) as follows:Let and , then the matrix inequality (30) is obtained.

3.3. The Controller Design Procedure

Based on the above analysis, we can achieve the following design procedure:(1)Select a cost function, adjust the weighting parameters , and get the state feedback gain , which satisfies the stability and tracking performance;(2)Select of the external disturbance;(3)Choose a low-pass filter that satisfies(4)Adjust and use the method in Theorem 5 to get a suitable , which satisfies the stability conditions (2) and (3) in Theorem 3.

4. Simulation Analysis

The Nanshan 26 m radio telescope (NSRT) is used to validate the effectiveness of the proposed approach. First, the identification experiment is carried out to obtain the model of NSRT. Then, we apply the proposed controller to the antenna servo system and analyze the performances of tracking, disturbance rejection, and robustness.

4.1. The Model of NSRT

To acquire the dynamic model of the NSRT, we use pseudorandom binary sequence (PRBS) as the input and the position encoder signal as the output. By analyzing the input data and the differential of the output data using the correlation analysis method and the eigensystem realization algorithm [38], we obtain the system parameter model.

The Bode diagrams of the data model and identification model are present in Figure 5. It shows that the frequency characteristics of the identified system model and the data model are consistent. The Bode plot reveals that the radio telescope demonstrates the desired rigidity characteristics at frequencies below 2 Hz. However, at frequencies above 2 Hz, the system exhibits noticeable flexibility, including resonance and antiresonance peaks. By analyzing the system model, it is discovered that the system has right half plane zeros, which further complicates the controller design. The antenna position loop model is formed by adding an integrator to the identification model (34). The state space equation of the position open-loop system is

4.2. Controller Design and Performance Analysis

Accurate pointing is a major function of radio telescopes, and now we mainly investigate its dynamic response characteristics when tracking a step signal and the antidisturbance performance after steady state. Based on the internal model principle, we get

To obtain the state feedback gain satisfying the tracking performance, the optimal tracking controller for the generalized controlled plant is designed by adjusting , , and . By trial and error, we obtain

By using the LQR control method, the state feedback gain is obtained.

Wind load is the main disturbance to the radio telescope. The wind spectrum exhibits significant energy values within the frequency range of 0-1 Hz. In order to assess the effectiveness of the EID estimator in mitigating disturbances, a composite disturbance consists of a step signal and various sinusoidal signals at different frequencies, as depicted in Figure 6, is introduced into the system, namely,

The maximum angular frequency of the disturbance is , then the first-order low-pass filter is selected as

The state observer gain is designed based on the dual system of (4). We adjust the parameter to achieve a state observer gain that satisfies the stability requirements using the method described in Theorem 5. The stability is test by examining the bode plots of and . Figure 7 displays the Bode plots of and obtained using different values of the parameter . When , we getand , which meets the stability condition in Theorem 3.

4.3. Performance Analysis

Assuming that the NSRT servo system is affected by matched and unmatched disturbances, the results of the simulation are shown in Figures 8 and 9, respectively.

At time t = 0 s, a step signal is applied to the system. Figures 8 and 9 show that regardless of whether the system includes the EID estimator, the tracking performance is the same, the system can track the reference input with almost no overshoot, and the settling time is 1.13 s. This indicates that the designed system has high tracking accuracy and excellent transient performance. After the system reaches a steady state, the disturbance is introduced to the system during t = 5–10 s. Since there is no model of the disturbance in the controller, the disturbance cannot be rejected to zero. For the matched disturbance, the maximum errors of the output caused by the disturbances are 0.0075 with the EID estimator and 0.0228 without the EID estimator. For the mismatched disturbance, they are 0.0113 and 0.0341, respectively. The system with the EID estimator exhibits better performance in rejecting disturbances compared to the system without the EID estimator.

Figure 10 shows the estimated disturbances of matched and unmatched disturbances as well as the actual disturbance. For the estimated disturbances of matched disturbance, there is a slight difference between them since there is in equation (10). For the estimated disturbances of mismatched disturbance, there is a large difference between them since the estimated disturbance is the estimation of the EID of mismatched disturbance.

To simplify the analysis, we focus on the worst-case scenario, in which the system is affected by a mismatched disturbance.

In order to demonstrate the advantage of the proposed method, the traditional PI controller and the ADRC controller are designed for the system. The method described in [4] is used to obtain the PI gains.

Since the relative order of the antenna system is two, the ADRC consists of a three-order LESO and a PD feedback control law. Based on the method in [], the parameters are , , and , and all poles of the LESO are chosen to be 24, which yielded the observer gain as follows:Also, a low-pass filter is added to the control input to guarantee the control performance, which is designed as

The results in Figure 11 demonstrate the system performance with three different control methods. Also, the tracking performance indices are presented in Table 1. Figure 11 and Table 1 show that when the system is controlled by the PI controller, the overshoot is about 23% and the settling time is about 2 s. And, for the ADRC controller, the overshoot is about 1% and the settling time is about 1.45 s. The above analysis indicates the proposed method has better tracking performance compared to the ADRC controller and the PI controller.

To evaluate the antidisturbance performance of the three methods, two performance indices are used, namely, the root mean square criterion (RMS) and the maximum error. The two indices are also presented in Table 1. We can get that the RMS values for the three methods are 0.0039, 0.0091, and 0.0127, while the maximum errors are 0.0128, 0.0317, and 0.0461, respectively. The RMS and maximum error value of the proposed method are 42.86% and 40.38% of the ADRC controller and 30.71% and 27.77% of the PI controller, respectively. This indicates when the EID-based controller applied to the NSRT servo system, there is a significant improvement in antidisturbance capability.

The antenna servo system is always subjected to parameter uncertainty. To verify the robustness of the proposed method, it is assumed that there are parameter uncertainties in the matrices , and the maximum perturbation is . By randomly selecting 5 uncertainties to verify the robustness, Figure 12 shows that the tracking and disturbance rejection performance of the NSRT servo system is not affected by the parameter perturbation, which indicates that the EID-based controller has certain robustness.

5. Conclusion

This paper introduces a composite control method to improve the tracking and disturbance suppression performance of the NSRT servo control system. The weighting parameters of the cost function with physical significance are investigated and the state feedback gain of the servo system is constructed using the LQR control method. For the nonminimum phase characteristics of the radio telescope system, the observer gain is obtained by using the LMI method, and the EID estimator is incorporated into the state feedback system to estimate and compensate for the disturbance. The simulation results show that the tracking performance of the proposed method is better than the ADRC controller and the traditional PI controller. Compared with the ADRC controller and the traditional PI controller, the RMS and the maximum disturbance-induced error in the output of the composite control method are smallest. These results indicate that the EID-based optimal control offers superior control performance.

Plans for future research are mainly in the following areas:(1)Study the uncertainty range of the antenna servo system and apply it to the design process of the controller to further improve the robustness of the whole system.(2)Combine the EID approach and other advanced control methods, such as MPC and SMC, to further improve the disturbance rejection performance and tracking accuracy of the radio telescope servo system.

Data Availability

The data used to support the findings of this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the National Key Research and Development Program of China under grant numbers 2021YFC2203603 and 2021YFC2203501. The research work was also partly supported by the Operation, Maintenance, and Upgrading Fund for Astronomical Telescopes and Facility Instruments, budgeted by the Ministry of Finance of China (MOF), and administrated by the Chinese Academy of Sciences (CAS).