Path-dependent options and randomization techniques
|
|
- Leslie York
- 5 years ago
- Views:
Transcription
1 Path-dependent options and randomization techniques Daniele Marazzina Joint work with G. Fusai and M. Marena Dipartimento di Matematica F. Brioschi Politecnico di Milano April 8, 2009
2 The Talk In this talk, we will show how we can price discretely monitored options Arithmetic Asian Barrier European Lookback Benchmark combining maturity randomization and Fourier transform techniques Numerical approach Numerical solution of a set of independent integral equations in the Lévy world
3 Discrete knock-and-out barrier option: as you hit you die Payoff at maturity (T = 1) if still alive Path-Dependent
4 Methods for pricing path-dependent options Barrier and Lookback Options Airoldi, Quantitative Finance, 2005 Perturbative Series Expansion Approximation of a generic probability function around another simpler pdf by matching all moments up to a generic order Chiarella, El-Hassan & Kucera, Computational Methods in Financial Engineering, 2008 Path integral approach Backward recursion functional equation solved by expanding the pricing functions in Fourier-Hermite series
5 Methods for pricing path-dependent options Barrier and Lookback Options - 2 Feng & Linetsky, Mathematical Finance, 2008 Fast Hilbert Transform Approach Apply the Hilbert transform and a discretization procedure using the Whittaker cardinal series Backward induction in the Fourier space, rather than in the state space Jackson, Jaimungal & Surkov, Journal of Computational Finance, 2008 Fourier Space Time-stepping method Apply the Fourier transform to the pricing PIDE to obtain a linear system of easily solvable ordinary differential equations Lord, Fang, Bervoets & Oosterlee, SIAM Journal on Scientific Computing, 2008 Convolution method Based on Fourier transforms combined with quadrature rules
6 Methods for pricing path-dependent options Asian Options Benhamou, Journal of Computational Finance, 2002 Convolution method Albrecher & Predota, Journal of Computational and Applied Mathematics, 2004 Moments approximation Approximations of the arithmetic option price based on the moments of the average Fusai & Meucci, Journal of Banking and Finance, 2008 Recursive integration Cerny & Kyriakou, to appear Backward price convolution They provide analytical upper bound for the pricing error due to truncations
7 CONTRACTS
8 The setup 1. Risk-free interest rate r 2. Asset having log-price x, i.e., x t = log(s t) 3. Derivative contract having maturity T N monitoring dates = T the time-distance between monitoring dates N strike K
9 Barrier and European Options For European and Barrier options, the pricing problem can be written in a recursive way (see Hull, 2003) C (x, j) = e r f (x, ξ; ) C (ξ, j + 1) dξ j = N 1,, 0 Ω C (x, N) = φ (x) where Ω depends on the contract f (x, ξ; ) the transition density φ (x) is the payoff C(x, j) is the option price at the j th monitoring date
10 1: Choice of Ω European Options: options that can only be exercised at the end of their lives Ω = R Barrier Options: options die/start to live if a lower or upper barrier is hit Ω = (ln(l), + ) for lower barrier Ω = (, ln(u)) for upper barrier Ω = (ln(l), ln(u)) for double barriers
11 2: Choice of the transition density: Lévy World Assumption The underlying asset evolves according to a generic Lévy process Thus the increments are i.i.d. the transition density f (x, ξ; ) from ξ at time t to x at time t + is of convolution type, i.e. (with an abuse of notation) f (x, ξ; ) = f (x ξ; ) Why Lévy Process? The volatility surface shows the inconsistence of the Black-Scholes model In general the transition density f is known only through its Fourier transform it is a computational advantage that f is of convolution type
12 The considered Lévy processes are G (Gaussian) CGMY (Carr-Geman-Madan-Yor) DE (Double Exponential) Lévy World 0 JD (Jump Diffusion or Merton) log return G CGMY Figure: Density of the log-returns For the CGMY, the DE and the JD processes tails are fatter than in the Gaussian case: we better model the occurrence of large (negative) returns Model (Parameters) Characteristic Exponent ψ (γ) G(σ) 1 2 σ2 γ 2 ) CGMY (C, G, M, Y ) C Γ ( Y ) ((M iγ) Y M Y + (G + iγ) Y G Y ( ) DE(σ, λ, p, η 1, η 2) 1 2 σ2 γ 2 + λ (1 p) η2 + p η 1 1 η 2 iγ η 1 +iγ ( ) JD(σ, α, λ, δ) 1 2 σ2 γ 2 + λ e iγα 1 2 γ2 δ 2 1 DE JD
13 Lookback Options Their settlement is based on the minimum or the maximum value of the underlying asset as registered during the lifetime of the options Consider fixed strike lookback put options, i.e., the payoff is ( ) + (K M N ) + = e k e m N where M N := min j=0,,n S j m N := min j=0,,n x j x j is the asset log-price at the j th monitoring date, i.e., x j = log(s j ) k is the logarithm of the strike K
14 Assuming x 0 > k, the lookback put with fixed strike has price k e rt e u P(m N u) du Problem: P(m N u) =???
15 P(m N u) = 1 v(x 0 u, N) where v(, ) satisfies the following recursion v(z, j) = + 0 v(z, 0) = 1 (z>0) f (z ξ; )v(ξ, j 1)dξ j = 1,, N Solving the recursion, we obtain the distribution of the minimum Once this distribution is computed, we have to perform a numerical integration to obtain the option price
16 Asian Options The payoff of an arithmetic Asian option depends on the arithmetic average of the values of the underlying asset at the monitoring dates Payoff of a Call Option ( ) + Fixed-Strike Asian Options: AN K N+1 ( Floating-Strike Asian Options: S N A N N+1 where A N N j=0 S j ) +
17 Fixed-Strike Asian Options We define a random variable Z such that A j+1 = A j + S j e Z j+1, j = 0,, N 1 Letting L N e Z N and defining recursively the following quantities we obtain L j e Z j (1 + L j+1 ), j = N 1,, 1 A N S 0 (1 + L 1) The density of L 1 turns out to be the key variable in pricing fixed strike Asian options Option Price e rt + where f B1 is the density of B 1 log(l 1) ( ) + S0 T + 1 (1 + ex ) K f B1 (x) dx
18 Recursion Problem: How to compute f B1? with a recursion The density of f Bn satisfies the recursion f Bj (x) = + f (x log (e y + 1)) f Bj+1 (y) dy, j = N 1,, 1 with initial condition f BN f, i.e., the transition density
19 Floating-Strike Asian Options The following recursion for the option price holds for j = N 1,, 0 C (S N, A N ; N) = ( S N A ) + N N C (S j, A j ; j) = e r f (s) C (S j e s, A j + S j e s ; j + 1) ds
20 Since the payoff is homogeneous of degree one in the stock price, we can write ( C (S j, A j ; j) = S j C 1, A ) j ; j S j If we set then the following recursion holds + g j (x) = e r f with initial condition 1 g j (x) C (1, x; j) ( log g N (x) = ( x )) y 1 ( 1 x ) + N + 1 x (y 1) 2 g j+1 (y) dy The floating strike option price is S 0g 0 (1)
21 Maturity Randomization
22 All the recursion can be written as v (x, j) = K (x, ξ; ) v (ξ, j 1) dξ Ω v (x, 0) = ψ (x) j = 1,, M for suitable choices of K, ψ, M and Ω In practice, to solve the pricing problem, we have to compute M nested integrals
23 The strategy The problem: recursive computation of an integral 1. Idea: let us randomize the option expiry according to a geometric distribution 2. Randomization transforms the recursive pricing problem into a set of integral equations 3. Let us look for numerical solutions of the integral equations 4. De-randomize the expiry and get the option price
24 Randomize the option expiry Let us consider a coin tossing game where q is the probability of getting tails If we get heads for the first time after j tails, we receive an option having j monitoring dates Therefore the expected gain W (x, q) is W (x, q) = (1 q)q j v (x, j) j=0
25 Randomize the option expiry Let us consider a coin tossing game where q is the probability of getting tails If we get heads for the first time after j tails, we receive an option having j monitoring dates Therefore the expected gain W (x, q) is W (x, q) = (1 q)q j v (x, j) j=0 The recursive equation becomes an integral equation W (x, q) = q K (x, ξ; ) W (ξ, q) dξ + (1 q)ψ (x) Ω
26 De-Randomization How to go back from W to v Complex inversion integral v (x, M) = 1 2πρ M 2π 0 ( W x, ρe ) 1u e 1Mu du Numerical approximation { v h (x, M)= 1 M 1 W (x, ρ)+( 1) M W (x, ρ)+2 ( 1) j Re 2Mρ M Solution of M + 1 independent integral equations j=1 ( ( W x, ρe )) } 1jπ M
27 The Algorithm 1. Compute the transition density f using FFT 2. For j = 0,, M, set q (j) = ρ exp ( 1jπ/M ), ρ = 10 γ/2m, γ = 8 3. Discretization of the integral equation Truncate the domain Ω to [a, b], if necessary Choose a quadrature rule and define the m nodes x i We obtain the linear systems (I m q(j)k md m) W m = Ψ m j = 0,, M where I m is the identity matrix of size m K m is the m m square matrix with elements K ij = K(x i, x j ; ) D m is the m m matrix which contains the weights associated to the quadrature formula Ψ i = (1 q)ψ (x i ) 4. Solve the M + 1 linear systems 5. Using the M + 1 solutions of the integral equation parametrized by q (j), reconstruct v h (x, M)
28 Computational Cost
29 How to improve performances? 1.EULER ACCELERATION 2.PRECONDITIONING TECHNIQUES 3.GRID COMPUTING
30 1.Euler Acceleration
31 Euler Acceleration Euler Summation is a convergence-acceleration technique well suited for evaluating alternating series Numerical Inversion Formula { v h (x, M)= 1 M 1 W (x, ρ)+( 1) M W (x, ρ)+2 ( 1) j Re 2Mρ M This formula can be rewritten as v h (x, M) = 1 ρ M j=0 j=1 ( ( W x, ρe M ( ( 1) j a j Re W (x, ρe )) 1jπ/M )) } 1jπ M
32 The idea of the Euler acceleration technique is to approximate the above sum with the following one ( ) v h (x, M) 1 m e m b 2 me ρ M ne +j j where b k = j=0 j=0 k ( ( 1) j a j Re W (x, ρe )) 1jπ/M with n e and m e suitably chosen, such that n e + m e < M Thus the number of linear systems to be solved is min{m, n e + m e} + 1 We set n e = 12 and m e = 10
33 2.Preconditioning Technique Barrier and Lookback Options K(, ) is of convolution type we deal with Toeplitz type matrices [Joint work with M. Ng]
34 Barrier and Lookback Options We consider Composite Rectangle, Trapezoidal and Simpson s formula (I qkd) W = Φ All the linear systems are solved with the GMRES iterative method Composite Rectangle K is a Toeplitz matrix D is the identity matrix I qkd is a Toeplitz matrix Composite Trapezoidal and Simpson s they are higher-order formulas K is a Toeplitz matrix I qkd is not a Toeplitz matrix Toeplitz Matrix Each descending diagonal from left to right is constant Multiplication of a m m Toeplitz matrix to a vector requires O(m log m) operations m m Toeplitz matrix store only 2m 1 elements
35 Preconditioner A preconditioner P of a matrix A is a matrix such that P 1 A has a smaller condition number than A Preconditioners are useful when using an iterative method to solve a large linear system Ax = b, since the rate of convergence for most iterative linear solvers degrades as the condition number of a matrix increases: instead of solving the original linear system, one solve P 1 Ax = P 1 b Lin, Ng & Chan Preconditioner Based on Circulant Matrix Already Inverted Lin, Ng & Chan
36 GMRES iterations GMRES Iterations & Execution Time (ETime - in sec): m = 3001 No Preconditioner CR CT CS M Iterations ETime Iterations ETime Iterations ETime Preconditioner CR CT CS M Iterations ETime Iterations ETime Iterations ETime Computational Cost: O(min{M, n e + m e}m log m)
37 2.Preconditioning Technique Asian Options K(, ) is no more of convolution type Reichel algorithm
38 Asian Options We consider Gauss-Chebyshew quadrature rules for the following linear system (I qkd) W = Ψ In order to speed up the solution of the above linear system, we consider an algorithm due to Reichel, who proposes a fast solution method for the one-dimensional Fredholm integral equation Reichel proves that the matrix K can be well approximated by a matrix K of rank much smaller than m Thus he proposes an iterative algorithm with a suitable preconditioner (based on K ) having a computational cost of O(m 2 ) arithmetic operations
39 Fixed Strike Asian Option Price. M = 50 Reichel GMRes GE m Price CPU Price CPU Price CPU Benchmark price: Cérny and Kyriakou
40 How to improve performances? 3.Grid Computing
41 Grid Solving independent integral equations Grid Computing Solve linear systems on different computers Grid Computing is a special type of parallel computing which relies on complete computers connected to a network Our grid consists of six Intel Pentium systems, each equipped with 4GB of RAM and Intel Core 2 Quad Q6600 (2400MHz) processor Floating Strike Asian Option - CPU Time (in seconds) N m = 1000 m = 2000 m = 3000 m = 4000 Serial Grid
42 NUMERICAL RESULTS
43 Barrier Options Double Barrier Options and Double Exponential Distribution: N = 252 Logarithmic Scales Benchmark Price: Feng & Linetsky, Rectangle Trapezoidal Simpson 10 3 Pointwise Error m
44 Lookback Options Lookback Options and DE Model - Prices N m = 1001 m = 2001 m = 3001 m = 4001 CI Composite Rectangle Composite Trapezoidal Composite Simpson s Parameters: r = , T = 1, S 0 = 100, and K = 100 DE: σ = , λ = , p = , η 1 = and η 2 = m >> 4001
45 Fixed Strike Asian Options Fixed Strike Asian Options: a comparison with Cérny & Kyriakou Model Parameters m Price Benchmark G σ = CGMY C = G = M = Y = NIG δ = α = β = Parameters: r = 0.04, T = 1, S 0 = 100, M = 50 and K = 100
46 Floating Strike Asian Options Floating Strike Asian Options: a comparison with Lim (2002) Model m Price Benchmark G Parameters: σ = 0.2, r = 0.1, T = 182, S0 = 100, M = 91 and K =
47 Conclusions We have presented a method for pricing options intuitive easy to be implemented accurate can be used with different products and processes
48 Conclusions We have presented a method for pricing options intuitive easy to be implemented accurate can be used with different products and processes References G.Fusai, D.M., M.Marena, Option Pricing, Maturity Randomization and Grid Computing, IPDPS2008 Proceedings G.Fusai, D.M., M.Marena, M.Ng, Randomization and Preconditioning Techniques for Option Pricing, submitted G.Fusai, D.M., M.Marena, Pricing Discretely Monitored Asian Options Using Maturity Randomization, in preparation
49 Discrete and continuous monitoring Path-Dependent Options One of the most important characteristic is the monitoring frequency of the underlying asset With discrete monitoring, the updating conditions occur only at prefixed dates, e.g. the end of the day/week/month The available analytical formula in the Black-Scholes setting assume continuous monitoring Approximations for the discrete monitoring case are given in Broadie, Glasserman & Kou (1999) Back
50 Lookback Options and Composite Rectangle Formula: Prices m DE CGMY ETime N = 50 N = 100 N = 50 N = 100 N = 50 N = CT CT CS CS Back
51 Lin, Ng & Chan Preconditioner We have to solve the linear system where: I is the identity matrix K is a complex Toeplitz matrix (I + KD)W = Φ D is a diagonal matrix - for example, for the trapezoidal rule D = diag([ 1 2, 1, 1, 1, 1,, 1, 1 2 ])
52 We consider the equivalent preconditioned linear system P 1 (I + KD)W = P 1 Φ Let C be a circulant preconditioners for K, then (I + C) 1 = I + P where the eigenvalues of the circulant matrix P are given by λ 1 + λ and λ are the eigenvalues of C Hence, the preconditioned linear system becomes (I + PD)(I + KD)W = (I + PD)Φ P 1 = I + PD can be computed easily (the preconditioner is already inverted) Back
Pricing Bermudan options in Lévy process models
1 1 Dept. of Industrial & Enterprise Systems Engineering University of Illinois at Urbana-Champaign Joint with Xiong Lin Bachelier Finance Society 6th Congress 6/24/2010 American options in Lévy models
More informationNUMERICAL METHODS OF PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS FOR OPTION PRICE
Trends in Mathematics - New Series Information Center for Mathematical Sciences Volume 13, Number 1, 011, pages 1 5 NUMERICAL METHODS OF PARTIAL INTEGRO-DIFFERENTIAL EQUATIONS FOR OPTION PRICE YONGHOON
More informationTime-changed Brownian motion and option pricing
Time-changed Brownian motion and option pricing Peter Hieber Chair of Mathematical Finance, TU Munich 6th AMaMeF Warsaw, June 13th 2013 Partially joint with Marcos Escobar (RU Toronto), Matthias Scherer
More informationFast Pricing and Calculation of Sensitivities of OTM European Options Under Lévy Processes
Fast Pricing and Calculation of Sensitivities of OTM European Options Under Lévy Processes Sergei Levendorskĭi Jiayao Xie Department of Mathematics University of Leicester Toronto, June 24, 2010 Levendorskĭi
More informationFourier Space Time-stepping Method for Option Pricing with Lévy Processes
FST method Extensions Indifference pricing Fourier Space Time-stepping Method for Option Pricing with Lévy Processes Vladimir Surkov University of Toronto Computational Methods in Finance Conference University
More informationHilbert transform approach for pricing Bermudan options in Lévy models
Hilbert transform approach for pricing Bermudan options in Lévy models 1 1 Dept. of Industrial & Enterprise Systems Engineering University of Illinois at Urbana-Champaign Joint with Xiong Lin Spectral
More informationA Highly Efficient Shannon Wavelet Inverse Fourier Technique for Pricing European Options
A Highly Efficient Shannon Wavelet Inverse Fourier Technique for Pricing European Options Luis Ortiz-Gracia Centre de Recerca Matemàtica (joint work with Cornelis W. Oosterlee, CWI) Models and Numerics
More informationValuation of Equity / FX Instruments
Technical Paper: Valuation of Equity / FX Instruments MathConsult GmbH Altenberger Straße 69 A-4040 Linz, Austria 14 th October, 2009 1 Vanilla Equity Option 1.1 Introduction A vanilla equity option is
More informationIntegrated structural approach to Counterparty Credit Risk with dependent jumps
1/29 Integrated structural approach to Counterparty Credit Risk with dependent jumps, Gianluca Fusai, Daniele Marazzina Cass Business School, Università Piemonte Orientale, Politecnico Milano September
More informationPricing Discretely Monitored Asian Options under Lévy Processes
?? pages 1 Pricing Discretely Monitored Asian Options under Lévy Processes Gianluca Fusai Dipartimento SEMEQ Università degli Studi del Piemonte Orientale Associate Fellow, Financial Options Research Centre,
More informationImplementing Models in Quantitative Finance: Methods and Cases
Gianluca Fusai Andrea Roncoroni Implementing Models in Quantitative Finance: Methods and Cases vl Springer Contents Introduction xv Parti Methods 1 Static Monte Carlo 3 1.1 Motivation and Issues 3 1.1.1
More informationPricing European Options by Stable Fourier-Cosine Series Expansions
Pricing European Options by Stable Fourier-Cosine Series Expansions arxiv:171.886v2 [q-fin.cp] 8 Jan 217 Chunfa Wang Fiance School of Zhejiang University of Finance and Economics, Hangzhou, China, cfwang@zufe.edu.cn
More informationExact Sampling of Jump-Diffusion Processes
1 Exact Sampling of Jump-Diffusion Processes and Dmitry Smelov Management Science & Engineering Stanford University Exact Sampling of Jump-Diffusion Processes 2 Jump-Diffusion Processes Ubiquitous in finance
More informationFinancial Mathematics and Supercomputing
GPU acceleration in early-exercise option valuation Álvaro Leitao and Cornelis W. Oosterlee Financial Mathematics and Supercomputing A Coruña - September 26, 2018 Á. Leitao & Kees Oosterlee SGBM on GPU
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT 08-11 Pricing Early-Exercise and Discrete Barrier Options by Fourier-Cosine Series Expansions F. Fang and C.W. Oosterlee ISSN 1389-6520 Reports of the Department of
More informationUSC Math. Finance April 22, Path-dependent Option Valuation under Jump-diffusion Processes. Alan L. Lewis
USC Math. Finance April 22, 23 Path-dependent Option Valuation under Jump-diffusion Processes Alan L. Lewis These overheads to be posted at www.optioncity.net (Publications) Topics Why jump-diffusion models?
More informationA Continuity Correction under Jump-Diffusion Models with Applications in Finance
A Continuity Correction under Jump-Diffusion Models with Applications in Finance Cheng-Der Fuh 1, Sheng-Feng Luo 2 and Ju-Fang Yen 3 1 Institute of Statistical Science, Academia Sinica, and Graduate Institute
More informationFE610 Stochastic Calculus for Financial Engineers. Stevens Institute of Technology
FE610 Stochastic Calculus for Financial Engineers Lecture 13. The Black-Scholes PDE Steve Yang Stevens Institute of Technology 04/25/2013 Outline 1 The Black-Scholes PDE 2 PDEs in Asset Pricing 3 Exotic
More informationDELFT UNIVERSITY OF TECHNOLOGY
DELFT UNIVERSITY OF TECHNOLOGY REPORT 10-03 Acceleration of Option Pricing Technique on Graphics Processing Units Bowen Zhang Cornelis. W. Oosterlee ISSN 1389-6520 Reports of the Department of Applied
More informationTangent Lévy Models. Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford.
Tangent Lévy Models Sergey Nadtochiy (joint work with René Carmona) Oxford-Man Institute of Quantitative Finance University of Oxford June 24, 2010 6th World Congress of the Bachelier Finance Society Sergey
More informationExtended Libor Models and Their Calibration
Extended Libor Models and Their Calibration Denis Belomestny Weierstraß Institute Berlin Vienna, 16 November 2007 Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Vienna, 16 November
More informationsymmys.com 3.2 Projection of the invariants to the investment horizon
122 3 Modeling the market In the swaption world the underlying rate (3.57) has a bounded range and thus it does not display the explosive pattern typical of a stock price. Therefore the swaption prices
More informationA distributed Laplace transform algorithm for European options
A distributed Laplace transform algorithm for European options 1 1 A. J. Davies, M. E. Honnor, C.-H. Lai, A. K. Parrott & S. Rout 1 Department of Physics, Astronomy and Mathematics, University of Hertfordshire,
More informationImplied Lévy Volatility
Joint work with José Manuel Corcuera, Peter Leoni and Wim Schoutens July 15, 2009 - Eurandom 1 2 The Black-Scholes model The Lévy models 3 4 5 6 7 Delta Hedging at versus at Implied Black-Scholes Volatility
More informationAsian options and meropmorphic Lévy processes
Asian options and meropmorphic Lévy processes July 9 13, 2013 Outline Introduction Existing pricing approaches and the exponential functional I q Methodology overview Meromorphic Lévy processes Theoretical
More informationOption Pricing and Calibration with Time-changed Lévy processes
Option Pricing and Calibration with Time-changed Lévy processes Yan Wang and Kevin Zhang Warwick Business School 12th Feb. 2013 Objectives 1. How to find a perfect model that captures essential features
More informationValuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility model
Valuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility model 1(23) Valuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility
More informationStochastic Grid Bundling Method
Stochastic Grid Bundling Method GPU Acceleration Delft University of Technology - Centrum Wiskunde & Informatica Álvaro Leitao Rodríguez and Cornelis W. Oosterlee London - December 17, 2015 A. Leitao &
More informationComputational Finance
Path Dependent Options Computational Finance School of Mathematics 2018 The Random Walk One of the main assumption of the Black-Scholes framework is that the underlying stock price follows a random walk
More informationPricing and Hedging of Commodity Derivatives using the Fast Fourier Transform
Pricing and Hedging of Commodity Derivatives using the Fast Fourier Transform Vladimir Surkov vladimir.surkov@utoronto.ca Department of Statistical and Actuarial Sciences, University of Western Ontario
More informationLocal vs Non-local Forward Equations for Option Pricing
Local vs Non-local Forward Equations for Option Pricing Rama Cont Yu Gu Abstract When the underlying asset is a continuous martingale, call option prices solve the Dupire equation, a forward parabolic
More informationExtended Libor Models and Their Calibration
Extended Libor Models and Their Calibration Denis Belomestny Weierstraß Institute Berlin Haindorf, 7 Februar 2008 Denis Belomestny (WIAS) Extended Libor Models and Their Calibration Haindorf, 7 Februar
More informationFractional PDE Approach for Numerical Solution of Some Jump-Diffusion Models
Fractional PDE Approach for Numerical Solution of Some Jump-Diffusion Models Andrey Itkin 1 1 HAP Capital and Rutgers University, New Jersey Math Finance and PDE Conference, New Brunswick 2009 A.Itkin
More informationThe stochastic calculus
Gdansk A schedule of the lecture Stochastic differential equations Ito calculus, Ito process Ornstein - Uhlenbeck (OU) process Heston model Stopping time for OU process Stochastic differential equations
More informationCity, University of London Institutional Repository
City Research Online City, University of London Institutional Repository Citation: Fusai, G., Germano, G. and Marazzina, D. (2016). Spitzer identity, Wiener-Hopf factorization and pricing of discretely
More informationThe data-driven COS method
The data-driven COS method Á. Leitao, C. W. Oosterlee, L. Ortiz-Gracia and S. M. Bohte Delft University of Technology - Centrum Wiskunde & Informatica Reading group, March 13, 2017 Reading group, March
More informationBarrier options. In options only come into being if S t reaches B for some 0 t T, at which point they become an ordinary option.
Barrier options A typical barrier option contract changes if the asset hits a specified level, the barrier. Barrier options are therefore path-dependent. Out options expire worthless if S t reaches the
More informationFast narrow bounds on the value of Asian options
Fast narrow bounds on the value of Asian options G. W. P. Thompson Centre for Financial Research, Judge Institute of Management, University of Cambridge Abstract We consider the problem of finding bounds
More informationCS 774 Project: Fall 2009 Version: November 27, 2009
CS 774 Project: Fall 2009 Version: November 27, 2009 Instructors: Peter Forsyth, paforsyt@uwaterloo.ca Office Hours: Tues: 4:00-5:00; Thurs: 11:00-12:00 Lectures:MWF 3:30-4:20 MC2036 Office: DC3631 CS
More informationPricing of some exotic options with N IG-Lévy input
Pricing of some exotic options with N IG-Lévy input Sebastian Rasmus, Søren Asmussen 2 and Magnus Wiktorsson Center for Mathematical Sciences, University of Lund, Box 8, 22 00 Lund, Sweden {rasmus,magnusw}@maths.lth.se
More informationFast and accurate pricing of discretely monitored barrier options by numerical path integration
Comput Econ (27 3:143 151 DOI 1.17/s1614-7-991-5 Fast and accurate pricing of discretely monitored barrier options by numerical path integration Christian Skaug Arvid Naess Received: 23 December 25 / Accepted:
More informationApplication of an Interval Backward Finite Difference Method for Solving the One-Dimensional Heat Conduction Problem
Application of an Interval Backward Finite Difference Method for Solving the One-Dimensional Heat Conduction Problem Malgorzata A. Jankowska 1, Andrzej Marciniak 2 and Tomasz Hoffmann 2 1 Poznan University
More informationAMS Spring Western Sectional Meeting April 4, 2004 (USC campus) Two-sided barrier problems with jump-diffusions Alan L. Lewis
AMS Spring Western Sectional Meeting April 4, 2004 (USC campus) Two-sided barrier problems with jump-diffusions Alan L. Lewis No preprint yet, but overheads to be posted at www.optioncity.net (Publications)
More informationPricing Early-exercise options
Pricing Early-exercise options GPU Acceleration of SGBM method Delft University of Technology - Centrum Wiskunde & Informatica Álvaro Leitao Rodríguez and Cornelis W. Oosterlee Lausanne - December 4, 2016
More informationPricing Algorithms for financial derivatives on baskets modeled by Lévy copulas
Pricing Algorithms for financial derivatives on baskets modeled by Lévy copulas Christoph Winter, ETH Zurich, Seminar for Applied Mathematics École Polytechnique, Paris, September 6 8, 26 Introduction
More informationThe data-driven COS method
The data-driven COS method Á. Leitao, C. W. Oosterlee, L. Ortiz-Gracia and S. M. Bohte Delft University of Technology - Centrum Wiskunde & Informatica CMMSE 2017, July 6, 2017 Álvaro Leitao (CWI & TUDelft)
More informationStochastic Volatility (Working Draft I)
Stochastic Volatility (Working Draft I) Paul J. Atzberger General comments or corrections should be sent to: paulatz@cims.nyu.edu 1 Introduction When using the Black-Scholes-Merton model to price derivative
More informationNumerical Solution of Two Asset Jump Diffusion Models for Option Valuation
Numerical Solution of Two Asset Jump Diffusion Models for Option Valuation Simon S. Clift and Peter A. Forsyth Original: December 5, 2005 Revised: January 31, 2007 Abstract Under the assumption that two
More informationOption pricing in the stochastic volatility model of Barndorff-Nielsen and Shephard
Option pricing in the stochastic volatility model of Barndorff-Nielsen and Shephard Indifference pricing and the minimal entropy martingale measure Fred Espen Benth Centre of Mathematics for Applications
More information17 Lookback Options: A Discrete Problem *
17 Lookback Options: A Discrete Problem * Key words: PDE, Monte Carlo simulation, transform method, discrete monitoring, exotic option The number of exotic options traded on the market has dramatically
More informationEFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS
Commun. Korean Math. Soc. 23 (2008), No. 2, pp. 285 294 EFFICIENT MONTE CARLO ALGORITHM FOR PRICING BARRIER OPTIONS Kyoung-Sook Moon Reprinted from the Communications of the Korean Mathematical Society
More informationM5MF6. Advanced Methods in Derivatives Pricing
Course: Setter: M5MF6 Dr Antoine Jacquier MSc EXAMINATIONS IN MATHEMATICS AND FINANCE DEPARTMENT OF MATHEMATICS April 2016 M5MF6 Advanced Methods in Derivatives Pricing Setter s signature...........................................
More information2.1 Mathematical Basis: Risk-Neutral Pricing
Chapter Monte-Carlo Simulation.1 Mathematical Basis: Risk-Neutral Pricing Suppose that F T is the payoff at T for a European-type derivative f. Then the price at times t before T is given by f t = e r(t
More informationModel Estimation. Liuren Wu. Fall, Zicklin School of Business, Baruch College. Liuren Wu Model Estimation Option Pricing, Fall, / 16
Model Estimation Liuren Wu Zicklin School of Business, Baruch College Fall, 2007 Liuren Wu Model Estimation Option Pricing, Fall, 2007 1 / 16 Outline 1 Statistical dynamics 2 Risk-neutral dynamics 3 Joint
More information"Pricing Exotic Options using Strong Convergence Properties
Fourth Oxford / Princeton Workshop on Financial Mathematics "Pricing Exotic Options using Strong Convergence Properties Klaus E. Schmitz Abe schmitz@maths.ox.ac.uk www.maths.ox.ac.uk/~schmitz Prof. Mike
More informationEquity correlations implied by index options: estimation and model uncertainty analysis
1/18 : estimation and model analysis, EDHEC Business School (joint work with Rama COT) Modeling and managing financial risks Paris, 10 13 January 2011 2/18 Outline 1 2 of multi-asset models Solution to
More informationMATH3075/3975 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS
MATH307/37 FINANCIAL MATHEMATICS TUTORIAL PROBLEMS School of Mathematics and Statistics Semester, 04 Tutorial problems should be used to test your mathematical skills and understanding of the lecture material.
More informationOption Pricing for a Stochastic-Volatility Jump-Diffusion Model
Option Pricing for a Stochastic-Volatility Jump-Diffusion Model Guoqing Yan and Floyd B. Hanson Department of Mathematics, Statistics, and Computer Science University of Illinois at Chicago Conference
More informationNumerical valuation for option pricing under jump-diffusion models by finite differences
Numerical valuation for option pricing under jump-diffusion models by finite differences YongHoon Kwon Younhee Lee Department of Mathematics Pohang University of Science and Technology June 23, 2010 Table
More informationValuation of European Call Option via Inverse Fourier Transform
ISSN 2255-9094 (online) ISSN 2255-9086 (print) December 2017, vol. 20, pp. 91 96 doi: 10.1515/itms-2017-0016 https://www.degruyter.com/view/j/itms Valuation of European Call Option via Inverse Fourier
More informationSolving the Black-Scholes Equation
Solving the Black-Scholes Equation An Undergraduate Introduction to Financial Mathematics J. Robert Buchanan 2010 Initial Value Problem for the European Call rf = F t + rsf S + 1 2 σ2 S 2 F SS for (S,
More informationBarrier option pricing for assets with Markov-modulated dividends
Barrier option pricing for assets with Markov-modulated dividends Giuseppe Di Graziano Statistical Laboratory University of Cambridge L.C.G Rogers Statistical Laboratory University of Cambridge 19th April
More informationFrom Characteristic Functions and Fourier Transforms to PDFs/CDFs and Option Prices
From Characteristic Functions and Fourier Transforms to PDFs/CDFs and Option Prices Liuren Wu Zicklin School of Business, Baruch College Fall, 2007 Liuren Wu Fourier Transforms Option Pricing, Fall, 2007
More informationNormal Inverse Gaussian (NIG) Process
With Applications in Mathematical Finance The Mathematical and Computational Finance Laboratory - Lunch at the Lab March 26, 2009 1 Limitations of Gaussian Driven Processes Background and Definition IG
More information1.1 Basic Financial Derivatives: Forward Contracts and Options
Chapter 1 Preliminaries 1.1 Basic Financial Derivatives: Forward Contracts and Options A derivative is a financial instrument whose value depends on the values of other, more basic underlying variables
More informationPricing Variance Swaps under Stochastic Volatility Model with Regime Switching - Discrete Observations Case
Pricing Variance Swaps under Stochastic Volatility Model with Regime Switching - Discrete Observations Case Guang-Hua Lian Collaboration with Robert Elliott University of Adelaide Feb. 2, 2011 Robert Elliott,
More informationRichardson Extrapolation Techniques for the Pricing of American-style Options
Richardson Extrapolation Techniques for the Pricing of American-style Options June 1, 2005 Abstract Richardson Extrapolation Techniques for the Pricing of American-style Options In this paper we re-examine
More informationChapter 5 Finite Difference Methods. Math6911 W07, HM Zhu
Chapter 5 Finite Difference Methods Math69 W07, HM Zhu References. Chapters 5 and 9, Brandimarte. Section 7.8, Hull 3. Chapter 7, Numerical analysis, Burden and Faires Outline Finite difference (FD) approximation
More informationOption Pricing for Discrete Hedging and Non-Gaussian Processes
Option Pricing for Discrete Hedging and Non-Gaussian Processes Kellogg College University of Oxford A thesis submitted in partial fulfillment of the requirements for the MSc in Mathematical Finance November
More informationFinal Exam Key, JDEP 384H, Spring 2006
Final Exam Key, JDEP 384H, Spring 2006 Due Date for Exam: Thursday, May 4, 12:00 noon. Instructions: Show your work and give reasons for your answers. Write out your solutions neatly and completely. There
More informationSolving the Black-Scholes Equation
Solving the Black-Scholes Equation An Undergraduate Introduction to Financial Mathematics J. Robert Buchanan 2014 Initial Value Problem for the European Call The main objective of this lesson is solving
More informationFractional Black - Scholes Equation
Chapter 6 Fractional Black - Scholes Equation 6.1 Introduction The pricing of options is a central problem in quantitative finance. It is both a theoretical and practical problem since the use of options
More informationThe Binomial Model. Chapter 3
Chapter 3 The Binomial Model In Chapter 1 the linear derivatives were considered. They were priced with static replication and payo tables. For the non-linear derivatives in Chapter 2 this will not work
More informationHints on Some of the Exercises
Hints on Some of the Exercises of the book R. Seydel: Tools for Computational Finance. Springer, 00/004/006/009/01. Preparatory Remarks: Some of the hints suggest ideas that may simplify solving the exercises
More informationBarrier Option. 2 of 33 3/13/2014
FPGA-based Reconfigurable Computing for Pricing Multi-Asset Barrier Options RAHUL SRIDHARAN, GEORGE COOKE, KENNETH HILL, HERMAN LAM, ALAN GEORGE, SAAHPC '12, PROCEEDINGS OF THE 2012 SYMPOSIUM ON APPLICATION
More informationMath 623 (IOE 623), Winter 2008: Final exam
Math 623 (IOE 623), Winter 2008: Final exam Name: Student ID: This is a closed book exam. You may bring up to ten one sided A4 pages of notes to the exam. You may also use a calculator but not its memory
More informationSparse Wavelet Methods for Option Pricing under Lévy Stochastic Volatility models
Sparse Wavelet Methods for Option Pricing under Lévy Stochastic Volatility models Norbert Hilber Seminar of Applied Mathematics ETH Zürich Workshop on Financial Modeling with Jump Processes p. 1/18 Outline
More informationNumerical schemes for SDEs
Lecture 5 Numerical schemes for SDEs Lecture Notes by Jan Palczewski Computational Finance p. 1 A Stochastic Differential Equation (SDE) is an object of the following type dx t = a(t,x t )dt + b(t,x t
More informationPricing Barrier Options under Local Volatility
Abstract Pricing Barrier Options under Local Volatility Artur Sepp Mail: artursepp@hotmail.com, Web: www.hot.ee/seppar 16 November 2002 We study pricing under the local volatility. Our research is mainly
More informationOptimal Hedging of Variance Derivatives. John Crosby. Centre for Economic and Financial Studies, Department of Economics, Glasgow University
Optimal Hedging of Variance Derivatives John Crosby Centre for Economic and Financial Studies, Department of Economics, Glasgow University Presentation at Baruch College, in New York, 16th November 2010
More informationSTOCHASTIC VOLATILITY MODELS: CALIBRATION, PRICING AND HEDGING. Warrick Poklewski-Koziell
STOCHASTIC VOLATILITY MODELS: CALIBRATION, PRICING AND HEDGING by Warrick Poklewski-Koziell Programme in Advanced Mathematics of Finance School of Computational and Applied Mathematics University of the
More informationTEST OF BOUNDED LOG-NORMAL PROCESS FOR OPTIONS PRICING
TEST OF BOUNDED LOG-NORMAL PROCESS FOR OPTIONS PRICING Semih Yön 1, Cafer Erhan Bozdağ 2 1,2 Department of Industrial Engineering, Istanbul Technical University, Macka Besiktas, 34367 Turkey Abstract.
More informationMonte Carlo Simulation of Stochastic Processes
Monte Carlo Simulation of Stochastic Processes Last update: January 10th, 2004. In this section is presented the steps to perform the simulation of the main stochastic processes used in real options applications,
More informationA NOTE ON NUMERICAL SOLUTION OF A LINEAR BLACK-SCHOLES MODEL
GANIT J. Bangladesh Math. Soc. (ISSN 1606-3694) Vol. 33 (2013) 103-115 A NOTE ON NUMERICAL SOLUTION OF A LINEAR BLACK-SCHOLES MODEL Md. Kazi Salah Uddin 1*, Mostak Ahmed 2 and Samir Kumar Bhowmik 1,3 1
More informationChapter 3: Black-Scholes Equation and Its Numerical Evaluation
Chapter 3: Black-Scholes Equation and Its Numerical Evaluation 3.1 Itô Integral 3.1.1 Convergence in the Mean and Stieltjes Integral Definition 3.1 (Convergence in the Mean) A sequence {X n } n ln of random
More informationSimulating Stochastic Differential Equations
IEOR E4603: Monte-Carlo Simulation c 2017 by Martin Haugh Columbia University Simulating Stochastic Differential Equations In these lecture notes we discuss the simulation of stochastic differential equations
More informationFRAME AND FOURIER METHODS FOR EXOTIC OPTION PRICING AND HEDGING
FRAME AND FOURIER METHODS FOR EXOTIC OPTION PRICING AND HEDGING A Thesis Presented to The Academic Faculty by Justin Lars Kirkby In Partial Fulfillment of the Requirements for the Degree Doctor of Philosophy
More informationQueens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Fall 2017 Instructor: Dr. Sateesh Mane.
Queens College, CUNY, Department of Computer Science Computational Finance CSCI 365 / 765 Fall 217 Instructor: Dr. Sateesh Mane c Sateesh R. Mane 217 13 Lecture 13 November 15, 217 Derivation of the Black-Scholes-Merton
More informationTwo-dimensional COS method
Two-dimensional COS method Marjon Ruijter Winterschool Lunteren 22 January 2013 1/29 Introduction PhD student since October 2010 Prof.dr.ir. C.W. Oosterlee). CWI national research center for mathematics
More informationFinancial Risk Management
Financial Risk Management Professor: Thierry Roncalli Evry University Assistant: Enareta Kurtbegu Evry University Tutorial exercices #4 1 Correlation and copulas 1. The bivariate Gaussian copula is given
More informationLecture 4. Finite difference and finite element methods
Finite difference and finite element methods Lecture 4 Outline Black-Scholes equation From expectation to PDE Goal: compute the value of European option with payoff g which is the conditional expectation
More informationChapter 14. The Multi-Underlying Black-Scholes Model and Correlation
Chapter 4 The Multi-Underlying Black-Scholes Model and Correlation So far we have discussed single asset options, the payoff function depended only on one underlying. Now we want to allow multiple underlyings.
More informationSADDLEPOINT APPROXIMATIONS TO OPTION PRICES 1. By L. C. G. Rogers and O. Zane University of Bath and First Chicago NBD
The Annals of Applied Probability 1999, Vol. 9, No. 2, 493 53 SADDLEPOINT APPROXIMATIONS TO OPTION PRICES 1 By L. C. G. Rogers and O. Zane University of Bath and First Chicago NBD The use of saddlepoint
More informationMath 416/516: Stochastic Simulation
Math 416/516: Stochastic Simulation Haijun Li lih@math.wsu.edu Department of Mathematics Washington State University Week 13 Haijun Li Math 416/516: Stochastic Simulation Week 13 1 / 28 Outline 1 Simulation
More informationKing s College London
King s College London University Of London This paper is part of an examination of the College counting towards the award of a degree. Examinations are governed by the College Regulations under the authority
More informationarxiv: v1 [math.st] 18 Sep 2018
Gram Charlier and Edgeworth expansion for sample variance arxiv:809.06668v [math.st] 8 Sep 08 Eric Benhamou,* A.I. SQUARE CONNECT, 35 Boulevard d Inkermann 900 Neuilly sur Seine, France and LAMSADE, Universit
More information1 The continuous time limit
Derivative Securities, Courant Institute, Fall 2008 http://www.math.nyu.edu/faculty/goodman/teaching/derivsec08/index.html Jonathan Goodman and Keith Lewis Supplementary notes and comments, Section 3 1
More informationPricing Dynamic Guaranteed Funds Under a Double Exponential. Jump Diffusion Process. Chuang-Chang Chang, Ya-Hui Lien and Min-Hung Tsay
Pricing Dynamic Guaranteed Funds Under a Double Exponential Jump Diffusion Process Chuang-Chang Chang, Ya-Hui Lien and Min-Hung Tsay ABSTRACT This paper complements the extant literature to evaluate the
More informationLarge Deviations and Stochastic Volatility with Jumps: Asymptotic Implied Volatility for Affine Models
Large Deviations and Stochastic Volatility with Jumps: TU Berlin with A. Jaquier and A. Mijatović (Imperial College London) SIAM conference on Financial Mathematics, Minneapolis, MN July 10, 2012 Implied
More informationUsing radial basis functions for option pricing
Using radial basis functions for option pricing Elisabeth Larsson Division of Scientific Computing Department of Information Technology Uppsala University Actuarial Mathematics Workshop, March 19, 2013,
More information