if line[0] <> linekick:
# print line
fields = line.strip().split(separator)
for i, number in enumerate(fields):
number = number.replace(",", ".")
data[i].append(number)
for i in xrange(len(dtype)):
data[i] = cast[dtype[i]](data[i])
return N.rec.array(data, dtype=dtype)
mydescr = N.dtype([("Zeit", "int32"), ("Ofen", "float"), ("Mitte", "float"), ("Halb", "float"), ("Ausen", "float")])
myrecarray = read_array("Stahlblock_150_C.ASC", mydescr)
timeline = myrecarray["Zeit"]
Temperatur = savitzky_golay(myrecarray["Mitte"], window_size=7, order=3)
Temperatur_Ofen = savitzky_golay(myrecarray["Ofen"], window_size=7, order=3)
Temperatur_Ofen_2 = savitzky_golay(myrecarray["Ausen"], window_size=17, order=6)
Temperatur_Ofen_3 = savitzky_golay(myrecarray["Halb"], window_size=17, order=6)
# Ofen = savitzky_golay(myrecarray['Ofen'], window_size=17, order=6)
# p.plot(timeline,Temperatur,lw=2)
# p.plot(timeline,75*1E-8*(Temperatur_Ofen-Temperatur)**4,lw=3)
# p.plot(timeline,Temperatur_Ofen_2-Temperatur,lw=2)
# p.plot(timeline,150-Temperatur,lw=2)
p.plot(timeline, Temperatur_Ofen_2)
p.plot(timeline, Temperatur_Ofen)
p.plot(timeline, Temperatur)
myrecarray = read_array("Stahlblock_150_C2.ASC", mydescr)
timeline = myrecarray["Zeit"]
from matplotlib.mlab import normpdf
import pylab as p
from scipy import *
from plot_XKCD import *
import pylab
from AuswertH import *
# eigene Unterprogramme
from filters import savitzky_golay
mydescr = N.dtype([('Zeit', 'int32'), ('Ofen', 'float'), ('Mitte', 'float'), ('Halb', 'float'),('Ausen', 'float')])
myrecarray = read_array('M121106A.ASC', mydescr,pointch=True)
timeline = myrecarray['Zeit']
Temperatur = savitzky_golay(myrecarray['Mitte'], window_size=17, order=6)
Temperatur1 = savitzky_golay(myrecarray['Halb'], window_size=17, order=6)
Temperatur2 = savitzky_golay(myrecarray['Ausen'], window_size=17, order=6)
Ofen = savitzky_golay(myrecarray['Ofen'], window_size=17, order=6)
ax = p.axes()
np.random.seed(0)
#x = np.linspace(0, 10, 100)
ax.plot(timeline, Ofen, 'b', lw=1, label='Ofen')
ax.plot(timeline, Temperatur1, 'r', lw=2, label='Probe')
#ax.plot(x, -np.cos(x) * np.exp(-0.1 * (x - 5) ** 2), 'r', lw=1, label='damped cosine')
ax.set_title('Zeitlicher Temperaturverlauf')
ax.set_xlabel('Zeit ')
ax.set_ylabel('Temperatur')
for line in open(filename, 'r'):
if line[0]<>linekick:
# print line
fields = line.strip().split(separator)
for i, number in enumerate(fields):
number=number.replace(',', '.')
data[i].append(number)
for i in xrange(len(dtype)):
data[i] = cast[dtype[i]](data[i])
return N.rec.array(data, dtype=dtype)
mydescr = N.dtype([('Zeit', 'int32'), ('Ofen', 'float'), ('Mitte', 'float'), ('Halb', 'float'),('Ausen', 'float')])
myrecarray = read_array('M081021A.ASC', mydescr)
timeline = myrecarray['Zeit']
Temperatur = savitzky_golay(myrecarray['Mitte'], window_size=17, order=6)
Ofen = savitzky_golay(myrecarray['Ofen'], window_size=17, order=6)
mydescr = N.dtype([('Zeit', 'int32'), ('Ofen', 'float'), ('Mitte', 'float'), ('Halb', 'float'),('Ausen', 'float')])
myrecarray2 = read_array('M081023A.ASC', mydescr)
timeline2 = myrecarray2['Zeit']
Temperatur2 = savitzky_golay(myrecarray2['Mitte'], window_size=17, order=6)
Ofen2 = savitzky_golay(myrecarray2['Ofen'], window_size=17, order=6)
timeline_d = N.diff(timeline)
Temperatur_d = N.diff(Temperatur)
timeline_d2 = N.diff(timeline2)
Temperatur_d2 = N.diff(Temperatur2)
