ISSN ONLINE(2278-8875) PRINT (2320-3765)
M.S. Das1, S. K. Singh2
|
Related article at Pubmed, Scholar Google |
Visit for more related articles at International Journal of Advanced Research in Electrical, Electronics and Instrumentation Engineering
In this paper unknown input observer using projection operator method is used to estimate five states of the fifth order lateral axis model of L-1011 system in discrete domain. Discrete reduced order UI observer is used to estimate the states of the system in noiseless and noisy environment. Kalman filter has been designed to estimate the state of L-1011 system when noise is considered. Finally error between the actual and estimated states of the system by UI observer and Kalman filter i.e performance of both are compared. It has been shown that Kalman filter is good estimator than UI observer in noisy environment in respect of estimation. The simulations have been done using MATLAB and SIMULINK Toolbox & the results of the simulations have been displayed in this paper.
Keywords |
Unknown input observer (UIO), projection operator method, Estimated state, Estimated error, Kalman filter. |
INTRODUCTION |
Observer and Kalman filter is widely used in estimation of different linear and non linear systems. Many authors have designed unknown input observer in different approach and Kalman filter for estimation of linear and non linear systems. Over the past three decades, many significant works have been carried out on the construction of observers for linear and nonlinear systems in the control systems literature[1-15].The observer was first proposed and developed by Luenberger in the early sixties of the last century [1][2][4].Full order and reduced order observer was introduced to estimate the states of the linear and non linear systems[12][16].Sliding mode observer were proposed for uncertain systems[7][13][14]. Later, Sliding mode observer was introduced for fault detection and isolation [14-15]. Observers for systems with unknown inputs has an important role in robust model-based fault detection [14][15].Alexander Stotsky and Ilya Kolmanovsky designed an observer to estimate the unknown input from available state measurements in automotive control application [18]. Talel. Bessaoudi, Karim Khemiri, Faycal. Ben Hmida and Moncef. Gossa [19] estimate the unknown input and states of a linear discrete time systems using recursive least-square approach.Estimation of states and unknown input of a nonlinear communication system has been addressed [20].Kalyana C. Veluvolu and Soh Yeng Chai designed a high gain observer with multiple sliding mode for state and unknown input estimations[21]. Simultaneous estimation [22] of states and unknown input for a class of nonlinear systems has been proposed by Q. P. Ha and H. Trinh. Thierry Floquet and Jean-Pierre Barbot designed a state and unknown input delayed estimator for discrete-time linear systems [23]. Abhijit Banerjee and Prof. G. Das [24] used the reduced order observer [26] to estimate the unknown input of a linear time invariant system. Ashis De, Abhijit. Banerjee and Prof. G. Das [25] estimate the unknown input of an LTI system using full order observer constructed by the method of generalized matrix inverse [27]. |
In 1960, R.E. Kalman published his famous paper [28] describing a recursive solution to the discrete-data linear filtering problem. Since then,the Kalman filter has been the subject of extensive research and application, particularly in the area of autonomous or assisted navigation.Many researchers are using different Kalman filter like Discrete Kalman filter Extended Kalman filter[31] ,Unscented Kalman for estimation of tracking[31-32], military system, flight control system, Image and video tracking in a real time system etc.The filter is said to be very powerful as it can estimate the state of a process even when the precise model of the system is unknown. A Kalman filter is simply an optimal recursive data processing algorithm that blends all available information, including measurement outputs, prior knowledge about the system and measuring sensors, to estimate the state variables in such a manner that the error is statistically minimized In practice, linear equation system with white Gaussian noises is commonly taken as the standard model of a Kalman filter. |
Stefen Hui and Stanislaw H. Zak [17] designed both unknown input full order and reduced order observer using projection operator method in continuous domain. They did not estimate the states of the L-1101system in noisy environment. In this paper a discrete reduced order UI observer using projection operator has been designed to estimate the states of L-1011 system in noisy environment. Also discrete Kalman filter is used to estimate the states. Finally accuracy in estimation by discrete Kalman filter and observer is compared. |
II. SYSTEM’S MODEL AND DESCRIPTION |
![]() |
III.DISCRETISED MODEL OF REDUCED ORDER UNKNOWN INPUT OBSERVER |
At first the continuous fifth order lateral axis model of L-1101 system and reduced order unknown input observer are discretized .This discrete UI observer is used to estimate the state of discrete system ( fifth order lateral axis model of L-1011 system).Firstly it is applied on the system to estimate the state in noisy and without noisy condition. Here it has been seen that Sampling time is taken 0.01 sec during discretisation .Dynamic Model of the discretised L-1011 system and reduced order unknown input observer using projection operator is given below respectively . |
X(k+1) = Ad x(k) + Bd u(k) and Y (k) = Cd x (k) + Dd u (k); Where |
![]() |
IV. SIMULATION RESULT (WITHOUT NOISE) |
Actual state and estimated state by UI observer using projection operator method [1] and error are shown below |
![]() |
![]() |
![]() |
V.ESTIMATION BY REDUCED ORDER UI OBSERVER WHEN NOISE IS INTRODUCED IN THE SYSTEM |
All systems are affected by unwanted noise that changes the actual output or actual state. Here it is assumed that the noise that is introduced with the system (Fifth order lateral axis model of L-1011) is zero mean white Gaussian in nature.UIO is applied on noisy system to estimate state of the system. |
VII. SIMULATION RESULT (IN NOISY ENVIRONMENT) |
Actual state of the noisy system estimated state by reduced order observer using projection operator method in discrete domain has been shown in below |
![]() |
![]() |
VII. STATE ESTIMATION BY DISCRETE KALMAN FILTER |
The Kalman filter is a mathematical method named after Rudolf E. Kalman.The Kalman filter produces estimates of the true values of measurements and their associated calculated values by predicting a value, estimating the uncertainty of the predicted value, and computing a weighted average of the predicted value and the measured value. Weight is chosen such that uncertainty is minimum. It is an essential part of development part of space and military technology. The Kalman filter is very important tool which is used to estimate the state of noisy linear dynamic system, described in state space form. It provides a recursive solution to the linear optimal filtering problem. It can be applied in stationary and non stationary environments. The solution is recursive in that each updated estimate of the state is computed from the previous estimate and the new input data, so only the previous estimate requires memory. In this case it is essential to store only the previous estimate instead of storing entire past state.Kalman filter is computationally more efficient than computing the estimate directly from the entire past observed data at each step of the filtering process.The state is the least amount of data on the past behavior of the system that is needed to predict its future behavior.Typically, the state xk is unknown. To estimate it, we use a set of observed data, denoted by the vector yk.Process equation can be written as |
x(k+1)= Φ(k) x(k) + Bu(k) + wk |
Where Φk is the state transition matrix of nxn dimension taking the state xk from time k to time k + 1. The process noise wk is assumed to be additive, white,and Gaussian, with zero mean wk ~ N(0,Q), Where Q process noise covarinace matrix. and with process noise covariance matrix defined by E(wn* wT k) = Qk when n=k otherwise zero. Measurement equation can be written as |
Zk=Hxk+vk |
where, vk is the measurement noise of m x 1 vector which is denoted by vk~N(0,R) where R measurement noise covariance matrix. xk = n×1 process state vector at time tk. zk = m × 1 measurement at time tk and C = Output matrix. Kalman Filter estimate state in two phase one is predict and another is update .Predict Equations and update equations are given below. In predict phase it uses the state estimate from the previous time step to produce an estimate of the state at the current time step. |
Predicted state (a priori state) : ïÃâ¬Ã xˆk =A ˆk ïÃâ¬ÃÂ1 x +B u k-1 And predicted(a priori state) error covariance Pk= APk-1 AÃÂÃÂ+Q |
The current a priori prediction is combined with current observation information to refine the state estimate. This improved estimate is termed as posteriori state estimate. Typically, the two phases alternate, with the prediction advancing the state until the next scheduled observation and the update incorporating the observation.Updates equations are given below.. |
![]() |
VIII. ESTIMATION BY KALMAN FILTER |
The discrete Kalman filter is applied on noisy fifth order lateral axis model o L-1011 system to estimate the state of the system. Here it is assumed zero mean white Gaussian noise is mixed with system and output equation .It means process noise Wk and measurement noise Vk all are zero mean white Gaussian noise. So process noise Wk and measurement noise Vk can be represented as Wk ~N(0,Qk) and Vk ~N(0,Rk) Where Qk is process noise covariance matrix and Rk measurement noise covariance matrix.It is assumed that process noise covariance matrix Qk and measurement noise covariance matrix Rk is identity matrix. The value of Qk and Rk is given below. |
![]() |
IX. SIMULATION RESULTS |
![]() |
![]() |
X.COMPARISON OF ESTIMATED ERROR BETWEEN KALMAN FILTER AND UI OBSERVER (PROJECTION OPERATOR METHOD) |
In this thesis the comparison of Kalman Filter and Unknown Input Observer is shown. These Two techniques ie discrete Kalman Filter and discretised Unknown Input Observer (using projection operator approach) are applied on a same noisy system to estimate the state of the system. These two techniques are applied on the fifth-order lateral axis model of an L-1011 fixed wing aircraft model. Here it is considered that zero mean white Gaussian noise is introduced to the system .First Unknown Input observer and secondly Kalman Filter is applied on same the noisy system to estimate the state and hence estimated error is calculated for each cases .Estimated state error by UI observer and The Kalman Filter is given below. There are five state in the system i.e five error plot is shown below. Five states are as follows x1=Φ , bank angle (rad), x2=yaw rate (rad/sec), x3=roll rate (rad/sec) , x4=sideslip-angle (rad) and x5=washed out filter state. |
XI.SIMULATION RESULT |
Actual state of a noisy system , estimated state by UI observer using projection operator method and estimated state by Kalman filter in noisy environment are shown in the following figure. |
![]() |
![]() |
![]() |
COMPARISON OF ERROR OF KALMAN FILTER AND UI OBSERVER |
Errors in estimation of states by Kalman Filter and UI reduced order observer are compared . Comparison of errors are shown in the following figure. |
![]() |
![]() |
RMS ERROR COMPARISON |
Finally I have compared the RMS error of UI observer and RMS error of Kalman Filter and error covariance. Simulated outputs are shown in the following figures. |
![]() |
![]() |
XII.CONCLUSION |
This paper summarizes a number of conclusions that can be drawn from the work . Reduced order unknown input observer is used here to estimate the state of the fifth order lateral axis model of an aircraft model of L-1011 in continuous and discrete domain. Also it is seen that UIO is capable of estimating state of a system in noisy environment in discrete domain. Here initial value of UIO observer q(0) is assumed to be zero.It is seen that estimated error using UIO using projection operator method is low in noise free system.It is a good estimator for noise free system. But it gives poor performance in estimation in noisy environment. Finally discrete Kalman filter is used to estimate the state of the fifth order lateral axis model of an aircraft model of L-1011 system in noisy environment. Kalman filter is good minimum variance estimator and it gives satisfactory result in noisy environment.It is easy to conclude from the error plot ( in fig. 21 to fig.25) and RMSE plot( in fig.26 to fig.30) that Kalman fillter is more efficient estimator than unknown input observer in noisy field. There is some limitation of using Kalman filter .We should know mean and correlation of the process noise Wk and measurement noise vk at each time instant.We should know process noise covariance Q , measurement noise covariance R and error covariance P also. If it is not given then estimation by Kalman filter is not possible correctly. In this paper process noise Wk and measurement noise vk is assumed to be zero mean white Gaussian noise .The value of P,Q,R is assumed as identity matrix in this paper.We have applied Kalman filter and UI observer to estimate state for the linear system only not for non-linear system. In future Extended Kalman Filter ,Unscented Kalman Filter ,UI observer ,Sliding mode observer, particle filter and Disturbance observer will be applied to estimate states of non-linear systems and compared their performance in estimation. |
ACKNOWLEDGEMENT |
I would like to thank to all professors of the department of Electrical Engineering, Kolkata, Jadavpur University and particularly to my guide, my Professor, Dr. Smita Sadhu . I am greatly indebted to her for his valuable advice, suggestion, kind co-operation, liberal help and guidance extended throughout the preparation of my work. It is my proud privilege to thank Professor Dr. T.K Ghoshal for encouragingto work. |
References |
|