Buscar

Aula 01 CALCULO NUMÉRICO



Continue navegando


Prévia do material em texto

CÁLCULO NUMÉRICO 
Aula 1 – Introdução ao Programa de 
Computação Numérica 
PROF. Reinaldo Camargo 
Reinaldo.camargo@estacio.br 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
EMENTA 
 Introdução ao programa de computação numérica. 
 Erros. 
 Zeros de funções. 
 Resolução de sistemas de equações lineares. 
 Aproximação. 
 Integração numérica. 
 Resolução de Equações Diferenciais Ordinárias. 
 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
UNIDADE 1 - Introdução ao Programa de Computação Numérica (PCN) 
 1.1 Operações Aritméticas (Escalares, Vetores e Matrizes) 
 1.2 Gráficos 
1.3 Programação Estruturada 
1.4 Funções 
 
UNIDADE 2 - Teoria dos Erros 
 2.1.Origem - Erros de Arredondamento e Truncamento 
 2.2.Propagação de Erros 
 
UNIDADE 3 - Solução de Equações Transcendentes e Polinomiais - 
Raízes de equações 
 3.1.Introdução - Identificação de Métodos Numéricos 
 3.2.Métodos de Intervalos (Bisseção e Falsa Posição) 
3.3.Métodos de Aproximação (Newton-Raphson e Ponto Fixo) 
 3.4.Método da Secante 
 
CONTEÚDOS 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
CONTEÚDOS 
UNIDADE 4 - Sistemas de Equações Lineares 
4.1.Métodos Diretos (Gauss-Jordan e Decomposição LU) 
 4.2.Métodos Iterativos (Gauss-Jacobi e Gauss-Seidel) 
4.3.Implementação e Aplicação com apoio de PCN 
 
 UNIDADE 5 - Aproximação de Funções 
5.1.Interpolação Polinomial (Métodos de Lagrange e Newton) 
5.2.Ajuste de Funções (Funções Polinomiais e Linearizáveis) 
 5.3.Implementação e Aplicação com apoio de PCN 
 
UNIDADE 6 - Integração Numérica 
 6.1 Métodos de Interpolação (Newton-Cotes: Retângulos, Trapézios, Regras de Simpson) 
6.2. Métodos de extrapolação: Método de Romberg 
6.3 Implementação e Aplicação com apoio de PCN 
 
UNIDADE 7 - Resolução de Equações Diferenciais Ordinárias de 1a. Ordem. 
 7.1 Modelos Dinâmicos; Métodos de Euler 
 7.2 Método de Runge-Kutta (1a. e 2a. Ordens) 
7.3 Implementação e Aplicação com apoio de PCN 
 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
PROCEDIMENTOS DE AVALIAÇÃO 
• AV1, AV2, AV3 
• Média >= 6,0 
• 2 Av’s com nota >= 4 para não reprovar direto na disciplina 
• Trabalhos individuais e em grupo 
 
 
• BARROSO, Leônidas Conceição et al. Cálculo numérico: 
com aplicações. 2. ed. São Paulo: Harbra, 1987. 
• ARENALES, Selma Helena de Vasconcelos; DAREZZO, 
Artur. Cálculo numérico: aprendizagem com apoio de 
software. São Paulo: Thomson Learning, 2008. 
• RUGGIERO, Marcia A. Gomes; LOPES, Vera Lucia da 
Rocha. Cálculo numérico: aspectos teóricos e 
computacionais. 2. ed. São Paulo: Pearson, 2006. 
 
 
BIBLIOGRAFIA BÁSICA 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
PROCEDIMENTOS DE AVALIAÇÃO 
• AV1, AV2, AV3 
• Média >= 6,0 
• 2 Av’s com nota >= 4 para não reprovar direto na disciplina 
• Trabalhos individuais e em grupo 
 
 
• BARROSO, Leônidas Conceição 
 et al. Cálculo numérico: com 
 aplicações. 2. ed. São Paulo: 
 Harbra, 1987. 
 
 
BIBLIOGRAFIA BÁSICA 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
PROCEDIMENTOS DE AVALIAÇÃO 
• ARENALES, Selma Helena de Vasconcelos; DAREZZO, 
Artur. Cálculo numérico: aprendizagem com apoio de 
software. São Paulo: Thomson Learning, 2008. 
• RUGGIERO, Marcia A. Gomes; LOPES, Vera Lucia da 
Rocha. Cálculo numérico: aspectos teóricos e 
computacionais. 2. ed. São Paulo: Pearson, 2006. 
 
 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
RESUMO DESTA AULA 
 Identificar e executar as operações 
aritméticas: 
 Escalares; 
 Vetores; 
 Matrizes; 
 Identificar os tipos de funções e seus 
respectivos gráficos; 
 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
VETORES – REPRESENTAÇÃO GRÁFICA 
y 
x 
),( bav

a 
b 
y 
x 
),,( cbav

a 
b 
c 
z jbia  .. v kcjbia

... v
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
OPERAÇÕES COM VETORES. 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
• MULTIPLICAÇÃO POR UM ESCALAR: Seja o vetor v (a,b,c) e 
o escalar real a. O vetor a.v é dado por (a.a, a.b, a.c) 
 
Ex.Se o vetor v é (1,2), o vetor 5.v será (5,10) 
 
 
y 
x 
)2,1(v

1 
2 
y 
x 
)10,5(v

5 
10 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
OPERAÇÕES COM VETORES. 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
• ADIÇÃO: Sejam os vetores v (a,b,c) e u (d,e,f). O vetor 
soma u + v = v + u = (a+d, b+e, c+f). = prop comutativa 
 
Graficamente, temos que: 
 
 v
u

vu


v

u

vu


Polígono Paralelogramo 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
OPERAÇÕES COM MATRIZES. 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
Considere uma tabela com m linhas e n colunas em que 
cada elemento que ocupa a “i-ésima” linha e a “j-ésima” 
coluna é denominado aij 











333231
232221
131211
aaa
aaa
aaa
A
 
 
Matriz com 3 linhas e 3 colunas. Matriz quadrada de ordem 3. 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
OPERAÇÕES COM MATRIZES. 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
Multiplicação de uma matriz A por um escalar real a. Seja a 
matriz Am x n. O produto de a por A, isto é, a.A é igual à 
multiplicação de cada elemento aij por a. 











333231
232221
131211
aaa
aaa
aaa
A











333231
232221
131211
...
...
...
.
aaa
aaa
aaa
A
aaa
aaa
aaa
a
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
OPERAÇÕES COM MATRIZES. 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
Adição de matrizes – para que esteja definida entre duas 
matrizes Am x n e Bp x q é necessário que m = p e n = q. Para 
encontrar a matriz C = A + B, basta adicionar os elementos 
respectivos, isto é, cij = aij + bij 



































333332323131
232322222121
131312121111
333231
232221
131211
333231
232221
131211
bababa
bababa
bababa
bbb
bbb
bbb
aaa
aaa
aaa
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
OPERAÇÕES COM MATRIZES. 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
Produto de matrizes – para que esteja definida entre duas 
matrizes Am x n e Bp x q é necessário que o número de colunas 
da primeira matriz seja igual ao número de linhas da segunda 
matriz, isto é, n = p. A matriz produto terá o número de 
linhas da primeira e o número de colunas da segunda, isto é, 
m linhas e q colunas. Observe: 
A3x4 . B4x5 = C3x5 
 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
OPERAÇÕES COM MATRIZES. 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
Produto de matrizes – Uma vez que esta operação esteja 
definida, cada elemento cij seráformado pela multiplicação 
dos elementos da linha i da matriz A pelos correspondentes 
elementos da coluna j da matriz B. 
Exemplo. 
 
c32 = a31.b12 + a32.b22+a33.b32 
244241
3231
2221
1211
233231
2221
1211
34434241
333231
232221
131211
.
x
x
x
cc
cc
cc
cc
bb
bb
bb
aaa
aaa
aaa
aaa







































AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
OPERAÇÕES COM MATRIZES. 
Dado as matrizes A = (aij)4x3, e B = (bij)3x4, pergunta-se: 
É possível calcular o produto AxB? Qual seria o resultado? 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
OPERAÇÕES COM MATRIZES. 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
PROPRIEDADES: 
• Em regra, Se A.B = B.A diz-se que A e B 
comutam; 
• A.I = A, I – matriz identidade 
• A.(B+C) = A.B + A.C - distributiva à esquerda 
• (B+C).A = B.A + C.A - distributiva à direita 
• A.(B.C) = (A.B).C – associativa 
• A.0 = 0, sendo 0, a matriz nula. 
 
ABBA .. 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
APLICANDO O CONHECIMENTO – EX 1 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
(PETROBRÁS - engenheiro) Sejam os vetores u = (1,2), v = (-2,5) e w = (x,y) do R2. Para que w = 3u – v, devemos ter x + y igual a: 
 
(PETROBRÁS - Engenheiro) 
Sejam os vetores u = (1,2), v = (-2,5) e w = (x,y) do R2. Para 
que w = 3u – v, devemos ter x + y igual a? 
SOLUÇÃO: 
• Multiplicação de um escalar por um vetor: 
 3u = 3.(1,2) = (3,6) 
• Adição/subtração de vetores: 
 3u – v = (3,6) – (-2,5) = (5,1) = w = (x,y) 
• Por comparação, x = 5 e y = 1. Logo, x + y = 6 
 
 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
APLICANDO O CONHECIMENTO – EX 2 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
Considere as seguintes matrizes: 
M = (mij)2x3, N = (nij)axb, P = (pij)cx4, Q = (qij)dxe. 
Para que seja possível determinar M+N, MxP e P-Q, quais 
os valores de a, b, c, d, e ? 
SOLUÇÃO: 
• ADIÇÃO: M2x3 + Naxb  a = 2 e b = 3 
• MULTIPLICAÇÃO: M2x3.Pcx4  c = 3 
• SUBTRAÇÃO: Pcx4 – Qdxe  e = 4 e c =d = 3 
 
 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
FUNCÕES ELEMENTARES 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
Suponha dois conjuntos A e B. Diz-se que f: A  B é uma 
função se para todo elemento x  A existe um único 
elemento y  B. Observe. 
 
 
 
 
 
 
f: x  x+3 
domínio 
• 1 
• 2 
• 3 
• 4 
• 5 
• 6 
• 7 
f 
Imagem 
Contradomínio 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
FUNCÕES ELEMENTARES 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
Graficamente, podemos identificar se uma curva é uma 
função traçando retas verticais. Se as retas cortarem em 
apenas um único ponto a curva, é uma função. 
x 
y 
É função 
y 
Não é função 
x 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
FUNCÕES ELEMENTARES - FUNÇÕES POLINOMIAIS 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
• f(x) = an.x
n + an-1.x
n-1 + an-2.x
n-2 + ...+ a2.x
2 + a1.x +
 a0 
• As raízes reais são os valores de x para os quais y é nulo, 
ou seja, a interseção do gráfico com o eixo x. 
 
 
 
 
 
•P(X) = 4x5 + 2x4 – 6 P(1) = 0 
x 
y 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
TEOREMA DE BOLZANO Ex. (-1,7) 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
Considere um intervalo (a,b) do domínio da função f(x). 
• Se f(a).f(b) > 0, existe um número par de raízes reais no 
intervalo (a,b); 
• Se f(a).f(b) < 0, existe um número ímpar de raízes reais no 
intervalo (a,b). 
 
 
 
 
 x 
y 
a 
b 
f(a) 
f(b) 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
APLICANDO O CONHECIMENTO – EX 3 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
Seja a função polinomial f(x) = 2x3 - 12x2 -3x + 8. Mostre que 
existe ao menos uma raiz real no intervalo (0, 1) da equação 
f(x) = 0. 
SOLUÇÃO: 
 f(0) = 2.(0)3 – 12.(0)2 -3.(0) + 8 = 8 
 f(1) = 2.(1)3 – 12.(1)2 -3.(1) + 8 = -5 
Pelo Teorema de Bolzano, como f(0).f(1) < 0, podemos inferir 
que existe um número ímpar de raízes reais no intervalo 
(0,1). 
 
 
 
 
 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
FUNÇÕES CRESCENTE E DECRESCENTE 
• Se x2 > x1  f(x2) > f(x1) diz-se que a função é 
estritamente crescente; 
• Se x2 > x1  f(x2) < f(x1) diz-se que a função é 
estritamente decrescente. 
 
x 
y 
x1 x2 
f(x2) 
f(x1) 
x 
y 
x1 x2 
f(x1) 
f(x2) 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
FUNÇÕES ELEMENTARES. 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
• Função afim / linear: y = a.x + b 
• Função quadrática: y = a.x2 + b.x + c 
• Função exponencial: y = ax 
• Função Logarítmica: y = Logb(x) 
• Função seno: y = sen(x) 
• Função co-seno: y = cos(x) 
• Função tangente: y = tg(x) 
 
 
 
 
 
 
 
 
 
 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
GRÁFICOS – FUNÇÃO AFIM 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
y 
x 
y 
x 
Crescente / a > 0 Decrescente / a < 0 
b b 
x = -b/a 
x = -b/a 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
GRÁFICOS – FUNÇÃO DO 20 GRAU 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
y 
x 
a > 0 
c 
Raízes reais 
Vértice 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
y 
x 
a < 0 
c 
Raízes reais 
Vértice 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
GRÁFICOS – FUNÇÃO EXPONENCIAL 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
y 
x 
a > 1 
(0,1) 
Crescente 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
y 
x 
0 <a < 1 
Decrescente 
(0,1) 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
GRÁFICOS – FUNÇÃO SENO 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
y 
x 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
1 
-1 
p/2 p 
3p/2 
2p 
Período principal : 2p 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
GRÁFICOS – FUNÇÃO CO-SENO 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
y 
x 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
1 
-1 
p/2 p 
3p/2 
2p 
Período principal : 2p 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
GRÁFICOS – FUNÇÃO TANGENTE 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
y 
x 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
p/2p 3p/2 -p/2 
Período principal : p 
-p 
-3p/2 
AULA 1: INTRODUÇÃO AO PROGRAMA DE COMPUTAÇÃO NUMÉRICA 
CÁLCULO NUMÉRICO 
RESUMINDO 
Nesta aula vocês estudaram: 
 As operações aritméticas: 
 Escalares; 
 Vetores; 
 Matrizes; 
 Os tipos de funções e seus respectivos 
gráficos.