Difference between revisions of "Linear Systems"
(→Homework) 
(→Errata) 

(42 intermediate revisions by the same user not shown)  
Line 5:  Line 5:  
{{chaptertable left}}  {{chaptertable left}}  
== Textbook Contents ==  == Textbook Contents ==  
−  {{am05pdf  +  {{am05pdfam08linsys28Sep12Linear Systems}} 
−  * 1.  +  * 1. Basic Definitions 
−  * 2.  +  * 2. The Matrix Exponential 
−  * 3.  +  * 3. Input/Output Response 
−  * 4.  +  * 4. Linearization 
−  * 5  +  * 5. Further Reading 
−  +  * Exercises  
−  *  +  
{{chaptertable right}}  {{chaptertable right}}  
Line 22:  Line 21:  
== Supplemental Information ==  == Supplemental Information ==  
* [[#Frequently Asked QuestionsFrequently Asked Questions]]  * [[#Frequently Asked QuestionsFrequently Asked Questions]]  
+  * [[#ErrataErrata]]  
* Wikipedia entries: [http://en.wikipedia.org/wiki/Linear_system linear system], [http://en.wikipedia.org/wiki/LTI_system_theory LTI system theory]  * Wikipedia entries: [http://en.wikipedia.org/wiki/Linear_system linear system], [http://en.wikipedia.org/wiki/LTI_system_theory LTI system theory]  
* [[#Additional InformationAdditional Information]]  * [[#Additional InformationAdditional Information]]  
Line 29:  Line 29:  
This chapter introduces the analysis tools for linear input/output systems.  This chapter introduces the analysis tools for linear input/output systems.  
<ol>  <ol>  
−  <li> <p>A ''linear system'' is one in which the output is jointly linear in the intitial condition for the system and the input to the system. In particular, a linear system has the property that if we apply an input <  +  <li> <p>A ''linear system'' is one in which the output is jointly linear in the intitial condition for the system and the input to the system. In particular, a linear system has the property that if we apply an input <amsmath>u(t) = \alpha u_1(t) + \beta u_2(t)</amsmath> with zero initial condition, the corresponding output will be <amsmath>y(t) = \alpha y_1(t) + \beta y_2(t)</amsmath>, where <amsmath>y_i</amsmath> is the output associated with the input <amsmath>u_i</amsmath>. This propery is called linear ''superposition''.</p></li> 
−  <li><p>A differential equation of the form  +  <li><p><span id=SISO>A differential equation of the form 
+  <center><amsmath>  
+  \begin{aligned}  
+  \dot x &= A x + B u &\quad x &\in R^n, u \in R \\  
+  y &= C x + D u & y &\in R  
+  \end{aligned}  
+  </amsmath></center>  
+  is a ''singleinput, singleoutput'' (SISO) ''linear differential equation''</span>. Its solution can be written in terms of the ''[[matrix exponential]]''  
<center>  <center>  
−  +  <amsmath>  
−  +  
−  +  
−  +  
−  +  
−  +  
−  +  
−  +  
−  +  
−  +  
−  +  
−  +  
−  +  
−  +  
e^{At}=I+At+\frac{1}{2}A^2t^2+\frac{1}{3!}A^3t^3+\cdots=\sum_{k=0}^\infty \frac{1}{k!} A^kt^k.  e^{At}=I+At+\frac{1}{2}A^2t^2+\frac{1}{3!}A^3t^3+\cdots=\sum_{k=0}^\infty \frac{1}{k!} A^kt^k.  
−  </  +  </amsmath> 
</center>  </center>  
The solution to the differential equation is given by the ''convolution equation''  The solution to the differential equation is given by the ''convolution equation''  
−  <center><  +  <center><amsmath> 
y(t) = C e^{At} x(0) + \int_0^t C e^{A(t\tau)}B u(\tau) d\tau + Du(t).  y(t) = C e^{At} x(0) + \int_0^t C e^{A(t\tau)}B u(\tau) d\tau + Du(t).  
−  </  +  </amsmath></center> 
−  </p>  +  </p></li> 
−  <li><p>A linear system <center><  +  <li><p>A linear system <center><amsmath>\dot x = A x</amsmath></center> is [[Dynamic Behavior#stabilityasymptotically stable]] if and only if all eigenvalues of <amsmath>A</amsmath> all have strictly negative real part and is {{unstable}} if any eigenvalue of <amsmath>A</amsmath> has strictly positive real part. For systems with eigenvalues having zero realpart, stability is determined by using the ''[http://en.wikipedia.org/wiki/Jordan_normal_form Jordan normal form]'' associated with the matrix. A system with eigenvalues that have no strictly positive real part is stable if and only if the Jordan block corresponding to each eigenvalue with zero part is a scalar (1x1) block.</p></li> 
−  <li><p>The input/output response of a (stable) linear system contains a transient region portion, which eventually decays to zero, and a steady state portion, which persists over time. Two special responses are the ''step response'', which is the output corresponding to an step input applied at <  +  <li><p>The input/output response of a (stable) linear system contains a transient region portion, which eventually decays to zero, and a steady state portion, which persists over time. Two special responses are the ''step response'', which is the output corresponding to an step input applied at <amsmath>t = 0</amsmath> and the ''frequency response'', which is the response of the system to a sinusoidal input at a given frequency.</p></li> 
−  * The ''steady state value'', <  +  
−  * The ''rise time'', <  +  <li><p> 
−  * The ''overshoot'', <  +  The step response is characterized by the following parameters: 
−  * The ''settling time'', <  +  * The ''steady state value'', <amsmath>y_{ss}</amsmath>, of a step response is the final level of the output, assuming it converges. 
−  </p>  +  * The ''rise time'', <amsmath>T_r</amsmath>, is the amount of time required for the signal to go from 10% of its final value to 90% of its final value. 
+  * The ''overshoot'', <amsmath>M_p</amsmath>, is the percentage of the initial value by which the signal initially rises above the final value.  
+  * The ''settling time'', <amsmath>T_s</amsmath>, is the amount of time required for the signal to stay within 5% of its final value for all future times.  
+  </p></li>  
+  
+  <li><p>The frequency response is given by  
+  <center><amsmath>  
+  y(t) = \underbrace{C e^{At} \Bigl( x(0)(s IA)^{1} B  
+  \Bigr)}_{\text{transient}} +  
+  \underbrace{\Bigl( D+C(s IA)^{1}B \Bigr) e^{s  
+  t}}_{\text{steady state}},  
+  </amsmath></center>  
+  where <amsmath>\cos\omega t = \frac{1}{2} \left(e^{j \omega t} + e^{j \omega t}\right)</amsmath> and <amsmath>s = j \omega</amsmath>. The gain and phase of the frequency response are given by  
+  <center><amsmath>  
+  \text{gain}(\omega) = \frac{A_y}{A_u} = M \qquad  
+  \text{phase}(\omega) = \phi  \psi = \theta.  
+  </amsmath></center>  
+  </p></li>  
<li><p>A nonlinear system of the form  <li><p>A nonlinear system of the form  
+  <center><amsmath>  
+  \begin{aligned}  
+  \dot x &= f(x, u) &\quad x &\in R^n, u \in R \\  
+  y &= h(x, u) & y &\in R  
+  \end{aligned}  
+  </amsmath></center>  
+  is a singleinput, singleoutput (SISO) nonlinear system. It can be linearized about an equibrium point <amsmath>x = x_e</amsmath>, <amsmath>u = u_e</amsmath>, <amsmath>y = y_e</amsmath> by defining new variables  
<center>  <center>  
−  +  <amsmath>z = x  x_e \qquad v = u  u_e \qquad w = y  h(x_e, u_e)</amsmath>.  
−  +  
−  +  
−  +  
−  +  
−  +  
−  +  
−  +  
−  +  
−  +  
−  +  
−  +  
−  +  
−  +  
−  +  
−  +  
</center>  </center>  
The dynamics of the system near the equilibrium point can then be approximated by the linear system  The dynamics of the system near the equilibrium point can then be approximated by the linear system  
−  <center>  +  <center><amsmath> 
−  {  +  \begin{aligned} 
−  +  \dot x &= A x + B u \\  
−  +  y &= C x + D u  
−  +  \end{aligned}  
−  +  </amsmath></center>  
−  +  
−  </center>  +  
where  where  
−  <center>  +  <center><amsmath> 
−  {  +  \begin{aligned} 
−  +  A &=\left.\frac{\partial f(x, u)}{\partial x}\right_{x_e, u_e} &\quad  
−  +  B &= \left.\frac{\partial f(x, u)}{\partial u}\right_{x_e, u_e} \\  
−  +  C &= \left.\frac{\partial h(x, u)}{\partial x}\right_{x_e, u_e} &\quad  
−  +  D &= \left.\frac{\partial y(x, u)}{\partial u}\right_{x_e, u_e}  
−  +  \end{aligned}  
−  +  </amsmath></center>  
−  +  
−  +  
−  +  
−  </center>  +  
The equilibrium point for a nonlinear system is locally asymptotically stable if the real part of the eigenvalues of the linearization about that equilibrium point have strictly negative real part.  The equilibrium point for a nonlinear system is locally asymptotically stable if the real part of the eigenvalues of the linearization about that equilibrium point have strictly negative real part.  
−  </p>  +  </p></li> 
</ol>  </ol>  
−  == Exercises ==  +  {{chaptertable begin}} 
+  {{chaptertable left}}  
+  
+  <!  
+  The following exercises cover some of the topics introduced in this chapter. Exercises marked with a * appear in the printed text.  
+  == Additional Exercises ==  
<ncl>Linear Systems Exercises</ncl>  <ncl>Linear Systems Exercises</ncl>  
+  >  
+  == Frequently Asked Questions ==  
+  <ncl>Linear Systems FAQ</ncl>  
+  == Errata ==  
+  <ncl>Linear Systems errata v2.11b</ncl>  
+  * [[:Category:Linear Systems errataFull list of errata starting from first printing]]  
+  * {{submitbug}}  
+  Minor corrections:  
+  * After equation (5.31), the parameter <amsmath>c</amsmath> listed in the text should be <amsmath>b_g</amsmath>  
+  {{chaptertable right}}  
+  
+  == MATLAB code ==  
+  The following MATLAB scripts are available for producing figures that appear in this chapter.  
+  * Figure 5.1: {{matlabfilelinsyslinearity.m}}, {{matlabfile.springmass.m}}  
+  * Figure 5.2: {{matlabfilelinsyspiecewise.m}}, {{matlabfile.springmass.m}}  
+  * Figure 5.5: {{matlabfilelinsysmodes.m}}  
+  * Figure 5.6: {{matlabfilelinsyspulses.m}}  
+  * Figure 5.8: {{matlabfilelinsystransient.m}}, {{matlabfile.springmass.m}}  
+  * Figure 5.9, 5.11: {{matlabfilelinsysresponse.m}}  
+  * Figure 5.10: {{matlabfilelinsyscompartment_stepresp.m}}  
+  * Figure 5.12: {{matlabfilelinsysopamp_bandbode.m}}  
+  * Figure 5.13: {{matlabfilelinsysafm_freqresp.m}}, {{matlabfileafmafm_data.m}}  
+  * Figure 5.14: {{matlabfilelinsyscruisepi_lin_nl_ch5.m}}, {{matlabfile.cruise_conpar.m}}, {{matlabfile.cruise_lin.m}}  
+  See the [[softwaresoftware page]] for more information on how to run these scripts.  
−  
−  
== Additional Information ==  == Additional Information ==  
−  
* [[Other Time Domain Specifications]]  * [[Other Time Domain Specifications]]  
+  {{chaptertable end}} 
Latest revision as of 22:45, 4 November 2012
Prev: Dynamic Behavior  Chapter 5  Linear Systems  Next: State Feedback 
Previous chapters have focused on the dynamics of a system with relatively little attention to the inputs and outputs. This chapter gives an introduction to input/output behavior for linear systems and shows how a nonlinear system can be approximated near an equilibrium point by a linear model.
Textbook Contents

Lecture MaterialsSupplemental Information

Chapter Summary
This chapter introduces the analysis tools for linear input/output systems.

A linear system is one in which the output is jointly linear in the intitial condition for the system and the input to the system. In particular, a linear system has the property that if we apply an input with zero initial condition, the corresponding output will be , where is the output associated with the input . This propery is called linear superposition.
A differential equation of the form
is a singleinput, singleoutput (SISO) linear differential equation. Its solution can be written in terms of the matrix exponential
The solution to the differential equation is given by the convolution equation
A linear system
is asymptotically stable if and only if all eigenvalues of all have strictly negative real part and is unstable if any eigenvalue of has strictly positive real part. For systems with eigenvalues having zero realpart, stability is determined by using the Jordan normal form associated with the matrix. A system with eigenvalues that have no strictly positive real part is stable if and only if the Jordan block corresponding to each eigenvalue with zero part is a scalar (1x1) block. The input/output response of a (stable) linear system contains a transient region portion, which eventually decays to zero, and a steady state portion, which persists over time. Two special responses are the step response, which is the output corresponding to an step input applied at and the frequency response, which is the response of the system to a sinusoidal input at a given frequency.
The step response is characterized by the following parameters:
 The steady state value, , of a step response is the final level of the output, assuming it converges.
 The rise time, , is the amount of time required for the signal to go from 10% of its final value to 90% of its final value.
 The overshoot, , is the percentage of the initial value by which the signal initially rises above the final value.
 The settling time, , is the amount of time required for the signal to stay within 5% of its final value for all future times.
The frequency response is given by
where and . The gain and phase of the frequency response are given by
A nonlinear system of the form
is a singleinput, singleoutput (SISO) nonlinear system. It can be linearized about an equibrium point , , by defining new variables
The dynamics of the system near the equilibrium point can then be approximated by the linear system
where
The equilibrium point for a nonlinear system is locally asymptotically stable if the real part of the eigenvalues of the linearization about that equilibrium point have strictly negative real part.
Frequently Asked Questions
Errata
Minor corrections: 
MATLAB codeThe following MATLAB scripts are available for producing figures that appear in this chapter.
See the software page for more information on how to run these scripts. Additional Information 