Chapter 3 Boundary Value Problem
|
|
- Scott Barker
- 5 years ago
- Views:
Transcription
1 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 us consider a genearal ODE of he form x (n) = f(, x, x, x,,x (n 1) ), [a, b] (3.1) A = a and = b he soluion is supposed o saisfy r 1 (x(a)x (a),,x (n 1) (a), x(b)x (b),,x (n 1) (b)) = 0, (3.2). r n (x(a)x (a),,x (n 1) (a), x(b)x (b),,x (n 1) (b)) = 0. The resuling problem (3.1) (3.2) is called a wo poin boundary value problem [8]. In order o be useful in applicaions, a BVP (3.1) (3.2) should be well posed. This means ha given he inpu o he problem here exiss a unique soluion, which depends coninuously on he inpu. However, quesions of exisence and uniqueness for BVPs are much more difficul han for IVPs and here is no general heory. 3.1 Single shooing mehods Linear shooing mehod Consider a linear wo-poin second-order BVP of he form x () = p()x ()+q()x()+r(), [a, b] (3.3) wih x(a) = α, x(b) = β. 37
2 The main idea of he mehod is o reduce he soluion of he BVP (3.3) o he soluion of an iniial value problem [10, 6]. Namely, le us consider wo special IVPs for wo funcions u() and v(). Suppose ha u() is a soluion of he IVP u () = p()u ()+q()u()+r(), u(a) = α, u (a) = 0 and v() is he unique soluion o he IVP v () = p()v ()+q()v(), v(a) = 0, v (a) = 1. Then he linear combinaion x() = u()+cv(), c = cons. (3.4) is a soluion o BVP (3.3). The unknown consan c can be found from he boundary condiion on he righ end of he ime inerval, i.e., x(b) = u(b)+cv(b) = β c = β u(b). v(b) Tha is, if v(b) 0 he unique soluion of (3.3) reads x() = u()+ β u(b) v(). v(b) Example 1 Le us solve a BVP [6] x () = x () 2 x()+1, (3.5) 1+ 2 x(0) = 1.25, x(1) = over he ime inerval [0, 4] using he linear shooing mehod (3.4). According o Eq. (3.4) he soluion of his equaion has he form x() = u() 0.95+u(4) v(), v(4) where u() and v() are soluions of wo IVPs and u () = v () = u () u()+1, u(0) = 1.25, u (0) = v () v(), v(0) = 0, v (0) = 1.
3 4 3 u() v() x() 2 1 u,v,x 0 1 Fig. 3.1 Numerical soluion of Eq. (3.5) over he inerval [0, 4] by he linear shooing mehod (3.4) Numerical soluion of he problem 3.5 as well as boh duncions u() and v() are presened on Fig Single shooing for general BVP For a general BVP for a second-order ODE, he simple shooing mehod is saed as follows: Le x () = f(, x(), x ()), [a, b] (3.6) x(a) = α, x(b) = β. be he BVP in quesion and le x(, s) denoe he soluion of he IVP x () = f(, x(), x ()), [a, b] (3.7) x(a) = α, x (a) = s, where s is a parameer ha can be varied. The IVP (3.7) is solved wih differen values of s wih, e.g., RK4 mehod ill he boundary condiion on he righ side x(b) = β becomes fulfilled. As menioned above, he soluion x(, s) of (3.7) depends on he parameer s. Le us define a funcion F(s) := x(b, s) β. If he BVP (3.6) has a soluion, hen he funcion F(s) has a roo, which is jus he value of he slope x (a) giving he soluion x() of he BVP in quesion. The zeros of F(s) can be found wih, e.g., Newon s mehod [7]. The Newon s mehod is probably he bes known mehod for finding numerical approximaions o he zeroes of a real-valued funcion. The idea of he mehod is
4 o use he firs few erms of he Taylor series of a funcion F(s) in he viciniy of a suspeced roo, i.e., F(s n + h) = F(s n )+F (s n )h+o(h 2 ). where s n is a n h approximaion of he roo. Now if one insers h = s s n, one obains F(s) = F(s n )+F (s n )(s s n ). As he nex approximaion s n+1 o he roo we choose he zero of his funcion, i.e., F(s n+1 ) = F(s n )+F (s n )(s n+1 s n ) = 0 s n+1 = s n F(s n) F (s n ). (3.8) The derivaive F (s n ) can be calculaed using he forward difference formula F (s n ) = F(s n + δ s) F(s n ) δ s where δ s is small. Noice ha his procedure can be unsable near a horizonal asympoe or a local exremum. Example 1 Consider a simple nonlinear BVP [10] x () = 3 2 x()2, (3.9) x(0) = 4, x(1) = 1 over he inerval [0, 1] and le us solve i numerically wih he single shooing mehod discussed above. Firs of all we define a corresponding IVP x () = 3 2 x()2 x(0) = 4, x (0) = s over [0, 1] and solve i for differen values of s, e.g., s [ 100, 0] wih he classical RK4 mehod. The resul of calculaion is presened on Fig. 3.2 (a). One can see, ha he funcion F(s) = x(, s) 1 admis wo zeros, depiced on Fig. 3.2 (a) as green poins. In order o find hem we use he Newon s mehod, discussed above. The mehod gives an approximaion o boh zeros of he funcion F(s): s = { 35.8, 8.0}, which give he righ slope x (0). Boh soluions, corresponding o wo differen values of s are presened on Fig. 3.2 (b).
5 (a) (b) F(,s) x Fig. 3.2 Numerical soluion of BVP (3.9) wih single shooing mehod. (a) The Funcion F(s) = x(, s) 1 is presened. Green poins depic wo zeros of his funcion, which can be found wih Newon s mehod. (b) Two soluions of (3.9) corresponding o wo differen values of parameer s (he red line corresponds o s = 35.8, whereas he blue one o s = 8.0). Example 2 Le us consider a linear eigenvalue problem of he form x + λ x = 0, x(0) = x(1) = 0, x (0) = 1 (3.10) over [0, 1] wih he simple shooing mehod. The exac soluion is λ = n 2 π 2, n N. In order o apply he simple shooing mehod we consider a corresponding IVP of he firs order wih addiional equaion for he unknown funcion λ(): wih x = y, y = λ x, λ = 0 x(0) = 0, x (0) = 1, λ(0) = s. where s is a free shooing parameer. Here we choose s = {0.5, 50, 100}. Resuls of he shooing wih hese iniial parameers are shown on Fig One can see, ha numerical soluions correspond o firs hree eigenvalues λ = {π 2,(2π) 2,(3π) 2 }. Example 3 Consider a nonlinear BVP of he fourh order [8]
6 λ=0.5 λ=50 λ= x Fig. 3.3 Numerical soluions of Eq. (3.10) over he inerval [0, 1] by single shooing mehod. Firs hree eigenfuncions, corresponding o eigenvalues λ = {π 2,(2π) 2,(3π) 2 } are presened wih x (4) () (1+ 2 )x () 2 + 5x() 2 = 0, [0, 1] (3.11) x(0) = 1, x (0) = 0, x (1) = 2, x (1) = 3. Our goal is o solve his equaion wih he simple shooing mehod. To his end, firs we rewrie he equaion as a sysem of four ODE s of he firs order: x 1 = x 2, x 2 = x 3, x 1 (0) = 1, x 3 (1) = 2, x 3 = x 4, x 2 (0) = 0, x 4 (1) = 3, x 4 = (1+ 2 )x 2 3 5x2 1. As he second sep we consider correspondig IVP x 1 = x 2, x 2 = x 3, x 1 (0) = 1, x 3 (1) = p, x 3 = x 4, x 2 (0) = 0, x 4 (1) = q, x 4 = (1+ 2 )x 2 3 5x2 1 wih wo free shooing parameers p and q. The soluion of his IVP fulfilles following wo requiremens: F 1 (p, q) : = x 3 (1, p, q)+2 = 0, F 2 (p, q) : = x 4 (1, p, q)+3 = 0. Tha is, a sysem of nonlinear algebraic equaions should be solved o find (p, q). The zeros of he sysem can be found wih he Newon s mehod (3.8). In his case he ieraion sep reads
7 x 0.92 Fig. 3.4 Numerical soluions of (3.11) over he inerval [0, 1] by single shooing mehod. Parameers are: p = q = 0.05, he ime sep h = 0.025, iniial shooing parameers (p 0,q 0 ) = (0, 0) where s = (p, q) T, F = (F 1, F 2 ) T and is a Jacobian of he sysem and s i+1 = s i F(s i) DF(s i ) DF(s i ) = ( F1 F 1 p q F 2 F 2 p q ) F i p = F i(p+ p, q) F i (p, q), p F i q = F i(p, q+ q) F i (p, q), q where i = 1, 2 and p, q are given values. Numerical soluion of he problem in quesion is presened on Fig Finie difference Mehod One way o solve a given BVP over he ime inerval [a, b] numerically is o approximae he problem in quesion by finie differences [8, 10, 6]. We form a pariion of he domain [a, b] using mesh poins a = 0, 1,..., N = b, where i = a+ih, h = b a, i = 0, 1,... N. N Difference quoien approximaions for derivaives can be used o solve BVP in quesion [10, 6]. In paricular, using a Taylor expansion in he viciniy of he poin
8 j, for he firs derivaive one obains a forward difference x ( i ) = x( i+1) x( i ) h In a similar way one ges a backward difference x ( i ) = x( i) x( i 1 ) h +O(h). (3.12) +O(h). (3.13) We can combine hese wo approaches and derive a cenral difference, which yields a more accurae approximaion: x ( i ) = x( i+1) x( i 1 ) +O(h 2 ). (3.14) 2h The second derivaive x ( i ) can be found in he same way using he linear combinaion of differen Taylor expansions. For example, a cenral difference reads x ( i ) = x( i+1) 2x( i )+x( i 1 ) h 2 +O(h 2 ). (3.15) Finie Difference for linear BVP Le us consider a linear BVP of he second order (3.3) x = p()x ()+q()x()+r(), [a, b], x(a) = α, x(b) = β. and inroduce he noaion x( i ) = x i, p( i ) = p i, q( i ) = q i and r( i ) = r i. Then, using Eq. (3.14) and Eq. (3.15) one can rewrie Eq. (3.3) as a difference equaion x 0 = α, x i+1 2x i + x i 1 x i+1 x i 1 h 2 = p i 2h x N = β. + q i x i + r i, i = 1,..., N 1, Now we can muliply boh sides of he second equaion wih h 2 and collec erms, involving x i 1, x i and x i+1. As resul we ge a sysem of linear equaions ( 1+ h ) ( 2 p i x i 1 (2+h 2 q i )x i + 1 h ) 2 p i x i+1 = h 2 r i, i = 1, 2,... N 1. or, in marix noaion or, more precisely Ax = b, (3.16)
9 (2+h 2 q 1 ) 1 h 2 p h 2 p 2 (2+h 2 q 2 ) 1 h 2 p h 2 p 3 (2+h 2 q 3 ) 1 h 2 p h 2 p N h p N 1 (2+h 2 q N 1 ) where ( ) ( h γ 1 = α 2 p 1 + 1, γ N = β 1 h ) 2 p N 1. Our goal is o find unknown vecor x. To his end we should inver he marix A. This marix has a band srucure and is ridiagonal. For marices of his kind a ridiagonal marix algorihm (TDMA), also known als Thomas algorihm can be used (see Appendix A for deails). x 1 x 2 x 3 x N 1 h 2 r 1 γ 1 h 2 r 2 = h 2 r 3, h 2 r N 1 γ N Example Solve a linear BVP [8] x () (1+ 2 )x() = 1, (3.17) x( 1) = x(1) = 0 over [ 1, 1] wih finie difference mehod. Firs we inroduce discree se of nodes i = 1 + ih wih given ime sep h. According o noaions used in previous secion, p() = 0, q() = (1+ 2 ), r() = 1, α = β = 0. Hence, he linear sysem (3.16) we are ineresed in reads (2+h 2 q 1 ) x 1 h 2 1 (2+h 2 q 2 ) x 2 h (2+h 2 q 3 ) x = h (2+h 2 q N 1 ) x N 1 h 2 The numerical soluion of he problem in quesion is presened on Fig Finie difference for linear eigenvalue problems Consider a Surm-Liouville problem of he form x ()+q()x() = λ v()x(), (3.18) over [a, b] wih
10 x Fig. 3.5 Numerical soluions of (3.17) over he inerval [ 1, 1] by finie difference mehod x(a) = 0, x(b) = 0. Inroducing noaion x i := x( i ), q i := q( i ), v i := v( i ), we can wrie down a difference equaion for Eq. (3.18) x 0 = 0 x i+1 2x i + x i 1 h 2 + q i x i λ v i x i = 0, i = 1,... N 1, x N = 0. If v i 0 for all i we can rewrie he difference equaion above as an eigenvalue problem (A λ I)x = 0 (3.19) for a ridiagonal marix A 2 + q h v 1 v h 2 v q h 2 v 2 h v 2 v h 2 v A = 0 + q h 2 v 3 h v 3 v h 2 v h 2 v N q h 2 v N 1 h 2 3 v N 1 v N 1 and a vecor x = (x 1, x 2,..., x N 1 ) T.
Chapter 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 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 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 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 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 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 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 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 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 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 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 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 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 informationVariational Iteration Method for Solving System of Fractional Order Ordinary Differential Equations
IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 1, Issue 6 Ver. II (Nov - Dec. 214), PP 48-54 Variaional Ieraion Mehod for Solving Sysem of Fracional Order Ordinary Differenial
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 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 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 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 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 informationOrdinary Differential Equations
Lecure 22 Ordinary Differenial Equaions Course Coordinaor: Dr. Suresh A. Karha, Associae Professor, Deparmen of Civil Engineering, IIT Guwahai. In naure, mos of he phenomena ha can be mahemaically described
More informationThen. 1 The eigenvalues of A are inside R = n i=1 R i. 2 Union of any k circles not intersecting the other (n k)
Ger sgorin Circle Chaper 9 Approimaing Eigenvalues Per-Olof Persson persson@berkeley.edu Deparmen of Mahemaics Universiy of California, Berkeley Mah 128B Numerical Analysis (Ger sgorin Circle) Le A be
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 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 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 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 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 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 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 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 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 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 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 informationEchocardiography Project and Finite Fourier Series
Echocardiography Projec and Finie Fourier Series 1 U M An echocardiagram is a plo of how a porion of he hear moves as he funcion of ime over he one or more hearbea cycles If he hearbea repeas iself every
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 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 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 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 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 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 information14 Autoregressive Moving Average Models
14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class
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 informationRecursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems
8 Froniers in Signal Processing, Vol. 1, No. 1, July 217 hps://dx.doi.org/1.2266/fsp.217.112 Recursive Leas-Squares Fixed-Inerval Smooher Using Covariance Informaion based on Innovaion Approach in Linear
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 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 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 informationOpen loop vs Closed Loop. Example: Open Loop. Example: Feedforward Control. Advanced Control I
Open loop vs Closed Loop Advanced I Moor Command Movemen Overview Open Loop vs Closed Loop Some examples Useful Open Loop lers Dynamical sysems CPG (biologically inspired ), Force Fields Feedback conrol
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 informationThe following report makes use of the process from Chapter 2 in Dr. Cumming s thesis.
Zaleski 1 Joseph Zaleski Mah 451H Final Repor Conformal Mapping Mehods and ZST Hele Shaw Flow Inroducion The Hele Shaw problem has been sudied using linear sabiliy analysis and numerical mehods, bu a novel
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 informationGENERALIZATION OF THE FORMULA OF FAA DI BRUNO FOR A COMPOSITE FUNCTION WITH A VECTOR ARGUMENT
Inerna J Mah & Mah Sci Vol 4, No 7 000) 48 49 S0670000970 Hindawi Publishing Corp GENERALIZATION OF THE FORMULA OF FAA DI BRUNO FOR A COMPOSITE FUNCTION WITH A VECTOR ARGUMENT RUMEN L MISHKOV Received
More informationComparison between the Discrete and Continuous Time Models
Comparison beween e Discree and Coninuous Time Models D. Sulsky June 21, 2012 1 Discree o Coninuous Recall e discree ime model Î = AIS Ŝ = S Î. Tese equaions ell us ow e populaion canges from one day o
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 informationIn this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should
Cambridge Universiy Press 978--36-60033-7 Cambridge Inernaional AS and A Level Mahemaics: Mechanics Coursebook Excerp More Informaion Chaper The moion of projeciles In his chaper he model of free moion
More informationTHE BELLMAN PRINCIPLE OF OPTIMALITY
THE BELLMAN PRINCIPLE OF OPTIMALITY IOANID ROSU As I undersand, here are wo approaches o dynamic opimizaion: he Ponrjagin Hamilonian) approach, and he Bellman approach. I saw several clear discussions
More informationMathcad Lecture #8 In-class Worksheet Curve Fitting and Interpolation
Mahcad Lecure #8 In-class Workshee Curve Fiing and Inerpolaion A he end of his lecure, you will be able o: explain he difference beween curve fiing and inerpolaion decide wheher curve fiing or inerpolaion
More information= ( ) ) or a system of differential equations with continuous parametrization (T = R
XIII. DIFFERENCE AND DIFFERENTIAL EQUATIONS Ofen funcions, or a sysem of funcion, are paramerized in erms of some variable, usually denoed as and inerpreed as ime. The variable is wrien as a funcion of
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 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 informationFractional Method of Characteristics for Fractional Partial Differential Equations
Fracional Mehod of Characerisics for Fracional Parial Differenial Equaions Guo-cheng Wu* Modern Teile Insiue, Donghua Universiy, 188 Yan-an ilu Road, Shanghai 51, PR China Absrac The mehod of characerisics
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 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 informationApplication of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing
Applicaion of a Sochasic-Fuzzy Approach o Modeling Opimal Discree Time Dynamical Sysems by Using Large Scale Daa Processing AA WALASZE-BABISZEWSA Deparmen of Compuer Engineering Opole Universiy of Technology
More information( ) a system of differential equations with continuous parametrization ( T = R + These look like, respectively:
XIII. DIFFERENCE AND DIFFERENTIAL EQUATIONS Ofen funcions, or a sysem of funcion, are paramerized in erms of some variable, usually denoed as and inerpreed as ime. The variable is wrien as a funcion of
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 informationPENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD
PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD HAN XIAO 1. Penalized Leas Squares Lasso solves he following opimizaion problem, ˆβ lasso = arg max β R p+1 1 N y i β 0 N x ij β j β j (1.1) for some 0.
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 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 informationMATH 4330/5330, Fourier Analysis Section 6, Proof of Fourier s Theorem for Pointwise Convergence
MATH 433/533, Fourier Analysis Secion 6, Proof of Fourier s Theorem for Poinwise Convergence Firs, some commens abou inegraing periodic funcions. If g is a periodic funcion, g(x + ) g(x) for all real x,
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 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 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 informationSOLUTIONS TO ECE 3084
SOLUTIONS TO ECE 384 PROBLEM 2.. For each sysem below, specify wheher or no i is: (i) memoryless; (ii) causal; (iii) inverible; (iv) linear; (v) ime invarian; Explain your reasoning. If he propery is no
More informationOrdinary differential equations. Phys 750 Lecture 7
Ordinary differenial equaions Phys 750 Lecure 7 Ordinary Differenial Equaions Mos physical laws are expressed as differenial equaions These come in hree flavours: iniial-value problems boundary-value problems
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 informationGuest Lectures for Dr. MacFarlane s EE3350 Part Deux
Gues Lecures for Dr. MacFarlane s EE3350 Par Deux Michael Plane Mon., 08-30-2010 Wrie name in corner. Poin ou his is a review, so I will go faser. Remind hem o go lisen o online lecure abou geing an A
More informationIntroduction to Probability and Statistics Slides 4 Chapter 4
Inroducion o Probabiliy and Saisics Slides 4 Chaper 4 Ammar M. Sarhan, asarhan@mahsa.dal.ca Deparmen of Mahemaics and Saisics, Dalhousie Universiy Fall Semeser 8 Dr. Ammar Sarhan Chaper 4 Coninuous Random
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 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 informationarxiv: v1 [math.ca] 15 Nov 2016
arxiv:6.599v [mah.ca] 5 Nov 26 Counerexamples on Jumarie s hree basic fracional calculus formulae for non-differeniable coninuous funcions Cheng-shi Liu Deparmen of Mahemaics Norheas Peroleum Universiy
More informationRobust estimation based on the first- and third-moment restrictions of the power transformation model
h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,
More informationEconomics 8105 Macroeconomic Theory Recitation 6
Economics 8105 Macroeconomic Theory Reciaion 6 Conor Ryan Ocober 11h, 2016 Ouline: Opimal Taxaion wih Governmen Invesmen 1 Governmen Expendiure in Producion In hese noes we will examine a model in which
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 informationDistance Between Two Ellipses in 3D
Disance Beween Two Ellipses in 3D David Eberly Magic Sofware 6006 Meadow Run Cour Chapel Hill, NC 27516 eberly@magic-sofware.com 1 Inroducion An ellipse in 3D is represened by a cener C, uni lengh axes
More informationMost Probable Phase Portraits of Stochastic Differential Equations and Its Numerical Simulation
Mos Probable Phase Porrais of Sochasic Differenial Equaions and Is Numerical Simulaion Bing Yang, Zhu Zeng and Ling Wang 3 School of Mahemaics and Saisics, Huazhong Universiy of Science and Technology,
More informationVehicle Arrival Models : Headway
Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where
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 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 informationTHE BERNOULLI NUMBERS. t k. = lim. = lim = 1, d t B 1 = lim. 1+e t te t = lim t 0 (e t 1) 2. = lim = 1 2.
THE BERNOULLI NUMBERS The Bernoulli numbers are defined here by he exponenial generaing funcion ( e The firs one is easy o compue: (2 and (3 B 0 lim 0 e lim, 0 e ( d B lim 0 d e +e e lim 0 (e 2 lim 0 2(e
More informationApproximation Algorithms for Unique Games via Orthogonal Separators
Approximaion Algorihms for Unique Games via Orhogonal Separaors Lecure noes by Konsanin Makarychev. Lecure noes are based on he papers [CMM06a, CMM06b, LM4]. Unique Games In hese lecure noes, we define
More informationSliding Mode Extremum Seeking Control for Linear Quadratic Dynamic Game
Sliding Mode Exremum Seeking Conrol for Linear Quadraic Dynamic Game Yaodong Pan and Ümi Özgüner ITS Research Group, AIST Tsukuba Eas Namiki --, Tsukuba-shi,Ibaraki-ken 5-856, Japan e-mail: pan.yaodong@ais.go.jp
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 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 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 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 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 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 informationUnsteady Flow Problems
School of Mechanical Aerospace and Civil Engineering Unseady Flow Problems T. J. Craf George Begg Building, C41 TPFE MSc CFD-1 Reading: J. Ferziger, M. Peric, Compuaional Mehods for Fluid Dynamics H.K.
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 informationT L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB
Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal
More information