1.2-BELLHOP-手册-和-程序日志课件

1、海洋声场计算指南BELLHOP手册和用户指南Michael B. Porter1-31-2011Heat, Light, and Sound Research, Inc.热、光和声音研究公司概 要特征：波束追踪结构基于几何和物理扩展规律导出的一个特别简单的算法。输入：应用了高斯形状、草帽形状等几种波束。它允许顶端、底端边界（海面测高和海底测深）和声速剖面随距离变化。也允许以描述声源指向性和边界介质地声属性的文件作为输入附件，还可接纳顶端和底端反射系数文件。输出：BELLHOP可以生成各种有用的输出，包括传播损失、本征声线、到达声线、接收时间序列。计算机语言：BELLHOP采用Fortran、M

2、atlab和Python语言实现。计算机平台：并可应用于Mac、Windows和Linux等多种操作系统。1.1 输入1、最简单、也最典型的是只提供一个文件，它被称为环境文件，包括声速剖面和海底的信息。2、如果海底随距离变化，则需添加一个以距离-深度来定义“水深”的测深文件。3、如果海洋声速随距离变化，还需添加在规则网格上标记声速的SSP文件。4、如果想以反射系数来描述海底，则需提供以角度-反射系数来定义“反射特征”的海底反射系数文件。5、海面情况与此类似。6、如果声源有指向性，则需提供以角度-幅度来定义“声源指向性特征”的文件。I、BELLHOP依靠在主环境文件中设置选项来读取以上各类附属文

3、件。II、有绘图程序（(plotssp、plotbty、plotbrc等）来展示每一个输入文件。1.2 输出BELLHOP依靠在主环境文件中设置选项来生成各种各样的输出。A、通常从设置“声线追踪”选项开始，并用plotray程序来展现。B、如果对传播损失（单位强度声波经传播后的强度），就生成“渲染（shade）”文件。并用plotshd绘制二维表面图，用plottlr 和 plottld分别绘制深度和距离上的切片。C、想得到单频声波经传播后的强度和整个时间序列，就选择“到达”计算。用plotarr绘制信道的脉冲响应序列。经卷积求和后，用plotts绘制声源或接收器的时间序列。2.1 算例120

4、11年原文档算例。环境文件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回传输入数据的打印文件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距离轴单位为米。如果更喜欢千米，那么就简单地设置MATLAB全局变量：global units;units = km;依据声线是触碰单侧边界还是触碰两侧边界，在绘图中使用不同颜色来表达。2.2 案例22018.8

26、.3算例。环境文件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；程序运行（通关）提醒：回显信息最后绘图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）基本算例环境文件：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”就执行一次所谓的“相干”TL计算。简单地将该选项第一个字母改为“S”或“I”，我们就分别执行“半相干”和“非相干”TL计算。5 指向性声源均匀介质、自由空间、全向点源的环境文件 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声学-弹性半空间。需要另起如表4a所描述的表面参数行 “SURFACE-LINE”：表面高度z 表面cp 表面cs 表面密度 表面 /FdB/(kmHz)，F指与频率相关（Freq. dependent）。

43、OPTIONS2A声学-弹性半空间（ACOUSTO-ELASTIC）。需要另起一行写明半空间参数。参见版块（6a）的说明。ZB CPB CSB RHOB APB ASB选项A通常用于海底建模。OPTIONS2C相干TL计算声源具有波束图案的环境文件：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)拓展计算空间自由空间的环境文件： 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环境文件：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环境文件：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程序运行时，会产生非常奇妙的结果；应对之策：逐一排查原因，最后落脚到TL的矩阵结构。应用多种方法之后，发现子程序RAMtoSHD.m中其实已经有提示。按照提示修改如下，把该语句注释部分的注释符号（%）去掉即可，这样就得到了图（8）。6.2 绘制单波束Dickins海山-单波束环境文件：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