Wintersemester 2023/24
© 2023 Prof. Dr. Rüdiger W. Braun
Beispiel: Beim Experiment ``vierfacher Wurf einer fairen Münze'' gibt es 16 Elementarereignisse ('a': Adler, 'z': Zahl) \begin{array}{cccc}
(a,a,a,a) & (a,a,a,z) & (a,a,z,a) & (a,a,z,z) \\
(a,z,a,a) & (a,z,a,z) & (a,z,z,a) & (a,z,z,z) \\
(z,a,a,a) & (z,a,a,z) & (z,a,z,a) & (z,a,z,z) \\
(z,z,a,a) & (z,z,a,z) & (z,z,z,a) & (z,z,z,z)
\end{array}
Die Aussage "Die Münze ist fair" bedeutet, dass jedes dieser Elementarereignisse die Wahrscheinlichkeit $1/16$ besitzt
\begin{equation*}
A = \{ (a,a,z,z), (a,z,a,z), (a,z,z,a),
(z,a,a,z), (z,a,z,a), (z,z,a,a) \}
\end{equation*}\Omega &=
\begin{aligned}[t]
\{ &(1,1), (1,2), (1,3), (1,4), (1,5), (1,6), \\
&(2,1), (2,2), (2,3), (2,4), (2,5), (2,6), \\
&(3,1), (3,2), (3,3), (3,4), (3,5), (3,6), \\
&(4,1), (4,2), (4,3), (4,4), (4,5), (4,6), \\
&(5,1), (5,2), (5,3), (5,4), (5,5), (5,6), \\
&(6,1), (6,2), (6,3), (6,4), (6,5), (6,6) \}
\end{aligned}
\end{align*}A = \{(1,1), (2,2), (3,3), (4,4), (5,5), (6,6) \}
\end{equation*}B = \{(3,4), (4,3) \}
\end{equation*}Aus einfachen Ereignissen werden komplexere aufgebaut.
Der Durchschnitt $A \cap B$ zweier Ereignisse $A$ und $B$ besteht aus allen Elementarereignissen, die sowohl zu $A$ als auch zu $B$ gehören.
Die Vereinigung $A \cup B$ zweier Ereignisse $A$ und $B$ besteht aus allen Elementarereignissen, die zu mindestens einem der beiden Ereignisse gehören.
Die Differenz $A \setminus B$ zweier Ereignisse $A$ und $B$ besteht aus allen Elementarereignissen, die zu $A$, aber nicht zu $B$ gehören.
Das Komplement $A^c$ eines Ereignisses $A$ besteht aus allen Elementarereignissen, die nicht zu $A$ gehören.
Die Menge aller Elementarereignisse ist der Ereignisraum. Seine Teilmengen heißen (Zufalls)-Ereignisse. Die Mengenlehre dient uns als Sprechweise, Ereignisse kurz und zweifelsfrei zu beschreiben.
verbal | mathematisch |
---|---|
Ereignisse $A$ und $B$ treffen ein | $ A \cap B $ |
Ereignis $A$ oder Ereignis $B$ trifft ein | $ A \cup B $ |
Ereignis $A$ trifft nicht ein | $ A^c $ |
Ereignis $A$ trifft ein, Ereignis $B$ aber nicht | $ A \setminus B $ |
unmögliches Ereignis | $ \emptyset $ |
sicheres Ereignis (= Ereignisraum) | $ \Omega $ |
Elementarereignis $\omega$ gehört zu $A$ | $ \omega \in A $ |
Elementarereignis $\omega$ gehört nicht zu $A$ | $ \omega \notin A $ |
alle Elementarereignisse von $A$ gehören zu $B$ | $ A \subseteq B $ |
$\Omega = \text{"Wurf eines Würfels"} = \{1,2,3,4,5,6\}$, $A = \text{"ungerade Zahl gewürfelt"} = \{1, 3, 5\}$ und $B = \text{"Zahl kleiner $4$ gewürfelt"} = \{1, 2, 3\}$
A \subseteq \Omega_1 \qquad B \subseteq \Omega_2
\end{equation*}A \times B
= \{ (a,b) \mid a \in A,\, b \in B \}
\end{equation*}$\color{blue} A = \{a_1, a_2, a_3\}$
$\color{green} B = \{b_1, b_2, b_3, b_4\} $
\begin{array}{c|cccc} & \color{green} b_1 & \color{green} b_2 & \color{green} b_3 & \color{green} b_4 \\\hline \color{blue} a_1 & (\color{blue} a_1, \color{green} b_1) & (\color{blue} a_1, \color{green} b_2) & (\color{blue} a_1, \color{green} b_3) & (\color{blue} a_1, \color{green} b_4) \\ \color{blue} a_2 & (\color{blue} a_2, \color{green} b_1) & (\color{blue} a_2, \color{green} b_2) & (\color{blue} a_2, \color{green} b_3) & (\color{blue} a_2, \color{green} b_4) \\ \color{blue} a_3 & (\color{blue} a_3, \color{green} b_1) & (\color{blue} a_3, \color{green} b_2) & (\color{blue} a_3, \color{green} b_3) & (\color{blue} a_3, \color{green} b_4) \end{array}Zweifacher Wurf eines Würfels \begin{align*} \Omega &= \{ 1,2,3,4,5,6\}^2 \\ &= \begin{aligned}[t] \{ &(1,1), (1,2), (1,3), (1,4), (1,5), (1,6), \\ &(2,1), (2,2), (2,3), (2,4), (2,5), (2,6), \\ &(3,1), (3,2), (3,3), (3,4), (3,5), (3,6), \\ &(4,1), (4,2), (4,3), (4,4), (4,5), (4,6), \\ &(5,1), (5,2), (5,3), (5,4), (5,5), (5,6), \\ &(6,1), (6,2), (6,3), (6,4), (6,5), (6,6) \} \end{aligned} \end{align*}
Für jedes Ereignis $A$ gebe es eine Zahl $P(A)$ mit
Dann ist $P$ eine Wahrscheinlichkeitsverteilung auf $\Omega$, und $(\Omega, P)$ ist ein wahrscheinlichkeitstheoretisches Modell des Zufallsexperiments
import scipy
scipy.special.factorial(70)
1.197857166996989e+100
scipy
: Bibliothek für wissenschaftliches Rechnenscipy.special
mit speziellen Funktionenscipy.stats
mit statistischen FunktionenEs gibt auch Bibliotheken, die $49!$ genau ausrechnen können
$$ 49! = 608281864034267560872252163321295376887552831379210240000000000 $$49 \cdot 48 \cdot 47 \cdot 46 \cdot 45 \cdot 44
= 10 068 347 520
\end{equation*}nPr
\cdot 4 \cdot 3 \cdot 2 \cdot 1}
= 13 983 816
$$49 \\ 6
\end{pmatrix} $nCr
n \\k
\end{pmatrix}
= \frac{n!}{k! \cdot (n-k)!}
= \frac{n \cdot (n-1) \cdots (n-k+2) \cdot (n-k+1)} {k \cdot
(k-1) \cdots 2 \cdot 1} $$
ist die Anzahl der möglichen Auswahlen von $k$ Objekten aus $n$-Objekten.n \\ k
\end{pmatrix}
$ heißt *Binomialkoeffizient.* Man sagt "$n$ über $k$".scipy.special.binom(49, 6)
13983816.0
Wie viele Möglichkeiten gibt es, zwei Ecken eines Quadrats auszuwählen?
Die Antwort ist $$ \begin{pmatrix}
4 \\ 2
\end{pmatrix}
= \frac{4 \cdot 3}{2 \cdot 1} = 6 $$
10 \\ 3
\end{pmatrix}
= \frac{10 \cdot 9 \cdot 8}{3 \cdot 2 \cdot 1}
= 120 $$
10 \\ 7
\end{pmatrix}
= \frac{10 \cdot 9 \cdot 8 \cdot 7 \cdot 6 \cdot 5 \cdot 4}{7
\cdot 6 \cdot 5 \cdot 4 \cdot 3 \cdot 2 \cdot 1}
= 120 $$
scipy.special.binom(10, 7)
120.0
scipy.special.binom(100, 33)
2.9469242702254086e+26
Für jedes $n$ und jedes $k \le n$ gelten
$$ \begin{pmatrix} n \\ 0 \end{pmatrix} = 1 $$ $$ \begin{pmatrix} n \\ 1 \end{pmatrix} = n $$ * $$ \begin{pmatrix} n \ k
\end{pmatrix}
=
\begin{pmatrix}
n \\ n-k
\end{pmatrix} $$
Wahrscheinlichkeitstheorie | Experiment |
---|---|
Zufallsvariable $X$ | Messvorrichtung |
Ereignisraum $\Omega$ | Menge aller möglichen Versuchsabläufe |
Elementarereignis $\omega$ | beobachteter Versuchsablauf |
Wert $X(\omega)$ | beobachteter Messwert |
$X$ eine Zufallsvariable auf $\Omega$. Wir schreiben zur Abkürzung (hierbei sind $a$ und $b$ irgendwelche Zahlen): \begin{align*} \{ X=a \} &= \{ \text{alle Elementarereignisse $ \omega $, für die $ X(\omega) = a $} \} \\ \{X \le a \} &= \{ \text{alle Elementarereignisse $ \omega $, für die $ X(\omega) \le a $} \} \\ \{ a < X \le b \} & = \{ \text{alle Elementarereignisse $ \omega $, für die $ a < X(\omega) \le b $} \} \end{align*}
\begin{array}{c|cccc}
n & 0 & 1 & 2 & 3 \\\hline
\rule{0pt}{2.25ex} P_X(n) & \frac18 & \frac38 & \frac38 & \frac18
\end{array}
B_{n,p}(k) =
\begin{pmatrix}
n \\ k
\end{pmatrix}
\cdot p^k \cdot (1-p)^{n-k} $$
5 \\ 2
\end{pmatrix} \left( \frac16 \right)^2 \left( 1 - \frac16 \right)^3
= 10 \cdot \frac1{36} \cdot
\frac{125}{216} = 0.1608 $$e$: Erfolg, $m$: Misserfolg, $q = 1 - p$
\begin{array}{ccc} P(eemmm) & = \color{blue} p \cdot \color{blue} p \cdot \color{green} q \cdot \color{green} q \cdot \color{green} q & = \color{blue} p^2 \cdot \color{green} q^3 \\ P(ememm) & = \color{blue} p \cdot \color{green} q \cdot \color{blue} p \cdot \color{green} q \cdot \color{green} q & = \color{blue} p^2 \cdot \color{green} q^3 \\ P(emmem) &= \color{blue} p \cdot \color{green} q \cdot \color{green} q \cdot \color{blue} p \cdot \color{green} q &= \color{blue} p^2 \cdot \color{green} q^3 \\ P(emmme) &= \color{blue} p \cdot \color{green} q \cdot \color{green} q \cdot \color{green} q \cdot \color{blue} p &= \color{blue} p^2 \cdot \color{green} q^3 \\ P(meemm) &= \color{green} q \cdot \color{blue} p \cdot \color{blue} p \cdot \color{green} q \cdot \color{green} q &= \color{blue} p^2 \cdot \color{green} q^3 \\ P(memem) &= \color{green} q \cdot \color{blue} p \cdot \color{green} q \cdot \color{blue} p \cdot \color{green} q &= \color{blue} p^2 \cdot \color{green} q^3 \\ P(memme) &= \color{green} q \cdot \color{blue} p \cdot \color{green} q \cdot \color{green} q \cdot \color{blue} p &= \color{blue} p^2 \cdot \color{green} q^3 \\ P(mmeem) &= \color{green} q \cdot \color{green} q \cdot \color{blue} p \cdot \color{blue} p \cdot \color{green} q &= \color{blue} p^2 \cdot \color{green} q^3 \\ P(mmeme) &= \color{green} q \cdot \color{green} q \cdot \color{blue} p \cdot \color{green} q \cdot \color{blue} p &= \color{blue} p^2 \cdot \color{green} q^3 \\ P(mmmee) &= \color{green} q \cdot \color{green} q \cdot \color{green} q \cdot \color{blue} p \cdot \color{blue} p &= \color{blue} p^2 \cdot \color{green} q^3 \end{array}$P(A)$ ist dann die Summe, $P(A) = 10 \cdot p^2 \cdot q^3$.
from scipy import stats
P = stats.binom(5, 1/6)
P.pmf(2)
0.16075102880658423
pmf
probability mass distribution
P = stats.binom(50, 0.01)
P.pmf(1)
0.3055586197664328
P.pmf(2)
0.07561804226543027
P = stats.binom(50, 0.99)
P.pmf(49)
0.30555861976643284