#ms.ple("torque", M, dM) #ms.pve("deflection force", Dt, dDt) #ms.ps("linreg vs graph", 0.0217, 0.0004, Dt, dDt) # 1.2 Deflecting force determination by period time measurement Jd = md / 2.0 * rd**2 dJd = rd * ms.sqrt((0.5 * rd * dmd)**2 + (md * drd)**2) Dp = 4.0 * ms.pi**2 * Jd / (T2**2 - T1**2) dDp = 4.0 * ms.pi**2 / (T2**2 - T1**2) * ms.sqrt(dJd**2 + ((2.0 * Jd * T1 * dT1)**2 + (2.0 * Jd * T2 * dT2)**2) / (T2**2 - T1**2)**2) print() ms.pve("T1", T1, dT1, False) ms.pve("T2", T2, dT2) ms.pve("Jd", Jd, dJd) ms.pve("Dp", Dp, dDp) # Deviation between 1.1 and 1.2 ms.ps("Deviation", Dt, Dp, dDt, dDp, False) D = ms.mean_value([Dt, Dp]) dD = 1 / 2 * ms.sqrt(dDt**2 + dDp**2) ms.pve("Mean Value", D, dD) # 2 Verification of the theorem of Steiner a2 = [] da2 = [] J = [] dJ = []
return deg + mins / 60.0 lamHg = [690.7, 623.4, 579.1, 577.0, 546.1, 499.2, 491.6, 435.8, 407.8, 404.7] phiHg = [mtd(181, 56), mtd(180, 34), mtd(180, 18), mtd(180, 16), mtd(179, 59), mtd(179, 30), mtd(179, 26), mtd(178, 31), mtd(177, 55), mtd(177, 52)] dphiHg = [mtd(0, 5)] * len(phiHg) phiHe = [mtd(180, 49), mtd(180, 22), mtd(179, 34), mtd(179, 28), mtd(179, 9), mtd(178, 48)] dphiHe = [mtd(0, 5)] * len(phiHe) phiH = [mtd(180, 45), mtd(179, 21), mtd(178, 29), mtd(178, 0)] dphiH = [mtd(0, 5)] * len(phiH) for i in range(len(phiHg)): ms.pv("lambda[" + str(i) + "]", lamHg[i]) print() for i in range(len(phiHg)): ms.pve("deltaHg(lambda)[" + str(i) + "]", phiHg[i], dphiHg[i]) print() for i in range(len(phiHe)): ms.pve("deltaHe(lambda)[" + str(i) + "]", phiHe[i], dphiHe[i]) print() for i in range(len(phiH)): ms.pve("deltaH(lambda)[" + str(i) + "]", phiH[i], dphiH[i]) for i in range(len(phiHg)): phiHg[i] *= m.pi / 180.0 for i in range(len(phiHe)): phiHe[i] *= m.pi / 180.0 for i in range(len(phiH)): phiH[i] *= m.pi / 180.0
import math as m import measure as ms # Gasthermometer p0 = 1013.25 pl = 1020.2 dpl = 0.2 p = [916, 945, 991, 1027, 1061, 1092, 1125, 1163, 198, 1226, 1245] dp = [1, 5, 6, 5, 5, 5, 3, 2, 5, 4, 2] pc = 653 dpc = 2 pn = 255 dpn = 1 pnb = p[len(p) - 1] * p0 / pl dpnb = p0 / pl * m.sqrt(dp[len(p) - 1]**2 + p[len(p) - 1]**2 / pl**2 * dpl**2) ms.pve("pnb", pnb, dpnb)
mw = 497.61e-3 dmw = 0.13e-3 mw -= mk dmw = m.sqrt(dmw**2 + dmk**2) Tw = 55.8 dTw = 0.2 Tl = 24.1 dTl = 0.1 Tm = 53.25 # 53.19 dTm = 0.27 # 0.4 W = mw * cw * (Tw - Tm) / (Tm - Tl) dW = m.sqrt((Tw - Tm)**2 * ((cw * dmw)**2 + (mw * dcw)**2) + (mw * cw)**2 * (dTw**2 + (((Tl - Tw) * dTm)**2 + ((Tw - Tm) * dTl)**2) / (Tm - Tl)**2)) / (Tm - Tl) ms.pve("mw", mw, dmw) ms.pve("Wasserwert", W, dW) # T=100°C R = 8.3144598 Mp = [207.2e-3, 26.9815385e-3, 12.011e-3] p = 1024.6 dp = 2.0 mp = [501.773e-3, 131.08e-3, 125.05e-3] dmp = [0.06e-3, 0.04e-3, 0.02e-3] mw = [612.38e-3, 600.03e-3, 615.13e-3] dmw = [0.04e-3, 0.09e-3, 0.05e-3] Tw = [27.4, 24.7, 23.4] dTw = [0.2, 0.2, 0.2] Tp = [100 + 0.0276 * (p - p0) for i in range(len(mp))]
#print(f) #print(d_f) #print(sigma) # Auflösungsvermögen Mikroskop s = 41.5 ds = 1.0 b = 250 f = 40 lam = 550e-6 beta = b / f - 1.0 a = [0.2, 0.5, 0.4] da = ms.std_dev_m(a) a = ms.mean_value(a) alpha = m.atan(a / (2 * s)) dalpha = 1.0 / m.sqrt(4 * s**2 + a**2) * m.sqrt(da**2 + (a / s)**2 * ds**2) g = 1 / (2 * beta) gMin = 0.61 * lam / m.sin(alpha) dgMin = abs(0.61 * lam * (m.cos(alpha) / (m.sin(alpha)**2)) * dalpha) ms.pve("a", a, da) ms.pv("beta", beta) ms.pve("alpha", alpha * 180 / m.pi, dalpha * 180 / m.pi) ms.pv("g", g) ms.pve("gMin", gMin, dgMin) ms.pv("sigma", abs(gMin - g) / dgMin)