{"id":6285,"date":"2020-06-10T18:53:15","date_gmt":"2020-06-10T23:53:15","guid":{"rendered":"http:\/\/blog.espol.edu.ec\/analisisnumerico\/?p=6285"},"modified":"2026-04-05T19:45:18","modified_gmt":"2026-04-06T00:45:18","slug":"s1eva2009tii_t1-movimiento-particula-plano","status":"publish","type":"post","link":"https:\/\/blog.espol.edu.ec\/algoritmos101\/mn-s1eva10\/s1eva2009tii_t1-movimiento-particula-plano\/","title":{"rendered":"s1Eva2009TII_T1 Movimiento de part\u00edcula en plano"},"content":{"rendered":"\n<p><strong>Ejercicio<\/strong>: <a href=\"https:\/\/blog.espol.edu.ec\/algoritmos101\/mn-1eva10\/1eva2009tii_t1-movimiento-particula-plano\/\" data-type=\"post\" data-id=\"439\">1Eva2009TII_T1 Movimiento de part\u00edcula en plano<\/a><\/p>\n\n\n\n<figure class=\"wp-block-image alignright size-full is-resized\"><img loading=\"lazy\" decoding=\"async\" width=\"640\" height=\"480\" src=\"http:\/\/blog.espol.edu.ec\/algoritmos101\/files\/2020\/06\/s1EIIT20019T1MovPuntoPlano01.png\" alt=\"s1EIIT20019T1 Mov Punto Plano gr\u00e1fica 01\" class=\"wp-image-18714\" style=\"width:350px\" \/><\/figure>\n\n\n\n<h2 class=\"wp-block-heading\">a. Planteamiento del problema<\/h2>\n\n\n\n<p>Las ecuaciones expresan las trayectorias de dos part\u00edculas,<\/p>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> x(t) = 3 \\sin ^{3}(t)-1<\/span>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> y(t) = 4 \\sin (t)\\cos (t)<\/span>\n\n\n\n<p>que para que se encuentren o choquen, sus coordenadas deber\u00edan ser iguales.<\/p>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> x(t) = y(t) <\/span>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> 3 \\sin ^{3}(t)-1 = 4 \\sin (t)\\cos (t)<\/span>\n\n\n\n<p>Se reordena la expresi\u00f3n, de la forma f(t) = 0 para usarla en el algoritmo de b\u00fasqueda de ra\u00edces.<\/p>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> 3 \\sin ^{3}(t)-1 - 4 \\sin (t)\\cos (t) = 0<\/span>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> f(t) = 3 \\sin ^{3}(t)-1 - 4 \\sin (t)\\cos (t)<\/span>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\" \/>\n\n\n\n<h2 class=\"wp-block-heading\">b. Intervalo de b\u00fasqueda de ra\u00edz<\/h2>\n\n\n\n<p>Como la variable independiente es tiempo, el evento a buscar se supone sucede en tiempos positivos t&gt;=0, por lo que el valor inicial a la izquierda del intervalo ser\u00e1 a=0<\/p>\n\n\n\n<p>Para el caso de b, a la derecha, se usa lo indicado en el enunciado para la pregunta del literal b, donde se indica <strong>t<\/strong> \u2208 [0, \u03c0\/2], por lo que b = \u03c0\/2<\/p>\n\n\n\n<p class=\"has-text-align-center\">[0, \u03c0\/2]<\/p>\n\n\n\n<p>verificando que exista cambio de signo entre f(a) y f(b)<\/p>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> f(0) = 3 \\sin ^{3}(0)-1 - 4 \\sin (0)\\cos (0) = -1<\/span>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> f(\\pi \/2) = 3 \\sin ^{3}(\\pi \/2)-1 - 4 \\sin (\\pi \/2)\\cos (\\pi \/2)<\/span>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> = 3 (1)^3-1 - 4 (1)(0) = 2<\/span>\n\n\n\n<p>con lo que se comprueba que al existir cambio de signo [-,+], debe existir una ra\u00edz en el intervalo.<\/p>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\" \/>\n\n\n\n<h2 class=\"wp-block-heading\">c. M\u00e9todo de Falsa Posici\u00f3n<\/h2>\n\n\n\n<h3 class=\"wp-block-heading\">Desarrollo anal\u00edtico con l\u00e1piz y papel<\/h3>\n\n\n\n<p>Se trata de mostrar los pasos en al menos tres iteraciones del m\u00e9todo, usando las siguientes expresiones:<\/p>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> f(t) = 3 \\sin ^{3}(t)-1 - 4 \\sin (t)\\cos (t)<\/span>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> c = b - f(b) \\frac{a-b}{f(a)-f(b)} <\/span>\n\n\n\n<p class=\"has-text-align-center\">[0, \u03c0\/2]<\/p>\n\n\n\n<h4 class=\"wp-block-heading\">iteraci\u00f3n 1<\/h4>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> a = 0 , b = \\pi\/2 <\/span>\n\n\n\n<p>tomando los datos al validar los puntos extremos<\/p>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> f(0) = -1 <\/span>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> f(\\pi \/2) = 2 <\/span>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> c = \\pi\/2 - 2 \\frac{0-\\pi\/2}{-1-2} = \\pi\/6 <\/span>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> f(\\pi\/6) = 3 \\sin ^{3}(\\pi\/6)-1 - 4 \\sin (\\pi\/6)\\cos (\\pi\/6)<\/span>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> = -2.3570 <\/span>\n\n\n\n<p>el signo de f(c) es el mismo que f(a), se ajusta el lado izquierdo<\/p>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> tramo = |c-a| = |\\pi \/6 - 0| = \\pi\/6 <\/span>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> a = c = \\pi\/6 , b = \\pi\/2<\/span>\n\n\n\n<h4 class=\"wp-block-heading\">iteraci\u00f3n 2<\/h4>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> a = \\pi\/6 , b = \\pi\/2 <\/span>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> f(\\pi\/6) = -2.3570 <\/span>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> f(\\pi \/2) = 2 <\/span>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> c = \\pi\/2 - 2 \\frac{\\pi\/6-\\pi\/2}{-1-2} = 1.0901 <\/span>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> f(1.0901) = 3 \\sin ^{3}(1.0901)-1 - 4 \\sin (1.0901)\\cos (1.0901)<\/span>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> = -0.5486 <\/span>\n\n\n\n<p>el signo de f(c) es el mismo que f(a), se ajusta el lado izquierdo<\/p>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> tramo = |c-a| = | 1.0901-\\pi\/6| = 1.0722 <\/span>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> a = c = 1.0901 , b = \\pi\/2<\/span>\n\n\n\n<h4 class=\"wp-block-heading\">iteraci\u00f3n 3<\/h4>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> a = 1.0901 , b = \\pi\/2 <\/span>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> f(1.0901) = -0.5486 <\/span>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> f(\\pi \/2) = 2 <\/span>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> c = \\pi\/2 - 2 \\frac{-0.5486-\\pi\/2}{-0.5486-2} = 1.19358 <\/span>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> f(1.19358) = 3 \\sin ^{3}(1.19358)-1 - 4 \\sin (1.19358)\\cos (1.19358)<\/span>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> = 0.0409 <\/span>\n\n\n\n<p>el signo de f(c) es el mismo que f(b), se ajusta el lado derecho<\/p>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> tramo = |b-c| = | \\pi\/2- 1.19358| = 0.3772 <\/span>\n\n\n<span class=\"wp-katex-eq katex-display\" data-display=\"true\"> a = 1.0901 , b = 1.19358<\/span>\n\n\n\n<hr class=\"wp-block-separator has-alpha-channel-opacity\" \/>\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=\"M\u00e9todos de la Bisecci\u00f3n y Falsa Posici\u00f3n con Python - Ejercicio. Movimiento de part\u00edcula\" width=\"500\" height=\"281\" src=\"https:\/\/www.youtube.com\/embed\/8tIqNQrh9v0?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<hr class=\"wp-block-separator has-alpha-channel-opacity\" \/>\n\n\n\n<h2 class=\"wp-block-heading\">Algoritmo con Python<\/h2>\n\n\n\n<p>Los par\u00e1metros aplicados en el algoritmo son los desarrollados en el planteamiento del problema e intervalo de b\u00fasqueda, con lo que se obtiene los siguientes resultados:<\/p>\n\n\n\n<pre class=\"wp-block-code alignwide\"><code>m\u00e9todo de la Posici\u00f3n Falsa \ni &#091;'a', 'c', 'b'] &#091;'f(a)', 'f(c)', 'f(b)']\n   tramo\n0 &#091;0.       0.523599 1.570796] &#091;-1.       -2.357051  2.      ]\n    0.5235987755982989\n1 &#091;0.523599 1.090105 1.570796] &#091;-2.357051 -0.548671  2.      ]\n    0.5665065528830355\n2 &#091;1.090105 1.193587 1.570796] &#091;-0.548671  0.040977  2.      ]\n    0.37720921125460705\n3 &#091;1.090105 1.186396 1.193587] &#091;-0.548671 -0.000573  0.040977]\n    0.09629045311092965\n4 &#091;1.186396 1.186495 1.193587] &#091;-5.731531e-04 -2.335321e-07  4.097678e-02]\n    9.919955449055884e-05\nraiz:   1.1864949811467547<\/code><\/pre>\n\n\n\n<figure class=\"wp-block-image size-full\"><img loading=\"lazy\" decoding=\"async\" width=\"640\" height=\"480\" src=\"http:\/\/blog.espol.edu.ec\/algoritmos101\/files\/2020\/06\/s1EIIT20019T1MovPuntoPlano01.png\" alt=\"s1EIIT20019T1 Mov Punto Plano gr\u00e1fica 01\" class=\"wp-image-18714\" \/><\/figure>\n\n\n\n<p>las instrucciones en Python son:<\/p>\n\n\n<div class=\"wp-block-syntaxhighlighter-code alignwide\"><pre class=\"brush: python; title: ; notranslate\" title=\"\">\n# 1Eva_IIT2009_T1 Movimiento de part\u00edcula en plano\nimport numpy as np\n\n#INGRESO\nxt = lambda t: 3*(np.sin(t)**3)-1\nyt = lambda t: 4*np.sin(t)*np.cos(t)\nfx = lambda t: 3*(np.sin(t)**3)-1 - 4*np.sin(t)*np.cos(t) \n\na = 0\nb = np.pi\/2\ntolera = 0.001\n\n# intervalo para gr\u00e1fica\nLa = a\nLb = b\nmuestras = 21\n\n# PROCEDIMIENTO\ntabla = &#x5B;]\ntramo = abs(b-a)\nfa = fx(a)\nfb = fx(b)\nwhile not(tramo&lt;=tolera):\n    c = b - fb*(a-b)\/(fa-fb)\n    fc = fx(c)\n    unafila = &#x5B;a,c,b,fa,fc,fb,tramo]\n    cambio = np.sign(fa)*np.sign(fc)\n    if cambio&gt;0:\n        tramo = abs(c-a)\n        a = c\n        fa = fc\n    else:\n        tramo = abs(b-c)\n        b = c\n        fb = fc\n    unafila&#x5B;6]=tramo\n    tabla.append(unafila)\ntabla = np.array(tabla)\nntabla = len(tabla)\n \n# SALIDA\nnp.set_printoptions(precision=6)\nprint('m\u00e9todo de la Posici\u00f3n Falsa ')\nprint('i', &#x5B;'a','c','b'],&#x5B; 'f(a)', 'f(c)','f(b)'])\nprint('  ','tramo')\nfor i in range(0,ntabla,1):\n    print(i, tabla&#x5B;i,0:3],tabla&#x5B;i,3:6])\n    print('   ', tabla&#x5B;i,6])\n \nprint('raiz:  ',c)\n\n# GRAFICA -------------------\nimport matplotlib.pyplot as plt\nti = np.linspace(La,Lb,muestras)\nxi = xt(ti)\nyi = yt(ti)\nfi = fx(ti)\n\nplt.plot(ti,xi, label='x(t)')\nplt.plot(ti,yi, label='y(t)')\nplt.plot(ti,fi, label='f(t)')\nplt.plot(c,fx(c),'ro')\nplt.axhline(0, color='green')\nplt.axvline(c, color='magenta')\nplt.legend()\nplt.xlabel('t')\nplt.title('M\u00e9todo de Falsa Posici\u00f3n')\nplt.show()\n<\/pre><\/div>","protected":false},"excerpt":{"rendered":"<p>Ejercicio: 1Eva2009TII_T1 Movimiento de part\u00edcula en plano a. Planteamiento del problema Las ecuaciones expresan las trayectorias de dos part\u00edculas, que para que se encuentren o choquen, sus coordenadas deber\u00edan ser iguales. Se reordena la expresi\u00f3n, de la forma f(t) = 0 para usarla en el algoritmo de b\u00fasqueda de ra\u00edces. b. Intervalo de b\u00fasqueda de [&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":[44],"tags":[58,54],"class_list":["post-6285","post","type-post","status-publish","format-standard","hentry","category-mn-s1eva10","tag-ejemplos-python","tag-mnumericos"],"_links":{"self":[{"href":"https:\/\/blog.espol.edu.ec\/algoritmos101\/wp-json\/wp\/v2\/posts\/6285","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=6285"}],"version-history":[{"count":6,"href":"https:\/\/blog.espol.edu.ec\/algoritmos101\/wp-json\/wp\/v2\/posts\/6285\/revisions"}],"predecessor-version":[{"id":23793,"href":"https:\/\/blog.espol.edu.ec\/algoritmos101\/wp-json\/wp\/v2\/posts\/6285\/revisions\/23793"}],"wp:attachment":[{"href":"https:\/\/blog.espol.edu.ec\/algoritmos101\/wp-json\/wp\/v2\/media?parent=6285"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/blog.espol.edu.ec\/algoritmos101\/wp-json\/wp\/v2\/categories?post=6285"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/blog.espol.edu.ec\/algoritmos101\/wp-json\/wp\/v2\/tags?post=6285"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}