EDO de 2ª Ordem com Coeficientes Constantes

Ver Todos os Exercícios

Teoria

Introdução

Basicamente, esse capítulo e os próximos são a aplicação das Transformadas de Laplace, então se você ficou com alguma dúvida volta lá para dar uma revisada!

Você já viu a definição da Transformada de Laplace, como calculá-la e várias de suas propriedades. Mas para que isso tudo? Para chegarmos aqui: vamos usar todos esses conceitos para resolver EDO’S.

Transformada de Laplace de derivadas

Antes de continuar, é bom lembrar bem das transformadas de derivadas. Vimos anteriormente que a transformada de uma derivada de ordem n é:

L f n t = s n L f t - s n - 1 f 0 - s n - 2 f ' 0 - … - s f n - 2 0 - f n - 1 ( 0 )

Traduzindo isso aí, a primeira coisa que vai aparecer é um s elevado a ordem da derivada e que multiplica a transformada da função:

s n L f t

Depois, todos os termos terão o sinal de menos que pode sumir na hora de substituir os valores.

- s n - 1 f 0 - s n - 2 f ' 0 - … - s f n - 2 0 - f n - 1 ( 0 )

O grau do s vai diminuindo um a um a cada termo até chegar ao grau zero e a ordem de f 0 vai aumentando um a um a cada termo até chegar ao grau n - 1 .

Não vai esquecer isso na prova, hein?! Para ordem 2 , vai ficar:

L f ' ' t = s 2 L f t - s f 0 - f ' ( 0 )

E ordem 1 :

L f ' t = s L f t - f ( 0 )

Repare que precisamos saber os valores iniciais da função para calcular.

Entendeu tudo isso? Vamos seguir. Não entendeu? Volte no tipo de transformadas de derivadas que está tudo bem mais suave por lá.

EDO e Laplace

Vamos ver um exemplo prático desse método, para você entender bem como ele funciona.

Exemplo: Resolva o problema de valor inicial:

y ' ' - y ' - 2 y = 0

y 0 = 1

y ' 0 = 0

Provavelmente você já resolveu esse tipo de problema pelo método da equação característica, agora vamos ver como resolver por Laplace.

Passo 1: Aplicar a Transformada de Laplace à equação.

Lembra que a Transformada de Laplace é linear? Temos

L y ' ' - y ' - 2 y = L y ' ' - L y ' - 2 L y = 0

Agora vamos aplicar aquela propriedade da Transformada da derivada, lembra dela?

L y ' = s L y - y 0

L y ' ' = s 2 L y - s y 0 - y ´ ( 0 )

Substituindo na equação, temos:

s 2 L y s - s y 0 - y ´ 0 - s L y - y 0 - 2 L y = 0

Passo 2: Substituir os valores iniciais e isolar L y

Como y 0 = 1   e   y ' 0 = 0 , temos:

s 2 L y - s - s L y - 1 - 2 L y = 0

s 2 L y - s - s L y + 1 - 2 L y = 0

L y s 2 - s - 2 = s - 1

L y = s - 1 s 2 - s - 2

Passo 3: Agora vamos decompor a expressão que encontramos em frações parciais.

O tipo que trata das transformadas tradicionais fala bem de frações parciais, se você não entender esse passo, volte rapidinho nele ;-)

As raízes de ( s 2 - s - 2 ) são 2 e - 1 , portanto, podemos escrever s 2 - s - 2 = ( s - 2 ) ( s + 1 ) .

s - 1 s 2 - s - 2 = a s - 2 + b s + 1

Temos que juntar tudo agora para achar essas constantes a e b :

s - 1 s 2 - s - 2 = a s + 1 + b ( s - 2 ) s - 2 ( s + 1 ) = a s + a + b s - 2 b s 2 - s - 2

Juntando os termos pelos graus de s :

L y = s - 1 s 2 - s - 2 = s a + b + ( a - 2 b ) s 2 - s - 2

Como o denominador já é o mesmo, temos:

s - 1 = s a + b + ( a - 2 b )

Portanto:

a + b = 1         ⟹         a = 1 - b a - 2 b = - 1                                                        

1 - b - 2 b = - 1         ⟹             b = 2 / 3

a = 1 - 2 / 3         ⟹             a = 1 / 3

  1. Agora que já temos os coeficientes a e b , basta reescrever as frações que criamos.

L y = s - 1 s 2 - s - 2 = 1 / 3 s - 2 + 2 / 3 s + 1

Passo 4: Agora que já temos as frações parciais, precisamos reconhecer a função y pela sua transformada. Ou seja, achar a transformada inversa!

Lembra que L e a t = 1 / ( s - a ) ? Essa transformada é muito parecida com cada uma das frações parciais, não é? Vamos, então, reescrever essas frações para reconhecer a função.

L y = 1 3 ⋅ 1 s - 2 + 2 3 ⋅ 1 s + 1

L y = 1 3 L e 2 t + 2 3 L e - t

L y = L 1 3 e 2 t + 2 3 e - t

Portanto, y é exatamente aquela expressão dentro dos colchetes da transformada:

y ( t ) = 1 3 e 2 t + 2 3 e - t

Lembre-se que y é uma função de t , não de s ! Não se confunda com isso, ok?

Talvez você esteja se perguntando: por que resolver EDO por Laplace? O método antigo já não estava bom?

Realmente, em certos casos resolver por Laplace dá muito mais trabalho! Mas Laplace também tem suas vantagens.

Uma delas é que conseguimos trabalhar com EDO’s não-homogêneas muito facilmente (se você não lembra, não tem problema, EDO não-homogênea é aquela que não tem “ 0 ” do outro lado do “ = ”. Tem uma função f ( t ) qualquer perdida).

Nesse caso, a única diferença é que temos que calcular a transformada dessa função também f ( t ) , quando formos aplicar Laplace na EDO. Outra vantagem é que não é necessário substituir as condições iniciais depois, ao longo do método elas já serão substituídas (lembra que na transformada da derivada, é necessário conhecer os valores iniciais?).

Resolver esse tipo de EDO não homogênea pelo método da equação característica dá muito mais trabalho! Se você já aprendeu isso, deve se lembrar.

Então, vamos resumir o que vimos para organizar melhor as ideias!

Passo a Passo para EDO por Laplace:

  1. Aplicar Laplace à EDO;
  2. Isolar L y , encontrando uma expressão para a transformada de y (não se esqueça de substituir os valores iniciais que o problema deu);
  3. Decompor a expressão de L y em frações parciais;
  4. Reconhecer funções por suas transformadas (ver e “ajustar” cada fração parcial separadamente).

Basicamente todos os exercícios de EDO resolvidos por Laplace vão seguir esse mesmo padrão. A parte mais complicada será reconhecer a função y depois que você encontrar sua transformada L y , às vezes você vai precisar usar um pouco de algebrismo. Por isso, é uma boa você treinar bastante essa “criatividade” fazendo exercícios!

EDO e Laplace

Exercícios Resolvidos

Exercício Resolvido #1

Resolva o problema de valor inicial

y ' ' + 4 y = 2 t ,     y 0 = 1 ,     y ' 0 = 0

Passo 1

Vamos usar a Transformada de Laplace para resolver esse PVI, beleza? Para isso, aplicamos a transformada de Laplace em todo mundo.

L y ( t ) = Y ( s )

L y ' ' t = s 2 Y s - s y 0 - y ' 0

Agregando os valores iniciais dados, temos

L y ' ' t = s 2 Y s - s

Além disso, para o termo do lado direito, sabemos que

L 2 t = 2 s 2

Passo 2

Substituindo na equação, temos:

y ' ' + 4 y = 2 t

L [ y ' ' ] + L [ 4 y ] = L [ 2 t ]

s 2 Y s - s + 4 Y s = 2 s 2

Reordenando:

Y s s 2 + 4 = s 3 + 2 s 2

Y s = s 3 + 2 s 2 s 2 + 4

Passo 3

Agora, vamos usar frações parciais para achar a inversa de Y ( s ) , isto é, o y ( t ) .

s 3 + 2 s 2 s 2 + 4 = A s + B s 2 + C s + D s 2 + 4

s 3 + 2 s 2 s 2 + 4 = A s + B s 2 + 4 + C s + D s 2 s 2 s 2 + 4

O numerador da direita, então, se torna:

A s + B s 2 + 4 + C s + D s 2 = A s 3 + B s 2 + 4 A s + 4 B + C s 3 + D s 2

A s + B s 2 + 4 + C s + D s 2 = ( A + C ) s 3 + ( B + D ) s 2 + 4 A s + 4 B

E, comparando termo a termo, para que a expressão se iguale a s 3 + 2 :

A + C = 1 B + D = 0 4 A = 0 4 B = 2

Assim, concluímos que A = 0 , B = 1 / 2 , C = 1 e D = - 1 / 2 . Logo:

s 3 + 2 s 2 s 2 + 4 = 1 / 2 s 2 + s - 1 / 2 s 2 + 4

s 3 + 2 s 2 s 2 + 4 = 1 / 2 s 2 + s s 2 + 4 - 1 / 2 s 2 + 4

Passo 4

Por fim, vamos reconhecer as inversas...

Y ( s ) = 1 / 2 s 2 + s s 2 + 4 - 1 / 2 s 2 + 4

Como y t =   L - 1 Y s , vem:

y t = L - 1 1 2 . 1 s 2 + L - 1 s s 2 + 4 - L - 1 1 4 . 2 s 2 + 4

y t = t 2 + cos ⁡ ( 2 t ) - s e n   ( 2 t ) 4

Resposta

y t = t 2 + cos ⁡ ( 2 t ) - s e n   ( 2 t ) 4

Exercício Resolvido #2

Considere o seguinte PVI:

x ' ' + 6 x ' + 18 x = 25 e t ,   x 0 = 1 ,   x ' 0 = 4

Resolva usando transformada de Laplace.

Passo 1

Vamos lembrar as transformadas das derivadas:

L x ' s = s L x s - x 0 = s L x s - 1

L x ' ' s = s 2 L x s - s x 0 - x ´ 0 = s 2 L x s - s - 4

Em breve usaremos isso!

Vamos aplicar a transformada em ambos os lados da EDO:

L x ' ' + 6 x ' + 18 x = 25 e t L x ' ' + 6 L x ' + 18 L x = L [ 25 e t ]

E a transformada de uma exponencial:

L 25 e t = 25 s - 1

Substituindo as transformadas:

s 2 L x - s - 4 + 6 s L x - 6 + 18 L x = 25 s - 1

Passo 2

Agrupando tudo que tem L [ x ] e colocando o resto para outro lado da equação:

s 2 + 6 s + 18 L x = 25 s - 1 + s + 10 = 25 + s 2 + 10 s - s - 10 s - 1 = s 2 + 9 s + 15 s - 1

Assim:

L x = s 2 + 9 s + 15 s - 1 s 2 + 6 s + 18

Passo 3

Vamos achar as raízes de s 2 - 6 s + 18 :

∆ = 36 - 4.18 = - 36 < 0

Eita, não deu galera! As raízes vão dar números complexos! Então vamos só completar quadrado perfeito:

s 2 + 6 s + 18 = s 2 + 6 s + 9 + 9 = s + 3 2 + 9

Substituindo isso na EDO:

L x = s 2 + 9 s + 15 s - 1 s + 3 2 + 9

Passo 4

Vamos aplicar frações parciais :

s 2 + 9 s + 15 s - 1 s + 3 2 + 9 = A s - 1 + B s + C s 2 + 6 s + 18

Galera, eu voltei a s 2 + 6 s + 18 ao invés de s + 3 2 + 9 para fazer o mmc!

A s - 1 + B s + C s 2 + 6 s + 18 = A s 2 + 6 A s + 18 A + B s 2 - B s + C s - C s - 1 s + 3 2 + 9

Assim:

s 2 + 9 s + 15 s - 1 s + 3 2 + 9 = A + B s 2 + 6 A - B + C s + 18 A - C s - 1 s + 3 2 + 9

Igualando os coeficientes:

A + B = 1 → B = 1 - A 6 A - B + C = 9 → 6 A - 1 + A + C = 9 → 7 A + C = 10 18 A - C = 15

Analisando a segunda e a terceira equação:

7 A + C = 10

18 A - C = 15

Somando esses 2 termos:

25 A = 25

A = 1

Assim:

C = 10 - 7 = 3

E

B = 1 - 1 = 0

Logo ficou:

s 2 + 9 s + 15 s - 1 s + 3 2 + 9 = 1 s - 1 + 3 s + 3 2 + 9

Então finalmente:

x = L - 1 s 2 + 9 s + 15 s - 1 s + 3 2 + 9 = L - 1 1 s - 1 + 3 s + 3 2 + 9 = L - 1 1 s - 1 + L - 1 3 s + 3 2 + 9

Passo 5

Então vamos resolver cada um:

L - 1 1 s - 1 = e t

L - 1 3 s + 3 2 + 9 = ? ?

Esse s + 3 representa que a função s está deslocada em - 3 , o macete é esquecer ela por um momentinho e ficar só com: L - 1 3 s 2 + 9 = sen ⁡ ( 3 t )

Beleza, agora lembre que era está deslocada, isso representa uma exponencial, em que a constante do expoente é o deslocamento da função s , tá lá na nossa teoria ;), então fica assim:

L - 1 3 s + 3 2 + 9 = e - 3 t sen ⁡ ( 3 t )

E finalmente a solução:

x = L - 1 1 s - 1 + L - 1 3 s + 3 2 + 9 = e t + e - 3 t sen ⁡ ( 3 t )

Resposta

x = e t + e - 3 t sen ⁡ ( 3 t )

Exercício Resolvido #3

Considere o PVI:

2 y ' ' + 5 y ' + 2 y = cos ⁡ ( 3 t ) ,   y 0 = 1 ,   y ' 0 = 4

a) Ache a transformada de Laplace da solução deste PVI.

b) Expresse a solução particular y ( t ) do PVI em função de uma convolução

Passo 1

a)

Utilizando a propriedade de Laplace, transformada da derivada:

L y ' s = s L y s - y 0

L y ' ' s = s 2 L y s - s y 0 - y ´ ( 0 )

Com as condições iniciais

L y ' = s L y - 1

L y ' ' = s 2 L y - s - 4

Aplicando a transformada na EDO:

L 2 y ' ' + 5 y ' + 2 y = L [ cos ⁡ 3 t ]

Mas sabemos que:

L 2 y ' ' + 5 y ' + 2 y = 2 L y ' ' + 5 L y ' + 2 L [ y ]

Substituindo as transformadas:

L 2 y ' ' + 5 y ' + 2 y = 2 s 2 L y - s - 4 + 5 s L y - 1 + 2 L [ y ]

Arrumando:

L 2 y ' ' + 5 y ' + 2 y = L y 2 s 2 + 5 s + 2 - 2 s - 13

Passo 2

Agora o lado direito:

L cos ⁡ 3 t = s s 2 + 9

Assim a EDO fica:

2 s 2 + 5 s + 2 L y s - 2 s - 13 = s s 2 + 9

Jogando os termos sem Laplace para o lado direito:

2 s 2 + 5 s + 2 L y s = s s 2 + 9 + 2 s + 13

Vou dividir tudo por 2 :

s 2 + 5 2 s + 1 L y s = 1 2 s s 2 + 9 + s + 13 2

Colocando L [ y ] ( s ) em evidência:

L y s = 1 2 s s 2 + 9 s 2 + 5 2 s + 1 + s + 13 2 s 2 + 5 2 s + 1

Eita que doidera, né?! Mas a letra a já foi!

Passo 3

b) Agora pra letra b , vamos pegar o resultado da a e arrumar de novo:

L y s = 1 2 s s 2 + 9 s 2 + 5 2 s + 1 + s + 13 2 s 2 + 5 2 s + 1

L y s = 1 s 2 + 5 2 s + 1 1 2 s s 2 + 9 + s + 13 2

Porque isso ajuda a gente? Bem, porque podemos usar a propriedade distributiva da integral de convolução!

f * g 1 + g 2 = f * g 1 + f * g 2

L f * g 1 + g 2 = L f ∙ L g 1 + g 2

L f * g 1 + g 2 = L f ∙ L g 1 + L g 2

Dessa forma, temos que:

L f t = 1 s 2 + 5 2 s + 1

L g 1 t = 1 2 s s 2 + 9

L g 2 t = s + 13 2

Mas opa, dois desses aí são tranquilos, vamos ter:

L - 1 1 2 s s 2 + 9 = g 1 t = 1 2 cos ⁡ ( 3 t )

L - 1 s + 13 2 = g 2 t = L - 1 s + L - 1 13 2 = δ ' t + 13 2 δ t

O que sobrou a gente tenta achar a transformada inversa também

L - 1 1 s 2 + 5 2 s + 1 = L - 1 1 s + 1 2 s + 2

Utilizando frações parciais

1 s + 1 2 s + 2 = A s + 1 2 + B s + 2

Assim obteremos o seguinte sistema

A + B = 0 2 A + B 2 = 1

E encontraremos que

A = 2 3   ,     B = - 2 3    

E assim,

L - 1 1 s + 1 2 s + 2 = L - 1 2 3 s + 1 2 - 2 3 ( s + 2 ) = 2 3   L - 1 1 s + 1 2 - 1 s + 2

E desses dois dentro dos colchetes nós sabemos a transformada inversa.

Com isso podemos concluir que:

L - 1 1 s 2 + 5 2 s + 1 = 2 3   L - 1 1 s + 1 2 - 1 s + 2 = 2 3 e - t 2 - 2 3 e - 2 t

Chegamos ao fim, vamos ter:

f t = 2 3 e - t 2 - 2 3 e - 2 t

g t = g 1 t + g 2 t = 1 2 cos ⁡ ( 3 t ) + δ ' t + 13 2 δ t

E dá uma olhada nisso: Tá vendo que temos um produto de transformadas? E nós sabemos quais funções deram origem a essas transformadas! Então o que temos que fazer é só escrever na forma da definição de convolução igual aqui em baixo!

L y s = L f * g 1 + g 2 = L f * g

y p = f * g = ∫ 0 t f τ   g t - τ d τ

y p = ∫ 0 t 2 3 e - τ 2 - 2 3 e - 2 τ 1 2 cos ⁡ ( 3 ( t - τ ) ) + δ ' ( t - τ ) + 13 2 δ ( t - τ ) d τ

y p = ∫ 0 t 1 3 e - τ 2 - 1 3 e - 2 τ cos ⁡ ( 3 ( t - τ ) ) + 2 δ ' ( t - τ ) + 13 δ ( t - τ ) d τ

y p = 1 3 ∫ 0 t e - τ 2 - e - 2 τ cos ⁡ ( 3 t - 3 τ ) + 2 δ ' ( t - τ ) + 13 δ ( t - τ ) d τ

Pronto!

Resposta

a )   L y s = 1 2 s s 2 + 9 s 2 + 5 2 s + 1 + s + 13 2 s 2 + 5 2 s + 1

b )   y p = 1 3 ∫ 0 t e - τ 2 - e - 2 τ cos ⁡ ( 3 t - 3 τ ) + 2 δ ' ( t - τ ) + 13 δ ( t - τ ) d τ

Exercício Resolvido #4

Use transformadas de Laplace para resolver o seguinte problema do valor inicial

y ' + y = - t sen ⁡ t

Com y 0 = 0

Passo 1

Vamos aplicar transformada de Laplace na equação

L y ' + y = L - t sen ⁡ t

L y ' + L y = L - t sen ⁡ t

Lembrando que

L y ' = s L y - y 0 = s L y

Vamos fazer

L y = Y s

Para facilitar nossa vida, a EDO fica

s Y s + Y s = - L t sen ⁡ t

Passo 2

Vamos fazer agora

L t sen ⁡ t

Lembrando que

L t n f t = - 1 n d n d s n L f t

E

L sen ⁡ a t = a s 2 + a 2

Vamos ter

L t sen ⁡ t = - 1 1 d d s L sen ⁡ t = 2 s ( s 2 + 1 ) 2

Passo 3

A EDO ficou

s Y s + Y s = - 2 s ( 1 + s 2 ) 2

Isolando Y s vamos ter

Y s = - 2 s s 2 + 1 2 . 1 s + 1 = - 2 s s + 1 ( s 2 + 1 ) 2

Agora a gente vai ter um trabalhinho porque vamos precisar separar isso em frações parciais:

Y s = - 2 s s + 1 ( s 2 + 1 ) 2 = A s + 1   + B s + C s 2 + 1 + D s + E s 2 + 1 2

Essa forma de separar as frações parciais é ótima porque temos 5 variáveis ( A até E ) e vamos montar um sistema com 5 equações (de s 4 até s 0 ). E lembre-se, o numerador em frações parciais sempre tem que ser pelo menos um grau a menos que o denominador. Então colocando tudo no mesmo denominador temos:

- 2 s s + 1 ( s 2 + 1 ) 2 = A s 2 + 1 2 s + 1   + B s + C s + 1 s 2 + 1 s 2 + 1 + D s + E s + 1 s 2 + 1 2

- 2 s = A s 4 + 2 s 2 + 1 + B s + C s 3 + s 2 + s + 1 + D s 2 + E s + D s + E

- 2 s = A s 4 + 2 A s 2 + A + B s 4 + B s 3 + B s 2 + B s + C s 3 + C s 2 + C s + C + D s 2 + E s + D s + E

0 s 4 + 0 s 3 + 0 s 2 + - 2 s + 0 = A + B s 4 + B + C s 3 + 2 A + B + C + D s 2 + B + C + E + D s + A + C + E

Igualando os índices montamos o seguinte sistema:

A + B = 0 B + C = 0 2 A + B + C + D = 0 B + C + E + D = - 2 A + C + E = 0

Da linha 2 , conseguimos simplificar a linha 3 e 4 :

A + B = 0 B + C = 0 2 A + D = 0 E + D = - 2 A + C + E = 0

Fazendo linha 1 menos linha 2 achamos A = C . Substituindo todos os C por A temos:

A + B = 0 B + A = 0 2 A + D = 0 E + D = - 2 A + A + E = 0

A + B = 0 2 A + D = 0 E + D = - 2 2 A + E = 0

Reduzimos em uma equação porque a primeira e segunda linha eram iguais. Do sistema acima, podemos fazer linha 2 menos linha 4 pra obter D = E e trocando D por E na terceira linha, ficamos com:

E + E = - 2

2 E = - 2

E = - 1 = D

E daí conseguimos pegar os valores das constantes que faltavam. Por exemplo, do ultimo sistema, linha 4 , conseguimos obter:

2 A + E = 0

2 A - 1 = 0

A = 1 2 = C

E por fim da linha 1 do último sistema:

A + B = 0

1 2 + B = 0

B = - 1 2

Passo 4

Do resultado do passo anterior, já que:

Y s = A s + 1   + D s + E s 2 + 1 2 + B s + C s 2 + 1

A = 1 2 B = - 1 2 C = 1 2 D = - 1 E = - 1

Vamos ter:

Y s = 1 2 1 s + 1 - s + 1 s 2 + 1 2 - 1 2 ⋅ s - 1 s 2 + 1

Assim

y t = L - 1 Y ( s ) = 1 2 L - 1 1 s + 1 - L - 1 s + 1 s 2 + 1 2 - 1 2 L - 1 s - 1 s 2 + 1

A primeira é a transformada de 1 deslocada em 1 unidade.

L - 1 1 s + 1 = e - t

A segunda é dada por

L - 1 s + 1 s 2 + 1 2 = L - 1 s s 2 + 1 2 + L - 1 1 s 2 + 1 2

L - 1 s + 1 s 2 + 1 2 = 1 2 t s e n t + 1 2 s e n t - t c o s t

A terceira, a gente separa em duas:

L - 1 1 - s s 2 + 1 = L - 1 1 s 2 + 1 - L - 1 s s 2 + 1

L - 1 1 - s s 2 + 1 = s e n t - cos ⁡ t

Que são as transformadas do cosseno e seno respectivamente. Assim juntando tudo isso:

y t = L - 1 Y ( s ) = 1 2 e - t - 1 2 t s e n t + 1 2 s e n t - t c o s t - 1 2 s e n t - cos ⁡ t    

Simplificando,

y t = L - 1 Y ( s ) = 1 2 e - t - t s e n t - 2 s e n t - cos ⁡ t

Resposta

y t = L - 1 Y ( s ) = 1 2 e - t - t s e n t - 2 s e n t - cos ⁡ t