Posts
学习与思考
-
Python拟合吸附等温线并计算吸附焓
Python, Scipy, Matplotlib, Csv, Pandas ·前一阵子帮人用C++写了一个用以计算吸附焓的简单小程序，仅实现了得到吸附等温线方程后计算吸附焓的功能。最近重拾Python语言，在尝试了Python的各类库之后对其爱不释手。遂尝试着将这个程序用Python完善了一下，实现了数据读入，拟合，计算吸附焓一系列流程。
作为一个新手，写程序的过程依然很揪心。好在换用了Anaconda后的Jupyter Notebook十分友好，可以在每个cell中直接运行查看结果进行调试。在经历了不断地与错误信息斗争后，终于磕磕绊绊地将想要的功能实现了。所以在这里记录一下学习到的一些语法知识以及程序目前的缺陷。
程序内容
#拟合BET数据中的吸附曲线并计算吸附焓 #Version 2.00 #以实现对数据的读取，转化单位和langmuir方程拟合并做图 #对BET软件直接生成的xls文件尚不能直接打开，需要另存为xlsx格式 #Author: lewisbase #Date: 2019.04.08 import csv import math import pandas as pd import numpy as np import matplotlib.pyplot as plt from matplotlib import ticker from matplotlib import rcParams from scipy.optimize import leastsq GRID=30000 R=8.314510e-3 ############################################################################################### def Langmuir(c): b=c return y-qmax*b*x/(1+b*x) def Langmuirfit(x): return qmax*b*x/(1+b*x) def Freundlich(f):...