MNEaula07

Information about MNEaula07

Published on December 28, 2007

Author: Arkwright26

Source: authorstream.com

Content

Métodos Numéricos e Estatísticos:  Métodos Numéricos e Estatísticos Prof. Marcone Jamilson Freitas Souza Aula 7: Métodos numéricos para equações diferenciais 1a ordem Passos múltiplos 2a ordem Equações diferenciais de 1a ordem:  Equações diferenciais de 1a ordem Métodos numéricos são usados quando não é possível obter uma solução geral, ou a forma dela é tão complicada que seu uso não é prático. Uma equação diferencial de 1a ordem tem a forma , e em geral podemos escrevê-la como: Problema do valor inicial - uma equação diferencial - uma condição que deve ser satisfeita pela solução Slide3:  Os métodos que estudaremos partem da idéia de que o espaço da variável independente (x) pode ser discretizado, formando uma rede x0 x1= x0+h x2= x1+h....... h é o passo . O valor da função em cada ponto da rede é calculado a partir de expansões em série de Taylor. Slide4:  Método de Euler ou Euler-Cauchy O valor de y para um passo h é dado pela expansão: Como em geral h é pequeno, suprimimos os termos de ordem O(h2): h2, h3, ..... Resultando na aproximação Slide5:  O que resulta no processo iterativo A omissão dos termos de ordem superior a 2 causa erros de truncagem (que podem ocorrer junto a erros de arredondamento). Slide6:  Exemplo: passo h=0,2 O erro não é (em geral) conhecido. Podemos estimá-lo utilizando um passo h´=2h Método de Euler melhorado (2a ordem):  Método de Euler melhorado (2a ordem) Método chamado de preditor-corretor. Slide8:  Exemplo: o mesmo visto anteriormente Método de Runge-Kutta (4a ordem):  Método de Runge-Kutta (4a ordem) Se f(x,y) não depender de y, o método reduz-se à regra de integração de Simpson Slide10:  Comparação entre os métodos Slide11:  Qual o valor mais adequado para o passo h? Se a função f varia muito com y, então h deve ser pequeno, para evitar erros de truncagem. Em geral, adota-se a “proposta” de que h  h/2 se K  0,05 h  2h se 0,01  K h não muda se 0,05  K  0,01 Estimativa de erro: Métodos para eq. dif. de segunda ordem:  Métodos para eq. dif. de segunda ordem P.V.I. Novamente o problema é obter os valores de yn e yn´ para a seqüência x1 = x0 + h; x2 = x0 + 2h; ... Começamos mais uma vez pelas expansões em série de Taylor da função e de sua derivada: Slide13:  O método mais simples consiste em desprezar os termos em derivadas de ordem y´´´ ou superiores 1o passo: 2o passo: Runge-Kutta-Nyström:  Runge-Kutta-Nyström Valores iniciais: x0, y0, y0´ passo h Saída Equações diferenciais parciais:  Equações diferenciais parciais Uma equação é dita quasilinear se for linear nas derivadas mais altas: Equações de diferenças para Eq. de Laplace e Poisson:  Equações de diferenças para Eq. de Laplace e Poisson Laplace Poisson Vamos ver o caso mais simples em duas dimensões (x e y): Slide17:  (x-h,y) (x,y) (x+h,y) h h k k (x,y-k) (x,y+k) Slide18:  Para as derivadas segundas, desprezando os termos O(h4), temos Juntando as aproximações das derivadas primeiras e segundas, fazendo h=k, obtemos a equação de diferenças correspondente à equação de Poisson: Slide19:  Para f(x,y) = 0 temos a equação de Laplace. h é chamado de o comprimento da malha (mesh size). Equações elípticas - em geral - devem levar em conta problemas de contorno (condições previamente definidas numa dada fronteira - espacial, por exemplo). Casos mais comuns: Dirichlet: se u é definido na fronteira C Neumann: se un=u/n (derivada na direção normal) é definida na fronteira. Para resolver o problema, é necessário criar uma malha.: nós da rede ou da malha (Pij) Fronteira C Exemplo:  Exemplo Uma placa de 12 cm de lado tem suas bordas mantidas às temperaturas mostradas na figura. Quais os valores das temperaturas no interior da placa? Será escolhido um comprimento h = 4 cm. 12 x y 12 u=0 u=100 u=100 u=100 R u=0 u=100 u=100 P02 P10 P20 P01 P11 P21 P12 Slide21:  A equação de transferência de calor é ut = c2(uxx+uyy) Para o regime estacionário ut = 0, a equação se reduz à de Laplace uxx+uyy = 0 Para cada ponto da malha, temos a seguinte equação: ui+1,j + ui-1,j + ui,j+1 + ui,j-1 -4 ui,j = 0 P11: - 4u11 + u21 + u01 + u12 + u10 = 0 - 4u11 + u21 + 100 + u12 + 100 = 0 - 4u11 + u21 + u12 = - 200 ui+1,j ui-1,j ui,j+1 ui,j-1 ui,j Slide22:  - 4u11 + u21 + u12 = -200 u11 - 4u21 + u22 = -200 u11 - 4u12 + u22 = -100 u21 +u12 - 4u22 = -100 Dando como resultados u11 = u21 = 87,5 (88,1) u12 = u22 = 62,5 (61,9)

Related presentations


Other presentations created by Arkwright26

transportation
07. 11. 2007
0 views

transportation

wonderful world
19. 06. 2007
0 views

wonderful world

2006911155950435
28. 04. 2008
0 views

2006911155950435

dietrich
17. 04. 2008
0 views

dietrich

ME Individual DM JG 2006
16. 04. 2008
0 views

ME Individual DM JG 2006

H106n
14. 04. 2008
0 views

H106n

DM GlobalFDI Movements240306
13. 04. 2008
0 views

DM GlobalFDI Movements240306

may30
10. 04. 2008
0 views

may30

Ulad using crop residues
09. 04. 2008
0 views

Ulad using crop residues

coral reef and climate change
07. 04. 2008
0 views

coral reef and climate change

Ian Brinkley DtF 07 06
30. 03. 2008
0 views

Ian Brinkley DtF 07 06

Temperature
14. 02. 2008
0 views

Temperature

AP Review 1400 1800
20. 02. 2008
0 views

AP Review 1400 1800

New Sony
03. 10. 2007
0 views

New Sony

Literary Vocabulary Rhyme
10. 10. 2007
0 views

Literary Vocabulary Rhyme

Chapter1McMurry
13. 10. 2007
0 views

Chapter1McMurry

FinanceTransition
16. 10. 2007
0 views

FinanceTransition

rexcor baker
15. 10. 2007
0 views

rexcor baker

kakande
28. 11. 2007
0 views

kakande

bernsteintwo
16. 10. 2007
0 views

bernsteintwo

What Is Internal Control
29. 10. 2007
0 views

What Is Internal Control

11 40 063
07. 11. 2007
0 views

11 40 063

Il Nazismo
14. 11. 2007
0 views

Il Nazismo

kryukov 20041004
12. 10. 2007
0 views

kryukov 20041004

1015 1
19. 11. 2007
0 views

1015 1

AI 120 Examples
17. 10. 2007
0 views

AI 120 Examples

galaxy physics
01. 12. 2007
0 views

galaxy physics

Qualitative tools
29. 11. 2007
0 views

Qualitative tools

pannebecker
03. 01. 2008
0 views

pannebecker

infectious
05. 01. 2008
0 views

infectious

b e flows
07. 01. 2008
0 views

b e flows

CSI pres
07. 10. 2007
0 views

CSI pres

1A Quality of our Water
02. 01. 2008
0 views

1A Quality of our Water

OPVII AldusEquity
01. 10. 2007
0 views

OPVII AldusEquity

wheat 1
04. 10. 2007
0 views

wheat 1

Lexical Semantics II
21. 11. 2007
0 views

Lexical Semantics II

Forklift Standard 12 14 99
27. 02. 2008
0 views

Forklift Standard 12 14 99

manuel scott powerpoint
25. 03. 2008
0 views

manuel scott powerpoint

subspace
19. 06. 2007
0 views

subspace

skos ecoterm 2006
19. 06. 2007
0 views

skos ecoterm 2006

services
19. 06. 2007
0 views

services

Working with Automatic PGA
19. 06. 2007
0 views

Working with Automatic PGA

wider context
19. 06. 2007
0 views

wider context

weinberg wfi
19. 06. 2007
0 views

weinberg wfi

VS Mod Presentation
19. 06. 2007
0 views

VS Mod Presentation

Unicode from a distance
19. 06. 2007
0 views

Unicode from a distance

Unicode AndIndia
19. 06. 2007
0 views

Unicode AndIndia

tunable abw
19. 06. 2007
0 views

tunable abw

Tsunefum Mizuno sep14 05
19. 06. 2007
0 views

Tsunefum Mizuno sep14 05

tlstut
19. 06. 2007
0 views

tlstut

synergy redesign demo
19. 06. 2007
0 views

synergy redesign demo

acadien
19. 06. 2007
0 views

acadien

y report
19. 06. 2007
0 views

y report

Tom Worthington
19. 06. 2007
0 views

Tom Worthington

Millennials
14. 07. 2007
0 views

Millennials

vienna a6
19. 06. 2007
0 views

vienna a6

unit armorer sustainment
28. 02. 2008
0 views

unit armorer sustainment

SCI1010 C2
13. 11. 2007
0 views

SCI1010 C2

Slides 2006 fin year web3
19. 06. 2007
0 views

Slides 2006 fin year web3

OUR SCAVENGER HUNT edited
16. 11. 2007
0 views

OUR SCAVENGER HUNT edited

ImplicationsResearch
03. 01. 2008
0 views

ImplicationsResearch

Jonh Roberts
31. 07. 2007
0 views

Jonh Roberts

wstechnology
19. 06. 2007
0 views

wstechnology

DiapoAnglaisdÃf
23. 10. 2007
0 views

DiapoAnglaisdÃf

QM chip
15. 10. 2007
0 views

QM chip

zend talk
19. 06. 2007
0 views

zend talk

seminarpresent
24. 02. 2008
0 views

seminarpresent

High School Counsellor Session
23. 11. 2007
0 views

High School Counsellor Session

xml cop feb05
19. 06. 2007
0 views

xml cop feb05

Value of Org RWG
19. 06. 2007
0 views

Value of Org RWG

Promo wkshp Downes
13. 03. 2008
0 views

Promo wkshp Downes

yw
17. 10. 2007
0 views

yw

WP1b
15. 10. 2007
0 views

WP1b

Regency Traffic111508 3 1
11. 03. 2008
0 views

Regency Traffic111508 3 1

Superstar
19. 06. 2007
0 views

Superstar

BUS 400
05. 10. 2007
0 views

BUS 400

950321
11. 10. 2007
0 views

950321

perry presentation
04. 03. 2008
0 views

perry presentation

vergados 1
20. 11. 2007
0 views

vergados 1

vlad
19. 06. 2007
0 views

vlad

Darstellung des HH AZM
15. 11. 2007
0 views

Darstellung des HH AZM