Ecuación diferencial lineal

Un ejemplo de ello.

En matemáticas, una ecuación diferencial lineal es aquella ecuación diferencial cuyas soluciones pueden obtenerse mediante combinaciones lineales de otras soluciones. Estas últimas pueden ser ordinarias (EDOs) o en derivadas parciales (EDPs). Las soluciones a las ecuaciones diferenciales lineales cuando son homogéneas forman un espacio vectorial, a diferencia de las ecuaciones diferenciales no lineales.

Definición

La ecuación

$\phi (x^{(n)},...,x',x,t)=0$

se llama lineal cuando la función $\phi$ es lineal a las variables $x^{(k)}$. ​

Introducción

Una ecuación diferencial lineal tiene la forma:

$Ly=f$

donde el operador diferencial L es un operador lineal, y es la función incógnita o desconocida (una función que podría ser dependiente del tiempo y(t)), y del lado derecho f es una función conocida de la misma naturaleza que y (denominada término de excitación). Para una función dependiente del tiempo se puede escribir la ecuación más detalladamente como:

$Ly(t)=f(t)$

y también se puede usar la notación con corchetes:

$L[y(t)]=f(t)$

El operador lineal L puede ser de la siguiente forma:​

$L_{n}(y)\equiv a_{n}(t){\frac {d^{n}y}{dt^{n}}}+a_{n-1}(t){\frac {d^{n-1}y}{dt^{n-1}}}+\cdots +a_{1}(t){\frac {dy}{dt}}+a_{0}(t)y$

o sino:

$L_{n}(y)\equiv \sum _{k=0}^{n}a_{k}(x)D^{k}(y)$

La condición de linealidad sobre L se da mientras no aparezcan productos de la función desconocida consigo misma, ni con ninguna de sus derivadas. Es conveninente reescribir esta ecuación en donde la forma del operador es:

$L_{n}(y)\equiv \left[\,a_{n}(t)D^{n}+a_{n-1}(t)D^{n-1}+\cdots +a_{1}(t)D+a_{0}(t)\right]y$

donde D es el operador diferencial d/dt (es decir, Dy = y' , D²y = y",... ), y a_k son funciones conocidas. Se dice que la ecuación tiene un orden n, si es el índice más alto de la derivada de y.

Estructura del espacio de soluciones

Si $f(t)$ es idénticamente nula la ecuación se denomina homogénea. La solución general de la ecuación homogénea viene dada por todas las combinaciones lineales de tantas soluciones linealmente independientes como el orden de dicha ecuación; es decir, que sus soluciones forman un espacio vectorial de dimensión $n$. Para comprobar la independencia lineal de un conjunto de $n$ soluciones podemos calcular su wronskiano. Si el wronskiano no se anula en el intervalo de definición de las soluciones, estas son linealmente independientes.

Si $f(t)\neq 0$ la ecuación se denomina no homogénea. La solución general de la ecuación no homogénea viene dada por

$\,y(t)=y_{p}(t)+y_{g}(t)$,

donde $y_{g}(t)$ es la solución general de la ecuación homogénea asociada e $y_{p}(t)$ es una solución particular a la ecuación no homogénea. Es decir, que sus soluciones forman un espacio afín de dimensión $n$.

En el caso de la solución particular, existen varios métodos para encontrala, entre ellos, el método de variación de los parámetros.

Si fijamos ciertas condiciones iniciales, tenemos garantizada la existencia y unicidad de solución local por el Teorema de Picard-Lindelöf siempre que las funciones $a_{k}$ y $f(t)$ sean continuas y acotadas en un entorno de los valores iniciales. Si, además, tenemos que $a_{k}$ y $f$ son continuas y acotadas en todo el espacio, tendriamos garantizada la existencia y unicidad global de las soluciones en todo el espacio.

Ecuación lineal de primer orden

Las Ecuaciones diferenciales de primer orden se caracterizan por ser de la forma:

${\begin{cases}y'+p(x)y=q(x)\\y(x_{0})=y_{0}\end{cases}}$

Donde $p(x)$ y $q(x)$ son funciones continuas en un intervalo abierto $(a,b)\subseteq \mathbb {R}$, y el valor inicial es $y(x_{0})=y_{0}$.

La solución de esta ecuación viene dada por:

$y(x)=e^{-\int _{x_{0}}^{x}p(s)ds}\left[y_{0}+\int _{x_{0}}^{x}q(s)e^{\int _{x_{0}}^{s}\!p(t)dt}ds\right]$

Ecuaciones lineales de orden n

Del mismo modo que se ha definido la ecuación diferencial lineal de primer orden podemos definir una ecuación diferencial de orden n como:

$y^{(n)}+A_{1}(x)y^{(n-1)}+\dots +A_{n}(x)y=R(x)$

Donde la derivada mayor que aparece es de orden n-ésimo.

Resolución caso general

Esta ecuación se dice que es lineal si la función incógnita o sus derivadas no están multiplicadas por sí mismas o si tampoco aparecen en forma de funciones compuestas (por ejemplo, sen(y) ). Una ecuación diferencial lineal de orden superior puede atacarse convirtiéndola en un sistema de n ecuaciones diferenciales de primer orden. Para hacer esto se definen las n funciones incógnita adicionales dadas por:

$Y_{0}(x):=y(x),\quad Y_{k}(x):={\frac {d^{k}y}{dx^{k}}}$

Puesto que:

$Y_{k}(x)={\frac {dY_{k-1}}{dx}}\ {\text{con}}\ k\leq n-1,\quad Y_{n}(x):=-A_{1}(x)Y_{n-1}(x)-\dots -A_{n}(x)Y_{0}(x)+R(x)$

El sistema de ecuaciones diferenciales puede escribirse en forma de ecuación matricial como:

${\begin{bmatrix}Y_{0}'\\Y_{1}'\\\dots \\Y_{n}'\end{bmatrix}}={\begin{bmatrix}0&1&0&\dots &0\\0&0&1&\dots &0\\\dots &&&\dots \\-A_{n}&-A_{n-1}&-A_{n-2}&\dots &-A_{1}\end{bmatrix}}{\begin{bmatrix}Y_{0}\\Y_{1}\\\dots \\Y_{n}\end{bmatrix}}$

Resolución con coeficientes constantes

La resolución de ecuaciones y sistemas de ecuaciones diferenciales se simplifica mucho si las ecuaciones son de coeficientes constantes. En el caso de una ecuación de primer orden la búsqueda de un factor integrante nos lleva en la mayoría de los casos a una ecuación en derivadas parciales. Si la ecuación es de orden superior, a no ser que sea una ecuación de Euler o similar, tendremos que proponer una solución que no viene dada, en general, por funciones elementales. En estos casos los métodos preferidos (sin contar el cálculo numérico) son los que emplean series de potencias o series de Fourier. En el caso de los sistemas, si la matriz del sistema es de coeficientes constantes podemos resolver el sistema usando el método de los valores propios ya que en ese caso la matriz resultante de la reducción de la ecuación a un sistema de primer orden es constante y puede encontrarse fácilmente su solución calculando la exponencial de la matriz del sistema.

Para estudiar otros métodos de encontrar la solución aparte de la exponenciación de matrices consideraremos una ecuación del tipo:

$y^{(n)}+a_{1}y^{(n-1)}+\dots +a_{n}y+b=0$

Donde $a_{k}(k=0,1,\dots ,n)\in \mathbb {R}$ son coeficientes constantes conocidos. Observemos que la derivada n-ésima va acompañada por el coeficiente unidad. Definimos el polinomio característico de la ecuación como

$r^{n}+a_{1}r^{n-1}+\dots +a_{n}=0$

que es una ecuación algebraica de orden n. Se demuestra que si hallamos las n raíces $\lambda _{n}$ del polinomio característico la solución de la ecuación homogénea:

$y(x)=y_{1}(x)+\dots +y_{n}(x)$

Al calcular las raíces $\lambda _{n}$ del polinomio característico pueden darse los siguientes casos:

• Raíces reales distintas: En este caso la solución viene dada directamente por $y(x)=y_{1}(x)+\dots +y_{n}(x)$, donde $y_{k}(x)=C_{k}e^{\lambda _{k}x}$, siendo $C_{k}$ constantes de integración.
• Raíces reales repetidas: Ilustraremos este caso con un ejemplo; sea una ecuación de segundo orden con coeficientes constantes cuyo polinomio característico tiene la raíz $\lambda _{i}$ doble. En este caso no podemos expresar la solución como $y(x)=2Ce^{\lambda x}$, ya que si lo hacemos de este modo tenemos una información redundante. En este caso particular la solución de la ecuación es $y(x)=C_{1}e^{\lambda x}+C_{2}xe^{\lambda x}$. En general, en una ecuación de orden n, si una raíz $\lambda _{0}$ aparece repetida q veces la solución parcial asociada a ella es:

$y(x)=\sum _{j=1}^{q}C_{j}x^{j-1}e^{\lambda _{0}x}$

• Raíces complejas: Si las raíces son del tipo $\lambda _{k}=a_{k}+b_{k}i$ debemos expresar la solución como combinación lineal de senos, cosenos y exponenciales en la forma

$y_{k}(x)=e^{a_{k}x}[\cos(b_{k}x)+\sin(b_{k}x)]$

Si las raíces complejas conjugadas están repetidas q veces, la ecuación es del tipo

$y_{k}(x)=e^{a_{k}x}\left[\sum _{j=1}^{q}C_{j}x^{j-1}\cos(b_{k}x)+\sum _{j=1}^{q}C_{j}x^{j-1}\sin(b_{k}x)\right]$

Una vez resuelto el problema homogéneo podemos atacar el problema completo. Para tener la solución del problema completo debemos sumar una solución particular a la solución homogénea ya obtenida:

$y(x)=y_{p}(x)+y_{h}(x)=y_{p}(x)+\sum _{k=1}^{n}C_{k}e^{\lambda _{k}x}$

Para hallar $\!y_{p}(x)$ empleamos el método de la conjetura razonable, consistente en analizar el término inhomogéneo de la ecuación y proponer funciones del mismo tipo como solución. Nótese que no es necesario que $\!b$ sea un coeficiente constante.

Ejemplos

• Tenemos $b=x$. Proponemos $y_{p}(x)=Ax+B$ (polinomio de primer orden). Las constantes $A$ y $B$ quedan determinadas tras aplicar los requerimientos de la ecuación a la solución particular (derivar n veces, multiplicar por $a_{n}$ coeficientes constantes, etc.).
• Tenemos $b=\cos(2x)$. Proponemos $y_{p}(x)=A\cos(2x)+B\sin(2x)$. Las constantes $A$ y $B$ se determinan como en el ejemplo 1.

Véase también

• Ecuaciones diferenciales de primer orden
• Ecuaciones diferenciales ordinarias
• Ecuaciones diferenciales

Referencias

Bibliografía

• Zill, Dennis G. (2006). Ecuaciones diferenciales con aplicaciones (2ª edición). Grupo Editorial Iberoamérica. 
• Aranda Iriarte, José Ignacio (2008). Apuntes de ecuaciones diferenciales I. Universidad Complutense de Madrid. Archivado desde el original el 15 de febrero de 2017. Consultado el 4 de abril de 2016. 
• Birkhoff, Garrett; Rota, Gian-Carlo (1978). Ordinary Differential Equations (en inglés). New York: John Wiley and Sons, Inc. ISBN 0-471-07411-X. 
• Gershenfeld, Neil (1999). The Nature of Mathematical Modeling (en inglés). Cambridge, UK.: Cambridge University Press. ISBN 978-0-521-57095-4. 
• Robinson, James C. (2004). An Introduction to Ordinary Differential Equations (en inglés). Cambridge, UK.: Cambridge University Press. ISBN 0-521-82650-0. 

• «Ordinary Differential Equations». Handbook of Differential Equations: Ordinary Differential Equations. 2008. ISSN 1874-5725. doi:10.1016/s1874-5725(08)x8001-0. Consultado el 8 de mayo de 2022. 

• «Équations Différentielles Ordinaires, Vincent Millot». 

Enlaces externos

• Soluciones exactas de Ecuaciones Diferenciales Ordinarias

• Datos: Q1129902