Atividade 01

Disciplina: Métodos Numéricos

Instruções:

• A lista deve ser realizada à mão usando folhas A4 ou de caderno. Não há necessidade
de escrever os enunciados das questões;

• Essa lista pode ser realizada INDIVIDUALMENTE ou em DUPLA, com o nome dos
alunos “escritos” no topo de cada folha. Mesmo fazendo em dupla, CADA aluno
deverá enviar seu trabalho;

• A lista deverá ser digitalizada (escaneada ou fotografada) em um único arquivo no

formato PDF antes de enviar pelo CANVAS. Há aplicativos de celulares que fazem isto,
por exemplo, Tiny Scanner, CamScanner e Office Lens.

• Se o tamanho da lista escaneada ficar muito grande, há programas online que

ajudam a reduzir o tamanho de arquivo PDF (smallpdf.com/pt/compressor-de-pdf).

• Atividades entregue FORA de TAREFA/TESTE no CANVAS, não serão aceitas. Ou seja,

atividades enviadas pelo SGA, por e-mail, por links de armazenamento em Nuvens ou
pela caixa de entrada do Canvas NÃO serão avaliadas;

• Atividades ilegíveis ou impossibilitados de abrir, serão zeradas;

• E, é claro, atividades IGUAIS serão zeradas tanto de quem COPIOU e de quem DEIXOU
copiar;

• Bom aprendizado.

1|Página
 Veja o que pode e o que não pode ao enviar o arquivo
Assim não: tudo digitado Assim não: nome digitado e com recortes

Assim não: questões recortada É assim que deve ser: nome escrito e sem
recortes

Assunto: unidade 1

• Sistemas Lineares: métodos diretos;

• Sistemas Lineares: métodos iterativos

ATENÇÃO

• Onde estiver escrito “à mão”, você deverá fazer a mão.

• Onde estiver escrito (Computer-Aided Software), você poderá usar alguma ferramenta
computacional como o Excel, Matlab, Scilab, Octave, R, Python ou outro.
• No caso de usar a ferramenta tecnológica, basta dar um print da tela do computador e colocar
como anexo no final na lista;

2|Página
 LISTA 01 - CÁLCULO NUMERICO - PUC MINAS
Aluno 01________________________________ Aluno 02 ______________________________

QUESTÃO 1) (à mão)

Resolva o sistema linear a seguir mantendo quatro casas decimais e, logo em seguida, dê o vetor
resíduos

RESPOSTA:

3|Página
 LISTA 01 - CÁLCULO NUMERICO - PUC MINAS
Aluno 01________________________________ Aluno 02 ______________________________

QUESTÃO 2) (à mão)

Resolva o sistema abaixo usando o método da pivotação parcial (ver ANEXO A) mantendo quatro
casas decimais

−2𝑥1 + 𝑥2 + 𝑥4 = −2
4𝑥 + 3𝑥2 + 3𝑥3 + 𝑥4 = 23
{ 1
−8𝑥1 + 7𝑥2 + 9𝑥3 + 5𝑥4 = −1
6𝑥1 − 7𝑥2 + 10𝑥3 + 8𝑥4 = −18

RESPOSTA:

4|Página
 LISTA 01 - CÁLCULO NUMERICO - PUC MINAS
Aluno 01________________________________ Aluno 02 ______________________________

QUESTÃO 3) (usando CAS)

Dê a solução do sistema abaixo usando o programa Octave (ver ANEXO B)

−2𝑥1 + 𝑥2 + 𝑥4 + 3𝑥5 − 2𝑥6 = 8.32

4𝑥1 + 3𝑥2 + 3𝑥3 + 𝑥4 + 2𝑥5 + 0.5𝑥6 = 5.95
−8𝑥1 + 7𝑥2 + 9𝑥3 + 5𝑥4 + 2.5𝑥5 − 3𝑥6 = 2.66
6𝑥1 − 7𝑥2 + 10𝑥3 + 8𝑥4 = 48.53
−𝑥1 + 2𝑥2 − 3𝑥3 + 𝑥4 + 0.5𝑥5 + 1.5𝑥6 = −15.65
{3𝑥3 + 5𝑥4 − 2𝑥5 + 4𝑥6 = −7.29

[0.48458 ; -3.66932 ; 2.65861 ; -0.83110 ; 4.11721 ; -0.71897]

5|Página
 LISTA 01 - CÁLCULO NUMERICO - PUC MINAS
Aluno 01________________________________ Aluno 02 ______________________________

QUESTÃO 4) (usando CAS)

Dê o vetor resíduos para a solução x = [-4, 2, 0.5, -3] apresentada para o sistema linear a seguir.

+0.531𝑥1 +0.742𝑥2 −1.251𝑥3 −2.104𝑥4 = −5.10

+7.428𝑥1 −1.161𝑥2 +0.399𝑥3 +1.516𝑥4 = 13.31
{
+2.126𝑥1 +3.071𝑥2 +0.455𝑥3 −2.706𝑥4 = 1.55
−1.007𝑥1 +6.418𝑥2 −1.442𝑥3 +0.774𝑥4 = 16.77

6|Página
 LISTA 01 - CÁLCULO NUMERICO - PUC MINAS
Aluno 01________________________________ Aluno 02 ______________________________

QUESTÃO 5) Resolva o sistema abaixo usando o método de Jacobi-Richardson com solução inicial
igual a x0 = (2, 2, 2).

𝑥1 − 2𝑥2 + 4𝑥3 = 14.2

{2𝑥1 + 𝑥2 − 𝑥3 = −4.1
𝑥1 + 3𝑥2 + 2𝑥3 = 16.7

a) Pelo critério de convergência, é garantido que o sistema acima irá convergir para algum
resultado?
b) (à mão) Mostre os cálculos para as três primeiras iterações.
c) (usando CAS) Continue até a 10ª iteração. Você acredita que os resultados estão
convergindo para uma solução do sistema?
d) (usando CAS) Em alguns sistemas lineares é possível rearranjar as linhas e/ou colunas de
forma que o critério de convergência para o método iterativo seja satisfeito. Faças as
modificações necessárias no sistema linear acima e o resolva (usando o Excel) usando uma
tolerância de 0,001.

Solução --------------------------------------------------------------------------------------------------------------------
a). x = [-75.5 39.1 46.275]t
b). x = [41463.59 33411.13 -30669.21]t
c). solução = -0.80 ; 5 ; 2.5

7|Página
 LISTA 01 - CÁLCULO NUMERICO - PUC MINAS
Aluno 01________________________________ Aluno 02 ______________________________

QUESTÃO 6) (usando CAS)

Use a função triangsup que está no ANEXO C para resolver o sistema

2𝑥1 + 4𝑥2 + 𝑥3 − 𝑥4 + 2𝑥5 + 𝑥6 = 34

7𝑥2 + 𝑥3 + 2𝑥4 + 6𝑥5 + 2𝑥6 = 69
4𝑥3 − 3𝑥4 − 4𝑥5 + 9𝑥6 = −23
𝑥4 + 3𝑥5 + 2𝑥6 = 35
2𝑥5 − 6𝑥6 = −4
{ 3𝑥6 = 12
RESPOSTAS
a) xs = [24, 8, -12, -2]t
b) xs =[3, 2, -7, -3, 10, 4]

8|Página
 LISTA 01 - CÁLCULO NUMERICO - PUC MINAS
Aluno 01________________________________ Aluno 02 ______________________________

QUESTÃO 7) (usando CAS)

No circuito elétrico abaixo temos duas fontes de tensões (10V e 4V) e alguns resistores. O objetivo
é obter as correntes elétricas i1, i2 e i3 mostradas em cada malha do circuito.

Pela Lei de Kirchoff para Tensões (LKT), a soma algébrica das tensões em qualquer caminho fechado
(laço ou malha) é zero. Pela Lei de Ohm, sabe-se que a tensão V em resistor com resistência R e
corrente 𝐼 passando por ele, é dada por 𝑉 = 𝑅𝐼. Aplicando essas duas leis no circuito acima, temos
o sistema linear a seguir para obter as correntes i1, i2 e i3

Determine valor das correntes i1, i2 e i3 usando o método iterativo de Jacobi com tolerância de
0,0001. Se preferir, troque i1, i2 e i3 por x1, x2 e x3.

RESPOSTA:
A solução que deverá ser próxima da solução exata [2, 1, 1] amperes

9|Página
 LISTA 01 - CÁLCULO NUMERICO - PUC MINAS
Aluno 01________________________________ Aluno 02 ______________________________

QUESTÃO 8) Aplicação em TRELIÇAS (usar CAS)

Uma treliça é uma estrutura resultante das combinações de diversas barras entrelaçadas. Em uma
treliça planar as barras se conectam por nós. Considerando que um sistema estático, é possível
calcular as tensões que ocorrem nas barras com base nas forças externas e de apoio sobre o sistema.
As barras podem estar em estados de compressão ou tração, e as estimativas de tensões obtidas
podem ser usadas para avaliar a eficácia da estrutura frente a esforços aplicados ao sistema.

Para o cálculo das tensões considera-se que cada nó estático tem resultante de forças nulo. Ao
analisar o sistema completo obtém-se geralmente um sistema linear com 2 vezes os números de
nós, sendo que para cada nó consideramos duas equações resultante vertical nula e resultante
horizontal nula.

A figura abaixo é um exemplo de treliça com 4 nós, 5 barras, 1 apoio fixo (símbolo triangular com
retângulo cinza abaixo) e 1 apoio vertical com movimento livre na horizontal (símbolo triangular
com círculos abaixo) e uma força de 10KN no nó 3.

(www.ime.usp.br/~map3121/2017/map3121/programas/EP1-trelica.pdf)

As forças atuantes em cada nó estão mostradas na figura a seguir. As forças obliquas (f1 e f4) são
substituídas pelas suas correspondentes componentes horizontais e verticais.

10 | P á g i n a
 LISTA 01 - CÁLCULO NUMERICO - PUC MINAS
Aluno 01________________________________ Aluno 02 ______________________________

Usando as leis da estática, é possível obter o sistema linear a seguir que permite encontrar as forças.
Note que esse sistema é esparso com muitos elementos nulos na matriz.

(usando CAS) Usando o método de Gauss-Seidel com tolerância de 0,01 e começando com solução
(F1, F2, F3, f4, f5, f6, f7, f8) = (0, 0, 0, 0, 0, 0, 0, 0) dê os valores aproximados das forças.

Resposta:

11 | P á g i n a
 LISTA 01 - CÁLCULO NUMERICO - PUC MINAS
Aluno 01________________________________ Aluno 02 ______________________________

ANEXOS

ANEXO A – Escalonamento de Gauss om pivotação parcial

A eliminação gaussiana com pivotamento parcial consiste em fazer trocas de linhas de forma que o
pivô a ser escolhido seja o elemento de maior módulo em cada coluna. Assim, evita-se que tenha
pivô nulo e os efeitos dos erros de arredondamentos serão menores nesse método

Como exemplo, vamos resolver o sistema a seguir usando o escalonamento de Gauss com
pivotamento parcial

Na 1ª coluna, o maior módulo está na 2ª linha

(a21 = 2) e terceira linha (a32 = 1). Podemos
tanto a 2ª ou 3ª linha para permutar. Foi
escolhida a 2ª linha

𝐿1 ↔ 𝐿2

𝐿2 ← 𝐿2 + (−1/2)𝐿1
𝐿3 ← 𝐿3 + (−1)𝐿1

Na 2ª coluna, o maior módulo

está na 3ª linha com a31 = 1 (a
𝐿2 ↔ 𝐿3
1ª linha foi descartada, já foi
usada).vamos permutar as
linhas 2 e 3

𝐿3 ← 𝐿3 + (−1/2)𝐿2

12 | P á g i n a
 LISTA 01 - CÁLCULO NUMERICO - PUC MINAS
Aluno 01________________________________ Aluno 02 ______________________________

ANEXO B - Resolvendo sistema linear usando escalonamento

2𝑥1 + 4𝑥2 − 2𝑥3 + 4𝑥4 = 46

4𝑥 − 𝑥2 + 3𝑥3 − 5𝑥4 = −21
{ 1
𝑥1 + 2𝑥2 + 2𝑥3 + 3𝑥4 = 36
2𝑥1 + 4𝑥2 − 2𝑥3 − 3𝑥4 = −3

** SUGESTÃO: Use Baixe o programa em https://octave.org/ ou usar uma

versão on-line disponível em https://octave-online.net/

# matriz com os coeficientes

A = [2, 4, -2, 4 ; 4, -1, 3, -5 ; 1, 2, 2, 3 ; 2, 4, -2, -3];

# vetor com os termos independentes

b = [46;-21;36;-3];

# matriz M aumentada
M = [A b]

# rref dá a matriz escalonada reduzida com a solução na última

# coluna (método de Gauss-Jordan)
rref(M)

# a solução xs é a 5ª coluna de rref(M)

xs = rref(M)(:,5)

# calculando o vetor residuos

r = b - A*xs

13 | P á g i n a
 LISTA 01 - CÁLCULO NUMERICO - PUC MINAS
Aluno 01________________________________ Aluno 02 ______________________________

ANEXO C – função para resolver sistema triangular

function [sol] = triangsup(x,b)

# Resolucao de sistema triangular superior 4x4

# ENTRADA
# # x = matriz com os coeficientes
# # b = vetor dos termos independentes
# SAIDA
# # sol = vetor solucao

s(4) = b(4) / x(4,4);

s(3) = (b(3) - ( x(3,4:4) * s(4:4)' )) / x(3,3);
s(2) = (b(2) - ( x(2,3:4) * s(3:4)' )) / x(2,2);
s(1) = (b(1) - ( x(1,2:4) * s(2:4)' )) / x(1,1);
sol = s;

endfunction

OBS: O apostrofo (') usado acima é para transpor o vetor/matriz. Se você “copia” os comandos acima e “colar” no
Matlab/Octave/Scilab, costuma esse apóstrofo (') e o sinal de menos (“-“) saírem errados, nesse caso, o programa dá
problemas. Se isso ocorrer, basta apagar e digitar novamente eles novamente dentro do programa Octave.

Usando a função triangsup()acima para resolver o sistema

2𝑥1 + 3𝑥2 + 6𝑥3 + 𝑥4 = −2

8𝑥2 + 4𝑥3 + 6𝑥4 = 4
{
−𝑥3 + 3𝑥4 = 6
2𝑥4 = −4

# Entrada dos dados

clc # limpa a tela
clear # reinicializa as variáveis
A = [2 3 6 1 ; 0 8 4 6 ; 0 0 -1 3 ; 0 0 0 2] # entra com matriz A
b = [-2 ; 4 ; 6 ; -4] # entra com vetor termos independentes
triangsup(A, b) # usar a funcao criada

14 | P á g i n a