Integral de función sinc sobre intervalo acotado

Estoy tratando de determinar si la integral

0 1 pecado X X d X
puede calcularse analíticamente.

Soy consciente de la definición de la función integral del seno. Si ( X ) , pero no he podido encontrar una referencia que indique si tal integral se puede calcular sin recurrir a métodos numéricos.

Los métodos habituales empleados para calcular la integral impropia sobre la línea real, como mediante la integración de contorno o definiendo

GRAMO ( t ) = 0 1 pecado X X mi t X d X
no parecen trabajar aquí.

¿Puede alguien darme una pista aquí?

Qué hay de malo en Si ? Sí, debe calcularse utilizando "métodos numéricos". Pero también todo lo demás. pecado se calcula mediante una serie de potencias. π se calcula con una serie infinita. No veo cómo son diferentes a algo como Si .

Respuestas (3)

Dado que el rango de integración está restringido a [ 0 , 1 ] puede hacer uso de Taylor Series para la función seno:

pecado ( X ) = k = 0 + ( 1 ) k X 2 k + 1 ( 2 k + 1 ) !

Por eso

k = 0 + ( 1 ) k ( 2 k + 1 ) ! 0 1 X 2 k + 1 X   d X

La integración es trivial y te da

k = 0 + ( 1 ) k ( 2 k + 1 ) ! 1 2 k + 1

Entonces puedes ir ad libitum con numérico.

Note que esa serie no es más que la expresión de la serie de \text{Si}(1) (SineIntegral).

El valor numérico es, en aras de la exhaustividad, 0.946083 ( . . . )

El valor de la serie, para algunos valores de k es:

k = 2 0.94611
k = 5 0.946083

Se necesitan pocos términos para obtener una alta precisión numérica.

0 1 pecado X X d X = 0 1 ( 1 X 2 3 ! + X 4 5 ! . . . ) d X = 1 1 3 3 ! + 1 5 5 ! . . .

Usando la respuesta de @Turing, podría ser divertido saber la cantidad de términos que se agregarán para una precisión determinada.

Escribiendo

k = 0 + ( 1 ) k ( 2 k + 1 ) ( 2 k + 1 ) ! = k = 0 pag ( 1 ) k ( 2 k + 1 ) ( 2 k + 1 ) ! + k = pag + 1 + ( 1 ) k ( 2 k + 1 ) ( 2 k + 1 ) !
como es una serie alterna, buscamos pag tal que
1 ( 2 pag + 3 ) ( 2 pag + 3 ) ! ϵ ( 2 pag + 3 ) ( 2 pag + 3 ) ! 1 ϵ
Planteamos simplemente el problema (ya que buscamos una aproximación) a
( 2 pag + 4 ) ! 1 ϵ

Si echas un vistazo aquí , encontrarás una excelente aproximación del factorial inverso. Aplicado al presente caso, esto daría

pag registro ( 2 π ϵ ) 2 W ( registro ( 2 π ϵ ) mi ) 9 4
dónde W ( . ) es la función de Lambert.

Dado que el argumento es bastante grande, puede utilizar la aproximación

W ( t ) = L 1 L 2 + L 2 L 1 + ( L 2 2 ) L 2 2 L 1 2 + ( 2 L 2 2 9 L 2 + 6 ) L 2 6 L 1 3 + . . .
dónde L 1 = registro ( t ) y L 2 = registro ( L 1 ) y, seguro, usarás pag .

Para ϵ = 10 10 , esto daría pag = 4.59018 es decir 6 . Informática, esto daría

1 13   13 ! = 1.235 × 10 11