Financial Maths Notes.
Financial Maths Notes.
Preface
The aim of these notes is to present (some of) the fundamental features of
financial mathematics in a rigorous way but avoiding stochastic calculus (stochastic
differential equations.) This means that time is discrete, and the continuous case
is considered as the limiting case when the time intervals → 0.
These notes are not intended to be a stand alone text on financial mathemat-
ics; rather, they are intended to be a mathematical supplement to some text on
financial economics. A good such text book is John C. Hull’s: Options, Futures,
and Other Derivatives (Prentice Hall.)
Harald Lang
Lecture Note 1
Assume that we want to buy a quantity of coffee beans with delivery in nine
months. However, we are concerned about what the (spot) price of coffee beans
might be then, so we draw up a contract where we agree on the price today. There
are now at least three ways in which we can arrange the payment: 1) we pay now,
in advance. We call this price the present price of coffee beans with delivery in
nine months time, and denote it by P . Note that this is completely different from
the spot price of coffee beans, i.e., the price of coffee for immediate delivery. 2)
we pay when the coffee is delivered, i.e., in nine months time. This price is the
forward price, which we denote by G. 3) we might enter a futures contract with
delivery in nine months time. A futures contract works as follows:
Let us denote the days from today to delivery by the numbers 0, 1, . . . , T , so
that day 0 is today, and day T is the day of delivery. Each day n a futures price
Fn of coffee beans with delivery day T is noted, and FT equals the prevailing spot
price of coffee beans day T . The futures price Fj is not known until day j; it will
depend on how the coffee bean crop is doing, how the weather has been up to that
day and the weather prospects up till day T , the expected demand for coffee, and
so on. One can at any day enter a futures contract, and there is no charge for doing
so. The long holder of the contract will each day j receive the amount Fj − Fj−1
(which may be negative, in which he has to pay the corresponding amount,) so if
I enter a futures contract at day 0, I will day one receive F1 − F0 , day two F2 − F1
and so on, and day T , the day of delivery, FT − FT −1 . The total amount I receive
is thus FT − F0 . There is no actual delivery of coffee beans, but if I at day T buy
the beans at the spot price FT , I pay FT , get my coffee beans and cash the amount
FT − F0 from the futures contract. In total, I receive my beans, and pay F0 , and
since F0 is known already day zero, the futures contract works somewhat like a
forward contract. The difference is that the value FT − F0 is paid out successively
during the time up to delivery rather than at the time of delivery.
The simplest of these three contracts is the one when we pay in advance, at
least if the good that is delivered is non-pecuniary, since in that case the interest
does not play a part. For futures contrats, the interest rate clearly plays a part,
since the return of the contract is spread out over time.
We will derive some book-keeping relations between the present prices, for-
ward prices and futures prices, but first we need some interest rate securities.
1
Lecture Note 1
Theorem
The following relations hold:
a) P , G and F0 are linear functions, i.e., if X and Y are random payments
made at time T , then for any constants a and b P (T ) [aX + bY ] = aP (T ) [X] +
bP (T ) [Y ], and similarly for G and F0 .
(T )
b) G(T ) [1] = 1, F0 [1] = 1 and P (T ) [1] = ZT
c) P (T ) [X] = ZT G(T ) [X]
(T )
d) P (T ) [XeR(0,T ) ] = F0 [X]
(T )
e) P (T ) [X] = F0 [Xe−R(0,T ) ]
Proof
The proof relies on an assumption of the model: the law of one price. It
means that there can not be two contracts that both yield the same payoff X at
time T , but have different prices today. Indeed, if there were two such contracts,
we would buy the cheaper and sell the more expensive, and make a profit today,
and have no further cash flows in the future. But so would everyone else, and this
is inconsistent with market equilibrium. In Lecture Note 4 we will extend this
model assumption somewhat.
To prove c), note that if we take a long position on a forward contract on X
and at the same time a long position of G(T ) [X] zero coupon bonds, then we have
a portfolio which costs ZT G(T ) [X] today, and yields the income X at time T . By
the law of one price, it hence must be that c) holds.
To prove d), consider the following strategy: Deposit F0 on the money market
account, and take er1 long positions on the futures contract on X for delivery at
time T .
2
Lecture Note 1
The next day the total balance is then F0 er1 + er1 (F1 − F0 ) = F1 er1 . Deposit
this on the money market account, and increase the futures position to er1 +r2
contracts.
The next day, day 2, the total balance is then F1 er1 +r2 + er1 +r2 (F2 − F1 ) =
F2 er1 +r2 . Deposit this on the money market account.
And so on, up to day T when the total balance is FT er1 +···+rT = Xer1 +···+rT .
In this way we have a strategy which is equivalent to a contract where we pay F0
day zero, and receive the value FT er1 +···+rT = Xer1 +···+rT day T . This proves d).
Since relation d) is true for any random variable X whose outcome is known
day T , we may replace X by Xe−R(0,T ) in that relation. This proves e).
It is now easy to prove b). The fact that P (T ) [1] = ZT is simply the definition
of ZT . The relation G(T ) [1] = 1 now follows from c) with X = 1. In order to
(T )
prove that F0 [1] = 1, note that by the definition of money market account, the
price needed to be paid day zero in order to receive eR(0,T ) day T is 1; hence
(T )
1 = P (T ) [eR(0,T ) ]. The relation F0 [1] = 1 now follows from d) with X = 1.
Finally, to prove a), note that if we buy a contracts which cost P (T ) [X] day
zero and gives the payoff X day T , and b contracts that gives payoff Y , then we
have a portfolio that costs aP (T ) [X] + bP (T ) [Y ] day zero and yields the payoff
aX + bY day T ; hence P (T ) [aX + bY ] = aP (T ) [X] + bP (T ) [Y ]. The other two
relations now follow immediately employing c) and d). This completes the proof.
Hence
(T ) (T ) (T )
ZT G(T ) [X] = P (T ) [X] = F0 [Xe−R(0,T ) ] = F0 [X]e−R(0,T ) = F0 [X]ZT
and hence
(T )
G(T ) [X] = F0 [X].
3
Lecture Note 1
Corollary
If interest rates are deterministic, the forward price and the futures price
(T )
coincide: G(T ) [X] = F0 [X]
The equality of forward- and futures prices does not in general hold if interest
rates are random, though. To see this, we show as an example that if eR(0,T ) is
(T )
random, then F0 [eR(0,T ) ] > G(T ) [eR(0,T ) ].
Indeed, note that the function y = x1 is convex for x > 0. This implies that
1
its graph lies over its tangent. Let y = m + k(x − m) be the tangent line through
1 1 1
the point (m, m ). Then x ≥ m + k(x − m) with equality only for x = m (we
consider only positive values of x.) Now use this with x = e−R(0,T ) and m = ZT .
We then have
eR(0,T ) ≥ ZT−1 + k(e−R(0,T ) − ZT )
where the equality holds only for one particular value of R(0, T ). In the absence
of arbitrage (we will come back to this in Lecture Note 4), the futures price of the
value of the left hand side is greater than the futures price of the value of the right
hand side, i.e.,
(T ) (T )
F0 [eR(0,T ) ] > ZT−1 + k(F0 [e−R(0,T ) ] − ZT )
(T )
But F0 [e−R(0,T ) ] = P (T ) [1] = ZT , so the parenthesis following k is equal to zero,
hence
(T )
F0 [eR(0,T ) ] > ZT−1
(T )
so G(T ) [eR(0,T ) ] = ZT−1 , and it follows that F0 [eR(0,T ) ] > G(T ) [eR(0,T ) ].
In general, if X is positively correlated with the interest rate, then the futures
price tends to be higher than the forward price, and vice versa.
4
Lecture Note 1
could make an arbitrage by buying the share today, borrow for the cost and take
a short position on a forward contract. There is then no net payment today, and
none at T (deliver the share, collect the delivery price G of the forward and pay
the S0 erT = G for the loan.) But it would give the trader the dividend of the
share for free, for this dividend goes to the holder of the share, not the holder of
the forward. Likewise, the holder of the share might have the possibility of taking
part in the annual meeting of the company, so there might be a convenience yield.
5
Lecture Note 2
Example 1.
Consider a share of a stock which costs S0 today, and which gives a known
dividend amount d in t years, and whose (random) spot price at time T > t is ST .
Assume that there are no other dividends or other convenience yield during the
time up to T . What is the forward price G on this stock for delivery at time T ?
Assume that we buy the stock today, and sell it at time T . The present value
of the dividend is Zt d and the present value of the income ST at time T is ZT G.
Hence we have the relation
S 0 = Zt d + ZT G
from which we can solve for G
Example 2.
Consider a share of a stock which costs S0 now, and which gives a known
dividend yield d St i t years, where St is the spot price immediately before the
dividend is payed out. Let the (random) spot price at time T > t be ST . Assume
that there are no other dividends or other convenience yield during the time up
to T . What is the forward price G on this stock for delivery at time T ?
Consider the strategy of buying the stock now, and sell it at time t immedi-
ately before the dividend is paid out. With the notation of Lecture Note 1, we
have the relation
P (t) (St ) = S0 (1)
Consider now the strategy of buying the stock now, cash the dividend at time
t, and eventually sell the stock at time T . The present value of the dividend is
d P (t) (St ) and the present value of the income ST at time T is ZT G. Hence we
have the relation
S0 = d P (t) (St ) + ZT G
If we combine with (1) we get
(1 − d)S0 = ZT G
Example 3.
With the same setting as in example 2, assume that there are dividend yield
payments at several points in time t1 , . . . , tn , each time with the amount d Stj . As
in the above example, we can buy the stock today and sell it just before the first
dividend is paid out, so the relation
6
Lecture Note 2
holds. For any k = 2, 3, . . . , n we can buy the stock at time tk−1 immediately before
the payment of the dividend, collect the dividend, and sell the stock immediately
before the dividend is paid out at time tk . The price of this strategy today is zero,
so we have
where we have used (2) to obtain the last equality. Hence, by the theorem of
Lecture Note 1, we have the relation
ZT G = (1 − d)n S0 (3)
Example 4.
We now consider the setting in example 3, but with a continuous dividend
yield ρ, i.e., for any small interval in time (t, t + δt) we get the dividend ρSt δt.
If we divide the time interval (0, T ) into a large number n of intervals of length
δt = T /n, we see from (3) that
ZT G = (1 − ρ δt)n S0
Example 5.
Assume we want to buy a foreign currency in t years time at an exchange
rate, the forward rate, agreed upon today. Assume that the interest on the foreign
currency is ρ and the domestic rate is r per year. Let X0 be the exchange rate
now (one unit of foreign currency costs X0 in domestic currency,) and let Xt be
the (random) exchange rate as of time t. Let G be the forward exchange rate.
Consider now the strategy: buy one unit of foreign currency today, buy foreign
zero coupon bonds for the amount, so that we have eρt worth of bonds at time
t when we sell the bonds. Since the exchange rate at that time is Xt , we get
Xt eρt in domestic currency. Since the price we have paid today is X0 we have the
relation
¡ ¢
X0 = P (t) Xt eρt = e−rt G eρt = G e(ρ−r)t
7
Lecture Note 2
From this we can solve for f . The interest rate f is the forward rate from t to T .
The floating rate between tj and tk is the zero coupon rate that prevails between
these two points in time. The amount that A pays at time tk is thus
Lk · (1/Z(tk−1 , tk ) − 1), where Z(tk−1 , tk ) of course is the price of the zero coupon
bond at time tk−1 that matures at tk . Note that this price is unknown today but
known at time tk−1 . The total amount that B will receive, and A will pay is thus
random.
On the other hand, party B pays A a fixed amount c at each of the times
t1 , . . . , tn . The question is what c ought to be in order to make this deal “fair”.
Consider the following strategy: Do nothing today, but at time tk−1 buy
1/Z(tk−1 , tk ) zero coupon bonds that mature at time tk . This gives a cash flow of
−1 at time tk−1 and 1/Z(tk−1 , tk ) at time tk The price today of this cash flow is
zero, so
n
X
PAB = Lk (Ztk−1 − Ztk )
1
8
Lecture Note 2
so we can calculate the fair value of c by solving the equation we get by setting
PBA = PAB .
9
Lecture Note 3
Lemma
If we choose the coefficients βi such that Cov(Fti , e) = 0 for i = 1, . . . , n, then
the variance Var(e) is minimised.
Proof
Assume that we have chosen βi such that Cov(Fti , e) = 0 for i = 1, . . . , n, and
consider any other choise of coefficients:
n
X
S= γi Fti + f
i=1
10
Lecture Note 3
Q.E.D.
Theorem
The set of coefficients βi that minimises the variance Var(e) of the residual is
the solution to the system
n
X
Cov(Ftj , Fti )βi = Cov(Ftj , S) j = 1, . . . , n.
i=1
Proof
The condition that Cov(Ftj , e) = 0 means that
n
X
0= Cov(Ftj , S − βi Fti )
i=1
Xn
= Cov(Ftj , S) − βi Cov(Ftj , Fti )
i=1
If we regard S and Fti as any random variables, then the coefficients βi are called
the regression coefficients of S onto Ft1 , . . . , Ftn ; in the context here they are the
optimal hedge ratios when the futures Ft1 , . . . , Ftn are used to hedge S: for each
unit of volume of S we should use futures Fti corresponding to βi units of volume
in the hedge.
Let us consider the case n = 1. In this case we have that
Cov(S, Ft )
β=
Var(Ft )
2 Cov2 (S, Ft )
Var(S) = β Var(Ft ) + Var(e) = + Var(e)
Var(Ft )
= ρ2 (S, Ft ) Var(S) + Var(e)
p
where ρ(S, Ft ) = Cov(S, Ft )/ Var(S) Var(Ft ) is the correlation coefficient be-
tween S and Ft . ¡ ¢
Solving for Var(e) we get the plesant relation Var(e) = Var(S) 1 − ρ2 (S, Ft )
p
which means that the standard deviation of the hedged position is 1 − ρ2 (S, Ft )
times that of the unhedged position.
11
Lecture Note 4
Let us consider a two period market model of contracts where one party, the “long”
holder, pays a sum P ,—the present price—today when the contract is drawn up,
and the other part delivers some value X at some later date t. A contract is
defined by the random payoff X and the present price P ; we describe it by the
pair (P, X). The total cash flow for the long holder at time t is thus X, which
is random, whereas the cash flow today is −P and is deterministic. There is no
other cash flow in relation to entering the contract.
The market consists of a set of contracts, and we assume that we can com-
pose arbitrary portfolios of contracts, i.e., if we have contracts 1, . . . , n with payoffs
X1 , . . . , Xn and present prices P1 , . . . , Pn , then we can compose a portfolio con-
sisting of λi units of contract i, i = 1, . . . , n where λ1 , . . . , λn are arbitrary real
numbers. I.e., we assume that we can take a short position in any contract, and
ignore anyPdivisibility problems. The total Pnpresent price of such a portfolio is
n
of course 1 λi Pi and the total payoff is 1 λi Xi . We call also such portfolios
“contracts”, so the set of contracts, defined as pairs (P, X), constitute a linear
space.
We assume that the market is arbitrage free in the sense that the law of one
price prevails (Lecture Note 1,) but also in the sense that
1. If X ≥ 0, then P (X) ≥ 0.
2. If X = 0 almost surely (abbr. “a.s.”, i.e., with probability 1) then P (X) = 0.
3. If X ≥ 0 and P (X) = 0, then X = 0 a.s.
Theorem
If the sample space is finite and the market is arbitrage free in the sense given
above, then there exists a random variable U such that U > 0, and for any
payoff X it holds that P (X) = E [XU ].
Proof
Let Ω = {ω1 , . . . , ωn , . . . , ωN } be the sample space and without loss of gen-
erality we may assume that ωk have positive probabilities pk for k ≤ n and zero
probability for k > n. ¡ ¢
We associate with each contract (P, X) the vector − P, X(ω1 ), . . . , X(ωn )
in Rn+1 i.e., we ignore X:s values on events with probability zero. The set of
such vectors constitute a linear subspace V of Rn+1 . We now prove that there is
a vector q̄ which is orthogonal to V all of whose coordinates are positive.
Let K be the subset of Rn+1 : K = {ū = (u0 , . . . , un ) ∈ Rn+1 | u0 +· · ·+un =
1 and ui ≥ 0 for all i}. Obviously K and V have no vector in common; indeed, it
12
Lecture Note 4
is easy to see that any such common vector would represent an arbitrage. Now let
q̄ be the vector of shortest Euclidean length such that q̄ = ū − v̄ for some vectors
ū and v̄ in K and V respectively. The fact that such a vector exists needs to be
proved, however, we will skip that proof. We write q̄ = ū∗ − v̄ ∗ where ū∗ ∈ K and
v̄ ∗ ∈ V.
Now note that for any t ∈ [0, 1] and any ū ∈ K, v̄ ∈ V , the vector tū + (1 −
t)ū ∈ K and tv̄ + (1 − t)v̄ ∗ ∈ V, hence |(tū + (1 − t)ū∗ ) − (tv̄ + (1 − t)v̄ ∗ )| as a
∗
for all contracts, and by scaling the qi :s we can arrange that q0 is equal to one,
hence
Xn
P = qi X(ωi ) (1)
1
We now define U (ωi ) = pqii for i = 1, . . . , n and U (ωi ) = arbitrarily positive for
i = n + 1, . . . , N . U is thus a positive random variable, and it is easy to verify
that it satisfies the conditions of the theorem. Q.E.D.
Note that by definition, the zero coupon price Zt = P (1) = E [U ]. Define the
random variable Q by Q = Zt−1 U , then Q > 0, E [Q] = 1 and
P (X) = Zt E [XQ]
G(X) = E [XQ].
13
Lecture Note 5
Black’s Model
Let us compute the price of a European derivative on some underlying asset
with value X at maturity t. We assume that there already is a forward contract
on the market, or that we can compute the forward price, and denote the current
forward price for delivery at t by G. Black’s model assumes that the value X has
a log-normal probability distribution:
√
X = Aeσ tz
where z ∈ N (0, 1) (1)
1 2
Hence A = G(X)eσλt− 2 σ t
which we substitute into (1):
1 2
√
X = G(X)eσλt− 2 σ t+σ t z
(2)
14
Lecture Note 5
Finally, if G denotes the forward price of the underlying asset, the present
price of the derivative is:
¡ 1 2
√ ¢¤
p = Zt E [f Ge− 2 σ t+σ t w where w ∈ N (0, 1) (3)
15
Lecture Note 5
E [X] = Geλσt
which means that the expected rate of return on a forward contract on X is λσ.
The expected rate of return is thus proportional to the standard deviation σ, and
λ is called the market price of risk.
The situation when λ > 0 is called backwardation, and when λ < 0 the
situation is known as contango.
16
Lecture Note 6
where Pt (y) is the value at time t, if the zero coupon interest rate at that time
equals y, which may differ from y0 .
The idea of the concepts yield (y) and duration (D) is to define a number D
for a coupon paying bond, or a portfolio of bonds, such that the relations (2) and
(3) still hold approximately. Hence, let P0 be the present value of a bond, or a
portfolio of bonds:
Xn
P0 = Zti ci
i=1
where Zs is the price today of a zero coupon bond maturing at time s and ci is
the payment received at time ti . Now define the yield y0 by the relation
n
X
P0 = e−y0 ti ci = P0 (y0 ) (4)
i=1
Here (20 ) follows from the definition of P̂D and (30 ) holds if the yield is y at time
t, as is seen as follows:
n
X n
X
−y(ti −t) yt
Pt (y) = e ci = e e−yti ci = eyt P0 (y) = P̂D e−y(D−t)
i=1 i=1
17
Lecture Note 6
Summary
A bond, or a portfolio of bonds, has a yield y0 , defined by (4), and a duration
D, defined by (5). At any time t before any coupon or other payments have been
paid out, the value Pt (y) of the asset is approximately equal to
and y is the yield prevailing at that time. In particular, the value of the asset at
time t = D is P̂D = P0 ey0 D , and hence independent of any changes in the yield to
a first order approximation.
Example
Consider a portfolio of bonds that gives the payment 1’000 after one year,
1’000 after two years and 2’000 after three years. Assume that Z1 = 0.945, Z2 =
0.890 and Z3 = 0.830. The present value of the portfolio is then
so the duration is D = 7’716.5 years = 2.208 years. The value of the portfolio
3’495
at time t, if the yield at that time is y and t is less than one year, can thus be
approximated by
Pt (y) = 3’995 e−y(2.208−t)
18
Lecture Note 6
When we arrive at time t it may be that the prevailing yield y at that time
is not the same as yF . Let y be the prevailing yield at time t. Then the value of
the security at that time is
n
X
Pt = Pt (y) = ci e−y(ti −t)
i=1
Pt = PF e−DF y
Summary
An interest rate security that after time t gives the deterministic payments
c1 , c2 , . . . , cn at times t1 < t2 < · · · < tn has a forward yield yF , defined by (6),
and a forward duration DF defined by (7). The value of the asset as of time t,
0 < t < t1 is then, to a first order approximation
Example
In the previous example, assume that Z1.5 = 0.914. The forward price G1.5
of the portfolio for delivery in one and a half years is obtained from the relation
(note that the first payment of 1’000 has already been paid out, and hence does
not contribute to the forward value:)
19
Lecture Note 6
Now
−G01,5 (yF ) = 0.5 e−0.5yF 1’000 + 1.5 e−1.5yF 2’000 = 3’215.8
The minus sign is there for notational convenience later on; note that z and −z
have the same distribution. Hence, the present value Pt at time t of the asset is
√ √
Pt = PF e−DF (α−σ t z)
= A eDF σ tz
where A = PF e−DF α .
Note that this is exactly the same specification as (1) in Lecture Note 5 with
σ replaced by DF σ. We √ can thus use Black’s formula for European options of
Lecture Note 5, where σ t in that formula stands for the product of the standard
deviation of the forward yield and the forward duration of the underlying asset.
Portfolio Immunising
Assume that we have portfolio P of bonds whose duration is DP . The duration
reflects its sensitivity to changes in the yield, so we might well want to change
that without making any further investments. One way to do this is to add a bond
futures to the portfolio (or possibly short such a futures.)
Consider again the security discussed under “Forward Yield and Forward Du-
ration” above. We will look at a futures contract on the value Pt . Using the no-
tation of Lecture Note 1, the futures price is F0 = F0 (Pt ) = F0 (PF e−DF y ), where
the yield y at time t is random as seen from today. Assume now that the present
yield y0 of the portfolio P changes to y0 + δy. The question arises what impact this
has on y. If y = (current yield) + (random variable independent of current yield),
then y should simply be replaced by y + δy. The new futures price should then be
20
Lecture Note 6
the marking-to-market is added to the portfolio, so after the shift the value of the
total portfolio is
1 0 F0
Dtot = − Ptot (y0 ) = DP + DF
P P
21
Lecture Note 7
In Lecture Note 4 we introduced the notion that a forward price can be expressed as
an expectation: G(t) (X) = E [XQ] for some random variable Q > with E [Q] = 1.
In probability theory, what we have done is called a change of probability measure.
If Q is as above, then it is convenient to introduce a notation for the expected
value E [XQ], for example
E∗ [X] = E [XQ]
and say that E∗ [X] is the expected value of X with respect to a new probability
measure Pr∗ , where the probability of an event A with this new probability measure
is defined by
Pr∗ (A) = E∗ [IA ] = E [IA Q]
(t)
G0 (X) = E(t) [X]
Choise of Numeraire
The forward probability measure introduced above is a very important exam-
ple of a risk adjusted probability measure, but one can also use other numeraires.
Let N be some random value whose outcome becomes known at some date t, and
assume that N > 0. For instance, N may be the value of one barrel of crude oil
at time t. Now define the random variable QN by
1
QN = (t)
UN
P0 (N )
22
Lecture Note 7
The verification is left as an exercise for the reader. If we denote by EN [·] the
expectation associated with QN , then
(t) (t) £X ¤
P0 (X) = P0 (N ) EN (1)
N
Here N is the numeraire. For example, if X is the price of jet fuel at time t and
the numeraire N is the price of crude oil at time t, the idea could be that we want
to calculate the present price for jet fuel to be delivered at time t, but rather than
creating a model for the price of jet fuel, we might find it easier to come up with a
reasonable model for how the market perceives the £ ratio
¤ X/N of jet fuel to crude
N X
oil. We can then from that model compute E N , whereas the present price
(t)
P0 (N ) can be directly observed on the market (since there is a futures market
for crude oil.)
If we choose N = one unit of some currency , i.e., the numeraire asset is a
zero coupon bond, then we have
(t) (t)
P0 (X) = P0 (1) EN [X] = Zt EN [X]
(t)
i.e., EN = G0 = E(t) , the forward mesaure with maturity date t.
Theorem 1
Assume that t < T and that the zero coupon bond price Z(t, T ) is known at
the present time 0. Then, for any random variable X whose outcome is known
at time t,
E(t) [X] = E(T ) [X]
Proof:
Consider the following strategy: Enter a contract forward contract on X
maturing at t (so G(t) = E(t) [X],) invest the net payment X − G(t) in zero coupon
bonds maturing at T . In this way we have constructed a contract which pays
Z(t, T )−1 (X − G(t) ) at time T and costs 0 today. Hence
23
Lecture Note 7
holds. We then say that Y is the conditional expectation at time t of the variable
X. It is easy to see that the variable Y is essentially uniquely determined by X,
i.e., if Y1 and Y2 both satisfy E [IA X] = E [IA Yi ] for all A, then Y1 = Y2 a.s.
The conditional expectation is denoted Et [X], i.e., Y = Et [X]. Note in par-
ticular the law of iterated expectations
£ ¤
E Et [X] = E [X] (3)
Martingale Prices
Consider a contract that gives the random payoff X at time T . The forward
(T )
price Gt , 0 ≤ t < T , of this contract at time t is a random variable whose
outcome is determined at time t.
24
Lecture Note 7
Theorem 2
(T )
For any t, 0 ≤ t ≤ T , the forward price Gt equals the conditional expectation
(T ) (T )
Gt = Et [X]
Hence, the forward process has the martingale property (4) w.r.t. the forward
measure, i.e.,
(T ) (T )
Gj = Ej [Gj+1 ]
Proof
Consider the following strategy: Let A be any event whose outcome occurs at
time t. Wait until time t and then, if A has occurred, enter a forward contract on
(T )
X maturing at T with forward price Gt , but if A has not occurred, do nothing.
(T )
We have thus constructed a contract which gives payoff IA (X − Gt ) at time T
which costs 0 today. Hence
(T ) (T )
ZT E(T ) [IA (X − Gt )] = 0 i.e., E(T ) [IA X] = E(T ) [IA Gt ]
Since this is true for any event A whose outcome is known at time t, we have by
the definition of conditional expectation
(T ) (T )
Gt = Et [X]
Q.E.D.
Remark
From this we deduce (c.f. Lecture Note 1) that the present price pt at time t
of the contract yielding X at time T is
(T )
pt = Z(t, T )Et [X]
where Z(t, T ) is the price at time t of a zero coupon bond maturing at time T .
Especially, if the interest rate is deterministic and equal to r, then the present
prices satisfy
pj = e−r ∆t E∗j [pj+1 ] (5)
25
Lecture Note 8
In Black’s model, there are only two relevant points in time: the time when the
contract is written, and the time when it matures. In many cases we have to
consider also what happens in between these two points in time.
First some notation. Time is discrete, t = t0 , . . . , tn = T , and we let for
notational simplicity t0 = 0 and all time spells tk − tk−1 = ∆t be equally long;
hence tk = k∆t. Let X be some random value whose outcome becomes known at
time T ; we are interested in evaluating the prices of derivatives of X.
Black-Scholes Dynamics
The Black-Scholes assumption about the random behaviour of X is that
n
X = Aeσ Σ1 zj (1)
where A is some constant, σ a parameter called the volatility, and the zj :s are
stochastically independent normal
√ N(0, ∆t) random variables (∆t is the variance;
the standard deviation is thus ∆t) where the outcome of zj occurs at time tj .
This is the model for the price dynamics under the true probability measure.
But in order to compute prices, we want the dynamics under the forward measure.
for some constants µ and λ. In order for E [Q] = 1 to be satisfied, we must choose
1 1
µ = − ∆t n λ2 = − λ2 T.
2 2
for a normally distributed random variable y ∈ N(0, v), as is easily checked. Our
suggestion thus looks like
2
1
T −λ Σn
Q = e− 2 λ 1 zi (4)
Now we prove
26
Lecture Note 8
1 n 2 n
E [h(z1 , . . . , zn )e− 2 Σ1 λj ∆t−Σ1 λj zj ] = E [h(w1 − λ1 ∆t, . . . , wn − λn ∆t)]
Proof
This is just a substitution of variables in a multiple integral:
1 2
E [h(z1 , . . . , zn )e− 2 Σλj ∆t−Σλj zj ]
Z Z
1 − 12 Σλ2j ∆t−Σλj xj
= . . . h(x 1 , . . . , x n )e ×
(2π∆t)n/2
1 n 2
× e− 2 ∆t Σ1 xj dx1 . . . dxn
Z Z
1 1 2
= n/2
. . . h(x1 , . . . , xn )e− 2 ∆t Σ(xj +λj ∆t) dx1 . . . dxn
(2π∆t)
= [change of variables: xj + λj ∆t = yj ]
Z Z
1 1 2
= n/2
. . . h(y1 − λ1 ∆t, . . . , yn − λn ∆t)e− 2 ∆t Σyj dy1 . . . dyn
(2π∆t)
= E [h(w1 − λ1 ∆t, . . . , wn − λn ∆t)]
Q.E.D.
Remark
Girsanov’s theorem generally refers to the more general case when ∆t → 0 so
the sum is replaced by an integral (“Itô integral”.)
We can now employ Girsanov’s theorem in our case. In our case λj = λ, and
the Girsanov transformation (4) hence defines a new measure, which will be our
forward measure E(T ) , such that under this measure (1) can be written
n
X = A eσ Σ1 (wj −λ∆t)
n (5)
= B eσ Σ1 wj
for some constant B and where wj are independent N(0, ∆t)-variables under the
(T ) (T ) 1 2
forward measure. Since G0 = E(T ) [X], it follows that B = G0 e− 2 σ T .
27
Lecture Note 8
Binomial Approximation
In almost all cases when the whole price path has to be taken into account, it is
necessary to use some numerical procedure to calculate the price of the derivative.
A useful numerical procedure is to approximate the specification (5) by a√binomial
tree. In (5), we replace wj by a binary variable bj which takes the value −√ ∆t with
probability (under the forward probability measure) 0.5 and the value + ∆t with
probability 0.5. For large values of n, and small ∆t, this is a good approximation—
in fact, as n → ∞ and ∆t → 0 the price of a derivative calculated from the binomial
tree (see next Lecture Note) will converge to the theoretical value it would have
from the specification (5).
The binomial specification is thus
n
X = B eσ Σ1 b j
√
Note that E(T ) [eσbj ] = cosh(σ ∆t). We know from the previous Lecture Note
(T )
that Gk (X) = Ek [X], so
k √
Gk = B eσ Σ1 bj coshn−k (σ ∆t)
Combining this relation with the same for Gk+1 , we get the relation
eσ bk+1 ¡ √
Gk+1 = Gk √ = Gk 1 ± ε) where ε = tanh(σ ∆t)
cosh(σ ∆t)
and the plus and minus sign occur with probability (under the forward probability
measure) 0.5 each.
Remarks
It is worthwhile to note that (5) is identical to (1) except for the constant
A; the volatility σ is the same. The parameter λ that appears in Q has a similar
interpretation as in Lecture Note 5: E [X] = G0 eλσT and is called the market
price of risk.
28
Lecture Note 9
G0 G0 u G0 u2 G0 u3 G0 u4 ···
G0 d G0 du G0 du G0 du3
2
···
G 0 d2 G0 d2 u G0 d2 u2 ···
G0 d3 G0 d3 u ···
G0 d4 ···
In this diagram, each column represents an instant of time, and moving from one
instant to the next (one step right) means that the forward price either moves “up”
or “down”. In this diagram “up” means going one step to the right on the same
line, and “down” going one step to the right on the line below (this configuration
is convenient in a spread sheet.) Under the forward measure the probabilities of
moving up and down are both 21 . The specification of u and d are as follows:
½
u=1+ε √
where ε = tanh(σ ∆t)
d=1−ε
where ∆t is the time spell between columns, and σ is the volatility of the forward
price.
Here the prices in the rightmost column are known, since they represent the value
of the option at maturity when the underlying futures price is that in the corre-
sponding position of the tree of futures prices. Then option prices are calculated
backwards:
¡ ¢
p̂j,k = e−r∆t 0.5 pj+1,k + 0.5 pj+1,k+1 (1)
pj,k = max[p̂j,k , value if exercised] (2)
29
Lecture Note 9
Indeed, p̂j,k = e−r∆t E∗ [pj+1,• ], which is just another way of writing (1), is the
value if the option is not exercised (see (5) of Lecture Note 7,) and the highest
value of this and the value if exercised is the present value. In order to calculate
the value if exercised, we employ the original tree of futures prices, of course. The
end result p is the price of the option today.
The price of the option is then calculated in the same way as above.
It is now easy to build the tree of stock prices: First compute G0 from (4), then
construct the tree for Gj , and from these prices, compute Sj from (3).
30
Lecture Note 9
For times j before i we can compute the stock price in accordance with example
2 of Lecture Note 2:
Sj (1 − d) = e−r(n−j)∆t Gj
In particular,
G0 = (1 − d)S0 ern∆t (7)
It is now easy to build the tree of stock prices: First compute G0 from (7), then
construct the tree for Gj , and from these prices, compute Sj from (5) and (6).
31
Lecture Note 10
We have assumed interest rates to be for the most part deterministic. In order to
study interest rate derivatives, or other situations where interest rates are assumed
to be random, we need a different probability measure than the forward measure.
One reason for this is that when interest rates are random, forward measures for
different maturity dates are not the same. Rather than having zero coupon bonds
as numeraires (see Lecture Note 7) we use the money market account (see Lecture
Note 1) as a numeraire. Hence, let the numeraire N be N = eR(0,t) where we have
used the notation introduced in Lecture Note 1 and 7. Then (see (1) of Lecture
Note 7)
(t) (t) b b
P0 (X) = P0 (eR(0,t) ) E[Xe −R(0,t)
] = E[Xe −R(0,t)
]
b (t) (t)
E[X] = P0 (XeR(0,t) ) = F0 (X) (1)
so it seems natural to call the probability measure associated with this numeraire
the futures measure. In the litterature it is often called the Equivalent Martingale
Measure (Hull uses the term “traditional risk neutral measure.”)
It is important to note that the forward measures E(T ) differ for different
maturities T when interest rates are random. This is not the case for the futures
measure, though. Indeed, we will now prove two theorems on the futures measure:
Theorem 1
The futures measure as defined above is independent of the date of maturity T
b (n) [X] =
in the following sense: if the outcome of X is known at time tn then E
b (m) [X] if m > n, where we denote by E
E b (k) the futures measure for contracts
maturing at time tk .
Theorem 2
The futures prices {Fj } have the martingale property w.r.t. the futures mea-
sure: Fj = Eb t [Fk ] for j < k. In particular, the futures price F0 is the expected
j
value w.r.t. the futures measure of X, the spot price at delivery.
32
Lecture Note 10
Proof of Theorem 1
Let P be the present price of the value XeR(0,n) to be delivered at time tn .
We can convert this contract to one where instead the value XeR(0,m) is delivered
at time tm by simply depositing the payoff XeR(0,n) in the money market account
up to time tm . The present price is of course the same. Using the theorem in
Lecture Note 1, we have the following equalities from the two contracts:
(tn ) (tm )
P = P (tn ) [XeR(0,n) ] = F0 [X] and P = P (tm ) [XeR(0,m) ] = F0 [X]
Q.E.D.
Proof of theorem 2
Consider the following strategy: Let A be any event whose outcome is known
at time tj < tk . At time tk−1 enter a long position of eR(0,k) futures contracts if A
has occurred, otherwise do nothing. At tk collect IA (Fk −Fk−1 )eR(0,k) (which may
be negative) and close the contract. This gives the payment IA (Fk − Fk−1 )eR(0,k)
at time tk and no other cash flow. The present price is zero, hence
b A (Fk − Fk−1 )]
0 = P (tk ) [IA (Fk − Fk−1 )eR(0,k) ] = F0 [IA (Fk − Fk−1 )] = E[I
Hence,
b A Fk ] = E[I
E[I b A Fk−1 ]
b A Fj ] = E[I
E[I b A Fk ]
Since A can be any event whose outcome is known at time tj , this means that (see
“Conditional Expectations and Martingales” in Lecture Note 7)
Theorem 3
The following formula for the present prices pj as of time tj for the asset X
to be delivered at time T holds:
b t [pk e−R(j,k) ],
pj = E j < k.
j
In particular,
b t [Xe−R(j,T ) ]
pj = E j
33
Lecture Note 10
Proof
With the notation of Lecture Note 1,
(T ) b t [Xe−R(j,T ) ]
p j = Fj [Xe−R(j,T ) ] = E j
and likewise
b t [Xe−R(k,T ) ]
pk = E k
Q.E.D.
34
Lecture Note 11
We divide time into short time intervals t0 , . . . , tn = T , tk −tk−1 = ∆t. Ideally, the
points in time tk should coincide with the times of settlement of futures contracts.
At tk−1 there is an interest rate rk ∆t prevailing from tk−1 to tk . This interest
rate is random, but its value is known at time tk−1 . The Ho-Lee model of the
interest rate is:
√
rk = ak + σ ∆t (z1 + · · · + zk−1 )
where ak are some numbers, the factor σ is the volatility of the short interest rate.
The zj :s are independent N(0,1)-variables and the outcome of each zj occurs at
time tj . This is under the true probability measure.
Under the futures measure, the model looks the same, with σ maintained, but
with new terms ak ; the reasoning is the same as in Lecture Note 8 (“remarks”.)
Thus, under the futures measure
√
rk = θk + σ ∆t (z1 + · · · + zk−1 )
We now normalise ∆t = 1. This means that rk is the one period interest rate—it
is proportional to ∆t whereas σ is proportional to ∆t3/2 . Now the model reads
Let Ztk be the price of a zero coupon bond maturing at tk with face value 1. Then,
according to Theorem 3, Lecture Note 10,
b e−r1 −···−rk ]
Ztk = E[1
In order for the model (1) to correctly represent the current term structure, we
must thus have
i.e.,
σ2
((k−1)2 +···+12 )
Ztk = e−θ1 −···−θk e 2
35
Lecture Note 11
b t [e−rt+1 −···−rT ]
Z(t, T ) = E
ZT b − σ2 (t2 +···+(T −1)2 )−σ(T −t)(z1 +···+zt )
= Et [e 2
Zt
× e−σ((T −t−1)zt+1 +(T −t−2)zt+2 +···+1 zT −1 ) ]
ZT − σ2 (t2 +···+(T −1)2 ) −σ(T −t)(z1 +···+zt )
= e 2 e
Zt
×E b t [e−σ((T −t−1)zt+1 +(T −t−2)zt+2 +···+1 zT −1 ) ]
ZT −σ(T −t)(z1 +···+zt ) − σ2 (t2 +···+(T −1)2 ) σ2 ((T −t−1)2 +···+12 )
= e e 2 e2
Zt
σ2 2 σ2
− (t + · · · + (T − 1)2 ) + ((T − t − 1)2 + · · · + 12 )
2 2
T −t−1
σ2 X
=− (t + j)2 − j 2
2 j=0
T −t−1
σ2 X
=− t (2j + t)
2 j=0
σ2
= − (T − 1)(T − t)t
2
Hence,
ZT − σ2 (T −1)(T −t)t −σ(T −t)(z1 +···+zt )
Z(t, T ) = e 2 e (3)
Zt
36
Lecture Note 11
regardless of the interest rate model. We leave the proof to the reader. The futures
price F0 is
b ZT e− σ2
2
b
F0 = E[Z(t, T )] = E[ (T −1)(T −t)t −σ(T −t)(z1 +···+zt )
e ]
Zt
σ2
= Ge− 2 (T −1)(T −t)t b −σ(T −t)(z1 +···+zt ) ]
E[e
σ2 σ2 σ2
(T −t)2 t
= Ge− 2 (T −1)(T −t)t
e 2 = Ge− 2 (T −t)(t−1)t
σ2 σ2
(T −t)t2
F0 = Ge− 2 (T −t)(t−∆t)t
−→ Ge− 2 when ∆t → 0.
We see that in contrast to the situation when interest rates are deterministic, the
forward price G and the futures price F0 differ. Since the price of the underlying
asset (the bond) is negatively correlated with the interest rate, the futures price
σ2
(T −t)t2
is lower; in the Ho-Lee model by factor e− 2 .
37
Lecture Note 11
Theorem
The variable zj is N (−σ(t − j), 1) under the forward measure with maturity
at time t ≥ j. Thus we may write zj = wj − σ(t − j) where wj ∈ N (0, 1) under
the forward (t) measure.
But ZZTt = G0 , the forward price of the underlying bond, hence, (also with no
normalisation of time,)
£ ¡ ¢¤
p = Zt E(t) F Z(t, T )
Z ∞
Zt ¡ σ2 2
√ ¢ 1 2
=√ F G0 e− 2 (T −t) t eσ(T −t) t x e− 2 x dx
2π −∞
Note that the same formula may be obtained by Black’s model for bond options
(Lecture Note 6.)
38
Lecture Note 12
k
X
rk = θk + σ bj k = 1, . . . , n
2
where {bk } are binomial random variables which takes the values ±1, each with
probability 0.5 under the futures measure; they are thus assumed to be statistically
independent. The outcome of the variable bk occurs at time k − 1. Since we don’t
need the dynamics under the true probability measure, we model the dynamics
from the start under the futures measure. However, as in Lecture Note 11, σ
represents the volatility of the interest rate under either measure, so it can be
3
measured from real data; σ is proportional to ∆t 2 just as in Lecture Note 11.
We will choose the terms θj :s such that the model is consistent with the
current term structure—it will be close to, but not exactly, the same as in the
Normal distribution case.
Let us compute the price Z4 at time t0 of a zero coupon bond maturing at
time t4 with face value 1:
³Z ´ £ ¡ ¢¤ k
X
k−1
rk = ln + ln cosh (k − 1)σ + σ bj k = 1, . . . , n
Zk 2
(
1
bj = 1 with probability 2
where 1
bj = −1 with probability 2
Once we have the parameters of the model, we can price interest rate deriva-
tives in the binomial tree. We show the procedure by an example, where we want
39
Lecture Note 12
t2 t3
84.794 91.119
88.254 92.960
91.856 94.838
96.754
The value of the option can now also be obtained by backward recursion:
t0 t1 t2
2.309 1.050 0
3.853 2.254
5.856
The value of the option is thus 2.309. It is also easy to price other more exotic
derivatives in this binomial tree model.
40