基于MATLAB的简谐振动合成图形的动态演示135181224

【作者单位】北京林业大学理学院100083【通讯作者】霍虎E-mail：[e-mail]作者简介：谷晓琪，女(回)，1991年，黑龙江省鸡西市，学生，本科，主要研究方向：电子信息科学与技术。地址：北京市清华东路35号北京林业大学理学院，E-mail：[e-mail]；周修文，女，1992年，山东省青岛市，学生，本科，主要研究方向：电子信息科学与技术。地址：北京市清华东路35号北京林业大学理学院，E-mail：[e-mail]；田振国，男，1992年，湖北，学生，本科，主要研究方向：电子信息科学与技术。地址：北京市清华东路35号北京林业大学理学院，E-mail：[e-mail]。基于MATLAB的简谐振动合成图形的动态演示谷晓琪周修文霍虎*田振国（北京林业大学理学院，北京100083）摘要：利用MATLAB编程制作GUI用户界面，通过用户输入简谐振动的各个参数，得到一维、二维简谐振动合成的图形，并研究不同参数对简谐振动合成的影响，观察合成图像的变化。在此基础上，利用MATLAB实现三维简谐振动合成的仿真实验。【关键字】MATLAB；简谐振动；合成；动态演示【中图分类号】O4-39【文献标识码】A简谐振动的合成是力学的主要内容之一，在大学物理的学习中已经对多个同频率简谐振动的合成进行了较为详细的讨论,但研究仅局限于两、三个不同频率简谐振动在一、二维坐标中的合成情况，三维无法体现。利用MATLAB绘出多不同频率简谐振动在一、二及三维坐标中合成的波形及轨迹,并依据这些波形与轨迹,探讨多个不同频率简谐振动的合成规律。1.研究意义利用MATLABGUI软件的制作，对一维、二维及三维简谐振动的授课方式及仿真实验提供了新的方法。在传统教学的过程当中，简谐振动的合成通过示波器来观察，其合成图像往往与理论值有较大偏差，其主要存在的问题有以下几个方面：（1）课堂时间紧迫，图像合成演示的直观显示多要在实验课当中进行操作，不利于学习的及时记忆和直观了解。（2）电路布线不合理引起的交叉干扰、电感漏磁容易引起合成图像失真【1】。（3）示波器无法合成及演示三维简谐振动的合成图像。（4）示波器探头种类繁多，对于精确的理论图形的演示难以企及，且探头可提供测试需要的保真度往往较低。（5）对于简谐振动合成的计算比较复杂，示波器无法高精度的实现图形的模拟合成。针对以上问题制作的MATLABGUI简谐振动合成的程序，能帮助授课教师在课堂上直观的演示简谐振动的合成，直观对比理论测量与实际测量。2.基本原理2.1一维简谐振动合成原理2.1.1多个一维同频率简谐振动的合成设质点在x方向上同时参与n个同频率简谐振动，振动方程【2】为：)1()cos(lxllwtAx其中，nl....2,1其中xlA为振幅，w为角频率，为初相位，则质点的合振动为：)2()cos(1xxnllwtAxx式中，21121sin)cos(nllxlnllxlxAAA)cossinarctan(11nllxlnllxlxAA分别为和振动的振幅和初相位，由（1）式可知多个一维同频率简谐振动可合成为一个同频率的简谐振动，其轨迹是余弦(或正弦)曲线。利用MATLAB进行合成演示，上述推证可得到证实。2.1.2多个一维不同频率简谐振动的合成一般情况下，多个不同频率简谐振动的合振动不再是简谐振动，而是复杂的运动。利用MATLAB进一步研究可知，多个一维频率比为有理数简谐振动的合振动虽然复杂但具有周期性，而多个一维频率比为无理数简谐振动的合振动则既复杂又无周期性【3】。由图1可知，多个一维同振幅、同相位频率相差不大简谐振动的合成，其结果形成多个大小不一的拍[4]。进一步研究可知，多个一维同振幅、同相位频率相差不大简谐振动的合振动是这些简谐振动两两合成的幅度减小的拍的叠加，其结果形成1n个大小不一的拍——多拍现象。其中，主拍的拍幅很大(为单个简谐振动振幅的n倍)，而次拍的拍幅则比较小【5】。2.2二维简谐振动的合成—李萨茹图形2.2.1相互垂直同频率简谐振动的合成当一个质点同时参与两个不同方向的振动时，一般情况下质点将做平面曲线运动，其运动轨迹的形状将由两个分振动的周期、振幅和它们的相位差决定[6]。沿两个振动的方向分别建立x，y轴，并以质点的平衡位置作为坐标原点，则这两个分振动可分别表示为：)3()cos(11wtAx)4()cos(22wtAy在t时刻，质点的位置可由坐标x，y确定。上述方程是以时间t作为参变量的运动轨迹的参数方程，从中消去t，便得轨迹方程：)5)((sin)cos(21221221222212AAxyAyAx此式是椭圆方程，它表示两个相互垂直且同频率的简谐振动合成的轨迹是椭圆。随着相位差值的不同，合成椭圆的形状也不同【7】。2.2.2相互垂直不同频率简谐振动的合成如果两个相互垂直的振动频率不相同，它们的合运动比较复杂，若随意选取两种分振动，可以看到合成轨迹是不稳定的，而且没有规律可循。（1）两振动的频率有很小的差异，可近似看成同频率振动的合成，不过相位差在缓慢地变化，在范围内由直线变成椭圆再变成直线等。（2）如果两振动频率相差较大，但有简单整数比，则合成运动有稳定封闭的运动轨迹，称为李萨如图形，其形状与互相垂直的分振动的角频率之比、各自初相位以及初相位差都有关系[5]。图2为两振动的频率有简单的整数比，合成运动为李萨如图形。运行中可以看到以质点运动轨迹的方式呈现的动画。图1一维多个不同频率简谐振动合成2.3三维简谐振动的合成设分别沿z、y和z向的频率不同的简谐振动的表达式为：)cos(111twAx)cos(222twAy)cos(333twAz此方程组就是合运动轨迹的参数方程。设x、y和z向振动满足频率比为整数比：2121nnvv3131nnvv其中1n，2n为不可约的整数，1n和3n为不可约的整数，那么存在：1nTT其中1T是x方向的振动周期，n是1n和2n的最小公倍数，T就是(1)式描述的空间曲线的参数周期，也就是上述曲线为闭合曲线。以x、y和z为坐标的空间点在时间T内完成闭合曲线的一次扫描，然后重复扫描，这就是三维李萨如曲线。如果不满足(2)式，那么(1)式描述的曲线不是闭合曲线【8】。利用MATLAB进行三维简谐振动的合成，得到图3。可看出，当分振动为简单整数比时，三维合成的图像也有类似李萨茹图像的情况。3.GUI界面及仿真流程利用MATLAB制作的用户界面人工操作，用户通过主页面选择所需功能，进入各功能页面后完成各函数读取目标，完成函数识别与图像的合成。整个仿真实验的流程如图4。4.结束语该实验通过MATLAB制作，利用MATLAB作为仿真平台，动态显示多个不同频率的简谐振动、互相垂直的不同频率的两个简谐振动、三维简谐振动的合成，让学生能在计算机上形象直观的认识并掌握三维的简谐振动的合成图象，尤其是多个、稍大频率的李萨如图形的合成，得出大频率李萨如图形仿真实验图象。参考文献：[1]百度文库.示波器知识百问[EB/OL].[2011-1-12].[2]马文蔚等.物理学[M].北京市:高等教育初版社2001:图2一维多个不同频率简谐振动合成图3三维简谐振动的合成图4仿真流程图设置个控件属性编辑各方向的简谐振动开始扫描各输入参数的调用comet[9]函数进行合成送用axes显示图像结束134-136.[3]蓝海江.多个简谐振动的合成[J].广西科学院学报,2009,1:22-25.[4]王宁星.拍图形的对称性[J].广州大学学报，2002，1(5)：89-91.[5]吴俊，张毓麟，晏世雷等.三个简谐振动合成的拍现象和实验仪[J].物理与工程，2008(03):38-39.[6]刘慧，王玉连，李勇.MATLAB在简谐振动的合成教学中的应用[J].科技信息博士·专家论坛，2010(20)：434-436.[7]熊德永，马慧.互相垂直同频率简谐振动的合成轨迹方程的推导[J].安徽科技学院学报，2012(02)：27[8]康文秀.三维李萨如图[J].物理与工程，2004(01):11-12.[9]张国辉.基于MATLAB的函数作图[J].湖南城市学院学报，2009，18(1):35-37.StudyoncomposographofsimpleharmonicvibrationusingMATLABGUXiao-qi,ZHOUXiu-wen,Huohu*,TIANZhen-guo(Collegeofscience,BeijingForestryUniversity,Beijing100083,China)Abstract:BasedonMATLABprogramming,producingtheGUIuserinterface.Userinputsparametersandcouldgetone-dimensional,two-dimensionalharmonicsynthesisofgraphics,andtostudytheeffectofdiffe-rentparametersonthecompositionofsimpleharmonicvibrationeffects,observationofsyntheticimagechanges.Onthisbasis,usingMATLABtoachievethree-dimensionalcompositionofsimpleharmonicvibra-tionsimulationexperiment..Keywords:MATLAB;harmonicvibration;synthesis;dynamicdemonstration全文总字数:3308;谷晓琪GuXiaoqi出生年：1991年出生地：黑龙江省鸡西市工作单位：北京林业大学；通信地址：北京市海淀区清华东路35号北京林业大学；邮编：100083；联系电话：[phone]；E-mail：[e-mail]周修文ZhouXiuwen出生年：1992年出生地：山东省青岛市工作单位：北京林业大学；通信地址：北京市海淀区清华东路35号北京林业大学；邮编：100083；联系电话：[phone]；E-mail：[e-mail]田振国TianZhenguo出生年：1992年出生地：湖北省工作单位：北京林业大学；通信地址：北京市海淀区清华东路35号北京林业大学；邮编：100083；联系电话：[phone]；E-mail：[e-mail]通讯作者霍虎HuoHu出生年：1972年出生地：北京市工作单位：北京林业大学；通信地址：北京市海淀区清华东路35号北京林业大学；邮编：100083；联系电话：[phone]；E-mail：[e-mail]