Ecuación en diferencia lineal

En matemáticas, se define una ecuación en diferencias lineal o relación de recurrencia lineal como una sucesión $\{a_n\}$ definida en función de elementos anteriores de esa misma sucesión para todo n no negativo y sus términos únicamente pueden tener grado 0 y 1 para cumplir que sea lineal.

Se define una ecuación en diferencias lineal de grado k como:

$a_n={\alpha }_1{a}_{n-1}+{\alpha }_2{a}_{n-2}+...+{\alpha }_k{a}_{n-k}+F(n)\mathrm{\forall }{\alpha }_i\in ℝ\wedge {\alpha }_k\ne 0$,

donde $F(n)$ puede ser una constante, un polinomio o una función exponencial en función de n.

Si $F(n)=0\mathrm{\forall }n\in {\mathbb{N}}_{\mathrm{𝟘}}$, es decir, que $F(n)$ sea la constante 0 se tratará de una ecuación lineal en diferencias homogénea.

Sus aplicaciones están relacionadas principalmente con el tratamiento de sistemas dinámicos donde las relaciones de recurrencia puedan moldear de manera adecuada este tipo de sistema.

Resolución

Debido a que el tratamiento de las relaciones de recurrencias puede ser costoso computacionalmente, es adecuado expresar las ecuaciones en diferencias lineales únicamente en función de n. Para ello existe un método sistemático para "resolverlas", tanto para las homogéneas y no homogéneas. Además, se deben tener k condiciones iniciales para poder dar una respuesta única a la ecuación. Las condiciones iniciales serían elementos $a_i$ de la sucesión de los que ya conocemos su valor. En la mayoría de casos los calcularemos para lo valores más sencillos si no los conocemos, es decir, $i=0,1,..$ .

Ecuación en diferencias lineal homogénea

En primer lugar hay que obtener el polinomio característico de la relación de recurrencia, que tendría la siguiente forma:

$r^k-{\alpha }_1{r}^{k-1}-{\alpha }_2{r}^{k-2}-...-{\alpha }_k=0$

Al ser un polinomio de grado k por el teorema fundamental del álgebra tendrá k raíces complejas. En este artículo solo se tratará siendo todas las raíces reales. Además, tendríamos dos casos dependiendo de si hay multiplicidad de raíces, es decir, que la cualquier raíz se repita dos o más veces.

Raíces únicas

Obtenemos el conjunto de soluciones $\{r_1,r_2,...,r_k\}$ y podemos expresar $a_n$ de la siguiente manera:

$a_n={\delta }_1{r_1}^n+{\delta }_2{r_2}^n+...+{\delta }_k{r_k}^n$ donde ${\delta }_i$ son contantes desconocidas cuyo valor hallaremos usando las condiciones iniciales.

Las condiciones iniciales serían: $\{a_{{\gamma }_1}={\beta }_1,a_{{\gamma }_2}={\beta }_2,...,a_{{\gamma }_k}={\beta }_k\}{\gamma }_i\in {\mathbb{N}}_{\mathrm{𝟘}},{\beta }_i\in ℝ$

Planteamos un sistema de ecuaciones lineal con las condiciones iniciales:

$\{\begin{array}{c}{\beta }_1={\delta }_1{r_1}^{{\gamma }_1}+{\delta }_2{r_2}^{{\gamma }_1}+...+{\delta }_k{r_k}^{{\gamma }_1}\\ {\beta }_2={\delta }_1{r_1}^{{\gamma }_2}+{\delta }_2{r_2}^{{\gamma }_2}+...+{\delta }_k{r_k}^{{\gamma }_2}\\ ....\\ ....\\ {\beta }_k={\delta }_1{r_1}^{{\gamma }_k}+{\delta }_2{r_2}^{{\gamma }_k}+...+{\delta }_k{r_k}^{{\gamma }_k}\end{array}$

Una vez ya obtenidas las soluciones las sustituimos por las constantes y obtenemos el resultado final.

Ejemplo

Usaremos este procedimiento para resolver la sucesión Fibonacci:

La sucesión Fibonacci se define cómo $a_n=a_{n-1}+a_{n-2}$ y sus condiciones iniciales son $\{a_0=0,a_1=1\}$.

Obtenemos el polinomio característico:

$r^2-r-1=0$

Las soluciones son $\{\varphi ,1-\varphi \}$ donde $\varphi$ es la razón áurea. Y nos queda:

$a_n={\delta }_1{\varphi }^n+{\delta }_2(1-\varphi {)}^n$

Resolvemos y el resultado es:

${\delta }_1=1/\sqrt{5}$ ${\delta }_2=-1/\sqrt{5}$

$a_n=\frac{{\varphi }^n-(1-\varphi {)}^n}{\sqrt{5}}$

Raíces múltiples

En el caso de que el polinomio característico tenga raíces múltiples no se puede formalizar la ecuación de la misma manera, la formalización necesaria sería:

Sean $\{r_1,r_2,..,r_t\}$ las raíces del polinomio característico y con multiplicidades $\{m_1,m_2,..,m_t\}/m_1+m_2+...+m_t=k$.

Tenemos:

$a_n=({\alpha }_{1,1}+{\alpha }_{1,2}n+...+{\alpha }_{1,m_1}{n}^{m_1-1})r_1^n+({\alpha }_{2,1}+{\alpha }_{2,2}n+...+{\alpha }_{2,m_2}{n}^{m_2-1})r_2^n+...+({\alpha }_{t,1}+{\alpha }_{t,2}n+...+{\alpha }_{t,m_t}{n}^{m_t-1})r_t^n$

para $\mathrm{\forall }n\in {\mathbb{N}}_{\mathrm{𝟘}}$, $\mathrm{\forall }{\alpha }_{i,j}\in ℝ,1\le i\le t,1\le j\le m_i$

De la misma manera plantearemos un sistema de ecuaciones lineales para resolverlo usando las condiciones iniciales de la ecuación en diferencias lineal.

Ejemplo

Tenemos la recurrencia $a_n=4a_{n-1}-4a_{n-2}$ y condiciones iniciales $\{a_0=8,a_1=20\}$.

El polinomio característico sería:

$r^2-4r+4=0$

Resolvemos y obtenemos 2 con multiplicidad de 2. Entonces, la solución sería de la forma:

$a_n=({\alpha }_1+{\alpha }_{2}n)2^n$

Resolvemos el sistema de ecuaciones lineales:

$\{\begin{array}{c}8={\alpha }_1\\ 20=2({\alpha }_1+{\alpha }_2)\end{array}$

Las soluciones de ese sistema de ecuaciones son $\{{\alpha }_1=8,{\alpha }_2=2\}$ y la solución a la recurrencia es $a_n=(8+2n)2^n$.

Ecuación en diferencias lineal no homogénea

Cuando $F(n)$ es una constante distinta de 0, $a_n$ se formula como $a_n^{(p)}+a_n^{(h)}$ donde $a_n^{(p)}$ es una solución de la relación de recurrencia lineal original y $a_n^{(h)}$una solución a la relación de recurrencia original sin el término $F(n)$. No existe un esquema general para obtener una solución de $a_n^{(p)}$, pero en el caso de que $F(n)$ sea un polinomio de grado k puede ser expresado por $\sum _{i=0}^k{\beta }_i{n}^i$ donde ${\beta }_i$ son constantes a determinar y en el caso de ser una función exponencial ${\beta }^n$como $\delta {\beta }^n$ donde $\delta$ es también una constante a determinar. Una vez resuelto $a_n^{(p)}$ se resuelve $a_n^{(h)}$ con su polinomio característico y se suman, después procedemos a determinar las constantes de $a_n^{(h)}$ para terminar de resolverlo usando las condiciones iniciales y obtendríamos el resultado final.

Ejemplo

Tenemos la recurrencia $a_n=5a_{n-1}+5n$ y condición inicial $\{a_1=4\}$.

Resolvemos $a_n^{(h)}$ = $5a_{n-1}^{(h)}$ y nos queda $a_n^{(h)}$ = $\alpha 5^n$ donde $\alpha$ es una constante a determinar.

Para $a_n^{(p)}$ al ser $F(n)$ un polinomio de grado 1 la solución sería $a_n^{(p)}=\gamma n+\delta$ donde $\gamma$ y $\delta$ son constantes a determinar.

Sustituimos en $a_n$:

$\gamma n+\delta =5(\gamma (n-1)+\delta )+5n$

Resolvemos y nos queda $\gamma =-5/4$ y $\delta =-25/16$ y por tanto $a_n$ es:

$a_n=a_n^{(p)}+a_n^{(h)}=-\frac{5}{4}n-\frac{25}{16}+\alpha 5^n$

Solo nos queda determinar $\alpha$ usando la condición inicial:

$4=-\frac{5}{4}-\frac{25}{16}+5\alpha ,\alpha =\frac{109}{80}$

Y la solución sería:

$a_n=-\frac{5}{4}n-\frac{25}{16}+\frac{109}{80}{5}^n$