Ejemplo de integración numérica

De Laplace

1 Enunciado

Una partícula se mueve a lo largo de una recta, siendo su velocidad (en el SI) como función del tiempo, la dada por la gráfica

La partícula parte de $s = 0$ .

Aprovechando los puntos en que la curva cruza la cuadrícula, calcule aproximadamente la posición en que se encontrará la partícula en $t=16\,\mathrm{s}$ .
Calcule el valor exacto de esta posición, sabiendo que la ley para la velocidad es

$v = \sqrt{1+16t-t^2}$

¿Cuál es el error relativo cometido en el apartado anterior?

2 Solución numérica

Si conocemos la velocidad instantánea, podemos hallar el desplazamiento a base de sumar los desplazamientos diferenciales

$\Delta s = \int_0^T v(t)\,\mathrm{d}t$

Gráficamente, esto equivale a hallar el área bajo la curva de la velocidad frente al tiempo.

Examinando la figura, podemos ver que la curva corta la cuadrícula en los puntos de coordenadas enteras (en el SI) (0,1), (1,4), (4,7), (7,8), (9,8), (12,7), (15,4) y (16,1):

El área bajo la curva puede entonces ser aproximada por una suma de trapecios, siendo la fórmula para el área de cada uno

$\Delta S = \frac{b+B}{2}h = \frac{v(t_1)+v(t_2)}{2}\Delta t$

Aplicando esto a nuestro caso obtenemos el desplazamiento aproximado

$\Delta s = \left(\frac{1+4}{2}\cdot 1+\frac{4+7}{2}\cdot 3+\frac{7+8}{2}\cdot 3+\frac{8+8}{2}\cdot 2 + \frac{8+7}{2}\cdot 3+\frac{7+4}{2}\cdot 3+\frac{4+1}{2}\cdot 1\right)\,\mathrm{m} = 99\,\mathrm{m}$

Existen otras formas numéricas de aproximar el resultado, según como se aproxime la curva. En este caso, se puede hacer coincidir con muy buena aproximación por media circunferencia de radio $\sqrt{65}$ , con lo que el área valdría

$\Delta s \simeq \frac{65 \pi}{2}\,\mathrm{m}=102.10\,\mathrm{m}$

3 Solución analítica

Integrando analíticamente la función obtenemos, midiendo todas las magnitudes en el SI

$\Delta s =\int_0^{16} \sqrt{1+16t-t^2}\,\mathrm{d}t$

Para resolver esta integral, primero completamos cuadrados, sumando y restando 64

$\Delta s =\int_0^{16} \sqrt{65-(t-8)^2}\,\mathrm{d}t$

Haciendo un cambio de variable

$u = t -8\qquad\Rightarrow\qquad \Delta s = \int_{-8}^8 \sqrt{65-u^2}\,\mathrm{d}u$

Con un nuevo cambio de variable

$u = \sqrt{65}\,\mathrm{sen}(\varphi)\qquad\Rightarrow\qquad \Delta s = 65 \int_{-\varphi_0}^{\varphi_0} \cos^2(\varphi)\,\mathrm{d}\varphi$

siendo

$\varphi_0=\mathrm{arcsen}\left(\frac{8}{\sqrt{65}}\right)=\,\mathrm{arctg}(8)$

La integral trigonométrica se calcula de forma sencilla

$\Delta s = 65 \int_{-\varphi_0}^{\varphi_0} \cos^2(\varphi)\,\mathrm{d}\varphi = \frac{65}{2}\int_{-\varphi_0}^{\varphi_0}(1+\cos(2\varphi))\mathrm{d}\varphi=65\left(\varphi_0+\frac{1}{2}\mathrm{sen}(2\varphi_0)\right)$

Sustituyendo el valor de $\varphi_0$