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=t1​lncc0​​
其中,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​=k1​lnc/2c0​​=kln2​

一级反应有如下特点

  1. kkk的数值与浓度无关,其量纲为时间的负一次幂,常用单位是s−1s^{-1}s−1、min−1min^{-1}min−1、h−1h^{-1}h−1
  2. 半衰期与反应的起始浓度无关(在计算中被消去了)
  3. 以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=t1​lncc0​​=t1​lnαt​−α∞​α0​−α∞​​


求解思想

该实验分为两个部分,

  1. 旋光度αt\alpha_tαt​和时间ttt的测量;
  2. 反应完全结束时旋光度α∞\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×物理化学实验:蔗糖水解反应速率常数测定相关推荐

  1. Matlab×物理化学实验:乙醇水溶液偏摩尔体积测定

    MATLAB×物理化学实验:乙醇水溶液偏摩尔体积测定 摘要 理论依据 一.比重瓶法测密度 二.多组分液体的体积 加和公式 求解思想 核心语句: 绘图相关语句: 程序实现 程序内容 范例 输入数据 调用 ...

  2. 计算机辅助物理化学实验 唐典勇课后答案,计算机辅助物理化学实验(第2版)...

    第1章  绪论 1.1  物理化学实验的目的和要求 1.1.1  学习物理化学实验课程的目的 1.1.2  实验前的准备 1.1.3  实验注意事项 1.1.4  实验报告要求 1.1.5  实验室规 ...

  3. 计算机量子化学计算焓变,物理化学实验

    绪论 1.1 物理化学实验目的与要求 1.2 误差及数据处理 1.3 计算机处理物理化学实验数据的方法 1.4 学生实验规则与安全 实验内容 2.1 热力学 实验1 恒温槽控温精度的研究 实验2 燃烧 ...

  4. 2021物理化学实验4:液体饱和蒸气压的测定

    #作者:王日睿 #中国科学技术大学生命科学学院 #2021.11.20 #物理化学实验:实验04 液体饱和蒸气压的测定 # Jupyter labimport numpy as np import m ...

  5. 2021物理化学实验3:燃烧热的测定

    #作者:王日睿 #中国科学技术大学生命科学学院 #2021.11.12 #物理化学实验:实验03 燃烧热的测定 # Jupyter labimport numpy as np import matpl ...

  6. matlab物理实验数据处理,利用Matlab对物理实验数据进行处理

    利用Matlab对物理实验数据进行处理 李凌云;王海军;韩之乃 [期刊名称]<南京晓庄学院学报> [年(卷),期]2005(021)005 [摘要]本文以热电偶的定标实验.二极管伏安特性实 ...

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

    #作者:王日睿 #中国科学技术大学生命科学学院 #2021.12.1 #物理化学实验:实验04 液体饱和蒸气压的测定 # Jupyter lab# 1. 绘制工作曲线 # 原始数据: propanol ...

  8. matlab操作入门实验报告,MATLAB基本操作实验报告.doc

    MATLAB基本操作实验报告MATLAB基本操作实验报告 南昌航空大学 数学与信息科学学院 实 验 报 告 课程名称: 实验名称: 实验类型: 验证性 综合性□ 设计性□ 实验室名称: 班级学号: 学 ...

  9. matlab数学实验 胡良剑,matlab数学实验习题全部答案(胡良剑)

    matlab数学实验习题全部答案(胡良剑) 1数学实验答案%Page20,ex1 (5) 等于[exp(1),exp(2);exp(3),exp(4)] (7) 3=1*3, 8=2*4 (8) a ...

最新文章

  1. Solr学习总结(二)Solr的安装与配置
  2. python数据挖掘工程师待遇_数据挖掘工程师岗位的具体职责
  3. 摘自《Java工程师成神之路》2018修订版,自我勉励
  4. Oh my God, Swagger API文档竟然可以这样写?
  5. 中公学python要多久_学完风变编程我才知道:一天的工作,五分钟就能完成!
  6. Java Protected 解读
  7. import com.google.common.* 出错,找不到
  8. android bin目录是什么,android bin目录下的.ap_是神马文件?
  9. sqlplus中设置在屏幕中上不打印出输出
  10. C# Winform关于控件TabControl闪烁的问题
  11. linux shell通过i2cget命令获取I2C器件寄存器的值
  12. Open Robotics 十周年(2012-2022)
  13. oracle10.2.0.4 dbca,10.2.0.4 DBCA problem :Error securing Database Control,...
  14. 查看电脑无线网 wifi密码
  15. java远控_利用Java实现远程控制
  16. opencv给图片加滑动条调节窗口大小
  17. 促销活动的2个环节与7个关键点
  18. 字符编码笔记:ASCII,Unicode 和 UTF-8(转帖、留着自己学习)
  19. ABAP subroutine 的定义和使用
  20. ANSA二次开发——BCGUI简介(3)

热门文章

  1. 使用opencv和C++实现多焦距图像合成
  2. PyQt5最全73 布局之addStretch设置布局的伸缩量
  3. 试题 算法训练 24点
  4. [AHK]输入法状态提示,中文状态提示“中”,英文状态提示“EN”[转]
  5. Spark Sql 聚合
  6. Manifest merger failed : uses-sdk:minSdkVersion 19 cannot be smaller than version 21 declared in lib
  7. 全面解析四大主流音频技术
  8. Python飞机大战(究极迫害版)
  9. mysql数据表分析统计记录
  10. 一个公众号绑定多个商户支付