2021物理化学实验8：乙酸乙酯皂化反应动力学研究

#作者：王日睿 
#中国科学技术大学生命科学学院
#2021.12.23
#物理化学实验：实验08 乙酸乙酯皂化反应动力学研究
# Jupyter lab
 
import pandas as pd
import numpy as np
import matplotlib.pyplot as plt
plt.rcParams['font.sans-serif'] = ['SimHei']  # 用来正常显示中文标签
plt.rcParams['axes.unicode_minus'] = False  # 用来正常显示负号

# 乙酸乙酯初始时刻的电导率L_0
# 30度下：
L_0_1 = [1633, 1633, 1632]
L_0_1 = sum(L_0_1)/3
# 35度下：
L_0_2 = [1742, 1742, 1742]
L_0_2 = sum(L_0_2)/3

# 手动把数据复制到txt文件了，省点时间
L_t_30 = np.loadtxt('30.txt')
L_t_35 = np.loadtxt('35.txt')

# 实际反应温度
T_30 = np.loadtxt('Temperature30.txt')
T_35 = np.loadtxt('Temperature35.txt')

# 实际反应温度：
T1 = sum(T_30)/len(T_30)
T2 = sum(T_35)/len(T_35)
print("第一组实际的反应温度平均值为%.2f"%T1)
print("第二组实际的反应温度平均值为%.2f"%T2)

>第一组实际的反应温度平均值为29.19
>第二组实际的反应温度平均值为34.17

# 反应时间（采样间隔为10s）
t0 = 10
time1 = np.arange(0, L_t_30.shape[0]*t0, t0)
time2 = np.arange(0, L_t_35.shape[0]*t0, t0)

# 用于拟合的y值
y1 = [(L_0_1-L_t_30[i])/time1[i] for i in range(time1.shape[0])]
y2 =