Integral der Sinc-Funktion über ein begrenztes Intervall

Ich versuche festzustellen, ob das Integral

0 1 Sünde X X D X
kann analytisch berechnet werden.

Die Definition der Sinusintegralfunktion ist mir bekannt Si ( X ) , aber ich konnte keine Referenz finden, die besagt, ob ein solches Integral berechnet werden kann, ohne auf numerische Methoden zurückzugreifen.

Die üblichen Methoden zur Berechnung des uneigentlichen Integrals über der realen Linie, z. B. durch Konturintegration oder durch Definieren

G ( T ) = 0 1 Sünde X X e T X D X
scheinen hier nicht zu funktionieren.

Kann mir hier jemand einen Tipp geben?

Was ist falsch mit Si ? Ja, es muss mit "numerischen Methoden" berechnet werden. Aber alles andere auch. Sünde wird mit einer Potenzreihe berechnet. π wird mit einer unendlichen Reihe berechnet. Ich sehe nicht, wie sie sich von so etwas unterscheiden Si .

Antworten (3)

Da beschränkt sich der Integrationsbereich auf [ 0 , 1 ] Sie können die Taylor-Reihe für die Sinusfunktion verwenden:

Sünde ( X ) = k = 0 + ( 1 ) k X 2 k + 1 ( 2 k + 1 ) !

Somit

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

Die Integration ist trivial und gibt Ihnen

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

Dann können Sie mit numerisch nach Belieben gehen.

Beachten Sie, dass diese Reihe nichts anderes ist als der Ausdruck der Reihe von \text{Si}(1) (SineIntegral).

Der Zahlenwert ist der Vollständigkeit halber 0,946083 ( . . . )

Der Wert der Reihe, für einige Werte von k Ist:

k = 2 0,94611
k = 5 0,946083

Es braucht wenige Terme, um eine hohe numerische Genauigkeit zu erhalten.

0 1 Sünde X X D X = 0 1 ( 1 X 2 3 ! + X 4 5 ! . . . ) D X = 1 1 3 3 ! + 1 5 5 ! . . .

Mit der Antwort von @Turing könnte es interessant sein, die Anzahl der Begriffe zu kennen, die für eine bestimmte Genauigkeit hinzugefügt werden müssen.

Schreiben

k = 0 + ( 1 ) k ( 2 k + 1 ) ( 2 k + 1 ) ! = k = 0 P ( 1 ) k ( 2 k + 1 ) ( 2 k + 1 ) ! + k = P + 1 + ( 1 ) k ( 2 k + 1 ) ( 2 k + 1 ) !
Da es sich um alle alternierenden Reihen handelt, suchen wir nach P so dass
1 ( 2 P + 3 ) ( 2 P + 3 ) ! ϵ ( 2 P + 3 ) ( 2 P + 3 ) ! 1 ϵ
Wir gehen einfach dem Problem (da wir nach einer Annäherung suchen) zu
( 2 P + 4 ) ! 1 ϵ

Wenn Sie hier nachsehen , finden Sie eine hervorragende Annäherung an die inverse Fakultät. Übertragen auf den vorliegenden Fall ergäbe sich dies

P Protokoll ( 2 π ϵ ) 2 W ( Protokoll ( 2 π ϵ ) e ) 9 4
Wo W ( . ) ist die Lambert-Funktion.

Da das Argument ziemlich groß ist, können Sie die Näherung verwenden

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 + . . .
Wo L 1 = Protokoll ( T ) Und L 2 = Protokoll ( L 1 ) und sicher werden Sie verwenden P .

Für ϵ = 10 10 , würde dies geben P = 4.59018 das heißt 6 . Rechnen würde das geben

1 13   13 ! = 1.235 × 10 11