# -*- coding: utf-8 -*-
import matplotlib.pyplot as plt
import numpy as np
import sympy as sy
import ocd
messung251=ocd.open_csv("messung251.csv")
messung252=ocd.open_csv("messung252.csv")
messung253=ocd.open_csv("messung253.csv")
messung254=ocd.open_csv("messung254.csv")
U=sy.Symbol("U")
I=sy.Symbol("I")
R=sy.Symbol("R")
var=[U,I,R]
def plot253():
plt.figure()
(a,b,Sa,Sb,Sy)=ocd.plot_var(I,U,var,messung232,True)
plt.xlabel("I in A")
plt.ylabel("U in V")
plt.title(ur"Spannungsstabilisiertes Netzgerät")
print "a:",a,"+-",Sa
print "b:",b,"+-",Sb
plt.xlim(0,1.4)
plt.show()
def plot252():
plt.figure()
(a,b,Sa,Sb,Sy)=ocd.plot_var(I,U,var,messung252,True)
plt.xlabel("I in A")
plt.ylabel("U in V")
plt.title(ur"Bleiakku")
# -*- coding: utf-8 -*-
import matplotlib.pyplot as plt
import numpy as np
import sympy as sy
import ocd
messung21=ocd.open_csv("messung21.csv")
messung22=ocd.open_csv("messung22.csv")
messung231=ocd.open_csv("messung231.csv")
messung2321=ocd.open_csv("messung2321.csv")
messung2322=ocd.open_csv("messung2322.csv")
U=sy.Symbol("U")
I=sy.Symbol("I")
var=[U,I]
def plot21():
plt.figure(figsize=(16,12))
(a,b,Sa,Sb,Sy)=ocd.plot_var(I,U,var,messung21,True)
ocd.plot_var(I,U,var,messung22)
plt.xlabel("I in A")
plt.ylabel("U in V")
plt.title(ur"Ohmscher Widerstand und Glühbirne")
print "a:",a,"+-",Sa
print "b:",b,"+-",Sb
plt.xlim(0,0.3)
plt.ylim(0,2.1)
plt.show()
def plot231a():
plt.figure(figsize=(16,12))
ocd.plot_var(I,U,var,messung231)
plt.grid(True)
# -*- coding: utf-8 -*-
import sympy as sy
import ocd
import matplotlib.pyplot as plt
import numpy as np
############ Daten der Spulenwicklungen ##############
messung21=ocd.open_csv("21.csv")
messung21[2].Sx+=np.std(messung21[2].x[2:])
messung22=ocd.open_csv("22.csv")
Uz = sy.Symbol("Uz")
UA = sy.Symbol("UA")
N = sy.Symbol("N")
alpha = sy.Symbol("alpha")
var=[Uz,UA,N,alpha]
def plot21():
fig=plt.figure(figsize=(16,12))
t=ocd.plot_var(Uz,N,var,messung21,False)
#(a,b,Sa,Sb,Sy)=t
#textstr=ur"$a= (%.3f\pm %.3f)V$ $b= (%.3f\pm %.3f)\frac{V}{A}$"%(a,Sa,b,Sb)
#plt.title(ur"Bestimmung des ohmschen Widerstandes der Primärspule",fontsize=16, family= "monospace")
#props = dict(boxstyle="round", facecolor="wheat", alpha=0.5)
#plt.text(0.05, 0.9, textstr, fontsize=14, verticalalignment="top", bbox=props)
plt.xlabel(ur"Uz in V")
plt.ylabel(ur"N")
plt.title(ur"Messung 2.1 Zählrohrcharakteristik")
plt.grid(True)
plt.show()
def plot22():
fig=plt.figure(figsize=(16,12))
import ocd
Werte=ocd.open_csv("versuch80/tabelle2.csv")
s_unten=sy.Symbol("s_unten")
s_oben=sy.Symbol("s_oben")
r_unten=sy.Symbol("r_unten")
r_oben=sy.Symbol("r_oben")
I=sy.Symbol("I")
U=sy.Symbol("U")
alpha_offset=sy.Symbol("alpha_offset")
alpha=sy.Symbol("alpha")
B=0.78*10**(-3)*I
import scipy.constants as c
wz=c.e/c.m_e*B/(2*sy.pi)
test=ocd.eval_expr(wz,[s_unten,s_oben,r_oben,I,U,alpha_offset,alpha],Werte,"test")
print str(test)
b = np.linspace(0, 5, 100)
b2 = np.arange(0, 5, 1)
# plt.plot(b, normal_k(b), 'g--', linewidth=1)
plt.plot(b, normal_k_min(b), "g--", linewidth=1)
# plt.plot(b, normal_k_max(b), 'g--', linewidth=1)
# plt.plot(b, poisson(b), 'r--', linewidth=1)
plt.plot(b, poisson_min(b), "r--", linewidth=1)
# plt.plot(b, poisson_max(b), 'r--', linewidth=1)
plt.scatter(b2, normal_k(b2), marker="x", c="b")
plt.scatter(b2, poisson(b2), marker="x", c="b")
plt.title("Messung 3")
plt.xlabel("Anzahl der Pulse")
plt.ylabel("Wahrscheinlichkeit")
plt.grid(True)
plt.xlim(0, 5)
plt.ylim(0, 0.7)
plt.show()
messung1 = ocd.open_csv("messung1.csv")
n1 = messung1[0].x.magnitude
t1 = messung1[1].x.magnitude
plot_messung(n1, t1, 50, 0.1, 32, 1.5)
messung2 = ocd.open_csv("messung2.csv")
n2 = messung2[1].x.magnitude
t2 = messung2[0].x.magnitude
plot_messung(n2, t2, 12, 0.2, 13, 0.5)
plot_messung3()
# -*- coding: utf-8 -*-
import sympy as sy
import ocd
import matplotlib.pyplot as plt
import numpy as np
Werte=ocd.open_csv("messung.csv")
max_werte=[k.max() for k in Werte]
alpha=sy.Symbol("alpha")
U=sy.Symbol("U")
U_rot=sy.Symbol("U_rot")
U_gruen=sy.Symbol("U_gruen")
U_max=sy.Symbol("U_max")
U_rot_max=sy.Symbol("U_rot_max")
U_gruen_max=sy.Symbol("U_gruen_max")
var=[alpha,U,U_rot,U_gruen,U_max,U_rot_max,U_gruen_max]
cos2=sy.cos(alpha)**2
plt.figure()
t1=ocd.plot_var(alpha,U/U_max,var,Werte+max_werte,False)
t2=ocd.plot_var(alpha,U_rot/U_rot_max,var,Werte+max_werte,False)
t3=ocd.plot_var(alpha,U_gruen/U_gruen_max,var,Werte+max_werte,False)
x=np.linspace(-1.58,1.58,1000)
plt.plot(x,np.cos(x)**2)
#print "a: %.3f b: %.3f Sa: %.3f Sb: %.3f Sy: %.3f"%t1
#print "a: %.3f b: %.3f Sa: %.3f Sb: %.3f Sy: %.3f"%t2
#print "a: %.4f b: %.3f Sa: %.3f Sb: %.3f Sy: %.3f"%t3
plt.xlabel(ur"$\alpha$")
plt.ylabel(ur"$\frac{I}{I}$")
plt.title(ur"Bestätigung des Gesetzes $I=\cos(\alpha)^3 I$")
plt.grid(True)
plt.show()
import sympy as sy
import ocd
Werte=ocd.open_csv("tables/table1.csv")
t=sy.Symbol("t")
s1=sy.Symbol("s1")
s2=sy.Symbol("s2")
q=sy.Symbol("q")
f=(s1/t + s2/t )*q*sy.cos((s2/s1))
test=ocd.eval_expr(f,[t,s1,s2,q],Werte,"test")
print str(test)
# -*- coding: utf-8 -*-
import ocd
import numpy as np
import sympy as sy
import matplotlib.pyplot as plt
messung211=ocd.open_csv("211.csv")
messung212=ocd.open_csv("212.csv")
f = sy.Symbol("f")
Ue = sy.Symbol("Ue")
Ua = sy.Symbol("Ua")
var=[f,Ue,Ua]
f_=np.linspace(0,120000,1000)
def plot_ver():
fig=plt.figure(figsize=(16,12))
ocd.plot_var(f,Ua/Ue,var,messung212)
t=ocd.plot_var(f,Ua/Ue,var,messung211,True)
(a,b,Sa,Sb,Sy)=t
textstr=ur"$a= (%.3f\pm %.3f)V$ $b= (%.3f\pm %.3f)\frac{V}{A}$"%(a,Sa,b,Sb)
props = dict(boxstyle="round", facecolor="wheat", alpha=0.5)
plt.text(0.05, 0.9, textstr, fontsize=14, verticalalignment="top", bbox=props)
plt.ylim(0,100)
plt.xlabel("f in Hz")
plt.ylabel(ur"Linearer Verstärker")
plt.show()
##### Frequenzabhängiger Verstärker #####
messung213=ocd.open_csv("213.csv")
R1 = 4700
R2 = 470000
R3 = 47000
# -*- coding: utf-8 -*-
import sympy as sy
import ocd
import matplotlib.pyplot as plt
import numpy as np
############ Daten der Spulenwicklungen ##############
messung211=ocd.open_csv("211.csv")
U1 = sy.Symbol("U1")
U2 = sy.Symbol("U2")
I1 = sy.Symbol("I1")
I2 = sy.Symbol("I2")
var=[U1,U2,I1,I2]
def plot_primaer():
fig=plt.figure(figsize=(16,12))
t=ocd.plot_var(I1,U1,var,messung211,True)
(a,b,Sa,Sb,Sy)=t
textstr=ur"$a= (%.3f\pm %.3f)V$ $b= (%.3f\pm %.3f)\frac{V}{A}$"%(a,Sa,b,Sb)
plt.title(ur"Bestimmung des ohmschen Widerstandes der Primärspule",fontsize=16, family= "monospace")
props = dict(boxstyle="round", facecolor="wheat", alpha=0.5)
plt.text(0.05, 0.9, textstr, fontsize=14, verticalalignment="top", bbox=props)
plt.xlabel(ur"I in A")
plt.ylabel(ur"U in V")
plt.grid(True)
plt.show()
def plot_sekundaer():
fig=plt.figure(figsize=(16,12))
t=ocd.plot_var(I2,U2,var,messung211,True)
(a,b,Sa,Sb,Sy)=t
textstr=ur"$a= (%.3f\pm %.3f)V$ $b= (%.3f\pm %.3f)\frac{V}{A}$"%(a,Sa,b,Sb)
import ocd
import sympy as sy
import matplotlib.pyplot as plt
import numpy as np
ex1=ocd.open_csv("ex1.csv")
ex2=ocd.open_csv("ex2.csv")
ex3=ocd.open_csv("ex3.csv")
ex4=ocd.open_csv("ex4.csv")
ex5=ocd.open_csv("ex5.csv")
offsets=ocd.open_csv("offsets.csv")
e_=sy.Symbol("e_")
g_=sy.Symbol("g_")
xb=sy.Symbol("xb")
xs=sy.Symbol("xs")
xa=sy.Symbol("xa")
xl=sy.Symbol("xl")
xk=sy.Symbol("xk")
var=[e_,g_,xb,xs,xa,xl,xk]
e = e_ + xa + xs
g = g_ + xa + xb
b = e - g
beta = b / g
f = g*b/ (g+b)
def plot_gbe(ex,f_re):
plt.figure()
beta_r=np.linspace(0.2,5,1000)
