[水声学原理-课件]-深海负梯度和负跃层 - 虚源法声场计算与波形预报仿真程序.docxVIP

程序 1： %%虚源法声场计算%% %对应水声学原理第四章 4.5.3 节的部分内容 clear all close all clc c=1500; %声速f=25; %声源频率 w=2*pi*f;%角频率k=w/c; %波数z0=100; %声源深度H=200; %波导深度 L=10000; %水平最远距离step=1000; %运算点数step2=200; %运算点数x=0:L/(step-1):L; z=0:H/(step2-1):H; m=200; %反射叠加次数 pp=zeros(step2,step);%初始化for n=0:1:m zn1=(2*H*n+z0-z).*ones(1,step); %四条本征声线对应虚源与接收点深度的差值zn2=(2*H*(n+1)-z0-z).*ones(1,step); zn3=(2*H*n+z0+z).*ones(1,step); zn4=(2*H*(n+1)-z0+z).*ones(1,step); r1=sqrt((ones(step2,1)*x).^2+zn1.^2);%四条本征声线对应虚源与接收点之间的传播距离r2=sqrt((ones(step2,1)*x).^2+zn2.^2); r3=sqrt((ones(step2,1)*x).^2+zn3.^2); r4=sqrt((ones(step2,1)*x).^2+zn4