last updated: September 19, 2024

Stochastic Processes and Financial Mathematics
(part two)

\(\newcommand{\footnotename}{footnote}\) \(\def \LWRfootnote {1}\) \(\newcommand {\footnote }[2][\LWRfootnote ]{{}^{\mathrm {#1}}}\) \(\newcommand {\footnotemark }[1][\LWRfootnote ]{{}^{\mathrm {#1}}}\) \(\let \LWRorighspace \hspace \) \(\renewcommand {\hspace }{\ifstar \LWRorighspace \LWRorighspace }\) \(\newcommand {\mathnormal }[1]{{#1}}\) \(\newcommand \ensuremath [1]{#1}\) \(\newcommand {\LWRframebox }[2][]{\fbox {#2}} \newcommand {\framebox }[1][]{\LWRframebox } \) \(\newcommand {\setlength }[2]{}\) \(\newcommand {\addtolength }[2]{}\) \(\newcommand {\setcounter }[2]{}\) \(\newcommand {\addtocounter }[2]{}\) \(\newcommand {\arabic }[1]{}\) \(\newcommand {\number }[1]{}\) \(\newcommand {\noalign }[1]{\text {#1}\notag \\}\) \(\newcommand {\cline }[1]{}\) \(\newcommand {\directlua }[1]{\text {(directlua)}}\) \(\newcommand {\luatexdirectlua }[1]{\text {(directlua)}}\) \(\newcommand {\protect }{}\) \(\def \LWRabsorbnumber #1 {}\) \(\def \LWRabsorbquotenumber "#1 {}\) \(\newcommand {\LWRabsorboption }[1][]{}\) \(\newcommand {\LWRabsorbtwooptions }[1][]{\LWRabsorboption }\) \(\def \mathchar {\ifnextchar "\LWRabsorbquotenumber \LWRabsorbnumber }\) \(\def \mathcode #1={\mathchar }\) \(\let \delcode \mathcode \) \(\let \delimiter \mathchar \) \(\def \oe {\unicode {x0153}}\) \(\def \OE {\unicode {x0152}}\) \(\def \ae {\unicode {x00E6}}\) \(\def \AE {\unicode {x00C6}}\) \(\def \aa {\unicode {x00E5}}\) \(\def \AA {\unicode {x00C5}}\) \(\def \o {\unicode {x00F8}}\) \(\def \O {\unicode {x00D8}}\) \(\def \l {\unicode {x0142}}\) \(\def \L {\unicode {x0141}}\) \(\def \ss {\unicode {x00DF}}\) \(\def \SS {\unicode {x1E9E}}\) \(\def \dag {\unicode {x2020}}\) \(\def \ddag {\unicode {x2021}}\) \(\def \P {\unicode {x00B6}}\) \(\def \copyright {\unicode {x00A9}}\) \(\def \pounds {\unicode {x00A3}}\) \(\let \LWRref \ref \) \(\renewcommand {\ref }{\ifstar \LWRref \LWRref }\) \( \newcommand {\multicolumn }[3]{#3}\) \(\require {textcomp}\) \(\newcommand {\intertext }[1]{\text {#1}\notag \\}\) \(\let \Hat \hat \) \(\let \Check \check \) \(\let \Tilde \tilde \) \(\let \Acute \acute \) \(\let \Grave \grave \) \(\let \Dot \dot \) \(\let \Ddot \ddot \) \(\let \Breve \breve \) \(\let \Bar \bar \) \(\let \Vec \vec \) \(\DeclareMathOperator {\var }{var}\) \(\DeclareMathOperator {\cov }{cov}\) \(\DeclareMathOperator {\indeg }{deg_{in}}\) \(\DeclareMathOperator {\outdeg }{deg_{out}}\) \(\newcommand {\nN }{n \in \mathbb {N}}\) \(\newcommand {\Br }{{\cal B}(\R )}\) \(\newcommand {\F }{{\cal F}}\) \(\newcommand {\ds }{\displaystyle }\) \(\newcommand {\st }{\stackrel {d}{=}}\) \(\newcommand {\uc }{\stackrel {uc}{\rightarrow }}\) \(\newcommand {\la }{\langle }\) \(\newcommand {\ra }{\rangle }\) \(\newcommand {\li }{\liminf _{n \rightarrow \infty }}\) \(\newcommand {\ls }{\limsup _{n \rightarrow \infty }}\) \(\newcommand {\limn }{\lim _{n \rightarrow \infty }}\) \(\def \ra {\Rightarrow }\) \(\def \to {\rightarrow }\) \(\def \iff {\Leftrightarrow }\) \(\def \sw {\subseteq }\) \(\def \wt {\widetilde }\) \(\def \mc {\mathcal }\) \(\def \mb {\mathbb }\) \(\def \sc {\setminus }\) \(\def \v {\textbf }\) \(\def \p {\partial }\) \(\def \E {\mb {E}}\) \(\def \P {\mb {P}}\) \(\def \R {\mb {R}}\) \(\def \C {\mb {C}}\) \(\def \N {\mb {N}}\) \(\def \Q {\mb {Q}}\) \(\def \Z {\mb {Z}}\) \(\def \B {\mb {B}}\) \(\def \~{\sim }\) \(\def \-{\,;\,}\) \(\def \|{\,|\,}\) \(\def \qed {$\blacksquare $}\) \(\def \1{\unicode {x1D7D9}}\) \(\def \cadlag {c\`{a}dl\`{a}g}\) \(\def \p {\partial }\) \(\def \l {\left }\) \(\def \r {\right }\) \(\def \F {\mc {F}}\) \(\def \G {\mc {G}}\) \(\def \H {\mc {H}}\) \(\def \Om {\Omega }\) \(\def \om {\omega }\) \(\def \Vega {\mc {V}}\)

13.2 Geometric Brownian motion

In this section we focus on a particular SDE, namely

\begin{equation} \label {eq:gbm} dX_t=\alpha X_t\, dt+\sigma X_t\, dB_t \end{equation}

where \(\alpha \in \R \) and \(\sigma \geq 0\) are deterministic constants. The parameter \(\alpha \) is known as the drift, and \(\sigma \) is known as the volatility. Equation (13.8) will be important to us because it will be our next step in establishing better models for stock prices (to be continued, in Section 15.1).

The solution to equation (13.8), which we will shortly show exists, is known as geometric Brownian motion.

The key step to solving (13.8) is to work with the logarithm of \(X\). With this aim in mind, we assume (for now) that there is a solution \(X\) that is strictly positive, and consider the process

\[Z_t=\log X_t.\]

Of course our assumption may not be true - but if such a solution does exist we hope to (do some calculations and with \(Z\) and) find an explicit formula for it, at which point we can go back and check we really do have a solution.

  • Remark 13.2.1 Taking logarithms is a natural idea to try. To see this, consider the special case \(\sigma =0\), where we find ourselves back in the world of differential equations: \(x(t)=\int _0^t \alpha x(u)\,du\). The fundamental theorem of calculus gives \(\frac {dx}{dt}=\alpha x\). We can solve this equation by considering \(z=\log x\), obtaining

    \[\frac {dz}{dt}=\frac {dz}{dx}\frac {dx}{dt}=\frac {1}{x}\alpha x=\alpha \]

    Thus \(z(t)=\alpha t+C\) and \(x(t)=e^{z(t)}=C'e^{\alpha t}\).

Using Ito’s formula, with \(Z_t=\log X_t\) (i.e. \(f(t,x)=\log x\)) we obtain

\begin{align*} dZ_t&=\l (0+\alpha X_t\frac {1}{X_t}+\frac 12(\sigma X_t)^2\frac {-1}{X_t^2}\r )\,dt + \sigma X_t \frac {1}{X_t}\,dB_t\\ &=\l (\alpha -\frac 12\sigma ^2\r )\,dt+\sigma \,dB_t. \end{align*} In integral form this gives us

\begin{align*} Z_t&=Z_0+\int _0^t\alpha -\frac {1}{2}\sigma ^2\,du+\int _0^t\sigma \,dB_u\\ &=Z_0+\l (\alpha -\frac 12\sigma ^2\r )t+\sigma B_t. \end{align*} Since \(Z_t=\log X_t\), raising both sides to the power \(e\) gives us

\begin{equation} \label {eq:gbm_soln} X_t=X_0\exp \l ((\alpha -\tfrac 12\sigma ^2)t+\sigma B_t\r ). \end{equation}

As we said, this was all based on the assumption that a (strictly positive) solution exists. But, now we have found the formula (13.9), we can go back and check that it does really give us a solution. This part is left for you, see exercise 13.12.

For future use, applying (13.9) at times \(t\) and \(T\), we obtain that

\begin{equation} \label {eq:gbm_soln_tT} X_T=X_t\exp \l ((\alpha -\tfrac 12\sigma ^2)(T-t)+\sigma (B_T-B_t)\r ). \end{equation}