弗兰克—赫兹实验的数据处理方法设计.doc

弗兰克赫兹实验的数据处理方法设计系别： 物理系专业： 物理学班级： 2009级B班学生姓名：刘立宏 学号：20091711233指导教师姓名：王蕴杰 职称：讲师 最后完成时间：2013年5月弗兰克赫兹实验的数据处理方法设计 刘立宏（青海师范大学 物理系 810008）摘要:本文利用计算机中的Matlab编程技术，并用插值的方法，绘制出了弗兰克赫兹实验数据的曲线,大大提高了实验结果的精度。关键词：弗兰克赫兹实验 插值 拟合 MatlabAbstract: In this paper, using the Matlab programming technology in the computer, and the method of interpolation, drawed the Hertz experiment data of the Frank curve greatly,improved the accuracy of experimental results.Keywords: Frank - Hertz experiment fitting interpolation Matlab0 引言 在计算机应用普及的今天，计算机技术广泛应用于科学计算，数据处理等方面,将计算机与有关的实验观测仪器相结合，可对实验数据进行记录、整理、加工、分析和绘制图表， 大大减小了繁杂的数据处理过程,使实验结果更直观,准确。本文通过对弗兰克赫兹管的实验数据进行分析处理，并在Matlab程序里利用插值法对数据进行分析处理，给出了弗兰克赫兹管实验时的最佳参数1 弗兰克赫兹实验1.1理论背景根据波尔理论，原子具有分立的能量,，又称为能级。正常状态下原子既不辐射也不吸收能量，当原子内电子受激发从低能级跃迁到高能级时，原子就处于受激发状态。原子在能及之间跃迁时，发射或吸收的能量满足普朗克公式 弗兰克赫兹实验便是利用电子与原子的碰撞，使其实现原子状态的跃迁，即当 时，即可实现跃迁。若原子吸收能量从基态跃迁到第一激发态，则称为第一激发电位。1.2 弗兰克赫兹实验实验原理如图1所示，在弗兰克赫兹管中，K是发射电子的热阴极；是控制阴极发射电流的第一栅极；为第二栅极；A是板极；在板极A和第二栅极之间加有拒斥电压;电子通过第一栅极进入区域，与管中的氩原子发生碰撞。当加速电压较小时，损失的能量很小，电子与原子发生弹性碰撞。当加速电压超过第一激发电位时，电子与原子发生非弹性碰撞，碰撞后电子将动能部分或全部传给被测气体，使其跃迁到第一激发态，而电子流则显著下降，由于原子内部的能量是量子化的，因此板极的电流会随加速的增加而有规律性地上下浮动变化。2弗兰克赫兹实验数据及处理2.1 matlab简介在Matlab中可以将数值分析，矩阵计算，科学数据可视化以及非线性动态系统的建模和仿真等诸多强大功能集成在一个易于使用的可视化环境当中。用户也可以将自己编写的实用程序导入Matlab函数库中来进行调用，本实验在对数据进行绘图处理时，就是用Matlab程序插值的方法将曲线绘制出来的。2.2插值的原理插值法是实用的数值方法，是函数逼近的重要方法。在生产和科学实验中，自变量x与因变量y的函数y = f(x)的关系式有时不能直接写出表达式的关系式，而只能得到函数在若干个点的函数值或导数值。当要求知道观测点之外的函数值时，需要估计函数值在该点的值。如何根据观测点的值，构造一个比较简单的函数y=(x)，使函数在观测点的值等于已知的数值或导数值的值。用简单函数y=(x)在点x处的值来估计未知函数y=f(x)在x点的值。寻找这样的函数(x)，办法很多。(x)可以是一个代数多项式，或是三角多项式，也可以是有理分式；(x)可以是任意光滑（任意阶导数连续续）的函数或是分段函数。函数的类不同，自然也有不同的逼近效果。2.3 弗兰克赫兹实验数据及处理2.3.1 弗兰克赫兹实验数据加速电压范围为14.5V-72V，测量时电压间隔取0.5V，对应的板级电流如下：0.2740 0.2900 0.3060 0.3190 0.3330 0.3480 0.3610 0.3710 0.3790 0.3820 0.3820 0.3790 0.3680 0.3620 0.3510 0.3520 0.3780 0.4330 0.5000 0.5670 0.6250 0.6800 0.7310 0.7710 0.8160 0.8470 0.8680 0.8830 0.8860 0.8750 0.8490 0.8060 0.7470 0.6720 0.5890 0.5020 0.4470 0.4780 0.5700 0.6770 0.7810 0.8990 0.9720 1.0660 1.1470 1.2260 1.2950 1.3470 1.3860 1.4100 1.4130 1.3930 1.3930 1.3410 1.2620 1.1470 0.9950 0.8120 0.6240 0.4990 0.5170 0.6290 0.7640 0.9100 1.0490 1.1920 1.2750 1.5070 1.6180 1.7050 1.7810 1.8350 1.8710 1.8780 1.8530 1.7930 1.6930 1.5500 1.3590 1.1270 0.8780 0.6950 0.6690 0.7620 0.8940 1.0430 1.1940 1.3530 1.5140 1.6590 1.8020 1.9350 2.0480 2.1450 2.2250 2.2810 2.3170 2.3050 2.2710 2.1980 2.0830 1.9190 1.7100 1.4680 1.2210 1.0440 1.0020 1.0570 1.1600 1.2960 1.4470 1.6070 1.7740 1.9370 2.0970 2.2430 2.3700 2.4890 2.5850 2.6590 2.7080 2.7220 2.7110 2.6620 2.5720 2.4440 2.27302.3.2处理程序及结果利用Matlab编制处理程序如下：clc;close all;x_max=max(u1); x_min=min(u1); %对实验数据进行拟合m=9; %对实验数据进行拟合时选用的阶数 x_max=max(u1); %确定函数图象x轴长度 x_min=min(u1); %确定函数图象x轴起始点 x=x_min:0.1:x_max; %取计算拟合曲线各点斜率的间隔 N=length(x)-1; y=interp1(u1,u2,x,spline);figure set(gca,FontSize,13); plot(u1,u2,*,x,z,r,linewidth,1.4) ,grid % 绘制曲线拟合后的图形 str3=fontsize13fontname黑体波谷序号N; str4=fontsize13fontname黑体加速电压U_G2K(单位：V); xlabel(str3) ylabel(str4) legend(fontsize13fontname黑体数据点,fontsize13fontname黑体线性拟合曲线,4) title(fontsize13fontname黑体弗兰克-赫兹谷值线性拟合曲线) hold on利用程序段处理后显示结果如图2所示。 图 2 弗兰克赫兹实验曲线3 结论 利用Matlab编程技术对弗兰克赫兹实验数据用最小二乘法进行了线性拟合，并绘制出了其实验曲线，并得到氩原子的第一激发电位为11.61 V，；详细的分析表明该