Ecuación diferencial ordinaria

En matemáticas, una ecuación diferencial ordinaria (comúnmente abreviada "EDO") es un tipo de ecuación diferencial caracterizada porque la variable dependiente está en función de una variable independiente; este hecho las distingue de las Ecuaciones en derivadas parciales. Estas ecuaciones son importantes en diversas áreas de estudio como la geometría, mecánica y astronomía, además de muchas otras aplicaciones.

Se ha dedicado mucho estudio a la resolución de este tipo de ecuaciones, estando casi completamente desarrollada la teoría para ecuaciones lineales. Sin embargo la mayoría de las ecuaciones diferenciales interesantes son no-lineales, a las cuales en la mayoría de los casos no se les puede encontrar una solución exacta.

Introducción

Si F es esta relación o función, la ecuación diferencial ordinaria (EDO) es

$\ F(x,y,y',y'',\dots ,y^{(n)})=0$

La ecuación diferencial lineal más general, de orden n está dada por:

$\ a_{n}(t)y^{(n)}+a_{n-1}(t)y^{(n-1)}+\ldots +a_{1}(t)y'+a_{0}(t)y=g(t)$

Donde los $a_{i}$ representan funciones dependientes de t.

Una solución tipo de esta ecuación será una "familia" de curvas o funciones del tipo $x=x(t)$ que verifica la ecuación.

Tipos de EDOs y forma de resolución

Existen diversos tipos de ecuaciones diferenciales ordinarias, cada una con una forma de resolución distinta; para clasificarlas, hay que hacer la diferencia entre ecuaciones diferenciales de primer orden y ecuaciones de orden superior (ya que las primeras son, por lo general, de más fácil resolución).

Existencia y unicidad de soluciones

El teorema de Peano-Picard garantiza la existencia de una solución y su unicidad para toda ecuación diferencial ordinaria lineal con coeficientes continuos en un intervalo tiene solución única en dicho intervalo. Para el caso de ecuaciones diferenciales no-lineales no existen resultados análogos al de Peano-Picard.

Soluciones analíticas

Existen métodos de resolución generales para ecuaciones diferenciales ordinarias lineales que permiten encontrar soluciones analíticas. En particular si los coeficientes de la ecuación lineal son constantes o periódicos la solución es casi siempre fácil de construir. Para coeficientes no constantes o no periódicos, pero que son desarrollables en serie de Taylor o serie de Laurent es aplicable con ciertas restricciones el método de Frobenius. Otra posibilidad es reducir una ecuación diferencial lineal de orden n a un sistema de n ecuaciones diferenciales lineales de primer orden.

Para las ecuaciones diferenciales ordinarias no-lineales no existen métodos generales.

Soluciones numéricas

Algunos de los métodos de solución numérica de ecuaciones diferenciales son el método de Runge-Kutta, los métodos multipaso y los métodos de extrapolación.

Ecuaciones diferenciales ordinarias de primer orden

Una ecuación diferencial de primer orden con la condición inicial se expresa de la siguiente forma:

$[L]=\left\{{\begin{array}{*{20}c}{{\cfrac {dy}{dt}}=f(t,y)}\\{y(t_{0})=y_{0}}\\\end{array}}\right.$

Donde $y(t_{0})=y_{0}\,$ es la condición inicial.

Entre los tipos de EDOs de primer orden se encuentran:​

Ecuación de variables separables

Son EDOs de la forma:

${\frac {dy}{dt}}=f(t,y)$

En donde es posible "despejar" todos los términos con la variable dependiente en función de la variable independiente, quedando ahora la ecuación:

$\ g(y)dy=h(t)dt$

En donde se procede integrando ambos miembros de la ecuación

$\int g(y)dy=\int h(t)dt$

De donde es posible obtener la solución

Ecuación exacta

Véase también: Ecuación diferencial exacta

Una ecuación de la forma:

$M(x,y)\,dx+N(x,y)\,dy=0,\,\!$

se dice exacta si existe una función F que cumpla:

${\frac {\partial F}{\partial x}}(x,y)=M$

y

${\frac {\partial F}{\partial y}}(x,y)=N.$

Su solución es entonces:

$F(x,y)=C.\,$

• Ecuación de Coeficientes Homogéneos (llamada comúnmente homogénea).

Ecuación lineal

Véase también: Ecuación diferencial lineal

Una ecuación diferencial es lineal si presenta la forma:

$y'+P(x)y=Q(x)\,$

Y que tienen por solución:

$y(x)=e^{-\int P(x)dx}\left(C+\int Q(x)e^{\int P(x)dx}dx\right)$

Como se puede apreciar, esta ecuación es una ecuación diferencial de Bernoulli, con n=0.

Ecuación de Bernoulli

Véase también: Ecuación diferencial de Bernoulli

Una ecuación diferencial de Bernoulli, que es a su vez una generalización de la ecuación diferencial lineal, fue formulada por Jakob Bernoulli y resuelta por su hermano, Johann Bernoulli y presenta la forma:

$y'+P(x)y=Q(x)y^{n}\,$

En la cual, si se hace la sustitución $z=y^{1-n}$, la ecuación se transforma en una ecuación lineal con z como variable dependiente, resolviéndose de manera análoga.

Ecuación de Riccati

Véase también: Ecuación diferencial de Riccati

Una ecuación diferencial tiene la forma de la introducida por Jacobo Francesco Riccati cuando presenta la estructura:

$y'(x)+P(x)y^{2}+Q(x)y+R(x)=0\,$

Para resolverla, se debe hacer la sustitución $y=y_{p}+{\frac {1}{z}}$, donde $y_{p}$ es una solución particular cualquiera de la ecuación.

Ecuación de Lagrange

Véase también: Ecuación diferencial de Lagrange

Una ecuación diferencial de Lagrange presenta la forma:

$y=g(y')x+f(y')\,$

Resolviéndose con la sustitución $y'=p$, obteniéndose una solución general y una solución particular.

Ecuación de Clairaut

Véase también: Ecuación diferencial de Clairaut

Una ecuación diferencial de Clairaut, llamada así en honor a Alexis-Claude Clairaut,tiene la forma:

$y=xy'+f(y')\,$

Como se puede apreciar, esta ecuación es una forma particular de la ecuación diferencial de Lagrange, con $g(y')=y'$, por lo cual, su resolución es análoga a la anterior.

Ecuaciones diferenciales ordinarias de segundo orden

Muchos problemas físicos importantes tanto en mecánica como en electromagnetismo conllevan la resolución de ecuaciones diferenciales de segundo orden.

Ecuación lineal con coeficientes constantes

La ecuación diferencial de segundo orden con coeficientes constantes tiene la forma:

$a{\frac {d^{2}y}{dx^{2}}}+b{\frac {dy}{dx}}+cy=f(x)$

La resolución de esta ecuación depende de las raíces del polinomio característico:

$a\lambda ^{2}+b\lambda +c=0\,$

En función de como sean las raíces de dicho polinomio se distinguen tres casos posibles:

• Caso 1: dos raíces reales y distintas $(\lambda _{1}\neq \lambda _{2})\,$, en este caso la solución general tiene la forma:

$y(x)=C_{1}e^{\lambda _{1}x}+C_{2}e^{\lambda _{2}x}+{\frac {e^{\lambda _{1}x}}{\lambda _{1}-\lambda _{2}}}\int _{x_{0}}^{x}e^{-\lambda _{1}u}f(u)du+{\frac {e^{\lambda _{2}x}}{\lambda _{2}-\lambda _{1}}}\int _{x_{0}}^{x}e^{-\lambda _{2}u}f(u)du$

• Caso 2: dos raíces reales e iguales $(\lambda _{1}=\lambda _{2})\,$, en este caso la solución general tiene la forma:

$y(x)=C_{1}e^{\lambda _{1}x}+C_{2}xe^{\lambda _{1}x}+xe^{\lambda _{1}x}\int _{x_{0}}^{x}e^{-\lambda _{1}u}f(u)du-e^{\lambda _{1}x}\int _{x_{0}}^{x}xe^{-\lambda _{2}u}f(u)du$

• Caso 3: dos raíces complejas conjugadas $(\lambda _{1}=p+qi,\lambda _{2}=p-qi)\,$, en este caso la solución general tiene la forma:

$y(x)=e^{px}(C_{1}\cos qx+C_{2}\sin qx)+{\frac {e^{px}\sin qx}{q}}\int _{x_{0}}^{x}e^{-pu}f(u)\cos qu\ du-{\frac {e^{px}\cos qx}{q}}\int _{x_{0}}^{x}xe^{-pu}f(u)\sin qu\ du$

El último término de esta última ecuación está relacionado con la integral de Duhamel.

Ecuación diferencial de Euler o de Cauchy

Esta ecuación tiene la forma:

$x^{2}{\frac {d^{2}y}{dx^{2}}}+ax{\frac {dy}{dx}}+by=g(x)$

Y puede resolverse mediante haciendo el cambio de variable $x=e^{t}\,$ que reduce la ecuación anterior a una ecuación de coeficientes constantes resoluble por los métodos de la sección anterior:

${\frac {d^{2}{\bar {y}}}{dt^{2}}}+(a-1){\frac {d{\bar {y}}}{dx}}+b{\bar {y}}=g(e^{t}),\qquad {\bar {y}}(t)=y(e^{t})$

Ecuaciones de Bessel

La ecuación diferencial de Bessel, aparece con frecuencia en la resolución del problema de Dirichlet en coordenadas cilíndricas. Dicha ecuación tiene la forma:

$x^{2}{\frac {d^{2}y}{dx^{2}}}+x{\frac {dy}{dx}}+(x^{2}-n^{2})y=0$

Esta ecuación es resoluble mediante las llamadas funciones de Bessel:

$y(x)=C_{1}J_{n}(x)+C_{2}Y_{n}(x)\,$

Además de esta ecuación existe otra ecuación resoluble mediante funciones de Bessel. La ecuación diferencial de Bessel modificada, aparece con frecuencia en la resolución del problema de Dirichlet en coordenadas cilíndricas. Dicha ecuación tiene la forma:

$x^{2}{\frac {d^{2}y}{dx^{2}}}+(2p+1)x{\frac {dy}{dx}}+(\alpha x^{2r}+\beta ^{2})y=0$

Cuya solución viene dada por:

$y(x)=x^{-p}\left[C_{1}J_{q/r}\left({\frac {\alpha }{r}}x^{r}\right)+C_{2}Y_{q/r}\left({\frac {\alpha }{r}}x^{r}\right)\right],\qquad q:={\sqrt {p^{2}-\beta ^{2}}}$

Ecuación de Legendre

La ecuación diferencial de Legendre, aparece con frecuencia en la resolución del problema de Dirichlet en coordenadas esféricas. La ecuación tiene la forma:

$(1-x^{2}){\frac {d^{2}y}{dx^{2}}}-2x{\frac {dy}{dx}}+n(n+1)y=0$

Cuando n es un entero una de las dos soluciones independientes que conforman la solución general de la ecuación anterior es el polinomio de Legendre de grado n:

$P_{n}(x)={\frac {1}{2^{n}n!}}\ {\frac {d^{n}}{dx^{n}}}(x^{2}-1)^{n}$

Las solución general puede expresarse en la forma:

$y(x)=C_{1}U_{n}(x)+C_{2}V_{n}(x)\,$, o bien, $=y(x)={\bar {C}}_{1}P_{n}(x)+{\bar {C}}_{2}Q_{n}(x)$

Donde:

${\begin{cases}U_{n}(x)=1-{\cfrac {n(n+1)}{2!}}x^{2}+{\cfrac {n(n-1)(n+1)(n+3}{4!}}x^{4}-\ldots \\V_{n}(x)=x-{\cfrac {(n-1)(n+2)}{3!}}x^{3}+{\cfrac {(n-1)(n-3)(n+2)(n+4}{5!}}x^{5}-\ldots \end{cases}}$

$P_{n}(x)={\begin{cases}U_{n}(x)/U_{n}(1)&n=0,2,4,\ldots \\V_{n}(x)/V_{n}(1)&n=1,3,5,\ldots \end{cases}}$, y $Q_{n}(x)={\begin{cases}V_{n}(x)U_{n}(1)&n=0,2,4,\ldots \\-U_{n}(x)V_{n}(1)&n=1,3,5,\ldots \end{cases}}$

Ecuaciones diferenciales ordinarias de orden superior

Ecuación lineal de orden n con coeficientes constantes

La ecuación diferencial lineal de orden n con coeficientes constantes es dela siguiente forma:

$\ a_{n}y^{(n)}+a_{n-1}y^{(n-1)}+\ldots +a_{1}y'+a_{0}y=g(t)$

Donde los términos $a_{i}\,$ representan constantes $\forall i\in \mathbb {N}$ En el caso homogéneo cuando el segundo miembro es idénticamente nulo, las soluciones de esta ecuación se pueden obtener a partir de la raíces del polinomio característico de la ecuación:

$\ a_{n}\lambda ^{n}+a_{n-1}\lambda ^{n-1}+\ldots +a_{1}\lambda +a_{0}=0$

En el caso de que todas las raíces sean diferentes la solución viene dada por:

$y(x)=C_{1}e^{\lambda _{1}x}+\ldots +C_{n}e^{\lambda _{n}x}=\sum _{i=1}^{n}C_{i}e^{\lambda _{i}x}$

En el caso de que existan varias raíces repetidas, siendo $m_{i}$ la multiplicidad de la raíz i-ésima, la solución es de la forma:

$y(x)=\sum _{i=1}^{k}\left(C_{i,0}+C_{i,1}x+\ldots +C_{i,m_{i}-1}x^{m_{i}-1}\right)e^{\lambda _{i}x}=\sum _{i=1}^{k}\left(\sum _{j=0}^{m_{i}-1}C_{i,j}x^{j}\right)e^{\lambda _{i}x},\quad k\leq n,\ \sum _{j=1}^{k}m_{j}=n$

Las multiplicidades de cada raíz son el exponente de la siguiente descomposición:

$a_{n}(\lambda -\lambda _{1})^{m_{1}}\ldots (\lambda -\lambda _{k})^{m_{k}}=a_{n}\lambda ^{n}+a_{n-1}\lambda ^{n-1}+\ldots +a_{1}\lambda +a_{0}=0$

Bibliografía

• José Ignacio Aranda Iriarte (2007). Apuntes de ecuaciones diferenciales I. Universidad Complutense de Madrid.

Véase también

• Ecuación diferencial
• Ecuaciones en derivadas parciales
• Teorema de Picard

Enlaces externos

• Soluciones exactas de Ecuaciones Diferenciales Ordinarias
• Métodos de resolución de Ecuaciones Diferenciales Ordinarias