wiki:u/erica/truncerror

Version 3 (modified by Erica Kaminski, 11 years ago) ( diff )

Taylor series and finite differences

To formally analyse the error associated with a finite difference approximation for a derivative, one expands each possible term of the approximation in a Taylor series (TS) so that Dapprox = Taylor series. You then would take the difference between the approximation (LHS) and the derivative term from the TS (on RHS), and you are left with the 'error'. As h gets small, the biggest term of this series will dominate the error. You now have an explicit coefficient for E ~ O(hn).

Replacing derivatives in differential equations with finite difference approximations replaces the original diff-eq problem with a system of algebraic equations. To analyse the error in this case, we need the concepts of "local truncation error" at each grid point, as well as some handle on the "global error" over the entire mesh.

Local truncation error (LTE)

LTE is a measure of how well a solution to a particular ODE satisfies the finite difference approximation of that ODE. Say we discretize the following equation:

The concept of LTE is illustrated with the following example. Suppose we replace the 2nd derivative of the 1D Poisson equation,

with the centered (2nd order, as can be verified by above discussion) finite different approximation,

(which is easily interpreted as the finite difference version of the derivative of the derivative). With these, and replacing the functions with discrete versions of the functions, we have

and discretize the functions over some finite number of grid points:

with the finite difference version

Note: See TracWiki for help on using the wiki.