Lecture 2 – Linear Systems

1 EE392m - Spring 2005 GorinevskyControl Engineering2-1 Lecture 2 Linear SystemsThis Lecture : EE263 material recap + some controls motivation Continuous time (physics) Linear state space model Transfer functions Black-box models; frequency domain analysis LinearizationEE392m - Spring 2005 GorinevskyControl Engineering2-2 Modeling and Analysis This Lecture considers Linear models. More detail on modeling in Lecture 7 Simulation: computing state evolution and output signal Stability: does the solution diverge after some time? Approximate Linear modelsSystem Model[internal states]Input signal(control)Output signal(observation)

2 State evolutionEE392m - Spring 2005 GorinevskyControl Engineering2-3 Linear Models Model is a mathematical representations of a system Models allow simulating the system Models can be used for conceptual analysis Models are never exact Linear models Have simple structure Can be analyzed using powerful mathematical tools Can be matched against real data using known procedures Many simple physics models are Linear They are just models, not the real Systems EE392m - Spring 2005 GorinevskyControl Engineering2-4 State space model Generic state space model is described by ODEs , physics-based system model state vector x observation vector y control vectoru),(),,(txgytuxfdtdx== state evolution observationExample: F16 Longitudinal Model.

3 180081 8201095210152910021105421705802328281093 112342 + == + =+ + =&&&&zx e Veu = Control and Simulationby Stevens and Lewisx1- velocity V [ft/sec] x2- angle of attack [rad] x3- pitch angle [rad] x4- pitch rate q [rad/sec] e- elevator deflection [deg]EE392m - Spring 2005 GorinevskyControl Engineering2-5 Linear state space model Linear Time Invariant (LTI) state space model: Can be integrated analytically or numerically (simulation) Can be well analyzed: stability, responseCxyBuAxdtdx=+= state evolution observationsExample: F16 Longitudinal Model []xyuxdtdx = + = - Spring 2005 GorinevskyControl Engineering2-6 Integrating Linear autonomous system Matrix exponential Can be computed in Matlab as expm(A) The definition corresponds to integrating the ODE by Euler methodAxdtdx=()ttttAIAt +=/0lim)exp(())()()()(TxtAIttAxTxtTx += + +Example: >>A = [ ; 0 ;0 0 0 1; 0 ].

4 >> expm(A)ans = ())()exp()()(TxtAnTxtAItnTxn += +EE392m - Spring 2005 GorinevskyControl Engineering2-7 Integrating Linear autonomous system Initial condition response0)0(,xxAxdtdx==0)exp()(xAttx=Exa mple: %Take A, B, C from the F16 example>> k = >> G = A + k*B*C;>> x0 = [ ; ;0; ];>> for j = 1:length(t); x(:,j)=expm(G*t(j))*x0; end; MEEE392m - Spring 2005 GorinevskyControl Engineering2-8 Eigenvalues and Stability Consider eigenvalues of the matrix A Suppose Ahas all different and nonzero eigenvalues, then The system solution is exponentially stable if If Ahas eigenvalues with multiplicity more than 1, things are a bit more complicated: Jordan blocks, polynomials in t Still the condition of exponentially stability is 0)det().

5 Eig(}{= =AIAj 1}diag{ =VVAj 0Re<j 1}diag{)exp( =VeVAttj Example: % take A from % the F16 example>> eig(A)ans = + - <j EE392m - Spring 2005 GorinevskyControl Engineering2-9 Input-output models Black-box models describe system Pas an operatorPxuinput signalyoutput signalinternal state (hidden) Historically (50 years ago) Black-box models EE State-space models ME, AAEE392m - Spring 2005 GorinevskyControl Engineering2-10 Linear system (input-output) Linearity Linear Time-Invariant Systems - LTI )()(11 yuP)()(TyTuP )()()()(2121 + + byaybuauP)()

6 (22 yuP PutytEE392m - Spring 2005 GorinevskyControl Engineering2-11 Convolution representation Convolution integral Impulse response Step response: u = 1 for t > 0 =tduthty )()()()()()()(thtyttu= = = =ttdhdthtg00)()()( )()(tgdtdth=uhy*=signal processing notationEE392m - Spring 2005 GorinevskyControl Engineering2-12 Impulse Response for State Space Model Impulse response for the state x system impulse response {1)0()(tuBtAxtx + BtuBtx= )(BAttx)exp()(=CxyBuAxdtdx=+= state evolution observationExample: >> A = [ 0; 0;0 0];>> B = [ ; ; 0];>> C = [0, 0, 1].)}

7 Exp()(=EE392m - Spring 2005 GorinevskyControl Engineering2-13 Formal transfer function Rational transfer function = IIR (Infinite Impulse Response) model Broad class of input-output Linear models Differentiation operator Formal transfer function rational function of s For a causal system m nsdtd ubdtudbdtudbyadtydadtydannnnnmmmmm111211 1121+ + +++=+++KKusDsNusHy = =)()()(11)(++++=mmmasasasNK11)(++++=nnnb sbsbsDKEE392m - Spring 2005 GorinevskyControl Engineering2-14 Poles, Impulse Response ExpandThen Quasi-polynomial impulse response see a textbook Example:usDsNy =)()(0)(=sN0)(=sD- zeros- polesKMKM pspssD)()()(11 =KupssNupssNyKMKKM ++ =)()()()(111K()()tpMKMKtpMMKKK ectcectcth1,11,1,111,1111)(+ + ++++++=KKKudtyd=22usy =2100)(xtvth+=Transfer function:Impulse response.)

8 EE392m - Spring 2005 GorinevskyControl Engineering2-15 Transfer Function for State Space Characteristic polynomial Poles are the same as eigenvalues of the state-space matrix A For stability we need Re pk= Re k< 0()()BAsICsHuBAsIy11)( = =Poles eigenvalues()0det= AsI()0det)(= =AsIsNCxyBuAxsx=+= Formal transfer function for a state space model sdtd EE392m - Spring 2005 GorinevskyControl Engineering2-16 Laplace transform Laplace integral transform: Laplace transform of the convolution integral yields Transfer function: function of complex variable s analytical in a right half-plane Res a for a stable system a 0 for an IIR model =0)()(dtethsHst)( )()( susHsy=)( sxsdtdx dtetxsxtxst = 0)()( )(Re sImags)(sH)()()(sDsNsH=EE392m - Spring 2005 GorinevskyControl Engineering2-17 Frequency decomposition Sinusoids are eigenfunctions of an LTI system .

9 LTIP lant Frequency domain analysisusHy)(= = tikktikkkeiHuyeu )(tie tieiH )(()tititieiedtdes = EE392m - Spring 2005 GorinevskyControl Engineering2-18 Frequency domain description = = deuiHydeuutiyti43421)(~)(~)(21)(~21)(~ uetiuPacket of sinusoids)(~ yetiPacket of sinusoids)( iHy Frequency domain analysis)( )()()(~ iudetudetuuisstti= === Fourier transform numerical analysis Laplace transform complex analysis u(t) = 0, for t < 0 EE392m - Spring 2005 GorinevskyControl Engineering2-19 Continuous Systems in frequency domain Fourier transform Inverse Fourier transform I/O impulse response model Transfer function system frequency response)(~)()(~)()()()()()(~21)()()(~0 uiHydtethsHdtuthtydextxdtetxxstttiti== === ],[],[ EE392m - Spring 2005 GorinevskyControl Engineering2-20 Frequency domain description Bode plots.

10 TitieiHyeu )(== )( =ssH Bode DiagramFrequency (rad/sec) Phase (deg) Magnitude (dB) -505101510-210-1100-180-135-90-450 |H|is often measured in dB [dB] = 20 log10M)(arg)()()( iHiHM==EE392m - Spring 2005 GorinevskyControl Engineering2-21 Model Approximation Model structure physics, computational Determine parameters from data Step/impulse responses are close the input/output models are close Example fit step response Linearization of nonlinear model EE392m - Spring 2005 GorinevskyControl Engineering2-22 Black-box model from data Linear black-box model can be determined from the data, , step response data, or frequency response Example problem.