Maxima

Binomialverteilung

Wir gehen von einem Zufallsexperiment (Basisexperiment oder auch Bernoulliexperiment) aus. Ein derartiges Experiment hat 2 mögliche Ausgänge:

  • A ... das Ereignis tritt ein
  • ¬A ... das Ereignis tritt nicht ein
  • P(A) = p ... Wahrscheinlichkeit für den Eintritt von A
  • P(¬A) = 1 - p ... Wahrscheinlichkeit, dass A nicht eintritt
Beispiel: A ... man würfelt einen Sechser, P(A) = 1/6, P(¬A) = 1 - 1/6 = 5/6

Dieses Experiment wird nun n-mal wiederholt - wir betrachten, wie oft sich das Ereignis A wiederholt und die damit zusammenhängende Wahrscheinlichkeit.
Beispiel: Wir würfeln 20 mal, wie oft kommt ein Sechser - wie hoch ist die Wahrscheinlichkeit, dass bei 20 Würfen 3 Sechser dabei sind etc.

Welche Funktionen bietet Maxima für derartige Berechnungen - welche Funktionen kann man sich leicht zurechtzimmern:

  • binomial(n,k) ... eingebaute Maxima-Funktion für n!/(k!.(n-k)!)
  • binv(n,p,k) ... selbst definierte Funktion für die Anzahl P(k), mit der eine bestimmte Eigenschaft mit Wahrscheinlichkeit p k-mal auftritt. Sie wird über die Formel für "Ziehen mit Zurücklegen" bestimmt:
    Binomialverteilung 01
  • binvert(n,p,a,b) ... meist wollen wir Bereiche wissen - "mindestens k-mal" oder "k zwischen a und b" - dies lässt sich (selbst definiert) wie folgt ausdrücken:
    Binomialverteilung 02

Binomialverteilung 03

Beispiel

Bei einem Spielautomaten gewinnt man mit der Wahrscheinlichkeit 0.4. Mit welcher Wahrscheinlichkeit gewinnt man bei 100 Spielen:

* mindestens 50
* höchstens 30
* mindestens 30 - aber höchstens 40 Spiele

Binomialverteilung 04

Wir wollen die Binomialverteilung durch eine geeignete Normalverteilung approximieren und anschließend beide Verteilungen grafisch darstellen.

Binomialverteilung 05

Darstellung:

des dritten Falls - mindestens 30, höchstens 40 Spiele

Binomialverteilung 06

Tip: Maxima verlangt bei einigen Befehlen wie bars, rectangle, ... keine klassische Liste der Form [Element1, Element2, ...] (wobei Element1 usw. wieder spezielle Listen sind), sondern die Form rectangle(Element1, Element2, ...) - rectangle([Element1, ...]) funktioniert nicht :-(
Hier muss man zu einem Trick greifen - mit Hilfe von

apply(draw2d,
  append([... Optionen ...],
         makeliste(rectangle(...))
  )
)	    
kann man Maxima "überlisten" und dem Befehle rectangle auch direkt eine Liste übergeben.

Binomialverteilung 07

MAXIMA-Datei


 
zurueck
Links:

© PH-NOe, letzte Änderung am 14. Juni 2008, erstellt von Walter Wegscheider