Matlab software tools for model identification and data analysis 10/11/2017 Prof. Marcello Farina
|
|
- Ashley Wood
- 5 years ago
- Views:
Transcription
1 Matlab software tools for model identification and data analysis 10/11/2017 Prof. Marcello Farina Model Identification and Data Analysis (Academic year ) Prof. Sergio Bittanti
2 Outline Data generation and analysis Definition and generation of time series and dynamic systems Analysis of systems Analysis of time series (sampled covariance and spectral analysis) Optimal predictors Identification Parameter identification Structure selection Model validation Exercizes
3 Time series and systems White noise >> N=100; % number of samples >> lambda=1; % standard deviation >> mu=1; % expected value I method >> u=mu+lambda*randn(n,1); II method >> u=wgn(m,n,p); % zero-mean signal M: number of signals P: signal power (dbw) (10log10(lambda^2)) Sum of a WGN and a generic signal >> u=awgn(x,snr,sigpower); X: original signal SRN: signal-to-noise ratio (db) SIGPOWER: power of signal X (dbw) (10log10( (0))) - if SIGPOWER= measured, AWGN estimates the power of signal X.
4 Time series and systems Definition of dynamic linear discrete-time systems: external representation (transfer function LTI format) >> system = tf(num,den,ts) num: numerator den: denominator Ts: sampling time (-1 default) >> system=tf([1 1],[1-0.8],1) Transfer function: z z - 0.8
5 Time series and systems Step response of a LTI system >> T=[0:Ts:Tfinal]; >> [y,t]=step(system,tfinal); % step response >> step(system,tfinal) % plot Response of a LTI system to a general input signal >> T=[0:Ts:Tfinal]; >> u1=sin(0.1*t); % e.g., sinusoidal input >> u2=1+randn(t,1); % e.g., white gaussian noise >> [y,t]=lsim(system,u1+u2,t);
6 Time series and systems The SYSTEM IDENTIFICATION TOOLBOX uses different formats for models (detto IDPOLY format) and signals (detto IDDATA format) Definition of IDDATA objects (i.e., inputs&outputs of a system in a single object): >> data=iddata(y,u,ts); u: input data y: output data Ts: sampling time
7 Time series and systems Definition of IDPOLY objects: >> model_idpoly = idpoly(a,b,c,d,f,noisevariance,ts) A,B,C,D,F are vectors (i.e., containing coefficients of the polynomials in incresing order of the exponent of z -1 ) Model_idpoly: A(q) y(t) = [B(q)/F(q)] u(t-nk) + [C(q)/D(q)] e(t) (If ARMAX: F=[], D=[]) The delay nk is defined in polynomial B(q). >> conversion from LTI format: Model_idpoly = idpoly(model_lti) model_lti is a LTI model (defined, e.g., using tf, zpk, ss) Response of an IDPOLY system to a generic input signal (vector u) >> y=sim(model_idpoly,u); -> It does not simulate the noise component, but just the deterministic one
8 Time series and systems Analysis of systems (both for LTI and IDPOLY formats): Zeros and poles: >> [Z,P,K]=zpkdata(model) Z = [-1] % positions of the zeros P = [0.8000] % positions of the poles K = 1 transfer constant Bode diagrams: >> [A,phi] = bode(model,w) W: values of the angular frequency w where to compute A(w) and phi(w) A: amplitude (A(w)) Phi: phase (phi(w))
9 Time series and systems Bode diagrams: >> [A,phi] = bode(model,[0.1 1]) A(:,:,1) = A(:,:,2) = phi(:,:,1) = phi(:,:,2) = >> bode(model) >> % draws the Bode plots (amplitude db- and phase)
10 Time series and systems Analysis of time series - m=mean(x) % Average or mean value. Y=detrend(X, constant ) % Remove constant trend from data sets. Y=detrend(X, linear ) % Remove linear trend from data sets. Covariance function ϒ( ) >> gamma=covf(y,10) % tau=0,1,,10-1 >> plot([0:9],gamma) Spectrum (!) periodogram method >> periodogram(y); % plots the periodogram >> [Pxx,w]=periodogram(y); >> plot(w/pi,10*log10(pxx)) Other methods are also available (e.g., fast fourier transform, see «help spectrum») Warning: covf computed the correlation function! TO OBTAIN THE COVARIANCE FUNCTION THE SIGNAL MUST HAVE ZERO MEAN Warning: no regularization technique is applied: the estimator is not consistent!
11 Predictors
12 Predictors
13 Predictors >> [ypred, x0_est, pred_model] = predict(model, data, k) INPUTS: model idpoly object data iddata object k prediction steps (integer number) OUTPUT ypred y_hat(t+1/t), predicted output, iddata object x0_est estimate of the initial condition (for state_space model) pred_model prediction model (discrete-time state-space model): >> pred_model_tf = tf(pred_model) % to obtain the transfer functions From input "y1" to output "y1": From input "u1" to output "y1 : It includes the prediction delay It includes the control delay
14 - EXAMPLE 1 (ARMAX model): Predictors >> Model_idpoly=idpoly([1 0.3],[ ],[1 0.5],[],[],1,1) Discrete-time IDPOLY model: A(q)y(t) = B(q)u(t) + C(q)e(t) A(q) = q^-1 B(q) = q^-1 C(q) = q^-1 tau = 1 % control delay >> [yd,xoe,modpred]=predict(model_idpoly,data,1); >> tf(modpred) From input "y1" to output "y1": 0.2 z^ z^-1 From input "u1" to output "y1": z^ z^ z^-1
15 Predictors - EXAMPLE 2 (OE model): >> Model_idpoly=idpoly([1 0.3],[ ],[1 0.3],[],[],1,1) Discrete-time IDPOLY model: A(q)y(t) = B(q)u(t) + C(q)e(t) A(q) = q^-1 B(q) = q^-1 C(q) = q^-1 tau = 1 % control delay >> [yd,xoe,modpred]=predict(model_idpoly,data,1); >> tf(modpred) From input "y1" to output "y1": 0 From input "u1" to output "y1": z^ z^ z^-1
16 Predictors - EXAMPLE 3 (ARX model): >> Model_idpoly=idpoly([1 0.3],[ ],[1],[],[],1,1) Discrete-time IDPOLY model: A(q)y(t) = B(q)u(t) + C(q)e(t) A(q) = q^-1 B(q) = q^-1 C(q) = 1 tau = 1 % control delay >> [yd,xoe,modpred]=predict(model_idpoly,data,1); >> tf(modpred) From input "y1" to output "y1": -0.3 z^-1 From input "u1" to output "y1": z^ z^-2
17 Prediction error Predictors >> Epsilon = pe(model_idpoly,data) % k-steps prediction error ARDERSON TEST FOR THE PREDICTION ERROR (99% confidence reguion) >> e=resid(model_idpoly,data) inputs model: IDPOLY format data: IDDATA format (input-output data) outputs e: IDDATA ([prediction error,input]) If the outputs are not specified, a plot is shown: - Empirical correlation function (empirical - sampled) of the prediction error - 99% confidence region -> if more than 99% of the samples of the normalized correlation function lie in the confidence region: positive results of the whiteness Test!
18 Example: Predictors >> A_idpoly=[1 0.3]; >> B_idpoly=[0 1 1/5]; %tau=1 >> C_idpoly=[1 1/2]; >> Model_idpoly = idpoly(a_idpoly,b_idpoly,c_idpoly,[],[],sigma,1); >>[yd,xoe,modpred]=predict(model_idpoly,data,1); >> resid(model_idpoly,data) >> data2=iddata(randn(n,2),1); >> resid(model_idpoly,data2) Correlation function of residuals. Output y Correlation function of residuals. Output y1 lag Cross corr. function between input u1 and residuals from output y
19 Identification >> systemidentification % start the MATLAB identification toolbox GUI
20 1.IMPORT DATA Identification
21 Identification 2.PREPROCESS DATA NEW (DETRENDED) DATA
22 3.ESTIMATE MODEL SELECT MODEL STRUCTURE Identification
23 Identification 3.1. MODEL STRUCTURES ARX(na,nb,k): A(z)y(t) = B(z)u(t) + e(t), where A(z)=1+ a_1 z^ a_na z^-na B(z)=b_1 z^-k b_nb z^-(k+nb-1) y(t)=-a_1 y(t-1)-...-a_na y(t-na) + b_1 u(t-k)+...+b_nb u(t-n_b-k+1)+e(t) ARMAX(na,nb,nc,k): A(z)y(t) = B(z)u(t) + C(z)e(t), where A(z)=1+ a_1 z^ a_na z^-na B(z)=b_1 z^-k b_nb z^-(k+nb-1) C(z)=1+ c_1 z^ c_nc z^-nc y(t)=-a_1 y(t-1)-...-a_na y(t-na) + b_1 u(t-k)+...+b_nb u(t-n_b-k+1) +e(t)+c_1 e(t-1) c_nc e(t-nc)
24 Identification 3.1. MODEL STRUCTURES OE(nb,nf,k): y(t) = B(z)/F(z)u(t) + e(t), where B(z)=b_1 z^-k b_nb z^-(k+nb-1) F(z)=1+ f_1 z^ f_nf z^-nf y(t)=-f_1 y(t-1)-...-f_na y(t-nf) + b_1 u(t-k)+...+b_nb u(t-n_b-k+1) +e(t)+f_1 e(t-1) f_nf e(t-nf) BOX-JENKINS(nb,nc,nd,nf,k): y(t) = B(z)/F(z)u(t) + C(z)/D(z)e(t) STATE-SPACE (n: order) [PEM or subspace identification] x(t+1)=ax(t)+bu(t)+ke(t) y(t)=cx(t)+du(t)+e(t)
25 Identification 3.2. EXTERNAL REPRESENTATION MODEL STRUCTURES BOX-JENKINS WHY ARE THESE CLASSES OF INTEREST? OUTPUT ERROR (OE) ARX ARMAX
26 Identification ARX(na,nb,k): y(t)=-a_1 y(t-1)-...-a_na y(t-na) + b_1 u(t-k)+...+b_nb u(t-n_b-k+1)+e(t) OPTIMAL 1-step PREDICTOR: yp(t-1)=-a_1 y(t-1)-...-a_na y(t-na) + b_1 u(t-k)+...+b_nb u(t-n_b-k+1) least square method (COMPUTATIONALLY LIGHTWEIGHT). IF THE SYSTEM DOES NOT BELONG TO THIS MODEL STRUCTURE: BAD RESULTS ( WRONG POLES AND ZEROS)! In our example, we estimate from data an ARX model [1 2 1] (with N=10000 data): Discrete-time IDPOLY model: A(q)y(t) = B(q)u(t) + e(t) A(q) = ( ) q^-1 B(q) = q^ ( ) q^-2 Estimated using ARX from data set edat Loss function and FPE Sampling interval: 1
27 OE(nb,nf,k): y(t) = B(z)/F(z)u(t) + e(t) OPTIMAL PREDICTOR: yp(t) = B(z)/F(z)u(t) Identification The OPTIMAL PREDICTOR is the SIMULATOR (i.e., it does not need data). The simulation with OE is also said Simulation Error Minimization. Computationally intensive It is able to well capture the dominant dynamics (poles and zeros) even if the model structure is not correct! In our example, we estimate from data a OE model [2 1 1] (with N=10000 data): Discrete-time IDPOLY model: y(t) = [B(q)/F(q)]u(t) + e(t) B(q) = q^ ( ) q^-2 F(q) = ( ) q^-1 Estimated using PEM using SearchMethod = Auto from data set z Loss function and FPE Sampling interval: 1
28 Identification In our example, we estimate from data a ARMAX model [ ] (with N=10000 data) i.e., the right model structure: Discrete-time IDPOLY model: A(q)y(t) = B(q)u(t) + C(q)e(t) A(q) = ( ) q^-1 B(q) = q^ ( ) q^-2 C(q) = ( ) q^-1 Estimated using PEM using SearchMethod = Auto from data set z Loss function and FPE Sampling interval: 1 Both with OE model class and with ARMAX model class (the structure of the data generator) we are able to correctly estimate the transfer function u->y
29 Identification 3.ESTIMATE MODEL IDENTIFICATION OF THE (CONTINUOUS-TIME) TRANSFER FUNCTION
30 Identification 3.ESTIMATE MODEL IDENTIFICATION OF MODELS: -NONLINEAR ARX - HAMMERSTEIN WIENER
31 4.Compare and validate models Identification
32 Identification For further details, see
33 Identification Other related Matlab functions (2 Identification toolbox) Parametric model estimation ar - AR-models of signals using various approaches. armax - Prediction error estimate of an ARMAX model. arx - LS-estimate of ARX-models. bj - Prediction error estimate of a Box-Jenkins model. ivar - IV-estimates for the AR-part of a scalar time series (with instrumental variable method). iv4 - Approximately optimal IV-estimates for ARX-models. n4sid - State-space model estimation using a sub-space method. nlarx - Prediction error estimate of a nonlinear ARX model. nlhw - Prediction error estimate of a Hammerstein-Wiener model. oe - Prediction error estimate of an output-error model. pem - Prediction error estimate of a general model.
34 Identification Model structure selection aic - Compute Akaike's information criterion. fpe - Compute final prediction error criterion. arxstruc - Loss functions for families of ARX-models. selstruc - Select model structures according to various criteria. idss/setstruc - Set the structure matrices for idss objects. struc - Generate typical structure matrices for ARXSTRUC. For more details please see: >> help ident
35 Exercizes
36 Exercizes
37 Exercizes
38 Exercizes
39 Exercizes Model Identification Consider the ARMAX system y(t)=0.1 y(t-1)+2u(t-1)-0.8 u(t-2)+e(t)+0.3e(t-1) Where e(t) is a white noise with zero mean and unitary variance. 1. Generate a 1000-samples realization of the output signal y(t), when the input u(t) is a persistently exciting signal (e.g., a white noise), and where y(0)=0;. 2. Identify, using the identification toolbox: 1. an ARX(1,1,1) model 2. an ARX(1,2,1) model 3. an ARMAX(1,2,2,1) model 4. an ARMAX(1,1,2,1) model 5. an OE(2,1,1) model 3. Compare and analyze the results (possibly add validation data to validate the results).
40 Exercizes
41 Exercizes
42 Exercizes
43 Exercizes
Matlab 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 informationComputer Exercise 1 Estimation and Model Validation
Lund University Time Series Analysis Mathematical Statistics Fall 2018 Centre for Mathematical Sciences Computer Exercise 1 Estimation and Model Validation This computer exercise treats identification,
More informationIdentification of a best thermal formula and model for oil and winding 0f power transformers using prediction methods
Identification of a best thermal formula and model for and 0f power transformers using prediction methods Kourosh Mousavi Takami Jafar Mahmoudi Mälardalen University, Västerås, sweden Kourosh.mousavi.takami@mdh.se
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 informationA Mathematica Toolbox for Signals, Models and Identification
The International Federation of Automatic Control A Mathematica Toolbox for Signals, Models and Identification Håkan Hjalmarsson Jonas Sjöberg ACCESS Linnaeus Center, Electrical Engineering, KTH Royal
More informationEL1820 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 informationModel Identification and Validation for a Heating System using MATLAB System Identification Toolbox
IOP Conference Series: Materials Science and Engineering OPEN ACCESS Model Identification and Validation for a Heating System using MATLAB System Identification Toolbox To cite this article: Muhammad Junaid
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 informationData-Driven Modelling of Natural Gas Dehydrators for Dew Point Determination
International Journal of Oil, Gas and Coal Engineering 2017; 5(3): 27-33 http://www.sciencepublishinggroup.com/j/ogce doi: 10.11648/j.ogce.20170503.11 ISSN: 2376-7669(Print); ISSN: 2376-7677(Online) Data-Driven
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 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 information12. Prediction Error Methods (PEM)
12. Prediction Error Methods (PEM) EE531 (Semester II, 2010) description optimal prediction Kalman filter statistical results computational aspects 12-1 Description idea: determine the model parameter
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 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 informationOutline 2(42) Sysid Course VT An Overview. Data from Gripen 4(42) An Introductory Example 2,530 3(42)
Outline 2(42) Sysid Course T1 2016 An Overview. Automatic Control, SY, Linköpings Universitet An Umbrella Contribution for the aterial in the Course The classic, conventional System dentification Setup
More informationTHERE are two types of configurations [1] in the
Linear Identification of a Steam Generation Plant Magdi S. Mahmoud Member, IAEG Abstract The paper examines the development of models of steam generation plant using linear identification techniques. The
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 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 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 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 informationStudy of Time Series and Development of System Identification Model for Agarwada Raingauge Station
Study of Time Series and Development of System Identification Model for Agarwada Raingauge Station N.A. Bhatia 1 and T.M.V.Suryanarayana 2 1 Teaching Assistant, 2 Assistant Professor, Water Resources Engineering
More informationRozwiązanie zagadnienia odwrotnego wyznaczania sił obciąŝających konstrukcje w czasie eksploatacji
Rozwiązanie zagadnienia odwrotnego wyznaczania sił obciąŝających konstrukcje w czasie eksploatacji Tadeusz Uhl Piotr Czop Krzysztof Mendrok Faculty of Mechanical Engineering and Robotics Department of
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 informationChap 4. State-Space Solutions and
Chap 4. State-Space Solutions and Realizations Outlines 1. Introduction 2. Solution of LTI State Equation 3. Equivalent State Equations 4. Realizations 5. Solution of Linear Time-Varying (LTV) Equations
More informationMODEL PREDICTIVE CONTROL and optimization
MODEL PREDICTIVE CONTROL and optimization Lecture notes Model Predictive Control PhD., Associate professor David Di Ruscio System and Control Engineering Department of Technology Telemark University College
More informationForecasting of ATM cash demand
Forecasting of ATM cash demand Ayush Maheshwari Mathematics and scientific computing Department of Mathematics and Statistics Indian Institute of Technology (IIT), Kanpur ayushiitk@gmail.com Project guide:
More informationOutline. What Can Regularization Offer for Estimation of Dynamical Systems? State-of-the-Art System Identification
Outline What Can Regularization Offer for Estimation of Dynamical Systems? with Tianshi Chen Preamble: The classic, conventional System Identification Setup Bias Variance, Model Size Selection Regularization
More informationIntroduction to system identification
Introduction to system identification Jan Swevers July 2006 0-0 Introduction to system identification 1 Contents of this lecture What is system identification Time vs. frequency domain identification Discrete
More informationModeling of Surface EMG Signals using System Identification Techniques
Modeling of Surface EMG Signals using System Identification Techniques Vishnu R S PG Scholar, Dept. of Electrical and Electronics Engg. Mar Baselios College of Engineering and Technology Thiruvananthapuram,
More informationEECE Adaptive Control
EECE 574 - Adaptive Control Recursive Identification Algorithms Guy Dumont Department of Electrical and Computer Engineering University of British Columbia January 2012 Guy Dumont (UBC EECE) EECE 574 -
More informationCentre for Mathematical Sciences HT 2017 Mathematical Statistics
Lund University Stationary stochastic processes Centre for Mathematical Sciences HT 2017 Mathematical Statistics Computer exercise 3 in Stationary stochastic processes, HT 17. The purpose of this exercise
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 informationCONTINUOUS TIME D=0 ZOH D 0 D=0 FOH D 0
IDENTIFICATION ASPECTS OF INTER- SAMPLE INPUT BEHAVIOR T. ANDERSSON, P. PUCAR and L. LJUNG University of Linkoping, Department of Electrical Engineering, S-581 83 Linkoping, Sweden Abstract. In this contribution
More informationData-driven methods in application to flood defence systems monitoring and analysis Pyayt, A.
UvA-DARE (Digital Academic Repository) Data-driven methods in application to flood defence systems monitoring and analysis Pyayt, A. Link to publication Citation for published version (APA): Pyayt, A.
More informationPractical Spectral Estimation
Digital Signal Processing/F.G. Meyer Lecture 4 Copyright 2015 François G. Meyer. All Rights Reserved. Practical Spectral Estimation 1 Introduction The goal of spectral estimation is to estimate how the
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 informationÜbersetzungshilfe / Translation aid (English) To be returned at the end of the exam!
Prüfung Regelungstechnik I (Control Systems I) Prof. Dr. Lino Guzzella 3.. 24 Übersetzungshilfe / Translation aid (English) To be returned at the end of the exam! Do not mark up this translation aid -
More informationJ. Sjöberg et al. (1995):Non-linear Black-Box Modeling in System Identification: a Unified Overview, Automatica, Vol. 31, 12, str
Dynamic Systems Identification Part - Nonlinear systems Reference: J. Sjöberg et al. (995):Non-linear Black-Box Modeling in System Identification: a Unified Overview, Automatica, Vol. 3,, str. 69-74. Systems
More informationResearch Article Experimental Parametric Identification of a Flexible Beam Using Piezoelectric Sensors and Actuators
Shock and Vibration, Article ID 71814, 5 pages http://dx.doi.org/1.1155/214/71814 Research Article Experimental Parametric Identification of a Flexible Beam Using Piezoelectric Sensors and Actuators Sajad
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 informationTopics in Undergraduate Control Systems Design
Topics in Undergraduate Control Systems Design João P. Hespanha April 9, 2006 Disclaimer: This is an early draft and probably contains many typos. Comments and information about typos are welcome. Please
More informationAdvanced Process Control Tutorial Problem Set 2 Development of Control Relevant Models through System Identification
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
More informationUnivariate ARIMA Models
Univariate ARIMA Models ARIMA Model Building Steps: Identification: Using graphs, statistics, ACFs and PACFs, transformations, etc. to achieve stationary and tentatively identify patterns and model components.
More informationSGN Advanced Signal Processing: Lecture 8 Parameter estimation for AR and MA models. Model order selection
SG 21006 Advanced Signal Processing: Lecture 8 Parameter estimation for AR and MA models. Model order selection Ioan Tabus Department of Signal Processing Tampere University of Technology Finland 1 / 28
More informationSound Listener s perception
Inversion of Loudspeaker Dynamics by Polynomial LQ Feedforward Control Mikael Sternad, Mathias Johansson and Jonas Rutstrom Abstract- Loudspeakers always introduce linear and nonlinear distortions in a
More informationA6523 Modeling, Inference, and Mining Jim Cordes, Cornell University
A6523 Modeling, Inference, and Mining Jim Cordes, Cornell University Lecture 19 Modeling Topics plan: Modeling (linear/non- linear least squares) Bayesian inference Bayesian approaches to spectral esbmabon;
More informationECE 636: Systems identification
ECE 636: Systems identification Lectures 7 8 onparametric identification (continued) Important distributions: chi square, t distribution, F distribution Sampling distributions ib i Sample mean If the variance
More informationf-domain expression for the limit model Combine: 5.12 Approximate Modelling What can be said about H(q, θ) G(q, θ ) H(q, θ ) with
5.2 Approximate Modelling What can be said about if S / M, and even G / G? G(q, ) H(q, ) f-domain expression for the limit model Combine: with ε(t, ) =H(q, ) [y(t) G(q, )u(t)] y(t) =G (q)u(t) v(t) We know
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 informationWhat can regularization offer for estimation of dynamical systems?
8 6 4 6 What can regularization offer for estimation of dynamical systems? Lennart Ljung Tianshi Chen Division of Automatic Control, Department of Electrical Engineering, Linköping University, SE-58 83
More information6. Methods for Rational Spectra It is assumed that signals have rational spectra m k= m
6. Methods for Rational Spectra It is assumed that signals have rational spectra m k= m φ(ω) = γ ke jωk n k= n ρ, (23) jωk ke where γ k = γ k and ρ k = ρ k. Any continuous PSD can be approximated arbitrary
More informationLecture Note #7 (Chap.11)
System Modeling and Identification Lecture Note #7 (Chap.) CBE 702 Korea University Prof. Dae Ryoo Yang Chap. Real-time Identification Real-time identification Supervision and tracing of time varying parameters
More informationLTI Approximations of Slightly Nonlinear Systems: Some Intriguing Examples
LTI Approximations of Slightly Nonlinear Systems: Some Intriguing Examples Martin Enqvist, Lennart Ljung Division of Automatic Control Department of Electrical Engineering Linköpings universitet, SE-581
More informationMatlab Controller Design. 1. Control system toolbox 2. Functions for model analysis 3. Linear system simulation 4. Biochemical reactor linearization
Matlab Controller Design. Control system toolbox 2. Functions for model analysis 3. Linear system simulation 4. Biochemical reactor linearization Control System Toolbox Provides algorithms and tools for
More informationGENERALIZED MODEL OF THE PRESSURE EVOLUTION WITHIN THE BRAKING SYSTEM OF A VEHICLE
1. Marin MARINESCU,. Radu VILAU GENERALIZED MODEL OF THE PRESSURE EVOLUTION WITHIN THE BRAKING SYSTEM OF A VEHICLE 1-. MILITARY TECHNICAL ACADEMY, BUCHAREST, ROMANIA ABSTRACT: There are lots of methods
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 informationContinuous Order Identification of PHWR Models Under Step-back for the Design of Hyper-damped Power Tracking Controller with Enhanced Reactor Safety
Continuous Order Identification of PHWR Models Under Step-back for the Design of Hyper-damped Power Tracking Controller with Enhanced Reactor Safety Saptarshi Das a,b, Sumit Mukherjee b,c, Shantanu Das
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 informationModeling and Predicting Chaotic Time Series
Chapter 14 Modeling and Predicting Chaotic Time Series To understand the behavior of a dynamical system in terms of some meaningful parameters we seek the appropriate mathematical model that captures the
More informationThermal error compensation for a high precision lathe
IWMF214, 9 th INTERNATIONAL WORKSHOP ON MICROFACTORIES OCTOBER 5-8, 214, HONOLULU, U.S.A. / 1 Thermal error compensation for a high precision lathe Byung-Sub Kim # and Jong-Kweon Park Department of Ultra
More informationCHAPTER 2 RANDOM PROCESSES IN DISCRETE TIME
CHAPTER 2 RANDOM PROCESSES IN DISCRETE TIME Shri Mata Vaishno Devi University, (SMVDU), 2013 Page 13 CHAPTER 2 RANDOM PROCESSES IN DISCRETE TIME When characterizing or modeling a random variable, estimates
More informationForward and Inverse Modelling Approaches for Prediction of Light Stimulus from Electrophysiological Response in Plants
Forward and Inverse Modelling Approaches for Prediction of Light Stimulus from Electrophysiological Response in Plants Shre Kumar Chatterjee a, Sanmitra Ghosh a, Saptarshi Das a, *, Veronica Manzella b,c,
More informationModeling and Identification of Dynamic Systems (vimmd312, 2018)
Modeling and Identification of Dynamic Systems (vimmd312, 2018) Textbook background of the curriculum taken. In parenthesis: material not reviewed due to time shortage, but which is suggested to be read
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 informationClassic Time Series Analysis
Classic Time Series Analysis Concepts and Definitions Let Y be a random number with PDF f Y t ~f,t Define t =E[Y t ] m(t) is known as the trend Define the autocovariance t, s =COV [Y t,y s ] =E[ Y t t
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 informationBounding the parameters of linear systems with stability constraints
010 American Control Conference Marriott Waterfront, Baltimore, MD, USA June 30-July 0, 010 WeC17 Bounding the parameters of linear systems with stability constraints V Cerone, D Piga, D Regruto Abstract
More informationFrequency methods for the analysis of feedback systems. Lecture 6. Loop analysis of feedback systems. Nyquist approach to study stability
Lecture 6. Loop analysis of feedback systems 1. Motivation 2. Graphical representation of frequency response: Bode and Nyquist curves 3. Nyquist stability theorem 4. Stability margins Frequency methods
More informationComparison of different seat-to-head transfer functions for vibrational comfort monitoring of car passengers
Comparison of different seat-to-head transfer functions for vibrational comfort monitoring of car passengers Daniele Carnevale 1, Ettore Pennestrì 2, Pier Paolo Valentini 2, Fabrizio Scirè Ingastone 2,
More informationTutorial lecture 2: System identification
Tutorial lecture 2: System identification Data driven modeling: Find a good model from noisy data. Model class: Set of all a priori feasible candidate systems Identification procedure: Attach a system
More informationTransfer function for the seismic signal recorded in solid and fractured rock surrounding deep level mining excavations
Transfer function for the seismic signal recorded in solid and fractured rock surrounding deep level mining excavations by A. Cichowicz*, A.M. Milev, and R.J. Durrheim Synopsis The obective of this work
More informationA Log-Frequency Approach to the Identification of the Wiener-Hammerstein Model
A Log-Frequency Approach to the Identification of the Wiener-Hammerstein Model The MIT Faculty has made this article openly available Please share how this access benefits you Your story matters Citation
More informationUNIVERSITY OF CALGARY. Nonlinear Closed-Loop System Identification In The Presence Of. Non-stationary Noise Sources. Ibrahim Aljamaan A THESIS
UNIVERSITY OF CALGARY Nonlinear Closed-Loop System Identification In The Presence Of Non-stationary Noise Sources by Ibrahim Aljamaan A THESIS SUBMITTED TO THE FACULTY OF GRADUATE STUDIES IN PARTIAL FULFILMENT
More informationIdentification, Model Validation and Control. Lennart Ljung, Linköping
Identification, Model Validation and Control Lennart Ljung, Linköping Acknowledgment: Useful discussions with U Forssell and H Hjalmarsson 1 Outline 1. Introduction 2. System Identification (in closed
More informationFinal Exam January 31, Solutions
Final Exam January 31, 014 Signals & Systems (151-0575-01) Prof. R. D Andrea & P. Reist Solutions Exam Duration: Number of Problems: Total Points: Permitted aids: Important: 150 minutes 7 problems 50 points
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 information5 Autoregressive-Moving-Average Modeling
5 Autoregressive-Moving-Average Modeling 5. Purpose. Autoregressive-moving-average (ARMA models are mathematical models of the persistence, or autocorrelation, in a time series. ARMA models are widely
More informationEE 602 TERM PAPER PRESENTATION Richa Tripathi Mounika Boppudi FOURIER SERIES BASED MODEL FOR STATISTICAL SIGNAL PROCESSING - CHONG YUNG CHI
EE 602 TERM PAPER PRESENTATION Richa Tripathi Mounika Boppudi FOURIER SERIES BASED MODEL FOR STATISTICAL SIGNAL PROCESSING - CHONG YUNG CHI ABSTRACT The goal of the paper is to present a parametric Fourier
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 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 informationLecture 6: Deterministic Self-Tuning Regulators
Lecture 6: Deterministic Self-Tuning Regulators Feedback Control Design for Nominal Plant Model via Pole Placement Indirect Self-Tuning Regulators Direct Self-Tuning Regulators c Anton Shiriaev. May, 2007.
More informationClosed-Loop Identification of Fractionalorder Models using FOMCON Toolbox for MATLAB
Closed-Loop Identification of Fractionalorder Models using FOMCON Toolbox for MATLAB Aleksei Tepljakov, Eduard Petlenkov, Juri Belikov October 7, 2014 Motivation, Contribution, and Outline The problem
More informationLateral dynamics of a SUV on deformable surfaces by system identification. Part II. Models reconstruction
IOP Conference Series: Materials Science and Engineering PAPER OPEN ACCESS Lateral dynamics of a SUV on deformable surfaces by system identification. Part II. Models reconstruction To cite this article:
More information3. ESTIMATION OF SIGNALS USING A LEAST SQUARES TECHNIQUE
3. ESTIMATION OF SIGNALS USING A LEAST SQUARES TECHNIQUE 3.0 INTRODUCTION The purpose of this chapter is to introduce estimators shortly. More elaborated courses on System Identification, which are given
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 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 informationcovariance function, 174 probability structure of; Yule-Walker equations, 174 Moving average process, fluctuations, 5-6, 175 probability structure of
Index* The Statistical Analysis of Time Series by T. W. Anderson Copyright 1971 John Wiley & Sons, Inc. Aliasing, 387-388 Autoregressive {continued) Amplitude, 4, 94 case of first-order, 174 Associated
More information17. Example SAS Commands for Analysis of a Classic Split-Plot Experiment 17. 1
17 Example SAS Commands for Analysis of a Classic SplitPlot Experiment 17 1 DELIMITED options nocenter nonumber nodate ls80; Format SCREEN OUTPUT proc import datafile"c:\data\simulatedsplitplotdatatxt"
More informationResampling techniques for statistical modeling
Resampling techniques for statistical modeling Gianluca Bontempi Département d Informatique Boulevard de Triomphe - CP 212 http://www.ulb.ac.be/di Resampling techniques p.1/33 Beyond the empirical error
More informationSTABILITY OF CLOSED-LOOP CONTOL SYSTEMS
CHBE320 LECTURE X STABILITY OF CLOSED-LOOP CONTOL SYSTEMS Professor Dae Ryook Yang Spring 2018 Dept. of Chemical and Biological Engineering 10-1 Road Map of the Lecture X Stability of closed-loop control
More informationIdentification of MIMO linear models: introduction to subspace methods
Identification of MIMO linear models: introduction to subspace methods Marco Lovera Dipartimento di Scienze e Tecnologie Aerospaziali Politecnico di Milano marco.lovera@polimi.it State space identification
More information9. Model Selection. statistical models. overview of model selection. information criteria. goodness-of-fit measures
FE661 - Statistical Methods for Financial Engineering 9. Model Selection Jitkomut Songsiri statistical models overview of model selection information criteria goodness-of-fit measures 9-1 Statistical models
More informationStatistics of Stochastic Processes
Prof. Dr. J. Franke All of Statistics 4.1 Statistics of Stochastic Processes discrete time: sequence of r.v...., X 1, X 0, X 1, X 2,... X t R d in general. Here: d = 1. continuous time: random function
More informationControl Systems I. Lecture 6: Poles and Zeros. Readings: Emilio Frazzoli. Institute for Dynamic Systems and Control D-MAVT ETH Zürich
Control Systems I Lecture 6: Poles and Zeros Readings: Emilio Frazzoli Institute for Dynamic Systems and Control D-MAVT ETH Zürich October 27, 2017 E. Frazzoli (ETH) Lecture 6: Control Systems I 27/10/2017
More informationPiezoelectric Transducer Modeling: with System Identification (SI) Method
World Academy of Science, Engineering and Technology 5 8 Piezoelectric Transducer Modeling: with System Identification (SI) Method Nora Taghavi, and Ali Sadr Abstract System identification is the process
More informationSystem identification - L17-19
Stochastic Adaptive Control (02421) www.imm.dtu.dk/courses/02421 Niels Kjølstad Poulsen Build. 303B, room 016 Section for Dynamical Systems Dept. of Applied Mathematics and Computer Science The Technical
More informationSPEECH ANALYSIS AND SYNTHESIS
16 Chapter 2 SPEECH ANALYSIS AND SYNTHESIS 2.1 INTRODUCTION: Speech signal analysis is used to characterize the spectral information of an input speech signal. Speech signal analysis [52-53] techniques
More informationPart III Spectrum Estimation
ECE79-4 Part III Part III Spectrum Estimation 3. Parametric Methods for Spectral Estimation Electrical & Computer Engineering North Carolina State University Acnowledgment: ECE79-4 slides were adapted
More informationChapter 18 The STATESPACE Procedure
Chapter 18 The STATESPACE Procedure Chapter Table of Contents OVERVIEW...999 GETTING STARTED...1002 AutomaticStateSpaceModelSelection...1003 SpecifyingtheStateSpaceModel...1010 SYNTAX...1013 FunctionalSummary...1013
More informationDirect continuous-time approaches to system identification. Overview and benefits for practical applications
Direct continuous-time approaches to system identification. Overview and benefits for practical applications Hugues GARNIER University of Lorraine, CNRS Research Centre for Automatic Control (CRAN Nancy,
More information