Hammerstein System Identification by a Semi-Parametric Method
|
|
- Percival Garrett
- 5 years ago
- Views:
Transcription
1 Hammerstein System Identification by a Semi-arametric ethod Grzegorz zy # # Institute of Engineering Cybernetics, Wroclaw University of echnology, ul. Janiszewsiego 11/17, Wroclaw, oland, grmz@ict.pwr.wroc.pl ABSAC A semi-parametric algorithm for identification of Hammerstein systems in the presence of correlated noise is proposed. he procedure is based on the non-parametric ernel regression estimator and the standard least squares. he advantages of the method in comparison with the standard non-parametric approach are discussed. Limit properties of the proposed estimator are studied, and the simulation results are presented. Keywords: onlinear systems, non-parametric regression, system identification 1. IODUCIO he nonlinear, dynamic system modeling has been widely discussed in the literature for the last three decades. umerous applications (in signal processing, telecommunication systems, designing of adaptive filters, pattern recognition, automation, biological systems (see [1]-[3])) show necessity of the research of more effective identification algorithms. he oldest identification techniques based on the Volterra and Wiener functional expansions gave very complicated algorithms. In 8 s the another approach appeared, basing on the assumption that the nonlinear dynamic system can be represented by interconnections of linear dynamic and static nonlinear elements [4]. he most of authors proposed the parametric algorithms without rigorous convergence proofs, restricted to the non-linearity of a polynomial form only. on-parametric methods were first applied by Greblici and awla ([5]-[8]) and based on the estimation of the regression function (ernel estimator, orthogonal series estimator, and recently on the wavelet approach [9]). Convergence of the proposed estimators was proved with the poor restrictions on the nonlinearity. he main disadvantage of these methods is non-analytic form of the solution. An idea of our algorithm is similar to that in [1], where the parameter estimator for the simple, static, linear element constructed by non-parametric estimation was analyzed. First the evaluations of the unmeasureable signal are computed by a non-parametric method, to be exploited in the parametric (least-squares) estimation of the unnown non-linearity. In section 2 the studied Hammerstein system (Fig. 1) is presented, and the identification problem is defined. ext, a two-stage identification algorithm is proposed in section 3, and its theoretical limit properties are analyzed in section 4. Simulation examples for the sample nonlinear system are shown in section 5.
2 2. SAEE OF HE OBLE 2.1. System under consideration he system consists of a static non-linearity µ () followed by a linear FI filter with the impulse response coefficients { λ p} p =. he system output is disturbed by a random, correlated noise process z. Only the input x and the output y of the overall system are accessible for measurements. x ε { ω r } w µ ( ) { λ p } p = r = v z y Fig.1. Hammerstein system with the correlated noise process. he system is described by the equations w which yelds ( x ) = µ, v = λ w p p y, z = ω r ε r ( x ) p p r= r= = λ µ + ω ε r r, y = v + z (1) 2.2. Assumptions and comments We assume the following: (A1) the random input sequence { x } (A2) the disturbance { } ε = = is iid, and is bounded: x x < is independent of inputs, zero mean, bounded, stationary white noise, i.e. Εε =, Εε ε j = for j, and ε < ε < (A3) the non-linearity () µ x = ϕ x a *, where µ x is of the form () () = ( a a a S ) () = ( 1(), 2(),..., S ()), where fs () x a * 1, 2,..., is an unnown true parameter vector, and ϕ x f x f x f x, s=1...s is a set of a priori nown, linearly independent functions, such that fs() x f < for x x (cf. (A1))
3 (A4) the correlated noise process { } filter, and ω r r = < z = (A() stationary process) is generated by a linear finite impulse response (A5) for clarity of exposition we additionally assume that µ () = and λ Under (A1)-(A5) we have: (C1) from (A2) and (A4) we conclude that ( z z <, where z = ε ω r r = ) { } z is zero-mean and bounded (C2) from (A1) and (A3) the vector process ϕ( x ) = Denoting W = ( w1, w2,..., w) and = ( ϕ( x1) ϕ( x2) ϕ( x )) W =Φ a *. From (C2) one can infer that the matrix Φ Φ a * = ( ) Φ Φ Φ is persistently exciting of order S,,..., we get Φ W is nonsingular, and 2.3. Identification tas he aim of the identification is to recover the unnown parameter vector a * of the nonlinearity using the input-output measurements {( x, y )} =1 (2) from the whole system, obtained in the experiment. We emphasize that the non-linearity output w cannot be directly measured. 3. SEI-AAEIC ALGOIH Observe that (cf. (1)) ( ) ( ) y = λ µ x + λ µ x + z p p 1 and () () Ε () x = λ µ x + µ x λ p Hence, including assumption (A5), we get () () µ () x = () x () where x () = { y x = x} 1 = Εµ x λ and further p Ε is the regression function. Basing on this observation we propose the following identification algorithm: { } Stage 1. Exploiting ( x, y ) { } w =1 =1 1 estimate the values of the unmeasureable signal by the non-parametric method (ernel regression estimate): µ ( ) ( ) (), where ( ) w! =! x =! x!,! x yk x x i = i i h K x x i ( ) i h( ) and K () is a ernel function, h( ) is a smoothing parameter (see[5]). (3)
4 Stage 2. Using the pairs {( x w ) },!, =1 a where W! ( w!,, w!,,..., w!, ) = 1 2. = compute the least-squares parameter estimate: ( ) Φ Φ 1 Φ! (4) W 4. LII OEIES Introducing the error of the non-parametric estimator as η = w! w =! x! w =! µ x µ x ( ( ) ()) ( ) ( ),, the parametric part of the algorithm can be interpreted as the identification of a non-linear, static element µ () x ϕ () x = a * corrupted by the disturbance η, (Fig. 2) from the noisy input-output measurements {( x w ) },!, =1. x!, µ ( ) w η, w Fig.2. Illustration of the second stage of the algorithm. Denoting Θ = ( η1, η2, η, ) (2) and (4)): a ( ) W ( ),,..., we obtain W! 1 = W + Θ which yields (from = Φ Φ Φ + Φ Φ Φ Θ = a + *, where = 1 1 Φ Φ Φ Θ ( ) = 1 1 = ϕ ϕ ϕ η 1, 2, he matrix 1, is invertible (see (C2)) and has finite elements (see (A3)), so (, ) well defined. aing use of the results in [5], we realize that w!, η, provided that () in probability, i.e. lim (, ) 1 w is and η > δ = for each δ > and each =1..., K fulfils standard conditions for ernel functions [5] and ( ) appropriately selected, such that h( ) and h( ) heorem 1. Assume that (A1)-(A6) hold and ( η δ). h is α, > c for some α > (in τ stage 1). Let = and τ >. If ατ > 1 then a a * in probability as.
5 roof. Each element 2, [] s of the vector 2, fulfils the condition 1 1 f = η η η, = UB [] s f ( x ) f ( x ) 2, s, s, UB ( > f ) = > δ η δ η > δ,, α ( η > δ) = τ ( ) α 1 ατ, c c = c If ατ > 1 then 1 ατ as. hus UB * a s a s in probability for each s=1...s. [] [] in probability, and 5. EXEIEAL ESULS he results of the experiments under the sample nonlinear system are presented in Fig. 3. It illustrates consistency property and shows the solutions for =1 and for =5 measurements. 3 E a - a * 2 E()= % a * y =1 =21 x ~ U(-1,1) ε ~ U(-2,2) µ (x) System specification: y () = (,,sin ) 2 ( ) ϕ x x x x a * = 21,, 2 = 2, λ, λ 1 = 1, λ 2 = 2, ω, ω 1, ω 2 Used ernel function Ku () h( ) = 5. τ 5. =5 =63 x ~ U(-1,1) ~ U(-2,2) ε = e u 2 µ (x) ϕ (x)*a 21 measurements x -4-8 ϕ (x)*a 63 measurements Fig.3. Experimental results for a sample nonlinear system identification x
6 6. COCLUSIOS he proposed method has the following advantages in comparison to the standard nonparametric approach: the resulting model has the analytic form, which is very important in practice we can estimate µ ( x ) not having any measurements from neighborhood of x However certain a priori nowledge about the non-linearity is needed. amely, the nonlinearity has to be nown with accuracy to the parameters but is not assumed to be of a polynomial type, which was commonly assumed earlier (eg. [4]). o increase numerical efficiency, both stages of the algorithm can be computed in recursive fashion. he recursive version of (3) in stage 1 can be found in [8]. he recursive version of (4) in stage 2 is the standard recursive least-squares algorithm. Experimental examinations over the systems containing non-typical non-linearities confirmed good behavior of the algorithm. he algorithm is simple to compute and can be well applied in many fields. EFEECES [1] W.X. Zheng, On a Least-Squares-Based Algorithm for Identification of Stochastic Linear Systems, IEEE ransactions on Signal rocessing, Vol. 46, o. 6, 1998, pp [2] G.Giunta, G. Jacovitti, A. eri, Bandpass onlinear System Identification by Higher Cross Correlation, IEEE ransactions on Signal rocessing, 29, 1991, pp [3] S.S. Hayin, Adaptive filter theory (third edition), rentice-hall, [4] S.A. Billings, S.Y.Fahouri, Identification of Systems Containing Linear Dynamic and Static onlinear Elements, Automatica, Vol. 18, o. 1, 1982, pp [5] W. Greblici,. awla, Hammerstein System Identification by on-parametric egression Estimation, International Journal of Control, Vol. 45, o. 1, 1987, pp [6] W. Greblici,. awla, onparametric Identification of Hammerstein Systems, IEEE ransactions on Information heory, Vol. 35, o. 2, 1989, pp [7] W. Greblici,. awla, Cascade on-linear System Identification by a on-arametric ethod, International Journal of System Science, Vol. 25, o. 1, 1994, pp [8] W. Greblici,. awla, ecursive onparametric Identification of Hammerstein Systems, J.Franlin Inst., Vol. 326, o. 4, [9]. awla, Z. Hasiewicz, onlinear System Identification by the Haar ultiresolution Analysis, IEEE ransactions on Circuits and Systems, Vol. 45, o. 9, 1998, pp [1] J.A. Cristobal Cristobal,. Feraldo oca, W. Gonzales anteiga, A Class of Linear egression arameter Estimators Constructed by onparametric Estimation, he Annals of Statistics, Vol. 15, o. 2, 1987, pp
APPLICATION OF INSTRUMENTAL VARIABLE METHOD TO THE IDENTIFICATION OF HAMMERSTEIN-WIENER SYSTEMS
APPLICAIO OF ISRUMEAL VARIABLE MEHOD O HE IDEIFICAIO OF HAMMERSEI-WIEER SYSEMS GREGOR MYK Wrocław University of echnology, Institute of Engineering Cybernetics, ul. Janisewsiego /7, 57 Wrocław, Poland,
More informationChapter 2 Hammerstein Systems
Chapter 2 Hammerstein Systems Abstract A discrete-time cascade of a static nonlinearity followed by a linear dynamics, i.e. the Hammerstein system, is presented. Its equivalence (from the proposed identification
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 informationFrom Volterra series through block-oriented approach to Volterra series
From Volterra series through block-oriented approach to Volterra series Pawel Wachel and Przemyslaw Sliwinski Department of Control Systems and Mechatronics Wroclaw University of Technology European Research
More informationCramér-Rao Bounds for Estimation of Linear System Noise Covariances
Journal of Mechanical Engineering and Automation (): 6- DOI: 593/jjmea Cramér-Rao Bounds for Estimation of Linear System oise Covariances Peter Matiso * Vladimír Havlena Czech echnical University in Prague
More informationAn Iterative Algorithm for the Subspace Identification of SISO Hammerstein Systems
An Iterative Algorithm for the Subspace Identification of SISO Hammerstein Systems Kian Jalaleddini R. E. Kearney Department of Biomedical Engineering, McGill University, 3775 University, Montréal, Québec
More informationExpressions for the covariance matrix of covariance data
Expressions for the covariance matrix of covariance data Torsten Söderström Division of Systems and Control, Department of Information Technology, Uppsala University, P O Box 337, SE-7505 Uppsala, Sweden
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 informationELEC system identification workshop. Subspace methods
1 / 33 ELEC system identification workshop Subspace methods Ivan Markovsky 2 / 33 Plan 1. Behavioral approach 2. Subspace methods 3. Optimization methods 3 / 33 Outline Exact modeling Algorithms 4 / 33
More informationNONLINEAR ECHO CANCELLATION FOR HANDS-FREE SPEAKERPHONES. Bryan S. Nollett and Douglas L. Jones
NONLINEAR ECHO CANCELLATION FOR HANDS-FREE SPEAKERPHONES Bryan S. Nollett and Douglas L. Jones Coordinated Science Laboratory University of Illinois at Urbana-Champaign 1308 W. Main St. Urbana, IL 61801
More informationPrediction-based adaptive control of a class of discrete-time nonlinear systems with nonlinear growth rate
www.scichina.com info.scichina.com www.springerlin.com Prediction-based adaptive control of a class of discrete-time nonlinear systems with nonlinear growth rate WEI Chen & CHEN ZongJi School of Automation
More informationParameter Bounds for Discrete-Time Hammerstein Models With Bounded Output Errors
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 48, NO. 10, OCTOBER 003 1855 Parameter Bounds for Discrete-Time Hammerstein Models With Bounded Output Errors V. Cerone and D. Regruto Abstract In this note,
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 informationSystem Identification: From Data to Model
1 : From Data to Model With Applications to Aircraft Modeling Division of Automatic Control Linköping University Sweden Dynamic Systems 2 A Dynamic system has an output response y that depends on (all)
More informationAccurate Model Identification for Non-Invertible MIMO Sandwich Block-Oriented Processes
Extended abstract for 006 AIChE Annual Meeting, Topic: Process Modeling and Identification, Session: 10B0 Accurate Model Identification for Non-Invertible MIMO Sandwich Block-Oriented Processes Swee-Teng
More information( ) ( = ) = ( ) ( ) ( )
( ) Vρ C st s T t 0 wc Ti s T s Q s (8) K T ( s) Q ( s) + Ti ( s) (0) τs+ τs+ V ρ K and τ wc w T (s)g (s)q (s) + G (s)t(s) i G and G are transfer functions and independent of the inputs, Q and T i. Note
More informationChapter 4 Neural Networks in System Identification
Chapter 4 Neural Networks in System Identification Gábor HORVÁTH Department of Measurement and Information Systems Budapest University of Technology and Economics Magyar tudósok körútja 2, 52 Budapest,
More informationSemiparametric Identification of Wiener Systems Using a Single Harmonic Input and Retrospective Cost Optimization
American Control Conference Marriott Waterfront, Baltimore, MD, USA June 3-July, ThC7. Semiparametric Identification of Wiener Systems Using a Single Harmonic Input and Retrospective Cost Optimization
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 informationHand Written Digit Recognition using Kalman Filter
International Journal of Electronics and Communication Engineering. ISSN 0974-2166 Volume 5, Number 4 (2012), pp. 425-434 International Research Publication House http://www.irphouse.com Hand Written Digit
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 informationRecursive Least Squares for an Entropy Regularized MSE Cost Function
Recursive Least Squares for an Entropy Regularized MSE Cost Function Deniz Erdogmus, Yadunandana N. Rao, Jose C. Principe Oscar Fontenla-Romero, Amparo Alonso-Betanzos Electrical Eng. Dept., University
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 informationNonlinear System Identification Using Maximum Likelihood Estimation
ISS (Print) : 30 3765 ISS (Online): 78 8875 (An ISO 397: 007 Certified Organizatio Vol., Special Issue, December 03 onlinear System Identification Using Maximum Lielihood Estimation Siny Paul, Bindu Elias
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 informationA STATE-SPACE APPROACH FOR THE ANALYSIS OF WAVE AND DIFFUSION FIELDS
ICASSP 2015 A STATE-SPACE APPROACH FOR THE ANALYSIS OF WAVE AND DIFFUSION FIELDS Stefano Maranò Donat Fäh Hans-Andrea Loeliger ETH Zurich, Swiss Seismological Service, 8092 Zürich ETH Zurich, Dept. Information
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 informationQ3. Derive the Wiener filter (non-causal) for a stationary process with given spectral characteristics
Q3. Derive the Wiener filter (non-causal) for a stationary process with given spectral characteristics McMaster University 1 Background x(t) z(t) WF h(t) x(t) n(t) Figure 1: Signal Flow Diagram Goal: filter
More information! # % & () +,.&/ 01),. &, / &
! # % & () +,.&/ ),. &, / & 2 A NEW METHOD FOR THE DESIGN OF ENERGY TRANSFER FILTERS Xiaofeng Wu, Z Q Lang and S. A. Billings Department of Automatic Control and Systems Engineering The University of Sheffield
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 informationLinearization problem. The simplest example
Linear Systems Lecture 3 1 problem Consider a non-linear time-invariant system of the form ( ẋ(t f x(t u(t y(t g ( x(t u(t (1 such that x R n u R m y R p and Slide 1 A: f(xu f(xu g(xu and g(xu exist and
More informationIMPROVEMENTS IN MODAL PARAMETER EXTRACTION THROUGH POST-PROCESSING FREQUENCY RESPONSE FUNCTION ESTIMATES
IMPROVEMENTS IN MODAL PARAMETER EXTRACTION THROUGH POST-PROCESSING FREQUENCY RESPONSE FUNCTION ESTIMATES Bere M. Gur Prof. Christopher Niezreci Prof. Peter Avitabile Structural Dynamics and Acoustic Systems
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 informationA Canonical Genetic Algorithm for Blind Inversion of Linear Channels
A Canonical Genetic Algorithm for Blind Inversion of Linear Channels Fernando Rojas, Jordi Solé-Casals, Enric Monte-Moreno 3, Carlos G. Puntonet and Alberto Prieto Computer Architecture and Technology
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 informationNON-LINEAR PARAMETER ESTIMATION USING VOLTERRA AND WIENER THEORIES
Journal of Sound and Vibration (1999) 221(5), 85 821 Article No. jsvi.1998.1984, available online at http://www.idealibrary.com on NON-LINEAR PARAMETER ESTIMATION USING VOLTERRA AND WIENER THEORIES Department
More informationDevelopment of Nonlinear Black Box Models using Orthonormal Basis Filters: A Review*
Development of Nonlinear Black Box Models using Orthonormal Basis Filters: A Review* Sachin C. atwardhan Abstract Over the last two decades, there has been a growing interest in the use of orthonormal
More informationA technique for simultaneous parameter identification and measurement calibration for overhead transmission lines
A technique for simultaneous parameter identification and measurement calibration for overhead transmission lines PAVEL HERING University of West Bohemia Department of cybernetics Univerzitni 8, 36 4 Pilsen
More informationNonlinear State Estimation Methods Overview and Application to PET Polymerization
epartment of Biochemical and Chemical Engineering Process ynamics Group () onlinear State Estimation Methods Overview and Polymerization Paul Appelhaus Ralf Gesthuisen Stefan Krämer Sebastian Engell The
More informationInput-output data sets for development and benchmarking in nonlinear identification
Input-output data sets for development and benchmarking in nonlinear identification Torbjörn Wigren Division of Systems and Control, Department of Information Technology, Uppsala University, PO box 337,
More informationIdentification of Nonlinear Dynamic Systems with Multiple Inputs and Single Output using discrete-time Volterra Type Equations
Identification of Nonlinear Dnamic Sstems with Multiple Inputs and Single Output using discrete-time Volterra Tpe Equations Thomas Treichl, Stefan Hofmann, Dierk Schröder Institute for Electrical Drive
More informationA SQUARE ROOT ALGORITHM FOR SET THEORETIC STATE ESTIMATION
A SQUARE ROO ALGORIHM FOR SE HEOREIC SAE ESIMAION U D Hanebec Institute of Automatic Control Engineering echnische Universität München 80290 München, Germany fax: +49-89-289-28340 e-mail: UweHanebec@ieeeorg
More informationECE 636: Systems identification
ECE 636: Systems identification Lectures 3 4 Random variables/signals (continued) Random/stochastic vectors Random signals and linear systems Random signals in the frequency domain υ ε x S z + y Experimental
More informationCONE CRUSHER MODEL IDENTIFICATION USING BLOCK-ORIENTED SYSTEMSWITH ORTHONORMAL BASIS FUNCTIONS
CONE CRUSHER MODEL IDENIFICAION USING BLOCK-ORIENED SYSEMSWIH ORHONORMAL BASIS FUNCIONS Oleksii Mykhailenko 1 1 Department of Power Supply and Energy Management, State institution of higher education «KryvyiRih
More informationCambridge University Press Nonparametric System Identification Wlodzimierz Greblicki and Miroslaw Pawlak Excerpt More information
1 Introduction System identification, as a particular process of statistical inference, exploits two types of information. The first is experiment; the other, called a priori, is known before making any
More informationLayered Polynomial Filter Structures
INFORMATICA, 2002, Vol. 13, No. 1, 23 36 23 2002 Institute of Mathematics and Informatics, Vilnius Layered Polynomial Filter Structures Kazys KAZLAUSKAS, Jaunius KAZLAUSKAS Institute of Mathematics and
More information2 Introduction of Discrete-Time Systems
2 Introduction of Discrete-Time Systems This chapter concerns an important subclass of discrete-time systems, which are the linear and time-invariant systems excited by Gaussian distributed stochastic
More informationWavelet-Based Nonparametric Modeling of Hierarchical Functions in Colon Carcinogenesis
Wavelet-Based Nonparametric Modeling of Hierarchical Functions in Colon Carcinogenesis Jeffrey S. Morris University of Texas, MD Anderson Cancer Center Joint wor with Marina Vannucci, Philip J. Brown,
More informationPARAMETER ESTIMATION AND ORDER SELECTION FOR LINEAR REGRESSION PROBLEMS. Yngve Selén and Erik G. Larsson
PARAMETER ESTIMATION AND ORDER SELECTION FOR LINEAR REGRESSION PROBLEMS Yngve Selén and Eri G Larsson Dept of Information Technology Uppsala University, PO Box 337 SE-71 Uppsala, Sweden email: yngveselen@ituuse
More informationA Kind of Frequency Subspace Identification Method with Time Delay and Its Application in Temperature Modeling of Ceramic Shuttle Kiln
American Journal of Computer Science and Technology 218; 1(4): 85 http://www.sciencepublishinggroup.com/j/ajcst doi: 1.1148/j.ajcst.21814.12 ISSN: 2412X (Print); ISSN: 24111 (Online) A Kind of Frequency
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 informationCALCULATION METHOD FOR NONLINEAR DYNAMIC LEAST-ABSOLUTE DEVIATIONS ESTIMATOR
J. Japan Statist. Soc. Vol. 3 No. 200 39 5 CALCULAION MEHOD FOR NONLINEAR DYNAMIC LEAS-ABSOLUE DEVIAIONS ESIMAOR Kohtaro Hitomi * and Masato Kagihara ** In a nonlinear dynamic model, the consistency and
More informationImpulsive Noise Filtering In Biomedical Signals With Application of New Myriad Filter
BIOSIGAL 21 Impulsive oise Filtering In Biomedical Signals With Application of ew Myriad Filter Tomasz Pander 1 1 Division of Biomedical Electronics, Institute of Electronics, Silesian University of Technology,
More informationAn Exact Solution to Support Vector Mixture
An Exact Solution to Support Vector Mixture Monjed Ezzeddinne, icolas Lefebvre, and Régis Lengellé Abstract This paper presents a new version of the SVM mixture algorithm initially proposed by Kwo for
More informationStochastic Collocation Methods for Polynomial Chaos: Analysis and Applications
Stochastic Collocation Methods for Polynomial Chaos: Analysis and Applications Dongbin Xiu Department of Mathematics, Purdue University Support: AFOSR FA955-8-1-353 (Computational Math) SF CAREER DMS-64535
More informationOnline Appendix. j=1. φ T (ω j ) vec (EI T (ω j ) f θ0 (ω j )). vec (EI T (ω) f θ0 (ω)) = O T β+1/2) = o(1), M 1. M T (s) exp ( isω)
Online Appendix Proof of Lemma A.. he proof uses similar arguments as in Dunsmuir 979), but allowing for weak identification and selecting a subset of frequencies using W ω). It consists of two steps.
More informationBayesian and Monte Carlo change-point detection
Bayesian and Monte Carlo change-point detection ROMN CMEJL, PVEL SOVK, MIROSLV SRUPL, JN UHLIR Department Circuit heory Czech echnical University in Prague echnica, 66 7 Prague 6 CZECH REPUBLIC bstract:
More informationNonparametric inference for ergodic, stationary time series.
G. Morvai, S. Yakowitz, and L. Györfi: Nonparametric inference for ergodic, stationary time series. Ann. Statist. 24 (1996), no. 1, 370 379. Abstract The setting is a stationary, ergodic time series. The
More informationIssues on quantile autoregression
Issues on quantile autoregression Jianqing Fan and Yingying Fan We congratulate Koenker and Xiao on their interesting and important contribution to the quantile autoregression (QAR). The paper provides
More informationAnalysis of Finite Wordlength Effects
Analysis of Finite Wordlength Effects Ideally, the system parameters along with the signal variables have infinite precision taing any value between and In practice, they can tae only discrete values within
More informationDiscretization of Continuous Linear Anisochronic Models
Discretization of Continuous Linear Anisochronic Models Milan Hofreiter Czech echnical University in Prague Faculty of Mechanical Engineering echnicka 4, 66 7 Prague 6 CZECH REPUBLIC E-mail: hofreite@fsid.cvut.cz
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 informationChapter 6: Nonparametric Time- and Frequency-Domain Methods. Problems presented by Uwe
System Identification written by L. Ljung, Prentice Hall PTR, 1999 Chapter 6: Nonparametric Time- and Frequency-Domain Methods Problems presented by Uwe System Identification Problems Chapter 6 p. 1/33
More informationRoot-MUSIC Time Delay Estimation Based on Propagator Method Bin Ba, Yun Long Wang, Na E Zheng & Han Ying Hu
International Conference on Automation, Mechanical Control and Computational Engineering (AMCCE 15) Root-MUSIC ime Delay Estimation Based on ropagator Method Bin Ba, Yun Long Wang, Na E Zheng & an Ying
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 information14 - Gaussian Stochastic Processes
14-1 Gaussian Stochastic Processes S. Lall, Stanford 211.2.24.1 14 - Gaussian Stochastic Processes Linear systems driven by IID noise Evolution of mean and covariance Example: mass-spring system Steady-state
More informationSubspace Identification for Nonlinear Systems That are Linear in Unmeasured States 1
Subspace Identification for Nonlinear Systems hat are Linear in nmeasured States 1 Seth L Lacy and Dennis S Bernstein 2 Abstract In this paper we apply subspace methods to the identification of a class
More informationIntroduction to Systems Analysis and Decision Making Prepared by: Jakub Tomczak
Introduction to Systems Analysis and Decision Making Prepared by: Jakub Tomczak 1 Introduction. Random variables During the course we are interested in reasoning about considered phenomenon. In other words,
More informationNon-stationary extremal eigenvalue approximations in iterative solutions of linear systems and estimators for relative error
on-stationary extremal eigenvalue approximations in iterative solutions of linear systems and estimators for relative error Divya Anand Subba and Murugesan Venatapathi* Supercomputer Education and Research
More informationTAKEHOME FINAL EXAM e iω e 2iω e iω e 2iω
ECO 513 Spring 2015 TAKEHOME FINAL EXAM (1) Suppose the univariate stochastic process y is ARMA(2,2) of the following form: y t = 1.6974y t 1.9604y t 2 + ε t 1.6628ε t 1 +.9216ε t 2, (1) where ε is i.i.d.
More informationA nonparametric test for seasonal unit roots
Robert M. Kunst robert.kunst@univie.ac.at University of Vienna and Institute for Advanced Studies Vienna To be presented in Innsbruck November 7, 2007 Abstract We consider a nonparametric test for the
More informationSingular Kalman Filtering: New Aspects Based on an Alternative System Theory
Singular Kalman Filtering: New Aspects Based on an Alternative System heory René Arnošt and Pavel Žampa Department of Cybernetics & Cybernetic Systems Research Center University of West Bohemia Univerzitní
More informationPower Amplifier Linearization Using Multi- Stage Digital Predistortion Based On Indirect Learning Architecture
Power Amplifier Linearization Using Multi- Stage Digital Predistortion Based On Indirect Learning Architecture Sreenath S 1, Bibin Jose 2, Dr. G Ramachandra Reddy 3 Student, SENSE, VIT University, Vellore,
More informationSet-Membership Identification of Wiener models with noninvertible nonlinearity
Set-Membership Identification of Wiener models with noninvertible nonlinearity V. Cerone, Dipartimento di Automatica e Informatica Politecnico di Torino (Italy) e-mail: vito.cerone@polito.it, diego.regruto@polito.it
More informationClosed-loop Identification of Hammerstein Systems Using Iterative Instrumental Variables
Proceedings of the 18th World Congress The International Federation of Automatic Control Closed-loop Identification of Hammerstein Systems Using Iterative Instrumental Variables Younghee Han and Raymond
More informationNON-LINEAR CONTROL OF OUTPUT PROBABILITY DENSITY FUNCTION FOR LINEAR ARMAX SYSTEMS
Control 4, University of Bath, UK, September 4 ID-83 NON-LINEAR CONTROL OF OUTPUT PROBABILITY DENSITY FUNCTION FOR LINEAR ARMAX SYSTEMS H. Yue, H. Wang Control Systems Centre, University of Manchester
More informationIdentification of Nonlinear Systems Using Neural Networks and Polynomial Models
A.Janczak Identification of Nonlinear Systems Using Neural Networks and Polynomial Models A Block-Oriented Approach With 79 Figures and 22 Tables Springer Contents Symbols and notation XI 1 Introduction
More informationSparse Least Mean Square Algorithm for Estimation of Truncated Volterra Kernels
Sparse Least Mean Square Algorithm for Estimation of Truncated Volterra Kernels Bijit Kumar Das 1, Mrityunjoy Chakraborty 2 Department of Electronics and Electrical Communication Engineering Indian Institute
More informationChapter 2 Wiener Filtering
Chapter 2 Wiener Filtering Abstract Before moving to the actual adaptive filtering problem, we need to solve the optimum linear filtering problem (particularly, in the mean-square-error sense). We start
More informationMultivariate ARMA Processes
LECTURE 8 Multivariate ARMA Processes A vector y(t) of n elements is said to follow an n-variate ARMA process of orders p and q if it satisfies the equation (1) A 0 y(t) + A 1 y(t 1) + + A p y(t p) = M
More informationRefined Instrumental Variable Methods for Identifying Hammerstein Models Operating in Closed Loop
Refined Instrumental Variable Methods for Identifying Hammerstein Models Operating in Closed Loop V. Laurain, M. Gilson, H. Garnier Abstract This article presents an instrumental variable method dedicated
More informationA. Poulimenos, M. Spiridonakos, and S. Fassois
PARAMETRIC TIME-DOMAIN METHODS FOR NON-STATIONARY RANDOM VIBRATION IDENTIFICATION AND ANALYSIS: AN OVERVIEW AND COMPARISON A. Poulimenos, M. Spiridonakos, and S. Fassois DEPARTMENT OF MECHANICAL & AERONAUTICAL
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 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 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 informationThe Rationale for Second Level Adaptation
The Rationale for Second Level Adaptation Kumpati S. Narendra, Yu Wang and Wei Chen Center for Systems Science, Yale University arxiv:1510.04989v1 [cs.sy] 16 Oct 2015 Abstract Recently, a new approach
More informationBlind channel deconvolution of real world signals using source separation techniques
Blind channel deconvolution of real world signals using source separation techniques Jordi Solé-Casals 1, Enric Monte-Moreno 2 1 Signal Processing Group, University of Vic, Sagrada Família 7, 08500, Vic
More informationParameter-dependent normal forms for bifurcation equations of dynamical systems
Parameter-dependent normal forms for bifurcation equations of dynamical systems Angelo uongo Dipartento di Ingegneria delle Strutture delle Acque e del Terreno Università dell Aquila Monteluco Roio 67040
More informationGradient-Adaptive Algorithms for Minimum Phase - All Pass Decomposition of an FIR System
1 Gradient-Adaptive Algorithms for Minimum Phase - All Pass Decomposition of an FIR System Mar F. Flanagan, Member, IEEE, Michael McLaughlin, and Anthony D. Fagan, Member, IEEE Abstract Adaptive algorithms
More informationNONLINEAR SYSTEM IDENTIFICATION VIA WAVELET EXPANSIONS. Z. Hasiewicz* M. Pawlak**
Copyright IFAC System Identification Santa Barbara, California, USA, 2000 NONLINEAR SYSTEM IDENTIFICATION VIA WAVELET EXPANSIONS Z. Hasiewicz* M. Pawlak** * Institute of Engineering Cybernetics, Wroclaw
More informationWavelet Methods for Time Series Analysis. Part IV: Wavelet-Based Decorrelation of Time Series
Wavelet Methods for Time Series Analysis Part IV: Wavelet-Based Decorrelation of Time Series DWT well-suited for decorrelating certain time series, including ones generated from a fractionally differenced
More informationarxiv: v1 [cs.sy] 22 Aug 2017
Parametric identification of parallel Wiener-Hammerstein systems Maarten Schoukens, Anna Marconato, Rik Pintelon, Gerd Vandersteen, Yves Rolain Vrije Universiteit Brussel (VUB), Dept. ELEC, Pleinlaan 2,
More informationIdentification of Unknown Functions in Dynamic Systems
Identification of Unknown Functions in Dynamic Systems Stéphane Binczak, Eric Busvelle, Jean-Paul Gauthier and Sabir Jacquir Abstract We consider the problem of representing a complex process by a simple
More informationMaximum Likelihood Diffusive Source Localization Based on Binary Observations
Maximum Lielihood Diffusive Source Localization Based on Binary Observations Yoav Levinboo and an F. Wong Wireless Information Networing Group, University of Florida Gainesville, Florida 32611-6130, USA
More informationSYSTEM RECONSTRUCTION FROM SELECTED HOS REGIONS. Haralambos Pozidis and Athina P. Petropulu. Drexel University, Philadelphia, PA 19104
SYSTEM RECOSTRUCTIO FROM SELECTED HOS REGIOS Haralambos Pozidis and Athina P. Petropulu Electrical and Computer Engineering Department Drexel University, Philadelphia, PA 94 Tel. (25) 895-2358 Fax. (25)
More informationADSP ADSP ADSP ADSP. Advanced Digital Signal Processing (18-792) Spring Fall Semester, Department of Electrical and Computer Engineering
Advanced Digital Signal rocessing (18-792) Spring Fall Semester, 201 2012 Department of Electrical and Computer Engineering ROBLEM SET 8 Issued: 10/26/18 Due: 11/2/18 Note: This problem set is due Friday,
More informationFeedback Control of Turbulent Wall Flows
Feedback Control of Turbulent Wall Flows Dipartimento di Ingegneria Aerospaziale Politecnico di Milano Outline Introduction Standard approach Wiener-Hopf approach Conclusions Drag reduction A model problem:
More informationState Space Representation of Gaussian Processes
State Space Representation of Gaussian Processes Simo Särkkä Department of Biomedical Engineering and Computational Science (BECS) Aalto University, Espoo, Finland June 12th, 2013 Simo Särkkä (Aalto University)
More informationStochastic Optimization with Inequality Constraints Using Simultaneous Perturbations and Penalty Functions
International Journal of Control Vol. 00, No. 00, January 2007, 1 10 Stochastic Optimization with Inequality Constraints Using Simultaneous Perturbations and Penalty Functions I-JENG WANG and JAMES C.
More informationAnalyticity and Derivatives of Entropy Rate for HMP s
Analyticity and Derivatives of Entropy Rate for HMP s Guangyue Han University of Hong Kong Brian Marcus University of British Columbia 1 Hidden Markov chains A hidden Markov chain Z = {Z i } is a stochastic
More information