Two-dimensional COS method
|
|
- Amelia Leonard
- 5 years ago
- Views:
Transcription
1 Two-dimensional COS method Marjon Ruijter Winterschool Lunteren 22 January /29
2 Introduction PhD student since October 2010 Prof.dr.ir. C.W. Oosterlee). CWI national research center for mathematics and computer science. CPB Netherlands Bureau for Economic Policy Analysis. Impact of climate change on investments and policy decisions. Financial mathematics. 2/29
3 Stochastic optimization - 2D Climate-economics problem Two stochastic processes: temperature and capital. Goal: maximize expected utility [ T ] vt, x 1, x 2 )= max E e ρs t) UC s ))ds F t. 1) {a s,c s} t 3/29
4 Stochastic optimization - 2D Climate-economics problem Two stochastic processes: temperature and capital. Goal: maximize expected utility [ T ] vt, x 1, x 2 )= max E e ρs t) UC s ))ds F t. 1) {a s,c s} t Rainbow option pricing problem Two stochastic processes: asset price 1 and asset price 2. Goal: maximize expected profit [ ] vt, x 1, x 2 )= max E e rτ t) gxτ 1, Xτ 2 ) F t. 2) τ [t,t ] 3/29
5 Stochastic optimization - 2D Climate-economics problem Two stochastic processes: temperature and capital. Goal: maximize expected utility [ T ] vt, x 1, x 2 )= max E e ρs t) UC s ))ds F t. 1) {a s,c s} t Rainbow option pricing problem Two stochastic processes: asset price 1 and asset price 2. Goal: maximize expected profit [ ] vt, x 1, x 2 )= max E e rτ t) gxτ 1, Xτ 2 ) F t. 2) τ [t,t ] The math is almost) the same 3/29
6 Table of Contents 1 Rainbow options 2 2D-COS formula 3 European options 4 Bermudan options 5 Heston stochastic volatility model 6 Conclusion 4/29
7 Financial mathematics In financial markets, traders deal in assets and options. The payoff of an option depends on the value of the underlying asset prices). Asset price X t is stochastic. 5/29
8 Financial mathematics In financial markets, traders deal in assets and options. The payoff of an option depends on the value of the underlying asset prices). Asset price X t is stochastic. Payoff call: g call x) =maxx K, 0) 3) Payoff put: g put x) =maxk x, 0) 4) 30 call 30 put gx) gx) x x Figure 1: Payoff call and put option, strike price K = 100 1D). 5/29
9 2 correlated asset prices Two stochastic asset price processes, X 1 t and X 2 t. 6/29
10 2 correlated asset prices Two stochastic asset price processes, Xt 1 and Xt 2. For example, correlated geometric Brownian motions: dxt 1 = μ 1Xt 1 dt + σ 1Xt 1 dw t 1, 5) dxt 2 = μ 2Xt 2 dt + σ 2Xt 2 dw t 2, 6) with dwt 1 dwt 2 = ρdt. 130 Correlated Asset Paths Asset Price Time to Expiry 6/29
11 Payoff rainbow options Basket option: weighted sum or average of different assets, e.g., gx 1, x 2 )=max 1 2 x 1 + x 2 ) K, 0 ). 7) gx 1,x 2 ) x x 1 Figure 2: Basket option. 7/29
12 Payoff rainbow options Basket option: weighted sum or average of different assets, e.g., gx 1, x 2 )=max 1 2 x 1 + x 2 ) K, 0 ). 7) gx 1,x 2 ) 10 5 gx 1,x 2 ) x x x x Figure 2: Basket option. Call on maximum option: Figure 3: Call on max. option. gx 1, x 2 )=maxmaxx 1, x 2 ) K, 0). 8) 7/29
13 European, American, and Bermudan-style European-style: you buy the option now, wait until terminal time T, then the option may be exercised. American-style: may be exercised at any time before the terminal time T. Bermudan-style: fixed exercise dates t m m =1,...,M) atwhich you can exercise the option. Financial mathematics: efficient computation of option price. vt 0, x 0 )=e rδt E [vt, X T )]. 9) 8/29
14 COS method Based on Fourier cosine series expansions. 9/29
15 COS method Based on Fourier cosine series expansions. Pricing financial and real options: European options F. Fang, C.W. Oosterlee, 2008), Bermudan and American options F. Fang, C.W. Oosterlee, 2009), Swing options, which are frequently used in energy markets B. Zhang, C.W. Oosterlee, 2010), Asian-style options B. Zhang, C.W. Oosterlee, 2011), Optimal dike height, M.J. Ruijter, master thesis, 2010),... 9/29
16 Fourier-cosine series expansion of function hx) on[a, b]: hx) = H k cos kπ x a ), x [a, b], 10) with coefficients k=0 H k = 2 b a hy)cos kπ y a ) dy. 11) 10 / 29
17 Fourier-cosine series expansion of function hx) on[a, b]: N 1 ĥx) = H k cos kπ x a ), x [a, b], 10) with coefficients k=0 H k = 2 b a hy)cos kπ y a ) dy. 11) 10 / 29
18 Fourier-cosine series expansion of function hx) on[a, b]: N 1 ĥx) = H k cos kπ x a ), x [a, b], 10) with coefficients k=0 H k = 2 b a hy)cos kπ y a ) dy. 11) N=2 60 hx) x 10 / 29
19 Fourier-cosine series expansion of function hx) on[a, b]: N 1 ĥx) = H k cos kπ x a ), x [a, b], 10) with coefficients k=0 H k = 2 b a hy)cos kπ y a ) dy. 11) N=3 60 hx) x 10 / 29
20 Fourier-cosine series expansion of function hx) on[a, b]: N 1 ĥx) = H k cos kπ x a ), x [a, b], 10) with coefficients k=0 H k = 2 b a hy)cos kπ y a ) dy. 11) N=4 60 hx) x 10 / 29
21 Fourier-cosine series expansion of function hx) on[a, b]: N 1 ĥx) = H k cos kπ x a ), x [a, b], 10) with coefficients k=0 H k = 2 b a hy)cos kπ y a ) dy. 11) N=6 60 hx) x 10 / 29
22 Fourier-cosine series expansion of function hx) on[a, b]: N 1 ĥx) = H k cos kπ x a ), x [a, b], 10) with coefficients k=0 H k = 2 b a hy)cos kπ y a ) dy. 11) N=10 60 hx) x 10 / 29
23 Fourier-cosine series expansion of function hx) on[a, b]: N 1 ĥx) = H k cos kπ x a ), x [a, b], 10) with coefficients k=0 H k = 2 b a hy)cos kπ y a ) dy. 11) N=50 60 hx) x 10 / 29
24 1D-COS formula We use the COS formula to approximate expectations. vt, x) =E [vt, X T )] = vt, y)f y x)dy R 11 / 29
25 1D-COS formula We use the COS formula to approximate expectations. vt, x) =E [vt, X T )] = vt, y)f y x)dy b a R vt, y)f y x)dy 11 / 29
26 1D-COS formula We use the COS formula to approximate expectations. vt, x) =E [vt, X T )] = vt, y)f y x)dy b a = 2 R vt, y)f y x)dy k=0 V k T )F k x), with coefficients V k T ):= 2 b vt, y)cos kπ y a a F k x) := 2 b f y x)cos kπ y a a Exponential convergence in N for smooth density f y x). ) dy, 12) ) dy. 13) 11 / 29
27 1D-COS formula We use the COS formula to approximate expectations. vt, x) =E [vt, X T )] = vt, y)f y x)dy b a 2 R vt, y)f y x)dy N 1 k=0 V k T )F k x), with coefficients V k T ):= 2 b vt, y)cos kπ y a a F k x) := 2 b f y x)cos kπ y a a Exponential convergence in N for smooth density f y x). ) dy, 12) ) dy. 13) 11 / 29
28 Fourier coefficients F k : F k x) = b 2 2 a R f y x)cos kπ y a ) dy f y x)cos kπ y a ) dy 14) 12 / 29
29 Fourier coefficients F k : F k x) = = b f y x)cos 2 a 2 f y x)cos R 2 Re f y x)exp R kπ y a ) dy kπ y a ) dy ikπ y a ) ) dy 14) 12 / 29
30 Fourier coefficients F k : 2 b F k x) = f y x)cos kπ y a ) dy a 2 f y x)cos kπ y a ) dy R 2 = Re f y x)exp ikπ y a ) ) dy R 2 kπ ) ) = Re ϕ x e ikπ a b a 14) Characteristic function of random variable Y : ϕu) =E [exp iuy )] = expiuy)f y)dy. 15) For many asset price processes the characteristic function is available. 12 / 29 R
31 Fourier coefficients F k : 2 b F k x) = f y x)cos kπ y a ) dy a 2 f y x)cos kπ y a ) dy R 2 = Re f y x)exp ikπ y a ) ) dy R 2 kπ ) ) = Re ϕ x e ikπ a b a ) ) 2 kπ = Re x a ikπ ϕ levy e b a. 14) Characteristic function of random variable Y : ϕu) =E [exp iuy )] = expiuy)f y)dy. 15) For many asset price processes the characteristic function is available. 12 / 29 R
32 2D-COS formula In 1D: In 2D: vt, x) =e rδt E [vt, X T )] N 1 2 e rδt V k T )F k x). 16) vt, x) =e rδt E [vt, X T )] k=0
33 2D-COS formula In 1D: In 2D: vt, x) =e rδt E [vt, X T )] N 1 2 e rδt V k T )F k x). 16) vt, x) =e rδt E [vt, X T )] k=0
34 2D-COS formula In 1D: In 2D: vt, x) =e rδt E [vt, X T )] N 1 2 e rδt V k T )F k x). 16) vt, x) =e rδt E [vt, X T )] = e rδt vt, y)f y x)dy b 1 a 1 2 R b 2 a 2 2 k=0 N 1 1 N 2 1 e rδt k 1 =0 k 2 =0 This can be extended to higher dimensions. V k1,k 2 T )F k1,k 2 x). 17) 13 / 29
35 Coefficients - most difficult part. ) ) F k1,k 2 x) 2 2 b 1 a 1 b 2 a 2 f y x)cos k 1 π R y1 a1 b 1 a 1 cos k 2 π y2 a2 b 2 a 2 dy 1 dy )
36 Coefficients - most difficult part. ) ) F k1,k 2 x) 2 2 b 1 a 1 b 2 a 2 f y x)cos k 1 π R y1 a1 b 1 a 1 cos k 2 π y2 a2 b 2 a 2 dy 1 dy ) We use the following goniometric relation: 2cosα)cosβ) =cosα + β)+cosα β). 19)
37 Coefficients - most difficult part. ) ) F k1,k 2 x) 2 2 b 1 a 1 b 2 a 2 f y x)cos k 1 π R y1 a1 b 1 a 1 cos k 2 π y2 a2 b 2 a 2 dy 1 dy ) We use the following goniometric relation: 2cosα)cosβ) =cosα + β)+cosα β). 19) Then 2F k1,k 2 x) =F + k 1,k 2 x)+f k 1,k 2 x), 20) where ) F ± k 1,k 2 x) := 2 2 b 1 a 1 b 2 a 2 f y x)cos k 1 π R y1 a1 b 1 a 1 ± k 2 π y2 a2 b 2 a 2 dy 1 dy )
38 Coefficients - most difficult part. ) ) F k1,k 2 x) 2 2 b 1 a 1 b 2 a 2 f y x)cos k 1 π R y1 a1 b 1 a 1 cos k 2 π y2 a2 b 2 a 2 dy 1 dy ) We use the following goniometric relation: Then where F ± k 1,k 2 x) := 2 b 1 a 1 2 2cosα)cosβ) =cosα + β)+cosα β). 19) = 2 b 1 a 1 2 b 2 a 2 Re 2F k1,k 2 x) =F + k 1,k 2 x)+f k 1,k 2 x), 20) b 2 a 2 f y x)cos R 2 R 2 f y x)exp ) k 1 π y1 a1 b 1 a 1 ± k 2 π y2 a2 b 2 a 2 dy 1 dy 2 ) ) ik 1 π y1 a1 b 1 a 1 ± ik 2 π y2 ±a2 b 2 a 2 dy 1 dy 2 21)
39 Coefficients - most difficult part. ) ) F k1,k 2 x) 2 2 b 1 a 1 b 2 a 2 f y x)cos k 1 π R y1 a1 b 1 a 1 cos k 2 π y2 a2 b 2 a 2 dy 1 dy ) We use the following goniometric relation: Then where F ± k 1,k 2 x) := 2 b 1 a 1 2 2cosα)cosβ) =cosα + β)+cosα β). 19) = 2 2 b 1 a 1 b 2 a 2 Re = 2 b 1 a 1 2 b 2 a 2 Re 2F k1,k 2 x) =F + k 1,k 2 x)+f k 1,k 2 x), 20) b 2 a 2 f y x)cos R 2 ϕ ) k 1 π y1 a1 b 1 a 1 ± k 2 π y2 a2 b 2 a 2 dy 1 dy 2 ik 1 π y1 a1 b 1 a 1 f y x)exp R 2 ) k1π b 1 a 1, ± k2π x b 2 a 2 exp ) ) ± ik 2 π y2 ±a2 b 2 a 2 dy 1 dy 2 ik 1 π a1 b 1 a 1 ik 2 π a2 21) b 2 a 2 ))
40 Approximate the terminal coefficients V k1,k 2 T ) with DCTs. Take Q max[n 1, N 2 ] grid-points and y n i i := a i +n i )b i a i Q and Δy i := b i a i, i =1, 2. 22) Q 15 / 29
41 Approximate the terminal coefficients V k1,k 2 T ) with DCTs. Take Q max[n 1, N 2 ] grid-points and y n i i := a i +n i )b i a i Q and The midpoint-rule integration gives us V k1,k 2 T ) = 2 b 1 a 1 2 b 2 a 2 b2 Q 1 Q 1 n 1 =0 n 2 =0 a 2 b1 a 1 gy n 1 1, y n 2 2 )cos Δy i := b i a i, i =1, 2. 22) Q ) ) gy)cos k 1 π y 1 a 1 b 1 a 1 cos k 2 π y 2 a 2 b 2 a 2 dy 1 dy 2 k 1 π 2n 1+1 2Q ) cos k 2 π 2n 2+1 2Q ) b1 a 1 Q b 2 a 2 Q. The above 2D-DCT can be calculated efficiently by, for example, MATLAB s function dct2. 15 / 29
42 Results - European options Results geometric basket call under correlated geometric Brownian motion, vt 0, x 0 )= Table 1: N 2 = N 1 ). N Error 7.60e e e e e-13 CPU ms) log 10 error N 1 16 / 29
43 Jump-diffusion process The log-jump-diffusion process dst i =r λe[ej i 1])St i dt + σ ist i dw t i + S t i ej i 1)dq t, 23) with q t a Poisson process with intenstity λ, andj =J 1, J 2 ) bivariate normally distributed jumps. Figure 4: Density recovery ˆf X T x 0 ). 17 / 29
44 Jump-diffusion process The log-jump-diffusion process dst i =r λe[ej i 1])St i dt + σ ist i dw t i + S t i ej i 1)dq t, 23) with q t a Poisson process with intenstity λ, andj =J 1, J 2 ) bivariate normally distributed jumps. Table 2: Put-on-min option values ˆvt 0, x 0 )N 1 = N 2 = 125), they correspond to value in Clift, 2008). Figure 4: Density recovery ˆf X T x 0 ). S 2 0 S / 29
45 Spark spread option - 3D 3-dimensional GBM S t =[St power, St gas, S CO 2 t ]. Spark spread: net revenue from selling power. Payoff : gs T )=Ωmax ) S power T α g S gas T αco 2 S CO 2 T K, / 29
46 Spark spread option - 3D 3-dimensional GBM S t =[St power, St gas, S CO 2 t ]. Spark spread: net revenue from selling power. Payoff : gs T )=Ωmax ) S power T α g S gas T αco 2 S CO 2 T K, / 29
47 Spark spread option - 3D 3-dimensional GBM S t =[St power, St gas, S CO 2 t ]. Spark spread: net revenue from selling power. Payoff : gs T )=Ωmax ) S power T α g S gas T αco 2 S CO 2 T K, 0. vt 0, S 0 ) CPU time s) N N Q n/a n/a Calculation of the option s Greeks is straightforward. 18 / 29
48 Bermudan options Bermudan option: fixed exercise dates t m m =1,...,M) atwhich you can either exercise the option or continue. Option value in 1D: vt m, x) =max[gx), ct m, x)]. 24) 19 / 29
49 Bermudan options Bermudan option: fixed exercise dates t m m =1,...,M) atwhich you can either exercise the option or continue. Option value in 1D: vt m, x) =max[gx), ct m, x)]. 24) Coefficients V k at time t m : V k t m ):= 2 b a b a vt m, y)cos kπ y a b a ) dy. 25) 19 / 29
50 Bermudan options Bermudan option: fixed exercise dates t m m =1,...,M) atwhich you can either exercise the option or continue. Option value in 1D: vt m, x) =max[gx), ct m, x)]. 24) Coefficients V k at time t m : V k t m ):= 2 b a b a vt m, y)cos kπ y a b a ) dy. 25) Asset price X t Continue 85 Exercise Time t m 19 / 29
51 Bermudan option - 2D Coefficients V k1,k 2 at time t m : V k1,k 2 t m ):= with b2 b1 a 2 a 1 vt m, y)cos ) ) k 1 π y 1 a 1 b 1 a 1 cos k 2 π y 2 a 2 b 2 a 2 dy 1 dy 2 vt m, x) =max[gx), ct m, x)]. 26) 20 / 29
52 Bermudan option - 2D Coefficients V k1,k 2 at time t m : V k1,k 2 t m ):= with b2 b1 a 2 a 1 vt m, y)cos ) ) k 1 π y 1 a 1 b 1 a 1 cos k 2 π y 2 a 2 b 2 a 2 dy 1 dy 2 vt m, x) =max[gx), ct m, x)]. 26) 6 Left: Optimal exercise domains blue) and continuation domains green) at initial time t 0. J rectangular sub-domains. y y 1 20 / 29
53 Equidistant vs. non-equidistant grid x x 1 Figure 5: Equidistant grid. 21 / 29
54 Equidistant vs. non-equidistant grid x 2 x x x 1 Figure 5: Equidistant grid. Figure 6: Non-equidistant grid. 21 / 29
55 Recursive recovery V k1,k 2 t m )= b2 b1 a 2 = p + q a 1 vt m, y)cos k 1 π y 1 a 1 gy)cos G p C q ct m, y)cos b 1 a 1 ) cos k 1 π y 1 a 1 b 1 a 1 ) cos k 2 π y 2 a 2 b 2 a 2 ) dy 1 dy 2 k 2 π y 2 a 2 b 2 a 2 ) dy ) ) k 1 π y 1 a 1 b 1 a 1 cos k 2 π y 2 a 2 b 2 a 2 dy 22 / 29
56 Recursive recovery V k1,k 2 t m )= b2 b1 a 2 = p + q a 1 vt m, y)cos k 1 π y 1 a 1 gy)cos G p C q ĉt m, y)cos b 1 a 1 ) cos k 1 π y 1 a 1 b 1 a 1 ) cos k 2 π y 2 a 2 b 2 a 2 ) dy 1 dy 2 k 2 π y 2 a 2 b 2 a 2 ) dy ) ) k 1 π y 1 a 1 b 1 a 1 cos k 2 π y 2 a 2 b 2 a 2 dy 22 / 29
57 Recursive recovery V k1,k 2 t m )= b2 b1 a 2 = p + q a 1 vt m, y)cos k 1 π y 1 a 1 gy)cos G p C q ĉt m, y)cos b 1 a 1 ) cos k 1 π y 1 a 1 b 1 a 1 ) cos k 2 π y 2 a 2 b 2 a 2 ) dy 1 dy 2 k 2 π y 2 a 2 b 2 a 2 ) dy ) ) k 1 π y 1 a 1 b 1 a 1 cos k 2 π y 2 a 2 b 2 a 2 dy The resulting matrix-vector products Mu can be computed efficiently by a Fourier-based algorithm. The computation time achieved is ON log 2 N). 22 / 29
58 Algorithm 2D-COS method for pricing Bermudan rainbow options Initialisation: Calculate coefficients V k1,k 2 t M ). Main loop to recover ˆV t m ): For m = M 1to1: Determine the optimal continuation regions C q and early-exercise regions G p. Compute ˆV t m ) with the help of the FFT algorithm. Final step: Compute ˆvt 0, x 0 ) by inserting ˆV k1,k 2 t 1 ) into COS formula. 23 / 29
59 Results Bermudan option Error geometric basket option under GBM. N 1 is number of terms in series expansion N 2 = N 1 ). J is number of sub-domains. 24 / 29
60 Results Bermudan option Error geometric basket option under GBM. N 1 is number of terms in series expansion N 2 = N 1 ). J is number of sub-domains. log 10 error N J / 29
61 Results Bermudan option Error geometric basket option under GBM. N 1 is number of terms in series expansion N 2 = N 1 ). J is number of sub-domains. Call on maximum option log 10 error N J Table 3: Results call-on-max GBM). Andersen 2004) Shashi S 0 2D-COS Binomial SGM Direct ) ) ) 24 / 29
62 Heston model X t represents the asset price process and ν t is the variance process, with dwt 1 dwt 2 = ρdt) dx t =r 1 2 ν t)x t dt + ν t X t dwt 1, dν t = κ ν ν t )dt + η ν t dwt 2. 27) 25 / 29
63 Heston model X t represents the asset price process and ν t is the variance process, with dwt 1 dwt 2 = ρdt) dx t =r 1 2 ν t)x t dt + ν t X t dwt 1, dν t = κ ν ν t )dt + η ν t dwt 2. 27) 25 / 29
64 Heston model X t represents the asset price process and ν t is the variance process, with dwt 1 dwt 2 = ρdt) dx t =r 1 2 ν t)x t dt + ν t X t dwt 1, dν t = κ ν ν t )dt + η ν t dwt 2. 27) The variance process remains strictly positive if the Feller condition is satisfied, 2κū/η 2 1:=q Feller 0, otherwise it may reach zero. 25 / 29
65 European option, with Bermudan framework M = 12 time steps. 26 / 29
66 European option, with Bermudan framework M = 12 time steps. Feller satisfied q Feller =0.98), v = N e e e e e e e e-6 N e e e e-6 26 / 29
67 European option, with Bermudan framework M = 12 time steps. Feller satisfied q Feller =0.98), v = N e e e e e e e e-6 N e e e e-6 Feller not satisfied q Feller = 0.47), v = N e e e e e e e e-5 N e e e e-5 26 / 29
68 Feller not satisfied q Feller = 0.84), v = N e e e e e e e e-2 N e e e e-2 27 / 29
69 Bermudan put option N 1 Feller satisfied q Feller =0.98) N / 29
70 Bermudan put option N 1 Feller satisfied q Feller =0.98) N N 1 Feller not satisfied q Feller = 0.47) N / 29
71 Summary and conclusion COS method is based on Fourier-cosine series expansion. Exponential convergence for f C. Can be extended to higher dimensions for pricing rainbow options. Experiments with financial and spark options. Heston stochastic volatility model. Future research Higher dimensions, Gibbs phenomenon and filters, Asian options. Contact: m.j.ruijter@cwi.nl 29 / 29
Time-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 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 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 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 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 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 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 informationImportance Sampling for Option Pricing. Steven R. Dunbar. Put Options. Monte Carlo Method. Importance. Sampling. Examples.
for for January 25, 2016 1 / 26 Outline for 1 2 3 4 2 / 26 Put Option for A put option is the right to sell an asset at an established price at a certain time. The established price is the strike price,
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 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 informationMonte Carlo Simulations
Monte Carlo Simulations Lecture 1 December 7, 2014 Outline Monte Carlo Methods Monte Carlo methods simulate the random behavior underlying the financial models Remember: When pricing you must simulate
More informationNEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 MAS3904. Stochastic Financial Modelling. Time allowed: 2 hours
NEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 Stochastic Financial Modelling Time allowed: 2 hours Candidates should attempt all questions. Marks for each question
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 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 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 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 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 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 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 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 informationRough Heston models: Pricing, hedging and microstructural foundations
Rough Heston models: Pricing, hedging and microstructural foundations Omar El Euch 1, Jim Gatheral 2 and Mathieu Rosenbaum 1 1 École Polytechnique, 2 City University of New York 7 November 2017 O. El Euch,
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 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 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 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 informationEfficient Valuation of Equity-Indexed Annuities Under Lévy Processes Using Fourier-Cosine Series
Efficient Valuation of Equity-Indexed Annuities Under Lévy Processes Using Fourier-Cosine Series Geng Deng, PhD, CFA, FRM Tim Dulaney, PhD, FRM, Craig McCann, PhD, CFA Mike Yan, PhD, FRM April 1, 14 Abstract
More informationTo link to this article:
This article was downloaded by: [Centrum Wiskunde & Informatica] On: 24 July 2012, At: 02:56 Publisher: Taylor & Francis Informa Ltd Registered in England and Wales Registered Number: 1072954 Registered
More informationApplication of Moment Expansion Method to Option Square Root Model
Application of Moment Expansion Method to Option Square Root Model Yun Zhou Advisor: Professor Steve Heston University of Maryland May 5, 2009 1 / 19 Motivation Black-Scholes Model successfully explain
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 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 informationNumerical Methods for Pricing Energy Derivatives, including Swing Options, in the Presence of Jumps
Numerical Methods for Pricing Energy Derivatives, including Swing Options, in the Presence of Jumps, Senior Quantitative Analyst Motivation: Swing Options An electricity or gas SUPPLIER needs to be capable,
More informationAnumericalalgorithm for general HJB equations : a jump-constrained BSDE approach
Anumericalalgorithm for general HJB equations : a jump-constrained BSDE approach Nicolas Langrené Univ. Paris Diderot - Sorbonne Paris Cité, LPMA, FiME Joint work with Idris Kharroubi (Paris Dauphine),
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 informationLinearity-Generating Processes, Unspanned Stochastic Volatility, and Interest-Rate Option Pricing
Linearity-Generating Processes, Unspanned Stochastic Volatility, and Interest-Rate Option Pricing Liuren Wu, Baruch College Joint work with Peter Carr and Xavier Gabaix at New York University Board of
More informationApplication of Stochastic Calculus to Price a Quanto Spread
Application of Stochastic Calculus to Price a Quanto Spread Christopher Ting http://www.mysmu.edu/faculty/christophert/ Algorithmic Quantitative Finance July 15, 2017 Christopher Ting July 15, 2017 1/33
More information3.1 Itô s Lemma for Continuous Stochastic Variables
Lecture 3 Log Normal Distribution 3.1 Itô s Lemma for Continuous Stochastic Variables Mathematical Finance is about pricing (or valuing) financial contracts, and in particular those contracts which depend
More informationNumerical treatment of stochastic control problems by Fourier-cosine series expansions
Delft University of Technology Faculty of Electrical Engineering, Mathematics and Computer Science Delft Institute of Applied Mathematics Numerical treatment of stochastic control problems by Fourier-cosine
More informationLecture 2: Rough Heston models: Pricing and hedging
Lecture 2: Rough Heston models: Pricing and hedging Mathieu Rosenbaum École Polytechnique European Summer School in Financial Mathematics, Dresden 217 29 August 217 Mathieu Rosenbaum Rough Heston models
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 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 informationMASM006 UNIVERSITY OF EXETER SCHOOL OF ENGINEERING, COMPUTER SCIENCE AND MATHEMATICS MATHEMATICAL SCIENCES FINANCIAL MATHEMATICS.
MASM006 UNIVERSITY OF EXETER SCHOOL OF ENGINEERING, COMPUTER SCIENCE AND MATHEMATICS MATHEMATICAL SCIENCES FINANCIAL MATHEMATICS May/June 2006 Time allowed: 2 HOURS. Examiner: Dr N.P. Byott This is a CLOSED
More informationModelling Credit Spread Behaviour. FIRST Credit, Insurance and Risk. Angelo Arvanitis, Jon Gregory, Jean-Paul Laurent
Modelling Credit Spread Behaviour Insurance and Angelo Arvanitis, Jon Gregory, Jean-Paul Laurent ICBI Counterparty & Default Forum 29 September 1999, Paris Overview Part I Need for Credit Models Part II
More informationPricing and Modelling in Electricity Markets
Pricing and Modelling in Electricity Markets Ben Hambly Mathematical Institute University of Oxford Pricing and Modelling in Electricity Markets p. 1 Electricity prices Over the past 20 years a number
More informationSupplementary Appendix to The Risk Premia Embedded in Index Options
Supplementary Appendix to The Risk Premia Embedded in Index Options Torben G. Andersen Nicola Fusari Viktor Todorov December 214 Contents A The Non-Linear Factor Structure of Option Surfaces 2 B Additional
More informationLecture Note 8 of Bus 41202, Spring 2017: Stochastic Diffusion Equation & Option Pricing
Lecture Note 8 of Bus 41202, Spring 2017: Stochastic Diffusion Equation & Option Pricing We shall go over this note quickly due to time constraints. Key concept: Ito s lemma Stock Options: A contract giving
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 informationIEOR E4703: Monte-Carlo Simulation
IEOR E4703: Monte-Carlo Simulation Generating Random Variables and Stochastic Processes Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com
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 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 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 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 informationPRICING EUROPEAN OPTIONS BASED ON
INSTITUTO SUPERIOR DE CIÊNCIAS DO TRABALHO E DA EMPRESA FACULDADE DE CIÊNCIAS DA UNVIVERSIDADE DE LISBOA DEPARTMENT OF FINANCE DEPARTMENT OF MATHEMATICS PRICING EUROPEAN OPTIONS BASED ON FOURIER-COSINE
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 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 informationChapter 3 Common Families of Distributions. Definition 3.4.1: A family of pmfs or pdfs is called exponential family if it can be expressed as
Lecture 0 on BST 63: Statistical Theory I Kui Zhang, 09/9/008 Review for the previous lecture Definition: Several continuous distributions, including uniform, gamma, normal, Beta, Cauchy, double exponential
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 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 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 informationHeston Stochastic Local Volatility Model
Heston Stochastic Local Volatility Model Klaus Spanderen 1 R/Finance 2016 University of Illinois, Chicago May 20-21, 2016 1 Joint work with Johannes Göttker-Schnetmann Klaus Spanderen Heston Stochastic
More informationMathematics of Finance Final Preparation December 19. To be thoroughly prepared for the final exam, you should
Mathematics of Finance Final Preparation December 19 To be thoroughly prepared for the final exam, you should 1. know how to do the homework problems. 2. be able to provide (correct and complete!) definitions
More informationIMPA Commodities Course : Forward Price Models
IMPA Commodities Course : Forward Price Models Sebastian Jaimungal sebastian.jaimungal@utoronto.ca Department of Statistics and Mathematical Finance Program, University of Toronto, Toronto, Canada http://www.utstat.utoronto.ca/sjaimung
More informationPartial differential approach for continuous models. Closed form pricing formulas for discretely monitored models
Advanced Topics in Derivative Pricing Models Topic 3 - Derivatives with averaging style payoffs 3.1 Pricing models of Asian options Partial differential approach for continuous models Closed form pricing
More informationDefinition Pricing Risk management Second generation barrier options. Barrier Options. Arfima Financial Solutions
Arfima Financial Solutions Contents Definition 1 Definition 2 3 4 Contenido Definition 1 Definition 2 3 4 Definition Definition: A barrier option is an option on the underlying asset that is activated
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 informationAn Analytical Approximation for Pricing VWAP Options
.... An Analytical Approximation for Pricing VWAP Options Hideharu Funahashi and Masaaki Kijima Graduate School of Social Sciences, Tokyo Metropolitan University September 4, 215 Kijima (TMU Pricing of
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 informationMLEMVD: A R Package for Maximum Likelihood Estimation of Multivariate Diffusion Models
MLEMVD: A R Package for Maximum Likelihood Estimation of Multivariate Diffusion Models Matthew Dixon and Tao Wu 1 Illinois Institute of Technology May 19th 2017 1 https://papers.ssrn.com/sol3/papers.cfm?abstract
More informationCONTINUOUS TIME PRICING AND TRADING: A REVIEW, WITH SOME EXTRA PIECES
CONTINUOUS TIME PRICING AND TRADING: A REVIEW, WITH SOME EXTRA PIECES THE SOURCE OF A PRICE IS ALWAYS A TRADING STRATEGY SPECIAL CASES WHERE TRADING STRATEGY IS INDEPENDENT OF PROBABILITY MEASURE COMPLETENESS,
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 informationLattice (Binomial Trees) Version 1.2
Lattice (Binomial Trees) Version 1. 1 Introduction This plug-in implements different binomial trees approximations for pricing contingent claims and allows Fairmat to use some of the most popular binomial
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 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 informationWKB Method for Swaption Smile
WKB Method for Swaption Smile Andrew Lesniewski BNP Paribas New York February 7 2002 Abstract We study a three-parameter stochastic volatility model originally proposed by P. Hagan for the forward swap
More informationStochastic Dynamical Systems and SDE s. An Informal Introduction
Stochastic Dynamical Systems and SDE s An Informal Introduction Olav Kallenberg Graduate Student Seminar, April 18, 2012 1 / 33 2 / 33 Simple recursion: Deterministic system, discrete time x n+1 = f (x
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 informationThe Valuation of Bermudan Guaranteed Return Contracts
The Valuation of Bermudan Guaranteed Return Contracts Steven Simon 1 November 2003 1 K.U.Leuven and Ente Luigi Einaudi Abstract A guaranteed or minimum return can be found in different financial products,
More informationEfficient Pricing of European-Style Asian Options under Exponential Lévy Processes Based on Fourier Cosine Expansions
SIAM J. FINANCIAL MATH. Vol. 4, pp. 399 426 c 2013 Society for Industrial and Applied Mathematics Efficient Pricing of European-Style Asian Options under Exponential Lévy Processes Based on Fourier Cosine
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 informationPricing and hedging with rough-heston models
Pricing and hedging with rough-heston models Omar El Euch, Mathieu Rosenbaum Ecole Polytechnique 1 January 216 El Euch, Rosenbaum Pricing and hedging with rough-heston models 1 Table of contents Introduction
More informationIntroduction Credit risk
A structural credit risk model with a reduced-form default trigger Applications to finance and insurance Mathieu Boudreault, M.Sc.,., F.S.A. Ph.D. Candidate, HEC Montréal Montréal, Québec Introduction
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 informationThe Pennsylvania State University. The Graduate School. Department of Industrial Engineering AMERICAN-ASIAN OPTION PRICING BASED ON MONTE CARLO
The Pennsylvania State University The Graduate School Department of Industrial Engineering AMERICAN-ASIAN OPTION PRICING BASED ON MONTE CARLO SIMULATION METHOD A Thesis in Industrial Engineering and Operations
More informationContinous time models and realized variance: Simulations
Continous time models and realized variance: Simulations Asger Lunde Professor Department of Economics and Business Aarhus University September 26, 2016 Continuous-time Stochastic Process: SDEs Building
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 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 informationEconomics has never been a science - and it is even less now than a few years ago. Paul Samuelson. Funeral by funeral, theory advances Paul Samuelson
Economics has never been a science - and it is even less now than a few years ago. Paul Samuelson Funeral by funeral, theory advances Paul Samuelson Economics is extremely useful as a form of employment
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 informationWrite legibly. Unreadable answers are worthless.
MMF 2021 Final Exam 1 December 2016. This is a closed-book exam: no books, no notes, no calculators, no phones, no tablets, no computers (of any kind) allowed. Do NOT turn this page over until you are
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 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 informationMonte Carlo Based Numerical Pricing of Multiple Strike-Reset Options
Monte Carlo Based Numerical Pricing of Multiple Strike-Reset Options Stavros Christodoulou Linacre College University of Oxford MSc Thesis Trinity 2011 Contents List of figures ii Introduction 2 1 Strike
More informationUser s Guide for the Matlab Library Implementing Closed Form MLE for Diffusions
User s Guide for the Matlab Library Implementing Closed Form MLE for Diffusions Yacine Aït-Sahalia Department of Economics and Bendheim Center for Finance Princeton University and NBER This Version: July
More informationMarkovian Projection, Heston Model and Pricing of European Basket Optio
Markovian Projection, Heston Model and Pricing of European Basket Options with Smile July 7, 2009 European Options on geometric baskets European Options on arithmetic baskets Pricing via Moment Matching
More informationStatistica Sinica Preprint No: SS R2
0 Statistica Sinica Preprint No: SS-2016-0434R2 Title Multi-asset empirical martingale price estimators derivatives Manuscript ID SS-2016-0434R2 URL http://www.stat.sinica.edu.tw/statistica/ DOI 10.5705/ss.202016.0434
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 informationLimit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies
Limit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies George Tauchen Duke University Viktor Todorov Northwestern University 2013 Motivation
More informationOptimized Least-squares Monte Carlo (OLSM) for Measuring Counterparty Credit Exposure of American-style Options
Optimized Least-squares Monte Carlo (OLSM) for Measuring Counterparty Credit Exposure of American-style Options Kin Hung (Felix) Kan 1 Greg Frank 3 Victor Mozgin 3 Mark Reesor 2 1 Department of Applied
More informationEstimating the Greeks
IEOR E4703: Monte-Carlo Simulation Columbia University Estimating the Greeks c 207 by Martin Haugh In these lecture notes we discuss the use of Monte-Carlo simulation for the estimation of sensitivities
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 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 information