Keywords
Kalman filter, Twin Rotor MIMO system, analytical redundancy, neural network, Luenberger observer, soft sensing
This article is included in the Manipal Academy of Higher Education gateway.
Kalman filter, Twin Rotor MIMO system, analytical redundancy, neural network, Luenberger observer, soft sensing
Study of aerospace systems has always been a subject of interest by many researchers, engineers, and technical students. However, it is practically very difficult to analyse these aerospace systems, so a replica is designed to understand the behaviour of the actual systems. One such system, which replicates the behaviour of a helicopter system, is a TRMS. The TRMS has two of the three movements of a helicopter, pitch angle and yaw angle.1 The TRMS position is controlled by the rotor speed. The two input variables to the TRMS are voltage to the main rotor and tail rotor, and the outputs are the pitch and yaw angle as shown in Figure 1. Excitation to these motors is given by a controller based on a set point given by user. Desired control action can be achieved only by accurate measurements of pitch and yaw.
Several researchers have reported work on computing the yaw and pitch positions using different techniques.1–5 Different algorithms are incorporated to maintain and retain the stability of the TRMS. Controllers like Proportional Integral Derivative (PID), fuzzy PID, sliding mode controller, fuzzy sliding mode controllers and estimators like Luenberger and Kalman are incorporated on the TRMS.2 Rahideh and Shaheed3 discussed the design of a model predictive controller for a TRMS. The controller is simulated based on the state model using motor armature current and angular speed as inputs. Jahed and Farrokhi4 discussed the design of a fuzzy based robust control for a TRMS using angular speed given by tachometer in a simulation platform, and Tao et al.5 discussed the fuzzy sliding and integral sliding controller design for a TRMS. In Rohith,6 a new control law was proposed for the design of sliding motor controllers, which mitigate the chattering problem with the gain variation and thereby guarantees faster system response and robustness. Design of an auto tuning based PID controller with fractional-order reference model approximation for a DC rotor in a TRMS model is discussed in Alagoz et al.7 The data regarding the process variable is derived from the angular position of pitch and yaw. In Patel and Janardhanan,8 a tuning method based on the Moore Skelboe Algorithm for a PID controller was presented. This was effective in finding the optimal PID values for the given initial range and applicable for some of the higher-order Linear Time Invariant (LTI) systems as well as for all first and second order linear time invariant systems. Further the method was effective in stabilizing unstable systems. Halim and Ismail9 presented a PID controller tuning using tree physiology optimization, which was based on the tree growth concept whose simulation results showed better results compared to other tuning methods in Single Input Single Output (SISO) and MIMO problems. Rao et al.,10 reported design of an observer using robust PID controller logic with H∞ observer to obtain the stable output in TRMS with sensor and actuator failure. Netto et al.,11 reported an Adaptive PID controller to cancel the effect of cross coupling between the tail rotor and main rotor when operating simultaneously in a TRMS. Adaptive linear quadratic regulator design for stable system operating at a single reference point with six tail and main rotors is reported by Faisal and Omar Waleed.12 Ghellab et al.,13 reported an adaptive radial basis function neural network with a dynamic terminal sliding mode control with cross coupling between tail and main subsystem for tracking the set point in the presence of wind gust and other external disturbance.
In Sleimi et al.14 a linear time-varying controller was designed using a differential flatness property leading to a two Degree of Freedom (DOF) controller for which the system must be in canonical controllable form with no need to define its dynamics. In Panda et al.15 the proposed control strategy used an adaptive back stepping controller implemented on a Twin Rotor Multi Input Multi Output System. It provided an explicit relationship between the saturation bound of the input signal and upper bounds of tracking errors, uncertainties, and disturbances. Mondal and Dey16 presented the development of a two DOF control system design providing an additional degree of freedom depending on the nature of the plant and loop compensators. The design methodology can be implemented for integer as well as non-integer order plants with better tracking and loop robustness.
Neural network based controller design using model inversion control for a twin rotor MIMO system was reported in Rahideh et al.17 Design of a differential evolution based neural network model to control the TRMS with data of angular positions was reported by Subudhi and Jena.18 Pratap and Purwar,19 reported implementation of a neuro-adaptive robust back stepping controller for TRMS. Many researchers have worked on the TRMS, but the design technique of most controllers is based on the data received from the sensor. If the sensor data is erroneous the complete system fails. To overcome the above problem, this paper attempts to design a soft sensor or observer for computing the actual yaw and pitch data to the controller. In order to start the design of a soft sensor it is essential to understand the TRMS model.
The TRMS dynamics are given in Rahideh et al.,20 Sun and Song21 and Ahmad et al.22 It was found that the system was a nonlinear and coupled one with considerable error in the measurements. The method based on the Kalman filter is considered as it provides for updating the estimates based on errors, and also takes into account the process and the measurement noise uncertainties to give out the best estimate. Model dynamics of the system are known and are used for estimation. The Kalman filter is widely used as a state estimator and is used in various fields. For example, the Kalman filter is used to estimate track soil parameters,23 for tracking objects in videos based on color features and mean shift,24 for estimation of flow rate and isolation of occurrence of fault,25 for estimation of wind speed and comparison with the estimation from Takagi-Sugeno observers.26
The TRMS is a prototype device used to understand the dynamics of a helicopter system. It consists of two rotor fans to operate the device in two Degrees of Freedom (2-DoF). For controlling/stabilizing the system performance based on a desired set point, it is essential to sense the actual pitch and yaw positions. The presence of faults during data transmission to the controller lead to the controller failing to take the desired action and thus the system destabilizes. Hence, a technique is proposed with the objective of designing an estimator so that it can be used in case of any contingencies. The Kalman filter technique is used to design the estimators for estimating the pitch and yaw position of the TRMS. In order to verify the effectiveness of the algorithm, it is applied to the system in real life and compared with the actual trajectory of the response. The performance of the designed Kalman filter for estimating angular pitch and yaw information is compared with the Luenberger observer and neural network estimator to evaluate its performance. The outline of the work carried out is shown in Figure 2.
The twin rotor MIMO system is a prototype of a helicopter propeller system. It is shown in Figure 3. The system was a non-linear MIMO system with significant cross-coupling. The angle of attack of the rotors was fixed and the aerodynamic forces were controlled by varying the speed of the motors. Significant cross-coupling was observed between the actions of the rotors, with each rotor influencing both angle positions. Two propellers were driven by DC motors controlled by their supply voltages. The measured signals were positioned on the beam in the space i.e. two position angles which are measured by rotary optical encoders mounted on each of the rotor shafts. The measured output was obtained from the rotary optical encoders. It was calibrated for the inertial axis using the loop up table approach.19, 20
The TRMS has two modes of operation, namely, 1-DOF control and 2-DOF control. In 1-DOF control, the pitch and the yaw are controlled individually and independently, whereas in 2-DOF control the pitch and the yaw are controlled simultaneously in a coupled system. The block diagram of the proposed scheme is shown in Figure 4. Here the coupled TRMS was controlled by two separate PID controllers for the pitch and the yaw respectively. The controllers were tuned to maintain the angles at the desired positions or set point, except in the event of some contingency such that the pitch or the yaw failed to maintain the desired position. The transient output of the sensors was used to design an estimator to find the angles where models were generated using a system identification process. A decoupled pitch and yaw model was used for designing the estimator. In this case pitch and yaw angles were independently calculated using the Kalman filter.
The Kalman filter is developed from the Bayesian filter. It was initially used to extract the signal from the noisy output data of sensors, and/or actuators.27–29 Over time, it has been used as a state estimator also. It is a stochastic estimator and requires explicit modeling of the process noise and the measurement noise in addition to the system model. The Kalman filter estimation is a two-step process. Initially it generates an estimate from the knowledge of the system dynamics, which is embedded in the system model along with the noise model. This is called the apriori estimate. Once the measurement output is available, the apriori estimate is updated to the posteriori estimate. The second stage involves the Kalman gain which is altered in every step based on an optimization problem. The Kalman gain was the main feature of this estimator and it decided if the estimate derived much of its information from the measurement output or the system dynamics. Hence, the Kalman gain was also updated on every step such that the error between the actual output and the estimate was reduced i.e. the residue reduced to zero. The noise covariance matrix was also updated at each step based on the error. Here the Markovian model was used where the state was said to be Markovian but the measurement was usually not Markovian. Also, we considered the process noise and measurement noise as Gaussian and they were not correlated with each other. It is also called the Gauss-Markov model. This allows nonlinearity. The mean was assumed to be zero for both measurement and process noises. Covariance was chosen in this case such that the covariance of measurement noise was slightly less than the covariance of the process noise. In general, the probability distribution of the state is given by equation (1)
where was determined by which was the probability distribution of the process noise. The system state was modeled as a linear combination of the previous states along with the input U and the process noise W, as given in equation (2)
For the considered TRMS, the matrices for the decoupled pitch model were found to be:
The matrices for the decoupled yaw model were found to be:
The measurement obtained was modeled as a linear combination of the system states and measurement noise V as in equation (5)
In the TRMS model, the measured output i.e. pitch and yaw, were related to states directly through output matrices S i.e. , . Hence for pitch, . For Yaw,. Let Q and R be processed noise covariance and measurement noise covariance respectively. The apriori error covariance P is given by equation (6).
Once the measurement was obtained, the Kalman gain was computed and the apriori state estimates and the apriori error covariance were updated as follows:
The posteriori error covariance
Here was the estimate of the state at k given the measurement at k and the apriori estimate at time k. Also, it should be noted that the Kalman filter was a one step ahead predictor. The noise covariance matrices R and Q were the most complex matrices to compute, therefore the idea was to start with an initial estimate of identity matrices for both R and Q and to change them based on the convergence between the actual and the estimated outputs. The Kalman filter estimated the states. However, the objective was to estimate the pitch and the yaw which were the output of the TRMS. Hence, the output was obtained from equation (10).
and the output residue (error) was given by equation (11)
where was the actual value of pitch and yaw
The Kalman filter provided a stochastic method of estimation, whereas the Luenberger observer and neural network provided a deterministic method of estimation of the states. The estimation of yaw and pitch carried out with these deterministic methods had a poorer performance due to their inability to incorporate uncertainties that emerge inherently in any real systems, such as the uncertainties due to modeling errors or sensor output errors. However, these estimates can be feasible if the error margin is not very stringent.
Luenberger based observers are widely used for numerous applications and modified according to the type of the system.27 The extended Luenberger observer and adaptive Luenberger observer are the latest class of these observers.27 The Luenberger observer is used for estimation of battery charge for electric vehicles,30 flux in motor drives,31 motorcycle dynamics,32 sensor-less speed estimation.33 The scheme of this observer is shown in Figure 5. The following were the expressions for the observer design from Figure 5.
State equations of the observer:
From equation 17, it is understood that if the eigenvalues were all negative, the estimated state vector error, ex, would decay to zero. The design then consisted of solving for the values of L to yield a desired characteristic equation.
Then the selection of eigenvalues for the observer was carried out to achieve a required closed loop response. These eigenvalues determined a characteristic equation that was made equal to equation 19 to solve for L. The poles for the observer for both pitch and yaw were chosen to be at -1, -2, and -3.
The Luenberger gain L for pitch was:
The Luenberger gain L for yaw was:
The estimation based on the neural network was based on a time series correlation between the input and the output, also called the targets. This assumed a black body model, where the dynamics of the system were not explicitly parameterized and the forecasting was done only from the input output relationship of the system. Here the neural network was trained to replicate the behavior of the system by using the given set of input output pair of data. The Levenberg-Marquardt algorithm with a back-propagation network having one hidden layer with a size of 10 neurons, which was used to predict the output i.e. the pitch and the yaw from the input. The neural network used for the present work is shown in Figure 6. The ‘nftool’ of MATLAB version R2014a,34 provided the platform to train the neural network and this was used in the current work. The regression graph for the trained neural network is shown in Figure 7 for pitch and Figure 8 for yaw systems.
A soft sensor for computing the yaw and pitch in TRMS using three different techniques, a neural network, Kalman filter, and a Luenberger observer, was designed. It was subjected to tests in real life. Tests were conducted by applying input and disturbance. The response of the proposed sensing system using the Luenberger observer is represented in Figures 9-12. The output obtained for the Kalman filter is shown in Figures 13 and 14. Similarly, the output obtained from the neural network is shown in Figures 15 and 16. From the responses it was found that the proposed sensing technique was able to track the pitch and yaw positions accurately in a practical system. The performance measures IAE, ISE, and ITAE were used to quantitatively compare the outputs obtained from the Luenberger observer, Kalman filter, and neural network estimator. The results of these for the pitch and yaw measurement are shown in Figures 17 and 18, respectively. This method provided a way to classify the errors that occurred in different stages of the system operation and provided a means to judge the accuracy of estimation. The estimator having the highest number of errors was not accurate. The steady state of the system occurred at around time = 50 seconds, and the errors at that instant for all the estimators were tabulated for pitch and yaw measurements in Table 1.35
Twin rotor multi input systems are fundamental for any aerial system, analysis of any property on the TRMS makes it important when designing any system further. The TRMS contains a motor system along with a sensor system for testing. Sensors are an integral part of the TRMS for providing actual information to the controller. Failure of the sensors would lead to errors in control action and thus lead to system failure. An attempt was made in this work to design an observer-based system, which would function accurately even with sensor faults. The effectiveness of the Kalman filter algorithm for estimation of pitch and yaw angular positions was verified in real life on the TRMS with external disturbances. The estimate was updated at every sampling time to predict the yaw and pitch for angular positions for the given input and measurement data. The performance of the Kalman filter was compared with that of the neural network and Luenberger observer. From the results of IAE, ISE, and ITAE, the error was least when using the Kalman filter followed by the Luenberger observer and lastly, the neural network. Hence, the estimation by the Kalman filter was more accurate during both transient and steady states.
Open Science Framework: Extended data for ‘Design of a soft sensing technique for measuring pitch and yaw angular positions for a Twin Rotor MIMI System’, https://doi.org/10.17605/OSF.IO/52V9D.35
This project contains the following extended data:
• Supplementary Video 1: Experimental model of a soft sensor design technique for estimation of pitch and yaw angular positions of a Twin Rotor MIMO System (TRMS)
• Supplementary Data 1: Simulink file used to carry out the real time experimentation of TRMS system.
Data are available under the terms of Creative Commons Zero “No rights reserved” data waiver (CC0 1.0 Public domain dedication).
Views | Downloads | |
---|---|---|
F1000Research | - | - |
PubMed Central
Data from PMC are received and updated monthly.
|
- | - |
Is the work clearly and accurately presented and does it cite the current literature?
Yes
Is the study design appropriate and is the work technically sound?
Yes
Are sufficient details of methods and analysis provided to allow replication by others?
Yes
If applicable, is the statistical analysis and its interpretation appropriate?
Partly
Are all the source data underlying the results available to ensure full reproducibility?
Yes
Are the conclusions drawn adequately supported by the results?
Partly
Competing Interests: No competing interests were disclosed.
Reviewer Expertise: Solar energy, Energy modelling, Energy systems
Is the work clearly and accurately presented and does it cite the current literature?
Partly
Is the study design appropriate and is the work technically sound?
Yes
Are sufficient details of methods and analysis provided to allow replication by others?
Partly
If applicable, is the statistical analysis and its interpretation appropriate?
Not applicable
Are all the source data underlying the results available to ensure full reproducibility?
Yes
Are the conclusions drawn adequately supported by the results?
Yes
References
1. Tastemirov A, Lecchini-Visintini A, Morales-Viviescas R: Complete dynamic model of the Twin Rotor MIMO System (TRMS) with experimental validation. Control Engineering Practice. 2017; 66: 89-98 Publisher Full TextCompeting Interests: No competing interests were disclosed.
Reviewer Expertise: Control, system identification, mechatronics systems, automotive
Alongside their report, reviewers assign a status to the article:
Invited Reviewers | ||
---|---|---|
1 | 2 | |
Version 2 (revision) 25 Nov 21 |
read | read |
Version 1 05 May 21 |
read | read |
Provide sufficient details of any financial or non-financial competing interests to enable users to assess whether your comments might lead a reasonable person to question your impartiality. Consider the following examples, but note that this is not an exhaustive list:
Sign up for content alerts and receive a weekly or monthly email with all newly published articles
Already registered? Sign in
The email address should be the one you originally registered with F1000.
You registered with F1000 via Google, so we cannot reset your password.
To sign in, please click here.
If you still need help with your Google account password, please click here.
You registered with F1000 via Facebook, so we cannot reset your password.
To sign in, please click here.
If you still need help with your Facebook account password, please click here.
If your email address is registered with us, we will email you instructions to reset your password.
If you think you should have received this email but it has not arrived, please check your spam filters and/or contact for further assistance.
Comments on this article Comments (0)