CURSO 2023-2024

Departamento de Fı́sica Interdisciplinar / Grado en Fı́sica - 61041088 Métodos Matemáticos I

Tema 6 Planos de fase de sistemas lineales bidimensionales

Problema 1
Consideremos el problema de valores iniciales para x1 (t) y x2 (t) dado por
dx1


 = −2x1 + x2
dt , x1 (0) = 1, x2 (0) = 0.
 dx2 = −5x1 + 2x2

dt
1. Obtener la solución del problema.

2. Hacer un esquema del diagrama de fases donde se muestre el punto crı́tico y la trayectoria que
corresponde a la solución particular obtenida antes. ¿De qué tipo es el punto crı́tico?
Solución
El sistema es lineal y de coeficientes constantes, por lo que el problema de valores iniciales propuesto
tiene solución única definida en R. Para resolverlo, es útil escribir el sistema en forma matricial:
   
d⃗x −2 1 x1
= A⃗x, donde A = es la matriz del sistema y ⃗x = .
dt −5 2 x2

1. La solución del problema vendrá dada por ⃗x(t) = c1 ⃗x1 (t)+c2 ⃗x2 (t), donde {⃗x1 , ⃗x2 } es un conjunto
fundamental de soluciones del sistema de ecuaciones diferenciales y c1 y c2 son constantes reales
que determinaremos a partir de las condiciones iniciales prescritas en t = 0.
Empecemos viendo si la matriz A es diagonalizable o no, para lo cual calculamos primero los
valores propios de A; si I denota la matriz identidad de orden dos,

det (A − λI) = 0 ⇒ λ2 + 1 = 0 ⇒ λ1 = i, λ2 = −i.

Puesto que los valores propios son complejos conjugados, la matriz es diagonalizable (en C); por
tanto, podemos encontrar dos vectores propios (complejos) linealmente independientes con los
que construir la solución general.
No obstante, en este caso de valores propios complejos, no es necesario buscar los dos vectores
propios independientes. Para obtener la solución general, basta considerar uno (cualesquiera) de
los valores propios, por ejemplo λ1 = i, y calcular un vector propio asociado a dicho autovalor,
⃗u1 . De esta forma, dos vectores solución linealmente independientes son:
   
⃗x1 (t) = Re ⃗u1 eλ1 t , ⃗x2 (t) = Im ⃗u1 eλ1 t .

Ası́, solo falta encontrar un autovector asociado al autovalor λ1 = i, resolviendo (A − λ1 I) ⃗u1 = ⃗0:
 
⃗ 1
(A − iI) ⃗u1 = 0 ⇒ ⃗u1 = (por ejemplo).
2+i

Entonces, sin más que operar, obtenemos:

   
1 cos(t) + i sen(t)
⃗u1 eλ1 t = eit = ,
2+i [2 cos(t) − sen(t)] + [2 sen(t) + cos(t)] i

1
 y, por lo tanto, los dos vectores solución son:
   
cos(t) sen(t)
⃗x1 (t) = , ⃗x2 (t) = .
2 cos(t) − sen(t) 2 sen(t) + cos(t)

En consecuencia, la solución general del sistema está dada por:

   
cos(t) sen(t)
⃗x(t) = c1 + c2 , c1 , c2 ∈ R,
2 cos(t) − sen(t) 2 sen(t) + cos(t)

o bien,

x1 (t) =c1 cos(t) + c2 sen(t),

x2 (t) =c1 [2 cos(t) − sen(t)] + c2 [2 sen(t) + cos(t)] , c1 , c2 ∈ R.

Solo nos falta determinar las constantes c1 y c2 , para lo cual usamos las condiciones iniciales:

x1 (0) = 1 ⇒ c1 = 1, x2 (0) = 0 ⇒ 2c1 + c2 = 0,

y, al resolver el sistema para c1 y c2 , se obtiene inmediatamente que c1 = 1 y c2 = −2.

En total, la solución del problema está dada por

x1 (t) = cos(t) − 2 sen(t), x2 (t) = −5 sen(t).

2. Los puntos crı́ticos se obtienen al resolver el sistema de ecuaciones


−2x1 + x2 = 0
, x1 = 0, x2 = 0,
−5x1 + 2x2 = 0

esto es, el sistema tiene un único punto crı́tico, C(0, 0), localizado en el origen de coordenadas
(esto ya era evidente, ya que la matriz A del sistema es regular, esto es, det(A) ̸= 0).
x2

x1

La naturaleza y estabilidad del punto crı́tico se deduce inmediatamente a partir de los valores
propios de la matriz A, los cuales se han obtenido en el apartado anterior: λ1 = i y λ2 = −i.
Puesto que los dos valores propios son imaginarios puros, se deduce inmediatamente que el punto
crı́tico C(0, 0) es un centro, y por tanto, es un punto crı́tico estable.

2
 La figura muestra el diagrama de fases del sistema, donde se observa el punto crı́tico del sistema,
C(0, 0), y la trayectoria que corresponde a la solución particular obtenida antes (obsérvese que
pasa por el punto (1, 0). Evidentemente, la trayectoria es una elipse. Y, además, se recorre en
sentido horario conforme aumenta el valor de t, lo cual se deduce fácilmente evaluando (en
las ecuaciones del sistema) el signo de dx1 /dt y dx2 /dt en cualquier punto (lo que indica el
crecimiento o decrecimiento de x1 y x2 ); por ejemplo, en el punto (1, 0), tenemos dx1 /dt = −2 < 0
(esto indica que x1 decrece en dicho punto) y dx2 /dt = −5 (indicando que x2 también decrece
en dicho punto), deduciendo entonces que la trayectoria se recorre en sentido horario.
NOTA: en el problema no se pedı́a el diagrama de fases exacto; en este problema lo importante
es mostrar que la trayectoria es una elipse centrada en el punto crı́tico (no es necesario hacer el
dibujo exacto) indicando el sentido de recorrido de la misma.

Problema 2
Considerar el sistema de ecuaciones diferenciales
 ′
x1 (t) = 2x1 (t) − x2 (t)
′ .
x2 (t) = 9x1 (t) − 4x2 (t)

1. Hallar la solución general, usando técnicas matriciales, del sistema.

2. Para cualquier solución particular, calcular su lı́mite cuando t tiende a +∞. ¿Es el punto (0, 0)
estable o inestable?

Solución
Escribamos el sistema en forma matricial:
   
d⃗x 2 −1 x1
= A⃗x, donde A = es la matriz del sistema y ⃗x = .
dt 9 −4 x2

1. La solución general del sistema está dada por ⃗x(t) = c1 ⃗x1 (t) + c2 ⃗x2 (t), donde {⃗x1 , ⃗x2 } es un
conjunto fundamental de soluciones del sistema y c1 y c2 son constantes reales.
Estudiamos si la matriz A es diagonalizable o no, para lo cual calcularemos sus valores propios
y vectores propios.

Valores propios: denotando por I la matriz identidad de orden dos,

det (A − λI) = 0 ⇒ (λ + 1)2 = 0 ⇒ λ = −1 (doble).

Vectores propios: al resolver (A − λI) ⃗u = ⃗0 para λ = −1 (con multiplicidad m = 2), se

obtiene que el subespacio vectorial propio asociado, V (λ = −1), tiene dimensión igual a
uno,
 
⃗ 1
(A + I) ⃗u = 0 ⇒ ⃗u = , o cualquier otro no nulo proporcional a él.
3

Como dim [V (λ = −1)] = 1 < m = 2, la matriz A no es diagonalizable, por lo que el vector propio
⃗u solo nos proporciona uno de los dos vectores solución linealmente independientes, digamos ⃗x1 :
 
1
λt
⃗x1 = ⃗u e = e−t ;
3

Ası́, el segundo vector solución ⃗x2 se debe obtener a partir de un vector propio generalizado.
En efecto, el defecto d del autovalor λ = −1 es d = m − dim [V (λ)] = 2 − 1 = 1, por lo que
podremos encontrar un vector propio generalizado, ⃗v1 ̸= ⃗0. Y más aún, la segunda solución ⃗x2
la buscaremos de la forma ⃗x2 = (t⃗v2 + ⃗v1 ) eλt , donde ⃗v2 ̸= ⃗0 es un vector propio ordinario que
se obtendrá una vez conocido ⃗v1 .

3
 El vector propio generalizado ⃗v1 se obtiene buscando un vector no nulo cualquiera que
pertenezca al núcleo de (A − λI)2 y que no dependa linealmente de ⃗u:
 
2 1
(A + I) ⃗v1 = ⃗0, con ⃗v1 ∈ / V (λ), ⃗v1 ̸= ⃗0 ⇒ ⃗v1 = (por ejemplo).
0

El vector propio ⃗v2 se obtendrá, entonces, de la siguiente forma:

 
3
⃗v2 = (A − λI) ⃗v1 = (A + I) ⃗v1 = ∈ V (λ).
9

Con esto, se construye la segunda solución linealmente independiente ⃗x2 :

 
1 + 3t
λt
⃗x2 = (t⃗v2 + ⃗v1 ) e = e−t .
9t

En consecuencia, la solución general del sistema está dada por:

   
1 −t 1 + 3t
⃗x(t) = c1 e + c2 e−t , c1 , c2 ∈ R,
3 9t

o bien,

x1 (t) =c1 e−t + c2 (1 + 3t) e−t ,

x2 (t) =3c1 e−t + 9c2 te−t , c1 , c2 ∈ R.

2. El lı́mite cuando t → +∞ se deduce inmediatamente a partir de la solución general obtenida

antes. En efecto, en dicho lı́mite, se tiene e−t → 0 y te−t → 0. En consecuencia, para cualquier
solución particular que se considere, tendremos lı́mt→+∞ x1 (t) = 0 y lı́mt→+∞ x2 (t) = 0.
Observación. La cuestión podrı́a haberse resuelto con facilidad de otra forma. En efecto, el
sistema tiene un único punto crı́tico, C(0, 0), (¿por qué?). Y ya que los dos valores propios
de la matriz A (en nuestro caso, iguales) son reales y negativos, el punto crı́tico C(0, 0) es
asintóticamente estable. Ası́, en el lı́mite t → +∞, cualquier trayectoria tenderá al punto (0, 0),
de donde se obtiene el resultado obtenido antes. Este comportamiento asintótico se observa en
la figura, donde se muestra el diagrama de fases del sistema (no se pedı́a en el problema)
x2

x1

4
Problema 3
Dado el sistema de ecuaciones diferenciales
 ′
x = −x/2 − y/2
.
y ′ = y/2

1. Determinar su solución general usando técnicas matriciales (autovalores y autovectores).

2. Representar el diagrama de fases del sistema, y clasificar la estabilidad de sus puntos crı́ticos.

Solución
Escribamos el sistema en forma matricial:
     
dx/dt x −1/2 −1/2
=A , donde A = es la matriz del sistema.
dy/dt y 0 1/2

1. La solución general del sistema está dada por

 
x(t)
= c1⃗s1 (t) + c2⃗s2 (t), c1 , c2 ∈ R,
y(t)

donde {⃗s1 , ⃗s2 } es un conjunto fundamental de soluciones del sistema.

Estudiamos si la matriz A es diagonalizable o no, para lo cual calculamos sus valores propios y
vectores propios.

Valores propios de A: obsérvese que la matriz A es triangular, por lo que sus valores propios
son justo los elementos de la diagonal principal, λ1 = −1/2 y λ2 = 1/2; y como λ1 ̸= λ2 ,
deducimos inmediatamente que la matriz A es diagonalizable.
Vectores propios de A: resolvemos (A − λI) ⃗u = ⃗0 para λ = λ1 y λ = λ2 ; y ya que A es
diagonalizable, cada uno de los subespacios propios, V (λ1 ) y V (λ2 ), tendrán dimensión uno.
Para λ = λ1 = −1/2, obtenemos:
 
⃗ 1
(A + (1/2)I) ⃗u1 = 0 ⇒ ⃗u1 = , o cualquier otro no nulo proporcional a él.
0

Para λ = λ2 = 1/2, obtenemos:

 
−1
(A − (1/2)I) ⃗u2 = ⃗0 ⇒ ⃗u2 = , u otro no nulo proporcional a él.
2

Puesto que la matriz A es diagonalizable, los vectores propios nos proporcionan ya dos vectores
solución linealmente independientes, ⃗s1 y ⃗s2 :
   
1 −1
⃗s1 = ⃗u1 eλ1 t = e−t/2 , ⃗s2 = ⃗u2 eλ2 t = et/2 .
0 2

Por tanto, la solución general del sistema es,

     
x(t) 1 −t/2 −1
= c1 e + c2 et/2 , c1 , c2 ∈ R,
y(t) 0 2

o bien,

x(t) =c1 e−t/2 − c2 et/2 ,

y(t) =2c2 et/2 , c1 , c2 ∈ R.

5
2. El sistema de ecuaciones diferenciales dado es un sistema autónomo lineal y homogéneo, y la
matriz del sistema, A, es regular (det(A) ̸= 0). De aquı́, se deduce inmediatamente que el sistema
tiene un único punto crı́tico, C(0, 0), localizado en el origen de coordenadas.
La naturaleza y estabilidad del punto crı́tico se deducen a partir de los valores propios de la
matriz A, obtenidos en el apartado anterior: λ1 = −1/2 y λ2 = 1/2. Puesto que los dos valores
propios son reales, distintos y de signo contrario, se deduce inmediatamente que el punto crı́tico
C(0, 0) es un punto de silla, y por tanto, es un punto crı́tico inestable.
En la figura se representa el diagrama de fases correspondiente a este sistema, donde se aprecia
el punto crı́tico del sistema, C(0, 0), y se muestran algunas trayectorias representativas (cada
una con un color distinto).
y

El diagrama de fases del sistema se dibuja de una forma muy sencilla teniendo en cuenta que
el sistema tiene un único punto crı́tico, C(0, 0), y el cual es un punto de silla. De esta forma,
tendremos cuatro trayectorias rectilı́neas tendiendo al punto crı́tico (dos de las cuales tienden al
punto crı́tico cuando t → +∞ y las otras dos cuando t → −∞), y entre las cuales habrá infinitas
trayectorias de tipo hiperbólico alejándose del punto crı́tico cuando |t| → ∞.
En primer lugar, dibujamos las cuatro trayectorias rectilı́neas (las cuales tienen la dirección de
los vectores propios):

Dos de ellas (la roja y la azul) tienen la dirección del vector propio asociado a λ1 = −1/2, es
decir, ⃗u1 = (1, 0)T (donde T indica traspuesta) y tienden al punto crı́tico cuando t → +∞;
en particular, la roja (resp. la azul) corresponde a considerar en la solución general c1 > 0
y c2 = 0 (resp. c1 < 0 y c2 = 0).
Las otras dos (la verde y la amarilla) tienen la dirección del vector propio asociado a
λ2 = 1/2, es decir, ⃗u2 = (−1, 2)T y tienden al punto crı́tico cuando t → −∞; en particular,
la verde (resp. la amarilla) corresponde a considerar en la solución c1 = 0 y c2 > 0 (resp.
c1 = 0 y c2 < 0).

Todas las demás trayectorias no son rectilı́neas, y se obtendrán al considerar en la solución

general valores no nulos de ambas constantes c1 y c2 (por ejemplo, la rosa se obtiene con c1 < 0
y c2 > 0, y similarmente todas las demás).
NOTA: En este problema no se pedı́a el diagrama de fases exacto, sólo se pedı́a un esquema
aproximado del mismo.

6
Problema 4
Dado el sistema de ecuaciones diferenciales
 ′ √
x = −3x/2
√ − 3y/2
.
y ′ = − 3x/2 − 5y/2

1. Determinar su solución general.

2. Representar el diagrama de fases del sistema, y clasificar la estabilidad de sus puntos crı́ticos.

Solución
Escribamos el sistema en forma matricial:
     √ 
dx/dt x −3/2 − 3/2
√
=A , donde A = es la matriz del sistema.
dy/dt y − 3/2 −5/2

1. La solución general del sistema está dada por

 
x(t)
= c1⃗s1 (t) + c2⃗s2 (t), c1 , c2 ∈ R,
y(t)

donde {⃗s1 , ⃗s2 } es un conjunto fundamental de soluciones del sistema.

Como la matriz A es una matriz real y simétrica, es diagonalizable; por tanto, existen dos vectores
propios de A linealmente independientes, con los cuales construiremos ⃗s1 y ⃗s2 .

Valores propios de A: si I es la matriz identidad de orden dos,

det (A − λI) = 0 ⇒ λ2 + 4λ + 3 = 0 ⇒ λ1 = −1, λ2 = −3.

Vectores propios de A: resolvemos (A − λI) ⃗u = ⃗0 para λ = λ1 y λ = λ2 , de lo cual

obtendremos dos vectores propios linealmente independientes (y, además, ortogonales).
Para λ = λ1 = −1, obtenemos:
 √ 
⃗ − 3
(A + I) ⃗u1 = 0 ⇒ ⃗u1 = , o cualquier otro no nulo proporcional a él.
1

Para λ = λ2 = −3, obtenemos:

 
(A + 3I) ⃗u2 = ⃗0 ⇒ ⃗u2 = √1 , u otro no nulo proporcional a él.
3

Los dos vectores propios calculados arriba nos proporcionan ya dos vectores solución linealmente
independientes, ⃗s1 y ⃗s2 :
 √   
− 3 −t 1
λ1 t
⃗s1 = ⃗u1 e = e , ⃗s2 = ⃗u2 e = √
λ2 t
e−3t .
1 3

Por tanto, la solución general del sistema es,

   √   
x(t) − 3 −t 1
= c1 e + c2 √ e−3t , c1 , c2 ∈ R,
y(t) 1 3

o bien,
√
x(t) = − 3c1 e−t + c2 e−3t ,
√
y(t) =c1 e−t + 3c2 e−3t , c1 , c2 ∈ R.

7
2. El sistema dado es un sistema autónomo lineal y homogéneo, y la matriz del sistema, A, es
regular (det(A) ̸= 0). De aquı́, se deduce que el sistema tiene un único punto crı́tico, C(0, 0),
localizado en el origen de coordenadas.
La naturaleza y estabilidad del punto crı́tico se deduce a partir de los valores propios de la matriz
A, los cuales se han obtenido en el apartado anterior: λ1 = −1 y λ2 = −3. En efecto, estamos
en el caso en el que λ2 < λ1 < 0; como los dos valores propios son reales, distintos y del mismo
signo, deducimos que el punto crı́tico C(0, 0) es un nodo; además, como tienen signo negativo,
se deduce que el nodo es asintóticamente estable.
En la figura se representa el diagrama de fases del sistema, donde se observa el punto crı́tico
C(0, 0) y se muestran algunas trayectorias representativas (cada una con un color distinto).
y

Para obtenerlo, tenemos en cuenta que el sistema tiene un único punto crı́tico, C(0, 0), y que
es un nodo; y ya que es asintóticamente estable, todas las trayectorias tienden al punto crı́tico
cuando t → +∞ (ver figura). De entre todas estas trayectorias, hay exactamente cuatro que son
rectilı́neas (las cuales tienen la dirección de los vectores propios); y entre ellas, habra una familia
de infinitas trayectorias de forma aproximadamente semiparabólica.
Primero, dibujamos las cuatro trayectorias rectilı́neas:

Dos de ellas (la roja y la azul) tienen la dirección del vector propio asociado a λ1 = −1, es
√
T
decir, ⃗u1 = − 3, 1 ; en particular, la roja (resp. la azul) corresponde a considerar en la
solución c1 > 0 y c2 = 0 (resp. c1 < 0 y c2 = 0).
las otras dos (la amarilla y la verde) tienen la dirección del vector propio asociado a λ2 = −3,
√
T
⃗u2 = 1, 3 ; en particular, la amarilla (resp. la verde) corresponde a considerar en la
solución general c1 = 0 y c2 > 0 (resp. c1 = 0 y c2 < 0).

Todas las demás trayectorias no son rectilı́neas, y se obtienen al considerar en la solución general
valores no nulos de las constantes c1 y c2 (por ejemplo, la rosa o la gris se obtienen con c1 < 0
y c2 > 0, y similarmente todas las demás). Puesto que λ2 < λ1 < 0, todas estas trayectorias no
rectilı́neas tienden al punto crı́tico (cuando t → +∞) según la dirección dada por el vector propio
asociado a λ1 = −1 (es decir, ⃗u1 ) y tienden al infinito (cuando t → −∞) según la dirección dada
por el vector propio asociado a λ2 = −3 (es decir, ⃗u2 ), como puede verse en la figura.
NOTA: En el problema sólo se pedı́a un esquema aproximado del diagrama de fases.

8
Problema 5
Dado el sistema de ecuaciones diferenciales
 ′
 x = 2x + y
y ′ = −x + 4y .
 ′
z = −x + y + 3z

1. Determinar la solución general del mismo usando técnicas matriciales

2. Discutir la establidad del sistema, exponiendo si es estable o inestable y qué significa que sea
estable o inestable según el caso.

Solución
Escribimos el sistema en forma matricial,
     
dx/dt x 2 1 0
 dy/dt  = A  y  , donde A =  −1 4 0  es la matriz del sistema.
dz/dt z −1 1 3

1. La solución general del sistema está dada por

 
x(t)
 y(t)  = c1⃗s1 (t) + c2⃗s2 (t) + c3⃗s3 (t),
z(t)

donde {⃗s1 , ⃗s2 , ⃗s3 } es un conjunto fundamental de soluciones del sistema de ecuaciones diferen-
ciales y c1 , c2 y c3 son constantes reales.
Como siempre, empezamos estudiando si la matriz A es diagonalizable o no, para lo cual calcu-
lamos sus valores propios y los correspondientes vectores propios.

Valores propios de A: se calculan en la forma usual, si ahora I es la matriz identidad de

orden tres,

det (A − λI) = 0 ⇒ (λ − 3)3 = 0 ⇒ λ=3 (triple).

Vectores propios de A: hay que resolver (A − λI) ⃗u = ⃗0 para λ = 3 (con ı́ndice de multipli-
cidad m = 3); al resolver, obtenemos que el subespacio vectorial propio asociado, V (λ = 3),
tiene dimensión dos, esto es, podemos encontrar dos vectores linealmente independientes:
   
0 1
⃗
(A − 3I) ⃗u = 0 ⇒ ⃗u1 =  0  , ⃗u2 =  1  (por ejemplo).
1 0

Como dimV (λ = 3) = 2 < m = 3, es claro que la matriz A no es diagonalizable, por lo

que los vectores propios ⃗u1 y ⃗u2 solo pueden proporcionarnos dos vectores solución linealmente
independientes:    
0 1
λt 3t λt
⃗s1 = ⃗u1 e =  0 e ,
 ⃗s2 = ⃗u2 e =  1  e3t .
1 0
Por tanto, el tercer vector solución ⃗s3 se debe obtener a partir de un vector propio generalizado.
En efecto, el defecto d del autovalor λ = 3 está dado por d = m − dim [V (λ)] = 3 − 2 = 1, por lo
que podremos encontrar un vector propio generalizado, ⃗v1 ̸= ⃗0. Además, la tercera solución ⃗s3
la buscaremos de la forma ⃗s3 = (t⃗v2 + ⃗v1 ) eλt , donde ⃗v2 ̸= ⃗0 es un vector propio ordinario que se
obtendrá una vez conocido ⃗v1 .

9
 El vector propio generalizado ⃗v1 se obtiene determinando un vector no nulo cualquiera que
pertenezca al núcleo de (A − λI)2 y que no dependa linealmente de ⃗u1 y de ⃗u2 :
 
1
(A − 3I)2 ⃗v1 = ⃗0, con ⃗v1 ∈/ V (λ = 3), ⃗v1 ̸= ⃗0 ⇒ ⃗v1 =  0  (por ejemplo).
0

Y el vector propio ⃗v2 se obtendrá ası́:

 
−1
⃗v2 = (A − λI) ⃗v1 = (A − 3I) ⃗v1 =  −1  ∈ V (λ = 3).
−1

Con esto, se construye la tercera solución linealmente independiente ⃗s3 :

 
1−t
⃗s3 = (t⃗v2 + ⃗v1 ) eλt =  −t  e3t .
−t

En total, la solución general del sistema es,

       
x(t) 0 1 1−t
 y(t)  = c1  0  e3t + c2  1  e3t + c3  −t  e3t , c1 , c2 , c3 ∈ R,
z(t) 1 0 −t

o bien,

x(t) =c2 e3t + c3 (1 − t) e3t ,

y(t) =c2 e3t − c3 te3t , c1 , c2 , c3 ∈ R.
3t 3t
z(t) =c1 e − c3 te ,

2. El sistema de ecuaciones diferenciales dado en el problema es un sistema autónomo lineal, ho-

mogéneo y tridimensional. Para poder responder a la cuestion, obtengamos los puntos crı́ticos y
estudiemos su estabilidad.
Los puntos crı́ticos se obtienen, por definición, resolviendo el sistema

 2x + y = 0
−x + 4y = 0 , x = 0, y = 0, z = 0,
−x + y + 3z = 0


es decir, el sistema tiene un único punto crı́tico, C(0, 0, 0), localizado en el origen de coordenadas
(esto ya era evidente, puesto que la matriz A del sistema dado es regular, esto es, verifica que
det(A) ̸= 0).
Para estudiar la estabilidad de este sistema, estudiamos el comportamiento asintótico de las
soluciones (que ya las conocemos del apartado anterior) cuando t → +∞. Es muy sencillo
comprobar que, en el lı́mite t → +∞, cualquier solución particular distinta de la trivial diverge
del punto crı́tico C(0, 0, 0) (téngase en cuenta que e3t y te3t tienden a infinito en el lı́mite
t → +∞). En este sentido, concluimos que el sistema de ecuaciones diferenciales dado en este
problema es inestable.

10
Problema 6
Dado el sistema de ecuaciones diferenciales
dx1


 = −x1 + 92 x2
dt .
 dx2 = −x1 + 2x2

dt
1. Determinar su solución general usando técnicas matriciales.
2. Representar el diagrama de fases del sistema, y clasificar la estabilidad de sus puntos crı́ticos.
Solución
Escribamos el sistema en forma matricial:
   
d⃗x −1 9/2 x1
= A⃗x, donde A = es la matriz del sistema y ⃗x = .
dt −1 2 x2
1. La solución general del sistema está dada por ⃗x(t) = c1 ⃗x1 (t) + c2 ⃗x2 (t), donde {⃗x1 , ⃗x2 } es un con-
junto fundamental de soluciones del sistema de ecuaciones diferenciales y c1 y c2 son constantes.
Veamos primero si la matriz A es diagonalizable o no. calculando los valores propios y vectores
propios de A. Los valores propios de A se obtienen de la forma usual:
5 1 3 1 3
det (A − λI) = 0 ⇒ λ2 − λ +
=0 ⇒ λ1 = + i, λ2 = − i,
2 2 2 2 2
donde I denota la matriz identidad de orden dos. Como los valores propios son complejos con-
jugados, la matriz es diagonalizable (en C); por tanto, es seguro que podremos encontrar dos
vectores propios (complejos) linealmente independientes con los que construir la solución general.
No obstante, por estar en el caso de valores propios complejos, no necesitamos obtener los dos
vectores propios independientes para construir la solución general. Es más sencillo considerar
uno cualquiera de los valores propios, por ejemplo λ1 = 1/2 + 3i/2, y calcular un vector propio
asociado a dicho autovalor, ⃗u1 . Entonces, dos vectores solución linealmente independientes ⃗x1 (t)
y ⃗x2 (t) estarán dados por
   
⃗x1 (t) = Re ⃗u1 eλ1 t , ⃗x2 (t) = Im ⃗u1 eλ1 t .

Ası́, solo falta encontrar un vector propio asociado al autovalor λ1 = 1/2 + 3i/2:
!
3 3
2 − 2 i
(A − λ1 I) ⃗u1 = ⃗0 ⇒ ⃗u1 = (por ejemplo).
1
Entonces, tras operar, obtenemos:
!
 !
3 3 3 3
+ 32 sen 3 3 3
− 32 cos 3




− i 1 3 cos 2t 2 t
+ 2 sen 2 t
2t i
⃗u1 eλ1 t = 2 2
e( 2 + 2 i)t = 2
3 3 et/2 ,
1 cos 2 t + i sen 2 t

y, por lo tanto, los dos vectores solución son:

3 3

3 3

! 3 3

!
− 32 cos 3


2 cos 2 t + 2
sen 2 t sen 2t 2t
⃗x1 (t) = et/2 , ⃗x2 (t) = 2
et/2 .
cos 32 t sen 23 t

En consecuencia, la solución general del sistema está dada por:

3 3

3 3

! 3 3

3 3

!
2 cos 2 t + 2
sen 2 t t/2 2 sen 2 t − 2
cos 2t
⃗x(t) = c1 e + c2 et/2 , c1 , c2 ∈ R,
cos 32 t sen 23 t
o bien,
         
3 t/2 3 3 3 t/2 3 3
x1 (t) = c1 e cos t + sen t + c2 e sen t − cos t ,
2 2 2 2 2 2
   
t/2 3 t/2 3
x2 (t) =c1 e cos t + c2 e sen t , c1 , c2 ∈ R.
2 2

11
 2. El sistema es autónomo lineal y homogéneo, y la matriz del sistema, A, es regular (det(A) ̸= 0).
De aquı́, deducimos que el sistema tiene un único punto crı́tico, C(0, 0), el cual está localizado
en el origen de coordenadas.
La naturaleza y estabilidad del punto crı́tico se deduce a partir de los valores propios de la matriz
A, ya obtenidos anteriormente: λ1 = 1/2 + 3i/2 y λ2 = 1/2 − 3i/2. Como los dos valores propios
son complejos conjugados con parte real positiva, se deduce inmediatamente que el punto crı́tico
C(0, 0) es un foco (inestable).
La figura muestra el diagrama de fases del sistema, donde se observa el punto crı́tico del sistema,
C(0, 0), y algunas trayectorias representativas. Como se esperaba, las trayectorias tienden (en
modo espiral) al punto crı́tico C(0, 0) cuando t → −∞; además, se recorren en sentido horario
conforme aumenta el valor de t, lo cual se puede deducir fácilmente observando el signo de dx1 /dt
y dx2 /dt en cualquier punto del plano.

x2

x1

NOTA: En el problema no se pedı́a el diagrama de fases exacto, sólo un esquema aproximado

del mismo.

Problema 7
Hacer una representación gráfica aproximada de la solución del problema de valores iniciales dado por
dx
= 1 − x3 , x(0) = 2.
dt
Solución
La ecuación diferencial dada es una ecuación de la forma dx/dt = f (x), donde la función f (x) está
dada por f (x) = 1−x3 = (1−x)(x2 +x+1) (ver Figura (A)). Puesto que f (x) es una función continua
y con derivada continua en R, el problema propuesto tiene solución única en algún entorno de t = 0.
Al tratarse de una ecuación autónoma, podemos hacer un esquema gráfico de la solución sin necesidad
de calcularla explı́citamente. Para ello, basta con obtener los puntos crı́ticos de la ecuación y el
diagrama de fases.
Por definición, los puntos crı́ticos de la ecuación salen al resolver f (x) = 0. Veamos:

(1 − x) x2 + x + 1 = 0


f (x) = 0 ⇒ ⇒ x = 1,

con lo que la ecuación dada tiene un único punto crı́tico, localizado en el punto x = 1.

12
Para obtener el diagrama de fases (notar que en realidad es un eje), basta estudiar el signo de f (x).
Para x < 1, tenemos f (x) > 0 (ver Figura (A)) y, por tanto, dx/dt > 0; esto implica que x(t) es
una función estrictamente creciente.
Para x > 1, tenemos f (x) < 0 (ver Figura (A)) y, entonces, dx/dt < 0; es decir, x(t) es una
función estrictamente decreciente.
De lo anterior, se obtiene inmediatamente el diagrama de fases de la ecuación (ver Figura (B)); y
deducimos también que el punto crı́tico x = 1 es un punto asintóticamente estable.

dx/dt

(A) f (x) (B) x

La gráfica de la solución, x(t), del problema de valores iniciales propuesto se muestra en la Figura (C),
y se dibuja ya de forma sencilla teniendo en cuenta lo siguiente.
Puesto que la condición inicial es x(0) = 2 > 1, del diagrama de fases deducimos que x(t) será
una función estrictamente decreciente; y, además, podemos garantizar que lı́mt→+∞ x(t) = 1.
Es más, puesto que d2 x/dt2 > 0 para todo t ∈ R (ver justo debajo), se deduce que la curva es
cóncava hacia arriba. En efecto, derivando en la propia ecuación diferencial obtenemos d2 x/dt2 :
d2 x
 
d dx d df dx
2
= = [f (x(t))] = = f ′ (x)f (x);
dt dt dt dt dx dt
y puesto que f (x) = 1 − x3 y f ′ (x) = −3x2 son funciones estrictamente negativas en todo el
intervalo x > 1, es claro que d2 x/dt2 > 0 en todo su dominio como habı́amos adelantado.

(C)

13
Observación. la Figura (C) ha sido obtenida resolviendo la ecuación de forma numérica. Nótese, que
la solución diverge para un valor finito de t (en particular, para t ≈ −0,1318). En el problema solo se
pedı́a un esquema aproximado de la gráfica.

14