MN PPF 21jan

Classificação:
Departamento de Matemática da Universidade de Aveiro

Métodos Numéricos 2019/2020 (1¯o Semestre)
Prova Prática Final - 21/01/2020 Duração: 2h
NOME: No :
CURSO: Turma Prática:
Todas as respostas devem incluir os comandos MatLab usados, as justificações e os resultados pedidos.
1. Sejam as funções cos(x) e sin(x), e pn (x) e qn (x) os respetivos polinómios de MacLaurin de grau n ∈ N.
(a) Utilize a função taylor do MatLab para determinar os polinómios p8 (x) e q9 (x).
Resposta (Cálculos efectuados, comandos usados e polinómios p8 (x) e q9 (x).)
√ √
(b) Tendo em conta que cos(π/6) = 3/2√e sin(π/3) = 3/2, utilize dois processos de aproximação
para calcular um valor aproximado de 3, um recorrendo ao polinómio p8 e outro ao polinómio q9 .
√
Resposta (Cálculos efectuados, comandos usados e as duas aproximações de 3 em formato long.)
√
(c) Compare as duas aproximações obtidas na alı́nea anterior com o valor de 3 dado pelo MatLab e
diga qual delas é a mais precisa. O que terá contribuı́do para essa maior precisão? Justifique.
Resposta (Cálculos efectuados, comandos usados, aproximação mais precisa e justificação.)
2. Considere o sistema linear Ax = b, onde
   
0 −1 −1 4 −100
 −1 4 −1 
0  e b =  −100  .
 
A=  −1 4 0 −1   100 
4 −1 −1 0 100
(a) Resolva o sistema pelo método da fatorização LU .
Resposta (Cálculos efectuados, comandos usados, factorizações e solução x no formato short.)
(b) Suponha que A é perturbada pela matriz δA = 10−3 A, tal que Ap = A + δA. Aplicando o
método de eliminação de Gauss com escolha parcial de pivot determine a solução xp do sistema
||A − Ap ||1 ||x − xp ||1
perturbado Ap xp = b, calcule e . Compare os valores destas perturbações e
||A||1 ||x||1
conclua, justificando, se o sistema Ax = b é bem ou mal condicionado.
||A−Ap ||1 ||x−xp ||1
Resposta (Comandos usados, xp , ||A||1 e ||x||1 , no formato short; conclusão justificada.)
(c) Obtenha um sistema W x = s equivalente ao sistema dado e para o qual esteja garantida a
convergência dos métodos iterativos de Jacobi e Gauss-Seidel. Considere como aproximação
T
inicial x(0) =

1 1 1 1 e compare a rapidez de convergência desses métodos,
determinando, para cada um deles, o número de iterações necessárias para se obter o valor de
x̃, tal que, ||x − x̃||∞ < 2−48 .
Resposta (Comandos usados, justificação do sistema W x = s, justificação do critério de paragem,
número de iterações e x̃ no formato long, para cada método.)
3. A tabela apresenta a velocidade v de queda de um paraquedista em função do tempo t:
t (s) 1 3 5 7 20 25
v (m/s) 8.0 23.1 30.9 39.4 80.0 95.5
(a) Determine a parábola p2 (t) = a2 t2 + a1 t + a0 , com a2 , a1 , a0 ∈ IR, que melhor se ajusta aos dados da
tabela no sentido dos mı́nimos quadrados e, a partir daı́, estime o tempo t∗ para o qual a velocidade
atingiu o valor de 55 m/s.
Resposta (Comandos usados, expressão de p2 (t) e t∗ em formato short)
(b) Indique a sequência de comandos necessários para representar num gráfico em MatLab os pontos
tabelados e a parábola p2 (t) para t ∈ [1, 25]. Calcule em seguida o erro quadrático médio da
aproximação p2 (t) aos dados da tabela.
Resposta (Sequência de comandos usados e erro quadrático médio em formato short)
4. Considere a função f (x) = sin(x) + 3 cos(x) − 2 e seja α a raiz de f (x) = 0 no intervalo I = [−1, −0.5].
(a) Obtenha uma estimativa de α usando interpolação cúbica inversa, com valores nodais igualmente
espaçados em I.
Resposta (Cálculos efectuados e comandos usados; Justificação da possibilidade de usar interpolação
inversa; Estimativa de α no formato short.)
(b) Mostre que são satisfeitas em I as condições suficientes para a convergência do método de Newton.
Em seguida, aplique-o para fornecer uma aproximação xk de α, tal que |xk − xk−1 | < 0.5 × 10−2 .
Escolha x0 de forma adequada.
Resposta (Cálculos efectuados e comandos usados; Justificações; Valor de x0 , número de iterações
e aproximação pedida no formato short.)
Cotações:
1.(a) 1.(b) 1.(c) 2.(a) 2.(b) 2.(c) 3.(a) 3.(b) 4.(a) 4.(b)
1.5 1.0 1.5 2.0 2.0 2.0 2.5 2.5 2.0 3.0

MN PPF 21jan

Enviado por

Direitos autorais:

Formatos disponíveis

MN PPF 21jan

Enviado por

Dados do documento

Título original

Direitos autorais

Formatos disponíveis

Compartilhar este documento

Compartilhar ou incorporar documento

Opções de compartilhamento

Você considera este documento útil?

Este conteúdo é inapropriado?

Direitos autorais:

Formatos disponíveis

MN PPF 21jan

Enviado por

Direitos autorais:

Formatos disponíveis

Classificação:

Departamento de Matemática da Universidade de Aveiro

CURSO: Turma Prática:

Você também pode gostar