Design of Composite Feedback and Feedforward Control Law for Aircraft Inertially Stabilized Platforms
The design of the control systems of the inertially stabilized platforms (ISPs) as part of airborne equipment for the majority of aircraft has its peculiarity. The presence of rate gyros in the inertial measurement unit gives the possibility to measure the rotation rate of the ISP base, which is the main disturbance interfering with the ISP accuracy. Inclusion of the feedforward disturbance gain in the control law with the simplest PI feedback significantly improves the accuracy of stabilization by the invariance theory. A combination of feedback and feedforward controllers produces a synergetic effect, thus, improving ISP accuracy. This article deals with the design of the airborne ISP control systems consisting of two stages: the parametric optimization of the PI feedback control based on composite “performance-robustness” criterion and the augmentation of the obtained system with feedforward gain. To prove the efficiency of the proposed control laws, the simulation of the ISP was undertaken. We have used a simulation of the heading-hold system of the commuter aircraft Beaver and the yaw rate output of this closed-loop system we have used as a source of the disturbance. The results of modeling proved the efficiency of the proposed design method.
Nowadays, the inertially stabilized platforms (ISPs) are the most powerful means of stabilizing high-precision optical and optical-electronic equipment on moving vehicles. They may be used as the onboard equipment of various vehicles: aerial, terrestrial, and marine [1, 2]. Their applications for piloted aircraft and unmanned aerial vehicles (UAVs) are significantly important, as they could be the main part of commercial payload  as well as an important part of navigation equipment, especially for UAV landing , stabilization, and trajectory tracking . In these cases, UAVs are unable to perform their flight missions without inertially stabilized visual equipment. The application of ISPs in aerial vehicles has its peculiarities that distinguish it from terrestrial and marine applications.
Firstly, the ISPs in aerial vehicles must be lighter, smaller, cheaper, and more power-efficient than their terrestrial and marine counterparts. This imposes further restrictions on their design parameters, particularly, on the computational power of the on-board computers. These restrictions are absent in the terrestrial and marine ISPs . Taking this fact into account, ISP control laws must be as simple as possible to implement them in the airborne computers.
Secondly, practically every aircraft is equipped with the Strapdown Inertial Navigation System (SINS), which is the kernel of the aircraft navigation system . This fact gives the possibility of direct measurement of the aircraft’s body rotation rate, which is the main source of the ISP disturbance. It is possible to design the ISP control laws based on feedback and feedforward (or combined) structure as far as the disturbance is measurable . Such an approach was intensively developed in the former Soviet literature as “invariance theory” . In American literature, it is known as the “disturbance rejection” . The simplest verbal description of this approach is given in . Some applications of this approach in other areas are given in [11, 12]. As it is known [1, 2], the entire ISP consists of two main contours: high-bandwidth rate inner loop for image stabilization and a lower-bandwidth outer loop for some target pointing or tracking. The first one is frequently used as an independent system for manual target tracking, and this article deals with the control law design for such a system. Aside from this introduction, the structure of the article consists of the second item, where the problem statement and the ISP linearized mathematical model are considered; the third item deals with the design of feedback and feedforward control law based on robustness and invariance theories; the fourth item deals with the simulation of the designed control system in the Simulink toolbox; and the fifth item is the conclusion.
2. Problem Statement and Mathematical Model of ISP
The ISP consists of the biaxial gimbals (for the control of the azimuth and elevation angles) driven by two electric motors. We will consider the design of the control system for the azimuth axis only because the elevation angle control is similar. Figure 1 represents the block diagram of ISP combined control system, consisting of a feedback controller “Fb.con.” and feedforward controller “Ff.con.”.
Other elements in Figure 1 stand for as follows: “Dist.mod.” is the disturbance model, “ISP model” represents a model of the DC motor with the azimuth frame of gimbals, are the ISP output signal and the ISP command signal in the tracking mode, respectively, is the disturbance signal, and is the pulse width of the PWM transducer, which is the control input for ISP. Here and further, we will assume the linear model of the pulse-width modulator. So, in this case, there are two channels of disturbance propagation. The first one is through the disturbance input of the controlled plant, and the second one is through the control input. For the sake of brevity, we will concisely describe the disturbance model, which is the series connection of two systems. The first one is the turbulence Dryden model [13, 14] of the lateral turbulent wind, and the second one is the model of the closed-loop heading angle stabilization system for commuter aircraft Beaver. We have used only the yaw rate output of this system as the source of the disturbance for ISP. As far as, in our case, this disturbance can be directly measured, it is possible to obtain the rejection of the measured disturbance by feedforward controller, and suppressing the fluctuations of the output variable by the feedback controller. Under [8–10], the feedback controller is used for the suppression of all possible measurable and nonmeasurable disturbances, and the feedforward controller is used for main disturbance rejection.
Now, we can formulate the problem statement for the feedback controller as follows. Let the dynamics of the controlled plant be described by the standard state-space model . where are the state propagation, control input, disturbance input, and observation matrices, respectively. Other variables in (1) are control input , disturbance , and output . The state vector consists of components , where is the torque of the DC electric motor. The most relevant approach to solve the disturbance suppression problem for the system (1) is the minimization of the bounded [16, 17]. In the general case, this gain is defined as follows [16, 17]: where are corresponding weighting matrices. The paper  describes the procedure of the minimization for static output feedback only. It was shown in  that for the more general case of static and dynamic output feedback, it is possible to use parametric optimization procedure for minimization of the and simultaneously -norm of the transfer function between the disturbance and ISP output where is the possible minimal value , thus guaranteeing certain robustness of system [16, 17]. Controllers, based on results [16, 17], must have the structure of the static output feedback and also the observer-based dynamic output feedback . In our case, it is possible to use simple controllers’ proportional-integral type recommended in [1, 2], which can be designed by a less sophisticated parametric optimization procedure. For achieving this goal, it is necessary to define the ISP mathematical model using concrete numerical parameters. Analysis of its dynamic properties gives the possibility to make demands on the optimization procedure.
For considered ISP complied with initial specifications, we have chosen the DC servomotor EDU-02  with reducer ratio . Total moment of inertia of servomotor loaded with gimbals equals . Other parameters have the following values: time constant of electric motor , viscous friction coefficient . Coulomb friction moment equals , where . The speed control of the servomotor uses a PWM transducer with maximal pulse width 2.5 msec. Suppose that the model of PWM transducer is linear, then the transfer function of servomotor from this transducer to the motor torque will have the following form where Taking into account these facts, it is possible to write down the following equations of dynamics [1, 2, 6]. where is the viscous friction moment, and . Here, is the rotation rate of the motor with the motionless base, and is the vehicle rotation rate (the yaw rate in our case). The last one is the main disturbance of this control system. We suppose that high-frequency vibrations of the vehicle body are suppressed by the vibration isolation means . One of the difficulties of the control system design is the nonlinear term in the first equation of system (1) due to Coulomb friction . As far as the disturbance and the rotation rate of gimbals are the random processes, it is necessary to use the random-input describing function [19, 20] for the statistical linearization of the nonlinear term . As it is known [19, 20], the coefficients of linearization depend on the ratio , where and are the mean value and the root mean square (rms) of the input random process , respectively. Due to the ISP’s main goal, the control system must maintain the rotation rate of the gimbals frame near zero in the steady-state mode; therefore, it can be possible to accept as the first approximation that . In this case, it will be enough to use only an equivalent statistical gain of nonlinear elements respective to the root mean squares , where and are the root mean squares of output and input signals of this element. Taking into account this circumstance, the equivalent statistical gain of a nonlinear element equals .
The root mean square depends on the type of aircraft used as the ISP carrier and its flight conditions. As it was mentioned above, we have chosen the small commuter aircraft Beaver as a carrier. All aerodynamic characteristics, which are necessary to create its mathematical model, are circumstantially presented in  along with the Dryden model of atmospheric turbulence (see also ). The yaw rate of this carrier in the automatic heading-hold flight mode in the stochastically disturbing atmosphere was accepted as the model of disturbance . The Dryden filter was created for the following flight conditions: altitude , true airspeed , and rms of instant wind velocity . And its output was applied to the correspondent input of the mathematical model of the closed-loop heading angle stabilization system. The root mean square of the yaw rate, in this case, obtained as a result of system modeling was . Substituting this value in (6), we can estimate . Using as the state vector, as the control input, and as the disturbance, we can rewrite system (5) in the standard Cauchy form (1), or
Substituting numerical values of coefficients in (7), we obtain model of the open loop system
Note that two complex eigenvalues of the state propagation matrix in (8) have very large eigenfrequency: 779 rad/sec (124 Hz). This fact must be taken into account in the procedure of control law design. As it was stated in [1, 2], the most popular control law for such controlled plants is the simplest PI controller. The transient processes in the closed-loop systems with such types of controlled plants have very high-frequency and low-frequency modes; therefore, such systems from a mathematical point of view can be classified as “stiff” systems. Differentiating the transient processes in real time requires a very high computational burden. Therefore, PID controllers are not well suited for these situations. Taking into account these considerations, we can specify the initial general problem statement for ISP control law design as follows. (1)Design of ISP feedback control law for the disturbance suppression based on the bounded approach. As far as the type of the PI feedback was determined from the very beginning [1, 2], it is necessary to design a parametric optimization procedure for this simplest control law.(2)Design of the ISP feedforward control law for the closed-loop system, determined at the previous step. This procedure based on the invariance theory will be described below.(3)Simulation of the designed system in order to evaluate the efficiency of the proposed control law.
2.1. Design of ISP Control Law
As it was stated before, we have chosen the PI control law for ISP feedback loop as
Therefore, the goal of the optimization procedure is the determination of the proportional and integral gains as a result of running some optimization procedure.
From item 1 of the problem statement, it is obvious that, first of all, we have to find some optimal gains in the control law (9) without feedforward, i.e., the gain of “Ff.con” (see Figure 1) must be equal to zero at this stage. Following aforementioned the minimization approach and expressions (2) and (3), we proposed cost function for optimization procedure including and with weighting coefficients allowing to achieve the trade-off between robustness and performance of control system. Also, sometimes, it is necessary to impose the limitations on the values of the control gains in order to avoid the saturation of the amplifiers of digital to analog converters (like PWM, for example) . Eventually, the composite performance-robustness index based on the above-mentioned considerations looks like where are weighting coefficients, and is the argument of the frequency response characteristic. We have to use a weighting transfer function in order to estimate norms in the realistic frequency bandwidth (no more than 50 Hz), so this transfer function is chosen as the second-order filter with eigenfrequency 314 rad/sec and unit damping ratio. Weighting coefficients were accepted as . These values were obtained by trial and error method, which is the standard approach for many optimization problems [6, 22]. PF is the penalty function, which is necessary to preserve closed-loop system stability during executing the optimization procedure on the one hand and to limit the absolute values of the gains in (6) in order to avoid the aforementioned saturation. It consists of three parts [6, 22]
The first one as a function of minimum distance to the imaginary axis could be graphically shown in Figure 2 and defined over the area for its first border (see Figure 2) as follows [6, 22]. where is a very large value (for example ), is the distance between the maximum and the minimum real parts of the poles. defines minimum stability margin of the closed-loop system.
The penalty function for the second border, which defines the closed-loop system bandwidth, looks like (12). The third component PF3 has to restrict too large values of gains in (6), and it was chosen in the following form.
This component plays the same role as the weighting matrix in expression (2) for standard minimization procedure . It is useful to mention that restriction imposed on controller (9) gains and by cost function (13) constrains the bandwidth of the closed-loop system. Eventually, we can formulate the optimization procedure as follows. where is the admissible domain in the controller parameters plane determined by PF, and is the performance-robustness index defined by (10), which is minimized by the function “fminsearch.m” in MATLAB. As it is known, this function uses Nelder-Mead minimization procedure. After running this procedure, we obtained the following parameters of the controller : . The values of norms are equal to . These norms demonstrate a good ability to suppress external disturbance [16, 17, 23]. The small value of the norm of the transfer function shows high degree of robustness of ISP control. Note that these values are upper estimates, and simulation demonstrates even higher efficiency of the controller (6). The eigenvalues of the closed-loop system equal . Then, this system is the very “stiff” one, as we noted before. Poles of the plant (5) substantiate the high-frequency complex poles of the closed-loop system, and low-frequency real pole appears due to the integral gain of controller.
Now, we will determine the gain of the feedforward controller “Ff.con.” (see Figure 1). Under [8, 9], the invariance condition concludes in the following: two aforementioned channels of the disturbance propagation must be equal and have opposite signs. Applying this condition to the system (7), one can obtain
Introducing the pseudoinverse matrix , it is possible to find gain as follows
2.2. Simulation of ISP with and without Feedforward Gain
The first stage of this process is the estimation of efficiency of ISP disturbance suppression in stabilization mode with a zero reference signal. Figure 3(a) represents the disturbance plot, and Figure 3(b) represents the plot of output ISP signal or the rotation rate of the azimuth frame of the gimbals. The rms of the input disturbance equals ; meanwhile, the rms of the ISP rotation rate is much less: . Therefore, the efficiency of this ISP is good enough, because it suppresses the disturbance times. However, including the feedforward gain in the control law much more improves ISP efficiency for suppressing input disturbance. Figure 4(a) represents the same plot of disturbance , and Figure 4(b) represents the plot of the rotation rate of the ISP azimuth frame of gimbals with feedforward controller. The rms of input disturbance equals the same , but the rms of ISP rotation rate is much smaller: . The input disturbance is suppressed times, so the efficiency of ISP increases at more than 71 times by the simplest mean. The experimentally obtained plot demonstrates a sharp peak in the proximity of point (, ). So, using feedforward/feedback control leads to improvement of stabilization accuracy.
The second stage of the modeling is the estimation of the ISP efficiency in the tracking mode with the nonzero reference signal, which is accepted as step function having magnitude 0.25 rad/sec (14.325 deg/sec). The transient processes, in this case, are shown in Figures 5 and 6. Plots in Figure 5 demonstrate the fast response of the ISP control systems (Figure 5(a)) and the absence of saturation of the control input (Figure 5(b)). Plots in Figure 6 show the ISP rotation rate error in the tracking mode without feedforward (Figure 6(a)) and with feedforward (Figure 6(b)). It is obvious that including feedforward gain practically eliminates the jitter of the target image in the field of sight of the stabilized camera in the tracking mode, thus essentially improving the ISP performance.
The design of effective but simple control laws for the airborne ISP was one of the initial goals of this research. However, further progress in the microelectronics opens up the possibility to use more effective ISP control laws. The new studies and papers will be devoted to the possibility of using other types of controllers for disturbance rejection and comparing them with the represented controller.
The most widespread sensors applicable to the designed controller are MEMS gyroscopes. For more expensive applications, it is possible to use Coriolis vibratory gyroscopes. The choice of actuators depends on a mass of loading. For the small mass, it is convenient to use moment electric motors; in other cases, it is necessary to use electric motors with reducers. In both cases, electric motors of direct current are used.
Usage of the rate gyroscopes in the inertial measurement unit (IMU) of the airborne navigation equipment gives an additional source to increase the ISP accuracy because it gives the possibility to include the feedforward disturbance gain in the control law following the main principles of the invariance theory.
The design of a simple PI control law based on the parametric robust optimization provides a quite acceptable accuracy and robustness of the ISP feedback control system.
Augmenting PI feedback control law with feedforward disturbance gain significantly improves the accuracy of inertial platform stabilization to such degree, which is inaccessible for feedback only even when the more sophisticated control law is used .
The results of the simulation of the ISP control systems with and without the feedforward disturbance gain prove the efficiency of the proposed algorithms for the stabilization of the target image in the camera’s field of view.
The paper “Design of Composite Feedback and Feedforward Control Law for Aircraft Inertially Stabilized Platforms” data used to support the findings of this study are included within the article in the list of references.
Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this paper.
The authors are grateful to the Ministry of Education and Science of Ukraine. This study was carried out in the frames of the research work “Methodology of Designing High-Precision Stabilized Platforms for Aircraft Observation Equipment using the Perspective Inertial Technologies” supported by the grant of the Ministry.
J. M. Hilkert, “Inertially stabilized platform technology,” IEEE Control Systems Magazine, vol. 1, pp. 26–46, 2008.View at: Google Scholar
E. Altug, J. P. Ostrowsky, and C. P. Taylor, “Control of a quadrotor helicopter using dual camera visual feedback,” The International Journal of Robotics Research, vol. 24, no. 5, pp. 329–341, 2016.View at: Google Scholar
L. R. Garsia Carillo, E. Rondon, A. Sanchez, A. Dzul, and R. Lozano, “Stabilization and trajectory tracking of a quad-rotor using vision,” Journal of Intelligent Robot Systems, vol. 61, pp. 103–118, 2011.View at: Google Scholar
A. A. Tunik and O. A. Suchchenko, “Usage of vector parametric pptimization for robust stabilization of ground vehicles information-measuring devices,” Proceedings of the National Aviation University, vol. 57, no. 4, pp. 23–32, 2013.View at: Google Scholar
V. B. Larin, A. A. Tunik, and S. I. Ilnytska, Some Algorithms for Unmanned Aerial Vehicles Navigation Systems, Outskirts Press, Inc., Colorado (USA), 2019.
B. N. Petrov and A. I. Kukhtenko, “Theory of invariant systems design, Chapter I in the book,” in Modern Methods of Automatic Control Systems Design, B. N. Petrov, V. V. Solodovnikov, and Y. I. Topcheev, Eds., p. 703, Moscow, Mashinostroenie, 1967.View at: Google Scholar
G. Tian and Z. Gao, “From poncelet’s invariance principle to active disturbance rejection,” in 2009 American Control Conference, pp. 2451–2457, Hyatt Regency Riverfront, St. Louis, MO, USA, June 2009.View at: Google Scholar
O. Katsuhiko, Modern Control Engineering, Pearson, 5th edition, 2009.
Y. Luo, W. Ren, Y. Huang et al., “Feedforward control based on error and disturbance observation for the CCD and fiber-optic gyroscope-based mobile optoelectronic tracking system,” Electronics, vol. 7, no. 223, p. 18, 2018.View at: Google Scholar
M. Kara-Mohamed, A. Lanzon, and W. P. Heath, “Feedforward/feedback multivariable control design for high speed nanopositioning,” in Proceedings of 2014 European Control Conference (ECC), pp. 1939–1944, Strasbourg, France, June 2014.View at: Google Scholar
M. O. Rauw, The flight dynamics and control toolbox, Math Works Company, 2000.
D. McLean, Automatic Flight Control Systems, Prentice Hall Inc., 1990.
R. C. Dorf and R. H. Bishop, “Modern Control Systems,” in Modern Control Systems, p. 1032, Pearson, 13th edition, 2017.View at: Google Scholar
V. B. Larin and A. A. Tunik, “Flight stabilization and exogenous uncertain disturbance suppression via static and dynamic output feedback,” in Proceedings of 13th St. Petersburg International Conference on Integrated Navigation Systems, pp. 27–35, St. Petersburg, Russia, May 2006.View at: Google Scholar
“Electrical Motor EDM-20,” http://tehnika.agroserver.ru/zapasnye-chasti/elektrodvigatel-edm-20-64010.htm.View at: Google Scholar
I. E. Kazakov and B. G. Dostupov, Statistical Dynamics of Nonlinear Automatic Systems, Fizmatgiz, Moscow, 1962, 332 p.
A. Gelb, J. F. Kasper Jr., R. A. Nash Jr., C. F. Price, and A. A. Sutherland Jr., Applied Optimal Estimation, MIT Press, Cambridge, MA, 2001.
E. J. Davison and I. J. Ferguson, “The design of controllers for the multivariable robust servomechanism problem using parameter optimization methods,” IEEE Transaction on Automatic Control, vol. AC-26, no. 1, pp. 93–110, 1981.View at: Google Scholar
A. A. Tunik and T. A. Galaguz, “Robust stabilization and nominal performance of the flight control system for small UAV,” Applied and Computational Mathematics, vol. 3, no. 1, pp. 34–45, 2004.View at: Google Scholar