MS211C/F - C Alculo Num Erico PROVA 1 (12/05/2022) Nome: RA

Fazer download em pdf ou txt
Fazer download em pdf ou txt
Você está na página 1de 2

MS211C/F – Cálculo Numérico PROVA 1 (12/05/2022)

Nome: RA: ~~~~~~

Justifique e explique todos os passos de cada resolução.


Respostas puramente numérica não serão consideradas.

13 2
1. Seja p : R → R, dada por p(x) = −1 + 3
x − 2x4 .

(a) Partindo de x0 = 1, mostre que a sequência exata de iterandos gerada pelo método de Newton
não converge a um zero de p.
(b) Um estudante, duvidando do afirmado no item acima, resolveu fazer um experimento numérico
no Octave. As primeiras 21 iterações de Newton resultaram nos dados exibidos na tabela a
seguir.
k xk f (xk ) k xk f (xk ) k xk f (xk )
1 −1.00000 1.3333 · 100 8 1.00062 1.3337 · 100 15 1.67714 −4.6349 · 100
2 1.00000 1.3333 · 100 9 −1.02894 1.3460 · 100 16 1.47740 −1.0700 · 100
3 −1.00000 1.3333 · 100 10 5.61611 −1.8539 · 103 17 1.39505 −1.4173 · 10−1
4 1.00000 1.3333 · 100 11 4.26129 −5.8178 · 102 18 1.38033 −4.0848 · 10−3
5 −1.00000 1.3333 · 100 12 3.26184 −1.8130 · 102 19 1.37988 −3.7517 · 10−6
6 1.00000 1.3333 · 100 13 2.53481 −5.5725 · 101 20 1.37988 −3.1761 · 10−12
7 −1.00001 1.3333 · 100 14 2.02039 −1.6637 · 101 21 1.37988 1.7764 · 10−15

Usando um Sistema de Ponto Flutuante com 4 dı́gitos significativos, compute a primeira


iteração do método de Newton. Explique o que pode ter acontecido no experimento acima,
que foi conduzido em precisão dupla.

2 +2x
2. Seja f (x) = x − e−x .

(a) Prove que f tem pelo menos um zero. Exiba um intervalo finito que contenha um único zero
de f e tenha comprimento menor ou igual a 0.5.
(b) Aplique um dos métodos visto nesta disciplina para estimar o valor do zero de f , com pelo
menos 4 casas decimais corretas.
   
2 4 1 5
3. Considere o sistema Ax = b, com A =  5 −2 3  e b =  −13 .
1 6 4 0
Abaixo está exibida a sequência de operações elementares utilizadas no escalonamento de A, com
pivoteamento parcial.

     
`2 2 4 1 5 −2.0 3.0 5 −2.0 3.0
`1  5 −2 3  ∼ `2 − 0.4`1  2 4.0 1.0  ∼ `3  0 4.8 −0.2 
1 6 4 `3 − 0.2`1 1 6.0 4.0 `2 0 6.4 3.4
   
5 −2.0 3.0 5 −2.0 3.0
∼  0 6.4 3.4  ∼  0 6.4 3.4 
`3 − 0.75`2 0 4.8 −0.2 0 0 −2.75

(a) Usando estes resultados, exiba as matrizes L, U e P , que correspondem à decomposição LU


de A. Justifique.
(b) Resolva o sistema linear Ax = b pela fatoração LU , aproveitando os cálculos já fornecidos.

4. Sejam A ∈ Rn×n e b ∈ Rn , e considere o problema de aproximar a solução do sistema linear


Ax = b.

(a) Caso a matriz A fosse a seguir, seria possı́vel assegurar, de alguma forma, a convergência
da sequência gerada pelo método de Jacobi? Se sim, você pode dizer alguma coisa sobre a
velocidade com que o erro reduziria a cada passo?
 
10 3 1 5
2 5 1 1
A=
1
.
4 0 9
3 1 8 1
(b) Escreva um algoritmo para o método de Gauss-Seidel. Deixe claro quais outras informações
iniciais, além de A e b, seu algoritmo necessita. Qual o custo computacional, em termos do
número de operações aritméticas de ponto flutuante, de cada iteração do algoritmo?

Você também pode gostar