实验二-分层介质中的声线传播路径研究.docxVIP

实验二 分层介质中的声线传播路径研究一、实验目的1、了解海水介质的垂直分层特性；2、掌握Snell折射定律和声线弯曲规律；3、掌握基本的声线传播轨迹的搜索方法；二、实验仪器 计算机三、实验原理在本课程的研究中，一般只考虑海水介质的垂直分层的特性，即声速（折射率）不随水平方向变化，仅是海水深度的函数。若令x、y代表水平坐标，z代表垂直坐标，则在分层介质中：c(x, y, z) = c(z) （1）n(x, y, z) = n(z) （2）除了超远程声传播问题外，海水介质的分层模型是对实际海洋介质的一种接近的理想模型。本课程中所获得的结果，都是在分层介质模型下得到的。3.1 Snell折射定律和声线弯曲已知射线声学所遵循的Snell定律为： （3）α为声传播方向与水平坐标ox的夹角，称为掠射角；c为该处声速。α0和c0为声线出射处的对应值。若α0和声线的垂直分布c(z)已知，可以按照Snell定律求出海洋中任意深度处声线的掠射角，也就确定了任意深度处声波传播方向。不同的起始掠射角α0求出不同的α值，对应不同的声线轨迹。负梯度下声速随深度增加而下降，掠射角α随深度增加而增加，声线弯向海底；正梯度下声速随深度增加而增加，掠射角α随深度增加而减小，声线弯向海面。3.2 声线轨迹二维平面问题的声线曲率的表达式为： （4）（4）式中，θ为声线入射角，即声线长度微元ds与垂直轴线z的夹角；c为微元ds处的声速。由Snell定律可知，当起始值给定时（α0和c0），比值cosα/c等于常数。对于恒定声速梯度而言c=c0(1+az)，a为相对声速梯度（a为常数），则dc/dz =c0a=常数，即恒定声速梯度情况下声线曲率到处相等，轨迹是圆弧。3.3 声线轨迹方程对于恒定声速梯度c=c0(1+az)，a=1/c0·dc/dz m。若声线在海面以掠射角α0=0（θ=π/2）出射，该声线曲率半径等于： （5）进一步可写出以α0=0掠射的声线轨迹方程： （6）对于声源在海面以任意掠射角α1出射的声线，同样可给出其声线轨迹方程： （7）3.4 声线经过的水平距离若声源位于(x, z1)处，接收点位于(x, z)处，声速按c=c(z)分布，可以从下列积分求出声线经过的水平距离： （8）根据Snell定律可导出，其中n(z)=c(z1)/c(z)。因而水平距离： （9）对于恒定声速梯度情况，轨迹是一圆弧，直接从声线轨迹图来求水平距离x更为方便。3.5 声线经过的水平距离声线因深度而异，声线经过微元ds距离所需要的时间dt=ds/c，声线从z1声度传播到z深度所需要的时间等于： （10）根据Snell定律，，其中n=c(z1)/ c(z)= c/ c，则： （11）（11）式为计算声线传播时间的一般表示式。在恒定声速梯度下，从Snell定律可求得，代入公式（10）中得到： （12）α1和α分别为深度z1和z处的声线掠射角。对（12）式积分可得： （13）四、实验内容 画出仿真条件下，从声源出发到接收点处三条主要声线的传播规律，三条声线分别直达声信号、海面一次反射信号、海底一次反射信号。仿真条件1：仿真条件2：实验步骤mk：function y=mk(x1,y1,x2,y2,x) k=(y2-y1)/(x2-x1); y=k*(x-x2)+y2;Endmp：function [m,n]=mp(a,b) syms x y; eq1=(y-b)/(x-a)-100; eq2=sqrt(x^2+(y-50)^2)-sqrt(a^2+(b-50)^2); s=solve(eq1,eq2,x,y); m=double(s.x); n=double(s.y); endclear;clc;axis([0 1000 0 50]);x=0:1000;y1=abs(mk(0,-20,1000,10,x));y2=abs(mk(0,-20,1000,40,x));y3=mk(0,80,1000,10,x);y4=mk(0,80,1000,40,x);y5=mk(0,20,1000,10,x);y6=mk(0,20,1000,40,x);for i=1:1001 if y3(i)50 y3(i)=100-y3(i); end if y4(i)50 y4(i)=100-y4(i); endendy=[y1;y2;y3;y4;y5;y6];for j=1:6 plot(x,y(j,:));hold onendset(gca,xaxislocation,top,ydir,reverse); 图1clear;clc;axis([0 1000 0 50]);x=0:1000;y0=mk(0,5