Cálculo Diferencial e Integral II: Integrales impropias del $2^{do}$ tipo

Por Miguel Ángel Rodríguez García

Introducción

En la sección anterior vimos las integrales impropias del primer tipo en el cual son integrales de la forma:

$$\int_{-\infty}^{b}f(x)dx$$

$$\int_{a}^{\infty}f(x)dx$$

$$\int_{-\infty}^{\infty}f(x)dx$$

En esta sección veremos las integrales impropias del segundo tipo que son integrales donde la función $f(x)$ no está definida en todo el intervalo, en algún punto dentro de un intervalo o en los extremos del intervalo.

Integrales impropias del $2^{do}$ tipo

Las integrales impropias del segundo tipo son integrales de la forma:

$$\int_{a}^{b}f(x)dx$$

Donde $f(x)$ es una función positiva y continua en un intervalo finito $[a ,b)$, por lo que se tiene una asíntota vertical en $b$, es decir, la función $f(x)$ tiene una discontinuidad en $x=b$.

Análogamente, si tenemos el intervalo $(a, b]$ se tiene una asíntota vertical en $a$, es decir, la función $f(x)$ tiene una discontinuidad en $x=a$.

Si tenemos que la función se indefine en un punto o varios puntos dentro del intervalo $[a, b]$ con $a<c<b$, entonces las integrales, que también son integrales impropias del segundo tipo, son integrales de funciones que se vuelven infinitas en un punto o varios puntos dentro de un intervalo, para solucionar este problema, veamos la siguiente definición:

Definición: Sea $f(x)$ una función continua en un intervalo $[a, b)$ y con una discontinuidad infinita en $b$ entonces definimos:

$$\int_{a}^{b}f(x)dx=\lim_{x \to b^{-}}\int_{a}^{x}f(t)dt \tag{1}$$

Definición: Análogamente, al caso anterior, sea $f(x)$ una función continua en un intervalo $(a, b]$ y con una discontinuidad infinita en $a$ entonces definimos:

$$\int_{a}^{b}f(x)dx=\lim_{x \to a^{+}}\int_{x}^{b}f(t)dt \tag{2}$$

Definición: Si $f(x)$ es continua en $[a, c) \cup (c,b]$, con $a<c<b$, entonces definimos:

$$\int_{a}^{b}f(x)dx=\int_{a}^{c}f(x)dx+\int_{c}^{b}f(x)dx \tag{2}$$

En estos casos, si la integral de la función $f(x)$ converge a un valor $L$ entonces se dice que la integral es convergente y converge al valor $L$, en caso contrario, se dice que la integral de la función $f(x)$ diverge.

Veamos unos ejemplos.

Ejemplos

  • $$\int_{0}^{1}\frac{1}{\sqrt[3]{x}}dx$$

Vemos que la función es discontinua en $x=0$, por la definición $(2)$, podemos calcular la integral como:

$$\int_{0}^{1}\frac{1}{\sqrt[3]{x}}dx=\lim_{x \to 0^{+}}\int_{x}^{1}\frac{1}{\sqrt[3]{t}}dt=\lim_{x \to 0^{+}}\int_{x}^{1}t^{-1/3}dt$$

$$=\lim_{x \to 0^{+}}(\frac{3}{2})t^{\frac{2}{3}}\bigg|_{x}^{1}=\lim_{x \to 0^{+}}(\frac{3}{2}-\frac{3}{2}x^{\frac{2}{3}})=\frac{3}{2}$$

$$\therefore \int_{0}^{1}\frac{1}{\sqrt[3]{x}}dx \space \space converge \space a \space \space \frac{3}{2}$$

  • $$\int_{2}^{5}\frac{1}{\sqrt{x-2}}dx$$

Vemos que la función es discontinua en $x=2$, por la definición $(2)$, calculamos la integral como:

$$\int_{2}^{5}\frac{1}{\sqrt{x-2}}dx=\lim_{x \to 2^{+}}\int_{x}^{5}\frac{1}{\sqrt{t-2}}dt=\lim_{x \to 2^{+}}2\sqrt{t-2}\bigg|_{x}^{5}=\lim_{x \to 2^{+}}2(\sqrt{3}-\sqrt{x-2})=2\sqrt{3}$$

$$\therefore \int_{2}^{5}\frac{1}{\sqrt{x-2}}dx=2\sqrt{3}$$

Las integrales impropias del tercer tipo son integrales que mezcla los dos tipos anteriores de integrales impropias, es decir, son integrales impropias tanto del primer tipo como del segundo tipo, por lo que para resolver este tipo de integrales del tercer tipo, se utiliza las mismas estrategias para resolver las integrales del primer y segundo tipo.

Ejemplo

  • $$\int_{-\infty}^{\infty}\frac{1}{x^{2}}dx$$

Vemos en los límites de integración que es una integral de primer tipo, pero la función tiene una discontinuidad en $x=0$, por lo que esta integral es de tercer tipo, determinamos su convergencia con las definiciones de las integrales de primer y segundo tipo como sigue:

$$\int_{-\infty}^{\infty}\frac{1}{x^{2}}dx=\lim_{x \to -\infty}\int_{x}^{0}\frac{1}{t^{2}}dt+\lim_{x \to \infty}\int_{0}^{\infty}\frac{1}{t^{2}}dt$$

$$=\lim_{x \to -\infty}\left ( \lim_{t \to 0^{-}} \int_{x}^{t} \frac{1}{v^{2}}dv \right )+\lim_{x \to \infty}\left ( \lim_{t \to 0^{+}} \int_{t}^{x} \frac{1}{v^{2}}dv \right )$$

$$=\lim_{x \to -\infty}\left ( \lim_{t \to 0^{-}} \frac{-1}{v}\bigg{|}^{t}_{x} \right)+\lim_{x \to \infty}\left ( \lim_{t \to 0^{+}} \frac{-1}{v}\bigg{|}^{t}_{x} \right )$$

$$=\lim_{x \to -\infty}\left ( \lim_{t \to 0^{-}} \left ( \frac{-1}{t}+\frac{1}{x} \right ) \right)+\lim_{x \to \infty}\left ( \lim_{t \to 0^{+}} \left ( \frac{-1}{t}+\frac{1}{x} \right ) \right )$$

Sabemos que $\lim_{t \to 0} \frac{1}{t}$ se indetermina.

$$\therefore \int_{-\infty}^{\infty}\frac{1}{x^{2}}dx \space \space diverge$$

Tarea moral

Los siguientes ejercicios no son para evaluación, pero son ejercicios para que practiques lo aprendido que te ayudaran en el desarrollo del entendimiento del tema, por lo que te invitamos a resolver los siguientes ejercicios propuestos relacionados con el tema visto.

Determine la convergencia de las siguientes integrales:

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. $$\int_{0}^{2}\frac{1}{x^{3}}dx$$
  2. $$\int_{0}^{\infty}\frac{dx}{\sqrt{x}(x+1)}$$
  3. $$\int_{0}^{1}\frac{1}{1-x}dx$$
  4. $$\int_{0}^{3}\frac{1}{(x-1)^{2/3}}dx$$
  5. $$\int_{0}^{1}ln(x)dx$$

Más adelante…

En esta sección vimos las integrales del segundo tipo y tercer tipo, así, como algunos ejemplos, por lo que terminamos de ver los tipos de integrales impropias, en la siguiente sección veremos algunos criterios de convergencia importante para las integrales impropias.

Entradas relacionadas

Cálculo Diferencial e Integral II: Integrales impropias del $1^{er}$ tipo

Por Miguel Ángel Rodríguez García

Introducción

En la sección anterior vimos el teorema del valor medio para integrales, en esta sección veremos las integrales impropias de primer tipo.

Al introducir el concepto de integral definida se exigió que las funciones estuvieran definidas en intervalos cerrados y que la integral de esas funciones en ese intervalo este definida. En esta entrada se suprimen esas restricciones y veremos integrales del tipo:

$$\int_{a}^{\infty }f(x)dx$$

$$\int_{0}^{1}\frac{dx}{\sqrt{1-x}}$$

Obsérvese que en la primera integral el límite de integración se escribe el símbolo de infinito y en la segunda integral para el punto $x=1$ el integrando no está definido en $1$, por lo que veremos las definiciones siguientes.

Integrales impropias del $1^{er}$ tipo

Definición. Sea $f$ una función continua definida en $[a, \infty)$ entonces definimos:

$$\int_{a}^{\infty }f(x)dx=\lim_{x \to \infty }\int_{a}^{x}f(t)dt \tag{1}$$

A $\int_{a}^{\infty }f(x)dx$ se le llama la integral impropia del $1^{er}$ tipo de la función $f$ de $a$ hasta $\infty$.

Definición. Si $\lim_{x \to \infty }\int_{a}^{x}f(t)dt$ es un numero real $L$ se dice que la integral es convergente y converge al valor $L$.

En cambio, si $\lim_{x \to \infty }\int_{a}^{x}f(t)dt$ da como resultado $\infty$ o $-\infty$, es decir, la integral diverge, entonces se dice que la integral diverge a $\infty$ o $-\infty$.

Análogamente, se puede dar la misma definición para cuando el límite de integración inferior tiende a $- \infty$.

Definición. Sea $f$ una función continua definida en $(-\infty, b]$ entonces definimos:

$$\int_{-\infty}^{ b}f(x)dx=\lim_{x \to -\infty }\int_{x}^{b}f(t)dt \tag{2}$$

Podemos tener integrales impropias de una función $f(x)$, tal que, los límites de integración van de $-\infty$ a $\infty$, en este caso, definimos lo siguiente:

Definición. Sea una función continua en $(-\infty, \infty)$ entonces:

$$\int_{-\infty}^{\infty }f(x)dx=\int_{-\infty}^{a}f(x)dx+\int_{a}^{\infty}f(x)dx \tag{3}$$

Veamos unos ejemplos.

Ejemplos

  • Calcula, si es posible, la integral $\int_{1}^{\infty }\frac{dx}{x^{2}}$.

Usamos la definición $(1)$, así:

$$\int_{1}^{\infty }\frac{dx}{x^{2}}=\lim_{x \to \infty }\int_{1}^{x}\frac{1}{t^{2}}dt$$

$$\Rightarrow =\lim_{x \to \infty }\int_{1}^{x}{t^{-2}}dt=\lim_{x \to \infty }\left [ (-1)t^{-1} \right ] \bigg|_{1}^{x}=\lim_{x \to \infty }(-x^{-1}-(-1)^{-1})=-\lim_{x \to \infty }\frac{1}{x}+\lim_{x \to \infty }1$$

Sabemos que:

$$\lim_{x \to \infty }\frac{1}{x}=0$$

Entonces:

$$\int_{1}^{\infty }\frac{dx}{x^{2}}=1$$

  • Calcula la siguiente integral impropia $\int_{-\infty}^{\infty} e^{x-e^{x}}dx$.

Por definición $(3)$, se tiene que:

$$\int_{-\infty}^{\infty} e^{x-e^{x}}dx=\int_{-\infty}^{a} e^{x-e^{x}}dx+\int_{\infty}^{a} e^{x-e^{x}}dx$$

Usamos ahora las definiciones $(1)$ y $(2)$ como:

$$=\lim_{x \to -\infty} \int_{x}^{a}e^{t-e^{t}}dt+\lim_{x \to \infty} \int_{a}^{x}e^{t-e^{t}}dt$$

Para integrar esta función solo utilizamos el método de cambio de variable, para esto, sea $u=e^{t}$, entonces:

$$\int e^{t-e^{t}}dt=\int e^{-u}dt=-e^{-u}=-e^{-e^{t}}$$

Así, la integral impropia se resuelve como:

$$\lim_{x \to -\infty} \int_{x}^{a}e^{t-e^{t}}dt+\lim_{x \to \infty} \int_{a}^{x}e^{t-e^{t}}dt=\lim_{x \to -\infty}-e^{-e^{x}}+\lim_{x \to \infty}-e^{-e^{x}}=1+0=1$$

Por tanto, la integral converge a 1.

Veamos el teorema siguiente que nos dice para que casos la función $ \frac{dx}{x^{s}} $ converge:

Teorema: $\forall \space s >1$ la integral:

$$\int_{1}^{\infty } \frac{dx}{x^{s}}=\frac{1}{s-1}$$

Es decir, $\frac{1}{x^{s}}$ converge. Sin embargo, si $s\leq 1$ la integral: $$\int_{1}^{\infty } \frac{dx}{x^{s}} \space diverge \space a \space \infty$$

Demostración:

Veamos la demostración por casos.

Sea $s\neq 1$, entonces por definición $(1)$ se tiene que:

$$\int_{1}^{\infty } \frac{dx}{x^{s}}=\lim_{x \to \infty}\int_{1}^{x} \frac{dt}{t^{s}}=\lim_{x \to \infty}\frac{t^{-s+1}}{-s+1}\bigg|_{1}^{x}= \lim_{x \to \infty}\left ( \frac{x^{-s+1}}{-s+1}-\frac{1^{-s+1}}{-s+1} \right )$$

  • Si $s>1 \Rightarrow -s<-1 \Rightarrow -s+1<0$, entonces tenemos que:

$$\lim_{x \to \infty}\frac{x^{-s+1}}{-s+1}=0$$

$$\therefore \int_{1}^{\infty}\frac{dx}{x^{s}}=\lim_{x \to \infty}\left (-\frac{1^{-s+1}}{-s+1} \right )=\frac{-1}{-s+1}=\frac{1}{s-1}$$

  • Si $s<1 \Rightarrow -s>-1 \Rightarrow -s+1>0$, entonces:

$$\lim_{x \to \infty}\frac{x^{-s+1}}{-s+1} \to \infty$$

$$\therefore \int_{1}^{\infty}\frac{dx}{x^{s}} \space Diverge$$

  • Si $s=1$, entonces:

$$\int_{1}^{\infty}\frac{dx}{x^{s}}=\int_{1}^{\infty}\frac{dx}{x}=\lim_{x \to \infty}\int_{1}^{\infty}\frac{dt}{t}=\lim_{x \to \infty}\left ( ln(t) \right )\bigg|_{1}^{x}=\lim_{x \to \infty}(ln(x)-ln(1))=\lim_{x \to \infty}ln(x)-0\to \infty$$

$$\therefore \int_{1}^{\infty}\frac{dx}{x^{s}} \space Diverge$$

$$\therefore \int_{1}^{\infty } \frac{dx}{x^{s}}=\frac{1}{s-1}$$

Converge para $s>1$ y diverge para $s\leq 1$.

$\square$

Veamos el siguiente ejemplo.

Ejemplo

  • $$\int_{1}^{\infty}\frac{dx}{x^{3/2}}$$

Vemos que del integrando podemos usar el teorema visto anteriormente donde $s=\frac{3}{2}>1$ por lo que, en ese caso, tenemos que:

$$\int_{1}^{\infty}\frac{dx}{x^{3/2}}=\frac{1}{\frac{3}{2}-1}=\frac{1}{\frac{1}{2}}=2$$

Tarea moral

Los siguientes ejercicios no son para evaluación, pero son ejercicios para que practiques lo aprendido que te ayudaran en el desarrollo del entendimiento del tema, por lo que te invitamos a resolver los siguientes ejercicios propuestos relacionados con el tema visto.

Calcule las siguiente integrales.

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. $$\int_{0}^{\infty}e^{-x}dx$$
  2. $$\int_{-\infty}^{0}sin(x)dx$$
  3. $$\int_{-\infty}^{\infty}\frac{dx}{x^{2}+1}$$
  4. $$\int_{0}^{\infty}\frac{e^{x}}{1+e^{2x}}dx$$
  5. $$\int_{1}^{\infty}(1-x)e^{-x}dx$$

Más adelante…

En esta sección vimos integrales impropias del $1^{er}$ tipo que son integrales en donde se integra en un intervalo infinito y se necesita saber el área bajo la curva de una función $f(x)$, es decir, en intervalos no acotados, en la siguiente sección veremos integrales impropias del $2^{do}$ tipo que son integrales impropias en donde la discontinuidad de la función $f(x)$ no está definida en algún punto o todo intervalo en $(a,b)$.

Entradas relacionadas

Cálculo Diferencial e Integral II: Teorema del valor medio para integrales

Por Miguel Ángel Rodríguez García

Introducción

En las secciones anteriores vimos algunos métodos numéricos de integración que se utilizan para dar solución a la integral de funciones, en esta sección veremos el teorema del valor medio para integrales.

Teorema del valor medio para integrales

El teorema del valor medio es una consecuencia del teorema de valor medio para la derivada y el teorema fundamental del Cálculo [Hipervinculo: Calculo II-Teorema fundamental del calculo], geométricamente significa que para funciones no negativas y continuas en un intervalo $[a, b]$ existe un valor $c$ en el mismo intervalo, tal que, el rectángulo con base $[a, b]$ y altura $f(c)$ tiene la misma área que la región bajo la gráfica de $f$ en el intervalo $[a, b]$ o lo que es lo mismo decir, alcanza su valor promedio en al menos un punto $c$ (ver figura 1).

Figura 1: Teorema del valor medio para integrales.

Enunciamos el siguiente teorema:

Teorema del valor medio para integrales

Sea función continua $f(x)$ en un intervalo $[a, b]$ entonces existe $c \space \epsilon \space [a, b]$, tal que:

$$\int_{a}^{b}f(x)dx=f(c)(b-a)$$

Demostración:

Como $c \space \epsilon \space [a, b]$ supongamos sin perdida de generalidad que $a<c<b$ entonces:

$$\int_{a}^{b} f(x)dx=\int_{a}^{c} f(x)dx+\int_{c}^{b} f(x)dx$$

Por las propiedades de la integral [Hipervinculo: Calculo II-Propiedades de la integral], la función $f$ es integrable en $[a, c]$ y $[c, b]$.

Ahora por el teorema del valor extremo sabemos que $f$ alcanza sus valores máximo y mínimo absolutos en el intervalo, el cual se denotan como $M$ y $m$ respectivamente, así, sabemos que:

$$\int_{a}^{b} mdx\leq \int_{a}^{b} f(x)dx\leq \int_{a}^{b} Mdx$$

Las integrales de la izquierda y derecha se pueden evaluar fácilmente:

$$m(b-a)\leq \int_{a}^{b} f(x)dx\leq M(b-a) \tag{1}$$

Por otro lado, como $c \space \epsilon \space [a, b]$ entonces: $m\leq f(c) \leq M$ para alguna $c \space \epsilon \space [a, b]$.

Si $m$ y $M$ son infinitesimalmente pequeños, entonces $m=f(c)=M$, por lo que en $(1)$:

$$f(c)(b-a)\leq \int_{a}^{b} f(x)dx\leq f(c)(b-a)$$

$$\Rightarrow \int_{a}^{b}f(x)dx=f(c)(b-a)$$

$\square$

Veamos un ejemplo.

Ejemplo

  • Determine el valor promedio de la función $f(x)=1+x^{2}$ en el intervalo $[-1,2]$.

Vemos que $a=-1$ y $b=2$, para calcular el valor promedio de la función $f(x)$ utilizamos el teorema del valor medio como sigue:

$$f_{prom}=f(c)=\frac{1}{b-a}\int_{a}^{b}f(x)dx$$

Así tenemos que:

$$f(c)=\frac{1}{2-(-1)}\int_{-1}^{2}(1+x^{2})dx=\frac{1}{3}\left [ x+\frac{x^{3}}{3} \right ]\bigg{|}_{-1}^{2}$$

$$=\frac{1}{3}\left [ 2+\frac{2^{3}}{3}-(-1)-\frac{(-1)^{3}}{3} \right ]=2$$

Vemos que $f(c)=2$, evaluamos en la función $f(x)$ el valor $c$ para encontrar su valor:

$$f(c)=1+c^{2}=2 \Rightarrow c=\pm 1$$

Sucede que en este caso hay dos números $c=1$ y $c=-1$ que toman el valor medio de la función $f(x)$.

Tarea moral

Los siguientes ejercicios no son para evaluación, pero son ejercicios para que practiques lo aprendido que te ayudaran en el desarrollo del entendimiento del tema, por lo que te invitamos a resolver los siguientes ejercicios propuestos relacionados con el tema visto.

Halle el valor promedio de las siguientes funciones en el intervalo indicado.

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. $f(x)=4x-x^{2}$, $[0,4]$
  2. $f(x)=4-x$, $[0,3]$
  3. $f(x)=3x^{2}-2x$, $[1,4]$
  • Determine los números $b$ talque el valor promedio de la función $f(x)=2+6x-3x^{2}$ en el intervalo $[0,b]$ sea igual a 3.
  • Demuestre que la velocidad promedio de un automóvil en un intervalo de tiempo $[t_{1}, t_{2}]$ es la misma que el promedio de sus velocidades.

Más adelante…

En esta sección vimos el teorema del valor intermedio aplicado a las integrales, en las siguientes secciones veremos las integrales impropias, es decir, integrales en donde se evalúa una función dentro de un intervalo que tiende a infinito o casos en donde la integral de una función se evalúa en todo $\mathbb{R}$.

Entradas relacionadas

Cálculo Diferencial e Integral II: Métodos Numéricos de Integración – Regla de Simpson

Por Miguel Ángel Rodríguez García

Introducción

En la sección anterior vimos dos métodos numéricos de integración: el método del punto medio y el método del trapecio. Otra regla de aproximación numérica a las integrales se llama regla de Simpson, el cual consiste en usar parábolas (como se muestra en la figura $1$) en lugar de segmentos de rectas para aproximarse a una curva.

Método de la regla de Simpson

Comencemos deduciendo la regla de Simpson.

Sea una curva dada por $f(x)$ en el plano en un intervalo $[a, b]$, dividimos el intervalo $[a, b]$ en $n$ subintervalos de igual longitud dado como:

$$\Delta x=\frac{b-a}{n}$$

En el que esta vez se requiere que $n$ sea un número par.

Figura 1: Regla de Simpson con aproximaciones parabólicas a la función $f(x)$.

La ecuación de una parábola está dada como:

$$y=Ax^{2}+Bx+C \tag{1}$$

Por lo que su área en el intervalo $[-h, h]$ es:

$$Área=\int_{-h}^{h}(Ax^{2}+Bx+C)dx=\left [ A\frac{x^{3}}{3}+B\frac{x^{2}}{2}+Cx+D \right ]\bigg{|}_{-h}^{h}$$

$$=\left [ A\frac{h^{3}}{3}+B\frac{h^{2}}{2}+Ch+D \right ]- \left [ A\frac{(-h)^{3}}{3}+B\frac{ (-h) ^{2}}{2}+C (-h) +D \right ]$$

$$=\left [ A\frac{h^{3}}{3}+B\frac{h^{2}}{2}+Ch+D \right ]+\left [ A\frac{h^{3}}{3}-B\frac{h^{2}}{2}+Ch-D \right ] $$

$$=\frac{2Ah^{3}}{3}+2Ch=h\frac{(2Ah^{2}+6C)}{3} \tag{2}$$

De la figura $1$ vemos que una de las curvas pasa por los puntos $(-h, y_{0})$, $(0, y_{1})$ y $(h, y_{2})$, evaluando estos puntos en la ecuación cuadrática $(1)$ se obtiene lo siguiente:

$$y_{0}=Ah^{2}-Bh+C$$

$$y_{1}=C$$

$$y_{2}=Ah^{2}+Bh+C$$

Si sumamos estas relaciones como:

$$y_{0}+4y_{1}+y_{2}= Ah^{2}-Bh+C +4C+ Ah^{2}+Bh+C= 2Ah^{2}+6C $$

Podemos expresar el área $(2)$ en términos de $y_{0}$, $y_{1}$ y $y_{2}$, como:

$$A_{1}=\frac{h}{3}(y_{0}+4y_{1}+y_{2})$$

Que es el área debajo de la parábola que pasa por los puntos $(x_{0}=-h, y_{0})$, $(x_{1}=0, y_{1})$ y $(x_{2}=h, y_{2})$, imaginemos que la segunda parábola intercepta en los puntos: $(x_{2}, y_{2})$, $(x_{3}, y_{3})$ y $(x_{4}, y_{4})$ entonces el área de esta segunda parábola es:

$$A_{2}=\frac{h}{3}(y_{2}+4y_{3}+y_{4})$$

Si sumamos todas las áreas hasta un n-esima parábola que se aproxima a la función $f(x)$, tendremos que el área total es:

$$\int_{a}^{b}f(x)dx\approx S_{n}=\frac{h}{3}(y_{0}+4y_{1}+y_{2})+\frac{h}{3}(y_{2}+4y_{3}+y_{4})+…+\frac{h}{3}(y_{n-2}+4y_{n-1}+y_{n})$$

$$=\frac{h}{3}(y_{0}+4y_{1}+2y_{2}+4y_{3}+…+2y_{n-2}+4y_{n-1}+y_{n})$$

Vemos que hay un patrón en los coeficientes:

$$1, \space 4, \space2, \space4, \space2 \space…. \space 2, \space4, \space1$$

Por lo que la regla de Simpson se define como:

$$\int_{a}^{b}f(x)dx\approx \frac{\Delta x}{3}(y_{0}+4y_{1}+y_{2})+\frac{\Delta x}{3}(y_{2}+4y_{3}+y_{4})+…+\frac{\Delta x}{3}(y_{n-2}+4y_{n-1}+y_{n})$$

$$=\frac{\Delta x}{3}(y_{0}+4y_{1}+2y_{2}+4y_{3}+…+2y_{n-2}+4y_{n-1}+y_{n}) \tag{3}$$

Con $\Delta x=\frac{b-a}{n}$, $n$ un número par, y los puntos $x_{i}$ los calculamos como:

$$x_{0}=a$$

$$x_{1}=a+\Delta x$$

$$…..$$

$$x_{n-1}=a+(n-1)\Delta x$$

$$x_{n}=b \tag{4}$$

Cota de error para la regla de Simpson

Para la estimación de la cota de error en la regla de Simpson, suponga que $|f^{4}(x)|\leq K$ para $a\leq x\leq b$ con $|f^{4}(x)|$ el valor absoluto de la cuarta derivada de la función. Si $E_{s}$ es el error relacionado con la regla de Simpson, entonces la cota de error para la regla de Simpson es:

$$E_{s}\leq\frac{K(b-a)^{5}}{180n^{4}}$$

Veamos un ejemplo.

Ejemplo

  • Usar la regla de Simpson para aproximar la integral $\int_{1}^{2}\frac{1}{x}dx$ con $n=10$.

Tenemos que $n=10$, $a=1$ y $b=2$ lo que implica que $\Delta x=\frac{b-a}{n}=0.1$.

Por la regla de Simpson $(3)$ y calculando los puntos $x_{i}$ $(4)$ tenemos que:

$$\int_{1}^{2}\frac{1}{x}dx\approx\ S_{10}=\frac{\Delta x}{3}\left [ f(1)+4f(1.1)+2f(1.2)+…+2f(1.8)+4f(1.9)+f(2) \right ]$$

$$=\frac{0.1}{3}\left [ \frac{1}{1}+\frac{4}{1.1}+\frac{2}{1.2}+\frac{4}{1.3}+\frac{2}{1.4}+\frac{4}{1.5}+\frac{2}{1.6}+\frac{4}{1.7}+\frac{2}{1.8}+\frac{4}{1.9}+\frac{1}{2}+ \right ]\approx 0.693150$$

Comparando este resultado con lo obtenido con la regla del punto medio y regla del trapecio, la regla de Simpson nos da una aproximación mucho mejor respecto a estos dos métodos, pues resulta que la regla de Simpson son promedios ponderados de la regla del punto medio y regla del trapecio, se puede demostrar que:

$$S_{2n}=\frac{1}{3}T_{n}+\frac{2}{3}M_{n}$$

Tarea moral

Los siguientes ejercicios no son para evaluación, pero son ejercicios para que practiques lo aprendido que te ayudaran en el desarrollo del entendimiento del tema, por lo que te invito a resolver los siguientes ejercicios propuestos relacionados con el tema visto.

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. Demuestre que: $S_{2n}=\frac{1}{3}T_{n}+\frac{2}{3}M_{n}$
  2. ¿Qué tan grande debe de ser n para que al utiliza la regla de Simpson al aproximar la integral $\int_{1}^{2}\frac{1}{x}dx$, sea exacta hasta dentro de 0.0001?
    1. Use la regla de Simpson con n=10 para aproximar la integral $\int_{0}^{1} e^{x^{2}}dx$
    2. Estime el error con esta aproximación
    1. Estimar la integral con n=3: $\int_{0}^{2} x^{3}dx$
    2. En este caso, ¿la regla de Simpson es exacta? ¿Porque?

Más adelante…

En esta sección vimos la regla de Simpson que consiste en otro método de aproximación numérica para las integrales por medio de parábolas y que es este método es un promedio ponderado de los métodos del punto medio y del trapecio. Aunque existen más métodos numéricos para aproximar integrales, solo veremos estos métodos. En la siguiente sección veremos el teorema del valor medio para las integrales.

Entradas relacionadas

Ecuaciones Diferenciales I: Sistemas de ecuaciones diferenciales

Por Omar González Franco

El conocimiento de las matemáticas añade vigor a la mente,
la libera del prejuicio, credulidad y superstición.
– John Arbuthnot

Introducción

¡Bienvenidos a la tercera unidad del curso de Ecuaciones Diferenciales I!.

En esta unidad estudiaremos los sistemas de ecuaciones diferenciales lineales de primer orden.

En la unidad 1 de este curso estudiamos el sistema Depredador – Presa, en nuestro análisis el modelo matemático determinado fue el siguiente sistema de ecuaciones diferenciales.

\begin{align*}
\dfrac{dC}{dt} &= aC(t) -bC(t)Z(t) \\
\dfrac{dZ}{dt} &= -cZ(t) + dC(t)Z(t)
\end{align*}

Puedes revisar la entrada correspondiente para recordar que representa cada una de las variables y constantes.

Este sistema fue nuestro primer ejemplo de un sistema de ecuaciones diferenciales y en esta unidad nuestro propósito será desarrollar distintos métodos que nos permitan resolver sistemas de hasta $n > 2$ ecuaciones diferenciales acopladas.

Es importante mencionar que a lo largo de esta unidad usaremos un enfoque matricial, por lo que es recomendable tener presente, al menos, la teoría básica sobre matrices y sus operaciones y propiedades vistas en el curso de Álgebra Lineal I.

En esta entrada comenzaremos por definir los que es un sistema de ecuaciones diferenciales, sus propiedades y veremos cómo es que la notación matricial nos puede ayudar.

¡Comencemos!

Sistemas de ecuaciones diferenciales lineales de primer orden

En esta unidad, a menos que indiquemos lo contrario, la variable independiente se denotará por $t$, mientras que las variables dependientes de $t$ por

$$y_{1} = y_{1}(t), \hspace{0.5cm} y_{2} = y_{2}(t), \hspace{0.5cm} \cdots, \hspace{0.5cm} y_{n} = y_{n}(t)$$

y las funciones $F_{i}$, $i = 1, 2, 3, \cdots, n$ son funciones con valores reales que dependen de las $n + 1$ variables en un intervalo $\delta$.

Notación: Para mayor comodidad, en esta unidad usaremos la notación de prima para la derivada.

$$\dfrac{dy}{dt} = y^{\prime}(t) \label{2} \tag{2}$$

Con esta notación el sistema de ecuaciones (\ref{1}) se puede escribir de la siguiente manera.

\begin{align*}
y_{1}^{\prime}(t) &= F_{1}(t, y_{1}, y_{2}, \cdots, y_{n}) \\
y_{2}^{\prime}(t) &= F_{2}(t, y_{1}, y_{2}, \cdots, y_{n}) \\
&\vdots \\
y_{n}^{\prime}(t) &= F_{n}(t, y_{1}, y_{2}, \cdots, y_{n}) \label{3} \tag{3}
\end{align*}

En el sistema lineal (\ref{5}) se supone que los coeficientes $a_{ij}(t)$, así como las funciones $g_{i}(t)$, $i, j = \{1, 2, 3, \cdots, n \}$ son continuas en un intervalo común $\delta$.

Ejemplo: El sistema de ecuaciones diferenciales

\begin{align*}
y_{1}^{\prime}(t) &= -3y_{1} + 4y_{2} -9y_{3} \\
y_{2}^{\prime}(t) &= 6y_{1} -y_{2} \\
y_{3}^{\prime}(t) &= 10y_{1} + 4y_{2} + 3y_{3}
\end{align*}

es un sistema lineal de primer orden compuesto por tres ecuaciones diferenciales lineales de primer orden cada una.

Notación: Si el sistema es de dos o tres ecuaciones diferenciales denotaremos por $x(t), y(t)$ o $x(t), y(t)$, $z(t)$ a las variables dependientes de $t$, respectivamente.

Considerando esta notación, el sistema del ejemplo anterior se puede escribir de la siguiente manera.

\begin{align*}
x^{\prime}(t) &= -3x + 4y -9z\\
y^{\prime}(t) &= 6x -y \\
z^{\prime}(t) &= 10x + 4y + 3z
\end{align*}

Problema de valores iniciales

Es posible demostrar la existencia y unicidad de soluciones de sistemas tanto lineales como no lineales (caso general) y de soluciones a sistemas lineales homogéneos y no homogéneos (casos particulares), sin embargo las demostraciones de estos teoremas suelen ser bastantes extensas y complejas para nosotros en estos momentos, ya que requieren de herramientas matemáticas que aún desconocemos. A continuación enunciamos el teorema de existencia y unicidad para el caso general y para el caso lineal homogéneo.

En este teorema la región $R$ se construye con el producto cartesiano de los intervalos abiertos en los que $t_{0} \in \delta$, $b_{1} \in \delta_{1}$, $b_{2} \in \delta_{2}$, $\cdots$, $b_{n} \in \delta_{n}$, así $(t_{0}, b_{1}, b_{2}, \cdots, b_{n}) \in R$.

Para el caso particular de sistemas lineales homogéneos, el teorema de existencia y unicidad se puede enunciar de la siguiente forma.

Como mencionamos antes, es complejo demostrar estos teoremas, sin embargo más adelante en esta unidad los retomaremos y los justificaremos. Por ahora hay que tener en cuenta que para el caso general se requiere de volver a algunos de los conceptos vistos para demostrar el teorema de existencia y unicidad de Picard – Lindelöf de la primera unidad y para los casos particulares ¡la definición de exponencial de una matriz nos ayudará a demostrarlos!.

Ahora veamos la utilidad de la notación matricial.

Sistemas lineales de primer orden en forma matricial

Daremos por hecho que se conocen las operaciones y propiedades básicas de las matrices, así como algunas propiedades de espacios vectoriales vistas en el curso de Álgebra Lineal I.

Definamos las siguientes matrices de funciones.

$$\mathbf{Y}(t) = \begin{pmatrix}
y_{1}(t) \\ y_{2}(t) \\ \vdots \\ y_{n}(t)
\end{pmatrix} \hspace{1cm} \Rightarrow \hspace{1cm} \mathbf{Y^{\prime}}(t) = \begin{pmatrix}
y_{1}^{\prime}(t) \\ y_{2}^{\prime}(t) \\ \vdots \\ y_{n}^{\prime}(t)
\end{pmatrix} $$

y

$$\mathbf{A}(t) = \begin{pmatrix}
a_{11}(t) & a_{12}(t) & \cdots & a_{1n}(t) \\
a_{21}(t) & a_{22}(t) & \cdots & a_{2n}(t) \\
\vdots & & & \vdots \\
a_{n1}(t) & a_{n2}(t) & \cdots & a_{nn}(t)
\end{pmatrix}, \hspace{1cm}
\mathbf{G}(t) = \begin{pmatrix}
g_{1}(t) \\ g_{2}(t) \\ \vdots \\ g_{n}(t)
\end{pmatrix}$$

Usando estas matrices, el sistema de ecuaciones diferenciales lineales de primer orden (\ref{5}) se puede escribir de la siguiente manera.

$$\begin{pmatrix}
y_{1}^{\prime}(t) \\ y_{2}^{\prime}(t) \\ \vdots \\ y_{n}^{\prime}(t)
\end{pmatrix} = \begin{pmatrix}
a_{11}(t) & a_{12}(t) & \cdots & a_{1n}(t) \\
a_{21}(t) & a_{22}(t) & \cdots & a_{2n}(t) \\
\vdots & & & \vdots \\
a_{n1}(t) & a_{n2}(t) & \cdots & a_{nn}(t)
\end{pmatrix} \begin{pmatrix}
y_{1}(t) \\ y_{2}(t) \\ \vdots \\ y_{n}(t)
\end{pmatrix} + \begin{pmatrix}
g_{1}(t) \\ g_{2}(t) \\ \vdots \\ g_{n}(t)
\end{pmatrix} \label{8} \tag{8}$$

o bien,

$$\mathbf{Y^{\prime}} = \mathbf{AY} + \mathbf{G} \label{9} \tag{9}$$

Si el sistema es homogéneo, entonces escribimos

$$\mathbf{Y^{\prime}} = \mathbf{AY} \label{10} \tag{10}$$

La solución de un sistema lineal la podemos definir como sigue.

Usando la notación matricial, un PVI se puede escribir de la siguiente manera.

El teorema de existencia y unicidad para el caso lineal se puede enunciar de la siguiente forma.

Verifica que el sistema de ecuaciones diferenciales usado como ejemplo al inicio de la entrada se puede escribir en notación matricial de la siguiente forma.

$$\mathbf{Y}^{\prime} = \begin{pmatrix}
-3 & 4 & -9 \\ 6 & -1 & 0 \\ 10 & 4 & 3
\end{pmatrix} \mathbf{Y}$$

Veamos un ejemplo más.

Ejemplo: Escribir el siguiente sistema lineal en forma matricial.

\begin{align*}
x^{\prime}(t) &= x -y + z + t + 1 \\
y^{\prime}(t) &= 2x + y -z -3t^{2} \\
z^{\prime}(t) &= x + y + z + t^{2} -t + 2
\end{align*}

Solución: Primero escribamos cada lado de las ecuaciones en una matriz.

$$\begin{pmatrix}
x^{\prime}(t) \\ y^{\prime}(t) \\ z^{\prime}(t)
\end{pmatrix} = \begin{pmatrix}
x -y + z + t -1 \\ 2x + y -z -3t^{2} \\ x + y + z + t^{2} -t + 2
\end{pmatrix}$$

La matriz derecha la separamos en dos, una que contenga a las variables dependientes y otra a la variable independiente.

$$\begin{pmatrix}
x^{\prime}(t) \\ y^{\prime}(t) \\ z^{\prime}(t)
\end{pmatrix} = \begin{pmatrix}
x -y + z \\ 2x + y -z \\ x + y + z
\end{pmatrix} + \begin{pmatrix}
t -1 \\ -3t^{2} \\ t^{2} -t + 2
\end{pmatrix}$$

Finalmente podemos escribir

$$\begin{pmatrix}
x^{\prime}(t) \\ y^{\prime}(t) \\ z^{\prime}(t)
\end{pmatrix} = \begin{pmatrix}
1 & -1 & 1 \\ 2 & 1 & -1 \\ 1 & 1 & 1
\end{pmatrix} \begin{pmatrix}
x \\ y \\ z
\end{pmatrix} + \begin{pmatrix}
t -1 \\ -3t^{2} \\ t^{2} -t + 2
\end{pmatrix}$$

O bien,

$$\mathbf{Y}^{\prime} = \begin{pmatrix}
1 & -1 & 1 \\ 2 & 1 & -1 \\ 1 & 1 & 1
\end{pmatrix} \mathbf{Y} + \begin{pmatrix}
t -1 \\ -3t^{2} \\ t^{2} -t + 2
\end{pmatrix}$$

Donde,

$$\mathbf{A} = \begin{pmatrix} 1 & -1 & 1 \\ 2 & 1 & -1 \\ 1 & 1 & 1 \end{pmatrix} \hspace{1cm} y \hspace{1cm} \mathbf{G}(t) = \begin{pmatrix} t -1 \\ -3t^{2} \\ t^{2} -t + 2 \end{pmatrix}$$

$\square$

Usando la notación matricial verifiquemos que un vector solución en efecto es solución de un sistema lineal.

Ejemplo: Probar que el vector

$$\mathbf{Y} = \begin{pmatrix}
5 \cos(t) \\ 3 \cos(t) -\sin(t)
\end{pmatrix}e^{t}$$

es solución del sistema lineal

$$\begin{pmatrix}
x^{\prime}(t) \\ y^{\prime}(t)
\end{pmatrix} = \begin{pmatrix}
-2 & 5 \\ -2 & 4
\end{pmatrix} \begin{pmatrix}
x(t) \\ y(t)
\end{pmatrix}$$

Solución: El vector dado es

$$\mathbf{Y} = \begin{pmatrix}
x(t) \\ y(t)
\end{pmatrix} = \begin{pmatrix}
5e^{t} \cos(t) \\ 3e^{t} \cos(t) -e^{t} \sin(t)
\end{pmatrix}$$

Por una lado, derivemos el vector

$$\mathbf{Y}^{\prime} = \begin{pmatrix}
x^{\prime}(t) \\ y^{\prime}(t)
\end{pmatrix} = \begin{pmatrix}
5e^{t} \cos(t) -5e^{t} \sin(t) \\ 3e^{t} \cos(t) -3e^{t} \sin(t) -e^{t} \sin(t) -e^{t} \cos(t)
\end{pmatrix}$$

Esto es,

$$\mathbf{Y}^{\prime} = \begin{pmatrix}
5 \cos(t) -5 \sin(t) \\ 2 \cos(t) -4 \sin(t)
\end{pmatrix} e^{t}$$

Por otro lado, sustituyamos los valores de $x(t)$ y $y(t)$ en el sistema y veamos si se obtiene el mismo resultado.

$$\mathbf{Y}^{\prime} = \begin{pmatrix}
-2 & 5 \\ -2 & 4
\end{pmatrix} \begin{pmatrix}
5e^{t} \cos(t) \\ 3e^{t} \cos(t) -e^{t} \sin(t)
\end{pmatrix} = \begin{pmatrix}
-10e^{t} \cos(t) + 15e^{t} \cos(t) -5e^{t} \sin(t) \\ -10e^{t} \cos(t) + 12e^{t} \cos(t) -4e^{t} \sin(t)
\end{pmatrix}$$

Esto es,

$$\mathbf{Y}^{\prime} = \begin{pmatrix}
5 \cos(t) -5 \sin(t) \\ 2 \cos(t) -4 \sin(t)
\end{pmatrix} e^{t}$$

Como el resultado es el mismo concluimos que, en efecto, el vector $\mathbf{Y}$ es solución del sistema lineal dado.

$\square$

Para concluir con esta entrada veamos un resultado interesante que nos conecta con la unidad anterior.

¡Una ecuación diferencial de orden $n \geq 2$ lineal puede ser reescrita como un sistema lineal de $n$ ecuaciones de primer orden!.

Reducción de una ecuación de orden $n$ a un sistema de ecuaciones

Consideremos una ecuación diferencial lineal de orden $n$.

$$a_{n}(x) \dfrac{d^{n}y}{dx^{n}} + a_{n -1}(x) \dfrac{d^{n -1}y}{dx^{n -1}} + \cdots + a_{1}(x) \dfrac{dy}{dx} + a_{0}(x)y = g(x) \label{13} \tag{13}$$

Para adaptar este ejercicio a la notación que estamos usando en esta entrada tomemos a $x = x(t)$ como la variable dependiente de $t$ y dividamos toda la ecuación por $a_{n}(t) \neq 0$, tal que se obtenga la siguiente ecuación de orden $n$.

$$\dfrac{dx^{n}}{dt^{n}} + b_{1}(t) \dfrac{d^{n -1}x}{dt^{n -1}} + \cdots + b_{n -2}(t) \dfrac{d^{2}x}{dt^{2}} + b_{n -1}(t) \dfrac{dx}{dt} + b_{n}(t)x = g(t) \label{14} \tag{14}$$

Ahora realicemos las siguientes definiciones.

$$y_{1} = x, \hspace{1cm} y_{2} = \dfrac{dx}{dt}, \hspace{1cm} y_{3} = \dfrac{d^{2}x}{dt^{2}}, \hspace{1cm} \cdots, \hspace{1cm} y_{n} = \dfrac{d^{n -1}x}{dt^{n -1}} \label{15} \tag{15}$$

y notemos que

$$y^{\prime}_{1} = \dfrac{dx}{dt}, \hspace{1cm} y^{\prime}_{2} = \dfrac{d^{2}x}{dt^{2}}, \hspace{1cm} y^{\prime}_{3} = \dfrac{d^{3}x}{dt^{3}}, \hspace{1cm} \cdots, \hspace{1cm} y^{\prime}_{n -1} = \dfrac{d^{n -1}x}{dt^{n -1}} \label{16} \tag{16}$$

De los resultados (\ref{15}) y (\ref{16}) obtenemos que

$$y^{\prime}_{1} = y_{2}, \hspace{1cm} y^{\prime}_{2} = y_{3}, \hspace{1cm} y^{\prime}_{3} = y_{4}, \hspace{1cm} \cdots, \hspace{1cm} y^{\prime}_{n -1} = y_{n} \label{17} \tag{17}$$

Para obtener $y^{\prime}_{n}$ sólo despejamos de la ecuación diferencial (\ref{14}).

$$y^{\prime}_{n} = \dfrac{d^{n}x}{dt^{n}} = g(t) -b_{1}(t) \dfrac{d^{n -1}x}{dt^{n -1}} -\cdots -b_{n -2}(t) \dfrac{d^{2}x}{dt^{2}} -b_{n -1}(t) \dfrac{dx}{dt} -b_{n}(t)x$$

Si usamos (\ref{15}) podemos escribir

$$y^{\prime}_{n} = g(t) -b_{1}(t)y_{n} -\cdots -b_{n -2}(t)y_{3} -b_{n -1}(t)y_{2} -b_{n}(t)y_{1} \label{18} \tag{18}$$

Con estos resultados nos damos cuenta que hemos formado un sistema lineal de $n$ ecuaciones diferenciales.

\begin{align*}
y^{\prime}_{1} &= y_{2} \\
y^{\prime}_{2} &= y_{3} \\
y^{\prime}_{3} &= y_{4} \\
&\vdots \\
y^{\prime}_{n -1} &= y_{n} \\
y^{\prime}_{n} &= g(t) -b_{1}(t)y_{n} -\cdots -b_{n -2}(t)y_{3} -b_{n -1}(t)y_{2} -b_{n}(t)y_{1}
\end{align*}

Usando la notación matricial obtenemos finalmente que

$$\begin{pmatrix}
y^{\prime}_{1}(t) \\ y^{\prime}_{2}(t) \\ \vdots \\ y^{\prime}_{n -1}(t) \\ y^{\prime}_{n}(t)
\end{pmatrix} = \begin{pmatrix}
0 & 1 & 0 & \cdots & 0 \\ 0 & 0 & 1 & \cdots & 0 \\ \vdots & \vdots & \vdots & \ddots & \vdots \\ 0 & 0 & 0 & \cdots & 1 \\ -b_{n}(t) & -b_{n-1}(t) & -b_{n-2}(t) & \cdots & -b_{1}(t)
\end{pmatrix} \begin{pmatrix}
y_{1}(t) \\ y_{2}(t) \\ \vdots \\ y_{n -1}(t) \\ y_{n}(t)
\end{pmatrix} + \begin{pmatrix}
0 \\ 0 \\ \vdots \\ 0 \\ g(t)
\end{pmatrix}$$

Esto por supuesto trae muchas ventajas, ya que en ocasiones será mucho más sencillo resolver un sistema de $n$ ecuaciones con los métodos que veremos más adelante que intentar resolver la ecuación de orden $n$ con los métodos desarrollados en la unidad anterior.

Para que quede más claro el procedimiento anterior realicemos un ejemplo.

Ejemplo: Escribir la ecuación diferencial de orden $n = 4$

$$\dfrac{d^{4}x}{dt^{4}} + 12 \dfrac{d^{3}x}{dt^{3}} -5 \dfrac{d^{2}x}{dt^{2}} + 8x = 2 \cos(t)$$

en un sistema lineal usando notación matricial.

Solución: Aplicamos las definiciones de (\ref{15}) y (\ref{16}).

$$y_{1} = x, \hspace{1cm} y_{2} = \dfrac{dx}{dt} = y^{\prime}_{1}, \hspace{1cm} y_{3} = \dfrac{d^{2}x}{dt^{2}} = y^{\prime}_{2} \hspace{1cm} y \hspace{1cm} y_{4} = \dfrac{d^{3}x}{dt^{3}} = y^{\prime}_{3}$$

Y de la ecuación diferencial obtenemos que

$$\dfrac{d^{4}x}{dt^{4}} = 2 \cos(t) -12y_{4} + 5y_{3} -8y_{1} = y^{\prime}_{4}$$

El sistema que se forma, es

\begin{align*}
y^{\prime}_{1} &= y_{2} \\
y^{\prime}_{2} &= y_{3} \\
y^{\prime}_{3} &= y_{4} \\
y^{\prime}_{4} &= 2 \cos(t) -12y_{4} + 5y_{3} -8y_{1}
\end{align*}

Por lo tanto, la ecuación diferencial de orden $4$ es equivalente al sistema lineal de $4$ ecuaciones diferenciales

$$\begin{pmatrix}
y^{\prime}_{1}(t) \\ y^{\prime}_{2}(t) \\ y^{\prime}_{3}(t) \\ y^{\prime}_{4}(t)
\end{pmatrix} = \begin{pmatrix}
0 & 1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \\ -8 & 0 & 5 & -12
\end{pmatrix} \begin{pmatrix}
y_{1} \\ y_{2} \\ y_{3} \\ y_{4}
\end{pmatrix} + \begin{pmatrix}
0 \\ 0 \\ 0 \\ 2 \cos (t)
\end{pmatrix}$$

$\square$

Hemos concluido con esta entrada.

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. Escribir los siguientes sistemas lineales en forma matricial.
  • $\begin{align*}
    x^{\prime}(t) &= 3x -5y \\
    y^{\prime}(t) &= 4x + 8y
    \end{align*}$
  • $\begin{align*}
    x^{\prime}(t) &= -3x + 4y + e^{-t} \sin(2t) \\
    y^{\prime}(t) &= 5x + 9z + 4e^{-t} \cos(2t) \\
    z^{\prime}(t) &= y + 6z -e^{-t}
    \end{align*}$
  1. Reescribir los siguientes sistemas lineales sin el uso de matrices.
  • $\mathbf{Y}^{\prime} = \begin{pmatrix}
    7 & 5 & -9 \\ 4 & 1 & 1 \\ 0 & -2 & 3 \\
    \end{pmatrix} \mathbf{Y} + \begin{pmatrix}
    0 \\ 2 \\ 1
    \end{pmatrix} e^{5t} -\begin{pmatrix}
    8 \\ 0 \\ 3
    \end{pmatrix} e^{-2t}$
  • $\begin{pmatrix}
    x^{\prime}(t) \\ y^{\prime}(t) \\ z^{\prime}(t)
    \end{pmatrix} = \begin{pmatrix}
    1 & -1 & 2 \\ 3 & -4 & 1 \\ -2 & 5 & 6
    \end{pmatrix} \begin{pmatrix}
    x \\ y \\ z
    \end{pmatrix} + \begin{pmatrix}
    1 \\ 2 \\ 2
    \end{pmatrix} e^{-t} -\begin{pmatrix}
    3 \\ -1 \\ 1
    \end{pmatrix} t$
  1. Probar que el vector dado $\mathbf{Y}$ es solución del sistema lineal correspondiente.
  • $\mathbf{Y}^{\prime} = \begin{pmatrix}
    2 & 1 \\ -1 & 0
    \end{pmatrix} \mathbf{Y}, \hspace{1cm} \mathbf{Y} = \begin{pmatrix}
    1 \\ 3
    \end{pmatrix} e^{t} + \begin{pmatrix}
    4 \\ -4
    \end{pmatrix} te^{t}$
  • $\mathbf{Y}^{\prime} = \begin{pmatrix}
    1 & 0 & 1 \\ 1 & 1 & 0 \\ -2 & 0 & -1
    \end{pmatrix} \mathbf{Y}, \hspace{1cm} \mathbf{Y} = \begin{pmatrix}
    \sin(t) \\ -\dfrac{1}{2} \sin(t) -\dfrac{1}{2} \cos(t) \\ -\sin(t) + \cos(t)
    \end{pmatrix}$
  1. Escribir las siguientes ecuaciones diferenciales de orden superior en un sistema lineal usando notación matricial.
  • $\dfrac{d^{4}x}{dt^{4}} -10 \dfrac{d^{3}x}{dt^{3}} + 35 \dfrac{d^{2}x}{dt^{2}} -50 \dfrac{dx}{dt} + 24x = 0$
  • $\dfrac{d^{4}x}{dt^{4}} -4 \dfrac{d^{3}x}{dt^{3}} + 8 \dfrac{d^{2}x}{dt^{2}} -8 \dfrac{dx}{dt} + 4x = 8 \sin (2t)$

Más adelante…

Nos hemos introducido en los sistemas lineales de primer orden, en la siguiente entrada estudiaremos las propiedades de las soluciones de estos sistemas de manera muy similar que en el caso de las ecuaciones diferenciales de orden superior.

Veremos que mucho de lo visto en la unidad anterior aparecerá nuevamente, pues conceptos como dependencia e independencia lineal, conjunto fundamental de soluciones, Wronskiano, principio de superposición, entre otros, volverán a aparecer, sólo habrá que adaptarlos a los sistemas lineales.

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»