弦振动规律及MATLAB数据处理

1、弦振动形成驻波的规律和数据的MATLAB处理苗锟1，黄育红2, 本论文受陕西师范大学“大学生开放性实验基金项目”资助。 作者简介：黄育红（1980-），陕西宝鸡人，陕西师范大学物理学院讲师，从事物理实验教学；现为西北工业大学在读博士，主要从事半导体物理理论计算。，李康1，宋翠1（1陕西师范大学理工科基础教学部，陕西 西安，700062；2陕西师范大学物理学与信息技术学院，陕西 西安，710062）摘 要 本文通过对大量实验数据进行分析，同时采用MATLAB软件对数据进行系统处理，作出了数据间的关系图。借助软件使物理分析简单化，形象化，进而使实验者有更多的精力关注其内在的物理思想。并在此基础上，

2、分析总结出弦振动形成驻波时波长与张力以及波长与振动频率的规律，提出检验判断实验数据的方法，可以使实验者在做实验的过程中随时判断数据的正确与否。关键词 弦振动；驻波；MATLAB中图分类号：O4-34 文献标识码：A1 引言“弦振动形成驻波”一直是高等学校普通物理实验中的传统力学实验之一，它是研究波的形成、传播和干涉的一个重要途径。在固定均匀弦振动实验中，利用XZDYB型均匀弦振动仪，通有交流电的金属弦线在磁场中受到安培力的作用能够振动，在满足一定的条件下可形成驻波，通过观察驻波波形，测量相关的物理量，可测量弦线的线密度和横波传播的波长和波速，进而可以了解驻波传播的规律。有关学者在这方面做了不少

3、讨论。以前都是电动音叉受振动产生驻波的实验装置，现在改用均匀固定弦振动仪，使实验操作简化，物理现象明显，但是仪器仍然有些不尽人意的地方。另外这个实验数据多而繁杂，很容易引起混淆，且计算量大，在实验数据处理方面，到目前为止，还没有人提及对实验数据进行实时判断的具体方法，大量的实验数据极易引起混淆。为了更清晰明确整个实验的体系和思路，作者系统研究弦振动产生驻波的各种情况，在取得大量实验数据的同时，用MATLAB软件处理实验数据并作图，使多而杂的实验数据所内含的物理规律清晰而形象化，并在此基础上，提出一种能实时检验判断数据的方法，使后来的实验者能随时判断自己测量数据的正确性，也为下一步的实验测量奠定

4、了基础。 2 弦振动形成驻波的原理驻波是两个振幅相同、频率相同的相干波在同一直线上沿相反方向传播的叠加，固定均匀弦振动在介质中的传播满足一定的条件就可形成驻波，如图1所示。弦振动是最直观、最简单的驻波模型。正向传播的波为：反向传播的波为：式中x为质点位置坐标；t为时间；A为振幅；f为频率；w=2pf称为圆频率；l为波长；k=2p/l称为波矢； n=lf为波的传播速度。两列波叠加的结果，任一点x的合成振动为： 由上式可见，入射波与反射波合成后，弦上各点都以同一频率作简谐振动，它们的振幅为，即驻波的振幅与时间t无关，与质点的位置x有关。ABxxxt=0t=T/4t=T/2000/2波腹波节图1 相

5、反方向传播的两列波形成驻波的示意图波节处振幅为零，即通过简单计算可求得相邻两波节之间的距离为：波腹处的质点振幅为最大，可计算出相邻两波腹间的距离也是半个波长。因此，在驻波实验中，只要测得相邻两波节（或相邻两波腹）间的距离，就能确定该波的波长。由于固定弦的两端是用劈尖支住的，故两端处的点不会动，必为波节，所以，只有当弦线的两个固定端之间的距离（弦长）等于半波长的整数倍时，才能形成驻波，这就是均匀弦振动产生驻波的条件，其数学表示式为： 由此可得沿弦线传播的横波波长为： (1)式中的为弦线上驻波的段数。波动理论指出，弦线中横波的传播速度为 (2) 式中为弦线中的张力，为单位长度弦线的质量(即弦线的线

6、密度)。根据波速、频率及波长的普遍关系式，可得 (3) 由(2)式和(3)式可得 (4)由此可见，当给定、时，频率只有满足上式的要求才能在弦上形成驻波。同理，当用外力（例如流过金属弦线上的交变电流在磁场中受到交变安培力的作用）去驱动弦线振动时，外力的频率必须与这个频率一致，才会促使弦振动的传播形成驻波。实验数据记录及MATLAB处理调整好实验装置，实验中分别测量频率一定和张力一定条件下相关的物理量，并利用MATLAB 处理实验数据，绘制各物理量之间的曲线图，使它们之间的关系更加明确和清晰化。（西安市的重力加速度值为）3.1 在频率一定的条件下，改变张力的大小，测量弦上横波的传播速度。选取频率，

7、张力T由砝码挂在弦线的一端产生（需考虑托盘质量），以20.0g砝码为起点，逐次增加5.0g直至140.0g为止。同时调节相应的弦长，使弦线上出现、两个驻波段。记录相应的、值，如表1所示，做测量曲线图，如图1所示。由(3)式计算弦上横波速度的测量值。表1 频率一定时实验数据测量表砝码质量张力/N2029.72044.3064.8020.5025.8064.5038.7039.852534.72044.4064.2019.8036.0079.4043.4041.503039.72037.4060.8023.4026.3571.8045.4546.133544.72035.5060.1024.602

8、0.0269.4149.3949.304049.72039.6065.8826.2824.8276.5551.7352.154554.72040.9068.0027.1030.1084.3854.2854.245059.72039.6568.0528.4022.4079.3856.9856.895564.72039.3269.1029.7820.1479.2059.0659.316069.72039.8070.6530.8526.9088.6061.7061.706574.72036.7868.6931.917.5071.2363.7363.787079.72034.0067.3033.301

9、0.0075.8065.8066.207584.72034.7468.8234.0810.0078.1568.1568.168089.72034.5869.5935.018.1077.9069.8069.918594.72036.3072.1135.819.0081.0872.0871.859099.72035.0571.6036.559.0083.1074.1073.6095104.72036.0073.9537.955.4581.7576.3076.10100109.72031.8570.6838.835.0583.0077.9577.81105114.72032.4271.5539.13

10、3.8683.2079.3478.80110119.72034.4774.6240.153.7084.8981.1980.75115124.72023.8065.1041.305.5188.7283.2182.91120129.72027.0069.3042.303.2888.4085.1284.86125134.72029.0071.3042.302.5088.5086.0085.30130139.72029.6372.0642.433.5591.3387.7886.32135144.72034.1077.4543.353.5290.2086.6886.69140149.72028.4972

11、.5244.032.1590.3988.2488.15图 1 测量曲线（注：其中包括砝码的质量和托盘的质量）程序：Datax=29.72，34.72，39.72，44.72，49.72，54.72，59.72，64.72，69.72，74.72，79.72，84.72，89.72，94.72，99.72，104.72，109.72，114.72，119.72，124.72，129.72，134.72，139.72，144.72，149.72; Datay=39.85，41.50，46.13，49.30，52.15，54.24，56.89，59.31，61.70，63.78，66.20，68.1

12、6，69.91，71.85，73.60，76.10，77.81，78.80，80.75，82.91，84.86，85.30，86.32，86.69，88.15;plot (Datax(1:25),Datay(1:25),k*-),xlabel(张力T/N), ylabel(平均波长/cm)3.2 在张力一定的条件下，改变频率，测量弦线上横波的传播速度.将40.0 g砝码挂在弦线的一端，频率以60为起点，逐次增加5直至175为止，调节相应的弦长，仍使弦长出现、两个驻波段。记录相应的、值，如表2所示，做测量曲线图，如图2所示。由(3)式计算出弦上横波传播的速度.表2 张力一定时实验数据测量表603

13、3.4076.2042.804.4594.3089.8587.736531.7271.9040.185.5088.0182.5181.447032.9269.2936.375.0080.9075.9074.327530.0064.5034.505.0076.2971.2970.158044.0076.6032.606.0072.6866.6865.948540.0073.5033.506.0068.5662.5664.789040.0068.7828.7810.0068.2058.2057.889545.0072.1027.1010.0066.5056.5055.3510037.0062.90

14、25.9020.0071.9051.9051.8510540.0064.4824.4825.0074.0049.0048.9811041.0064.7023.7025.0071.7046.7047.0511540.0062.5022.5030.0074.7044.7044.8512040.0061.5021.5022.2065.2743.0743.0412539.8060.5020.7022.1164.1242.0141.7113038.6058.4119.8123.3863.6140.2339.9313539.3158.6219.3125.0963.8638.7738.7014040.285

15、8.9118.6327.2564.5337.2837.2714539.0056.4917.4927.1963.1735.9835.4815039.3857.0217.6427.3461.4034.0634.6715539.5056.0116.5129.2862.3833.1033.0616040.2656.1015.8431.8963.9032.0131.8516546.0061.8415.8433.3565.1431.7931.7417044.5159.7815.2734.3865.0330.6530.6017544.1058.8814.7832.8562.5929.7429.65图 2 测

16、量曲线程序：Datax=60，65，70，75，80，85，90，95，100，105，110，115，120，125，130，135，140，145，150，155，160，165，170，175; Datay=87.73，81.44，74.32，70.15，65.94，64.78，57.88，55.35，51.85，48.98，47.05，44.85，43.04，41.71，39.93，38.7，37.27，35.48，34.67，33.06，31.85，31.74，30.6，29.65;plot (Datax(1:24),Datay(1:24),k*-), axis(40,180,20,

17、90),xlabel(频率/Hz), ylabel(平均波长/cm)弦振动形成驻波的规律由(4)式 分析可知： 因铜丝线密度是一定值，也可以设计实验测量其大小。当频率固定不变时，张力越大波长越大，这也可由表1的计算数据和图1看出。此时，对表1中的数据横向比较，即在某一特定张力下，理论上也就一定了，两段驻波对应弦长之比：=2：1，推而广之，若实验条件允许，能测出段驻波，则各段驻波对应弦长之比为：2：1。 当张力固定不变时，频率越大波长越小，即某一段驻波下，频率越大，测得对应弦长越小。由表2和图2可以很直观的看到这个规律。对表2中的数据横向比较，即在某一固定频率下，理论上两段驻波对应弦长之比应为：

18、=2：1，当然在实验过程中，由于各方面的误差，弦长之比近似为2：1，这也可以侧面检验我们实验数据的可靠性。在学生掌握以上规律后，实验进行过程中若出现明显不符合此规律的数据时，即可很快判断引起该粗大误差的原因，分析导致此现象的深层次原因。该方法能有效避免因操作不当而记录的错误数据，通过提炼总结理论中的规律，也培养了学生理论联系实际，站在理论高度分析问题的能力。实践证明此法行之有效。检验判断实验数据的方法判断实验数据的有效性可通过以下几种方式进行：(1) 列表法。使用数据表检验数据，表中未列出的可根据内插法推出。(2) 作图法。通过简单运算，检验在一定的张力或频率下，所得的波长与图中对应的数据的差

19、距，判断所测数据是否合理。(3) 根据总结出的弦振动形成驻波的规律。实际在做实验的时候，前两种方法实施起来比较麻烦，要实时初步判断数据的合理性与准确性，建议采用第三种方法，在误差范围内，该方法简单易行，可操作性强，能快速剔除粗差。 完善仪器方面的建议1. 托盘的重量对张力的影响。建议把托盘制成标准的质量，如8g或10g等，方便张力的准确测量和计算。2. 增大振幅。可以通过适当增大电流、提高磁场强度、减小金属丝电阻等方式，增大驻波振动的振幅，从而使驻波的现象更明显，使实验者更容易的观察判断出最大振幅的位置。3. 稳定驻波波形。设法使电流更加稳定，使磁场更加均匀，并减弱托盘的摆动幅度，以达到短时间

20、内使驻波波形稳定的目的。结束语本文通过对大量实验数据进行分析，总结出弦振动时驻波波长与张力，波长与振动频率的规律，提出了检验判断实验数据的方法，可以使实验者在做实验的过程中随时判断数据的正确与否，同时利用MATLAB软件处理数据和作图，使实验规律更加形象化，物理思想更加明确化。 参考文献：1 吴俊林，袁胜利，刘志存，鲁百佐基础物理实验 M陕西：陕西师范大学出版社，2005，78802 钞曦旭，杨万民，唐纯青MATLAB及其在大学物理课程中的应用M西安：陕西师范大学出版社，20063 余凤梅，潘玲珠，冯显灿判定“弦振动研究”实验数据的两种方法J仲恺农业技术学院学报，2005，18(3)：444T

21、he Law of Standing Wave During String Vibration and Data Processing with MATLABMiao Kun1, Huang Yuhong2, Li Kang1, Song Cui1( 1. Department of the Basic Teaching of Science and Engineering ,Shaanxi Normal University, Xian Shaanxi, 700062;2. College of Physics and Information Technology, Shaanxi Norm

22、al University, Xian Shaanxi, 710062)Abstract: In this paper, MATLAB software system is used to process a large number of experimental data. The physical analysis becomes simple and visual by means of this software, thus it can make experimenters pay more attention on the inner physical thought. Furt

23、hermore, the law of the wavelength and tension and the law of the wavelength and frequency are analyzed and summarized. The method to judge the validity of experimental data is putforward, which can help experimenters to judge the data in the process of conducting experiments at any time.Key words:

24、String vibration, Standing wave, MATLAB附录资料：matlab画二次曲面一、螺旋线1.静态螺旋线a=0:0.1:20*pi;h=plot3(a.*cos(a),a.*sin(a),2.*a,b,linewidth,2);axis(-50,50,-50,50,0,150);grid onset(h,erasemode,none,markersize,22);xlabel(x轴);ylabel(y轴);zlabel(z轴);title(静态螺旋线); 2.动态螺旋线t=0:0.1:10*pi;i=1;h=plot3(sin(t(i),cos(t(i),t(i)

25、,*,erasemode,none);grid onaxis(-2 2 -2 2 0 35)for i=2:length(t) set(h,xdata,sin(t(i),ydata,cos(t(i),zdata,t(i); drawnow pause(0.01)endtitle(动态螺旋线);(图略) 3.圆柱螺旋线t=0:0.1:10*pi;x=r.*cos(t);y=r.*sin(t);z=t;plot3(x,y,z,h,linewidth,2);grid onaxis(square)xlabel(x轴);ylabel(y轴);zlabel(z轴);title(圆柱螺旋线) 二、旋转抛物面

26、b=0:0.2:2*pi;X,Y=meshgrid(-6:0.1:6);Z=(X.2+Y.2)./4;meshc(X,Y,Z);axis(square)xlabel(x轴);ylabel(y轴);zlabel(z轴);shading flat;title(旋转抛物面)或直接用：ezsurfc(X.2+Y.2)./4) 三、椭圆柱面load clownezsurf(2*cos(u),4*sin(u),v,0,2*pi,0,2*pi)view(-105,40) %视角处理shading interp %灯光处理colormap(map) %颜色处理grid on %添加网格线axis equal

27、%使x,y轴比例一致xlabel(x轴);ylabel(y轴);zlabel(z轴);shading flat;title(椭圆柱面) %添加标题四、椭圆抛物面b=0:0.2:2*pi;X,Y=meshgrid(-6:0.1:6);Z=X.2./9+Y.2./4;meshc(X,Y,Z);axis(square)xlabel(x轴);ylabel(y轴);zlabel(z轴);shading flat;title(椭圆抛物面)或直接用：ezsurfc(X.2./9+Y.2./4)五、双叶双曲面ezsurf(8*tan(u)*cos(v),8.*tan(u)*sin(v),2.*sec(u),-

28、pi./2,3*pi./2,0,2*pi)axis equalgrid onaxis squarexlabel(x轴);ylabel(y轴);zlabel(z轴);shading flat;title(双叶双曲面)六、双曲柱面load clownezsurf(2*sec(u),2*tan(u),v,-pi/2,pi/2,-3*pi,3*pi)hold on %在原来的图上继续作图ezsurf(2*sec(u),2*tan(u),v,pi/2,3*pi/2,-3*pi,3*pi)colormap(map)shading interpview(-15,30)axis equalgrid onaxi

29、s equalxlabel(x轴);ylabel(y轴);zlabel(z轴);shading flat;title(双曲柱面)七、双曲抛物面（马鞍面）X,Y=meshgrid(-7:0.1:7);Z=X.2./8-Y.2./6;meshc(X,Y,Z);view(85,20)axis(square)xlabel(x轴);ylabel(y轴);zlabel(z轴);shading flat;title(双曲抛物面)或直接用：ezsurfc(X.2./8-Y.2./6) 八、抛物柱面X,Y=meshgrid(-7:0.1:7);Z=Y.2./8;h=mesh(Z);rotate(h,1 0 1,

30、180) %旋转处理%axis(-8,8,-8,8,-2,6);axis(square)xlabel(x轴);ylabel(y轴);zlabel(z轴);shading flat;title(抛物柱面)或直接用：ezsurfc(Y.2./8) 九、环面ezmesh(5+2*cos(u)*cos(v),(5+2*cos(u)*sin(v),2*sin(u),0,2*pi,0,2*pi)axis equalgrid onxlabel(x轴);ylabel(y轴);zlabel(z轴);shading flat;title(环面)十、椭球ezsurfc(5*cos(u)*sin(v),(3*sin(u)*sin(v),4*cos(v),0,2*pi,0,2*pi)axis equalgrid onxlabel(x轴);ylabel(y轴);