Linear dynamic filtering with noisy input and output 1
|
|
- Verity Franklin
- 5 years ago
- Views:
Transcription
1 Katholieke Universiteit Leuven Departement Elektrotechniek ESAT-SISTA/TR Linear dynamic filtering with noisy input and output 1 Ivan Markovsky and Bart De Moor 2 2 November 22 Published in the proc of 13th IFAC Symposium on System Identification (SYSID-23) pages This report is available by anonymous ftp from ftpesatkuleuvenacbe in the directory pub/sista/markovsky/reports/2-191psgz 2 KULeuven, Dept of Electrical Engineering (ESAT), Research group SCD (SISTA), Kasteelpark Arenberg 1, 31 Leuven-Heverlee, Belgium, Tel 32/16/ , Fax 32/16/ , WWW: ivanmarkovsky@esatkuleuvenacbe Research Council KUL: GOA Mefisto 666, several PhD/postdoc & fellow grants; Flemish Government: - FWO: PhD/postdoc grants, projects, G2499 (multilinear algebra), G472 (support vector machines), G1972 (power islands), G1413 (Identification and cryptography), G4913 (control for intensive care glycemia), G123 (QIT), research communities (ICCoS, ANMMM); - AWI: Bil Int Collaboration Hungary/Poland; - IWT: PhD Grants, Soft4s (softsensors), Belgian Federal Government: DWTC (IUAP IV 2 ( ) and IUAP V 22 (22 26), PODO II (CP/4: TMS and Sustainibility); EU: CAGE; ERNSI; Eureka 263 IMPACT; Eureka 2419 FliTE; Contract Research/agreements: Data4s, Electrabel, Elia, LMS, IPCOS, VIB;
2 Abstract We establish the equivalence between the optimal least-squares state estimator for a linear time-invariant dynamic system with noise corrupted input and output, and an appropriately modified Kalman filter The approach used is algebraic and the result shows that the noisy input/output filtering problem is not fundamentally different from the classical Kalman filtering problem The result is illustrated with a simulation example
3 LINEAR DYNAMIC FILTERING WITH NOISY INPUT AND OUTPUT Ivan Markovsky and Bart De Moor ESAT, SCD-SISTA, KULeuven, Kasteelpark Arenberg 1, B-31 Leuven-Heverlee, Belgium Tel: , Fax: Abstract: We establish the equivalence between the optimal least-squares state estimator for a linear time-invariant dynamic system with noise corrupted input and output, and an appropriately modified Kalman filter The approach used is algebraic and the result shows that the noisy input/output filtering problem is not fundamentally different from the classical Kalman filtering problem The result is illustrated with a simulation example Keywords: dynamic errors-in-variables model, Kalman filtering, optimal smoothing, total least squares 1 INTRODUCTION Optimal least-squares state estimation for linear dynamical systems is a well developed topic with many practical applications A central result is the Kalman filter In the discrete-time case, one considers the model x(t + 1) = Ax(t)+Bu(t)+w(t), y(t) = Cx(t)+Du(t)+v(t), (1) with x() = x and t =,1, Here u(t) R m, y(t) R l, and x(t) R n are the input, output, and state vectors at time instant t When the model matrices A, B, C, and D are known, the state x depends linearly on the input u and the output y, so that it can be estimated from the input/output (I/O) data via the leastsquares method The unknown process noise w and the measurement noise v play the role of equation errors in (1) Even though v and w are unknown, we assume that they are zero mean, white, Gaussian noises, with covariance matrix E w(t) v(t) ] w T (t + τ) v T (t + τ)] = Q(t) S(t) S T δ(τ) (t) R(t) The importance of the Kalman filter is that it solves the least-squares problem recursively and in real time Model (1), however, is in a certain sense asymmetric The process noise w acts as an unobserved input while the measurement noise v represents measurement error on y However, u is assumed to be noiseless! In the paper, we pose and answer the question: How should we modify the Kalman filter when both the input and the output of the system are noisy? The paradigm of treating the input and the output on an equal footing leads to the behavioral approach (Polderman and Willems, 1998) A derivation of the Kalman filter in the behavioral context is given in (Fagnani and Willems, 1997) In this paper, we consider the deterministic discretetime LTI state-space system x(t + 1) = Ax(t)+Bu(t), x() = x, y(t) = Cx(t)+Du(t), t =,1, together with the measurement errors model (2) u d (t) = u(t)+ũ(t), = y(t)+ỹ(t) (3) and refer to the model (2) together with the measurement errors model (3) as the noisy I/O model (see Figure 1) We note that (2 3) is an errors-in-variables (EIV) model in the sense of, eg, (Söderström et
4 al, 22; Zheng, 22), when the problem would be to estimate A, B, C, and D But since the purpose of this paper is to estimate x, knowing A, B, C, and D, we call it a noisy I/O model Kalman filter Section 6 further establishes the optimal estimates of the true input/output signals In Section 7, we confirm the results on an example and in Section 8, give conclusions x u ũ u d (A,B,C,D) Fig 1 Block scheme of the noisy I/O model A full analysis of a static EIV estimation problem, also known as the total least squares (TLS) problem, is given in the monograph (Van Huffel and Vandewalle, 1991) The dynamic equivalent of the TLS problem, see (Aoki and Yue, 197; De Moor and Roorda, 1994), is a system identification problem: given noisy I/O measurements, find the model (2) The dynamic TLS problem can be expressed as a static TLS problem with Toeplitz (or Hankel) structured data matrices The structured total least squares problem is treated in (De Moor, 1993; Lemmerling, 1999) Different type of estimation problem, in the EIV context, occurs when the model is exactly known and the state vector has to be estimated from the noisy I/O measurements We refer to this latter problem as the noisy I/O state estimation problem In (Guidorzi et al, 23; Diversi et al, 23), the noisy I/O state estimation problem has been considered, using the language of transfer functions The authors there claim that it is fundamentally different problem from the Kalman filtering problem and describe new recursive algorithms for its solution In this paper, we prove that the noisy I/O state estimation problem is not fundamentally different from the Kalman filtering problem Its solution boils down to the solution of a sequence of linear least-squares problems with a special structure coming from the state space equation (1) Our approach is linear algebraic We represent the model over a finite time horizon as a set of linear equations and apply standard linear algebra techniques for its analysis The optimal solution is shown to be a Kalman filter with correlated process and measurement noises The continuous-time version of the noisy I/O state estimation problem is treated in (Markovsky et al, 22), where a completion of squares approach is used and the solution is also shown to be a Kalman filter-type system In Section 2, we define the smoothing and filtering noisy I/O state estimation problems Two explicit block solutions of the smoothing problem are derived in Section 3 They are weighted least-squares problems In Section 4, we transform the noisy I/O model in the form (1) and define the Kalman filter for the resulting system as the modified Kalman filter In Section 5, we prove the deterministic equivalence of the estimates of the noisy I/O filter and the modified y d y ỹ 2 PROBLEM FORMULATION An on-line, recursive estimation procedure can be realized by a causal dynamical system called a filter Such a filter operates on previous and current measurements and produces an estimate of the to-beestimated signal for the current moment of time The problem of a filter synthesis is referred to as a filtering problem A smoothing problem is an estimation problem, in which on the basis of the available measurements, an estimate for the to-be-estimated signal is produced for the whole (past) period of observation We introduce some notation used in the rest of the paper A signal variable, without time index, denotes the vector obtained by stacking one over another the signal samples for the consecutive time instances For example, over the time horizon,1,,t f 1, the vector of the consecutive input samples u is defined as u := u T () u T (t f 1)] T and the vector of the consecutive state samples is x := x T () x T (t f 1) x T (t f )] T For a time indexed matrix sequence {V(t)} t f 1 t=, we denote by V, without time argument, the block matrix V := blk diag ( V(),,V(t f 1) ) Definition 1 (Optimal noisy I/O smoothing problem) Assume that the measurement errors ũ and ỹ are random, centered, normal, uncorrelated, and white with known covariance matrices cov ( ũ(t) ) =: Vũ(t), cov ( ỹ(t) ) =: Vỹ(t), (4) and that the initial condition x is unknown Then, given the matrices A, B, C, D, the optimal noisy I/O smoothing problem is defined as 1/2 û ] min Vũ ud 2 û,ŷ, ˆx Vỹ ŷ y d (5) ˆx(t + 1) = A ˆx(t)+Bû(t) st ŷ(t) = C ˆx(t)+Dû(t) for t =,1,,t f 1 The optimal smoothed state estimate ˆx(,t f ) is the solution of (5) Under the normality assumption for the noises, ˆx(,t f ) is the minimum variance estimate, the maximum likelihood estimate, and the conditional expectation estimate of the state x The equivalence is well known in the Kalman filter case (Willems, 22; Bryson and Ho, 1975; Anderson and Moore, 1979), and is shown in the noisy I/O case in (Guidorzi et al, 23)
5 Definition 2 (Optimal noisy I/O filtering problem) Given the model (2 3), satisfying the assumptions of Definition 1, the optimal noisy I/O filtering problem is to find a dynamical system, z(t + 1) = A f (t)z(t)+b f (t), (6) ˆx(t) = C f (t)z(t)+d f (t), such that ˆx(t) = ˆx(t,t + 1), where ˆx( ) is the solution of (6), ie, the optimal filtered state estimate, and ˆx(,t + 1) is the optimal smoothed state estimate with a time horizon t SMOOTHING BY BLOCK PROCESSING In this section, we write the optimal noisy I/O smoothing problem (5) as a weighted least-squares (WLS) problem This representation is used as a conceptual tool for the analysis and not as a means to carry out the actual computations needed for the estimation We represent the I/O dynamics of the system (2), over the time horizon,,t f 1, explicitly as where Γ := C CA CA t f 2 y = Γx + Tu, (7) H() H(1) H(), and T := H(t f 1) H(t f 2) H() The matrix Γ is an extended observability matrix and T is a Toeplitz matrix formed from the Markov parameters H() = D, H(t) = CA t 1 B, t = 1,,t f 1 Using (7), we see that the optimal noisy I/O smoothing problem (5) is a weighted least-squares problem 1 ( ) min Vũ 2 ud I 2 ˆx,û (8) Vỹ y d Γ T]ˆx û Alternatively, we represent the input/state/output dynamics of the system, over the time horizon,,t f 1, as y() y(1) y(t f 1) C A I x() C x(1) = A I x(2) + C x(t f ) A I } {{ } A D B D u() + B u(1) D u(t f 1) B } {{ } B Substituting y d ỹ for y and u d ũ for u (see (3)), we have y d () y d (1) u d () u d (1) +B = y d (t f 1) u d (t f 1) }{{} ȳ d I ũ() ũ(1) I ỹ() = A x+b + ỹ(1), ũ(t f 1) I ỹ(t f 1) } {{ } C or with the definition of the new variables ȳ d +Bu d = A x+bũ+cỹ (9) Using (9), the optimal noisy I/O smoothing problem is equivalent to the following problem 1/2 min Vũ u 2 ˆx, u, y Vỹ y (1) u st ȳ d +Bu d = A ˆx+B C], y which solution is given in Section 5 The solutions (8) and (1) of the noisy I/O smoothing problem are not recursive and thus not practical for large data sets See (Markovsky et al, 22) for recursive solution of an equivalent continuous-time problem In the discrete-time case, the recursive solution is given by two time-varying filters; one running backward in time and one running forward in time The forward recursion is defined by a time-varying Riccati equation The backward filter produces the optimal smoothed state estimate 4 THE MODIFIED KALMAN FILTER We convert the noisy I/O model (2 3) in the form (1) by substituting u d (t) ũ for u(t) and ỹ for y(t) (see (3)) in (2) x(t + 1) = Ax(t)+Bu d (t) Bũ(t) = Cx(t)+Du d (t) Dũ(t)+ỹ(t) and define (fake) process noise w and measurement noise v by w := Bũ and v := Dũ+ỹ The resulting system x(t + 1) = Ax(t)+Bu d (t)+w(t) = Cx(t)+Du d (t)+v(t) is in the form (1), where Q(t) S(t) B Vu (t) S T = (t) R(t) D I (11) ] T B V y (t) D I
6 We call the Kalman filter of the modified system (11), ie, the system z(t + 1) = A KF (t)z(t)+b KF (t), (12) ˆx(t) = C KF (t)z(t)+d KF (t), where A KF (t) = ( A K(t)C ), B KF (t) = B K(t)D, K(t) ], and C KF (t) = I P(t)C T( CP(t)C T + R(t) ) 1 C, D KF (t) = P(t)C T( CP(t)C T + R(t) ) 1 D I ], K(t) = ( AP(t)C T + S(t) )( CP(t)C T + R(t) ) 1, P(t + 1) = AP(t)A T ( AP(t)C T + S(t) ) ( CP(t)C T + R(t) ) 1( AP(t)C T + S(t) ) T + Q(t), the modified Kalman filter It recursively solves (9) for the last block entry of the unknown x The solution is in the sense of the weighted least-squares problem 1/2 min Ve ê 2 st ȳ d +Bu d = A ˆx+ê (13) ˆx,ê The variable ê accounts for the cumulative noise ũ ] e := B C] ỹ added to the equation and the covariance matrix of e is V e = BVũB T +CVỹC T When the measurement noise covariances Vũ(t) and Vũ(t) does not depend on t, one can replace the timevarying Kalman filter with the (suboptimal) timeinvariant filter, obtained by replacing P(t) in (12) with the steady-state solution P of the algebraic Riccati equation P = APA T ( A PC T +S )( C PC T +R ) 1 ( A PC T +S ) T +Q In the following section, we investigate the relation between the modified Kalman filter (12) state estimate and the noisy I/O filter state estimate 5 EQUIVALENCE OF THE MODIFIED KALMAN FILTER AND THE NOISY I/O FILTER Consider the linear system of equations (9) and the two solution methods (1) and (13) Denote ũ ] Vũ z := ȳ d +Bu d, δ :=, V ỹ δ :=, Vỹ u ˆδ :=, and D := B C] y We want to find a relation between the solutions of the following problems 1/2 min Ve ê 2 2 st z = A ˆx+ê, (14) ˆx,ê where V e = DV δ D T and min ˆx, ˆδ V 1/2 δ ˆδ 2 2 st z = A ˆx+D ˆδ (15) The first problem is a weighted least-squares problem and its solution is ˆx KF = (A T Ve 1 A ) 1 A T Ve 1 z (16) The noisy I/O estimation problem (15) is a minimumnorm type problem and its solution is (see Lemma 3 in the Appendix) ˆx = I]V e A A ] 1 z, (17) where A is a matrix which columns form a basis for the orthogonal complement of the range space of A We transform (16) and (17) by the change of variables Then and A := V 1/2 e A and z := V 1/2 e z ˆx KF = ( A T A ) 1 A T z, (18) ˆx = I] A A ] 1 z (19) ((19) follows from the identity A = Ve 1/2 A ) Now the question of the solutions equivalence is answered by Theorem 4, see the Appendix, which states that ˆx KF = ˆx Thus the two solutions are deterministically equal and the noisy I/O filtering problem is solved by the modified Kalman filter (12), ie A f = A KF, B f = B KF, C f = C KF, and D f = D KF 6 OPTIMAL ESTIMATION OF THE TRUE INPUT/OUTPUT SIGNALS Up to now we were interested in the optimal filtering in the sense of state estimation In this section, we show how the optimal estimates of the input and the output can be derived from the modified Kalman filter The solutions ê and ˆδ of (14) and (15), respectively, satisfy the following relation ê = V e A V e A A ] 1 z = DV δ D T A V e A A ] 1 z = D ˆδ This implies that the state estimate ˆx, the one-stepahead prediction z(t + 1), and the optimal input estimate û satisfy the equation z(t + 1) = A ˆx(t)+Bû(t) (2) Then we can find û exactly from ˆx and z(t + 1), obtained from the modified Kalman filter (12) In fact, (2) and the Kalman filter equations imply that û(t) = E(t)z(t)+F(t), where E(t) := VũD T ( CP(t)C T + R(t) ) 1 C and F(t) := I VũD T ( CP(t)C T + R(t) ) 1 D, VũD T( CP(t)C T + R(t) ) ] 1 The optimal output estimate is ŷ(t) = C ˆx(t) + Dû(t)
7 7 NUMERICAL EXAMPLE In this section, we verify numerically the equivalence of the solutions established in Section 5 The particular system, we use, is A = ], B = 1 ], C = ], and D = 5381 The time horizon is t f = 1, the initial state is x =, and the input signal is a normal white noise sequence with unit variance The input and the output of the system are corrupted by independent, centered, normal, white noises with variances var(ũ(t)) = 4 and var(ỹ(t)) = 4 for all t The estimate of the noisy I/O filter is computed directly from the definition, ie, we solve a sequence of optimal smoothing problems with increasing timehorizon Every smoothing problem is a weighted leastsquares problem that is solved explicitly according to (8) The last block entries of the obtained sequence of solutions form the noisy I/O filter state estimate We compare the noisy I/O filter estimate with the estimate of the modified Kalman filter (12) The experiment is carried out in MATLAB The state estimate ˆx KF obtained by the modified Kalman filter is up to the numerical errors equal to the state estimate ˆx f obtained by the noisy I/O filter, ˆx KF ˆx f = 57723e 15 This is the desired numerical verification of the theoretical result of the paper The absolute errors of estimation ˆx x 2, û u 2, ŷ y 2 for all estimation methods, discussed in the paper is given in Table 1 Table 1 Comparison of the absolute errors of the state, input, and output estimates for all methods and the noisy data (MKF modified Kalman filter) Method ˆx x 2 û u 2 ŷ y 2 optimal smoothing optimal filtering time-varying MKF time-invariant MKF noisy data CONCLUSIONS We considered optimal noisy I/O estimation problems for discrete-time LTI systems The filtering problem is solved via a modified Kalman filter The equivalence between the optimal noisy I/O filter and the modified Kalman filter is proven algebraically using explicit state-space representation of the system ACKNOWLEDGEMENTS Ivan Markovsky is a research assistant and Dr Bart De Moor is a full professor at the Katholieke Universiteit Leuven, Belgium Our research is supported by Research Council KUL: GOA Mefisto 666, several PhD/postdoc & fellow grants; Flemish Government: FWO: PhD/postdoc grants, projects, G2499 (multilinear algebra), G472 (support vector machines), G1972 (power islands), G1413 (Identification and cryptography), G4913 (control for intensive care glycemia), G123 (QIT), research communities (ICCoS, ANMMM); AWI: Bil Int Collaboration Hungary/Poland; IWT: PhD Grants, Soft4s (softsensors), Belgian Federal Government: DWTC (IUAP IV 2 ( ) and IUAP V 22 (22 26), PODO II (CP/4: TMS and Sustainibility); EU: CAGE; ERNSI; Eureka 263 IMPACT; Eureka 2419 FliTE; Contract Research/agreements: Data4s, Electrabel, Elia, LMS, IPCOS, VIB 9 REFERENCES Anderson, B D O and J B Moore (1979) Optimal Filtering Prentice Hall Aoki, M and P C Yue (197) On certain convergence questions in system identification SIAM J Control 8(2), Bryson, A and Y Ho (1975) Applied Optimal Control Hemisphere, Washington, DC De Moor, B (1993) Structured total least squares and L 2 approximation problems Lin Alg and Its Appl , De Moor, B and B Roorda (1994) L 2 -optimal linear system identification structured total least squares for SISO systems In: In the proceedings of the CDC pp Diversi, R, R Guidorzi and U Soverini (23) Algorithms for optimal errors-in-variables filtering Systems & Control Letters 48, 1 13 Fagnani, F and J C Willems (1997) Deterministic Kalman filtering in a behavioral framework Systems & Control Letters 32, Guidorzi, R, R Diversi and U Soverini (23) Optimal errors-in-variables filtering Automatica 39, Lemmerling, P (1999) Structured total least squares: analysis, algorithms and applications PhD thesis ESAT/SISTA, KU Leuven Markovsky, I, J C Willems and B De Moor (22) Continuous-time errors-in-variables filtering In: Proc of the Conference on Decision and Control pp Meyer, C D (2) Matrix Analysis and Applied Linear Algebra SIAM Polderman, J W and J C Willems (1998) Introduction to mathematical systems theory Springer- Verlag Söderström, T, U Soverini and K Mahata (22) Perspectives on errors-in-variables estimation for dynamic systems Signal Processing 82, Van Huffel, S and J Vandewalle (1991) The total least squares problem: Computational aspects and analysis SIAM, Philadelphia Willems, J C (22) Deterministic Kalman filtering J of Econometrics, to appear
8 Zheng, W X (22) A bias correction method for identification of linear dynamic errors-invariables models IEEE Trans on Aut Control 47(7), Appendix A SOLUTION OF THE OPTIMIZATION PROBLEM (15) AND PROOF OF THE STATE ESTIMATES EQUIVALENCE Lemma 3 Assuming that A is full rank and V δ is positive definite The solution of the minimum-norm type problem (15) is ˆδ Vδ D = T A ] V ˆx I e A A ] 1 z PROOF The Lagrangian of (15) is L( ˆx, ˆδ,λ) = ˆδ T V 1 δ ˆδ + λ T (A ˆx+Dδ z) The first order optimality condition L 1 = 2V ˆδ +D T λ = δ δ ˆδ = 1 2 V δd T λ, (A1) L ˆx = A T λ = λ = 2A λ, (A2) L λ = A ˆx+D ˆδ z = z = A ˆx+D ˆδ (A3) is a necessary and sufficient condition for a global minimum The matrix A is any matrix which columns form a basis for the range space of A T Substituting (A1) and (A2) in (A3), we have ] λ DV δ D T A A ] = z ˆx Using V e = DV δ D T and the assumption that A is full rank, the result follows Theorem 4 For a full rank matrix A R m n, with m greater then n, 1, (A T A) 1 A T = n (m n) I n ] A A] (A4) and the right-hand side is 1 I n ] A A = In ] U 2 U 1 Σ 1 V T] 1 U12 U 11 Σ 1 V T 1 = I n ] U 21 Σ 1 V T (A5) U 22 To find explicitly the inverse matrix in (A5), we use the formula for inverse of a block matrix (Meyer, 2, p123) B C D E where ] 1 = B 1 + B 1 CS 1 DB 1 B 1 CS 1 S 1 DB 1 S = E DB 1 C S 1 ], is the Schur complement of B in the matrix B C D E] For the block matrix in (A5), the Schur complement of U 12 is Then I n ] S = U 21 Σ 1 V T U 22 U 1 12 U 11Σ 1 V T = (U 21 U 22 U 1 12 U 11)Σ 1 V T A A] 1 = VΣ 1 1 (U 21 U 22 U 1 12 U 11) 1 U 22 U 1 12 I] (A6) Because of the orthogonality of U, U T 1 U 2 = U T 11U 12 +U T 21U 22 = Then U 22 U 1 12 = U T 21 U T 11 (A7) (U 21 U 22 U 1 12 U 11) 1 = (U 21 +U T 21 U T 11 U 11) 1 = ( U21 T (U 21 T U 21 +U11 T U 11) ) 1 = U T }{{} 21 (A8) U1 T U 1=I Substitution of the expressions of (A7) and (A8) into (A6) establishes the identity where A is a matrix which columns form a basis for the orthogonal complement of the range space of A, ie, Range(A ) = Null(A T ), and rank(a ) = m n PROOF In the proof, we use the SVD of the matrix A A = UΣV T Σ1 = U 1 U 2 ] V T = U 1 Σ 1 V T Partition U as follows n m n U = U 1 U 2 ] = U11 U 12 U 21 U 22 m n n The matrix U 2 satisfies Range(U 2 ) = Null(A T ) and rank(u 2 ) = m n, so it serves as a particular A The left-hand side of the desired identity (A4) is (A T A) 1 A T = VΣ 1 1 U 1 T = VΣ 1 1 U T 11 U21 T ],
Structured weighted low rank approximation 1
Departement Elektrotechniek ESAT-SISTA/TR 03-04 Structured weighted low rank approximation 1 Mieke Schuermans, Philippe Lemmerling and Sabine Van Huffel 2 January 2003 Accepted for publication in Numerical
More informationBlock-row Hankel Weighted Low Rank Approximation 1
Katholieke Universiteit Leuven Departement Elektrotechniek ESAT-SISTA/TR 03-105 Block-row Hankel Weighted Low Rank Approximation 1 Mieke Schuermans, Philippe Lemmerling and Sabine Van Huffel 2 July 2003
More informationInformation, Covariance and Square-Root Filtering in the Presence of Unknown Inputs 1
Katholiee Universiteit Leuven Departement Eletrotechnie ESAT-SISTA/TR 06-156 Information, Covariance and Square-Root Filtering in the Presence of Unnown Inputs 1 Steven Gillijns and Bart De Moor 2 October
More informationCompactly supported RBF kernels for sparsifying the Gram matrix in LS-SVM regression models
Compactly supported RBF kernels for sparsifying the Gram matrix in LS-SVM regression models B. Hamers, J.A.K. Suykens, B. De Moor K.U.Leuven, ESAT-SCD/SISTA, Kasteelpark Arenberg, B-3 Leuven, Belgium {bart.hamers,johan.suykens}@esat.kuleuven.ac.be
More informationOn Weighted Structured Total Least Squares
On Weighted Structured Total Least Squares Ivan Markovsky and Sabine Van Huffel KU Leuven, ESAT-SCD, Kasteelpark Arenberg 10, B-3001 Leuven, Belgium {ivanmarkovsky, sabinevanhuffel}@esatkuleuvenacbe wwwesatkuleuvenacbe/~imarkovs
More informationData Assimilation in 2D Magneto-Hydrodynamics Systems 1
Katholieke Universiteit Leuven Departement Elektrotechniek ESAT-SISTA/TR 0-7 Data Assimilation in D Magneto-Hydrodynamics Systems Oscar Barrero Mendoza, Dennis S. Bernstein, and Bart L.R. De Moor. September
More informationx 104
Departement Elektrotechniek ESAT-SISTA/TR 98-3 Identication of the circulant modulated Poisson process: a time domain approach Katrien De Cock, Tony Van Gestel and Bart De Moor 2 April 998 Submitted for
More informationSUBSPACE IDENTIFICATION METHODS
SUBSPACE IDENTIFICATION METHODS Katrien De Cock, Bart De Moor, KULeuven, Department of Electrical Engineering ESAT SCD, Kasteelpark Arenberg 0, B 300 Leuven, Belgium, tel: +32-6-32709, fax: +32-6-32970,
More informationStructured Matrices and Solving Multivariate Polynomial Equations
Structured Matrices and Solving Multivariate Polynomial Equations Philippe Dreesen Kim Batselier Bart De Moor KU Leuven ESAT/SCD, Kasteelpark Arenberg 10, B-3001 Leuven, Belgium. Structured Matrix Days,
More informationComparison of identification algorithms on the database DAISY
Proceedings of the 17th International Symposium on Mathematical Theory of Networks and Systems, Kyoto, Japan, July 24-28, 26 FrP11.2 Comparison of identification algorithms on the database DAISY Ivan Markovsky,
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 informationData assimilation for magnetohydrodynamics systems
Journal of Computational and Applied Mathematics 89 (00) 9 www.elsevier.com/locate/cam Data assimilation for magnetohydrodynamics systems O. Barrero Mendoza a,, B. De Moor a, D.S. Bernstein b a Department
More informationEquivalence of state representations for hidden Markov models
Systems & Control Letters 57 (2008) 410 419 www.elsevier.com/locate/sysconle Equivalence of state representations for hidden Markov models Bart Vanluyten, Jan C. Willems, Bart De Moor K.U.Leuven, ESAT/SCD(SISTA),
More informationDepartement Elektrotechniek ESAT-SISTA/TR 98- Stochastic System Identication for ATM Network Trac Models: a Time Domain Approach Katrien De Cock and Bart De Moor April 998 Accepted for publication in roceedings
More informationDynamic measurement: application of system identification in metrology
1 / 25 Dynamic measurement: application of system identification in metrology Ivan Markovsky Dynamic measurement takes into account the dynamical properties of the sensor 2 / 25 model of sensor as dynamical
More informationJoint Regression and Linear Combination of Time Series for Optimal Prediction
Joint Regression and Linear Combination of Time Series for Optimal Prediction Dries Geebelen 1, Kim Batselier 1, Philippe Dreesen 1, Marco Signoretto 1, Johan Suykens 1, Bart De Moor 1, Joos Vandewalle
More informationA missing data approach to data-driven filtering and control
1 A missing data approach to data-driven filtering and control Ivan Markovsky Abstract In filtering, control, and other mathematical engineering areas it is common to use a model-based approach, which
More informationLow-rank approximation and its applications for data fitting
Low-rank approximation and its applications for data fitting Ivan Markovsky K.U.Leuven, ESAT-SISTA A line fitting example b 6 4 2 0 data points Classical problem: Fit the points d 1 = [ 0 6 ], d2 = [ 1
More informationImproved initial approximation for errors-in-variables system identification
Improved initial approximation for errors-in-variables system identification Konstantin Usevich Abstract Errors-in-variables system identification can be posed and solved as a Hankel structured low-rank
More informationarxiv:quant-ph/ v2 17 Feb 2004
Graphical description of the action of local Clifford transformations on graph states Maarten Van den Nest Jeroen Dehaene and Bart De Moor Katholieke Universiteit Leuven ESAT-SCD Belgium. (Dated: February
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 informationLecture 2. Linear Systems
Lecture 2. Linear Systems Ivan Papusha CDS270 2: Mathematical Methods in Control and System Engineering April 6, 2015 1 / 31 Logistics hw1 due this Wed, Apr 8 hw due every Wed in class, or my mailbox on
More informationMaximum Likelihood Estimation and Polynomial System Solving
Maximum Likelihood Estimation and Polynomial System Solving Kim Batselier Philippe Dreesen Bart De Moor Department of Electrical Engineering (ESAT), SCD, Katholieke Universiteit Leuven /IBBT-KULeuven Future
More informationOn max-algebraic models for transportation networks
K.U.Leuven Department of Electrical Engineering (ESAT) SISTA Technical report 98-00 On max-algebraic models for transportation networks R. de Vries, B. De Schutter, and B. De Moor If you want to cite this
More informationOptimal trac light control for a single intersection Bart De Schutter and Bart De Moor y ESAT-SISTA, KU Leuven, Kardinaal Mercierlaan 94, B-3 Leuven (
Katholieke Universiteit Leuven Departement Elektrotechniek ESAT-SISTA/TR 97- Optimal trac light control for a single intersection Bart De Schutter and Bart De Moor Proceedings of the 997 International
More informationRECURSIVE SUBSPACE IDENTIFICATION IN THE LEAST SQUARES FRAMEWORK
RECURSIVE SUBSPACE IDENTIFICATION IN THE LEAST SQUARES FRAMEWORK TRNKA PAVEL AND HAVLENA VLADIMÍR Dept of Control Engineering, Czech Technical University, Technická 2, 166 27 Praha, Czech Republic mail:
More informationAdvancing density waves and phase transitions in a velocity dependent randomization traffic cellular automaton
Katholieke Universiteit Leuven Departement Elektrotechniek ESAT-SCD (SISTA) / TR 3- Advancing density waves and phase transitions in a velocity dependent randomization traffic cellular automaton Sven Maerivoet
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 informationDepartement Elektrotechniek ESAT-SISTA/TR About the choice of State Space Basis in Combined. Deterministic-Stochastic Subspace Identication 1
Katholieke Universiteit Leuven Departement Elektrotechniek ESAT-SISTA/TR 994-24 About the choice of State Space asis in Combined Deterministic-Stochastic Subspace Identication Peter Van Overschee and art
More informationA unified framework for EIV identification methods in the presence of mutually correlated noises
Preprints of the 19th World Congress The International Federation of Automatic Control A unified framework for EIV identification methods in the presence of mutually correlated noises Torsten Söderström
More informationThe characteristic equation and minimal state space realization of SISO systems in the max algebra
KULeuven Department of Electrical Engineering (ESAT) SISTA Technical report 93-57a The characteristic equation and minimal state space realization of SISO systems in the max algebra B De Schutter and B
More informationON THE REGULARIZATION OF CANONICAL CORRELATION ANALYSIS. Tijl De Bie. Bart De Moor
ON THE REGULARIZATION OF CANONICAL CORRELATION ANALYSIS Tijl De Bie Katholiee Universiteit Leuven ESAT-SCD Kasteelpar Arenberg 0 300 Leuven tijl.debie@esat.uleuven.ac.be Bart De Moor Katholiee Universiteit
More informationTHROUGHOUT the last few decades, the field of linear
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL 50, NO 10, OCTOBER 2005 1509 Subspace Identification of Hammerstein Systems Using Least Squares Support Vector Machines Ivan Goethals, Kristiaan Pelckmans, Johan
More informationData-driven signal processing
1 / 35 Data-driven signal processing Ivan Markovsky 2 / 35 Modern signal processing is model-based 1. system identification prior information model structure 2. model-based design identification data parameter
More informationPrimal-Dual Monotone Kernel Regression
Primal-Dual Monotone Kernel Regression K. Pelckmans, M. Espinoza, J. De Brabanter, J.A.K. Suykens, B. De Moor K.U. Leuven, ESAT-SCD-SISTA Kasteelpark Arenberg B-3 Leuven (Heverlee, Belgium Tel: 32/6/32
More informationEstimating Polynomial Structures from Radar Data
Estimating Polynomial Structures from Radar Data Christian Lundquist, Umut Orguner and Fredrik Gustafsson Department of Electrical Engineering Linköping University Linköping, Sweden {lundquist, umut, fredrik}@isy.liu.se
More information14 th IFAC Symposium on System Identification, Newcastle, Australia, 2006
1 th IFAC Symposium on System Identification, Newcastle, Australia, 6 THE REDUCED RANK TRANSFORM SQUARE ROOT FILTER FOR DATA ASSIMILATION S. Gillijns D. S. Bernstein B. De Moor Department of Electrical
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 informationMultirate MVC Design and Control Performance Assessment: a Data Driven Subspace Approach*
Multirate MVC Design and Control Performance Assessment: a Data Driven Subspace Approach* Xiaorui Wang Department of Electrical and Computer Engineering University of Alberta Edmonton, AB, Canada T6G 2V4
More informationAPPROXIMATE REALIZATION OF VALVE DYNAMICS WITH TIME DELAY
APPROXIMATE REALIZATION OF VALVE DYNAMICS WITH TIME DELAY Jan van Helvoirt,,1 Okko Bosgra, Bram de Jager Maarten Steinbuch Control Systems Technology Group, Mechanical Engineering Department, Technische
More informationLocal Convergence of Sequential Convex Programming for Nonconvex Optimization
Local Convergence of Sequential Convex Programming for Nonconvex Optimization Quoc Tran Dinh and Moritz Diehl *Department of Electrical Engineering (SCD-ESAT) and OPTEC, Katholieke Universiteit Leuven,
More informationThe QR decomposition and the singular value decomposition in the symmetrized max-plus algebra
K.U.Leuven Department of Electrical Engineering (ESAT) SISTA Technical report 96-70 The QR decomposition and the singular value decomposition in the symmetrized max-plus algebra B. De Schutter and B. De
More informationOptimization-Based Control
Optimization-Based Control Richard M. Murray Control and Dynamical Systems California Institute of Technology DRAFT v1.7a, 19 February 2008 c California Institute of Technology All rights reserved. This
More informationMatrix factorization and minimal state space realization in the max-plus algebra
KULeuven Department of Electrical Engineering (ESAT) SISTA Technical report 96-69 Matrix factorization and minimal state space realization in the max-plus algebra B De Schutter and B De Moor If you want
More informationSYSTEMTEORI - KALMAN FILTER VS LQ CONTROL
SYSTEMTEORI - KALMAN FILTER VS LQ CONTROL 1. Optimal regulator with noisy measurement Consider the following system: ẋ = Ax + Bu + w, x(0) = x 0 where w(t) is white noise with Ew(t) = 0, and x 0 is a stochastic
More informationECONOMETRIC METHODS II: TIME SERIES LECTURE NOTES ON THE KALMAN FILTER. The Kalman Filter. We will be concerned with state space systems of the form
ECONOMETRIC METHODS II: TIME SERIES LECTURE NOTES ON THE KALMAN FILTER KRISTOFFER P. NIMARK The Kalman Filter We will be concerned with state space systems of the form X t = A t X t 1 + C t u t 0.1 Z t
More informationThe Linear Dynamic Complementarity Problem is a special case of the Extended Linear Complementarity Problem B. De Schutter 1 B. De Moor ESAT-SISTA, K.
Katholieke Universiteit Leuven Departement Elektrotechniek ESAT-SISTA/TR 9-1 The Linear Dynamic Complementarity Problem is a special case of the Extended Linear Complementarity Problem 1 Bart De Schutter
More informationProblem set 5: SVD, Orthogonal projections, etc.
Problem set 5: SVD, Orthogonal projections, etc. February 21, 2017 1 SVD 1. Work out again the SVD theorem done in the class: If A is a real m n matrix then here exist orthogonal matrices such that where
More informationNonlinear State Estimation! Extended Kalman Filters!
Nonlinear State Estimation! Extended Kalman Filters! Robert Stengel! Optimal Control and Estimation, MAE 546! Princeton University, 2017!! Deformation of the probability distribution!! Neighboring-optimal
More informationFIR Filters for Stationary State Space Signal Models
Proceedings of the 17th World Congress The International Federation of Automatic Control FIR Filters for Stationary State Space Signal Models Jung Hun Park Wook Hyun Kwon School of Electrical Engineering
More informationOn Spectral Factorization and Riccati Equations for Time-Varying Systems in Discrete Time
On Spectral Factorization and Riccati Equations for Time-Varying Systems in Discrete Time Alle-Jan van der Veen and Michel Verhaegen Delft University of Technology Department of Electrical Engineering
More informationELEC system identification workshop. Behavioral approach
1 / 33 ELEC system identification workshop Behavioral approach Ivan Markovsky 2 / 33 The course consists of lectures and exercises Session 1: behavioral approach to data modeling Session 2: subspace identification
More informationTransductively Learning from Positive Examples Only
Transductively Learning from Positive Examples Only Kristiaan Pelckmans and Johan A.K. Suykens K.U.Leuven - ESAT - SCD/SISTA, Kasteelpark Arenberg 10, B-3001 Leuven, Belgium Abstract. This paper considers
More informationSVD-based optimal ltering with applications to noise reduction in speech signals Simon Doclo ESAT - SISTA, Katholieke Universiteit Leuven Kardinaal Me
Departement Elektrotechniek ESAT-SISTA/TR 999- SVD-based Optimal Filtering with Applications to Noise Reduction in Speech Signals Simon Doclo, Marc Moonen April, 999 Internal report This report is available
More informationDesign of FIR Smoother Using Covariance Information for Estimating Signal at Start Time in Linear Continuous Systems
Systems Science and Applied Mathematics Vol. 1 No. 3 2016 pp. 29-37 http://www.aiscience.org/journal/ssam Design of FIR Smoother Using Covariance Information for Estimating Signal at Start Time in Linear
More informationRiccati difference equations to non linear extended Kalman filter constraints
International Journal of Scientific & Engineering Research Volume 3, Issue 12, December-2012 1 Riccati difference equations to non linear extended Kalman filter constraints Abstract Elizabeth.S 1 & Jothilakshmi.R
More informationBasic Concepts in Data Reconciliation. Chapter 6: Steady-State Data Reconciliation with Model Uncertainties
Chapter 6: Steady-State Data with Model Uncertainties CHAPTER 6 Steady-State Data with Model Uncertainties 6.1 Models with Uncertainties In the previous chapters, the models employed in the DR were considered
More informationErrors-in-variables identification through covariance matching: Analysis of a colored measurement noise case
008 American Control Conference Westin Seattle Hotel Seattle Washington USA June -3 008 WeB8.4 Errors-in-variables identification through covariance matching: Analysis of a colored measurement noise case
More informationBlind Identification of FIR Systems and Deconvolution of White Input Sequences
Blind Identification of FIR Systems and Deconvolution of White Input Sequences U. SOVERINI, P. CASTALDI, R. DIVERSI and R. GUIDORZI Dipartimento di Elettronica, Informatica e Sistemistica Università di
More informationOptimal traffic light control for a single intersection
KULeuven Department of Electrical Engineering (ESAT) SISTA Technical report 97- Optimal traffic light control for a single intersection B De Schutter and B De Moor If you want to cite this report, please
More informationKatholieke Universiteit Leuven
Katholieke Universiteit Leuven Departement Elektrotechniek ESA-SISA/R 24-167 An instrumental variable method for adaptive feedback cancellation in hearing aids 1 Ann Spriet 2, Ian Proudler 3,Marc Moonen
More informationImplementation of the regularized structured total least squares algorithms for blind image deblurring
Linear Algebra and its Applications 391 (2004) 203 221 wwwelseviercom/locate/laa Implementation of the regularized structured total least squares algorithms for blind image deblurring N Mastronardi a,,1,
More informationOn instrumental variable-based methods for errors-in-variables model identification
On instrumental variable-based methods for errors-in-variables model identification Stéphane Thil, Marion Gilson, Hugues Garnier To cite this version: Stéphane Thil, Marion Gilson, Hugues Garnier. On instrumental
More informationA Behavioral Approach to GNSS Positioning and DOP Determination
A Behavioral Approach to GNSS Positioning and DOP Determination Department of Communications and Guidance Engineering National Taiwan Ocean University, Keelung, TAIWAN Phone: +886-96654, FAX: +886--463349
More informationarxiv: v1 [math.ag] 10 Oct 2011
Are diverging CP components always nearly proportional? Alwin Stegeman and Lieven De Lathauwer November, 8 Abstract arxiv:.988v [math.ag] Oct Fitting a Candecomp/Parafac (CP) decomposition (also nown as
More informationLessons in Estimation Theory for Signal Processing, Communications, and Control
Lessons in Estimation Theory for Signal Processing, Communications, and Control Jerry M. Mendel Department of Electrical Engineering University of Southern California Los Angeles, California PRENTICE HALL
More informationControl of traffic with anticipative ramp metering
Van den Berg, Bellemans, De Schutter, Hellendoorn 1 Control of traffic with anticipative ramp metering Submission date: Word count: 5033 words + (4 figures)*(250 words) = 6033 words Authors: M. van den
More informationCALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems. CDS 110b
CALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems CDS 110b R. M. Murray Kalman Filters 14 January 2007 Reading: This set of lectures provides a brief introduction to Kalman filtering, following
More informationSteady State Kalman Filter
Steady State Kalman Filter Infinite Horizon LQ Control: ẋ = Ax + Bu R positive definite, Q = Q T 2Q 1 2. (A, B) stabilizable, (A, Q 1 2) detectable. Solve for the positive (semi-) definite P in the ARE:
More informationSystem Identification by Nuclear Norm Minimization
Dept. of Information Engineering University of Pisa (Italy) System Identification by Nuclear Norm Minimization eng. Sergio Grammatico grammatico.sergio@gmail.com Class of Identification of Uncertain Systems
More informationLecture 19 Observability and state estimation
EE263 Autumn 2007-08 Stephen Boyd Lecture 19 Observability and state estimation state estimation discrete-time observability observability controllability duality observers for noiseless case continuous-time
More informationPerformance Analysis of an Adaptive Algorithm for DOA Estimation
Performance Analysis of an Adaptive Algorithm for DOA Estimation Assimakis K. Leros and Vassilios C. Moussas Abstract This paper presents an adaptive approach to the problem of estimating the direction
More informationVirtual Reference Feedback Tuning for non-linear systems
Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference 25 Seville, Spain, December 2-5, 25 ThA9.6 Virtual Reference Feedback Tuning for non-linear systems
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 informationPrediction, filtering and smoothing using LSCR: State estimation algorithms with guaranteed confidence sets
2 5th IEEE Conference on Decision and Control and European Control Conference (CDC-ECC) Orlando, FL, USA, December 2-5, 2 Prediction, filtering and smoothing using LSCR: State estimation algorithms with
More informationTHE GEOMETRY OF MULTIVARIATE POLYNOMIAL DIVISION AND ELIMINATION
THE GEOMETRY OF MULTIVARIATE POLYNOMIAL DIVISION AND ELIMINATION KIM BATSELIER, PHILIPPE DREESEN, AND BART DE MOOR Abstract. Multivariate polynomials are usually discussed in the framework of algebraic
More informationCONTROL SYSTEMS, ROBOTICS AND AUTOMATION - Vol. XIII - Nonlinear Observers - A. J. Krener
NONLINEAR OBSERVERS A. J. Krener University of California, Davis, CA, USA Keywords: nonlinear observer, state estimation, nonlinear filtering, observability, high gain observers, minimum energy estimation,
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 informationSequential State Estimation (Crassidas and Junkins, Chapter 5)
Sequential State Estimation (Crassidas and Junkins, Chapter 5) Please read: 5.1, 5.3-5.6 5.3 The Discrete-Time Kalman Filter The discrete-time Kalman filter is used when the dynamics and measurements are
More informationA Perturbation Analysis using Second Order Cone Programming for Robust Kernel Based Regression
A Perturbation Analysis using Second Order Cone Programg for Robust Kernel Based Regression Tillmann Falc, Marcelo Espinoza, Johan A K Suyens, Bart De Moor KU Leuven, ESAT-SCD-SISTA, Kasteelpar Arenberg
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 informationFor extensive reviews of the total least squares (TLS)
Advanced Review Ivan Markovsky, 1 Diana M. Sima, and Sabine Van Huffel Recent advances in total least squares approaches for solving various errorsin-variables modeling problems are reviewed, with emphasis
More informationA note on the characteristic equation in the max-plus algebra
K.U.Leuven Department of Electrical Engineering (ESAT) SISTA Technical report 96-30 A note on the characteristic equation in the max-plus algebra B. De Schutter and B. De Moor If you want to cite this
More informationObservability and state estimation
EE263 Autumn 2015 S Boyd and S Lall Observability and state estimation state estimation discrete-time observability observability controllability duality observers for noiseless case continuous-time observability
More informationSparsity in system identification and data-driven control
1 / 40 Sparsity in system identification and data-driven control Ivan Markovsky This signal is not sparse in the "time domain" 2 / 40 But it is sparse in the "frequency domain" (it is weighted sum of six
More informationSIMON FRASER UNIVERSITY School of Engineering Science
SIMON FRASER UNIVERSITY School of Engineering Science Course Outline ENSC 810-3 Digital Signal Processing Calendar Description This course covers advanced digital signal processing techniques. The main
More informationSubspace Identification A Markov Parameter Approach
Subspace Identification A Markov Parameter Approach Nelson LC Chui JM Maciejowski Cambridge University Engineering Deptartment Cambridge, CB2 1PZ, England 22 December 1998 Technical report CUED/F-INFENG/TR337
More informationCALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems. CDS 110b
CALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems CDS 110b R. M. Murray Kalman Filters 25 January 2006 Reading: This set of lectures provides a brief introduction to Kalman filtering, following
More informationDepartement Elektrotechniek ESAT-SISTA/TR Dynamical System Prediction: a Lie algebraic approach for a novel. neural architecture 1
Katholieke Universiteit Leuven Departement Elektrotechniek ESAT-SISTA/TR 1995-47 Dynamical System Prediction: a Lie algebraic approach for a novel neural architecture 1 Yves Moreau and Joos Vandewalle
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 informationStatistics 910, #15 1. Kalman Filter
Statistics 910, #15 1 Overview 1. Summary of Kalman filter 2. Derivations 3. ARMA likelihoods 4. Recursions for the variance Kalman Filter Summary of Kalman filter Simplifications To make the derivations
More informationX t = a t + r t, (7.1)
Chapter 7 State Space Models 71 Introduction State Space models, developed over the past 10 20 years, are alternative models for time series They include both the ARIMA models of Chapters 3 6 and the Classical
More informationThe norms can also be characterized in terms of Riccati inequalities.
9 Analysis of stability and H norms Consider the causal, linear, time-invariant system ẋ(t = Ax(t + Bu(t y(t = Cx(t Denote the transfer function G(s := C (si A 1 B. Theorem 85 The following statements
More informationStatistical and Adaptive Signal Processing
r Statistical and Adaptive Signal Processing Spectral Estimation, Signal Modeling, Adaptive Filtering and Array Processing Dimitris G. Manolakis Massachusetts Institute of Technology Lincoln Laboratory
More informationCANONICAL LOSSLESS STATE-SPACE SYSTEMS: STAIRCASE FORMS AND THE SCHUR ALGORITHM
CANONICAL LOSSLESS STATE-SPACE SYSTEMS: STAIRCASE FORMS AND THE SCHUR ALGORITHM Ralf L.M. Peeters Bernard Hanzon Martine Olivi Dept. Mathematics, Universiteit Maastricht, P.O. Box 616, 6200 MD Maastricht,
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 informationThe Kalman Filter ImPr Talk
The Kalman Filter ImPr Talk Ged Ridgway Centre for Medical Image Computing November, 2006 Outline What is the Kalman Filter? State Space Models Kalman Filter Overview Bayesian Updating of Estimates Kalman
More informationUnbiased minimum variance estimation for systems with unknown exogenous inputs
Unbiased minimum variance estimation for systems with unknown exogenous inputs Mohamed Darouach, Michel Zasadzinski To cite this version: Mohamed Darouach, Michel Zasadzinski. Unbiased minimum variance
More informationREGLERTEKNIK AUTOMATIC CONTROL LINKÖPING
Recursive Least Squares and Accelerated Convergence in Stochastic Approximation Schemes Lennart Ljung Department of Electrical Engineering Linkping University, S-581 83 Linkping, Sweden WWW: http://www.control.isy.liu.se
More informationSubspace angles and distances between ARMA models
ÔÖØÑÒØ ÐØÖÓØÒ Ë̹ËÁËÌ»ÌÊ ß ËÙ Ô ÒÐ Ò ØÒ ØÛÒ ÊÅ ÑÓÐ ÃØÖÒ Ó Ò ÖØ ÅÓÓÖ ÅÖ ÈÙÐ Ò ÈÖÓÒ Ó Ø ÓÙÖØÒØ ÁÒØÖÒØÓÒÐ ËÝÑÔÓ ÙÑ Ó ÅØÑØÐ ÌÓÖÝ Ó ÆØÛÓÖ Ò ËÝ ØÑ ÅÌÆË µ ÈÖÔÒÒ ÖÒ ÂÙÒ ß Ì ÖÔÓÖØ ÚÐÐ Ý ÒÓÒÝÑÓÙ ØÔ ÖÓÑ ØÔº غÙÐÙÚÒºº
More information