Fourier Space Time-stepping Method for Option Pricing with Lévy Processes
|
|
- Adelia Gibbs
- 5 years ago
- Views:
Transcription
1 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 of Waterloo July 27, 2007 Joint work with Ken Jackson and Sebastian Jaimungal, University of Toronto 1 / 32
2 FST method Extensions Indifference pricing 1 Fourier Space Time-stepping method Infinitesimal generator and characteristic exponent Method derivation Numerical results 2 Extensions Multi-asset options Regime switching 3 Indifference pricing Optimal investment problems Application of FST to solution of HJB equations 2 / 32
3 FST method Extensions Indifference pricing Characteristic exponent Method derivation Numerical results The Option Pricing Problem Option payoff is given by ϕ(s) Stock price follows an exponential Lévy model: S(t) = S(0)e µt+x (t), X (t) is a Lévy process 3 / 32
4 FST method Extensions Indifference pricing Characteristic exponent Method derivation Numerical results The Option Pricing Problem Option payoff is given by ϕ(s) Stock price follows an exponential Lévy model: S(t) = S(0)e µt+x (t), X (t) is a Lévy process Generalizing PIDE for Lévy processes { t v + Lv = 0 v(t, x) = ϕ(s(0) e x ) where L is the infinitesimal generator of the Lévy process: Lf = γ x f + σ2 2 [ xxf + f (x + y) f (x) y1 { y <1} x f (x) ] ν(dy) R/{0} 3 / 32
5 FST method Extensions Indifference pricing Characteristic exponent Method derivation Numerical results Finite Difference Methods for Option Pricing Alternating Direction Implicit-FFT - Andersen and Andreasen (2000) Implicit-Explicit (IMEX) - Cont and Tankov (2004) IMEX Runge-Kutta - Briani, Natalini, and Russo (2004) Fixed Point Iteration - d Halluin, Forsyth, and Vetzal (2005) 4 / 32
6 FST method Extensions Indifference pricing Characteristic exponent Method derivation Numerical results Finite Difference Methods for Option Pricing Alternating Direction Implicit-FFT - Andersen and Andreasen (2000) Implicit-Explicit (IMEX) - Cont and Tankov (2004) IMEX Runge-Kutta - Briani, Natalini, and Russo (2004) Fixed Point Iteration - d Halluin, Forsyth, and Vetzal (2005) Common Features: Treat the integral term explicitly to avoid solving a dense system of linear equations. Use the Fast Fourier Transform (FFT) to speed up the computation of the integral term (which can be regarded as a convolution) 4 / 32
7 FST method Extensions Indifference pricing Characteristic exponent Method derivation Numerical results Finite Difference Methods for Option Pricing Alternating Direction Implicit-FFT - Andersen and Andreasen (2000) Implicit-Explicit (IMEX) - Cont and Tankov (2004) IMEX Runge-Kutta - Briani, Natalini, and Russo (2004) Fixed Point Iteration - d Halluin, Forsyth, and Vetzal (2005) Common Features: Treat the integral term explicitly to avoid solving a dense system of linear equations. Use the Fast Fourier Transform (FFT) to speed up the computation of the integral term (which can be regarded as a convolution) Drawbacks: Diffusive and integral terms are treated asymmetrically Large jump are truncated and small jumps approximated by diffusion Difficult to extend to higher dimensions 4 / 32
8 FST method Extensions Indifference pricing Characteristic exponent Method derivation Numerical results Infinitesimal Generator and Characteristic Exponent The characteristic exponent of a Lévy-Khinchin representation can be factored from a Fourier transform of the operator (Sato 1999) F[Lv](τ, ω) = {iγω σ2 ω 2 } + [e iωx 1 iωy1 { ω <1} ]ν(dy) F[v](τ, ω) 2 = ψ(ω)f[v](τ, ω) Model Characteristic Exponent ψ(ω) Black-Scholes-Merton iµω σ2 ω 2 2 Merton Jump-Diffusion iµω σ2 ω λ(e i µω σ2 ω 2 /2 1) Variance Gamma 1 κ log(1 iµκω + σ2 κω 2 2 ) CGMY CΓ( Y ) [ (M iω) Y M Y + (G +iω) Y G Y ] 5 / 32
9 FST method Extensions Indifference pricing Characteristic exponent Method derivation Numerical results Numerical Method Derivation Apply the Fourier transform to the pricing PIDE { t F[v](t, ω) + Ψ(ω)F[v](t, ω) = 0, F[v](T, ω) = F[ϕ](ω) 6 / 32
10 FST method Extensions Indifference pricing Characteristic exponent Method derivation Numerical results Numerical Method Derivation Apply the Fourier transform to the pricing PIDE { t F[v](t, ω) + Ψ(ω)F[v](t, ω) = 0, F[v](T, ω) = F[ϕ](ω) Resulting ODE has explicit solution F[v](t 1, ω) = F[v](t 2, ω) e (t2 t1)ψ(ω) 6 / 32
11 FST method Extensions Indifference pricing Characteristic exponent Method derivation Numerical results Numerical Method Derivation Apply the Fourier transform to the pricing PIDE { t F[v](t, ω) + Ψ(ω)F[v](t, ω) = 0, F[v](T, ω) = F[ϕ](ω) Resulting ODE has explicit solution F[v](t 1, ω) = F[v](t 2, ω) e (t2 t1)ψ(ω) Apply the inverse Fourier transform v(t 1, x) = F 1 { F[v](t 2, ω) e (t2 t1)ψ(ω)} (x) 6 / 32
12 FST method Extensions Indifference pricing Characteristic exponent Method derivation Numerical results Numerical Method Derivation Apply the Fourier transform to the pricing PIDE { t F[v](t, ω) + Ψ(ω)F[v](t, ω) = 0, F[v](T, ω) = F[ϕ](ω) Resulting ODE has explicit solution F[v](t 1, ω) = F[v](t 2, ω) e (t2 t1)ψ(ω) Apply the inverse Fourier transform v(t 1, x) = F 1 { F[v](t 2, ω) e (t2 t1)ψ(ω)} (x) Fourier Space Time-stepping (FST) method v n 1 = FFT 1 [FFT[v n ] e Ψ t ] 6 / 32
13 FST method Extensions Indifference pricing Characteristic exponent Method derivation Numerical results European Call Option N Value Change log 2 Ratio CPU-Time Option: European call S = 1.0, K = 1.0, T = 0.2 Model: Kou jump-diffusion σ = 0.2, λ = 0.2, p = 0.5, η = 3, η + = 2, r = 0.0 Quoted Price: Almendral and Oosterlee (2005) 7 / 32
14 FST method Extensions Indifference pricing Characteristic exponent Method derivation Numerical results American Put Option N M Value Change log 2 Ratio CPU-Time Option: American put S = 90.0, K = 98.0, T = 0.25 Model: CGMY C = 0.42, G = 4.37, M = 191.2, Y = , r = 0.06 Quoted Price: Forsyth, Wan, and Wang (2006) 8 / 32
15 FST method Extensions Indifference pricing Characteristic exponent Method derivation Numerical results American Put Option N M Value Change log 2 Ratio CPU-Time Option: American put S = 10.0, K = 10.0, T = 0.25 Model: CGMY C = 1.0, G = 8.8, M = 9.2, Y = 1.8, r = / 32
16 FST method Extensions Indifference pricing Characteristic exponent Method derivation Numerical results Up-and-Out Barrier Call Option N M Value Change log 2 Ratio CPU-Time Option: Up-and-Out Barrier Call S = 100.0, K = 100.0, B = 110, T = 1.0 Model: Black-Scholes-Merton σ = 0.15, r = 0.05, q = 0.02 Closed-Form Price: Hull (2005) 10 / 32
17 FST method Extensions Indifference pricing Multi-asset options Regime switching Multi-asset options Pricing PIDE { ( t + L) v = 0, v(t, x) = ϕ(s(0) e x ) 11 / 32
18 FST method Extensions Indifference pricing Multi-asset options Regime switching Multi-asset options Pricing PIDE { ( t + L) v = 0, v(t, x) = ϕ(s(0) e x ) Multi-dimensional FST method v n 1 = FFT 1 [FFT[v n ] e Ψ t ] 11 / 32
19 FST method Extensions Indifference pricing Multi-asset options Regime switching Multi-asset options Pricing PIDE { ( t + L) v = 0, v(t, x) = ϕ(s(0) e x ) Multi-dimensional FST method v n 1 = FFT 1 [FFT[v n ] e Ψ t ] Jackson, Jaimungal and Surkov (2007) discuss application of FST to pricing of spread, American spread and catastrophe equity put options 11 / 32
20 FST method Extensions Indifference pricing Multi-asset options Regime switching Catastrophe Equity Put (CatEPut) In the event of large (catastrophic) losses U, the insurer receives a put option on its own stock ϕ(s(t ), L T ) = 1 LT >U (K S(T )) + 12 / 32
21 FST method Extensions Indifference pricing Multi-asset options Regime switching Catastrophe Equity Put (CatEPut) Presence of losses drives the share value down, not an independent jump process S(t) = S(0) exp { α L(t) + γ t + σ W t } L(t) = N(t) l i n=1 13 / 32
22 FST method Extensions Indifference pricing Multi-asset options Regime switching Catastrophe Equity Put (CatEPut) Presence of losses drives the share value down, not an independent jump process S(t) = S(0) exp { α L(t) + γ t + σ W t } L(t) = N(t) l i n=1 When losses are modeled as a Gamma r.v., characteristic exponent is [ ] ( Ψ(ω 1, ω 2 ) = i γ ω σ2 ω1 2 + λ 1 i( αω 1 + ω 2 ) v ) m 2 m ) v 1 13 / 32
23 FST method Extensions Indifference pricing Multi-asset options Regime switching European CatEPut 14 / 32
24 FST method Extensions Indifference pricing Multi-asset options Regime switching American CatEPut 15 / 32
25 FST method Extensions Indifference pricing Multi-asset options Regime switching Regime Switching Let K := {1,..., K} denote the possible hidden states of the world, driven by a continuous time Markov chain Z t. The transition probability from state k at time t 1 to state l at time t 2 is given by P t1t2 kl = Q(Z t2 = l Z t1 = k) = (exp{(t 2 t 1 )A}) kl where A is the Markov chain generator Within state k, log-stock follows Lévy model k 16 / 32
26 FST method Extensions Indifference pricing Multi-asset options Regime switching Regime Switching Pricing PIDE { [ t + ( A kk + L (k))] v(x, k, t) + j k A jk v(x, j, t) = 0 v(x, k, T ) = ϕ(s(0)e x ) 17 / 32
27 FST method Extensions Indifference pricing Multi-asset options Regime switching Regime Switching Pricing PIDE { [ t + ( A kk + L (k))] v(x, k, t) + j k A jk v(x, j, t) = 0 v(x, k, T ) = ϕ(s(0)e x ) FST Method where v n 1 = FFT 1 [FFT[v n ] e Ψ t ] { Akk + Ψ Ψ(ω) kl = (k) (ω), k = l A kl, k l No time-stepping required for European options 17 / 32
28 FST method Extensions Indifference pricing Multi-asset options Regime switching American put option with 3 regimes N M Value Change log 2 Ratio CPU-Time Option: American put S = 100.0, K = 100.0, T = 1.0 Model: Merton jump-diffusion σ = 0.15, µ = 0.5, σ = 0.45, r = 0.05, q = 0.02, λ [0.3, 0.5, 0.7], p = [0.2, 0.3, 0.5], A = [ 0.8, 0.6, 0.2; 0.2, 1, 0.8; 0.1, 0.3, 0.4] 18 / 32
29 FST method Extensions Indifference pricing Multi-asset options Regime switching Exercise boundary for American put option with 3 regimes 19 / 32
30 FST method Extensions Indifference pricing Multi-asset options Regime switching FST Method Strengths Stable and robust, even for options with discontinuous payoffs Easily extendable to various stochastic processes and no loss of performance for infinite activity processes Can be applied to multi-dimensional and regime-switching problems in a natural manner Computationally efficient Computational cost is O(MNlogN) while the error is O( x 2 + t 2 ) European options priced in a single time-step Bermudan style options do not require time-stepping between monitoring dates 20 / 32
31 FST method Extensions Indifference pricing Optimal investment problems Application of FST to HJB Indifference Pricing In incomplete markets (non-traded assets, transaction costs, portfolio constraints, etc.) perfect replication is impossible Investor can still maximize the expected utility of wealth through dynamic trading The price of a claim is the initial wealth forgone so that the investor is no worse off in expected utility terms at maturity The framework incorporates wealth dependence, non-linear pricing and risk-aversion 21 / 32
32 FST method Extensions Indifference pricing Optimal investment problems Application of FST to HJB The Modelling Framework Utility function is a strictly increasing and concave function ranking investor s preferences of wealth Popular choices include CRRA model U(x) = x 1 γ 1 γ or CARA model U(x) = 1 1 γ e γx 22 / 32
33 FST method Extensions Indifference pricing Optimal investment problems Application of FST to HJB The Modelling Framework Utility function is a strictly increasing and concave function ranking investor s preferences of wealth Popular choices include CRRA model U(x) = x 1 γ 1 γ or CARA model U(x) = 1 1 γ e γx An economic agent over a fixed-time horizon attempts to optimally allocate his investment between risky (S t ) and risk-free (B t ) assets ds t = µs t dt + σs t dw t dx t = π t X t ds t S t + (1 π t )X t db t B t = (π t (µ r) + r)x t dt + π t σx t dw t where π t is the share of wealth allocated in stocks 22 / 32
34 FST method Extensions Indifference pricing Optimal investment problems Application of FST to HJB Two Optimal Control Problems Investor maximizes the expected utility of wealth at time horizon T, given initial wealth x at t: V 0 (t, x) = sup E [U(XT π X t = x)] π t 23 / 32
35 FST method Extensions Indifference pricing Optimal investment problems Application of FST to HJB Two Optimal Control Problems Investor maximizes the expected utility of wealth at time horizon T, given initial wealth x at t: V 0 (t, x) = sup E [U(XT π X t = x)] π t In addition to initial endowment x, investor receives k derivative contracts with payoff C(S T ): V k (t, x, s) = sup E [U(XT π + k C(S T ) X t = x, S t = s)] π t 23 / 32
36 FST method Extensions Indifference pricing Optimal investment problems Application of FST to HJB Two Optimal Control Problems Investor maximizes the expected utility of wealth at time horizon T, given initial wealth x at t: V 0 (t, x) = sup E [U(XT π X t = x)] π t In addition to initial endowment x, investor receives k derivative contracts with payoff C(S T ): V k (t, x, s) = sup E [U(XT π + k C(S T ) X t = x, S t = s)] π t Indifference Pricing Principle Indifference buy price pbuy k (s) and sell price pk sell (s) satisfy V 0 (t, x) = V k (t, x p k buy (s), s) V 0 (t, x) = V k (t, x + p k sell(s), s) 23 / 32
37 FST method Extensions Indifference pricing Optimal investment problems Application of FST to HJB Optimal Investment Problem Optimal value HJB equation { t V 0 (t, x) + sup π t {A π V 0 (t, x)} = 0 where V 0 (T, x) = U(x) A π = (π t (µ r) + r)x x π2 t σ 2 x 2 xx 24 / 32
38 FST method Extensions Indifference pricing Optimal investment problems Application of FST to HJB Optimal Investment Problem Optimal value HJB equation { t V 0 (t, x) + sup π t {A π V 0 (t, x)} = 0 where V 0 (T, x) = U(x) A π = (π t (µ r) + r)x x π2 t σ 2 x 2 xx Example: CARA utility γ = 0.8; GBM µ = 0.1, σ = 0.2, r = / 32
39 FST method Extensions Indifference pricing Optimal investment problems Application of FST to HJB Optimal Expected Utility 25 / 32
40 FST method Extensions Indifference pricing Optimal investment problems Application of FST to HJB Optimal Investment Problem with Option Optimal value 2D HJB equation { t V k (t, x, s) + sup π t {A π V k (t, x, s)} = 0 where V k (T, x, s) = U(x + k C(s)) A π = (π t (µ r) + r)x x π2 t σ 2 x 2 xx + µs s σ2 ss ss + πσ 2 xs xs 26 / 32
41 FST method Extensions Indifference pricing Optimal investment problems Application of FST to HJB Optimal Investment Problem with Option Optimal value 2D HJB equation { t V k (t, x, s) + sup π t {A π V k (t, x, s)} = 0 where V k (T, x, s) = U(x + k C(s)) A π = (π t (µ r) + r)x x π2 t σ 2 x 2 xx + µs s σ2 ss ss + πσ 2 xs xs Example continued: European put K = 2, T = 5; k = 2 26 / 32
42 FST method Extensions Indifference pricing Optimal investment problems Application of FST to HJB Optimal Expected Utility (t=5) 27 / 32
43 FST method Extensions Indifference pricing Optimal investment problems Application of FST to HJB Optimal Expected Utility (t=0) 28 / 32
44 FST method Extensions Indifference pricing Optimal investment problems Application of FST to HJB Indifference Buy Price 29 / 32
45 FST method Extensions Indifference pricing Optimal investment problems Application of FST to HJB Application of FST to solution of HJB equations General Approach: Fix π t (x) over the time-step [t 1, t 2 ] and solve the resulting PIDE Iterate to converge to the optimal policy 30 / 32
46 FST method Extensions Indifference pricing Optimal investment problems Application of FST to HJB Application of FST to solution of HJB equations General Approach: Fix π t (x) over the time-step [t 1, t 2 ] and solve the resulting PIDE Iterate to converge to the optimal policy Optimization Algorithm 1: Fix π t (x) = π t over the entire space F[A π ] has an analytic form and the resulting ODE can be solved explicitly (just like in option pricing) Inefficient if optimal policy is uniformly distributed in space 30 / 32
47 FST method Extensions Indifference pricing Optimal investment problems Application of FST to HJB Application of FST to solution of HJB equations General Approach: Fix π t (x) over the time-step [t 1, t 2 ] and solve the resulting PIDE Iterate to converge to the optimal policy Optimization Algorithm 1: Fix π t (x) = π t over the entire space F[A π ] has an analytic form and the resulting ODE can be solved explicitly (just like in option pricing) Inefficient if optimal policy is uniformly distributed in space Optimization Algorithm 2: Let π t (x) vary in space F[A π ] involves convolutions of F[A π ] and F[π(x)], F[π 2 (x)] Can apply a policy iteration approach of Wang and Forsyth (2006) Working in Fourier space we can solve an Integral HJB 30 / 32
48 FST method Extensions Indifference pricing Future Work Exotic, multi-asset options Mean reverting processes in energy markets HJB equations arising from optimal control problems Efficient policy iteration algorithm Optimal control with jump-diffusions 31 / 32
49 FST method Extensions Indifference pricing Thank You! 32 / 32
Pricing 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 informationOptimal investments under dynamic performance critria. Lecture IV
Optimal investments under dynamic performance critria Lecture IV 1 Utility-based measurement of performance 2 Deterministic environment Utility traits u(x, t) : x wealth and t time Monotonicity u x (x,
More informationControl Improvement for Jump-Diffusion Processes with Applications to Finance
Control Improvement for Jump-Diffusion Processes with Applications to Finance Nicole Bäuerle joint work with Ulrich Rieder Toronto, June 2010 Outline Motivation: MDPs Controlled Jump-Diffusion Processes
More informationMulti-period mean variance asset allocation: Is it bad to win the lottery?
Multi-period mean variance asset allocation: Is it bad to win the lottery? Peter Forsyth 1 D.M. Dang 1 1 Cheriton School of Computer Science University of Waterloo Guangzhou, July 28, 2014 1 / 29 The Basic
More informationBACHELIER FINANCE SOCIETY. 4 th World Congress Tokyo, Investments and forward utilities. Thaleia Zariphopoulou The University of Texas at Austin
BACHELIER FINANCE SOCIETY 4 th World Congress Tokyo, 26 Investments and forward utilities Thaleia Zariphopoulou The University of Texas at Austin 1 Topics Utility-based measurement of performance Utilities
More informationPDE Methods for Option Pricing under Jump Diffusion Processes
PDE Methods for Option Pricing under Jump Diffusion Processes Prof Kevin Parrott University of Greenwich November 2009 Typeset by FoilTEX Summary Merton jump diffusion American options Levy Processes -
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 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 informationHigh Frequency Trading in a Regime-switching Model. Yoontae Jeon
High Frequency Trading in a Regime-switching Model by Yoontae Jeon A thesis submitted in conformity with the requirements for the degree of Master of Science Graduate Department of Mathematics University
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 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 informationAll Investors are Risk-averse Expected Utility Maximizers. Carole Bernard (UW), Jit Seng Chen (GGY) and Steven Vanduffel (Vrije Universiteit Brussel)
All Investors are Risk-averse Expected Utility Maximizers Carole Bernard (UW), Jit Seng Chen (GGY) and Steven Vanduffel (Vrije Universiteit Brussel) First Name: Waterloo, April 2013. Last Name: UW ID #:
More informationOptimal asset allocation under forward performance criteria Oberwolfach, February 2007
Optimal asset allocation under forward performance criteria Oberwolfach, February 2007 Thaleia Zariphopoulou The University of Texas at Austin 1 References Indifference valuation in binomial models (with
More informationIncorporating Managerial Cash-Flow Estimates and Risk Aversion to Value Real Options Projects. The Fields Institute for Mathematical Sciences
Incorporating Managerial Cash-Flow Estimates and Risk Aversion to Value Real Options Projects The Fields Institute for Mathematical Sciences Sebastian Jaimungal sebastian.jaimungal@utoronto.ca Yuri Lawryshyn
More informationUtility Indifference Pricing and Dynamic Programming Algorithm
Chapter 8 Utility Indifference ricing and Dynamic rogramming Algorithm In the Black-Scholes framework, we can perfectly replicate an option s payoff. However, it may not be true beyond the Black-Scholes
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 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 informationSPDE and portfolio choice (joint work with M. Musiela) Princeton University. Thaleia Zariphopoulou The University of Texas at Austin
SPDE and portfolio choice (joint work with M. Musiela) Princeton University November 2007 Thaleia Zariphopoulou The University of Texas at Austin 1 Performance measurement of investment strategies 2 Market
More informationWeak Reflection Principle and Static Hedging of Barrier Options
Weak Reflection Principle and Static Hedging of Barrier Options Sergey Nadtochiy Department of Mathematics University of Michigan Apr 2013 Fields Quantitative Finance Seminar Fields Institute, Toronto
More informationRemarks on American Options for Jump Diffusions
Remarks on American Options for Jump Diffusions Hao Xing University of Michigan joint work with Erhan Bayraktar, University of Michigan Department of Systems Engineering & Engineering Management Chinese
More informationAMH4 - ADVANCED OPTION PRICING. Contents
AMH4 - ADVANCED OPTION PRICING ANDREW TULLOCH Contents 1. Theory of Option Pricing 2 2. Black-Scholes PDE Method 4 3. Martingale method 4 4. Monte Carlo methods 5 4.1. Method of antithetic variances 5
More informationAll Investors are Risk-averse Expected Utility Maximizers
All Investors are Risk-averse Expected Utility Maximizers Carole Bernard (UW), Jit Seng Chen (GGY) and Steven Vanduffel (Vrije Universiteit Brussel) AFFI, Lyon, May 2013. Carole Bernard All Investors are
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 informationForward Dynamic Utility
Forward Dynamic Utility El Karoui Nicole & M RAD Mohamed UnivParis VI / École Polytechnique,CMAP elkaroui@cmapx.polytechnique.fr with the financial support of the "Fondation du Risque" and the Fédération
More informationGeneralized Multi-Factor Commodity Spot Price Modeling through Dynamic Cournot Resource Extraction Models
Generalized Multi-Factor Commodity Spot Price Modeling through Dynamic Cournot Resource Extraction Models Bilkan Erkmen (joint work with Michael Coulon) Workshop on Stochastic Games, Equilibrium, and Applications
More informationIlliquidity, Credit risk and Merton s model
Illiquidity, Credit risk and Merton s model (joint work with J. Dong and L. Korobenko) A. Deniz Sezer University of Calgary April 28, 2016 Merton s model of corporate debt A corporate bond is a contingent
More informationCredit Risk using Time Changed Brownian Motions
Credit Risk using Time Changed Brownian Motions Tom Hurd Mathematics and Statistics McMaster University Joint work with Alexey Kuznetsov (New Brunswick) and Zhuowei Zhou (Mac) 2nd Princeton Credit Conference
More informationUnified Credit-Equity Modeling
Unified Credit-Equity Modeling Rafael Mendoza-Arriaga Based on joint research with: Vadim Linetsky and Peter Carr The University of Texas at Austin McCombs School of Business (IROM) Recent Advancements
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 informationOptimal Investment for Worst-Case Crash Scenarios
Optimal Investment for Worst-Case Crash Scenarios A Martingale Approach Frank Thomas Seifried Department of Mathematics, University of Kaiserslautern June 23, 2010 (Bachelier 2010) Worst-Case Portfolio
More informationOptimal Investment with Deferred Capital Gains Taxes
Optimal Investment with Deferred Capital Gains Taxes A Simple Martingale Method Approach Frank Thomas Seifried University of Kaiserslautern March 20, 2009 F. Seifried (Kaiserslautern) Deferred Capital
More informationFourier Methods for Pricing Early-Exercise Options Under Lévy Dynamics
Fourier Methods for Pricing Early-Exercise Options Under Lévy Dynamics by Tolulope Rhoda Fadina Thesis presented in partial fulfilment of the requirements for the degree of Master of Science in Mathematics
More informationSome mathematical results for Black-Scholes-type equations for financial derivatives
1 Some mathematical results for Black-Scholes-type equations for financial derivatives Ansgar Jüngel Vienna University of Technology www.jungel.at.vu (joint work with Bertram Düring, University of Mainz)
More informationMSc Financial Engineering CHRISTMAS ASSIGNMENT: MERTON S JUMP-DIFFUSION MODEL. To be handed in by monday January 28, 2013
MSc Financial Engineering 2012-13 CHRISTMAS ASSIGNMENT: MERTON S JUMP-DIFFUSION MODEL To be handed in by monday January 28, 2013 Department EMS, Birkbeck Introduction The assignment consists of Reading
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 informationSpot and forward dynamic utilities. and their associated pricing systems. Thaleia Zariphopoulou. UT, Austin
Spot and forward dynamic utilities and their associated pricing systems Thaleia Zariphopoulou UT, Austin 1 Joint work with Marek Musiela (BNP Paribas, London) References A valuation algorithm for indifference
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 informationAnalytical Option Pricing under an Asymmetrically Displaced Double Gamma Jump-Diffusion Model
Analytical Option Pricing under an Asymmetrically Displaced Double Gamma Jump-Diffusion Model Advances in Computational Economics and Finance Univerity of Zürich, Switzerland Matthias Thul 1 Ally Quan
More informationIn chapter 5, we approximated the Black-Scholes model
Chapter 7 The Black-Scholes Equation In chapter 5, we approximated the Black-Scholes model ds t /S t = µ dt + σ dx t 7.1) with a suitable Binomial model and were able to derive a pricing formula for option
More informationLocal Variance Gamma Option Pricing Model
Local Variance Gamma Option Pricing Model Peter Carr at Courant Institute/Morgan Stanley Joint work with Liuren Wu June 11, 2010 Carr (MS/NYU) Local Variance Gamma June 11, 2010 1 / 29 1 Automated Option
More information13.3 A Stochastic Production Planning Model
13.3. A Stochastic Production Planning Model 347 From (13.9), we can formally write (dx t ) = f (dt) + G (dz t ) + fgdz t dt, (13.3) dx t dt = f(dt) + Gdz t dt. (13.33) The exact meaning of these expressions
More informationStock Loan Valuation Under Brownian-Motion Based and Markov Chain Stock Models
Stock Loan Valuation Under Brownian-Motion Based and Markov Chain Stock Models David Prager 1 1 Associate Professor of Mathematics Anderson University (SC) Based on joint work with Professor Qing Zhang,
More informationHedging Under Jump Diffusions with Transaction Costs. Peter Forsyth, Shannon Kennedy, Ken Vetzal University of Waterloo
Hedging Under Jump Diffusions with Transaction Costs Peter Forsyth, Shannon Kennedy, Ken Vetzal University of Waterloo Computational Finance Workshop, Shanghai, July 4, 2008 Overview Overview Single factor
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 informationPricing 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 informationAn Explicit Example of a Shadow Price Process with Stochastic Investment Opportunity Set
An Explicit Example of a Shadow Price Process with Stochastic Investment Opportunity Set Christoph Czichowsky Faculty of Mathematics University of Vienna SIAM FM 12 New Developments in Optimal Portfolio
More informationRegression estimation in continuous time with a view towards pricing Bermudan options
with a view towards pricing Bermudan options Tagung des SFB 649 Ökonomisches Risiko in Motzen 04.-06.06.2009 Financial engineering in times of financial crisis Derivate... süßes Gift für die Spekulanten
More informationOptimal Trade Execution: Mean Variance or Mean Quadratic Variation?
Optimal Trade Execution: Mean Variance or Mean Quadratic Variation? Peter Forsyth 1 S. Tse 2 H. Windcliff 2 S. Kennedy 2 1 Cheriton School of Computer Science University of Waterloo 2 Morgan Stanley New
More informationPath-dependent options and randomization techniques
Path-dependent options and randomization techniques Daniele Marazzina Joint work with G. Fusai and M. Marena Dipartimento di Matematica F. Brioschi Politecnico di Milano daniele.marazzina@polimi.it April
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 informationLecture 1: Lévy processes
Lecture 1: Lévy processes A. E. Kyprianou Department of Mathematical Sciences, University of Bath 1/ 22 Lévy processes 2/ 22 Lévy processes A process X = {X t : t 0} defined on a probability space (Ω,
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 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 informationTime-consistent mean-variance portfolio optimization: a numerical impulse control approach
1 2 3 Time-consistent mean-variance portfolio optimization: a numerical impulse control approach Pieter Van Staden Duy-Minh Dang Peter A. Forsyth 4 5 6 7 8 9 1 11 12 13 14 15 16 17 18 19 2 21 22 23 Abstract
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 informationCalibration Lecture 4: LSV and Model Uncertainty
Calibration Lecture 4: LSV and Model Uncertainty March 2017 Recap: Heston model Recall the Heston stochastic volatility model ds t = rs t dt + Y t S t dw 1 t, dy t = κ(θ Y t ) dt + ξ Y t dw 2 t, where
More informationOn Asymptotic Power Utility-Based Pricing and Hedging
On Asymptotic Power Utility-Based Pricing and Hedging Johannes Muhle-Karbe ETH Zürich Joint work with Jan Kallsen and Richard Vierthauer LUH Kolloquium, 21.11.2013, Hannover Outline Introduction Asymptotic
More informationBIRKBECK (University of London) MSc EXAMINATION FOR INTERNAL STUDENTS MSc FINANCIAL ENGINEERING DEPARTMENT OF ECONOMICS, MATHEMATICS AND STATIS- TICS
BIRKBECK (University of London) MSc EXAMINATION FOR INTERNAL STUDENTS MSc FINANCIAL ENGINEERING DEPARTMENT OF ECONOMICS, MATHEMATICS AND STATIS- TICS PRICING EMMS014S7 Tuesday, May 31 2011, 10:00am-13.15pm
More informationA discretionary stopping problem with applications to the optimal timing of investment decisions.
A discretionary stopping problem with applications to the optimal timing of investment decisions. Timothy Johnson Department of Mathematics King s College London The Strand London WC2R 2LS, UK Tuesday,
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 informationOn Using Shadow Prices in Portfolio optimization with Transaction Costs
On Using Shadow Prices in Portfolio optimization with Transaction Costs Johannes Muhle-Karbe Universität Wien Joint work with Jan Kallsen Universidad de Murcia 12.03.2010 Outline The Merton problem The
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 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 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 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 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 informationRisk Neutral Valuation
copyright 2012 Christian Fries 1 / 51 Risk Neutral Valuation Christian Fries Version 2.2 http://www.christian-fries.de/finmath April 19-20, 2012 copyright 2012 Christian Fries 2 / 51 Outline Notation Differential
More informationOn the valuation of barrier and American options in local volatility models with jumps
On the valuation of barrier and American options in local volatility models with jumps A thesis presented for the degree of Doctor of Philosophy of Imperial College London by Bjorn Eriksson Department
More informationOptimal Securitization via Impulse Control
Optimal Securitization via Impulse Control Rüdiger Frey (joint work with Roland C. Seydel) Mathematisches Institut Universität Leipzig and MPI MIS Leipzig Bachelier Finance Society, June 21 (1) Optimal
More informationOptimal Execution: II. Trade Optimal Execution
Optimal Execution: II. Trade Optimal Execution René Carmona Bendheim Center for Finance Department of Operations Research & Financial Engineering Princeton University Purdue June 21, 212 Optimal Execution
More informationMean-Quadratic Variation Portfolio Optimization: A desirable alternative to Time-consistent Mean-Variance Optimization?
1 2 3 4 Mean-Quadratic Variation Portfolio Optimization: A desirable alternative to Time-consistent Mean-Variance Optimization? Pieter M. van Staden Duy-Minh Dang Peter A. Forsyth October 24, 218 5 6 7
More informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Simulating Stochastic Differential Equations Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
More informationReal Options and Game Theory in Incomplete Markets
Real Options and Game Theory in Incomplete Markets M. Grasselli Mathematics and Statistics McMaster University IMPA - June 28, 2006 Strategic Decision Making Suppose we want to assign monetary values to
More informationA High-order Front-tracking Finite Difference Method for Pricing American Options under Jump-Diffusion Models
A High-order Front-tracking Finite Difference Method for Pricing American Options under Jump-Diffusion Models Jari Toivanen Abstract A free boundary formulation is considered for the price of American
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 informationChapter 15: Jump Processes and Incomplete Markets. 1 Jumps as One Explanation of Incomplete Markets
Chapter 5: Jump Processes and Incomplete Markets Jumps as One Explanation of Incomplete Markets It is easy to argue that Brownian motion paths cannot model actual stock price movements properly in reality,
More information2.1 Mean-variance Analysis: Single-period Model
Chapter Portfolio Selection The theory of option pricing is a theory of deterministic returns: we hedge our option with the underlying to eliminate risk, and our resulting risk-free portfolio then earns
More informationHedging Credit Derivatives in Intensity Based Models
Hedging Credit Derivatives in Intensity Based Models PETER CARR Head of Quantitative Financial Research, Bloomberg LP, New York Director of the Masters Program in Math Finance, Courant Institute, NYU Stanford
More informationPricing in markets modeled by general processes with independent increments
Pricing in markets modeled by general processes with independent increments Tom Hurd Financial Mathematics at McMaster www.phimac.org Thanks to Tahir Choulli and Shui Feng Financial Mathematics Seminar
More informationOptimal Asset Allocation with Stochastic Interest Rates in Regime-switching Models
Optimal Asset Allocation with Stochastic Interest Rates in Regime-switching Models Ruihua Liu Department of Mathematics University of Dayton, Ohio Joint Work With Cheng Ye and Dan Ren To appear in International
More informationEconomathematics. Problem Sheet 1. Zbigniew Palmowski. Ws 2 dw s = 1 t
Economathematics Problem Sheet 1 Zbigniew Palmowski 1. Calculate Ee X where X is a gaussian random variable with mean µ and volatility σ >.. Verify that where W is a Wiener process. Ws dw s = 1 3 W t 3
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 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 informationEfficient Static Replication of European Options under Exponential Lévy Models
CIRJE-F-539 Efficient Static Replication of European Options under Exponential Lévy Models Akihiko Takahashi University of Tokyo Akira Yamazaki Mizuho-DL Financial Technology Co., Ltd. January 8 CIRJE
More informationOption Pricing for a Stochastic-Volatility Jump-Diffusion Model with Log-Uniform Jump-Amplitudes
Option Pricing for a Stochastic-Volatility Jump-Diffusion Model with Log-Uniform Jump-Amplitudes Floyd B. Hanson and Guoqing Yan Department of Mathematics, Statistics, and Computer Science University of
More informationPortability, salary and asset price risk: a continuous-time expected utility comparison of DB and DC pension plans
Portability, salary and asset price risk: a continuous-time expected utility comparison of DB and DC pension plans An Chen University of Ulm joint with Filip Uzelac (University of Bonn) Seminar at SWUFE,
More informationAsset Pricing Models with Underlying Time-varying Lévy Processes
Asset Pricing Models with Underlying Time-varying Lévy Processes Stochastics & Computational Finance 2015 Xuecan CUI Jang SCHILTZ University of Luxembourg July 9, 2015 Xuecan CUI, Jang SCHILTZ University
More informationPolynomial processes in stochastic portofolio theory
Polynomial processes in stochastic portofolio theory Christa Cuchiero University of Vienna 9 th Bachelier World Congress July 15, 2016 Christa Cuchiero (University of Vienna) Polynomial processes in SPT
More informationComputational Methods in Finance
Chapman & Hall/CRC FINANCIAL MATHEMATICS SERIES Computational Methods in Finance AM Hirsa Ltfi) CRC Press VV^ J Taylor & Francis Group Boca Raton London New York CRC Press is an imprint of the Taylor &
More informationLikelihood Estimation of Jump-Diffusions
Likelihood Estimation of Jump-Diffusions Extensions from Diffusions to Jump-Diffusions, Implementation with Automatic Differentiation, and Applications Berent Ånund Strømnes Lunde DEPARTMENT OF MATHEMATICS
More informationSample Path Large Deviations and Optimal Importance Sampling for Stochastic Volatility Models
Sample Path Large Deviations and Optimal Importance Sampling for Stochastic Volatility Models Scott Robertson Carnegie Mellon University scottrob@andrew.cmu.edu http://www.math.cmu.edu/users/scottrob June
More informationThe investment game in incomplete markets
The investment game in incomplete markets M. R. Grasselli Mathematics and Statistics McMaster University Pisa, May 23, 2008 Strategic decision making We are interested in assigning monetary values to strategic
More informationStochastic Partial Differential Equations and Portfolio Choice. Crete, May Thaleia Zariphopoulou
Stochastic Partial Differential Equations and Portfolio Choice Crete, May 2011 Thaleia Zariphopoulou Oxford-Man Institute and Mathematical Institute University of Oxford and Mathematics and IROM, The University
More informationShort-time-to-expiry expansion for a digital European put option under the CEV model. November 1, 2017
Short-time-to-expiry expansion for a digital European put option under the CEV model November 1, 2017 Abstract In this paper I present a short-time-to-expiry asymptotic series expansion for a digital European
More informationStochastic modelling of electricity markets Pricing Forwards and Swaps
Stochastic modelling of electricity markets Pricing Forwards and Swaps Jhonny Gonzalez School of Mathematics The University of Manchester Magical books project August 23, 2012 Clip for this slide Pricing
More informationSTOCHASTIC VOLATILITY AND OPTION PRICING
STOCHASTIC VOLATILITY AND OPTION PRICING Daniel Dufresne Centre for Actuarial Studies University of Melbourne November 29 (To appear in Risks and Rewards, the Society of Actuaries Investment Section Newsletter)
More informationA new PDE approach for pricing arithmetic average Asian options
A new PDE approach for pricing arithmetic average Asian options Jan Večeř Department of Mathematical Sciences, Carnegie Mellon University, Pittsburgh, PA 15213. Email: vecer@andrew.cmu.edu. May 15, 21
More informationThe investment game in incomplete markets.
The investment game in incomplete markets. M. R. Grasselli Mathematics and Statistics McMaster University RIO 27 Buzios, October 24, 27 Successes and imitations of Real Options Real options accurately
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 informationDynamic Relative Valuation
Dynamic Relative Valuation Liuren Wu, Baruch College Joint work with Peter Carr from Morgan Stanley October 15, 2013 Liuren Wu (Baruch) Dynamic Relative Valuation 10/15/2013 1 / 20 The standard approach
More information