1、重 庆 师 范 大 学物 理 与 电 子 工 程 学 院学 年 论 文专 业 ： 物理学（师范） 年 级 ： 2012级 姓 名 ： 尚俊霖 学 号 ： 20120511346 导 师 ： 毋志民 二 一 四 年 九 月 二 十 二 日用牛顿环测透镜曲率半径的数据处理方法尚 俊 霖 20120511346 物 理 学 （ 师 范 ） 2012级摘要：牛顿环实验是大学物理实验中非常重要的实验，以单色平行光投射到牛顿环装置上，则由空气膜上下表面反射的光波将互相干涉，形成的干涉条纹为博的等厚各点的轨迹，这种干涉是一种等厚干涉。处理该实验的测量数据常采用逐差法，最小二乘法，加权平均以及其它方法。通过介
2、绍用逐差法、线性回归法、加权平均法处理牛顿环测透镜曲率半径数据的方法和过程。比较三种实验数据处理方法的优缺点,其中加权平均法既考虑了如何克服实验的系统误差,又能按照处理原则去对待非等精度测量,且建立在数理统计理论基础上。该方法主要是比较相应的权,进而求出加权平均值,得出加权平均法为牛顿环实验数据处理的最佳方法，但加权平均法中要计算的数据较多，公式较多，较传统的方法要复杂的多。故探讨如何简化加权平均法，利用 Matlab 软件进行处理。关键字：牛顿环实验、加权平均法、非等精度实验数据的处理、干涉条纹 、加权平均值 引言：牛顿环是牛顿在 1675 年所做的著名实验。牛顿环是等厚干涉的一种，它在光学
3、计量、基本物理量测量等方面有广泛的应用：用牛顿环测定光波的波长、透镜曲率半径，用牛顿环来检验磨制透镜的质量等。处理牛顿环的测量数据常采用逐差法，最小二乘法，加权平均以及其它方法，因为该实验测量是非等精度的测量，逐差法可以很好的克服实验的系统误差，但是没有按照数据的处理原则去对待非等精度的测量，该方法虽然被广泛的使用，但是处理的结果并不是很理想，最小二乘法回避了非等精度性的困难，但是它没有考虑实验的系统误差，所以用该方法来处理实验得到的数据也不是很理想的。基于以上两中方法的利弊，本实验可以采用加权平均法来处理实验数据，它既考虑了如何消除实验的系统误差，又按照了数据处理原则来处理非等精度的测量，因
4、此它是处理牛顿环实验数据的较为理想的方法。1 实验原理牛顿环仪是由待测平凸透镜 L(曲率半径约为 200700 cm)和磨光的平玻璃板P 叠合装在金属框架 F 中构成(如图一)，框架上有三个螺旋 H，用以调节 L 和 P之间的接触，改变干涉环纹的形状和位置。调节 H 时，螺旋不可旋的过紧，以免接触压力过大引起透镜弹性形变，甚至损坏透镜。如图二所示将曲率半径很大的平凸透镜的凸面放在一光学平面玻璃上，在透镜和平面之间形成空气膜，以平行单色光垂直照射时，经空气膜层上，下两表面反射的两束光发生干涉，在空气膜上表面出现一组干涉条纹。干涉条纹是以接触点 O 为圆心的一系列同心圆环，称为牛顿环。在图中，设
5、r 为牛顿环某环的半径，e 为与该环对应的空气膜层的厚度。考虑到光在空气膜下表面反射的光，是从光疏介质(空气)入射到光密介质(玻璃)，有半波损失，而在空气膜上表面反射的光，是从光密介质入射到光疏介质，无半波损失。所以在空气膜上、下表面反射的两束反射光的光程差为(1)2e在直角三角形 AOC 中，有2)(rR从而得图二测量原理图一 牛顿环仪eRr2考虑到 eR, e 跟 R 相比可以略去，即(2)r2代入（1）式，可得到2Rr根据干涉相长和干涉相消的条件可得明环半径为(3)2)1(Rkr),21(k暗纹半径为(4)kr ),0(必须指出，由于干涉条纹有一定宽度，上式中的 r 是第 K 级牛顿环的
6、条纹中心到圆环中心的距离，根据二光相干时光强分布的理论计算可知，各级牛顿环的条纹宽度并不相同，K 愈小，即距离环中心愈近，条纹愈粗，中心为一圆形暗斑。将（4）式加以变换，可得(5)krR/2显然，只要测出第 K 级暗纹的半径 rk，由已知波长 即可根据上式算出曲率半径 R。但由于接触压力引起的弹性形变使接触部位不是一个点而是一个小圆面， 圆环中心为一暗斑，使得中心难以找准，这样，干涉级数 k 和第 k 级暗纹半径 rk都难以测准，另外，接触面镜面上可能有微小灰尘存在，会引起附加光程差，这会给测量带来系统误差。设 Dm、D n分别是第 m 级和第 n 级暗环的直径，由式(4)可得到：R42RDn
7、42(6)(/)(2nm(7)rr不难证明用上面的式子计算曲率半径 R 可消除前面所述因素的影响，而且, 即使中心 O 未能找准，测量的 D 不是直径而是牛顿环的弦长，也不产生原理性误差。本实验已知钠光波长 =5 8 9.3 nm，测出 Dm和 Dn代入(6)式即可求出 R。2 牛顿环干涉条纹的特点（1） 干涉图样是以接触点为圆心的一组明、暗相间的同心圆环，有半波损失时，中间为一暗斑。（2） 从中心向外，条纹级数越来越高，条纹的间隔越来越密。（3）用白光照射将形成彩色光谱，对每一级光谱，红色的在外圈，紫色的在内圈。（4） 增大透镜与平板玻璃间的距离，膜的等厚线向中心收缩，则干涉圆环也向中心收缩
8、（内陷） ，膜厚每改变 ，条纹就向外冒出（扩张）或向中心内陷一条。2/n3 实验数据处理：定义 为 个相邻牛顿换直径平方差的测量精度，由误差传22mkD递推倒可以得到： )(02. 222 mkmk D由此可见，当 取不同的值的时候， 也不同，所以该实验是非等k k精度的测量，应该用加权平均法处理该实验的数据，不妨令，由于本实验的 ，且 相应的权重2mkkDX 20,20kkX所 以 kX为，112202kkkDa所以加权平均值为： ，kkaX)(很明显， 都是非等精度的测量值，其对应的权值分别为：,.,21nX，所以 的均方差误差为：na.,21 k，kkkanXXS)1(2由此可见， ，m
9、DRkk4)(2XSk4最终的结果为： R3.1 加权平均法: R3=R3R3=(1 271.10.2) mm 下面对这三种数据处理方法进行检验,选择最优的数据处理方法,检验方法较多,现选择采用 t 分布检验9: t=x1-x2(n1-1)2 1+(n2-1)22(1/n1+1/n2) 式中:n1 和 n2 分别为凸透镜球面的上、下两面的折射率,由于凸透镜球面周围都为空气薄膜,故 n1=n2,则令 =n1+n2-2=2(n-1), 从而有: t=（-R2)/(2R1+2R2)方法 1 与方法 2 比较计算,可得:t1=0.350;方法 2 与方法 3 比较计算,可得:t2=0.340。 若取显
10、著水平 =10%,则置信率 p=90%,=18,查 t 分布表可得10t=1.734, 则|t1|=0.3541.734,|t2|=0.3401.734。 若取 = 50 %,则 p=50%,=18,查表得 t=0.688,则|t1|=0.3540.688,|t2|=0.3400.688 。 通过上面分析可以看出三种数据处理方法有如下特点: (1) 逐差法主要是围绕如何克服实验的系统误差来进行的,是建立在算术计算的基础上,但并不满足非等精度测量实验数据处理的条件,而牛顿环干涉实验是非等精度测量,故逐差法对于牛顿环实验来说并不是一种理想的数据处理方法。(2) 线性回归法主要是为了避免非等精度测量
11、的困难,但未考虑该次实验中的系统误差,所以线性回归法对于牛顿环实验来说也不是理想的数据处理方法。 (3) 加权平均法既考虑了如何克服实验的系统误差,又能按照处理原则去对待非等精度测量,且建立在数理统计理论基础上,所以加权平均法是处理牛顿环实验数据的最佳方法。3.2 实验数据处理 程序MATLBfunction DP(D)L=D(1,:);n=numel(L);c=0.0002;X=0;m=20,k=2,h=5893*power(10,-7);S3=0;S4=0;for i=1:n;a(i)=D(1,i)-D(2,i);b(i)=D(3,i)-D(4,i);x(i)=a(i)2-b(i)2;y(
12、i)=a(i)2+b(i)2;p(i)=1/(c*y(i); endfor i=1:n;X=X+(p(i)*x(i)/sum(p(:);endfor i=i:n;S3=S3+p(i)*(x(i)-X)2;S4=S4+p(i);endS1=(S3/(S4*(n-1)0.5;R=X/(4*m*h*1000),S=k*S1/(4*m*h*1000),3.3 数据处理举例环的级数（k） 45 44 43 42 41 40左环的位置/mm 38.353 38.271 38.176 38.075 37.994 37.901右环的位置/mm 22.206 22.383 22.531 22.644 22.73
13、2 22.834环的级数（k-m） 25 24 23 22 21 20左环的位置/mm 36.321 36.256 36.142 36.014 35.885 35.785右环的位置/mm 24.164 24.406 24.529 24.656 24.781 24.912表 1：用牛顿环测量透镜的曲率半径实验数据应用上述 程序处理该实验数据得到的结果为：MATLBm = 20k = 2R = 2.3502S = 0.0381其中 为该实验所采用的级差， 为不确定度的扩展系数， 为侧得的透镜kR的曲率半径的平均值， 为不确定度。S所以该透镜的曲率半径为：)2()(038.5.2 kmSR ：扩 展
14、 系 数 为单 位 ：3.4 程序的使用效果该程序使用起来比较简单，做完实验后只需要将实验数据按表 1 的形式填入到表格中，然后将表格中的数据直接导入到 中，运行程序就可MATLB以了，方便且易行。结束语：本文对牛顿环实验数据加权平均法进行分析。逐差法在牛顿环干涉实验中是一种常用的实验处理方法,其原理简单且便于理解,对它的实验原理不用再做过多的叙述,但由于逐差法不满足非等精度测量实验数据的条件,而牛顿环干涉实验就是一种非等精度测量,故该方法对于牛顿环干涉实验并不是一种理想的实验处理方法;线性回归法先利用数值插值法对实验数据进行处理,再利用最小二乘法将实验数据拟合成一条直线函数,最后用Matlab软件计算出线性拟合系数B及相关系数r,进而算出凸透镜的曲率半径R和测量的相对不确定度;加权平均值法主要是比较相应的权,进而求出加权平均值,利用Matlab软件处理较为方便,在优化模型中应用较广。经过分析与讨论可知应用加权平均值法为牛顿环实验数据处理的最佳方法。参考文献：大学物理实验 姚安居 吴庆州 编 中国矿业大学出版社光学新概念物理教程赵凯华主编 高等教育出版社 MATLAB基 础 应 用 简 明 教 程 张 平 主 编 北京航空航天大学出版社
