Matemática Numérica III

5 Integração 5.1 Integração Autoadaptativa Referências

Ajude a manter o site livre, gratuito e sem propagandas. Colabore!

5.2 Integrais múltiplas

Vamos trabalhar com métodos para a computação de integrais múltiplas

\int\int_{R}f(x,y)\,dA.

(5.13)

Em uma região retangular $A=[a,b]\times[c,d]$ , podemos reescrevê-la como uma integral iterada

\int\int_{R}f(x,y)\,dA=\int_{a}^{b}\int_{c}^{d}f(x,y)\,dy\,dx.

(5.14)

5.2.1 Regras de Newton-Cotes

Regra do Trapézio

A Regra do Trapézio²²2Notas de Aula - Matemática Numérica. nos fornece

\int_{c}^{d}f(x,y)\,dy=\frac{h_{y}}{2}\left[f(x,c)+f(x,d)\right]-\frac{h_{y}^{% 3}}{12}f^{\prime\prime}(x,\eta)

(5.15)

com $h_{y}=(d-c)$ e $\eta\in(c,d)$ . De forma iterada, temos

	$\displaystyle\int_{a}^{b}\int_{c}^{d}f(x,y)\,dy\,dx$	$\displaystyle=\frac{h_{y}}{2}\int_{a}^{b}f(x,c)\,dx+\frac{h_{y}}{2}\int_{a}^{b% }f(x,d)\,dx$		(5.16)
		$\displaystyle-\frac{h_{y}^{3}}{12}\int_{a}^{b}f^{\prime\prime}(x,\eta)\,dx.$		(5.17)

Então, à exceção do termo do erro, aplicamos a Regra do Trapézio para as integrais em $x$ . Obtemos

$\displaystyle\int_{a}^{b}\int_{c}^{d}f(x,y)\,dy\,dx$	$\displaystyle=\frac{h_{y}}{2}\frac{h_{x}}{2}\left[f(a,c)+f(b,c)\right]$	(5.18)
	$\displaystyle+\frac{h_{y}}{2}\frac{h_{x}}{2}\left[f(a,d)+f(b,d)\right]$	(5.19)
	$\displaystyle-\frac{h_{y}}{2}\frac{h_{x}^{3}}{12}f^{\prime\prime}(\mu^{\prime}% ,c)$	(5.20)
	$\displaystyle-\frac{h_{y}}{2}\frac{h_{x}^{3}}{12}f^{\prime\prime}(\mu^{\prime% \prime},d)$	(5.21)
	$\displaystyle-\frac{h_{y}^{3}}{12}\int_{a}^{b}f^{\prime\prime}(x,\eta)\,dx,$	(5.22)

com $h_{x}=(b-a)$ , $\mu^{\prime},\mu^{\prime\prime}\in(a,b)$ . Pelos Teorema do Valor Intermediário e pelo Teorema do Valor Médio, podemos ver que o erro é $O(h_{x}h_{y}^{3}+h_{x}^{3}h_{y})$ . Por fim, obtemos a Regra do Trapézio para Integrais Iteradas

	$\displaystyle\int_{a}^{b}\int_{c}^{d}f(x,y)\,dy\,dx$	$\displaystyle=\frac{h_{y}}{2}\frac{h_{x}}{2}\left[f(a,c)+f(b,c)+f(b,d)+f(a,d)\right]$		(5.23)
		$\displaystyle+O(h_{x}h_{y}^{3}+h_{x}^{3}h_{y}).$		(5.24)

Exemplo 5.2.1.

A Regra do Trapézio fornece

\int_{1.5}^{2}\int_{1}^{1.5}\ln(x+2y)\,dy\,dx\approx 0.36.

(5.25)

Verifique!

Regra de Simpson

A Regra do Simpson³³3Notas de Aula - Matemática Numérica. nos fornece

	$\displaystyle\int_{c}^{d}f(x,y)\,dy$	$\displaystyle=\frac{h_{y}}{3}\left[f(x,y_{1})+4f(x,y_{2})+f(x,y_{3})\right]$		(5.26)
		$\displaystyle-\frac{h_{y}^{5}}{90}f^{(4)}(x,\eta)$		(5.27)

com $h_{y}=(d-c)/2$ , $y_{j}=(j-1)h_{y}$ , $j=1,2,3$ , e $\eta\in(c,d)$ . De forma iterada, temos

	$\displaystyle\int_{a}^{b}\int_{c}^{d}f(x,y)\,dy\,dx$	$\displaystyle=\frac{h_{y}}{3}\left[\int_{a}^{b}f(x,y_{1})\,dx+4\int_{a}^{b}f(x% ,y_{2})\,dx+\int_{a}^{b}f(x,y_{3})\,dx\right]$		(5.28)
		$\displaystyle-\frac{h_{y}^{5}}{90}\int_{a}^{b}f^{(4)}(x,\eta)\,dx$		(5.29)

Então, à exceção do termo do erro, aplicamos a Regra de Simpson para as integrais em $x$ . Obtemos

$\displaystyle\int_{a}^{b}\int_{c}^{d}f(x,y)\,dy\,dx$	$\displaystyle=\frac{h_{x}h_{y}}{9}\left[f(x_{1},y_{1})+4f(x_{2},y_{1})+f(x_{3}% ,y_{1})\right]$	(5.30)
	$\displaystyle+\frac{4h_{x}h_{y}}{9}\left[f(x_{1},y_{2})+4f(x_{2},y_{2})+f(x_{3% },y_{2})\right]$	(5.31)
	$\displaystyle+\frac{h_{x}h_{y}}{9}\left[f(x_{1},y_{3})+4f(x_{2},y_{3})+f(x_{3}% ,y_{3})\right]$	(5.32)
	$\displaystyle-\frac{h_{x}^{5}h_{y}}{270}f^{(4)}(\mu_{1},y_{1})$	(5.33)
	$\displaystyle-\frac{4h_{x}^{5}h_{y}}{270}f^{(4)}(\mu_{2},y_{2})$	(5.34)
	$\displaystyle-\frac{h_{x}^{5}h_{y}}{270}f^{(4)}(\mu_{3},y_{3})$	(5.35)
	$\displaystyle-\frac{h_{y}^{5}}{90}\int_{a}^{b}f^{(4)}(x,\eta)\,dx$	(5.36)

com $h_{x}=(b-a)/2$ , $\mu_{1},\mu_{2},\mu_{3}\in(a,b)$ . Pelos Teorema do Valor Intermediário e Teorema do Valor Médio, podemos ver que o erro é $O(h_{x}h_{y}^{5}+h_{x}^{5}h_{y})$ . Por fim, obtemos a Regra de Simpson para Integrais Iteradas

$\displaystyle\int_{a}^{b}\int_{c}^{d}f(x,y)\,dy\,dx$	$\displaystyle=\frac{h_{x}h_{y}}{9}\left[f(x_{1},y_{1})+4f(x_{2},y_{1})+f(x_{3}% ,y_{1})\right]$	(5.37)
	$\displaystyle+\frac{4h_{x}h_{y}}{9}\left[f(x_{1},y_{2})+4f(x_{2},y_{2})+f(x_{3% },y_{2})\right]$	(5.38)
	$\displaystyle+\frac{h_{x}h_{y}}{9}\left[f(x_{1},y_{3})+4f(x_{2},y_{3})+f(x_{3}% ,y_{3})\right]$	(5.39)
	$\displaystyle+O(h_{x}h_{y}^{5}+h_{x}^{5}h_{y}).$	(5.40)

Exemplo 5.2.2.

A Regra de Simpson fornece

\int_{1.5}^{2}\int_{1}^{1.5}\ln(x+2y)\,dy\,dx\approx 0.361003.

(5.41)

Verifique!

5.2.2 Regras Compostas de Newton-Cotes

A ideia é particionar a região de integração em células e o resultado da integração é a soma da aplicação da regra de quadratura em cada uma das células.

Regra Composta do Trapézio

Para uma região retangular $R=[a,b]\times[c,d]$ , vamos construir a malha

M=\{c_{k}=[x_{i},x_{i+1}]\times[y_{j},y_{j+1}]:\leavevmode\nobreak\ k=i+(j-1)n% _{x}\},

(5.42)

onde $x_{i}=(i-1)h_{x}$ , $h_{x}=(b-a)/n_{x}$ , $i=1,2,\dotsc,n_{x}+1$ e $y_{j}=(j-1)h_{y}$ , $h_{y}=(d-c)/n_{y}$ , $j=1,2,\dotsc,n_{y}+1$ . Consulte a Figura 5.1.

Refer to caption — Figura 5.1: Representação da malha para a Regra Composta do Trapézio.

Aplicando a ideia, temos

	$\displaystyle\int_{a}^{b}\int_{c}^{d}f(x,y)\,dy\,dx$	$\displaystyle=\sum_{k=1}^{n_{x}n_{y}}\int\int_{C_{k}}f(x,y)\,dy\,dx$		(5.43)
		$\displaystyle=\sum_{i=1}^{n_{x}}\sum_{j=1}^{n_{y}}\int_{x_{i}}^{x_{i+1}}\int_{% y_{j}}^{y_{j+1}}f(x,y)\,dy\,dx$		(5.44)

Em cada integral em $C_{k}$ , aplicamos a Regra do Trapézio, segue

	$\displaystyle\int_{x_{i}}^{x_{i+1}}\int_{y_{j}}^{y_{j+1}}f(x,y)\,dy\,dx$	$\displaystyle\approx\frac{h_{y}}{2}\frac{h_{x}}{2}\left[f(x_{i},y_{j})+f(x_{i+% 1},y_{j})\right.$		(5.45)
		$\displaystyle+\left.f(x_{i+1},y_{j+1})+f(x_{i+1},y_{j})\right]$		(5.46)

Observamos que nos conjuntos de nodos (marcados em azul na Figura 5.1)

	$\displaystyle\{(i,j):\leavevmode\nobreak\ i=2,\dotsc,n_{x},j=1\text{ ou }j=n_{% y}+1\},$		(5.47)
	$\displaystyle\{(i,j):\leavevmode\nobreak\ i=1\text{ ou }i=n_{x}+1,j=2,\dotsc,n% _{y}\}$		(5.48)

a função integranda será avaliada $2$ vezes. Já, em todos os nodos internos, $i=2,\dotsc,n_{x}$ , $j=2,\dotsc,n_{y}$ , a função será avaliada $4$ vezes. Com isso, chegamos à Regra Composta do Trapézio

$\displaystyle\int_{a}^{b}\int_{c}^{d}f(x,y)\,dy\,dx$	$\displaystyle=\frac{h_{x}h_{y}}{4}\left[f(x_{1},y_{1})+f(x_{n_{x}+1},y_{1})\right.$
	$\displaystyle\left.+f(x_{n_{x}+1},y_{n_{y}+1})+f(x_{1},y_{n_{y}+1})\right]$
	$\displaystyle+\frac{h_{x}h_{y}}{2}\sum_{i=2}^{n_{x}}\left[f(x_{i},y_{1})+f(x_{% i},y_{n_{y}+1})\right]$
	$\displaystyle+\frac{h_{x}h_{y}}{2}\sum_{j=2}^{n_{y}}\left[f(x_{1},y_{j})+f(x_{% n_{x}+1},y_{j})\right]$
	$\displaystyle+h_{x}h_{y}\sum_{i=2}^{n_{x}}\sum_{j=2}^{n_{y}}f(x_{i},y_{j})$
	$\displaystyle+O(h_{x}^{2}+h_{y}^{2})$	(5.49)

Observamos que esta é uma quadratura de $(n_{x}+1)(n_{y}+1)$ nodos.

E. 5.2.1.

Verifique a aplicação da Regra Composta do Trapézio para computar

\int_{1.5}^{2}\int_{1}^{1.5}\ln(x+2y)\,dy\,dx.

(5.50)

Regra Composta de Simpson

Aqui, vamos construir uma malha

M=\left\{C_{k}=[x_{i},x_{i+2}]\times[y_{j},y_{j+2}]:\leavevmode\nobreak\ i=1,3% ,\dotsc,n_{x}-1,j=1,3,\dotsc,n_{y}-1,\leavevmode\nobreak\ \right\},

(5.51)

onde $x_{i}=(i-1)h_{x}$ , $h_{x}=(b-a)/n_{x}$ , $i=1,2,\dotsc,n_{x}+1$ e $y_{j}=(j-1)h_{y}$ , $h_{y}=(d-c)/n_{y}$ , $j=1,2,\dotsc,n_{y}+1$ . Com $n_{x},n_{y}\leq 2$ números pares. Consulte a Figura 5.2.

A Regra Composta de Simpson para Integrais Iteradas fica

$\displaystyle\int_{a}^{b}\int_{c}^{d}f(x,y)\,dx\,dy$	$\displaystyle=\frac{h_{x}h_{y}}{9}\left\{f(x_{1},y_{1})+f(x_{n_{x}+1},y_{1})\right.$
	$\displaystyle+f(x_{1},y_{n_{y}+1})+f(x_{n_{x}+1},y_{n_{y}+1})$
	$\displaystyle+2\sum_{i=1}^{n_{x}/2-1}\left[f(x_{2i+1},y_{1})+f(x_{2i+1},y_{n_{% y}+1})\right]$
	$\displaystyle+2\sum_{j=1}^{n_{y}/2-1}\left[f(x_{1},y_{2j+1})+f(x_{n_{x}+1},y_{% 2j+1})\right]$
	$\displaystyle+4\sum_{i=1}^{n_{x}/2}\left[f(x_{2i},y_{1})+f(x_{2i},y_{n_{y}+1})\right]$
	$\displaystyle+4\sum_{j=1}^{n_{y}/2}\left[f(x_{1},y_{2j})+f(x_{n_{x}+1},y_{2j})\right]$
	$\displaystyle+4\sum_{i=1}^{n_{x}/2-1}\sum_{j=1}^{n_{y}/2-1}f(x_{2j+1},y_{2j+1})$
	$\displaystyle+8\sum_{i=1}^{n_{x}/2-1}\sum_{j=1}^{n_{y}/2}f(x_{2i+1},y_{2j})$
	$\displaystyle+8\sum_{i=1}^{n_{x}/2}\sum_{j=1}^{n_{y}/2-1}f(x_{2i},y_{2j+1})$
	$\displaystyle\left.+16\sum_{i=1}^{n_{x}/2}\sum_{j=1}^{n_{y}/2}f(x_{2i},y_{2j})\right\}$
	$\displaystyle+O(h_{x}^{4}+h_{y}^{4}).$	(5.52)

E. 5.2.2.

Verifique a aplicação da Regra Composta de Simpson para computar

\int_{1.5}^{2}\int_{1}^{1.5}\ln(x+2y)\,dy\,dx.

(5.53)

Envie seu comentário

Aproveito para agradecer a todas/os que de forma assídua ou esporádica contribuem enviando correções, sugestões e críticas!

Este texto é disponibilizado nos termos da Licença Creative Commons Atribuição-CompartilhaIgual 4.0 Internacional. Ícones e elementos gráficos podem estar sujeitos a condições adicionais.

Política de uso de dados