In this section, we briefly go over some problems for stochastic calculus, the counterpart of stochastic processes in continuous space. Since the basic definitions and theorems of Brownian motion and stochastic calculus are directly used as interview problems, we’ll simply integrate them into the problems instead of starting with an overview of definitions and theorems.
Brownian motion
A. Define and enumerate some properties of a Brownian motion?1
Solution: This is the most basic Brownian motion question. Interestingly, part of the definition, such as W(0)=0, and some properties are so obvious that we often fail to recite all the details.
A continuous stochastic process W(t),t≥0, is a Brownian motion if
W(0)=0;
The increments of the process W(t1)−W(0),W(t2)−W(t1),…,W(tn)−W(tn−1), ∀0≤t1≤t2≤⋯≤tn are independent;
Each of these increments is normally distributed with distribution W(ti+1)−W(ti)∼N(0,ti+1−ti).
Some of the important properties of Brownian motion are the following: continuous (no jumps); E[W(t)]=0; E[W(t)2]=t; W(t)∼N(0,t); martingale property E[W(t+s)∣W(t)]=W(t); cov(W(s),W(t))=s,∀0<s<t; and Markov property (in continuous space).
There are two other important martingales related to Brownian motion that are valuable tools in many applications.
Y(t)=W(t)2−t is a martingale.
Z(t)=exp{λW(t)−21λ2t}, where λ is any constant and W(t) is a Brownian motion, is a martingale. (Exponential martingale).
130
We’ll show a proof of the first martingale using Ito’s lemma in the next section. A sketch for the exponential martingale is the following:2
B. What is the correlation of a Brownian motion and its square?
Solution: The solution to this problem is surprisingly simple. At time t, Bt∼N(0,t), by symmetry, E[Bt]=0 and E[Bt3]=0. Applying the equation for covariance Cov(X,Y)=E[XY]−E[X]E[Y], we have Cov(Bt,Bt2)=E[Bt3]−E[Bt]E[Bt2]=0−0=0. So the correlation of a Brownian motion and its square is 0, too.
C. Let Bt be a Brownian motion. What is the probability that B1>0 and B2<0?
Solution: A standard solution takes advantage of the fact that B1∼N(0,1), and B2−B1 is independent of B1, which is again a normal distribution: B2−B1∼N(0,1). If B1=x>0, then for B2<0, we must have B2−B1<−x.
But do we really need the integration step? If we fully take advantage of the facts that B1 and B2−B1 are two IID N(0,1), the answer is no. Using conditional probability and independence, we can reformulate the equation as
This approach is better demonstrated in Figure 5.10. When we have B1>0 and B2−B1<−B1, which accounts for 1/8 of the density volume. (All 8 regions separated by x=0, y=0, y=x, and y=−x have the same density volume by symmetry.)
Figure 5.10 Probability density graph of (B1,B2−B1)
Stopping time/ first passage time
A. What is the mean of the stopping time for a Brownian motion to reach either -1 or 1?
Solution: As we have discussed, Bt2−t is martingale. It can be proved by applying Ito’s lemma:
So d(Bt2−t) has no drift term and is a martingale. Let T=min{t;Bt=1 or −1}. At continuous time and space, the following property still applies: A martingale stopped at 132 a stopping time is a martingale! So BT2−T is a martingale and E[BT2−T]=B02−0=0. The probability that Bt hits 1 or -1 is 1, so BT2=1⇒E[T]=E[BT2]=1.
B. Let W(t) be a standard Wiener process and τx (x>0) be the first passage time to level x (τx=min{t;W(t)=x}). What is the probability density function of τx and the expected value of τx?
Solution: This is a textbook problem that is elegantly solved using the reflection principle, so we will simply summarize the explanation. For any Wiener process paths that reach x before t (τx≤t), they have equal probability ending above x or below x at time t, P(τx≤t,W(t)≥x)=P(τx≤t,W(t)≤x). The explanation lies in the reflection principle. As shown in Figure 5.11, for each path that reaches x before t and is at a level y above x at time t, we can switch the sign of any move starting from τx and the reflected path will end at 2x−y that is below x at time t. For a standard Wiener process (Brownian motion), both paths have equal probability.
From part A, it’s easy to show that the expected stopping time to reach either α (α>0) or −β (β>0) is again E[N]=αβ. The expected first passage time to level x is
133
essentially the expected stopping time to reach either x or −∞ and E[τx]=x×∞=∞. Although we have P(τx≤∞)=2−2N(x/∞)=1, the expected value of τx is ∞!
Figure 5.11 Sample path of a standard Weiner process and its reflected path
C. Suppose that X is a Brownian motion with no drift, i.e. dX(t)=dW(t). If X starts at 0, what is the probability that X hits 3 before hitting -5? What if X has drift m, i.e. dX(t)=mdt+dW(t)?
Solution: A Brownian motion is a martingale. Let p3 be the probability that the Brownian motion hits 3 before -5. Since a martingale stopped at a stopping time is a martingale, we have 3P3+(−5)(1−P3)=0⇒P3=5/8. Similar to random walk, if we have stopping boundaries (α>0) and −β (β>0), the probability that it stops at α instead of −β is pα=β/(α+β). The expected stopping time to reach either α or −β is again E[N]=αβ.
When X has drift m, the process is no longer a martingale. Let P(t,x) be the probability that the process hits 3 before hitting -5 when X=x at time t. Although X is no longer a 134 martingale process, it is still a Markov process. So P(t,x)=P(x) is actually independent of t. Applying the Feynman-Kac equation4, we have
mPx(x)+1/2Pxx(x)=0for −5<x<3.
We also have boundary conditions that P(3)=1 and P(−5)=0.
mPx(x)+1/2Pxx(x)=0 is a homogeneous linear differential equation with two real roots: r1=0 and r2=−2m. So the general solution is P(x)=c1e0x+c2e−2mx=c1+c2e−2mx. Applying the boundary conditions, we have
A different and simpler approach takes advantage of the exponential martingale: Z(t)=exp{λW(t)−21λ2t}. Since W(t)=X(t)−mt, X(t)−mt is a Brownian motion as well. Applying the exponential martingale, we have E[exp(λ(X−mt)−21λ2t)]=1 for any constant λ. To remove the terms including time t, we can set λ=−2m and the equation becomes E[exp(−2mX)]=1. Since a martingale stopped at a stopping time is a martingale, we have P3exp(−2m×3)+(1−P3)exp(−2m×−5)=1⇒e10m−e−6me10m−1.
D. Suppose that X is a generalized Weiner process dX=dt+dW(t), where W(t) is a Brownian motion. What is the probability that X ever reaches -1?
Solution: To solve this problem, we again can use the equation E[exp(−2mX)]=1 from the previous problem with m=1. It may not be obvious since we only have one apparent boundary, -1. To apply the stopping time, we also need a corresponding positive boundary. To address this problem, we can simply use +∞ as the positive boundary and the equation becomes
Ito’s lemma is the stochastic counterpart of the chain rule in ordinary calculus. Let X(t) be an Ito process satisfying dX(t)=β(t,X)dt+γ(t,X)dW(t), and f(X(t),t) be a twice-differentiable function of X(t) and t. Then f(X(t),t) is an Ito process satisfying
A. Let Bt be a Brownian motion and Zt=tBt. What is the mean and variance of Zt? Is Zt a martingale process?
Solution: As a Brownian motion, Bt∼N(0,t), which is symmetric about 0. Since t is a constant at t, Zt=tBt is symmetric about 0 and has mean 0 and variance t×var(Bt)=t2. More exactly, Zt∼N(0,t2).
Although Zt has unconditional expected value 0, it is not a martingale. Applying Ito’s lemma to Zt=tBt, we have dZt=∂Bt∂ZtdBt+∂t∂Ztdt+21×∂Bt2∂2Ztdt=21t−1/2Btdt+tdBt. For all the cases that Bt=0, which has probability 1, the drift term 21t−1/2Btdt is not zero.5 Hence, the process Zt=tBt is not a martingale process.
B. Let W(t) be a Brownian motion. Is W(t)3 a martingale process?
136
Solution: Applying Ito’s lemma to f(W(t),t)=W(t)3, we have ∂W(t)∂f=3W(t)2, ∂t∂f=0, ∂W(t)2∂2f=6W(t), and df(W(t),t)=3W(t)dt+3W(t)2dW(t). So again for the cases W(t)=0, which has probability 1, the drift term is not zero. Hence, W(t)3 is not a martingale process.