Integral ∫∞−∞exp(−x2)1+x4dx∫−∞∞exp⁡(−x2)1+x4dx\int_{-\infty}^\infty \frac{\exp{(-x^2)}} {1+x^4}x

Tengo una integral que me gustaría algunos consejos sobre cómo encontrar. Es:

Exp ( X 2 ) 1 + X 4 d X

Tengo algo de experiencia con integrales de contorno, así que intenté usar una integral de contorno donde el contorno es un semicírculo de radio R en la mitad superior del plano complejo. Creo que la integral a lo largo del arco se desvanece como R así que utilicé el teorema del residuo y obtengo un valor de π 8 ( porque ( 1 ) pecado ( 1 ) ) . ¡Pero esto no puede ser correcto ya que el valor es negativo mientras que el integrando siempre es positivo! No puedo averiguar qué es lo que va mal.

¿Alguien tiene una idea de cómo evaluar esta integral? La integración del contorno es el método preferido si es posible, pero también estoy abierto a otros métodos. Cualquier sugerencia es apreciada!

Bienvenido a MSE. Buena primera pregunta!
El problema con tu semicírculo es que si X = i y entonces Exp ( X 2 ) / ( 1 + X 4 ) = Exp ( y 2 ) / ( 1 + y 4 ) , que va a + como y + .
Solo una advertencia para cualquiera que intente esto. Parece que las integrales de Fresnel son inevitables. WolframAlpha dará una respuesta exacta en términos de Fresnel S y C . Intentar usar un método similar a esta publicación es una opción, pero conduce a resolver y + y = π / X , lo que de nuevo conduce a las integrales de Fresnel.
Mathematica da la respuesta
π ( ( 2 C ( 2 π ) 1 ) pecado ( 1 ) + ( 1 2 S ( 2 π ) ) porque ( 1 ) ) 2 .
Aquí, S y C son las funciones integrales de Fresnel mencionadas anteriormente.
fracciones parciales descompone la integral en dos partes 1 2 mi X pag ( X 2 ) 1 i X 2 d s + 1 2 mi X pag ( X 2 ) 1 + i X 2 d s
Maple evalúa esto en términos de la función Ira j y la función de Weber mi . Otra razón para pensar que no es "elemental".
Curiosamente parece que + Exp ( X k ) 1 + X 4 d X no converge siempre k es impar. Por lo que vale la pena + Exp ( X 4 ) 1 + X 4 d X = 2 mi Γ ( 5 4 ) Γ ( 3 4 , 1 ) ; dónde Γ ( X , y ) es la función gamma incompleta. Hay patrones similares para k>2 pares.
Gracias por agregar ese enlace, Brian. No pensé en usar la diferenciación debajo del signo integral. De hecho, obtengo su ODE de segundo orden y, por lo tanto, ¡veo por qué surgen las integrales de Fresnel!
Solo por curiosidad, ya que el contorno del semicírculo superior es una mala idea, sería interesante saber si hay algún contorno que permita evaluar tal integral.

Respuestas (2)

mi X 2 1 + X 4 d X = mi X 2 0 mi X 2 t pecado t d t d X = 0 pecado t mi ( 1 + t ) X 2 d X d t
= π 0 pecado t 1 + t d t = 1 + t = X 2 2 π 1 pecado ( X 2 1 ) d X
= 2 π porque 1 1 pecado ( X 2 ) d X 2 π pecado 1 1 porque ( X 2 ) d X
= π porque 1 1 2 S ( 2 π ) 2 π pecado 1 1 2 C ( 2 π ) 2
Dónde S ( X ) y C ( X ) son integrales de Fresnel .

¡Buena solución! El único truco no obvio es reconocer 1 / ( 1 + X 4 ) como 0 Exp ( X 2 t ) s i norte ( t ) d t
Tienes un talento único para resolver integrales.
¡Gracias por las buenas palabras!
Esto demuestra efectivamente la importancia de recordar las transformadas de Laplace al integrar :D

¡La sugerencia de Brian anterior merece una mención ya que pude resolver la integral de esa manera también!

Dejar

I ( t ) = Exp ( t X 2 ) 1 + X 4 d X
para un parámetro t 0 . Diferenciando da
I ( t ) = X 2 Exp ( t X 2 ) 1 + X 4 d X
y una segunda diferenciación da
I ( t ) = X 4 Exp ( t X 2 ) 1 + X 4 d X .
Por eso,
I ( t ) + I ( t ) = Exp ( t X 2 ) = π t
como Brian señaló.

Recordando mis notas de clase, ahora tenemos una EDO lineal no homogénea de segundo orden con coeficientes constantes. La solución se puede escribir

I ( t ) = I C ( t ) + I PAG ( t )
dónde I C ( t ) resuelve la EDO homogénea
I C ( t ) + I C ( t ) = 0.
Dejar I 1 ( t ) = pecado ( t ) y I 2 ( t ) = porque ( t ) sean las dos soluciones de la EDO homogénea. La solución especial I PAG ( t ) ahora se puede encontrar usando el método descrito aquí http://tutorial.math.lamar.edu/Classes/DE/VariationofParameters.aspx

Observando que el Wronskiano para I1 e I2 es -1, la solución particular es

I PAG ( t ) = pecado ( t ) 0 t π tu porque ( tu ) d tu porque ( t ) 0 t π tu pecado ( tu ) d tu
o
I PAG ( t ) = 2 π ( pecado ( t ) 0 t porque ( tu 2 ) d tu porque ( t ) 0 t pecado ( tu 2 ) d tu ) .
Poner todo esto junto da
I ( t ) = A pecado ( t ) + B porque ( t ) + π 2 { C ( 2 t π ) pecado ( t ) S ( 2 t π ) porque ( t ) } ,
dónde C ( X ) y S ( X ) son las integrales de coseno y seno de Fresnel respectivamente.

Las condiciones iniciales son

I ( 0 ) = d X 1 + X 4
y
I ( 0 ) = X 2 1 + X 4 d X .

Usando un contorno en forma de D en el plano complejo, estas integrales se pueden mostrar tanto como iguales π 2 . Por eso B = A = π 2 .

Esto finalmente da

I ( 1 ) = π porque ( 1 ) 1 2 S ( 2 π ) 2 π pecado ( 1 ) 1 2 C ( 2 π ) 2 .