Integral de Euler-Poisson. Detalles sobre métodos de cálculo



En el artículo, en detalle, hasta el más mínimo detalle, se consideran tres métodos para tomar la integral de Euler-Poisson. En uno de los métodos, se deriva una fórmula de reducción auxiliar. Para encontrar algunas integrales complejas, uno puede usar fórmulas de reducción que le permitan a uno bajar el grado del integrando y calcular las integrales correspondientes en un número finito de pasos.

Esta integral se toma de la función gaussiana: I= int limits 0inftyex2dx
Hay una forma matemática muy interesante. Para encontrar la integral original, primero busque el cuadrado de esta integral y luego saque la raíz del resultado. Por qué Sí, porque es mucho más fácil e indoloro ir a las coordenadas polares. Por lo tanto, considere el cuadrado de la integral gaussiana:
I2= int limits 0inftyex2dx int limits 0inftyey2dy= int limits 0infty int limits 0inftye left(x2+y2 right)dxdy
Vemos que obtenemos una doble integral de alguna función g left(x,y right)= exp left[ left(x2+y2 right) right] . Al final de esta integral de superficie se encuentra el elemento de área en el sistema de coordenadas cartesianas dS=dxdy .
Ahora pasemos al sistema de coordenadas polares:

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


Aquí debe tenerse en cuenta que r puede variar de 0 a + ∞, porque x variado dentro del mismo rango. Pero el ángulo φ varía de 0 a π / 2, que describe la región de integración en el primer trimestre del sistema de coordenadas cartesianas. Sustituyendo en la fuente, obtenemos:

 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


Debido a la simetría de la integral y el rango positivo de valores del integrando, podemos concluir que

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



¿Encontremos algunas soluciones más? Esto es interesante! :)

Considera la función g left(t right)= left(1+t right)et
Ahora recordemos las matemáticas escolares y realicemos un estudio simple de una función usando derivadas y límites. No es que consideremos límites complejos aquí (después de todo, no los pasan en la escuela), solo discutimos lo que sucederá con la función si su argumento tiende a cero o al infinito, por lo tanto, calcularemos el comportamiento asintótico, que siempre es muy importante en matemáticas. Esto es como una evaluación cualitativa de lo que está sucediendo.

 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 left[ beginarraylt<0 totet>0 tog left(t right) rmaumentat>0 totet<0 tog left(t right) rmdisminuye 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


Está limitado anteriormente por la unidad en el intervalo (-∞; + ∞) y cero en el intervalo [-1; + ∞).

Realizamos el siguiente cambio de variables t= pmx2
Y obtenemos:

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} \ right. $



En la primera desigualdad, restringimos la variación (0,1), y en la segunda, el intervalo (0; + ∞), elevamos ambas desigualdades a la potencia n, ya que las desigualdades con términos positivos pueden elevarse a cualquier grado positivo. Obtenemos:

\ 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}


Construyamos gráficas para n = 1 para demostrar las desigualdades.

Ahora intentamos integrar las desigualdades dentro de los límites indicados en los sistemas correspondientes. E inmediatamente combine todo en una desigualdad:

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


Nuevamente, si nos fijamos en los gráficos, entonces esta desigualdad es cierta.

Dado un pequeño reemplazo, es fácil ver que:

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


Es decir en esa gran desigualdad en el medio, tenemos la integral de Euler-Poisson, y ahora necesitamos encontrar las integrales que se encuentran en los límites de esta desigualdad.

Encuentra la integral desde el borde izquierdo:

\ 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}


Para calcularlo y evaluarlo, primero encontremos una integral general. Ahora te mostraré cómo derivar la fórmula de reducción (en matemáticas, por tales fórmulas significan bajar el grado) para una integral dada.

\ 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} y {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} \ right | _ \ alpha ^ \ beta + \ left ({n - 1} \ right) \ int \ limits_ \ alpha ^ \ beta {\ cos ^ {n - 2} t \ sin ^ 2 tdt} = \\ = \ left. {\ cos ^ {n - 1} t \ sin t} \ right | _ \ alpha ^ \ beta + \ left ({n - 1} \ right) \ int \ limits_ \ alpha ^ \ beta {\ cos ^ {n - 2} t \ left ({1 - \ cos ^ 2 t} \ right) dt} = \\ = \ left. {\ cos ^ {n - 1} t \ sin t} \ right | _ \ alpha ^ \ beta + \ left ({n - 1} \ right) \ 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 right| alpha beta+ left(n1 right) 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 right| alpha beta+ left(n1 right) int limits alpha beta cosn2tdtn int limits alpha beta cosntdt= left. cosn1t sint right| alpha beta+ left(n1 right) int limits alpha beta cosn2tdt int limits alpha beta cosntdt= frac1n left. cosn1t sint right| alpha beta+ fracn1n int limits alpha beta cosn2tdt endarray


Ahora, si usamos la fórmula de reducción consideramos la misma integral, pero con nuestros límites de 0 a π / 2, entonces podemos hacer algunas simplificaciones:

 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 left. cosn3t sint right| 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


Como vemos, puedes bajarlo al infinito (depende de n). Sin embargo, hay una sutileza. La fórmula cambia dependiendo de si n es un número par o no.
Para esto, consideramos dos casos.

 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 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 limits frac pi20 cos9tdt= frac8 cdot6 cdot4 cdot29 cdot7 cdot5 cdot3 int limits frac pi20 costdt= frac8 cdot6 cdot4 cdot29 cdot7 cdot5 cdot3alaizquierda. left( sint right) right| frac pi20== frac8 cdot6 cdot4 cdot29 cdot7 cdot5 cdot3 cdot1= frac left(n1 right)!!n!! endarray$


Donde esta n !! - doble factorial. ¡El factorial doble de n se denota por n! y se define como el producto de todos los números naturales en el intervalo [1, n] que tiene la misma paridad que n

Debido al hecho de que 2n + 1 es un número impar para cualquier valor de n, obtenemos para el límite izquierdo de nuestra desigualdad:

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


Encuentre la integral desde el borde derecho:
(Aquí usamos la misma fórmula de reducción que se demostró anteriormente)

\ begin {array} {l} \ int \ limits_0 ^ \ infty {\ frac {1} {{\ left ({1 + x ^ 2} \ right) ^ n}} dx = \ left [\ begin {array } {l} x = \ tan t \ to \ begin {array} {* {20} c} {x = 0 \ to t = 0} \\ {x = \ infty \ to t = \ frac {\ pi} {2}} \\ \ end {array} \\ dx = \ frac {{dt}} {{\ cos ^ 2 t}} \\ \ frac {1} {{1 + x ^ 2}}} \ \ frac {1} {{1 + \ tan ^ 2 t}} = \ cos ^ 2 t \\ \ end {array} \ right]} = \\ = \ int \ limits_0 ^ {\ frac {\ pi} {2} } {\ cos ^ {2n - 2} tdt} = \ left [{\ left ({2n - 2} \ right) - {\ rm {even}}} right] = \ frac {{\ left ({2n - 3} \ right) !!}} {{\ left ({2n - 2} \ right) !!}} \ cdot \ frac {\ pi} {2} \\ \ end {array}


Después de haber estimado los lados izquierdo y derecho de la desigualdad, hacemos algunas transformaciones para evaluar los límites de los lados izquierdo y derecho de la desigualdad, siempre que n tienda a ∞:

 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


Cuadrar ambos lados de la desigualdad:

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


Ahora hagamos una pequeña digresión. En 1655, John Wallis (un matemático inglés, uno de los precursores del análisis matemático) propuso una fórmula para determinar el número π. J. Wallis se acercó a ella, calculando el área de un círculo. Este producto converge extremadamente lento; por lo tanto, la fórmula de Wallis es de poca utilidad para el cálculo práctico del número π. Pero es genial para evaluar nuestra expresión :)

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


Ahora transformamos nuestra desigualdad para que podamos ver dónde sustituir la fórmula 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 left( left(2n1 right)!! right)2 right]<I2< frac1 mathop lim limitsn to infty left[ frac1n cdot frac left( left(2n2 right)!! derecha)2 left( left(2n3 right)!! right)2 right] 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


De la fórmula de Wallis se deduce que las expresiones izquierda y derecha tienden a π / 4 como n → ∞
Debido al hecho de que la función exp [-x²] es par, asumimos con seguridad que

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


Por primera vez, Euler calculó la integral gaussiana unidimensional en 1729, luego Poisson encontró una manera simple de calcularla. En este sentido, se llamó la integral de Euler - Poisson.

Intentemos calcular la integral gaussiana. Se puede escribir en diferentes formas. Después de todo, nada cambia el nombre de la variable por la cual se está llevando a cabo la integración.

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


Puede pasar de coordenadas cartesianas tridimensionales a esféricas y considerar el cubo de la integral de Gauss.

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


El jacobiano de esta transformación se puede calcular de la siguiente manera:

\ 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}} \ derecho | = \ left | {\ begin {array} {* {20} c} {\ sin \ theta \ cos \ varphi} & {r \ cos \ theta \ cos \ varphi} & {- r \ sin \ theta \ sin \ varphi} \\ {\ sin \ theta \ sin \ varphi} y {r \ cos \ theta \ sin \ varphi} y {r \ sin \ theta \ cos \ varphi} \\ {r \ cos \ theta} y {- r \ sin \ theta} & 0 \\ \ end {array}} \ right | = \\ = r ^ 2 \ sin \ theta \\ \ end {array}


 beginarraylI3= int limits  inftyinfty int limits  inftyinfty int limits  inftyinftyex2y2z2dxdydz= int limits2 pi0 int limits 0pi int limits 0inftyer2Jdrd thetad varphi== int limits2 pi0d varphi int limits 0pi sin thetad theta int limits 0inftyer2r2dr endarray


Calculamos las integrales secuencialmente, comenzando por la interna.

 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. left( cos theta right) right| 0pi= left( cos pi right) left( cos0 right)=1+1=2 int limits2 pi0d varphi= left. varphi right|2 pi0=2 pi endarray


Entonces como resultado obtenemos:

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


La integral de Euler-Poisson se usa a menudo en la teoría de la probabilidad.

Espero que para alguien el artículo sea útil y ayude a comprender algunas técnicas matemáticas :)

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


All Articles