Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t =...
Transcript of Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t =...
![Page 1: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/1.jpg)
Capítulo 2
![Page 2: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/2.jpg)
Seção 2.1: Equações lineares; Fator integrante
• Uma EDO de primeira ordem é da forma
onde f é linear na variável y. Alguns exemplos típicos deste tipo de equações (com coeficientes constantes) saõ,
ou equações com coeficientes variáveis:
),( ytfdtdy
=
)()( tgytpdtdy
=+
bayy +−=′
![Page 3: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/3.jpg)
Equações com coeficientes constantes
• Para uma equação de 1ª Ordem,
precisamos de técnicas de integração elementares:
Cat ekkeaby
Ctaaby
dtaaby
dy
aaby
dtdy
±=+=
+−=−
−=−
−=−
∫∫
,/
/ln/
//
,bayy +−=′
![Page 4: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/4.jpg)
Equações com coeficientes variáveis: Fatores Integrantes
• Consideremos a seguintes equação:
• Esta técnica consiste em multiplicar a ambos lados da equação por uma função μ(t), escolhida adequadamente, de tal forma que o processo de integração seja simples.
)()( tgytpdtdy
=+
![Page 5: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/5.jpg)
Exemplo 1: Fator integrante (1 de 2)
• Considere a equação:
• Multiplicando ambos lados da equação por μ(t), obtemos
• Devemos escolher μ(t) de forma que o lado esquerdo da equação seja a derivada de uma expressão conhecida. Lembramos a regra da derivada de um produto:
• Devemos selecionar μ(t) tal que
2/2 teyy =+′
[ ] ydt
tddtdytyt
dtd )()()( μμμ +=
tettt 2)()(2)( =⇒=′ μμμ
)()(2)( 2/ teytdtdyt t μμμ =+
![Page 6: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/6.jpg)
Exemplo 1: Solução geral (2 de 2)
• Multipliamos μ(t) = e2t, a ambos lados da equação:
[ ]
tt
tt
tt
ttt
t
t
Ceey
Ceye
eyedtd
eyedtdye
etytdtdyt
eyy
22/
2/52
2/52
2/522
2/
2/
52
52
2
)()(2)(
2
−+=
+=
=
=+
=+
=+′
μμμ
![Page 7: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/7.jpg)
Fator integrante: Lado direito sendo uma função de t
• Consideremos o seguinte procedimento de solução:
[ ]
atatat
atat
atat
atatat
Cedttgeey
dttgeye
tgeyedtd
tgeyaedtdye
tgtytadtdyt
tgayy
−− +=
=
=
=+
=+
=+′
∫∫
)(
)(
)(
)(
)()()()(
)(
μμμ
![Page 8: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/8.jpg)
Exemplo 2: Solução geral (1 de 2)
• Consideremos a resolução de
usando a fórmula:
• Integrando por partes,
• Assim
tyy −=+′ 551
5/5/5/ )5()( tttatatat CedtteeCedttgeey −−−− +−=+= ∫∫
[ ]5/5/
5/5/5/
5/5/5/
550
5525
5)5(
tt
ttt
ttt
tee
dtetee
dttedtedtte
−=
−−=
−=−
∫∫∫∫
( ) 5/5/5/5/5/ 550550 ttttt CetCeteeey −−− +−=+−=
![Page 9: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/9.jpg)
Exemplo 2: Gráficos (2 de 2)
• O gráfico da esquerda mostra o campo de direções com várias curvas integrais.
• O gráfico da direita mostra várias soluções e uma particular (em vermelho) cujo gráfico contem o ponto (0,50).
5/550551 tCetytyy −+−=⇒−+−=′
![Page 10: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/10.jpg)
Exemplo 3: Solução geral (1 de 2)
• Podemos resolver
usando a fórmula:
• Integrando por partes,
• Então
tyy −=−′ 551
5/5/5/ )5()( tttatatat CedtteeCedttgeey +−=+= ∫∫ −−−
[ ]5/
5/5/5/
5/5/5/
5
5525
5)5(
t
ttt
ttt
te
dtetee
dttedtedtte
−
−−−
−−−
=
+−−−=
−=−
∫∫∫∫
[ ] 5/5/5/5/ 55 tttt CetCeteey +=+= −
![Page 11: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/11.jpg)
Exemplo 4: Gráficos (2 de 2)
• Curvas dadas pelo campo de direções.
• O gráfico a direita mostra várias curvas integrais, uma solução particular (em vermelho) cujo ponto inicial sobre o eixo y separa soluções que crescem positivamente e negativamente quando t→∞.
5/555/ tCetytyy +=⇒−=−′
![Page 12: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/12.jpg)
Método do fator integrante para uma equação linear de 1ª ordem.
• Seja
• Multiplicamos ambos os lados por μ(t),
• μ(t) deve ser tal que μ'(t) = p(t)μ(t). Se segue que
)()( tgytpy =+′
[ ] yttpdtdytyt
dtd )()()()( μμμ +=
)()()()()( ttgyttpdtdyt μμμ =+
![Page 13: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/13.jpg)
Método do fator integrante para uma equação linear de 1ª ordem.
• Escolhemos μ(t) tal que μ'(t) = p(t)μ(t). • Assumindo que μ(t) > 0, temos que
• Com k = 0, temos
assim μ(t) > 0 como desejado.
ktdtpttdtpttd
+=⇒= ∫∫∫ )()(ln)()()( μ
μμ
,)( )( tdtpet ∫=μ
![Page 14: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/14.jpg)
Solução Geral
• Resumindo:
• Logo
( )
( ) ( )
( ) ( ) ( ) ( ) ( ), onde ( ) p t dt
y p t y g tdyt p t t y t g t t edt
μ μ μ μ ∫
′ + =
+ = =
[ ]
( )
( ) ( ) ( )
( ) ( ) ( )
( ) ( ), ( )
( )p t dt
d t y t g tdt
t y t g t dt c
t g t dt cy onde t e
t
μ μ
μ μ
μμ
μ∫
=
= +
+= =
∫∫
![Page 15: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/15.jpg)
Exemplo 4: Solução geral (1 de 3)
• Resolver
levamos para a forma padrão:
• Logo
e então
( ) ,21,52 2 ==−′ ytyyt
2 5 , para 0y y t tt
′ − = ≠
222
2
2ln55
1
51
)(
)()(CtttCdt
tt
t
Ctdtt
t
Cdttgty +=⎥⎦
⎤⎢⎣⎡ +=
+=
+= ∫
∫∫μ
μ
2
1lnln22
)( 1)( 2
teeeet ttdt
tdttp===∫=∫=
⎟⎠⎞
⎜⎝⎛
−−μ
![Page 16: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/16.jpg)
Exemplo 4: Solução particular (2 de 3)
• Usando a condição inicial y(1) = 2 e a solução geral
obtemos
ou
22 2ln52)1( tttyCy +=⇒==
,ln5 22 Cttty +=
( )5/2ln5 2 += tty
![Page 17: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/17.jpg)
Exemplo 4: Gráficos (3 de 3)
• Os gráficos abaixo mostram curvas integrais e a solução particular que passa pelo ponto inicial (1,2).
( )2
2 2
2 2
IVP: 2 5 , 1 2
Solução geral: 5 ln
Solução particular: 5 ln 2
ty y t y
y t t Ct
y t t t
′ − = =
= +
= +
![Page 18: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/18.jpg)
Seção 2.2: Equações Separáveis
• Consideraremos uma sub‐classe de equações lineares e não‐lineares de primeira ordem
• Que podemos rescrever
• Por exemplo, seja M(x,y) = ‐ f (x,y) w N(x,y) = 1.
• Se M for apenas função de x e N apenas função de y
• Estas equações se chamam separáveis.
0),(),( =+dxdyyxNyxM
0),(),( =+ dyyxNdxyxM
0)()( =+ dyyNdxxM
),( yxfdxdy
=
![Page 19: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/19.jpg)
Exemplo 1:
• Resolver :
• Separando variáveis
• A solução está definida implicitamente..
11
2
2
−+
=yx
dxdy
( ) ( )( ) ( )
Cxxyy
Cxxyy
dxxdyy
dxxdyy
++=−
++=−
+=−
+=−
∫∫
3331
31
11
11
33
33
22
22
![Page 20: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/20.jpg)
Exemplo 2: Soluções Implícitas e Explícitas (1 de 4)
• Resolver:
• Separando variáveis
• A equação define a solução implicitamente. Podemos, neste caso explicitar a solução:
( )12243 2
−++
=y
xxdxdy
( ) ( )( ) ( )
Cxxxyy
dxxxdyy
dxxxdyy
+++=−
++=−
++=−
∫∫222
24312
24312
232
2
2
( ) ( )
Cxxxy
CxxxyCxxxyy
+++±=
++++±=⇒=+++−−
2212
224420222
23
23232
![Page 21: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/21.jpg)
Exemplo 2: PVI (2 de 4)
• Se a solução deve satisfazer y(0) = ‐1 podemos usar a forma implícita para determinar C
• Subsituindo C temos
• Usando a e forma explícita
• Substituindo
411
221 23
=⇒±=−
+++±=
CC
Cxxxy
3)1(2)1(222
2
232
=⇒=−−−
+++=−
CCCxxxyy
3222 232 +++=− xxxyy
4221 23 +++−= xxxy
![Page 22: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/22.jpg)
Exemplo 2 : condição Inicial y(0) = 3 (3 de 4)
• Com a condição inicial y(0) = 3, devemos escolher a raiz positiva:
4221 23 ++++= xxxy
![Page 23: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/23.jpg)
Exemplo 2: Domínio (4 de 4)
As soluções de
Está dada por
Para uma representação dada em forma explícita
e aqui o domínio (‐2, ∞). Note que x = ‐2 leva a y = 1, o que faz que o produz uma tangente vertical.
Reciprocamente, o domínio da solução pode ser determinada pelas tangentes verticais do gráfico.
2 3 2
3 2
2 2 2 3 (implícitamente)
1 2 2 4 (explícitamente)
y y x x x
y x x x
− = + + +
= − + + +
( ) ( ) ( )( )2212221 22 ++−=+++−= xxxxxy
( ) 1)0(,12
243 2
−=−++
= yy
xxdxdy
![Page 24: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/24.jpg)
Exemplo 3: a solução implícita para um PVI (1 de 2)
• Considere o seguinte PVI:
• Separando variáveis obtemos
• Usando a condição inicial
1)0(,31
cos3 =
+=′ y
yxyy
Cxyy
xdxdyyy
xdxdyy
y
+=+
=⎟⎟⎠
⎞⎜⎜⎝
⎛+
=+
∫∫sinln
cos31
cos31
3
2
3
1sinln 3 +=+ xyy
![Page 25: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/25.jpg)
• A modelagem de sistemas físicos são realizados com freqüência usando equações diferenciais.
• Construção do modelo: Transformar o problema real em termos matemáticos. Isto é feito considerando as leis físicas que regem o processo. Uma equação diferencial é um modelo matemático aproximado do problema original.
• Análise do modelo: Resolvendo as equações ou obtendo uma compreensão qualitativa da solução. Podemos simplificar o problema real, porém conservando as características do problema original.
• Comparação experimental: Verificar a solução ou sugerir reformulação matemática do problema.
Modelagem
![Page 26: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/26.jpg)
Exemplo 1: Problemas de misturas (1 de 7)
• Inicialmente (t = 0), um tanque contém Q0 lb de sal em 100 gal de agua. Supondo que uma solução salina entra no tanque com uma concentração de ¼ lb de sal/gal com uma vazão de r gal/min, e deixa o tanque com igual valor. (a) Escreva o PVI que descreve o processo.
(b) Determine a quantidade de sal Q(t) que permanece no tanque para qualquer valor de t.
(c) Determine a quantidade de sal no tanque depois de um tempo grande.
(d) Se r = 3 e Q0 = 2QL , determine o tempo T após do qual a sal dentro do tanque é igual a 2% de QL .
![Page 27: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/27.jpg)
Exemplo 1: (a) O PVI (2 de 7)
• Em t = 0, o tanque contém Q0 lb de sal dissolvidas em 100 gal de água. Uma solução salina ingressa ao tanque com concentração ¼ lb de sal/gal , com vazão r gal/min, (V=cte).
• A mistura é homogênea
• Entrada: (1/4 lb salt/gal)(r gal/min) = (r/4) lb/min
• Saída: Se há Q(t) lbs de sal no tanque no instante t, então a concentração de sal é Q(t) lb/100 gal, e o fluxo mássico de saída é [Q(t)r/100] lb/min.
• Assim
/ massa sal entrando - massa sal saindodQ dt =
0)0(,1004
QQrQrdtdQ
=−=
![Page 28: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/28.jpg)
Exemplo 1: (b) Achando a solução Q(t) (3 de 7)
• Para achar Q(t) resolvemos o PVI
• Usamos a técnica de fator integrante:
ou
0)0(,4100
QQrrQdtdQ
==+
[ ][ ] 100/
0
100/100/100/100/
100/
100/
2525)(
25254
)(
)(
rt
rtrtrtrt
rt
rtat
eQtQ
CeCeedtreetQ
eet
−
−−−
−+=
+=+=⎥⎦
⎤⎢⎣
⎡=
==
∫
μ
( ) 100/0
100/125)( rtrt eQetQ −− +−=
![Page 29: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/29.jpg)
Exemplo 1: (c) Achando a solução assintótica QL (4 of 7)
• Para achar QL de sal Q(t) :
• Este resultado faz sentido ?
• O gráfico mostra a curva integral para r = 3 e vários valores de Q0.
[ ]( ) lb 252525lim)(lim 100/0 =−+== −
∞→∞→
rt
ttL eQtQQ
( ) 100/0
100/125)( rtrt eQetQ −− +−=
![Page 30: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/30.jpg)
Exemplo 1: (d) Achar T (5 de 7)
• Supor r = 3 e Q0 = 2QL . Para achar o tempo T depois que Q(t) é 2% de QL , observamos que Q0 = 2QL = 50 lb, então
• Depois, 2% de 25 lb são 0.5 lb,
[ ] trt eeQtQ 03.100/0 25252525)( −− +=−+=
min 4.13003.0
)02.0ln(03.0)02.0ln(
02.025255.25
03.0
03.0
≈−
=
−==
+=−
−
T
Te
eT
T
![Page 31: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/31.jpg)
Exemplo 2: (e) Achando a vazão (6 de 7)
• Para achar a vazão r requerida se T não ultrapassa 45 minutes, temos di item (d) que Q0 = 2QL = 50 lb, com
e a curva solução decresce de from 50 para 25.5.
• Assim
100/2525)( rtetQ −+=
gal/min 69.845.0
)02.0ln(45.)02.0ln(
02.025255.25
45.0
10045
≈−
=
−==
+=−
−
r
re
er
r
![Page 32: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/32.jpg)
Exemplo 2: Polução de uma represa (1 de 7)
• Considere uma represa contendo inicialmente 10 milhões de litros de água limpa. Um córrego despeja água contaminada com uma vazão de 5 milhões l/ano, sendo que o nível da represa permanece constante. A concentração c(t) de água contaminada varia com o tempo periodicamente segundo:
c(t) = 2 + sen 2t g/l(a) Construa um modelo matemático para esse processo para determinar a quantidade de material poluente Q(t).(b) Grafique a solução e descreva o efeito de mudar a concentração de entrada.
![Page 33: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/33.jpg)
Exemplo 3: (a) PVI (2 de 7)
• Um balanço de massa de material poluente leva a
• Entrada: (2 + sen 2t g/l)(5 x 106 l/ano) • Saida: Se há Q(t) g de material poluente no instante t, a
concentração é Q(t) g/107 l, e sai com uma de:
[Q(t) g/107 l][5 x 106 l/ano]
/ fluxo mássico entrada - fluxo mássico saidadQ dt =
![Page 34: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/34.jpg)
Exemplo 3: (a) PVI (3 de 7)
– entrada: (2 + sin 2t g/l)(5 x 106 l/ano) – saida: [Q(t) g/107 l][5 x 106 l/ano] = Q(t)/2 g/ano
• Temos
• Mudando de variáveis (re‐escalamento) : Seja q(t) = Q(t)/106.
( )( )6 ( ) 2 sen 2 5 x 10 , (0) 02
dQ Q tt Qdt
= + − =
( ) 10 5 sen 2 , (0) 02
dq q t t qdt
+ = + =
![Page 35: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/35.jpg)
Exemplo 3: (a) Resolução (4 de 7)
• Para resolver o PVI
usamos o método do fator integrante:
• Integrando,
( )∫ +=
==− dtteetq
eettt
tat
2sin510)(
)(2/2/
2/μ
0)0(,2sin5102/ =+=+′ qtqq
2/
2/2/2/2/
173002sin
17102cos
174020)(
2sin17102cos
174020)(
t
tttt
etttq
Cteteeetq
−
−
−+−=
⎥⎦⎤
⎢⎣⎡ ++−=
![Page 36: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/36.jpg)
Exemplo 3: (a) Integração por partes (5 de 7)
( )
Ctetetdte
Ctetetdte
Ctetetdte
tdtetete
tdtetete
tdtetetdte
ttt
ttt
ttt
ttt
ttt
ttt
++−=
++−=
++−=
⎥⎦⎤
⎢⎣⎡ −+−=
⎥⎦
⎤⎢⎣
⎡⎟⎠⎞
⎜⎝⎛ −+−=
⎥⎦⎤
⎢⎣⎡ +−=
∫
∫
∫
∫
∫
∫∫
2sin17102cos
17402sin5
2sin1722cos
1782sin
2sin812cos
212sin
1617
2sin1612sin
812cos
21
2sin412sin
21
412cos
21
2cos412cos
212sin
2/2/2/
2/2/2/
2/2/2/
2/2/2/
2/2/2/
2/2/2/
![Page 37: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/37.jpg)
Exemplo 3: (b) Análise da solução (6 de 7)
• O nosso PVI e a nossa soloução é
• O gráfico e o campo de direções
• Note que o termo exponencial é
importante para valores pequeno
de t some para t grande.
y = 20 seria a solução de equilíbrio
se no fosse pelo termo sen(2t).
2/
173002sin
17102cos
174020)(
0)0(,2sin51021
tetttq
qtqdtdq
−−+−=
=+=+
![Page 38: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/38.jpg)
Seção 2.4: Equações lineares e não lineares
• Uma EDO de 1ª ordem da forma y' = f (t,y), é linear se a f for linear em y, e não linear se a f for não linear em y.
• Exemplos: y' = ty ‐ et, y' = ty2. • Veremos que equações lineares e não lineares diferenciam‐
se de várias formas:– Os teoremas sobre existência e unicidade das soluções são
diferentes. – As soluções das equações lineares podem ser expressas na forma de
uma solução geral, o que usualmente não acontece com as não lineares.
– As equações lineares possuem soluções dadas em forma explícita; em geral não temos soluções dadas implicitamente para as não lineares.
• Para ambos tipos de equações os métodos numéricos e gráficos são importantes.
![Page 39: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/39.jpg)
Teorema 2.4.1
• Considere o PVI de 1ª ordem linear :
Se as funções p e g são contínuas no intervalo aberto (α, β) que contém o ponto t = t0,
então existe uma única solução y = φ(t) que satisfaz as condições iniciais para cada t em (α, β).
• Sem prova: Use:
0)0(),()( yytgytpdtdy
==+
∫=+
=∫ t
tdssp
t
t ett
ydttgty 00
)(0)( where,
)(
)()(μ
μ
μ
![Page 40: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/40.jpg)
Teorema 2.4.2
• Considere o PVI de 1ª ordem não linear:
• Supor que f e ∂f/∂y são contínuas no retângulo aberto (t, y) ∈ (α, β) x (γ, δ) contendi o ponto (t0, y0). Então no intervalo (t0 ‐ h, t0 + h) ⊆ (α, β) existe uma única solução y = φ(t) que satisfaz o PVI.
• As condições do teorema são suficiente porém não necessárias. A continuidade de f garante existência porém não unicidade de φ.
0)0(),,( yyytfdtdy
==
![Page 41: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/41.jpg)
Exemplo 1: Linear PVI
• Consideremos o PVI:
• A solução está definida para t > 0, o intervalo para o qual p(t) = ‐2/t é contínua.
• Se a condição inicial for y(‐1) = 2, a solução está dada pela mesma expressão, porém definida para t < 0.
• Está garantido que a solução é
única no correspondente intervalo.
( ) 222 2ln521,52 tttyytyyt +=⇒==−′
![Page 42: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/42.jpg)
Exemplo 2: PVI não linear (1 e 2)
• Consideremos o PVI não linear:
• As funções f e ∂f/∂y estão dadas por
e são contínuas a exceção dos pontos sobre a reta y = 1.
• Então podemos construir um retângulo aberto em torno de (0, ‐1) no qual f e ∂f/∂y são contínuas, tão grande como possível sem cobrir y = 1.
• Qual a largura do retângulo? Para t > ‐2,
( ) ( ),
12243),(,
12243),( 2
22
−++
−=∂∂
−++
=y
xxyxyf
yxxyxf
( ) 1)0(,12
243 2
−=−++
= yy
xxdxdy
4221 23 +++−= xxxy
![Page 43: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/43.jpg)
Exemplo 2: Mudando a condição inicial (2 de 2)
• O nosso problema
sendo
que são ontínuas exceto na reta y = 1.
• Se mudamos a condição incial para y(0) = 1, o Teorema 2.4.2 não é satisfeito. Resolvendo temos
• Então a solução não é única.
( ) ( ),
12243),(,
12243),( 2
22
−++
−=∂∂
−++
=y
xxyxyf
yxxyxf
( ) 1)0(,12
243 2
−=−++
= yy
xxdxdy
0,221 23 >++±= xxxxy
![Page 44: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/44.jpg)
Exemplo 3: PVI Não linear
• Considere o seguinte PVI não linear
• As funções f e ∂f/∂y são
• Assim f é contínua, porém ∂f/∂y não é em y = 0, e o Teorema 2.4.2 não é satisfeito. Existem soluções não únicas. Separando variáveis e resolvendo, temos
• If initial condition is not on t‐axis, then Theorem 2.4.2 does guarantee existence and uniqueness.
3/23/1
31),(,),( −=
∂∂
= yytyfyytf
( )00)0(,3/1 ≥==′ tyyy
0,32
23 2/3
3/23/1 ≥⎟⎠⎞
⎜⎝⎛±=⇒+=⇒=− ttyctydtdyy
![Page 45: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/45.jpg)
Exemplo 4: PVI Não‐linear
• Consideremos o PVI não‐linear
• As funções f e ∂f/∂y estão dadas por
• f e ∂f/∂y são contínuas em t = 0, assim temos que a solução existe e é única.
• Separando variáveis obtemos
• A solução está definida em (‐∞, 1). Notamos que a singularidade em t = 1 não está presente na apresentação do problema.
yytyfyytf 2),(,),( 2 =∂∂
=
1)0(,2 ==′ yyy
ty
ctyctydtdyy
−=⇒
+−
=⇒+=−⇒= −−
11112
![Page 46: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/46.jpg)
Intervalo de definição: Equação Linear
• Temos que a solução de
existe em uma vizinhança de t = t0 onde p e g são contínuas.
• Assíntotas verticais ou outras descontinuidades na solução podem aparecer apenas nas descontinudades de p ou g.
• Porém as soluções podem ser diferenciáveis em pontos de descontinuidade de p ou g.
0)0(),()( yytgytpy ==+′
![Page 47: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/47.jpg)
Intervalos de definição: Equações não‐lineares
• Os intervalos de definição onde a solução existe podem ser difíceis de determinar.
• A solução y = φ(t) existe sempre que (t,φ(t)) permaneça dentro da região retangular indicada anteriormente. Isto é o que determina o valor de h. Desde que φ(t) não é conhecida, pode ser impossível determinar esta região.
• Neste caso, o intervalo no qual a solução existe pode não ter uma relação simples com a função f na equação diferencial
y' = f (t, y), em contraste com equações lineares.
• Além disso, quaisquer singularidades na solução dependerão da condição inicial e da própria equação.
• Comparemos esses comentários nos seguintes exemplos.
![Page 48: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/48.jpg)
Soluções gerais
• Para uma equação linear de primeira ordem, é possível obter uma solução contendo uma constante arbitrária, a partira da qual todas as soluções podem ser obtidas especificando essa constante.
• Para equações não‐lineares a solução geral pode não existir. Isto é, ainda com uma solução contendo uma constante arbitrária pode haver outras soluções que não podem ser obtidas especificando valores para esta constante.
• Consideremos o Exemplo 4: A função y = 0 é uma solução da equação diferencial, porém não pode ser obtida especificando um valor para c na solução achada usando separação de variáveis:
ctyy
dtdy
+−
=⇒=12
![Page 49: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/49.jpg)
Soluções explícitas: Equações lineares
• O PVI
existe em uma vizinhança de t0 na qual p e g são contínuas e é única.
• A solução possui uma representação explícita,
0 0( )0( ) ( )
, onde ( ) ,( )
t
t
t
p s dstt g t dt y
y t et
μμ
μ
+ ∫= =∫
0)0(),()( yytgytpy ==+′
![Page 50: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/50.jpg)
Aproximação da solução explícita
• A representação explícita da solução pode ser obtida sempre que as integrais envolvidas possam ser resolvidas.
• Se as integrais não podem ser resovlidas devemos usar métodos numéricos.
∑∫
∫
=
Δ≈
∫=+
=
n
kkkk
t
t
dssp
t
t
ttgtdttgt
ett
Cdttgty
t
t
1
)(
)()()()(
)( where,)(
)()(
0
00
μμ
μμ
μ
![Page 51: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/51.jpg)
Soluções implícitas para equações não lineares
• Para equações não‐lineares as representações explícitas podem não existir.
• As vezes é possível obter uma solução que defina a solução explicitamente.
• As vezes, são necessários fazer uma tabela y vs t. Como em um exemplo anterior
1sinln1)0(,31
cos 33 +=+⇒=
+=′ xyyy
yxyy
![Page 52: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/52.jpg)
Seção 2.5: Dinâmica Populacional
• Consideraremos equações do tipo y' = f (y), chamadas equações autônomas, onde a variável independente t não aparece explicitamente.
• O objetivo imediato é utilizar métodos geométricos que nos permitem obter informação qualitativa da solução sem precisar resolver a equação.
• Exemplo ( Crescimento exponencial):
• Solução:
0, >=′ rryy
rteyy 0=
![Page 53: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/53.jpg)
Crescimento logístico
• O modelo exponencial y' = ry, cuja solução é y = ert, predice um crescimento ilimitado com taxa r > 0 (que independe da população).
• Supondo agora que a taxa de crescimento depende to tamanho da população, podemos substituir r por uma função h(y) para obter y' = h(y)y.
• A função h(y) selecionada deve ser tal que– h(y) ≅ r > 0 quando y for pequeno, – h(y) decresça quando y crescer, e– h(y) < 0 quando y for suficientemente grande.
Obtemos assim h(y) = r – a y, sendo a > 0.
• Esta equação é chamada a equação logística.( ) 0,, >−=′ aryayry
![Page 54: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/54.jpg)
A equação Logística
• Que pode ser escrita como
onde K = r/a. A constante r é chamada taxa de crescimento intrínseca, e K representa a capacidade crítica da população.
• O campo de direções com
r = 1 e K = 10
,1 yKyr
dtdy
⎟⎠⎞
⎜⎝⎛ −=
( ) 0,, >−=′ aryayry
![Page 55: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/55.jpg)
A equação Logística: Soluções de equilíbrio
• A equação logística
• Temos dois estado de equilíbrio:
• No campo de direções, (r = 1, K = 10), observamos o comportamento das soluções próximas
das soluções de equilíbrio:
y = 0 é instável
y = 10 é assintóticamente estável
0,,1 >⎟⎠⎞
⎜⎝⎛ −= Kry
Kyr
dtdy
Ktyty ==== )(,0)( 21 φφ
![Page 56: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/56.jpg)
Soluções de equilíbrio
• As soluções de equilíbrio para y' = f (y) podem ser obtidas determinando as raízes de f (y) = 0.
• As raízes de f (y)= 0 são chamados pontos críticos.
• Os pontos crítico de
são y = 0 e y = K.
• Pontos críticos são funções
constantes (soluções de equilíbrio).
yKyr
dtdy
⎟⎠⎞
⎜⎝⎛ −= 1
![Page 57: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/57.jpg)
A equação logísitica: Análise qualitativa e gráfica (1 de 7)
• Começamos graficando f (y) vs. y.
• No caso de crescimento logístico, isto significa fazer o estudo da curva z=f (y) .
yKyryf ⎟⎠⎞
⎜⎝⎛ −= 1)(
![Page 58: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/58.jpg)
Pontos críticos (2 de 7)
• As intersecções de f com o eixo y acontecem em (0,0) e (K,0), correspondem aos pontos críticos.
• O vértice da parábola é (K/2, rK/4).
[ ]
4221
2
202
11)(
1)(
rKKKKrKf
KyKyKr
Kyy
Kryf
yKyryf
set
=⎟⎠⎞
⎜⎝⎛⎟⎠⎞
⎜⎝⎛ −=⎟
⎠⎞
⎜⎝⎛
=⇒=−−=
⎥⎦
⎤⎢⎣
⎡⎟⎠⎞
⎜⎝⎛ −+⎟
⎠⎞
⎜⎝⎛−=′
⎟⎠⎞
⎜⎝⎛ −=
![Page 59: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/59.jpg)
Crescendo e decrescendo (3 de 7)
• dy/dt > 0 para 0 < y < K, e então y é uma função crescente de t (indicado com uma seta para a direita ao longo do eixo y para 0 < y < K).
• Análogamente, y é decrescente para y > K (indicado com uma seta para a esquerda para y > K).
• O eixo y é chamado linha de fase.
0,1 >⎟⎠⎞
⎜⎝⎛ −= ry
Kyr
dtdy
![Page 60: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/60.jpg)
Comportamento assintótico (4 de 7)
• Observamos que dy/dt ≅ 0 quando y ≅ 0 e/ou y ≅ K, assim ysofre um achatamento próximo desses valores e y pula um degrau quando y se desloca desde 0 ou desde K.
yKyr
dtdy
⎟⎠⎞
⎜⎝⎛ −= 1
![Page 61: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/61.jpg)
Concavidade (5 de 7)
• Para estudar a concavidade y(t), calculamos y'':
• O gráfico possui concavidade para cima quando f e f' têm o mesmo sinal, o que acontece quando 0 < y < K/2 e y > K.
• y tem concavidade para abaixo quando f e f' têm sinal oposto, o que acontece quando K/2 < y < K.
• Os pontos de inflexão aparecem na intersecção da curva solução y e a reta y = K/2.
)()()()( 2
2
yfyfdtdyyf
dtydyf
dtdy ′=′=⇒=
![Page 62: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/62.jpg)
Análise das curvas solução (6 de 7)
• Temos:– y cresce quando 0 < y < K.
– y decresce quando y > K.
– Declividade de y aproximadamente nula quando y ≅ 0 ou y ≅ K.
– Concavidade para cima quando 0 < y < K/2 fr y > K.
– Concavidade para abaixo quando K/2 < y < K.
– Ponto de inflexão quando y = K/2.
![Page 63: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/63.jpg)
Resumo (7 de 7)
• A constante K é chamada capacidade crítica o nível de saturação para a espécie.
• Observamos que o comportamento da solução é bem diferente da solução exponencial que é determinado pelo efeito do termo não linear na equação logística.
![Page 64: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/64.jpg)
Resolvendo a equação log´sitica (1 de 3)
• Assumindo que y ≠ 0 e y ≠ K, podemos escrever:
• Expandindo em frações parciais,
• Assim temos
• Integrando
( ) rdtyKy
dy=
−1
rdtdyKy
Ky
=⎟⎟⎠
⎞⎜⎜⎝
⎛−
+/1
/11
CrtKyy +=−− 1lnln
( ) ( ) KyABKyBAyyB
KyA
yKy==⇒−+=⇒+
−=
−,111
111
![Page 65: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/65.jpg)
Resolvendo … (2 de 3)
• Temos:
• Se 0 < y0 < K, logo 0 < y < K e então
• ou
CrtKyy +=⎟⎠⎞
⎜⎝⎛ −− 1lnln
CrtKyy +=−− 1lnln
( )0
00 0
ln1 1 1
ou , (0)
rt C rt
rt
y y yrt C e cey K y K y K
y Ky onde y yy K y e
+
−
⎡ ⎤= + ⇔ = ⇔ =⎢ ⎥− − −⎣ ⎦
= =+ −
![Page 66: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/66.jpg)
Resolvendo …(3 de 3)
• Temos:
para 0 < y0 < K.
• Podemos mostrar que esta solução é também verdadeira para y0 > K.
• Finalmente a nossa solução é
( ) rteyKyKyy −−+
=00
0
( ) rteyKyKyy −−+
=00
0
![Page 67: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/67.jpg)
Comportamento assintótico
• Temos que
• Usamos limite para ver o comportamento assintótico:
• Concluímos que a solução de equilíbrio y(t) = K é assintóticamente estável, entanto que s solução y(t) = 0 é instável.
• A única forma de garantir que a solução permaneça próxima de zero é fazer y0 = 0.
( ) rteyKyKyy −−+
=00
0
( ) KyKy
eyKyKyy
trttt==
−+=
∞→−∞→∞→0
0
00
0 limlimlim
![Page 68: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/68.jpg)
Exemplo: Peixes (1 de 2)
• Seja y a biomassa (em kg) de população de peixes em t, com r = 0.71/ano e K = 80.5 x 106 kg. Se y0 = 0.25K, ache
(a) biomassa 2 anos depois
(b) o tempo τ de tal forma que y(τ) = 0.75K.
(a) Re‐escalando
então
e assim
( ) rteyKyKyy −−+
=00
0
( ) rteKyKyKy
Ky
−−+=
00
0
1
5797.075.025.025.0)2(
)2)(71.0( ≈+
= −eKy
kg 107.465797.0)2( 6×≈≈ Ky
![Page 69: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/69.jpg)
Exemplo: Peixes (2 de 2)
(b) Ache o tempo τ para o qual y(τ) = 0.75K.
( )
( )
( ) ( )
( )
0
0 0
0 0 0
0 0 0
0 0
0 0
0.751
0.75 1
0.75 0.75 10.25
0.75 1 3 1
1 0.25ln 3.095 anos0.71 3 0.75
r
r
r
r
y Ky K y K e
y y yeK K K
y K y K e y Ky K y Ke
y K y K
τ
τ
τ
τ
τ
−
−
−
−
=+ −
⎡ ⎤⎛ ⎞+ − =⎜ ⎟⎢ ⎥⎝ ⎠⎣ ⎦+ − =
= =− −
⎛ ⎞−= ≈⎜ ⎟⎜ ⎟
⎝ ⎠
( ) rteKyKyKy
Ky
−−+=
00
0
1
![Page 70: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/70.jpg)
Suporte crítico (1 de 2)
• Considere a equação logísitica modificada:
• O gráfico do lado direito f (y) é
0,1 >⎟⎠⎞
⎜⎝⎛ −−= ry
Tyr
dtdy
![Page 71: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/71.jpg)
Capacidade crítica: análise qualitativo e solução (2 de 2)
• Fazendo uma análise semelhante ao do caso logístico obtemos o gráfico da curva solução.
• T é um valor crítico para y0, no qual a população desaparece o cresce ilimitadamente, dependendo do lado em que fica y0em relação a T.
• Pode se mostrado que a solução para a equação crítica
é 0,1 >⎟⎠⎞
⎜⎝⎛ −−= ry
Tyr
dtdy
( ) rteyTyTyy
00
0
−+=
![Page 72: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/72.jpg)
Crescimento logístico com capacidade crítica (1 de 2)
• Para evitar crescimento ilimitado para y > T como antes, consideremos a equação logística modificada :
• O gráfico do lado direito f (y) é
1 1 , 0 e 0dy y yr y r T Kdt T K
⎛ ⎞⎛ ⎞= − − − > < <⎜ ⎟⎜ ⎟⎝ ⎠⎝ ⎠
![Page 73: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/73.jpg)
Crescimento logístico com capacidade crítica (2 de 2)
• T é a capacidade crítica para y0, em que a população some ou cresce para o valor K, dependendo onde está y0 em relação a T.
• K é o nível da capacidade crítica.
• Note que : y = 0 e y = K são soluções de equilíbrio estável,
e y = T é uma solução de equilíbrio instável.
![Page 74: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/74.jpg)
Seção 2.6: Equações exatas & Fatores Integrantes
• Consideremos uma EDO de primeira ordem da forma
• Supor que existe uma função ψ tal que
e tal que ψ(x,y) = c define y = φ(x) implicitamente. Logo
assim a equação original se converte em
• Então ψ(x,y) = c define uma solução implicitamente. • Neste caso, a EDO é exata.
0),(),( =′+ yyxNyxM
),(),(),,(),( yxNyxyxMyx yx == ψψ
[ ])(,),(),( xxdxd
dxdy
yxyyxNyxM φψψψ
=∂∂
+∂∂
=′+
[ ] 0)(, =xxdxd φψ
![Page 75: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/75.jpg)
• Suporque uma EDO pode ser escrita na forma
onde as funções M, N, My and Nx são todas contínuas em uma região retangular R: (x, y) ∈ (α, β ) x (γ, δ ). Então a Eq. (1) é uma equação diferencial exata sse
• Isto é, existe uma função ψ que satisfazendo as condições
sse M e N satisfazem a Eq. (2).
)1(0),(),( =′+ yyxNyxM
)2(),(),,(),( RyxyxNyxM xy ∈∀=
)3(),(),(),,(),( yxNyxyxMyx yx == ψψ
Teorema
![Page 76: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/76.jpg)
Exemplo 1: A equação exata (1 de 4)
Consideremos a EDO
sendo
temos
Usando o teorema anterior
Então
0)4()4(4
4=′−++⇔
−+
−= yyxyxyxyx
dxdy
yxyxNyxyxM −=+= 4),(,4),(
( , ) 4 ( , ) DO é exatay xM x y N x y E= = ⇒
yxyxyxyx yx −=+= 4),(,4),( ψψ
( ) )(4214),(),( 2 yCxyxdxyxdxyxyx x ++=+== ∫∫ψψ
![Page 77: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/77.jpg)
Exemplo 1: Solução (2 de 4)
Temos que
e
Se segue que
então
A solução está implicitamente por
yxyxyxyx yx −=+= 4),(,4),( ψψ
( ) )(4214),(),( 2 yCxyxdxyxdxyxyx x ++=+== ∫∫ψψ
kyyCyyCyCxyxyxy +−=⇒−=′⇒′+=−= 2
21)()()(44),(ψ
kyxyxyx +−+= 22
214
21),(ψ
cyxyx =−+ 22 8
![Page 78: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/78.jpg)
Exemplo 1: Campos de direção e curvas solução (3 de 4)
A EDO e as soluções estão dadas por
cyxyxyyxyxyxyx
dxdy
=−+⇒=′−++⇔−+
−= 22 80)4()4(4
4
![Page 79: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/79.jpg)
Exemplo 1: soluções explícitas e gráficos (4 de 4)
• A nossa solução está dada implicitamente por.
• Neste caso podemos explicitar y:cxxycxxyy +±=⇒=−−− 222 17408
cyxyx =−+ 22 8
![Page 80: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/80.jpg)
Exemplo 2: Solução exata (1 de 3)
Considere a EDO
Temos
e então
logo
0)1(sin)2cos( 2 =′−+++ yexxxexy yy
1sin),(,2cos),( 2 −+=+= yy exxyxNxexyyxM
( , ) cos 2 ( , ) é exatayy xM x y x xe N x y= + = ⇒
1sin),(,2cos),( 2 −+==+== yy
yx exxNyxxexyMyx ψψ
( ) )(sin2cos),(),( 2 yCexxydxxexydxyxyx yyx ++=+== ∫∫ψψ
![Page 81: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/81.jpg)
Exemplo 2: Solução (2 de 3)
Temos
Se segue que
logo
E a solução é
1sin),(,2cos),( 2 −+==+== yy
yx exxNyxxexyMyx ψψ
( ) )(sin2cos),(),( 2 yCexxydxxexydxyxyx yyx ++=+== ∫∫ψψ
kyyCyC
yCexxexxyx yyy
+−=⇒−=′⇒
′++=−+=
)(1)(
)(sin1sin),( 22ψ
kyexxyyx y +−+= 2sin),(ψ
cyexxy y =−+ 2sin
![Page 82: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/82.jpg)
Exemplo 2: Campos de direção e curvas solução (3 de 3)
A equação diferencial e a solução está dada por
cyexxyyexxxexy
y
yy
=−+
=′−+++2
2
sin,0)1(sin)2cos(
![Page 83: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/83.jpg)
Exemplo 3: A equação não exata (1 de 3)
Consideremos a EDO
Temos que
e então
• Para mostrar que o método não funciona escolhamos uma função ψ tal que
• Logo
0)2()3( 32 =′+++ yxxyyxy
32 2),(,3),( xxyyxNyxyyxM +=+=
2( , ) 3 2 2 3 ( , ) não é exatay xM x y x y y x N x y= + ≠ + = ⇒
32 2),(,3),( xxyNyxyxyMyx yx +==+== ψψ
( ) )(2/33),(),( 222 yCxyyxdxyxydxyxyx x ++=+== ∫∫ψψ
![Page 84: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/84.jpg)
Exemplo 3: A equação não exata (2 de 3)
Vemos que
e
• Então
• Assim, não existe uma tal função ψ. Entretanto, se procedemos (incorretamente) como antes obtemos
que não é a solução da EDO.
32 2),(,3),( xxyNyxyxyMyx yx +==+== ψψ
( ) )(2/33),(),( 222 yCxyyxdxyxydxyxyx x ++=+== ∫∫ψψ
kyxyxyCxxyC
yCxyxxxyyxy
+−=⇒−=′⇒
′++=+=
2/3)(2/3)(
)(22/32),(
23??
23?
23ψ
cxyyx =+ 23
![Page 85: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/85.jpg)
Exemplo 3: Gráficos (3 de 3)
A EDO e a solução dada implicitamente são
é solução ?
cxyyxyxxyyxy
=+
=′+++23
32 ,0)2()3(
![Page 86: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/86.jpg)
• As vezes é possível converte que não é não exata em uma exata multiplicando‐a por uma função, um fator integrante adequado μ(x,y):
• Para que esta equação seja exata devem ser satisfeitas
• Esta EDO pode ser difícil de resolver. Se μ for apenas função de x, então μy = 0 e assim podemos resolver
deixando o lado direito como uma função de x. Similarmente se μ for uma função de apenas de y.
Fatores integrantes
0),(),(),(),(0),(),(
=′+=′+
yyxNyxyxMyxyyxNyxM
μμ
( ) ( ) ( ) 0=−+−⇔= μμμμμ xyxyxy NMNMNM
,μμN
NMdxd xy −=
![Page 87: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/87.jpg)
Exemplo 4: A equação não exata
Consideremos a seguinte EDO.
com um fator integrante adequado resolvemos
• Multiplicando a equação diferencial por μ, obtemos uma equação exata
cuja solução está dada por
0)()3( 22 =′+++ yxyxyxy
xxxdx
dN
NMdxd xy =⇒=⇔
−= )(μμμμμ
,0)()3( 2322 =′+++ yyxxxyyx
cyxyx =+ 223
21
![Page 88: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/88.jpg)
Seção 2.7: Aproximação numérica; o método de Euler
• Seja o PVI
• Sabemos que se f e ∂f /∂y são continuas, então o PVI possui uma única solução y = φ(t) em alguma vizinhança de t0.
• Quando a equação diferencial for linear, separável ou exata podemos achar uma solução analítica.
• Entretanto, as soluções para a grande maioria de equações diferenciais não podem ser obtidas com estes métodos.
• Assim devemos encontrar formas alternativas para resolver este problemas.
00 )(),,( ytyytfdtdy
==
![Page 89: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/89.jpg)
Campos de direção
• Para o problema de valor inicial
podemos esboçar o campo de direções e visualizamos o comportamento das soluções. Este procedimento é relativamente simples e funciona ainda para equações complexas. Porém, não conduzem a soluções quantitativas as quais possam ser avaliadas ou comparadas.
,)(),,( 00 ytyytfy ==′
![Page 90: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/90.jpg)
Métodos numéricos
• Para o PVI
uma possibilidade consiste em calcular valores aproximadas da solução y = φ(t) em um conjunto de pontos t.
• Idealmente, a solução aproximada contém erros, que deveriam estar limitados para garantir um nível de exatidão.
• Para isso introduzimos o métodos das tangentes ou o método de Euler.
,)(),,( 00 ytyytfy ==′
![Page 91: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/91.jpg)
O Método de Euler
• Começamos por encontrar uma primeira aproximação de
para y = φ(t) em t0.
• A solução passa por (t0, y0) com coeficiente angular
f(t0,y0). A reta tangente no ponto inicial é
• A reta tangente é uma boa aproximação para a curva solução em um pequeno intervalo
• Se t1 estiver próximo de t0,
podemos aproximar φ(t1) por
( )( )0000 , ttytfyy −+=
,)(),,( 00 ytyytfy ==′
( )( )010001 , ttytfyy −+=
![Page 92: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/92.jpg)
A fórmula de Euler
• Para um ponto t2 perto de t1, aproximamos φ(t2) usando uma reta que passa por (t1, y1) com coeficiente angular f(t1,y1):
• Assim criamos uma sequência yn de aproximações para φ(tn):
onde fn = f(tn,yn).
• Para um passo h = tn – tn‐1, obtemos a fórmula de Euler
( )( )
( )nnnnn ttfyy
ttfyyttfyy
−⋅+=
−⋅+=−⋅+=
++ 11
12112
01001
( )( )121112 , ttytfyy −+=
…,2,1,0,1 =+=+ nhfyy nnn
![Page 93: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/93.jpg)
Aproximação de Euler
• Para graficar a solução aproximada de Euler, graficamos os pontos
(t0, y0), (t1, y1),…, (tn, yn), e então unimos os pontos com segmento s de reta.
( ) ( )1 1 , ,n n n n n n n ny y f t t onde f f t y+ += + ⋅ − =
![Page 94: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/94.jpg)
Exemplo 1: Método de Euler (1 de 3)
• Seja o PVI
usamos o método de Euler com h = 0.1 para obter a solução aproximada em t = 0.1, 0.2, 0.3, 0.4.
( )( )( )( )( )( )( )( )( ) 80.3)1.0(88.22.08.988.2
88.2)1.0(94.12.08.994.194.1)1.0(98.2.08.998.
98.)1.0(8.90
334
223
112
001
≈−+=⋅+=≈−+=⋅+=
≈−+=⋅+==+=⋅+=
hfyyhfyyhfyyhfyy
0)0(,2.08.9 =−=′ yyy
![Page 95: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/95.jpg)
Exemplo 1: Solução Exata (2 de 3)
• Podemos achar a solução exata :
( )
( )t
Ct
eyky
ekkey
Cty
dty
dyyy
yyy
2.0
2.0
149491)0(
,49
2.049ln
2.049
492.00)0(,2.08.9
−
−
−=⇒
−=⇒=±=+=
+−=−
−=−
−−=′=−=′
![Page 96: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/96.jpg)
Exemplo 1: Análise do erro (3 de 3)
• Analisando a tabela embaixo, observamos que os erros são pequenos. Os erros de arredondamento são pequenos e a solução exata é aproximadamente é linear em [0, 0.4].
t Exact y Approx y Error % Rel Error0.00 0 0.00 0.00 0.000.10 0.97 0.98 -0.01 -1.030.20 1.92 1.94 -0.02 -1.040.30 2.85 2.88 -0.03 -1.050.40 3.77 3.8 -0.03 -0.80
Erro relativo Porcentual 100exato aprox
exato
y yy−
= ×
![Page 97: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/97.jpg)
Exemplo 2: Método de Euler (1 de 3)
• Para o PVI
podemos usar o Método de Euler com h = 0.1 para aproximar a solução em t = 1, 2, 3, e 4.
• Solução exata:
( )( )( )( ) 15.4)1.0()15.3)(2(3.0415.3
15.3)1.0()31.2)(2(2.0431.231.2)1.0()6.1)(2(1.046.1
6.1)1.0()1)(2(041
334
223
112
001
≈+−+=⋅+=≈+−+=⋅+=
=+−+=⋅+==+−+=⋅+=
hfyyhfyyhfyyhfyy
1)0(,24 =+−=′ yyty
tety 2
411
21
47
++−=
![Page 98: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/98.jpg)
Exemplo 2: Análise do erro (2 de 3)
• As primeiras dez aproximações de Euler são mostradas na tabela embaixo.
• Os erro são pequenos no início, depois são catrastroficos, o que sugere uma solução não linear.
t Exato y Aprox y Erro % Rel Error0,00 1,00 1,00 0,00 0,000,10 1,66 1,60 0,06 3,550,20 2,45 2,31 0,14 5,810,30 3,41 3,15 0,26 7,590,40 4,57 4,15 0,42 9,140,50 5,98 5,34 0,63 10,580,60 7,68 6,76 0,92 11,960,70 9,75 8,45 1,30 13,310,80 12,27 10,47 1,80 14,640,90 15,34 12,89 2,45 15,961,00 19,07 15,78 3,29 17,27
t Exato y Aprox y Erro % Erro Rel0,00 1,00 1,00 0,00 0,001,00 19,07 15,78 3,29 17,272,00 149,39 104,68 44,72 29,933,00 1109,18 652,53 456,64 41,174,00 8197,88 4042,12 4155,76 50,69
2
Solução exata:7 1 114 2 4
ty t e= − + +
![Page 99: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/99.jpg)
Exemplo 2: Análise do erro & Graficos (3 de 3)
• Os gráficos mostram a solução exata (vermelho) e a aproximada de Euler (azul).
t Exata y Aprox y Erro % Erro rel0,00 1,00 1,00 0,00 0,001,00 19,07 15,78 3,29 17,272,00 149,39 104,68 44,72 29,933,00 1109,18 652,53 456,64 41,174,00 8197,88 4042,12 4155,76 50,69
2
Solução exata:7 1 114 2 4
ty t e= − + +
![Page 100: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/100.jpg)
Análise do erro (1 de 4)
• Se f e ∂f /∂y são continuas, então o PVI
possui uma solução y = φ(t) em uma vizinhança de t0.
• φ é o membro de uma família uni‐paramétrica de soluções que satisfaz φ(t0) = y0.
00 )(),,( ytyytfy ==′
![Page 101: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/101.jpg)
Análise do erro (2 de 4)
• A primeira etapa de Euler usa a reta tangente de φ no ponto (t0, y0) para aproximar φ(t1) que chamamos y1.
• O ponto (t1, y1) não pertence ao gráfico de φ, desde que y1é uma aproximação de φ(t1).
• Na iteração seguinte o método de Euler não usa a reta tangente de φ, senão uma curva próxima φ1 que passa pelo ponto (t1, y1).
• Assim o método de Euler usa
uma sequência de retas tangentes
de uma sequência de soluções
φ, φ1, φ2,… da equação diferencial
![Page 102: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/102.jpg)
Exemplo do análise do erro: Familia convergente (3 de 4)
• Desde que o Método de Euler usa as retas tangentes de uma sequência de diferentes soluções, a exatidão depende após vários passos do comportamento da solução em (tn, yn), n = 1, 2, 3, … . Consideremos o seguinte PVI
• Observemos o gráfico embaixo; não interessa quais soluções estamos aproximando com retas tangentes desde que todas as soluções ficam próximas quando t aumenta.
2/326)(1)0(,23 ttt eetyyyey −−− −−==⇒=−+=′ φ
![Page 103: Capítulo 2 - UFSCmtm.ufsc.br/~daniel/eqdifor/cap2.pdf · Exemplo 1: (a) O PVI (2 de 7) • Em t = 0, o tanque contém Q. 0. lb de sal dissolvidas em 100 gal de água. Uma solução](https://reader033.fdocumentos.tips/reader033/viewer/2022051916/6007faf82a0c2d6f1f6c68a5/html5/thumbnails/103.jpg)
Exemplo do análise do erro: Familia divergente (4 de 4)
• Consideremos agora:
• A família de curvas solução são divergentes. A cada passo do método de Euler as retas tangentes seguem uma solução diferente ao passo anterior, de tal forma que as aproximações se afastam cada vez mais na medida que t cresce.
4112471)0(,24 2tetyyyty ++−=⇒=+−=′