大学物理实验教学中基于MATLAB的液体表面张力系数测定的实验数据处理.docVIP

大学物理实验教学中基于MATLAB的液体表面张力系数测定的实验数据处理摘要：液体表面张力系数的测定是大学物理实验教学中的一个重要实验。本文将利用MATLAB编程来处理实验数据，得到液体表面张力系数，和传统数据处理方法相比，用MATLAB处理实验数据能有效避免手工处理所带来的误差，适合在实验教学中使用，而且能提高大学本科学生利用计算机语言编程处理实验数据的能力，可以取得良好的教学效果。关键词：大学物理实验教学；MATLAB软件；表面张力系数中图分类号：G642.3 文献标志码：A 文章编号：1674-9324（2015）33-0237-02液体表面张力系数的测定是大学物理实验教学中的一个重要实验。液体表面张力系数是反映液体性质的一个非常重要的物理量，我们在众多自然科学领域中有着重要的实际应用，比如我们在工农业生产上的浮选技术和液体输运技术等都要对表面张力进行研究。对液体表面张力系数进行测量有多种方法，拉脱法、毛细管法和液滴测重法等，大学物理实验中一般采用拉脱法测液体表面张力系数。[1，2]液体表面张力系数的测定作为基本的大学物理实验之一，它的实验数据处理方法很重要，但是一般方法数据处理多存在不方便核查校对，耗时费力，同时也无法找出误差较大甚至错误的数据点等问题，而MATLAB计算机软件是把符号运算、数学计算和图形处理、计算机语言编程等功能集中于一身的计算机语言，它经常被应用于科学研究中。[3-6]运用MATLAB计算机语言来处理大学物理实验数据我们只需编写一些简单的像通常算术的简单程序，运行后就可得到我们设想的结果，我们运用MATLAB语言编程既克服了最小二乘法人工计算量大的不足之处，又使麻烦无味的数学计算转变成一种简洁的可视化操作流程，而且能非常精确地标记大学物理实验的数据点和绘制出实验数据的拟合曲线。[7，8]本文利用MATLAB语言编程来处理实验数据，得到液体表面张力系数，和传统数据处理方法相比，用MATLAB处理实验数据能有效避免手工处理所带来的误差，适合在实验教学中使用，而且能提高学生利用计算机语言编程处理实验数据的能力。一、液体表面张力系数的测定的实验原理液体表面张力系数的测定是大学物理实验教学中的一个重要实验。我们在本实验用到的器材有硅压阻力敏传感器、铁架台、铝合金吊片、吊盘、0.5g片码、玻璃器皿等。表面张力系数为表示液体表面性质的一个重要物理参考量，其中拉脱法是测算表面张力系数的最常用方法。通过测已有周长的金属片从待测的液体的表面脱落时需要的拉力，来测算出液体的表面张力系数的物理实验方法，我们称之为拉脱法。正如图1所示，假设在液体表面层中有一分界线MN，这两部分表面层中的分子存在相互作用引力（即表面张力）1和2，它们大小相等、方向相反、垂直于MN并沿液体表面分别作用在表面层相互接触部分。表面张力正比于表面分界线的长度f=σl，σ为比例系数，为表面张力特征的一个物理量，称为液体的表面张力系数。实验中，测量一铝合金吊片从待测液体表面拉脱时需要的力。拉脱前瞬间，可以近似地认为圆形吊片脱离液体表面前的瞬间受力为F1=σ?π（D1+D2）+mg，F1是传感器挂钩受到的拉力，D1、D2分别为圆形吊片底面的内外直径，mg为吊片的重力。拉脱后传感器挂钩受拉力为F2=mg。硅压阻力敏传感器由弹性梁和贴在梁上的传感器芯片组成，其中芯片由四个硅扩散电阻集成一个非平衡电桥。当外界压力作用于金属梁上时，在压力作用下，电桥失去平衡，此时将有电压信号输出，输出电压大小与所加外力成正比。即U=kF，k是硅压阻力敏传感器的灵敏度，其值可由实验定标测出，F是外力大小，U是传感器输出电压大小。拉脱前后输出电压的改变量ΔU可表示为ΔU=U1-U2=kF1-kF2=kσπ（D1+D2）。液体表面张力系数可表示为σ=二、基于MATLAB的液体表面张力系数测定的实验数据处理我们所测的液体为自来水。首先，金属环的外径和内径经过三次测量然后求平均值，得到外径D1=3.491cm和内径D2=3.323cm。然后测量硅压阻力敏传感器电压和重力的对应值。然后利用MATLAB语言编程的最小二乘法拟合实验数据，求出传感器灵敏度，程序编写如下：%测量传感器灵敏度系数k的实验数据处理clearg=9.8M=[0.500 1.000 1.500 2.000 2.500 3.000 3.500]；%重力的实验值U=[17.5 35.0 52.5 70.0 88.2 105.6 122.8]；%电压的测量值p=polyfit（M，U，1）；%用最小二乘法数值拟合的方程i=linspace（0.0，4.0，100）；%在0.0-4.0之间等间隔地取100个点z=polyval（p，i）；%作100个点的拟合运算plot（M，U，*，i，z，b）%以‘*’标实验值，以蓝色画拟合曲线x