Baixe o app para aproveitar ainda mais
Prévia do material em texto
UNIVERSIDADE FEDERAL DE MINAS GERAIS INSTITUTO DE CIEˆNCIAS EXATAS DEPARTAMENTO DE MATEMA´TICA EQUAC¸O˜ES DIFERENCIAIS A HORA´RIO: 11:10 a`s 12:50 - 09/09/2004 1a. Avaliac¸a˜o 1. Suponha que uma sala contendo 45 m3 de ar puro e´ invadida por fumantes que introduzem 0, 05 g/m3 de mono´xido de carbono a uma taxa de 0, 1 m3/min. Permite- se que a mistura (ar puro e mono´xido de carbono), bem circulada, saia da sala na mesma taxa de 0, 1 m3/min. (a) Encontre uma fo´rmula para a quantidade de mono´xido de carbono presente em cada instante t. (b) A concentrac¸a˜o de 1 g/m3 e´ atingida? (c) Fac¸a um esboc¸o do gra´fico da quantidade de mono´xido de carbono em func¸a˜o do tempo. Link para a soluc¸a˜o. 2. Considere o seguinte problema: dy dt = f(y) = (ey − 1)(y − 3), y(0) = y0. (1) (a) Esboce o gra´fico de f(y) em func¸a˜o de y, determine os pontos de equil´ıbrio e classifique cada um deles como assintoticamente esta´vel ou insta´vel. Justifique. (b) Esboce algumas soluc¸o˜es para diferentes valores de y0. Link para a soluc¸a˜o. 3. (a) Encontrar a soluc¸a˜o do problema de valor inicial dy dx = 2x− 2y y + 2x y(2) = −4 + √ 23 (b) Determinar o intervalo de validade da soluc¸a˜o. Link para a soluc¸a˜o. 1 Soluc¸a˜o 1. (a) O problema pode ser modelado pelo seguinte problema de valor inicial dQ dt = 0, 05 · 0, 1− 0, 1Q 45 . Q(0) = 0 dQ dt + 1 450 Q = 5 · 10−3 Multiplicando-se a equac¸a˜o pelo fator integrante µ(t) = e � 1 450 dt = e t 450 obtemos d dt ( e t 450Q ) = 5 · 10−3e t450 . Integrando-se obtemos e t 450Q(t) = 2, 25 e t 450 + C ⇒ Q(t) = 2, 25 + Ce− t450 . Substituindo t = 0 e Q = 0: C = −2, 25. Substituindo o valor de C encontrado: Q(t) = 2, 25− 2, 25e− t450 . (b) c(t) = Q(t) V = 0, 05− 0, 05e− t450 e nunca atingira´ 1 g/m3. (c) 0 500 1000 1500 0 0.5 1 1.5 2 2.5 t Q Link para a pro´xima questa˜o. 2 2. (a) −1 −0.5 0 0.5 1 1.5 2 2.5 3 3.5 4 −10 0 10 20 30 40 50 60 y y’=f(y) Os pontos de equil´ıbrio sa˜o as ra´ızes de f(y) = (ey − 1)(y − 3), ou seja, y1 = 0 e y2 = 3. i. y1 = 0 e´ ponto de equil´ıbrio esta´vel pois para valores de y pro´ximos de y1 = 0 temos • y′ = f(y) > 0, para y < y1 = 0 • y′ = f(y) < 0, para y > y1 = 0. O que implica que se y0 = y(0) e´ pro´ximo de y1 = 0 a soluc¸a˜o correspon- dente y(t) esta´ se aproximando de y1 = 0, quando t cresce. ii. y2 = 3 e´ ponto de equil´ıbrio insta´vel pois para valores de y pro´ximos de y2 = 3 temos • y′ = f(y) < 0, para y < y2 = 3 • y′ = f(y) > 0, para y > y2 = 3. O que implica que se y(0) e´ pro´ximo de y2 = 3 a soluc¸a˜o correspondente y(t) esta´ se afastando de y2 = 3, quando t cresce. 3 (b) −4 −3 −2 −1 0 1 2 3 4 −4 −3 −2 −1 0 1 2 3 4 t y Link para a pro´xima questa˜o. 4 3. (a) 2x− 2y + (−y − 2x)dy dx = 0 M = 2x− 2y, N = −y − 2x ∂M ∂y = −2, ∂N ∂x = −2 ∂M ∂y = ∂N ∂x ⇒ A equac¸a˜o e´ exata! ψ(x, y) = ∫ Mdx = x2 − 2xy + h(y) N = −y − 2x = ∂ψ ∂y = −2x+ h′(y) h′(y) = −y ⇒ h(y) = −1 2 y2 + C1 A soluc¸a˜o geral da equac¸a˜o e´ dada implicitamente por x2 − 2xy − 1 2 y2 = C Substituindo-se x = 2 e y = −4 + √ 23 na soluc¸a˜o acima 4− 4(−4 + √ 23)− 1 2 (−4 + √ 23)2 = 1 2 = C. Logo a soluc¸a˜o do problema de valor inicial e´ dada implicitamente por x2 − 2xy − 1 2 y2 = 1 2 (b) dx dy = 0 ⇔ y + 2x = 0 ⇔ y = −2x Substituindo-se na soluc¸a˜o x2 + 4x2 − 2x2 = 1 2 ⇔ 3x2 = 1 2 ⇔ x = ± √ 1 6 Como x = 2 > √ 1 6 , enta˜o a soluc¸a˜o e´ va´lida para x > √ 1 6 . 5
Compartilhar