2021物理化学实验6：电池电动势的测定及其应用

#作者：王日睿 
#中国科学技术大学生命科学学院
#2021.12.9
#物理化学实验：实验06 电池电动势的测定及其应用
# Jupyter lab

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

## raw data
E_30 = [457.759, 457.573, 457.420]
E_35 = [453.517, 453.453, 453.451]
E_40 = [451.238, 450.889, 450.566]
E_45 = [447.726, 447.288, 447.137]
E_cu = [427.891, 428.046, 428.400]

E_t = [sum(E_30)/3, sum(E_35)/3, sum(E_40)/3, sum(E_45)/3]

T = [30,35,40,45]

T = [t+273.15 for t in T]
T = np.array(T)

from scipy import optimize

def func(x,k,b):
    return k*x+b

k, b = optimize.curve_fit(func,T,E_t)[0]

X=np.arange(T.min()-3,T.max()+3,0.1)
Y=func(X,k,b)

def get_R_square(x, y):
    y_mean = np.mean(y)
    y_predict = [k * e + b for e in x]
    TSS = [(x-y_mean)**2 for x in y]
    RSS &