Direct multiple shooting method

Direct multiple shooting method

In the area of mathematics known as numerical ordinary differential equations, the direct multiple shooting method is a numerical method for the solution of boundary value problems. The method divides the interval over which a solution is sought into several smaller intervals, solves an initial value problem in each of the smaller intervals, and imposes additional matching conditions to form a solution on the whole interval. The method constitutes a significant improvement in distribution of nonlinearity and numerical stability over single shooting methods.

Contents

Single shooting methods

Shooting methods can be used to solve boundary value problems (BVP) like

 y'(t) = f(t, y(t)), \quad y(t_a) = y_a, \quad y(t_b) = y_b,

in which the time points ta and tb are known but the initial and terminal values ya and yb are unknown and sought for.

Single shooting methods proceed as follows. Let y(t; t0, y0) denote the solution of the initial value problem (IVP)

 y'(t) = f(t, y(t)), \quad y(t_0) = y_0.

Define the function F(p) as the difference between y(tb; p) and the specified boundary value yb: F(p) = y(tb; p) − yb. Then for every solution (ya, yb) of the boundary value problem we have ya=y0 while yb corresponds to a root of F. This root be solved by any root-finding method given that certain method-dependent prerequisites are satisfied. This often will require initial guesses to ya and yb. Typically, analytic root finding is impossible and iterative methods such as Newton's method are used for this task.

The application of single shooting for the numerical solution of boundary value problems suffers from several drawbacks.

  • For a given initial value y0 the solution of the IVP obviously must exist on the interval [ta,tb] so that we can evaluate the function F whose root is sought.

For highly nonlinear or unstable ODEs, this requires the initial guess y0 to be extremely close to an actual but unknown solution ya. Initial values that are chosen slightly off the true solution may lead to singularities or breakdown of the ODE solver method. Choosing such solutions is inevitable in an iterative root-finding method, however.

  • Finite precision numerics may make it impossible at all to find initial values that allow for the solution of the ODE on the whole time interval.
  • The nonlinearity of the ODE effectively becomes a nonlinearity of F, and requires a root-finding technique capable of solving nonlinear systems. Such methods typically converge slower as nonlinearities become more severe. The boundary value problem solver's performance suffers from this.
  • Even stable and well-conditioned ODEs may make for unstable and ill-conditioned BVPs. A slight alteration of the initial value guess y0 may generate an extremely large step in the ODEs solution y(tb; ta, y0) and thus in the values of the function F whose root is sought. Non-analytic root-finding methods can seldom cope with this behaviour.

Multiple shooting

A direct multiple shooting method partitions the interval [ta, tb] by introducing additional grid points

 t_a = t_0 < t_1 < \cdots < t_N = t_b .

The method starts by guessing somehow the values of y at all grid points tk with 0 ≤ kN − 1. Denote these guesses by yk. Let y(t; tk, yk) denote the solution emanating from the kth grid point, that is, the solution of the initial value problem

 y'(t) = f(t, y(t)), \quad y(t_k) = y_k.

All these solutions can be pieced together to form a continuous trajectory if the values y match at the grid points. Thus, solutions of the boundary value problem correspond to solutions of the following system of N equations:

 \begin{align}
& y(t_1; t_0, y_0) = y_1 \\
& \qquad\qquad\vdots \\
& y(t_{N-1}; t_{N-2}, y_{N-2}) = y_{N-1} \\
& y(t_N; t_{N-1}, y_{N-1}) = y_b.
\end{align}

The central N−2 equations are the matching conditions, and the first and last equations are the conditions y(ta) = ya and y(tb) = yb from the boundary value problem. The multiple shooting method solves the boundary value problem by solving this system of equations. Typically, a modification of the Newton's method is used for the latter task.

Bock's Direct Multiple Shooting Method

References


Wikimedia Foundation. 2010.

Нужно сделать НИР?

Look at other dictionaries:

  • Shooting method — 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… …   Wikipedia

  • Direct-Shift Gearbox — Transmission types Manual Sequential manual Non synchronous Preselector Automatic Manumatic Semi automatic Electrohydraulic Dual …   Wikipedia

  • List of numerical analysis topics — This is a list of numerical analysis topics, by Wikipedia page. Contents 1 General 2 Error 3 Elementary and special functions 4 Numerical linear algebra …   Wikipedia

  • List of mathematics articles (D) — NOTOC D D distribution D module D D Agostino s K squared test D Alembert Euler condition D Alembert operator D Alembert s formula D Alembert s paradox D Alembert s principle Dagger category Dagger compact category Dagger symmetric monoidal… …   Wikipedia

  • Optimal control — theory, an extension of the calculus of variations, is a mathematical optimization method for deriving control policies. The method is largely due to the work of Lev Pontryagin and his collaborators in the Soviet Union[1] and Richard Bellman in… …   Wikipedia

  • Model predictive control — Model Predictive Control, or MPC, is an advanced method of process control that has been in use in the process industries such as chemical plants and oil refineries since the 1980s. Model predictive controllers rely on dynamic models of the… …   Wikipedia

  • Law, Crime, and Law Enforcement — ▪ 2006 Introduction Trials of former heads of state, U.S. Supreme Court rulings on eminent domain and the death penalty, and high profile cases against former executives of large corporations were leading legal and criminal issues in 2005.… …   Universalium

  • Foley, Jack — (1891–1967)    While his first name may not elicit recognition in movie fans, undoubtedly most have heard of the technique which bears his surname: Foley, or the act of performing sound effects during postproduction to match the action of the… …   The Encyclopedia of Stanley Kubrick

  • AR-15 — The AR 15 comes in many sizes and has many options, depending on the manufacturer. The part shown bottom center is the lower receiver, which under US law is the component legally considered the firearm . Type …   Wikipedia

  • Cinema of Hungary — List of Hungarian films 1901 1947 1948 1989 1990 present Actors • Directors • Producers European cinema …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”