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⁷⁰⁷⁰endnote: ⁷⁰Notas 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⁷¹⁷¹endnote: ⁷¹Notas 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)

Notas

1 Siméon Denis Poisson, 1781 - 1840, matemático francês. Fonte: Wikipédia:Siméon Denis Poisson.
2 Siméon Denis Poisson, 1781 - 1840, matemático francês. Fonte: Wikipédia:Siméon Denis Poisson.
3 Pierre-Simon Laplace, 1749 - 1827, matemático francês. Fonte: Wikipédia: Pierre-Simon Laplace.
4 Siméon Denis Poisson, 1781 - 1840, matemático francês. Fonte: Wikipédia:Siméon Denis Poisson.
5 Siméon Denis Poisson, 1781 - 1840, matemático francês. Fonte: Wikipédia:Siméon Denis Poisson.
6 Siméon Denis Poisson, 1781 - 1840, matemático francês. Fonte: Wikipédia:Siméon Denis Poisson.
7 Siméon Denis Poisson, 1781 - 1840, matemático francês. Fonte: Wikipédia:Siméon Denis Poisson.
8 Siméon Denis Poisson, 1781 - 1840, matemático francês. Fonte: Wikipédia:Siméon Denis Poisson.
9 Siméon Denis Poisson, 1781 - 1840, matemático francês. Fonte: Wikipédia:Siméon Denis Poisson.
10 CSR é mais eficiente em muitos casos.
11 Siméon Denis Poisson, 1781 - 1840, matemático francês. Fonte: Wikipédia:Siméon Denis Poisson.
12 Siméon Denis Poisson, 1781 - 1840, matemático francês. Fonte: Wikipédia:Siméon Denis Poisson.
13 Siméon Denis Poisson, 1781 - 1840, matemático francês. Fonte: Wikipédia:Siméon Denis Poisson.
14 Semyon Aronovich Gershgorin, 1901 - 1933, matemático russo-soviético. Fonte: Wikipedia.
15 Siméon Denis Poisson, 1781 - 1840, matemático francês. Fonte: Wikipédia:Siméon Denis Poisson.
16 Siméon Denis Poisson, 1781 - 1840, matemático francês. Fonte: Wikipédia:Siméon Denis Poisson.
17 Siméon Denis Poisson, 1781 - 1840, matemático francês. Fonte: Wikipédia:Siméon Denis Poisson.
18 Siméon Denis Poisson, 1781 - 1840, matemático francês. Fonte: Wikipédia:Siméon Denis Poisson.
19 Siméon Denis Poisson, 1781 - 1840, matemático francês. Fonte: Wikipédia:Siméon Denis Poisson.
20 Magnus R. Hestenes, 1906 - 1991, matemático americano estadunidense. Fonte: Wikipedia.
21 Eduard Stiefel, 1909 - 1978, matemático suíço. Fonte: Wikipédia.
22 Siméon Denis Poisson, 1781 - 1840, matemático francês. Fonte: Wikipédia:Siméon Denis Poisson.
23 Georgi Iwanowitsch Petrov, 1912 - 1987, engenheiro soviético. Fonte: Wikipedia.
24 Boris Galerkin, 1871 - 1945, engenheiro e matemático soviético. Fonte: Wikipédia.
25 Desenvolvido por Yousef Saad e H. Schultz, 1986. Fonte: Wikipedia.
26 Alexei Nikolajewitsch Krylov, 1863 - 1945, engenheiro e matemático russo. Fonte: Wikipédia.
27 Georgi Iwanowitsch Petrov, 1912 - 1987, engenheiro soviético. Fonte: Wikipedia.
28 Boris Galerkin, 1871 - 1945, engenheiro e matemático soviético. Fonte: Wikipédia.
29 Jørgen Pedersen Gram, 1850 - 1916, matemático dinamarquês. Fonte: Wikipédia.
30 Erhard Schmidt, 1876 - 1959, matemático alemão. Fonte: Wikipédia.
31 Walter Edwin Arnoldi, 1917 - 1995, engenheiro americano estadunidense. Fonte: Wikipédia.
32 Karl Adolf Hessenberg, 1904 - 1959, engenheiro e matemático alemão. Fonte: Wikipédia.
33 Alston Scott Householder, 1904 - 1993, matemático americano estadunidense. Fonte: Wikipédia.
34 Carl Gustav Jakob Jacobi, 1804 - 1851, matemático alemão. Fonte: Wikipédia: Carl Gustav Jakob Jacobi.
35 Isaac Newton, 1642 - 1727, matemático, físico, astrônomo, teólogo e autor inglês. Fonte: Wikipédia: Isaac Newton.
36 Rudolf Otto Sigismund Lipschitz, 1832 - 1903, matemático alemão. Fonte: Wikipédia.
37 Jan Burgers, 1895 - 1981, físico neerlandês. Fonte: Wikipédia: Jan Burgers.
38 Johann Peter Gustav Lejeune Dirichlet, 1805 - 1859, matemático alemão. Fonte: Wikipédia: Johann Peter Gustav Lejeune Dirichlet.
39 Erich Hans Rothe, 1895 - 1988, matemático alemão. Fonte: Wikipédia.
40 Leonhard Paul Euler, 1707-1783, matemático e físico suíço. Fonte: Wikipédia: Ronald Fisher.
41 Isaac Newton, 1642 - 1727, matemático, físico, astrônomo, teólogo e autor inglês. Fonte: Wikipédia: Isaac Newton.
42 Carl Gustav Jakob Jacobi, 1804 - 1851, matemático alemão. Fonte: Wikipédia: Carl Gustav Jakob Jacobi.
43 Carl Gustav Jakob Jacobi, 1804 - 1851, matemático alemão. Fonte: Wikipédia: Carl Gustav Jakob Jacobi.
44 Isaac Newton, 1642 - 1727, matemático, físico, astrônomo, teólogo e autor inglês. Fonte: Wikipédia: Isaac Newton.
45 Rudolf Otto Sigismund Lipschitz, 1832 - 1903, matemático alemão. Fonte: Wikipédia.
46 Jan Burgers, 1895 - 1981, físico neerlandês. Fonte: Wikipédia: Jan Burgers.
47 Johann Peter Gustav Lejeune Dirichlet, 1805 - 1859, matemático alemão. Fonte: Wikipédia: Johann Peter Gustav Lejeune Dirichlet.
48 Erich Hans Rothe, 1895 - 1988, matemático alemão. Fonte: Wikipédia.
49 Leonhard Paul Euler, 1707-1783, matemático e físico suíço. Fonte: Wikipédia: Ronald Fisher.
50 Isaac Newton, 1642 - 1727, matemático, físico, astrônomo, teólogo e autor inglês. Fonte: Wikipédia: Isaac Newton.
51 Carl Gustav Jakob Jacobi, 1804 - 1851, matemático alemão. Fonte: Wikipédia: Carl Gustav Jakob Jacobi.
52 Brook Taylor, 1685 - 1731, matemático britânico. Fonte: Wikipédia:Brook Taylor.
53 Philip Wolfe, 1927 - 2016, matemático estadunidense. Fonte: Wikipédia.
54 Howard Harry Rosenbrock, 1920 - 2010, engenheiro britânico. Fonte: Wikipedia: Howard Harry Rosenbrock.
55 Leopold Kronecker, 1923 - 1891, matemático alemão. Fonte: Wikipédia: Leopold Kronecker.
56 Howard Harry Rosenbrock, 1920 - 2010, engenheiro britânico. Fonte: Wikipedia: Howard Harry Rosenbrock.
57 Isaac Newton, 1642 - 1727, matemático, físico, astrônomo, teólogo e autor inglês. Fonte: Wikipédia: Isaac Newton.
58 Alexei Nikolajewitsch Krylov, 1863 - 1945, engenheiro e matemático russo. Fonte: Wikipédia.
59 Howard Harry Rosenbrock, 1920 - 2010, engenheiro britânico. Fonte: Wikipedia: Howard Harry Rosenbrock.
60 Howard Harry Rosenbrock, 1920 - 2010, engenheiro britânico. Fonte: Wikipedia: Howard Harry Rosenbrock.
61 Uma matriz $A\in\mathbb{C}^{n\times n}$ é diagonalizável, quando tem $n$ autovetores linearmente independentes.
62 John William Strutt, 3º Barão Rayleigh, 1842 - 1919, físico e matemático britânico. Fonte: Wikipédia: John William Strutt.
63 A afirmação segue por indução matemática.
64 Issai Schur, 1875 - 1941, matemático russo-alemão. Fonte: Wikipédia.
65 Uma matriz $U$ é dita unitária quando $U^{-1}=U^{H}$ .
66 Uma matriz $Q$ é dita ortogonal quando $Q^{T}Q=I$ .
67 Karl Adolf Hessenberg, 1904 - 1959, engenheiro e matemático alemão. Fonte: Wikipédia.
68 James Wallace Givens Jr., 1910 - 1993, matemático estadunidense. Fonte: Wikipédia.
69 Consulte mais sobre a Regra de Simpson em Seção 10.1 Regras de Newton-Cotes.
70 Notas de Aula - Matemática Numérica.
71 Notas de Aula - Matemática Numérica.

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