Universidad Nacional de
Ingeniería
Fundación Extensus
Maestría en Ingeniería Estructural y Sísmica.
Métodos Numéricos
Tarea # 05
Docente: Dra. María Penkova
Julio, 2017
Universidad Nacional de
Ingeniería
Fundación Extensus
1. El objetivo de este problema es comparar las aproximaciones por diferencias finitas de
segundo orden exactas hacia delante, atrás y centradas, de la primera derivada de una
función con el valor real de la derivada. Esto se hará para
a) Use el cálculo para determinar el valor correcto de la derivada en x=2
b) Para evaluar las aproximaciones por diferencias finitas centradas, comience con x
= 0.5, asi para la primera evaluación, los valores de x para la aproximación por
diferencias centradas será x= 2 ±0.5 o x= 1.5 y 2.5. entonces disminuya en pasos de
0.01 hacia abajo hasta un valor mínimo de ∆x=0.01
CENTRADA
x+ x-
i x ΔX ΔX ΔX f(x) f(x+ΔX) f(x-ΔX) f'(x) error
1 2 0.5 2.5 1.5 52.59815003 145.913159 18.5855369 127.327622 17.682048
2 2 0.49 2.49 1.51 52.59815003 142.984382 18.9812917 126.533765 16.948329
3 2 0.48 2.48 1.52 52.59815003 140.113796 19.3852432 125.758909 16.232171
4 2 0.47 2.47 1.53 52.59815003 137.30025 19.7975572 125.002864 15.533400
5 2 0.46 2.46 1.54 52.59815003 134.542613 20.2184024 124.265447 14.851845
6 2 0.45 2.45 1.55 52.59815003 131.83978 20.6479513 123.546476 14.187339
7 2 0.44 2.44 1.56 52.59815003 129.190664 21.0863796 122.845778 13.539721
8 2 0.43 2.43 1.57 52.59815003 126.594202 21.5338669 122.163181 12.908834
9 2 0.42 2.42 1.58 52.59815003 124.049352 21.9905959 121.498519 12.294523
10 2 0.41 2.41 1.59 52.59815003 121.555091 22.4567536 120.851631 11.696639
11 2 0.4 2.4 1.6 52.59815003 119.110418 22.9325302 120.222359 11.115037
12 2 0.39 2.39 1.61 52.59815003 116.71435 23.4181202 119.610551 10.549576
13 2 0.38 2.38 1.62 52.59815003 114.365926 23.9137217 119.016058 10.000118
14 2 0.37 2.37 1.63 52.59815003 112.064202 24.4195371 118.438736 9.466531
15 2 0.36 2.36 1.64 52.59815003 109.808253 24.9357727 117.878444 8.948683
16 2 0.35 2.35 1.65 52.59815003 107.597172 25.4626389 117.335048 8.446451
Universidad Nacional de
Ingeniería
Fundación Extensus
17 2 0.34 2.34 1.66 52.59815003 105.430073 26.0003506 116.808415 7.959713
18 2 0.33 2.33 1.67 52.59815003 103.306082 26.5491267 116.298417 7.488350
19 2 0.32 2.32 1.68 52.59815003 101.224348 27.1091909 115.804932 7.032248
20 2 0.31 2.31 1.69 52.59815003 99.1840321 27.6807711 115.327840 6.591298
21 2 0.3 2.3 1.7 52.59815003 97.1843156 28.2641 114.867026 6.165392
22 2 0.29 2.29 1.71 52.59815003 95.2243942 28.859415 114.422378 5.754428
23 2 0.28 2.28 1.72 52.59815003 93.3034798 29.4669582 113.993789 5.358306
24 2 0.27 2.27 1.73 52.59815003 91.4208001 30.0869765 113.581155 4.976931
25 2 0.26 2.26 1.74 52.59815003 89.575598 30.7197221 113.184377 4.610210
26 2 0.25 2.25 1.75 52.59815003 87.7671313 31.365452 112.803359 4.258056
27 2 0.24 2.24 1.76 52.59815003 85.9946727 32.0244285 112.438009 3.920382
28 2 0.23 2.23 1.77 52.59815003 84.2575091 32.6969192 112.088239 3.597109
29 2 0.22 2.22 1.78 52.59815003 82.5549417 33.3831971 111.753965 3.288158
30 2 0.21 2.21 1.79 52.59815003 80.8862854 34.0835408 111.435106 2.993453
31 2 0.2 2.2 1.8 52.59815003 79.2508687 34.7982344 111.131586 2.712926
32 2 0.19 2.19 1.81 52.59815003 77.6480334 35.5275678 110.843330 2.446507
33 2 0.18 2.18 1.82 52.59815003 76.0771344 36.2718367 110.570271 2.194134
34 2 0.17 2.17 1.83 52.59815003 74.5375393 37.0313429 110.312343 1.955744
35 2 0.16 2.16 1.84 52.59815003 73.0286283 37.8063941 110.069482 1.731281
36 2 0.15 2.15 1.85 52.59815003 71.5497937 38.5973044 109.841631 1.520691
37 2 0.14 2.14 1.86 52.59815003 70.10044 39.4043941 109.628735 1.323923
38 2 0.13 2.13 1.87 52.59815003 68.6799835 40.2279902 109.430743 1.140929
39 2 0.12 2.12 1.88 52.59815003 67.2878518 41.068426 109.247608 0.971667
40 2 0.11 2.11 1.89 52.59815003 65.9234843 41.9260417 109.079284 0.816095
41 2 0.1 2.1 1.9 52.59815003 64.586331 42.8011845 108.925733 0.674175
42 2 0.09 2.09 1.91 52.59815003 63.2758532 43.6942083 108.786916 0.545874
43 2 0.08 2.08 1.92 52.59815003 61.9915226 44.6054744 108.662801 0.431162
44 2 0.07 2.07 1.93 52.59815003 60.7328214 45.5353514 108.553358 0.330009
45 2 0.06 2.06 1.94 52.59815003 59.4992423 46.4842151 108.458560 0.242393
46 2 0.05 2.05 1.95 52.59815003 58.2902876 47.4524491 108.378385 0.168291
47 2 0.04 2.04 1.96 52.59815003 57.1054698 48.4404448 108.312813 0.107687
48 2 0.03 2.03 1.97 52.59815003 55.9443111 49.4486013 108.261830 0.060565
49 2 0.02 2.02 1.98 52.59815003 54.8063428 50.4773259 108.225421 0.026915
50 2 0.01 2.01 1.99 52.59815003 53.6911058 51.5270342 108.203580 0.006728
c) Repita el inciso b) para las diferencias de segundo orden hacia adelante y hacia
atrás. (Observe que esto se puede hacer al mismo tiempo que la diferencia
centrada se calcula en el lazo)
Universidad Nacional de
Ingeniería
Fundación Extensus
Adelante Atrás
x+ x-
i x ΔX ΔX ΔX f(x) f(x+ΔX) f(x-ΔX) f'(x) error f'(x) error
1 2 0.5 2.5 1.5 52.598 145.913 18.586 186.630018 72.492052 68.025226 37.127955
2 2 0.49 2.49 1.51 52.598 142.984 18.981 184.461697 70.487990 68.605833 36.591331
3 2 0.48 2.48 1.52 52.598 140.114 19.385 182.324262 68.512474 69.193556 36.048131
4 2 0.47 2.47 1.53 52.598 137.300 19.798 180.217233 66.565061 69.788495 35.498261
5 2 0.46 2.46 1.54 52.598 134.543 20.218 178.140137 64.645313 70.390756 34.941624
6 2 0.45 2.45 1.55 52.598 131.840 20.648 176.092510 62.752802 71.000442 34.378124
7 2 0.44 2.44 1.56 52.598 129.191 21.086 174.073895 60.887105 71.617660 33.807663
8 2 0.43 2.43 1.57 52.598 126.594 21.534 172.083842 59.047807 72.242519 33.230139
9 2 0.42 2.42 1.58 52.598 124.049 21.991 170.121909 57.234498 72.875129 32.645452
10 2 0.41 2.41 1.59 52.598 121.555 22.457 168.187660 55.446776 73.515601 32.053498
11 2 0.4 2.4 1.6 52.598 119.110 22.933 166.280669 53.684247 74.164050 31.454172
12 2 0.39 2.39 1.61 52.598 116.714 23.418 164.400513 51.946520 74.820589 30.847368
13 2 0.38 2.38 1.62 52.598 114.366 23.914 162.546779 50.233214 75.485338 30.232977
14 2 0.37 2.37 1.63 52.598 112.064 24.420 160.719059 48.543951 76.158413 29.610889
15 2 0.36 2.36 1.64 52.598 109.808 24.936 158.916952 46.878361 76.839937 28.980994
16 2 0.35 2.35 1.65 52.598 107.597 25.463 157.140064 45.236079 77.530032 28.343177
17 2 0.34 2.34 1.66 52.598 105.430 26.000 155.388007 43.616748 78.228822 27.697322
18 2 0.33 2.33 1.67 52.598 103.306 26.549 153.660400 42.020014 78.936434 27.043315
19 2 0.32 2.32 1.68 52.598 101.224 27.109 151.956867 40.445530 79.652997 26.381034
20 2 0.31 2.31 1.69 52.598 99.184 27.681 150.277039 38.892956 80.378642 25.710360
21 2 0.3 2.3 1.7 52.598 97.184 28.264 148.620552 37.361954 81.113500 25.031170
22 2 0.29 2.29 1.71 52.598 95.224 28.859 146.987049 35.852195 81.857707 24.343340
23 2 0.28 2.28 1.72 52.598 93.303 29.467 145.376178 34.363354 82.611400 23.646743
24 2 0.27 2.27 1.73 52.598 91.421 30.087 143.787593 32.895111 83.374717 22.941250
25 2 0.26 2.26 1.74 52.598 89.576 30.720 142.220954 31.447151 84.147800 22.226731
26 2 0.25 2.25 1.75 52.598 87.767 31.365 140.675925 30.019164 84.930792 21.503053
27 2 0.24 2.24 1.76 52.598 85.995 32.024 139.152178 28.610847 85.723840 20.770082
28 2 0.23 2.23 1.77 52.598 84.258 32.697 137.649387 27.221899 86.527091 20.027681
29 2 0.22 2.22 1.78 52.598 82.555 33.383 136.167235 25.852025 87.340695 19.275710
30 2 0.21 2.21 1.79 52.598 80.886 34.084 134.705406 24.500936 88.164806 18.514029
31 2 0.2 2.2 1.8 52.598 79.251 34.798 133.263593 23.168346 88.999578 17.742494
32 2 0.19 2.19 1.81 52.598 77.648 35.528 131.841491 21.853974 89.845170 16.960959
33 2 0.18 2.18 1.82 52.598 76.077 36.272 130.438802 20.557544 90.701741 16.169277
34 2 0.17 2.17 1.83 52.598 74.538 37.031 129.055231 19.278784 91.569454 15.367296
35 2 0.16 2.16 1.84 52.598 73.029 37.806 127.690489 18.017427 92.448475 14.554865
36 2 0.15 2.15 1.85 52.598 71.550 38.597 126.344291 16.773209 93.338971 13.731827
37 2 0.14 2.14 1.86 52.598 70.100 39.404 125.016357 15.545871 94.241114 12.898025
38 2 0.13 2.13 1.87 52.598 68.680 40.228 123.706411 14.335158 95.155076 12.053300
Universidad Nacional de
Ingeniería
Fundación Extensus
39 2 0.12 2.12 1.88 52.598 67.288 41.068 122.414182 13.140821 96.081034 11.197487
40 2 0.11 2.11 1.89 52.598 65.923 41.926 121.139402 11.962611 97.019166 10.330421
41 2 0.1 2.1 1.9 52.598 64.586 42.801 119.881810 10.800286 97.969655 9.451936
42 2 0.09 2.09 1.91 52.598 63.276 43.694 118.641146 9.653608 98.932686 8.561859
43 2 0.08 2.08 1.92 52.598 61.992 44.605 117.417157 8.522340 99.908445 7.660017
44 2 0.07 2.07 1.93 52.598 60.733 45.535 116.209592 7.406253 100.897124 6.746235
45 2 0.06 2.06 1.94 52.598 59.499 46.484 115.018204 6.305117 101.898916 5.820332
46 2 0.05 2.05 1.95 52.598 58.290 47.452 113.842751 5.218710 102.914019 4.882128
47 2 0.04 2.04 1.96 52.598 57.105 48.440 112.682995 4.146810 103.942631 3.931436
48 2 0.03 2.03 1.97 52.598 55.944 49.449 111.538702 3.089201 104.984958 2.968070
49 2 0.02 2.02 1.98 52.598 54.806 50.477 110.409639 2.045669 106.041204 1.991839
50 2 0.01 2.01 1.99 52.598 53.691 51.527 109.295579 1.016005 107.111581 1.002548
d) Grafique los resultados de b y c versus x. Para efectos de comparación, incluya el
resultado exacto de la gráfica.
Grafico
200.000000
180.000000
160.000000
140.000000
120.000000
Diferencia Centrada
100.000000
Hacia Adelante
80.000000 Hacia Atras
60.000000
40.000000
20.000000
0.000000
0 0.1 0.2 0.3 0.4 0.5 0.6
Universidad Nacional de
Ingeniería
Fundación Extensus
2. Use los datos siguientes para encontrar la velocidad y aceleración en t=10 segundos:
Emplee los métodos de diferencias finitas correctas de segundo orden
a) Centradas
Velocidad
Aceleración
b) Hacia adelante
Velocidad
Aceleración
Universidad Nacional de
Ingeniería
Fundación Extensus
c) Hacia atrás
Velocidad
Aceleración
3. Evalúe la integral siguiente:
a) En forma analítica
b) Con una sola aplicación de la regla del trapecio
c) Con aplicación múltiple de la regla del trapecio, con n=2 y 4
Universidad Nacional de
Ingeniería
Fundación Extensus
Universidad Nacional de
Ingeniería
Fundación Extensus
d) Con una sola aplicación de la regla de Simpson 1/3
e) Con una sola aplicación de la regla de Simpson 1/3 con n=4
Universidad Nacional de
Ingeniería
Fundación Extensus
f) Con una sola aplicación de la regla de Simpson 3/8
g) Con la aplicación múltiple de la regla de Simpson 1/3 con n = 5.
4. Evalúe la integral de los datos tabulados a continuación, con
Universidad Nacional de
Ingeniería
Fundación Extensus
a) La regla del trapecio
b) Las reglas de Simpson
5. Determine la distancia recorrida para los datos siguientes:
a) Use la regla del trapecio
b) La mejor combinación de las reglas del trapecio y de Simpson
Universidad Nacional de
Ingeniería
Fundación Extensus
c) La integración analítica de polinomios de segundo y tercer orden, determinados
por regresión
Obteniendo la gráfica por regresión polinómica de 2do orden
Grafica de aceleracion
y = -0.1117x2 + 1.3526x + 3.4447
9
8
7
6
5
4
3
2
1
0
0 2 4 6 8 10 12
Regresión polinómica de 3er orden
Universidad Nacional de
Ingeniería
Fundación Extensus
Grafica de aceleracion
y = -0.0165x3 + 0.1625x2 + 0.0788x + 4.8547
9
8
7
6
5
4
3
2
1
0
0 2 4 6 8 10 12
6. Las áreas (A) de la sección transversal de una corriente se requieren para varias tareas
de la ingeniería de recursos hidráulicos, como el pronóstico del escurrimiento y el
diseño de presas. A menos que se disponga de dispositivos electrónicos muy avanzados
para obtener perfiles continuos del fondo del canal, el ingeniero debe basarse en
mediciones discretas de la profundidad para calcular A. El área de la sección transversal
de un canal se calcula con:
𝐵𝐵
𝐴𝐴𝑐𝑐 = � 𝐻𝐻(𝑦𝑦)𝑑𝑑𝑑𝑑
0
donde B = ancho total del canal (m), H = profundidad (m), y y = distancia desde uno
de los márgenes (m).
En forma similar, el flujo promedio Q (m3/s) se calcula por medio de:
𝐵𝐵
𝑄𝑄 = � 𝑈𝑈𝑈𝑈(𝑦𝑦)𝐻𝐻(𝑦𝑦)𝑑𝑑𝑑𝑑
0
donde U = velocidad del agua (m/s). Use estas relaciones y algún método numérico
para determinar Ac y Q, para los datos siguientes:
Universidad Nacional de
Ingeniería
Fundación Extensus
Utilizando la regla del trapecio:
Universidad Nacional de
Ingeniería
Fundación Extensus
7. Una viga de 11m está sujeta a una carga, y la fuerza cortante sigue la ecuación:
Donde V es la fuerza cortante y x es la distancia a lo largo de la viga. Se sabe que V=dM/dx,
y M es el momento flexionante. La integración conduce a la relación:
Si M0 es cero y x=11, calcule M con el empleo de:
Universidad Nacional de
Ingeniería
Fundación Extensus
d) Integración analítica
e) Aplicación múltiple de la regla del trapecio
Universidad Nacional de
Ingeniería
Fundación Extensus
El valor usando la regla múltiple del trapecio es 277.75
f) Aplicación múltiple de las reglas de Simpson.
8. Sea la integral
a) Emplee la integración de Romberg para evaluar con una cantidad de es = 0.5%. Debe
presentar sus resultados en la forma de una tabla como se mostró en las lecciones.
La solución analítica
Para el 1er intervalo con n=1b
Para el intervalo con n = 2
Universidad Nacional de
Ingeniería
Fundación Extensus
Para el intervalo con n = 4
Para el intervalo con n = 8
• Usando el método de Romberg
Para n=1,2
Para n = 2,4
Se muestran los resultados en la tabla
1er 2da 3ra 4ta
1.34377
1.81556 1.97283
1.91172 1.94377 1.94184
1.93310 1.94023 1.93999 1.93996
b) Obtenga una estimación de la integral con el uso de la fórmula de Gauss-Legendre con
cinco puntos.
Universidad Nacional de
Ingeniería
Fundación Extensus
Realizando un cambio de variable antes de la integración
Sustituyendo en la ecuación inicial tenemos:
Entonces:
Usando la fórmula de gauss-Legendre para 5 puntos:
Donde:
Universidad Nacional de
Ingeniería
Fundación Extensus
9. La cantidad de masa transportada por un tubo durante cierto periodo de tiempo se
calcula con
donde M = masa (mg), t1 = tiempo inicial (min), t2 = tiempo final (min), Q(t) = tasa de
flujo (m3/min), y c(t) = concentración (mg/m3). Las representaciones funcionales
siguientes de definen las variaciones temporales en el flujo y la concentración:
Determine la masa transportada entre t1 = 2 min y t2 = 8 min, con integración de
Romberg para una tolerancia de 0.1%
Universidad Nacional de
Ingeniería
Fundación Extensus
Universidad Nacional de
Ingeniería
Fundación Extensus
