Docsity
Docsity

Prepare-se para as provas
Prepare-se para as provas

Estude fácil! Tem muito documento disponível na Docsity


Ganhe pontos para baixar
Ganhe pontos para baixar

Ganhe pontos ajudando outros esrudantes ou compre um plano Premium


Guias e Dicas
Guias e Dicas

Equações Diferenciais, Notas de estudo de Engenharia Elétrica

Apostila do Professor Altemir sobre equações diferenciais, EDO, EDP

Tipologia: Notas de estudo

2010
Em oferta
30 Pontos
Discount

Oferta por tempo limitado


Compartilhado em 26/02/2010

ana-burlan-11
ana-burlan-11 🇧🇷

1 documento

Pré-visualização parcial do texto

Baixe Equações Diferenciais e outras Notas de estudo em PDF para Engenharia Elétrica, somente na Docsity! 1 Ministério da Educação Universidade Tecnológica Federal do Paraná - Campus Curitiba Departamento Acadêmico de Matemática EQUAÇÕES DIFERENCIAIS Notas de aula Professor: Altemir José Borges Curitiba Agosto de 2006 2 EQUAÇÕES DIFERENCIAIS Definição: Chama-se equação diferencial à equação que possui as derivadas ou diferenciais de uma ou mais variáveis dependentes, em relação a uma ou mais variáveis livres. Exemplos: a) 13 −= x dx dy b) xey dx dy dx yd 5 2 2 6127 =+− c) x dx dy dx yd cos5 43 2 2 =⎟ ⎠ ⎞ ⎜ ⎝ ⎛−⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ d) xyz y zx x z 3= ∂ ∂ − ∂ ∂ Classificação: A equação será chamada de ordinária se as variáveis dependentes forem função de uma única variável livre, caso contrário, serão chamadas de equações diferenciais parciais. As equações dos exemplos a, b e c anteriores são equações diferenciais ordinárias e a equação do exemplo d é uma equação diferencial parcial. Ordem: Chama-se ordem de uma equação diferencial à ordem da derivada de maior ordem. As equações a) e d) são de primeira ordem, já os exemplos b) e c) são de segunda ordem. Grau: Grau é o maior expoente da derivada de maior ordem. As equações a, b e d são de primeiro grau e o exemplo c é do terceiro grau. Solução: É uma função que quando substituída na equação diferencial a transforma numa identidade. As soluções podem ser: solução geral, particular ou singular. Chama-se solução geral à família de curvas integrais que verifica a equação diferencial e possui constantes arbitrárias. Chama-se solução particular de uma equação diferencial à solução obtida a partir da solução geral impondo condições iniciais ou de contorno. Geralmente as condições iniciais serão dadas para o instante inicial, já as condições de contorno aparecem quando nas equações de ordem superior os valores da função e de suas derivadas são dadas em pontos distintos. Por exemplo: Resolver a equação diferencial ordinária (EDO) xyy 6'''5 −=+ , sujeita às condições iniciais y(0) =2 e y’(0) = 3, ou resolver a EDO xyy 6'''5 −=+ , sujeita às condições de contorno y(0)=2 e y’(1)=3. Chama-se solução singular de uma equação diferencial à envoltória1 da família de curvas integrais. Teorema da existência: A equação ),( yxg dx dy = admite solução se: • g(x,y) é contínua e unívoca em uma região D de pontos (x,y). • yg ∂∂ existe e é contínua em todos os pontos de D. ____________________ 1 Envoltória de uma família de curvas é a uma curva tangente a todas as curvas da família. 5 5. 01)1( 222 =−−− dyxdxyx 6. 0)1( =−− ydxdyx 7. 2 2 1 1 x y dx dy + + = 8. x dx dy 5sen= 9. 03 =+ dyedx x 10. 6)1( +=+ x dx dyx 11. yxy 4'= 12. 2 3 x y dx dy = 13. x yx dy dx + = 1 22 14. yxe dx dy 23 += 15. 0)2()4( 22 =+−+ dxxyxdyyxy 16. xdxdyxy =+ )1(2 17. 21ln ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ += x y dy dxxy 18. dyxxdxe y )cos1(sen)1( +=+− , com y(0)=0 19. dxyxydy 2/12 )1(4 += , com y(0)=1 20. )1(4 2 += x dy dx , com 1)4( = πx 21. xyyyx −='2 , com y(-1)=-1 22. 2)( y dx dyee xx =+ − 23. 2pp dt dp −= 24. xyyx1 dx dy +++= 25. 2y(0) com ,2y2xxy dx dy =−−+= 26. 0sen)1(cos =++ − ydyeydx x , com 4)0( π=y 27. 1 10 2 + = xdx dy , com y(0)=0 28. )cos( yx dx dy += (Dica: Faça x+y=t) 29. 2)1(' ++= yxy (Dica observe o ex. 28) 6 30. )(' 2 yxtgy += (Dica observe o ex. 28) 31. 322' +−+= xyy (Dica observe o ex. 28) 32. Encontre as soluções singulares da equação dydxyx =− 21 RESPOSTAS 1. Cyxx =−− 2 3 2 2. C y x = 3. Cyln3)x4(2x424 3 =−−+−− 4. Cyx =+− coscos 5. Cy x x =−+ arcsen1 6. )1( −= xCy 7. Cx Cxy − + = 1 8. Cxy +−= 5cos 5 1 9. Cey x += −3 3 1 10. Cxxy +++= )1ln(5 11. 4Cxy = 12. Cxy += −− 12 2 13. Cxxyxx +=+− 3)ln(33 14. Cee xy +=− − 32 23 15. )4(2 22 xCy +=+ 16. Cxxy ++−= )1ln(2 17. Cyyyxxx +++=− ln2 29 ln 3 233 18. ( )( ) 41cos1 =++ yex 19. 221 22 +=+ xy 20. ( )434 π−= ytgx 21. ( )xexy 11+−= 22. Cetgy x +=− −− )(11 23. tCe p p = −1 24. C 2 xx)y1ln( 2 ++++ 25. 1e3y x22 2x −= − 26. 22sec)1( =+ ye x 27. arctgxy 10= 28. Cxyxyxec +=+−+ )cot()(cos 29. )(1 Cxtgxy ++−−= 30. Cyxxy =++− )(2sen22 31. 2)()32(4 Cxxy +=+− 32. y=1 ou y=-1 2° TIPO: EQUAÇÕES HOMOGÊNEAS Definição: A função definida por z=f(x,y) será uma função homogênea de grau m se tivermos f(λx,λy)= λmf(x,y). Exemplos: a) f(x,y)=2x3+5xy2 é homogênea de grau 3, pois f(λx,λy)=2(λx)3+5λx.(λy)2=λ3f(x,y). b) f(x,y)=yex/y é homogênea de grau 1, pois f(λx,λy)=λyeλx/λy=λf(x,y). Definição: A equação M(x,y)dx+N(x,y)dy=0 será chamada de equação diferencial homogênea se M e N forem funções homogêneas de mesmo grau. Resolução: 7 Se Mdx + Ndy = 0 for uma equação diferencial homogênea, então ela poderá ser escrita da forma ⎟ ⎠ ⎞ ⎜ ⎝ ⎛= x yf dx dy , onde a mudança de variáveis x yt = irá separar as variáveis. Exemplo: Determine a solução de 06)32( 22 =−− xydydxyx , sujeita à condição inicial y(1)=1/3. Como as funções M(x,y)=2x2-3y2 e N(x,y)=-6xy são funções homogêneas de grau 2, então a equação dada é homogênea. Fazendo x yt = , ou y=x.t (1) e diferenciando, teremos dy=x.dt+t.dx (2). Substituindo (1) e (2) na equação dada vem: 0)...(.6))(32( 22 =+−− dtxdxtxtxdxxtx 0)..(.6)32( 222 =+−− dtxdxttxdxtx 0..6)632( 22 =−−− dttxdxtt 0..6)92( 2 =−− dttxdxt Separando as variáveis, resulta: 0 92 .6 2 =− − t dtt x dx . Integrando teremos Ctx =−+ )92ln(ln3 2 Eliminando os logaritmos Ctx =− )92.( 23 Voltando para as variáveis x e y: C x yx = ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎣ ⎡ ⎟ ⎠ ⎞ ⎜ ⎝ ⎛− 2 3 92. Cxyx =− 23 92 Impondo a condição inicial y(1)=1/3, teremos a solução particular: 192 23 =− xyx Resolva as seguintes equações: 1) 0xydy2dx)yx( 22 =−− 2) 0dy)y4x(dx)yx2( =+−− 3) 0xydydx)yx( 22 =−+ 4) 0xydy2dx)y3x( 22 =+− , com y=1 e x=2 5) 0xdydx)yx( =+− 6) 0dy)x2y(xdx =−+ 7) 0dyxdx)yxy( 22 =−+ 8) xy xy dx dy + − = 9) 0dy)xyx(ydx =++− 10) dy)yx3(ydxx2 332 += 11) y x x y dx dy += 12) y x2 ye4x dy dxy − += 10 30) ( ) ( ) 0dy5y6xdx3y4x =−−−−− 31) ( ) ( ) 0dy1y3x9dx2yx3 =+−++− RESPOSTAS 1. Cxy3x 23 =− 2. Cy4xy2x2 22 =−− 3. 2x2y2 Cex = 4. 8 x31 x y −= 5. Cxyxx =+ln 6. )()ln()( yxCyyxyx −+=−− 7. Cyxyx =+ ln 8. Cx ytgyx =++ − )(2)ln( 122 9. 2)(ln4 Cyyx −= 10. 2339 )( yxCy += 11. Cyx y += ln2)( 2 12. Cye yx += ln82 13. Cx yx =)cos( 14. xyeCxxy 2=+ 15. 333 x8xlnx3y =+ 16. 2 2/12/3 yyxx =+ 17. 1 -exln x y = 18. 0x-yxlnxx ylnx4 =++ 19. 2 3 2 32/12 3 523ln3 xyyxxx =++ 20. 0xyln)yx( =++ 21. 2)y x--2(1lny 2 1 += 22. Cy4x2yxy6x2 22 =+−−− 23. Cy2yxy6x2 22 =++− 24. )3yx(C)1yx( 3 −+=−− 25. C)3y5x10ln(4y15x5 =−−+− 26. C)7y3x2ln(9y3x3 +−+−=+ 27. C)3y3x3ln(2yx3 =+−−++ 28. 3)y-(1)-y(x 3 +=+ xC 29. C4x-8y5)8yln(4x =+++ 30. 2)-3y-(1)-2y-(x 2 xC= 31. 1)2y-ln(6xC6y2x +−=++ 4° TIPO: EQUAÇÕES EXATAS Forma : A equação Mdx+Ndy=0 será uma equação diferencial exata , quando existir uma função f(x,y)=C tal que df=Mdx+Ndy = 0 ou se a relação x N y M ∂ ∂ = ∂ ∂ for verdadeira. Resolução: Dada a equação diferencial exata Mdx+Ndy=0 (1) e seja z=f(x,y)=C sua solução, cuja diferencial dada por dy y fdx x fdz ∂ ∂ + ∂ ∂ = (2). Então, comparando (1) e (2) teremos: ),( yxM x f = ∂ ∂ (3) e ),( yxN y f = ∂ ∂ (4). Para obtermos a sua solução z=f(x,y) deveremos integrar, por exemplo,a expressão (3), em relação à variável x, da qual teremos ∫ += )(),(),( ygdxyxMyxf (5). Derivando parcialmente (5) em relação à y teremos: )(' ),( yg y dxyxM y f + ∂ ∂ = ∂ ∂ ∫ (6). 11 Igualando (6) e (4) resulta: ),()(' ),( yxNyg y dxyxM =+ ∂ ∂∫ . Isolando g’(y) e integrando em relação a y acharemos 1 ),( ),()( Cdy y dxyxM yxNyg +⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎝ ⎛ ∂ ∂ −= ∫ ∫ (7). Substituindo (7) em (5) teremos a solução geral da equação exata, que é Cdy y dxyxM yxNdxyxMyxf =⎟ ⎟ ⎠ ⎞ ⎜ ⎜ ⎝ ⎛ ∂ ∂ −+= ∫ ∫∫ ),( ),(),(),( . Exemplo: Resolver a seguinte equação diferencial 0)542()23( 2 =+−++ dyyxdxyx . Inicialmente vamos verificar a que modelo esta equação pertence. i. Ela não é de variáveis separáveis porque temos soma das variáveis x e y, ii. Ela não é homogênea porque os coeficientes das diferenciais não são funções homogêneas, iii. Para verificarmos se a equação é exata vamos utilizar a relação x N y M ∂ ∂ = ∂ ∂ . 2)23( 2 = ∂ +∂ = ∂ ∂ y yx y M 2)542( = ∂ +−∂ = ∂ ∂ x yx x N Como a condição x N y M ∂ ∂ = ∂ ∂ é verificada temos que a equação é exata. A solução f(x,y)=C verifica 0= ∂ ∂ + ∂ ∂ = dy y fdx x fdf , assim comparando com a equação dada teremos ),( yxM x f = ∂ ∂ ou yx x f 23 2 += ∂ ∂ , que integrado parcialmente em relação a x resulta )(23 ygyxxf ++= . Comparando ),( yxN y f = ∂ ∂ teremos 542)('2 +−=+ yxygx . Logo 54)(' +−= yyg que integrado nos fornece yyyg 52)( 2 +−= . Daí a solução f(x,y)=C fica: Cyyyxx =+−+ 522 23 Resolver as seguintes equações diferenciais: 1) 02)( 22 =−− xydydxyx 2) 0)23()12( =−+−+− dyyxdxyx 3) 0)2( =−+ dyyxedxe yy 4) 0)cos2()( 23 =+++ dyyxydxyx 5) 0]12)cos([])cos([ =++++ dy y xxyxdx x yxyy 6) 0)73()12( =++− dyydxx 7) 0)84()45( 3 =−++ dyyxdxyx 12 8) 0)42()32( 22 =++− dyyxdxxy 9) 0dy)ylnx(dx)xln4yx3( 32 =−+− 10) 0)cos23()( 223 =++−− dyxyxydxxsenxyy 11) xy xy xey2 ye2 dx dy − + = 12) 0)3()154( 2423 =−++−− dyxyxdxyxyx 13) 1)1(,0)12()( 22 ==−+++ ydyxxydxyx 14) 2)1(,0)146()524( =−=−++−+ ydyxydxxy 15) .03131 =⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ +−+⎟ ⎠ ⎞ ⎜ ⎝ ⎛ +− dyx y dxy x 16) 0 91 1 23 2 32 =+⎟ ⎠ ⎞ ⎜ ⎝ ⎛ + − yx dy dx x yx 17) 0coscos)( =+− ydyxdxsenxsenytgx 18) xyx dx dyyx 44)221( 32 +=−− 19) eydyyxysenxdxxyxxy ==+−+−− )0(,0)ln2()23cos( 322 20) 262 xyxe dx dyx x +−= 5° TIPO: EQUAÇÕES REDUTÍVEIS A EXATAS Na equação Mdx+Ndy=0, quando as derivadas parciais y M ∂ ∂ e x N ∂ ∂ diferirem, muitas vezes pode-se determinar um fator integrante que irá transformar a equação dada numa equação exata. Vejamos o exemplo: Resolver a equação 02)( 2 =+− xdydxxy . Primeiramente, é sempre importante verificar a que modelo esta equação pertence: i. Ela não é de variáveis separáveis porque temos soma das variáveis. ii. Ela não é homogênea porque os coeficientes das diferencias são polinômios que não têm os mesmos graus. iii. Para verificarmos se a equação é exata vamos utilizar a relação x N y M ∂ ∂ = ∂ ∂ . Como x N y M ∂ ∂ ≠ ∂ ∂ , pois 1)( 2 = ∂ −∂ = ∂ ∂ y xy y M e 2)2( = ∂ ∂ = ∂ ∂ x x x N a equação também não é exata. Agora vamos determinar um fator integrante, isto é, um fator que ao se multiplicar ambos os membros da equação a transforme em exata. Seja ),( yxλ este fator integrante. Impondo que 0),(2),()( 2 =+− dyyxxdxyxxy λλ seja exata, teremos: [ ] [ ] x yxx y yxxy ∂ ∂ = ∂ −∂ ),(2),()( 2 λλ 15 Voltando para (4) determinaremos CQdxet Pdx +∫= ∫ . . Assim, resulta ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ +∫∫= ∫ − CQdxeey PdxPdx . que é a solução geral da equação linear. 2. Fator de integração O fator ∫= dxxP e )( λ transformará a equação (1) numa equação diferencial exata, isto é: Escrevendo (1) com diferenciais, vem ( ) 0=−+ dxQPydy . Quando multiplicada pelo fator integrante λ, resultará na equação exata 0.. =⎟ ⎠ ⎞⎜ ⎝ ⎛ ∫−∫+∫ dxeQyePdye PdxPdxPdx . Resolva as seguintes equações diferenciais: 1. 2−=− x x y dx dy 2. senxytgx dx dy =− 3. 0cot =−+ x gx x y dx dy 4. 0cos)1( =−−+ ydxdysenyx 5. arctgxy dx dyx =++ )1( 2 6. y dx dy 5= 7. 4123 =+ y dx dy 8. xey dx dy 3=+ 9. 223' xyxy =+ 10. 1'2 =+xyyx 11. ( ) 024 2 =++ ydxdyyx 12. ( )dxyxsenxxdy −= 13. ( ) 01 =++ ye dx dye xx 14. 1cos =+ ysenx dx dyx 15. xxy dx dyx −=+ 34 16. ( ) xeyxx dx dyx =++ 22 17. ( ) 01coscos 32 =−+ dxxyxsenxdy 18. ( ) 02 =−++ dyyexxyydx y (dica escreva dx/dy) 19. ( ) xeyx dx dyx 313 −=++ 20. ( ) 04 6 =+− dyyxydx 21. xx x ee ey dx dy − − + − =+ 21 22. ( ) 022 2 =−++ dyyxyxydx 23. θ=θ+ θ cossecr d dr 24. ( ) xyy dx dyx 4852 2 −−=+ 25. 205 =+ y dx dy , com y(0)=2 26. ERi dt diL =+ , sendo L, R e E constantes, com i(0)=io 27. ( ) xytgxy 2cos' =+ , com y(0)=-1 28. )50( −= Tk dt dT , com T(0)=200 29. ( ) xy dx dyx ln1 =++ , sendo y(1)=10 30. ( ) 022 =+− y dx dyxx , com y(3)=6 31. xy y dx dy − = , sendo y(5)=2 32. Encontre uma solução contínua satisfazendo )(xfy dx dy =+ , em que ⎩ ⎨ ⎧ > ≤≤ = 1 xse 0, 1x0 se ,1 )(xf e a condição y(0)=0 16 RESPOSTAS 1) )ln2( Cxxxy +−= 2) ) 2 (sec 2 Cxsenxy += 3) x Csenxy += )ln( 4) )2sec2)(sec( Cytgyyytgyx ++−+= 5) arctgxCearctgxy −+−= 1 6) xcey 5= 7) xcey 4 3 1 −+= 8) xx ceey −+= 3 4 1 9) 3 3 1 xcey −+= 10) 11 ln −− += cxxxy 11) 2 12 5 4 −+−= cyyx 12) x c x senxxy ++−= cos 13) 1+ = xe cy 14) xcsenxy cos.+= 15) 43 5 1 7 1 −+−= cxxxy 16) xx e x ce x y −+= 222 1 17) ecxcxy cos.sec += 18) yyyy e y ce y e y ex −++−= 224 1 2 1 2 1 19) xx e x cey 33 −− += 20) 462 cyyx += 21) xxxx ceeeey −−− ++= )ln( 22) 21 ye y c y x −+= 23) crtg +−=+ θθθθ cos)(sec 24) 4)2(63 5 + + + = x c x y 25) xey 524 −−= 26) LRto eREiREti /)/(/)( −−+= 27) xxsenxy coscos. −= 28) ktetT 15050)( += 29) 21ln)1( +−=+ xxxyx 30) 2 2 − = x xy 31) y yx 8 2 1 += 32) ⎩ ⎨ ⎧ >− ≤≤− = − − 1 xse ,)1( 1x0 se ,1 x x ee e y 7° TIPO: EQUAÇÕES DE BERNOULLI Conceito: As equações da forma nyxQyxP dx dy )()( =+ (1) com 1≠n , onde P e Q são funções de x ou constantes, são chamadas de equações de Bernoulli. Resolução: Para resolvermos a equação de Bernoulli iremos transformá-la numa equação linear multiplicando ambos os membros de (1) por y-n, o que implicará em )()( 1 xQyxP dx dyy nn =+ −− (2). 17 Em (2), chamando ty n =−1 , obteremos )().( 1 1 xQtxP dx dt n =+⋅ − que escrita como )().1().().1( xQntxPn dx dt −=−+ representa uma equação linear. Como exemplo da equação de Bernoulli, podemos citar um modelo empírico usado para a determinação do peso de peixes, que é a equação de Von Bertalanffly, 3/2pp dt dp αβ =+ , onde p é peso de cada peixe em função do tempo t, α é a constante de anabolismo, isto é, a taxa de síntese de massa por unidade de superfície do peixe e β é a constante de catabolismo, representando a taxa de diminuição da massa por unidade de massa. Resolva as seguintes equações de Bernoulli: 1. 232 xy x y dx dy =− 2. 32 xyxy dx dy =− 3. 33 yxxy dx dy =+ 4. 2 1 y y dx dyx =+ 5. ( )13 −= xyy dx dy 6. xyy dx dyx =+ 22 7. 1/2y(1) com ,32 42 ==− yxy dx dyx 8. 33 yxy dx dyx =+ 9. yxy xdx dy += 4 10. 02 2 =+− xy dx dyxy 11. 222 y x y dx dy =+ 12. dxyyxdy )1( 2 += 13. 22 )1( xyxy dx dyx +=− Respostas: 1. Cx xy + − = 4 2 3 4 2. Ce ey x x + − = 2 2 2 2 2 2 3. ( ) 2/12 21 −++= xCexy 4. Cxyx =− 333 5. xCe y x 33 1 3 1 =+−− 6. Cx y x =− ln 20 APLICAÇÕES DAS EQUAÇÕES DE 1a ORDEM E 1o GRAU 1. Determine a equação das curvas que possuem a subnormal constante. 2. Determine a equação das curvas que possuem a subtangente constante. 3. Nos problemas a seguir determine as trajetórias ortogonais de cada família de curvas dadas: a. cxy = b. 2cxy = c. 122 =+ ycx d. xcey −= e. 32 cxy = f. cx xy + = 1 g. cxyx 42 22 =+ h. θcos2cr = i. θ2sen2 cr = 4. Encontre as curvas das trajetórias ortogonais de yceyx =+ , que passam por P(0,5). 5. Um investidor aplica determinada quantia que triplica em 30 meses. Em quanto tempo essa quantia estará quadruplicada, supondo que o aumento é proporcional ao capital existente a cada instante? 6. Sabe-se que a população de uma certa comunidade cresce a uma taxa proporcional ao número de pessoas presentes em qualquer instante. Se a população duplicou em 5 anos, quando ela triplicará? 7. Suponha que a população da comunidade do problema 6 anterior seja 10.000 após 3 anos. Qual era a população inicial? Qual será a população em 10 anos? 8. A população de bactérias em uma cultura cresce a uma taxa proporcional ao número de bactérias presentes em qualquer tempo. Após 3 horas, observa-se que há 400 bactérias presentes. Após 10 horas existem 2000. Qual era o número inicial de bactérias ? 9. O isótopo radioativo de chumbo, Pb-209, decresce a uma taxa proporcional à quantidade presente em qualquer tempo. Sua meia-vida é 3,3 horas. Se 1 grama de chumbo está presente inicialmente, quanto tempo levará para 90% de chumbo desaparecer? 10. Quando um raio de luz vertical passa através de uma substância transparente, a taxa na qual sua intensidade I decresce é proporcional a I(t), em que t representa a espessura do meio (em metros). No mar a intensidade a 3 m abaixo da superfície é de 25% da intensidade inicial Io do raio incidente. Qual é a intensidade do raio a 15m abaixo da superfície? 11. Segundo a Lei de Newton, a velocidade de resfriamento de um corpo no ar é proporcional à diferença entre a temperatura do corpo e a temperatura do ar. Se a temperatura do ar é 20oC e o corpo se resfria em 20 minutos de 100oC para 60oC, dentro de quanto tempo sua temperatura descerá para 30oC? 12. Um termômetro é retirado de uma sala, em que a temperatura é 70ºF, e colocado no lado fora onde a temperatura é 10ºF. Após 0,5 minuto o termômetro marcava 50ºF. Qual será a temperatura marcada pelo termômetro no instante t=1 minuto? Quanto levará para marcar 15ºF? 13. Um indivíduo é encontrado morto em seu escritório pela secretária que liga imediatamente para a polícia. Quando a polícia chega, 2 horas depois da chamada, examina o cadáver e o ambiente tirando os seguintes dados. A temperatura do escritório era de 20oC, o cadáver inicialmente tinha uma temperatura de 35oC. Uma hora depois medindo novamente a temperatura do corpo obteve 34.2oC. O investigador, supondo que a temperatura de uma pessoa 21 viva é de 36.5oC, prende a secretária. Por que?. No dia seguinte o advogado da secretária a liberta, alegando o que? 14. Em um depósito há 100l de uma solução aquosa que contém 10kg de sal. Joga- se água neste depósito com uma velocidade de 3l/min ao mesmo tempo em que, através de um orifício desse tanque, a mistura escoa com uma velocidade de 2l/min. A mistura se mantém homogênea por agitação. Que quantidade de sal haverá no tanque 1h depois de iniciada a operação 15. Inicialmente, 50 gramas de sal são dissolvidos em um tanque contendo 300 litros de água. Uma solução salina é bombeada para dentro do tanque a uma taxa de 3 litros por minuto e a solução bem misturada é então drenada na mesma taxa. Se a concentração da solução que entra é 2 gramas por litro, determine a quantidade de sal no tanque em qualquer instante. Quantas gramas de sal estão presentes após 50 minutos? E após um longo tempo? 16. Um tanque contém 500 litros de água pura. Uma solução salina contendo 2g de sal por litro é bombeada para dentro do tanque a uma taxa de 5 litros por minuto. A mistura é drenada à mesma taxa. Encontre a quantidade de gramas de sal no tanque em qualquer instante. 17. Suponha que um estudante infectado com um vírus da gripe retorne a uma faculdade isolada no campus onde se encontra 1000 estudantes. Presumindo que a taxa na qual o vírus se espalha é proporcional não somente à quantidade de alunos infectados, mas também à quantidade de alunos não infectados, determine o número de alunos infectados após 6 dias se ainda é observado que depois de 4 dias x(4)= 50. 18. Uma lancha se desloca numa lagoa com uma velocidade de 10m/s. Em dado instante seu motor é desligado, com isso a lancha sofre uma redução de velocidade proporcional à velocidade instantânea. Sabendo que ao final de 5 segundos sua velocidade é de 8m/s, qual será o tempo necessário para que a lancha adquira velocidade de 1m/s? 19. Um bote está sendo rebocado a uma velocidade de 12nós(6,17m/s). No instante em que o cabo do reboque é largado, um homem no bote começa a remar, no sentido do movimento com uma força de 10N. Sabendo que o peso do homem e do bote é 200N e que a resistência ao deslocamento, em N, é de 2.6v, sendo v a velocidade em m/s, achar a velocidade do bote no fim de 30 segundos. 20. Uma bateria de 12 volts é conectada a um circuito em série no qual a indutância é de 0.5 Henry e a resistência 10 ohms. Determine a corrente i se a corrente inicial é zero. 21. Achar a equação da curva que passa pelo ponto P(5,6), conhecendo-se a declividade de sua tangente num ponto qualquer y x dx dy 3 2 = . 22. Achar a equação da curva cuja subtangente seja o dobro da abscissa do ponto de contato. 23. Achar a equação da curva cuja subtangente num ponto P(x,y) seja igual à ordenada de P. 24. Uma curva dada passa pelos pontos (0,0) e (3,9). Achar a sua equação sabendo que a mesma tem a propriedade de dividir o retângulo formado pelos eixos coordenados e pelas retas paralelas a estes, tomadas por um ponto P(x,y), em duas partes, sendo a área de uma dela o triplo da outra. 25. Achar a equação da família de curvas em que a subnormal, num ponto P(x,y) seja igual à abscissa desse ponto. 22 26. Um marca passo, como indicado na figura abaixo, consiste em uma bateria, um capacitor e o coração como resistor. Quando a chave S está em P, o capacitor C é carregado; quando S está em Q, o capacitor R descarregado, enviando um impulso elétrico ao coração. Durante esse tempo, a voltagem E aplicada ao coração é dada por 21 t, 1 ttE RCdt dE <<−= , onde R e C são constantes. Determine E(t) se E(t1)=E0. (É claro que a chave é aberta e fechada periodicamente para simular o batimento cardíaco natural.) 27. Em março de 1987 a população mundial atingiu cinco bilhões, e estava crescendo à taxa de 380 mil pessoas por dia. Assumindo-se taxas de natalidade e mortalidade constantes, para quando se deve esperar uma população mundial de 10 bilhões de pessoas. 28. É um fato da física que os elementos radioativos se desintegram espontaneamente em um processo chamado decaimento radioativo. Os experimentos têm mostrado que a taxa de desintegração é proporcional à quantidade de elemento presente. Sabe-se que a meia-vida específica do carbono-14 radioativo está em torno de 5730 anos. Em 1988, o Vaticano autorizou o Museu Britânico a datar a relíquia de pano conhecida como o Sudário de Turim, possivelmente o sudário de Jesus de Nazaré. Este pano, que apareceu em 1356, contém o negativo da imagem de um corpo humano que se acreditava no mundo inteiro ser o de Jesus. O relatório do Museu mostrou que as fibras no pano continham entre 92 e 93% do carbono-14 original. Use esta informação para estimar a idade do sudário. 29. Ache uma curva do plano xy que passa pelo ponto P(0,3) e cuja reta tangente em um ponto qualquer tem inclinação 2x/y2. 30. Uma bala de massa m=3.56x10-3kg é disparada para cima com uma velocidade inicial vo=988m/s, e torna-se mais lenta pela força da gravidade e uma força de resistência do ar de kv2, sendo k=7.3x10-6kg/m. Determine a altura máxima atingida pela bala.(Considere g=9,8m/s2) Coração R C Q S P E0 25 Na envoltória y = g(x), C é uma função de x e y, ou seja, C=C(x,y)≠ constante. Um ponto P(x,y) pertencente à envoltória também satisfaz a equação F(x, y, C(x,y))=0, pois pertence a certa curva da família. Neste ponto P(x,y), Ee dx dy dx dy ⎟ ⎠ ⎞ ⎜ ⎝ ⎛=⎟ ⎠ ⎞ ⎜ ⎝ ⎛ , onde : edx dy ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ é a declividade da reta tangente à envolvida e; Edx dy ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ é a declividade da reta tangente à envoltória E Derivando F(x, y, C(x,y))=0 em relação a x, vem: 0..... = ∂ ∂ ∂ ∂ ∂ ∂ + ∂ ∂ ∂ ∂ + ∂ ∂ ∂ ∂ + ∂ ∂ ∂ ∂ x y y C C F x C C F x y y F x x x F (1) Nas envolvidas, como C= constante, vem de (1): 0 ,0. ≠ ∂ ∂ ∂ ∂ ∂ ∂ −=⇒= ∂ ∂ + ∂ ∂ y F y F x F dx dy dx dy y F x F . Na envoltória, como em qualquer ponto P (x,y) Ee dx dy dx dy ⎟ ⎠ ⎞ ⎜ ⎝ ⎛=⎟ ⎠ ⎞ ⎜ ⎝ ⎛ , vem de (1) que: 0.0... =⎟⎟ ⎠ ⎞ ⎜⎜ ⎝ ⎛ ∂ ∂ ∂ ∂ + ∂ ∂ ∂ ∂ ⇒= ∂ ∂ ∂ ∂ ∂ ∂ + ∂ ∂ ∂ ∂ x y y C x C C F x y y C C F x C C F . Como C = C(x,y) ≠ constante, vem que 0= ∂ ∂ C F . Daí, a equação da envoltória é dada resolvendo-se o seguinte sistema: ⎪⎩ ⎪ ⎨ ⎧ = ∂ ∂ 0 C F 0=y))C(x, y, F(x, . EXERCÍCIOS: 1) Dar a envoltória das seguintes famílias de curvas, onde α é o parâmetro. Represente num mesmo sistema cartesiano as curvas integrais e sua envoltória: a) α α 1.4 2 += xy b) 0).2.(2 222 =++++ αα yyx 02) Determinar a envoltória da família de retas que forma com os semi-eixos positivos um triângulo de área constante igual a 20. Resposta: 1) a) y3 = 27x b) x2 + 4y =0 2) x.y=10 e e e E F(x,y,C(x,y))=0 26 Solução singular de uma equação diferencial: Conceito: A solução singular de uma equação diferencial é uma solução que satisfaz a equação, mas não é uma de suas soluções particulares. Geometricamente, a solução singular é representada pela envoltória das curvas integrais, quando esta envoltória existe. Isto decorre do fato de que em cada ponto (x0, y0) da envoltória, o coeficiente angular da reta tangente à envoltória e à curva integral corresponde a dx dy0 . Assim, os elementos x0, y0 e dx dy0 em cada ponto da envoltória satisfazem a equação diferencial F(x,y, dx dy )=0, uma vez que são sempre elementos de uma linha integral. EXERCÍCIOS: 01) Encontre a solução singular da equação dydxyx =− 21. . Represente geometricamente a solução geral e a singular num mesmo sistema cartesiano. 02) Obter a solução geral e singular das seguintes equações: a) 12 2 2 =+⎟ ⎠ ⎞ ⎜ ⎝ ⎛ y dx dyy b) y - x. dx dy = 2 ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ dx dy c) y = 2 ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ dx dy - x. dx dy + 2 2x d) y = x. dx dy - ln dx dy e) y = y.(y´)2 +2xy´ Resposta: 01) y= sen( 2 2x +C) e y = 1± 02) a) (x-C)2 + y2 = 1 e y = 1± b) y = Cx + C2 e y = - 4 2x c) y = 2 2x +Cx + C2 e y = 4 2x d) y = Cx – lnC e y = 1+lnx e) CxCy 44 22 −= e como solução singular o ponto P(0,0). 27 EQUAÇÕES DIFERENCIAIS DE 1ª ORDEM E GRAU DIFERENTE DE 1: EQUAÇÕES DE CLAIRAUT Conceito: São as equações da forma ⎟ ⎠ ⎞ ⎜ ⎝ ⎛+= dx dyf dx dyxy . Resolução: Chamando p dx dy = a equação de Clairaut fica ( )pfxpy += . Derivando a equação anterior em relação a x, teremos: dx dppfp dx dpx dx dy )('1. ++= ( ) 0)(' =+ pfx dx dp Logo p=C e a solução geral será: ( )CfCxy += Derivando a solução geral parcialmente em relação ao parâmetro C, teremos 0)(' =+ Cfx , que é a condição para obtermos a solução singular. Resolva as seguintes equações e obtenha uma solução singular: 1. 'ln1' yxyy −+= 2. 3 ⎟ ⎠ ⎞ ⎜ ⎝ ⎛−= dx dy dx dyxy 3. '' yeyxy =− 4. ( )2' 2 1' yxyy += 5. 0 2 =+−⎟ ⎠ ⎞ ⎜ ⎝ ⎛ y dx dyx dx dy 6. ( )2'3' yxyy =− 7. 01 23 =+⎟ ⎠ ⎞ ⎜ ⎝ ⎛−⎟ ⎠ ⎞ ⎜ ⎝ ⎛ dx dyy dx dyx 8. 045 =+⎟ ⎠ ⎞ ⎜ ⎝ ⎛ +− y dx dyx dx dy 9. ( ) 2'' −−= yxyy 10. 2 1 ⎟ ⎠ ⎞ ⎜ ⎝ ⎛++= dx dy dx dyxy Aplicações: 11. Achar a curva, em que a soma dos segmentos determinados sobre os eixos cartesianos pela reta tangente seja igual a k. 12. Achar a curva, em que o produto dos segmentos determinados sobre os eixos cartesianos pela reta tangente seja igual a k. 30 EQUAÇÕES LINEARES DE ORDEM SUPERIOR Tipos especiais de equações de 2ª ordem: 1º) Equação do tipo: )(2 2 xf dx yd = Solução: )(2 2 xf dx yd = dxxf dx dydxf dx dy dx d )()( =⎟ ⎠ ⎞ ⎜ ⎝ ⎛⇒=⎟ ⎠ ⎞ ⎜ ⎝ ⎛⇒ . Integrando ambos os membros, vem: dx dy = ∫ + 1)( Cdxxf dy = [ ]dxCdxxf∫ + 1)( y = [ ] 21)( CdxCdxxf ++∫ ∫ Ex: Resolva a equação 0762 2 =+− x dx yd 2º) Equação do tipo ),(2 2 dx dyxf dx yd = : Faz-se )( , xppp dx dy == , vem: dx dp dx yd =2 2 . Assim, tem-se dx dp ),( pxf= , que é uma equação de primeira ordem em relação a p, cuja solução geral desta equação é p =F(x, C1). Como p= dx dy , vem: ∫ +=⇒=⇒= 2111 ),(),(),( CdxCxFydxCxFdyCxFdx dy Ex.: Resolva as equações: a) (1+x) 02 2 =+ dx dy dx yd b) xeyy 6´´´ =− 3º) Equação do tipo )(2 2 yf dx yd = : Faz-se )( , yppp dx dy == , donde vem: dy dpp dx dy dy dp dx dp dx yd === .2 2 . Como [ ] 12 )(2)()()( Cdyyfpdyyfpdpdyyfpdpyfdy dpp +=⇒=⇒=⇒= ∫∫ ∫ . Daí vem: 31 [ ] [ ] [ ] 111 2 )(2 )(2)(2 Cdyyf dydxCdyyf dx dyCdyyf dx dy +± =⇒+±=⇒+=⎟ ⎠ ⎞ ⎜ ⎝ ⎛ ∫ ∫∫ que é uma equação de variáveis separadas em x e y. Ex.: Resolva a equação y´´+9y = 0 Ex: Uma partícula de massa m se desloca ao longo do eixo dos x atraída por outra, situada na origem, com a força F = -4mx-3, sendo x > 0. Determinar a equação do movimento, sabendo-se que para t =0 se tem x = 2 e a velocidade v = - 3 . 4º) Equação do tipo ),(2 2 dx dyyf dx yd = : Procedendo de modo análogo ao anterior, a equação se reduz a ),( pyf dy dpp = . Resolvendo-a em relação a p e substituindo pelo seu valor dx dy , obtém-se uma equação de variáveis separadas. Ex.: Resolver a equação y.y´´ - y2.y´=(y´)2 Equações lineares de ordem superior Forma: Equações diferenciais lineares de ordem superior são as equações da forma ByA dx dyA dx ydA dx ydA dx ydA n n nn n n =+++++ − − − 012 2 21 1 1 L (1), onde Ai e B são constantes ou funções de x, com i = 0 ... n. Quando B=0 diremos que a equação é linear homogênea. Resolução: Iremos inicialmente resolver as equações lineares homogêneas de coeficientes constantes. Observe que se fizermos An=...=A2=0 teremos uma equação linear de primeira ordem cuja solução particular pode ser da forma rxey = . Impondo que tal solução seja também uma solução particular da equação linear homogênea de coeficientes constantes, teremos a equação polinomial 001 2 2 1 1 =+++++ − − ArArArArA n n n n L , chamada de equação característica. Em relação à equação característica podemos ter três casos a considerar: i. Todas as raízes da equação característica são reais e distintas Sejam nrrrr ,..., , , 321 as raízes reais e distintas da equação característica, então a solução geral será dada por: xr n xrxrxr neCeCeCeCy +⋅⋅⋅+++= 321 321 ii. A equação característica tem raízes complexas Sejam bjar +=1 e bjar −=2 as raízes complexas da equação característica 001 2 2 =++ ArArA , proveniente da equação linear de segunda ordem 0012 2 2 =++ yAdx dyA dx ydA , então a solução geral será dada por: 32 ( )senbxCbxCey ax 21 cos += iii. A equação característica tem raízes múltiplas Sejam 21 rr = raízes múltiplas da equação característica 001 2 2 =++ ArArA , proveniente da equação linear de segunda ordem ByA dx dyA dx ydA =++ 012 2 2 , então a solução geral será dada por: xrxr xeCeCy 11 21 += EXERCÍCIOS: Encontre a solução geral para cada equação dada: 1. 0'y"y4 =+ 2. 0y36"y =− 3. 0y9"y =+ 4. 0y6'y"y =−− 5. 0y16'y8"y =++ 6. 0y5'y3"y =−+ 7. 0y2'y5"y12 =−− 8. 0y5'y4"y =+− 9. 0y'y2"y3 =++ 10. 0'y5"y4'''y =−− 11. 0y'''y =− 12. 0y9'y3''y5'''y =++− 13. 0y2''y'''y =−+ 14. 0y'y3''y3'''y =+++ 15. 0 dx yd dx yd dx yd 2 2 3 3 4 4 =++ 16. 0y9 dx yd24 dx yd16 2 2 4 4 =++ Resolva as seguintes equações sujeita às condições indicadas: 17. -2(0)y' e 2y(0) ,0y16''y ===+ 18. 3(0)y' e 0y(0) ,0y5'y6''y ===++ 19. 0(0)y' e -1y(0) ,05'2''2 ===+− yyy 20. 0(0)y'y(0) ,02''' ===++ yyy 21. 1(1)y' e 0y(1) ,02'3'' ===+− yyy 22. -7(0)'y' e 1(0)y' e 0y(0) ,0'36''12''' ====++ yyy Respostas: 1. 4/x21 eccy −+= 2. x62 x6 1 ececy += − 3. x3sencx3coscy 21 += 4. x22 x3 1 ececy −+= 5. x42 x4 1 xececy −− += 35 25. 0(0) x'e 0 x(0),sen22 2 ===+ tFx dt xd o ωω 26. 0 2 y' e 0 2 y ,2sencos'' =⎟ ⎠ ⎞ ⎜ ⎝ ⎛=⎟ ⎠ ⎞ ⎜ ⎝ ⎛−=+ ππxxyy 27. 2 9-(0)'y' e 2 5(0)y' , 2 1y(0) ,40242'''2''' 5 ===+−=+− xx eeyyy Respostas 1. 32 ++= −− xx BeAey 2. 5 3 5 655 +++= xBxeAey xx 3. 2/74222 +−++= −− xxBxAey xx 4. 9 2 52)62()62( −−−+= +−+− xxBeAey xx 5. 633 −+= − xx BeAey 6. 73 316 73 3cos6) 2 3 2 3cos(2/1 xsenxxBsenxAey −+++= 7. 4 55cos senxxBsenxAy ++= 8. 8/2/2/cos 2/2/2/ xxx xeDsenxxCBeAey ++++= − 9. 3 84 xxx xeBeAey −+= 10. ( )1 2 2 +++++= xxx exDxeCeBxAy 11. xexxxBsenxAy 32 )3/444(33cos −+−++= 12. xBeAy x 3++= 13. 2/22/2/ 2 112 xxx exBxeAey +++= 14. xxxBsenxAy 2cos 4 322cos −+= 15. xsenxxxxBAsenxy 2 1cos 2 cos 2 +−+= 16. 4 2)22cos( xsenxexBsenxAey x x ++= 17. 25 2cos9 25 212 2 cos xxsenxBxeAey xx −+−+= −− 18. 3737 cos6 4 2 6 senxxxCeBxAy x +−−++= 19. 3 23 3 2 x xxx exxeCxBxeAey −−−++= 20. 32coscos 2 −−++++= xxDxsenxxCxBsenxxAy 36 21. 3 2cos44cos xxBAsenxy +++= 22. 2/122 −= xseny 23. xxey x 303200200 25/ +−+−= − 24. xxx esenxexey 422 79cos10 −−− ++−= 25. ttFtsenFx oo ω ω ω ω cos 22 2 −= 26. 3 2 246 cos xsenxsenxsenxxy ++−−= π 27. 2/12291111 52 xxxx eexxxeey +−++−= MÉTODO DA VARIAÇÃO DE PARÂMETROS (LAGRANGE) Vamos desenvolver o método inicialmente para uma equação linear de segunda ordem ByA dx dyA dx yd =++ 012 2 (1). A solução característica de (1) é dada por 2211 yCyCyc += e a solução particular será dada por 2211 yuyuy p += , onde 21 e uu são funções que serão determinadas pela resolução do sistema: ⎩ ⎨ ⎧ =+ =+ Byuyu yuyu 2211 2211 '''' 0'' EXERCÍCIOS: Resolva as seguintes equações diferenciais pelo método da variação de parâmetros: 1. xyy sec" =+ 2. x eyyy x =+− '2" 3. senx yy 1" =+ 4. xgyy 3cot9" =+ 5. 2242" xexyy =− 6. senxyy =+" 7. xyy 2cos" =+ 8. xyy cosh" =− 9. xeyy x cos4" =− 10. xe yyy + =++ 1 12'3" 11. xseneyyy =++ 2'3" 12. xyy 3sec29'' =+ 13. xexyyy 2'2'' −=+− 14. xyy 2sen4'' =+ Respostas: 1. )ln(cos.coscos xxxsenxBsenxxAy +++= 2. xxeeBxAy xx ln)( ++= 3. )ln(.cos.cos senxsenxxxBsenxxAy +−+= 37 4. ⎟ ⎠ ⎞ ⎜ ⎝ ⎛⋅++= 2 3ln 9 333cos xtgxsenxBsenxAy 5. 222 xxx eBeAey ++= − 6. 2 coscos xxBsenxxAy −+= 7. xBsenxxAy 2cos 6 1 2 1cos −++= 8. 244 xsenhxBeAexexeBeAey xx xx xx ++=−++= − − − 9. )cos2( 10 22 xsenxeBeAey x xx −++= − 10. )1ln()( 22 xxxxx eeeBeAey ++++= −−−− 11. xxxx seneeBeAey 22 −−− −+= 12. )3ln(cos)3(cos 9 23 3 23cos3 xxxsenxBxAseny +++= 13. )ln1( xeBxeAey xxx +−+= 14. )21( 8 12cos2 xxsenxBxAseny −++= MÉTODO DO OPERADOR DERIVADA Conceito: Dada uma função definida por y=f(x), chama-se operador derivada, denotado por D, a dx dD = , 2 2 2 dx dD = , 3 3 3 dx dD = , ... Propriedades: Sejam u=u(x) e v =v(x): P1. D(u+v)=Du+Dv P2. D(a.u)=a.Du, a ∈ℜ P3. Dm(Dnu)=Dm+nu, com m∈ℜ e n∈ℜ. P4. O operador direto uaDuuaD .)( −=− , a ∈ℜ. P5. O operador inverso ∫ −=− dxueeuaD axax ..1 , a ∈ℜ. Exemplo: Resolver a equação ( ) xeyDD 32 65 =+− , utilizando o operador inverso. ( ) xeyDD 32 65 =+− xeyDD 3)3)(2( =−− xe D yD 3 2 1)3( − =− ∫ −=− dxeeeyD xxx 322 .)3( ( )CeeyD xx +=− .)3( 2 xx CeeyD 23)3( +=− 40 EQUAÇÃO DE EULER-CAUCHY A equação de Euler-Cauchy tem a seguinte forma: ByA dx dybaxA dx ydbaxA dx ydbaxA n n n n =+++++++ 012 2 2 2 )()()( L , onde A0, A1, ..., An, a e b são constantes. Para resolver tal equação faremos teabax .=+ , que irá eliminar os coeficientes variáveis. EXERCÍCIOS: Resolver as seguintes equações diferenciais: 1. xy dx dyx dx ydx 612)12(2)12( 2 2 2 =−+−+ 2. 0122 2 2 2 =−+ y dx dyx dx ydx 3. xy dx dyx dx ydx ln 2 2 2 =+− 4. xexyxyyx 42 23'3'' =+− 5. xy dx dyx dx ydx 322 2 2 2 =+− 6. 02'2''3''' 23 =+−+ yxyyxyx 7. xxxyxyyx 3ln2'2''' 23 +=−+ 8. )1ln(6')1(18'')1(9''')1( 23 xyyxyxyx +=++++++ 9. 4(1)y' e 0y(1) com ,0'3''2 ===+ xyyx 10. 2(1)y' e 1y(1) com ,0'''2 ===++ yxyyx Resolva as seguintes equações diferenciais por desenvolvimento em série: 11. 02 =−− xy dx dyx 12. 0' 2 =−− xexyxy 13. 0''')1( 2 =−++ yxyyx Respostas 1. 4 1 16 36 2 12 12 2 3 + + −⎟ ⎠ ⎞ ⎜ ⎝ ⎛ ++ + = xxB x Ay 2. 43 −+= BxAxy 3. xxBxAxy ln2ln +++= 4. xx xeexBxAxy 22 23 −++= 5. xxBxAxy ln32 −+= 6. 2ln −++= CxxBxAxy 41 7. [ ] xxxxxxCsenxBxAxy ln3 2 ln)(ln)cos(ln 2 2 +−+++= 8. 36 11 6 )1ln( )1()1(1 32 − + + + + + + + = x x C x B x Ay 9. 2 22 x y −= 10. )(ln2)cos(ln xsenxy += 11. 2xAxy += 12. xxeAxy += 13. ... 82 42 1 +−++= x AxAxAAy ooo APLICAÇÕES 1. Molas Um corpo de massa m é conectado a uma mola de comprimento l e constante elástica k, provocando um deslocamento s na mola, atingindo o equilíbrio. Após o equilíbrio, se a massa for deslocada de uma distância x e solta, teremos um movimento harmônico simples. Pela 2ª lei de Newton maF = . Como 2 2 dt xda = teremos: mgkxks dt xdm +−−=2 2 Mas como na posição de equilíbrio mg=ks, vem: kx dt xdm −=2 2 , (1) sujeito às condições iniciais x(0)=x0 e x’(0)=x1. Resolvendo, teremos a equação do movimento. Obs.: Quando tivermos uma força de resistência ao movimento, devida ao meio ambiente, por exemplo, vamos supor que esta força seja proporcional à velocidade. Assim a equação (1) acima ficará: s l x equilíbrio Posição inicial mg K(s+x) 42 dt dxkx dt xdm α−−=2 2 , onde α é uma constante de proporcionalidade. 2. Deformação em vigas horizontais Dada uma viga simplesmente apoiada de comprimento (vão) 2l, sujeita a uma carga uniformemente distribuída q. Para determinar as reações de apoio, poderemos associar a carga uniformemente distribuída a uma carga concentrada equivalente, aplicada no centro de gravidade da carga uniforme. Aplicando as equações de equilíbrio da Estática ( )∑ ∑∑ === 0M e 0V , 0H chegaremos a qlRR BA == , onde H, V e M são as componentes horizontais, verticais e momentos estáticos, respectivamente. Para a determinação da equação dos momentos, tomaremos uma seção S, qualquer, na estrutura. Chegando a: 222 .. 2qxqlxxqxxRM AS −=−= . Sabemos da Mecânica que M R EI = , onde E é o módulo de elasticidade, I é o momento de inércia da seção transversal, R é o raio de curvatura da linha elástica. Do Cálculo Diferencial, sabemos que 2 2 2/32 1 dx yd dx dy R ⎥ ⎥ ⎦ ⎤ ⎢ ⎢ ⎣ ⎡ ⎟ ⎠ ⎞ ⎜ ⎝ ⎛+ = . Como a inclinação da linha elástica é muito pequena, podemos impor que 0= dx dy , chegando a EI M dx yd =2 2 , que sujeita as condições de contorno y(0)=0 e y’(l)=0, nos dará a equação da linha elástica. x R A S l l R A R B 2ql 2l 45 2) ⎪ ⎪ ⎪ ⎩ ⎪⎪ ⎪ ⎨ ⎧ −−++= −+−−= +−++= − − − 10 cos 10 17 5 2 42 5 cos4 5 12 5 4 2 cos5 2 2 2 usenuCeBeAez usenuCeBey usenuCeBeAex uuu uu uuu 3) ⎪⎩ ⎪ ⎨ ⎧ +−++−= −+= x x esenxBAxBAz eBsenxxAy 2)3(cos)3( 2 cos 4) ⎪ ⎪ ⎩ ⎪⎪ ⎨ ⎧ +−−= +++= − − −− 130 61 130 cos33 3 4 )cos8( 65 1 5 3 53 senxxBeAez xsenxBeAey x x xx 5) ⎪ ⎪ ⎩ ⎪⎪ ⎨ ⎧ +−−= +−−++= − − 36 11 182 3 54 11 189 23 3 2 23 32 xxCeBey xxxCeBeAz x x xx 6) )6( 5 2y , )( 5 4 2323 tttt eeeex −− −=−= 7) 10 }3)(3cos){( 3 1y , 420 1133cos ttt etsenBAtBAeetBsentAx +++−=−−+= − EQUAÇÕES DIFERENCIAIS PARCIAIS Conceitos: São as equações diferenciais que possuem derivadas parciais de uma função de várias variáveis. A maior ordem da derivada que aparece na equação diferencial é chamada de ordem da equação diferencial parcial. Com respeito às soluções de uma equação diferencial parcial devemos citar as soluções: Solução geral que é aquela que possui funções arbitrárias, a solução completa que possui constantes arbitrárias e a solução singular que é a envoltória da família de superfícies correspondentes à solução completa. Usualmente, nas equações diferenciais parciais que possuam derivadas parciais da função z=f(x,y), denota-se p x z = ∂ ∂ e q y z = ∂ ∂ , ou seja, a equação xy y zyz x zzx 2= ∂ ∂ + ∂ ∂ pode ser escrita da forma xyyzqzxp 2=+ . As equações da forma RqQpP =+ .. são chamadas de equações lineares, onde P=P(x,y,z), Q=Q(x,y,z) e R=R(x,y,z) Determinação da solução geral: Nos casos particulares das equações lineares RqQpP =+ .. , onde P=0 ou Q=0 a solução geral é facilmente determinada por integração, vejamos os exemplos: 46 a) 34 −+= ∂ ∂ yx x z terá solução geral )(32 2 yfxxyxz +−+= b) 34 −+= ∂ ∂ yx y z terá solução geral )(3 2 4 2 xfyyxyz +−+= EXERCÍCIOS: Determine a solução geral das equações diferenciais parciais: 1. 0=+ ypx 2. zyxxp 22 ++= 3. 0=− xqy 4. xzyxp −=− 5. xz x z x z 12652 2 =+ ∂ ∂ − ∂ ∂ 6. xez x z x z =− ∂ ∂ − ∂ ∂ 542 2 7. 22 2 yx yx z += ∂∂ ∂ 8. 2 2 2xy yx z = ∂∂ ∂ Respostas: 1. )( 2 2 y y xz φ+−= 2. yxyxz −−= )(2φ 3. )( 2 2 x x yz φ+= 4. ⎥⎦ ⎤ ⎢⎣ ⎡ +−−= )(ln xx x yxz φ 5. 3 52).().( 32 2 1 +++= xeyeyz xx φφ 6. 8 ).().( 521 x xx eeyeyz −+= − φφ 7. )()( 33 21 33 yxxyyxz φφ +++= 8. )()( 3 21 32 yxyxz φφ ++= Nos casos gerais poderemos empregar o método de Lagrange, que consiste na resolução do sistema R dz Q dy P dx == , cujas soluções são u=u(x,y,z)=a e v=v(x,y,z)=b e as relações 0),( =vuφ ou )(vu φ= ou ainda )(uv φ= serão soluções gerais da equação diferencial linear, desde que pelo menos u ou v tenham a variável z. Exemplos: Determine a solução geral das seguintes equações diferenciais parciais: 1) zqypx 232 =− Na comparação com a equação linear vemos que xP 2= , yQ 3−= e zR 2= , que substituído no sistema de Lagrange R dz Q dy P dx == , resulta z dz y dy x dx 232 = − = . De y dy x dx 32 − = obtemos ayx =23 e de z dz y dy 23 = − teremos byz =23 Assim uma solução geral pode ser )( 2323 yxyz φ= 2) 0=+ xqyp 47 Substituindo no sistema de Lagrange yP = , xQ = e 0=R , teremos: 0 dz x dy y dx == De x dy y dx = obtemos ayx =− 22 e de 0 dz x dy = teremos bz = , logo: )( 22 yxz −= φ é uma solução geral. 3) zqzypxyx =−++− )2()( O sistema auxiliar é dado por z dz zy dy zyx dx = − = +− 2 De z dz zy dy = −2 vem a equação linear 12 −=− z y dz dy cuja solução é a zz y =− 1 2 Para determinarmos uma segunda equação diferencial a partir do sistema auxiliar, vamos aplicar propriedades das proporções, assim: yx yxd zyzyx dydx z dz zy dy zyx dx + + = −++− + == − = +− )( 22 , de onde obteremos: yx yxd z dz + + = )( b yx z = + Logo ⎟ ⎠ ⎞ ⎜ ⎝ ⎛ −= + zz y yx z 1 2φ é uma solução geral. EXERCÍCIOS: Determine a solução geral das equações diferenciais parciais: 1. 132 =+ qp 2. yxzqxzpy 222 =− 3. xyt t zt y zy t zx = ∂ ∂ + ∂ ∂ + ∂ ∂ 4. zqp =+ 5. 243 =+ qp 6. zyqxp =+− 7. xyyzqxzp =+ 8. 222 zqypx =+ 9. xyzxqyp 2=− 10. 1cos.. =+ xqsenxp 11. zq y xp x y =+ 2 3 2 3 12. zyxxyqpyx 322 )(2)( +=++ Respostas: 1. 0)23,2( =−− yxzxφ 2. )( 3322 yxzy +=+ φ 3. 0)3,/,/( =− zxytytyxφ 4. )( yxez y −= φ 5. )43(23 xyxz −+= φ 6. )(xyxz φ= 7. )( 2zxyxy −= φ 8. )/1/1( zxxyyx −=− φ 9. )(. 22 2 yxez x += φ 10. ⎥⎦ ⎤ ⎢⎣ ⎡ −=− ) 2 ln()ln( xtgzysenx φ
Docsity logo



Copyright © 2024 Ladybird Srl - Via Leonardo da Vinci 16, 10126, Torino, Italy - VAT 10816460017 - All rights reserved