Você está aqui: Matemática > Equações Diferenciais Ordinárias > Transformada de Laplace (3).
Formato Livro
08/04/2011
Equações Diferenciais Ordinárias
7: Transformada de Laplace (3)

5.4 Equações Diferenciais com Entradas Descontínuas

A transformada de Laplace é apropriada para o tratamento de uma equação diferencial com um termo não homogêneo, também denominado a função de entrada, descontínua, como já mencionado.


Exemplo 5.18. Considere a equação diferencial

y''+y'+54y=g(t),y(0)=0,y'(0)=0,

onde g(t)=1-uπ(t), uma função descontínua. Transformando a equação obtemos

s2Y(s)+sY(s)+54Y(s)=L{1}-L{uπ(t)}=1s(1-eπs),

de onde se pode concluir que

Y(s)=1-eπss(s2+s+54).

Vamos denotar H(s) a transformada de h(t), onde

H(s)=1s(s2+s+54),h(t)=L-1{H(s)}.

Então, com esta notação

Y(s)=(1-eπs)H(s)

e, pelo teorema 1 da seção anterior

y(t)=L-1{H(s)-eπsH(s)}=h(t)-uπ(t)h(t-π).

Resta então encontrar L-1{H(s)}. Para fazer isto escrevemos H em termos das frações parciais

H(s)=as+bs+cs2+s+54

o que pode ser conseguido se as constantes são a=4/5,b=-4/5   e   c=-4/5. Logo

H(s)=451s-45s+1s2+s+54.

Usando as transformadas conhecidas e suas propriedades encontramos

h(t)=45-45(cost+12sent) e-t/2

e, portanto,

y(t)={ 45-45(cost+12sent) e-t/2, t<π -45(1+eπ/2) (cost+12sent) e-t/2, tπ.

Exemplo 5.19 Veremos como utilizar a transformada acima para resolver uma equação diferencial que tem como função de entrada uma distribuição delta. Considere o problema

y''+y=4δ (t-2π),y (0)=1,y' (0)=0.

Este problema pode ser visto como a equação de movimento de uma partícula presa a uma mola com constante elástica unitária e em meio sem atrito, submetida a uma força de curta duração, como uma colisão, no momento t=2π. Ela parte afastada de uma unidade de distância da posição de equilíbrio no instante t=0 e tem velocidade inicial nula. Transformando toda a equação temos

s2Y(s)-s+Y(s)=4e-2πs

ou seja

Y(s)=ss2+1+4e-2πss2+1.

Utilizando os teoremas e transformadas inversas pertinentes temos

y(t)=cost+4sen(t-2π)u2π(t).

Como o seno é periódica com período T=2π temos sen(t-2π)=sen(t) e

y(t)={ cost, 0t2π, cost+4sent, t2π.


Figura 6: Exemplo 5.19

O gráfico desta solução mostra que a mola oscila com amplitude 1 até o instante t=2π. Depois disto passa a oscilar com maior amplitude, devido ao impulso sofrido.

5.4.1 Exercícios

Encontre a solução dos problemas de valor inicial

1. y''+y=f(t), y(0)=0,y'(0)=1onde f(t)={ 1, 0t<π/2 0, tπ/2

2. y''+2y'+2y=f(t), y(0)=0,y'(0)=1onde f(t)={ 1, πt<2π 0, 0t<πe t2π

3. y''+4y=sen t-u2π(t) sen (t-2π), y(0)=0,y'(0)=0

4. y''+4y=sen t+uπ(t)sen (t-π), y(0)=0,y'(0)=0

5. y''+2y'+y=f(t), y(0)=1,y'(0)=0onde f(t)={ 1, 0t<1 0, t1


Encontre a solução dos problemas de valor inicial involvendo deltas de Dirac:

6. y''+2y'+2y=δ(t-π), y(0)=1,y'(0)=0

7. y''+4y=δ(t-π)-δ(t-2π), y(0)=0,y'(0)=0

8. y''+2y'+y=δ(t)+u2π(t), y(0)=0,y'(0)=1

9. y''-y=2δ(t-1), y(0)=1,y'(0)=0

10. y''+2y'+3y=sen t+δ(t-π), y(0)=0,y'(0)=1

11. y''+ω2y=δ(t-π/ω), y(0)=1,y'(0)=1


Algumas Soluções:

1. y=1-cost+sen t-uπ/2(t)(1-sent)

3. y=16[1-u2π(t)](2sent-sen2t)

5. y=1-u1(t)[1-e-(t-1)-(t-1)e-(t-1)]

6. y=etcost+e-tsen t-uπ(t)e-(t-π)sen t

7. y=12uπ(t)sen 2t-12u2π(t)sen 2t

9. y=cosht+2u1(t)senh(t-1)

11. y=cosωt-ω-1uπ/ω(t)sen ωt

5.5 A Integral de Convolução

Pode ocorrer, na solução de uma equação diferencial, que a solução transformada Y (s) seja dada em termos do produto de duas transformadas conhecidas. Se F(s)=L {f(t) } e G(s)=L {g(t) } queremos saber qual é a transformada inversa do produto H (s)=F(s)G(s),ou seja L-1{F(s)G(s)}. O teorema da convolução nos fornece a transformada inversa deste produto.


Teorema 1: Seja f(t) e g(t) duas funções que possuem transformadas de Laplace, e F(s)=L{f(t)},G(s)=L{g(t)} as suas transformadas. Se H(s)=F(s)G(s) é o produto das transformadas então sua transformada inversa, h(t)=L-1{H(s)} é

h(t )=0tf (t'-τ)g (τ)dτ= 0tf (τ)g(t'-τ )dτ.

h(t) é a convolução de f e g e as integrais são chamadas de integrais de convolução.


Demonstração: por definição temos as expressões para F e G,

F(s)=0f(u)e-sudu,G(s)=0g(v)e-svdv,
H(s)=F(s)G(s)=0g(v)e-svdv 0f(u)e-sudu =0g(v)[0f(u)e-s(u+v)du]dv.

Para avaliar esta última integral introduzimos a variável u = t + v. Neste caso t = v quando u = 0 e du = dt o que permite escrever

H(s)=0 g(v) [ vf(t-v) e-stdt]  dv.

Admitindo que a ordem de integração pode ser invertida, integramos primeiro em v. Observe que os limites de integração, inicialmente v<t<, passam a ser 0<v<t (veja figura 6) e temos

H(s)=0e -st[0tf(t-v)g(v)dv]dt.

Figura 6: Teorema da convolução

Concluimos, portanto, que H(s) é a transformada da função h(t), dentro dos colchetes, ou seja H(s)=L{h(t)} onde

h(t)=0tf(t-v)g(v)dv.


é usada também a seguinte notação para a integral de convolução

h(t)=0tf(t-τ)g(τ)dτ=(fg)(t).

Exemplo 5.20: Vamos usar o teorema da convolução para encontrar a transformada inversa de

F(s)=1(s2+a2)2,a=constante.

Como conhecemos a transformada inversa

L-1 {1s2+a2 }=1asenat,

podemos encontrar

I=L-1 {1s2+a21 s2+a2 }=1a20t sen(at)sena (t-τ)dτ.

Para calcular esta integral podemos usar a relação trigonométrica

senAsenB= [cos(a-b)+cos(a+b)]/2

ou, no caso presente

senaτsena(t-τ)=12[cosa(2τ-t)-cosat],

para obter

I=12a2[0tcosa(2τ-t)dτ-cosat0tdτ]=12a3(senat-atcosat).

A primeira integral pode ser avaliada por meio de uma substituição de variáveis, u=a(2τ-t), du=2a dτ.

5.5.1 Exercícios

1. Determine as transformadas inversas:

a. F(s)=1s4(s2+1) b. F(s)=s(s+1)(s2+4)
c. F(s)=1(s+1)3(s2+4) d. F(s)=G(s)s2+1

Encontre a solução dos problemas de valor inicial

2. y''+ω2y=g(t), y(0)=0,y'(0)=1

3. y''+2y'+y=sen at, y(0)=0,y'(0)=0

4. 4y''+4y'+17y=g(t), y(0)=0,y'(0)=0

5. y''+y'+5/4y=1-uπ(t), y(0)=1,y'(0)=-3

6. y''+4y'+4y=g(t), y(0)=2,y'(0)=-3


Prove as seguintes propriedades da integral de convolução

7. fg=gf

8. f(g1+g2)=fg1+fg2

9. f(gh)=(fg)h


Algumas Soluções:

1a. f(t)=160t(t-τ)3sen τdτ

1b. f(t)=0te-(t-τ)cos2τdτ

1c. f(t)=120t(t-τ)e-(t-τ)sen 2τdτ

1d. f(t)=0tsen (t-τ)g(τ)dτ

2. y=1ωsen ωt+1ω0tsen ω(t-τ)g(τ)dτ

3. y=0te-(t-τ)sen (t-τ)sen aτdτ

4. y=180te-(t-τ)/2sen 2(t-τ)g(τ)dτ

5. y=e-t/2cost-12e-t/2sen t+0te-(t-τ)/2sen (t-τ)[1-uπ(τ)]dτ

5.6 Algumas Transformadas de Laplace e Propriedades

Listamos nesta seção um resumo das principais transformadas e as propriedades mais importantes.

Para n inteiro positivo:

L{1}=1s     L{coshat}=ss2-a2 L{t}=1s2     L{senhat}=as2-a2 L{tn}=n!sn+1     L{teat}=1(s-a)2 L{eat}=1s-a     L{tneat}=n!(s-a)n+1 L{senat}=as2+a2     L{ua(t)}=e-ass L{cosat}=ss2+a2     L{δ(t-a)}=e-sa


Denotando F(s)=L{f(t)}:

L{eatf(t)}=F(s-a) L-1{F(s-a)}=eatf(t)
L{ua(t)f(t-a)}=e-asF(s) L{e-asF(s)}=ua(t)f(t-a).

Se f é função periódica, f(t+T)=f(t) então

L{f(t)}=11-e-sT0Tf(t)e-stdt.


Transformada das derivadas:

L{f(t)}=sL{f(t)}-f(0),

L{f''(t)}=s2L{f(t)}-sf(0)-f'(0).


Integral de Convolução:

h(t)=L-1[F(s)G(s)]=0tf(t-τ)g(τ)dτ.


Leia no site: | Devemos Acreditar na Ciência? | Hipótese, Modelo e Teoria em Física | Cosmologia - Estrutura do Universo | História da Pessoa com Deficiência |