#作者:王日睿
#中国科学技术大学生命科学学院
#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