声传播理论模型.docVIP

声道中的声线绘制与频散方程数值求解 作者：王翰卓 模型内容 给定大洋某处深度上离散采集的声速值，某深度处存在声速的极小值点，因为存在声道。认为在大洋一定范围内声速在水平方向上没有变化。解决以下问题- 画出连续的声速剖面曲线 当在适当深度上放置声源，使得当声波掠射角在之间时，所有的声线都会上下发生反转而不发生海底或海面的反射。应用射线声线的模型，画出其在给定距离上的声线。 将射线声线与波动声学相联系，在只存在声线反转的情况下，给定声波频率，应用频散方程，求解特征波数。 算法分析 1、给定深度-声速曲线上的离散的点，可应用分段样条插值的方法得到深度上更多点的声速信息，同时保证了曲线的平滑。 2、算法如下所示 （1）输入传播距离，设置声源深度。 （2）对深度进行离散，对每个声层上的声速进行线性插值。 （3）计算声源处的声速值，用找出每条声线翻转处的声速，从声源处分开上下两层，分别用线性插值找出不同声线的上下翻转深度。 （4）对每条声线进行计算。先求出声线每层分层处的掠射角，在计算其所走的水平距离；每个角度值记录其从出发到其到达第一次翻转位置处的每层的水平位移。 （5）根据同一深度上某一声线的略射角相同，即声线的轴对称性，得到声线从第一次反转点到回到声源高度位置上的水平位移信息，最终得到一个跨度上的水平位移信息。 （6）进行画图。与的声线形状是一致的，故的上层可以用的上层画出，的下层可以用的的下层画出；在利用每条声线的周期性，可将其延拓到指定水平距离。 3、问题分析 波矢的方向即射线的某点的切线，水平方向的波矢具有不变性，一定时，由于，由声源极小值处发射的声线在向上下传播的过程中由于变小，相应与的夹角变小，直到即达到射线的反转点。 为了使得所有的声线都反转，而且海底声速大于海面声速，故；同时因为，所以。越小，反转点离声源越远。 当使得成立时，其为本征波数，显然，一定时，越小，越大，且积分上下限变长，简正波号数增大。增大时，简正波的最大号数也会增加。 编程的算法如下： 给定，计算分层后的声速以及波数。 计算声速最小值，给出的范围。 用计算最大的简正波号数。 用二分法在中找的根，积分的计算选用复化辛普生公式，最终确定每一号水平波数。 源代码及运行结果 1、声线绘制 clc,clear all,format long T=tic; depth1=[0.0,150.0,305.0,533.0,610.0,680.0];depth2=[762.0,1372.0,1829.0,3048.0,4000.0]; depth=[depth1 depth2]; c1=[1507.2,1498.1,1491.7,1480.7,1478.9,1478.0];c2=[1478.6,1483.2,1488.6,1507.5,1523.0]; c=[c1 c2]; z0=700;%声源位置 c0=interp1(depth,c,z0); x0=input(input the distance of propagation); z1=0:700;z2=700:4000;z=[z1 z2];%深度上进行分层 figure(1) plot(spline(depth,c,z),z,r);%样条差值后画出声速-深度曲线 set(gca,yDir,reverse);set(gca,XAxisLocation,top);%y轴反转，x轴取有效部分 title(sound speed section); xlabel(speed(m/s));ylabel(depth(m)); grid on; c1n=interp1([depth1,z0],[c1,c0],z1);%线性差值，求封层后的声速 c2n=interp1([z0,depth2],[c0,c2],z2); theta=[-pi/18:pi/180:-pi/180,pi/180:pi/180:pi/18];%入射角度 cr=c0./cos(theta); reverse_up=interp1(c1n,z1,cr(1:10));%计算上下翻转深度 reverse_down=interp1(c2n,z2,cr(11:20)); x=zeros(20,2500);%存储每条声线的离散水平分量 number=ceil([700-reverse_up,reverse_down-700])+1;%每条声线记录的基本x点的个数,同时也是声线的高度 x(:,1)=0; for i=1:20 if i=10 beta=acos(c1n(701:-1:701-number(i)+2)*cos(theta(i))/c0);%记录掠射角 beta(length(beta