Ecuaciones lineales de segundo orden · Tema 22

Raíces reales, repetidas y complejas

Tres tipos de raíces características producen tres formas de solución y tres comportamientos dinámicos reconocibles.

01 · Un cálculo, tres caminos

El discriminante decide la forma de la solución

Para la ecuación lineal homogénea

a y″+b y′+c y=0,    a≠0

la propuesta y=ert conduce a a r²+b r+c=0. El discriminante

Δ=b²−4ac

determina cuántas raíces reales existen y, por lo tanto, cómo se construyen los dos modos independientes que necesita la solución general.

Mapa de las soluciones según el discriminante.
CondiciónRaícesSolución general real
Δ>0r₁≠r₂ realesC₁er₁t+C₂er₂t
Δ=0r real doble(C₁+C₂t)ert
Δ<0α±βieαt(C₁cos βt+C₂sin βt)

02 · Primer régimen

Dos raíces reales distintas: dos exponenciales

Si Δ>0, la fórmula cuadrática produce dos números reales diferentes:

r₁=(−b+√Δ)/(2a),    r₂=(−b−√Δ)/(2a)

Cada raíz genera una solución, y como r₁≠r₂, las exponenciales son independientes:

y(t)=C₁er₁t+C₂er₂t

Por ejemplo, para y″+y′−6y=0,

r²+r−6=(r−2)(r+3)=0
y=C₁e2t+C₂e−3t

Un modo crece y el otro decae. Salvo que C₁=0, el modo e2t termina dominando cuando t→∞.

03 · Segundo régimen

Una raíz repetida necesita un segundo modo nuevo

Si Δ=0, la única raíz es r=−b/(2a). La función ert aporta un modo, pero escribirla dos veces no crea dos constantes efectivas:

C₁ert+C₂ert=(C₁+C₂)ert

El segundo modo independiente es t ert. La familia correcta es

y(t)=(C₁+C₂t)ert

Para y″−4y′+4y=0, la característica es (r−2)²=0, de modo que

y=(C₁+C₂t)e2t

El factor t no es decorativo: recupera el grado de libertad perdido por la coincidencia de las raíces.

04 · Por qué aparece t

La raíz repetida puede verse como un límite

Supongamos dos raíces cercanas, r y r+h. Una combinación normalizada de sus modos es

(e(r+h)t−ert)/h

Al hacer h→0, esta diferencia cociente se convierte en la derivada respecto del parámetro r:

∂ert/∂r=t ert

Así, cuando dos modos exponenciales se fusionan, el modo tert sobrevive como dirección independiente. El wronskiano lo confirma:

W(ert,tert)=e2rt≠0

Esta idea de multiplicar por una potencia de t reaparece en ecuaciones de orden superior cuando una raíz tiene multiplicidad mayor.

05 · Tercer régimen

Raíces complejas: crecimiento y oscilación

Si Δ<0, las raíces forman un par conjugado:

r=α±βi,    α=−b/(2a),   β=√(−Δ)/(2|a|)

Los modos complejos e(α±βi)t se transforman mediante la identidad de Euler, eiθ=cos θ+i sin θ. Una base real es

eαtcos(βt),    eαtsin(βt)

Por tanto, la solución general real queda

y=eαt(C₁cos(βt)+C₂sin(βt))

Para y″+2y′+5y=0, las raíces son −1±2i y

y=e−t(C₁cos(2t)+C₂sin(2t))

06 · Leer α y β

La parte real es la envolvente; la imaginaria, la frecuencia

Parte real α

Amplitud exponencial

Si α<0, la envolvente decae; si α=0, la amplitud se conserva; si α>0, crece.

Parte imaginaria β

Frecuencia angular

Determina el ritmo de oscilación. El período es T=2π/β.

La combinación trigonométrica también puede escribirse como una sola función con fase:

C₁cos(βt)+C₂sin(βt)=A cos(βt−φ)
A=√(C₁²+C₂²)

Esta forma separa la amplitud inicial A de la fase φ. La envolvente completa es ±Aeαt.

07 · Seleccionar una trayectoria

Cómo aplicar y(0)=y₀, y′(0)=v₀ en cada caso

Las constantes cambian de interpretación según el tipo de raíz. Para datos en t=0:

Constantes determinadas por el estado inicial.
CasoPrimera constanteSegunda constante
r₁≠r₂C₁=(v₀−r₂y₀)/(r₁−r₂)C₂=(r₁y₀−v₀)/(r₁−r₂)
r dobleC₁=y₀C₂=v₀−ry₀
α±βiC₁=y₀C₂=(v₀−αy₀)/β

Estas fórmulas son atajos derivados de evaluar y e y′ en cero. Para un instante t₀≠0, conviene centrar las expresiones en τ=t−t₀ o resolver el sistema algebraico correspondiente.

08 · Tres PVI breves

La misma información inicial, tres respuestas diferentes

Raíces distintas

Para y″−y′−2y=0, y(0)=3, y′(0)=0, las raíces son 2 y −1. Las fórmulas dan C₁=1 y C₂=2:

y=e2t+2e−t

Raíz repetida

Para y″+2y′+y=0 con los mismos datos, r=−1, C₁=3 y C₂=3:

y=(3+3t)e−t

Raíces complejas

Para y″+2y′+5y=0, también con y(0)=3, y′(0)=0, tenemos α=−1, β=2, C₁=3 y C₂=3/2:

y=e−t(3cos(2t)+1,5sin(2t))

09 · Implementación robusta

Cerca de Δ=0, la clasificación requiere tolerancia

En aritmética exacta, el signo de Δ decide el caso. En una computadora, coeficientes aproximados pueden producir un valor diminuto como −2×10−15 cuando matemáticamente se esperaba cero.

Conviene comparar el discriminante con una tolerancia relativa a la escala:

ε=ε₀·max(1,b²,|4ac|)

Entonces se clasifica como repetida si |Δ|≤ε. También hay que evitar fórmulas con r₁−r₂ en el denominador cuando las raíces están casi fusionadas: aunque sean correctas algebraicamente, pueden amplificar errores de redondeo.

La clasificación numérica es una decisión computacional; no cambia el concepto matemático, pero sí la estabilidad de la implementación.

10 · Exploración interactiva

Laboratorio: atravesar los tres regímenes

Modificá p y q en y″+py′+qy=0. El laboratorio clasifica las raíces, construye la solución que cumple y(0)=y₀, y′(0)=v₀ y grafica su evolución.

EXPERIMENTO 22

Tres formas de solución

y″ + 2,00y′ + 5,00y = 0 · Δ=−16,00

Los resultados numéricos aparecen debajo.
TipoComplejas conjugadas
Raíces−1,00 ± 2,00i
ConstantesC₁=2,00; C₂=1,00
Valor en t=2y(2) ≈ −0,28

Oscilación con envolvente decreciente e⁻ᵗ.

La curva dorada es la solución del PVI y las líneas discontinuas muestran su envolvente cuando las raíces son complejas.

Preguntas para explorar

  1. Usá p=−1, q=−2. ¿Qué modos reales aparecen y cuál domina?
  2. Elegí p=2, q=1. ¿Qué raíz se repite y qué factor adicional aparece?
  3. Con q=5, cambiá el signo de p. ¿La envolvente crece o decrece?
  4. Mantené la ecuación y modificá los datos iniciales. ¿Cambian las raíces?

11 · Representación en código

Un solucionador exacto por casos

Clasificar las raíces con una tolerancia escalada

function clasificarRaices(a, b, c, epsilonBase = 1e-12) {
  if (![a, b, c].every(Number.isFinite) || a === 0) {
    throw new TypeError("Coeficientes inválidos");
  }
  const delta = b ** 2 - 4 * a * c;
  const tolerancia = epsilonBase * Math.max(1, b ** 2, Math.abs(4 * a * c));

  if (delta > tolerancia) {
    const raiz = Math.sqrt(delta);
    return { tipo: "distintas", r1: (-b + raiz) / (2 * a), r2: (-b - raiz) / (2 * a) };
  }
  if (delta < -tolerancia) {
    return { tipo: "complejas", alfa: -b / (2 * a), beta: Math.sqrt(-delta) / (2 * Math.abs(a)) };
  }
  return { tipo: "repetida", r: -b / (2 * a) };
}

console.log(clasificarRaices(1, 2, 5));

Determinar constantes desde el estado inicial

function ajustarCondiciones(raices, y0, v0) {
  if (raices.tipo === "distintas") {
    const { r1, r2 } = raices;
    return {
      c1: (v0 - r2 * y0) / (r1 - r2),
      c2: (r1 * y0 - v0) / (r1 - r2)
    };
  }
  if (raices.tipo === "repetida") {
    return { c1: y0, c2: v0 - raices.r * y0 };
  }
  return {
    c1: y0,
    c2: (v0 - raices.alfa * y0) / raices.beta
  };
}

const raices = clasificarRaices(1, 2, 5);
console.log(ajustarCondiciones(raices, 2, 0));

Evaluar cualquiera de las tres soluciones

function evaluarSolucion(raices, constantes, t) {
  const { c1, c2 } = constantes;
  if (raices.tipo === "distintas") {
    return c1 * Math.exp(raices.r1 * t)
      + c2 * Math.exp(raices.r2 * t);
  }
  if (raices.tipo === "repetida") {
    return (c1 + c2 * t) * Math.exp(raices.r * t);
  }
  return Math.exp(raices.alfa * t) * (
    c1 * Math.cos(raices.beta * t)
    + c2 * Math.sin(raices.beta * t)
  );
}

const constantes = ajustarCondiciones(raices, 2, 0);
console.log(evaluarSolucion(raices, constantes, 0));
console.log(evaluarSolucion(raices, constantes, 2));

12 · Comprensión

Confusiones frecuentes

«Una raíz doble aporta dos exponenciales iguales»

Dos copias de la misma función no son independientes. El segundo modo es tert.

«Las raíces complejas producen necesariamente soluciones complejas»

Al combinar el par conjugado se obtiene la familia real con seno y coseno.

«β controla el crecimiento»

La parte real α controla la envolvente; β controla la frecuencia angular.

«Si Δ<0, no hay solución»

No hay raíces reales, pero sí raíces complejas y soluciones reales oscilatorias.

«C₁ y C₂ siempre tienen la misma interpretación»

Multiplican bases diferentes en cada caso y sus fórmulas dependen de esas bases.

«En código basta comprobar delta === 0»

Los errores de representación hacen necesaria una comparación con tolerancia.

13 · Práctica guiada

Ejercicios de comprensión

Ejercicio 1: raíces distintas

Resolvé y″−5y′+6y=0.

Ver solución

r²−5r+6=(r−2)(r−3). Las raíces son 2 y 3, por lo que y=C₁e2t+C₂e3t.

Ejercicio 2: raíz repetida

Resolvé y″+6y′+9y=0.

Ver solución

(r+3)²=0. La raíz doble es −3 y y=(C₁+C₂t)e−3t.

Ejercicio 3: raíces complejas

Resolvé y″−4y′+13y=0.

Ver solución

Las raíces son 2±3i. Entonces y=e2t(C₁cos(3t)+C₂sin(3t)).

Ejercicio 4: PVI con raíz doble

Resolvé y″−2y′+y=0, y(0)=2, y′(0)=5.

Ver solución razonada

La raíz doble es r=1. C₁=2 y C₂=v₀−rC₁=5−2=3. Por tanto, y=(2+3t)et.

Ejercicio 5: PVI oscilatorio

Resolvé y″+4y=0, y(0)=3, y′(0)=−2.

Ver solución

α=0, β=2. Así, C₁=3 y C₂=−2/2=−1: y=3cos(2t)−sin(2t).

Ejercicio 6: interpretar raíces

Una ecuación tiene raíces −0,5±4i. Describí su comportamiento y calculá el período.

Ver solución razonada

La solución oscila dentro de una envolvente proporcional a e−0,5t, por lo que decae. Su frecuencia angular es 4 y el período es T=2π/4=π/2.

14 · Síntesis

Ideas para recordar

  • El signo de Δ=b²−4ac determina el tipo de raíces características.
  • Dos raíces reales distintas generan dos exponenciales independientes.
  • Una raíz doble requiere los modos ert y tert.
  • Las raíces α±βi producen eαtcos(βt) y eαtsin(βt).
  • La parte real controla crecimiento o decaimiento; la imaginaria controla oscilación.
  • En todos los casos deben existir dos modos linealmente independientes.
  • Las condiciones iniciales determinan los coeficientes de esos modos.
  • Cerca de Δ=0, una implementación debe utilizar tolerancias.
  • Las raíces permiten anticipar el comportamiento antes de calcular las constantes.

En el siguiente tema agregaremos fuerzas externas y separaremos la respuesta en solución complementaria y solución particular.