APPROXIMATE SOLUTION OF A SYSTEM OF LINEAR EQUATIONS WITH RANDOM PERTURBATIONS
|
|
- Ursula Weaver
- 5 years ago
- Views:
Transcription
1 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 Sciences, Brunel University, U.K. Abstract This work suggests a way of finding an approximate solution to a system of linear equations of the form AX = b, A = A + with a known square matrix A, a known vector b and a small structured random perturbation. Under certain realistic assumptions about the structure and the smallness of perturbations, it is shown that a uniformly convergent sequence of approximations to the random variable X = A b can be obtained in terms of polynomials of the random perturbations. Possible applications of this method in control theory are suggested. Numerical examples demonstrate the applicability of this method. 1 Introduction Many estimation and control problems rely on solution of linear least-squares problems or uncertain matrix inversion problems of the form AX = b. When the nominal data {A, b} are subject to uncertainty and/or disturbances, the performance of the optimal estimator may degrade appreciably. Similar problems arise while computing frequency response of a system with uncertain poles. Numerous approaches to alleviate this problem have been suggested. Min-max approach to least-squares problem under bounded data uncertainties has been studied in [1] and []. This approach formulates and solves problems of the form min X max AX b A A, b B where A and B are appropriate bounded uncertainty sets. Data-fitting under bounded uncertainties has corresponding author been considered in a more general setting in [3]. Applications of min-max based least-squares solutions have been reported in state estimation in general dynamical systems [4] and in water networks [5]. If statistical information is available about the uncertainty, the min-max based approach described above may be too conservative for use in certain applications. If the joint probability distribution function of all the uncertain parameters is known, the joint probability distribution of A b may be completely characterised using results in [6]. However, this level of information is rarely available. In this paper, the focus is to obtain computationally simpler estimates with a good average-case accuracy possibly at the expense of worst-case accuracy using substantially less information than the full probability density. To this end, this paper considers a sequence of uniform approximations to the true solution which is a random variable of the system of equations. The statistics of approximants is much simpler to compute than that of the true solution. The rest of the paper is organised as follows. The next section gives the exact assumptions and the problem formulation for the case when A is a square matrix and is invertible almost surely. Section 3 discusses the main results of this paper. Section 4 describes the case when the random perturbations have a special affine structure. Section 5 describes applications to analysis of uncertainty in frequency response of linear systems. Section 6 demonstrates these results with numerical examples. Problem formulation Notation in the paper is standard. Let R m n respectively, C m n denote the space of m n real
2 resp. complex matrices. Let R n denote the space of n vectors. I denotes identity matrix; its size is determined by the context. Ai, j denotes i, j th element of matrix A. Vectors will be denoted by boldface characters. In the problems considered in this paper, a system of linear equations of the following form is assumed to be given: AX = b, A = A +, A Q m m, b R m 1 where Q = R or C depending on context and is a real matrix-valued random variable satisfying P A < 1 = 1. Here denotes the maximum singular value and PΦ denotes the probability of occurence of event Φ. It is seen that the random perturbation is assumed to be small relative to the nominal value in a welldefined sense. Specific instances of the uncertainty set will be considered in section 4. 3 Main result Theorem 1 Suppose the given system of equations 1 satisfies and suppose that the inverse of A exists with probability 1. Define X = A b, X N = N A i A 3 with A = I. Then X N X with probability 1. Proof : The proof rests on the following standard result in probability see, e.g., [7], theorem 7.4: If N P X N X > ɛ < for all ɛ >, then X N X with probability 1. The proof is based on deriving an upper bound on P X N X > ɛ and then showing that summation of these upper bounds over N is finite. 1 Since P A < 1 = 1, the following power series expansion holds with probability 1 see, e.g. [8], chapter 5: A + = I + A = A A i A 1 From proof of [7], theorem 7.4, it is easy to see that analogous result holds for vector-valued random variables. so that P X = A i A = 1. Suppose that P A α = 1 for some α < 1. Also, let A = β. For a given ɛ >, let N ɛ be the smallest integer such that αnɛ+1 β 1 α ɛ. Now, using definition of X N, XNɛ > ɛ { = P i=n ɛ A i } A > ɛ =, since, with probability 1, 4 { } A i A i=n ɛ Next, for any ɛ >, N= N ɛ XN > ɛ N= since N ɛ completes the proof. αnɛ+1 β 1 α ɛ. 5 XN > ɛ <, 6 log ɛ1 α αβ logα is finite for any ɛ >. This The above result shows a way of obtaining a series of random vectors which converges with probability 1 to the true solution of the system of equations. The expected value of X N is given by N EX N = E A i A. This is a function of the first N moments of random variable and may thus be found using a limited information about the uncertainty. In the next section, with a special structure is considered, which yields particularly simple low order approximations. 4 Affinely parameterised perturbations Consider the system 1 with which satisfies and has the following structure: = a i A i, A i R m m, 7
3 and a i are scalar random variables with Ea i =, i = 1,,..., k. The matrices A i are in general, sparse matrices, with 1 s in places where the random parameter a i has an impact on the nominal entry in the matrix A and zeros at all other elements numerical examples will make this point clear. The assumption that a i are zero-mean random variables is perfectly reasonable and implies that the mean is absorbed in A. It is possible to obtain useful low order approximations may be obtained with a very small number of terms for this structure, as explained below. 1. The first order approximation X 1 = I a i A A i A is a linear combination of zero mean random variables a i and is of independent interest. If the variance of a i are known and if N is sufficiently large, this expression may be used alongwith central limit theorem to build approximate confidence intervals around each element of the nominal th order solution A. If N is small but the upper and the lower bounds on a i are known, Hoeffding s inequality [9] may be used to obtain potentially conservative approximate error bounds. See [1] for a recent application of Hoeffding s inequalities in engineering. Alternatively, X 1 may simply be used to test the sensitivity of the solution of a system of linear equations to small perturbations in certain entries in the matrix without having to solve the perturbed system.. If a i are uncorrelated, it is easy to show that EX = I + Ea i A A i A. Note that the computation for EX does not require sampling the distribution of a i. For a large number of uncertain parameters, the cost-saving to build an accurate estimate of the average value of A b may be significant. Simulation experiments indicate that higher order approximations indeed tend to be more accurate than the nominal th order solution A. This is demonstrated in section 6 through an example. If first 4 moments of a i are known, it is possible to obtain covariance matrix of X. The expression for covariance is straighforward if tedious and is ommitted. It is necessary here to comment on the relationship between the order of approximation and the size of perturbations. The bound on the relative size of perturbation A will not be known in general. However, this bound does not appear in X N itself. If X N for some N 1 differs substantially by more than 1%, say from X = A, it may indicate that the small perturbations condition may be violated and the data is not reliable enough to build an estimate. In general, it may be seen from the proof of theorem 1 that X N X X A N+1 1 α holds with probability 1. The choice of N then depends on the trade-off between tractability and a priori knowledge of the size of uncertainty. Finally, the results presented here may be trivially extended to the case when b = b + δ b for a random perturbation δ b which is uncorrelated to. 5 Applications The technique presented in this paper is generic and may have a wide variety of applications. Two applications are described in some detail here. The focus of the first application is to map the available information about uncertainty in parameters of a linear system into corresponding uncertainty in its frequency response. The second application relates to perturbation of a transfer function of a controller and its effect on the closed-loop transfer functions. 5.1 Uncertainty in frequency response of a linear system due to parametric perturbations If the parametric uncertainty in a model of a linear system is in the form of a covariance matrix of parameter estimates, there are numerous other ways of mapping it into frequency response; see [11] and references therein. However, the parameters may be obtained by physical knowledge of the underlying dynamics, e.g. by knowing the values of nomial motor resistance and inductance in a DC drive. The motivation of this application is the latter situation, when the parameters may be known within certain tolerances around their nominal values.
4 Consider a continuous time linear, shift-invariant system in a standard notation, dx = Axt + But, dt y k = Cxt + Dut where A = A +, A is a constant real matrix and is a matrix-valued real random variable satisfying P jωi A < 1 = 1, ω. 8 The perturbation may represent the uncertainty in physical parameters which yield the poles of the system. Following the steps of theorem 1, the first order approximation of the uncertain frequency response P jω = CjωI A B + D at a frequency ω is given by P 1 = CjωI A B + D CjωI A jωi A B. 9 Unlike P, P 1 is affine in uncertainty. If the distribution of this uncertainty is known, approximate confidence intervals of pointwise frequency response of the system may be easily obtained from Uncertainty in closed-loop frequency response due to controller perturbations Consider a closed-loop described by the following equations: y = P u + w u = Cy. 1 Here, P is the plant, C is the controller, w is a disturbance acting on the output and u and y are the plant input and plant output respectively. The plant and the controller are assumed to be matrix valued linear operators while the signals are assumed to be vector-valued. It may be easily shown that the transfer function from w to [ y T u ] T T is given by [ [ ] y I P C = u] CI P C w. 11 Suppose the controller is given by C = C + where C is a known transfer function and is a zero mean random perturbation e.g. accounting for tolerances of nominal values of electrical or mechanical components used in hardware implementation. It is of interest to find the effect of this perturbation on the frequency response of the transfer function matrix. For notational simplicity, let S = I P C and let I P C I P C T = CI P C, T = C I P C. It is assumed that P S P jω < 1 = 1 holds at each frequency ω of interest. Using the first order approximation for I P C and then retaining only the first order terms in in the resulting expression yields after some elementary manipulation [ ] S T T + P S I C P. 1 S This expression may be used to find pointwise confidence intervals for the frequency response of the transfer matrix to perturbations in the controller transfer function. 6 Numerical Examples Results of some numerical experiments related to better average case accuracy of higher order approximations are presented here. Specific applications, such as the ones described in the last section, are not discussed due to space constraints. Consider a simple linear system of equations A + X = b with, [ ] [ 1 1 A =, b =, = a 3 3 5] 1 A 1 + a A, 1 A 1 =, A = 1 and a i P. a 1 is assumed to be uniformly distributed in [.4, +.4] and a is assumed to be uniformly distributed in [.6, +.6]. This gives a variation of ±% around nominal values of the corresponding two entries in A. 1 samples of each a 1, a are generated and expected value of solution E X = [ ] T is computed. The nominal solution in this case is A = [ ] T. Using the result described earlier, the expected value of second order approximation is computed as EX = A + Ea 1 A A 1 + Ea A A A = [ ] T.
5 Note that first order moments are zero and the corresponding term need not be computed. It is seen that using a simple second order approximation with a modest extra computation yields a clear improvement in the accuracy of solution in this case. As another example, consider an iterative numerical experiment. An i th iteration proceeds as follows. A random matrix A of size and a random vector of size 1 is generated. Two random variables are considered: a 1 is normally distributed with Ea 1 = and Ea 1 =. A 1, and a is uniformly distributed with Ea = and Ea =.1 A,. For 1 realisations of a 1 and a, 1 A = A + a 1 + a 1 is computed and the mean solution E X = EA b is computed. Care taken was to ensure that the matrix A does not become too illconditioned. The solutions EX and X = A are computed. The following quantities are taken as measure of errors: L 1 = EX E X, and L = X E X. This entire experiment is repeated 1 times. It was found that L 1 < L held for 94 out of 1 times. In the remaining instances the assumption A < 1 was violated more than once. It is thought that this presents a convincing case for using a higher order approximation when better average case accuracy is required. 7 Conclusion A method to build a sequence of random vectors which converges almost surely to the solution of a system of linear equations under random perturbation is suggested. Probabilistic description of N th vector in this sequence may serve as an approximate description of the actual solution. When information about higher order moments for random perturbations is available, it may be used to build an accurate approximation to the expected value of the solution as well as to the confidence intervals around this expected value. The extension of these results to least squares problems with non-square matrices and examples of suggested applications will be reported elsewhere. References [1] S. Chandrasekaran, G. Golub, M. Gu, and A. Sayed, Parameter estimation in the presence of bounded data uncertaities, SIAM J. Matrix Anal. and Appl., vol. 19, pp. 35 5, [] L. E. Ghaoui and H. Lebret, Robust solutions to least-squares problems with uncertain data, SIAM J. Matrix Anal. Appl., vol. 18, pp , [3] G. Watson, Data fitting problem with bounded uncertainties in the data, SIAM J. Matrix Anal. Appl., vol., pp , 1. [4] A. Sayed, A framework for state space estimation with uncertain models, IEEE Trans. Automat. Contr., vol. 46, pp , 1. [5] A. Nagar and R. Powell, LFT/SDP approach to the uncertainty analysis for state estimation of water distribution systems, IEE Proceedings, vol. 149, pp ,. [6] J.Feinberg, On the universality of the probability distribution of the product B X of random matrices. [7] G. Grimmett and D. Stirzaker, Probability and Random Processes. Oxford University Press, 1. [8] R. Horn and C. Johnson, Matrix Analysis. Cambridge University Press, [9] W. Hoeffding, Probability inequalities for sums of bounded random variables, Amer. Statistical Asso. Journal, no. 3, pp. 13 3, [1] F. Paganini, A set-based approach for white noise modeling, IEEE Trans. Automat. Contr., vol. 41, pp , [11] X. Bombois, M. Gevers, G. Scorletti, and B. Anderson, Robustness analysis tools for an uncertainty set obtained by prediction error identification, Automatica, vol. 37, pp , 1.
Semidefinite and Second Order Cone Programming Seminar Fall 2012 Project: Robust Optimization and its Application of Robust Portfolio Optimization
Semidefinite and Second Order Cone Programming Seminar Fall 2012 Project: Robust Optimization and its Application of Robust Portfolio Optimization Instructor: Farid Alizadeh Author: Ai Kagawa 12/12/2012
More informationLikelihood Bounds for Constrained Estimation with Uncertainty
Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference 5 Seville, Spain, December -5, 5 WeC4. Likelihood Bounds for Constrained Estimation with Uncertainty
More informationHandout 8: Dealing with Data Uncertainty
MFE 5100: Optimization 2015 16 First Term Handout 8: Dealing with Data Uncertainty Instructor: Anthony Man Cho So December 1, 2015 1 Introduction Conic linear programming CLP, and in particular, semidefinite
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 informationFurther Results on Model Structure Validation for Closed Loop System Identification
Advances in Wireless Communications and etworks 7; 3(5: 57-66 http://www.sciencepublishinggroup.com/j/awcn doi:.648/j.awcn.735. Further esults on Model Structure Validation for Closed Loop System Identification
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 informationThM06-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 informationRobust Tuning of Power System Stabilizers Using Coefficient Diagram Method
International Journal of Electrical Engineering. ISSN 0974-2158 Volume 7, Number 2 (2014), pp. 257-270 International Research Publication House http://www.irphouse.com Robust Tuning of Power System Stabilizers
More informationOn the Power of Robust Solutions in Two-Stage Stochastic and Adaptive Optimization Problems
MATHEMATICS OF OPERATIONS RESEARCH Vol. 35, No., May 010, pp. 84 305 issn 0364-765X eissn 156-5471 10 350 084 informs doi 10.187/moor.1090.0440 010 INFORMS On the Power of Robust Solutions in Two-Stage
More informationECEN 420 LINEAR CONTROL SYSTEMS. Lecture 6 Mathematical Representation of Physical Systems II 1/67
1/67 ECEN 420 LINEAR CONTROL SYSTEMS Lecture 6 Mathematical Representation of Physical Systems II State Variable Models for Dynamic Systems u 1 u 2 u ṙ. Internal Variables x 1, x 2 x n y 1 y 2. y m Figure
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 informationUncertainty and Robustness for SISO Systems
Uncertainty and Robustness for SISO Systems ELEC 571L Robust Multivariable Control prepared by: Greg Stewart Outline Nature of uncertainty (models and signals). Physical sources of model uncertainty. Mathematical
More informationSOLVING POWER AND OPTIMAL POWER FLOW PROBLEMS IN THE PRESENCE OF UNCERTAINTY BY AFFINE ARITHMETIC
SOLVING POWER AND OPTIMAL POWER FLOW PROBLEMS IN THE PRESENCE OF UNCERTAINTY BY AFFINE ARITHMETIC Alfredo Vaccaro RTSI 2015 - September 16-18, 2015, Torino, Italy RESEARCH MOTIVATIONS Power Flow (PF) and
More informationOn the Power of Robust Solutions in Two-Stage Stochastic and Adaptive Optimization Problems
MATHEMATICS OF OPERATIONS RESEARCH Vol. xx, No. x, Xxxxxxx 00x, pp. xxx xxx ISSN 0364-765X EISSN 156-5471 0x xx0x 0xxx informs DOI 10.187/moor.xxxx.xxxx c 00x INFORMS On the Power of Robust Solutions in
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 informationStatic Output Feedback Stabilisation with H Performance for a Class of Plants
Static Output Feedback Stabilisation with H Performance for a Class of Plants E. Prempain and I. Postlethwaite Control and Instrumentation Research, Department of Engineering, University of Leicester,
More informationChance Constrained Input Design
5th IEEE Conference on Decision and Control and European Control Conference (CDC-ECC) Orlando, FL, USA, December -5, Chance Constrained Input Design Cristian R. Rojas, Dimitrios Katselis, Håkan Hjalmarsson,
More informationH Estimation. Speaker : R.Lakshminarayanan Guide : Prof. K.Giridhar. H Estimation p.1/34
H Estimation Speaker : R.Lakshminarayanan Guide : Prof. K.Giridhar H Estimation p.1/34 H Motivation The Kalman and Wiener Filters minimize the mean squared error between the true value and estimated values
More informationDenis ARZELIER arzelier
COURSE ON LMI OPTIMIZATION WITH APPLICATIONS IN CONTROL PART II.2 LMIs IN SYSTEMS CONTROL STATE-SPACE METHODS PERFORMANCE ANALYSIS and SYNTHESIS Denis ARZELIER www.laas.fr/ arzelier arzelier@laas.fr 15
More informationLong-Run Covariability
Long-Run Covariability Ulrich K. Müller and Mark W. Watson Princeton University October 2016 Motivation Study the long-run covariability/relationship between economic variables great ratios, long-run Phillips
More informationEstimation and Optimization: Gaps and Bridges. MURI Meeting June 20, Laurent El Ghaoui. UC Berkeley EECS
MURI Meeting June 20, 2001 Estimation and Optimization: Gaps and Bridges Laurent El Ghaoui EECS UC Berkeley 1 goals currently, estimation (of model parameters) and optimization (of decision variables)
More informationRobust linear optimization under general norms
Operations Research Letters 3 (004) 50 56 Operations Research Letters www.elsevier.com/locate/dsw Robust linear optimization under general norms Dimitris Bertsimas a; ;, Dessislava Pachamanova b, Melvyn
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 informationA Gain-Based Lower Bound Algorithm for Real and Mixed µ Problems
A Gain-Based Lower Bound Algorithm for Real and Mixed µ Problems Pete Seiler, Gary Balas, and Andrew Packard Abstract In this paper we present a new lower bound algorithm for real and mixed µ problems.
More informationBackward Error Estimation
Backward Error Estimation S. Chandrasekaran E. Gomez Y. Karant K. E. Schubert Abstract Estimation of unknowns in the presence of noise and uncertainty is an active area of study, because no method handles
More informationLearning Gaussian Process Models from Uncertain Data
Learning Gaussian Process Models from Uncertain Data Patrick Dallaire, Camille Besse, and Brahim Chaib-draa DAMAS Laboratory, Computer Science & Software Engineering Department, Laval University, Canada
More informationComparison of four state observer design algorithms for MIMO system
Archives of Control Sciences Volume 23(LIX), 2013 No. 2, pages 131 144 Comparison of four state observer design algorithms for MIMO system VINODH KUMAR. E, JOVITHA JEROME and S. AYYAPPAN A state observer
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 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 informationCHAPTER 5 ROBUSTNESS ANALYSIS OF THE CONTROLLER
114 CHAPTER 5 ROBUSTNESS ANALYSIS OF THE CONTROLLER 5.1 INTRODUCTION Robust control is a branch of control theory that explicitly deals with uncertainty in its approach to controller design. It also refers
More informationPERFORMANCE ANALYSIS OF CLOSED LOOP SYSTEM WITH A TAILOR MADE PARAMETERIZATION. Jianhong Wang, Hong Jiang and Yonghong Zhu
International Journal of Innovative Computing, Information and Control ICIC International c 208 ISSN 349-498 Volume 4, Number, February 208 pp. 8 96 PERFORMANCE ANALYSIS OF CLOSED LOOP SYSTEM WITH A TAILOR
More informationReview of Probability Theory
Review of Probability Theory Arian Maleki and Tom Do Stanford University Probability theory is the study of uncertainty Through this class, we will be relying on concepts from probability theory for deriving
More informationConvex relaxations of chance constrained optimization problems
Convex relaxations of chance constrained optimization problems Shabbir Ahmed School of Industrial & Systems Engineering, Georgia Institute of Technology, 765 Ferst Drive, Atlanta, GA 30332. May 12, 2011
More informationPlant friendly input design for system identification in closed loop
Plant friendly input design for system identification in closed loop Sridharakumar Narasimhan, Xavier Bombois Dept. of Chemical Engineering, IIT Madras, Chennai, India. (e-mail: sridharkrn@iitm.ac.in.
More informationLecture 6. Chapter 8: Robust Stability and Performance Analysis for MIMO Systems. Eugenio Schuster.
Lecture 6 Chapter 8: Robust Stability and Performance Analysis for MIMO Systems Eugenio Schuster schuster@lehigh.edu Mechanical Engineering and Mechanics Lehigh University Lecture 6 p. 1/73 6.1 General
More informationROBUST STABILITY AND PERFORMANCE ANALYSIS* [8 # ]
ROBUST STABILITY AND PERFORMANCE ANALYSIS* [8 # ] General control configuration with uncertainty [8.1] For our robustness analysis we use a system representation in which the uncertain perturbations are
More informationStructured Uncertainty and Robust Performance
Structured Uncertainty and Robust Performance ELEC 571L Robust Multivariable Control prepared by: Greg Stewart Devron Profile Control Solutions Outline Structured uncertainty: motivating example. Structured
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 informationOPTIMAL EXPERIMENT DESIGN IN CLOSED LOOP. KTH, Signals, Sensors and Systems, S Stockholm, Sweden.
OPTIMAL EXPERIMENT DESIGN IN CLOSED LOOP Henrik Jansson Håkan Hjalmarsson KTH, Signals, Sensors and Systems, S-00 44 Stockholm, Sweden. henrik.jansson@s3.kth.se Abstract: In this contribution we extend
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 informationControl Systems Design
ELEC4410 Control Systems Design Lecture 3, Part 2: Introduction to Affine Parametrisation School of Electrical Engineering and Computer Science Lecture 3, Part 2: Affine Parametrisation p. 1/29 Outline
More informationGaussian processes. Chuong B. Do (updated by Honglak Lee) November 22, 2008
Gaussian processes Chuong B Do (updated by Honglak Lee) November 22, 2008 Many of the classical machine learning algorithms that we talked about during the first half of this course fit the following pattern:
More informationPARAMETER DEPENDENT H CONTROLLER DESIGN BY FINITE DIMENSIONAL LMI OPTIMIZATION: APPLICATION TO TRADE-OFF DEPENDENT CONTROL
PARAMETER DEPEDET H COTROLLER DESIG BY FIITE DIMESIOAL LMI OPTIMIZATIO: APPLICATIO TO TRADE-OFF DEPEDET COTROL M Dinh, G Scorletti V Fromion E Magarotto GREYC Equipe Automatique, ISMRA 6 boulevard du Maréchal
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 informationAdaptive Dual Control
Adaptive Dual Control Björn Wittenmark Department of Automatic Control, Lund Institute of Technology Box 118, S-221 00 Lund, Sweden email: bjorn@control.lth.se Keywords: Dual control, stochastic control,
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 informationHowever, reliability analysis is not limited to calculation of the probability of failure.
Probabilistic Analysis probabilistic analysis methods, including the first and second-order reliability methods, Monte Carlo simulation, Importance sampling, Latin Hypercube sampling, and stochastic expansions
More informationSYNTHESIS OF ROBUST DISCRETE-TIME SYSTEMS BASED ON COMPARISON WITH STOCHASTIC MODEL 1. P. V. Pakshin, S. G. Soloviev
SYNTHESIS OF ROBUST DISCRETE-TIME SYSTEMS BASED ON COMPARISON WITH STOCHASTIC MODEL 1 P. V. Pakshin, S. G. Soloviev Nizhny Novgorod State Technical University at Arzamas, 19, Kalinina ul., Arzamas, 607227,
More informationLMI Based Model Order Reduction Considering the Minimum Phase Characteristic of the System
LMI Based Model Order Reduction Considering the Minimum Phase Characteristic of the System Gholamreza Khademi, Haniyeh Mohammadi, and Maryam Dehghani School of Electrical and Computer Engineering Shiraz
More informationThe ϵ-capacity of a gain matrix and tolerable disturbances: Discrete-time perturbed linear systems
IOSR Journal of Mathematics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 11, Issue 3 Ver. IV (May - Jun. 2015), PP 52-62 www.iosrjournals.org The ϵ-capacity of a gain matrix and tolerable disturbances:
More informationTen years of progress in Identification for Control. Outline
Ten years of progress in Identification for Control Design and Optimization of Restricted Complexity Controllers Grenoble Workshop, 15-16 January, 2003 Michel Gevers CESAME - UCL, Louvain-la-Neuve, Belgium
More informationRobust Optimization for Risk Control in Enterprise-wide Optimization
Robust Optimization for Risk Control in Enterprise-wide Optimization Juan Pablo Vielma Department of Industrial Engineering University of Pittsburgh EWO Seminar, 011 Pittsburgh, PA Uncertainty in Optimization
More information2nd Symposium on System, Structure and Control, Oaxaca, 2004
263 2nd Symposium on System, Structure and Control, Oaxaca, 2004 A PROJECTIVE ALGORITHM FOR STATIC OUTPUT FEEDBACK STABILIZATION Kaiyang Yang, Robert Orsi and John B. Moore Department of Systems Engineering,
More informationMIMO analysis: loop-at-a-time
MIMO robustness MIMO analysis: loop-at-a-time y 1 y 2 P (s) + + K 2 (s) r 1 r 2 K 1 (s) Plant: P (s) = 1 s 2 + α 2 s α 2 α(s + 1) α(s + 1) s α 2. (take α = 10 in the following numerical analysis) Controller:
More informationAlgorithms for Uncertainty Quantification
Algorithms for Uncertainty Quantification Lecture 9: Sensitivity Analysis ST 2018 Tobias Neckel Scientific Computing in Computer Science TUM Repetition of Previous Lecture Sparse grids in Uncertainty Quantification
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 informationSimultaneous State and Fault Estimation for Descriptor Systems using an Augmented PD Observer
Preprints of the 19th World Congress The International Federation of Automatic Control Simultaneous State and Fault Estimation for Descriptor Systems using an Augmented PD Observer Fengming Shi*, Ron J.
More informationRobust Pareto Design of GMDH-type Neural Networks for Systems with Probabilistic Uncertainties
. Hybrid GMDH-type algorithms and neural networks Robust Pareto Design of GMDH-type eural etworks for Systems with Probabilistic Uncertainties. ariman-zadeh, F. Kalantary, A. Jamali, F. Ebrahimi Department
More informationModel Order Reduction of Interval Systems using Mihailov Criterion and Factor Division Method
Model Order Reduction of Interval Systems using Mihailov Criterion Factor Division Method D. Kranthi Kumar Research Scholar S. K. Nagar Professor J. P. Tiwari Professor ABSTRACT This paper presents a mixed
More informationCONTROL SYSTEMS, ROBOTICS AND AUTOMATION Vol. XI Stochastic Stability - H.J. Kushner
STOCHASTIC STABILITY H.J. Kushner Applied Mathematics, Brown University, Providence, RI, USA. Keywords: stability, stochastic stability, random perturbations, Markov systems, robustness, perturbed systems,
More informationOptimal Control of Linear Systems with Stochastic Parameters for Variance Suppression: The Finite Time Case
Optimal Control of Linear Systems with Stochastic Parameters for Variance Suppression: The inite Time Case Kenji ujimoto Soraki Ogawa Yuhei Ota Makishi Nakayama Nagoya University, Department of Mechanical
More informationStructural Reliability
Structural Reliability Thuong Van DANG May 28, 2018 1 / 41 2 / 41 Introduction to Structural Reliability Concept of Limit State and Reliability Review of Probability Theory First Order Second Moment Method
More informationTHE INVERSE FUNCTION THEOREM
THE INVERSE FUNCTION THEOREM W. PATRICK HOOPER The implicit function theorem is the following result: Theorem 1. Let f be a C 1 function from a neighborhood of a point a R n into R n. Suppose A = Df(a)
More informationState feedback gain scheduling for linear systems with time-varying parameters
State feedback gain scheduling for linear systems with time-varying parameters Vinícius F. Montagner and Pedro L. D. Peres Abstract This paper addresses the problem of parameter dependent state feedback
More informationOn Identification of Cascade Systems 1
On Identification of Cascade Systems 1 Bo Wahlberg Håkan Hjalmarsson Jonas Mårtensson Automatic Control and ACCESS, School of Electrical Engineering, KTH, SE-100 44 Stockholm, Sweden. (bo.wahlberg@ee.kth.se
More informationROBUST BLIND CALIBRATION VIA TOTAL LEAST SQUARES
ROBUST BLIND CALIBRATION VIA TOTAL LEAST SQUARES John Lipor Laura Balzano University of Michigan, Ann Arbor Department of Electrical and Computer Engineering {lipor,girasole}@umich.edu ABSTRACT This paper
More informationThe Multivariate Gaussian Distribution [DRAFT]
The Multivariate Gaussian Distribution DRAFT David S. Rosenberg Abstract This is a collection of a few key and standard results about multivariate Gaussian distributions. I have not included many proofs,
More informationSTRUCTURE MATTERS: Some Notes on High Gain Observer Design for Nonlinear Systems. Int. Conf. on Systems, Analysis and Automatic Control 2012
Faculty of Electrical and Computer Engineering Institute of Control Theory STRUCTURE MATTERS: Some Notes on High Gain Observer Design for Nonlinear Systems Klaus Röbenack Int. Conf. on Systems, Analysis
More informationFINITE HORIZON ROBUST MODEL PREDICTIVE CONTROL USING LINEAR MATRIX INEQUALITIES. Danlei Chu, Tongwen Chen, Horacio J. Marquez
FINITE HORIZON ROBUST MODEL PREDICTIVE CONTROL USING LINEAR MATRIX INEQUALITIES Danlei Chu Tongwen Chen Horacio J Marquez Department of Electrical and Computer Engineering University of Alberta Edmonton
More informationsc Control Systems Design Q.1, Sem.1, Ac. Yr. 2010/11
sc46 - Control Systems Design Q Sem Ac Yr / Mock Exam originally given November 5 9 Notes: Please be reminded that only an A4 paper with formulas may be used during the exam no other material is to be
More informationROBUST ESTIMATOR FOR MULTIPLE INLIER STRUCTURES
ROBUST ESTIMATOR FOR MULTIPLE INLIER STRUCTURES Xiang Yang (1) and Peter Meer (2) (1) Dept. of Mechanical and Aerospace Engineering (2) Dept. of Electrical and Computer Engineering Rutgers University,
More informationLinear Matrix Inequalities in Robust Control. Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University MTNS 2002
Linear Matrix Inequalities in Robust Control Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University MTNS 2002 Objective A brief introduction to LMI techniques for Robust Control Emphasis on
More informationVariable Learning Rate LMS Based Linear Adaptive Inverse Control *
ISSN 746-7659, England, UK Journal of Information and Computing Science Vol., No. 3, 6, pp. 39-48 Variable Learning Rate LMS Based Linear Adaptive Inverse Control * Shuying ie, Chengjin Zhang School of
More informationGraph and Controller Design for Disturbance Attenuation in Consensus Networks
203 3th International Conference on Control, Automation and Systems (ICCAS 203) Oct. 20-23, 203 in Kimdaejung Convention Center, Gwangju, Korea Graph and Controller Design for Disturbance Attenuation in
More informationHybrid System Identification via Sparse Polynomial Optimization
2010 American Control Conference Marriott Waterfront, Baltimore, MD, USA June 30-July 02, 2010 WeA046 Hybrid System Identification via Sparse Polynomial Optimization Chao Feng, Constantino M Lagoa and
More informationModel Uncertainty and Robust Stability for Multivariable Systems
Model Uncertainty and Robust Stability for Multivariable Systems ELEC 571L Robust Multivariable Control prepared by: Greg Stewart Devron Profile Control Solutions Outline Representing model uncertainty.
More informationChapter 2. Classical Control System Design. Dutch Institute of Systems and Control
Chapter 2 Classical Control System Design Overview Ch. 2. 2. Classical control system design Introduction Introduction Steady-state Steady-state errors errors Type Type k k systems systems Integral Integral
More informationScalable robust hypothesis tests using graphical models
Scalable robust hypothesis tests using graphical models Umamahesh Srinivas ipal Group Meeting October 22, 2010 Binary hypothesis testing problem Random vector x = (x 1,...,x n ) R n generated from either
More informationJ. Sadeghi E. Patelli M. de Angelis
J. Sadeghi E. Patelli Institute for Risk and, Department of Engineering, University of Liverpool, United Kingdom 8th International Workshop on Reliable Computing, Computing with Confidence University of
More informationTHE STABLE EMBEDDING PROBLEM
THE STABLE EMBEDDING PROBLEM R. Zavala Yoé C. Praagman H.L. Trentelman Department of Econometrics, University of Groningen, P.O. Box 800, 9700 AV Groningen, The Netherlands Research Institute for Mathematics
More informationInterval solutions for interval algebraic equations
Mathematics and Computers in Simulation 66 (2004) 207 217 Interval solutions for interval algebraic equations B.T. Polyak, S.A. Nazin Institute of Control Sciences, Russian Academy of Sciences, 65 Profsoyuznaya
More informationHere represents the impulse (or delta) function. is an diagonal matrix of intensities, and is an diagonal matrix of intensities.
19 KALMAN FILTER 19.1 Introduction In the previous section, we derived the linear quadratic regulator as an optimal solution for the fullstate feedback control problem. The inherent assumption was that
More informationA Spectral Approach to Linear Bayesian Updating
A Spectral Approach to Linear Bayesian Updating Oliver Pajonk 1,2, Bojana V. Rosic 1, Alexander Litvinenko 1, and Hermann G. Matthies 1 1 Institute of Scientific Computing, TU Braunschweig, Germany 2 SPT
More informationChapter 8 SENSITIVITY ANALYSIS OF INTERVAL SYSTEM BY ROOT LOCUS APPROACH
Chapter 8 SENSITIVITY ANALYSIS OF INTERVAL SYSTEM BY ROOT LOCUS APPROACH SENSITIVITY ANALYSIS OF INTERVAL SYSTEM BY ROOT LOCUS APPROACH 8.1 INTRODUCTION Practical systems are difficult to model, as this
More informationUTILIZING PRIOR KNOWLEDGE IN ROBUST OPTIMAL EXPERIMENT DESIGN. EE & CS, The University of Newcastle, Australia EE, Technion, Israel.
UTILIZING PRIOR KNOWLEDGE IN ROBUST OPTIMAL EXPERIMENT DESIGN Graham C. Goodwin James S. Welsh Arie Feuer Milan Depich EE & CS, The University of Newcastle, Australia 38. EE, Technion, Israel. Abstract:
More informationcomponent risk analysis
273: Urban Systems Modeling Lec. 3 component risk analysis instructor: Matteo Pozzi 273: Urban Systems Modeling Lec. 3 component reliability outline risk analysis for components uncertain demand and uncertain
More informationBasics of Uncertainty Analysis
Basics of Uncertainty Analysis Chapter Six Basics of Uncertainty Analysis 6.1 Introduction As shown in Fig. 6.1, analysis models are used to predict the performances or behaviors of a product under design.
More informationHandout 6: Some Applications of Conic Linear Programming
ENGG 550: Foundations of Optimization 08 9 First Term Handout 6: Some Applications of Conic Linear Programming Instructor: Anthony Man Cho So November, 08 Introduction Conic linear programming CLP, and
More informationTwo Results About The Matrix Exponential
Two Results About The Matrix Exponential Hongguo Xu Abstract Two results about the matrix exponential are given. One is to characterize the matrices A which satisfy e A e AH = e AH e A, another is about
More informationA Characterization of the Hurwitz Stability of Metzler Matrices
29 American Control Conference Hyatt Regency Riverfront, St Louis, MO, USA June -2, 29 WeC52 A Characterization of the Hurwitz Stability of Metzler Matrices Kumpati S Narendra and Robert Shorten 2 Abstract
More informationState estimation of uncertain multiple model with unknown inputs
State estimation of uncertain multiple model with unknown inputs Abdelkader Akhenak, Mohammed Chadli, Didier Maquin and José Ragot Centre de Recherche en Automatique de Nancy, CNRS UMR 79 Institut National
More informationLINEAR QUADRATIC OPTIMAL CONTROL BASED ON DYNAMIC COMPENSATION. Received October 2010; revised March 2011
International Journal of Innovative Computing, Information and Control ICIC International c 22 ISSN 349-498 Volume 8, Number 5(B), May 22 pp. 3743 3754 LINEAR QUADRATIC OPTIMAL CONTROL BASED ON DYNAMIC
More informationRobust Principal Component Analysis
ELE 538B: Mathematics of High-Dimensional Data Robust Principal Component Analysis Yuxin Chen Princeton University, Fall 2018 Disentangling sparse and low-rank matrices Suppose we are given a matrix M
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 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 informationRegularization Parameter Estimation for Least Squares: A Newton method using the χ 2 -distribution
Regularization Parameter Estimation for Least Squares: A Newton method using the χ 2 -distribution Rosemary Renaut, Jodi Mead Arizona State and Boise State September 2007 Renaut and Mead (ASU/Boise) Scalar
More informationEFFICIENT SHAPE OPTIMIZATION USING POLYNOMIAL CHAOS EXPANSION AND LOCAL SENSITIVITIES
9 th ASCE Specialty Conference on Probabilistic Mechanics and Structural Reliability EFFICIENT SHAPE OPTIMIZATION USING POLYNOMIAL CHAOS EXPANSION AND LOCAL SENSITIVITIES Nam H. Kim and Haoyu Wang University
More informationKey words. conjugate gradients, normwise backward error, incremental norm estimation.
Proceedings of ALGORITMY 2016 pp. 323 332 ON ERROR ESTIMATION IN THE CONJUGATE GRADIENT METHOD: NORMWISE BACKWARD ERROR PETR TICHÝ Abstract. Using an idea of Duff and Vömel [BIT, 42 (2002), pp. 300 322
More informationME 132, Fall 2017, UC Berkeley, A. Packard 334 # 6 # 7 # 13 # 15 # 14
ME 132, Fall 2017, UC Berkeley, A. Packard 334 30.3 Fall 2017 Final # 1 # 2 # 3 # 4 # 5 # 6 # 7 # 8 NAME 20 15 20 15 15 18 15 20 # 9 # 10 # 11 # 12 # 13 # 14 # 15 # 16 18 12 12 15 12 20 18 15 Facts: 1.
More informationGeneralized Information Reuse for Optimization Under Uncertainty with Non-Sample Average Estimators
Generalized Information Reuse for Optimization Under Uncertainty with Non-Sample Average Estimators Laurence W Cook, Jerome P Jarrett, Karen E Willcox June 14, 2018 Abstract In optimization under uncertainty
More information