Konturintegration zur Bewertung eines reellen Integrals

Nov 07 2020

Ich bewerte dieses Integral:

$$\int_{-\infty}^{\infty} \frac{\sin x}{x(x^2+1)^2}\,dx$$

mit der Formel

$$\int_{-\infty}^{\infty} f(x) \sin(sx) dx = 2\pi \sum\text{Re } \text{Res}[f(z) e^{isz}]$$

wobei die Summe über den Resten in der oberen Halbebene liegt.

Da also die einzigen zwei Singularitäten innerhalb der oberen Halbebene liegen $z = 0$ und $z=i$, Ich habe das gefunden

$$\begin{align} 2\pi \sum\text{Re } \text{Res}\left(f(z) e^{isz}\right) &= 2 \pi \left(\text{Re } \text{Res}_{z= 0}\left[\frac{1}{z (z^2+1)^2} e^{i z}\right] + \text{Re } \text{Res}_{z= i}\left[\frac{1}{z (z^2+1)^2} e^{i z}\right]\right) \\\\ &=2 \pi \left(1 + \frac{-3}{4e}\right) \end{align}$$

Ich bin mir ziemlich sicher, dass ich die beiden Reste richtig berechnet habe, da in mathematica

Residue[E^(I z)/(z (z^2 + 1)^2), {z, 0}]

ist $1$ und

Residue[E^(I z)/(z (z^2 + 1)^2), {z, I}]

ist $\frac{-3}{4e}$

Aber das Integral bewerten

Integrate[Sin[x]/(x (x^2 + 1)^2), {x, -Infinity, Infinity}]

mathematica gibt $\pi - \frac{3 \pi}{2e}$.

Ich frage mich, ob dies daran liegt, dass ich irgendwo etwas falsch gemacht habe oder ob es daran liegt, dass mathematica die falsche Antwort gibt.

Dankeschön!

Antworten

1 MarkViola Nov 07 2020 at 04:06

Der Einschluss des Rückstands bei $z=0$das ist nicht richtig. Wir beginnen vielmehr mit dem Schreiben

$$\int_{-\infty}^\infty \frac{\sin(x)}{x(x^2+1)^2}\,dx=\text{Im}\left(\text{PV}\int_{-\infty}^\infty \frac{e^{ix}}{x(x^2+1)^2}\,dx\right)$$

wo der Cauchy-Hauptwert gegeben ist durch

$$\text{PV}\int_{-\infty}^\infty \frac{e^{ix}}{x(x^2+1)^2}\,dx=\lim_{\varepsilon\to 0^+}\int_{|x|>\varepsilon}\frac{e^{ix}}{x(x^2+1)^2}\,dx$$


Als nächstes bewegen wir uns zur komplexen Ebene. Ler$R>1$, $\varepsilon>0$, und $C$ sei die Kontur in der oberen Halbebene, die als $(i)$ die geraden Wege von $-R$ zu $-\varepsilon$ und von $\varepsilon$ zu $R$, $(ii)$ der halbkreisförmige Bogen zentriert bei $z=0$ mit Radius $\varepsilon$ von $-\varepsilon$ zu $\varepsilon$, und $(iii)$ der halbkreisförmige Bogen zentriert bei $z=0$ mit Radius $R$ von $R$ zu $-R$. Beachten Sie, dass$z=0$ ist aus dem durch begrenzten Innenbereich ausgeschlossen $C$.

Dann können wir schreiben

$$\begin{align} \oint_{C}\frac{e^{iz}}{z(z^2+1)^2}\,dz&=\text{PV}\int_{-\infty}^\infty \frac{e^{ix}}{x(x^2+1)^2}\,dx\\\\ &+\int_{\pi}^0 \frac{e^{i\varepsilon e^{i\phi}}}{\varepsilon e^{i\phi}((\varepsilon e^{i\phi})^2+1)^2}\,i\varepsilon e^{i\phi}\,d\phi\\\\ &+\int_0^{\pi} \frac{e^{iR e^{i\phi}}}{R e^{i\phi}((R e^{i\phi})^2+1)^2}\,iR e^{i\phi}\,d\phi\tag1 \end{align}$$

Wie $R\to \infty$, das letzte Integral auf der rechten Seite von $(1)$ nähert sich $0$.

Wie $\varepsilon\to0^+$, das zweite Integral auf der rechten Seite von $(1)$ nähert sich $-i\pi$.

Schon seit $C$ hat die ausgeschlossen $z=0$ist der einzige betroffene Rückstand bei $z=i$. Daher finden wir

$$\text{PV}\int_{-\infty}^\infty \frac{e^{ix}}{x(x^2+1)^2}\,dx=i\pi +\text{Res}\left(\frac{e^{iz}}{z(z^2+1)^2}, z=i\right)\tag2$$

Berechnen Sie nun den Rückstand bei $z=i$ und nimm den Imaginärteil beider Seiten von $(2)$. Kannst du jetzt fertig werden?

MostafaAyaz Nov 07 2020 at 04:05

Schon seit $z=0$ ist ein einzelner Pol von $f(z)$ an der Grenze des oberen Halbkreises sollte es multipliziert werden $\pi$ lieber als $2\pi$.

FelixMarin Nov 07 2020 at 10:09

$\newcommand{\bbx}[1]{\,\bbox[15px,border:1px groove navy]{\displaystyle{#1}}\,} \newcommand{\braces}[1]{\left\lbrace\,{#1}\,\right\rbrace} \newcommand{\bracks}[1]{\left\lbrack\,{#1}\,\right\rbrack} \newcommand{\dd}{\mathrm{d}} \newcommand{\ds}[1]{\displaystyle{#1}} \newcommand{\expo}[1]{\,\mathrm{e}^{#1}\,} \newcommand{\ic}{\mathrm{i}} \newcommand{\mc}[1]{\mathcal{#1}} \newcommand{\mrm}[1]{\mathrm{#1}} \newcommand{\on}[1]{\operatorname{#1}} \newcommand{\pars}[1]{\left(\,{#1}\,\right)} \newcommand{\partiald}[3][]{\frac{\partial^{#1} #2}{\partial #3^{#1}}} \newcommand{\root}[2][]{\,\sqrt[#1]{\,{#2}\,}\,} \newcommand{\totald}[3][]{\frac{\mathrm{d}^{#1} #2}{\mathrm{d} #3^{#1}}} \newcommand{\verts}[1]{\left\vert\,{#1}\,\right\vert}$ $\ds{\bbox[5px,#ffd]{}}$


\begin{align} &\bbox[5px,#ffd]{\int_{-\infty}^{\infty}{\sin\pars{x} \over x\pars{x^{2} + 1}^{2}}\,\dd x} = \Im\int_{-\infty}^{\infty}{\expo{\ic x} - 1\over x\pars{x^{2} + 1}^{2}}\,\dd x \\[5mm] = &\ \Im\braces{2\pi\ic\,\lim_{x \to \ic}\,\totald{}{x} \bracks{\pars{x - \ic}^{2}\,{\expo{\ic x} - 1\over x\pars{x^{2} + 1}^{2}}}} \\[5mm] = &\ 2\pi\,\Re\braces{\lim_{x \to \ic}\,\totald{}{x} \bracks{{\expo{\ic x} - 1\over x\pars{x + \ic}^{2}}}} \\[5mm] = &\ 2\pi\,\Re\bracks{\lim_{x \to \ic}\ {\expo{\ic x}\pars{\ic x^{2} - 4x - \ic} + 3x + \ic \over x^{2}\,\pars{x + \ic}^{3}}} \\[5mm] = &\ \bbx{\pi - {3\pi \over 2\expo{}}} \approx 1.4080 \\ & \end{align}