Advanced
Robust Dissipative Filtering for Polytopic Uncertain Singular Systems in Continuous & Discrete Time
Robust Dissipative Filtering for Polytopic Uncertain Singular Systems in Continuous & Discrete Time
The Transactions of The Korean Institute of Electrical Engineers. 2014. Jul, 63(7): 994-1000
Copyright © 2014, The Korean Institute of Electrical Engineers
  • Received : May 22, 2014
  • Accepted : June 16, 2014
  • Published : July 01, 2014
Download
PDF
e-PUB
PubReader
PPT
Export by style
Article
Author
Metrics
Cited by
TagCloud
About the Authors
데니 진
Dept. of Electronic Engineering, Sun Moon University, Korea
종해 김
Dept. of Electronic Engineering, Sun Moon University, Korea
도창 오
Corresponding Author : Dept. of Medical IT Engineering, Konyang University, Korea E-mail :docoh@konyang.ac.kr

Abstract
This paper considers the problem of continuous and discrete time dissipative filter design method for the singular systems with polytopic uncertainties. Two bounded real lemmas (BRL) for the robust filters of dissipative singular systems along with the polytopic uncertainties are proposed on the basis of Lyapunov criterion. The sufficient conditions for both continuous and discrete time dissipative filter design methods are derived using the obtained BRL by linear matrix inequality (LMI) approach. Finally, the validities of the proposed methods are shown by numerical examples.
Keywords
1. 서 론
Dissipative System takes on an important role when considering a flow of any entity. Matter such as liquid, gas, or energy could be considered as a dissipative system [1] , as described and modeled by many scientists. In the early 70's, one of the vikings for dissipative system, J. C. Willems [1] , wrote the general theory about dissipativeness. J. C. Willems [1] used state space model and dissipativeness and defined an inequality involving storage function and the supply function [2] . However, the concept of dissipativeness steered in a different direction. Ilya Prigogine, a well known pioneer in dissipative system was the first to establish the dissipative structure for thermodynamics [3] . Though Prigogine70's fundamental concept is based on chemistry and physics, nevertheless, the theory of dissipative structure is treated as an open system model with extended capacity to transformation and revolutionize based on the same concept, the flow of entity. Progogine70's dissipative structure concluded that the order in an open system can be maintained only in a non-equilibrium condition [4] , which regards to the need of a continuous exchange of energy and resources in an open system with the environment in order to be able to renew itself. However, dissipative system under control theory has developed to be somewhat opposite to thermodynamic70's model; a model in which the dissipative structure is considered to be operated under an equilibrium state. The extensive study of this topic redefines its significants, as mechanical and electrical engineers continue to strive deeper and expand its role in various applications.
Nevertheless, dissipative system is not only being researched under the bases of control theory, other branches of study such as Physics, specifically quantum mechanics, model the system as quantum dissipative systems, where the dissipative systems in quantum mechanics studies the system in microscopic level that focuses on the behaviour of the system [5] .
The system dynamics of dissipative system is quite diverse ever since the model was introduced by Prigogine [4] . The topic of study for dissipative system branches of into a wide range of theories and contexts, the branches of study also becomes narrow, complex, and specific. In terms of control theory, it is important to continue to explore and pursue the impact of dissipative system, for its contribution meets no bound nor a single form which can describes its essence.
Regardless with the perspective of dissipative system, the dissipative property, for which it can be considered when dealing with a flow of any entity, is handle uniquely by the system modeller as desired. Thus, there is no generalization nor a "right way" to model dissipative system, it act as an add-on toward a pre-existing model, which one can only judge the model's efficiency, simplicity and performance without the lost of generality.
The models and considering systems of dissipative systems usually deals with unique cases, such as Yan and Liu [6] or Kaneko et al. [7]. Though the cases are different, modelling and the system dynamics are unlikely to be the same, it is difficult to compare or judge one's model with another. However, even with a similar objective and considering systems, by using different methods and techniques, it would still be difficult to compare. A similar paper back in 2010 [8] , with similar objective and considering system, it is noticeable that their method is basic and lack of conservativeness.
In this paper, we will present a more efficient and compact method, to retain assurance of the system performance in spite of model inaccuracies and changes [9] , and the notion of dissipativity being the boundedness of free supply that can be extracted from a system [1] is expressed in a compact filter form. The advantage of having dissipative properties in such filter form allows us to relate and handle many mechanical and electrical engineering applications [10] , where energy is of interest. Hence, the search for a compacted approach to combine these properties into one useful tool begins here.
2. Continuous-time filter design
Consider the continuous-time singular systems,
PPT Slide
Lager Image
where x ( t )∈ R n is the state vector, y ( t )∈ R m is the measurement vector, z ( t )∈ Rq is the signal to be estimated, w ( t )∈ Rq is the input signal, and E is singular matrix that satisfies the rank ( E )= r n . All systems matrices have appropriate dimensions.
Definition 1. [9 - 10] Linear uncertain system is a convex polytopic set defined as δ ( x ) = A (α) x ,and is said to be quadratically stable in the uncertain box of
PPT Slide
Lager Image
, given that αi = 1 and αi ≥ 0 , if and only if there exists a matrix P = PT >0 such that
PPT Slide
Lager Image
.
By Definition 1, the considering system (1) will now be represented with polytopic uncertainty properties in the following form,
PPT Slide
Lager Image
where i is the natural number index.
The objective is to design the full-order linear filter with polytopic uncertainty, defined in Definition 1, in the following is the state-space realization form
PPT Slide
Lager Image
where Af Rn × n , Bf Rn × m, C f Rp × n , and Df Rp × m . In the following sections, time t will omit the for simplicity.
- 2.1 Robust filtering for continuous-time
By combining the considering system (2) with the full-order linear filter (3), where the augmented state vector is defined as
PPT Slide
Lager Image
= [ xT
PPT Slide
Lager Image
T ] T and the estimation error as
PPT Slide
Lager Image
= z-ẑ , we can obtain the following filtering error singular system as
PPT Slide
Lager Image
Taking the difference between the estimated signal with the linear filter signal, the following solution is obtained:
PPT Slide
Lager Image
By defining some terms as follows
PPT Slide
Lager Image
thus, the filtering error dynamics is expressed by
PPT Slide
Lager Image
Using the Lyapunov stability criterion, a sufficient condition for the stability of an equilibrium system can be determined without solving the state equation. Let V (
PPT Slide
Lager Image
) be a continuously differentiable vector function. The system is asymptotically stable if there exists a positive-definite function V (
PPT Slide
Lager Image
) : P Rn [11] , where
PPT Slide
Lager Image
such that
PPT Slide
Lager Image
PPT Slide
Lager Image
Therefore (8) is said to be negative definite, where
PPT Slide
Lager Image
is the derivative along the trajectory of filtering error dynamics (6).
Definition 2. [12] The filtering error dynamic system expressed by (5) is said to be passive with a dissipation η if
PPT Slide
Lager Image
holds all trajectories with x (0) = 0 and all t ≥ 0 .
To add on the dissipativity characteristics to this linear filter, the descriptor system must satisfy the dissipation performance inequality by Definition 2, expressed as Jp . η is the dissipative constant that controls the overall performance of the system, wT
PPT Slide
Lager Image
is the supply rate of the system, and wTw is the weighting factor of the dissipation rate.
Theorem 1. Consider the continuous-time system (2). For a given positive constant η that is less than 1, if there exist a positive-definite matrix of P and R∈Rnxn matrix satisfying ETR=0 such that
PPT Slide
Lager Image
PPT Slide
Lager Image
is quadratically stable with optimum system performance and dissipativeness. The robust filter variables are defined as the following
PPT Slide
Lager Image
Proof. An auxiliary term introduced to (8) to prevent any rank deficiency, and to maintain the system to be negative definite in the final solution as shown in the following [13]
PPT Slide
Lager Image
PPT Slide
Lager Image
is a constant matrix yet to be determined by the LMI condition simulator, and
PPT Slide
Lager Image
through R that satisfies
PPT Slide
Lager Image
The robust filter variables can be determined from the following inequality,
PPT Slide
Lager Image
With the appropriate manipulation of (12), the following condition is derived,
PPT Slide
Lager Image
PPT Slide
Lager Image
To further derive a sufficient condition, some terms are defined as follows:
PPT Slide
Lager Image
In order to separate out the variables within the same term in Ʌ to construct the LMI, similarity transformation technique [14] will be performed, and the following terms are defined as
PPT Slide
Lager Image
The result of ΓɅΓT < 0 can be changed to the following form
PPT Slide
Lager Image
PPT Slide
Lager Image
Where * in (16) denotes the symmetric terms in the matrix, and for simplicity, terms such as Φ +ΦT are represented as < Φ> . We define the following new variables as
PPT Slide
Lager Image
in order to avoid nonlinear terms within the LMI of (16). Therefore, by substituting in (17), we can obtained the LMI condition of (9).
Example 1. Consider the following continuous-time systems matrices of (2) with polytopic uncertainty of two vertices that satisfies Theorem 1:
PPT Slide
Lager Image
All four possible configurations of the two vertices are implemented for simulation. Using the command “[tmin,aopt]=feasp(LMISYS);” , a feasible solution of the system is obtained, with the feasibility radius saturation at 0.005%. the system variables are as follows:
PPT Slide
Lager Image
By using the transfer matrix identity [15] , the relationship between
PPT Slide
Lager Image
and [ Af, Bf, Cf, Df ] is derived to be
PPT Slide
Lager Image
Therefore,
PPT Slide
Lager Image
Thus, the resulting filter guarantees the system's stability and optimal performance with the given dissipative variable η =0.8 .
3. Discrete-time filter design
Consider the discrete-time singular systems,
PPT Slide
Lager Image
where x ( k )∈ Rn is the state vector, y ( k )∈ Rm is the measurement vector, z ( k )∈ Rp is the signal to be estimated, w ( k )∈ Rq is the input signal, and E is singular matrix that satisfies the rank ( E )= r n . All systems matrices have appropriate dimensions.
By Definition 1, the considering system (18) will now be represented with polytopic uncertainty properties in the following form,
PPT Slide
Lager Image
where i is the natural number index.
With the objective same as continuous-time, the design the full-order linear filter with polytopic uncertainty, defined by Definition 1, is in the following state-space realization form
PPT Slide
Lager Image
where Af Rn× n , Bf Rn×m , Cf Rp×n , Df Rp×m .
- 3.1 Robust filtering for discrete-time
Similar to continuous-time, system (19) and (20) are combined, where the augmented state vector is defined as
PPT Slide
Lager Image
( k ) = [xT ( k )
PPT Slide
Lager Image
T ( k )] T and the estimation error as
PPT Slide
Lager Image
( k )= z ( k )- ( k ) , we can obtain the following filtering error descriptor system as
PPT Slide
Lager Image
and by taking the difference between the estimated signal with the linear filter signal we will obtain the following
PPT Slide
Lager Image
Once again, some terms will be defined as follows
PPT Slide
Lager Image
then the filtering error dynamics is expressed by
PPT Slide
Lager Image
By Lyapunov stability criterion, V (
PPT Slide
Lager Image
) is defined as the infinite sum vector function, where
PPT Slide
Lager Image
such that ΔV < 0 . Therefore
PPT Slide
Lager Image
Similar to continuous-time, definition 2 is extended to discrete-time, the filtering error dynamic system expressed by (23) is said to be passive with a dissipation η if
PPT Slide
Lager Image
( w ( k ) T
PPT Slide
Lager Image
( k )- ηw ( k ) T w ( k ) ≥ 0 holds for all trajectories with x (0)=0 and k ≥ 0 [12] .
To join the dissipative characteristics with the linear filter, the singular system must satisfy the dissipation performance inequality defined by definition 2, and it is expressed below as Jp .
PPT Slide
Lager Image
where η is the dissipation variable that controls the overall performance of the system, w ( k ) T
PPT Slide
Lager Image
is the supply rate of the system, and w ( k ) Tw ( k ) is the weighting factor of the dissipation rate.
Theorem 2. Consider the discrete-time system (19). For a given positive constant η that is less than 1, if there exist a positive-definite matrix of P and R∈R n × n matrix satisfying ETR=0 such that
PPT Slide
Lager Image
PPT Slide
Lager Image
is quadratically stable with optimum system performance and dissipativeness. The robust filter variables are defined as the following
PPT Slide
Lager Image
Proof. Since the approach is different than continuous-time, minor steps are highlighted to show the manipulation difference for both cases. Likewise to continuous-time, an auxiliary term is introduced to prevent any rank deficiency, and to maintain the system negative definite properties in the final solution as shown in the following [13]
PPT Slide
Lager Image
where
PPT Slide
Lager Image
is the constant matrix to be determined and
PPT Slide
Lager Image
can be selected satisfy
PPT Slide
Lager Image
. The linear filter variables can be determined by combining all three properties matrix inequality: filtering error dynamic system, dissipativity, and auxiliary term.
PPT Slide
Lager Image
With the appropriate manipulation and substitution of (29), the matrix inequality can be expressed as
PPT Slide
Lager Image
Note: all function of k is omitted from here on for simplicity.
By reforming (30) in matrix inequality form and apply the first technique, Schur complement [16] , we will obtain the following matrix inequality,
PPT Slide
Lager Image
where (31) is defined as follows
PPT Slide
Lager Image
Since the matrix inequality (31) becomes nonlinear from simply matrix substitutions and expansions, the stability check using LMI solver in Matlab cannot be performed and will lead to an unsolvable matrix inequality if we proceed further. Therefore, a second technique is introduced, called the slack variable denote by H [16 - 17] . To derive the sufficient condition from (30), we will define the following as
PPT Slide
Lager Image
There exist matrix H where H + HT >0 , thus, X 4 + X 4 T >0 , which implies X 4 is invertible.
PPT Slide
Lager Image
The * denotes the symmetric term in a symmetric matrix. For simplicity, terms such as Φ + ΦT are represented as < Φ >
Although the slack variable allows us to proceed on with the matrix manipulation, however, a similar nonlinear result will arise as illustrated above,
PPT Slide
Lager Image
. The LMI solver still cannot be used until a third technique is applied. Similar transformation will be performing on Ʌ to separate out multiple variables in one term.
PPT Slide
Lager Image
The results of ΓΛΓT < 0 ,along with its internal matrices expansion, the final form is obtained as (27).
Example 2 . Consider the following discrete-time system:
PPT Slide
Lager Image
that satisfies Theorem 2. Once again, all possible configurations of the two vertices are used for this simulation. Thus, by using the same Matlab command “[tmin,aopt]=feasp(LMISYS);” , a feasible solution is obtained for this system, with the feasibility radius saturation at 0%. The variables are obtained as follows:
PPT Slide
Lager Image
By using the transfer matrix identity [16] , the relationship between
PPT Slide
Lager Image
and [ Af,Bf,Cf,Df ] is derived to be
PPT Slide
Lager Image
The resulting filter guarantees the system's stability and optimal performance with the given dissipative variable ƞ =0.8.
4. Conclusion
A sufficient condition has been satisfied for this robust dissipative filtering system. Formulated as a compact set of linear matrix inequality, the filtering variable has been determined in continuous and discrete time case. The filter is constructed with properties of robustness and dissipativity. Filtering variables can be easily obtained by using a LMI solver.
Hence, this dissipative filtering technique raises significant improvement and uses for application, in both continuous and discrete time, that deals with the flow of energy due to its compactness and robustness. The dissipative filtering technique guarantees maximum performance of flow rate and stability of the system. Though this may not be the most generic form, however, the proposed method offers multi-purpose and flexibility within its design. Future study, research, and extension may continue along this path for deeper and wider use of dissipativeness.
Acknowledgements
이 논문은 2010년도 정부(교육과학기술부)의 재원으로 한국연구재단의 지원을 받아 수행된 연구임. (No. 2010-0025670)
BIO
진 데 니 (陳 敏 英)
2007년 캐나다 요크 대학 졸업. 2007년~2009년 캐나다 Qauntec Geoscience Inc.근무, 2013년~현재 선문대학교 전자공학과 석사과정
E-mail : chandanny@sunmoon.ac.kr
김 종 해 (金 鍾 海)
1993년 경북대학교 전자공학과 졸업. 1998년 동 대학원 전자공학과 졸업(공박). 1998년~2002년 경북대학교 센서기술연구소 전임연구원. 2000년~2001년 일본 오사카대학 객원연구원. 2010년~2011년 미국 조지아텍 방문연구원. 2002년~현재 선문대학교 전자공학과 교수.
Tel : 041-530-2352
E-mail : kjhae@sunmoon.ac.kr
오 도 창 (吳 道 昌)
1991년 경북대학교 전자공학과 졸업. 1997년 동 대학원 전자공학과 졸업(공박). 1997년 2월~1997년 8월 창원대학교 공과대학 국책교수. 2007년~2008년 미국 Univ. of Florida 방문교수. 1997년 8월~현재 건양대학교 의공학부 교수.
Tel : 041-730-5180
E-mail : docoh@konyang.ac.kr
References
Willems J. C. 2007 "Dissipative Dynamical Systems" European Journal of Control 134 - 151
Willems J. C. 1972 "Dissipative Dynamical Systems Part I: General Theory" The Archive for Rational Mechanics and Analysis 45 (5) 321 - 351    DOI : 10.1007/BF00276493
Prigogine I. , Jantsch E. , Waddington C.H. 1976 Evolution and Consciousness Addison-Wesley London "Order through fluctuation: self-organization and social system" 93 - 133
Nicolis G. , Prigogine I. 1977 Self-organization in nonequilibrium systems Wiley New York
Sieberer L. M. , Huber S. D. , Altman E. , Diehi S. 2013 "Dynamical Critical Phenomena in Driven-Dissipative Systems" Phys. Rev. Lett. arXiv.org 110 195301 -    DOI : 10.1103/PhysRevLett.110.195301
Yan L. , Lin B. 2010 "Stabilization with Optimal Performance for Dissipative Discrete-Time Impulsive Hybrid Systems" Advances in Difference Equations 2010 278240 -    DOI : 10.1186/1687-1847-2010-278240
Kaneko O. , Rapisarda P. , Takaba K. 2005 "Totally dissipative systems" Systems & Control Letters 54 705 - 711    DOI : 10.1016/j.sysconle.2004.11.009
LV L. , Duan G. , Wu A. 2010 "Robust dissipative filtering for continuous-time polytopic uncertain time-delay systems" Journal of Harbin Institute of Technology (New Series) 17 (4)
Dorf Richard C. , Bishop Robert H. 1998 Modern Control Systems 8th Edition Addison Wesley Longman, Inc. 689 - 752
Barb Florin Dan , Ben-tel Aharon , Nemirovski Arkadi 2003 Robust Dissipativity of Interval Uncertain Linear Systems SIAM Journal on Control and Optimization 41 (6) 1661 - 1695    DOI : 10.1137/S0363012901398174
Bao J. , Lee P.L. 2007 Dissipativity and Passivity, Process Control: The Passive Systems Approach, 253 Springer 77 -
Grman Ĺubomír , et al. 2005 "Robust stability conditions for polytopic systems" International Journal of Systems Science et al. 36 (15) 961 - 973    DOI : 10.1080/00207720500389592
Kim J.H. 2012 "Discrete Passive Filtering of Singular Systems with Time-varying Delay" Applied Mathematical Sciences 6 (81) 4047 - 4055
Boyd S. , Ghaoui L. E. , Feron E. , Balakrishnan V. 1994 Linear Matrix Inequalities in System and Control Theory, Studies in Applied Mathematics "Analysis of LDIs: Input/Output properties : Dissipativity" 93 - 94
Brogan William L. 1974 Modern Control Theory Third Edition Prentice-Hall Intl., Inc.
Boyd S. , Vandenberghe L. 2004 "Convex Optimization" Cambridge University Press 131 - 132
de Oliveria M.C. 1999 "A new discrete-time robust stability condition" Systems & Control Letters 37 261 - 265    DOI : 10.1016/S0167-6911(99)00035-3