Let ( α, β be the eigenvector associated with the eigenvalue λ i
|
|
- Arleen Barber
- 6 years ago
- Views:
Transcription
1 ENGI 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 eigenvalues (wih λ < λ < 0 ): Two linearly independen soluions are x, y e λ, e λ and x, y e λ λ α α, e ( ) ( ) ( ) ( ) The general soluion is ( x, y ) ( c e λ ) c e λ, c e λ λ α + α + c e ( ) ( ) One can see ha x y lim, 0, 0. All orbis herefore erminae a he criical poin a he origin. The sysem is asympoically sable. If boh arbirary consans are zero, hen we have he rivial soluion (x y 0 for all ). If one of he arbirary consans is zero (say c ), hen x, y cα e λ, ce λ y x α which is a sraigh line hrough he origin, of slope α. [The siuaion is similar if c is zero.] ( ) ( ) We herefore obain sraigh-line rajecories ending a he singular poin, when exacly one of he arbirary consans is zero.
2 ENGI Sabiliy Analysis (Linear) Page 4.6 If neiher arbirary consan is zero, hen λ λ ( λλ) ( λλ) y c e + ce c + ce ce + c x λ λ ( λλ) ( λλ) cα e + c α e cα + c α e cα e + c α Because λ < λ < 0, lim and lim y x y x ( λ λ ) + c e c lim ( λ λ ) α c e c α + α ( λ λ ) + ( λλ) α α + α c c e lim c c e All orbis herefore come in from infiniy parallel o he line y x. α All orbis share he same angen a he origin, y x. α We obain a sable node ha is also asympoically sable. [The case illusraed here is α, α 3,,, λ 5, λ 0, which is + 3 generaed from A 4.] Case of real, disinc, posiive eigenvalues (wih λ > λ > 0 ): The analysis leads o he same phase space, excep ha he arrows are reversed. The resul is an unsable node.
3 ENGI Sabiliy Analysis (Linear) Page 4.7 Case of real, disinc eigenvalues of opposie sign (wih λ < 0 < λ): The general soluion is ( x, y ) ( c e λ ) c e λ, c e λ λ α + α + c e x y x y ( ) ( ( ) ( )) λ < 0 < λ lim, and lim, do no exis (infinie), (wih he excepion of he orbi for c 0). All orbis (excep c 0) herefore move away from he criical poin a he origin. The sysem is unsable. If boh arbirary consans are zero, hen we have he rivial soluion (x y 0 for all ). If one of he arbirary consans is zero (say c ), hen λ λ x, y cα e, ce y x α which is a sraigh line hrough he origin, of slope α. [The siuaion is similar if c is zero.] ( ) ( ) We herefore obain sraigh-line rajecories when one of he arbirary consans is zero. One of hem (c 0) ends a he singular poin while he oher begins here. If neiher arbirary consan is zero, hen λ λ ( λλ) ( λλ) y c e + ce c + ce ce + c x λ λ ( λλ) ( λλ) cα e + c α e cα + cα e cα e + c α Because λ < 0 < λ, lim and lim y x y x ( λ λ ) + c e c lim ( λ λ ) α c e c α + α ( λ λ ) + ( λλ) α α + α c c e lim c c e
4 ENGI Sabiliy Analysis (Linear) Page 4.8 All orbis herefore share he same asympoes, y x (incoming) and α y x (ougoing). α We obain a saddle poin, which is an unsable criical poin. [The case illusraed here is α, α 3,,, λ + 5, λ 5, which is generaed from A ]
5 ENGI Sabiliy Analysis (Linear) Page 4.9 Case of real, equal, negaive eigenvalues ( λ λ < 0) and b c 0: The sysem is uncoupled: dx d ax dy d dy and equal eigenvalues now require a d λ. x, y λ λ c e, c e. ( ) ( ) The general soluion is λ 0 lim ( x, y ) (, ) and lim ( x, y ) ( 0 ) <,0. All orbis herefore erminae a he criical poin a he origin. The sysem is asympoically sable. If boh arbirary consans are zero, hen we have he rivial soluion (x y 0 for all ). y c c 0 x c and c 0, c 0 x 0 The orbis are sraigh lines ending a he criical poin a he origin. The criical poin is an asympoically sable sar-shaped node. Addiional Noe: The eigenvalues of any riangular marix are he diagonal enries of ha marix: a b The characerisic equaion of A 0 d de A λ I 0 is ( ) a λ b 0 d λ ( a λ)( d λ) 0 λ a o r d
6 ENGI Sabiliy Analysis (Linear) Page 4.0 Case of real, equal, negaive eigenvalues ( λ λ < 0) and b, c no boh zero: The characerisic equaion λ ( a d ) λ ( ad bc) ( ) ( ) ( ) + + has he discriminan a + d 4 ad bc a d + 4bc 0. a+ d The soluion of he characerisic equaion simplifies o λ. x, y c c e λ λ α + α, c + c e. ( ) The general soluion is ( ) ( ) ( ) λ 0 lim x, y, and lim x, y 0 ( ) ( ) 0 ( ) ( ) <,0. All orbis herefore erminae a he criical poin a he origin. The sysem is asympoically sable. If boh arbirary consans are zero, hen we have he rivial soluion (x y 0 for all ). y c + c If c 0, hen as ± x cα + cα α All orbis (excep for c 0) herefore come in from infiniy parallel o he line y x, which is also a angen a he origin. I can be shown ha α when c 0, so ha he rajecories for c 0 and c 0 are boh y x. α α α Neiher eigenvalue can be zero, oherwise (0, 0) is no he only criical poin (as shown on page 4.4).
7 ENGI Sabiliy Analysis (Linear) Page 4. Case of real, equal, posiive eigenvalues ( λ λ > 0 ) The analysis leads o he same phase planes as in he case of real equal negaive eigenvalues, bu he signs of he arrows are reversed and he resul is an unsable node. Case of complex conjugae pair of eigenvalues wih negaive real par The eigenvalues (roos of he characerisic equaion) are λ a + jb, λ a jb, ( a<0 ). The general soluion has he form x c( Acosb Asin b) c( Asin b Acosb) + + e y c( Bcosb Bsin b) c( Bsin b Bcosb) + + e Using he definiions ( ) + ( + ), B ( c B c B ) + ( c B + c B ) A ca ca ca ca ca + ca ca ca cb + cb cb cb cos α, sinα, cos, sin A A B B he general soluion can be wrien more compacly as a a x, y A e cos b + α, B e cos b + a a ( ) ( ) ( ) ( ) ( ) ( ) ( ) ( ) a< 0 lim x, y, and lim x, y 0,0. π If x 0 hen b n n π If y 0 hen b+ + nπ n y Bcos( b+ ) x Acos( b+ α ) y is periodic, wih period π. x b The orbis spiral in o he origin. + α + π ( ) ( ) We have an asympoically sable spiral, also known as a sable focus.
8 ENGI Sabiliy Analysis (Linear) Page 4. Case of complex conjugae pair of eigenvalues wih posiive real par The analysis leads o he same phase planes as in he case of negaive real par, bu he signs of he arrows are reversed and he resul is an unsable focus. Case of complex conjugae pair of eigenvalues wih zero real par (pure imaginary) The eigenvalues (roos of he characerisic equaion) are λ jb, λ + jb. The general soluion has he compac form ( x, y ) ( Acos ( b+ α ), Bcos( b+ )) π If α 0 and, hen x y ( x, y ) ( Acos b, Bsinb) A + B so ha he orbis are ellipses, cenred on he criical poin a he origin. This is a sable cenre. Oher choices of α and also lead o concenric ses of ellipses, bu roaed wih respec o he coordinaes axes. Noe ha his is he only case of a sable criical poin ha is no asympoically sable.
9 ENGI Sabiliy Analysis (Linear) Page 4.3 Summary for he Linear Sysem dx dy ax + by, d d cx + dy, abcd,,, consans Characerisic equaion: λ a + d λ + ad bc 0 Discriminan ( ) ( ) ( ) ( ) ( ) ( ) D a+ d 4 ad bc a d + 4bc a b Roos of characerisic equaion ( eigenvalues of A ): c d Cases: λ ( ) a+ d ± D a + d D oher condiion λ Type of poin a + d < 0 D > 0 ad bc > 0 real, disinc Sable negaive node a + d < 0 D 0 b c 0 real, equal Sable negaive sar shape a + d < 0 D 0 b, c no boh 0 real, equal Sable negaive node a + d < 0 D < 0 complex Sable conjugae pair focus [spiral] a + d 0 D < 0 Pure Sable imaginary pair cenre a + d > 0 D > 0 ad bc > 0 real, disinc Unsable posiive node (any) D > 0 ad bc < 0 real, disinc Unsable opposie signs saddle poin a + d > 0 D 0 b c 0 real, equal Unsable posiive sar shape a + d > 0 D 0 b, c no boh 0 real, equal Unsable posiive node a + d > 0 D < 0 complex Unsable conjugae pair focus [spiral] Noe ha ad bc de A and ha a + d he race of he marix A. In brief, if he real pars of boh eigenvalues are negaive (or boh zero), hen he origin is sable. Oherwise i is unsable. [See also he example a "
10 ENGI Sabiliy Analysis (Linear) Page 4.4 Example Find he naure of he criical poin of he sysem dx dy 4x 3 y, 5x4y d d and find he general soluion. a b 4 3 The coefficien marix is A c d 5 4. race(a) a + d 4 + ( 4) 0 D (a d) + 4bc (4 + 4) + 4( 3)(5) > de A < D > 0 and ad bc < 0 λ are real wih opposie signs and he criical poin is a saddle poin (unsable). Solving he sysem: ( a+ d) ± D 0 ± 4 λ ± x, y cα e + c α e, c e + c e ( ) ( ) α where is he eigenvecor associaed wih he eigenvalue λ α and is he eigenvecor associaed wih he eigenvalue λ +. To find he eigenvecors, find non-zero soluions o he equaion a λ b α 0 c d λ 0 A λ : 4+ 3 α 5 3 α Any non-zero choice such ha 5α 3 0 will provide an eigenvecor. α 3 Selec. 5
11 ENGI Sabiliy Analysis (Linear) Page 4.5 Example (coninued) A λ +: 4 3 α 3 3 α Any non-zero choice such ha α 0 will provide an eigenvecor. α Selec. The general soluion is ( x, y ) ( 3 ce, 5ce ) ]. [I is simple o check ha (4x 3y, 5x 4y) is indeed equal o ( x, y) Also noe ha y 5ce y ( ) 5 y( ) lim ( c 0) and lim c 0 x 3ce x 3 + x so ha all orbis for which boh c and c are non-zero share he same asympoes, 3y 5x (which is he incoming orbi, when c 0) and y x (which is he ougoing orbi, when c 0). ( ) A few represenaive orbis and he wo asympoes are ploed in his phase space diagram:
12 ENGI Sabiliy Analysis (Linear) Page 4.6 Example Find he naure of he criical poin of he sysem dx dy x + y, d d x y and find he general soluion. a b The coefficien marix is A c d. race(a) a + d + 4 < 0. D (a d) + 4bc ( + ) > 0 λ are real, disinc and negaive and de A ad bc 4 3 > 0 he criical poin is a sable node. Solving he sysem: ( a+ d) ± D 4 ± 4 λ ± 3,, 3, 3 x y cα e + c α e c e + c e ( ) ( ) α where is he eigenvecor associaed wih he eigenvalue λ 3 α and is he eigenvecor associaed wih he eigenvalue λ. To find he eigenvecors, find non-zero soluions o he equaion a λ b α 0 c d λ 0 A λ 3: + 3 α α Any non-zero choice such ha α + 0 will provide an eigenvecor. α Selec. A λ : + α α Any non-zero choice such ha α + 0 will provide an eigenvecor. α Selec.
13 ENGI Sabiliy Analysis (Linear) Page 4.7 Example (coninued) The general soluion is ( ) ( 3, 3 x y ce, ce ) ]. [I is simple o check ha ( x + y, x y) is indeed equal o ( x, y) Also noe ha y 3 ce ce + c c x 3 ce ce + c c y y lim ( c 0) and lim c 0 x + x and ( ) ( ) ( 3 3 x y ce ce ) ( ) lim, lim, 0, 0 so ha all orbis for which boh c and c are non-zero come in from a direcion parallel o y x (which is he orbi when c 0) and share he same angen a he origin, y x (which is he orbi when c 0). A few represenaive orbis and he common angen are ploed in his phase space diagram:
After 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 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 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 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 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 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 informationSolutions for homework 12
y Soluions for homework Secion Nonlinear sysems: The linearizaion of a nonlinear sysem Consider he sysem y y y y y (i) Skech he nullclines Use a disincive marking for each nullcline so hey can be disinguished
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 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 informationPredator - Prey Model Trajectories and the nonlinear conservation law
Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories
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 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 informationChapter #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 informationBifurcation Analysis of a Stage-Structured Prey-Predator System with Discrete and Continuous Delays
Applied Mahemaics 4 59-64 hp://dx.doi.org/.46/am..4744 Published Online July (hp://www.scirp.org/ournal/am) Bifurcaion Analysis of a Sage-Srucured Prey-Predaor Sysem wih Discree and Coninuous Delays Shunyi
More informationExercises: Similarity Transformation
Exercises: Similariy Transformaion Problem. Diagonalize he following marix: A [ 2 4 Soluion. Marix A has wo eigenvalues λ 3 and λ 2 2. Since (i) A is a 2 2 marix and (ii) i has 2 disinc eigenvalues, we
More informationStability and Bifurcation in a Neural Network Model with Two Delays
Inernaional Mahemaical Forum, Vol. 6, 11, no. 35, 175-1731 Sabiliy and Bifurcaion in a Neural Nework Model wih Two Delays GuangPing Hu and XiaoLing Li School of Mahemaics and Physics, Nanjing Universiy
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 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 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 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 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 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 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 motions of the celt on a horizontal plane with viscous friction
The h Join Inernaional Conference on Mulibody Sysem Dynamics June 8, 18, Lisboa, Porugal The moions of he cel on a horizonal plane wih viscous fricion Maria A. Munisyna 1 1 Moscow Insiue of Physics and
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 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 informationApplication 5.4 Defective Eigenvalues and Generalized Eigenvectors
Applicaion 5.4 Defecive Eigenvalues and Generalized Eigenvecors The goal of his applicaion is he soluion of he linear sysems like where he coefficien marix is he exoic 5-by-5 marix x = Ax, (1) 9 11 21
More informationChapter #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 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 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 informationFor example, the comb filter generated from. ( ) has a transfer function. e ) has L notches at ω = (2k+1)π/L and L peaks at ω = 2π k/l,
Comb Filers The simple filers discussed so far are characeried eiher by a single passband and/or a single sopband There are applicaions where filers wih muliple passbands and sopbands are required The
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 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 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 informationMath 1. Two-Hours Exam December 10, 2017.
Mah. Two-Hours Exam December, 7. JE/JKL..7 Problem resar;wih(linearalgebra): Given he inhomogeneous linear sysem of equaions lign:=x-*x+3*x3=a^+*a-3; lign x x 3 x3 = a a 3 lign:=x+*x-*x3=a^+3; lign x x
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 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 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 informationHamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t
M ah 5 2 7 Fall 2 0 0 9 L ecure 1 0 O c. 7, 2 0 0 9 Hamilon- J acobi Equaion: Explici Formulas In his lecure we ry o apply he mehod of characerisics o he Hamilon-Jacobi equaion: u + H D u, x = 0 in R n
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 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 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 information4. Advanced Stability Theory
Applied Nonlinear Conrol Nguyen an ien - 4 4 Advanced Sabiliy heory he objecive of his chaper is o presen sabiliy analysis for non-auonomous sysems 41 Conceps of Sabiliy for Non-Auonomous Sysems Equilibrium
More informationCHAPTER 2 Signals And Spectra
CHAPER Signals And Specra Properies of Signals and Noise In communicaion sysems he received waveform is usually caegorized ino he desired par conaining he informaion, and he undesired par. he desired par
More informationThe equation to any straight line can be expressed in the form:
Sring Graphs Par 1 Answers 1 TI-Nspire Invesigaion Suden min Aims Deermine a series of equaions of sraigh lines o form a paern similar o ha formed by he cables on he Jerusalem Chords Bridge. Deermine he
More informationV L. DT s D T s t. Figure 1: Buck-boost converter: inductor current i(t) in the continuous conduction mode.
ECE 445 Analysis and Design of Power Elecronic Circuis Problem Se 7 Soluions Problem PS7.1 Erickson, Problem 5.1 Soluion (a) Firs, recall he operaion of he buck-boos converer in he coninuous conducion
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 information3, so θ = arccos
Mahemaics 210 Professor Alan H Sein Monday, Ocober 1, 2007 SOLUTIONS This problem se is worh 50 poins 1 Find he angle beween he vecors (2, 7, 3) and (5, 2, 4) Soluion: Le θ be he angle (2, 7, 3) (5, 2,
More informationWEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x
WEEK-3 Reciaion PHYS 131 Ch. 3: FOC 1, 3, 4, 6, 14. Problems 9, 37, 41 & 71 and Ch. 4: FOC 1, 3, 5, 8. Problems 3, 5 & 16. Feb 8, 018 Ch. 3: FOC 1, 3, 4, 6, 14. 1. (a) The horizonal componen of he projecile
More information4.5 Constant Acceleration
4.5 Consan Acceleraion v() v() = v 0 + a a() a a() = a v 0 Area = a (a) (b) Figure 4.8 Consan acceleraion: (a) velociy, (b) acceleraion When he x -componen of he velociy is a linear funcion (Figure 4.8(a)),
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 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 informationChapter Three Systems of Linear Differential Equations
Chaper Three Sysems of Linear Differenial Equaions In his chaper we are going o consier sysems of firs orer orinary ifferenial equaions. These are sysems of he form x a x a x a n x n x a x a x a n x n
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 informationOutline of Topics. Analysis of ODE models with MATLAB. What will we learn from this lecture. Aim of analysis: Why such analysis matters?
of Topics wih MATLAB Shan He School for Compuaional Science Universi of Birmingham Module 6-3836: Compuaional Modelling wih MATLAB Wha will we learn from his lecure Aim of analsis: Aim of analsis. Some
More informationRobotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1.
Roboics I April 11, 017 Exercise 1 he kinemaics of a 3R spaial robo is specified by he Denavi-Harenberg parameers in ab 1 i α i d i a i θ i 1 π/ L 1 0 1 0 0 L 3 0 0 L 3 3 able 1: able of DH parameers of
More information1 st order ODE Initial Condition
Mah-33 Chapers 1-1 s Order ODE Sepember 1, 17 1 1 s order ODE Iniial Condiion f, sandard form LINEAR NON-LINEAR,, p g differenial form M x dx N x d differenial form is equivalen o a pair of differenial
More informationTHE 2-BODY PROBLEM. FIGURE 1. A pair of ellipses sharing a common focus. (c,b) c+a ROBERT J. VANDERBEI
THE 2-BODY PROBLEM ROBERT J. VANDERBEI ABSTRACT. In his shor noe, we show ha a pair of ellipses wih a common focus is a soluion o he 2-body problem. INTRODUCTION. Solving he 2-body problem from scrach
More informationarxiv: v1 [math.gm] 4 Nov 2018
Unpredicable Soluions of Linear Differenial Equaions Mara Akhme 1,, Mehme Onur Fen 2, Madina Tleubergenova 3,4, Akylbek Zhamanshin 3,4 1 Deparmen of Mahemaics, Middle Eas Technical Universiy, 06800, Ankara,
More informationSection 2.2 Charge and Current 2.6 b) The current direction is designated as the direction of the movement of positive charges.
Chaper Soluions Secion. Inroducion. Curren source. Volage source. esisor.4 Capacior.5 Inducor Secion. Charge and Curren.6 b) The curren direcion is designaed as he direcion of he movemen of posiive charges..7
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 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 informationMath 315: Linear Algebra Solutions to Assignment 6
Mah 35: Linear Algebra s o Assignmen 6 # Which of he following ses of vecors are bases for R 2? {2,, 3, }, {4,, 7, 8}, {,,, 3}, {3, 9, 4, 2}. Explain your answer. To generae he whole R 2, wo linearly independen
More informationLinear Dynamic Models
Linear Dnamic Models and Forecasing Reference aricle: Ineracions beween he muliplier analsis and he principle of acceleraion Ouline. The sae space ssem as an approach o working wih ssems of difference
More informationTests of Nonlinear Resonse Theory. We compare the results of direct NEMD simulation against Kawasaki and TTCF for 2- particle colour conductivity.
ess of Nonlinear Resonse heory We compare he resuls of direc NEMD simulaion agains Kawasaki and CF for 2- paricle colour conduciviy. 0.0100 0.0080 0.0060 Direc CF, BK and RK J x 0.0040 0.0020 DIR KAW CF
More informationDIFFERENTIAL GEOMETRY HW 5
DIFFERENTIAL GEOMETRY HW 5 CLAY SHONKWILER 3. Le M be a complee Riemannian manifold wih non-posiive secional curvaure. Prove ha d exp p v w w, for all p M, all v T p M and all w T v T p M. Proof. Le γ
More informationPhysics 221 Fall 2008 Homework #2 Solutions Ch. 2 Due Tues, Sept 9, 2008
Physics 221 Fall 28 Homework #2 Soluions Ch. 2 Due Tues, Sep 9, 28 2.1 A paricle moving along he x-axis moves direcly from posiion x =. m a ime =. s o posiion x = 1. m by ime = 1. s, and hen moves direcly
More informationOptimality Conditions for Unconstrained Problems
62 CHAPTER 6 Opimaliy Condiions for Unconsrained Problems 1 Unconsrained Opimizaion 11 Exisence Consider he problem of minimizing he funcion f : R n R where f is coninuous on all of R n : P min f(x) x
More information2. Nonlinear Conservation Law Equations
. Nonlinear Conservaion Law Equaions One of he clear lessons learned over recen years in sudying nonlinear parial differenial equaions is ha i is generally no wise o ry o aack a general class of nonlinear
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 informationDynamical Systems, Phase Plan Analysis- Linearization Method
Journal of Mulidisciplinary Engineering Science Technology (JMEST) ISSN: 458-943 Vol. 3 Issue 8, Augus - 16 Dynamical Sysems, Phase Plan Analysis- Linearizaion Mehod Alfred DACI Deparmen of Mahemaics,
More informationTraveling Waves. Chapter Introduction
Chaper 4 Traveling Waves 4.1 Inroducion To dae, we have considered oscillaions, i.e., periodic, ofen harmonic, variaions of a physical characerisic of a sysem. The sysem a one ime is indisinguishable from
More informationThe Asymptotic Behavior of Nonoscillatory Solutions of Some Nonlinear Dynamic Equations on Time Scales
Advances in Dynamical Sysems and Applicaions. ISSN 0973-5321 Volume 1 Number 1 (2006, pp. 103 112 c Research India Publicaions hp://www.ripublicaion.com/adsa.hm The Asympoic Behavior of Nonoscillaory Soluions
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 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 informationEE363 homework 1 solutions
EE363 Prof. S. Boyd EE363 homework 1 soluions 1. LQR for a riple accumulaor. We consider he sysem x +1 = Ax + Bu, y = Cx, wih 1 1 A = 1 1, B =, C = [ 1 ]. 1 1 This sysem has ransfer funcion H(z) = (z 1)
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 informationCHAPTER 12 DIRECT CURRENT CIRCUITS
CHAPTER 12 DIRECT CURRENT CIUITS DIRECT CURRENT CIUITS 257 12.1 RESISTORS IN SERIES AND IN PARALLEL When wo resisors are conneced ogeher as shown in Figure 12.1 we said ha hey are conneced in series. As
More information1 Solutions to selected problems
1 Soluions o seleced problems 1. Le A B R n. Show ha in A in B bu in general bd A bd B. Soluion. Le x in A. Then here is ɛ > 0 such ha B ɛ (x) A B. This shows x in B. If A = [0, 1] and B = [0, 2], hen
More informationYimin Math Centre. 4 Unit Math Homework for Year 12 (Worked Answers) 4.1 Further Geometric Properties of the Ellipse and Hyperbola...
4 Uni Mah Homework for Year 12 (Worked Answers) Suden Name: Grade: Dae: Score: Table of conens 4 Topic 2 Conics (Par 4) 1 4.1 Furher Geomeric Properies of he Ellipse and Hyperbola............... 1 4.2
More informationLinear Quadratic Regulator (LQR) - State Feedback Design
Linear Quadrai Regulaor (LQR) - Sae Feedbak Design A sysem is expressed in sae variable form as x = Ax + Bu n m wih x( ) R, u( ) R and he iniial ondiion x() = x A he sabilizaion problem using sae variable
More informationAsymptotic instability of nonlinear differential equations
Elecronic Journal of Differenial Equaions, Vol. 1997(1997), No. 16, pp. 1 7. ISSN: 172-6691. URL: hp://ejde.mah.sw.edu or hp://ejde.mah.un.edu fp (login: fp) 147.26.13.11 or 129.12.3.113 Asympoic insabiliy
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 information15. Vector Valued Functions
1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,
More information4.6 One Dimensional Kinematics and Integration
4.6 One Dimensional Kinemaics and Inegraion When he acceleraion a( of an objec is a non-consan funcion of ime, we would like o deermine he ime dependence of he posiion funcion x( and he x -componen of
More informationSignal and System (Chapter 3. Continuous-Time Systems)
Signal and Sysem (Chaper 3. Coninuous-Time Sysems) Prof. Kwang-Chun Ho kwangho@hansung.ac.kr Tel: 0-760-453 Fax:0-760-4435 1 Dep. Elecronics and Informaion Eng. 1 Nodes, Branches, Loops A nework wih b
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 informationMEI Mechanics 1 General motion. Section 1: Using calculus
Soluions o Exercise MEI Mechanics General moion Secion : Using calculus. s 4 v a 6 4 4 When =, v 4 a 6 4 6. (i) When = 0, s = -, so he iniial displacemen = - m. s v 4 When = 0, v = so he iniial velociy
More informationA Necessary and Sufficient Condition for the Solutions of a Functional Differential Equation to Be Oscillatory or Tend to Zero
JOURNAL OF MAEMAICAL ANALYSIS AND APPLICAIONS 24, 7887 1997 ARICLE NO. AY965143 A Necessary and Sufficien Condiion for he Soluions of a Funcional Differenial Equaion o Be Oscillaory or end o Zero Piambar
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 informationModule 2 F c i k c s la l w a s o s f dif di fusi s o i n
Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms
More information18 Biological models with discrete time
8 Biological models wih discree ime The mos imporan applicaions, however, may be pedagogical. The elegan body of mahemaical heory peraining o linear sysems (Fourier analysis, orhogonal funcions, and so
More informationTHE SINE INTEGRAL. x dt t
THE SINE INTEGRAL As one learns in elemenary calculus, he limi of sin(/ as vanishes is uniy. Furhermore he funcion is even and has an infinie number of zeros locaed a ±n for n1,,3 Is plo looks like his-
More information!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)
"#"$%&#'()"#&'(*%)+,&',-)./)1-*) #$%&'()*+,&',-.%,/)*+,-&1*#$)()5*6$+$%*,7&*-'-&1*(,-&*6&,7.$%$+*&%'(*8$&',-,%'-&1*(,-&*6&,79*(&,%: ;..,*&1$&$.$%&'()*1$$.,'&',-9*(&,%)?%*,('&5
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 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 informationnon -negative cone Population dynamics motivates the study of linear models whose coefficient matrices are non-negative or positive.
LECTURE 3 Linear/Nonnegaive Marix Models x ( = Px ( A= m m marix, x= m vecor Linear sysems of difference equaions arise in several difference conexs: Linear approximaions (linearizaion Perurbaion analysis
More informationModule 4: Time Response of discrete time systems Lecture Note 2
Module 4: Time Response of discree ime sysems Lecure Noe 2 1 Prooype second order sysem The sudy of a second order sysem is imporan because many higher order sysem can be approimaed by a second order model
More informationHamilton Jacobi equations
Hamilon Jacobi equaions Inoducion o PDE The rigorous suff from Evans, mosly. We discuss firs u + H( u = 0, (1 where H(p is convex, and superlinear a infiniy, H(p lim p p = + This by comes by inegraion
More informationPhysics 20 Lesson 5 Graphical Analysis Acceleration
Physics 2 Lesson 5 Graphical Analysis Acceleraion I. Insananeous Velociy From our previous work wih consan speed and consan velociy, we know ha he slope of a posiion-ime graph is equal o he velociy of
More information