基于Matlab用最大公约数法求电子电量.pdf.pdfVIP

第 25卷 第 2期 大 学 物 理 实 验 Vo1．25No．2 2012年 4月 PHYSICAL EXPERIMENT OFCOLIEGE Apr．2012 文章编号：1007—2934(2012)02—0063—03 基于 Matlab用最大公约数法求 电子 电量 王玉梅 (陕西理工学院，湖北 汉中 723003) 摘 要：用Matlab程序辅助，在允许一定误差存在的情况下，给出了求正实数 的最大公约数的函 数，并把其运用于密立根油滴实验的数据处理中，求出了最小电荷电量。 关 键 词：密立根油滴实验；最大公约数法；Matlab；最小电荷电量 中图分类号：O411—2 文献标志码：A 密立根油滴实验是证 明电荷电量量子化的基 本实验 ，在近代物理学中具有极其重要的地位 。 I 基本算法 其构思巧妙、方法简便，目前仍是各高等院校理工 科学生必做的实验。从实验设计至今，实验仪器 1．1 求两个整数最大公约数的基本思想 不断更新、改进 ，测量也变得更加方便。但是在数 整数a、b组成数组p===[a，b]，用其中的大数 据处理问题上，始终没有找到一个很好的方法[1]。 max(p)减去小数 min(p)，得到差 cha，若 cha／ 李波欣L2介绍的方法处理结果固然好，但还是利 rain(p)小于一定的数值 1／max(p)，则 min(p) 用了电子电量的理论值在先 (事先规定所求结果 即为最大公约数；否则令 a—cha，b— rain(p)，重 在 1．550001～1．650001之间)，还是存在倒推法 复上述步骤，直至输 出最大公约数。用此思想编 意识。丁红星等 介绍的平均值正整数法依然不 写的程序名为 wgys(a，b)，用它可以得到任意两 能完全排除人为干扰，另外处理起来还是 比较繁 个正整数的最大公约数。 琐的。李娟_4]等介绍的是统计法处理数据，不适 1．2 求两个小数最大公约数的基本思想 合个人处理。宋五~,llE]介绍的方法倒是思路新 允许一定误差存在的情况下，求小数的最大 颖、巧妙，值得学习、借鉴。大家公认的获得基本 公约数，难在如何允许一定的误差存在，如何控制 电荷量的方法是最大公约数法，但是最大公约数 误差。笔者从数据量子化的特点开始考虑，又考 目前只能求整数的，求取小数的最大公约数还没 虑到实验观察到的油滴所带电量必为最小电荷 电 有现成的函数。王广涛_6等介绍的方法就有最大 量的简单整数倍 n／m(n、m均取值于 1～2O)，因 公约数的思想，只是不太容易理解 ，且其 中阈值的 此制定了下述思想来求最小电荷 电量。 取法不太清楚。本文在求整数最大公约数思想的 实数 a、b(用科学计数法表示，取 同一幂次， 基础上，本着合理处理密立根油滴数据的想法，用 a、b即为两个实数的去掉幂次的那一部分，a、b中 Matlab程序_7]辅助，在允许一定误差存在的情况 的最小值处于区间 l～10中)组成数组 P一[a， 下，给出了求正实数的最大公约数的函数，并把其 b]，用其中的小数 min(p)除以大数 max(p)，所 运用于密立根油滴实验的数据处理中，求出了最 得商 sh减去 n／m得到差cha，cha除以sh和n／m 小电荷 电量 。此程序不仅可 以大大简化运算过 两个数中的最大值得到判据 pan，先令 n一1，m一 程 ，而且其思想 内涵可 以让大家真正体会到电荷