拟牛顿法matlab程序_牛顿环实验的数据处理改进及图像分析

牛顿环是光的一种薄膜干涉图样&＃xff0c;最早是由牛顿于1675年首先观察得到的&＃xff0c;当用单色光照射牛顿环装置时&＃xff0c;在视觉上为一系列明暗相间的同心圆环。这些圆圈的距离不同&＃xff0c;随着到中心点的距离增加而逐渐变窄&＃xff0c;它们是由球面上与平面上反射的光线相互干涉而形成的干涉条纹。以往实验室通常使用的逐差法处理实验数据&＃xff0c;过程较为繁杂&＃xff0c;其计算结果也不够准确。在牛顿环实验中获得的某组数据通过逐差法计算牛顿环透镜曲率&＃xff0c;其相对扩展不确定度为0.82%&＃xff0c;由于“逐差”先两两相减使得自由度减半&＃xff0c;即重复测量或组合测量时计算实验标准偏差所用的独立残差个数减半&＃xff0c;未能充分利用测得的有效数据信息&＃xff0c;这不符合实验设计随机化原则的思路。本文通过线性拟合法处理实验数据&＃xff0c;并利用GUI建立一个可进行数值分析以及图像动态分析的交互式界面&＃xff0c;不仅可以得到更精确的实验结果&＃xff0c;其相对扩展不确定度为0.65%&＃xff0c;还能结合实际选择牛顿环实验参量以观察到不同的实验现象&＃xff0c;简便易行&＃xff0c;应用范围广。

1 牛顿环实验理论

牛顿环装置如图1所示&＃xff0c;牛顿环干涉的明暗条纹条件满足&＃xff1a;

得到牛顿环暗纹公式[1]

由于单色光的波长λ已知&＃xff0c;测量得到第k级的暗环半径&＃xff0c;可以计算得到牛顿环凸透镜的曲率半径。由于中心切点处的压力形变或残留粉尘作用&＃xff0c;实际方程式由上式修正为

(1)

b0指的是由于凸透镜中心切点处的压力形变或残留粉尘作用&＃xff0c;使得牛顿环的暗环公式不再是这样的正比例关系&＃xff0c;而是需要在右边加上一个常数&＃xff0c;这个常数反映的是压力形变或残留粉尘作用所带来的影响&＃xff0c;即b0。

图1 牛顿环装置示意图e指的是牛顿环装置中的凸透镜和底面平板之间的厚度距离

2 传统实验数据处理方法——逐差法

考虑到凸透镜和平面镜表面的尘埃和凹陷&＃xff0c;以及波长误差、平晶不平度、部分余弦(投影)误差、被测面偏离理想球面等误差使得显微镜观察下的牛顿环中心将会是一个圆斑。于是改为测量牛顿环的各暗环直径即可避免这种误差。若求得了第m级及第n级的暗环直径Dm,Dn&＃xff0c;记m-n&＃61;25即可得到计算公式&＃xff1a;

(2)

此时已知单色光源为钠光灯&＃xff0c;其波长λ&＃61;589.3nm&＃xff0c;通过逐差法即可计算得到曲率半径R。具体的逐差法计算结果见表2所示。(其中XL,XR分别表示某暗环中心在牛顿环左侧和右侧的位置)

表1 牛顿环各级暗环半径位置记录表

表2 逐差法处理数据结果

将数据代入式(2)&＃xff0c;计算透镜曲率半径值为&＃xff1a;R&＃61;875.34mm。

测量结果的B类不确定度主要由两个分量组成&＃xff0c;即仪器误差限和瞄准误差限。因变量(第k级暗环直径平方值)的扩展不确定度[2]

(3)

其中UDka即瞄准误差限&＃xff0c;通常是直径测量不确定度的主要因素之一[3]。

先考虑到仪器误差限带来的B类不确定度分量&＃xff0c;实验使用的是含有滚珠丝杆和1级螺纹副的测量显微镜&＃xff0c;其误差限ΔINS≈0.05,由B类仪器误差限分量的相对不确定度公式得到B类相对不确定度第一项分量&＃xff1a;

再考虑到瞄准误差限带来的B类不确定度分量&＃xff0c;由B类瞄准误差限分量的相对不确定度公式得到B类相对不确定度第二项分量&＃xff1a;

下面计算逐差法的B类方和根&＃xff0c;使用扩展不确定度的计算公式[2]&＃xff1a;

其中t为自由度ν&＃61;n-1&＃61;4&＃xff0c;在置信水平0.95时的学生氏t分布的值&＃xff0c;满足下列公式

(4)

将ν&＃61;4代入式(4)得到该t因子为&＃xff1a;t&＃61;2.7764。b1为直径平方逐差值(表2中一列5个实测数据)的平均值为直径平方逐差值与m-n比值的实验标准差&＃xff1a;0.00576。Ub1为逐差法B类不确定度分量&＃xff1a;

计算B类方和根为

B类相对不确定度合成为

此外&＃xff0c;符合1-α&＃61;0.95的A类相对扩展不确定度为tsb1/b1&＃61;0.8%

3 实验数据处理的改进技术和方法

考虑到“逐差”先两两相减使得自由度减半&＃xff0c;因而未能充分利用测得的有效数据信息。并且其操作计算步骤不具有简便性。由式(1)知&＃xff1a;

将离散的牛顿环曲率半径平方值与暗环级数画在一个坐标图中&＃xff0c;如图2所示。

图2 牛顿环暗环半径平方与暗环级数散点图

由图2知二者近似地成一次关系。对这些离散的点进行线性拟合所得曲线的斜率代表凸透镜曲率半径R和单色入射光的波长λ之乘积。因此牛顿环凸透镜曲率半径R亦可通过线性拟合法计算得出。

4 线性拟合处理实验数据

根据上述计算原理计算得到最后的拟合结果&＃xff0c;通过自制的GUI界面直观方便地呈现&＃xff0c;如图3所示。

图3 牛顿环曲率半径计算GUI界面

计算透镜曲率半径值为&＃xff1a;R&＃61;873.15mm。

针对瞄准误差限带来的B类不确定度分量&＃xff0c;由线性拟合斜率扩展不确定度的计算公式[2]&＃xff1a;

选取作为线性拟合因变量yi&＃xff0c;各环数(55,50&＃xff0c;…10)作为自变量xi。实验共拟合了10组数据&＃xff0c;故令线性拟合组数n&＃61;10;式(5)中t为自由度ν&＃61;n-2&＃61;8&＃xff0c;在置信水平0.95时的学生氏t分布的值&＃xff0c;同样满足式(4)。将ν&＃61;8代入得到该t因子为&＃xff1a;t&＃61;2.3060。sy为因变量的标准差。满足下式&＃xff1a;

sb1为拟合直线斜率b1的标准差。满足下式&＃xff1a;

xi为各环数为环数均值并将式(3)代入&＃xff0c;计算B类扩展不确定度为&＃xff1a;

考虑到仪器误差限对应的平均直径的B类标准差为显著小于因变量标准差sy&＃xff0c;可以忽略不计之[2]。此外&＃xff0c;符合1-α&＃61;0.95的A类相对扩展不确定度为tsb1/b1&＃61;0.63%

5 实验结果分析对比

对两种处理得到结果的相对不确定度进行分析比较&＃xff0c;由相对不确定度的计算公式&＃xff1a;

(5)

将和所分别代表的逐差法和线性拟合法对应的扩展不确定度以及量值斜率b1代入式(5)得到表3。由表3知&＃xff0c;通过两种方法计算相对扩展不确定度&＃xff0c;线性拟合法明显小于逐差法。由于逐差法两两相减&＃xff0c;自由度减半&＃xff0c;未能充分利用测量所得数据的有效信息。根2分布“总体标准差的散布宽度与自由度成反单调关系” 的规律&＃xff0c;自由度减半原则上将使结果散布范围增大&＃xff0c;显著减小测量准确度[2]。而线性拟合法在增加自由度的基础上其相对扩展不确定度更小&＃xff0c;可信度更高。因此线性拟合法的结果更为精确。

表3 逐差法和线性拟合处理数据相对扩展不确定度的对比表

6 牛顿环图像动态分析模拟

除了上述数值分析的模拟之外&＃xff0c;对牛顿环图像的分析也可通过该GUI界面进行模拟。利用计算机进行牛顿环图像的仿真处理[4]&＃xff0c; 仅考虑通过对干涉光强的模拟而绘制图像&＃xff0c;并未结合实际情况中入射单色光波长与凸透镜曲率半径的取值&＃xff0c;造成图像分析较为片面。考虑到实际选取的入射光波长有一定的取值范围&＃xff0c;该界面提供两种波长的选取方式。可以模拟实验室进行牛顿环实验的常用单色光&＃xff0c;如钠灯、汞灯&＃xff0c;通过弹出式菜单即可选取&＃xff1b;亦可通过滑动条任意选取位于可见光范围内的波长。移动各滑动条就可以改变凸透镜的曲率半径取值。

牛顿环图像的模拟选取波长为420nm的单色光照射&＃xff0c;凸透镜的曲率半径为520mm&＃xff0c;得到牛顿环图像如图4所示。实测暗环第5环、第15环、第25环的直径依次是&＃xff1a;2.087mm、3.621mm、4.669mm。

图4 初始牛顿环图像

图5 改变曲率半径参量后的牛顿环图像

控制波长不变&＃xff0c;增大凸透镜的曲率半径至1325.5mm&＃xff0c;得到牛顿环图像如图5所示。实测暗环第5环、第15环、第25环的直径依次是&＃xff1a;3.335mm、5.778mm、7.457mm。

控制曲率半径不变&＃xff0c;增加单色光的波长至739.6nm&＃xff0c;得到的牛顿环图像如图6所示。实测暗环第5环、第15环、第25环的直径依次是&＃xff1a;2.772mm、4.804mm、6.196mm。

图6 改变单色光波长参量后的牛顿环图像

可见由GUI界面所呈现的牛顿环图像得到得以下结论均与实验结果相符&＃xff1a;

(1) 若控制单色光波长参量不变&＃xff0c;牛顿环各环半径将随着凸透镜曲率半径的增大而增大。

(2) 若控制凸透镜曲率半径参量不变&＃xff0c;牛顿环各环半径将随着单色光波长的增大而增大。

(3) 实验室单色入射光的波长以及牛顿环装置的凸透镜曲率半径这两个参量的取值不可过小&＃xff0c;否则产生的牛顿环的各环半径太小&＃xff0c;由光强公可知&＃xff0c;光强在空间中的明暗分布过于紧密&＃xff0c;且牛顿环具有内疏外密的特点&＃xff0c;在显微镜的有效放大率内可能会造成牛顿环混杂&＃xff0c;使得实验观察困难&＃xff0c;观察者也会产生用眼疲劳。

7 结语

采用线性拟合法从相对不确定度的相关理论处理牛顿环的实验数据&＃xff0c;得到结果的相对扩展不确定度为0.65%&＃xff0c;更精确且符合实验设计随机化原则的思路。并且采用计算机程序进行实验的数值与图像模拟&＃xff0c;将实验室中繁杂的测量、计算、数据处理过程通过简单的GUI界面表示&＃xff0c;排除实验过程中遇到的包括透镜与平面镜之间的压力形变以及调整实验参量进行动态分析较困难的问题。通过计算机模拟与GUI界面的结合&＃xff0c;提高实验的准确程度&＃xff0c;便于总结牛顿环的实验规律&＃xff0c;可以得到更广泛的应用。

致谢: 本文在实验数据的不确定度分析方面得到了清华大学物理实验中心朱鹤年教授的悉心指导&＃xff0c;在此表示对朱教授的真诚感谢。

本文在选题、研究和修改过程中得到了西南交通大学物理实验中心姜向东教授的支持和帮助&＃xff0c;在此表示对姜老师的真诚感谢。

参考文献

[1] 张晓&＃xff0c;王莉.大学物理学[M].2版.北京:高等教育出版社,2014.

[2] 朱鹤年.新概念物理实验测量引论[M].北京:高等教育出版社&＃xff0c;2007.

[3] 考洛米佐夫.干涉仪的理论基础及应用[M].北京:技术标准出版社&＃xff0c;1982.

[4] 牛顿环干涉实验的MATLAB模拟[EB/OL].[2018-07-01].https://wenku.baidu.com/view/da6c174e5bcfa1c7aa 00b52acfc789eb172d9e05.html.

[5] 叶玉堂&＃xff0c;肖峻等.光学教程[M].北京:清华大学出版社,2011.

作者简介: 姜向东&＃xff0c;西南交通大学教授&＃xff0c;研究方向为光电传感器、物理电子学&＃xff0c;xdjiang&＃64;home.swjtu.edu.cn。

引文格式: 朱晓梅&＃xff0c;向伟铭&＃xff0c;姜向东. 牛顿环实验的数据处理改进及图像分析[J]. 物理与工程&＃xff0c;2019&＃xff0c;29(3)&＃xff1a;55-58,62.

END

更多精彩文章请点击下面“蓝字”标题查看&＃xff1a;

• 关于举办首届“全国慕课教育创新大会暨高校在线开放课程联盟联席会年会”的通知

• 关于发展中国大学物理教育MOOC联盟成员的通知

• 黄昆&＃xff1a;永远怀念吴大猷老师

• 黄昆&＃xff1a;我的研究生涯

• 葛惟昆&＃xff1a;香港科技大学授予荣誉理学博士的颂词

• 朱邦芬院士&＃xff1a;一代宗师黄昆

• 朱邦芬院士&＃xff1a;纪念中国半导体物理及固体物理奠基人黄昆先生

• 朱邦芬院士&＃xff1a;“减负”误区及我国科学教育面临的挑战 

• 朱邦芬院士&＃xff1a;从基础科学班到清华学堂物理班

• 朱邦芬院士&＃xff1a;对培养一流拔尖创新人才的思考

• 王炳燊、葛惟昆&＃xff1a;玻恩-黄昆方程的历史回顾——纪念黄昆诞辰100周年

• 贾惟义、葛惟昆&＃xff1a;在黄昆先生指导下研究固体中的晶格振动、光色散及电声子耦合——纪念黄昆先生诞辰100周年

• 李学潜教授&＃xff1a;物理是一种文化

• 李学潜教授&＃xff1a;如何帮助物理系学生迈过从高三到大一这个坎

• 李焱、杨宏&＃xff1a;从超短光到超短超强光的突破——解读获2018年诺贝尔物理学奖的啁啾脉冲放大技术

• 王雯宇&＃xff1a;光和引力波专题 Ⅰ ——广义相对性原理、光速不变原理及引力论

• 王雯宇&＃xff1a;光和引力波专题Ⅱ——电磁及引力介质理论

• 穆良柱&＃xff1a;什么是物理及物理文化&＃xff1f;

• 罗洪刚等&＃xff1a;“课程群”组织体系: 对物理学基础课程教学体系的思考

• 邓崇林&＃xff1a;创新量纲分析重导玻尔模型能级公式

• 物理教育的春天来了&＃xff01;

• 吴国祯教授&＃xff1a;我的国外研究生经历印象——应清华大学物理系“基科班20年·学堂班10年纪念活动”而写

• 《物理与工程》2019年第1期目录

• 陈佳洱&＃xff0c;赵凯华&＃xff0c;王殖东&＃xff1a;面向21世纪&＃xff0c;急待重建我国的工科物理教育

• 王亚愚教授&＃xff1a;清华物理系本科人才培养理念与实践

• 葛惟昆教授&＃xff1a;关于中外人才培养的几点思考

• 王青教授&＃xff1a;小班教学与翻转课堂&＃xff1a;《费曼物理学Ⅱ》的10年教学实践——纪念费曼先生百年诞辰

• 安宇教授&＃xff1a;为什么传统的课堂讲授模式需要改变

• 刘玉鑫教授&＃xff1a;关于本科生物理基础课程教学和教材编著的一些思考

• 沈乾若&＃xff1a;重创理科教育的美加课程改革

• Henderson C&＃xff1a;美国研究基金支持下的物理教育研究及其对高等物理教育的影响

《物理与工程》期刊是专注于物理教育教学研究的学术期刊&＃xff0c;是中国科技核心期刊&＃xff0c;1981年创刊&＃xff0c;欢迎踊跃投稿&＃xff0c;期刊投审稿采编平台&＃xff1a;

http://gkwl.cbpt.cnki.net

欢迎关注

《物理与工程》微信公众号