1、海洋聲場計算指南BELLHOP手冊和用戶指南Michael B. Porter1-31-2011Heat, Light, and Sound Research, Inc.熱、光和聲音研究公司概 要特征：波束追蹤結(jié)構(gòu)基于幾何和物理擴展規(guī)律導(dǎo)出的一個特別簡單的算法。輸入：應(yīng)用了高斯形狀、草帽形狀等幾種波束。它允許頂端、底端邊界（海面測高和海底測深）和聲速剖面隨距離變化。也允許以描述聲源指向性和邊界介質(zhì)地聲屬性的文件作為輸入附件，還可接納頂端和底端反射系數(shù)文件。輸出：BELLHOP可以生成各種有用的輸出，包括傳播損失、本征聲線、到達聲線、接收時間序列。計算機語言：BELLHOP采用Fortran、M

2、atlab和Python語言實現(xiàn)。計算機平臺：并可應(yīng)用于Mac、Windows和Linux等多種操作系統(tǒng)。1.1 輸入1、最簡單、也最典型的是只提供一個文件，它被稱為環(huán)境文件，包括聲速剖面和海底的信息。2、如果海底隨距離變化，則需添加一個以距離-深度來定義“水深”的測深文件。3、如果海洋聲速隨距離變化，還需添加在規(guī)則網(wǎng)格上標記聲速的SSP文件。4、如果想以反射系數(shù)來描述海底，則需提供以角度-反射系數(shù)來定義“反射特征”的海底反射系數(shù)文件。5、海面情況與此類似。6、如果聲源有指向性，則需提供以角度-幅度來定義“聲源指向性特征”的文件。I、BELLHOP依靠在主環(huán)境文件中設(shè)置選項來讀取以上各類附屬文

3、件。II、有繪圖程序（(plotssp、plotbty、plotbrc等）來展示每一個輸入文件。1.2 輸出BELLHOP依靠在主環(huán)境文件中設(shè)置選項來生成各種各樣的輸出。A、通常從設(shè)置“聲線追蹤”選項開始，并用plotray程序來展現(xiàn)。B、如果對傳播損失（單位強度聲波經(jīng)傳播后的強度），就生成“渲染（shade）”文件。并用plotshd繪制二維表面圖，用plottlr 和 plottld分別繪制深度和距離上的切片。C、想得到單頻聲波經(jīng)傳播后的強度和整個時間序列，就選擇“到達”計算。用plotarr繪制信道的脈沖響應(yīng)序列。經(jīng)卷積求和后，用plotts繪制聲源或接收器的時間序列。2.1 算例120

4、11年原文檔算例。環(huán)境文件MunkB_ray_2011.env1Munk profile ! TITLE 250.0 ! FREQ (Hz) 31 ! NMEDIA 4SVF ! SSPOPT (Analytic or C-linear interpolation) 551 0.0 5000.0 ! DEPTH of bottom (m) 6 0.0 1548.52 / 7 200.0 1530.29 / 8 250.0 1526.69 / 9 400.0 1517.78 / 10 600.0 1509.49 / 11 800.0 1504.30 / 121000.0 1501.38 / 13

5、1200.0 1500.14 / 141400.0 1500.12 / 151600.0 1501.02 / 161800.0 1502.57 / 172000.0 1504.62 / 182200.0 1507.02 / 192400.0 1509.69 / 202600.0 1512.55 / 212800.0 1515.56 / 223000.0 1518.67 / 233200.0 1521.85 / 243400.0 1525.10 / 253600.0 1528.38 / 263800.0 1531.70 / 274000.0 1535.04 / 284200.0 1538.39

6、/ 294400.0 1541.76 / 304600.0 1545.14 / 314800.0 1548.52 / 325000.0 1551.91 / 33A 0.0 345000.0 1600.00 0.0 1.0 / 351 ! NSD 361000.0 / ! SD(1:NSD) (m) 3751 ! NRD 380.0 5000.0 / ! RD(1:NRD) (m) 391001 ! NR 400.0 100.0 / ! R(1:NR ) (km) 41R ! R/C/I/S 4241 ! NBeams 43-20.0 20.0 / ! ALPHA1,2 (degrees) 44

7、0.0 5500.0 101.0 ! STEP (m), ZBOX (m), RBOX (km) 第一步：繪制聲速剖面Matlab命令的語法是：plotssp(MunkB_ray_2011) 或者 plotssp MunkB_ray_2011需要提供配合的其他matlab文件如下：read_env.m；read_env_core.m；topbot.m；：crci.m；readsdrd.m；readvector.m；回顯信息最后繪圖MunkB_ray_2011.env-matlab命令行窗口輸出-Munk profileFrequency = 50 Hz Number of media = 1

8、C-Linear approximation to SSP Attenuation units: dB/mkHz VACUUM z (m) alphaR (m/s) betaR rho (g/cm3) alphaI betaI ( Number of pts = 51 Roughness = 0.00 Depth = 5000.00 ) 0.00 1548.52 0.00 1.00 0.0000 0.0000 200.00 1530.29 0.00 1.00 0.0000 0.0000 250.00 1526.69 0.00 1.00 0.0000 0.0000 400.00 1517.78

9、0.00 1.00 0.0000 0.0000 600.00 1509.49 0.00 1.00 0.0000 0.0000 800.00 1504.30 0.00 1.00 0.0000 0.0000 1000.00 1501.38 0.00 1.00 0.0000 0.0000 1200.00 1500.14 0.00 1.00 0.0000 0.0000 1400.00 1500.12 0.00 1.00 0.0000 0.0000 1600.00 1501.02 0.00 1.00 0.0000 0.0000 1800.00 1502.57 0.00 1.00 0.0000 0.000

10、0 2000.00 1504.62 0.00 1.00 0.0000 0.0000 2200.00 1507.02 0.00 1.00 0.0000 0.0000 2400.00 1509.69 0.00 1.00 0.0000 0.0000 2600.00 1512.55 0.00 1.00 0.0000 0.0000 2800.00 1515.56 0.00 1.00 0.0000 0.0000 3000.00 1518.67 0.00 1.00 0.0000 0.0000 3200.00 1521.85 0.00 1.00 0.0000 0.0000 3400.00 1525.10 0.

11、00 1.00 0.0000 0.0000 3600.00 1528.38 0.00 1.00 0.0000 0.0000 3800.00 1531.70 0.00 1.00 0.0000 0.0000 4000.00 1535.04 0.00 1.00 0.0000 0.0000 4200.00 1538.39 0.00 1.00 0.0000 0.0000 4400.00 1541.76 0.00 1.00 0.0000 0.0000 4600.00 1545.14 0.00 1.00 0.0000 0.0000 4800.00 1548.52 0.00 1.00 0.0000 0.000

12、0 5000.00 1551.91 0.00 1.00 0.0000 0.0000 Number of pts = 51 ACOUSTO-ELASTIC half-space 5000.00 1600.00 0.00 1.00 0.0000 0.0000 cLow = 1.0 m/s cHigh = m/sRMax = 1000.000000 km - Number of source depths = 51 Source depths (m) 0.00 . 5000.00 Producing source depths by interpolation between sd(1) and s

13、d(2)- Number of receivers depths = 1001 Receiver depths (m) 0.00 . 100.00 Producing receiver depths by interpolation between rd(1) and rd(2)第二步：計算聲場并檢查輸入文件運行BELLHOP的MATLAB命令是：bellhop( MunkB_ray_2011 )；或者bellhop MunkB_ray_2011 查看回饋文件*.prt回傳輸入數(shù)據(jù)的打印文件MunkB_ray_2011.prt BELLHOP/BELLHOP3D BELLHOP- Munk p

14、rofile frequency = 50.00 Hz Dummy parameter NMedia = 1 Spline approximation to SSP Attenuation units: dB/mkHz VACUUM Depth = 5000.00 m Sound speed profile: z (m) alphaR (m/s) betaR rho (g/cm3) alphaI betaI 0.00 1548.52 0.00 1.00 0.0000 0.0000 200.00 1530.29 0.00 1.00 0.0000 0.0000 250.00 1526.69 0.0

15、0 1.00 0.0000 0.0000 400.00 1517.78 0.00 1.00 0.0000 0.0000 600.00 1509.49 0.00 1.00 0.0000 0.0000 800.00 1504.30 0.00 1.00 0.0000 0.0000 1000.00 1501.38 0.00 1.00 0.0000 0.0000 1200.00 1500.14 0.00 1.00 0.0000 0.0000 1400.00 1500.12 0.00 1.00 0.0000 0.0000 1600.00 1501.02 0.00 1.00 0.0000 0.0000 18

16、00.00 1502.57 0.00 1.00 0.0000 0.0000 2000.00 1504.62 0.00 1.00 0.0000 0.0000 2200.00 1507.02 0.00 1.00 0.0000 0.0000 2400.00 1509.69 0.00 1.00 0.0000 0.0000 2600.00 1512.55 0.00 1.00 0.0000 0.0000 2800.00 1515.56 0.00 1.00 0.0000 0.0000 3000.00 1518.67 0.00 1.00 0.0000 0.0000 3200.00 1521.85 0.00 1

17、.00 0.0000 0.0000 3400.00 1525.10 0.00 1.00 0.0000 0.0000 3600.00 1528.38 0.00 1.00 0.0000 0.0000 3800.00 1531.70 0.00 1.00 0.0000 0.0000 4000.00 1535.04 0.00 1.00 0.0000 0.0000 4200.00 1538.39 0.00 1.00 0.0000 0.0000 4400.00 1541.76 0.00 1.00 0.0000 0.0000 4600.00 1545.14 0.00 1.00 0.0000 0.0000 48

18、00.00 1548.52 0.00 1.00 0.0000 0.0000 5000.00 1551.91 0.00 1.00 0.0000 0.0000 ( RMS roughness = 0.00 ) ACOUSTO-ELASTIC half-space 5000.00 1600.00 0.00 1.00 0.0000 0.0000 _ Number of source depths = 1 Source depths (m) 1000.00 Number of receiver depths = 51 Receiver depths (m) 0.00000 100.000 200.000

19、 300.000 400.000 500.000 600.000 700.000 800.000 900.000 1000.00 1100.00 1200.00 1300.00 1400.00 1500.00 1600.00 1700.00 1800.00 1900.00 2000.00 . 5000.00000 Number of receiver ranges = 1001 Receiver ranges (km) 0.00000 0.100000 0.200000 0.300000 0.400000 0.500000 0.600000 0.700000 0.800000 0.900000

20、 1.00000 1.10000 1.20000 1.30000 1.40000 1.50000 1.60000 1.70000 1.80000 1.90000 2.00000 . 100.000000 Ray trace run Geometric hat beams in Cartesian coordinates Point source (cylindrical coordinates) Rectilinear receiver grid: Receivers at rr( : ) x rd( : ) Number of beams in elevation = 41 Beam tak

21、e-off angles (degrees) -20.0000 -19.0000 -18.0000 -17.0000 -16.0000 -15.0000 -14.0000 -13.0000 -12.0000 -11.0000 -10.0000 -9.00000 -8.00000 -7.00000 -6.00000 -5.00000 -4.00000 -3.00000 -2.00000 -1.00000 0.00000 . 20.000000000000000 _ Step length, deltas = 0.0000000000000000 m Maximum ray depth, Box%

22、z = 5500.0000000000000 m Maximum ray range, Box%r = 101000.00000000000 m Step length, deltas = 500.00000000000000 m (automatically selected)Tracing beam 1 -20.00Tracing beam 2 -19.00Tracing beam 3 -18.00Tracing beam 4 -17.00Tracing beam 5 -16.00Tracing beam 6 -15.00Tracing beam 7 -14.00Tracing beam

23、8 -13.00Tracing beam 9 -12.00Tracing beam 10 -11.00Tracing beam 11 -10.00Tracing beam 12 -9.00Tracing beam 13 -8.00Tracing beam 14 -7.00Tracing beam 15 -6.00Tracing beam 16 -5.00Tracing beam 17 -4.00Tracing beam 18 -3.00Tracing beam 19 -2.00Tracing beam 20 -1.00Tracing beam 21 0.00Tracing beam 22 1.

24、00Tracing beam 23 2.00Tracing beam 24 3.00Tracing beam 25 4.00Tracing beam 26 5.00Tracing beam 27 6.00Tracing beam 28 7.00Tracing beam 29 8.00Tracing beam 30 9.00Tracing beam 31 10.00Tracing beam 32 11.00Tracing beam 33 12.00Tracing beam 34 13.00Tracing beam 35 14.00Tracing beam 36 15.00Tracing beam

25、 37 16.00Tracing beam 38 17.00Tracing beam 39 18.00Tracing beam 40 19.00Tracing beam 41 20.00 CPU Time = 0.109 s第三步：繪制聲線軌跡MATLAB命令繪制聲線：Figure； plotray( MunkB_ray_2011 )；或者 plotray MunkB_ray_2011距離軸單位為米。如果更喜歡千米，那么就簡單地設(shè)置MATLAB全局變量：global units;units = km;依據(jù)聲線是觸碰單側(cè)邊界還是觸碰兩側(cè)邊界，在繪圖中使用不同顏色來表達。2.2 案例22018.8

26、.3算例。環(huán)境文件MunkB_ray.envMunk profile ! TITLE50.0 ! FREQ (Hz)1 ! NMEDIACVF ! SSPOPT (Analytic or C-linear interpolation)51 0.0 5000.0 ! DEPTH of bottom (m) 0.0 1548.52 / 200.0 1530.29 / 250.0 1526.69 / 400.0 1517.78 / 600.0 1509.49 / 800.0 1504.30 / 1000.0 1501.38 / 1200.0 1500.14 / 1400.0 1500.12 / 16

27、00.0 1501.02 / 1800.0 1502.57 / 2000.0 1504.62 / 2200.0 1507.02 / 2400.0 1509.69 / 2600.0 1512.55 / 2800.0 1515.56 / 3000.0 1518.67 / 3200.0 1521.85 / 3400.0 1525.10 / 3600.0 1528.38 / 3800.0 1531.70 / 4000.0 1535.04 / 4200.0 1538.39 / 4400.0 1541.76 / 4600.0 1545.14 / 4800.0 1548.52 / 5000.0 1551.9

28、1 /A 0.0 5000.0 1600.00 0.0 1.0 /2 ! NSD1000.0 4000 / ! SD(1:NSD) (m)51 ! NRD0.0 5000.0 / ! RD(1:NRD) (m)1001 ! NR0.0 100.0 / ! R(1:NR ) (km)R ! R/C/I/S41 ! NBeams-20.0 20.0 / ! ALPHA1,2 (degrees)0.0 5500.0 25 101.0 ! STEP (m), ZBOX (m), RBOX (km)第一步：繪制聲速剖面plotssp(MunkB_ray)所需其他matlab文件：read_env.m；r

29、ead_env_core.m；topbot.m；crci.m；readsdrd.m；readvector.m；程序運行（通關(guān)）提醒：回顯信息最后繪圖MunkB_ray.env-matlab命令行窗口輸出-Munk profileFrequency = 50 Hz Number of media = 1 C-Linear approximation to SSP Attenuation units: dB/mkHz VACUUM z (m) alphaR (m/s) betaR rho (g/cm3) alphaI betaI ( Number of pts = 51 Roughness = 0

30、.00 Depth = 5000.00 ) 0.00 1548.52 0.00 1.00 0.0000 0.0000 200.00 1530.29 0.00 1.00 0.0000 0.0000 250.00 1526.69 0.00 1.00 0.0000 0.0000 400.00 1517.78 0.00 1.00 0.0000 0.0000 600.00 1509.49 0.00 1.00 0.0000 0.0000 800.00 1504.30 0.00 1.00 0.0000 0.0000 1000.00 1501.38 0.00 1.00 0.0000 0.0000 1200.0

31、0 1500.14 0.00 1.00 0.0000 0.0000 1400.00 1500.12 0.00 1.00 0.0000 0.0000 1600.00 1501.02 0.00 1.00 0.0000 0.0000 1800.00 1502.57 0.00 1.00 0.0000 0.0000 2000.00 1504.62 0.00 1.00 0.0000 0.0000 2200.00 1507.02 0.00 1.00 0.0000 0.0000 2400.00 1509.69 0.00 1.00 0.0000 0.0000 2600.00 1512.55 0.00 1.00

32、0.0000 0.0000 2800.00 1515.56 0.00 1.00 0.0000 0.0000 3000.00 1518.67 0.00 1.00 0.0000 0.0000 3200.00 1521.85 0.00 1.00 0.0000 0.0000 3400.00 1525.10 0.00 1.00 0.0000 0.0000 3600.00 1528.38 0.00 1.00 0.0000 0.0000 3800.00 1531.70 0.00 1.00 0.0000 0.0000 4000.00 1535.04 0.00 1.00 0.0000 0.0000 4200.0

33、0 1538.39 0.00 1.00 0.0000 0.0000 4400.00 1541.76 0.00 1.00 0.0000 0.0000 4600.00 1545.14 0.00 1.00 0.0000 0.0000 4800.00 1548.52 0.00 1.00 0.0000 0.0000 5000.00 1551.91 0.00 1.00 0.0000 0.0000 Number of pts = 51 ACOUSTO-ELASTIC half-space 5000.00 1600.00 0.00 1.00 0.0000 0.0000 cLow = 2.0 m/s cHigh

34、 = m/sRMax = 1000.000000 km - Number of source depths = 51 Source depths (m) 0.00 . 5000.00 Producing source depths by interpolation between sd(1) and sd(2)- Number of receivers depths = 1001 Receiver depths (m) 0.00 . 100.00 Producing receiver depths by interpolation between rd(1) and rd(2)第二步：計算聲場

35、計算文件：atWin10_2018_7atbinbellhop.exebellhop( MunkB_ray )第三步：繪制聲線軌跡Matlab命令：global units;units = km;Figure;plotray( MunkB_ray )3 繪制本征聲線MunkB_eigenray.envMunk profile! TITLE50.0! FREQ (Hz)1! NMEDIACVF! SSPOPT (Analytic or C-linear interpolation)51 0.0 5000.0! DEPTH of bottom (m) 0.0 1548.52 / 200.0 153

36、0.29 / 250.0 1526.69 / 400.0 1517.78 / 600.0 1509.49 / 800.0 1504.30 / 1000.0 1501.38 / 1200.0 1500.14 / 1400.0 1500.12 / 1600.0 1501.02 / 1800.0 1502.57 / 2000.0 1504.62 / 2200.0 1507.02 / 2400.0 1509.69 / 2600.0 1512.55 / 2800.0 1515.56 / 3000.0 1518.67 / 3200.0 1521.85 / 3400.0 1525.10 / 3600.0 1

37、528.38 / 3800.0 1531.70 / 4000.0 1535.04 / 4200.0 1538.39 / 4400.0 1541.76 / 4600.0 1545.14 / 4800.0 1548.52 / 5000.0 1551.91 /A 0.0 5000.0 1600.00 0.0 1.0 /1! NSD1000.0 /! SD(1:NSD) (m)1! NRD800.0 /! RD(1:NRD) (m)1! NR100.0 /! R(1:NR ) (km)E! R/C/I/S5001! NBeams-25.0 25.0 / ! ALPHA1,2 (degrees)0.0

38、5500.0 101.0! STEP (m), ZBOX (m), RBOX (km)4 傳播損失4.1 相干傳播損失（TL）基本算例環(huán)境文件：MunkB_Coh.envMunk profile, coherent! TITLE50.0! FREQ (Hz)1! NMEDIACVW! SSPOPT (Analytic or C-linear interpolation)51 0.0 5000.0! DEPTH of bottom (m) 0.0 1548.52 / 200.0 1530.29 / 250.0 1526.69 / 400.0 1517.78 / 600.0 1509.49 / 8

39、00.0 1504.30 / 1000.0 1501.38 / 1200.0 1500.14 / 1400.0 1500.12 / 1600.0 1501.02 / 1800.0 1502.57 / 2000.0 1504.62 / 2200.0 1507.02 / 2400.0 1509.69 / 2600.0 1512.55 / 2800.0 1515.56 / 3000.0 1518.67 / 3200.0 1521.85 / 3400.0 1525.10 / 3600.0 1528.38 / 3800.0 1531.70 / 4000.0 1535.04 / 4200.0 1538.3

40、9 / 4400.0 1541.76 / 4600.0 1545.14 / 4800.0 1548.52 / 5000.0 1551.91 /A 0.0 5000.0 1600.00 0.0 1.8 0.8 /1! NSD1000.0 /! SD(1:NSD) (m)201! NRD0.0 5000.0 /! RD(1:NRD) (m)501! NR0.0 100.0 /! R(1:NR ) (km)C ! R/C/I/S0! NBEAMS-20.3 20.3 /! ALPHA1, 2 (degrees)0.0 5500.0 101.0! STEP (m), ZBOX (m), RBOX (k

41、m)小貼士：4.2 相干、半相干和非相干傳播損失（TL）Runtype=“C”就執(zhí)行一次所謂的“相干”TL計算。簡單地將該選項第一個字母改為“S”或“I”，我們就分別執(zhí)行“半相干”和“非相干”TL計算。5 指向性聲源均勻介質(zhì)、自由空間、全向點源的環(huán)境文件 omni.envPoint source in free space100.0 ! Frequency (Hz)1CAF-10000.0 1500.0 0.0 1.0 /1500 0.0 10000.0-10000.0 1500.0 0.0 1.0 / 10000.0 /A 0.0/1 ! NSD0.0 / ! SD(1:NSD)501 !

42、NRD-5000.0 5000.0 / ! RD(1:NRD)501! NR-10.0 10.0 /! R(1:NR ) (km)C ! Run type: Ray/Coh/Inc/Sem0 ! NBEAMS-180 180 /! ALPHA1,2 (degrees)0.0 10001.0 10.0! STEP (m), ZBOX (m), RBOX (km)OPTIONS1CC-線性插值A(chǔ)聲學-彈性半空間。需要另起如表4a所描述的表面參數(shù)行 “SURFACE-LINE”：表面高度z 表面cp 表面cs 表面密度 表面 /FdB/(kmHz)，F(xiàn)指與頻率相關(guān)（Freq. dependent）。

43、OPTIONS2A聲學-彈性半空間（ACOUSTO-ELASTIC）。需要另起一行寫明半空間參數(shù)。參見版塊（6a）的說明。ZB CPB CSB RHOB APB ASB選項A通常用于海底建模。OPTIONS2C相干TL計算聲源具有波束圖案的環(huán)境文件：shaded.envShaded point source in free space100.0 ! Frequency (Hz)1CAF0.0 1500.0 0.0 1.0 /1500 0.0 5000.0-5000.0 1500.0 0.0 1.0 / 5000.0 /A, 0.0/1! NSD0.0 /! SD(1:NSD)501! NRD-

44、5000.0 5000.0 /! RD(1:NRD)501! NR-10.0 10.0 /! R(1:NR ) (km)C * ! Run type: Ray/Coh/Inc/Sem361! NBEAMS-180 180 /! ALPHA1,2 (degrees)0.0 5001.0 10.0! STEP (m), ZBOX (m), RBOX (km)聲源波束指向性圖文件： shaded.sbp37-180 10-170 -10-160 0-150 -20-140 -10-130 -30-120 -20-110 -40-100 -30-90 -50-80 -30-70 -40-60 -20-

45、50 -30-40 -10-30 -20-20 0-10 -10 0 10 10 -10 20 0 30 -20 40 -10 50 -30 60 -20 70 -40 80 -30 90 -50100 -30110 -40120 -20130 -30140 -10150 -20160 0170 -10180 10聲源波束模式圖：Matlab命令：polar(a(:,1)/180*pi,a(:,2)拓展計算空間自由空間的環(huán)境文件： omni_extend.env（范圍擴展到 100km100km）Point source in free space100.0 ! Frequency (Hz)1

46、CAF-100000.0 1500.0 0.0 1.0 /1500 0.0 100000.0-100000.0 1500.0 0.0 1.0 / 100000.0 /A 0.0/1! NSD0.0 / ! SD(1:NSD)5001! NRD-50000.0 50000.0 / ! RD(1:NRD)5001! NR-100.0 100.0 /! R(1:NR ) (km)C ! Run type: Ray/Coh/Inc/Sem0 ! NBEAMS-180 180 /! ALPHA1,2 (degrees)0.0 100001.0 100.0! STEP (m), ZBOX (m), RBO

47、X (km)6 邊界隨距離變化6.1 分段線性邊界：Dickins海山分段線性邊界-Dickins海山：DickinsB.btyL5 0 3000 10 3000 20 500 30 3000100 3000環(huán)境文件：DickinsB.envDickins seamount; unit density! TITLE230.0! FREQ (Hz)1! NMEDIACVW! SSPOPT (Analytic or C-linear interpolation)525 0.0 3000.0! DEPTH of bottom (m) 0 1476.7 / 38 1476.7 / 50 1472.6

48、/ 70 1468.8 / 100 1467.2 / 140 1471.6 / 160 1473.6 / 170 1473.6 / 200 1472.7 / 215 1472.2 / 250 1471.6 / 300 1471.6 / 370 1472.0 / 450 1472.7 / 500 1473.1 / 700 1474.9 / 900 1477.0 /1000 1478.1 /1250 1480.7 /1500 1483.8 /2000 1490.5 /2500 1498.3 /3000 1506.5 /A 0.03000.0 1550.0 0.0 1. 0.5 /1! NSD 18

49、.0 /! SD(1:NSD) (m)601! NRD 0.0 3000.0 /! RD(1:NRD) (m)1001! NR 0.0 100.0 /! R(1:NR) (km)CB! R/C/I/S0! NBEAMS-45.0 45.0 / ! ALPHA1,2 (degrees)0.0 3100.0 101.0 ! STEP (m), ZBOX (m), RBOX (km)需要輔助的matlab文件：readbty.m環(huán)境文件：DickinsB.envDickins seamount! TITLE230.0! FREQ (Hz)1! NMEDIACVW ! SSPOPT (Analytic

50、 or C-linear interpolation)525 0.0 3000.0! DEPTH of bottom (m) 0 1476.7 / 38 1476.7 / 50 1472.6 / 70 1468.8 / 100 1467.2 / 140 1471.6 / 160 1473.6 / 170 1473.6 / 200 1472.7 / 215 1472.2 / 250 1471.6 / 300 1471.6 / 370 1472.0 / 450 1472.7 / 500 1473.1 / 700 1474.9 / 900 1477.0 /1000 1478.1 /1250 1480

51、.7 /1500 1483.8 /2000 1490.5 /2500 1498.3 /3000 1506.5 /A* 0.03000.0 1550.0 0.0 1.5 0.5 /1! NSD 18.0 /! SD(1:NSD) (m)201! NRD 0.0 3000.0 /! RD(1:NRD) (m)1001! NR 0.0 100.0 /! R(1:NR) (km)CB! R/C/I/S0! NBEAMS-89.0 89.0 / ! ALPHA1,2 (degrees)0.0 3100.0 101.0 ! STEP (m), ZBOX (m), RBOX (km)Dickins海山場景下

52、的傳播損失a、BELLHOP， b、RAMMatlab 腳本小貼士：注意：ram.m程序運行時，會產(chǎn)生非常奇妙的結(jié)果；應(yīng)對之策：逐一排查原因，最后落腳到TL的矩陣結(jié)構(gòu)。應(yīng)用多種方法之后，發(fā)現(xiàn)子程序RAMtoSHD.m中其實已經(jīng)有提示。按照提示修改如下，把該語句注釋部分的注釋符號（%）去掉即可，這樣就得到了圖（8）。6.2 繪制單波束Dickins海山-單波束環(huán)境文件：DickinsB_oneBeam.envDickins seamount! TITLE230.0 ! FREQ (Hz)1! NMEDIACVW I! SSPOPT (Analytic or C-linear interpolation)525 0.0 3000.0! DEPTH of bottom (m) 0 14