Matlab用于物理实验常见数据处理.pdfVIP

维普资讯 第 20卷 第 4期 大 学 物 理 实 硷 Vo1．2O No．4 2007年 J2月 出版 PHYSICALEXPERIMENrOFCOU正GE Dec． 2oo7 文章编号 ：1007—2934(2007)04—0085—02 Matlab用于物理实验常见数据处理 张世 功 (贵州大学，贵阳，550003；同济大学 ，上海 ，200092) 摘 要 本文利用 madab强大的计算功能对常用 的物理实验数据处理方法进行结合处 理，剔除异常数据，并对不同计算方法进行比较，将结果显示出来。 关键词 Matlab；数据处理 ；键盘 中图分类号 ：04—39 文献标识码 ：A 物理实验的数据处理方法一般有 ：逐差法 、作图法、最小二乘法等。在物理教学中，一 般将某种实验的数据用某种处理方法对 学生进行教学或深化。但各处理方法各有优缺 点：逐差法可能会产生有较大 的误差 ；作 图法有人的主观因素，误差也不好标定 ；而最小二 乘法虽然准确 ，但计算量偏大，错误数据会影响相关系数。这三种方法如果配合起来使用 可以发现究竟用那种方法 比较好，从而选用恰 当的实验数据处理方法。 1 逐差法 逐差法的原理已经非常清楚，这里不在赘述。但是实验数据的个数却对逐差法有影 响，是偶数对数据 (rl=2m)时，我们可以很容 易算 出来斜率 ；当是奇数对数据 (rl=2q一1) 时，需要将中间数据(n=q)数据对舍弃，比如 11对(rl=11，q=6，nl=n／2=5．5数据 ，将第 6对数据舍弃，这样可以减小误差 ：假如 X=1，2，3…I1，很简单， 与 相 比较 ，虽 然结果差不多，但是误差却明显不同，假如分子 的误差都是 0．1，分母没有误差 ，前者结果 的误差是 0．017，而后者是 0．02，相差 0．03，相对误差 17％!所以奇数对数据做逐差法时， 绝对不是随便舍弃最前或最后的数据 !当然实验数据的增多会减小这种误差。 2 最小二乘法 最小二乘法使拟合的线性方程与每一个数据点的均方根差最小 ，就充分地利用到了 每一个数据 ，但计算量太大 ，如果交给计算机去做 ，就省事多了。但是如果有一组数据的 偏差太大 ，或者说是一组错误的数据不能发现 ，也会影响拟合的线性方程。 我们先让计算机显示数据点．然后最小二乘法计算，显示拟合的直线，直观在图像 中 收稿 日期 ：2007—08—13 ～ 85 — 维普资讯 观察数据 ，确定数据是否有问题 ，决定修改(或删除)该数据点，再次进行计算 ，就得到较好 的拟合线性方程 。 Mdab里有一个命令叫做 keyb0ard，它可 以中断程序、激活键盘、等待命令，就可以修 改数据了，修改完后 ，利用 retum命令回到程序继续运行 。从而实现再次计算 。 键盘修改部分是在matlab的命令窗口进行的，注意函数窗13与命令窗13的切换。另 外输入字符 “y，，后将删除某组数据 ，数据对将 由n变为 n一1了!否则就出错了。这样新 数据(剔除了一组数据)将再次进行最小二乘法计算 ，利用hold。n命令将两次计算的图合 在同一个图片里 ，如图 1。 图1中，X：[1234567]，y=，明显 由于x(7)=7；y(7)=12与其他 点偏离太远，拟合的直线 (图1中的虚线)比较别扭，然后 keyboard命令里的x(7)：[]；y 7)：[]；删除了这个点，拟合的直线(图 1中的实线)就 比原先的好 了。而计算结果可以 在命令窗口看到。 圈 1 最小二乘法翻除掉错误点后再次拟合 作图法是人为地画条曲线，均匀地从数据点中穿过，使数据点均匀