Quelltext der Seite Kurs:Modellierung und Numerische Methoden von Finanzderivaten/2 Binomialmethode

== 2.1 Binomialbäume ==
Wir betrachten einen sehr einfachen Finanzmarkt, in dem mit einem Bond, einer Aktie oder einer (europäischen) Call-Option nur zu den Zeitpunkten <math>t = 0</math> und <math>t = \Delta t</math> gehandelt werden kann. Dabei werden folgende Annahmen gemacht:

* Der ''Bond'' <math>B_t</math> habe zur Zeit <math>t = 0</math> den Wert <math>B_0 = 1</math>. Der risikofreie Zinssatz betrage <math>r > 0</math> für Guthaben und für Kredite; es werde kontinuierlich verzinst, d.h. <math>B_{\Delta t} = e^{r\Delta t}</math>.

* Die ''Aktie'' habe zur Zeit <math>t = 0</math> den Wert <math>S_0 = S</math>. Zur Zeit <math>t = \Delta t</math> gebe es zwei Möglichkeiten:
** Der Kurs der Aktie betrage mit Wahrscheinlichkeit <math>q > 0</math>: <math>S_{\Delta t} = uS</math> (<math>u</math> steht für einen ”up”-Zustand) bzw.
** Der Kurs betrage mit Wahrscheinlichkeit <math>1 - q > 0</math>: <math>S_{\Delta t} = dS</math> (<math>d</math> steht für ”down”-Zustand). Es sei <math>u > d > 0</math>.

* Die ''Call-Option'' habe den Ausübungspreis <math>K</math> und die Verfallszeit <math>\Delta t</math>.
* Der Finanzmarkt sei ''arbitragefrei'' und erlaube ''Aktienleerverkäufe'' (sog. ''short selling''). Das bedeutet, dass man Aktien noch nicht besitzen muss, sondern später liefern kann. Weiter fallen keine Transaktionskosten für Kauf/Verkauf von Wertpapieren an. Es erfolgen zwischenzeitlich keine Dividendenzahlungen.

Die vierte Voraussetzung impliziert <math>d \le e^{r\Delta t} \le u</math>. Gilt nämlich <math>e^{r\Delta t} > u</math>, so kann der Kauf von Bonds durch Aktienleerverkäufe finanziert und Arbitrage erzielt werden. Falls aber <math>e^{r\Delta t} < d</math> ist, so besteht eine Arbitrage-Möglichkeit durch den Kauf von Aktien, der durch Kredite finanziert wird.

Zur Zeit <math>\Delta t</math> beträgt der Wert <math>C_{\Delta t}</math> der Call-Option entweder <math>C_u := (uS - K)^+</math> oder <math>C_d := (dS - K)^+</math>. Wir wollen den Wert <math>C_0</math> der Call-Option zur Zeit <math>t = 0</math> bestimmen. Dazu verwenden wir das '''Duplikationsprinzip'''. Wir sichern die Call-Option dadurch ab, dass wir <math>c_1</math> Anteile des Bonds und <math>c_2</math> Anteile der Aktie kaufen bzw. verkaufen. Wir suchen also ein zur Call-Option äquivalentes Portfolio, so dass gilt 
<div align="center"><math>c_1 \cdot B_t + c_2 \cdot S_t = C_t, \quad t = 0,</math> bzw. <math>t = \Delta t</math>.</div>
Die unbekannten Größen <math>c_1, c_2</math> bestimmen wir aus dieser Gleichung zum Zeitpunkt <math>t = \Delta t</math>:
<div align="center"><math>c_1e^{r\Delta t} + c_2u \cdot S = C_u,</math></div>
<div align="center"><math>c_1e^{r\Delta t} + c_2d \cdot S = C_d.</math></div>
Als Lösung ergibt sich
<div align="center"><math>c_1 = \frac{uC_d - dC_u}{(u - d)e^{r\Delta t}}, \quad c_2 = \frac{C_u - C_d}{(u - d)S}.</math></div>
Wegen <math>d < u</math> folgt <math>c_1 \le 0</math> und <math>c_2 \ge 0</math>, d. h. es muss zum Aktienkauf stets ein Kredit aufgenommen werden. Nunmehr können wir die Optionsprämie <math>C_0</math> bestimmen:
<div align="center">(2.1) <math>C_0 = c_1 \cdot e^0 + c_2  S = e^{-r\Delta t} (pC_u + (1 - p)C_d)</math> mit <math>p = \frac{e^{r\Delta t} - d}{u - d}</math>.</div>
Aus der Ungleichungskette <math>d \le e^{r\Delta t} \le u</math> folgt <math>0 \le p \le 1</math>.

Die Optionsprämie kann als ''diskontierter Erwartungswert'' bzgl. der Wahrscheinlichkeit <math>p</math> interpretiert werden. Definieren wir nämlich den Erwartungswert einer Zufallsfunktion <math>X</math>, die nur die beiden (diskreten) Zustände <math>X = X_u</math> und <math>X = X_d</math> mit Wahrscheinlichkeit <math>p</math> bzw. <math>1 - p</math> annehmen kann, durch
<div align="center"><math>E_p(X) := pX_u + (1 - p)X_d,</math></div>
so lautet die Formel (2.1) einfach
<div align="center"><math>C_0 = e^{-r\Delta t}E_p (S_{\Delta t} - K)^+.</math></div>
Der Optionspreis hängt nicht von der Wahrscheinlichkeit <math>q</math> ab; dies ist verständlich, da wir ''alle'' Pfade im Binomialbaum betrachten. Wegen der Relation
<div align="center"><math>puS + (1 - p)dS = \frac{e^{r\Delta t} - d}{u - d} uS + \frac{u - e^{r\Delta t}}{u - d} dS = e^{r\Delta t} S</math></div>
können wir <math>p</math> als die ''risikoneutrale Wahrscheinlichkeit'' interpretieren, da der erwartete Kurs des Basiswertes mit der Wahrscheinlichkeit <math>p</math> des Up-Zustandes gleich dem Erlös aus dem risikofreien Bond ist.

Wir verallgemeinern nun diese Idee, indem wir einen <math>n</math>-Perioden-Finanzmarkt betrachten. In jeder Periode der Länge <math>\Delta t</math> kann sich der Aktienkurs mit der Wahrscheinlichkeit <math>q</math> um den Faktor <math>u</math> und mit Wahrscheinlichkeit <math>1 - q</math> um den Faktor <math>d</math> ändern. Ein solches Preismodell heißt in der Literatur auch Ross-Rubinstein-Modell [2]. Seien also <math>n \in \mathbb N</math> und <math>T = n\Delta t</math> die Verfallszeit. Dann ist <math>S^n_k := u^kd^{n-k}S</math> der Wert der Aktie zur Zeit <math>T</math> bei <math>k</math> Up-Zuständen und <math>n - k</math> Down-Zuständen.

Sei zunächst <math>n = 2</math>. Dann berechnet sich der Optionspreis zur Zeit <math>t = 0</math> analog zum Ein-Perioden-Modell aus
<div align="center"><math>C_0 = e^{-r\Delta t} (pC_u + (1 - p)C_d),</math></div>
wobei <math>C_u</math> und <math>C_d</math> gegeben sind durch
<div align="center"><math>C_u = e^{-r\Delta t} (pC_{uu} + (1 - p)C_{ud}), \quad C_d = e^{-r\Delta t} (pC_{ud} + (1 - p)C_{dd}),</math></div>
und <math>C_{uu}, C_{ud}</math> und <math>C_{dd}</math> sind gegeben. Daher ergibt sich im Zwei-Perioden-Modell
<div align="center"><math>C_0 = e^{-r\Delta t} (p^2C_{uu} + 2p(1 - p)C_{ud} + (1 - p)^2C_{dd}).</math></div>
Nun kann auf die Verallgemeinerung von <math>n</math> Perioden geschlossen werden:
<div align="center">(2.2) <math>C_0 = e^{-nr\Delta t} \sum^n_{k=0} \frac{n!}{k!(n - k)!} p^k(1 - p)^{n-k} (S^n_k -K)^+.</math></div>
Unter Verwendung der Binomialkoeffizienten
<div align="center"><math>\binom n k = \frac{n!}{k!(n - k)!}</math></div>
und dem Erwartungswert einer Zufallsfunktion <math>X</math>, die mit Wahrscheinlichkeit
<div align="center"><math>\frac{n!}{k!(n - k)!} p^k(1 - p)^{n-k}</math></div>
den Wert <math>X = X_k, k = 0, \ldots, n</math> annehmen kann, durch
<div align="center">(2.3) <math>E_p(X) = \sum^n_{k=0} \binom n k p^k(1 - p)^{n-k}X_k,</math></div>
so erhalten wir die Beziehung
<div align="center"><math>C_0 = e^{-rT} E_p ((S^n - K)^+).</math></div>
Der Wert
<div align="center"><math>\binom n k p^k(1 - p)^{n-k}</math></div>
ist die Wahrscheinlichkeit, dass der Aktienkurs zur Zeit <math>T = n \Delta t</math> den Wert <math>u^kd^{n- k}S</math> besitzt.

Wir können die Formel (2.3) auch anders interpretieren. Mit dem minimalen Index
<div align="center"><math>m = \min_{0 \le k \le n} \{u^kd^{n-k} S - K\} \ge 0</math></div>
können wir (2.3) schreiben als
<div align="center"><math>C_0 = S \sum^n_{k = m} \binom n k \left( pue^{-r\Delta t} \right)^k \left( (1 - p)de^{-r\Delta t} \right)^{n - k} - Ke^{-rT} \sum^n_{k = m} \binom n k p^k (1 - p)^{n-k}.</math></div>
Sei <math>p' = pue^{-r\Delta t}</math>. Dann folgt aus der Definition von <math>p</math> die Beziehung <math>pu+(1-p)d = e^{r\Delta t}</math> oder <math>1 - p' = (1 - p)de^{-r\Delta t}</math> und daher
<div align="center">(2.4) <math>C_0 =S \sum^n_{k = m} \binom n k \left( p' \right)^k (1 - p')^{n - k} - Ke^{-rT} \sum^n_{k = m} \binom n k p^k (1 - p)^{n-k}</math></div>
<div align="center">(2.5) <math>= S \Phi(m, p') - Ke^{-rT} \Phi(m, p),</math></div>
wobei <math>\Phi</math> die Verteilungsfunktion der Binomialverteilung ist:
<div align="center"><math>\Phi(m, p) = \sum^n_{k = m} \binom n k p^k (1 - p)^{n-k}.</math></div>
Die Gleichung (2.5) kann als ''diskrete Black-Scholes-Formel'' interpretiert werden (s.u.).

== 2.2 Die Binomial-Optionspreis-Berechnung mit Mathematica ==
Wir bestimmen den auf der Binomial-Methode beruhenden Optionspreis mit dem Formelmanipulationssystem '''''Mathemtica''''' und wählen als Beispiel für den Aktien-Basispreis <math>S = 50</math>. Wir nehmen an, dass der Aktienpreis in der kommenden Periode entweder um 10% steigt oder um 3% fällt; Der (feste) Zinssatz sei <math>r = 6\%</math>. Aus diesen Werten bestimmen wir die Wahrscheinlichkeiten und schließlich den Optionspreis der Aktie.
 Solve[{55*p + 48.5*q == 50, p + q == 1/1.06}, {p, q}]
 a = Simplify[Solve[{p*up + q*down == 1, p + q == 1/R}, {p, q}]]

Weiter bestimmen wir den (Zustands-) Optionspreis über <math>n</math> Perioden:
 Clear[statePrices]
 statePrices[up_, down_, R_] := Solve[{p*up + q*down == 1, p + q == 1/R}, {p, q}][[1]]
 Clear[binomialCall]
 binomialCall[s_, x_, n_] := Sum[p^j*q^(n - j)*Binomial[n, j]*Max[s*p^j*down^(n - j) - x, 0], {j, 0, n}]
 statePrices[up, down, R];

Nun sind die Parameter der Binomial-Verteilung für ein <math>n</math>-Perioden-Modell zu bestimmen. Dabei verwenden wir die asymptotische Beziehung und setzen für <math>S_t/S</math> die entsprechenden Werte ein:
<div align="center">(2.6) <math>\log \left( \frac{S_t}{S} \right) \sim \mathcal N(\mu t, \sigma \sqrt{t}),</math></div>
<div align="center"><math>\log \left( \frac{Su^jd^{n - j}}{S} \right) = j \log u + (n - j) \log d = j \log(u/d) + n \log(d),</math></div>
<div align="center">(2.7) <math>\operatorname{E}\, \left\{ \left( \frac{Su^jd^{n - j}}{S} \right) \right\} = \operatorname{E}(j) \log(u/d) + n \log(d) = np \log(u/d) + n \log d,</math></div>
<div align="center">(2.8) <math>\operatorname{Var}\, \left\{ \left( \frac{Su^jd^{n - j}}{S} \right) \right\} = np (1 - p) \left[ \log \left( \frac{u}{d} \right) \right]^2.</math></div>
Die Beziehung für die Varianz folgt aus der Tatsache, dass in jeder Periode die Gleichung
<div align="center"><math>(u - d)^2 \left( p(1 - p)^2 + (1 - p)p^2 \right) = (u - d)^2 p (1 - p)</math></div>
gilt. Damit gelten für den ''Erwartungswert'' <math>\mu</math> und für die ''Varianz (Volatilität)'' <math>\sigma</math> die
beiden Gleichungen
<div align="center"><math>np \log(u/d) + n \log d = \mu, \quad np (1 - p) \left[ \log \left( \frac{u}{d} \right) \right]^2 = \sigma^2.</math></div>

Wir realisieren die Formeln in '''''Mathematica''''' und erhalten mit der Anweisungsfolge zunächst allgemeine Formeln, in die konkrete Werte eingesetzt werden: Somit sind folgende Anweisungen zu realisieren:
 Clear[u, d, mu]
 a = Solve[{(p*Log[u/d] + Log[d])*n == mu, n*p*(1 - p)*Log[u/d]^2 == var}/.p->1/2, {u, d}]
 a[[4]]/.(mu->0.12, var->0.2^2, p->0.5, n->100)

Die Werte <math>p = 0.5, \mu = 0.12, \sigma^2 = 0.2^2</math> (d.h. var), <math>n = 100</math> liefern folgenden Plot:
 p = 0.5; 
 mu = 0.12;
 var = 0.2^2;
 n = 100;
 aa = 
  ListPlot[
   Table[{u^j*d^(n - j), p^j*(1 - p)^(n - j)*Binomial[n, j]}, {j, 0, n}]/.a[[4]], 
   Joined->True, PlotRange->A11]

==== Bemerkung: ====
Das ist nicht die einzige Lösungsmenge, die gegen die Log-Normal-Verteilung konvergiert. Cox, Ross, Rubinstein [2] benutzen die Werte
<div align="center"><math>u = \exp \left( \sigma \sqrt{1/n} \right), \quad d = \frac{1}{u} = \exp \left( - \sigma \sqrt{1/n} \right), \quad p = \frac{1}{2} + \frac{1}{2} \frac{\mu}{\sigma} \sqrt{1/n}.</math></div>
Diese Werte lösen die Gleichungen zwar nicht exakt, sie konvergieren aber für <math>n \to \infty</math> gegen die exakten Werte.

== 2.3 Brownsche Bewegung und ein Aktienkursmodell ==
Wir stellen zunächst einige Begriffe aus der Stochastik zusammen.

Wir bezeichnen mit <math>(\Omega, \mathcal F, P)</math> einen Wahrscheinlichkeitsraum. <math>\Omega</math> bezeichnet hier im einfachsten Fall eine endliche (oder - etwas komplizierter - eine abzählbare) Menge, die Menge der Elementarereignisse, <math>\mathcal F</math> ist die Klasse aller Teilmengen von <math>\Omega</math> (d.h. die Klasse von Ereignissen), <math>P(A)</math> wird die Wahrscheinlichkeit des Ereignisses <math>A \in \mathcal F</math> genannt.
==== Bemerkung: ====
(a) Die Klasse <math>\mathcal F</math> ist eine <math>\sigma</math>-Algebra, d.h.
# <math>A \in \mathcal F</math> impliziert <math>\overline A = \Omega \setminus A \in \mathcal F</math>,
# <math>A_k \in \mathcal F, k = 1, 2, \ldots</math> impliziert <math>\bigcup A_k \in \mathcal F</math> und <math>\bigcap A_k \in \mathcal F</math>.
(b) Die numerische, auf <math>\mathcal F</math> definierte Funktion <math>P</math> ist
# positiv semidefinit, d.h. <math>P(A) \ge 0 \quad \forall A \in \mathcal F</math>,
# normalisiert, d.h. <math>P(\Omega) = 1</math>,
# <math>\sigma</math>-additiv, d.h. <math>P\left( \bigcup_k A_k \right) = \sum_k P(A_k)</math> für <math>A_k, k \in \mathbb N</math>, paarweise disjunkt.

Im allgemeinen Fall, wo <math>\Omega</math> z. B. überabzählbar ist, gilt die Forderung, dass die Klasse <math>\mathcal F</math> eine <math>\sigma</math>-Algebra ist, ebenfalls. Zur Definition und für Beispiele informiere man sich in Grundlagenwerken der Stochastik, siehe z. B. [10], oder man vgl. die Aussagen im Kurs Stochastik.

=== Definition 2.1 ===
:'''''(1)''' Ein Wahrscheinlichkeitsraum ist ein Tripel <math>(\Omega, \mathcal F, P)</math>, bestehend aus einer Menge <math>\Omega</math>, einer <math>\sigma</math>-Algebra <math>\mathcal F</math> und einem Wahrscheinlichkeitsmaß <math>P</math>.
:'''''(2)''' Eine Funktion <math>X: \Omega \to \mathbb R</math> auf einem Wahrscheinlichkeitsraum <math>(\Omega, \mathcal F, P)</math> heißt Zufallsvariable genau dann, wenn für alle Borelmengen <math>B \in B(\mathbb R)</math> gilt: <math>X^{-1}(B) \in \mathcal F</math>.

==== Bemerkung: ====
Die Menge <math>B(\mathbb R)</math> aller Borelmengen ist die kleinste <math>\sigma</math>-Algebra, die alle offenen Mengen von <math>\mathbb R</math> enthält.

=== Definition 2.2 ===
:''Sei <math>(\Omega, \mathcal F, P)</math> ein Wahrscheinlichkeitsraum. Die Zufallsvariablen <math>X_1, \ldots, X_n</math> heißen unabhängig, falls <math>\forall A_k \in B(\mathbb R), k = 1, \ldots, n</math> gilt
<div align="center"><math>P \Bigl( \{X_1 \in A_1\} \cap \ldots \cap \{X_n \in A_n\} \Bigr) = P \Bigr( \{X_1 \in A_1\} \Bigr) \cdot \ldots \cdot P \Bigl( \{X_n \in A_n\} \Bigr).</math></div>

Wir gehen nun zu stetigen Zufallsvariablen über. Unter gewissen Voraussetzungen können Erwartungswert und Varianz mittels eines Riemann-Integrals berechnet werden.

=== Definition 2.3 ===
:''1. Sei <math>(\Omega, \mathcal F, P)</math> ein Wahrscheinlichkeitsraum und <math>X</math> eine Zufallsvariable. Dann ist der Erwartungswert <math>E(X)</math> definiert durch
<div align="center"><math>E(X) := \int\limits_\Omega X(\omega)\, dP, \quad \omega \in \Omega.</math></div>
:''(Das ist das Lebesgue-Integral der summierbaren Funktion <math>X: \Omega \to \mathbb R</math> bzgl. des Maßes <math>P</math>.)
:''2. Die Varianz <math>\operatorname{Var}(X)</math> ist definiert durch
<div align="center"><math>\operatorname{Var}(X) := E \Bigl( \bigl( X - E(X) \bigr)^2 \Bigr) = E(X^2) - (E(X))^2,</math></div>
:''falls die entsprechenden Integrale existieren. Die Standardabweichung ist gegeben durch <math>\sigma := \sqrt{\operatorname{Var}(X)}</math>.
:''3. Die Kovarianz zweier Zufallsvariabler X, Y</math> ist definiert durch
<div align="center"><math>\operatorname{Cov}(X, Y) := E \Bigl( (X - E(X))(Y - E(Y )) \Bigr) = E(XY) - E(X)E(Y),</math></div>
:''falls <math>E(X), E(XY)</math> und <math>E(Y)</math> existieren.

Für unabhängige Zufallsvariable <math>X</math> und <math>Y</math> gilt <math>E(XY) = E(X)E(Y)</math> und insbesondere <math>\operatorname{Cov}(X, Y) = 0</math>. Daher kann die Kovarianz als ein Maß für die Abhängigkeit von <math>X</math> und <math>Y</math> betrachtet werden. Eine weitere Konsequenz aus der Unabhängigkeit von <math>X</math> und <math>Y</math> ist die Additivität der Varianzen:
<div align="center"><math>(2.9) \operatorname{Var}(X + Y ) = E(X^2) + 2E(X)E(Y) + E(Y^2) \Bigl( E(X) + E(Y) \Bigr)^2 = \operatorname{Var}(X) + \operatorname{Var}(Y).</math></div>

Die Integrale werden i. a. über ihre Darstellung der Zufallsvariablen mittels charakteristischer Funktionen definiert. Sind die Funktionen Riemann-integrabel, so erhalten wir folgende Darstellungen:

=== Definition 2.4 ===
:''Sei <math>X: \Omega \to \mathbb R</math> eine Zufallsvariable auf einem Wahrscheinlichkeitsraum <math>(\Omega, \mathcal F, P)</math>.
:''1. Die Verteilungsfunktion <math>F</math> von <math>X</math> ist definiert durch
<div align="center"><math>F(x) = P(X < x) = P(\{y \in \Omega: X(y) < x\}).</math></div>
:''Die Dichtefunktion <math>f</math> von <math>X</math> (falls sie existiert) ist gegeben durch
<div align="center">''<math>F(x) = \int\limits_{- \infty}^x f(s)\, ds</math> oder <math>F' = f</math>.</div>

:''2. Der Erwartungswert und die Varianz einer Zufallsvariablen <math>X</math>, die eine Dichtefunktion besitzt, sind gegeben durch
<div align="center"><math>E(X) = \int\limits_\mathbb{R} xf(x)\, dx, \quad \operatorname{Var}(X) = E(X^2) - E(X)^2 = \int\limits_\mathbb{R} (x - E(X))^2f(x)\, dx,</math></div>
:''falls die entsprechenden Integrale existieren.

Wir geben folgende Beispiele an:

'''(1)''' Betrachten <math>\Omega = \mathbb R, \mu \in \mathbb R</math> und <math>\sigma > 0</math>. Eine Zufallsvariable <math>X</math> mit der Dichtefunktion
<div align="center"><math>f(x) = \frac{1}{\sqrt{2\pi \sigma^2}} \exp \left( - \frac{(x - \mu)^2}{2\sigma^2} \right), \quad x \in \mathbb R</math></div>
heißt ''normalverteilt'' oder <math>N(\mu, \sigma^2)</math>-verteilt. Der Erwartungswert und die Varianz lauten <math>E(X) = \mu</math> und <math>\operatorname{Var}(X) = \sigma^2</math>. Eine <math>N(0, 1)</math>-verteilte Zufallsvariable nennen wir ''standardnormalverteilt''. Für die Verteilungsfunktion <math>F</math> einer standardnormalverteilten Zufallsvariable schreiben wir auch
<div align="center"><math>\Phi(x) = \frac{1}{\sqrt{2\pi}} \int\limits_{- \infty}^x e^{-s^2/2}\, ds.</math></div>
Eine Zufallsvariable <math>X</math>, deren Logarithmus normalverteilt ist, d. h. für die gilt, dass <math>\ln X</math> normalvertelit ist, heißt ''lognormalverteilt''.

'''(2)''' Wir betrachten den Kurs einer Aktie in einem zeit-diskreten Finanzmarkt. Damit kann der Aktienkurs um einen konstanten Faktor mit Wahrscheinlichkeit <math>q</math> steigen (”up”) und mit Wahrscheinlichkeit <math>1 - q</math> fallen (”down”). Die Zufallsvariable <math>X</math> beschreibe die Anzahl der ”up”-Zustände in einem <math>n</math>-Perioden-Finanzmarkt. Dann ist der Wertebereich von <math>X</math> die Menge <math>\{0, 1, 2, \ldots, n\}</math>, und die Wahrscheinlichkeit, dass der Aktienkurs <math>k</math>-mal steigt, beträgt
<div align="center"><math>P(X = k) = \binom n k q^k (1 - q)^{n - k}.</math></div>
Wir sagen, dass <math>X</math> ''binomialverteilt'' oder <math>B(n, q)</math>-verteilt ist. Der Aktienkurs lautet dann (in den gleichen Notationen, wie im 1. Teil angegeben)
<div align="center"><math>S^n_X = u^X d^{n - X} S.</math></div>
Die Verteilungsfunktion ist dann gegeben durch
<div align="center"><math>F(m) = P(X \le m) = \sum^m_{k=0} P(X = k) = \sum^m_{k=0} \binom n k q^k (1 - q)^{n - k}.</math></div>

Ersetzen wir in der Definition der Verteilungsfunktion das Integral durch eine Summe,
so erkennen wir, dass die Dichtefunktion gerade <math>f(k) = P(X = k)</math> ist. Damit sind der
Erwartungswert und die Varianz gegeben durch
<div align="center"><math>E(X) = \sum^n_{k=0} k f(k) = nq, \quad \operatorname{Var}(X) = \sum^n_{k=0} (k - E(X))^2 f(k) = nq(1 - q).</math></div>
Damit können wir die diskrete Formel (2.5) zur Berechnung der Optionsprämie <math>C_0</math> in der folgenden Form schreiben:
<div align="center"><math>C_0 = S P(X_{p'} \ge m) - Ke^{-rT}P(X_p \ge m),</math></div>
wobei <math>X_{p'}</math> bzw. <math>X_p</math> jetzt <math>B(n, p')</math>- bzw. <math>B(n, p)</math>-verteilte Zufallsvariablen sind.

'''Stochastische Prozesse und Brownsche Bewegung.''' Der zeitstetige Aktienkurs <math>S_t</math> ist eine Zufallsvariable, d.h. wir sollten eigentlich <math>S(\omega, t)</math> schreiben, wobei <math>\omega</math> ein Element des zugrundeliegenden Wahrscheinlichkeitsraumes ist. Man lässt jedoch i. a. das Argument <math>\omega</math> weg. Die Funktion <math>\omega \to S(\omega, t)</math> ist eine Zufallsvariable für festes <math>t</math>. Welche Regularität besitzt nun <math>t \to S(\omega, t)</math> für festes <math>\omega</math>? Die Klärung dieses Zusammenhangs führt auf den Begriff des stochastischen Prozesses.

=== Definition 2.5 ===
:''Ein (stetiger) stochastischer Prozess <math>X_t, t \in [0,\infty)</math>, ist eine Familie von Zufallsvariablen <math>X: \Omega \times [0,\infty) \to \mathbb R</math>, wobei <math>t \to X(\omega, t)</math> stetig ist für (fast alle) <math>\omega \in \Omega</math>.
:''Wir schreiben <math>X_t = X(t) = X(\cdot, t)</math>, d.h. <math>X_t</math> ist eine Zufallsvariable.

Ein stochastischer Prozess ist also eine funktionenwertige Zufallsvariable. Ein spezieller stochastischer Prozess ist durch die sog. ''Brownsche Bewegung'' gegeben: <br />
Sei <math>Z_k</math> eine Folge von unabhängigen Zufallsvariablen, die mit jeweils mit Wahrscheinlichkeit <math>\frac{1}{2}</math> die Werte <math>+\sqrt{\Delta t}</math> oder <math>-\sqrt{\Delta t}</math> annehmen und es sei weiter
<div align="center"><math>X^n_t = \sum^n_{k = 1} Z_k</math>, wobei <math>n \Delta t = t</math></div>
ist. Dann gilt für alle <math>k</math> für Erwartungswert und Varianz <math>E(Z_k) = 0, \operatorname{Var}(Z_k) = \Delta t</math> und wegen (2.9)
<div align="center"><math>E(X^n_t) = 0, Var(X^n_t) = \Bigl( E(X^n_t)^2 \Bigr) = nE(Z^2_1) = n \operatorname{Var}(Z_1) = n\Delta t = t.</math></div>
Man kann zeigen, dass <math>X^n_t</math> für <math>n \to \infty</math> (und damit für <math>\Delta t \to 0</math>, so dass <math>n\Delta t = t</math> gilt) konvergiert. Den Grenzwert nennen wir die Brownsche Bewegung <math>W_t</math>.
=== Satz 2.1 (Wiener) ===
:''Es gibt einen stetigen stochastischen Prozess <math>W_t</math> mit den folgenden Eigenschaften:
:# <math>W_0 = 0,</math>
:# ''Für alle <math>0 \le s \le t</math> gilt: <math>W_t - W_s</math> ist <math>N(0, t - s)</math>-verteilt.
:# ''Für alle <math>0 \le r \le u \le s \le t</math> gilt: W_t - W_s</math> und <math>W_u - W_r</math> sind unabhängig.

Der stochastische Prozess <math>W_t</math> aus Satz 2.1 wird ''Wiener Prozess'' oder ''Brownsche Bewegung'' genannt. Aus Satz 2.1, (2) folgt unmittelbar, dass <math>W_t</math> <math>N(0, t)</math>-verteilt ist und dass gilt
<div align="center">(2.10) <math>E \Bigl( (W_t - W_s)^2 \Bigr) = t - s</math></div>
Der Beweis von Satz 2.1 ist aufwändig.

'''Ein Aktienkursmodell.''' Mit Hilfe des Wiener-Prozesses können wir ein Modell für die Entwicklung von Aktienkursen angeben. Wir betrachten zunächst einen Bond <math>B_t</math> mit dem risikofreien Zinssatz <math>r \ge 0</math>. Dann gilt <math>B_t = B_0e^{rt}</math> oder
<div align="center"><math>\ln B_t = \ln B_0 + r \cdot t.</math></div>
Das legt den folgenden Ansatz für den Aktienkurs <math>S_t</math> nahe:
<div align="center"><math>\ln S_t = \ln S_0 + b \cdot t + \text{Zufallseinfluss}.</math></div>
Da in einem arbitragefreien Markt die Anlage einer Aktie die gleiche Rendite wie ein Bond bringen soll, ersetzen wir <math>r \cdot t</math> durch <math>b \cdot t</math>. ein solches Verhalten kann dem (längerfristigen) Verlauf von Aktienkursen entnommen werden. Er verändert sich mit dem Driftterm und wird durch zufällige Schwankungen überlagert.

Für den Zufall nehmen wir an, dass er ohne Tendenz ist, d.h. der Erwartungswert der (zufälligen) Schwankungen sei Null; er hänge jedoch von der Zeit ab. Das wird durch den Ansatz
<div align="center"><math>\ln S_t = \ln S_0 + b_t + \sigma W_t</math></div>
erfüllt. Man definiert üblicherweise <math>\mu = b + \frac{1}{2} \sigma^2</math>, das impliziert
<div align="center">(2.11) <math>\ln S_t = \ln S_0 + \left( \mu - \frac{1}{2} \sigma^2 \right) t + \sigma W_t.</math></div>
Der Aktienkurs ist dann gegeben durch
<div align="center">(2.12) <math>S_t = S_0 \exp \left( \mu t + \sigma W_t - \frac{1}{2} \sigma^2t \right).</math></div>
Man nennt <math>S_t</math> eine ''geometrische Brownsche Bewegung''; <math>S_t</math> ist lognormalverteilt. Bei der Analyse der Black-Scholes-Formel kommen wir auf die Beziehungen zurück.

Für den Erwartungswert und die Varianz von <math>S_t</math> gelten nun folgende Aussagen:

=== Lemma 2.1 ===
:''Sei <math>S_t</math> eine geometrische Brownsche Bewegung. Dann gilt
<div align="center">(2.13) <math>E(S_t) = S_0e^{\mu t},</math></div>
<div align="center">(2.14) <math>\operatorname{Var}(S_t) = S^2_0 \exp(2\mu t) \Bigl( e^{\sigma^2 t} - 1 \Bigr).</math></div>

==== Beweis: ====
Wir verwenden die Eigenschaft, dass für eine normalverteilte Zufallsgröße <math>X</math> gilt, dass <math>\exp(X)</math> logarithmisch normalverteilt (<math>\log</math>-normalverteilt) ist. Da <math>W_t</math> <math>N(0, t)</math>-verteilt ist, folgt
<div align="center"><math>E(e^{\sigma W_t}) = \frac{1}{\sqrt{2\pi} t} \int\limits_\mathbb{R} e^{\sigma x} e^{- \frac{x^2}{2t}}\, dx = \frac{1}{\sqrt{2\pi} t} e^\frac{\sigma^2 t}{2} \int\limits_\mathbb{R} e^{- \frac{(x - \sigma t)^2}{2t}}\, dx = \exp \left( \frac{\sigma^2 t}{2} \right),</math></div>
also
<div align="center"><math>E(S_t) = S_0 \exp \left( \mu t - \frac{\sigma^2 t}{2} \right) E(e^{\sigma W_t}) = S_0 e^{\mu t}</math></div>
und
<div align="center"><math>\operatorname{Var}(S_t) = E(S^2_t) - E(S_t)^2 = S^2_0 \exp((2\mu - \sigma^2) t) E(e^{2\sigma W_t}) - S^2_0 e^{2\mu t} = S^2_0 e^{2\mu t} (e^{\sigma^2 t} - 1).</math></div>

q.e.d. 

Der Aktienkurs ist nach Gleichung (2.14) ein Produkt aus dem mittleren Kurs <math>S_0e^{\mu t}</math> und einer zufälligen Schwankung <math>\exp \left( \sigma W_t - \frac{1}{2} \sigma^2 t \right)</math> um diesen mittleren Kurs. Die Formeln (2.14) benötigen wir für die Betrachtungen der Binomialmethode im nächsten Abschnitt.

== 2.4 Der Algorithmus der Binomialmethode ==
Wir sind nun in der Lage, einen ersten numerischen Algorithmus vorzustellen, um den Preis einer europäischen oder einer amerikanischen Option approximativ zu bestimmen: <br />
Wir zerlegen das Zeitintervall <math>[0, T]</math> in <math>N</math> Zeitschritte der Länge <math>\Delta t = T/N</math> und suchen Approximationen <math>S_k</math> von <math>S(t_k)</math> zu den Zeitpunkten <math>t_k = k\Delta t, k = 0, 1, \ldots,N</math>. Wir treffen folgende Annahmen:
:Der Kurs <math>S_k</math> nimmt nach Ablauf der Zeit <math>\Delta t</math> entweder mit Wahrscheinlichkeit <math>p</math> den Wert <math>S_{k+1} = uS_k</math> (”up”) oder mit Wahrscheinlichkeit <math>1 - p</math> den Wert <math>S_{k+1} = dS_k</math> (”down”) an.
:Die erwartete Rendite im Zeitraum <math>\Delta t</math> entspricht dem risikolosen Zinssatz, d. h. wir setzen <math>\mu = r</math> in (2.14):
<div align="center">(2.15) <math>E(S(t_{k+1})) = S(t_k) e^{r\Delta t}, \quad \operatorname{Var}(S(t_{k+1})) = \Bigl( S(t_k) \Bigr)^2 e^{2r \Delta t} \Bigl( e^{\sigma^2 \Delta t} - 1 \Bigr).</math></div>
:Für den Optionspreis <math>V(t_k)</math> gelte analog
<div align="center">(2.16) <math>E(V(t_{k+1})) = V(t_k) e^{r\Delta t}.</math></div>
:Für den Kauf oder Verkauf von Wertpapieren fallen keine Transaktionskosten an, es erfolgen keine Dividendenzahlungen.

Die Parameter <math>u, d, p</math> sind (zunächst) unbekannt.

Die Grundidee der Binomialmethode ist es, die Erwartungswerte und Varianzen des (zeit-)kontinuierlichen und des (zeit-)diskreten Modells gleichzusetzen. Für das diskrete Modell gilt
<div align="center">(2.17) <math>E(S_{k+1}) = p \cdot uS_k + (1 - p) \cdot dS_k,</math></div>
<div align="center">(2.18) <math>\operatorname{Var}(S_{k+1}) = E(S^2_{k+1}) - \Bigl( E(S_{k+1}) \Bigr)^2</math></div>
<div align="center">(2.19) <math>= p(uS_k)^2 + (1 - p)(dS_k)^2 - \Bigl( puS_k + (1 - p)dS_k \Bigr)^2.</math></div>

Ersetzen wir in (2.15) <math>S(t_k)</math> durch die Approximation <math>S_k</math> und verwenden wir (2.19), so erhalten wir
<div align="center"><math>S_ke^{r\Delta t} = puSk + (1 - p)dS_k,</math></div>
<div align="center"><math>S^2_k e^{2r \Delta t} \Bigl( e^{\sigma^2 \Delta t} - 1 \Bigr) = p(uS_k)^2 + (1 - p)(dS_k)^2 - \Bigl( puS_k + (1 - p)dS_k \Bigr)^2.</math></div>
Die erste Gleichung ist äquivalent zu
<div align="center">(2.20) <math>e^{r \Delta t} = pu + (1 - p)d.</math></div>
Wir setzen diese Beziehung in die zweite Gleichung ein und erhalten
<div align="center"><math>e^{2r \Delta t} \Bigl( e^{\sigma^2 \Delta t} - 1 \Bigr) = pu^2 + (1 - p)d^2 - \Bigl( e^{r \Delta t} \Bigr)^2</math></div>
oder die dazu äquivalente Gleichung
<div align="center">(2.21) <math>e^{2r \Delta t+\sigma^2 \Delta t} = pu^2 + (1 - p)d^2.</math></div>
Damit haben wir mit (2.20) und (2.21) zwei Gleichungen für die drei Unbekannten <math>u, d, p</math>. Wir wählen zusätzlich dazu eine Symmetriebeziehung für einen Kursanstieg bzw. einen Kursabfall:
<div align="center">(2.22) <math>u \cdot d = 1.</math></div>
Die Wahl von <math>p = \frac{1}{2}</math> wäre auch möglich. Damit ist noch das nichtlineare Gleichungssystem (2.20) bis (2.22) zu lösen:
<div align="center"><math>u = \beta + \sqrt{\beta^2 - 1}, \quad d = \beta - \sqrt{\beta^2 - 1}, \quad p = \frac{e^{r \Delta t} - d}{u - d},</math></div>
wobei für <math>\beta</math> gilt
<div align="center"><math>\beta = \frac{1}{2} \Bigl( e^{- r \Delta t} + e^{(r + \sigma^2) \Delta t} \Bigr).</math></div>
Die Binomialmethode wird nun in drei Schritten realisiert: <br />
Sei <math>S_0</math> der Aktienkurs zur Zeit <math>t = 0</math> und seien <math>S_{lk} = u^ld^{k - l} S_0</math> für <math>k = 0, \ldots, N</math> und <math>l = 0, \ldots, k</math> die möglichen Aktienkurse zur Zeit <math>t_k</math>.

:'''1. Schritt:''' Initialisierung des Binomialbaumes. Berechne (für den Fall europäischer oder amerikanischer Optionen) <math>\forall l = 0, \ldots, N</math>:
<div align="center"><math>S_{lN} = u^ld^{N - l} S_0.</math></div>

:'''2. Schritt:''' Berechne die Optionswerte <math>V</math>. Für <math>t = T</math> ist <math>V(S, T)</math> durch die Endbedingung bekannt. Zu berechnen ist folglich für <math>l = 0, \ldots, N</math>:
<div align="center"><math>V_{lN} = \begin{cases} (S_{lN} - K)^+ & \text{(Europ.) Call} \\ (K - S_{lN}) & \text{(Europ.) Put} \end{cases}.</math></div>
:Bei amerikanischen Optionen ist hier zu überprüfen, ob ein vorzeitiger Kauf/Verkauf günstiger ist.
Den dritten Schritt der Berechnung der Optionspreise <math>V_{lk}</math> nehmen wir wie folgt vor:<br />
Wir können (2.20) mit Hilfe der Definition von <math>S_{lk}</math> auch in der Form
<div align="center"><math>S_{lk} e^{r \Delta t} = puS_{lk} + (1 - p)dS_{lk} = pS_{l+1,k+1} + (1 - p)S_{l,k+1}</math></div>
schreiben. Ersetzen wir in (2.16) den zeit-kontinuierlichen Optionspreis durch den zeitdiskreten, d. h. <math>V_k e^{r \Delta t} = E(V_{k+1})</math>, so folgt für den Aktienkurs die Gleichung
<div align="center"><math>V_{lk} e^{r \Delta t} = pV_{l+1,k+1} + (1 - p)V_{l,k+1}.</math></div>

:'''3. Schritt:''' (''Rückwärtsiteration'')
:Für alle <math>k = N-1, \ldots, 0</math> und für <math>l = 0, \ldots, k</math> bestimme man im Falle ''europäischer Optionen''
<div align="center"><math>V_{lk} = e^{- r \Delta t} \Bigl( pV_{l+1,k+1} + (1 - p)V_{l,k+1} \Bigr)</math></div>
:und im Falle ''amerikanischer Optionen''
<div align="center"><math>\tilde V_{lk} = e^{- r \Delta t} \Bigl( pV_{l+1,k+1} + (1 - p)V_{l,k+1} \Bigr)</math></div>
<div align="center"><math>V_{lk} = \max \{V_{lN}, \tilde V_{lk}\}.</math></div>

Der Wert <math>V_{00}</math> ist eine Approximation der Optionsprämie <math>V(S_0, 0)</math>.

Die angegebenen Schritte können mittels '''''Mathematica''''' realisiert werden. Wie bei dem nachfolgenden Matlab-Programm wählen wir den Basiswert <math>S_0 = 5</math>, die Zinsrate <math>r = 4\%</math>, die Volatilität <math>\sigma = 0.3</math>, die Periodendauer <math>T = 1</math> und greifen z. B. auf ein Fünf-Perioden-Modell zurück.

Wir realisieren die ''payoff''-Funktion
<div align="center"><math>\max \{S_T - X, 0\}, \quad \max \{X - S_T, 0\}</math></div>
in '''''Mathematica''''' durch die Funktionen
 Clear[payoffCall, payoffPut];
 payoffCall[s_] := Max[s - X, 0];
 payoffPut[s_] := Max[X - s, 0];
in den benötigten Varianten für Call- und Put-Option. Das geschieht durch jeweiliges Vertauschen der beiden Summanden. Wir bilden weiter den Binomialbaum und berechnen mit Hilfe einer Prozedur die Optionswerte.
 Clear[r, n, S0, h1, up, down, p]
 {* up=1.1; down=0.95; *}
 r = 0.04;
 sigma = 0.3; 
 n = 4; 
 S0 = 5; 
 h1 = 1/2*(Exp[-r/(n + 1)] + Exp[(r + sigma^2)/(n + 1)]);
 up = h1 + Sqrt[h1^2 - 1];
 down = 1/up;
 p0 = (Exp[r/(n + 1)] - down)/(up - down);
 Print["up = ", up]
 Print["down = ", down]

 stock = Table[5*up^(j - 1)*down^(i - j), {i, 1, 5}, {j, 1, i}];
 MatrixForm[stock]
 
 solution = Solve[{p*(1 + r) + q*(1 + r) == 1, up*p + down*q == 1}{p, q}]
Die Zinsrate gilt für jede Periode (was hier mit <math>T = 1</math> zusammenfällt!). Der gesamte Baum für die errechneten Werte <math>p</math> und <math>q</math> hat dann folgende Gestalt:
 n = 4;
 p = solution[[1, 1, 2]];
 q = solution[[1, 2, 2]];
 statePrices = Table[p^(j - 1)*q^(i - j), {i, 1, n + 1}, {j, 1, i}];
 MatrixForm[statePrices]
mit dem (aktuellen, d. h. dem gegenwärtigen) Optionspreis für <math>X = 6</math> von
 X = 6;
 statePrices[[n + 1]] . payoffPut /@ stock[[n + 1]]
 statePrices[[n + 1]] . Map[payoffPut, stock[[n + 1]]];
 {* --> alternativer Befehl *}
Ein anderer Weg zur Berechnung des Optionspreises ist die Konstruktion eines intermediären Binomial-Baumes, wo auch Dividenden-Ausschüttungen, Änderungen von Zinssätzen oder Volatilitäten usw. eingebaut werden können, durch eine Prozedur, die
sog. ''Europäische Optionen'' berücksichtigt.
{| class="wikitable sortable"
|- class="hintergrundfarbe5"
! N  || V(5,0) || CPU-Zeit [sec]
|-
|10
| 1.096558
|< 0.001
|-
| 100
| 1.092326
| 0.06
|-
| 500
| 1.094202 
| 1.52
|-
| 1000
| 1.094364
| 5.78
|-
| 2000 
| 1.094361 
| 24.8
|}
Die Übereinstimmung mit dem weiter unten durch eine MATLAB-Prozedur berechneten Put-Wert ist (bei Beachtung der jeweiligen Unterteilungszahl in Perioden) gegeben.

Als weiteres Beispiel ist ein Programm (in Matlab) zur Berechnung einer europäischen Put-Option mit <math>K = 6, r = 0.04, \sigma = 0.3, T = 1</math> angegeben. Mit Hilfe der (weiter unten betrachteten) Black-Scholes-Gleichung erhält man als ”exakte” Lösung <math>V(5, 0) = 1.09435</math>; für die Binomialmethode sind für verschiedene Zerlegungen N die zeit-diskreten Lösungen in der Tabelle oben angegeben.

==== Bemerkung: ====
1. In Abschnitt 2.1 wurde begründet, dass die Arbitrage-Freiheit des Marktes die Ungleichungskette <math>d \le e^{r\Delta t} \le u</math> nach sich zieht. Diese Ungleichungekette ist erfüllt, falls <math>\sigma \ge 0</math> ist, denn es gilt
<div align="center"><math>\beta - \sqrt{\beta^2 - 1} \le e^{r\Delta t} \le \beta + \sqrt{\beta^2 - 1} \Leftrightarrow |e^{r\Delta t} - \beta| \le \sqrt{\beta^2 - 1} \Leftrightarrow (e^{r\Delta t} - \beta)^2 \le \beta^2 - 1 \Leftrightarrow \beta \ge \frac{1}{2} (e^{r\Delta t} + e^{- r\Delta t}).</math></div>
Die letzte Ungleichung ist für alls <math>r \ge 0, \sigma \ge 0, \Delta t > 0</math> erfüllt, was in praktischen Fällen zutreffend ist.

2. Eine alternative Binomialmethode erhalten wir durch Lösung der Gleichungen (2.20), (2.21) und <math>p = \frac{1}{2}</math>. Dieses Vorgehen liefert (Man versuche, den Beweis selbständig zu führen!)
<div align="center"><math>d = e^{r\Delta t} \left( 1 - \sqrt{e^{\sigma^2 \Delta t} - 1} \right), \quad u = e^{r\Delta t} \left( 1 +\sqrt{e^{\sigma^2 \Delta t} - 1} \right), \quad p = \frac{1}{2}.</math></div>
In diesem Fall können die Parameter <math>r, \sigma, \Delta t</math> nicht beliebig sein, da sonst <math>d < 0</math> möglich ist. Dies wird verhindert, wenn wir den Zeitschritt klein genug wählen, nämlich
<div align="center"><math>\Delta t \le \frac{\ln 2}{\sigma^2}.</math></div>
Hier findet sich ein ''Stabilitätsproblem'' der Numerik wieder, das wir aber momentan nicht näher beleuchten wollen.

3. Die beschriebene Binomialmethode basiert darauf, Erwartungswerte und Varianzen des zeit-kontinuierlichen und des zeit-diskreten Modells gleich zu setzen und weiter <math>u \cdot d = 1</math> oder <math>p = \frac{1}{2}</math> vorauszusetzen. Eine weitere Variante erhält man, indem nicht nur die Erwartungswerte und die Varianzen, sondern auch die ''dritten Momente'', d. h. Ausdrücke <math>E(S_k)^3</math>, gleichgesetzt werden. Das führt auf die Gleichungen (2.20), (2.21) und auf
<div align="center"><math>p(S_k)^3 + (1 - p)(dS_k)^3 = S^3_k e^{3(r+\sigma^2) \Delta t}.</math></div>
Die Lösung lautet dann
<div align="center"><math>d = \frac{Q}{2} e^{r \Delta t} \left( Q + 1 - \sqrt{Q^2 + 2Q - 3} \right), \quad u = \frac{Q}{2} e^{r \Delta t} \left( Q + 1 + \sqrt{Q^2 + 2Q - 3} \right),</math></div>
<div align="center"><math>p = \frac{e^{r \Delta t} - d}{u - d}, \quad Q = e^{\sigma \Delta t}.</math></div>
In diesem Falle gilt <math>d > 0</math> und <math>p > 0</math> für alle <math>r, \sigma, \Delta t</math>.

4. In dem Falle von Dividendenzahlungen auf den Basiswert (etwa zur Zeit <math>t < T</math>), fällt der Kurs sprunghaft um den Ausschüttungsbetrag. Dies kann modelliert werden, indem die Werte von <math>S</math> im Binomialbaum zur Zeit <math>t_k</math> entsprechend vermindert werden.

5. Die Binomialmethode kann zur Trinomialmethode erweitert werden, indem zu jedem Zeitpunkt <math>t_k</math> drei Änderungsmöglichkeiten für den Kurs des Basiswertes mit Wahrscheinlichkeit <math>p, q, r</math> mit <math>p + q + r = 1</math> zugelassen werden. Details findet man in [13].

 ''% Binomialmethode fuer einen europäischen Put''
 ''% Input''
 K = 6; S0 = 5; r = 0.04; sigma = 0.3; T = 1; N = 99;
 dt = T/(N+1);
 beta = 0.5*('''exp'''(-r*dt) + '''exp'''((r + sigma^2)*dt));
 u = beta + '''sqrt'''(beta^2 - 1);
 d = 1/u;
 p = ('''exp'''(r*dt) - d)/(u - d);
 ''% 1. Schritt''
 '''for''' j=1:N + 1
   S(j,N+1) = S0*u^(j - 1)*d^(N - j+1);
 '''end'''
 ''% 2. Schritt''
 '''for''' j=1:N + 1
   V(j,N+1) = '''max'''([K - S(j,N+1) 0]);
 '''end'''
 ''% 3. Schritt''
 e = '''exp'''(- r*dt);
 '''for''' i=N - 1:1
   '''for''' j=1:i
     V(j,i) = e*(p*V(j + 1, i + 1) + (1 - p)*V(j, i + 1));
   '''end'''
 '''end'''
 ''% output''
 '''fprintf'''(’V(%f,0) = %f\n\n’, S0, V(1,1))

== 2.5 Ein Übergang vom Binomialbaum zur Black-Scholes-Formel ==
In Abschnitt 2.1 wurde der Preis einer europäischen Call-Option zur Zeit <math>t = 0</math> aus einem <math>n</math>-Perioden-Binomialbaum hergeleitet:
<div align="center">(2.24) <math>C_0 = SP(X_{p'} \ge m) - Ke^{- rT} P(X_p \ge m),</math></div>
wobei <math>X_{p'}</math> und <math>X_p</math> binomialverteilte Zufallsgrößen mit <math>B(n, p')</math> und <math>B(n, p)</math> sind. Weiter ist <math>m = \min_{0 \le k \le n} \{u^k d^{n-k} S - K \ge 0\}</math> und
<div align="center"><math>p = \frac{e^{r \Delta t} - d}{u - d}, \quad p' = pue^{-r \Delta t}.</math></div>
Wir wollen nun den Grenzwert von <math>C_0</math> für <math>n \to \infty</math> untersuchen. Wir setzen dafür <math>T = n \Delta t, u > 1</math> und <math>ud = 1</math>. Weiter sei eine Zahl \sigma > 0 definiert durch
<div align="center"><math>u = exp(\sigma \sqrt{\Delta t}).</math></div>
Der Grenzübergang beruht auf dem zentralen Grenzwertsatz, den wir in folgender Variante benutzen:

=== Satz 2.2 ===
:''Sei <math>X_n, n \in \mathbb N</math> eine Folge <math>B(n, p)</math>-verteilter Zufallsvariable auf einem Wahrscheinlichkeitsraum. Dann gilt
<div align="center"><math>\lim_{n \to \infty} P \left( \frac{X_n - np}{\sqrt{np(1 - p)}} \le x \right) = \Phi(x) = \frac{1}{\sqrt{2 \pi}} \int\limits^x_{- \infty} e^{\frac{- s^2}{2}}\, ds.</math></div>

Wir beweisen, dass im zeit-kontinuierlichen Grenzfall <math>\Delta t \to 0</math> der diskrete Optionspreis (2.24) gegen die zeit-kontinuierliche Black-Scholes-Formel konvergiert.

=== Satz 2.3 ===
:''Es gelte <math>T = n \Delta t, u = \exp(\sigma \sqrt{\Delta t}), d = \exp(- \sigma \sqrt{\Delta t})</math>. Dann folgt
<div align="center"><math>\lim_{\Delta t \to 0} C_0 = S \Phi(d_1) - Ke^{- rT} \Phi(d_2),</math></div>
:''wobei der Grenzwert <math>\Delta t \to 0</math> bei <math>T = n\Delta t</math> gerade <math>n \to \infty</math> entspricht und
<div align="center"><math>d_{1/2} = \ln(S/K) + \frac{(r \pm \sigma^2/2)T}{\sigma \sqrt{T}}.</math></div>

==== Beweis: ====
Es genügt zu zeigen, dass
<div align="center"><math>\lim_{\Delta t \to 0} P(X_{p'} m) = \Phi(d_1), \quad \lim_{\Delta t \to 0} P(X_{p} m) = \Phi(d_2)</math></div>
Wir beweisen exemplarisch den zweiten Grenzwert. Wir zeigen zuerst, dass
<div align="center">(2.25) <math>\lim_{\Delta t \to 0} np(1 - p) \left( \ln \frac{u}{d} \right)^2 = \sigma^2 T,</math></div>
<div align="center">(2.26) <math>\lim_{\Delta t \to 0} n \left( p \ln \frac{u}{d} + \ln d \right) = \left( r - \frac{\sigma^2}{2} \right) T.</math></div>
Sowohl <math>u</math> als auch <math>d</math> hängen von <math>\Delta t</math> ab; somit ist auch <math>p</math> eine Funktion von <math>\Delta t</math>. Durch Anwendung der Regel von l’Hospital auf <math>p = (\exp(r\Delta t) - d)/(u - d), u = 1/d = \exp(\sigma \sqrt{\Delta t})</math> folgt
<div align="center"><math>\lim_{\Delta t \to 0} p = \lim_{\Delta t \to 0} \frac{e^{r \Delta t} - e^{- \sigma \sqrt{\Delta t}}}{e^{\sigma \sqrt{\Delta t}} - e^{- \sigma \sqrt{\Delta t}}} = \lim_{\Delta t \to 0} \frac{r e^{r \Delta t} + \frac{\sigma}{2\sqrt{\Delta t}} e^{- \sigma \sqrt{\Delta t}}}{\frac{\sigma}{2\sqrt{\Delta t}} e^{\sigma \sqrt{\Delta t}} + \frac{\sigma}{2\sqrt{\Delta t}} e^{- \sigma \sqrt{\Delta t}}} = \frac{1}{2}.</math></div>
Eine weitere Rechnung (- die etwas mühsam ist -) zeigt, dass
<div align="center"><math>\lim_{\Delta t \to 0} \frac{2p - 1}{\sqrt{\Delta t}} = \frac{r}{\sigma} - \frac{\sigma}{2}.</math></div>
Wir schließen aus diesen beiden Beziehungen (2.25), denn
<div align="center"><math>\lim_{\Delta t \to 0} np(1 - p) \left( \ln \frac{u}{d} \right)^2 = \lim_{\Delta t \to 0} \frac{T}{\Delta t} p(1 - p) \left( 2 \sigma \sqrt{\Delta t} \right)^2 = \lim_{\Delta t \to 0} 4p(1 - p) \sigma^2 T = \sigma^2 T,</math></div>
<div align="center"><math>\lim_{\Delta t \to 0} n \left( p \ln \frac{u}{d} + \ln d \right) = \lim_{\Delta t \to 0} \frac{T}{\sqrt{\Delta t}} (2p - 1) \sigma = \left( r - \frac{\sigma^2}{2} \right) T.</math></div>
Um den zentralen Grenzwertsatz anwenden zu können, formulieren wir die Wahrscheinlichkeit für <math>X_p</math> um:
<div align="center"><math>P(X_p \ge m) = 1 - P(X_p < m) = 1 - P \left( \frac{X_p - np}{\sqrt{np(1 - p)}} < \frac{m - np}{\sqrt{np(1 - p)}} \right)</math></div>
Nach Definition von <math>m</math> gilt
<div align="center"><math>m \ln u + (n - m) \ln d - \ln \frac{K}{S} \Longleftrightarrow m \ge - \frac{\ln \frac{S}{K} + n \ln d}{\ln \frac{u}{d}}</math></div>
und es exisitiert eine Zahl <math>0 \le \alpha < 1</math>, so dass
<div align="center"><math>m = - \frac{\ln \frac{S}{K} + n \ln d}{\ln \frac{u}{d}} + \alpha.</math></div>
Damit ergibt sich
<div align="center"><math>P(Xp \ge m) = 1 - P \left( \frac{X_p - np}{\sqrt{np(1 - p)}} < \frac{- \ln \frac{S}{K} - n \left( p \ln \frac{u}{d} + \ln d \right) + \alpha \ln \frac{u}{d}}{\ln \frac{u}{d} \sqrt{np(1 - p)}} \right).</math></div>
Wegen (2.25) und <math>\alpha \ln(u/d) = 2\alpha \sigma \sqrt{\Delta t}</math> erhalten wir den Grenzwert für <math>\Delta t \to 0</math>
<div align="center"><math>P(X_p \ge m) \stackrel{(\Delta t \to 0)}{\longrightarrow} 1 - \Phi \left( \frac{- \ln(S/K) - (r - \sigma^2/2)T}{\sigma \sqrt{T}} \right),</math></div>
woraus wegen <math>1 - \Phi(-x) = \Phi(x)</math>
<div align="center"><math>P(X_p \ge m) \stackrel{(\Delta t \to 0)}{\longrightarrow} \Phi \left( \frac{\ln(S/K) + (r - \sigma^2/2)T}{\sigma \sqrt{T}} \right) = \Phi(d_2)</math></div>
und damit die Behauptung folgt. 

q.e.d.