Basic numerical procedures
Download
1 / 38

Basic Numerical Procedures - PowerPoint PPT Presentation


  • 153 Views
  • Updated On :

Basic Numerical Procedures. Chapter 19. Tree Approaches to Derivatives Valuation. Trees Monte Carlo simulation Finite difference methods. Binomial Trees. Binomial trees are frequently used to approximate the movements in the price of a stock or other asset

loader
I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
capcha
Download Presentation

PowerPoint Slideshow about 'Basic Numerical Procedures' - nuru


An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript
Basic numerical procedures l.jpg

Basic Numerical Procedures

Chapter 19

Zheng Zhenlong, Dept of Finance,XMU


Tree approaches to derivatives valuation l.jpg
Tree Approaches to Derivatives Valuation

  • Trees

  • Monte Carlo simulation

  • Finite difference methods

Zheng Zhenlong, Dept of Finance,XMU


Binomial trees l.jpg
Binomial Trees

  • Binomial trees are frequently used to approximate the movements in the price of a stock or other asset

  • In each small interval of time the stock price is assumed to move up by a proportional amount u or to move down by a proportional amount d

Zheng Zhenlong, Dept of Finance,XMU


Movements in time d t figure 19 1 page 400 l.jpg

Su

p

S

1 – p

Sd

Movements in Time Dt(Figure 19.1, page 400)

Zheng Zhenlong, Dept of Finance,XMU


1 tree parameters for asset paying a dividend yield of q l.jpg
1. Tree Parameters for asset paying a dividend yield of q

Parameters p, u, and d are chosen so that the tree gives correct values for the mean & variance of the stock price changes in a risk-neutral world

Mean: e(r-q)Dt= pu + (1– p )d

Variance: s2Dt = pu2 + (1– p )d2 – e2(r-q)Dt

A further condition often imposed is u = 1/ d

Zheng Zhenlong, Dept of Finance,XMU


2 tree parameters for asset paying a dividend yield of q equations 19 4 to 19 7 l.jpg
2. Tree Parameters for asset paying a dividend yield of q(Equations 19.4 to 19.7)

When Dt is small a solution to the equations is

Zheng Zhenlong, Dept of Finance,XMU


The complete tree figure 19 2 page 402 l.jpg

S0u3

S0u2

S0u

S0u

S0

S0

S0

S0d

S0d

S0d2

S0d3

The Complete Tree(Figure 19.2, page 402)

S0u4

S0u2

S0d 2

S0d4

Zheng Zhenlong, Dept of Finance,XMU


Backwards induction l.jpg
Backwards Induction

  • We know the value of the option at the final nodes

  • We work back through the tree using risk-neutral valuation to calculate the value of the option at each node, testing for early exercise when appropriate

Zheng Zhenlong, Dept of Finance,XMU


Example put option example 19 1 page 402 l.jpg
Example: Put Option(Example 19.1, page 402)

S0= 50; K = 50; r =10%; s = 40%;

T = 5 months = 0.4167;

Dt = 1 month = 0.0833

The parameters imply

u = 1.1224; d = 0.8909;

a = 1.0084; p = 0.5073

Zheng Zhenlong, Dept of Finance,XMU


Example continued figure 19 3 page 403 l.jpg
Example (continued)Figure 19.3, page 403

Zheng Zhenlong, Dept of Finance,XMU


Calculation of delta l.jpg
Calculation of Delta

Delta is calculated from the nodes at time Dt

Zheng Zhenlong, Dept of Finance,XMU


Calculation of gamma l.jpg
Calculation of Gamma

Gamma is calculated from the nodes at time 2Dt

Zheng Zhenlong, Dept of Finance,XMU


Calculation of theta l.jpg
Calculation of Theta

Theta is calculated from the central nodes at times 0 and 2Dt

Zheng Zhenlong, Dept of Finance,XMU


Calculation of vega l.jpg
Calculation of Vega

  • We can proceed as follows

  • Construct a new tree with a volatility of 41% instead of 40%.

  • Value of option is 4.62

  • Vega is

Zheng Zhenlong, Dept of Finance,XMU


Trees for options on indices currencies and futures contracts l.jpg
Trees for Options on Indices, Currencies and Futures Contracts

As with Black-Scholes:

  • For options on stock indices,qequals the dividend yield on the index

  • For options on a foreign currency, qequals the foreign risk-free rate

  • For options on futures contracts q = r

Zheng Zhenlong, Dept of Finance,XMU


Binomial tree for dividend paying stock l.jpg
Binomial Tree for Dividend Paying Stock Contracts

  • Procedure:

    • Draw the tree for the stock price less the present value of the dividends

    • Create a new tree by adding the present value of the dividends at each node

  • This ensures that the tree recombines and makes assumptions similar to those when the Black-Scholes model is used

Zheng Zhenlong, Dept of Finance,XMU


Extensions of tree approach l.jpg
Extensions of Tree Approach Contracts

  • Time dependent interest rates(用远期利率)

  • The control variate technique

    fA+fbs-fE

Zheng Zhenlong, Dept of Finance,XMU


Alternative binomial tree section 19 4 page 414 l.jpg
Alternative Binomial Tree Contracts(Section 19.4, page 414)

Instead of setting u = 1/d we can set each of the 2 probabilities to 0.5 and

Zheng Zhenlong, Dept of Finance,XMU


Trinomial tree page 409 l.jpg

Su Contracts

pu

pm

S

S

pd

Sd

Trinomial Tree (Page 409)

Zheng Zhenlong, Dept of Finance,XMU


Time dependent parameters in a binomial tree page 409 l.jpg
Time Dependent Parameters in a Binomial Tree Contracts(page 409)

  • Making r or q a function of time does not affect the geometry of the tree. The probabilities on the tree become functions of time.

  • We can make s a function of time by making the lengths of the time steps inversely proportional to the variance rate.

Zheng Zhenlong, Dept of Finance,XMU


Monte carlo simulation l.jpg
Monte Carlo Simulation Contracts

When used to value European stock options, Monte Carlo simulation involves the following steps:

1. Simulate 1 path for the stock price in a risk neutral world

2. Calculate the payoff from the stock option

3. Repeat steps 1 and 2 many times to get many sample payoff

4. Calculate mean payoff

5. Discount mean payoff at risk free rate to get an estimate of the value of the option

Zheng Zhenlong, Dept of Finance,XMU


Sampling stock price movements equations 19 13 and 19 14 page 419 l.jpg
Sampling Stock Price Movements Contracts(Equations 19.13 and 19.14, page 419)

  • In a risk neutral world the process for a stock price is

  • We can simulate a path by choosing time steps of length Dt and using the discrete version of this

    where e is a random sample from f(0,1)

Zheng Zhenlong, Dept of Finance,XMU


A more accurate approach equation 19 15 page 420 l.jpg
A More Accurate Approach Contracts(Equation 19.15, page 420)

Zheng Zhenlong, Dept of Finance,XMU


Extensions l.jpg
Extensions Contracts

When a derivative depends on several underlying variables we can simulate paths for each of them in a risk-neutral world to calculate the values for the derivative

Zheng Zhenlong, Dept of Finance,XMU


Sampling from normal distribution page 422 l.jpg
Sampling from Normal Distribution Contracts(Page 422)

  • One simple way to obtain a sample

    from f(0,1) is to generate 12 random numbers between 0.0 & 1.0, take the sum, and subtract 6.0

  • In Excel =NORMSINV(RAND()) gives a random sample from f(0,1)

Zheng Zhenlong, Dept of Finance,XMU


To obtain 2 correlated normal samples l.jpg
To Obtain 2 Correlated Normal Samples Contracts

  • Obtain independent normal samples x1 and x2 and set

  • A procedure known as Cholesky’s decomposition when samples are required from more than two normal variables (see page 422)

Zheng Zhenlong, Dept of Finance,XMU


Standard errors in monte carlo simulation l.jpg
Standard Errors in Monte Carlo Simulation Contracts

The standard error of the estimate of the option price is the standard deviation of the discounted payoffs given by the simulation trials divided by the square root of the number of observations.

Zheng Zhenlong, Dept of Finance,XMU


Application of monte carlo simulation l.jpg
Application of Monte Carlo Simulation Contracts

  • Monte Carlo simulation can deal with path dependent options, options dependent on several underlying state variables, and options with complex payoffs

  • It cannot easily deal with American-style options

Zheng Zhenlong, Dept of Finance,XMU


Determining greek letters l.jpg
Determining Greek Letters Contracts

For D:

1. Make a small change to asset price

2. Carry out the simulation again using the same random number streams

3. Estimate D as the change in the option price divided by the change in the asset price

Proceed in a similar manner for other Greek letters

Zheng Zhenlong, Dept of Finance,XMU


Variance reduction techniques l.jpg
Variance Reduction Techniques Contracts

  • Antithetic variable technique(对偶变量技术)

  • Control variate technique

  • Importance sampling(虚值部分不用抽)

  • Stratified sampling(分层抽样,抽样次数已知)

  • Moment matching

  • Using quasi-random sequences(平衡抽样)

Zheng Zhenlong, Dept of Finance,XMU


Sampling through the tree l.jpg
Sampling Through the Tree Contracts

Instead of sampling from the stochastic process we can sample paths randomly through a binomial or trinomial tree to value a derivative

Zheng Zhenlong, Dept of Finance,XMU


Finite difference methods l.jpg
Finite Difference Methods Contracts

  • Finite difference methods aim to represent the differential equation in the form of a difference equation

  • We form a grid by considering equally spaced time values and stock price values

  • Define ƒi,jas the value of ƒ at time iDt when the stock price is jDS

Zheng Zhenlong, Dept of Finance,XMU


Finite difference methods continued l.jpg
Finite Difference Methods Contracts(continued)

Zheng Zhenlong, Dept of Finance,XMU


Implicit finite difference method equation 19 25 page 429 l.jpg
Implicit Finite Difference Method Contracts(Equation 19.25, page 429)

Zheng Zhenlong, Dept of Finance,XMU


Explicit finite difference method page 422 428 l.jpg
Explicit Finite Difference Method Contracts(page 422-428)

Zheng Zhenlong, Dept of Finance,XMU


Implicit vs explicit finite difference method l.jpg
Implicit vs Explicit Finite Difference Method Contracts

  • The explicit finite difference method is equivalent to the trinomial tree approach

  • The implicit finite difference method is equivalent to a multinomial tree approach

Zheng Zhenlong, Dept of Finance,XMU


Implicit vs explicit finite difference methods figure 19 16 page 433 l.jpg

ƒ Contractsi , j

ƒi +1, j

ƒi +1, j –1

Implicit vs Explicit Finite Difference Methods(Figure 19.16, page 433)

ƒi +1, j +1

ƒi , j +1

ƒi , j

ƒi +1, j

ƒi , j –1

Implicit Method

Explicit Method

Zheng Zhenlong, Dept of Finance,XMU


Other points on finite difference methods l.jpg
Other Points on Finite Difference Methods Contracts

  • It is better to have ln Srather thanSas the underlying variable

  • Improvements over the basic implicit and explicit methods:

    • Hopscotch method

    • Crank-Nicolson method(显性与隐性方法的平均)

Zheng Zhenlong, Dept of Finance,XMU