Figure 1:
Campos direccionales y soluciones
Las soluciones de la ecuación diferencial de primer orden dy/dx = f (x, y) pueden ser representadas gráficamente.
Dado un punto P (x, y), la ecuación diferencial da la pendiente de la recta tangente a la solución en el punto P .
Si m es una tal pendiente entonces la ecuación diferencial dice que
dy
m= = f (P ) = f (x, y).
dx P
La idea de un campo de direcciones o de pendientes es similar a la de un campo vectorial, donde f(x, y), es
representado por un vector en el punto (x, y). En el caso de una ecuación diferencial, f(x, y) da la pendiente de
la tangente a la curva en el punto (x, y) y se representa dicha tangente mediante un pequeño segmento de recta
en dicho punto. Si se representan varios de estos segmentos, se obtiene una representación visual del conjunto
de las soluciones de la ecuación diferencial.
En el siguiente gráfico se indica el campo de direcciones de la ecuación diferencial x′ = x2 − t, así como
también algunas de las soluciones de la ecuación diferencial.
EJEMPLOS
1. Trazar el campo direccional asociado a la ecuación diferencial y ′ = 1.
Solución. En este caso y ′ = f(x, y) es equivalente a y ′ = 1; es decir, que f (x, y) = 1 cualquiera que sea el
punto (x, y) del plano. Trazaremos entonces pequeños segmentos de pendiente 1 en un buen número de
puntos del plano. Se tiene entonces el gráfico siguiente:
2. Trazar el campo direccional asociado a la ecuación diferencial y ′ = x.
Solución. En este caso y′ = f (x, y) es equivalente a y ′ = x; es decir, que f (x, y) = x. Es decir que la
pendiente tangente a una curva solución en un punto cualquiera (x, y) , es igual a la abscisa del punto.
1
Figure 2: Campo dureccional de y ′ = 1 y akgunas curvas solución.
Construyamos una tabla de valores y dibujemos el campo direccional:
y y′ = y
−3 −3
−2 −2
−1 −1
0 0
1 1
2 2
3 3
3. Trazar el campo direccional asociado a la ecuación diferencial y ′ = y.
Solución. En este caso y ′ = f (x, y) es equivalente a y′ = y; es decir, que f(x, y) = y. Es decir que la
pendiente tangente a una curva solución en un punto cualquiera (x, y) , es igual a la ordenada del punto.
Construyamos una tabla de valores y dibujemos el campo direccional:
Trazaremos entonces pequeños segmentos de pendiente 1 en un buen número de puntos del plano. Se
tiene entonces el gráfico siguiente:
4. Trazar el campo direccional asociado a la ecuación diferencial y ′ = x + y.
Solución. En este caso y ′ = f(x, y) es equivalente a y′ = x + y; es decir, que f (x, y) = x + y. Es decir
que la pendiente tangente a una curva solución en un punto cualquiera (x, y) , es igual a la suma de las
coordenadas del punto.
Construyamos una tabla de valores y dibujemos el campo direccional:
x+y
5. Trazar el campo direccional asociado a la ecuación diferencial y ′ = .
x−y
x+y x+y
Solución. En este caso y ′ = f(x, y) es equivalente a y ′ = ; es decir, que f(x, y) = . Es decir
x−y x−y
que la pendiente tangente a una curva solución en un punto cualquiera (x, y) , es igual a la suma de las
coordenadas del punto dividida para la diferencia de las mismas.
2
Figure 3:
Construyamos una tabla de valores y dibujemos el campo direccional:
Observación. Trabajo tedioso, sinembargo existe otra herramienta que facilita este trabajo y son las
denominadas isoclinas.
6. En el siguiente gráfico se indica el campo de direcciones de la ecuación diferencial x′ = x2 − t, así como
también algunas de las soluciones de la ecuación diferencial.
Isoclinas
Definition 1 Curvas en las que la pendiente de la recta tangente a las soluciones es constante. En general,
para una EDO de orden 1, las isoclinas vienen dadas por F (x, y) = c.
Ejemplos
1. Consideremos nuevamente la ecuación diferencial y′ = x + y.
Solución. En este caso y ′ = f(x, y) es equivalente a y′ = x + y; es decir, que f (x, y) = x + y. Es decir
que la pendiente tangente a una curva solución en un punto cualquiera (x, y) , es igual a la suma de las
coordenadas del punto.
Construyamos una tabla de valores y dibujemos el campo direccional:
3
Figure 4:
Las isoclinas están dadas por x + y = k; o lo que es lo mismo, y = k − x, mismas que constituyen rectas
paralelas de pendiente −1.
k y =k−x
−3 y = −3 − x
−2 y = −2 − x
−1 y = −1 − x
0 y = −x
1 y =1−x
2 y =2−x
3 y =3−x
2. Consideremos la ecuación diferencial y′ = x2 + y 2 . Las isoclinas están dadas por x2 √+ y2 = C. En este
caso las isoclinas son círculos concéntricos de centro el origen de coordenadas y radio C, con C > 0.
x+y
3. Determinar las isoclinas de la ecuación diferencial y′ = .
x−y
4. y′ = y − x es una ecuación lineal. Las isoclinas son muy sencillas en este caso, y − x = m son rectas de
pendiente 1 y ordenada en el origen m, es decir la pendiente de las soluciones que pasan por ellas. El
dibujo de algunas isoclinas y soluciones es el que se muestra en la siguiente figura.
(x − 1) y
5. Dada la ecuación diferencial xy ′ + (1 − x) y = 0, las isoclinas son ahora = m, o si despejamos
x
y:
mx
y=
x−1
curvas que pasan por (0, 0) y tienen una asíntota vertical en el punto x = 1.
La pendiente en el origen no está definida. De la expresión obtenida en primer lugar, uno puede deducir
que la recta x = 1 es una isoclina de pendiente igual a cero. Asímismo, x = 0 es una isoclina de la ecuación
dx
( ) con pendiente igual a cero (o si se quiere de la ecuación inicial con pendiente ∞). El dibujo de las
dy
isoclinas y soluciones se muestra a continuación.
4
Figure 5:
Otras características de las soluciones pueden obtenerse de la propia ecuación, Por ejemplo, el conjunto de
puntos de inflexión está contenido en el conjunto de puntos que anulan a la segunda derivada y puede obtenerse
de la ecuación sin más que derivar:
∂f ∂f
y ′′ = + f (x, y) = 0.
∂x ∂y
EJEMPLOS
1. Dada la ecuación diferencial y ′ = y2 − x, las isoclinas de esta ecuación son las curvas y2 − x = m, es decir
paráboLas con eje horizontal. La derivada segunda de una solución es:
y ′′ = 2y′ y − 1 = 2y y 2 − x − 1
y por tanto, la ecuación de los posibles puntos de inflexión es la expresion anterior igualada a cero, es
decir;
1
x = y2 − .
2y
En la siguiente figura pueden verse las isoclinas, curva de puntos de inflexión y algunas soluciones.
2. Ejercicio. Representar las isoclinas de la ecuación diferencial ordinaria y ′ = 2x − y. ¿Qué tipo de curvas
son dichas isoclinas? Representar las isoclinas correspondientes a m = 0 y m = 2.¿Qué particularidad
tiene la correspondiente a m = 2?
3. Dada la ecuación diferencial y ′ = x2 − y, las curvas en las que la pendiente de la recta tangente a la
solución se mantiene constante son de la forma x2 − y = c, es decir, parábolas. En la Figura 2.3 puede
5
Figure 6: Campo direccional de y ′ = y y algunas curvas solución.
observarse el campo de pendientes de la ecuación diferencial.
Campo de pendientes y′ = x2 − y
Nos gustaría estudiar dónde están los puntos de inflexión de las soluciones. Sabemos que los puntos de
6
Figure 7: Campo direccional de y ′ = y así como algunas curvas solución.
inflexión aparecen cuando y ′′ = 0 Vemos entonces que
y′′ = 2x − y′ =⇒ y′′ = 2x − x2 + y
y′′ = 0 =⇒ 2x − x2 + y = 0 =⇒ y = x (x − 2) .
Los puntos de inflexión están en la parábola y = x (x − 2)). Veamos cómo hallar la solución de esta
ecuación diferencial.
y′ + y = x2 =⇒ ex y ′ + ex y = x2 ex =⇒ (ex y)′ = x2 ex
=⇒ ex y = x2 ex dx + C =⇒ y = e−x x2 ex dx + C .
4. Ejercicio. Construir el campo de direcciones y las curvas de nivel de la edo y ′ = sen(x) + y.
2x − y
5. Ejercicio resuelto. Dibujemos aproximadamente las soluciones de y ′ = .
x−y
2x − y 2−C
Solución. Trazamos algunas isoclinas = C, de donde y = x (rectas que pasan por el origen),
x−y 1−C
para diferentes valores de C y sobre cada una de ellas dibujamos algunos segmentos de pendiente C.
• C = 0 =⇒ y = 2x (segmentos horizontales: posibles máximos y mínimos de las soluciones)
3
• C = 1 =⇒ x = 0; C = −1 =⇒ y = x; ...
2
Una vez que sabemos que las isoclinas son rectas y = mx (es trivial ver que esto sucede en toda ecuación
homogenea) es más cómodo dibujar la recta de pendiente m que uno quiera y trazar sobre ella segmentos
2−m
de pendiente C = f (x, mx) = :
1−m
3
m = 0 =⇒ C = 2; m = 1 =⇒ K = +∞; m = −1 =⇒ C =
2
Las curvas tangentes a estos segmentos parecen ser cerradas (o tal vez espirales poco abiertas).
Podemos resolver la ecuación y comprobar (el ejemplo es poco práctico). Hay dos formas de hacerlo:
mirándola como ecuación homogénea o como exacta:
y
2− x
y′ = y o (2x − y) + (y − x) y′ = 0.
1− x
7
Figure 8: Campo direcional de y′ = x + y así como algunas curvas solución.
Por los dos caminos se llega a y2 −2xy +2x2 = C con lo que las soluciones son elipses. Con más propiedad,
√ √ √ √
cada una de ellas define de hecho dos soluciones en − C, C : y = x + C − x2 , y = x − C − x2
√ √ √
funciones definidas en − C, C no derivables en x = ± C.
1 1
Dada la ecuación diferencial y ′ = x−2y, las isoclinas son y = (x − C) (rectas de pendiente ). Dibujamos
2 2
1 1 1 1 3
las de C = −1, − 12 , 0, 1
2, 1 y 3
2 (que cortan respectivamente x = 0 en y = , , 0, − , − y −
2 4 4 2 4
).
1
Si C = , la recta y los segmentos trazados sobre ella tienen la misma pendiente y por tanto la isoclina
2
8
x+y
Figure 9: Campo direccional de y ′ = .
x−y
es solución de la ecuación (por ser tangente al campo de direcciones).
Podemos, también en este caso, resolver la ecuación (que es lineal) y comprobar. Bastará sumar la solución
x 1
general de la homogénea a la particular ya encontrada: y = − − + Ce−2x (a lo mismo llegaríamos con
2 4
la fórmula y = Ce−2x + e−2x xe2x dx).
EJEMPLO
3xy + 2y2
y′ = .
x2 + xy
Es una ecuación homogénea (tanto el numerador como el denominador son homogéneos de grado 2), con
3z + 2z 2 z (z + 2)
f (1, z) = =⇒ f (1, z) − z = .
1+z z+1
Las soluciones lineales son por tanto y = 0 y y = −2x. Las demás soluciones están dadas por
z+1 1 1 1 1
log |x| + C = dz = + dz? log |z (z + 2)| ,
z (z + 2) 2 z z+2 2
de donde
z (z + 2) = Cx2 .
y
con K = ±e2C constante no nula. Resolviendo esta ecuación cuadrática en z y sustituyendo z = obtenemos
x
finalmente
y = −x ± x 1 + Cx2 .
1
Esta solución está definida para todo x si C > 0, y para |x| < si C < 0. En este caso f es singular en
|C|
las rectas x = 0 y x + y = 0,´por lo que el abierto U estará contenido enuna de las cuatro regiones abiertas
determinadas por la intersección de estas rectas. En cada una de estas regiones el signo del radical en la fórmula
anterior está bien determinado: por ejemplo, si x < 0 y x + y > 0 deberá tomarse el signo −. Finalmente,
nótese que haciendo C = 0 en la fórmula obtenemos las dos soluciones lineales.
Sin embargo esta estrategia no siempre produce los frutos deseados. Sin ir más lejos, el problema
y 2 + (y ′ )2 = 0
y(0) = 0
no tiene solución y
y ′ = 3y2/3
y(0) = 0
tiene al menos dos soluciones dadas por y(x) = 0 y y(x) = x3 , Se verá más adelante bajo qué condiciones los
problemas de existencia y unicidad tienen asociados una única solución.
9
Figure 10:
Envolvente de una familia de curvas
Definition 2 (Envolvente). Una curva envolvente es aquella que “toca” a todas las curvas de una familia y
es tangente en los puntos de contacto.
Envolvente de una familia de curvas
Vamos a analizar como calcular la curva envolvente a una familia de curvas.
Método para hallar la envolvente: Como ya hemos visto, en general una familia de curvas viene dada por
una expresión de la forma F (x; y; c) = 0. Para hallar la envolvente a una familia de curvas basta con eliminar
c del sistema de ecuaciones
F (x; y; c) = 0
∂
F (x; y; c) = 0
∂c
donde F(x; y; c) = 0 define la familia de curvas de la cual queremos hallar la envolvente.
10
Figure 11:
Pongamos en práctica lo aprendido con un ejemplo:
EJEMPLO: Supongamos que tenemos un cañón antiaéreo en el origen de coordenadas que dispara un
proyectil con una velocidad inicial V . El ángulo θ en el que dispara el cañón es variable. Sabemos que la curva
que describe el proyectil es una parábola. El objetivo de este problema es hallar la zona en la que un avión
podría volar sin ser alcanzado por un proyectil. Es sencillo darse cuenta de que la zona de peligro viene descrita
por la que queda bajo la curva envolvente a la familia de parábolas que pueden describir los proyectiles lanzados
(Ver Figura 2.9).
Envolvente a una familia de parbolas
En primer lugar hallaremos la familia de curvas, tenemos en principio:
Movimiento horizontal: x(t) = V cos(θ) × t
g
Movimiento vertical: y(t) = V sen(θ) × t − t2 , donde t es el tiempo y g es la aceleración de la gravedad.
2
2V sen(θ)
y (tmax ) = 0 =⇒ tmax = .
g
Obtenemos así la ecuación de la posición:
g
σ (t) = (x (t) , y (t)) = V cos θ × t, t V sen θ − t , donde t ∈ [0; tmax ]
2
Despejando t e igualando términos se obtiene la ecuación de la trayectoria:
g
tan θ × x − x2 − y = 0
2V cos2 θ
2
1
Usando que 1 + tan2 θ = simplificamos la ecuación anterior:
cos2 θ
g 1 + tan2 θ 2
tan θ × x − x − y = 0.
2V 2
11
Figure 12: isoclinas de la edo y ′ = x2 + y 2 .
Hemos obtenido la familia de parábolas F (x; y; θ) = 0 donde
g tan2 θ + 1 2
F (x; y; θ) = [tan θ] x − x −y
2V 2
Para facilitar los cálculos llamamos c = tan θ y calculamos
∂ g
F (x; y; c) = x − 2 cx2
∂c V
∂ V2
Notamos que ∂c F (x; y; c) = 0 ⇐⇒ c =.
gx
Sustituyendo y simplificando llegamos a que la envolvente a nuestra familia de parábolas es
V2 g 2
y= − x .
2g 2V 2
12
Figure 13: Campo direccional, isoclinas y algunas curvas solución de y′ = x2 + y2 .
Figure 14:
Figure 15:
13
Figure 16:
Figure 17:
Figure 18:
14
Figure 19:
Figure 20:
15