9 The LQR Problem Revisited
|
|
- Evangeline Jennings
- 5 years ago
- Views:
Transcription
1 9 he LQR Problem Revisited Problem: Compute a state feedback controller u(t) = Kx(t) that stabilizes the closed loop system and minimizes J := lim E [ z(t) z(t) ] t for the LI system Assumptions: a) D zu D zu 0; b) (A, B u ) stabilizable; ẋ(t) = Ax(t) + B u u(t) + B w w(t), x(0) = 0 z(t) = C z x(t) + D zu u(t). c) w(t) is a Gaussian white noise vector with zero mean and covariance W 0. Solution: As before, restate the cost function with the help of the Gramian J = trace ( ) PB w WBw where (A + B u K) P + P(A + B u K) + (C z + D zu ) (C z + D zu ) = 0. Use the Comparison Lemma to conclude that for any K such that (A + B u K) is Hurwitz then Y P = trace ( Y B w WB w) J where Y satisfies the inequality (A + B u K) Y + Y (A + B u K) + (C z + D zu K) (C z + D zu K) 0. MAE 280 B 124 Maurício de Oliveira
2 9.1 LQR revisited: solution by congruence + change-of-variables Recipe I (congruence + change-of-variables): First apply Schur complement (A + B u K) Y + Y (A + B u K) + (C z + D zu K) (C z + D zu K) 0 [ ] (A + Bu K) Y + Y (A + B u K) Cz + K Dzu 0 C z + D zu K I Under the assumption that (A, B u ) is stabilizable and Dzu D zu 0 we have that Y P 0 so that we can apply the congruence transformation [ ][ ][ ] Y 1 0 (A + Bu K) Y + Y (A + B u K) Cz + K Dzu Y I C z + D zu K I 0 I [ (A + Bu K)Y 1 + Y 1 (A + B u K) Y 1 Cz + Y ] 1 K Dzu C z Y 1 + D zu KY 1 0 I Apply change-of-variables X := Y 1, L := KY 1 to obtain the LMI [ AX + XA + B u L + L Bu XCz + ] L Dzu 0 C z X + D zu L I As for the cost function trace ( X 1 B w WB w) = trace ( Y Bw WB w) J Now introduce an extra variable Z such that herefore Z B wx 1 B w trace(zw) = trace ( B w X 1 B w W ) J Using Schur Complement Z B wx 1 B w [ ] Z B w 0. X B w MAE 280 B 125 Maurício de Oliveira
3 9.1.1 Summary: LQR revisited (first form) he state feedback controller u(t) = Kx(t), K = LX 1 where X S n, L R m n and Z S r solve the SDP min X S n,l R m n,z S r s.t. trace(zw) [ ] AX + XA + B u L + L Bu XCz + L Dzu 0, C z X + D zu L I [ ] Z B w 0, X 0. X B w stabilizes the closed loop system and minimizes the cost function J := lim E [ z(t) z(t) ] t for the LI system under the assumptions a) D zud zu 0; b) (A, B u ) stabilizable; ẋ(t) = Ax(t) + B u u(t) + B w w(t), x(0) = 0 z(t) = C z x(t) + D zu u(t) c) w(t) is a Gaussian white noise vector with zero mean and covariance W 0. MAE 280 B 126 Maurício de Oliveira
4 9.2 Elimination Lemma Let x R n, Q S n, B R m n and B R p n such that rank(b) < n and rank(c) < n. he following statements are equivalent: i) B QB 0 and C QC 0 ii) µ R : Q µb B 0, and Q µc C 0. iii) X R p m : Q + C X B + B X C 0. Remark: Proofs can be found in Skelton, Iwasaki and Grigoriadis (heorem , p. 29.) or Boyd, El Gahoui, Feron & Balakrishnan (p. 32). None as nice as the next one :). Proof: ii) and iii) i) is as in Finsler s Lemma, just multiply by B and C. i) ii) use Finsler s Lemma to show that i) implies the existence of scalars µ 1 and µ 2 such that Q µ 1 B B 0, Q µ 2 C C 0. hen µ = max{µ 1, µ 2 } is such that ii) holds for some µ. ii) iii) requires an explicit construction of X. MAE 280 B 127 Maurício de Oliveira
5 9.2.1 Construction of X Assume that ii) holds. If µ 0 then X = 0 solves the problem. So assume µ > 0. Now let ρ 2 = µ so that Define hen Q ρ 2 B B 0, Q ρ 2 C C 0. Φ := Q ρ 2 B B ρ 2 C C Φ + ρ 2 C C 0, Φ + ρ 2 B B 0. Also notice that Φ 0 because Φ := 1 2 [( Q 2ρ 2 B B ) + ( Q 2ρ 2 C C )] and 2ρ 2 > ρ 2 = µ. herefore, using Schur Complement [ ] Φ ρc Φ + ρ 2 C C 0 0. ρc I Schur Complement one more time [ ] Φ ρc 0 ρc I I ρ 2 CΦ 1 C 0. Likewise Φ + ρ 2 B B 0 I ρ 2 BΦ 1 B 0. We will now prove that X = ρ 2 Y Y = ρ 2 C Φ 1 B satisfies iii). With that in mind notice that [ ] I ρ 2 B Φ 1 B Y ρ 2 B Φ 1 C Y ρ 2 C Φ 1 B I ρ 2 C Φ 1 C 0. MAE 280 B 128 Maurício de Oliveira
6 Using Schur Complement this is equivalent to Φ ρ B ρ C ρ B I Y 0 ρ C Y I Using Schur Complement one more time [ ] [ ] Φ + ρ 2 C C ρ B + ρ C Y Q ρ ρ B + ρy C I + Y = 2 B B ρ B + ρ C Y Y ρ B + ρy C I + Y 0 Y and once again Q ρ 2 B B + ρ 2 (B + C Y)(I Y Y) 1 (B + Y C) 0 and I Y Y 0. Using the Lemma of the Matrix Inverse 2 so that (I Y Y) 1 = I + Y (I YY ) 1 Y 0 Q ρ 2 B B + ρ 2 (B + C Y)(I Y Y) 1 (B + Y C) = Q + ρ 2 C YB + ρ 2 B Y C + ρ 2 C YY C hence iii) holds because + ρ 2 (B + C Y)Y (I Y Y) 1 Y(B + Y C) ρ 2 C YY C + ρ 2 (B + C Y)Y (I Y Y) 1 Y(B + Y C) 0. 2 Verify that [ I + Y (I YY ) 1 Y ] (I Y Y) = I Y Y + Y (I YY ) 1 Y(I Y Y) = I Y Y + Y (I YY ) 1 (I YY )Y = I Y Y + Y Y = I. MAE 280 B 129 Maurício de Oliveira
7 9.3 LQR revisited: solution by Elimination Lemma Recipe II (existence conditions): As before we start with [ ] (A + Bu K) Y 0, Y + Y (A + B u K) Cz + K Dzu 0 C z + D zu K I Note that we are using strict inequalities! Now rewrite [ ] [ ] [ ] A Y + Y A Cz Y Bu [ ] I + C z I D zu }{{} K I 0 + }{{} 0 }{{}}{{} X C }{{} Q B C We need a more general version of Finsler s Lemma. K }{{} X Computing B and C under no assumptions [ B u Y Dzu] 0. }{{} B First compute C = [ I 0 ] = [ 0 I ]. he problem of computing a basis for 0 = Bx = [ [ ] Bu Y Dzu] [ ] x = B u Dzu Y 0 x 0 I is more involved. As before, because Y is nonsingular [ ] [ B Y 0 u Dzu] y = 0, y = 0 I x 0. Let [ E1 E 2 ] := [ ] Bu Dzu so that B = [ ] Y 1 0 [B u 0 I D zu [ ] ] Y 1 E = 1. E 2 MAE 280 B 130 Maurício de Oliveira
8 his choice produce the inequalities and 0 B QB = [ ] [ ][ ] E1 Y 1 E2 A Y + Y A Cz Y 1 E 1 C z I E 2 = [ ] [ ][ ] E1 E2 Y 1 A + AY 1 Y 1 Cz E1 C z Y 1. I 0 C QC = [ 0 I ] [ ][ ] A Y + Y A Cz 0 C z I I = I which is trivially satisfied. he first inequality, along with the cost inequality derived previously can be linearized with the change-of-variables X := Y 1. E 2 MAE 280 B 131 Maurício de Oliveira
9 9.3.2 Summary: LQR revisited (second form) he optimal state feedback controller u(t) = Kx(t) which can be computed from the solution to the SDP in the variables X S n, Z S r solve the SDP min X S n,z S r s.t. trace(zw) [ ] [ E 1 E2 AX + XA XCz C z X I [ ] Z B w 0, X 0. X B w ][ E1 E 2 ] 0, where [ E1 E 2 ] = [ ] Bu Dzu stabilizes the closed loop system and minimizes the cost function J := lim E [ z(t) z(t) ] t for the LI system under the assumptions a) D zu D zu 0; b) (A, B u ) stabilizable; ẋ(t) = Ax(t) + B u u(t) + B w w(t), x(0) = 0 z(t) = C z x(t) + D zu u(t) c) w(t) is a Gaussian white noise vector with zero mean and covariance W 0. MAE 280 B 132 Maurício de Oliveira
10 9.3.3 Computing B under the assumption D zu D zu 0 Proceed as before but now notice that if Dzu D zu 0 then [ ] [ ] E1 I 0 = E 2 D zu (DzuD zu ) 1 Bu (Dzu) = [ ] Bu Dzu is a possible choice for E 1 and E 2 so that, using X := Y 1 we have where 0 B QB = [ ] [ ][ ] E1 E2 AY 1 + Y 1 A + Y 1 Cz E1 C z Y 1 I E 2 [ ] I 0 = D zu (Dzu D zu) 1 Bu (Dzu [ ) ][ ] AX + XA XCz I 0 C z X I D zu (Dzu D zu) 1 Bu (Dzu [ ] ) I 0 = D zu (DzuD zu ) 1 Bu (Dzu) [ ] AX + XÃ XCz (Dzu) C z X + D zu (Dzu D zu) 1 Bu (Dzu [ÃX ) ] + X Ã B u (Dzu D zu) 1 Bu XCz (D zu ) = (D zu ) C z X (D zu ) (D zu ) Ã := A B u (D zud zu ) 1 D zuc z. MAE 280 B 133 Maurício de Oliveira
11 9.3.4 Summary: LQR revisited (third form) he optimal state feedback controller u(t) = Kx(t) which can be computed from the solution to the SDP in the variables X S n, Z S r solve the SDP min X S n,z S r s.t. trace(zw) ] [ÃX + X à B u (Dzu D zu) 1 Bu XCz (D zu ) (Dzu) C z X (Dzu) (Dzu) 0, [ ] Z B w 0, X 0. X B w where à := A B u (D zud zu ) 1 D zuc z, stabilizes the closed loop system and minimizes the cost function J := lim E [ z(t) z(t) ] t for the LI system under the assumptions a) D zu D zu 0; b) (A, B u ) stabilizable; ẋ(t) = Ax(t) + B u u(t) + B w w(t), x(0) = 0 z(t) = C z x(t) + D zu u(t) c) w(t) is a Gaussian white noise vector with zero mean and covariance W 0. MAE 280 B 134 Maurício de Oliveira
12 9.4 Example: satellite in circular orbit r u 1 u 2 m θ e Satellite of mass m with thrust in the radial direction u 1 and in the tangential direction u 2. Continuing... m( r r θ 2 ) = u 1 km r 2 + w 1, m(2ṙ θ + r θ) = u 2 + w 2, where w 1 and w 2 are independent white noise disturbances with variances δ 1 and δ 2. As before, putting in state space and linearizing around r ) ( ) x(t) = ωt (ū1 0, ū = w1 = 0, w = = 0, ω = k/ r ū 2 w 3 2 ω one obtains the linearized system ẋ = 3 ω 2 2 r ω 2 ω/ r 0 Problem: Given x+ 1/m 0 0 1/(m r) ( u1 u 2 ) + 1/m 0 0 1/(m r) m = 100 kg, r = R km, k = GM where G N m 2 /kg 2 is the universal gravitational constant, and M kg and R km are the mass and radius of the earth. If the variances δ 1 and δ 2 are 0.1N find solutions to the LQR control problem where Q = I, using u 2 only first, then using u 1 and u 2. R = ρi, MAE 280 B 135 Maurício de Oliveira ( w1 w 2 ).
13 % MAE 280 B - Linear Control Design % Mauricio de Oliveira % % LMI LQR Control % m = 100; % 100 kg r = 300E3; % 300 km R = 6.37E6; % ˆ3 km G = 6.673E-11; % N mˆ2/kgˆ2 M = 5.98E24; % ˆ24 kg k = G * M; % gravitational force constant w = sqrt(k/((r+r)ˆ3)); % angular velocity (rad/s) v = w * (R + r); % "ground" velocity (m/s) % linearized system matrices A = [ 1 0; 0 1; 3*wˆ2 2*(r+R)*w; -2*w/(r+R) 0]; Bu = [; ; 1/m 0; 0 1/(m*r)]; Bw = [; ; 1/m 0; 0 1/(m*r)]; % noise variances W = 0.1 * eye(2) W = e e-01 % scale = diag([1 r 1 r]) = % similarity transformation At = * A / At = e e e e e-04 0 But = * Bu But = e e-02 Bwt = * Bw Bwt = e e-02 MAE 280 B 136 Maurício de Oliveira
14 % optimal LQR state feedback control (using u2 only) n = size(at,1); m = 1; % LQR solution Cz = [eye(n); zeros(m,n)]; Dzu = [zeros(n,m); eye(m)]; [K,X,S] = lqr(at, But(:,2), Cz * Cz, Dzu * Dzu); Klqr = - K Klqr = e e e e+01 % LMI solution % declare variables X = sdpvar(n,n, symmetric ); Z = sdpvar(m,m, symmetric ); L = sdpvar(m,n); % declare LMIs LMI1 = [At*X+X*At +But(:,2)*L+L *But(:,2), X*Cz +L *Dzu ;... Cz*X+Dzu*L, -eye(m+n)]; LMI2 = [Z, Bwt(:,2) ; Bwt(:,2), X]; LMI3 = X; LMI = set(lmi1 < 0) + set(lmi2 > 0) + set(lmi3 > 0); options = sdpsettings( solver, sedumi ); solution =solvesdp(lmi,trace(w(2,2)*z),options) SeDuMi 1.1R3 by AdvOL, 2006 and Jos F. Sturm, Alg = 2: xz-corrector, theta = 0.250, beta = eqs m = 15, order n = 19, dim = 123, blocks = 4 nnz(a) = , nnz(ada) = 217, nnz(l) = 116 it : b*y gap delta rate t/tp* t/td* feas cg cg prec 0 : 4.49E : -6.05E E E+00 2 : -1.66E E E-01 3 : -3.52E E E-01 4 : -1.59E E E-01 5 : -2.29E E E-01 6 : -3.45E E E-01 7 : -4.71E E E-01 8 : -5.65E E E-01 9 : -7.25E E E : -8.33E E E : -1.07E E E : -1.23E E E : -1.40E E E : -1.52E E E : -1.63E E E : -1.70E E E : -1.76E E E : -1.79E E E : -1.82E E E : -1.84E E E-07 MAE 280 B 137 Maurício de Oliveira
15 21 : -1.84E E E : -1.85E E E : -1.85E E E : -1.85E E E : -1.85E E E : -1.85E E E : -1.85E E E : -1.85E E E-10 iter seconds digits c*x b*y Inf e e-02 Ax-b = 8.8e-10, [Ay-c]_+ = 3.0E-10, x = 5.0e+03, y = 2.4e-01 Detailed timing (sec) Pre IPM Post 1.500E E E-02 Max-norms: b = e-01, c = 1, Cholesky add =0, skip = 0, L.L = solution = yalmiptime: e-01 solvertime: e-01 info: Numerical problems (SeDuMi-1.1) problem: 4 dimacs: [7.9553e e e e-08] % Construct K and check closed loop stability Klqr Klqr = e e e e+01 K = double(l) / double(x) K = e e e e+01 Acl = At + But(:,2)*K; eig(acl) ans = e e-02i e e-02i e e-05 % optimal LQR state feedback control (using u1 and u2) n = size(at,1); m = 2; % LQR solution Cz = [eye(n); zeros(m,n)]; Dzu = [zeros(n,m); eye(m)]; [K,X,S] = lqr(at, But, Cz * Cz, Dzu * Dzu); Klqr = - K Klqr = e e e e e e e e+01 MAE 280 B 138 Maurício de Oliveira
16 % LMI solution % declare variables X = sdpvar(n,n, symmetric ); Z = sdpvar(m,m, symmetric ); L = sdpvar(m,n); % declare LMIs LMI1 = [At*X+X*At +But*L+L *But, X*Cz +L *Dzu ;... Cz*X+Dzu*L, -eye(m+n)]; LMI2 = [Z, Bwt ; Bwt, X]; LMI = set(lmi1 < 0) + set(lmi2 > 0); options = sdpsettings( solver, sedumi ); solution =solvesdp(lmi,trace(w*z),options) SeDuMi 1.1R3 by AdvOL, 2006 and Jos F. Sturm, Alg = 2: xz-corrector, theta = 0.250, beta = eqs m = 21, order n = 17, dim = 137, blocks = 3 nnz(a) = , nnz(ada) = 393, nnz(l) = 207 it : b*y gap delta rate t/tp* t/td* feas cg cg prec 0 : 5.86E : -1.21E E E+00 2 : -3.20E E E-01 3 : -5.95E E E-01 4 : -3.86E E E-01 5 : -5.46E E E-01 6 : -7.75E E E-01 7 : -1.09E E E-03 8 : -1.38E E E-04 9 : -1.74E E E : -2.02E E E : -2.38E E E : -2.60E E E : -2.78E E E : -2.84E E E : -2.86E E E : -2.86E E E : -2.86E E E-11 iter seconds digits c*x b*y Inf e e-02 Ax-b = 4.4e-11, [Ay-c]_+ = 2.0E-11, x = 1.6e+02, y = 2.9e-01 Detailed timing (sec) Pre IPM Post 1.000E E E-02 Max-norms: b = e-01, c = 1, Cholesky add =0, skip = 0, L.L = solution = yalmiptime: e-02 solvertime: e-01 info: No problems detected (SeDuMi-1.1) problem: 0 dimacs: [3.9717e e e e-10] MAE 280 B 139 Maurício de Oliveira
17 % Construct K and check closed loop stability Klqr Klqr = e e e e e e e e+01 K = double(l) / double(x) K = e e e e e e e e+01 Acl = At + But*K; eig(acl) ans = e e-02i e e-02i e e-02i e e-02i diary off MAE 280 B 140 Maurício de Oliveira
18 Let s solve it again with the second form. % MAE 280 B - Linear Control Design % Mauricio de Oliveira % % LMI LQR Control % m = 100; % 100 kg r = 300E3; % 300 km R = 6.37E6; % ˆ3 km G = 6.673E-11; % N mˆ2/kgˆ2 M = 5.98E24; % ˆ24 kg k = G * M; % gravitational force constant w = sqrt(k/((r+r)ˆ3)); % angular velocity (rad/s) v = w * (R + r); % "ground" velocity (m/s) % linearized system matrices A = [ 1 0; 0 1; 3*wˆ2 2*(r+R)*w; -2*w/(r+R) 0]; Bu = [; ; 1/m 0; 0 1/(m*r)]; Bw = [; ; 1/m 0; 0 1/(m*r)]; % noise variances W = 0.1 * eye(2) W = e e-01 % scale = diag([1 r 1 r]) = % similarity transformation At = * A / At = e e e e e-04 0 But = * Bu But = e e-02 Bwt = * Bw Bwt = e e-02 MAE 280 B 141 Maurício de Oliveira
19 % optimal LQR state feedback control (using u1 and u2) n = size(at,1); m = size(but,2); % LQR solution Cz = [eye(n); zeros(m,n)]; Dzu = [zeros(n,m); eye(m)]; [K,X,S] = lqr(at, But, Cz * Cz, Dzu * Dzu); Klqr = - K Klqr = e e e e e e e e+01 % LMI solution (second form) E = null([but Dzu ]); % declare variables X = sdpvar(n,n, symmetric ); Z = sdpvar(m,m, symmetric ); % declare LMIs LMI1 = E *[At*X+X*At, X*Cz ;... Cz*X, -eye(m+n)]*e; LMI2 = [Z, Bwt ; Bwt, X]; LMI = set(lmi1 < 0) + set(lmi2 > 0); options = sdpsettings( solver, sedumi ); solution =solvesdp(lmi,trace(w*z),options) SeDuMi 1.1R3 by AdvOL, 2006 and Jos F. Sturm, Alg = 2: xz-corrector, theta = 0.250, beta = eqs m = 13, order n = 15, dim = 101, blocks = 3 nnz(a) = , nnz(ada) = 169, nnz(l) = 91 it : b*y gap delta rate t/tp* t/td* feas cg cg prec 0 : 4.81E : -1.24E E E+00 2 : -2.72E E E-01 3 : -4.83E E E-01 4 : -3.67E E E-01 5 : -5.44E E E-01 6 : -7.84E E E-02 7 : -1.09E E E-04 8 : -1.39E E E-04 9 : -1.76E E E : -2.05E E E : -2.42E E E : -2.63E E E : -2.79E E E : -2.85E E E : -2.86E E E : -2.86E E E-10 iter seconds digits c*x b*y MAE 280 B 142 Maurício de Oliveira
20 Inf e e-02 Ax-b = 6.5e-10, [Ay-c]_+ = 3.1E-10, x = 1.6e+02, y = 2.9e-01 Detailed timing (sec) Pre IPM Post 1.400E E E-02 Max-norms: b = e-01, c = 1, Cholesky add =0, skip = 0, L.L = solution = yalmiptime: e-01 solvertime: e-01 info: No problems detected (SeDuMi-1.1) problem: 0 dimacs: [5.9141e e e e-09] % Construct K and check closed loop stability Y = inv(double(x)); QQ = [Y*At+At *Y, Cz ;... Cz, -eye(m+n)]; BB = [But *Y Dzu ]; CC = [eye(n) zeros(n,m+n)]; mu = 1e8 mu = max(eig(qq - mu*bb *BB)) ans = e+10 max(eig(qq - mu*cc *CC)) ans = e+00 Phi = QQ - mu*bb *BB - mu*cc *CC; XX = muˆ2 * CC * inv(phi) * BB ; % Check closed loop stability Klqr Klqr = e e e e e e e e+01 Klmi = XX Klmi = e e e e e e e e+01 Acl = At + But*Klmi; eig(acl) ans = e e-02i e e-02i e e-02i e e-02i diary off MAE 280 B 143 Maurício de Oliveira
8 A First Glimpse on Design with LMIs
8 A First Glimpse on Design with LMIs 8.1 Conceptual Design Problem Given a linear time invariant system design a linear time invariant controller or filter so as to guarantee some closed loop indices
More information4 Optimal State Estimation
4 Optimal State Estimation Consider the LTI system ẋ(t) = Ax(t)+B w w(t), y(t) = C y x(t)+d yw w(t), z(t) = C z x(t) Problem: Compute (Â,F) such that the output of the state estimator ˆx(t) = ˆx(t) Fy(t),
More information1 Controllability and Observability
1 Controllability and Observability 1.1 Linear Time-Invariant (LTI) Systems State-space: Dimensions: Notation Transfer function: ẋ = Ax+Bu, x() = x, y = Cx+Du. x R n, u R m, y R p. Note that H(s) is always
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 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 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 informationLecture 10: Linear Matrix Inequalities Dr.-Ing. Sudchai Boonto
Dr-Ing Sudchai Boonto Department of Control System and Instrumentation Engineering King Mongkuts Unniversity of Technology Thonburi Thailand Linear Matrix Inequalities A linear matrix inequality (LMI)
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 informationAppendix A Solving Linear Matrix Inequality (LMI) Problems
Appendix A Solving Linear Matrix Inequality (LMI) Problems In this section, we present a brief introduction about linear matrix inequalities which have been used extensively to solve the FDI problems described
More informationLOW ORDER H CONTROLLER DESIGN: AN LMI APPROACH
LOW ORDER H CONROLLER DESIGN: AN LMI APPROACH Guisheng Zhai, Shinichi Murao, Naoki Koyama, Masaharu Yoshida Faculty of Systems Engineering, Wakayama University, Wakayama 640-8510, Japan Email: zhai@sys.wakayama-u.ac.jp
More informationEE363 homework 8 solutions
EE363 Prof. S. Boyd EE363 homework 8 solutions 1. Lyapunov condition for passivity. The system described by ẋ = f(x, u), y = g(x), x() =, with u(t), y(t) R m, is said to be passive if t u(τ) T y(τ) dτ
More information2 The Linear Quadratic Regulator (LQR)
2 The Linear Quadratic Regulator (LQR) Problem: Compute a state feedback controller u(t) = Kx(t) that stabilizes the closed loop system and minimizes J := 0 x(t) T Qx(t)+u(t) T Ru(t)dt where x and u are
More informationH 2 Optimal State Feedback Control Synthesis. Raktim Bhattacharya Aerospace Engineering, Texas A&M University
H 2 Optimal State Feedback Control Synthesis Raktim Bhattacharya Aerospace Engineering, Texas A&M University Motivation Motivation w(t) u(t) G K y(t) z(t) w(t) are exogenous signals reference, process
More informationON STABILITY TESTS FOR LINEAR SYSTEMS
Copyright 22 IFAC 15th Triennial World Congress Barcelona Spain ON STABILITY TESTS FOR LINEAR SYSTEMS Maurício C. de Oliveira 1 Robert E. Skelton Department of Telematics School of Electrical and Computer
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 information7.1 Linear Systems Stability Consider the Continuous-Time (CT) Linear Time-Invariant (LTI) system
7 Stability 7.1 Linear Systems Stability Consider the Continuous-Time (CT) Linear Time-Invariant (LTI) system ẋ(t) = A x(t), x(0) = x 0, A R n n, x 0 R n. (14) The origin x = 0 is a globally asymptotically
More informationPoluqenie rezulьtatov v vide. line inyx matreqnyh neravenstv dl robastnosti. Tradicionna Xkola Upravlenie, Informaci i Optimizaci
Poluqenie rezulьtatov v vide line inyx matreqnyh neravenstv dl robastnosti Dmitri i Жanoviq Konovalov Dimitri i Poselь Dmitry Peaucelle Dimitri Peaucelle Tradicionna Xkola Upravlenie, Informaci i Optimizaci
More informationStatic Output Feedback Stabilisation with H Performance for a Class of Plants
Static Output Feedback Stabilisation with H Performance for a Class of Plants E. Prempain and I. Postlethwaite Control and Instrumentation Research, Department of Engineering, University of Leicester,
More informationFrom Convex Optimization to Linear Matrix Inequalities
Dep. of Information Engineering University of Pisa (Italy) From Convex Optimization to Linear Matrix Inequalities eng. Sergio Grammatico grammatico.sergio@gmail.com Class of Identification of Uncertain
More informationH 2 Suboptimal Estimation and Control for Nonnegative
Proceedings of the 2007 American Control Conference Marriott Marquis Hotel at Times Square New York City, USA, July 11-13, 2007 FrC20.3 H 2 Suboptimal Estimation and Control for Nonnegative Dynamical Systems
More informationState Estimation with Finite Signal-to-Noise Models
State Estimation with Finite Signal-to-Noise Models Weiwei Li and Robert E. Skelton Department of Mechanical and Aerospace Engineering University of California San Diego, La Jolla, CA 9293-411 wwli@mechanics.ucsd.edu
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 informationRobust Multi-Objective Control for Linear Systems
Robust Multi-Objective Control for Linear Systems Elements of theory and ROMULOC toolbox Dimitri PEAUCELLE & Denis ARZELIER LAAS-CNRS, Toulouse, FRANCE Part of the OLOCEP project (includes GloptiPoly)
More informationModern Optimal Control
Modern Optimal Control Matthew M. Peet Arizona State University Lecture 22: H 2, LQG and LGR Conclusion To solve the H -optimal state-feedback problem, we solve min γ such that γ,x 1,Y 1,A n,b n,c n,d
More informationSemidefinite Programming Duality and Linear Time-invariant Systems
Semidefinite Programming Duality and Linear Time-invariant Systems Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University 2 July 2004 Workshop on Linear Matrix Inequalities in Control LAAS-CNRS,
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 14: LMIs for Robust Control in the LF Framework ypes of Uncertainty In this Lecture, we will cover Unstructured,
More informationRobust and Optimal Control, Spring 2015
Robust and Optimal Control, Spring 2015 Instructor: Prof. Masayuki Fujita (S5-303B) D. Linear Matrix Inequality D.1 Convex Optimization D.2 Linear Matrix Inequality(LMI) D.3 Control Design and LMI Formulation
More informationECE 680 Linear Matrix Inequalities
ECE 680 Linear Matrix Inequalities Stan Żak School of Electrical and Computer Engineering Purdue University zak@purdue.edu October 11, 2017 Stan Żak (Purdue University) ECE 680 Linear Matrix Inequalities
More informationLMI based output-feedback controllers: γ-optimal versus linear quadratic.
Proceedings of the 17th World Congress he International Federation of Automatic Control Seoul Korea July 6-11 28 LMI based output-feedback controllers: γ-optimal versus linear quadratic. Dmitry V. Balandin
More informationSYNTHESIS OF ROBUST DISCRETE-TIME SYSTEMS BASED ON COMPARISON WITH STOCHASTIC MODEL 1. P. V. Pakshin, S. G. Soloviev
SYNTHESIS OF ROBUST DISCRETE-TIME SYSTEMS BASED ON COMPARISON WITH STOCHASTIC MODEL 1 P. V. Pakshin, S. G. Soloviev Nizhny Novgorod State Technical University at Arzamas, 19, Kalinina ul., Arzamas, 607227,
More informationIntroduction to Nonlinear Control Lecture # 4 Passivity
p. 1/6 Introduction to Nonlinear Control Lecture # 4 Passivity È p. 2/6 Memoryless Functions ¹ y È Ý Ù È È È È u (b) µ power inflow = uy Resistor is passive if uy 0 p. 3/6 y y y u u u (a) (b) (c) Passive
More informationIntroduction to Linear Matrix Inequalities (LMIs)
ECE 680 Fall 2017 Introduction to Linear Matrix Inequalities (LMIs) by Stanislaw H. Żak and Guisheng Zhai October 18, 2017 Henry Ford offered customers the Model T Ford in any color the customer wants,
More informationRobust Control 5 Nominal Controller Design Continued
Robust Control 5 Nominal Controller Design Continued Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 4/14/2003 Outline he LQR Problem A Generalization to LQR Min-Max
More informationRiccati Equations and Inequalities in Robust Control
Riccati Equations and Inequalities in Robust Control Lianhao Yin Gabriel Ingesson Martin Karlsson Optimal Control LP4 2014 June 10, 2014 Lianhao Yin Gabriel Ingesson Martin Karlsson (LTH) H control problem
More informationAn LMI Optimization Approach for Structured Linear Controllers
An LMI Optimization Approach for Structured Linear Controllers Jeongheon Han* and Robert E. Skelton Structural Systems and Control Laboratory Department of Mechanical & Aerospace Engineering University
More informationRobust Control 9 Design Summary & Examples
Robust Control 9 Design Summary & Examples Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University 5/6/003 Outline he H Problem Solution via γ-iteration Robust stability via
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 informationLinear Matrix Inequality (LMI)
Linear Matrix Inequality (LMI) A linear matrix inequality is an expression of the form where F (x) F 0 + x 1 F 1 + + x m F m > 0 (1) x = (x 1,, x m ) R m, F 0,, F m are real symmetric matrices, and the
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 informationPreliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012
Instructions Preliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012 The exam consists of four problems, each having multiple parts. You should attempt to solve all four problems. 1.
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 informationLecture 15: H Control Synthesis
c A. Shiriaev/L. Freidovich. March 12, 2010. Optimal Control for Linear Systems: Lecture 15 p. 1/14 Lecture 15: H Control Synthesis Example c A. Shiriaev/L. Freidovich. March 12, 2010. Optimal Control
More information1 Some Facts on Symmetric Matrices
1 Some Facts on Symmetric Matrices Definition: Matrix A is symmetric if A = A T. Theorem: Any symmetric matrix 1) has only real eigenvalues; 2) is always iagonalizable; 3) has orthogonal eigenvectors.
More informationEL2520 Control Theory and Practice
EL2520 Control Theory and Practice Lecture 8: Linear quadratic control Mikael Johansson School of Electrical Engineering KTH, Stockholm, Sweden Linear quadratic control Allows to compute the controller
More informationRobust Anti-Windup Controller Synthesis: A Mixed H 2 /H Setting
Robust Anti-Windup Controller Synthesis: A Mixed H /H Setting ADDISON RIOS-BOLIVAR Departamento de Sistemas de Control Universidad de Los Andes Av. ulio Febres, Mérida 511 VENEZUELA SOLBEN GODOY Postgrado
More informationA New Strategy to the Multi-Objective Control of Linear Systems
Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference 25 Seville, Spain, December 12-15, 25 TuC8.6 A New Strategy to the Multi-Objective Control of Linear
More informationResearch Article An Equivalent LMI Representation of Bounded Real Lemma for Continuous-Time Systems
Hindawi Publishing Corporation Journal of Inequalities and Applications Volume 28, Article ID 67295, 8 pages doi:1.1155/28/67295 Research Article An Equivalent LMI Representation of Bounded Real Lemma
More informationEE363 homework 7 solutions
EE363 Prof. S. Boyd EE363 homework 7 solutions 1. Gain margin for a linear quadratic regulator. Let K be the optimal state feedback gain for the LQR problem with system ẋ = Ax + Bu, state cost matrix Q,
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 informationLinear Matrix Inequality Approach
1 MEM800-007 Chapter 4a Linear Matrix Inequality Approach Reference: Linear matrix inequalities in system and control theory, Stephen Boyd et al. SIAM 1994. Linear Matrix Inequality (LMI) approach have
More informationH 2 and H 1 cost estimates for time-invariant uncertain
INT. J. CONTROL, 00, VOL. 75, NO. 9, ±79 Extended H and H systems norm characterizations and controller parametrizations for discrete-time M. C. DE OLIVEIRAy*, J. C. GEROMELy and J. BERNUSSOUz This paper
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 informationController Coefficient Truncation Using Lyapunov Performance Certificate
Controller Coefficient Truncation Using Lyapunov Performance Certificate Joëlle Skaf Stephen Boyd Information Systems Laboratory Electrical Engineering Department Stanford University European Control Conference,
More informationIntroduction to linear matrix inequalities Wojciech Paszke
Introduction to linear matrix inequalities Wojciech Paszke Institute of Control and Computation Engineering, University of Zielona Góra, Poland e-mail: W.Paszke@issi.uz.zgora.pl Outline Introduction to
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 informationOutput Stabilization of Time-Varying Input Delay System using Interval Observer Technique
Output Stabilization of Time-Varying Input Delay System using Interval Observer Technique Andrey Polyakov a, Denis Efimov a, Wilfrid Perruquetti a,b and Jean-Pierre Richard a,b a - NON-A, INRIA Lille Nord-Europe
More informationTopic # Feedback Control Systems
Topic #17 16.31 Feedback Control Systems Deterministic LQR Optimal control and the Riccati equation Weight Selection Fall 2007 16.31 17 1 Linear Quadratic Regulator (LQR) Have seen the solutions to the
More informationEconomic sensor/actuator selection and its application to flexible structure control
To appear in Proc. SPIE Int. Soc. Opt. Eng. 2004 Economic sensor/actuator selection and its application to flexible structure control Robert E. Skelton a and Faming Li a a Department of Mechanical and
More informationConvex Optimization Problems. Prof. Daniel P. Palomar
Conve Optimization Problems Prof. Daniel P. Palomar The Hong Kong University of Science and Technology (HKUST) MAFS6010R- Portfolio Optimization with R MSc in Financial Mathematics Fall 2018-19, HKUST,
More informationCOVARIANCE CALCULATION FOR FLOATING POINT STATE SPACE REALIZATIONS
opyright 2002 IFA 15th riennial World ongress, Barcelona, Spain OVARIANE ALULAION FOR FLOAING POIN SAE SPAE REALIZAIONS Sangho Ko,1 Robert R Bitmead,1 Department of Mechanical and Aerospace Engineering,
More informationAn Introduction to Linear Matrix Inequalities. Raktim Bhattacharya Aerospace Engineering, Texas A&M University
An Introduction to Linear Matrix Inequalities Raktim Bhattacharya Aerospace Engineering, Texas A&M University Linear Matrix Inequalities What are they? Inequalities involving matrix variables Matrix variables
More informationLinear Matrix Inequalities in Robust Control. Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University MTNS 2002
Linear Matrix Inequalities in Robust Control Venkataramanan (Ragu) Balakrishnan School of ECE, Purdue University MTNS 2002 Objective A brief introduction to LMI techniques for Robust Control Emphasis on
More information1 / 21 Perturbation of system dynamics and the covariance completion problem Armin Zare Joint work with: Mihailo R. Jovanović Tryphon T. Georgiou 55th
1 / 21 Perturbation of system dynamics and the covariance completion problem Armin Zare Joint work with: Mihailo R. Jovanović Tryphon T. Georgiou 55th IEEE Conference on Decision and Control, Las Vegas,
More informationANALYSIS OF DISCRETE-TIME H2 GUARANTEED COST PERFORMANCE
ANALYSIS OF DISCREE-IME H2 GUARANEED COS PERFORMANCE Richard Conway Computer Mechanics Laboratory Department of Mechanical Engineering University of California at Berkeley Berkeley, California 94720-740
More informationNonlinear Control Lecture 7: Passivity
Nonlinear Control Lecture 7: Passivity Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2011 Farzaneh Abdollahi Nonlinear Control Lecture 7 1/26 Passivity
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 informationModern Optimal Control
Modern Optimal Control Matthew M. Peet Arizona State University Lecture 21: Optimal Output Feedback Control connection is called the (lower) star-product of P and Optimal Output Feedback ansformation (LFT).
More informationMassachusetts Institute of Technology. Department of Electrical Engineering and Computer Science : MULTIVARIABLE CONTROL SYSTEMS by A.
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Q-Parameterization 1 This lecture introduces the so-called
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 informationROBUST ANALYSIS WITH LINEAR MATRIX INEQUALITIES AND POLYNOMIAL MATRICES. Didier HENRION henrion
GRADUATE COURSE ON POLYNOMIAL METHODS FOR ROBUST CONTROL PART IV.1 ROBUST ANALYSIS WITH LINEAR MATRIX INEQUALITIES AND POLYNOMIAL MATRICES Didier HENRION www.laas.fr/ henrion henrion@laas.fr Airbus assembly
More informationLMIs for Observability and Observer Design
LMIs for Observability and Observer Design Matthew M. Peet Arizona State University Lecture 06: LMIs for Observability and Observer Design Observability Consider a system with no input: ẋ(t) = Ax(t), x(0)
More informationOutline. Linear Matrix Inequalities in Control. Outline. System Interconnection. j _jst. ]Bt Bjj. Generalized plant framework
Outline Linear Matrix Inequalities in Control Carsten Scherer and Siep Weiland 7th Elgersburg School on Mathematical Systems heory Class 3 1 Single-Objective Synthesis Setup State-Feedback Output-Feedback
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 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 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 informationStabilization of persistently excited linear systems
Stabilization of persistently excited linear systems Yacine Chitour Laboratoire des signaux et systèmes & Université Paris-Sud, Orsay Exposé LJLL Paris, 28/9/2012 Stabilization & intermittent control Consider
More informationSemidefinite Programming Basics and Applications
Semidefinite Programming Basics and Applications Ray Pörn, principal lecturer Åbo Akademi University Novia University of Applied Sciences Content What is semidefinite programming (SDP)? How to represent
More informationAlgebra C Numerical Linear Algebra Sample Exam Problems
Algebra C Numerical Linear Algebra Sample Exam Problems Notation. Denote by V a finite-dimensional Hilbert space with inner product (, ) and corresponding norm. The abbreviation SPD is used for symmetric
More informationFRTN10 Exercise 12. Synthesis by Convex Optimization
FRTN Exercise 2. 2. We want to design a controller C for the stable SISO process P as shown in Figre 2. sing the Yola parametrization and convex optimization. To do this, the control loop mst first be
More informationLecture 5 Linear Quadratic Stochastic Control
EE363 Winter 2008-09 Lecture 5 Linear Quadratic Stochastic Control linear-quadratic stochastic control problem solution via dynamic programming 5 1 Linear stochastic system linear dynamical system, over
More informationSTABILITY AND STABILIZATION OF A CLASS OF NONLINEAR SYSTEMS WITH SATURATING ACTUATORS. Eugênio B. Castelan,1 Sophie Tarbouriech Isabelle Queinnec
STABILITY AND STABILIZATION OF A CLASS OF NONLINEAR SYSTEMS WITH SATURATING ACTUATORS Eugênio B. Castelan,1 Sophie Tarbouriech Isabelle Queinnec DAS-CTC-UFSC P.O. Box 476, 88040-900 Florianópolis, SC,
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 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 informationConstrained Nonlinear Optimization Algorithms
Department of Industrial Engineering and Management Sciences Northwestern University waechter@iems.northwestern.edu Institute for Mathematics and its Applications University of Minnesota August 4, 2016
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 informationQuadratic Stability of Dynamical Systems. Raktim Bhattacharya Aerospace Engineering, Texas A&M University
.. Quadratic Stability of Dynamical Systems Raktim Bhattacharya Aerospace Engineering, Texas A&M University Quadratic Lyapunov Functions Quadratic Stability Dynamical system is quadratically stable if
More informationSLIDING SURFACE MATCHED CONDITION IN SLIDING MODE CONTROL
Asian Journal of Control, Vol. 9, No. 3, pp. 0-0, September 007 1 SLIDING SURFACE MACHED CONDIION IN SLIDING MODE CONROL Ji Xiang, Hongye Su, Jian Chu, and Wei Wei -Brief Paper- ABSRAC In this paper, the
More informationTaylor expansions for the HJB equation associated with a bilinear control problem
Taylor expansions for the HJB equation associated with a bilinear control problem Tobias Breiten, Karl Kunisch and Laurent Pfeiffer University of Graz, Austria Rome, June 217 Motivation dragged Brownian
More information6.241 Dynamic Systems and Control
6.41 Dynamic Systems and Control Lecture 5: H Synthesis Emilio Frazzoli Aeronautics and Astronautics Massachusetts Institute of Technology May 11, 011 E. Frazzoli (MIT) Lecture 5: H Synthesis May 11, 011
More informationLinear-Quadratic Optimal Control: Full-State Feedback
Chapter 4 Linear-Quadratic Optimal Control: Full-State Feedback 1 Linear quadratic optimization is a basic method for designing controllers for linear (and often nonlinear) dynamical systems and is actually
More informationROBUST STABILITY TEST FOR UNCERTAIN DISCRETE-TIME SYSTEMS: A DESCRIPTOR SYSTEM APPROACH
Latin American Applied Research 41: 359-364(211) ROBUS SABILIY ES FOR UNCERAIN DISCREE-IME SYSEMS: A DESCRIPOR SYSEM APPROACH W. ZHANG,, H. SU, Y. LIANG, and Z. HAN Engineering raining Center, Shanghai
More informationControl, Stabilization and Numerics for Partial Differential Equations
Paris-Sud, Orsay, December 06 Control, Stabilization and Numerics for Partial Differential Equations Enrique Zuazua Universidad Autónoma 28049 Madrid, Spain enrique.zuazua@uam.es http://www.uam.es/enrique.zuazua
More informationLMI Methods in Optimal and Robust Control
LMI Methods in Optimal and Robust Control Matthew M. Peet Arizona State University Lecture 4: LMIs for State-Space Internal Stability Solving the Equations Find the output given the input State-Space:
More informationFrom the multivariable Nyquist plot, if 1/k < or 1/k > 5.36, or 0.19 < k < 2.7, then the closed loop system is stable.
1. (a) Multivariable Nyquist Plot (diagonal uncertainty case) dash=neg freq solid=pos freq 2.5 2 1.5 1.5 Im -.5-1 -1.5-2 -2.5 1 2 3 4 5 Re From the multivariable Nyquist plot, if 1/k < -.37 or 1/k > 5.36,
More informationInput/output delay approach to robust sampled-data H control
Systems & Control Letters 54 (5) 71 8 www.elsevier.com/locate/sysconle Input/output delay approach to robust sampled-data H control E. Fridman, U. Shaked, V. Suplin Department of Electrical Engineering-Systems,
More informationProblem Set 8 Solutions 1
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6245: MULTIVARIABLE CONTROL SYSTEMS by A Megretski Problem Set 8 Solutions 1 The problem set deals with H-Infinity
More informationLinear Matrix Inequalities vs Convex Sets
Linear Matrix Inequalities vs Convex Sets with admiration and friendship for Eduardo Harry Dym Math Dept Weitzman Inst. Damon Hay Math Dept Florida Texas Igor Klep Math Dept Everywhere in Solvenia Scott
More informationOPTIMAL CONTROL SYSTEMS
SYSTEMS MIN-MAX Harry G. Kwatny Department of Mechanical Engineering & Mechanics Drexel University OUTLINE MIN-MAX CONTROL Problem Definition HJB Equation Example GAME THEORY Differential Games Isaacs
More informationTo appear in IEEE Trans. on Automatic Control Revised 12/31/97. Output Feedback
o appear in IEEE rans. on Automatic Control Revised 12/31/97 he Design of Strictly Positive Real Systems Using Constant Output Feedback C.-H. Huang P.A. Ioannou y J. Maroulas z M.G. Safonov x Abstract
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 information