FACULTAD CS.
FISICAS Y MATEMATICAS
UNIVERSIDAD DE CHILE
MA26A Ecuaciones Diferenciales Ordinarias. Semestre 2007-2
Profesor: Axel Osses Auxiliares: Nicol
as Carre
no, Jorge Lemus.
Pauta Control 3
22 de Octubre de 2007
P1.- (i) Sean n N \ {0} y 6= 0. Use la transformada de Laplace para probar la igualdad:
tn cos t =
n n1
t
sen t .
Solucion: Se aplica transformada de Laplace a la expresion de la izquierda y se usa la propiedad de la
convolucion:
L[tn cos t](s) = L[tn ](s)L[cos t](s).
Estas u
ltimas dos transformadas son conocidas, con lo que se obtiene:
L[tn cos t](s) =
n!
s
n (n 1)!
=
sn+1 s2 + w2
sn s2 + 2
Reconociendo el producto de transformadas se obtiene:
L[tn cos t](s) =
n!
sn+1
s2
s
n
= L[tn1 sen t](s),
2
+w
de donde usando el Teorema de Lerch se obtiene:
tn cos t =
n n1
t
sen t .
(ii) Resuelva la siguiente ecuacion diferencial usando transformada de Laplace:
(
sen t 0 t <
00
0
ty 2y + ty = f (t), y(0) = 0, y
= 0, f (t) =
2
cos t t .
Lo primero es escribir la funcion f en una lnea, para lo que se usa la funcion escalon o Heaviside.
Recordando que L[Ha f (t a)] = eas F (s) y que cos(t ) = cos(t), sen(t ) = sen(t), se tiene:
f (t) = sen t + H (t)[ cos(t) + sen(t)] = sen t + H (t)[cos(t ) sen(t )],
aplicando transformada de Laplace:
L[f (t)](s) =
1
es s
es
+
.
s2 + 1 s2 + 1 s2 + 1
Ahora, llamando F (s) = L[y](s), al aplicar transformada de Laplace al lado izquierdo de la ecuacion se
obtiene:
d
d
[s2 F (s) sy(0) y 0 (0)] 2[sF (s) y(0)] [F (s)]
ds
ds
= 2sF (s) s2 F 0 (s) 2sF (s) F 0 (s) = (s2 + 1)F 0 (s) 4sF (s).
Juntando las 2 expresiones anteriores, se obtiene:
(s2 + 1)F 0 (s) + 4sF (s) =
F 0 (s) +
s2
es s
es
1
2
+ 2
.
+1 s +1 s +1
4s
1
es s
es
F (s) = 2
2
+ 2
.
2
2
+1
(s + 1)
(s + 1)
(s + 1)2
s2
Esto es una EDO de primer orden para F (s), que resolvemos con factor integrante. e
La ecuacion queda entonces:
[F (s)(s2 + 1)2 ]0 = 1 + (1 s)es ,
1
1 + s
e
,
F (s)(s2 + 1)2 = s ses +
2
F (s) =
(s2
4s
ds
s2 +1
= (s2 + 1)2 .
s
ses
es
+ (1 + ) 2 2
.
2
2
2
+ 1)
(s + 1)
(s + 1)2
Aplicando antitransformada y aplicando propiededes de traslacion y convolucion se obtiene:
y(t) = (sen() cos())(t)
1
1+
H (t)(sen() cos())(t ) +
H (t)(cos() cos())(t ).
P2.- Sea A(t) una matriz de n n funciones continuas para t R. Una matriz fundamental no canonica es una
matriz W (t) tal que:
(L)
W 0 = A(t)W, con W (0) invertible.
(i) Pruebe que W (t) existe y que si es la matriz fundamental canonica entonces
W (t) = (t)W (0).
Deduzca de ello que W (t) es invertible para todo t R.
Solucion: Ver que W (t) existe es una consecuencia del TEU (Teorema de Existencia y Unicidad). Para ello
denotemos por wi (t) a la i-esima columna de W (t), y consideremos el sistema lineal para cada i = 1 . . . n:
wi0 (t) =
wi (0) =
A(t)wi (t)
wi
donde {wi }ni=1 son l.i. Por el TEU cada sistema anterior tiene solucion, y es u
nica dada la condicion inicial.
As, es facil ver que W (t) satisface (SL). Sea ahora (t) = W (t) (t)W (0). Derivando:
0 (t) = W 0 (t) 0 (t)W (0)
0 (t) = A(t)W (t) A(t)(t)W (0)
0 (t) = A(t)(t),
o sea, satisface un sistema lineal, que llamammos (SL0 ). Ademas, (0) = W (0) (0)W (0) = 0. Por
el TEU, (SL0 ) tiene solucion u
nica, pero la funcion identicamente cero es solucion de (SL0 ), por lo que
podemos concluir que (t) = 0 t, es decir W (t) = (t) t. Por u
ltimo, sabemos que una matriz
fundamental es invertible si y solo si es invertible en un punto. Como W (t) es una matriz fundamental y
es invertible en 0, concluimos que W (t) es invertible t.
t
e
et
(ii) Sabiendo que W (t) =
es la solucion de (L), encuentre (t) y A(t).
et et
Compruebe que 0 (0) = A(0) para verificar sus calculos.
1 1
Solucion: Calculemos (t). De la parte anterior se tiene que (t) = W (t)W 1 (0). Como W (0) =
1 1
1/2 1/2
W (0)1 =
, se tiene:
1/2 1/2
(t)
et
=
et
t
et
et
e +et
2
et et
2
1/2 1/2
1/2 1/2
!
t
et e
2
et +et
2
De la misma forma calculemos A(t) = W 0 (t)W (t)1 :
A(t) =
et
et
0
1
2
et
et
1
0
et
2
et
2
e2
et
2
P3.- Una red genica esta formada de tres genes a, b y c. La actividad del gen a aumenta proporcionalmente a la
actividad del gen b y se indica por una flecha que va de b a a con una constante de proporcionalidad positiva
(b activa a a). Por otro lado, la actividad del gen c disminuye proporcionalmente a la actividad del gen a y se
dibuja un segmento terminado en una barra desde a hacia c con una constante de proporcionalidad negativa
(a inhibe a c). Y as sucesivamente. Puede tambien haber auto-inhibicion como ocurre con el gen c.
a
1
+1
+1
c
3
3
(i) Escriba el sistema de ecuaciones lineales en las variables xa , xb y xc representando las actividades de cada
gen. Encuentre la matriz A del sistema y su polinomio caracterstico.
Solucion: El sistema de ecuaciones vienen dado por:
x0a
x0b
x0c
= xb
= xc
= xa 3xb 3xc
o matricialmente
0
xa
xb
xc
0
= 0
1
1
0
3
0
xa
1 xb
3
xc
Encontremos el polinomio caracterstico:
1 3
0
1
3
(( + 3) + 3) 1
=
=
3 32 3 1
( + 1)3
(ii) Encuentre los valores y vectores propios (o vectores propios generalizados) de A. A partir de esto encuentre
la descomposicion diagonal (o en forma de Jordan) de A.
nico valor propio (con multiplicidad 3). Busquemos
Solucion: Se puede ver de la parte anterior que 1 es el u
los vectores propios:
1
0
1
1
1
3
0
v1
1 v2
2
v3
0
= 0
0
De donde
v1
v3
=
=
v2
v2
1
Fijando v2 = 1, obtenemos el primer vector propio: ~v1 = 1.
1
Encontremos ahora un vector propio generalizado:
1
0
1
0
v1
1 v2 =
2
v3
1
1
3
v1
v3
=
=
1
1
1
1 v2
1 v2
1
Fijando v2 = 0, ~v2 = 0 .
1
Como todava no completamos una base del espacio, buscamos un nuevo vector propio generalizado:
1
0
1
1
1
3
v1
0
1 v2 =
2
v3
v1
v3
1
Fijando v3 = 0, ~v3 = 0.
0
1
1 1
1
Si definimos P = 1 0 0 y J = 0
1 1 0
0
entonces:
1
0
1
= 1 v2
= v2
1
1
0
0
1 , la descomposicion en forma de Jordan es
1
A = P JP 1
(iii) Establezca la forma general que tiene la actividad de la red si no hay influencias externas a ella y muestre
que las actividades tienden a cero.
Solucion: Si no hay influencias externas, se debe resolver el sistema homogeneo para encontrar la forma
general de la red. Sabemos que la solucion viene dada por:
x(t) = eAt x0
xa
con x(t) = xb , y x0 una condicion inicial. Explcitamente, la solucion es:
xc
xa
xb =
xc
1
1
1
1 1
et
0 0 0
1 0
0
tet
et
0
t2 t
2e
t ~
c,
te
et
~c = P 1 x0
Cuando t , todos los terminos de la matriz anterior tienden a cero, con lo que se concluye que todas
las actividades son nulas en el lmite.
(iv) Pruebe que la red anterior es equivalente en la base propia a la red siguiente:
+1
+1
donde debe determinar que combinaciones de las actividades originales representan las nuevas actividades
de d, e y f . Finalmente, a partir de esta nueva red, diga cual de estas tres actividades converge en principio
mas lentamente a cero. (Nota: le puede servir definir las variables xd , xe y xf ).
Solucion: El sistema para esta nueva red es:
x0f
x0e
=
=
xf + xe
xe + xd
x0d
xd
o en forma matricial
0
xf
xe
xd
1
= 0
0
1
1
0
xf
0
1 xe
xd
1
xf
xa
Llamando y = xe y x(t) = xb .
xd
xc
El sistema es entonces:
y0
Jy
Por otro lado, el sistema original:
x0
x0
P 1 x0
=
=
=
Ax
P JP 1 x
JP 1 x
es decir, el nuevo sistema es equivalente a expresar las variables originales en otra base, mediante el cambio
y = P 1 x.
De la relacion x = P y, se puede deducir que
xa
xb
xc
= xf + xe + xd
= xf
= xf xe .
Ademas, para las nuevas variables y se tiene:
et
y= 0
0
tet
et
0
t2 t
2e
tet ~c1 ,
t
~c1 = y0 .
Dependiendo de las condiciones iniciales, se tendra que existira alguna componente yi , i = 1, 2, 3, tal que
t > t0 , yi (t) yj (t), i 6= j. Esa componente sera la que converge mas lentamente a cero.
