SYLLABUS Previous: 1.4 Numerical discretization
Up: 1.4 Numerical discretization
Next: 1.4.2 Sampling on a
One of the most important properties of a numerical approximation based on a finite number of data points is that the approximation converges to the exact value as more information is gathered (exercise 1.01). This convergence can be defined locally for any arbitrary point in space and time (more restrictive) or by monitoring a global quantity (more permissive).
The convergence rate can be estimated experimentally from a geometric sequence of approximations , by successively refining the numerical resolution , for example by doubling the resolution a couple of times when
Knowing the convergence rate , it is possible to further improve the accuracy of the numerical approximation by using a so-called Richardson extrapolation
where a more accurate solution is obtained from an evaluation and its refinement .
Let's take an example and approximate the function on a grid in , simply by sampling the value of the function in the middle of intervals at . The approximation for the value at the origin converges to zero provided that , which is best visualized in a convergence study using a lin-lin plot (figure 1.4.1#fig.1, left).
The mid-point rule (sect.3.3) can be used to approximate an integral by computing the area of small boxes in order to evaluate the global quantity with . Global convergence is achieved for ; because of the weak singularity, the convergence rate drops from the expected from a quadrature of smooth functions to when .
|Copyright © Lifelong-learners at 09:16:54, January 18th, 2018|