Advanced Process Control Tutorial Problem Set 2 Development of Control Relevant Models through System Identification
|
|
- Gwendoline Burns
- 6 years ago
- Views:
Transcription
1 Advanced Process Control Tutorial Problem Set 2 Development of Control Relevant Models through System Identification 1. Consider the time series x(k) = β 1 + β 2 k + w(k) where β 1 and β 2 are known constants and w(k) is a white noise process with variance σ 2. (a) Show that the mean of the moving average process y(k) = 1 p x(k j) 2p + 1 is β 1 + β 2 k. Is x(k) a stationary process? j= p (b) Find a transformation that produces a stationary process starting from x(k). (Hint: Consider transformation using backward difference operator, i.e. z(k) = (1 q 1 )x(k)) 2. Show that autocovariance function r(s, t) = E (v(s) µ v (s)) (v(k) µ v (k))] = E v(s)v(k)] µ v (s)µ v (k) where E v(s)] = µ v (s). 3. For a moving average process of the form x(k) = (1/2)w(k 2) + w(k 1) + 2w(k) (1/2)w(k + 1) where w(k) are independent with zero means and variance σ 2 w, determine the autocovariance and autocorrelation functions as a function of lag τ = s k. 4. Estimate the autocorrelation of the finite sequence u = {1, 2, 3, 4, 5, 6}. Comment on the relationship between r u,u (τ) and r u,u ( τ). 5. If h = {1, 2, 3, 4} and u = {5, 6, 7, 8}, estimate the cross-correlation r hu. 6. Consider two series x(k) = w(k) y(k) = w(k) θw(k 1) + u(k) where w(k) and u(k) are independent zero mean white noise sequences with variances σ 2 and λ 2, respectively, and θ is a unspecified constant. 1
2 (a) Express the autocorrelation function ρ y (τ) of sequence {y(k)} for τ = ±1, ±2,..as a function of σ 2,λ 2, and θ. (b) Determine cross-correlation function ρ x,y (τ) relating {x(k)} and {y(k)}. (c) Show that {x(k)} and {y(k)} are jointly stationary. (Series with constant means and aucovariance and cross covariance functions depending only on τ are said to be jointly stationary). 7. Consider a moving average process v(k) = e(k) + c 1 e(k 1) + c 2 e(k 2) (1) where {e(k)} is a zero mean white noise process with variance λ 2. Show that stochastic process {v(k)} has zero mean and auto-correlation R v (0) = E v(k), v(k)] = (1 + c c 2 2)λ 2 (2) R v (1) = E v(k), v(k 1)] = (c 1 + c 1 c 2 )λ 2 (3) R v (2) = E v(k), v(k 2)] = c 2 λ 2 (4) R v (k) = 0 for k > 2 (5) Note that {v(k)} is a typical example of colored noise. 8. Consider ARX model of the form y(k) = ay(k 1) + bu(k 1) + e(k) (6) It is desired to estimate the model parameters (a, b) using measurement data set {y(k) : k = 0, 1,...N} collected from an experiment in which input sequence {u(k) : k = 0, 1,...N} was injected into a system. (a) Show that the least square estimate of parameters generated from input-output data is given by y(k 1) 2 y(k 1)u(k 1) y(k 1)u(k 1) u(k 1) 2 where all summations are from k = 1 to N. ] â b ] = y(k)y(k 1) y(k)u(k 1) (7) (b) When data length is large (i.e. N ), show that equation (7) is equivalent to ] ] E y(k 1) 2 ] E y(k 1)u(k 1)] â E y(k)y(k 1)] = E y(k 1)u(k 1)] E u(k 1) 2 ] b E y(k)u(k 1)] (8) 2
3 or ] R y (0) R yu (0) R yu (0) R u (0) â b ] = R y(1) R yu (1) where R y (τ) represents auto-correlation function and R yu (τ) represents crosscorrelation function. (c) Defining regressor vector ] T ϕ(k) = y(k 1) u(k 1) (10) θ = â b ] T (11) show that equation (7) can be written as (9) E ϕ(k)ϕ(k) T ] θ = E ϕ(k)y(k)] (12) Hint: Show that Ω T Ω = E ϕ(k)ϕ(k) ] T Ω T Y = E ϕ(k)y(k)] ϕ(1) T y(1) ϕ(2) T Ω =... Y = y(2)... ϕ(n) T y(n) 9. Generalize the results of the previous for a general ARX model of the form y(k) = a 1 y(k 1)... a n y(k 1) + b 1 u(k 1) b n u(k n) + e(k) (13) 10. Model conversions (a) Consider OE model of the form y(k) = 2q 1 u(k) + v(k) 1 0.6q 1 Using long division, convert the model into the following form y(k) = h 1 u(k 1) h n u(k n) + v(k) where n is selected such that h i < 0.01 are neglected. How many terms are required and what can you say about h n as n increases? The resulting model is called finite impulse response model (FIR) and h i are called as impulse response coeffi cients (why?). 3
4 (b) Consider OE model of the form y(k) = 2q 1 u(k) + v(k) 1 1.5q 1 Can you find FIR model for this system? Justify your answer. (c) Consider AR model of the form v(k) = 1 e(k) 1 0.5q 1 where {e(k)} is a zero mean white noise signal with unit variance. Using long division, convert the model into moving average (MA) form y(k) = e(k) + h 1 e(k 1) h n e(k n) n is selected such that h i < 0.01 are neglected. (d) Consider AR model of the form v(k) = 1 (1 0.5q 1 )(1 0.25q 1 ) e(k) Using long division, convert the model into moving average (MA) form. (e) Consider AR model of the form v(k) = 1 (1 q 1 ) e(k) Using long division, is it possible to convert the model into moving average (MA) form? 11. Consider process governed by FIR equation of the form y(k) = h 1 u(k 1) + h 2 u(k 2) + e(k) (14) where {e(k)} is a sequence of independent normal N(0, λ) random variables. (a) Determine estimates of (ĥ1, ĥ2) when input signal {u(k)} is step input introduced at k = 0. (b) Make same investigation as part(a) when the input signal {u(k)} is a white noise with unit variance. 4
5 12. Consider data generated by the discrete time system System : y(k) = h 1 u(k 1) + h 2 u(k 2) + e(k) (15) where {e(k)} is a sequence of independent normal N(0, 1) random variables. Assume that parameter h of the model Model : y(k) = hu(k) (16) is determined by least square. (a) Determine estimates obtained for large observation sets when the input u(k) is a step function. (This is a simple illustration of the problem of fitting a low order model to a data generated by a complex system. The result obtained will critically depend on the character of the input signal.) (b) Make same investigation as part (a) when the input signal is white noise with unit variance. 13. Consider FIR model of the form y(k) = h 1 u(k 1) h N u(k N) + v(k) (17) show that least square estimates of impulse response coeffi cients are given by equation (12) where ] T ϕ(k) = u(k 1)... u(k N) (18) ] T θ = ĥ1... ĥ N (19) In other words, generalize results of Problem 8 to a general FIR model 14. If it is desired to identify parameters of FIR model (17), taking clues from the previous problem, what is the requirement on rank of matrix E ϕ(k)ϕ(k) ] T? This condition is called as persistency of excitation. 15. For a FIR model, show that parameter estimates are unbiased if {v(k)} is a zero mean sequence. 16. Consider discrete time system given by equation (6) where the input signal {u(k)} and noise {e(k)} are sequences on independent random variables with zero mean and standard deviations σ and λ, respectively.. Determine the covariance of parameter estimates obtained for large observation sets. 5
6 17. Consider discrete time system given by equation y(k) = a 0 y(k 1) + b 0 u(k 1) + e(k) + c 0 e(k 1) (20) where the input signal {u(k)} and noise {e(k)} are sequences on independent random variables with zero mean and standard deviations σ and λ, respectively. Assume that a model of the form y(k) = ay(k 1) + bu(k 1) + ε(k) (21) are estimated by least squares. Determine the asymptotic values of the estimates when (a) {u(k)} is a zero mean white noise process with standard deviation σ (b) {u(k)} is step input of magnitude σ (c) In particular, compare the estimated values (â, b) with the true values (a 0, b 0 ) for the following system a 0 = 0.8 ; b 0 = 1 ; c 0 = 0.5 (22) for the cases (a) σ = 1, λ = 0.1 (b) σ = 1, λ =. By comparing the estimates for cases (a) and (b) with true values, what can you conclude about the effect of signal to noise ration (σ 2 /λ 2 ) on the parameter estimates? 18. Consider a discrete time model v(k) = a + b k + e(k) (23) where {e(k)} is a sequence of independent normal N(0, λ) random variables. Determine least square estimates of model parameters and covariance of the estimates. Discuss behavior of the estimates as the number of data points increases. 19. Consider data generated by y(k) = b + e(k) ; k = 1, 2,...N (24) where {e(k) : k = 1, 3, 4,...} is a sequence of independent random variables. Furthermore, assume that there is a large error at k = 2, i.e., e(2) = A where A is a large number. Determine the estimate obtained and discuss how it depends on A.(This is a simple example that shows how sensitive the least square estimate is with respect to occasional large errors.) 6
7 20. Suppose that we wish to identify a plant that is operating in closed loop as follows Plant dynamics : y(k) = ay(k 1) + bu(k 1) + ε(k) (25) Feedback control law : u(k) = βy(k) (26) where {e(k)} is a sequence of independent normal N(0, λ) random variables. (a) Show that we cannot identify parameters (a, b) from observations of y and u, even when β is known. (b) Assume that an external independent perturbation was introduced in input signal u(k) = βy(k) + r(k) (27) where {r(k)} is a sequence of independent normal N(0, σ) random variables. Show that it is now possible to recover estimates of open loop model parameters using the closed loop data. (Note: Here {r(k)} has been taken as a zero mean white noise sequence to simplify the analysis. In practice, an independent PRBS signal is added to manipulated input to make the model parameters identifiable in closed loop conditions.) 21. The English mathematician Richardson has proposed the following simple model for the arms race between two countries x(k + 1) = ax(k) + by(k) + f (28) y(k + 1) = cx(k) + dy(k) + g (29) where x(k) and y(k) are yearly expenditures on arms of the two nations and (a, b, c, d, f, g) are model parameters. The following data has been obtained from World Armaments and Disarmaments Year Book 1982Determine the parameters of the model by least squares and investigate the stability of the model. 22. Consider an ARMA model of the form which is equivalent to y(k) = ay(k 1) + e(k) + ce(k 1) (30) y(k) = H(q)e(k) = 1 + cq 1 e(k) (31) 1 + aq 1 {e(k)} is a sequence of independent normal N(0, λ) random variables. Develop 1 step ahead predictor for ŷ(k + 1 k), which uses only the current and the past measurements of y. 7
8 23. Consider an ARMAX model of the form which is equivalent to y(k) = ay(k 1) + bu(k 1) + e(k) + ce(k 1) (32) y(k) = G(q)u(k) + H(q)e(k) = bq cq 1 u(k) + e(k) (33) 1 + aq aq 1 {e(k)} is a sequence of independent normal N(0, λ) random variables. Develop 1-step ahead predictor for ŷ(k + 1 k), which uses only the current and the past measurements of y. 24. Consider Box-Jenkin s model y(k) = G(q)u(k) + H(q)e(k) Derive one step prediction G(q) = q + b q + a H(q) = q + c q + d ŷ(k k 1) = H(q)] 1 G(q)u(k) + 1 (H(q)) 1] y(k) y(k) = ŷ(k k 1) + e(k) and express dynamics of ŷ(k k 1) as a time domain difference equation. 25. Consider moving average (MA) process y(k) = H(q)e(k) (34) H(q) = 1 1.1q q 2 (35) 8
9 Compute H 1 (q) as an infinite expansion by long division and develop an auto-regressive model of the form e(k) = H 1 (q)y(k) (36) This model facilitates estimation of noise e(k) based on current and past measurements of y(k). 26. Given an ARMAX model of the form y(k) = B(q) C(q) u(k) + A(q) A(q) e(k) = Rearrange this model as 0.1q q 1 0.2q 1 u(k) + e(k) (37) q 1 y(k) = C 1 (q)b(q) C 1 (q)a(q) u(k) + 1 e(k) (38) C 1 (q)a(q) Compute C 1 (q) as an infinite expansion by long division and truncate the expansion after finite number of terms when coeffi cients become small, i.e. C 1 T (q) 1 + c 1q c n q n (39) Using this truncated C 1 T (q), express the model in ARX form y(k) = B(q) Ã(q) u(k) + 1 e(k) (40) Ã(q) Ã(q) = C 1 T (q)a(q) ; B(q) = C 1 T (q)b(q) (41) This simple calculation will illustrate how a low order ARMAX model can be approximated by a high order ARX model. 27. Consider transfer functions G 1 (q) = G 2 (q) = H(q) = (q 0.5)(q + 0.5) (q 1)(q 2 1.5q + 0.7) (q 0.2)(q + 0.2) (q 1)(q 2 1.5q + 0.7) (q 0.8) (q 2 1.5q + 0.7) Derive state-space realization using observable canonical form for the following systems (cases (a) to (d)) (a) y(k) = G 1 (q)u 1 (k) + v(k) 9
10 (b) y(k) = G 1 (q)u 1 (k) + G 2 (q)u 2 (k) + v(k) (c) y(k) = G 1 (q)u 1 (k) + H(q)e(k) (d) y(k) = G 1 (q)u(k) + G 2 (q)u(k) + H(q)e(k) (e) Given that sequence {e(k)} is a zero mean white noise sequence with standard deviation equal to 0.5, express the resulting state space models for cases (c) and (d) in the form x(k + 1) = Φx(k) + Γu(k) + w(k) y(k) = Cx(k) + e(k) and estimate covariance of white noise sequence {w(k)}. (f) Derive state-space realization using controllable canonical form for case (a). 28. Derive a state realization for ] y 1 (k) 1 = y 2 (k) (q 2 1.5q + 0.8) q q 1.5 q 0.5 q ] u 1 (k) u 2 (k) ] + v 1 (k) v 2 (k) ] in controllable and observable canonical forms. 29. A system is represented by G(s) = 3 (s + 4)(s + 1) (a) Derive continuous time state-space realizations dx = Ax + Bu ; y = Cx dt in (a) Controllable canonical form (b) Observable canonical form (b) Convert each of the continuous time state space models into discrete state space form x(k + 1) = Φx(k) + Bu(k) ; y(k) = Cx(k) Is canonical structure in continuous time preserved after discretization? Show that both the discrete realizations have identical transfer function G(q). (c) If canonical structures are not preserved after discretization, derive discrete state realizations in (a) Controllable canonical form (b) Observable canonical form starting from G(q). 30. Derive a state realizations for ] y 1 (t) 1 = y 2 (t) (s 2 + 3s + 2) s s 2 s 3 s + 2 ] u 1 (t) u 2 (t) ] in controllable and observable canonical forms. 10
EL1820 Modeling of Dynamical Systems
EL1820 Modeling of Dynamical Systems Lecture 9 - Parameter estimation in linear models Model structures Parameter estimation via prediction error minimization Properties of the estimate: bias and variance
More informationEL1820 Modeling of Dynamical Systems
EL1820 Modeling of Dynamical Systems Lecture 10 - System identification as a model building tool Experiment design Examination and prefiltering of data Model structure selection Model validation Lecture
More informationEECE Adaptive Control
EECE 574 - Adaptive Control Basics of System Identification Guy Dumont Department of Electrical and Computer Engineering University of British Columbia January 2010 Guy Dumont (UBC) EECE574 - Basics of
More informationNonlinear System Identification Using MLP Dr.-Ing. Sudchai Boonto
Dr-Ing Sudchai Boonto Department of Control System and Instrumentation Engineering King Mongkut s Unniversity of Technology Thonburi Thailand Nonlinear System Identification Given a data set Z N = {y(k),
More informationIdentification of ARX, OE, FIR models with the least squares method
Identification of ARX, OE, FIR models with the least squares method CHEM-E7145 Advanced Process Control Methods Lecture 2 Contents Identification of ARX model with the least squares minimizing the equation
More informationAnalysis of Discrete-Time Systems
TU Berlin Discrete-Time Control Systems 1 Analysis of Discrete-Time Systems Overview Stability Sensitivity and Robustness Controllability, Reachability, Observability, and Detectabiliy TU Berlin Discrete-Time
More informationCONTROL SYSTEMS, ROBOTICS, AND AUTOMATION - Vol. V - Prediction Error Methods - Torsten Söderström
PREDICTIO ERROR METHODS Torsten Söderström Department of Systems and Control, Information Technology, Uppsala University, Uppsala, Sweden Keywords: prediction error method, optimal prediction, identifiability,
More informationAnalysis of Discrete-Time Systems
TU Berlin Discrete-Time Control Systems TU Berlin Discrete-Time Control Systems 2 Stability Definitions We define stability first with respect to changes in the initial conditions Analysis of Discrete-Time
More informationControl Design. Lecture 9: State Feedback and Observers. Two Classes of Control Problems. State Feedback: Problem Formulation
Lecture 9: State Feedback and s [IFAC PB Ch 9] State Feedback s Disturbance Estimation & Integral Action Control Design Many factors to consider, for example: Attenuation of load disturbances Reduction
More informationData mining for system identi cation
LERTEKNIK REG AU T O MA RO T IC C O N T L applications to process André Carvalho Bittencourt Automatic Control, Linköping Sweden Outline 1. Problem formulation. Theoretical guiding principles Modeling
More informationModelling and Control of Dynamic Systems. Stability of Linear Systems. Sven Laur University of Tartu
Modelling and Control of Dynamic Systems Stability of Linear Systems Sven Laur University of Tartu Motivating Example Naive open-loop control r[k] Controller Ĉ[z] u[k] ε 1 [k] System Ĝ[z] y[k] ε 2 [k]
More informationOn Input Design for System Identification
On Input Design for System Identification Input Design Using Markov Chains CHIARA BRIGHENTI Masters Degree Project Stockholm, Sweden March 2009 XR-EE-RT 2009:002 Abstract When system identification methods
More informationWe are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors
We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists 3,800 116,000 120M Open access books available International authors and editors Downloads Our
More informationEE538 Final Exam Fall 2007 Mon, Dec 10, 8-10 am RHPH 127 Dec. 10, Cover Sheet
EE538 Final Exam Fall 2007 Mon, Dec 10, 8-10 am RHPH 127 Dec. 10, 2007 Cover Sheet Test Duration: 120 minutes. Open Book but Closed Notes. Calculators allowed!! This test contains five problems. Each of
More informationAbstract. 1 Introduction
Relationship between state-spaceand input-output models via observer Markov parameters M.Q. Phan,' R. W. Longman* ^Department of Mechanical and Aerospace Engineering, Princeton University, Princeton, NJ
More informationLECTURES 2-3 : Stochastic Processes, Autocorrelation function. Stationarity.
LECTURES 2-3 : Stochastic Processes, Autocorrelation function. Stationarity. Important points of Lecture 1: A time series {X t } is a series of observations taken sequentially over time: x t is an observation
More informationIdentification of Linear Systems
Identification of Linear Systems Johan Schoukens http://homepages.vub.ac.be/~jschouk Vrije Universiteit Brussel Department INDI /67 Basic goal Built a parametric model for a linear dynamic system from
More informationModule 9: State Feedback Control Design Lecture Note 1
Module 9: State Feedback Control Design Lecture Note 1 The design techniques described in the preceding lectures are based on the transfer function of a system. In this lecture we would discuss the state
More informationTime Series Analysis
Time Series Analysis hm@imm.dtu.dk Informatics and Mathematical Modelling Technical University of Denmark DK-2800 Kgs. Lyngby 1 Outline of the lecture Input-Output systems The z-transform important issues
More informationStochastic Process II Dr.-Ing. Sudchai Boonto
Dr-Ing Sudchai Boonto Department of Control System and Instrumentation Engineering King Mongkuts Unniversity of Technology Thonburi Thailand Random process Consider a random experiment specified by the
More informationOn the Equivalence of OKID and Time Series Identification for Markov-Parameter Estimation
On the Equivalence of OKID and Time Series Identification for Markov-Parameter Estimation P V Albuquerque, M Holzel, and D S Bernstein April 5, 2009 Abstract We show the equivalence of Observer/Kalman
More informationIV. Covariance Analysis
IV. Covariance Analysis Autocovariance Remember that when a stochastic process has time values that are interdependent, then we can characterize that interdependency by computing the autocovariance function.
More informationModule 3. Descriptive Time Series Statistics and Introduction to Time Series Models
Module 3 Descriptive Time Series Statistics and Introduction to Time Series Models Class notes for Statistics 451: Applied Time Series Iowa State University Copyright 2015 W Q Meeker November 11, 2015
More information6.435, System Identification
SET 6 System Identification 6.435 Parametrized model structures One-step predictor Identifiability Munther A. Dahleh 1 Models of LTI Systems A complete model u = input y = output e = noise (with PDF).
More information1 Linear Difference Equations
ARMA Handout Jialin Yu 1 Linear Difference Equations First order systems Let {ε t } t=1 denote an input sequence and {y t} t=1 sequence generated by denote an output y t = φy t 1 + ε t t = 1, 2,... with
More informationCDS Final Exam
CDS 22 - Final Exam Instructor: Danielle C. Tarraf December 4, 2007 INSTRUCTIONS : Please read carefully! () Description & duration of the exam: The exam consists of 6 problems. You have a total of 24
More informationAutomatic Control 2. Model reduction. Prof. Alberto Bemporad. University of Trento. Academic year
Lecture: Automatic Control 2 Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011 1 / 17 Lecture:
More informationSystem Identification
System Identification Arun K. Tangirala Department of Chemical Engineering IIT Madras July 26, 2013 Module 6 Lecture 1 Arun K. Tangirala System Identification July 26, 2013 1 Objectives of this Module
More informationEEE582 Homework Problems
EEE582 Homework Problems HW. Write a state-space realization of the linearized model for the cruise control system around speeds v = 4 (Section.3, http://tsakalis.faculty.asu.edu/notes/models.pdf). Use
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 informationRural/Urban Migration: The Dynamics of Eigenvectors
* Analysis of the Dynamic Structure of a System * Rural/Urban Migration: The Dynamics of Eigenvectors EGR 326 April 11, 2019 1. Develop the system model and create the Matlab/Simulink model 2. Plot and
More informationRoss Bettinger, Analytical Consultant, Seattle, WA
ABSTRACT DYNAMIC REGRESSION IN ARIMA MODELING Ross Bettinger, Analytical Consultant, Seattle, WA Box-Jenkins time series models that contain exogenous predictor variables are called dynamic regression
More informationEcon 623 Econometrics II Topic 2: Stationary Time Series
1 Introduction Econ 623 Econometrics II Topic 2: Stationary Time Series In the regression model we can model the error term as an autoregression AR(1) process. That is, we can use the past value of the
More informationLecture Note #6 (Chap.10)
System Modeling and Identification Lecture Note #6 (Chap.) CBE 7 Korea University Prof. Dae Ryook Yang Chap. Model Approximation Model approximation Simplification, approximation and order reduction of
More informationLecture 1: Introduction to System Modeling and Control. Introduction Basic Definitions Different Model Types System Identification
Lecture 1: Introduction to System Modeling and Control Introduction Basic Definitions Different Model Types System Identification What is Mathematical Model? A set of mathematical equations (e.g., differential
More informationIdentification in closed-loop, MISO identification, practical issues of identification
Identification in closed-loop, MISO identification, practical issues of identification CHEM-E7145 Advanced Process Control Methods Lecture 4 Contents Identification in practice Identification in closed-loop
More informationLinear Model Under General Variance Structure: Autocorrelation
Linear Model Under General Variance Structure: Autocorrelation A Definition of Autocorrelation In this section, we consider another special case of the model Y = X β + e, or y t = x t β + e t, t = 1,..,.
More information18.S096 Problem Set 4 Fall 2013 Time Series Due Date: 10/15/2013
18.S096 Problem Set 4 Fall 2013 Time Series Due Date: 10/15/2013 1. Covariance Stationary AR(2) Processes Suppose the discrete-time stochastic process {X t } follows a secondorder auto-regressive process
More informationControl Systems Lab - SC4070 System Identification and Linearization
Control Systems Lab - SC4070 System Identification and Linearization Dr. Manuel Mazo Jr. Delft Center for Systems and Control (TU Delft) m.mazo@tudelft.nl Tel.:015-2788131 TU Delft, February 13, 2015 (slides
More informationExam in Automatic Control II Reglerteknik II 5hp (1RT495)
Exam in Automatic Control II Reglerteknik II 5hp (1RT495) Date: August 4, 018 Venue: Bergsbrunnagatan 15 sal Responsible teacher: Hans Rosth. Aiding material: Calculator, mathematical handbooks, textbooks
More informationy k = ( ) x k + v k. w q wk i 0 0 wk
Four telling examples of Kalman Filters Example : Signal plus noise Measurement of a bandpass signal, center frequency.2 rad/sec buried in highpass noise. Dig out the quadrature part of the signal while
More informationApplication of Modified Multi Model Predictive Control Algorithm to Fluid Catalytic Cracking Unit
Application of Modified Multi Model Predictive Control Algorithm to Fluid Catalytic Cracking Unit Nafay H. Rehman 1, Neelam Verma 2 Student 1, Asst. Professor 2 Department of Electrical and Electronics
More informationWeek 5 Quantitative Analysis of Financial Markets Characterizing Cycles
Week 5 Quantitative Analysis of Financial Markets Characterizing Cycles Christopher Ting http://www.mysmu.edu/faculty/christophert/ Christopher Ting : christopherting@smu.edu.sg : 6828 0364 : LKCSB 5036
More informationIndustrial Model Predictive Control
Industrial Model Predictive Control Emil Schultz Christensen Kongens Lyngby 2013 DTU Compute-M.Sc.-2013-49 Technical University of Denmark DTU Compute Matematiktovet, Building 303B, DK-2800 Kongens Lyngby,
More informationLecture 7: Discrete-time Models. Modeling of Physical Systems. Preprocessing Experimental Data.
ISS0031 Modeling and Identification Lecture 7: Discrete-time Models. Modeling of Physical Systems. Preprocessing Experimental Data. Aleksei Tepljakov, Ph.D. October 21, 2015 Discrete-time Transfer Functions
More informationIDENTIFICATION OF A TWO-INPUT SYSTEM: VARIANCE ANALYSIS
IDENTIFICATION OF A TWO-INPUT SYSTEM: VARIANCE ANALYSIS M Gevers,1 L Mišković,2 D Bonvin A Karimi Center for Systems Engineering and Applied Mechanics (CESAME) Université Catholique de Louvain B-1348 Louvain-la-Neuve,
More informationAn Exponentially Weighted Moving Average Method for Identification and Monitoring of Stochastic Systems
Ind. Eng. Chem. Res. 2008, 47, 8239 8249 8239 PROCESS DESIGN AND CONTROL An Exponentially Weighted Moving Average Method for Identification and Monitoring of Stochastic Systems Shyh-Hong Hwang,* Ho-Tsen
More informationMoving Average (MA) representations
Moving Average (MA) representations The moving average representation of order M has the following form v[k] = MX c n e[k n]+e[k] (16) n=1 whose transfer function operator form is MX v[k] =H(q 1 )e[k],
More information13. Power Spectrum. For a deterministic signal x(t), the spectrum is well defined: If represents its Fourier transform, i.e., if.
For a deterministic signal x(t), the spectrum is well defined: If represents its Fourier transform, i.e., if jt X ( ) = xte ( ) dt, (3-) then X ( ) represents its energy spectrum. his follows from Parseval
More informationUniversity of Oxford. Statistical Methods Autocorrelation. Identification and Estimation
University of Oxford Statistical Methods Autocorrelation Identification and Estimation Dr. Órlaith Burke Michaelmas Term, 2011 Department of Statistics, 1 South Parks Road, Oxford OX1 3TG Contents 1 Model
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 informationLecture 5: Recurrent Neural Networks
1/25 Lecture 5: Recurrent Neural Networks Nima Mohajerin University of Waterloo WAVE Lab nima.mohajerin@uwaterloo.ca July 4, 2017 2/25 Overview 1 Recap 2 RNN Architectures for Learning Long Term Dependencies
More informationAn All-Interaction Matrix Approach to Linear and Bilinear System Identification
An All-Interaction Matrix Approach to Linear and Bilinear System Identification Minh Q. Phan, Francesco Vicario, Richard W. Longman, and Raimondo Betti Abstract This paper is a brief introduction to the
More informationTime Series Analysis. James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY
Time Series Analysis James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY & Contents PREFACE xiii 1 1.1. 1.2. Difference Equations First-Order Difference Equations 1 /?th-order Difference
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 informationOptimal control and estimation
Automatic Control 2 Optimal control and estimation Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011
More information6.435, System Identification
System Identification 6.435 SET 3 Nonparametric Identification Munther A. Dahleh 1 Nonparametric Methods for System ID Time domain methods Impulse response Step response Correlation analysis / time Frequency
More informationProblem Set 2: Box-Jenkins methodology
Problem Set : Box-Jenkins methodology 1) For an AR1) process we have: γ0) = σ ε 1 φ σ ε γ0) = 1 φ Hence, For a MA1) process, p lim R = φ γ0) = 1 + θ )σ ε σ ε 1 = γ0) 1 + θ Therefore, p lim R = 1 1 1 +
More informationAN IDENTIFICATION ALGORITHM FOR ARMAX SYSTEMS
AN IDENTIFICATION ALGORITHM FOR ARMAX SYSTEMS First the X, then the AR, finally the MA Jan C. Willems, K.U. Leuven Workshop on Observation and Estimation Ben Gurion University, July 3, 2004 p./2 Joint
More informationDynamic and Time Series Modeling for Process Control
Dynamic and ime Series Modeling for Process Control Sachin C. Patwardhan Dept. of Chemical Engineering Why Mathematical Modeling? Key Component of All Advanced Monitoring, Control and Optimization Schemes
More informationClosed-Loop Identification of Unstable Systems Using Noncausal FIR Models
23 American Control Conference (ACC) Washington, DC, USA, June 7-9, 23 Closed-Loop Identification of Unstable Systems Using Noncausal FIR Models Khaled Aljanaideh, Benjamin J. Coffer, and Dennis S. Bernstein
More information4F3 - Predictive Control
4F3 Predictive Control - Lecture 2 p 1/23 4F3 - Predictive Control Lecture 2 - Unconstrained Predictive Control Jan Maciejowski jmm@engcamacuk 4F3 Predictive Control - Lecture 2 p 2/23 References Predictive
More information4F3 - Predictive Control
4F3 Predictive Control - Lecture 3 p 1/21 4F3 - Predictive Control Lecture 3 - Predictive Control with Constraints Jan Maciejowski jmm@engcamacuk 4F3 Predictive Control - Lecture 3 p 2/21 Constraints on
More informationTIME SERIES ANALYSIS. Forecasting and Control. Wiley. Fifth Edition GWILYM M. JENKINS GEORGE E. P. BOX GREGORY C. REINSEL GRETA M.
TIME SERIES ANALYSIS Forecasting and Control Fifth Edition GEORGE E. P. BOX GWILYM M. JENKINS GREGORY C. REINSEL GRETA M. LJUNG Wiley CONTENTS PREFACE TO THE FIFTH EDITION PREFACE TO THE FOURTH EDITION
More informationInfinite Horizon LQ. Given continuous-time state equation. Find the control function u(t) to minimize
Infinite Horizon LQ Given continuous-time state equation x = Ax + Bu Find the control function ut) to minimize J = 1 " # [ x T t)qxt) + u T t)rut)] dt 2 0 Q $ 0, R > 0 and symmetric Solution is obtained
More informationDiscrete time processes
Discrete time processes Predictions are difficult. Especially about the future Mark Twain. Florian Herzog 2013 Modeling observed data When we model observed (realized) data, we encounter usually the following
More informationLinear Systems. Manfred Morari Melanie Zeilinger. Institut für Automatik, ETH Zürich Institute for Dynamic Systems and Control, ETH Zürich
Linear Systems Manfred Morari Melanie Zeilinger Institut für Automatik, ETH Zürich Institute for Dynamic Systems and Control, ETH Zürich Spring Semester 2016 Linear Systems M. Morari, M. Zeilinger - Spring
More information1 Introduction to Generalized Least Squares
ECONOMICS 7344, Spring 2017 Bent E. Sørensen April 12, 2017 1 Introduction to Generalized Least Squares Consider the model Y = Xβ + ɛ, where the N K matrix of regressors X is fixed, independent of the
More informationTime Series Analysis
Time Series Analysis Christopher Ting http://mysmu.edu.sg/faculty/christophert/ christopherting@smu.edu.sg Quantitative Finance Singapore Management University March 3, 2017 Christopher Ting Week 9 March
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 informationAn Algorithm for Finding Process Identification Intervals from Normal Operating Data
Processes 015, 3, 357-383; doi:10.3390/pr300357 OPEN ACCESS processes ISSN 7-9717 www.mdpi.com/journal/processes Article An Algorithm for Finding Process Identification Intervals from Normal Operating
More informationEECE Adaptive Control
EECE 574 - Adaptive Control Recursive Identification in Closed-Loop and Adaptive Control Guy Dumont Department of Electrical and Computer Engineering University of British Columbia January 2010 Guy Dumont
More informationOnline monitoring of MPC disturbance models using closed-loop data
Online monitoring of MPC disturbance models using closed-loop data Brian J. Odelson and James B. Rawlings Department of Chemical Engineering University of Wisconsin-Madison Online Optimization Based Identification
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 informationText of Intelligent Control. M. Yamakita
Text of Intelligent Control M. Yamakita 1 0.1 Sampled data system and its mathematical representation 0.1.1 Discretization of continuous linear time invariant (LTI)systems We consider a system represented
More informationEconometría 2: Análisis de series de Tiempo
Econometría 2: Análisis de series de Tiempo Karoll GOMEZ kgomezp@unal.edu.co http://karollgomez.wordpress.com Segundo semestre 2016 II. Basic definitions A time series is a set of observations X t, each
More informationLecture 4a: ARMA Model
Lecture 4a: ARMA Model 1 2 Big Picture Most often our goal is to find a statistical model to describe real time series (estimation), and then predict the future (forecasting) One particularly popular model
More informationProcess Dynamics & Control LECTURE 1: INTRODUCTION OF MODEL PREDICTIVE CONTROL A Multivariable Control Technique for the Process Industry
Process Dynamics & Control LECTURE 1: INTRODUCTION OF MODEL PREDICTIVE CONTROL A Multivariable Control Technique for the Process Industry Jong Min Lee Chemical and Biological Engineering Seoul National
More informationErrors-in-variables identification through covariance matching: Analysis of a colored measurement noise case
008 American Control Conference Westin Seattle Hotel Seattle Washington USA June -3 008 WeB8.4 Errors-in-variables identification through covariance matching: Analysis of a colored measurement noise case
More informationProbability Space. J. McNames Portland State University ECE 538/638 Stochastic Signals Ver
Stochastic Signals Overview Definitions Second order statistics Stationarity and ergodicity Random signal variability Power spectral density Linear systems with stationary inputs Random signal memory Correlation
More informationQuis custodiet ipsos custodes?
Quis custodiet ipsos custodes? James B. Rawlings, Megan Zagrobelny, Luo Ji Dept. of Chemical and Biological Engineering, Univ. of Wisconsin-Madison, WI, USA IFAC Conference on Nonlinear Model Predictive
More informationA summary of Modeling and Simulation
A summary of Modeling and Simulation Text-book: Modeling of dynamic systems Lennart Ljung and Torkel Glad Content What re Models for systems and signals? Basic concepts Types of models How to build a model
More informationENSC327 Communications Systems 19: Random Processes. Jie Liang School of Engineering Science Simon Fraser University
ENSC327 Communications Systems 19: Random Processes Jie Liang School of Engineering Science Simon Fraser University 1 Outline Random processes Stationary random processes Autocorrelation of random processes
More informationwhere r n = dn+1 x(t)
Random Variables Overview Probability Random variables Transforms of pdfs Moments and cumulants Useful distributions Random vectors Linear transformations of random vectors The multivariate normal distribution
More informationEmpirical Market Microstructure Analysis (EMMA)
Empirical Market Microstructure Analysis (EMMA) Lecture 3: Statistical Building Blocks and Econometric Basics Prof. Dr. Michael Stein michael.stein@vwl.uni-freiburg.de Albert-Ludwigs-University of Freiburg
More informationStatistical signal processing
Statistical signal processing Short overview of the fundamentals Outline Random variables Random processes Stationarity Ergodicity Spectral analysis Random variable and processes Intuition: A random variable
More informationStochastic Processes. A stochastic process is a function of two variables:
Stochastic Processes Stochastic: from Greek stochastikos, proceeding by guesswork, literally, skillful in aiming. A stochastic process is simply a collection of random variables labelled by some parameter:
More informationXIII Simpósio Brasileiro de Automação Inteligente Porto Alegre RS, 1 o 4 de Outubro de 2017 ALGORITHM-AIDED IDENTIFICATION USING HISTORIC PROCESS DATA
ALGORITHM-AIDED IDENTIFICATION USING HISTORIC PROCESS DATA José Genario de Oliveira Jr, Claudio Garcia Departamento de Engenharia de Telecomunicações e Controle Escola Politécnica da Universidade de São
More informationSome Time-Series Models
Some Time-Series Models Outline 1. Stochastic processes and their properties 2. Stationary processes 3. Some properties of the autocorrelation function 4. Some useful models Purely random processes, random
More informationELEG 3143 Probability & Stochastic Process Ch. 6 Stochastic Process
Department of Electrical Engineering University of Arkansas ELEG 3143 Probability & Stochastic Process Ch. 6 Stochastic Process Dr. Jingxian Wu wuj@uark.edu OUTLINE 2 Definition of stochastic process (random
More informationESSE Mid-Term Test 2017 Tuesday 17 October :30-09:45
ESSE 4020 3.0 - Mid-Term Test 207 Tuesday 7 October 207. 08:30-09:45 Symbols have their usual meanings. All questions are worth 0 marks, although some are more difficult than others. Answer as many questions
More informationModel structure. Lecture Note #3 (Chap.6) Identification of time series model. ARMAX Models and Difference Equations
System Modeling and Identification Lecture ote #3 (Chap.6) CHBE 70 Korea University Prof. Dae Ryoo Yang Model structure ime series Multivariable time series x [ ] x x xm Multidimensional time series (temporal+spatial)
More informationS deals with the problem of building a mathematical model for
JOURNAL OF GUIDANCE, CONTROL, AND DYNAMICS Vol. 18, NO. 4, July-August 1995 Identification of Linear Stochastic Systems Through Projection Filters Chung-Wen Chen* Geophysical and Environmental Research
More informationTime Series Analysis. James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY
Time Series Analysis James D. Hamilton PRINCETON UNIVERSITY PRESS PRINCETON, NEW JERSEY PREFACE xiii 1 Difference Equations 1.1. First-Order Difference Equations 1 1.2. pth-order Difference Equations 7
More informationRegression with time series
Regression with time series Class Notes Manuel Arellano February 22, 2018 1 Classical regression model with time series Model and assumptions The basic assumption is E y t x 1,, x T = E y t x t = x tβ
More informationMatlab software tools for model identification and data analysis 11/12/2015 Prof. Marcello Farina
Matlab software tools for model identification and data analysis 11/12/2015 Prof. Marcello Farina Model Identification and Data Analysis (Academic year 2015-2016) Prof. Sergio Bittanti Outline Data generation
More informationTIME DELAY TEMPERATURE CONTROL WITH IMC AND CLOSED-LOOP IDENTIFICATION
TIME DELAY TEMPERATURE CONTROL WITH IMC AND CLOSED-LOOP IDENTIFICATION Naoto ABE, Katsuyoshi KARAKAWA and Hiroyuki ICHIHARA Department of Mechanical Engineering Informatics, Meiji University Kawasaki 4-857,
More informationChapter 4: Models for Stationary Time Series
Chapter 4: Models for Stationary Time Series Now we will introduce some useful parametric models for time series that are stationary processes. We begin by defining the General Linear Process. Let {Y t
More informationSAMPLE SOLUTION TO EXAM in MAS501 Control Systems 2 Autumn 2015
FACULTY OF ENGINEERING AND SCIENCE SAMPLE SOLUTION TO EXAM in MAS501 Control Systems 2 Autumn 2015 Lecturer: Michael Ruderman Problem 1: Frequency-domain analysis and control design (15 pt) Given is a
More informationState Estimation of Linear and Nonlinear Dynamic Systems
State Estimation of Linear and Nonlinear Dynamic Systems Part I: Linear Systems with Gaussian Noise James B. Rawlings and Fernando V. Lima Department of Chemical and Biological Engineering University of
More information