Numerical Methods for Partial Differential Equations
|
|
- Leo Lloyd
- 6 years ago
- Views:
Transcription
1 Numerical Methods for Partial Differetial Equatios CAAM 45 Sprig 005 Lecture 4 -step time-steppig methods: stability, accuracy Ruge-Kutta Methods, Istructor: Tim Warburto
2 Today Recall AB stability regios ad start up issues Group aalysis of the Leap-Frog scheme Oe-step methods Example Ruge-Kutta methods: Modified Euler Geeral family of d order RK methods Heu s 3 rd Order Method The 4 th Order Ruge-Kutta method Jameso-Schmidt-Turkel Liear Absolute Stability Regios for the d order family RK Global error aalysis for geeral -step methods (stops slightly short of a full covergece aalysis) Warig o usefuless of global error estimate Discussio o AB v. RK Embedded lower order RK schemes useful for a posteriori error estimates. CAAM 45 Sprig 005
3 Recall: AB v. AB3 v. AB4 These are the margis of absolute stability for the AB methods: Startig with the yellow AB (Euler-Forward) we see that as the order of accuracy goes up the stability regio shriks. i.e. we see that to use the higher order accurate AB scheme we are required to take more time steps. Q) how may more? CAAM 45 Sprig 005
4 AB: AB: AB3: Recall: Requiremets Startig Requiremets u0 = u( 0) u ( ) + = u + dt f u u = u( dt) u = u( 0) 0 3 u = u + dt f u f u u = u dt + 0 ( 0) u = u dt u = u dt u = u + f f + f ( ) + solutio level for start solutio levels for start 3 solutio levels for start CAAM 45 Sprig 005
5 cot So as we take higher order versio of the AB scheme we also eed to provide iitial values at more ad more levels. For a problem where we do ot kow the solutio at more tha the iitial coditio we may have to: Use AB with small dt to get the secod restart level Use AB with small dt to get the third restart level March o usig AB3 started with the three levels achieved above. AB AB AB3 CAAM 45 Sprig 005
6 Recall: Derivatio of AB Schemes The AB schemes were motivated by cosiderig the exactly time itegrated ODE: Which we approximated by usig a p th order polyomial iterpolatio of the fuctio f + u t u t f u t dt + = + t + p t t t u t u t I f u t dt + + CAAM 45 Sprig 005
7 Leap Frog Scheme We could also have started the itegral at: t + Ad used the mid poit rule: Which suggests the leap frog scheme: t t + u t u t f u t dt = + + u t u t dtf u t + u u dtf u + = + CAAM 45 Sprig 005
8 Voluteer Exercise u u dtf u + = + ) accuracy: what is the local trucatio error? ) stability: what is the maifold of absolute liear stability (try aalytically) i the u=dt*mu plae? a) what is the regio of absolute liear-stability? CAAM 45 Sprig 005
9 cot u u dtf u + = + 3) How may startig values are required? 4) Do we have covergece? 5) What is the global order of accuracy? 6) Whe is this a good method? CAAM 45 Sprig 005
10 Oe Step Methods Give the difficulties iheret i startig the higher order AB schemes we are ecouraged to look for oe-step methods which oly require to evaluate i.e. u u + u = u + dt Φ + u, t; dt Euler-Forward is a oe-step method: u = u + dtf u Φ + u, t ; dt : = f u We will cosider the oe-step Ruge-Kutta methods. For itroductory details see: A itroductio to umerical aalysis, Suli ad Mayers,. (p37) ad o Trefethe p75- Gustafsso,Kreiss ad Oliger p4- CAAM 45 Sprig 005
11 Ruge-Kutta Methods The Ruge-Kutta are a family of oe-step methods. They cosist of s stages (i.e. require s evaluatios of f) They will be p th order accurate, for some p. They are self startig!!!. CAAM 45 Sprig 005
12 Example Ruge-Kutta Method (Modified Euler) a = dtf ( u, t ) Modified Euler: a dt b = dtf u, t + + u+ = u + b Note how we oly eed oe startig value. We ca also reiterpret this through itermediate values: dt uˆ = u + f ( u, t ) u+ = u + dtf ( uˆ, t+ / ) This looks like a half step to approximate the miditerval u ad the a full step. This is a -stage, d order, sigle step method. CAAM 45 Sprig 005
13 Liear Stability Aalysis As before we assume that f is liear i u ad idepedet of time The scheme becomes (for some give mu): dt uˆ = u + f u t u = u + dtf u t (, ) ( ˆ, ) + + / dt uˆ = u + µ u u = u + dtµ uˆ + Which we simplify (elimiate the uhat variable): dt uˆ = u + µ u ( dt ) u = u + µ dtµ u + + u u = u + dtµ uˆ + CAAM 45 Sprig 005
14 cot We gather all terms o the right had side: ( dtµ ) u+ = + dtµ + u [ Note: the bracketed term is exactly the first 3 terms of the Taylor series for exp(dt*mu), more o that later ] We also ote for the umerical solutio to be bouded, ad the scheme stable, we require: ( dtµ ) + dtµ + CAAM 45 Sprig 005
15 cot The stability regio is the set of u=mu*dt i the complex plae such that: ν + ν + The maifold of margial stability ca be foud (as i the liear multistep methods) by fixig the multiplier to be of uit magitude ad lookig for the correspodig values of u which produce this multiplier. i.e. for each theta fid u such that ν + ν + = i e θ CAAM 45 Sprig 005
16 cot We ca maually fid the roots of this quadratic: ν + ν + = i e θ To obtai a parameterized represetatio of the maifold of margial stability: ( i ) ν = ± e θ CAAM 45 Sprig 005
17 Plottig Stability Regio for Modified Euler ( i ) ν = ± e θ CAAM 45 Sprig 005
18 Checkig Modified Euler at the Imagiary Axis As before we wish to check how much of the imagiary axis is icluded iside the regio of absolute stability. Here we plot the real part of the + root CAAM 45 Sprig 005
19 Is the Imagiary Axis i the Stability Regio? We ca aalytically zoom i by choosig u=i*alpha (i.e. o the imagiary axis). We the check the magitude of the multiplier: 4 ν α α α + ν + = + iα = α + = + 4 So we kow that the oly poit o the imagiary axis with multiplier magitude bouded above by is the origi. Modified Euler is ot suitable for the advectio equatio. CAAM 45 Sprig 005
20 Geeral stage RK family Cosider the four parameter family of RK schemes of the form: k = (, ) ( β, α ) f u t k = f u + dtk t + dt u = u + dt ak + bk + where we will determie the parameters (a,b,alpha,beta) by cosideratio of accuracy. [ Euler-Forward is i this family with a=,b=0 CAAM 45 Sprig 005
21 cot The sigle step operator i this case is: k = (, ) ( β, α ) (, ) β f u t k = f u + dtk t + dt u = u + dt ak + bk + u = u + Φ u t + where Φ u, t = af u, t + bf u + dtf u, t, t + αdt CAAM 45 Sprig 005
22 cot We ow perform a trucatio aalysis, similar to that performed for the liear multistep methods. We will use the followig fact: du dt = (, t) f u t d u d f f du f f (, ) dt dt t u dt t u 3 d u d f f f... 3 dt dt t u = f u t t = + = + f = + = f f f f f f f t t u u t u u t u f f f f CAAM 45 Sprig 005
23 cot (accuracy) We expad Phi i terms of powers of dt usig the bivariate Taylor s expasio ( ) β where: ( α ) Φ u t, t = af u, t + bf u + dtf u, t, t + dt f + f f = af + b αdt + βdtf + + O dt t u ( αdt) f f ( βdtf ) f ( αdt)( βdt) f + +! t t u! u (, ) f = f u t t 3 CAAM 45 Sprig 005
24 cot We costruct the local trucatio error as: ( ) (, ) T = u t + dt u t dtφ u t t dt dt = dtf + ft + ffu + ftt + ftu f + fuu f + fu ft + fu f 3! ( α dt) ( βdtf ) dtaf + b f + αdtf + βdtff + f + αβdt ff + f + O dt 4 t u tt tu uu Now we choose a,b,alpha,beta to miimize the local trucatio error. Note we use subidexig to represet partial derivatives. CAAM 45 Sprig 005
25 cot Cosider terms which are the same order i dt i the local trucatio error: dt dt T = dtf + ( ft + ffu ) + ( ftt + ftu f + fuu f + fu ( ft + fu f )) 3! αdt βdtf dtaf + b f + αdtft + βdtffu + ftt + αβdt fftu + fuu + O dt Coditio : a b = 0 4 Coditio : ( ft + ffu ) b( αdtft + βdtffu ) = 0 f bα = bβ = Uder these coditios, the trucatio is order 3 so the method is d order accurate. It is ot possible to further elimiate the dt^3 terms by adjustig the parameters. CAAM 45 Sprig 005
26 Case: No Explicit t Depedece i f ( u ( t ), t ) bf u βdtf ( u) Φ = + ( βdtf ) f f = bf + βdtf + O dt + u! u 3 du d u f d u f f ( ), 3 = f u t = f = f + f dt dt u dt u u ( ; ) T = u t u t dtφ u t dt + ( dt) ( 3 ) = dt + + ( fuu f + fu f ) dt b + + f + O d 3! dt dt β f ff 4 u f βdtff u uu t b =, β = It is easier to geeralize to higher order RK i this case whe there is o explicit time depedece i f. CAAM 45 Sprig 005
27 Secod Example Ruge-Kutta: Heu s Third Order Formula Traditioal versio (, ) a = dtf u t a dt b = dtf u, t b dt c = dtf u +, t u+ = u + ( a + 3c) 4 I terms of itermediate variables: dt uˆ = u + f ( u, t ) 3 dt uˆ = u + f ( uˆ, t+ /3 ) 3 u = u + f u t + f u t 4 ( (, ) 3 ( ˆ, )) + + /3 This is a 3 rd order, 3 stage sigle step explicit Ruge-Kutta method. CAAM 45 Sprig 005
28 Agai Let s Check the Stability Regio dt uˆ = u + f ( u, t ) 3 dt uˆ = u + f ( uˆ, t+ /3 ) 3 u = u + f u t + f u t 4 ( (, ) 3 ( ˆ, )) + + /3 With f=mu*u reduces to a sigle level recursio with a very familiar multiplier: dt uˆ = u + µ u 3 dt dt uˆ = u + µ u + µ u 3 3 dt dt dt u = u + µ u 3µ u µ u µ u dt dt dt = u + µ 3µ µ µ u ( µ dt) ( µ dt) = + µ dt u CAAM 45 Sprig 005
29 Stability of Heu s 3 rd Order Method Each margially stable mu*dt is such that the multiplier is of magitude, i.e. 3 ν ν + ν + + = 6 This traces a curve i the u=mu*dt complex plae. Sice we are short o time we ca plot this usig Matlab s roots fuctio i e θ CAAM 45 Sprig 005
30 Stability Regio for RK (s=p) rk3 rk rk4 CAAM 45 Sprig 005
31 Heu s Method ad The Imagiary Axis This time we cosider poits o the imagiary axis which are close to the origi: ν = iα 3 α α + iα i 6 3 α α = + α α α = + 36 rk3 Ad this is bouded above by if α 3.73 CAAM 45 Sprig 005
32 Observatio o RK liear stability For the s th order, s stage RK we see that the stability regio grows with icreasig s: Cosequetly we ca take a larger time step (dt) as the order of the RK scheme icrease. O the dow side, we require more evaluatios of f CAAM 45 Sprig 005
33 Popular 4 th Order Ruge-Kutta Formula Four stages: a = dtf u, t b = dtf u + a /, t c = dtf u + b /, t d = dtf u + c, t + + / + / u+ = u + a + b + c + d 6 see: of miimum umber of stages to achieve p th order. p76 for details CAAM 45 Sprig 005
34 Imagiary Axis (agai) With the obvious multiplier we obtai: ν = iα 3 4 ν ν ν + ν = α α α α α + iα i + = For stability we require: α α α i.e. α.83 rk4 CAAM 45 Sprig 005
35 Imagiary Axis Stability Summary.83 for the 4 th Order Ruge-Kutta method.73 for Heu s 3 rd Order Method 0 for modified Euler CAAM 45 Sprig 005
36 Boudig the Global Error i Terms of the Local Trucatio Error Theorem: Cosider the geeral oe-step method u = u + dt Φ + u, t; dt ad we assume that Phi is Lipschitz cotiuous with respect to the first argumet (with costat L Φ ) i.e. for (, ; ) (, ; ) { } u, t, v, t D = u, t : t t t, u u C we have: Φ u t dt Φ v t dt L u v 0 max 0 The assumig u ( t ) u ( t ) C = N it follows that 0,,.., T u u t e, = 0,,..., N where T = max T ( L Φ( t ) ) t 0 L 0 N Φ Φ CAAM 45 Sprig 005
37 cot Proof: we use the defiitio of the local trucatio error: T = ( u( t + dt) u( t )) dtφ( u( t ), t ) to costruct the error equatio: we use the Lipschitz cotiuity of Phi: tidyig: ( ) u t + u + = u t u + dt Φ u t, t Φ u, t + T u t u u t u + dtl u t u + T + + Φ u t u dtl u t u T Φ + CAAM 45 Sprig 005
38 ( ) proof cot u t u + dtl u t u + T + + Φ { } ( ) ( ) + dtl + dtl u t u + T + T + ( dtl ) u( t ) u T ( dtl ) m= T + m= 0 m Φ Φ Φ 0 0 m Φ m= 0 ( dtl ) Φ m m= { } m= max ( ) { max } m dtlφ Tm dtlφ Tm + = m m m= 0 Φ T ( + ) T ( + ) dtl + dtlφ e dtl dtl T dtl Φ Φ ( L ) Φ t+ t0 e Φ Φ m dtl CAAM 45 Sprig 005
39 proof summary We ow have the global error estimate: T u t u e ( L ) Φ t+ t0 + + dtlφ Broadly speakig this implies that if the local trucatio error is h^{p+} the the error at a give time step will scale as O(h^p): p ( + ) + u t u O h Covergece follows uder restrictios o the ODE which guaratee existace of a uique C solutio ad stable choice of dt. CAAM 45 Sprig 005
40 Warig About Global Error Estimate It should be oted that the error estimate is of almost zero practical use. T u t u e ( L ) Φ t+ t0 + + dtlφ I the full covergece aalysis we pick a fial time t ad we will see that expoetial term agai. Covergece is guarateed but the costat ca be extraordiarily large for fiite time: L Φ ( L ) Φ t t0 e CAAM 45 Sprig 005
41 A Posteriori Error Estimate There are examples of RK methods which have embedded lower order schemes. i.e. after oe full RK time step, for some versios it is possible to use a secod set of coefficiets to recostruct a lower order approximatio. Thus we ca compute the differece betwee the two differet approximatios to estimate the local trucatio error committed over the time step. google: ruge kutta embedded Numerical recipes i C: CAAM 45 Sprig 005
42 My Favorite s Stage Ruge-Kutta Method There is a s stage Ruge-Kutta method of particular simplicity due to Jameso-Schmidt- Turkel, which is of iterest whe there is o explicit time depedece for f u for m=0:s- + u dt u = u + f u s m ed u = = u CAAM 45 Sprig 005
43 RK v. AB Whe should we use RK ad whe should we use AB? rk3 rk rk4 CAAM 45 Sprig 005
44 Class Cacelled o 0/7/05 There will be o class o Thursday 0/7/05 The homework due for that class will be due the followig Thursday 0/4/05 CAAM 45 Sprig 005
Numerical Methods for Partial Differential Equations
Numerical Methods for Partial Differetial Equatios CAAM 452 Sprig 2005 Lecture 9 Istructor: Tim Warburto Today Chage of pla I will go through i detail how to solve homework 3 step by step. Q1) Build a
More information3. Z Transform. Recall that the Fourier transform (FT) of a DT signal xn [ ] is ( ) [ ] = In order for the FT to exist in the finite magnitude sense,
3. Z Trasform Referece: Etire Chapter 3 of text. Recall that the Fourier trasform (FT) of a DT sigal x [ ] is ω ( ) [ ] X e = j jω k = xe I order for the FT to exist i the fiite magitude sese, S = x [
More information6.003 Homework #3 Solutions
6.00 Homework # Solutios Problems. Complex umbers a. Evaluate the real ad imagiary parts of j j. π/ Real part = Imagiary part = 0 e Euler s formula says that j = e jπ/, so jπ/ j π/ j j = e = e. Thus the
More informationMath 113, Calculus II Winter 2007 Final Exam Solutions
Math, Calculus II Witer 7 Fial Exam Solutios (5 poits) Use the limit defiitio of the defiite itegral ad the sum formulas to compute x x + dx The check your aswer usig the Evaluatio Theorem Solutio: I this
More information1 Approximating Integrals using Taylor Polynomials
Seughee Ye Ma 8: Week 7 Nov Week 7 Summary This week, we will lear how we ca approximate itegrals usig Taylor series ad umerical methods. Topics Page Approximatig Itegrals usig Taylor Polyomials. Defiitios................................................
More informationNICK DUFRESNE. 1 1 p(x). To determine some formulas for the generating function of the Schröder numbers, r(x) = a(x) =
AN INTRODUCTION TO SCHRÖDER AND UNKNOWN NUMBERS NICK DUFRESNE Abstract. I this article we will itroduce two types of lattice paths, Schröder paths ad Ukow paths. We will examie differet properties of each,
More informationMath 312 Lecture Notes One Dimensional Maps
Math 312 Lecture Notes Oe Dimesioal Maps Warre Weckesser Departmet of Mathematics Colgate Uiversity 21-23 February 25 A Example We begi with the simplest model of populatio growth. Suppose, for example,
More informationSection A assesses the Units Numerical Analysis 1 and 2 Section B assesses the Unit Mathematics for Applied Mathematics
X0/70 NATIONAL QUALIFICATIONS 005 MONDAY, MAY.00 PM 4.00 PM APPLIED MATHEMATICS ADVANCED HIGHER Numerical Aalysis Read carefully. Calculators may be used i this paper.. Cadidates should aswer all questios.
More informationMIDTERM 3 CALCULUS 2. Monday, December 3, :15 PM to 6:45 PM. Name PRACTICE EXAM SOLUTIONS
MIDTERM 3 CALCULUS MATH 300 FALL 08 Moday, December 3, 08 5:5 PM to 6:45 PM Name PRACTICE EXAM S Please aswer all of the questios, ad show your work. You must explai your aswers to get credit. You will
More informationChapter 10: Power Series
Chapter : Power Series 57 Chapter Overview: Power Series The reaso series are part of a Calculus course is that there are fuctios which caot be itegrated. All power series, though, ca be itegrated because
More informationIntroduction to Signals and Systems, Part V: Lecture Summary
EEL33: Discrete-Time Sigals ad Systems Itroductio to Sigals ad Systems, Part V: Lecture Summary Itroductio to Sigals ad Systems, Part V: Lecture Summary So far we have oly looked at examples of o-recursive
More informationComplex Analysis Spring 2001 Homework I Solution
Complex Aalysis Sprig 2001 Homework I Solutio 1. Coway, Chapter 1, sectio 3, problem 3. Describe the set of poits satisfyig the equatio z a z + a = 2c, where c > 0 ad a R. To begi, we see from the triagle
More informationMath 210A Homework 1
Math 0A Homework Edward Burkard Exercise. a) State the defiitio of a aalytic fuctio. b) What are the relatioships betwee aalytic fuctios ad the Cauchy-Riema equatios? Solutio. a) A fuctio f : G C is called
More informationA NEW CLASS OF 2-STEP RATIONAL MULTISTEP METHODS
Jural Karya Asli Loreka Ahli Matematik Vol. No. (010) page 6-9. Jural Karya Asli Loreka Ahli Matematik A NEW CLASS OF -STEP RATIONAL MULTISTEP METHODS 1 Nazeeruddi Yaacob Teh Yua Yig Norma Alias 1 Departmet
More information1 6 = 1 6 = + Factorials and Euler s Gamma function
Royal Holloway Uiversity of Lodo Departmet of Physics Factorials ad Euler s Gamma fuctio Itroductio The is a self-cotaied part of the course dealig, essetially, with the factorial fuctio ad its geeralizatio
More informationTaylor expansion: Show that the TE of f(x)= sin(x) around. sin(x) = x - + 3! 5! L 7 & 8: MHD/ZAH
Taylor epasio: Let ƒ() be a ifiitely differetiable real fuctio. A ay poit i the eighbourhood of 0, the fuctio ƒ() ca be represeted by a power series of the followig form: X 0 f(a) f() f() ( ) f( ) ( )
More informationSOLUTIONS TO EXAM 3. Solution: Note that this defines two convergent geometric series with respective radii r 1 = 2/5 < 1 and r 2 = 1/5 < 1.
SOLUTIONS TO EXAM 3 Problem Fid the sum of the followig series 2 + ( ) 5 5 2 5 3 25 2 2 This series diverges Solutio: Note that this defies two coverget geometric series with respective radii r 2/5 < ad
More informationL 5 & 6: RelHydro/Basel. f(x)= ( ) f( ) ( ) ( ) ( ) n! 1! 2! 3! If the TE of f(x)= sin(x) around x 0 is: sin(x) = x - 3! 5!
aylor epasio: Let ƒ() be a ifiitely differetiable real fuctio. At ay poit i the eighbourhood of =0, the fuctio ca be represeted as a power series of the followig form: X 0 f(a) f() ƒ() f()= ( ) f( ) (
More informationChapter 4. Fourier Series
Chapter 4. Fourier Series At this poit we are ready to ow cosider the caoical equatios. Cosider, for eample the heat equatio u t = u, < (4.) subject to u(, ) = si, u(, t) = u(, t) =. (4.) Here,
More information4x 2. (n+1) x 3 n+1. = lim. 4x 2 n+1 n3 n. n 4x 2 = lim = 3
Exam Problems (x. Give the series (, fid the values of x for which this power series coverges. Also =0 state clearly what the radius of covergece is. We start by settig up the Ratio Test: x ( x x ( x x
More informationApply change-of-basis formula to rewrite x as a linear combination of eigenvectors v j.
Eigevalue-Eigevector Istructor: Nam Su Wag eigemcd Ay vector i real Euclidea space of dimesio ca be uiquely epressed as a liear combiatio of liearly idepedet vectors (ie, basis) g j, j,,, α g α g α g α
More informationChapter 7: The z-transform. Chih-Wei Liu
Chapter 7: The -Trasform Chih-Wei Liu Outlie Itroductio The -Trasform Properties of the Regio of Covergece Properties of the -Trasform Iversio of the -Trasform The Trasfer Fuctio Causality ad Stability
More informationMATH 10550, EXAM 3 SOLUTIONS
MATH 155, EXAM 3 SOLUTIONS 1. I fidig a approximate solutio to the equatio x 3 +x 4 = usig Newto s method with iitial approximatio x 1 = 1, what is x? Solutio. Recall that x +1 = x f(x ) f (x ). Hece,
More informationCOMM 602: Digital Signal Processing
COMM 60: Digital Sigal Processig Lecture 4 -Properties of LTIS Usig Z-Trasform -Iverse Z-Trasform Properties of LTIS Usig Z-Trasform Properties of LTIS Usig Z-Trasform -ve +ve Properties of LTIS Usig Z-Trasform
More informationPolynomials with Rational Roots that Differ by a Non-zero Constant. Generalities
Polyomials with Ratioal Roots that Differ by a No-zero Costat Philip Gibbs The problem of fidig two polyomials P(x) ad Q(x) of a give degree i a sigle variable x that have all ratioal roots ad differ by
More information9.3 The INTEGRAL TEST; p-series
Lecture 9.3 & 9.4 Math 0B Nguye of 6 Istructor s Versio 9.3 The INTEGRAL TEST; p-series I this ad the followig sectio, you will study several covergece tests that apply to series with positive terms. Note
More informationChapter 2: Numerical Methods
Chapter : Numerical Methods. Some Numerical Methods for st Order ODEs I this sectio, a summar of essetial features of umerical methods related to solutios of ordiar differetial equatios is give. I geeral,
More informationCS321. Numerical Analysis and Computing
CS Numerical Aalysis ad Computig Lecture Locatig Roots o Equatios Proessor Ju Zhag Departmet o Computer Sciece Uiversity o Ketucky Leigto KY 456-6 September 8 5 What is the Root May physical system ca
More informationZeros of Polynomials
Math 160 www.timetodare.com 4.5 4.6 Zeros of Polyomials I these sectios we will study polyomials algebraically. Most of our work will be cocered with fidig the solutios of polyomial equatios of ay degree
More informationRecurrence Relations
Recurrece Relatios Aalysis of recursive algorithms, such as: it factorial (it ) { if (==0) retur ; else retur ( * factorial(-)); } Let t be the umber of multiplicatios eeded to calculate factorial(). The
More informationME 501A Seminar in Engineering Analysis Page 1
Accurac, Stabilit ad Sstems of Equatios November 0, 07 Numerical Solutios of Ordiar Differetial Equatios Accurac, Stabilit ad Sstems of Equatios Larr Caretto Mecaical Egieerig 0AB Semiar i Egieerig Aalsis
More informationMa 530 Introduction to Power Series
Ma 530 Itroductio to Power Series Please ote that there is material o power series at Visual Calculus. Some of this material was used as part of the presetatio of the topics that follow. What is a Power
More informationCS537. Numerical Analysis and Computing
CS57 Numerical Aalysis ad Computig Lecture Locatig Roots o Equatios Proessor Ju Zhag Departmet o Computer Sciece Uiversity o Ketucky Leigto KY 456-6 Jauary 9 9 What is the Root May physical system ca be
More informationZ - Transform. It offers the techniques for digital filter design and frequency analysis of digital signals.
Z - Trasform The -trasform is a very importat tool i describig ad aalyig digital systems. It offers the techiques for digital filter desig ad frequecy aalysis of digital sigals. Defiitio of -trasform:
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.436J/15.085J Fall 2008 Lecture 19 11/17/2008 LAWS OF LARGE NUMBERS II THE STRONG LAW OF LARGE NUMBERS
MASSACHUSTTS INSTITUT OF TCHNOLOGY 6.436J/5.085J Fall 2008 Lecture 9 /7/2008 LAWS OF LARG NUMBRS II Cotets. The strog law of large umbers 2. The Cheroff boud TH STRONG LAW OF LARG NUMBRS While the weak
More informationMath 113 Exam 3 Practice
Math Exam Practice Exam 4 will cover.-., 0. ad 0.. Note that eve though. was tested i exam, questios from that sectios may also be o this exam. For practice problems o., refer to the last review. This
More informationTERMWISE DERIVATIVES OF COMPLEX FUNCTIONS
TERMWISE DERIVATIVES OF COMPLEX FUNCTIONS This writeup proves a result that has as oe cosequece that ay complex power series ca be differetiated term-by-term withi its disk of covergece The result has
More informationKinetics of Complex Reactions
Kietics of Complex Reactios by Flick Colema Departmet of Chemistry Wellesley College Wellesley MA 28 wcolema@wellesley.edu Copyright Flick Colema 996. All rights reserved. You are welcome to use this documet
More informationENGI Series Page 6-01
ENGI 3425 6 Series Page 6-01 6. Series Cotets: 6.01 Sequeces; geeral term, limits, covergece 6.02 Series; summatio otatio, covergece, divergece test 6.03 Stadard Series; telescopig series, geometric series,
More informationx a x a Lecture 2 Series (See Chapter 1 in Boas)
Lecture Series (See Chapter i Boas) A basic ad very powerful (if pedestria, recall we are lazy AD smart) way to solve ay differetial (or itegral) equatio is via a series expasio of the correspodig solutio
More informationTEACHER CERTIFICATION STUDY GUIDE
COMPETENCY 1. ALGEBRA SKILL 1.1 1.1a. ALGEBRAIC STRUCTURES Kow why the real ad complex umbers are each a field, ad that particular rigs are ot fields (e.g., itegers, polyomial rigs, matrix rigs) Algebra
More informationProblem Cosider the curve give parametrically as x = si t ad y = + cos t for» t» ß: (a) Describe the path this traverses: Where does it start (whe t =
Mathematics Summer Wilso Fial Exam August 8, ANSWERS Problem 1 (a) Fid the solutio to y +x y = e x x that satisfies y() = 5 : This is already i the form we used for a first order liear differetial equatio,
More informationThe z-transform. 7.1 Introduction. 7.2 The z-transform Derivation of the z-transform: x[n] = z n LTI system, h[n] z = re j
The -Trasform 7. Itroductio Geeralie the complex siusoidal represetatio offered by DTFT to a represetatio of complex expoetial sigals. Obtai more geeral characteristics for discrete-time LTI systems. 7.
More informationThe Method of Least Squares. To understand least squares fitting of data.
The Method of Least Squares KEY WORDS Curve fittig, least square GOAL To uderstad least squares fittig of data To uderstad the least squares solutio of icosistet systems of liear equatios 1 Motivatio Curve
More information3.2 Properties of Division 3.3 Zeros of Polynomials 3.4 Complex and Rational Zeros of Polynomials
Math 60 www.timetodare.com 3. Properties of Divisio 3.3 Zeros of Polyomials 3.4 Complex ad Ratioal Zeros of Polyomials I these sectios we will study polyomials algebraically. Most of our work will be cocered
More informationPRELIM PROBLEM SOLUTIONS
PRELIM PROBLEM SOLUTIONS THE GRAD STUDENTS + KEN Cotets. Complex Aalysis Practice Problems 2. 2. Real Aalysis Practice Problems 2. 4 3. Algebra Practice Problems 2. 8. Complex Aalysis Practice Problems
More information(c) Write, but do not evaluate, an integral expression for the volume of the solid generated when R is
Calculus BC Fial Review Name: Revised 7 EXAM Date: Tuesday, May 9 Remiders:. Put ew batteries i your calculator. Make sure your calculator is i RADIAN mode.. Get a good ight s sleep. Eat breakfast. Brig:
More informationWe are mainly going to be concerned with power series in x, such as. (x)} converges - that is, lims N n
Review of Power Series, Power Series Solutios A power series i x - a is a ifiite series of the form c (x a) =c +c (x a)+(x a) +... We also call this a power series cetered at a. Ex. (x+) is cetered at
More information6.3 Testing Series With Positive Terms
6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial
More informationIntroduction to Optimization Techniques. How to Solve Equations
Itroductio to Optimizatio Techiques How to Solve Equatios Iterative Methods of Optimizatio Iterative methods of optimizatio Solutio of the oliear equatios resultig form a optimizatio problem is usually
More informationINFINITE SEQUENCES AND SERIES
11 INFINITE SEQUENCES AND SERIES INFINITE SEQUENCES AND SERIES 11.4 The Compariso Tests I this sectio, we will lear: How to fid the value of a series by comparig it with a kow series. COMPARISON TESTS
More informationLecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting
Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would
More informationTime-Domain Representations of LTI Systems
2.1 Itroductio Objectives: 1. Impulse resposes of LTI systems 2. Liear costat-coefficiets differetial or differece equatios of LTI systems 3. Bloc diagram represetatios of LTI systems 4. State-variable
More informationPC5215 Numerical Recipes with Applications - Review Problems
PC55 Numerical Recipes with Applicatios - Review Problems Give the IEEE 754 sigle precisio bit patter (biary or he format) of the followig umbers: 0 0 05 00 0 00 Note that it has 8 bits for the epoet,
More informationNumerical Methods for Ordinary Differential Equations
Numerical Methods for Ordiary Differetial Equatios Braislav K. Nikolić Departmet of Physics ad Astroomy, Uiversity of Delaware, U.S.A. PHYS 460/660: Computatioal Methods of Physics http://www.physics.udel.edu/~bikolic/teachig/phys660/phys660.html
More informationChapter 9: Numerical Differentiation
178 Chapter 9: Numerical Differetiatio Numerical Differetiatio Formulatio of equatios for physical problems ofte ivolve derivatives (rate-of-chage quatities, such as velocity ad acceleratio). Numerical
More informationMa 4121: Introduction to Lebesgue Integration Solutions to Homework Assignment 5
Ma 42: Itroductio to Lebesgue Itegratio Solutios to Homework Assigmet 5 Prof. Wickerhauser Due Thursday, April th, 23 Please retur your solutios to the istructor by the ed of class o the due date. You
More informationlim za n n = z lim a n n.
Lecture 6 Sequeces ad Series Defiitio 1 By a sequece i a set A, we mea a mappig f : N A. It is customary to deote a sequece f by {s } where, s := f(). A sequece {z } of (complex) umbers is said to be coverget
More informationAppendix: The Laplace Transform
Appedix: The Laplace Trasform The Laplace trasform is a powerful method that ca be used to solve differetial equatio, ad other mathematical problems. Its stregth lies i the fact that it allows the trasformatio
More informationECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations
ECE-S352 Itroductio to Digital Sigal Processig Lecture 3A Direct Solutio of Differece Equatios Discrete Time Systems Described by Differece Equatios Uit impulse (sample) respose h() of a DT system allows
More informationCHAPTER 10 INFINITE SEQUENCES AND SERIES
CHAPTER 10 INFINITE SEQUENCES AND SERIES 10.1 Sequeces 10.2 Ifiite Series 10.3 The Itegral Tests 10.4 Compariso Tests 10.5 The Ratio ad Root Tests 10.6 Alteratig Series: Absolute ad Coditioal Covergece
More information1.3 Convergence Theorems of Fourier Series. k k k k. N N k 1. With this in mind, we state (without proof) the convergence of Fourier series.
.3 Covergece Theorems of Fourier Series I this sectio, we preset the covergece of Fourier series. A ifiite sum is, by defiitio, a limit of partial sums, that is, a cos( kx) b si( kx) lim a cos( kx) b si(
More informationSequences and Series of Functions
Chapter 6 Sequeces ad Series of Fuctios 6.1. Covergece of a Sequece of Fuctios Poitwise Covergece. Defiitio 6.1. Let, for each N, fuctio f : A R be defied. If, for each x A, the sequece (f (x)) coverges
More informationLecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting
Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would
More informationMath 128A: Homework 1 Solutions
Math 8A: Homework Solutios Due: Jue. Determie the limits of the followig sequeces as. a) a = +. lim a + = lim =. b) a = + ). c) a = si4 +6) +. lim a = lim = lim + ) [ + ) ] = [ e ] = e 6. Observe that
More informationx x x Using a second Taylor polynomial with remainder, find the best constant C so that for x 0,
Math Activity 9( Due with Fial Eam) Usig first ad secod Taylor polyomials with remaider, show that for, 8 Usig a secod Taylor polyomial with remaider, fid the best costat C so that for, C 9 The th Derivative
More informationLinear Differential Equations of Higher Order Basic Theory: Initial-Value Problems d y d y dy
Liear Differetial Equatios of Higher Order Basic Theory: Iitial-Value Problems d y d y dy Solve: a( ) + a ( )... a ( ) a0( ) y g( ) + + + = d d d ( ) Subject to: y( 0) = y0, y ( 0) = y,..., y ( 0) = y
More informationfrom definition we note that for sequences which are zero for n < 0, X[z] involves only negative powers of z.
We ote that for the past four examples we have expressed the -trasform both as a ratio of polyomials i ad as a ratio of polyomials i -. The questio is how does oe kow which oe to use? [] X ] from defiitio
More information( θ. sup θ Θ f X (x θ) = L. sup Pr (Λ (X) < c) = α. x : Λ (x) = sup θ H 0. sup θ Θ f X (x θ) = ) < c. NH : θ 1 = θ 2 against AH : θ 1 θ 2
82 CHAPTER 4. MAXIMUM IKEIHOOD ESTIMATION Defiitio: et X be a radom sample with joit p.m/d.f. f X x θ. The geeralised likelihood ratio test g.l.r.t. of the NH : θ H 0 agaist the alterative AH : θ H 1,
More informationMath 2784 (or 2794W) University of Connecticut
ORDERS OF GROWTH PAT SMITH Math 2784 (or 2794W) Uiversity of Coecticut Date: Mar. 2, 22. ORDERS OF GROWTH. Itroductio Gaiig a ituitive feel for the relative growth of fuctios is importat if you really
More informationA widely used display of protein shapes is based on the coordinates of the alpha carbons - - C α
Nice plottig of proteis: I A widely used display of protei shapes is based o the coordiates of the alpha carbos - - C α -s. The coordiates of the C α -s are coected by a cotiuous curve that roughly follows
More informationMath 220B Final Exam Solutions March 18, 2002
Math 0B Fial Exam Solutios March 18, 00 1. (1 poits) (a) (6 poits) Fid the Gree s fuctio for the tilted half-plae {(x 1, x ) R : x 1 + x > 0}. For x (x 1, x ), y (y 1, y ), express your Gree s fuctio G(x,
More informationNotes 8 Singularities
ECE 6382 Fall 27 David R. Jackso Notes 8 Sigularities Notes are from D. R. Wilto, Dept. of ECE Sigularity A poit s is a sigularity of the fuctio f () if the fuctio is ot aalytic at s. (The fuctio does
More informationTHE SOLUTION OF NONLINEAR EQUATIONS f( x ) = 0.
THE SOLUTION OF NONLINEAR EQUATIONS f( ) = 0. Noliear Equatio Solvers Bracketig. Graphical. Aalytical Ope Methods Bisectio False Positio (Regula-Falsi) Fied poit iteratio Newto Raphso Secat The root of
More informationHonors Calculus Homework 13 Solutions, due 12/8/5
Hoors Calculus Homework Solutios, due /8/5 Questio Let a regio R i the plae be bouded by the curves y = 5 ad = 5y y. Sketch the regio R. The two curves meet where both equatios hold at oce, so where: y
More information(b) What is the probability that a particle reaches the upper boundary n before the lower boundary m?
MATH 529 The Boudary Problem The drukard s walk (or boudary problem) is oe of the most famous problems i the theory of radom walks. Oe versio of the problem is described as follows: Suppose a particle
More informationSIDDHARTH GROUP OF INSTITUTIONS :: PUTTUR Siddharth Nagar, Narayanavanam Road QUESTION BANK (DESCRIPTIVE)
QUESTION BANK 8 SIDDHARTH GROUP OF INSTITUTIONS :: PUTTUR Siddharth Nagar, Narayaavaam Road 5758 QUESTION BANK (DESCRIPTIVE Subject with Code : (6HS6 Course & Brach: B.Tech AG Year & Sem: II-B.Tech& I-Sem
More informationThe axial dispersion model for tubular reactors at steady state can be described by the following equations: dc dz R n cn = 0 (1) (2) 1 d 2 c.
5.4 Applicatio of Perturbatio Methods to the Dispersio Model for Tubular Reactors The axial dispersio model for tubular reactors at steady state ca be described by the followig equatios: d c Pe dz z =
More informationmx bx kx F t. dt IR I LI V t, Q LQ RQ V t,
Lecture 5 omplex Variables II (Applicatios i Physics) (See hapter i Boas) To see why complex variables are so useful cosider first the (liear) mechaics of a sigle particle described by Newto s equatio
More informationRiemann Sums y = f (x)
Riema Sums Recall that we have previously discussed the area problem I its simplest form we ca state it this way: The Area Problem Let f be a cotiuous, o-egative fuctio o the closed iterval [a, b] Fid
More information6 Integers Modulo n. integer k can be written as k = qn + r, with q,r, 0 r b. So any integer.
6 Itegers Modulo I Example 2.3(e), we have defied the cogruece of two itegers a,b with respect to a modulus. Let us recall that a b (mod ) meas a b. We have proved that cogruece is a equivalece relatio
More information*X203/701* X203/701. APPLIED MATHEMATICS ADVANCED HIGHER Numerical Analysis. Read carefully
X0/70 NATIONAL QUALIFICATIONS 006 MONDAY, MAY.00 PM.00 PM APPLIED MATHEMATICS ADVANCED HIGHER Numerical Aalysis Read carefully. Calculators may be used i this paper.. Cadidates should aswer all questios.
More informationMAT 271 Project: Partial Fractions for certain rational functions
MAT 7 Project: Partial Fractios for certai ratioal fuctios Prerequisite kowledge: partial fractios from MAT 7, a very good commad of factorig ad complex umbers from Precalculus. To complete this project,
More informationAnalysis of Algorithms. Introduction. Contents
Itroductio The focus of this module is mathematical aspects of algorithms. Our mai focus is aalysis of algorithms, which meas evaluatig efficiecy of algorithms by aalytical ad mathematical methods. We
More information(a) (b) All real numbers. (c) All real numbers. (d) None. to show the. (a) 3. (b) [ 7, 1) (c) ( 7, 1) (d) At x = 7. (a) (b)
Chapter 0 Review 597. E; a ( + )( + ) + + S S + S + + + + + + S lim + l. D; a diverges by the Itegral l k Test sice d lim [(l ) ], so k l ( ) does ot coverge absolutely. But it coverges by the Alteratig
More informationAssignment 1 : Real Numbers, Sequences. for n 1. Show that (x n ) converges. Further, by observing that x n+2 + x n+1
Assigmet : Real Numbers, Sequeces. Let A be a o-empty subset of R ad α R. Show that α = supa if ad oly if α is ot a upper boud of A but α + is a upper boud of A for every N. 2. Let y (, ) ad x (, ). Evaluate
More informationInfinite Sequences and Series
Chapter 6 Ifiite Sequeces ad Series 6.1 Ifiite Sequeces 6.1.1 Elemetary Cocepts Simply speakig, a sequece is a ordered list of umbers writte: {a 1, a 2, a 3,...a, a +1,...} where the elemets a i represet
More informationUnit 4: Polynomial and Rational Functions
48 Uit 4: Polyomial ad Ratioal Fuctios Polyomial Fuctios A polyomial fuctio y px ( ) is a fuctio of the form p( x) ax + a x + a x +... + ax + ax+ a 1 1 1 0 where a, a 1,..., a, a1, a0are real costats ad
More informationA collocation method for singular integral equations with cosecant kernel via Semi-trigonometric interpolation
Iteratioal Joural of Mathematics Research. ISSN 0976-5840 Volume 9 Number 1 (017) pp. 45-51 Iteratioal Research Publicatio House http://www.irphouse.com A collocatio method for sigular itegral equatios
More informationSeunghee Ye Ma 8: Week 5 Oct 28
Week 5 Summary I Sectio, we go over the Mea Value Theorem ad its applicatios. I Sectio 2, we will recap what we have covered so far this term. Topics Page Mea Value Theorem. Applicatios of the Mea Value
More informationb i u x i U a i j u x i u x j
M ath 5 2 7 Fall 2 0 0 9 L ecture 1 9 N ov. 1 6, 2 0 0 9 ) S ecod- Order Elliptic Equatios: Weak S olutios 1. Defiitios. I this ad the followig two lectures we will study the boudary value problem Here
More informationPolynomial Functions and Their Graphs
Polyomial Fuctios ad Their Graphs I this sectio we begi the study of fuctios defied by polyomial expressios. Polyomial ad ratioal fuctios are the most commo fuctios used to model data, ad are used extesively
More informationIn algebra one spends much time finding common denominators and thus simplifying rational expressions. For example:
74 The Method of Partial Fractios I algebra oe speds much time fidig commo deomiators ad thus simplifyig ratioal epressios For eample: + + + 6 5 + = + = = + + + + + ( )( ) 5 It may the seem odd to be watig
More informationChapter 7 z-transform
Chapter 7 -Trasform Itroductio Trasform Uilateral Trasform Properties Uilateral Trasform Iversio of Uilateral Trasform Determiig the Frequecy Respose from Poles ad Zeros Itroductio Role i Discrete-Time
More informationDefinition of z-transform.
- Trasforms Frequecy domai represetatios of discretetime sigals ad LTI discrete-time systems are made possible with the use of DTFT. However ot all discrete-time sigals e.g. uit step sequece are guarateed
More informationMath 21B-B - Homework Set 2
Math B-B - Homework Set Sectio 5.:. a) lim P k= c k c k ) x k, where P is a partitio of [, 5. x x ) dx b) lim P k= 4 ck x k, where P is a partitio of [,. 4 x dx c) lim P k= ta c k ) x k, where P is a partitio
More informationMath 61CM - Solutions to homework 3
Math 6CM - Solutios to homework 3 Cédric De Groote October 2 th, 208 Problem : Let F be a field, m 0 a fixed oegative iteger ad let V = {a 0 + a x + + a m x m a 0,, a m F} be the vector space cosistig
More informationMTH Assignment 1 : Real Numbers, Sequences
MTH -26 Assigmet : Real Numbers, Sequeces. Fid the supremum of the set { m m+ : N, m Z}. 2. Let A be a o-empty subset of R ad α R. Show that α = supa if ad oly if α is ot a upper boud of A but α + is a
More informationj=1 dz Res(f, z j ) = 1 d k 1 dz k 1 (z c)k f(z) Res(f, c) = lim z c (k 1)! Res g, c = f(c) g (c)
Problem. Compute the itegrals C r d for Z, where C r = ad r >. Recall that C r has the couter-clockwise orietatio. Solutio: We will use the idue Theorem to solve this oe. We could istead use other (perhaps
More informationn 3 ln n n ln n is convergent by p-series for p = 2 > 1. n2 Therefore we can apply Limit Comparison Test to determine lutely convergent.
06 微甲 0-04 06-0 班期中考解答和評分標準. ( poits) Determie whether the series is absolutely coverget, coditioally coverget, or diverget. Please state the tests which you use. (a) ( poits) (b) ( poits) (c) ( poits)
More information