Rate of convergence

From Wikipedia, the free encyclopedia

In numerical analysis (a branch of mathematics), the speed at which a convergent sequence approaches its limit is called the rate of convergence. Although strictly speaking, a limit does not give information about any finite first part of the sequence, this concept is of practical importance if we deal with a sequence of successive approximations for an iterative method, as then typically fewer iterations are needed to yield a useful approximation if the rate of convergence is higher. This may even make the difference between needing ten or a million iterations.

Series acceleration is a collection of techniques for improving the rate of convergence of a series. Such acceleration is commonly accomplished with sequence transformations.

Contents

[edit] Basic definition

Suppose that the sequence {xk} converges to the number ξ.

We say that this sequence converges linearly to ξ, if there exists a number μ ∈ (0, 1) such that

\lim_k \frac{|x_{k+1}-\xi|}{|x_k-\xi|} = \mu.

The number μ is called the rate of convergence.

If the above holds with μ = 0, then the sequence is said to converge superlinearly. One says that the sequence converges sublinearly if it converges, but (1) does not hold for any μ < 1.

The next definition is used to distinguish superlinear rates of convergence. We say that the sequence converges with order q for q > 1 to ξ if

\lim_k \frac{|x_{k+1}-\xi|}{|x_k-\xi|^q} = \mu \mbox{ with } \mu > 0

In particular, convergence with order 2 is called quadratic convergence, and convergence with order 3 is called cubic convergence.

[edit] Extended definition

The drawback of the above definitions (1) and (2) is that these do not catch some sequences which still converge reasonably fast, but whose "speed" is variable, such as the sequence {bk} below. Therefore, the definition of rate of convergence is sometimes extended as follows.

Under the new definition, the sequence {xk} converges with at least order q if there exists a sequence {εk} such that

|x_k - \xi| \le \varepsilon_k \quad\mbox{for all } k,

and the sequence {εk} converges to zero with order q according to the above "simple" definition.

[edit] Acceleration of convergence

Many methods exist to increase the rate of convergence of a given sequence, i.e. to transform a given sequence into one converging faster to the same limit. Such techniques are in general known as "series acceleration". The goal of the transformed sequence is to be much less "expensive" to calculate than the original sequence. One example of series acceleration is Aitken's delta-squared process.

[edit] Examples

Consider the following sequences:

a_0 = 1 ,\, a_1 = \frac12 ,\, a_2 = \frac14 ,\, a_3 = \frac18 ,\, a_4 = \frac1{16} ,\, a_5 = \frac1{32} ,\, \ldots ,\, a_k = \frac1{2^k} ,\, \ldots
b_0 = 1 ,\, b_1 = 1 ,\, b_2 = \frac14 ,\, b_3 = \frac14 ,\, b_4 = \frac1{16} ,\, b_5 = \frac1{16} ,\, \ldots ,\, b_k = \frac1{4^{\operatorname{floor}(k/2)}} ,\, \ldots
c_0 = \frac12 ,\, c_1 = \frac14 ,\, c_2 = \frac1{16} ,\, c_3 = \frac1{256} ,\, c_4 = \frac1{65536} ,\, \ldots ,\, c_k = \frac1{2^{2^k}} ,\, \ldots
d_0 = 1 ,\, d_1 = \frac12 ,\, d_2 = \frac13 ,\, d_3 = \frac14 ,\, d_4 = \frac15 ,\, d_5 = \frac16 ,\, \ldots ,\, d_k = \frac1{k+1} ,\, \ldots

The sequence {ak} converges linearly to 0 with rate 1/2. More generally, the sequence k converges linearly with rate μ if |μ| < 1. The sequence {bk} also converges linearly to 0 with rate 1/2 under the extended definition, but not under the simple definition. The sequence {ck} converges superlinearly. In fact, it is quadratically convergent. Finally, the sequence {dk} converges sublinearly.

[edit] References

The simple definition is used in

  • Michelle Schatzman (2002), Numerical analysis: a mathematical introduction, Clarendon Press, Oxford. ISBN 0-19-850279-6.

The extended definition is used in

  • Kendell A. Atkinson (1988), An introduction to numerical analysis (2nd ed.), John Wiley and Sons. ISBN 0-471-50023-2.
  • Walter Gautschi (1997), Numerical analysis: an introduction, Birkhäuser, Boston.
  • Endre Süli and David Mayers (2003), An introduction to numerical analysis, Cambridge University Press. ISBN 0-521-00794-1.