Marsangriff (Wahrscheinlichkeit,

8

Angenommen, die Erde wurde von Marsraumschiffen angegriffen und wir haben Raketen, die gegen die Raumschiffe abgefeuert werden können. Die Wahrscheinlichkeit, jedes Raumschiff von jeder Rakete zu treffen und zu zerstören, beträgt (unabhängig vom Rest).m = k n n pnm=knnp

Wie groß ist die Wahrscheinlichkeit, alle Raumschiffe zu zerstören, wenn wir alle Raketen gleichzeitig loslassen, aber jede Rakete zufällig ein Raumschiff auswählt?

will198
quelle

Antworten:

14

Ein äquivalentes Modell für diesen Prozess besteht darin, zuerst die Raumschiffe in eine Flasche zu füllen. Setzen Sie die Anzahl der zerstörten Schiffe auf Null. Zählen Sie die Raketen . Um festzustellen, auf welches Schiff die Rakete zielt , schütteln Sie die Flasche gut und ziehen Sie zufällig ein Schiff aus der Flasche. Markieren Sie es mit der Wahrscheinlichkeit als zerstört; Andernfalls ändern Sie keine der Markierungen. Wenn es ursprünglich intakt war und jetzt als zerstört markiert wurde, erhöhen Sie die Anzahl der zerstörten Schiffe. Geben Sie dieses Schiff wieder in die Flasche und wiederholen Sie den Vorgang.1 , 2 , , m i pn1,2,,mip

Diese beschreibt eine Markov - Kette auf den Zählungen das wird durchlaufen werden Iterationen. Nachdem Schiffe zerstört wurden, ist die Chance, dass ein anderes zerstört wird (wodurch ein Übergang von Zustand zu Zustand ) die Chance, ein unzerstörtes Schiff (von dem es ) auszuwählen, mal die Chance, dies zu zerstören Schiff (das ist ). Das ist,m i i i + 1 n - i p0,1,,nmiii+1nip

Pr(ii+1)=ninp.

Ansonsten bleibt die Kette im Zustand . Der Anfangszustand ist . Das Interesse konzentriert sich auf die Chance, in Zustand ist nach Iterationen.i = 0 n mii=0nm

Die Übergangsmatrix dieser Wahrscheinlichkeiten, wobei die Wahrscheinlichkeit für den Übergang von nach , lässt sich leicht diagonalisieren:P i j i jPPijij

P=(1pp00001n1npn1np00001n2npn2np000011np1np00001)=V(100000npn00000n2pn00000n(n1)pn00000nnpn)V1

wo

V=((n0)(n1)(n2)(nn1)(nn)(n10)(n11)(n12)(n1n1)0(10)(11)000(00)0000)

ist Pascals Dreieck. Das Inverse ist leicht zu finden

V1=(0000(00)000(11)(10)00(22)(21)(20)0(n1n1)(1)n1+2(n12)(1)n1+1(n12)(1)n1+0(n10)(n0)(n1)(1)n+2(n2)(1)n+1(n2)(1)n+0(n0))

Lassen Sie diese zentrale (diagonale) Matrix schreiben , so dassΛ j j = n - j pΛ

Λjj=njpn.

Die Matrix für Iterationenm

(*)Pm=(VΛV1)m=VΛmV1

und natürlich

(Λm)jj=Λjjm=(njpn)m.

Wenn wir die Multiplikation in , finden wir

(**)(Pm)0n=j=0min(m,n)(1)j(nj)(njpn)m.

Dies ist die Chance, nach dem Start in Zustand im Zustand . Es ist Null für und danach ist es mal ein Polynom vom Grad (mit Termen ungleich Null von Grad bis ), was bedeutet, dass keine weitere Vereinfachung möglich erscheint. Wenn jedoch größer ist (ungefähr bis oder so), können die Potenzen in der Summe durch Exponentiale angenähert werden.n0m=0,1,,n1pnmn0mnn/p1020

(njpn)m=(1jpn)m(emp/n)j,

was wiederum über den Binomialsatz summiert werden kann und ergibt

(Pm)0n(1emp/n)n.

(Wenn und beide groß sind, kann dies weiter als angenähert werden .)mp/nnexp(emp)

Zur Veranschaulichung zeigt diese Grafik die korrekten Werte in Blau und die Näherung in Rot für wobei und . Die Unterschiede betragen höchstens wenige Prozent.m100n=5p=1/3

Zahl

Die Annäherung kann verwendet werden, um ein zu schätzen , das wahrscheinlich alle Schiffe auslöscht. Wenn Sie möchten, dass dies mindestens eine Chance gibt, wählen Sie som1εm

  1. mp/n ist groß und

  2. mn(log(n)log(ε))/p .

Dies wird aus einem Taylorreihenausdruck erster Ordnung für die Approximation erhalten. Angenommen, wir möchten eine Chance haben, alle Schiffe im Beispiel der Abbildung auszulöschen. Dann ist und95%ε=0.05

m5(log(5)log(0.05))/(1/3)=69.

Beachten Sie, dass nicht besonders groß ist, aber es kommt dahin: Die Annäherung könnte in Ordnung sein. Tatsächlich beträgt die ungefähre Chance während die korrekte Chance beträgt . 95,07 % 95,77 %mp/n=69(1/3)/5=4.695.07%95.77%


Dies ist eine modifizierte Coupon Sammler Problem , bei dem jeder Coupon, der nur eine gefunden wird , hat Chance, nützlich. Die hier verwendete Methode erzeugt die gesamte Verteilung der zerstörten Schiffe für jedes : Überprüfen Sie einfach die erste Zeile von .m P mpmPm

whuber
quelle
2
Ich habe es auch getan: Anscheinend verallgemeinert die Verteilungsfunktion einer "Fisher Verteilung". Siehe stats.stackexchange.com/questions/203410 . Eine Kopie von Fischers Originalpapier (1929) ist als PDF unter hekyll.services.adelaide.edu.au/dspace/bitstream/2440/15201/1/… verfügbar . ξ()ξ
whuber
1
Neben dem Problem des Bildchensammler dieses Problem betrifft auch die Belegungsproblem und die Annäherung von Exponentialgrößen bezieht sich auf die Poissonization hier getan math.stackexchange.com/a/631534/466748 (Ping @ Ben da er viele dieser Probleme hat)
Sextus Empiricus
@Martijn Vielen Dank, dass Sie darauf hingewiesen haben: Dieselbe Markov-Kette erscheint im Belegungsproblem.
whuber
1
Ich stelle mir vor, Sie könnten dies auch als multinomiale Verteilung betrachten, bei der Kugeln über Zellen mit der ersten Wahrscheinlichkeit und der letzten Zellwahrscheinlichkeit . Drücken Sie das Problem dann eher wie ein kombinatorisches Problem aus und verwenden Sie Stirling-Zahlen für Gleichung ** (und möglicherweise können ihre Eigenschaften verwendet werden, um zu derselben Annäherung zu gelangen). (nicht, dass es besser oder möglicherweise sogar schlechter ist, aber es gibt der Lösung nur einen anderen Blickwinkel)n + 1 n p / n 1 - pknn+1np/n1p
Sextus Empiricus