Numerical Methods in Computational Finance. Daniel J. Duffy

Чтение книги онлайн.

Читать онлайн книгу Numerical Methods in Computational Finance - Daniel J. Duffy страница 39

Numerical Methods in Computational Finance - Daniel J. Duffy

Скачать книгу

behaves for large values of time; the answer depends on a relationship between the drift and diffusion parameters:

StartLayout 1st Row 1st Column normal i right-parenthesis 2nd Column mu greater-than one half sigma squared comma limit Underscript t right-arrow infinity Endscripts xi left-parenthesis t right-parenthesis equals infinity almost surely 2nd Row 1st Column ii right-parenthesis 2nd Column mu less-than one half sigma squared comma limit Underscript t right-arrow infinity Endscripts xi left-parenthesis t right-parenthesis equals 0 almost surely 3rd Row 1st Column iii right-parenthesis 2nd Column mu equals one half sigma squared comma xi left-parenthesis t right-parenthesis fluctuates between arbitrary large and 4th Row 1st Column Blank 2nd Column arbitrary small values as t right-arrow infinity almost surely period EndLayout

      In general it is not possible to find an exact solution, and in these cases we must resort to numerical approximation techniques.

      Equation (3.17) is a one-factor equation because there is only one dependent variable (namely xi left-parenthesis t right-parenthesis) to be modelled. It is possible to define equations with several dependent variables. The prototypical non-linear stochastic differential equation is given by the system:

      where:

StartLayout 1st Row 1st Column mu 2nd Column colon left-bracket 0 comma upper T right-bracket times normal double struck upper R Superscript n Baseline right-arrow normal double struck upper R Superscript n Baseline left-parenthesis vector right-parenthesis 2nd Row 1st Column sigma 2nd Column colon left-bracket 0 comma upper T right-bracket times normal double struck upper R Superscript n Baseline right-arrow normal double struck upper R Superscript n times m Baseline left-parenthesis matrix right-parenthesis 3rd Row 1st Column upper X 2nd Column colon left-bracket 0 comma upper T right-bracket right-arrow normal double struck upper R Superscript n Baseline left-parenthesis vector right-parenthesis 4th Row 1st Column upper W 2nd Column colon left-bracket 0 comma upper T right-bracket right-arrow normal double struck upper R Superscript m Baseline left-parenthesis vector right-parenthesis period EndLayout StartLayout 1st Row mu equals mu left-parenthesis t comma upper X right-parenthesis 2nd Row sigma equals sigma left-parenthesis t comma upper X right-parenthesis period EndLayout

      This is a generalisation of Equation (3.17). Thus, instead of scalars this system employs vectors for the solution, drift and random number terms while the diffusion term is a rectangular matrix.

      Existence and uniqueness theorems for the solution of the SDE system (3.20) are similar to those in the one-factor case. For example, theorem 5.2.1 in Øksendal (1998) addresses these issues. We discuss SDEs in more detail in Chapter 13.

      In this section we introduce a class of one-step methods to approximate the solution of ODE system (3.1).

      The first step is to replace continuous time by discrete time. To this end, we divide the interval [0, T] into a number of subintervals. We define upper N plus 1 mesh points as follows:

0 equals t 0 less-than t 1 less-than ellipsis less-than t Subscript n Baseline less-than t Subscript n plus 1 Baseline less-than ellipsis less-than t Subscript upper N Baseline equals upper T period

      In this case we define a set of subintervals left-parenthesis t Subscript n Baseline comma t Subscript n plus 1 Baseline right-parenthesis of size normal upper Delta t Subscript n Baseline identical-to t Subscript n plus 1 Baseline minus t Subscript n, 0 less-than-or-equal-to n less-than-or-equal-to upper N minus 1.

      In general, we speak of a non-uniform mesh when the sizes of the subintervals are not necessarily the same. However, in this book we consider in the main a class of finite difference schemes where the N subintervals have the same length (we then speak of a uniform mesh), namely normal upper Delta t equals upper T slash upper N. The variable h equals upper T slash upper N is also used to denote the uniform mesh size.

      In general, we define y Subscript n to be the approximate solution at time t Subscript n and we write the functional dependence of y Subscript n plus 1 on t Subscript n Baseline comma y Subscript n Baseline and h by:

      where normal upper Phi is called the increment function. For example, in the case of the explicit Euler method, this function is:

normal upper Phi left-parenthesis t comma y comma h right-parenthesis equals f left-parenthesis t comma y right-parenthesis period

Скачать книгу