4 votos

Integral con la delta de Dirac (a mí o a wolfram mathematica?)

Traté de calcular con Wolfram Mathematica la siguiente integral $$I=\int_0^\pi\int_{-\infty}^\infty x e^{-jx\cos(\theta-\varphi)}f(\alpha\cos(\theta-\psi))\mathrm \, dx\mathrm \, d\theta$$ donde $-\pi\leq\varphi,\psi\leq\pi$.

Suponiendo que el interior de la intergral es, por definición, la derivada de la función delta de Dirac que obtengo: $$I=2\pi j\int_0^\pi\delta'(\cos(\theta-\varphi))f(\alpha\cos(\theta-\psi)) \, \mathrm \, d\theta$$ Wolfram Mathematica me dice que $$I_\mathrm{wolfram}=0$$ Entonces traté de hacerlo a mano. Si yo uso la definición de una función delta: $$ \begin{eqnarray} I_\mathrm{me}=\delta^{'} (\cos(\theta-\varphi))&=&\left[ \sum_i \frac{\delta(\theta-\varphi-\frac{\pi}{2}-i\pi)}{|\sin(\frac{\pi}{2}+i\pi)|}\right]^{'}=\sum_i \delta^{'}\left(\theta-\left(\varphi+\frac{\pi}{2}+i\pi\right)\right) \end{eqnarray} $$ $$ \begin{eqnarray} I_\mathrm{me}&=&2\pi j\sum_i\int_0^\pi\delta'\left(\theta-\left(\varphi+\frac{\pi}{2}+i\pi\right)\right)f(\alpha\cos(\theta-\psi)) \, \mathrm d\theta=\\ &=&-2\pi j\int_0^\pi\delta\left(\theta-\left(\varphi+\frac{\pi}{2}\right)\right)f'(\alpha\cos(\theta-\psi)) \, \mathrm d\theta=\\ &=&-2\pi j f'(\alpha\cos(\varphi+\frac{\pi}{2}-\psi))u(\pi -2 \varphi ) u(2 \varphi +\pi ) \end{eqnarray} $$ donde $u(\cdot)$ es una función de Heaviside. Y es que en realidad no $0$. Así que ¿dónde está el truco?

1voto

Felix Marin Puntos 32763

\begin{align} &{\rm I}\left(\varphi,\psi\right) \equiv \int_{0}^{\pi}{\rm d}\theta\ \ {\rm f}\left(\alpha\cos\left(\theta - \psi\right)\right) \int_{-\infty}^{\infty}{\rm d}x\ \ x\, {\rm e}^{-{\rm i}x\cos\left(\theta - \varphi\right)} \\[3mm]&= \int_{0}^{\pi}{\rm d}\theta\ \ {\rm f}\left(\alpha\cos\left(\theta - \psi\right)\right)\left\lbrack% {2\pi \over {\rm i}\sin\left(\theta - \psi\right)}\, {\partial \over \partial\theta}\ \overbrace{\quad% \int_{-\infty}^{\infty}{{\rm d}x \over 2\pi}\, {\rm e}^{-{\rm i}x\cos\left(\theta - \varphi\right)}\quad} ^{\delta\left(\cos\left(\theta - \varphi\right)\right)} \right\rbrack \\[3mm]&= 2\pi{\rm i}\int_{0}^{\pi}{\rm d}\theta\, \delta\left(\cos\left(\theta - \varphi\right)\right)\, {\partial \over \partial\theta}\left\lbrack% {{\rm f}\left(\alpha\cos\left(\theta - \psi\right)\right) \over \sin\left(\theta - \psi\right)} \right\rbrack = 2\pi{\rm i}\int_{0}^{\pi}\!\!\!{\rm d}\theta\, \delta\left(\cos\left(\theta - \varphi\right)\right)\, {\partial{\rm F}\left(\alpha,\theta - \psi\right) \over \partial\theta} \end{align}

$\displaystyle{% \mbox{donde}\quad {\rm F}\left(\alpha,\phi\right) \equiv {{\rm f}\left(\alpha\cos\left(\phi\right)\right) \más \sin\left(\phi\right)}} $

\begin{align} \begin{array}{rcl} \sin\left(\theta_{\ell}\left(\varphi\right) - \varphi\right) & = & \left(-1\right)^{\ell} \\ \sin\left(\theta_{\ell}\left(\varphi\right) - \psi\right) & = & \left(-1\right)^{\ell}\cos\left(\varphi - \psi\right) \\ \cos\left(\theta_{\ell}\left(\varphi\right) - \psi\right) & = & \left(-1\right)^{\ell + 1}\sin\left(\varphi - \psi\right) \end--&\begin{align} {\rm I}\left(\varphi,\psi\right) &= -2\pi{\rm i}\sum_{\ell = -\infty}^{\infty} {\partial{\rm F}\left(\alpha,\theta_{\ell}\left(\varphi\right) - \psi\right) \over \partial\psi} \int_{0}^{\pi}{\rm d}\theta\, \delta\left(\theta - \theta_{\ell}\left(\varphi\right)\right) \\[3mm]&= -2\pi{\rm i}\!\!\!\!\!\!\!\!\!\! \sum_{% \ell = -\infty \atop {0\ <\ \theta_{\ell}\left(\varphi\right)\ <\ \pi\vphantom{\Huge A^{A}}}} ^{\infty}\!\!\!\!\!\!\!\!\!\! {\partial{\rm F}\left(\alpha,\theta_{\ell}\left(\varphi\right) - \psi\right) \over \partial\psi} = -2\pi{\rm i}\!\!\!\!\!\!\!\!\!\! \sum_{% \ell = -\infty \atop {\left\vert\ell + \varphi/\pi\right\vert\ <\ 1/2\vphantom{\Huge A}}} ^{\infty}\!\!\!\!\!\!\!\!\!\! {\partial{\rm F}\left(\alpha,\varphi - \psi + \ell\pi + \pi/2\right) \over \partial\psi} \end\begin{array}{rcl} \cos\left(\phi + \ell\pi + {\pi \over 2}\right) & = & -\sin\left(\ell\pi + \phi\right) = \left(-1\right)^{\ell + 1}\sin\left(\phi\right) \\ \sin\left(\phi + \ell\pi + {\pi \over 2}\right) & = & \phantom{-}\cos\left(\ell\pi + \phi\right) = \left(-1\right)^{\ell\phantom{+1}}\cos\left(\phi\right) \end\begin{align} ---&-------------------------------------- \\ {\rm I}\left(\varphi,\psi\right) &= \color{#ff0000}{-2\pi{\rm i}\!\!\!\!\!\!\!\!\!\! \sum_{% \ell = -\infty \atop {\left\vert\ell + \varphi/\pi\right\vert\ <\ 1/2\vphantom{\Huge A}}} ^{\infty}\!\!\!\!\!\!\!\!\!\!\left(-1\right)^{\ell}\, {\partial \over \partial\psi}\left\lbrack% {{\rm f}\left(\left(-1\right)^{\ell + 1}\alpha \sin\left(\varphi - \psi\right)\right) \over \cos\left(\varphi - \psi\right)} \right\rbrack = {\rm I}_{1}\left(\phi,\psi\right) + {\rm I}_{2}\left(\phi,\psi\right)} \\[5mm] {\rm I}_{1}\left(\varphi,\psi\right) &= \color{#ff0000}{-2\pi{\rm i}\alpha\!\!\!\!\!\!\!\!\!\! \sum_{% \ell = -\infty \atop {\left\vert\ell + \varphi/\pi\right\vert\ <\ 1/2\vphantom{\Huge A}}} ^{\infty}\!\!\!\!\!\!\!\!\!\! {\rm f}'\left(\vphantom{\LARGE A}\left(-1\right)^{\ell + 1}\alpha \sin\left(\varphi - \psi\right)\right)\sec\left(\varphi - \psi\right)} \\[3mm] {\rm I}_{2}\left(\varphi,\psi\right) &= \color{#ff0000}{2\pi{\rm i}\!\!\!\!\!\!\!\!\!\! \sum_{% \ell = -\infty \atop {\left\vert\ell + \varphi/\pi\right\vert\ <\ 1/2\vphantom{\Huge A}}} ^{\infty}\!\!\!\!\!\!\!\!\!\!\left(-1\right)^{\ell} {\rm f}\left(\vphantom{\LARGE A}\left(-1\right)^{\ell + 1}\alpha \sin\left(\varphi - \psi\right)\right)\sec\left(\varphi - \psi\right) \tan\left(\varphi - \psi\right)} \end------------------- \\[3mm] {\rm I}\left(\varphi,\psi\right) &= -2\pi{\rm i}\,{\partial \over \partial\psi} \int_{0}^{\pi}{\rm d}\theta\, {\rm F}\left(\alpha,\theta - \psi\right) \delta\left(\cos\left(\theta - \varphi\right)\right) \\[3mm]&= -2\pi{\rm i}\,{\partial \over \partial\psi} \int_{0}^{\pi}{\rm d}\theta\, {\rm F}\left(\alpha,\theta - \psi\right) \sum_{\ell = -\infty}^{\infty} {\delta\left(\theta - \theta_{\ell}\left(\varphi\right)\right) \over \left\vert -\sin\left(\theta_{\ell}\left(\varphi\right) - \varphi\right) \right\vert} \end{align} donde $$ \theta_{\ell}\left(\varphi\right) \equiv \varphi + \left(2\ell + 1\right)\pi/2. \quad\mbox{Aviso de que}\quad \left\vert% -#-#-{array}\right. $$ Entonces -#-#-{align}

$$ \left\lbrace% -#-#-{array}\right. $$

-#-#-{align}

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X