import Praktikum
import numpy as np
from pylab import *


font = {'family' : 'DejaVu Sans',
        'weight' : 'normal',
        'size'   : 14}

rc('font', **font)


tgll=[(7.18-0.37)/10,(7.52-0.71)/10,(7.49-0.68)/10,(7.83-1.01)/10]
tggg=[(7.62-0.82)/10,(7.62-0.82)/10]
tss=[(8.25-2.63)*2,(8.2-2.67)*2,(8.30-2.65)*2,(8.25-2.67)*2]

tgl=np.mean(tgll)
#ss=0
#n=len(tgll)
#for j in range(0,n):
#    ss=ss+(tgll[j]-tgl)**2
#otgl=np.sqrt(ss/((n-1)*n))

#tgg=np.mean(tggg)
#ss=0
#n=len(tggg)
#for j in range(0,n):
#    ss=ss+(tggg[j]-tgg)**2
#otgg=np.sqrt(ss/((n-1)*n))

tgg=0.608
otgg=0.01/np.sqrt(6)
otgl=0.01/np.sqrt(6)

ts=np.mean(tss)
ss=0
n=len(tss)
for j in range(0,n):
    ss=ss+(tss[j]-ts)**2
ots=np.sqrt(ss/((n-1)*n))

fgl=1/tgl
fgg=1/tgg
fs=1/ts
ofgl=otgl/(tgl**2)
ofgg=otgg/(tgg**2)
ofs=ots/(ts**2)

k=(fgg**2-fgl**2)/(fgg**2+fgl**2)

ok=4/((fgl**2+fgg**2)**2)*np.sqrt(ofgl**2*fgg**2+ofgg**2*fgl**2)

r=0.6
C=0.0000047
L=0.0022

t=ts/(2*np.pi)*(np.pi/2-np.arctan(k/r*np.sqrt(L/C)))

print(tgl)
print(tgg)
print(ts)
print(otgl)
print(ots)
print(k)
print(ok)
print(t)


data = Praktikum.lese_lab_datei('lab/s1.lab')
t = data[:,1]
ub = data[:,2]
ua = data[:,3]

voltageError = [0.01]*len(ub)
offset = Praktikum.gewichtetes_mittel(np.array(ub), np.array(voltageError))[0]
ub = ub - offset

offset = Praktikum.gewichtetes_mittel(np.array(ua), np.array(voltageError))[0]
ua = ua - offset

figure()
plot(t,ub,',')
plot(t,ua,',')
points = np.linspace(0,0.02,1000)
plot(points, [0 for xi in points])