El libro interactivo de “Ecuaciones diferenciales” se pensó para que el estudiante tenga una herramienta práctica en el estudio de las ecuaciones sin perder el rigor que ésta conlleva, el libro contiene teoría, ejemplos resueltos paso a paso y ejercicios propuestos con su respectiva respuesta, pero también, tiene escenas interactivas donde el lector puede cambiar los enunciados y visualizar las diferentes gráficas e interactuar con ellas, también cuenta con escenas para que el estudiante repase conceptos de asignaturas previas como matemáticas operativas, cálculo diferencial, cálculo integral y álgebra lineal. Recuerde que para tener un buen desempeño en esta asignatura se debe manejar muy bien los conceptos previos adquiridos en las asignaturas mencionadas anteriormente.
El texto incluye 40 escenas interactivas, diseñadas en DescartesJS y GeoGebra, para que los estudiantes y el lector puedan tener un acercamiento más práctico con las ecuaciones, también se incluyen 26 vídeos de YouTube creados por el mismo autor del libro para facilitar el desarrollo y mejor entendimiento de los diferentes contenidos.
El libro tiene cinco capítulos, introducción a las ecuaciones diferenciales, ecuaciones diferenciales de primer orden, ecuaciones diferenciales de orden superior, transformada de Laplace y series de Fourier. Con esto se puede comenzar aprendiendo a clasificar las ecuaciones para luego resolverlas de acuerdo a las características de cada una de ellas.
Para evitar, al máximo, la dependencia con la conectividad en la red, las expresiones matemáticas se han construido recurriendo al API de KATEX.
9
Capítulo I
INTRODUCCIÓN A LAS ECUACIONES DIFERENCIALES
1.1 Ecuaciones Diferenciales y modelos matemáticos
Las diferentes leyes que rigen el universo están descritas en lenguaje matemático. Cuando se deben describir fenómenos estáticos, el álgebra es suficiente para resolverlos, pero la mayoría de dichos fenómenos naturales involucran cambios que dependen del tiempo, es decir, implican una razón de cambio. Cuando un proceso describe el comportamiento de algunos fenómenos de la vida real en términos matemáticos recibe el nombre de modelo matemático y se utiliza para comprender las variables que lo conforman
Luego de plantear un modelo, es necesario resolverlo y confrontar la solución para saber si es razonable con los datos experimentales o con los hechos conocidos. Las hipótesis que se plantean implican una razón de cambio de una o más variables independientes con respecto a una o más variables dependientes, dicho modelo matemático implica una ecuación diferencial (ED)
Las ED tienen una importancia fundamental en las matemáticas debido a que son aplicables a casi todas las ramas del saber, como: geometría, física, química y biología. Éstas permiten identificar y aplicar técnicas y conceptos matemáticos que admiten obtener e interpretar la solución de una situación problema en la que intervienen relaciones entre tasas o razones de cambio de variables definidas en diversos contextos.
Algunas aplicaciones son: movimiento rectilíneo y caída libre, deformación de resortes y vigas, mezclas de soluciones salinas, circuitos en serie, desintegración radiactiva, fechado de fósiles con carbono 14, crecimiento y decrecimiento poblacional, propagación de virus y enfermedades, cambio en las temperaturas y reacciones químicas entre otros.
13
1.2 Definición de ecuación Diferencial
Definición de ecuación diferencial: una ecuación que contiene derivadas de una o más variables respecto a una o más variables independientes, se dice que es una ecuación diferencial.
Hay dos formas de expresar ecuaciones diferenciales, las dos se van a utilizar en este libro.
a) dxdy+4y=1
b) y′′′−5y′′+4y=0
NOTACIÓN:
Notación de Leibniz: Utiliza los símbolos dy∧dx para representar incrementos infinitamente pequeños. Se representa como sigue:
dxdy,dx2d2y,dx3d3y…dxndny
Notación de Primas: La notación prima se usa para denotar el número de la derivada siendo esta usada hasta tres prima (′′′), de ahí en adelante se utiliza notación numérica: la cuarta derivada se denota y(4). La n-ésima derivada se escribe como y(n). Se representa como sigue:
y′,y′′y′′′y(4)…y(n)
14
la notación de Leibniz llamada así en honor del filósofo y matemático alemán del siglo XVII Gottfried Wilhelm Leibniz, utiliza los símbolos dx y dy para representar incrementos infinitamente pequeños de x∧y respectivamente, al igual que Δx∧Δy representan incrementos finitos x∧y.
Leibniz comenzó a usar el carácter ∫ basándose en el carácter de la palabra latina summa, que escribió como una "s" alargada. Este uso apareció por primera vez públicamente en su artículo "De Geometría", publicado en Acta Eruditorum de junio de 1686
La prima es un signo gráfico que se usa en matemáticas, artes, ciencias y unidades de medida. En trigonometría, la prima sencilla indica la precisión de minutos y la doble los segundos al referirse a medidas angulares, cuando se aplica la nomenclatura sexagesimal. Por ejemplo: 12018′5′′ denota 12 grados, 18 minutos y 5 segundos.
En el caso de las funciones matemáticas la prima sencilla indica la primera derivada de la función, la doble prima indica consecuentemente la segunda derivada y así sucesivamente, de doble prima en adelante son derivadas de orden superior.
En este libro vamos a utilizar las dos notaciones, pero es importante resaltar que aunque la notación de primas es más práctica y rápida de escribir, la notación de Leibniz tiene una ventaja sobre la notación de primas y es que muestra claramente ambas variables, la dependiente y la independiente.
15
1.3 Clasificación de las Ecuaciones Diferenciales
Se clasifican de acuerdo a tres propiedades.
1.3.1 Según el tipo:
Ecuación Diferencial Ordinaria (EDO): es la ecuación que contiene derivadas respecto a una sola variable independiente. (Son las que se estudiarán en este libro).
dxdy+y=senx
y′′′+5y′−3y=0
"x" Variable independiente.
"y" Variable dependiente
Ecuación Diferencial Parcial (EDP): es la que contiene las derivadas parciales respecto a dos o más variables independientes.
∂y∂u=∂x∂v
"x" ∧ "y" Variables independientes.
"u" ∧ "v" Variables dependientes
Nota:Como todas las ecuaciones diferenciales que vamos a trabajar en este libro son ordinarias las vamos a llamar ED
16
1.3.2 Según el orden:
El orden de una ED es el de la derivada de mayor orden que hay en la ecuación.
dx2d2y+5y=senx ED de segundo orden
y′′′+5y′′−3y′+y=0 ED de tercer orden
dx4d4y+dx2d2y+y=ex ED de cuarto orden
y′+2y=tanx ED de primer orden
1.3.3 Según la linealidad:
Se dice que una ED de n-ésimo orden es lineal si F es lineal en la variable dependiente, esto significa que una ED es lineal si tiene la forma
- la variable dependiente ∧ todas sus derivadas son de primer grado.
- Los coeficientes de an,an−1,...,a1,a0∧ las funciones dependen a los más de la variable independiente.
xy′′−4y′=x2 ED lineal de segundo orden
17
La ecuación es lineal porque la variable que tiene exponente es la variable independiente.
x3y′′′+xy′−y=ln∣y∣ ED no lineal de tercer orden
La ecuación no es lineal porque la función logaritmo depende de la variable dependiente
dx4d4y+17(dx2d2y)2+16y=0 ED no lineal de cuarto orden
La ecuación no es lineal porque la segunda derivada está elevada a un exponente
y′=y−y2 ED no lineal de primer orden
La ecuación no es lineal porque la variable dependiente está elevada a un exponente
En general es más fácil identificar cuándo una ED es no lineal, que cuándo es lineal.
Practiquemos
En la siguiente escena interactiva, adaptada de Plantillas con descartes_JS, debes arrastrar la flecha desde la ED lado izquierdo hasta el texto del lado derecho para identificar el orden y la linealidad de cada ecuación diferencial. Es importante aprender a identificar las ED según el orden y la linealidad, para aplicar el método adecuado de solución que veremos más adelante.
18
1.4 Solución de una ecuación Diferencial
Definición: Cualquier función f definida en un intervalo I∧ que tiene al menos n derivadas continuas en f, las cuales cuando se sustituyen en una ED ordinaria de n-ésimo orden reducen la ecuación a una identidad, se dice que es una solución de la ecuación en el intervalo.
Intervalo de definición: La solución de una ecuación diferencial debe estar definida en un intervalo. Dicho intervalo se conoce como intervalo de definición, intervalo de existencia o dominio de la solución.
19
1.4.1 Tipos de soluciones
Solución trivial: Es la solución de una ED que es idéntica a cero en un intervalo I
y=0
Puede ocurrir que una ecuación de orden superior tenga una o varias soluciones triviales, lo que no puede ocurrir es que todas las soluciones sean triviales
Solución explícita: Es la solución en la que la variable dependiente se expresa sólo en términos de la variable independiente ∧ las constantes.
La función y=4+21e−5x es solución explícita de la ED
dxdy+5y=20
En el ejemplo (1) de la página 21 se verificará que y=4+21e−5x es solución de la ED dxdy+5y=20
Solución implícita: Se dice que una función G(x,y)=0 es una solución implícita de una ED en un intervalo I, suponiendo que existe al menos una función f que satisface la relación así como la ED en I
La función x2+y2=4 es solución implícita de la ED dxdy=−yx en el intervalo −2<x<2
20
Ejemplo 1. Verificación de una solución
Verifique que la función y=4+ce−5x es una solución de la ED dxdy+5y=20
Solución
Como la ecuación es de primer orden, se debe derivar una vez la función.
dxdy=−5ce−5x
Reemplazando la función ∧ la derivada en la ED se tiene
−5ce−5x+5(4+ce−5x)=20
Aplicando propiedad distributiva
−5ce−5x+20+5ce−5x=20
Agrupando términos semejantes queda 20=20
Al demostrar que el lado derecho de la ecuación es igual al lado izquierdo, se verifica que la función
y=4+ce−5x es solución de la ED dxdy+5y=20
21
Ejemplo 2. Verificación de una solución
Verifique que la función y=xe−3x es una solución de la ED y′′+6y′+9y=0
Solución
Como la ecuación es de segundo orden, se debe derivar dos veces la función, tenga en cuenta que se debe derivar como un producto.
y=xe−3x
y′=−3xe−3x+e−3x
y′′=9xe−3x−3e−3x−3e−3x
Agrupando términos semejantes
y′′=9xe−3x−6e−3x
Ahora reemplazamos la función y las dos derivadas en la ED
9xe−3x−6e−3x+6(−3xe−3x+e−3x)+9(xe−3x)=0
Aplicando propiedad distributiva
9xe−3x−6e−3x−18xe−3x+6e−3x+9xe−3x=0
Agrupando términos semejantes
22
0=0
Al demostrar que el lado derecho de la ecuación es igual al lado izquierdo, se verifica que la función
y=xe−3x es solución de la ED y′′+6y′+9y=0
Ejemplo 3. Comprobación de una solución
Verificar si la función y=xsen(ln∣x∣) es solución explícita de la ED x2y′′+xy′+2y=0
Solución
Como la ecuación es de segundo orden, se debe derivar dos veces la función, tenga en cuenta que se debe derivar como un producto.
y=xsen(ln∣x∣)
y′=xcos(ln∣x∣)x1+sen(ln∣x∣)(1)
y′=cos(ln∣x∣)+sen(ln∣x∣)
y′′=−sen(ln∣x∣)x1+cos(ln∣x∣)x1
23
Ahora se reemplaza la función y sus derivadas en la ED de segundo orden.
Como no se cumple la identidad, se puede afirmar que
y=xsen(ln∣x∣)No es la solución de la ED x2y′′+xy′+2y=0
Ejemplo 4. Comprobación de una solución implícita.
Demuestre que la relación x2+y2=4 es solución implícita de la ED dxdy=−yx en el intervalo (−2,2)
24
Solución
Se deriva dos veces la función
2x+2ydxdy=0 despejando se obtiene dxdy=−yx
Además, resolviendo x2+y2=4 para y en términos de x se obtiene.
y=±4−x2
las dos funciones y=4−x2∧y=−4−x2 satisfacen la relación x2+y2=4∧ son las soluciones explícitas definidas en el intervalo (−2,2)
Las dos funciones, tanto la positiva como la negativa, satisfacen la relación ∧ son la solución implícita, definidas en el intervalo abierto (−2,2)
Ejemplo 5. Comprobación de una solución implícita.
Demuestre que la función x2y−tanx+y2=0 es
solución implícita de la ED (2xy−sec2x)dx+(x2+2y)dy=0
25
Solución
Primero debemos organizar la ED, dividiéndola toda por dx
(2xy−sec2x)+(x2+2y)dxdy=0
Ahora despejamos el dxdy obtenemos
dxdy=x2+2ysec2x−2xy
Ahora vamos a derivar la función x2y−tanx+y2=0 aplicando derivación implícita.
x2dxdy+2xy−sec2x+2ydxdy=0
Agrupando los términos con dxdy al lado izquierdo y los otros términos al lado derecho, obtenemos
x2dxdy+2ydxdy=−2xy+sec2x
Sacando factor común dxdy
dxdy(x2+2y)=−2xy+sec2x
Despejando dxdy obtenemos
26
dxdy=x2+2ysec2x−2xy
Como los dos resultados son iguales, podemos decir que la función x2y−tanx+y2=0 es solución implícita de la ED
(2xy−sec2x)dx+(x2+2y)dy=0
1.4.2 Familias de soluciones
Familias de soluciones: una ED puede tener una cantidad α de soluciones que corresponden a las elecciones ilimitadas de parámetros. Todas estas familias se llaman familias paramétricas de soluciones. El parámetro está directamente relacionado con el orden de la ED.
Si la ED es de primer orden, la solución es uniparamétrica. Hay un parámetro c que puede tomar diferentes valores.
Si la ED es de segundo orden, la solución es biparamétrica. Hay dos parámetros c1∧c2 que pueden tomar diferentes valores.
Si la ED es de tercer orden, la solución es triparamétrica. Hay tres parámetros c1, c2∧c3 que pueden tomar diferentes valores.
Si la ED es de orden n, la solución es n-paramétrica. Hay n parámetros c1, c2, c3…, cn que pueden tomar diferentes valores.
27
Solución general: Es la función que satisface a la ED ∧ que contiene una o más constantes arbitrarias dependiendo del orden de la ecuación diferencial y que fueron obtenidas de las sucesivas integraciones.
La función y=c1xcos(ln∣x∣)+c2xsen(ln∣x∣)
Es una solución general de la ED
x2y′′−xy′+2y=0
La función y=c1e2x+c2xe2x
Es una solución general de la ED
y′′−4y′+4y=0
Solución particular: Es la función que satisface a la ED ∧ cuyas constantes arbitrarias toman un valor específico.
La función y=52xcos(ln∣x∣)−31xsen(ln∣x∣)
Es una solución particular de la ED
x2y′′−xy′+2y=0
La función y=53e2x−21xe2x
Es una solución particular de la ED
y′′−4y′+4y=0
28
En el siguiente video, realizado por Jaime H. Ramírez Rios, podrás observar un ejemplo de comprobación de una solución explícita.
En el siguiente video, realizado por Jaime H. Ramírez Rios, podrás observar un ejemplo de comprobación de una solución explícita.
29
1.4.3 Curva solución
Curva solución: La gráfica de la solución de una ED se llama curva solución. Puesto que la función es derivable y continua en el intervalo de definición I. Es probable que la gráfica de la ED sea diferente a la gráfica de la solución, por lo cual el dominio de cada una no necesariamente debe ser igual.
Geométricamente, la solución general representa una familia de curvas. La solución general es una familia de gráficas que van cambiando dependiendo del valor del parámetro c, la solución particular es una de las curvas de la familia.
Geométricamente se puede hallar la solución de una ED y′=f(x,y) dada. En cada punto (x,y) del plano xy, el valor de f(x,y) determina una pendiente m=f(x,y), esto significa que la solución de una ED es simplemente una función derivable cuya gráfica y=y(x) tiene su “pendiente correcta” en cada punto (x,y(x)).
Una curva solución de la ED y′=f(x,y) es claramente una curva en el plano xy cuya línea tangente en cada punto (x,y) tiene pendiente m=f(x,y)
En la siguiente escena interactiva podrás interactuar con la solución y=ln∣∣t2+1∣∣+c de la ecuación diferencial dtdy=t2+12t.
En la siguiente escena interactiva, diseñada por Jorge L. Herrera, podrás ver cómo cambia la gráfica de una ED al cambiar el valor del parámetro c.
30
En la siguiente escena interactiva, diseñada por Adrián Rocha Ichante, podrás ver cómo cambia la gráfica de la ED al cambiar el valor del parámetro c. La solución es una familia de soluciones.
31
1.5 Problemas de valor inicial
A menudo es interesante resolver una ecuación diferencial de primer orden de la forma dxdy=f(x,y), sujeta a la condición inicial y(x0)=y0
Ejemplo 1. PVI de primer orden
y=cex es una familia uniparamétrica de soluciones de la ED de primer orden y′=y en el intervalo (−∞,∞). Determine una solución del problema de valor inicial PVI de primer orden que consiste en esta ecuación diferencial ∧ la condición inícial y(0)=5
Solución
y=cex
Se reemplaza la condición inicial y(0)=5. Recuerde que la condición inicial es y(x0)=y0, es decir, y=5 cuando x=0
5=ce0
El valor del parámetro es:
c=5
reemplazando en la solución general, obtenemos la solución particular.
32
y=5ex
Ejemplo 2. PVI de segundo orden
y=c1e−x+c2e2x es una familia biparamétrica de soluciones de la ED de segundo orden y′′−y′−2y=0
Determinar c1∧c2 de modo que se cumplan las condiciones iníciales y(0)=2∧y′(0)=−3
Solución
y=c1e−x+c2e2x es la solución general de la ED y′′−y′−2y=0
Para hallar los valores de c1∧c2 y expresar la solución particular se debe reemplazar la condición y(0)=2 en la solución de la siguiente forma:
2=c1e0+c2e0. Recuerde que e0=1, entonces queda una ecuación con dos incógnitas 2=c1+c2 y la llamamos ecuación (1)
Para resolver esta ecuación se necesita tener otra ecuación con dos incógnitas, entonces se deriva la solución y se reemplaza la condición y′(0)=−3 de la siguiente forma:
Derivando la solución obtenemos y′=−c1e−x+2c2e2x
33
Reemplazando la condición inicial −3=−c1e0+2c2e0 se llega a la ecuación −3=−c1+2c2 y la llamamos ecuación (2)
Con las ecuaciones (1)∧(2) se puede formar un sistema de ecuaciones
{c1+c2=2−c1+2c2=−3
Este sistema se puede resolver por igualación, sustitución o eliminación.
Teniendo en cuenta las dos ecuaciones resultantes, lo vamos a resolver por eliminación
c1+c2=2
−c1+2c2=−3
3c2=−1
Obtenemos c2=−31
Reemplazando en la ecuación 2=c1+c2 llegamos a
2=c1−31
Despejando c1 se obtiene
c1=37
34
Entonces la solución particular es
y=37e−x−31e2x
Ejemplo 3. Intervalo de definición de una solución
y=x2+c1 es una familia uniparamétrica de soluciones de la ecuación diferencial de Primer orden y′+2xy2=0. Determine una solución del PVI y(1)=4 de primer orden que consiste en esta ecuación diferencial y las condiciones iníciales dadas. Dé el intervalo en el cual está definida la solución.
Solución
y=x2+c1 es la solución general de la ED y′+2xy2=0
Para hallar el valor de c se debe reemplazar la condición y(1)=4 en la solución de la siguiente forma:
4=12+c1. Despejando obtenemos. 4+4c=1
Despejando el parámetro c obtenemos c=−43
35
La solución particular es y=x2−431 el dominio es R−{±23}
El intervalo de definición está dado por.
(−∞12−3)U(2−3,23)∪(23,∞)
1.5.1 Campos direccionales
A partir del capítulo 2 se hará un análisis de los diferentes métodos para resolver ED, para lo cual, es necesario identificar primero el tipo de ecuación para saber cuál o cuáles métodos aplicar, pero antes de proceder a dicho análisis, nos detendremos en la interpretación geométrica de las ED y sus soluciones.
Desde el punto de vista geométrico hay un método para obtener soluciones aproximadas de la ED y′=f(x,) a través de cada grupo representativo de puntos (x,y) en el plano se obtiene un segmento lineal con pendiente m=f(x,y), todos los segmentos lineales constituyen un campo de pendientes o un campo direccional
La terna (x,y,y′) determina la dirección de una recta que pasa por el punto (x,y). El conjunto de los segmentos de estas rectas es la representación geométrica del campo direccional
Una Isóclina es el conjunto de los puntos del plano en donde las rectas tangentes a las gráficas de las soluciones de la ecuación diferencial tienen la misma pendiente. Estos puntos son aquellos que satisfacen y′=f(x,y)=c
36
En las siguientes tres escenas interactivas podres ver el comportamiento de los campos direccionales y las Isóclinas de diferentes ED.
En la primera escena interactiva podras ver el comportamiento del campo de direcciones de una ED de primer orden dxdy=−yx y la gráfica solución para diferentes valores del parámetro (radio de circunferencia).
En la segunda escena interactiva podras ver el comportamiento del campo de direcciones, densidad y longitud de una ED de primer orden dxdy=x, pero puedes cambiar la ED ingresando otra función diferente, además, puedes visualizar soluciones particulares para diferentes valores del parámetro.
En la tercera escena interactiva podrás interactuar con las Isóclinas que están formadas por la solución de la ED
y′=x+y−1.
Al mover el deslizador, podrás ver las pendientes de las rectas tangentes de acuerdo con la recta
k=x+y+1.
Las flechas indican cómo va cambiando el campo de direcciones de la ED
En la siguiente escena interactiva, diseñada por Jorge Olivares Funes, podrás ver el campo direccional de una ED y cómo cambia el radio de la circunferencia.
37
En la siguiente escena interactiva, diseñada por Soed Torres, podrás ver el campo direccional de una ED.
38
En la siguiente escena interactiva, diseñada por Diego Sandoval, podrás ver las Isóclinas de una ED y cómo cambian a medida que se cambia el valor del parámetro.
1.5.2 Existencia y unicidad
Cuando se considera un problema de valor inicial (PVI) se deben considerar dos preguntas:
¿Existe la solución del problema? Si la respuesta es afirmativa se debe preguntar ¿La solución es única?
Dada una ED dxdy=f(x,y)
Donde f(x,y) está definida en una región rectangular R que contiene al punto (x0,y0)
39
Si f(x,y) satisface las condiciones:
1.f(x,y) es continua en R,
2.∂y∂f es continua en R
Entonces existe un intervalo I con centro en x0∧ existe una y sólo una función y=g(x) definida en el intervalo I que satisface la condición inicial y(x0)=y0
En otras palabras, las condiciones para la existencia de la solución son:
• Continuidad de f(x,y) en R,
• Acotamiento de f(x,y) por R,
Y las condiciones para la unicidad son:
• Continuidad de f(x,y)∧∂y∂f en R,
• Acotamiento de f(x,y)∧∂y∂f por R,
Estas condiciones son suficientes pero no necesarias, porque puede existir una solución que satisface y(x0)=y0, pero que no cumple la condición 1, o la condición 2, o ninguna de las dos.
A continuación vamos a ver un ejemplo y luego a definir la existencia y unicidad de soluciones.
40
Ejemplo 1. Existencia y unicidad
Cada una de las funciones y=0∧y=x2 satisface la ED dxdy=2y∧ la condición inicial y(0)=0
Solución
El PVI dxdy=2y, y(0)=0 tiene al menos dos soluciones. como se muestra en la figura 1.
El estudiante puede verificar que las dos funciones son solución de la ED, además, las gráficas de las dos soluciones pasan por el mismo punto (0,0)
Gráfica 1 Dos soluciones del mismo PVI.
41
Existencia y unicidad de soluciones: Suponga que tanto la función f(x,y)∧ su derivada Dy∧(x,y) son continuas en algún rectángulo R en el plano xy que contiene el punto (a,b) en su interior. Entonces, para algún intervalo abierto I conteniendo el punto a, el problema de valor inicial
dxdy=f(x,y)
Tiene una y sólo una solución que está definida en el intervalo I
Ejemplo 2. Existencia y unicidad
Verificar si la función y=−1+cxx cumple la CI. y(0)=1
Solución
Reemplazando las condiciones iniciales y=1 cuando x=0 en la función, obtenemos
1=−1+c(0)0 nos queda la ecuación 1=0 lo cual es un absurdo.
Esto quiere decir que es imposible que dicha solución cumpla las condiciones iniciales
42
Ejemplo 3. Existencia y unicidad
Demuestre que la ED dxdy+y=0 tiene solución única para la CI. y(1)=3
Solución
Vamos a demostrar que la ED tiene solución única sin resolverla, para esto debemos despejar la ED.
dxdy=−y
La función f(x,y)=−y es un polinomio, por lo cual es continua en todo el plano R2
Ahora se deriva parcialmente la función respecto a y
∂y∂f=−1
Es una función constante, por lo cual es continua en todo el plano R2
Ambas funciones son continuas en cualquier rectángulo que contenga el punto (1,3), como se cumplen las condiciones 1∧2 en un intervalo con centro en x0 entonces el PVI tiene solución única
43
Ejemplo 4. Existencia y unicidad
Determine una región del plano xy para el que la ED dxdy=xy tiene solución única, cuyas gráficas pasen por un punto (x0,y0) en la región
Solución
f(x,y)=xy
Una raíz par no puede ser negativa, pero como está en el numerador puede ser 0
La función es continúa para x≥0∧y≥0 o para x≤0∧y≤0
Ahora se deriva parcialmente la función respecto a y
∂y∂f=y21x21= 21yx
Una raíz par no puede ser negativa, pero como y está en el denominador no puede ser 0
La función es continúa para x>0∧y>0 o para x < 0∧y < 0
La función es continúa en la intercepción de los dos rectángulos es decir para x>0∧y>0 o para x < 0∧y < 0
44
1.6 Repasando trigonometría y derivadas
Como acabas de ver en los contenidos vistos hasta ahora es muy importante que manejes bien los conceptos de trigonometría básica y de derivadas, por eso vamos a repasar un poco dichos contenidos, si ya los manejas bien, pasa directamente a la página 53 para practicar lo visto hasta el momento.
A continuación encontrarás una tabla con las identidades trigonométricas básicas y pitagóricas para que uses al resolver ejercicios que contengan funciones trigonométricas.
Tabla 1. Identidades trigonométricas
45
Ahora vamos a repasar derivadas. Recuerda que una derivada se puede resolver utilizando el concepto de límite, pero como ya viste la asignatura cálculo diferencial y aprendiste a derivar utilizando el concepto de límite y nuestro interés es sólo derivar de forma correcta, vamos a utilizar la siguiente tabla que contiene las derivadas básicas, la regla de la suma y de la diferencia, la regla del producto y la regla del cociente, además las derivadas de las funciones hiperbólicas las cuales vamos a utilizar cuando resolvamos ED de orden superior.
Tabla 2. Derivadas
En la siguiente escena interactiva, resuelve los ejercicios propuestos que corresponden a la derivada de la suma y la diferencia de funciones.
46
Para comenzar debes recordar las siguientes propiedades
dxd(u±v)=u′±v′
dxd(x)n=nxn−1
Debes resolver en tu cuaderno o libreta de apuntes cada ejercicio y después verificar que tu solución corresponda con la de cada ejemplo, luego da click en el botón otro ejemplo para resolver más ejercicios hasta que consideres que ya recuerdas bien el concepto de derivada.
En la siguiente escena interactiva, Que se encuentra en la página 112 del libro Cálculo diferencial interactivo, podrás repasar derivadas de polinomios
En la siguiente escena interactiva, resuelve los ejercicios propuestos que corresponden a la derivada de un producto.
47
Selecciona la opción apropiada y dale verificar. Si deseas ver el paso a paso, pulsa el botón Ver solución. Luego selecciona siguiente para ver otro ejercicio
Recuerda que es muy importante derivar los productos de forma adecuada para verificar soluciones de ecuaciones diferenciales tanto de primer orden, como las de orden superior que se verán en el próximo capítulo.
En la siguiente escena interactiva, que se encuentra en la página 120 del libro Cálculo diferencial interactivo, podrás repasar derivadas de un producto
48
En la siguiente escena interactiva, resuelve los ejercicios propuestos que corresponden a la derivada de un cociente, verifica que tu solución corresponda con la de cada ejercicio, luego selecciona siguiente para ver otro ejemplo.
Cuando se deriva un cociente es muy importante tener en cuenta el orden de la derivada. La fórmula de la derivada del producto y del cociente las encuentras en la tabla de la página 45
En la siguiente escena interactiva, que se encuentra en la página 119 del libro Cálculo diferencial interactivo, podrás repasar derivadas de un cociente
49
En la siguiente escena interactiva, resuelve cada ejercicio y selecciona la respuesta correcta, dale verificar para comprobar si es correcta, luego selecciona siguiente para ver otro ejercicio.
Podrás comprobar tu destreza para derivar un cociente en el cual hay polinomios, funciones trigonométricas, logaritmos y exponenciales.
En la siguiente escena interactiva, que se encuentra en la página 124 del libro Cálculo diferencial interactivo, podrás repasar derivadas de un cociente
50
Los siguientes ejercicios fueron diseñados por Miguel Ángel Cabezón Ochoa. En esta escena interactiva puedes realizar ejercicios aplicando regla de la cadena con diferentes funciones. Selecciona el tipo de derivada a practicar, haciendo clic en una de las opciones de la columna izquierda. La última opción muestra diferentes funciones de forma aleatoria.
Tienes a disposición una gran cantidad de ejercicios. Te sugerimos realizarlos en tu cuaderno o libreta de apuntes, cuando hayas terminado puedes verificarlos con la solución del interactivo.
Así que mucho ánimo y ponte a practicar.
En la siguiente escena interactiva, que se encuentra en la página 132 del libro Cálculo diferencial interactivo, podrás repasar derivadas utilizando regla de la cadena
51
Practiquemos
En la siguiente escena interactiva, adaptada de Plantillas con descartes_JS, encontrarás una sopa de letras con 10 palabras relacionadas a las definiciones que acabas de estudiar en las páginas anteriores sobre clasificación de las ecuaciones diferenciales, tipos de soluciones y problemas de valor inicial entre otros. Recuerda que es muy importante comprender bien los conceptos teóricos para resolver de forma adecuada los ejercicios.
52
1.6.1 Ejercicios y respuestas del capítulo 1
53
54
Capítulo 2
ECUACIONES DIFERENCIALES DE PRIMER ORDEN
56
2.1 Ecuaciones de variables separables
Definición: se dice que una E.D. de primer orden de la forma dxdy=h(y)g(x) es separable o que tiene variables separables.
Método de solución La solución de la ED es de la forma
∫h(y)dy=∫g(x)dx+c
Nota: No es necesario emplear dos constantes de integración al resolver la ED, porque si escribimos ∫h(y)dy+c1=∫g(x)dx+c2 la diferencia entre c2−c1 da como resultado una constante c
Ejemplo 1: Solución de una ED separable
Resolver la ED xsenxe−ydx−ydy=0
Solución
xsenxe−ydx−ydy=0
Primero debemos tener un diferencial a cada lado de la ecuación
57
xsenxe−ydx=ydy
Ahora separamos variables
xsenxdx=e−yydy
Pasamos Euler al numerador aplicando propiedades de potenciación
xsenxdx = yeydy
Los diferenciales indican que se debe integrar
∫xsenxdx = ∫yeydy
Ahora se integra cada lado aplicando integración por partes
∫xsenxdx
u=xdv=senxdxdu=dxv=−cosx
∫xsenxdx=−xcosx+senx+c
∫yeydy
u=ydv=eydydu=dyv=ey
∫yeydy=ye−y−e−y
Uniendo los resultados de las dos integrales obtenemos la solución de la ED
−xcosx+senx+c = ye−y−e−y
58
Factorizando obtenemos la solución implícita
e−y(y−1)=−xcosx+senx+c
Ejemplo 2: Separando variables
Resolver la ED xtanxdx−ycosxdy=0
Solución
xtanxdx−ycosxdy=0
Separamos variables
cosxxtanxdx=ydy
La integral no es por sustitución, entonces se debe aplicar identidades para integrar por sustitución.
Recuerde que tanx=cosxsenx
xcos2xsenxdx=ydy
Los diferenciales indican que se debe integrar
∫xcos2xsenxdx=∫ydy
La integral del lado izquierdo es por partes, primero integramos el dv
59
∫cos2xsenxdx
u=cosxdu=−senxdx
−∫u2du=−∫u−2du=u1=cosx1=secx
∫xcos2xsenxdx
u=xdv=cos2xsenx
du=dxv=secx
∫xcos2xsenx=xsecx−∫secxdx
Organizando la integral
xsecx−∫secx=∫ydy
Integrando a ambos lados obtenemos la solución general
Nota La integral de secante está en la tabla de la página 70
xsecx−ln∣secx+tanx∣+c=2y2
Despejando la y obtenemos la solución explícita
y=2xsecx−2ln∣secx+tanx∣+c
Recuerde que 2c=c
60
Ejemplo 3: Quitando logaritmos
Resolver la ED xdxdy=4y
Solución
Antes de resolver la ED es importante anotar que cuando sea posible se deben quitar los logaritmos
xdxdy=4y
Subimos el diferencial
xdy=4ydx
Separamos variables
ydy=4xdx
Los diferenciales indican que se debe integrar
∫ydy=4∫xdx
Integrando obtenemos
ln∣y∣=4ln∣x∣+c
Aplicamos la propiedad de logaritmos cln∣A∣=ln∣A∣c
ln∣y∣=ln∣x∣4+c
Recuerde que eln∣x∣=x
61
eln∣y∣=eln∣x∣4+c
Aplicando propiedades de potenciación
eln∣y∣=eln∣x∣4⋅ec
Teniendo en la cuenta que am.am=am+n∧ec=c
Obtenemos la solución explícita
y=cx4
Ejemplo 4: Problema de valor inicial
Resolver el PVI dxdy=−yx Sujeto a la condición y(4)=3
Solución
dxdy=−yx
Recuerde que para integrar no pueden haber diferenciales en el denominador.
ydy=−xdx
Los diferenciales indican que se debe integrar
∫ydy=−∫xdx
Integrando obtenemos
62
2y2=−2x2+c
despejando la constante
2y2+2x2=c
Simplificando
y2+x2=c
Reemplazando condiciones iniciales
32+42=c
El valor de la c
c=25
La solución particular es
Nota: Al multiplicar 2 por una constante c el resultado es c
y2+x2=25
En el ejemplo anterior solucionamos una ecuación diferencial por variables separables y en la solución general reemplazamos las condiciones iniciales para hallar el valor de la c y obtener una solución particular. noten que el valor del parámetro "c" varía de acuerdo a las condiciones iniciales
La ecuación del ejemplo 4 la vamos a graficar en la siguiente escena interactiva teniendo en la cuenta el PVI y podrás observar cómo cambia la gráfica al cambiar el valor del parámetro. La c varía desde 0.1 hasta 25
63
Las ED al igual que las ecuaciones algebraicas se pueden graficar en el plano cartesiano, algunas ED son muy sencillas de graficar, pero hay otras muy complejas.
En la siguiente escena interactiva, diseñada por Jorge Olivares Funes, podrás ver cómo cambia la gráfica del ejemplo 4 al cambiar el valor del parámetro. Al resolver la ED el parámetro se convierte en el radio de la circunferencia.
Ejemplo 5: Problema de valor inicial
Resolver la ED dxdy=−2y
Sujeta a la condición y(0)=−5
64
Solución
dxdy=−2y
Separando variables
ydy=−2dx
Los diferenciales indican que hay que integrar
∫ydy=−2∫dx
Integrando obtenemos
ln∣y∣=−2x+c
Quitando el logaritmo
y=ce−2x
Reemplazando condiciones iniciales
−5=ce0
El valor del parámetro es
c=−5
La solución particular es
y=−5e−2x
En la siguiente escena interactiva puedes ver la gráfica del ejemplo 4. Podrás ver en la escena interactiva, cómo cambia la gráfica de la ED al cambiar el valor del parámetro en un intervalo de [−15,15], el valor de la x lo puedes cambiar de [0,100] y el valor de la k lo puedes cambiar de [−5,5]
65
Al interactuar con los deslizadores podrás observar cómo cambia la gráfica, también podrás observar la ecuación diferencial y veras la solución particular de dicha ecuación
En la siguiente escena interactiva, diseñada por Jorge Olivares Funes, podrás ver cómo cambia la gráfica del ejemplo 4 al cambiar el valor de la constante.
En las siguientes dos escenas interactiva podrás resolver ED por variables separables, en la primera debes tener en cuenta que el resultado se da con el parámetro despejado.
En la primera escena interactiva debes ingresar la función f(x) (función con términos en la variable x) en la primera entrada y la función g(y) (función con términos en la variable y) en la segunda entrada, luego vas a ver la ED igualada a cero y la solución genberal
66
En la segunda escena interactiva encuentras una ED resuelta en dos pasos, los pasos los puedes ver interactuando con el deslizador, luego puedes practicar con diferentes ED cada una con su respectiva respuesta.
En la siguiente escena interactiva, diseñada por Ravinder Kumar, podrás resolver ED por variables separables y visualizar su gráfica.
En la siguiente escena interactiva, diseñada por Ravinder Kumar, podrás resolver ED por variables separables. Con el botón "Practice" podrás ver diferentes ejercicios y con el botón "Show" visualizarás el resultado.
En el vídeo podrás observar un ejemplo de ecuaciones diferenciales por variables separables.
67
En el siguiente video, realizado por Jaime H. Ramírez Rios, podrás observar un ejemplo de variables separables.
68
2.1.1 Tabla de Integrales
Como acabas de ver en los contenidos vistos hasta ahora es muy importante que manejes bien las integrales, por eso a continuación encontrarás una tabla con algunas de las integrales que encontrarás en los diferentes ecuaciones diferenciales que tendrás que resolver en este libro tanto en la segunda unidad como en las posteriores.
A continuación encontrarás una tabla con integrales para que uses al resolver ejercicios.
Tabla 3. Integrales
69
2.1.2 Repasando integrales
Como acabas de ver para resolver ED es muy importante que manejes bien las diferentes técnicas de integración, por eso vamos a repasar un poco integrales, si ya los manejas bien, pasa directamente a la página 76 para resolver ED por variables separables.
En la siguiente escena interactiva, que se encuentra en la página 40 del libro Integrando con Paco, podrás repasar integrales directas
70
En la siguiente escena interactiva vamos a resolver integrales por sustitución, recuerda que si una integral no es directa, se debe examinar si una función es la derivada de la otra, para llevarla a la forma ∫f(u)du
Puedes escoger el tipo de integral que quieres resolver, luego hazlo en tu cuaderno y verifica tu respuesta, luego da clikc en el botón otro ejercicio para realizar uno diferente.
En la siguiente escena interactiva, que se encuentra en la página 77 del libro Integrando con Paco, podrás repasar integrales por sustitución
71
La única forma de volverse experto en algún contenido, es practicando, y esto se logra resolviendo muchos ejercicios, por eso aquí tienes más integrales que se resuelven por sustitución para que sigas practicando en caso de que todavía tengas dificultad con dichas integrales.
Puedes escoger el tipo de integral que quieres resolver, luego hazlo en tu cuaderno y verifica tu respuesta, luego da clikc en el botón otro ejercicio para realizar uno diferente.
En la siguiente escena interactiva, que se encuentra en la página 70 del libro Integrando con Paco, podrás repasar integrales por partes
72
Cuando una integral no es directa, ni por sustitución, se verifica si es por partes. En la siguiente escena interactiva vamos a resolver integrales utilizando integración por partes
Puedes escoger el tipo de integral que quieres resolver, luego hazlo en tu cuaderno y verifica tu respuesta, luego da clikc en el botón otro ejercicio para realizar uno diferente.
Aquí encontrarás combinaciones de: algebraicas con exponencial, con logaritmo natural y con coseno; integrales logarítmicas y coseno con exponencial, las cuales son cíclicas.
Recuerda la fórmula de integración por partes ∫udv=uv−∫vdu
En la siguiente escena interactiva, que se encuentra en la página 106 del libro Integrando con Paco, podrás repasar integrales por partes
73
La técnica de integración por partes, fuera de que es muy utilizada, a veces presenta alguna complicación al momento de iniciar a resolverla, por eso aquí encontrarás otra escena interactiva muy amigable para que sigas repasando dicha integral
Puedes escoger el tipo de integral que quieres resolver, identificas u∧dv, luego buscas con los botones la respectiva derivada e integral, armas la solución y le das verificar para comprobar la respuesta.
Aquí encontrarás combinaciones diferentes a las vistas en la escena interactiva anterior.
Recuerda la fórmula de integración por partes ∫udv=uv−∫vdu
En la siguiente escena interactiva, que se encuentra en la página 105 del libro Integrando con Paco, podrás repasar integrales directas
74
2.1.3 Ejercicios y respuestas de la sección 2.1
75
2.2 Ecuaciones lineales
Definición: una ecuación diferencial de primer orden de la forma
a1(x)dxdy+a0(x)y=g(x)(E1)
se dice que es una ecuación lineal en la variable dependiente y
Se dice que la ecuación lineal es homogénea cuando g(x)=0; si no es no homogénea.
Al dividir ambos lados de la ecuación (E1) por el coeficiente a1(x) se obtiene la ecuación
a1(x)a1(x)dxdy+a1(x)a0(x)y=a1(x)g(x)
Donde a1(x)a0(x) lo vamos a llamar p(x)∧a1(x)g(x) lo vamos a llamar f(x)
Reemplazando estos valores en la ecuación (E1) obtenemos la ecuación que vamos a llamar Forma estándar
dxdy+p(x)y=f(x)
76
Esta ecuación tiene la propiedad de que su solución es la suma de dos soluciones
y=yc+yp
Donde yc es una solución de la ecuación homogénea ∧yp es una solución particular de la ecuación no homogénea.
Método de solución
- Convertir la ecuación lineal (E1) a la forma estándar.
- Identifique a p(x) para determinar el factor integrante e∫p(x)dx
- Multiplique la forma estándar por el factor integrante, inmediatamente la ecuación líneal se transforma en variables separables, por lo cual se debe integrar ambos lados de la ecuación.
Propiedad
Las ecuaciones lineales tienen la propiedad que el lado izquierdo de la ecuación resultante es la derivada del producto del factor integrante por la variable dependiente, esto es
dxd[e∫p(x)dxy]=e∫p(x)dxf(x)
Reemplazando el factor integrante e∫p(x)dx por F.I en la ecuación y despejando, llegamos a la ecuación.
F.Iy=∫F.If(x)dx
77
Cuando se aplica la propiedad, se resuelven de forma más sencilla las ecuaciones lineales, pues sólo se debe integrar el lado derecho de la ED.
Ejemplo 1: ED lineal homogénea
Resolver la ED dxdy=−2y
Solución Esta ED la resolvimos por variables separables, ejemplo 5, página 65. Ahora la vamos a resolver como una ecuación lineal.
dxdy=−2y
Pasando −2y al lado izquierdo de la ecuación nos queda una ED lineal homogénea.
dxdy+2y=0
Como la ecuación está en forma estándar identificamos p(x) y hallamos el factor integrante
p(x)=2F.I=e2x
Multiplicando la forma estándar por el F.I
e2xdxdy+2e2xy=0
Tenemos una ecuación de variables separables
78
∫ydy=−2∫dx
Ahora integramos ambos lados de la ecuación
ln∣y∣=−2x+c
Eliminando logaritmos llegamos a la solución
y=ce−2x
Ejemplo 2: ED lineal no homogénea
Resolver la ED dxdy=−2y+10
Solución Esta ED la resolvimos homogénea en el ejemplo 1. Ahora la vamos a resolver no homogénea, vamos a obtener la solución complementaria yc del ejemplo anterior más la solución particular yp asociada a la ED no homogénea.
dxdy=−2y+10
Pasando −2y al lado izquierdo de la ecuación, obtenemos una ecuación lineal no homogénea.
dxdy+2y=10
ecuación en forma estándar, identificamos p(x)∧ y hallamos el F.I
79
p(x)=2F.I=e2x
Multiplicando la forma estándar por el F.I
e2xdxdy+2e2xy=10e2x
Ahora tenemos una ecuación de variables separables.
e2xdxdy=−2e2xy+10e2x
Sacando factor común para separar variables.
e2xdxdy=−2e2x(y−5)
ahora pasamos el dx al lado derecho y separamos variables
y−5dy=−2dx
Los diferenciales indican que hay que integrar ambos lados de la ED
∫y−5dy=−2∫dx
Integrando
ln∣y−5∣=−2x+c
Quitando el logaritmo
eln∣y−5∣=e−2x+c
La ecuación queda de la forma
y−5=ce−2x
Despejando
y=5+ce−2x
80
La respuesta contiene la solución complementaria yc=ce−2x la cual obtuvimos en el ejemplo 1 y también la solución particular yp=5 que corresponde a la solución de la ED no homogénea, por lo cual la solución general es la suma de dos soluciones y=yc+yp
Cuando se resuelven ED lineales de primer orden, no es necesario hacer la distinción en la solución entre yc∧yp, pero cuando se resuelvan las ED de orden superior si es necesario, pues cada una se resuelve de forma diferente.
Otra forma de resolver una ED lineal es aplicando la propiedad
F.Iy=∫F.If(x)dx
Observe que al usar la propiedad sólo se debe integrar el lado derecho de la ecuación y luego despejar la variable dependiente y esto reduce las operaciones matemáticas para llegar a la solución de la ED.
Se debe llevar la ED a la forma estándar, identificar p(x) para hallar el factor integrante. En lugar de multiplicar toda la forma estándar por el factor integrante para resolver la ecuación por variables separables, se aplica la propiedad F.Iy=∫F.If(x)dx
Vamos a resolver nuevamente el ejercicio anterior aplicando la propiedad
e2xy=∫10e2xdx
Integrando
e2xy=5e2x+c
Despejando y
81
y=e2x5e2x+e2xc
La solución es
y=5+ce−2x
Nota: Todas las ecuaciones diferenciales lineales se pueden resolver aplicando la propiedad F.Iy=∫F.If(x)dx. Observe que aplicando la propiedad se reduce el método de solución, pues sólo se debe integrar el lado derecho de la ecuación.
Ejemplo 3: ED lineal no homogénea
Resolver la ED xdxdy−y=x2senx
Solución Esta ED no está en forma estándar. Se debe dividir cada término por x para llevarla a la forma estándar.
dxdy−x1y=xsenx
Como la ecuación está en forma estándar identificamos p(x)∧ hallamos el factor integrante.
p(x)=x−1F.I=x−1
Ahora aplicamos la propiedad para resolverla
82
x−1y=∫x−1xsenxdx
Aplicando potenciación para unir términos semejantes.
x−1y=∫senxdx
Integrando el lado derecho de la ecuación
x−1y=−cosx+c
Despejamos la variable dependiente y
y=x−1−cosx+x−1c
La solución es
y=−xcosx+cx
Ejemplo 4: ED lineal no homogénea
Resolver la ED dxdy+ny=enx
Solución La ED está en forma estándar.
dxdy+ny=enx
Como la ecuación está en forma estándar, identificamos p(x) y hallamos el factor integrante
83
p(x)=nF.I=enx
Aplicamos la propiedad para resolver la ecuación lineal
enxy=∫enxenxdx
Ahora aplicamos potenciación para agrupar términos semejantes.
enxy=∫e2nxdx
Integrando el lado derecho. Recuerde que esta integral es por sustitución donde u=nx
enxy=2ne2nx+c
Despejando la y
y=2nenxe2nx+enxc
La solución es
y=2nenx+ce−nx
En la siguiente escena interactiva, podrás ver la gráfica del ejemplo 4, pero también podrás interactuar con los valores de n para visualizar diferentes ecuaciones diferenciales y observar cómo cambia la gráfica a medida que se cambia el valor de n o el valor del parámetro c, también podrás observar la solución particular para cada ED
Al cambiar el valor de n en la escena interactiva, la gráfica varia en la inclinación de la curva
84
Al cambiar el valor de c en la escena interáctiva, cambia la forma de la gráfica
En la siguiente escena interactiva, diseñada por Jorge Olivares Funes, podrás ver cómo cambia la gráfica del ejemplo 4 al cambiar los valores de n y de la constante.
Es recomendable tomar tu cuaderno de apuntes y resolver varios ejercicios cuando cambies los valores de n o de c y luego compruebes tus resultados con los de la escena interactiva
En la página siguiente encontrarás una escena donde puedes ingresar el p(x) y el f(x) para resolver diferentes ED
En la siguiente escena interactiva, diseñada por Ravinder Kumar, podrás Resolver ecuaciones lineales.
85
En la siguiente escena interactiva, diseñada por Ravinder Kumar, podrás Resolver ecuaciones lineales homogéneas y no homogéneas y visualizar su gráfica.
86
En el siguiente video, realizado por Jaime H. Ramírez Rios, podrás observar un ejemplo de ecuación lineal.
En la siguiente escena interactiva, diseñada por Diego Sandoval, podrás resolver ED lineales y visualizar su campo de direcciones.
87
2.2.1 Ejercicios y respuestas de la sección 2.2
88
2.3 Ecuaciones exactas
Definición:Una ED de la forma M(x,y)dx+N(x,y)dy es una diferencial exacta en una región R del plano xy si corresponde a la diferencial de alguna función f(x,y)
Una ecuación diferencial de primer orden de la forma
M(x,y)dx+N(x,y)dy=0
es una ecuación exacta si la expresión del lado izquierdo es una diferencial exacta
Teorema. Criterio para una diferencial exacta: Sean continuas M(x,y)dx∧N(x,y)dy, con derivadas parciales continuas en una región rectangular R, definida por a<x<b,c<y<d. Entonces la condición necesaria y suficiente para que M(x,y)dx+N(x,y)dy sea una diferencial exacta es que:
∂y∂M=∂x∂N
Método de solución
- Dada una función de la forma M(x,y)dx+N(x,y)dy=0 se determina si es exacta. En caso afirmativo, existe una función para la cual
89
∂x∂f=M(x,y)
- Determinamos f si integramos M(x,y) con respecto a x, manteniendo a y constante
f(x,y)=∫M(x,y)dx+g(y)(1)
Donde g(y) es la constante de integración.
- Luego derivamos la ecuación (1) con respecto a y∧ suponemos que
∂y∂f=N(x,y)
- Por último integramos con respecto a y∧ sustituimos en la ecuación (1)
- La solución implícita de la ecuación es
f(x,y)=c
Es importante anotar que este método de solución se puede invertir comenzando con ∂y∂f=N(x,y) y el resultado es el mismo.
90
Ejemplo 1: ED exacta
Determinar si la ED 2y(y−1)dx+x(2y−1)dy=0 es exacta, si lo es resuélvala.
Solución
2y(y−1)dx+x(2y−1)dy=0
Primero debemos aplicar la propiedad distributiva
(2y2−2y)dx+(2xy−x)dy=0
Determinamos si es exacta. Derivamos parcialmente M con respecto a y dejando x constante
∂y∂M=4y−2
Derivamos parcialmente N con respecto a x dejando y constante
∂x∂N=2y−1
Como ∂y∂M=∂x∂N
Podemos concluir que la ED no es exacta, por tanto, no hay que resolverla
91
Ejemplo 2: ED exacta
Determinar si la ED 2xydx+(x2−1)dy=0 es exacta, si lo es resuélvala.
Solución
2xydx+(x2−1)dy=0
Determinamos si es exacta. Derivamos parcialmente M con respecto a y dejando x constante
∂y∂M=2x
Derivamos parcialmente N con respecto a x dejando y constante
∂x∂N=2x
Como ∂y∂M=∂x∂N, entonces, la ecuación es exacta.
2xydx+(x2−1)dy=0
Asumimos ∂x∂f=M(x,y).
∂f=(2xy)∂x
Los diferenciales indican que se debe integrar
92
∫∂f=∫2xy∂x
Como vamos a integrar con respecto a x tanto el 2 como la y son constantes
∫∂f=2y∫x∂x
Integramos
f(x,y)=x2y+g(y)(1)
Derivamos la ecuación (1) con respecto a y dejando x constante
x2+g′(y)
Ahora igualamos este resultado con N.
x2+g′(y)=x2−1
Simplificamos e integramos con respecto a y
g(y)=−y
Reemplazando este valor en la ecuación (1) y como f(x,y)=c obtenemos
c=x2y−y
En la siguiente escena interactiva, diseñada por Victor Cruz, podrás resolver ED exactas y visualizar su gráfica.
93
En la siguiente escena interactiva, diseñada por Víctor Cruz, podrás resolver ED exactas ingresando M(x,y)∧N(x,y)
94
Ejemplo 3: Resolver el PVI
Determinar si la ED es exacta, si lo es resuélvala.
(y2cosx−3x2y−2x)dx+(2ysenx−x3+ln∣y∣)dy=0
Sujeta a la condición inicial y(0)=e
Solución
Determinamos si es exacta. Derivamos parcialmente M con respecto a y dejando x constante
2ycosx−3x2
Ahora derivamos parcialmente N con respecto a x dejando y constante
2ycosx−3x2
Como ∂y∂M=∂x∂N, entonces, la ecuación es exacta.
Asumimos ∂x∂f=M(x,y)
∂x∂f=y2cosx−3x2y−2x
Despejamos e integramos a ambos lados.
95
∫∂f=∫(y2cosx−3x2y−2x)∂x
Sacamos la constante de cada integral
∫∂f=y2∫cosx∂x−3y∫x2∂x−2∫x∂x
Integramos
f(x,y)=y2senx−3y3x3−22x2+g(y)
Simplificando obtenemos
f(x,y)=y2senx−x3y−x2+g(y)(1)
Derivamos la ecuación (1) con respecto a y dejando x constante
2ysenx−x3+g′(y)
Ahora igualamos este resultado con N.
2ysenx−x3+g′(y)=2ysenx−x3+ln∣y∣
Ahora simplificamos
g′(y)=ln∣y∣
Integramos con respecto a y. Se integra por partes
u=ln∣y∣dv=dy
du=y1dyv=y
96
∫g′(y)∂y=yln∣y∣−∫∂y
El valor g(y) es
g(y)=yln∣y∣−y
Reemplazando este valor en la ecuación (1) y teniendo en cuenta que f(x,y)=c obtenemos la solución general
c=y2senx−x3y−x2+yln∣y∣−y
Ahora debemos reemplazar las condiciones iniciales
y(0)=e para obtener la solución particular
c=e2sen0−(0)3e−(0)2+eln∣e∣−e
Recuerde que ln∣e∣=1
c=e−e
El valor del parámetro es c=0
Reemplazando este valor en la solución general
0=y2senx−x3y−x2+yln∣y∣−y
Obtenemos la solución particular es
0=y2senx−x3y−x2+yln∣y∣−y
97
Ejemplo 4: ED exacta
Determinar si la ED es exacta, si lo es resuélvala.
(1+ln∣x∣+xy)dx=(1−ln∣x∣)dy
Solución
Antes de determinar si la ecuación es o no exacta se debe organizar para que cumpla la estructura M(x,y)dx+N(x,y)dy=0 que vimos en la definición.
Se puede pasar cualquiera de los dos términos al otro lado, teniendo en cuenta cual es M y cual es N Para resolver este ejercicio pasamos el término del lado derecho al lado izquierdo
(1+ln∣x∣+xy)dx−(1−ln∣x∣)dy=0
Los dos términos deben ser positivos, por lo cual sacamos factor común −1 del segundo término
(1+ln∣x∣+xy)dx+(−1+ln∣x∣)dy=0
Ahora si podemos determinar si la ED es exacta. Derivamos parcialmente M con respecto a y dejando x constante. Cuando en la ED no aparecen los dos términos al lado izquierdo o están negativos, se debe ajustar antes de comprobar si es exacta.
98
x1
Derivamos parcialmente N con respecto a x dejando y constante
x1
Como ∂y∂M=∂x∂N, entonces, la ecuación es exacta.
Asumimos ∂x∂f=M(x,y)
∫∂f=∫(1+ln∣x∣+xy)∂x
Separamos las integrales
∫∂f=∫∂x+∫ln∣x∣∂x+y∫x1∂x
Integramos M con respecto a x dejando y constante.
f(x,y)=x+xln∣x∣−x+yln∣x∣+g(y)(1)
Agrupando términos semejantes
f(x,y)=xln∣x∣+yln∣x∣+g(y)(1)
Derivamos la ecuación (1) con respecto a y dejando x constante
ln∣x∣+g′(y)
Ahora igualamos este resultado con N.
99
ln∣x∣+g′(y)=−1+ln∣x∣
Simplificamos
g′(y)=−1
Integramos con respecto a y
g(y)=−y
Reemplazando este valor en la ecuación (1)
Teniendo en cuenta que f(x,y)=c obtenemos la solución general
c=xln∣x∣+yln∣x∣−y
En la siguiente escena interactiva, diseñada por Ravinder Kumar, podrás resolver ED exactas y visualizar su gráfica.
100
En el siguiente video, realizado por Jaime H. Ramírez Rios, podrás observar un ejemplo de ecuación exacta.
En la siguiente escena interactiva, diseñada por Ravinder Kumar, podrás resolver ED exactas
101
2.3.1 Ejercicios y respuestas de la sección 2.3
102
2.4 Ecuaciones transformables a exactas
Muchas ED que no son exactas se pueden transformar en exactas usando un factor integrante. luego de transformar la ED se utiliza el método de solución visto en las ecuaciones exactas. Si la ecuación no se puede transformar en exacta se resuelve por otro método que veremos más adelante.
Para transformar una ecuación no exacta en exacta se utiliza la siguiente sustitución.
p(x)=N(x,y)∂y∂M−ax∂N
Si p(x) es función sólo de la variable x, entonces el factor integrante e∫p(x)dx la transforma en exacta
Si p(x) depende de dos variables, entonces se utiliza la siguiente sustitución
p(y)=M(x,y)∂x∂N−ay∂M
Si p(y) es función sólo de la variable y, entonces el factor integrante e∫p(y)dy la transforma en exacta
103
Si una de las dos sustituciones queda en términos de una sola variable, entonces el factor integrante la transforma en exacta. es importante anotar que no todas las ED de primer orden tienen factor integrante y la sustitución tampoco permite encontrar todos los factores integrantes posibles de la ED. Sólo permite encontrar factores integrantes univariables
Ejemplo 1: ED reducible a exacta
Determinar el factor integrante para resolver la ED
xydx+(2x2+3y2−20)dy=0
Solución
Determinamos si es exacta. Derivamos parcialmente M con respecto a y dejando x constante
x
Ahora derivamos parcialmente N con respecto a x dejando y constante
4x
Como ∂y∂M=∂x∂N, entonces, la ecuación no es exacta.
104
Realizamos la sustitución p(x)=N(x,y)∂y∂M−ax∂N
2x2+3y2−20x−4x=2x2+3y2−20−3x
Como el resultado depende de dos variables buscamos p(y)
xy4x−x=xy3x=y3
Como el resultado depende de una variable, utilizamos p(y) para hallar el factor integrante
F.I=e∫y3dy=y3
Ahora multiplicamos la ED no exacta por el F.I
[xydx+(2x2+3y2−20)dy=0]y3
xy4dx+(2x2y3+3y5−20y3)dy=0
Derivamos M con respecto a y dejando x constante
4xy3
Derivamos N con respecto a x dejando y constante
4xy3
Como ∂y∂M=∂x∂N, entonces, la ecuación es exacta.
105
Asumimos ∂x∂f=M(x,y)
∂f=(xy4)∂x
Los diferenciales indican que se debe integrar
∫∂f=∫(xy4)∂x
Integramos M con respecto a x dejando y constante
f(x,y)=2x2y4+g(y)(1)
Derivamos la ecuación (1) con respecto a y dejando x constante
2x2y3+g′(y)
Igualamos este resultado con N.
2x2y3+g′(y)=2x2y3+3y5−20y3
Simplificando obtenemos
g′(y)=3y5−20y3
Integramos con respecto a y
∫g′(y)∂y=∫(3y5−20y3)∂y
Integral directa
g(y)=63y6−420y4
Simplificando
106
g(y)=2y6−5y4
Reemplazando este valor en la ecuación (1) y teniendo en cuenta que f(x,y)=c obtenemos la solución general
c=2x2y4+2y6−5y4
Ejemplo 2: ED reducible a exacta
Determinar el factor integrante para resolver la ED
(cos2y−senx)dx+(−2tanxsen2y)dy=0
Solución
Determinamos si es exacta. Derivamos parcialmente M con respecto a y dejando x constante
−2sen2y
Ahora derivamos parcialmente N con respecto a x dejando y constante
−2sec2xsen2y
107
Como ∂y∂M=∂x∂N, entonces, la ecuación no es exacta.
Como el resultado depende de una variable, utilizamos p(x) para hallar el factor integrante
F.I=e∫−tanxdx=cosx
Ahora multiplicamos la ED no exacta por el F.I
[(cos2y−senx)dx+(−2tanxsen2y)dy=0]cosx
(cosxcos2y−senxcosx)dx+(−2senxsen2y)=0
Tenga en cuenta que se aplicó la identidad cosxtanx=senx
Derivamos M con respecto a y dejando x constante
108
−2cosxsen2y
Derivamos N con respecto a x dejando y constante
−2cosxsen2y
.
Como ∂y∂M=∂x∂N, entonces, la ecuación es exacta.
Asumimos ∂x∂f=M(x,y)
∂f=(cosxcos2y−senxcosx)∂x
Los diferenciales indican que se debe integrar
∫∂f=∫(cosxcos2y−senxcosx)∂x
Integramos M con respecto a x dejando y constante
f(x,y)=senxcos2y+2cos2x+g(y)(1)
Derivamos la ecuación (1) con respecto a y dejando x constante
−2senxsen2y+g′(y)
Ahora igualamos este resultado con N.
−2senxsen2y+g′(y)=−2senxsen2y
109
Al simplificar vemos que g′(y)=0. Cuando integramos el 0 el resultado es 0+c como la constante sólo se tiene en cuenta en la respuesta podemos decir que:
g(y)=0
Reemplazando este valor en la ecuación (1)∧ teniendo en cuenta que f(x,y)=c obtenemos la solución general
c=senxcos2y+2cos2x
Ejemplo 3: ED reducible a exacta
Determinar el factor integrante para resolver la ED
(2xyln∣y∣)dx+(x2+y2y2+1)dy=0
Solución
Determinamos si es exacta. Derivamos parcialmente M con respecto a y dejando x constante
2x+2xln∣y∣
Derivamos parcialmente N con respecto a x dejando y constante
110
2x
Como ∂y∂M=∂x∂N, entonces, la ecuación no es exacta.
Realizamos la sustitución p(x)=N(x,y)∂y∂M−ax∂N
x2+y2y2+12x+2xln∣y∣−2x
Simplificando
x2+y2y2+12xln∣y∣
Como el resultado depende de dos variables, utilizamos la sustitución p(y)
2xyln∣y∣2x−2x−2xln∣y∣
Simplificando
2xyln∣y∣−2xln∣y∣=y−1
Hallamos el factor integrante
F.I=e∫y−1dy=y−1
Ahora multiplicamos la ED no exacta por el F.I
(2xln∣y∣)dx+(x2y−1+yy2+1)dy=0
Derivamos M con respecto a y dejando x constante
111
2xy−1
Derivamos N con respecto a x dejando y constante
2xy−1
Como ∂y∂M=∂x∂N, entonces, la ecuación es exacta.
Asumimos ∂x∂f=M(x,y)
∂f=(2xln∣y∣)∂x
Los diferenciales indican que se debe integrar
∫∂f=∫2xln∣y∣∂x
Integramos M con respecto a x dejando y constante
f(x,y)=x2ln∣y∣+g(y)(1)
Derivamos la ecuación (1) con respecto a y dejando x constante
x2y−1+g′(y)
Ahora igualamos este resultado con N.
x2y−1+g′(y)=x2y−1+yy2+1
Simplificando
112
g′(y)=yy2+1
Ahora integramos esta ecuación con respecto a y.
∫g′(y)∂y=∫yy2+1∂y
Integral por sustitución
u=y2+1∧du=2ydy
Reemplazando en la integral
g′(y)=21∫u21du
Integrando obtenemos
g(y)=31∫u23
Reemplazamos el valor de u en la respuesta
g(y)=31(y2+1)23
Entonces g(y) es:
g(y)=31(y2+1)3
Reemplazando este valor en la ecuación (1)∧ teniendo en cuenta que f(x,y)=c obtenemos la solución general
c=x2ln∣y∣+31(y2+1)3
113
Ejemplo 4: Resolver el PVI
Determinar el factor integrante para resolver la ED
xtanxdx−ycosxdy=0
Sujeta a la condición inicial y(0)=2
Solución
Debemos organizar la ED para que cumpla la estructura M(x,y)dx+N(x,y)dy=0
xtanxdx+(−ycosx)dy=0
Determinamos si es exacta. Derivamos parcialmente M con respecto a y dejando x constante
0
Ahora derivamos parcialmente N con respecto a x dejando y constante
ysenx
Como ∂y∂M=∂x∂N, entonces, la ecuación no es exacta.
Buscamos un factor integrante para convertirla en exacta
114
Realizamos la sustitución p(x)=N(x,y)∂y∂M−ax∂N
−ycosx0−ysenx
Simplificando
cosxsenx
hallamos el factor integrante
F.I=e∫tanxdx=secx
La integral de secante está en la página 70
Ahora multiplicamos la ED no exacta por el F.I
[xtanxdx+(−ycosx)dy=0]secx
(xsecxtanx)dx+(−y)dy=0
Tenga en cuenta que se aplicó la identidad cosxsecx=1
Derivamos M con respecto a y dejando x constante
0
Derivamos N con respecto a x dejando y constante
0
Como ∂y∂M=∂x∂N, entonces, la ecuación es exacta.
115
Asumimos ∂x∂f=M(x,y)
∫∂f=∫xsecxtanx∂x
Integramos M con respecto a x dejando y constante.
u=xdv=secxtanxdxdu=dxv=secx
Aplicando la fórmula de integración por partes
∫∂f=xsecx−∫secxdx
Integrando la secx obtenemos
f(x,y)=xsecx−ln∣secx+tanx∣+g(y)(1)
Derivamos la ecuación (1) con respecto a y dejando x constante y la igualamos con N
g′(y)=−y
Integramos con respecto a y.
∫g′(y)∂y=∫−y∂y
Integral directa
g(y)=−2y2
Reemplazando en la ecuación (1)∧ teniendo en cuenta que f(x,y)=c
c=xsecx−ln∣secx+tanx∣−2y2
116
Ahora reemplazamos las condiciones iniciales y(0)=2
c=0sec0−ln∣sec0+tan0∣−222c=−2
Entonces la solución particular de la ED es
−2=xsecx−ln∣secx+tanx∣−2y2
En el siguiente video, realizado por Jaime H. Ramírez Rios, podrás observar un ejemplo de una ecuación reducible a exacta.
117
2.4.1 Ejercicios y respuestas de la sección 2.4
118
2.5 Ecuaciones homogéneas
Si una ED tiene la forma M(x,y)dx+N(x,y)dy=0
Tiene la siguiente propiedad
a) M(tx,ty)=tnM(x,y)
b) N(tx,ty)=tnN(x,y)
Entonces se dice que tiene coeficientes homogéneos o que es una ecuación diferencial homogénea
Nota: Una ecuación diferencial homogénea siempre se puede llevar a una ecuación diferencial de variables separables por medio de una sustitución algebraica adecuada
Definición: Se dice que f(x,y) es una función homogénea de grado n si para algún número R
f(tx,ty)=tnf(x,y)
Método de solución
- Si f(x,y) es una función homogénea de grado n, es posible hacer una sustitución con el fin de llevarla a una E.D. de variables separables.
- La sustitución que normalmente se hace es:
119
y=ux
Recuerde que, con esta sustitución, la ecuación se vuelve variables separables, sí la integral es muy difícil o imposible de resolverse utiliza la siguiente sustitución.
x=vy
Ejemplo 1: ED homogénea
Resolver la ED
xdxdy=(y+xexy)
Solución
Primero se determina si la ecuación es homogénea
txdxdy=(ty+txetxty)
Se simplifica ∧ se saca factor común
txdxdy=t(y+xexy)
120
Se determina que la ecuación es homogénea de grado 1. Ahora se hace el cambio de variable
y=ux derivando dy=udx+xdu
xdxudx+xdu=(ux+xeu)
Al hacer la sustitución la ecuación se vuelve variables separables, subimos dx y sacamos factor común x.
x(udx+xdu)=x(u+eu)dx
Simplificamos la x∧ hacemos distributiva
udx+xdu=udx+eudx
Agrupamos los dx a un lado ∧ el du al otro lado
xdu=udx−udx+eudx
Simplificamos ∧ separamos variables
eudu=xdx
Los diferenciales indican que se debe integrar
∫e−udu=∫xdx
integrando obtenemos
−e−u=ln∣x∣+c
Sustituimos nuevamente la u
121
−e−xy=ln∣x∣+c
Multiplicamos por −1 para dar la solución
e−xy=−ln∣x∣+c
Ejemplo 2: ED homogénea
Resolver la ED
(y2+yx)dx−x2dy=0
Solución
Primero se determina si la ecuación es homogénea
(t2y2+tytx)dx−t2x2dy=0
Se multiplica el parámetro
(t2y2+t2yx)dx−t2x2dy=0
Se saca factor común
t2(y2+yx)dx−t2x2dy=0
Se determina que la ecuación es homogénea de grado 2. Ahora se hace el cambio de variable
122
y=ux derivando dy=udx+xdu
((ux)2+uxx)dx−x2(udx+xdu)=0
Para resolver la ecuación se saca factor común x2∧ pasamos el segundo término a la derecha.
x2(u2+u)dx=x2(udx+xdu)
Simplificamos x2∧ hacemos distributiva
u2dx+udx=udx+xdu
Se deben separar los diferenciales a diferente lado para poder separar variables
u2dx=xdu
Ahora se separan variables
xdx=u2du
Los diferenciales indican que se debe integrar
∫xdx=∫u−2du
integrando obtenemos
ln∣x∣+c=−u−1
Sustituimos nuevamente la u
ln∣x∣+c=−yx
despejando la y, queda la solución
123
y=ln∣x∣+c−x
Ejemplo 3: Resolver el PVI
Resolver la ED
(y+xcotxy)dx−xdy=0
Solución
Primero se determina si la ecuación es homogénea
(ty+txcottxty)dx−txdy=0
Se simplifica ∧ se saca factor común
t(y+xcotxy)dx−txdy=0
Como los dos términos tienen el mismo factor t, se determina que la ecuación es homogénea de grado 1. Ahora se hace el cambio de variable
y=ux derivando dy=udx+xdu
(ux+xcotu)dx−x(udx+xdu)=0
124
Al hacer la sustitución la ecuación se vuelve variables separables, ahora sacamos factor común x∧ pasamos el segundo término a la derecha.
x(u+cotu)dx=x(udx+xdu)
Simplificamos x∧ hacemos distributiva
udx+cotudx=udx+xdu
separamos diferenciales
udx+cotudx−udx=xdu
Simplificamos ∧ separamos variables
xdx=cotudu
Los diferenciales indican que se debe integrar ∧ se aplica la identidad cotu1=tanu
∫xdx=∫tanudu
integrando obtenemos
ln∣x∣+c=ln∣secu∣
Aplicamos axioma de igualdad
eln∣x∣+c=eln∣secu∣
Quitando logaritmos. Recuerde que eln∣x∣=x
cx=secu
Sustituyendo la u
125
cx=secxy
Axioma de igualdad para despejar y
sec−1(cx)=sec−1secxy
Operando
sec−1(cx)=xy
despejando la y queda la solución
y=xsec−1(cx)
Ejemplo 4: Resolver la ED homogénea
Resolver la ED
2x2ydx=(3x3+y3)dy
Solución
Primero se determina si la ecuación es homogénea
2t2x2tydx=(3t3x3+t3y3)dy
Operando y sacando factor común
2t3x2ydx=t3(3x3+y3)dy
126
Como los dos términos tienen el mismo factor t3, se determina que la ecuación es homogénea de grado 3. Ahora hacemos el cambio de variable
y=ux derivando dy=udx+xdu
2x2uxdx=(3x3+(ux)3)(udx+xdu)
Al hacer la sustitución la ecuación se vuelve variables separables, ahora sacamos factor común x3
2x3udx=x3(3+u3)(udx+xdu)
Simplificamos x3∧ hacemos distributiva
2udx=3udx+3xdu+u4dx+xu3du
separamos diferenciales
−udx−u4dx=3xdu+xu3du
Sacamos factor común
−dx(u+u4)=xdu(3+u3)
Separando variables
−xdx=u+u43+u3du
Los diferenciales indican que se debe integrar
−∫xdx=∫u+u43+u3du
127
La integral del lado derecho es muy difícil de resolver, pues no se puede hacer por sustitución ∧ por fracciones parciales queda un término irreductible, cuando esto ocurre, se debe volver a comenzar el ejercicio pero utilizando la otra sustitución x=vy
x=vy derivando dx=vdy+ydv
2(vy)2y(vdy+ydv)=(3(vy)3+y3)dy
Operando tenemos
2v2y3(vdy+ydv)=(3v3y3+y3)dy
Sacamos factor común y3
2v2y3(vdy+ydv)=y3(3v3+1)dy
Simplificando ∧ haciendo distributiva
2v3dy+2v2ydv=3v3dy+dy
separamos diferenciales y agrupamos términos semejantes
2v2ydv=v3dy+dy
Sacamos factor común
2v2ydv=dy(v3+1)
Separando variables
v3+12v2dv=ydy
Los diferenciales indican que se debe integrar
128
∫v3+12v2dv=∫ydy
La integral se resuelve por sustitución
Es necesario hacer un cambio de variable para resolver la integral del lado izquierdo cambiamos v3+1 por u
u=v3+1⇒du=3v2dv⇒3du=v2dv
32∫vdv=∫ydy
Integrando obtenemos
32ln∣u∣+c=ln∣y∣
Se organiza el término
2ln∣∣v3+1∣∣+c=3ln∣y∣
Aplicamos propiedades de logaritmos
ln∣∣v3+1∣∣2+c=ln∣y∣3
Aplicamos axioma de igualdad
eln∣v3+1∣2+c=eln∣y∣3
Quitando logaritmos
c(v3+1)2=y3
Sustituyendo la v
c(y3x3+1)2=y3
Operando fraccionarios
129
c(y3x3+y3)2=y3
Aplicando propiedades de potenciación
(y3)2C(x3+y3)2=y3
Pasando y6 a multiplicar obtenemos la solución
y9=c(x3+y3)2
En el siguiente video, realizado por Jaime H. Ramírez Rios, podrás observar un ejemplo de ecuación homogénea.
En la siguiente escena interactiva, diseñada por Ravinder Kumar, podrás resolver ecuaciones homogéneas.
130
En la siguiente escena interactiva, diseñada por Ravinder Kumar, podrás resolver ecuaciones homogéneas.
131
2.5.1 Ejercicios y respuestas de la sección 2.5
132
2.6 ecuación de Bernoulli
Algunas veces es necesario hacer un cambio de variable para transformar una ecuación diferencial no lineal en lineal. Lo primero que se debe hacer es llevar la ecuación a la forma estándar, lego hacer una sustitución. Si la ecuación presenta la forma:
dxdy+p(x)y=f(x)yn (1)
Donde n es cualquier número real, es la ecuación de Bernoulli.
Nota: Si n=0n=1 la ecuación es lineal.
Si n=1 la sustitución
u=y1−n
Reduce cualquier ecuación de la forma (1) en una ecuación lineal.
La ecuación diferencial de Bernoulli es una ecuación diferencial ordinaria de primer orden, formulada por Jacob Bernoulli. Esta ecuación fue transformada, por Gottfried Leibniz en 1693 y por Johann Bernoulli en 1697, es una ecuación diferencial lineal de primer orden, mediante la sustitución u=y1−n
133
Ejemplo 1: Resolver la ED de Bernoulli
Resolver la ED
xdxdy+y=x2y2
Solución
Se determina el valor de n para realizar la sustitución. En este ejercicio n=2
u=y1−2⇒u=y−1
Despejando la y
y=u−1
Derivando la y
dy=−u−2du
Reemplazando la y∧ el dy en la ecuación de Bernoulli
xdx−u−2du+u−1=x2(u−1)2
Después de hacer el cambio de variable, la ecuación se transforma en lineal, ahora se debe llevar a la forma estándar
−u−2xdx−u−2xdu+−u−2xu−1=−u−2xx2u−2
Simplificando
134
dxdu−x1u=−x
La ecuación es lineal en la variable dependiente u∧ además está en forma estándar. Ahora identificamos p(x) para hallar el factor integrante.
p(x)=−x1F.I=x−1
Ahora aplicamos la propiedad para resolver la ED lineal, multiplicando la forma estándar por el F.I
x−1u=−∫x−1xdx
Ahora debemos integrar el lado derecho de la ecuación
x−1u=−x+c
Despejamos la variable dependiente u
u=−x−1x+x−1c
Sustituyendo u∧ aplicando potenciación
y1=−x2+cx
Despejando y se llega a la solución
y=−x2+cx1
135
Ejemplo 2: Resolver la ED de Bernoulli
Resolver la ED
dxdy=y(xy3−1)
Solución
Se aplica propiedad distributiva
dxdy=xy4−y
Se organiza la ecuación para que la variable lineal quede al lado izquierdo
dxdy+y=xy4
La ecuación ya está en la forma de la ecuación de Bernoulli
dxdy+p(x)y=f(x)yn, ahora se determina el valor de n para realizar la sustitución. En este ejercicio n=4
u=y1−4⇒u=y−3
Despejando la y
y=u−31
Derivando la y
136
dy=−31u−34du
Reemplazando la y∧ el dy en la ecuación de Bernoulli
dx−31u−34du+u−31=xu−34
Después de hacer el cambio de variable, la ecuación se transforma en lineal, ahora se debe llevar a la forma estándar
integrando el lado derecho de forma directa, obtenemos
x6u=5−9x5+c
Despejando la u y simplificando
141
u=5x−9+x6c
Sustituyendo u
y−3=5x−9+x6c
Reemplazando valores iniciales
(21)−3=5(1)−9+(1)6c
Operando
8=5−9+c
Despejando obtenemos el valor de c
c=549
Reemplazando en la respuesta obtenemos
y−3=5x−9+5x649
Operando el lado derecho
y−3=5x6−9x5+49
Bajamos la y para volver el exponente positivo
y31=5x6−9x5+49
Despejando la y∧ sacando raíz cúbica obtenemos la siguiente respuesta
y=3−9x5+495x6
142
En la siguiente escena interactiva, diseñada por Raul Enrique Escobar Caro, podrás resolver la ecuación de Bernoulli.
En el siguiente video, realizado por Jaime H. Ramírez Rios, podrás observar un ejemplo de ecuación de Bernoulli.
143
2.6.1 Ejercicios y respuestas de la sección 2.6
144
2.7 Ejercicios de modelado
Las ED tienen una importancia fundamental en las matemáticas debido a que es aplicable a casi todas las ramas del saber, entre ellas están: la geometría, la física, la química y la biología. Éstas permiten identificar y aplicar técnicas y conceptos matemáticos que admiten obtener e interpretar la solución de una situación problema en la que intervienen relaciones entre tasas o razones de cambio de variables definidas en diversos contextos.
Algunas de las aplicaciones de las ED de primer orden son: el crecimiento y decrecimiento poblacional, la desintegración radiactiva, el fechado de fósiles con carbono 14, el cambio en las temperaturas y los circuitos en serie, entre otros.
Las ED se utilizan para describir los comportamientos de algunos fenómenos de la vida real en términos matemáticos, dicho proceso recibe el nombre de modelo matemático y se utiliza para comprender las variables que lo conforman, la relación entre ellas, además de predecir el comportamiento de dicho fenómeno.
Cuando dicho modelo lo debe resolver un científico es necesario conocer la fuerza de gravedad que ejerce la tierra sobre dicho cuerpo o la temperatura de ebullición del agua dependiendo de la altitud donde se realice el experimento, más para efectos de este libro, el objetivo se reduce a resolver modelos de baja resolución por lo cual la fuerza de gravedad que se utiliza es de 9.8seg2m en el sistema MKS o 32seg2pies en el sistema inglés y la temperatura de ebullición del agua que se utilizará es de 100°C.
Es por esto que para terminar el capítulo se introduce el modelado de problemas a partir de las leyes físicas que los rigen y que conducen a EDO.
145
2.7.1 Crecimiento y decrecimiento
En 1798 el economista inglés Thomas Malthus intentó modelar el crecimiento de la población humana por medio de las matemáticas, su teoría dice que la rapidez con que varía la población de un país es proporcional al número de personas presentes en ese momento. En términos matemáticos, si p(t) denota la población en función del tiempo, su hipótesis se puede expresar cómo
dtdpαp
Recuerde que en matemáticas, si dos cantidades x∧y son proporcionales, se escribe xαy significando que x es un múltiplo constante de y, es decir, x=ky, por tanto, la ecuación dtdpαp se puede reescribir como
dtdp=kp
Donde k es una constante de proporcionalidad.
Podemos observar que la ecuación contiene la derivada de una variable dependiente (población) con respecto a una variable independiente (tiempo), a esta ecuación se le conoce como ED y al resolverla teniendo en cuenta las condiciones iniciales del modelo matemático podemos predecir la cantidad de personas presentes en cualquier momento t. Aunque es un modelo simple bastante exacto, puede fallar si no se tienen en cuenta algunos factores como, por ejemplo, migración, inmigración y catástrofes naturales entre otros. Cuando se habla de seres vivos se puede usar para modelar el crecimiento de pequeñas poblaciones en intervalos cortos de tiempo.
146
Observemos que esta es una ED de primer orden, al organizarla nos queda
dtdp−kp=0
Esta es una ED lineal y está en forma estándar. Identificamos p(t) y hallamos el factor integrante
p(t)=−kF.I=e−kt
Aplicando la propiedad
e−ktp=∫0dx
Integrando
e−ktp=c
despejando la y
p=e−ktc
La solución es
p=cekt
Esta es la solución de la ED para resolver ejercicios de crecimiento y decrecimiento poblacional. Esta ecuación es la misma que se utiliza para ejercicios de descomposición radiactiva, carbono 14 (C-14) e interés compuesto entre otros, por lo cual para efectos de simplificación de variables cuando no estemos trabajando ejercicios de crecimiento y decrecimiento poblacional, la ecuación la vamos a trabajar con la variable x
147
x=cekt
Ejemplo 1: Crecimiento poblacional
La población de Colombia crece a una tasa proporcional a la población presente en el tiempo (t). En el año 2005 el censo arrojó que había 41′468.384 habitantes, en 2018 el número de habitantes era de 48′258.494. ¿Cuál será el número de personas en 2030? Y ¿Cuándo habrá 100′000.000 de habitantes en Colombia?
Solución
Condiciones iniciales del problema. En 2005, el cual va a ser el tiempo inicial t=0 el número de habitantes era de 41′468.384 y en 2018, o sea, 13 años más tarde el número de habitantes era de 48′258.494, las otras dos condiciones son las preguntas.
Condiciones iniciales.
t=0p=41′468.384
t=13p=48′258.494
t=12p=?
t=?p=100′000.000
148
Vamos a comenzar reemplazando el primer dato p=41′468.384 cuando t=0 en la ecuación p=cekt(1).
Note que la ecuación (1) tiene cuatro incógnitas
41′468.384=cek(0)⇒c=41′468.384
Reemplazando este valor en la ecuación (1) queda de la siguiente forma p=41′468.384ekt(2).
Note que la ecuación (2) ya tiene sólo tres incógnitas
Ahora reemplazamos p=48′258.494 cuando t=13 en la ecuación (2)
48′258.494=41′468.384ek(13) despejando queda
41′468.38448′258.494=e13k
Aplicamos logaritmo natural a ambos lados de la ecuación
ln∣∣41′468.38448′258.494∣∣=ln∣∣e13k∣∣
0.1516405482=13k despejando obtenemos el valor de k
k=0.01166465755
Reemplazando este valor en la ecuación (2) obtenemos
p=41′468.384e0.01166465755t(3).
149
Ya podemos responder la primera pregunta. Cuál será la población en 2030, es decir, en t=25
p=41′468.384e0.01166465755(25)⇒p=55′509.143
La población en Colombia en 2030 será de 55′509.143 habitantes
Ahora vamos a responder la segunda pregunta que dice, ¿Cuándo habrá 100′000.000 de habitantes en Colombia?
100′000.000=41′468.384e0.01166465755(t)⇒t=75.46
A mediados del año 2093 se espera que la población en Colombia sea de 100′000.000 de habitantes.
Ejemplo 2: Datado con carbono
En una cueva de Sudáfrica se halló un cráneo humano junto con los restos de una hoguera. Los arqueólogos creen que la edad del cráneo sea igual a la edad de la hoguera. Se determinó que sólo queda el 2% de la cantidad original de C−14 en los restos de la madera. Estime la edad del cráneo.
150
En 1955, Willard Frank Libby publicó una obra titulada Radiocarbon datingel en la cual habla sobre la teoría del fechado con radiocarbono (carbono 14) que se produce en la atmósfera por acción de la radiación cósmica sobre el nitrógeno.
El carbono 14 es un isótopo radioactivo inestable que se encuentra en la atmósfera, al oxidarse se transforma en dióxido de carbono que absorben las plantas y que se transmiten a otros organismos cuando éstas son consumidas. Este intercambio constante de carbono entre la atmósfera y los seres vivos se ve interrumpida cuando éstos mueren y es cuando comienza a desintegrarse en fracciones constantes, en los laboratorios especializados miden la cantidad de carbono 14 que se conservó de la muestra al paso del tiempo.
Solución
Para resolver el ejercicio es necesario plantear las condiciones iniciales del problema, teniendo en cuenta que la vida media del C−14 es 5.600 años. Se sabe que inicialmente el cráneo tenía el 100% de c−14. Siempre que un ejercicios nos hable de porcentaje lo vamos a dividir por cien.
Condiciones iniciales.
t=0p=1
t=5600p=0.5
t=?p=0.02
151
Vamos a comenzar reemplazando el primer dato p=1 cuando t=0 en la ecuación p=cekt(1).
1=cek(0)⇒c=1
Reemplazando este valor en la ecuación (1) queda de la siguiente forma p=ekt(2).
Ahora reemplazamos p=0.5 cuando t=5600 en la ecuación (2)
0.5=ek(5600)
Aplicamos logaritmo natural a ambos lados de la ecuación
ln∣0.5∣=ln∣∣e5600k∣∣
−0.6931471806=5600k despejando obtenemos el valor de k
k=−0.0001237762
Reemplazando este valor en la ecuación (2) obtenemos
p=e−0.0001237762t(3).
Ahora vamos a calcular la edad del cráneo, teniendo en cuenta que al momento de encontrar los restos de la madera sólo quedaba el 2% de la cantidad original de C−14.
0.02=e−0.0001237762t⇒t=31.607,1.
La edad aproximada del cráneo es de 31.607,6 años.
152
Ejemplo 3: Crecimiento poblacional
La población de cierta ciudad aumenta proporcionalmente al número de habitantes que hay en un momento dado en ella. Si después de 5 años, la población se ha triplicado y después de 8 años es de 45.000 habitantes, hallar el número de ciudadanos que había inicialmente.
Solución
Se deben plantear las condiciones iniciales del problema. Como la incógnita es la población inicial, la vamos a llamar p0 para poder trabajar la ecuación. Note que cuando el enunciado dice que a los 5 años la población se ha triplicado, se representa 3po
Condiciones iniciales.
t=0p=p0
t=5p=3p0
t=8p=45.000
Reemplazamos el primer dato p=p0 cuando t=0 en la ecuación p=cekt(1).
p0=cek(0)⇒c=p0
Reemplazando este valor en la ecuación (1) queda de la siguiente forma p=p0ekt(2).
153
Ahora reemplazamos p=3p0 cuando t=5 en la ecuación (2)
3p0=p0ek(5) despejando queda
p03p0=e5k
Aplicamos logaritmo natural a ambos lados de la ecuación
ln∣3∣=ln∣∣e5k∣∣
1.098612289=5k despejando obtenemos el valor de k
k=0.2197224577
Reemplazando este valor en la ecuación (2) obtenemos
p=p0e0.2197224577t(3).
Ahora vamos a calcular la cantidad inicial de habitantes
45000=p0e0.2197224577(8)
Despejando p0
p0=e0.2197224577(8)45000.
La cantidad inicial de habitantes en la ciudad es 7.759 habitantes
154
Ejemplo 4: Descomposición radiactiva
Se ha detectado que el 0.5% de una sustancia radioactiva desaparece en 12 años. ¿Qué porcentaje desaparecerá en 1.000 años? ¿Cuál es la vida media de dicha sustancia?
Solución
Condiciones iniciales.
t=0
t=12
t=?
t=1000
x=1
x=0.995
x=0.5
x=?
Reemplazamos el primer dato x=1 cuando t=0 en la ecuación p=cekt(1).
1=cek(0)⇒c=1
Reemplazando este valor en la ecuación (1) queda de la siguiente forma x=ekt(2).
Ahora reemplazamos x=0.995 cuando t=12 en la ecuación (2)
155
0.995=ek(12) aplicamos logaritmo natural a ambos lados de la ecuación
ln∣0.995∣=ln∣∣e12k∣∣
Despejando obtenemos k=−0.00041771181
Reemplazando este valor en la ecuación (2) obtenemos
x=e−0.00041771181t(3).
Ahora vamos a calcular qué porcentaje desaparecerá en 1.000 años
x=e−0.00041771181(1000)
x=65.85 Esto quiere decir que a los mil años quedará el 65.85%
En 1.000 años desaparecerá el 34.15%
Ahora vamos a calcular la vida media de la sustancia radiactiva
0.5=e−0.00041771181(t).
ln∣0.5∣=ln∣∣e−0.00041771181(t)∣∣
t=−0.00041771181−0.6931471806
La vida media de la sustancia radiactiva es 1659.4 años
156
Ejemplo 5: Interés compuesto
Previo al nacimiento de su primer hijo, una pareja deposita 5.000 dólares en una cuenta que paga el 8% de interés compuesto. Los pagos de interés son acumulables al capital. ¿Cuánto habrá en la cuenta en el dieciochoavo cumpleaños del niño?
Solución
En los ejercicios de modelado la k es la constante de crecimiento o decrecimiento del ejercicio, es la que determina a qué velocidad está creciendo o decreciendo una población o se descompone el carbono 14 o se desintegra determinado material radiactivo, en los ejercicios de interés compuesto el interés determina la rapidez de crecimiento del capital, es decir, que la k es el interés compuesto, para este ejercicio k=0.08
Condiciones iniciales.
t=0
t=18
x=5000
x=?
Reemplazamos el primer dato x=5000 cuando t=0 en la ecuación p=cekt(1).
5000=cek(0)⇒c=5000
157
Reemplazando este valor en la ecuación (1) queda de la siguiente forma x=5000ekt(2).
Ahora reemplazamos k=0.08∧t=18 en la ecuación (2)
x=5000e(0.08)(18)x=21103,47.
Al dieciochoavo cumpleaños del niño habrá en la cuenta 21103,47 dólares
En la página 199 puedes observar un vídeo de aplicación de crecimiento y decrecimiento.
En la siguiente escena interactiva, podrás resolver ejercicios de crecimiento poblacional, debes ingresar las condiciones iniciales del ejercicio.
158
2.7.2 Ejercicios y respuestas de la sección 2.7.1
159
2.7.3 Temperaturas
La ley de enfriamiento de Newton dice que en un cuerpo que se está enfriando, la rapidez con que la temperatura T cambia, es proporcional a la diferencia entre la temperatura del cuerpo y la temperatura constante Tm del medio que lo rodea, esto es:
dtdT=k(T−Tm)
Podemos ver que esta es una ED de variables separables. Separando variables tenemos
(T−Tm)dT=kdt
Integrando ambos lados de la ED. El lado izquierdo se integra por sustitución.
ln∣T−Tm∣=kt+c
Quitando el logaritmo
eln∣T−Tm∣=ekt+c
T−Tm=cekt
Despejando la temperatura obtenemos la solución
T=Tm+cekt
160
La solución es muy similar a la de crecimiento y decrecimiento, se diferencia en que tenemos un dato más que es la temperatura del medio Tm, por esto los ejercicios se resuelven de forma similar.
Ejemplo 1: Temperatura
Una taza de café se sirve a una temperatura de 70∘C en un cuarto cuya temperatura es de 17∘C. Dos minutos más tarde la temperatura del café es de 60∘C. ¿Después de cuánto tiempo la taza de café alcanza la temperatura de 20∘C?
Solución
Se plantean las condiciones iniciales del problema.
Condiciones iniciales.
t=0
t=2
t=?
T=70∘C
T=60∘C
T=20∘C
Tm=17∘C
Vamos a comenzar reemplazando el primer dato T=70∘C cuando t=0 en la ecuación T=Tm+cekt(1).
161
70=17+cek(0)⇒c=53
Reemplazando este valor en la ecuación (1) queda de la siguiente forma T=17+53ekt(2).
Ahora reemplazamos T=60 cuando t=2 en la ecuación (2)
60=17+53ek(2) despejando queda 5343=ek(2)
Aplicamos logaritmo natural a ambos lados de la ecuación
ln∣∣5343∣∣=ln∣∣e2k∣∣
−0.2090917979=2k despejando obtenemos el valor de k
k=−0.1045458989
Reemplazando este valor en la ecuación (2) obtenemos
T=17+53e−0.1045458989t(3).
Ahora vamos a calcular cuánto tiempo se demora la taza de café para llegar a 20°C
20=17+53e−0.1045458989t⇒t=27,46.
La taza de café se demora 27,46 minutos para llegar a 20∘C
162
Ejemplo 2: Temperatura
Un termómetro que indica 70∘F se coloca en un horno precalentado a temperatura constante. A través de una ventana de vidrio del horno, un observador registra que la temperatura es de 110∘F después de medio minuto y de 145∘F después de un minuto. ¿A qué temperatura está el horno?
Solución
Se plantean las condiciones iniciales del problema.
Condiciones iniciales.
t=0
t=0.5
t=1
T=70∘F
T=110∘F
T=145∘F
Tm=?
Vamos a comenzar reemplazando el primer dato T=70∘C cuando t=0 en la ecuación T=Tm+cekt(1).
70=Tm+cek(0)⇒c=70−Tm
163
Reemplazando este valor en la ecuación (1) queda de la siguiente forma T=Tm+(70−Tm)ekt(2).
Ahora reemplazamos T=110 cuando t=0.5 en la ecuación (2)
110=Tm+(70−Tm)ek(0.5)(3) Nos queda una ecuación con dos incógnitas, como no se puede resolver una ecuación con dos incógnitas, entonces, reemplazamos la otra condición para tener una segunda ecuación.
145=Tm+(70−Tm)ek(1)(4)
Despejando Euler de la ecuación (3)e0.5k=70−Tm110−Tm(5)
Despejando Euler de la ecuación (4)ek=70−Tm145−Tm(6)
No podemos igualar las ecuaciones (5)∧(6)
Aplicamos axioma de igualdad a la ecuación (5) elevando al cuadrado a ambos lados
Cuando se resuelven ejercicios de crecimiento y decrecimiento o de temperaturas y piden población inicial, temperatura inicial o temperatura del medio, siempre se debe plantear un sistema de ecuaciones para poder resolver el ejercicio.
165
Ejemplo 3: Temperatura
Un termómetro se lleva del interior de una habitación donde la temperatura del aire es de 80∘F al exterior donde la temperatura del aire es de 15∘F. después de un minuto el termómetro indica 55∘F. ¿Cuál es la temperatura del termómetro en t=2 minutosS? ¿Cuánto tiempo le tomará alcanzar 19∘F
Solución Se plantean las condiciones iniciales del problema.
Condiciones iniciales.
t=0
t=1
t=2
t=?
T=80∘F
T=55∘F
T=?
T=19∘F
Tm=15∘F
Vamos a comenzar reemplazando el primer dato T=80∘C cuando t=0 en la ecuación T=Tm+cekt(1).
80=15+cek(0)⇒c=65
166
Reemplazando este valor en la ecuación (1) queda de la siguiente forma T=15+65ekt(2).
Ahora reemplazamos T=55 cuando t=1 en la ecuación (2)
55=15+65ek(1) despejando queda 6540=ek
Aplicamos logaritmo natural a ambos lados de la ecuación
ln∣∣6540∣∣=ln∣∣ek∣∣ obtenemos el valor de k
k=−0.4855078
Reemplazando este valor en la ecuación (2) obtenemos
T=15+65e−0.4855078t(3).
Ahora vamos a calcular la temperartura a los 2 minutos
T=15+65e(−0.4855078)(2)⇒T=39.61.
La temperatura a los dos minutos es 39.61∘F
Ahora vamos a calcular cuánto tiempo se demora el termómetro para llegar a 19∘F
19=15+65e−0.4855078t $
ln∣∣654∣∣=ln∣∣e−0.4855078t∣∣⇒t=5,74
167
El termómetro se demora 5,74 minutos para llegar a 19∘F
Ejemplo 4: Temperatura
Se sumerge una barra de acero en un líquido que está a una temperatura de 110∘C, la barra está a 32∘C, un segundo después su temperatura aumenta 4∘C. ¿Cuál será la temperatura de la barra a los 10 segundos? ¿Cuánto tiempo le tomará alcanzar 100∘C
Solución
Tenga en cuenta que la temperatura de la barra aumenta 4∘C en un segundo, es decir en t=1 la temperatura es 36∘C.
Condiciones iniciales.
t=0
t=1
t=10
t=?
T=32∘C
T=36∘C
T=?
T=100∘C
Tm=110∘C
168
Vamos a comenzar reemplazando el primer dato T=32∘C cuando t=0 en la ecuación T=Tm+cekt(1).
32=110+cek(0)⇒c=−78
Reemplazando este valor en la ecuación (1) queda de la siguiente forma T=110−78ekt(2).
Ahora reemplazamos T=36 cuando t=1 en la ecuación (2)
36=110−78ek(1) despejando queda −78−74=ek
Aplicamos logaritmo natural a ambos lados de la ecuación
ln∣∣7874∣∣=ln∣∣ek∣∣ obtenemos el valor de k
k=−0.05264373
Reemplazando este valor en la ecuación (2) obtenemos
T=110−78e−0.05264373t(3).
Ahora vamos a calcular la temperatura a los 10 segundos
T=110−78e(−0.05264373)(10)⇒T=63.92.
La temperatura a los 10 segundos es 63.92∘C
Ahora vamos a calcular cuánto tiempo se demora el termómetro para llegar a 100∘C
169
100=110−78e−0.05264373t despejando queda
−78−10=e−0.05264373t
Aplicamos logaritmo natural a ambos lados de la ecuación
ln∣∣7810∣∣=ln∣∣e−0.05264373t∣∣⇒t=39.01
El termómetro se demora 39.01 segundos para llegar a 100∘F
Ejemplo 5: Temperatura
Se encuentra el cuerpo de un hombre una mañana fría en una habitación cuya temperatura es 15∘C, un detective mide la temperatura del cuerpo y el termómetro indica 36∘C, una hora después vuelve a medir la temperatura del cuerpo y el termómetro indica 35.2∘C. En qué momento ocurrió el asesinato.
[Sugerencia: La temperatura normal del cuerpo es de 37∘C].
Solución
Condiciones iniciales.
170
t=0
t=t
t=t+1
T=37∘C
T=33∘C
T=32.1
Tm=15∘C
Vamos a comenzar reemplazando el primer dato T=37∘C cuando t=0 en la ecuación T=Tm+cekt(1).
37=15+cek(0)⇒c=22
La ecuación (1) queda de la siguiente forma T=15+22ekt(2).
Ahora reemplazamos T=33 cuando t=t en la ecuación (2)
33=15+22ek(t) despejando queda 2218=ekt
Aplicamos logaritmo natural a ambos lados de la ecuación
ln∣∣2218∣∣=ln∣∣ekt∣∣ obtenemos una ecuación con dos incógnitas$
kt=−0.20067079(3) Nos queda una ecuación con dos incógnitas, como no se puede resolver, reemplazamos la otra condición para tener una segunda ecuación.
Aplicamos logaritmo natural a ambos lados de la ecuación
171
ln∣∣2217.1∣∣=ln∣∣ek(t+1)∣∣ obtenemos
k(t+1)=−0.25196398(4)
Despejamos k de las ecuaciones (3)
k=t−0.20067079(5)
Despejamos k de las ecuaciones (4)
k=t+1−0.25196398(6)
Igualamos (5)∧(6)
t−0.20067069=t+1−0.25196398
Operando obtenemos
−0.20067069t−0.20067069=−0.25196398t
Agrupando términos semejantes
0.05129329t=0.20067069 despejando t=3.912
3.912 horas equivale a 235 minutos
El asesinato ocurrió a las 8:06 am
En la página 199 puedes observar un vídeo de aplicación de temperaturas.
172
2.7.4 Ejercicios y respuestas de la sección 2.7.3
173
2.7.5 Mezclas
las soluciones químicas, entendiéndose estas como un sistema homogéneo donde hay mayor proporción de soluto que de solvente, conllevan la necesidad de comprender el concepto de concentración, el cual no es más que la relación entre la cantidad de soluto y la cantidad de solvente. Las ED son una herramienta para comprender el comportamiento de las mezclas de dos o más soluciones en un recipiente bajo ciertas condiciones
Se tomará la sal como soluto para definir la ED, pero es importante resaltar que es igualmente aplicable a todo tipo de soluciones homogéneas, pues la ecuación mide la variación de la cantidad de soluto presente en la mezcla en cualquier momento t
Al mezclar dos soluciones salinas de distintas concentraciones surge una ecuación diferencial de primer orden, que define la cantidad de sal contenida en la mezcla. A(t) denota la cantidad de sal (medida en libras) en el tanque en el tiempo t, entonces la razón con la que A(t) cambia es una razón neta:
dtdA=R1−R2
La rapidez está dada por
R1= Velocidad de entrada por concentración
R2= Velocidad de salida por concentración
La rapidez con la que la sal entra es minlbs
174
Ejemplo 1: Mezclas
Se disuelve inicialmente 50 libras de sal en un tanque con 300 galones de agua. Se bombea Salmuera al tanque a razón 3mingal y luego la solución adecuadamente mezclada se bombea fuera del tanque a razón de 3mingal. Si la concentración de la solución que entra es de 2gallib, determine la cantidad de sal que hay en el tanque en un instante cualquiera. ¿Cuánta sal hay después de 50 min? ¿Cuánta después de un largo tiempo?
Solución
La ilustración ayuda a comprender lo que está pasando en el tanque.
En el tanque está entrando salmuera a razón de 3mingal con una concentración de 2gallib. Es decir, R1=3mingal.2gallib=6minlib
Inicialmente en el tanque hay 300 galones de salmuera, entonces para hallar R2 se plantea la ecuación R2=re.V+(re−rs)tA.
Reemplazando datos en la ecuación.
175
R2=3mingal.300+(3−3)tAgallib=100Aminlib.
Reemplazando estos datos en la ED dtdA=R1−R2 nos queda
dtdA=6−100A organizando la ED
dtdA+1001A=6
Tenemos una ED lineal no homogénea y está en forma estándar. Identificamos p(t)∧ hallamos el factor integrante.
p(t)=1001F.I=e100t
Aplicamos la propiedad para resolverla.
F.Iy=∫F.If(x)dx
e100tA=∫6e100tdt
Integrando el lado derecho
e100tA=600e100t+c
Despejamos la variable dependiente A
A=600+ce−100t
Ahora reemplazamos condiciones iniciales A(0)=50 para hallar la c
50=600+ce0
El valor de la constante es
176
c=−550
La solución es
A(t)=600−550e−100t
Teniendo la ecuación de concentración de sal, podemos responder las dos preguntas
Para saber la cantidad de sal a los 50 minutos se reemplaza, la t=50
A(50)=600−550e−10050
A(50)=266.4
La cantidad de sal que hay a los 50 minutos es 266.4 libras
Ahora necesitamos saber la cantidad de sal después de un largo tiempo, es decir, t=∞
A(∞)=600−550e−100∞
A(∞)=600
La cantidad de sal que hay despues de un largo tiempo es 600 libras
177
Ejemplo 2: Mezclas
Un tanque contiene 450 litros de agua en el que se disuelven 30 gramos de sal: Una salmuera que contiene 3lgr se bombea al tanque con una intensidad de 6minl, la solución adecuadamente mezclada se bombea hacia afuera con una intensidad de 8minl. Determine la cantidad de libras de sal que hay en el tanque en cualquier instante. ¿Cuándo se vacía el tanque?
Solución
En el tanque está entrando salmuera a razón de 6minl con una
concentración de 3lgr. Es decir, R1=6minl.3lgr=18mingr
Inicialmente en el tanque hay 450 litros de salmuera, entonces para hallar R2 se plantea la ecuación R2=re.V+(re−rs)tA.
Reemplazando datos en la ecuación.
R2=8minl.450+(6−8)tAgallib=8450−2tAmingr.
Reemplazando estos datos en la ED dtdA=R1−R2 nos queda
dtdA=18−4225−tA organizando la ED
178
dtdA+225−t4A=18
Tenemos una ED lineal no homogénea y está en forma estándar. Identificamos p(t)∧ hallamos el factor integrante.
p(t)=225−t4F.I=(225−t)−4
Ahora aplicamos la propiedad para resolverla
(225−t)−4A=∫18(225−t)−4dt
Integrando el lado derecho de la ecuación
(225−t)−4A=6(225−t)−3+c
Despejamos la variable dependiente A
A=6(225−t)+c(225−t)4
Ahora reemplazamos condiciones iniciales A(0)=30 para hallar la c
30=6(225)+c(225)4
El valor de la constante es
c=−(225)41320
La solución es
A=6(225−t)−1320(225225−t)4
179
La ecuación que nos ayuda a determinar la cantidad de sal que hay en el tanque en cualquier momento es
A(t)=6(225−t)−1320(225225−t)4
Ahora necesitamos saber cuándo se vacía el tanque. Como la intensidad de salida es mayor que la intensidad de entrada, es lógico que llega un momento en que el tanque estará vacío
Utilizamos la condición 450+(6−8)t. entonces es necesario igualar esta ecuación a cero y despejar t
El tanque se vacía en t=225 minutos
Ejemplo 3: Mezclas
Un estanque contiene 100m3 de agua contaminada. Con el propósito de descontaminarlo se introduce agua limpia a razón de 2minm3 y el agua contaminada (uniformemente mezclada) se deja salir del estanque a la misma razón. ¿Qué porcentaje de contaminantes se habrá eliminado después de 1 h? ¿Qué tiempo debe transcurrir para que los contaminantes disminuyan en un 90%?
Solución
180
La ilustración ayuda a comprender lo que está pasando en el tanque.
En el tanque está entrando agua limpia a razón de 2minm3 con una concentración de 0 contaminación. Es decir, R1=2minm3.0=0
Inicialmente en el tanque hay 100m3 de agua contaminada, entonces para hallar R2 se plantea la ecuación R2=re.V+(re−rs)tA.
Reemplazando datos en la ecuación.
R2=2.100+(2−2)tA=50Aminm3.
Reemplazando estos datos en la ED dtdA=R1−R2 nos queda
dtdA=0−50A organizando la ED
dtdA+501A=0
ED lineal homogénea en forma estándar. Identificamos p(t)
p(t)=501F.I=e50t
Ahora aplicamos la propiedad para resolverla
181
e50tA=∫0dt
Integrando obtenemos
e50tA=c
Despejamos la variable dependiente A
A=ce50t
Ahora reemplazamos condiciones iniciales A(0)=A0 para hallar la c
A0=ce500
El valor de la constante es
c=A0
La solución es
A(t)=A0e50t
Teniendo la ecuación de concentración de contaminante, podemos responder las dos preguntas
Para saber el porcentaje de contaminante que se habrá eliminado después de 60 minutos
A(t)=A0e5060
La cantidad de contaminante que hay después de 1 hora es 0.3012A0
182
Ahora vamos a calcular qué porcentaje es esta cantidad de A0
A0A0−0.3012A0=A0A0(1−0.3012)=0.6989
Después de una hora se habrá eliminado el 69.89% del contaminante del estanque
Ejemplo 4: Mezclas
El aire del interior de un pequeño cuarto con dimensiones de 12x8x8 metros contiene 3% de monóxido de carbono. Empezando en t=0, se sopla aire fresco que no contiene monóxido de carbono hacía el interior del cuarto a razón de 100minm3. Si el aire del cuarto sale al exterior a través de una abertura a la misma velocidad.
¿Cuándo tendrá el aire del interior del cuarto 0.01% de monóxido de carbono?
Solución
La cantidad de monóxido de carbono que hay inicialmente en el cuarto es el 3% del volumen V=768m3
A(0)=100768x3=23.04
183
En el cuarto entra aire limpio a razón de 100minm3. Es decir, R1=100minm3.0=0
Inicialmente en el cuarto hay 768m3 litros de salmuera, entonces para hallar R2 se plantea la ecuación R2=re.V+(re−rs)tA.
Reemplazando datos en la ecuación.
R2=100.768+(100−100)tA=768100A=19225A.
Reemplazando estos datos en la ED dtdA=R1−R2 nos queda
dtdA=0−19225A organizando la ED
dtdA+19225A=0
ED lineal homogénea en forma estándar. Identificamos p(t)∧ hallamos el factor integrante.
p(t)=19225F.I=e19225t
Ahora aplicamos la propiedad para resolverla
e19225tA=∫0dt
Integrando el lado derecho de la ecuación
e19225tA=c
Despejamos A
184
A=ce−19225t
Ahora reemplazamos condiciones iniciales A(0)=23.04 para hallar la c
23.04=ce−19225(0)
El valor de la constante es
c=23.04
La solución es
La ecuación que determina la concentración de monóxido de carbono en función del tiempo es:
A=23.04e−19225t
Ahora podemos responder la pregunta. ¿Cuándo tendrá el aire del interior del cuarto 0.01% de monóxido de carbono?
0.0768=23.04e−19225t
ln∣∣23.040.0768∣∣=ln∣∣e−19225t∣∣ obtenemos
t=43.805
En un tiempo de t=43.8 minutos el aire del interior del cuarto tendrá 0.01% de monóxido de carbono
En la página 200 puedes observar un vídeo de aplicación de mezclas.
185
2.7.6 Ejercicios y respuestas de la sección 2.7.5
186
2.7.7 Circuitos en serie
Circuitos LR: Cuando un circuito en serie sólo contiene un inductor (L) y un resistor (R), la segunda ley de Kirchhoff establece que:
la suma de las caídas de voltaje a través de un inductor Ldtdi y de un resistor (iR) es igual al voltaje aplicado E(t) al circuito
Circuito LR en serie
Ldtdi+Ri=E(t)(1)
Circuitos RC: Cuando un circuito en serie sólo contiene un resistor (R) y un capacitor (C), la caída de voltaje a través de un capacitor de capacitancia (C) es
cQ(t) donde q es la carga en el capacitor.
Por la segunda ley de kirchhoff Ri+Ci+q=E(t)(2)
La corriente (i) y la carga (q) se relacionan mediante idtdq.
Así la ecuación (2) se transforma en
187
Circuito RC en serie
Rdtdq+C1q=E(t)(3)
Ejemplo 1: Circuito en serie LR
Se conecta una bateria de 20 Voltios a un circuito en serie LR con 0.2 henry de inductancia y 30 Ohms de resistencia. Determine la corriente i(t) sí i(0)=0.
Solución
Recordemos que
Ldtdi+Ri=E(t)
Reemplazando valores tenemos 0.2dtdi+30i=20
Es una ED lineal pero no está en forma estándar.
188
LLevandola a la forma estándar, obtenemos.
dtdi+150i=100
Identificamos p(t) y hallamos el FI.
p(t)=150F.I=e150t
Ahora aplicamos la propiedad para resolverla
e150ti=∫100e150tdt
Integrando el lado derecho de la ecuación
e150ti=150100e150t+c
Despejamos la variable dependiente i
i=32+ce−150t
Ahora reemplazamos condiciones iniciales i(0)=0 para hallar la c
0=32+ce0
El valor de la constante es
c=−32
La corriente en función del tiempo es
i(t)=32−32e−150t
189
Ejemplo 2: Circuito en serie RC
Se aplica una fuerza electromotriz de 100 Voltios a un circuito en serie RC, donde la resistencia es de 200 Ohms y la capacitancia es de 10−4 farads. Determine la carga q(t) en el capacitor, Si q(0)=0. Encuentre la corriente i(t)
Solución
Recordemos que
Rdtdq+C1q=E(t)
Reemplazando valores tenemos
200dtdq+10−41q=100
Es una ED lineal pero no está en forma estándar.
LLevandola a la forma estándar, obtenemos.
dtdq+50q=21
Identificamos p(t) y hallamos el FI.
190
p(t)=50F.I=e50t
Ahora aplicamos la propiedad para resolverla
e50tq=∫21e50tdt
Integrando el lado derecho de la ecuación
e50tq=1001e50t+c
Despejamos la variable dependiente q
q=1001+ce−50t
Ahora reemplazamos condiciones iniciales q(0)=0 para hallar la c
c=−1001
La carga en función del tiempo es
q(t)=1001−1001e−50t
La corriente es la derivada de la carga, por lo cual debemos derivar la respuesta anterior
q′(t)=10050e−50t
Entonces la corriente en función del tiempo es:
i(t)=21e−50t
191
Ejemplo 3: Circuito en serie LR
Se aplica una fuerza electromotriz de 450 voltios a un circuito LR con 0.25 henry de inductancia y 90 ohms de resistencia. Determine la corriente i(t) sí i(0)=0. Cuando t tiende al infinito.
Solución
Recordemos que Ldtdi+Ri=E(t)
Reemplazando valores tenemos 0.25dtdi+90i=450
Es una ED lineal pero no está en forma estándar.
LLevandola a la forma estándar, obtenemos.
dtdi+360i=1800
Identificamos p(t) y hallamos el FI.
p(t)=360F.I=e360t
Ahora aplicamos la propiedad para resolverla
192
e360ti=∫1800e360tdt
Integrando el lado derecho de la ecuación
e360ti=3601800e360t+c
Despejamos la variable dependiente i
i=e360t5e360t+e360tc
Simplificando
i=5+ce−360t
Ahora reemplazamos condiciones iniciales i(0)=0 para hallar la c
0=5+ce0
El valor de la constante es
c=−5
La corriente en función del tiempo es
i(t)=5−5e−360t
La corriente cuando t tiende a infinito es
i(t)=5−5e(−360)(∞)
i(∞)=5
193
Ejemplo 4: Circuito en serie RC
Se aplica una fuerza electromotriz de 200cos2t Voltios a un circuito en serie RC, donde la resistencia es de 50 Ohms y la capacitancia es de 10−2 farads. Determine la carga q(t) en el capacitor, Si q(0)=0. Encuentre la corriente i(t)
Solución
Recordemos que
Rdtdq+C1q=E(t)
Reemplazando valores tenemos
50dtdq+10−21q=200cos2t
50dtdq+100q=200cos2t
Es una ED lineal pero no está en forma estándar.
194
LLevándola a la forma estándar, obtenemos.
dtdq+2q=4cos2t
Identificamos p(t) y hallamos el FI.
p(t)=2F.I=e2t
Ahora aplicamos la propiedad para resolverla
e2tq=4∫e2tcos2tdt
Utilizamos el siguiente resultado de la tabla de integrales
∫eaxcosbxdx=a2+b2eax(acosbx+bsenbx)+c
∫e2tcos2tdt=4+4e2t(2cos2t+2sen2t)+c
4∫e2tcos2tdt=4(8e2t)2(cos2t+sen2t)+c
Simplificando
4∫e2tcos2tdt=e2t(cos2t+sen2t)+c
Reemplazando este resultado en la ED
e2tq=e2t(cos2t+sen2t)+c
Despejamos la variable dependiente q
195
q=e2te2t(cos2t+sen2t)+e2tc
Simplificando
q=cos2t+sen2t+ce−2t
Ahora reemplazamos condiciones iniciales q(0)=0 para hallar la c
0=cos2(0)+sen2(0)+ce−2(0)
El valor de la constante es
c=−1
La carga en función del tiempo es
q=cos2t+sen2t−e−2t
Para calcular la corriente recuerde que I=dtdq
I=−2sen2t+2cos2t+2e−2t
La corriente en función del tiempo es
I=−2sen2t+2cos2t+2e−2t
En la página 200 puedes observar un vídeo de aplicación de circuitos.
196
En el siguiente video, realizado por Jaime H. Ramírez Rios, podrás observar un ejemplo de aplicación de crecimiento.
En el siguiente video, realizado por Jaime H. Ramírez Rios, podrás observar un ejemplo de aplicación de temperatura.
197
En el siguiente video, realizado por Jaime H. Ramírez Rios, podrás observar un ejemplo de aplicación de mezclas.
En el siguiente video, realizado por Jaime H. Ramírez Rios, podrás observar un ejemplo de aplicación de circuitos.
198
2.7.8 Ejercicios y respuestas de la sección 2.7.7
199
2.7.9 Ejercicios y respuestas del capítulo 2
200
201
202
Capítulo 3
ECUACIONES DIFERENCIALES DE ORDEN SUPERIOR
204
3.1 Principio de superposición
En la segunda unidad se estudiaron ED de primer orden. En esta unidad se estudiarán las ecuaciones de orden superior n≥2 comenzando con las ecuaciones lineales.
Una ED lineal de orden superior es de la forma:
any(n)+an−1y(n−1)+…+a2y′′+a1y′+a0y=g(x)
Cuando g(x)=0 se dice que la ED es homogénea, en caso contrario es no homogénea.
Soluciones
La función y=h(x) se llama solución de la ED si está definida y es derivable n veces en algún intervalo de tal manera que al sustituirla en la ecuación junto con sus derivadas se obtenga una identidad.
Teorema 1. Principio de superposición; ecuaciones homogéneas. Sean y1,y2,y3,…yn soluciones de la ecuación homogénea de n-ésimo orden en un intervalo I. entonces la combinación lineal
y=c1y1(x)+c2y2(x)+…+cnyn(x)
Donde c1,c2,…,cn son constantes arbitrarias, también es una solución en el intervalo.
205
Ejemplo 1. Principio de superposición
Usando el principio de superposición, probar que las funciones
y1=c1e−x, y2=c2xe−x son soluciones de la ED y′′+2y′+y=0
Solución
Escribimos la solución como una combinación lineal de las soluciones.
Como se cumple la identidad, se prueba por el principio de superposición que las funciones
y1=c1e−x∧y2=c2xe−x son solución de la ED y′′+2y′+y=0
206
3.1.1 Ejercicios y respuestas de la sección 3.1
207
3.2 Dependencia e independencia lineal
Definición se dice que un conjunto de funciones f1(x),f2(x),…,fn(x) es linealmente dependiente en un intervalo I si existen constantes c1,c2,…,cn no todas ceros, tales que
c1f1(x)+c2f2(x)+…+cnfn(x)=0
Para toda x en el intervalo. Si el conjunto de funciones no es linealmente dependiente en el intervalo, se dice que es linealmente independiente.
Definición del Wronskiano Suponga que cada una de las funciones f1(x),f2(x),…,fn(x) tiene al menos n−1 derivadas. El determinante
Donde las primas denotan derivadas, se llama el wronskiano de las funciones.
208
Teorema 2. Criterio para soluciones linealmente independientes. Sean y1,y2,y3,…ynn soluciones de la ecuación diferencial lineal homogénea de n-ésimo orden en un intervalo I. El conjunto de soluciones es linealmente independiente en I si y sólo si W(y1,y2,y3,…yn)=0 para toda x en el intervalo.
Ejemplo 1. Dependencia e independencia lineal
Verificar si el conjunto de funciones f1(x)=e3x,f2(x)=5e3x es linealmente dependiente (LD) o linealmente independiente (LI) en el intervalo (−∞,∞)
Solución
Para determinar si el conjunto de funciones es LI o LD se deben analizar las dos funciones.
Si multiplicamos f1 por 5, obtenemos f2 , es decir, 5f1=f2
Entonces f1 es multiplo escalar de f2
Si dos funciones son multiplo escalar, son LD
Las funciones son linealmente dependientes.
209
Otra forma de resolver el ejercicio es aplicando el Wronskian de las funciones.
Verificar si el conjunto de funciones f1(x)=x,f2(x)=x2∧f3(x)=3x2−2x es LD o LI en el intervalo (−∞,∞)
Solución
Para determinar si el conjunto de funciones es LI o LD se deben analizar las tres funciones.
La función f3(x) se puede escribir como una combinación lineal de las otras dos funciones f3(x)=3f2(x)−2f1(x), entonces,
211
Las funciones son linealmente dependientes.
Nota un conjunto de funciones f1(x),f2(x),…,fn(x) es LD en un intervalo, si por lo menos una función se puede expresar como una combinación lineal de las otras funciones.
Ahora lo vamos a resolver aplicando el Wronskian de las funciones.
En la página 231 puedes observar un vídeo de dependencia e independencia lineal
Como W=0 las funciones son LI.
213
Repasando Determinantes
Asociado a cada matriz cuadrada A hay un número llamado determinante de A denotado como “detA” o también por ∣A∣ (no confundir con el valor absoluto). Los determinantes nos proporcionan un método para el cálculo de la matriz inversa (en caso de existir) y un criterio para estudiar si una matriz es o no invertible. Sus aplicaciones son múltiples en todas las ramas de las ciencias que tratan problemas lineales en los que necesariamente aparecen matrices y por tanto, determinantes. Sólo se puede calcular el determinante si la matriz es cuadrada o de orden n.
En la escena interactiva que aparece a continuación se explica, el cálculo del determinante de matrices de orden 2 y 3. Es importante que repase el concepto antes de resolver ejercicios.
214
3.2.1 Ejercicios y respuestas de la sección 3.2
215
3.3 Conjunto fundamental de soluciones
Definición: Conjunto fundamental de soluciones Cualquier conjunto y1,y2,y3,…yn de n soluciones linealmente independientes de la ecuación diferencial lineal homogénea de n-ésimo orden en un intervalo I es un conjunto fundamental de soluciones (CFS) en el intervalo.
Teorema 3: solución general; ecuaciones homogéneas. Sea y1,y2,y3,…yn un CFS de la ED lineal homogénea de n-ésimo orden en un intervalo I. Entonces la solución general de la ecuación en el intervalo es:
y=c1y1(x)+c2y2(x)+…+cnyn(x)
Donde c1,c2,…,cn, son constantes arbitrarias
Ejemplo 1. Conjunto fundamental de soluciones
Compruebe que las funciones e−3x,e4x forman un CFS de la ecuación diferencial y′′−y′−12y=0 en el intervalo (−∞,∞).
216
Forme la solución general y demuestre que son la solución de la ED.
Solución
Para demostrar que las funciones forman un CFS, se debe demostrar que las funciones son LI.
Al demostrar que el lado derecho de la ecuación es igual al lado izquierdo, se verifica que la función
y=c1cos(2ln∣x∣)+c2sen(2ln∣x∣) es solución de la ED x2y′′+xy′+4y=0
Cuando un conjunto de funciones es LD, significa que esas funciones no son la solución de ninguna ED, sólo se verifica si son la solución de un ED cuando son LI, es decir, si forman un CFS.
220
3.3.1 Ejercicios y respuestas de la sección 3.3
221
3.4 Reducción de orden
En la sección anterior vimos que la solución general de una ecuación diferencial lineal homogénea de segundo orden
a2y′′+a1y′+a0y=0(E1)
Es una combinación lineal y=c1y1(x)+c2y2(x) donde y1∧y2 son soluciones linealmente independientes en un intervalo I.
En esta sección vamos a ver un método para determinar dichas soluciones cuando los coeficientes de la ED en (E1) son constantes
REDUCCIÓN DE ORDEN Suponga que y1 denota una solución no trivial de a2y′′+a1y′+a0y=0 y que y1 se define en un intervalo I.
Se busca una segunda solución y2 tal que y1∧y2 sean un conjunto LI en I.
CASO GENERAL dividiendo la ecuación (E1) por a2(x) queda la ecuación en la forma estándar.
y′′+P(x)y′+Q(x)y=0
Para encontrar la segunda solución se reemplaza en la ecuación:
222
y2(x)=y1(x)∫[y1(x)]2e−∫p(x)dxdx(E2)
Ejemplo 1. Reducción de orden
La función y1=e−3x es una solución de la ED y′′−y′−12y=0. Use reducción de orden para encontrar una segunda solución y2. Forme la solución General.
Solución
Este ejercicio lo resolvimos como un CFS en la página 180, demostramos que las soluciones e−3x,e4x forman un CFS, luego demostramos que son la solución de la ED y′′−y′−12y=0. Ahora vamos a comenzar con la solución e−3x y utilizando reducción de orden vamos a hallar la segunda solución e4x
Como la ED está en forma estándar identificamos p(x) que es el término que acompaña a y′
p(x)=−1
El p(x) y la solución y1=e−3x los reemplazamos en la ecuación (E2) para hallar la segunda solución y2
223
y2=e−3x∫[e−3x]2e∫dxdx
Integramos el numerado y operamos el denominador
y2=e−3x∫e−6xedxdx
Aplicando potenciación
y2=e−3x∫e7xdx
Integrando
y2=e−3x(71e7x)
Operando
y2=71e4x
Solución general
y=c1e−3x+c2e4x
Ejemplo 2. Reducción de orden
La función y1=sen3x es una solución de la ED y′′+9y=0. Use reducción de orden para encontrar una segunda solución y2. Forme la solución General.
224
Solución
En la ED no hay y′, eso quiere decir que el número que acompaña este término es 0
p(x)=0
El p(x) y la solución y1=sen3x los reemplazamos en la ecuación (E2) para hallar la segunda solución y2
y2=sen3x∫[sen3x]2e∫0dxdx
La integral queda de la forma
y2=sen3x∫sen23x1dx
Identidad
y2=sen3x∫csc23xdx
Integrando
y2=sen3x(−31cot3x)
Identidad
y2=−31sen3x(sen3xcos3x)
Simplificando
y2=−31cos3x
225
Solución general
y=c1sen3x+c2cos3x
Ejemplo 3. Reducción de orden
La función y1=e32x es una solución de la ED 9y′′−12y′+4y=0. Use reducción de orden para encontrar una segunda solución y2.
Forme la solución General.
Solución
A diferencia de los ejemplos anteriores, esta ED no está en forma estándar, entonces debemos dividir toda la ecuación por 9 antes de identificar p(x)
99y′′−912y′+94y=90
Simplificando
y′′−34y′+94y=0
Identificamos p(x)
p(x)=34
226
El p(x) y la solución y1=e32x los reemplazamos en la ecuación (E2) para hallar la segunda solución y2
y2=e32x∫[e32x]2e∫34xdxdx
La integral queda de la forma
y2=e32x∫[e34x]e34xdx
Simplificando
y2=e32x∫dx
Integrando
y2=e32xx
La segunda solución es
y2=xe32x
Solución general
y=c1e32x+c2xe32x
227
Ejemplo 4. Reducción de orden
La función y1=xcos(ln∣x∣) es una solución de la ED x2y′′−xy′+2y=0. Use reducción de orden para encontrar una segunda solución y2. Forme la solución General.
Solución
Como la ED no está en forma estándar, debemos dividir toda la ecuación por x2 para identificar p(x)
x2x2y′′−x2xy′+x22y=x20
Simplificando
y′′−x1y′+x22y=0
Identificamos p(x)
p(x)=−x1
El p(x) y la solución y1=xcos(ln∣x∣) los reemplazamos en la ecuación (E2) para hallar la segunda solución y2
y2=xcos(ln∣x∣)∫[xcos(ln∣x∣)]2e∫x1dxdx
La integral queda de la forma
228
y2=xcos(ln∣x∣)∫x2cos2(ln∣x∣)xdx
Simplificando y aplicando identidades
y2=xcos(ln∣x∣)∫xsec2(ln∣x∣)dx
Esta integral se resuelve por sustitución, siendo u=(ln∣x∣)∧du=x1dx
y2=xcos(ln∣x∣)∫sec2udu
Integrando
y2=xcos(ln∣x∣)tanu
Identidad
y2=xcos(ln∣x∣)(cosusenu)
Sustituyendo u
y2=xcos(ln∣x∣)(cos(ln∣x∣)sen(ln∣x∣))
Simplificando obtenemos la segunda solución
y2=xsen(ln∣x∣)
Solución general
y=c1xcos(ln∣x∣)+c2xsen(ln∣x∣)
229
En el siguiente video, realizado por Jaime H. Ramírez Rios, podrás observar un ejemplo de dependencia e independencia lineal.
En el siguiente video, realizado por Jaime H. Ramírez Rios, podrás observar un ejemplo de reducción de orden.
230
3.4.1 Ejercicios y respuestas de la sección 3.4
231
3.5 Ecuaciones lineales homogéneas con coeficientes constantes
Se trata de determinar si existen soluciones exponenciales en (−∞,∞) de las ecuaciones lineales homogéneas de orden superior del tipo
any(n)+an−1y(n−1)+…+a2y′′+a1y′+a0y=0
Donde los coeficientes an,an−1,…,a2,a1,a0 son constantes reales ∧an=0
Método de solución: comenzamos para el caso especial de la ecuación de segundo orden
a2y′′+a1y′+a0y=0
La única función elemental no trivial cuya derivada es una constante múltiple de sí misma es la función exponencial emx, entonces la solución es de la forma y=emx
Vamos a reemplazar la solución y=emx y sus derivadas en la ED de segundo orden y vamos a cambiar a2,a1,a0, por a,b,c puesto que sólo necesitamos tres constantes
y=emxy′=memxy′′=m2emx
am2emx+bmemx+cemx=0
Factor común
emx(am2+bm+c)=0
Simplificando
232
am2+bm+c=0
Esta solución se conoce como ecuación auxiliar o ecuación característica. Existen tres formas de la solución general que corresponden a los tres casos:
Tabla 4. Soluciones ED homogéneas
Ejemplo 1. Coeficientes constantes
Resolver la ED 2y′′−5y′−3y=0.
Solución
Primero se lleva la ED a la ecuación auxiliar
233
2m2−5m−3=0
Se determinan los coeficientes
a=2,b=−5,c=−3
Se reemplazan en la ecuación cuadrática
m=2(2)−(−5)±(−5)2−4(2)(−3)
Operando
m=45±49
Como el discriminante es 49>0 obtenemos dos respuestas diferentes
m=45+7∧m=45−7
Simplificando
m1=3m2=−21
Las dos soluciones se reemplazan en las soluciones de la ED de la tabla 4
y1=e3xy2=e−21x
Solución general
y=c1e3x+c2e−21x
234
Ejemplo 2. Coeficientes constantes
Resolver la ED y′′+6y′+9y=0.
Solución
Primero se lleva la ED a la ecuación auxiliar
m2+6m+9=0
Se factoriza
m=2(1)−(6)±(6)2−4(1)(9)
Operando
m=2−6±0
Como el discriminante es 0 obtenemos dos respuestas iguales
m1=−3m2=−3
Las dos soluciones son
y1=e−3xy2=xe−3
Solución general
y=c1e−3x+c2xe−3x
235
Ejemplo 3. Coeficientes constantes
Resolver la ED y′′+4y′+7y=0.
Solución
Primero se lleva la ED a la ecuación auxiliar
m2+4m+7=0
Se factoriza
m=2(1)−(4)±(4)2−4(1)(7)
Operando
m=2−4±−12
Como el discriminante es −12<0 Multiplicamos por i2=−1
m=2−4±12i2
Obtenemos dos respuestas complejas
m=2−4±23i
Obtenemos la solución de la ecuación auxiliar
m=−2±3i
Los valores de α∧β son
236
α=−2β=3
Las soluciones utilizando la tabla 4 son
y1=e−2xcos3x
y2=e−2xsen3x
Solución general
y=c1e−2xcos3x+c2e−2xsen3x
Cuando las respuestas son complejas no se tiene en cuenta los signos más y menos, pues la función coseno es una función par y la función seno es impar, por tanto, no se colocan los signos, quedan de forma implícita, tampoco se coloca la i que indica imaginarios y que la respuesta es compleja.
Para resolver la ecuación auxiliar no es necesario utilizar la ecuación cuadrática, se pueden usar todas las técnicas de factorización
Ejemplo 4. Coeficientes constantes
Resolver la ED y′′−25y=0.
Solución
Primero se lleva la ED a la ecuación auxiliar
237
m2−25=0
Se factoriza como una diferencia de cuadrados
(m+5)(m−5)=0
Obtenemos las soluciones de la ecuación auxiliar
m1=5m2=−5
Las dos soluciones son
y1=e5xy2=e−5x
Solución general
y=c1e5x+c2e−5x
Ejemplo 5. Coeficientes constantes
Resolver la ED y′′+25y=0.
Solución
Primero se lleva la ED a la ecuación auxiliar
m2+25=0
La suma de cuadrados no se puede factorizar, entonces despejamos la variable
238
m=±−25
Obtenemos la solución de la ecuación auxiliar
m=0±5i
Los valores de α∧β son
α=0β=5
Las dos soluciones son
y1=e0cos5xy2=e0sen5x
Solución general
y=c1cos5x+c2sen5x
Ejemplo 6. Resolver el PVI
Resolver la ED y′′+y=0 Sujeta a y(3π)=0∧y′(3π)=2
Solución
Primero se lleva la ED a la ecuación auxiliar
m2+1=0
Despejamos la variable pues no es factorizable
239
m=±−1
Obtenemos la solución de la ecuación auxiliar
m=0±i
Los valores de α∧β son
α=0β=1
Las dos soluciones son
y1=e0cosxy2=e0senx
Solución general
y=c1cosx+c2senx
Como la ED tiene condiciones iniciales, debemos pasar de la solución general a la solución particular, para esto, reemplazamos las condiciones iniciales en la solución general
0=c1cos(3π)+c2sen(3π)
Evaluando obtenemos la ecuación (1)
0=21c1+23c2
Derivamos la solución general
y′=−c1senx+c2cosx
Reemplazamos condiciones iniciales
2=−c1sen(3π)+c2cos(3π)
Evaluando obtenemos la ecuación (2)
240
2=−23c1+21c2
Tenemos un sistema de ecuaciones.
⎩⎨⎧0=21c1+23c22=−23c1+21c2
Multiplicamos la ecuación (1) por 3 y resolvemos el sistema de ecuaciones por eliminación.
⎩⎨⎧0=23c1+23c22=−23c1+21c2
2=2c2
Obtenemos c2=1
Reemplazando este valor en la ecuación (1) llegamos a
0=21c1+23
Despejando obtenemos c1=−3
Entonces la solución particular es
y=−3cosx+senx
241
Ecuaciones de orden superior: para resolver una ecuación diferencial de n-ésimo orden de la forma
anm(n)+an−1m(n−1)+…+a2m2+a1m+a0=0
Si todas las raíces son reales y distintas, entonces se tiene la solución:
y=c1em1x+c2em2x+…+ccemnx
Si es de multiplicidad k, entonces se tiene la solución:
y=c1em1x+c2xem1x+c2x2em1x+…+ckxk−1em1x
Una ED de orden superior puede tener una combinación de soluciones diferentes, repetidas y complejas, para expresar la solución se utiliza la tabla 4
Ejemplo 7. Coeficientes constantes
Resolver la ED y′′′+3y′′−4y=0
Solución
Primero se lleva la ED a la ecuación auxiliar
m3+3m2−4=0
para factorizarla se hace división sintética
242
Con la división sintética obtenemos una raíz m=1 y con el residuo obtenemos una ecuación de segundo grado.
(m−1)(m2+4m+4)
Factorizando el trinomio cuadrado perfecto
(m−1)(m+2)2
Tenemos tres soluciones, una diferente y otra de multiplicidad dos.
y1=exy2=e−2xy3=xe−2x
Solución general
y=c1ex+c2e−2x+c3xe−2x
Ejemplo 8. Coeficientes constantes
Resolver la ED dx4d4y+2dx2d2y+y=0
Solución
Primero se lleva la ED a la ecuación auxiliar
243
m4+2m2+1=0
Es un polinomio de cuarto grado y se factoriza como un trinomio cuadrado perfecto
(m2+1)2=0
El exponente indica que la solución es de multiplicidad dos. Se soluciona la ecuación de segundo grado
m2+1=0
Al despejar y sacar raíz cuadrada obtenemos
m=±−i⇒m=±i
Tenemos una raíz compleja, buscamos en la tabla 4 y obtenemos dos soluciones
y1=e0cosxy2=e0senx
Como la solución es de multiplicidad dos
y1=cosxy2=senxy3=xcosxy4=xsenx
Solución general
y=c1cosx+c2senx+c3xcosx+c4xsenx
En la siguiente escena interactiva, podrás resolver ED de segundo orden por medio de coeficientes constantes.
244
En la siguiente escena interactiva, diseñada por Marco Peñaloza, podrás resolver ED de segundo orden homogéneas.
245
En el siguiente video, realizado por Jaime H. Ramírez Rios, podrás observar un ejemplo de coeficientes constantes.
En el siguiente video, realizado por Jaime H. Ramírez Rios, podrás observar un ejemplo de coeficientes constantes.
246
3.5.1 Ejercicios y respuestas de la sección 3.5
247
3.6 Ecuaciones lineales no homogéneas
Para resolver una ecuación diferencial lineal no homogénea
any(n)+an−1y(n−1)+…+a2y′′+a1y′+a0y=g(x)
Se debe hacer dos cosas
• Encontrar la función complementaria yc de la ecuación homogénea asociada
• Encontrar alguna solución particular yp de la ecuación no homogénea asociada
3.7 Coeficientes indeterminados
Es un camino directo y sencillo para determinar una solución de una ecuación lineal no homogénea con coeficientes constantes. El método consiste en hacer una conjetura acerca de la forma del yp, motivada por las clases de funciones que conforman g(x). El método general se limita a ED lineales donde
• Los coeficientes ai,i=0,1,2...,n son constantes y
• La función g(x) es una constante k, una función polinomial, una función exponencial eαx una función seno o cosenosenβx,cosβx o combinaciones de estas funciones
El método de coeficientes indeterminados no es aplicable a ecuaciones cuando la función g(x) corresponde a: ln∣x∣, x1, tanx, sen−1x, etc. Las ecuaciones en las que g(x) es una de estas funciones se consideran en la sección 3.8
248
$$
$$
Tabla 5. Soluciones particulares de prueba
Ejemplo 1. Coeficientes indeterminados
Resolver la ED y′′+4y′−2y=2x2−3x+6
Solución
Primero se debe resolver la ED homogénea
y′′+4y′−2y=0
Se lleva a la ecuación auxiliar
249
m2+4m−2=0
Se factoriza.
m=2−4±24
Obtenemos la solución de la ecuación auxiliar
m=−2±6
La solución de la ecuación homogénea es
yc=c1e(−2+6)x+c2e(−2−6)x
Factorizando
Ahora vamos a calcular yp, para eso buscamos en la tabla de la página 243 el yp de prueba para un polinomio de grado 2
yp=Ax2+Bx+C
Ahora vamos a reemplazar el yp y sus derivadas en la ED para calcular las constantes
yp′=2Ax+B
yp′′=2A
2A+4(2Ax+B)−2(Ax2+Bx+C)=2x2−3x+6
Se aplica propiedad distributiva
2A+8Ax+4B−2Ax2−2Bx−2C=2x2−3x+6
Con esta información vamos a construir un sistema de ecuaciones de tres ecuaciones con tres incógnitas
250
⎩⎨⎧2=−2A−3=8A−2B6=2A+4B−2C
Resolviendo el sistema de ecuaciones obtenemos los siguientes valores
A=−1B=−25c=−9
La solución particular es:
yp=−x2−25x−9
La solución general de la ED es:
y=c1e(−2+6)x+c2e(−2−6)x−x2−25x−9
Ejemplo 2. Coeficientes indeterminados
Resolver la ED y′′+9y=2sen3x
Solución
Primero se debe resolver la ED homogénea
y′′+9y=0
Se lleva a la ecuación auxiliar
251
m2+9=0
Despejamos la m
m=±3i
La solución de la ecuación homogénea es
yc=c1cos3x+c2sen3x
Para hallar el yp, buscamos en la tabla de la página 250 para la función seno
yp=Acos3x+Bsen3x
Derivamos y reemplazamos en la ED.
yp′=−3Asen3x+3Bcos3x
yp′′=−9Acos3x−9Bsen3x
−9Acos3x−9Bsen3x+9(Acos3x+Bsen3x)=2sen3x
Se aplica propiedad distributiva
−9Acos3x−9Bsen3x+9Acos3x+9Bsen3x)=2sen3x
Agrupando términos semejantes obtenemos
0=2sen3x
El resultado es contradictorio, lo que demuestra que yp no es el adecuado. Esto ocurre porque la solución complementaria yc ya tiene sen3x, entonces para solucionar el problema se debe multiplicar todo el yp por x
252
yp=Axcos3x+Bxsen3x
Ahora vamos a reemplazar el yp y sus derivadas en la ED para calcular las constantes
Con esta información vamos a construir un sistema de ecuaciones de dos ecuaciones con dos incógnitas
253
{2=−6A0=6B
Los valores de las incógnitas son:
A=−61B=0
La solución particular es:
yp=−61xcos3x
La solución general de la ED es:
y=c1cos3x+c2sen3x−61xcos3x
Ejemplo 3. Coeficientes indeterminados
Resolver la ED y′′−2y′−3y=4x−5+6xe2x
Solución
Primero se debe resolver la ED homogénea
y′′−2y′−3y=0
Se lleva a la ecuación auxiliar
254
m2−2m−3=0
Se factoriza
(m−3)(m+1)
Las raíces del polinomio son
m=3∧m=−1
La solución de la ecuación homogénea es
yc=c1e3x+c2e−x
Ahora vamos a calcular yp
Es importante anotar que aparecen dos yp, un yp1 que corresponde a un polinomio de grado 1 y un yp2 que corresponde a la función exponencial combinada con un polinomio de primer grado
El procedimiento utilizado para resolver una E.D. lineal de primer orden es también aplicable para resolver una E.D. de orden superior. Para adaptar el método de variación de parámetros a una ecuación diferencial de segundo orden
a2(x)y′′+a1(x)y′+a0(x)y=g(x)
Se debe llevar la ecuación a la forma estándar
y′′+P(x)y′+Q(x)y=f(x)
P(x)Q(x)∧f(x) con continuas en I. Se halla yc, la solución general de la ecuación homogénea. La solución particular para la ecuación lineal de segundo orden tiene la forma
yp=u1(x)y1(x)+u2(x)y2(x)
Donde y1∧y2 forman un conjunto fundamental de soluciones en I de la forma homogénea. Como la ecuación busca determinar dos funciones desconocidas u1∧u2 y se cuenta con una sola ecuación, se deriva dos veces yp y se sustituye en la ecuación en forma estándar obteniendo las ecuaciones:
y1u1+y2u2=0
y1′u1+y2′u2=f(x)
Este sistema de dos ecuaciones con dos incógnitas se puede resolver por la regla de Cramer
Las funciones u1∧u2 se obtienen integrando u1′∧u2′.
Ejemplo 1. Variación de parámetros
Resolver la ED y′′+2y′+y=x2e−x por variación de parámetros
Solución
La ecuación auxiliar es m2+2m+1=0
Factorizando (m+1)2=0
Las raíces son m=−1 raíz de multiplicidad dos
Solución complementaria yc=c1e−x+c2xe−x
Calculamos el Wronskiano W=[e−x−e−xxe−x−xe−x+e−x]
W=−xe−2x+e−2x−(−xe−2x)=−xe−2x+e−2x+−xe−2x
W=e−2x
280
Ahora calculamos W1=[0x2e−xx2e−x−xe−x+e−x]
W1=0−(x3e−2x)=−x3e−2x
Ahora calculamos W2=[e−x−e−x0x2e−x]
W2=x2e−2x−0=x2e−2x
Entonces u1′=e−2x−x3e−2x para calcular u1 se integra
∫−x3dx⇒u1=−4x4
u2′=e−2xx2e−2x para calcular u2 se integra
∫x2dx⇒u2=3x3 Recuerde que
yp=u1y1+u2y2
Entonces yp=−4x4e−x+3x3xe−x
Operando yp=−4x4e−x+3x4e−x=12x4e−x
y=c1e−x+c2xe−x+12x4e−x
281
No es necesario verificar si hay respuestas repetidas, el método le coloca el exponente a las respuestas repetidas. Observe que la respuesta del yp tiene exponente x4, significa que estaba repetida.
Ejemplo 2. Variación de parámetros
Resolver la ED y′′+25y=sen5x por variación de parámetros.
Solución
La ecuación auxiliar es m2+25=0
Como no se puede factorizar, se despeja m2=−25
Las raíces son m=±5i raíz compleja
Solución complementaria yc=c1cos5x+c2sen5x
Calculamos el Wronskiano W=[cos5x−5sen5xsen5x5cos5x]
W=5cos5x−(−5sen5x)=5cos5x+5sen5x
W=5
Ahora calculamos W1=[0sen5xsen5x5cos5x]
W1=0−(sen25x)=−sen25x
282
Ahora calculamos W2=[cos5x−5sen5x0sen5x]
W2=cos5xsen5x−0=sen5xcos5x
Entonces u1′=5−sen25x para calcular u1 se integra
Recuerde la identidad sen2x=21−sen(2x)
−51∫21−sen(10x)
⇒u1=−10x+1001cos10x
u2′=5sen5xcos5x para calcular u2 se integra por sustitución
Donde los coeficientes an,an−1,...,a2,a1,a0 son constantes se conoce como ecuación de Cauchy-Euler. Tiene como característica que el grado de los coeficientes nominales coincide con el orden de la ecuación. La ecuación de segundo orden tiene la forma:
ax2y′′+bxy′+cy=0
Método de solución: La función elemental no trivial cuya derivada es una constante múltiple de sí misma es la función exponencial xm, entonces la solución es de la forma y=xm
Tabla 6. Soluciones ecuación de Cauchy homogéneas
291
Ejemplo 1. ecuación homogénea
Resolver la ecuación de Cauchy Euler x2y′′+xy′−y=0
Solución
Vamos a reemplazar la solución y=xm y sus derivadas en la ED de segundo orden
y=xmy′=mxm−1y′′=m(m−1)xm−2
x2m(m−1)xm−2+xmxm−1−xm=0
Agrupamos las x aplicando propiedades de potenciación
m(m−1)xm+mxm−xm=0
Factor común
xm(m(m−1)+m−1)=0
Propiedad distributiva
m2−m+m−1=0
Agrupando términos semejantes
m2−1=0
Factorizando obtenemos
m1=1∧m2=−1
La solución general es
292
y=c1x+c2x−1
Ejemplo 2. ecuación homogénea
Resolver la ecuación de Cauchy Euler x2y′′+3xy′+y=0
Solución
Vamos a reemplazar la solución y=xm y sus derivadas en la ED de segundo orden
y=xmy′=mxm−1y′′=m(m−1)xm−2
x2m(m−1)xm−2+3xmxm−1+xm=0
Agrupamos las x aplicando propiedades de potenciación
m(m−1)xm+3mxm+xm=0
Factor común
xm(m(m−1)+3m+1)=0
Se pasa xm a dividir al lado derecho y se aplica propiedad distributiva
m2+2m+1=0
Factorizando
293
(m+1)2=0
Tenemos dos raíces iguales. m=−1. La solución general es
y=c1x−1+c2x−1ln∣x∣
Ejemplo 3. ecuación homogénea
Resolver la ecuación de Cauchy Euler x2y′′−xy′+2y=0
Solución
Vamos a reemplazar la solución y=xm y sus derivadas en la ED de segundo orden
y=xmy′=mxm−1y′′=m(m−1)xm−2
x2m(m−1)xm−2−xmxm−1+2xm=0
Agrupamos las x aplicando propiedades de potenciación
m(m−1)xm−mxm+2xm=0
Factor común
xm(m(m−1)−m+2)=0
294
Se pasa xm a dividir al lado derecho y se aplica propiedad distributiva
m2−2m+2=0
Como el polinomio no es factorizable, usamos la ecuación cuadrática
m=22±4−8
Obtenemos la solución de la ecuación auxiliar
m=22±−4
Raíz negativa, entonces obtenemos solución compleja
m=1±i
Solución general
y=c1xcos(ln∣x∣)+c2xsen(ln∣x∣)
Ejemplo 4. ecuación homogénea
Resolver la ecuación de Cauchy Euler x3y′′′−6y=0
Solución
295
Vamos a reemplazar la solución y=xm y sus derivadas en la ED de segundo orden
y=xmy′=mxm−1y′′=m(m−1)xm−2
y′′′=m(m−1)(m−2)xm−3
x3m(m−1)(m−2)xm−3−6xm=0
Agrupamos las x aplicando propiedades de potenciación
m(m−1)(m−2)xm−6xm=0
Factor común
xm(m(m−1)(m−2)−6)=0
Se pasa xm a dividir al lado derecho, se aplica propiedad distributiva y se agrupan términos semejantes
m3−3m2+2m−6=0
Con la división sintética obtenemos una raíz m=3 y con el residuo obtenemos una ecuación de segundo grado.
(m−3)(m2+2)=0
Se despeja el segundo término para hallar las tres raíces de la ecuación auxiliar
296
m=3∧m=±2
Solución general
y=c1x3+c2cos(2ln∣x∣)+c3sen(2ln∣x∣)
Ejemplo 5. ecuación no homogénea
Resolver la ecuación de Cauchy Euler x2y′′−3xy′+3y=2x4ex
Solución
Vamos a reemplazar la solución y=xm y sus derivadas en la ED de segundo orden
y=xmy′=mxm−1y′′=m(m−1)xm−2
x2m(m−1)xm−2−3xmxm−1+3xm=0
Agrupamos las x aplicando propiedades de potenciación
m(m−1)xm−3mxm+3xm=0
Factor común
xm(m(m−1)−3m+3)=0
Se pasa xm a dividir al lado derecho y se aplica distributiva
297
m2−4m+3=0
Factorizando
(m−3)(m−1)=0
Las dos raíces son
m=3∧m=1
La solución complementaria es:
yc=c1x3+c2x
Calculamos el Wronskiano W=[x33x2x1]
W=x3−3x3
W=−2x3
para calcular W1∧W2 debemos llevar la ED a la forma estándar
y′′−x3y′+x23y=2x2ex⇒f(x)=2x2ex
Ahora calculamos W1=[02x2exx1]
W1=0−(2x3ex)=−2x3ex
Ahora calculamos W2=[x33x202x2ex]
W2=2x5ex−0=2x5ex
298
Entonces u1′=−2x3−2x3ex para calcular u1 se integra
∫exdx⇒u1=ex
u2′=−2x32x5ex para calcular u2 se integra
∫−x2e2dx⇒ Se integra por partes
u2=−x2ex+2xex−2ex
Recuerde que yp=u1y1+u2y2
Entonces yp=(ex)(x3)+(−x2ex+2xex−2ex)(x)
yp=x3ex−x3ex+2x2ex−2xex
Agrupando términos semejantes obtenemos
yp=2x2ex−2xex
La solución general es:
299
y=c1x3+c2x+2x2ex−2xex
Ejemplo 6. ecuación no homogénea
Resolver la ecuación de Cauchy Euler x2y′′+8xy′+10y=x−1ln∣x∣
Solución
Vamos a reemplazar la solución y=xm y sus derivadas en la ED de segundo orden
y=xmy′=mxm−1y′′=m(m−1)xm−2
x2m(m−1)xm−2+8xmxm−1+10xm=0
Agrupamos las x aplicando propiedades de potenciación
m(m−1)xm+8mxm+10xm=0
Factor común
xm(m(m−1)+8m+10)=0
Se pasa xm a dividir al lado derecho y se aplica distributiva
m2+7m+10=0
Factorizando
300
(m+5)(m+2)=0
Las dos raíces son
m=−5∧m=−2
La solución complementaria es:
yc=c1x−5+c2x−2
Calculamos el Wronskiano W=[x−5−5x−6x−2−2x−3]
W=−2x−8+5x−8
W=3x−8
para calcular W1∧W2 debemos llevar la ED a la forma estándar
y′′+x8y′+x210y=x2x−1ln∣x∣⇒f(x)=x−3ln∣x∣
Ahora calculamos W1=[0x−3ln∣x∣x−2−2x−3]
W1=0−(x−5ln∣x∣)=−x−5ln∣x∣
Ahora calculamos W2=[x−5−5x−60x−3ln∣x∣]
W2=x−8ln∣x∣−0=x−8ln∣x∣
Entonces u1′=3x−8−x−5ln∣x∣ para calcular u1 se integra
301
∫3−x3ln∣x∣dx
Se integra por partes
u=ln∣x∣dv=−3x3dx
du=x1dxv=−12x4
Aplicando la fórmula de integración por partes
u1=−12x4ln∣x∣+121∫x3dx
u1 es
u1=−12x4ln∣x∣+481x4
Ahora calculamos u2
u2′=3x−8x−8ln∣x∣ para calcular u2 se integra
∫3ln∣x∣dx
Se integra por partes
u=3ln∣x∣dv=dx
du=3x1dxv=x
Aplicando la fórmula de integración por partes
302
u2=3xln∣x∣−31∫dx
u2=3xln∣x∣−3x
Ahora calculamos yp
yp=u1y1+u2y2
yp=(−12x4ln∣x∣+48x4)x−5+(3xln∣x∣−3x)x−2
Aplicando distributiva
yp=−12x−1ln∣x∣+48x−1+3x−1ln∣x∣−3x−1
Agrupando términos semejantes obtenemos
yp=−4x−1ln∣x∣−165x−1
La solución general es:
y=c1x−5+c2x−2−4x−1ln∣x∣−165x−1
Ahora vamos a resolver una ED de Cauchy Euler con valores iniciales
303
Ejemplo 7. Ecuación no homogénea
Resolver la ecuación de Cauchy Euler x2y′′−5xy′+8y=8x6
Sujeta a la condición y(21)=0∧y′(21)=0
Solución
Reemplazamos la solución y=xm y sus derivadas en la ED.
y=xmy′=mxm−1y′′=m(m−1)xm−2
x2m(m−1)xm−2−5xmxm−1+8xm=0
Agrupamos las x aplicando propiedades de potenciación
m(m−1)xm−5mxm+8xm=0
Factor común
xm(m(m−1)−5m+8)=0
Se pasa xm a dividir al lado derecho y se aplica distributiva
m2−6m+8=0
Factorizando
(m−2)(m−4)=0
Las dos raíces son
304
m=2∧m=4
La solución complementaria es:
yc=c1x2+c2x4
Calculamos el Wronskiano W=[x22xx44x3]
W=4x5−2x5
W=2x5
para calcular W1∧W2 debemos llevar la ED a la forma estándar
y′′−x5y′+x28y=8x4⇒f(x)=8x4
Ahora calculamos W1=[08x4x44x3]
W1=0−(8x8)=−8x8
Ahora calculamos W2=[x22x08x4]
W2=8x6−0=8x6
Entonces u1′=2x5−8x8 para calcular u1 se integra
∫−4x3dx
305
u1=−x4
Ahora calculamos u2
u2′=2x58x6 para calcular u2 se integra
∫4xdx
u2=2x2
Ahora calculamos yp
yp=u1y1+u2y2
yp=(−x4)(x2)+(2x2)(x4)
yp=−x6+2x6
Agrupando términos semejantes obtenemos
yp=x6
La solución general es:
y=c1x2+c2x4+x6
Reemplazamos las condiciones iniciales en la solución general
0=41c1+161c2+641
ecuación (1)
Derivamos la solución general para reemplazar la otra condición
306
y′=2c1x+4c2x3+6x5
Reemplazamos C.I
0=c1+21c2+163
ecuación (2). Resolvemos el sistema de ecuaciones
⎩⎨⎧0=41c1+161c2+6410=c1+21c2+163
Multiplicamos la ecuación (2) por −41 y resolvemos el sistema de ecuaciones por eliminación.
⎩⎨⎧−641=41c1+161c2643=−41c1−81c2
321=−161c2
Obtenemos c2=−21
Reemplazando en la ecuación (1) obtenemos c1=161
La solución particular es:
y=161x2−21x4+x6
307
En el siguiente video, realizado por Jaime H. Ramírez Rios, podrás observar un ejemplo de Cauchy Euler homogénea.
En el siguiente video, realizado por Jaime H. Ramírez Rios, podrás observar un ejemplo de Cauchy Euler no homogénea.
308
3.10.1 Ejercicios y respuestas de la sección 3.10
309
310
Capítulo 4
TRANSFORMADA DE LAPLACE
312
4.1 Transformada Integral
En este capítulo se examina un tipo especial de transformada integral llamada transformada de Laplace. Además de tener la propiedad de linealidad, la transformada de Laplace tiene muchas otras propiedades interesantes que la hacen muy útil para resolver problemas lineales con valores iniciales.
La Transformada de Laplace corresponde a un operador T(f) definido por la integral
f(t)=∫0∞k(s,t)f(t)dt(E1)
Con f(t) definido para T≥0
Transforma una función f de la variable t en una función F de la variable s.
La integral (E1) se define como un límite
∫0∞k(s,t)f(t)dt=b→∞lim∫0bk(s,t)f(t)dt
Si el límite existe se dice que la integral converge, en caso contrario diverge. En general el límite existirá sólo para ciertos valores de la variable s.
La función k(s,t) en la ecuación (E1) se llama Kernel o Núcleo de la transformada. La elección de k(s,t)=e−st como el núcleo nos proporciona una integral especialmente importante, la cual vamos a definir a continuación.
313
Definición: sea f una función definida por t≥0 entonces se dice que la integral
L{f(t)}=∫0∞e−stf(t)dt(E2)
Es la transformada de Laplace de f, siempre que la integral converja.
Cuando la integral de la ecuación (E2) converge el resultado es una función de s. En el análisis general se usa una letra minúscula para denotar la función que se transforma y la letra mayúscula correspondiente para denotar su transformada, por ejemplo:
L{f(t)}=F(s)L{g(t)}=G(s)L{y(t)}=Y(s)
Ejemplo 1. Transformada de Laplace
Evaluar la Transformada de Laplace de k
Solución
L{k}=∫0∞ke−stdt
Como no se puede evaluar una integral con un integrando infinito, transformamos la integral en un límite
314
b→∞lim∫0bke−stdt
Sacamos la constante k
b→∞limk∫0be−stdt
Se integra por sustitución
b→∞lim−[se−st]0b
Evaluamos la integral
b→∞lim−se−sb+se0
Evaluamos el límite
b→∞lim−se−sb+se0=s1
Concluimos que la Transformada de Laplace de la función constante k es:
s1
Este resultado es válido para s>0
Ahora vamos a evaluar la Transformada de Laplace de la función identidad t
315
Ejemplo 2. Transformada de Laplace
Evaluar la Transformada de Laplace de t
Solución
L{t}=∫0∞te−stdt
Como no se puede evaluar una integral con un integrando infinito, transformamos la integral en un límite. También se debe tener en cuenta que esta integral es diferente a la anterior porque la t es una variable
b→∞lim∫0tte−stdt
Se integra por partes
u=tdv=e−st
du=dtv=−se−st
Aplicando la fórmula de integración por partes
b→∞lim−ste−st+s1∫0te−stdt
Se integra por sustitución
b→∞lim[−ste−st−s21e−st]0b
Evaluamos la integral
316
b→∞lim[−sbe−sb−s21e−sb]−[−s0e0−s21e0]
b→∞lim[−sbe−sb−s21e−sb]+s21
Evaluamos el límite. Se debe notar que el primer término es una indeterminación, pues multiplicamos ∞ por 0. Por tanto, debemos aplicar el teorema de L´Hopital.
Recuerde que para aplicar el teorema debemos tener un cociente.
b→∞lim−esbb=−s2esb1=0
Al evaluar el límite nos da 0
Concluimos que la Transformada de Laplace de la función identidad t es:
s21
Este resultado es válido para s>0
Con los dos ejemplos anteriores se establece la generalización de la Transformada de Laplace mediante el siguiente teorema llamado transformadas básicas. Se sobreentiende la restricción de s para la convergencia de cada Transformada de Laplace.
317
Tabla 7. Transformada de funciones básicas
Mediante integrales básicas como sustitución y partes vistas en cálculo integral el estudiante puede verificar las transformadas básicas que están en la tabla 7
La transformada de Laplace corresponde a un operador lineal.
L{αf(t)±βg(t)}
Podemos aplicar las mismas propiedades de las integrales porque la transformada de Laplace es una integral, entonces podemos dividir la función anterior en dos transformadas
L{αf(t)}±L{βg(t)}
Ahora sacamos las constantes de las transformadas
αL{f(t)}±βL{g(t)}
318
Ejemplo 3. Transformada básica
Evaluar la Transformada L{5t+2−4sen5t}
Solución
Como la transformada de Laplace es un operador lineal podemos separar la función en tres transformadas.
L{5t}+L{2}−L{4sen5t}
Sacamos las constantes de cada transformada
5L{t}+2L{1}−4L{sen5t}
Ahora aplicamos las transformadas básicas del teorema 1, aplicamos las propiedades 2, 1 y 5
5s21+2s1−4s2+525
Operando los términos obtenemos:
s25+s2−s2+2520
Estos resultados también se pueden obtener integrando cada término
319
Ejemplo 4. Transformada básica
Evaluar la Transformada L{3t4+3cos4t−8e−3t+2senh5t}
Solución
Separamos la función en cuatro transformadas, esto se puede hacer cuando los términos estén separados por sumas o restas
L{3t4}+L{3cos4t}−L{8e−3t}+L{2senh5t}
Ahora sacamos la constante de cada transformadas
3L{t4}+3L{cos4t}−8L{e−3t}+2L{senh5t}
Ahora aplicamos las transformadas básicas del teorema 1, aplicamos las propiedades 3, 4, 6 y 8
3s4+14!+3s2+42s−8s−(−3)1+2s2−525
Operando los términos obtenemos:
s574+s2+163s−s+38+s2−2510
320
Ejemplo 5. Transformada básica
Evaluar la Transformada L{cos2t}
Solución
La transformada de Laplace de cos2t no está en las transformadas básicas, entonces es necesario aplicar identidades para transformar la función
Recuerde la identidad cos2t=21+cos2t
L{cos2t}=L{21+cos2t}=
L{21}+21L{cos2t}
Ahora aplicamos las transformadas básicas del teorema 1, aplicamos las propiedades 1, y 6
Operando los términos obtenemos:
2s1+2(s2+4)s
321
Ejemplo 6. Transformada de una función continua por tramos
Evalúe L{f(t)} donde f(t)={00≤t=32t≥3
Solución
La función que se muestra en la figura, es continua por tramos y de orden exponencial para t>0, la función f se define en dos tramos, la L{f(t)} se expresa como la suma de dos integrales
L{f(t)}=∫030e−stdt+∫3∞2e−stdt
Sólo se evalua la segunda integral, la primera es 0
b→∞lim∫3b2e−stdt
Se saca la constante de la integral
b→∞lim2∫3be−stdt
Integrando
322
b→∞lim[−s2e−st]3b
Evaluamos la integral
b→∞lim[−s2e−sb]−[−s2e−3s]
Evaluando el límite obtenemos
s2e−3s
En el siguiente video, realizado por Jaime H. Ramírez Rios, podrás observar un ejemplo de Transformada de Laplace por tramos.
323
4.1.1 Ejercicios y respuestas de la sección 4.1
324
4.2 Transformada Inversa
El problema inverso:
Si F(s) representa la transformada de Laplace de una función f(s), es decir, L{f(t)}=F(s). Se dice entonces que f(t) es la transformada de Laplace inversa de F(s) y se escribe:
f(t)=L−1{f(s)}
Al igual que la transformada de Laplace básica, existe un teorema para generalizar las transformadas básicas inversas
Tabla 8. Transformada de funciones inversas
La transformada de Laplace inversa corresponde a un operador lineal.
L−1{αf(t)±βg(t)}
325
Podemos aplicar las mismas propiedades de las integrales porque la transformada de Laplace es una integral, entonces podemos dividir la función anterior en dos transformadas
L−1{αf(t)}±L−1{βg(t)}
Ahora sacamos las constantes de las transformadas
αL−1{f(t)}±βL−1{g(t)}
Ejemplo 1. Transformada inversa
Evalúe L−1{s2−163s}
Solución
Como la transformada inversa es un operador lineal, primero se debe sacar la constante de la transformada y descomponer el 16
3L−1{s2−(4)2s}
La función corresponde a la propiedad 8 de la tabla 8, donde k=4
cosh4t
En ocasiones se debe ajustar la función
326
Ejemplo 2. Transformada inversa
Evalúe L−1{s51}
Solución
Corresponde a la propiedad 3 de la tabla 8, pero se debe ajustar.
n+1=5, entonces n=4, por tanto, la función se debe multiplicar y dividir por 4!
4!1L−1{s54!}
241t4
Ejemplo 3. Transformada inversa
Evalúe L−1{s2+71}
Solución
327
Corresponde a la propiedad 5 de la tabla 8, pero se debe ajustar.
k2=7, entonces k=7, por tanto, la función se debe multiplicar y dividir por 7
71L−1{s2+77}
71sen7t
Ejemplo 4. Transformada inversa
Evalúe L−1{s2+5s+4}
Solución
Se debe separar en dos transformadas y sacar la constante
L−1{s2+5s}+4L−1{s2+51}
Corresponde a las propiedades 6 y 5 respectivamente pero se debe ajustar el segundo término´para que coincida con la propiedad 5. Se puede apreciar que k2=5, por tanto, k=5. Se debe ajustar multiplicando y dividiendo por 5
328
L−1{s2+5s}+54L−1{s2+55}
cos5t+545sen5t
Ejemplo 4. Transformada inversa
Evalúe L−1{2s3+3s2−2ss2+2s−1}
Solución
La función no corresponde a ninguna de las propiedades de la tabla 8, se debe separar el denominador
L−1{s(2s−1)(s+2)s2+2s−1}
Se aplica fracciones parciales. Se factoriza el denominador.
s(2s−1)(s+2)s2+2s−1
El denominador tiene tres factores lineales distintos, por lo cual corresponde al caso 1 de fracciones parciales y se organiza de la siguiente manera:
329
El denominador tiene tres factores lineales distintos, por lo cual corresponde al caso 1 de fracciones parciales.
s(2s−1)(s+2)s2+2s−1=sA+2s−1B+s+2C
Esta identidad es válida para todos los valores de s excepto s=0,s=1/2∧s=−2. Para hallar los valores de A, B∧C, multiplicamos ambos lados de esta ecuación por s(2s−1)(s+2)
s2+2s−1=A(2s−1)(s+2)+Bs(s+2)+Cs(2s−1)
Ahora reemplazamos esos valores en la ecuación para hallar A, B∧C
Cuando s=0
−1=A(−1)(2)
A=21
Cuando s=1/2
41=B(21)(25)
B=51
Cuando s=−2
−1=C(−2)(−5)C=−101
330
Con los valores de la constante procedemos a evaluar la transformada inversa
L−1{s(2s−1)(s+2)s2+2s−1}=
21L−1{s1}+51L−1{(2s−1)1}−101L−1{(s+2)1}
Las transformadas corresponden a la propiedad 1 y 4 de la tabla 8, pero el segundo término hay que ajustarlo antes de aplicar las propiedades
La función no corresponde a ninguna de las propiedades de la tabla 8, se debe hacer fracciones parciales para separar el denominador.
s(s2+5)1
El denominador tiene dos factores, uno lineal y uno cuadrático distintos, por lo cual corresponde al caso 3 de fracciones parciales y se organiza de la siguiente manera:
s(s2+5)1=sA+s2+5Bs+C
Esta identidad es válida para todos los valores de s excepto s=0. Para hallar los valores de A, B∧C, multiplicamos ambos lados de esta ecuación por s(s2+5) y obtenemos
1=A(s2+5)+(Bs+C)(s)
Para hallar los valores de las constantes debemos hacer distributiva y armar un sistema de ecuaciones
1=As2+5A+Bs2+Cs
Con esta información vamos a construir un sistema de ecuaciones de tres ecuaciones con tres incógnitas
⎩⎨⎧0=A+B0=C1=5A
Resolviendo el sistema de ecuaciones obtenemos los siguientes valores
332
A=51B=−51c=0
Con los valores de la constante procedemos a evaluar la transformada inversa
L−1{s3+5s1}=51L−1{s1}−51L−1{s2+5s}
Las transformadas corresponden a la propiedad 1 y 6 de la tabla 8
51−51cos5t
En el siguiente video, realizado por Jaime H. Ramírez Rios, podrás observar un ejemplo de Transformada de Laplace inversa.
333
4.2.1 Repasando fracciones parciales
En la escena interactiva que aparece a continuación se explica, el caso 1 de fracciones parcielas. Es importante que repase el concepto antes de resolver ejercicios.
334
4.2.2 Ejercicios y respuestas de la sección 4.2
335
4.3 Transformada de derivadas
TEOREMA 3. Transformada de una derivada: Si f,f′,f′′,…,f(n−1) son continuas en [0,∞) y son de orden exponencial y si f(n)(t) es continua por tramos en [0,∞), entonces
Lfn(t)=snF(s)−s(n−1)f(0)−s(n−2)f′(0)−…−f(n−1)(0)
Donde F(s)=Lf(t)
Solución de problemas con valores iniciales
Para resolver una ecuación diferencial lineal con coeficientes constantes sujeta a condiciones iniciales se utiliza:
L{y′}=sY(s)−y(0)
L{y′′}=s2Y(s)−sy(0)−y′(0)
L{y′′′}=s3Y(s)−s2y(0)−sy′(0)−y′′(0)
la transformada de Laplace de una ecuación diferencial lineal con coeficientes constantes se convierte en una ecuación algebraica en Y(s).
En el ejemplo siguiente se ilustra el método para resolver una ED, así como la descomposición en fracciones parciales.
336
Ejemplo 1. Resolver el PVI
Resolver la Ed y′+2y=t sujeta a la condición inicial y(0)=−1
Solución
Para resolver la ED se reemplaza el teorema tres al lado izquierdo y el teorema uno al lado derecho
sY(s)−y(0)+2Y(s)=s21
Se reemplazan las condiciones iniciales
sY(s)−(−1)+2Y(s)=s21
Se saca factor común Y(s)
Y(s)(s+2)+1=s21
Se despeja el factor común
Y(s)(s+2)=s21−1
Se opera el lado derecho
Y(s)(s+2)=s21−s2
337
Se despeja Y(s)
Y(s)=s2(s+2)1−s2
La solución y(t) del PVI original es y(t)=L−1{Y(s)}, donde la transformada inversa se hace término a término, por tanto, es necesario descomponer la expresión del lado derecho en fracciones parciales para poder aplicar el teorema dos.
s2(s+2)1−s2=sA+s2B+s+2C
Esta identidad es válida para todos los valores de s excepto s=0∧s=−2. Para hallar los valores de A, B∧C, multiplicamos ambos lados de esta ecuación por s2(s+2) y obtenemos
1−s2=As(s+2)+B(s+2)+C(s2)
Para hallar los valores de las constantes debemos hacer distributiva y armar un sistema de ecuaciones
1−s2=As2+2As+Bs+2B+Cs2
Con esta información vamos a construir un sistema de ecuaciones de tres ecuaciones con tres incógnitas
⎩⎨⎧−1=A+C2A+B=01=2B
Resolviendo el sistema de ecuaciones obtenemos los siguientes valores
338
A=−41B=21c=−43
Con los valores de la constante procedemos a evaluar la transformada inversa
Y(s)=−41L−1{s1}+21L−1{s21}−43L−1{s+21}
Aplicando la transformada inversa obtenemos la solución de la ED
y(t)=−41+21t−43e−2t
Ejemplo 2. Resolver el PVI
Resolver la ED y′′+3y′=0 sujeta a la condición inicial y(0)=1∧y′(0)=−1
Solución
Para resolver la ED se reemplaza el teorema tres al lado izquierdo
s2Y(s)−sy(0)−y′(0)+3sY(s)−3y(0)=0
Se reemplazan las condiciones iniciales
339
s2Y(s)−s+1+3sY(s)−3=0
Se saca factor común Y(s)
Y(s)(s2+3s)−s−2=0
Se despeja el factor común
Y(s)(s2+3s)=s+2
Se despeja Y(s) y se saca factor común
Y(s)=s(s+3)s+2
La solución y(t) del PVI original es y(t)=L−1{Y(s)}, donde la transformada inversa se hace término a término, por tanto, es necesario descomponer la expresión del lado derecho en fracciones parciales para poder aplicar el teorema dos.
s(s+3)s+2=sA+s+3B
Esta identidad es válida para todos los valores de s excepto s=0∧s=−3. Para hallar los valores de A∧B, multiplicamos ambos lados de esta ecuación por s(s+3) y obtenemos
s+2=A(s+3)+Bs
Ahora reemplazamos esos valores en la ecuación para hallar A∧B
Cuando s=0
340
2=A(3)
A=32
Cuando s=−3
−1=B(−3)
B=31
Con los valores de la constante procedemos a evaluar la transformada inversa
Y(s)=32L−1{s1}+31L−1{s+31}
Aplicando la transformada inversa obtenemos la solución de la ED
y(t)=32+31e−3t
Ejemplo 3. Resolver el PVI
Resolver la ED y′′−3y′+2y=e−4t sujeta a la condición inicial y(0)=1∧y′(0)=5
341
Solución
Para resolver la ED se reemplaza el teorema tres al lado izquierdo y el teorema uno al lado derecho
s2Y(s)−sy(0)−y′(0)−3sY(s)+3y(0)+2Y(s)=s+41
Se reemplazan las condiciones iniciales
s2Y(s)−s−5−3sY(s)+3+2Y(s)=s+41
Se saca factor común Y(s)
Y(s)(s2−3s+2)−s−2=s+41
Se despeja el factor común
Y(s)(s2−3s+2)=s+41+s+2
Se opera el lado derecho
Y(s)(s2−3s+2)=s+4s2+6s+9
Se despeja Y(s) y se factoriza el trinomio
Y(s)=(s+4)(s−2)(s−1)s2+6s+9
Se descompone el lado derecho en fracciones parciales para poder aplicar el teorema dos.
342
(s+4)(s−2)(s−1)s2+6s+9=s+4A+s−2B+s−1C
Esta identidad es válida para todos los valores de s excepto s=−4,s=2∧s=1. Para hallar los valores de A, B∧C, multiplicamos ambos lados de esta ecuación por (s+4)(s−2)(s−1) y obtenemos
s2+6s+9=A(s−2)(s−1)+B(s+4)(s−1)+C(s+4)(s−2)
Cuando s=−4
1=A(−6)(−5)
A=301
Cuando s=2
25=B(6)(1)
B=625
Cuando s=1
16=C(5)(−1)
C=−516
Con los valores de la constante procedemos a evaluar la transformada inversa
343
Y(s)=301L−1{s+41}+625L−1{s−21}
−516L−1{s−11}
Aplicando la transformada inversa obtenemos la solución de la ED
y(t)=301e−4t+625e2t−516et
Ejemplo 4. Resolver el PVI
Resolver la Ed 2y´´+y′−y=t+1 sujeta a la condición inicial y(0)=1∧y′(0)=0
Solución
Para resolver la ED se reemplaza el teorema tres al lado izquierdo y el teorema uno al lado derecho
2s2Y(s)−2sy(0)−2y′(0)+sY(s)−y(0)−Y(s)=s21+s1
Se reemplazan las condiciones iniciales
2s2Y(s)−2s(1)−2(0)+sY(s)−(1)−Y(s)=s21+s1
344
2s2Y(s)−2s+sY(s)−1−Y(s)=s21+s1
Se saca factor común Y(s)
Y(s)(2s2+s−1)−2s−1=s21+s
Se pasan al lado derecho el término que no dependen de Y(s)
Y(s)(2s2+s−1)=s21+s+2s+1
Se opera el lado derecho
Y(s)(2s2+s−1)=s22s3+s2+s+1
Se despeja Y(s) y se factoriza el trinomio
Y(s)=s2(2s−1)(s+1)2s3+s2+s+1
Se descompone el lado derecho en fracciones parciales
s2(2s−1)(s+1)2s3+s2+s+1=sA+s2B+2s−1C+s+1D
Esta identidad es válida para todos los valores de s excepto s=−21,0,1. Para hallar los valores de A,B,C∧D, multiplicamos ambos lados de esta ecuación por s2(2s−1)(s+1) y obtenemos
Aplicando la transformada inversa obtenemos la solución de la ED
y(t)=−2−t+38e2t+31e−t
En la siguiente escena interactiva, diseñada por Carlos Olvera, podrás resolver ED usando la Transofrmada de Laplace.
347
En el siguiente video, realizado por Jaime H. Ramírez Rios, podrás observar un ejemplo de Transformada de derivada.
348
4.3.1 Ejercicios y respuestas de la sección 4.3
349
4.4 Traslación en el eje "s"
Este teorema facilita encontrar transformadas sin resolver la integral, basta con recorrer la función. Gráficamente se vería así:
Si se considera s una variable real, entonces la gráfica de F(s−a) es la gráfica de F(s) desplazada en el eje s por la cantidad ∣a∣. Si a>0, la gráfica de F(s) se desplaza a unidades a la derecha, mientras que si a<0, la gráfica se desplaza a unidades a la izquierda.
Evaluar transformadas tales como {e−2tt4} o {e3tcosh5t} se puede hacer de forma directa siempre que se conozca {t4} o {cosh5t}. si se conoce la Transformada de Laplace de una función f(t), es posible calcular la Transformada de Laplace de un múltiplo exponencial de f, es decir, {e−atf(t)}, sólo con trasladar o desplazar, F(s) a F(s−a). El resultado se conoce como primer teorema de traslación.
350
Primer teorema de traslación: Si Lf(t)=F(s) y a es cualquier número real, entonces
L{eatf(t)}=F(s−a)
Ejemplo 1. Uso del primer teorema de traslación
Evalúe L{e−2tt4}
Solución
Se utiliza la transformada básica (teorema 1)
El denominador de la transformada de Euler se coloca a la derecha
L{e−2tt4}=L{t4∣s→s+2}
Luego se aplica el teorema 1 a la función y donde va la s se coloca la traslación
(s+2)54!
351
Ejemplo 2. Uso del primer teorema de traslación
Evalúe L{e3tcosh5t}
Solución
Se utiliza la transformada básica (teorema 1)
El denominador de la transformada de Euler se coloca a la derecha
L{e3tcosh5t}=L{cosh5t∣s→s−3}
Luego se aplica el teorema 1 a la función, cosh5t
L{e3tcosh5t}=L{s2−25s∣s→s−3}
Para terminar se debe colocar la traslación que está al lado derecho, en la función, donde esté la s, tanto en el numerador como en el denominador
(s−3)2−25s−3
352
Ejemplo 3. Uso del primer teorema de traslación
Evalúe L{e−4tsen3t}
Solución
Se utiliza la transformada básica (teorema 1)
El denominador de la transformada de Euler se coloca a la derecha
L{e−4tsen3t}=L{sen3t∣s→s+4}
Luego se aplica el teorema 1 a la función sen3t
L{e−4tsen3t}=L{s2+93∣s→s+4}
Donde está la s se coloca la traslación
(s+4)2+93
Forma inversa del teorema:
L−1{F(s−a)}=L−1{F(s)∣s→s−a}
353
Ejemplo 1. Forma inversa del teorema
Evalúe L−1{(s−3)25}
Solución
Se utiliza la transformada inversa (teorema 2)
Se saca la constante de la transformada y la traslación se coloca a un lado
L−1{(s−3)25}=5L−1{s21∣s→s−3}
Se aplica el teorema 2 a la función y a la traslación
5e3tt
Ejemplo 2. Forma inversa del teorema
Evalúe L−1{(s−3)2+162s+5}
354
Solución
Antes de aplicar el teorema dos se debe ajustar el numerador, se busca que la traslación del denominador también esté en el numerador, por ese motivo se suma y se resta 6
L−1{(s−3)2+162s−6+6+5}
Luego se saca factor común y se suma el 6 y el 5
L−1{(s−3)2+162(s−3)+11}
Se separan las dos traslaciones y se sacan las constantes
2L−1{(s−3)2+16(s−3)}+11L−1{(s−3)2+161}
La traslación se coloca a un lado
2L−1{s2+16s∣s→s−3}+11L−1{s2+161∣s→s−3}
Antes de aplicar el teorema 2 se debe ajustar el segundo termino.
2L−1{s2+16s∣s→s−3}+411L−1{s2+164∣s→s−3}
Se aplica el teorema 2 a la función y a la traslación
2e3tcos4t+411e3tsen4t
355
Ejemplo 3. Resolver el PVI
Resolver la Ed y′′−y′=e−2tcos2t sujeta a la condición inicial y(0)=0∧y′(0)=1
Solución
Se utiliza la transformada de la primera derivada (teorema 3) para el lado izquierdo y la transformada básica (teorema 1) para el lado derecho como sigue
s2Y(s)−sy(0)−y′(0)−sY(s)+y(0)=(s+2)2+4s+2
Se reemplazan las condiciones iníciales
s2Y(s)−s−sY(s)+1=(s+2)2+4s+2
factor común Y(s) y se opera el denominador del lado derecho
Y(s)(s2−s)−s+1=s2+4s+8s+2
Se pasan los terminos que no están en el factor comun al lado derecho
Y(s)(s2−s)=s2+4s+8s+2+s−1
356
Se operan los fraccionarios del lado derecho
Y(s)(s2−s)=s2+4s+8s3+3s2+5s−6
Se despeja Y(s)
Y(s)=s(s−1)(s2+4s+8)s3+3s2+5s−6
Se descompone el lado derecho en fracciones parciales
Esta identidad es válida para todos los valores de s excepto s=0,1. Para hallar los valores de A,B,C∧D, multiplicamos ambos lados de esta ecuación por s(s−1)(s2+4s+8) y obtenemos
s3+3s2+5s−6=A(s−1)(s2+4s+8)+Bs(s2+4s+
8)+(Cs+D)s(s−1)
Propiedad distributiva y se construye un sistema de ecuaciones
s3+3s2+5s−6=As3+3As2+4As−8A+Bs3+4Bs2+
8Bs+Cs3−Cs2+Ds2−Ds
A+B+C3A+4B−C+D4A+8B−D−8A====135−6
357
Resolviendo el sistema de ecuaciones obtenemos:
A=43B=133C=521D=132
Con los valores de la constante procedemos a evaluar la transformada inversa
Y(s)=43L−1{s1}+133L−1{s−11}+
521L−1{s2+4s+8s+8}
El tercer término es necesario factorizar el denominador, se ajusta el trinomio cuadrado perfecto y luego se factoriza
Y(s)=43L−1{s1}+133L−1{s−11}+
521L−1{(s+2)2+4s+8}
Ahora se ajusta el numerador
Y(s)=43L−1{s1}+133L−1{s−11}+
521L−1{(s+2)2+4s+2+6}
Se divide el tercer término en dos y se sacan las constantes
358
Y(s)=43L−1{s1}+133L−1{s−11}+
521L−1{(s+2)2+4s+2}+526L−1{(s+2)2+41}
Se coloca la traslación a un lado
Y(s)=43L−1{s1}+133L−1{s−11}+
521L−1{s2+4s∣s→s+2}+263L−1{s2+41∣s→s+2}
Antes de aplicar el teorema 2 se debe ajustar el último término.
Y(s)=43L−1{s1}+133L−1{s−11}+
521L−1{s2+4s∣s→s+2}+(26)23L−1{s2+42∣s→s+2}
Se aplica el teorema 2 a la función y a la traslación
y(t)=43+133et+521e−2tcos2t+523e−2tsen2t
En el siguiente video, realizado por Jaime H. Ramírez Rios, podrás observar un ejemplo de Traslación en el eje s.
359
En el siguiente video, realizado por Jaime H. Ramírez Rios, podrás observar un ejemplo de Traslación en el eje s.
360
4.4.1 Ejercicios y respuestas de la sección 4.4
361
4.5 Traslación en el eje "t"
Definición: Función escalón unitario La función escalón unitario μ(t−a) se define como:
μ(t−a)={0,1,0≤t<at≥a
La función escalón unitario es útil para representar funciones definidas a trozos. Gráficamente se vería así:
Ejemplo 1. Función escalón
suponer f(t)={g(t),h(t),0≤t<at≥a escribir como escalones unitarios
362
Solución
La función escalón está dada por:
f(t)=g(t)+[h(t)−g(t)]μ(t−a)
Aplicando propiedad distributiva queda de la forma:
f(t)=g(t)+h(t)μ(t−a)−g(t)μ(t−a)
Ejemplo 2. Función escalón
suponer f(t)=⎩⎨⎧g(t),h(t),0,0≤t<aa≤t<bt≥b escribir como escalones unitarios
Solución
La función escalón está dada por:
f(t)=g(t)+[h(t)−g(t)]μ(t−a)+[0−h(t)]μ(t−b)
Aplicando propiedad distributiva queda de la forma:
f(t)=g(t)+h(t)μ(t−a)−g(t)μ(t−a)−h(t)μ(t−b)
363
Ejemplo 3. Función escalón
Exprese f(t)=(20t0SiSi0≤t<5t≥5 en términos de la función escalón unitario
Solución
La función escalón está dada por:
f(t)=20t+[0−20t]μ(t−5)
Aplicando propiedad distributiva queda de la forma:
f(t)=20t−20tμ(t−5)
Teorema: Segundo teorema de traslación. Si F(s)=Lf(t) y a>0 entonces
L{f(t−a)μ(t−a)=e−asF(s)
Ahora vamos a evaluar la transformada de Laplace con traslación en el eje t. es importante recordar que se utilizan los tres teoremas vistos desde el comienzo de la unidad
364
Ejemplo 4
evaluar L{(t−1)μ(t−1)}
Solución
La función cumple con el segundo teorema, entonces se aplica la transformada de Laplace básica
s21e−s
Ejemplo 5
evaluar L{(3t+1)μ(t−2)}
Solución
se debe ajustar la función para que cumpla el segundo teorema de traslación
L{(3t−6+6+1)μ(t−2)}
L{(3t−6+7)μ(t−2)}
365
Factor común y propiedad distributiva
L{3(t−2)μ(t−2)+7μ(t−2)}
Ahora se aplica la transformada de Laplace básica
s23e−2s+s7e−2s
Ejemplo 6
evaluar L{costμ(t−π)}
Solución
se debe ajustar la función para que cumpla el segundo teorema de traslación
L{cos(t−π)μ(t−π)}
Como no se puede restar y sumar dentro del ángulo porque vuelve a dar lo mismo, se aplica la identidad de suma y resta de ángulos del coseno
L{(costcosπ+sentsenπ)μ(t−π)}
Evaluando obtenemos
L{−costμ(t−π)}
366
Ahora se aplica la transformada de Laplace básica
−s2+1se−πs
Forma inversa del teorema:Si f(t)=L−1F(s) la forma inversa del teorema anterior con a>0, es
L−1{e−asF(s)}=f(t−a)μ(t−a)
Ejemplo 7
evaluar L−1{s−41e−2s}
Solución
Se resuelve con la tranformada de laplace inversa para {s−41}=e4t
Ahora donde está la t se debe colocar la traslación
e−4(t−2)μ(t−2)
367
Ejemplo 8
evaluar L−1{s2+9se−2πs}
Solución
Tranformada de laplace inversa para {s2+9s}=cos3t
Ahora donde está la t se debe colocar la traslación
L−1{s2+9se−2πs}=cos3(t−2π)μ(t−2π)
Aplicando la identidad para el coseno
cos3(t−2π)μ(t−2π)=cos3tcos32π+sen3tsen32π
Evaluando obtenemos
cos3(t−2π)μ(t−2π)=cos3t(0)+sen3t(−1)
La solución es:
−sen3tμ(t−2π)
368
Ejemplo 9
Evalúe la función f(t)={2−20≤t<3t≥3
Solución
se escribe f(t) como función escalón unitario
f(t)=2+[−2−2]μ(t−3)
f(t)=2−4μ(t−3)
Se aplica la transformada de Laplace básica
s2−s4e−3s
Ejemplo 10
Evalúe la función f(t)=⎩⎨⎧0,sent,0≤t<23πt≥23π
369
Solución
Se escribe f(t) como función escalón unitario
f(t)=0+[sent−0]μ(t−23π)
Propiedad distributiva
f(t)=sentμ(t−23π)
Se ajusta el ángulo
f(t)=sen(t−23π)μ(t−23π)
Aplicando la identidad para el seno
sen(t−23π)μ(t−23π)=sentcos23π−sen23πcost
Evaluando
sen(t−23π)μ(t−23π)=costμ(t−23π)
Se aplica la transformada de Laplace básica
s2+1se−23πs
370
Ejemplo 11
Resolver el PVI y′′−5y′+6y=f(t) sujeta a la condición inicial
y(0)=0∧y′(0)=1 donde f(t)={100≤t<1t≥1
Solución
se escribe f(t) como función escalón unitario
f(t)=1−μ(t−1)
Vamos a resolver el PVI
y′′−5y′+6y=1−μ(t−1)
Se utiliza la transformada de la primera derivada (teorema 3) para el lado izquierdo y la transformada básica (teorema 1) para el lado derecho como sigue
s2Y(s)−sy(0)−y′(0)−5sY(s)+5y(0)+6Y(s)=s1−s1e−s
Se reemplazan las condiciones iníciales
s2Y(s)−s−1−5sY(s)+5+6Y(s)=s1−s1e−s
factor común Y(s) y se despeja
371
Y(s)(s2−5s+6)=s1+s−4−s1e−s
Se operan los términos semejantes al lado derecho
Y(s)(s2−5s+6)=ss2−4s+1−s1e−s
Se despeja Y(s) y se factoriza el termino del lado izquierdo
Y(s)=s(s−2)(s−3)s2−4s+1−s(s−2)(s−3)1e−s
Para aplicar la transformada inversa se debe hacer fracciones parciales a los dos terminos por separado.
Se descompone el primer término en fracciones parciales
s(s−2)(s−3)s2−4s+1=sA+(s−2)B+(s−3)C
Esta identidad es válida para todos los valores de s excepto s=0,2∧3
Para hallar los valores de A,B∧C, multiplicamos ambos lados de esta ecuación por s(s−2)(s−3) y obtenemos
s2−4s+1=A(s−2)(s−3)+Bs(s−3)+Cs(s−2)
Reemplazando las restricciones del dominio en la ecuación y obtenemos los valores de las constantes
A=61B=23C=−32
372
Se descompone el segundo término en fracciones parciales, no vamos a tener en cuenta el signo menos y e−s para hacer la fracción parcial, sólo se vuelven a utilizar en la respuesta
s(s−2)(s−3)1=sA+(s−2)B+(s−3)C
Esta identidad es válida para todos los valores de s excepto s=0,2∧3
Para hallar los valores de A,B∧C, multiplicamos ambos lados de esta ecuación por s(s−2)(s−3) y obtenemos
1=A(s−2)(s−3)+Bs(s−3)+Cs(s−2)
Reemplazando las restricciones en la ecuación obtenemos:
Resolver el PVI y′′+4y=sentμ(t−2π) sujeta a la condición inicial y(0)=1∧y′(0)=0
Solución
Antes de evaluar el escalón se debe ajustar el ángulo
L{sentμ(t−2π)}=L{sen(t−2π)μ(t−2π)}
Se aplica la identidad de suma y resta de ángulos del seno
L{(sentcos2π−costsenπ)μ(t−π)}=L{sentμ(t−2π)}
L{sentμ(t−2π)}=s2+11e−2πs
Se utiliza la transformada básica (teorema 1) para el lado derecho y se iguala al resultado que obtuvimos del escalón
s2Y(s)−sy(0)−y′(0)+4Y(s)=s2+11e−2πs
Se reemplazan las condiciones iníciales
s2Y(s)−s+4Y(s)=s2+11e−2πs
Se saca factor común Y(s) y se despeja
374
Y(s)(s2+4)=s+s2+11e−2πs
Se despeja Y(s)
Y(s)=s2+4s+(s2+4)(s2+1)1e−2πs
Se descompone el último término en fracciones parciales
(s2+4)(s2+1)1=s2+4As+B+s2+1Cs+D
Para hallar los valores de A,B∧C, multiplicamos ambos lados de esta ecuación por (s2+4)(s2+1) y obtenemos
1=(As+B)(s2+1)+(Cs+D)(s2+4)
Aplicamos propiedad sitributiva
1=As3+As+Bs2+B+Cs3+4Cs+Ds2+4D
Construimos un sistema de ecuaciones
A+CB+DA+4CB+4D====0001
Para resolver el sistema de ecuaciones se puede igualar la ecuación 1 con la 3 y la ecuación 2 con la 4 y hacerlo por el método de eliminación.
Se multiplica la ecuación 1 por (−1) y la sumamos a la ecuación 3
375
−A−C=0
A+4C=0
3C=0
Obtenemos el valor de la CC=0
Reemplazando en la ecuación 4 obtenemos A=0
Ahora se multiplica la ecuación 2 por (−1) y la sumamos a la ecuación 4
−B−D=0
B+4D=1
3D=1
Obtenemos el valor de la DD=31
Reemplazando en la ecuación 1 obtenemos B=−31
Con los valores de la constante procedemos a evaluar la transformada inversa
Y(s)=L−1{s2+4s}+
(−31L−1{s2+41}+31L−1{s2+11})e−2πs
376
Antes de aplicar el teorema 2$$ se debe ajustar el término de la mitad.
Y(s)=L−1{s2+4s}+
(−(3)21L−1{s2+42}+31L−1{s2+11})e−2πs
Aplicando la transformada inversa obtenemos:
y(t)=cos2t+[−61sen2(t−2π)+31sen(t−2π)]μ(t−2π)
En el siguiente video, realizado por Jaime H. Ramírez Rios, podrás observar un ejemplo de Traslación en el eje t.
377
4.5.1 Ejercicios y respuestas de la sección 4.5
378
4.6 Derivada de una transformada
Teorema: Derivadas de transformadas
Si F(s)=Lf(t)∧n=1,2,3,… entonces
L{tnf(t)}=(−1)ndsndnF(s)
Ejemplo 1
Evaluar L{tcosh2t)}
Solución
Para aplicar derivada de una transformada debe estar la variable "t" multiplicando cualquier función. Se determina el valor de la n para saber cuántas veces se debe derivar. En este caso n=1
L{tcosh2t)}=(−1)dsds2−4s
L{tcosh2t)}=(−1)⋅(s2−4)2(s2−4)(1)−(s)(2s)
Propiedad distributiva
379
L{tcosh2t)}=(−1)⋅(s2−4)2s2−4−2s2
Agrupando términos semejantes
(s2−4)2s2+4
Ejemplo 2
Evaluar L{t2sen3t)}
Solución
n=2 se debe derivar dos veces
L{t2sen3t)}=(−1)2ds2d2s2+93
Se realiza la primer derivada y se escribe de la siguiente forma
Como hay un producto y aparece Euler también lo podemos resolver como traslación en el eje s
L{e−5tt3}=L{t3∣s→s+5}
L{e−5tt3}=(s+5)46
382
Ejemplo 4
Evaluar L{te2tsen6t}
Solución
n=1 se debe derivar una vez y se hace traslación en el eje s
L{te2tsen6t}=(−1)dsd(s−2)2+366
L{te2tsen6t}=−[(s−2)2+36]2−6[2(s−2)]
La solución es:
[(s−2)2+36]212(s−2)
Ejemplo 5
Evaluar L{tsenkt}
383
Solución
n=1 se debe derivar una vez y se hace traslación en el eje s
L{tsenkt}=(−1)1dsds2+k2k
L{tsenkt}=−(s2+k2)2−k(2s)
La solución es:
(s2+k2)22ks
Ejemplo 6
Resolver el PVI y′′+16y=cos4t sujeta a la condición inicial y(0)=0∧y′(0)=1
Solución
Se utiliza la transformada de la segunda derivada (teorema 3) al lado izquierdo y la transformada básica (teorema 1) al lado derecho
s2Y(s)−sy(0)−y′(0)+16Y(s)=s2+16s
Se reemplazan las condiciones iniciales
384
s2Y(s)−1+16Y(s)=s2+16s
Factor común Y(s) y se despeja
Y(s)(s2+16)=s2+16s+1
Se opera el lado derecho
Y(s)(s2+16)=s2+16s2+s+16
Se despeja Y(s)
Y(s)=(s2+16)2s2+s+16
Se descompone el lado derecho en fracciones parciales
(s2+16)2s2+s+16=s2+16As+B+(s2+16)2Cs+D
Se multiplica ambos lados de la ecuación por (s2+16)2 para quitar denominadores
s2+s+16=(As+B)(s2+16)+Cs+D
Propiedad distributiva
s2+s+16=As3+16As+Bs2+16B+Cs+D
Se arma un sistema de cuatro ecuaciones con cuatro incógnitas para hallar el valor de las constantes
385
AB16A+C16B+D====01116
Los valores de las constantes son:
A=0B=1C=1D=0
Y(s)=L−1{s2+161}+L−1{(s2+16)2s}
Se debe ajustar el primer término y el segundo término antes de aplicar la transformada inversa, para eso se utiliza el resultado del ejemplo 5 para ajustar el segundo término
Y(s)=41L−1{s2+164}+81L−1{(s2+16)28s}
Aplicando la transformada inversa obtenemos la solución
y(t)=41sen4t+81tsen4t
Ejemplo 7
Resolver el PVI y′+y=tsent sujeta a la condición inicial y(0)=0
386
Solución
Al lado derecho hay una traslación en el eje s
tsent=(−1)1dsds2+11
tsent=−(s2+1)2−2s
tsent=(s2+1)22s
Se utiliza la transformada de la segunda derivada (teorema 3) al lado izquierdo y se iguala con el resultado de la traslación
sY(s)−y(0)+Y(s)=(s2+1)22s
Se reemplazan las condiciones iniciales y se saca factor común
Y(s)(s+1)=(s2+1)22s
Se despeja Y(s)
Y(s)=(s+1)(s2+1)22s
Se descompone el lado derecho en fracciones parciales
(s+1)(s2+1)22s=s+1A+s2+1Bs+C+(s2+1)2Ds+E
387
Se multiplica ambos lados de la ecuación por (s+1)(s2+1)2 para quitar denominadores
Aplicando la transformada inversa obtenemos la solución
y(t)=−21e−t+21cost+21tsent−21tcost
389
4.6.1 Ejercicios y respuestas de la sección 4.6
390
4.7 Transformada de integrales
La transformada de Laplace de una solución desconocida de una ED, puede reconocerse algunas veces como el producto de las transformadas de dos funciones conocidas.
Ejemplo 1
resolver el PVI y′′+y=cost sujeta a la condición inicial y(0)=0∧y′(0)=0
Solución
Reemplazando el teorema 1 y el teorema 2
s2Y(s)−sy(0)−y′(0)+Y(s)=s2+1s
Se reemplazan las condiciones iniciales
s2Y(s)+Y(s)=s2+1s
Factor común Y(s)
Y(s)(s2+1)=s2+1s
Se despeja Y(s)
Y(s)=(s2+1)2s
391
Se descompone el lado derecho en fracciones parciales
(s2+1)2s=s2+1As+B+(s2+1)2Cs+D
Se multiplica ambos lados de la ecuación por (s2+1)2 para quitar denominadores
s=(As+B)(s2+1)+Cs+D
Propiedad distributiva
s=As3+As+Bs2+B+Cs+D
Se arma un sistema de cuatro ecuaciones con cuatro incógnitas para hallar el valor de las constantes
ABA+CB+D====0010
Los valores de las constantes son:
A=0B=0C=1D=0E=1
Y(s)=L−1{(s2+1)2s}
No hay una transformada que nos ayude a resolver el ejercicio. Debe haber una forma de combinar las 2 funciones sent∧cost para obtener y(t) cuya transformada sea el producto de sus transformadas.
392
Pero y(t) no es el producto de cost por sent debido a que:
L{costsent}L{21sen2t}
Evaluando la transformada
s2+41=(s2+1)2s
Se puede concluir que:
L{costsent}=L{cost}L{sent}
El teorema que nos ayuda a resolver el problema indica que la función:
h(t)=∫0tf(τ)g(t−τ)dτ
Tiene la propiedad
L{h(t)}=H(s)=F(s)⋅G(s)
La función de t definida como la integral
∫0tf(τ)g(t−τ)dτ
Depende sólo de f∧g, y se conoce como la convolución de f∧g.
Se representa como f∗g, de tal forma que su transformada es el producto de las transformadas de f∧g.
393
Convolución:
si las funciones f∧g son continuas por tramos en el intervalo [0,∞) entonces un producto especial, denotado por f∗g se define mediante la integral
f∗g=∫0tf(τ)g(t−τ)dτ
Y se llama convolución de f∧g. La convención de f∗g es una función de t
Teorema de convolución
si f(t)∧g(t) son funciones continuas por tramos en el intervalo [0,∞) y de orden exponencial, entonces
L{f∗g}=L{f(t)}L{g(t)}=F(s)G(s)
Para terminar el ejemplo 1 debemos tener en cuenta que:
La convolución de sent∧cost es
sent∗cost=∫0tsenτcos(t−τ)dτ
Aplicando la identidad de producto a suma
senmx.cosnx=21(sen[(m−n)x]+sen[(m+n)x])
394
∫0tsenτcos(t−τ)dτ=21∫0tsen(τ−t+τ)+sen(τ+t−τ)dτ
Agrupando términos semejantes
21∫0tsen(2τ−t)+sen(t)dτ
Integrando obtenemos
21[−21cos(2τ−t)+sen(t)τ]0t
Reemplazando los límites de integración
21[(−21cos(t)+sen(t)t)−(−21cos(−t)+sen(t)0)]
Ahora evaluamos
21[−21cos(t)+tsent+21cos(t)]
Agrupando términos semejantes obtenemos la solución
21tsent
La solución del ejemplo 1 es:
y(t)=21tsent
395
Ejemplo 2
Evaluar L{2∗t2}
Solución
L{2∗t2}=s2⋅s32!
s44
Ejemplo 3
Evaluar L{e−2t∗senh3t}
Solución
L{e−2t∗senh3t}=s+21⋅s2−93
(s+2)(s2−9)3
396
Ejemplo 4
Evaluar L{∫0tcosτdτ}
Solución
La función f(τ)=1 y la función g(t−τ)=cosτ por tanto, se aplica el teorema de convolución
L{∫0tcosτdτ}=s1⋅s2+1s
Al multiplicar los dos resultados obtenemos la solución:
s2+11
Ejemplo 5
Evaluar L{t∫0tsenτdτ}
Solución
La función f(τ)=1 y la función g(t−τ)=senτ, se aplica el teorema de convolución, pero afuera de la integral hay una t, por
397
tanto, también se debe aplicar la derivada de una transformada
L{t∫0tsenτdτ}=s1⋅(−1)dsds2+11
L{t∫0tsenτdτ}=s1⋅(−1)⋅(s2+1)2−2s
Al multiplicar signos y simplificar obtenemos:
(s2+1)22
Ejemplo 6
Hallar la convolución de f(t)=sen2t∧g(t)=et
Solución
Para resolver el ejercicio debemos aplicar el teorema de convolución integrando las dos funciones
sen2t∗et=∫0tet−τsen2τdτ
Aplicando propiedades de potenciación
∫0tet⋅e−τsen2τdτ=et∫0te−τsen2τdτ
398
Recuerde la integral
cos2t=∫eaxsenbxdx=a2+b2eax(asenbx−bcosbx)+c
Reemplazamos los valores teniendo en cuenta que: a=−1∧b=2, obtenemos:
et∫0te−τsen2τdτ=et(5e−τ(−sen2τ−2cos2τ)]0t
Ahora debemos evaluar la integral reemplazando los límites de integración
et[(5e−t(−sen2t−2cos2t))−(5e0(−sen0−2cos0))]
Evaluando obtenemos
et[(5e−t(−sen2t−2cos2t))−(51(−2))]
Operando el último término
et[(5e−t(−sen2t−2cos2t))+52]
Aplicando la propiedad distributiva obtenemos la solución
52et−52sen2t−52cos2t
399
4.7.1 Ejercicios y respuestas de la sección 4.7
400
4.7.2 Ejercicios y respuestas del capítulo 4
401
402
Capítulo 5
SERIES DE FOURIER
404
5.1 Funciones ortogonales
Las series e integrales de Fourier establecen un tema clásico del análisis matemático. Aparecen en el siglo XVIII como resultado del estudio de las vibraciones de una cuerda, las series de Fourier han contribuido al desarrollo de los conceptos básicos del análisis –función, integral, serie, convergencia y se han obtenido por los trabajos de varios matemáticos sobre series trigonométricas.
Las series de Fourier son series de términos coseno y seno y surgen en la tarea práctica de representar funciones periódicas generales. Como aplicación constituyen una herramienta muy importante en la solución de problemas en los que intervienen ecuaciones diferenciales ordinarias y parciales.
Otra propiedad esencial de las funciones es generalizar el concepto de ortogonalidad, el concepto de producto interno de vectores pierde su interpretación geométrica en las funciones. Pero dicho concepto se ha generalizado y es común considerar una función como un vector. Entonces se puede decir que dos funciones distintas son ortogonales cuando su producto interno es cero, este producto interno es en realidad una integral definida. Este concepto se aplica en movimientos amortiguados de una masa “m” en un resorte y en oscilaciones amortiguadas forzadas.
Producto interno de funciones:
El producto interno de dos funciones f1∧f2 en un intervalo [a,b] es el número.
(f1,f2)=∫abf1(x)f2(x)dx
405
Dos funciones f1∧f2 son ortogonales en un intervalo [a,b] sí.
(f1,f2)=∫abf1(x)f2(x)dx=0
Un conjunto de funciones de valor real {φ0(x),φ1(x),φ2(x)…} se dice que es ortogonal en un intervalo [a,b] si
(φm,φn)=∫abφm(x)φn(x)dx=0m=n
Ejemplo 1
demostrar que las funciones f1(x)=x2∧f2(x)=x3 son ortogonales en el intervalo [−1,1]
Solución
Se debe demostrar que el producto interno entre las dos funciones es cero
406
(f1,f2)=∫−11x2x3dx
Obtenemos la integral
∫−11x5dx=6x6∣∣−11
Evaluando los límites
61−61=0
Como el producto interno de las funciones es cero, queda demostrado que las funciones son ortogonales
Ejemplo 2
demostrar que las funciones f1(x)=cosx∧f2(x)=sen2x son ortogonales en el intervalo [0,π]
Solución
(f1,f2)=∫0πcosxsen2xdx
Se integra por sustitución
u=senxdu=cosxdx
407
∫0πcosxsen2xdx=∫0πu2du=3u3=3sen3x∣∣0π
Evaluando obtenemos
3sen3π−3sen30=0−0=0
Las funciones son ortogonales
Conjuntos ortonormales: la norma o longitud ∥u∥ de un vector u, se puede expresar en términos del producto interno. La expresión (u,u)=∥∥u2∥∥ se llama norma cuadrada, por lo que la norma es ∥u∥=(u,u). De igual modo la norma cuadrada de una función φ0 es ∥φ(x)∥2=(φn,φn) y así la norma o su longitud generalizada es ∥φ(x)∥=(φn,φn). En otras palabras, la norma cuadrada y la norma de una función φn en un conjunto ortogonal {φn(x)} son, respectivamente
∥φ(x)∥2=∫abφn2(x)dx
∥φ(x)∥=∫abφn2(x)dx
Ejemplo 3
Demuestre que el conjunto {1,cosx,cos2x,...} es ortogonal en el intervalo [−π,π] y encuentre las normas de cada función.
408
Solución
Se debe demostrar que el producto interno de las funciones {φ0(x),φ1(x),φ2(x)...} es cero, donde:
φ0(x)=1,φ1(x)=cosx,φ2(x)=cos2x,φn(x)=cosnx...
Primero se demuestra que φ0(x),φn(x)=0, es decir, el 1 con todos los cosenos