2021物理化学实验5：双液系的气液平衡相图

#作者：王日睿 
#中国科学技术大学生命科学学院
#2021.12.1
#物理化学实验：实验04 液体饱和蒸气压的测定
# Jupyter lab

# 1. 绘制工作曲线
# 原始数据：
propanold = [0, 0.5145,1.0230,1.5360,2.0142,2.5632,3.0618,3.5026,4.0367,4.5147, 5] # 正丙醇质量
total = [5, 4.9609,5.0349,4.9806,5.1029,5.1980,5.0913,4.9578,5.0469,4.9968, 5] # 总质量
reflac1 = [1.3325, 1.3423,1.3504,1.3572,1.3620,1.3672,1.3723,1.3767,1.3794,1.3820, 1.3831]
reflac2 = [1.3325, 1.3423,1.3504,1.3571,1.3620,1.3672,1.3723,1.3767,1.3794,1.3820, 1.3831]

# 初步处理
reflac = [(reflac1[i]+reflac2[i])/2 for i in range(len(reflac1))]
concentrate = [propanold[i]/total[i] for i in range(len(total))]

# 拟合：
# 选择2-4次的函数拟合
from scipy import optimize
def fun_2(x, a, b, c):
    return a*x**2 +b*x +c

def fun_3(x, a, b, c,d):
    return a*x**3 +b*x**2 +c*x+d

def fun_4(x, a,b,c,d,e):
    return a*x**4 +b*x**3 +c*x**2+d*x+e

parameters_2 = optimize.curve_fit(fun_2, reflac, conce