Intégrale d'Euler-Poisson. Détails sur les méthodes de calcul



Dans l'article, en détail, jusqu'au moindre détail, trois méthodes de prise de l'intégrale d'Euler-Poisson sont considérées. Dans l'une des méthodes, une formule de réduction auxiliaire est dérivée. Pour trouver des intégrales complexes, on peut utiliser des formules de réduction qui permettent d'abaisser le degré de l'intégrande et de calculer les intégrales correspondantes en un nombre fini d'étapes.

Cette intégrale est tirée de la fonction gaussienne: I= int limits 0inftyex2dx
Il existe une méthode mathématique très intéressante. Pour trouver l'intégrale d'origine, recherchez d'abord le carré de cette intégrale, puis prenez la racine du résultat. Pourquoi? Oui, car il est tellement plus facile et indolore d'aller aux coordonnées polaires. Par conséquent, considérons le carré de l'intégrale gaussienne:
I2= int limits 0inftyex2dx int limits 0inftyey2dy= int limits 0infty int limits 0inftye left(x2+y2 right)dxdy
Nous voyons que nous obtenons une double intégrale d'une fonction g left(x,y right)= exp left[ left(x2+y2 right) right] . À la fin de cette intégrale de surface se trouve l'élément de surface dans le système de coordonnées cartésiennes dS=dxdy .
Passons maintenant au système de coordonnées polaires:

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


Il convient de noter ici que r peut varier de 0 à + ∞, car x variait dans la même plage. Mais l'angle φ varie de 0 à π / 2, qui décrivent la région d'intégration dans le premier quart du système de coordonnées cartésiennes. En substituant à la source, on obtient:

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


En raison de la symétrie de l'intégrale et de la plage positive de valeurs de l'intégrande, nous pouvons conclure que

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



Trouvons-nous d'autres solutions? C'est intéressant! :)

Considérez la fonction g left(t right)= left(1+t right)et
Rappelons maintenant les mathématiques de l'école et réalisons une étude simple d'une fonction en utilisant des dérivées et des limites. Ce n'est pas que nous considérerons des limites complexes ici (après tout, ils ne les dépassent pas à l'école), nous discutons simplement de ce qui arrivera à la fonction si son argument tend à zéro ou à l'infini, ainsi nous allons estimer le comportement asymptotique, qui est toujours très important en mathématiques. C'est comme une évaluation qualitative de ce qui se passe.

 beginarraylg left(t right)= left(1+t right)etg left(t right)=et left(1+t right)et=tetg left(t right)=0 tot=0 gauche[ beginarraylt<0 totet>0 tog left(t right) rmaugmentet>0 totet<0 àg gauche(t droite) rmdiminue endarray droite.g left(0 right)= left(1+0 right)e0=1g left(1 right)= left(11 droite)e gauche(1 droite)=0g gauche( infty droite)= gauche(1+ infty droite)e infty=0 endarray


Il est délimité ci-dessus par l'unité sur l'intervalle (-∞; + ∞) et zéro sur l'intervalle [-1; + ∞).

Nous effectuons le changement de variables suivant t= pmx2
Et nous obtenons:

t = \ pm x ^ 2 \ to \ 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} \ droite.



Dans la première inégalité, nous limitons la variation (0,1) et dans la seconde, l'intervalle (0; + ∞), nous élevons les deux inégalités à la puissance n, car les inégalités avec des termes positifs peuvent être augmentées à n'importe quel degré positif. Nous obtenons:

\ 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} \ droite) ^ n}} \\ x> 1 \\ \ end {array} \ droite.} \\ \ end {array}


Construisons des graphes pour n = 1 pour démontrer les inégalités

Nous essayons maintenant d'intégrer les inégalités dans les limites indiquées dans les systèmes correspondants. Et combinez immédiatement tout en une seule inégalité:

 int limits10 gauche(1x2 droite)ndx< int limits10enx2dx< int limits 0inftyenx2dx< int limits 0infty frac1 left(1+x2 right)ndx


Encore une fois, si vous regardez les graphiques, cette inégalité est vraie.

Étant donné un petit remplacement, il est facile de voir que:

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


C'est-à-dire dans cette grande inégalité au milieu, nous avons l'intégrale d'Euler-Poisson, et maintenant nous devons trouver les intégrales qui se trouvent aux frontières de cette inégalité.

Trouvez l'intégrale à partir de la bordure gauche:

\ begin {array} {l} \ int \ limits_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 \ to t = \ arcsin 1 = \ frac {\ pi} {2} \\ x = 0 \ to t = \ arcsin 0 = 0 \\ \ end {array} \\ \ end {array}} \ right] = \\ = \ int \ limits_0 ^ {\ frac {\ pi} {2}} {\ cos ^ {2n} t \ cdot \ cos tdt} = \ int \ limits_0 ^ {{ \ frac {\ pi} {2}} {\ cos ^ {2n + 1} tdt} \\ \ end {array}


Pour le calculer et l'évaluer, trouvons d'abord une intégrale générale. Je vais maintenant vous montrer comment dériver la formule de réduction (en mathématiques, par de telles formules, elles signifient abaisser le degré) pour une intégrale donnée.

\ begin {array} {l} \ int \ limits_ \ alpha ^ \ beta {\ cos ^ n tdt} = \ int \ limits_ \ alpha ^ \ beta {\ cos ^ {n - 1} t \ cos tdt} = \ int \ limits_ \ 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} \ droite | _ \ alpha ^ \ beta + \ gauche ({n - 1} \ droite) \ int \ limits_ \ alpha ^ \ beta {\ cos ^ {n - 2} t \ sin ^ 2 tdt} = \\ = \ gauche. {\ cos ^ {n - 1} t \ sin t} \ droite | _ \ alpha ^ \ beta + \ gauche ({n - 1} \ droite) \ int \ limits_ \ alpha ^ \ beta {\ cos ^ {n - 2} t \ gauche ({1 - \ cos ^ 2 t} \ droite) dt} = \\ = \ gauche. {\ cos ^ {n - 1} t \ sin t} \ droite | _ \ alpha ^ \ beta + \ gauche ({n - 1} \ droite) \ int \ limits_ \ alpha ^ \ beta {\ cos ^ {n - 2} tdt} - \ left ({n - 1} \ right) \ int \ limits_ \ alpha ^ \ beta {\ cos ^ n tdt} \\ \ end {array}


 beginarrayl int limits alpha beta cosntdt= left. cosn1t sint droite| alpha beta+ gauche(n1 droite) int limits alpha beta cosn2tdt left(n1 right) int limits alpha beta cosntdt int limits alpha beta cosntdt+ left(n1 right) int limits alpha beta cosntdt= left. cosn1t sint droite| alpha beta+ gauche(n1 droite) int limits alpha beta cosn2tdtn int limits alpha beta cosntdt= left. cosn1t sint droite| alpha beta+ gauche(n1 droite) int limits alpha beta cosn2tdt int limits alpha beta cosntdt= frac1n left. cosn1t sint right| alpha beta+ fracn1n int limits alpha beta cosn2tdt endarray


Maintenant, si nous utilisons la formule de réduction, nous considérons la même intégrale, mais avec nos limites de 0 à π / 2, alors nous pouvons faire quelques simplifications:

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


Comme nous le voyons, vous pouvez l'abaisser à l'infini (dépend de n). Cependant, il y a une subtilité. La formule change selon que n est un nombre pair ou non.
Pour cela, nous considérons deux cas.

 beginarrayln=10: int limits frac pi20 cos10tdt= frac9 cdot7 cdot5 cdot310 cdot8 cdot6 cdot4 int limits frac pi20 cos2tdt= frac9 cdot7 cdot5 cdot310 cdot8 cdot6 cdot4 int limits frac pi20 left( frac12+ frac12 cos2t right)dt== frac9 cdot7 cdot5 cdot310 cdot8 cdot6 cdot4 gauche. 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 àdroite)!!n!! cdot frac pi2 endarray

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


Où est n !! - double factorielle. La double factorielle de n est notée n !! et est défini comme le produit de tous les nombres naturels dans l'intervalle [1, n] ayant la même parité que n

Du fait que 2n + 1 est un nombre impair pour toute valeur de n, nous obtenons pour la frontière gauche de notre inégalité:

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


Trouvez l'intégrale à partir de la bordure droite:
(ici, nous utilisons la même formule de réduction qui a été prouvée précédemment)

 beginarrayl int limits 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 limits frac pi20 cos2n2tdt= left[ left(2n2 right) rmevenright]= frac left(2n3 droite)!! gauche(2n2 droite)!! cdot frac pi2 endarray


Après avoir estimé les côtés gauche et droit de l'inégalité, nous effectuons quelques transformations pour évaluer les limites des côtés gauche et droit de l'inégalité, à condition que n tend vers ∞:

\ begin {array} {l} \ frac {{\ left ({2n} \ right) !!}} {{\ left ({2n + 1} \ right) !!}} <\ frac {1} { {\ sqrt n}} \ cdot I <\ frac {{\ left ({2n - 3} \ right) !!}} {{\ left ({2n - 2} \ right) !!}} \ cdot \ frac {\ pi} {2} \\ \ sqrt n \ cdot \ frac {{\ left ({2n} \ right) !!}} {{\ left ({2n + 1} \ right) !!}} <I <\ sqrt n \ cdot \ frac {{\ gauche ({2n - 3} \ droite) !!}} {{\ gauche ({2n - 2} \ droite) !!}}} cdot \ frac {\ pi} {2} \\ \ end {array}


Carré des deux côtés de l'inégalité:

n cdot frac gauche( gauche(2n droite)!! droite)2 gauche( gauche(2n+1 droite)!! droite)2<I2<n cdot frac gauche( gauche(2n3 droite)!! droite)2 gauche( gauche(2n2 droite)!! droite)2 cdot frac pi24


Faisons maintenant une petite digression. En 1655, John Wallis (un mathématicien anglais, l'un des précurseurs de l'analyse mathématique) a proposé une formule pour déterminer le nombre π. J. Wallis s'approcha d'elle, calculant l'aire d'un cercle. Ce produit converge extrêmement lentement; par conséquent, la formule de Wallis est peu utile pour le calcul pratique du nombre π. Mais c'est génial pour évaluer notre expression :)

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


Maintenant, nous transformons nos inégalités afin de voir où substituer la formule de Wallis:

 beginarrayl fracn2 left(2n+1 right)2 cdot frac1n cdot frac left( left(2n right)!! right)2 left( left(2n1 right)!! right)2<I2< frac1 frac1n cdot frac left( left(2n2 right)!! right)2 left( left(2n3 right)!! right)2 cdot frac pi24 mathop lim limitsn to infty left[ fracn2 left(2n+1 right)2 right] cdot mathop lim limitsn to infty left[ frac1n cdot frac left( left(2n right)!! right)2 gauche( gauche(2n1 droite)!! droite)2 droite]<I2< frac1 mathop lim limitsn to infty left[ frac1n cdot frac left( left(2n2 right)!! droite)2 gauche( gauche(2n3 droite)!! droite)2 droite] cdot frac pi24 frac14 cdot pi<I2< frac1 pi cdot frac pi24 to frac pi4<I2< frac pi4I2= frac pi4 toI= frac sqrt pi2 endarray


Il résulte de la formule de Wallis que les expressions gauche et droite ont tendance à π / 4 comme n → ∞
Étant donné que la fonction exp [-x²] est paire, nous supposons en toute sécurité que

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


Pour la première fois, l'intégrale gaussienne unidimensionnelle a été calculée en 1729 par Euler, puis Poisson a trouvé un moyen simple de la calculer. À cet égard, il a été appelé l'intégrale d'Euler-Poisson.

Essayons de calculer l'intégrale gaussienne. Il peut être écrit sous différentes formes. Après tout, rien ne change le nom de la variable par laquelle l'intégration a lieu.

 beginarraylI= int limits 0inftyex2dxI= int limits  inftyinftyex2dx= int limits  inftyinftyey2dy= int limits  inftyinftyez2dz endarray


Vous pouvez passer des coordonnées cartésiennes tridimensionnelles aux coordonnées sphériques et considérer le cube de l'intégrale de Gauss.

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


Le jacobien de cette transformation peut être calculé comme suit:

\ begin {array} {l} J = \ left | {\ begin {array} {* {20} c} {\ frac {{\ partial x}} {{\ partial r}}} & {\ frac {{\ partial x}} {{\ partial \ theta}} } & {\ frac {{\ partial x}} {{\ partial \ varphi}}} \\ {\ frac {{\ partial y}} {{\ partial r}}} & {\ frac {{\ partial y }} {{\ partial \ theta}}} & {\ frac {{\ partial y}} {{\ partial \ varphi}}} \\ {\ frac {{\ partial z}} {{\ partial r}} } & {\ frac {{\ partial z}} {{\ partial \ theta}}} & {\ frac {{\ partial z}} {{\ partial \ varphi}}} \\ \ end {array}} \ à droite | = \ gauche | {\ 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 \ thêta} & 0 \\ \ end {array}} \ right | = \\ = r ^ 2 \ sin \ theta \\ \ end {array}


 beginarraylI3= int limits  inftyinfty int limits  inftyinfty int limits  infty inftyex2y2z2dxdydz= int limits2 pi0 int limits 0pi int limits 0inftyer2Jdrd thetad varphi== int limits2 pi0d varphi int limits 0pi sin thetad thêta int limits 0inftyer2r2dr endarray


Nous calculons les intégrales séquentiellement, en commençant par l'intérieure.

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


Ensuite, nous obtenons:

 beginarraylI3=2 pi cdot2 cdot fracI4 toI3= piI toI2= pi toI= sqrt piI= int limits  inftyinftyex2dx= sqrt pi endarray


L'intégrale d'Euler-Poisson est souvent utilisée dans la théorie des probabilités.

J'espère que pour quelqu'un l'article sera utile et aidera à comprendre certaines techniques mathématiques :)

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


All Articles