EDITAL DA DISCIPLINA DE1 CALCULO NUM ERICO · EDITAL DA DISCIPLINA DE1 CALCULO NUM ERICO Prof. Dr. Catalunha Atualizado via LATEX em 26 de Outubro de 2015 as 15:03 1Este EDITAL cont

EDITAL DA DISCIPLINA DE1

CALCULO NUMERICO

Prof. Dr. Catalunha

Atualizado via LATEX em 26 de Outubro de 2015 as 15:03

1Este EDITAL contem textos de minha autoria e outros retirados das bibliografias indicadas, outextos correlatos no assunto, sempre que possıvel citadas as fontes. Tais notas nao excluem a consultaao conteudo na integra da bibliografia original citada e sao apenas uma forma de guia de conteudo paraestudos.

Conteudo

I Ementa e Conteudo Programatico 4

II Conteudo da disciplina 7

1 Erros 8

1.1 Exercıcios Complementares gerais do Capıtulo . . . . . . . . . . . . . . . . . . . . 9

2 Sistemas Lineares 10

2.1 Substituicao Retroativa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12

2.2 Metodos Diretos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12

2.2.1 Metodo de Gauss . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12

2.2.2 Metodo de Jordan . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14

2.2.3 Metodo da Pivotacao Completa . . . . . . . . . . . . . . . . . . . . . . . . 14

2.3 Metodos iterativos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14

2.3.1 Metodo de Jacobi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15

2.3.2 Metodo de Gauss-Seidel . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15

2.3.3 Criterios de convergencia . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18


3 Ajuste de Modelo 21


4 Interpolacao 36

4.1 Interpolacao por ajuste de modelo . . . . . . . . . . . . . . . . . . . . . . . . . . . 36

4.2 Interpolacao por lagrange . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38

4.3 Interpolacao por diferenca dividida . . . . . . . . . . . . . . . . . . . . . . . . . . 39


5 Zero de Funcao 44

5.1 Metodo da Bissecao . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44

5.2 Metodo de Pegaso . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44

5.3 Metodo de Newton . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44

5.4 Comparacao dos metodos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 45


2

Prof. Dr. Catalunha - EDITAL Calculo Numerico

6 Integracao 466.1 Metodo de Newton-Cottes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46

6.1.1 Trapezio . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 466.1.2 Simpson I . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 466.1.3 Simpson II . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 466.1.4 Extrapolacao Richardson . . . . . . . . . . . . . . . . . . . . . . . . . . . . 476.1.5 Integral dupla . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 47

6.2 Quadratura Gaussiana . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 476.3 Exercıcios Complementares gerais do Capıtulo . . . . . . . . . . . . . . . . . . . . 47

7 Equacoes Diferenciais Ordinarias 497.1 Metodos de Runge-Kutta . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 497.2 Exercıcios Complementares gerais do Capıtulo . . . . . . . . . . . . . . . . . . . . 49

III Administracao da disciplina 50

8 Dicas de como estudar 518.1 Rotina . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 518.2 Morto . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 518.3 Faminto . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 528.4 Anti-social . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 528.5 Arsenal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 52

9 Procedimentos de avaliacao 539.1 Conteudo e avaliacoes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 539.2 Estrutura do trabalho . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54

9.2.1 Para casa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 549.2.2 Em sala . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54

9.3 Estrutura do Projeto . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 559.4 Estrutura da prova: . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56

9.4.1 Preparacao da sala para a avaliacao . . . . . . . . . . . . . . . . . . . . . . 569.4.2 A prova . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 57

9.5 Outras orientacoes: . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58

10 Elaboracao das Tarefas 6010.1 Tarefa exemplo . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60

10.1.1 Exercıcios Propostos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60

11 Formulario 63

Atualizado via LATEX em 26 de Outubro de 2015 as 15:03 3

Parte I

Ementa e Conteudo Programatico

4


Ementa de disciplina do curso de Engenharia Ambiental conforme Plano Polıtico Pegagogicoaprovado em 2007.

DISCIPLINA:Calculo NumericoCH Total CH Teorica CH Pratica Creditos60 60 0 4

PERIODO: PRE-REQUISITOS:3 Calculo IIOBJETIVO:Capacitar os alunos para a aplicacao de tecnicas de calculo numerico para resolucao de proble-mas de engenharia e utilizacao de ferramentas computacionais auxiliares.

CONTEUDO BASICO:O Metodo numerico. Utilizacao de softwares e hardwares matematicos. Precisao e erro. Sis-temas de equacoes lineares, Ajuste de curvas. Interpolacao. Raızes de equacoes. Integracaonumerica. Solucao numerica de equacoes diferenciais ordinarias.METODOLOGIA DE ENSINO:O ensino sera ministrado de forma expositiva em sala de aula, utilizando os recursos audiovisuaisdisponıveis, com consulta ao material bibliografico e debates sobre o tema. A disciplina seraadministrada utilizando todos os recursos disponıveis no sistema Moodle.

PROCEDIMENTOS DE AVALIACAO:Trabalhos e provas.

BIBLIOGRAFIA BASICA:* Barroso, Leonidas Conceicao et al.. Calculo Numerico (com aplicacoes) / Leonidas ConceicaoBarroso, Magali Maria de Araujo Barroso, Frederico Ferreira Campos, filho, Marcio Luiz Buntede Carvalho, Mirian Lourenco Maia. 2 ed. Sao Paulo. editora HARBRA ltda, 1987.* Ruggiero, Marcia A. Gomes. Calculo numerico: aspectos teoricos e computacionais / MarciaA. Gomes Ruggiero, Vera Lucia da Rocha Lopes. 2 ed.Sao Paulo: MAKRON Books. 1996.* Claudio, Dalcidio Moraes. Calculo numerico computacional, teoria e pratica: Sao Paulo.Atlas, 1989. 464p.BIBLIOGRAFIA COMPLEMENTAR:* Boulos, Paulo; Abud, Zara Issa, Calculo diferencial e Integral, Volume 2, Sao Paulo, PearsonEducation do Brasil, 2002.* Bassanezi, Rodney Carlos. Ensino-aprendizagem com modelagem matematica: uma novaestrategia / Rodney Carlos Bassanezi. 2 ed. Sao Paulo, Contexto, 2004* Sperandio, Decio; Mendes, Joao Teixeira; Silva, Luiz Henry Monken e Silva. Calculo Numericoe Aplicacoes. Sao Paulo. Prentice Hall (Pearson). 2007. 368p.



Na Tabela 1 temos as datas conforme calendario do periodo vigente. Na Tabela 2 o encontroque sera realizado na referida data.

Data15-05-201522-05-201529-05-201505-06-201512-06-201519-06-201526-06-201517-07-201507-08-2015P124-07-201531-07-201514-08-201521-08-201528-08-201504-09-201511-09-201518-09-201525-09-2015P2?-10-201502-10-2015

Tabela 1: Datas

Encontro Conteudo Planejado1 Apresentacao da disciplina. Erros2 Sistema Linear - Gauss, Pivotacao, Jordan.3 Sistema Linear - Jacobi, Newton.4 Ajuste de Modelo - Multiplo, Polinominal5 Ajuste de Modelo - Transformacao6 Interpolacao - Ajuste, Lagrange7 Interpolacao - Diferenca Dividida8 Revisao9 Prova 0110 Zero de Funcao - Bissecao, Newton11 Zero de Funcao - Pegaso12 Integracao - Trapezio, Simpson13 Integracao - Richardson14 Integracao - Dupla, Quadratura15 EDO - Euler, Runge-Kutta16 EDO - Euler, Runge-Kutta17 Revisao18 Prova 0219 Reposicao de Provas20 Prova de Recuperacao

Tabela 2: Conteudo

Legenda: RA: indica data de reposicao de aula. P1 e P2: indicam Prova 1 e Prova 2. PR:indica data de prova de recuperacao. RP: indica reposicao de prova.Observacao: Horario de atendimento aos alunos sera na quinta de 08h as 11h30, Bloco II, Sala10, Ramal 8229 ou Telefone 32328229.

6 Atualizado via LATEX em 26 de Outubro de 2015 as 15:03

Parte II

Conteudo da disciplina

7

Capıtulo 1

Erros

Teoria

Desenvolver o conteudo selecionado da seguinte bibliografia: [Barroso, 1987, Capıtulo 1.1, Pag. 1a 3.]

Desenvolver o conteudo selecionado da seguinte bibliografia: [Freitas, 2000, Capıtulo 2.1, Pag.17.].



Exercıcios Resolvidos

1. Fonte: [Freitas, 2000, Pag. 18, Exemplo do item 2.2.3, adaptado] (exercicios/Sergio018item2.2.3)A diferenca entre um valor exata x e sua aproximacao x foi de x = 100, x = 100.1. Bemcomo y = 0.0006 e y = 0.0004.Pede-se:

(a) Calcule o erro absoluto e relativo no calculo das aproximacoes.Resposta: Concluimos que a aproximacao de x e melhor, pois o erro relativo naaproximacao foi menor.Solucao:Aplicando as equacoes temos

ex =| 100− 100.1 | (1.1)

ex = 0.1 (1.2)

ex =| 0.0006− 0.0004 | (1.3)

ey = 0.0002 (1.4)

Ex =0.1

100.1(1.5)

Ex = 0.000999 (1.6)

Ey =0.0002

0.0004(1.7)

Ey = 0.5 (1.8)

8


Exercıcios Basicos

1.1 Exercıcios Complementares gerais do Capıtulo

Veja lista de exercıcios em cada bibliografia apresentada nas teorias. Alem destes apresento outrosexercıcios conforme lista a seguir.


Capıtulo 2

Sistemas Lineares

Introducao

Desenvolver o conteudo selecionado da seguinte bibliografia: [Barroso, 1987], Capıtulo 2.1, Pag.17 a 19.

Veja Figura 2.1

Figura 2.1: Grafico ilustrativo

Se estas equacoes {y = 1 +xy = 5 −x (2.1)

Forem representadas na forma de um sistema linear, como em:{y −x = 1y +x = 5

(2.2)

10


Ou {x1 −x2 = 1x1 +x2 = 5

(2.3)

A solucao seria:

1 octave :1> A=[1 ,−1;1 ,1 ]2 A =34 1 −15 1 167 octave :2> b = [ 1 ; 5 ]8 b =9

10 111 51213 octave :3> A∗∗−1∗b14 ans =1516 317 21819 octave :4>

Caso multipliquemos uma das equacoes por um numero a solucao do sistema nao se altera.Veja Figura 2.2 e equacao multiplicada.{

y −x = 110y +10x = 50

(2.4)

Outro exemplo usando

1 # c o n s i d e r e x =[1 ,2 ,3 ] e y =[1 ,2 ,1 ]2 >>> import numpy as np3 >>> np . matrix ( [ [ 3 , 6 , 1 4 ] , [ 6 , 1 4 , 3 6 ] , [ 1 4 , 3 6 , 9 8 ] ] )4 matrix ( [ [ 3 , 6 , 1 4 ] ,5 [ 6 , 14 , 3 6 ] ,6 [ 1 4 , 36 , 9 8 ] ] )7 >>> X=np . matrix ( [ [ 3 , 6 , 1 4 ] , [ 6 , 1 4 , 3 6 ] , [ 1 4 , 3 6 , 9 8 ] ] )8 >>> Y=np . matrix ( [ [ 4 ] , [ 8 ] , [ 1 8 ] ] )9 >>> Y

10 matrix ( [ [ 4 ] ,11 [ 8 ] ,12 [ 1 8 ] ] )13 >>> X∗∗−1∗Y14 matrix ( [ [ − 2 . ] ,15 [ 4 . ] ,16 [ − 1 . ] ] )



2.1 Substituicao Retroativa

Teoria

Desenvolver o conteudo selecionado da seguinte bibliografia: [Barroso, 1987], Capıtulo 2.1 a 2.1.6,Pag. 20 a 27


Exercıcios Basicos

2.2 Metodos Diretos

Teoria

Desenvolver o conteudo selecionado da seguinte bibliografia: [Barroso, 1987], Capıtulo 2.2, Pag.27 a 49. Abrangendo os itens:

Figura 2.2: Grafico ilustrativo




Exercıcios Basicos

2.2.1 Metodo de Gauss

Teoria


1. Com (n-1) passos um Sistema Linear ax=b e transformado num sistema triangular equi-valente, o qual se resolve facilmente por substituicao. Um destes metodos e o de Gauss.Conforme descrito neste exemplo.

Considere o seguinte sistema linear de equacoes:2x1 +3x2 −x3 = 54x1 +4x2 −3x3 = 32x1 −3x2 +x3 = −1

(2.5)

Escrevendo este sistema linear numa forma de matriz aumentada, M, teremos:

M =

2 3 −1 | 54 4 −3 | 32 −3 1 | −1

=(A | b

)(2.6)

A escolha do pivo sera sempre aii de Mk. Os multiplicadores sao mki = −

akijpakipjp

. Em que

ip e a linha pivotal; jp e a coluna pivotal; k e a etapa. Sendo as linhas calculadas comoLk+1i = Lki + Lkip ∗m

ki e sucessivamente teremos:

Etapa 0

Como a011 = 2, os multiplicadores serao m02 = −4

2= −2 e m0

3 = −2

2= −1 temos as linhas

atualizadas:

L11 = L0

1 (2.7)

L12 = L0

2 + L01 ∗m0

2 (2.8)

L13 = L0

3 + L01 ∗m0

3 (2.9)

(2.10)

Aplicando os valores nas equacoes temos:

M1 =

2 3 −1 | 54 + (2 ∗ −2) 4 + (3 ∗ −2) −3 + (−1 ∗ −2) | 3 + (5 ∗ −2)2 + (2 ∗ −1) −3 + (3 ∗ −1) 1 + (−1 ∗ −1) | −1 + (5 ∗ −1)

(2.11)

Calculando teremos:

M1 =

2 3 −1 | 50 −2 −1 | −70 −6 2 | −6

(2.12)



Etapa 1

Como a122 = −2, os multiplicadores serao m13 = −−6

−2= −3 e as linhas atualizadas serao:

L21 = L0

1 (2.13)

L22 = L1

2 (2.14)

L23 = L1

3 + L12 ∗m1

3 (2.15)

(2.16)

Aplicando os valores nas equacoes temos:

M1 =

2 3 −1 | 50 −2 −1 | −70 −6 + (−2 ∗ −3) 2 + (−1 ∗ −3) | −6 + (−7 ∗ −3)

(2.17)

Calculando teremos:

M1 =

2 3 −1 | 50 −2 −1 | −70 0 5 | 15

(2.18)

Retornando a matriz para forma de sistema linear teremos o sistema triangular desejado.

2x1 +3x2 −x3 = 5

−2x2 −1x3 = −75x3 = 15

(2.19)

Basta resolver o sistema pelo metodo de substituicao retroativa.

x1 =5− (3 ∗ 2− 1 ∗ 3)

2= 1

x2 =−7− (−1 ∗ 3)

−2= 2

x3 =15

5= 3

(2.20)

Resultando em x1 = 1, x2 = 2, x3 = 3.



Exercıcios Basicos

2.2.2 Metodo de Jordan

Teoria


Exercıcios Basicos

2.2.3 Metodo da Pivotacao Completa

Teoria


Exercıcios Basicos

2.3 Metodos iterativos

Teoria

Desenvolver o conteudo selecionado da seguinte bibliografia: [Barroso, 1987], Capıtulo 2.3, Pag.49 a 72. [Buffoni, 2000], Pag 20 a 23. Abrangendo os itens:

2.3.1 Metodo de Jacobi

Teoria


Exercıcios Basicos

2.3.2 Metodo de Gauss-Seidel

Teoria


1. Considere o seguinte sistema linear de equacoes:2x1 +3x2 −x3 = 54x1 +4x2 −3x3 = 32x1 −3x2 +x3 = −1

(2.21)

Isolando as equacoes teremos: Considere o seguinte sistema linear de equacoes:x1 =

5− (3x2 − x3)2

x2 =3− (4x1 − 3x3)

4

x3 =−1− (2x1 − 3x2)

1

(2.22)



Considere que os valores inciais de iteracao sao:x1 = 0x2 = 0x3 = 0

(2.23)

Aplicando o processo iterativo teremos:k x1 x2 x3

0 0 0 0

15 − (3 ∗ 0 − 0)

2=

5

2

3 − (4 ∗ 52

− 3 ∗ 0)

4=

−7

4

−1 − (2 ∗ 52

− 3 ∗ −74

)

1= −11.25

2

5 − (3 ∗−7

4− (−11.25))

2= −0.5

3 − (4 ∗ −0.5 − 3 ∗ −11.25)

4= −7.1875

−1 − (2 ∗ −0.5 − 3 ∗ −7.1875)

1= −21.562

3 ... ... ...

(2.24)

Veja as equacoes do octave para ajudar nos calculos.

1 func t i on s l g a u s s s e i d e l ( )2 f1=@( x2 , x3 ) (5−(3∗x2−x3 ) ) /2 ;3 f 2=@( x1 , x3 ) (3−(4∗x1−3∗x3 ) ) /4 ;4 f 3=@( x1 , x2 ) (−1−(2∗x1−3∗x2 ) ) /1 ;5 x1=0;6 x2=0;7 x3=0;8 [ 0 , x1 , x2 , x3 ]9 f o r i =1:10

10 x1=f1 ( x2 , x3 ) ;11 x2=f2 ( x1 , x3 ) ;12 x3=f3 ( x1 , x2 ) ;13 [ i , x1 , x2 , x3 ]14 endfor15 endfunct ion

2. Considere o seguinte sistema linear de equacoes:x1 +4x2 +52x3 = 57

27x1 +110x2 −3x3 = 13422x1 +2x2 +14x3 = 38

(2.25)

Veja que desta forma o sistema nao converge entao precisamos reordenar as equacoes:22x1 +2x2 +14x3 = 3827x1 +110x2 −3x3 = 134x1 +4x2 +52x3 = 57

(2.26)

Isolando as equacoes teremos: Considere o seguinte sistema linear de equacoes:x1 =

38− (2x2 + 14x3)

22

x2 =134− (27x1 − 3x3)

110

x3 =57− (x1 + 4x2)

52

(2.27)



Considere que os valores inciais de iteracao sao:x1 = 0x2 = 0x3 = 0

(2.28)

Aplicando o processo iterativo teremos:

k x1 x2 x3

0 0 0 0

138 − (2 ∗ 0 + 14 ∗ 0)

22= 1.7272

134 − (27 ∗ 1.7272 − 3 ∗ 0)

110= 0.7942

57 − (1.7272 + 4 ∗ 0.7942)

52= 1.0018

238 − (2 ∗ 0.7942 + 14 ∗ 1.0018)

22= 1.0175

134 − (27 ∗ 1.0175 − 3 ∗ 1.0018)

110= 0.9957

57 − (1.0175 + 4 ∗ 0.9957)

52= 0.9999

338 − (2 ∗ 0.9957 + 14 ∗ 0.9999)

22= 1.0004

134 − (27 ∗ 1.0004 − 3 ∗ 0.9999)

110= 0.9999

57 − (1.0004 + 4 ∗ 0.9999)

52= 1.0000

3 ... ... ...

(2.29)

Veja as equacoes do octave para ajudar nos calculos.

1 octave :7> A=[22 2 14 ;27 110 −3;1 4 52 ]2 A =34 22 2 145 27 110 −36 1 4 527 octave :10> b = [ 3 8 ; 1 3 4 ; 5 7 ]8 b =9

10 3811 13412 5713 octave :11> A∗∗−1∗b14 ans =15 1.0000016 1.0000017 1.000001819 func t i on s l g a u s s s e i d e l 2 ( )20 f1=@( x2 , x3 ) (38−(2∗x2+14∗x3 ) ) /22 ;21 f2=@( x1 , x3 ) (134−(27∗x1−3∗x3 ) ) /110 ;22 f3=@( x1 , x2 ) (57−(x1+4∗x2 ) ) /52 ;23 x1=0;24 x2=0;25 x3=0;26 [ 0 , x1 , x2 , x3 ]27 f o r i =1:1028 x1=f1 ( x2 , x3 ) ;29 x2=f2 ( x1 , x3 ) ;30 x3=f3 ( x1 , x2 ) ;



31 [ i , x1 , x2 , x3 ]32 endfor33 endfunct ion343536 octave :12> s l g a u s s s e i d e l 237 ans =3839 0 0 0 04041 ans =4243 1.00000 1.72727 0.79421 1.001844445 ans =4647 2.00000 1.01753 0.99575 0.999994849 ans =5051 3.00000 1.00039 0.99990 1.000005253 ans =5455 4.00000 1.00001 1.00000 1.000005657 ans =5859 5.00000 1.00000 1.00000 1.00000

Exercıcios Basicos

2.3.3 Criterios de convergencia

Teoria


Exercıcios Basicos



1. Fonte: SistemaLinear99 (exercicios/SistemaLinear99) Resolva cada um dos sistemas deequacoes lineares dados a seguir (I,II,...) conforme pede-se:



I

3x1 4x2 −5x3 x4 = −10x2 x3 −2x4 = −1

4x3 −5x4 = 32x4 = 2

II

2x1 −1x2 3x3 5x4 = −76x1 −3x2 12x3 11x4 = 44x1 −1x2 10x3 8x4 = 4

−2x2 −8x3 10x4 = −60

III

x1 4x2 52x3 = 5727x1 110x2 −3x3 = 13422x1 2x2 14x3 = 38

IV

3.2x1 x2 2x3 = 8.2−1x1 1.5x3 −2.4x4 = 2.844.1x1 2.5x2 2x4 = 13.6x1 2.8x4 = 4.72

V5x y z = 53x 4y z = 63x 3y 6z = 0

VI−4x1 10x2 = 195x1 3x2 = 15

VII3x1 0 x3 = 3x1 −x2 0 = 13x1 x2 2x3 = 9

VIIIx1 0 0 = 12x1 5x2 0 = 23x1 6x2 4x3 = 3

IXx1 = 1x1 x2 = −1x1 x2 x3 = 3x1 x2 x3 x4 = 3



Xx1 x2 x3 x4 = 4x1 x2 x3 = 3x1 x2 = 2x1 = 1

XIx1 −3x2 x3 = 6

4x2 −x3 = 3x3 = 4

XIIx1 x2 2x3 = 42x1 −x2 −x3 = 0x1 −x2 −x3 = −1

XIII10x1 2x2 6x3 = 28x1 10x2 9x3 = 72x1 −7x2 −10x3 = −17

XIV−4x1 10x2 = 195x1 3x2 = 15

.Pede-se:

(a) Metodo de Gauss.

1 ##################################################################2 (%i 2 ) a l g s y s ( [ 5∗ x+y+z =5,3∗x+4∗y+z =6,3∗x+3∗y+6∗z =0 ] , [ x , y , z ] ) ;3 (%o2 ) [ [ x = 1 , y = 1 , z = − 1 ] ]4 ##################################################################5 octave :2> a=[3 0 1 ;1 −1 0 ; 3 1 2 ]6 a =78 3 0 19 1 −1 0

10 3 1 21112 octave :3> b = [ 3 ; 1 ; 9 ]13 b =1415 316 117 91819 octave :4> a∗∗−1∗b20 ans =2122 −2.0000



23 −3.000024 9 .000025 ##################################################################

(b) Pivotacao Completa

(c) Jordan

(d) Jacobi, considerando o vetor inicial x(0) = [0, 0, 0, 0]T com k ≥ 3 ou Erro ≤ 10−2

(e) Gauss-Seidel, considerando o vetor inicial x(0) = [0, 0, 0, 0]T com k ≥ 3 ou Erro ≤10−2


Capıtulo 3

Ajuste de Modelo

Teoria

Uma regressao ou curva de tendencia pode ser o primeiro passo para uma modelagem. Umarelacao funcional, obtida atraves de uma ajuste dos dados, propicia condicoes para a elaboracaode hipoteses que levam a formulacao dos modelos.

Os modelos sao relacoes funcionais que incorporam as particularidades do fenomeno analisado.Um reta ou curva ajustada nao pode ser considerado um modelo matematico para uma deter-

minada situacao. Neste caso, a reta ou curva simplesmente descreve uma tendencia dos fatos nointervalo pesquisado.

Mesmo que uma curva possa fazer alguma previsao de futuros valores para o fenomeno estu-dado, ainda assim, tal formulacao nao poderia ser considerado um modelo matematico do fenomenoenquanto seus parametros nao tiverem algum significado biologico, quımico, fısico ou probabilısticocom o fenomeno em estudo !

Contudo o processo de ajuste de curvas e um dos mais importantes passos para o treinamentoem modelagem e entendimento do comportamento de um fenomeno que nao seja aleatorio.

Desenvolver conteudo da seguinte bibliografia: [Barroso, 1987, Paginas 323 ate 356], e [Rodney, 2006,Paginas 54 ate 85].

Uma seguencia de etapas importante para qualquer ajuste de modelo pode ser vista a seguir,sendo que estas etapas devem ser seguidas nos exercicios deste item.

1. O grafico dos pontos originais.

2. O modelo transformado e seus coeficientes.

3. A representacao da modelo transformado e pontos transformados num grafico.

4. Os coeficientes da equacao de regressao do modelo transformado

5. O coeficiente de ajustamento dos dados transformados.

6. O modelo original e seus coeficientes, obtido da relacao com os coeficientes do modelo trans-formado.

7. A representacao da modelo original e pontos originais num grafico.

8. O coeficiente de ajustamento dos dados originais.

22


9. O ajustamento por processo nao-linear, R ou Gnuplot.

10. O coeficiente de ajustamento pelo processo nao-linear.

11. A representacao num grafico dos pontos originais, modelo por mınimos quadrados e modelopor ajuste nao-linear.


Em cada situacao a seguir aplicando as etapas anteriores ajuste um modelo aos dados:

1. Fonte: [Barroso, 1987, Pag. 340, Exercıcio 7.7]. (exercicios/Barroso340Exc7.7) Considereos pontos da Tabela 3.1.

Tabela 3.1: Dados de Campox 0.1 1.5 3.3 4.5 5.0y 5.9 8.8 12 19.8 21.5

Pede-se:

(a) Ajustar os pontos ao modelo y = aebx.Resposta: Resposta: y = 5.7014e0.26283x

Solucao:E sempre bom plotar o grafico dos pontos para ter uma nocao nas manipulacoes ma-tematicas, veja Figura 3.1.

Figura 3.1: Modelo y = abx



Efetuando as transformacoes teremos:

y = aebx (3.1)

ln y = ln a+ ln ebx (3.2)

ln y = ln a+ bx ln e (3.3)

yt = at + bx (3.4)

Ficando os pontos na forma, veja Tabela3.2. Como os pontos originais foram alterados,

Tabela 3.2: Pontos Transformadosx 0.1 1.5 3.3 4.5 5.0y 5.9 8.8 12 19.8 21.5yt 1.7750 2.1748 2.4849 2.9857 3.0681

e importante ver os pontos transformados, veja Figura 3.2.

Figura 3.2: Modelo yt = at+ bx

A somas para montagem da matriz normal serao:

1 octave :1> x = [ 0 . 1 , 1 . 5 , 3 . 3 , 4 . 5 , 5 . 0 ]2 x =3 0.10000 1.50000 3.30000 4.50000 5.000004 octave :2> y =[5 .9 , 8 . 8 , 12 , 19 . 8 , 2 1 . 5 ]5 y =6 5.9000 8 .8000 12.0000 19.8000 21.50007 octave :3> yt=log ( y )8 yt =9 1.7750 2 .1748 2 .4849 2 .9857 3 .0681



10 octave :4> mx=[ l ength ( x ) ,sum( x ) ; sum( x ) ,sum( x . ˆ 2 ) ]11 mx =12 5.0000 14.400013 14.4000 58.400014 octave :5> my=[sum( yt ) ; sum( x .∗ yt ) ]15 my =16 12.48817 40 .416

Podemos ver estas somas na forma de planilha, como a seguir, Tabela 3.3,

Tabela 3.3: Pontos Transformadosi x y yt x2 xyt (y − y)2 y2

1 0.10000 5.9000 1.7750 0.010000 0.17750 0.0021868 348102 1.50000 8.8000 2.1748 2.250000 3.26213 0.1178681 774403 3.30000 12.0000 2.4849 10.890000 8.20019 2.4729728 1440004 4.50000 19.8000 2.9857 20.250000 13.43557 1.4272960 3920405 5.00000 21.5000 3.0681 25.000000 15.34026 0.0793819 462250∑

14.4000∑

68∑

12.488∑

58.4000∑

40.416∑

4, 0997∑

1110.5

Ficando a matriz normal na forma:[5.0000 14.400014.4000 58.4000

]∗[atb

]=

[12.48840.416

](3.5)

Este sistema linear deve ser resolvido conforme capıtulo 2, Sistemas Lineares, por qual-quer metodo pedido.

Antes de continuar e importante saber que a matriz normal, anterior, pode ser obtidatambem por manipulacao matricial.

X =

1 x11 x12 · · · x1p1 x21 x22 · · · x2p...

......

. . ....

1 xn1 xn2 · · · xnp

β =

β0β1...βp

Y =

y1y2...yn

(3.6)

Matriz de soma x = X ′X (3.7)

Matriz de soma y = X ′Y (3.8)

Como a seguir.

1 octave :13> mx= [ [ 1 ; 1 ; 1 ; 1 ; 1 ] , x ’ ]2 mx =3 1.00000 0.100004 1.00000 1.500005 1.00000 3.300006 1.00000 4.500007 1.00000 5.000008 octave :15> my=yt ’



9 my =10 1.775011 2 .174812 2 .484913 2 .985714 3 .068115 octave :16> (mx’∗mx)16 ans =17 5.0000 14.400018 14.4000 58.400019 octave :17> mx’∗my20 ans =21 12 .48822 40 .416

A solucao do sistema linear foi resolvido de forma matricial, para fins de simplificacao.Ficando:

1 octave :6> mc=mx∗∗−1∗my2 mc =3 1.740714 0.26283

Potanto mc e a matriz de coeficiente do modelo transformado, e como as relacoes paraa sao:

at = ln a (3.9)

eat = eln a (3.10)

eat = a (3.11)

(3.12)

Precisamos aplica exp() para obter os valores originais, ficando o b inalterado.

1 octave :8> a=exp (mc(1 ) )2 a = 5.70143 octave :9> b=mc(2)4 b = 0.26283

Entao o valor dos coeficientes originais sao a = 5.7014 e b = 0.26283 ficando o modeloy = 5.7014e0.26283x conforme Figura 3.3

(b) Calcule o grau de ajustamento r2 do modelo aos dados.Resposta: Resposta: 97.8% ou 0.97793.

Solucao:O grau de ajustamento mede o grau de explicacao do modelo aos pontos originais. Or2 tem o seguinte modelo em uma de suas formas:

R2 = 1−∑

(yi − yi)2∑(y2i )−

(∑yi)2

n

(3.13)



Uma das informacoes importantes para calculo do r2 e o valor de yi, veja Tabela 3.4

1 octave :20> ye =5.7014∗ e . ˆ ( 0 . 2 6 2 8 3 . ∗ x )2 ye =34 5.8532 8 .4567 13.5726 18.6053 21.2183

Tabela 3.4: Pontos Estimados pelo modelox 0.1 1.5 3.3 4.5 5.0y 5.9 8.8 12 19.8 21.5yi 5.8532 8.4567 13.5726 18.6053 21.2183

Os valores intermediarios pode ser vistos na Tabela 3.3. Aplicando no modelo r2 temos:

1 octave :21> 1−(sum ( ( y−ye ) . ˆ 2 ) ) /(sum( y . ˆ 2 )−(sum( y ) ˆ2) /( l ength ( x ) ) )2 ans = 0.97793

Pontanto o modelo explica 97% dos dados originais.

(c) Utilizando o gnuplot ajuste o modelo aos dados.Resposta: Resposta: .

Solucao:Para ajustar os dados a um modelo utilizando o gnuplot basta criar o arquivo de pontos,conforme a seguir:

Arquivo 3.1: exercicios/Barroso340Exc7.7/grafico.pts

Figura 3.3: Modelo y = aebx



1 #x ; y ;2 0 .1 5 .93 1 .5 8 .84 3 .3 125 4 .5 19 .86 5 21 .5

O script do gnuplot para ajustamento pode ser como a seguir:

Arquivo 3.2: exercicios/Barroso340Exc7.7/graficoC1.txt

1 # Comentario2 reset3 set term pngca i ro4 set output ’ g ra f i coC1 . png ’5 set grid6 set key out s id e c en te r bottom t i t l e ’ Legenda : ’7 set t i t l e ”Modelo Gnuplot”8 set xlabel ”x”9 set ylabel ”y”

10 f ( x )=a∗exp(b∗x )11 f i t f ( x ) ” g r a f i c o . pts ” using ( $1 ) : ( $2 ) v ia a , b12 plot ” g r a f i c o . pts ” using ( $1 ) : ( $2 ) t i t l e ”Pontos O r i g i n a i s ” ,

\13 f ( x ) t i t l e ”Modelo Gnuplot”

Apresentando o grafico conforme Figura 3.4.

Figura 3.4: Ajustando modelo a dados com Gnuplot

O gnuplot tambem fornece uma arquivo, chamado ”fit.log”de resultados do ajustamento



com alguma analise estatıstica. Os parametros do modelo sao f(x) = 5.41218e0.277491x

conforme segue no arquivo.

Arquivo 3.3: exercicios/Barroso340Exc7.7/fit.log

123 ∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗∗

4 Mon Feb 4 0 9 : 1 3 : 02 2013567 FIT : data read from ” g r a f i c o . pts ” us ing ( $1 ) : ( $2 )8 format = x : z9 #datapo int s = 5

10 r e s i d u a l s are weighted equa l l y ( un i t weight )1112 func t i on used f o r f i t t i n g : f ( x )13 f i t t e d parameters i n i t i a l i z e d with cur rent v a r i a b l e va lue s14151617 I t e r a t i o n 018 WSSR : 21307.4 d e l t a (WSSR) /WSSR : 019 d e l t a (WSSR) : 0 l i m i t f o r s topping : 1e−0520 lambda : 275 .1752122 i n i t i a l s e t o f f r e e parameter va lue s2324 a = 125 b = 12627 After 8 i t e r a t i o n s the f i t converged .28 f i n a l sum of squares o f r e s i d u a l s : 3 .6864929 r e l . change during l a s t i t e r a t i o n : −5.55842e−073031 degree s o f freedom (FIT NDF) : 332 rms o f r e s i d u a l s (FIT STDFIT) = s q r t (WSSR/ ndf ) : 1 .1085333 var iance o f r e s i d u a l s ( reduced ch i square ) = WSSR/ ndf : 1 .228833435 Fina l s e t o f parameters Asymptotic Standard Error36 ======================= ==========================3738 a = 5.41218 +/− 0 .6779 (12.53%)39 b = 0.277491 +/− 0.02867 (10.33%)404142 c o r r e l a t i o n matrix o f the f i t parameters :4344 a b



45 a 1 .00046 b −0.964 1 .000

1 >>> x=np . matrix ( [ 0 . 1 , 1 . 5 , 3 . 3 , 4 . 5 , 5 ] )2 >>> y=np . matrix ( [ 5 . 9 , 8 . 8 , 1 2 . 0 , 1 9 . 8 , 2 1 . 5 ] )3 >>> v1=np . matrix (np . ones (5 ) )4 >>> X=np . concatenate ( ( v1 .T, x .T) ,1 )5 >>> Y=y .T6 >>> Y=np . l og (Y)7 >>> (X.T∗X) ∗∗−1∗(X.T∗Y)8 matrix ( [ [ 1 . 74071419 ] ,9 [ 0 . 2 6 2 8 3 1 5 7 ] ] )

Exercıcios Basicos



Em cada situacao a seguir aplicando as etapas anteriores ajuste um modelo aos dados:

1. FAHL et al, 1982, estudando as caracterısticas fisiologicas de tres cultivares de mandioca emque x=massa seca de raızes de mandioca, ; y=massa seca total da planta, . Tabela 3.5 .

Tabela 3.5: Dados para analisex y x y x y

478,8 116,1 375 63,1 492,5 168,1695,4 231,9 565,1 192,9 676,7 312,21095,1 462,2 816,7 338,8 908,1 479,61301,2 615,7 965 504,5 1141,8 600,61463,1 764,8 1120,4 589,8 1304,3 786,5

Pede-se:

(a) Ajuste os dados da Tabela 3.5 ao modelo y = a+ bx.

2. PEREIRA et al, 1982, estudando as relacoes radiometricas em tres cultivares de mandioca,em que x = radiacao refletida do infravermelho proximo, ; y=radiacao incidente do infra-vermelho proximo,. Tabela 3.6 .Pede-se:

(a) Ajuste os dados da Tabela 3.6 ao modelo y = a+ bx.

3. QUAGGIO et al, 1985, estudando metodos de determinacao da necessidade de calagem, emque x = (H+Al), meq/100 ; y=potencial hidrogenionico,. Tabela 3.7 .Pede-se:



Tabela 3.6: Dados para analise

x y x y x y237 77 461 161 292 105286 98 489 161 223 77335 126 482 175 161 63384 147 432 154419 154 391 133


x y x y x y x y6,43 2,8 5,96 6,2 6,79 1,4 6,37 2,86,66 2 5,58 5,5 6,42 2,9 5,59 5,95,7 8 4,93 13,6 6,7 1,9 4,91 11,66,66 2,3 5,24 8,4 6,05 3,8 4,86 17,55,63 7,2 4,76 11,8 5,61 6,3 4,62 21,55,85 5 4,22 28,9 6,05 4,54,89 11,4 6,43 2,4 5,79 5,3

(a) Ajuste os dados da Tabela 3.7 ao modelo y = axb.

4. CHANG et al, 1982, estudando a resistencia do corte do colmo de cana-de-acucar, em quex = energia armazenada,kgf.cm; y=angulo do pendulo, graus.. Tabela 3.8 .


x y x y5,5 3 55 4020 8 60 4329 13 67 6339 17 80 8644 27

Pede-se:

(a) Ajuste os dados da Tabela 3.8 ao modelo y = abx.

5. PAES DE CAMARGO et al, 1982, estudando a construcao de um tensiometro simples deleitura direta, em que x = tensao,mb; y=altura da camara, mm.. Tabela 3.9 .Pede-se:

(a) Ajuste os dados da Tabela 3.9 ao modelo y =ax

b+ x. A solucao do sistema linear deve

ser por Gauss.




x 9 12 30 42 57 102 147 210 290y 217 291 439 515 603 681 716 746 755

(b) Na solucao do sistema linear do item anterior utilize o metodo da pivotacao completae depois comprove os resultados com um bom metodo iterativo.

(c) Obtenha o valor da altura da camara para uma tensao igual a sua idade. Use inter-polacao utilizando diferenca dividida ou lagrange, para um polinomio interpolador degrau 3.

6. QUAGGIO et al, 1985, estudando calagem para a sucessao batata-triticale-milho, em que y= porcentagem de tuberculos graudos,mb; x=calcio no solo, meq/100. Tabela 3.10 .


x y x y x y0,1 65,9 0,2 68,4 0,3 64,91 92,7 0,7 89,1 1,1 87

1,6 93 0,8 93,2 1,2 89,62,8 94,8 1,4 95,5 1,6 93,23,3 95,5 1,9 94,8 1,8 91,7

Pede-se:

(a) Ajuste os dados da Tabela 3.10 ao modelo y = ab

x.

7. QUAGGIO et al, 1985, estudando metodos de determinacao da necessidade de calagem, emque x = (H+Al), meq/100 ; y=potencial hidrogenionico, . Tabela 3.11 .


x y x y x y x y6,43 2,8 5,96 6,2 6,79 1,4 6,37 2,86,66 2 5,58 5,5 6,42 2,9 5,59 5,95,7 8 4,93 13,6 6,7 1,9 4,91 11,66,66 2,3 5,24 8,4 6,05 3,8 4,86 17,55,63 7,2 4,76 11,8 5,61 6,3 4,62 21,55,85 5 4,22 28,9 6,05 4,54,89 11,4 6,43 2,4 5,79 5,3

Pede-se:



(a) Ajuste os dados da Tabela 3.11 ao modelo y = ea+bx.

8. FAHL et al, 1982, estudando as caracterısticas fisiologicas de tres cultivares de mandioca,em que y = ındice de area foliar; x=dias apos o plantio.. Tabela 3.12 .


x y60 0,5182 1,41116 2,97143 3,54172 3,35214 1,83249 0,48

Pede-se:

(a) Ajuste os dados da Tabela 3.12 ao modelo y = a+ bx+ cx2.

9. Fonte: Ajuste09 (exercicios/Ajuste09) MACHADO et al, 1982, estudando as caracterısticasfisiologicas de quadro variedades de milho, em que y = ındice de area foliar; x=massa secada planta, . Tabela 3.13 .


x y24,6 0,45128,4 1,46292,5 2,5624,3 2,99857,5 2,731041,2 2,11191,4 1,5

Pede-se:

(a) Ajuste os dados da Tabela 3.13 ao modelo y = a + b√x + cx. Use o processo dos

Mınimos Quadrados.Resposta: y = −1.3965268 + 0.3729929 ∗

√(x)− 0.0081978 ∗ x

Solucao:

1 octave :33> x2 x =3



4 24 .600 128.400 292.500 624.300 857.500 1041.2001191.400

56 octave :34> y7 y =89 0.45000 1.46000 2.50000 2.99000 2.73000 2.10000

1.50000101112 octave :50> x1=s q r t ( x )13 x1 =1415 4.9598 11.3314 17.1026 24.9860 29.2831 32.2676

34.51671617 octave :51> x2=x18 x2 =1920 24 .600 128.400 292.500 624.300 857.500 1041.200

1191.4002122 octave :52> X=[ ones (7 , 1 ) , x1 ’ , x2 ’ ]23 X =2425 1.0000 4 .9598 24.600026 1 .0000 11.3314 128.400027 1 .0000 17.1026 292.500028 1 .0000 24.9860 624.300029 1 .0000 29.2831 857.500030 1 .0000 32.2676 1041.200031 1 .0000 34.5167 1191.4000323334 octave :40> Y=[y ’ ]35 Y =3637 0.4500038 1.4600039 2.5000040 2.9900041 2.7300042 2.1000043 1.5000044 octave :54> X’∗X45 ans =4647 7.0000 e+00 1.5445 e+02 4.1599 e+0348 1.5445 e+02 4.1599 e+03 1.2201 e+05



49 4 .1599 e+03 1.2201 e+05 3.7312 e+065051 octave :55> X’∗Y52 ans =5354 13 .73055 335.72056 9111.0365758 octave :42> (X’∗X) ∗∗−1∗(X’∗Y)59 ans =6061 −1.396526862 0.372992963 −0.00819786465 octave :45> f=@( x ) −1.3965268+0.3729929∗ s q r t ( x ) −0.0081978∗x66 f =6768 @( x ) −1.3965268 + 0.3729929 ∗ s q r t ( x ) − 0.0081978 ∗ x6970 octave :46> f ( 2 4 . 6 )71 ans = 0.2517972 octave :57> f ( 1 2 8 . 4 )73 ans = 1.777474 octave :60> f ( x )75 ans =7677 0.25179 1.77740 2.58478 2.80519 2.49625 2.10352

1.71108

10. MACHADO et al, 1982, estudando as caracterısticas fisiologicas de quadro variedades demilho, em que x = dias apos o plantio; y=massa seca da planta, . Tabela 3.14 .


x y36 245,950 526,264 1102,978 1346,792 1619,8106 1775,8120 1775

Pede-se:



(a) Ajuste os dados da Tabela 3.14 ao modelo y = ea+bx+cx2.

11. Considere os dados da Tabela 3.15 .


x y22 1536 114,150 292,864 665,278 826,292 843,6106 782,5120 722

Pede-se:

(a) Ajuste os dados da Tabela 3.15 ao modelo y = αxβe−γx.


Capıtulo 4

Interpolacao

Introducao

Desenvolver o conteudo selecionado da seguinte bibliografia: [Barroso, 1987], Pag. 151 a 153.

4.1 Interpolacao por ajuste de modelo

Teoria

Desenvolver o conteudo selecionado da seguinte bibliografia: [Barroso, 1987], Capıtulo 4.3, Pag.153-155, 159.

Desenvolver o conteudo selecionado da seguinte bibliografia: [Barroso, 1987], Capıtulo 4.4,Pag. 159 a 161.

Observacoes basicas na aplicacao do metodo:

1. Usa um Polinomio Interpolador, Pg de qualquer grau n >= g >= n sendo a quantidade depontos n.

2. Resulta em uma equacao.


1. Fonte: [Barroso, 1987, Pag. 167, Exemplo 4.9]. (exercicios/Barroso167Exp4.9) Considere ospontos da Tabela 10.1.

Tabela 4.1: Dados de Campoi 0 1 2 3x 0.0 0.2 0.4 0.5y 0.000 2.008 4.064 5.125

Pede-se:

(a) Utilizando ajuste de modelo determine o valor de x = 0.3 por um polinomio de 3 grau.Resposta: Resposta: 3.027

Solucao:

37


O interpolador de lagrange fornece o mesmo resultado que um polinomio de mesmograu para o valor desejado. Provando iremos estimar este polinomio usando o octave,para montar as matrizes normais e calcular a resolucao do sistema, mas voce deveresolver o sistema por um processo tipo Gauss,etc .

1 octave :1> x = [ 0 , 0 . 2 , 0 . 4 , 0 . 5 ]2 x =3 0.00000 0.20000 0.40000 0.500004 octave :2> y = [ 0 , 2 . 0 0 8 , 4 . 0 6 4 , 5 . 1 2 5 ]5 y =6 0.00000 2.00800 4.06400 5.125007 octave :7> mx=[ l ength ( x ) ,sum( x ) ,sum( x . ˆ 2 ) ,sum( x . ˆ 3 ) ; sum( x ) ,sum( x . ˆ 2 )

,8 sum( x . ˆ 3 ) ,sum( x . ˆ 4 ) ; sum( x . ˆ 2 ) ,sum( x . ˆ 3 ) ,sum( x . ˆ 4 ) ,sum( x . ˆ 5 ) ; sum( x

. ˆ 3 )9 ,sum( x . ˆ 4 ) ,sum( x . ˆ 5 ) ,sum( x . ˆ 6 ) ]

10 mx =11 4.000000 1.100000 0.450000 0.19700012 1.100000 0.450000 0.197000 0.08970013 0.450000 0.197000 0.089700 0.04181014 0.197000 0.089700 0.041810 0.01978515 octave :8> my=[sum( y ) ; sum( x .∗ y ) ; sum( x . ˆ 2 . ∗ y ) ; sum( x . ˆ 3 . ∗ y ) ]16 my =17 11.1970018 4.5897019 2.0118120 0.9167921 octave :11> c=mx∗∗−1∗my22 c =23 2.1285 e−1524 1 .0000 e+0125 −4.5830e−1326 1 .0000 e+00

Ficando o modelo y = 2.1285e− 15 + 1.0000e+ 01x− 4.5830e− 13x21.0000e+ 00x3 eaplicando o valor de x = 0.3 temos:

1 octave :10> xd = [ 1 ; 0 . 3 ; 0 . 3 ˆ 2 ; 0 . 3 ˆ 3 ]2 xd =34 1.0000005 0.3000006 0.0900007 0.0270008 octave :12> sum( c .∗ xd )9 ans = 3.0270

Teremos o mesmo resultado, P3(0.3) = 3.027, que o interpolador de lagrange.



Exercıcios Basicos

4.2 Interpolacao por lagrange

Teoria

Desenvolver o conteudo selecionado da seguinte bibliografia: [Barroso, 1987], Capıtulo 4.5, Pag.164-165, 167-170 e [Freitas, 2000], Capıtulo 7.3.1, Pag. 139 a 140].

Observacoes basicas na aplicacao do metodo:

1. Usa um Polinomio Interpolador de grau n− 1 sendo a quantidade de pontos n.

2. Resulta apenas em valor numerico




Pede-se:

(a) Utilizando o polinomio interpolador de Lagrange calcule o valor de P3(0.3)Resposta: Resposta: 3.027

Solucao:

O modelo generico do interpolador de Lagrange e:

Pn(x) =n∑i=0

yi ∗n∏j=0j 6=i

x− xjxi − xj

(4.1)

Para nosso problema temos n = 3 e x = 0.3, ou seja estimar um polinomio interpolador



de grau 3 para x = 0.3. Montando o somatorio e produtorio temos:

Pn(x) = y0∗x− x1x0 − x1

∗ x− x2x0 − x2

∗ x− x3x0 − x3

+ y1∗x− x0x1 − x0

∗ x− x2x1 − x2

∗ x− x3x1 − x3

+ y2∗x− x0x2 − x0

∗ x− x1x2 − x1

∗ x− x3x2 − x3

+ y3∗x− x0x3 − x0

∗ x− x1x3 − x1

∗ x− x2x3 − x2

(4.2)

Aplicando os numeros nas variaveis teremos as seguintes fracoes:

P3(0.3) = 0.000∗ 0.3− 0.2

0.0− 0.2∗ 0.3− 0.4

0.0− 0.4∗ 0.3− 0.5

0.0− 0.5

+ 2.008∗ 0.3− 0.0

0.2− 0.0∗ 0.3− 0.4

0.2− 0.4∗ 0.3− 0.5

0.2− 0.5

+ 4.064∗ 0.3− 0.0

0.4− 0.0∗ 0.3− 0.2

0.4− 0.2∗ 0.3− 0.5

0.4− 0.5

+ 5.125∗ 0.3− 0.0

0.5− 0.0∗ 0.3− 0.2

0.5− 0.2∗ 0.3− 0.4

0.5− 0.4

(4.3)

Resolvendo as fracoes teremos as seguintes numeros:

P3(0.3) = 0.000∗ −0.5∗ 0.25∗ 0.4

+ 2.008∗ 1.5∗ 0.5∗ 0.6666

+ 4.064∗ 0.75∗ 0.5∗ 2

+ 5.125∗ 0.6∗ 0.3333∗ −1

(4.4)

Efetuando os calculos temos:P3(0.3) = 3.027 (4.5)

Exercıcios Basicos

4.3 Interpolacao por diferenca dividida

Teoria

Desenvolver o conteudo selecionado da seguinte bibliografia: [Barroso, 1987], Capıtulo 4.6, Pag.175-189. Observacoes basicas na aplicacao do metodo:



1. Usa um Polinomio Interpolador de grau n− 1 sendo a quantidade de pontos n.

2. Resulta apenas em valor numerico




Pede-se:

(a) Utilizando o polinomio interpolador Newton calcule o valor de P3(0.3)Resposta: Resposta: 3.027

Solucao:O modelo generico do interpolador de Newton e:

Pn(x) = y0 +n∑i=1

Niy0 ∗i−1∏j=0

(x− xj) (4.6)

Para nosso problema temos n = 3 e x = 0.3, ou seja estimar um polinomio interpoladorde newtom de grau 3 para x = 0.3. Veja que temos que fazer a diferenca dividida aten = 3. Calculando primeiro a diferenca dividida, temos:

i xi yi Nyi N2yi N3yi

0 0.0 0.000y1 − y0x1 − x0

Ny1 − Ny0x2 − x0

N2y1 − N2y0x3 − x0

1 0.2 2.008y2 − y1x2 − x1

Ny2 − Ny1x3 − x1

2 0.4 4.064y3 − y2x3 − x2

3 0.5 5.125

(4.7)



Substituindo as equacoes pelos numeros teremos:

i xi yi Nyi N2yi N3yi

0 0.0 0.0002.008− 0.000

0.2− 0.0

Ny1 − Ny00.4− 0.0

N2y1 − N2y00.5− 0.0

1 0.2 2.0084.064− 2.008

0.4− 0.2

Ny2 − Ny10.5− 0.2

2 0.4 4.0645.125− 4.064

0.5− 0.4

3 0.5 5.125

(4.8)

Calculando os valores das fracoes temos:

i xi yi Nyi N2yi N3yi0 0.0 0.000 10.04 0.6 11 0.2 2.008 10.28 1.12 0.4 4.064 10.613 0.5 5.125

(4.9)

De posse das diferencas divididas podemos aplicar no modelo original de interpolacaode Newton.

Pn(x) = y0+ Ny0∗ (x− x0)

+ N2y0∗ (x− x0)∗ (x− x1)

+ N3y0∗ (x− x0)∗ (x− x1) (x− x2)

(4.10)

Substituindo as equacoes pelos numeros teremos:

P3(0.3) = 0.000+ 10.040∗ (0.3− 0.0)

+ 0.6∗ (0.3− 0.0)∗ (0.3− 0.2)

+ 1∗ (0.3− 0.0)∗ (0.3− 0.2)∗ (0.3− 0.4)

(4.11)

Calculando as somas e produtos encontramos P3(0.3) = 3.027.

Para este procedimento podemos usar tambem uma tabela para facilitar diferencas,somantorios e produtorios da ultima parte da equacao.

x 0.3 0.3 0.3xi 0 0.2 0.4

x− xi 0.3 0.1 −0.1∏i−1j=0(x− xj) 0.3 (0.3 ∗ 0.1) = 0.03 (0.3 ∗ 0.1 ∗ −0.1) = −0.003

(4.12)



Desta forma fica simples utilizar os resultados das Tabelas 4.9 e 4.12 para aplicar nomodelo principal Eq. 4.6.

P3(0.3) = 0.000+ 10.040∗ 0.3

+ 0.6∗ 0.03

+ 1∗ −0.003

(4.13)

Calculando as somas e produtos encontramos P3(0.3) = 3.027 como anteriormente.

Exercıcios Basicos



1. Fonte: [Barroso, 1987, Pag. 152, Exemplo 4.1] (exercicios/Barroso152exp4.1) Em Belo Ho-rizonte em 1950 haviam 352724 hab, em 1960 haviam 683908 hab, em 1970 haviam 1235030hab e em 1980 haviam 1814990 hab.Pede-se:

(a) Calcular o numero aproximado de habitantes de Belo Horizonte em 1975, com basenum polinomial de 2o grau.

2. Fonte: [Barroso, 1987, Pag. 159, Exercıcio 4.3.3.1] (exercicios/Barroso159exc1) Data afuncao f(x) = 10x4 + 2x+ 1 com os valores de f(0.1) e f(0.2) determinar P1(0.15).Pede-se:

(a) Compare os valores funcao original com a funcao e valores interpolados.

3. [Barroso, 1987, Pag. 159, Exercıcio 4.3.3.3]

4. Fonte: [Barroso, 1987, Pag. 164, Exercıcio 4.4.3.3] (exercicios/Barroso164exc3) Data afuncao f(x) = 10x4 +2x+1 com os valores de f(0.1) e f(0.2) e f(0.3) determinar P2(0.15).Pede-se:

(a) Compare os valores funcao original com a funcao e valores interpolados.













15. [Barroso, 1987, Pag. 201, Exercıcio 4.9.3]













Capıtulo 5

Zero de Funcao

Introducao

Desenvolver o conteudo selecionado da seguinte bibliografia: [Barroso, 1987], Capıtulo 3.1, Pag.83-86.

5.1 Metodo da Bissecao

Teoria

Desenvolver o conteudo selecionado da seguinte bibliografia: [Barroso, 1987], Capıtulo 3.4, Pag.106-107, 109-110.


Exercıcios Basicos

5.2 Metodo de Pegaso

Teoria



Exercıcios Basicos

5.3 Metodo de Newton

Teoria


45



Exercıcios Basicos

5.4 Comparacao dos metodos

Teoria



Exercıcios Basicos



Plote as funcoes caso necessario para determinar um intervalo aproximado para as iteracoes.

1. Fonte: [Barroso, 1987, Pag. 129, Exemplo 3.29]. (exercicios/Barroso129Exp3.29)



4. Fonte: [Barroso, 1987, Pag. 120, Exemplo 3.25] (exercicios/Barroso120Exp3.25)






Capıtulo 6

Integracao

Introducao

Desenvolver o conteudo selecionado da seguinte bibliografia: [Barroso, 1987], Capıtulo 5.1, Pag.205-206. [Freitas, 2000], Capıtulo 8, Pag. 163

6.1 Metodo de Newton-Cottes

6.1.1 Trapezio

Teoria

Desenvolver o conteudo selecionado da seguinte bibliografia: [Barroso, 1987], Capıtulo 5.2, Pag.206-209. [Freitas, 2000], Capıtulo 8.1, Pag. 164-168.


Exercıcios Basicos

6.1.2 Simpson I

Teoria

Desenvolver o conteudo selecionado da seguinte bibliografia: [Freitas, 2000], Capıtulo 8.2, Pag.169-173. [Barroso, 1987], Capıtulo 5.3,5.4; Pag. 214-232.


Exercıcios Basicos

6.1.3 Simpson II

Teoria

Desenvolver o conteudo selecionado da seguinte bibliografia: [Freitas, 2000], Capıtulo 8.2, Pag.169-173. [Barroso, 1987], Capıtulo 5.3,5.4; Pag. 214-232.

47



Exercıcios Basicos

6.1.4 Extrapolacao Richardson

Teoria



Exercıcios Basicos

6.1.5 Integral dupla

Teoria

Desenvolver o conteudo selecionado da seguinte bibliografia: [Barroso, 1987], Capıtulo 5.6, Pag.243.


Exercıcios Basicos

6.2 Quadratura Gaussiana

Teoria



Exercıcios Basicos



Nao precisa calcular o erro cometido, nos exercıcios a seguir.



























Capıtulo 7

Equacoes Diferenciais Ordinarias

Introducao

Desenvolver o conteudo selecionado da seguinte bibliografia: [Barroso, 1987], Capıtulo 6.1 a 6.1.2,Pag. 275-279.

Desenvolver o conteudo selecionado da seguinte bibliografia: [Barroso, 1987], Capıtulo 6.1.3,Pag. 279-283.

7.1 Metodos de Runge-Kutta

Teoria



Exercıcios Basicos



Aplique apenas o metodo de Runge-Kutta de 4a ordem.




50

Parte III

Administracao da disciplina

51

Capıtulo 8

Dicas de como estudar

Parece ate um tanto ousado um topico como este, mas por incrıvel que pareca existem algumasorientacoes que sao importantes de como se deve proceder para ter sucesso no estudo de umconhecimento qualquer, como no do presente curso.

Voce ao ler estas dicas pode dizer, mas isto nao serve pra mim. No entanto de tudo que voceler neste topico alguma coisa lhe servira de auxılio no seu estudo, e ja valeu a pena todo o esforcodeste texto.

Produtividade, e quando precisamos fazer muito em pouco tempo/espaco. Para isto precisamosaos poucos construir uma base solida para que a cada etapa possamos ir aumentando nossaprodutividade nos estudos. Nem sempre ir direto ao resumo resolve o problema, ele e um guiasobre o que sabemos de um assunto, se nao sabemos nada o resumo tem pouca validade.

8.1 Rotina

E muito importante voce estabelecer uma rotina de estudos para cada materia, lendo sempre oconteudo basico e nao apenas os resumos ou apostilas. E sem medo

Separe tempo para cada coisa: famılia, igreja, saude e estudos. Viu como os estudos vem porultimo, pois se vc parar para estudar e algum dos demais estiverem em falta, seu organismo/-conciencia comecara a cobrar a conta e vc nao tera paz para produzir com qualidade.

Tirar alguns momentos para programar um descanso/lazer e muito importante. Ficar dias edias trabalhando nao tras nenhum retorno a medio ou longo prazo.

8.2 Morto

Estudar cansado nao rende nada, com 30 minutos vc ja esta dormindo. Entao esteja sempredescansado antes das atividades de estudo. Nao adiante dormir a tarde toda e estudar ate as 4 damadrugada, o dia seguinte cobra a conta. Apos se acatalunhaentar ou fazer uma atividade fısicao organismo cobra um descanso, permita se relaxar apos o almoco por alguns minutos, uns 20sao um bom referencial. Aprenda a dormir numa rede, ela e confortavel, leve, higienica e facil detransportar.

52


8.3 Faminto

Com fome ficamos muito incomodados ao parar para estudar, entao faca um lanche antes deestudar, coisa leve. E ruim tambem comer enquanto estuda. O organismo nao gosta desta ati-vidades paralelas com a acatalunhaentacao. Acabar de almocar e estudar tambem e ruim, poisapos as refeicoes precisamos dar ao organismo oportunidade de processar a acatalunhaentacaoe isto consome energia e por isto ficamos mais cansados apos uma boa refeicao, importante daroportunidade, com um cochilo, para o organismo trabalhar em paz. Evite ouvir musica apos oalmoco, feche os olhos e deixe o organismo aproveitar a energia para processar o almoco.

8.4 Anti-social

Nao adianta vc querer estudar ligado no facebook, ou outra rede social. A mente precisa seconcentrar, de dedicacao para entender o conteudo a ser assimilado. Apos estudar algumas horas,vc para um descanso e ai sim se atualize, mas aquele conceito de conectado a todo momento,cobra seu preco na falta de rendimento nos estudos, vc parace que estudou horas, mas na verdadea mente aprendeu pouco coisa. Ai vc se irrita com o conteudo pois fica horas ’estudando’ e naoaprende, mas veja que vc nao estudou, ficou apenas olhando para o conteudo e concentrado naconversa da amiga sobre as mais recentes fofocas do dia. Ser anti-social no momento de estudare prova de equilıbrio e sucesso em algo mais importante, assimilar o conteudo com qualidade.

Ouvir musica e outro fator questionado. Muitos dizem, estudo apenas com musica, nao sei aocerto ate que ponto nossa mente consegue assimilar ambas as atividades com qualidade. Sempreexiste perda de alguma. E a musica e um lazer, entao reserve tempo de qualidade para ela e naodurante seu estudo. Igualmente para a televisao.

8.5 Arsenal

Antes de estudar, veja se vc possui todo material necessario para entendimento daquele conteudo.Nao perca tempo baixando, via internet, toneladas de tutoriais ou pegando dezenas de livros nabiblioteca, escolha alguns e tente entender a teoria partindo deles.

Estudar de apenas uma fonte tambem e ruim, pois nem todos os livros tem o mesmo foco oulinguagem, para isto tenha duas a quatro fontes boas como base de seus estudos.

Construa seus resumos bem organizadamente. Os resumos sao a prova de que sua menteorganizou o entendimento. Apos isto resolva os exercıcios somente baseando nos resumos. Antesda prova estude pela teoria e reforce o conhecimento nos resumos.


Capıtulo 9

Procedimentos de avaliacao

Com o intuito de esclarecer e evitar alguns problemas sobre o procedimento de avaliacao, seguemalgumas regras que DEVEM ser seguidas no presente EDITAL.

9.1 Conteudo e avaliacoes

1. O conteudo conforme listado neste edital sera ministrado continuamente ao longo do semes-tre.

2. A avaliacao do desempenho do aluno dar-se-a por trabalhos de casa, trabalhos em sala, 1projeto e 1 prova.

3. Nas tarefas ou Prova cada letra (a,b,...) da questao (1,2,...) vale 1 (um) ponto ou mais.Somando todas a pontuacao de todas as letras acertadas e dividindo pelo total de letrastemos sua pontuacao total. Este valor sera aplicado ao peso da atividade. Veja Tabela 9.1e Equacoes.

Tabela 9.1: Valor das questoes numa provaProva/Tarefa Letra Valor

1

a 1b 1c 1d 1

2a 1b 1

3 a 1Total 7/7=1

4. Sera utilizado planilha eletronica, para calculo de notas e medias, sem interferencia subjetivanas notas ou casas decimais. No diario a nota 1 corresponde a nota da etapa 1, idem paranota 2. O lancamento da nota sera conforme padrao daquele sistema, uma casa decimal,arredondados conforme planilha.

54


5. As tarefas de Casa e Sala tem juntas peso 6, o projeto tem peso 2 e a prova peso 2. Nacorrecao as notas valem zero ou um para cada letra. O projeto tem avaliacao subjetiva dezero a um.

NotaF inal =(TCasa+ TSala) ∗ Peso+NotaProj ∗ Peso+NotaProva ∗ Peso∑

(PesoTarefas+ PesoProj + PesoProva)

NotaDiario = NotaF inal ∗ 10

9.2 Estrutura do trabalho

6. Cada tarefa tem um valor especıfico e a soma de todas as tarefas corretas divida pela somatotal de tarefas valera de zero a um. Multiplando este valor pelo peso total das tarefas.Somente serao pontuadas tarefas entregues completas.

7. O sistema de gestao de tarefas sera informado via moodle. Nele serao postados e corrigidostodos trabalhos.

8. Sera considerada uma situacao, exercicio, onde se pede varios calculos, um por letra.

9. O tarefas serao reunidas por capıtulo/conteudo. Durante todo o capıtulo quaisquer per-guntas podem ser feitas sobre o mesmo. O aluno deve ficar atento para a publicacao dastarefas.

10. Caso na tarefa em sala o aluno force a saıda igual a solicitada sem os calculos corretos omesmo perdera todas as notas de todas as tarefas ja conquistadas.

9.2.1 Para casa

11. Cada tarefa tem prazo de entrega definido com data e hora inicial e data e hora final. Estesprazos sao improrrogaveis.

12. Cada tarefa tem um tempo limite de resolucao em horas e minutos. Estes prazos sao im-prorrogaveis.

13. Cada tarefa tem um limite de tentativas. Estes limites sao improrrogaveis.

9.2.2 Em sala

14. As tarefas que forem aplicados em sala de aula, terao nota de 100% se entregues duranteo prazo estabelecido em sala. E nota de 50% se aplicados fora de sala, caso haja falta doaluno naquela aula. Todos as tarefas aplicados em sala serao avisados na aula corrente.

15. As tarefas aplicados em sala serao um dos aplicados para casa. O sorteio sera feito aleatoriopara cada aluno.



9.3 Estrutura do Projeto

16. O projeto substitui a prova 01. O aluno pode abandonar o projeto e realizar a prova.

17. O projeto sera a pesquisa e desenvolvimento completo de um exercıcio de calculo numeropor capıtulo. Devendo ser entregue no formado a seguir:

18. Nao sera aceito projetos simulares.

19. A correcao do trabalho e totalmente subjetivo a criterio do professor. Valendo de zero a um.Multiplicando assim pelo peso do projeto.

20. O arquivo com o relatorio, somente texto puro, deve ser entregue no sistema de tarefas coma seguinte estrutura. Substitua os ”[]”pelo que se pede:

1 ############################################################2 Fonte :3 −−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−45 [ informe onde vc obteve o exerc ı c i o . Pode s e r um l i v r o ( informe nome

completo do l i v r o , autor , pagina e exerc ı c i o ) ou s i t e ( c i t e caminhocompleto ) ]

67 ############################################################8 Introdu c ao te o r i c a :9 −−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−

1011 [ Apresente uma d e s c r i c ao do contexto do problema abordado na s i t u a c ao . ]1213 ############################################################14 D e f i n i c ao da s i t u a c ao espec ı f i c a :15 −−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−1617 [ I s o l e uma s i t u a c ao a s e r r e s o l v i d a ]1819 ############################################################20 O que se pede da s i t u a c ao :21 −−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−2223 [ Determine o que se de s e j a obter da s i t u a c ao em d i f e r e n t e s l e t r a s ]2425 ############################################################26 Resposta :27 −−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−2829 [ Apenas a r e spo s ta s imp le s e d i r e t a ]3031 ############################################################32 Desenvolvimento :33 −−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−−



3435 [O desenvolvimento deta lhado e s c r i t o da s egu in t e forma : ]3637 Exemplo :38 y=a+bx39 a=140 b=241 x=342 y=1+2∗343 y=744 ]

9.4 Estrutura da prova:

9.4.1 Preparacao da sala para a avaliacao

21. A avaliacao sera realizada em sala de aula comum. Podendo haver modificacao de sala pararealizacao da prova, devendo o aluno ficar atendo ao moodle onde sera publicada a nova salapara avaliacao.

22. A carteira onde o aluno devera se colocar dentro da sala de aula sera decidida pelo professor.Durante a prova o aluno podera ser modificado de carteira sem previo aviso.

23. Os alunos devem organizar a sala de forma que as carteiras formem circulos, podendo ser maisde um circulo na sala. Distribuıda de forma equidistante entre as carteiras dentro da sala.Cada circulo deve ter o mesmo numero de carteiras. As carteiras posicionadas proximasas paredes poderao permanecer com o respectivo lado encostado na mesma. Sempre quepossıvel deve haver uma carteira vazia entre dois alunos.

24. Para inıcio da prova os alunos devem portar apenas roupa pessoal, oculos (caso necessariopara auxilio na visao), borracha, grafite, caneta, regua, calculadora (autorizada na disci-plina). Os demais objetos devem ser deixados dentro da bolsa na frente da sala.

25. Caso o aluno esteja portando blusa, a mesma deve estar sendo usada, senao deve deixar nafrente da sala. A calculadora nao deve permitir leitura de textos longos ou imagens, sendoque o aluno devera consultar o professor sobre se sua calculadora atende a este requisito parao mesmo realizar a prova. A realizacao da prova preve necessariamente o uso de calculadorae portanto o aluno nao deve presumir que tera condicoes de realizar a prova sem portar uma.

26. Qualquer outro objeto eletronico, celular, relogio, etc, devera ser deixado dentro da bolsa edesligado. O aluno nao tera acesso a sua bolsa durante a realizacao da prova, caso o alarmeou outros sinais sonoros acontecam durante a prova o aluno nao podera ter acesso a seuobjeto para manuseio. O aluno deve tambem certificar de que nao esta de posse de nenhumoutro papel, material ou objeto senao os autorizados. Acarretara o cancelamento da provado aluno, mesmo que o aluno esteja portanto os objetos especificados neste item, quer sejapor esquecimento ou outro motivo.



27. Como alimento durante a prova somente sera autorizado o uso de agua em recipiente trans-parente. Portanto o aluno devera efetuar um adequado desjejum antes do inıcio da prova.Que tera duracao de 200 minutos.

28. O acesso ao banheiro esta restrito a casos necessarios. Sendo a liberacao feita pelo professore nao pela vontade propria do aluno. Devendo o aluno sempre cientificar o professor de suasaıda e chegada.

29. Todo material necessario para a realizacao da prova, fora os citados acima serao entreguespelo professor e deverao ser devolvidos no final da prova. Nao ficando com o aluno nenhummaterial usada durante estas atividades. As folhas de rascunho serao ecatalunhainados nomomento do recolhimento da prova e as folhas de desenvolvimento devidamente identificadasserao anexadas a prova.

30. O aluno que estiver portanto qualquer material bibliografico referente a disciplina, fora oentregue pelo professor, provocara o cancelamento imediato de sua prova.

9.4.2 A prova

31. O aluno dispoe de 1 encontro, 200 minutos, para resolucao da prova. Para isto o aluno devecontrolar individualmente o tempo de resolucao de cada questao, bem como o tempo total.

32. O aluno recebera uma folha tipo A4, ou mais, onde constam alguns questoes e o que se pedede cada um. Nesta folha nao deve ser feita nenhuma resolucao.

33. Cada questao devera ser resolvida individualmente em folha tipo A4, fornecida para de-senvolvimento. O aluno deve usar uma linha dupla para separacao entre cada questao. Oaluno deve evitar uso de duas colunas na folha ou quando fizer, usar uma linha vertical paraseparar o desenvolvimento. Devendo evitar o uso de linhas diagonais. Veja Capıtulo 10 paraentendimento do procedimento de elaboracao de tarefas.

34. A presenca do aluno na prova se da pela assinatura da folha de presenca.

35. A resolucao de uma questao tem duas partes: (1) O aluno devera efetuar o desenvolvimentoa lapis ou caneta. (2) A resposta somente a caneta, deve ser um texto claro onde consta ovalor numerico encontrado e a respectiva unidade quando for o caso.

36. A prova consta basicamente de desenvolvimento de qualquer problema matematico praticoutilizando conhecimento teorico constante na bibliografia e ministrado em sala. E nao ape-nas dos exercıcios entregues pelo aluno ou constantes na lista de exercıcios resolvidos oupropostos ou presentes na bibliografia.

37. As operacoes basicas da calculadora compreendem soma, subtracao, multiplicacao, divisao,operacoes trigonometricas, operacoes de potencia, operacoes logarıtmicas e somatorio delista de numeros. As demais operacoes e desenvolvimentos devem ser apresentados, como detalhamento solicitado, como parte da resolucao da questao. A equacoes ou modelosmatematicos devem apresentar todos valores envolvidos nas suas variaveis e parametros, e



podem ser resolvidas diretamente na calculadora. O uso de funcoes avancadas, nao constan-tes da lista anterior, podera ser usada apenas para conferencia das resposta e nao constaracomo resposta do exercıcio.

38. A resolucao da questao, e a obtencao da nota da mesma, se da pela elaboracao completae conjunta de todos os calculos em cada passo, apresentando-os de forma clara, coerente,didatica, explicativa, organizada, legıvel e detalhada. E da resposta coerente com o desenvol-vimento. Nao serao consideradas fracoes de calculo ou desenvolvimento como atendimento anota da questao. A nota compoe-se da coerencia de todos os calculos envolvidos no exercıcioe nao somente da resposta final. E o aluno nao deve presumir que a escrita de algarismossoltos na folha, ou pedacos desconectos de desenvolvimento, irao justificar o desenvolvimentoe a resposta ao exercıcio. A divergencia entre desenvolvimento e resposta anula a questao.

39. Deverao ser utilizada 4 casas decimais. Com numeros diferentes de zero, em notacao ci-entıfica, para todos os calculos. O erro relativo maximo permitido entre o valor do ga-barito e o valor fornecido pelo aluno e de 1%. com base na equacao ErroRelativo =V alorGabarito− V alorAluno

V alorAluno∗ 100.

40. Em situacoes de iteracoes quando o erro nao for informado, sera adotado o erro padrao de0.01 entre o valor estimado e a iteracao anterior.

41. A plotagem ou construcao de graficos e obrigatoria para qualquer modelo discutido na dis-ciplina, bem como aplicacoes de conhecimentos de calculo e fısica, ja vistos em outras disci-plinas do curso. Sendo que as derivadas de modelos polinomiais e obrigatorio o aluno sabere as demais derivadas sao fornecidas.

42. O professor recebera a prova pessoalmente de cada aluno, somente se todas as informacoessolicitadas na prova estiverem preenchidas. Podendo o aluno sair de sala somente aposliberacao do professor, o abandono da prova em sala caracteriza desistencia da mesma e seraarrolada testemunha para tais fatos.

9.5 Outras orientacoes:

43. O exame final, que envolve toda a materia, tem os mesmos procedimentos de uma etapa,sem a atividade de trabalho.

44. A turma devera eleger um representante que sera o responsavel por organizar as duvidas/necessidadesgerais da turma e me consultar pessoalmente, via email ou celular para que possamos decidirsobre questoes de natureza geral. Isto nao inclui as atividades de conteudo ou duvidas damateria, trata-se apenas de assuntos gerais como trabalhos, provas, reposicoes, bem comooutras.

45. As provas ficarao com o professor, a disposicao do aluno, para revisoes de nota em data aser definida com o professor. Apos este perıodo estas atividades serao arquivadas. Exceto aprova de recuperacao, que ficara na Secretaria Academica do Curso.



46. O professor a todo momento podera revisar as notas lancadas, ate o fechamento do diariono final do semestre.

47. O aluno devera assinar lista de presenca em sala de aula e no atendimento a alunos noperıodo programado, como comprovacao de sua presenca naquela atividade.

48. As regras previstas no capıtulo Elaboracao de Tarefas,Capıtulo 10 serao adotadas comoreferencia para todas as tarefas: sejam elas provas, trabalhos, projetos, questoes, exercıciosetc; desenvolvidas durante o curso. Somente sendo considerados para correcao as tarefasque estiverem de acordo com estas regras.

49. A administracao da disciplina sera feita somente via moodle ou link relacionado neste.

50. A prova de recuperacao, que envolve toda a materia, tem os mesmos procedimentos de umaetapa, sem a atividade de trabalho.

51. A solicitacao para reposicao de provas seguira os tramites da Secretaria Academica do Curso,SAC, e sera analisadas de acordo com cada caso pelo professor. Sendo marcada uma unicadata para realizacao destas atividades.

52. Todas as atividades da disciplina como aulas, trabalhos e provas, estao previamente agen-dados e sao apresentados aos alunos no primeiro dia de aula. E constam do EDITAL dadisciplina.

53. A mudanca de qualquer compromisso requer concordancia unanime dos alunos, disponibi-lidade de sala e disponibilidade do professor. O pedido de alteracao de deve ser feito porescrito e em tempo habil para estas mudancas.

54. As provas podem ser alteradas apenas uma vez da data da inicialmente prevista no inıciodo curso.

55. As notas serao divulgadas via moodle, logo apos a realizacao das correcoes que poderaocorrer mesmo depois da realizacao da prova seguinte. Sendo que as provas que estiveremcom dificuldade de interpretacao serao corrigidas por ultimo. O aluno nao deve presumirque deve ter acesso a nota anterior para realizacao de uma prova seguinte na agenda deavaliacoes.

Quaisquer outros procedimentos nao previstos neste texto serao resolvidos pelo professor eatualizados neste EDITAL.


Capıtulo 10

Elaboracao das Tarefas

Uma disciplina e composta basicamente da ministracao/apresentacao de conteudo teorico e de-senvolvimento de atividades para fixacao e ampliacao destes conteudos.

O Moodle referencia qualquer atividade da disciplina (exercıcios, trabalhos, provas) comotarefa, o que sera seguido nestas orientacoes.

Para desenvolvimento das tarefas sera utilizado formulario de equacoes, material de escrita ecalculadora.

Uma tarefa sao situacoes a resolver oriundas de diversas fontes; numeradas assim 1,2,... . Eo que se deseja deste exercıcio no item ”Pede-se”numeradas assim (a),(b),... . Que devem serdesenvolvida no relatorio manual individualmente em cada letra conforme orientacoes a seguir.Seguindo sempre de uma resposta que coaduna com o desenvolvimento.

10.1 Tarefa exemplo

Considere que foi proposto uma tarefa ao aluno de nome Fulano Ferreira da Silva. Para CalculoNumerico deve ser feito desenvolvimento manual com resultados e formatacao similares a umdesenvolvimento digital apresentado.

A referencia a tarefa no item ”Exercıcios Propostos”de cada capıtulo e apresentada da seguinteforma:

10.1.1 Exercıcios Propostos

1. Considere os valores da Tabela 10.1 de um experimento onde y=f(x).

Tabela 10.1: Experimentox y

1.3 23.4 5.25.1 3.86.8 6.18 5.8

Pede-se:

61


(a) Esboce um grafico dos dados

(b) Ajuste um modelo polinomial de 1 grau, para os pontos dados.

(c) Calcule a raiz deste modelo.

(d) Integre o modelo num intervalo qualquer.

2. Veja na Bibliografia [Barroso, 1987, Pag. 335, Exemplo 7.5].

Veja que esta tarefa tem 2 exercicios. O primeiro constante na lista de exercıcios propostos.O segundo aponta para determinado livro usado na disciplina.

A estrutura de desenvolvimento deve seguir a orientacao conforme abaixo para que, quando eufor consulta-lo terei um entendimento claro do que voce quis fazer, de quais ferramentas utilizoue das manipulacoes matematicas necessarias ao desenvolvimento do exercıcio. Vendo claramenteo desenvolvimento e a resposta.

Evite usar coluna dupla na folha A4. Mas necessitando faze-lo utilize uma linha vertical paraseparar os conteudos.

Considere o retangulo a seguir como a primeira folha A4 utilizada na resolucao da tarefa.

1 Prova 01 de Calc . Num. − Fulano F e r r e i r a da S i l v a − f o l h a 1 de 22 ############################################################3 E x e r c i c i o 01 − LETRA A4 ############################################################5 DESENVOLVIMENTO:6 Neste item voc e desenvo lve toda a r e s o l u c ao da r e f e r i d a l e t r a do exerc ı c i o7 RESPOSTA:8 Ao f i n a l voc e DEVE e s c r e v e r um texto c la ro , a caneta , como re spo s ta da

l e t r a em desenvolvimento , c i tando os v a l o r e s nume r i c o s , v a r i a v e i s eunidades nece s s a r i a s para composi c ao da r e spo s ta .

9 ############################################################10 E x e r c i c i o 01 − LETRA B11 ############################################################12 DESENVOLVIMENTO:13 . . . . .14 RESPOSTA:15 . . .

Considere o retangulo a seguir como a segunda folha A4 utilizada na resolucao da tarefa.

1 Calc . Num. − Prova 02 − Fulano F e r r e i r a da S i l v a − f o l h a 2−22 ############################################################3 cont inua c ao do E x e r c i c i o 1 − Letra B.4 ############################################################5 DESENVOLVIMENTO:6 . . . . .7 RESPOSTA:8 . . .9

10 ############################################################11 E x e r c i c i o 02 − LETRA A12 ############################################################



13 DESENVOLVIMENTO:14 . . . . .15 RESPOSTA:16 . . .

Quaisquer outras orientacoes neste item serao atualizadas e informadas aos alunos.


Capıtulo 11

Formulario

O uso do formulario e de extrema importancia em todas as atividades, inclusive nas provas, entaoestude pelo livro mas resolva os exercıcios usando apenas o formulario.

A estrutura do mesmo e apresentada em anexo.

64

CÁLCULO NUMÉRICO1) Erros

e x=∣x−x∣ , E x=ex

∣x∣

2) Sistemas Lineares:

a11 x1a12 x2...a1n x n=b1

a21 x1a22 x2...a2n xn=b2

...an1 x1an2 x2...ann xn=bn

(1) ou ∑j=1

n

a ij x j=bi ou

Ax=b (2) com i=j=1,2,...,n; i=linha; j=coluna

2.1) Métodos diretos:

Etapa k ; pivô: a i p j p

k ; multiplicador: i≠i p ; mik=−

aij p

k

a i p j p

k;

linha: Lik1=Li

kLi p

k∗mik

Para Gauss: pivô é a i p j p onde i p= j p transformando a Li onde

ii p e Li p fica inalterada.

Para Pivotação Completa: pivô é a i p j p onde máx∣aij≠0∣

transformando a Li onde i≠i p e Li p fica inalterada.

Para Jordam: pivô é a i p j p onde i p= j p transformando a Li

onde i≠i p

2.2) Métodos iterativos:É condição suficiente para que a iteração de ambos os métodos

convirja se ∣a ii∣> ∑j=1, j≠i

n

∣aij∣ para i=1,2,...,n.

Para ambos os métodos, considere o sistema linear (1). Explicitando

x1 na 1ª equação, x2 na 2ª, …., tem-se:Para Jacobi:

x1(k + 1)

=b1−(a12 x2

(k )+ ...+ a1n xn

(k))

a11

x2(k+ 1)

=b2−(a21 x1

(k )+ ...+ a2n x n

(k ))

a22

...

xn(k+ 1)

=bn−(an1 x1

(k)+ an2 x2

(k )+ ...+ an ,n−1 xn−1

( k))

ann

Para Gauss-Seidel:

x1(k+ 1)

=b1−(a12 x2

(k )+ ...+ a1n xn

(k))

a11

x 2(k+ 1)

=b2−(a21 x1

(k+ 1)+ ...+ a2n xn

(k))

a22

...

xn(k+ 1)

=bn−(an1 x1

(k+ 1)+ an2 x2(k+ 1)+ ...+ an , n−1 xn−1

(k+ 1))

ann

Para ambos os métodos, os valores iniciais são xn

(0)

qualquer. E o

critério de parada é até k> Iterações ou

máx1≤i≤n∣x ik + 1

−x ik∣≤Erro

3) Ajuste de curvasy=c0+c1 x1+c2 x2+...+cg x g

[n ∑ x1 i ∑ x2 i ... ∑ x gi

∑ x1 i ∑ x1 i x1i ∑ x 1i x2 i ... ∑ x1 i xgi

∑ x2 i ∑ x 2 i x1i ∑ x 2 i x2 i ... ∑ x2i x gi

... ... ... ... ...

∑ xgi ∑ xgi x1 i ∑ x gi x2 i ... ∑ xgi xgi

]∗[c0

c1

c2

...cg

]=[∑ y i

∑ x 1 i y i

∑ x 2 i yi

...

∑ x gi y i

][1 x11 x21 ... xg 1

1 x12 x22 ... xg 2

1 x13 x23 ... xg 3

... ... ... ... ...1 x1n x2n ... x gn

]∗[c0

c1

c2

...cg

]=[y1

y2

y3

...yn

]c=( X ' X )

−1( X ' Y )

R2=1−

∑ ( y i−~y i)

2

∑ y i2−

(∑ y i)2

n ou

R2=c ' X ' Y −n( ym)

2

Y ' Y −n ( ym)2

4) Interpolação

4.1) Interpolação Lagrangeana

Pn x =∑i=0

n

y i∏j=0j≠i

n x− x j

x i−x j 4.2) Formula de Newton para Diferença Dividida

1 y i=

y i1− y i

x i1−x i

n y i=

n−1 y i1−n−1 y i

x in−x i

Pnx = yo∑i=1

n

i y o∏

j=0

i−1

x−x j

5) Zeros de funções

5.1) Método da bisseção

f x contínua no intervalo [a ,b] e f (a)∗ f (b)< 0 . Faz-

se x=ab

2. Se f a ∗ f x 0 então intervalo será

[a ,b= x ] . Senão, se f (x )∗ f (b)< 0 intervalo será

[a= x , b ] . Parada ∣b−a∣≤erro ou

iterações≥ln( b−a

erro )ln 2

−1

Formulário de Cálculo Numérico. Prof. Dr. Catalunha - Versão atualizada em 26/03/2015 às 11:06 hs Folha 1 de 2

5.2) Método Pégaso ou Secante

f x é contínua em [a , b] e f (a)∗ f (b)< 0 . Faz-se

x=b−f bb−a

f b− f a . Se f x ∗ f b0 então faz a=b e

f a = f b . Senão, se f x ∗ f b0 então faz a=a e

f a =f a f b

f b f x . Em ambos os casos faz b=x e

f b= f x na próxima iteração.

5.3) Método de Newton

f x é contínua no intervalo [a , b] e f (a)∗ f (b)< 0 . É

condição suficiente para a convergência do método de Newton que:

f ´ x e f ´ ´ x sejam não nulas e preservem o sinal em

a , b e x0 seja tal que f x0∗ f ´ ´ x00 . Então faz

xn1= xn−f xn

f ´ xn

6) Integração

Considere ∫a

b

f x dx e h=b−a

n e x i=a+ hi e

i=0,1 , ... , n

6.1) Regra do Trapézio

I =h2 y02y1...2yn−1 yn

6.2) 1ª regra de Simpson (n=múltiplo de 2)

I =h3 y04y12y24y32y4...2yn−24yn−1 yn

6.3) 2ª regra de Simpson(n=múltiplo de 3)

I=3h8

y03y13y22y33y 43y52y6...3yn−23yn−1 yn

6.4) Richardson

I =I 2n1

p I 2−I 1

n2p−n1

p; p=2 se trapézio e p=4 se simpson. n2>n1

6.5) Integral dupla

I =∫a

b

dx∫c

d

f (x , y )dy=kx ky∑i=0

i=nx

∑j=0

j=ny

((cx i cy j)∗ f ( x i , y j))

kx=hx2

se trapézio; kx=hx3

se 1 reg. de simpson; kx=3hx8

se 2 reg. de simpson; idem para ky.

6.6) Gauss

I =∑i=0

n−1

Ai F t i e

F t =12b−a f 1

2b−a t

12ba sendo Ai e t i

coeficientes conforme tabela a seguir:n i t_i A_i

1 0 0 2

2 1;0 ±0,57735027 1

3 0;12

±0,774596670

5/98/9

4 0;12;3

±0,86113631±0,33998104

0,347854840,65214516

7) EDO

7.1) Runge-Kutta de 4ª Ordem

y j1= y jh6K 12K 22K 3K 4

x j1=x jh

K 1= f x j , y j

K 2= f x jh2

, y jh2

K1

K 3= f x jh2

, y jh2

K2

K 4= f x jh , y jhK3

8) Observações e revisõesPropriedades logarítmicas

ln ab=ln a ln b ln ab=ln a – ln b ln ab=b ln a

ln a=log e a ln e=1

Formulário de Cálculo Numérico. Prof. Dr. Catalunha - Versão atualizada em 26/03/2015 às 11:06 hs Folha 2 de 2

Bibliografia

[Barroso, 1987] Barroso, L. C. (1987). Calculo Numerico (com Aplicacoes) 2a Edicao. EditoraHarbra ltda.

[Buffoni, 2000] Buffoni, S. S. d. O. (2000). Metodos Numericos. Apostila de Introducao aosMetodos Numericos, Edicao 2002, UFF.

[Freitas, 2000] Freitas, S. R. (2000). Metodos Numericos. Edicao 12-01-2000, UFMS.

[Rodney, 2006] Rodney, C. B. (2006). Ensino-aprendizagem com modelagem matematica: umanova estrategia. 3a edicao. Editora Contexto. 2006. Sao Paulo. ISBN 85-7244-207-3.

67

Documents

EDITAL DA DISCIPLINA DE1 CALCULO NUM ERICO · EDITAL DA DISCIPLINA DE1 CALCULO NUM ERICO Prof. Dr. Catalunha Atualizado via LATEX em 26 de Outubro de 2015 as 15:03 1Este EDITAL cont