EXACT QUANTIFICATION OF VARIANCE ERROR
|
|
- Julius Dalton
- 6 years ago
- Views:
Transcription
1 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. FAX: Department of Sensors, Signals and Systems Automatic Control), The Royal Institute of Technology, S Stockholm, Sweden, FAX: Abstract: This paper establishes a method for quantifying variance error in cases where the input spectral density has a rational factorisation. Compared to previous work which has involved asymptotic-in-model-order techniques and yielded expressions which are only approximate for finite orders, the quantifications provided here are exact for finite model order, although they still only apply asymptotically in the observed data length. The key principle employed here is the use of a reproducing kernel in order to characterise the model class, and this implies a certain geometric-type approach to the problem. Keywords: Parameter Estimation, Maximum Likelihood Estimators.. ITRODUCTIO There has been much work over the last several years directed at quantifying the noise induced error variance error) in frequency response estimates L.Ljung and Z.D.Yuan 985, L.Ljung 985, P.M.J. Van den Hof et al. 995, Wahlberg 99, Wahlberg 994, inness et al. 999). In these contributions, in the interests of tractable analysis, a common theme has been to consider the limiting value of the variance error as the model order grows, and then assume that this same value can be approximately used for finite model order. This approach leaves open the question of the accuracy of the ensuing approximate quantification, and indeed recent works P.M.J. Van den Hof et al. 995, Wahlberg 99, Wahlberg 994, inness et al. 999) have specifically addressed the issue by using certain orthonormal basis techniques designed to maximise the finite model order accuracy. The paper here takes a different approach, and shows how in some circumstances it is possible to quantify the variance error in a manner that is exact even for finite model order, although it is still asymptotic in the data length. The key technique used to achieve this is a geometric one, in which recognition of certain subspace invariants reproducing kernels) provides a means for the exact quantification. The importance of this issue of non-asymptotic in model order variance quantification has been recognised by other authors Xie and Ljung 000), who have employed an This work was supported by the Australian Research Council and the Centre for Integrated Dynamics and Control CIDAC). It was completed while the first author was visiting the Department of Sensors, Signals and Systems Automatic Control), The Royal Institute of Technology, Stockholm, Sweden, analysis technique that is quite different to that of this work. The penultimate section of this paper profiles the results of this approach versus the one taken here.. PROBLEM SETTIG The estimation problem considered here is one in which observed input-output data y t, u t is generated according to y t = Gq)u t + e t. Here e t is a zero-mean white noise sequence that satisfies Ee t = σ. As well, Gq) is assumed to be an n th order rational transfer function with poles at ξ 0,, ξ n, and u t is taken to be a second order stationary process with associated spectral density 0 < Φ u ω) <. For the purposes of estimating the dynamics Gq), it is supposed that the following parameterised model structure is used where y t = Gq, θ)u t + e t ) n Gq, θ) = A q) θ k q k, Aq) = n q ξ k ). Here, the fixed poles ξ 0,, ξ m are meant to be chosen via a-priori knowledge Wahlberg 99, Wahlberg 994, Heuberger et al. 995, P.M.J. Van den Hof et al. 995), and the special case of ξ k = 0 renders ) as the common FIR model structure. With this in mind, an estimate θ of θ [θ 0,, θ n ] T based on observations of y t and u t may be found via a least squares approach:
2 θ = arg min V θ) ) θ V θ) y t Gq, θ)u t ). t= Furthermore, since Gq) is assumed to be n th order, then there exists a true parameter vector θ such that Gq, θ ) = Gq) and hence with ψ t = A q)λq)u t, Λq) [q n,, q, ] T 3) it is clear that the expression θ θ = φ t φ T t t= ) φ t e t 4) t= is a quantification of the parameter space estimation error. Furthermore, since the estimated frequency response Ge jω, θ ) depends linearly on the parameter estimate θ according to Ge jω, θ ) = A e jω )Λ T e jω ) θ then the frequency response estimation error may be quantified from 4) as Ge jω, θ ) Ge jω )=A e jω )Λ T e jω ) ) φ t φ T t φ t e t 5) t= t= Of course, this expression 5) depends on the exact noise realisation e t which cannot be known, so in order to quantify the estimation error Ge jω, θ ) Ge jω ) it is usual to consider its average over the ensemble of possible noise realisations. That is, consideration is given to ), which by 5) and assuming openloop data collection, may be computed as where ) = σ R Λ e jω )R Λejω ) Ae jω ) 6) φ t φ T t t= and represents the operation of conjugate transpose. ow, while 6) gives an exact quantification of the noise induced estimation error, it does not expose very much about how experiment design choices or other factors) might affect it. In consideration of this, the approach pioneered in works such as Whittle 953, Hannan and icholls 977, L.Ljung and Z.D.Yuan 985) is to notice that according to the assumptions made on the input R lim R Φu = T n A 7) Λe jλ )Λ e jλ ) Φ uλ) Ae jλ dλ. 8) ) Here, the notation T n Φ u ) is chosen to denote a Toeplitz matrix, of size n n, that is completely characterised by the spectral density function Φ u. In fact, a well known feature of such matrices is that for large n, then ) A T n T n Φu A Φ u where the approximation sign means that the Hilbert Schmidt matrix norm Golub and Loan 989) between the two quantities can be made arbitrarily small for arbitrarily large model roder n Grenander and Szegö 958). In this case, for the same large n it is inviting to approximate ) from 6) as Finally, note that ) σ Ae jω ) A Λ e jω )T n Λe jω ). 9) nπ n Λ e jω )T n Φ u A Φ u ) Λe jω ) = Λ e jω )Λe jλ ) Aejλ ) Φ u λ) τ= dλ = τ ) c τ e jωτ 0) where the c τ are the Fourier co-efficients of A /Φ u : c τ Ae jλ ) Φ u λ) e jλτ dλ. ) Since the Fourier series in 0) can be expected to converge to the function A /Φ u determining its coefficients via ), it is arguable via 9), 0) that for large model order n ) σ n Ae jω ) Aejω ) Φ u ω) = n σ Φ u ω). ) In contrast to 6), this approximate expression ) very clearly exposes the key factors that contribute to the noise induced estimation error. As such, its use has become a fundamental part of the science of System Identification Ljung 999) since it was first presented for the FIR case ξ k = 0) in L.Ljung and Z.D.Yuan 985) and then extended to much more general model structures in L.Ljung 985, P.M.J. Van den Hof et al. 995, Wahlberg 99, Wahlberg 994). However, as argued in inness et al. 999, inness and Hjalmarsson 00a, inness and Hjalmarsson 00b), this reliance can be misplaced. For example, the Achillee s Heel of ) is that it clearly depends on the model order
3 n being large in order for the approximating steps to be accurate. In recognition of this, the work inness et al. 999) has suggested that a quantification that can be more reliable for low model order n is ) σ n Φ u ω) ξ k e jω ξ k 3) ote that as pointed out in inness et al. 999, Xie and Ljung 000), the greater the proportion of the poles ξ k located at the origin, the closer the approximation 3) is to the seminal one ), and when all the poles are at the origin, then ) becomes a special case of 3). evertheless, the approximation 3) still depends on high model order n in order to obtain high accuracy. In recognition of this last point, the purpose of this paper is to derive an expression for lim ) that is exact for finite model order n. 3. MATHEMATICAL BACKGROUD The key tool used in this paper is what might be called a geometric one, and it depends on the idea of a Reproducing Kernel. To explain this, note that any frequency response of interest can be assumed to lie in a certain Hilbert Space H, where any geometric ideas such as orthogonality are obtained from the inner product between two functions f, g H which is defined as f, g µ fλ)gλ) µλ) dλ 4) with µλ) being some positive definite weighting function. ow, suppose that f is an element of an n-dimensional subspace X n defined by certain orthonormal basis elements ϕ 0,, ϕ n as f X n Span ϕ 0,, ϕ n. Then it is a consequence of the Riesz Representation TheoremRiesz and Sz.-agy 955, Rudin 966) that for any fixed ω, a further element K n λ, ω) X n exists such that fω) = fλ), K n λ, ω) µ f X n. 5) This reproducing kernel will be vital to the arguments of this paper, so it is important to present its fundamental features. Firstly it is Hermitian Symmetric in that, since for any fixed λ the kernel K n ω, λ) X n, then integration in the inner product is with respect to the common variable) K n ω, λ) = K n σ, λ), K m σ, ω) µ = K n σ, ω), K n σ, λ) µ = K n λ, ω). This then implies that K n ω, λ) is the unique element in X n that has the property 5), since if another function H n ω, σ) also satisfied 5), then it would hold that H n ω, λ) = H n ω, λ) = H n σ, λ), K n σ, ω) µ = K n σ, ω), H n σ, λ) µ = K n λ, ω). Finally, although the reproducing kernel is unique, there may of course) be many different ways of expressing it. One obvious one uses the orthonormal basis ϕ k for X n to express the quantity as n K n λ, ω) = ϕ k λ)ϕ k ω). 6) This formulation can be verified by noting that if X n f = τ c τϕ τ for some constants c τ, then n n c τ ϕ τ λ), ϕ k λ)ϕ k ω) = τ=0 n n c τ τ=0 ϕ k ω) ϕ τ λ), ϕ k λ) µ = n c τ ϕ τ ω) = fω). 7) τ=0 4. MAI RESULT With these preliminary ideas in mind, then as just argued via 6) and 9) lim ) = σ Ae jω ) Λ e jω )Tn Φu A Λe jω ). 8) The challenge then is to quantify this quadratic form in an exact manner, as opposed to the previously profiled approximation schemes based on asymptotic analysis. This ambition may be achieved by the following Theorem, which is the main result of the paper. Theorem 4.. Suppose that Φ u ω) has a finite dimensional rational spectral factorisation Φ u ω) = κ Ψe jω )Ψe jω ), Then Ψz) ν t=0 µ z βt z α t Ae jω ) Λ e jω )Tn Φu A Λe jω ) ). = n κ ϕ k e jω ) 9) where ϕ k z) is such that Ψz)ϕ 0 z),, Ψz)ϕ n z) is an orthonormal set of functions satisfying
4 X n Span Ψz)ϕ 0 z),, Ψz)ϕ n z) = Span Ψz) zn Az),, Ψz). 0) Az) Proof: ote that as a simple consequence of the definition 8) Ψe jω ) Ae jω ) Λ ne jω )Tn Φu A Λ n e jλ ) κ Ψejλ )Ψe jλ ) Ae jλ )Ae jλ ) Λ ne jλ ) dλ = Ψejω ) Ae jω ) Λ ne jω ). Therefore, the reproducing kernel K n λ, ω) for the space X n defined in 0) is given by the following quadratic form: K n λ, ω) = κ Ψejλ )Ψe jω ) Ae jλ )Ae jω ) Λ e jω )Tn Φu A Λe jλ ). ) However, as established in 7), if the set Ψz)ϕ 0 z),, Ψz)ϕ n z) is an orthonormal basis for X n then n K n λ, ω)= Ψe jλ )ϕ k e jλ )Ψe jω )ϕ k e jω ) n =Ψe jλ )Ψe jω ) ϕ k e jλ )ϕ k e jω ). ) Equating ) with ) when λ = ω then gives 9). On account of 8), this result immediately gives the following exact in n) formulation for variance error lim ) = σ n κ ϕ k e jω ) 3) which implies the following error quantification, which is not asymptotic in model order n: ) σ n κ ϕ k e jω ). 4) Of course, for this to produce an explicit quantification it is necessary to compute the basis Ψz)ϕ k z) for the space X n defined in 0). The following examples will demonstrate how this may be achieved analytically. 5. FIRST ORDER EXAMPLE In order to illustrate the preceding ideas, consider the simplest and lowest order example possible in which the observed data is obtained from the first order system ξ0 y t = u t + e t q ξ 0 where the spectral factor Ψz) of Φ u is of the form Ψz) = z β z α. Then an elementary calculation shows that Ψz) z ξ 0 ) = + αξ 0) + β ) βα + ξ 0 ) α ) αξ 0 ) ξ 0 ) and hence the unit norm version Ψz)ϕ 0 z) of Ψz)z ξ 0 ) is [ α ] ) αξ 0 ) Ψz)ϕ 0 z) = + αξ 0 ) + β ) βα + ξ 0 ) Therefore, according to 3) Ψz) z ξ 0 ). lim VarGejω, θ n ) = σ κ α ) ξ 0 ) αξ 0 ) + αξ 0 ) + β ) βα + ξ 0 ) e jω ξ 0. 5) For the case of α = 0.5, β = 0, ξ 0 = 0.9, κ =, σ = e 4 the results of using the ensuing quantification VarGe jω, θ n ) σ α ) ξ 0 ) αξ 0 ) κ + αξ 0 ) + β ) βα + ξ 0 ) e jω ξ 0. 6) are shown in figure. ote that the expression 6) is es Estimate variability vs exact quantification and asymptotic approximation Sample Variability Exact Quantification Asymptotic Approximation Fig.. Variability of Fixed Denominator Estimate - First Order Example. True variability vs. theoretically derived approximations. Solid line is Monte Carlo estimate of true variability, dashed line is quantification 6), dash-dot line is asymptotic-based approximation ). sentially) exact, which provides one example validating the claim of this paper that it is possible to derive accurate variance error quantifications that apply for arbitrary low model order. Here, and in what follows, the true variability is estimated via the sample average over 0000 different realisations of input and measurement noise, and is illustrated as a solid line. The exact quantification, such as given by 6) is shown as a dashed line, and the well known approximation ) is shown as a dash-dot line.
5 6. SECOD ORDER EXAMPLE Suppose we increase the model complexity in the previous example to that of Gq) = ξ 0) ξ ) q ξ 0 )q ξ ). In order to compute the asymptotic variance associated with this model structure, it is necessary to compute the orthonormal basis Ψϕ 0, Ψϕ that spans the same space as Ψz)z ξ 0 ), Ψz)z ξ ). In fact, this is quite difficult, primarily because of the restriction that all elements in this space have the same zeroes as Ψz). evertheless it can be achieved via the well known Gramm Schmidt procedure. Specifically, we have already established that [ Ψz)ϕ 0 z) = α ) αξ 0 ) + αξ 0 ) + β ) βα + ξ 0 ) Ψz) z ξ 0 ) ] is of unit norm, and hence a basis element fz) that is orthonormal to Ψz)ϕ 0 z) and also such that SpanΨz)ϕ 0, f = SpanΨz)z ξ 0 ), Ψz ξ ) is given by fz)=ψz) z ξ ) Ψz) z ξ ), ϕ 0Ψ Ψz)ϕ 0 z). While this looks innocent enough, it evaluates to where fz) = z β)γz + δ) z ξ 0 )z ξ )z α) ξ 0 )β α)αβ ) γ [ + αξ 0 ) + β ) βα + ξ 0 )]αξ ) ξ 0, δ + and hence, f = ξ 0 )β α)αβ )ξ [ + αξ 0 ) + β ) βα + ξ 0 )]αξ ) 7) ξ 0 β)ξ 0 γ + δ) βξ 0 )γ + δξ 0 ) ξ 0 ξ )ξ 0 α) ξ0 ) ξ 0ξ ) αξ 0 ) + ξ β)ξ γ + δ) βξ )γ + δξ ) ξ ξ 0 )ξ α) ξ ) ξ 0ξ ) αξ ) + α β)αγ + δ) βα)γ + δα) α ξ 0 )α ξ ) ξ 0 α) ξ α) α ). This is surprisingly complicated! Furthermore, it only applies for the case when all of ξ 0, ξ and α are real valued and distinct. evertheless, with the definitions K 0 α ) αξ 0 ) ξ 0 ) + αξ 0 ) + β ) βα + ξ 0 ), K f 8) it permits the computation in this situation for the second order case of VarGe jω, θ n ) σ [ κ K 0 e jω ξ 0 + K γe jω + δ ] e jω ξ 0 e jω ξ 9) The performance of this quantification, for the same experiment conditions as for the first order case and with ξ = 0.8, is shown in figure, and again it is seen to be essentially exact Estimate variability vs exact quantification and asymptotic approximation Sample Variability Exact Quantification Asymptotic Approximation Fig.. Variability of Fixed Denominator Estimate: True variability vs. theoretically derived approximations. Solid line is Monte Carlo estimate of true variability, dashed line is approximation 9), dash-dot line is asymptotic-based approximation ) 7. COMPARISO WITH OTHER APPROACHES The recent paper Xie and Ljung 000) has also tackled this problem of providing variance error quantification that is exact for finite model order. In that work, a completely different analysis approach is taken. Furthermore, some specific assumptions are made there that are different to those of this paper. Specifically, in Xie and Ljung 000) it is required that ) The spectral factor Ψz) is strictly all-pole ie. u t is a strictly auto-regressive time series); ) The model order chosen must be greater than or equal to the order n of the underlying dynamics plus the order ν of the spectral factor Ψz); 3) In the model structure ), all fixed poles ξ k for n < k n + ν are set to zero. With these conditions in mind, the quantification proposed in Xie and Ljung 000) is lim ) = σ n ) ξ k ν Φ u ω) e jω ξ k + α k e jω α k. 30) Clearly, it is important to reconcile this expression with that derived here. For this purpose, consider the first order example of 5 in which case, since Ψz) is also first order in that example, to satisfy the above requirements of Xie and Ljung 000) a second order model of the form Gq, θ) = θ 0 + θ q qq ξ 0 )
6 needs to be employed. In this case, substituting ξ = 0, β = 0 into the analysis of 6 leads to a quantification of ) σ κ K0 + K γe jω δ e jω ξ 0 3) with γ, δ, K 0 and K given by 7), 8). On the other hand, according to 30), the analysis of Xie and Ljung 000) gives a quantification of: ) σ ξ 0 ) e jω α + α ) e jω ξ 0 κ e jω ξ 0 3) The performance of both these quantifications for the case of ξ 0 = 0.9, α = 0.5 is shown in figure 3. Clearly, the quantification 3) seems to be essentially exact, and is also identical to that of 3) derived in Xie and Ljung 000), so that for the scenario where 3) applies, the results of this paper are consistent with it Estimate variability vs exact quantification and asymptotic approximation Sample Variability Exact Quantification Xie Ljung Quantification Fig. 3. Variability of Fixed Denominator Estimate - Second Order Case. True variability vs. theoretically derived approximations. Solid line is Monte Carlo estimate of true variability, dashed line is approximation 3), dash-dot line is Xie Ljung quantification 3) 8. COCLUSIOS The good news presented here is that in cases of rational input spectral density, it is possible to obtain variance quantifications that apply for arbitrarily low model order. The bad news is that the analytical manipulations required to compute the necessary orthonormal basis lead to very complicated expressions; this happens because some zeroes as well as all the poles are constrained in the required orthonormal spanning elements. This implies that some sort of numerical procedure may need to be developed to compute these quantities. This work is related to work by other authors, but there are important differences. In particular, the results here are derived via new reproducing kernel techniques that avoid the need for the detailed residue-based computations employed in Xie and Ljung 000). This approach allows for the inclusion of more general classes of input spectra and also avoids the requirement of model order being larger than the underlying true dynamics. However, because of this generality the quantifications here are more complicated than those of Xie and Ljung 000). 9. REFERECES Golub, G. and Charles Van Loan 989). Matrix Computations. Johns Hopkins University Press. Grenander, U. and G. Szegö 958). Toeplitz Forms and their Applications. University of California Press, Berkeley. Hannan, E. and D.F. icholls 977). The estimation of the prediction error variance. Journal of the American Statistical Association 7360, part ), Heuberger, P., P.M.J. Van den Hof and O.H. Bosgra 995). A generalized orthonormal basis for linear dynamical systems. IEEE Transactions on Automatic Control AC-403), Ljung, L. 999). System Identification: Theory for the User, nd edition). Prentice-Hall, Inc.. ew Jersey. L.Ljung 985). Asymptotic variance expressions for identified black-box transfer function models. IEEE Transactions on Automatic Control AC-309), L.Ljung and Z.D.Yuan 985). Asymptotic properties of black-box identification of transfer functions. IEEE Transactions on Automatic Control 306), inness, B. and Hakan Hjalmarsson 00a). The quantification of variance error part I: Fundamental principles. Submitted to IEEE Transactions on Automatic Control. inness, B. and Hakan Hjalmarsson 00b). The quantification of variance error part II: General models and applications. Submitted to IEEE Transactions on Automatic Control. inness, B., Håkan Hjalmarsson and Fredrik Gustafsson 999). The fundamental role of general orthonormal bases in system identification. IEEE Transactions on Automatic Control 447), P.M.J. Van den Hof, P.S.C. Heuberger and J. Bokor 995). System identification with generalized orthonormal basis functions. Automatica 3), Riesz, F. and B. Sz.-agy 955). Functional Analysis. Fredrick Ungar. ew York. Rudin, W. 966). Real and Complex Analysis. McGraw Hill. ew York. Wahlberg, B. 99). System identification using Laguerre models. IEEE Transactions on Automatic Control AC-365), Wahlberg, B. 994). System identification using Kautz models. IEEE Transactions on Automatic Control AC-396), Whittle, P. 953). Estimation and information in stationary time series. Arkiv För Matematik 3), Xie, L.-L. and Lennart Ljung 000). Asymptotic variance expressions for estimated frequency functions. Technical Report LiTH-ISY-R-77. Department of Electrical Engineering, Linköping University,To appear, IEEE Trans. Automatic Control, 00.
The Asymptotic CRLB for the Spectrum of ARMA Processes
The Asymptotic CRLB for the Spectrum of ARMA Processes Brett Ninness Abstract This paper addresses the issue of quantifying the frequency domain accuracy of ARMA spectral estimates as dictated by the Cramér
More informationCover page. : On-line damage identication using model based orthonormal. functions. Author : Raymond A. de Callafon
Cover page Title : On-line damage identication using model based orthonormal functions Author : Raymond A. de Callafon ABSTRACT In this paper, a new on-line damage identication method is proposed for monitoring
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 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 informationTHIS paper studies the input design problem in system identification.
1534 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL. 50, NO. 10, OCTOBER 2005 Input Design Via LMIs Admitting Frequency-Wise Model Specifications in Confidence Regions Henrik Jansson Håkan Hjalmarsson, Member,
More informationOPTIMAL EXPERIMENT DESIGN IN CLOSED LOOP. KTH, Signals, Sensors and Systems, S Stockholm, Sweden.
OPTIMAL EXPERIMENT DESIGN IN CLOSED LOOP Henrik Jansson Håkan Hjalmarsson KTH, Signals, Sensors and Systems, S-00 44 Stockholm, Sweden. henrik.jansson@s3.kth.se Abstract: In this contribution we extend
More informationFinite-time experiment design with multisines
Proceedings of the 7th World Congress The International Federation of Automatic Control Seoul, Korea, July 6-, 8 Finite-time experiment design with multisines X. Bombois M. Barenthin P.M.J. Van den Hof
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 informationLocal Modelling with A Priori Known Bounds Using Direct Weight Optimization
Local Modelling with A Priori Known Bounds Using Direct Weight Optimization Jacob Roll, Alexander azin, Lennart Ljung Division of Automatic Control Department of Electrical Engineering Linköpings universitet,
More informationProblem Description We remark that this classication, however, does not exhaust the possibilities to assess the model quality; see e.g., Ljung and Guo
Non-Stationary Stochastic Embedding for Transfer Function Graham C. Goodwin eegcg@cc.newcastle.edu.au Estimation Julio H. Braslavsky julio@ee.newcastle.edu.au Department of Electrical and Computer Engineering
More informationFurther Results on Model Structure Validation for Closed Loop System Identification
Advances in Wireless Communications and etworks 7; 3(5: 57-66 http://www.sciencepublishinggroup.com/j/awcn doi:.648/j.awcn.735. Further esults on Model Structure Validation for Closed Loop System Identification
More informationRobust Expansion of Uncertain Volterra Kernels into Orthonormal Series
21 American Control Conference Marriott Waterfront, Baltimore, MD, USA June 3-July 2, 21 FrA22.4 Robust Expansion of Uncertain Volterra Kernels into Orthonormal Series Alex da Rosa 1, Ricardo J. G. B.
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 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 informationIDENTIFICATION OF A TWO-INPUT SYSTEM: VARIANCE ANALYSIS
IDENTIFICATION OF A TWO-INPUT SYSTEM: VARIANCE ANALYSIS M Gevers,1 L Mišković,2 D Bonvin A Karimi Center for Systems Engineering and Applied Mechanics (CESAME) Université Catholique de Louvain B-1348 Louvain-la-Neuve,
More informationROBUST FREQUENCY DOMAIN ARMA MODELLING. Jonas Gillberg Fredrik Gustafsson Rik Pintelon
ROBUST FREQUENCY DOMAIN ARMA MODELLING Jonas Gillerg Fredrik Gustafsson Rik Pintelon Department of Electrical Engineering, Linköping University SE-581 83 Linköping, Sweden Email: gillerg@isy.liu.se, fredrik@isy.liu.se
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 informationTHE IDENTIFICATION of models on the basis of frequency
656 IEEE TRANSACTIONS ON AUTOMATIC CONTROL, VOL 43, NO 5, MAY 1998 Frequency Domain Identification Generalized Orthonormal Basis Functions Douwe K de Vries and Paul M J Van den Hof, Senior Member, IEEE
More informationCHEAPEST IDENTIFICATION EXPERIMENT WITH GUARANTEED ACCURACY IN THE PRESENCE OF UNDERMODELING
CHEAPEST IDENTIFICATION EXPERIMENT WITH GUARANTEED ACCURACY IN THE PRESENCE OF UNDERMODELING Xavier Bombois, Marion Gilson Delft Center for Systems and Control, Delft University of Technology, Mekelweg
More informationDirection of Arrival Estimation: Subspace Methods. Bhaskar D Rao University of California, San Diego
Direction of Arrival Estimation: Subspace Methods Bhaskar D Rao University of California, San Diego Email: brao@ucsdedu Reference Books and Papers 1 Optimum Array Processing, H L Van Trees 2 Stoica, P,
More informationSDP APPROXIMATION OF THE HALF DELAY AND THE DESIGN OF HILBERT PAIRS. Bogdan Dumitrescu
SDP APPROXIMATION OF THE HALF DELAY AND THE DESIGN OF HILBERT PAIRS Bogdan Dumitrescu Tampere International Center for Signal Processing Tampere University of Technology P.O.Box 553, 3311 Tampere, FINLAND
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 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 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 informationDefinition 1. A set V is a vector space over the scalar field F {R, C} iff. there are two operations defined on V, called vector addition
6 Vector Spaces with Inned Product Basis and Dimension Section Objective(s): Vector Spaces and Subspaces Linear (In)dependence Basis and Dimension Inner Product 6 Vector Spaces and Subspaces Definition
More informationPart 1a: Inner product, Orthogonality, Vector/Matrix norm
Part 1a: Inner product, Orthogonality, Vector/Matrix norm September 19, 2018 Numerical Linear Algebra Part 1a September 19, 2018 1 / 16 1. Inner product on a linear space V over the number field F A map,
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 informationNONLINEAR STRUCTURE IDENTIFICATION WITH LINEAR LEAST SQUARES AND ANOVA. Ingela Lind
OLIEAR STRUCTURE IDETIFICATIO WITH LIEAR LEAST SQUARES AD AOVA Ingela Lind Division of Automatic Control, Department of Electrical Engineering, Linköping University, SE-58 83 Linköping, Sweden E-mail:
More informationReview of similarity transformation and Singular Value Decomposition
Review of similarity transformation and Singular Value Decomposition Nasser M Abbasi Applied Mathematics Department, California State University, Fullerton July 8 7 page compiled on June 9, 5 at 9:5pm
More informationA NOVEL METHOD TO DERIVE EXPLICIT KLT KERNEL FOR AR(1) PROCESS. Mustafa U. Torun and Ali N. Akansu
A NOVEL METHOD TO DERIVE EXPLICIT KLT KERNEL FOR AR() PROCESS Mustafa U. Torun and Ali N. Akansu Department of Electrical and Computer Engineering New Jersey Institute of Technology University Heights,
More informationA general transform theory of rational orthonormal basis function expansions
A general transform theory of rational orthonormal basis function expansions homas J. de Hoog Peter S.C. Heuberger Paul M.J. Van den Hof Signals, Systems and Control Group, Department of Applied Physics
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 informationf-domain expression for the limit model Combine: 5.12 Approximate Modelling What can be said about H(q, θ) G(q, θ ) H(q, θ ) with
5.2 Approximate Modelling What can be said about if S / M, and even G / G? G(q, ) H(q, ) f-domain expression for the limit model Combine: with ε(t, ) =H(q, ) [y(t) G(q, )u(t)] y(t) =G (q)u(t) v(t) We know
More informationElementary linear algebra
Chapter 1 Elementary linear algebra 1.1 Vector spaces Vector spaces owe their importance to the fact that so many models arising in the solutions of specific problems turn out to be vector spaces. The
More informationMATH Linear Algebra
MATH 304 - Linear Algebra In the previous note we learned an important algorithm to produce orthogonal sequences of vectors called the Gramm-Schmidt orthogonalization process. Gramm-Schmidt orthogonalization
More informationGENERALISED FOURIER AND TOEPLITZ RESULTS FOR RATIONAL ORTHONORMAL BASES
GENERALSED FOURER AND TOEPLTZ RESULTS FOR RATONAL ORTHONORMAL BASES BRETT NNNESS*, HaKAN HJALMARSSONt, AND FREDRK GUSTAFSSONt Abstract. This paper provides a generalisation of certain classical Fourier
More informationA NEW APPROACH TO MIXED H 2 /H OPTIMAL PI/PID CONTROLLER DESIGN
Copyright 2002 IFAC 15th Triennial World Congress, Barcelona, Spain A NEW APPROACH TO MIXED H 2 /H OPTIMAL PI/PID CONTROLLER DESIGN Chyi Hwang,1 Chun-Yen Hsiao Department of Chemical Engineering National
More informationCOMPLEX CONSTRAINED CRB AND ITS APPLICATION TO SEMI-BLIND MIMO AND OFDM CHANNEL ESTIMATION. Aditya K. Jagannatham and Bhaskar D.
COMPLEX CONSTRAINED CRB AND ITS APPLICATION TO SEMI-BLIND MIMO AND OFDM CHANNEL ESTIMATION Aditya K Jagannatham and Bhaskar D Rao University of California, SanDiego 9500 Gilman Drive, La Jolla, CA 92093-0407
More informationDegenerate Perturbation Theory. 1 General framework and strategy
Physics G6037 Professor Christ 12/22/2015 Degenerate Perturbation Theory The treatment of degenerate perturbation theory presented in class is written out here in detail. The appendix presents the underlying
More informationPublished in: Proceedings of the 16th IFAC Symposium on System Identification (SYSID 2012), July , Brussels, Belgium
Dynamic network structure identification with prediction error methods - basic examples Dankers, AG; Van den Hof, PMJ; Heuberger, PSC; Bombois, X Published in: Proceedings of the 16th IFAC Symposium on
More informationBCT Lecture 3. Lukas Vacha.
BCT Lecture 3 Lukas Vacha vachal@utia.cas.cz Stationarity and Unit Root Testing Why do we need to test for Non-Stationarity? The stationarity or otherwise of a series can strongly influence its behaviour
More informationModel Error Concepts in Identification for Control
Model Error Concepts in Identification for Control Wolfgang Reinelt, Andrea Garulli, Lennart Ljung, Julio H. Braslavsky and Antonio Vicino Department of Electrical Engineering Linköping University, 581
More informationBOUNDARY CONDITIONS FOR SYMMETRIC BANDED TOEPLITZ MATRICES: AN APPLICATION TO TIME SERIES ANALYSIS
BOUNDARY CONDITIONS FOR SYMMETRIC BANDED TOEPLITZ MATRICES: AN APPLICATION TO TIME SERIES ANALYSIS Alessandra Luati Dip. Scienze Statistiche University of Bologna Tommaso Proietti S.E.F. e ME. Q. University
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 informationCITY UNIVERSITY LONDON. BEng (Hons) in Electrical and Electronic Engineering PART 2 EXAMINATION. ENGINEERING MATHEMATICS 2 (resit) EX2003
No: CITY UNIVERSITY LONDON BEng (Hons) in Electrical and Electronic Engineering PART 2 EXAMINATION ENGINEERING MATHEMATICS 2 (resit) EX2003 Date: August 2004 Time: 3 hours Attempt Five out of EIGHT questions
More informationDevelopment of Nonlinear Black Box Models using Orthonormal Basis Filters: A Review*
Development of Nonlinear Black Box Models using Orthonormal Basis Filters: A Review* Sachin C. atwardhan Abstract Over the last two decades, there has been a growing interest in the use of orthonormal
More informationLeast-Squares Performance of Analog Product Codes
Copyright 004 IEEE Published in the Proceedings of the Asilomar Conference on Signals, Systems and Computers, 7-0 ovember 004, Pacific Grove, California, USA Least-Squares Performance of Analog Product
More informationADAPTIVE FILTER THEORY
ADAPTIVE FILTER THEORY Fourth Edition Simon Haykin Communications Research Laboratory McMaster University Hamilton, Ontario, Canada Front ice Hall PRENTICE HALL Upper Saddle River, New Jersey 07458 Preface
More informationHANKEL-NORM BASED INTERACTION MEASURE FOR INPUT-OUTPUT PAIRING
Copyright 2002 IFAC 15th Triennial World Congress, Barcelona, Spain HANKEL-NORM BASED INTERACTION MEASURE FOR INPUT-OUTPUT PAIRING Björn Wittenmark Department of Automatic Control Lund Institute of Technology
More informationPolicy Gradient Reinforcement Learning for Robotics
Policy Gradient Reinforcement Learning for Robotics Michael C. Koval mkoval@cs.rutgers.edu Michael L. Littman mlittman@cs.rutgers.edu May 9, 211 1 Introduction Learning in an environment with a continuous
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationQuantification of frequency domain error bounds with guaranteed confidence level in Prediction Error Identification
Quantification of frequency domain error bounds with guaranteed confidence level in Prediction Error Identification X. Bombois (1), B.D.O. Anderson (2), M. Gevers (3) (1) Delft Center for Systems and Control,
More informationWWW: WWW:
Comparing Different Approaches to Model Error Modeling in Robust Identification Wolfgang Reinelt #, Andrea Garulli and Lennart Ljung # # Department of Electrical Engineering Linköping University, 581 83
More informationWe are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists. International authors and editors
We are IntechOpen, the world s leading publisher of Open Access books Built by scientists, for scientists 3,800 116,000 120M Open access books available International authors and editors Downloads Our
More informationDiscrete Simulation of Power Law Noise
Discrete Simulation of Power Law Noise Neil Ashby 1,2 1 University of Colorado, Boulder, CO 80309-0390 USA 2 National Institute of Standards and Technology, Boulder, CO 80305 USA ashby@boulder.nist.gov
More informationAdvanced Digital Signal Processing -Introduction
Advanced Digital Signal Processing -Introduction LECTURE-2 1 AP9211- ADVANCED DIGITAL SIGNAL PROCESSING UNIT I DISCRETE RANDOM SIGNAL PROCESSING Discrete Random Processes- Ensemble Averages, Stationary
More information,, rectilinear,, spherical,, cylindrical. (6.1)
Lecture 6 Review of Vectors Physics in more than one dimension (See Chapter 3 in Boas, but we try to take a more general approach and in a slightly different order) Recall that in the previous two lectures
More informationGaussian Processes. 1. Basic Notions
Gaussian Processes 1. Basic Notions Let T be a set, and X : {X } T a stochastic process, defined on a suitable probability space (Ω P), that is indexed by T. Definition 1.1. We say that X is a Gaussian
More informationSUMMARY OF MATH 1600
SUMMARY OF MATH 1600 Note: The following list is intended as a study guide for the final exam. It is a continuation of the study guide for the midterm. It does not claim to be a comprehensive list. You
More informationUNIFORMLY MOST POWERFUL CYCLIC PERMUTATION INVARIANT DETECTION FOR DISCRETE-TIME SIGNALS
UNIFORMLY MOST POWERFUL CYCLIC PERMUTATION INVARIANT DETECTION FOR DISCRETE-TIME SIGNALS F. C. Nicolls and G. de Jager Department of Electrical Engineering, University of Cape Town Rondebosch 77, South
More informationBAYESIAN INFERENCE FOR MULTI-LINE SPECTRA IN LINEAR SENSOR ARRAY
BAYESIA IFERECE FOR MULTI-LIE SPECTRA I LIEAR SESOR ARRAY Viet Hung Tran Wenwu Wang Yuhui Luo Jonathon Chambers CVSSP, University of Surrey, UK ational Physical Laboratory (PL), UK ewcastle University,
More informationCLOSED-FORM FREQUENCY ESTIMATION USING SECOND-ORDER NOTCH FILTERS. S.M. Savaresi, S. Bittanti, H.C. So*
CLOSED-FORM FREQUECY ESTIMATIO USIG SECOD-ORDER OTCH FILTERS S.M. Savaresi S. Bittanti H.C. So* Dipartimento di Elettronica e Informazione Politecnico di Milano Piazza L. da Vinci Milano ITALY. * Department
More informationRobustness in Experiment Design
IEEE TRANSACTIONS ON AUTOMATIC CONTROL Robustness in Experiment Design Cristian R. Rojas, Juan C. Agüero, Member, IEEE, James S. Welsh, Member, IEEE, Graham C. Goodwin, Fellow, IEEE and Arie Feuer, Fellow,
More informationDETECTION theory deals primarily with techniques for
ADVANCED SIGNAL PROCESSING SE Optimum Detection of Deterministic and Random Signals Stefan Tertinek Graz University of Technology turtle@sbox.tugraz.at Abstract This paper introduces various methods for
More informationMath 350 Fall 2011 Notes about inner product spaces. In this notes we state and prove some important properties of inner product spaces.
Math 350 Fall 2011 Notes about inner product spaces In this notes we state and prove some important properties of inner product spaces. First, recall the dot product on R n : if x, y R n, say x = (x 1,...,
More informationCHAPTER VIII HILBERT SPACES
CHAPTER VIII HILBERT SPACES DEFINITION Let X and Y be two complex vector spaces. A map T : X Y is called a conjugate-linear transformation if it is a reallinear transformation from X into Y, and if T (λx)
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 informationLinear Algebra: Matrix Eigenvalue Problems
CHAPTER8 Linear Algebra: Matrix Eigenvalue Problems Chapter 8 p1 A matrix eigenvalue problem considers the vector equation (1) Ax = λx. 8.0 Linear Algebra: Matrix Eigenvalue Problems Here A is a given
More informationFinite-dimensional spaces. C n is the space of n-tuples x = (x 1,..., x n ) of complex numbers. It is a Hilbert space with the inner product
Chapter 4 Hilbert Spaces 4.1 Inner Product Spaces Inner Product Space. A complex vector space E is called an inner product space (or a pre-hilbert space, or a unitary space) if there is a mapping (, )
More informationDiscrete-Time H Gaussian Filter
Proceedings of the 17th World Congress The International Federation of Automatic Control Discrete-Time H Gaussian Filter Ali Tahmasebi and Xiang Chen Department of Electrical and Computer Engineering,
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 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 informationBayes spaces: use of improper priors and distances between densities
Bayes spaces: use of improper priors and distances between densities J. J. Egozcue 1, V. Pawlowsky-Glahn 2, R. Tolosana-Delgado 1, M. I. Ortego 1 and G. van den Boogaart 3 1 Universidad Politécnica de
More informationMathematical Methods wk 1: Vectors
Mathematical Methods wk : Vectors John Magorrian, magog@thphysoxacuk These are work-in-progress notes for the second-year course on mathematical methods The most up-to-date version is available from http://www-thphysphysicsoxacuk/people/johnmagorrian/mm
More informationMathematical Methods wk 1: Vectors
Mathematical Methods wk : Vectors John Magorrian, magog@thphysoxacuk These are work-in-progress notes for the second-year course on mathematical methods The most up-to-date version is available from http://www-thphysphysicsoxacuk/people/johnmagorrian/mm
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 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 informationThe Laplace Transform
The Laplace Transform Generalizing the Fourier Transform The CTFT expresses a time-domain signal as a linear combination of complex sinusoids of the form e jωt. In the generalization of the CTFT to the
More informationVectors To begin, let us describe an element of the state space as a point with numerical coordinates, that is x 1. x 2. x =
Linear Algebra Review Vectors To begin, let us describe an element of the state space as a point with numerical coordinates, that is x 1 x x = 2. x n Vectors of up to three dimensions are easy to diagram.
More informationMATH 205C: STATIONARY PHASE LEMMA
MATH 205C: STATIONARY PHASE LEMMA For ω, consider an integral of the form I(ω) = e iωf(x) u(x) dx, where u Cc (R n ) complex valued, with support in a compact set K, and f C (R n ) real valued. Thus, I(ω)
More informationLinear Algebra. Min Yan
Linear Algebra Min Yan January 2, 2018 2 Contents 1 Vector Space 7 1.1 Definition................................. 7 1.1.1 Axioms of Vector Space..................... 7 1.1.2 Consequence of Axiom......................
More informationCOMPLEX HERMITE POLYNOMIALS: FROM THE SEMI-CIRCULAR LAW TO THE CIRCULAR LAW
Serials Publications www.serialspublications.com OMPLEX HERMITE POLYOMIALS: FROM THE SEMI-IRULAR LAW TO THE IRULAR LAW MIHEL LEDOUX Abstract. We study asymptotics of orthogonal polynomial measures of the
More informationTRUNCATED TOEPLITZ OPERATORS ON FINITE DIMENSIONAL SPACES
TRUNCATED TOEPLITZ OPERATORS ON FINITE DIMENSIONAL SPACES JOSEPH A. CIMA, WILLIAM T. ROSS, AND WARREN R. WOGEN Abstract. In this paper, we study the matrix representations of compressions of Toeplitz operators
More informationAsymptotic Geometry. of polynomials. Steve Zelditch Department of Mathematics Johns Hopkins University. Joint Work with Pavel Bleher Bernie Shiffman
Asymptotic Geometry of polynomials Steve Zelditch Department of Mathematics Johns Hopkins University Joint Work with Pavel Bleher Bernie Shiffman 1 Statistical algebraic geometry We are interested in asymptotic
More informationBLIND SEPARATION OF TEMPORALLY CORRELATED SOURCES USING A QUASI MAXIMUM LIKELIHOOD APPROACH
BLID SEPARATIO OF TEMPORALLY CORRELATED SOURCES USIG A QUASI MAXIMUM LIKELIHOOD APPROACH Shahram HOSSEII, Christian JUTTE Laboratoire des Images et des Signaux (LIS, Avenue Félix Viallet Grenoble, France.
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 informationJoint Direction-of-Arrival and Order Estimation in Compressed Sensing using Angles between Subspaces
Aalborg Universitet Joint Direction-of-Arrival and Order Estimation in Compressed Sensing using Angles between Subspaces Christensen, Mads Græsbøll; Nielsen, Jesper Kjær Published in: I E E E / S P Workshop
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 informationThe Laplace Transform
The Laplace Transform Syllabus ECE 316, Spring 2015 Final Grades Homework (6 problems per week): 25% Exams (midterm and final): 50% (25:25) Random Quiz: 25% Textbook M. Roberts, Signals and Systems, 2nd
More informationOn the Behavior of Information Theoretic Criteria for Model Order Selection
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 49, NO. 8, AUGUST 2001 1689 On the Behavior of Information Theoretic Criteria for Model Order Selection Athanasios P. Liavas, Member, IEEE, and Phillip A. Regalia,
More informationFinding eigenvalues for matrices acting on subspaces
Finding eigenvalues for matrices acting on subspaces Jakeniah Christiansen Department of Mathematics and Statistics Calvin College Grand Rapids, MI 49546 Faculty advisor: Prof Todd Kapitula Department
More informationInner product spaces. Layers of structure:
Inner product spaces Layers of structure: vector space normed linear space inner product space The abstract definition of an inner product, which we will see very shortly, is simple (and by itself is pretty
More informationChance Constrained Input Design
5th IEEE Conference on Decision and Control and European Control Conference (CDC-ECC) Orlando, FL, USA, December -5, Chance Constrained Input Design Cristian R. Rojas, Dimitrios Katselis, Håkan Hjalmarsson,
More informationMULTIPLE-CHANNEL DETECTION IN ACTIVE SENSING. Kaitlyn Beaudet and Douglas Cochran
MULTIPLE-CHANNEL DETECTION IN ACTIVE SENSING Kaitlyn Beaudet and Douglas Cochran School of Electrical, Computer and Energy Engineering Arizona State University, Tempe AZ 85287-576 USA ABSTRACT The problem
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 informationHILBERT SPACES AND THE RADON-NIKODYM THEOREM. where the bar in the first equation denotes complex conjugation. In either case, for any x V define
HILBERT SPACES AND THE RADON-NIKODYM THEOREM STEVEN P. LALLEY 1. DEFINITIONS Definition 1. A real inner product space is a real vector space V together with a symmetric, bilinear, positive-definite mapping,
More informationVector spaces. DS-GA 1013 / MATH-GA 2824 Optimization-based Data Analysis.
Vector spaces DS-GA 1013 / MATH-GA 2824 Optimization-based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_fall17/index.html Carlos Fernandez-Granda Vector space Consists of: A set V A scalar
More informationBernstein-Szegö Inequalities in Reproducing Kernel Hilbert Spaces ABSTRACT 1. INTRODUCTION
Malaysian Journal of Mathematical Sciences 6(2): 25-36 (202) Bernstein-Szegö Inequalities in Reproducing Kernel Hilbert Spaces Noli N. Reyes and Rosalio G. Artes Institute of Mathematics, University of
More informationMathematics Department Stanford University Math 61CM/DM Inner products
Mathematics Department Stanford University Math 61CM/DM Inner products Recall the definition of an inner product space; see Appendix A.8 of the textbook. Definition 1 An inner product space V is a vector
More informationCommutants of Finite Blaschke Product. Multiplication Operators on Hilbert Spaces of Analytic Functions
Commutants of Finite Blaschke Product Multiplication Operators on Hilbert Spaces of Analytic Functions Carl C. Cowen IUPUI (Indiana University Purdue University Indianapolis) Universidad de Zaragoza, 5
More information