Homec4science
Diffusion yearII_reports (master) Browse

rSTICAZZI/labo/C1/errors.py
384c0a276062master

Last Change
Raw File

errors.py
View Options

import math
Tp = 955 + 273.0
Tm = 19 + 273.0
dTp = 5.0
dTm = 0.05
Vp = 184.45
Vm = 10.57
dVp = 21.36
dVm = 2.47
c = 2.5
a = [15, 17, 18, 20, 23, 25]
da = [3, 2, 1, 0.5, 0.5, 0.5]
F = [1.2, 0.55, 0.21, 0.16, 0.1, 0.06]
dF = [0.1, 0.05, 0.01, 0.01, 0.01, 0.01]
w = [252.6, 318.4, 373.2, 387.95, 396.6, 406.85]
dw = [1.7, 4.7, 0.1, 0.05, 2.3, 1.85]
R = 0.04
dR = 0.0005
U = 14.53
dU = 0.005
I = 14.04
dI = 0.005
# adjustments
for i in range(len(F)):
a[i] *= math.pi / 180.0
da[i] *= math.pi / 180.0 # rad
w[i] /= 60.0 # 1/s
dw[i] /= 60.0
def flux_2(U, I):
return U * I
def errFlux_2(U, I, dU, dI):
return I * dU + U * dI
def rend_emp(W, U, I):
return W / flux_2(U, I)
def errRend_emp(W, U, I, dW, dU, dI):
return dW / flux_2(U, I) + W * errFlux_2(U, I, dU, dI) / flux_2(U, I)**2
def errV(h, d, dh, dd):
return math.pi * ((d/2)**2 * dh + math.pi/2 * h * d * dd)
def errSigma(Tp, Tm, dTp, dTm):
return dTm / Tp + Tm / Tp**2 * dTp
def errPhi(Vp, Vm, dVp, dVm):
return dVp / Vp + dVm / Vm
def sigma(Tp, Tm):
return 1 - (Tm / Tp)
def phi(Vp, Vm):
return math.log(Vp/Vm)
def rend(Tp, Tm, Vp, Vm, c):
return 1.0 / (1 / sigma(Tp, Tm) + c / phi(Vp, Vm))
def errRend(Tp, Tm, dTp, dTm, Vp, Vm, dVp, dVm, c):
return rend(Tp, Tm, Vp, Vm, c)**2 * \
(errSigma(Tp, Tm, dTp, dTm) / sigma(Tp, Tm)**2 + c * errPhi(Vp, Vm, dVp, dVm) / phi(Vp, Vm)**2)
def puiss(F, R, w, a):
return F * R * w * math.cos(a)
def errPuiss(F, R, w, a, dF, dR, dw, da):
cos = math.cos(a)
return R * cos * w * dF + \
F * cos * w * dR + \
F * R * cos * dw + \
F * R * w * math.sin(a) * da
eta = rend(Tp, Tm, Vp, Vm, c)
d_eta = errRend(Tp, Tm, dTp, dTm, Vp, Vm, dVp, dVm, c)
print "eta_th = ", eta, " += ", d_eta
P = []
dP = []
Eta = []
dEta = []
for i in range(len(F)):
P.append(puiss(F[i], R, w[i], a[i]))
dP.append(errPuiss(F[i], R, w[i], a[i], dF[i], dR, dw[i], da[i]))
Eta.append(rend_emp(P[i], U, I))
dEta.append(errRend_emp(P[i], U, I, dP[i], dU, dI))
print "puissance =", P[i], " +- ", dP[i]
print "Eta =", Eta[i], " +- ", dEta[i]
out_table = open("mesure.dat", 'w')
out_table.write("%s %s %s %s \n" % ("w", "dw", "Eta", "dEta"))
for i in range(len(P)):
out_table.write("%f %f %f %f \n" % (w[i], dw[i], Eta[i], dEta[i]))
out_table.close()
c4science · Help