Verteilung eines Polynoms zweiten Grades einer Gaußschen Zufallsvariablen

8

Ich möchte berechnen

P(Y=aX2+bX+c<0)

wobei . Ich kann es ganz einfach mit Monte Carlo machen. Ich wurde jedoch gebeten, das analytische PDF von und dann zu berechnenf Y ( y ) Y.XN(0,σ)fY(y)Y

I=0fY(y)dy

Ich denke, wird so sein, dass nur numerisch berechnet werden kann. Da es sich jedoch um ein univariates Integral handelt, stehen numerische Methoden zur Verfügung, um es mit sehr hoher Genauigkeit zu berechnen. Gibt es einen (relativ einfachen) Ausdruck für , damit ich eine numerische Integration durchführen kann? Oder gibt es neben Monte Carlo (was meiner Meinung nach der sinnvollste Ansatz ist) eine andere Möglichkeit, zu berechnen ? I f Y ( y ) I.fY(y)IfY(y)I

DeltaIV
quelle
1
Haben Sie haben , um die pdf von finden zuerst und dann über die negative Real Linie integrieren oder können Sie mit der Methode von mpiktas wies darauf hin, das die pdf von vermeidet finden ? Y.YY
Dilip Sarwate
@ DilipSarwate, danke für die Frage. Ich wurde ausdrücklich gebeten, 1. und 2. über . Eine Antwort, die genau das tut, wäre also großartig. Andererseits kann ich darauf hinweisen, dass die Anfrage nicht zumutbar ist und dass ich bereits zwei sehr nette Methoden (MC und @mpiktas) habe, die gut funktionieren. Somit ist die Antwort auf Ihre Frage: Ich habe nicht unbedingt müssen (ich bin nicht gefeuert zu werden, wenn ich nicht), aber ich würde sicher zu schätzen in der Lage zu tun (also noch ein weitere Diskussion mit dem Anforderer zu vermeiden) . [ - , 0 ]fY(y)[,0]
DeltaIV
3
OK hier geht. Es ist zu beachten, dass was in Form des Standard-Gaußschen CDF Verwendung ausgedrückt werden kann die in der Antwort von @ mpkitas beschriebene Methode. Nehmen der Ableitung WRT, wird dann geben Sie die pdf . Sagen Sie Ihrem Anforderer außerdem, dass Sie das PDF nicht explizit integrieren müssen, um zu finden, da dessen Wert Sie bereits ermittelt haben. Φ ( ) y f Y ( y ) I I = F Y ( 0 )FY(y)=P{Yy}=P{aX2+bX+cy0}Φ()yfY(y)II=FY(0)
Dilip Sarwate
@ DilipSarwate fantastisch! Mit anderen Worten, und in der Antwort von mpkitas werden zu Funktionen von , und dann wende ich einfach die Kettenregel zur Ableitung an. Vielen Dank! x 2 yx1x2y
DeltaIV

Antworten:

16

Man beachte, dassax2+bx+c=a(xx1)(xx2) , wobei und x 2 Wurzeln des Polynoms a x 2 + b x + c sind . Wir müssen annehmen, dass x 1 und x 2 real und ungleich sind, andernfalls ist die fragliche Wahrscheinlichkeit trivial null oder eins.x1x2ax2+bx+cx1x2

Wir haben zwei Fälle.

  1. , dann ist P ( a X 2 + b X + c < 0 ) = P ( x 1 < X < x 2 ) .a>0P(aX2+bX+c<0)=P(x1<X<x2)

  2. a<0P(aX2+bX+c<0)=P(X<x1X>x2)=1P(x1<X<x2).

X

mpiktas
quelle
Ausgezeichnet! Und natürlich haben wir analytische Ausdrücke für und x 2 (Wurzeln einer quadratischen Gleichung), also ist alles sehr einfach. Vielen Dank! PS natürlich sind x 1 und x 2 in meinem Fall immer real und verschieden. Ich habe vergessen, das anzugeben. x1x2x1x2
DeltaIV
Sehr clevere Lösung. Ich wollte gerade vorschlagen, die Verteilung von abzuleiten, war aber völlig unnötig. Gut gemacht! Y
Ramhiser
@ JohnA.Ramey Siehe die Kommentare zur Hauptfrage.
Dilip Sarwate