水声学原理4.1.5表面声道-临近海面的水下点源声场仿真程序.pdf

%邻近海面的水下点源声场的仿真 %对应水声学原理第四章4.1 节的内容 close all clear all clc c=1500;%水中声速 f=1000;%声源频率 w=2*pi*f;%角频率 k=w/c;%水中波数 z1=50;%声源深度 Lmin=20;%最小接收水平距离 Lmax=8000;%最远接收水平距离 z=30;%接收深度 llim=2*k*z*z1/pi;%临界距离 r=Lmin:1:Lmax; R1=Lmin:1:llim; p1=2./R1.*abs(sin(k*z1*z./R1))+10^-5;%菲涅尔干涉区声场 R2=llim:1:Lmax; p2=2*k*z*z1./(R2.^2)+10^-5;%夫琅和费区声场 rep=p2(1)/p1(length(p1));%修正系数避免出现跳变 p2=p2/rep; P1=-20*log10(p1); P2=-20*log10(p2); plot(R1,P1);hold on;%绘图 plot(R2,P2);hold on; title('临近海面水下点源声场传播损失') xlabel('Range/m') ylabel('TL/dB') axis ij