Matlab软件在牛顿环实验数据处理方法上的创新.doc

Matlab软件在牛顿环实验数据处理方法上的创新 徐少刚 夏雪琴 大学物理实验学A12船舶2 （浙江海洋学院，浙江 舟山 316022） 摘要：本文运用了用Matlab软件的功能函数对牛顿环测量平凸透镜的曲率半径实验的数据进行处理，简化了数据处理过程，提高了实验结果的准确度；同时也提高了我们大学生应用软件的能力，激发了学习兴趣，提高了教学质量。 关键词：牛顿环；Matlab软件；最小二乘法；数据处理；创新 牛顿环实验是一个古老而又非常重要的光学实验，一般用它来测量平凸透镜的曲率半径 十分方便，这个实验在我校已开设许多年了。目前数据处理方法主要有：逐差法和最小二乘法，用人工计算既繁琐又篇幅大，占用了我们学生较多的学习时间，并且往往运算容易出差错，得不到十分理想的实验结果，这样必然会影响我们学生的学习积极性和教学效果。随着计算机软件的迅速发展，如用Matlab软件的功能函数进行最小二乘法曲线拟合，得到平凸透镜的曲率半径，而且还能进行误差理论分析，一方面简化了数据处理过程，提高了实验结果的准确度；另一方面提高了我们大学生应用软件的能力，并且对实验数据处理方法上的一种创新，激发了我们大学生对学习大学物理实验的兴趣，提高了教学质量。 1实验原理[1] 在一块平面玻璃上安放上一焦距很大的平凸透镜，使其凸面与平面相接触，在接触点附近就形成一层空气膜。当用一平行的准单色光垂直照射时，在空气膜上表面反射的光束和下表面反射的光束在膜上表面相遇相干，形成以接触点为圆心的明暗相间的环状干涉图样，称为牛顿环，其如图1。 由干涉条件知，当时，干涉条纹为暗条纹。于是： 可改为 用Matlab软件作图并处理实验数据，作出～的关系图。 令 ； ； ， 可得到拟合直线方程 ，并从其斜率ｋ中求出平凸透镜曲率半径的值。 2实验数据记录 已知钠光的波长，仪器误差　 表1透镜曲率半径测量数据 环数 环数 K K 55 32.394 21.538 30 31.036 22.918 50 32.112 21.786 25 30.702 23.338 45 31.902 22.062 20 30.318 23.602 40 31.638 22.314 15 29.918 24.026 35 31.338 22.600 10 29.410 24.506 3 Matlab软件处理实验数据[2] 3.1得到与的数据表 令，，，用Matlab软件编写程序计算： x=[55 50 45 40 35 30 25 20 15 10] c1=[0.032394 0.032112 0.031902 0.031638 0.031338 0.031036 0.030702 0.030318 0.029918 0.02941] % XL的值 c2=[0.021538 0.021786 0.022062 0.022314 0.0226 0.022918 0.023338 0.023602 0.024026 0.024506] % XR的值 c3=c1-c2 %DK=（XL- XR）的值 y=c3.^2 % y-DK2的值 表2由与之间关系得到与的数据 10 15 20 25 30 35 40 45 50 55 24 34.7 45.1 54.2 65.9 76.4 86.9 96.8 107 118 3.2用polyfit功能函数作与的拟合图 3.2.1用Matlab软件最小二乘法公式求出斜率k、相关系数 用最小二乘法求出曲线的斜率k、相关系数， ， Matlab程序： a=polyfit(x,y,1) 结果得： k=0.2079e-5,b=0.3307e-5 3.2.2用plot函数绘与的拟合图2 Matlab程序： x1=10:5:55;k=2, y1=a(1)*x1+a(2); plot(x,y,’*’,x1,y1,’k’),xlabel(‘P’),ylabel(‘I’),grid on, Matlab拟合图见图2 所以拟合直线方程为: 4 Matlb软件最小二乘法误差分析计算 4.1 Matlab程序： yn=a(1)*x+a(2),xn=(y-a(2))/a(1), sigmayn=sqrt((yn-y)*(