Práctica 1

Medina González Paola Ximena

Rayo Romero Jaqueline
Vargas Silva Daniela Alejandra
Villagomez Cuevas Héctor Alonso

Asignatura: Métodos Numéricos

Ingeniería en Biotecnología y Química

Periodo Febrero - Julio 2025

Fecha de entrega: 21 de Febrero del 2025

Docente: Lazcano González Roberto

 Problema 1
Es conocido que el periodo T de un péndulo simple de longitud viene dado por 𝑇 = 2π 𝐿/𝑔
siendo g la gravedad. Se pretende medir la aceleración de la gravedad g usando un
péndulo. Aproximar el error relativo cometido si las medidas efectuadas son L=1 ±0,001
metros y T =2± 0,01 segundos.

Procedimiento
Utilizando la fórmula del periodo del péndulo se despeja para obtener la fórmula de la
gravedad
𝑇 = 2π 𝐿/𝑔
2 2 2 2 2 2
𝑇 = (2π 𝐿/𝑔) → 𝑇 = 4π 𝐿/𝑔 → 𝑔 * 𝑇 = 4π * 𝐿
2 2
𝑔 = 4π * 𝐿/𝑇
Una vez obtenida la fórmula, se sustituyen los valores de L=1 y T=2 en la fórmula que da
2
como resultado =9.87 𝑚/𝑠 =g
Ahora, para obtener el valor relativo se utiliza la fórmula de propagación de errores relativos
proporcionada previamente

En la que como primer paso se tiene que obtener las derivadas parciales de cada uno.
La derivada de la gravedad con respecto a la longitud es dada por:
2 2 2
∂𝑔 4π 𝐿 4π ∂𝑔 ∂𝑔 4π
∂𝐿
= 2 → 2 ∂𝐿
=𝐿 → ∂𝐿
= 2
𝑇 𝑇 𝑇
La derivada de la gravedad con respecto al periodo es dada por:
2
∂𝑔 4π 𝐿 2 ∂𝑔 1 2 ∂𝑔 −2 2 ∂𝑔 −3 2 2
∂𝑇
= 2 → 4π 𝐿 ∂𝑇
= 2 → 4π 𝐿 ∂𝑇
=𝑇 → 4π 𝐿 ∂𝑇
=− 2𝑇 → 4π 𝐿 *− 3
𝑇 𝑇 𝑇
2
∂𝑔 8π
∂𝑇
= − 3
𝑇
Una vez teniendo las derivadas parciales solo faltaría calcular el valor relativo de cada
componente.La fórmula básica es
𝑣𝑎𝑙𝑜𝑟 𝑎𝑏𝑠.
ϵ= 𝑣𝑎𝑙𝑜𝑟 𝑟𝑒𝑎𝑙
0.001 0.01
ϵ𝐿 = 1
= 0. 001 ϵ𝑇 = 2
= 0. 005
Conociendo estos valores ya podemos utilizar la fórmula de propagación de errores
relativos, quedaría así:
𝐿 ∂𝑔 𝑇 ∂𝑔
ϵ𝑔 = || 𝑔 ∂𝐿 ||ϵ𝐿 + || 𝑔 ∂𝑇 ||ϵ𝑇
| 𝐿 4π2 | |𝑇 2
8π |
ϵ𝑔 = | 𝑔 2 |ϵ𝐿 + | 𝑔 − 3 |ϵ𝑇
| 𝑇 | | 𝑇 |
| 1 4π2 | | 2 2
8π |
ϵ𝑔 = | 9.87 2 |0. 001 + | 9.87 − 3 |0. 005
| 2 | | 2 |
ϵ𝑔 = |⋍0. 99|0. 001 + |⋍1. 99|0. 005
ϵ𝑔 = 0. 001 + 2 * 0. 005
ϵ𝑔 = 0. 011
El error relativo cometido fue de 0.011
 Problema 2
Lo resaltado en color rojo es uso y apoyo de IA
La velocidad límite v de una gota esférica de radio R compuesta de un fluido de viscosidad
µ 1 y densidad ρ 1 desplazándose bajo la acción de la gravedad constante g dentro de otro
fluido de viscosidad µ 2 ≠ 0 y densidad ρ 2 viene dada por:
2
2𝑅 (ρ1−ρ2)(µ1+µ2)
±𝑣 = 3µ2(2µ2+3µ1)

(i) Si se conocen de manera exacta los siguientes datos en el sistema internacional ρ1 = 10

, ρ2=15 y µ2 = 0, 3, pero sin embargo µ1 = 0, 2 ≠ 0, 1 y 𝑅 = 1, 5 ≠ 0, 2, aproximar el error
absoluto con el que obtiene la velocidad v.

(ii) Suponiendo que solo se tienen errores en la medida de R, averiguar si el problema de

evaluar v está mal condicionado para alguna relación entre las variables del problema.

Procedimiento
El problema nos solicita observar como los errores afectan la velocidad. Lo que se dedujo
fue el utilizar las derivadas parciales, básicamente estas nos indican cuánto cambia la
velocidad (v) cuando cambian R y µ 1

​ Fórmula general para la incertidumbre en una función v(𝑥1, 𝑥2...) es:

𝑛
| ∂𝑣 |
∆𝑣 ≈ ∑ | ∂𝑥 |∆𝑥𝑖
𝑖=1| 𝑖 |
en términos de interpretación de las diferenciales para cambios pequeños en las
variables (𝑑𝑥𝑖 ≈ ∆𝑥𝑖) obtenemos la “versión” usada en el problema
∂𝑣 | ∂𝑣 |
∆𝑣 ≈ || ∂𝑅 ||∆𝑅 +| ∂µ |∆µ1
| 1|
Datos:
∂𝑣
●​ || ∂𝑅 || se refiere a la derivada parcial de “v” con respecto a la R
| ∂𝑣 |
●​ | ∂µ |es la derivada parcial de la “v” con respecto a µ 1
| 1|
●​ ∆𝑅 y ∆µ1 son los errores en la edición de R y µ 1
​ Valores dados:
​ R= 1.5 con error en ∆𝑅= ± 0. 2
​ µ1= 0.2 con error en ∆µ1 =± 0. 2
​ ρ1 = 10
​ ρ2 = 0. 3
​ µ2 = 0. 3
​ 𝑔 = 9. 81

∂𝑣 | ∂𝑣 |
Lo que se entendió es que se deben calcular primero || ∂𝑅 || y | ∂µ |
| 1|
Derivada parcial respecto a la R
2
2𝑅 (ρ1−ρ2)(µ1+µ2)
𝑣= 3µ2(2µ2+3µ1)

∂𝑣
Sacamos || ∂𝑅 || (la derivada de R) de la ecuación:
𝑑𝑣 2
𝑑𝑅
= (2𝑅 𝑔(ρ1 − ρ2)(µ1 + µ2)) = 4𝑅𝑔(ρ1 − ρ2)(µ1 + µ2)
El denominador es constante respecto a R y no tiene R así que se queda igual
4𝑅𝑔(ρ1−ρ2)(µ1+µ2)
| ∂𝑣 | =
| ∂𝑅 | 3µ2(2µ2+3µ1)

Sustituyendo valores queda:

| ∂𝑣 | = 4(1.5)(9,81)(10−0.3)(0.2+0.3) = 264. 33
| ∂𝑅 | 3(0.3)(2(0.3)+3(0.2))
| ∂𝑣 | ≈ 264. 33
| ∂𝑅 |

Derivada respecto a µ1
| ∂𝑣 |
Sacamos | ∂µ | de la ecuación:
| 1|
Se busco en IA solo para verificar la derivada ya que nos revolvimos (se espera que se
pueda entender este desglose de procedimiento):
Derivando el numerador y denominador separados queda:
Numerador
2
2𝑅 𝑔(ρ1 − ρ2)(µ1 + µ2)
𝑑𝑣 2
µ1
= 2𝑅 𝑔(ρ1 − ρ2)
Denominador
3µ2(2µ2 + 3µ1)
𝑑𝑣
µ1
(3µ2(2µ2 + 3µ1) = 9µ2
(Aquí ya no se sabía bien cómo incorporar ambas derivadas, se buscó que lo más
conveniente era aplicar la regla del cociente)
2 2
𝑑𝑣 (2𝑅 𝑔(ρ1−ρ2))(3µ2(2µ2+3µ1))−(2𝑅 𝑔(ρ1−ρ2)(µ1+µ2))(9µ2)
µ1
= 2
3µ2(2µ2+3µ1)

Factorización y simplificación
2
| ∂𝑣 | 6𝑅 𝑔(ρ1−ρ2)µ2(µ2+3µ1)
| ∂µ | = − 2
| 1| 3µ2(2µ2+3µ1)

Sustituyendo

| ∂𝑣 | 2
6(1.5) (9.81)(9.7)(0.3)(0.3+3(0.2))
| ∂µ | = − 2 = − 98. 2
| 1| (3(0.3)(2(0.3)+3(0.2))
| ∂𝑣 |
| ∂µ | ≈− 98. 2
| 1|

Se calcula el error ya final:

 ∂𝑣 | ∂𝑣 |
Haciendo una resta con la fórmula ∆𝑣 ≈ || ∂𝑅 ||∆𝑅 +| ∂µ |∆µ1
| 1|

∆𝑣 ≈ |264. 33|(0. 2) + |− 98. 2|(0. 2)

∆𝑣 ≈ 264. 33(0. 2) + (98. 2)(0. 2) = 52. 86 + 19. 64 = 72. 5
∆𝑣 ≈ 72. 5

Finalmente la incertidumbre o el error en la velocidad es de 72.5 a consecuencia de los

errores en la medición de R y µ1

Conclusiones
Fue un análisis en dónde se analizaron los errores en las mediciones, utilizando técnicas de
errores y como calcular derivadas parciales, con el propósito de evaluar la incertidumbre en
la determinación de la aceleración de la gravedad con la ayuda de un péndulo y en la
velocidad límite de una gota esférica en un fluido. Esto destacando la propagación de
errores en las mediciones experimentales y reforzando la necesidad de minimizar la
incertidumbre en variables con mayor alteración en los cálculos para así mejorar la
precisión y disminuir el rango de error por medio de un buen método numérico.