Prévia do material em texto
DESCRIÇÃO Aplicação do conceito de equações diferenciais. PROPÓSITO Definir as equações diferenciais e resolver as equações diferenciais de primeira ordem. PREPARAÇÃO Antes de iniciar o conteúdo deste tema, tenha em mãos papel, caneta, uma calculadora científica ou a calculadora de seu smartphone ou computador. OBJETIVOS MÓDULO 1 Descrever os conceitos iniciais de uma equação diferencial MÓDULO 2 Classificar as equações diferenciais MÓDULO 3 Calcular equações diferenciais de primeira ordem MÓDULO 4 Reconhecer situações possíveis de serem modeladas por equações diferenciais de primeira ordem EQUAÇÕES DIFERENCIAIS DE PRIMEIRA ORDEM MÓDULO 1 Descrever os conceitos iniciais de uma equação diferencial CONCEITOS INICIAIS DA EQUAÇÃO DIFERENCIAL EQUAÇÕES DIFERENCIAIS Alguns problemas, em diversas áreas da Ciência e da Engenharia, podem ser modelados por uma equação que envolve uma variável e suas taxas de variação representadas por suas derivadas. Essas equações são denominadas equações diferenciais e têm uma função matemática como solução. As equações diferenciais representam um ramo importantíssimo da Matemática e tem diversas aplicações práticas. Neste módulo, serão apresentados os conceitos iniciais da equação diferencial. CONCEITOS INICIAIS A busca pela solução de um problema, no ramo da Engenharia ou da Ciência, pode ser feita por meio da modelagem de uma equação matemática. Em outras palavras, busca-se obter a solução do problema por meio da resolução de uma equação. Você já solucionou diversas vezes equações algébricas que envolviam variáveis e suas funções. No entanto, grande variedade desses problemas serão modelados por uma equação que representa o relacionamento entre a variável estudada e as suas taxas de variação. EQUAÇÃO ALGÉBRICA VERSUS EQUAÇÃO DIFERENCIAL AS EQUAÇÕES QUE RELACIONAM UMA VARIÁVEL E SUAS DERIVADAS SÃO DENOMINADAS EQUAÇÕES DIFERENCIAIS — EXPRESSÃO UTILIZADA DESDE 1676, QUANDO FOI CRIADA PELO MATEMÁTICO LEIBNIZ. A diferença entre uma equação algébrica e uma equação diferencial é que esta última envolve a(s) derivada(s) ou a(s) derivadas parciais de determinada variável. A VARIÁVEL PARA A QUAL SE BUSCA A SOLUÇÃO RECEBE O NOME DE VARIÁVEL DEPENDENTE OU INCÓGNITA DA EQUAÇÃO. As derivadas são obtidas pela taxa de variação da variável dependente em relação a uma ou mais variáveis independentes da equação. Lembre-se de que: Variável Independente: Entrada da equação Variável Dependente: Saída da equação, isto é, a solução desejada a ser obtida A solução de uma equação diferencial, caso exista, será uma função matemática que representa de que modo a variável estudada — dependente — irá depender de uma ou mais variáveis independentes. VARIÁVEL INDEPENDENTE VOCÊ SE LEMBRA DE QUANDO ESTUDOU CINEMÁTICA NA DISCIPLINA DE FÍSICA? Na ocasião, os problemas para os quais se desejava obter a posição de um objeto — representada pela variável s — em relação a uma variável independente — o tempo (t). Desse modo, era preciso obter a função s(t), que seria uma função matemática que solucionaria o modelo estudado. Em outras palavras, s(t) representa como a posição varia em função do tempo. Relacionávamos, para isso, a posição do objeto com as suas taxas de variação com o tempo (t), que eram a velocidade (v) — primeira derivada de s em função de t — e a aceleração (a) — segunda derivada de s em função de t. Portanto, modelávamos o problema por uma equação diferencial. EXEMPLO Poderia ser uma equação do tipo: � Atenção! Para visualização completa da equação utilize a rolagem horizontal com m, n e p sendo números reais. No entanto, s (t) = m + n v (t) + p a(t) v (t) = (t) = s′(t) ds dt e , de forma que: � Atenção! Para visualização completa da equação utilize a rolagem horizontal com m, n e p sendo números reais. Repare que a equação envolve a variável independente s e suas derivadas s’ e s’’, constituindo uma equação diferencial. Vejamos, a seguir, um exemplo com números. Nesse caso, um problema poderia ser modelado por: , PARA � Atenção! Para visualização completa da equação utilize a rolagem horizontal Na equação apresentada, os coeficientes m, n e p são constantes reais. No entanto, os coeficientes de uma equação diferencial também podem ser uma função matemática que dependa da variável independente. Por exemplo: , PARA � Atenção! Para visualização completa da equação utilize a rolagem horizontal Vimos, até aqui, que a equação diferencial envolve a variável e suas derivadas, e apresenta os seguintes elementos: Uma variável dependente Uma ou mais variáveis independentes E seus coeficientes TIPOS DE EQUAÇÕES DIFERENCIAIS a (t) = (t) = v′ (t) = s′′(t) dv dt s (t) = m + n s′ (t) + p s′′ (t) s (t) = 2 + s′ (t) + 4 s′′ (t) t ≥ 0 s (t) = e−t + 2s′ (t) − s′′ (t) 1 t + 1 t ≥ 0 Existem dois tipos de equações diferenciais — as equações diferenciais ordinárias (EDO) e as equações diferenciais parciais (EDP). Vamos conhecê-las a seguir: EQUAÇÕES DIFERENCIAIS ORDINÁRIAS – EDO As equações diferenciais podem envolver uma variável dependente que dependa apenas de uma variável independente, por exemplo . Observe que f é a variável dependente ou incógnita, x é a variável independente e os coeficientes dessa equação serão o 1, 2 e 4x. A variável dependente é normalmente representada por um símbolo, e não por uma função. Vejamos: Na qual y = f(x). Tal equação terá como solução uma função f(x) que depende apenas da variável x. Esse tipo de equação diferencial é denominado de equação diferencial ordinária ou (EDO). Vejamos outro exemplo de EDO já estudada em Física. Lembre-se da Mecânica e das Leis de Newton. Você estudou que a aceleração de um objeto, de massa m, dependia da força aplicada nele. Vamos imaginar uma força h conhecida que depende do tempo(t) aplicado a esse objeto. Com isso, teremos: como Conhecendo-se a função h(t), pode-se obter a função s(t), que representa a posição do objeto de massa m em relação ao tempo. Dessa forma, resolve-se a equação diferencial dada, por exemplo: , Repare que, na equação, aparece derivada de segunda ordem da variável s, constituindo uma equação diferencial. Além disso, existe apenas uma variável independente — a variável t —, de modo que se trata de uma EDO. EQUAÇÕES DIFERENCIAIS PARCIAIS (EDP) Outro tipo de equação diferencial é quando a variável dependente depende de mais de uma variável independente. Nesse caso, na equação, aparecerão derivadas parciais de diversas ordens. A solução da equação será uma função escalar que dependerá das variáveis envolvidas. A equação, a seguir, muito utilizada no eletromagnetismo, exemplifica esse tipo de equação: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Repare que a variável f depende das variáveis x,y e z. A solução, então, será uma função f(x,y,z), que apresenta a dependência de f com as variáveis independentes. Trata-se do tipo de equação denominada equação diferencial parcial ou EDP. Nos próximos módulos, estudaremos apenas as equações do tipo ordinárias (EDO). Porém, é fundamental saber reconhecer se uma equação é diferencial ou não, e de que tipo. f (x) + 2f ′ (x) = 4x y + 2y′ = 4x h (t) = ma(t) a (t) = s′′(t) s′′ (t) = h(t) 1 m s′′ = 4 cos(2t) t ≥ 0 (x, y, z) + (x, y, z) + (x, y, z) = 0 ∂2f ∂x2 ∂2f ∂y2 ∂2f ∂z2 javascript:void(0) javascript:void(0) CLASSIFIQUE AS EQUAÇÕES DIFERENCIAIS, A SEGUIR, EM ALGÉBRICA, EDP OU EDO: A) B) C) D) RESOLUÇÃO A equação da letra A é uma equação diferencial ordinária (EDO), pois apresenta a variável y relacionada com as suas derivadas e com a variável x. Por isso, y é a variável dependente que depende apenas de uma variável independente x. A equação da letra B é uma equação algébrica, pois não apresenta nenhuma derivada. Nessa equação, definimos y dependendo de x, ou x dependendo de y. Depende de qual das duas é conhecida no problema. A equação da letra C é uma equação diferencial, pois envolve derivadas em seustermos. Como as derivadas que aparecem são derivadas parciais, será uma equação diferencial parcial (EDP). Nesse caso, z será a variável dependente, e x e y serão as variáveis independentes. A equação da letra D, por fim, também é uma EDO, pois relaciona as derivadas da variável dependente m com a variável independente t. SOLUÇÃO DE UMA EQUAÇÃO DIFERENCIAL COMO OBTER A SOLUÇÃO DE UMA EQUAÇÃO DIFERENCIAL? RESPOSTA Não existe um método único para resolver todos os tipos de equações diferenciais. Há, no entanto, alguns métodos de grande abrangência, isto é, que resolvem grande número de equações (veremos esses métodos em módulos posteriores). LEMBRE-SE DE QUE A SOLUÇÃO DE UMA EQUAÇÃO DIFERENCIAL, SE EXISTIR, SERÁ UMA FUNÇÃO QUE ATENDE A EQUAÇÃO APRESENTADA. A solução será uma ou mais famílias de funções. A família de funções que atende a equação diferencial é denominada solução geral da equação diferencial. Para determinar uma função específica, isto é, uma solução particular, são necessárias algumas informações adicionais, que denominaremos condições iniciais ou condições de contorno. A descoberta de uma solução particular, às vezes, é denominada problema de valor inicial. 3y + y′ − 2 cos (x) = 4 3y + x2– cos(x) = 2 z − 2 + = x + y ∂z ∂x ∂2z ∂x∂y = 4 + 2t − d2m dt2 dm dt EXEMPLO Seja a equação Para essa equação, uma possível solução seria s(t) = cos (2t). Ainda iremos aprender a determinar a solução dessa equação, mas já podemos verificar se a função s(t) apresentada é ou não solução da equação diferencial. Como sabemos, a solução será a função que satisfará a equação dada, de modo que: Se, então , então e Substituindo, verifica-se que , satisfazendo a EDO apresentada. Dessa maneira, s(t) = cos (2t) é solução da equação diferencial. REPARE QUE NÃO APENAS COS (2T) SERÁ SOLUÇÃO, MAS TODA FAMÍLIA DE FUNÇÃO DO TIPO , COM K REAL. OBTENHA A PRIMEIRA E SEGUNDA DERIVADAS E VEJA QUE A FAMÍLIA DE FUNÇÕES SATISFAZ A EQUAÇÃO. ESSA FAMÍLIA DE FUNÇÕES SERÁ DENOMINADA SOLUÇÃO GERAL DA EQUAÇÃO DIFERENCIAL. Para um valor de k, por exemplo k = 3, tendo-se s(t) = cos (2t) + 3, temos uma solução particular. A própria função s(t) = cos (2t) é uma outra solução particular. Vamos testar, agora, se a função é ou não solução da equação diferencial dada. Se , então e . Assim, não satisfaz a equação, já que: � Atenção! Para visualização completa da equação utilize a rolagem horizontal s′′ = s′ + 2 sen (2t) − 4 cos (2t) s(t) = cos(2t) s′(t) = −2sen(2t) s′′(t) = −4 cos(2t) s′′ = s′ + 2 sen (2t) − 4 cos (2t) s(t) = cos(2t) + k s(t) = t2 + 2 s(t) = t2 + 2 s′ (t) = 2t s′′ (t) = 2 s′′ ≠ s′ + 2 sen (2t) − 4 cos (2t) Logo, não é a solução da equação dada. VERIFIQUE SE A FUNÇÃO , COM K REAL, É UMA SOLUÇÃO DA EQUAÇÃO DIFERENCIAL . RESOLUÇÃO Usando as regras de derivação, temos: se , então: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Substituindo, temos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Comprovando que é solução geral para equação diferencial dada. DETERMINE A SOLUÇÃO PARTICULAR QUE ATENDE A EQUAÇÃO DIFERENCIAL y = cosx (senx − k) y′ + tg (x) y − cos2x = 0 y = cosx (senx − k) y′ = (−senx) (senx − k) + cosx (cosx) = −sen2x + ksenx + cos2x y′ + tg (x) y − cos2x = −sen2x + ksenx + cos2x + tgx cosx (senx − k) − cos2x = = −sen2x + ksenx + cos2x + senx (senx − k) − cos2x = = −sen2x + ksenx + cos2x + sen2x − ksenx − cos2x = 0 y = cosx (senx − k) y′ + tg (x) y − cos2x = 0 E A CONDIÇÃO INICIAL DE PARA . RESOLUÇÃO No exemplo anterior, vimos que a solução geral da equação será , com k real. Quando . � Atenção! Para visualização completa da equação utilize a rolagem horizontal Porém, para x = 0 temos y = – 1, portanto – 1 = – k, de modo que k = 1. Assim, a solução particular será . VERIFIQUE SE A FUNÇÃO ESCALAR É UMA SOLUÇÃO PARTICULAR DA EQUAÇÃO DIFERENCIAL . RESOLUÇÃO Precisamos testar se a função u(x,t) dada satisfaz a equação. Repare que, agora, temos uma EDP. Precisamos lembrar de como fazer derivadas parciais. Derivamos em função da variável mantendo todas as demais como constantes. SE � Atenção! Para visualização completa da equação utilize a rolagem horizontal Agora, obteremos a derivada parcial de segunda ordem y =– 1 x = 0 y = cosx (senx − k) x = 0 → y = cos0 (sen0 − k) = 1 (0 − k) = −k y = cosx (senx − 1) u (x, t) = e−4tcos(2x) − = 0 ∂2u ∂x2 ∂u ∂t u (x, t) = e−4tcos (2x) → = (−4) e−4tcos (2x) ∂u ∂t ∂2u ∂x2 SE � Atenção! Para visualização completa da equação utilize a rolagem horizontal SE � Atenção! Para visualização completa da equação utilize a rolagem horizontal Substituindo, temos: SE � Atenção! Para visualização completa da equação utilize a rolagem horizontal Assim, satisfaz a equação, provando que é uma solução da equação diferencial. TEORIA NA PRÁTICA Seja um objeto de massa m, medida em kg, presa na extremidade de uma mola com constante elástica k, medida em N/m. Considere que o ponto de equilíbrio da mola se encontre na origem, isto é, x = 0. Quando comprimimos ou esticamos a mola de x, medido em metros, o objeto ficará sujeito a uma força elástica dada por Pela Lei de Newton, essa força será igual à massa do objeto vezes a aceleração. Com isso, teremos: SE u (x, t) = e−4tcos (2x) → = (−2) e−4tsen (2x) ∂u ∂x = (−2) e−4tsen (2x) → = (−2) (2) e−4tcos (2x) = (−4)e−4tcos (2x) ∂u ∂x ∂2u ∂x2 − = (−4) e−4tcos (2x) − (−4) e−4tcos (2x) = 0 ∂2u ∂x2 ∂u ∂t (x, t) = e−4tcos (2x) →F = −kx. →F = m→a = m d2x dt2 � Atenção! Para visualização completa da equação utilize a rolagem horizontal Assim, Verifique se a função , C real, é solução da equação diferencial que representa o movimento de um objeto de 5kg preso em uma mola de constante elástica 400 N/m. Determine a solução particular sabendo que para t = 0 a posição do objeto vale x = 0. RESOLUÇÃO SOLUÇÃO DE EQUAÇÃO DIFERENCIAL MÃO NA MASSA 1. MARQUE A ALTERNATIVA QUE APRESENTA UMA EQUAÇÃO QUE NÃO É UMA EQUAÇÃO DIFERENCIAL: A) � Atenção! Para visualização completa da equação utilize a rolagem horizontal B) � Atenção! Para visualização completa da equação utilize a rolagem horizontal C) � Atenção! Para visualização completa da equação utilize a rolagem horizontal D) � Atenção! Para visualização completa da equação utilize a rolagem horizontal E) m = −kx d2x dt2 x(t) = sen (√80 t + C) = (3 − y) dy dx 4 x + = cos(x + y) ∂2v ∂x2 ∂2v ∂y2 s′′ + cos (x)s′ = 4x u + ln v = (u + v) � Atenção! Para visualização completa da equação utilize a rolagem horizontal 2. MARQUE A ALTERNATIVA QUE APRESENTA UMA EDP: A) � Atenção! Para visualização completa da equação utilize a rolagem horizontal B) � Atenção! Para visualização completa da equação utilize a rolagem horizontal C) � Atenção! Para visualização completa da equação utilize a rolagem horizontal D) � Atenção! Para visualização completa da equação utilize a rolagem horizontal E) � Atenção! Para visualização completa da equação utilize a rolagem horizontal 3. SEJA A EQUAÇÃO DIFERENCIAL . MARQUE A ALTERNATIVA QUE APRESENTA UMA SENTENÇA ERRADA EM RELAÇÃO A ESSA EQUAÇÃO A) Trata-se de uma equação diferencial ordinária. B) A variável independente é a variável u. C) A variável dependente é a variável v. D) A família de funções , com k real, é uma solução geral dessa equação. E) A função é uma solução particular dessa equação. = −2x d2x dt2 4 = 3t − 2x dx dt + − ysen(x + y + z) = 0 ∂2y ∂x2 ∂y ∂z m′′ + ln (n)m = 4 x + y + z = sen(x + y + z) cos (v) = 2u v − 1 u − 3v = u2 dv du v (u) = k u3 − u2 v (u) = u2 4. SABENDO QUE UMA DAS SOLUÇÕES GERAIS DA EQUAÇÃO DIFERENCIAL É A FUNÇÃO , COM K REAL, DETERMINE A SOLUÇÃO PARTICULAR DESSA FAMÍLIA DE SOLUÇÃO QUE ATENDA A CONDIÇÃO INICIAL S(0) = 2. A) B) C) D) E) 5. MARQUE A ALTERNATIVA QUE APRESENTA UMA SOLUÇÃO GERAL PARA A EQUAÇÃO DIFERENCIAL : A) , com k real � Atenção! Para visualização completada equação utilize a rolagem horizontal B) , com k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal C) , com k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal D) , com k real s′ − s + u = 0 s (u) = 1 − keu + u s (u) = 2u + 1 s (u) = 1 + eu + u s (u) = 1 − eu + u s (u) = eu + u s (u) = 1 + 3eu + u y − 2x + xy′ = 0 y = cos(2x) + 3k y = 3k ln(x) y = x − k x y = k 2x � Atenção! Para visualização completa da equação utilize a rolagem horizontal E) , com k � Atenção! Para visualização completa da equação utilize a rolagem horizontal 6. SABE-SE QUE A FAMÍLIA DE FUNÇÕES , M E N REAIS ATENÇÃO! PARA VISUALIZAÇÃO COMPLETA DA EQUAÇÃO UTILIZE A ROLAGEM HORIZONTAL É UMA SOLUÇÃO GERAL PARA EQUAÇÃO DIFERENCIAL . SABENDO QUE V = 1, PARA QUANDO , E QUE V = 4, PARA QUANDO , DETERMINE A ALTERNATIVA QUE APRESENTA UMA SOLUÇÃO PARTICULAR DESSA SOLUÇÃO GERAL QUE ATENDA AS CONDIÇÕES DE CONTORNO DADAS. A) � Atenção! Para visualização completa da equação utilize a rolagem horizontal B) � Atenção! Para visualização completa da equação utilize a rolagem horizontal C) � Atenção! Para visualização completa da equação utilize a rolagem horizontal D) � Atenção! Para visualização completa da equação utilize a rolagem horizontal E) � Atenção! Para visualização completa da equação utilize a rolagem horizontal y = kx + 1 x v (x, y) = m ln (x2 + y2) + n + = 0 ∂2v ∂x2 ∂2v ∂y2 x2 + y2 = 1 x2 + y2 = e v (x, y) = 2 ln (x2 + y2) v (x, y) = 2 exp (x2 + y2) v (x, y) = 3 ln (x2 + y2) + 1 v (x, y) = ln (x2 + y2) v (x, y) = 2 ln (x2 + y2) − 1 GABARITO 1. Marque a alternativa que apresenta uma equação que não é uma equação diferencial: A alternativa "D " está correta. Como já vimos, uma equação diferencial é aquela que envolve as derivadas da variável analisada. Na equação da letra A, verifica-se que a variável dependente é y e a variável independente é x. Na equação, aparece a derivada de primeira ordem de y em relação a x, sendo, então, uma EDO. Na equação da letra B, verifica-se que a variável dependente é v, e as variáveis independentes são x e y. Na equação, aparecem as derivadas parciais de segunda ordem de v, sendo, portanto, uma EDP. Na equação da letra C, verifica-se que a variável dependente é s e a variável independente é x. Na equação, aparece a derivada de primeira e segunda ordem de s em relação a x, sendo, assim, uma EDO. Na equação da letra D, não aparece nenhuma derivada, de modo que se trata de uma equação algébrica que envolve as variáveis u e v. Sendo assim, a resposta da questão. Na equação da letra E, verifica-se que a variável dependente é x e a variável independente é t. Na equação, aparece a derivada de segunda ordem, de x em relação a t, sendo, portanto, uma EDO. 2. Marque a alternativa que apresenta uma EDP: A alternativa "B " está correta. Como já estudamos, uma equação diferencial é aquela que envolve as derivadas da variável analisada. EDP é a equação diferencial na qual a variável dependente depende de mais de uma variável independente. Repara que apenas a equação da letra B apresenta uma variável y que depende de duas variáveis x e z. Como aparecem, nessa equação, as derivadas parciais, trata-se de uma equação diferencial parcial, sendo a resposta da questão. As equações da letra A e da letra C são EDOs. As equações da letra D e letra E são equações algébricas. Atenção, apesar da equação da letra D apresentar uma variável que depende de outras duas, ela não tem nenhuma derivada entre seus termos, não sendo uma equação diferencial. 3. Seja a equação diferencial . Marque a alternativa que apresenta uma sentença ERRADA em relação a essa equação A alternativa "E " está correta. Veja a solução no vídeo abaixo: 4. Sabendo que uma das soluções gerais da equação diferencial é a função , com k real, determine a solução particular dessa família de solução que atenda a condição inicial s(0) = 2. A alternativa "B " está correta. Inicialmente, como exercício, você pode testar e confirmar que a solução geral satisfaz a equação diferencial dada. Vejamos: Se u − 3v = u2 dv du s′ − s + u = 0 s (u) = 1 − keu + u � Atenção! Para visualização completa da equação utilize a rolagem horizontal Desse modo: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Para que a solução particular seja parte da família da solução geral, ela deve ser obtida por meio de um valor da constante k. Observe que as equações das letras A e E não podem ser obtidas da solução geral . Isto é, para nenhum valor de k é possível partir da equação geral e obter essas equações. Agora, vamos analisar qual a solução particular da família em que se tem s = 2 quando u = 0. � Atenção! Para visualização completa da equação utilize a rolagem horizontal Desse modo, Logo, a solução particular será � Atenção! Para visualização completa da equação utilize a rolagem horizontal Que está na letra B. 5. Marque a alternativa que apresenta uma solução geral para a equação diferencial : A alternativa "C " está correta. É preciso testar todas as alternativas e verificar se satisfazem ou não a equação apresentada no enunciado. Usando as regras de derivação, temos: LETRA A Se � Atenção! Para visualização completa da equação utilize a rolagem horizontal Assim, � Atenção! Para visualização completa da equação utilize a rolagem horizontal Não é a solução LETRA B Se s (u) = 1 − keu + u → s′ = −keu + 1 s′ − s + u = −keu + 1 − (1 − keu + u) + u = 0 s (u) = 1 − keu + u s (u) = 1 − keu + u s (u) = 1 − keu + u → s (0) = 1 − ke0 + 0 = 1 − k = 2 k = 1– 2 =– 1 s (u) = 1 − (−1) eu + u = 1 + eu + u y − 2x + xy′ = 0 y = cos(2x) + 3k → y′ = −2sen(2x) y − 2x + xy′ = cos (2x) + 3k − 2x + x(−2sen (2x)) ≠ 0 y = 3kln(x) → y′ = 3k x � Atenção! Para visualização completa da equação utilize a rolagem horizontal Assim, � Atenção! Para visualização completa da equação utilize a rolagem horizontal Não é a solução LETRA C Se � Atenção! Para visualização completa da equação utilize a rolagem horizontal Então � Atenção! Para visualização completa da equação utilize a rolagem horizontal Assim, � Atenção! Para visualização completa da equação utilize a rolagem horizontal Desse modo, é uma solução geral da equação dada, sendo a resposta da questão. LETRA D Se � Atenção! Para visualização completa da equação utilize a rolagem horizontal Então � Atenção! Para visualização completa da equação utilize a rolagem horizontal Assim, � Atenção! Para visualização completa da equação utilize a rolagem horizontal Não é a solução LETRA E Se � Atenção! Para visualização completa da equação utilize a rolagem horizontal y − 2x + xy′ = 3k ln(x) − 2x + x ≠ 0 3k x y = x − k x y′ = 1 − k (−1) = 1 + 1 x2 k x2 y − 2x + xy′ = x − − 2x + x(1 + ) = x − − 2x + x + = 0k x k x2 k x k x y = x − k x y = k 2x y′ = k (−1) 1 2x2 y − 2x + xy′ = − 2x + x( ) = − 2x − = −2x ≠ 0k 2x −k 2x2 k 2x k 2x y = kx + 1 x Então � Atenção! Para visualização completa da equação utilize a rolagem horizontal Assim, � Atenção! Para visualização completa da equação utilize a rolagem horizontal Não é a solução 6. Sabe-se que a família de funções , m e n reais Atenção! Para visualização completa da equação utilize a rolagem horizontal é uma solução geral para equação diferencial . Sabendo que v = 1, para quando , e que v = 4, para quando , determine a alternativa que apresenta uma solução particular dessa solução geral que atenda as condições de contorno dadas. A alternativa "C " está correta. Inicialmente, como exercício, podemos testar e confirmar que a solução geral satisfaz a equação diferencial dada. Para que a solução particular seja parte da família da solução geral, ela deve ser obtida por um valor da constante k. Com isso, eliminamos a alternativa B. Precisamos testar se a função v(x,y) dada satisfaz a equação. Agora, vamos analisar quala solução particular da família em que se tem v = 1 para quando e v = 4 para quando � Atenção! Para visualização completa da equação utilize a rolagem horizontal Quando � Atenção! Para visualização completa da equação utilize a rolagem horizontal Portanto � Atenção! Para visualização completa da equação utilize a rolagem horizontal y′ = k − 1 x2 y − 2x + xy′ = kx + − 2x + x(k − ) = kx + − 2x + kx − ≠ 0k x k x2 k x k x v (x, y) = m ln (x2 + y2) + n + = 0 ∂2v ∂x2 ∂2v ∂y2 x2 + y2 = 1 x2 + y2 = e v (x, y) = m ln (x2 + y2) + n x2 + y2 = 1 x2 + y2 = e v (x, y) = m ln (x2 + y2) + n x2 + y2 = 1 → v = m ln (1) + n = m.0 + n = n = 1 n = 1e v (x, y) = m ln (x2 + y2) + 1 Quando � Atenção! Para visualização completa da equação utilize a rolagem horizontal Assim, m=3 Logo, a resposta está na letra C, com . GABARITO VERIFICANDO O APRENDIZADO 1. MARQUE A ALTERNATIVA QUE APRESENTA UMA EQUAÇÃO QUE NÃO É UMA EQUAÇÃO DIFERENCIAL: A) B) C) D) E) 2. MARQUE A ALTERNATIVA QUE APRESENTA UMA SOLUÇÃO PARA EQUAÇÃO DIFERENCIAL . A) B) C) x2 + y2 = e → v = m ln (e) + 1 = m.1 + 1 = m + 1 = 4 v (x, y) = 3 ln (x2 + y2) + 1 y = 3x − y dx dy s′′ + 2s = 3u − 2 = 3t + 1 + 2s ∂v ∂t 3x + y = 2z − 1 − 2v = u ∂2x ∂u2 ∂2x ∂v2 y′′ − 2y′ + y = 0 y = xex y = 2 cos(x) D) E) GABARITO 1. Marque a alternativa que apresenta uma equação que NÃO é uma equação diferencial: A alternativa "D " está correta. Uma equação diferencial é aquela que envolve as derivadas da variável analisada. Na equação da letra A, verifica-se que a variável dependente é x e a variável independente é y. Na equação, aparece a derivada de primeira ordem de x em relação a y, sendo, então, uma EDO. Na equação da letra B, verifica-se que a variável dependente é s e a variável independente é u. Na equação, aparece a derivada de segunda ordem de s em relação a u, sendo também uma EDO. Na equação da letra C, verifica-se que a variável dependente é v e as variáveis independentes são t e s. Na equação, aparece a derivada parcial de v em relação a t, sendo, assim, uma EDP. A equação da letra D não apresenta nenhuma derivada, de modo que a resposta da questão é uma equação algébrica. Na equação da letra E, por fim, verifica-se que a variável dependente é x e as variáveis independentes são u e v. A equação apresenta derivadas parciais de x em relação a u e v, sendo uma EDP. Portanto, a resposta correta é a letra D. 2. Marque a alternativa que apresenta uma solução para equação diferencial . A alternativa "A " está correta. Para verificar qual a função é solução da equação, necessitamos verificar se a função satisfaz a equação. LETRA A Se � Atenção! Para visualização completa da equação utilize a rolagem horizontal Substituindo, temos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Trata-se, dessa forma, da solução da equação e a resposta da questão. LETRA B Se � Atenção! Para visualização completa da equação utilize a rolagem horizontal Substituindo, temos: y = xe−x y = sen(x) y = e−x y′′ − 2y′ + y = 0 y = xex → y′ = ex + xex → y′′ = ex + ex + xex = 2ex + xex y′′ − 2y′ + y = 2ex + xex − 2 (ex + xex) + xex = 0 y = 2 cos (x) → y′ = −2 sen (x) → y′′ = −2 cos(x) � Atenção! Para visualização completa da equação utilize a rolagem horizontal Não é a solução da equação. LETRA C Se � Atenção! Para visualização completa da equação utilize a rolagem horizontal Substituindo, temos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Não é a solução da equação. LETRA D Se � Atenção! Para visualização completa da equação utilize a rolagem horizontal Substituindo, temos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Não é a solução da equação. LETRA E Se � Atenção! Para visualização completa da equação utilize a rolagem horizontal Substituindo, temos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Não é a solução da equação. MÓDULO 2 Classificar as equações diferenciais CLASSIFICAÇÃO DAS EQUAÇÕES DIFERENCIAIS y′′ − 2y′ + y = −2 cos(x) − 2(−2 sen(x)) + 2 cos(x) = −4 sen(x) ≠ 0 y = xe−x → y′ = e−x − xe−x → y′′ = −e−x − e−x + xe−x = −2e−x + xe−x y′′ − 2y′ + y = −2e−x + xe−x − 2 (e−x − xex) + xe−x ≠ 0 y = sen (x) → y′ = cos(x) → y′′ = − sen(x) y′′ − 2y′ + y = −sen(x) − 2 cos(x) + sen(x) = −2 cos(x) ≠ 0 y = e−x → y′ = −e−x → y′′ = e−x y′′ − 2y′ + y = e−x − 2(−e −x ) + e−x = 4e−x TIPOS DE CLASSIFICAÇÃO DAS EQUAÇÕES DIFERENCIAIS No módulo anterior, vimos os conceitos iniciais da equação diferencial, e conhecemos as equações diferenciais ordinária e parcial. Neste módulo, veremos que as equações diferenciais podem ser classificadas de diversas formas, e conheceremos algumas delas, principalmente a classificação quanto à ordem e ao grau. VOCÊ JÁ SABE DISTINGUIR UMA EQUAÇÃO DIFERENCIAL ORDINÁRIA DE UMA EQUAÇÃO DIFERENCIAL PARCIAL, CERTO? Equação ordinária A incógnita depende apenas de uma variável independente, e aparecem, na equação, as derivadas dessa incógnita em relação à variável independente, em suas diversas ordens. Equação parcial A incógnita a ser descoberta depende de duas ou mais variáveis independentes e, na equação, aparecem derivadas parciais dessa incógnita em relação às variáveis independentes. ATENÇÃO Lembre-se de que a solução da EDO — caso exista — é uma função real, e solução da EDP — caso exista — é uma função escalar. As equações diferenciais podem ser classificadas de diversas formas. A seguir, veremos a classificação da equação diferencial quanto à ordem, ao grau e à linearidade. ORDEM DA EQUAÇÃO DIFERENCIAL As derivadas ou derivadas parciais que aparecem na equação diferencial podem ser de diversas ordens: DERIVAÇÃO DE ORDEM SUPERIOR A derivada de uma função também é uma função. Dessa forma, a função derivada também pode possuir uma derivada. Nesse sentido, a derivação de uma função derivada é denominada derivação de ordem superior. A ordem de uma derivada corresponde ao número de vezes que derivamos a função. Suponha a função y = x6 + 5 A primeira derivada de y, ou derivada de primeira ordem, será A segunda derivada de y, ou derivada de segunda ordem, será a derivada da primeira derivada, de modo que: , e assim sucessivamente. A derivada de ordem n, portanto, será a derivada da derivada de ordem n – 1 da função. Representamos a derivada de ordem superior n – n inteiro positivo – por f(n)(x) ou D(n)f(x). Utilizando a notação de Leibniz, representaremos a derivada de ordem n por: � Atenção! Para visualização completa da equação utilize a rolagem horizontal As derivadas parciais de ordem superior seguem raciocínio semelhante. DERIVAÇÃO PARCIAL DE ORDEM SUPERIOR A derivada parcial de uma função escalar também é uma função escalar. Por serem funções escalares, também podemos determinar as suas derivadas parciais em relação às variáveis independentes. A derivada parcial de uma função que já é derivada parcial de uma função é denominada derivada parcial de segunda ordem. Se repetirmos o processo, teremos as derivadas parciais de terceira, quarta, quinta, até a enésima ordem. Essas derivadas parciais são conhecidas como derivadas parciais de ordem superior. Por exemplo, seja , então: E � Atenção! Para visualização completa da equação utilize a rolagem horizontal Vamos, agora, determinar as derivadas parciais de segunda ordem, isto é, a derivada parcial da função escalar y′ = = 6x5. dy dx y′′ = = = 30x4 dy′ dx d2y dx2 = ( ) = ( ( )) dny dx n d d dn−1y dxn−1 d dx d dx dn−2y dxn−2 f(x, y) = 8x2y3 fx (x, y) = 16xy 3 fy (x, y) = 24x 2y2 fx (x, y) = 16xy 3 � Atenção! Para visualização completa da equação utilize a rolagem horizontal Relembre a notação: OU OU � Atenção! Para visualização completa da equação utilize a rolagem horizontal As funções e são denominadas derivadas parciais de primeira ordem da função . Já as funções , fx (x, y) = 16xy 3 →= 16y3 ∂fx ∂x fx (x, y) = 16xy 3 → = 48xy2 ∂fx ∂y = ( ) = = 16y3∂fx ∂x ∂ ∂x ∂f ∂x ∂2f ∂x2 (fx)x = fxx = 16y 3 = ( ) = = 48xy2∂fx ∂y ∂ ∂y ∂f ∂x ∂2f ∂y∂x (fx)y = fxy = 48xy 2 fx(x, y) fy(x, y) f(x, y) fxx(x, y) fxy(x, y) , e são as derivadas de segunda ordem da função f(x,y). Podemos repetir esse processo sucessivamente. SAIBA MAIS A ordem de uma equação diferencial será a ordem da mais alta derivada ou da derivada parcial da função incógnita que aparece na equação. Para entender o conceito na prática, vamos estudar os seguintes exemplos: DETERMINE A ORDEM DAS SEGUINTES EQUAÇÕES DIFERENCIAIS: A) � Atenção! Para visualização completa da equação utilize a rolagem horizontal B) � Atenção! Para visualização completa da equação utilize a rolagem horizontal C) � Atenção! Para visualização completa da equação utilize a rolagem horizontal D) fyx(x, y) fyy(x, y) s′ − 2s + 5x4 = cosx z2 − − 4 = − ln(z) d3y dz d2y dz2 − xyz = 2 ∂3m ∂z3 uvw − = ∂u ∂v ∂2u ∂v∂w � Atenção! Para visualização completa da equação utilize a rolagem horizontal RESOLUÇÃO a) Temos, aqui, uma equação diferencial ordinária na qual a função incógnita s — ou variável dependente — está relacionada apenas com uma variável dependente x. Analisando a equação, observa-se que a derivada de mais alta ordem é a primeira derivada, de forma que a ordem dessa EDO é 1. b) Nesse caso, temos uma equação diferencial ordinária na qual a função incógnita y está relacionada apenas com uma variável dependente z. A equação apresenta mais de uma derivada, mas a derivada de mais alta ordem é a terceira derivada de y em função de z. Com isso, a ordem dessa EDO é 3. c) Aqui temos uma equação diferencial parcial, na qual a variável dependente m depende das variáveis x,y e z. Observe que a equação e a derivada parcial de mais alta ordem é uma derivada parcial de terceira ordem. Desse modo, essa EDP apresenta ordem 3 d) A equação da letra D é novamente uma equação diferencial parcial em que a variável u depende das variáveis v e w. Nessa equação, temos uma derivada parcial de primeira ordem e uma derivada parcial de ordem 2. Desse modo, a EDP terá ordem 2. FORMA PADRÃO Toda vez que o coeficiente que multiplica a derivada de mais alta ordem da equação diferencial for um, dizemos que a equação diferencial está na sua forma padrão. No exemplo anterior, as equações diferenciais da letra A e da letra C estão na forma padrão. Já as equações da letra B e D não estão, mas podem ser colocadas. Vejamos: , multiplicando-se todos os coeficientes por ( – 1). � Atenção! Para visualização completa da equação utilize a rolagem horizontal , dividindo-se todos os coeficientes por 3. � Atenção! Para visualização completa da equação utilize a rolagem horizontal AS EDOS SEMPRE PODEM SER COLOCADAS NA SUA FORMA PADRÃO. GRAU DE UMA EQUAÇÃO DIFERENCIAL s′ − 2s + 5x4 = cosx z2 − − 4 = − ln(z) d3y dz d2y dz2 − xyz = 2 ∂3m ∂z3 uvw − 6 = 3 ∂u ∂v ∂2u ∂v∂w z2 − − 4 = − ln (z) → −z2 + + 4 = − + ln(z) d3y dz d2y dz2 d3y dz d2y dz2 uvw − 6 = 3 → uvw − 2 = ∂u ∂v ∂2u ∂v∂w 1 3 ∂u ∂v ∂2u ∂v∂w Uma equação diferencial também é classificada quanto ao seu grau. Nesse caso, o grau será o expoente da derivada mais alta que existe na equação diferencial. ATENÇÃO É preciso observar a derivada de mais alta ordem para definir a ordem da equação diferencial. Para isso, basta analisar essa derivada e ver o expoente a que ela está submetida. Dessa forma, iremos definir seu grau. Se olharmos o exemplo anterior, veremos que todas as equações diferenciais apresentadas têm grau 1, já que as derivadas de maior ordem estão sempre elevadas ao expoente um. Para entender melhor, vamos estudar alguns exemplos. DETERMINE A ORDEM E O GRAU DAS EQUAÇÕES DIFERENCIAIS APRESENTADAS A) � Atenção! Para visualização completa da equação utilize a rolagem horizontal B) � Atenção! Para visualização completa da equação utilize a rolagem horizontal C) � Atenção! Para visualização completa da equação utilize a rolagem horizontal RESOLUÇÃO a) A equação diferencial apresentada é uma equação diferencial ordinária na qual a função incógnita y depende da variável dependente x. Analisando a equação, observa-se que a derivada de mais alta ordem é a de segunda derivada, de modo que a ordem dessa EDO é 2. A derivada está elevada ao expoente três, portanto, o grau dessa EDO vale 3. b) y′ − 3x(y′′) 3 − x4 = 3 ( ) 3 − ( ) 2 = s5 ds du d3s du3 ( ) 2 + = 3 ∂3z ∂u∂v∂w ∂2z ∂w2 y′ − 3x(y′′) 3 − x4 = 3 Nesse caso, também temos uma equação diferencial ordinária, cuja variável dependente s depende da variável u. Observe que derivada de mais alta ordem é uma derivada de terceira ordem e está elevada ao expoente 2. Com isso, essa EDO é de ordem 3 e grau 2. c) Por fim, temos uma equação diferencial parcial, com a incógnita z dependendo das variáveis independentes u,v e w. A derivada de maior grau é a derivada de ordem 3. Como o termo está elevado ao expoente 2, a EDP tem ordem 3 e grau 2. LINEARIDADE DA EQUAÇÃO DIFERENCIALO Outra classificação para as equações diferenciais diz respeito à sua linearidade, isto é, se a equação diferencial é linear ou não linear. A equação diferencial linear ocorre nos seguintes casos: A variável dependente e suas derivadas só podem aparecer na forma simples, isto é, elevadas ao expoente um. Os coeficientes da equação diferencial, isto é, os termos que multiplicam a incógnita ou suas derivadas, só podem depender da(s) variável(is) independente(s) ou serem números reais. Já a equação diferencial não linear ocorre nas seguintes situações: A variável dependente e suas derivadas aparecerem elevadas a um expoente numérico diferente de 1. Se aparecer algum produto entre a variável dependente e suas derivadas, ou das derivadas das incógnitas entre si. Se aparecer, como coeficiente, uma função da variável dependente ou de suas derivadas. É importante enfatizar que as variáveis independentes, nas equações diferenciais lineares, podem aparecer como funções ou com expoentes diferentes da unidade. ATENÇÃO Como a incógnita e suas derivadas estão sempre elevadas ao expoente 1, toda equação diferencial linear, obrigatoriamente, apresentará um grau 1. Vamos estudar os exemplos a seguir para compreender melhor. DETERMINE SE AS EQUAÇÕES DIFERENCIAIS A SEGUIR SÃO OU NÃO LINEARES. A) ( ) 3 − ( ) 2 = s5 ds du d3s du3 ( ) 2 + = 3 ∂3z ∂u∂v∂w ∂2z ∂w2 � Atenção! Para visualização completa da equação utilize a rolagem horizontal B) � Atenção! Para visualização completa da equação utilize a rolagem horizontal C) � Atenção! Para visualização completa da equação utilize a rolagem horizontal D) � Atenção! Para visualização completa da equação utilize a rolagem horizontal RESOLUÇÃO a) � Atenção! Para visualização completa da equação utilize a rolagem horizontal A equação diferencial é uma EDO. Apesar de os coeficientes dependerem apenas da variável independente x ou serem número, na equação, tem-se uma derivada da incógnita y com expoente 3. Desse modo, a equação é não linear. b) � Atenção! Para visualização completa da equação utilize a rolagem horizontal A equação diferencial é uma EDO. Os coeficientes dependem apenas da variável independente u ou são números. O expoente da incógnita v e de suas derivadas vale 1, de forma que a equação diferencial é linear. c) � Atenção! Para visualização completa da equação utilize a rolagem horizontal y′ − 3x(y′′) 3 − x4 = 3 sen (x) y − 3y′ + 4x2 = 0 − xyz = 2 ∂3m ∂z3 yy′ − lnu = y′′ y′ − 3x(y′′) 3 − x4 = 3 sen (u) v − 3v′ + 4u2 = 0 8z − xyz = 2 ∂3m ∂z3 A equação diferencial é uma EDP. Os coeficientes dependem apenas das variáveis independentes x,y e z. O expoente da derivada parcial da variável m, que existe na equação, é a unidade. Com isso, a EDP é linear. d) � Atenção! Para visualização completa da equação utilize a rolagem horizontal A equação diferencial é uma EDO. Apesar daincógnita y e de suas derivadas estarem elevadas ao expoente 1, existe um produto entre y e a derivada de y. Logo, a equação é não linear. EXPRESSÃO DAS EDO LINEARES As equações diferenciais ordinárias lineares de ordem n podem e sempre serão expressas da forma � Atenção! Para visualização completa da equação utilize a rolagem horizontal Nesse caso, a função incógnita s depende da variável independente e , com 0 ≤ j ≤ n, são os coeficientes da equação. Se for impossível transformar uma EDO na forma apresentada, ela não será uma equação diferencial linear. Isso também é verdade para as equações diferenciais parciais. ELAS SERÃO LINEARES SE FOR POSSÍVEL COLOCAR DA FORMA SEMELHANTE A UM POLINÔMIO. Vejamos o exemplo, a seguir, para uma EDP de segunda ordem, em que a variável z depende das variáveis x e y: � Atenção! Para visualização completa da equação utilize a rolagem horizontal HOMOGENEIDADE E COEFICIENTES Uma equação linear pode ser classificada como: HOMOGÊNEA yy′ − lnu = y′′ fn (u) s (n) + fn−1 (u) s (n−1) + fn−2 (u) s (n−2) + … + f1 (u) s ′ + fo (u) s = g(u) u fj (u) f (x, y) zxx + g (x, y) zxy + h (x, y) zyy + p (x, y) zx + q (x, y) zy + m (x, y) z = w(x, y) NÃO HOMOGÊNEA AS EQUAÇÕES HOMOGÊNEAS SÃO AQUELAS QUE NÃO APRESENTAM TERMOS SEM A VARIÁVEL DEPENDENTE OU SEM UMA DE SUAS DERIVADAS. NOS EXEMPLOS ANTERIORES, AS EQUAÇÕES HOMOGÊNEAS APRESENTAM OU IGUAIS A ZERO. CASO CONTRÁRIO, A EQUAÇÃO SERÁ NÃO HOMOGÊNEA, E O TERMO OU SERÁ DENOMINADO TERMO NÃO HOMOGÊNEO ATENÇÃO SE OS COEFICIENTES DA EQUAÇÃO DIFERENCIAL QUE MULTIPLICAM A VARIÁVEL DEPENDENTE OU SUAS DERIVADAS FOREM TODOS NÚMEROS REAIS, DIZEMOS QUE A EQUAÇÃO TEM COEFICIENTES CONSTANTES. A EQUAÇÃO DIFERENCIAL SERÁ DE COEFICIENTES VARIÁVEIS CASO ALGUM DESSES TERMOS DEPENDA DA(S) VARIÁVEL(IS) INDEPENDENTE(S). Classifique as equações diferenciais lineares, a seguir, quanto à homogeneidade e indique se são ou não de coeficientes constantes. A) � Atenção! Para visualização completa da equação utilize a rolagem horizontal B) � Atenção! Para visualização completa da equação utilize a rolagem horizontal C) g(u) w(x, y) g(u) w(x, y) s′ − 2s + 5x4 = cosx − 4y = d3y dz d2y dz2 � Atenção! Para visualização completa da equação utilize a rolagem horizontal D) � Atenção! Para visualização completa da equação utilize a rolagem horizontal RESOLUÇÃO � Atenção! Para visualização completa da equação utilize a rolagem horizontal Observe que os coeficientes que multiplicam a variável dependente e suas derivadas são numéricos, de modo que a equação diferencial é de coeficientes constantes. O termo independente vale cos x, logo, é diferente de zero, e a equação não homogênea. Podemos reescrever a equação na forma: � Atenção! Para visualização completa da equação utilize a rolagem horizontal B) � Atenção! Para visualização completa da equação utilize a rolagem horizontal A EDO linear tem todos os coeficientes numéricos e não tem nenhum termo que independe de y e de suas derivadas. Portanto, trata-se de uma equação diferencial de coeficientes constantes e homogênea. Podemos reescrever a equação na forma: � Atenção! Para visualização completa da equação utilize a rolagem horizontal C) 2z + m = 2 ∂3m ∂z3 sen (x) y − 3y′ + 4x2y = 0 a)s′ − 2s + 5x4 = cosx s′ − 2s = cosx − 5x4 − 4y = d3y dz d2y dz2 y′′′ − y′′ − 4y = 0 2z + sen w − = 2m ∂3m ∂z3 ∂m ∂w � Atenção! Para visualização completa da equação utilize a rolagem horizontal A equação diferencial linear é uma EDP com coeficientes que dependem das variáveis independentes z e w. Além disso, apresenta um termo que independe da variável dependente m, de modo que é uma equação de coeficientes variáveis e não homogênea. Reescrevendo a equação na forma polinomial, temos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal D) � Atenção! Para visualização completa da equação utilize a rolagem horizontal A EDO linear apresenta coeficientes que são funções da variável independente x, de forma que constitui uma equação de coeficientes variáveis. Não existe termo que independa da variável y e de suas derivadas, logo, é uma equação homogênea. Ordenando a equação, ficamos com: � Atenção! Para visualização completa da equação utilize a rolagem horizontal As equações diferenciais lineares, assim como as equações algébricas, são mais fáceis de serem resolvidas do que as não lineares. No próximo módulo, analisaremos a resolução de alguns tipos de equações diferenciais de primeira ordem. TEORIA NA PRÁTICA Sabemos que não existe um método único para solucionar qualquer equação diferencial. Existem métodos que são apropriados para grupos de equações diferenciais. Desse modo, é importante classificar uma equação para definir o método de solução que pode ser utilizado na prática. Escolha o método para solucionar a equação apresentada a seguir: Classifique quanto à ordem e ao grau, indique se são lineares ou não lineares, homogênea ou não homogênea e, por fim, se têm coeficientes constantes ou variáveis. � Atenção! Para visualização completa da equação utilize a rolagem horizontal RESOLUÇÃO CLASSIFICAÇÃO DA EQUAÇÃO DIFERENCIAL 2z − − 2m = −sen w ∂3m ∂z3 ∂m ∂w sen (x) y − 3y′ + 4x2y = 0 3y′ − (sen(x) + 4x2)y = 0 y′′ + 2x y′ − 3, 4 y = x3 + ex MÃO NA MASSA 1. MARQUE A ALTERNATIVA QUE APRESENTA UMA EQUAÇÃO DIFERENCIAL DE TERCEIRA ORDEM: A) � Atenção! Para visualização completa da equação utilize a rolagem horizontal B) � Atenção! Para visualização completa da equação utilize a rolagem horizontal C) � Atenção! Para visualização completa da equação utilize a rolagem horizontal D) � Atenção! Para visualização completa da equação utilize a rolagem horizontal E) � Atenção! Para visualização completa da equação utilize a rolagem horizontal 2. MARQUE A ALTERNATIVA QUE APRESENTA UMA EQUAÇÃO DIFERENCIAL DE GRAU DOIS: A) � Atenção! Para visualização completa da equação utilize a rolagem horizontal B) = (3 − y) dy dx 4 x + ( ) 2 = cos(x + y) ∂2v ∂x2 ∂2v ∂y2 s′′ + cos (x)s′ = 4x ( ) 2 + 4 = x ∂3y ∂u∂v∂w ∂y ∂w = −2x d4x dt4 3 = (3 − y2) dy dx 4 x � Atenção! Para visualização completa da equação utilize a rolagem horizontal C) � Atenção! Para visualização completa da equação utilize a rolagem horizontal D) � Atenção! Para visualização completa da equação utilize a rolagem horizontal E) � Atenção! Para visualização completa da equação utilize a rolagem horizontal 3. MARQUE A ALTERNATIVA QUE APRESENTA UMA EQUAÇÃO DIFERENCIAL LINEAR HOMOGÊNEA A) � Atenção! Para visualização completa da equação utilize a rolagem horizontal B) � Atenção! Para visualização completa da equação utilize a rolagem horizontal C) � Atenção! Para visualização completa da equação utilize a rolagem horizontal D) � Atenção! Para visualização completa da equação utilize a rolagem horizontal E) � Atenção! Para visualização completa da equação utilize a rolagem horizontal 4. QUAL É A ORDEM E O GRAU DA EQUAÇÃO DIFERENCIAL ? + ( ) 2 − ln (xy) = 2 ∂3v ∂x3 ∂2v ∂y2 (s′′) 2 + cos (x)(s′) 3 = 4x ( ) 3 + 4 = x ∂3y ∂u∂v∂w ∂y ∂w = −2x d4x dt4 + ( ) 2 − ln (xy) = 2 ∂3v ∂x3 ∂2v ∂y2 s′′ + cos (x) s′ = 4x 2xy′– cos(x)y′′ + 4ey = 6x + 3v = 5u du dv d2u dv2 uv − v′ + v′′ (2v′ + 1) = 0 2(u′′′)2 − 2v5u + (u′′)3 = 4 ATENÇÃO! PARA VISUALIZAÇÃO COMPLETA DA EQUAÇÃO UTILIZE A ROLAGEM HORIZONTAL A) Ordem 2 e grau 3 B) Ordem 3 e grau 5 C) Ordem 1 e grau 1 D) Ordem 1 e grau 2 E) Ordem 3 e grau 2 5. SEJA A EQUAÇÃO DIFERENCIAL , MARQUE A ALTERNATIVA QUE APRESENTA UMA CLASSIFICAÇÃO VERDADEIRA EM RELAÇÃO A ESSA EQUAÇÃO. A) Linear, coeficientes constantes e homogênea B) Não linear, coeficientes variáveis e homogênea C) Linear, coeficientes variáveis e não homogênea D) Não linear, coeficientes constantes e homogênea E) Linear, coeficientes variáveis e homogênea 6. SEJA A EQUAÇÃO DIFERENCIAL F(X,Y) Y’’ – 2Y’ + Y = X G(X,Y). MARQUE A ALTERNATIVA QUE APRESENTA VALORESPARA F(X,Y) E G(X,Y) DE FORMA QUE A EQUAÇÃO DIFERENCIAL SEJA DE SEGUNDA ORDEM, LINEAR E HOMOGÊNEA. A) f(x,y) = 0 e g(x,y) = y B) f(x,y) = x e g(x,y) = y C) f(x,y) = y e g(x,y) = 0 D) f(x,y) = x e g(x,y) = x E) f(x,y) = y e g(x,y) = x GABARITO 1. Marque a alternativa que apresenta uma equação diferencial de terceira ordem: A alternativa "D " está correta. Como vimos, uma equação diferencial terá ordem 3 se a derivada de maior ordem for de terceira ordem. A alternativa A apresenta uma equação diferencial de primeira ordem; as alternativas B e C apresentam equações diferenciais de segunda ordem; a alternativa E apresenta uma equação diferencial de quarta ordem. Portanto, a resposta é a alternativa D, que apresenta uma derivada parcial de ordem 3 como derivada de maior ordem. 2. Marque a alternativa que apresenta uma equação diferencial de grau dois: A alternativa "C " está correta. Como vimos, uma equação diferencial terá grau dois se a derivada de maior ordem estiver elevada a um expoente 2. Na equação da letra A, a ordem é um, e a derivada de ordem 1 está elevada ao expoente unitário, sendo grau também um. A ordem da equação diferencial da letra B 2y′ − cos (x)y′′ + 3y = 5 vale três, pois é a derivada de mais alta ordem. O grau dessa equação também é um. A derivada de mais alta ordem na alternativa da letra C é de segunda ordem. A derivada está elevada a um expoente dois. Desse modo, é uma equação diferencial de ordem 2 e grau 2, sendo a resposta correta. Na letra D, existe uma equação de ordem três e grau três. Por fim, na letra E, a equação é de ordem quatro e grau um. 3. Marque a alternativa que apresenta uma equação diferencial linear homogênea A alternativa "D " está correta. A equação, para ser linear, deve conter coeficientes que independam da variável independente. Assim, a alternativa E não é linear. Para ser linear, tanto a incógnita quanto as derivadas devem ter expoente elevado a um. Com isso, as alternativas A e C não apresentam equações lineares. Já as alternativas B e D são equações lineares, mas como a alternativa B apresenta um termo que independe da incógnita, a equação é não homogênea. A equação diferencial da letra D — resposta da questão — não apresenta um termo que independe da função incógnita, de modo que é a uma equação linear homogênea. 4. Qual é a ordem e o grau da equação diferencial ? Atenção! Para visualização completa da equação utilize a rolagem horizontal A alternativa "E " está correta. A ordem da equação diferencial é a ordem da derivada de maior ordem na equação. Analisando a equação, vemos que existe uma derivada de terceira ordem e uma de segunda ordem, de forma que a ordem vale 3. O grau é o expoente desta derivada, que está elevada ao expoente dois na equação, de modo que o grau da equação diferencial vale 2. Então, a letra E é a alternativa correta. 5. Seja a equação diferencial , marque a alternativa que apresenta uma classificação verdadeira em relação a essa equação. A alternativa "C " está correta. A variável dependente é y, e a independente é a variável x. Como todos os coeficientes independem de y e todas as incógnitas e as derivadas têm expoente um, trata-se de uma equação diferencial linear. Como existe um coeficiente que é uma função de x, trata-se de uma equação de coeficientes variáveis. Por fim, existe um termo que independe de y, sendo uma equação não homogênea. Logo, a alternativa que apresenta as características linear, coeficientes variáveis e não homogênea é a alternativa C. 6. Seja a equação diferencial f(x,y) y’’ – 2y’ + y = x g(x,y). Marque a alternativa que apresenta valores para f(x,y) e g(x,y) de forma que a equação diferencial seja de segunda ordem, linear e homogênea. A alternativa "B " está correta. Veja a solução no vídeo abaixo CLASSIFICAÇÃO DE EQUAÇÕES DIFERENCIAIS GABARITO 2(u′′′)2 − 2v5u + (u′′)3 = 4 2y′ − cos (x)y′′ + 3y = 5 VERIFICANDO O APRENDIZADO 1. MARQUE A ALTERNATIVA QUE APRESENTA UMA EQUAÇÃO DE ORDEM TRÊS E GRAU DOIS. A) � Atenção! Para visualização completa da equação utilize a rolagem horizontal B) � Atenção! Para visualização completa da equação utilize a rolagem horizontal C) � Atenção! Para visualização completa da equação utilize a rolagem horizontal D) � Atenção! Para visualização completa da equação utilize a rolagem horizontal E) � Atenção! Para visualização completa da equação utilize a rolagem horizontal 2. SEJA A EQUAÇÃO DIFERENCIAL ATENÇÃO! PARA VISUALIZAÇÃO COMPLETA DA EQUAÇÃO UTILIZE A ROLAGEM HORIZONTAL MARQUE A ALTERNATIVA QUE APRESENTA UMA CLASSIFICAÇÃO VERDADEIRA EM RELAÇÃO A ESSA EQUAÇÃO. A) Ordinária, linear e não homogênea B) Ordinária, de grau dois e ordem um C) Parcial, linear e homogênea D) Parcial, de coeficientes constantes e de segunda ordem E) Ordinária, linear de grau dois GABARITO x′′′ − yx = 4 − 3v2 = m + 1 ∂2m ∂v2 + (cos(u) ) 2 = 4v dv du d3v du3 ( ) 3 = 3uv ∂2y ∂u∂v (s′′) 2 + cos (x)(s′) 3 = 4x euv − ( ) 2 + ln(u) = dv du d2v du2 1. Marque a alternativa que apresenta uma equação de ordem três e grau dois. A alternativa "C " está correta. Na letra A, há uma equação diferencial de terceira ordem, pois é a ordem da mais alta derivada da variável dependente x. Essa derivada está elevada à unidade, de forma que seu grau é um. A equação da letra B é uma equação diferencial de segunda ordem e, como a derivada de mais alta ordem está elevada ao expoente um, seu grau também vale um. A letra C, resposta da questão, apresenta a derivada de mais alta ordem de terceira ordem e elevada ao expoente dois. Assim é uma equação de ordem 3 e grau 2. Na letra D, a derivada de mais alta ordem é de segunda ordem e está elevada ao expoente três. Com isso, temos uma equação de ordem dois e grau três. Por fim, uma equação de ordem dois com grau 2, pois há uma derivada de segunda ordem elevada ao expoente dois. A resposta correta, então, é a letra C. 2. Seja a equação diferencial Atenção! Para visualização completa da equação utilize a rolagem horizontal Marque a alternativa que apresenta uma classificação verdadeira em relação a essa equação. A alternativa "A " está correta. Comentário: Analisando a equação, verificamos que a variável dependente é v, e a independente, a variável u. A equação tem derivadas de v em relação a um, sendo uma equação diferencial ordinária. A derivada de mais alta ordem é de segunda ordem e está elevada a um expoente unitário, de modo que se trata de uma equação de ordem dois e grau um. A equação é linear, já que os coeficientes dependem apenas de u, e não há termo que independa de v, logo, é uma equação homogênea. Por fim, como os coeficientes são funções de u, é uma equação de coeficientes variáveis. Com isso, a única alternativa verdadeira é a letra A. MÓDULO 3 Calcular equações diferenciais de primeira ordem MÉTODOS DE RESOLUÇÃO DE EQUAÇÃO DIFERENCIAL DE PRIMEIRA ORDEM euv − ( ) 2 + ln(u) = dv du d2v du2 SOLUÇÃO DAS EQUAÇÕES DIFERENCIAIS DE PRIMEIRA ORDEM Nos módulos anteriores, estudamos os conceitos iniciais das equações diferenciais e analisamos como verificar se uma função é ou não solução de uma equação diferencial. Vimos que não existe um método geral para resolver todos os tipos de equação diferencial. Existem, no entanto, alguns métodos com grande abrangência que resolvem alguns tipos de equações diferenciais. Neste módulo, estudaremos alguns dos métodos para resolução de equações diferenciais ordinárias de primeira ordem. Vale lembrar que a solução de uma equação diferencial ordinária será uma função real que satisfaz a equação dada. PROBLEMA DE VALOR INICIAL CHAMAREMOS DE SOLUÇÃO GERAL A FAMÍLIA DE FUNÇÕES QUE SATISFAZ A EQUAÇÃO, E DE SOLUÇÃO PARTICULAR UMA FUNÇÃO QUE FAZ PARTE DESSA SOLUÇÃO GERAL E ATENDE ÀS CONDIÇÕES INICIAIS FORNECIDAS. ESSE TIPO DE PROBLEMA SERÁ DENOMINADO PROBLEMA DE VALOR INICIAL. Veremos, aqui, a solução de equações diferenciais de primeira ordem, portanto a derivada de maior ordem que apareceráserá a derivada de primeira ordem, elevada ao expoente 1. A EDO de primeira ordem será do tipo: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Mesmo para o caso das EDOs de primeira ordem, não existe um método único que abrange uma solução para todas as EDOs. Analisaremos alguns métodos de resolução, e todos os que serão apresentados se aplicam para equações diferenciais de primeira ordem que podem ser colocadas na forma: � Atenção! Para visualização completa da equação utilize a rolagem horizontal A função f(x,y) será uma função que pode depender da variável independente y e/ou da variável dependente x. MÉTODO DIRETO O primeiro método é denominado de método direto e sua solução é mais simples. Ela ocorrerá quando for possível isolar a derivada de y em um lado e, do outro, uma função que dependa apenas da variável independente x. Para a solução direta, basta usar, depois, a integração direta para obter a solução. Em outras palavras, quando . y′ + a (x, y) y = b(x, y) y′ = f(x, y) f(x, y) = f(x) � Atenção! Para visualização completa da equação utilize a rolagem horizontal Vejamos o exemplo: OBTENHA A SOLUÇÃO GERAL DA EQUAÇÃO DIFERENCIAL RESOLUÇÃO Para obter o valor da função incógnita y, devemos tentar isolar e obter o valor da derivada da variável � Atenção! Para visualização completa da equação utilize a rolagem horizontal Então: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Basta, agora, integrar ambos os lados, de forma que a solução geral da equação diferencial será y′ = f (x) → dy = f (x) dx → y = ∫ f (x) dx 4x3 − y′ − 2x2 + 1 = 0. 4x3 − y′ − 2x2 + 1 = 0 → y′ = 4x3 − 2x2 + 1 = 4x3 − 2x2 + 1 dy dx dy = (4x3 − 2x2 + 1) dx ∫ dy = ∫ (4x3 − 2x2 + 1)dx y = x4 − x3 + x + k 2 3 , com k número real. OBTENHA UMA SOLUÇÃO PARTICULAR DA EQUAÇÃO DIFERENCIAL , TAL QUE, PARA ATENDER A CONDIÇÃO X = 1, O VALOR DE Y VALE 2. RESOLUÇÃO Obtivemos a solução geral , com k número real. Agora, iremos substituir a condição inicial e obter a solução particular: Para x = 1 → No entanto, y = 2 para x = 1, de forma que A solução particular será SOLUÇÃO PARTICULAR Infelizmente, nem sempre será possível realizar a solução da equação geral por meio de uma manipulação matemática, como apresentado no exemplo. Por isso, torna-se necessário definir métodos que podem ser empregados na obtenção da solução de um conjunto específico de equações. Analisaremos três métodos: EQUAÇÕES SEPARÁVEIS EQUAÇÃO EXATA EQUAÇÕES LINEARES 4x3 − y′ − 2x2 + 1 = 0 y = x4 − x3 + x + k 2 3 y = 14 − 13 + 1 + k = + k 2 3 4 3 + k = 2 → k = 4 3 2 3 = x4 − x3 + x + 2 3 2 3 ANTES DE ANALISARMOS TAIS MÉTODOS, PRECISAMOS NOS PERGUNTAR: COMO GARANTIR QUE UMA EDO DE PRIMEIRA ORDEM TENHA SEMPRE UMA SOLUÇÃO PARTICULAR PARA O PROBLEMA E, QUANDO ESSA SOLUÇÃO EXISTIR, QUE ELA SEJA ÚNICA? RESPOSTA Vejamos o teorema: Seja uma EDO do tipo , em que f(x,y) é contínua em um retângulo R contendo o ponto (x0,y0), com y(x0) = y0. Essa EDO terá sempre uma solução no retângulo R que contém (x0,y0) e essa solução será única, sempre que for contínua em R. Retornando ao nosso exemplo, para EDO , repare que é contínua para todo (x,y) em . Além disso, é contínua em todo . Portanto, essa EDO terá sempre uma solução para todos os problemas de valor inicial, e a solução será única em cada caso. Nos tópicos seguintes, estaremos focados em determinar a solução de uma EDO de primeira ordem para alguns tipos de EDOs de primeira ordem. EDO SEPARÁVEIS DE PRIMEIRA ORDEM Estamos estudando a EDO de primeira ordem: � Atenção! Para visualização completa da equação utilize a rolagem horizontal y′ = f(x, y) ∂f ∂y y′ = 4x3 − 2x2 + 1 f(x, y) = 4x3 − 2x2 + 1 R2 = 0 ∂f ∂y R2 y′ = f(x, y) A equação será do tipo separável quando for possível transformar a função f(x,y) em um produto h(x)g(y), isto é, em um produto de uma função, que só depende da variável independente, por outra que só depende da variável dependentes. Resolvendo, temos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal A equação é possível para um intervalo onde g(y) é diferente de zero. Desse modo, para solucionar a equação vamos aplicar a integral definida em ambos os lados e obter a solução geral: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Apesar de cada integral definida apresentar uma constante de integração, as duas constantes podem ser substituídas apenas por uma constante de integração k. RELEMBRANDO Atenção às regras de integração, já que elas serão bastante usadas na solução de equações diferenciais. A solução geral pode ser uma função implícita relacionado y com o x, não sendo possível, por vezes, explicitar a função. Vejamos alguns exemplos desse método a seguir: UMA POPULAÇÃO DE BACTÉRIA P, MEDIDA EM MILHÕES, CRESCE EM RELAÇÃO AO TEMPO T, MEDIDO EM HORAS, POR MEIO DO SEGUINTE MODELO y′ = h (x) g (y) y′ = h (x) g (y) → = h (x) g (y) → = h (x) dx dy dx dy g (y) ∫ dy = ∫ h (x) dx + k, com k real 1 g(y) = 2(p − 1) dp dt � Atenção! Para visualização completa da equação utilize a rolagem horizontal Determine a quantidade de bactérias que haverá após 5 horas, sabendo que, para t = 0, o valor de p é de 2 milhão de bactérias. RESOLUÇÃO Precisamos, inicialmente, solucionar a equação diferencial de primeira ordem para descobrir como p varia com o tempo. Repare que temos uma EDO separável, uma vez que: , para � Atenção! Para visualização completa da equação utilize a rolagem horizontal Assim: , com � Atenção! Para visualização completa da equação utilize a rolagem horizontal Como p é maior ou igual a 2, então, p – 1 é positivo, de modo que � Atenção! Para visualização completa da equação utilize a rolagem horizontal Com isso, a solução geral da EDO será , com k real, mas para t = 0, p = 2. � Atenção! Para visualização completa da equação utilize a rolagem horizontal Já a solução particular será . Para obter a população para t = 5 horas, temos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Observe que, no exemplo, a expressão obtida no método levava em conta que . Como, no exemplo, o valor de p ≥ 2, não se perdeu nenhuma solução possível. Pode acontecer em alguns casos, no entanto, que g(y) = 0 forneça solução para equação diferencial analisada. Nesse caso, diz-se que as soluções são singulares, já que não podem ser obtidas pelo método empregado. Vejamos o exemplo: Obtenha a solução da equação = 2 (p − 1) → dp = 2 dt dp dt 1 p − 1 p ≠ 1 ∫ dp = ∫ 2 dt → ln |p − 1| = 2t + k1 p − 1 k ∈ R |p − 1| = p − 1 ln |p − 1| = ln (p − 1) =2t + k → p − 1 = exp(2t + k) p = 1 + exp(2t + k) p = 1 + exp (2.0 + k) = 2 exp (k) = 2 − 1 = 1 → k = 0 p = 1 + exp(2t) p (5) = 1 + exp (2.5) = 1 + exp (10) ≈ 22.028 milhões p ≠ 1 RESOLUÇÃO Trata-se de uma EDO de primeira ordem. A EDO é do tipo separável, uma vez que: , com � Atenção! Para visualização completa da equação utilize a rolagem horizontal Dessa maneira: , com real. , com k real e para . � Atenção! Para visualização completa da equação utilize a rolagem horizontal Vemos, com isso, que essa solução geral não engloba o caso de y = 0. Se fizermos y = 0 na EDO, observaremos que é solução para a EDO, pois é atendida para a função y = 0. A função y = 0 seria uma solução singular para a EDO. Assim, a solução da EDO seria , com k real , ou y = 0. ATENÇÃO Às vezes, com uma substituição de variável, podemos transformar uma EDO que não é do tipo separável para uma EDO separável. Vejamos um exemplo desse caso: DETERMINE A FUNÇÃO U, QUE É SOLUÇÃO DA EQUAÇÃO ATENÇÃO! PARA VISUALIZAÇÃO COMPLETA DA EQUAÇÃO UTILIZE A ROLAGEM HORIZONTAL E QUE PARA V = 0 TENHA U = 1. y′ = x3y2 y′ = x3y2 → = x3y2 → dy = x3dx dy dx 1 y2 ≠ 0 ∫ dy = ∫ x3dx + k → − = x4 + k1 y2 1 y 1 4 y = − = − 1 x4 + k1 4 4 x4 + 4k y ≠ 0 y′ = x3y2 y = − = − 1 x4 + k1 4 4 x4 +4k u′ = 2v + 2u + v2 + u2 + 4v + 4u + 2uv + 8 RESOLUÇÃO Temos uma EDO de primeira ordem que relaciona a incógnita u com a variável independente v. Inicialmente, não é possível se transformar a EDO para uma EDO separável, isto é: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Manipulando a equação, no entanto, teremos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Se fizemos , � Atenção! Para visualização completa da equação utilize a rolagem horizontal Derivando w em relação à variável independente v, teremos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Substituindo Então , que é uma EDO separável. Com isso: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Retornando a variável original, podemos obter uma equação implícita, que relaciona u com v da seguinte forma: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Para v = 0, temos u = 1 � Atenção! Para visualização completa da equação utilize a rolagem horizontal Assim sendo, a função u que atende às condições do problema, em relação a variável v, é obtida pela equação: u′ = f (u) g(v) u′ = 2v + 2u + v2 + u2 + 4v + 4u + 2uv + 8 u′ = (2v + 2u + 4) + (v2 + u2 + 4 + 4v + 4u + 2uv) u′ = 2 (v + u + 2) + (v + u + 2)2 w = v + u + 2 u′ = 2 (v + u + 2) + (v + u + 2) 2 → u′ = 2w + w2 w = v + u + 2 → w′ = 1 + u′ w′ − 1 = u′ = 2w + w2 w′ = 1 + 2w + w2 w′ = = 1 + 2w + w2 → dw = dv dw dv 1 1 + 2w + w2 ∫ dw = ∫ dv + k1 1 + 2w + w2 ∫ dw = ∫ dv + k1 (1 + w)2 − = v + k 1 1 + w − = v + k → v + = − k 1 1 + (v + u + 2) 1 v + u + 3 0 + = − k → k = − 1 0 + 1 + 3 1 4 � Atenção! Para visualização completa da equação utilize a rolagem horizontal EDO EXATAS DE PRIMEIRA ORDEM Uma equação diferencial de primeira ordem será considerada exata em um domínio S se puder ser colocada na seguinte forma: e existir uma função M(x,y), definida em S tal que: E para todos os pontos (x,y) da região S. � Atenção! Para visualização completa da equação utilize a rolagem horizontal Observe que u(x,y) e v(x,y) são duas funções escalares de x e y, e serão dadas pela equação diferencial. Uma forma de verificar se a equação é exata é analisando: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Vejamos: VERIFICAÇÃO DE EQUAÇÃO EXATA OBSERVE SE A EQUAÇÃO DIFERENCIAL v + = 1 v + u + 3 1 4 u (x, y) + v (x, y) y′ = 0 (x, y) = u(x, y) ∂M ∂x (x, y) = v(x, y) ∂M ∂y = = ∂u ∂y ∂v ∂x ∂2M ∂x∂y 2 (x2 + y) y′ = −4xy − x2 É UMA EQUAÇÃO DIFERENCIAL EXATA. RESOLUÇÃO Colocando a equação na forma � Atenção! Para visualização completa da equação utilize a rolagem horizontal Portanto: e � Atenção! Para visualização completa da equação utilize a rolagem horizontal Obtendo e Retornando a solução da equação de primeira ordem exata, vamos definir a sua solução. Se a equação diferencial for exata, então a solução da equação diferencial será dada pela equação implícita , com k real. Sendo a função a função tal que � Atenção! Para visualização completa da equação utilize a rolagem horizontal Vamos provar essa afirmativa: Lembrando de o estudo da função escalar, temos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Substituindo, teremos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Mas u (x, y) + v (x, y) y′ = 0 2 (x2 + y) y′ = −4xy − x2 → 2 (x2 + y) y′ + (4xy + x2) = 0 u (x, y) = 4xy + x2 v (x, y) = 2 (x2 + y) (x, y) = 4x ∂u ∂y (x, y) = 4x ∂v ∂x u(x, y) + v(x, y)y′ = 0 M(x, y) = k M(x, y) (x, yt) = u(x, y)e (x, y) = v(x, y) ∂M ∂x ∂M ∂y M(x, y) = k → = 0 dM dx = + = 0 dM dx ∂M ∂x dx dx ∂M ∂y dy dx + = 0 → + y′ = 0 ∂M ∂x dx dx ∂M ∂y dy dx ∂M ∂x ∂M ∂y e , então: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Dessa forma, satisfazemos a equação diferencial apresentada. Precisamos obter a função M(x,y) para solucionar a equação diferencial exata. A função M(x,y) pode ser obtida por meio da integração. Isso significa que, se desejarmos obter a função M(x,y), basta resolvermos: E Logo, obtemos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Lembre-se de que, como estamos integrando em relação a x, e M(x,y) é uma função escalar que depende de x e y, a constante de integração será uma função de y, representada por g(y). Derivando, parcialmente, M(x,y) em relação a y, teremos: (x, y) = u(x, y) ∂M ∂x (x, y) = v(x, y) ∂M ∂y + y′ = 0 → u (x, y) + v (x, y) y′ = 0 ∂M ∂x ∂M ∂y (x, y) = u(x, y) ∂M ∂x (x, y) = v(x, y) ∂M ∂y (x, y) = u (x, y) → M = ∫ u (x, y) dx + g (y) ∂M ∂x (x, y) = (∫ u (x, y) dx + g (y)) = (∫ u (x, y) dx) + g′(y)∂M ∂y ∂ ∂y ∂ ∂y No entanto: De modo que: Substituindo, obtemos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal O equacionamento parece complicado, mas basta seguir os passos para obter a solução. Repare que a solução da equação será, normalmente, uma equação implícita. Veja: EXEMPLO DETERMINE A SOLUÇÃO DA EQUAÇÃO DIFERENCIAL RESOLUÇÃO No exemplo anterior, provamos que se trata de uma equação diferencial exata de primeira ordem, com: e � Atenção! Para visualização completa da equação utilize a rolagem horizontal Assim: g′ (y) = (x, y) − (∫ u (x, y) dx)∂M ∂y ∂ ∂y (x, y) = v (x, y) → g′ (y) = v (x, y) − (∫ u (x, y) dx)∂M ∂y ∂ ∂y g (y) = ∫ [v (x, y) − (∫ u (x, y) dx)] dy ∂ ∂y M(x, y) = ∫ u (x, y) dx + ∫ [v (x, y) − (∫ u (x, y) dx)] dy∂ ∂y 2 (x2 + y) y′ = −4xy − x2 u (x, y) = 4xy + x2 v (x, y) = 2 (x2 + y) � Atenção! Para visualização completa da equação utilize a rolagem horizontal Vamos, agora, obter a derivada parcial de M em relação a y: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Como: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Então: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Com isso, substituindo: � Atenção! Para visualização completa da equação utilize a rolagem horizontal A solução da equação diferencial será dada pela equação implícita , k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal O valor de k será obtido por uma condição inicial. Por exemplo, considere que y(0) = 2, de modo que: � Atenção! Para visualização completa da equação utilize a rolagem horizontal A solução particular será . EDO LINEAR DE PRIMEIRA ORDEM O terceiro método que estudaremos neste módulo é apropriado para as equações diferenciais lineares. A equação ordinária linear de primeira ordem pode ser sempre expressa na forma: (x, y) = u (x, y) = 4xy + x2 → M = ∫ (4xy + x2) dx + g(y) = 2x2y + x 3 + g(y) ∂M ∂x 1 3 (x, y) = (2x2y + x 3 + g(y)) = 2x2 + g′(y)∂M ∂y ∂ ∂y 1 3 (x, y) = v (x, y) = 2 (x2 + y) ∂M ∂y 2x2 + g′ (y) = 2x2 + g′ (y) → g′ (y) = 2y → g (y) = ∫ 2ydy = y2 M = 2x2y + x 3 + g (y) = 2x2y + x 3 y2 + y2 1 3 1 3 2x2y + x 3 y2 + y2 = k 1 3 2.02.2 + .0 3 22 + 22 = k → k = 4 1 3 2x2y + x 3 y2 + y2 = 4 1 3 y′ + u (x) y = v(x) � Atenção! Para visualização completa da equação utilize a rolagem horizontal Na qual u(x) e v(x) são duas funções da variável independente x, dadas na equação diferencial. Vamos considerar uma função auxiliar dada por P(x) y. Derivando essa função por meio da regra do produto, teremos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Vamos partir da equação diferencial linear e multiplicar os dois lados pela função P(x): � Atenção! Para visualização completa da equação utilize a rolagem horizontal Compare as duas equações. Se fizemos , igualaremos o lado esquerdo da segunda equação com a primeira. Isto é, se , então . Integrando ambos os lados, teremos: [P (x) y]′ = P ′ (x) y + P (x) y′ p(x)y′ + P(x)u (x) y = P(x)v(x) P′ (x) = P(x)u (x) P ′ (x) = P(x)u (x) [P (x) y] ′ = P(x)v(x) ∫ [P (x) y]′dx = ∫ P (x) v (x) dx + k , k real P (x) y = ∫ P (x) v (x) dx + k , k real y = (∫ P (x) v (x) dx + k ) , k real 1 P(x) � Atenção! Para visualização completa da equação utilize a rolagem horizontal Conseguiremos, então, obter a solução da equação diferencial linear de primeira ordem se escolhermos corretamente a função P(x). A função P(x) é denominada fator integrante. CÁLCULO DE P(X) Vamos, agora, ver como calculamos P(X). Retornando, , então: Assim: logo: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Como necessitamos apenas de um P(x) que atenda à condição, não precisamos incluir a constante de integração. Portanto, para obter o fator integrante, basta resolver: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Após obter P(x), devemos obter a equação geral da EDO linear de primeira ordem: P ′ (x) = P(x)u (x) = u(x) → ∫ = ∫ u (x) dx P ′(x) P(x) dP P ln |P(x)| = ∫ u (x) dx + C, C real P (x) = exp(∫ u (x) dx + C) , C real P (x) = exp(∫ u (x) dx) , K REAL � Atenção! Para visualização completa da equação utilize a rolagem horizontal ETAPAS Agora que já entendemos o procedimento, vamos listar o passo a passo para não esquecermos nenhuma etapa do processo: ETAPA A Colocar a equação diferencial linear na sua forma padrão, isto é, o coeficiente que multiplica y’ deve ser um. ETAPA B Obter o fator integrante, usando o valor de u(x) da equação, pela fórmula: � Atenção! Para visualização completa da equação utilize a rolagem horizontal ETAPA C c) Obter o valor da integral a seguir, usando v(x) da equação diferencial dada: � Atenção! Para visualização completa da equação utilize a rolagem horizontal ETAPA D Determinar a solução geral da equação diferencial pela fórmula: � Atenção! Para visualização completa da equação utilize a rolagem horizontal A SOLUÇÃO PARTICULAR SERÁ OBTIDA ACHANDO-SE O VALOR DE K REAL POR MEIO DA CONDIÇÃO INICIAL. Vejamos um exemplo a seguir: Vamos obter a equação geral da equação diferencial . RESOLUÇÃO Vamos colocar a equação diferencial linear na forma padrão: y = (∫ P (x) v (x) dx + k ) 1 P(x) P (x) = exp(∫ u (x) dx) ∫ P (x) v (x) dx y = (∫ P (x) v (x) dx + k ) , k real1 P(x) 2y′ − 2y − 4ex = 0 javascript:void(0) javascript:void(0) javascript:void(0) javascript:void(0) � Atenção! Para visualização completa da equação utilize a rolagem horizontal Então, . Agora, devemos obter o fator integrante: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Resolvendo a integral, teremos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal De forma que: , k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal Com isso, a solução geral será: , k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal Vamos obter a solução que atenda à equação diferencial e que, para x = 0, tenhamos y = 1. RESOLUÇÃO No exemplo anterior, obtivemos a equação geral , k real, fazendo � Atenção! Para visualização completa da equação utilize a rolagem horizontal Desse modo, pela condição fornecida no enunciado, y = 1 = k. Logo, a solução particular será: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Para o caso em que u(x) e v(x) forem número reais, isto é, a equação diferencial linear seja de coeficientes constantes, a resolução da equação vai ser um pouco mais simplificada. Vamos supor que u(x) = a e v(x) = b, com a e b reais. Então, teremos Assim: y′ − y = 2ex u(x) =– 1ev(x) = 2ex P (x) = exp(∫ u (x) dx) = exp(∫ (−1) dx) = exp(−x) ∫ P (x) v (x) dx = ∫ e−x 2ex dx = ∫ 2 dx = 2x y = (∫ P (x) v (x) dx + k ) = (2x + k ) = ex (2x + k)1 P(x) 1 e−x y = ex(2x + k) 2y′ − 2y − 4ex = 0 y = ex(2x + k) x = 0 → y = e0 (2.0 + k) = 1 (0 + k) = k y = ex(2x + 1) y′ + ay = b � Atenção! Para visualização completa da equação utilize a rolagem horizontal Logo: , k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal Observe que, se as funções u(x) e v(x) forem contínuas em um domínio S, podemos garantir que a solução obtida pelo procedimento descrito sempre existirá e será única. TEORIA NA PRÁTICA O modelo logístico é mais sofisticado do que o modelo de crescimento exponencial para se trabalhar com crescimentos populacionais. Seja P(t) o tamanho de uma população em um instante t. Define-se a equação diferencial logística como: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Na qual C e K são constante reais. A constante K é denominada de capacidade de suporte. Use o método da equação diferencial separável para obter a solução da equação diferencial logística: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Com a condição inicial de que P(0) = 200. RESOLUÇÃO Precisamos resolver a equação: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Tal equação vale para todo P diferente de zero e de 1000. P (x) = exp(∫ u (x) dx) = exp(∫ a dx) = eax y = (∫ P (x) v (x) dx + k ) = (∫ b eaxdx + k )1 P(x) 1 eax y = (b ∫ eaxdx + k ) = ( eax + k )1 eax 1 eax b a y = + ke−ax b a = CP (1 − )dP dt P K = 0, 05 P (1 − )dP dt P 1000 = 0, 05P (1 − ) → = 0, 05 dtdP dt P 1000 dP P (1 − )P 1000 dP = ( + ) dP = 0, 05dt1000 P(1000 − P) 1 P 1 1000 − P ∫ dP + ∫ dP = ∫ 0, 05 dt1 P 1 1000 − P ln |P | − ln |1000 − P | = 0, 05t + k, k real ln ∣ ∣ ∣ ∣ ∣ ∣ = 0, 05t + k P 1000 − P ∣ ∣ ∣ ∣ ∣ ∣ = exp (0, 05t + k) = exp (0, 05t) exp(k) P 1000 − P � Atenção! Para visualização completa da equação utilize a rolagem horizontal De modo que: � Atenção! Para visualização completa da equação utilize a rolagem horizontal em que , com k real. Resolvendo, teremos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Mas � Atenção! Para visualização completa da equação utilize a rolagem horizontal A solução particular será: � Atenção! Para visualização completa da equação utilize a rolagem horizontal MÃO NA MASSA 1. OBTENHA A SOLUÇÃO DA EQUAÇÃO DIFERENCIAL . A) � Atenção! Para visualização completa da equação utilize a rolagem horizontal B) � Atenção! Para visualização completa da equação utilize a rolagem horizontal C) � Atenção! Para visualização completa da equação utilize a rolagem horizontal D) = Be0,05t P 1000 − P B = ±ek P = 1000Be0,05t − BPe0,05t P(t) = = , com B real 1000Be0,05t 1 + Be0,05t 1000B B + e−0,05t P (0) = = = 200 → 1000B = 200B + 200 → B = = 1000B B + e−0,05.0 1000B B + 1 200 800 1 4 P (t) = = 250 + e−0,05t1 4 1000 1 + 4e−0,05t s′ − cosu + 3u2 + 1 = 0 s = u3 + cosu + k , k real s = u3 − senu − u + k , k real s = −u3 + cosu − 1 + k , k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal E) � Atenção! Para visualização completa da equação utilize a rolagem horizontal 2. SEJA A EQUAÇÃO ATENÇÃO! PARA VISUALIZAÇÃO COMPLETA DA EQUAÇÃO UTILIZE A ROLAGEM HORIZONTAL MARQUE A ALTERNATIVA FALSA RELACIONADA A ESSA EQUAÇÃO DIFERENCIAL. A) É uma equação exata. B) É uma equação não linear. C) É uma equação separável. D) É uma equação de primeira ordem. E) É uma equação de coeficientes variáveis. 3. DESCUBRA A ALTERNATIVA QUE APRESENTA A SOLUÇÃO PARTICULAR PARA EQUAÇÃO DIFERENCIAL , SABENDO QUE, PARA X = 1, O VALOR DE Y VALE 2. A) � Atenção! Para visualização completa da equação utilize a rolagem horizontal B) � Atenção! Para visualização completa da equação utilize a rolagem horizontal C) � Atenção! Para visualização completa da equação utilize a rolagem horizontal D) � Atenção! Para visualização completa da equação utilize a rolagem horizontal E) � Atenção! Para visualização completa da equação utilize a rolagem horizontal s = −u3 + senu − u + k , k real s = −u3 + senu − u +k , k real y′x3 + 3y − 6xy2y′ = 2y3 − 3x2y − 3xy′ y2y′ − x2 − 1 = 0 y = x + 1 y3 = x3 + 3x + 4 xy2 = 2y + 1 y3 = 3x + 5 x2 + y2 = 5 4. DETERMINE A SOLUÇÃO PARTICULAR PARA A EQUAÇÃO DIFERENCIAL ATENÇÃO! PARA VISUALIZAÇÃO COMPLETA DA EQUAÇÃO UTILIZE A ROLAGEM HORIZONTAL E QUE ATENDA À CONDIÇÃO INICIAL DE Y = 1 PARA X = 1. A) � Atenção! Para visualização completa da equação utilize a rolagem horizontal B) � Atenção! Para visualização completa da equação utilize a rolagem horizontal C) � Atenção! Para visualização completa da equação utilize a rolagem horizontal D) � Atenção! Para visualização completa da equação utilize a rolagem horizontal E) � Atenção! Para visualização completa da equação utilize a rolagem horizontal 5. OBTENHA A SOLUÇÃO GERAL DA EQUAÇÃO DIFERENCIAL ATENÇÃO! PARA VISUALIZAÇÃO COMPLETA DA EQUAÇÃO UTILIZE A ROLAGEM HORIZONTAL A) , k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal B) , k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal C) , k real y′x3 + 3y − 6xy2y′ = 2y3 − 3x2y − 3xy′ 3yx + x3y − 2xy3 − 2 = 0 3yx + x3y − 2xy2 − 1 = 0 2yx + x3y + 2y3 − 2 = 0 yx − 2xy3 − 2 = 0 3yx + x3y − 2xy3 + 4 = 0 8x3y + 2y′ − 16x3 = 0 y = ke−x 4 y = 2 + ke−x 4 y = 2x + ke−x 4 � Atenção! Para visualização completa da equação utilize a rolagem horizontal D) , k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal E) , k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal 6. OBTENHA A SOLUÇÃO DA EQUAÇÃO DIFERENCIAL ATENÇÃO! PARA VISUALIZAÇÃO COMPLETA DA EQUAÇÃO UTILIZE A ROLAGEM HORIZONTAL COM X PERTENCENTE AO INTERVALO QUE ATENDA Y(0)=2 A) , k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal B) , K real � Atenção! Para visualização completa da equação utilize a rolagem horizontal C) , k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal D) , k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal E) , k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal y = x3 + kex 3 y = k + 2ex 4 y′ − xsen 2x + y tgx = 0 (− , )π 2 π 2 y = 2 cosx + 2xcos2x + cosx y = 2senx cosx − 2xcos2x + 2cosx y = 2senx cosx − 2cos2x + 2senx y = 2senx cosx + 2xcos2x − 2cosx y = 2senx−2xcosx − 2 GABARITO 1. Obtenha a solução da equação diferencial . A alternativa "D " está correta. Analisando a equação, verificamos que é possível resolvê-la de forma direta. Vejamos: , k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal A alternativa correta é, portanto, a letra D. 2. Seja a equação Atenção! Para visualização completa da equação utilize a rolagem horizontal Marque a alternativa falsa relacionada a essa equação diferencial. A alternativa "C " está correta. Organizando a equação diferencial dada, temos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Repare que a derivada de mais alto grau é de primeira ordem, de modo que a EDO é de primeira ordem. Como os coeficientes dependem de x e y, trata-se de uma equação diferencial com coeficientes variáveis. Por fim, é uma equação não linear, já que aparece um produto de y e y’. Vamos verificar, agora, se a equação diferencial é exata. � Atenção! Para visualização completa da equação utilize a rolagem horizontal Como , trata-se de uma equação diferencial exata. Não é possível colocar a EDO na forma y’=g(x)h(y), assim a equação é não separável, e a alternativa C apresenta uma afirmativa falsa. 3. Descubra a alternativa que apresenta a solução particular para equação diferencial , sabendo que, para x = 1, o valor de y vale 2. s′ − cosu + 3u2 + 1 = 0 s′ − cosu + 3u2 + 1 = 0 → s′ = cosu − 1 − 3u2 = cosu − 1 − 3u2 → ds = (cosu − 1 − 3u2) du ds du ∫ ds = ∫ (cosu − 1 − 3u2) du → s = senu − u − u3 + k y′x3 + 3y − 6xy2y′ = 2y3 − 3x2y − 3xy′ y′x3 + 3y − 6xy2y′ = 2y3 − 3x2y − 3xy′ y′ (x3 − 6xy2 + 3x) = 2y3 − 3y − 3x2y y′ (x3 − 6xy2 + 3x) + (3y + 3x2y − 2y3) = 0 u (x, y) = 3y + 3x2y − 2y3 → = 3 + 3x2 − 6y2 ∂u ∂y v (x, y) = x3 − 6xy2 + 3x → = 3x2 − 6y2 + 3 ∂v ∂x = ∂u ∂y ∂v ∂x y2y′ − x2 − 1 = 0 A alternativa "B " está correta. Analisando a equação, verificamos que se trata de uma equação separável: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Resolvemos passando tudo o que depende do y para um lado e tudo o que depende do x para o outro lado. Vejamos: , C real � Atenção! Para visualização completa da equação utilize a rolagem horizontal Para � Atenção! Para visualização completa da equação utilize a rolagem horizontal Então � Atenção! Para visualização completa da equação utilize a rolagem horizontal Indicado na alternativa B. 4. Determine a solução particular para a equação diferencial Atenção! Para visualização completa da equação utilize a rolagem horizontal e que atenda à condição inicial de y = 1 para x = 1. A alternativa "A " está correta. Veja a solução no vídeo abaixo: 5. Obtenha a solução geral da equação diferencial Atenção! Para visualização completa da equação utilize a rolagem horizontal A alternativa "B " está correta. Organizando a equação diferencial dada, verificamos que se trata de uma equação linear: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Então y2y′ = x2 + 1 → y′ = = g (y)h(x) x2 + 1 y2 y2y′ = x2 + 1 → y2dy = (x2 + 1)dx ∫ y2dy = ∫ (x2 + 1)dx + k → y3 = x3 + x + k1 3 1 3 y3 = x3 + 3x + C x = 1 → y3 = 1 + 3 + C = 23 = 8 → C = 4 y3 = x3 + 3x + 4 y′x3 + 3y − 6xy2y′ = 2y3 − 3x2y − 3xy′ 8x3y + 2y′ − 16x3 = 0 8x3y + 2y′ − 16x3 = 0 → y′ + 4x3y = 68 u(x) = 4x3ev(x) = 8x3 � Atenção! Para visualização completa da equação utilize a rolagem horizontal Agora, temos de obter o fator integrante: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Resolvendo a integral, teremos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Para resolver a integral, vamos fazer , k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal A opção correta é a letra B. 6. Obtenha a solução da equação diferencial Atenção! Para visualização completa da equação utilize a rolagem horizontal Com x pertencente ao intervalo que atenda y(0)=2 A alternativa "B " está correta. Veja a solução no vídeo abaixo: RESOLUÇÃO DE EQUAÇÃO DIFERENCIAL LINEAR DE PRIMEIRA ORDEM GABARITO P (x) = exp(∫ u (x) dx) = exp(∫ 4x3dx) = exp(x4) P (x) = exp(∫ u (x) dx) = exp(∫ 4x3dx) = exp(x4) u = x4 → du = 4x3 2 ∫ 4x3 ex4dx = 2 ∫ eudu = 2eu = 2ex4 y = (∫ P (x) v (x) dx + k ) = (2ex4 + k ) = 2 + ke−x41 P(x) 1 ex 4 y′ − xsen 2x + y tgx = 0 (− , )π 2 π 2 VERIFICANDO O APRENDIZADO 1. MARQUE A ALTERNATIVA QUE APRESENTA A SOLUÇÃO GERAL PARA A EQUAÇÃO DIFERENCIAL ATENÇÃO! PARA VISUALIZAÇÃO COMPLETA DA EQUAÇÃO UTILIZE A ROLAGEM HORIZONTAL A) k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal B) , k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal C) , k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal D) , k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal E) , k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal 2. MARQUE A ALTERNATIVA QUE APRESENTA A SOLUÇÃO PARTICULAR PARA A EQUAÇÃO DIFERENCIAL ATENÇÃO! PARA VISUALIZAÇÃO COMPLETA DA EQUAÇÃO UTILIZE A ROLAGEM HORIZONTAL SABENDO QUE Y(0) = 1. A) � Atenção! Para visualização completa da equação utilize a rolagem horizontal 6x2 − 2ex + 2y′ = 0 y = ln(x) − 2x3 + k y = e2x + k y = x4 + k y = ex − x3 + k y = 2ex − 4x2 + k y′ (3y − 7x − 1) + (2x − 7y + 3) = 0 x2 − xy + 6x + 3y2 − 2y − 1 = 0 B) � Atenção!Para visualização completa da equação utilize a rolagem horizontal C) � Atenção! Para visualização completa da equação utilize a rolagem horizontal D) � Atenção! Para visualização completa da equação utilize a rolagem horizontal E) � Atenção! Para visualização completa da equação utilize a rolagem horizontal GABARITO 1. Marque a alternativa que apresenta a solução geral para a equação diferencial Atenção! Para visualização completa da equação utilize a rolagem horizontal A alternativa "D " está correta. , k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal 2. Marque a alternativa que apresenta a solução particular para a equação diferencial Atenção! Para visualização completa da equação utilize a rolagem horizontal Sabendo que y(0) = 1. A alternativa "B " está correta. Repare que se trata de uma equação exata, pois, 2x2 − 4xy + 6x + 3y2 − 2y − 1 = 0 2x2 − xy + 6x + 3y2 + 2y − 4 = 0 2x2 − 2xy + 6x + 3y2 + 2y − 3 = 0 2x2 − xy − x + y2 − 2y − 1 = 0 6x2 − 2ex + 2y′ = 0 6x2 − 2ex + 2y′ = 0 y′ = = ex − 3x2 dy dx dy = (ex − 3x2) dx ∫ dy = ∫ (ex − 3x2) dx + k → y = ex − x3 + k y′ (3y − 7x − 1) + (2x − 7y + 3) = 0 = ∂u ∂y ∂v ∂x u (x, y) = 2x − 7y + 3 → = −7 ∂u ∂y � Atenção! Para visualização completa da equação utilize a rolagem horizontal Resolvendo a equação diferencial pela solução geral M(x,y) = C, C real. Lembre-se de que e � Atenção! Para visualização completa da equação utilize a rolagem horizontal Vamos, agora, obter a derivada parcial de M em relação a y: � Atenção! Para visualização completa da equação utilize a rolagem horizontal No entanto: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Então: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Desse modo, substituindo, teremos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal A solução da equação diferencial será dada pela equação implícita: , C real � Atenção! Para visualização completa da equação utilize a rolagem horizontal Como y(0) = 1, então: � Atenção! Para visualização completa da equação utilize a rolagem horizontal A solução particular será: � Atenção! Para visualização completa da equação utilize a rolagem horizontal v (x, y) = 3y − 7x − 1 → = −7 ∂v ∂x v (x, y) = 3y − 7x − 1 → = −7 ∂v ∂x (x, y) = v(x, y) ∂M ∂y (x, y) = u (x, y) = 2x − 7y + 3 → M = ∫ (2x − 7y + 3) dx + g(y)∂M ∂x M = ∫ (2x − 7y + 3) dx + g (y) = x2 − 7xy + 3x + g(y) (x, y) = (x2 − 7xy + 3x + g(y)) = −7x + g′(y) ∂M ∂y ∂ ∂y (x, y) = v (x, y) = 3y − 7x − 1 ∂M ∂y −7x + g′ (y) = 3y − 7x − 1 → gprime (y) = 3y − 1 g (y) = ∫ (3y − 1) dy = y2 − y3 2 M = x2 − 7xy + 3x + g (y) = x2 − 7xy + 3x + y2 − y 3 2 M(x, y) = C → x2 − 7xy + 3x + y2 − y = C 3 2 02 − 7.0.1 + 3.0 + 12 − 1 = C → C = 3 2 1 2 x2 − 7xy + 3x + y2 − y = 3 2 1 2 Ou � Atenção! Para visualização completa da equação utilize a rolagem horizontal MÓDULO 4 Reconhecer situações possíveis de serem modeladas por equações diferenciais de primeira ordem MODELAGEM POR MEIO DA EQUAÇÃO DIFERENCIAL DE PRIMEIRA ORDEM MODELAGEM POR EQUAÇÕES DIFERENCIAIS DE PRIMEIRA ORDEM Diversas são as situações — seja nas Ciências, seja na Engenharia — que envolvem uma incógnita e suas taxas de variação representadas por suas derivadas. Esses problemas podem ser modelados por uma equação diferencial. Estudaremos, neste módulo, algumas situações que podem ser modeladas por uma equação diferencial de primeira ordem. Ao analisarmos uma situação prática, buscamos modelar o problema com os dados e as observações realizadas, de forma a equacionar um modelo para buscar a solução. APÓS O LEVANTAMENTO DOS DADOS E DO RELACIONAMENTO ENTRE ELES, MONTAMOS UMA OU MAIS EQUAÇÕES, QUE DEVEM SER RESOLVIDAS. AS EQUAÇÕES DE MODELAMENTO, EM GERAL, SERÃO EQUAÇÕES DIFERENCIAIS, QUANDO AS RELAÇÕES DEPENDEM DA INCÓGNITA A SER OBTIDA E DE SUAS TAXAS DE VARIAÇÃO. Vejamos algumas situações que podem ser modeladas por meio de uma equação diferencial de primeira ordem. QUEDA LIVRE DE UM OBJETO SUJEITO À RESISTÊNCIA DO AR A segunda Lei de Newton relaciona a força, em N, que age em um objeto de massa m, em kg, e sua aceleração em m/s², da seguinte forma: 2x2 − 4xy + 6x + 3y2 − 2y − 1 = 0 � Atenção! Para visualização completa da equação utilize a rolagem horizontal Considerando uma trajetória retilínea, podemos retirar os vetores e trabalhar apenas com os módulos. Devemos lembrar, no entanto, que a velocidade é a primeira derivada da posição em relação ao tempo, já a aceleração é a primeira derivada da velocidade em relação ao tempo. Dessa forma, a aceleração será a segunda derivada da posição pelo tempo: � Atenção! Para visualização completa da equação utilize a rolagem horizontal No caso de um objeto em queda livre, se desprezarmos a resistência do ar, o objeto estará sujeito apenas ao seu peso, e a aceleração será constante e igual à aceleração da gravidade. Nesse caso, não será necessária uma equação diferencial para modelar o problema. Na prática, o ar resiste ao movimento de queda livre, com uma força que é proporcional à sua velocidade, de modo que � Atenção! Para visualização completa da equação utilize a rolagem horizontal K é uma constante de proporcionalidade determinada experimentalmente. Com isso, um objeto em queda livre de massa m, medida em kg, estará sujeito ao peso, que empurra o objeto para baixo, e à resistência do ar, contrária ao peso. Vejamos: em que g é aceleração da gravidade. � Atenção! Para visualização completa da equação utilize a rolagem horizontal Conseguimos modelar o problema da posição em relação ao tempo dessa forma. Porém, temos uma equação diferencial de segunda ordem, que ainda não sabemos resolver. Nesse caso, podemos modelar a velocidade com o tempo: →F = m→a a = = dv dt d2s dt2 Far = Kv = K ds dt FR = P − Far = ma mg − K = m ds dt d2s dt2 � Atenção! Para visualização completa da equação utilize a rolagem horizontal Temos, agora, uma equação diferencial linear de primeira ordem, que já sabemos resolver. Resolvendo a equação diferencial da velocidade pelo método para equação linear, obteremos a seguinte solução: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Após obter a velocidade, usamos a relação: De forma que: Então: � Atenção! Para visualização completa da equação utilize a rolagem horizontal UM OBJETO COM MASSA DE 1 KG ESTÁ EM QUEDA LIVRE EM UM AMBIENTE NO QUAL A CONSTANTE DE PROPORCIONALIDADE DA RESISTÊNCIA DO AR É DE 0,8 NS²/M. O OBJETO SAI DO REPOUSO. DETERMINE A VELOCIDADE MÁXIMA OBTIDA PELO OBJETO, CONSIDERANDO A ACELERAÇÃO DA GRAVIDADE COMO 10 M/S². RESOLUÇÃO O modelo de queda livre será dado pela equação que relaciona a velocidade com o tempo: mg − Kv = m dv dt v (t) = (1 − e− t) m/s mg K K m v = → s = ∫ t 0 v dt ds dt s (t) = ∫ t 0 (1 − e− t) dt mg K K m s (t) = [t − (1 − e− t)]m mg K m k K m mg − Kv = m → + 0, 8v = 10 dv dt dv dt � Atenção! Para visualização completa da equação utilize a rolagem horizontal Trata-se de uma equação diferencial linear com . Então a(t) = 0,8 e b(t) = 10. Agora, temos de obter o fator integrante: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Resolvendo a integral, teremos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Então: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Como � Atenção! Para visualização completa da equação utilize a rolagem horizontal Assim , k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal Quando t tende ao infinito, a exponencial tenderá a zero, e a velocidade chega no seu valor máximo de 12,5 m/s². TAXA DE CRESCIMENTO POPULACIONAL EXPONENCIAL A modelagem do crescimento de uma população — por exemplo, de insetos ou bactérias — é necessária, muitas vezes, em Biologia. Uma forma simplesde modelar o crescimento de uma população é utilizando um modelo exponencial, no qual a taxa de variação da população com o tempo é proporcional ao tamanho da população. Desse modo, seja P(t) a população em um instante de tempo t: em que k é a taxa líquida de aumento da população. � Atenção! Para visualização completa da equação utilize a rolagem horizontal + a (t) v = b (t) dv dt P (t) = exp(∫ a (t) dt) = exp(∫ 0, 8 dt) = e0,8t ∫ P (t) b (t) dt = ∫ e0,8t10 dt = e0,8t = 12, 5e0,8t100 8 v(t) = (∫ P (t) b (t) dt + k ) = (12, 5e0,8t + k )1 P(t) 1 e0,8t v(t) = 12, 5 + ke−0,8tm/s , k real v(0) = 0 → 0 = 12, 5 + ke−0 → k = − 12, 5 v (t) = 12, 5 ( 1 − e−0,8t)m/s = kP dP dt O modelo pode ser empregado para o decrescimento exponencial se k for negativa. Se observarmos a equação do modelo, trata-se de uma equação diferencial linear, de coeficientes constantes e homogênea. Podemos resolver a equação pelo método estudado no módulo anterior. Assim, a solução do modelo será obtida pela resolução dessa equação diferencial linear. Para isso, use o método e você obterá a solução geral: , C REAL Para obter o C, é preciso conhecer uma condição de contorno. Se conhecermos o valor da polução P para t = 0, denominando P0, teremos: Se conhecermos o valor da polução P para t = T, denominando PT, teremos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal ATENÇÃO O modelo de crescimento populacional na forma exponencial é chamado de Lei de Malthus. Vale saber que esse modelo se torna menos preciso para valores grandes de tempo. DECAIMENTO EXPONENCIAL Vejamos, agora, um modelo para um decaimento exponencial, aplicado em um decaimento radioativo, por exemplo. Considere m(t) a massa remanescente de uma substância radioativa em um instante t. Nesse caso, a taxa de decaimento será dada por: em que k é a constante negativa denominada constante de decaimento. � Atenção! Para visualização completa da equação utilize a rolagem horizontal A vida média de um material é definida como o tempo necessário para que metade da massa do material radioativo decaia. P (t) = Cekt P (t) = P0e kt P (t) = PTe k(t−T ) (t) = k m dm dt Determine a vida média de um material sabendo que a constante de decaimento é de 1,4 10-4 por ano. Considere o decaimento exponencial RESOLUÇÃO A equação diferencial que modela o problema será: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Resolvendo a equação � Atenção! Para visualização completa da equação utilize a rolagem horizontal Teremos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Desejamos calcular a vida média do material, que seria o instante em que sua massa chega à metade. Considerando A0 como t = 0 e tv o tempo vida média, temos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Então � Atenção! Para visualização completa da equação utilize a rolagem horizontal MISTURA DE SOLUÇÕES Vamos conhecer agora um problema de química modelado por uma equação diferencial. Um problema de mistura de soluções — por exemplo, a salmoura — está relacionado com um recipiente de capacidade fixa no qual se mistura uma substância, como o sal, em um líquido, como a água. A solução, em dada concentração, entra no recipiente a uma taxa fixa. A mistura, bem agitada no tanque, sai do recipiente também com uma taxa fixa, que pode ser diferente da taxa de entrada (t) = k m dm dt dm = kdt → ∫ dm = ∫ kdt → lnm = kt + C, C real1 m ′ m m = exp (kt + C) = exp (kt) exp (C) = Aekt , A real t = 0 → m = A0 t = tv → m = A0 2 = A0e ktv → ektv = → ktv = ln 2 →tv = ln 2 = ln2 = 4951 anos A0 2 1 2 1 k 1 1, 4 104 Imagem: Jorge Luís Rodrigues Pedreira de Cerqueira A taxa de variação do sal com o tempo será dada por . Essa taxa será determinada pela diferença entre a taxa de entrada, TE, e a taxa de saída, TS, do sal no tanque, ou seja: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Vamos entender esse processo por meio de um exemplo numérico: Seja um recipiente com 5000L de água e 20 kg de sal, inicialmente. Insere-se, no recipiente, uma solução (água salgada) com uma concentração de 0,05 kg de sal por litro de água, a uma taxa fixa de 50L/min. A solução é misturada completamente e tem uma saída do tanque com uma taxa de 50 L/min. Determine a quantidade máxima de sal que permanece no recipiente. RESOLUÇÃO Nosso problema é calcular quanto de sal permanece no tanque depois de certo tempo. Seja s(t) a quantidade de sal, em kg, depois de t minutos. Para t = 0, teremos apenas a quantidade de sal na solução inicial – 20 kg. A taxa de variação do sal com o tempo será dada por . Essa taxa será dada pela diferença entre a taxa de entrada, TE, e a taxa de saída, TS, do sal no tanque, ou seja: � Atenção! Para visualização completa da equação utilize a rolagem horizontal A taxa de entrada, em nosso exemplo, seria dada por 0,05 kg/L vezes 50L/min, de modo que TE = 2,50 kg/min. ds dt = TE − TS ds dt ds dt = TE − TS ds dt Como a taxa de saída é similar à taxa de entrada — no exemplo, de 50L/min —, o recipiente sempre fica com sua capacidade fixa, de 5000L. Considera-se que o sal não aumenta o volume da água. A taxa de saída do sal será de s(t)/5000 kg/L, que mede a quantidade de sal no recipiente pelo volume total, vezes a vazão de saída de 50L/min. Assim, TS = s(t)/100 kg/min � Atenção! Para visualização completa da equação utilize a rolagem horizontal Observe que se trata de uma equação diferencial separável e linear. Com isso, podemos solucionar a equação pelos métodos estudados no módulo anterior: , C REAL � Atenção! Para visualização completa da equação utilize a rolagem horizontal Como t = 0, temos s = 20 kg � Atenção! Para visualização completa da equação utilize a rolagem horizontal De modo que: = 2, 50 − = ds dt s 100 250 − s 100 = dt → ∫ = ∫ dt + C, C realds 250 − s 1 100 ds 250 − s 1 100 −ln |250 − s| = + C t 100 −ln |250 − 20| = + C → C = − ln 230 0 100 −ln |250 − s| = − ln 230 t 100 � Atenção! Para visualização completa da equação utilize a rolagem horizontal Então: , COM T EM MINUTOS � Atenção! Para visualização completa da equação utilize a rolagem horizontal VEJA QUE PODEMOS OBTER A QUANTIDADE DE SAL PARA QUALQUER INSTANTE T. ALÉM DISSO, PODEMOS DETERMINAR A MÁXIMA QUANTIDADE DE SAL NO RECIPIENTE. CONFORME T TENDE PARA INFINITO, A EXPONENCIAL TENDE A ZERO, DE FORMA QUE SMAX = 250KG. CIRCUITOS ELÉTRICOS RL OU RC A resolução de circuitos elétricos com uma resistência e um capacitor em série, denominado RC, ou de um resistor e um indutor em série, denominado RL, podem ser solucionados por meio de uma equação diferencial de primeira ordem. Vejamos: Imagem: Jorge Luís Rodrigues Pedreira de Cerqueira Usando a lei das malhas para o circuito RL, podemos escrever a seguinte equação diferencial: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Conhecendo-se a tensão da fonte v(t) e os componentes do circuito, obtém-se a função i(t), que fornece a corrente elétrica com o tempo. Trata-se de uma equação linear com coeficientes constantes e não homogênea: 250 − s = exp(ln 230 − ) = 230exp(− ) t 100 t 100 s(t) = 250 − 230exp(− )t 100 v(t) − Ri(t) − L = 0 di(t) dt A tensão no indutor pode ser obtida por . � Atenção! Para visualização completa da equação utilize a rolagem horizontal Para o caso do circuito RC, usa-se a lei dos nós do circuito elétrico. O modelo do circuito é obtido pela equação: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Com esse modelo, dada a tensão da fonte v(t) e os elementos do circuito, obtém-se a dependência da tensão no capacitor com o tempo, vc(t). Para se obter a corrente i(t), pode-se fazer . Seja um circuito RL em série com resistência de 15 Ω e indutor de 5 H. A tensão é fornecida por uma fonte contínua de 50 V, que éligada em t = 0s. Determine a corrente após 1 s e a corrente limite do circuito. RESOLUÇÃO Conforme estudamos, o modelo utilizado será dado pela equação diferencial: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Substituindo os dados do problema, temos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Trata-se de uma equação diferencial linear do tipo � Atenção! Para visualização completa da equação utilize a rolagem horizontal Então, a(t) = 3 e b(t) = 10. Iremos, agora, obter o fator integrante: + i (t) = v(t) di(t) dt R L 1 L vL (t) = L di(t) dt = C v(t) − vc(t) R dvc(t) dt i (t) = C dvc(t) dt + i (t) = v(t) di(t) dt R L 1 L + i (t) = 50 → + 3i (t) = 10 di(t) dt 15 5 1 5 di(t) dt + a (t) i (t) = b (t) di (t) dt P (t) = exp(∫ a (t) dt) = exp(∫ 3 dt) = e3t � Atenção! Para visualização completa da equação utilize a rolagem horizontal Resolvendo a integral, temos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Então: , k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal Como circuito é ligado em t = 0s, então i(0) = 0, de modo que: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Então: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Calculando � Atenção! Para visualização completa da equação utilize a rolagem horizontal Quando temos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal JUROS COMPOSTOS CONTINUAMENTE A modelagem também pode ser aplicada na matemática financeira. Considere uma aplicação financeira com um valor inicial A0 aplicado a uma taxa de juros anuais representada por r. O juro do investimento será composto n vezes ao ano. Assim, em cada período de composição, a taxa de juros será dada por uma taxa de . Após período de n vezes o período de composição t, teremos um capital calculado por: ∫ P (t) b (t) dt = ∫ e3t10 dt = e3t 10 3 i(t) = (∫ P (t) b (t) dt + k ) = ( e 3t + k )1 P(t) 1 e3t 10 3 i(t) = + ke−3tA , k real 10 3 0 = + k → k = − 10 3 10 3 i (t) = (1 − e−3t)A 10 3 i(1) = (1 − e−3.1) = 3, 16 A 10 3 t → ∞ lim t→∞ (t) = (1 − 0) = A 10 3 10 3 r n � Atenção! Para visualização completa da equação utilize a rolagem horizontal Por exemplo, seja um investimento inicial de R$ 1.000,00 com uma taxa de 5% ao ano, mas composto de diversas formas. Após 2 anos de investimento, o capital será: COMPOSIÇÃO ANUAL Em um período de dois anos, teremos dois instantes de composição de juros com uma taxa de 5%: � Atenção! Para visualização completa da equação utilize a rolagem horizontal COMPOSIÇÃO SEMESTRAL Em um período de dois anos, teremos quatro instantes de composição de juros com taxa de : � Atenção! Para visualização completa da equação utilize a rolagem horizontal COMPOSIÇÃO MENSAL Em um período de dois anos, teremos 24 instantes de composição de juros com taxa de : � Atenção! Para visualização completa da equação utilize a rolagem horizontal Vamos, agora, considerar que n tende ao infinito. Nesse caso, seria denominado juro composto continuamente da seguinte forma: � Atenção! Para visualização completa da equação utilize a rolagem horizontal O limite dentro dos colchetes é um limite fundamental que vale e, de modo que: A = A0(1 + ) ntr n A = 1.000(1 + ) 2 = 1.102, 50 0, 05 1 5% = 2, 5% 1 2 A = 1.000(1 + ) 4 = 1.103, 81 0, 05 2 5% 1 12 A = 1.000(1 + ) 24 = 1.104, 94 0, 05 12 A = lim n→∞ A0(1 + ) nt = A0 lim n→∞ [(1 + ) ] rt = A0[ lim n→∞ (1 + ) ] rt r n r n n r r n n r javascript:void(0) javascript:void(0) javascript:void(0) � Atenção! Para visualização completa da equação utilize a rolagem horizontal Derivando ambos os lados, em relação ao tempo, obtemos um modelo com taxas de juros compostas continuamente: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Repare que o modelo é uma equação do primeiro grau que diz que a taxa de crescimento do investimento é proporcional ao seu valor. Um capital de R$ 10.000,00 é investido em uma aplicação com taxa de juros de 5% ao ano composta continuamente. Determine a equação que modele o problema e o valor de capital, após dois anos de aplicação. RESOLUÇÃO Conforme estudado o modelo será dado por: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Substituindo os dados do problema, temos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Resolvendo a equação, teremos � Atenção! Para visualização completa da equação utilize a rolagem horizontal No entanto, A(0) = 10.000 = exp (k) A = A0e rt = A0re rt = rA(t) dA dt = A0re rt = rA(t) dA dt = 0, 05 A(t) dA dt = 0, 05 A → = 0, 05 dt → lnA = 0, 05t + k → A (t) = exp (0, 05t) exp(k) dA dt dA A � Atenção! Para visualização completa da equação utilize a rolagem horizontal De forma que: � Atenção! Para visualização completa da equação utilize a rolagem horizontal TRAJETÓRIAS ORTOGONAIS A UMA FAMÍLIA DE CURVAS Considere uma família de curvas no R2. Nesse caso, podemos obter uma trajetória ortogonal, em qualquer ponto pertencente a essas curvas, seguindo o seguinte conceito: Por meio das equações das curvas, obtemos a relação entre y e x Derivamos y em relação a x, obtendo o coeficiente angular, isto é, a inclinação da reta tangente às curvas em qualquer ponto. A inclinação da trajetória ortogonal é obtida pela troca do sinal e pela inversão do coeficiente angular da reta tangente, obtendo-se uma equação diferencial a ser solucionada. Vejamos a situação a seguir: Considere uma família de curvas dadas pela equação , k real. Determine as trajetórias ortogonais a essa família de curvas dadas. RESOLUÇÃO A equação da família de curvas fornece a relação entre y e x, isto é, . Derivando em relação a x, ambos os lados, usando a derivação implícita, temos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal A (t) = 10.000e0,05t A(2) = 10.000e0,05.2 = 10.000e0,1 = 11.051, 71 (x0, y0) y = kx2 y = kx2 = 2kx dy dx Temos de obter uma inclinação que independa de k. Como , temos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Então � Atenção! Para visualização completa da equação utilize a rolagem horizontal Precisamos, agora, resolver a equação diferencial � Atenção! Para visualização completa da equação utilize a rolagem horizontal Usando o método da equação separável, teremos: , C REAL � Atenção! Para visualização completa da equação utilize a rolagem horizontal Observe que são equações de uma família de elipses centradas na origem. TEORIA NA PRÁTICA Você deseja fazer uma aplicação bancária com um capital de R$ 1.000,00. Seu gerente ofereceu um investimento novo que fornecia uma taxa de juros anual de 5% com composição contínua. y = kx2 → k = y x2 = 2 x = 2 dy dx y x2 y x mt = = 2 → mort = = − dy dx y x dy dx x 2y y′ = − → 2yy′ = −x → 2y dy = −xdx x 2y ∫ 2y dy = ∫ −x dx → y2 = − x2 + C → + y2 = C 1 2 x2 2 Para ajudar na sua decisão, o gerente informou que se você investisse com uma composição continua, após dez anos, seu capital seria dez por cento maior, se a composição fosse mensal. Use seus conhecimentos de cálculo para verificar a veracidade da informação dada pelo gerente. RESOLUÇÃO Temos, portanto, um investimento de juros anual de 5%. Se aplicarmos com composição de juros mensal durante cinco anos, teremos 10 x 12 = 120 instantes de composição com juros de : � Atenção! Para visualização completa da equação utilize a rolagem horizontal Para uma aplicação de juros contínua, deve-se calcular uma equação diferencial do tipo: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Substituindo os dados do problema: � Atenção! Para visualização completa da equação utilize a rolagem horizontalResolvendo a equação: � Atenção! Para visualização completa da equação utilize a rolagem horizontal No entanto, A(0) = 1.000 = exp (k) � Atenção! Para visualização completa da equação utilize a rolagem horizontal Com isso, após dez anos, teremos � Atenção! Para visualização completa da equação utilize a rolagem horizontal O valor obtido não é 10% maior do que a composição mensal, que obteve um capital de 1.647,00. MÃO NA MASSA 1. UM CRESCIMENTO POPULACIONAL É MODELADO POR UM CRESCIMENTO EXPONENCIAL. SABE-SE QUE, PARA T = 0, A POPULAÇÃO É DE 1.000 ESPÉCIES E, PARA T = 2 ANOS, É DE 1000E4 ESPÉCIES. DETERMINE A POPULAÇÃO PARA UM INSTANTE DE TEMPO DE 5 ANOS: A) 5% 1 12 A = 1.000(1 + ) 120 = 1.647, 00 0, 05 12 = rA(t) dA dt = 0, 05 A(t) dA dt = 0, 05 A → = 0, 05 dt → lnA = 0, 05t + k → A (t) = exp (0, 05t) exp(k) dA dt dA A A (t) = 1.000e0,05t A(10) = 1.000e0,05.10 = 1.000e0,5 = 1.648, 72 1000 e10 B) C) D) E) 2. UM CAPITAL A0 FOI INVESTIDO EM T = 0 COM UMA TAXA ANUAL DE 10% ANUAIS. ESSE INVESTIMENTO TEVE COMPOSIÇÃO DE JUROS DE FORMA CONTÍNUA. DETERMINE O VALOR DE A0 SABENDO QUE, DOIS ANOS APÓS O INÍCIO DO INVESTIMENTO, O CAPITAL ERA DE 2.000E0,2 A) R$ 1.000,00 B) R$ 2.000,00 C) R$ 3.000,00 D) R$ 4.000,00 E) R$ 5.000,00 3. SEJA UMA FAMÍLIA DE ELIPSES CENTRADAS NA ORIGEM DADAS PELA EQUAÇÃO , COM K REAL. DETERMINE AS TRAJETÓRIAS ORTOGONAIS A ESSA FAMÍLIA DE CURVAS DADAS. A) , k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal B) , k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal C) , k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal D) , k real 2000 e10 3000 e10 4000 e10 5000 e10 x2 + = k y2 2 x = eKy y = Kx2 x = Ky2 x2 + y2 = K � Atenção! Para visualização completa da equação utilize a rolagem horizontal E) , k real � Atenção! Para visualização completa da equação utilize a rolagem horizontal 4. SEJA UM RECIPIENTE COM, INICIALMENTE, 10.000L DE ÁGUA E 50KG DE SAL. INSERE-SE, NO RECIPIENTE, UMA SOLUÇÃO (ÁGUA SALGADA), COM UMA CONCENTRAÇÃO DE 0,1KG DE SAL POR LITRO DE ÁGUA, A UMA TAXA FIXA DE 50L/MIN. A SOLUÇÃO É MISTURADA COMPLETAMENTE E TEM UMA SAÍDA DO TANQUE COM UMA TAXA DE 50L/MIN. DETERMINE A QUANTIDADE MÁXIMA DE SAL QUE PERMANECE NO RECIPIENTE. A) 250 kg B) 500 kg C) 750 kg D) 1.000 kg E) 1.250 kg 5. UM OBJETO CAI EM QUEDA LIVRE A PARTIR DO REPOUSO A UMA ALTURA DE 100M DO SOLO. O OBJETO TEM UMA MASSA DE 5KG. CONSIDERE A CONSTANTE DE RESISTÊNCIA DO AR DE 1 NS²/M E A ACELERAÇÃO DA GRAVIDADE IGUAL A 10M/S². DETERMINE A QUE DISTÂNCIA DO SOLO, APROXIMADAMENTE, O OBJETO ESTARÁ PARA UM INSTANTE DE 4S APÓS O INÍCIO DA QUEDA. A) 8,4 m B) 37,67 m C) 25,7 m D) 32,4 m E) 41,3 m 6. SEJA UM CIRCUITO RC EM SÉRIE COM RESISTÊNCIA DE E CAPACITOR DE 5 F. A TENSÃO É FORNECIDA POR UMA FONTE CONTÍNUA DE 50 V, QUE É LIGADA EM T = 0S. DETERMINE A CORRENTE APÓS 10 S. A) B) C) D) x = Ky3 + 1 10Ω 5e−0,1A 15e−0,5A 15e−0,2A E) GABARITO 1. Um crescimento populacional é modelado por um crescimento exponencial. Sabe-se que, para t = 0, a população é de 1.000 espécies e, para t = 2 anos, é de 1000e4 espécies. Determine a população para um instante de tempo de 5 anos: A alternativa "A " está correta. Como o modelo segue uma modelagem de crescimento exponencial, temos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Então: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Para � Atenção! Para visualização completa da equação utilize a rolagem horizontal Para � Atenção! Para visualização completa da equação utilize a rolagem horizontal Então, Assim, � Atenção! Para visualização completa da equação utilize a rolagem horizontal 2. Um capital A0 foi investido em t = 0 com uma taxa anual de 10% anuais. Esse investimento teve composição de juros de forma contínua. Determine o valor de A0 sabendo que, dois anos após o início do investimento, o capital era de 2.000e0,2 A alternativa "B " está correta. Conforme estudado, o modelo será dado por: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Substituindo os dados do problema, temos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Resolvendo a equação, teremos: 5e−0,5A 5e−0,2A = k P dP dt = kdt → ∫ dP = ∫ k dt + CdP P 1 P lnP = kt + C → P = exp (kt) exp (C) = Aekt, A real t = 0 → P = 1000 = Aet.0 = A → A = 1000 t = 2 → P = 1000ek.2 = 1000e4 → 2k = 4 → k = 2 P (t) = 1000e2t t = 5 → P = 1000 e2.5 = 1000 e10 = rA(t) dA dt = 0, 1 A(t) dA dt � Atenção! Para visualização completa da equação utilize a rolagem horizontal Mas A(0) = A0 = exp (k) � Atenção! Para visualização completa da equação utilize a rolagem horizontal De forma que: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Então A0 = R$ 2.000,00 3. Seja uma família de elipses centradas na origem dadas pela equação , com k real. Determine as trajetórias ortogonais a essa família de curvas dadas. A alternativa "C " está correta. A equação da família de curvas fornece a relação entre y e x, isto é, . Derivando em relação a x, em ambos os lados, usando a derivação implícita, temos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Então � Atenção! Para visualização completa da equação utilize a rolagem horizontal É preciso resolver a equação diferencial � Atenção! Para visualização completa da equação utilize a rolagem horizontal Usando o método da equação separável, temos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Para facilitar a manipulação matemática, vamos transformar a constante C = ln m, com m real: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Então: = 0, 1 A → = 0, 1 dt → lnA = 0, 1t + k → A (t) = exp (0, 1t) exp(k) dA dt dA A A (t) = A0e 0,1t A(2) = A0e 0,1.2 = A0e 0,2 = 2.000 e0,2 x2 + = k y2 2 x2 + = k y2 2 2x + 2y = 0 → 2x + y = 0 → y = −2x dx dx 1 2 dy dx dy dx dy dx mt = = − → mort = = dy dx 2x y dy dx y 2x y′ = → 2xy′ = y y 2x = → ∫ dy = ∫ dx → ln |y| = ln |x| + C → ln y2 = ln |x| + Cy′ 2y 1 x 1 2y 1 x 1 2 ln y2 = ln |x| + C = ln |x| + lnm y2 = m |x| → x = y2 1 m � Atenção! Para visualização completa da equação utilize a rolagem horizontal Fazendo , K real , K real, que corresponde à alternativa C, que são equações de parábolas que passam na origem. 4. Seja um recipiente com, inicialmente, 10.000L de água e 50kg de sal. Insere-se, no recipiente, uma solução (água salgada), com uma concentração de 0,1kg de sal por litro de água, a uma taxa fixa de 50L/min. A solução é misturada completamente e tem uma saída do tanque com uma taxa de 50L/min. Determine a quantidade máxima de sal que permanece no recipiente. A alternativa "D " está correta. O modelo utilizado será A taxa de entrada seria dada por 0,1kg/L vezes 50L/min, então TE = 5,0 kg/min O recipiente sempre fica com sua capacidade fixa, de 10.000L com uma taxa de saída de 50L/min. Assim a taxa de saída do sal será de s(t)/10.000kg/L, que mede a quantidade de sal no recipiente pelo volume total, vezes a vazão de saída de 50L/min. Assim TS = s(t)/200 kg/min � Atenção! Para visualização completa da equação utilize a rolagem horizontal Observe que se trata de uma equação diferencial separável e linear, de forma que podemos solucionar a equação pelos métodos estudados no módulo anterior. Vejamos: , C real , C real � Atenção! Para visualização completa da equação utilize a rolagem horizontal Como t = 0 se tem s = 50kg � Atenção! Para visualização completa da equação utilize a rolagem horizontal Assim: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Então: , com t em minutos = K 1 m → x = Ky2 = TE − TS ds dt = 5, 0 − = ds dt s 200 1000 − s 200 = dt →∫ = ∫ dt + Cds 1000 − s 1 200 ds 1000 − s 1 200 −ln |1000 − s| = + C t 200 −ln |1000 − 50| = + C → C = − ln 950 0 200 −ln |1000 − s| = − ln 950 t 200 1000 − s = exp(ln 950 − ) = 950exp(− )t 200 t 200 s(t) = 1000 − 950exp(− )t 200 � Atenção! Para visualização completa da equação utilize a rolagem horizontal Conforme t tende para infinito, a exponencial tende a zero, de modo que smax = 1.000kg. 5. Um objeto cai em queda livre a partir do repouso a uma altura de 100m do solo. O objeto tem uma massa de 5kg. Considere a constante de resistência do ar de 1 Ns²/m e a aceleração da gravidade igual a 10m/s². Determine a que distância do solo, aproximadamente, o objeto estará para um instante de 4s após o início da queda. A alternativa "B " está correta. Veja a solução no vídeo abaixo: QUEDA LIVRE - APLICAÇÃO EQUAÇÃO DIFERENCIAL PRIMEIRA ORDEM 6. Seja um circuito RC em série com resistência de e capacitor de 5 F. A tensão é fornecida por uma fonte contínua de 50 V, que é ligada em t = 0s. Determine a corrente após 10 s. A alternativa "E " está correta. Veja a solução no vídeo abaixo: CIRCUITO RC - APLICAÇÃO EQUAÇÃO DIFERENCIAL PRIMEIRA ORDEM GABARITO VERIFICANDO O APRENDIZADO 1. UM CRESCIMENTO POPULACIONAL É MODELADO POR UM CRESCIMENTO EXPONENCIAL. SABE-SE QUE, PARA T = 5 ANOS, A POPULAÇÃO É DE 8.000 ESPÉCIES E, PARA T = 0 ANOS, É DE 1.000 ESPÉCIES. DETERMINE O MODELO DO CRESCIMENTO POPULACIONAL COM O TEMPO. APROXIME LN8=2, CASO NECESSITE. A) 10Ω � Atenção! Para visualização completa da equação utilize a rolagem horizontal B) � Atenção! Para visualização completa da equação utilize a rolagem horizontal C) � Atenção! Para visualização completa da equação utilize a rolagem horizontal D) � Atenção! Para visualização completa da equação utilize a rolagem horizontal E) � Atenção! Para visualização completa da equação utilize a rolagem horizontal 2. SEJA UM RECIPIENTE COM, INICIALMENTE, 1.000L DE ÁGUA E 50KG DE SAL. INSERE-SE, NO RECIPIENTE, UMA SOLUÇÃO (ÁGUA SALGADA), COM UMA CONCENTRAÇÃO DE 0,5KG DE SAL POR LITRO DE ÁGUA, A UMA TAXA FIXA DE 20L/MIN. ESSA SOLUÇÃO É MISTURADA COMPLETAMENTE E TEM UMA SAÍDA DO TANQUE COM UMA TAXA DE 20L/MIN. DETERMINE A QUANTIDADE MÁXIMA DE SAL QUE PERMANECE NO RECIPIENTE A) 200kg B) 300kg C) 500kg D) 700kg E) 900kg GABARITO 1. Um crescimento populacional é modelado por um crescimento exponencial. Sabe-se que, para t = 5 anos, a população é de 8.000 espécies e, para t = 0 anos, é de 1.000 espécies. Determine o modelo do crescimento populacional com o tempo. Aproxime ln8=2, caso necessite. A alternativa "D " está correta. Como o modelo segue uma modelagem de crescimento exponencial, temos: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Então: P(t) = 1.000 e t 4 5 P(t) = 500 e t 2 5 P(t) = 1.000 e− t 3 5 P(t) = 1.000 e t 2 5 P(t) = 500 e t 1 5 = k P dP dt = kdt → ∫ dP = ∫ k dt + CdP P 1 P , A real � Atenção! Para visualização completa da equação utilize a rolagem horizontal Para � Atenção! Para visualização completa da equação utilize a rolagem horizontal Para � Atenção! Para visualização completa da equação utilize a rolagem horizontal Então, A alternativa correta é a letra D. 2. Seja um recipiente com, inicialmente, 1.000L de água e 50kg de sal. Insere-se, no recipiente, uma solução (água salgada), com uma concentração de 0,5kg de sal por litro de água, a uma taxa fixa de 20L/min. Essa solução é misturada completamente e tem uma saída do tanque com uma taxa de 20L/min. Determine a quantidade máxima de sal que permanece no recipiente A alternativa "C " está correta. O modelo utilizado será A taxa de entrada seria dada por 0,5 kg/L vezes 20L/min, então, TE = 10,0 kg/min. Como a taxa de saída é similar à taxa de entrada, de 20L/min, o recipiente sempre fica com sua capacidade fixa, de 1.000L. Desse modo, a taxa de saída do sal será de s(t)/1.000 kg/L, que mede a quantidade de sal no recipiente pelo volume total, vezes a vazão de saída de 20L/min. Assim, TS = s(t)/50 kg/min: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Observe que se trata de uma equação diferencial separável e linear, de forma que podemos solucionar a equação pelos métodos estudados no módulo anterior. Vejamos: , C real , C real � Atenção! Para visualização completa da equação utilize a rolagem horizontal Como t = 0, temos s = 50kg � Atenção! Para visualização completa da equação utilize a rolagem horizontal lnP = kt + C → P = exp (kt) exp (C) = Aekt t = 0 → P = 1.000 = Aek.0 → A = 1000 t = 5 → P = 8.000= A ek.5 = 1000 e5k → e5k = 8 5k = ln8 → k = ln 8 = 1 5 2 5 P(t) = 1.000 e t 2 5 = TE − TS ds dt = 10, 0 − = ds dt s 50 500 − s 50 = dt → ∫ = ∫ dt + Cds 500 − s 1 50 ds 500 − s 1 50 −ln |500 − s| = + C t 50 −ln |500 − 50| = + C → C = − ln 450 0 50 Assim: � Atenção! Para visualização completa da equação utilize a rolagem horizontal Então: , com t em minutos � Atenção! Para visualização completa da equação utilize a rolagem horizontal Conforme t tende para infinito, a exponencial tende a zero, de modo que smax = 500kg. A alternativa correta é, portanto, a letra C. CONCLUSÃO CONSIDERAÇÕES FINAIS Neste tema, conhecermos e aplicamos o conceito das equações diferenciais de primeira ordem. Estudamos, inicialmente, os conceitos iniciais e a classificação da equação diferencial, principalmente quanto ao grau e à ordem. Vimos, em seguida, alguns métodos de resolução de equações diferenciais de primeira ordem. Por fim, aplicamos a resolução das equações diferenciais de primeira ordem na modelagem de alguns problemas em diversas áreas. AVALIAÇÃO DO TEMA: REFERÊNCIAS ÇENGEL, Y.; PAUL III, W. J. Equações diferenciais. Porto Alegre: Mc Graw Hill Education, 2012. GUIDORIZZI, H. L. Cálculo. Vol. 4. 5 ed. São Paulo: LTC, 2013. STEWART, J. Cálculo. Vol. 2. 5 ed. São Paulo: Thomson Learning, 2008. THOMAS, G. B. Cálculo. Vol. 1. 12 ed. São Paulo: Pearson, 2013. EXPLORE+ Para saber mais sobre os assuntos tratados neste tema, pesquise: Técnicas de solução de sistemas de equações diferenciais e algébricas: aplicação em sistemas de energia elétrica −ln |500 − s| = − ln 450 t 50 500 − s = exp(ln 450 − ) = 450exp(− )t 50 t 50 s(t) = 500 − 450exp(− )t 50 CONTEUDISTA Jorge Luís Rodrigues Pedreira de Cerqueira CURRÍCULO LATTES javascript:void(0);