Poincaré–Lindstedt method

In perturbation theory, the Poincaré–Lindstedt method or Lindstedt–Poincaré method is a technique for uniformly approximating periodic solutions to ordinary differential equations, when regular perturbation approaches fail. The method removes secular terms—terms growing without bound—arising in the straightforward application of perturbation theory to weakly nonlinear problems with finite oscillatory solutions.

The method is named after Henri Poincaré, and Anders Lindstedt.

Example: the Duffing equation

The undamped, unforced Duffing equation is given by

:$\ddot + x + \varepsilon\, x^3 = 0\,$

for ''t'' > 0, with 0 < ''ε'' ≪ 1.J. David Logan. ''Applied Mathematics'', Second Edition, John Wiley & Sons, 1997. .

Consider initial conditions

:$x(0) = 1,\,$ $\dot x(0) = 0.\,$

A perturbation-series solution of the form ''x''(''t'') = ''x''₀(''t'') + ''ε'' ''x''₁(''t'') + …  is sought. The first two terms of the series are

:x(t) = \cos(t) + \varepsilon \left \tfrac\, \left( \cos(3t) - \cos(t) \right) - \tfrac\, t\, \sin(t) \right+ \cdots.\,

This approximation grows without bound in time, which is inconsistent with the physical system that the equation models.The Duffing equation has an invariant energy $\scriptstyle E=\tfrac12\, \dot^2 + \tfrac12\, x^2 + \tfrac14\, \varepsilon\, x^4$ = constant, as can be seen by multiplying the Duffing equation with $\scriptstyle \dot$ and integrating with respect to time ''t''. For the example considered, from its initial conditions, is found: ''E'' = ½ + ¼ ''ε''. The term responsible for this unbounded growth, called the secular term, is $t\sin(t)$. The Poincaré–Lindstedt method allows for the creation of an approximation that is accurate for all time, as follows.

In addition to expressing the solution itself as an asymptotic series, form another series with which to scale time ''t'':

:$\tau = \omega t,\,$ where $\omega = \omega_0 + \varepsilon \omega_1 + \cdots.\,$

For convenience, take ''ω''₀ = 1 because the leading order of the solution's angular frequency is 1. Then the original problem becomes

:$\omega^2\, x''(\tau) + x(\tau) + \varepsilon\, x^3(\tau) = 0\,$

with the same initial conditions. Now search for a solution of the form ''x''(''τ'') = ''x''₀(''τ'') + ''ε'' ''x''₁(''τ'') + … . The following solutions for the zeroth and first order problem in ''ε'' are obtained:

:$\begin x_0 &= \cos(\tau) \\ \text x_1 &= \tfrac\, \left(\cos(3\tau)-\cos(\tau)\right) + \left( \omega_1 - \tfrac \right)\, \tau\, \sin(\tau). \end$

So the secular term can be removed through the choice: ''ω''₁ = . Higher orders of accuracy can be obtained by continuing the perturbation analysis along this way. As of now, the approximation—correct up to first order in ''ε''—is

:
x(t) \approx \cos\Bigl(\left(1 + \tfrac\, \varepsilon \right)\, t \Bigr)
+ \tfrac\, \varepsilon\, \left cos\Bigl( 3 \left(1 + \tfrac\,\varepsilon\, \right)\, t \Bigr)-\cos\Bigl(\left(1 + \tfrac\,\varepsilon\, \right)\, t \Bigr)\right \,

References and notes

{{DEFAULTSORT:Poincare-Lindstedt method
Perturbation theory