Apostila - Métodos Numéricos
Apostila - Métodos Numéricos
Apostila - Métodos Numéricos
1.1 Base 10
1.2 Base 2
Um mesmo número pode ser representado em diferentes bases, vamos nesse curso apenas
trabalhar com a base 2 (binária) e base 10.
Exemplo 1.3.1
3
Capítulo 1. Conversão de Base Numérica (Unidade 01)
Erros (Unidade 02 e 03)
Exemplo 1.3.2
13 2
1 6 2
Resolução: Então, (13)10 = (1101)2
0 3 2
1 1
1.4 Exercícios
1. Transformar os números: x1 = 34, x2 = 125 e x3 = 323, que estão na base 10, para base
2.
3. Em uma Universidade a senha de acesso à rede wi- é composta por 4 dígitos que es-
tão na base 10, porém por motivos de segurança esses números são apresentados aos
alunos na forma binária separados por duplas. A senha segue o formato ABCD onde,
AB = (1010111)2 e CD = (10110)2 . Podemos armar que a senha wi- é:
(b) 8722
(c) 8272
Respostas
1. x1 = (34)10 = (100010)2 x2 = (125)10 = (1111101)2 x3 = 323 = (101000011)2
3. (b)
Erro de Truncamento
É quando simplesmente desprezamos uma parte da mantissa e camos com o que restou.
Exemplo: π = 3, 14159265... truncando o valor de π com 4 casas decimais caria π = 3, 1415.
Erro de Arredondamento
É quando é feito a "aproximação"do número obtido por excesso (para cima), por falta (para
baixo) ou simétrico (mais próximo). Assim, em linhas gerais, para arredondar um número para
o número de máquina mais próximo, na base 10, devemos apenas observar o primeiro dígito a
ser descartado. Se este dígito é menor que 5, deixamos os dígitos inalterados; e se é maior ou
igual a 5, devemos somar 1 ao último dígito que restou.
Exemplo: O número 0,66666..., arrendondando cará 0,66667
A diferença entre o valor aproximado e o valor exato de um número pode ser medida pelo
erro absoluto ou pelo erro relativo.
|x − x0 | |x − x0 |
Erro Absoluto : EA = |x − x0 | Erro Relativo : ER = ou ER =
|x0 | |x|
Exemplo 1.5.1
Resolução: EA = |0, 003 − 0, 002| = 0, 001 O que pode parecer um valor pequeno de erro.
Agora vamos comparar com o erro relativo.
Exemplo 1.5.2
|101000 − 10000|
ER = = 0, 01 → 1%, neste caso temos um erro relativamente pequeno.
|10000|
1.6 Exercícios
n=0
n!
n 6= 0 e 0! = 1. Calcule o erro absoluto e o erro relativo do número e = 2, 718282, em
5
1
relação à sua aproximação dada por:
X
n=0
n!
Respostas
|31415 − 31416|
1. a) 31415 m2 b) 31416 m2 c) ER = = 0, 000031831
|31416|
2. f (1.25) = −1, 843200; f (1.248) = −1, 843703 ER = 0, 0002728
3. ER = 0, 019 4. EA = ER = 0, 000594
Exemplo 2.1.1
3x + 2x2 − 4x3 = 1
1
Resolva por escalonamento o sistema x1 + x2 + 2x3 = 2
4x + 3x + 2x = 3
1 2 3
0 0 4 0
7
Capítulo 2. Resolução de Sistemas Lineares (Unidade 04)
3x + 2x2 − 4x3 = 1
1
Agora vamos reescrever o sistema: 1/3x2 + 10/3x3 = 5/3
= 04x3
Resolvendo o sistema da terceira linha para primeira, chegamos em: (x1 , x2 , x3 ) = (−3, 5, 0)
2.2 Exercícios
x1 − x 2 + 2x3 − x4 = −8
2x − 2x
+ 3x3 − 3x4 = −20
1 2
2. Escalonar e resolver o sistema
x1 + x2 + x3 = −2
x − x
+ 4x3 + 3x4 = 4
1 2
3x + y + 2z = −1
3. Considere o seguinte sistema x + y = 2 .
2x + 2y −
z = 1
(a) x + z = −1 (c) y + z = −2
(b) z − y = −4 (d) x + y = 0
x + y + z = 7
6. Determine x + y · z para 2x + y − z = 0
x − 2y
+ 2z = 2
2x + y − z = 0
7. Use a resolução do sistema x − y + z = 3 , para se calcular cos(x + y) + z .
3x − y
+ 2z = 6
(a) 0 Vídeo
(b) 2
(c) −1
(d) 1
3x + 3y + z = 7
8. Resolver, usando o método de Eliminação de Gauss, o sistema 2x + 2y − z = 3
x − y + 5z = 5
Respostas
1. (a) (x1 , x2 , x3 , x4 ) = (−1, 2, 0, 1) (b) (x, y, z) = (1, 1, −1)
2. (x1 , x2 , x3 , x4 ) = (−7, 3, 2, 2)
3. x + y − z = −1
4. (c) 7. (b)
5. (a)
6. x + y · z = 11.13 8. (d)
Teorema 3.1.1
Seja f uma função contínua no intervalo [a,b], tal que
f(a).f(b) < 0. Então essa função possui pelo menos uma raíz real nesse intervalo
[a,b].
Considere o intervalo [a,b] para o qual f (a).f (b) < 0. No método da bisseccão calculamos o
a+b
valor da função f (x) no ponto médio: x1 = . Portanto teremos com isso três possibilidades:
2
1. f (x1 ) = 0, então x1 é a raíz. Nada a se fazer!!
Vídeo
2. f (a).f (x1 ) < 0, então f (x) tem um zero (raíz) em
[a, x1 ]. Repete-se o processo.
10
Capítulo 3. Zeros de Funções Reais (Unidade 05)
A repetição desse método é chamado iteração. A quantidade n de iterações é dada por:
log(b − a) − log(E)
n>
log(2)
onde E é o erro.
Exemplo 3.2.1
Determinar um valor aproximado da raiz quadrada de 5. Com erro menor ou igual a 0,01.
√
Portanto, 5 = 2, 2421875
Exemplo 3.2.2
2. Estime o valor de uma das raízes da função f (x) = x5 + x − 1 no intervalo [0,1] com erro
de 0,01.
3. Calcule araíz positiva da função f (x) = x2 + ln(x) com erro menor igual a 0,01 no
1
intervalo , 1 .
2
Vídeo
4. Calcule a raíz positiva da função
f (x) = x − cos(x) com erro menor igual a 10−1 .
Obs.: Construa o gráco com duas funções separadas para
estimar o intervalo.
5. Calcular a raíz da função f (x) = x3 − 10 com erro menor igual a 0,1 no intervalo [2,3].
1
6. Calcule uma raíz real de f (x) = 3x − cos(x) com erro de 10 −2
no intervalo 0, .
2
Respostas
1. x7 = 2,648438 4. x4 = 0,656250
2. x7 = 0,757813 5. x4 = 2,187500
3. x6 = 0,648438 6. x6 = 0,320312
Exemplo 3.4.1
√
Estimar o valor de 3, pelo método de Newton com erro de 0,01.(Usar 6 casas decimais)
Resolução: Usando a função f (x) = x2 − 3, no intervalo [1, 2], vamos usar o número 2
como valor de x0 . Devemos ter em mãos a derivada de f . f 0 (x) = 2x.
f (x0 ) f (2) 1
x 1 = x0 − = 2− =2− x1 = 1.75
f 0 (x0 ) f 0 (2) 4
O critério de parada será quando |xk − xk+1 | < e , então temos, |x3 − x2 | < 0.01
Exemplo 3.4.2
Exemplo 3.4.3
Em engenharia ambiental, a seguinte equação pode ser usada para calcular o nível de
concentração de oxigênio c num rio, em função da distância x, medida a partir d local de
descarga de poluentes:
c(x) = 10 − 20(e−0,2x − e−0,75x )
Calcule, usando o método de Newton, a distância para a qual o nível de oxigênio desce
para o valor 5. Utilize para aproximação inicial o valor x0 = 1.
A derivada desta nova função será f 0 (x) = 4.e−0,2x − 15.e−0,75x . Agora aplicamos na fórmula
de iteração.
f (x0 ) f (x1 )
x1 = x0 − = 0.494228 x2 = x1 − = 0.596414
f 0 (x0 ) f 0 (x1 )
f (x2 ) f (x3 )
x3 = x2 − = 0.602337 x4 = x3 − = 0.602355
f 0 (x2 ) f 0 (x3 )
2. Em um determinado
! circuito elétrico!a corrente pode ser calculada pela seguinte expressão:
π π
I = 5.e−5t . sin + 5. sin 5t − , onde t é tempo em segundos. Determine um valor
4 4
de t para qual o valor da intensidade da corrente é nula. Usar t0 = 1s e 4 casas decimais.
3. Dispõe-se de uma lâmina retangular de 10 cm × 16 cm, para construir uma caixa re-
tangular sem tampa, cortando um quadrado de igual tamanho em cada uma das quinas.
Estimar um valor para o lado do quadrado de tal forma que o volume de caixa seja de
100 cm3 . Utilizar o método de Newton-Raphson com precisão de 0,001 cm e valor inicial
de x0 = 1.
Vídeo
6. Faça a estimativa da raíz usando o método de Newton para f (x) = ex − 4x2 com erro
< 10−1 . Considere x0 = 1.
7. Sabemos que f (x) = x3 − 3x + 1 possui uma raíz no intervalo [0,1]. Estime, pelo método
de Newton-Raphson, o valor desta raíz. Considere ≤ 0, 002.
8. Calcular pelo menos uma raiz real das equações abaixo, com ≤ 10−3 , usando o método
de Newton.
a) f (x) = 2x − sen(x) + 4
b) f (x) = ex − tg(x)
c) f (x) = 10x + x3 + 2
9. Dada a função f (x) = e−x + x2 − 2. Estime o valor da raíz com erro ≤ 0, 05, no intervalo
[1,2].
10. Usando o método de Newton-Raphson, com erro inferior a 10−2 , determinar uma raiz das
seguintes equações: Obs.: Colocar a calculadora
em radianos para as trigonométricas.
3
a) 2x − tg(x) = 0 No intervalo 1;
2
b) 5x3 + x2 − 12x + 4 = 0 No intervalo [0; 1]
7
c) sen(x) − ex = 0 No ponto x0 =
2
Respostas
1. a) x = 1, 218636 m b) x = 0, 345974 m
2. t = 0, 7881
3. x = 0, 839019 cm ou x ≈ 0, 84 cm
4. x = 4, 095462
5. θ = 1, 275698 rad
6. x = 0, 714806
7. x = 0, 347296
9. x = 1, 315974
17
Capítulo 4. Interpolação Polinomial (Unidade 06)
Exemplo 4.1.1
Dados os pares de pontos (−1, 15) , (0, 8) e (3, −1) determinar o polinômio de interpolação
para a função denida por esse conjunto de pares de pontos.
a − a1 + a2 = 15
0
a0 = 8 ⇒ a0 = 8, a1 = −6 e a2 = 1
a + 3.a + 9.a = −1
0 1 2
∴ p(x) = 8 − 6x + x2
4.2 Exercícios
Respostas
5 x
(a) p(x) = −
3 3
(b) p(x) = −0.5x2 + 0.5x + 2
i=0
onde:
n
Y (x − xi )
L(xi ) =
i=0, i6=j
(xj − xi )
Exemplo 4.3.1
(x − x1 ).(x − x2 ) (x − 0).(x − 2) x2 − 2x
L0 (x) = = =
(x0 − x1 ).(x0 − x2 ) (−1 − 0).(−1 − 2) 3
(x − x0 ).(x − x2 ) (x − (−1)).(x − 2) x2 − x − 2
L1 (x) = = =
(x1 − x0 ).(x1 − x2 ) (0 − (−1)).(0 − 2) −2
(x − x0 ).(x − x1 ) (x − (−1)).(x − 0) x2 + x
L2 (x) = = =
(x2 − x0 ).(x2 − x1 ) (2 − (−1)).(2 − 0) 6
x2 − 2x x2 − x − 2 x2 + x
p(x) = f (x0 ).L0 (x)+f (x1 ).L1 (x)+f (x2 ).L2 (x) = (1). +(−1). +(1). =
3 −2 6
p(x) = x2 − x − 1
Exemplo 4.3.2
1
Encontre o polinômio que interpola f (x) = nos pontos x0 = 2 x1 = 2.5 x2 = 4.
x2
1 1 1
Resolução: f (x0 ) = = 0.25, f (x1 ) = = 0.16 e f (x2 ) = = 0.0625
4 6.25 16
(x − x1 ).(x − x2 ) (x − 2.5).(x − 4)
L0 (x) = = = x2 − 6.5x + 10
(x0 − x1 ).(x0 − x2 ) (2 − 2.5).(2 − 4)
4.4 Exercícios-Lagrange
2. Obtenha o polinômio interpolador de Lagrange para a função f (x) = Ln(x) nos pontos
x0 = 1, x1 = 1.2 e x2 = 1.5.
x −1 0 3
f (x) 15 8 −1
5. Considere a tabela:
x 1 3 4 5
f (x) 0 6 24 60
6. Usar o método de Lagrange para determinar o polinômio interpolador a partir dos dados
fornecidos na tabela a seguir.
Respostas
11 2 13
1. p(x) = x − x−2
24 12
x f (x)
2. p(x) = −0, 335x2 + 1, 6485x − 1, 3135
0 0
0,2 2,008 3. p(x) = 0.5x2 − 0.5x − 1
6. p(x) = x3 + 10x
O método consiste em interpolar um conjunto de pontos por meio de um polinômio que usa
diferenças divididas como coecientes.
Polinômio de Newton
Exemplo 4.5.1
x −1 0 3
f (x) 15 8 −1
Determine o polinômio de interpolação na forma de Newton.
Resolução:
x f (x) ∆1 ∆2
-1 15
8 − 15
= −7
0 − (−1)
−3 − (−7)
0 8 =1
3 − (−1)
−1 − 8
= −3
3−0
3 -1
∴ p(x) = x2 − 6x + 8
Exemplo 4.5.2
Considere a tabela:
x 1 3 4 5
f (x) 0 6 24 60
Resolução:
x (x0 ) ∆1 ∆2 ∆3
x0 = 1 0
6−0
=3
3−1
18 − 3
x1 = 3 6 =5
4−1
24 − 6 9−5
= 18 =1
4−3 5−1
36 − 18
x2 = 4 24 =9
5−3
60 − 24
= 36
5−4
x3 = 5 60
∴ p(x) = x3 − 3x2 + 2x
b) f (4.5) = (4.5)3 − 3.(4.5)2 + 2(4.5) = 91, 125 − 60, 75 + 9 = 39, 375 ∴ f (4, 5) = 39, 375
Vídeo
x -2 -1 1 2
f (x) 0 1 -1 0
x 0 1 1,5 2,5
f (x) 1.0 0.5 0.4 0.286
5. A tabela abaixo apresenta a população dos EUA (em milhões) de 1940 a 1980.
Tempo (s) 1 3 5 7 20
Velocidade (cm/s) 800 2310 3090 3940 8000
Respostas
Integral numericamente uma função y = f (x) num dado intervalo [a, b] é integrar um
polinômio Pn (x) que aproxime f (x) no dado intervalo.
b
(b − a)
Z
f (x)dx ≈ [f (a) + f (b)]
a 2
Uma maneira de se calcular usando a regra dos trapézios é subdividindo o intervalo [a, b]
em n subintervalos de amplitude iguais a h e a cada subintervalo aplica-se a regra.
Vídeo
Z b
h
f (x) ≈ . y0 + 2y1 + 2y2 + ... + 2yn−1 + yn
a 2
26
Capítulo 5. Integração Numérica (Unidade 07)
Exemplo 5.1.1
Determine o volume do sólido gerado pela rotação da função f (x) = sen(x) + 2 em torno
do eixo x no intervalo de 0 ≤ x ≤ 2π .
Resolução: Z b
Lembrando que a fórmula do volume é V = π. [f (x)]2 dx, então o volume pedido é
a
Z 2π
V = π. [sen(x) + 2]2 dx.
0
2π − 0 π
Para o intervalo dado vamos escolher subdividí-lo em 4 partes. O que ca h = =
4 2
obs.: A função usada para o método dos trapézios será g(x) = [f (x)]2 = (sen(x) + 2)2 .
i xi f (xi )
π π
1 f (x1 ) = (sen + 2)2 = 9
2 2
2π
π2
Z
2 π/2
V = π. [sen(x) + 2] dx = π. .(4 + 2.(9 + 4 + 1) + 4) = .36 = 9π 2 u.v
0 2 4
Exemplo 5.1.2
Z 3.6
1
Calcular a integral dx utilizando a regra dos Trapézios composta e subdividindo o
3.0 x
intervalo em 6 partes.
Resolução:
b−a 3.6 − 3.0 0.6
Primeiro calcularemos o incremento dos subintervalos. h = = = = 0.1
n 6 6
Vamos construir uma tabela para saber o valor das bases dos trapézios.
i xi yi = f (xi )
0 x0 = 3.0 y0 = 0.333333
1 x1 = 3.1 y1 = 0.322581
2 x2 = 3.2 y2 = 0.312500
3 x3 = 3.3 y3 = 0.303030
4 x4 = 3.4 y4 = 0.294118
5 x5 = 3.5 y5 = 0.285714
6 x6 = 3.6 y6 = 0.277778
Z 3.6
1 h
dx = y0 + 2y1 + 2y2 + 2y3 + 2y4 + 2y5 + y6 = 0, 182350
3.0 x 2
Vídeo
Z 6
1. Calcule usando a regra dos Trapézios o valor de ex .cos(x) dx. Sub-
2
dividindo em 5 partes.
Z 2
1
2. Calcule a integral aproximando seu valor pela regra dos Trapézios dx com
0 x2 +4
n = 8.
3. Dada a tabela:
Z 0,8
Determine o valor de x.ex dx usando a regra dos Trapézios.
0
Z 5
x
4. Use a regra dos Trapézios para calcular √ dx com n = 10.
2 x2+4
Respostas
1. 158,347063
2. 0,392374
3. 0,5649
4. 2,555605