Métodos de primer orden · Tema 16

Ecuaciones de Bernoulli

Cómo una potencia no lineal puede ocultar una ecuación lineal y qué precauciones exige volver desde la variable transformada a la solución original.

01 · Reconocer la familia

Una no linealidad con estructura aprovechable

Una ecuación de Bernoulli tiene la forma

y′+p(t)y=q(t)yn

donde n es constante. Si n≠0 y n≠1, la potencia hace que la ecuación no sea lineal en y, pero una sustitución específica la convierte en una ecuación lineal.

Reconocimiento de ecuaciones de Bernoulli.
Ecuación¿Bernoulli?Exponente
y′+ty=ety³Sín=3
y′+(2/t)y=t√ySí, donde la raíz esté definidan=1/2
y′+y=sin(t)Ya es linealn=0
y′+y²=tNo en esta formaLa potencia está en el lado lineal esperado.

02 · Los exponentes excepcionales

n=0 y n=1 ya son casos lineales

La sustitución de Bernoulli está diseñada para n≠0,1. Los exponentes excluidos no son problemas: producen ecuaciones lineales directamente.

n = 0

y′+py=q

Como y⁰=1, queda la ecuación lineal no homogénea estudiada en el tema anterior.

n = 1

y′+(p−q)y=0

El término derecho se mueve al izquierdo y queda una ecuación lineal homogénea.

También hay que considerar q(t)=0: la ecuación es homogénea lineal sin importar el valor escrito de n.

03 · Derivar la sustitución

Por qué v=y^(1−n) elimina la potencia

En una región donde las potencias involucradas están definidas y y≠0, dividimos la ecuación por yn:

y−ny′+p(t)y1−n=q(t)

Definimos

v=y1−n

Por la regla de la cadena:

v′=(1−n)y−ny′

Multiplicando la ecuación dividida por 1−n y sustituyendo:

v′+(1−n)p(t)v=(1−n)q(t)

Esta ecuación es lineal en v. Después de resolverla, debemos reconstruir y mediante

y=v1/(1−n)

04 · Método operativo

Resolver sin perder información

  1. Normalizar e identificar.

    Escribir y′+p(t)y=q(t)yⁿ y registrar n.

  2. Revisar soluciones especiales.

    Comprobar y=0 y cualquier restricción de las potencias antes de dividir.

  3. Transformar.

    Usar v=y1−n y v′=(1−n)y−ny′.

  4. Resolver la lineal.

    Aplicar factor integrante a v′+(1−n)pv=(1−n)q.

  5. Volver a y.

    Invertir la potencia, seleccionar ramas y aplicar la condición inicial.

  6. Declarar el intervalo.

    Excluir ceros de denominadores, radicandos inválidos y puntos donde la transformación falla.

  7. Comprobar en la ecuación original.

    No validar únicamente la ecuación transformada.

05 · Ejemplo completo

Una ecuación cuadrática que se vuelve lineal

Consideremos

y′+y=y²

Es Bernoulli con p=1, q=1 y n=2. Antes de dividir observamos que y=0 es solución. Para y≠0, usamos v=y−1=1/y. La ecuación transformada es:

v′−v=−1

Un factor integrante es e−t:

(e−tv)′=−e−t

Integrando:

e−tv=e−t+C  →  v=1+Cet

Volvemos a la variable original:

y(t)=1/[1+Cet]

Además de esta familia, debemos conservar y=0. El valor C=0 produce el otro equilibrio y=1.

06 · Condición inicial y dinámica

El dato inicial determina decaimiento o explosión

Para y(0)=y₀≠0,

C=(1−y₀)/y₀
Comportamiento de y′=y(y−1).
DatoConstante CComportamiento futuro
0<y₀<1PositivaLa solución decrece hacia cero.
y₀=1CeroPermanece en el equilibrio inestable y=1.
y₀>1NegativaEl denominador se anula y la solución explota en tiempo finito.

Para y₀>1, el tiempo de explosión se obtiene de 1+Cet=0:

t*=ln[y₀/(y₀−1)]

La ecuación transformada v=1/y atraviesa cero sin dificultad, pero ese cero corresponde a |y|→∞. Invertir la transformación revela el límite real de la solución original.

07 · Coeficientes variables

Bernoulli conduce al método lineal completo

Resolvamos, para t>0,

y′+(1/t)y=ty²

Con n=2 y v=1/y:

v′−(1/t)v=−t

El factor integrante lineal es μ=e∫−1/t dt=1/t. Entonces:

(v/t)′=−1

Integramos y reconstruimos:

v/t=−t+C  →  v=Ct−t²
y=1/(Ct−t²)

También y=0 satisface la ecuación original y no aparece al usar v=1/y. Cada solución no nula debe restringirse a un intervalo positivo que no contenga ceros del denominador.

08 · Ramas y dominios

Invertir una potencia puede introducir decisiones

La reconstrucción y=v1/(1−n) no siempre define una única función real para cualquier valor de v.

Exponente entero negativo
Puede introducir denominadores; los ceros de v se convierten en singularidades de y.
Raíz par
Requiere v≥0 y puede producir ramas y=±√v según la transformación.
Potencia fraccionaria original
La propia ecuación puede imponer y≥0 u otra restricción real.
Dato y₀=0
La división por yⁿ puede excluirlo; debe estudiarse directamente.

Ejemplo con n=1/2

Si v=√y, entonces reconstruir da y=v², pero la sustitución original exige v≥0. Una solución lineal para v que se vuelve negativa ya no representa esa rama de √y.

09 · Una conexión importante

La ecuación logística también es Bernoulli

El modelo logístico puede reorganizarse:

P′=rP−(r/K)P²
P′−rP=−(r/K)P²

Es una ecuación de Bernoulli con n=2, p=−r y q=−r/K. La sustitución v=1/P produce:

v′+rv=r/K

Esta ecuación lineal tiene equilibrio v=1/K. Al invertir, recuperamos la solución logística. La conexión explica por qué una ecuación poblacional no lineal puede resolverse con herramientas lineales.

10 · Exploración interactiva

Laboratorio: la transformación anticipa el destino

Explorá y′+y=y² con distintos valores positivos de y₀. El panel muestra simultáneamente y(t) y la variable lineal v(t)=1/y(t)=1+Cet. Cuando v llega a cero, y explota.

EXPERIMENTO 16

Variable original y transformada

y(t) = 1 / [1 + 1,00e^t]

Los resultados numéricos aparecen debajo.
Variable inicialv₀ = 2,00
Constante linealC = 1,00
Valor observadoy(1,00) = 0,27
Intervalo futuroSin explosión futura

Como 0<y₀<1, la solución decrece hacia cero.

La curva dorada representa y(t), la violeta representa v(t)=1/y(t), y una línea coral marca el equilibrio y=1.

Preguntas para explorar

  1. Elegí y₀=1. ¿Qué valores toman C, v e y?
  2. Probá y₀>1. ¿Cómo se relaciona el cruce de v=0 con la asíntota de y?
  3. Acercá y₀ a 1 desde arriba. ¿Qué ocurre con el tiempo de explosión?
  4. Para 0<y₀<1, ¿por qué v crece mientras y decrece?

11 · Representación en código

Resolver el PVI y declarar su intervalo

function resolverBernoulli(yInicial) {
  if (!(yInicial > 0)) throw new RangeError("Este ejemplo requiere y₀ > 0");
  const constante = (1 - yInicial) / yInicial;
  const explosion = yInicial > 1
    ? Math.log(yInicial / (yInicial - 1))
    : Infinity;

  return {
    constante,
    explosion,
    evaluar(t) {
      if (t >= explosion) throw new RangeError("t fuera del intervalo máximo");
      return 1 / (1 + constante * Math.exp(t));
    }
  };
}

const solucion = resolverBernoulli(0.5);
console.log(solucion.constante);
console.log(solucion.evaluar(1));

Comprobar la transformación lineal

const yInicial = 0.5;
const C = (1 - yInicial) / yInicial;
const v = t => 1 + C * Math.exp(t);
const derivadaV = t => C * Math.exp(t);
const puntos = [0, 0.5, 1, 2];

const cumpleLineal = puntos.every(t =>
  Math.abs(derivadaV(t) - v(t) + 1) < 1e-12
);

console.log(cumpleLineal);
console.log(v(0) === 1 / yInicial);

Transformar metadatos de una Bernoulli

function transformarBernoulli(p, q, n) {
  if (n === 0 || n === 1) {
    throw new RangeError("n=0 y n=1 ya son ecuaciones lineales");
  }
  return {
    coeficienteLineal(t) {
      return (1 - n) * p(t);
    },
    terminoIndependiente(t) {
      return (1 - n) * q(t);
    },
    reconstruir(v) {
      return v ** (1 / (1 - n));
    }
  };
}

const transformada = transformarBernoulli(t => 1, t => 1, 2);
console.log(transformada.coeficienteLineal(0));
console.log(transformada.terminoIndependiente(0));
console.log(transformada.reconstruir(2));

12 · Comprensión

Confusiones frecuentes

«Toda ecuación con yⁿ es de Bernoulli»

La ecuación completa debe organizarse como y′+p(t)y=q(t)yⁿ. La ubicación y los demás términos importan.

«La ecuación transformada usa el mismo p y q»

Los nuevos coeficientes son (1−n)p y (1−n)q.

«Puedo dividir por yⁿ sin revisar y=0»

La división puede eliminar una solución válida. Hay que comprobarla en la ecuación original.

«Resolver v resuelve automáticamente y»

Falta invertir la potencia, elegir ramas reales y excluir ceros o signos incompatibles.

«Si v es global, y también lo es»

Un cero de v=1/y corresponde a una explosión de y.

«n=0 o n=1 requieren la sustitución»

Esos valores ya producen ecuaciones lineales y se resuelven directamente.

13 · Práctica guiada

Ejercicios de transformación

Ejercicio 1: reconocer

Identificá p(t), q(t) y n en y′−2ty=t³y⁴.

Ver solución

p(t)=−2t, q(t)=t³ y n=4. La sustitución es v=y−3.

Ejercicio 2: obtener la lineal

Transformá y′+3y=6y².

Ver solución

Con n=2 y v=1/y, resulta v′−3v=−6.

Ejercicio 3: resolver

Resolvé y′+2y=2y² para soluciones no nulas.

Ver solución razonada

v=1/y produce v′−2v=−2. Entonces v=1+Ce2t y y=1/[1+Ce2t]. Además, y=0 es solución.

Ejercicio 4: condición inicial

Para y′+y=y², y(0)=1/3, encontrá la solución.

Ver solución

C=(1−1/3)/(1/3)=2, por lo que y=1/(1+2et).

Ejercicio 5: explosión

Con la misma ecuación y y(0)=2, calculá el extremo derecho del intervalo máximo.

Ver solución

t*=ln[2/(2−1)]=ln 2. La solución existe hacia adelante solo para t<ln 2.

Ejercicio 6: recuperar una solución

¿Qué solución se pierde al transformar la logística con v=1/P?

Ver solución

Se pierde P=0, porque 1/P no está definida allí. Debe agregarse después de comprobarla en la ecuación original.

14 · Síntesis

Ideas para recordar

  • La forma de Bernoulli es y′+p(t)y=q(t)yⁿ.
  • Para n=0 o n=1, la ecuación ya es lineal.
  • La sustitución es v=y1−n.
  • La ecuación transformada es v′+(1−n)pv=(1−n)q.
  • Dividir por yⁿ puede perder la solución y=0.
  • La reconstrucción puede exigir ramas, signos y restricciones.
  • Un cero de la variable transformada puede representar una explosión de la original.
  • La solución debe comprobarse en la ecuación original.
  • El modelo logístico es un caso importante de Bernoulli con n=2.

En el próximo tema estudiaremos ecuaciones exactas, donde la clave será reconocer una diferencial total y reconstruir una función potencial.