Enunciado

Considere o problema de valor inicial

y ' ' + 0.1 y ' + y = f ( t ) ,   y ( 0 ) = 0 ,   y ' ( 0 ) = 0

em que f(t) é a mesma que no Problema 19.

  1. Desenhe o gráfico da solução. Use um valor de n suficientemente grande e um intervalo de tempo longo o bastante para que a parte transiente da solução torne-se desprezível e o estado estacionário apareça claramente.
  2. Estime a amplitude e a frequência da parte correspondente ao estado estacionário da solução.
  3. Compare os resultados do item (b) com os da Seção 3.8 para um oscilador sob a ação de uma força senoidal.

Passo 1

Na letra a:

Como a EDO é linear, a transformação de Laplace pode ser aplicada para resolvê-lo. A transformada de Laplace de uma função y(t) é definida aqui como:

Y ( s )   =   ℒ { y ( t ) } = ∫ 0 ∞ e - s t y ( t )   d t

Consequentemente, a primeira e a segunda derivadas se transformam da seguinte maneira.

ℒ { d y d t } = s Y ( s )   -   y ( 0 )

ℒ { d ² y d t ² } = s ² Y ( s )   -   s y ( 0 ) - y ' ( 0 )

Aplique a conversão Laplace nos dois lados da EDO:

ℒ { y ' ' + 0.1 y ' + y } = ℒ { u 0 ( t )   + 2 ∑ k = 1 n ( - 1 ) k u k π ( t ) }

ℒ { y ' ' } + ℒ { 0.1 y ' } + ℒ { y } = ℒ { u 0 ( t ) }   + 2 ∑ k = 1 n ( - 1 ) k ℒ { u k π ( t ) }

[ s ² Y ( s )   -   s y ( 0 ) - y ' ( 0 ) ] + 0.1 [ s Y ( s ) - y ( 0 ) ] + Y ( s ) = ∫ 0 ∞ e - s t u 0 ( t ) d t   + ∑ k = 1 n ( - 1 ) k ∫ 0 ∞ e - s t u k π ( t ) d t

Passo 2

Aplique as condições iniciais, y ( 0 ) = 0 e y ' ( 0 ) = 0 .

( s 2 + 1 10 s + 1 ) Y ( s ) = 1 s + 2 ∑ k = 1 n ( - 1 ) k ( e - k π s s )

Como resultado da aplicação da transformada de Laplace, a EDO reduziu a uma equação algébrica para Y, (a solução transformada):

Y ( s ) = [ 1 s - s + 1 20 ( s + 1 20 ) 2 + 399 400 - 1 399 399 20 ( s + 1 20 ) 2 + 399 400 ]  

+ 2 ∑ k = 1 n ( - 1 ) k [ 1 s - s + 1 20 ( s + 1 20 ) 2 + 399 400 - 1 399 399 20 ( s + 1 20 ) 2 + 399 400 ] e - k π s

Passo 3

Agora, basta resolver esta equação e achar sua inversa!

y ( t ) = ℒ - 1 { Y ( s ) } =

= ℒ - 1 { [ 1 s - s + 1 20 ( s + 1 20 ) 2 + 399 400 - 1 399 399 20 ( s + 1 20 ) 2 + 399 400 ]   + 2 ∑ k = 1 n ( - 1 ) k [ 1 s - s + 1 20 ( s + 1 20 ) 2 + 399 400 - 1 399 399 20 ( s + 1 20 ) 2 + 399 400 ] e - k π s }

Chegando em:

y ( t ) =   1 - e - t 20 c o s ( 399 20 t ) - 1 399 e - t / 20 s e n ( 399 20 t )

+ 2 ∑ k = 1 n ( - 1 ) k { 1 - e - ( t - k π ) / 20 c o s [ 399 20 ( t - k π ) ] - 1 399 e - ( t - k π ) 20 s e n [ 399 20 ( t - k π ) ] } H ( t - k π )

Passo 4

Abaixo estão vários gráficos de y (t) versus t para vários valores de n.

Passo 5

Na letra b:

A amplitude da parte estável da solução é de cerca de 12,6 e o período é T = 2 π , o que significa que a frequência angular é ω = 2 π / T = 1 . O problema correspondente com um termo de força senoidal foi considerado na Seção 3.8. Observe que o custo se assemelha mais a f(t).

z ' ' + 0.1 z ' + z = c o s ( t ) ,             z ( 0 ) = 0 ,   z ' ( 0 ) = 0

Passo 6

A solução para esta EDO está plotada abaixo

A amplitude da parte estável da solução aqui é de cerca de 9,9, e o período é T = 2 π o que significa que a frequência angular é ω = 2 π / T = 1 . Essa diferença de amplitude se deve ao fato de que a força do cosseno varia em força com o tempo, enquanto f(t) fornece pulsos com força máxima.

Resposta

  1. y ( t ) =   1 - e - t 20 c o s ( 399 20 t ) - 1 399 e - t / 20 s e n ( 399 20 t ) + 2 ∑ k = 1 n ( - 1 ) k { 1 - e - ( t - k π ) / 20 c o s [ 399 20 ( t - k π ) ] - 1 399 e - ( t - k π ) 20 s e n [ 399 20 ( t - k π ) ] } H ( t - k π )
  2. z ' ' + 0.1 z ' + z = c o s ( t ) ,             z ( 0 ) = 0 ,   z ' ( 0 ) = 0
  3. A amplitude da parte estável da solução aqui é de cerca de 9,9, e o período é T = 2 π o que significa que a frequência angular é ω = 2 π / T = 1 . Essa diferença de amplitude se deve ao fato de que a força do cosseno varia em força com o tempo, enquanto f(t) fornece pulsos com força máxima.
Ver Outros Exercícios desse livro

Exercícios de Livros Relacionados

Encontre a solução do problema de valor inicial dado.Desenhe os gráficos da solução e do termo não homogêneo; explique a relação entre eles. y ' '+ y=f ( t ) ;y ( 0 )=0 ,y ' ( 0 ) = 1 ; f ( t )=1 ,0 ≤
Encontre a solução do problema de valor inicial dado.Desenhe os gráficos da solução e do termo não homogêneo; explique a relação entre eles. y ' '+y ' + 5 4 y = t - u π / 2 ( t )( t - π 2 ) ;y ( 0 )=0
Encontre uma expressão envolvendo u c ( t ) para uma função f cujo gráfico é uma rampacrescente de zero em t = t 0 até o valor h em t = t 0 + k
Encontre a solução do problema de valor inicial dado.Desenhe os gráficos da solução e do termo não homogêneo; explique a relação entre eles. y ' ' + 3 y ' + 2 y = f ( t ) ;y ( 0 ) = 0 ,y ' ( 0 ) = 0 ;
Encontre a solução do problema de valor inicial dado.Desenhe os gráficos da solução e do termo não homogêneo; explique a relação entre eles. y ' ' + 3 y ' + 2 y = u 2 ( t ) ;y ( 0 ) = 0 ,y ' ( 0 ) = 1