
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 0inftye−x2dxEs 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 0Inftye−x2dx int Limits 0Inftye−y2dy= int Limits 0Infty int limit 0inftye− left(x2+y2 right)dxdyWir 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 0inftye−r2rd varphidr= int limit frac pi20d varphi int limit 0inftye−r2rdr= int limit frac pi20d varphi int limit 0inftye−r2 frac12d left(r2 right)== frac12 int limit frac pi20d varphi left( left.−e−r2 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 − inftyinftye−x2dx=2 int limit 0inftye−x2dx=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)e−tErinnern 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)e−tg′ left(t right)=e−t− left(1+t right)e−t=−te−tg′ left(t right)=0 bist=0 left[ beginarraylt<0 to−te−t>0 tog left(t right)− rmerhöhtt>0 to−te−t<0 bisg left(t right)− rmverringert endarray right.g left(0 right)= left(1+0 right)e−0=1g left(−1 right)= left(1−1 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= pmx2Und 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(1−x2 right)ndx< int limit10e−nx2dx< int limit 0inftye−nx2dx< 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 0inftye−nx2dx= left[ beginarraylp= sqrtnxp2=nx2 fracdp sqrtn=dx endarray right]= frac1 sqrtn int limit 0inftye−p2dp= 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. cosn−1t sint rechts| alpha beta+ links(n−1 rechts) int limit alpha beta cosn−2tdt− left(n−1 right) int limit alpha beta cosntdt int border alpha beta cosntdt+ left(n−1 right) int border alpha beta cosntdt= left. cosn−1t sint rechts| alpha beta+ links(n−1 rechts) int limit alpha beta cosn−2tdtn int limit alpha beta cosntdt= left. cosn−1t sint rechts| alpha beta+ links(n−1 rechts) int limit alpha beta cosn−2tdt int limit alpha beta cosntdt= frac1n left. cosn−1t sint right| alpha beta+ fracn−1n int limit alpha beta cosn−2tdt 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. cosn−1t sint right| frac pi20+ fracn−1n int border frac pi20 cosn−2tdt= left[ frac1n left. cosn−1t sint right| frac pi20=0 right]== fracn−1n int limit frac pi20 cosn−2tdt= fracn−1n left( frac1n−2 left. cosn−3t sint right| frac pi20+ fracn−3n−2 int border frac pi20 cosn−4tdt right)== fracn−1n left( fracn−3n−2 int border frac pi20 cosn−4tdt right)= fracn−1n left( fracn−3n−2 left( fracn−5)n−4 int border frac pi20 cosn−6tdt right) right)== fracn−1n left( fracn−3n−2 left( fracn−5n−4 left( fracn−7n−6 int border frac pi20 cosn−8tdt 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(n−1 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(n−1 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 beginarray∗20cx=0 tot=0x= infty tot= frac pi2 endarraydx= fracdt cos2t frac11+x2= frac11+ tan2t= cos2t endarray right]== int limit frac pi20 cos2n−2tdt= left[ left(2n−2 right)− rmevenright]= frac left(2n−3 right)!! left(2n−2 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(2n−3 right)!! left(2n−2 right)!! cdot frac pi2 sqrtn cdot frac left(2n right)!! left(2n+1 right)!!<I< sqrtn cdot frac left(2n−3 right)!! left(2n−2 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(2n−3 right)!! right)2 left( left(2n−2 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(2n−1 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 − inftyinftye−x2dx=2 int limit 0inftye−x2dx=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 0inftye−x2dxI= int border − inftyinftye−x2dx= int Grenzen − InftyInftye−y2dy= int Grenzen − InftyInftye−z2dz 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 − Inftyinftye−x2−y2−z2dxdydz= int limit2 pi0 int limit 0pi int limit 0inftye−r2Jdrd thetad varphi== int begrenzt2 pi0d varphi int Grenzen 0pi sin thetad theta int limit 0inftye−r2r2dr endarray
Wir berechnen die Integrale nacheinander, beginnend mit dem inneren.
beginarrayl int border 0inftye−r2r2dr= left[ beginarraylu=r todu=drdv=re−r2dr bisv= intre−r2dr= frac12 inte−r2dr2=− frac12e−r2 endarray right]== left. left(− frac12re−r2 right) right| 0infty+ frac12 int limit 0inftye−r2dr= frac12 int limit 0inftye−r2dr= 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 − inftyinftye−x2dx= 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 :)