Archivo de la etiqueta: discriminante

Ecuaciones Diferenciales I: Ecuaciones diferenciales homogéneas con coeficientes constantes

Por Omar González Franco

La esencia de las matemáticas no es hacer las cosas simples complicadas,
sino hacer las cosas complicadas simples.
– S. Gudder

Introducción

Continuando con nuestro desarrollo de métodos de resolución de ecuaciones diferenciales de orden superior, en particular de segundo orden, en esta entrada estudiaremos un método aplicado sólo a ecuaciones diferenciales lineales homogéneas de segundo orden con coeficientes constantes, es decir, de la forma

$$a \dfrac{d^{2}y}{dx^{2}} + b \dfrac{dy}{dx} + cy = 0$$

Con $a, b$ y $c$ contantes.

Antes de comenzar motivemos el análisis que desarrollaremos a lo largo de la entrada considerando primero las ecuaciones diferenciales de primer orden.

Ecuaciones diferenciales lineales homogéneas de primer orden con coeficientes constantes

En la primera unidad estudiamos las ecuaciones diferenciales lineales homogéneas de la forma

$$a_{1}(x) \dfrac{dy}{dx} + a_{0}(x) y = 0 \label{1} \tag{1}$$

Si $a_{1}(x) \neq 0$ podemos definir

$$P(x) = \dfrac{a_{0}(x)}{a_{1}(x)}$$

y reescribir la ecuación anterior en su forma canónica como

$$\dfrac{dy}{dx} + P(x) y = 0 \label{2} \tag{2}$$

Cuando estudiamos su método de resolución concluimos que la solución general de la ecuación (\ref{2}) es

$$y(x) = k e^{-\int{P(x) dx}} \label{3} \tag{3}$$

Consideremos ahora el caso en el que $a_{1}(x) = a \neq 0$ y $a_{0}(x) = b$ son números constantes, con esto la ecuación (\ref{1}) se puede escribir como

$$a \dfrac{dy}{dx} + b y = 0 \label{4} \tag{4}$$

Si definimos la constante $c = \dfrac{b}{a}$ podremos escribir la ecuación anterior en su forma canónica como

$$\dfrac{dy}{dx} + cy = 0 \label{5} \tag{5}$$

En la unidad anterior vimos que esta ecuación se puede resolver ya sea por variables separables o con ayuda de un factor integrante, sin embargo hay un método de resolución que sólo implica hacer un poco de álgebra.

Definiendo $k = -c$ la ecuación (\ref{5}) la podemos escribir de la siguiente manera.

$$\dfrac{dy}{dx} = ky \label{6} \tag{6}$$

La solución se puede intuir rápidamente, buscamos una función $y$, tal que su derivada sea igual a ella misma multiplicada por una constante, la función que satisface esto es

$$y(x) = e^{kx} \label{7} \tag{7}$$

Observamos que al derivarla recuperamos la ecuación (\ref{6}).

$$\dfrac{dy}{dx} = k e^{kx} = ky$$

Ahora que conocemos la solución, sustituyamos la función (\ref{7}) y su derivada en la ecuación (\ref{4}).

$$a k e^{kx} + b e^{kx} = 0$$

Factoricemos la función exponencial.

$$e^{kx} (ak + b) = 0$$

Como $e^{kx} \neq 0$, $\forall x \in \mathbb{R}$, entonces necesariamente

$$ak + b = 0 \label{8} \tag{8}$$

De donde

$$k = -\dfrac{b}{a}$$

Sustituyendo en la solución (\ref{7}), se tiene

$$y(x) = e^{-bx/a}$$

Por lo tanto, la solución general de la ecuación diferencial lineal homogénea de primer orden con coeficientes constantes (\ref{4}) es

$$y(x) = c_{1}e^{-bx/a} \label{9} \tag{9}$$

Ejemplo: Determinar la solución general de la ecuación diferencial

$$8 \dfrac{dy}{dx} + 16y = 0$$

Solución: Al ser una ecuación con coeficientes constantes, sabemos que la solución es de la forma (\ref{7}). Sustituimos $y$ y su derivada $\dfrac{dy}{dx}$ en la ecuación diferencial.

\begin{align*}
8ke^{kx} + 16e^{kx} &= 0 \\
e^{kx}(8k + 16) &= 0 \\
8k + 16 &= 0
\end{align*}

De la última relación despejamos a $k$, obteniendo

$$k = -\dfrac{16}{8} = -2$$

Una solución de la ecuación diferencial es

$$y(x) = e^{-2x}$$

Y por tanto su solución general en el intervalo $\delta = (-\infty, \infty)$ es

$$y(x) = c_{1} e^{-2x}$$

$\square$

Lo interesante es que esta idea de soluciones exponenciales se puede extender a ecuaciones diferenciales homogéneas de orden superior

$$a_{n} \dfrac{d^{n}y}{dx^{n}} + a_{n -1} \dfrac{d^{n -1}y}{dx^{n -1}} + \cdots + a_{1} \dfrac{dy}{dx} + a_{0} y = 0 \label{10} \tag{10}$$

donde los coeficientes $a_{i}$, $i = 0, 1, 2, …, n$, son constantes reales y $a_{n} \neq 0$.

Este análisis nos sirvió para motivar el siguiente método de resolución de ecuaciones diferenciales de segundo orden.

Ecuaciones diferenciales lineales homogéneas de segundo orden con coeficientes constantes

La ecuación diferencial que queremos resolver es de la forma

$$a \dfrac{d^{2}y}{dx^{2}} + b \dfrac{dy}{dx} + cy = 0 \label{11} \tag{11}$$

Con $a, b$ y $c$ constantes. Igual que antes, lo que intentamos es encontrar una solución de la forma

$$y(x) = e^{kx}$$

La primera y segunda derivada están dadas de la siguiente forma, respectivamente.

$$\dfrac{dy}{dx} = k e^{kx} \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}y}{dx^{2}} = k^{2} e^{kx}$$

Sustituimos en la ecuación (\ref{11}).

\begin{align*}
a(k^{2} e^{kx}) + b(ke^{kx}) + c(e^{kx}) &= 0 \\
e^{kx} (ak^{2} + bk + c) &= 0
\end{align*}

Como $e^{kx} \neq 0, \forall x \in \mathbb{R}$, entonces necesariamente

$$ak^{2} + bk + c = 0 \label{12} \tag{12}$$

El problema se ha reducido a obtener las raíces de esta ecuación, dicha ecuación tiene un nombre particular.

La dos raíces de la ecuación auxiliar son

$$k_{1} = \dfrac{-b + \sqrt{b^{2} -4ac}}{2a} \hspace{1cm} y \hspace{1cm} k_{2} = \dfrac{-b -\sqrt{b^{2} -4ac}}{2a} \label{13} \tag{13}$$

El discriminante puede ser positivo, cero o negativo, en cada caso ocurre lo siguiente.

  • Si $\Delta > 0$, entonces $k_{1}$ y $k_{2}$ son reales y distintos.
  • Si $\Delta = 0$, entonces $k_{1}$ y $k_{2}$ son reales e iguales, y
  • Si $\Delta < 0$, entonces $k_{1}$ y $k_{2}$ son números conjugados complejos.

Estudiemos cada caso y veamos el tipo de solución que se obtiene en cada uno.

Caso 1: Discriminante positivo

La solución de la ecuación (\ref{11}) es de la forma

$$y = e^{kx}$$

Si el discriminante es positivo, entonces al resolver la ecuación auxiliar obtendremos dos raíces reales y distintas $k_{1}$ y $k_{2}$, de manera que se tendrán dos soluciones de la ecuación diferencial (\ref{11}), dichas soluciones son

$$y_{1}(x) = e^{k_{1} x} \hspace{1cm} y \hspace{1cm} y_{2}(x) = e^{k_{2} x}$$

Notemos lo siguiente.

\begin{align*}
W(y_{1}, y_{2}) &= \begin{vmatrix}
e^{k_{1} x} & e^{k_{2} x} \\
k_{1} e^{k_{1} x} & k_{2} e^{k_{2} x} \end{vmatrix} \\
&= k_{2} e^{k_{2} x} e^{k_{1} x} -k_{1} e^{k_{1} x} e^{k_{2} x} \\
&= (k_{2} -k_{1}) e^{(k_{1} + k_{2}) x} \neq 0
\end{align*}

Como el Wronskiano es distinto de cero $\forall x \in \mathbb{R}$, entonces ambas funciones son linealmente independientes en $\mathbb{R}$, por tanto forman un conjunto fundamental de soluciones y la solución general de la ecuación diferencial (\ref{11}), en el caso en el que el discriminante es positivo, es

$$y(x) = c_{1} e^{k_{1} x} + c_{2} e^{k_{2} x} \label{15} \tag{15}$$

Con

$$k_{1} = \dfrac{-b + \sqrt{b^{2} -4ac}}{2a} \hspace{1cm} y \hspace{1cm} k_{2} = \dfrac{-b -\sqrt{b^{2} -4ac}}{2a}$$

Caso 2: Discriminante igual a cero

En este caso como

$$\Delta = b^{2} -4ac = 0$$

entonces,

$$k_{1} = k_{2} = -\dfrac{b}{2a}$$

De esta manera sólo obtendremos una solución exponencial.

$$y_{1}(x) = e^{k_{1} x} = e^{k_{2} x} = e^{-bx/2a} \label{16} \tag{16}$$

Para obtener la segunda solución vamos a aplicar el método de reducción de orden visto en la entrada anterior, en donde obtuvimos que una segunda solución linealmente independiente es de la forma

$$y_{2}(x) = y_{1}(x) \int{\dfrac{e^{-\int{P(x) dx}}}{y_{1}^{2}(x)} dx} \label{17} \tag{17}$$

En este caso de coeficientes constantes si la ecuación (\ref{11}) la dividimos por la constante $a \neq 0$ obtenemos

$$\dfrac{d^{2}y}{dx^{2}} + \dfrac{b}{a} \dfrac{dy}{dx} + \dfrac{c}{a}y = 0$$

Si definimos las constantes

$$P = \dfrac{b}{a} \hspace{1cm} y \hspace{1cm} Q = \dfrac{c}{a}$$

la ecuación anterior se puede reescribir como

$$\dfrac{d^{2}y}{dx^{2}} + P \dfrac{dy}{dx} + Qy = 0 \label{18} \tag{18}$$

Sustituyamos $P(x) = P$ y $y_{1}(x)$ en la segunda solución (\ref{17}).

\begin{align*}
y_{2}(x) &= e^{-bx/2a} \int{\dfrac{e^{-\int{\frac{b}{a} dx}}}{\left( e^{-bx/2a} \right)^{2}} dx} \\
&= e^{-bx/2a} \int{\dfrac{e^{-bx/a}}{e^{-bx/a}} dx} \\
&= e^{-bx/2a} \int{dx} \\
&= x e^{-bx/2a}
\end{align*}

Con esto, la segunda solución de la ecuación diferencial (\ref{11}), para este segundo caso es

$$y_{2}(x) = xe^{k_{1}x} = xe^{k_{2}x} = x e^{-bx/2a} \label{19} \tag{19}$$

Usando (\ref{16}) y (\ref{19}), notamos que

\begin{align*}
W(y_{1}, y_{2}) &= \begin{vmatrix}
e^{k_{1} x} & xe^{k_{1} x} \\
k_{1} e^{k_{1} x} & e^{k_{1} x} + x k_{1} e^{k_{1} x} \end{vmatrix} \\
&= e^{k_{1} x}(e^{k_{1} x} + x k_{1} e^{k_{1} x}) -x e^{k_{1} x} (k_{1} e^{k_{1} x}) \\
&= e^{2k_{1} x} + x k_{1} e^{2k_{1} x} -x k_{1} e^{2k_{1} x} \\
&= e^{2k_{1} x} \neq 0
\end{align*}

Como el Wronskiano es distinto de cero $\forall x \in \mathbb{R}$, entonces ambas funciones son linealmente independientes en $\mathbb{R}$, por tanto forman un conjunto fundamental de soluciones y la solución general de la ecuación diferencial (\ref{11}), en el caso en el que el discriminante es cero, es

$$y(x) = c_{1}e^{k_{1}x} + c_{2} xe^{k_{1}x} \label{20} \tag{20}$$

Con $k_{1} = -\dfrac{b}{2a}$.

Revisemos el último caso.

Caso 3: Discriminante negativo

Sabemos que

$$k_{1} = \dfrac{-b + \sqrt{b^{2} -4ac}}{2a} \hspace{1cm} y \hspace{1cm} k_{2} = \dfrac{-b -\sqrt{b^{2} -4ac}}{2a}$$

Si el discriminante es menor a cero definimos

$$-w = b^{2} -4ac \label{21} \tag{21}$$

Con $w \in \mathbb{R}^{+}$, tal que

$$k_{1} = \dfrac{-b + \sqrt{-w}}{2a} = -\dfrac{b}{2a} + i\dfrac{\sqrt{w}}{2a} \label{22} \tag{22}$$

y

$$k_{2} = \dfrac{-b -\sqrt{-w}}{2a} = -\dfrac{b}{2a} -i\dfrac{\sqrt{w}}{2a} \label{23} \tag{23}$$

Donde $i^{2} = -1$, definimos

$$\alpha = -\dfrac{b}{2a} \hspace{1cm} y \hspace{1cm} \beta = \dfrac{\sqrt{w}}{2a} \label{24} \tag{24}$$

De esta forma las raíces de la ecuación auxiliar son los números complejos conjugados

$$k_{1} = \alpha + i\beta \hspace{1cm} y \hspace{1cm} k_{2} = \alpha -i\beta \label{25} \tag{25}$$

Donde $\alpha$ y $\beta$ son números reales positivos. Dado que ambas raíces son distintas, similar al caso 1, tendremos dos soluciones dadas por

$$y_{1}(x) = e^{(\alpha + i\beta) x} \hspace{1cm} y \hspace{1cm} y_{2}(x) = e^{(\alpha -i\beta) x} \label{26} \tag{26}$$

De tarea moral demuestra que

$$W(y_{1}, y_{2}) = -2i \beta e^{2 \alpha x} \neq 0 \label{27} \tag{27}$$

y por tanto $\{ y_{1}, y_{2} \}$ forma un conjunto fundamental de soluciones, es así que la solución general está dada por la combinación lineal

$$y(x) = c_{1}e^{(a + i\beta) x} + c_{2}e^{(a -i\beta) x} \label{28} \tag{28}$$

A pesar de que la función anterior es la solución de la ecuación diferencial (\ref{11}), es común trabajar con una solución que contenga funciones reales en lugar de exponenciales complejas. Con el propósito de hallar una expresión real para la solución $y(x)$, vamos a considerar la formula de Euler

$$e^{i\theta} = \cos \theta + i\sin \theta \label{29} \tag{29}$$

con $\theta \in \mathbb{R}$. De esta formula y considerando las propiedades

$$\cos (-\beta x) = \cos (\beta x) \hspace{1cm} y \hspace{1cm} \sin (-\beta x) = -\sin (\beta x) \label{30} \tag{30}$$

es que podemos escribir las siguientes expresiones.

$$e^{i\beta x} = \cos (\beta x) + i \sin (\beta x) \hspace{1cm} y \hspace{1cm} e^{-i\beta x} = \cos (\beta x) -i \sin (\beta x) \label{31} \tag{31}$$

De las ecuaciones anteriores observamos que se cumplen las siguientes relaciones.

$$e^{i\beta x} + e^{-i\beta x} = 2 \cos (\beta x) \hspace{1cm} y \hspace{1cm} e^{i\beta x} -e^{-i\beta x} = 2i \sin (\beta x) \label{32} \tag{32}$$

Estas ecuaciones nos servirán mas adelante ya que, dada la solución (\ref{28}), si ocurre que $c_{1} = c_{2} = 1$, entonces se obtiene la solución

\begin{align*}
y_{1}(x) &= e^{(a + i\beta )x} + e^{(a -i\beta )x} \\
&= e^{ax} (e^{i\beta x} + e^{-i\beta x}) \\
&= 2 e^{ax} \cos (\beta x)
\end{align*}

Y si ocurre que $c_{1} = 1$ y $c_{2} = -1$, entonces se obtiene la solución

\begin{align*}
y_{2}(x) &= e^{(a + i\beta )x} -e^{(a -i\beta )x} \\
&= e^{ax}(e^{i\beta x} -e^{-i\beta x}) \\
&= 2ie^{ax}\sin (\beta x)
\end{align*}

Con estos resultados vemos que las funciones

$$g(x) = e^{\alpha x} \cos(\beta x) \hspace{1cm} y \hspace{1cm} h(x) = e^{\alpha x} \sin(\beta x) \label{33} \tag{33}$$

son ahora funciones reales y además de ello son soluciones de la ecuación diferencial (\ref{11}).

En la primer entrada de esta segunda unidad demostramos que un múltiplo constante

$$y(x) = c y_{1}(x)$$

de una solución $y_{1}$ de una ecuación diferencial lineal homogénea es también una solución, usando este resultado es que podemos asegurar que las funciones $C_{1}g(x)$ y $C_{2}h(x)$ son también solución. De tarea moral muestra que

$$W(C_{1} e^{\alpha x} \cos (\beta x), C_{2} e^{\alpha x} \sin (\beta x)) = C_{1}C_{2} \beta e^{2 \alpha x} \neq 0 \label{34} \tag{34}$$

Es decir, el Wronskiano de las soluciones es distinto de cero. De esta manera podemos afirmar que ambas funciones forman un conjunto fundamental de soluciones en cierto intervalo $\delta$ y, por lo tanto, podemos concluir que la solución general real de la ecuación diferencial (\ref{11}) en el caso en el que el discriminante es negativo es

$$y(x) = C_{1} e^{ax} \cos (\beta x) + C_{2} e^{ax} \sin (\beta x) = e^{ax} (C_{1} \cos (\beta x) + C_{2}\sin (\beta x)) \label{35} \tag{35}$$

Ejemplos

Realicemos una serie de ejemplos en los que tengamos que identificar a que caso pertenecen las ecuaciones diferenciales y así poder obtener su solución.

Ejemplo: Resolver el siguiente problema con valores iniciales.

$$\dfrac{d^{2}y}{dx^{2}} -4\dfrac{dy}{dx} -5y = 0; \hspace{1cm} y(1) = 0, \hspace{0.5cm} y^{\prime}(1) = 2$$

Solución: Consideremos la solución $y = e^{kx}$ y sus primeras dos derivadas

$$\dfrac{dy}{dx}= ke^{kx} \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}y}{dx^{2}}= k^{2}e^{kx}$$

Sustituimos en la ecuación diferencial.

$$k^{2} e^{kx} -4ke^{kx} -5e^{kx} = e^{kx}(k^{2} -4k -5) = 0$$

Como $e^{kx} \neq 0, \forall x \in \mathbb{R}$, entonces la ecuación auxiliar es

$$k^{2} -4k -5 = 0$$

Resolviendo para $k$, se tiene

$$k= \dfrac{4\pm \sqrt{16+20}}{2}= \dfrac{4\pm 6}{2}$$

De donde $k_{1} = 5$ y $k_{2} = -1$. Como ambas raíces son reales y distintas (y $\Delta = 36 > 0$), entonces la ecuación pertenece al caso 1 por lo que podemos afirmar que la solución es de la forma (\ref{15}).

$$y(x) = c_{1} e^{5x} + c_{2} e^{-x}$$

La derivada es

$$\dfrac{dy}{dx} = 5 c_{1} e^{5x} -c_{2}e^{-x}$$

Apliquemos las condiciones iniciales.

$$y(1) = c_{1}e^{5} + c_{2}e^{-1} = 0 \hspace{1cm} y \hspace{1cm} y^{\prime}(1) = 5 c_{1}e^{5} -c_{2}e^{-1} = 2$$

Resolviendo el sistema de ecuaciones obtendremos que

$$c_{1} = \dfrac{e^{-5}}{3} \hspace{1cm} y \hspace{1cm} c_{2} = -\dfrac{e^{1}}{3}$$

Por lo tanto, la solución particular del PVI es

$$y(x) = \dfrac{1}{3}e^{5(x -1)} -\dfrac{1}{3}e^{1 -x}$$

$\square$

Ejemplo: Resolver el siguiente problema con valores en la frontera.

$$\dfrac{d^{2}y}{dx^{2}} -10\dfrac{dy}{dx} + 25y = 0; \hspace{1cm} y(0) = 1, \hspace{0.5cm} y(1) = 0$$

Solución: Consideramos nuevamente la solución $y = e^{kx}$ y sus derivadas

$$\dfrac{dy}{dx}= ke^{kx} \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}y}{dx^{2}}= k^{2}e^{kx}$$

Sustituimos en la ecuación diferencial.

$$k^{2}e^{kx} -10ke^{kx} + 25e^{kx} = e^{kx}(k^{2} -10k + 25) = 0$$

Como $e^{kx} \neq 0, \forall x \in \mathbb{R}$, entonces la ecuación auxiliar es

$$k^{2} -10k +25 = 0$$

Resolviendo para $k$, se tiene

$$k = \dfrac{10 \pm \sqrt{100 -100}}{2} = \dfrac{10}{2} = 5$$

Notamos que $k_{1} = k_{2} = 5$, es decir, son raíces reales e iguales ($\Delta = 0$), por lo que estamos situados en el segundo caso y la solución está dada por la ecuación (\ref{20}).

$$y(x) = c_{1} e^{5x} + c_{2}x e^{5x}$$

Apliquemos las condiciones en la frontera.

$$y(0) = c_{1} = 1 \hspace{1cm} y \hspace{1cm} y(1) = c_{1}e^{5} + c_{2}e^{5} = 0$$

Obtenemos que $c_{1} = 1$ y $c_{2} = -1$. Por lo tanto, la solución particular del PVF es

$$y(x) = e^{5x} -x e^{5x}$$

Y además es única al no depender de ningún parámetro libre.

$\square$

Ejemplo: Obtener la solución general de la ecuación diferencial

$$\dfrac{d^{2}y}{dx^{2}} + 2 \dfrac{dy}{dx} + 3y = 0$$

Solución: Considerando la solución $y= e^{kx}$, calculando su primera y segunda derivada y sustituyendo en la ecuación diferencial obtendremos que la ecuación auxiliar es

$$k^{2} + 2k + 3 = 0$$

Resolviendo para $k$, tenemos

$$k = \dfrac{-2 \pm \sqrt{4 -12}}{2} = \dfrac{-2 \pm \sqrt{-8}}{2} = \dfrac{-2 \pm 2i\sqrt{2}}{2} = -1 \pm i \sqrt{2}$$

Las raíces son

$$k_{1} = -1 + i \sqrt{2} \hspace{1cm} y \hspace{1cm} k_{2} = -1 -i \sqrt{2}$$

Así mismo, identificamos que $\alpha = -1$ y $\beta = \sqrt{2}$. Al tratarse de raíces complejas notamos que corresponde al caso 3 y su solución compleja esta dada, de acuerdo a la ecuación (\ref{28}), como

$$y(x) = c_{1} e^{(-1 + i \sqrt{2})x} + c_{2} e^{(-1 -i \sqrt{2})x}$$

Sin embargo, una solución real es mucho más práctica de forma que si consideramos la ecuación (\ref{35}) concluimos que la solución general real de la ecuación diferencial es

$$y(x) = e^{-x}(C_{1} \cos (\sqrt{2} x) + C_{2} \sin (\sqrt{2} x))$$

$\square$

Ejemplo: Obtener la solución general de la ecuación diferencial

$$2\dfrac{d^{2}y}{dx^{2}} + 2\dfrac{dy}{dx} + y = 0$$

Solución: Considerando la solución $y= e^{kx}$, calculando su primera y segunda derivada y sustituyendo en la ecuación es como obtendremos que la ecuación auxiliar es

$$2k^{2} + 2k + 1 = 0$$

Resolvamos para $k$.

$$k = \dfrac{-2 \pm \sqrt{4 -4}}{4} = -\dfrac{2}{4}$$

En este caso

$$k_{1} = k_{2} = -\dfrac{1}{2}$$

es decir, las raíces son reales e iguales, de manera que estamos en el caso dos y por lo tanto la solución general es

$$y(x) = c_{1}e^{-x/2} + c_{2}x e^{-x/2}$$

$\square$

Ejemplo: Obtener la solución general de la ecuación diferencial

$$\dfrac{d^{2}y}{dx^{2}} -\dfrac{dy}{dx} -6y = 0$$

Solución: Repitiendo el mismo paso de siempre obtenemos que la ecuación auxiliar es

$$k^{2} -k -6 = 0$$

Resolviendo para $k$, tenemos

$$k = \dfrac{1 \pm \sqrt{1 + 24}}{2} = \dfrac{1 \pm 5}{2}$$

Las raíces son $k_{1} = 3$ y $k_{2} = -2$. Al ser reales y distintas notamos que estamos en circunstancias del caso uno y por tanto la solución general es

$$y(x) = c_{1}e^{3x} + c_{2}e^{-2x}$$

$\square$

Uno más

Ejemplo: Obtener la solución general de la ecuación diferencial

$$\dfrac{d^{2}y}{dx^{2}} + 8\dfrac{dy}{dx} + 16y = 0$$

Solución: La ecuación auxiliar en este caso es

$$k^{2} + 8k + 16= 0$$

Las raíces se obtienen de hacer

$$k = \dfrac{8 \pm \sqrt{64 -64}}{2} = \dfrac{8}{2} = 4$$

Las raíces son $k_{1} = k_{2} = 4$, al ser reales e iguales concluimos que la solución general es

$$y(x) = c_{1}e^{4x} + c_{2}x e^{4x}$$

$\square$

Revisemos dos ecuaciones con una forma particular.

Dos ecuaciones particulares

Las ecuaciones diferenciales

$$\dfrac{d^{2}y}{dx^{2}} + m^{2}y = 0 \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}y}{dx^{2}} -m^{2}y = 0$$

Con $m \in \mathbb{R}$, son importantes en matemáticas aplicadas. Vamos a obtener la forma de la solución de cada una de ellas.

Comencemos con la ecuación diferencial

$$\dfrac{d^{2}y}{dx^{2}} + m^{2}y = 0 \label{36} \tag{36}$$

Como es una ecuación con coeficientes constantes, entonces la solución es de la forma $y(x) = e^{kx}$, así la ecuación auxiliar en este caso es

$$k^{2} + m^{2} = 0$$

cuyas raíces son complejas y están dadas por

$$k_{1} = im \hspace{1cm} y \hspace{1cm} k_{2} = -im$$

De donde $\alpha = 0$ y $\beta = m$. Por lo tanto, de la ecuación (\ref{35}) concluimos que la solución de la ecuación (\ref{36}) es

$$y(x) = c_{1} \cos (mx) + c_{2}\sin (mx) \label{37} \tag{37}$$

Ejemplo: Obtener la solución general de la ecuación diferencial

$$\dfrac{d^{2}y}{dx^{2}} + 9y = 0$$

Solución: La ecuación a resolver es de la forma (\ref{36}) con $m = 3$. La ecuación auxiliar en este caso es

$$k^{2} + 9 = 0$$

De donde $k_{1} = i3$ y $k_{2} = -i3$, es decir, $\alpha = 0, \beta = m = 3$. De acuerdo a la ecuación (\ref{37}), la solución general de la ecuación diferencial es

$$y(x) = c_{1} \cos (3x) + c_{2} \sin (3x)$$

$\square$

Para el caso de la ecuación diferencial

$$\dfrac{d^{2}y}{dx^{2}} -m^{2}y = 0 \label{38} \tag{38}$$

La ecuación auxiliar es

$$k^{2} -m^{2} = 0$$

cuyas raíces son

$$k_{1} = m \hspace{1cm} y \hspace{1cm} k_{2} = -m$$

En este caso las raíces son reales y distintas, entonces la ecuación pertenece al caso 1 y por tanto la solución general está dada por la ecuación (\ref{15}), esto es

$$y(x) = c_{1} e^{mx} + c_{2} e^{-mx} \label{39} \tag{39}$$

Lo interesante de la ecuación (\ref{38}) es que si en su solución (\ref{39}) se elige $c_{1} = c_{2} = \dfrac{1}{2}$, se tiene

$$y_{1}(x) = \dfrac{1}{2} \left( e^{mx} + e^{-mx} \right) = \cosh (mx) \label{40} \tag{40}$$

Y si se elige $c_{1} = \dfrac{1}{2}$ y $c_{2} = -\dfrac{1}{2}$, entonces

$$y_{2}(x) = \dfrac{1}{2} \left( e^{mx} -e^{-mx} \right) = \sinh (mx) \label{41} \tag{41}$$

Se puede comprobar que $W(y_{1}, y_{2}) \neq 0$, lo que prueba que ambas soluciones son linealmente independientes en algún intervalo $\delta$. Este análisis nos permite establecer una forma alterna de la solución de la ecuación (\ref{38}), dada por

$$y(x) = c_{1} \cosh (mx) + c_{2} \sinh (mx) \label{42} \tag{42}$$

Ejemplo: Resolver el problema con valores iniciales

$$\dfrac{d^{2}y}{dx^{2}} -3y = 0; \hspace{1cm} y(0) = 1, \hspace{0.5cm} y^{\prime}(0) = 5$$

Usando primero la forma de la solución dada en (\ref{39}) y posteriormente resolverlo de nuevo usando la forma de la solución dada en (\ref{42}).

Solución: La ecuación a resolver es de la forma (\ref{38}) con $m = \sqrt{3}$. Considerando la solución $y = e^{kx}$, obtenemos que la ecuación auxiliar es

$$k^{2} -3 = 0$$

De donde $k_{1} = m = \sqrt{3}$ y $k_{2} = -m = -\sqrt{3}$. Usando la ecuación (\ref{39}) concluimos que la solución de la ecuación diferencial es

$$y(x) = c_{1} e^{\sqrt{3} x} + c_{2} e^{-\sqrt{3} x}$$

La derivada de la solución es

$$\dfrac{dy}{dx} = \sqrt{3} c_{1} e^{\sqrt{3} x} -\sqrt{3} c_{2} e^{-\sqrt{3} x}$$

Apliquemos las condiciones iniciales, por un lado

$$y(0) = c_{1} + c_{2} = 1 \hspace{1cm} \Rightarrow \hspace{1cm} c_{2} = 1 -c_{1}$$

por otro lado,

$$y^{\prime}(0) = \sqrt{3} c_{1} -\sqrt{3} c_{2} = 5$$

Si sustituimos el valor de $c_{2}$, tenemos

\begin{align*}
\sqrt{3} c_{1} -\sqrt{3} (1 -c_{1}) &= 5 \\
{\sqrt{3}}(c_{1} -1 +c_{1}) &= 5 \\
2 c_{1} -1 &= \dfrac{5}{\sqrt{3}} \\
2c_{1} &= \dfrac{5}{\sqrt{3}} +1 \\
c_{1} &= \dfrac{5}{2\sqrt{3}} + \dfrac{1}{2}
\end{align*}

Sustituyendo en $c_{2}$, se tiene

$$c_{2} = 1 -\left( \dfrac{5}{2\sqrt{3}} -\dfrac{1}{2} \right) = \dfrac{1}{2} -\dfrac{5}{2\sqrt{3}}$$

Por lo tanto, la solución particular del PVI es

$$y(x) = \dfrac{1}{2} \left( 1 + \dfrac{5}{\sqrt{3}} \right) e^{\sqrt{3} x} + \dfrac{1}{2} \left(1 -\dfrac{5}{\sqrt{3}} \right) e^{-\sqrt{3} x}$$

Si por otro lado consideramos la forma de la solución (\ref{42}), obtenemos que

$$y(x) = c_{1} \cosh(\sqrt{3} x) + c_{2} \sinh(\sqrt{3} x)$$

y su derivada

$$\dfrac{dy}{dx} = \sqrt{3} c_{1} \sinh(\sqrt{3} x) + \sqrt{3} c_{2} \cosh(\sqrt{3} x)$$

Aplicando las condiciones iniciales, se tiene

$$y(0) = c_{1} = 1 \hspace{1cm} y \hspace{1cm} y^{\prime}(0) = {\sqrt{3}c_{2}} = 5$$

De donde $c_{1} = 1$ y $c_{2} = \dfrac{5}{\sqrt{3}}$. Por lo tanto, la solución alterna del PVI es

$$y(x) = \cosh (\sqrt{3} x) + \dfrac{5}{\sqrt{3}} \sinh(\sqrt{3} x)$$

¿Qué relación hay entre las constantes $c_{1}$ y $c_{2}$ obtenidas en cada caso?.

$\square$

Para concluir la entrada revisemos brevemente estos mismos resultados para el caso de las ecuaciones de orden superior a dos.

Ecuaciones de orden superior

Es posible aplicar éste método de resolución para el caso de las ecuaciones diferenciales lineales homogéneas de orden superior a dos con coeficientes constantes, sin embargo desarrollar esta teoría desde cero puede ser muy complejo. Ahora que hemos sido muy detallados en el desarrollo de este método para el caso de las ecuaciones de segundo orden, vamos sólo a enunciar los posibles resultados para las ecuaciones de orden superior.

Es importante recordar que estamos estudiando las ecuaciones homogéneas y lineales con coeficientes constantes por lo que una ecuación de orden $n > 2$ con estas características tiene la siguiente forma.

$$a_{n} \dfrac{d^{n}y}{dx^{n}} + a_{n -1} \dfrac{d^{n -1}y}{dx^{n -1}} + \cdots + a_{2} \dfrac{d^{2}y}{dx^{2}} + a_{1} \dfrac{dy}{dx} + a_{0}y = 0 \label{43} \tag{43}$$

Donde $a_{i}, i = 0, 1, \cdots, n$ son constantes.

Igual que antes, se considera la solución $y(x) = e^{kx}$. Al sustituir esta función y las correspondientes derivadas en la ecuación (\ref{43}) obtendremos la siguiente ecuación auxiliar.

$$a_{n} k^{n} + a_{n -1} k^{n -1} + \cdots + a_{2} k^{2} + a_{1} k + a_{0} =0 \label{44} \tag{44}$$

Esta ecuación tendrá $n$ raíces, estas raíces, como sabemos, pueden ser reales o complejas, iguales o distintas. La solución general para caso está dada de la siguiente manera.

  • Si las raíces son reales y distintas, la solución estará dada por

$$y(x) = c_{1} e^{k_{1} x} + c_{2} e^{k_{2} x} + \cdots + c_{n} e^{k_{n}x} \label{45} \tag{45}$$

  • Si las raíces son reales e iguales, la solución estará dada por

$$y(x) = e^{kx}(c_{1} + c_{2}x + c_{3}x^{2} + \cdots + c_{n}x^{n -1}) \label{46} \tag{46}$$

En el caso de orden superior es posible tener raíces reales y de ellas que unas sean iguales y otras diferentes, en este caso se usan los dos puntos anteriores respectivamente.

Por ejemplo, supongamos que de una ecuación de sexto orden se obtienen seis raíces, tales que

$$k_{1} \neq k_{2} = k_{3} \neq k_{4} \hspace{1cm} y \hspace{1cm} k_{1} \neq k_{4} = k_{5} = k_{6}$$

Entonces la solución estaría dada por

$$y(x) = c_{1} e^{k_{1} x} + c_{2} e^{k_{2} x} + c_{3}x e^{k_{3} x} + c_{4} e^{k_{4} x} + c_{5}x e^{k_{5} x} + c_{6}x^{2} e^{k_{6} x}$$

Ya que, el par de raíces $k_{1} \neq k_{2}$ al ser diferentes genera la solución $c_{1} e^{k_{1} x} + c_{2} e^{k_{2} x}$, el par de raíces $k_{2} = k_{3}$ al ser iguales genera la solución $c_{2} e^{k_{2} x} + c_{3}x e^{k_{3} x}$, el par de raíces $k_{1} \neq k_{4}$ genera la solución $c_{1} e^{k_{1} x} + c_{4} e^{k_{4} x}$ y finalmente las raíces $k_{4} = k_{5} = k_{6}$ genera la solución $c_{4} e^{k_{4} x} + c_{5}x e^{k_{5} x} + c_{6}x^{2} e^{k_{6} x}$, esto de acuerdo a los dos puntos anteriores (\ref{45}) y (\ref{46}).

Finalmente,

  • Si las raíces son complejas, para cada par conjugado la solución es: $$y(x) = e^{\alpha x}(C_{1} \cos (\beta x) + C_{2} \sin (\beta x))$$ Si hay otro par igual, la función $$y(x) = e^{\alpha x}x (C_{1} \cos (\beta x) + C_{2} \sin (\beta x))$$ será solución y si hay otro par igual, la función $$y(x) = e^{\alpha x}x^{2} (C_{1} \cos (\beta x) + C_{2} \sin (\beta x))$$ será solución y así sucesivamente.

Ejemplo: Resolver el siguiente problema con valores iniciales.

$$\dfrac{d^{3}y}{dx^{3}} -7 \dfrac{d^{2}y}{dx^{2}} + 4\dfrac{dy}{dx} + 12y = 0$$

Con $y(0) = 1$, $y^{\prime}(0) = 0$ y $y^{\prime \prime}(0) = 36$.

Solución: Proponiendo la solución $y = e^{kx}$ la ecuación auxiliar que se obtiene es

$$k^{3} -7k^{2} +4k +12 = 0$$

Factorizando esta ecuación, se obtiene

$$(k + 1)(k -2)(k -6) = 0$$

De donde $k_{1} = -1$, $k_{2} = 2$ y $k_{3} = 6$. Como las raíces son reales y diferentes, entonces la solución de la ecuación diferencial es de la forma (\ref{45}).

$$y(x) = c_{1} e^{-x} + c_{2} e^{2x} + c_{3} e^{6x}$$

Para aplicar las condiciones iniciales calculemos la primera y segunda derivada de la solución.

$$\dfrac{dy}{dx} = -c_{1} e^{-x} + 2c_{2} e^{2x} + 6c_{3} e^{6x} \hspace{1cm} y \hspace{1cm} \dfrac{d^{2}y}{dx^{2}} = c_{1} e^{-x} + 4c_{2} e^{2x} + 36c_{3} e^{6x}$$

Apliquemos las condiciones iniciales.

$$y(0) = c_{1} e^{0} + c_{2} e^{0} + c_{3} e^{0} = c_{1} + c_{2} + c_{3} = 1$$

$$y^{\prime}(0) = -c_{1} e^{0} + 2c_{2} e^{0} + 6c_{3} e^{0} = -c_{1} + 2c_{2} + 6c_{3} = 0$$

$$y^{\prime \prime}(0) = c_{1} e^{0} + 4c_{2} e^{0} + 36c_{3} e^{0} = c_{1} + 4c_{2} + 36c_{3} = 36$$

El sistema que se obtiene es

\begin{align*}
c_{1} + c_{2} + c_{3} &= 1 \\
-c_{1} + 2c_{2} + 6c_{3} &= 0 \\
c_{1} + 4c_{2} + 36c_{3} &= 36 \\
\end{align*}

Resolviendo el sistema obtendremos que

$$c_{1} = \dfrac{16}{7}, \hspace{1cm} c_{2} = -\dfrac{5}{2} \hspace{1cm} y \hspace{1cm} c_{3} = \dfrac{17}{14}$$

Por lo tanto, la solución particular de la ecuación diferencial de orden 3 es

$$y(x) = \dfrac{16}{7} e^{-x} -\dfrac{5}{2} e^{2x} + \dfrac{17}{14} e^{6x}$$

$\square$

La dificultad de resolver ecuaciones de orden mayor a 2 realmente radica en que se vuelve más complicado encontrar las raíces de la ecuación auxiliar y resolver el sistema de ecuaciones que se genere en problemas con valores iniciales o con valores en la frontera. Un método para factorizar la ecuación auxiliar y obtener las raíces puede ser la división sintética, así como el método de Gauss – Jordan para obtener las soluciones de un sistema de ecuaciones lineales.

Tarea moral

Los siguientes ejercicios no forman parte de la evaluación del curso, pero servirán para entender mucho mejor los conceptos vistos en esta entrada, así como temas posteriores.

  1. Obtener la solución general de las siguientes ecuaciones diferenciales con coeficientes constantes.
  • $4 \dfrac{d^{2}y}{dx^{2}} + \dfrac{dy}{dx} = 0$
  • $3 \dfrac{d^{2}y}{dx^{2}} + 2 \dfrac{dy}{dx} + y = 0$
  • $\dfrac{d^{2}y}{dx^{2}} -4 \dfrac{dy}{dx} + 13y = 0$
  • $\dfrac{d^{2}y}{dx^{2}} -2 \sqrt{3} \dfrac{dy}{dx} + 3y = 0$
  • $\dfrac{d^{2}y}{dx^{2}} + 4 \dfrac{dy}{dx} -y = 0$
  • $\dfrac{d^{2}y}{dx^{2}} -4 \dfrac{dy}{dx} + 5y = 0$
  • $\dfrac{d^{2}y}{dx^{2}} -6 \dfrac{dy}{dx} + 13y = 0$
  • $\dfrac{d^{2}y}{dx^{2}} -\dfrac{4}{3} \dfrac{dy}{dx} + \dfrac{4}{9} y = 0$
  1. Resolver los siguientes problemas con valores iniciales.
  • $\dfrac{d^{2}y}{dx^{2}} + y = 0; \hspace{1cm} y \left( \dfrac{\pi}{3} \right) = 0, \hspace{0.5cm} y^{\prime} \left( \dfrac{\pi}{3} \right) = 2$
  • $\dfrac{d^{3}y}{dx^{3}} + 2 \dfrac{d^{2}y}{dx^{2}} -5 \dfrac{dy}{dx} -6y = 0; \hspace{1cm} y(0) = y^{\prime}(0) = 0, \hspace{0.5cm} y^{\prime \prime} (0) = 1$
  1. Resolver los siguientes problemas con valores en la frontera.
  • $\dfrac{d^{2}y}{dx^{2}} + 4y = 0; \hspace{1cm} y(0) = 0, \hspace{0.5cm} y(\pi) = 0$
  • $\dfrac{d^{2}y}{dx^{2}} -2 \dfrac{dy}{dx} + 2y = 0; \hspace{1cm} y(0) = 1, \hspace{0.5cm} y(\pi) = 1$
  1. Resolver el siguiente problema con valores iniciales usando primero la forma de la solución dada en (\ref{39}) y posteriormente resolverlo de nuevo usando la forma de la solución dada en (\ref{42}).
  • $\dfrac{d^{2}y}{dx^{2}} -y = 0; \hspace{1cm} y(0) = 1, \hspace{0.5cm} y^{\prime}(1) = 0$

Más adelante…

Ahora sabemos resolver ecuaciones diferenciales lineales y homogéneas de orden superior con coeficientes constantes, en la siguiente entrada estudiaremos este mismo tipo de ecuaciones, pero en el caso no homogéneo.

Para resolver ecuaciones no homogéneas existen dos métodos, el primero de ellos se denomina coeficientes indeterminados, este método suele no ser tan complicado, sin embargo esta limitado a ciertas funciones $g(x)$. Un segundo método se conoce como variación de parámetros siendo una extensión del método visto en la unidad anterior y que sirve para cualquier función $g(x)$, es decir, es un método general. En la siguiente entrada comenzaremos presentando el primer método.

Entradas relacionadas

Agradecimientos

Trabajo realizado con el apoyo del Programa UNAM-DGAPA-PAPIME PE104522 «Hacia una modalidad a distancia de la Licenciatura en Matemáticas de la FC-UNAM – Etapa 2»

Álgebra Superior II: Desigualdades de polinomios reales

Por Leonardo Ignacio Martínez Sandoval

Introducción

En la entrada anterior mostramos el teorema de factorización para polinomios con coeficientes reales. Lo que haremos ahora es ver que podemos aplicarlo en la resolución de desigualdades de polinomios en $\mathbb{R}[x]$. El objetivo es que, al final de la entrada, entendamos cómo se pueden resolver problemas como los siguientes:

Problema 1. Determina todos los números $x$ en $\mathbb{R}$ para los cuales $$x^6-12x^4-49x^2-30 > 3x^5-48x^3-51x+6.$$

Problema 2. Determina todos los números $x$ en $\mathbb{R}$ para los cuales $$\frac{1}{x}>x^3-x^2+1.$$

Antes de hablar de resolución de desigualdades de polinomios, veremos una forma alternativa de factorizar en $\mathbb{R}[x]$ usando potencias.

Teorema de factorización de polinomios reales con potencias

De acuerdo al teorema de factorización en $\mathbb{R}[x]$, un polinomio $p(x)$ se puede factorizar de manera única en factores lineales y factores cuadráticos con discriminante negativo. De ser necesario, podemos agrupar los factores lineales iguales y reordenarlos para llegar a una factorización de la forma $$a(x-r_1)^{\alpha_1}\cdots(x-r_m)^{\alpha_m}(x^2-b_1x+c_1)\cdots (x^2-b_{n}x+c_{n}),$$ en donde:

  • $a$ es un real distinto de cero,
  • $\alpha_1,\ldots,\alpha_m$ y $n$ son enteros positivos tales que $2n+\sum_{i=1}^m \alpha_i$ es igual al grado de $p(x)$,
  • para cada $i$ en $\{1,\ldots,m\}$ se tiene que $r_i$ es raíz real de $p(x)$ y $r_1<r_2<\ldots<r_m$
  • para cada $j$ en $ \{1,\ldots,n\}$ se tiene que $b_j,c_j$ son reales tales que $b_j^2-4c_j<0$.

Observa que los $r_i$ son ahora distintos y que están ordenados como $r_1<\ldots<r_m$. De aquí, obtenemos que $(x-r_i)^{\alpha_i}$ es la mayor potencia del factor lineal $x-r_i$ que divide a $p(x)$. Este número $\alpha_i$ se usa frecuentemente, y merece una definición por separado.

Definición. Sea $p(x)$ un polinomio en $\mathbb{R}[x]$ y $r$ una raíz de $p(x)$. La multiplicidad de $r$ como raíz de $p(x)$ es el mayor entero $\alpha$ tal que $$(x-r)^\alpha \mid p(x).$$ Decimos también que $r$ es una raíz de multiplicidad $\alpha$.

Ejemplo. El polinomio $k(x)=x^4-x^3-3x^2+5x-2$ se factoriza como $(x-1)^3(x+2)$. Así, la multiplicidad de $1$ como raíz de $k(x)$ es $3$. Además, $-2$ es una raíz de $k(x)$ de multiplicidad $1$.

$\triangle$

Después hablaremos de una forma práctica en la que podemos encontrar la multiplicidad de una raíz, cuando hablemos de continuidad de polinomios y sus derivadas.

Desigualdades de polinomios reales factorizados

Supongamos que tenemos un polinomio $p(x)$ no constante en $\mathbb{R}[x]$ para el cual conocemos su factorización en la forma $$a(x-r_1)^{\alpha_1}\cdots(x-r_m)^{\alpha_m}(x^2-b_1x+c_1)\cdots (x^2-b_{n}x+c_{n}),$$ y que queremos determinar para qué valores reales $r$ se cumple que $$p(r)>0.$$

Daremos por cierto el siguiente resultado, que demostraremos cuando hablemos de continuidad de polinomios.

Proposición. Las evaluaciones en reales de un polinomio cuadrático y mónico en $\mathbb{R}[x]$ de discriminante negativo, siempre son positivas.

Lo que nos dice este resultado es que, para fines de la desigualdad que queremos resolver, podemos ignorar los factores cuadráticos en la factorización de $p(x)$ pues

$$a(x-r_1)^{\alpha_1}\cdots(x-r_m)^{\alpha_m}(x^2-b_1x+c_1)\cdots (x^2-b_{n}x+c_{n})$$ y $$a(x-r_1)^{\alpha_1}\cdots(x-r_m)^{\alpha_m}$$ tienen el mismo signo.

Por la miasma razón, podemos ignorar aquellos factores lineales con exponente par, y de los de exponente impar, digamos $(x-r)^{2\beta +1}$ obtenemos una desigualdad equivalente si los remplazamos por exponente $1$, pues $(x-r)^{2\beta}$ es positivo y por lo tanto no cambia el signo de la desigualdad si lo ignoramos.

En resumen, cuando estamos resolviendo una desigualdad del estilo $p(x)>0$ podemos, sin cambiar el conjunto solución, reducirla a una de la forma $$q(x):=a(x-r_1)(x-r_2)\ldots(x-r_m)>0.$$ La observación clave para resolver desigualdades de este estilo está resumida en el siguiente resultado.

Proposición. Tomemos un polinomio $q(x)$ en $\mathbb{R}[x]$ de la forma $$q(x)=a(x-r_1)(x-r_2)\ldots(x-r_m)$$ con $r_1<\ldots<r_m$ reales.

Si $m$ es par:

  • Para reales $r$ en la unión de intervalos $$(-\infty,r_1)\cup(r_2,r_3)\cup\ldots \cup (r_{m-2},r_{m-1})\cup (r_m,\infty),$$ la evaluación $q(r)$ tiene el mismo signo que $a$
  • Para reales $r$ en la unión de intervalos $$(r_1,r_2)\cup(r_3,r_4)\cup\ldots \cup (r_{m-3},r_{m-2})\cup (r_{m-1},r_m),$$ la evaluación $q(r)$ tiene signo distinto al de $a$.

Si $m$ es impar:

  • Para reales $r$ en la unión de intervalos $$(r_1,r_2)\cup(r_3,r_4)\cup\ldots \cup (r_{m-2},r_{m-1})\cup (r_m,\infty),$$ la evaluación $q(r)$ tiene el mismo signo que $a$.
  • Para reales $r$ en la unión de intervalos $$(-\infty,r_1)\cup(r_2,r_3)\cup\ldots \cup (r_{m-3},r_{m-2})\cup (r_{m-1},r_m),$$ la evaluación $q(r)$ tiene signo distinto al de $a$.

Demostración. El producto $(r-r_1)(r-r_2)\ldots(r-r_m)$ es positivo si y sólo si tiene una cantidad par de factores negativos. Si $r>r_m$, todos los factores son positivos, y por lo tanto $q(r)$ tiene el mismo signo que $a$ cuando $r$ está en el intervalo $(r_m,\infty)$.

Cada que movemos $r$ de derecha a izquierda y cruzamos un valor $r_i$, cambia el signo de exactamente uno de los factores, y por lo tanto la paridad de la cantidad de factores negativos. El resultado se sigue de hacer el análisis de casos correspondiente.

$\square$

Veamos cómo podemos utilizar esta técnica para resolver desigualdades polinomiales que involucran a un polinomio que ya está factorizado en irreducibles.

Problema 1. Determina para qué valores reales $x$ se tiene que $$-2(x-5)^7(x+8)^4(x+2)^3(x+10)(x^2-x+2)^3$$ es positivo.

Solución. Por la discusión anterior, podemos ignorar el polinomio cuadrático del final, pues es irreducible. También podemos ignorar los factores lineales con potencia par, y podemos remplazar las potencias impares por unos. Así, basta con encontrar los valores reales de $x$ para los cuales $$q(x)=-2(x-5)(x+2)(x+10)$$ es positivo. Tenemos $3$ factores, así que estamos en el caso de $m$ impar en la proposición.

Las tres raíces, en orden, son $-10, -2, 5$. Por la proposición, para $x$ en la unión de intervalos $$(-\infty,-10)\cup (-2,5)$$ se tiene que $q(x)$ tiene signo distinto al de $a=-2$ y por lo tanto es positivo. Para $x$ en el conjunto $$(-10,-2)\cup (5,\infty)$$ se tiene que $q(x)$ tiene signo igual al de $a=-2$, y por lo tanto es negativo. De esta forma, la respuesta es el conjunto $$(-\infty,-10)\cup (-2,5).$$

Puedes dar clic aquí para ver en GeoGebra las gráfica de $q(x)$ y del polinomio original, y verificar que tienen el mismo signo en los mismos intervalos.

$\triangle$

Si estamos resolviendo una desigualdad y el valor de $a$ en la factorización es positivo, es un poco más práctico ignorarlo desde el principio, pues no afecta a la desigualdad.

Problema 2. Determina para qué valores reales $x$ se tiene que $$7(x+7)^{13}(x+2)^{31}(x-5)^{18}(x^2+1)$$ es positivo.

Solución. Tras las cancelaciones correspondientes, obtenemos la desigualdad equivalente $$(x+7)(x+2)>0.$$

Las raíces del polinomio que aparece son $-7$ y $-2$. De acuerdo a la proposición, estamos en el caso con $m$ par. De esta forma, la expresión es negativa en el intervalo $(-7,-2)$ y es positiva en la unión de intervalos $$(-\infty,-7)\cup (-2,\infty).$$

$\triangle$

Otras desigualdades de polinomios y manipulaciones algebraicas

Si tenemos otras expresiones polinomiales, también podemos resolverlas con ideas similares, solo que a veces se tienen que hacer algunas manipulaciones previas para llevar la desigualdad a una de la forma $p(x)>0$.

Problema. Determina todos los números $x$ en $\mathbb{R}$ para los cuales $$x^6-12x^4-49x^2-30 > 3x^5-48x^3-51x+6.$$

Solución. El problema es equivalente a encontrar los reales $x$ para los cuales $$x^6-3x^5+12x^4+48x^3-29x^2+51x-36>0.$$ El polinomio del lado izquierdo se puede factorizar como $(x-3)^2(x-1)(x+4)(x^2+1)$, así que obtenemos el problema equivalente $$(x-3)^2(x-1)(x+4)(x^2+1)>0,$$ que ya sabemos resolver. El resto de la solución queda como tarea moral.

Puedes ver la gráfica del polinomio $$(x-3)^2(x-1)(x+4)(x^2+1)$$ en GeoGebra si das clic aquí.

$\triangle$

Tener cuidado al multiplicar por denominadores

Hay que tener cuidado al realizar algunas manipulaciones algebraicas, pues pueden cambiar el signo de la desigualdad que estamos estudiando. Veamos un ejemplo donde sucede esto.

Problema. Determina todos los números $x$ en $\mathbb{R}$ para los cuales $$\frac{1}{x}>x^3-x^2+1.$$

Solución. La expresión no está definida en $x=0$, pues se anula un denominador. Supongamos entonces que $x\neq 0$, y recordémoslo al expresar la solución final. Vamos a multiplicar la desigualdad por $x$, pero tenemos que hacer casos.

Si $x>0$, entonces el signo de desigualdad no se altera y obtenemos la desigualdad equivalente $$0>x^4-x^3+x-1=(x-1)(x+1)(x^2-x+1).$$ El factor cuadrático es irreducible y lo podemos ignorar. Si estuviéramos trabajando en todo $\mathbb{R}$, el conjunto solución sería el intervalo $(-1,1)$. Sin embargo, tenemos que restringir este conjunto solución sólo al caso en el que estamos, es decir, $x>0$. Así, para este caso sólo los reales en $(0,1)$ son solución.

Si $x<0$, entonces el signo de la desigualdad sí se altera, y entonces obtenemos la desigualdad equivalente $$0<x^4-x^3+x-1=(x-1)(x+1)(x^2-x+1).$$ De nuevo podemos ignorar el factor cuadrático. La desigualdad tiene solución en todo $\mathbb{R}$ al conjunto $(-\infty,-1)\cup (1,\infty)$, pero en este caso debemos limitarlo adicionalmente con la restricción $x<0$. De este modo, las soluciones para este caso están en el intervalo $(-\infty,-1)$.

Ahora sí, juntando ambos casos, tenemos que el conjunto solución final es $$(-\infty,-1)\cup(0,1).$$

Puedes ver la gráfica en GeoGebra de $\frac{1}{x}-x^3+x^2-1$ dando clic aquí. Ahí puedes verificar que esta expresión es positiva exactamente en el conjunto que encontramos.

$\triangle$

Más adelante…

Como queda claro, resulta ser útil tener un polinomio en su forma factorizada para resolver desigualdades de polinomios reales. En los ejemplos que dimos en esta entrada, se dieron las factorizaciones de los polinomios involucrados. En el resto del curso veremos herramientas que nos permitirán encontrar la factorización de un polinomio o, lo que es parecido, encontrar sus raíces:

  • Veremos propiedades de continuidad de polinomios para mostrar la existencia de raíces para polinomios reales en ciertos intervalos.
  • El teorema del factor nos dice que si $r$ es raíz de $p(x)$, entonces $x-r$ divide a $p(x)$. Sin embargo, no nos dice cuál es la multiplicidad de $r$. Veremos que la derivada de un polinomio nos puede ayudar a determinar eso.
  • También veremos el criterio de la raíz racional, que nos permite enlistar todos los cantidatos a ser raíces racionales de un polinomio $p(x)$ con coeficientes racionales.
  • Finalmente, veremos que para los polinomios de grado $3$ y $4$ hay formas de obtener sus raíces de forma explícita, mediante las fórmulas de Cardano y de Ferrari.

Tarea moral

A continuación hay algunos ejercicios para que practiques los conceptos vistos en esta entrada. Te será de mucha utilidad intentarlos para entender más la teoría vista.

  1. Completa la solución del problema enunciado en la sección de manipulaciones algebraicas.
  2. Encuentra el conjunto solución de números reales $x$ tales que $$(x+1)(x+2)^2(x+3)^3(x+4)^4>0.$$
  3. Determina las soluciones reales a la desigualdad $$\frac{x-1}{x+2}>\frac{x+2}{x-1}.$$ Ten cuidado con los signos. Verifica tu respuesta en este enlace de GeoGebra, que muestra la gráfica de $f(x)=\frac{x-1}{x+2}-\frac{x+2}{x-1}$.
  4. Realiza las gráficas de otros polinomios de la entrada en GeoGebra para verificar las soluciones dadas a las desigualdades de polinomios.
  5. Revisa esta entrada, en donde se hablan de aplicaciones de desigualdades polinomiales para un problema de un concurso de matemáticas.

Entradas relacionadas

Agradecimientos

Trabajo realizado con el apoyo del Programa UNAM-DGAPA-PAPIME PE104522 «Hacia una modalidad a distancia de la Licenciatura en Matemáticas de la FC-UNAM – Etapa 2»