Métodos Numéricos Ingenierı́a Quı́mica/Mecánica CURSO 2021/2022 Laboratorio

Práctica 4: Diferenciación Numérica

Objetivo: aplicar las fórmulas de dos y tres puntos para aproximar las derivadas de funciones que
se pueden plantear en problemas de ingenierı́a.

Si los datos que tenemos son experimentales, es decir, vienen dados en una tabla, simplemente
habrá que aplicar las fórmulas adecuadas con el tamaño de paso que nos proporciona la tabla:

Ejemplo: La distancia s(t) recorrida por un móvil se muestra en la siguiente tabla

t (s) 8 9 10
s (m) 17.453 21.460 25.752
(a) Aproxima la velocidad v(9) mediante las fórmulas de la diferencia progresiva, regresiva y
diferencias centradas.
(b) Aproxima la aceleración en t = 9 s.
Podemos introducir los datos de la tabla o bien en una matriz 3×2, o más simple en dos vectores,
uno con los datos de la variable independiente, t en este ejemplo, y otro con los valores de la
variable dependiente, s. El cálculo de las aproximaciones de la velocidad y aceleración en t = 9
se refleja en el siguiente fichero de comandos

1 % % Ejemplo de Derivación numérica

2 % A partir de una tabla de valores calcular velocidad y aceleración en t=9
3 clc;clear
4 t=8:10; % o bien t=[8,9,10], serı́a prescindible
5 h=1; % tama~
no de paso
6 s=[17.453,21.460, 25.752];
7 % cálculo velocidad: derivación aproximada en t=9
8 vprog=(s(3)-s(2))/h % diferencia progresiva
9 vreg=(s(2)-s(1))/h % diferencia regresiva
10 vcen=(s(3)-s(1))/(2*h) % diferencias centradas
11 % cálculo de aceleración: derivada segunda aproximada en t=9
12 acelaprox=(s(3)-2*s(2)+s(1))/(h^2) % derivada segunda aprox.

Observa que basta con la definición del tamaño de paso, h = 1, para aplicar las fórmulas de
derivación y que el vector correspondiente a la variable independiente, t, no serı́a necesario
definirlo. Si los valores de la variable independiente no fueran equiespaciados sı́ serı́a necesaria
su definición.

En el caso de tener que aproximar las derivadas de una función dada veamos como podemos
actuar

19
Métodos Numéricos Ingenierı́a Quı́mica/Mecánica CURSO 2021/2022 Laboratorio

Ejemplo: La distancia s(t) recorrida por un móvil viene dada por

s(t) = −70 + 7t + 70e−t/10

(a) Aproxima la velocidad v(9) mediante la fórmulas de las diferencias centradas con tamaños
de paso h = 10−k k = 1, 2, . . . , 10
(b) Evalúa los errores cometidos.

En este caso será preciso definir la función, s(t), y su derivada, s0 (t) para comparar con los
valores obtenidos por la fórmula de diferenciación. El fichero de instrucciones serı́a

1 % % Ejemplo derivación numérica

2 % a partir de la función dada
3 clc,clear
4 format long
5 % definición de la función espacio, s(t),y su derivada
6 s6=@(t) -70+7*t+70*exp(-t/10);
7 der1s6=@(t) 7-7*exp(-t/10);
8

9 % % Cálculo aproximado de s’(9) con distintos valores de h y error cometido

10 % diferencia centrada
11 kmax=10;
12 k=1:kmax;hk=10.^(-k);
13 dercen=(s6(9+hk)-s6(9-hk))./(2*hk); % diferencias centradas
14 erdercen=abs(dercen-der1s6(9)); % error cometido
15 disp(’h aprox. error’)
16 [hk’ dercen’ erdercen’]

Nota: MatLab trabaja con variables matriciales, por tanto la programación de las diferencias
centradas, en este ejemplo, puede hacerse para el rango completo de valores de tamaño de paso
con una única sentencia.

20
Métodos Numéricos Ingenierı́a Quı́mica/Mecánica CURSO 2021/2022 Laboratorio

Práctica 4

Ejercicio 1 Sea s(t) la posición (metros) de un paracaidı́sta, a partir del momento en que se abre
su paracaı́das.
381 −16t/25 381(16t − 5)
s(t) = e + t≥0
80 400
(a) Aproxima la velocidad a los t = 15 segundos, mediante las fórmulas de las diferencias
progresiva y regresiva con tamaños de paso h = 10−k k = 1, 2, . . . , 10.
(b) Evalúa los errores cometidos comparando con el valor real. ¿A partir de que h cambia la
tendencia decreciente y los errores comienzan a aumentar?
(c) Acota el error de truncamiento de la aproximación de la derivada mediante las fórmulas de
las diferencias progresiva y regresiva para los tamaños de paso h = 10−k k = 1, 2, . . . , 10
(utiliza el intervalo [14.9, 15.1]).
(d) Halla el tamaño de paso, H1, que minimiza el error de truncamiento de la fórmula de
diferencias centradas. Calcula la aproximación y error cometido para el tamaño de paso
hallado, H1.
(e) Aproxima la aceleración a los t = 15 segundos, mediante la fórmula de la derivada segunda
con tamaños de paso h = 10−k k = 1, 2, . . . , 10 y calcula los errores cometidos.
(f) Acota el error de truncamiento de la aproximación de la derivada segunda para los tamaños
de paso h = 10−k k = 1, 2, . . . , 10 (utiliza el intervalo [14.9, 15.1]).
(g) Halla el tamaño de paso, H2, que minimiza el error de truncamiento de la fórmula de
la derivada segunda. Calcula la aproximación y error cometido para el tamaño de paso
hallado, H2.

Ejercicio 2 Para un circuito con voltaje E(t) e inductancia L = 0.06 henrios, la primera ley de
Kirchoff expresa la relación
dI
E=L +RI
dt
siendo R = 4 Ohmios la resistencia del circuito e I la intensidad de corriente. Se obtuvieron
las siguientes medidas de la intensidad I en varios instantes de tiempo

t (s) 1 1.03 1.06 1.09 1.12 1.15 1.18 1.21 1.24 1.27 1.3
I (A) 6.41 6.22 6.7 5.73 5.54 5.25 5.0 4.56 4.27 3.92 3.59

(a) Aproxima el voltaje E(t) en los instantes de tiempo de la tabla utilizando las fórmulas de
tres puntos adecuadas en cada instante.
(b) Aproxima el voltaje en t = 1.054 s utilizando interpolación cuadrática sobre los resultados
obtenidos en el apartado anterior.

7
Métodos Numéricos Ingenierı́a Quı́mica/Mecánica CURSO 2021/2022 Laboratorio

Ejercicio 3- La ley de Newton para el enfriamiento dice que la tasa de enfriamiento es proporcional
a la diferencia de temperaturas del cuerpo y del medio. Se expresa como:
dT
= k(T − Ta )
dt
donde T es la temperatura del cuerpo (◦ C), Ta es la temperatura del medio (◦ C) y k es una
constante de proporcionalidad. Una bola de cobre se calienta hasta una temperatura de 100o C
en el instante t = 0, se sumerge en agua que se mantiene a Ta = 30◦ C constante, la temperatura
de la bola cambia, ası́

tiempo, (min) 0 2 4 6 8 10 12 14 16 18 20
◦
T ( C) 100 78 62 52.8 45.7 40.8 37.5 35.1 33.5 32.4 31.7

dT
(a) Utiliza la diferenciación numérica de orden O(h2 ) para determinar en cada valor del
dt
tiempo.
dT
(b) Dibuja frente a T − Ta para los datos del intervalo de tiempo [0, 20], y estima el valor
dt
de k.

Nota: para estimar el valor de k realiza un ajuste lineal a los puntos dibujados:

Primero abre la ventana gráfica de Matlab con los puntos dibujados: Open in figure window
En la lı́nea del menú elige: Tools/Basic Fitting
En la ventana de ajuste elige: linear/Show equation,5 Significant Digits.
La pendiente de la recta será la aproximación de k buscada.