%?******************************************************
%?FFP
%?******************************************************
clear?all

%?***?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??=?500;			%?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（?‘Depth?（m）‘?）;?ylabel（?‘Sound?Speed?（m/s）‘?）

%?******************************************************
%?set?up?finite?difference?matrix
%?******************************************************

v?=?omega?*?omega?./?（?c?.*?c

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件???????3379??2009-07-15?20:59??海洋聲傳播程序\ffp.m

?????文件???????2826??2009-07-15?20:59??海洋聲傳播程序\modes.m

?????文件???????2922??2009-07-15?20:59??海洋聲傳播程序\pe.m

?????文件????????463??2009-07-15?20:58??海洋聲傳播程序\rayf.m

?????文件????????942??2009-07-15?20:58??海洋聲傳播程序\rays.m

?????目錄??????????0??2010-11-12?16:27??海洋聲傳播程序

-----------?---------??----------?-----??----

????????????????10532????????????????????6