简单PE模型(Munk声速剖面)-程序员宅基地

% ******************************************************
% PE
% ******************************************************

% *** set up Munk SSP ***

zs  = 1000.0;            % source depth
f   = 50; omega = 2 * pi * f    % frequency in Hertz
eps = 0.00737; c0  = 1500;

d   = 5000;            % bottom depth
nz  = 250;            % number of finite-difference points
h   = d / nz; h2  = h * h;    % mesh spacing
z   = linspace( 0, d, nz );    % grid coordinates

x = 2 * ( z - 1300 ) / 1300;  
c = c0 * ( 1 + eps * ( x - 1 + exp( -x ) ) );

plot( z, c ); view( 90, 90 );
xlabel( '