ECUACIONES EN DIFERENCIAS FINITAS

Luis Puello Garca

Jose

UNIVERSIDAD INDUSTRIAL DE SANTANDER

FACULTAD DE CIENCIAS

ESCUELA DE MATEMATICAS
BUCARAMANGA
2008

ECUACIONES EN DIFERENCIAS FINITAS

Luis Puello Garca

Jose
Monografa presentada como requisito para optar al
ttulo de Licenciado en Matematicas

Director

n Barbosa
Gildardo Guzma
Magister en Matematicas

UNIVERSIDAD INDUSTRIAL DE SANTANDER

FACULTAD DE CIENCIAS

ESCUELA DE MATEMATICAS
BUCARAMANGA
2008

Dedicatoria

La vida esta llena de peque

nas cosas que nos van llevando a ser felices, alcanzar
las metas es una de ellas y luchar por obtenerlas es el mejor de los esfuerzos

Agradecimientos

Agradezco muy especialmente a: Dios por alcanzar esta meta en mi vida, a mi madre
Elba Garca, a mi ta Cleotilde Meja y a mis Hermanos Robinson, Lilian y Ernesto que
me apoyaron y confiaron en mi, al profesor Gildardo Guzman por su colaboracion para la
realizacion de esta monografa, a mis amigos y compa
neros de carrera Francisco Gutierrez,

oscary Avila,
Oscar Madiedo, Arturo Castro, Jaiver Hernandez y Felix Paez por su compa
na
y motivacion, y a mi novia Liliana Ramos por su colaboracion en los momentos difciles.

TITLE: EQUATIONS IN FINITES DIFFERENCES*

AUTHOR: JOSE LUIS PUELLO GARCIA**
KEY WORDS: DIFFERENCE OPERATOR, INDEFINITE SUMMATION, OPERATOR, DIFFERENCE EQUATIONS, HOMOGENEOUS EQUATION,GENERAL SOLUTION, GENERAL PROPERTIES.
DESCRIPTION:
In this essay a bibliographic revision it is made, about difference equations in special the
linears. In the first chapter contains the study of the main operators involved in the calculus and their essential properties. They are also studied also some excellent functions that
allow to simplify the notation, as well as to obtain formulates important for the calculus in
differences, and the indefinite summation, as the inverse process to the calculus in differences. In the second chapter a special development it is made (properties and methods of
analytic solution and graph) to the lineal equations in finite difference of first order with
some applications to the economy and likeness with the continuous calculus. It is also considered some equations that do not fulfill the definition of lineal equations of first order,
but with a variable change they can be taken to equation of that type or as some equations
with variable coefficients, to reduce it or to convert it to constant coefficients. The Third
chapter a similar study it is made (chapter two) for lineal equations in differences of more
order where the properties and the methods are generalized for any order. In each chapter
illustrative examples of the properties are given that are presented.

*
**

Monograph
BARBOSA.
Faculty of sciences. Mathematics school. Director: GILDARDO GUZMAN

TITULO: ECUACIONES EN DIFERENCIAS FINITAS*

LUIS PUELLO GARCIA***
AUTOR: JOSE
PALABRAS CLAVES: Operador diferencia, suma indefinida, operador de traslacion,
ecuaciones en diferencias, ecuacion homogenea, solucion general, propiedades generales.
RESUMEN
En el presente trabajo se hace una revision bibliografica sobre los aspectos generales de
las ecuaciones en diferencias finitas en especial las lineales. En el primer captulo contiene el
estudio de los principales operadores involucrados en el calculo y sus propiedades esenciales.
Se estudiaran asimismo algunas funciones relevantes que permiten simplificar la notacion,
as como obtener formulas importantes para el calculo en diferencias, y la sumas indefinidas,
como el proceso inverso al calculo en diferencias. En el segundo captulo se hace un desarrollo especial (propiedades y metodos de solucion analtica y grafica) a las ecuaciones lineales
en diferencia finitas de primer orden con algunas aplicaciones a la economa y semejanzas
con el calculo continuo. Tambien se considera algunas ecuaciones que no cumplen con la
definicion de ecuaciones lineales de primer orden, pero con un cambio de variable se pueden
llevar a ecuacion de ese tipo o como algunas ecuaciones con coeficientes variables, reducirla
o convertirla a coeficientes constantes. El Tercer y u
ltimo captulo se hace un estudio similar
(capitulo dos) para ecuaciones lineales en diferencias de orden mayor donde se generalizan
las propiedades y los metodos para cualquier orden. En cada uno de los captulos se realizan
ejemplos ilustrativos a cada una de las propiedades que se presentan.

*
***

Monografa
BARBOSA.
Facultad de Ciencias. Escuela de matematicas. Director: GILDARDO GUZMAN

Indice general

INTRODUCCION

1. DIFERENCIAS FINITAS
1.1. Primera diferencia de una funcion .
1.2. Diferencia de Orden Superior . . . .
1.3. Operador de Traslacion E . . . . . .
1.4. La Funcion Factorial . . . . . . . . .
1.5. Igualdad de los Operadores . . . . . .
1.6. Suma Indefinida: El Operador 1 .
1.6.1. Aplicacion al Calculo Integral

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

3
3
5
9
11
15
20
27

.
.
.
.
.
.
.
.

.
.
.
.
.
.
.
.

28
28
31
33
36
40
43
50
52

SUPERIOR.
. . . . . . . . .
. . . . . . . . .
. . . . . . . . .
. . . . . . . . .
. . . . . . . . .
. . . . . . . . .
. . . . . . . . .

.
.
.
.
.
.
.

54
54
55
56
59
64
64
67

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

.
.
.
.
.
.
.

LINEAL EN DIFERENCIAS DE PRIMER ORDEN.

2. ECUACION
2.1. Introduccion . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.2. Solucion de una Ecuacion Lineal en Diferencia de Primer Orden . . . .
2.3. Propiedades de la Ecuacion Homogenea . . . . . . . . . . . . . . . . . .
2.4. Solucion General de una Ecuacion no Homogenea. . . . . . . . . . . . .
2.5. Ecuacion Lineal en Diferencia con Coeficientes Constantes . . . . . . .
2.6. El Metodo de Coeficientes Indeterminados . . . . . . . . . . . . . . . .
2.7. Solucion Geometrica . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.8. Solucion de Ecuaciones no Lineales . . . . . . . . . . . . . . . . . . . .
LINEAL EN DIFERENCIAS DE ORDEN
3. ECUACION
3.1. Introduccion . . . . . . . . . . . . . . . . . . . . . . . . . .
3.2. Existencia de una Solucion . . . . . . . . . . . . . . . . . .
3.3. Propiedades de la Ecuacion Homogenea . . . . . . . . . . .
3.4. Ecuaciones Homogeneas con Coeficientes Constantes . . .
3.5. Solucion General de la Ecuacion Lineal no Homogenea . .
3.5.1. Variacion de Parametros . . . . . . . . . . . . . . .
3.5.2. Metodo de coeficiente indeterminados . . . . . . . .
i

.
.
.
.
.
.
.
.
.
.
.
.
.
.
.

.
.
.
.
.
.
.

BIBLIOGRAFIA.

69

ii


INTRODUCCION

Las ecuaciones en DIFERENCIAS FINITAS juegan un papel muy importante en el

estudio de problemas en los cuales algunas variables asumen solo valores discretos, y de
aqu su influencia en el desarrollo de los metodos numericos relacionados con diferentes
areas como la de las ecuaciones diferenciales, la probabilidad y en general problemas donde
existen relaciones secuenciales entre variables. Se constituyen entonces en una herramienta
adecuada para discretizar problemas en variable continua y hoy en da mucho mas con la
abundancia de computadoras cada vez mas potentes.
El calculo en diferencias se ocupa a grosso modo del estudio de propiedades de las funciones reales de variable real a partir de las diferencias sucesivas entre los valores que toman
dichas funciones en un conjunto determinado de puntos del dominio. Antes del nacimiento
del calculo infinitesimal las funciones que se manejaban tenan como dominio un conjunto
numerable de puntos y el estudio de sus propiedades solo se entenda a partir del calculo
en diferencias finitas. El origen del calculo en diferencias se le puede atribuir a B. Taylor,
aunque el primero que resolvio cuestiones avanzadas fue J. Stirling. No obstante, el primer
tratado sobre el calculo en diferencias aparece en la obra de L. Euler, que introduce por
primera vez el smbolo para expresar las diferencias.
Dedicaremos esta monografa al estudio de los aspectos generales del ecuaciones lineales en
diferencias. Lo dividiremos en tres captulos, en el primero abordaremos el estudio de los
principales operadores involucrados en el calculo y sus propiedades esenciales. Estudiaremos asimismo algunas funciones relevantes que permiten simplificar la notacion, as como
obtener formulas importantes para el calculo en diferencias, y la sumas indefinida como
el proceso inverso al calculo en diferencias, todos estos aspectos se desarrollaran tomando
como referencias los libros [3] y [4]. De este ultimo libro se tomara un ejemplo que haga
analoga entre el operador suma indefinida y la antiderivada como proceso inverso del operador diferencial. En el segundo y tercer capitulo se proporcionara un desarrollo bastante
completo (propiedades y metodos de solucion) de las ecuaciones lineales en diferencia finitas
y algunas aplicaciones hacia la economa, para el desarrollo de estos temarios se tomaran
1

como referencia los textos [1] y [2].

CAPITULO 1

DIFERENCIAS FINITAS

En este captulo se presenta una vision general de las diferencias finitas y algunas de sus
propiedades, que aparecen en situaciones en que interesa estudiar los cambios en los valores
de una funcion, al efectuar variaciones en el dominio de esta. Ademas estas diferencias
aparecen naturalmente en el estudio de las ecuaciones en diferencias finitas, objeto de la
presente monografa.
Ejemplos de diferencias finitas se dan en diversos contextos, por ejemplo, cuando se estudian
sucesiones numericas. Al considerar una sucesion (xn ) cuyos primeros terminos vienen dados
por 0,2,6,12,20,30,42,56 ,. . . y asumiendo que la construccion de sus terminos sigue una
regularidad que puede intuirse de los terminos dados, puede observarse que x2 x1 = 2,
x3 x2 = 4, x4 x3 = 6, x5 x4 = 8, x6 x5 = 10, x7 x6 = 12, x8 x7 = 14. Siguiendo
este patron una sucesion como la indicada tiene la propiedad que xn+1 xn = 2n, para
todo n
umero entero mayor o igual que uno. El lado izquierdo de la u
ltima igualdad, es una
diferencia finita de primer orden.

1.1.

Primera diferencia de una funci

on

Definici
on 1.1.1. Dadas una funcion f : S R, con S R y una constante h tales
que siempre que x pertenezca al dominio de f , x + h tambien; se define la funcion primera
diferencia de f como la funcion f tal que
(1.1)

f (x) = f (x + h) f (x).

Ejemplo 1.1.1. La primera diferencia de la funcion f definida por f (x) = ax + b , permite

una interpretaci
on sencilla, es la funcion cuyo valor en x es igual al cambio producido en
el valor de f , cuando en el dominio se produce un incremento h. Entonces f (x) = ah
define una funcion constante, y si h = 1, la primera diferencia coincide con el valor a de la
pendiente de la recta.
3

Ejemplo 1.1.2. Sea u(x) = xm , m N m 1. Hallamos la primera diferencia teniendo

en cuenta que
m1
X
m
m
a b = (a b)
amk1 bk .
k=0

Entonces

[u(x)] = (xm ) = (x h)m xm =

= (x + h x)
=h

m1
X

m1
X

(x + h)mk1 xk =

k=0

(x + h)mk1 xk .

k=0

De esta expresion puede observarse que la primera diferencia de una potencia de x, es un

polinomio de grado menor en una unidad.
Una propiedad muy sencilla pero muy u
til al calcular diferencias de ciertas funciones como
las polinomicas, es la siguiente:
Teorema 1.1.1. Si f, g son funciones y a, b n`
umeros reales entonces
(1.2)

[af1 + bf2 ] (x) = af1 (x) + bf2 (x) .

Demostraci
on.
[af1 + bf2 ](x) = [af1 (x + h) + bf2 (x + h)] [af1 (x) + bf2 (x)] =
= [af1 (x + h) af1 (x)] + [bf2 (x + h) bf2 (x)] =
= a[f1 (x + h) f1 (x)] + b[f2 (x + h) f2 (x)] =
= af1 (x) + bf2 (x).

Ejemplo 1.1.3. Hallar la primera diferencia de la funcion u (t) = 3t2 5t.
Sean f (t) = t2 y g(t) = t entonces u(t) = [3f + 5g](x). Aplicando el Teorema 1.1.1 se
obtiene
[3f + 5g](x) = 3t2 5t =
= 3[(t + h)2 t] 5[t + h t] =
= 3[t2 + 2ht + h2 t2 ] 5[h] =
= 3t2 6ht + 3h2 3t2 + 5h =
= 5h + 6ht + 3h2 .
Por induccion matematica, el Teorema 1.1.1 puede extenderse a n funciones y n constantes,
como se indica a continuacion.

Corolario 1.1.1. sean f1 , f2 , , fn funciones y c1 , c2 , , cn constantes arbitrarias. Entonces para cualquier entero positivo n se cumple que
(1.3) [c1 f1 (x) + c2 f1 (x) + + cn f1 (x)] = c1 f1 (x) + c2 f2 (x) + + cn fn (x).
Demostraci
on. Procediendo por induccion se tiene: Para n = 2, seg
un el teorema precedente el resultado es valido.
Si el resultado se cumple para n = k, entonces tambien es cierto para n = k + 1. De hechor,
si
[c1 f1 (x) + c2 f2 (x) + + ck fk (x)] = c1 f1 (x) + c2 f2 (x) + + ck fk (x)
entonces
[c1 f1 (x) + c2 f1 (x) + + ck fk+1 (x)] = c1 f1 (x) + c2 f2 (x) + + ck+1 fk+1 (x).
c1 f1 (x) + + ck+1 fk+1 (x) = [c1 f1 (x) + + ck fk (x)] + ck+1 fk+1 (x)
= [c1 f1 (x) + + ck fk (x)] + ck+1 fk+1 (x)
= [c1 f1 (x) + + ck fk (x) + ck+1 fk+1 (x)].

Observaci
on 1. Del Teorema anterior se podra preguntar que condiciones debe tener

P
P
P
fi (x) para que se cumpla la siguiente igualdad
fi (x) = fi (x).?
i=1

i=1

i=1

Ejemplo 1.1.4. Primera diferencia del polinomio

P (x) = a1 + a1 x + + am xm ,
es de grado m 1.
Aplicando el Corolario 1.1.1 y el Ejemplo 1.1.2 se obtiene que
P (x) = a1 1 + a2 x + + am xm =
= 0 + ha2

0
X
k=0

(x + h)

1k k

x + ha3

1
X

(x + h)

2k1 k

k=0

x + + ham

m1
X

(x + h)m1k xk .

k=0

Por lo tanto se obtiene un grado menor en el polinomio P (x).

1.2.

Diferencia de Orden Superior

Como se establecio en la definicion precedente la primera diferencia de una funcion dada,

produce una nueva funcion. La primera diferencia de esta nueva funcion se llama diferencia
de orden dos y en general se puede definir de forma inductiva para cualquier n
umero natural
n, la diferencia de orden n.
5

Definici
on 1.2.1. La diferencia de orden cero de una funcion se define como
0 f (t) = f (t) .

(1.4)

Se suele identificar 0 con I el operador identidad, esto es

0 f = If = f.

(1.5)

Definici
on 1.2.2. La diferencia de orden n, n de una funcion f se puede definir como


n f (x) = n1 f (x) ,
para todo n
umero natural n 1.

Ejemplo 1.2.1. Sea f (x) = x3 + 1 entonces aplicando la definici

on (1.1.1) tenemos :
f (x) = f (x + h) f (x) = (x + h)3 + 1 x2 1 = 3x2 h + 3xh2 + h3 .
Aplicado nuevamente la definici
on 1.1.1 y el teorema 1.1.1 a la funcion diferencia tenemos:
(f )(x) = (3x2 h + 3xh2 + h3 ) =
= 3hx2 + 3h2 x + h3 =
= 3h[(x + h)2 x2 ] + 3h2 [x + h x] + 0 =
= 3h[(x + h)2 x2 ] + 3h2 h =
= 3h[(x + h)2 x2 ] + 3h3 =
= 3h[x2 + 2xh + h2 x2 ] + 3h3 =
= 3h[2xh + h2 ] + 3h3 =
= 6xh2 + 6h3 .
Por la Definicion 1.2.2 se tiene que
(f )(x) = 2 f (x) = 6xh2 + 6h3 ,
que es la segunda diferencia. Si calculamos la tercera diferencia tenemos que
((f (x))) = (2 f (x)) = 3 f (x) = [6xh2 + 6h3 ] =
= [6xh2 ] + [6h3 ] =
= 6h2 [x] + [6h3 ] =
= 6h3 .
Entonces la tercera diferencia es constante y por consiguiente, 4 f (x) = 6h3 6h3 = 0.
Anotemos que n f (x) = n1 [f (x)]. Inductivamente se tiene que
[f (x)] = 0 [f (x)].
6

Si
k f (x) = k1 [f (x)];
k+1 f (x) = [k f (x)] = [k1 (f (x))] = k [f (x)].
La primera diferencia de una funcion tambien tiene formulaciones de tipo general para productos y cocientes de funciones, las cuales tienen bastante similitud con las ya conocidas
para la derivacion de productos y cocientes. A continuacion se formulan y demuestran estos
resultados.
Teorema 1.2.1. Dadas las funciones f y g entonces:
1. [f g] (x) = f (x) g (x) + g (x) f (x + h) .
 
f
f (x) g (x) f (x) g (x)
2.
si g (x) g (x + h) 6= 0.
(x) =
g
g (x) g (x + h)
Demostraci
on.
1.
[f g](x) = f (x + h)g(x + h) f (x)g(x) =
= f (x + h)g(x + h) f (x + h)g(x) + f (x + h)g(x) f (x)g(x) =
= f (x + h)[g(x + h) g(x)] + g(x)[f (x + h) f (x)] =
= f (x + h)g(x) + g(x)f (x) = f (x)g(x) + g(x)f (x + h).
2.
 
f (x + h) f (x)
f

=
(x) =
g
g(x + h)
g(x)
=

g(x)f (x + h) g(x + h)f (x)

=
g(x + h)g(x)

g(x)f (x + h) f (x)g(x) + f (x)g(x) g(x + h)f (x)]

=
g(x + h)g(x)

g(x)[f (x + h) f (x)] + f (x)[g(x) g(x + h)]

=
g(x + h)g(x)

g(x)[f (x + h) f (x)] f (x)[g(x + h) g(x)]

=
g(x + h)g(x)

g(x)f (x) f (x)g(x)

.
g(x)g(x + h)
7


Ejemplo 1.2.2. Calcular la primera diferencia u(x) = 5x(1 x).
Sean f (x) = 5x y g(x) = (1 x) aplicando el teorema 1.2.1 la parte 1 se obtiene
u(x) = [f.g](x) = [5(x + h) 5x][1 x] + [(1 x h) 1 + x][5(x + h)]
= (5x + 5h 5x)(1 x) + (h)(5x + 5h)
= 5h 5hx 5hx 5h2
= 5h 10hx 5h2 .
x
+1
Sean f (x) = x y g(x) = x2 + 1, aplicando el Teorema 1.2.1 la parte 2 se obtiene
 
f
(x)(x2 + 1) x(x2 + 1)
v(x) =
=
(x) =
g
[(x + h)2 + 1][x2 + 1]
(x + h x)(x2 + 1) [(x + h)2 + 1 (x2 + 1)]
=
=
(x2 + 2hx + h2 + 1)(x2 + 1)
h(x2 + 1) x(x2 + 2hx + h2 + 1 x2 1)
=
=
(x2 + 2hx + h2 + 1)(x2 + 1)
h(x2 + 1) x(2hx + h2 )
= 2
=
(x + 2hx + h2 + 1)(x2 + 1)
hx2 + h 2hx2 + xh2
=
= 2
(x + 2hx + h2 + 1)(x2 + 1)
h hx2 + xh2
= 2
.
(x + 2hx + h2 + 1)(x2 + 1)
Ejemplo 1.2.3. Calcular la primera diferencia v(x) =

x2

Por induccion matematica se probara a continuacion que la diferencia de orden n es lineal.

Teorema 1.2.2. n [c1 f1 (x) + c2 f2 (x)] = c1 n f1 (x) + c2 n f2 (x) n = 0, 1, 2, 3
Demostraci
on. Procediendo por induccion, para n = 0 se tiene:
0 [c1 f1 (x) + c2 f2 (x)] = I[c1 f1 (x) + c2 f2 (x)] = c1 f1 (x) + c2 f2 (x),
por la definicion de 0 f = If = f .
Asumamos que para n = k el resultado es verdadero y demostremos que para n = k + 1
tambien lo es
k+1 [c1 f1 (x) + c2 f2 (x)] = hk [c1 f1 (x) + c2 f2 (x)]i = [c1 k f1 (x) + c2 k f2 (x)].
Aplicando el teorema 1.1.1 se tiene
hk [c1 f1 (x) + c2 f2 (x)]i = c1 k f1 (x) + c2 k f2 (x) = c1 k+1 f1 (x) + c2 k+1 f2 (x)
luego queda demostrado para n = k + 1, lo cual implica que se cumple para todo n.
8

1.3.

Operador de Traslaci
on E

Introduciremos una nueva funcion que facilitara los calculos, ya que el operador aplicado
a una funcion f requiere de dos operaciones para conseguir f (x). Si nosotros cambiamos
de f (x) a f (x + h) y substraemos f (x) del resultado, resulta ser sumamente u
til introducir
un smbolo especial para la primera operacion o funcion, que llamaremos traslacion y se
denotara por la letra E
Definici
on 1.3.1. Sean f una funcion y h cualquier constante tales que x y x + h pertenezcan al dominio de f , se define operador de traslacion de f como la funcion Ef, cuyo valor
en x es
Ef (x) = f (x + h).
Observando podemos escribir
f (x) = Ef (x) f (x) = Ef (x) If (x) = (E I)f (x).
El operador E, como , puede aplicarse mas de una vez a una funcion. Por ejemplo:
E 2 f (x) = E[Ef (x)] = Ef (x + h) = f (x + 2h),
E 3 f (x) = E[E 2 (x)] = Ef (x + 2h) = f (x + 3h),
y en general, si se define
(1.6)

E 0 f (x) = If (x) = f (x).

Se tiene la siguiente definicion para un entero n 1 .

Definici
on 1.3.2. Para un entero n 0 se define
E n f (x) = E(E n1 f (x))

n = 1, 2,

A continuacion demostraremos un resultado que nos permite hallar E n de una forma mas
simple, y a la vez nos va permitir demostrar algunos teoremas mas fuertes.
Lema 1.3.1.
E n f (x) = f (x + nh)

para todo n 1.

Demostraci
on.
Haciendo induccion sobre n, para n = 1 se tiene que E 1 f (x) = f (x + h) y esto es verdadero por la definicion.
Asumamos que
E k f (x) = f (x + kh).
Aplicando el operador E a E k f, se tiene:
E[E k f (x)] = E[f (x + kh),
9

usando la Definicion 1.3.1 y 1.3.2 se obtiene que

E k+1 f (x) = f (x + kh + h)
= f (x + [k + 1]h),
entonces el resultado es verdadero para todo n 1.

Igual que para n , el operador E n es lineal.

Teorema 1.3.1 (Propiedad de linealidad).
E n [c1 f1 (x) + c2 f2 (x)] = c1 E n f1 (x) + c1 E n f1 (x), c1 , c2 R, n = 0, 1, 2, 3
Demostraci
on.
Apliquemos induccion sobre n. Para n = 0 es verdadero por que
E 0 [c1 f1 + c2 f2 ](x) = I[c1 f1 + c1 f1 ](x) = c1 f1 (x) + c1 f1 (x)
y
c1 E 0 f1 (x) + c1 E o f1 (x) = c1 If1 (x) + c1 If1 (x) = c1 f1 (x) + c1 f1 (x)
por la definicion E 0 = I.
Supongamos que es verdadero para n = k y demostremos para n = k + 1.
Si
E k [c1 f1 (x) + c2 f2 (x)] = c1 E k f1 (x) + c2 E k f2 (x)
aplicando el operador E, se tiene
E(E k [c1 f1 (x) + c2 f2 (x)]) = E[c1 E k f1 (x) + c2 E k f2 (x)],
aplicando las definiciones y el Lema 1.3.1
E k+1 [c1 f1 (x) + c2 f2 (x)] = E(c1 f1 (x + kh) + c2 f2 (x + kh)
= c1 f1 (x + kh + h) + c2 f2 (x + kh + h)
= c1 f1 (x + h(k + 1)) + c2 f2 (x + k(k + 1))
= c1 E k+1 f1 (x) + c2 E k+1 f2 (x).
Lo cual demuestra que se cumple para n = k + 1.

Teorema 1.3.2. Sea f y g dos funciones. Entonces

(1.7)

[f g](x) = Ef (x) g(x) + g(x) f (x).

Demostraci
on. Por el teorema 1.2.1 parte 1 y la definicion del operador E se tiene
[f g] (x) = f (x + h) g (x) + f (x) g (x + h)
= Ef (x)g (x) + f (x) g(x).

10

Teorema 1.3.3 (Formula para la diferencia del cociente para dos funciones).
Sean f y g dos funciones entonces
 
f
g(x)f (x) f (x)g(x)
si g(x)Eg(x) 6= 0.
(x) =

g
g(x)Eg(x)
Demostraci
on. Por el teorema 1.2.1 parte 2 y la definicion del operador E se tiene
 
f (x) g (x) f (x) g (x)
f
(x) =

g
g (x) g (x + h)
f (x) g (x) f (x) g (x)
=
.
g (x) Eg (x)


1.4.

La Funci
on Factorial

Una funcion fundamental para calcular diferencias de polinomios es la factorial que se define
continuacion.
Definici
on 1.4.1. La funcion factorial de orden n(n = 1, 2, 3 ) es aquella cuyo valor en
x esta dado por
x(n) = x(x h)(x 2h) [x (n 1)h]

para n = 1, 2, . . . y x(0) = 1

Ejemplo 1.4.1.
x(1) = x

(x + h)(1) = x + h

x(2) = x(x h)

(x + h)(2) = (x + h)x

x(3) = x(x h)(x 2h)

(x + h)(3) = (x + h)(x)(x h)

De modo que, si h = 2,
3(1) = 3
4(1) = 4

3(2) = 3 1 = 3

4(2) = 4 2 = 8

3(3) = 3 1 1 = 3
4(3) = 4 2 0 = 0

Si h = 1
x(1) = x
x(2) = x(x 1) = x2 x

x(3) = x(x 1)(x 2) = x3 3x2 + 2x

En el siguiente teorema se estudia la primera diferencia de la funcion factorial.

11

Teorema 1.4.1. Si n es un entero positivo, entonces

x(n) = nhx(n1) .
Demostraci
on. Para n 1
x(n) = (x + h)(n) x(n) .
Teniendo en cuenta
(x + h)(n) = (x + h)(x)(x h) [x (n 2)h] =
= (x + h)x(n1) ,

y que
x(n) = x(x h)(x 2h) [x (n 2)h][x (n 1)h]
= [x (n 1)h]x(n1) ,

se sigue que
x(n) = (x + h)x(n1) [x (n 1)h]x(n1) =
= [(x + h) [x (n 1)h]]x(n1) =

= [x + h x + nh 1h]x(n1) =
= nhx(n1) .


Puede observase que excepto por una constante, la primera diferencia de la funcion factorial
es una funcion factorial de orden menor en una unidad que la funcion original. Entonces se
tiene el siguiente resultado para la diferencia de orden n de la funcion factorial de orden n.
Teorema 1.4.2. Si n es un entero positivo, entonces
n x(n) = n!hn
Demostraci
on. Por induccion matematica.
Para n = 1, 1 x(1) = h.
Suponga que k x(k) = k!hk . Entonces
k+1 x(k+1) = k [x(k+1) ] =
= k [(k + 1)hx(k) ] =
= (k + 1)hk x(k) =
= (k + 1)hk!hk =
= (k + 1)!hk+1 .

12

Teorema 1.4.3. sea f un polinomio de grado n, defenida por

a0 + a1 x + a2 x 2 + + a n x n ,

(1.8)

con a0 , a1 , , an constantes arbitrarias y an 6= 0. Entonces la enesima diferencia de f es

una funcion constante y las subsiguientes diferencias son iguales a cero, es decir,
p f (x) = 0

(1.9)

si p > n.

Demostraci
on. Por (1.8)
(1.10)

f (x) = a0 1 + a1 x + a2 x2 + + an xn .

Si n es un entero positivo,usando el teorema del binomio obtenemos

xm = (c + h)m xm =

m(m 1) m2 2
x
h + + hm x m =
2
m(m

1)
xm2 h2 + + hm .
= mxm1 h +
2

= xm + mxm1 h +

Es decir, aplicando a xm produce un n

umero finito de terminos con xm1 como la potencia
mas alta de x. Usando (1.10) se obtiene una suma de terminos. Cada constante multiplicada
por una potencia de x y la maxima potencia de x aparece n1. En otras palabras a aplicando
a un polinomio de grado n resulta un polinomio de grado (n 1). Por un razonamiento
similar 2 f produce un polinomio de grado n2y porconsiguiente si p > n p f (x) = 0. 
Si f (x) = xn con n N, f se puede escribir en terminos de funciones factoriales.
Ejemplo 1.4.2. Si n = 1, x(1) = x; si n = 2 x(2) = x(x h) = x2 xh y x2 = x(2) + x1 h.
Multiplicando la igualdad procedente por x entonces x3 = xx(2) + xx(1) h.
Como
x(3) = x(x h)(x 2h) =
= x(2) (x 2h) =

= xx(2) 2hx(2) ,
luego

x3 = x(3) + 2hx(2) + hx2 =

= x(3) + 2hx(2) + h[x(2) + x(1) h] =
= x(3) + 3hx(2) + h2 x(1) .
Se observa que x3 es igual x(3) mas una constante por x(2) mas una constante por x(1) .
Se probara en el caso general, que existen constantes A1 , A2 , A3 , . . . , An1 tales que xn =
x(n) +An1 x(n1) +An2 x(n2) + +, A2 x(2) +A1 x(1) y por consiguiente xn se puede expresar
en terminos de funciones factoriales.
13

Teorema 1.4.4. Si n N existen constantes A1 , A2 , A3 , . . . , An1 tales que

xn = x(n) + An1 x(n1) + An2 x(n2) + + A2 x(2) + A1 x(1) .
Demostraci
on.
Por induccion.
Si n = 2, x(2) = x(x h) = x2 xh entonces x2 = x(2) + x(1) h.
As A1 = h.
Suponganse que el resultado se cumple para k. Entonces existen constantes A1 , A2 , A3 , . . . , Ak1
tales que
xk = x(k) + Ak1 x(k1) + Ak2 x(k2) + + A2 x(2) + A1 x(1) .

Multiplicando la igualdad anterior por x se sigue que

xxk = xx(k) + Ak1 xx(k1) + Ak2 xx(k2) + + A2 xx(2) + A1 xx(1) .

Como
x(k+1) = x(x h)(x 2h) (x (k 1)h)(x kh)
= x(k) [x kh],

entonces
entonces

xx(k) = x(k+1) khx(k) ,

para todo k = 0, 1, . . .

xk+1 = x(k+1) khx(k) + Ak1 [x(k) + (k 1)hx(k1) ] + Ak2 [x(k1) + (k 2)x(k2) ] + +

+ A2 [x(3) + 2hx(2) ] + A1 [x(2) + hx(1) ] =

= x(k+1) + x(k) [kh + Ak1 ] + x(k1) [Ak1 (k 1)h + Ak2 ] + x(k2) [Ak2 (k 2)h + Ak3 ]+
+ + x(2) [A2 2h + A1 ] + A1 hx(1) .

Claramente se ve que lo que esta entre parentesis son constantes, por lo cual queda demostrado el teorema.

Teorema 1.4.5. Si f es un polinomio de grado n, f (x) = a0 + a1 x + a2 x2 + + an xn
entonces
n f (x) = n!hn an .
Demostraci
on. Tenemos que f (x) = a0 + a1 x + a2 x2 + + an xn y aplicando la diferencia
enesima tenemos
n f (x) = n a0 + n a1 x + n a2 x2 + + n an xn .
Por el teorema 1.4.3 tenemos que
n f (x) = n an xn ,
y por el teorema 1.4.2 se tiene que
n f (x) = n an xn = n!hn an .

14

El siguiente ejemplo que se presenta a continuacion es semejante al teorema de Taylor para

derivadas.
Ejemplo 1.4.3. Considere h = 1 y f (x) = x(n) + An1 x(n1) + An2 x(n2) + + A2 x(2) +
A1 x(1) + A0 .
f (x) = A1 + 2A2 x(1) + 3A3 x(2) + + nx(n1)

2 f (x) = 2 1 A2 + 3 2A3 x(1) + + n(n 1)x(n2)

3 f (x) = 3 2 1 A2 + + n(n 1)(n 2)x(n3)

..
.
n f (x) = n!x(0) .
Si x = 0, como 0(1) = 0(2) = 0n = 0 entonces

f (0) = A0 , f (0) = A1 , 2 f (0) = 2!A2 , 3 f (0) = 3!A3 , . . . , n f (0) = n!,

luego
(1.11)

f (x) = f (0) + f (0)x(1) + 2 f (0)

x(2)
x(3)
xn
+ 3 f (0)
+ + n f (0).
2!
3!
n!

Tambien se puede representar en terminos de coeficientes binomiales, ya que

 
n
n!
1 2 . . . (n k) (n k 1) . . . n
=
=
=
k
k!(n k)
k!(1 2 . . . n k)
n (n 1) . . . (n k + 2)(n k + 1)
=
=
k!
x(k)
.
=
k!
Por consiguiente
(1.12)

 
 
 
 
x
x
x
x
2
3
f (x) = f (0) +
f (0) +
f (0) +
f (0) + +
n f (0).
1
2
3
n

La ecuaciones (1.11) y (1.12) se llaman formulas de interpolacion de Newton, y se usan

para hallar aproximaciones de funciones para cualquier n
umero real.

1.5.

Igualdad de los Operadores

La ecuacion (1.3) se puede reescribir de la siguiente manera:

(1.13)

f (x) = Ef (x) If (x).

Seg
un esta ecuacion se obtiene el mismo resultado al aplicar el operador a cualquier
funcion, como aplicar el operador (E I). Casos como estos son muy comunes, se introduce
una terminologa especial para describirlos.
15

Definici
on 1.5.1. Dos operadores, O1 y O2 , son iguales si tienen un dominio com
un D y
para todo f D, O1 f = O2 f .
As (1.13) puede escribirse como una equivalencia
(1.14)

=EI

De igual manera despejando E se tiene que E = + I. Lo cual es equivalente a

Ef (x) = f (x) + If (x)
Teorema 1.5.1. Los operadores y E conmutan, es decir,
(1.15)

E = E.

Demostraci
on. Lo que vamos demostrar es que Ef (x) y Ef (x) son iguales.
Entonces
(Ef (x)) = f (x + h)
(1.16)

= f (x + 2h) f (x + h),

(1.17)

(E)f (x) = E[f (x)] = E[f (x + h) f (x)]

= f (x + 2h) f (x + h).

por (1.16) y (1.17) se tiene lo que se quera demostrar.

Teorema 1.5.2. Los operadores y E cumplen la ley de los exponentes, es decir, para n
y m enteros no negativos,
(1.18)
(1.19)
(1.20)

m n = n m = m+n ,
E m E n = E n E m = E m+n ,
m E n = E n m .

Demostraci
on. Para demostrar (1.18) se procede por induccion. Para n = 0 es verdadero
porque
m 0 f (x) = m [0 f (x)] = m If (x) = m f (x) = Im f (x) = 0 m f (x) = m+0 f (x).
luego
m 0 = 0 m = m+0 .
Supongamos que es verdadero para n = k y demostremos para n = k + 1.
Si
m k+1 f (x) = m [k+1 f (x)] = m 1 [k f (x)] = 1 m [k f (x)] =
= 1 k [m f (x)] = k+1 [m f (x)] = k+1 m f (x)
16

y
k+1 [m f (x)] = 1 k [m f (x)] = 1 [m+k f (x)] = m+k+1 f (x) = m+(k+1) f (x)
Por transitividad se concluye que
(1.21)

m k+1 f (x) = k+1 m f (x) = m+(k+1) f (x).

Ahora para demostrar (1.19) utilizaremos el lema 1.3.1

(1.22)

E m (E n f (x)) = E m f (x + nh) = f (x + nh + mh) = f (x + mh + nh),

y
(1.23)

(E n E m )f (x) = E n (E m f (x)) = E n f (x + mh) = f (x + mh + nh).

Por (1.22) y (1.23) se obtiene

(1.24)

(E m E n )f (x) = (E n E m )f (x).

De la Ecuacion (1.23) se tiene que

(1.25)

(E n E m )f (x) = E m+n f (x).

Haciendo transitividad en (1.24) y (1.25) se cumple lo pedido.

Para demostrar (1.20), Apliquemos induccion sobre n = 0. Para n = 0 es verdadero porque
(1.26)

E m 0 f (x) = E m If (x) = E m f (x),

(1.27)

0 E m = 0 f (x + hm) = If (x + hm) = f (x + hm).

luego por (1.26) y (1.27) se tiene E m 0 = 0 E m .

Asumamos que
(E m k )f (x) = (k E m )f (x).
Aplicando el operador , se tiene
E m [k f (x)] = k [E m f (x)],
E m [k f (x)] = k [E m f (x)],
(E m k+1 )f (x) = (k+1 E m )f (x).
Entonces se cumple para todo n 0, luego
E m n = n E m .

A continuacion presentaremos un resultado que nos permite escribir las diferencias en terminos del operador E y viceversa, es decir, como una combinacion lineal de E o como de .
17

Teorema 1.5.3. Si n es un entero positivo, entonces

 
n
P
n
n
i.) f (x) =
(1)nk E k f (x).
k
k=0
 
n
P
n
n
ii.) E f (x) =
k f (x).
k=0 k

Demostraci
on. i.) Sea f (x) cualquiera funcion. Entonces la primera diferencia f (x) es
f (x) = Ef (x) If (x). Por por Definicion 1.5.1 son equivalente, es decir
(E I).
As, para n N se tiene que

n (E I)n .

Aplicando el Teorema del binomio obtenemos

n  
X
n
n
(1 I)nk E k

k
k=0


n
X
n
(1)nk E k I nk .

k
k=0

Pero E k I nk E k , por lo tanto

n  
X
n

(1)nk E k .

f (x)

n  
X
n

(1)nk E k f (x).

Entonces
n

k=0

k=0

ii.) De igual forma como se hizo en lo anterior

E n ( + I)n .
Aplicando el Teorema del binomio a lo anterior se obtiene
n  
X
n
n
k I nk .
E
k
k=0
Pero k I nk = k , por lo tanto

n  
X
n

k .

E f (x)

n  
X
n

k f (x).

Entonces
n

k=0

k=0

18

A continuacion presentamos algunos ejemplos de aplicacion del anterior teorema

Ejemplo 1.5.1. Sea y(x) = cx , mostrar que
(1.28)

n x

c =c

n  
X
n

k=0

(1)nk ckh ,

con paso fijo o incremento h.

Aplicando el Teorema 1.5.3 parte i) a n cx obtenemos que
(1.29)

n  
X
n
c =
(1)nk E k cx .
k
k=0
n x

Por el Lema 1.3.1 se tiene que

E k ck = cx+hk = cx ckh .

(1.30)

Reemplazando (1.30) en (1.29) se tiene que

(1.31)

n x

c =

n  
X
n

k=0

(1)nk cx ckh .

Aplicando la propiedad homogenea de la sumatoria a (1.31) se obtiene que

n x

c =c

n  
X
n

k=0

(1)nk ckh .

Ejemplo 1.5.2. Sea la funcion y(x) = xm con m N. Vamos derivar la siguiente formula
X n
n m
(1)nk (x + kh)m ,
x =
k
k=0
con incremento h.
Aplicando el Teorema 1.5.3 parte i) a n xm obtenemos que
X n 
n m
(1.32)
x =
(1)nk E k xm ,
k
k=0
por el Lema 1.3.1 se tiene que
E k xm = (x + kh)m .

(1.33)

Reemplazando (1.33) en (1.32) se obtiene que

n m

x =

n  
X
n
k=0

(1)nk (x + kh)m .

El Teorema 1.5.3 puede ser usado para desarrollar problemas de interpolacion.

19

Ejemplo 1.5.3. Supongamos que se tiene la siguiente informaci

on referente a la funcion
y(x):
1. Es un polinomio de grado 3.
2. y(0) = 4, y(1) = 6, y(3) = 4.
3. Paso h = 1.
Estimar el valor de y(4).
Usando la informacion del punto 2 podemos construir la siguiente tabla.
x
0
1
2
3
4

y(x) y(x) 2 y(x) 3 y(x)

4
2
0
6
6
2
6
8
4
4
?

Por el Lema 1.3.1 Tenemos que

E k [y(x)] = y(x + hk).
Tomando x = 0, k = 4 y aplicando el Teorema 1.5.3 se obtiene que
4  
X
4
4
k y(0) =
E [y(0)] = y(4) =
k
k=0

= y(0) + 4y(0) + 62 y(0) + 43 y(0) + 4 y(0).

Como y(x) es un polinomio de grado 3 entonces 4 y(x) = 0.

Luego
y(4) = 4 + 4(2) + 6(0) + 4(6) + 0 = 12.

1.6.

Suma Indefinida: El Operador 1

En esta seccion se trata de introducir el operador de suma indefinida que puede entenderse
como el operador inverso al de la primera diferencia.
El nombre de suma indefinida puede intuirse de la siguiente observacion donde se asume
que las funciones consideradas son de variable entera no negativa. Si f es una funcion cuyo
n1
P
valor en n es f (n) , si la funcion F (n) =
f (k) . Tomando como paso h = 1,
k=1

F (n) = F (n + 1) F (n) =
=

n
X
k=1

f (k)

= f (n) +

n1
X
k=1

n1
X

f (k) =

k=1

f (k)

= f (n) + 0 = f (n).
20

n1
X
k=1

f (k) =

Puede notarse entonces, que dada la funcion f , F es una funcion cuya primera diferencia
es la funcion f. Escribiendo F (n) = 1 f (n), se tiene entonces que (1 f (n)) = f (n).
El comentario precedente induce la siguiente definicion.
Definici
on 1.6.1. Dada la funcion f, una funcion F cuya primera diferencia es f, sera llamada suma indefinida de f , y se denotar
a por F = 1 f.
En otros terminos si F y f son funciones tales que F (x) = f (x), F se llama suma de f,
o equivalentemente una suma de f es una soluci
on de la ecuaci
on F (x) = f (x).
El problema de encontrar una suma indefinida de f es el problema inverso de encontrar la
primera diferencia de F. Veamos el siguiente ejemplo.
Ejemplo 1.6.1. Una suma indefinida de la funcion f definida por f (x) = 3, es la funcion
F definida por F (x) = x + 4 , ya que considerando como paso h = 3 se tiene que
F (x) = F (x + 3) F (x) = (x + 3) + 4 (x + 4)
= x + 3 + 4 x 3 = 3.
Otra suma indefinida para la misma funcion es la definida por F (x) = x + c, ya que
tomando como paso h = 3, se tiene que
F (x) = F (x + 3) F (x) = (x + 3) + c (x + c)
= x + 3 + c x c = 3.
De lo anterior se puede observar, que si se consideran diferencias con un paso fijo h, una
suma indefinida de la funcion constante definida por f (x) = h, es la funcion F definida por
F (x) = x y de manera mas general una suma indefinida para f , es la funcion definida por
F (x) = x + c donde c es un n
umero real.
x(n+1)
+ C, con paso h 6= 0.
h(n + 1)
x(2)
= x(1) , lo cual implica que
dividiendo por 2h se tiene
2h
x(2)
1 (1)
x =
+ C.
2h

Ejemplo 1.6.2. Probemos que 1 x(n) =

Como x(2) = 2hx(1)

De igual forma
x(3)
= x(2)
3h
De forma general tenemos

entonces 1 x(2) =

1 x(n) =

x(3)
+ C.
3h

x(n+1)
+ C.
h(n + 1)

2x
Ejemplo 1.6.3. Hallemos que 1 2x = h1 + 4, con paso fijo h R.
2
Como
2x = 2x+h 2x = 2x (2h 1).
21

Entonces
(2h
Por consiguiente

1
2x = 2x .
1)

2x
+ C.
(2h 1)
Donde C es una constante cualesquiera, para este caso se toma C = 4.
Luego
2x
+ 4.
1 2x = h
(2 1)
1 2x =

Cuando f es una sucesion y se toma como paso h = 1 se tendra que una suma indefinida
n1
P
de f tambien es la funcion F definida por F (n) =
f (k) +c, donde c es una constante
k=1

real.
De una manera mas amplia, si q es una funcion periodica de periodo h, se tiene que q(x) =
0 y entonces una suma indefinida para la funcion f definida por f (x) = h, es la funcion
definida por F (x) = x + q(x). Se puede concluir entonces que si F es una suma indefinida
para una funcion f entonces para toda funcion periodica q de periodo h, G(x) = F (x) +
g(x) define una suma indefinida de f. Observese de otro lado que dada una funcion g y
considerando un paso fijo h, si g(x) = 0, entonces g(x + h) g(x) = 0 y por lo cual g es
una funcion periodica de periodo h. En consecuencia si se tienen dos sumas indefinidas para
la misma funcion su diferencia sera una funcion periodica de periodo h. Los comentarios
precedente permiten formular el siguiente teorema.
Teorema 1.6.1. Sean F1 y F2 sumas indefinidas de una funcion f con el mismo periodo
h, entonces su diferencia (F1 F2 ) es una funcion con periodo h,
Demostraci
on. Si F1 y F2 son sumas indefinidas de f, F1 (x) = f (x) y F2 (x) = f (x).
De aqu
F1 (x) F2 (x) = [F1 (x) F2 (x)] = 0.
Por el comentario mencionado F1 F2 es periodica de periodo h.

Este teorema nos lleva concluir que en general 1 6= 1 .

Si F es una suma indefinida de f , entonces F = f . De un lado 1 [F (x)] = F (x) + c es
una funcion cuya diferencia es F (x). [1 F (x)] = F (x) ya que 1 F (x) es una funcion
cuya diferencia es F. De esta forma:
1 [F (x)] [1 F (x)] = c y 1 6= 1 .

El operador de suma indefinida tambien cumple la propiedad de linealidad, la cual se demuestra a continuacion.
Teorema 1.6.2. Si F1 y F2 son sumas indefinidas de f1 y f2 . respectivamente con el mismo
paso, y si c1 y c2 son constante arbitrarias, entonces c1 F1 + c2 F2 es una suma indefinida de
la funcion c1 f1 + c2 f2 , es decir,
(1.34)

1 (c1 f1 + c2 f2 ) = c1 1 F1 + c2 1 F2 .
22

Demostraci
on. Tenemos que F1 = f1 y F2 = f2 . Debemos demostrar que la diferencia
de c1 F1 + c2 F2 es igual c1 f1 + c2 f2 . Usando (1.2) se obtiene
(c1 F1 + c2 F2 ) = c1 F1 + c2 F2
= c1 f1 + c2 f2 .
Esto muestra claramente que c1 F1 + c2 F2 es una funcion cuya diferencia es c1 f1 + c2 f2 y por
consiguiente:
1 [c1 F1 + c2 F2 ] = c1 f1 + c2 f2 .

Ejemplo 1.6.4. Hallemos
1 [x3 x2 x + 1].
Sabemos que
x2 = x(2) + hx(1) ,
x3 = x(3) + 3hx(2) + h2 x(1) .
Pasando x3 x2 x + 1 a funciones factoriales obtenemos
x(3) + (3h 1)x(2) + (h2 h 1)x(1) + x(0) ,
entonces aplicando la suma indefinida y el Teorema 1.6.2 se obtiene
1 [x(3) + (3h 1)x(2) + (h2 h 1)x(1) + x(0) ] = 1 x(3) + (3h 1)1 x(2) +
+ (h2 h 1)1 x(1) + 1 x(0) =

x(3)
x(2)
x(1)
x(4)
+ C1 + (3h 1)
+ C2 + (h2 h 1)
+ C3 +
+ C4 =
4h
3h
2h
h
x(3)
x(2) x(1)
x(4)
+ (3h 1)
+ (h2 h 1)
+
+ C1 + C2 + C3 + C4 ,
=
4h
3h
2h
h

haciendo C = C1 + C2 + C3 + C4 siguen siendo constantes.

Luego
1 [x(3) + (3h 1)x(2) + (h2 h 1)x(1) + x(0) ] =

x(3)
x(4)
+ (3h 1)
+
4h
3h
x(2) x(1)
+
+ C.
+ (h2 h 1)
2h
h

La formula (1.34) puede ser extendida a una combinacion lineal de n funciones, es decir,
(1.35)

n
X

ck fk =

k=1

n
X
k=1

la cual se puede demostrar por induccion.

23

ck 1 fk ,

Observaci
on 2. An
alogamente a la Observaci
on, que condiciones debe tener

P
P
que se cumpla la siguiente igualdad
1 fi (x) = 1
fi (x)?.
i=1

i=1

i=1

fi para

Teorema 1.6.3. Si F es una suma indefinida para F . Entonces

1 f (x + a) = F (x + a) + C,

para alguna constante C

Demostraci
on. Por hipotesis tenemos que F (z) = f (z) y haciendo z = x + a.
[F (z) + C] = F (z) + C = f (z) + 0 = f (z).
Aplicando la definicion de suma indefinida obtenemos
1 f (z) = F (z) + C,
luego
1 f (x + a) = F (x + a) + C.

En el siguiente teorema se estudia la suma indefinida de un producto, la cual conduce a la
llamada sumacion por partes.
Teorema 1.6.4 (Sumacion por partes). Sea u y v dos funciones, entonces
1 [v(x) u(x)] = u(x)v(x) 1 [Eu(x) v(x)].
Demostraci
on. Como por Teorema 1.2.1 parte 1
[u(x) v(x)] = Eu(x) v(x) + v(x) u(x).
Aplicando la definicion de suma indefinida se obtiene que
1 [Eu(x) v(x) + v(x) u(x)] = u(x) v(x)
y por el teorema 1.6.2 se concluye que
1 [Eu(x) v(x)] + 1 [v(x) u(x)] = u(x) v(x).
Despejando
1 [v(x) u(x)] = u(x) v(x) + 1 [Eu(x) v(x)].

Ejemplo 1.6.5. Por medio de la sumacion por partes vamos mostrar lo siguiente:
1 [x2x ] =

h2x+h
x2x

+ p(x),
2k 1 (2h 1)2

con paso fijo h R.

Sean v(x) = x, u(x) = 2x . Antes de aplicar de teorema de sumacion por partes, calculemos
24

1 cx .
Sea
cx = cx+h cx = cx (ch 1).
Entonces

 x 
c
cx
= h
.
c = h
c 1
c 1
x

Por consiguiente

1 cx =

(1.36)

cx
+ C.
ch 1

Reemplazando c = 2 en (1.36) se obtiene

1 2x =

2x
+ C.
2h 1

2x
. Entonces aplicando teorema 1.6.3
2h 1
 x+h 
1 h2
+C =

2h 1
h
h
1 [2x+h ] + C =
(2 1)
h
2x+h
h
h
+C =
(2 1) (2 1)
h2x+h
h
+ C.
(2 1)2

Como u(x) = 2x entonces 1 2x = u(x) =

x

[x2 ] =
1

=
=
=

x2x
2h 1
x2x
2h 1
x2x
2h 1
x2x
2h 1

Teorema 1.6.5 (Suma de series finitas). Sea f una funcion con dominio en los enteros
positivos y sea F cualquier suma indefinida de f . Entonces
n
X
k=1

f (k) = F (n + 1) F (1).

Demostraci
on. Como f tiene dominio en los naturales tenemos paso h = 1.
Por hipotesis tenemos que F (k) = f (k). Entonces
n
X
k=1

f (k) =

n
X

F (k) =

k=1

n
X
k=1

F (k + 1) F (k),

utilizando la propiedad telescopica en la expresion anterior, se obtiene que

n
X
k=1

f (k) = F (n + 1) F (1).


Ejemplo 1.6.6. Aplicando el Teorema anterior vamos obtener las siguientes formulas:
25

1.

n
P

k=

k=1

2.

n
P

n(n + 1)
.
2

ark1 =

k=1

3.

n
P

k=1

a arn
1r

(r 6= 1).

(2k 1) = n2 .

Tomando h = 1
1. Sea F (k) = f (k) = k, entonces
1 k = 1 k (1) =

k(k 1)
k (2)
= Y (k) =
.
2h
2h

Aplicando el teorema 1.6.5 (suma de series finitas) se obtiene

n
X

k=

k=1

k(k 1)
(n + 1)(n + 1 1)
n(n + 1)
=
0=
.
2
2
2

2. Sea Y (k) = y(k) = ark1 entonces

1 ark1 =

ark1
= Y (k).
rh 1

Luego
n
X

ark1 =

n
X
ark1
k=1

k=1

r1

arn+11 ar11

=
r1
r1

arn
a
arn
a

=
+
=
r1 r1
1r 1r
a
arn
a arn
=

=
.
1r 1r
1r
=

3. Sea Y (k) = y(k) = 2k 1 = 2k (1) k (0) entonces

1 [2k (1) k (0) ] =

2k (2)
k (1) = k (2) k (1) = k(k1)k = k 2 kk = k 2 2k = Y (x),
2

entonces
n
X
(2k 1) = k 2 2k = [(n + 1)2 2(n + 1)] [1 2] = n2 + 2n + 1 2n 2 + 1 = n2 .
k=1

Luego

n
X
k=1

(2k 1) = n2 .
26

1.6.1.

Aplicaci
on al C
alculo Integral

Existe en el calculo integral situaciones donde es imposible hallar el valor exacto de una
integral, ya que es necesario conocer una antiderivada de la funcion. Sin embargo, a veces
es difcil o imposible, hallar estas antiderivadas.
Por ejemplo evaluar con exactitud las siguientes integrales
Z 1
Z 1
x2
1 + x3 dx.
e dx,
0

Rb

Supongamos que queremos hallar la a f (x)dx. Dividamos el intervalo [a, b] en n subinter(b a)

valos de igual longitud x =
y sea a = x0 < x1 < x2 < < xn1 < xn puntos
n
intermedios del intervalo, Podemos escribir la siguiente aproximacion discreta de nuestra
integral
Z b
n
X
f (xi )xi .
f (x)dx
a

i=1

Aplicando el Teorema Suma de series finitas obtenemos

Z b
f (x)dx 1 f (xn )xn
a

En ocasiones es difcil de calcular suma indefinida de una funcion, una forma de evitar esto,
es sustituir la funcion f por otra funcion g que la aproxime y la cual facilite los calculos.
2

Por ejemplo, conocido el desarrollo de la serie de Taylor ex =

k=0

x2k
k!

para todo x R,

como la serie converge uniformemente, se deduce que:

Z 1
Z 1X

Z 1 2k

n
X
X
X
x2k
1
1
x
x2
e dx =
dx =
dx =

k!
(2k + 1)k!
(2k + 1)k!
0
0 k=0 k!
k=0 0
k=0
k=0

donde se puede hallar una aproximacion hasta termino enesimo por medio del operador
suma indefinida, ahora:


Z 1
1
x2
1
.
e dx
(2n + 1)n!
0
Se puede observar que este operador suma indefinida es difcil de calcular.

P
1
Considerando desarrollo de Taylor de funcion
=
xk , se deduce que
x 1 k=0
Z 2
Z 2X

Z 2

X
X
1
2k+1 1
k
k
x dx =
x dx =
dx =
.
k
+
1
1 x1
1
1 k=0
k=0
k=0
Tomando suma de los primeros n terminos, tenemos: tenemos
 n+1

Z 2
2
1
1
1
dx
.
n+1
1 x1
Lo cual es una suma indefinida facil de calcular.

Para el calculo diferencial, en Ejemplo 2.1.3 se hara una relacion entre el operador diferencial
y diferencia.
27

CAPITULO 2

LINEAL EN DIFERENCIAS
ECUACION
DE PRIMER ORDEN.

2.1.

Introduccion

Se comienza el estudio de las ecuaciones en diferencias con el analisis detallado de un tipo

especial de ecuacion, llamado ecuacion lineal en diferencias de primer orden. Expondremos
algunos ejemplos iniciales que permitiran ilustrar que es una ecuacion en diferencias de
primer orden y como una ecuacion de este tipo puede originarse.
Ejemplo 2.1.1. Este ejemplo est
a relacionado con el llamado interes simple e interes compuesto. Una cierta cantidad de dinero gana interes simple a una tasa de interes i en un
periodo determinado, cuando el capital obtenido en un periodo dado, es igual al capital en
el periodo anterior, mas el interes ganado durante el periodo por el capital inicial invertido.
Si el capital inicial se denota por C, y el capital en el periodo k se denota por Ck , entonces
Ck+1 = Ck + iC, k = 0, 1, , donde C0 = C.
El interes compuesto es obtenido cuando los intereses son agregados al capital. Si i representa la rata o tasa de interes por periodo de conversion (por ejemplo si se tiene una tasa anual
del 8 % capitalizable semestralmente i = 0,04, si la capitalizaci
on se hace trimestralmente
i = 0,02 etc) entonces el capital Ck+1 , acumulado en el periodo k + 1 de conversion es igual
al capital Ck obtenido en el periodo k mas iCk . Entonces Ck+1 = Ck + iCk , k = 0, 1, ,
donde C0 = C.
Las ecuaciones Ck+1 = Ck + iC, y Ck+1 = Ck + iCk , donde k = 0, 1, 2, son ejemplos de
ecuaciones en diferencias de primer orden.
Ejemplo 2.1.2. En este ejemplo se muestra la relaci
on entre la ecuaciones en diferencia de
ex
, esta
primer orden y las series de potencias. Considerando la funcion f definida por
1x
28

funcion tiene una expansion en series de potencias por

an xn , la cual es absolutamente

n=0

convergente en el intervalo (1, 1).

Al multiplicar por 1 x se obtiene

X
X
X
xn
n
n
e = (1 x)
an x =
(an an1 )x = a0 +
(an an1 )x =
n!
n=0
n=0
n=1
n=1
x

X
n=0

an x

an x

n+1

X
n=0

n=0

an x

an1 x = a0 +

n=1

X
n=1

(an an1 )x =

X
xn
n=1

n!

De aqu se concluye que

an an1 =
La ecuaci
on an an1 =

1
n!

n = 1, 2, . . . y a0 = 1.

1
es una ecuaci
on en diferencias de primer orden.
n!

Ejemplo 2.1.3. En este ejemplo se muestra la relaci

on cercana entre el operador de diferencia finita y el operador derivada, no es sorprendente que se puedan establecer conexiones entre ecuaciones en diferencias y ecuaciones diferenciales. Este ejemplo muestra esta
situacion.
Considere la ecuaci
on diferencial
(2.1)

dy(x)
= Ay(x) + B, para a x b,
dx

A y B son constantes con A 6= 0. Si y(a) = y0 , resolver el problema de valor inicial implica

encontrar una funcion derivable que en x0 = a tenga el valor presente y0 .
Para pasar de una ecuaci
on (2.1) a una ecuaci
on en diferencias se divide el intervalo [a, b]
ba
. Esto produce (n+1) partes x0 = a, x1 =
en n subintervalos, considerando el paso h =
n
x0 + h, x2 = x0 + 2h, . . . , xn = x0 + nh = b, encontrando yk = y(xk ) = y(x0 + kh).
Recordando que
y(x + h) y(x)
dy(x)
= lm
,
0
dx

parece razonable reemplazar la ecuaci

on (2.1) por
yk
= Ayk + B,
h

k = 0, 1, 2, . . . , n 1

o
yk+1 yk = Ahyk + Bh.
Despejando yk+1 , se tiene que
yk+1 = (1 + hA)yk + Bh

k = 0, 1, 2, . . . , n 1 con y0 = y(a).

La ecuaci
on yk+1 = (a + hA)yk + Bh es una ecuaci
on lineal en diferencia de primer orden
con coeficientes constantes.
29

Despues de haber estudiado lo ejemplos anteriores formularemos a continuacion la definicion

de una ecuacion en diferencia de primer orden.
Definici
on 2.1.1. Una ecuaci
on en diferencias de la forma
(2.2)

yk+1 + a (k) yk = b(k),

donde a y b son funciones definidas sobre un conjunto S, se llama ecuaci

on lineal en
diferencias de primer orden, sobre el conjunto S.
Definici
on 2.1.2. Cuando en (2.2), b(k) es la funcion nula la ecuaci
on llama ecuaci
on
homogenea.
Ejemplo 2.1.4. Ahora retomemos los ejemplos anteriores y miremos cuales de ellas son
ecuaciones lineales de primer orden, seg
un la definici
on:
1. La relacionada con interes simple, llevada a la forma de la ecuaci
on (2.2) es
Ck+1 Ck = iC

k = 0, 1, 2, . . .

donde a(k) = 1 y b(k) = iC, por lo tanto es una ecuaci

on lineal y no homogenea ya
que b(k) 6= 0.
2. La relacionada con el interes compuesto, llevada a la forma de la ecuaci
on (2.2),
conduce a
Ck+1 (i + 1)Ck = 0,

k = 0, 1, 2, . . .

donde a(k) = (i + 1) y b(k) = 0, por lo tanto es una ecuaci

on lineal homogenea, ya
que b(k) = 0.
3. En el Ejemplo 2.1.2, la ecuaci
on se puede llevar a la forma de la ecuaci
on (2.2)
obteniendo
1
k = 0, 1, 2, . . . .
an an1 =
n!
1
on lineal no homogenea, ya que
Donde a(k) = 1 y b(k) = , lo cual es una ecuaci
n!
b(k) 6= 0.
4. En el Ejemplo 2.1.3, la ecuaci
on se puede llevar a la forma de la ecuaci
on (2.2),
obteniendo lo siguiente:
yk+1 (a + hA)yk = B.
Donde a(k) = (a + hA) y b(k) = B, lo cual es una ecuaci
on no homogenea, ya que
b(k) 6= 0.

30

2.2.

Soluci
on de una Ecuaci
on Lineal en Diferencia de
Primer Orden

Ya definimos que es una ecuacion lineal en diferencias de primer orden y algunos situaciones
donde se pueden presentar este tipo de ecuaciones, ahora vamos a mirar como se comportan
y cual es su solucion general.
La idea de solucion para una ecuacion en diferencias sobre un conjunto de S es basicamente
la misma que se tiene para una ecuacion algebraica, o para una ecuacion diferencial como
se puede ver en los siguientes ejemplos.
Ejemplo 2.2.1. Tomemos la ecuaci
on en diferencias
yk+1 2yk = 0, k = 0, 1, 2, .
Claramente podemos ver que es una ecuaci
on lineal homogenea. Probemos que al reemplazar
k
yk = 2 k = 0, 1, 2, , en la ecuaci
on se obtiene una identidad.
yk+1 2yk = 2k+1 2(2k )
= 2k+1 2k+1
= 0.

Se dice que la funcion y definida por yk = y(k) = 2k es una soluci

on. Tambien podemos
k
tomar yk = 3 2 k = 0, 1, 2, , y reemplazando tenemos
yk+1 2yk = 3 2k+1 2(3 2k )
= 3 2k+1 3 2k+1
= 0.

En general se puede verificar que yk = C2k k = 0, 1, 2, , con C R, es una soluci

on y
que para C1 , C2 en R
yk = C1 2k + C2 2k ,
define una soluci
on.
yk+1 2yk = (C1 2k+1 + C2 2k+1 ) 2(C1 2k + C2 2k )

= (C1 2k+1 + C2 2k+1 ) (C1 2k+1 + C2 2k+1 )

= 0.

Lo anterior nos muestra que una ecuaci

on en diferencias puede tener muchas soluciones
sobre un conjunto especifico de valores.
Ejemplo 2.2.2. Miremos que la ecuaci
on
(2.3)

yk =

C
1 + Ck

k = 0, 1, 2,
31

es una soluci
on de la ecuaci
on
(2.4)

yk+1 =

yk
.
1 + yk

De la ecuaci
on (2.3) obtenemos
yk+1 =

C
.
1 + C(k + 1)

Ahora reemplazando la ecuaci

on (2.3) en (2.4) tenemos

yk+1

C
C
yk
C
1 + Ck
.
= 1 + Ck =
=
=
C
1 + C(k + 1)
1 + yk
1 + C(k + 1)
1+
1 + Ck
1 + Ck

Luego hace verdadera la igualdad.

Tomando estos ejemplos como soporte podemos definir a continuacion que es una solucion
de una ecuacion en diferencias finitas.
Definici
on 2.2.1. Una funcion f es una soluci
on de la ecuaci
on en diferencia
yk+1 = a(k)yk + b(k)
en un conjunto S si
fk+1 = a(k)fk + b(k),

para todo k S.

Ejemplo 2.2.3. La funcion y dada por

(2.5)

yk = 1

2
k

k = 1, 2, 3,

es una soluci
on de la ecuaci
on diferencias de primer orden
(2.6)

(k + 1)yk+1 + kyk = 2k 3

k = 1, 2, 3,

Para probar lo anterior se sustituyen los valores y en (2.6) para todo entero positivo k. De
(2.5) nosotros tenemos
(2.7)

yk+1 = 1

2
.
k+1

Reemplazando (2.5) y (2.7) en (2.6) se obtiene





2
2
(k + 1) 1
+k 1
= 2k 3,
k+1
k
simplificando al lado izquierdo tenemos
(k + 1) 2 + k 2 = 2k 3.
Lo cual implica yk satisface la ecuaci
on en diferencia (2.6).
32

Ejemplo 2.2.4. La funcion y dada por la

(2.8)

yk =

k(k 1)
+ C,
2

donde C es cualquier constante arbitraria, es soluci

on de la ecuaci
on en diferencias de
primer orden
(2.9)

yk+1 yk = k

k = 1, 2, 3,

Para probar lo anterior se sustituye los valores y en (2.9). De (2.8) se conoce que
(2.10)

yk+1 =

k(k + 1)
+ C.
2

Reemplazando (2.8) y (2.10) en (2.9) se obtiene



k(k + 1)
k(k 1)
k2 k
k2 k
+C
+C =
+ +C
+ C = k,
2
2
2
2
2
2
Simplificando el lado izquierdo se tiene
2k
k2 + k k2 + k
+C C =
= k,
2
2
lo cual demuestra que y es soluci
on de la ecuaci
on en diferencias (2.9).
En la seccion siguiente se estudia la ecuacion homogenea, la cual es fundamental a la hora
obtener soluciones para la ecuacion no homogenea.

2.3.

Propiedades de la Ecuaci
on Homog
enea

Denotaremos R como el conjunto de todas la sucesiones reales, el cual es un espacio

vectorial bajo las operaciones naturales de suma de sucesiones y multiplicaciones de escalares reales. A continuacion vamos demostrar que el conjunto solucion s de la ecuacion
homogenea es un subespacio R .
Teorema 2.3.1. El conjunto soluci
on de la ecuaci
on lineal homogenea en diferencias de
primer orden es un subespacio del espacio vectorial R .
Demostraci
on.
Se va demostrar que si u, v s y R, entonces (u + v) s y (u) s .
i) Como u, v s ,

uk+1 a(k)uk = 0 y vk+1 a(k)vk = 0.

Sumando las dos igualdades tenemos

(uk+1 + vk+1 ) a(k)(uk + vk ) = 0,
lo cual nos dice que u + v tambien es un elemento de s .
33

ii) Como u s entonces

uk+1 a(k)uk = 0,

multiplicando la igualad anterior por se tiene

[uk+1 a(k)uk ] = uk+1 a(k)uk = 0 = 0,
lo cual implica que u pertenece s .
Por i) y ii) se obtiene que s es subespacio de R .

Como consecuencia del teorema anterior se cumple el siguiente principio de superposicion

Corolario 2.3.1. Si u y v son dos soluciones de la ecuaci
on
yk+1 = a(k)yk
y c1 , c2 R, entonces c1 u + c2 v tambien es soluci
on.
Teorema 2.3.2 (Existencia y unicidad). Dada la ecuaci
on lineal homogenea en diferencia
de primer orden
(2.11)

yk+1 = a(k)yk con a(k) 6= 0 para todo k

y la constante m0 , existe una u

nica soluci
on, f con la condici
on f (0) = m0 .
Demostraci
on.
Se observa que una solucion f debe satisfacer f (1) = a(0)f (0), f (2) = a(1)f (1) = a(1)a(0)f (0);
f (3) = a(2)f (2) = a(2)a(1)a(0)f (0). Esto nos lleva a proponer y demostrar que la funcion
definida por*
f (k) = f (0)

k1
Y

k = 1, 2, . . .

a(j),

j=0

es una solucion. En efecto:

f (k + 1) = f (0)

k
Y

"

a(j) = f (0)

j=0

k1
Y

a(j) a(k) = a(k)f (k).

j=0

Para tener bien definida la expresi

on siguiente en k = 0 hay que suponer
1
Y

a(j) = 1.

j=0

En adelante se adoptara lo siguiente

(
a(r)a(r + 1) a(s.),
a(j) =
1
,
j=r
s
Y

para garantizar su validez para todo k 0.

34

sr
s<r

Si f (0) = m0 , f es una solucion de la ecuacion con la condiciones requerida.

De otro lado si g es una solucion,
a(k)g(k)
g(k + 1)
= k
,
k
Q
Q
a(j)
a(j)
j=0

entonces

j=0

g(k)
g(k + 1)
g(k + 1)
a(k)g(k)
g(k)

k1
= 0,
k1
k
k
k1
Q
Q
Q
Q
Q

a(k)
a(j)
a(j)
a(j)
a(j)
a(j)
j=0

j=0

j=0

j=0

Entonces existe una constante C tal que g(k) = C

k1
Q

j=0

a(j). De aqu se sigue que si g(0) = m0

j=0

entonces g = f y se tiene la unicidad.

Teorema 2.3.3. la dimension del subespacio de soluciones de la ecuaci

on lineal homogeneas
de primer orden es 1.
Demostraci
on. Sea
H :R s
7 H() =

n1
Y

a(k).

k=0

Probemos que H es una transformacion lineal biyectiva.

i) Sea x1 , x2 R entonces
H(x1 + x2 ) = (x1 + x2 )

n1
Y

a(k) = x1

n1
Y

a(k) + x2

a(k) = H(x1 ) + H(x2 )

k=0

k=0

k=0

n1
Y

ii) Sea , x1 R entonces

H(x1 ) = x1

n1
Y

a(k) = H(x1 ).

k=0

iii) H es inyectiva ya que si

H() = 0 =

k1
Y
j=0

como a(j) 6= 0 entonces = 0.

35

a(j),

iv) H es sobreyectiva ya que si g(k) = C

k1
Q
j=0

a(j) con C R es una solucion de la ecuacion

homogenea, entonces g es imagen de C.

Se concluye por i) y ii) que H es una transformacion lineal biyectiva. Por iii) y iv) dim R =
dim s = 1.

Ejemplo 2.3.1. Resolver la ecuaci
on con condici
on inicial
yk+1 = kyk ,

y(0) = 2.

Como a(k) = k, entonces la soluci

on general es
yk = y0

k1
Y
j=0

j = C(k 1)! k = 1, 2, 3, . . .

Teniendo en cuenta la condici

on inicial, la soluci
on es
yk = 2

k1
Y
j=0

j = 2(k 1)!.

Ejemplo 2.3.2. Retomando el Ejemplo 2.1.1 de interes compuesto hallemos la soluci

on
general de esta ecuaci
on.
Sea
Ck+1 = Ck (i + 1),
donde a = (i + 1). Entonces
Ck = C0

k veces

k1
Y
j=0

}|
{
z
(i + 1) = C0 (1 + i) (1 + i) (1 + i) . . . (1 + i) = C0 (1 + i)k .

Luego la soluci
on es:
Ck = C0 (1 + i)k .

2.4.

Soluci
on General de una Ecuaci
on no Homog
enea.

El Teorema que se enunciara a continuacion nos permitira ver la relacion que hay entre la
solucion particular de la ecuacion lineal completa y la solucion general de la homogenea
asociada a esta.
Teorema 2.4.1. Dada la ecuaci
on lineal no homogenea
yk+1 = a (k) yk + b(k).
Si f y g son soluciones de esta ecuaci
on, entonces la funcion s = f g es soluci
on de la
ecuaci
on homogenea asociada.
36

Demostraci
on. Por hipotesis
f (k + 1) a (k) f (k) = b(k)
y
g(k + 1) a (k) g(k) = b(k),
entonces
s(k + 1) = f (k + 1) g(k + 1) = a (k) f (k) + b(k) a (k) g(k) b(k) =
= a(k)[f (k) g(k)] = a(k)s(k).

El teorema muestra que si se conoce una solucion particular g, de la ecuacion no homogenea,
cualquier solucion f de ella, es la suma de esta solucion y una solucion de la ecuacion
homogenea asociada.
Ejemplo 2.4.1. Tomando la ecuaci
on
yk+1 yk = 1,
donde facilmente se puede ver que yk = k es soluci
on.
Ahora hallemos la soluci
on general de la ecuaci
on homogenea asociada, la cual es:
yk = y0

k1
Y

1 = y0 .

j=0

Miremos que s(k) = k + y0 es tambien una soluci

on. Reemplazando en yk+1 yk = 1
obtenemos que
y0 + k + 1 y0 k = 1,
lo cual cumple la igualdad.
Teorema 2.4.2. La soluci
on general de la ecuaci
on lineal no homogenea
yk+1 = a(k)yk+1 + b(k)
tiene la forma

# k1
X b(j)
,
yk = A
a(j) +
a(j)
j
Q
j=0
j=0
j=1
a(i)
k1
Y

"k1
Y

donde A es una constante.

Demostraci
on. Dividiendo ambos lados de
yk+1 = a(k)yk+1 + b(k),
37

i=0

por

k
Q

a(j) se tiene que

j=0

yk+1
k
Q

a(j)

j=0

entonces

a(k)yk
b(k)
a(k)yk
b(k)
+ k
=
,
+ k
k
k1
Q
Q
Q
Q
a(j)
a(j)
a(j)
a(k)
a(j)

j=0

j=0

yk+1
k
Q

a(j)

j=0

Por consiguiente

yk
k1
Q

a(j)

j=0

k1
Q

b(k).

a(j)

yk
yk+1
1

b(k).
=

k
k1
k
Q
Q
Q
Q
k1

a(j)
a(j)
a(j)
a(j)

j=0

yk

1
k
Q

j=0

Entonces

j=0

j=0

j=0

j=0

es una suma indefinida de

a(j)

j=0

j=0

b(k)
. Entonces por le teorema se series
k
Q
a(j)

j=0

finitas existe una constante A tal que:

k1
X
b(j)
=A+
.
j
k1
Q
Q
j=1
a(j)
a(i)

yk

j=0

De aqu

i=0

# k1
X b(j)
a(j) +
a(j)
yk = A
.
j
Q
j=1
j=0
j=0
a(i)
"k1
Y

k1
Y

i=0

Puede observase que la estructura de la solucion es la suma de una solucion de la ecuacion

homogenea y una particular de la no homogenea.
Ejemplo 2.4.2. Resolver
1. yk+1 = yk + (k + 1)2 . En este ejemplo a(k) = 1 y b(k) = (k + 1)2 .
Entonces :

yk = A

k1
Y

1+

j=1

=A+

k1
X
j=1

k1
X
(j + 1)2

1
=
j
Q

j=1
j=1
i

k1
Y

i=1

k1

X
(j + 1)2
(j + 1)2 .
=A+
1
j=1
38

Ahora aplicando propiedades de la sumatoria se tiene

yk = A +

k1
X

(j + 1) = A +

j +

k1
X

2j +

k1
X

1=

j=1

j=1

j=1

j=1

=A+

k1
X

2k 3 + 3k 2 + k 6
k(k 1)(2k 1)
+ k(k + 1) + (k 1) = A +
.
6
6

Luego

2k 3 + 3k 2 + k 6
yk = A +
.
6

2. yk+1 = (k a)yk . En este ejemplo a(k) = (k a) y b(k) = 0. Entonces poe el Teorema

2.4.2 la soluci
on general es:

k1

X
0

a(j)
a(j) +
yk = A
=
j
Q

j=1
j=0
j=0
(i a)
k1
Y

k1
Y

=A

i=1

k1
Y
j=0

(j a) + 0 = A

k1
Y
j=0

(a j) =

= A[a(a 1)(a 2) (a (k 1))] = Aa(k) .

3. Dado yk+1 = (k a)(k b)yk . Se observa que a(k) = (k a)(k b) y b(k) = 0.
Entonces la soluci
on general es:

k1
k1
X
Y
yk = A (k a)(k b) + (j a)(j b)
j
Q
j=0
j=0
j=1
k1
Y

=A

i=0

k1
Y
j=0

(i a)(i b)

(a j)(b j) =

= A[a(a 1)(a 2) (a hk 1i)] [b(b 1)(b 2)((b hk 1))] =

= Aa(k) b(k) .

4. Dado yk+1 = (k a1 )(k a2 ) (k ak )yk entonces la soluci

on general se deriva al
observar los casos 3 y 4, luego se obtiene que
(k) (k)

(k)

(k)

yk+1 = Ca1 a2 ak1 ak .

39

2.5.

Ecuaci
on Lineal en Diferencia con Coeficientes
Constantes

En esta seccion hallaremos la solucion general de un caso particular de la ecuacion no

homogenea, donde a(k) y b(k) son constantes.
Teorema 2.5.1. Dada la ecuaci
on lineal en diferencia de primer orden
yk+1 = ayk+1 + b,
con a y b constantes. Entonces la soluci
on general tiene la forma

k 1 ak
Aa +
b,
(a 6= 1)
yk =
1a

A + kb,
(a = 1)

Demostraci
on. Aplicando el teorema 2.4.2 tenemos

(2.12)

" k1  
#
k1
X
X 1 (j+1)
1

k
k
yk = A
.
a+
a
b = Aa + a b
k
a
j=0 Q

j=0
j=0
j=0
a
k1
Y

k1
Y

j=0

Aplicando en la sumatoria de la ecuacion anterior el teorema de suma series finitas de igual

forma como se hizo en el ejercicio 1.6.6 numeral 2 tenemos
k1  j+1
X
1

(2.13)

j=0

1 ak
.
ak (1 a)

Reemplazando (2.13) en (2.12) obtenemos



1 ak
1 ak
k
k
k
=
Aa
+
b.
yk = Aa + a b
ak (1 a)
1a
Si a = 1 tenemos

" k1 #
k1
X
X
1

yk = A
1+
1
b
=
A
+
b
1 = A + kb.

k
j=0 Q

j=0
j=0
j=0
1
k1
Y

k1
Y

j=0

Ejemplo 2.5.1. El llamado fenomeno de telara

na en economa, proporciona ejemplos en
los cuales se puede aplicar las ecuaciones en diferencias de primer orden. Consideremos la
siguiente situacion: Los agricultores deciden, sobre los precios del presente a
no, el area que
van a sembrar de cuerdo al producto.
40

Anticipando que el nivel de los precios se mantendr

a si los precios del a
no son altos, los
agricultores tienden a sembrar grandes areas. El siguiente a
no cuando la cosecha es recogida
y llevada al mercado, la oferta supera a la demanda, los precios caen, entonces los cultivadores sembraran menos area. Cuando la siguiente cosecha es recogida la oferta podra ser
menos que la demanda y los precios se incrementan, entonces los agricultores siembran mas
y el proceso continua.
Para el analisis de este proceso estudiaremos tres funciones para k = 0, 1, 2, . . . Como S,
donde S(k) es el n
umero de unidades ofrecidas en el periodo k, D es la funcion para la cual
D(k) es el numero de unidades demandadas en el periodo k y la funcion P, tal que P (k) es
el precio por unidades en el periodo k. Se hacen los siguientes supuestos:
i. La demanda es una funcion del precio D(k) = f (P (k)), es decir, la cantidad demandada
depende del precio.
ii. La oferta S es funcion del precio S(k +1) = g(P (k)), la cantidad ofrecida en un periodo
dado es funcion del precio en el periodo anterior.
iii. El precio del mercado esta determinado por la disponibilidad de oferta y las transacciones que se presentan en el precio en el cual la cantidad demandada son iguales, es
decir, P es determinado como la soluci
on de la ecuaci
on S(k) = D(K)
Si P (0) esta dada f y g son conocidas se puede calcular, S(1) = g(P (0)) y D(1) se puede
calcular puesto que D(1) = S(1) entonces D(1) = f (P (1)) = S(1) y as P (1) podra obtenerse de la ecuaci
on f (P (1)) = S(1). Este proceso podra repetirse para obtener P2 y P3 . . . ,
etc. Se estudiara la sucesion P0 , P1 , P2 , . . . .
Si asumimos que f y g son polinomios de primer grado, es decir que la gr
afica de la demanda y la oferta en funcion de P son rectas; obtendremos una ecuaci
on de diferencias lineal
de primer orden para P, ya que
Dk = mPk + b,

m>0yb>0

Sk+1 = mP
b k+C

m
b > 0 y C > 0.

La pendiente de la demanda se toma negativa ya que el incremento de una unidad en el

precio produce un decrecimiento igual a uno en la demanda y un incremento de una unidad
en el precio produce un incremento m
b en la oferta.
Por el supesto iii) S(k + 1) = D(k + 1) entonces

(2.14)

Despejando Pk+1 se obtiene

mPk + b = mP
b k+1 + C.

(2.15)

Pk+1 =

Haciendo A =

m
bC
Pk +
.
m
b
m
b

bC
m
yB=
, se obtiene que
m
b
m
b
Pk+1 = APk + .

A < 0 la sucesion {Pk } oscila.

41

El siguiente ejemplo mostraremos otra forma de obtener la solucion general de una ecuacion
en diferencia lineal de primer orden
Ejemplo 2.5.2. Dada la ecuaci
on lineal de primer orden
(2.16)

yk+1 = a(k)yk+1 + b(k) a(k) 6= 0.

Supongamos que yk es la soluci

on general, la cual se escribe de la siguiente forma
(2.17)

y k = uk v k ,

ya que el producto de dos funciones u y v puede ser determinado. Entonces

yk+1 ak yk = uk+1 vk+1 ak uk vk =
= uk+1 vk+1 + uk+1 vk uk+1 vk ak uk vk =
= uk+1 (vk + vk ) ak uk vk =
= vk (uk+1 ak uk ) + uk+1 vk = b(k).
Ahora tomando
(2.18)

uk+1 ak uk = 0

y
(2.19)

uk+1 vk = b(k),

entonces yk es una soluci

on de la ecuaci
on (2.16). La funcion u es una soluci
on de la
ecuaci
on homogenea (2.18), la cual tiene la siguiente soluci
on
u0 = C1 ,

uk = C1

k1
Y
j=0

a(j) k = 1, 2, 3,

y C1 R.

Despejando vk de (2.19) se obtiene

vk =

b(k)
,
uk+1

ya que uk 6= 0, por que a(j) 6= 0 para todo j.

Aplicando suma indefinida obtenemos lo siguiente
(2.20)

vk = 1

b(k)
+ C2 .
uk+1

Ahora hallemos nuevamente la soluci

on general de la ecuaci
on lineal en diferencias de coeficientes constantes
yk+1 = ayk + b,

42

donde a(k) = a y B(k) = b. Falta hallar v se usa la ecuaci

on (2.20) y se obtiene lo siguiente
vk =

b
ak+1

b
= 1
a

 k
1
,
a

aplicando la ecuaci
on (1.36) a

Entonces

 k
 k
1
a

1
,
1

= 1a a

k,

 k
b
1
,
vk = 1 a a

bk + C ,
2

si a = 1.

si a 6= 1
si a = 1.

Luego la soluci
on de yk+1 = ayk + b, esta dada por

C C a k + b ,
1 2
y k = uk v k =
1a
bk + C C ,
1

si a 6= 1

si a 6= 1
si a = 1.

Haciendo A = C1 C2 se obtiene que

Aak + b ,
y k = uk v k =
1a
A + bk,

si a 6= 1
si a = 1.

Se obtiene la misma soluci

on del Teorema 2.5.1.

2.6.

El M
etodo de Coeficientes Indeterminados

A continuacion presentaremos otro metodo que nos permitira hallar la solucion particular
de la ecuacion lineal en diferencia no homogenea
(2.21)

yk+1 = yk + b(k).

El metodo de coeficiente indeterminado consiste en hacer una suposicion de la forma general

de yp (k) cuando la funcion b(k) es bastante sencilla. A continuacion se presentaran algunos
casos:
1) b(k) es un polinomio
Comencemos con el caso donde b(k) es un polinomio de primer grado para despues
generalizar para polinomios de cualquier grado. Sea b(k) = b0 + b1 k, entonces la ecuacion
(2.2) tiene la forma
(2.22)

yk+1 ayk = b0 + b1 k,
43

b 6= 0.

Dado que el segundo miembro de la ecuacion (2.22) es un polinomio de primer grado,

entonces yk+1 ayk debe ser tambien un polinomio de primer grado. Luego se puede
conjeturar la existencia de una solucion particular que sea tambien un polinomio de
primer grado:
(2.23)

yp (k) = 0 + 1 k,

reemplazando (2.23) en (2.22) se tiene

0 + 1 (k + 1) a(0 + 1 k) = b0 + b1 k.
Agrupando terminos semejante se obtiene
(0 + 1 a0 ) + (1 a1 )k = b0 + b1 k.
Igualando los coeficientes de potencias iguales de k:
(
0 + 1 a0 = b0 ,
1 a1 = b1 ,
se producen los siguientes valores
b0 b1
b0 b1
=
,
2
(1 a)(1 a)
(1 a)(1 a)3
b1
,
1 =
1a
0 =

siempre que a 6= 1.
Si a = 1 la conjetura anterior ha fallado y no existe ninguna solucion de la forma
propuesta. La razon esta en que, en este caso, yk+1 yk es un polinomio de grado cero
y el segundo miembro de (2.22) es de primer grado. Lo cual se propone para yp (k) un
polinomio de grado dos. Pero como a = 1, la solucion general de la ecuacion homogenea
es yh (k) = C, el termino constante de yp (k) es absorbido por la solucion de la ecuacion
homogenea. Propongamos la solucion particular de la forma
yp (k) = 0 k + 1 k 2 .

(2.24)

Sustituyendo (2.24) en (2.22) con a = 1, resulta

0 (k + 1) + 1 (k + 1)2 (0 + 1 k 2 ) = b0 + b1 k.
Simplificando se obtiene
(0 + 1 ) + 21 k = b0 + b1 k.
Igualando coeficientes

0 1 = b0
21 = b1 ,
44

y as
0 = b0
1 =

b1
,
2

b1
.
2

Luego la ecuacion (2.22) tiene solucion particular de la forma:

b
b1
b1

0
k, si a 6= 1
+
2
1a
yp (k) = 1 a (1 a)

(b0 b1 )k + b1 k 2 ,
si a = 1
2
2
Lo anterior se puede generalizar para ecuaciones con termino b(k) = b0 + b1 k + +
bn k n , bn 6= 0, la ecuacion
(2.25)

yn+1 ayk = b0 + b1 k + + bn k n

cuya solucion general es de la forma

(2.26)

(
0 + 1 k + + n k n ,
si a 6= 1
yp (k) =
n
k(0 + 1 k + + n k ), si a = 1

Los coeficientes 0 + 1 k + + n k n se determinan sustituyendo (2.26) en la ecuacion

(2.28) e igualando el coeficiente que resulta en el primer miembro para cada potencia de
k con el correspondiente bk del segundo miembro.
Ejemplo 2.6.1. Resolvamos, por el metodo de coeficientes indeterminados, la siguiente
ecuaci
on
(2.27)

yk = yk1 + k 2 .

Se puede comprobar que la soluci

on que satisface y(0) = 0 es
1
1
1
yk = k 3 + k 2 + k.
3
2
6
La soluci
on general de la ecuaci
on homogenea asociada es
yh (k) = C.
Dado que b(k) es un polinomio de grado 2, y a = 1, lo adecuado es buscar una soluci
on
de la forma
(2.28)

yp (k) = k(0 + 1 k + 2 k 2 ) = (0 k + 1 k 2 + 2 k 3 ).

Sustituyendo (2.28) en la ecuaci

on(2.27)
0 k + 1 k 2 + 2 k 2 [0 (k 1) + 1 (k 1)2 + 2 (k 1)2 ] = k 2 .
45

resulta,

0 1 + 2 = 0
21 32 = 0

32 = 1.

(2.29)

Resolviendo el sistema de ecuaciones se obtiene que

1
0 = ,
6

1 =

1
2

1
y 2 = .
3

Entonces la soluci
on particular es
1
1
1
yp (k) = k + k 2 + k 3 .
6
2
3

(2.30)

Luego, la soluci
on general es de la ecuaci
on completa (2.27) esta dada por
1
1
1
y(k) = yh (k) + yp (k) = C + k + k 2 + k 3 .
6
2
3
Ahora reemplazando las condiciones iniciales obtenemos
1
1
1
y(0) = C + (0) + (0)2 + (0)3 = 0
6
2
3
se obtiene la constante C = 0, luego se comprueba que la soluci
on que satisface y(0) = 0
es
1
1
1
yk = k + k 2 + k 3 .
6
2
3

(2.31)

Tomando la ecuaci
on (2.27) y haciendo iteraciones con y0 = 0 tenemos

(2.32)

y1 = y0 + 12 = 12 ,
y2 = y1 + 22 = 12 + 22 ,
y3 = y2 + 32 = 12 + 22 + 32 ,
..
.
yn = yn1 + n2 = 12 + 22 + + n2 ,

k = 1.
k = 2.
k = 3.
k = n.

Lo cual se puede demostrar por induccion.

Ahora haciendo k = n en (2.31) tenemos yn = 16 n + 21 n2 + 13 n3 y reemplazando en (2.32)
se obtiene
1
1
1
1
3
1
n
1
n(n + 1)(2n + 1)
n + n2 + n3 = n(1 + n + n2 ) = (n + )(n + 1) =
=
6
2
3
3
2
2
3
2
6
=12 + + n2 ,
lo cual implica que
n(n + 1)(2n + 1)
= 12 + 22 + + n2
6
46

es la formula conocida para la suma de cuadrados de numeros enteros. De esta manera

se puede establecer la validez para otras formulas conocidas como
1
13 + 23 + + n3 = n2 (n + 1)2 = (1 + 2 + + n)2 ,
2
1
14 + 24 + + n4 = n(n + 1)(2n + 1)(3n2 + 3n 1).
30
2) Sea b(k) = (0 + 1 k + + k n n ) k .
Supongamos ahora que el termino b(k) es el producto de una potencia por un polinomio,
o sea de la forma
b(k) = Pn (k) k , 6= 0
la ecuacion
(2.33)

yk+1 = yk + Pn (k) k = yk + (0 + 1 k + + k n n ) k .

Dividiendo por k+1 en la ecuacion (2.33), obtenemos

yk+1
yk
Pn (k) k
1 yk
Pn (k)
.
=
a
+
=a
+
k+1
k+1
k+1
k

yk
Haciendo un cambio de variable en (2.34), donde uk = k , resulta

(2.34)

(2.35)

uk+1 =

1
a
uk + Pn (k).

Las ecuaciones (2.21) y (2.35) son equivalentes en el sentido de que a cada solucion yk
de (2.21) le corresponde una y solo una solucion u(k) de (2.35). Entonces la ecuacion
(2.33) tiene una solucion particular de la forma
(
k (0 + 1 k + + n k n ),
si 6= a
(2.36)
yp (k) =
k k(0 + 1 k + + n k n ), si = a.
De nuevo los coeficientes 0 , 1 , n se determinan sustituyendo (2.36) en la ecuacion
(2.33) e igualando el coeficiente que resulta en el primer miembro para cada potencia de
k con el correspondiente bk del segundo miembro.
El siguiente ejemplo nos va permitir pasar de una ecuacion lineal en diferencia cuando a(k) y b(k) no son constantes a una ecuacion de la forma (2.21).
Ejemplo 2.6.2. Dada la constante real 6= 1, resolver la ecuaci
on
(2.37)
Multiplicando
(2.38)

yk+1 =

k+2
yk + (k + 2)k
k+1

1
ambos lados la ecuaci
on (2.37), se obtiene
k+2
1
yk+1
=
yk + k .
(k + 2)
k+1
47

Haciendo u(k) =

yk
y sustituyendo en (2.38), se obtiene que
k+1
uk+1 = uk + k .

(2.39)

Como 6= 1 y a = 1 en la ecuaci
on (2.39) entonces 6= a, se toma una soluci
on de la
forma
up (k) = 0 k .
Reemplazando en (2.39) resulta que
0 k+1 0 k = k = 0 k 0 k = k ,
entonces
0 0 = 1,

y as 0 =

1
y la soluci
on particular de (2.39) es
( 1)
up (k) =

1
k .
1

Por lo tanto la soluci

on particular de (2.37) es
yp (k) =
y su soluci
on general es
y(k) = C(k + 1) +

1
(k + 1)k
1
1
(k + 1)k ,
1

C R.

3) Sea b(k) = a0 cosk + b0 sen k. Nuevamente la idea consiste en buscar una solucion
particular del mismo tipo, es decir,
(2.40)

yp (k) = A0 cosk + B0 sen k.

Sustituyendo en el primer termino de la ecuacion (2.21)

yk+1 ayk = A0 cos(k + ) + B0 sen(k + ) aA0 cosk aB0 sen k
= A0 (cos k cos sen k sen ) + B0 (sen k cos + cos k sen )
aA0 cos k aB0 sen
= (A0 cos + B0 sen aA0 ) cos k + (A0 sen + B0 cos aB0 ) sen k
e igualando al segundo termino resulta
(A0 cos +B0 sen aA0 ) cos k+(A0 sen +B0 cos aB0 ) sen k = a0 cosk+b0 sen k.
Para que estas funciones trigonometricas coincidan para todo k > 0, han de ser iguales
los coeficientes de cos k y sen k; as, se deduce que
(
A0 cos + B0 sen aA0 = a0
A0 sen + B0 cos aB0 = b0 .
48

Este sistema si tiene solucion, determina los valores de A0 y B0 . En el caso de que no

exista solucion del sistema, habra que buscar una solucion particular de la forma
yp (k) = A0 k cos k + B0 k sen k,
esto sucede cuando alguno de los sumandos de la expresion (2.40) es solucion de ecuacion
homogenea.
Ejemplo 2.6.3. Dada la ecuaci
on lineal en diferencia de primer orden
k
,
2
hallemos la soluci
on homogenea asociada. Como a(k) = 1 se tiene que
yk+1 = yk + sen

yh (k) = y0

k1
Y

1 = y0 .

j=0

Supongamos que la soluci

on particular tiene la forma
k
k
+ B0 sen
,
2
2
entonces sustituyendo en la ecuaci
on obtenemos
yp (k) = A0 cos

A0 cos

(k + 1)
(k + 1)
k
k
k
+ B0 sen
A0 cos
B0 sen
= sen
.
2
2
2
2
2

Dado que
k
(k + 1)
= sen
,
2
2
k
(k + 1)
= cos
,
sen
2
2
cos

resulta
A0 sen

k
k
k
k
k
+ B0 cos
A0 cos
B0 sen
= sen
.
2
2
2
2
2

Simplificando tenemos
(A0 B0 ) sen
de lo cual se deduce

k
k
k
+ (B0 A0 ) cos
= sen
,
2
2
2

(
A0 B0 = 1
B0 A0 = 0.

on particular es
As A0 = 12 y B0 = 32 , por lo que una soluci
1
k 3
k
yp (k) = cos
+ sen
,
2
2
2
2
y la soluci
on general de la ecuaci
on completa viene dada por:
yk = y0

k 3
k
1
cos
+ sen
.
2
2
2
2
49

2.7.

Soluci
on Geom
etrica

Ya vimos en secciones anteriores como resolver una ecuacion lineal en diferencia de primer
orden, para responder preguntas especificas de estas soluciones. Este metodo grafico o geometrico nos permite obtener informacion cualitativas de la solucion general de ecuaciones
difciles de resolver por metodos especficos.
Este nuevo metodo esta basado en utilizar la grafica de la funcion f (x) para generar el
siguiente punto que toma la trayectoria de la solucion y luego proyectarlo con ayuda de la
bisectriz y = x, en el eje horizontal; aplicando reiteradamente este procedimiento se generan
los sucesivos puntos de la trayectoria. Veamoslo con mas detenimiento el siguiente ejemplo.
Ejemplo 2.7.1. Retomemos nuevamente la ecuaci
on de interes compuesto con un ahorro
A al final de cad periodo, entonces
Ck+1 = (1 + i)Ck + A,
la cual se puede representar por la recta o ecuaci
on
y = f (x) = (1 + i)x + A.
Comencemos con un punto inicial x0 = 0, el siguiente punto de la trayectoria es x1 =
f (x0 ) = A, est es, x1 = A es el punto del eje vertical tal que (x0 , x1 ) = (0.A) pertenece la
gr
afica y = f (x). Proyectando dicho punto en el eje horizontal con ayuda de la recta y = x,
podemos aplicar de nuevo la funcion f (x) y obtener x2 = f (x1 ) = (1 + i)A + A = A(2 + i))
y proyectarlo nuevamente por medio de la ecuaci
on x = y. Hacemos lo mismo para obtener
x3 y as sucesivamente.
Claramente se ve que el valor siguiente crece de forma proporcional, por lo tanto su soluci
on
gr
afica es la siguiente:
y = Ck+1
Ck+1 = (1 + i)Ck + A

y=x

x0

x1

x2

x = Ck

Ente caso como la pendiente de la ecuaci

on es positiva la sucesion que se forma es creciente,
en caso contrario tendramos una sucesion decreciente y las flechas de la gr
afica iran en
sentido contrario.
50

Ejemplo 2.7.2. Aqu vamos retomar el Ejemplo 2.5.1 relacionado con la oferta y la demanda. Veamos sus posibles soluciones gr
aficas. Partiremos desde la ecuaci
on (2.14) y la
separemos en dos ecuaciones equivalentes
y = mx + b,

m>0yb>0

y
y = mx
b + C,

m
b > 0 y C > 0,

y se precede de la misma forma que le ejemplo anterior pero con la u

nica diferencia que no se
proyecta sobre la ecuaci
on y = x, de esta forma se obtiene la gr
afica referente P (k)(precios)
y demanda o oferta. A continuaci
on se presentaran los diferentes comportamientos que
puede tomar A en la ecuaci
on (2.15).
i) Si 1 < A < 0 la recta P (k + 1) = AP (k) + b es decreciente y genera una sucesion
puntos que converge al punto de interseccion entre la dos rectas.
D(k) 0 S(k + 1)

P1

P3P5

P4 P2

P0

P (k)

ii) Si A = 1, todas las soluciones no constates son peri

odicas de periodo dos.
P (k + 1)

P2

P0

iii) Si A < 1 se generan oscilaciones no acotadas

51

P (k)

P (k + 1)

P1

2.8.

P0

P2

P (k)

Soluci
on de Ecuaciones no Lineales

En este seccion vamos considerarar algunas ecuaciones que no cumplen la Definicion 2.1.1
de ecuaciones lineales de primer orden, pero con un cambio de variable la podemos llevar
a una ecuacion de ese tipo. Otra forma de utilizar el metodo que vamos desarrollar, es
llevar algunas ecuaciones con coeficientes variables, es reducirla o convertirla a coeficientes
constantes.
Ejemplo 2.8.1. Sea la ecuaci
on yk yk+1 = kyk yk+1 . Se observa que yk = 0 es una soluci
on
de la ecuaci
on. En este caso, esta soluci
on es trivial. No obstante, podemos encontrar la
soluciones no triviales de la forma siguiente:
Dividamos le ecuaci
on por yk yk+1 , con la cual
1
yk+1

1
= k.
yk

on de primer orden
Haciendo un cambio de variable uk = y1k , tenemos la siguiente ecuaci
lineal
uk+1 uk = k,
cuya soluci
on general es uk = C + k, volviendo a la variable original tenemos yk =
As, tenemos todas las soluciones de la ecuaci
on lineal, incluyendo la trivial.

1
.
C+k

El siguiente ejemplo nos permite llevar una ecuacion de coeficiente variables a una de coeficientes constantes, dichas soluciones se dejamos al lector.
2
Ejemplo 2.8.2. i) considere la ecuaci
on yk+1
yk2 = 1. Haciendo uk = yk2 nos queda la
siguiente ecuaci
on lineal de primer orden

uk+1 uk = 1.

ii) Considerando la ecuaci

on yk+1 = (1+ 3 yk )3 . Haciendo un cambio de variable uk = 3 yk
tenemos la siguiente ecuaci
on
q
p
3
3
uk+1 = (1 + uk ) 3 u3k+1 = 3 (1 + uk )3 uk+1 = 1 + uk
52

luego tenemos la siguiente ecuaci

on lineal primer orden
uk+1 uk = 1.
Nota: Este metodo de cambio de variable se puede aplicar para ecuaciones de orden superior
no lineales y llevarlas a ecuaciones lineales de orden superior que se va desarrollar en el
siguiente capitulo.

53

CAPITULO 3

LINEAL EN DIFERENCIAS
ECUACION
DE ORDEN SUPERIOR.

3.1.

Introducci
on

En el capitulo anterior se presto atencion a las ecuaciones en diferencias de primer orden,

haciendo enfasis en las ecuaciones lineales, para las cuales se estudio la estructura de las
soluciones, algunos metodos para calcular soluciones, aspectos geometricos de estas soluciones y se presentaron algunas situaciones que conducan a este tipo de ecuaciones.
En el presente capitulo se hara un estudio similar para ecuaciones lineales en diferencias de
orden mayor o igual que dos. Para efectos de calculo tomaremos ecuaciones de orden dos y
tres.
Comenzaremos con la sucesion de Fibonacci, quec se puede plantear de la siguiente forma
yn+2 = yn+1 + yn . La cual es una ecuacion en diferencia finita de orden dos.
Damos a continuacion la definicion de una ecuacion lineal en diferencia de orden n.
Definici
on 3.1.1. Una ecuaci
on en diferencia de la forma
(3.1)

yk+n + a1 (k)yk+n1 + + an1 (k)yk+1 + an (k)yk = b(k).

donde ai , b; para i = 0, 1, 2, . . . son funciones que depende k N, donde an (k) 6= 0, es una

ecuaci
on lineal en diferencias de orden n.
En el caso especial en que b(k) = 0, para todo k, la ecuaci
on se denomina homogenea, y
cuando las funciones ai son constantes, la ecuaci
on se llama lineal con coeficientes constantes.
54

Ejemplo 3.1.1. Un ejemplo de una ecuaci

on de orden tres es
yk+3 + 3yk+2 + 9yk = cos k.
Aqu
a1 (k) = 3,

3.2.

a2 (k) = 9,

b(k) = cos k.

Existencia de una Soluci

on

Estableceremos condiciones bajo las cuales una ecuacion en diferencias de orden mayor que
uno tiene solucion, El teorema de existencia es importante, ya que por ejemplo
yk+2 kyk+1 yk = 0 k = 0, 1, 2,
no tiene solucion. Si y0 = 0 y y2 = 1, pues que y2 y0 = 1 0 = 1. Por otro lado si
y0 = y2 = 1; entonces y2 y0 = 0; y como para k = 1 y3 y2 y1 = y3 y2 y1 ; para
satisfacer la ecuacion se tendra que y3 dependa de y1 y as la ecuacion tendra infinitas
soluciones.
Teorema 3.2.1 (Existencia y unicidad). Dada la ecuaci
on lineal
yk+n + a1 (k)yk+n1 + a2 (k)yk+n2 + + an1 (k)yk+1 + an (k)yk = b(k),

para k = 0, 1, 2,

y las constantes 1 , 2 , . . . , n1 , existe una u

nica soluci
on tal que y0 = 0 , y1 = 1 , . . . , yn1 =
n1 .
Demostraci
on. Como se conocen
y0 = 0 , y1 = 1 , y2 = 2 , . . . , yn1 = n1 ,
al reemplazarlo en la ecuacion, haciendo k = 0, se obtiene yn , ya que
yn + a0 (0)yn1 + a2 (0)yn2 + an1 (0)y1 + an (k)y0 = b(0).
Ahora suponganse que se conocen los terminos de la sucesion hasta yj para alg
un j n.
Mostraremos que el valor de yj+1 puede construirse.
Reemplazando en la ecuacion para k = j + 1 n se obtiene
yj+1 = b(j + 1 n) a1 (j + 1 n)yj +an (j + 1 n)yj+1n .
Como por hipotesis los terminos de la derecha son conocidos se conoce yj+1 . De esta forma
se concluye inductivamente que la ecuacion tiene solucion. Ademas puede observase que
cualquier otra sucesion con valores iniciales iguales coincide con la solucion encontrada. 

55

3.3.

Propiedades de la Ecuaci
on Homog
enea

La ecuacion homogenea juega un papel central en la b

usqueda de la solucion de la ecuacion
no homogenea. Dos hechos indican la importancia de la ecuacion homogenea en la teora de
ecuaciones lineales. El primero es que la diferencia entre dos soluciones de la ecuacion no
homogenea es una solucion de la homogenea asociada y en segundo lugar que si dos solucion
de la ecuacion homogenea se suman sigue siendo una solucion de la ecuacion homogenea.
Por consiguiente se estudian aspectos muy importantes de la ecuacion homogenea.
Teorema 3.3.1. Suponganse que las sucesiones, {wk } y {zk } son soluciones de la ecuaci
on
lineal homogenea de orden n. Entonces para n
umeros reales 1 y 2 , 1 wk +2 zk es soluci
on
de la homogenea
Demostraci
on. Tomando yk = 1 wk + 2 zk y reemplazando en la ecuacion lineal homogenea
yk+n + a1 (k)yk+n1 + a2 (k)yk+n2 + + an1 (k)yk+1 + an (k)yk = 0 para k = 0, 1, 2,
tenemos
1 wk+n + 2 zk+n + a1 (k)[1 wk+n1 + 2 zk+n1 ] + + an (k)[1 wk + 2 zk ] =
1 [wk+n + a1 (k)wk+n1 + + an (k)wk ] + 2 [zk+n + a2 (k)zk+n1 + + a2 zk ].
Como {wk } y {zk } son soluciones de la ecuacion homogenea, luego se obtiene que yk =
1 wk + 2 zk es una solucion de la ecuacion lineal homogenea.

El teorema anterior nos muestra que puede existir un conjunto muy grande de soluciones
para una ecuacion en diferencia de orden n. Esto nos lleva a preguntarnos si es posible encontrar un conjunto especial de soluciones de tal forma que pueda ser usada para representar
todas las soluciones. La repuesta esta a esta pregunta nos lleva a formular la definicion de
soluciones independientes.
Definici
on 3.3.1. Dada la ecuaci
on homogenea de orden n, las soluciones {yk } y {ybk } se
llaman independientes si dado que yk + ybk = 0 para todo k entonces = = 0. De lo
contrario son linealmente dependientes.
Ejemplo 3.3.1.

yk+2 yk+1 yk = 0
yk = 1, ybk = k.

La soluci
on {yk } y {ybk } son independientes ya que si yk + ybk = 0 para todo k, entonces
y1 + yb1 = 0,
y2 + yb2 = 0,

entonces + = 0.
entonces + 2 = 0.

De lo ultimo se tiene = = 2 y as = 0 y = 0.
56

El siguiente teorema muestra que toda ecuacion lineal homogenea de orden n tiene n soluciones independientes.
Teorema 3.3.2. Toda ecuaci
on lineal homogenea de orden n tiene n soluciones independientes.
Demostraci
on. Sea la ecuacion homogenea de orden n
yk+n + a1 (k)yk+n1 + + an1 (k)yk+1 + an (k)yk = 0.
Por el Teorema de existencia y unicidad se tiene lo siguiente:
1) Si y 1 (k) es la u
nica solucion tal que
1
y01 = 1, y11 = 0, . . . , yn1
= 0,

2) si y 2 (k) es la u
nica solucion tal que
2
= 0,
y02 = 0, y12 = 1, . . . , yn1

..
.
n) si y n (k) es u
nica solucion tal que
n
y0n = 0, y1n = 0, . . . , yn1
= 1.

Las soluciones {y 1 (k)}, {y 2 (k)}, , {y n (k)} son linealmente independientes, ya que si

C1 y 1 (k) + C2 y 2 (k) + + Cn y n (k) = 0,
tomando cada problema de condicion inicial y haciendo k = 0, 1, 2, . . . , n 1 obtenemos
1 y 1 (0) + 2 y 2 (0) + + n y n (0) = 0 entonces 1 y 1 (0) = 0. Por lo tanto 1 = 0.
1 y 2 (1) + 2 y 2 (1) + + n y n (1) = 0 entonces 2 y 2 (1) = 0. Por lo tanto 2 = 0.
As sucesivamente resulta que 1 = 2 = = n = 0.

A continuacion mostraremos que toda solucion de la ecuacion homogenea se puede escribir
en terminos soluciones independientes.
Teorema 3.3.3. Sea {zk } una soluci
on de la ecuaci
on lineal homogenea en diferencias
finitas de orden n, y y1 (k), y2 (k), . . . , yn (k) las soluciones obtenidas en el teorema precedente.
Entonces existen constantes 1 , 2 , . . . n , tales que zk = 1 y1 (k) + 2 y2 (k) + . . . + n yn (k).
Demostraci
on. Sea y1 (k), y2 (k), . . . , yn (k) soluciones de la ecuacion lineal homogenea por
el Teorema 3.3.1 existen escalares C1 , C2 , . . . , Cn tal que
yk = C1 y1 (k) + C2 y2 (k) + . . . + Cn yn (k),
es solucion de ecuacion homogenea. Como {z} es otra solucion de la ecuacion homogenea
tenemos el siguiente conjunto de n + 1 de soluciones y1 (k), y2 (k), . . . , yn (k), z(k) linealmente
57

dependiente por el Teorema 3.3.2 que solo se admiten n soluciones linealmente independiente. Aplicando la definicion de dependencia lineal, existen escalares C1 , C2 , . . . , Cn , Cn+1 no
todos ceros tal que
C1 y1 (k) + C2 y2 (k) + . . . + Cn yn (k) + Cn+1 z(k) = 0.
Pos ser y1 (k), y2 (k), . . . , yn (k) linealmente independientes Cn+1 6= 0, entonces
z(k) =
Tomando i =

C1
C2
Cn
y1 (k) +
y2 (k) + . . . +
yn (k).
Cn+1
Cn+1
Cn+1

yi (k)
, i {1, 2, . . . n} entonces
Cn+1
z(k) = 1 y1 (k) + 1 y2 (k) + . . . + n yn (k),

lo cual es una solucion particular de la ecuacion homogenea lineal.


Tenemos seg
un los teoremas precedentes que el conjunto solucion de la ecuacion lineal en
diferencias de orden n es un espacio vectorial de dimension n. y por consiguiente, cualquier
solucion no solamente puede expresarse como combinacion lineal de las soluciones basicas
dadas, sino que puede expresarse en terminos de cualquier conjunto independiente de n
soluciones de la ecuacion.
Ejemplo 3.3.2. Estudiemos la ecuaci
on yk+2 3yk+1 + 2yk = 0. Puede verificarse que
yk = 1 y zk = 2k son soluciones ya que
13+2=0 y

2k+2 3 2k+1 + 2 2k = 2k [22 3 2 + 2] = 0.

Estas soluciones son independientes ya que si tenemos el siguiente sistema de ecuaciones

(
(
C1 yk + C2 zk ,
C1 + C2 2k ,

C1 yk+1 + C2 zk+1 .
C1 + C2 2k+1 .
Despejando C1 de la primera y reemplazando en la segunda se tiene
C2 (2k + 2k+1 ) = 0.
Como 2k + 2k+1 6= 0 para k = 1, 2, 3, . . . entonces C2 = 0. Reemplazando el valor de C2 en
la primera se obtiene
C1 + 0 2k = 0.
Luego C1 = 0.
Entonces toda soluci
on de esta ecuaci
on, es de la forma wk = 1 yk + 2 2k para alguna
pareja 1 , 2 de n
umeros reales.

58

3.4.

Ecuaciones Homog
eneas con Coeficientes
Constantes

En la seccion 3.3 vimos que una solucion de la ecuacion lineal en diferencias finitas de orden
n es una combinacion lineal de n soluciones particulares linealmente independientes, pero
no dijimos como encontrar tal solucion. Explicaremos ahora como encontrar de manera
explcita y directa, n soluciones linealmente independientes de una ecuacion lineal en diferencias finitas de orden n con coeficientes constantes. A continuacion presentaremos algunos
ejemplos de ecuaciones homogeneas con coeficientes constantes.
Ejemplo 3.4.1. Considere las ecuaciones lineales en diferencia
yk+1 yk = 0,
yk+2 3yk+1 + 2yk = 0,
yk+3 yk = 0.
Son ecuaciones con coeficientes constantes de orden 1, 2 y 3 respectivamente.
Busquemos una solucion sencilla de la ecuacion homogenea y observemos lo que sucede al
reemplazar yk = mk en la ecuacion lineal homogenea
(3.2)

yk+n + a1 yk+n1 + a2 yk+n2 + + an1 yk+1 + (k)yk = 0,

con a1 , a2 , . . . , an constantes.
mk+n + a1 mk+n1 + a2 mk+n2 + + an1 mk+1 + an mk = 0.

Factorizando tenemos

mk [mn + a1 mn1 + a2 mn2 + + an1 m + an ] = 0,

entonces como mk nunca es cero, yk = mk sera una solucion de la ecuacion (3.2) cuando m
sea una raz de la ecuacion
mn + a1 mn1 + a2 mn2 + + an1 m + an = 0.
Esto nos lleva plantear la siguiente definicion.
Definici
on 3.4.1. Llamaremos ecuaci
on o polinomio caracterstico asociado a la ecuaci
on
(3.2) a
p(m) = mn + a1 mn1 + + an .
Ejemplo 3.4.2. Sea la ecuaci
on lineal homogenea con coeficientes constante yk+2 3yk+1 +
2yk = 0. Su polinomio caracterstico es :
p(m) = m2 3m + 2 = 0.
Este polinomio se puede descomponer en dos factores
m2 3m + 2 = (m 2)(m 1),
entonces sus races son: m = 2 y m = 1.
59

Ejemplo 3.4.3. Sea la ecuaci

on lineal homogenea con coeficientes constantes yk+2 +yk = 0.
Su polinomio caracterstico es el siguiente:
m2 + 1 = 0 o

m2 = 1,

lo que implica que tiene dos soluciones complejas diferentes m = i y m = i.

Entonces, nuestro problema se reduce a la solucion de una ecuacion puramente algebraica.
Por el teorema fundamental del algebra, todo polinomio de grado n tiene n races no todas
necesariamente distintas ni reales como se vio en los dos ejemplos anteriores, por tal la
solucion de una ecuacion en diferencia presenta diferentes casos, que a continuacion desarrollaremos.
Teorema 3.4.1. Sea la ecuaci
on lineal en diferencia, homogenea con coeficientes constantes
yk+n + a1 yk+n1 + a2 yk+n2 + + an1 yk+1 + (k)yk = 0,
y su polinomio asociado p(m) = mn + a1 mn1 + + an .
a) Si m0 es raz de p(m), entonces, mk0 es soluci
on de la ecuaci
on, homogenea.
b) Si a + bi es raz de p(m), entonces rk sen k, rk cos k son soluciones de la ecuaci
on
homogenea siendo r el modulo y el argumento de a + bi.
Demostraci
on.
a) Sean
yk+n + a1 yk+n1 + + an1 yk+2 + an yk
y
p(m0 ) = mn0 + a1 mn1
+ + an .
0
Sea yk = mk0 entonces
m0k+n + a1 m0k+n1 + + an1 mk+2
+ an mk0 = mk0 (mk0 + a1 mn1
+ + an1 m20 + an ) =
0
0
= mk0 p(m0 ) = 0,

por lo tanto yk = mk0 es solucion de la ecuacion homogenea

b) Tomando la parte a) tenemos que si a + bi es solucion del polinomio caracterstico p(m)
entonces (a + bi)k es solucion de la ecuacion homogenea en diferencias.

b
Sea r = a2 + b2 que es el modulo de a + bi y = arc tg el argumento.
a
Por el Teorema de Demoivre se tiene que
(a + bi) = r(cos + i sen ) y (a + bi)k = rk (cos k + i sen k) = rk cos k + irk sen k.
Como las races de un polinomio vienen en pares conjugados entonces a bi tambien es
solucion del polinomio caracterstico por lo tanto
(a bi)k = rk (cos k i sen k) = rk cos k irk sen k,
60

tambien es solucion de la ecuacion homogenea en diferencia.

Tomando
1
1
rk cos k = [rk cos k + irk sen k] + [rk cos k irk sen k]
2
2
1
1
= (a + bi)k + (a bi)k
2
2
y
1 k
1
[r cos k + irk sen k] [rk cos k irk sen k]
2i
2i
1
1
= (a + bi)k (a bi)k .
2i
2i

rk sen k =

Por lo tanto, rk sen k y rk cos k son soluciones de la ecuacion homogenea.

Del teorema anterior podemos concluir que, si se tienen n-races reales distintas de la
ecuacion caracterstica, la solucion ecuacion homogenea tiene la siguiente forma
yk = C1 mk1 + C2 mk2 + + Cn mkn ,
ya que el conjunto de soluciones es linealmente independiente.
De igual forma para races complejas se tiene
yk = C1 (a + ib)k + C2 (a ib)k =rk cos k + irk sen k + rk cos k irk sen k =
=(C1 + C2 )rk sen k + (C1 C2 )i cos k =

=c1 rk sen k + c2 cos k,

donde c1 = (C1 + C2 ) y c2 = (C1 C2 )i.

Ejemplo 3.4.4. Retomando los Ejemplos 3.4.2 y 3.4.3 se tiene las siguientes soluciones:
Para el Ejemplo 3.4.2
yk = C1 2k + C2 ,
y para el Ejemplo 3.4.3

yk = C3 cos k + C4 sen k.
2
2
Teorema 3.4.2. Sea la ecuaci
on en diferencias homogenea
a) Si m0 es raz con multiplicidad p, entonces
mk0 , kmk0 , k 2 mk0 , . . . k p1 mk0 ,
son soluciones de la ecuaci
on de diferencias.

61

b) Si a + bi, a bi son complejos conjugados con multiplicidad p, entonces

rk cos k,

krk cos k, . . . ,

k p1 rk cos k,

rk sen k,

krk sen k, . . . ,

k p1 rk sen k,

son soluciones de la ecuaci

on homogenea en diferencias.
Demostraci
on.
a) Utilizando el operador de traslacion podemos escribir la ecuacion
yk+n + a1 yk+n1 + + an1 yk+2 + an yk = 0,
como sigue
[E n + a1 E n1 + + an1 E + an ]yk = 0,

como m0 es una raz de multiplicidad p, lo anterior se puede expresar como

(E m0 )p yk = 0.

(3.3)

Tenemos que y1 = mk0 es una solucion y hagamos y(k) = uk y1 = uk mk0 , entonces

k+1
(E m0 )uk mk0 = E(uk mk0 ) uk mk+1
= uk+1 mk+1
uk mk+1
= mk+1
0
0
0
0 (uk+1 uk ) = m0 uk .

Entonces
p
(E m0 )p yk = mk+1
0 uk ,

para cualquier sucesion uk . entonces yk es solucion de la ecuacion (3.3) si y solo si

k uk = 0, luego
u(x) = 1 + k + k 2 + + k p1 ,

en particular mk0 , kmk0 , k 2 mk0 , . . . , k p1 mk0 .

b) Usando la parte a) tenemos

(a ib)k , k(a ib)k , k 2 (a ib)k , . . . k p1 (a ib)k .

Entonces
(r[cos i sen ])k , k(r[cos i sen ])k , k 2 (r[cos i sen )k , . . . k p1 (r[cos i sen ])k .
Realizando calculos obtenemos
rk cos k irk sen k, krk cos k ikrk sen k, k 2 rk cos ik 2 rk sen ,
, . . . k p1 rk cos k ik p1 rk sen k.

En particular se tiene
rk cos k,

krk cos k, . . . ,

k p1 rk cos k,

rk sen k,

krk sen k, . . . ,

k p1 rk sen k.

62


Ejemplo 3.4.5. Resolvamos un problema con condiciones iniciales
yk+2 + 2yk+1 + yk = 0.
y0 = 5,
y1 = 3.
Observemos que la ecuaci
on caracterstica
m2 + 2m + 1 = (m + 1)2 = 0,
tiene races iguales m1 = m2 = 1, es decir tiene multiplicidad p = 2. De aqu la soluci
on
general proporcionada por el Teorema 3.4.2 parte a) es
yk = c1 (1)k + c2 k(1)k .
Reemplazando las condiciones iniciales producen las siguientes ecuaciones
y0 = c1 (1)0 + c2 0 (1)0 = c1 = 5,
y1 = c1 (1)1 + c2 1 (1)1 = c1 c2 = 3,
lo cual implica que c1 = 5, y c2 = 2. As la soluci
on deseada del problema con condiciones
iniciales es
yk = 5(1)k 2k(1)k .
Ejemplo 3.4.6. La ecuaci
on caracterstica de la ecuaci
on lineal homogenea en diferencia
finita
yk+3 + yk+1 10yk = 0,
es la ecuaci
on cubica
m3 + r 10
la cual tiene races complejas 1 2i y una raz real 2. De aqu la soluci
on general proporcionada por el teorema 3.4.1 es
yk = c1 2k + (c2 cos k + c3 sen k).
Ejemplo 3.4.7. Sea Pn la poblaci
on n-enesima generaci
on de una especie de bacterias.
Suponga que la poblacion en la n-enesima generaci
on depende de las dos generaciones anteriores:
Pn+2 rPn+1 sPn = 0
donde r y s son medidas de la importancia relativa de las dos generaciones.
La ecuaci
on caracterstica es m2 rm s = 0 con races

r r2 + 4s
r + r2 + 4s
, m2 =
m1 =
2
2
63

Si r2 > 4, las races son reales diferentes y la soluci

on general es
Pn = c1 mk1 + c2 mk2 .
Si |m1 | < 1 y |m2 | < 1, entonces Pn tiende a cero cuando k tiende a infinito, pero |m1 | o
|m2 | es mayor que 1, entonces Pn tiende a infinito. Si r2 = 4s(s < 0), entonces la soluci
on
general es
 r k
 r k
Pn = c1
+ c2 k
,
2
2
en cuyo caso Pn tiende a cero si |r| < 2, pero si |r| 2, la soluci
on general es
k

donde = tan1
manera oscilatoria.

3.5.

Pn = (s) 2 (c1 cos k + +c2 k),

!
r
4s
. Cuando s < 1, la soluci
on Pn tiende a cero de una
r2
r

Soluci
on General de la Ecuaci
on Lineal no
Homog
enea con Coeficientes Constantes

En la secciones anteriores presentamos como hallar la solucion de una ecuacion homogenea

en diferencia de orden n. Igual que las ecuaciones de primer grado hay una relacion entre
esta solucion y la particular de la no homogenea. El Teorema 2.4.1 se puede generalizar para
ecuaciones de grado n donde su demostracion es similar a la demostrado, para ecuaciones
de orden uno.
A continuacion presentaremos dos formas de encontrar la solucion de una ecuaciones no
homogenea. El primer metodo que presentaremos es el de variacion de parametros y el
segundo de coeficientes indeterminados que es similar a expuesto para solucionar ecuaciones
de primer grado lo cual nos limitaremos a presentar algunos ejemplos para ecuaciones de
orden n.

3.5.1.

Variaci
on de Par
ametros

Supongamos la solucion de la ecuacion homogenea asociada a la ecuacion lineal en diferencia

de orden n
yk+n + a1 yk+n1 + a2 yk+n2 + + an1 yk+1 + an (k)yk = b(k),
tiene la siguiente forma
yH (k) = C1 y1 (k) + C2 y2 (k) + + Cn yn (k),
donde yi (k) i = 1, 2, 3, . . . , n son soluciones linealmente independientes.
A partir de esto, buscaremos una solucion completa de la siguiente forma:
(3.4)

yH (k) = C1 (k)y1 (k) + C2 (k)y2 (k) + + Cn (k)yn (k),

64

donde Ci (k) i = 1, 2, 3, . . . , n estan funcion de k.

Por otro lado
(3.5)

yH (k + 1) = C1 (k + 1)y1 (k + 1) + C2 (k + 1)y2 (k + 1) + + Cn (k + 1)yn (k + 1)

sumando y retando Ci (k + 1) a 3.5 tenemos

yH (k + 1) = [C1 (k) + C1 (k + 1) C1 (k + 1)]y1 (k + 1) + [C2 (k) + C2 (k + 1) C2 (k + 1)]y2 (k + 1)+
+ + [Cn (k) + Cn (k + 1) Cn (k + 1)]yn (k + 1) =
= [C1 (k) + C1 (k)]y1 (k + 1) + [C2 (k) + C2 (k)]y2 (k + 1)+
+ + [Cn (k) + Cn (k)]yn (k + 1)
= [C1 (k)y1 (k + 1) + C2 (k)y2 (k + 1) + + Cn (k)yn (k + 1)] + [C1 (k)y1 (k + 1)+
+ C2 (k)y2 (k + 1) + + Cn (k)yn (k + 1)].
la cual se puede escribir como
(3.6)

yH (k + 1) = C1 (k)y1 (k + 1) + C2 (k)y2 (k + 1) + + Cn (k)yn (k + 1)

siempre y cuando se tome Ci (k), i = 1, 2, 3, . . . , n de tal forma que

C1 (k)y1 (k + 1) + C2 (k)y2 (k + 1) + + Cn (k)yn (k + 1) = 0.
Analogamente, sumando y restando Ci (k) para i = 1, 2, . . . , n a la ecuacion 3.6 se puede
escribir de la siguiente forma
yH (k + 2) = C1 (k)y1 (k + 2) + C2 (k)y2 (k + 2) + + Cn (k)yn (k + 2),
siempre y cuando se tome Ci (k), i = 1, 2, 3, . . . , n de tal manera que
(3.7)

C1 (k)y1 (k + 2) + C2 (k)y2 (k + 2) + + Cn (k)yn (k + 2) = 0.

sucesivamente podemos escribir

yH (k + n 1) = C1 (k)y1 (k + n 1) + C2 (k)y2 (k + n 1) + + Cn (k)yn (k + n 1)
si elegimos Ci (k), i = 1, 2, 3, . . . , n que verifique
(3.8)

C1 (k)y1 (k + n 1) + C2 (k)y2 (k + n 1) + + Cn (k)yn (k + n 1) = 0

y finalmente
yH (k + n) = C1 (k)y1 (k + n) + C2 (k)y2 (k + n) + + Cn (k)yn (k + n)
si elegimos Ci (k), i = 1, 2, 3, . . . , n tal que
(3.9)

C1 (k)y1 (k + n) + C2 (k)y2 (k + n) + + Cn (k)yn (k + n) = b(k).

Asi tenemos un sistema de n ecuaciones (3.7),(3.8),. . . ,(3.9) linealmente independientes el

cual tiene solucion por el Teorema 3.3.3. Calculando los valores de Ci (k) podemos obtener
los Ci (k) por medio del operador suma indefinida 1 .
65

Ejemplo 3.5.1. Hallemos la soluci

on de la ecuaci
on
yk+2 5yk+1 + 6yk = k + 2.
Su polgono caracterstico
r2 5r + 6 = (r 2)(r 3),
y su respectivas soluciones son r1 = 2 y r2 = 3, luego la soluci
on homogenea asociada es
yk = c1 2k + c2 3k .
Tomando lo anterior tenemos que
yk = c1 (k)2k + c2 (k)3k ,
lo cual implica que
(3.10)

c1 (k)2k+1 + c2 (k)3k+1 = 0,

(3.11)

c1 (k)2k+2 + c2 (k)3k+2 = k + 2

Despejando Ci (k) en las ecuaciones (3.10) y (3.11) e igualando obtenemos

k + 2 3k+2
3k+1
c
(k)
=
k+2 c2 (k),
2
2k+1
2k+2
2
despejando se tiene
(3.12)

c2 (k) =

k+2
3k+1

c1 (k) =

k+2
.
2k+1

y
(3.13)

Aplicando la suma indefinida y el Teorema sumacion por partes en (3.12) y 3.13) se obtiene
c1 (k) =

4k + 7
k+2
+ D1 y c2 (k) = 3 k+1 + D2 .
k+2
2
2 3

Entonces la soluci
on general es:
yk =

k+2 k
4k + 7
k + 2 4k + 7
2 + 3 k+1 3k + 2k D1 + 3k D2 =
+
+ 2k D1 + 3k D2 ,
k+2
2
2 3
4
27

luego
yk = 2k D1 + 3k D2 +

k + 2 4k + 7
+
.
4
27

Ahora hallemos en el siguiente ejemplo una formula general para resolver ecuaciones homogeneas lineales de grado 2 utilizando el metodo de variacion de parametros.

66

Ejemplo 3.5.2. Sea la ecuaci

on
(3.14)

yk+2 + a1 (k)yk+1 + a2 (k)yk = bk

a2 (k) 6= 0,

supongamos que la soluci

on general es de la forma
yk = ck y1 (k) + dk y2 (k).
donde y1 (k) e y2 (k) son soluciones linealmente independiente.
Utilizando el metodo de variaci
on parametros resulta lo siguiente
ck y1 (k) + dk y2 (k) = 0,
ck y1 (k) + dk y2 (k) = bk .
El determinante de este sistema que lo denotaremos detk+1 (y1 , y2 ) es diferentes de cero,
as que
bk y2 (k + 2)
bk y1 (k + 1)
, dk =
.
detk+1 (y1 , y2 )
detk+1 (y1 , y2 )
Aplicando el operador de suma indefinida tenemos que
ck =

k1
X
bn y1 (n + 1)
,
ck = c0
det
(y
,
y
)
n+1
1
2
n=0

k1
X
bn y2 (n + 2)
dk = d0 +
,
det
(y
,
y
)
n+1
1
2
n=0

entonces la soluci
on general de la ecuaci
on (3.14) est
a dada por
k1
k1
X
X
bn y1 (n + 1)
bn y2 (n + 1)
yk = co y1 (k) + d0 y2 (k) y1 (k)
+ y2 (k)
.
detn+1 (y1 , y2 )
detn+1 (y1 , y2 )
n=0
n=0

3.5.2.

M
etodo de coeficiente indeterminados

Este metodo es analogo al que se desarrollo para ecuaciones de lineales de primer grado,
es esta seccion lo que se pretende es dar algunos ejemplos para ilustrarlo a ecuaciones de
grado superior.
Ejemplo 3.5.3. Sea la ecuaci
on
yk+2 + 5yk+1 + 6yk = k + 2,
la soluci
on homogenea asociada es
yh (k) = 2k C1 + 3k C2 .
Para hallar la soluci
on particular supongamos que este tiene siguiente forma
yp (k) = Ak + C,
entonces calculemos los valores A y C, reemplazando la soluci
on particular en la ecuaci
on
completa tenemos
A(k + 2) + C 5[A(k + 1) + C] + 6(Ak + C) = 2Ak 3A + 2C = k 2 + k,
67

lo cual implica que A =

Entonces

1
2

y C = 43 .
3
1
yp (k) = k + ,
2
4

luego la soluci
on general es
1
3
yk = 2k C1 + 3k C2 + k + .
2
4
Ejemplo 3.5.4. Sea la ecuaci
on yk+2 3yk+1 + yk = sen(k), la soluci
on de la ecuaci
on
homogenea asociada es:
yp (k) = C1 + C2 2k .
Supongamos que la soluci
on particular de la ecuaci
on completa tiene la siguiente forma
yp (k) = A sen(k) + B cos(k),
reemplazando en la ecuaci
on completa tenemos
7A sen(k) + 6B cos(k) = sen(K),
entonces las constantes A = 17 y B = 0. Lo caul implica que
yp (k) =

1
sen(k).
7

Luego la soluci
on general es:
yk = C1 + C2 2k +

1
sen(k).
7

Ejemplo 3.5.5. Sea la ecuaci

on yk+3 4yk+2 +3yk+1 = 2k . La soluci
on homogenea asociada
es:
yh (k) = 2k C1 + C2 .
Ahora supongamos que la soluci
on particular es de la forma
yp (k) = Ak2k .
Reemplazando en la ecuaci
on completa tenemos
A(k + 3)2k+3 4A(k + 2)2k+2 + 2A(k + 1)2k+1 = 2k ,
por lo tanto el valor de A =

1
.
14

Reemplazando en la soluci
on particular tenemos que
yp (k) =

1 k
k2 .
14

Luego la soluci
on general es
yk = 2k C1 + C2 +

68

1 k
k2 .
14

Bibliografa

[1] FERNANDEZ
PEREZ,
Carlos y VEGAS H. Francisco Jos
e. Ecuaciones Diferenciales y en Diferencia International Thomson Editores, 2003.

[2] GONZALEZ,
C. y BARRIOS, J. A. Analisis discreto en Economa Empresa International Thomson Editores, Espa
na, 2000.
[3] TAKAHASHI, Takehito. Ecuaciones en Diferencias con Aplicaciones Grupo Editorial Iberoamerica , Mexico, 1990.
[4] GOLDBERG, Samuel. Introduction to difference Equations Caurier Dover Publications, New York, 1986.

69