= 0 otherwise. Eu(n) = 0 and Eu(n)u(m) = δ n m
|
|
- Sherman French
- 5 years ago
- Views:
Transcription
1 A-AE 567 Final Homework Spring 212 You will need Matlab and Simulink. You work must be neat and easy to read. Clearly, identify your answers in a box. You will loose points for poorly written work. You must work alone. Exam is due in a box outside my office ARMS 321 by 12: pm Friday May 4. NAME: 1
2 1 Problem Consider the joint density function f x,y (x, y) determined by f x,y (x, y) = cx if x 2 y 1 and x 1 = otherwise where c is a constant. Find: The constant c. The conditional expectation ĝ(y) = E(x y = y). Compute the error E(x ĝ(y)) 2. Find a polynomial p(y) = α + βy + γy 2 of degree at most two which is the best approximation of x, that is, E(x p(y)) 2 = inf{e(x p(y)) 2 : p is a polynomial of degree at most 2} Compute the error E(x p(y)) 2. 2 Problem Let u(n) be a discrete time white noise process, that is, Eu(n) = and Eu(n)u(m) = δ n m for all integers n and m where δ n is the Kronecker delta. Let θ be a uniform random variable over the interval [, 2π]. Consider the random process defined by y(n) = 2 sin(n + θ) + u(n) 2
3 where n is an integer. Moreover, assume that θ and u(n) are independent random variables for all integers n. Consider the prediction problem: σ 2 = min{e y(n) α 1 y(n 1) α 2 y(n 2) α 3 y(n 3) α 4 y(n 4) 2 : α j R}. The minimum is taken over all constants {α j } 4 1. (i) Find Ey(n)y(m) where n and m are integers. (ii) Find the optimal constants { α j } 4 1 and the error σ such that 4 σ 2 = E y(n) α j y(n j) 2. (iii) What is your guess for the error σ in the following optimization problem: σ 2 = min{e y(n) α j y(n j) 2 }. j=1 The theory behind this problem is beyond the scope of the course. However, you should be able to guess at σ. j=1 3 Problem Consider the discrete time system x(n + 1) = Ax(n) + Bu(n) and y(n) = Cx(n) + Dv(n) where u and v are independent Gaussian white noise processes and the initial condition is x() =. The system is third order. Moreover, the matrices {A, B, C, D} and the output {y(j)} 1 are given on the course webpage in the file a567 data12.mat. Let M n = span{y(j)} n. Find the following estimates which are vectors in R 3. P M1 x(9). 3
4 P M1 x(1) P M1 x(11) P M1 x(13) 4 Problem Consider the state space system x(n + 1) = Ax(n) + u(n) and y(n) = Cx(n) + v(n) (4.1) where u(n) and v(n) are mean zero Gaussian random process which are independent to the initial condition x(). Moreover, assume that E u(n) [ ] u(m) v(m) = R 11 R 12 δ n m. v(n) R 21 R 22 Here δ = 1 and for k nonzero, δ k =. Let M n = span{y(k)} n and x(n) = P Mn 1 x(n) denote the optimal state estimate. Let Q n = E x(n) x(n) = E (x(n) x(n)) (x(n) x(n)) be the error covariance matrix. Find the Kalman filter for the state space system in (4.1). To be precise, find a recursive estimate for the optimal state x(n) and a recursive formula for the error covariance Q n. (i) The recursion for the optimal state is x(n + 1) = (ii) The recursion for the error covariance is Q n+1 = (iii) Find x(n n) = P Mn x(n) in term of x(n), Q n and the data. (iv) Find x(n 1 n) = P Mn x(n 1) in term of x(n), Q n and the data. 4
5 5 Problem: The inverted pendulum problem Part of this problem is taken from the web site: Control tutorials for Matlab, Consider the classical inverted pendulum problem, consisting of a cart of mass m 1 with an inverted pendulum with inertia J on the cart. The linearized equations of motion in state space form are given by ẋ 1 x ẍ (J+m 2l 2 )ζ m 2 2 gl2 ẋ = q q + φ 1 φ φ m 2lζ m 2 gl(m 1 +m 2 ) q q φ y = x 1 ẋ v 1 1 φ 1 v 2 2 φ q = J(m 1 + m 2 ) + m 1 m 2 l 2. J+m 2 l 2 q m 2 l q 1 w + 4 u 1 u Here x is the position of the cart and φ is the angle of the inverted pendulum from the vertical position. The control force on the cart is denoted by w. Moreover, u 1, u 2, v 1 and v 2 are all independent white noise processes. Notice that u 1 is a random force on the cart and u 2 is a random torque on the pendulum. The notation is m 1 =.5kg is the mass of the cart m 2 =.2kg is the mass of the pendulum ζ =.1N/m/sec is the friction of the cart l =.3m is the distance to the pendulum s center of mass 5
6 J =.6kgm 2 is the inertia of the pendulum w is the force applied to the cart x is the position of the cart φ is the angle in radians of the pendulum from the vertical position. Assume that all the initial conditions are zero. Design a feedback controller w = K x based on the steady state Kalman filter such that the cart moves not more that one meter from the origin ( x 1) and the angle moves no more than ±2 degrees from ] the vertical position ( φ.35). Your state feedback gain K = [k 1 k 2 k 3 k 4 must satisfy k j 3 for j = 1, 2, 3, 4. Simulate your controller in Simulink for 3 seconds. Hand the graphs from your Simulink program for: (i) The position x and the estimated position of the cart x on the same graph. (ii) The velocity ẋ and the estimated velocity of the cart x on the same graph. (iii) The angle φ and the estimated angle of the pendulum φ on the same graph. (iv) The angular velocity φ and the estimated angular velocity of the pendulum φ on the same graph. (v) Hand in your gain K. On the band limited white nose generators, set both the noise power and sample time equal to.1. The seed for u 1, u 2, v 1 and v 2 is respectively, 23341, 23342, and
7 6 Problem Consider the third order state space system ẋ = Ax and y = Cx 1 [ ] A = 1 and C = Find the error σ and the optimal initial condition x to solve the following optimization problem σ 2 = e 2t Ce At x { } 2 dt = inf e 2t Ce At x 2 dt : x C 3 7 A review of steady state Kalman filtering Let us review the continuous time steady state Kalman filter. Throughout we assume that the system is time invariant, that is, ẋ = Ax + Bu and y = Cx + Dv (7.1) where A is an operator on X and B maps U into X while C maps X into Y and D maps V into Y. The operators {A, B, C, D} are all fixed and do not depend upon time t. We also assume that the operator D is onto Y, or equivalently, DD is invertible. Let Q(t) be the solution to the Riccati differential equation: Q = AQ + QA + BB QC (DD ) 1 CQ subject to Q() =. (7.2) The following is a classical result in linear quadratic regulator theory and Riccati differential equations. THEOREM 7.1 Consider the linear time invariant system {A, B, C, D} where D is onto. Then there exists a solution Q(t) to the Riccati differential equation in (7.2) 7
8 for all t. In other words, the Riccati differential equation in (7.2) does not have a finite escape time. Moreover, the following holds. (i) The solution {Q(t) : t } forms an increasing sequence of positive operators, that is, Q(σ) Q(t) for all σ t. (ii) If the pair {C, A} is observable, then Q(t) converges to a positive operator P as t tends to infinity, that is, P = lim t Q(t). (7.3) In this case, P is a positive solution for the algebraic Riccati equation = AP + P A + BB P C (DD ) 1 CP. (7.4) (iii) If {A, B, C, D} is controllable and observable, then P is strictly positive. (iv) If {A, B, C, D} is controllable and observable, then P is the only strictly positive solution to the algebraic Riccati equation (7.4). To be precise, if Ξ is a strictly positive solution to (7.4), then Ξ = P = lim t Q(t). (v) If {A, B, C, D} is controllable and observable, then A P C (DD ) 1 C is a continuous time stable operator, that is, all the eigenvalues of A P C (DD ) 1 C are contained in the open left half plane. 7.1 The steady state Kalman filter with control Assume that {A, B, C, D} is a controllable and observable system. Let P be the unique positive solution to the algebraic Riccati equation in (7.4). By passing limits in the Kalman filter, we arrive at the steady state Kalman filter defined by x = (A LC) x + Ly L = P C (DD ) 1. (7.5) 8
9 Theorem 7.1 guarantees that A LC is a continuous time stable operator. Moreover, L is called the Kalman gain or Kalman observer. The steady state Kalman filter provides an estimate x(t) for the state x(t) given the past {y(σ) : σ < t}. The Kalman filter converges to the steady state Kalman filter. In other words, the steady state Kalman filter converges to the optimal state estimator as t tends to infinity. Assume that {A, B, C, D} is a controllable and observable system. Consider the Kalman filtering problem with a control added to the state equation, that is, ẋ = Ax + Bu + B 1 w and y = Cx + Dv. (7.6) Here w is a known input or control with values in W, and B 1 is an operator mapping W into X. As before, u and v are independent white noise processes, which are also independent to the initial condition x(). In this case, the steady state Kalman filter is given by x = (A LC) x + Ly + B 1 w L = P C (DD ) 1. (7.7) As expected, P is the unique positive solution to the algebraic Riccati equation (7.4) in Theorem 7.1. Finally, it is noted that (7.7) reduces to the steady state Kalman filter in (7.5) when B 1 =. 9
10 simouty To Workspace3 output y simoutx1 To Workspace1 Band Limited White Noise2 x = Ax+Bu y = Cx+Du State Space em x = Ax+Bu y = Cx+Du em x1 c* u Kalman Filter d Matrix Gain1 x2 Band Limited White Noise1 Gain simoutx2 To Workspace2 Sine Wave Figure 1: The Simulink block diagram for w(t) = sin(t). 8 An example of a sinusoid input Consider the stable state space system given by ẋ1 = 1 x 1 + u + 1 sin(t) ẋ 2 x 2 y = x 1 x 2 + v. (8.1) Notice that { 1, 2} are the eigenvalues for the 2 2 state space matrix in (8.1). As expected, u and v are independent white noise processes. For simplicity we have assumed that all the initial conditions are zero. The known input w = sin(t). Using [ ] tr [ ] P = are(a, C C/(d 2 ), B B ) in Matlab with d = 1, B = 1 and C = 1 1 (the transpose of a matrix is denoted by tr), we arrived at.818 P =.31 and L = The eigenvalues for P are {.817,.166}. Since the eigenvalues for P are small, we expect the Kalman filter to perform rather well. The Simulink block diagram to 1
11 5 The plot of the output y 2 The plot of x 1 and hatx time time (a) The output y feed into the Kalman filter. (b) The state space x 1 and x The plot of x 2 and hatx time (c) The state space x 2 and x 2 Figure 2: The plots of y, x 1, x 1, x 2 and x 2. implement this Kalman filter is given in Figure 1. In this Simulink program, we used.1 for both the noise power and sample time for our white noise generators. The seed for the band limited white noise generator 1 is 23342, while the seed for the band limited white noise generator 2 is The output y for the state space model (8.1) is given in Figure 2, Part (a). The graphs for x 1 and its state estimate x 1 are presented in Figure 2, Part (b). Finally, the graphs for x 2 and its corresponding state estimate x 2 are given in Figure 2, Part (c). 11
12 9 State feedback and the Kalman filter As before, assume that {A, B, C, D} is a controllable and observable system. The state equations for the Kalman filtering problem with a control w are given by ẋ = Ax + Bu + B 1 w and y = Cx + Dv. (9.1) Here w is a known input or control with values in W, and B 1 is an operator mapping W into X. As before, u and v are independent white noise processes. Recall that the steady state Kalman filter is given by x = (A LC) x + Ly + B 1 w L = P C (DD ) 1 = AP + P A + BB P C (DD ) 1 CP where P >. (9.2) Now let us use the estimated state x to design a state feedback controller K. To this end, compute a state feedback controller K mapping X into W such that A B 1 K is stable. One can use pole placement or the linear quadratic regulator to determine a state feedback gain K satisfying a specified performance criteria. Since x is the estimated state, the feedback control is given by w = K x. Substituting w = K x into ẋ = Ax + Bu + B 1 w and (9.2) with y = Cx + Dv yields ẋ = A B 1K x x + B u. (9.3) x LC A LC B 1 K LD v So the feedback system is precisely a state space system driven by the white noise [ tr. process u v] Recall that two state space systems {A 2, B 2 } and {A 3, B 3 } are called similar if there exists an invertible transformation Ξ such that A 3 = ΞA 2 Ξ 1 and B 3 = ΞB 2. If {A 2, B 2 } and {A 3, B 3 } are similar, then A 2 and A 3 have the same eigenvalues. In particular, A 2 is stable if and only if A 3 is stable. 12
13 We claim that the state space system (9.3) is similar to the following state space system ẋ = A B 1K B 1 K x x + B u. (9.4) x A LC B LD v The error in estimation is denoted by x = x x. By construction A B 1 K and A LC are both stable. Since the state space matrix in (9.4) is upper triangular with A B 1 K and A LC on the diagonal, the state space system in (9.4) is stable. Because the two state space matrices in (9.3) and (9.4) are similar, the state space operator T = A LC B 1K on X (9.5) A LC B 1 K X in (9.3) must also be stable. In fact, the eigenvalues of T are precisely the eigenvalues of A B 1 K union the eigenvalues of A LC. To show that (9.3) and (9.4) are similar, consider the operator determined by Ξ = I on X. I I X Recall that x = x x is the error in estimation. Now observe that Ξ x x = I x x = x x. I I In other words, Ξ maps [ x ] tr [ tr. x into x x] The operator Ξ is invertible. In fact, Ξ 1 = Ξ. To verify that (9.3) is similar to (9.4), notice that ΞT Ξ 1 = I A B 1K Ξ I I LC A LC B 1 K = A B 1K I = A B 1K B 1 K A LC LC A I I A LC Ξ B = I B = B. LD I I LD B LD 13
14 This yields the state space system in (9.4) and proves our claim. 1 An example of state feedback control Consider the state space system ẋ1 = 1 x 1 + u + 2 w(t) ẋ 2 x 2 y = x 1 + x 2 + v. (1.1) The eigenvalues for the 2 2 state space matrix in (1.1) are {1, 2}. The state space system in (1.1) is unstable. Using P = are(a, C C/(d 2 ), B B ) in Matlab with [ ] tr [ ] d = 1, B = 1 and C = 1 1, we arrived at.5 P =.5 and L = The eigenvalues for the matrix P are {.4418, }. By employing the linear [ tr, quadratic regulator in Matlab K = lqr(a, B 1, eye(2), 1/4) with B 1 = 2 1] we obtained K = [ ] In this case, the matrix T in (9.5) is given by T = A B 1K =. LC A LC B 1 K The eigenvalues for T are { 1, 3.568±.484i, }. As expected, T is stable. The Simulink block diagram to implement this state feedback control w = K x with the Kalman filter to estimate x is given in Figure 3. In this Simulink program, we 14
15 used.1 for both the noise power and sample time for our white noise generators. Moreover, the seed for the band limited white noise generator 1 is while the seed for the band limited white noise generator 2 is The output y for the state space model (1.1) is given in Figure 4, Part (a). Figure 4, Part (b) graphs x 1 and its state estimate x 1, while Figure 4, Part (c) graphs x 2 and its state estimate x 2. Finally, it is noted that our estimate x 1 for the state space x 1 is very good, that is, the graphs for x 1 and x 1 are practically identical. 15
16 simout1 To Workspace1 output y simout2 To Workspace2 Band Limited White Noise2 x = Ax+Bu y = Cx+Du State Space em x = Ax+Bu y = Cx+Du em x1 c* u Kalman FIlter d Matrix Gain1 x2 Band Limited White Noise1 Gain simout k* u To Workspace Feedback Figure 3: The Simulink block diagram for w = K x. 16
17 3 The plot of the output y 2 The plot of x 1 and hatx time time (a) The output y feed into the Kalman filter. (b) The state space x 1 and x The plot of x 2 and hatx time (c) The state space x 2 and x 2 Figure 4: The plots of y, x 1, x 1, x 2 and x State feedback and the discrete Kalman filter Assume that {A, B, C, D} is a controllable and observable discrete time system. The state equations for the discrete time Kalman filtering problem with a control w are given by x(n + 1) = Ax(n) + Bu(n) + B 1 w(n) and y(n) = Cx(n) + Dv(n). (11.1) 17
18 Here w is a known input or control with values in W, and B 1 is an operator mapping W into X. As before, u and v are independent white noise processes. Recall that the steady state Kalman filter is given by x(n + 1) = (A LC) x(n) + Ly(n) + B 1 w(n) L = AP C (CP C + DD ) 1 (11.2) P = AP A + BB AP C (CP C + DD ) 1 CP A where P >. Finally, A LC is stable, that is all the eigenvalues of A LC are contained in the open unit disc. Now let us use the estimated state x to design a state feedback controller K. To this end, compute a state feedback controller K mapping X into W such that A B 1 K is stable. One can use pole placement or the discrete time linear quadratic regulator to determine a state feedback gain K satisfying a specified performance criteria. Since x is the steady state estimate, the feedback control is given by w = K x. Substituting w(n) = K x(n) into x(n + 1) = Ax(n) + Bu(n) + B 1 w(n) and (11.2) with y = Cx + Dv yields x(n + 1) = A B 1K x(n) + B u(n). (11.3) x(n + 1) LC A LC B 1 K x(n) LD v(n) So the feedback system is precisely a state space system driven by the white noise [ tr. process u v] Recall that two state space systems {A 2, B 2 } and {A 3, B 3 } are called similar if there exists an invertible transformation Ξ such that A 3 = ΞA 2 Ξ 1 and B 3 = ΞB 2. If {A 2, B 2 } and {A 3, B 3 } are similar, then A 2 and A 3 have the same eigenvalues. In particular, A 2 is stable if and only if A 3 is stable. 18
19 We claim that the state space system (11.3) is similar to the following state space system: x(n + 1) = A B 1K B 1 K x(n) + B u(n). (11.4) x(n + 1) A LC x(n) B LD v(n) The error in estimation is denoted by x = x x. By construction A B 1 K and A LC are both stable. Since the state space matrix in (11.4) is upper triangular with A B 1 K and A LC on the diagonal, the state space system in (11.4) is stable. Because the two state space matrices in (11.3) and (11.4) are similar, the state space operator T = A LC B 1K on X (11.5) A LC B 1 K X in (11.3) must also be stable. In fact, the eigenvalues of T are precisely the eigenvalues of A B 1 K union the eigenvalues of A LC. To show that (11.3) and (11.4) are similar, consider the operator determined by Ξ = I on X. I I X Recall that x = x x is the error in estimation. Now observe that Ξ x x = I x x = x x. I I [ ] tr [ tr. In other words, Ξ maps x x into x x] The operator Ξ is invertible. In fact, Ξ 1 = Ξ. To verify that (11.3) is similar to (11.4), notice that ΞT Ξ 1 = I A B 1K Ξ I I LC A LC B 1 K = A B 1K I = A B 1K B 1 K A LC LC A I I A LC Ξ B = I B = B. LD I I LD B LD 19
20 This yields the state space system in (11.4) and proves our claim. 12 An example of state feedback control Consider the state space system: x 1(n + 1) = 1 x 1(n) + u + 2 w(t) x 2 (n + 1) 2 3 x 2 (n) 1 1 y = x 1 + x 2 + v. (12.1) The eigenvalues for the 2 2 state space matrix in (12.1) are {.5616, }. The discrete state space system in (12.1) is unstable. Using P = dare(a, C, B B, d d ) [ ] tr [ ] in Matlab with d = 1, B = 1 and C = 1 1, we arrived at.649 P = and L = The eigenvalues for the matrix P are {.121, }. By employing the discrete time [ tr, linear quadratic regulator in Matlab K = dlqr(a, B 1, eye(2), 1/8) with B 1 = 2 1] we obtained K = [ ] In this case, the matrix T in (11.5) is given by T = A B 1K =. LC A LC B 1 K The eigenvalues for T are {.5738,.344,.1545,.2456}. As expected, T is stable. The Simulink block diagram to implement this state feedback control w = K x with the steady state Kalman filter to estimate x is given in Figure 5. In this Simulink 2
21 program, we used.1 for the sample time of the white noise or random number generators. Moreover, we set the seed equal to 2 and 22. Furthermore, we presented a Simulink program based on the equations in (11.3). The output y for the state space model (12.1) is given in Figure 6, Part (a). Figure 6, Part (b) graphs x 1 and its steady state estimate x 1, while Figure 6, Part (c) graphs x 2 and its steady state estimate x 2. Finally, it is noted that our steady state estimate x 1 for the state space x 1 is very good, that is, the graphs for x 1 and x 1 are practically identical. Let M n = span{y(j) : j n}. Recall that x(n) = P Mn 1 x(n) is the optimal state estimate given the past {y(j)} n 1. Now consider the state estimate x(n n) = P Mn x(n). Clearly, x(n n) is a better estimate of the state x(n) than x(n). According to the class notes in steady state x(n n) ( I P C (CP C + DD ) 1 C ) x(n) + P C (CP C + DD ) 1 y(n). This steady state estimate of x(n n) is also included in our Simulink program as the third variable going into the respective scopes. As expected, x(n n) is a slightly better estimate of the state x(n) than x(n); see Figure 6. 21
22 K*u Gain y output K*u Gain2 Random Number Random Number1 y(n)=cx(n)+du(n) x(n+1)=ax(n)+bu(n) Discrete State Space K*u Gain1 first state second state Figure 5: The Simulink block diagram for w = K x. 22
23 Time offset: Time offset: (a) The output y. (b) The state x 1, x 1 and x 1 (n n) Time offset: (c) The state x 2, x 2 and x 2 (n n) Figure 6: The plots of y, x 1, x 1, x 2 and x 2. 23
24 13 The Kalman filter and state feedback control Assume that {A, B, C, D} is a controllable and observable discrete time system. The state equations for the discrete time Kalman filtering problem with a control w are given by x(n + 1) = Ax(n) + Bu(n) + B 1 w(n) and y(n) = Cx(n) + Dv(n). (13.1) Here w is a known input or control with values in W, and B 1 is an operator mapping W into X. As before, u and v are independent white noise processes. Recall that the Kalman filter is given by x(n + 1) = (A L n C) x(n) + L n y(n) + B 1 w(n) L n = AQ n C (CQ n C + DD ) 1 (13.2) Q n+1 = AQ n A + BB A n C (C n C + DD ) 1 C n A. The initial condition Q = Ex()x(). Finally, if Q =, then Q n is increasing (Q n Q n+1 for all n), and Q n converges to P, the unique positive solution to the algebraic Riccati equation in (11.2). Now let us use the estimated state x to design a state feedback controller K. To this end, compute a state feedback controller K mapping X into W such that A B 1 K is stable. One can use pole placement or the discrete time linear quadratic regulator to determine a state feedback gain K satisfying a specified performance criteria. Since x is the estimated state, the feedback control is given by w = K x. Substituting w(n) = K x(n) into x(n + 1) = Ax(n) + Bu(n) + B 1 w(n) and (13.2) with y = Cx + Dv yields x(n + 1) = A B 1K x(n) + B u(n). (13.3) x(n + 1) L n C A L n C B 1 K x(n) LD v(n) 24
25 So the feedback system is precisely a state space system driven by the white noise [ tr. process u v] By repeating our previous arguments, with x = x x, we see that the system in (13.3) is similar to the following time varying state space system x(n + 1) = A B 1K B 1 K x(n) + B u(n). (13.4) x(n + 1) A L n C x(n) B L n D v(n) The example in (12.1) revisited. system in (12.1), that is, x 1(n + 1) = x 2 (n + 1) Now let us apply the Kalman filter to the 1 x 1(n) + u + 2 w(t) 2 3 x 2 (n) 1 1 y = x 1 + x 2 + v. (13.5) The Matlab commands we used to simulate the state space system in (13.3) are given by a = [, 1; 2, 3]; b = [; 1]; b1 = [2; 1]; c = [1, 1]; d = 1; q{1} = zeros(2, 2); for j = 1 : 1; q{j+1} = a q{j} a +b b a q{j} c inv(c q{j} c +d d ) c q{j} a ; end for j = 1 : 1; L{j} = a q{j} c inv(c q{j} c + d d ); end p = dare(a, c, b b, d d ); norm(p q{1}) k = dlqr(a, b1, eye(2), 1/8) for j = 1 : 1; t{j} = [a, b1 k; L{j} c, a L{j} c b1 k]; B{j} = [b, [; ]; [; ], L{j} d]; end for j = 1 : 1; u1{j} = randn(2, 1); end x1 = zeros(4, 1); for j = 1 : 1; x{j + 1} = t{j} x{j} + B{j} u1{j}; end 25
26 for j = 1 : 1; x1(j) = x{j}(1); x2(j) = x{j}(2); x3(j) = x{j}(3); x4(j) = x{j}(4); end Next we compute the steady state response in (11.3). L s = a p c inv(c p c + d d ); t s = [a, b1 k; L s c, a L s c b1 k]; B s = [b, [; ]; [; ], L s d]; x s {1} = zeros(4, 1); for j = 1 : 1; x s {j + 1} = t s x s {j} + B s u1{j}; end for j = 1 : 1; x s1 (j) = x s {j}(1); x s2 (j) = x s {j}(2); x s3 (j) = x s {j}(3); x s4 (j) = x s {j}(4); end The stairs command was used to generate the graphs. Notice that we computed the state x which comes from the system using Q n, and the steady state x s using P. In this example, P Q 1 where Q =. Figure 7 Parts (a) and (b) compare x 1 with the steady state x 1s, and respectively, x 2 with the steady state x 2s. Because Q n yields a more accurate estimate of the state than the steady state, in the beginning x 1 x 1s and x 2 x 2s. Since the feedback controller is trying to make the state as small as possible, using Q n yields a better result for about 1 steps. Then the system converges to the steady state response. Finally, Figure 7, Parts (c) and (d) compares x and x. The estimates are quite good. 26
27 The first state with Q n and P The second state with Q n and P (a) states x 1 and x 1s. (b) The states x 2 and x 2s. 8 The first state and its estimate 6 The second state and its estimate (c) The state x 1 and x 1. (d) The state x 2 and x 2 Figure 7: The plots of y, x 1, x 1, x 2 and x 2. 27
x(n + 1) = Ax(n) and y(n) = Cx(n) + 2v(n) and C = x(0) = ξ 1 ξ 2 Ex(0)x(0) = I
A-AE 567 Final Homework Spring 213 You will need Matlab and Simulink. You work must be neat and easy to read. Clearly, identify your answers in a box. You will loose points for poorly written work. You
More informationThe Control of an Inverted Pendulum
The Control of an Inverted Pendulum AAE 364L This experiment is devoted to the inverted pendulum. Clearly, the inverted pendulum will fall without any control. We will design a controller to balance the
More informationRECURSIVE ESTIMATION AND KALMAN FILTERING
Chapter 3 RECURSIVE ESTIMATION AND KALMAN FILTERING 3. The Discrete Time Kalman Filter Consider the following estimation problem. Given the stochastic system with x k+ = Ax k + Gw k (3.) y k = Cx k + Hv
More informationState space control for the Two degrees of freedom Helicopter
State space control for the Two degrees of freedom Helicopter AAE364L In this Lab we will use state space methods to design a controller to fly the two degrees of freedom helicopter. 1 The state space
More informationLecture 9. Introduction to Kalman Filtering. Linear Quadratic Gaussian Control (LQG) G. Hovland 2004
MER42 Advanced Control Lecture 9 Introduction to Kalman Filtering Linear Quadratic Gaussian Control (LQG) G. Hovland 24 Announcement No tutorials on hursday mornings 8-9am I will be present in all practical
More informationHere represents the impulse (or delta) function. is an diagonal matrix of intensities, and is an diagonal matrix of intensities.
19 KALMAN FILTER 19.1 Introduction In the previous section, we derived the linear quadratic regulator as an optimal solution for the fullstate feedback control problem. The inherent assumption was that
More informationThe Control of an Inverted Pendulum
The Control of an Inverted Pendulum AAE 364L This experiment is devoted to the inverted pendulum. Clearly, the inverted pendulum will fall without any control. We will design a controller to balance the
More informationCALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems. CDS 110b
CALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems CDS 110b R. M. Murray Kalman Filters 14 January 2007 Reading: This set of lectures provides a brief introduction to Kalman filtering, following
More informationLab 6a: Pole Placement for the Inverted Pendulum
Lab 6a: Pole Placement for the Inverted Pendulum Idiot. Above her head was the only stable place in the cosmos, the only refuge from the damnation of the Panta Rei, and she guessed it was the Pendulum
More informationState Regulator. Advanced Control. design of controllers using pole placement and LQ design rules
Advanced Control State Regulator Scope design of controllers using pole placement and LQ design rules Keywords pole placement, optimal control, LQ regulator, weighting matrixes Prerequisites Contact state
More informationOptimization-Based Control
Optimization-Based Control Richard M. Murray Control and Dynamical Systems California Institute of Technology DRAFT v1.7a, 19 February 2008 c California Institute of Technology All rights reserved. This
More informationCALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems
CDS 101 1. For each of the following linear systems, determine whether the origin is asymptotically stable and, if so, plot the step response and frequency response for the system. If there are multiple
More informationState Feedback MAE 433 Spring 2012 Lab 7
State Feedback MAE 433 Spring 1 Lab 7 Prof. C. Rowley and M. Littman AIs: Brandt Belson, onathan Tu Princeton University April 4-7, 1 1 Overview This lab addresses the control of an inverted pendulum balanced
More informationLaboratory 11 Control Systems Laboratory ECE3557. State Feedback Controller for Position Control of a Flexible Joint
Laboratory 11 State Feedback Controller for Position Control of a Flexible Joint 11.1 Objective The objective of this laboratory is to design a full state feedback controller for endpoint position control
More informationOptimization-Based Control
Optimization-Based Control Richard M. Murray Control and Dynamical Systems California Institute of Technology DRAFT v2.1a, January 3, 2010 c California Institute of Technology All rights reserved. This
More informationCONTROL DESIGN FOR SET POINT TRACKING
Chapter 5 CONTROL DESIGN FOR SET POINT TRACKING In this chapter, we extend the pole placement, observer-based output feedback design to solve tracking problems. By tracking we mean that the output is commanded
More informationCALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems. CDS 110b
CALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems CDS 110b R. M. Murray Kalman Filters 25 January 2006 Reading: This set of lectures provides a brief introduction to Kalman filtering, following
More 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 informationOptimization-Based Control
Optimization-Based Control Richard M. Murray Control and Dynamical Systems California Institute of Technology DRAFT v1.7b, 23 February 2008 c California Institute of Technology All rights reserved. This
More informationSteady State Kalman Filter
Steady State Kalman Filter Infinite Horizon LQ Control: ẋ = Ax + Bu R positive definite, Q = Q T 2Q 1 2. (A, B) stabilizable, (A, Q 1 2) detectable. Solve for the positive (semi-) definite P in the ARE:
More informationLQR, Kalman Filter, and LQG. Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin
LQR, Kalman Filter, and LQG Postgraduate Course, M.Sc. Electrical Engineering Department College of Engineering University of Salahaddin May 2015 Linear Quadratic Regulator (LQR) Consider a linear system
More informationSAMPLE SOLUTION TO EXAM in MAS501 Control Systems 2 Autumn 2015
FACULTY OF ENGINEERING AND SCIENCE SAMPLE SOLUTION TO EXAM in MAS501 Control Systems 2 Autumn 2015 Lecturer: Michael Ruderman Problem 1: Frequency-domain analysis and control design (15 pt) Given is a
More informationAutomatic Control II Computer exercise 3. LQG Design
Uppsala University Information Technology Systems and Control HN,FS,KN 2000-10 Last revised by HR August 16, 2017 Automatic Control II Computer exercise 3 LQG Design Preparations: Read Chapters 5 and 9
More informationMath 504 (Fall 2011) 1. (*) Consider the matrices
Math 504 (Fall 2011) Instructor: Emre Mengi Study Guide for Weeks 11-14 This homework concerns the following topics. Basic definitions and facts about eigenvalues and eigenvectors (Trefethen&Bau, Lecture
More informationEL 625 Lecture 10. Pole Placement and Observer Design. ẋ = Ax (1)
EL 625 Lecture 0 EL 625 Lecture 0 Pole Placement and Observer Design Pole Placement Consider the system ẋ Ax () The solution to this system is x(t) e At x(0) (2) If the eigenvalues of A all lie in the
More information6 OUTPUT FEEDBACK DESIGN
6 OUTPUT FEEDBACK DESIGN When the whole sate vector is not available for feedback, i.e, we can measure only y = Cx. 6.1 Review of observer design Recall from the first class in linear systems that a simple
More informationControls Problems for Qualifying Exam - Spring 2014
Controls Problems for Qualifying Exam - Spring 2014 Problem 1 Consider the system block diagram given in Figure 1. Find the overall transfer function T(s) = C(s)/R(s). Note that this transfer function
More informationLecture 19 IIR Filters
Lecture 19 IIR Filters Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/5/10 1 General IIR Difference Equation IIR system: infinite-impulse response system The most general class
More informationECEn 483 / ME 431 Case Studies. Randal W. Beard Brigham Young University
ECEn 483 / ME 431 Case Studies Randal W. Beard Brigham Young University Updated: December 2, 2014 ii Contents 1 Single Link Robot Arm 1 2 Pendulum on a Cart 9 3 Satellite Attitude Control 17 4 UUV Roll
More informationChapter 9 Observers, Model-based Controllers 9. Introduction In here we deal with the general case where only a subset of the states, or linear combin
Lectures on Dynamic Systems and Control Mohammed Dahleh Munther A. Dahleh George Verghese Department of Electrical Engineering and Computer Science Massachuasetts Institute of Technology c Chapter 9 Observers,
More informationA Crash Course on Kalman Filtering
A Crash Course on Kalman Filtering Dan Simon Cleveland State University Fall 2014 1 / 64 Outline Linear Systems Probability State Means and Covariances Least Squares Estimation The Kalman Filter Unknown
More informationECE-320: Linear Control Systems Homework 8. 1) For one of the rectilinear systems in lab, I found the following state variable representations:
ECE-30: Linear Control Systems Homework 8 Due: Thursday May 6, 00 at the beginning of class ) For one of the rectilinear systems in lab, I found the following state variable representations: 0 0 q q+ 74.805.6469
More informationControl Systems Design
ELEC4410 Control Systems Design Lecture 18: State Feedback Tracking and State Estimation Julio H. Braslavsky julio@ee.newcastle.edu.au School of Electrical Engineering and Computer Science Lecture 18:
More informationReview: control, feedback, etc. Today s topic: state-space models of systems; linearization
Plan of the Lecture Review: control, feedback, etc Today s topic: state-space models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered
More informationThe control of a gantry
The control of a gantry AAE 364L In this experiment we will design a controller for a gantry or crane. Without a controller the pendulum of crane will swing for a long time. The idea is to use control
More informationInverted Pendulum: State-Space Methods for Controller Design
1 de 12 18/10/2015 22:45 Tips Effects TIPS ABOUT BASICS HARDWARE INDEX NEXT INTRODUCTION CRUISE CONTROL MOTOR SPEED MOTOR POSITION SYSTEM MODELING ANALYSIS Inverted Pendulum: State-Space Methods for Controller
More informationLecture 2 and 3: Controllability of DT-LTI systems
1 Lecture 2 and 3: Controllability of DT-LTI systems Spring 2013 - EE 194, Advanced Control (Prof Khan) January 23 (Wed) and 28 (Mon), 2013 I LTI SYSTEMS Recall that continuous-time LTI systems can be
More informationChapter 3. LQ, LQG and Control System Design. Dutch Institute of Systems and Control
Chapter 3 LQ, LQG and Control System H 2 Design Overview LQ optimization state feedback LQG optimization output feedback H 2 optimization non-stochastic version of LQG Application to feedback system design
More informationROBOTICS Laboratory Problem 02
ROBOTICS 2015-2016 Laboratory Problem 02 Basilio Bona DAUIN PoliTo Problem formulation The planar system illustrated in Figure 1 consists of a cart C sliding with or without friction along the horizontal
More informationState Feedback and State Estimators Linear System Theory and Design, Chapter 8.
1 Linear System Theory and Design, http://zitompul.wordpress.com 2 0 1 4 State Estimator In previous section, we have discussed the state feedback, based on the assumption that all state variables are
More informationState Feedback Controller for Position Control of a Flexible Link
Laboratory 12 Control Systems Laboratory ECE3557 Laboratory 12 State Feedback Controller for Position Control of a Flexible Link 12.1 Objective The objective of this laboratory is to design a full state
More informationEEE582 Homework Problems
EEE582 Homework Problems HW. Write a state-space realization of the linearized model for the cruise control system around speeds v = 4 (Section.3, http://tsakalis.faculty.asu.edu/notes/models.pdf). Use
More informationCDS 101: Lecture 4.1 Linear Systems
CDS : Lecture 4. Linear Systems Richard M. Murray 8 October 4 Goals: Describe linear system models: properties, eamples, and tools Characterize stability and performance of linear systems in terms of eigenvalues
More informationOptimal control and estimation
Automatic Control 2 Optimal control and estimation Prof. Alberto Bemporad University of Trento Academic year 2010-2011 Prof. Alberto Bemporad (University of Trento) Automatic Control 2 Academic year 2010-2011
More informationSYSTEMTEORI - KALMAN FILTER VS LQ CONTROL
SYSTEMTEORI - KALMAN FILTER VS LQ CONTROL 1. Optimal regulator with noisy measurement Consider the following system: ẋ = Ax + Bu + w, x(0) = x 0 where w(t) is white noise with Ew(t) = 0, and x 0 is a stochastic
More informationẋ n = f n (x 1,...,x n,u 1,...,u m ) (5) y 1 = g 1 (x 1,...,x n,u 1,...,u m ) (6) y p = g p (x 1,...,x n,u 1,...,u m ) (7)
EEE582 Topical Outline A.A. Rodriguez Fall 2007 GWC 352, 965-3712 The following represents a detailed topical outline of the course. It attempts to highlight most of the key concepts to be covered and
More information1 Steady State Error (30 pts)
Professor Fearing EECS C28/ME C34 Problem Set Fall 2 Steady State Error (3 pts) Given the following continuous time (CT) system ] ẋ = A x + B u = x + 2 7 ] u(t), y = ] x () a) Given error e(t) = r(t) y(t)
More informationTHE REACTION WHEEL PENDULUM
THE REACTION WHEEL PENDULUM By Ana Navarro Yu-Han Sun Final Report for ECE 486, Control Systems, Fall 2013 TA: Dan Soberal 16 December 2013 Thursday 3-6pm Contents 1. Introduction... 1 1.1 Sensors (Encoders)...
More informationIntroduction to Modern Control MT 2016
CDT Autonomous and Intelligent Machines & Systems Introduction to Modern Control MT 2016 Alessandro Abate Lecture 2 First-order ordinary differential equations (ODE) Solution of a linear ODE Hints to nonlinear
More informationComputationally, diagonal matrices are the easiest to work with. With this idea in mind, we introduce similarity:
Diagonalization We have seen that diagonal and triangular matrices are much easier to work with than are most matrices For example, determinants and eigenvalues are easy to compute, and multiplication
More informationPredictive Control of Gyroscopic-Force Actuators for Mechanical Vibration Damping
ARC Centre of Excellence for Complex Dynamic Systems and Control, pp 1 15 Predictive Control of Gyroscopic-Force Actuators for Mechanical Vibration Damping Tristan Perez 1, 2 Joris B Termaat 3 1 School
More informationProfessor Fearing EE C128 / ME C134 Problem Set 10 Solution Fall 2010 Jansen Sheng and Wenjie Chen, UC Berkeley
Professor Fearing EE C28 / ME C34 Problem Set Solution Fall 2 Jansen Sheng and Wenjie Chen, UC Berkeley. (5 pts) Final Value Given the following continuous time (CT) system ẋ = Ax+Bu = 5 9 7 x+ u(t), y
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 informationLinearization problem. The simplest example
Linear Systems Lecture 3 1 problem Consider a non-linear time-invariant system of the form ( ẋ(t f x(t u(t y(t g ( x(t u(t (1 such that x R n u R m y R p and Slide 1 A: f(xu f(xu g(xu and g(xu exist and
More informationThe Hilbert Space of Random Variables
The Hilbert Space of Random Variables Electrical Engineering 126 (UC Berkeley) Spring 2018 1 Outline Fix a probability space and consider the set H := {X : X is a real-valued random variable with E[X 2
More informationLab 5a: Pole Placement for the Inverted Pendulum
Lab 5a: Pole Placement for the Inverted Pendulum November 1, 2011 1 Purpose The objective of this lab is to achieve simultaneous control of both the angular position of the pendulum and horizontal position
More informationComputer Problem 1: SIE Guidance, Navigation, and Control
Computer Problem 1: SIE 39 - Guidance, Navigation, and Control Roger Skjetne March 12, 23 1 Problem 1 (DSRV) We have the model: m Zẇ Z q ẇ Mẇ I y M q q + ẋ U cos θ + w sin θ ż U sin θ + w cos θ θ q Zw
More informationCDS 101/110: Lecture 3.1 Linear Systems
CDS /: Lecture 3. Linear Systems Goals for Today: Describe and motivate linear system models: Summarize properties, examples, and tools Joel Burdick (substituting for Richard Murray) jwb@robotics.caltech.edu,
More informationMachine Learning 4771
Machine Learning 4771 Instructor: ony Jebara Kalman Filtering Linear Dynamical Systems and Kalman Filtering Structure from Motion Linear Dynamical Systems Audio: x=pitch y=acoustic waveform Vision: x=object
More informationLecture 7 and 8. Fall EE 105, Feedback Control Systems (Prof. Khan) September 30 and October 05, 2015
1 Lecture 7 and 8 Fall 2015 - EE 105, Feedback Control Systems (Prof Khan) September 30 and October 05, 2015 I CONTROLLABILITY OF AN DT-LTI SYSTEM IN k TIME-STEPS The DT-LTI system is given by the following
More informationChap. 3. Controlled Systems, Controllability
Chap. 3. Controlled Systems, Controllability 1. Controllability of Linear Systems 1.1. Kalman s Criterion Consider the linear system ẋ = Ax + Bu where x R n : state vector and u R m : input vector. A :
More informationRepresent this system in terms of a block diagram consisting only of. g From Newton s law: 2 : θ sin θ 9 θ ` T
Exercise (Block diagram decomposition). Consider a system P that maps each input to the solutions of 9 4 ` 3 9 Represent this system in terms of a block diagram consisting only of integrator systems, represented
More informationEE226a - Summary of Lecture 13 and 14 Kalman Filter: Convergence
1 EE226a - Summary of Lecture 13 and 14 Kalman Filter: Convergence Jean Walrand I. SUMMARY Here are the key ideas and results of this important topic. Section II reviews Kalman Filter. A system is observable
More informationFinal Exam Solutions
EE55: Linear Systems Final Exam SIST, ShanghaiTech Final Exam Solutions Course: Linear Systems Teacher: Prof. Boris Houska Duration: 85min YOUR NAME: (type in English letters) I Introduction This exam
More informationMSMS Matlab Problem 02
MSMS 2014-2015 Matlab Problem 02 Basilio Bona DAUIN PoliTo Problem formulation The planar system illustrated in Figure 1 consists of a cart C sliding with friction along the horizontal rail; the cart supports
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 informationCDS 110b: Lecture 2-1 Linear Quadratic Regulators
CDS 110b: Lecture 2-1 Linear Quadratic Regulators Richard M. Murray 11 January 2006 Goals: Derive the linear quadratic regulator and demonstrate its use Reading: Friedland, Chapter 9 (different derivation,
More informationControl Systems Design, SC4026. SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft
Control Systems Design, SC4026 SC4026 Fall 2009, dr. A. Abate, DCSC, TU Delft Lecture 4 Controllability (a.k.a. Reachability) vs Observability Algebraic Tests (Kalman rank condition & Hautus test) A few
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 informationDenis ARZELIER arzelier
COURSE ON LMI OPTIMIZATION WITH APPLICATIONS IN CONTROL PART II.2 LMIs IN SYSTEMS CONTROL STATE-SPACE METHODS PERFORMANCE ANALYSIS and SYNTHESIS Denis ARZELIER www.laas.fr/ arzelier arzelier@laas.fr 15
More informationEC Control Engineering Quiz II IIT Madras
EC34 - Control Engineering Quiz II IIT Madras Linear algebra Find the eigenvalues and eigenvectors of A, A, A and A + 4I Find the eigenvalues and eigenvectors of the following matrices: (a) cos θ sin θ
More informationD(s) G(s) A control system design definition
R E Compensation D(s) U Plant G(s) Y Figure 7. A control system design definition x x x 2 x 2 U 2 s s 7 2 Y Figure 7.2 A block diagram representing Eq. (7.) in control form z U 2 s z Y 4 z 2 s z 2 3 Figure
More informationCALIFORNIA INSTITUTE OF TECHNOLOGY Control and Dynamical Systems
CDS 101 1. Åström and Murray, Exercise 1.3 2. Åström and Murray, Exercise 1.4 3. Åström and Murray, Exercise 2.6, parts (a) and (b) CDS 110a 1. Åström and Murray, Exercise 1.4 2. Åström and Murray, Exercise
More informationLecture 8. Applications
Lecture 8. Applications Ivan Papusha CDS270 2: Mathematical Methods in Control and System Engineering May 8, 205 / 3 Logistics hw7 due this Wed, May 20 do an easy problem or CYOA hw8 (design problem) will
More informationControl Systems Design, SC4026. SC4026 Fall 2010, dr. A. Abate, DCSC, TU Delft
Control Systems Design, SC4026 SC4026 Fall 2010, dr. A. Abate, DCSC, TU Delft Lecture 4 Controllability (a.k.a. Reachability) and Observability Algebraic Tests (Kalman rank condition & Hautus test) A few
More information1 Lyapunov theory of stability
M.Kawski, APM 581 Diff Equns Intro to Lyapunov theory. November 15, 29 1 1 Lyapunov theory of stability Introduction. Lyapunov s second (or direct) method provides tools for studying (asymptotic) stability
More informationKalman Filter and Parameter Identification. Florian Herzog
Kalman Filter and Parameter Identification Florian Herzog 2013 Continuous-time Kalman Filter In this chapter, we shall use stochastic processes with independent increments w 1 (.) and w 2 (.) at the input
More informationSuppose that we have a specific single stage dynamic system governed by the following equation:
Dynamic Optimisation Discrete Dynamic Systems A single stage example Suppose that we have a specific single stage dynamic system governed by the following equation: x 1 = ax 0 + bu 0, x 0 = x i (1) where
More informationLinear Systems. Manfred Morari Melanie Zeilinger. Institut für Automatik, ETH Zürich Institute for Dynamic Systems and Control, ETH Zürich
Linear Systems Manfred Morari Melanie Zeilinger Institut für Automatik, ETH Zürich Institute for Dynamic Systems and Control, ETH Zürich Spring Semester 2016 Linear Systems M. Morari, M. Zeilinger - Spring
More informationEE363 homework 2 solutions
EE363 Prof. S. Boyd EE363 homework 2 solutions. Derivative of matrix inverse. Suppose that X : R R n n, and that X(t is invertible. Show that ( d d dt X(t = X(t dt X(t X(t. Hint: differentiate X(tX(t =
More informationDr Ian R. Manchester Dr Ian R. Manchester AMME 3500 : Review
Week Date Content Notes 1 6 Mar Introduction 2 13 Mar Frequency Domain Modelling 3 20 Mar Transient Performance and the s-plane 4 27 Mar Block Diagrams Assign 1 Due 5 3 Apr Feedback System Characteristics
More informationTime-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 2015
Time-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 15 Asymptotic approach from time-varying to constant gains Elimination of cross weighting
More information1 Kalman Filter Introduction
1 Kalman Filter Introduction You should first read Chapter 1 of Stochastic models, estimation, and control: Volume 1 by Peter S. Maybec (available here). 1.1 Explanation of Equations (1-3) and (1-4) Equation
More informationTime-Invariant Linear Quadratic Regulators!
Time-Invariant Linear Quadratic Regulators Robert Stengel Optimal Control and Estimation MAE 546 Princeton University, 17 Asymptotic approach from time-varying to constant gains Elimination of cross weighting
More informationHomework Solution # 3
ECSE 644 Optimal Control Feb, 4 Due: Feb 17, 4 (Tuesday) Homework Solution # 3 1 (5%) Consider the discrete nonlinear control system in Homework # For the optimal control and trajectory that you have found
More informationState Feedback and State Estimators Linear System Theory and Design, Chapter 8.
1 Linear System Theory and Design, http://zitompul.wordpress.com 2 0 1 4 2 Homework 7: State Estimators (a) For the same system as discussed in previous slides, design another closed-loop state estimator,
More informationMath Assignment 3 - Linear Algebra
Math 216 - Assignment 3 - Linear Algebra Due: Tuesday, March 27. Nothing accepted after Thursday, March 29. This is worth 15 points. 10% points off for being late. You may work by yourself or in pairs.
More informationLinear-Quadratic-Gaussian (LQG) Controllers and Kalman Filters
Linear-Quadratic-Gaussian (LQG) Controllers and Kalman Filters Emo Todorov Applied Mathematics and Computer Science & Engineering University of Washington Winter 204 Emo Todorov (UW) AMATH/CSE 579, Winter
More informationSignals and Systems. Problem Set: The z-transform and DT Fourier Transform
Signals and Systems Problem Set: The z-transform and DT Fourier Transform Updated: October 9, 7 Problem Set Problem - Transfer functions in MATLAB A discrete-time, causal LTI system is described by the
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 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 informationFINAL EXAM SOLUTIONS, MATH 123
FINAL EXAM SOLUTIONS, MATH 23. Find the eigenvalues of the matrix ( 9 4 3 ) So λ = or 6. = λ 9 4 3 λ = ( λ)( 3 λ) + 36 = λ 2 7λ + 6 = (λ 6)(λ ) 2. Compute the matrix inverse: ( ) 3 3 = 3 4 ( 4/3 ) 3. Let
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 informationUniversity of Toronto Department of Electrical and Computer Engineering ECE410F Control Systems Problem Set #3 Solutions = Q o = CA.
University of Toronto Department of Electrical and Computer Engineering ECE41F Control Systems Problem Set #3 Solutions 1. The observability matrix is Q o C CA 5 6 3 34. Since det(q o ), the matrix is
More informationCourse Outline. Higher Order Poles: Example. Higher Order Poles. Amme 3500 : System Dynamics & Control. State Space Design. 1 G(s) = s(s + 2)(s +10)
Amme 35 : System Dynamics Control State Space Design Course Outline Week Date Content Assignment Notes 1 1 Mar Introduction 2 8 Mar Frequency Domain Modelling 3 15 Mar Transient Performance and the s-plane
More informationSimple Harmonic Motion Concept Questions
Simple Harmonic Motion Concept Questions Question 1 Which of the following functions x(t) has a second derivative which is proportional to the negative of the function d x! " x? dt 1 1. x( t ) = at. x(
More informationLecture 12. Upcoming labs: Final Exam on 12/21/2015 (Monday)10:30-12:30
289 Upcoming labs: Lecture 12 Lab 20: Internal model control (finish up) Lab 22: Force or Torque control experiments [Integrative] (2-3 sessions) Final Exam on 12/21/2015 (Monday)10:30-12:30 Today: Recap
More informationFeedback Control part 2
Overview Feedback Control part EGR 36 April 19, 017 Concepts from EGR 0 Open- and closed-loop control Everything before chapter 7 are open-loop systems Transient response Design criteria Translate criteria
More informationRobotics. Control Theory. Marc Toussaint U Stuttgart
Robotics Control Theory Topics in control theory, optimal control, HJB equation, infinite horizon case, Linear-Quadratic optimal control, Riccati equations (differential, algebraic, discrete-time), controllability,
More informationLecture 9 Nonlinear Control Design
Lecture 9 Nonlinear Control Design Exact-linearization Lyapunov-based design Lab 2 Adaptive control Sliding modes control Literature: [Khalil, ch.s 13, 14.1,14.2] and [Glad-Ljung,ch.17] Course Outline
More information