Shooting method: Difference between revisions
No edit summary |
|||
Line 25: | Line 25: | ||
:<math> f(t, y(t), y'(t))=p(t)y'(t)+q(t)y(t)+r(t). \, </math> |
:<math> f(t, y(t), y'(t))=p(t)y'(t)+q(t)y(t)+r(t). \, </math> |
||
In this case, the solution to the boundary value problem is usually given by: |
In this case, the solution to the boundary value problem is usually given by: |
||
:<math>y(t) = y_{(1)}(t)+\frac{ |
:<math>y(t) = y_{(1)}(t)+\frac{y_{(1)}(t)-y_{(1)}(t_1)}{y_{(2)}(t_1)}y_{(2)}(t)</math> |
||
where <math>y_{(1)}(t)</math> is the solution to the initial value problem: |
where <math>y_{(1)}(t)</math> is the solution to the initial value problem: |
||
:<math>y_{(1)}''(t) = p(t)y_{(1)}'(t)+q(t)y_{(1)}(t)+r(t),\quad y_{(1)}(t_0) = y_0, \quad y_{(1)}'(t_0) = 0, </math> |
:<math>y_{(1)}''(t) = p(t)y_{(1)}'(t)+q(t)y_{(1)}(t)+r(t),\quad y_{(1)}(t_0) = y_0, \quad y_{(1)}'(t_0) = 0, </math> |
Revision as of 07:36, 29 December 2014
In numerical analysis, the shooting method is a method for solving a boundary value problem by reducing it to the solution of an initial value problem. The following exposition may be clarified by this illustration of the shooting method.
For a boundary value problem of a second-order ordinary differential equation, the method is stated as follows. Let
be the boundary value problem. Let y(t; a) denote the solution of the initial value problem
Define the function F(a) as the difference between y(t1; a) and the specified boundary value y1.
If F has a root a then obviously the solution y(t; a) of the corresponding initial value problem is also a solution of the boundary value problem. On the other hand, if the boundary value problem has a solution y(t), then y(t) is also the unique solution y(t; a) of the initial value problem where a = y'(t0), thus a is a root of F.
The usual methods for finding roots may be employed here, such as the bisection method or Newton's method.
Linear shooting method
The boundary value problem is linear if f has the form
In this case, the solution to the boundary value problem is usually given by:
where is the solution to the initial value problem:
and is the solution to the initial value problem:
See the proof for the precise condition under which this result holds.
Example
A boundary value problem is given as follows by Stoer and Burlisch (Section 7.3.1).
was solved for s = −1, −2, −3, ..., −100, and F(s) = w(1;s) − 1 plotted in the first figure. Inspecting the plot of F, we see that there are roots near −8 and −36. Some trajectories of w(t;s) are shown in the second figure.
Solutions of the initial value problem were computed by using the LSODE algorithm, as implemented in the mathematics package GNU Octave.
Stoer and Bulirsch state that there are two solutions, which can be found by algebraic methods. These correspond to the initial conditions w′(0) = −8 and w′(0) = −35.9 (approximately).
See also
References
- Josef Stoer and Roland Bulirsch. Introduction to Numerical Analysis. New York: Springer-Verlag, 1980. (See Section 7.3.)
- Press, WH; Teukolsky, SA; Vetterling, WT; Flannery, BP (2007). "Section 18.1. The Shooting Method". Numerical Recipes: The Art of Scientific Computing (3rd ed.). New York: Cambridge University Press. ISBN 978-0-521-88068-8.
External links
- Brief Description of ODEPACK (at Netlib; contains LSODE)
- Shooting method of solving boundary value problems – Notes, PPT, Maple, Mathcad, Matlab, Mathematica at Holistic Numerical Methods Institute [1]
- Shooting Method for Boundary Value Problems