Ordinary Differential Equations
|
|
- Wendy Gaines
- 5 years ago
- Views:
Transcription
1 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 y n y = a y a y 1 1 n = + + yn an 1y1... anny n n n (5.1) and he linear n-h-order differenial equaion (n > 1) wih consan coefficiens x + a x a x + a x= ( n) ( n 1) n 1 1 (5.) We discuss suiable mehods o compue soluions and he fundamenal marix. In paricular, we compue he fundamenal marix for all planar sysems (n = ) and classify hese sysems according o heir properies in phase space. We also inroduce he noion of sabiliy for linear differenial equaions Planar sysems 1. Equaion of small oscillaions. We sudy he linear second-order differenial equaion x+ a x + a x= (5.3) 1 for real coefficiens a, a 1. This equaion arises from problems in mechanics and elecriciy. Oscillaions of a spring. A mass m coupled o a spring is submied o a linear force. Newon s equaion of moion reads mx = Dx where D is a posiive consan. I is of he form (5.3) wih a = D/m and a 1 =. If fricion is aken ino accoun hen a 1 is a posiive consan. Small oscillaions of he mahemaical pendulum. We consider a weighless rod of lengh l aached a one end and bearing a poin mass m a he oher end. If x denoes he angle by which he rod deviaes from he verical hen ml x = mgl sin x. For small deviaions sin(x) x and we ge an equaion of he form (5.3) wih a = g/l and a 1 =. 1
2 Again, if we ake ino accoun a force due o fricion ha is proporional o he velociy hen a = a 1 is a posiive consan. Invered pendulum. If we urn he pendulum upside down, hen he sign of he force due o graviaion changes sign and herefore a = - g/l. Elecrical circuis. The charge y of a condenser of capaciy C in a closed elecrical circui wih induciviy L and resisance R saisfies he second order differenial equaion 1 Ly + Ry + C y = which is of he form (5.3) wih a = 1/LC and a 1 = R/L Transformaion o a firs-order sysem. The second-order differenial equaion (5.3) is equivalen o a wo-componen linear firs-order sysem wih marix A given by 1 A = a a1 (5.4) I is convenien o rewrie A using he race and he deerminan or is eigenvalues. Then where 1 1 A = = de + ( A) r( A) λ1λ λ1 λ (5.5) 4de r A r A A λ1 = + 4de r A r A A λ = (5.6) are he (possibly complex) eigenvalues of A. If here is no fricion (no resisance), hen r(a) = which means ha he sysem is area preserving in phase space. 1.. Equaion in phase space. The equaion (5.3) describes a moion in he wodimensional phase space wih coordinaes x, y = x. The marix A defines he direcion field, which is a composiion of a roaion and a dilaion of he plane. The orbi of a x (), y (), or if we consider y as a funcion of x he soluion is he se of poins { } se of poins {x, y(x)} where y is (a leas formally) a soluion of he non-linear firs-order differenial equaion dy ( x ) = a x a1 = de A x + r A dx yx yx
3 which is of homogeneous ype and can be solved (a leas in principle) by elemenary inegraion mehods (see ch..3). Anoher imporan quaniy defined on he phase space is he energy of he sysem. The energy funcional associaed o (5.3) is defined by (, ) 1 1 E x x = x + a x (5.7) If x() is a soluion of (5.3) hen d E ( x, x ) = x ( x + a x ) = a x r 1 = ( A ) x (5.8) d If here is no fricion ( a 1 = -r(a) = ) he energy of a soluion remains consan and he soluion follows he lines of consan energy in phase space. In his case he sysem is called conservaive. Since he equaion is linear he energy defines a quadraic form and he orbis are eiher ellipses (de(a) = a > ), hyperbolas (de(a) = a < ), or sraigh lines (de(a) = a = ). In he presence of fricion in a mechanical sysem or elecrical resisance in an elecrical circui ( a 1 = -r(a) > ) he energy is decreasing Explici soluions of he second order equaion. According o he resuls of he previous chaper (heorem 4.7) he general complex valued soluion is of he form () x = cexp λ+ cexp λ, λ λ x = cexp λ + cexp λ, λ = λ λ 1 1 (5.9) The corresponding real soluions are lised in he following able. λ 1, λ real, r(a) 4de(A) λ 1 λ r(a) > 4de(A) λ 1 = λ = λ, r(a) = 4de(A) () = ( λ ) + ( λ ) = ( λ ) + ( x c1exp 1 cexp λ 1, λ complex conjugae, r(a) < 4de(A) x c exp c exp λ) 1 λ 1 = µ iω, λ = µ + iω λ 1 = iω, λ = iω, r(a) = () = exp( µ ) 1sin ( ω ) + cos( ω ) or 1 ( ω ) ( () = exp( µ ) sin ( ω + δ) x = Csin ( ω+ δ ) x c c x C x = c sin + c cos ω) or The invered pendulum belongs o he firs class since de(a) = a <. In all oher models de(a) = a > and damping -a 1 = r(a) For weak damping i.e. small r(a) here are wo complex conjugae, criical damping corresponds o a double real eigenvalue If boh eigenvalues are real and differen he sysem is called overdamped. 3
4 1.4. Compuing exp(a). A complex valued fundamenal marix of he associaed firs order sysem is easily compued if he eigenvalues are differen. The marix A defined in (5.4) resp. (5.5) has eigenvecors 1 1 v1 =, v = λ1 λ Therefore a fundamenal marix Y() is given by Since we ge ( λ1) exp( λ) exp exp Y() = ( v1exp( λ1) vexp( λ) ) = λ1exp λ1 λ λ 1 λv1 λ1v ( λ λ1) =, v v1= ( λ λ1) 1 (5.1) exp ( A) If λ 1 = µ iω, λ = µ + iω (5.11) reads λ1 λ λ λ1 1 λe λ1e e e = λ λ λλ( e e ) λ e λe λ λ λ λ (5.11) exp ( A) µ e ωcos ω µ sin ω sin ω = ω sin ( ω )( µ + ω ) ωcos( ω ) + µ sin ( ω ) (5.1) To compue exp(a) in case λ 1 = λ = λ we le λ 1 λ = λ. We obain exp 1 λ = λ 1 + λ (5.13) λ ( A) e The direc compuaion goes as follows: By heorem 4.7. he funcions exp(λ) and exp(λ) are soluions of he second order equaion. The corresponding soluion vecors are λ 1 λ 1 = e, v = e v λ 1+ λ from which we easily consruc a fundamenal marix. 4
5 . General firs order sysems and heir classificaion. Le A be an arbirary real marix. Then A is eiher similar o he marix given in (5.5) or A is a muliple of he uni marix E. In he firs case all resuls are similar o he resuls for he second order equaion. To classify he sysems i is sufficien o consider heir normal forms, which are one of he following real marices: λ1 λ 1 µ ω,, λ λ ω µ (5.14) where we used he noaions inroduced above. Diagonal case λ 1 < λ < λ 1 < < λ < λ 1 < λ Sable node Saddle poin Unsable node λ 1 = λ < λ 1 < = λ < λ 1 = λ Sable sar Marginally sable Unsable sar 5
6 Jordan case λ 1 = λ = λ < λ1 = λ = λ = λ 1 = λ = λ = Sable node Marginally unsable Unsable node Complex case µ <, ω > µ =, ω > µ >, ω > Sable focus Cener Unsable focus.1. Sabiliy. We precise he meaning of sabiliy and insabiliy of linear equaions. We noe ha y() = is always a saionary soluion of a linear homogeneous sysem. We define sabiliy and insabiliy as a propery of soluions wih iniial condiions close o he zero soluion. We say ha y() = is sable if all soluions remain bounded, asympoically sable if all soluions end o zero as ends o infiniy, and unsable if here is an unbounded soluion of he linear equaion. We have seen ha y() = is sable if here is no eigenvalue wih posiive real par ( wih excepion of he marginally unsable Jordan block), asympoically sable if all eigenvalues have negaive real par, and unsable if here is an eigenvalue wih posiive real par or a marginally unsable Jordan block. A more appropriae formulaion of sabiliy of a saionary soluion, which also exends o nonlinear sysems is he concep of Lyapunov sabiliy. Definiion: The saionary soluion y() = is called sable or Lyapunov sable if for every ε > here is a δ > such ha for every y wih y < δ he soluion y() wih 6
7 iniial condiion y() = y exiss for all > and saisfies y() < ε for all >. The saionary soluion is called asympoically sable if in addiion all soluions saisfy lim y =. As an applicaion of he linear resul we presen he following resul for nonlinear planar auonomous sysems. Theorem 1: Le A M (, ) and g:, g( y) = O( y ) a differeniable funcion. Consider he differenial equaion y = Ay+ g( y) (5.15) If all eigenvalues of A have negaive real pars hen he saionary soluion y = of (5.15) is asympoically sable. Proof: The proof of he sabiliy resul will be a beauiful applicaion of Gronwall-s inequaliy. If y() is a soluion of (5.15) hen i saisfies he inegral equaion () = exp + exp ( ( )) ( ) y A y A s g y s ds (5.16) Le. be a norm on R. Since all eigenvalues of A have negaive real pars here exis posiive consans a, C such ha a exp( A ) x Ce x for any x R. By he assumpio n on g here is a posiive consan K such ha g( y) Therefore, if y r, hen K y if y r. g( y) Kr y Laer we shall choose r sufficienly small such ha he consan Kr will be small. Taking he norm on boh sides of (5.16) and using hese esimaes we obain he following inegral inequaliy which holds for all such ha y rin [,]: () a + ( ) a s y Ce y CKr e y() s ds 7
8 We se z () = y () e a. Then z() saisfies he inegral inequaliy + CKr By Gronwall s inequaliy z () Cy e. z() C y CKr z( s) ds We choose r sufficienly small such ha CKr = a/. We ge If we choose y such ha y y () rfor all > concluding he proof. y () Cy e a r C, hen he above inequaliy shows in paricular ha Remark: I is clear ha heorem 1 also holds for higher dimensional auonomous sysems. For he second order differenial equaion (5.3) he rivial soluion x = is asympoically sable if de(a) = a > and -a 1 = r(a) <. In his case he energy funcional defined in (5.7) is posiive definie and sricly decreasing for nonrivial soluions. Funcionals on phase space having his propery will be called Lyapunov funcionals. The exisence of an appropriae Lyapunov funcional will imply he sabiliy resul. 3. Applicaions An example from pharmacokineics. The reacion equaion k1 k 1 A k1 A describes a basic model for sudying he dynamics of a drug in he body where A 1 and A represen he drug concenraion in issue and blood, respecively. The rae k is due o removal from he body by he kidneys. The corresponding sysem of differenial equaions is hen y = k y k y = y y k y k y k (5.17) Hence A k k 1 1 = k1 k1 k (5.18) 8
9 and r(a) <, de(a) = k 1 k >. The rivial soluion is asympoically sable: he drug will compleely disappear as one expecs. 3.. The mahemaical pendulum wih fricion. The equaion for he mahemaical pendulum wih fricion is given by 1 x + γ x = gl sin ( x) (5.19) By heorem 1 he saionary soluion x = is asympoically sable. 5.. How o solve linear equaions wih consan coefficiens We consider he differenial equaion y = Ay for a real n n marix A. While i is sraighforward o find all soluions if all eigenvalues o f A have algebraic mulipliciy one (see heorem 4.5) he compuaion of soluions in he case of muliple eigenvalues is a lile bi more involved. 1. A soluion algorihm. We describe a rouine o compue a fundamenal sysem of any linear differenial equaions wih consan coefficiens. ( sep 1) Compue he characerisic polynomial P A(σ) and compue is zeros which are he eigenvalues of A. (sep ) For simple eigenvalues λ find he corresponding eigenvecors w. consruc he soluion as in heorem 4.5, ha is se y = e λ wand ake real and imaginary par in case o f complex conjugae eigenvalues λ, λ. (sep 3)Le λ be an eigenvalue wih m ulipliciy k > 1. Firs of all, solve for corresponding eigenvecors, i.e. solve ( A λe) w= The number of linear independen soluions n 1 is equal o he dimension of he nullspace (or kernel) of A λe. The soluions of he differenial equaion are consruced as in (sep ) yielding n 1 linearly independen soluions for he eigenvalue λ. (sep 4) If n 1 = k, hen you have found all soluions for his eigenvalue. If n 1 < k, compue he nullspace of (A λe), i.e. solve 9
10 ( A λe) w= The number of linear independen soluions n is equal o he dimension of he nullspace (or kernel) of (A λe). Take he n n 1 linear independen soluions, which do no belong o he nullspace of A λe, i.e. A λe w. The corresponding soluion y() of he differenial equaion is given by y() = exp(a)w. I can be explicily compued as follows: T exp( A) w = exp( ( λe + A λe) ) w = exp( λ ) exp ( λ ) ( E A E w e λ = E+ ( λ ) ( λ ) = + A E + 1 A E +... w ( λ ) λ λ e w e A E w he las equaliy is a consequence of ( A λe) w= imaginary par in case of comp lex conjugae eigenvalues λ, λ. n ) for all n > 1. Again ake real and (sep 5) If n = k, hen you have found all soluions for his eigenvalue. If n < k, hen compue he nullspace o f (A λe) 3 and proceed as before. (sep 6) Repea he procedure for he following powers of A λe unil he dimension of he nullspace and hence he number of independen soluion of he differenial equaion equals k.. Examples. We presen a few explici compuaions of soluion in he case of muliple eigenvalues..1. Solve he Cauchy problem y = 1 1 y, y() = 1 3 (5.) If A denoes he marix in (5.) hen P A (λ) = (λ 1) 3 and herefore λ = 1 is a riple eigenvalue of A. Le {e 1, e, e 3 } denoe he sandard orhonormal basis of R 3. We compue 1
11 1 A E = 1 We see ha he nullspace of A E is spanned by e 1. Hence y 1 () = exp() e 1 solves (5.). We compue (A E). ( A E) = The nullspace of (A E) is spanned by e 1 and e. Only e yields a new soluion of he differenial equaion. According o he general procedure described above i is given by y () = ee + e A E e = ee + e e We noe ha since (A E) 3 = is nullspace is R 3. Take e 3 as he hird vecor spanning R 3. The hird linear independen soluion is y3() = ee3 + e ( A E) e3+ e ( A E) e3 = e + e 1 + e 1 We have compued he fundamenal marix 1 Y() ( y (), y (), y () ) e 1 1 = 1 3 = Since Y() = E we have Y() = exp(a). The soluion of he Cauchy problem is herefore () y = e 1 = e + 3 (5.1) Solve he Cauchy problem 11
12 y = y, y() = 1 (5.) If A denoes he marix in (5.) hen P A (λ) = λ 3 and herefore λ = is a riple eigenvalue of A. The nullspace of B is spanned by he vecors 1 1 w1 = 1, w = 1 Hence y 1 () = w 1 and y () = w solve (5.). We noe ha since A = is nullspace is R 3. Take w 3 = w = (,1,) as he hird vecor spanning R 3. The hird linear independen soluion is 1 y3() = e + Ae = = 1+ We have compued he fundamenal marix 1 1 Y() = ( y1(), y(), y3() ) = We easily check ha exp A ( 1+ = y3 y1, y3, y y1 + y3 ) = 1+ 1 The soluion of he Cauchy problem is herefore y () = 1+ = (5.3) 1
13 .3. Consider he equaions of moion for wo coupled pendulums x = ax k x y y = ay k( y x) (5.4) where a and k are posiive consans. The marix of he corresponding firs-order sysem T for he vecor x, xyy,, is given by 1 a k A = k 1 k a k i has degenerae purely imaginary eigenvalues ± i a, ± i a+ k and can be solved by he procedure described above. In his case, however, i is more convenien o decouple he equaions of moion by considering he equaions for he sum and he difference of x and y, which are given by ( y) x+ y = a x+ x y = a+ k ( x y) (5.5) which are easily solved. 13
Let ( α, β 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 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 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 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 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 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 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 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 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 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 informationAn Introduction to Malliavin calculus and its applications
An Inroducion o Malliavin calculus and is applicaions Lecure 5: Smoohness of he densiy and Hörmander s heorem David Nualar Deparmen of Mahemaics Kansas Universiy Universiy of Wyoming Summer School 214
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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 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 informationDISCRETE GRONWALL LEMMA AND APPLICATIONS
DISCRETE GRONWALL LEMMA AND APPLICATIONS JOHN M. HOLTE MAA NORTH CENTRAL SECTION MEETING AT UND 24 OCTOBER 29 Gronwall s lemma saes an inequaliy ha is useful in he heory of differenial equaions. Here is
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 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 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 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 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 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 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 informationKEY. Math 334 Midterm III Winter 2008 section 002 Instructor: Scott Glasgow
KEY Mah 334 Miderm III Winer 008 secion 00 Insrucor: Sco Glasgow Please do NOT wrie on his exam. No credi will be given for such work. Raher wrie in a blue book, or on your own paper, preferably engineering
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 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 informationCERTAIN CLASSES OF SOLUTIONS OF LAGERSTROM EQUATIONS
SARAJEVO JOURNAL OF MATHEMATICS Vol.10 (22 (2014, 67 76 DOI: 10.5644/SJM.10.1.09 CERTAIN CLASSES OF SOLUTIONS OF LAGERSTROM EQUATIONS ALMA OMERSPAHIĆ AND VAHIDIN HADŽIABDIĆ Absrac. This paper presens sufficien
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 informationOn a Fractional Stochastic Landau-Ginzburg Equation
Applied Mahemaical Sciences, Vol. 4, 1, no. 7, 317-35 On a Fracional Sochasic Landau-Ginzburg Equaion Nguyen Tien Dung Deparmen of Mahemaics, FPT Universiy 15B Pham Hung Sree, Hanoi, Vienam dungn@fp.edu.vn
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 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 informationAnn. Funct. Anal. 2 (2011), no. 2, A nnals of F unctional A nalysis ISSN: (electronic) URL:
Ann. Func. Anal. 2 2011, no. 2, 34 41 A nnals of F uncional A nalysis ISSN: 2008-8752 elecronic URL: www.emis.de/journals/afa/ CLASSIFICAION OF POSIIVE SOLUIONS OF NONLINEAR SYSEMS OF VOLERRA INEGRAL EQUAIONS
More informationHamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:
M ah 5 7 Fall 9 L ecure O c. 4, 9 ) Hamilon- J acobi Equaion: Weak S oluion We coninue he sudy of he Hamilon-Jacobi equaion: We have shown ha u + H D u) = R n, ) ; u = g R n { = }. ). In general we canno
More informationUndetermined coefficients for local fractional differential equations
Available online a www.isr-publicaions.com/jmcs J. Mah. Compuer Sci. 16 (2016), 140 146 Research Aricle Undeermined coefficiens for local fracional differenial equaions Roshdi Khalil a,, Mohammed Al Horani
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 informationKEY. Math 334 Midterm III Fall 2008 sections 001 and 003 Instructor: Scott Glasgow
KEY Mah 334 Miderm III Fall 28 secions and 3 Insrucor: Sco Glasgow Please do NOT wrie on his exam. No credi will be given for such work. Raher wrie in a blue book, or on your own paper, preferably engineering
More informationPOSITIVE SOLUTIONS OF NEUTRAL DELAY DIFFERENTIAL EQUATION
Novi Sad J. Mah. Vol. 32, No. 2, 2002, 95-108 95 POSITIVE SOLUTIONS OF NEUTRAL DELAY DIFFERENTIAL EQUATION Hajnalka Péics 1, János Karsai 2 Absrac. We consider he scalar nonauonomous neural delay differenial
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 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 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 informationStructural Dynamics and Earthquake Engineering
Srucural Dynamics and Earhquae Engineering Course 1 Inroducion. Single degree of freedom sysems: Equaions of moion, problem saemen, soluion mehods. Course noes are available for download a hp://www.c.up.ro/users/aurelsraan/
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 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 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 information8. Basic RL and RC Circuits
8. Basic L and C Circuis This chaper deals wih he soluions of he responses of L and C circuis The analysis of C and L circuis leads o a linear differenial equaion This chaper covers he following opics
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 informationMATH 5720: Gradient Methods Hung Phan, UMass Lowell October 4, 2018
MATH 5720: Gradien Mehods Hung Phan, UMass Lowell Ocober 4, 208 Descen Direcion Mehods Consider he problem min { f(x) x R n}. The general descen direcions mehod is x k+ = x k + k d k where x k is he curren
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 informationMatrix Versions of Some Refinements of the Arithmetic-Geometric Mean Inequality
Marix Versions of Some Refinemens of he Arihmeic-Geomeric Mean Inequaliy Bao Qi Feng and Andrew Tonge Absrac. We esablish marix versions of refinemens due o Alzer ], Carwrigh and Field 4], and Mercer 5]
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 informationChapter 1 Fundamental Concepts
Chaper 1 Fundamenal Conceps 1 Signals A signal is a paern of variaion of a physical quaniy, ofen as a funcion of ime (bu also space, disance, posiion, ec). These quaniies are usually he independen variables
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 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 informationPOSITIVE AND MONOTONE SYSTEMS IN A PARTIALLY ORDERED SPACE
Urainian Mahemaical Journal, Vol. 55, No. 2, 2003 POSITIVE AND MONOTONE SYSTEMS IN A PARTIALLY ORDERED SPACE A. G. Mazo UDC 517.983.27 We invesigae properies of posiive and monoone differenial sysems wih
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 information6. Stochastic calculus with jump processes
A) Trading sraegies (1/3) Marke wih d asses S = (S 1,, S d ) A rading sraegy can be modelled wih a vecor φ describing he quaniies invesed in each asse a each insan : φ = (φ 1,, φ d ) The value a of a porfolio
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 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 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 informationLecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.
Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in
More informationResearch Article Existence and Uniqueness of Periodic Solution for Nonlinear Second-Order Ordinary Differential Equations
Hindawi Publishing Corporaion Boundary Value Problems Volume 11, Aricle ID 19156, 11 pages doi:1.1155/11/19156 Research Aricle Exisence and Uniqueness of Periodic Soluion for Nonlinear Second-Order Ordinary
More informationOn Oscillation of a Generalized Logistic Equation with Several Delays
Journal of Mahemaical Analysis and Applicaions 253, 389 45 (21) doi:1.16/jmaa.2.714, available online a hp://www.idealibrary.com on On Oscillaion of a Generalized Logisic Equaion wih Several Delays Leonid
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 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 informationThe Existence, Uniqueness and Stability of Almost Periodic Solutions for Riccati Differential Equation
ISSN 1749-3889 (prin), 1749-3897 (online) Inernaional Journal of Nonlinear Science Vol.5(2008) No.1,pp.58-64 The Exisence, Uniqueness and Sailiy of Almos Periodic Soluions for Riccai Differenial Equaion
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 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 informationProperties Of Solutions To A Generalized Liénard Equation With Forcing Term
Applied Mahemaics E-Noes, 8(28), 4-44 c ISSN 67-25 Available free a mirror sies of hp://www.mah.nhu.edu.w/ amen/ Properies Of Soluions To A Generalized Liénard Equaion Wih Forcing Term Allan Kroopnick
More informationTheory of! Partial Differential Equations-I!
hp://users.wpi.edu/~grear/me61.hml! Ouline! Theory o! Parial Dierenial Equaions-I! Gréar Tryggvason! Spring 010! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and
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 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 informationModal identification of structures from roving input data by means of maximum likelihood estimation of the state space model
Modal idenificaion of srucures from roving inpu daa by means of maximum likelihood esimaion of he sae space model J. Cara, J. Juan, E. Alarcón Absrac The usual way o perform a forced vibraion es is o fix
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 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 informationSolutions from Chapter 9.1 and 9.2
Soluions from Chaper 9 and 92 Secion 9 Problem # This basically boils down o an exercise in he chain rule from calculus We are looking for soluions of he form: u( x) = f( k x c) where k x R 3 and k is
More informationTheory of! Partial Differential Equations!
hp://www.nd.edu/~gryggva/cfd-course/! Ouline! Theory o! Parial Dierenial Equaions! Gréar Tryggvason! Spring 011! Basic Properies o PDE!! Quasi-linear Firs Order Equaions! - Characerisics! - Linear and
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 informationSingle and Double Pendulum Models
Single and Double Pendulum Models Mah 596 Projec Summary Spring 2016 Jarod Har 1 Overview Differen ypes of pendulums are used o model many phenomena in various disciplines. In paricular, single and double
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 informationSolution of Integro-Differential Equations by Using ELzaki Transform
Global Journal of Mahemaical Sciences: Theory and Pracical. Volume, Number (), pp. - Inernaional Research Publicaion House hp://www.irphouse.com Soluion of Inegro-Differenial Equaions by Using ELzaki Transform
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 informationElementary Differential Equations and Boundary Value Problems
Elemenar Differenial Equaions and Boundar Value Problems Boce. & DiPrima 9 h Ediion Chaper 1: Inroducion 1006003 คณ ตศาสตร ว ศวกรรม 3 สาขาว ชาว ศวกรรมคอมพ วเตอร ป การศ กษา 1/2555 ผศ.ดร.อร ญญา ผศ.ดร.สมศ
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 informationClass Meeting # 10: Introduction to the Wave Equation
MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion
More informationLecture 20: Riccati Equations and Least Squares Feedback Control
34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he
More informationOn linear ODEs with a time singularity of the first kind and unsmooth inhomogeneity
Rachůnková e al. RESEARCH On linear ODEs wih a ime singulariy of he firs kind and unsmooh inhomogeneiy Irena Rachůnková, Svaoslav Saněk, Jana Vampolová and Ewa B Weinmüller 2* * Correspondence: ewa.weinmueller@uwien.ac.a
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 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 information