import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
import sys
sys.path.insert(1, '../Base')
from reg_lin import reg_lin_w as rl
import mpl_config as mpl
mpl.inicio(4)#, [7., 5.])
D = [pd.read_csv("CC-{}.csv".format(i), sep=';', decimal=',') for i in range(1,4)]
I = np.asarray([[D[j]["I{}".format(i)].to_numpy() for i in range(1,5)] for j in range(0,3)])
T = np.asarray([[D[j]["T{}".format(i)].to_numpy() for i in range(1,5)] for j in range(0,3)])
LnI = np.asarray([[D[j]["LnI{}".format(i)].to_numpy() for i in range(1,5)] for j in range(0,3)])
sLnI = np.asarray([[D[j]["sLnI{}".format(i)].to_numpy() for i in range(1,5)] for j in range(0,3)])
#Selector de experimento
exp = 2 #Número de experiencia (0-2)
rep = 3 #Número de repetición (0-3)
c = ["royalblue", "mediumseagreen", "sandybrown", "tomato"]
for r in range(4):
x = T[exp,r][~np.isnan(T[exp,r])]
y = LnI[exp,r][~np.isnan(LnI[exp,r])]
sy = sLnI[exp,r][~np.isnan(sLnI[exp,r])]
a, b = rl(x,y,sy)[0:2]
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
import sys
sys.path.insert(1, '../Base')
import mpl_config as mpl
mpl.inicio(3)
#Constantes
pm = 4 * np.pi * 10**(-7)
i = 2.55
n = 154
r = 0.2
#Calculo campo magnético
B = lambda pm, i, n, r, a, z: ((pm * i * n) / (2 * r)) * ((1 / (1 + ((z - (a/2)) / (r))**2)**(3/2)) + (1 / (1 + ((z + (a/2)) / (r))**2)**(3/2)))
#Leer datos
d1 = pd.read_csv("BH-1.csv", sep=';', decimal=',')
z1 = d1["z"]
Be1 = d1["Bexp"]
d2 = pd.read_csv("BH-2.csv", sep=';', decimal=',')
z2 = d2["z"]
Be2 = d2["Bexp"]
d3 = pd.read_csv("BH-3.csv", sep=';', decimal=',')
z3 = d3["z"]
Be3 = d3["Bexp"]
#Curvas teóricas
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
import uncertainties as unc
import sys
sys.path.insert(1, '../Base')
from reg_lin import reg_lin_b as rl
from varias_medidas import tratamiento_datos as tdatos
import mpl_config as mpl
mpl.inicio(3, [3, 2.2])
#Datos
m1 = 0.20089
m2 = 0.26079
m3 = 0.30167
L = 0.1 #Diafragma
sm = 0.00001
st = 0.001
sl = 0.001
d = [
pd.read_csv("LN_MRU_M{}.csv".format(i), sep=';', decimal=',')
for i in range(1, 4)
]
S = np.array([d[j]["S"] for j in range(0, 3)])
T = np.array([d[j]["T"] for j in range(0, 3)])
T1 = np.array([d[j]["T1"] for j in range(0, 3)])
T2 = np.array([d[j]["T2"] for j in range(0, 3)])
import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
import uncertainties as unc
import uncertainties.unumpy as unp
import sys
sys.path.insert(1, '../Base')
from varias_medidas import tratamiento_datos as tdatos
from reg_lin import reg_lin as rl
import mpl_config as mpl
mpl.inicio(1)
#CONSTANTES GLOBALES
D = 0.02167
sD = 0.00048
sb = 0.001 * 2
#FUNCIONES
def m_inercia_exp(d):
#Tratamiento Datos
T = np.zeros(len(d))
sT = np.zeros(len(d))
for i in range(len(d)):
T[i], sT[i] = tdatos(d[i], sb)
print("---\nT{} = {:.2u}".format(i * 3, unc.ufloat(T[i], sT[i])))
#Momento Inercia