Collocation method

In mathematics, a collocation method is a method for the numerical solution of ordinary differential equation and partial differential equations and integral equations. The idea is to choose a finite-dimensional space of candidate solutions (usually, polynomials up to a certain degree) and a number of points in the domain (called collocation points), and to select that solution which satisfies the given equation at the collocation points.

Ordinary differential equations

Suppose that the ordinary differential equation

$y'(t)=f(t,y(t)),\quad y(t_{0})=y_{0},$

is to be solved over the interval [t₀, t₀+h]. Denote the collocation points by 0 ≤ c₁< c₂< … < c_n ≤ 1.

The corresponding (polynomial) collocation method approximates the solution y by the polynomial p of degree n which satisfies the initial condition p(t₀) = y₀, and the differential equation p'(t) = f(t,p(t)) at all points t = t₀ + c_kh where k = 1, …, n. This gives n + 1 conditions, which matches the n + 1 parameters needed to specify a polynomial of degree n.

All these collocation methods are in fact implicit Runge–Kutta methods. However, not all Runge–Kutta methods are collocation methods.

Example

Pick, as an example, the two collocation points $c_{1}=0$ and $c_{2}=1$ (so $n=2$). The collocation conditions are

$p(t_{0})=y_{0},\,$

$p'(t_{0})=f(t_{0},p(t_{0})),\,$

$p'(t_{0}+h)=f(t_{0}+h,p(t_{0}+h)).\,$

There are three conditions, so p should be a polynomial of degree 2. Write p in the form

$p(t)=\alpha (t-t_{0})^{2}+\beta (t-t_{0})+\gamma \,$

to simplify the computations. Then the collocation conditions can be solved to give the coefficients

${\begin{aligned}\alpha &={\frac {1}{2h}}{\Big (}f(t_{0}+h,p(t_{0}+h))-f(t_{0},p(t_{0})){\Big )},\\\beta &=f(t_{0},p(t_{0})),\\\gamma &=y_{0}.\end{aligned}}$

The collocation method is now given (implicitly) by

$y_{1}=p(t_{0}+h)=y_{0}+{\frac {1}{2}}h{\Big (}f(t_{0}+h,y_{1})+f(t_{0},y_{0}){\Big )},\,$

where $y_{1}=p(t_{0}+h)$ is the approximate solution at $t=t_{0}+h$.

This method is known as the "trapezoidal rule." Indeed, this method can also be derived by rewriting the differential equation as

$y(t)=y(t_{0})+\int _{t_{0}}^{t}f(\tau ,y(\tau ))\,{\textrm {d}}\tau ,\,$

and approximating the integral on the right-hand side by the trapezoidal rule for integrals.

References

• Ernst Hairer, Syvert Nørsett and Gerhard Wanner, Solving ordinary differential equations I: Nonstiff problems, second edition, Springer Verlag, Berlin, 1993. ISBN 3-540-56670-8.
• Arieh Iserles, A First Course in the Numerical Analysis of Differential Equations, Cambridge University Press, 1996. ISBN 0-521-55376-8 (hardback), ISBN 0-521-55655-4 (paperback).