Stochastic Processes and Advanced Mathematical Finance. Stochastic Differential Equations and the Euler-Maruyama Method
|
|
- Jocelyn Paul
- 6 years ago
- Views:
Transcription
1 Steven R. Dunbar Department of Mathematics 203 Avery Hall University of Nebraska-Lincoln Lincoln, NE Voice: Fax: Stochastic Processes and Advanced Mathematical Finance Stochastic Differential Equations and the Euler-Maruyama Method Rating Mathematically Mature: proofs. may contain mathematics beyond calculus with 1
2 Section Starter Question Explain how to use a slope-field diagram to solve the ordinary differential equation dx dt = x. How would you turn that process into an algorithm to numerically compute an approximate solution without a diagram? Key Concepts 1. We can numerically simulate the solution to stochastic differential equations with an analog to Euler s method, called the Euler-Maruyama (EM) method. Vocabulary 1. A stochastic differential equation is a mathematical equation relating a stochastic process to its local deterministic and random components. The goal is to extend the relation to find the stochastic process. Under mild conditions on the relationship, and with a specifying initial condition, solutions of stochastic differential equations exist and are unique. 2. The Euler-Maruyama (EM) method is a numerical method for simulating the solutions of a stochastic differential equation based on 2
3 the definition of the Itô stochastic integral: Given dx(t) = G(X(t)) dt +H(X(t)) dw (t), X(t 0 ) = X 0, and a step size dt, we approximate and simulate with X j = X j 1 + G(X j 1 ) dt +H(X j 1 )(W (t j 1 + dt) W (t j 1 )) 3. Extensions and variants of standard Brownian motion defined through stochastic differential equations are Brownian motion with drift, scaled Brownian motion, and geometric Brownian motion. Mathematical Ideas Stochastic Differential Equations: Symbolically The straight line segment is the building block of differential calculus. The basic idea behind differential calculus is that differentiable functions, no matter how difficult their global behavior, are locally approximated by straight line segments. In particular, this is the idea behind Euler s method for approximating differentiable functions defined by differential equations. We know that rescaling ( zooming in on) Brownian motion does not produce a straight line, it produces another image of Brownian motion. This self-similarity is ideal for an infinitesimal building block, for instance, we could build global Brownian motion out of lots of local chunks of Brownian motion. This suggests we could build other stochastic processes out of suitably scaled Brownian motion. In addition, if we include straight line segments we can overlay the behavior of differentiable functions onto the stochastic processes as well. Thus, straight line segments and chunks of Brownian motion are the building blocks of stochastic calculus. With stochastic differential calculus, we can build new stochastic processes. We do this by specifying how to build the new stochastic processes 3
4 locally from our base deterministic function, the straight line and our base stochastic process, standard Brownian motion. We write the local change in value of the stochastic process over a time interval of (infinitesimal) length dt as dx = G(X(t)) dt +H(X(t)) dw (t), X(t 0 ) = X 0. (1) Note that we are not allowed to write dx dt = G(X(t)) + H(X(t)) dw dt, X(t 0) = X 0 since standard Brownian motion is nowhere differentiable with probability 1. Actually, the informal stochastic differential equation (1) is a compact way of writing a rigorously defined, equivalent implicit Itô integral equation. Since we do not have the required rigor, we will approach the stochastic differential equation intuitively. The stochastic differential equation says the initial point (t 0, X 0 ) is specified, perhaps with X 0 a random variable with a given distribution. A deterministic component at each point has a slope determined through G at that point. In addition, some random perturbation affects the evolution of the process. The random perturbation is normally distributed with mean 0. The variance of the random perturbation is (H(X(t))) 2 at (t, X(t)). This is a simple expression of a Stochastic Differential Equation (SDE) which determines a stochastic process, just as an Ordinary Differential Equation (ODE) determines a differentiable function. We extend the process with the incremental change information and repeat. This is an expression in words of the Euler-Maruyama method for numerically simulating the stochastic differential expression. Example. A very simple stochastic differential equation is dx = r dt + dw, X(0) = b with r a constant. Take a deterministic initial condition to be X(0) = b. The new process is the stochastic extension of the differential equation expression of a straight line. The new stochastic process X is drifting or trending at constant rate r with a random variation due to Brownian motion perturbations around that trend. We will later show explicitly that the solution of this SDE is X(t) = b + rt + W (t) although it is seems intuitively clear that this should be the process. We will call this Brownian motion with drift. 4
5 Example. Another very simple stochastic differential equation is dx = σ dw, X(0) = b This stochastic differential equation says that the process is evolving as a multiple of standard Brownian motion. The solution may be easily guessed as X(t) = σw (t) which has variance σ 2 t on increments of length t. Sometimes the new process is called Brownian motion (in contrast to standard Brownian motion which has variance t on increments of length t). We combine the previous two examples to consider dx = r dt +σ dw, X(0) = b which has solution X(t) = b+rt+σw (t), a multiple of Brownian motion with drift r started at b. Sometimes this extension of standard Brownian motion is called Brownian motion. Some authors consider this process directly instead of the more special case we considered in the previous chapter. Example. The next simplest and first non-trivial differential equation is dx = X dw. Here the differential equation says that process is evolving like Brownian motion with a variance which is the square of the process value. When the process is small, the variance is small, when the process is large, the variance is large. Expressing the stochastic differential equation as dx /X = dw we may say that the relative change acts like standard Brownian motion. The resulting stochastic process is called geometric Brownian motion and it will figure extensively later as a model of security prices. Example. The next simplest differential equation is dx = rx dt +σx dw, X(0) = b. Here the stochastic differential equation says that the growth of the process at a point is proportional to the process value, with a random perturbation proportional to the process value. Again looking ahead, we could write the differential equation as dx /X = r dt +σ dw and interpret it to say the relative rate of increase is proportional to the time observed together with a random perturbation like a Brownian increment corresponding to the length of time. We will show later that the analytic expression for the stochastic process defined by this SDE is b exp((r 1 2 σ2 )t + σw (t)). 5
6 Stochastic Differential Equations: Numerically The sample path that the Euler-Maruyama method produces numerically is the analog of using the Euler method. The formula for the Euler-Maruyama (EM) method is based on the definition of the Itô stochastic integral: X j = X j 1 + G(X j 1 ) dt +H(X j 1 )(W (t j 1 + dt) W (t j 1 )), t j = t j 1 + dt. Note that the initial conditions X 0 and t 0 set the starting point. In this text, we use coin-flipping sequences of an appropriate length scaled to create an approximation to W (t) just as in the section Approximation to Brownian Motion. The coin-flipping sequences emphasize the discrete nature of the simulations with an easily constructed random process. This is consistent with the approach of this text which always uses coin-flipping sequences to create random processes. Note that since the increments W (t j 1 + dt) W (t j 1 ) are independent and identically distributed, we will use independent coin-flip sequences to generate the approximation of the increments. The EM method could use independent normal random variates directly to obtain the increments W (t j 1 + dt) W (t j 1 ). Using independent normal random variates directly would be easier and more efficient. The exercises modify the example scripts to use independent normal random variates directly. Then dw = W (t j 1 + dt) W (t j 1 ) = W ( dt) ŴN( dt) = Ŵ (N dt) N = (N dt) dtŵ. N dt The first equality above is the definition of an increment, the second equality means the random variables W (t j 1 + dt) W (t j 1 ) and W ( dt) have the same distribution because of the definition of standard Brownian motion which specifies that increments with equal length are normally distributed with variance equal to the increment length. The approximate equality occurs because of the approximation of Brownian motion by coin-flipping sequences. We generate the approximations using a random number generator, but we could as well use actual coin-flipping. In Table 1 the generation of the sequences is not recorded, only the summed and scaled (independently 6
7 j t j X j 2X j dt dw X j dw 2X j dt + X j + X j dw 2X j dt + X j dw Figure 1: Simulation with the Euler-Maruyama method of a process defined by a stochastic differential equation. sampled) outcomes. For convenience, take dt = 1/10, N = 100, so we need Ŵ (100 (1/10))/ 100 = T 10 /10. Then to obtain the entries in the column labeled dw in the table we flip a coin 10 times and record T 10 /10. Take r = 2, b = 1, and σ = 1, so we simulate the solution of dx = 2X dt +X dw, X(0) = 1. A computer program can produce such a table with the step size made much smaller, presumably resulting in better approximation properties. In fact, it is possible to consider kinds of convergence for the EM method comparable to the Strong Law of Large Numbers and the Weak Law of Large Numbers. See the Problems for examples. Discussion The numerical approximation procedure using coin-flipping makes it clear that the Euler-Maruyama method generates a random process. The value of the process depends on the time value and the coin-flip sequence. Each generation of an approximation will be different because the coin-flip sequence is different. The Euler-Maruyama method generates a stochastic process 7
8 path approximation. To derive distributions and statistics about the process requires generating multiple paths, see the Problems for examples. This shows that stochastic differential equations provide a way to define new stochastic processes. This is analogous to the notion that ordinary differential equations define new functions to study and use. In fact, one approach to developing calculus and the analysis of functions is to start with differential equations, use the Euler method to define approximations of solutions, and then to develop a theory to handle the passage to continuous variables. This approach is especially useful for a mathematical modeling viewpoint since the model often uses differential equations. This text follows the approach of starting with stochastic differential equations to describe a situation and numerically defining new stochastic processes to model the situation. At certain points, we appeal to more rigorous mathematical theory to justify the modeling and approximation. One important justification asserts that if we write a stochastic differential equation, then solutions exist and the stochastic differential equation always yields the same process under equivalent conditions. The Existence-Uniqueness Theorem shows that under reasonable modeling conditions stochastic differential equations do indeed satisfy this requirement. Theorem 1 (Existence-Uniqueness). For the stochastic differential equation assume dx = G(t, X(t)) dt +H(t, X(t)) dw (t), X(t 0 ) = X 0 1. Both G(t, x) and H(t, x) are continuous on (t, x) [t 0, T ] R. 2. The coefficient functions G and H satisfy a Lipschitz condition: G(t, x) G(t, y) + H(t, x) H(t, y) K x y. 3. The coefficient functions G and H satisfy a growth condition in the second variable for all t [t 0, T ] and x R. G(t, x) 2 + H(t, x) 2 K(1 + x 2 ) 8
9 Then the stochastic differential equation has a strong solution on [t 0, T ] that is continuous with probability 1 and sup E [ X 2 (t) ] < t [t 0,T ] and for each given Wiener process W (t), the corresponding strong solutions are pathwise unique which means that if X and Y are two strong solutions, then [ ] P sup X(t) Y (t) = 0 t [t 0,T ] = 1. See [4] for a precise definition of strong solution but essentially it means that for each given Wiener process W (t) we can generate a solution to the SDE. Note that the coefficient functions here are two-variable functions of both time t and location x, which is more general than the functions considered in equation (1). The restrictions on the functions G(t, x) and H(t, x), especially the continuity condition, can be considerably relaxed and the theorem will still remain true. Sources This section is adapted from: An Algorithmic Introduction to the Numerical Simulation of Stochastic Differential Equations, by Desmond J. Higham, in SIAM Review, Vol. 43, No. 3, pp , 2001 and Financial Calculus: An introduction to derivative pricing by M. Baxter, and A. Rennie, Cambridge University Press, 1996, pages The Existence-Uniqueness Theorem is adapted from An Introduction to Stochastic Processes with Applications to Biology, by L. J. S. Allen, Pearson Prentice-Hall, 2003, pages and Numerical Solution of Stochastic Differential Equations, by Peter Kloeden and Eckhard Platen, Springer Verlag, 1992, pages
10 Algorithms, Scripts, Simulations Algorithm The scripts apply the EM method to simulate the solution of dx = rx dt +σx dw, X(0) = b. The parameters N for the number of steps in the EM method, T for the ending time, and stochastic differential equation parameters r and s are set. Find the time step and initialize the arrays holding the time steps and the solution simulation. Using M = 30N create a piecewise linear function ŴM(t) using the approximation scripts in Approximation to Brownian Motion. The parameter 30 is chosen according to the rule of thumb in the DeMoivre-Laplace Central Limit Theorem, Central Limit Theorem. Then each time increment will have 30 coin flips, sufficient according to the rule of thumb to guarantee a scaled sum which is appropriately normally distributed. Then loop over the number of steps using the EM algorithm X j = X j 1 + G(X j 1 ) dt +H(X j 1 )(W (t j 1 + dt) W (t j 1 )) and plot the resulting simulation. Scripts Geogebra GeoGebra R R script for stochasticsdes.r 1 r <- -1 # growth / decay rate 2 sigma <- 0.5 # relative standard deviation 3 b <- 3 # initial condition 4 5 M <- 100 # number of steps for EM method to take 6 T <- 1 # maximum time 7 h <- T/ M # time step 8 t <- seq ( length =M+1, from =0, by= h) # t is the vector [0 1h 2h 3h... Nh] 10
11 9 X <- array (0, c( M +1) ) # place to store locations N <- 30*( M +1) # number of steps for the Brownian Motion approx p < S <- array (0, c(n +1) ) 15 rw <- cumsum ( 2 *( runif (N) <= p) -1 ) 16 S [2:( N +1) ] <- rw WcaretN <- function ( z) { 19 Delta <- T/N # add 1 since arrays are 1- based 22 prior = floor (z/ Delta ) subsequent = ceiling ( z/ Delta ) retval <- sqrt ( Delta )*(S[ prior ] + ((z/ Delta +1) - prior )*(S[ subsequent ] - S[ prior ])) 26 } X [1] <- b 29 for (i in 1:M) { 30 X[i +1] <- X[i]+r*X[i]*h+ sigma *X[i]*( WcaretN (t[i]+h)- WcaretN (t[i])) 31 } plot (t,x,"l", xlim =c(0, T), ylim =c(x[1] - exp ( abs (r)*t +1), X [1]+ exp ( abs (r)*t +1) )) 34 title ( main = paste ("r = ", r, " sigma = ", sigma, " steps =", M)) 35 Octave Octave script for stochasticdes.m 1 r = -1; # growth / decay rate 2 sigma = 0.5; # \ sigma 3 b = 3; # initial value 4 5 M =100; # number of steps for EM method to take 6 global T =1; # maximum time 7 h=t/m; # time step 8 t =(0: h:t); # t is the vector [0 1h 2h 3h... Nh] 11
12 9 X= zeros ( size (t));# prepare place to store locations global N = 30*( M +1) ; # number of steps for the Brownian 12 # Motion approximation 13 global S; 14 p = 1/ 2; 15 S = zeros (N+1, 1); 16 S (2: N +1) = cumsum ( 2 * ( rand (N,1) <=p) - 1); function retval = WcaretN ( z) 19 global N; 20 global T; 21 global S; 22 Delta = T/N; # add 1 since arrays are 1- based 25 prior = floor (z/ Delta ) + 1; 26 subsequent = ceil ( z/ Delta ) + 1; retval = sqrt ( Delta )*(S( prior ) + ((z/ Delta +1) - prior ).*( S( subsequent )-S( prior ))); endfunction X (1) = b; # initial height at t = 0 33 for i =1: M # start taking steps 34 X(i +1) =X(i)+r*X(i)*h+ sigma *X(i)*( WcaretN (t(i)+h)- WcaretN (t(i)) ); 35 end ; plot ( t, X) # plot more permanently 38 axis ([0 T X (1) -exp ( abs (r)*t +1) X (1) + exp ( abs (r)*t +1) ]); # set axis limits 39 grid on; 40 title (["r = ", num2str (r), ", sigma = ", num2str ( sigma ), ", steps = ", num2str (M)] ); 41 Perl Perl PDL script for stochasticdes.pl 1 # #!/ usr / bin / pdl 2 use PDL :: NiceSlice ; 3 12
13 4 $r = -1; # growth / decay rate 5 $sigma = 0. 5; # standard deviation 6 $b = 3; # initial value 7 8 $M = 100; # number of steps for EM method to take 9 $T = 1; # maximum time 10 $h = $T/ $M; # time step 11 $t = zeros ( $M + 1 ) - > xlinvals ( 0, $T ); # vector of [0, 1h, 2h, 3h... Mh] 12 $X = zeros ( $M + 1 ); $N = 30*( $M +1) ; # number of steps for the Brownian Motion $p = 0. 5; 17 $S = zeros ( $N + 1 ); # the random walk 18 $S ( 1 : $N ).= cumusumover ( 2 * ( random ( $N) <= $p ) - 1 ); # function WcaretN interpolating random walk 21 sub WcaretN { 22 my $x = 23 $Delta = $T / $N; $prior = floor ( $x / $Delta ); 26 $subsequent = ceil ( $x / $Delta ); $retval = 29 sqrt ( $Delta ) 30 * ( $S ( $prior ) 31 + ( ( $x / $Delta ) - $prior ) 32 * ( $S ( $subsequent ) - $S ( $prior ) ) ); 33 } $X (0).= $b; # initial value at t = 0 36 for ( $i =0; $i <= $M -1; $i ++) { # start taking steps 37 $X($i +1).= $X($i)+$r*$X($i)*$h+ $sigma *$X($i)*( WcaretN ($t($i)+$h)-wcaretn ($t($i)) ); 38 } # file output to use with external plotting programming 41 # such as gnuplot, R, octave, etc. 42 # Start gnuplot, then from gnuplot prompt 43 # plot " stochasticdes. dat " with lines 13
14 44 45 open ( F, ">stochasticdes. dat " ) die " cannot write : $! "; 46 foreach $j ( 0.. $M ) { 47 print F $t -> range ( [$j] ), " ", $X -> range ( [$j] ), "\ n"; 48 } 49 close (F); 50 SciPy Scientific Python script for stochasticdes.py 1 import scipy 2 3 r = -1. # growth / decay rate 4 sigma = 0.5 # standard deviation 5 b = 3. # initial value 6 7 M = 100 # number of steps for EM method to take 8 T = 1. # maximum time, note type real 9 h = T/M # time step 10 t = scipy. linspace (0,T,M +1) # vector of [0, 1h, 2h, 3h... Nh] 11 X = scipy. zeros ( M + 1 ) N = 30*( M +1) # number of steps for the Brownian Motion # the random walk 16 p = S = scipy. zeros (N +1) 18 S [1: N +1] = scipy. cumsum ( 2*( scipy. random. random (N) <= p ) - 1 ) def WcaretN (x): 21 Delta = T/ N # T real coerces Delta real 22 prior = scipy. floor (x/ Delta ). astype ( int ) 23 subsequent = scipy. ceil (x/ Delta ). astype ( int ) 24 return scipy. sqrt ( Delta )*(S[ prior ] + (x/ Delta - prior )*(S[ subsequent ] - S[ prior ])) X [0] = b # iniital value at t = 0 27 for i in range (0,M): 14
15 28 X[i +1] = X[i]+r*X[i]*h+ sigma *X[i ]*( WcaretN (t[i]+h)- WcaretN (t[i]) ) # optional file output to use with external plotting programming 31 # such as gnuplot, R, octave, etc. 32 # Start gnuplot, then from gnuplot prompt 33 # plot " stochasticdes. dat " with lines 34 f = open ( stochasticdes. dat, w ) 35 for j in range (0,M +1) : 36 f. write ( str (t[j])+ + str (X[j])+ \n ); f. close () Problems to Work for Understanding 1. Graph an approximation of a multiple of Brownian motion with drift with parameters b = 2, r = 1/2 and σ = 2 on the interval [0, 5] in two ways. (a) Flip a coin 25 times, recording whether it comes up Heads or Tails each time, Scoring Y i = +1 for each Heads and Y i = 1 for each flip, also keep track of the accumulated sum T n = n i=1 T i for i = Using N = 5 compute the rescaled approximation Ŵ 5 (t) = (1/ 5)T 5t at the values t = 0, 1/5, 2/5, 3/5,... 24/5, 5 on [0, 5]. Finally compute and graph the value of X(t) = b + rt + σŵ5(t). (b) Using the same values of Ŵ5(t) as approximations for W (dt) compute the values of the solution of the stochastic differential equation dx = r dt +σ dw, X(0) = b on the interval [0, 5]. 2. Repeat the previous problem with parameters b = 2, r = 1/2 and σ = 2. 15
16 3. Repeat the previous problem with parameters b = 2, r = 1/2 and σ = Modify the scripts to use normally distributed random deivates and then simulate the solution of the stochastic differential equation dx(t) = X(t) dt +2X(t) dx on the interval [0, 1] with initial condition X(0) = 1 and step size t = 1/ Modify the scripts to use normally distributed random deviates and then simulate the solution of the stochastic differential equation dx(t) = tx(t) dt +2X(t) dx on the interval [0, 1] with initial condition X(0) = 1 and step size t = 1/10. Note the difference with the previous problem, now the multiplier of the dt term is a function of time. 6. Modify the scripts to use general functions G(X) and H(x) and then apply it to the SDE for the (a) Ornstein-Uhlenbeck process: dx = θ(µ X(t)) dt +σ dw (t), X(t 0 ) = X 0. (b) Cox-Ingersoll-Ross process: dx = θ(µ X(t)) dt +σ (X(t)) dw (t), X(t 0 ) = X 0. (c) the modified Cox-Ingersoll-Ross process: dx = θx(t) dt +θ (X(t)) dw (t), X(t 0 ) = X Write a program with parameters r, σ, b, T and N (so dt = T/N) that computes and graphs the approximation of the solution of the stochastic differential equation dx(t) = rx(t) dt +σx(t) dx with X(0) = b on the interval [0, T ]. Apply the program to the stochastic differential equation with r = 2, σ = 1, b = 1, and N = 2 6, 2 7, 2 8 on the interval [0, 1]. 16
17 8. Generalize the program from the previous problem to include a parameter M for the number of sample paths computed. Then using this program on the interval [0, 1] with M = 1000, and N = 2 8 compute E [ X n X(1) ], where X(1) = be r 1 2 σ2 +σŵ2 8 (1). 9. Using the program from the previous problem with M = 1000 and N = 2 5, 2 6, 2 7, 2 8, 2 9 compute E [ X N X(1) ], where X(1) = be r 1 2 σ2 +σŵ2 9 (1). Then for the 5 values of N, make a log-log plot of E [ X N X(1) ] on the vertical axis against t = 1/N on the horizontal axis. Using the slope of the resulting best-fit line experimentally determine the order of convergence γ so that E [ X N X(1) ] C( t) γ. Reading Suggestion: References [1] Linda J. S. Allen. An Introduction to Stochastic Processes with Applications to biology. Pearson Prentice-Hall, [2] M. Baxter and A. Rennie. Financial Calculus: An introduction to derivative pricing. Cambridge University Press, HG 6024 A2W554. [3] Desmond J. Higham. An algorithmic introduction to numerical simulation of stochastic differential equations. SIAM Review, 43(3): , [4] P. Kloeden and E. Platen. Numerical Solution of Stochastic Differential Equations, volume 23 of Stochastic Modelling and Applied Probability. Springer,
18 Outside Readings and Links: 1. Matlab program files for Stochastic Differential Equations offers a number of MATLAB routines for stochastic differential equations. I check all the information on each page for correctness and typographical errors. Nevertheless, some errors may occur and I would be grateful if you would alert me to such errors. I make every reasonable effort to present current and accurate information for public use, however I do not guarantee the accuracy or timeliness of information on this website. Your use of the information from this website is strictly voluntary and at your risk. I have checked the links to external sites for usefulness. Links to external websites are provided as a convenience. I do not endorse, control, monitor, or guarantee the information contained in any external website. I don t guarantee that the links are active at all times. Use the links here with the same caution as you would all information on the Internet. This website reflects the thoughts, interests and opinions of its author. They do not explicitly represent official positions or policies of my employer. Information on this website is subject to change without notice. Steve Dunbar s Home Page, to Steve Dunbar, sdunbar1 at unl dot edu Last modified: Processed from L A TEX source on August 2,
Stochastic Processes and Advanced Mathematical Finance. Multiperiod Binomial Tree Models
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 informationStochastic Processes and Advanced Mathematical Finance. Hitting Times and Ruin Probabilities
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 informationStochastic Processes and Advanced Mathematical Finance. A Stochastic Process Model of Cash Management
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 informationStochastic Processes and Advanced Mathematical Finance. Single Period Binomial Models
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 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 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 information23 Stochastic Ordinary Differential Equations with Examples from Finance
23 Stochastic Ordinary Differential Equations with Examples from Finance Scraping Financial Data from the Web The MATLAB/Octave yahoo function below returns daily open, high, low, close, and adjusted close
More informationAMH4 - ADVANCED OPTION PRICING. Contents
AMH4 - ADVANCED OPTION PRICING ANDREW TULLOCH Contents 1. Theory of Option Pricing 2 2. Black-Scholes PDE Method 4 3. Martingale method 4 4. Monte Carlo methods 5 4.1. Method of antithetic variances 5
More 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 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 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 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 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 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 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 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 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 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 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 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 informationBarrier Options Pricing in Uncertain Financial Market
Barrier Options Pricing in Uncertain Financial Market Jianqiang Xu, Jin Peng Institute of Uncertain Systems, Huanggang Normal University, Hubei 438, China College of Mathematics and Science, Shanghai Normal
More informationA Numerical Approach to the Estimation of Search Effort in a Search for a Moving Object
Proceedings of the 1. Conference on Applied Mathematics and Computation Dubrovnik, Croatia, September 13 18, 1999 pp. 129 136 A Numerical Approach to the Estimation of Search Effort in a Search for a Moving
More informationStochastic Calculus - An Introduction
Stochastic Calculus - An Introduction M. Kazim Khan Kent State University. UET, Taxila August 15-16, 17 Outline 1 From R.W. to B.M. B.M. 3 Stochastic Integration 4 Ito s Formula 5 Recap Random Walk Consider
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 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 informationAnalysing multi-level Monte Carlo for options with non-globally Lipschitz payoff
Finance Stoch 2009 13: 403 413 DOI 10.1007/s00780-009-0092-1 Analysing multi-level Monte Carlo for options with non-globally Lipschitz payoff Michael B. Giles Desmond J. Higham Xuerong Mao Received: 1
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 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 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 informationModule 4: Monte Carlo path simulation
Module 4: Monte Carlo path simulation Prof. Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute Module 4: Monte Carlo p. 1 SDE Path Simulation In Module 2, looked at the case
More informationOrnstein-Uhlenbeck Theory
Beatrice Byukusenge Department of Technomathematics Lappeenranta University of technology January 31, 2012 Definition of a stochastic process Let (Ω,F,P) be a probability space. A stochastic process is
More informationMonte Carlo Methods for Uncertainty Quantification
Monte Carlo Methods for Uncertainty Quantification Abdul-Lateef Haji-Ali Based on slides by: Mike Giles Mathematical Institute, University of Oxford Contemporary Numerical Techniques Haji-Ali (Oxford)
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 informationMath 489/Math 889 Stochastic Processes and Advanced Mathematical Finance Dunbar, Fall 2007
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 Math 489/Math 889 Stochastic
More informationContinuous Processes. Brownian motion Stochastic calculus Ito calculus
Continuous Processes Brownian motion Stochastic calculus Ito calculus Continuous Processes The binomial models are the building block for our realistic models. Three small-scale principles in continuous
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 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 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 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 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 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 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 informationReading: You should read Hull chapter 12 and perhaps the very first part of chapter 13.
FIN-40008 FINANCIAL INSTRUMENTS SPRING 2008 Asset Price Dynamics Introduction These notes give assumptions of asset price returns that are derived from the efficient markets hypothesis. Although a hypothesis,
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 informationBROWNIAN MOTION Antonella Basso, Martina Nardon
BROWNIAN MOTION Antonella Basso, Martina Nardon basso@unive.it, mnardon@unive.it Department of Applied Mathematics University Ca Foscari Venice Brownian motion p. 1 Brownian motion Brownian motion plays
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 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 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 information"Vibrato" Monte Carlo evaluation of Greeks
"Vibrato" Monte Carlo evaluation of Greeks (Smoking Adjoints: part 3) Mike Giles mike.giles@maths.ox.ac.uk Oxford University Mathematical Institute Oxford-Man Institute of Quantitative Finance MCQMC 2008,
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 informationMonte Carlo Simulation of Stochastic Processes
Monte Carlo Simulation of Stochastic Processes Last update: January 10th, 2004. In this section is presented the steps to perform the simulation of the main stochastic processes used in real options applications,
More informationComputer labs. May 10, A list of matlab tutorials can be found under
Computer labs May 10, 2018 A list of matlab tutorials can be found under http://snovit.math.umu.se/personal/cohen_david/teachlinks.html Task 1: The following MATLAB code generates (pseudo) uniform random
More informationMonte Carlo Methods in Option Pricing. UiO-STK4510 Autumn 2015
Monte Carlo Methods in Option Pricing UiO-STK4510 Autumn 015 The Basics of Monte Carlo Method Goal: Estimate the expectation θ = E[g(X)], where g is a measurable function and X is a random variable such
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 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 informationAD in Monte Carlo for finance
AD in Monte Carlo for finance Mike Giles giles@comlab.ox.ac.uk Oxford University Computing Laboratory AD & Monte Carlo p. 1/30 Overview overview of computational finance stochastic o.d.e. s Monte Carlo
More informationComputer Exercise 2 Simulation
Lund University with Lund Institute of Technology Valuation of Derivative Assets Centre for Mathematical Sciences, Mathematical Statistics Fall 2017 Computer Exercise 2 Simulation This lab deals with pricing
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 informationPricing Dynamic Solvency Insurance and Investment Fund Protection
Pricing Dynamic Solvency Insurance and Investment Fund Protection Hans U. Gerber and Gérard Pafumi Switzerland Abstract In the first part of the paper the surplus of a company is modelled by a Wiener process.
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 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 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 informationProbability. An intro for calculus students P= Figure 1: A normal integral
Probability An intro for calculus students.8.6.4.2 P=.87 2 3 4 Figure : A normal integral Suppose we flip a coin 2 times; what is the probability that we get more than 2 heads? Suppose we roll a six-sided
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 information1 Mathematics in a Pill 1.1 PROBABILITY SPACE AND RANDOM VARIABLES. A probability triple P consists of the following components:
1 Mathematics in a Pill The purpose of this chapter is to give a brief outline of the probability theory underlying the mathematics inside the book, and to introduce necessary notation and conventions
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 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 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 informationThe value of foresight
Philip Ernst Department of Statistics, Rice University Support from NSF-DMS-1811936 (co-pi F. Viens) and ONR-N00014-18-1-2192 gratefully acknowledged. IMA Financial and Economic Applications June 11, 2018
More informationFinancial Engineering MRM 8610 Spring 2015 (CRN 12477) Instructor Information. Class Information. Catalog Description. Textbooks
Instructor Information Financial Engineering MRM 8610 Spring 2015 (CRN 12477) Instructor: Daniel Bauer Office: Room 1126, Robinson College of Business (35 Broad Street) Office Hours: By appointment (just
More informationarxiv: v1 [math.pr] 15 Dec 2011
Parameter Estimation of Fiber Lay down in Nonwoven Production An Occupation Time Approach Wolfgang Bock, Thomas Götz, Uditha Prabhath Liyanage arxiv:2.355v [math.pr] 5 Dec 2 Dept. of Mathematics, University
More informationComputer Exercise 2 Simulation
Lund University with Lund Institute of Technology Valuation of Derivative Assets Centre for Mathematical Sciences, Mathematical Statistics Spring 2010 Computer Exercise 2 Simulation This lab deals with
More informationCash Accumulation Strategy based on Optimal Replication of Random Claims with Ordinary Integrals
arxiv:1711.1756v1 [q-fin.mf] 6 Nov 217 Cash Accumulation Strategy based on Optimal Replication of Random Claims with Ordinary Integrals Renko Siebols This paper presents a numerical model to solve the
More informationUsing of stochastic Ito and Stratonovich integrals derived security pricing
Using of stochastic Ito and Stratonovich integrals derived security pricing Laura Pânzar and Elena Corina Cipu Abstract We seek for good numerical approximations of solutions for stochastic differential
More informationApplications of Stochastic Processes in Asset Price Modeling
Applications of Stochastic Processes in Asset Price Modeling TJHSST Computer Systems Lab Senior Research Project 2008-2009 Preetam D Souza May 26, 2009 Abstract Stock market forecasting and asset price
More informationStochastic Calculus, Application of Real Analysis in Finance
, Application of Real Analysis in Finance Workshop for Young Mathematicians in Korea Seungkyu Lee Pohang University of Science and Technology August 4th, 2010 Contents 1 BINOMIAL ASSET PRICING MODEL Contents
More informationNumerical Simulation of Stochastic Differential Equations: Lecture 2, Part 2
Numerical Simulation of Stochastic Differential Equations: Lecture 2, Part 2 Des Higham Department of Mathematics University of Strathclyde Montreal, Feb. 2006 p.1/17 Lecture 2, Part 2: Mean Exit Times
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 informationConvergence Analysis of Monte Carlo Calibration of Financial Market Models
Analysis of Monte Carlo Calibration of Financial Market Models Christoph Käbe Universität Trier Workshop on PDE Constrained Optimization of Certain and Uncertain Processes June 03, 2009 Monte Carlo Calibration
More informationMath489/889 Stochastic Processes and Advanced Mathematical Finance Solutions to Practice Problems
Math489/889 Stochastic Processes and Advanced Mathematical Finance Solutions to Practice Problems Steve Dunbar No Due Date: Practice Only. Find the mode (the value of the independent variable with the
More informationParameter estimation of diffusion models from discrete observations
221 Parameter estimation of diffusion models from discrete observations Miljenko Huzak Abstract. A short review of diffusion parameter estimations methods from discrete observations is presented. The applicability
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 informationMONTE CARLO EXTENSIONS
MONTE CARLO EXTENSIONS School of Mathematics 2013 OUTLINE 1 REVIEW OUTLINE 1 REVIEW 2 EXTENSION TO MONTE CARLO OUTLINE 1 REVIEW 2 EXTENSION TO MONTE CARLO 3 SUMMARY MONTE CARLO SO FAR... Simple to program
More informationAdvanced computational methods X SDE Lecture 5
Advanced computational methods X071521-SDE Lecture 5 1 Other weak schemes Here, I list out some typical weak schemes. If you are interested in them, you can read in details. 1.1 Weak Taylor approximations
More informationUsing R for teaching financial mathematics and statistics
MSOR Connections Vol No Spring Term 20 Julian Stander and John Eales Using R for teaching financial mathematics and statistics Julian Stander School of Computing and Mathematics University of Plymouth
More informationVolatility Smiles and Yield Frowns
Volatility Smiles and Yield Frowns Peter Carr NYU CBOE Conference on Derivatives and Volatility, Chicago, Nov. 10, 2017 Peter Carr (NYU) Volatility Smiles and Yield Frowns 11/10/2017 1 / 33 Interest Rates
More informationA THREE-FACTOR CONVERGENCE MODEL OF INTEREST RATES
Proceedings of ALGORITMY 01 pp. 95 104 A THREE-FACTOR CONVERGENCE MODEL OF INTEREST RATES BEÁTA STEHLÍKOVÁ AND ZUZANA ZÍKOVÁ Abstract. A convergence model of interest rates explains the evolution of the
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 informationBluff Your Way Through Black-Scholes
Bluff our Way Through Black-Scholes Saurav Sen December 000 Contents What is Black-Scholes?.............................. 1 The Classical Black-Scholes Model....................... 1 Some Useful Background
More informationMarket Risk: FROM VALUE AT RISK TO STRESS TESTING. Agenda. Agenda (Cont.) Traditional Measures of Market Risk
Market Risk: FROM VALUE AT RISK TO STRESS TESTING Agenda The Notional Amount Approach Price Sensitivity Measure for Derivatives Weakness of the Greek Measure Define Value at Risk 1 Day to VaR to 10 Day
More informationAppendix G: Numerical Solution to ODEs
Appendix G: Numerical Solution to ODEs The numerical solution to any transient problem begins with the derivation of the governing differential equation, which allows the calculation of the rate of change
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 informationNumerical Simulation of Stochastic Differential Equations: Lecture 1, Part 1. Overview of Lecture 1, Part 1: Background Mater.
Numerical Simulation of Stochastic Differential Equations: Lecture, Part Des Higham Department of Mathematics University of Strathclyde Course Aim: Give an accessible intro. to SDEs and their numerical
More informationRISK-NEUTRAL VALUATION AND STATE SPACE FRAMEWORK. JEL Codes: C51, C61, C63, and G13
RISK-NEUTRAL VALUATION AND STATE SPACE FRAMEWORK JEL Codes: C51, C61, C63, and G13 Dr. Ramaprasad Bhar School of Banking and Finance The University of New South Wales Sydney 2052, AUSTRALIA Fax. +61 2
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 informationParameter estimation in SDE:s
Lund University Faculty of Engineering Statistics in Finance Centre for Mathematical Sciences, Mathematical Statistics HT 2011 Parameter estimation in SDE:s This computer exercise concerns some estimation
More informationOrnstein-Uhlenbeck Processes. Michael Orlitzky
Ornstein-Uhlenbeck Processes Introduction Goal. To introduce a new financial dervative. No fun. I m bad at following directions. The derivatives based on Geometric Brownian Motion don t model reality anyway.
More informationStochastic Processes and Advanced Mathematical Finance. Limitations of the Black-Scholes Model
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 6 Analyzing Accumulated Change: Integrals in Action
Chapter 6 Analyzing Accumulated Change: Integrals in Action 6. Streams in Business and Biology You will find Excel very helpful when dealing with streams that are accumulated over finite intervals. Finding
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 information