Chapter #1 EEE8013 EEE3001. Linear Controller Design and State Space Analysis
|
|
- Karen Blankenship
- 5 years ago
- Views:
Transcription
1 Chaper EEE83 EEE3 Chaper # EEE83 EEE3 Linear Conroller Design and Sae Space Analysis Ordinary Differenial Equaions.... Inroducion.... Firs Order ODEs Second Order ODEs General Maerial Roos are real and unequal Roos are Complex (and hence no equal) Roos are real and equal Tuorial Exercise I... 9 Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk /9
2 Chaper EEE83 EEE3 Ordinary Differenial Equaions. Inroducion To undersand he propeies (dynamics) of a sysem, we can model (represen) i using differenial equaions (DEs). The response/behaviour of he sysem is found by solving he DEs. In our cases, he DE is an Ordinary DE (ODE), i.e. no a paial derivaive. The main purpose of his Chaper is o learn how o solve firs and second order ODEs in he ime domain. This will serve as a building block o model and sudy more complicaed sysems. Our ulimae goal is o conrol he sysem when i does no show a saisfacory behaviour. Effecively, his will be done by modifying he ODE. Noe for EEE83 sudens: There are foonoes hroughou he noes, which is assessed maerial!. Firs Order ODEs The general form of a firs order ODE is: d f x, () where x, Analyical soluion: Explici formula for x() (a soluion which can be found using various mehods) which saisfies f ( x, ) d The proper noaion is x() and no x bu we drop he brackes in order o simplify he presenaion. Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk /9
3 Chaper EEE83 EEE3 Example.: Prove ha x e and x e are soluions of 3x d. de 3x 3e 3e 3e d d de 3x 3e 3e 3e d d Obviously here are infinie soluions o an ODE and for ha reason he found soluion is called he General Soluion of he ODE. Firs order Iniial Value Problem : f ( x, ), x x An iniial value problem is an ODE wih an iniial condiion, hence we do no find he general soluion bu he Specific Soluion ha passes hrough x a =. d Analyical soluion: Explici formula for x() which saisfies f ( x, ) and d passes hrough x when. Example.: Prove ha x e is a soluion, while x e is no a soluion of 3, xx d Boh expressions ( x e and x e ) saisfy he 3x bu a = d x e x x e x 3 clc, clear all, syms, x=exp(-3*); =diff(x,); isequal(,-3*x) x=-*exp(-3*); =diff(x,); isequal(,-3*x) 3 clc, clear all, syms, x=exp(-3*); x=-*exp(-3*); x=; x_=double(subs(x,,)); x_=double(subs(x,,)); isequal(x,x_), isequal(x,x_), Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 3/9
4 Chaper EEE83 EEE3 For ha reason some books use a differen symbol for he specific soluion:,, x. You mus be clear abou he difference beween an ODE and he soluion o an IVP! From now on we will jus sudy IVP unless oherwise explicily menioned. Linear Firs Order ODEs A linear s order ODE is given by: x' b x c, a a ax' bx c, a Non auonomous Auonomous () wih a,b,c and a. In engineering books he mos common form of () is (since a ): x' kx u (3) wih k,u Noe: We say ha u is he inpu o our sysem ha is represened by (3) The soluion of (3) (using he inegraing facor) is given by: k k k x e x e e u d Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 4/9
5 Chaper EEE83 EEE3 k The erm e x is called ransien response, while comes k k e e u d from he inpu signal u. If we assume ha u is consan: k k k x e x e e ud x e k x u e k k Hence: lim x u u/ k, k k, k Thus we say ha if k> he sysem is sable (and he soluion converges exponenially a u/k) while if k< he sysem is unsable (and he soluion diverges exponenially o, ). Example.3: u= and k= & 5, x= x e,limx,as >.5 Transien Toal.5 Transien Toal Inpu componen Inpu componen clc, clear all, syms x(), =diff(x); dsolve(+*x, x()==) Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 5/9
6 Chaper EEE83 EEE3 Example.4: u= and k=- & 5, x= Toal Transien Inpu componen.5 x Toal Transien Inpu componen.5.5 Example.5: u= and k=5, x= & Toal Transien Inpu componen 4 3 Toal Transien Inpu componen Example.6: u=- & and k=5, x=.5 Transien Inpu componen Toal Toal Transien Inpu componen clc, clear all, close all, syms, x=; k=5; =; u=; =:.:; x_x=exp(-k*)*x; x_u=exp(-k*)*in(exp(k*)*u,,); x_x_=double(subs(x_x,,)); x_u_=double(subs(x_u,,)); hold on, plo(,x_x_), plo(,x_u_), plo(,x_u_+x_x_) Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 6/9
7 Chaper EEE83 EEE3 Commens: In real sysems we canno have a sae (say he speed of a mass-spring sysem) ha becomes infinie, obviously he sysem will be desroyed when x ges o a high value. For he dynamics (seling ime, sabiliy ) of he sysem we should only focus on he homogenous ODE: x' k x 3. Second Order ODEs 3. General Maerial A second order ODE has as a general form: f x', x, d (4) A linear nd order ODE is given by: x'' A x' B x u, Non auonomous x'' Ax' Bx u, Auonomous (5) And again we focus on auonomous homogeneous sysems: x'' A x' B x (6) Again we define as an analyical soluion of (6) an expression ha saisfies i. Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 7/9
8 Chaper EEE83 EEE3 Example.7: Given x'' x' 3x prove ha soluions: e '' e ' 3 e 9e 6e 3e e '' e ' 3 e x e and x e are wo e e 3e 6 Assume ha you have soluions for a nd order ODE x and x (we will see laer how o ge hese wo soluions), hen: x A x B x x A x B x obviously I can muliply hese wo equaions wih arbirary consans: Cx CA x CB x Cx CA x CB x and now I can add hem and collec similar erms: Cx C x '' A Cx C x ' B Cx C x Common Term Common Term Common Term which means ha Cx Cx is also a soluion of he ODE. (i.e. he linear combinaion of x and x) 6 clc, clear all, close all, syms x() Dx=diff(x); Dx=diff(x,); ODE=Dx-*Dx-3*x; subs(ode, x, exp(-)), subs(ode, x, exp(3*)) Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 8/9
9 Chaper EEE83 EEE3 Example.8: Given x'' x' 3x prove ha e e '' e e ' 3 e e e e e e e e e e 6e 4e 3e 6e 9e 6e 3e e 4e 6e x e e is a soluion: 7 Now, he quesion is, if we have x and x, can ALL oher soluions of he ODE, be expressed as a linear combinaion of x and x? So assume a hird soluion : A B '' ' Now, he quesion can be wrien as, can we find consans C and C such as: Cx Cx ' C x ' C x ' This equaion can be seen as a by sysem wih unknowns C and C as: C x x x' x' C ' From linear algebra his sysem of equaions has a unique soluion if: 7 clc, clear all, close all, syms x() ; Dx=diff(x); Dx=diff(x,); ODE=Dx-*Dx-3*x; subs(ode, x, exp(3*)+exp(-)) Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 9/9
10 Chaper EEE83 EEE3 ' ' x x x x x x ' x x ' Noe: The marix W x, x of he ODE. ' ' x x x x is called he Wronskian 8 We also know from linear algebra ha he deerminan is no zero if: ' ' x x C x x So if he wo soluions x and x are linear independen (LI) hen ANY oher soluion can be described by he linear combinaion of x and x. So now we have o look for wo LI soluions for he nd order ODE. Example.9: Prove ha wo soluions of x'' x' 3x, x e are linear independen. x e and x x e e e e Wx, x W x' x' 3e e 3e e W e e e e e e e 9 8 From he Polish mahemaician Józef Maria Hoëne-Wroński 9 clc, clear all, close all, syms, x=exp(3*); x=exp(-); Dx=diff(x); Dx=diff(x); W=[x, x; Dx, Dx], de(w) Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk /9
11 Chaper EEE83 EEE3 Example.: Prove ha wo soluions of x'' x' 3x, x e are NOT linear independen. 3 x e and W x x x x e e, x' x' 3e 6e e e W e e 3e 6e Example.: For he ODE x'' x' 3x prove ha he soluion 3 x e e canno be wrien as any combinaion of x e and x e x Cx Cx e e Ce Ce C C e From his expression we have ha C C (and hence we have he erm e ) bu here is no erm e for e. Bu how can we find wo LI soluions? For homogeneous s order ODEs wih k u= he soluion was: x e C so we will ry a similar approach for nd order ODEs: x '' Ax' Bx, assume x e => x' re & x' ' r e => x '' Ax' Bx r e Are Be r Ar B (7) This is called he Characerisic Equaion (CE) and we have o check is roos: r A A 4B, hese are he Characerisic values or Eigenvalues. Noice ha we do NOT know wha is he value of r. Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk /9
12 Chaper EEE83 EEE3 3. Roos are real and unequal If A 4B he sysem is called Overdamped and he wo roos are r and r wih r r, r, r. Then soluions as: ' ' r x x e e x x re re e and r x e x e re e re are wo linear independen hence he general soluion is r Cx Cx Ce Ce (8) x If r and r < hen x and he sysem is sable. If r or r > hen x and he sysem is unsable. Example.: The CE of x'' x' 3x is r r3 which means 43 r 5 ha he wo roos are: r, r 6 5 x e e and hence he LI soluions are 6 x e e 5 6 This means ha he general soluion is x Ce Ce and hence he ODE is sable. The Wronskian is x x e e 6 5 e x x e e e e e e 5 6 ' ' 5 6 If he iniial condiion is x, x' hen: roos([ 3]) clc, clear all, close all, syms x() Dx=diff(x,); Dx=diff(x,); ODE=Dx+*Dx+3*x; dsolve(ode) Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk /9
13 Chaper EEE83 EEE3 C C C 6 x 6e 5e 5C6C C Roos are Complex (and hence no equal) If A 4B hen he sysem is called Underdamped and he wo roos are r abj and r r abj wih r r, r, r. Then a bj and x e e are wo linear independen soluions as x e e a bj e abj e a bj e a bj e e a bj e abj abj abj e abj e a a a a abje abje e abjabje bj a bj a bj a bj a bj Hence he general soluion is x Cx Cx Ce Ce (9) bu remember ha C and C are complex now variables such as x. Example.3: The CE of x'' x' 5x is r r5 which means 6 4j r j ha he wo roos are: r, j r j j x e e and hence he LI soluions are j x e e 3 clc, clear all, close all, syms x() Dx=diff(x,); Dx=diff(x,); ODE=Dx+*Dx+3*x; dsolve(ode, x()==, Dx()==) Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 3/9
14 Chaper EEE83 EEE3 j This means ha he general soluion is x Ce C e ODE is sable. The Wronskian is x x e j j j x ' x' j e j e j j j j j e e j e e j e j e j j e 4je If he iniial condiion is x jc e j and hence he j, x' hen: C C j C 4 jc C j x je je 4 4 j j 4 An alernaive approach is no o use x & x bu a linear combinaion of hem: y e e, y e e e e e e Noe ha re re re re Using Euler s formula: abj a e e cosb jsin b and hence: a bj a bj a a y e e e cosb jsinbcosb jsinb e cosb a bj a bj a a y e e e cosb jsinb cosb jsinb je sinb Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 4/9
15 Chaper EEE83 EEE3 As y and y are soluions so do y, y. So he general soluion when j we have complex roos is: a x e CcosbC sin b, CC, () Example.4: The CE of x'' x' 5x is r r5 which means 6 4j r j ha he wo roos are: r, j r j x e cos and hence he LI soluions are x e sin This means ha he general soluion is x e Ccos Csin and hence he ODE is sable. The Wronskian is x x e cos e sin e x ' x ' e cos e sin e sin e cos If he iniial condiion is x, x' hen: C C CC C.5 x e cos.5sin 3.3 Roos are real and equal If A 4Bhen he sysem is called Criically damped and he wo roos are r r r wih r. One soluion is x e bu how abou x? We can use x e and he general soluion: x C x C x Ce C e () Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 5/9
16 Chaper EEE83 EEE3 The Wronskian is: e e re e e e e e re e e e e e Example.5: The CE of x'' x' x is r r which means ha r he wo roos are: r, r x e and hence he LI soluions are x e This means ha he general soluion is x Ce Ce and hence he ODE is sable. The Wronskian is e e e e e e If he iniial condiion is x C C CC C x e e, x' hen: No assessed maerial To see why x e is he nd soluion go o he ODE and place e '' Ae ' Bxe r Ar B Since r is a double roo of he CE: a. So: e '' A e ' Bxe a r r Taking he ime derivaive w r: x e : r ArBa r r for some consan Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 6/9
17 Chaper EEE83 EEE3 '' ' d e d e d e d e a r r A B dr dr dr dr And as we can change he sequence of he differeniaion: d e ar r '' ' d e d e d e A B dr dr dr dr By using simple calculus: '' ' e Ae Beear r ear r '' ' d e d a r r e Ae Be ar r e dr dr '' ' By placing now where r=r: Which means ha e e re re e And hence e Ae Be e mus be a soluion of my ODE and: e re e e re re e re e x e is my second soluion. Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 7/9
18 Chaper EEE83 EEE3 Roo Space jb jb a a Sable Unsable Name Oscillaions? Componens of soluion Overdamped No Two exponenials: Criically damped No k k e, e, k k, Two exponenials: k k e, e, k Underdamped Yes One exponenial and one cosine e k cos,, k Undamped Yes one cosine cos Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 8/9
19 Chaper EEE83 EEE3 4. Tuorial Exercise I. By using he general form of he analyic soluion ry o predic he response of he following sysems. Your answer mus describe he sysem as sable/unsable, convergen o zero/nonzero value. Crosscheck your answer by solving he DE: d 5 6x, x, x, x d 5 6x, x, x, x d 5 6x, x, x, x d 5 6x, x, x, x d 3, x, x, x. Find he soluion of x 6x 5x, x, x 3. Briefly describe how he soluion behaves for hese iniial condiions. Draw a skech of he response. 3. Find he soluion of x x 6x, x, x. Briefly describe how he soluion behaves for hese iniial condiions. Draw a skech of he response. 4. Find he soluion of x x.5x, x, x / 3. Briefly describe how he soluion behaves for hese iniial condiions. Draw a skech of he response. 5. Find he Wronskian marices of he soluions of Q-5. Module Leader: Dr Damian Giaouris - damian.giaouris@ncl.ac.uk 9/9
Chapter #1 EEE8013 EEE3001. Linear Controller Design and State Space Analysis
Chaper EEE83 EEE3 Chaper # EEE83 EEE3 Linear Conroller Design and Sae Space Analysis Ordinary Differenial Equaions.... Inroducion.... Firs Order ODEs... 3. Second Order ODEs... 7 3. General Maerial...
More informationSome Basic Information about M-S-D Systems
Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,
More informationODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004
ODEs II, Lecure : Homogeneous Linear Sysems - I Mike Raugh March 8, 4 Inroducion. In he firs lecure we discussed a sysem of linear ODEs for modeling he excreion of lead from he human body, saw how o ransform
More informationChapter 2. First Order Scalar Equations
Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.
More informationSystem of Linear Differential Equations
Sysem of Linear Differenial Equaions In "Ordinary Differenial Equaions" we've learned how o solve a differenial equaion for a variable, such as: y'k5$e K2$x =0 solve DE yx = K 5 2 ek2 x C_C1 2$y''C7$y
More informationChapter 6. Systems of First Order Linear Differential Equations
Chaper 6 Sysems of Firs Order Linear Differenial Equaions We will only discuss firs order sysems However higher order sysems may be made ino firs order sysems by a rick shown below We will have a sligh
More informationt is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...
Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger
More informationMath Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.
Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need
More informationAfter the completion of this section the student. Theory of Linear Systems of ODEs. Autonomous Systems. Review Questions and Exercises
Chaper V ODE V.5 Sysems of Ordinary Differenial Equaions 45 V.5 SYSTEMS OF FIRST ORDER LINEAR ODEs Objecives: Afer he compleion of his secion he suden - should recall he definiion of a sysem of linear
More information10. State Space Methods
. Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he
More informationt + t sin t t cos t sin t. t cos t sin t dt t 2 = exp 2 log t log(t cos t sin t) = Multiplying by this factor and then integrating, we conclude that
ODEs, Homework #4 Soluions. Check ha y ( = is a soluion of he second-order ODE ( cos sin y + y sin y sin = 0 and hen use his fac o find all soluions of he ODE. When y =, we have y = and also y = 0, so
More informationSolutions of Sample Problems for Third In-Class Exam Math 246, Spring 2011, Professor David Levermore
Soluions of Sample Problems for Third In-Class Exam Mah 6, Spring, Professor David Levermore Compue he Laplace ransform of f e from is definiion Soluion The definiion of he Laplace ransform gives L[f]s
More informationLAPLACE TRANSFORM AND TRANSFER FUNCTION
CHBE320 LECTURE V LAPLACE TRANSFORM AND TRANSFER FUNCTION Professor Dae Ryook Yang Spring 2018 Dep. of Chemical and Biological Engineering 5-1 Road Map of he Lecure V Laplace Transform and Transfer funcions
More informationLet us start with a two dimensional case. We consider a vector ( x,
Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our
More informationSection 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients
Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous
More informationWall. x(t) f(t) x(t = 0) = x 0, t=0. which describes the motion of the mass in absence of any external forcing.
MECHANICS APPLICATIONS OF SECOND-ORDER ODES 7 Mechanics applicaions of second-order ODEs Second-order linear ODEs wih consan coefficiens arise in many physical applicaions. One physical sysems whose behaviour
More informationTwo Coupled Oscillators / Normal Modes
Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own
More informationChapter 3 Boundary Value Problem
Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le
More informationSolutions to Assignment 1
MA 2326 Differenial Equaions Insrucor: Peronela Radu Friday, February 8, 203 Soluions o Assignmen. Find he general soluions of he following ODEs: (a) 2 x = an x Soluion: I is a separable equaion as we
More informationAnnouncements: Warm-up Exercise:
Fri Apr 13 7.1 Sysems of differenial equaions - o model muli-componen sysems via comparmenal analysis hp//en.wikipedia.org/wiki/muli-comparmen_model Announcemens Warm-up Exercise Here's a relaively simple
More informationdt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.
Mah Rahman Exam Review Soluions () Consider he IVP: ( 4)y 3y + 4y = ; y(3) = 0, y (3) =. (a) Please deermine he longes inerval for which he IVP is guaraneed o have a unique soluion. Soluion: The disconinuiies
More informationMath 334 Fall 2011 Homework 11 Solutions
Dec. 2, 2 Mah 334 Fall 2 Homework Soluions Basic Problem. Transform he following iniial value problem ino an iniial value problem for a sysem: u + p()u + q() u g(), u() u, u () v. () Soluion. Le v u. Then
More informationMTH Feburary 2012 Final term PAPER SOLVED TODAY s Paper
MTH401 7 Feburary 01 Final erm PAPER SOLVED TODAY s Paper Toal Quesion: 5 Mcqz: 40 Subjecive quesion: 1 4 q of 5 marks 4 q of 3 marks 4 q of marks Guidelines: Prepare his file as I included all pas papers
More informationDifferential Equations
Mah 21 (Fall 29) Differenial Equaions Soluion #3 1. Find he paricular soluion of he following differenial equaion by variaion of parameer (a) y + y = csc (b) 2 y + y y = ln, > Soluion: (a) The corresponding
More informationMA 366 Review - Test # 1
MA 366 Review - Tes # 1 Fall 5 () Resuls from Calculus: differeniaion formulas, implici differeniaion, Chain Rule; inegraion formulas, inegraion b pars, parial fracions, oher inegraion echniques. (1) Order
More informationLet ( α, β be the eigenvector associated with the eigenvalue λ i
ENGI 940 4.05 - Sabiliy Analysis (Linear) Page 4.5 Le ( α, be he eigenvecor associaed wih he eigenvalue λ i of he coefficien i i) marix A Le c, c be arbirary consans. a b c d Case of real, disinc, negaive
More informationModule 3: The Damped Oscillator-II Lecture 3: The Damped Oscillator-II
Module 3: The Damped Oscillaor-II Lecure 3: The Damped Oscillaor-II 3. Over-damped Oscillaions. This refers o he siuaion where β > ω (3.) The wo roos are and α = β + α 2 = β β 2 ω 2 = (3.2) β 2 ω 2 = 2
More informationElectrical Circuits. 1. Circuit Laws. Tools Used in Lab 13 Series Circuits Damped Vibrations: Energy Van der Pol Circuit
V() R L C 513 Elecrical Circuis Tools Used in Lab 13 Series Circuis Damped Vibraions: Energy Van der Pol Circui A series circui wih an inducor, resisor, and capacior can be represened by Lq + Rq + 1, a
More informationMA 214 Calculus IV (Spring 2016) Section 2. Homework Assignment 1 Solutions
MA 14 Calculus IV (Spring 016) Secion Homework Assignmen 1 Soluions 1 Boyce and DiPrima, p 40, Problem 10 (c) Soluion: In sandard form he given firs-order linear ODE is: An inegraing facor is given by
More informationSimulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010
Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid
More information1 1 + x 2 dx. tan 1 (2) = ] ] x 3. Solution: Recall that the given integral is improper because. x 3. 1 x 3. dx = lim dx.
. Use Simpson s rule wih n 4 o esimae an () +. Soluion: Since we are using 4 seps, 4 Thus we have [ ( ) f() + 4f + f() + 4f 3 [ + 4 4 6 5 + + 4 4 3 + ] 5 [ + 6 6 5 + + 6 3 + ]. 5. Our funcion is f() +.
More informationu(x) = e x 2 y + 2 ) Integrate and solve for x (1 + x)y + y = cos x Answer: Divide both sides by 1 + x and solve for y. y = x y + cos x
. 1 Mah 211 Homework #3 February 2, 2001 2.4.3. y + (2/x)y = (cos x)/x 2 Answer: Compare y + (2/x) y = (cos x)/x 2 wih y = a(x)x + f(x)and noe ha a(x) = 2/x. Consequenly, an inegraing facor is found wih
More informationMath Final Exam Solutions
Mah 246 - Final Exam Soluions Friday, July h, 204 () Find explici soluions and give he inerval of definiion o he following iniial value problems (a) ( + 2 )y + 2y = e, y(0) = 0 Soluion: In normal form,
More informationMath 23 Spring Differential Equations. Final Exam Due Date: Tuesday, June 6, 5pm
Mah Spring 6 Differenial Equaions Final Exam Due Dae: Tuesday, June 6, 5pm Your name (please prin): Insrucions: This is an open book, open noes exam. You are free o use a calculaor or compuer o check your
More informationMath 333 Problem Set #2 Solution 14 February 2003
Mah 333 Problem Se #2 Soluion 14 February 2003 A1. Solve he iniial value problem dy dx = x2 + e 3x ; 2y 4 y(0) = 1. Soluion: This is separable; we wrie 2y 4 dy = x 2 + e x dx and inegrae o ge The iniial
More informationConcourse Math Spring 2012 Worked Examples: Matrix Methods for Solving Systems of 1st Order Linear Differential Equations
Concourse Mah 80 Spring 0 Worked Examples: Marix Mehods for Solving Sysems of s Order Linear Differenial Equaions The Main Idea: Given a sysem of s order linear differenial equaions d x d Ax wih iniial
More informationMATH 128A, SUMMER 2009, FINAL EXAM SOLUTION
MATH 28A, SUMME 2009, FINAL EXAM SOLUTION BENJAMIN JOHNSON () (8 poins) [Lagrange Inerpolaion] (a) (4 poins) Le f be a funcion defined a some real numbers x 0,..., x n. Give a defining equaion for he Lagrange
More informationOscillation of an Euler Cauchy Dynamic Equation S. Huff, G. Olumolode, N. Pennington, and A. Peterson
PROCEEDINGS OF THE FOURTH INTERNATIONAL CONFERENCE ON DYNAMICAL SYSTEMS AND DIFFERENTIAL EQUATIONS May 4 7, 00, Wilmingon, NC, USA pp 0 Oscillaion of an Euler Cauchy Dynamic Equaion S Huff, G Olumolode,
More informationdy dx = xey (a) y(0) = 2 (b) y(1) = 2.5 SOLUTION: See next page
Assignmen 1 MATH 2270 SOLUTION Please wrie ou complee soluions for each of he following 6 problems (one more will sill be added). You may, of course, consul wih your classmaes, he exbook or oher resources,
More informationME 391 Mechanical Engineering Analysis
Fall 04 ME 39 Mechanical Engineering Analsis Eam # Soluions Direcions: Open noes (including course web posings). No books, compuers, or phones. An calculaor is fair game. Problem Deermine he posiion of
More information23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes
Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals
More informationSecond-Order Differential Equations
WWW Problems and Soluions 3.1 Chaper 3 Second-Order Differenial Equaions Secion 3.1 Springs: Linear and Nonlinear Models www m Problem 3. (NonlinearSprings). A bod of mass m is aached o a wall b means
More information( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is
UNIT IMPULSE RESPONSE, UNIT STEP RESPONSE, STABILITY. Uni impulse funcion (Dirac dela funcion, dela funcion) rigorously defined is no sricly a funcion, bu disribuion (or measure), precise reamen requires
More informationVoltage/current relationship Stored Energy. RL / RC circuits Steady State / Transient response Natural / Step response
Review Capaciors/Inducors Volage/curren relaionship Sored Energy s Order Circuis RL / RC circuis Seady Sae / Transien response Naural / Sep response EE4 Summer 5: Lecure 5 Insrucor: Ocavian Florescu Lecure
More informationIntermediate Differential Equations Review and Basic Ideas
Inermediae Differenial Equaions Review and Basic Ideas John A. Burns Cener for Opimal Design And Conrol Inerdisciplinary Cener forappliedmahemaics Virginia Polyechnic Insiue and Sae Universiy Blacksburg,
More informationMEI STRUCTURED MATHEMATICS 4758
OXFORD CAMBRIDGE AND RSA EXAMINATIONS Advanced Subsidiary General Cerificae of Educaion Advanced General Cerificae of Educaion MEI STRUCTURED MATHEMATICS 4758 Differenial Equaions Thursday 5 JUNE 006 Afernoon
More informationSection 7.4 Modeling Changing Amplitude and Midline
488 Chaper 7 Secion 7.4 Modeling Changing Ampliude and Midline While sinusoidal funcions can model a variey of behaviors, i is ofen necessary o combine sinusoidal funcions wih linear and exponenial curves
More informationINDEX. Transient analysis 1 Initial Conditions 1
INDEX Secion Page Transien analysis 1 Iniial Condiions 1 Please inform me of your opinion of he relaive emphasis of he review maerial by simply making commens on his page and sending i o me a: Frank Mera
More information5.1 - Logarithms and Their Properties
Chaper 5 Logarihmic Funcions 5.1 - Logarihms and Their Properies Suppose ha a populaion grows according o he formula P 10, where P is he colony size a ime, in hours. When will he populaion be 2500? We
More informationMorning Time: 1 hour 30 minutes Additional materials (enclosed):
ADVANCED GCE 78/0 MATHEMATICS (MEI) Differenial Equaions THURSDAY JANUARY 008 Morning Time: hour 30 minues Addiional maerials (enclosed): None Addiional maerials (required): Answer Bookle (8 pages) Graph
More information2.7. Some common engineering functions. Introduction. Prerequisites. Learning Outcomes
Some common engineering funcions 2.7 Inroducion This secion provides a caalogue of some common funcions ofen used in Science and Engineering. These include polynomials, raional funcions, he modulus funcion
More information(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)
Phase Plane Analysis of Linear Sysems Adaped from Applied Nonlinear Conrol by Sloine and Li The general form of a linear second-order sysem is a c b d From and b bc d a Differeniaion of and hen subsiuion
More informationEXERCISES FOR SECTION 1.5
1.5 Exisence and Uniqueness of Soluions 43 20. 1 v c 21. 1 v c 1 2 4 6 8 10 1 2 2 4 6 8 10 Graph of approximae soluion obained using Euler s mehod wih = 0.1. Graph of approximae soluion obained using Euler
More information2.9 Modeling: Electric Circuits
SE. 2.9 Modeling: Elecric ircuis 93 2.9 Modeling: Elecric ircuis Designing good models is a ask he compuer canno do. Hence seing up models has become an imporan ask in modern applied mahemaics. The bes
More information1 birth rate γ (number of births per time interval) 2 death rate δ proportional to size of population
Scienific Comuing I Module : Poulaion Modelling Coninuous Models Michael Bader Par I ODE Models Lehrsuhl Informaik V Winer 7/ Discree vs. Coniuous Models d d = F,,...) ) =? discree model: coninuous model:
More informationY 0.4Y 0.45Y Y to a proper ARMA specification.
HG Jan 04 ECON 50 Exercises II - 0 Feb 04 (wih answers Exercise. Read secion 8 in lecure noes 3 (LN3 on he common facor problem in ARMA-processes. Consider he following process Y 0.4Y 0.45Y 0.5 ( where
More informationMath Week 15: Section 7.4, mass-spring systems. These are notes for Monday. There will also be course review notes for Tuesday, posted later.
Mah 50-004 Week 5: Secion 7.4, mass-spring sysems. These are noes for Monday. There will also be course review noes for Tuesday, posed laer. Mon Apr 3 7.4 mass-spring sysems. Announcemens: Warm up exercise:
More informationPhysics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle
Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,
More informationTHE WAVE EQUATION. part hand-in for week 9 b. Any dilation v(x, t) = u(λx, λt) of u(x, t) is also a solution (where λ is constant).
THE WAVE EQUATION 43. (S) Le u(x, ) be a soluion of he wave equaion u u xx = 0. Show ha Q43(a) (c) is a. Any ranslaion v(x, ) = u(x + x 0, + 0 ) of u(x, ) is also a soluion (where x 0, 0 are consans).
More informationMon Apr 9 EP 7.6 Convolutions and Laplace transforms. Announcements: Warm-up Exercise:
Mah 225-4 Week 3 April 9-3 EP 7.6 - convoluions; 6.-6.2 - eigenvalues, eigenvecors and diagonalizabiliy; 7. - sysems of differenial equaions. Mon Apr 9 EP 7.6 Convoluions and Laplace ransforms. Announcemens:
More informationComplete solutions to Exercise 14(b) 1. Very similar to EXAMPLE 4. We have same characteristic equation:
Soluions 4(b) Complee soluions o Exercise 4(b). Very similar o EXAMPE 4. We have same characerisic equaion: 5 i Ae = + Be By using he given iniial condiions we obain he simulaneous equaions A+ B= 0 5A
More informationMath 10B: Mock Mid II. April 13, 2016
Name: Soluions Mah 10B: Mock Mid II April 13, 016 1. ( poins) Sae, wih jusificaion, wheher he following saemens are rue or false. (a) If a 3 3 marix A saisfies A 3 A = 0, hen i canno be inverible. True.
More informationMA Study Guide #1
MA 66 Su Guide #1 (1) Special Tpes of Firs Order Equaions I. Firs Order Linear Equaion (FOL): + p() = g() Soluion : = 1 µ() [ ] µ()g() + C, where µ() = e p() II. Separable Equaion (SEP): dx = h(x) g()
More informationChapter Q1. We need to understand Classical wave first. 3/28/2004 H133 Spring
Chaper Q1 Inroducion o Quanum Mechanics End of 19 h Cenury only a few loose ends o wrap up. Led o Relaiviy which you learned abou las quarer Led o Quanum Mechanics (1920 s-30 s and beyond) Behavior of
More informationMath Wednesday March 3, , 4.3: First order systems of Differential Equations Why you should expect existence and uniqueness for the IVP
Mah 2280 Wednesda March 3, 200 4., 4.3: Firs order ssems of Differenial Equaions Wh ou should epec eisence and uniqueness for he IVP Eample: Consider he iniial value problem relaed o page 4 of his eserda
More informationd 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3
and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or
More informationMath 36. Rumbos Spring Solutions to Assignment #6. 1. Suppose the growth of a population is governed by the differential equation.
Mah 36. Rumbos Spring 1 1 Soluions o Assignmen #6 1. Suppose he growh of a populaion is governed by he differenial equaion where k is a posiive consan. d d = k (a Explain why his model predics ha he populaion
More informationOrdinary Differential Equations
Ordinary Differenial Equaions 5. Examples of linear differenial equaions and heir applicaions We consider some examples of sysems of linear differenial equaions wih consan coefficiens y = a y +... + a
More informationTopic Astable Circuits. Recall that an astable circuit has two unstable states;
Topic 2.2. Asable Circuis. Learning Objecives: A he end o his opic you will be able o; Recall ha an asable circui has wo unsable saes; Explain he operaion o a circui based on a Schmi inverer, and esimae
More informationCh.1. Group Work Units. Continuum Mechanics Course (MMC) - ETSECCPB - UPC
Ch.. Group Work Unis Coninuum Mechanics Course (MMC) - ETSECCPB - UPC Uni 2 Jusify wheher he following saemens are rue or false: a) Two sreamlines, corresponding o a same insan of ime, can never inersec
More informationOrdinary dierential equations
Chaper 5 Ordinary dierenial equaions Conens 5.1 Iniial value problem........................... 31 5. Forward Euler's mehod......................... 3 5.3 Runge-Kua mehods.......................... 36
More informationMath 2214 Sol Test 2B Spring 2015
Mah 14 Sol Tes B Sring 015 roblem 1: An objec weighing ounds sreches a verical sring 8 fee beond i naural lengh before coming o res a equilibrium The objec is ushed u 6 fee from i s equilibrium osiion
More informationADVANCED MATHEMATICS FOR ECONOMICS /2013 Sheet 3: Di erential equations
ADVANCED MATHEMATICS FOR ECONOMICS - /3 Shee 3: Di erenial equaions Check ha x() =± p ln(c( + )), where C is a posiive consan, is soluion of he ODE x () = Solve he following di erenial equaions: (a) x
More informationBoyce/DiPrima/Meade 11 th ed, Ch 3.1: 2 nd Order Linear Homogeneous Equations-Constant Coefficients
Boce/DiPrima/Meade h ed, Ch 3.: nd Order Linear Homogeneous Equaions-Consan Coefficiens Elemenar Differenial Equaions and Boundar Value Problems, h ediion, b William E. Boce, Richard C. DiPrima, and Doug
More informationMath 2142 Exam 1 Review Problems. x 2 + f (0) 3! for the 3rd Taylor polynomial at x = 0. To calculate the various quantities:
Mah 4 Eam Review Problems Problem. Calculae he 3rd Taylor polynomial for arcsin a =. Soluion. Le f() = arcsin. For his problem, we use he formula f() + f () + f ()! + f () 3! for he 3rd Taylor polynomial
More information2.3 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS
Andrei Tokmakoff, MIT Deparmen of Chemisry, 2/22/2007 2-17 2.3 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS The mahemaical formulaion of he dynamics of a quanum sysem is no unique. So far we have described
More informationENGI 9420 Engineering Analysis Assignment 2 Solutions
ENGI 940 Engineering Analysis Assignmen Soluions 0 Fall [Second order ODEs, Laplace ransforms; Secions.0-.09]. Use Laplace ransforms o solve he iniial value problem [0] dy y, y( 0) 4 d + [This was Quesion
More informationReview - Quiz # 1. 1 g(y) dy = f(x) dx. y x. = u, so that y = xu and dy. dx (Sometimes you may want to use the substitution x y
Review - Quiz # 1 (1) Solving Special Tpes of Firs Order Equaions I. Separable Equaions (SE). d = f() g() Mehod of Soluion : 1 g() d = f() (The soluions ma be given implicil b he above formula. Remember,
More informationChapter 8 The Complete Response of RL and RC Circuits
Chaper 8 The Complee Response of RL and RC Circuis Seoul Naional Universiy Deparmen of Elecrical and Compuer Engineering Wha is Firs Order Circuis? Circuis ha conain only one inducor or only one capacior
More informationMATH 2050 Assignment 9 Winter Do not need to hand in. 1. Find the determinant by reducing to triangular form for the following matrices.
MATH 2050 Assignmen 9 Winer 206 Do no need o hand in Noe ha he final exam also covers maerial afer HW8, including, for insance, calculaing deerminan by row operaions, eigenvalues and eigenvecors, similariy
More informationDynamic Effects of Feedback Control!
Dynamic Effecs of Feedback Conrol! Rober Sengel! Roboics and Inelligen Sysems MAE 345, Princeon Universiy, 2017 Inner, Middle, and Ouer Feedback Conrol Loops Sep Response of Linear, Time- Invarian (LTI)
More information6.2 Transforms of Derivatives and Integrals.
SEC. 6.2 Transforms of Derivaives and Inegrals. ODEs 2 3 33 39 23. Change of scale. If l( f ()) F(s) and c is any 33 45 APPLICATION OF s-shifting posiive consan, show ha l( f (c)) F(s>c)>c (Hin: In Probs.
More informationSection 3.8, Mechanical and Electrical Vibrations
Secion 3.8, Mechanical and Elecrical Vibraions Mechanical Unis in he U.S. Cusomary and Meric Sysems Disance Mass Time Force g (Earh) Uni U.S. Cusomary MKS Sysem CGS Sysem fee f slugs seconds sec pounds
More informationh[n] is the impulse response of the discrete-time system:
Definiion Examples Properies Memory Inveribiliy Causaliy Sabiliy Time Invariance Lineariy Sysems Fundamenals Overview Definiion of a Sysem x() h() y() x[n] h[n] Sysem: a process in which inpu signals are
More informationSignals and Systems Profs. Byron Yu and Pulkit Grover Fall Midterm 1 Solutions
8-90 Signals and Sysems Profs. Byron Yu and Pulki Grover Fall 07 Miderm Soluions Name: Andrew ID: Problem Score Max 0 8 4 6 5 0 6 0 7 8 9 0 6 Toal 00 Miderm Soluions. (0 poins) Deermine wheher he following
More informationt 2 B F x,t n dsdt t u x,t dxdt
Evoluion Equaions For 0, fixed, le U U0, where U denoes a bounded open se in R n.suppose ha U is filled wih a maerial in which a conaminan is being ranspored by various means including diffusion and convecion.
More informationMath 334 Test 1 KEY Spring 2010 Section: 001. Instructor: Scott Glasgow Dates: May 10 and 11.
1 Mah 334 Tes 1 KEY Spring 21 Secion: 1 Insrucor: Sco Glasgow Daes: Ma 1 and 11. Do NOT wrie on his problem saemen bookle, excep for our indicaion of following he honor code jus below. No credi will be
More informationSZG Macro 2011 Lecture 3: Dynamic Programming. SZG macro 2011 lecture 3 1
SZG Macro 2011 Lecure 3: Dynamic Programming SZG macro 2011 lecure 3 1 Background Our previous discussion of opimal consumpion over ime and of opimal capial accumulaion sugges sudying he general decision
More informationWeek 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem)
Week 1 Lecure Problems, 5 Wha if somehing oscillaes wih no obvious spring? Wha is ω? (problem se problem) Sar wih Try and ge o SHM form E. Full beer can in lake, oscillaing F = m & = ge rearrange: F =
More informationExam 1 Solutions. 1 Question 1. February 10, Part (A) 1.2 Part (B) To find equilibrium solutions, set P (t) = C = dp
Exam Soluions Februar 0, 05 Quesion. Par (A) To find equilibrium soluions, se P () = C = = 0. This implies: = P ( P ) P = P P P = P P = P ( + P ) = 0 The equilibrium soluion are hus P () = 0 and P () =..
More informationSolutionbank Edexcel AS and A Level Modular Mathematics
Page of 4 Soluionbank Edexcel AS and A Level Modular Mahemaics Exercise A, Quesion Quesion: Skech he graphs of (a) y = e x + (b) y = 4e x (c) y = e x 3 (d) y = 4 e x (e) y = 6 + 0e x (f) y = 00e x + 0
More informationLab 10: RC, RL, and RLC Circuits
Lab 10: RC, RL, and RLC Circuis In his experimen, we will invesigae he behavior of circuis conaining combinaions of resisors, capaciors, and inducors. We will sudy he way volages and currens change in
More informationRC, RL and RLC circuits
Name Dae Time o Complee h m Parner Course/ Secion / Grade RC, RL and RLC circuis Inroducion In his experimen we will invesigae he behavior of circuis conaining combinaions of resisors, capaciors, and inducors.
More informationBEng (Hons) Telecommunications. Examinations for / Semester 2
BEng (Hons) Telecommunicaions Cohor: BTEL/14/FT Examinaions for 2015-2016 / Semeser 2 MODULE: ELECTROMAGNETIC THEORY MODULE CODE: ASE2103 Duraion: 2 ½ Hours Insrucions o Candidaes: 1. Answer ALL 4 (FOUR)
More informationHOMEWORK # 2: MATH 211, SPRING Note: This is the last solution set where I will describe the MATLAB I used to make my pictures.
HOMEWORK # 2: MATH 2, SPRING 25 TJ HITCHMAN Noe: This is he las soluion se where I will describe he MATLAB I used o make my picures.. Exercises from he ex.. Chaper 2.. Problem 6. We are o show ha y() =
More informationChallenge Problems. DIS 203 and 210. March 6, (e 2) k. k(k + 2). k=1. f(x) = k(k + 2) = 1 x k
Challenge Problems DIS 03 and 0 March 6, 05 Choose one of he following problems, and work on i in your group. Your goal is o convince me ha your answer is correc. Even if your answer isn compleely correc,
More informatione a s a f t dt f t dt = p = p. t = a
Mah 225-4 Fri Apr 7 5, EP76 Today we finish discussing Laplace ransform echniques: Impulse forcing ("dela funcions")oday's noes Convoluion formulas o solve any inhomogeneous consan coefficien linear DE,
More information!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)
"#"$%&#'()"#&'(*%)+,&',-)./)1-*) #$%&'()*+,&',-.%,/)*+,-&1*#$)()5*6$+$%*,7&*-'-&1*(,-&*6&,7.$%$+*&%'(*8$&',-,%'-&1*(,-&*6&,79*(&,%: ;..,*&1$&$.$%&'()*1$$.,'&',-9*(&,%)?%*,('&5
More informationSecond Order Linear Differential Equations
Second Order Linear Differenial Equaions Second order linear equaions wih consan coefficiens; Fundamenal soluions; Wronskian; Exisence and Uniqueness of soluions; he characerisic equaion; soluions of homogeneous
More informationCHAPTER 6: FIRST-ORDER CIRCUITS
EEE5: CI CUI T THEOY CHAPTE 6: FIST-ODE CICUITS 6. Inroducion This chaper considers L and C circuis. Applying he Kirshoff s law o C and L circuis produces differenial equaions. The differenial equaions
More information