6 votos

Integrales múltiples de Fourier con la función Theta de Heaviside

Quiero evaluar la integral:

$$I=\int_{-\infty}^{\infty}dx_1 \int_{-\infty}^{\infty}dx_2 \ \Theta(x_1-x_2) \ e^{i(ax_1+bx_2)}$$ donde $\Theta(x)$ es la función de Heaviside.

Lo que estaba haciendo ahora era tomar la relación para $\Theta$ : $\Theta (x)=-\frac{1}{2\pi i}\int_{-\infty}^{\infty}d\tau \frac{1}{\tau + i\epsilon} e^{-ix\tau}$ y tengo: $$I= -\frac{1}{2\pi i}\int_{-\infty}^{\infty}dx_1 \int_{-\infty}^{\infty}dx_2\int_{-\infty}^{\infty}d\tau \ \frac{1}{\tau + i\epsilon}e^{-i(\tau-a)x_1}e^{-i(\tau+b)x_2}\\=2\pi i\int_{-\infty}^{\infty}d\tau\ \frac{1}{\tau + i\epsilon}\delta(\tau-a)\delta(\tau+b) =2\pi i\frac{\delta(a+b)}{a + i\epsilon}$$ No sabía si la integral era convergente y podía simplemente intercambiar las integrales, así que lo intenté de otra forma con $X=x_1+x_2$ y $x=x_1-x_2$ : $$I=\frac{1}{2}\int_{-\infty}^{\infty}dX \int_{-\infty}^{\infty}dx \ \Theta(x) \ e^{ia\frac{x+X}{2}+ib\frac{X-x}{2}} \\ =\pi \int_{-\infty}^{\infty}dx \ \Theta(x) e^{-2\pi i\frac{b-a}{4\pi}}\delta(\frac{b+a}{2})$$ Con la transformada de Fourier de la función Heaviside $\int_{-\infty}^{\infty}dk\ \Theta(k)e^{-2\pi i kx}=\frac{1}{2}(\delta(x)-\frac{i}{\pi k}) $ Me sale $$I=\pi \left(2\pi\delta(b-a)-\frac{4 i}{b-a}\right)\delta(a+b)=2\pi^2\delta(a)\delta(b)+2\pi i\frac{\delta(a+b)}{a}$$ Todavía no sé dónde está el $\delta(a)\delta(b)$ debe venir en el primer método. Cuando quiero comprobar que ahora e integrar $I$ en $a$ y $b$ Me sale desde la primera línea: $$\int_{-\infty}^{\infty}da \int_{-\infty}^{\infty}db \ I = \int_{-\infty}^{\infty}dx_1\int_{-\infty}^{\infty}dx_2 \Theta(x_1-x_2) \delta(x_1)\delta(x_2) \\ = \Theta(0)=\frac{1}{2}$$ y del segundo resultado: $$\int_{-\infty}^{\infty}da \int_{-\infty}^{\infty}db \ I=4\pi^2-\int_{-\infty}^{\infty}db \frac{1}{b}=-\infty$$ ¿En qué se equivocó? ¿Es correcta la integral?

EDIT: corregido el error en la derivación a causa del comentario.

2voto

Fredrik Puntos 26

Los productos de las distribuciones son un tema delicado bien conocido, véase por ejemplo ce Puesto de mathoverflow. A no ser que cada factor de distribución dependa de diferentes variables, como por ejemplo $\delta^3(\vec{r})=\delta(x)\delta(y)\delta(x)$ . Así que lo más seguro es escribir la integral de OP en dos variables manifiestamente independientes. Por ejemplo

$$ I(a_1,a_2) ~:=~ \iint_{\mathbb{R}^2} \! \mathrm{d}x^1 \mathrm{d}x^2~ e^{i (a_1 x^1+a_2 x^2)} \theta(x^1-x^2) $$ $$\tag{1}~=~ \iint_{\mathbb{R}^2} \! \mathrm{d}y^1 \mathrm{d}y^2 ~ e^{i (b_1 y^1 + b_2 y^2)} \theta(y^1)~=~I_1(b_1)~I_2(b_2) .$$

Aquí hemos definido una transformación lineal de coordenadas (con determinante unitario para evitar un factor jacobiano por simplicidad)

$$ \tag{2} \begin{bmatrix}y^1 \\ y^2 \end{bmatrix} ~=~\begin{bmatrix} \frac{1}{2} &-\frac{1}{2} \\ 1-t & t \end{bmatrix} \begin{bmatrix}x^1 \\ x^2 \end{bmatrix} $$

de manera que la nueva variable $y^1$ se convierte en el argumento de la función escalonada de Heaviside. Aquí $t\in\mathbb{R}$ es un parámetro libre. Es interesante observar cómo la integral (1) es independiente del valor de este parámetro libre $t$ . La transformación inversa es

$$ \tag{3} \begin{bmatrix}x^1 \\ x^2 \end{bmatrix} ~=~\begin{bmatrix} t&\frac{1}{2} \\ t-1 & \frac{1}{2} \end{bmatrix} \begin{bmatrix}y^1 \\ y^2 \end{bmatrix}. $$

Para que el argumento $a_1 x^1+a_2 x^2=b_1 y^1 + b_2 y^2$ de la exponencial (1) para que permanezca invariante de forma, definimos nuevas variables

$$ \tag{4} \begin{bmatrix}b_1 \\ b_2 \end{bmatrix} ~=~\begin{bmatrix} t& t-1 \\ \frac{1}{2}& \frac{1}{2} \end{bmatrix} \begin{bmatrix}a_1 \\ a_2 \end{bmatrix}. $$

El primer factor integral es

$$\tag{5} I_1(b_1)~:= \int_{\mathbb{R}} \! \mathrm{d}y^1~e^{i b_1 y^1}\theta(y^1) ~=~\frac{i}{b_1+i0^+}~=~P\frac{i}{b_1}+\pi\delta(b_1),$$

Véase el Fórmula Sokhotski-Plemelj . El segundo factor integral es

$$\tag{6} I_2(b_2)~:= \int_{\mathbb{R}} \! \mathrm{d}y^2~e^{i b_2 y^2} ~=~2\pi ~\delta(b_2).$$

Mientras que los dos factores (5) y (6) dependen claramente de $t$ , su producto [=integral de OP (1)] es independiente de $t$ :

$$ I(a_1,a_2)~=~I_1(b_1)~I_2(b_2)~=~\frac{2\pi i}{a_1+i0^+}\delta(a_1+a_2)$$ $$\tag{7} ~=~2\pi i ~P\frac{1}{a_1}\delta(a_1+a_2)+2\pi^2 \delta(a_1)\delta(a_2). $$

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