ROBUST FREQUENCY DOMAIN ARMA MODELLING. Jonas Gillberg Fredrik Gustafsson Rik Pintelon
|
|
- Baldwin Blair
- 5 years ago
- Views:
Transcription
1 ROBUST FREQUENCY DOMAIN ARMA MODELLING Jonas Gillerg Fredrik Gustafsson Rik Pintelon Department of Electrical Engineering, Linköping University SE Linköping, Sweden Department ELEC, Vrije Universiteit Brussel Pleinlaan 2, 15 Brussel, Belgium Astract: In this paper a method for the rejection of frequency domain outliers is proposed. The algorithm is ased on the work y Huer on M-estimators and the concept of influence function introduced y Hampel. The estimation takes place in the context of frequency domain continuous-time ARMA modelling, ut the method can e also e applied to the discrete time case.copyright c 26SYSID Keywords: roust estimation, outlier rejection, continuous-time, ARMA, CARMA, frequency Domain, influence function, M-estimator 1. INTRODUCTION A characteristic industrial signal processing prolem is viration analysis using tachometer measurements on rotating axles. Motivated y recent advances in system identification in the frequency domain Pintelon and Schoukens, 21; Ljung, 1999), a frequency domain approach for such analysis was outlined in an earlier paper Gillerg and Gustafsson, 25) y the authors and compared at a theoretical level to the time domain algorithm proposed in Persson, 22; Persson and Gustafsson, 21). Among other things, the new method involved virational modelling using a continuous-time AR model and the restriction of the frequency interval to e used in the estimation. Interpolation was also performed in order to counter the effect of non-uniform sampling. In the same paper the main specifications on a procedure aimed at high-sensitivity viration analysis were listed: 1) Being ased on parametric physical models of the viration. 2) Operate on short data atches in a prespecified speed interval where the data pass several quality checks. 3) Potential to reject wide and disturances that are non-interfering with the viration. 4) Potential to reject narrow and disturances that are interfering with the viration. The methods given in Persson, 22; Persson and Gustafsson, 21) and Gillerg and Gustafsson, 25) successfully solves the first three specifications, ut not the last one. The prolem of narrow-and disturances occurs in several applications. For instance in an automotive drive-line, see Figure 1, where the virations can e due to engine knock or gear shifts. In rootics there are virations from the load, just to mention a few. In this paper the focus will therefore e on roustness against these narrow and disturances.
2 where the model parameter vector θ includes the noise variance λ whose true value is σ 2 ). The transfer function G can e parameterized y θ in an aritrary way, for example y the conventional numerator and denominator parameters: Gp, θ) = Bp) Ap) Ap) = p n + a 1 p n 1 + a 2 p n a n Bp) = p m + 1 p m m θ = [a 1 a 2... a n m λ] T. 4) f [Hz] Fig. 1. Smoothed spectrum for virations extracted from ABS system signals 2. OUTLINE The outline of the paper will e the following. First in Section 3 the frequency domain continuous-time ARMA CARMA) modelling and identification approach will e descried. Then in Section 4 the identification criterion is modified in order to facilitate rejection of frequency domain outliers. The concept of M-estimators and the influence function are also introduced. In Section 5 the asymptotic variance and gross sensitivity of the estimator are calculated, and measures are taken in order to assure consistency in the case of correct modelling. A certain choice of criterion is also proved to make an optimal tradeoff etween decreased ias due to outliers and increased variance due to a more insensitive method. Finally in Section 6 the algorithm is illustrated y a numerical example. 3. FREQUENCY DOMAIN CARMA ESTIMATION In this paper we shall consider continuous-time ARMA models represented as yt) = G c p)et) 1) where et) is continuous time white noise such that E[et)] = E[et)es)] = σ 2 δt s) The operator p is here the differentiation operator. We assume that Gp) is strictly proper, so yt) itself does not have a white-noise component, ut is a well defined second order, stationary process. Its spectrum spectral density) can e written as Φ c ω) = σ 2 G c iω) 2 2) We shall consider a general model parameterization G c p, θ) 3) Let us define the truncated Fourier transformation of the continuous time output {yt) : t [, T ]} in expression 1) as Y T iω) = 1 T T yt)e iωt dt. A complicating element in a practical estimation procedure is that we do not have access to the entire continuous time realization of the output. Instead we have a finite numer of samples of yt) at time instances {t 1, t 2,..., t N }. This is an important issue which have een dealt with in the case of equidistant sampling in a previous paper y the first author Gillerg and Ljung, 25). The case of non-uniform sampling we refer the reader to a technical reports associated with two other papers sumitted to this conference Gillerg and Ljung, 25a) Gillerg and Ljung, 24a). However, if we have an estimate of the continuoustime Fourier transform, the continuous-time periodogram is ˆΦiω n ) = Y T iω n ) 2 Exp Φω n, θ ) 5) and the transform will e asymptotically independent at the frequencies ω k, k = 1,..., N ω where ω k = 2πk/T, k N. In this paper we will focus treat the continuous-time periodogram estimates as our measurements and will consider them to e completely uncorrelated. When an estimate of the power spectrum is availale a model can e identified y the following Maximum-Likelihood ML) procedure descried in Gillerg and Ljung, 24) and Gillerg, 24) ˆθ = arg min θ k=1 ˆΦiω k ) Φiω k, θ) + log Φiω k, θ). 6) where the asymptotic Fisher information at the true parameters θ is Jθ ) = Φ θ ω n, θ ) Φω n, θ ) Φ θ ω n, θ ) T Φω n, θ ). 7) This approach forms the asis for the work presented in the remaining part of the paper.
3 The method in 6) can also e descried as the θ which is the solution to the vector equation ) Φ θ ω n, θ) ˆΦωn ) Φω n, θ) Φω n, θ) 1 =. where Φ θ ω n, θ) is the gradient of the spectrum with respect to the parameters. 4. M-ESTIMATORS AND OUTLIER REJECTION If the model is the correct one, the maximumlikelihood method is known to e the est possile optimal) estimator with consistency and asymptotic efficiency Wald, 1949)Cramér, 146). If there, on the other hand, are unmodelled outliers present, the method might not e optimal. This may lead us to consider roust estimators which give up efficiency at the true model in exchange for reasonale performance if the model is not the truecasella and Berger, 22). Therefore we introduce the following method ΨˆΦ, θ) = Φ θ ω n, θ) ˆΦωn ) Φω n, θ) ) =. 8) inspired y the M-estimators introduced y Huer Huer, 1981). The unknown here is the function ψ which will e found using the influence function approach introduced y Hampel Hampel et al., 1986). Certain measures must however e taken in order to assure that the estimates are consistent when periodogram originate from the assumed model in 5). This implies that ψ ˆΦiωn ) Φiω n, θ ) ) ) df n ˆΦiωn ) = 9) ψ x) dg x) = 1) where we from now on define the distriutions F n Exp Φiω n, θ ) 11) G Exp 1 12) In this context we define the solution to equation 8) as a vector valued stochastic variale T F ) with a distriution dependent on F = F 1,..., F Nω ), the distriution of the periodogram ˆΦiω n ). The so called influence function introduced y Hampel is then defined as IF ˆΦ; T, F ) = lim h T 1 h)f + hδˆφ) T F ) h and measures the asymptotic ias caused y contamination in the data. Here is the multivariale measure which puts the mass 1 at ˆΦ δˆφ = ˆΦiω1 ),..., ˆΦωNω ). 5. INFLUENCE FUNCTION AND ASYMPTOTIC VARIANCE In the case of the estimator in 8) it is possile to show that see p.11 in Hampel et al., 1986)) IF ˆΦ; T, F ) = MΨ, T F )) 1 ΨˆΦ, T F )) where ) MΨ, F ) = ΨˆΦ, θ) θ=t F ) df ˆΦ). It is also possile to show see p.85 in Hampel et al., 1986)) that the asymptotic variance of the parameter estimates will e V T, F ) = IF ˆΦ; T, F )IF ˆΦ; T, F )df ˆΦ) =Mψ, F ) 1 Qψ, F )Mψ, F ) T where QΨ, F ) = ΨˆΦ, T F ))Ψ T ˆΦ, T F ))df ˆΦ) Here = and ) MΨ, F ) = ΨˆΦ, θ) Φ θ ω n, θ) ˆΦωn ) Φω n, θ) Φ θ ω n, θ ) Φ θ ω n, θ ) T Φω n, θ ) Φω n, θ ) + Φ θ ω n, θ ) Φ θ ω n, θ ) T Φω n, θ ) Φω n, θ ) df ˆΦ) θ=θ )) θ=θ df ˆΦ) Φ )) θ ω n, θ) ˆΦωn ) = Φω n, θ) θ=θ Φ θ ω ) ) n, θ ) ˆΦωn ) ψ Φω n, θ ) Φω n, θ ) ) ˆΦω n ) ˆΦωn ) Φω n, θ ) ψ. Φω n, θ ) Because of 9) and 1) the expected value of first expression on the right hand side of the equation aove will e zero. Therefore we will have MΨ, F ) = Jθ) xψ x)dgx). Since the periodogram at different frequencies are assumed independent we also have QΨ, F ) = ΨˆΦ, T F ))Ψ T ˆΦ, T F ))df ˆΦ) =Jθ) ψ 2 x)dgx) where J is the Fischer information in 7) and the asymptotic variance will e ψ 2 x)dgx) V T, F ) = xψ x)dgx) ) 2 J 1 θ)
4 5.1 Optimal B-Roust Estimator Assume now that we want to find a function ψ that limits the gross sensitivity measured as γut, F ) = sup{ IF ˆΦ; T, F ) } 13) ˆΦ where. means taking the asolute value of each vector component, while at the same time minimizes the trace of the variance T r V T, F ). In our case we will have ψ 2 x)dgx) T r V T, F ) =T r xψ x)dgx) ) 2 J 1 θ) and ψ 2 x)dgx) = xψ x)dgx) ) 2 T r J 1 θ) IF ˆΦ; T, F ) = MΨ, T G)) 1 ΨˆΦ, T F )) ) Jθ) 1 Φ = θ ω n, θ) xψ x)dgx) Φω n, θ) ψ ˆΦωn ) Φω n, θ) xψ x)dgx) Jθ) 1 Φ θ ω n, θ) Φω n, θ) 5.2 Computing a The parameter in 14) can e considered a tuning factor which is selected y the user in order to strike a alance etween the asymptotic ias and variance of the parameter estimates. The factor a is then selected such that the estimator is asymptotically uniased when the data is produced y the model. That is ψx)dgx) = [x 1 a] e x dx =. This means that three different cases have to e considered. 5.3 Case 1 The first, when 1 a which is trivial since yields =. [x 1 a] e x dx = = 16) According to Theorem 1 and Equation 2.4.1) on p. 122 in Hampel et al., 1986) a version of the so called Huer function ψx) =[x 1 a] 14) if < x 1 a = x 1 a if < x 1 a if x 1 a minimizes ψ 2 x)dgx) xψ x)dgx) ) 2 among all mappings ψ that satisfy ψx)dgx) = sup ψx) x xψ x)dgx) xψ x)dgx) Therefore, using the function in 14) in ) Φ θ ω n, θ) ˆΦωn ) = 15) Φω n, θ) will minimize while T r V T, F ) = sup ˆΦ ψ 2 x)dgx) xψ x)dgx) ) 2 T r J 1 θ) Nω Jθ) 1 Φ θ ω n,θ) Φω n,θ) IF ˆΦ; T, F ) xψ x)dgx) 5.4 Case 2 The second case is when < 1 a. Here [x 1 a] e x dx = 1+a+ + 1+a+ x 1 a)e x dx The first integral in this expression will e 1+a+ x 1 a)e x dx = [ e x xe x a) ] 1+a+ = [ a x)e x] 1+a+ = 1 )e e 1+a) a The second integral is 1+a+ e x dx = [ e x] 1+a+ = e e 1+a) x 1 a)e x dx. This means that given we have to chose a) such that 1 )e e 1+a) a + e e 1+a) = This means that in this case { e e 1+a) + a = 1 + a e e 1+a) a =
5 5.5 Case 3 In the third case we have 1 a and 1+a+ 1+a [x 1 a] e x dx = x 1 a)e x dx + 1+a 1+a+ The first integral will in this case e 1+a while the second is 1+a+ 1+a x 1 a)e x dx = e x dx e x dx = [ e x] 1+a = e e 1+a) x 1 a)e x dx. = [ a x)e x] 1+a+ 1+a = a 1 a )e e 1+a) a 1 a + )e e 1+a) = 1)e 1)e )e 1+a) The third integral is the same as previously 1+a+ This means that we have x 1 a)e x = e e 1+a) [x 1 a] e x dx = e e 1+a) + 1)e 1)e )e 1+a) + e e 1+a) = e e )e 1+a) = and the solution is a = ln e e 1 + a 1 a Fig. 2. The parameter a as a function of the parameter. The reakpoint etween the two solution occurs at.7968 The ojective is to estimate the parameters of the continuous-time spectrum Φiω) = λ iω) 2 + a 1 iω) + a ) where a 1 = 3 and a 2 = 2. The underlying continuous-time system is the autoregressive model σ yt) = p 2 et) 18) + a 1 p + a 2 where et) is continuous-time Gaussian white noise. In Figure 3 we have computed the standard deviation for the parameter estimates of the continuous-time model found in 18). Both the new and the original ML method have een used. Simulations have een performed using N MC = 1 Monte-Carlo runs for each value of. As expected the standard deviation of the roust method will increase as decrease, ut it is interesting to note that it only approximately doule for such a small value as =.1. In the estimation we have used frequencies w = { :.1 : 2π} 5.6 Equations σ a1.4.2 This means that a) is defined y the equations e e 1+a) + a = if 1 + a e 1+a = ln e e if 1 + a In Figure 2 the parameter a is presented as a function of the parameter. The reakpoint etween the two sets of solutions occurs when NUMERICAL EXAMPLE In this section we will compare the efficiency and roustness of the proposed estimation method in 8) and 14) with that of the ML approach in 6). σ a2 σ λ Fig. 3. The standard deviations of the parameter estimates. Roust dotted) and ML method dashed). In Figure 4 we have estimated the same parameters as in the example prolem aove. In order to simulate outliers, we have also introduced an
6 additive, random and exponentially distriuted disturance at 5 % of all frequencies, with a variance which is twice the maximum value of the power spectrum in 17). From the figure one can a 1 a 2 λ Fig. 4. The ias for the parameter estimates. True solid), ML method dashed) and roust dotted). see that the outliers cause ias in the parameter estimates when the ML method is utilized. This ias is then reduced y approximately 5 % with the use of the more roust criterion. 7. CONCLUSIONS In this paper a method for the rejection of frequency domain outliers is proposed. The algorithm is ased on the work y Huer on M- estimators and the concept of influence function introduced y Hampel. The estimation takes place in the context of frequency domain continuoustime ARMA modelling, ut the method can e also e applied to the discrete time case. It is also proved that a certain choice of criterion will produce an optimal tradeoff etween ias and variance under certain assumptions. Finally the method is illustrated y a numerical example which shows that the ias can e reduced y approximately 5%. REFERENCES Casella, G. and R.L. Berger 22). Statistical Inference. Duxury. Pacific Grove, CA. Cramér, H. 146). Mathematical Methods of Statistics. Princeton University Press. Princeton, NJ. Gillerg, J. 24). Methods for frequency domain estimation of continuous-time models. Lic. Dissertation Linköping Studies in Science and Technology Thesis No Department of Electrical Engineering, Linköping University. SE Linköping, Sweden. Gillerg, J. and F. Gustafsson 25). Frequencydomain continuous-time AR modelling using non-uniformly sampled measurements. In: Proceedings of IEEE Int Conf on Acoustics, Speech and Signal Processing in Philadelphia, USA. Gillerg, J. and L. Ljung 24a). Frequencydomain ARMA models: Interpolation and nonuniform sampling. Technical Report LiTH- ISY-R Department of Electrical Engineering, Linköping University. SE Linköping, Sweden. Gillerg, J. and L. Ljung 24). Frequencydomain ARMA models: Interpolation and nonuniform sampling. Technical Report LiTH- ISY-R Department of Electrical Engineering, Linköping University. SE Linköping, Sweden. Gillerg, J. and L. Ljung 25a). Frequencydomain ARMA models from non-uniformly sampled data. Technical Report LiTH-ISY-R Department of Electrical Engineering, Linköping University. SE Linköping, Sweden. Gillerg, J. and L. Ljung 25). Frequencydomain ARMA models from sampled data. In: Proc. 16th IFAC World Congress. Hampel, F. R., E.M. Ronchetti, P.J. Rousseeuw and W.A. Stael 1986). Roust Statistics: the approach ased on influence function. John Wiley & Sons. Huer, P. 1981). Roust Statistics. John Wiley & Sons. Ljung, L. 1999). System Identification: Theory for the User. Prentice-Hall. Upper Saddle River. Persson, N. 22). Event ased sampling with applicaion to spectral estimation. Lic. Dissertation Linköping Studies in Science and Technology Thesis No Department of Electrical Engineering, Linköping University. SE Linköping, Sweden. Persson, N. and F. Gustafsson 21). Event ased sampling with application to viration analysis in pneumatic tires. In: Proceedings of IEEE Int Conf on Acoustics, Speech and Signal Processing in Salt Lake City, USA. Pintelon, R. and J. Schoukens 21). System Identification - A Frequency Domain Approach. IEEE Press. Piscataway, NJ. Wald, A. 1949). Note on the consistency of the maximum likelihood estimate. Anti. Math. Statis. 2,
On 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 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 informationIN this paper, we consider the estimation of the frequency
Iterative Frequency Estimation y Interpolation on Fourier Coefficients Elias Aoutanios, MIEEE, Bernard Mulgrew, MIEEE Astract The estimation of the frequency of a complex exponential is a prolem that is
More informationNon-parametric identification
Non-parametric Non-parametric Transient Step-response using Spectral Transient Correlation Frequency function estimate Spectral System Identification, SSY230 Non-parametric 1 Non-parametric Transient Step-response
More informationMODELING AND IDENTIFICATION OF A MECHANICAL INDUSTRIAL MANIPULATOR 1
Copyright 22 IFAC 15th Triennial World Congress, Barcelona, Spain MODELING AND IDENTIFICATION OF A MECHANICAL INDUSTRIAL MANIPULATOR 1 M. Norrlöf F. Tjärnström M. Östring M. Aberger Department of Electrical
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 informationParameter Estimation in a Moving Horizon Perspective
Parameter Estimation in a Moving Horizon Perspective State and Parameter Estimation in Dynamical Systems Reglerteknik, ISY, Linköpings Universitet State and Parameter Estimation in Dynamical Systems OUTLINE
More informationThe Local Polynomial Method for nonparametric system identification: improvements and experimentation
The Local Polynomial Method for nonparametric system identification: improvements and experimentation Michel Gevers, Rik Pintelon and Johan Schoukens Abstract The Local Polynomial Method (LPM) is a recently
More informationOn Identification of Cascade Systems 1
On Identification of Cascade Systems 1 Bo Wahlberg Håkan Hjalmarsson Jonas Mårtensson Automatic Control and ACCESS, School of Electrical Engineering, KTH, SE-100 44 Stockholm, Sweden. (bo.wahlberg@ee.kth.se
More informationFrequency-Domain Identification of Continuous-Time ARMA Models from Sampled Data
Technical report from Automatic Control at Linköpings universitet Frequency-Domain Identification of Continuous-Time ARMA Models from Sampled Data Jonas Gillberg, Lennart Ljung Division of Automatic Control
More informationGTLS ALGORITHMS IN THE FREQUENCY DOMAIN SYSTEM IDENTIFICATION USING NOISE INFORMATION OUT OF A FINITE NUMBER OF REPEATED INDEPENDENT REALIZATIONS
GTLS ALGORITHMS IN THE REQUENCY DOMAIN SYSTEM IDENTIICATION USING NOISE INORMATION OUT O A INITE NUMBER O REPEATED INDEPENDENT REALIZATIONS Gerd Versteen *, Rik Pintelon Johan Schoukens Vrije Universiteit
More informationform the corresponding predictions (3) it would be tempting to use ^y(tj) =G(q )u(t)+(ih (q ))(y(t)g(q )u(t)) (6) and the associated prediction errors
Some Results on Identifying Linear Systems Using Frequency Domain Data Lennart Ljung Department of Electrical Engineering Linkoping University S-58 83 Linkoping, Sweden Abstract The usefulness of frequency
More informationFinQuiz Notes
Reading 9 A time series is any series of data that varies over time e.g. the quarterly sales for a company during the past five years or daily returns of a security. When assumptions of the regression
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 informationSpectrum Opportunity Detection with Weak and Correlated Signals
Specum Opportunity Detection with Weak and Correlated Signals Yao Xie Department of Elecical and Computer Engineering Duke University orth Carolina 775 Email: yaoxie@dukeedu David Siegmund Department of
More informationSELECTIVE ANGLE MEASUREMENTS FOR A 3D-AOA INSTRUMENTAL VARIABLE TMA ALGORITHM
SELECTIVE ANGLE MEASUREMENTS FOR A 3D-AOA INSTRUMENTAL VARIABLE TMA ALGORITHM Kutluyıl Doğançay Reza Arablouei School of Engineering, University of South Australia, Mawson Lakes, SA 595, Australia ABSTRACT
More informationAccuracy Analysis of Time-domain Maximum Likelihood Method and Sample Maximum Likelihood Method for Errors-in-Variables Identification
Proceedings of the 17th World Congress The International Federation of Automatic Control Seoul, Korea, July 6-11, 8 Accuracy Analysis of Time-domain Maximum Likelihood Method and Sample Maximum Likelihood
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 informationFUNDAMENTAL FILTERING LIMITATIONS IN LINEAR NON-GAUSSIAN SYSTEMS
FUNDAMENTAL FILTERING LIMITATIONS IN LINEAR NON-GAUSSIAN SYSTEMS Gustaf Hendeby Fredrik Gustafsson Division of Automatic Control Department of Electrical Engineering, Linköpings universitet, SE-58 83 Linköping,
More informationLarge Sample Properties of Estimators in the Classical Linear Regression Model
Large Sample Properties of Estimators in the Classical Linear Regression Model 7 October 004 A. Statement of the classical linear regression model The classical linear regression model can be written in
More informationTheory of Maximum Likelihood Estimation. Konstantin Kashin
Gov 2001 Section 5: Theory of Maximum Likelihood Estimation Konstantin Kashin February 28, 2013 Outline Introduction Likelihood Examples of MLE Variance of MLE Asymptotic Properties What is Statistical
More informationAUTOMATIC CONTROL COMMUNICATION SYSTEMS LINKÖPINGS UNIVERSITET. Questions AUTOMATIC CONTROL COMMUNICATION SYSTEMS LINKÖPINGS UNIVERSITET
The Problem Identification of Linear and onlinear Dynamical Systems Theme : Curve Fitting Division of Automatic Control Linköping University Sweden Data from Gripen Questions How do the control surface
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 informationarxiv: v1 [stat.ml] 10 Apr 2017
Integral Transforms from Finite Data: An Application of Gaussian Process Regression to Fourier Analysis arxiv:1704.02828v1 [stat.ml] 10 Apr 2017 Luca Amrogioni * and Eric Maris Radoud University, Donders
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 informationOn the estimation of the K parameter for the Rice fading distribution
On the estimation of the K parameter for the Rice fading distribution Ali Abdi, Student Member, IEEE, Cihan Tepedelenlioglu, Student Member, IEEE, Mostafa Kaveh, Fellow, IEEE, and Georgios Giannakis, Fellow,
More informationAn algorithm for robust fitting of autoregressive models Dimitris N. Politis
An algorithm for robust fitting of autoregressive models Dimitris N. Politis Abstract: An algorithm for robust fitting of AR models is given, based on a linear regression idea. The new method appears to
More informationMaximum Likelihood Modeling Of Orbits Of Nonlinear ODEs
Maimum Likelihood Modeling Of Orbits Of onlinear ODEs orsten Söderström, orbjörn Wigren and Emad Abd-Elrady Department of Systems and Control, Information echnology, Uppsala Univ., PO Bo 337, SE-75 05
More informationFRF parameter identification with arbitrary input sequence from noisy input output measurements
21st International Symposium on Mathematical Theory of Networks and Systems July 7-11, 214. FRF parameter identification with arbitrary input sequence from noisy input output measurements mberto Soverini
More informationCOMPARISON OF THE ESTIMATORS OF THE LOCATION AND SCALE PARAMETERS UNDER THE MIXTURE AND OUTLIER MODELS VIA SIMULATION
(REFEREED RESEARCH) COMPARISON OF THE ESTIMATORS OF THE LOCATION AND SCALE PARAMETERS UNDER THE MIXTURE AND OUTLIER MODELS VIA SIMULATION Hakan S. Sazak 1, *, Hülya Yılmaz 2 1 Ege University, Department
More informationModel Specification Testing in Nonparametric and Semiparametric Time Series Econometrics. Jiti Gao
Model Specification Testing in Nonparametric and Semiparametric Time Series Econometrics Jiti Gao Department of Statistics School of Mathematics and Statistics The University of Western Australia Crawley
More informationSliding Window Recursive Quadratic Optimization with Variable Regularization
11 American Control Conference on O'Farrell Street, San Francisco, CA, USA June 29 - July 1, 11 Sliding Window Recursive Quadratic Optimization with Variable Regularization Jesse B. Hoagg, Asad A. Ali,
More informationA NEW INFORMATION THEORETIC APPROACH TO ORDER ESTIMATION PROBLEM. Massachusetts Institute of Technology, Cambridge, MA 02139, U.S.A.
A EW IFORMATIO THEORETIC APPROACH TO ORDER ESTIMATIO PROBLEM Soosan Beheshti Munther A. Dahleh Massachusetts Institute of Technology, Cambridge, MA 0239, U.S.A. Abstract: We introduce a new method of model
More informationREGLERTEKNIK AUTOMATIC CONTROL LINKÖPING
Time-domain Identication of Dynamic Errors-in-variables Systems Using Periodic Excitation Signals Urban Forssell, Fredrik Gustafsson, Tomas McKelvey Department of Electrical Engineering Linkping University,
More informationSystem Identification Approach Applied to Drift Estimation.
Instrumentation and Measurement Technology Conference - IMTC 007 Warsaw, Poland, May 1-3, 007 System Identification Approach Applied to Drift Estimation. Frans Verbeyst 1,, Rik Pintelon 1, Yves Rolain
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 informationTime domain identification, frequency domain identification. Equivalencies! Differences?
Time domain identification, frequency domain identification. Equivalencies! Differences? J. Schoukens, R. Pintelon, and Y. Rolain Vrije Universiteit Brussel, Department ELEC, Pleinlaan, B5 Brussels, Belgium
More informationAn Adaptive Clarke Transform Based Estimator for the Frequency of Balanced and Unbalanced Three-Phase Power Systems
07 5th European Signal Processing Conference EUSIPCO) An Adaptive Clarke Transform Based Estimator for the Frequency of Balanced and Unalanced Three-Phase Power Systems Elias Aoutanios School of Electrical
More informationPerformance of an Adaptive Algorithm for Sinusoidal Disturbance Rejection in High Noise
Performance of an Adaptive Algorithm for Sinusoidal Disturbance Rejection in High Noise MarcBodson Department of Electrical Engineering University of Utah Salt Lake City, UT 842, U.S.A. (8) 58 859 bodson@ee.utah.edu
More informationEXACT QUANTIFICATION OF VARIANCE ERROR
Copyright 00 IFAC 5th Triennial World Congress, Barcelona, Spain EXACT QUATIFICATIO OF VARIACE ERROR Brett inness, Håkan Hjalmarsson Dept. of Elec. & Comp. Eng, Uni. ewcastle, Australia. email:brett@ee.newcastle.edu.au,
More informationEstimation of Dynamic Regression Models
University of Pavia 2007 Estimation of Dynamic Regression Models Eduardo Rossi University of Pavia Factorization of the density DGP: D t (x t χ t 1, d t ; Ψ) x t represent all the variables in the economy.
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 informationCHANGE DETECTION WITH UNKNOWN POST-CHANGE PARAMETER USING KIEFER-WOLFOWITZ METHOD
CHANGE DETECTION WITH UNKNOWN POST-CHANGE PARAMETER USING KIEFER-WOLFOWITZ METHOD Vijay Singamasetty, Navneeth Nair, Srikrishna Bhashyam and Arun Pachai Kannu Department of Electrical Engineering Indian
More informationFree Probability, Sample Covariance Matrices and Stochastic Eigen-Inference
Free Probability, Sample Covariance Matrices and Stochastic Eigen-Inference Alan Edelman Department of Mathematics, Computer Science and AI Laboratories. E-mail: edelman@math.mit.edu N. Raj Rao Deparment
More informationA Robust Strategy for Joint Data Reconciliation and Parameter Estimation
A Robust Strategy for Joint Data Reconciliation and Parameter Estimation Yen Yen Joe 1) 3), David Wang ), Chi Bun Ching 3), Arthur Tay 1), Weng Khuen Ho 1) and Jose Romagnoli ) * 1) Dept. of Electrical
More informationEE531 (Semester II, 2010) 6. Spectral analysis. power spectral density. periodogram analysis. window functions 6-1
6. Spectral analysis EE531 (Semester II, 2010) power spectral density periodogram analysis window functions 6-1 Wiener-Khinchin theorem: Power Spectral density if a process is wide-sense stationary, the
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 informationUTILIZING PRIOR KNOWLEDGE IN ROBUST OPTIMAL EXPERIMENT DESIGN. EE & CS, The University of Newcastle, Australia EE, Technion, Israel.
UTILIZING PRIOR KNOWLEDGE IN ROBUST OPTIMAL EXPERIMENT DESIGN Graham C. Goodwin James S. Welsh Arie Feuer Milan Depich EE & CS, The University of Newcastle, Australia 38. EE, Technion, Israel. Abstract:
More informationAsymptotic Nonequivalence of Nonparametric Experiments When the Smoothness Index is ½
University of Pennsylvania ScholarlyCommons Statistics Papers Wharton Faculty Research 1998 Asymptotic Nonequivalence of Nonparametric Experiments When the Smoothness Index is ½ Lawrence D. Brown University
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 informationMethods for Frequency Domain Estimation of Continuous-Time Models
Linköping Studies in Science and Technology Thesis No. 1133 Methods for Frequency Domain Estimation of Continuous-Time Models Jonas Gillberg REGLERTEKNIK AUTOMATIC CONTROL LINKÖPING Division of Automatic
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 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 informationECE531 Lecture 10b: Maximum Likelihood Estimation
ECE531 Lecture 10b: Maximum Likelihood Estimation D. Richard Brown III Worcester Polytechnic Institute 05-Apr-2011 Worcester Polytechnic Institute D. Richard Brown III 05-Apr-2011 1 / 23 Introduction So
More informationComputer Intensive Methods in Mathematical Statistics
Computer Intensive Methods in Mathematical Statistics Department of mathematics johawes@kth.se Lecture 16 Advanced topics in computational statistics 18 May 2017 Computer Intensive Methods (1) Plan of
More informationON THE COMPARISON OF BOUNDARY AND INTERIOR SUPPORT POINTS OF A RESPONSE SURFACE UNDER OPTIMALITY CRITERIA. Cross River State, Nigeria
ON THE COMPARISON OF BOUNDARY AND INTERIOR SUPPORT POINTS OF A RESPONSE SURFACE UNDER OPTIMALITY CRITERIA Thomas Adidaume Uge and Stephen Seastian Akpan, Department Of Mathematics/Statistics And Computer
More informationRandom Eigenvalue Problems Revisited
Random Eigenvalue Problems Revisited S Adhikari Department of Aerospace Engineering, University of Bristol, Bristol, U.K. Email: S.Adhikari@bristol.ac.uk URL: http://www.aer.bris.ac.uk/contact/academic/adhikari/home.html
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 informationChapter 4. Theory of Tests. 4.1 Introduction
Chapter 4 Theory of Tests 4.1 Introduction Parametric model: (X, B X, P θ ), P θ P = {P θ θ Θ} where Θ = H 0 +H 1 X = K +A : K: critical region = rejection region / A: acceptance region A decision rule
More informationA Few Notes on Fisher Information (WIP)
A Few Notes on Fisher Information (WIP) David Meyer dmm@{-4-5.net,uoregon.edu} Last update: April 30, 208 Definitions There are so many interesting things about Fisher Information and its theoretical properties
More informationFrequency estimation by DFT interpolation: A comparison of methods
Frequency estimation by DFT interpolation: A comparison of methods Bernd Bischl, Uwe Ligges, Claus Weihs March 5, 009 Abstract This article comments on a frequency estimator which was proposed by [6] and
More informationTIGHT BOUNDS FOR THE FIRST ORDER MARCUM Q-FUNCTION
TIGHT BOUNDS FOR THE FIRST ORDER MARCUM Q-FUNCTION Jiangping Wang and Dapeng Wu Department of Electrical and Computer Engineering University of Florida, Gainesville, FL 3611 Correspondence author: Prof.
More informationNew Developments in Tail-Equivalent Linearization method for Nonlinear Stochastic Dynamics
New Developments in Tail-Equivalent Linearization method for Nonlinear Stochastic Dynamics Armen Der Kiureghian President, American University of Armenia Taisei Professor of Civil Engineering Emeritus
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 informationFAST AND ACCURATE DIRECTION-OF-ARRIVAL ESTIMATION FOR A SINGLE SOURCE
Progress In Electromagnetics Research C, Vol. 6, 13 20, 2009 FAST AND ACCURATE DIRECTION-OF-ARRIVAL ESTIMATION FOR A SINGLE SOURCE Y. Wu School of Computer Science and Engineering Wuhan Institute of Technology
More informationERASMUS UNIVERSITY ROTTERDAM Information concerning the Entrance examination Mathematics level 2 for International Business Administration (IBA)
ERASMUS UNIVERSITY ROTTERDAM Information concerning the Entrance examination Mathematics level 2 for International Business Administration (IBA) General information Availale time: 2.5 hours (150 minutes).
More informationIntroduction to Estimation Methods for Time Series models Lecture 2
Introduction to Estimation Methods for Time Series models Lecture 2 Fulvio Corsi SNS Pisa Fulvio Corsi Introduction to Estimation () Methods for Time Series models Lecture 2 SNS Pisa 1 / 21 Estimators:
More informationStatistical Methods for Handling Incomplete Data Chapter 2: Likelihood-based approach
Statistical Methods for Handling Incomplete Data Chapter 2: Likelihood-based approach Jae-Kwang Kim Department of Statistics, Iowa State University Outline 1 Introduction 2 Observed likelihood 3 Mean Score
More informationEstimation, Detection, and Identification
Estimation, Detection, and Identification Graduate Course on the CMU/Portugal ECE PhD Program Spring 2008/2009 Chapter 5 Best Linear Unbiased Estimators Instructor: Prof. Paulo Jorge Oliveira pjcro @ isr.ist.utl.pt
More informationCONTROL SYSTEMS, ROBOTICS, AND AUTOMATION - Vol. V - Frequency Domain System Identification - J. Schoukens and R. Pintelon
COTROL SYSTEMS, ROBOTICS, AD AUTOMATIO - Vol. V - Frequency Domain System Identification - J. Schoukens and FREQUECY DOMAI SYSTEM IDETIFICATIO J. Schoukens and Department ELEC, Vrije Universiteit Brussel,
More informationPatterns of Scalable Bayesian Inference Background (Session 1)
Patterns of Scalable Bayesian Inference Background (Session 1) Jerónimo Arenas-García Universidad Carlos III de Madrid jeronimo.arenas@gmail.com June 14, 2017 1 / 15 Motivation. Bayesian Learning principles
More informationOutline. A Central Limit Theorem for Truncating Stochastic Algorithms
Outline A Central Limit Theorem for Truncating Stochastic Algorithms Jérôme Lelong http://cermics.enpc.fr/ lelong Tuesday September 5, 6 1 3 4 Jérôme Lelong (CERMICS) Tuesday September 5, 6 1 / 3 Jérôme
More informationSensitivity and Asymptotic Error Theory
Sensitivity and Asymptotic Error Theory H.T. Banks and Marie Davidian MA-ST 810 Fall, 2009 North Carolina State University Raleigh, NC 27695 Center for Quantitative Sciences in Biomedicine North Carolina
More information8.04 Spring 2013 March 12, 2013 Problem 1. (10 points) The Probability Current
Prolem Set 5 Solutions 8.04 Spring 03 March, 03 Prolem. (0 points) The Proaility Current We wish to prove that dp a = J(a, t) J(, t). () dt Since P a (t) is the proaility of finding the particle in the
More informationarxiv: v1 [cs.sy] 24 Oct 2016
Filter-based regularisation for impulse response modelling Anna Marconato 1,**, Maarten Schoukens 1, Johan Schoukens 1 1 Dept. ELEC, Vrije Universiteit Brussel, Pleinlaan 2, 1050 Brussels, Belgium ** anna.marconato@vub.ac.be
More informationEconometrics I, Estimation
Econometrics I, Estimation Department of Economics Stanford University September, 2008 Part I Parameter, Estimator, Estimate A parametric is a feature of the population. An estimator is a function of the
More informationAdaptive Rejection Sampling with fixed number of nodes
Adaptive Rejection Sampling with fixed number of nodes L. Martino, F. Louzada Institute of Mathematical Sciences and Computing, Universidade de São Paulo, Brazil. Abstract The adaptive rejection sampling
More informationROBUST MINIMUM DISTANCE NEYMAN-PEARSON DETECTION OF A WEAK SIGNAL IN NON-GAUSSIAN NOISE
17th European Signal Processing Conference EUSIPCO 2009) Glasgow, Scotland, August 24-28, 2009 ROBUST MIIMUM DISTACE EYMA-PEARSO DETECTIO OF A WEAK SIGAL I O-GAUSSIA OISE Georgy Shevlyakov, Kyungmin Lee,
More informationExample: Bipolar NRZ (non-return-to-zero) signaling
Baseand Data Transmission Data are sent without using a carrier signal Example: Bipolar NRZ (non-return-to-zero signaling is represented y is represented y T A -A T : it duration is represented y BT. Passand
More informationA NONLINEAR TRANSFORMATION APPROACH TO GLOBAL ADAPTIVE OUTPUT FEEDBACK CONTROL OF 3RD-ORDER UNCERTAIN NONLINEAR SYSTEMS
Copyright 00 IFAC 15th Triennial World Congress, Barcelona, Spain A NONLINEAR TRANSFORMATION APPROACH TO GLOBAL ADAPTIVE OUTPUT FEEDBACK CONTROL OF RD-ORDER UNCERTAIN NONLINEAR SYSTEMS Choon-Ki Ahn, Beom-Soo
More informationUsing Multiple Kernel-based Regularization for Linear System Identification
Using Multiple Kernel-based Regularization for Linear System Identification What are the Structure Issues in System Identification? with coworkers; see last slide Reglerteknik, ISY, Linköpings Universitet
More informationOptimal input design for nonlinear dynamical systems: a graph-theory approach
Optimal input design for nonlinear dynamical systems: a graph-theory approach Patricio E. Valenzuela Department of Automatic Control and ACCESS Linnaeus Centre KTH Royal Institute of Technology, Stockholm,
More informationAsymptotic inference for a nonstationary double ar(1) model
Asymptotic inference for a nonstationary double ar() model By SHIQING LING and DONG LI Department of Mathematics, Hong Kong University of Science and Technology, Hong Kong maling@ust.hk malidong@ust.hk
More informationStationary or Non-Stationary Random Excitation for Vibration-Based Structural Damage Detection? An exploratory study
Stationary or Non-Stationary Random Excitation for Vibration-Based Structural Damage Detection? An exploratory study Andriana S. GEORGANTOPOULOU & Spilios D. FASSOIS Stochastic Mechanical Systems & Automation
More informationRandom Matrix Eigenvalue Problems in Probabilistic Structural Mechanics
Random Matrix Eigenvalue Problems in Probabilistic Structural Mechanics S Adhikari Department of Aerospace Engineering, University of Bristol, Bristol, U.K. URL: http://www.aer.bris.ac.uk/contact/academic/adhikari/home.html
More informationUsing Hankel structured low-rank approximation for sparse signal recovery
Using Hankel structured low-rank approximation for sparse signal recovery Ivan Markovsky 1 and Pier Luigi Dragotti 2 Department ELEC Vrije Universiteit Brussel (VUB) Pleinlaan 2, Building K, B-1050 Brussels,
More informationIdentificationwithStochasticSamplingTimeJitter
IdentificationwithStochasticSamplingTimeJitter Frida Eng 1, Fredrik Gustafsson. Dept of EE, Linköpings universitet, SE-58183 Linköping, Sweden Abstract This work investigates how stochastic sampling jitter
More informationLIKELIHOOD-BASED ESTIMATION OF PERIODICITIES IN SYMBOLIC SEQUENCES. Dept. of Mathematical Statistics, Lund University, Sweden
LIKELIHOOD-BASED ESTIMATION OF PERIODICITIES IN SYMBOLIC SEQUENCES Stefan Ingi Adalbjörnsson, Johan Swärd, and Andreas Jaobsson Dept. of Mathematical Statistics, Lund University, Sweden ABSTRACT In this
More informationOn Stochastic Adaptive Control & its Applications. Bozenna Pasik-Duncan University of Kansas, USA
On Stochastic Adaptive Control & its Applications Bozenna Pasik-Duncan University of Kansas, USA ASEAS Workshop, AFOSR, 23-24 March, 2009 1. Motivation: Work in the 1970's 2. Adaptive Control of Continuous
More informationCROSS-VALIDATION OF CONTROLLED DYNAMIC MODELS: BAYESIAN APPROACH
CROSS-VALIDATION OF CONTROLLED DYNAMIC MODELS: BAYESIAN APPROACH Miroslav Kárný, Petr Nedoma,Václav Šmídl Institute of Information Theory and Automation AV ČR, P.O.Box 18, 180 00 Praha 8, Czech Republic
More information2 Statistical Estimation: Basic Concepts
Technion Israel Institute of Technology, Department of Electrical Engineering Estimation and Identification in Dynamical Systems (048825) Lecture Notes, Fall 2009, Prof. N. Shimkin 2 Statistical Estimation:
More informationLTI Systems, Additive Noise, and Order Estimation
LTI Systems, Additive oise, and Order Estimation Soosan Beheshti, Munther A. Dahleh Laboratory for Information and Decision Systems Department of Electrical Engineering and Computer Science Massachusetts
More 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 informationAkaike criterion: Kullback-Leibler discrepancy
Model choice. Akaike s criterion Akaike criterion: Kullback-Leibler discrepancy Given a family of probability densities {f ( ; ψ), ψ Ψ}, Kullback-Leibler s index of f ( ; ψ) relative to f ( ; θ) is (ψ
More informationSubspace-based Identification
of Infinite-dimensional Multivariable Systems from Frequency-response Data Department of Electrical and Electronics Engineering Anadolu University, Eskişehir, Turkey October 12, 2008 Outline 1 2 3 4 Noise-free
More informationExtremogram and ex-periodogram for heavy-tailed time series
Extremogram and ex-periodogram for heavy-tailed time series 1 Thomas Mikosch University of Copenhagen Joint work with Richard A. Davis (Columbia) and Yuwei Zhao (Ulm) 1 Zagreb, June 6, 2014 1 2 Extremal
More informationML Estimation of Process Noise Variance in Dynamic Systems
ML Estimation of Process Noise Variance in Dynamic Systems Patrik Axelsson, Umut Orguner, Fredrik Gustafsson and Mikael Norrlöf {axelsson,umut,fredrik,mino} @isy.liu.se Division of Automatic Control Department
More informationsine wave fit algorithm
TECHNICAL REPORT IR-S3-SB-9 1 Properties of the IEEE-STD-57 four parameter sine wave fit algorithm Peter Händel, Senior Member, IEEE Abstract The IEEE Standard 57 (IEEE-STD-57) provides algorithms for
More informationInformation geometry of mirror descent
Information geometry of mirror descent Geometric Science of Information Anthea Monod Department of Statistical Science Duke University Information Initiative at Duke G. Raskutti (UW Madison) and S. Mukherjee
More informationWavelet Transform and Its Applications to Acoustic Emission Analysis of Asphalt Cold Cracking
Wavelet Transform and Its Applications to Acoustic Emission Analysis of Asphalt Cold Cracking Nesvijski, E., Marasteanu, M. Department of Civil Engineering University of Minnesota Minneapolis, MN 55455
More information