Processing math: 100%
Pular para o conteúdo principal

Destaques

O que é mais rápido? Subir ou descer?

Suponha que uma bola de massa m seja lançada verticalmente para cima. Se desconsiderássemos a resistência do ar, teríamos que o tempo de subida seria exatamente igual ao tempo de descida, como já bem sabemos dos estudos básicos de cinemática do ensino médio. Mas e se a resistência for levada em conta, isso muda? Faça sua aposta!  Se considerarmos que a resistência do ar possa ser modelada como um fator proporcional à velocidade, ou seja, a força de resistência do ar seja de f_v = -bv , então podemos escrever a equação diferencial de acordo com a lei de Newton, da seguinte forma:  m\frac{dv}{dt} = -mg-bv
  Perceba que essa mesma equação é válida para o movimento de subida e descida. A fim de determinarmos a velocidade, que é a função que precisamos descobrir dada a equação diferencial anterior, resolvemos a equação diferencial  \frac{dv}{dt} +\frac{b}{m}v = -g
  A solução dessa equação é dada pela soma de duas soluções: a particular...

A integral da Função sinc{(x)}

Talvez você já tenha visto ou ouvido falar da função seno cardinal ou seno por argumento, ou talvez nem saiba que isso exista. Mas agora, se não sabia, vai saber. A função seno cardinal é definida como sinc{(x)}=\frac{ {sin{(x)}}}{x}
O gráfico dessa função é mostrado na figura 1. 


Figura1: A função seno cardinal

Como você pode imaginar, a função sinc{(x)} tende a zero conforme x\rightarrow\infty. No entanto, conforme x\rightarrow0, a função seno cardinal possui uma propriedade interessante. Na verdade, o limite da função seno cardinal quando x\rightarrow0 é um limite fundamental na matemática, representado por \lim _{ x\rightarrow 0 }{ \frac { \sin { (x) }  }{ x }  } =1
Como se pode ver no gráfico, a função tem um "pico" no ponto zero, e esse é o maior valor para a função seno cardinal. Os lóbulos restantes possuem sempre valores menores, e esse decaimento é "ditado" pela função f(x)=\frac{ {1}}{x}. Essa função representa um papel importante no estudo de sinais de comunicação por exemplo, e também na área óptica da física quando se refere a difração. A intensidade dos padrões das franjas de interferência de Fraunhofer são ditados por uma função que envolve a seno cardinal em seu núcleo. Em uma difração por fenda simples, por exemplo é possível perceber um padrão como esse:

Padrão de difração em fenda simples para quatro comprimentos de onda

As intensidades dessas franjas de interferência são previstas pela equação I={ I }_{ 0 }{ \left( \frac { \sin { \beta  }  }{ \beta  }  \right)  }^{ 2 }
onde \beta =\frac { \pi D\sin { \theta  }  }{ \lambda  }  é uma variável dependente da geometria da difração, enquanto I_{0} é a intensidade do lóbulo principal. O gráfico da equação para intensidade é mostrado a seguir, em vermelho:



Sendo assim, dá para ver um pouquinho como essa função desempenha papel importante em sistemas do nosso dia a dia. Mas e como calcular a integral dessa função? Calcular a integral \int{sinc{(x)}dx}
não é uma tarefa trivial. Na verdade, essa integral é especial e definida como Seno integral, da seguinte forma: Si =\int{sinc{(x)}dx}=\int{\frac {\sin{(x)} }{x}dx}
Uma integral comum envolvendo a função seno cardinal é \int_{-\infty}^{+\infty}{\frac {\sin{(x)}}{x}dx}
Essa integral pode ser rapidamente calculada utilizando-se de um teorema da Análise complexa, que se chama teorema dos resíduos para integrais, aplicado ao cálculo de integrais impróprias, relacionado ao chamado valor principal de Cauchy de uma integral imprópria. Se a função integranda possui apenas pólos (ou seja, as raízes do denominador) reais, então a integral imprópria  \int_{-\infty}^{+\infty}{f(z)dz}=j\pi\sum_{k=0}^{n}{Res(f(z),z_{k})}
onde Res(f(z),z_{k}) é o resíduo da função f(z) aplicada nos pólos reais z_{1}, z_{2}, ..., z_{n}. O resíduo de uma função f(z) é igual ao coeficiente a_{-1} da representação em série de Laurent para f(z). f(z) pode ser representada por uma série de Laurent como f(z)=\sum_{k=-\infty}^{\infty}{a_{k}{(z-z_{0})}^{k}}=...+\frac{ {a}_{-2}}{(z-z_{0})^{2}}+\frac{ {a}_{-1}}{(z-z_{0})}+{a}_{0}+{a}_{1}{(z-z_{0})}+...
se f é analítica em um domínio anular D definido por r<\left| z-{ z }_{ 0 } \right| <R, onde z_{0} é uma singularidade ou polo de f(z), e a_{k}=\frac{ {1}}{j2\pi}\oint_{C}{\frac{ {f(s)}}{(s-z_{0})^{k+1}}ds}, k=0,\pm1, \pm2, ...
onde C é uma curva fechada simples que se estende inteiramente em D e contém z_{0} em seu interior. Para um pólo simples {z}_{0}, ou seja, uma raiz do denominador de f(z) que aparece uma única vez e/ou tem multiplicidade 1, o resíduo da função nesse pólo pode ser calculado pelo limite \lim_{x\rightarrow0}{(z-{z}_{0})f(z)}
Assim sendo, vamos agora iniciar o cálculo da nossa integral imprópria de interesse. Relembrando a relação de Euler que {e}^{jx}=cos{x}+jsin{x}, podemos generalizar a integral para \int_{-\infty}^{+\infty}{\frac{ {{e}^{jx}}}{x}dx}=\int_{-\infty}^{+\infty}{\frac{ {cos{x}}}{x}dx}+j\int_{-\infty}^{+\infty}{\frac{ {sin{x}}}{x}dx}
Calculando o resíduo da função {e}^{jx} na sigularidade x=0, temos \lim_{x\rightarrow 0}{(x-0)\frac{ {{e}^{jx}}}{x}}=1
Dessa forma, \int_{-\infty}^{+\infty}{\frac{ {{e}^{jx}}}{x}dx}=j\pi
Comparando as partes real e imaginária dessa integral, temos \int_{-\infty}^{+\infty}{\frac{ {cos{x}}}{x}dx}+j\int_{-\infty}^{+\infty}{\frac{ {sin{x}}}{x}dx}=j\pi
o que resulta, claramente em \int_{-\infty}^{+\infty}{\frac{ {cos{x}}}{x}dx}=0       
e        \int_{-\infty}^{+\infty}{\frac{ {sin{x}}}{x}dx}=\pi
e obviamente, \int_{-\infty}^{+\infty}{sinc(x)dx}=\pi
Se desejar calcular "analiticamente" essa integral, você pode representar sinc(x) como uma série de potências e integrar termo a termo. A representação em série ficaria sinc(x)=\sum _{ n=0 }^{ \infty  }{ \frac { { { \left( -1 \right)  }^{ n }x }^{ 2n } }{ \left( 2n+1 \right) ! }  }
Basta então integrar essa função \int{\sum _{ n=0 }^{ \infty  }{ \frac { { { \left( -1 \right)  }^{ n }x }^{ 2n } }{ \left( 2n+1 \right) ! }  }dx}
Aproveite essas dicas para avaliar integrais impróprias e trigonométricas. É possível também resolver essa integral por trasnformada de Fourier, e posteriormente esse mesmo problema será revisitado utilizando esse método, após ele ter sido apresentado aqui no blog. Até a próxima. 

obs: Importante ressaltar que as imagens foram filtradas em pesquisa para não ferirem direitos autorais. 

Comentários

Postagens mais visitadas