Hilbert transform approach for pricing Bermudan options in Lévy models
|
|
- Curtis Barker
- 6 years ago
- Views:
Transcription
1 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 and cubature methods in finance and econometrics University of Leicester, 6/19/2009
2 Lévy models Introduction Lévy models Hilbert transform methods Better fit to empirical data, explain volatility smiles Finite activity jump diffusion models Merton 76 (normal jump diffusion), Kou 02 (double exponential jump diffusion) Infinite activity pure jump Lévy models Madan et al 90, 91, 98 (variance gamma), Barndorff-Nielsen 98 (normal inverse Gaussian), Carr, Geman, Madan & Yor 02 (CGMY), Eberlein, Keller & Prause 98 (Generalized hyperbolic) Books on Lévy processes Boyarchenko & Levendorskii 02, Cont & Tankov 04, Schoutens 03, Bertoin 96, Sato 99, Applebaum 04, Kyprianou 06
3 Lévy-Khintchine theorem Lévy models Hilbert transform methods Lévy process: independent and stationary increments Analytical tractability: characteristic functions available via Lévy-Khintchine Theorem φ t (ξ) = E[e iξxt ] = e tψ(ξ) where ψ is the characteristic exponent ψ(ξ) = 1 2 σ2 ξ 2 iµξ + (1 e iξx + iξx1 { x 1} )Π(dx) R
4 Fourier transform method Lévy models Hilbert transform methods Fourier transform method for European options: Carr & Madan 99 c(k) = 1 2π e αk e iξk Φ(ξ)dξ R where Φ(ξ) = e rt φ T (ξ i(α + 1)) (iξ + α)(iξ + α + 1) Extended in Lee 04: more general payoffs, exponentially decaying errors Bermudan options: vanilla, discrete barrier/lookback options, Hilbert transform method
5 Hilbert transform Introduction Lévy models Hilbert transform methods Hilbert transform of f L p (R), 1 p < Hf (x) = 1 π p.v. For any f L 1 (R) with ˆf L 1 (R) f (y) x y dy F(sgn f )(ξ) = ihˆf (ξ) F(1 (l, ) f )(ξ) = 1 2ˆf (ξ) + i 2 eiξl H(e iηl ˆf (η))(ξ)
6 Option pricing Introduction Lévy models Hilbert transform methods European vanilla options: E[(K S T ) + ] = KE[1 {ST <K}] E[S T 1 {ST <K}] P(X x) = 1 (,x) (y)p(y)dy R European style discrete barrier options: Feng & Linetsky 08 f j (x) = 1 (l,u) E tj,x[f j+1 (X tj+1 )] European style discrete lookbacks: Feng & Linetsky 09 M j X j = max(m j 1, X j ) X j = max(0, M j 1 X j 1 (X j X j 1 ))
7 Discrete Hilbert transform Lévy models Hilbert transform methods Discrete Hilbert transform with step size h > 0 H h f (x) = m= f (mh) 1 cos[π(x mh)/h], x R π(x mh)/h For f analytic in a horizontal strip {z C : I(z) < d} Hf H h f L (R) Ce πd/h πd(1 e πd/h ) Related to Whittaker cardinal series (sinc) expansion
8 Sinc expansion of analytic functions Lévy models Hilbert transform methods Whittaker cardinal series (sinc expansion) c(f, h)(x) = sin(π(x kh)/h) f (kh) π(x kh)/h For entire functions of exponential type π/h, sinc expansion is exact: c(f, h) = f For functions analytic in a strip {z C : I(z) < d}, Stenger 93 Ce πd/h f c(f, h) L πd(1 e 2πd/h )
9 Trapezoidal rule Introduction Lévy models Hilbert transform methods Take Hilbert transform on c(f, h), obtain discrete Hilbert transform Trapezoidal rule very accurate for f analytic in a strip {z C : I(z) < d} f (x)dx f (kh)h Ce 2πd/h 1 e 2πd/h R m= We use trapezoidal rule to compute Fourier inverse integral
10 Bermudan style vanilla options Backward induction in state space Discrete approximation Bermudan put: payoff G(S) = (K S) +, discrete monitoring T = {t 0, t 1,, t N } = {0,, 2,, N = T } Exponential Lévy model: X t a Lévy process in (Ω, F, F, P), start from ln(s 0 /K), equivalent martingale measure P is given Optimal stopping S t = Ke Xt V 0 (S 0 ) = sup E 0 [e rτ G(S τ )] τ
11 Backward induction in state space Backward induction in state space Discrete approximation Variable change x = ln(s/k), g(x) = G(Ke x ) = K(1 e x ) +, f 0 (x) = V 0 (Ke x ) Backward induction f N (x) = g(x) ( ) f j (x) = max g(x), e r E j,x [f j+1 (X (j+1) )], 0 j < N V 0 (S 0 ) = f 0 (ln(s 0 /K))
12 Direct implementation Backward induction in state space Discrete approximation Knowing the transition density p ( ) E j,x [f j+1 (X (j+1) )] = f j+1 (y)p (y x)dy When discretized, the convolution becomes a Toeplitz matrix vector multiplication, can be implemented using FFT (Eydeland 94 ) Density may not be available, polynomial convergence Double exponential fast Gauss transform (Broadie & Yamamoto 05 ), limited to the Black-Schoels-Merton model and Merton s model For general Lévy models: COS method Fang & Oosterlee 08 R
13 Dampening for integrability Backward induction in state space Discrete approximation For α > 0 (for puts), define f j α(x) = e αx f j (x), g α (x) = e αx g(x) Backward induction fα N (x) = g α (x) ( ) fα(x) j = max g α (x), e αx r E j,x [e αx (j+1) fα j+1 (X (j+1) )], V 0 (S 0 ) = e α ln(s 0/K) f 0 α (ln(s 0 /K))
14 Esscher transform Introduction Backward induction in state space Discrete approximation Define new measure P α via dp α dp F t = e αxt φ t (iα) Esscher transformed Lévy process is still a Lévy process φ α t (ξ) = φ t(ξ + iα) φ t (iα) Applying property E α s [Y t ] = E s [Y t Z t /Z s ] E j,x [e αx (j+1) fα j+1 (X (j+1) )] = e αx ψ(iα) E α j+1 j,x [fα (X (j+1) )]
15 Dampened backward induction Backward induction in state space Discrete approximation Dampened backward induction For 0 j < N f j α(x) = max f N α (x) = g α (x) ( ) g α (x), e (r+ψ(iα)) E α j+1 j,x [fα (X (j+1) )] = g α (x)1 (,x j )(x) +e (r+ψ(iα)) E α j,x [ f j+1 α (X (j+1) ) ] 1 [x j, )(x) where x j is the early exercise boundary at time j
16 Backward induction in state space Discrete approximation By convolution theorem, Fourier transform of E α [ j,x f j+1 α (X (j+1) ) ] = fα j+1 (y)p α (y x)dy is a product ˆf α j+1 (ξ)φ α ( ξ) R ˆf N α (ξ) = ĝ α (ξ) ( 1 ˆf α(ξ) j = F(g α 1 (,x j ))(ξ) + e (r+ψ(iα)) α 2ˆf j+1 (ξ)φ α ( ξ) + i ( ) ) 2 eiξx j H e iηx j ˆf α j+1 (η)φ α ( η) (ξ)
17 Early exercise boundary Backward induction in state space Discrete approximation Early exercise boundary solves g α (x) = e (r+ψ(iα)) E α j+1 j,x [fα (X (j+1) )] Using Fourier inverse representation g α (x) = 1 2π e (r+ψ(iα)) j+1 e iξxˆf α (ξ)φ α (ξ)dξ xn = K. To solve for x j, use Newton-Raphson, with starting point xj+1 R
18 Algorithm summarized Backward induction in state space Discrete approximation Start with Fourier transform of dampened payoff ˆf α N = ĝ α j+1 At time j, with ˆf α, compute early exercise boundary xj using Newton Raphson Compute ˆf α j from ˆf α j+1 and xj (Hilbert transform) With ˆf α 1, option value at time 0 ( fα 0 1 (x) = max g α (x), 2π e (r+ψ(iα)) e iξxˆf ) α 1 (ξ)φ α (ξ)dξ R
19 Discrete approximation Backward induction in state space Discrete approximation Need to repeatedly evaluate a Fourier inverse integral and Hϕ(ξ) Trapezoidal rule for Fourier inverse integral, truncate infinite series with truncation level M, computational cost O(M) Replace Hϕ by discrete Hilbert transform, truncate resulting infinite series Hϕ(ξ) M m= M 1 cos[π(ξ mh)/h] ϕ(mh) π(ξ mh)/h
20 Error estimate Introduction Backward induction in state space Discrete approximation Discretization error O(exp( πd/h) With φ t (ξ) exp( ct ξ ν ), truncation error is essentially Select h = h(m) according to O(exp( c(mh) ν )) h(m) = Total error: O(exp( CM ν 1+ν )) ( ) 1 πd 1+ν M ν 1+ν c
21 Toeplitz matrix vector multiplication Backward induction in state space Discrete approximation Evaluate Hψ(ξ) M m= M for ξ = Mh,, Mh 1 cos[π(ξ mh)/h] ψ(mh) π(ξ mh)/h Correspond to Toeplitz matrix vector multiplication FFT based method for such multiplications: O(M log(m)) Total computational cost of the method: O(NM log(m))
22 NIG model Introduction Backward induction in state space Discrete approximation Pricing Bermudan put option in NIG NIG process: time changed Brownian motion, with an inverse Gaussian process as the stochastic time Characteristic exponent ψ(ξ) = iµξ + δ NIG ( αnig 2 (β NIG + iξ) 2 αnig 2 β2 NIG ) φ t (ξ) has exponential tails with ν = 1, error estimate in M: O(e C M )
23 Bermudan put in the NIG model Backward induction in state space Discrete approximation Figure: T = 1, N = 252, S 0 = 100, K = 110, r = 0.1, q = 0, α NIG = 25, β NIG = 5, δ NIG = 0.5, implemented in Matlab using a Laptop with CPU 2GHz, RAM 1G, last node takes 2.9s
24 Bermudan barrier options Barrier options Lookback options Bermudan down-and-out put with lower barrier L, l = ln(l/k) Backward induction f N (x) = g(x)1 (l, ) (x) For 0 j < N, f j (x) = 1 (l, ) (x) max(g(x), e r E j,x [f j+1 (X (j+1) )]) = g(x)1 (l, ) (x)1 (,x j )(x) +e r E j,x [f j+1 (X (j+1) )]1 (l, ) (x)1 (x j,,)(x)
25 CGMY model Introduction Barrier options Lookback options Pricing Bermudan down-and-out put in CGMY Levy density Characteristic exponent π(x) = CeGx x 1+Y 1 {x<0} + Ce Mx x 1+Y 1 {x>0} ψ(ξ) = iµξ CΓ( Y )((M iξ) Y M Y + (G + iξ) Y G Y ) Characteristic function has exponential tails with ν = Y, error estimate in M: O(exp( CM Y /(1+Y ) ))
26 Barrier options Lookback options Bermudan DOP in the CGMY model Figure: T = 1, N = 252, S 0 = 100, K = 110, r = 0.1, q = 0, C = 3, G = 22, M = 28, Y = 0.8, implemented in Matlab using a Laptop with CPU 2GHz, RAM 1G, last node takes 2.0s
27 Barrier options Lookback options Bermudan floating strike lookbacks Standard backward induction involves two state variables: asset price, maximum asset price Can be reduced to one state variable, maximum asset price/asset price f j (y) = max(e y 1, e q E j,y [f j+1 (e Y (j+1) )] Double exponential fast Gauss transform method (Yamamoto 05 ), limited to BSM and Merton s models
28 Kou s model Introduction Barrier options Lookback options Pricing Bermudan floating strike lookback put Log asset price follows N t X t = µt+σb t + Z n, Z n pη 1 e η1x 1 {x>0} +(1 p)η 2 e η2x 1 {x<0} n=1 Characteristic exponent ψ(ξ) = iµξ σ2 ξ 2 + λ ( 1 pη 1 η 1 iξ (1 p)η ) 2 η 2 + iξ Characteristic function has exponential tails with ν = 2, error estimate O(exp( CM 2/3 ))
29 Barrier options Lookback options Bermudan floating strike lookback put in Kou s model Figure: T = 1, N = 252, S 0 = 100, r = 0.1, q = 0, σ = 0.25, λ = 3, p = 0.3, η 1 = 17, η 2 = 12, implemented in Matlab using a Laptop with CPU 2GHz, RAM 1G, last node takes 1.1s
30 Summary Introduction Barrier options Lookback options Hilbert transform method for pricing Bermudan style options in Lévy process models Very accurate with exponentially decay errors Fast with computational cost O(NM log(m)) European vanilla, barrier, lookback, defaultable bonds,, barrier, floating strike lookback, monte carlo simulation etc.
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 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 informationTHE HILBERT TRANSFORM AND ITS APPLICATIONS IN COMPUTATIONAL FINANCE XIONG LIN DISSERTATION
c 2010 Xiong Lin THE HILBET TANSFOM AND ITS APPLICATIONS IN COMPUTATIONAL FINANCE BY XIONG LIN DISSETATION Submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy in
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 informationUsing Lévy Processes to Model Return Innovations
Using Lévy Processes to Model Return Innovations Liuren Wu Zicklin School of Business, Baruch College Option Pricing Liuren Wu (Baruch) Lévy Processes Option Pricing 1 / 32 Outline 1 Lévy processes 2 Lévy
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 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 informationZisheng Chen Liming Feng. Xiong Lin. University of Illinois at Urbana-Champaign Department of Industrial and Enterprise Systems Engineering
Proceedings of the 2011 Winter Simulation Conference S. Jain, R. R. Creasey, J. Himmelspach, K. P. White, and M. Fu, eds. INVERSE TRANSFORM METHOD FOR SIMULATING LEVY PROCESSES AND DISCRETE ASIAN OPTIONS
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 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 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 informationTHE WIENER-HOPF FACTORIZATION AS A GENERAL METHOD FOR VALUATION OF AMERICAN AND BARRIER OPTIONS
THE WIENER-HOPF FACTORIZATION AS A GENERAL METHOD FOR VALUATION OF AMERICAN AND BARRIER OPTIONS SERGEI LEVENDORSKIǏ Department of Economics, The University of Texas at Austin, 1 University Station C31,
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 informationFractional differential equations for modelling financial processes with jumps
Hong Kong Baptist University HKBU Institutional Repository Open Access Theses and Dissertations Electronic Theses and Dissertations 9-12-2016 Fractional differential equations for modelling financial processes
More informationAdvantages of the Laplace transform approach in pricing first touch digital options in Lévy-driven models
Advantages of the Laplace transform approach in pricing first touch digital options in Lévy-driven models Oleg Kudryavtsev Department of Informatics Russian Customs Academy Rostov Branch Budennovskiy 20
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 informationChapter 1. Introduction and Preliminaries. 1.1 Motivation. The American put option problem
Chapter 1 Introduction and Preliminaries 1.1 Motivation The American put option problem The valuation of contingent claims has been a widely known topic in the theory of modern finance. Typical claims
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 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 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 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 informationEfficient valuation of exotic derivatives in Lévy models
Efficient valuation of exotic derivatives in models Ernst Eberlein and Antonis Papapantoleon Department of Mathematical Stochastics and Center for Data Analysis and Modeling (FDM) University of Freiburg
More informationPower Style Contracts Under Asymmetric Lévy Processes
MPRA Munich Personal RePEc Archive Power Style Contracts Under Asymmetric Lévy Processes José Fajardo FGV/EBAPE 31 May 2016 Online at https://mpra.ub.uni-muenchen.de/71813/ MPRA Paper No. 71813, posted
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 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 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 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 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 informationMachine Learning for Quantitative Finance
Machine Learning for Quantitative Finance Fast derivative pricing Sofie Reyners Joint work with Jan De Spiegeleer, Dilip Madan and Wim Schoutens Derivative pricing is time-consuming... Vanilla option pricing
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 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 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 informationQuadratic hedging in affine stochastic volatility models
Quadratic hedging in affine stochastic volatility models Jan Kallsen TU München Pittsburgh, February 20, 2006 (based on joint work with F. Hubalek, L. Krawczyk, A. Pauwels) 1 Hedging problem S t = S 0
More informationFactor Models for Option Pricing
Factor Models for Option Pricing Peter Carr Banc of America Securities 9 West 57th Street, 40th floor New York, NY 10019 Tel: 212-583-8529 email: pcarr@bofasecurities.com Dilip B. Madan Robert H. Smith
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 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 informationJump-type Lévy processes
Jump-type Lévy processes Ernst Eberlein Department of Mathematical Stochastics, University of Freiburg, Eckerstr. 1, 7914 Freiburg, Germany, eberlein@stochastik.uni-freiburg.de 1 Probabilistic structure
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 informationPricing exotic options under regime switching: A Fourier transform method
Pricing exotic options under regime switching: A Fourier transform method Peter Hieber Lehrstuhl für Finanzmathematik M13, Technische Universität München, Parkring 11, 85748 Garching-Hochbrück, Germany,
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 informationFast Hilbert transform algorithms for pricing discrete timer options under stochastic volatility models
Fast Hilbert transform algorithms for pricing discrete timer options under stochastic volatility models Pingping Zeng, Yue Kuen Kwok 1, and Wendong Zheng Department of Mathematics, Hong Kong 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 informationSato Processes in Finance
Sato Processes in Finance Dilip B. Madan Robert H. Smith School of Business Slovenia Summer School August 22-25 2011 Lbuljana, Slovenia OUTLINE 1. The impossibility of Lévy processes for the surface of
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 informationValuation of Volatility Derivatives. Jim Gatheral Global Derivatives & Risk Management 2005 Paris May 24, 2005
Valuation of Volatility Derivatives Jim Gatheral Global Derivatives & Risk Management 005 Paris May 4, 005 he opinions expressed in this presentation are those of the author alone, and do not necessarily
More informationA GENERAL FORMULA FOR OPTION PRICES IN A STOCHASTIC VOLATILITY MODEL. Stephen Chin and Daniel Dufresne. Centre for Actuarial Studies
A GENERAL FORMULA FOR OPTION PRICES IN A STOCHASTIC VOLATILITY MODEL Stephen Chin and Daniel Dufresne Centre for Actuarial Studies University of Melbourne Paper: http://mercury.ecom.unimelb.edu.au/site/actwww/wps2009/no181.pdf
More informationNumerical Algorithms for Pricing Discrete Variance and Volatility Derivatives under Time-changed Lévy Processes
Numerical Algorithms for Pricing Discrete Variance and Volatility Derivatives under Time-changed Lévy Processes WENDONG ZHENG, CHI HUNG YUEN & YUE KUEN KWOK 1 Department of Mathematics, Hong Kong University
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 informationASYMMETRICALLY TEMPERED STABLE DISTRIBUTIONS WITH APPLICATIONS TO FINANCE
PROBABILITY AND MATHEMATICAL STATISTICS Vol. 0, Fasc. 0 (0000), pp. 000 000 doi: ASYMMETRICALLY TEMPERED STABLE DISTRIBUTIONS WITH APPLICATIONS TO FINANCE A. A R E F I (ALLAMEH TABATABA I UNIVERSITY) AND
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 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 informationHedging of swaptions in a Lévy driven Heath-Jarrow-Morton framework
Hedging of swaptions in a Lévy driven Heath-Jarrow-Morton framework Kathrin Glau, Nele Vandaele, Michèle Vanmaele Bachelier Finance Society World Congress 2010 June 22-26, 2010 Nele Vandaele Hedging of
More informationOptimal Stopping for American Type Options
Optimal Stopping for Department of Mathematics Stockholm University Sweden E-mail: silvestrov@math.su.se ISI 2011, Dublin, 21-26 August 2011 Outline of communication Multivariate Modulated Markov price
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 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 informationPricing Foreign Equity Option with time-changed Lévy Process
Pricing Foreign Equity Option with time-changed Lévy Process Abstract. In this paper we propose a general foreign equity option pricing framework that unifies the vast foreign equity option pricing literature
More informationThe Heston model with stochastic interest rates and pricing options with Fourier-cosine expansion
The Heston model with stochastic interest rates and pricing options with Fourier-cosine expansion Kees Oosterlee 1,2 1 CWI, Center for Mathematics and Computer Science, Amsterdam, 2 Delft University of
More informationFair Valuation of Insurance Contracts under Lévy Process Specifications Preliminary Version
Fair Valuation of Insurance Contracts under Lévy Process Specifications Preliminary Version Rüdiger Kiesel, Thomas Liebmann, Stefan Kassberger University of Ulm and LSE June 8, 2005 Abstract The valuation
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 informationFast and accurate pricing of barrier options under Lévy processes
INSTITUT NATIONAL DE RECHERCHE EN INFORMATIQUE ET EN AUTOMATIQUE Fast and accurate pricing of barrier options under Lévy processes Oleg Kudryavtsev Sergei Levendorskiǐ N 6670 Octobre 2008 Thème NUM apport
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 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 informationLévy Processes. Antonis Papapantoleon. TU Berlin. Computational Methods in Finance MSc course, NTUA, Winter semester 2011/2012
Lévy Processes Antonis Papapantoleon TU Berlin Computational Methods in Finance MSc course, NTUA, Winter semester 2011/2012 Antonis Papapantoleon (TU Berlin) Lévy processes 1 / 41 Overview of the course
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 informationJOHN CROSBY, NOLWENN LE SAUX, AND ALEKSANDAR MIJATOVIĆ
APPROXIMATING LÉVY PROCESSES WITH A VIEW TO OPTION PRICING JOHN CROSBY, NOLWENN LE SAUX, AND ALEKSANDAR MIJATOVIĆ Abstract. We examine how to approximate a Lévy process by a hyperexponential jump-diffusion
More informationAN INTRODUCTION TO LÉVY PROCESSES WITH APPLICATIONS IN FINANCE
AN INTODUCTION TO LÉVY POCESSES WITH APPLICATIONS IN FINANCE ANTONIS PAPAPANTOLEON Abstract. These notes aim at introducing Lévy processes in an informal and intuitive way. Several important results about
More informationMgr. Jakub Petrásek 1. May 4, 2009
Dissertation Report - First Steps Petrásek 1 2 1 Department of Probability and Mathematical Statistics, Charles University email:petrasek@karlin.mff.cuni.cz 2 RSJ Invest a.s., Department of Probability
More informationPricing and Hedging of European Plain Vanilla Options under Jump Uncertainty
Pricing and Hedging of European Plain Vanilla Options under Jump Uncertainty by Olaf Menkens School of Mathematical Sciences Dublin City University (DCU) Financial Engineering Workshop Cass Business School,
More informationIMPLEMENTING THE SPECTRAL CALIBRATION OF EXPONENTIAL LÉVY MODELS
IMPLEMENTING THE SPECTRAL CALIBRATION OF EXPONENTIAL LÉVY MODELS DENIS BELOMESTNY AND MARKUS REISS 1. Introduction The aim of this report is to describe more precisely how the spectral calibration method
More informationSmall-time asymptotics of stopped Lévy bridges and simulation schemes with controlled bias
Small-time asymptotics of stopped Lévy bridges and simulation schemes with controlled bias José E. Figueroa-López 1 1 Department of Statistics Purdue University Computational Finance Seminar Purdue University
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 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 informationAnalysis of Fourier Transform Valuation Formulas and Applications
Analysis of Fourier Transform Formulas and Applications Ernst Eberlein Freiburg Institute for Advanced Studies (FRIAS) and Center for Data Analysis and Modeling (FDM) University of Freiburg (joint work
More informationABSTRACT. Professor Dilip B. Madan Department of Finance
ABSTRACT Title of dissertation: The Multivariate Variance Gamma Process and Its Applications in Multi-asset Option Pricing Jun Wang, Doctor of Philosophy, 2009 Dissertation directed by: Professor Dilip
More informationSkewness Premium with Lévy Processes
Skewness Premium with Lévy Processes José Fajardo and Ernesto Mordecki September 15, 2006 Abstract We study the skewness premium (SK) introduced by Bates (1991) in a general context using Lévy Processes.
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 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 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 informationShort-time asymptotics for ATM option prices under tempered stable processes
Short-time asymptotics for ATM option prices under tempered stable processes José E. Figueroa-López 1 1 Department of Statistics Purdue University Probability Seminar Purdue University Oct. 30, 2012 Joint
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 informationFourier, Wavelet and Monte Carlo Methods in Computational Finance
Fourier, Wavelet and Monte Carlo Methods in Computational Finance Kees Oosterlee 1,2 1 CWI, Amsterdam 2 Delft University of Technology, the Netherlands AANMPDE-9-16, 7/7/26 Kees Oosterlee (CWI, TU Delft)
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 informationA Modified Tempered Stable Distribution with Volatility Clustering
A Modified Tempered Stable Distribution with Volatility Clustering Y.S. Kim, S.T. Rachev, D.M. Chung, M.L. Bianchi Abstract We first introduce a new variant of the tempered stable distribution, named the
More informationStochastic volatility modeling in energy markets
Stochastic volatility modeling in energy markets Fred Espen Benth Centre of Mathematics for Applications (CMA) University of Oslo, Norway Joint work with Linda Vos, CMA Energy Finance Seminar, Essen 18
More informationPricing Variance Swaps on Time-Changed Lévy Processes
Pricing Variance Swaps on Time-Changed Lévy Processes ICBI Global Derivatives Volatility and Correlation Summit April 27, 2009 Peter Carr Bloomberg/ NYU Courant pcarr4@bloomberg.com Joint with Roger Lee
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 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 informationMARIANNA MOROZOVA IEIE SB RAS, Novosibirsk, Russia
MARIANNA MOROZOVA IEIE SB RAS, Novosibirsk, Russia 1 clue of ineffectiveness: BS prices are fair only in case of complete markets FORTS is clearly not complete (as log. returns are not Normal) Market prices
More informationThe Use of Importance Sampling to Speed Up Stochastic Volatility Simulations
The Use of Importance Sampling to Speed Up Stochastic Volatility Simulations Stan Stilger June 6, 1 Fouque and Tullie use importance sampling for variance reduction in stochastic volatility simulations.
More informationA Consistent Pricing Model for Index Options and Volatility Derivatives
A Consistent Pricing Model for Index Options and Volatility Derivatives 6th World Congress of the Bachelier Society Thomas Kokholm Finance Research Group Department of Business Studies Aarhus School of
More informationA Simulation Study of Bipower and Thresholded Realized Variations for High-Frequency Data
Washington University in St. Louis Washington University Open Scholarship Arts & Sciences Electronic Theses and Dissertations Arts & Sciences Spring 5-18-2018 A Simulation Study of Bipower and Thresholded
More informationPricing Asian Options in a Semimartingale Model
Pricing Asian Options in a Semimartingale Model Jan Večeř Columbia University, Department of Statistics, New York, NY 127, USA Kyoto University, Kyoto Institute of Economic Research, Financial Engineering
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 informationBasic Stochastic Processes
Basic Stochastic Processes Series Editor Jacques Janssen Basic Stochastic Processes Pierre Devolder Jacques Janssen Raimondo Manca First published 015 in Great Britain and the United States by ISTE Ltd
More informationUniversità degli Studi di Bergamo Facoltà di Economia
Università degli Studi di Bergamo Facoltà di Economia Dipartimento di Matematica, Statistica, Informatica e Applicazioni XX CICLO DOTTORATO DI RICERCA Metodi computazionali per le previsioni e decisioni
More informationSkewness in Lévy Markets
Skewness in Lévy Markets Ernesto Mordecki Universidad de la República, Montevideo, Uruguay Lecture IV. PASI - Guanajuato - June 2010 1 1 Joint work with José Fajardo Barbachan Outline Aim of the talk Understand
More informationAN ANALYTICALLY TRACTABLE UNCERTAIN VOLATILITY MODEL
AN ANALYTICALLY TRACTABLE UNCERTAIN VOLATILITY MODEL FABIO MERCURIO BANCA IMI, MILAN http://www.fabiomercurio.it 1 Stylized facts Traders use the Black-Scholes formula to price plain-vanilla options. An
More informationPricing Derivatives on Two Lévy-driven Stocks
Pricing Derivatives on Two Lévy-driven Stocks José Fajardo and Ernesto Mordecki March 21, 2003 Abstract The aim of this work is to study the pricing problem for derivatives depending on two stocks driven
More information