{"id":2358,"date":"2018-08-29T06:55:07","date_gmt":"2018-08-29T11:55:07","guid":{"rendered":"http:\/\/blog.espol.edu.ec\/matg1013\/?p=2358"},"modified":"2026-04-05T20:14:43","modified_gmt":"2026-04-06T01:14:43","slug":"s2eva2018ti_t4-dragado-acceso-maritimo","status":"publish","type":"post","link":"https:\/\/blog.espol.edu.ec\/algoritmos101\/mn-s2eva20\/s2eva2018ti_t4-dragado-acceso-maritimo\/","title":{"rendered":"s2Eva2018TI_T4 Dragado acceso mar\u00edtimo"},"content":{"rendered":"\n<p><em><strong>Ejercicio<\/strong><\/em>: <a href=\"https:\/\/blog.espol.edu.ec\/algoritmos101\/mn-2eva20\/2eva2018ti_t4-dragado-acceso-maritimo\/\" data-type=\"post\" data-id=\"2347\">2Eva2018TI_T4 Dragado acceso mar\u00edtimo<\/a><\/p>\n\n\n\n<p>a) La matriz para remover sedimentos se determina como la diferencia entre la profundidad y la matriz de batimetr\u00eda.<\/p>\n\n\n\n<p>Considere el signo de la profundidad para obtener el resultado:<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code>matriz remover sedimentos: \n&#091;&#091; 4.21  0.    0.96  0.    3.76  3.09]\n &#091; 2.15  0.11  2.05  3.77  0.    3.07]\n &#091; 0.    1.14  1.65  0.    1.62  1.35]\n &#091; 3.7   0.    0.59  2.33  0.    4.23]\n &#091; 0.    1.38  3.53  4.49  1.98  1.4 ]\n &#091; 0.    0.77  0.32  1.06  4.24  3.54]]<\/code><\/pre>\n\n\n\n<figure class=\"wp-block-embed is-type-video is-provider-youtube wp-block-embed-youtube wp-embed-aspect-16-9 wp-has-aspect-ratio\"><div class=\"wp-block-embed__wrapper\">\n<iframe loading=\"lazy\" title=\"Integral num\u00e9rico doble con Python. Ejercicio dragado acceso mar\u00edtimo\" width=\"500\" height=\"281\" src=\"https:\/\/www.youtube.com\/embed\/1kkzGVK4C-0?feature=oembed\" frameborder=\"0\" allow=\"accelerometer; autoplay; clipboard-write; encrypted-media; gyroscope; picture-in-picture; web-share\" referrerpolicy=\"strict-origin-when-cross-origin\" allowfullscreen><\/iframe>\n<\/div><\/figure>\n\n\n\n<p>se obtiene con la instrucciones:<\/p>\n\n\n<div class=\"wp-block-syntaxhighlighter-code \"><pre class=\"brush: python; title: ; notranslate\" title=\"\">\n# 2da Evaluaci\u00f3n I T\u00e9rmino 2018\n# Tema 4. canal acceso a Puertos de Guayaquil\nimport numpy as np\n\n# INGRESO\nprofcanal = 11\n\nxi = np.array(&#x5B; 0.,  50., 100., 150., 200., 250.])\nyi = np.array(&#x5B; 0., 100., 200., 300., 400., 500.])\n\nbatimetria = &#x5B;&#x5B; -6.79,-12.03,-10.04,-11.60, -7.24,-7.91],\n              &#x5B; -8.85,-10.89, -8.95, -7.23,-11.42,-7.93],\n              &#x5B;-11.90, -9.86, -9.35,-12.05, -9.38,-9.65],\n              &#x5B; -7.30,-11.55,-10.41, -8.67,-11.84,-6.77],\n              &#x5B;-12.17, -9.62, -7.47, -6.51, -9.02,-9.60],\n              &#x5B;-11.90,-10.23,-10.68, -9.94, -6.76,-7.46]]\n\nbatimetria = np.array(batimetria)\n# PROCEDIMIENTO\n&#x5B;n,m] = np.shape(batimetria)\n\n# Matriz remover sedimentos\nremover = batimetria + profcanal\nfor i in range(0,n,1):\n    for j in range(0,m,1):\n        if remover&#x5B;i,j]&lt;0:\n            remover&#x5B;i,j]=0\n# SALIDA\nprint('matriz remover sedimentos: ')\nprint(remover)\n<\/pre><\/div>\n\n\n<p>b) el volumen se calcula usando el algoritmo de Simpson primero por un eje, y luego con el resultado se contin\u00faa con el otro eje,<\/p>\n\n\n\n<p>Considere que existen 6 puntos, o 5 tramos integrar en cada eje.<\/p>\n\n\n\n<ul class=\"wp-block-list\">\n<li>Al usar Simpson de 1\/3 que usan tramos pares, faltar\u00eda integrar el \u00faltimo tramo.<\/li>\n\n\n\n<li>En el caso de Simpson de 3\/8 se requieren tramos m\u00faltiplos de 3, por lo que faltar\u00eda un tramo para volver a usar la f\u00f3rmula.<\/li>\n<\/ul>\n\n\n\n<p>La soluci\u00f3n por filas se implementa usando una combinaci\u00f3n de <strong>Simpson 3\/8<\/strong> para los puntos entre remover[i, 0:(2+1)] y <strong>Simpson 1\/3<\/strong> para los puntos entre remover[i, 3:(4+1)].<\/p>\n\n\n\n<p>Luego se completa el integral del otro eje con el resultado anterior, aplicando el mismo m\u00e9todo.<\/p>\n\n\n\n<p>Se obtuvieron los siguientes resultados:<\/p>\n\n\n\n<pre class=\"wp-block-code\"><code>Integral en eje x: \n&#091;435.1  346.5  287.44 255.58 590.54 440.52]\nVolumen:  199161.11111111115<\/code><\/pre>\n\n\n\n<p>que se obtiene usando las instrucciones a continuaci\u00f3n de las anteriores:<\/p>\n\n\n<div class=\"wp-block-syntaxhighlighter-code \"><pre class=\"brush: python; title: ; notranslate\" title=\"\">\n# literal b) ---------------------------\ndef simpson13(fi,h):\n    s13 = (h\/3)*(fi&#x5B;0] + 4*fi&#x5B;1] + fi&#x5B;2])\n    return(s13)\ndef simpson38(fi,h):\n    s38 = (3*h\/8)*(fi&#x5B;0] + 3*fi&#x5B;1] + 3*fi&#x5B;2]+ fi&#x5B;3])\n    return(s38)\n\nIntegralx = np.zeros(n,dtype = float)\n\nfor i in range(0,n,1):\n    hx = xi&#x5B;1]-xi&#x5B;0]\n    fi = remover&#x5B;i, 0:(0+4)]\n    s38 = simpson38(fi,hx)\n    fi = remover&#x5B;i, 3:(3+3)]\n    s13 = simpson13(fi,hx)\n    Integralx&#x5B;i] = s38 + s13\n\nhy = yi&#x5B;1] - yi&#x5B;0]\nfj = Integralx&#x5B;0:(0+4)]\ns38 = simpson38(fj,hy)\nfj = Integralx&#x5B;3:(3+3)]\ns13 = simpson13(fj,hy)\nvolumen = s38 + s13\n\n# Salida\nnp.set_printoptions(precision=2)\nprint('Integral en eje x: ')\nprint(Integralx)\nprint('Volumen: ', volumen)\n<\/pre><\/div>\n\n\n<p>Para el examen escrito, se requieren realizar al menos 3 iteraciones\/ filas para encontrar el integral.<\/p>\n","protected":false},"excerpt":{"rendered":"<p>Ejercicio: 2Eva2018TI_T4 Dragado acceso mar\u00edtimo a) La matriz para remover sedimentos se determina como la diferencia entre la profundidad y la matriz de batimetr\u00eda. Considere el signo de la profundidad para obtener el resultado: se obtiene con la instrucciones: b) el volumen se calcula usando el algoritmo de Simpson primero por un eje, y luego [&hellip;]<\/p>\n","protected":false},"author":8043,"featured_media":0,"comment_status":"closed","ping_status":"open","sticky":false,"template":"wp-custom-template-entrada-mn-ejemplo","format":"standard","meta":{"footnotes":""},"categories":[48],"tags":[58,54],"class_list":["post-2358","post","type-post","status-publish","format-standard","hentry","category-mn-s2eva20","tag-ejemplos-python","tag-mnumericos"],"_links":{"self":[{"href":"https:\/\/blog.espol.edu.ec\/algoritmos101\/wp-json\/wp\/v2\/posts\/2358","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/blog.espol.edu.ec\/algoritmos101\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/blog.espol.edu.ec\/algoritmos101\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/blog.espol.edu.ec\/algoritmos101\/wp-json\/wp\/v2\/users\/8043"}],"replies":[{"embeddable":true,"href":"https:\/\/blog.espol.edu.ec\/algoritmos101\/wp-json\/wp\/v2\/comments?post=2358"}],"version-history":[{"count":4,"href":"https:\/\/blog.espol.edu.ec\/algoritmos101\/wp-json\/wp\/v2\/posts\/2358\/revisions"}],"predecessor-version":[{"id":23864,"href":"https:\/\/blog.espol.edu.ec\/algoritmos101\/wp-json\/wp\/v2\/posts\/2358\/revisions\/23864"}],"wp:attachment":[{"href":"https:\/\/blog.espol.edu.ec\/algoritmos101\/wp-json\/wp\/v2\/media?parent=2358"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/blog.espol.edu.ec\/algoritmos101\/wp-json\/wp\/v2\/categories?post=2358"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/blog.espol.edu.ec\/algoritmos101\/wp-json\/wp\/v2\/tags?post=2358"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}