[水声学原理-课件]-声速剖面插值仿真程序说明文档.docxVIP

声速剖面插值 简介 本算例是对水声学原理第二章关于水下声速问题的部分内容进行仿真，利用 MATLAB 对声速进行插值。 基本原理 声速是影响声波在水中传播的最基本物理量。海水中声速的变化会导致声传播规律的改变，因此，精确的声速数据在理论研究和工程应用中都具有十分重要的意义。 本算例提供两种方法进行声速插值计算。第一种方法利用快速傅里叶插值的方法。利用 FFT 算法把测得的实验数据转换到变换域中，再通过补零的方法， 然后用更多点的傅里叶逆变换变换回来得到更多的数据，其结果相当于是对数据进行升采样。 第二种方法利用分段线性插值的方法。分段线性插值具有计算简单、稳定性好、收敛性好、各小段曲线在连接点上连续、且容易实现等多种优点。具体原理可以参阅各类数值计算的参考书籍，在此不作展开。 数值仿真 仿真参数：声速极小值 1500m/s；声道轴深度 1000m；同时可以调整傅里叶插值的点数，本例程中设置点数为原始数据长度的 50 倍。 仿真结果： 本例程使用 Munk 声速分布作为测量得到的数据。 图 1 利用快速傅里叶插值的方法得到的插值结果 图 2 利用分段线性插值的方法得到的插值结果 结论 (1)从仿真结果中可以看出利用分段线性插值得到的数据较好，由于 FFT 运算的特点，插值中会引入截断效应和混叠，导致插值结果的起伏，但有时傅里叶变换插值也不失为一种方法。 (2)实验测量的声速仅是在某些