This preview shows pages 1–3. Sign up to view the full content.
This preview has intentionally blurred sections. Sign up to view the full version.View Full Document
Unformatted text preview: Ep = diag( [ Ep(3,3), Ep(2,2), Ep(1,1)] ) E Ap Qp * Ep * inv(Qp) Q [Qc,Ec]=eig(Ac); Ac Qc * Ec * inv(Qc) Q [Qd,Ed]=eig(Ad); Ad Qd * Ed * inv(Qd) Q % transformation % T = [ 0 1 0; 0 0 0.6; 1 -1 0]; T T * Ap * inv(T) Ac A T * Bp Bc B Cp * inv(T) Cc % eigenvalue & eigenvector by null space calculation % q1 = [ -6/5; -3/5; 1 ]; q2 = [ -3/2; -3/10; 1]; q3 = [ -2; -1/5; 1 ]; q Q = [ q1 q2 q3 ]; Q Ad = inv(Q)*Ap*(Q) Td = inv(Q) T Bd = [ 3; -6; 3 ]; Cd = [ 1 1 1 ]; Dd = 0; D bigCp = [ Bp Ap*Bp Ap^2*Bp]; bigCc = [ Bc Ac*Bc Ac^2*Bc]; b T = bigCc*inv(bigCp) T T * Ap * inv(T) Ac A T * Bp Bc B Cp * inv(T) Cc C...
View Full Document
This note was uploaded on 07/21/2009 for the course EE 901-43400 taught by Professor Feng-lilian during the Fall '08 term at National Taiwan University.
- Fall '08