matlab声压,MATLAB模拟声速测量中声压的变化

声速测定广泛应用于测距、无损检测、液体测速、气体成分分析等生产生活领域,超声测量声速也就成为大学物理的一项重要实验内容.声波在介质中传播的能量损耗以及在介质界面反射的能量损耗是实验误差产生的主要原因[1-3],声波的强度变化体现在声压的变化上[4-8].本文首先利用驻波法测量声压变化,再利用MatLab软件对实验数据进行分析模拟[9-10],最大精度地测量声压变化,减少实验误差.实验结果表明,理论模拟分析与实验结果吻合度很高,此法可以应用于大学物理声速测量实验教学.1声速测量中声压的变化规律测试过程中,示波器给出的波形是声压波形.由于接收器受到超声波的声压作用,所以接收器能够接收到叠加波的信号.对于理想状态下的一维情况,声压与空气分子的位移关系为p=ksyx,(1)其中ks为媒质的弹性系数.则总的合成波y=y++y-=A+cos(x-2x)+A-cos(x+2x)=2A-cos2cosx+(A+-A-)cos(t-2x).(2)将(2)式代人(1)式,可得到各点处的声压波动方程为p=ksyx=-4A-sin2xcost+2(A+-A-)sin(t-2x)(3)仅在A+=A-时,合成声压波才是纯粹的驻波,p=-4A-sin2xcost.(4)驻波的波峰或者波谷与换能器距离呈周期性关系.当两列波实现共振时,驻波才有最大的振幅,声压最大,示波器上显示为波峰;当不共振时,显示波形随之减小,声压减小;当两换能器之间的距离改变了半个波长,示波器上就会连续出现相邻的两个波峰或者波谷.测试中任意两个波峰或者波谷所对应的距离移动仍为半个波长的整数倍,两个相邻的波峰之间对应的距离移动为半个波长.2测试中声压的变化规律采用驻波法测定声波波长,频率直接由信号发生器给出,可以最大精度地测量声压变化.由于声波在空气中衰减较大,随着距离增加,示波器显示波形振幅会越来越小,如图1所示.图1声压与换能器间距之间的关系3声压波曲线的计算机模拟与分析两换能器距离L与接收端换能器接收到的声压能量正比于声压振幅的平方A2,其关系推导为A2=A02e-2aL(1+r2)21-2r1r2e-2aLcos2ML+r12r22e-4aL,(5)式(5)中,a为空气的吸收系数,A0为入射波振幅,r1与r2分别为波在发射端和接收端的反射系数,M为波数.编写Matlab语句,通过数值模拟获得E-L关系曲线,见图2.声压振幅平方A2随着接收换能器位置L的变化而呈现周期性减弱,每间隔/2出现一次极大或极小值.声压振幅取极值的位置满足L=n/2(n=0,1,2,3,…)图2MATLAB模拟的E-L关系由驻波法得到的驻波波峰值点记录见表1,应用Matlab模拟获得声压波理论峰值点记录见表2,两者的拟合曲线如图3所示.表1声压波实验峰值点Ln/mm0.003.617.8412.4716.5220.8124.9429.2833.2037.42Un/mv10.007.425.644.924.023.162.142.081.981.72表2声压波理论峰值点Ln/mm0.004.208.4012.5016.4120.5124.6128.7332.8236.93Un/mv10.007.125.324.483.903.272.862.522.182.04图3实验与模拟峰值拟合曲线比较由表1的测试数据经拟合后,可获得测试方程f(x)=9.453*exp(-0.05341*x)即f(x)=9.453e-0.05341x.(5)由表2的Matlab模拟峰值点经拟合后,可获f(x)=9.221*exp(-0.05013*x),即f