http://libreria-universitaria.blogspot.

com

7
Aplicaciones

7.1 APlICACION A LAS ECUACIONES DIFERENCIALES

Muchas leyes de la fsica, qunica, biologa y economa se expresan en trminos de

ecuaciones diferenciales, es decir, ecuaciones que comprenden funciones y sus derivadas.
El pfopsito de esta seccin es ilustrar una manera en la que es posible aplicar el lgebra
lineal para resolver ciertos sistemas de ecuaciones diferenciales . Esta sec.cin no es amplia,
pero debe servir para convencer al lector respecto a que el lgebra lineal tiene aplicaciones
slidas .
Una de las ecuaciones diferenciales ms sencillas es

y' = ay (7.1)

en donde y ={(x) es una funcin desconocida que se debe determinar , y' = dy fdx es su
derivada y a es una constante . Como la mayora de las ecuaciones diferenciales, (7.1) tie-
ne una infmidad de soluciones,. las cuales son funciones, en este caso, de l forma

(7.2)

en donde c es una constante arbitraria. Cada funcin de este tipo es una solucin de y' = ay
dado que

y' = cae ax = ay

Inversamente, toda solucin de y' = ay debe ser una funcin de la forma ceO X (ejercicio
7), de modo qe (7.2) describe todas las soluciones de y' =ay. A (7.2) se le da el nombre
de solucin general de y' = ay.
A veces el problema fsico que genera una ecuacin diferencial impone alguna con-
dicin agregada que permite aislar una solucin particular a partir de la solucin general.
Por ejemplo, si se requiere que la solucin de y' = ay satisfaga la condicin agregada

y(O) =3 (7.3)

329
 http://libreria-universitaria.blogspot.com

330 APLICACIONES

es decir , y = 3 cuando x = O, entonces, al sustituir estos valores en la solucin general

y = cell x , se obtiene un valor para c, a saber

3 = ceo = c
Por tanto,
y = 3e ax

es la nica solucin de y ' = ay que satisface la condicin agregada . Una condicin, como
la (7.3), que especifica el valor de la solucin como la (73), que especifica el valor de la
solucin en un punto se conoce como condicin inicial, y el problema de resolver una
ecuacin diferencial sujeta a una condicin inicial se denomina problema con valor inicial.
En esta seccin se considera la resolucin de sistemas de ecuaciones diferenciales
que tienen la forma

/1 = a ll Y1 + a 12Y2 + ... + a 1"y"

y~ = a 21 Y1 + a 22Y2 + . .. + a 2"y"
(7.4)

en donde Y1 = 11 (x), Y2 = 12 (x), ... 'Y n = In (x) son funciones que van a determinarse y
las a son constantes . En notacin matricial , (7 .4) se puede escribir

0, ms brevemente,
l <lj
Y2


y~
ra11
a 21
..
.
a"l a""
"j
~2" lY1
a
.~:2j
.Y n

Y' =AY

Ejemplo 1

a) Escnbase el sistema

/ 1 =, 3Y1
y~ = -2Yz
y~ = 5Y3

en notacion matricial .
b) Resulvase el sistema.
c) Hllese una solucin del sistema que satisfaga las condiciones iniciales y 1 (O) = 1,
Y2 (O) = 4 YY3 (O) -
2- .
 http://libreria-universitaria.blogspot.com

APLlCACION A LAS ECUACION ES DIFERENCIALES 331

Solucin a).

o bien ,
nr
. ~~ = ~
O
-2
O ~W:l (7. 5)

~ [~ ~]
O
Y' -2
Y
O

b) Debido a que cada ecuacin comprende slo una funcin desconocida, se puede resol-
ver cada una por separado. De (7 .2) se obtiene

Y = ce 3x
Y2 = c 2 e - 2x
5x
Y3 = c 3e

o , en notacin matricial,

c) A partir de las condiciones iniciales dadas , se obtiene

1 = y (O) = ceo = c
4 = yz (O) = c2e o = C2

- 2 = Y3 (O) = C3 eo = C3

de modo que la solucin que satisface las condiciones iniciales es

o , en notacin matricial,

lI
3X

y ] e zx]
y = Y2 = 4e -
[ V _ 2e 5 x
_ 3

El sistema de este ejemplo fu e fcil de resolver porque cada ecuacin comprenda

slo una funcin desconocida, y ste fue el caso porque la matriz de coeficientes (7.5 )
para el sistema fue diagonal . Pero
, cmo manejar un sistema
 http://libreria-universitaria.blogspot.com

332 APLICACIONES

Y' = AY

en el cual la matriz A no sea diagonal? La idea es sencilla : Intntese hacer una sustitucin
para Y que conduzca a un nuevo sistema con una matriz diagoal de los coeficientes;
resulvase este nuevo sistema ms sencillo y , a continuacin , sese esta solucin para
determinar la del sistema original.
La clase de sustitucin que se tiene en mente es

YI = PIIU 1 + P12 UZ + ... + Pl"U"

YZ = P21 U I + PZZU Z + .. . + P2"U"
(7.6)

o , en notacin matricial ,

fPll lU lj'
PIZ

e~ ~::
f Yl j P22 P lI z
PI,,]
2"

P,,2 ~"II ~I"

o, ms brevemente,

Y = PU

En esta sustitucin, las Pi son constantes que deben determinarse de tal modo que el
nuevo sistema qUe comprenda las funciones desconocidas u 1 ,u z , ... , Un tenga una matriz
diagonal de los coeficientes . Se deja como ejercicio para el estudiante derivar cada una de
las ecuaciones en (7.6) Y deduCir que

Y' = PU'

Si se hacen las sustituciones Y =PU y y' = PU' en el sistema original

Y' = AY

y.si se supone que Pes inversible , se obtiene

P U' = A(PU)

o bien,

o bien ,

U'= D U
 http://libreria-universitaria.blogspot.com

APLICACION A LAS ECUACIONES DIFERENCIALES 333

en donde D =p- 1 AP. Ahora queda claro cmo elegir P; si se desea que la nueva matriz de
los coeficientes D sea diagonal, se debe elegir P como una matriz que diagonalice a A.
Esto sugiere el procedimiento que sigue para resolver un sistema

Y' = AY

con una matriz A de los coeficientes que sea diagonalizable:

Ejemplo 2

aj Resulvase

/ 1 = YI + Y2
Y~ = 4YI - 2Y2

b j Hllese la solucin que satisfaga las condiciones iniciales y 1 (O) = 1, Y2 (O) =6.
Solucin aj . La matriz de los coeficientes para el sistema es

Como se analiz en la seccin 6.2, A ser diagonalizada por cualquier matriz P cuyas co-
lumnas sean eigenvectores linealmente independientes de A. Dado que

/, - 1
det(1.J - A) = . - 1 1 = A2 + 1, - 6 = (1, + 3)(/. - 2)
1
-4 I~ +2

los eigenvalores de A son A. = 2, A. = -3. Por definicin,

x=
[xlJ
x 2

es un eigenvector de A correspondiente a A. si y slo si x es una solucin no trivial de

(/J - A)x = 0 , es decir, de

l. -
[ - 4
1 l.-+12J[xJ [OJO1
X2 -
 http://libreria-universitaria.blogspot.com

334 APLICACIONES

Si A. = 2, este ~ist,ema se convierte en

[-41 - 41J [xlJ = [OJ Xz

Al resolverlo. da

Xl = t, .\z = t

de modo que

Por tanto ,

es una base para el eigenespacio correspondiente a A. = 2. De manera anloga, el lector

puede demostrar que

es una base para el eigenespacio correspondiente a A. = -3. Por tanto,

diagonaliza a A y

Por tanto , la sustitucin

y = PU y Y' = PU'

conduce al nuevo " sistema diagonal "

U'1 = 2u 1

u~ = - 3u z

Con base en (7 .2 ), la solucin de este sistema es

= e ('2.\ .1
[{ 1

11 2 =
(' 1 (' 2.\'

('2 e Jx
U = [ e 1e - 3.\
2
J
 http://libreria-universitaria.blogspot.com

APLICACIONES A LAS ECUACIONES DIFERENCIALES 335

de modo que la ecuacin Y = PU proporciona como la solucin para Y

o bien,

y, = c,e zx - ,tcze- 3x
(7,7)
)'z = c,e h + c ze- 3x

, b) Si se sustituyen las condiciones iniciales dadas en (7,7), se obtiene

e, - ,tc 2 =
c , + C2 = 6

Al resolver este sistema'se btiene

ez =4

de modo que, por (7.7), la solucin que satisface las condiciones iniciales es

)' , = 2e h _ e - 3x
J'2 = 2e 2x + 4e - 3x

En esta seccin se ha supuesto que la matiiz de los coeficienfes de y' = A Y es dia-

gonalizable , Si este no es el caso , es necesario aplicar otros mtodos para resolver el sis-
tema, Los mtodos de este tipo se estudian en textos ms avanzados,

EJERCICIOS 7.1

1. a) Resuelva el sistema

/1 = .\"1 + 4.\"2
.\"~ = 2.\"1 + 3\"2 ,
b) Halle la solucin que satisfaga las condiciones iniciales y , (O) = O, Y2 (O) = O,
2. a) Resuelva el sistema
.r'1 = .\"1 + 3\"2
.I'~ = 4\'1 + 5\'1
b) Halle la solucin que satisfaga las condiciones y,(O) = 2, y; (O) = 1
3. a) Resuelva el sistema

Y'2 = - 2\" 1 + .\'2 '

,.-:, = -2.\"1 + .\'3

 http://libreria-universitaria.blogspot.com

336 APLICACIONES

b ) Halle la solucin que satisfaga las condiciones iniciales y (O) = -l. Y2 (O) =
1. Y3 (0)=0.

4. Resuelva el sistema

y'l = 4YI + 2Y2 + 2y,

y = 2YI + 4Y2 + 2Y3
Y3 = 2YI + 2Y2 + 4Y3
S. Resuelva la ecuacin diferencial y" - y' - 6y = O. (Sugerencia.) Haga Y1 =
Y Y2 = y' Y. a continuacin. demuestre que

/\ = Y2
y = y" = / + 6y = y'\ + 6YI = 6y\ + Y2
6. Resuelva la ecuacin diferencial y" , - 6y" + 11Y' - 6y = O. (Sugerencia. Haga
YI =Y.Y2 = y'.Y3 = y " Y. a continuacin. demuestre que

y'l = Y2
Y'2 = Y3
Y3 = 6YI - llY2 + 6Y3 '
7. Pruebe que toda solucin de y' =ay tiene la forma Y ::; ceax . (Sugerencia . Sea y =
f(x) una solucin y demuestre que f(x)e -x es constante.)

8. Pruebe que si A es diagonalizable y

y = r~:j
)1"

satisface y' = A Y. entonces cada y ; es una combinacin lineal de e A1x e A2x .. ..

eAnx en donde Al A2 ..... An son eigenvalores de A.

7.2 APLlCACION A PROBLEMAS DE APROXiMACION; SERIES DE

FOURIER
En muchas aplicaciones se tiene inters en encontrar la mejor aproximacin posible sobre
un intervalo. para una funcin f. por medio de otra funcin que pertenece a alguna clase
especificada; por ejemplo:

a). Halle la mejor aproximacin posible para eX sobre [O. 1l. por medio- de un
polinomio de la forma ao + a 1 x + a2 X2 .
b ) Halle la mejor aproximacin posible para sen 1TX sobre [-l. 1], por medio de
una funcin de la forma ao + a \ eX + a2 e?- X + a3 e3x .
e) Halle la mejor aproximacin posible para Ix I
sobre [O . 21Tl. por medio de una
funcin de la forma ao + a \ sen x + a2 sen 2x + b 1 COS X + b2 cos 2x.