2021 物理化学实验1：恒温槽的装配与性能测定

#作者：王日睿 
#中国科学技术大学生命科学学院
#2021.11.3
#物理化学实验：恒温槽的装配与性能测定

import numpy as np
import matplotlib.pyplot as plt
import os

def readfile(filename):
    x, y = [], [] 
    file = open(filename, 'r')
    count = 0
    for line in file.readlines():
        count = count + 1
        if count >= 4:
            x.append(int(line[0:4]))
            y.append(float(line[7:len(line)]))
    return x, y


FILE_PATH = "C:\\Users\\29070\\Desktop\\我的笔记\\物化实验\\恒温槽的装配与性能测定\\wangrirui-shiyang1\\wangrirui-shiyang1"
FILE_NAME = os.listdir(FILE_PATH)

time = []
temperature = [] # 按文件名依次储存在time[],和temperature[]里
for file in FILE_NAME:
    x, y = readfile(FILE_PATH + r'\\' + file)
    time.append(x)
    temperature.append(y)

# 作图,观察并切片，舍去无效数据
plt.figure(1)
fig, ax = plt.subplots(4, 2, figsize=(10, 10), dpi=400)

plt.rcParams['font.sans-serif'] = ['SimHei']  # 用来正常显示中文标签
plt.rcParams['axes.unicode_minus'] = False  # 用来正常显示负号

ax[0][0].set_xlabel("time/min")
ax[0][0].set_ylabel("temperature")
ax[0][0].set_title("DTC-125V")
ax[0][0].plot(time[0][30:], temperature[0][30:], color="r")
time[0] = time[0][30:]
temperature[0] = temperature[0][30:]

ax[0][1].set_xlabel("time/min")
ax[0][1].set_ylabel("temperature")
ax[0][1].set_title("DTC-200V")
ax[0][1].plot(time[1][100:1000], temperature[1][100:1000], color="b")
time[1] = time[1][100:1000]
temperature[1] = temperature[1][100:1000]

ax[1][0].set_xlabel("time/min")
ax[1][0].set_ylabel("temperature")
ax[1][0].set_title("继电器-125V")
ax[1][0].plot(time[2][80:],temperature[2][80:],color="g")
time[2] = time[2][80:]
temperature[2] = temperature[2][80:]

ax[1][1].set_xlabel("time/min")
ax[1][1].set_ylabel("temperature")
ax[1][1].set_title("继电器-200V")
ax[1][1].plot(time[3],temperature[3],color="y")

ax[2][0].set_xlabel("time/min")
ax[2][0].set_ylabel("temperature")
ax[2][0].set_title("通冷却水-DTC-200V")
ax[2][0].plot(time[4][760:],temperature[4