E N D
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas BC1419 - C´ alculo num´ erico - 3º (2016) UFABC 11 de Outubro de 2016
Exemplo: A ´ unica solu¸ c˜ ao do problema de contorno 1 x2y(x) +1 x3,1 < x < 2, y(1) = 1, y(2) =1 y00(x) = 2 ´ e y(x) =1 x. Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β.
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. Exemplo: A ´ unica solu¸ c˜ ao do problema de contorno 1 x2y(x) +1 x3,1 < x < 2, y(1) = 1, y(2) =1 y00(x) = 2 ´ e y(x) =1 x.
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. M´ etodo de diferen¸ cas finitas Fixado um n´ umero inteiro positivo n, estimar os valores da solu¸ c˜ ao y(t) nos pontos igualmente espa¸ cados a = x0< x1< x2< ... < xn= b, xi= x0+ ih,i = 1,2,...,n, h :=b−a n.
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. O vetor de aproxima¸ c˜ ao yhpor diferen¸ cas finitas ´ e a solu¸ c˜ ao do seguinte sitema linear de equa¸ c˜ oes ? ? −r(x1) +α −r(x2) 2 + h2q(x1) h2 −1 0 0 . . . 2 + h2q(x2) −1 −1 0 . . . ... ... ... ... 0 1 h2 × yh= ... ... ... ... ... 0 2 + h2q(xn−2) −1 −1 0 −1 0 −r(xn−2) −r(xn−1) + ? ? 2 + h2q(xn−1) 0 . . . β h2
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. Convergˆ encia do m´ etodo de diferen¸ cas finitas Sob algumas hip´ oteses, vale que y1− y(x1) y2− y(x2) ... yn−1− y(xn−1) , K n2 ||Eh||∞ ≤ Eh :=
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x
Projeto aplicado: resolu¸ c˜ ao de problemas de contorno por diferen¸ cas finitas Problema alvo: dadas duas fun¸ c˜ oes cont´ ınuas q e r, encontrar uma fun¸ c˜ ao suave y : [a,b] → R tal que y00(x) = q(x)y(x) + r(x), a < x < b, y(a) = α, y(b) = β. ?y(x) =1 ? 1 x2y(x) + 1 x3, y(1) = 1, y(2) =1 Exemplo: y00(x) = 2, x