10.4 Fórmula de inversão

A prova do Teorema de Unicidade que demos na seção anterior é muito instrutiva. Ela está baseada no fato de que, se duas distribuições são indistinguíveis quando testadas em funções trigonométricas, então elas são as mesmas.

Entretanto, aquela prova não é auto-contida pois usa a versão trigonométrica do Teorema de Weierstrass. Daremos nesta seção uma prova auto-contida, que terá seus pontos altos como transformada de Fourier inversa e cálculo da integral de Dirichlet sem utilizar Análise Complexa.

Enunciamos e provamos agora a fórmula de inversão, que tem o Teorema de Unicidade como corolário. Se bem que nem o enunciado nem as técnicas de prova serão usadas no restante desse livro, sua leitura envolve um uso mais engenhoso da integral de Lebesgue, o que pode causar certa satisfação.

Teorema 10.45 (Fórmula de inversão).

Sejam XX uma variável aleatória e a<ba<b dois pontos de continuidade de FXF_{X}. Então

FX(b)FX(a)=limc+12πabc+ceixtφX(t)dtdx.F_{X}(b)-F_{X}(a)=\lim_{c\to+\infty}\frac{1}{2\pi}\int_{a}^{b}\int_{-c}^{+c}e^% {-ixt}\varphi_{X}(t)\,\mathrm{d}t\,\mathrm{d}x.

Informalmente, se XX tem densidade suave e com cauda leve,

φX(t)=+eitxfX(x)dx\displaystyle\varphi_{X}(t)=\int_{-\infty}^{+\infty}e^{itx}f_{X}(x)\,\mathrm{d}x

cuja transformada de Fourier inversa1717 17 A transformada de Fourier inversa não será usada neste livro, é mencionada apenas para explicar de onde veio a fórmula de inversão. Em Análise, o sinal de menos no expoente aparece na transformada e não na sua inversa, e usa-se 2π\sqrt{2\pi} em ambas ao invés de 2π2\pi na inversa.  é dada por

fX(x)=12π+eixtφX(t)dt.f_{X}(x)=\frac{1}{2\pi}\int_{-\infty}^{+\infty}e^{-ixt}\varphi_{X}(t)\,\mathrm% {d}t.

Integrando em [a,b][a,b] obtemos

FX(b)FX(a)=abfX(x)dx=12πab+eixtφX(t)dtdx,F_{X}(b)-F_{X}(a)=\int_{a}^{b}f_{X}(x)\,\mathrm{d}x=\frac{1}{2\pi}\int_{a}^{b}% \int_{-\infty}^{+\infty}e^{-ixt}\varphi_{X}(t)\,\mathrm{d}t\,\mathrm{d}x,

explicando a fórmula de inversão. Havendo motivado a fórmula, passemos à demonstração, sem supor suavidade de nenhum tipo.

Demonstração.

Sejam <a<b<-\infty<a<b<\infty. Para cada c>0c>0, definimos

Ac=12πabc+ceixtφX(t)dtdx.A_{c}=\frac{1}{2\pi}\int_{a}^{b}\int_{-c}^{+c}e^{-ixt}\varphi_{X}(t)\,\mathrm{% d}t\,\mathrm{d}x.

Usando o Teorema de Fubini para inverter as integrais,

Ac=12π𝔼[c+cabeixteitXdxdt],A_{c}=\frac{1}{2\pi}\mathbb{E}\Big{[}\int_{-c}^{+c}\int_{a}^{b}e^{-ixt}e^{itX}% \,\mathrm{d}x\,\mathrm{d}t\Big{]},

que se justifica porque abcc|eitxeitX|dtdx=2c(ba)<\int_{a}^{b}\int_{-c}^{c}|e^{-itx}e^{itX}|\mathrm{d}t\mathrm{d}x=2c(b-a)<\infty.

A partir da expressão acima, definimos

gc(z)=c+cabeixteitzdxdt.\displaystyle g_{c}(z)=\int_{-c}^{+c}\int_{a}^{b}e^{-ixt}e^{itz}\,\mathrm{d}x% \,\mathrm{d}t.

Integrando, expandindo, usando que o cosseno é par e substituindo, obtemos

gc(z)=c(za)+c(za)senvvdvc(zb)+c(zb)senvvdv.\displaystyle g_{c}(z)=\int_{-c(z-a)}^{+c(z-a)}\frac{\mathop{\mathrm{sen}}% \nolimits v}{v}\,\mathrm{d}v-\int_{-c(z-b)}^{+c(z-b)}\frac{\mathop{\mathrm{sen% }}\nolimits v}{v}\,\mathrm{d}v.

Observe que a função h(r,s)=rssenvvdvh(r,s)=\int_{r}^{s}\frac{\mathop{\mathrm{sen}}\nolimits v}{v}\,\mathrm{d}v é limitada e tende a algum número positivo λ>0\lambda>0 quando s+s\to+\infty e rr\to-\infty. Mais abaixo vamos descobrir o valor desse limite λ\lambda.

De volta às funções gcg_{c}, separando em casos, obtemos

limc+gc(z)=2λ[𝟙(a,b)(z)+12𝟙{a,b}(z)].\lim_{c\to+\infty}g_{c}(z)=2\lambda\Big{[}\mathds{1}_{(a,b)}(z)+\tfrac{1}{2}% \mathds{1}_{\{a,b\}}(z)\Big{]}.

Como gcg_{c} é limitada, pelo Teorema da Convergência Dominada,

limc+Ac=limc+12π𝔼[gc(X)]=λπ(a<X<b)+12λπ(X{a,b}).\lim_{c\to+\infty}A_{c}=\lim_{c\to+\infty}\tfrac{1}{2\pi}\mathbb{E}[g_{c}(X)]=% \tfrac{\lambda}{\pi}\mathbb{P}(a<X<b)+\tfrac{1}{2}\tfrac{\lambda}{\pi}\mathbb{% P}(X\in\{a,b\}).

Repare que o enunciado do teorema é para pontos aa e bb onde FXF_{X} é contínua. Usando o fato de que c+csenvvdvπ\int_{-c}^{+c}\frac{\mathop{\mathrm{sen}}\nolimits v}{v}\mathrm{d}v\to\pi, isso completa a prova do teorema.

Provaremos agora que λ=π\lambda=\pi, isto é, c+csenvvdvπ\int_{-c}^{+c}\frac{\mathop{\mathrm{sen}}\nolimits v}{v}\mathrm{d}v\to\pi. Esta integral, condicionalmente convergente, é conhecida como integral de Dirichlet.

Consideramos o caso X𝒩(0,1)X\sim\mathcal{N}(0,1), que, convenientemente, tem como função característica um múltiplo da sua densidade. Pela fórmula de inversão,

λπ[FX(b)FX(a)]\displaystyle\tfrac{\lambda}{\pi}\big{[}F_{X}(b)-F_{X}(a)\big{]} =limc+12πabc+ceixtφX(t)dtdx.\displaystyle=\lim_{c\to+\infty}\frac{1}{2\pi}\int_{a}^{b}\int_{-c}^{+c}e^{-% ixt}\varphi_{X}(t)\,\mathrm{d}t\,\mathrm{d}x.
=limc+ab12π(c+ceitxet2/22πdt)dx\displaystyle=\lim_{c\to+\infty}\int_{a}^{b}\frac{1}{\sqrt{2\pi}}\Big{(}\int_{% -c}^{+c}e^{-itx}\frac{e^{-t^{2}/2}}{\sqrt{2\pi}}\,\mathrm{d}t\Big{)}\,\mathrm{% d}x
=ab12π(+eitxet2/22πdt)dx\displaystyle=\int_{a}^{b}\frac{1}{\sqrt{2\pi}}\Big{(}\int_{-\infty}^{+\infty}% e^{-itx}\frac{e^{-t^{2}/2}}{\sqrt{2\pi}}\,\mathrm{d}t\Big{)}\,\mathrm{d}x
=ab12πex2/2dx\displaystyle=\int_{a}^{b}\frac{1}{\sqrt{2\pi}}e^{-x^{2}/2}\,\mathrm{d}x
=FX(b)FX(a),\displaystyle=F_{X}(b)-F_{X}(a),

donde concluímos que λ=π\lambda=\pi. Nas igualdades acima usamos o Teorema da Convergência Dominada (a densidade é integrável e |eitx|=1|e^{-itx}|=1), e o valor da integral correspondente à função característica da normal.

Essa prova de que c+csenvvdvπ\int_{-c}^{+c}\frac{\mathop{\mathrm{sen}}\nolimits v}{v}\mathrm{d}v\to\pi ilustra a técnica de aceitar um resultado parcial com constante indeterminada e descobrir a constante após aplicar o argumento a um caso concreto. Outra prova de c+csenvvdvπ\int_{-c}^{+c}\frac{\mathop{\mathrm{sen}}\nolimits v}{v}\mathrm{d}v\to\pi, que tampouco usa Análise Complexa, consiste em escrever

0+csenvvdv\displaystyle\int_{0}^{+c}\tfrac{\mathop{\mathrm{sen}}\nolimits v}{v}\mathrm{d}v =0+c0+exvsenvdxdv=0+0+cexvsenvdvdx\displaystyle=\int_{0}^{+c}\int_{0}^{+\infty}e^{-xv}\mathop{\mathrm{sen}}% \nolimits v\,\mathrm{d}x\mathrm{d}v=\int_{0}^{+\infty}\int_{0}^{+c}e^{-xv}% \mathop{\mathrm{sen}}\nolimits v\,\mathrm{d}v\mathrm{d}x
=0+1ecx(cosc+xsenc)1+x2dx0+11+x2dx=π2\displaystyle=\int_{0}^{+\infty}\frac{1-e^{-cx}(\cos c+x\mathop{\mathrm{sen}}% \nolimits c)}{1+x^{2}}\mathrm{d}x\to\int_{0}^{+\infty}\frac{1}{1+x^{2}}\mathrm% {d}x=\frac{\pi}{2}

quando c+c\to+\infty. Comutar as integrais é justificado pelo Teorema de Fubini, pois 0c|senvv|dv1+[logc]+<\int_{0}^{c}|\frac{\mathop{\mathrm{sen}}\nolimits v}{v}|\mathrm{d}v\leqslant 1% +[\log c]^{+}<\infty, enquanto o limite é justificado pelo Teorema da Convergência Dominada pois o numerador é dominado por 1+1+e11+1+e^{-1} para todo x0x\geqslant 0 e c1c\geqslant 1.

Assim completamos a prova de que c+csenvvdvπ\int_{-c}^{+c}\frac{\mathop{\mathrm{sen}}\nolimits v}{v}\mathrm{d}v\to\pi e, com isso, provamos a fórmula de inversão. ∎