Archivo de la etiqueta: Wronskiano

Ecuaciones Diferenciales I: Soluciones a sistemas de ecuaciones diferenciales

Por Omar González Franco

Los errores y dificultades no resueltos en el pasado de las matemáticas
siempre han sido las oportunidades de su futuro.
– E. T. Bell

Introducción

En la entrada anterior vimos lo que es un sistema de ecuaciones diferenciales, en particular un sistema lineal de primer orden. Vimos también lo que es un problema de valores iniciales y establecimos la notación matricial.

Así mismo, vimos cómo es que una ecuación diferencial lineal de orden n se puede transformar en un sistema lineal de primer orden, esto tiene bastante ventaja ya que, una vez que veamos cómo resolver sistemas de ecuaciones diferenciales, muchas veces será más sencillo resolver el sistema que resolver la ecuación de orden n aplicando los métodos que ya conocemos.

En esta entrada estudiaremos las propiedades de las soluciones de los sistemas lineales de primer orden.

Cabe mencionar que mucho de lo que desarrollaremos en esta entrada es bastante similar a la teoría vista con las ecuaciones diferenciales de orden n, comenzando por la validez del principio de superposición.

A partir de ahora sólo usaremos la notación matricial y toda la teoría básica del álgebra lineal que éstas conllevan.

Soluciones de sistemas lineales de primer orden

Comencemos por estudiar el caso homogéneo. El sistema lineal de primer orden homogéneo es

(1)(y1y2yn)=(a11(t)a12(t)a1n(t)a21(t)a22(t)a2n(t)an1(t)an2(t)ann(t))(y1y2yn)

O bien,

(2)Y=AY

En la entrada anterior definimos la solución de un sistema de ecuaciones diferenciales en el intervalo δ como el conjunto de n funciones

(3)S0={y1(t),y2(t),,yn(t)}

definidas en δ y diferenciables en el mismo intervalo, tales que satisfacen simultáneamente las n ecuaciones diferenciables de un sistema lineal.

Las soluciones pueden ser escritas como el vector

(4)Y=(y1(t)y2(t)yn(t))

cuyos elementos son funciones derivables que satisfacen un sistema lineal en el intervalo δ.

En las siguientes definiciones y teoremas se supondrá que los coeficientes aij(t), i,j{1,2,3,,n} y ,para el caso no homogéneo, las funciones gi(t), son continuas en algún intervalo común δ.

Comencemos por mostrar que el principio de superposición también es valido para sistemas lineales.

Demostración: Consideremos la combinación lineal

Y=c1Y1+c2Y2++cmYm

con

Yi=(y1iy2iyni)

para i=1,2,,m. La derivada de Yi esta dada por

Yi=(y1iy2iyni)

Entonces la derivada de la combinación lineal es

Y=(c1y11+c2y12++cmy1mc1y21+c2y22++cmy2mc1yn1+c2yn2++cmynm)=c1(y11y21yn1)+c2(y12y22yn2)++cm(y1my2mynm)=c1Y1+c2Y2++cmYm

Como cada Yi, i=1,2,,m, es solución del sistema homogéneo (2) en δ, entonces

Yi=AYi

así

Y=c1(AY1)+c2(AY2)++cm(AYm)=A(c1Y1+c2Y2++cmYm)=AY

En donde se ha hecho uso de la propiedad distributiva de la matriz A y de la hipótesis (5). Por lo tanto, la combinación lineal

Y=c1Y1+c2Y2++cmYm

también es solución y los es en el mismo intervalo común δ ya que esta compuesta de soluciones definidas en dicho intervalo.

◻

Intenta hacer la demostración.

Realicemos un ejemplo.

Ejemplo: Probar que la combinación lineal

Y=c1Y1+c2Y2+c3Y3=c1(110)+c2(e2te2t0)+c3(00e3t)

es solución del sistema lineal

Y=(110110003)Y

Solución: Probemos que cada uno de los vectores de la combinación lineal es solución y usemos el principio de superposición.

Los vectores son

Y1=(110),Y2=(e2te2t0),Y3=(00e3t)

Por un lado, derivemos estos vectores.

Y1=(000),Y2=(2e2t2e2t0),Y3=(003e3t)

Por otro lado, sustituyamos cada uno de los vectores en el sistema lineal y usemos los resultados anteriores.

AY1=(110110003)(110)=(11110)=(000)=Y1

AY2=(110110003)(e2te2t0)=(e2t+e2te2t+e2t0)=(2e2t2e2t0)=Y2

y

AY3=(110110003)(00e3t)=(003e3t)=Y3

De esta manera queda mostrado que los tres vectores son solución, ya que satisfacen el sistema. Por el principio de superposición concluimos que la combinación lineal

Y=c1Y1+c2Y2+c3Y3=c1(110)+c2(e2te2t0)+c3(00e3t)

también es solución del sistema lineal.

◻

El principio de superposición nos indica que un sistema lineal puede tener más de una solución, sin embargo, similar al caso de ecuaciones diferenciales de orden n, buscamos soluciones que sean linealmente independientes entre sí. A continuación definimos la dependencia e independencia lineal de las soluciones en este contexto.

En la unidad anterior definimos una herramienta muy útil que, además de ayudarnos a resolver ecuaciones diferenciales de orden superior en algunos métodos, nos ayuda a determinar si un conjunto de soluciones es linealmente independiente, dicha herramienta es el Wronskiano, la definición en el caso de los sistemas lineales de primer orden, es la siguiente.

Se puede demostrar que si el Wronskiano es distinto de cero, entonces las soluciones son linealmente independientes, igual que antes, esto es conocido como el criterio para soluciones linealmente independientes. Para demostrar este hecho es conveniente recordar algunos resultados de álgebra que podremos usar en la demostración.

Recordemos que un sistema lineal de n ecuaciones con n incógnitas es un conjunto de ecuaciones

(9)b11u1+b12u2++b1nun=d1b21u1+b22u2++b2nun=d2bn1u1+bn2u2++bnnun=dn

Con bi,j y di, i,j{1,2,3,,n} números reales dados y ui, i=1,2,,n las incógnitas. Usando la notación matricial podemos escribir el sistema (9) como

(10)BU=D

con

B=(b11b12b1nb21b22b2nbn1bn2bnn),U=(u1u2un),D=(d1d2dn)

Los resultados que nos interesan son los siguientes.

Si D=0, el sistema (10) también recibe el nombre de sistema homogéneo.

Con estos resultados podemos demostrar el criterio para soluciones linealmente independientes que se enuncia a continuación.

Demostración:

) Por demostrar: W(Y1,Y2,,Yn)0.

Sea t0δ en el que W(t0)=0, en donde W(t0) denota al Wronskiano con cada vector solución evaluado en el punto t0.

W(t0)=W(Y1(t0),Y2(t0),,Yn(t0))

En una combinación de ambos teoremas de los resultados de álgebra podemos deducir que existen constantes c1,c2,,cn, no todos cero, tal que

(11)Y(t0)=c1Y1(t0)+c2Y2(t0)++cnYn(t0)=0

Lo que tenemos es un sistema lineal de n ecuaciones homogéneo con n incógnitas (sistema lineal en el contexto algebraico (10) con D=0, no sistema lineal de ecuaciones diferenciales), dichas incógnitas son las constantes ci, i=1,2,,n. La relación (11) se cumple debido a que si el Wronskiano es igual a cero, entonces es posible que el sistema no tenga solución trivial y mucho menos una solución única, esto lo deducimos de los teoremas de álgebra que establecimos.

Por otro lado, sabemos por hipótesis que los vectores Y1,Y2,,Yn son solución del sistema homogéneo (2) en el intervalo δ, por el principio de superposición sabemos también que la combinación lineal

Y(t)=c1Y1+c2Y2++cnYn

es solución de (2) en δ. Del resultado (11) y de la unicidad de la solución se deduce que Y(t)=0 para algún punto t=t0δ, es decir,

c1Y1+c2Y2++cnYn=0

Pero por hipótesis los vectores Y1,Y2,,Yn son linealmente independientes en δ, lo que implica que

c1=c2==cn=0

lo cual es una contradicción con lo que establecimos en (11). Por lo tanto, el Wronskiano tiene que ser distinto de cero, es decir

W(Y1,Y2,,Yn)0

) Por demostrar: S es linealmente independiente.

Este caso también lo demostraremos por contradicción. Supongamos que los vectores solución Y1,Y2,,Yn son linealmente dependientes en δ, esto implica que existen constantes c1,c2,,cn no todos cero, tal que

c1Y1+c2Y2++cnYn=0

Este sistema lo podemos escribir en la forma (9) como

c1y11+c2y12++cny1n=0c1y21+c2y22++cny2n=0c1yn1+c2yn2++cnynn=0

En donde las funciones yij, i,j{1,2,3,,n} son los coeficientes y las constantes ci, i=1,2,,n son las incógnitas. Debido a que las ci no son todas cero implica que el sistema no tiene solución trivial y por el segundo teorema de los resultados de álgebra concluimos que

W(Y1,Y2,,Yn)=0

Pero, por hipótesis

W(Y1,Y2,,Yn)0

lo cual es una contradicción y todo nace de considerar a S como un conjunto linealmente dependiente. Por lo tanto, el conjunto de soluciones

S={Y1,Y2,,Yn}

es linealmente independiente en δ.

◻

Un resultado interesante se enuncia a continuación.

Este resultado nos garantiza que si W0 para algún punto t0δ, entonces W0 para toda tδ y por el criterio anterior las soluciones serán linealmente independientes en ese intervalo.

El conjunto de soluciones linealmente independientes del sistema lineal (2) recibe un nombre especial.

El siguiente teorema nos garantiza la existencia de este conjunto.

El conjunto fundamental de soluciones está constituido por vectores que son linealmente independientes entre sí, con estos vectores es posible formar una matriz cuyas columnas están formadas con las entradas de dichos vectores, esta matriz tiene un nombre especial.

Un hecho interesante es que el determinante de la matriz fundamental de soluciones corresponde al Wronskiano.

(13)W(Y1,Y2,,Yn)=|M(t)|

Realicemos un ejemplo, para ello consideremos el sistema lineal del ejemplo anterior.

Ejemplo: Mostrar que las soluciones

Y1=(110),Y2=(e2te2t0),Y3=(00e3t)

del sistema lineal

Y=(110110003)Y

son linealmente independientes.

Solución: En el ejemplo anterior ya comprobamos que efectivamente son solución del sistema lineal dado. Para determinar si son linealmente independientes veamos si el Wronskiano es distinto de cero.

W(Y1,Y2,Y3)=|1e2t01e2t000e3t|=e5t+0+000(e5t)=2e5t0

Como W0, tR, entonces los vectores dados son linealmente independientes y por lo tanto forman un conjunto fundamental de soluciones en R.

S={(110),(e2te2t0),(00e3t)}

La matriz fundamental de soluciones es

M(t)=(1e2t01e2t000e3t)

◻

Un buen ejercicio sería mostrar que un conjunto de soluciones del sistema lineal homogéneo (2) forma un espacio vectorial, es relativamente sencillo probar cada una de las propiedades o axiomas que definen a un espacio vectorial. El resultado a demostrar de tarea moral es el siguiente.

Soluciones generales a sistemas lineales

Ahora que conocemos algunas propiedades de las soluciones de sistemas lineales, es momento de conocer la forma general de las soluciones de los sistemas lineales tanto homogéneos como no homogéneos.

Comencemos por enunciar el teorema que establece la forma de la solución general de un sistema lineal homogéneo (2).

Demostración: Sea Y(t) una solución arbitraria del sistema lineal homogéneo en el intervalo δ, sea t0δ y supongamos que

Y(t0)=(b1b2bn)=Y0

Es decir, la función Y(t) satisface el problema de valores iniciales Y=AY;Y(t0)=Y0.

Por otro lado, por el principio de superposición sabemos que la combinación lineal

Y^(t)=c1Y1(t)+c2Y2(t)++cnYn(t)

también es solución del sistema lineal Y=AY. Donde ci, i=1,2,,n son constantes arbitrarias y las Yi, i=1,2,,n son las soluciones del conjunto fundamental de soluciones del sistema lineal. Supongamos que

Y^(t0)=c1Y1(t0)+c2Y2(t0)++cnYn(t0)=Y0

Lo que tenemos es el siguiente sistema de n ecuaciones.

c1y11(t0)+c2y12(t0)++cny1n(t0)=b1c1y21(t0)+c2y22(t0)++cny2n(t0)=b2c1yn1(t0)+c2yn2(t0)++cnynn(t0)=bn

En donde las incógnitas son las contantes ci, i=1,2,,n. Como las funciones yij, i,j{1,2,3,,n} pertenecen a vectores del conjunto de soluciones, entonces sabemos que Y1,Y2,,Yn son linealmente independientes y por el criterio para soluciones linealmente independientes inferimos que W(t0)0, donde

W(t0)=W(Y1(t0),Y2(t0),,Yn(t0))

De los resultados de álgebra deducimos que el sistema de n ecuaciones tiene solución única, esto significa que existen constantes únicas c1,c2,,cn, tal que

c1Y1(t0)+c2Y2(t0)++cnYn(t0)=Y0

Esto nos indica que

Y^(t)=c1Y1(t)+c2Y2(t)++cnYn(t)

es solución del problema de valores iniciales. Por el teorema de existencia y unicidad para sistemas lineales homogéneas concluimos que Y(t)=Y^(t), es decir,

Y(t)=c1Y1(t)+c2Y2(t)++cnYn(t)

Como Y(t) es una solución arbitraria, entonces debe ser la solución general del sistema lineal homogéneo en δ.

◻

Para concluir la entrada estudiemos el caso no homogéneo.

Sistemas no homogéneos

El sistema lineal de primer orden no homogéneo es

(15)(y1(t)y2(t)yn(t))=(a11(t)a12(t)a1n(t)a21(t)a22(t)a2n(t)an1(t)an2(t)ann(t))(y1(t)y2(t)yn(t))+(g1(t)g2(t)gn(t))

O bien,

(16)Y=AY+G

El vector de funciones que satisface el sistema (16) es una solución y recibe un nombre.

A continuación se enuncia el teorema que nos muestra la forma general de la solución de un sistema lineal no homogéneo.

Demostración: Sea

Yp(t)=(y1py2pynp)

una solución particular de (16) y sean Y1(t),Y2(t),,Yn(t), n soluciones linealmente independientes del sistema homogéneo asociado Y=AY.

Sea Y(t) una solución arbitraria del sistema no homogéneo, notemos lo siguiente.

(Y(t)Yp(t))=Y(t)Yp(t)=(AY(t)+G)(AYp(t)+G)=A(Y(t)Yp(t))

Este resultado nos indica que Y(t)Yp(t) es solución del sistema homogéneo, eso significa que se puede escribir como

Y(t)Yp(t)=c1Y1(t)+c2Y2(t)++cnYn(t)

entonces, la solución Y tiene la forma

(19)Y(t)=c1Y1(t)+c2Y2(t)++cnYn(t)+Yp(t)

La solución Y(t), al ser cualquier solución del sistema lineal no homogéneo, podemos deducir que la solución general debe tener la forma (19), por lo que concluimos que Y(t) se trata de la solución general de (16).

Considerando la hipótesis (17) concluimos que la solución general del sistema lineal no homogéneo es

Y(t)=Yc(t)+Yp(t)

◻

Cuando estamos trabajando con un sistema lineal no homogéneo, la solución general del sistema lineal homogéneo asociado (17) recibe un nombre particular.

Concluyamos con un ejemplo.

Ejemplo: Probar que el vector

Yp=(14e2t+12te2tet+14e2t+12te2t12t2e3t)

es una solución particular del siguiente sistema lineal no homogéneo.

Y=(110110003)Y+(ete2tte3t)

Solución: Por un lado, derivemos el vector dado.

Yp=(12e2t+12e2t+te2tet+12e2t+12e2t+te2tte3t+32t2e3t)=(te2tet+e2t+te2tte3t+32t2e3t)

Por otro lado, sustituyamos directamente en el sistema al vector dado.

Y=(110110003)(14e2t+12te2tet+14e2t+12te2t12t2e3t)+(ete2tte3t)

Operando obtenemos lo siguiente.

Y=(14e2t+12te2tet+14e2t+12te2t+et14e2t+12te2tet+14e2t+12te2t+e2t32t2e3t+te3t)=(te2tet+e2t+te2tte3t+32t2e3t)

Los resultados obtenidos son los mismos, por lo tanto el vector Yp es solución del sistema.

En los ejemplos anteriores de esta entrada probamos que el conjunto fundamental de soluciones del sistema lineal homogéneo asociado

Y=(110110003)Y

esta constituido por los vectores linealmente independientes

Y1=(110),Y2=(e2te2t0),Y3=(00e3t)

de manera que la función complementaria es

Yc=c1(110)+c2(e2te2t0)+c3(00e3t)

Como la solución general es

Y=Yc+Yp

Entonces la solución general del sistema lineal no homogéneo es

Y(t)=c1(110)+c2(e2te2t0)+c3(00e3t)+(14e2t+12te2tet+14e2t+12te2t12t2e3t)

◻

Tarea moral

Los siguientes ejercicios no forman parte de la evaluación del curso, pero servirán para entender mucho mejor los conceptos vistos en esta entrada, así como temas posteriores.

  1. Los siguientes vectores son soluciones de un sistema lineal homogéneo Y=AY. Determinar si forman un conjunto fundamental de soluciones en R.
  • Y1=(11)et,Y2=(26)et+(88)tet
  • Y1=(1613),Y2=(121)e4t,Y3=(232)e3t
  1. Probar que el vector Yp es una solución particular del sistema lineal dado.
  • Y=(2134)Y(17)et,Yp=(11)et+(11)tet
  • Y=(123420610)Y+(143)sin(3t),Yp=(sin(3t)0cos(3t))
  1. Mostrar que la solución general de

    Y=(060101110)Y

    en el intervalo (,) es

    Y=c1(615)et+c2(311)e2t+c3(211)e3t
  1. Mostrar que la solución general de

    Y=(1111)Y+(11)t2+(46)t+(15)

    en el intervalo (,) es

    Y=c1(112)e2t+c2(11+2)e2t+(10)t2+(24)t+(10)
  1. Demostrar que el conjunto de soluciones del sistema lineal homogéneo Y=AY forma un espacio vectorial con la suma y el producto por escalares usuales de matrices.

Más adelante…

Ahora que conocemos lo que son los sistemas lineales de ecuaciones diferenciales y las propiedades de sus soluciones estamos casi listos para comenzar a desarrollar los distintos métodos de resolución, sin embargo, antes de ello es necesario definir una herramienta matemática que será de suma utilidad en el desarrollo posterior de esta unidad. Dicha herramienta es la exponencial de una matriz.

En la siguiente entrada definiremos lo que significa eAt, donde A es una matriz de n×n con componentes constantes y veremos como se relaciona con un sistema lineal Y=AY. Así mismo, profundizaremos en el concepto de matriz fundamental de soluciones.

Entradas relacionadas

Agradecimientos

Trabajo realizado con el apoyo del Programa UNAM-DGAPA-PAPIME PE104522 «Hacia una modalidad a distancia de la Licenciatura en Matemáticas de la FC-UNAM – Etapa 2»

Ecuaciones Diferenciales I: Ecuaciones lineales no homogéneas de segundo orden – Método de variación de parámetros

Por Omar González Franco

Las matemáticas son un lugar donde puedes hacer
cosas que no puedes hacer en el mundo real.
– Marcus du Sautoy

Introducción

Con lo que hemos estudiado en las dos últimas entradas somos capaces de resolver ecuaciones diferenciales lineales de segundo orden homogéneas y no homogéneas con coeficientes constantes, es decir, ecuaciones de la forma

Label '1' multiply defined

Con a,b,c constantes y g(x)=0 en el caso homogéneo o g(x)0 en el caso no homogéneo, en éste último caso aún estamos limitados a la forma que puede tener la función g, pues sabemos resolver las ecuaciones diferenciales en el caso en el que la función g es una constante, una función polinomial, una función exponencial, funciones seno o coseno, o una combinación entre ellas. La pregunta ahora es, ¿cómo resolver este tipo de ecuaciones para cualquier tipo de función g(x)?.

En esta entrada desarrollaremos un método que nos permite obtener la solución general independientemente de la forma que tenga la función g(x). A dicho método se le conoce como variación de parámetros.

El nombre de este método resulta familiar. En la unidad anterior desarrollamos éste método para el caso de las ecuaciones diferenciales lineales de primer orden como método alterno al método por factor integrante. Lo que haremos en esta entrada es una adaptación del método de variación de parámetros para el caso en el que las ecuaciones diferenciales son de orden superior, en particular, de segundo orden.

Variación de parámetros

Consideremos la ecuación diferencial

Label '2' multiply defined

Si a2(x)0 para toda x en el intervalo δ en el que está definida la solución, entonces podemos definir las funciones

P(x)=a1(x)a2(x),Q(x)=a0(x)a2(x)yg(x)=g0(x)a2(x)

de manera que la ecuación (2) la podemos escribir en su forma estándar como

Label '3' multiply defined

En el caso de primer orden se hizo la suposición de que la solución particular era de la forma

yp(x)=k(x)y1(x)=k(x)eP(x)dx

Manteniendo esta idea, en el caso de segundo orden se busca una solución de la forma

Label '4' multiply defined

Donde y1 y y2 forman un conjunto fundamental de soluciones en δ de la ecuación homogénea asociada a (3). Determinemos la primera y segunda derivada de yp para sustituir los resultados en la ecuación diferencial (3).

Label '5' multiply defined

y

Label '6' multiply defined

Sustituyendo en (3) y reorganizando los términos obtenemos lo siguiente.

k1[d2y1dx2+Pdy1dx+Qy1]+y1d2k1dx2+dk1dxdy1dx+k2[d2y2dx2+Pdy2dx+Qy2]+y2d2k2dx2+dk2dxdy2dx+P[y1dk1dx+y2dk2dx]+dy1dxdk1dx+dy2dxdk2dx=g(x)

Como y1 y y2 son soluciones de la ecuación homogénea asociada, entonces

d2y1dx2+Pdy1dx+Qy1=0yd2y2dx2+Pdy2dx+Qy2=0

y además notamos que

Label '7' multiply defined

y

Label '8' multiply defined

Considerando lo anterior la ecuación diferencial queda como

Label '9' multiply defined

Nuestro propósito es determinar a las funciones k1(x) y k2(x) de (4), esto implica que debemos formar un sistema con dos ecuaciones que debemos resolver para obtener dichas funciones. De acuerdo al resultado obtenido vamos a establecer la restricción de que las funciones k1 y k2 satisfacen la relación

Label '10' multiply defined

Considerando esto la ecuación se reduce a

Label '11' multiply defined

Las ecuaciones (10) y (11) corresponden al sistema de dos ecuaciones que debemos resolver.

Como podemos notar, es un sistema para determinar las derivadas de las funciones k1 y k2 y no las funciones mismas, esto implica que una vez que determinemos a las derivadas será necesario hacer una integración a cada una de ellas. Resolvamos el sistema.

Multipliquemos la ecuación (10) por dy2dx y la ecuación (11) por y2.

Label '12' multiply defined

Label '13' multiply defined

Si a la ecuación (12) le restamos la ecuación (13) obtenemos lo siguiente.

Label '14' multiply defined

Recordemos que el Wronskiano esta definido como

Label '15' multiply defined

Entonces la ecuación (14) la podemos escribir como

Label '16' multiply defined

Como y1 y y2 forman un conjunto fundamental de soluciones de la ecuación homogénea asociada, entonces

W(y1,y2)0

Así, de la ecuación (16) obtenemos que

Label '17' multiply defined

Hemos encontrado el valor de la derivada de la función k1(x), integrando obtenemos finalmente que

Label '18' multiply defined

En un proceso totalmente análogo, si multiplicamos a la ecuación (10) por dy1dx y a la ecuación (11) por y1 y realizamos los mismos pasos obtendremos la ecuación para la derivada de la función k2(x).

Label '19' multiply defined

Integrando obtendremos la función que buscamos

(20)k2(x)=y1(x)g(x)W(y1,y2)dx

Sustituyendo los resultados (18) y (20) en la solución particular (4) obtenemos finalmente la solución que buscábamos

(21)yp(x)=y1(x)y2(x)g(x)W(y1,y2)dx+y2(x)y1(x)g(x)W(y1,y2)dx

El resultado (21) corresponde a la solución particular de la ecuación diferencial (2) que, a diferencial del método de coeficientes indeterminados, se aplica para cualquier función g(x), aunque cabe mencionar que si la función g es muy compleja, entonces nos resultará, en algunas ocasiones, complicado resolver las integrales involucradas.

A lo largo del curso hemos motivado a no memorizar las formulas y en su lugar desarrollar el procedimiento del método correspondiente, sin embargo, en esta ocasión se trata de un método muy largo y complicado para usarse cada vez que se intente resolver una ecuación diferencial, por lo que se recomienda seguir los siguientes pasos.

  • Primero se determina la solución complementaria yc=c1y1+c2y2 de la ecuación diferencial homogénea asociada, esto nos permitirá determinar el conjunto fundamental de soluciones {y1,y2}.
  • Una vez conocido el conjunto fundamental de soluciones se procede a calcular el Wronskiano W(y1,y2).
  • Posteriormente se divide la ecuación diferencial por a2 para escribir la ecuación es su forma estándar (3) y así obtener la forma de la función g(x).
  • Se sustituyen los valores correspondientes en (18) y (20) para obtener las funciones k1 y k2 respectivamente.
  • Finalmente se sustituyen los resultados en la solución particular yp=k1y1+k2y2 y posteriormente en la solución general y=yc+yp

Cuando se calculan las integrales indefinidas (18) y (20) no es necesario considerar las constantes de integración. Para mostrar esto consideremos las constantes c3 y c4, tales que

y(x)=yc(x)+yp(x)=c1y1(x)+c2y2(x)+[k1(x)+c3]y1(x)+[k2(x)+c4]y2(x)=[c1+c3]y1(x)+[c2+c4]y2(x)+k1(x)y1(x)+k2(x)y2(x)=C1y1(x)+C2y2(x)+k1(x)y1(x)+k2(x)y2(x)

Es decir, las constantes de la solución complementaria contienen todas las constantes que puedan aparecer en el método.

Realicemos algunos ejemplos.

Ejemplo: Determinar la solución general de la ecuación diferencial

3d2ydx26dydx+6y=exsec(x)

Solución: El primer paso es obtener la solución complementaria. La ecuación auxiliar es

3k26k+6=0

De donde k1=1+i y k2=1i, identificamos que α=β=1, entonces la forma de la solución complementaria es

yc(x)=c1excos(x)+c2exsin(x)

El conjunto fundamental de soluciones esta conformado por las funciones

y1(x)=excos(x)yy2(x)=exsin(x)

La derivada de ambas soluciones son

dy1dx=excos(x)exsin(x)ydy2dx=exsin(x)+excos(x)

Usando estos resultados calculemos el Wronskiano.

W=|excos(x)exsin(x)excos(x)exsin(x)exsin(x)+excos(x)|=e2xcos(x)sin(x)+e2xcos2(x)e2xcos(x)sin(x)+e2xsin2(x)=e2x

El Wronskiano es

W(x)=e2x

¡Cuidado!, como en la ecuación diferencial la segunda derivada tiene un coeficiente, debemos dividir toda la ecuación por dicho coeficiente para obtener la forma estándar y así la función g(x). La ecuación en su forma estándar es

d2ydx22dydx+2y=exsec(x)3

En este caso la función g es

g(x)=exsec(x)3

Ahora que ya conocemos los valores que necesitábamos, recurrimos a las ecuaciones (18) y (20) para obtener las funciones que buscamos.

Para la función k1(x), tenemos lo siguiente.

k1(x)=y2(x)g(x)W(y1,y2)dx=(exsin(x))(exsec(x)3)e2xdx=e2xsin(x)sec(x)3e2xdx=13tan(x)dx=13ln|cos(x)|

La integral de la tangente es común. Por tanto, la función k1 es

k1(x)=13ln|cos(x)|

Para el caso de la función k2(x), tenemos lo siguiente.

k2(x)=y1(x)g(x)W(y1,y2)dx=(excos(x))(exsec(x)3)e2xdx=e2xcos(x)sec(x)3e2xdx=13dx=13x

La función k2 es

k2(x)=13x

Ya podemos establecer que la solución particular, de acuerdo a (4), es

yp(x)=13ln|cos(x)|[excos(x)]+13x[exsin(x)]

Por lo tanto, la solución general de la ecuación diferencial

3d2ydx26dydx+6y=exsec(x)

es

y(x)=c1excos(x)+c2exsin(x)+13excos(x)ln|cos(x)|+13xexsin(x)

◻

Con este ejemplo encontramos un buen momento para reflexionar y darnos cuenta de que ya hemos avanzado mucho, tan sólo observa el tipo de ecuación que acabamos de resolver y no sólo eso, observa que tan compleja es la solución general.

¡Sigamos adelante!.

Ejemplo: Determinar la solución general de la ecuación diferencial

d2ydx2+2dydx+y=exx

Solución: Como la ecuación ya está es su forma estándar la función g es

g(x)=exx

Determinemos la solución complementaria, la ecuación auxiliar es

k2+2k+1=0

De donde k1=k2=1, la multiplicidad de la solución nos indica que la forma de la solución complementaria es

yc(x)=c1ex+c2xex

El conjunto fundamental de soluciones esta conformado por las funciones

y1(x)=exyy2(x)=xex

Usando estas soluciones y sus derivadas calculemos el Wronskiano.

W=|exxexexexxex|=e2xxe2x+xe2x=e2x

El Wronskiano es

W(x)=e2x

Sustituyamos estos resultados directamente en la ecuación (21).

yp(x)=ex(xex)(exx)e2xdx+xex(ex)(exx)e2xdx=exdx+xexdxx=xex+xexln(x)

La solución particular es

yp(x)=xex+xexln(x)

Por lo tanto, la solución general de la ecuación diferencial es

y(x)=c1ex+c2xexxex+xexln(x)

◻

Un ejemplo más.

Ejemplo: Determinar la solución general de la ecuación diferencial

d2ydx2y=4x3ex

Solución: La función g es

g(x)=4x3ex

y la ecuación auxiliar es

k21=0

De donde k1=1 y k2=1. Entonces, la solución complementaria es

yc(x)=c1ex+c2ex

El conjunto fundamental de soluciones esta conformado por

y1(x)=exyy2(x)=ex

Usando estas funciones y sus derivadas calculemos el Wronskiano.

W=|exexexex|=2

El Wronskiano es W=2. Sustituyendo estos resultados directamente en la ecuación (21), obtenemos lo siguiente.

yp(x)=ex(ex)(4x3ex)2dx+ex(ex)(4x3ex)2dx=2exx3dx2exx3e2xdx=12x4ex2exx3e2xdx

La integral que nos falta se puede resolver por partes tomando u=x2 y v=e2x. Resolviendo la integral obtendremos lo siguiente.

x3e2xdx=12e2xx334e2xx2+34e2xx38e2x

Sustituyendo en la solución particular tenemos

yp(x)=12x4ex2ex(12e2xx334e2xx2+34e2xx38e2x)=12x4exx3ex+32x2ex32xex+34ex

Finalmente obtenemos como solución particular a la función

yp(x)=ex(12x4x3+32x232x+34)

Y por lo tanto, la solución general de la ecuación diferencial es

y(x)=c1ex+c2ex+ex(12x4x3+32x232x+34)

Este resultado es válido, sin embargo se puede simplificar más, ya que se puede reescribir a la solución como

y(x)=ex(c1+34)+c2ex+ex(12x4x3+32x232x)

y definir la constante C1=c1+34 para finalmente escribir la solución como

y(x)=C1ex+c2ex+ex(12x4x3+32x232x)

◻

Variación de parámetros en ecuaciones de orden superior

Este método se puede generalizar a ecuaciones de orden n aunque, por su puesto, los cálculos se vuelven más extensos.

A continuación mostraremos el panorama general para ecuaciones diferenciales de orden n y mostraremos los resultados para el caso n=3 que nos mostrará la forma en que aumenta la complejidad de los cálculos.

La ecuación de orden n es su forma estándar es

(22)dnydxn+Pn1(x)dn1ydxn1++P1(x)dydx+P0(x)y=g(x)

Si la solución complementaria de (22) es

(23)yc(x)=c1y1(x)+c2y2(x)++cnyn(x)

Entonces una solución particular debe ser

(24)yp(x)=k1(x)y1(x)+k2(x)y2(x)++kn(x)yn(x)

Análogo a las ecuaciones (10) y (11), las derivadas dkidx=ki con i=1,2,,n se determinan del sistema de n ecuaciones

y1k1+y2k2++ynkn=0y1k1+y2k2++ynkn=0(25)y1(n1)k1+y2(n1)k2++yn(n1)kn=g(x)

Al igual que el caso de segundo orden, las primeras n1 ecuaciones del sistema son suposiciones que se hacen para simplificar la ecuación resultante de sustituir la solución (24) en la ecuación (22).

Usando la regla de Cramer para resolver el sistema se obtiene que

(26)dkidx=WiW;i=1,2,,n

Donde W es el Wronskiano del conjunto fundamental {y1(x),y2(x),,yn(x)} y Wi es el determinante que se obtiene de remplazar la i-ésima columna del Wronskiano por la columna formada por el lado derecho de (25), es decir, la columna que consta de (0,0,,g(x)).

Para que quede más claro lo anterior, en el caso n=3 las dkidx, i=1,2,3 quedan como

(27)dk1dx=W1W,dk2dx=W2W,dk3dx=W3W

Donde

(28)W=|y1y2y3y1y2y3y1y2y3|

y

W1=|0y2y30y2y3g(x)y2y3|,W2=|y10y3y10y3y1g(x)y3|,W3=|y1y20y1y20y1y2g(x)|

Habrá que integrar las ecuaciones de (27) para obtener las funciones ki, i=1,2,3 y así obtener la solución particular

(29)yp(x)=k1(x)y1(x)+k2(x)y2(x)+k3(x)y3(x)

Notemos que usando esta notación, los resultados (17) y (19) del caso n=2 se pueden escribir como

(30)dk1dx=W1W=y2g(x)Wydk2dx=W2W=y1g(x)W

Donde

(31)W=|y1y2y1y2|,W1=|0y2g(x)y2|yW2=|y10y1g(x)|

Realicemos un ejemplo con una ecuación de orden 3.

Ejemplo: Resolver la ecuación diferencial de tercer orden

d3ydx3+dydx=tan(x)

Solución: La función g es

g(x)=tan(x)

y la ecuación auxiliar es

k3+k=k(k2+1)=0

De donde k1=0, k2=i y k3=i. Dos raíces son complejas conjugadas con α=0 y β=1. La primer raíz nos indica que la forma de una solución es

y1(x)=ek1x=1

mientras que las dos raíces restantes nos indican dos solución de la forma

y2(x)=eαxcos(βx)=cos(x)

y

y3(x)=eαxsin(βx)=sin(x)

Por lo tanto, la solución complementaria de la ecuación diferencial es

yc(x)=c1+c2cos(x)+c3sin(x)

Como vimos, el conjunto fundamental de soluciones es {y1=1,y2=cos(x),y3=sin(x)}, las primeras y segundas derivadas correspondientes son

dy1dx=0,dy2dx=sin(x),dy3dx=cos(x)

y

d2y1dx2=0,d2y2dx2=cos(x),d2y3dx2=sin(x)

Ahora calculemos los determinantes correspondientes, el primero de ellos es el Wronskiano

W=|1cos(x)sin(x)0sin(x)cos(x)0cos(x)sin(x)|=sin2(x)+cos2(x)=1

Para el resto de determinantes, tenemos

W1=|0cos(x)sin(x)0sin(x)cos(x)tan(x)cos(x)sin(x)|=tan(x)[cos2(x)+sin2(x)]=tan(x)

W2=|10sin(x)00cos(x)0tan(x)sin(x)|=cos(x)tan(x)=sin(x)

y

W3=|1cos(x)00sin(x)00cos(x)tan(x)|=sin(x)tan(x)=sin2(x)cos(x)

Sustituyendo estos resultados en (27), obtenemos

dk1dx=tan(x),dk2dx=sin(x),dk3dx=sin2(x)cos(x)

Procedemos a integrar cada ecuación (sin considerar constantes) para obtener las funciones que buscamos.

La primer integral es común,

k1(x)=tan(x)dx=ln|cos(x)|

La segunda integral es directa

k2(x)=sin(x)dx=cos(x)

Mientras que para la tercer integral si se requiere de un mayor cálculo.

k3(x)=sin2(x)cos(x)dx=1cos2(x)cos(x)dx=1cos(x)dx+cos(x)dx

Por un lado

1cos(x)dx=sec(x)dx=ln|tan(x)+sec(x)|

por otro lado,

cos(x)dx=sin(x)

entonces

k3(x)=ln|tan(x)+sec(x)|+sin(x)

Ahora que conocemos las funciones incógnita concluimos que la solución particular de la ecuación diferencial es

yp(x)=ln|cos(x)|(1)+cos(x)[cos(x)]+[ln|tan(x)+sec(x)|+sin(x)](sin(x))=ln|cos(x)|+cos2(x)sin(x)ln|tan(x)+sec(x)|+sin2(x)=ln|cos(x)|sin(x)ln|tan(x)+sec(x)|

Por lo tanto, la solución general de la ecuación diferencial de tercer orden es

y(x)=c1+c2cos(x)+c3sin(x)ln|cos(x)|sin(x)ln|tan(x)+sec(x)|

◻

Como podemos notar, los cálculos se hacen más extensos, sin embargo los pasos a seguir son los mismos para cualquier orden.

El método de variación de parámetros, a diferencia del método de coeficientes indeterminados, tiene la ventaja de que siempre produce una solución de la ecuación diferencial independientemente de la forma de la función g(x), siempre y cuando se pueda resolver la ecuación homogénea asociada. Además, el método de variación de parámetros es aplicable a ecuaciones diferenciales lineales con coeficientes variables.

Tarea moral

Los siguientes ejercicios no forman parte de la evaluación del curso, pero servirán para entender mucho mejor los conceptos vistos en esta entrada, así como temas posteriores.

  1. Obtener la solución general de las siguientes ecuaciones diferenciales.
  • d2ydx29dydx=18x2e9x
  • d2ydx2+9y=18exsin(x)
  • 4d2ydx24dydx+y=ex/21x2
  1. Resolver las siguientes ecuaciones diferenciales para las condiciones iniciales dadas.
  • d2ydx23dydx=12e4x(x+1);y(0)=0,y(0)=4
  • d2ydx2+dydx6y=10exsin(x);y(0)=217,y(0)=0
  1. Obtener la solución general de las siguientes ecuaciones diferenciales de tercer orden. Simplificar la forma de la solución redefiniendo las constantes.
  • d3ydx3+4dydx=sec(2x)
  • d3ydx3+d2ydx2dydxy=2ex

Más adelante…

Hemos concluido con el estudio de las ecuaciones diferenciales lineales de orden superior con coeficientes constantes.

Lo que sigue es estudiar este mismo tipo de ecuaciones, pero en el caso en el que los coeficientes no son constantes, es decir, son coeficientes variables. Estas ecuaciones suelen ser mucho más difícil de resolver, sin embargo existe un tipo de ecuación especial, conocida como ecuación de Cauchy – Euler, que contiene coeficientes variables, pero que su método de resolución es bastante similar a lo que hemos desarrollado en el caso de coeficientes constantes, pues su resolución involucra resolver una ecuación auxiliar. En la siguiente entrada estudiaremos dicha ecuación.

Entradas relacionadas

Agradecimientos

Trabajo realizado con el apoyo del Programa UNAM-DGAPA-PAPIME PE104522 «Hacia una modalidad a distancia de la Licenciatura en Matemáticas de la FC-UNAM – Etapa 2»

Ecuaciones Diferenciales I – Videos: Propiedades del conjunto de soluciones a un sistema lineal de ecuaciones de primer orden

Por Eduardo Vera Rosales

Introducción

En la entrada anterior comenzamos el estudio de los sistemas de ecuaciones diferenciales de primer orden x˙1=F1(t,x1,x2,,xn)x˙2=F2(t,x1,x2,,xn)x˙n=Fn(t,x1,x2,,xn) donde revisamos las principales definiciones y enunciamos el teorema de existencia y unicidad correspondiente a sistemas de primer orden y sus problemas de condición inicial. Es momento ahora de estudiar las principales propiedades que cumple el conjunto de soluciones a un sistema lineal de ecuaciones de primer orden, las cuales se comportan de una manera bastante similar al conjunto de soluciones a una ecuación de segundo orden lineal que revisamos en la unidad anterior.

Iniciaremos revisando al conjunto de soluciones al sistema lineal homogéneo X˙=AX el cual cumple el principio de superposición, es decir, si tenemos n soluciones, digamos X1(t),X2(t),,Xn(t), entonces cualquier combinación lineal de estas también lo será. Si recuerdas tus cursos de Álgebra Lineal, esta última propiedad nos dice que el conjunto de soluciones es cerrado bajo la suma y producto por escalar usuales definidos para matrices. Con estas operaciones, veremos que el conjunto de soluciones al sistema lineal homogéneo forma un espacio vectorial.

Posteriormente definiremos el Wronskiano de un subconjunto de soluciones al sistema lineal homogéneo, el cual es similar más no igual al Wronskiano que definimos para ecuaciones lineales de segundo orden. En la tarea moral demostrarás la relación que tienen estos dos Wronskianos.

Si hablamos del Wronskiano y del conjunto de soluciones como un espacio vectorial, debemos hablar también de dependencia e independencia lineal entre las soluciones al sistema. Además, demostraremos que si el Wronskiano no se anula entonces el subconjunto de soluciones es linealmente independiente. Además si lo último ocurre podremos expresar cualquier solución como una combinación lineal de las soluciones linealmente independientes. Con estos conceptos podremos definir a la matriz fundamental de soluciones del sistema, la cual revisaremos más a detalle en entradas posteriores.

Terminaremos revisando el caso no homogéneo X˙=AX+Q demostrando que su solución general será la suma de la solución general al sistema homogéneo y una solución particular al sistema no homogéneo.

El espacio vectorial del conjunto de soluciones a un sistema lineal homogéneo

En el primer video probamos el principio de superposición de soluciones al sistema lineal homogéneo. Además, vemos que el conjunto de soluciones al sistema forma un espacio vectorial con la suma y producto por escalar usuales para matrices.

El Wronskiano de un subconjunto de soluciones e independencia lineal

Definimos el Wronskiano de un subconjunto de soluciones al sistema lineal homogéneo, así como los conceptos de dependencia e independencia lineal de soluciones. Probamos un importante teorema que relaciona estos dos conceptos y nos dice cómo se ve la solución general al sistema. Finalizamos definiendo la matriz fundamental de soluciones del sistema.

Solución general al sistema lineal no homogéneo

Finalizamos la entrada demostrando que la solución general al sistema lineal no homogéneo es la suma de la solución general al sistema homogéneo y una solución particular al sistema no homogéneo.

Tarea moral

Los siguientes ejercicios no forman parte de la evaluación del curso, pero te servirán para entender mucho mejor los conceptos vistos en esta entrada, así como temas posteriores.

  • ¿El conjunto de soluciones a un sistema lineal no homogéneo forma un espacio vectorial con las operaciones usuales de matrices?
  • Prueba que X1(t)=(100);X2(t)=(t20);X3(t)=(t2t0) son linealmente independientes en R.
  • Sean X1(t),X2(t),,Xn(t) soluciones al sistema X˙=AX en el intervalo [a,b]. Demuestra que W[X1,X2,,Xn](t)=0t[a,b], ó W[X1,X2,,Xn](t)0t[a,b].
  • Considera el sistema lineal X˙=(0110)X. Prueba que X1(t)=(etet);X2(t)=(etet) son soluciones al sistema. Además prueba que son linealmente independientes en R y por lo tanto forma una matriz fundamental de soluciones al sistema.
  • Considera la ecuación y¨+p(t)y˙+q(t)y=0 y su sistema de ecuaciones correspondiente X˙=(01q(t)p(t))X. Prueba que si X1(t), X2(t) son soluciones linealmente independientes al sistema de ecuaciones, y si y1(t), y2(t) forman un conjunto fundamental de soluciones a la ecuación de segundo orden, entonces se satisface la identidad W[y1,y2](t)=cW[X1,X2](t) para alguna constante c0.

Más adelante

En la siguiente entrada comenzaremos a resolver algunos sistemas lineales bastante sencillos. El método que estudiaremos será el de eliminación de variables, el cual consiste en eliminar variables dependientes hasta quedarnos con una ecuación diferencial de orden superior. Resolviendo esta última ecuación podremos encontrar la solución general al sistema original. Este método funciona para sistemas lineales con coeficientes constantes.

¡Hasta la próxima!

Entradas relacionadas

Agradecimientos

Trabajo realizado con el apoyo del Programa UNAM-DGAPA-PAPIME PE104522 «Hacia una modalidad a distancia de la Licenciatura en Matemáticas de la FC-UNAM – Etapa 2»

Ecuaciones Diferenciales I: Soluciones a ecuaciones diferenciales de orden superior

Por Omar González Franco

Las matemáticas son la música de la razón.
– James Joseph Sylvester

Introducción

En la entrada anterior comenzamos a estudiar los problemas con valores iniciales (PVI) y problemas con valores en la frontera (PVF), ambos para el caso de las ecuaciones diferenciales lineales de orden superior. Vimos también que si y1,y2,,yk son k soluciones de una ecuación homogénea de n-ésimo orden en un intervalo δ, entonces la combinación lineal

Label '1' multiply defined

donde las ci, i=1,2,,k son constantes, también es solución en el intervalo δ, este resultado es conocido como principio de superposición y nuestro propósito en esta entrada es estudiar las propiedades de todas estas soluciones donde la independencia lineal de las funciones jugará un papel muy importante en la construcción del conjunto fundamental de soluciones y de la solución general.

Es importante tener presente el concepto de conjunto fundamental de soluciones presentado en la entrada anterior.

Soluciones a ecuaciones diferenciales lineales de orden superior

Estamos interesados en soluciones linealmente independientes de una ecuación diferencial lineal homogénea de orden superior de la forma

Label '2' multiply defined

Al intentar responder la pregunta de si el conjunto de n soluciones {y1,y2,,yn} de (2) es linealmente independiente podemos apelar directamente a la definición de independencia lineal, sin embargo esta pregunta se puede responder de una forma mecánica usando un determinante llamado el Wronskiano.

El Wronskiano es una herramienta que podemos utilizar para determinar si el conjunto de soluciones de la ecuación diferencial (2) es un conjunto linealmente independiente y la forma de hacerlo es a través del siguiente teorema conocido como criterio para soluciones linealmente independientes.

Este teorema nos dice que sólo basta mostrar que el Wronskiano es distinto de cero para garantizar que el conjunto de soluciones {y1,y2,,yn} es linealmente independiente y por tanto formará un conjunto fundamental de soluciones.

Al conjunto de soluciones linealmente independiente {y1,y2,,yn} de la ecuación (2) se le denomina fundamental porque, así como cualquier vector en R3 se puede expresar como una combinación lineal de los vectores linealmente independientes i^,j^ y k^, cualquier solución de una ecuación diferencial de la forma (2) se puede expresar como una combinación lineal de las n soluciones del conjunto fundamental, podemos decir que las soluciones {y1,y2,,yn} son los bloques básicos para la solución general de la ecuación.

En el siguiente teorema se enuncia la forma general de la solución de la ecuación diferencial (2).

Aterricemos estas ideas generales al caso de las ecuaciones diferenciales lineales homogéneas de segundo orden.

Ecuaciones lineales homogéneas de segundo orden

Una ecuación diferencial lineal homogénea de segundo orden es de la forma

Label '6' multiply defined

Sobre esta ecuación desarrollaremos la siguiente teoría. Primero definamos el Wronskiano para el caso n=2.

Ahora que conocemos la forma del Wronskiano para n=2, demostremos el teorema de la solución general para el caso n=2.

Demostración: Sea y(x) una solución de la ecuación diferencial (6) en el intervalo δ y sea x0δ, tal que

y(x0)=αydydx(x0)=β$

con α y β constantes. Supongamos que existen c1 y c2 constantes tales que

Label '8' multiply defined

y

Label '9' multiply defined

esto debido a que por hipótesis y1(x) y y2(x) son soluciones de la ecuación diferencial y por tanto la combinación lineal también lo será. Aplicando el teorema de existencia y unicidad obtenemos que la solución y(x) tiene que ser de la forma

y(x)=c1y1+c2y2

por lo que nuestro problema se reduce a demostrar que las constantes c1 y c2 existen.

Si multiplicamos a la ecuación (8) por dy2dx(x0) y a la ecuación (9) por y2(x0) obtenemos lo siguiente, respectivamente.

Label '10' multiply defined

y

Label '11' multiply defined

Restémosle a la ecuación (10) la ecuación (11).

Label '12' multiply defined

Sabemos que el Wronskiano, en x=x0, está definido como

Label '13' multiply defined

Por comodidad denotaremos a W(y1(x0),y2(x0)) como W(x0). Entonces la ecuación (12) se puede escribir de la siguiente manera.

Label '14' multiply defined

Debido a que por hipótesis W(y1,y2)0 para toda xδ, en particular lo es en x=x0, por tanto podemos despejar a la constante c1 y así obtener un valor para dicha constante lo que muestra su existencia.

Para obtener la expresión de c2 hacemos algo similar, multiplicamos a la ecuación (8) por dy1dx(x0) y a la ecuación (9) por y1(x0) y repetimos el mismo procedimiento demostrando así que existe un valor para la constante c2.

Como hemos encontrado valores para c1 y c2, entonces existen y por lo tanto la solución general a la ecuación (4) es

Label '15' multiply defined

◻

Ya hemos definido lo que es el conjunto fundamental de soluciones de una ecuación diferencial de orden n, para el caso n=2 lo podemos definir de la siguiente manera.

Así, si encontramos un conjunto fundamental de soluciones {y1(x),y2(x)}, entonces

W(y1,y2)0

para toda xδ y por tanto

y(x)=c1y1(x)+c2y2(x)

será la solución general de la ecuación diferencial (4).

Del criterio para soluciones linealmente independientes se puede hacer notar que cuando y1,y2,,yn son n soluciones de la ecuación diferencial (2) en un intervalo δ, el Wronskiano W(y1,y2,,yn) es siempre igual a cero o nunca es cero en todo δ. Vamos a demostrar este hecho para el caso n=2.

Demostración: Como y1(x) y y2(x) son soluciones de la ecuación (6), entonces

Label '16' multiply defined

y

Label '17' multiply defined

Si multiplicamos a la ecuación (16) por y2 y a la ecuación (17) por y1 obtenemos lo siguiente, respectivamente.

Label '18' multiply defined

y

Label '19' multiply defined

A la ecuación (19) vamos a restarle la ecuación (18).

Label '20' multiply defined

Sabemos que

W(y1,y2)=y1dy2dxy2dy1dx

y notemos lo siguiente

dWdx=ddx(y1dy2dxy2dy1dx)=dy1dxdy2dx+y1d2y2dx2dy2dxdy1dxy2d2y1dx2=y1d2y2dx2y2d2y1dx2

Es decir,

Label '21' multiply defined

En términos del Wronskiano la ecuación (20) se puede escribir como

Label '22' multiply defined

Como a2(x)0 para toda xδ, entonces podemos definir la función

P(x)=a1(x)a2(x)

tal que la ecuación (22) se pueda escribir como

Label '23' multiply defined

Esta resultado corresponde a una ecuación diferencial lineal homogénea de primer orden y ya sabemos que la solución es de la forma

W(x)=keP(x)dx

de manera que hay dos posibilidades:

  • Si k=0W(x)=0,xδ
  • Si k0W(x)0,xδ

◻

El criterio para soluciones linealmente independientes nos garantiza que si el Wronskiano es distinto de cero, entonces el conjunto de soluciones es linealmente independiente en δ, lo opuesto es cierto bajo ciertas condiciones, si el Wronskiano es igual a cero, entonces el conjunto de soluciones es linealmente dependiente. Demostremos este hecho.

Demostración: Por hipótesis

W(y1(x),y2(x))=0

xδ, es decir

Label '24' multiply defined

Consideremos el siguiente resultado.

Label '25' multiply defined

Donde hemos considerado la hipótesis y20. Si usamos la hipótesis (24) obtenemos que

ddx(y1y2)=0

xδ, integrando esta ecuación obtenemos que

y1y2=k

O bien,

Label '26' multiply defined

con k una constante. Esto demuestra que y1 y y2 son linealmente dependientes.

◻

Hay que tener muy presentes las hipótesis de este teorema, pues es posible que el Wronskiano sea cero aún cuando las funciones consideradas en un cierto intervalo sean linealmente independientes en él.

Como consecuencia del teorema anterior podemos establecer el criterio para soluciones linealmente independientes en el caso n=2.

Realicemos un ejemplo.

Ejemplo: En la entrada anterior de tarea moral tenías que verificar que las funciones

y1(x)=e3xyy2(x)=e4x

forman un conjunto fundamental de soluciones de la ecuación diferencial

d2ydx2dydx12y=0

en δ=(,). Demostremos esto mismo usando los teoremas vistos anteriormente.

Solución: Consideremos las soluciones

y1(x)=e3xyy2(x)=e4x

y sus respectivas derivadas

dy1dx=3e3xydy2dx=4e4x

Calculemos el Wronskiano.

W(y1,y2)=|e3xe4x3e3x4e4x|=e3x(4e4x)e4x(3e3x)=7ex0

Como

W(y1,y2)=7ex0

entonces {y1(x)=e3x,y2(x)=e4x} forma un conjunto fundamental de soluciones y la solución general está dada por

y(x)=c1e3x+c2e4x

◻

Con esto concluimos el estudio de algunas propiedades importantes de las soluciones a la ecuación diferencial lineal homogénea de orden superior, terminemos esta entrada con el estudio del caso no homogéneo.

Ecuaciones no homogéneas

La ecuación diferencial lineal no homogénea de n-ésimo orden es

Label '27' multiply defined

Nuestro objetivo es obtener la forma general de la solución de la ecuación no homogénea (27) y estudiar algunas propiedades de las soluciones.

Si recordamos al operador polinomial

Label '28' multiply defined

la definición anterior implica que

Label '29' multiply defined

Veamos el siguiente resultado.

Demostración: Sea y(x) la combinación lineal

y(x)=c1y1+c2y2++cnyn+yp

si aplicamos el operador polinomial, tenemos

L{y(x)}=L{c1y1+c2y2++cnyn+yp}=c1L{y1}+c2L{y2}++cnL{yn}+L{yp}=0+g(x)=g(x)

Ya que L{yi}=0 para cada i=1,2,,n por ser cada yi solución de la ecuación homogénea, mientras que L{yp}=g(x) por ser solución de la ecuación no homogénea. Entonces, como

L{y(x)}=g(x)

concluimos que la combinación lineal

y(x)=c1y1+c2y2++cnyn+yp

es solución de la ecuación diferencial no homogénea.

◻

¿Y qué ocurre si las soluciones y1,y2,,yn forman un conjunto fundamental de soluciones?. La respuesta es que la combinación lineal

y(x)=c1y1+c2y2++cnyn+yp

sería la solución general de la ecuación diferencial no homogénea (27). Demostremos este resultado.

Demostración: Sea y(x) la solución general de la ecuación no homogénea (27) y sea yp(x) una solución particular de la misma ecuación, ambas definidas en el intervalo δ, de manera que

L{y(x)}=L{yp(x)}=g(x)

con L el operador polinomial (28). Nuestro objetivo es encontrar la forma explícita de y(x).

Definamos la función

(32)h(x)=y(x)yp(x)

y notemos lo siguiente.

L{h(x)}=L{y(x)yp(x)}=L{y(x)}L{yp(x)}=g(x)g(x)=0

Esto es,

L{h(x)}=0

lo que significa que la función h(x) es solución de la ecuación homogénea (2) y por el teorema de la solución general de ecuaciones homogéneas podemos establecer que la función h(x) tiene la siguiente forma.

(33)h(x)=c1y1(x)+c2y2(x)++cnyn(x)

Con {y1,y2,,yn} un conjunto fundamental de soluciones. Sustituyendo (33) en (32) y despejando a la solución general y(x) obtenemos finalmente que

y(x)=c1y1(x)+c2y2(x)++cnyn(x)+yp(x)

que es lo que queríamos demostrar.

◻

La diferencia entre las soluciones (30) y (31) es que en (31) las yi,i=1,2,,n forman un conjunto fundamental de soluciones, es decir, son linealmente independientes entre sí, mientras que en (30) no necesariamente forman una conjunto fundamental y sin embargo, también son solución de la ecuación (27).

En el caso de las ecuaciones no homogéneas vemos que la solución general corresponde a la suma de la solución general de la ecuación homogénea asociada más una solución particular de la ecuación no homogénea. En este caso no homogéneo la solución general de la ecuación homogénea tiene un nombre particular.

Por tanto, resolver una ecuación lineal no homogénea implica resolver primero la ecuación homogénea asociada para obtener la función complementaria yc(x) y luego se encuentra una solución particular yp(x) de la ecuación no homogénea para finalmente sumarlas

(34)y(x)=yc(x)+yp(x)

Realicemos un ejemplo.

Ejemplo: Probar que la función

y(x)=c1e2x+c2xe2x+x2e2x+x2

definida en el intervalo δ=(,), es la solución general de la ecuación diferencial

d2ydx24dydx+4y=2e2x+4x12

Solución: Primero probemos que las funciones

y1=e2xyy2=xe2x

forman un conjunto fundamental de soluciones de la ecuación homogénea asociada

d2ydx24dydx+4y=0

Para ello veamos que y1 y y2 son soluciones de la ecuación homogénea y que son linealmente independientes, es decir, que W(y1,y2)0. Calculemos las derivadas.

y1=e2xdy1dx=2e2xd2y1dx2=4e2x

y2=xe2xdy2dx=e2x+2xe2xd2y2dx2=4e2x+4xe2x

De tarea moral muestra que ambas funciones son solución de la ecuación homogénea asociada, es decir, que

d2y1dx24dy1dx+4y1=0yd2y2dx24dy2dx+4y2=0

Ahora probemos que forman un conjunto fundamental de soluciones, para ello calculemos el Wronskiano.

W(y1,y2)=|e2xxe2x2e2xe2x+2xe2x|=e2x(e2x+2xe2x)xe2x(2e2x)=e4x0

Como W(y1,y2)0, xδ, por los teoremas vistos anteriormente concluimos que {y1=e2x,y2=xe2x} forma un conjunto fundamental de soluciones de la ecuación homogénea asociada y que la solución general a dicha ecuación es

yc(x)=c1e2x+c2xe2x

donde el subíndice c indica que es la función complementaria.

Ahora verifiquemos que la función

yp(x)=x2e2x+x2

es una solución particular de la ecuación no homogénea. Calculemos la primera y segunda derivada.

dypdx=2xe2x+2x2e2x+1

d2ypdx2=2e2x+8xe2x+4x2e2x

Sustituyamos en la ecuación diferencial.

d2ydx24dydx+4y=(2e2x+8xe2x+4x2e2x)4(2xe2x+2x2e2x+1)+4(x2e2x+x2)=2e2x+(8xe2x8xe2x)+(4x2e2x8x2e2x+4x2e2x)+4x12=2e2x+4x12

Esto es,

d2ydx24dydx+4y=2e2x+4x12

que justo corresponde a la ecuación diferencial no homogénea, por lo tanto, efectivamente yp es una solución particular.

Como {y1=e2x,y2=xe2x} es un conjunto fundamental de soluciones de la ecuación homogénea asociada y

yp(x)=x2e2x+x2

es una solución particular de la ecuación no homogénea, por el teorema de la solución general de ecuaciones no homogéneas concluimos que la función

y(x)=c1e2x+c2xe2x+x2e2x+x2

es la solución general de la ecuación no homogénea.

◻

Hay algo muy interesante que ocurre en el ejemplo anterior. Mostramos que la función

yp(x)=x2e2x+x2

es una solución particular de la ecuación no homogénea

d2ydx24dydx+4y=2e2x+4x12=g(x)

Sin embargo, si haces los cálculos correspondientes notarás que la función

yp1(x)=x2e2x

es una solución particular de la ecuación

d2ydx24dydx+4y=2e2x=g1(x)

mientras que la función

yp2(x)=x2

es una solución particular de la ecuación

d2ydx24dydx+4y=4x12=g2(x)

Así, si superponemos las soluciones particulares

yp(x)=yp1(x)+yp2(x)

obtenemos en la ecuación diferencial la superposición de la funciones

g(x)=g1(x)+g2(x)

Lo anterior es efecto del principio de superposición para ecuaciones no homogéneas.

Demostración: Sea L el operador polinomial (28) y sean ypi(x), i=1,2,,k, soluciones particulares de las ecuaciones no homogéneas

L{ypi(x)}=gi(x)

i=1,2,,k respectivamente. Definamos la función

yp(x)=yp1(x)+yp2(x)++ypk(x)

Nuestro objetivo es demostrar que la función yp(x) es una solución particular de la ecuación (37), es decir, que se cumple que

L{yp(x)}=g1(x)+g2(x)++gk(x)

En efecto

L{yp(x)}=L{yp1(x)+yp2(x)++ypk(x)}=L{yp1(x)}+L{yp2(x)}++L{ypk(x)}=g1(x)+g2(x)++gk(x)

Con esto queda probado que

yp(x)=yp1(x)+yp2(x)++ypk(x)

es solución de (37).

◻

Realicemos un último ejemplo.

Ejemplo: Probar que

  • yp1(x)=4x2 es solución particular de d2ydx23dydx+4y=16x2+24x8,
  • yp2(x)=e2x es solución particular de d2ydx23dydx+4y=2e2x,
  • yp3(x)=xex es solución particular de d2ydx23dydx+4y=2xexex.

y probar que la superposición

y(x)=yp1(x)+yp2(x)+yp3(x)=4x2+e2x+xex

es una solución de

d2ydx23dydx+4y=16x2+24x8+2e2x+2xexex

Solución: Sean

g1(x)=16x2+24x8,g2(x)=2e2xyg3(x)=2xexex

De tarea moral muestra que efectivamente,

d2yp1dx23dyp1dx+4yp1=g1(x)

d2yp2dx23dyp2dx+4yp2=g2(x)

d2yp3dx23dyp3dx+4yp3=g3(x)

Por el principio de superposición para ecuaciones no homogéneas sabemos que la función

y(x)=yp1(x)+yp2(x)+yp3(x)

es solución de la ecuación

d2ydx23dydx+4y=g1(x)+g2(x)+g3(x)

Por lo tanto, la función

y(x)=4x2+e2x+xex

es solución de la ecuación diferencial

d2ydx23dydx+4y=16x2+24x8+2e2x+2xexex

Si gustas puedes calcular la primera y segunda derivada de y(x) y verificar la ecuación anterior para asegurarte del resultado.

◻

Con esto concluimos nuestro estudio sobre algunas propiedades de las ecuaciones diferenciales lineales de orden superior. En la siguiente entrada conoceremos un primer método para resolver ecuaciones diferenciales de segundo orden.

Tarea moral

Los siguientes ejercicios no forman parte de la evaluación del curso, pero servirán para entender mucho mejor los conceptos vistos en esta entrada, así como temas posteriores.

  1. Dadas las soluciones de las siguientes ecuaciones diferenciales lineales homogéneas de orden superior en el intervalo dado, calcular el Wronskiano para determinar si es un conjunto fundamental de soluciones y en caso de serlo dar la solución general.
  • x3d3ydx3+6x2d2ydx2+4xdydx4y=0, con soluciones

y1=x,y2=1x2,y3=1x2ln(x);δ=(0,).

  • d4ydx4+d2ydx2=0, con soluciones

y1=1,y2=x,y3=cos(x),y4=sin(x);δ=(,).

  1. Dadas las soluciones de las siguientes ecuaciones diferenciales lineales no homogéneas de orden superior en el intervalo dado, probar que se trata de la solución general de la ecuación.
  • d2ydx27dydx+10y=24ex, con solución

y(x)=c1e2x+c2e5x+6ex;δ=(,).

  • 2x2d2ydx2+5xdydx+y=x2x, con solución

y(x)=c11x+c21x+115x216x;δ=(0,).

  1. Comprobar que las funciones yp1(x)=3e2xyyp2(x)=x2+3x son, respectivamente, soluciones particulares de d2ydx26dydx+5y=9e2xyd2ydx26dydx+5y=5x2+3x16
  1. Usando el ejercicio anterior, encontrar la solución particular de las siguientes ecuaciones.
  • d2ydx26dydx+5y=5x2+3x169e2x
  • d2ydx26dydx+5y=10x26x+32+e2x

Más adelante…

Ahora que ya conocemos algunas propiedades de las ecuaciones diferenciales de orden superior y sus soluciones, en particular de las ecuaciones lineales de segundo orden, es momento de comenzar a desarrollar los distintos métodos de resolución de estas ecuaciones diferenciales.

En la siguiente entrada comenzaremos con un método que permite reducir una ecuación de segundo orden en una ecuación de primer orden, de tal manera que podremos resolverla aplicando alguno de los métodos vistos en la unidad anterior. No es casualidad que dicho método se conozca como método de reducción de orden.

Entradas relacionadas

Agradecimientos

Trabajo realizado con el apoyo del Programa UNAM-DGAPA-PAPIME PE104522 «Hacia una modalidad a distancia de la Licenciatura en Matemáticas de la FC-UNAM – Etapa 2»

Ecuaciones Diferenciales I – Videos: Ecuaciones lineales homogéneas de segundo orden. Propiedades del conjunto de soluciones

Por Eduardo Vera Rosales

Introducción

Hola a todos. Después de haber estudiado ecuaciones diferenciales de primer orden, llegamos a la segunda unidad del curso donde analizaremos ecuaciones diferenciales de segundo orden. Dada la dificultad para resolver este tipo de ecuaciones, nos enfocaremos únicamente en las ecuaciones lineales de segundo orden, es decir, de la forma a0(t)d2ydt2+a1(t)dydt+a2(t)y=g(t).

En esta entrada comenzaremos con el caso de las ecuaciones homogéneas de segundo orden, es decir, cuando g(t) es la función constante cero en un intervalo (α,β). Estudiaremos la teoría de las soluciones a este tipo de ecuaciones antes de analizar las distintas técnicas para resolverlas. Debido a que el conjunto de soluciones a este tipo de ecuaciones se comportan de buena manera, podremos encontrar la solución general a la ecuación si previamente conocemos dos soluciones particulares que cumplan una hipótesis que daremos a conocer en el intervalo (α,β). Definiremos el Wronskiano y la independencia lineal de dos soluciones a una ecuación diferencial, y probaremos distintos teoremas y propiedades de las soluciones con base en estos conceptos.

¡Comencemos!

Ecuaciones lineales homogéneas de segundo orden, Teorema de existencia y unicidad y solución general

En este video damos una introducción a las ecuaciones diferenciales de segundo orden, y en particular, a las ecuaciones lineales de segundo orden. Enunciamos el teorema de existencia y unicidad para ecuaciones lineales de segundo orden, y comenzamos a desarrollar la teoría para encontrar la solución general a ecuaciones homogéneas.

Conjunto fundamental de soluciones y el Wronskiano

Continuando con la teoría de las soluciones a ecuaciones homogéneas de segundo orden, demostramos un par de teoremas que nos ayudan a encontrar la solución general a este tipo de ecuaciones. Además, definimos al conjunto fundamental de soluciones de la ecuación homogénea y el Wronskiano de dos funciones.

Independencia lineal de soluciones

En este último video definimos el concepto de independencia lineal de soluciones a la ecuación homogénea de segundo orden, y demostramos un teorema que nos da otra forma de encontrar un conjunto fundamental de soluciones a la ecuación diferencial homogénea.

Tarea moral

Los siguientes ejercicios no forman parte de la evaluación del curso, pero te servirán para entender mucho mejor los conceptos vistos en esta entrada, así como temas posteriores.

  • Prueba que y1(t)=sint y y2(t)=cost son soluciones a la ecuación diferencial d2ydt2+y=0. Posteriormente prueba que y(t)=k1sint+k2cost también es solución a la ecuación, donde k1, k2 son constantes.
  • Prueba que {sint,cost} es un conjunto fundamental de soluciones a la ecuación del ejercicio anterior. ¿En qué intervalo es el conjunto anterior un conjunto fundamental de soluciones?
  • Prueba que si p(t), q(t) son continuas en (α,β), y1(t), y2(t) son soluciones a la ecuación d2ydt2+p(t)dydt+q(t)y=0 en (α,β) y existe t0 en dicho intervalo, donde W[y1,y2](t0)0, entonces {y1(t),y2(t)} forman un conjunto fundamental de soluciones en (α,β).
  • Prueba que si p(t), q(t) son continuas en (α,β), entonces existe un conjunto fundamental de soluciones {y1(t),y2(t)} a la ecuación d2ydt2+p(t)dydt+q(t)y=0 en el mismo intervalo. (Hint: Toma un punto en el intervalo (α,β) y dos problemas de condición inicial adecuados de tal forma que puedas utilizar el teorema de existencia y unicidad y el Wronskiano para deducir el resultado).
  • Prueba que y1(t)=t|t|, y2(t)=t2 son linealmente independientes en [1,1] pero linealmente dependientes en [0,1]. Verifica que el Wronskiano se anula en R. ¿Pueden ser y1(t), y2(t) soluciones a la ecuación d2ydt2+p(t)dydt+q(t)y=0 en (1,1) si p y q son continuas en este intervalo?

Más adelante

En la próxima entrada conoceremos el método de reducción de orden, donde supondremos que ya conocemos una solución particular y1(t) a la ecuación lineal homogénea de segundo orden, y con ayuda de esta hallaremos una segunda solución y2(t) tal que forma un conjunto fundamental de soluciones junto con y1.

Entradas relacionadas

Agradecimientos

Trabajo realizado con el apoyo del Programa UNAM-DGAPA-PAPIME PE104522 «Hacia una modalidad a distancia de la Licenciatura en Matemáticas de la FC-UNAM – Etapa 2»