Math 416/516: Stochastic Simulation
|
|
- Brett Bryan
- 6 years ago
- Views:
Transcription
1 Math 416/516: Stochastic Simulation Haijun Li Department of Mathematics Washington State University Week 13 Haijun Li Math 416/516: Stochastic Simulation Week 13 1 / 28
2 Outline 1 Simulation of Itō Diffusions: The Euler Approximation 2 Second Order Approximation: The Milstein Scheme 3 A Summary: Monte Carlo Methods in Sciences and Engineering Haijun Li Math 416/516: Stochastic Simulation Week 13 2 / 28
3 Numerical Solution of Stochastic Differential Equations SDEs which admit an explicit solution are few exceptions. Therefore numerical techniques for the approximation of the solution to a SDE are often called for. One purpose is to visualize a variety of sample paths of the solution. A collection of such paths is called a scenario, which can be used for some kind of prediction of the stochastic process at future instants of time. A second objective is to achieve reasonable approximations to the distributional quantities (expectations, variances, covariance and higher-order moments) of the solution to a SDE. Only in a few cases one is able to give explicit formulas for these quantities, and even then they frequently involve special functions which have to be approximated numerically. Numerical solutions allow us to simulate as many sample paths as we want; they constitute the basis for Monte-Carlo techniques to obtain the distributional characteristics of Itō diffusions. Haijun Li Math 416/516: Stochastic Simulation Week 13 3 / 28
4 Euler Approximation Scheme Consider the SDE dx t = µ(x t )dt + σ(x t )db t, t [0, T ]. We assume that the coefficient functions µ(x) and σ(x) are Lipschitz continuous, and EX0 2 <, which guarantee the existence and uniqueness of a (strong) solution. 1 To approximate the solution, partition [0, T ] as follows, τ m : 0 = t 0 < t 1 < < t m 1 < t m = T, with i = t i t i 1, 1 i m, and mesh(τ m ) = max 1 i m i. Let i B = B ti B ti 1, 1 i m. 2 Define recursively, 1 i m, X ti = X ti 1 + µ(x ti 1 ) i + σ(x ti 1 ) i B, with initial value X 0. Haijun Li Math 416/516: Stochastic Simulation Week 13 4 / 28
5 Estimation via Euler Scheme Goal: To estimate v = E[h(X t, 0 t T )], where (X t, 0 t T ) is the solution of the SDE dx t = µ(x t )dt + σ(x t )db t, t [0, T ]. Generate n independent sample paths, 1 k n, X (k) t i = X (k) t i 1 + µ(x (k) t i 1 )(t i t i 1 ) + σ(x (k) t i 1 ) t i t i 1 Z (k) with X (k) 0 = X 0, where Z k 1,..., Z k m are i.i.d. standard normal random variables. Use the estimator ˆv n = 1 n approx of (X t,0 t T ) n {}}{ E[h( X (k) t 0, X (k) t 1,..., X (k) t m )] k=1 i, 1 i m Haijun Li Math 416/516: Stochastic Simulation Week 13 5 / 28
6 Example Consider the following SDE ds t = rs t dt + σs t db t where r is the interest rate and σ is the volatility. Use the Euler scheme: S (k) t 0 = S 0, S (k) t i = S (k) t i 1 + rs (k) t i 1 (t i t i 1 ) + σs (k) t ti i 1 t i 1 Z (k) The estimate of a European call option price: ˆv n = 1 n e rtm (S (k) t n m K ) + k=1 The estimate of an Asian call option price: ( ) + ˆv n = 1 n 1 m e rtm S (k) t n m i K. k=1 i=1 i, 1 i m Haijun Li Math 416/516: Stochastic Simulation Week 13 6 / 28
7 Discretization Error for v = E[h(X t, 0 t T )] Since we use a time-discretized version of Itō diffusion to approximate the continuous-time Itō diffusion, this would create bias, called the discritization error. That is, in general, ˆv n is a biased estimator. MSE (or L 2 -distance) = E(ˆv n v) 2 = var(ˆv n ) + [E(ˆv n ) v] 2 }{{} bias 2 Such discretization error cannot be eliminated by merely increasing the sample size n. The variance reduction techniques decreases the variance var(ˆv n ), but would not be able to eliminate the bias. Under some mild regularity conditions, when the time discretization becomes finer (m becomes larger), the discretization error converges to zero. Haijun Li Math 416/516: Stochastic Simulation Week 13 7 / 28
8 Convergence Rate of Numerical Solutions The Euler Approximation Rate is 0.5. Consider the SDE dx t = µ(x t )dt + σ(x t )db t, t [0, T ]. Assume that the coefficient functions µ(x) and σ(x) are Lipschitz continuous, and EX0 2 <. The Euler approximation converges strongly with order 0.5. That is, there exists a constant c > 0 such that E X T X (m) T c mesh(τ m ) 0.5. ( 0, as m ) One could use E sup 0 t T X t X (m) t as a more appropriate criteria to describe the pathwise closeness of X and X (m). But this quantity is more difficult to deal with theoretically. Haijun Li Math 416/516: Stochastic Simulation Week 13 8 / 28
9 160 CHAPTER 3. Numerical Solutions (dashed lines) VS Exact Solution LT ,! I 8,,,, t Figure : dx t = 0.01X t dt X t db t, X 0 = 1. Haijun Li Math 416/516: Stochastic Simulation Week 13 9 / 28
10 Example: Eliminating Discretization Error Consider the Itō diffusion X t that satisfies the SDE dx t = r dt + θ(t)db t, 0 t T, where r is a constant and θ(t) is a deterministic function of t. Rewrite: t X t = X 0 + r t + θ(s)db s 0 Using Itō s Isometry, we have t ( t ) θ(s)db s N 0, θ(s) 2 ds 0 0 When t = T, ( 1/2 T X T = X 0 + r T + θ(s) ds) 2 Z where Z has the standard normal distribution. Haijun Li Math 416/516: Stochastic Simulation Week / 28 0
11 Example: Eliminating Discretization Error (cont d) For a European-type call option, the estimate of the price: ˆv n = 1 n n k=1 e rt (X (k) T K ) + v = E[e rt (X T K ) + ], which does not involve discretization, and has no discretization error. For an Asian-type call option, the estimate of the price: ( ) + ˆv n = 1 n e rt 1 m X (k) t n m i K k=1 i=1 ( ) + ] v = E [e rt 1 T X s ds K, T 0 which still needs discretization. Haijun Li Math 416/516: Stochastic Simulation Week / 28
12 Milstein Approximation Scheme Consider the SDE dx t = µ(x t )dt + σ(x t )db t, t [0, T ]. In contrast to the first order approximation, the Milstein approximation exploits a so-called Taylor-Itō expansion that incorporates high order approximation. Heuristics: Apply the Itō lemma to the integrands µ(x s ) and σ(x s ) at each point t i 1 of discretization, and then estimate the higher order terms using separable multiple stochastic integrals. Taylor-Itō expansions involve multiple stochastic integrals. Their rigorous treatment requires a more advanced theory of the stochastic calculus. Haijun Li Math 416/516: Stochastic Simulation Week / 28
13 Milstein Approximation (Grigori Milstein, 1975) 1 To approximate the solution, partition [0, T ] as follows, τ m : 0 = t 0 < t 1 < < t m 1 < t m = T, with i = t i t i 1, 1 i m, and mesh(τ m ) = max 1 i m i. Let i B = B ti B ti 1, 1 i m. 2 With initial value X 0, define recursively, for 1 i m, X ti = X ti 1 + µ(x ti 1 ) i + σ(x ti 1 ) i B + 1 }{{} 2 σ(x t i 1 )σ (X ti 1 )[( i B) 2 i ]. }{{} Euler s first order approx Milstein s second order term Note that where Z N(0, 1). i B = t i t i 1 Z N(0, t i t i 1 ), ( i B) 2 i = (t i t i 1 )(Z 2 1) Haijun Li Math 416/516: Stochastic Simulation Week / 28
14 Estimation via Milstein Scheme Goal: To estimate v = E[h(X t, 0 t T )], where (X t, 0 t T ) is the solution of the SDE dx t = µ(x t )dt + σ(x t )db t, t [0, T ]. Generate n independent sample paths, 1 k n, X (k) t i = X (k) t i 1 + µ(x (k) t i 1 )(t i t i 1 ) + σ(x (k) t i 1 ) t i t i 1 Z (k) i + 1 (k) σ(x t 2 i 1 )σ (X (k) t i 1 )(t i t i 1 )[(Z (k) i ) 2 1], with X (k) 0 = X 0, where Z1 k,..., Z m k are i.i.d. standard normal random variables. Use the estimator ˆv n = 1 n approx of (X t,0 t T ) n {}}{ E[h( X (k) t 0, X (k) t 1,..., X (k) t m )] k=1 Haijun Li Math 416/516: Stochastic Simulation Week / 28
15 Convergence Rate of Milstein Scheme The Milstein Approximation Rate is 1.0. Consider the SDE dx t = µ(x t )dt + σ(x t )db t, t [0, T ]. Assume that the coefficient functions µ(x) and σ(x) are Lipschitz continuous, and EX0 2 <. The Euler approximation converges strongly with order 1.0. That is, there exists a constant c > 0 such that E X T X (m) T c mesh(τ m ) 1.0. ( 0, as m ) Haijun Li Math 416/516: Stochastic Simulation Week / 28
16 Euler (left 3.4. NUMERICAL column) SOLUTION VS Milstein (right column) 165 I I I 1C I n I Figure : dx t = 0.01X t dt + 2X t db t, X 0 = 1. Haijun Li Math 416/516: Stochastic Simulation Week / 28
17 The Lamperti Transform (John Lamperti, 1962, 1972) Consider the SDE: dx t = µ(x t )dt + σ(x t )db t, t [0, T ]. If σ(x) is constant (σ (x) = 0), then the Milstein scheme reduces to the Euler scheme. Example: Langevin Equation dx t = cx t dt + σdb t, t [0, T ]. Can we transform the SDE to a SDE with constant volatility? We standardize it! Use the Lamperti transform Y t = L(X t ), where L(x) = 1 σ(x) dx. The Lamperti transform is widely used in studies of Self-Similar phenomena and Fractals. Haijun Li Math 416/516: Stochastic Simulation Week / 28
18 Standardize SDEs via the Lamperti Transform Consider the SDE dx t = µ(x t )dt + σ(x t )db t, t [0, T ]. Take the derivatives on the Lamperti transform: Use Itō s Lemma, L (x) = 1 σ(x), L (x) = σ (x) σ 2 (x) dy t = dl(x t ) = L (X t )dx t L (X t )dt = [ ] µ(xt ) σ(x t ) σ (x) dt +db t. 2 Since X t = L 1 (Y t ), we have the standardized SDE for Y t : [ µ(l 1 (Y t )) dy t = σ(l 1 (Y t )) σ (L 1 ] (Y t )) dt + db t. 2 Haijun Li Math 416/516: Stochastic Simulation Week / 28
19 Example: Geometric Brownian Motion Consider the linear SDE: with initial value X 0 > 0. Use the Lamperti transform L(x) = x dx t = rx t dt + σx t db t, X 0 1 σs ds = 1 σ log x X 0, Y t = L(X t ) = 1 σ log X t X 0, X t = X 0 e σy t. Since dy t = ( r σ σ ) ( r dt + db t = Y t = 2 σ σ ) t + B t, 2 we obtain the geometric Brownian motion X t = X 0 e (r σ2 /2)t+σB t. Haijun Li Math 416/516: Stochastic Simulation Week / 28
20 Example: Cox-Ingersoll-Ross Process Consider the non-linear SDE: dx t = a(b X t )dt + σ X t db t, with initial value X 0 = 0. Here a, b, σ are positive constants satisfying that 2ab σ 2. Use the Lamperti transform L(x) = x 0 1 σ s ds = 2 x, Yt = L(X t ) = 2 Xt, X t = σ2 Yt 2 σ σ 4. Convert the SDE to a SDE with constant volatility dy t = ( 4ab σ 2 1 2σ 2 ay ) t dt + db t Y t 2 and then use the Euler scheme with no discretization error. Haijun Li Math 416/516: Stochastic Simulation Week / 28
21 Remark: Interest (Short) Rate Modeling The well-known Vasicek interest rate model (Oldrich Vasicek, Wells Fargo, 1977) describes the interest rate X t as an Ornstein-Uhlenbeck process. Haijun Li Math 416/516: Stochastic Simulation Week / 28
22 Remark: Interest Rate Modeling The basic Vasicek interest rate model was extended to the mean-reverting Itō diffusion. Haijun Li Math 416/516: Stochastic Simulation Week / 28
23 Remark: Interest Rate Modeling The Cox-Ingersoll-Ross model (John Cox, Jonathan Ingersoll and Stephen Ross, 1985; also see Lin Chen, US Federal Reserve, 1994) describes the evolution of interest rates, as an extension of the well-known Vasicek model. Haijun Li Math 416/516: Stochastic Simulation Week / 28
24 Monte Carlo vs Numerical Methods Once sample paths (or scenarios) of the solution of a SDE are obtained, they can be used to estimate the distributional quantities (expectations, variances, higher-order moments) of the solution. Since derivative prices are often written as expectations of underlying asset values, which are the solutions of SDEs, Monte Carlo method becomes an essential tool in the pricing of derivative securities and in risk management. Monte Carlo is generally not a competitive method for calculating univariate expectation. For example, the error in a trapezoidal rule for the integral of a d-dimensional twice continuously differentiable function is O(n 2/d ), which is in contrast to the standard error O(n 1/2 ) of the Monte Carlo method for the same problem. The performance degradation with increasing dimension is a characteristic of all deterministic integration methods, and thus Monte Carlo methods are attractive in high dimension. Haijun Li Math 416/516: Stochastic Simulation Week / 28
25 Black-Scholes Model Assume that the stock price follows a geometric Brownian motion: that is, S t = S 0 e (c σ2 /2)t+σB t, t 0, S t LogN(log S 0 + (c σ 2 /2)t, σ 2 t), where the drift c and volatility σ are the (physical) parameters of the stock price process. The risk-free interest rate is given by r. Under the risk-neutral probability measure, S t LogN(log S 0 + (r σ 2 /2)t, σ 2 t). Fundamental Theorem of Arbitrage-Free Pricing Estimate value at t = 0 with payoff X at maturity T via MC methods: v = E [e rt X] Haijun Li Math 416/516: Stochastic Simulation Week / 28
26 Typical Example: Discrete-time Asian Options Consider the payoff v = (X K ) +, where X = ( m j=1 X t j )/m for a fixed set of monitoring dates 0 = t 0 < t 1 < < t m = T. Goal: To find the value v = E(e rt (X K ) + ) where X tj = X tj 1 e (r 1 2 σ2 )(t j t j 1 )+σ t j t j 1 Z j, 1 j m. Typical Monte Carlo Algorithm for i = 1,..., n for j = 1,..., m generate the standard normal Z ij set X i (j) = X i (j 1)e (r 1 2 σ2 )(t j t j 1 )+σ t j t j 1 Z ij set X i = (X i (1) + + X i (m))/m set C i = e rt (X i K ) + set Ĉn = (C C n )/n. Haijun Li Math 416/516: Stochastic Simulation Week / 28
27 Efficiency of Simulation Estimators Ĉ n from above illustrative example is unbiased and asymptotically normal. More precisely, let s denote our computational budget, and τ denote the computational time needed for C i, then s[ Ĉ s/τ C] d N(0, σ 2 C τ), as s. In comparing unbiased estimators, we should prefer the one for which σc 2 τ is smallest. For finite but at least moderately large n, we can supplement the point estimate Ĉn with a (1 α)100% confidence interval Ĉ n ± t α/2,n 1 s C n, where s C is the sample standard deviation, and t α/2,n 1 is the upper 100(α/2)th percentage point of a t distribution with n 1 degrees of freedom. Haijun Li Math 416/516: Stochastic Simulation Week / 28
28 Be Aware of Your Estimation Biases! Bias frequently occurs in estimation via MC methods. For example, the bias can arise due to the following errors. 1 Model discretization error: For many models, exact sampling of the continuous-time dynamics is infeasible, some discretization approximation has to be used, resulting a bias. 2 Payoff discretization error: Discretization has to be used for the payoffs that are functionals of the underlying asset processes. 3 Nonlinear functions of means: In a compound option, the price of the first option depends on the price of the second option..., but these prices can only be estimated, resulting a bias. Haijun Li Math 416/516: Stochastic Simulation Week / 28
AMH4 - 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 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 informationA No-Arbitrage Theorem for Uncertain Stock Model
Fuzzy Optim Decis Making manuscript No (will be inserted by the editor) A No-Arbitrage Theorem for Uncertain Stock Model Kai Yao Received: date / Accepted: date Abstract Stock model is used to describe
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 informationNumerical schemes for SDEs
Lecture 5 Numerical schemes for SDEs Lecture Notes by Jan Palczewski Computational Finance p. 1 A Stochastic Differential Equation (SDE) is an object of the following type dx t = a(t,x t )dt + b(t,x t
More informationA new approach for scenario generation in risk management
A new approach for scenario generation in risk management Josef Teichmann TU Wien Vienna, March 2009 Scenario generators Scenarios of risk factors are needed for the daily risk analysis (1D and 10D ahead)
More information1.1 Basic Financial Derivatives: Forward Contracts and Options
Chapter 1 Preliminaries 1.1 Basic Financial Derivatives: Forward Contracts and Options A derivative is a financial instrument whose value depends on the values of other, more basic underlying variables
More information2.1 Mathematical Basis: Risk-Neutral Pricing
Chapter Monte-Carlo Simulation.1 Mathematical Basis: Risk-Neutral Pricing Suppose that F T is the payoff at T for a European-type derivative f. Then the price at times t before T is given by f t = e r(t
More 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 informationFrom Discrete Time to Continuous Time Modeling
From Discrete Time to Continuous Time Modeling Prof. S. Jaimungal, Department of Statistics, University of Toronto 2004 Arrow-Debreu Securities 2004 Prof. S. Jaimungal 2 Consider a simple one-period economy
More informationMODELLING 1-MONTH EURIBOR INTEREST RATE BY USING DIFFERENTIAL EQUATIONS WITH UNCERTAINTY
Applied Mathematical and Computational Sciences Volume 7, Issue 3, 015, Pages 37-50 015 Mili Publications MODELLING 1-MONTH EURIBOR INTEREST RATE BY USING DIFFERENTIAL EQUATIONS WITH UNCERTAINTY J. C.
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 informationOptimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing
Optimal Search for Parameters in Monte Carlo Simulation for Derivative Pricing Prof. Chuan-Ju Wang Department of Computer Science University of Taipei Joint work with Prof. Ming-Yang Kao March 28, 2014
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 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 informationAn Introduction to Stochastic Calculus
An Introduction to Stochastic Calculus Haijun Li lih@math.wsu.edu Department of Mathematics Washington State University Week 5 Haijun Li An Introduction to Stochastic Calculus Week 5 1 / 20 Outline 1 Martingales
More informationStochastic Differential Equations in Finance and Monte Carlo Simulations
Stochastic Differential Equations in Finance and Department of Statistics and Modelling Science University of Strathclyde Glasgow, G1 1XH China 2009 Outline Stochastic Modelling in Asset Prices 1 Stochastic
More informationMonte Carlo Methods for Uncertainty Quantification
Monte Carlo Methods for Uncertainty Quantification Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Mike Giles (Oxford) Monte Carlo methods 2 1 / 24 Lecture outline
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 informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 19 11/20/2013. Applications of Ito calculus to finance
MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.7J Fall 213 Lecture 19 11/2/213 Applications of Ito calculus to finance Content. 1. Trading strategies 2. Black-Scholes option pricing formula 1 Security
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 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 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 informationOption Pricing Models for European Options
Chapter 2 Option Pricing Models for European Options 2.1 Continuous-time Model: Black-Scholes Model 2.1.1 Black-Scholes Assumptions We list the assumptions that we make for most of this notes. 1. The underlying
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 informationNo-arbitrage theorem for multi-factor uncertain stock model with floating interest rate
Fuzzy Optim Decis Making 217 16:221 234 DOI 117/s17-16-9246-8 No-arbitrage theorem for multi-factor uncertain stock model with floating interest rate Xiaoyu Ji 1 Hua Ke 2 Published online: 17 May 216 Springer
More informationMonte Carlo Methods in Financial Engineering
Paul Glassennan Monte Carlo Methods in Financial Engineering With 99 Figures
More informationSTOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL
STOCHASTIC CALCULUS AND BLACK-SCHOLES MODEL YOUNGGEUN YOO Abstract. Ito s lemma is often used in Ito calculus to find the differentials of a stochastic process that depends on time. This paper will introduce
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 informationRandomness and Fractals
Randomness and Fractals Why do so many physicists become traders? Gregory F. Lawler Department of Mathematics Department of Statistics University of Chicago September 25, 2011 1 / 24 Mathematics and the
More informationMulti-Asset Options. A Numerical Study VILHELM NIKLASSON FRIDA TIVEDAL. Master s thesis in Engineering Mathematics and Computational Science
Multi-Asset Options A Numerical Study Master s thesis in Engineering Mathematics and Computational Science VILHELM NIKLASSON FRIDA TIVEDAL Department of Mathematical Sciences Chalmers University of Technology
More informationAn Introduction to Stochastic Calculus
An Introduction to Stochastic Calculus Haijun Li lih@math.wsu.edu Department of Mathematics Washington State University Week 2-3 Haijun Li An Introduction to Stochastic Calculus Week 2-3 1 / 24 Outline
More informationMarket interest-rate models
Market interest-rate models Marco Marchioro www.marchioro.org November 24 th, 2012 Market interest-rate models 1 Lecture Summary No-arbitrage models Detailed example: Hull-White Monte Carlo simulations
More informationTHE USE OF NUMERAIRES IN MULTI-DIMENSIONAL BLACK- SCHOLES PARTIAL DIFFERENTIAL EQUATIONS. Hyong-chol O *, Yong-hwa Ro **, Ning Wan*** 1.
THE USE OF NUMERAIRES IN MULTI-DIMENSIONAL BLACK- SCHOLES PARTIAL DIFFERENTIAL EQUATIONS Hyong-chol O *, Yong-hwa Ro **, Ning Wan*** Abstract The change of numeraire gives very important computational
More informationKing s College London
King s College London University Of London This paper is part of an examination of the College counting towards the award of a degree. Examinations are governed by the College Regulations under the authority
More informationAdvanced Stochastic Processes.
Advanced Stochastic Processes. David Gamarnik LECTURE 16 Applications of Ito calculus to finance Lecture outline Trading strategies Black Scholes option pricing formula 16.1. Security price processes,
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 informationOption Pricing Formula for Fuzzy Financial Market
Journal of Uncertain Systems Vol.2, No., pp.7-2, 28 Online at: www.jus.org.uk Option Pricing Formula for Fuzzy Financial Market Zhongfeng Qin, Xiang Li Department of Mathematical Sciences Tsinghua University,
More informationStochastic Modelling in Finance
in Finance Department of Mathematics and Statistics University of Strathclyde Glasgow, G1 1XH April 2010 Outline and Probability 1 and Probability 2 Linear modelling Nonlinear modelling 3 The Black Scholes
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 informationRisk, Return, and Ross Recovery
Risk, Return, and Ross Recovery Peter Carr and Jiming Yu Courant Institute, New York University September 13, 2012 Carr/Yu (NYU Courant) Risk, Return, and Ross Recovery September 13, 2012 1 / 30 P, Q,
More informationDr. Maddah ENMG 625 Financial Eng g II 10/16/06
Dr. Maddah ENMG 65 Financial Eng g II 10/16/06 Chapter 11 Models of Asset Dynamics () Random Walk A random process, z, is an additive process defined over times t 0, t 1,, t k, t k+1,, such that z( t )
More informationAdvanced Topics in Derivative Pricing Models. Topic 4 - Variance products and volatility derivatives
Advanced Topics in Derivative Pricing Models Topic 4 - Variance products and volatility derivatives 4.1 Volatility trading and replication of variance swaps 4.2 Volatility swaps 4.3 Pricing of discrete
More informationKing s College London
King s College London University Of London This paper is part of an examination of the College counting towards the award of a degree. Examinations are governed by the College Regulations under the authority
More information- 1 - **** d(lns) = (µ (1/2)σ 2 )dt + σdw t
- 1 - **** These answers indicate the solutions to the 2014 exam questions. Obviously you should plot graphs where I have simply described the key features. It is important when plotting graphs to label
More information1 The continuous time limit
Derivative Securities, Courant Institute, Fall 2008 http://www.math.nyu.edu/faculty/goodman/teaching/derivsec08/index.html Jonathan Goodman and Keith Lewis Supplementary notes and comments, Section 3 1
More 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 informationMFE/3F Questions Answer Key
MFE/3F Questions Download free full solutions from www.actuarialbrew.com, or purchase a hard copy from www.actexmadriver.com, or www.actuarialbookstore.com. Chapter 1 Put-Call Parity and Replication 1.01
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 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 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 Runge Kutta Methods with the Constant Elasticity of Variance (CEV) Diffusion Model for Pricing Option
Int. Journal of Math. Analysis, Vol. 8, 2014, no. 18, 849-856 HIKARI Ltd, www.m-hikari.com http://dx.doi.org/10.12988/ijma.2014.4381 Stochastic Runge Kutta Methods with the Constant Elasticity of Variance
More informationMFE/3F Questions Answer Key
MFE/3F Questions Download free full solutions from www.actuarialbrew.com, or purchase a hard copy from www.actexmadriver.com, or www.actuarialbookstore.com. Chapter 1 Put-Call Parity and Replication 1.01
More informationAnalytical formulas for local volatility model with stochastic. Mohammed Miri
Analytical formulas for local volatility model with stochastic rates Mohammed Miri Joint work with Eric Benhamou (Pricing Partners) and Emmanuel Gobet (Ecole Polytechnique Modeling and Managing Financial
More informationAmerican Option Pricing Formula for Uncertain Financial Market
American Option Pricing Formula for Uncertain Financial Market Xiaowei Chen Uncertainty Theory Laboratory, Department of Mathematical Sciences Tsinghua University, Beijing 184, China chenxw7@mailstsinghuaeducn
More informationDrunken Birds, Brownian Motion, and Other Random Fun
Drunken Birds, Brownian Motion, and Other Random Fun Michael Perlmutter Department of Mathematics Purdue University 1 M. Perlmutter(Purdue) Brownian Motion and Martingales Outline Review of Basic Probability
More informationAn Efficient Numerical Scheme for Simulation of Mean-reverting Square-root Diffusions
Journal of Numerical Mathematics and Stochastics,1 (1) : 45-55, 2009 http://www.jnmas.org/jnmas1-5.pdf JNM@S Euclidean Press, LLC Online: ISSN 2151-2302 An Efficient Numerical Scheme for Simulation of
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 informationProbability in Options Pricing
Probability in Options Pricing Mark Cohen and Luke Skon Kenyon College cohenmj@kenyon.edu December 14, 2012 Mark Cohen and Luke Skon (Kenyon college) Probability Presentation December 14, 2012 1 / 16 What
More informationModule 10:Application of stochastic processes in areas like finance Lecture 36:Black-Scholes Model. Stochastic Differential Equation.
Stochastic Differential Equation Consider. Moreover partition the interval into and define, where. Now by Rieman Integral we know that, where. Moreover. Using the fundamentals mentioned above we can easily
More informationMonte Carlo Methods in Structuring and Derivatives Pricing
Monte Carlo Methods in Structuring and Derivatives Pricing Prof. Manuela Pedio (guest) 20263 Advanced Tools for Risk Management and Pricing Spring 2017 Outline and objectives The basic Monte Carlo algorithm
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 informationComputing Greeks with Multilevel Monte Carlo Methods using Importance Sampling
Computing Greeks with Multilevel Monte Carlo Methods using Importance Sampling Supervisor - Dr Lukas Szpruch Candidate Number - 605148 Dissertation for MSc Mathematical & Computational Finance Trinity
More informationSlides for DN2281, KTH 1
Slides for DN2281, KTH 1 January 28, 2014 1 Based on the lecture notes Stochastic and Partial Differential Equations with Adapted Numerics, by J. Carlsson, K.-S. Moon, A. Szepessy, R. Tempone, G. Zouraris.
More informationLecture 11: Ito Calculus. Tuesday, October 23, 12
Lecture 11: Ito Calculus Continuous time models We start with the model from Chapter 3 log S j log S j 1 = µ t + p tz j Sum it over j: log S N log S 0 = NX µ t + NX p tzj j=1 j=1 Can we take the limit
More informationComputational Finance. Computational Finance p. 1
Computational Finance Computational Finance p. 1 Outline Binomial model: option pricing and optimal investment Monte Carlo techniques for pricing of options pricing of non-standard options improving accuracy
More informationForeign Exchange Derivative Pricing with Stochastic Correlation
Journal of Mathematical Finance, 06, 6, 887 899 http://www.scirp.org/journal/jmf ISSN Online: 6 44 ISSN Print: 6 434 Foreign Exchange Derivative Pricing with Stochastic Correlation Topilista Nabirye, Philip
More informationAn Overview of Volatility Derivatives and Recent Developments
An Overview of Volatility Derivatives and Recent Developments September 17th, 2013 Zhenyu Cui Math Club Colloquium Department of Mathematics Brooklyn College, CUNY Math Club Colloquium Volatility Derivatives
More informationMultilevel quasi-monte Carlo path simulation
Multilevel quasi-monte Carlo path simulation Michael B. Giles and Ben J. Waterhouse Lluís Antoni Jiménez Rugama January 22, 2014 Index 1 Introduction to MLMC Stochastic model Multilevel Monte Carlo Milstein
More informationValuation of Asian Option. Qi An Jingjing Guo
Valuation of Asian Option Qi An Jingjing Guo CONTENT Asian option Pricing Monte Carlo simulation Conclusion ASIAN OPTION Definition of Asian option always emphasizes the gist that the payoff depends on
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 informationNumerical Simulation of Stochastic Differential Equations: Lecture 1, Part 2. Integration For deterministic h : R R,
Numerical Simulation of Stochastic Differential Equations: Lecture, Part Des Higham Department of Mathematics University of Strathclyde Lecture, part : SDEs Ito stochastic integrals Ito SDEs Examples of
More informationMonte Carlo Methods. Prof. Mike Giles. Oxford University Mathematical Institute. Lecture 1 p. 1.
Monte Carlo Methods Prof. Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute Lecture 1 p. 1 Geometric Brownian Motion In the case of Geometric Brownian Motion ds t = rs t dt+σs
More informationConditional Full Support and No Arbitrage
Gen. Math. Notes, Vol. 32, No. 2, February 216, pp.54-64 ISSN 2219-7184; Copyright c ICSRS Publication, 216 www.i-csrs.org Available free online at http://www.geman.in Conditional Full Support and No Arbitrage
More informationContinuous Time Finance. Tomas Björk
Continuous Time Finance Tomas Björk 1 II Stochastic Calculus Tomas Björk 2 Typical Setup Take as given the market price process, S(t), of some underlying asset. S(t) = price, at t, per unit of underlying
More informationOptimal stopping problems for a Brownian motion with a disorder on a finite interval
Optimal stopping problems for a Brownian motion with a disorder on a finite interval A. N. Shiryaev M. V. Zhitlukhin arxiv:1212.379v1 [math.st] 15 Dec 212 December 18, 212 Abstract We consider optimal
More informationMath 623 (IOE 623), Winter 2008: Final exam
Math 623 (IOE 623), Winter 2008: Final exam Name: Student ID: This is a closed book exam. You may bring up to ten one sided A4 pages of notes to the exam. You may also use a calculator but not its memory
More informationCalibration of Interest Rates
WDS'12 Proceedings of Contributed Papers, Part I, 25 30, 2012. ISBN 978-80-7378-224-5 MATFYZPRESS Calibration of Interest Rates J. Černý Charles University, Faculty of Mathematics and Physics, Prague,
More informationRohini Kumar. Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque)
Small time asymptotics for fast mean-reverting stochastic volatility models Statistics and Applied Probability, UCSB (Joint work with J. Feng and J.-P. Fouque) March 11, 2011 Frontier Probability Days,
More informationModeling via Stochastic Processes in Finance
Modeling via Stochastic Processes in Finance Dimbinirina Ramarimbahoaka Department of Mathematics and Statistics University of Calgary AMAT 621 - Fall 2012 October 15, 2012 Question: What are appropriate
More informationFINANCIAL OPTION ANALYSIS HANDOUTS
FINANCIAL OPTION ANALYSIS HANDOUTS 1 2 FAIR PRICING There is a market for an object called S. The prevailing price today is S 0 = 100. At this price the object S can be bought or sold by anyone for any
More informationComputational Finance Improving Monte Carlo
Computational Finance Improving Monte Carlo School of Mathematics 2018 Monte Carlo so far... Simple to program and to understand Convergence is slow, extrapolation impossible. Forward looking method ideal
More informationIntroduction to Affine Processes. Applications to Mathematical Finance
and Its Applications to Mathematical Finance Department of Mathematical Science, KAIST Workshop for Young Mathematicians in Korea, 2010 Outline Motivation 1 Motivation 2 Preliminary : Stochastic Calculus
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 informationImplementing the HJM model by Monte Carlo Simulation
Implementing the HJM model by Monte Carlo Simulation A CQF Project - 2010 June Cohort Bob Flagg Email: bob@calcworks.net January 14, 2011 Abstract We discuss an implementation of the Heath-Jarrow-Morton
More informationOPTIMAL PORTFOLIO CONTROL WITH TRADING STRATEGIES OF FINITE
Proceedings of the 44th IEEE Conference on Decision and Control, and the European Control Conference 005 Seville, Spain, December 1-15, 005 WeA11.6 OPTIMAL PORTFOLIO CONTROL WITH TRADING STRATEGIES OF
More informationDynamic Hedging and PDE Valuation
Dynamic Hedging and PDE Valuation Dynamic Hedging and PDE Valuation 1/ 36 Introduction Asset prices are modeled as following di usion processes, permitting the possibility of continuous trading. This environment
More informationNumerical Methods for Stochastic Differential Equations with Applications to Finance
Numerical Methods for Stochastic Differential Equations with Applications to Finance Matilde Lopes Rosa Instituto Superior Técnico University of Lisbon, Portugal May 2016 Abstract The pricing of financial
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 informationCredit Risk : Firm Value Model
Credit Risk : Firm Value Model Prof. Dr. Svetlozar Rachev Institute for Statistics and Mathematical Economics University of Karlsruhe and Karlsruhe Institute of Technology (KIT) Prof. Dr. Svetlozar Rachev
More informationEuropean call option with inflation-linked strike
Mathematical Statistics Stockholm University European call option with inflation-linked strike Ola Hammarlid Research Report 2010:2 ISSN 1650-0377 Postal address: Mathematical Statistics Dept. of Mathematics
More informationTEST OF BOUNDED LOG-NORMAL PROCESS FOR OPTIONS PRICING
TEST OF BOUNDED LOG-NORMAL PROCESS FOR OPTIONS PRICING Semih Yön 1, Cafer Erhan Bozdağ 2 1,2 Department of Industrial Engineering, Istanbul Technical University, Macka Besiktas, 34367 Turkey Abstract.
More informationPractical example of an Economic Scenario Generator
Practical example of an Economic Scenario Generator Martin Schenk Actuarial & Insurance Solutions SAV 7 March 2014 Agenda Introduction Deterministic vs. stochastic approach Mathematical model Application
More informationFast and accurate pricing of discretely monitored barrier options by numerical path integration
Comput Econ (27 3:143 151 DOI 1.17/s1614-7-991-5 Fast and accurate pricing of discretely monitored barrier options by numerical path integration Christian Skaug Arvid Naess Received: 23 December 25 / Accepted:
More informationRisk Neutral Measures
CHPTER 4 Risk Neutral Measures Our aim in this section is to show how risk neutral measures can be used to price derivative securities. The key advantage is that under a risk neutral measure the discounted
More informationStochastic Processes and Advanced Mathematical Finance. Stochastic Differential Equations and the Euler-Maruyama Method
Steven R. Dunbar Department of Mathematics 203 Avery Hall University of Nebraska-Lincoln Lincoln, NE 68588-0130 http://www.math.unl.edu Voice: 402-472-3731 Fax: 402-472-8466 Stochastic Processes and Advanced
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 informationStochastic Volatility Effects on Defaultable Bonds
Stochastic Volatility Effects on Defaultable Bonds Jean-Pierre Fouque Ronnie Sircar Knut Sølna December 24; revised October 24, 25 Abstract We study the effect of introducing stochastic volatility in the
More informationSubject CT8 Financial Economics Core Technical Syllabus
Subject CT8 Financial Economics Core Technical Syllabus for the 2018 exams 1 June 2017 Aim The aim of the Financial Economics subject is to develop the necessary skills to construct asset liability models
More informationVariance Reduction for Monte Carlo Simulation in a Stochastic Volatility Environment
Variance Reduction for Monte Carlo Simulation in a Stochastic Volatility Environment Jean-Pierre Fouque Tracey Andrew Tullie December 11, 21 Abstract We propose a variance reduction method for Monte Carlo
More information