Coeficientes Constantes Não Homogênea-Variação dos Parâmetros

Ver Todos os Exercícios

Teoria

Introdução:

Faaaaaaala tu, tatu! Tudo bem? Meu anjo, a gente já viu um bocado de coisa de EDO superior né? Já vimos o caso homogêneo e o não homogêneo, mas ainda precisamos ver um outro caso não homogêneo.

Na real, esse é um método de resolver EDO superior não homogênea, beeeeeem similar ao que vimos em EDO de 2 ª ordem! Tá, é uma extensão, vai!

Agora vamos ver como resolver EDO’s de ordem superior com coeficientes constantes não homogênea. Elas têm a seguinte forma:

a n y n + a n - 1 y n - 1 + … a 3 y ' ' ' + a 2 y ' ' + a 1 y ' + a 0 y = f ( x )

Com todos a 0 , a 1 , a 2 … a n constantes.

Nesse caso, precisamos fazer como fazíamos lá em EDO de 2ª ordem não homogênea. Ou seja, resolver a EDO homogênea associada, e achar uma solução particular dessa EDO com f x . A solução vai ser:

y = y h + y p

Tu lembra como se resolve EDO não homogênea? Aqui uma colinha:

Encontrando a Solução Geral:

Agora vamos ver como encontrar a solução de EDO’s de ordem superior com coeficientes constantes não homogênea.

O processo será:

  1. Resolver a EDO como se fosse homogênea – solução homogênea: y h x
  2. Encontrar a solução particular: y p x
  3. Encontrar solução geral: y x = y h x + y p x

Bora começar com um exemplo? Se liga

Exemplo:

Considere a seguinte EDO:

y ' ' ' + 4 y ' = s e c ⁡ ( 2 x )

Vamos encontrar a solução da EDO homogênea associada:

y ' ' ' + 4 y ' = 0

A equação auxiliar é:

r 3 + 4 r = 0

Colocando r em evidência temos:

r r 2 + 4 = 0

Temos:

r 1 = 0

Ou

r 2 + 4 = 0

Então:

r 2 = - 4

Vamos ter:

r 2 = 2 i

r 3 = - 2 i

Portanto as raízes são:

r 1 = 0 ,   r 2 = + 2 i , r 3 = - 2 i

Portanto, a solução da equação homogênea associada é:

y h = c 1 + c 2 cos ⁡ 2 x + c 3 s e n ( 2 x )

Agora vamos encontrar a solução particular da equação não homogênea. Você lembra como fazíamos isso em EDO de 2 ª ordem? Dá uma olhadinha nisso aqui, a gente já volta pro nosso exemplo

Relembrando y p :

Agora vamos relembrar como fazíamos para encontrar y p para EDO’s de 2ª Ordem com não homogênea pelo Método da Variação dos Parâmetros. Só para recordar, a gente usava esse método quando f ( x ) não era um produto/soma das funções trigonométricas, exponenciais ou polinomiais.

A EDO era da forma:

y ' ' x + p x y ' x + q x y x = B x

Onde p , q , B são funções contínuas.

Nós supomos que:

y p x = f 1 x y 1 x + f 2 x y 2 ( x )

Onde y 1 ( x ) e y 2 ( x ) eram soluções da EDO homogênea associada.

Nós derivamos y p ( x ) , substituímos na EDO e após algumas contas achamos como calcular f 1 ( x ) e f 2 ( x ) :

f 1 x = ∫ y 2 x . B x M d x

f 2 x = ∫ - y 1 x . B x M d x

Onde M = y 1 ' y 2 ' y 1 y 2 = y 1 ' x . y 2 x - y 2 ' x y 1 x

Logo, a gente encontrava y p .

Vamos ver nessa seção que, o que fazíamos para EDO’s de ordem 2, é na verdade, um caso específico. Ou seja, vamos ver que na verdade o método da variação de parâmetros pode ser usado para EDO’s de ordem superiores. Vamos lá!!!

Voltando ao nosso caso...

Variação dos Parâmetros: Ordem Superior

Aqui vai ser um pouco parecido com antes, só que as contas vão ser um bocado maiores.

Imagina que temos aquela EDO não homogênea, a primeira coisa a fazer é achar a solução da EDO homogênea associada.

Com isso vamos achar y 1 ,   y 2 ,   y 3 ,   y 4 , … , y n , para a conta da variação de parâmetros vamos considerar essas soluções sem as constantes, por exemplo, antes tínhamos:

y 1 = C 1 e r x

Para a variação de parâmetros vamos considerar

y 1 = e r x

Beleza?

A nossa solução particular vai ser da seguinte forma

y p = y 1 u 1 + y 2 u 2 + y 3 u 3 + … + y n - 1 u n - 1 + y n u n

Precisamos achar esses u ,que nada mais são que funções de x . A conta é uma formulinha, um pouco trabalhosa, mas através dela conseguimos resolver o problema.

No nosso exemplo, o método supõe que:

y p = u 1 y 1 + u 2 y 2 + u 3 y 3

Sendo y 1 = 1 ,   y 2 = c o s ⁡ ( 2 x ) e y 3 = s e n ( 2 x )

Então:

y p = u 1 . 1 + u 2 cos ⁡ 2 x + u 3 s e n ( 2 x )

Beleza, agora pra gente chegar a algum lugar com isso, precisaremos introduzir umas novas fórmulinhas. Vamos primeiro considerar a matriz

W = y 1 y 2 y 3 … y n y 1 ' y 2 ' y 3 ' … y n ' y 1 ' ' y 2 ' ' y 3 ' ' … y n ' ' … … … … … y 1 n - 1 y 2 n - 1 y 3 n - 1 … y n n - 1

No nosso caso, essa matriz é:

W = y 1 y 2 y 3 y 1 ' y 2 ' y 3 ' y 1 " y 2 " y 3 "

Calculando as derivadas:

y 1 = 1

y 1 ' = 0

y 1 " = 0

Temos também:

y 2 = cos ⁡ 2 x

y 2 ' = - 2 s e n 2 x

y 2 " = - 4 c o s ⁡ ( 2 x )

Falta ainda:

y 3 = s e n 2 x

y 3 ' = 2 cos ⁡ 2 x

y 3 " = - 4 s e n ( 2 x )

Substituindo em W temos:

W = 1 cos ⁡ 2 x s e n 2 x 0 - 2 s e n 2 x 2 cos ⁡ 2 x 0 - 4 c o s ⁡ ( 2 x ) - 4 s e n ( 2 x )

Calculando W temos:

W = 1 - 2 s e n 2 x - 4 s e n 2 x - - 4 cos ⁡ 2 x 2 cos ⁡ 2 x - cos ⁡ 2 x 0 - 0 + s e n ( 2 x ) ( 0 - 0 )

W = 8 s e n 2 x 2 + 8 cos ⁡ 2 x 2 = 8

Pois s e n 2 x 2 + c o s 2 x 2 = 1

Além do W , vamos precisar saber da matriz W k , que vai ser a matriz W trocando a coluna k pela coluna

0 0 0 … 1

Por exemplo:

k = 1

W 1 = 0 y 2 y 3 … y n 0 y 2 ' y 3 ' … y n ' 0 y 2 ' ' y 3 ' ' … y n ' ' … … … … … 1 y 2 n - 1 y 3 n - 1 … y n n - 1

Mais um exemplo:

k = 2

W 2 = y 1 0 y 3 … y n y 1 ' 0 y 3 ' … y n ' y 1 ' ' 0 y 3 ' ' … y n ' ' … … … … … y 1 n - 1 1 y 3 n - 1 … y n n - 1

Outro exemplo:

k = n

W n = y 1 y 2 y 3 … 0 y 1 ' y 2 ' y 3 ' … 0 y 1 ' ' y 2 ' ' y 3 ' ' … 0 … … … … … y 1 n - 1 y 2 n - 1 y 3 n - 1 … 1

Deu pra entender né?!

No nosso caso,

Agora vamos calcular W 1 , W 2 e W 3 :

W 1 = 0 y 2 y 3 0 y 2 ' y 3 ' 1 y 2 " y 3 "

Substituindo temos:

W 1 = 0 cos ⁡ 2 x s e n 2 x 0 - 2 s e n 2 x 2 cos ⁡ 2 x 1 - 4 c o s ⁡ ( 2 x ) - 4 s e n ( 2 x )

W 1 = 0 - 2 s e n 2 x - 4 s e n 2 x - - 4 cos ⁡ 2 x 2 cos ⁡ 2 x - cos ⁡ 2 x 0 . - 4 s e n 2 x - ( 1 . 2 cos ⁡ 2 x + s e n 2 x 0 - 1 . - 2 s e n 2 x

W 1 = 0 + 2 cos ⁡ 2 x 2 + 2 s e n 2 x 2 = 2

Calculando W 2 :

W 2 = 1 0 s e n 2 x 0 0 2 cos ⁡ 2 x 0 1 - 4 s e n ( 2 x )

W 2 = 1 ( - 2 cos ⁡ 2 x ) - 0 0 - 0 + s e n 2 x 0 - 0 = - 2 c o s ⁡ ( 2 x )

Calculando W 3 :

W 3 = 1 cos ⁡ 2 x 0 0 - 2 s e n 2 x 0 0 - 4 c o s ⁡ ( 2 x ) 1

W 3 = 1 - 2 s e n ( 2 x ) - 0 - cos ⁡ 2 x 0 - 0 + 0 0 - 0 = - 2 s e n ( 2 x )

Tudo isso que estamos fazendo é pra chegar nos u k , com k = 1,2 , 3 … n , que vão ser encontrados resolvendo a seguinte conta

u k = ∫ f x det ⁡ W k det ⁡ W d x

Por exemplo, para k = 1 vamos ter:

u 1 = ∫ f x det ⁡ W 1 det ⁡ W d x

Para k = 2 :

u 2 = ∫ f x det ⁡ W 2 det ⁡ W d x

Para k = n :

u n = ∫ f x det ⁡ W n det ⁡ W d x

Esse k é importante para manter a ordem da parada e as contas darem certinho.

Sacou?! Agora vamos calcular u 1 ,   u 2 e u 3 do nosso caso:

Sendo f x = s e c ⁡ ( 2 x )

u 1 = ∫ f x det ⁡ W 1 det ⁡ W d x

Substituindo:

u 1 = ∫ sec ⁡ 2 x . 2 8 d x = 1 4 ∫ sec ⁡ 2 x d x = 1 4 ∫ sec ⁡ 2 x sec ⁡ 2 x + t g 2 x sec ⁡ 2 x + t g 2 x d x =

Eu apenas multipliquei em cima e em baixo por sec ⁡ 2 x + t g 2 x

Fazendo a distributiva temos:

u 1 = 1 4 ∫ sec 2 ⁡ 2 x + s e c ⁡ ( 2 x ) t g 2 x sec ⁡ 2 x + t g 2 x d x

Vamos fazer a substituição:

u = 2 x

d u = 2 d x

Fica:

u 1 = 1 8 ∫ sec 2 ⁡ u + s e c ⁡ ( u ) t g u sec ⁡ u + t g u d u

Vamos fazer a substituição:

v = sec ⁡ u + t g ( u )

d v = ( sec ⁡ u t g u + sec 2 ⁡ u ) d u    

Então:

u 1 = 1 8 ∫ sec 2 ⁡ u + sec ⁡ u t g u v d v sec 2 ⁡ u + sec ⁡ u t g u   =       1 8 ∫ d v v   = 1 8 ln ⁡ v      

u 1 = 1 8 ln ⁡ sec ⁡ u + t g u = 1 8 ln ⁡ sec ⁡ 2 x + t g 2 x  

Portanto:

u 1 = 1 8 ln ⁡ sec ⁡ 2 x + t g 2 x  

Vamos calcular u 2 :

u 2 = ∫ f x det ⁡ W 2 det ⁡ W d x

Substituindo:

u 2 = ∫ s e c ⁡ ( 2 x ) ( - 2 cos ⁡ 2 x ) 8 d x = 1 8 ∫ 1 . - 2 cos ⁡ 2 x c o s ⁡ ( 2 x ) d x = - 1 4 x

Portanto:

u 2 = - 1 4 x

Vamos calcular u 3 :

u 3 = ∫ f x det ⁡ W 3 det ⁡ W d x

Substituindo:

u 3 = ∫ s e c ⁡ ( 2 x ) ( - 2 sen ⁡ 2 x ) 8 d x = - 1   4 ∫ sen ⁡ 2 x c o s ⁡ ( 2 x ) d x

Vamos fazer a substituição:

u = 2 x

d u = 2 d x

Temos:

u 3 = - 1   8 ∫ sen ⁡ u c o s ⁡ ( u ) d u

Vamos fazer a substituição:

v = cos ⁡ u

d v = - s e n u d u

Fica:

u 3 = - 1   8 ∫ sen ⁡ u v d v - s e n ( u ) = 1 8 ∫ d v v = 1 8 ln ⁡ v = 1 8 l n ⁡ | c o s ⁡ ( u ) | = 1 8 l n ⁡ | c o s ⁡ ( 2 x ) |

Portanto:

u 3 = 1 8 l n ⁡ | c o s ⁡ ( 2 x ) |

Portanto a solução particular fica:

y p = u 1 . 1 + u 2 cos ⁡ 2 x + u 3 s e n ( 2 x )

y p = 1 8 ln ⁡ sec ⁡ 2 x + t g 2 x + - 1 4 x ( cos ⁡ 2 x )   + 1 8 l n ⁡ | c o s ⁡ ( 2 x ) | s e n ( 2 x )

Logo a solução geral será:

y = y h + y p

y = y h = c 1 + c 2 cos ⁡ 2 x + c 3 s e n 2 x + 1 8 ln ⁡ sec ⁡ 2 x + t g 2 x + - 1 4 x ( cos ⁡ 2 x )   + 1 8 l n ⁡ | c o s ⁡ ( 2 x ) | s e n ( 2 x )

Bora pros exercícios!!!

Relembrando y p :

Variação dos Parâmetros: Ordem Superior

Exercícios Resolvidos

Exercício Resolvido #1

Resolva pelo método da variação dos parâmetros a EDO a seguir e defina um intervalo no qual a solução geral seja válida:

y ' ' ' - 6 y " = x 2

Passo 1

Vamos começar encontrar a solução da EDO homogênea associada:

y ' ' ' - 6 y " = 0

A equação auxiliar é:

r 3 - 6 r 2 = 0

Passo 2

Agora vamos resolvê-la:

r 3 - 6 r 2 = 0

r 2 r - 6 = 0

As raízes são:

r 1 = 0

r 2 = 0

r 3 = 6

Passo 3

Agora que temos as raízes podemos encontrar a solução da homogênea:

y 1 = e 0 x = 1

y 2 = x e 0 x = x

y 3 = e 6 x

A solução da homogênea associada é:

y h = c 1 y 1 + c 2 y 2 + c 3 y 3

y h = c 1 + c 2 x + c 3 e 6 x

Passo 4

Agora vamos calcular a solução particular da EDO não homogênea:

y ' ' ' - 6 y " = x 2

O método da variação dos parâmetro supõe que :

y p = u 1 y 1 + u 2 y 2 + u 3 y 3

Substituindo:

y p = u 1 + u 2 x + u 3 e 6 x

Porém vimos que para encontrar u 1 . u 2 e u 3 precisamos calcular W ,   W 1 , W 2 e W 3

Passo 5

Vamos calcular W ,   W 1 , W 2 e W 3 :

W = y 1 y 2 y 3 y 1 ' y 2 ' y 3 ' y 1 " y 2 " y 3 "

Para calculá-lo, vamos calcular as derivadas necessárias:

De y 1 temos:

y 1 = 1

y 1 ' = 0

y 1 " = 0

De y 2 temos:

y 2 = x

y 2 ' = 1

y 2 " = 0

De y 3 temos:

y 3 = e 6 x

y 3 ' = 6 e 6 x

y 3 " = 36 e 6 x

Passo 6

Agora podemos encontrar:

W = y 1 y 2 y 3 y 1 ' y 2 ' y 3 ' y 1 " y 2 " y 3 "

Substituindo:

W = 1 x e 6 x 0 1 6 e 6 x 0 0 36 e 6 x

Como temos uma matriz triangular, o determinante é dado pelo produto dos elementos da diagonal principal:

det ⁡ W = 1 ⋅ 1 ⋅ 36 e 6 x = 36 e 6 x

Passo 7

Vamos calcular W 1 :

W 1 = 0 y 2 y 3 0 y 2 ' y 3 ' 1 y 2 " y 3 "

Substituindo:

W 1 = 0 x e 6 x 0 1 6 e 6 x 1 0 36 e 6 x

det ⁡ W 1 = 0 36 e 6 x - 0 - x 0 36 e 6 x - 6 e 6 x + e 6 x 0 - 1 = 6 x e 6 x - e 6 x = e 6 x ( 6 x - 1 )

Passo 8

Vamos calcular W 2 :

W 2 = y 1 0 y 3 y 1 ' 0 y 3 ' y 1 " 1 y 3 "

Substituindo:

W 2 = 1 0 e 6 x 0 0 6 e 6 x 0 1 36 e 6 x

det ⁡ W 2 = 1 ( 0 ) 36 e 6 x - 6 e 6 x   - 0 0 36 e 6 x - 0 6 e 6 x + e 6 x 0 - 0 = - 6 e 6 x

Passo 9

Vamos calcular W 3 :

W 3 = y 1 y 2 0 y 1 ' y 2 ' 0 y 1 " y 2 " 1

Substituindo:

W 3 = 1 x 0 0 1 0 0 0 1

det ⁡ W 3 = 1 ⋅ 1 ⋅ 1 = 1

Passo 10

Vamos calcular u 1 :

u 1 = ∫ f x det ⁡ W 1 det ⁡ W d x

Da EDO temos :

f x = x 2

Substituindo:

u 1 = ∫ x 2 e 6 x ( 6 x - 1 ) 36 e 6 x d x

Calculando a integral:

u 1 = 1 36 ∫ 6 x 3 - x 2 d x = 1 36 6 x 4 4 - x 3 3 = x 4 24 - x 3 108

Passo 11

Vamos calcular u 2 :

u 2 = ∫ f x det ⁡ W 2 det ⁡ W d x

Da EDO temos :

f x = x 2

Substituindo:

u 2 = ∫ x 2 ( - 6 e 6 x ) 36 e 6 x d x

Calculando a integral:

u 2 = - 1 6 ∫ x 2 d x = - 1 18 x 3

Passo 12

Vamos calcular u 3 :

u 3 = ∫ f x det ⁡ W 3 det ⁡ W d x

Da EDO temos :

f x = x 2

Substituindo:

u 3 = ∫ x 2 ( 1 ) 36 e 6 x d x

Calculando a integral:

u 3 = 1 36 ∫ x 2 e - 6 x d x

Vamos integrar por partes:

u = x 2

d v = e - 6 x d x

Encontramos:

d u = 2 x d x

v = - e - 6 x 6

Então temos:

∫ u d v = u v - ∫ v d u

u 3 = 1 36 ∫ x 2 e - 6 x d x = 1 36 - x 2 ( e - 6 x )   6 - ∫ - 2 x e - 6 x 6 d x

u 3 = - 1 216 x 2 e - 6 x - ∫ 2 x e - 6 x d x   ( i )

Vamos usar a integral por partes mais uma vez para calcular essa integral:

∫ 2 x e - 6 x d x

u = 2 x

d v = e - 6 x d x

Encontramos:

v = - e - 6 x 6

d u = 2 d x

A integral fica:

∫ 2 x e - 6 x d x = - 2 x e - 6 x 6 - ∫ - 2 e - 6 x 6 d x = - 1 3 x e - 6 x + e - 6 x 6

Substituindo me ( i ) :

u 3 = - 1 216 x 2 e - 6 x + 1 3 x e - 6 x + e - 6 x 6 = - e - 6 x 216 x 2 + x 3 + 1 18

Passo 13

Podemos então encontrar a solução particular:

y p = u 1 + u 2 x + u 3 e 6 x

Substituindo:

y p = x 4 24 - x 3 108 + x 5 24 - x 4 108 - 1 216 x 2 + x 3 + 1 18

Ajeitando:

y p = - 1 216 x 2 + x 3 + 1 18 - x 3 108 + x 4 1 24 - 1 108 + x 5 24

Passo 14

A solução geral vai ser:

y = y h + y p

y = c 1 + c 2 x + c 3 e 6 x - 1 216 x 2 + x 3 + 1 18 - x 3 108 + x 4 1 24 - 1 108 + x 5 24

Passo 15

Veja que a solução geral está definida para todo x   ∈ ( - ∞ , ∞ ) .

Resposta

y = c 1 + c 2 x + c 3 e 6 x - 1 216 x 2 + x 3 + 1 18 - x 3 108 + x 4 1 24 - 1 108 + x 5 24

x ∈ ( - ∞ , ∞ )

Exercício Resolvido #2

Use o método da variação dos parâmetros para determinar a solução geral da equação diferencial dada

y ' ' ' + y ' = tg ⁡ x

Com

0 < x < π

Passo 1

Bem aqui temos uma EDO que é não homogênea, para resolver essa EDO temos que achar a solução da EDO homogênea equivalente primeiro, depois vamos achar a solução particular.

A EDO homogênea equivalente vai ser

y ' ' ' + y ' = 0

Para resolvermos essa EDO temos que achar a equação característica dessa EDO

r 3 + r = 0

Passo 2

Vamos resolver essa equação agora

r 3 + r = 0

A primeira coisa que vamos fazer é colocar o r em evidência

r r 2 + 1 = 0

Assim vamos ter

r = 0         o u         r 2 + 1 = 0

E

r 1 = 0

r 2 = i

r 3 = - i

Passo 3

Agora que temos as raízes podemos ver como vão ser as soluções.

Repara que todas as raízes tem multiplicidade 1 , então as soluções vão ser

y 1 = C 1 e 0 x = C 1

y 2 + y 3 = e 0 x C 2 cos ⁡ x + C 3 sen ⁡ x = C 2 cos ⁡ x + C 3 sen ⁡ x

A nossa solução da homogênea vai ser

y h = y 1 + y 2 + y 3

y h = C 1 + C 2 cos ⁡ x + C 3 sen ⁡ x

Passo 4

Vamos agora para a solução particular dessa EDO.

y ' ' ' + y ' = tg ⁡ x

Com

0 < x < π

Ele disse para usarmos o método da variação de parâmetros para resolver isso. Mesmo que ele não tivesse dito como aquela função do lado direito é tg ⁡ x ela não se enquadra naquelas do coeficientes a determinar, então precisamos usar o método da variação de parâmetros mesmo.

A nossa solução particular vai ser

y p = y 1 u 1 + y 2 u 2 + y 3 u 3

Onde

y 1 = 1

y 2 = cos ⁡ x

y 3 = sen ⁡ x

Que são as funções que achamos da homogênea sem as constantes.

Precisamos achar u 1 , u 2 e u 3 , para isso temos que achar a matriz W que vimos na teoria, vai ser

W = y 1 y 2 y 3 y 1 ' y 2 ' y 3 ' y 1 ' ' y 2 ' ' y 3 ' ' = 1 cos ⁡ x sen ⁡ x 0 - sen ⁡ x cos ⁡ x 0 - cos ⁡ x - sen ⁡ x

Além disso devemos achar W 1 , W 2 e W 3

W 1 = 0 y 2 y 3 0 y 2 ' y 3 ' 1 y 2 ' ' y 3 ' ' = 0 cos ⁡ x sen ⁡ x 0 - sen ⁡ x cos ⁡ x 1 - cos ⁡ x - sen ⁡ x

W 2 = y 1 0 y 3 y 1 ' 0 y 3 ' y 1 ' ' 1 y 3 ' ' = 1 0 sen ⁡ x 0 0 cos ⁡ x 0 1 - sen ⁡ x

W 3 = y 1 y 2 0 y 1 ' y 2 ' 0 y 1 ' ' y 2 ' ' 1 = 1 cos ⁡ x 0 0 - sen ⁡ x 0 0 - cos ⁡ x 1

E nossas funções u vão ser

u k = ∫ tg ⁡ x det ⁡ W k det ⁡ W d x

Com k = 1,2 , 3

Passo 5

Vamos começar com u 1

u 1 = ∫ tg ⁡ x det ⁡ W 1 det ⁡ W d x

Vamos achar esses determinantes

det ⁡ W = 1 cos ⁡ x sen ⁡ x 0 - sen ⁡ x cos ⁡ x 0 - cos ⁡ x - sen ⁡ x = sen 2 ⁡ x + cos 2 ⁡ x = 1

det ⁡ W = 1

det ⁡ W 1 = 0 cos ⁡ x sen ⁡ x 0 - sen ⁡ x cos ⁡ x 1 - cos ⁡ x - sen ⁡ x = cos 2 ⁡ x + sen 2 ⁡ x = 1

det ⁡ W 1 = 1

Vamos ter

u 1 = ∫ tg ⁡ x d x = ln ⁡ sec ⁡ x

Passo 6

Vamos começar com u 2

u 2 = ∫ tg ⁡ x det ⁡ W 2 det ⁡ W d x

Vamos achar o det ⁡ W 2

det ⁡ W 2 = 1 0 sen ⁡ x 0 0 cos ⁡ x 0 1 - sen ⁡ x = - cos ⁡ x

det ⁡ W 2 = - cos ⁡ x

Vamos ter

u 2 = ∫ - tg ⁡ x cos ⁡ x d x = ∫ - sen ⁡ x d x = cos ⁡ x

Passo 7

Vamos começar com u 3

u 3 = ∫ tg ⁡ x det ⁡ W 3 det ⁡ W d x

Vamos achar o det ⁡ W 3

det ⁡ W 3 = 1 cos ⁡ x 0 0 - sen ⁡ x 0 0 - cos ⁡ x 1 = - sen ⁡ x

det ⁡ W 2 = - sen ⁡ x

Vamos ter

u 3 = - ∫ tg ⁡ x sen ⁡ x d x = - ∫ sen 2 ⁡ x cos ⁡ x d x = - ∫ 1 - cos 2 ⁡ x cos ⁡ x d x = - ∫ sec ⁡ x - cos ⁡ x d x

u 3 = ∫ - sec ⁡ x + cos ⁡ x d x = sen ⁡ x - ln ⁡ sec ⁡ x + tg ⁡ x

Passo 8

Assim vamos ter

u 1 = ln ⁡ sec ⁡ x

u 2 = cos ⁡ x

u 3 = sen ⁡ x - ln ⁡ sec ⁡ x + tg ⁡ x

E a nossa solução particular vai ser

y p = y 1 u 1 + y 2 u 2 + y 3 u 3

Onde

y 1 = 1

y 2 = cos ⁡ x

y 3 = sen ⁡ x

Vamos ter

y p = 1 . ln ⁡ sec ⁡ x + cos ⁡ x cos ⁡ x + sen ⁡ x sen ⁡ x - ln ⁡ sec ⁡ x + tg ⁡ x

y p = ln ⁡ sec ⁡ x + cos 2 ⁡ x + sen 2 ⁡ x - sen ⁡ x ln ⁡ sec ⁡ x + tg ⁡ x

y p = ln ⁡ sec ⁡ x - sen ⁡ x ln ⁡ sec ⁡ x + tg ⁡ x

Passo 9

A nossa solução vai ficar

y = y h + y p

y = C 1 + C 2 cos ⁡ x + C 3 sen ⁡ x + ln ⁡ sec ⁡ x - sen ⁡ x ln ⁡ sec ⁡ x + tg ⁡ x

Resposta

y = C 1 + C 2 cos ⁡ x + C 3 sen ⁡ x + ln ⁡ sec ⁡ x - sen ⁡ x ln ⁡ sec ⁡ x + tg ⁡ x

Exercício Resolvido #3

Resolva pelo método da variação dos parâmetros a EDO a seguir e defina um intervalo no qual a solução geral seja válida:

y ' ' ' - 2 y " - y ' + 2 y = e 3 x  

Passo 1

Vamos começar encontrar a solução da EDO homogênea associada:

y ' ' ' - 2 y " - y ' + 2 y = 0

A equação auxiliar é:

r 3 - 2 r 2 - r + 2 = 0

Passo 2

Agora vamos resolvê-la:

r 3 - 2 r 2 - r + 2 = 0

Colocando - r 2 em evidência:

- r 2 - r + 2 - r + 2 = 0

Colocando ( 2 - r ) em evidência:

2 - r - r 2 + 1 = 0

As raízes são:

r 1 = 2

r 2 = 1

r 3 = - 1

Passo 3

Agora que temos as raízes podemos encontrar a solução da homogênea:

y 1 = e 2 x

y 2 = e x

y 3 = e - x

A solução da homogênea associada é:

y h = c 1 y 1 + c 2 y 2 + c 3 y 3

y h = c 1 e 2 x + c 2 e x + c 3 e - x

Passo 4

Agora vamos calcular a solução particular da EDO não homogênea:

y ' ' ' - 2 y " - y ' + 2 y = e 2 x  

O método da variação dos parâmetro supõe que :

y p = u 1 y 1 + u 2 y 2 + u 3 y 3

Substituindo:

y p = u 1 e 2 x + u 2 e x + u 3 e - x

Porém vimos que para encontrar u 1 . u 2 e u 3 precisamos calcular W ,   W 1 , W 2 e W 3

Passo 5

Vamos calcular W ,   W 1 , W 2 e W 3 :

W = y 1 y 2 y 3 y 1 ' y 2 ' y 3 ' y 1 " y 2 " y 3 "

Para calculá-lo, vamos calcular as derivadas necessárias:

De y 1 temos:

y 1 = e 2 x

y 1 ' = 2 e 2 x

y 1 " = 4 e 2 x

De y 2 temos:

y 2 = e x

y 2 ' = e x

y 2 " = e x

De y 3 temos:

y 3 = e - x

y 3 ' = - e - x

y 3 " = e - x

Passo 6

Agora podemos encontrar:

W = y 1 y 2 y 3 y 1 ' y 2 ' y 3 ' y 1 " y 2 " y 3 "

Substituindo:

W = e 2 x e x e - x 2 e 2 x e x - e - x 4 e 2 x e x e - x

det ⁡ W = e 2 x 1 - ( - 1 ) - e x 2 e x + 4 e x + e - x 2 e 3 x - 4 e 3 x = 2 e 2 x - 6 e 2 x - 2 e 2 x = - 6 e 2 x

Passo 7

Vamos calcular W 1 :

W 1 = 0 y 2 y 3 0 y 2 ' y 3 ' 1 y 2 " y 3 "

Substituindo:

W 1 = 0 e x e - x 0 e x - e - x 1 e x e - x

det ⁡ W 1 = 0 1 - 1 - e x 0 + e - x + e - x - e x = - 1 - 1 = - 2

Passo 8

Vamos calcular W 2 :

W 2 = y 1 0 y 3 y 1 ' 0 y 3 ' y 1 " 1 y 3 "

Substituindo:

W 2 = e 2 x 0 e - x 2 e 2 x 0 - e - x 4 e 2 x 1 e - x

det ⁡ W 2 = e 2 x e - x   - 0 2 e x + 4 e x + e - x 2 e 2 x - 0 = e x + 2 e x = 3 e x

Passo 9

Vamos calcular W 3 :

W 3 = y 1 y 2 0 y 1 ' y 2 ' 0 y 1 " y 2 " 1

Substituindo:

W 3 = e 2 x e x 0 2 e 2 x e x 0 4 e 2 x e x 1

det ⁡ W 3 = e 2 x e x - 0   - e x 2 e 2 x - 0 + 0 2 e 3 x - 4 e 3 x = e 3 x - 2 e 3 x = - e 3 x

Passo 10

Vamos calcular u 1 :

u 1 = ∫ f x det ⁡ W 1 det ⁡ W d x

Da EDO temos :

f x = e 3 x

Substituindo:

u 1 = ∫ e 3 x ( - 2 ) - 6 e 2 x d x

Calculando a integral:

u 1 = 1 3 ∫ e x d x = e x 3

Passo 11

Vamos calcular u 2 :

u 2 = ∫ f x det ⁡ W 2 det ⁡ W d x

Da EDO temos :

f x = e 3 x

Substituindo:

u 2 = ∫ e 3 x 3 e x - 6 e 2 x d x

Calculando a integral:

u 2 = - 1 2 ∫ e 2 x d x = - 1 4 e 2 x

Passo 12

Vamos calcular u 3 :

u 3 = ∫ f x det ⁡ W 3 det ⁡ W d x

Da EDO temos :

f x = e 3 x

Substituindo:

u 3 = ∫ e 3 x - e 3 x - 6 e 2 x d x

Calculando a integral:

u 3 = 1 6 ∫ e 4 x d x

Encontramos:

u 3 = 1 6 ∫ e 4 x d x = 1 24 e 4 x

Passo 13

Podemos então encontrar a solução particular:

y p = u 1 e 2 x + u 2 e x + u 3 e - x

Substituindo:

y p = e x 3 e 2 x - 1 4 e 2 x e x + 1 24 e 4 x e - x

Ajeitando:

y p = e 3 x 3 - e 3 x 4 + e 3 x 24 = 3 e 3 x 24 = e 3 x 8

Passo 14

A solução geral vai ser:

y = y h + y p

y = c 1 e 2 x + c 2 e x + c 3 e - x + e 3 x 8

Passo 15

Veja que a solução geral está definida para todo x   ∈ ( - ∞ , ∞ ) .

Resposta

y = c 1 e 2 x + c 2 e x + c 3 e - x + e 3 x 8

x ∈ ( - ∞ , ∞ )

Exercícios de Livros Relacionados

14. Encontre uma fórmula envolvendo integrais para uma solução particular da equação diferencial y ' ' ' - y ' ' + y ' - y = g ( t )
Em cada um dos problemas de 9 a 12, encontre a solução do problema de valor inicial dado. Depois faça um gráfico da solução. 12. y ' ' ' - y ' = csc ⁡ t ; y π 2= 2 ,y ' α 2= 1 ,y ' ' α 2 = - 1
17. Encontre uma fórmula envolvendo integrais para uma solução particular da equação diferencial x 3 y ' ' ' - 3 x 2 y ' ' + 6 x y ' - 6 y = g t ,x0 Sugestão: Verifique se x , x 2 e x 3 são soluções d
15. Encontre uma fórmula envolvendo integrais para uma solução particular da equação diferencial y 4 - y = g ( t ) Sugestão: As funções sen ⁡ t , cos ⁡ t , senh ⁡ t e cosh ⁡ t formam um conjunto funda
16. Encontre uma fórmula envolvendo integrais para uma solução particular da equação diferencial y ' ' ' - 3 y ' ' + 3 y ' - y = g ( t ) Se g ( t ) = t - 2 e t , determine Y ( t ) .