Pat eee)Indice general 1. Ecuaciones Diferenciales Ordinarias de Primer Orden 1.1. Conceptos Bésicos : 1.2. Beuaciones Diferencisles de Variables Separables 1.8. Eouaciones que se Reducen a Varieble Separable 1.4, Ecuaciones Diferenciales Homogeneas y reducibles a ‘ Homogees Ecuaciones Diferenciales Exactas Factores Integrantes 1.7. Bounciones Diferencisles Lineales de Primer Orden... LLL, Método de Variacién de Pardmetras para una ecuacién de Primer orden 1.8. Aplicaciones a Circuitos Eléctricas RL-RO. . «. : Problemss Resueltos ...... 2.4 + Problemas Propuestos . . asset (a eta ee 2, Ecuaciones Diferenciales Lineales Ordinarias 28 2.1. Beuaciones Lineales Homogéneas de segundo orden . . . . one 23 2.2. Beunciones Homogéneas de Segundo Orden con Coeficientes Constantes . . 4 2.3, Rafces complejas de la Ecuacién Caracteristica . . Ss 25 23.1. Problema con Valor Inicial . . . . « . 9 eB 24, Rals Doble de la Beuacion Caracterfsticn . . . - ak ieee 26 2.5, Aplicaciones a Sistemas Mecéinicos y a Sistemas Elécticos, : 26 2.6. Ecuacién de Couchy- Euler ........ Deters 2.7. Método de Operadores Inversos ... . : 28, Conectores ....... 05 2.9, Beuaciones Diferenciaies Lineales Homogénees de Cualquier Orden - 2.9.1, ‘Teorema de Existencia y Unicidad ... . . 2.10. Aplicaciones a Cireuitos Eléetrioos .... +... 4. 2.10.1. Método Complejo para Obtener las soluciones Pasticulares 2.11. Método de Coeficientes Indeterminados ...... 2.12, Método de Variaciéu de Parémetros....... « ~~~ 249s4r-Beuaciones con. eoeficientes variables... Problemas Resueltos .. . . : |, Problemas Propuestos . . 3. Ecuaciones en diferencias finitas 8.1, Operador de Diferencia Finita Lineal... .- +. + 3.2. Ecuaciones en Diferencias finitas de Primer Orden 3.3. Evatt on Dikrrmci 9 Sermo Grade a8 Problemas Resueltos . . . Problemas Propuestos . . 4, Resolucién de las Ecuaciones Diferenciales por Serie de Potencia 41. Series de Potencia .... 66.55 4.1.1, Teorema fundamental 1... . Teoremas de la Sumatoria . . 418, Series Notables.......- : 2. Ecuaciones Diferenciales Importantes . . . . 4.2.1. Bouacién Diferencial de Laguerre . 42.2. Beuacién Diferencial de Hermite . . 42.3. Bcuacién Diferencial de LegendreCélcula Avanzado I FIEE-UNI 424 Método de Frobenius... 0... 00. 0c ever eects 76 Problemas Remeltos [Link] 0000s s sss eee 7 : 79 Problemas Propuestos erect Secbice Gusciepe a 5. Ecuacién de Bessel et 5.1, Funciones de Bessel de Primera Clase... . aorta 2 88 2 ie 5.2, Ecuacién Diferencial Bipergeometcica de Gauss - 89 5.21. La Ecuacién Hipergeométrica . 89 5.22. ‘Teoremas Elementales de la Funciéa Hipergeométrica 90 523. Repesenacin de varias funcons en tmins de furcinss Miprgeamtrcs 90 \< Broblomas Remseios = : z Bs + 92 Problemas Propuestas... . arbiaperocindeeers Periodic hab OF 6. Transformada de Laplace 98 6.1, Definicion de la Transformada de Laplace... 0... 20 e eee cece cece eee eee es 8 6.2. Representecién de In Transformade de Laplace... 2. sees eee e eee ee en ee eens 63. Condiciones Suet paral Bistenia den Tassformada de Laplace «600000201 98 84 Teoremas de la Transiomads de Laplace saewase = OD ‘Teorema de Linealidad a) Traslacién 2... 99 Cambio de Bacala .. . - as aan e ‘ransformade de Leplace de los derivedas « - 199 i Transfrmods de Laplace de intgrales « «« ang 7200 5. Multiplicaciéa por ¢* « 100 Division port... . = 101 Desplazamienta ent - 101 Funciones Periddicas - |» yo1 6.5. Algunas Funciones Eapecialos . . . 102 65.1. La Funciin Gomme |. 102 652. La Funcién Beta... 103, 653. Punciones de Bessel | - 103 6.54. Puncién de error. | 1038 Problemas Remeltos ...... 0... +104 Problemas Propuestos pre + 108 7. Transformada Inversa de Laplace 109 74. ‘Unleldad de le Trenetormede Inver de Leplooo « ‘ =igtiie OR 109 TAL. Teorema de Lerch . . 3 fice ws, cel 7, Meovemes de a Trensfrmads Ives de Laplace « eeeearyeed 109 Linealidad ... . : a 109 Traslacién .. .. - aera inonen gn os welll Segundo Teorema trasiacién « ee 10 Cambio de escala.... . suo ‘ransformada inversa de Laplace de las > 210 . Bacsormade inves de Laplace de ls Integral, ear ony peed . Multiplicacion por s . « EF waa um :uL 9. F :u2 73, Métodos para hallar Ia ‘Transformada de Laplace +112 73.1. Método de fracciones parciales . sapehg A 73:2. Método de las sories .... . « eave Problemas Resuellos..... +200 000 23 Problems Propuestos 5 Bans 216 8. Sistemas de Ecuaciones Diferenciales Lineales con Coeficientes Constantes uy 8.1 Bxponencal de una Matis y Solin de Hewson Dien -u7 Problemas Resueltos ..... . « ore he ses ee 29 Problemas Propuestes i fed eee lsGaleulo Avanzado I FIEE-UNI Anexo I: Operadores Inversos de Heaviside 124 Demostracion de los Operadores Inversos «6. +0 ee ee eee eee ee 126 Anexo Il: Tabla de Transformadas de Laplace at Anexo III: Examenes ciclo 2019-1 145Capitulo 1 Ecuaciones Diferenciales Ordinarias de Primer Orden 1.1. Conceptos Basicos Por ecuacién diferencial ordinaria se entiende una relacién en que intervienen una funcién no conocida y de = ¥ sus derivadas Por ejemplo: ¥ = cose ay y'+4y=0 2) yy + 267y" = (2? +2)y? 3) Son ecuaciones diferenciales ordinarias. El término ordinaria la distingue de una ecuacién diferencial parcial en la que se relacionan derivadas parciales de una funcién no conocida de dos o més varlables independientes. Por ejemplo : Uae + tyy = 0 6 uns ecuaclén diferencial parcial. En este capftulo solamente consideraremos ecuaciones diferenciales ordinarias. Se dice que una ecuacién diferencia! ordinaria es de orden n , si la derivada de mayor orden que contiene le ecuacién es a enésima derivada de y com respecto a 2. La nocién de orden de una ecuacién diferenciel conduce una elasifcacién util de las ecuaciones , de primer orden , segundo orden, ete. De este modo, (1.1)-es tna ecuacién diferencial de primer orden, (1.2) e de segundo orden y (1.8) es de tercer orden. En este capitulo consideraremos Ins ecuaciones de primer orden. Las ecuaciones de mayor orden se considerarén cen los capitulos 2 a 3. Una funeién: y=a(e) (ua) es una solucién de una ecuacién diferencial de primer orden, dada en algiin intervalo, digamos a <2 yur ») Celoutames y raemplezemes dy =2dutuds «) le ecuscion que queda en las variables x y u es de variables separables, la cual se resuelve por el metodo anterior y al final se vuelve a las variables originales reemplazando u por La ecuacién se puede poner en la forma ‘en euyo caso f(z,y) es de grado cero. Ejemplo 1 Resolver: (2+ v9) det (8y5 — ay?2) dy =0 (e+u!)do+8y? G2) dy =0 Hacemos u= (@+u)de+(u-2)du du _ stu a zy Wea) Sea yf = f(zy) La ect dienlal s homoge sf.) = (2) usu + wezite wi pue Zt te fe im» ae os hat perl ————(pauseton diferencia de-variible separable) ——~—-— Teena Integrando 1-v [Fepem lave ; acctan(y) — $n (140%) = Info| +O (4) Pecan at y= Cenen >)CaAloulo Avanzado I FIEE-UNI 1.5. Ecuaciones Diferenciales Exactas ‘Una ecuacién diferencial de primer orden M(x; y)de + N(zy)dy =0 (1.0) se dice que es exacta si el primer miembro es la diferencial exacta o total ou, Oe dun Feet ly (1.31) de alguna funciéa u(x;y). Entonces (1.10) puede escribirse: du ‘¢ integrando, obtenemos inmediatamente la solucién general de (1.5) en la forma wey) 0 (132) ‘Comparando (1.10) y (1.11) vemos que es exacta, si existe alguna funcién u(x;_y) tal que: Ou ou @ Bem 2 @ Ben 3) Supongames que M y N estén definidas y tienen primeres derivades parciales continuas en una regién en el ‘plano zy, limitada por tna curva, cerrada que no se intersecta a si misma. Entonces de (1.13) Oe Fy MO Be ay ~ Oxdy ‘Ox ~ yon y por la suposicién de continuidad, estas segundas derivadas son iguales. De esta manera aM _6N aye (1.14) Esta condicién no solamente es necesaria sino también suficiente para que Mdz+ Nady sea una diferencial total. Si (1.10) es exacta, la fimeién u(z;y) puede encontrarse por tanteos o en Ia forma sistemética siguiente. A partir de (1.13)(a) tenemos um [macs K) (15) ‘y se mantiene constante en esta integracién , k(y) desempetia el papel de ‘'constante’ de integraciéa. Para Asterina f(y), encontamnos 54, «partir de (1.15), wea (1.19)(0) para hallar $e negra. Ejemplo 1 Resolver: © + = SS ge yy S22 Sy? Solucién (? 2") de + 2rydy=0 6M Meayay-2 > oe on N@es)=2y > = aM _9N Recultando Se = Por definicién y el criterio de exactitud Is ecuacién diferencial separable es exacta (v= 27) dr + ey dy = dF,y) = % ar ar oz 7 7 Integrando: : Fa =ve— > +h) 8Céloulo Avanzado I FIBE-UNI Fm sey + 0) ray = 22y+H(y) > KG)=0 A Aye Fayavte-S4o | yusy Luego la solucién general seré pate : ve-Z4aec yh 4 e BE es ve-F ae 1.6. Factores Integrantes En algunas ocasiones , una ecuacién diferencial dada: M(aiyhde +N (a: y)dy = 0 ‘RO @6 Exacta, pero puede serlo si se multiplica por una funcién adecuada F(z;y) # 0. Esta funcién recibe el nombre de fictor integrante de (1.16) . Con algo de experiencia, los factares integrantes pueden encontrarse Por simple observacién. Con este objeto, se recomienda al estudiante recordar algunas diferenciales. En algunos casos especiales, pueden determinarse los factores integrantes de manera sistemética, tal como se verd en la seecién siguiente. Teorema Si (1.16) no es exacta.y tiene una solucién general u(x; y) = C, entonces existe un factor integrante (incluso un niimero infinito de factores). (1.18) Demostracién. Diferenciando u(2;y) = C, tenemos uy, OU, dan eet y= 0 ‘A pactir de esta fgualdad y de (1.16), ee deduce que ou ou nes Bae Ja que debe cumplirse idénticamente, Por lo tanto, debe existir una funcién F(z;y), tal que y entoness quay) Esto demuestra que F es un factor integrante de (1.16). Ademés, si multiplicamos (1.17) por una fancién H(u), 0 obtiene Ia expresién, H(w)F(Paz + Qdy) = H(u)au que también es exacta. Lo que significa que H(u)F(x;y) 3 otro factor integrante y como H es arbitraria, se sigue que existe un mimero infinito de factores integrantes para (1.16). Entonces: aN aM aN Rem GR Fe rn Sea u(a;y),un factor integrante uMdz+uNdy=0 + Mide+Nidy=0Clculo Avanzado FIEE-UNI Donde: (lat) _ 4 90a(u)) oz an aM 2 8 ay Bulle = Fe — hay = FOO ws 2y oe Bleecién de factor de integrantes de la forma u(z,y) = 2" y™ m=; usatayt) MIB Observacién: Existen otres factores integrantes que dependen de =e y , como a™yt;22mt1yntl;,_ Ejemplo Sen M(x, yd + N(x, y)dy = 0 iQué condicién se debe cumplir para obtener un factor de integracion de la forma u(z,y) = u(2ty + 2p)? Solucion Bligiendo » = 24y ay oe Eawy+y Batts Le contlciéa pedida serd 2M en Oe =e Meta) Maayan Ota) Sagun Ia condicion de factor de potencia, wv) ue) > cumple con la condicion Es um factor de integracton ule) =e! Pe 10Célculo Avanzado I FIEE-UNI 1.7. Ecuaciones Diferenciales Lineales de Primer Orden ‘Se dice que una ecuscién diferencial de primer orden es lineal s! puede escribirse en Ia forma ¥ +f@y=r(z) (4.18) El rasgo caracteristico de este ecuacién es que os lineal en y y y/', mientras que f y r pueden ser cuslesquiera funciones dadas de z Si r(z) = 0, se dice que la ecuacién es homogénea; en caso contrarlo se dice que es no homogénea.(Con r=0 para toda x). Deduciremos una ecuacién para encontrar la solucién general de (1.18) en algyin intervalo I, suponiendo que f Yr son continuas en I para la ecuacién homogénea . ¥+s(ay=0 (119) resulta muy sencill, Separando variables, tenemos 24 = ~ (2) y, por lo tanto tafp| = — f ‘Ha)de- +e" 0 200 vale SI" 5 (ome (120) Para resolver la ecuacién no homogénea (1.18), la eseribiremos en Ie forma (Ufy-r)dz+dy=0 emostraremos que es posible hallar un factor integrante F(2) que depende solamente de =. Si ese factor existe. F(@)(fy —r)dz + F(@)dy =0 debe ser exacts. Ahora, para esta ecuacién, la condicién (5) de Ie seeci6n 1.5 toma le forma Hev-ni=Z oven Sieparends mite Inkegrnde,citenemas w(e)= f sear de esta igualdad Fle)=eM*) donde wa) = f sate "Esto demuestia que F(z) es iii factor integrate We (1-18) ~Abora;al-multiplicar (1.18) por este factor; encon- aed ‘tramos Ay +f) Puesto que A = f, podemos escribir Zoe) =er ae Ahora podemos integrar ambos miembros y obtener whe [Prise Dividiendo ambos miembros entre e*, llegaros a In ecuacién deseada v2) 8 f Arae-+o] ; he fi fle)dz (21) representa la solucién general de (1.18) en Ia forma de una integral. (La seleccién de los valores de las cons- tantes de integracién no tiene importancia). uCalculo Avanzado I FIEE-UNI Ejemplo ‘scuaciones diferenciales lineales de primer orden Resolver (ain 2 — 22y?) do + 82°ydy = 0 ‘Solucién tomes aS usy > du=3ydy Reemplazando (ett—20u) de =4du=0 (t= ru) deedu=0 Sit=ng 4 met (te%—2u) det du = 0 Beuacién diferencial respecto a u(t) Factor de integracién ef ~## = a) = if wie tetat sc] uli) = [tet + el] + Ce u(t) = =te be +Cet yo = Pine +294 C2? 1.7.1. Método de Variacién de Pardmetros para una ecuacién de Primer orden Existe otra forma interesante de obtener Ia solucién general de la ecuacién diferencial lineal vtseware) (1.22) Hemos visto que 53. dee te Hoe i gs (2) = en Sie? (2.23) ‘es una solucién de la ecuacién lineal homogénea correspondiente. Trataremos'de determiner tuna funcién u(z) ‘alae 8 (a) =ufeo(z) os yf a ea yl (2.24) sea la solucién general de (1.22). Bste intento lo sugiere la forma de le solucién general p(x) de la ecuacién hhomogénea y consiste en sustituir el pardmetro c por wna Variable u(x). Por lo anterior, este método, debido a Lagrange, recibe el nombre de método por variacién de parémetros. - Puede generalizarse para ecuaciones de mayor orden donde adquiere una gran importancia. Nuestro caso presente 1 Por su sencillez, nos brinda una buens oportunidad para dar a conocer este mnétodo fundamental. Sustituyendo (1.24) en (1.22), obtenemos: Mews yah aye fe vis} se) vutu+fu)ar y's fy) = 9 ‘Puesto que v es una solueidn de la ecuacién homogénea, la ecuaclén anterior se reduce 8 ‘La integraciéa conduce a ‘Asf obtenemos ea ¥y de acuerdo con (1.23) vemos que este resultado es idénticn al (1.20) de la seccién anterior. 2DCAleulo Avanzade I FIEE-UNI 1 33 Ejemplo yepes) ye foo Resolver Ta ecuaci6n diferencial de Bernoulli ¥tprcy= fey ineR Solucién Fy + pong’ = fey Sea z=y v= (ny ty # = (= ny "[f@)y" ~ pe)o) 2 +(1=n)ple)z = f(@)(L-n) oT ef d-mpterae i fOr. (a) (1nd tC 1.8. Aplicaciones a Circuitos Eléctricos RL -RC Tas ecuaciones diferenciales lineales de primer orden tienen diverses aplicaciones en la fisica y en le préctica de la ingenlerfa. Para ilustrar lo anterior, consideremos algunos ejemplos modelo relacionados con los circuitos eléetricos. Estas y otras consideraciones Semejantes son importantes, debido a que ayudan a expresar situacio- nes fisicas en términos de relaciones matométicas. La transicién de la situacién fisica a un modelo matemético correspondiente siempre es el primer paso en las matematicas de ingenieria y este paso importante requiere ‘experiencia y entrenamiento que solamente puede adquirir considerando ejemplos tipicos de varios campos. El cireuito eléctrico més simple es un circuito serie que tiene una fuente de energfa eléctrica(fuerza electromo- ttiz) tal como generador o baterfs y un resistor que utiliza la energfa , por ejemplo uns bombilla eléctrica (Pig. 1). i) @) : EO, Si cerramos el interruptor, entonces fuiré una corriente eléctrica J a través del resistor, provocando una caida de voltaje, esto es, el potencial eléctrico en los dos extremes del resistor serén diferentes; esta diferencia de potencial o caida de voltaje se puede medir con un voltimetro. Los experimentos demuestran que se cumple la siguiente ley. "La caida de voltaje Bp a través de un resistor es proporcional a la corriente instanténea I, mateméticamente. Eg= RI (Ley de Ohm) (1.25) donde la constante de proporcionalidad R recibe el nombre de resistencia del resistor. La corriente J se mide en arperics, Ie resistencia R en ohmios y el voltaje Bq en voltios. Los otros dos elementos importantes en circuitos més complicactos, son los inductores y capacitores. Un inductor se opone a un cambio en la corriente y tiene un efecto de inercia, en electricidad, semejante al de la masa en Ja mecénfea; mas adelante consideraremos esta [Link] experimentos nmuestran que se cumple Is siguiente ley. La cafda de voltaje Bz a través de un inductor es proporcional a la razén de cambio instanténes de la corriente [Link] respecto al tiempo, mateméticamente. aeCalculo Avanzado I FIEE-UNI a matg (1.26) donde la constante de proporcionalidad 1 se llama inductancia y se mide en henrios; el tiempo t se mide en segundos. ‘Un capacitor es un elemento que almacena energie. Los experimentos muestran que se cumple la siguiente ley. La eaida de woltaje Bo a través de un capacitor es proporcional a la carga eléctrica instanténes Q que tiene el cepacitor, mateméticamente. 1 Fo= 5a donde C es la capacitancia y se mide en faradios; la carga Q se mide en columbios. Y puesto que aQ no=8 Entonces puede escribirse: Fo=h [evar 29 La corriente 1(!) que circula en un circuito, se puede determinar resolviendo Ia ecuacién (0 ecuaciones) que results de aplicar la segunds ley de Kirchoff. La sums algebraics de todas las cafdas do voltaje instanténess slrededor de una trayectoris cerrada en un cireulto eléctrico es cero o bien el voltaje aplicado en una trayectoria cerrada en un eireuito eléctrico es igusl a Ja suma de las cafdas de voltaje en el resto de la trayectoria Ejemplo (Circuito RL) = ‘Al eplicar Ia segunda Ley de Kirchoff en el *cireuito RL” de la 6i#2 ¥ Zon (1.25) y (1.26), obtenemos < et a id= BH) Pee (1.28) SFA lene Sets oe He dala eolucién general = cons, entonets ls ecuacl6n (4) do la seein 1.7 08 paiBi fg / +0} 1 = +00 1) =R/L 7.29) Bl ltimo término tiende a oxo cuando t tiene al infnito, de modo que I()tende al valor Iimite 22, desputs {etn tiempo sufcientemente largo, I, ser prcticamente constante. ~~ n= Ba) que corresponde a la condicién inicial 1(0) = ‘Caso B (Fuerza cloctromotriz periédies).Si E(t) = Easenut, ontonces por la ccuacién (1.28) de Is Seccién 1.8, la solucién general de (1.28) (de esta seceién) es [evnautar+e] -(@= R/L) Integrando por partes se obtiene 110 = Ceo" ay Fh (sent — ut eoeut) Se puede escribir: I) = oo Ey wh oye. ee 1.90) + grea eet 0), arcian(S) (30) El término exponencial tlende a cero cuando ¢ tiende al infinito. Esto significa que después de un tiempo six ficientemente largo, la corriente I(t) précticamente elecuta cscilaclones arménicas . y se presenta el éngulo de 14Célculo Avanzado I FIEE-UNI fase ¢ como une funcion de SE. Si L = 0 entonoes o = 0 y las oocilaciones de I(t) estén en fase con £'() Se dice que un sistems eléctrico (0 dinémico) se encuentra en estedo estacionario eusndo las variables que deseri- ben su comporiamiento son constanies o funciones periédicas del tiempo y se dice que su estado es transitorio (o estado no estacionario) cuando no se encuentra en el estado estacionario. Las variables correspondientes reciben, el nombre de "funciones del estado estacionario y funciones transitorias”, respectivamente, En el ejemplo 1, caso A, la funcién Eo/[Link] la funcién de estado estacionario o soluci6n del estado estacionario de (1.28) y en el caso B, le solucién del cstado estacionario esté representada por el iltimo término de (1.30). Antes de que el circuito alcance (précticamente) el estado estacionario, se encuentra el estado transitorio. Es evidente que este periodo transitorio o intermedio se presenta, debido a que los inductores y capacitores alma- conan energis y no pueden cambiar instantdneamente las corrientes, Situseioues semejantes surgen en. varios sistemas fisioos. Por ejemplo, si se opera el interruptor de un radiorreceptor provisto de tubos al vecio, del tipo de flamento, deberd transcurrir un intervalo de tiempo trausitorio durante el cusl los tubos cambian de frfos a calientes. Ejemplo (Circuito RC) c K———_ i) Aplicando la segunds leyde Kirchoff, y las eeuactones (1.25) y (1.27) al elrculto RC de la Figura, cbieneros la ecuacién: ard [1a 2) ee “i (st) Para librarnos de la integral, derivamos con respecto a t, entonces: . ae Z s eee es Y¥ de acuerdo con (1.28) de la Secci6n 1.7 , la solucién general de (1.32) es : eile fc HB ’ noqe%eg [oe Z +0) i: (133) Caso A(Fuerza Blectromotriz constante).- Si E es constante , entonces dE /dt = 0 y toma la siguiente forma T() =e-% (asa) Caso B (Fuerza Electromotriz, senoidal) - $i B(t) = Fysenut, entonces rr Gy 7 eB ocast Sustituyeado esta ecuseién en (1.33) ¢ integrando por partes, encontramos. H = ce" + eure +uRC sent) (135) Ee eal =e) WE + (wk) molt 15Célculo Avanzado T FIEE-UNI Donde ¢ = —1/wRC. Bl primer término decrece constantemente a medida que crece ty el tltimo término representa la corriente de estado estacionario, la cual es senoidal. La gréfica de I(t) es semejante a la de la Fig.45. Circuitos més complicados y la analogfa entre las vibraciones elSctricas y mecdnicas ce considerardn en relacién con las ecuaciones diferenciales de segundo orden que se estudian en la siguiente Seccién. Ejemplo 1 En el circuito mostrado, el interruptor $ se conecta al punto a en luego de t = seg. Se conecta al punto 6, calcular q(t), ¥é 2 0 = 0 (C esta inicialmente descargado) y Ry=10, Ry=50, C=0.2F; E(t) = 10ces(St)voltios Solucisn #0 [0.5] RA + ssa as EW) ‘2da ley de Kirchoff B= RK) +h BW = Ria () + 5a) 1 FO = RAO + Gale) A(t) +540) at) =e (/ 1Ocos(st}eP*at + 2) a(t) = sen (St) + cos( 5) + 0e-** Weos(5t) ent=0 g0)=0 + G=-i aft) = sen(5#) + cos(6e) + e-** 16Célculo Avanzado T FIEE-UNT mire c= Ida ley de kirchoft Jatt) + Rail) = 0 s(t) #9) =0 a(t) +l) =0 gate [f+] at) = Ket q ( Se l-e*P=Ket = (Ie) ene a= (Ine) eo Finalizando sen (5t) + 008(64) — ete [05 7] irCAélculo Avanzado I FIEE-UNI Problemas Resueltos 1, Resolver: Solucién Boy? + 4x (de + 1) y — (82 + 40? -1) Vemos por simple inspeccion vaste A yale? Reemplazando en la ecuscién dads Ata! = 4e(4e + 1)(u + 2) — 8e( + 2) — (82° + 427-1) ote = Bez? w+ Amu = 82 uae lee [fexel seae-+c] 2+ce% ptartce™™ 4 cambiamos de variable y calculamos ate Y= 2 Te (e*Inz — 22y?) de +32%y"dy = 0 a! ingds + 32%y4dy ~ 2ey"de = 0 eet aie Bayi 3 Resolver: (¢-ylny+ylnz)dz+2(Iny—Inz)dy=0 Solucién , z—ylny+ylns'+ 2 (Iny—Inz)y =0 a) faze an (H ey) av tun Hacemos y= or + y=av+u Integrando vinv—9 Infa| +e TaQ)faweve 8Célculo Avanzado I FIEE-UNI 4, (223 + 3ye — Tx) de — (30%y +24? — 8y) dy = 0 Solucién (22? + 8y? — 7) xd — (3x7 + 24? —8) ydy =0 zap > de=%ydy . tad + dt=2ede (45-7 £-(se4 22-8) E=0 5 dz _ 2432-7) og @ +a 8 Hecemst=uth A zautk + dt=du A deady e ge + Sut ao vem + vawen a k=1, h=2 Be (B=), _ 9d Te eG Integrando la ecuscién difercncial de variable separable Sarctan(A) — In (1+ 9?) = 2Inful +O. sarctan ( ) = (+ (23) - la [2-2] +0| 5. Resolver: (22 — 2y +1)? dy = (2 -y+1)'de Solucion dy a-ytl)\* a (Saas) Haciendot=2z—-y + =1-y¥ _ 4a Get Gant 7% weatne Ney) =2 1a ecuactén diferoncial no es exact Calculo de un factor integrante Siu=u(e) 1 (am =) = (SN Nae Jn oeCaleulo Avanzado 1 FIEF-UNT enna! S )* x7 Para u=u(s,y) See M(2,x)ér + N(z,u)dy uns ecuacion diferencial no exects , hacemos w= u(z) y 2 2(2,4) Para u= ult) ligtendo 2 = ay 1 rettando (2,9) = = Ae) = Ressvindo la ecutelén diferencia exact (st ae yismsdy oh yy , ae) _ ae aye ae 0 zie ee Jee Sarno 0) Establecer la ecuacion diferencia lineal de Ler order para V(t) 8) Hallar VQ), s(t) SoluciéaCéleulo Avanzado I FIEE-UNI R a(t) R a) = Fe i) oR ss v(t) E(t) f mae ache : iQ=itiz Usendo la 2da ley de Kirchoff = ene, — BQO=2RKt) + RAO (1.36) wiles. Rint} [a (odo= Rito sn i f ocr mr fav + a=ovw (2.38) De (1.36) ee) E(t) = RGix(t) + ia(t)) + R(t) E(t) = 2Rig(t) + Rial t) (1.39) De (131 on Riz +0 = RA) (2.40) De (1.40) en (1.99) BO) = 2(Rin() + v(0) + Rial) Blt) = 3Rin(t) + 20(¢) De (1.40) E(t) = 8RCv'(t) + 20(8) (0+ gage = 2 b)SIR=1 C=1/3 Bt) =c*sen(t) v'(t) + 2u(t) Saar Ef V(t) = —cos(i)e"* + Ke-™ EDO. Como C esta inicilmente descargado v(0) = Z4{0) = 0 pues g(0) = 0 Luego en el enterior ecuactén uf) =(1—cos(t))e™ ve2o De (1.40) en (1.3¢ mae Eid) = Ri(t) + Ria t v(t) E(t) = Rift) + ROW(e) +0(0) ty = BOQ RO - wf) q- eee Como R=1 , C=1/3 a = BW) -o)— FH) i(@) = i sen()— 3a ~ ext] ox a1CAlculo Avanzado I PIEE-UNI Problemas Propuestos & Resolver (1+2°y*) y + (ay—1)'a/ =0 2, Resolver da? — ay +9? byl (e? — av+4y?) =0 . Recolver (xe +14) de + (Bay + sen(u)) dy =0 4. Establecer que condicion debe cumplirse para que la ecuacion diferencial: M(z,y)dz + N(e,y)dy =0 42 — ay! Tooge un ator integrate de forme y = u (== 3P*) Resolver (2552 Sy?) de + (T—S2y?) dy =0 Resolver 3 (2? + y?)" (edy — ydz) + fr +22%(2? + y?)'] (ade + yay) =0 Resolver (p?sen(2) — ycos®(2)) dx — (sen(z)eos(z)) dy = 0 Hallar uns condicion necesaria para que la ecuscion diferencia] Mla,y)ds + N(e,9)dy = 0 ‘Tenga un factor integrante de le forma y = u (xy? + z%y*) Resolver (ay + 2ry(lny)* + ylny) de + (20*iny +2) dy =0 10, Resolver yf + Jy= ty? LL. Resolver (y? sen(2) — yeos?(z)) dz — (sen(z) cos(z)) dy = 0 12. Resolver yf = —2e~ + 2yer* — 22? 13, Resolver (16zy + 14y?) dz + (102? 4 212y) dy =0 14. Resolver x (a — 1) y'— (22+ ly ty? +22=0 15. Encontrar el factor de integracién F y resolver: ¢) dz +e!-Fdy =0 ¥) yeeileyyes— (22 - zsmn()) et ©) y+ ay)ér + 2edy =0 d) Demostrar que si Py Q, de (1) satisfacen le relacién oP _ 8Q ee ) Demostrar que F.= y es un factor integrante de yds + 2zdy = 0. A continuacién encontrar otscs factores integrantes de'esta. ecuscién. 16. Demostrar que no tiene importancia le seleccién dal valor de Ja constante de integracién en f fdz 17. Resolver los siguientes problemas con valor inicial 8) of ~ yoote = 22 — steotsy9(§) = +1 ) ¥— (Lt Dye etaiy(l)se-1 22D *kostte” Capitulo 2 Ecuaciones Diferenciales Lineales Ordinarias ‘Las ecuaciones diferenciales ordinarias pueden dividirse en dos grandes clases, las llamadas ecuaciones lineales y les ecuaciones no lineales. Mientras que las ecuaciones no lineales (de segundo grado y érdenes superiores) son dificles, las ecuaciones lineales son mucho més seneillas en muchos aspectos, porque varias propiedades de sus soluciones pueden considerarse de manera general, y, [Link] tanto, se tiene métodos ordinarios pare resolver muchas de estas ecusciones. En el presente capitulo consideramos las ecuaciones diferenciales lineales y sus aplicaciones. Estas ecuseiones juegan un papel muy importante on las matemétieas de ingonieria, por ejemplo, en relacién con les vibrociones mecénicas, cireuitos eléctronicos y redes. 2.1. Ecuaciones Lineales Homogéneas de segundo orden Electudiante ya conoce las ecuaciones diferenciales lineales de primer orden; ahora definiremos y consideraremos las ecuaciones lineales de segundo orden. Se dice que una ecuaciones diferencial de segundo orden es lineal, si puede escribirse en la forma y+ f@)y +9@)y=r(2) (a) Las caracteristices fundamentales de esta ecuscién os que os lineal en y y sus derivadas, mientras que fg, ¥ t pueden ser cualesquiera funciones dadas de x. Cualquier ecuacién de segundo orden que no puede escribirse en la forma (2.1) se dice que es no lineal. Por ejemplo, +4 =e Feene (22) x'y" + ry +(x? -1)y=0 (2.3) ‘son ecuaciones lineales, mientras que vty =0 we ‘son ecuaciones no lineales. Las ecuaciones diferenciales de segundo orden juegan un papel bésico en muchos y diferentes problemas de Ingenieria, Veremos que algunas de estas ecuaciones son muy sencillas porque sus soluciones son funciones elementales. Otras son més complicadas y sus soluciones son funciones especiales muy importantes, como las funciones Bessel y las hipergeométricas. Las consideraciones que presentaremos a continuacion demostrardn que la ecuactén (2.1) tiene caracteristicas algo més sencillas cuando ls funci6n r/z) del segundo miembro es idénticameate cero. Entonces la ecuacin (2.1) viene a ser w+ S(@)y' + (zy =0 (2.4) y entonces se dice que es homogénea. Si r(2) 70, se dice que la ecuacién (2.1) es no homogénea. Fstas nociones son semejantes 2 las presentadas én le Sec. 1.7 en relacin con las ecuaciones lineales de primer orden. Por ejemplo, (2.2) es una ecuaciéa no homogénea mientras que (2.3) es homogénea. 23,Célculo Avanzado FIEE-UNI Las funciones f y 9 de (2.1) y (2.4) reciben el nombre de coeficientes de las ecuaciones. En todas nuestras consideraciones supondremos que x varia en algtin rango arbitrariamente fjado, por ejemplo, en algin intervalo finito o sobre todo el eje x. Todas las supasiciones y afirmaciones se referirdn a ese rango fjo, cl cual no se especificaré necesariamente en cacia caso, Se dice que tna funcién v= $2) ¢[Link] solucién de una ecuacién éiferencial (lineal o no lineal) de segundo orden, sobre algiin intervalo (post- Dlgmente infinito), si #(2) esta definida y existe su segunda derivads en todo el intervalo y, edemés, I ecuacién se transforma en tna identidad cuando se sustituye la funcién y no especificada y sus derivedas por @ y sus derivades. ‘Teorema fundamental Si se multiplica por cualquier constante una soluci6n de Ia ecuacién diferencial lineal homogenea (2.4) , la fun- ci6n resultante también es una solucién de (2.4) . Si se suman dos soluciones de esa ecuacién, la suma resultante ‘también es tna solucién de esa ecuacién. ‘La demostracién se obtiene inmediatamente de los teoremas de linealidad de las derivadas, las cuales aseguran ‘que Is derivada de una constante por una funcién es igual a Ia constante por In derivada de la fumcion y In derivada de tna suma de funeiones es la suma de sus derivadas. En efecto, para probar la primera afirmacién del teorema, supongamos que ¢(2) es una solucién de (2.4) y demastraremes entonces que y = ef(2) también es una solucin. Si sustituimos y = efx) en (2.4), el primer ‘miembro viene a set of" + feb +906 = di" +16 + 99) Y como ¢ satidface (2.4), la expresin que se encuentra dentro de los paréntesis es igual a cero y queda demostrada la primera. afirmacién del teorema. La demostracién de la iltima afirmacién es igual de sencilla, por lo que se deje al lector. El estudiante debe conocer este importantisimo teorema tanto como saber su taba de rmultiplicar, pero no debe ‘olvidar que el teorema no se cumple para las ecuaciones lineales no homogéneas o para las ecuaciones no lineales. 2.2. Ecuaciones Homogéneas de Segundo Orden con Coeficientes Constantes Ahora consideraremos las ecuaciones diferenciales lineales de segundo orden con cosficientes constantes. Estas ecuaciones tienen importantes aplicaciones en la ingenier‘a, especialmente en relaci6n cn Ine vibraciones meedal- as y eléctricas ‘Empezaremos por las ecuaciones homogéneas, es decir, ecusciones de Is forma: y' tay +by=0 (25) donde a y b son constantes. ‘Supondremos que a y b son reales y que el dominio considerado para x es todo el eje x. ;Cémo podremas resolver Ja ecuacién (2.5) Recordemos que la solucién de la ecuacién lineal homogénea de primer orden con coeficientes constantes. vtky=0, ‘una funcién exponencial, a saber yoo Es natural supone que y=e* 6) puede ser una selucién de (25) si se escoge A apropiadamente, Sustituyendo en (2.5) la solucin (26) y sus derivadas wade yy" aXe obtenernos (arte =0. ‘De modo que (2.6) es tuna solucién de (2.8) , si Aes une solucién de la ecuacién euadrética M+ai+o=0 7 ryGélculo Avanzado T FIEE-UNi Esta ecuacién recibe el nombre de ecuscién caracteristics (o eouacién auxiliar de (2.5) . Sus rafces som: MaftetVe—B), =H ve-B) (28) Apayéindonos en el razonamiento seguido, se deduce que las funciones Der (29) m= A yw son soluciones de (2.5) . El estudiante puéde comprobar esta afirmacién sustituyendo (2.9) en (25) . De nuestros conocimientos de élgebra elemental sabemcs que, si a y b son reales, la ecuscién caracteristice puede tener Caso I dos rafces resles distintas Caso II dos raices complejas conjugadas Caso TIL una rafz (real) doble. Se discutirén separademente estos casos. 2.3. Raices complejas de la Ecuacién Caracterfstica 2.3.1. Problema con Valor Inicial * En Ia Sec 2.2 vimos que las soluciones do lo ecuacién lines! homogénes con coeficientes constantes. vitoyty=0 5 (ey breales) (2.10) procedimiento puramente algebraico,es decir, determinando las rafces A1, Ao de Ia ecuacién caracteristica. +ar4+b=0 (2.41) Yy estas soluciones con de la forma =e, = et (2.12) Si, £ Ag el cociente y:/y2 no es constante y, a partir del teorema I de la Sec. 2.3, se dice que estas soluciones ‘constituyen un sistema fimdamental para toda [Link] solucién general correspondiente es y= qe" toe" (213) Si estas rafces distintes son reales (Caso I de la Sec. 2.2 ), entonces las soluciones(2.12) gon [Link] amabargo,si di y Aa son complejos conjugadas (Caso TI de la See. 2.2), digamos M=ptig , A=p—iq (p,q reales, ¢ #0) eutonces las soluciones (2.12) son complejs, Z oni gle tae geez aletader = er y aplicando le segunda ecuacién de Buler, uy = et? = oF (cos qx + isingz) a yn = cP = * (cos gr —isingz) De estes 2 expresiones para yn y y2, obtenemos 1 Fn - a) =e oengs ‘Vemos que las dos funciones de los segundos miembros son reales y , 8 partir del teorema 1 de la Sec. 2.1 , se sigue ue son soluciones de la ecuacién diferencias (2.10) -Ademés, como su cociente no es constante ,son linealmente independientes en cualquier intervalo y , por lo tanto, constituyen un sistema fundamental de soluciones de la ‘ecuacién(2-10) en todo el eje x. La solucién general correspondiente es va) = oF (Acoses + Bscnez) donde A y B son constantes arbitrarias. 25CAleulo Avanzado I PIEE-UNI 2.4. Raiz Doble de la Ecuacién Caracteristica Para completar nuestra discusiém sobre les sohuciones de las ecuaciones diferencialeslineales homogéneas con coeficientes constantes, de la forma vita +ly=0 (24) ‘ahora consideraremos el caso en el que la ecuseién caracterfstica NM +aN +bA=0 (2.18) tiene una rats doble o repetida. Este caso de conoce con el nombre de caso [Link] presenta si, y solamente si, 1 discriminante de (2.18) es igual a cero: = 4=0 y, entonces, b= *S 7 Lae raz es —e/2, y, en principio, obtenemos solamente una solucién nae 3s (210) ‘Teorema 1. En al caso de una rats doble de (2.15) , otra solucién de (2.14) es ee (17) ‘Las dos soluciones (2.16) y (49 constituye una sistema fundamental y ls solucién general correspondiente es ,por Jo tanto, y= (ei toza)e* (2.18) 2.5. Aplicaciones a Sistemas Mecdnicos y a Sistemas Eléctricos Ahora consideramos una aplicacién muy importante de les ecuaciones diferenciales lineales homogéness con. coeficientes coustantes, en I [Link] problema se tomara de le mecdnica , pero veremos posteriormente ‘que en los circuitos eléctricos se tiene una analogia completa. La ecuacién caracteristica correspondiente es ‘Los rafees con =e oo om La forma de la solucton de ls ecuscin dependeré del amortiguamiento y como en Ia Sec 2.2, podemos distinguir Jos tres casos. Caso I. 2 > 4mk Raices reales y distintes (Sobreamortiguamiento) (Caso IL. c? < dmk Rafces complejas y conjugadas (Bajo amo: (Caso IIT. c? = dk Rats real doble (Amortiguanaiento eritico) 2.6. Ecuacién de Cauchy - Euler La Ecuacién de Cauchy o Eeuacién de Euler. zPy" tary +by=0 (a,b constantes) Gay ‘También puede resolverse sigulendo un método exclusivamente algebraico. Sustituyendo = (2.20) 26— Sewn rrr wr wr rw ww rr wr rr rrr ww wrrwewrwewrwruvurwvrwvwTweT wwe ee Calculo Avanzado I FIEE-UNI ‘yas [Link] la ecuacién (2.19) , se obtiene. 2mm —1)2"-* + oxma™—1 + be™ = 0 Dividiendo m2 +(a—1)m+b=0 centonces: y=as™ te” Ejemplo 1 ‘ Resolver: [2e-5)°Dt + 2(¢-3)D? +60] 2() = 241-2") Soluetén 2") =Die(t) + [(2e—9)°D4 4 ¢-3)D* + 6D + D4] a(t) = 28-41 [(2t - 5)*D* + (2 - 5)D? + 6D] a(t) = 2041 Maultiplicando por (2 — 5): [2-54 + (2t — 5)*D? + 6(¢ — 5)D] 2(¢) = (2+ + 1)(2t — 5) + Ecuacion de Cauchy -Buler Cambio de variable: (2t—5) =e" + 2t+1= e+ +6 y por Cauchy - Euler [2'D(D- (D2) -3) + PDD 1) +6- 20] alu) = Me" +6) fax(u): [16D —96D* + 1800 — s8D]x(u) = 0 reD 16r* — 9673 + 1807? — 88r = 0 Les soluciones son: m0, r2=0754 , ry =2623+064Ui y ry =2.623-0.644i > pu) = Or + Coe ™ + Gye? cos 0.644 + Cue? son 0.644 2.7. Método de Operadores Inversos _Laecuacién diferencial de n-fsimo orden se escribe dela siguiente manera PH an ay? tanay 2 ++ bay +0 = f(z) Hacemos: Dy=¥ A D¥ymy” Dy =a Reescribiendo la ecuaeién diferencial (D* + on-1D*"* +--+ + aD + a9) y = F(z) Pow = f(@) P(D) coincide con P(r) = 9" ag a7) + ag ar? bees bar tag Polinomio asociado a la ecuacién diferencial Fectorizando resulta, P(D) = (D—1:)(D—1)...(D tn) TEOREMA Ses la ecuacién diferencial P(D)y = 0 (D-n)(D- 12)... (Dra) y= 0 (2.21) 2CAleulo Avanzado I EIERUNE tiene por solucién el siguiente conjunto fe" ee, .. 7} jemmplo 1 Resolver: y/’— sy’ +6y=0 Sotucién (D?-5D+6)y=0 (D-3)(D-2y=0 + {} ta = Ac™ + Be ‘TEOREMA ‘Una solucion de Is ecuseién no homogé P(D)y = f(x) (D-n) (D129)... (D-tm) y= fle) ‘Esta dada por 1 % = Bey l®) Conocida como solucién particular TEOREMA, La solucién de le ecuscién: ¥f — ry = f(e) Esté dade por y=e* | 7 f(2)de Ejemplo 2 Resolver: y' ~8y = xe® Solucin Solucién homogénea y¥-3 (-3)=0 n=3 4 ye 3 made Solucién particular Por el teorema, Hae / e™ (xe) dr t=O" (ze 6") p= (e—De* Sumando yp y Yp ule) = Ae + BG — 1) 28CAlculo Avanzado I IEESUNT Bjemplo 3 Resolver: —4y = e% Solucién Solucién homogénea v-4y=0 (P-)=0 + (-2r+2)=0 n=2 > nae* as > in un = Ae™* + Be Hallando la soluci6n particular w= ae) mgt i 22) »- pla )] mimi xf .-2 (,22) wale free a DD wae fe (et)ae ae w= (3-4) Sumando yn ¥ Yp via) =Ae* + Be 40 G = Je TEOREMA (Desplazamiento exponencial) P(D) = ,D™ + by aD" ovet yr tb = PO = PD ay] i 1 ety we. Py ETH? Ejemplo 4 Hallar yy on la ecuacién diferencia: (D* -3D +2) y= 2e* Solucién — wary” L v ya. 2 (D+1)? -3(D +1) +2) Seen = DID) 29Calculo Avanzado I FIEE-UNI Recordando : intra ten.. t-r Entonces 5 Uae (L+ D+ D + D> +...) @) w=-<(5 414+ D4D?+ Tc) TEOREMA © pregene) = ppel= i * prepares) Ejemplo 5 Resolver: (D? — 2D +5) y = e*cos(2z) Solucion 1 = pT e) 1 Ops i aD FH) Or) 1 eon tem eyo) vp ‘sen(2z) ‘TEOREMAS 1 oxy) per PO 1s He" a Fepy emer +0)= rH ‘sen(ox +6) 1 » Ey mer += Brgy lan +) le) =e pelle) . Fate) =z Ejemplo 6 (Dy Poy pu Foe) Resolver: (D4 — 10D? +9) y = cos(2z +3) ‘alueion 1 = erate Hy 8) L 1 w= pe preg mee t9) 30Célculo Avanzado I = = Reemplazamos P(D*)= D249 + P(-a@)=-P+9 A a=? = Eprom ts) Heioates => Pcslaae +t won a v5 epee w= peels +8) Ejemplo 7 Resolver: (D? +3D ~4) y=sen(2z) cién ‘a a (22) = payee) = “a os 4D con(22) (D+1)(D-4) "tee a) a eens) = iw [D? -3D — 4] sen(2z) we iz [-4sen(2z) — 6cos(22) — Asen(2z)] ee +3 [asen(22) +3 008(22)] = P(—D)sen(ax) Po sen(ee) = PaAPa) © sh senfaz) = P(-D)sen(ez) PO) ‘P(ai)P(-ai) P(—D)senh(ax) z Se RC conten meet P(—D) cosh(az) Sy cosh (az) = SYP (—a) Drape) = aap) 1 \ r * Grane) = gaan (end) 2.8. Conectores 1. Diferenciador a(t) a(t) 31Calcul Avanzado 1 = FIEE-UNI 2. Amplificador por una constante a(t) ce k-a(t) 3. Retardo a(t) oo a(t—t) 4. Integrador 26) y(t) = [ oo §. Sumador y(t) + x(t) = aQ+2Q) 6. Multip! licadorCélculo Avanzado I FIEE-UNI 7, Anversor a} 8. Diferenciador Multiple x Cart mavSOS at Hap _Kosita” Bjemplo 1 ‘En el diagrama mostrado determine la sefial de entrada y(t) J (A, a B EA a A b (PD) Bertsen(2) *y(t) Solucién Sy" +10 +y = te~* sen(t) hallando la solucién homogénes Soluciém homogénea waste Bet Abora hallando la solucién particular (9D? + 10D +1) y= te~* sen(t) 1 L iain ree) wae ape yp same" [F,(z)cos(wz) + Ge(z)sen(wx)] WO}=9 ¥O=1 — y"0)=2 _ Solucién 2 = / QR Solucién homogénea 7, & r(r=Dr+1)=0 © (r-1Nr+1) 3 > mal v MAF Bel + Ce-* Solucion Particular wf) ene + Be BStHES ns Analizando miembro a miembro L(omultiplicidad 1) Por el 2do caso (r Por el 2do caso (rt) + yp = Be Up = Dze“* + BE* Reemplazando yp en la cevacién diferencial Resolviendo D=2 A £ = 2x2 4 Le pa 2ee* +50 Sumando ya ¥ up ve) = A+ BE +e? +208 + wo)=0 + A+B+cti=o ¥O=1 3 BC+ vW@)=2 + B+o-2=2Céleulo Avanzado I FIEE-UNI —y" = 122? +62 Solucién Solucién homogénea Part=0 + r(r-1)=0 i 1 > oma 3 wal m=0 + mae? yay + pear + mel + ys un = A+ Bot Ce* Hallando la solucién particular up = f(2) Ya Dsi+ E+ F Por el ler caso (r = 0) w=2" (Dz? +Er+F) comom=2 tp =2 (Ds? + Bo +P) ‘Reemplazando yp en la ecuacién diferencial oft — of = 1208 462 Resolviendo A 8S A vp = 23 (a? 32-18) Sumando va ¥ ¥p ye) =A+ Bes Crea (—2 Ejemplo 3 Resolver: of! +y/ = xc00(2) Solucién Solucién homogénea re Hallando la solucién particular ip = Prey 0s(2) + Qa) sen(2) tp = ef (Pte) c08(z) + Qa) sen(z)] Dondet=Oyu=1 Por el Ato caso (r 0) Up = (Cz + D)cos(2) + (Ex + F) sen(z) Reemplazando en la ecuacién diferencial Gant 2 ( 1 w=Célculo Avanzado T FIBE-UNI Sumando yp ¥ Yp ead coda) (F245 ve) = A+ Be + at5 sen(z) Ejemplo 4 “Resolver: y/’ +4y = sen(22) Solucién Sclucién homogénea Pee 3 (rr +2) =0 m= a cos(2z) + yy = c0s(22) m=O sen(2z) + yp = sen(2z) a = Acos(2z) + Bsen(22) Hiallando le solucién particular Up = 2" [P(x cos(ur) + (a) sen(w2)] Donde t=Oy w= 2 Por el dto caso (r 74 024) A (m= 1) tp = 2 [Coos(22) + Dsen(2x)] Reemplazando en la ecuscién diferencial Uf + yp = sea(22) Donde ‘Sumando vp ¥ vp v(e) = Acos(2x) + Bsen(2x) ~ Fcos(2x) TEOREMA (principio de superposicién) Sime va, a son goluclones particulares._.,. v' + Pla) + Qe = Pilz) v' + P(a)/ + Q(e)y = Ra(2) Respectivamente entonces Yrnti Es une solucién de: vu" + P(z)y! + O(r)y = Pi(z) + Pale) Ejemplo 5 Resolver: of + 4y = deos(2z) + Beos(2) +827 — Az Solucién Solucién homogénea Peano 4 (r-2i)(r42)=0 ta = Acos(22) + Bsen(22) Hallando la solucién perticular Ry(z) = 4oos(22) +" +4y = 4 cos(22) 40Gélculo Avanzade I FIBE-UNI vy = ae (Plz) cos(wr) + Q(z) sen(wz)] Donde t=O yu =2 o£2) Afm=2) tpt = #[Ccos(2e) + Den(2)] Donde C=0yD=1 pt = 2 [sen(22)] Para Ra(z) = 6cox(2) + "+ 4y = 6cos(z) wv" + 4y = 6cos(z) vp = ae [P(z) cos(uz) + Q(z) sen(uz)] Donde t=Oyw=1 Por el Ato caso (r #024) y (m=0) pa = Ccos(usz) + D sen (uz) ‘Reemplazando en Ja ecuactén diferencia! pa + dypa = Beos(x) Donde C=2y D=0 ‘upa = 2eos(wer) Pare Re(z)= 827-42 + y!" + dy = 8x? — 4x Por el ler caso (r #0) ups = Cx? +Dz+E ‘Reemplazando s+ Aups = 80? — Ae G=2 A De-l A Bet pg = 2-2-1 Por el principio de superpasicién: ip = 2o0n(2a) + 2008(2) +227 — 2-1 Sumando ya ¥ Yip u(x) = Acos(2z) + B sen{2x) + xsen(2x) + 2cos(z) + 2x” — 2-1 es) Ls eousciéa vw" + P(z)y’ + Q@)y=0 Se puede transformar en ult f(zu=0 Haciendo y = u(z)-v(z) Donde v(z) se elije apropiadamente DEMOSTRACION: yeu vavotw! ofl suv + Due! tw" Reemplazando en la ecuacién wt (i +E p(aye + a) at (r@ ei ra) % aLQUABIO Rostra” tA CAleulo Avanzado I FIEE-UNI zy A P@)+2v=0 Blemplo 6 oz) = eB SPO Resolver: ay"+2y' +ay=0 Solucién Por el teorema Hallando f(z) ‘Segtin el teorema Solucién homogéaea como y= uv Ejemplo 7 Resolver: y"+ 2zy/ —2y Solucién Como y=u-¥ a FO) = Fo + Po +20) fa)=1 w+ fejus0 vw’ tu=0 Pe1=0 + (-dr+)=0 m= eFcos(2) + ucos(z) uz=esen(z) + wsen(2) ne on gett pa, Gl =0 eget @@a-% vam ert SMe tlie 4 gat e=-5 + w'-Funo suse a2CAiculo Avanzado I FIES-UNT 2.12. Método de Variacién de Parametros El método de variacién de pardmetros como sigue. ¥" + P(e)y + Qe) = f(z) (239) ‘Sabemos que le. correspondiente ecuacion homogénea vo" + Pet! + Qv = (2.40) tiene una solucién general yp (2) de la forma. ' {uw} w= An + Bue El método consiste en reemplasar A y B por fumciones u(a) y u(z) que se determinan de maners que la funcién resultante = UM + UME (2.41) Derivando (2.41) , se obtiene = (ofan + my] + [ohva + maul = fein + viel + (oars + oh ‘Con ests condicin, la expresién de yf, se reduce a Ia forma Uy = ry} + av (2.42) Derivando esta funcién, se obtiene Uy = ah + mia] + [aus + vay] (2.43) ‘Sustituyendo (2.41) , (2.42) y (2.48) en (2.99) y agrupando separadamente los términas en u y en v, se obtiene jdamente [eich + visi] + [vine + veut] + Peay rith + vat) + Qe) [tits + vara] = fl) v1 [A + Peay + Qu] + 00 [ud + Payne + Qeyua] + ink +h = F(2) Y como y; ¥ ye son soluciones de la ecuacién homogénea (2.40) , esta expresién se reduce Oviatt = ( vinh + teh = Fle) ) Este es un sistema de dos ecuaciones lineales algebraicas que tienen como incégnitas las funclones u’ y v’. La ~Solucié[Link] obtiene aplicando Is regia de Cramer __ Ss 3+ fa BE (2.44) donde W = nuh vive en el wronskiano de yy y ye. Es evidente que W # 0, porque las funciones y: y y2 constituyen un sistema fundamental. Integrando las ecuaciones (2.44) , se tiene =- [Ma om [YE una f Mas, [we YY ya que r(x) es continua, estas integrales existen. Al sustituir estas expresiones para uy v, en (2.41) , se obtiene Ja solucién deseada de (2.39) ane we, tole) = 1 if We + nf a (2.45) ‘Obeérvese que si en (2.45) se dejan como arbitrarias a las constantes de integracién, entonces (2.45) representa, una solucién general de (2.39) . 8Célculo Avanzado I ‘FIEE-UNI Ejemplo 1 Resolver: yf'-+y=cse(z) Solucién Soluciéa homogénea v+y=0 4 41-04 n=akt =H nsetsen(is) + m=sem(z) 3 va —e™cos(b) > yx = cox(e) ‘wh = Asen(z) + B cos(z) Solucton particular Wp = v1c00(z) + vasen(z) Bn el método ¥{ con(x) + obsen(z) = 0 snultiplicando por coe(z) vj sen(xz) + v} cos(z) = ese(z) amultiplicando por — sen() ye a nee n v= cot(z) + v2 = In(sen(z)) Luego tenemos ‘Up = ~z.c08(z) + In(sen(z)) -sen(z) ‘Sumando yp ¥ Yp y= Asen(z) + B cos(a) ~ 7 00s(2) + In(een(z)) - sen(2) 2.12.1. Ecuaciones con coeficientes variables y+ On-1(z)y"t + bale)? Solucién fundamental de Ia homogénea thy! +b =S(z) fonsvasva,--- tn} Solucién particular Yo = Cala)pn + Cale) ya + Cal@)ys +--+ Ca(Z)Yn Por variecién de pardmetros Cia) + Ch(z)va + Co@)ys +--+ Cy(z)yn = 0 Cie + Hele + CE) +--+ Che), = 0 Cllelaa* + Osea? + Cladus be CLa)aa” = 52) Ejemplo 2 Resolver: (cos(2) — sen(x))y/” + 2sen(x)y’ — (sen(x) + c0s(2))y = e*(cos() — sen(z))? Soluetén Sumando coeficientes y por inspeceién nse por D’Alembert ywaucje® + yp =sen(z) n= Act + Boen(z) Solucién particular tp = vi(2)e* + va(2) sen(z) Método de variacién de parimetzos Ula)e* + C%(2) seals) = OL (z)e® + O5(2)c00(2) =e (c00(2) — sen(2)) 4Célculo Avanzado T Gi(z) = cos(z) Orla) = e* Up = C cos(z)e* + De® sen(z) Sumando tn ¥ vp y= Ae + Boen(e) + Ccos(a)e™ + De* sen(a) Bjemplo 3 Resolver: ay" ~ y/ — 4a%y = 1625037 Solucién Solucién homogénea ns Ae + Be? Solucién particular : ‘ tp = Cilz)e* + Coe-* Por variacién de parémetros * 7 Ol(z)e" + Oye" =0 i(e) (200%) + 04(2) (—20e-*") = 16:4" Ci(w) = 227 A Ons) = 2226 — p= 2atet + (2n%e* — o**) =? Sumando yp ¥ Yy wade + Bet 42ztet + (ater et") et FIEE-UNICéleulo Avanzado I FIEE-UNI Problemas Resueltos 1, Determine el WRONSKIANO para la ecuacion diferencial ¥"(2) + p@)y" + ae)’ + ry = oz) Solucién ua = {ya(=), voz), wole)} up = {Ci(z)s(z), Cale)u2(2), Ca(z)us(2)} Uy = Cin + Cah + Chae + Cath + Che + Cath + Cin + Chia + Chys =0 uy = Cah + Cou + Cash (2.46) y= Cink + Cau + Cith + Call + Cig + Cath + Oy + Oh + Ch =0 Ys = Cru + Caf + Cou 47) Up! = Chul + Cayh! + Caaf + Cays' + Cau + Cavs (248) Reemplazando (2.46), (2.47) y (2.48) en Ie ecuacién diferencial Cul + Cand + Cah = az) Cin + Cam + Cis = 0 Clash + Ch + Cink = 0 Chl + Cu + Cah =0 De la matris por Cramer: rome 29 ev | 2, Resolver Ja ecuacién diferencial a) (cosz)y" + (6cosx — 2senz)y/ + (8 cosx — 6senz)y = e* cos 2x + e~** sen 2c Solucién = (@sa)y uf = (een) + (cosz)y/ ut = (sen) (cos2)y + (c0s2)/" ~(6enz)¥ ul = (cosz)y" — 2(sen:2)y" — (cosz)y/Calculo Avanzado I FIBE-UNI ul + Ou! + 9u = e* cos2x +e“ sen 2x P+6r¢9=0 + un = exe 8? + egre ory = «e2)] ip F (acos2s + 4een24) +e (-f sen) yao [rowan ge ( € My = nt ty un = 0° + cane“ + & (Bc00 27 + deen 22) — 8, Hallar x(t) en la ecuacin diferencial x(t) ply im Bi Dy Be Solueion Primero, hallamos la solueién homogénea “2,(t)” @2!Y +162!" + 82! — Ba! — ar = fe Usando operador D 8D‘ + 16D° +8D? - 3D-3=0 Polinomio caracteristico = Soar tems enemy (eSNG AB Bem Ripe 1 neh nage 7 Kul = Cit + cx 4081 (cco (FE) + cin ()) Para la solucién particular, usaremos el método de coeficientes indeterminados X,(t) = te-*(AB + Be +Ct+D) Para el edleulo de cosficientes hallames 2/,2" 2!" y 2!¥ X's et(At + BE + OF + Dt) +e-\(4At + SBP +208 + D) X" = et(At" + Bt + C# + Dt) + e-'(12At? + 6BE +20) XM" m= —3e-*(4t? + 38H? + 2Ct + D) + e-*(24.4t +68) X™ = e"t(124t? + 6Bt +20) + e-"(244) a7Céleulo Avanzado I FIEE-UNI Resolviendo ‘Luego, la soluefén general es: 2(0) = Cre + Cae? +6 (Cy 00s 4. Resélver Iatecuscién diferencial 2) oy" —22y" 42x — dy =a Soluctén ‘La ecuacién diferencial es de Euler y la solucin sugerida es z [D(D -1)(D-2) — D(D -1) +2 -2y, [D402 +5D-2y =e Solucién homogénes y Polinomio carseterietico adr? 457-25 (2-1)? n=2 m=, nm=1 La solucién homogénea es ua(t) = Cie™ + Cae! + Cate! = Soluei6n particular: Método de operasiores inversos aacst=h() [aa VQ = unt ve Entonces la soluci6n general es a(t) = Cr+ (Ca — Det + (Cs + ett 8) oy! Ay + 4y = sen3e Solueién ~ Solucién homogénea y polinomio caracteristico y'-4/ +4y=0 Paat4=0 + (F-2)7=G n= a(t) = Cre + Caze™* ~ Solueién particular por método de cocficientes indeterminados ip: cen 3c 48CAleulo Avanzado I FIEE-UNI tp =e (Deos 32 + Esen32) Para el cfleulo de coeficientes hallamos yy ¥”” 12 (D ccs 2x + Bsen 2x) + o*(— 1e"* (D cos Su + E'sen3z) + 2e**(—3D sen Sa + 3E cos $x) + 2e**(-3D sen Sz +35 cos3z) +€*(=9D cos 3x — 9E sen 3z) Up(2) = —5e™*D cos3z — 12c™ sen 3x — 5e** Dsen 3x + 12e™*E cos3z . ‘Reemplazamos en la ecusciGn diferencial 82-4 9B eos22) —9e™ Deos Sx — 9e™* B sen 3z = e™ sen Sz 08 =14E=1/9, D=0 isle) = — Je ends ‘La solucién general es ult) = Oe + Cane — 2 sends 5, Resolver la ecuacién diferencial [(2—1)'D* + 2(32 — 1)8D* — 20(82 — 14D +-60(82 — DN] y= Se ED sent (ine ~1)) ividimos entre (2s — 1) [(e2—1)°D* +2(82—1)2D* — 2o(ae—1)D + 60]y = 2D ae en? (in(2a—1)) hhacemos t= In(3x— 1) ‘entonces [8:D(D —1)( — 2) +2-3D(D —1) — 20-20 +60] 2) = [= oat] [27D8 - 63D? — 24 + 60] x(t) = (1+ 2e~*) sen?(¢) Para el célculo de an(t) [a7D* — 63D? — 24 + 60] a(t) =0 eee 2048 6379 = 24r ¢60= n=1, n= a+) ae a v6) 2n(t) = Cre? + Ged HV + Cgedll-VOl= * Para hallar £,(¢) (B= nperD(D - 2) + 18D — eg2(0) = 7) sent) co0s(2t) a te 8(D — 1)(0D? — 12D — 20)z(¢) = E is *eteo(n] 10 «gp =x [eorctocm {f- ts 9}] 1 1 1s a) =— mapa) - IND =a OO + BDH 1 + — ID= HCI — aay He) 49Caleulo Avanzado I FIEE-UNL (8) = sg + peag Ht sen 2) + 52c0s(22) — ras =a +H) 6. Bn el esquema mostrado, determine y(t) , 20(t). Si P(t) = hie cos(t) , (0) = 22(0)=1 , 24(0) = 240) =2 k=3m=c (numéricamente). ‘ & } 3 m m LK FO a(t) A(t) Solucién Diagrama de eada bloque _#0 cy Mm L (x1 — 2) Hos (t)~ mlt)) eX) = mo may (t) + ci) + A(er(t) — 22(t)) =0 (2.49) _ 0) —k(a2 — 21) m ext) ma4{0) = —K(ca() ~ 23(0)) - 20. “mers(t) + ox (t) + k(za(t) — a1(¢)) = (2.50) (2.49) + (2.50) Hecemos &,(t) + za(t) =2(t) 40+ 240 = 20 +240 =" entonces m(z{(0) + a3(O) + efi (t) + 24()) = 0 m(2"(f)) +e2'(t) =0 2") +32!() = ecuacién caracteristiea 43r=0CAlculo Avanzado I FIEE-UNT por lo cusl se obtiene 2() = C1 +Cre* o Jee dates del pro! 2(0) = 21(0) + 22(0) 2=4 cAlculo de 2!(t) en t= 20) centonces se tiene a+ = Baa despejamos xo(t) en funcién de 21(¢) a(t) = 2 — fer a(t) (2.51) matey reai( +8 [rx - (F—Fer*—a9)] =0 maf{(t) +exh(t) +k [pox -2 a =0 mati teri) +2 (9 = [2 - Bate) + Sait tanto = 2-4 FtO+e1 +29 = 2-fem AA(O) + 3a (t) + 6xr(t) = 10— de“ + Pora la solucién homogénes : n= + |au)= Cre eos (4 + Para la solucion partionlar : [DP +3D+6] 21,() =10-4e-* i 2.0= pape usando los primeros términos del resto de la divisién L D?+3D+6 6 12 1. L 3) + syl= (3-8) 01- perapes tel 10 4 ge n,0-B- te 51CAlculo Avanzado I FIEE-UNI 2()=1 cBlculo de 24 (¢) Hl) =—oreH cos H()=2 nl) = De @51) 7 10 — ge" — aa(t) = erm 52Calculo Avanzado I FIBE-UNI Problemas Propuestos 1, Resolver 4a) (senz)(2’)® — 8(cos2)2'2” — (senz)z” — 2’senz=0 8) (ey? + day? + y)de + (oP? + 20%y + 2)dy =0 2, Resolver wy" + doy + (2? +2)y=1 3. Resolver (D° = 3D* +4? — 40? + 3D —1)y = 26 + 3senz +1 4. Determine la soluciéa de la ED, usendo D'alambers, 2y" —8ny/ — (4 824)y = 0 5. Resolver (tan*r)y/" + [(1 + sec} tan? zp + (sec? a)y = (1 + esc2) sec? 2 6. Resolver yO +8 4 16y = (2 — 1) cen 2x + e-* cos 2e 1., Resolver (cos z)2/" ~ 3(senz)2'2" — (cosz)z* — 9(cos2)2" + 8(senx)(z!)? + 3(cosz)2! = tet sent 8. Resolver 4: ay'42aty +y=2tt ‘9. Determine le solucién de (tsent)z! + 2(sen £ + toast + tseni)z! — 6(sent + cost + toost)x 10, Resolver les ecuaciones diferenciales 2) yy" = (ya? 8) 2(e—1)y" + (22° — 20? — a) <2 *e- 1)? Li, Resolver ‘9a? [(e— 2)*D* + 2(82 — 2)°D? — 20(82 — 2)*D + 1802 — 120] y= ( ) sen? (In(3a — 2) 3-2 12, En el circuito mostracdo, las corrientes y voltajes iniciales son nulos. Determine el voltaje instantaneo V(t). Si BQ) —tetsen(t) , M=20, Mall, Ra, Lait , Cute 13. En el cireuito mostrado, el interruptor “s” se cierra en ¢ = 0 y las corrientes iniciales son nulas. Si Ew iye ipa = 10? , B(e) = 6sen(3004). Hallar : aeCAleulo Avanzado 1 FIEE-UNI ) Las ecuaciones diferenciales lineales de 2do orden para i(t) ¢ éa(t). ») fA(€)in(t) y V(t) em el condensador. Ly £852 72 a + BW) G 1d. En el esquema mostredo, caleular 2(¢), si 2(0) = 2'(0) = 0 , F(t) = 100|sen(60)|. m=tkg , k=10% , c=20 m 15. Determine 2(t) en el diagrama adjunto. N~seg k L_, x(t) mp— Fo 2) 16. En el cireuito mostrado, los — Ose D3 — (4042)? (241) valores iniciales de corrientes y voltajes son mulos. Calcular la corriente instanténea i(t) a través de Is bobina. vite) 17. Resolver Lmigir , Omir , Bm 1m0eteen() (-1)? +1], =sen( ma 4 ) EyCAlculo Avanzado I FIEE-UNI 18. En el circuito mostrado, los condensadores estan inicialmente cargedos con un potencia Vo y los interrup- tores estan cerradas. Calcule los potenciales en los condensedores para i> 0. ¢ 19. Resolver ay! — 2oy' + (2—922)y = 2°(cosh(3z)Capitulo 3 Ecuaciones en diferencias finitas ‘Son ecuaciones de la forma (GnB" + dy1E™? + n_2E"~? + +--+ 028? + 0B" + a0) y = f(z) donde E + operador lineal de desplazamiento /E" f(z) = f(@-+n) a0} GaiX0 § y=ul2) L[E| = 0nE” + Ona E"~? + Ong B"2 ++---pagE* + a,B' +09 + LIE] = f(z) Ecuacién diferencial con diferencias finitas homogéneas :{(z) = 0, entonces (a, B" + ana 8"! + ong B™? +--+ + anB? + a1B" +09) y=0 (8.1) Solucién de una ecuacién diferencial: La solucién homogénea es de le forma y=r#, eR Byatt Bfy=1# a Ba (3.1) £1* [ogr" agent! +oqaat™2 4+ -b at? 4ayr? + 09] = 0, donde 7? £0 Polinomio Caracteristico cat” HOPE agar Et Ont Pott Donde hay n rafces ri,r2,rs,...Tn ¥ son linealmente independientes. Y= Curt + Card + Car§ bo + Care Solucién general yy, = Cir +CarZ Eeuaciones diferenciales con diferencias finitas de segundo orden Polinomio caracteristico: agr? + ayr + a9 = 0 Hay tres soluciones 1. Raices reales no repetidas 2, Raices reales repetidas m1 = ra Us 3. Raices complejas conjugadas r1=a+i8 ; ra=a un = Cr(a+48)* + Cola - 18)" 56CAlculo Avanzado I FIBE-UNI si AVR; oan (2) un = Ci(Re®)” + Ca(Re”)* ‘Un = OLR® [o0s (ix) + sen (iz)] + C2 R® [cos (iz) — sen (éz)] un = CiR* [A cos (9) + Been (0z)] Ejemplo 1 Con le condiciém f :<0,-+00 >—+ R hallar f(z) = 2 Solucién Anilisis discreto: fiat +R i wna L Li an {hae} Ejemplo 2 Con Ia condicién f+ R hallar f(z) =a Soluetén as fz uO +R neon? ni? = 0,1,4,9,.. 3.1. Operador de Diferencia Finita Lineal Af@) = f(z +1) - f(z) analizando las distintas diferencias y derivando Ay= vast Ua A? yn = B (Oya) = A (nti — Ye) = Atingt Aun Dyn = Unt2 — 2nga +90 ‘Tercera diferencia, M8 yn =A (07yn) A (Ynt2 = naa + Ye) Btn = Yost — Wns + Stet — Yo ky = A(A‘ yn) Sunt — tai = 4 nts —AYne2 +Yne — Bye =n? Orden de uns ecuscién: Es la diferencia entre el sub indice mayor y el sub {ndice menor K-ésima diferencia Uinta + Synth — Us ‘orden orden 2 Sunset tines = 2n orden orden L 37FIEE-UNI Un = Unt? — Wns + Ya, Yosa = (n+ 2)7-+2(n-42) =n? + 6048 Yat = (a+1)? +2(n-$1) =n? $4043 A Un = Unt ~ nat + Yn A? yn =2 Ejemplo 2 Si ya =e hallar A?y, Solucién Yai tein = 27 yn = Und — ais + Yn Ay, =e" (Pe 1) Ejemplo 3 Demuestre que tix = 2" —n os una solucién de la ecuacién: Yne2 — Bnet + 2y Soluciéa Yata = 24? — (n +2) ders = 24? — (n+ 1) Reemplazando en la ecuacién for? — (n+2)] — 3 [24 - (n+ 1] +2[2" mJ =. Demostrado que cumple con la coudicién Yat? = Syntt + 29m 3.2. Ecuaciones en Diferencias finitas de Primer Orden ayn toon 6 Va #0 tnt = Ata +B ‘Obtenemos una solucién a partir de la ecuseién anterior n= A+B = A(AptB)+B = My t+ AB+B w=A(42%y + AB+ B)+B= Ay +A°B+AB+B Ya = ATyot (Lt A+ A? +A +. FA) B 58