5.6 Multistep Methods (Cont'd) Example. Derive Adams-Bashforth Two-Step Explicit Method: Solve The IVP
5.6 Multistep Methods (Cont'd) Example. Derive Adams-Bashforth Two-Step Explicit Method: Solve The IVP
Derive Adams-Bashforth two-step explicit method: Solve the IVP: ( ) over Integrate ) ( Use ( approximate ( ). ( ) ( ( ))
( )
) where and the backward difference ( Consequently, Adams-Bashforth two-step explicit method is: ( Local Truncation Error. If ( ) solves the IVP ( ( the local truncation error is ( ) ( ) ( ) ( ) ( ) ) ) ) ( ( ) ) ( ) )
( ) )
and
))
( ))
NOTE: the local truncation error of a m-step explicit step is ( the local truncation error of a m-step implicit step is (
Predictor-Corrector Method Motivation: (1) Solve the IVP Solution: The three-step Adams-Moulton method is ( ) by the three-step Adams-Moulton method. ( ) ( ) can be solved by Newtons method. However, this can be quite computationally expensive. (2) combine explicit and implicit methods. 4th order Predictor-Corrector Method (we will combine 4th order Runge-Kutta method + 4th order 4-step explicit Adams-Bashforth method + 4th order three-step Adams-Moulton implicit method) Step 1: Use 4th order Runge-Kutta method to compute Step 2: For (a) Predictor sub-step. Use 4th order 4-step explicit Adams-Bashforth method to compute a predicated value ( ) ( ) ( ) ( )
(b) Correction sub-step. Use 4th order three-step Adams-Moulton implicit method to compute a correction (the approximation at time step) ( ) ( ) ( ) ( )
5.10 Stability Consistency and Convergence Definition. A one-step difference equation with local truncation error ( ) ( ) is said to be consistent if
Definition. A one-step difference equation is said to be convergent if ( ) where ( ) is the exact solution and is the approximate solution.
( ) ( ) ( ) ( ) be continuous and satisfy a Lipschitz Example. To solve . Let condition with Lipschitz constant . Show that Eulers method is consistent and convergent. Solution: ( ) ( ) Thus Eulers method is consistent. By Theorem 5.9, ( ) ( ) Thus Eulers method is convergent. The rate of convergence of Eulers method is ( ) [
( ( ) )
( )
Stability Motivation: How does round-off error affect approximation? To solve IVP method. Suppose is the round-off error associated with each step. ( Then ( ) ( )[
( )
( )
by Eulers
)
( )
. Here )
Stability: small changes in the initial conditions produce correspondingly small changes in the subsequent approximations. Convergence of One-Step Methods ( Theorem. Suppose the IVP form
( )
exists and ( {(
with constant
on
then
(1) The method is stable; (2) The method is convergent if and only if it is consistent: ( ) ( ) (3) If exists s.t. ( ) ( ) when then ( )
( )