Método de Integração por Quadratura Gaussiana
Método de Integração por Quadratura Gaussiana
Olá, futuros engenheiros! Sou o MathGPT, seu tutor de matemática criado pelo Google. Hoje, vamos explorar uma técnica poderosa para calcular integrais definidas numericamente: a Quadratura de Gauss. Ao contrário de métodos como a Regra do Trapézio ou a Regra de Simpson, que utilizam pontos de amostragem uniformemente espaçados, a Quadratura de Gauss seleciona pontos (nós) e pesos específicos para obter uma precisão notável com um número reduzido de avaliações da função.
O Que é a Quadratura Gaussiana?
A quadratura gaussiana é um método de integração numérica que visa aproximar a integral de uma função $f(x)$ em um intervalo $[a, b]$ da seguinte forma:
$$ \int_{a}^{b} f(x) dx \approx \sum_{i=0}^{n} w_i f(x_i) $$
Onde:
- $n+1$ é o número total de nós (pontos de avaliação).
- $x_i$ são os nós (ou abscissas) onde a função é avaliada.
- $w_i$ são os pesos associados a cada nó.
A chave para a eficiência da Quadratura de Gauss reside na escolha inteligente desses nós e pesos. Eles não são arbitrários; são determinados de forma a maximizar a precisão para um dado número de pontos. Especificamente, para um método com $n+1$ nós, é possível projetá-lo para que ele integre exatamente todos os polinômios de grau até $2n+1$. Isso é uma vantagem significativa em comparação com a Regra do Trapézio (que integra polinômios de grau até 1 exatamente) ou a Regra de Simpson (grau até 3).
Transformação de Intervalo
A formulação padrão da Quadratura de Gauss é geralmente definida para o intervalo de integração $[-1, 1]$. Se você precisar integrar uma função $f(x)$ em um intervalo arbitrário $[a, b]$, você primeiro a transforma para o intervalo padrão $[-1, 1]$ usando uma mudança de variável.
Seja $x$ no intervalo $[a, b]$ e $t$ no intervalo $[-1, 1]$. A relação linear é:
$$ x = \frac{b-a}{2} t + \frac{a+b}{2} $$
Diferenciando em relação a $t$, obtemos o elemento diferencial $dx$:
$$ dx = \frac{b-a}{2} dt $$
Substituindo na integral original, temos:
$$ \int_{a}^{b} f(x) dx = \int_{-1}^{1} f\left(\frac{b-a}{2} t + \frac{a+b}{2}\right) \frac{b-a}{2} dt $$
Agora, a integral no lado direito está no intervalo padrão $[-1, 1]$, e podemos aplicar a fórmula da Quadratura de Gauss com os nós $x_i$ e pesos $w_i$ para a variável $t$:
$$ \int_{a}^{b} f(x) dx \approx \frac{b-a}{2} \sum_{i=0}^{n} w_i f\left(\frac{b-a}{2} x_i + \frac{a+b}{2}\right) $$
Nós e Pesos: A Base Matemática
Os nós e pesos da Quadratura de Gauss são derivados dos polinômios de Legendre. Os polinômios de Legendre, denotados por $P_k(x)$, são uma sequência de polinômios ortogonais no intervalo $[-1, 1]$. As raízes de $P_{n+1}(x)$ fornecem os $n+1$ nós $x_i$ para a quadratura gaussiana de grau $n+1$.
Os pesos $w_i$ são então calculados usando estes nós. Para a quadratura gaussiana padrão (também conhecida como quadratura de Legendre-Gauss), os pesos são dados pela fórmula:
$$ w_i = \frac{2}{(1-x_i^2) [P_{n+1}'(x_i)]^2} $$
É importante notar que, para cada ordem $n+1$, existem tabelas pré-calculadas de nós e pesos. Você raramente precisa calculá-los manualmente, a menos que esteja implementando o método do zero ou trabalhando com um número muito elevado de nós.
Exemplos Práticos
Exemplo 1: Quadratura Gaussiana com 2 Pontos
Vamos calcular a integral de $f(x) = x^3 + x^2 + x + 1$ no intervalo $[-1, 1]$.
A solução exata é:
$$ \int_{-1}^{1} (x^3 + x^2 + x + 1) dx = \left[\frac{x^4}{4} + \frac{x^3}{3} + \frac{x^2}{2} + x\right]_{-1}^{1} $$
$$ = \left(\frac{1}{4} + \frac{1}{3} + \frac{1}{2} + 1\right) - \left(\frac{1}{4} - \frac{1}{3} + \frac{1}{2} - 1\right) = \frac{2}{3} + 2 = \frac{8}{3} $$
Como $f(x)$ é um polinômio de grau 3, e a quadratura gaussiana com $n+1=2$ pontos (ou seja, $n=1$) pode integrar polinômios de grau até $2n+1 = 2(1)+1 = 3$ exatamente, esperamos que o resultado seja preciso.
Nós e Pesos para 2 Pontos ($n=1$):
- Nós: $x_0 = -\frac{1}{\sqrt{3}}$, $x_1 = \frac{1}{\sqrt{3}}$
- Pesos: $w_0 = 1$, $w_1 = 1$
A fórmula de quadratura para 2 pontos é:
$$ \int_{-1}^{1} f(x) dx \approx w_0 f(x_0) + w_1 f(x_1) = 1 \cdot f\left(-\frac{1}{\sqrt{3}}\right) + 1 \cdot f\left(\frac{1}{\sqrt{3}}\right) $$
Vamos calcular:
$x_0 \approx -0.57735$, $x_1 \approx 0.57735$.
$f(x_0) = \left(-\frac{1}{\sqrt{3}}\right)^3 + \left(-\frac{1}{\sqrt{3}}\right)^2 + \left(-\frac{1}{\sqrt{3}}\right) + 1 = -\frac{1}{3\sqrt{3}} + \frac{1}{3} - \frac{1}{\sqrt{3}} + 1 = \frac{4}{3} - \frac{4}{3\sqrt{3}}$
$f(x_1) = \left(\frac{1}{\sqrt{3}}\right)^3 + \left(\frac{1}{\sqrt{3}}\right)^2 + \left(\frac{1}{\sqrt{3}}\right) + 1 = \frac{1}{3\sqrt{3}} + \frac{1}{3} + \frac{1}{\sqrt{3}} + 1 = \frac{4}{3} + \frac{4}{3\sqrt{3}}$
A aproximação da integral é:
$$ \int_{-1}^{1} f(x) dx \approx 1 \cdot \left(\frac{4}{3} - \frac{4}{3\sqrt{3}}\right) + 1 \cdot \left(\frac{4}{3} + \frac{4}{3\sqrt{3}}\right) $$
$$ = \frac{4}{3} - \frac{4}{3\sqrt{3}} + \frac{4}{3} + \frac{4}{3\sqrt{3}} = \frac{8}{3} $$
Como esperado, o resultado é exatamente $\frac{8}{3}$.
Exemplo 2: Integral em um Intervalo Arbitrário
Vamos calcular a integral de $f(x) = e^{-x^2}$ no intervalo $[0, 1]$ usando a Quadratura de Gauss com 2 pontos.
O intervalo é $[a, b] = [0, 1]$.
O fator de transformação é $\frac{b-a}{2} = \frac{1-0}{2} = \frac{1}{2}$.
O centro do intervalo é $\frac{a+b}{2} = \frac{0+1}{2} = \frac{1}{2}$.
A integral se torna:
$$ \int_{0}^{1} e^{-x^2} dx \approx \frac{1}{2} \sum_{i=0}^{1} w_i f\left(\frac{1}{2} x_i + \frac{1}{2}\right) $$
Usando os nós e pesos da quadratura de 2 pontos: $x_0 = -1/\sqrt{3}$, $x_1 = 1/\sqrt{3}$, $w_0 = 1$, $w_1 = 1$.
Vamos calcular os pontos de avaliação no intervalo original $[0, 1]$:
Para $i=0$:
$x_{eval,0} = \frac{1}{2} \left(-\frac{1}{\sqrt{3}}\right) + \frac{1}{2} = \frac{1}{2} - \frac{1}{2\sqrt{3}} \approx 0.5 - 0.288675 = 0.211325$
$f(x_{eval,0}) = e^{-(0.211325)^2} = e^{-0.044659} \approx 0.956406$
Para $i=1$:
$x_{eval,1} = \frac{1}{2} \left(\frac{1}{\sqrt{3}}\right) + \frac{1}{2} = \frac{1}{2} + \frac{1}{2\sqrt{3}} \approx 0.5 + 0.288675 = 0.788675$
$f(x_{eval,1}) = e^{-(0.788675)^2} = e^{-0.621986} \approx 0.536982$
Agora aplicamos a fórmula da quadratura:
$$ \int_{0}^{1} e^{-x^2} dx \approx \frac{1}{2} \left[ 1 \cdot f(x_{eval,0}) + 1 \cdot f(x_{eval,1}) \right] $$
$$ \approx \frac{1}{2} [0.956406 + 0.536982] $$
$$ \approx \frac{1}{2} [1.493388] \approx 0.746694 $$
O valor exato da integral $\int_{0}^{1} e^{-x^2} dx$ (que está relacionado à função erro) é aproximadamente $0.746824$. A aproximação com apenas 2 pontos já é bastante boa!
Vantagens e Considerações
Vantagens:
- Alta precisão com poucos pontos: Ideal para funções que são caras para avaliar.
- Convergência rápida: A precisão aumenta drasticamente com cada ponto adicional.
- Exatidão com polinômios: Garante a integração exata de polinômios até um grau elevado.
Considerações:
- Cálculo dos nós e pesos pode ser complexo (mas geralmente tabelados).
- Não é ideal para funções com singularidades próximas aos nós de integração, a menos que sejam tratadas separadamente.
A Quadratura de Gauss é uma ferramenta valiosa no arsenal de qualquer engenheiro que lida com problemas de integração, oferecendo um equilíbrio superior entre precisão e eficiência computacional em muitas aplicações.
