Matlab×物理化学实验:蔗糖水解反应速率常数测定
Matlab×物理化学实验:蔗糖水解反应速率常数测定
- 摘要
- 理论依据
- 一、蔗糖水解
- 化学反应式
- 半衰期
- 一级反应有如下特点
- 二、旋光度α\alphaα
- 求解思想
- 核心语句:
- 绘图相关语句:
- 程序实现
- 程序内容
- 范例
- 输入数据
- 调用函数
- 执行结果
- 写在最后
- 参考文献
版权声明:
如要转载,请与我联系。
摘要
对于一个化学反应,用普通的方法不太容易得知其反应进度和反应速率。该次举例的实验采用对反应溶液旋光度的测量,算出对应时间的旋光度,采用作图法画出直线,用直线斜率求得反应速率常数kkk。该文参照物理化学的知识原理,利用实验数据,应用Matlab求解,最后得出反应速率常数kkk。
关键词:乙醇水溶液、比重瓶法、偏摩尔体积
理论依据
一、蔗糖水解
化学反应式
蔗糖的分解反应属于一级反应
蔗糖能在酸性催化下能够水解为葡萄糖和果糖,方程式如下:
根据一级反应的特点,反应速率只与浓度有关,因此反应速率的公式可以写为:
v=−dcdt=ktv=\frac{\mathrm{-d}c}{\mathrm{d}t}=kt v=dt−dc=kt
对该公式积分并带入初始值可以得到:
k=1tlnc0ck=\frac{1}{t}ln\frac{c_0}{c} k=t1lncc0
其中,vvv是反应速率;kkk是反应速率常数;ttt是反应的时间;ccc是时间为ttt时的浓度,c0c0c0是起始浓度。
半衰期
当反应进行到一半,即反应物消耗一半,所需要的时间称为半衰期t1/2t_{1/2}t1/2,(要注意半衰期与起始浓度c0c_0c0无关)半衰期表达式如下:
t1/2=1klnc0c/2=ln2kt_{1/2}=\frac{1}{k}ln\frac{c_0}{c/2}=\frac{ln2}{k} t1/2=k1lnc/2c0=kln2
一级反应有如下特点
- kkk的数值与浓度无关,其量纲为时间的负一次幂,常用单位是s−1s^{-1}s−1、min−1min^{-1}min−1、h−1h^{-1}h−1
- 半衰期与反应的起始浓度无关(在计算中被消去了)
- 以lnclnclnc对ttt作图可以得到一条直线,该直线的斜率为−k-k−k,截距为lnc0lnc_0lnc0
因此,可由作图法可以得到直线斜率,以此计算反应常数kkk
二、旋光度α\alphaα
蔗糖及其水解产物葡萄糖和果糖都有不对称的碳原子,因此有旋光度。蔗糖和葡萄糖是右旋物质(旋光度α\alphaα>0),果糖是左旋物质(旋光度α\alphaα<0)。
由于果糖比旋光度的绝对值大于葡萄糖比旋光度的绝对值,因此在反应结束后,溶液的旋光度一定是负数。
我们可以通过测量反应溶液的旋光度,进而知晓该反应的反应进度和反应速率。
并且,在其他条件没有改变的情况下(如探测灯的波长λλλ、液层厚度lll、实验温度ttt),物质的旋光度α\alphaα与反应物质的浓度ccc成正比:
ct∼[αt−α∞]c_t\sim[\alpha _t-\alpha_\infty] ct∼[αt−α∞]
得到下式:
k=1tlnc0c=1tlnα0−α∞αt−α∞k=\frac{1}{t}ln\frac{c_0}{c}=\frac{1}{t}ln\frac{\alpha _0-\alpha_\infty}{\alpha _t-\alpha_\infty} k=t1lncc0=t1lnαt−α∞α0−α∞
求解思想
该实验分为两个部分,
- 旋光度αt\alpha_tαt和时间ttt的测量;
- 反应完全结束时旋光度α∞\alpha_\inftyα∞的测量
我们进行该试验之后,对数据进行处理,
将得到的ln(αt−α∞)−tln(\alpha_t-\alpha_\infty)-tln(αt−α∞)−t进行1次线性拟合,用直线斜率求反应速率kkk
核心语句:
polyfit(x,y,n) #最小二乘法对[x,y]拟合进行n次拟合
绘图相关语句:
plot(x,y)
lengend('String1','String2')
title('String')
xlabel('String')
ylabel('String')
gird on
程序实现
以下为matlab的.m函数文件,注意使用时候需要将文件名和函数名设置相同。
程序内容
%{filename:Rate_Hydrolysis.mauthor:Lyn#reocreate:2020.10.16#可以根据需要修改#如有二次发送的必要#请留下修改者和修改日期:)#感激您的真诚#
%}
%程序开始
function Rate_Hydrolysis(x,y)
p=polyfit(x,y,1); %直接进行拟合
disp('该直线方程为:')
ln_ataq=poly2str(p,'t')
q=polyder(p); %对多项式进行求导
k=-q; %通过斜率算得反应速率k
fprintf('斜率k=%0.6fmin^-1\n',k)
t12=0.6932/k; %求解半衰期t_1/2
fprintf('半衰期t1/2=%.6fmin\n',t12)
xr=0:0.5:60; %反应进行60min时基本完毕
yl=polyval(p,xr); %画出拟合曲线在0~50min的曲线
%画图的相关指令
plot(x,y,'bo',xr,yl,'r-')
title('蔗糖溶解')
xlabel('反应时间t/min')
ylabel('lnαt-aq')
grid on
%程序结束
范例
注:数据是我为举例随机顺序生成的,与真实数据差距巨大。为避免不必要的麻烦,请不要照搬!
以部分数据为例
时间t/mint/mint/min | ln(αt−α∞)ln(\alpha_t-\alpha_\infty)ln(αt−α∞) |
---|---|
21.12 | 1.90 |
23.65 | 1.77 |
26.20 | 1.65 |
28.34 | 1.58 |
30.67 | 1.45 |
34.23 | 1.30 |
35.45 | 1.25 |
36.53 | 1.20 |
36.23 | 1.18 |
输入数据
x=[21.12,23.65,26.20,28.34,30.67,34.23,35.45,36.53,36.23]
y=[1.90,1.77,1.65,1.58,1.45,1.30,1.25,1.20,1.18]
调用函数
Rate_Hydrolysis(x,y)
执行结果
>> Rate_Hydrolysis(x,y)
该直线方程为:
ln_ataq =' -0.04564 t + 2.857'
斜率k=0.045640min^-1
半衰期t1/2=15.188295min
写在最后
本次实例比较简单,用一次拟合即可达成。
如果有人看的话,祝你变得更强,下次见 : )
参考文献
·物理化学.上册/傅献彩等编.——5版——北京:高等教育出版社,2005.7
·物理化学实验(工科类专业用)/李晔,韦美菊主编.——北京:化学工业出版社,2013.3
·MATLAB与数学实验/艾冬梅等编著.——2版本.——北京:机械工业出版社,2014.6
Matlab×物理化学实验:蔗糖水解反应速率常数测定相关推荐
- Matlab×物理化学实验:乙醇水溶液偏摩尔体积测定
MATLAB×物理化学实验:乙醇水溶液偏摩尔体积测定 摘要 理论依据 一.比重瓶法测密度 二.多组分液体的体积 加和公式 求解思想 核心语句: 绘图相关语句: 程序实现 程序内容 范例 输入数据 调用 ...
- 计算机辅助物理化学实验 唐典勇课后答案,计算机辅助物理化学实验(第2版)...
第1章 绪论 1.1 物理化学实验的目的和要求 1.1.1 学习物理化学实验课程的目的 1.1.2 实验前的准备 1.1.3 实验注意事项 1.1.4 实验报告要求 1.1.5 实验室规 ...
- 计算机量子化学计算焓变,物理化学实验
绪论 1.1 物理化学实验目的与要求 1.2 误差及数据处理 1.3 计算机处理物理化学实验数据的方法 1.4 学生实验规则与安全 实验内容 2.1 热力学 实验1 恒温槽控温精度的研究 实验2 燃烧 ...
- 2021物理化学实验4:液体饱和蒸气压的测定
#作者:王日睿 #中国科学技术大学生命科学学院 #2021.11.20 #物理化学实验:实验04 液体饱和蒸气压的测定 # Jupyter labimport numpy as np import m ...
- 2021物理化学实验3:燃烧热的测定
#作者:王日睿 #中国科学技术大学生命科学学院 #2021.11.12 #物理化学实验:实验03 燃烧热的测定 # Jupyter labimport numpy as np import matpl ...
- matlab物理实验数据处理,利用Matlab对物理实验数据进行处理
利用Matlab对物理实验数据进行处理 李凌云;王海军;韩之乃 [期刊名称]<南京晓庄学院学报> [年(卷),期]2005(021)005 [摘要]本文以热电偶的定标实验.二极管伏安特性实 ...
- 2021物理化学实验5:双液系的气液平衡相图
#作者:王日睿 #中国科学技术大学生命科学学院 #2021.12.1 #物理化学实验:实验04 液体饱和蒸气压的测定 # Jupyter lab# 1. 绘制工作曲线 # 原始数据: propanol ...
- matlab操作入门实验报告,MATLAB基本操作实验报告.doc
MATLAB基本操作实验报告MATLAB基本操作实验报告 南昌航空大学 数学与信息科学学院 实 验 报 告 课程名称: 实验名称: 实验类型: 验证性 综合性□ 设计性□ 实验室名称: 班级学号: 学 ...
- matlab数学实验 胡良剑,matlab数学实验习题全部答案(胡良剑)
matlab数学实验习题全部答案(胡良剑) 1数学实验答案%Page20,ex1 (5) 等于[exp(1),exp(2);exp(3),exp(4)] (7) 3=1*3, 8=2*4 (8) a ...
最新文章
- Solr学习总结(二)Solr的安装与配置
- python数据挖掘工程师待遇_数据挖掘工程师岗位的具体职责
- 摘自《Java工程师成神之路》2018修订版,自我勉励
- Oh my God, Swagger API文档竟然可以这样写?
- 中公学python要多久_学完风变编程我才知道:一天的工作,五分钟就能完成!
- Java Protected 解读
- import com.google.common.* 出错,找不到
- android bin目录是什么,android bin目录下的.ap_是神马文件?
- sqlplus中设置在屏幕中上不打印出输出
- C# Winform关于控件TabControl闪烁的问题
- linux shell通过i2cget命令获取I2C器件寄存器的值
- Open Robotics 十周年(2012-2022)
- oracle10.2.0.4 dbca,10.2.0.4 DBCA problem :Error securing Database Control,...
- 查看电脑无线网 wifi密码
- java远控_利用Java实现远程控制
- opencv给图片加滑动条调节窗口大小
- 促销活动的2个环节与7个关键点
- 字符编码笔记:ASCII,Unicode 和 UTF-8(转帖、留着自己学习)
- ABAP subroutine 的定义和使用
- ANSA二次开发——BCGUI简介(3)
热门文章
- 使用opencv和C++实现多焦距图像合成
- PyQt5最全73 布局之addStretch设置布局的伸缩量
- 试题 算法训练 24点
- [AHK]输入法状态提示,中文状态提示“中”,英文状态提示“EN”[转]
- Spark Sql 聚合
- Manifest merger failed : uses-sdk:minSdkVersion 19 cannot be smaller than version 21 declared in lib
- 全面解析四大主流音频技术
- Python飞机大战(究极迫害版)
- mysql数据表分析统计记录
- 一个公众号绑定多个商户支付