[海洋中的声传播理论-课件]-简正波声场计算及时域波形预报仿真程序.docxVIP

程序 1： %理想边界（绝对软海底）平行平面层波导声场计算 %对应水声学原理第三章前两节针对波动声学的部分内容clear all close all clc; f=25;%声源频率 w=f*2*pi;%角频率 c=1500;%声速 k=w/c;%波数 z0=100; %声源深度 L=10000; %所要求解的距离H=200; %最大深度 step=1000; %运算点数pp=zeros(step,step);%声场初始化Tl=zeros(1,8000); lam=c/f;%波长 n=floor(H/(lam/2));%简正波阶数x=0:L/(step-1):L;%水平距离z=0:H/(step-1):H;%深度 for i=1:1:n kzn=i*pi/H; %本征值kn=sqrt(k^2-kzn^2); po=(j*pi*2/H)*sin(kzn*z0)*sin(kzn*z).*besselh(0,2,(kn*x));%声场解析解pp=po+pp;%求和 end; p0=exp(j*k);%声源归一化条件p=abs(pp)+10^-7;%避免奇异性TL=-20*log10(abs(p/p0)); figure pcolor(x,z,TL)%绘图shading interp; colormap(default) cmap = colormap; colormap(flipud(cm