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

Ich habe ein Integral, für das ich einen Ratschlag hätte, wie ich es finden kann. Es ist:

exp ( X 2 ) 1 + X 4 D X

Ich habe einige Erfahrung mit Konturintegralen, also habe ich versucht, ein Konturintegral zu verwenden, bei dem die Kontur ein Halbkreis mit dem Radius R in der oberen Hälfte der komplexen Ebene ist. Ich denke, das Integral entlang des Bogens verschwindet als R also habe ich den Residuensatz verwendet und bekomme einen Wert von π 8 ( cos ( 1 ) Sünde ( 1 ) ) . Das kann aber nicht stimmen, da der Wert negativ ist, während der Integrand immer positiv ist! Ich kann nicht herausfinden, was falsch läuft.

Hat jemand eine Idee, wie man dieses Integral auswerten kann? Die Konturintegration ist nach Möglichkeit die bevorzugte Methode, aber ich bin auch offen für andere Methoden. Jeder Vorschlag ist willkommen!

Willkommen bei MSE. Schöne erste Frage!
Das Problem mit Ihrem Halbkreis ist, dass if X = ich j Dann exp ( X 2 ) / ( 1 + X 4 ) = exp ( j 2 ) / ( 1 + j 4 ) , was geht + als j + .
Nur eine Warnung an alle, die das versuchen. Es sieht so aus, als ob die Fresnel-Integrale unvermeidlich sind. WolframAlpha wird eine genaue Antwort in Bezug auf den Fresnel geben S Und C . Der Versuch, eine ähnliche Methode wie in diesem Beitrag zu verwenden, ist eine Option, führt jedoch zur Lösung j + j = π / X , was wiederum zu den Fresnel-Integralen führt.
Mathematica gibt die Antwort
π ( ( 2 C ( 2 π ) 1 ) Sünde ( 1 ) + ( 1 2 S ( 2 π ) ) cos ( 1 ) ) 2 .
Hier, S Und C sind die oben erwähnten Fresnel-Integralfunktionen.
Partialbruch zerfällt ganzzahlig in zwei Teile 1 2 e X P ( X 2 ) 1 ich X 2 D S + 1 2 e X P ( X 2 ) 1 + ich X 2 D S
Maple wertet dies anhand der Anger-Funktion aus J und die Weber-Funktion E . Ein weiterer Grund zu glauben, dass es nicht "elementar" ist.
Merkwürdigerweise sieht es so aus + exp ( X k ) 1 + X 4 D X konvergiert nicht immer k ist ungerade. Für was es wert ist + exp ( X 4 ) 1 + X 4 D X = 2 e Γ ( 5 4 ) Γ ( 3 4 , 1 ) ; Wo Γ ( X , j ) ist die unvollständige Gammafunktion. Es gibt ähnliche Muster für gerade k>2.
Danke, dass du diesen Link hinzugefügt hast, Brian - ich habe nicht daran gedacht, die Differenzierung unter dem Integralzeichen zu verwenden. Ich bekomme tatsächlich Ihre ODE zweiter Ordnung und sehe daher, warum die Fresnel-Integrale auftauchen!
Nur aus Neugier - da die obere Halbkreiskontur eine schlechte Idee ist, wäre es interessant zu wissen, ob es eine Kontur gibt, die die Auswertung eines solchen Integrals ermöglichen würde.

Antworten (2)

e X 2 1 + X 4 D X = e X 2 0 e X 2 T Sünde T D T D X = 0 Sünde T e ( 1 + T ) X 2 D X D T
= π 0 Sünde T 1 + T D T = 1 + T = X 2 2 π 1 Sünde ( X 2 1 ) D X
= 2 π cos 1 1 Sünde ( X 2 ) D X 2 π Sünde 1 1 cos ( X 2 ) D X
= π cos 1 1 2 S ( 2 π ) 2 π Sünde 1 1 2 C ( 2 π ) 2
Wo S ( X ) Und C ( X ) sind Fresnel-Integrale .

Schöne Lösung! Der einzige nicht offensichtliche Trick ist das Erkennen 1 / ( 1 + X 4 ) als 0 exp ( X 2 T ) S ich N ( T ) D T
Sie haben ein einzigartiges Talent zum Lösen von Integralen.
Danke für die netten Worte!
Dies zeigt effektiv, wie wichtig es ist, sich an Laplace-Transformationen zu erinnern, während :D integriert wird

Brians obiger Vorschlag verdient eine Erwähnung, da ich das Integral auch auf diese Weise lösen konnte!

Lassen

ICH ( T ) = exp ( T X 2 ) 1 + X 4 D X
für einen Parameter T 0 . Differenzieren gibt
ICH ' ( T ) = X 2 exp ( T X 2 ) 1 + X 4 D X
und eine zweite Differenzierung ergibt
ICH ( T ) = X 4 exp ( T X 2 ) 1 + X 4 D X .
Somit,
ICH ( T ) + ICH ( T ) = exp ( T X 2 ) = π T
wie Brian betonte.

In Erinnerung an meine Vorlesungsunterlagen haben wir jetzt eine inhomogene lineare GDE zweiter Ordnung mit konstanten Koeffizienten. Die Lösung kann geschrieben werden

ICH ( T ) = ICH C ( T ) + ICH P ( T )
Wo ICH C ( T ) löst die homogene DGL
ICH C ( T ) + ICH C ( T ) = 0.
Lassen ICH 1 ( T ) = Sünde ( T ) Und ICH 2 ( T ) = cos ( T ) seien die beiden Lösungen der homogenen ODE. Die besondere Lösung ICH P ( T ) kann nun mit der hier beschriebenen Methode http://tutorial.math.lamar.edu/Classes/DE/VariationofParameters.aspx gefunden werden

Beachten Sie, dass der Wronskian für I1 und I2 -1 ist, die spezielle Lösung ist

ICH P ( T ) = Sünde ( T ) 0 T π u cos ( u ) D u cos ( T ) 0 T π u Sünde ( u ) D u
oder
ICH P ( T ) = 2 π ( Sünde ( T ) 0 T cos ( u 2 ) D u cos ( T ) 0 T Sünde ( u 2 ) D u ) .
Das alles zusammen ergibt
ICH ( T ) = A Sünde ( T ) + B cos ( T ) + π 2 { C ( 2 T π ) Sünde ( T ) S ( 2 T π ) cos ( T ) } ,
Wo C ( X ) Und S ( X ) sind die Fresnel-Kosinus- bzw. Sinus-Integrale.

Die Anfangsbedingungen sind

ICH ( 0 ) = D X 1 + X 4
Und
ICH ' ( 0 ) = X 2 1 + X 4 D X .

Unter Verwendung einer D-förmigen Kontur in der komplexen Ebene können diese Integrale als beide gleich gezeigt werden π 2 . Somit B = A = π 2 .

Das gibt endlich

ICH ( 1 ) = π cos ( 1 ) 1 2 S ( 2 π ) 2 π Sünde ( 1 ) 1 2 C ( 2 π ) 2 .