Cálculo Numérico: 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.