Euler-Poisson-Integral. Details zu Berechnungsmethoden



In dem Artikel werden bis ins kleinste Detail drei Methoden zur Verwendung des Euler-Poisson-Integrals betrachtet. Bei einer der Methoden wird eine Hilfsreduktionsformel abgeleitet. Um einige komplexe Integrale zu finden, kann man Reduktionsformeln verwenden, mit denen man den Grad des Integranden senken und die entsprechenden Integrale in einer endlichen Anzahl von Schritten berechnen kann.

Dieses Integral stammt aus der Gaußschen Funktion: I= int limit 0inftyex2dx
Es gibt einen sehr interessanten mathematischen Weg. Um das ursprüngliche Integral zu finden, suchen Sie zuerst nach dem Quadrat dieses Integrals und ziehen Sie dann die Wurzel aus dem Ergebnis. Warum? Ja, weil es so viel einfacher und schmerzloser ist, zu den Polarkoordinaten zu gehen. Betrachten Sie daher das Quadrat des Gaußschen Integrals:
I2= int Limits 0Inftyex2dx int Limits 0Inftyey2dy= int Limits 0Infty int limit 0inftye left(x2+y2 right)dxdy
Wir sehen, dass wir ein Doppelintegral einer Funktion erhalten g left(x,y right)= exp left[ left(x2+y2 right) right] . Am Ende dieses Oberflächenintegrals befindet sich das Flächenelement im kartesischen Koordinatensystem dS=dxdy .
Gehen wir nun zum Polarkoordinatensystem über:

 beginarrayldS=dxdy=rd varphi cdotdr left. beginarraylx=r cos varphiy=r sin varphi endarray right| bisx2 cos2 varphi+y2 sin2 varphi=r2 bisx2+y2=r2 endarray


Hierbei ist zu beachten, dass r von 0 bis + ∞ variieren kann, weil x variierte im gleichen Bereich. Der Winkel φ variiert jedoch von 0 bis π / 2, was den Integrationsbereich im ersten Viertel des kartesischen Koordinatensystems beschreibt. Wenn wir die Quelle einsetzen, erhalten wir:

 beginarraylI2= int border 0infty int border 0inftye left(x2+y2 right)dxdy= int limit frac pi20 int limit 0inftyer2rd varphidr= int limit frac pi20d varphi int limit 0inftyer2rdr= int limit frac pi20d varphi int limit 0inftyer2 frac12d left(r2 right)== frac12 int limit frac pi20d varphi left( left.er2 right| 0infty right)= frac12 int border frac pi20d varphi left(e infty left(e0 right) right)= frac12 int border frac pi20d varphi= frac12 left( left. varphi right| frac pi20 right)= frac pi4I2= frac pi4 toI= sqrt frac pi4= frac sqrt pi2 endarray


Aufgrund der Symmetrie des Integrals und des positiven Wertebereichs des Integranden können wir daraus schließen

 int limit  inftyinftyex2dx=2 int limit 0inftyex2dx=2 cdot frac sqrt pi2= sqrt pi



Lassen Sie uns weitere Lösungen finden? Das ist interessant! :) :)

Betrachten Sie die Funktion g left(t right)= left(1+t right)et
Erinnern wir uns nun an die Schulmathematik und führen eine einfache Untersuchung einer Funktion unter Verwendung von Ableitungen und Grenzen durch. Es ist nicht so, dass wir hier komplexe Grenzen berücksichtigen (schließlich werden sie in der Schule nicht überschritten), wir diskutieren nur, was mit der Funktion passieren wird, wenn ihr Argument gegen Null oder gegen Unendlich tendiert, und schätzen daher das asymptotische Verhalten, das in der Mathematik immer sehr wichtig ist. Dies ist wie eine qualitative Einschätzung dessen, was passiert.

 beginarraylg left(t right)= left(1+t right)etg left(t right)=et left(1+t right)et=tetg left(t right)=0 bist=0 left[ beginarraylt<0 totet>0 tog left(t right) rmerhöhtt>0 totet<0 bisg left(t right) rmverringert endarray right.g left(0 right)= left(1+0 right)e0=1g left(1 right)= left(11 right)e left(1 right)=0g left( infty right)= left(1+ infty right)e infty=0 endarray


Es ist oben durch die Einheit im Intervall (-∞; + ∞) und durch Null im Intervall [-1; + ∞) begrenzt.

Wir nehmen die folgende Änderung von Variablen vor t= pmx2
Und wir bekommen:

t = \ pm x ^ 2 \ nach \ left \ {\ begin {array} {l} 0 <\ left ({1 - x ^ 2} \ right) e ^ {x ^ 2} <1 \\ 0 < \ left ({1 + x ^ 2} \ right) e ^ {- x ^ 2} <1 \\ \ end {array} \ right. \ to \ left \ {\ begin {array} {l} 0 <\ left ({1 - x ^ 2} \ right) <e ^ {- x ^ 2} \\ 0 <e ^ {- x ^ 2} <\ frac {1} {{1 + x ^ 2}} \\ \ end {array} \ right.



In der ersten Ungleichung beschränken wir die Variation (0,1) und in der zweiten, dem Intervall (0; + ∞), erhöhen wir beide Ungleichungen auf die Potenz n, da Ungleichungen mit positiven Termen in jedem positiven Ausmaß erhöht werden können. Wir bekommen:

\ begin {array} {* {20} c} {\ left \ {\ begin {array} {l} \ left ({1 - x ^ 2} \ right) ^ n <e ^ {- nx ^ 2} \\ 0 <x <1 \\ \ end {array} \ right.} & {\ Left \ {\ begin {array} {l} e ^ {- nx ^ 2} <\ frac {1} {{\ left ({1 + x ^ 2} \ right) ^ n}} \\ x> 1 \\ \ end {array} \ right.} \\ \ end {array}


Konstruieren wir Graphen für n = 1, um die Ungleichungen zu demonstrieren

Nun versuchen wir, die Ungleichungen innerhalb der in den entsprechenden Systemen angegebenen Grenzen zu integrieren. Und kombinieren Sie sofort alles zu einer Ungleichung:

 int limit10 left(1x2 right)ndx< int limit10enx2dx< int limit 0inftyenx2dx< int limit 0infty frac1 left(1+x2 right)ndx


Wenn Sie sich die Grafiken ansehen, ist diese Ungleichung wieder wahr.

Bei einem kleinen Ersatz ist leicht zu erkennen, dass:

 int limit 0inftyenx2dx= left[ beginarraylp= sqrtnxp2=nx2 fracdp sqrtn=dx endarray right]= frac1 sqrtn int limit 0inftyep2dp= frac1 sqrtnI


Das heißt, In dieser großen Ungleichung in der Mitte haben wir das Euler-Poisson-Integral, und jetzt müssen wir die Integrale finden, die an den Grenzen dieser Ungleichung stehen.

Suchen Sie das Integral am linken Rand:

\ begin {array} {l} \ int \ border_0 ^ 1 {\ left ({1 - x ^ 2} \ right) ^ n dx} = \ left [{\ begin {array} {* {20} c} \ begin {array} {l} x = \ sin t \\ dx = \ cos tdt \\ 1 - x ^ 2 = 1 - \ sin ^ 2 t = \ cos ^ 2 t \\ \ end {array} & \ begin {array} {l} x = 1 \ bis t = \ arcsin 1 = \ frac {\ pi} {2} \\ x = 0 \ bis t = \ arcsin 0 = 0 \\ \ end {array} \\ \ end {array}} \ right] = \\ = \ int \ limit_0 ^ {\ frac {\ pi} {2}} {\ cos ^ {2n} t \ cdot \ cos tdt} = \ int \ limit_0 ^ { \ frac {\ pi} {2}} {\ cos ^ {2n + 1} tdt} \\ \ end {array}


Um es zu berechnen und zu bewerten, finden wir zunächst ein allgemeines Integral. Jetzt werde ich Ihnen zeigen, wie Sie die Reduktionsformel (in der Mathematik bedeuten solche Formeln das Verringern des Grades) für ein bestimmtes Integral ableiten.

\ begin {array} {l} \ int \ limit_ \ alpha ^ \ beta {\ cos ^ n tdt} = \ int \ limit_ \ alpha ^ \ beta {\ cos ^ {n - 1} t \ cos tdt} = \ int \ limit_ \ alpha ^ \ beta {\ cos ^ {n - 1} t \ cdot d \ left ({\ sin t} \ right)} = \\ = \ left [{\ begin {array} {* { 20} c} {u = \ cos ^ {n - 1} t} & {du = - \ left ({n - 1} \ right) \ cos ^ {n - 2} t \ sin tdt} \\ {dv = d \ left ({\ sin t} \ right)} & {v = \ sin t} \\ \ end {array}} \ right] = \\ = \ left. {\ cos ^ {n - 1} t \ sin t} \ rechts | _ \ alpha ^ \ beta + \ links ({n - 1} \ rechts) \ int \ limit_ \ alpha ^ \ beta {\ cos ^ {n - 2} t \ sin ^ 2 tdt} = \\ = \ left. {\ cos ^ {n - 1} t \ sin t} \ rechts | _ \ alpha ^ \ beta + \ links ({n - 1} \ rechts) \ int \ limit_ \ alpha ^ \ beta {\ cos ^ {n - 2} t \ left ({1 - \ cos ^ 2 t} \ right) dt} = \\ = \ left. {\ cos ^ {n - 1} t \ sin t} \ rechts | _ \ alpha ^ \ beta + \ links ({n - 1} \ rechts) \ int \ limit_ \ alpha ^ \ beta {\ cos ^ {n - 2} tdt} - \ left ({n - 1} \ right) \ int \ limit_ \ alpha ^ \ beta {\ cos ^ n tdt} \\ \ end {array}


 beginarrayl int limit alpha beta cosntdt= left. cosn1t sint rechts| alpha beta+ links(n1 rechts) int limit alpha beta cosn2tdt left(n1 right) int limit alpha beta cosntdt int border alpha beta cosntdt+ left(n1 right) int border alpha beta cosntdt= left. cosn1t sint rechts| alpha beta+ links(n1 rechts) int limit alpha beta cosn2tdtn int limit alpha beta cosntdt= left. cosn1t sint rechts| alpha beta+ links(n1 rechts) int limit alpha beta cosn2tdt int limit alpha beta cosntdt= frac1n left. cosn1t sint right| alpha beta+ fracn1n int limit alpha beta cosn2tdt endarray


Wenn wir nun die Reduktionsformel verwenden, betrachten wir dasselbe Integral, aber mit unseren Grenzen von 0 bis π / 2, können wir einige Vereinfachungen vornehmen:

 beginarrayl int border frac pi20 cosntdt= frac1n left. cosn1t sint right| frac pi20+ fracn1n int border frac pi20 cosn2tdt= left[ frac1n left. cosn1t sint right| frac pi20=0 right]== fracn1n int limit frac pi20 cosn2tdt= fracn1n left( frac1n2 left. cosn3t sint right| frac pi20+ fracn3n2 int border frac pi20 cosn4tdt right)== fracn1n left( fracn3n2 int border frac pi20 cosn4tdt right)= fracn1n left( fracn3n2 left( fracn5)n4 int border frac pi20 cosn6tdt right) right)== fracn1n left( fracn3n2 left( fracn5n4 left( fracn7n6 int border frac pi20 cosn8tdt right) right) right)=... endarray


Wie wir sehen, können Sie es auf unendlich senken (hängt von n ab). Es gibt jedoch eine Subtilität. Die Formel ändert sich abhängig davon, ob n eine gerade Zahl ist oder nicht.
Hierzu betrachten wir zwei Fälle.

 beginarrayln=10: int border frac pi20 cos10tdt= frac9 cdot7 cdot5 cdot310 cdot8 cdot6 cdot4 int limit frac pi20 cos2tdt= frac9 cdot7 cdot5 cdot310 cdot8 cdot6 cdot4 int limit frac pi20 left( frac12+ frac12 cos2t right)dt== frac9 cdot7 cdot5 cdot310 cdot8 cdot6 cdot4 left. left( frac12t+ frac12 sin2t right) right| frac pi20= frac9 cdot7 cdot5 cdot310 cdot8 cdot6 cdot4 cdot frac pi4= frac9 cdot7 cdot5 cdot3 cdot110 cdot8 cdot6 cdot4 cdot2 cdot frac pi2== frac left(n1 right)!!n!! cdot frac pi2 endarray

 beginarrayln=9: int limit frac pi20 cos9tdt= frac8 cdot6 cdot4 cdot29 cdot7 cdot5 cdot3 int limit frac pi20 costdt= frac8 cdot6 cdot4 cdot29 cdot7 cdot5 cdot3übrig. left( sint right) right| frac pi20== frac8 cdot6 cdot4 cdot29 cdot7 cdot5 cdot3 cdot1= frac left(n1 right)!!n!! endarray


Wo ist n !! - doppelte Fakultät. Die doppelte Fakultät von n wird mit n bezeichnet !! und ist definiert als das Produkt aller natürlichen Zahlen im Intervall [1, n] mit der gleichen Parität wie n

Aufgrund der Tatsache, dass 2n + 1 eine ungerade Zahl für jeden Wert von n ist, erhalten wir für die linke Grenze unserer Ungleichung:

 int limit frac pi20 cos2n+1tdt= frac left(2n right)!! left(2n+1 right)!!


Suchen Sie das Integral am rechten Rand:
(hier verwenden wir die gleiche Reduktionsformel, die zuvor bewiesen wurde)

 beginarrayl int border 0infty frac1 left(1+x2 right)ndx= left[ beginarraylx= tant to beginarray20cx=0 tot=0x= infty tot= frac pi2 endarraydx= fracdt cos2t frac11+x2= frac11+ tan2t= cos2t endarray right]== int limit frac pi20 cos2n2tdt= left[ left(2n2 right) rmevenright]= frac left(2n3 right)!! left(2n2 right)!! cdot frac pi2 endarray


Nachdem wir die linke und rechte Seite der Ungleichung geschätzt haben, führen wir einige Transformationen durch, um die Grenzen der linken und rechten Seite der Ungleichung zu bewerten, vorausgesetzt, n tendiert zu ∞:

 beginarrayl frac left(2n right)!! left(2n+1 right)!!< frac1 sqrtn cdotI< frac left(2n3 right)!! left(2n2 right)!! cdot frac pi2 sqrtn cdot frac left(2n right)!! left(2n+1 right)!!<I< sqrtn cdot frac left(2n3 right)!! left(2n2 right)!! cdot frac pi2 endarray


Quadrieren Sie beide Seiten der Ungleichung:

n cdot frac left( left(2n right)!! right)2 left( left(2n+1 right)!! right)2<I2<n cdot frac left( left(2n3 right)!! right)2 left( left(2n2 right)!! right)2 cdot frac pi24


Jetzt machen wir einen kleinen Exkurs. 1655 schlug John Wallis (ein englischer Mathematiker, einer der Vorläufer der mathematischen Analyse) eine Formel zur Bestimmung der Zahl π vor. J. Wallis kam zu ihr und berechnete die Fläche eines Kreises. Dieses Produkt konvergiert extrem langsam, daher ist die Wallis-Formel für die praktische Berechnung der Zahl π von geringem Nutzen. Aber um unseren Ausdruck zu bewerten, ist es großartig :)

 pi= mathop lim limitn to infty frac1n left[ frac left(2n right)!! left(2n1 right)!! right]2


Jetzt transformieren wir unsere Ungleichung, damit wir sehen können, wo die Wallis-Formel ersetzt werden kann:

\ begin {array} {l} \ frac {{n ^ 2}} {{\ left ({2n + 1} \ right) ^ 2}} \ cdot \ frac {1} {n} \ cdot \ frac { {\ left ({\ left ({2n} \ right) !!} \ right) ^ 2}} {{\ left ({\ left ({2n - 1} \ right) !!} \ right) ^ 2} } <I ^ 2 <\ frac {1} {{\ frac {1} {n} \ cdot \ frac {{\ left ({\ left ({2n - 2} \ right) !!} \ right) ^ 2 }} {{\ left ({\ left ({2n - 3} \ right) !!} \ right) ^ 2}}} \ cdot \ frac {{\ pi ^ 2}} {4} \\ \ mathop {\ lim} \ limit_ {n \ to \ infty} \ left [{\ frac {{n ^ 2}} {{\ left ({2n + 1} \ right) ^ 2}}} \ right] \ cdot \ mathop {\ lim} \ limit_ {n \ to \ infty} \ left [{\ frac {1} {n} \ cdot \ frac {{\ left ({\ left ({2n} \ right) !!} \ right ) ^ 2}} {{\ left ({\ left ({2n - 1} \ right) !!} \ right) ^ 2}}} \ right] <I ^ 2 <\ frac {1} {{\ mathop {\ lim} \ limit_ {n \ to \ infty} \ left [{\ frac {1} {n} \ cdot \ frac {{\ left ({\ left ({2n - 2} \ right) !!} \ rechts) ^ 2}} {{\ left ({\ left ({2n - 3} \ right) !!} \ right) ^ 2}}} \ right]}} \ cdot \ frac {{\ pi ^ 2} } {4} \\ \ frac {1} {4} \ cdot \ pi <I ^ 2 <\ frac {1} {\ pi} \ cdot \ frac {{\ pi ^ 2}} {4} \ to \ frac {\ pi} {4} <I ^ 2 <\ frac {\ pi} {4} \\ I ^ 2 = \ frac {\ pi} {4} \ to I = \ frac {{\ sqrt \ pi}} {2} \\ \ end {array}


Aus der Wallis-Formel folgt, dass sowohl der linke als auch der rechte Ausdruck zu π / 4 als n → ∞ tendieren
Aufgrund der Tatsache, dass die Funktion exp [-x²] gerade ist, gehen wir davon aus

 int limit  inftyinftyex2dx=2 int limit 0inftyex2dx=2 cdot frac sqrt pi2= sqrt pi


Zum ersten Mal wurde das eindimensionale Gaußsche Integral 1729 von Euler berechnet, dann fand Poisson einen einfachen Weg, es zu berechnen. In dieser Hinsicht wurde es das Euler-Poisson-Integral genannt.

Versuchen wir, das Gaußsche Integral zu berechnen. Es kann in verschiedenen Formen geschrieben werden. Schließlich ändert nichts den Namen der Variablen, unter der die Integration stattfindet.

 beginarraylI= int border 0inftyex2dxI= int border  inftyinftyex2dx= int Grenzen  InftyInftyey2dy= int Grenzen  InftyInftyez2dz endarray


Sie können von dreidimensionalen kartesischen zu sphärischen Koordinaten wechseln und den Würfel des Gauß-Integrals betrachten.

\ left \ {\ begin {array} {l} x = r \ sin \ theta \ cos \ varphi \\ y = r \ sin \ theta \ sin \ varphi \\ z = r \ cos \ theta \\ \ end {array} \ right. \ bis x ^ 2 + y ^ 2 + z ^ 2 = r ^ 2


Der Jacobi dieser Transformation kann wie folgt berechnet werden:

\ begin {array} {l} J = \ left | {\ begin {array} {* {20} c} {\ frac {{\ partielle x}} {{\ partielle r}}} & {\ frac {{\ partielle x}} {{\ partielle \ theta}} } & {\ frac {{\ partielles x}} {{\ partielles \ varphi}} \\ {\ frac {{\ partielles y}} {{\ partielles r}}} & {\ frac {{\ partielles y }} {{\ partielle \ theta}}} & {\ frac {{\ partielle y}} {{\ partielle \ varphi}}} \\ {\ frac {{\ partielle z}} {{\ partielle r}} } & {\ frac {{\ partielles z}} {{\ partielles \ theta}} & {\ frac {{\ partielles z}} {{\ partielles \ varphi}} \\ \ end {array}} \ rechts | = \ left | {\ begin {array} {* {20} c} {\ sin \ theta \ cos \ varphi} & {r \ cos \ theta \ cos \ varphi} & {- r \ sin \ theta \ sin \ varphi} \\ {\ sin \ theta \ sin \ varphi} & {r \ cos \ theta \ sin \ varphi} & {r \ sin \ theta \ cos \ varphi} \\ {r \ cos \ theta} & {- r \ sin \ Theta} & 0 \\ \ end {array}} \ right | = \\ = r ^ 2 \ sin \ theta \\ \ end {array}


 beginarraylI3= int Limits  InftyInfty int Limits  InftyInfty Int Limits  Inftyinftyex2y2z2dxdydz= int limit2 pi0 int limit 0pi int limit 0inftyer2Jdrd thetad varphi== int begrenzt2 pi0d varphi int Grenzen 0pi sin thetad theta int limit 0inftyer2r2dr endarray


Wir berechnen die Integrale nacheinander, beginnend mit dem inneren.

 beginarrayl int border 0inftyer2r2dr= left[ beginarraylu=r todu=drdv=rer2dr bisv= intrer2dr= frac12 inter2dr2= frac12er2 endarray right]== left. left( frac12rer2 right) right| 0infty+ frac12 int limit 0inftyer2dr= frac12 int limit 0inftyer2dr= frac12 cdot fracI2= fracI4 int limit 0pi sin thetad theta= left. left( cos theta right) right| 0pi= left( cos pi right) left( cos0 right)=1+1=2 int limit2 pi0d varphi= left. varphi right|2 pi0=2 pi endarray


Dann erhalten wir als Ergebnis:

 beginarraylI3=2 pi cdot2 cdot fracI4 bisI3= piI bisI2= pi bisI= sqrt piI= int limit  inftyinftyex2dx= sqrt pi endarray


Das Euler-Poisson-Integral wird häufig in der Wahrscheinlichkeitstheorie verwendet.

Ich hoffe, dass der Artikel für jemanden nützlich ist und hilft, einige mathematische Techniken zu verstehen :)

Source: https://habr.com/ru/post/de470553/


All Articles