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)