Scribd
Cargar
108 vistas4 páginas

Examen Pendulo Simple

Este documento presenta la simulación numérica de un péndulo simple utilizando dos métodos: el método de Euler (FPA) y el método de punto medio (LPA). El método de Euler se utiliza para calcular la posición y velocidad del péndulo en cada paso de tiempo. El método LPA es más preciso. Ambos métodos muestran que el período del péndulo es aproximadamente 0.33 segundos para el FPA y 0.31 segundos para el LPA.

Cargado por

Guillermo
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
108 vistas4 páginas

Examen Pendulo Simple

Este documento presenta la simulación numérica de un péndulo simple utilizando dos métodos: el método de Euler (FPA) y el método de punto medio (LPA). El método de Euler se utiliza para calcular la posición y velocidad del péndulo en cada paso de tiempo. El método LPA es más preciso. Ambos métodos muestran que el período del péndulo es aproximadamente 0.33 segundos para el FPA y 0.31 segundos para el LPA.

Cargado por

Guillermo
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd

4/5/2018 Examen pendulo simple

Péndulo simple
In [2]:

import numpy as np
import [Link] as plt
import math

In [3]:

x0 = [Link]/6 #Pasamos el ángulo a radianes.


w0 = 0
t0 = 0
l = 2.5
g = 10
N = 1000
dt = 0.01
Z = g/l
x = [Link](N+1)
x[0] = x0
w = [Link](N+1)
w[0] = w0
t = [Link](N+1)
t[0] = t0

FPA(Método de Euler)
In [4]:

for i in range (1,N+1):


t[i] = t[i-1] + dt
w[i] = w[i-1] + (-Z*[Link](x[i-1]))
x[i] = x[i-1] + w[i-1]*dt
print(w[i],x[i],t[i])

-2.0 0.523598775598 0.01


-4.0 0.503598775598 0.02
-5.93032259974 0.463598775598 0.03
-7.71900226792 0.404295549601 0.04
-9.29248707035 0.327105526922 0.05
-10.5777005932 0.234180656218 0.06
-11.5058849287 0.128403650286 0.07
-12.018089323 0.0133448009996 0.08
-12.0714669426 -0.10683609223 0.09
-11.6449350578 -0.227550761656 0.1
-10.7425666536 -0.344000112234 0.11
-9.39354449953 -0.45142577877 0.12
-7.64854877981 -0.545361223765 0.13
-5.57364098391 -0.621846711564 0.14
-3.2434923132 -0.677583121403 0.15
-0.735844751196 -0.710018044535 0.16
1.87154506961 -0.717376492047 0.17
4.5011852165 -0.698661041351 0.18
7.0739572885 -0.653649189186 0.19
9 5063070092 0 5829096163 0 2

[Link] 1/4
4/5/2018 Examen pendulo simple

In [19]:

T = 0.33 ##Del resultado obtenido deducimos que el periodo vale 0.33s


[Link](t,x)
[Link](0,2*T)
[Link](2.5,-2.5)
[Link]('Pendulo simple(FPA)')
[Link]('t/s')
[Link]('y/m')
[Link]()
##Como es normal en el método de Euler observamos como va perdiendo precisión en cada itera

LPA
In [6]:

t2 = [Link](N+1)
t2[0] = t0
w2 = [Link](N+1)
w2[0] = w0
x2 = [Link](N+1)
x2[0] = x0

[Link] 2/4
4/5/2018 Examen pendulo simple

In [7]:

for k in range (1,N+1):


t2[k] = t2[k-1] + dt
w2[k] = w2[k-1] + (-Z*[Link](x2[k-1]))
x2[k] = x2[k-1] + w2[k]*dt
print(w2[k],x2[k],t2[k])

-2.0 0.503598775598 0.01


-3.93032259974 0.464295549601 0.02
-5.72149474884 0.407080602113 0.03
-7.30521551902 0.334028446922 0.04
-8.61662139946 0.247862232928 0.05
-9.59794974521 0.151882735476 0.06
-10.2031475888 0.0498512595874 0.07
-10.4024700456 -0.0541734408685 0.08
-10.1858822573 -0.156032263442 0.09
-9.56428263716 -0.251675089813 0.1
-8.56817613026 -0.337356851116 0.11
-7.2441997735 -0.409798848851 0.12
-5.65050041344 -0.466303852985 0.13
-3.85214908598 -0.504825343845 0.14
-1.91753077263 -0.524000651571 0.15
0.0838612050362 -0.523162039521 0.16
2.0823481162 -0.502338558359 0.17
4.00825413146 -0.462256017045 0.18
5.79212807224 -0.404334736322 0.19
7 36575698338 0 330677166488 0 2
In [21]:

T2 = 0.31 ##En el caso del método LPA (que es mas preciso) deducimos de los valores obtenid
[Link](t2,x2)
[Link]()
[Link](0,2*T2)
[Link]('t/s')
[Link]('y/m')
[Link]('Pendulo simple(LPA)')

Out[21]:

<[Link] at 0x69adcc0>

[Link] 3/4
4/5/2018 Examen pendulo simple

[Link] 4/4

También podría gustarte