Copyright IFAC System Identification Santa Barbara, California, USA, 2000
|
|
- Dominick Richards
- 6 years ago
- Views:
Transcription
1 Copyright IFAC System Identification Santa Barbara, California, USA, 2000 ON CONSISTENCY AND MODEL VALIDATION FOR SYSTEMS WITH PARAMETER UNCERTAINTY S. Gugercin*,l, A.C. Antoulas *,1 * Department of Electrical and Computer Engineering Rice University Houston, Texas , USA Abstract: In this note we discuss the issue of consistency and model validation for the identification framework introduced in (Antoulas, 1997) and (Antoulas and Anderson, 1999). Keywords: Identification, parameter uncertainty, model validation, consistency. 1. INTRODUCTION To apply robust control, the existence of a nominal model together with some description of the uncertainty is assumed. Usually, uncertainty in the system is taken to be norm bounded. In the past decade system identification has investigated the problem of obtaining an uncertain family of models, where the uncertainty is norm bounded. One of the well-known approaches for this purpose is set membership identification; this approach characterizes the membership set of the parameters consistent with the data, the model structure, and the perturbation characteristics (see e.g. (M. K. H. Fan et al., 1991), (A.C. Antoulas et al., 1999)). In (A.C. Antoulas et al., 1999), measured inputoutput data with some a priori knowledge on the norm of the perturbation, are assumed. A family of uncertain plants is thus obtained with realparametric uncertainty. The family is given in terms of a linear fraction, the uncertainty being a real contraction. In this note we address the question: What is the consistency problem in this framework? Does there always exist an input function which will invalidate a given model? 1 Partially supported by the NSF Grant DMS For a given model how do the minimal energy input sequences required to invalidate it, look like both in the noise-free and noisy cases. The rest of this note is organized as follows. In Section 2, a summary of earlier results on the identification part is given. Section 3 summarizes the issue of consistency and model validation. We conclude with a numerical experiment in Section AN OVERVIEW OF THE IDENTIFICATION FRAMEWORK The results in this section are mainly based on (Antoulas, 1997), (Antoulas and Anderson, 1999), and (A.C. Antoulas et al., 1999). Let us consider the following input-output data records: &t= ~)t CII{2' t=o, 1,.-.,N generated by some discrete time dynamical system. We model the given input-output data as generated by a linear, time-invariant, discrete time system with an additive perturbation; consequently, the data is decomposable as follows: &~=wt+sjt = + Yt where wt comes from the nominal model (i.e. from the linear, time-invariant system) and &t is the 9t
2 perturbation. We introduce the input matrix U and output matrix Y: [ ] U0 ~1 ' " " Us--1?~1 U2 ' " " Us U :=.... E ]~(m+l) s ~m Urn+l "'" ~ra+s--1 YYl Yl " " " Ys-1 Yl Y~ Y~ N(~+i)x~ [ Y:=.. E Yn+l Yn+s-1 for appropriate indices m, n, s (n is the degree of the denominator, m the degree of the numerator, of the system generating the data, and s = N - n + 1). Define the data matrix Jv~: 34:= EN, r:=m+n+2 The data matrices 2Q and Ad are defined similarly. Assumption: throughout the paper it is assumed that some upper bound e on the 2-norm of the perturbation is a priori known. Thus = f14 + A4 with II.hdll2 <_ Consider a linear system with transfer function G(z) = p(a q(~), where degp(z) = m and degq(z) = n; this system will be described by means of its coefficient vector x = [-p0 - pl... pm qo ql --. q ]r (1) where Pi, qi are the coefficients of the monomiai z i in p(z), q(z) respectively. Given such a system x, its misfit with the data matrix 2~I is defined as follows:,(x, Jq)- IIxTJql]2 Ilxb Notice that as argued in (Antoulas, 1997), the true model x0 which generated the data satisfies #(xo, JQ) < e; we thus define the uncertain model family ~: F= { xe ]~; ~(x'yq)- ''xt2~4''2 - Our goal next is to give a characterization of the uncertain family 3 r. Assume that the data matrix has the following Singular Value Decomposition (SVD): J~ = W2V T : [Wr "'" Wl ] " V T where W E I~ ~x~, and V E R ~x~ are orthogonal matrices. A second assumption we will make is 0" 2 > ~: > 0"1; this assumption corresponds to the knowledge of the system order (e.g. by means of preprocessing of the data). Then x E ~ if and only if there exists a constant c such that x = cx~ and x~ = wl + a2(~2w2 + "'" Jr O~r(~rWr where S-4 i= 1,2,...,r, +-.. _< 1 3. CONSISTENCY AND INVALIDATION The problems of consistency and model validation have been studied extensively, starting with Ljung's book ( L. Ljung., 1987); more recent contributions can be found in (K. Poola et hi., 1994), (R.S. Smith and J.C. Doyle, 1992), (R.S. Smith, 1996). However, in our identification framework, these issues need to be re-examined. The main result of our identification scheme asserts that given a level of misfit e and an input u, all models x satisfying #(x, A)I) < e, are not invalidated by the current input u and the corresponding output y. Furthermore all models x satisfying #(x, M) > e, are in fact invalidated. Hence the issue of consistency in this framework can be phrased as follows: Given a certain model described by the model parameter x, does there exist an input function u for which this parameter value is invalidated?. In other words what is the intersection of ellipses obtained by letting u vary through all possible input sequences of a given length and given energy? How does u look like? Below we will address these issues. As shown in (A.C. Antoulas et hi., 1999) and (S. Gugercin, 2000), when we have noisy output measurements, for a given model ~, ~TjQ can be written as xtj~-~u~(xo, fc) + Jkf,~n(Xo, x) where b/=[u(0) u(1)... u(g- n + m) 1' N" = [n(0) n(1)... n(n) ], (xo,~) E ]~,x8 and :,~(Xo, 5) E I~ (N+1) 8 depend only on Xo and wherel=n-n+m+l ands =N-n+l. It is assumed that the system is initially rest (zero initial conditions.) Therefore, the system 2 is invalidated by this set of measurements if t[// + HI:nil > e[[~[[. To obtain all systems which can be invalidated by using some appropriate input, we search for the points on the boundary of this set. Namely, given a parameter vector, we formulate the following optimization problem: min[lut[ s.t. [[UI:+Af ~ll=e[]21[ (2) Assuming that the noise sequence Af has known characteristics (e.g., Gaussian or fixed upper bound), for the given system 2, the problem in (2) is a nonhomogenous quadratic optimization problem on the ball. The minimization problem
3 2 is examined in details in (A.C. Antoulas et al., 1999) and in (S. Gugercin, 2000). Here we briefly review those results, The problem is solved following the method developed in (D. C. Sorensen, 1997). Thus for each parameter vector with the given noise sequence Af, we compute the minimum energy/.4 required to invalidate it. Therefore for a given noise sequence Af, using inputs of given energy, not all systems can be invalidated. The set of models which can be invalidated depends on e and the energy of the inputs u. For fixed e, this set of invalidatable models decreases as the energy of either the input sequence or the initial conditions decreases. The minimal energy input 5/* which solves (2) is unique. Moreover, the true model x0 can never be invalidated using an input of finite energy. Let 5/~f be the solution for the above minimization problem in the noise-free case. Then u,:r = ~ll:~llv~ 71" 1 where lrl is the largest singular value of : and vl is the corresponding left singular vector. Finally, it can be shown that even for infinite length measurements, the set of plants which cannot be invalidated does not shrink to the true model. In this note, we mainly examine the shape and norm of the minimal energy 5/ which invalidates a given model 2. It is easy to prove that for infinite length measurements Z: T is symmetric, positive definite and Toeplitz. Hence, in the noise-free case for infinite length measurements, the minimal input U. I is given by the eigenvector of a symmetric positive definite toeplitz matrix corresponding to the largest eigenvalue. Let A(z) be the filter designed as follows: A(z) = N - n-~-m Z k=0 u~f(k) z-k, As stated in (J. Makhoul, 1981), the eigenvectors of a symmetric positive definite toeplitz matrix are either symmetric or skew-symmetric with respect to middle. Also, the location of the zeros of the filters whose coefficients are the elements of the eigenvectors corresponding to the minimum or maximum eigenvalues lie on the unit circle. Moreover as S. S. Reddi (1984) proves these zeros are distinct. Hence it follows that in the noisefree case for infinite length measurements the minimal energy input 5/*f is either symmetric or skew-symmetric with respect to middle. Also A(z) above has distinct zeros on unit circle. 4. EXAMPLES We consider a LTI first order dynamical system described by two parameters a and b: a(z) - b z--a The numerical values are assumed to be b = 1, a = 0.5, and number of data points N=100. The noise is additive white Gaussian Yt "~ N(0,0.05); e is chosen to be Minimal Norm Inputs We solve the minimization problem defined in (2) using the tools developed in (D. C. Sorensen, 1997). Experiments were run 100 times and then averaged. Let the solution be denoted by H*. Figure I shows the average minimum energy of H on a logarithmic scale, i.e. log10(115/*]1), required to invalidate a given model represented as a deviation from the true model. On that figure (~ stands for the deviation from the nominal pole and/~ for the deviation from the nominal gain. Notice that as the deviation from true model becomes smaller, the minimum energy ]lh*]l increases, and finally as 2-~xo, Ilu*ll-~cc. 05. ~- - -t. mm nocm O v dare,an from the mac mo~l I..-. : :-..! 15 i 2 a~pha Fig. 1. Minimal energy IlU*ll (log scale) necessary to invalidate given model Figure 2 shows the cross section of the 3 dimensional plot in Figure 1 for 115/11=1. It shows that for any choice of a unit energy input function, all models inside the area shown cannot be invalidated. Moreover, all the models 2, on the boundary of this area need unit energy to be invalidated. Now we examine the shape of the input sequences required to invalidate models &, which are on the boundary in Figure 2. We choose 4 models (see the figure) : &T = [ ] 2 T = [ ] 2[ = [ ] 24 ~ = [ ] beta
4 The ~t o( (amnl~t s Ihat c~,nol be,nvahdate<l by unit energy ' =*! O m,~alidati~ inputs for x2 0,,,,,,,,, o,i 0.~ t xt x : ~0 3O O0 ~0 -O2 X2 X4 o' / invahdaling inputs for x ~ Nobly ~.3 01 ~8 ~DI6 J i I I -04 ~Q t~ca 0 Fig. 2. Cross section of Figure I for ItUll = 1: Set of parameter values which cannot be invalidated For each of these chosen models, we compute the minimal energy input sequences L/* and L/*f required to invalidate them in the noisy and noisefree cases. The resulting plots are shown in Figure 3 for xl, x4 and in Figure 4 for x2, x3. It should be noticed that for each case the noisefree solutions L/~I (i) are given by either ± sin (~ i), or -4-(-1) -i sin(~i). Figure 3 and Figure 4 reveal that the noisy solution is a perturbation around the noise-free one; and when the noise-free solution alternates in sign, the resulting noisy solution has a complex structure. Fig. 4. Input sequences that invalidate x2 (above) and x3 (below). result illustrates the fact that even with infinite length measurements, not all the models can be invalidated. 10( 1 os 1.04 ~1.03 I.o2 ~n ~ U vs Oata L~gZh ~e~a~ ~put* for xl ~'i i : :! i : : i ~o 40 6O ao ice I Oala Leng~ (N) kwldk~at]ng inpots for x4 ] o,s.... i i i/ Fig. 5. Minimal energy IlU~N[I required to invalidate the given model x3 vs N, the length of the data. o -o.5 5. CONCLUSION ~t ( o 100 Fig. 3. Input sequences that invalidate xl (above) and x4 (below). 4.2 Norm of the Input vs Data Length In this example we examine the minimum energy of 5/, as a function of the length of the measurements. Let X3 be the model given in the previous example. We solve the minimization problem (2) for x3 as N varies from N=10,.-., 200. Let L/~N be the solution for a given N. Figure 5 shows IIL/~Nll corresponding to the chosen model x3 as N changes. Figure 5 indicates that even as N goes to infinity, ]IL/~NII does not go to zero; it stays constant with a finite nonzero value. This The set of models which can be invalidated using some measurement, depends on the misfit ( and the energy of the inputs u. For fixed ~, this set of invalidatable models decreases as the energy of either the input sequence or the initial conditions decreases. In the noise-free case for the infinite length measurements the minimal norm input L/~, I which is the eigenvector of a symmetric positive definite Toeplitz matrix is symmetric or skew symmetric with respect to middle. Moreover simulations illustrate that the minimal norm input N* for the noisy problem is a perturbation around ~f. 6. REFERENCES Antoulas, A.C. (1997). On modeling for robust control. Proc. SysId'97 pp
5 A.C.Antoulas and B.D.O. Anderson (1999). On the choice of inputs in identification for robust control, Automatica, 35, M.K.H. Fan, A.L. Tits, and J.C.Doyle (1991). Robustness in presence of mixed parametric uncertanity and unmodeled dynamics. IEEE Trans. Automat. Control, 4: 553:572. S. Gugercin. Model Validation and Consistency (2000). M. S. Thesis, Dept. of ECE, Rice University, Houston, TX. J. Makhoul. On the Eigenvectors of Symmetric Toeplitz Matrices (1981). IEEE Transactions on Acoustics, Speech, ans Signal Processing, Vol. ASSP-29, No.4, pp: L. Ljung (1987). System identification: Theory for the user. Prentice Hall. K. Pools, P. Khargonekar, A. Tikku, J. Krause, and K. Nagpal (1994). A time-domain approach to model validation. IEEE Trans. Automat. Control, TAC-39: S. S. Reddi. Eigenvector Properties of Toeplitz Matrices and Their Application to Spectral Analysis of Time Series (1984). Signal Processing, 7: R.S. Smith and J.C. Doyle (1992). Model validation: a connection between robust control and identification. IEEE Trans. Automat. Control, TAC-3~': R.S. Smith (1996). Continuous-time control model validation using finite experimental data. IEEE Trans. Automat. Control, TAC- 41: D. C. Sorensen (1997). Minimization of a Large- Scale Quadratic Function subject to a Spherical Constraint. SIAM J. Optim., Vol. 7, No. 1, pp S. Gugercin, H. Zhang, and A.C. Antoulas (1999). System identification for robust control. Submitted to Automatics.
ThM06-2. Coprime Factor Based Closed-Loop Model Validation Applied to a Flexible Structure
Proceedings of the 42nd IEEE Conference on Decision and Control Maui, Hawaii USA, December 2003 ThM06-2 Coprime Factor Based Closed-Loop Model Validation Applied to a Flexible Structure Marianne Crowder
More informationTHIS paper deals with robust control in the setup associated
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL 50, NO 10, OCTOBER 2005 1501 Control-Oriented Model Validation and Errors Quantification in the `1 Setup V F Sokolov Abstract A priori information required for
More informationLTI Systems, Additive Noise, and Order Estimation
LTI Systems, Additive oise, and Order Estimation Soosan Beheshti, Munther A. Dahleh Laboratory for Information and Decision Systems Department of Electrical Engineering and Computer Science Massachusetts
More informationUNCERTAINTY MODELING VIA FREQUENCY DOMAIN MODEL VALIDATION
AIAA 99-3959 UNCERTAINTY MODELING VIA FREQUENCY DOMAIN MODEL VALIDATION Martin R. Waszak, * NASA Langley Research Center, Hampton, Virginia Dominick Andrisani II, Purdue University, West Lafayette, Indiana
More informationRICE UNIVERSITY. System Identication for Robust Control. Huipin Zhang. A Thesis Submitted. in Partial Fulfillment of the. Requirements for the Degree
RICE UNIVERSITY System Identication for Robust Control by Huipin Zhang A Thesis Submitted in Partial Fulfillment of the Requirements for the Degree Master of Science Approved, Thesis Committee: Athanasios
More informationSparsity in system identification and data-driven control
1 / 40 Sparsity in system identification and data-driven control Ivan Markovsky This signal is not sparse in the "time domain" 2 / 40 But it is sparse in the "frequency domain" (it is weighted sum of six
More informationI. D. Landau, A. Karimi: A Course on Adaptive Control Adaptive Control. Part 9: Adaptive Control with Multiple Models and Switching
I. D. Landau, A. Karimi: A Course on Adaptive Control - 5 1 Adaptive Control Part 9: Adaptive Control with Multiple Models and Switching I. D. Landau, A. Karimi: A Course on Adaptive Control - 5 2 Outline
More informationStructured singular values for Bernoulli matrices
Int. J. Adv. Appl. Math. and Mech. 44 2017 41 49 ISSN: 2347-2529 IJAAMM Journal homepage: www.ijaamm.com International Journal of Advances in Applied Mathematics and Mechanics Structured singular values
More informationOn Moving Average Parameter Estimation
On Moving Average Parameter Estimation Niclas Sandgren and Petre Stoica Contact information: niclas.sandgren@it.uu.se, tel: +46 8 473392 Abstract Estimation of the autoregressive moving average (ARMA)
More informationRobust Loop Shaping Controller Design for Spectral Models by Quadratic Programming
Robust Loop Shaping Controller Design for Spectral Models by Quadratic Programming Gorka Galdos, Alireza Karimi and Roland Longchamp Abstract A quadratic programming approach is proposed to tune fixed-order
More informationSystem Identification by Nuclear Norm Minimization
Dept. of Information Engineering University of Pisa (Italy) System Identification by Nuclear Norm Minimization eng. Sergio Grammatico grammatico.sergio@gmail.com Class of Identification of Uncertain Systems
More informationA New Subspace Identification Method for Open and Closed Loop Data
A New Subspace Identification Method for Open and Closed Loop Data Magnus Jansson July 2005 IR S3 SB 0524 IFAC World Congress 2005 ROYAL INSTITUTE OF TECHNOLOGY Department of Signals, Sensors & Systems
More informationEIGENFILTERS FOR SIGNAL CANCELLATION. Sunil Bharitkar and Chris Kyriakakis
EIGENFILTERS FOR SIGNAL CANCELLATION Sunil Bharitkar and Chris Kyriakakis Immersive Audio Laboratory University of Southern California Los Angeles. CA 9. USA Phone:+1-13-7- Fax:+1-13-7-51, Email:ckyriak@imsc.edu.edu,bharitka@sipi.usc.edu
More informationAdaptive Channel Modeling for MIMO Wireless Communications
Adaptive Channel Modeling for MIMO Wireless Communications Chengjin Zhang Department of Electrical and Computer Engineering University of California, San Diego San Diego, CA 99- Email: zhangc@ucsdedu Robert
More informationModel reduction via tangential interpolation
Model reduction via tangential interpolation K. Gallivan, A. Vandendorpe and P. Van Dooren May 14, 2002 1 Introduction Although most of the theory presented in this paper holds for both continuous-time
More informationOptimal Polynomial Control for Discrete-Time Systems
1 Optimal Polynomial Control for Discrete-Time Systems Prof Guy Beale Electrical and Computer Engineering Department George Mason University Fairfax, Virginia Correspondence concerning this paper should
More informationRECURSIVE SUBSPACE IDENTIFICATION IN THE LEAST SQUARES FRAMEWORK
RECURSIVE SUBSPACE IDENTIFICATION IN THE LEAST SQUARES FRAMEWORK TRNKA PAVEL AND HAVLENA VLADIMÍR Dept of Control Engineering, Czech Technical University, Technická 2, 166 27 Praha, Czech Republic mail:
More informationROBUST CONSTRAINED REGULATORS FOR UNCERTAIN LINEAR SYSTEMS
ROBUST CONSTRAINED REGULATORS FOR UNCERTAIN LINEAR SYSTEMS Jean-Claude HENNET Eugênio B. CASTELAN Abstract The purpose of this paper is to combine several control requirements in the same regulator design
More informationA Subspace Approach to Estimation of. Measurements 1. Carlos E. Davila. Electrical Engineering Department, Southern Methodist University
EDICS category SP 1 A Subspace Approach to Estimation of Autoregressive Parameters From Noisy Measurements 1 Carlos E Davila Electrical Engineering Department, Southern Methodist University Dallas, Texas
More informationKrylov Techniques for Model Reduction of Second-Order Systems
Krylov Techniques for Model Reduction of Second-Order Systems A Vandendorpe and P Van Dooren February 4, 2004 Abstract The purpose of this paper is to present a Krylov technique for model reduction of
More informationPrediction, filtering and smoothing using LSCR: State estimation algorithms with guaranteed confidence sets
2 5th IEEE Conference on Decision and Control and European Control Conference (CDC-ECC) Orlando, FL, USA, December 2-5, 2 Prediction, filtering and smoothing using LSCR: State estimation algorithms with
More informationSubmitted to Electronics Letters. Indexing terms: Signal Processing, Adaptive Filters. The Combined LMS/F Algorithm Shao-Jen Lim and John G. Harris Co
Submitted to Electronics Letters. Indexing terms: Signal Processing, Adaptive Filters. The Combined LMS/F Algorithm Shao-Jen Lim and John G. Harris Computational Neuro-Engineering Laboratory University
More informationExploring Granger Causality for Time series via Wald Test on Estimated Models with Guaranteed Stability
Exploring Granger Causality for Time series via Wald Test on Estimated Models with Guaranteed Stability Nuntanut Raksasri Jitkomut Songsiri Department of Electrical Engineering, Faculty of Engineering,
More informationRobust Least-Squares Filtering With a Relative Entropy Constraint
Robust Least-Squares Filtering With a Relative Entropy Constraint Bernard C. Levy ECE Department University of California Davis, CA 95616 USA and Abstract Ramine Nikoukhah INRIA Domaine de Voluceau Rocquencourt
More informationThe Solvability Conditions for the Inverse Eigenvalue Problem of Hermitian and Generalized Skew-Hamiltonian Matrices and Its Approximation
The Solvability Conditions for the Inverse Eigenvalue Problem of Hermitian and Generalized Skew-Hamiltonian Matrices and Its Approximation Zheng-jian Bai Abstract In this paper, we first consider the inverse
More informationIterative Learning Control Analysis and Design I
Iterative Learning Control Analysis and Design I Electronics and Computer Science University of Southampton Southampton, SO17 1BJ, UK etar@ecs.soton.ac.uk http://www.ecs.soton.ac.uk/ Contents Basics Representations
More informationMathematics for Control Theory
Mathematics for Control Theory H 2 and H system norms Hanz Richter Mechanical Engineering Department Cleveland State University Reading materials We will use: Michael Green and David Limebeer, Linear Robust
More informationDESIGNING A KALMAN FILTER WHEN NO NOISE COVARIANCE INFORMATION IS AVAILABLE. Robert Bos,1 Xavier Bombois Paul M. J. Van den Hof
DESIGNING A KALMAN FILTER WHEN NO NOISE COVARIANCE INFORMATION IS AVAILABLE Robert Bos,1 Xavier Bombois Paul M. J. Van den Hof Delft Center for Systems and Control, Delft University of Technology, Mekelweg
More informationApproximate Bisimulations for Constrained Linear Systems
Approximate Bisimulations for Constrained Linear Systems Antoine Girard and George J Pappas Abstract In this paper, inspired by exact notions of bisimulation equivalence for discrete-event and continuous-time
More informationSubspace Identification
Chapter 10 Subspace Identification Given observations of m 1 input signals, and p 1 signals resulting from those when fed into a dynamical system under study, can we estimate the internal dynamics regulating
More informationNOISE ROBUST RELATIVE TRANSFER FUNCTION ESTIMATION. M. Schwab, P. Noll, and T. Sikora. Technical University Berlin, Germany Communication System Group
NOISE ROBUST RELATIVE TRANSFER FUNCTION ESTIMATION M. Schwab, P. Noll, and T. Sikora Technical University Berlin, Germany Communication System Group Einsteinufer 17, 1557 Berlin (Germany) {schwab noll
More informationPartial Eigenvalue Assignment in Linear Systems: Existence, Uniqueness and Numerical Solution
Partial Eigenvalue Assignment in Linear Systems: Existence, Uniqueness and Numerical Solution Biswa N. Datta, IEEE Fellow Department of Mathematics Northern Illinois University DeKalb, IL, 60115 USA e-mail:
More informationH-Infinity Controller Design for a Continuous Stirred Tank Reactor
International Journal of Electronic and Electrical Engineering. ISSN 974-2174 Volume 7, Number 8 (214), pp. 767-772 International Research Publication House http://www.irphouse.com H-Infinity Controller
More informationPOLYNOMIAL SINGULAR VALUES FOR NUMBER OF WIDEBAND SOURCES ESTIMATION AND PRINCIPAL COMPONENT ANALYSIS
POLYNOMIAL SINGULAR VALUES FOR NUMBER OF WIDEBAND SOURCES ESTIMATION AND PRINCIPAL COMPONENT ANALYSIS Russell H. Lambert RF and Advanced Mixed Signal Unit Broadcom Pasadena, CA USA russ@broadcom.com Marcel
More informationAuxiliary signal design for failure detection in uncertain systems
Auxiliary signal design for failure detection in uncertain systems R. Nikoukhah, S. L. Campbell and F. Delebecque Abstract An auxiliary signal is an input signal that enhances the identifiability of a
More informationObservability and state estimation
EE263 Autumn 2015 S Boyd and S Lall Observability and state estimation state estimation discrete-time observability observability controllability duality observers for noiseless case continuous-time observability
More informationNetwork Parameters of a Two-Port Filter
Contents Network Parameters of a Two-Port Filter S Parameters in the s-domain Relation between ε and ε R 2 2 Derivation of S 22 3 2 Network Parameters in the jω-domain 3 2 S Parameters in the jω-domain
More informationChapter 9 Observers, Model-based Controllers 9. Introduction In here we deal with the general case where only a subset of the states, or linear combin
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A. Dahleh George Verghese Department of Electrical Engineering and Computer Science Massachuasetts Institute of Technology c Chapter 9 Observers,
More informationMulti-Model Adaptive Regulation for a Family of Systems Containing Different Zero Structures
Preprints of the 19th World Congress The International Federation of Automatic Control Multi-Model Adaptive Regulation for a Family of Systems Containing Different Zero Structures Eric Peterson Harry G.
More informationMarcus Pantoja da Silva 1 and Celso Pascoli Bottura 2. Abstract: Nonlinear systems with time-varying uncertainties
A NEW PROPOSAL FOR H NORM CHARACTERIZATION AND THE OPTIMAL H CONTROL OF NONLINEAR SSTEMS WITH TIME-VARING UNCERTAINTIES WITH KNOWN NORM BOUND AND EXOGENOUS DISTURBANCES Marcus Pantoja da Silva 1 and Celso
More informationIyad T. Abu-Jeib and Thomas S. Shores
NEW ZEALAND JOURNAL OF MATHEMATICS Volume 3 (003, 9 ON PROPERTIES OF MATRIX I ( OF SINC METHODS Iyad T. Abu-Jeib and Thomas S. Shores (Received February 00 Abstract. In this paper, we study the determinant
More informationAn Exact Stability Analysis Test for Single-Parameter. Polynomially-Dependent Linear Systems
An Exact Stability Analysis Test for Single-Parameter Polynomially-Dependent Linear Systems P. Tsiotras and P.-A. Bliman Abstract We provide a new condition for testing the stability of a single-parameter,
More informationOn Solving Large Algebraic. Riccati Matrix Equations
International Mathematical Forum, 5, 2010, no. 33, 1637-1644 On Solving Large Algebraic Riccati Matrix Equations Amer Kaabi Department of Basic Science Khoramshahr Marine Science and Technology University
More informationHybrid Systems - Lecture n. 3 Lyapunov stability
OUTLINE Focus: stability of equilibrium point Hybrid Systems - Lecture n. 3 Lyapunov stability Maria Prandini DEI - Politecnico di Milano E-mail: prandini@elet.polimi.it continuous systems decribed by
More informationFundamental Limitations in Self-Sensing Magnetic Bearings when Modeled as Linear Periodic Systems
Fundamental Limitations in Self-Sensing Magnetic Bearings when Modeled as Linear Periodic Systems K.S. Peterson and R.H. Middleton and J.S. Freudenberg Abstract In Magnetic Bearing Measurement Configurations
More informationAn efficient algorithm to compute the real perturbation values of a matrix
An efficient algorithm to compute the real perturbation values of a matrix Simon Lam and Edward J. Davison 1 Abstract In this paper, an efficient algorithm is presented for solving the nonlinear 1-D optimization
More informationFunctional Analysis Review
Outline 9.520: Statistical Learning Theory and Applications February 8, 2010 Outline 1 2 3 4 Vector Space Outline A vector space is a set V with binary operations +: V V V and : R V V such that for all
More informationLecture 19 Observability and state estimation
EE263 Autumn 2007-08 Stephen Boyd Lecture 19 Observability and state estimation state estimation discrete-time observability observability controllability duality observers for noiseless case continuous-time
More informationGLOBAL ANALYSIS OF PIECEWISE LINEAR SYSTEMS USING IMPACT MAPS AND QUADRATIC SURFACE LYAPUNOV FUNCTIONS
GLOBAL ANALYSIS OF PIECEWISE LINEAR SYSTEMS USING IMPACT MAPS AND QUADRATIC SURFACE LYAPUNOV FUNCTIONS Jorge M. Gonçalves, Alexandre Megretski y, Munther A. Dahleh y California Institute of Technology
More informationAnalysis of Error Produced by Truncated SVD and Tikhonov Regularization Met hods *
Analysis of Error Produced by Truncated SVD and Tikhonov Regularization Met hods * Irina F. Gorodnitsky and Bhaskar D. Rao Dept. of Electrical and Computer Engineering University of California, San Diego
More informationTHE NUMERICAL EVALUATION OF THE MAXIMUM-LIKELIHOOD ESTIMATE OF A SUBSET OF MIXTURE PROPORTIONS*
SIAM J APPL MATH Vol 35, No 3, November 1978 1978 Society for Industrial and Applied Mathematics 0036-1399/78/3503-0002 $0100/0 THE NUMERICAL EVALUATION OF THE MAXIMUM-LIKELIHOOD ESTIMATE OF A SUBSET OF
More informationChapter 3. LQ, LQG and Control System Design. Dutch Institute of Systems and Control
Chapter 3 LQ, LQG and Control System H 2 Design Overview LQ optimization state feedback LQG optimization output feedback H 2 optimization non-stochastic version of LQG Application to feedback system design
More informationConvexity of the Joint Numerical Range
Convexity of the Joint Numerical Range Chi-Kwong Li and Yiu-Tung Poon October 26, 2004 Dedicated to Professor Yik-Hoi Au-Yeung on the occasion of his retirement. Abstract Let A = (A 1,..., A m ) be an
More informationH-INFINITY CONTROLLER DESIGN FOR A DC MOTOR MODEL WITH UNCERTAIN PARAMETERS
Engineering MECHANICS, Vol. 18, 211, No. 5/6, p. 271 279 271 H-INFINITY CONTROLLER DESIGN FOR A DC MOTOR MODEL WITH UNCERTAIN PARAMETERS Lukáš Březina*, Tomáš Březina** The proposed article deals with
More informationA NUMERICAL METHOD TO SOLVE A QUADRATIC CONSTRAINED MAXIMIZATION
A NUMERICAL METHOD TO SOLVE A QUADRATIC CONSTRAINED MAXIMIZATION ALIREZA ESNA ASHARI, RAMINE NIKOUKHAH, AND STEPHEN L. CAMPBELL Abstract. The problem of maximizing a quadratic function subject to an ellipsoidal
More informationIN THIS paper, we consider optimal control of vehicular platoons.
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 50, NO. 9, SEPTEMBER 2005 1307 On the Ill-Posedness of Certain Vehicular Platoon Control Problems Mihailo R. Jovanović, Member, IEEE, and Bassam Bamieh, Senior
More informationCramér-Rao Bounds for Estimation of Linear System Noise Covariances
Journal of Mechanical Engineering and Automation (): 6- DOI: 593/jjmea Cramér-Rao Bounds for Estimation of Linear System oise Covariances Peter Matiso * Vladimír Havlena Czech echnical University in Prague
More informationLecture 7 : Generalized Plant and LFT form Dr.-Ing. Sudchai Boonto Assistant Professor
Dr.-Ing. Sudchai Boonto Assistant Professor Department of Control System and Instrumentation Engineering King Mongkuts Unniversity of Technology Thonburi Thailand Linear Quadratic Gaussian The state space
More informationGeneralized Pencil Of Function Method
Generalized Pencil Of Function Method In recent times, the methodology of approximating a function by a sum of complex exponentials has found applications in many areas of electromagnetics For example
More informationTime-delay in high-gain observer based disturbance estimation
Time-delay in high-gain observer based disturbance estimation Xuewu Dai, Zhiwei Gao Usama Abou-Zayed Tim Breikin Control Systems Centre, The University of Manchester, Manchester UK, M6 1QD. Southwest University,
More informationAPPROXIMATE SOLUTION OF A SYSTEM OF LINEAR EQUATIONS WITH RANDOM PERTURBATIONS
APPROXIMATE SOLUTION OF A SYSTEM OF LINEAR EQUATIONS WITH RANDOM PERTURBATIONS P. Date paresh.date@brunel.ac.uk Center for Analysis of Risk and Optimisation Modelling Applications, Department of Mathematical
More informationAn Adaptive Sensor Array Using an Affine Combination of Two Filters
An Adaptive Sensor Array Using an Affine Combination of Two Filters Tõnu Trump Tallinn University of Technology Department of Radio and Telecommunication Engineering Ehitajate tee 5, 19086 Tallinn Estonia
More informationMULTIVARIABLE ZEROS OF STATE-SPACE SYSTEMS
Copyright F.L. Lewis All rights reserved Updated: Monday, September 9, 8 MULIVARIABLE ZEROS OF SAE-SPACE SYSEMS If a system has more than one input or output, it is called multi-input/multi-output (MIMO)
More informationRobust Stabilization of the Uncertain Linear Systems. Based on Descriptor Form Representation t. Toru ASAI* and Shinji HARA**
Robust Stabilization of the Uncertain Linear Systems Based on Descriptor Form Representation t Toru ASAI* and Shinji HARA** This paper proposes a necessary and sufficient condition for the quadratic stabilization
More informationControl Systems. Laplace domain analysis
Control Systems Laplace domain analysis L. Lanari outline introduce the Laplace unilateral transform define its properties show its advantages in turning ODEs to algebraic equations define an Input/Output
More informationTitle without the persistently exciting c. works must be obtained from the IEE
Title Exact convergence analysis of adapt without the persistently exciting c Author(s) Sakai, H; Yang, JM; Oka, T Citation IEEE TRANSACTIONS ON SIGNAL 55(5): 2077-2083 PROCESS Issue Date 2007-05 URL http://hdl.handle.net/2433/50544
More informationThe Essentials of Linear State-Space Systems
:or-' The Essentials of Linear State-Space Systems J. Dwight Aplevich GIFT OF THE ASIA FOUNDATION NOT FOR RE-SALE John Wiley & Sons, Inc New York Chichester Weinheim OAI HOC OUOC GIA HA N^l TRUNGTAMTHANCTINTHUVIIN
More informationTHIS paper studies the input design problem in system identification.
1534 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 50, NO. 10, OCTOBER 2005 Input Design Via LMIs Admitting Frequency-Wise Model Specifications in Confidence Regions Henrik Jansson Håkan Hjalmarsson, Member,
More informationStructured Stochastic Uncertainty
Structured Stochastic Uncertainty Bassam Bamieh Department of Mechanical Engineering University of California at Santa Barbara Santa Barbara, CA, 9306 bamieh@engineeringucsbedu Abstract We consider linear
More informationRobust fixed-order H Controller Design for Spectral Models by Convex Optimization
Robust fixed-order H Controller Design for Spectral Models by Convex Optimization Alireza Karimi, Gorka Galdos and Roland Longchamp Abstract A new approach for robust fixed-order H controller design by
More informationFloor Control (kn) Time (sec) Floor 5. Displacement (mm) Time (sec) Floor 5.
DECENTRALIZED ROBUST H CONTROL OF MECHANICAL STRUCTURES. Introduction L. Bakule and J. Böhm Institute of Information Theory and Automation Academy of Sciences of the Czech Republic The results contributed
More informationDynamic measurement: application of system identification in metrology
1 / 25 Dynamic measurement: application of system identification in metrology Ivan Markovsky Dynamic measurement takes into account the dynamical properties of the sensor 2 / 25 model of sensor as dynamical
More informationNOTICE WARNING CONCERNING COPYRIGHT RESTRICTIONS: The copyright law of the United States (title 17, U.S. Code) governs the making of photocopies or
NOTICE WARNING CONCERNING COPYRIGHT RESTRICTIONS: The copyright law of the United States (title 17, U.S. Code) governs the making of photocopies or other reproductions of copyrighted material. Any copying
More informationFundamental Limitations in Self-Sensing Magnetic Bearings when Modeled as Linear Periodic Systems
Proceedings of the 6 American Control Conference Minneapolis, Minnesota, USA, June 14-16, 6 FrA11.1 Fundamental Limitations in Self-Sensing Magnetic Bearings when Modeled as Linear Periodic Systems K.S.
More informationH 2 optimal model reduction - Wilson s conditions for the cross-gramian
H 2 optimal model reduction - Wilson s conditions for the cross-gramian Ha Binh Minh a, Carles Batlle b a School of Applied Mathematics and Informatics, Hanoi University of Science and Technology, Dai
More informationRobust Observer for Uncertain T S model of a Synchronous Machine
Recent Advances in Circuits Communications Signal Processing Robust Observer for Uncertain T S model of a Synchronous Machine OUAALINE Najat ELALAMI Noureddine Laboratory of Automation Computer Engineering
More informationEFFECTS OF ILL-CONDITIONED DATA ON LEAST SQUARES ADAPTIVE FILTERS. Gary A. Ybarra and S.T. Alexander
EFFECTS OF ILL-CONDITIONED DATA ON LEAST SQUARES ADAPTIVE FILTERS Gary A. Ybarra and S.T. Alexander Center for Communications and Signal Processing Electrical and Computer Engineering Department North
More informationReturn Difference Function and Closed-Loop Roots Single-Input/Single-Output Control Systems
Spectral Properties of Linear- Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 2018! Stability margins of single-input/singleoutput (SISO) systems! Characterizations
More informationSet-based adaptive estimation for a class of nonlinear systems with time-varying parameters
Preprints of the 8th IFAC Symposium on Advanced Control of Chemical Processes The International Federation of Automatic Control Furama Riverfront, Singapore, July -3, Set-based adaptive estimation for
More informationUpper and Lower Bounds of Frequency Interval Gramians for a Class of Perturbed Linear Systems Shaker, Hamid Reza
Aalborg Universitet Upper and Lower Bounds of Frequency Interval Gramians for a Class of Perturbed Linear Systems Shaker, Hamid Reza Published in: 7th IFAC Symposium on Robust Control Design DOI (link
More informationRiccati difference equations to non linear extended Kalman filter constraints
International Journal of Scientific & Engineering Research Volume 3, Issue 12, December-2012 1 Riccati difference equations to non linear extended Kalman filter constraints Abstract Elizabeth.S 1 & Jothilakshmi.R
More informationResearch Article An Equivalent LMI Representation of Bounded Real Lemma for Continuous-Time Systems
Hindawi Publishing Corporation Journal of Inequalities and Applications Volume 28, Article ID 67295, 8 pages doi:1.1155/28/67295 Research Article An Equivalent LMI Representation of Bounded Real Lemma
More informationCONTROL SYSTEMS, ROBOTICS AND AUTOMATION - Vol. XIII - Nonlinear Observers - A. J. Krener
NONLINEAR OBSERVERS A. J. Krener University of California, Davis, CA, USA Keywords: nonlinear observer, state estimation, nonlinear filtering, observability, high gain observers, minimum energy estimation,
More informationFixed Order H Controller for Quarter Car Active Suspension System
Fixed Order H Controller for Quarter Car Active Suspension System B. Erol, A. Delibaşı Abstract This paper presents an LMI based fixed-order controller design for quarter car active suspension system in
More informationComplex Variables...Review Problems (Residue Calculus Comments)...Fall Initial Draft
Complex Variables........Review Problems Residue Calculus Comments)........Fall 22 Initial Draft ) Show that the singular point of fz) is a pole; determine its order m and its residue B: a) e 2z )/z 4,
More informationTakagi Sugeno Fuzzy Sliding Mode Controller Design for a Class of Nonlinear System
Australian Journal of Basic and Applied Sciences, 7(7): 395-400, 2013 ISSN 1991-8178 Takagi Sugeno Fuzzy Sliding Mode Controller Design for a Class of Nonlinear System 1 Budiman Azzali Basir, 2 Mohammad
More informationADAPTIVE FILTER THEORY
ADAPTIVE FILTER THEORY Fourth Edition Simon Haykin Communications Research Laboratory McMaster University Hamilton, Ontario, Canada Front ice Hall PRENTICE HALL Upper Saddle River, New Jersey 07458 Preface
More informationDidier HENRION henrion
POLYNOMIAL METHODS FOR ROBUST CONTROL Didier HENRION www.laas.fr/ henrion henrion@laas.fr Laboratoire d Analyse et d Architecture des Systèmes Centre National de la Recherche Scientifique Université de
More informationProjection of state space realizations
Chapter 1 Projection of state space realizations Antoine Vandendorpe and Paul Van Dooren Department of Mathematical Engineering Université catholique de Louvain B-1348 Louvain-la-Neuve Belgium 1.0.1 Description
More informationFrequency-Domain Robust Control Toolbox
Frequency-Domain Robust Control Toolbox Alireza Karimi Abstract A new frequency-domain robust control toolbox is introduced and compared with some features of the robust control toolbox of Matlab. A summary
More informationLQR, Kalman Filter, and LQG. Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin
LQR, Kalman Filter, and LQG Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin May 2015 Linear Quadratic Regulator (LQR) Consider a linear system
More informationLyapunov Stability of Linear Predictor Feedback for Distributed Input Delays
IEEE TRANSACTIONS ON AUTOMATIC CONTROL VOL. 56 NO. 3 MARCH 2011 655 Lyapunov Stability of Linear Predictor Feedback for Distributed Input Delays Nikolaos Bekiaris-Liberis Miroslav Krstic In this case system
More informationDESIGN OF ROBUST CONTROL SYSTEM FOR THE PMS MOTOR
Journal of ELECTRICAL ENGINEERING, VOL 58, NO 6, 2007, 326 333 DESIGN OF ROBUST CONTROL SYSTEM FOR THE PMS MOTOR Ahmed Azaiz Youcef Ramdani Abdelkader Meroufel The field orientation control (FOC) consists
More informationc 21 w 2 c 22 s 2 c 12
Blind Adaptive Cross-Pole Interference Cancellation Using Fractionally-Spaced CMA Wonzoo Chung, John Treichler y, and C. Richard Johnson, Jr. wonzoo(johnson)@ee.cornell.edu y jrt@appsig.com School of Elec.
More informationSubspace-based Identification
of Infinite-dimensional Multivariable Systems from Frequency-response Data Department of Electrical and Electronics Engineering Anadolu University, Eskişehir, Turkey October 12, 2008 Outline 1 2 3 4 Noise-free
More informationIS NEGATIVE STEP SIZE LMS ALGORITHM STABLE OPERATION POSSIBLE?
IS NEGATIVE STEP SIZE LMS ALGORITHM STABLE OPERATION POSSIBLE? Dariusz Bismor Institute of Automatic Control, Silesian University of Technology, ul. Akademicka 16, 44-100 Gliwice, Poland, e-mail: Dariusz.Bismor@polsl.pl
More informationCover page. : On-line damage identication using model based orthonormal. functions. Author : Raymond A. de Callafon
Cover page Title : On-line damage identication using model based orthonormal functions Author : Raymond A. de Callafon ABSTRACT In this paper, a new on-line damage identication method is proposed for monitoring
More informationThe norms can also be characterized in terms of Riccati inequalities.
9 Analysis of stability and H norms Consider the causal, linear, time-invariant system ẋ(t = Ax(t + Bu(t y(t = Cx(t Denote the transfer function G(s := C (si A 1 B. Theorem 85 The following statements
More informationSubject: Optimal Control Assignment-1 (Related to Lecture notes 1-10)
Subject: Optimal Control Assignment- (Related to Lecture notes -). Design a oil mug, shown in fig., to hold as much oil possible. The height and radius of the mug should not be more than 6cm. The mug must
More informationFundamentals of Matrices
Maschinelles Lernen II Fundamentals of Matrices Christoph Sawade/Niels Landwehr/Blaine Nelson Tobias Scheffer Matrix Examples Recap: Data Linear Model: f i x = w i T x Let X = x x n be the data matrix
More information