Ecuaciones lineales a coeficientes constantes de orden n

Buscaremos soluciones de un sistema lineal a coeficientes constantes de orden n no homogéneo. Aún así, tendremos que añadir una restricción en el método que presentaremos.

Si nuestra EDO es: any(n)(x)+an1y(n1)(x)++a1y(x)+a0y(x)=f(x) (lineal a coeficientes constantes) tenemos que pedir que la función f(x) sea un polinomio, una exponencial, seno o coseno o cualquier combinación de éstas.

Es decir, ahora estaremos preparados para resolver por ejemplo: y+y=3cosx+e2x

Por el mismo motivo que en sistemas lineales, una solución general de esta ecuación es una suma de la solución general de la parte homogénea y una solución particular de la no homogénea.

Vamos a resolver la EDO por el método del polinomio anulador o de coeficientes indeterminados.

Supongamos que tenemos la EDO escrita anteriormente y f(x) una función que cumpla las condiciones que hemos pedido. Entonces:

  • Resolvemos la parte homogénea. De forma que obtenemos n soluciones linealmente independientes.

Ejemplo

En el ejemplo que hemos dado, la soluciones son:y1(x)=cosxy2(x)=sinx

  • Buscamos un polinomio que anule f(x). Esta operación consiste a encontrar un polinomio que sus coeficientes multipliquen las derivadas. Es decir:Q(D)=bnDn+bn1Dn1++b1D+b0Idwhere Dk significa derivar k veces la función que lo multiplica. Así,Q(D)f(x)=(bnDn+bn1Dn1++b1D+b0Id)f(x)==bnDnf(x)+bn1Dn1f(x)++b1Df(x)+b0Idf(x)==bnfn(x)+bn1fn1(x)++b1f(x)+b0f(x)=0Es decir, es como buscar qué EDO lineal y homogénea satisface f(x). Para hacerlo procedemos de forma inversa (de cuando encontramos soluciones en el caso homogéneo).

Ejemplo

En el ejemplo anterior, tenemos que encontrar un polinomio que anule f(x)=3cosx+e2x.

Procedamos de forma inversa que cuando encontramos soluciones, es decir: cosx proviene de λ=ie2x que proviene de λ=2.

Por lo tanto el polinomio anulador es: Q(D)=(D2+ID)(D2Id).

En efecto, Q(D)f(x)=(D2+Id)(D2Id)f(x)=(D32D2+D2Id)f(x)==f(x)ef(x)+f(x)2f(x)==3sinx+8e2x+6cosx8e2x3sinx+2e2x6cosx2e2x=0

  • Notemos que, introduciendo esta notación, nuestra EDO inicial se puede escribir como P(D)y(x)=f(x) , con P(D)=anDn+an1Dn1++a1D+a0Id Aplicando el polinomio Q(D) a la anterior igualdad tenemos: P(D)y(x)=f(x)Q(D)P(D)y(x)=Q(D)f(x)=0 y por lo tanto tenemos una nueva ecuación, pero homogénea de orden k (mayor que n). Entonces solucionamos este problema, obteniendo k funciones, de las cuales las n primeras son soluciones encontradas en (1).y(x)=C1y1(x)+C2y2(x)++Cnyn(x)+D1y~n+1(x)++Dky~k(x)

Ejemplo

En el ejemplo anterior, pues, tenemos Q(D)P(D)=(D2+Id)(D2+Id)(D2Id) que tiene por raíces (y por lo tanto por soluciones asociadas): λ=±i con multiplicidad 2 que da por soluciones cosx,sinx,xcosx,xsinx, λ=2 que da por solución e2x.

Por lo tanto tenemos que y(x)=C1cosx+C2sinx+D1xcosx+D2xsinx+D3e2x

  • Buscamos una solución particular de la ecuación no homogénea de la forma: yp(x)=D1y~n+1(x)++Dky~k(x) es decir, cogemos las soluciones que han aparecido en (3), que no teníamos en (1) y buscamos ciertos coeficientes para obtener la solución.

Ejemplo

En nuestro ejemplo, debemos buscar una solución particular de la forma: yp(x)=D1cosx.

Impongamos que sea solución: yp+yp=3cosx+e2xyp+yp=2D1sinxD1xcosx+2D2cosxD2xsinx+4D3e2x+D1xcosx+ +D2xsinx+D3e2x==2D1sinx+2D2cosx+5D3e2x Igualando coeficientes, obtenemos: D1=0D2=32D3=15 Por lo tanto, la solución particular es: yp(x)=32xsinx+15e2x

  • Finalmente, tenemos que la solución general de nuestra EDO no homogénea inicial es:y(x)=yh(x)+yp(x)

Ejemplo

Para acabar con nuestro ejemplo, tenemos que la solución general es: y(x)=c1cosx+c2sinx+32xsinx+15e2x