Model Reduction using a Frequency-Limited H 2 -Cost
|
|
- Tyrone Williamson
- 5 years ago
- Views:
Transcription
1 Technical report from Automatic Control at Linköpings universitet Model Reduction using a Frequency-Limited H 2 -Cost Daniel Petersson, Johan Löfberg Division of Automatic Control petersson@isy.liu.se, johanl@isy.liu.se 7th March 212 Report no.: LiTH-ISY-R-345 Submitted to the 51st IEEE Conference on Decision and Control, Maui, Hawaii Address: Department of Electrical Engineering Linköpings universitet SE Linköping, Sweden WWW: AUTOMATIC CONTROL REGLERTEKNIK LINKÖPINGS UNIVERSITET Technical reports from the Automatic Control group in Linköping are available from
2 Abstract We propose a method for model reduction on a given frequency range, without having to specify input and output filter weights. The method uses a nonlinear optimization approach where we formulate a H 2 like cost function which only takes the given frequency range into account. We derive a gradient of the proposed cost function which enables us to use off-the-shelf optimization software. Keywords: Optimization, Model Reduction
3 Model Reduction using a Frequency-Limited H 2 -Cost Daniel Petersson and Johan Löfberg Abstract We propose a method for model reduction on a given frequency range, without having to specify input and output filter weights. The method uses a nonlinear optimization approach where we formulate a H 2 like cost function which only takes the given frequency range into account. We derive a gradient of the proposed cost function which enables us to use off-the-shelf optimization software. 1 Introduction Given a linear time-invariant (lti dynamical model, ẋ(t = Ax(t + Bu(t, y(t = Cx(t + Du(t where A R n n, B R n m, C R p n and D R p m, the model reduction problem in this paper is to find a reduced order model ẋ r (t = A r x r (t + B r u(t, y r (t = C r x r (t + D r u(t, with A r R nr nr, B r R nr m, C r R p nr and D r R p m with n r < n, where this reduced order model describes the original model well in some metric. In this paper we are interested in a reduced model that describes the model well on a given frequency range. This is motivated by situations where the given model is valid only for a certain frequency range, for example as in Varga et al. [212] where models coming from aerodynamical and structural mechanics computations describing a flexible structure are only valid up to a certain frequency. For a review of model reduction approaches, both ordinary and frequencyweighted, see e.g. Gugercin and Antoulas [24] and Ghafoor and Sreeram [28]. Some of the most commonly used frequency-weighted methods, according to Gugercin and Antoulas [24], are Wang et al. [1999], Enns [1984] and D. Petersson and J. Löfberg are with the Division of Automatic Control, Department of Electrical Engineering, Linköpings universitet, SE Sweden {petersson@isy.liu.se,johanl@isy.liu.se} 1
4 Lin and Chiu [1992], which all use different balanced truncation approaches. In many of the frequency-weighted methods one need to specify input and output filter weights. In Gawronski and Juang [199] they introduce a method which does not need these weighting functions, by introducing frequency-limited Gramians. This method can be interpreted as using ideal low-, band- or highpass filters as weights. However, this method has the drawback of not always producing stable models. One approach to remedy this has been presented in Gugercin and Antoulas [24], where they introduce a modification of the method in Gawronski and Juang [199], and derive H bound for the error. Sahlan et al. [212] presents a modification of the method from Gawronski and Juang [199], however this method is only applicable to siso models. In this paper we propose a method, based on nonlinear optimization, that uses the frequency-limited Gramians introduced in Gawronski and Juang [199], and does not have the drawback of sometimes finding unstable models. We use these Gramians to construct a frequency-limited H 2 -norm which describes the cost function of the optimization problem. We derive a gradient of the proposed cost function which enables us to use off-the-shelve optimization software to solve the problem efficiently. 2 Frequency-Limited Model Reduction The method proposed in this paper is a model reduction method that, given a model G, finds a reduced order model, Ĝ, that is a good approximation on a given frequency interval, e.g. [, ω]. The objective is to minimize the error between the given model and the sought reduced model in a frequency-limited H 2 -norm, using the frequency-limited Gramians introduced in Gawronski and Juang [199]. We formulate the optimization problem ( where minimize Ĝ G Ĝ 2 = minimize E 2 H, (1 H 2,ω 2,ω Ĝ E 2 H = 1 ω 2,ω E(iν E(iνdν. (2 2π ω Assume that the system E is stable and described by E : { ẋ(t = AE x(t + B E u(t y(t = C E x(t + D E u(t (3 which will be denoted as [ AE B E : E C E D E ]. (4 Assuming G and Ĝ are represented as [ ] [ A B Â G :, C D Ĝ : ˆB Ĉ ˆD ], (5 2
5 then the error system can be realized in state space form as ( ( [ ] A AE B E : E = Â BˆB C E D E ( C Ĉ. (6 D ˆD This realization of the error system will later prove beneficial when rewriting the optimization problem. Throughout the paper we will assume that the given model is strictly proper, i.e., D =, or assume that D = ˆD. We will also assume that the given model is stable. To be able to calculate the frequency-limited H 2 -norm we need expressions for the frequency-limited Gramians. We start by defining the standard Gramians, see Zhou et al. [1996], in time and frequency domain, as P E = Q E = e A Eτ B E B T Ee AT E τ dτ = 1 2π e AT E τ C T EC E e A Eτ dτ = 1 2π H(νB E B T EH (νdν, H (νc T EC E H(νdν, (7a (7b where H(ω = (iωi A E 1 and H (ω denotes the conjugate transpose of H(ω. The controllability and observability Gramians satisfies, respectively, the Lyapunov equations A E P E + P E A T E + B E B T E =, A T EQ E + Q E A E + C T EC E =. (8a (8b The H 2 -norm of E can be expressed as E 2 H 2 = tr = 1 2π = tr = 1 2π C E e A Eτ B E B T Ee AT E τ C T Edτ C E H(νB E B T EH (νc T Edν = tr C E P E C T E B T Ee AT E τ C T EC E e A Eτ B E dτ B T EH (νc T EC E H(νB E dν = tr B T EQ E B E. (9a (9b (9c (9d Now we want to narrow the frequency band, from (, to ( ω, ω where ω <. To do this we define the frequency-limited Gramians, see Gawronski and Juang [199], as P E,ω = 1 ω H(νB E B T 2π EH (νdν, (1a ω Q E,ω = 1 ω H (νc T 2π EC E H(νdν, (1b ω where these Gramians satisfy the Lyapunov equations, see Gawronski and Juang [199], A E P E,ω + P E,ω A T E + S E,ω B E B T E + B E B T ES E,ω =, A T EQ E,ω + Q E,ω A E + S E,ωC T EC E + C T EC E S E,ω =, (11a (11b 3
6 with S E,ω = i 2π ln ( (A E + iωi(a E iωi 1 = i 2π ln f(a E. (12 Using the frequency-limited Gramians we can express the cost function of the optimization problem (1, i.e., the frequency-limited H 2 -norm for E, as E 2 H 2,ω = tr C EP E,ω C T E (13a = tr B T EQ E,ω B E. (13b Now we want to rewrite the cost function (13 to a more computationally tractable form. This is done by using the realization given in (6 and by partitioning the Gramians P E,ω and Q E,ω as ( ( Pω X ω Qω Y ω P E,ω =, Q ˆP E,ω =, (14 ω ˆQ ω X T ω Y T ω and S E,ω as ( Sω S E,ω =. (15 Ŝ ω P ω, Q ω, ˆP ω, ˆQ ω, X ω and Y ω satisfy, due to (11, the Sylvester and Lyapunov equations AP ω + P ω A T + S ω BB T + BB T S ω =, AX ω + X ω Â T + S ω B ˆB T + B ˆB T Ŝ ω =, Â ˆP ω + ˆP ω Â T + Ŝω ˆB ˆB T + ˆB ˆB T Ŝ ω =, A T Q ω + Q ω A + S ωc T C + C T CS ω =, A T Y ω + Y ω Â S ωc T Ĉ C T ĈŜω =, Â T ˆQω + ˆQ ω Â + Ŝ ωĉt Ĉ + ĈT ĈŜω =. (16a (16b (16c (16d (16e (16f Note that P ω and Q ω satisfy the Lyapunov equations for the frequency-limited controllability and observability Gramians for the given model, and ˆP ω and ˆQ ω satisfy the Lyapunov equations for the frequency-limited controllability and observability Gramians for the sought model. With the partitioning of P E,ω and Q E,ω it is possible to rewrite (13 in two alternative forms E 2 H = tr 2,ω (B T Q ω B + 2B T Y ω ˆB + ˆBT ˆQω ˆB, (17a E 2 H 2,ω = tr (CP ω C T 2CX ω Ĉ T + Ĉ ˆP ω Ĉ T. (17b Remark 1. We note that when solving the optimization problem (1 with the cost function (17 we do not need the term B T Q ω B or CP ω C T in the cost function as they do not depend on the optimization variables, Â, ˆB and Ĉ. 2.1 Gradient of the Cost Function An appealing feature of the proposed nonlinear optimization approach to solve the problem is that the equations (17 are differentiable in the system matrices, 4
7 Â, ˆB and Ĉ. In addition, the closed form expression obtained when differentiating the cost function is expressed in the given data (A, B and C, the optimization variables (Â, ˆB and Ĉ and solutions to equations in (16. To show this we start by differentiating with respect to ˆB and Ĉ. First we note that neither Q ω, Y ω nor ˆQ ω in equation (17a depends on ˆB which means that the equation is quadratic in ˆB. Analogous observations can be made with equation (17b and the variable Ĉ. Hence, the derivative of the cost function with respect ˆB and Ĉ becomes E 2 H 2,ω ˆB E 2 H 2,ω Ĉ ( = 2 ˆQω ˆB + Y T ω B, (18a = 2 (Ĉ ˆPω CX ω. (18b For the more complicated case of differentiating with respect to  we observe that ˆQ ω and Y ω do depend on Â, see the equations in (16. The calculations of this part of the gradient are lengthy and can be found in the appendix. The complete gradient becomes where W = i π E 2 H 2,ω  E 2 H 2,ω ˆB E 2 H 2,ω Ĉ ( =2 YωX T + ˆQ ω ˆP + W, (19a ( =2 ˆQω ˆB + Y T ω B, (19b =2 (Ĉ ˆPω CX ω, (19c T [ ( (I f(â L f(â, ĈT Ĉ ˆP T CX] ĈT (  iωi T, (2 with the function L(, being the Frechét derivative of the matrix logarithm, see Higham [28]. Remark 2. By using nonlinear optimization and the fact that the derivation of the gradient is done element-wise, it is possible to impose structure in the system matrices of the reduced model, e.g., to impose a block structure in the Â-matrix. Remark 3. By using addition/subtraction of two or more different frequencylimited Gramians it is possible to focus on one or more arbitrary frequency ranges, e.g., you can construct the frequency-limited controllability Gramian, P Ω, for the interval ω [ω 1, ω 2 ] [ω 3, ω 4 ] as P Ω = P ω2 P ω1 + P ω4 P ω3. Remark 4. The proposed method can, analogous to what is done in Petersson [21], easily be extended to a method for identifying lpv-models over a limited frequency domain. Remark 5. By supplying the cost function and its gradient in computationally efficient forms this method can be used in any off-the-shelf Quasi-Newton solver. Remark 6. By using a stable model, e.g., a model from a Hankel reduction, as an initial point in the optimization and using a line-search we can limit the search to stable models. 5
8 3 Numerical Examples In this section three examples are used to illustrate the applicability of the method and to compare it with other methods. In the examples we will use three different methods; Truncation of Hankel singular values (will be called Hankel, the method proposed in Gawronski and Juang [199] (called Gawronski and the proposed method (called Prop. method. Both the proposed method and the Gawronski method will take the limited frequency range into account, but Hankel will not. The reason that we use these methods for comparison is to have one standard method (Hankel and one method that takes the limited frequency range into account. The Gawronski method is a representative method among frequency-weighted methods, see Gugercin and Antoulas [24], with the benefit of not having to design weighting functions, but with the drawback that it cannot guarantee that the resulting model is stable. The proposed method uses a cost function which is non-convex, which makes it important to use a good initial point. For the examples presented here we have used the model obtained by the Hankel method as an initial point for the optimization in the proposed method. Example 1 (Small illustrative example. This example addresses a small model with four states. The model is two second order models in series, one with a resonance frequency at ω = 1 and the other at ω = 3. We will try to limit the frequency rage to ω [, 1.7] to try to only capture the first model. 1 9 G = G 1 G 2 = s 2 +.2s + 1 s 2 +.3s + 9. (21 The results from the different methods can be seen in Figures 1 and 2 and Table 1. As can be seen in the result we are successful in finding a good model for the first model with both the proposed method and the Gawronski method. All the reduced models are stable. The Hankel method captures the wrong resonance mode (from our perspective and fails completely in the lower frequency regions. The proposed method and the Gawronski method returns models essentially indiscernible. Table 1: Numerical results for Example 1 G Ĝ G Ĝ H 2,ω Re λ max H2,ω G H2,ω Hankel 1.765e+ 1.6e e-3 Gawronski 9.14e e e-2 Prop. method 8.512e e e-2 Example 2 (Example 1 in Gawronski and Juang [199]. In this example we reuse Example 1 from Gawronski and Juang [199]. The model, which is a spring-damper model with three masses, has six states and we will reduce the model to three states and limit the frequency interval to ω [, 1.3]. The results from the different methods can be seen in Figures 3 and 4 and Table 2. The proposed method and the Gawronski method are also in this example successful in finding low order models that approximates the given model on the given frequency range, and all the reduced models are stable. All the methods find the 6
9 Magnitude plot for the given and the reduced models 4 2 Magnitude (db Hankel Gawronski Prop. method True Frequency (rad/s Figure 1: Magnitude plot of the given and the reduced models in example 1 for ω [, 1.7]. The dashed black vertical line denotes ω = 1.7. The Hankel method captures the wrong resonance mode (from our perspective and fails completely in the lower frequency regions. The proposed method and the Gawronski method returns models essentially indiscernible Magnitude plot for the error systems Magnitude (db Hankel Gawronski Prop. method Frequency (rad/s Figure 2: Magnitude plot of the error systems in example 1 for ω [, 1.7]. The dashed black vertical line denotes ω =
10 resonance mode, but once again the Hankel method fails in the lower frequencies. The proposed method, the Gawronski method and the true model are indiscernible up to the limit frequency ω = 1.3. Magnitude plot for the given and the reduced models 4 2 Hankel Gawronski Prop. method True Magnitude (db Frequency (rad/s Figure 3: Magnitude plot of the given and the reduced models in example 2 for ω [, 1.3]. The black vertical line denotes ω = 1.3. All the methods find the resonance mode, but once again the Hankel method fails in the lower frequencies. The proposed method, the Gawronski method and the true model are indiscernible up to the limit frequency ω = 1.3. Table 2: Numerical results for Example 2 G Ĝ G Ĝ H 2,ω Re λ max H2,ω G H2,ω Hankel 1.62e e e-3 Gawronski 9.464e e e-3 Prop. method 3.926e e-3-3.8e-3 Example 3 (Aircraft example. The model in this example is a model with 22 states that describes the longitudinal motion of an aircraft, see Varga et al. [212]. We will reduce this model to five states and limit the frequency range to ω [, 15]. The results from the different methods can be seen in Figures 5 and 6 and Table 3. The proposed method and Gawronski method yields models with good agreement up to the limit frequency. Note though that the Gawronski method gives a model with a peculiar resonance mode at ω 11 rad/s, and the model turns out to be unstable. The Hankel method returns a model structurally different from the two other methods. In all three of the above examples we observe that the proposed method finds an at least as good model as the method proposed in Gawronski and 8
11 Magnitude plot for the error systems 1 1 Hankel Gawronski Prop. method 2 Magnitude (db Frequency (rad/s Figure 4: Magnitude plot of the error systems in example 2 for ω [, 1.3]. The dashed black vertical line denotes ω = 1.3. Table 3: Numerical results for Example 3 G Ĝ G Ĝ H 2,ω Re λ max H2,ω G H2,ω Hankel 3.7e e e+1 Gawronski 5.443e e-3 4.1e+3 Prop. method 5.899e e e+ Juang [199], but with the proposed model we can guarantee that the reduced model is stable and also impose structure in the system matrices. In the first two examples it takes less than one second and in the third example about seven seconds for the proposed method to find a reduced model. 4 Conclusions In this paper we have proposed a method, based on nonlinear optimization, that uses the frequency-limited Gramians introduced in Gawronski and Juang [199], and does not have the drawback of finding unstable models. We use these Gramians to construct a frequency-limited H 2 -norm which describes the cost function of the optimization problem. We derive a gradient of the proposed cost function which enables us to use off-the-shelve optimization software to solve the problem efficiently. The derivation of the method also enables us to impose structural constraints, e.g., upper triangular A-matrix, in the system matrices. The derivation follows closely the technique in Petersson [21] and it is easy to extend the method to identifying lpv-models. We also present three examples of different sizes and characteristics to show the applicability of the method. 9
12 Magnitude plot for the given and the reduced models Hankel Gawronski Prop. method True Magnitude (db Frequency (rad/s Figure 5: Magnitude plot of the given and the reduced models in example 3 for ω [, 15]. The dashed black vertical line denotes ω = 15. The proposed method and Gawronski method yields models with good agreement up to the limit frequency. Note though that the Gawronski method gives a model with a peculiar resonance mode at ω 11 rad/s, and the model turns out to be unstable. The Hankel method returns a model structurally different from the two other methods. References Dale F. Enns. Model reduction with balanced realizations: An error bound and a frequency weighted generalization. In Proceedings of the 23rd IEEE Conference on Decision and Control, pages , Las Vegas, USA, Wodek Gawronski and Jer-Nan Juang. Model reduction in limited time and frequency intervals. International Journal of Systems Science, 21(2: , 199. Abdul Ghafoor and Victor Sreeram. A survey/review of frequency-weighted balanced model reduction techniques. Journal of Dynamic Systems, Measurement and Control, 13:614, 28. Serkan Gugercin and Athanasios C. Antoulas. A survey of model reduction by balanced truncation and some new results. International Journal of Control, 77(8: , 24. Nicholas J. Higham. Functions of Matrices: Theory and Computation. SIAM, 28. Ching-An Lin and Tai-Yih Chiu. Model reduction via frequency weighted balanced realization. Control Theory and Advanced Technology, 8: ,
13 Magnitude plot for the error systems 5 4 Hankel Gawronski Prop. method 3 Magnitude (db Frequency (rad/s Figure 6: Magnitude plot of the error systems in example 3 for ω [, 15]. The dashed black vertical line denotes ω = 15. Daniel Petersson. Nonlinear optimization approaches to H 2 -norm based LPV modelling and control. Licentiate thesis no. 1453, Department of Electrical Engineering, Linköping University, 21. Shafishuhaza Sahlan, Abdul Ghafoor, and Victor Sreeram. A new method for the model reduction technique via a limited frequency interval impulse response gramian. Mathematical and Computer Modelling, 55(3-4:134 14, 212. Andreas Varga, Anders Hansson, and Guilhem Puyou, editors. Optimization Based Clearance of Flight Control Laws. Lecture Notes in Control and Information Science. Springer, 212. G. Wang, Victor Sreeram, and W. Q. Liu. A new frequency-weighted balanced truncation method and an error bound. IEEE Transactions on Automatic Control, 44(9: , Kemin Zhou, John C. Doyle, and Keith Glover. Robust and optimal control. Prentice-Hall, Inc., ISBN A Appendix In this appendix we present the differentiation of the cost function with respect to Â. The cost function is E 2 H = tr 2,ω (B T Q ω B + 2B T Y ω ˆB + ˆBT ˆQω ˆB (22 11
14 and by looking at the equations A T Y ω + Y ω  S ωc T Ĉ C T ĈŜω =,  T ˆQω + ˆQ ω  + Ŝ ωĉt Ĉ + ĈT ĈŜω =, (23a (23b we observe that ˆQ ω and Y ω depend on  which [ we need] to keep in mind when differentiating (22 with respect to Â. Hence, becomes where Yω in (23,  T ˆQ ω E 2 H 2,ω  [ ] ( 2 E H 2,ω = tr 2 ˆBB T Y ω +  ⠈B ˆB T ˆQ ω, (24 ij ij and ˆQ ω  T YT ω depend on  via the differentiated versions of the equations + YT ω A + ÂT Y T ω Ŝ ω a ij Ĉ T C =, ij (25a + ˆQ ω  + ÂT ˆQω + ˆQ ω  + Ŝ ω a ij Ĉ T Ĉ + ĈT Ĉ Ŝω a ij =. (25b To be able to substitute ˆQ ω and Yω tractable we use the following lemma. to something more computationally Lemma 1. If M and N satisfy the Sylvester equations then tr CN = tr DM. AM + MB + C =, NA + BN + D =, Studying Lemma 1, the two factors in front of Yω and ˆQ ω in (24 and the structure of the Lyapunov/Sylvester equations in (25, ÂT +  + = and  T + A + =, brings us to the conclusion that, to do the substitution, we need to solve two additional Lyapunov/Sylvester equations, namely AX + XÂT + B ˆB T =,  ˆP + ˆPÂT + ˆB ˆB T =. (26a (26b Note that ˆP in (26b happens to be the controllability Gramian for the reduced model. Rewriting (24 using Lemma 1, (25 and (26 leads to [ ] [ 2 E ( H 2,ω = 2 tr ÂT Y  â ωx T + ˆQ ω ˆP + Ŝ ω (ĈT Ĉ ij â ˆP ] ĈT CX. ij ij (27 What remains is to rewrite the second term in (27, which includes Ŝ ω. Recall the definition of Ŝω, Ŝ ω = i ( 2π ln ( + iωi(â iωi 1 (28 12
15 and differentiate with respect to an element in Â, i.e., a ij. This yields ( ( Ŝω = i a ij 2π L f(â f(â, = i a ij 2π L  f(â, (I f(â (A iωi 1 a ij (29 where L(A, E is the Frechét derivative of the matrix logarithm, see Higham [28], with L(A, E = 1 (t(a I + I 1 E (t(a I + I 1 dt, (3 f(â =( + iωi(â iωi 1. (31 The function L(A, E can be efficiently evaluated using the algorithm by Higham [28]. By substituting (29 into (27 and using (3 with the fact that we can interchange the tr-operator and the integral we obtain [ ] 2 E H 2,ω  where W = i π = 2 tr ij = 2 tr [ ( ÂT Y â ωx T + ˆQ ω ˆP + Ŝ ω (ĈT Ĉ ij â ˆP ] ĈT CX ij ( ( ÂT Y â ωx T + ˆQ ω ˆP ij ( [ L + i ( π tr  T ( T I â f(â ij f(â, ĈT Ĉ ˆP ] T ĈT CX (  iωi T = tr ( [ ( ÂT 2 Y â ωx T + ˆQ ω ˆP + W], (32 ij T [ ( (I f(â L f(â, ĈT Ĉ ˆP T CX] ĈT (  iωi T. (33 13
16
17 Avdelning, Institution Division, Department Datum Date Division of Automatic Control Department of Electrical Engineering Språk Language Rapporttyp Report category ISBN Svenska/Swedish Licentiatavhandling ISRN Engelska/English Examensarbete C-uppsats D-uppsats Övrig rapport Serietitel och serienummer Title of series, numbering ISSN URL för elektronisk version LiTH-ISY-R-345 Titel Title Model Reduction using a Frequency-Limited H 2 -Cost Författare Author Daniel Petersson, Johan Löfberg Sammanfattning Abstract We propose a method for model reduction on a given frequency range, without having to specify input and output filter weights. The method uses a nonlinear optimization approach where we formulate a H 2 like cost function which only takes the given frequency range into account. We derive a gradient of the proposed cost function which enables us to use off-the-shelf optimization software. Nyckelord Keywords Optimization, Model Reduction
An efficient implementation of gradient and Hessian calculations of the coefficients of the characteristic polynomial of I XY
Technical report from Automatic Control at Linköpings universitet An efficient implementation of gradient and Hessian calculations of the coefficients of the characteristic polynomial of I XY Daniel Ankelhed
More informationBlock diagonalization of matrix-valued sum-of-squares programs
Block diagonalization of matrix-valued sum-of-squares programs Johan Löfberg Division of Automatic Control Department of Electrical Engineering Linköpings universitet, SE-581 83 Linköping, Sweden WWW:
More informationModel Reduction for Unstable Systems
Model Reduction for Unstable Systems Klajdi Sinani Virginia Tech klajdi@vt.edu Advisor: Serkan Gugercin October 22, 2015 (VT) SIAM October 22, 2015 1 / 26 Overview 1 Introduction 2 Interpolatory Model
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 informationProperties and approximations of some matrix variate probability density functions
Technical report from Automatic Control at Linköpings universitet Properties and approximations of some matrix variate probability density functions Karl Granström, Umut Orguner Division of Automatic Control
More informationA relaxation of the strangeness index
echnical report from Automatic Control at Linköpings universitet A relaxation of the strangeness index Henrik idefelt, orkel Glad Division of Automatic Control E-mail: tidefelt@isy.liu.se, torkel@isy.liu.se
More informationUpper and Lower Bounds of Frequency Interval Gramians for a Class of Perturbed Linear Systems Shaker, Hamid Reza
Aalborg Universitet Upper and Lower Bounds of Frequency Interval Gramians for a Class of Perturbed Linear Systems Shaker, Hamid Reza Published in: 7th IFAC Symposium on Robust Control Design DOI (link
More informationRobust finite-frequency H2 analysis
Technical report from Automatic Control at Linköpings universitet Robust finite-frequency H2 analysis Alfio Masi, Ragnar Wallin, Andrea Garulli, Anders Hansson Division of Automatic Control E-mail: masi@dii.unisi.it,
More informationLinköping University Electronic Press
Linköping University Electronic Press Report Simulation Model of a 2 Degrees of Freedom Industrial Manipulator Patrik Axelsson Series: LiTH-ISY-R, ISSN 400-3902, No. 3020 ISRN: LiTH-ISY-R-3020 Available
More informationSystem analysis of a diesel engine with VGT and EGR
System analysis of a diesel engine with VGT and EGR Master s thesis performed in Vehicular Systems by Thomas Johansson Reg nr: LiTH-ISY-EX- -5/3714- -SE 9th October 25 System analysis of a diesel engine
More informationRobust Control Of A Flexible Manipulator Arm: A Benchmark Problem
Robust Control Of A Flexible Manipulator Arm: A Benchmark Problem Stig Moberg, Jonas Öhr Division of Automatic Control Department of Electrical Engineering Linköpings universitet, SE-581 83 Linköping,
More informationRobust Multivariable Control
Lecture 2 Anders Helmersson anders.helmersson@liu.se ISY/Reglerteknik Linköpings universitet Today s topics Today s topics Norms Today s topics Norms Representation of dynamic systems Today s topics Norms
More informationOn Consistency of Closed-loop Subspace Identifictaion with Innovation Estimation
Technical report from Automatic Control at Linköpings universitet On Consistency of Closed-loop Subspace Identictaion with Innovation Estimation Weilu Lin, S Joe Qin, Lennart Ljung Division of Automatic
More informationModel reduction for linear systems by balancing
Model reduction for linear systems by balancing Bart Besselink Jan C. Willems Center for Systems and Control Johann Bernoulli Institute for Mathematics and Computer Science University of Groningen, Groningen,
More informationLane departure detection for improved road geometry estimation
Lane departure detection for improved road geometry estimation Thomas B. Schön, Andreas Eidehall, Fredrik Gustafsson Division of Automatic Control Department of Electrical Engineering Linköpings universitet,
More informationH 2 optimal model reduction - Wilson s conditions for the cross-gramian
H 2 optimal model reduction - Wilson s conditions for the cross-gramian Ha Binh Minh a, Carles Batlle b a School of Applied Mathematics and Informatics, Hanoi University of Science and Technology, Dai
More informationH 2 -optimal model reduction of MIMO systems
H 2 -optimal model reduction of MIMO systems P. Van Dooren K. A. Gallivan P.-A. Absil Abstract We consider the problem of approximating a p m rational transfer function Hs of high degree by another p m
More informationRobust Heading Estimation Indoors
Technical report from Automatic Control at Linköpings universitet Robust Heading Estimation Indoors Jonas Callmer, David Törnqvist, Fredrik Gustafsson Division of Automatic Control E-mail: callmer@isy.liu.se,
More informationFREQUENCY-WEIGHTED MODEL REDUCTION METHOD WITH ERROR BOUNDS FOR 2-D SEPARABLE DENOMINATOR DISCRETE SYSTEMS
INERNAIONAL JOURNAL OF INFORMAION AND SYSEMS SCIENCES Volume 1, Number 2, Pages 105 119 c 2005 Institute for Scientific Computing and Information FREQUENCY-WEIGHED MODEL REDUCION MEHOD WIH ERROR BOUNDS
More informationImplementation of the GIW-PHD filter
Technical reort from Automatic Control at Linöings universitet Imlementation of the GIW-PHD filter Karl Granström, Umut Orguner Division of Automatic Control E-mail: arl@isy.liu.se, umut@isy.liu.se 28th
More informationBalanced Truncation 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 2004: MODEL REDUCTION Balanced Truncation This lecture introduces balanced truncation for LTI
More informationKrylov Techniques for Model Reduction of Second-Order Systems
Krylov Techniques for Model Reduction of Second-Order Systems A Vandendorpe and P Van Dooren February 4, 2004 Abstract The purpose of this paper is to present a Krylov technique for model reduction of
More informationIterative Rational Krylov Algorithm for Unstable Dynamical Systems and Generalized Coprime Factorizations
Iterative Rational Krylov Algorithm for Unstable Dynamical Systems and Generalized Coprime Factorizations Klajdi Sinani Thesis submitted to the Faculty of the Virginia Polytechnic Institute and State University
More informationBALANCING-RELATED MODEL REDUCTION FOR DATA-SPARSE SYSTEMS
BALANCING-RELATED Peter Benner Professur Mathematik in Industrie und Technik Fakultät für Mathematik Technische Universität Chemnitz Computational Methods with Applications Harrachov, 19 25 August 2007
More informationFluid flow dynamical model approximation and control
Fluid flow dynamical model approximation and control... a case-study on an open cavity flow C. Poussot-Vassal & D. Sipp Journée conjointe GT Contrôle de Décollement & GT MOSAR Frequency response of an
More informationDepartment of Physics, Chemistry and Biology
Department of Physics, Chemistry and Biology Master s Thesis Quantum Chaos On A Curved Surface John Wärnå LiTH-IFM-A-EX-8/7-SE Department of Physics, Chemistry and Biology Linköpings universitet, SE-58
More informationProjection of state space realizations
Chapter 1 Projection of state space realizations Antoine Vandendorpe and Paul Van Dooren Department of Mathematical Engineering Université catholique de Louvain B-1348 Louvain-la-Neuve Belgium 1.0.1 Description
More informationGramians based model reduction for hybrid switched systems
Gramians based model reduction for hybrid switched systems Y. Chahlaoui Younes.Chahlaoui@manchester.ac.uk Centre for Interdisciplinary Computational and Dynamical Analysis (CICADA) School of Mathematics
More informationA Weighting Method for Approximate Nonlinear System Identification
Technical report from Automatic Control at Linköpings universitet A Weighting Method for Approximate Nonlinear System Identification Martin Enqvist Division of Automatic Control E-mail: maren@isy.liu.se
More informationMODEL REDUCTION OF SWITCHED SYSTEMS BASED ON SWITCHING GENERALIZED GRAMIANS. Received May 2011; revised September 2011
International Journal of Innovative Computing, Information and Control ICIC International c 2012 ISSN 1349-4198 Volume 8, Number 7(B), July 2012 pp. 5025 5044 MODEL REDUCTION OF SWITCHED SYSTEMS BASED
More informationSingular Value Decomposition Based Model Order Reduction Techniques
Singular Value Decomposition Based Model Order Reduction Techniques by Ahmad Jazlan Bin Haja Mohideen A thesis submitted to the School of Electrical, Electronic and Computer Engineering in partial fulfilment
More informationRaktim Bhattacharya. . AERO 632: Design of Advance Flight Control System. Norms for Signals and Systems
. AERO 632: Design of Advance Flight Control System Norms for Signals and. Raktim Bhattacharya Laboratory For Uncertainty Quantification Aerospace Engineering, Texas A&M University. Norms for Signals ...
More informationInstitutionen för systemteknik
Institutionen för systemteknik Department of Electrical Engineering Examensarbete Multiple Platform Bias Error Estimation Examensarbete utfört i Reglerteknik vid Tekniska högskolan i Linköping av Åsa Wiklund
More informationExperimental Comparison of Methods for Multivariable Frequency Response Function Estimation
Technical report from Automatic Control at Linköpings universitet Experimental Comparison of Methods for Multivariable Frequency Response Function Estimation Erik Wernholt, Stig Moberg Division of Automatic
More informationFourier Model Reduction for Large-Scale Applications in Computational Fluid Dynamics
Fourier Model Reduction for Large-Scale Applications in Computational Fluid Dynamics K. Willcox and A. Megretski A new method, Fourier model reduction (FMR), for obtaining stable, accurate, low-order models
More informationFEL3210 Multivariable Feedback Control
FEL3210 Multivariable Feedback Control Lecture 8: Youla parametrization, LMIs, Model Reduction and Summary [Ch. 11-12] Elling W. Jacobsen, Automatic Control Lab, KTH Lecture 8: Youla, LMIs, Model Reduction
More informationModel reduction of interconnected systems
Model reduction of interconnected systems A Vandendorpe and P Van Dooren 1 Introduction Large scale linear systems are often composed of subsystems that interconnect to each other Instead of reducing the
More informationModel reduction via tangential interpolation
Model reduction via tangential interpolation K. Gallivan, A. Vandendorpe and P. Van Dooren May 14, 2002 1 Introduction Although most of the theory presented in this paper holds for both continuous-time
More informationControl Systems Design
ELEC4410 Control Systems Design Lecture 13: Stability Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 13: Stability p.1/20 Outline Input-Output
More information6.241 Dynamic Systems and Control
6.241 Dynamic Systems and Control Lecture 24: H2 Synthesis Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology May 4, 2011 E. Frazzoli (MIT) Lecture 24: H 2 Synthesis May
More informationWeighted balanced realization and model reduction for nonlinear systems
Weighted balanced realization and model reduction for nonlinear systems Daisuke Tsubakino and Kenji Fujimoto Abstract In this paper a weighted balanced realization and model reduction for nonlinear systems
More informationPosition Estimation and Modeling of a Flexible Industrial Robot
Position Estimation and Modeling of a Flexible Industrial Robot Rickard Karlsson, Mikael Norrlöf, Division of Automatic Control Department of Electrical Engineering Linköpings universitet, SE-581 83 Linköping,
More informationA new approach to lane guidance systems
A new approach to lane guidance systems Andreas Eidehall, Jochen Pohl, Fredrik Gustafsson Division of Automatic Control Department of Electrical Engineering Linköpings universitet, SE-581 83 Linköping,
More informationGramians of structured systems and an error bound for structure-preserving model reduction
Gramians of structured systems and an error bound for structure-preserving model reduction D.C. Sorensen and A.C. Antoulas e-mail:{sorensen@rice.edu, aca@rice.edu} September 3, 24 Abstract In this paper
More informationLMI Based Model Order Reduction Considering the Minimum Phase Characteristic of the System
LMI Based Model Order Reduction Considering the Minimum Phase Characteristic of the System Gholamreza Khademi, Haniyeh Mohammadi, and Maryam Dehghani School of Electrical and Computer Engineering Shiraz
More informationLINEAR QUADRATIC OPTIMAL CONTROL BASED ON DYNAMIC COMPENSATION. Received October 2010; revised March 2011
International Journal of Innovative Computing, Information and Control ICIC International c 22 ISSN 349-498 Volume 8, Number 5(B), May 22 pp. 3743 3754 LINEAR QUADRATIC OPTIMAL CONTROL BASED ON DYNAMIC
More informationNon-linear Dynamic Inversion
Control of Unmanned Aerial Vehicles using Non-linear Dynamic Inversion Master s Thesis Division of Automatic Control Department of Electrical Engineering Linköping University Mia Karlsson LiTH-ISY-EX-3-
More informationNorm invariant discretization for sampled-data fault detection
Automatica 41 (25 1633 1637 www.elsevier.com/locate/automatica Technical communique Norm invariant discretization for sampled-data fault detection Iman Izadi, Tongwen Chen, Qing Zhao Department of Electrical
More information1 Continuous-time Systems
Observability Completely controllable systems can be restructured by means of state feedback to have many desirable properties. But what if the state is not available for feedback? What if only the output
More informationDiscrete and continuous dynamic systems
Discrete and continuous dynamic systems Bounded input bounded output (BIBO) and asymptotic stability Continuous and discrete time linear time-invariant systems Katalin Hangos University of Pannonia Faculty
More information3 Gramians and Balanced Realizations
3 Gramians and Balanced Realizations In this lecture, we use an optimization approach to find suitable realizations for truncation and singular perturbation of G. It turns out that the recommended realizations
More informationApproximation of the Linearized Boussinesq Equations
Approximation of the Linearized Boussinesq Equations Alan Lattimer Advisors Jeffrey Borggaard Serkan Gugercin Department of Mathematics Virginia Tech SIAM Talk Series, Virginia Tech, April 22, 2014 Alan
More informationInstitutionen för systemteknik
Institutionen för systemteknik Department of Electrical Engineering Examensarbete Estimation of Inter-cell Interference in 3G Communication Systems Examensarbete utfört i Reglerteknik vid Tekniska högskolan
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 informationModel reduction using a frequency-limited H2- cost
Model reducton usng a frequency-lmted H2- cost Danel Petersson and Johan Löfberg Lnköpng Unversty Post Prnt N.B.: When ctng ths work, cte the orgnal artcle. Orgnal Publcaton: Danel Petersson and Johan
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 informationMulti-Model Adaptive Regulation for a Family of Systems Containing Different Zero Structures
Preprints of the 19th World Congress The International Federation of Automatic Control Multi-Model Adaptive Regulation for a Family of Systems Containing Different Zero Structures Eric Peterson Harry G.
More informationCME 345: MODEL REDUCTION
CME 345: MODEL REDUCTION Balanced Truncation Charbel Farhat & David Amsallem Stanford University cfarhat@stanford.edu These slides are based on the recommended textbook: A.C. Antoulas, Approximation of
More informationA Singular Value Decomposition Based Closed Loop Stability Preserving Controller Reduction Method
A Singular Value Decomposition Based Closed Loop Stability Preserving Controller Reduction Method Kin Cheong Sou and Anders Rantzer Abstract In this paper a controller reduction method which preserves
More informationOn Feedback Linearization for Robust Tracking Control of Flexible Joint Robots
Technical report from Automatic Control at Linköpings universitet On Feedback Linearization for Robust Tracking Control of Flexible Joint Robots Stig Moberg, Sven Hanssen Division of Automatic Control
More informationLarge-scale and infinite dimensional dynamical systems approximation
Large-scale and infinite dimensional dynamical systems approximation Igor PONTES DUFF PEREIRA Doctorant 3ème année - ONERA/DCSD Directeur de thèse: Charles POUSSOT-VASSAL Co-encadrant: Cédric SEREN 1 What
More informationModule 07 Controllability and Controller Design of Dynamical LTI Systems
Module 07 Controllability and Controller Design of Dynamical LTI Systems Ahmad F. Taha EE 5143: Linear Systems and Control Email: ahmad.taha@utsa.edu Webpage: http://engineering.utsa.edu/ataha October
More informationThe model reduction algorithm proposed is based on an iterative two-step LMI scheme. The convergence of the algorithm is not analyzed but examples sho
Model Reduction from an H 1 /LMI perspective A. Helmersson Department of Electrical Engineering Linkoping University S-581 8 Linkoping, Sweden tel: +6 1 816 fax: +6 1 86 email: andersh@isy.liu.se September
More informationME 234, Lyapunov and Riccati Problems. 1. This problem is to recall some facts and formulae you already know. e Aτ BB e A τ dτ
ME 234, Lyapunov and Riccati Problems. This problem is to recall some facts and formulae you already know. (a) Let A and B be matrices of appropriate dimension. Show that (A, B) is controllable if and
More information1. Find the solution of the following uncontrolled linear system. 2 α 1 1
Appendix B Revision Problems 1. Find the solution of the following uncontrolled linear system 0 1 1 ẋ = x, x(0) =. 2 3 1 Class test, August 1998 2. Given the linear system described by 2 α 1 1 ẋ = x +
More informationMODEL REDUCTION BY A CROSS-GRAMIAN APPROACH FOR DATA-SPARSE SYSTEMS
MODEL REDUCTION BY A CROSS-GRAMIAN APPROACH FOR DATA-SPARSE SYSTEMS Ulrike Baur joint work with Peter Benner Mathematics in Industry and Technology Faculty of Mathematics Chemnitz University of Technology
More informationA Novel Integral-Based Event Triggering Control for Linear Time-Invariant Systems
53rd IEEE Conference on Decision and Control December 15-17, 2014. Los Angeles, California, USA A Novel Integral-Based Event Triggering Control for Linear Time-Invariant Systems Seyed Hossein Mousavi 1,
More informationMixed Parametric/Unstructured LFT Modelling for Robust Controller Design
2 American Control Conference on O'Farrell Street, San Francisco, CA, USA June 29 - July, 2 Mixed Parametric/Unstructured LFT Modelling for Robust Controller Design Harald Pfifer and Simon Hecker Abstract
More informationProblem Set 5 Solutions 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Problem Set 5 Solutions The problem set deals with Hankel
More informationSiMpLIfy: A Toolbox for Structured Model Reduction
SiMpLIfy: A Toolbox for Structured Model Reduction Martin Biel, Farhad Farokhi, and Henrik Sandberg ACCESS Linnaeus Center, KTH Royal Institute of Technology Department of Electrical and Electronic Engineering,
More informationarxiv: v1 [math.oc] 17 Oct 2014
SiMpLIfy: A Toolbox for Structured Model Reduction Martin Biel, Farhad Farokhi, and Henrik Sandberg arxiv:1414613v1 [mathoc] 17 Oct 214 Abstract In this paper, we present a toolbox for structured model
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 information16.31 Fall 2005 Lecture Presentation Mon 31-Oct-05 ver 1.1
16.31 Fall 2005 Lecture Presentation Mon 31-Oct-05 ver 1.1 Charles P. Coleman October 31, 2005 1 / 40 : Controllability Tests Observability Tests LEARNING OUTCOMES: Perform controllability tests Perform
More informationModel reduction of large-scale dynamical systems
Model reduction of large-scale dynamical systems Lecture III: Krylov approximation and rational interpolation Thanos Antoulas Rice University and Jacobs University email: aca@rice.edu URL: www.ece.rice.edu/
More informationROBUST PASSIVE OBSERVER-BASED CONTROL FOR A CLASS OF SINGULAR SYSTEMS
INTERNATIONAL JOURNAL OF INFORMATON AND SYSTEMS SCIENCES Volume 5 Number 3-4 Pages 480 487 c 2009 Institute for Scientific Computing and Information ROBUST PASSIVE OBSERVER-BASED CONTROL FOR A CLASS OF
More informationLinear System Theory
Linear System Theory Wonhee Kim Chapter 6: Controllability & Observability Chapter 7: Minimal Realizations May 2, 217 1 / 31 Recap State space equation Linear Algebra Solutions of LTI and LTV system Stability
More informationFrequency Domain Versus Time Domain Methods in System Identification Revisited
Technical report from Automatic Control at Linköpings universitet Frequency Domain Versus Time Domain Methods in System Identification Revisited Lennart Ljung Division of Automatic Control E-mail: ljung@isy.liu.se
More informationMarcus Pantoja da Silva 1 and Celso Pascoli Bottura 2. Abstract: Nonlinear systems with time-varying uncertainties
A NEW PROPOSAL FOR H NORM CHARACTERIZATION AND THE OPTIMAL H CONTROL OF NONLINEAR SSTEMS WITH TIME-VARING UNCERTAINTIES WITH KNOWN NORM BOUND AND EXOGENOUS DISTURBANCES Marcus Pantoja da Silva 1 and Celso
More informationarxiv: v1 [math.na] 12 Jan 2012
Interpolatory Weighted-H 2 Model Reduction Branimir Anić, Christopher A. Beattie, Serkan Gugercin, and Athanasios C. Antoulas Preprint submitted to Automatica, 2 January 22 arxiv:2.2592v [math.na] 2 Jan
More informationInfluence of Laser Radar Sensor Parameters on Range Measurement and Shape Fitting Uncertainties
Influence of Laser Radar Sensor Parameters on Range Measurement and Shape Fitting Uncertainties Christina Grönwall, Ove Steinvall, Fredrik Gustafsson, Tomas Chevalier Division of Automatic Control Department
More informationObservability. Dynamic Systems. Lecture 2 Observability. Observability, continuous time: Observability, discrete time: = h (2) (x, u, u)
Observability Dynamic Systems Lecture 2 Observability Continuous time model: Discrete time model: ẋ(t) = f (x(t), u(t)), y(t) = h(x(t), u(t)) x(t + 1) = f (x(t), u(t)), y(t) = h(x(t)) Reglerteknik, ISY,
More informationDECENTRALIZED CONTROL DESIGN USING LMI MODEL REDUCTION
Journal of ELECTRICAL ENGINEERING, VOL. 58, NO. 6, 2007, 307 312 DECENTRALIZED CONTROL DESIGN USING LMI MODEL REDUCTION Szabolcs Dorák Danica Rosinová Decentralized control design approach based on partial
More informationModern Optimal Control
Modern Optimal Control Matthew M. Peet Arizona State University Lecture 19: Stabilization via LMIs Optimization Optimization can be posed in functional form: min x F objective function : inequality constraints
More informationModel reduction of large-scale systems by least squares
Model reduction of large-scale systems by least squares Serkan Gugercin Department of Mathematics, Virginia Tech, Blacksburg, VA, USA gugercin@mathvtedu Athanasios C Antoulas Department of Electrical and
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 informationKrylov-Subspace Based Model Reduction of Nonlinear Circuit Models Using Bilinear and Quadratic-Linear Approximations
Krylov-Subspace Based Model Reduction of Nonlinear Circuit Models Using Bilinear and Quadratic-Linear Approximations Peter Benner and Tobias Breiten Abstract We discuss Krylov-subspace based model reduction
More informationStability preserving post-processing methods applied in the Loewner framework
Ion Victor Gosea and Athanasios C. Antoulas (Jacobs University Bremen and Rice University May 11, 2016 Houston 1 / 20 Stability preserving post-processing methods applied in the Loewner framework Ion Victor
More informationInstitutionen för systemteknik
Institutionen för systemteknik Department of Electrical Engineering Examensarbete Load flow control and optimization of Banverket s 132 kv 16 2/3 Hz high voltage grid Examensarbete utfört i Reglerteknik
More informationBalancing of Lossless and Passive Systems
Balancing of Lossless and Passive Systems Arjan van der Schaft Abstract Different balancing techniques are applied to lossless nonlinear systems, with open-loop balancing applied to their scattering representation.
More informationGain Scheduled Missile Control Using Robust Loop Shaping Examensarbete utfört i Reglerteknik vid Tekniska Högskolan i Linköping av Henrik Johansson Reg nr: LiTH-ISY-EX-3291-22 Gain Scheduled Missile Control
More informationL 2 -optimal model reduction for unstable systems using iterative rational
L 2 -optimal model reduction for unstable systems using iterative rational Krylov algorithm Caleb agruder, Christopher Beattie, Serkan Gugercin Abstract Unstable dynamical systems can be viewed from a
More informationHankel Optimal Model Reduction 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.242, Fall 2004: MODEL REDUCTION Hankel Optimal Model Reduction 1 This lecture covers both the theory and
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science : Dynamic Systems Spring 2011
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science 6.4: Dynamic Systems Spring Homework Solutions Exercise 3. a) We are given the single input LTI system: [
More informationInstitutionen för datavetenskap Department of Computer and Information Science
Institutionen för datavetenskap Department of Computer and Information Science Final thesis Conditional Planning for Troubleshooting and Repair in a Partially Observable Environment by Håkan Warnquist
More informationInstitutionen för systemteknik
Institutionen för systemteknik Department of Electrical Engineering Examensarbete Modeling and Identification of the Gantry-Tau Parallel Kinematic Machine Examensarbete utfört i Reglerteknik vid Tekniska
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 informationThe Important State Coordinates of a Nonlinear System
The Important State Coordinates of a Nonlinear System Arthur J. Krener 1 University of California, Davis, CA and Naval Postgraduate School, Monterey, CA ajkrener@ucdavis.edu Summary. We offer an alternative
More informationECEN 420 LINEAR CONTROL SYSTEMS. Lecture 6 Mathematical Representation of Physical Systems II 1/67
1/67 ECEN 420 LINEAR CONTROL SYSTEMS Lecture 6 Mathematical Representation of Physical Systems II State Variable Models for Dynamic Systems u 1 u 2 u ṙ. Internal Variables x 1, x 2 x n y 1 y 2. y m Figure
More informationBUMPLESS SWITCHING CONTROLLERS. William A. Wolovich and Alan B. Arehart 1. December 27, Abstract
BUMPLESS SWITCHING CONTROLLERS William A. Wolovich and Alan B. Arehart 1 December 7, 1995 Abstract This paper outlines the design of bumpless switching controllers that can be used to stabilize MIMO plants
More informationControl Systems Design
ELEC4410 Control Systems Design Lecture 14: Controllability Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 14: Controllability p.1/23 Outline
More information