Bellhop 海底地形起伏条件下的传播特性

news2024/12/24 8:54:49

文章目录

  • 前言
  • 一、预备内容
  • 二、水平海底波导(水平海底)
    • 1、海底水平的深海波导中的声线
      • ①、环境文件
      • ②、Matlab 命令
      • ③、执行结果
    • 2、海底水平的深海波导中的本征声线
      • ①、环境文件
      • ②、Matlab 命令
      • ③、执行结果
    • 3、海底水平的深海波导中的相干传播损失
      • ①、环境文件
      • ②、Matlab 命令
      • ③、执行结果
    • 4、到达接收器处的脉冲响应
      • ①、环境文件
      • ②、Matlab 命令
      • ③、执行结果
  • 三、高斯海山(变化的海底)
    • 1、海底形状文件
    • 2、环境文件
    • 3、Matlab 代码
    • 3、执行结果
  • 四、变化的边界(变化的海底和海面)
    • 1、海面形状文件
    • 2、海底形状文件
    • 3、环境文件
    • 4、Matlab 代码
    • 5、执行结果
  • 五、分段线性边界(Dickins 海山)
    • 1、海底形状文件
    • 2、环境文件
    • 3、Matlab 代码
    • 4、执行结果
  • 六、曲线拟合边界(抛物线海底)
    • 1、海底形状文件
    • 2、环境文件
    • 3、Matlab 代码
    • 4、执行结果
  • 七、资源自取
  • 总结


前言

由于水下声信道课程大作业的需要,因此本节专门研究海底地形起伏条件下的声传播特性


一、预备内容

在所有实例中,我们均采用 Munk 深海声速剖面,界于 0 到 5000m 深度之间,声源频率 50Hz,位于 1000m 深度,声线步距 100m,声线70根,出射角扇面为 -13°~ 13°,海底声速1600m/s,海底密度 1.8g/cm3,海底衰减系数 0.8 dB/λ,海洋环境如下图所示:
在这里插入图片描述

应用的深海场景示意图

二、水平海底波导(水平海底)

注:波导含义:海洋波导可以理解为声波在海洋中传播的通道,比如经典的深海 Munk 剖面,声波会被限制在一定范围内的通道远距离传播,这个通道就可以理解为海洋波导

我们以从 0 到 101km 之间简单的声线轨迹(计算)作为开始。

下面是 Matlab 代码,分别绘制了海底水平的深海波导中的声线轨迹海底水平的深海波导中的本征声线海底水平的深海波导中的相干传播损失到达声线(脉冲响应)

clc;    clear;  close all;
global units;	units = 'km';
bellhop flatwav_R   %   Runtype = 'R'
figure; plotray flatwav_R	% 海底水平的深海波导中的声线
bellhop flatwav_E   %   Runtype = 'E'
figure; plotray flatwav_E	% 海底水平的深海波导中的本征声线
bellhop flatwav_C   %   Runtype = 'C'
figure; plotshd flatwav_C.shd	% 海底水平的深海波导中的相干传播损失
bellhop flatwav_A   %   Runtype = 'A'
% [ Arr, Pos ] = read_arrivals_asc( ARRFile, Narrmx )
[ Arr, Pos ] = read_arrivals_asc( 'flatwav_A.arr' );
% plotarr( filename, irr, ird, isd )
plotarr(  'flatwav_A.arr', 1, 1, 1 )	% 到达接收器处的脉冲响应

接下来我们会对上面代码分别进行讲解

1、海底水平的深海波导中的声线

①、环境文件

有关环境文件的具体讲解可以参考我之前的博客-> Bellhop 从入门到上手

flatwav_R.env

'Munk profile/Flat waveguide/Ray trace'  	% TITLE
50.0 										% FREQ (Hz)
1											% NMEDIA
'SVW'										% SSP-TOP-WATER-OPT
51 0.0 5000.0								% NMESH SIGMA Z(NSSP)
0.0 1548.52 /								% Z() CP() CS() RHO() AP() AS()
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 /
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.39 /
4400.0 1541.76 /
4600.0 1545.14 /
4800.0 1548.52 /
5000.0 1551.91 /
'A' 0.0										% BOTOPT SIGMA
5000.0 1600.00 0.0 1.8 .8 0.0				% ZB CPB CSB RHOB APB ASB
1 											% NSD
1000.0 /									% SD(1:NSD) (m)
1 											% NRD
1000.0 /									% RD(1:NRD) (m)
1 											% NRR
101.0 /										% RR(1:NRR ) (km)
'R'											% OPTION: 'R/E/C/A/I/S
70 											% NBEAMS ISINGLE
-13.0 13.0 /								% ALPHA(1:NBEAMS) (°)
100.0 5500.0 102.0							% STEP (m) ZBOX (m) RBOX (km)

②、Matlab 命令

bellhop flatwav_R   %   Runtype = 'R'
figure; plotray flatwav_R	% 海底水平的深海波导中的声线

一旦 Bellhop 完成计算,可以查验是否创建了两个文件:第一个是flatwav_R.prt,包含波导特性、出射角数目、计算时间等相关的综合信息;第二个是 flatwav_R.ray,它是包含射线坐标的 ASCII 码文件,并可用 M-文件 plotray.m 来绘进行图。

③、执行结果

在这里插入图片描述

海底水平的深海波导中的声线轨迹

2、海底水平的深海波导中的本征声线

将 OPTIONS3(1) = ‘R’ 改为 OPTIONS3(1) = ‘E’,和前面同样地运行程序,可以得到本征声线。

①、环境文件

将上面 flatwav_R.env 第 41 行的 ‘R’ 改为 ‘E’ 即可

'Munk profile/Flat waveguide/Eigenrays'
50.0
1
'SVW'
51 0.0 5000.0
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 /
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.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.8 .8 0.0
1 
1000.0 /
1 
1000.0 /
1 
101.0 /
'E'
70 
-13.0 13.0 /
100.0 5500.0 102.0

②、Matlab 命令

bellhop flatwav_E   %   Runtype = 'E'
figure; plotray flatwav_E

③、执行结果

在这里插入图片描述

海底水平的深海波导中的本征声线

3、海底水平的深海波导中的相干传播损失

计算相干传播损失需要对输入文件做点小修改。首先,设置OPTIONS3(1) = ‘C’;其次,需要考虑是在 501×501 点的距离-深度矩形网格上计算传播损失。最后,我们设置 nbeams = 0,让 Bellhop 自行决定所需的射线根数。

①、环境文件

flatwav_C.env

'Munk profile/Flat waveguide/Coherent transmission loss'
50.0
1
'SVW'
51 0.0 5000.0
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 /
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.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.8 .0 .0 /
1
1000.0 /
501
0.0 5000.0 /
501
0.0 101.0 /
'C'
0
-14.0 14.0 /
100.0 5500.0 102.0

②、Matlab 命令

bellhop flatwav_C   %   Runtype = 'C'
figure; plotshd flatwav_C.shd

运行 Bellhop,我们得到名为 flatwav_C.shd 的二进制文件,它实际上包含了经过相干计算的声压。我们用 M-文件 plotshd.m 来绘制传播损失图。

③、执行结果

在这里插入图片描述

海底水平的深海波导中的相干传播损失

4、到达接收器处的脉冲响应

修改 flatwav_C.env 设置 OPTIONS3(1)=‘A’

①、环境文件

flatwav_A.env

'Munk profile/Flat waveguide/Arrive'
50.0
1
'SVW'
51 0.0 5000.0
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 /
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.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.8 .0 .0 /
1
1000.0 /
1
1000.0 /
1
101.0 /
'A'
101
-14.0 14.0 /
100.0 5500.0 102.0

②、Matlab 命令

bellhop flatwav_A   %   Runtype = 'A'
% [ Arr, Pos ] = read_arrivals_asc( ARRFile, Narrmx )
[ Arr, Pos ] = read_arrivals_asc( 'flatwav_A.arr' );
% plotarr( filename, irr, ird, isd )
plotarr(  'flatwav_A.arr', 1, 1, 1 )

运行 Bellhop 之后,得到名为 flatwav_A.arr 的 ascii 码文件,其中包含到达接收器位置的射线的振幅和传播时间(我们仅指一个位置点,该模型能够计算到达阵列块中所有阵列点的声线的传播时间和振幅)。包含在 *.arr 文件中的数据可以使用 M-文件 read_arrivals_asc.m 读取。

③、执行结果

在这里插入图片描述

到达接收器处的脉冲响应

三、高斯海山(变化的海底)

本小节描述针对不平海底应用 Bellhop 计算声线的算例。应用高斯函数生成一个理想化的海山, 并写入文件 seamount.bty。

1、海底形状文件

seamount.bty

'L'					% 插值类型
101					% 点数
0 4997.16 			% r( ) z( )
1.01 4997.15 
2.02 4997.15 
3.03 4997.14 
4.04 4997.12 
5.05 4997.1 
6.06 4997.07 
7.07 4997.01 
8.08 4996.93 
9.09 4996.81 
10.1 4996.64 
11.11 4996.38 
12.12 4996.01 
13.13 4995.48 
14.14 4994.74 
15.15 4993.7 
16.16 4992.26 
17.17 4990.3 
18.18 4987.66 
19.19 4984.12 
20.2 4979.45 
21.21 4973.35 
22.22 4965.48 
23.23 4955.43 
24.24 4942.75 
25.25 4926.94 
26.26 4907.46 
27.27 4883.73 
28.28 4855.18 
29.29 4821.26 
30.3 4781.43 
31.31 4735.28 
32.32 4682.48 
33.33 4622.87 
34.34 4556.49 
35.35 4483.61 
36.36 4404.74 
37.37 4320.69 
38.38 4232.57 
39.39 4141.73 
40.4 4049.81 
41.41 3958.67 
42.42 3870.3 
43.43 3786.84 
44.44 3710.38 
45.45 3642.97 
46.46 3586.5 
47.47 3542.59 
48.48 3512.54 
49.49 3497.24 
50.5 3497.16 
51.51 3512.31 
52.52 3542.22 
53.53 3586 
54.54 3642.36 
55.55 3709.66 
56.56 3786.04 
57.57 3869.45 
58.58 3957.77 
59.59 4048.9 
60.6 4140.82 
61.61 4231.68 
62.62 4319.84 
63.63 4403.93 
64.64 4482.85 
65.65 4555.8 
66.66 4622.25 
67.67 4681.93 
68.68 4734.79 
69.69 4781.01 
70.7 4820.89 
71.71 4854.87 
72.72 4883.47 
73.73 4907.24 
74.74 4926.77 
75.75 4942.61 
76.76 4955.32 
77.77 4965.39 
78.78 4973.28 
79.79 4979.4 
80.8 4984.08 
81.81 4987.63 
82.82 4990.28 
83.83 4992.25 
84.84 4993.69 
85.85 4994.73 
86.86 4995.48 
87.87 4996.01 
88.88 4996.38 
89.89 4996.64 
90.9 4996.81 
91.91 4996.93 
92.92 4997.01 
93.93 4997.07 
94.94 4997.1 
95.95 4997.12 
96.96 4997.14 
97.97 4997.15 
98.98 4997.15 
99.99 4997.16 
101 4997.16 

2、环境文件

seamount_R.env

'Munk profile/Sea Mountain/Ray trace'		% TITLE
50.0										% FREQ (Hz)
1											% NMEDIA
'SVW'										% SSP-TOP-WATER-OPT
51 0.0 5000.0								% NMESH SIGMA Z(NSSP)
0.0 1548.52 /								% Z() CP() CS() RHO() AP() AS()
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 /
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.39 /
4400.0 1541.76 /
4600.0 1545.14 /
4800.0 1548.52 /
5000.0 1551.91 /
'A*' 0.0									% BOTOPT SIGMA
5000.0 1600.00 0.0 1.8 .0 0.0				% ZB CPB CSB RHOB APB ASB
1 											% NSD
1000.0 /									% SD(1:NSD) (m)
1 											% NRD
1000.0 /									% RD(1:NRD) (m)
1 											% NRR
101.0 /										% RR(1:NRR ) (km)
'R'											% OPTION: 'R/E/C/A/I/S'
71 											% NBEAMS ISINGLE
-14.0 14.0 /								% ALPHA(1:NBEAMS) (°)
100.0 5500.0 102.0							% STEP (m) ZBOX (m) RBOX (km)

3、Matlab 代码

下面是 Matlab 代码,分别绘制了高斯海山深海波导中的声线轨迹高斯海山的深海波导中的本征声线高斯海山的深海波导中的相干传播损失

clc; clear all;  close all;
global units;	units = 'km';

a=5;    sigma=1; 
x=linspace(0,10.1,101);
y=(1/((sqrt(2*pi))*sigma))*exp(-((x-a).^2)/(2*sigma.^2));
y = 4997.1624 - y / max(y) * 1500; 

fid = fopen('seamount.bty','wt');
fprintf(fid,'%1s%1s%1s\n',char(39),'L',char(39));
fprintf(fid,'%3d\n',length(y));
for mi = 1 : length(y)
    fprintf(fid,'%g %g \n',x(mi)*10,y(mi));
end
fclose(fid);

subplot(321);  bellhop('seamount_R');
plotray('seamount_R');ylim([0 5000])
hold on;    grid on; 
plot(x*1e1,y,'b','LineWidth',1.5); 

subplot(323);  bellhop('seamount_E');
plotray('seamount_E');ylim([0 5000])
hold on;    grid on; 
plot(x*1e1,y,'b','LineWidth',1.5); 

subplot(325);  bellhop('seamount_C');
plotshd('seamount_C.shd');ylim([0 5000])
hold on;    
plot(x*1e1,y,'y','LineWidth',1.5); 

subplot(322);  bellhop('seamount_R');
plotray('seamount_R');ylim([0 5000])
plotbty('seamount_R');  grid on; 

subplot(324);  bellhop('seamount_E');
plotray('seamount_E');ylim([0 5000])
plotbty('seamount_E');  grid on; 

subplot(326);  bellhop('seamount_C');
plotshd('seamount_C.shd');ylim([0 5000])
plotbty('seamount_C');

仿照之前的例子,将 OPTIONS3(1) = 'R’改为 OPTIONS3(1) = 'E’和 OPTIONS3(1) =‘C’,就能够分别计算得到本征声线和相干传播损失,因此这里其他环境文件不再一一列举。

3、执行结果

在这里插入图片描述
上图从上到下,从左到右依次为:高斯海山的深海波导中的声线轨迹高斯海山的深海波导中的本征声线高斯海山的深海波导中的相干传播损失
注:左图和右图的区别在于右边图使用已有的 .bty 文件进行绘制

四、变化的边界(变化的海底和海面)

Bellhop 不仅能够处理变化的海底,还能够同时处理海面、海底都起伏变化的场景。

1、海面形状文件

波浪海面形状文件 varbounds.ati

'L'							% 插值类型
101							% 点数
0 41.2215 					% r( ) z( )
1.01 69.0983 
2.02 100 
3.03 130.902 
4.04 158.779 
5.05 180.902 
6.06 195.106 
7.07 200 
8.08 195.106 
9.09 180.902 
10.1 158.779 
11.11 130.902 
12.12 100 
13.13 69.0983 
14.14 41.2215 
15.15 19.0983 
16.16 4.89435 
17.17 0 
18.18 4.89435 
19.19 19.0983 
20.2 41.2215 
21.21 69.0983 
22.22 100 
23.23 130.902 
24.24 158.779 
25.25 180.902 
26.26 195.106 
27.27 200 
28.28 195.106 
29.29 180.902 
30.3 158.779 
31.31 130.902 
32.32 100 
33.33 69.0983 
34.34 41.2215 
35.35 19.0983 
36.36 4.89435 
37.37 0 
38.38 4.89435 
39.39 19.0983 
40.4 41.2215 
41.41 69.0983 
42.42 100 
43.43 130.902 
44.44 158.779 
45.45 180.902 
46.46 195.106 
47.47 200 
48.48 195.106 
49.49 180.902 
50.5 158.779 
51.51 130.902 
52.52 100 
53.53 69.0983 
54.54 41.2215 
55.55 19.0983 
56.56 4.89435 
57.57 0 
58.58 4.89435 
59.59 19.0983 
60.6 41.2215 
61.61 69.0983 
62.62 100 
63.63 130.902 
64.64 158.779 
65.65 180.902 
66.66 195.106 
67.67 200 
68.68 195.106 
69.69 180.902 
70.7 158.779 
71.71 130.902 
72.72 100 
73.73 69.0983 
74.74 41.2215 
75.75 19.0983 
76.76 4.89435 
77.77 0 
78.78 4.89435 
79.79 19.0983 
80.8 41.2215 
81.81 69.0983 
82.82 100 
83.83 130.902 
84.84 158.779 
85.85 180.902 
86.86 195.106 
87.87 200 
88.88 195.106 
89.89 180.902 
90.9 158.779 
91.91 130.902 
92.92 100 
93.93 69.0983 
94.94 41.2215 
95.95 19.0983 
96.96 4.89435 
97.97 0 
98.98 4.89435 
99.99 19.0983 
101 41.2215 

2、海底形状文件

varbounds_R.bty

'L'
101
0 4997.16 
1.01 4997.15 
2.02 4997.15 
3.03 4997.14 
4.04 4997.12 
5.05 4997.1 
6.06 4997.07 
7.07 4997.01 
8.08 4996.93 
9.09 4996.81 
10.1 4996.64 
11.11 4996.38 
12.12 4996.01 
13.13 4995.48 
14.14 4994.74 
15.15 4993.7 
16.16 4992.26 
17.17 4990.3 
18.18 4987.66 
19.19 4984.12 
20.2 4979.45 
21.21 4973.35 
22.22 4965.48 
23.23 4955.43 
24.24 4942.75 
25.25 4926.94 
26.26 4907.46 
27.27 4883.73 
28.28 4855.18 
29.29 4821.26 
30.3 4781.43 
31.31 4735.28 
32.32 4682.48 
33.33 4622.87 
34.34 4556.49 
35.35 4483.61 
36.36 4404.74 
37.37 4320.69 
38.38 4232.57 
39.39 4141.73 
40.4 4049.81 
41.41 3958.67 
42.42 3870.3 
43.43 3786.84 
44.44 3710.38 
45.45 3642.97 
46.46 3586.5 
47.47 3542.59 
48.48 3512.54 
49.49 3497.24 
50.5 3497.16 
51.51 3512.31 
52.52 3542.22 
53.53 3586 
54.54 3642.36 
55.55 3709.66 
56.56 3786.04 
57.57 3869.45 
58.58 3957.77 
59.59 4048.9 
60.6 4140.82 
61.61 4231.68 
62.62 4319.84 
63.63 4403.93 
64.64 4482.85 
65.65 4555.8 
66.66 4622.25 
67.67 4681.93 
68.68 4734.79 
69.69 4781.01 
70.7 4820.89 
71.71 4854.87 
72.72 4883.47 
73.73 4907.24 
74.74 4926.77 
75.75 4942.61 
76.76 4955.32 
77.77 4965.39 
78.78 4973.28 
79.79 4979.4 
80.8 4984.08 
81.81 4987.63 
82.82 4990.28 
83.83 4992.25 
84.84 4993.69 
85.85 4994.73 
86.86 4995.48 
87.87 4996.01 
88.88 4996.38 
89.89 4996.64 
90.9 4996.81 
91.91 4996.93 
92.92 4997.01 
93.93 4997.07 
94.94 4997.1 
95.95 4997.12 
96.96 4997.14 
97.97 4997.15 
98.98 4997.15 
99.99 4997.16 
101 4997.16 

3、环境文件

varbounds_R.env

'Munk profile/Variable boundaries/Ray trace'	% TITLE
50.0											% FREQ (Hz)
1												% NMEDIA
'SVW *'											% SSP-TOP-WATER-OPT
51 0.0 5000.0									% NMESH SIGMA Z(NSSP)
0.0 1548.52 /									% Z() CP() CS() RHO() AP() AS()
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 /
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.39 /
4400.0 1541.76 /
4600.0 1545.14 /
4800.0 1548.52 /
5000.0 1551.91 /
'A*' 0.0										% BOTOPT SIGMA
5000.0 1600.00 0.0 1.8 .0 0.0					% ZB CPB CSB RHOB APB ASB
1 												% NSD
1000.0 /										% SD(1:NSD) (m)
1												% NRD
1000.0 /										% RD(1:NRD) (m)
1												% NRR
101.0 /											% RR(1:NRR ) (km)
'R'												% OPTION: 'R/E/C/A/I/S'
71 												% NBEAMS ISINGLE
-14.0 14.0 /									% ALPHA(1:NBEAMS) (°)
100.0 5500.0 102.0								% STEP (m) ZBOX (m) RBOX (km)

4、Matlab 代码

下面是 Matlab 代码,分别绘制了波浪海面和高斯海山的深海波导中声速剖面声线轨迹本征声线传播损失

clc;    clear all;  close all;
global units;	units = 'km';
%=======  Sea Surface  ============
xs = linspace(0,10*pi,101);
ys = 100 + sin(xs - pi/5) * 100;
xs = xs / max(xs) *101;
%=======  Write the Sea Surface file ============
% fid = fopen('varbounds_R.ati','wt');
% fprintf(fid,'%1s%1s%1s\n',char(39),'L',char(39));
% fprintf(fid,'%3d\n',length(ys));
% for mi = 1 : length(ys)
%     fprintf(fid,'%g %g \n',xs(mi),ys(mi));
% end
% fclose(fid);
%=======  Sea Bottom  ============
a = 5;    sigma = 1; 
x = linspace( 0,10.1,101 );
y = (1/((sqrt(2*pi)) * sigma)) * exp(-((x-a).^2)/(2*sigma.^2));
y = 4997.1624 - y / max(y) * 1500; 
%=======  Write the Sea Bottom file ============
% fid = fopen('varbounds_R.bty','wt');
% fprintf(fid,'%1s%1s%1s\n',char(39),'L',char(39));
% fprintf(fid,'%3d\n',length(y));
% for mi = 1 : length(y)
%     fprintf(fid,'%g %g \n',x(mi)*10,y(mi));
% end
% fclose(fid);
%=======  Calculating and Plotting  ============
subplot(3,6,2.5:6); bellhop('varbounds_R');
plotray('varbounds_R');ylim([0 5000])
hold on;    grid on; 
plotati('varbounds_E');plotbty('varbounds_E')
% plot(x*1e4,y,'b','LineWidth',1.5); 
% plot(xs*1e3,ys,'b','LineWidth',1.5);  

subplot(3,6,8.5:12); bellhop('varbounds_E');
plotray('varbounds_E');ylim([0 5000])
hold on;    grid on; 
plotati('varbounds_E');plotbty('varbounds_E');  xlabel('');
% plot(x*1e4,y,'b','LineWidth',1.5); 
% plot(xs*1e3,ys,'b','LineWidth',1.5);  

subplot(3,6,14.5:18); bellhop('varbounds_C');
plotshd('varbounds_C.shd');ylim([0 5000])
hold on;    grid on; 
plotati('varbounds_E');plotbty('varbounds_E')
% plot(x*1e4,y,'y','LineWidth',1.5); 
% plot(xs*1e3,ys,'y','LineWidth',1.5);  

%=======  Other Plottings  ============
subplot(3,6,[1 7 13]);plotssp('varbounds_C.env')

仿照之前的例子,将 OPTIONS3(1) = 'C’改为 OPTIONS3(1) = 'E’和 OPTIONS3(1) =‘R’,就能够分别计算得到本征声线和相干传播损失,因此这里其他环境文件不再一一列举。

5、执行结果

在这里插入图片描述
从左到右,从上到下分别绘制了波浪海面和高斯海山的深海波导中声速剖面声线轨迹本征声线传播损失

五、分段线性边界(Dickins 海山)

Dickins 海山场景:先将海底测深理想化为 3000m 深度的平面,在距离 20km 的地方,升起一座海山,山顶延伸至海深 500m 处,海山模拟成三角形,底边 20km。

1、海底形状文件

DickinsB.bty

'L'
5
  0 3000
 10 3000
 20  500
 30 3000
100 3000
  • 第 1 行指明“分段线性拟合”设为“L”,“曲线拟合”设为“C”。因为我们打算模拟一座三角形的海山,所以选择“分段线性拟合”。
  • 第 2 行是测深点的数目,后续各行是以“距离-深度”来定义的“海底测深”

2、环境文件

DickinsB.env

'Dickins seamount'	! TITLE
230.0			! FREQ (Hz)
1			! NMEDIA
'CVW'			! SSPOPT (Analytic 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.7 /
1500 1483.8 /
2000 1490.5 /
2500 1498.3 /
3000 1506.5 /
'A*' 0.0
3000.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/S'
0		       		! NBEAMS
-89.0 89.0 /		        ! ALPHA1,2 (degrees)
0.0  3100.0  101.0		! STEP (m), ZBOX (m), RBOX (km)
  • 第 4 行:我们看到顶端选项(针对声速剖面、海面以上半空间、海洋水体的设定)为“CVW”。“C”表示我们希望以“分段线性”方式对声速 c 进行插值。“V”是将海洋表面设置为“标准真空”的选项。“W”表示将衰减单位设定
    为 dB/λ(波长),这是用于环境描述的单位。
  • 第 5 行:水深设定为 3000m。在这种海底测深随距离变化的场景中,该深度值不是直接使用的,不过,因为它是海洋声速剖面定义的一部分,因此该深度值必须深到足以覆盖测深文件中随距离而变化的最深深度。
  • 第 36 行中:我们看到“海底选项”设置为“A*”。“A”表示海底是“声学弹性”半空间,这是典型的设置。这里引入“*”,作为标志,是告诉 BELLHOP 要读取一个附加的测深文件“Dickins.bty”,来标识海深是随距离变化的。

3、Matlab 代码

clc;    clear all; % close all
global units ;  units = 'km';

bellhop( 'DickinsB' )
plotshd( 'DickinsB.shd', 2, 1, 1 )
caxisrev( [ 70 120 ] )
plotbty 'DickinsB'	% 绘制海深曲线

ram
plotshd( 'RAM.shd.mat', 2, 1, 2 )
caxisrev( [ 70 120 ] )
plotbty 'DickinsB'

% bellhop DickinsB_oneBeam
% plotshd  DickinsB_oneBeam.shd
% caxisrev( [ 70 120 ] )
% plotbty 'DickinsB'
TL = load('tl.line');
plot(TL(:,1)/1000,TL(:,2));
axis ij;

4、执行结果

在这里插入图片描述

Dickins 海山场景下的传播损失,从上到下依次为应用 Bellhop 、 应用RAM(抛物线方程模型)

上半图板是 BELLHOP 的计算结果,下半图板是我们当做参考解的 RAM PE 的计算结果。

两者的一致性是令人满意的;不过,人为的海山尖顶导致了大量的能量衍射。通过在不连续的测深点附近插入额外的测深点,这种情况可以得到进一步改善。

六、曲线拟合边界(抛物线海底)

我们考虑了 McGirr 等人所描述的抛物线型海底剖面,其深度由 D ( r ) = 500 1 + 4 r D(r) = 500 \sqrt{1+4r} D(r)=5001+4r 给出,其中距离 r 单位为公里,深度 D 单位为米。声源位于焦点并作为原点,这是最理想的情况。

1、海底形状文件

ParaBot.bty

'C' 
1001 
 -0.250000 0.000000 
 -0.249975 5.000000 
 -0.249900 10.000000 
 -0.249775 15.000000 
 -0.249600 20.000000 
 -0.249375 25.000000 
 -0.249100 30.000000 
 -0.248775 35.000000 
 -0.248400 40.000000 
 -0.247975 45.000000 
 -0.247500 50.000000 
 -0.246975 55.000000 
 -0.246400 60.000000 
 -0.245775 65.000000 
 -0.245100 70.000000 
 -0.244375 75.000000 
 -0.243600 80.000000 
 -0.242775 85.000000 
 -0.241900 90.000000 
 -0.240975 95.000000 
 -0.240000 100.000000 
 -0.238975 105.000000 
 -0.237900 110.000000 
 -0.236775 115.000000 
 -0.235600 120.000000 
 -0.234375 125.000000 
 -0.233100 130.000000 
 -0.231775 135.000000 
 -0.230400 140.000000 
 -0.228975 145.000000 
 -0.227500 150.000000 
 -0.225975 155.000000 
 -0.224400 160.000000 
 -0.222775 165.000000 
 -0.221100 170.000000 
 -0.219375 175.000000 
 -0.217600 180.000000 
 -0.215775 185.000000 
 -0.213900 190.000000 
 -0.211975 195.000000 
 -0.210000 200.000000 
 -0.207975 205.000000 
 -0.205900 210.000000 
 -0.203775 215.000000 
 -0.201600 220.000000 
 -0.199375 225.000000 
 -0.197100 230.000000 
 -0.194775 235.000000 
 -0.192400 240.000000 
 -0.189975 245.000000 
 -0.187500 250.000000 
 -0.184975 255.000000 
 -0.182400 260.000000 
 -0.179775 265.000000 
 -0.177100 270.000000 
 -0.174375 275.000000 
 -0.171600 280.000000 
 -0.168775 285.000000 
 -0.165900 290.000000 
 -0.162975 295.000000 
 -0.160000 300.000000 
 -0.156975 305.000000 
 -0.153900 310.000000 
 -0.150775 315.000000 
 -0.147600 320.000000 
 -0.144375 325.000000 
 -0.141100 330.000000 
 -0.137775 335.000000 
 -0.134400 340.000000 
 -0.130975 345.000000 
 -0.127500 350.000000 
 -0.123975 355.000000 
 -0.120400 360.000000 
 -0.116775 365.000000 
 -0.113100 370.000000 
 -0.109375 375.000000 
 -0.105600 380.000000 
 -0.101775 385.000000 
 -0.097900 390.000000 
 -0.093975 395.000000 
 -0.090000 400.000000 
 -0.085975 405.000000 
 -0.081900 410.000000 
 -0.077775 415.000000 
 -0.073600 420.000000 
 -0.069375 425.000000 
 -0.065100 430.000000 
 -0.060775 435.000000 
 -0.056400 440.000000 
 -0.051975 445.000000 
 -0.047500 450.000000 
 -0.042975 455.000000 
 -0.038400 460.000000 
 -0.033775 465.000000 
 -0.029100 470.000000 
 -0.024375 475.000000 
 -0.019600 480.000000 
 -0.014775 485.000000 
 -0.009900 490.000000 
 -0.004975 495.000000 
 0.000000 500.000000 
 0.005025 505.000000 
 0.010100 510.000000 
 0.015225 515.000000 
 0.020400 520.000000 
 0.025625 525.000000 
 0.030900 530.000000 
 0.036225 535.000000 
 0.041600 540.000000 
 0.047025 545.000000 
 0.052500 550.000000 
 0.058025 555.000000 
 0.063600 560.000000 
 0.069225 565.000000 
 0.074900 570.000000 
 0.080625 575.000000 
 0.086400 580.000000 
 0.092225 585.000000 
 0.098100 590.000000 
 0.104025 595.000000 
 0.110000 600.000000 
 0.116025 605.000000 
 0.122100 610.000000 
 0.128225 615.000000 
 0.134400 620.000000 
 0.140625 625.000000 
 0.146900 630.000000 
 0.153225 635.000000 
 0.159600 640.000000 
 0.166025 645.000000 
 0.172500 650.000000 
 0.179025 655.000000 
 0.185600 660.000000 
 0.192225 665.000000 
 0.198900 670.000000 
 0.205625 675.000000 
 0.212400 680.000000 
 0.219225 685.000000 
 0.226100 690.000000 
 0.233025 695.000000 
 0.240000 700.000000 
 0.247025 705.000000 
 0.254100 710.000000 
 0.261225 715.000000 
 0.268400 720.000000 
 0.275625 725.000000 
 0.282900 730.000000 
 0.290225 735.000000 
 0.297600 740.000000 
 0.305025 745.000000 
 0.312500 750.000000 
 0.320025 755.000000 
 0.327600 760.000000 
 0.335225 765.000000 
 0.342900 770.000000 
 0.350625 775.000000 
 0.358400 780.000000 
 0.366225 785.000000 
 0.374100 790.000000 
 0.382025 795.000000 
 0.390000 800.000000 
 0.398025 805.000000 
 0.406100 810.000000 
 0.414225 815.000000 
 0.422400 820.000000 
 0.430625 825.000000 
 0.438900 830.000000 
 0.447225 835.000000 
 0.455600 840.000000 
 0.464025 845.000000 
 0.472500 850.000000 
 0.481025 855.000000 
 0.489600 860.000000 
 0.498225 865.000000 
 0.506900 870.000000 
 0.515625 875.000000 
 0.524400 880.000000 
 0.533225 885.000000 
 0.542100 890.000000 
 0.551025 895.000000 
 0.560000 900.000000 
 0.569025 905.000000 
 0.578100 910.000000 
 0.587225 915.000000 
 0.596400 920.000000 
 0.605625 925.000000 
 0.614900 930.000000 
 0.624225 935.000000 
 0.633600 940.000000 
 0.643025 945.000000 
 0.652500 950.000000 
 0.662025 955.000000 
 0.671600 960.000000 
 0.681225 965.000000 
 0.690900 970.000000 
 0.700625 975.000000 
 0.710400 980.000000 
 0.720225 985.000000 
 0.730100 990.000000 
 0.740025 995.000000 
 0.750000 1000.000000 
 0.760025 1005.000000 
 0.770100 1010.000000 
 0.780225 1015.000000 
 0.790400 1020.000000 
 0.800625 1025.000000 
 0.810900 1030.000000 
 0.821225 1035.000000 
 0.831600 1040.000000 
 0.842025 1045.000000 
 0.852500 1050.000000 
 0.863025 1055.000000 
 0.873600 1060.000000 
 0.884225 1065.000000 
 0.894900 1070.000000 
 0.905625 1075.000000 
 0.916400 1080.000000 
 0.927225 1085.000000 
 0.938100 1090.000000 
 0.949025 1095.000000 
 0.960000 1100.000000 
 0.971025 1105.000000 
 0.982100 1110.000000 
 0.993225 1115.000000 
 1.004400 1120.000000 
 1.015625 1125.000000 
 1.026900 1130.000000 
 1.038225 1135.000000 
 1.049600 1140.000000 
 1.061025 1145.000000 
 1.072500 1150.000000 
 1.084025 1155.000000 
 1.095600 1160.000000 
 1.107225 1165.000000 
 1.118900 1170.000000 
 1.130625 1175.000000 
 1.142400 1180.000000 
 1.154225 1185.000000 
 1.166100 1190.000000 
 1.178025 1195.000000 
 1.190000 1200.000000 
 1.202025 1205.000000 
 1.214100 1210.000000 
 1.226225 1215.000000 
 1.238400 1220.000000 
 1.250625 1225.000000 
 1.262900 1230.000000 
 1.275225 1235.000000 
 1.287600 1240.000000 
 1.300025 1245.000000 
 1.312500 1250.000000 
 1.325025 1255.000000 
 1.337600 1260.000000 
 1.350225 1265.000000 
 1.362900 1270.000000 
 1.375625 1275.000000 
 1.388400 1280.000000 
 1.401225 1285.000000 
 1.414100 1290.000000 
 1.427025 1295.000000 
 1.440000 1300.000000 
 1.453025 1305.000000 
 1.466100 1310.000000 
 1.479225 1315.000000 
 1.492400 1320.000000 
 1.505625 1325.000000 
 1.518900 1330.000000 
 1.532225 1335.000000 
 1.545600 1340.000000 
 1.559025 1345.000000 
 1.572500 1350.000000 
 1.586025 1355.000000 
 1.599600 1360.000000 
 1.613225 1365.000000 
 1.626900 1370.000000 
 1.640625 1375.000000 
 1.654400 1380.000000 
 1.668225 1385.000000 
 1.682100 1390.000000 
 1.696025 1395.000000 
 1.710000 1400.000000 
 1.724025 1405.000000 
 1.738100 1410.000000 
 1.752225 1415.000000 
 1.766400 1420.000000 
 1.780625 1425.000000 
 1.794900 1430.000000 
 1.809225 1435.000000 
 1.823600 1440.000000 
 1.838025 1445.000000 
 1.852500 1450.000000 
 1.867025 1455.000000 
 1.881600 1460.000000 
 1.896225 1465.000000 
 1.910900 1470.000000 
 1.925625 1475.000000 
 1.940400 1480.000000 
 1.955225 1485.000000 
 1.970100 1490.000000 
 1.985025 1495.000000 
 2.000000 1500.000000 
 2.015025 1505.000000 
 2.030100 1510.000000 
 2.045225 1515.000000 
 2.060400 1520.000000 
 2.075625 1525.000000 
 2.090900 1530.000000 
 2.106225 1535.000000 
 2.121600 1540.000000 
 2.137025 1545.000000 
 2.152500 1550.000000 
 2.168025 1555.000000 
 2.183600 1560.000000 
 2.199225 1565.000000 
 2.214900 1570.000000 
 2.230625 1575.000000 
 2.246400 1580.000000 
 2.262225 1585.000000 
 2.278100 1590.000000 
 2.294025 1595.000000 
 2.310000 1600.000000 
 2.326025 1605.000000 
 2.342100 1610.000000 
 2.358225 1615.000000 
 2.374400 1620.000000 
 2.390625 1625.000000 
 2.406900 1630.000000 
 2.423225 1635.000000 
 2.439600 1640.000000 
 2.456025 1645.000000 
 2.472500 1650.000000 
 2.489025 1655.000000 
 2.505600 1660.000000 
 2.522225 1665.000000 
 2.538900 1670.000000 
 2.555625 1675.000000 
 2.572400 1680.000000 
 2.589225 1685.000000 
 2.606100 1690.000000 
 2.623025 1695.000000 
 2.640000 1700.000000 
 2.657025 1705.000000 
 2.674100 1710.000000 
 2.691225 1715.000000 
 2.708400 1720.000000 
 2.725625 1725.000000 
 2.742900 1730.000000 
 2.760225 1735.000000 
 2.777600 1740.000000 
 2.795025 1745.000000 
 2.812500 1750.000000 
 2.830025 1755.000000 
 2.847600 1760.000000 
 2.865225 1765.000000 
 2.882900 1770.000000 
 2.900625 1775.000000 
 2.918400 1780.000000 
 2.936225 1785.000000 
 2.954100 1790.000000 
 2.972025 1795.000000 
 2.990000 1800.000000 
 3.008025 1805.000000 
 3.026100 1810.000000 
 3.044225 1815.000000 
 3.062400 1820.000000 
 3.080625 1825.000000 
 3.098900 1830.000000 
 3.117225 1835.000000 
 3.135600 1840.000000 
 3.154025 1845.000000 
 3.172500 1850.000000 
 3.191025 1855.000000 
 3.209600 1860.000000 
 3.228225 1865.000000 
 3.246900 1870.000000 
 3.265625 1875.000000 
 3.284400 1880.000000 
 3.303225 1885.000000 
 3.322100 1890.000000 
 3.341025 1895.000000 
 3.360000 1900.000000 
 3.379025 1905.000000 
 3.398100 1910.000000 
 3.417225 1915.000000 
 3.436400 1920.000000 
 3.455625 1925.000000 
 3.474900 1930.000000 
 3.494225 1935.000000 
 3.513600 1940.000000 
 3.533025 1945.000000 
 3.552500 1950.000000 
 3.572025 1955.000000 
 3.591600 1960.000000 
 3.611225 1965.000000 
 3.630900 1970.000000 
 3.650625 1975.000000 
 3.670400 1980.000000 
 3.690225 1985.000000 
 3.710100 1990.000000 
 3.730025 1995.000000 
 3.750000 2000.000000 
 3.770025 2005.000000 
 3.790100 2010.000000 
 3.810225 2015.000000 
 3.830400 2020.000000 
 3.850625 2025.000000 
 3.870900 2030.000000 
 3.891225 2035.000000 
 3.911600 2040.000000 
 3.932025 2045.000000 
 3.952500 2050.000000 
 3.973025 2055.000000 
 3.993600 2060.000000 
 4.014225 2065.000000 
 4.034900 2070.000000 
 4.055625 2075.000000 
 4.076400 2080.000000 
 4.097225 2085.000000 
 4.118100 2090.000000 
 4.139025 2095.000000 
 4.160000 2100.000000 
 4.181025 2105.000000 
 4.202100 2110.000000 
 4.223225 2115.000000 
 4.244400 2120.000000 
 4.265625 2125.000000 
 4.286900 2130.000000 
 4.308225 2135.000000 
 4.329600 2140.000000 
 4.351025 2145.000000 
 4.372500 2150.000000 
 4.394025 2155.000000 
 4.415600 2160.000000 
 4.437225 2165.000000 
 4.458900 2170.000000 
 4.480625 2175.000000 
 4.502400 2180.000000 
 4.524225 2185.000000 
 4.546100 2190.000000 
 4.568025 2195.000000 
 4.590000 2200.000000 
 4.612025 2205.000000 
 4.634100 2210.000000 
 4.656225 2215.000000 
 4.678400 2220.000000 
 4.700625 2225.000000 
 4.722900 2230.000000 
 4.745225 2235.000000 
 4.767600 2240.000000 
 4.790025 2245.000000 
 4.812500 2250.000000 
 4.835025 2255.000000 
 4.857600 2260.000000 
 4.880225 2265.000000 
 4.902900 2270.000000 
 4.925625 2275.000000 
 4.948400 2280.000000 
 4.971225 2285.000000 
 4.994100 2290.000000 
 5.017025 2295.000000 
 5.040000 2300.000000 
 5.063025 2305.000000 
 5.086100 2310.000000 
 5.109225 2315.000000 
 5.132400 2320.000000 
 5.155625 2325.000000 
 5.178900 2330.000000 
 5.202225 2335.000000 
 5.225600 2340.000000 
 5.249025 2345.000000 
 5.272500 2350.000000 
 5.296025 2355.000000 
 5.319600 2360.000000 
 5.343225 2365.000000 
 5.366900 2370.000000 
 5.390625 2375.000000 
 5.414400 2380.000000 
 5.438225 2385.000000 
 5.462100 2390.000000 
 5.486025 2395.000000 
 5.510000 2400.000000 
 5.534025 2405.000000 
 5.558100 2410.000000 
 5.582225 2415.000000 
 5.606400 2420.000000 
 5.630625 2425.000000 
 5.654900 2430.000000 
 5.679225 2435.000000 
 5.703600 2440.000000 
 5.728025 2445.000000 
 5.752500 2450.000000 
 5.777025 2455.000000 
 5.801600 2460.000000 
 5.826225 2465.000000 
 5.850900 2470.000000 
 5.875625 2475.000000 
 5.900400 2480.000000 
 5.925225 2485.000000 
 5.950100 2490.000000 
 5.975025 2495.000000 
 6.000000 2500.000000 
 6.025025 2505.000000 
 6.050100 2510.000000 
 6.075225 2515.000000 
 6.100400 2520.000000 
 6.125625 2525.000000 
 6.150900 2530.000000 
 6.176225 2535.000000 
 6.201600 2540.000000 
 6.227025 2545.000000 
 6.252500 2550.000000 
 6.278025 2555.000000 
 6.303600 2560.000000 
 6.329225 2565.000000 
 6.354900 2570.000000 
 6.380625 2575.000000 
 6.406400 2580.000000 
 6.432225 2585.000000 
 6.458100 2590.000000 
 6.484025 2595.000000 
 6.510000 2600.000000 
 6.536025 2605.000000 
 6.562100 2610.000000 
 6.588225 2615.000000 
 6.614400 2620.000000 
 6.640625 2625.000000 
 6.666900 2630.000000 
 6.693225 2635.000000 
 6.719600 2640.000000 
 6.746025 2645.000000 
 6.772500 2650.000000 
 6.799025 2655.000000 
 6.825600 2660.000000 
 6.852225 2665.000000 
 6.878900 2670.000000 
 6.905625 2675.000000 
 6.932400 2680.000000 
 6.959225 2685.000000 
 6.986100 2690.000000 
 7.013025 2695.000000 
 7.040000 2700.000000 
 7.067025 2705.000000 
 7.094100 2710.000000 
 7.121225 2715.000000 
 7.148400 2720.000000 
 7.175625 2725.000000 
 7.202900 2730.000000 
 7.230225 2735.000000 
 7.257600 2740.000000 
 7.285025 2745.000000 
 7.312500 2750.000000 
 7.340025 2755.000000 
 7.367600 2760.000000 
 7.395225 2765.000000 
 7.422900 2770.000000 
 7.450625 2775.000000 
 7.478400 2780.000000 
 7.506225 2785.000000 
 7.534100 2790.000000 
 7.562025 2795.000000 
 7.590000 2800.000000 
 7.618025 2805.000000 
 7.646100 2810.000000 
 7.674225 2815.000000 
 7.702400 2820.000000 
 7.730625 2825.000000 
 7.758900 2830.000000 
 7.787225 2835.000000 
 7.815600 2840.000000 
 7.844025 2845.000000 
 7.872500 2850.000000 
 7.901025 2855.000000 
 7.929600 2860.000000 
 7.958225 2865.000000 
 7.986900 2870.000000 
 8.015625 2875.000000 
 8.044400 2880.000000 
 8.073225 2885.000000 
 8.102100 2890.000000 
 8.131025 2895.000000 
 8.160000 2900.000000 
 8.189025 2905.000000 
 8.218100 2910.000000 
 8.247225 2915.000000 
 8.276400 2920.000000 
 8.305625 2925.000000 
 8.334900 2930.000000 
 8.364225 2935.000000 
 8.393600 2940.000000 
 8.423025 2945.000000 
 8.452500 2950.000000 
 8.482025 2955.000000 
 8.511600 2960.000000 
 8.541225 2965.000000 
 8.570900 2970.000000 
 8.600625 2975.000000 
 8.630400 2980.000000 
 8.660225 2985.000000 
 8.690100 2990.000000 
 8.720025 2995.000000 
 8.750000 3000.000000 
 8.780025 3005.000000 
 8.810100 3010.000000 
 8.840225 3015.000000 
 8.870400 3020.000000 
 8.900625 3025.000000 
 8.930900 3030.000000 
 8.961225 3035.000000 
 8.991600 3040.000000 
 9.022025 3045.000000 
 9.052500 3050.000000 
 9.083025 3055.000000 
 9.113600 3060.000000 
 9.144225 3065.000000 
 9.174900 3070.000000 
 9.205625 3075.000000 
 9.236400 3080.000000 
 9.267225 3085.000000 
 9.298100 3090.000000 
 9.329025 3095.000000 
 9.360000 3100.000000 
 9.391025 3105.000000 
 9.422100 3110.000000 
 9.453225 3115.000000 
 9.484400 3120.000000 
 9.515625 3125.000000 
 9.546900 3130.000000 
 9.578225 3135.000000 
 9.609600 3140.000000 
 9.641025 3145.000000 
 9.672500 3150.000000 
 9.704025 3155.000000 
 9.735600 3160.000000 
 9.767225 3165.000000 
 9.798900 3170.000000 
 9.830625 3175.000000 
 9.862400 3180.000000 
 9.894225 3185.000000 
 9.926100 3190.000000 
 9.958025 3195.000000 
 9.990000 3200.000000 
 10.022025 3205.000000 
 10.054100 3210.000000 
 10.086225 3215.000000 
 10.118400 3220.000000 
 10.150625 3225.000000 
 10.182900 3230.000000 
 10.215225 3235.000000 
 10.247600 3240.000000 
 10.280025 3245.000000 
 10.312500 3250.000000 
 10.345025 3255.000000 
 10.377600 3260.000000 
 10.410225 3265.000000 
 10.442900 3270.000000 
 10.475625 3275.000000 
 10.508400 3280.000000 
 10.541225 3285.000000 
 10.574100 3290.000000 
 10.607025 3295.000000 
 10.640000 3300.000000 
 10.673025 3305.000000 
 10.706100 3310.000000 
 10.739225 3315.000000 
 10.772400 3320.000000 
 10.805625 3325.000000 
 10.838900 3330.000000 
 10.872225 3335.000000 
 10.905600 3340.000000 
 10.939025 3345.000000 
 10.972500 3350.000000 
 11.006025 3355.000000 
 11.039600 3360.000000 
 11.073225 3365.000000 
 11.106900 3370.000000 
 11.140625 3375.000000 
 11.174400 3380.000000 
 11.208225 3385.000000 
 11.242100 3390.000000 
 11.276025 3395.000000 
 11.310000 3400.000000 
 11.344025 3405.000000 
 11.378100 3410.000000 
 11.412225 3415.000000 
 11.446400 3420.000000 
 11.480625 3425.000000 
 11.514900 3430.000000 
 11.549225 3435.000000 
 11.583600 3440.000000 
 11.618025 3445.000000 
 11.652500 3450.000000 
 11.687025 3455.000000 
 11.721600 3460.000000 
 11.756225 3465.000000 
 11.790900 3470.000000 
 11.825625 3475.000000 
 11.860400 3480.000000 
 11.895225 3485.000000 
 11.930100 3490.000000 
 11.965025 3495.000000 
 12.000000 3500.000000 
 12.035025 3505.000000 
 12.070100 3510.000000 
 12.105225 3515.000000 
 12.140400 3520.000000 
 12.175625 3525.000000 
 12.210900 3530.000000 
 12.246225 3535.000000 
 12.281600 3540.000000 
 12.317025 3545.000000 
 12.352500 3550.000000 
 12.388025 3555.000000 
 12.423600 3560.000000 
 12.459225 3565.000000 
 12.494900 3570.000000 
 12.530625 3575.000000 
 12.566400 3580.000000 
 12.602225 3585.000000 
 12.638100 3590.000000 
 12.674025 3595.000000 
 12.710000 3600.000000 
 12.746025 3605.000000 
 12.782100 3610.000000 
 12.818225 3615.000000 
 12.854400 3620.000000 
 12.890625 3625.000000 
 12.926900 3630.000000 
 12.963225 3635.000000 
 12.999600 3640.000000 
 13.036025 3645.000000 
 13.072500 3650.000000 
 13.109025 3655.000000 
 13.145600 3660.000000 
 13.182225 3665.000000 
 13.218900 3670.000000 
 13.255625 3675.000000 
 13.292400 3680.000000 
 13.329225 3685.000000 
 13.366100 3690.000000 
 13.403025 3695.000000 
 13.440000 3700.000000 
 13.477025 3705.000000 
 13.514100 3710.000000 
 13.551225 3715.000000 
 13.588400 3720.000000 
 13.625625 3725.000000 
 13.662900 3730.000000 
 13.700225 3735.000000 
 13.737600 3740.000000 
 13.775025 3745.000000 
 13.812500 3750.000000 
 13.850025 3755.000000 
 13.887600 3760.000000 
 13.925225 3765.000000 
 13.962900 3770.000000 
 14.000625 3775.000000 
 14.038400 3780.000000 
 14.076225 3785.000000 
 14.114100 3790.000000 
 14.152025 3795.000000 
 14.190000 3800.000000 
 14.228025 3805.000000 
 14.266100 3810.000000 
 14.304225 3815.000000 
 14.342400 3820.000000 
 14.380625 3825.000000 
 14.418900 3830.000000 
 14.457225 3835.000000 
 14.495600 3840.000000 
 14.534025 3845.000000 
 14.572500 3850.000000 
 14.611025 3855.000000 
 14.649600 3860.000000 
 14.688225 3865.000000 
 14.726900 3870.000000 
 14.765625 3875.000000 
 14.804400 3880.000000 
 14.843225 3885.000000 
 14.882100 3890.000000 
 14.921025 3895.000000 
 14.960000 3900.000000 
 14.999025 3905.000000 
 15.038100 3910.000000 
 15.077225 3915.000000 
 15.116400 3920.000000 
 15.155625 3925.000000 
 15.194900 3930.000000 
 15.234225 3935.000000 
 15.273600 3940.000000 
 15.313025 3945.000000 
 15.352500 3950.000000 
 15.392025 3955.000000 
 15.431600 3960.000000 
 15.471225 3965.000000 
 15.510900 3970.000000 
 15.550625 3975.000000 
 15.590400 3980.000000 
 15.630225 3985.000000 
 15.670100 3990.000000 
 15.710025 3995.000000 
 15.750000 4000.000000 
 15.790025 4005.000000 
 15.830100 4010.000000 
 15.870225 4015.000000 
 15.910400 4020.000000 
 15.950625 4025.000000 
 15.990900 4030.000000 
 16.031225 4035.000000 
 16.071600 4040.000000 
 16.112025 4045.000000 
 16.152500 4050.000000 
 16.193025 4055.000000 
 16.233600 4060.000000 
 16.274225 4065.000000 
 16.314900 4070.000000 
 16.355625 4075.000000 
 16.396400 4080.000000 
 16.437225 4085.000000 
 16.478100 4090.000000 
 16.519025 4095.000000 
 16.560000 4100.000000 
 16.601025 4105.000000 
 16.642100 4110.000000 
 16.683225 4115.000000 
 16.724400 4120.000000 
 16.765625 4125.000000 
 16.806900 4130.000000 
 16.848225 4135.000000 
 16.889600 4140.000000 
 16.931025 4145.000000 
 16.972500 4150.000000 
 17.014025 4155.000000 
 17.055600 4160.000000 
 17.097225 4165.000000 
 17.138900 4170.000000 
 17.180625 4175.000000 
 17.222400 4180.000000 
 17.264225 4185.000000 
 17.306100 4190.000000 
 17.348025 4195.000000 
 17.390000 4200.000000 
 17.432025 4205.000000 
 17.474100 4210.000000 
 17.516225 4215.000000 
 17.558400 4220.000000 
 17.600625 4225.000000 
 17.642900 4230.000000 
 17.685225 4235.000000 
 17.727600 4240.000000 
 17.770025 4245.000000 
 17.812500 4250.000000 
 17.855025 4255.000000 
 17.897600 4260.000000 
 17.940225 4265.000000 
 17.982900 4270.000000 
 18.025625 4275.000000 
 18.068400 4280.000000 
 18.111225 4285.000000 
 18.154100 4290.000000 
 18.197025 4295.000000 
 18.240000 4300.000000 
 18.283025 4305.000000 
 18.326100 4310.000000 
 18.369225 4315.000000 
 18.412400 4320.000000 
 18.455625 4325.000000 
 18.498900 4330.000000 
 18.542225 4335.000000 
 18.585600 4340.000000 
 18.629025 4345.000000 
 18.672500 4350.000000 
 18.716025 4355.000000 
 18.759600 4360.000000 
 18.803225 4365.000000 
 18.846900 4370.000000 
 18.890625 4375.000000 
 18.934400 4380.000000 
 18.978225 4385.000000 
 19.022100 4390.000000 
 19.066025 4395.000000 
 19.110000 4400.000000 
 19.154025 4405.000000 
 19.198100 4410.000000 
 19.242225 4415.000000 
 19.286400 4420.000000 
 19.330625 4425.000000 
 19.374900 4430.000000 
 19.419225 4435.000000 
 19.463600 4440.000000 
 19.508025 4445.000000 
 19.552500 4450.000000 
 19.597025 4455.000000 
 19.641600 4460.000000 
 19.686225 4465.000000 
 19.730900 4470.000000 
 19.775625 4475.000000 
 19.820400 4480.000000 
 19.865225 4485.000000 
 19.910100 4490.000000 
 19.955025 4495.000000 
 20.000000 4500.000000 
 20.045025 4505.000000 
 20.090100 4510.000000 
 20.135225 4515.000000 
 20.180400 4520.000000 
 20.225625 4525.000000 
 20.270900 4530.000000 
 20.316225 4535.000000 
 20.361600 4540.000000 
 20.407025 4545.000000 
 20.452500 4550.000000 
 20.498025 4555.000000 
 20.543600 4560.000000 
 20.589225 4565.000000 
 20.634900 4570.000000 
 20.680625 4575.000000 
 20.726400 4580.000000 
 20.772225 4585.000000 
 20.818100 4590.000000 
 20.864025 4595.000000 
 20.910000 4600.000000 
 20.956025 4605.000000 
 21.002100 4610.000000 
 21.048225 4615.000000 
 21.094400 4620.000000 
 21.140625 4625.000000 
 21.186900 4630.000000 
 21.233225 4635.000000 
 21.279600 4640.000000 
 21.326025 4645.000000 
 21.372500 4650.000000 
 21.419025 4655.000000 
 21.465600 4660.000000 
 21.512225 4665.000000 
 21.558900 4670.000000 
 21.605625 4675.000000 
 21.652400 4680.000000 
 21.699225 4685.000000 
 21.746100 4690.000000 
 21.793025 4695.000000 
 21.840000 4700.000000 
 21.887025 4705.000000 
 21.934100 4710.000000 
 21.981225 4715.000000 
 22.028400 4720.000000 
 22.075625 4725.000000 
 22.122900 4730.000000 
 22.170225 4735.000000 
 22.217600 4740.000000 
 22.265025 4745.000000 
 22.312500 4750.000000 
 22.360025 4755.000000 
 22.407600 4760.000000 
 22.455225 4765.000000 
 22.502900 4770.000000 
 22.550625 4775.000000 
 22.598400 4780.000000 
 22.646225 4785.000000 
 22.694100 4790.000000 
 22.742025 4795.000000 
 22.790000 4800.000000 
 22.838025 4805.000000 
 22.886100 4810.000000 
 22.934225 4815.000000 
 22.982400 4820.000000 
 23.030625 4825.000000 
 23.078900 4830.000000 
 23.127225 4835.000000 
 23.175600 4840.000000 
 23.224025 4845.000000 
 23.272500 4850.000000 
 23.321025 4855.000000 
 23.369600 4860.000000 
 23.418225 4865.000000 
 23.466900 4870.000000 
 23.515625 4875.000000 
 23.564400 4880.000000 
 23.613225 4885.000000 
 23.662100 4890.000000 
 23.711025 4895.000000 
 23.760000 4900.000000 
 23.809025 4905.000000 
 23.858100 4910.000000 
 23.907225 4915.000000 
 23.956400 4920.000000 
 24.005625 4925.000000 
 24.054900 4930.000000 
 24.104225 4935.000000 
 24.153600 4940.000000 
 24.203025 4945.000000 
 24.252500 4950.000000 
 24.302025 4955.000000 
 24.351600 4960.000000 
 24.401225 4965.000000 
 24.450900 4970.000000 
 24.500625 4975.000000 
 24.550400 4980.000000 
 24.600225 4985.000000 
 24.650100 4990.000000 
 24.700025 4995.000000 
 24.750000 5000.000000 

2、环境文件

'Parabolic bottom profile'	! TITLE
10.0,			! FREQ (Hz)
1,			! NMEDIA
'CVW *',			! SSPOPT (Analytic or C-linear interpolation)
525 0.0 5100.0		! DEPTH of bottom (m)
 -5100 1500.0 /
  5100 1500.0 /
'A~' 0.0
5100.0 5000.0 0 1.0 0.0 0.0 /
1				! NSD, 
0.0 /			! SD(1:NSD) (m)
201				! NRD
    0.0 5000.0 /		! RD(1:NRD) (m)
501				! NR,   
  0.0  20.0 /			! R(1:NR) (km)
'R'	  			! 'R/C/I/S'
50		       		! NBEAMS
-89.0 89.0 /		        ! ALPHA1,2 (degrees)
0  5100.0  25.1,		! STEP (m), ZBOX (m), RBOX (km)

3、Matlab 代码

clc;    clear all;  close all
global units;	units = 'km';
% linear boundary interpolation
make_bdry( 'L' )
% the rays:
bellhop ParaBot
subplot(211); plotray ParaBot;  axis( [ 0 20 -5000 5000 ] )
plotati 'ParaBot'   % superimpose an altimetry plot
plotbty 'ParaBot'   % superimpose a bathymetry plot

bellhop ParaBothl
subplot(212); plotray ParaBothl;    axis( [ 0 20 0 5000 ] )
plotbty 'ParaBot'   

% linear boundary interpolation
make_bdry( 'C' )
bellhop ParaBot
figure; subplot(311); plotray ParaBot;  axis( [ 0 20 -5000 5000 ] )
plotati 'ParaBot';  plotbty 'ParaBot'

bellhop ParaBothc
subplot(312); plotray ParaBothc;    axis( [ 0 20 0 5000 ] )
plotbty 'ParaBot'

bellhop ParaBothcC
subplot(313); plotshd ParaBothcC.shd;   axis( [ 0 20 0 5000 ] )
plotbty 'ParaBot'

4、执行结果

从底部反射的声线应该与表面平行,就像手电筒中的反射面所产生的均匀光束一样。当我们超出 20 公里的距离时,声线的位置对底部各块小面元的倾斜非常敏感。我们每增加 25m 海深就进行一次测深采样,并将“测深插值”选项设置为“分段线性插值”,就得到下图所示的声线轨迹。下图中不规则的声线结构明显揭露了该方法的缺陷
在这里插入图片描述

分段线性拟合边界反射的声线

为了获得更加平滑的声线轨迹,我们只需将测深文件内的第一个字母更改为“C”,即设置成“曲线拟合”选项,就可实现。这种改进的边界插值提供了一组完美的平行射线(在眼睛所能辨识的范围内)。
在这里插入图片描述

曲线拟合边界插值的抛物线测深剖面;从上到下依次是声线轨迹、声线轨迹、相干传播损失。

七、资源自取

CSDN 链接:https://download.csdn.net/download/qq_41839588/87855407?spm=1001.2014.3001.5503


总结

这就是文章的全部了,案例中详细讲解了水平海底、变化的海底、Dickins 海山以及抛物线海底场景下的声传播特性。


我的qq:2442391036,欢迎交流!


本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/611686.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

TDEngine3.0 环境安装、配置及使用经验总结

TDEngine3.0 环境安装、配置及使用经验总结 一、TDengine 介绍二、TDengine的下载三、TDengine Server安装及配置3.1 安装3.2 taos的参数配置3.3 启动3.4 taosAdapter 四、TDengine Client 安装4.1 linux客户端安装4.2 windows客户端安装 五、TDEngine3.x的使用总结 一、TDengi…

minhook探究

参考:https://github.com/TsudaKageyu/minhook minhook是windows平台上支持x86/x64的hook库,git上的自我介绍说是“mininalistic",其简约并不简单。在接口的设计,hook的兼容性等方面,还是值得我们初学者解决的。熟悉inline …

Hadoop之HDFS概述

Hadoop概述之HDFS HDFS架构概述优缺点HDFS架构HDFS文件块大小HDFS的shell命令HDFS读写流程写数据流程 HDFS读数据流程NameNode 和 SecondaryNameNode工作机制DataNode工作机制DataNode数据完整性如何保证 端口名称Hadoop2.xHadoop3.xNameNode内部通信端口8020/9000NameNode HTT…

两天搞定计算机专业毕业设计,附源码

两天搞定计算机专业毕业设计,附源码 适用者毕设专业 使用要求具备基本Unity 基本操作小白即可,无需编码 博主诉求快乐毕业 点赞 关注 收藏 资源说明Free资源太多了,看截图目录就知道了 适用者 毕设专业 鄙人也是计算机狗一只,会…

软考A计划-电子商务设计师-复习要点

点击跳转专栏>Unity3D特效百例点击跳转专栏>案例项目实战源码点击跳转专栏>游戏脚本-辅助自动化点击跳转专栏>Android控件全解手册点击跳转专栏>Scratch编程案例 👉关于作者 专注于Android/Unity和各种游戏开发技巧,以及各种资源分享&am…

Ampere 又放大招,推出自研192 核AmpereOne 系列处理器,已投产

作者 | 伍杏玲 近日,Ampere Computing 发布2023年度战略和产品路线图,并推出全新的AmpereOne系列处理器,拥有多达 192 个单线程 Ampere 核,内核数量为业界最高。这是第一款基于 Ampere 新自研核的产品,由 Ampere 自有…

java--正则表达式

一、作用 作用一:校验字符串是否满足规则 作用二:在一段文本中查找满足要求的内容 二、符号含义 1、字符类(只匹配一个字符) 符号含义[abc]只能是a,b或c中一个[^abc]除了a,b,c之外的任何字符[a-zA-Z]a到z A到Z[a-d[m-p]]a到d&…

弱网测试,Network Link Conditioner你知多少

网络环境的好坏,有时会让你的产品带给用户完全不同的体验,作为开发者,在开发项目过程中,我们需要进行对于网络环境的调试。Mac环境下模拟慢速网络可以使用苹果官方提供的工具: Network Link Conditioner 具体操作步骤…

数字化时代,公司如何成为数据驱动组织

当前,数据要素和数字经济提出了数据在生产过程中发挥的重要作用。其中最热点的话题包括数据资产、数据价值、数据驱动和数字化转型。如果数据是一种资产,那么它应该为公司及其利益相关者创造价值。那么如何获从数据中获得得不同类型的价值以维持公司的竞…

Python获取当当平台商品数据信息可视化效果展示

前言 嗨喽&#xff0c;大家好呀~这里是爱看美女的茜茜呐 环境使用: 版 本&#xff1a; python 3.8 编辑器&#xff1a;pycharm 2021.2 jupyter notebook 模块使用: 采集 requests >>> pip install requests 数据请求 csv <表格文件> 内置模块 保存数据 …

图书购物商城 图书后台管理系统

图书购物商城 图书后台管理系统 这个图书购物商城系统是一个基于JSP、Servlet和MySQL技术开发的综合性系统&#xff0c;它包括前台和后台功能。该系统旨在为用户提供一个便捷的购书平台&#xff0c;同时为管理员提供管理图书和订单的功能。 前台部分是用户使用的界面&#xf…

MQTT与EMQ

文章目录 1 MQTT协议与EMQ中间件1.1 物联网消息协议MQTT1.1.1 什么是MQTT1.1.2 MQTT相关概念1.1.3 消息服务质量QoS——信息的可靠投递1.1.3.1 QoS0——消息服务质量为0&#xff0c;消息发送至多一次1.1.3.2 QoS1——消息发送至少一次1.1.3.3 QoS2——消息发送仅一次1.1.3.4 不…

app渗透-常见问题及绕过

app渗透-常见问题及绕过 6.app常见问题和绕过前言6.1反代理操作前言6.1.1判断6.1.2实例演示-探探6.1.3绕过1-r0capture6.1.4绕过2-proxifier6.1.5绕过3-小黄鸟 6.2证书校验前言6.2.1判断6.2.2浏览器校验和解决6.2.3桡过证书单项校验-xp框架6.2.3绕过证书双向校验 6.app常见问题…

TongWeb8适配JakartaEE应用

历史&#xff1a; 2017年 Oracle将Java EE&#xff08;Java SE还自己保留&#xff09;交给开源组织&#xff0c;Eclipse基金会接手。但Oracle不允许开源组织使用Java名号&#xff0c;所以Jakarta EE名称于2018.02.26应运而生。 版本 发布日期 焦点说明 Java EE 8 2017.08 …

30天网络安全从入门到精通?

前言 毫无疑问&#xff0c;网络安全是当下最具潜力的编程方向之一。对于许多未曾涉足计算机编程的领域「小白」来说&#xff0c;深入地掌握网络安全看似是一件十分困难的事。至于一个月能不能学会网络安全&#xff0c;这个要看个人&#xff0c;对于时间管理不是很高的&#xf…

花3个月面过阿里测开岗,拿个25K不过分吧?

计算机专业&#xff0c;代码能力一般&#xff0c;之前有过两段实习以及一个学校项目经历。第一份实习是大二暑期在深圳的一家互联网公司做前端开发&#xff0c;第二份实习由于大三暑假回国的时间比较短&#xff08;小于两个月&#xff09;&#xff0c;于是找的实习是在一家初创…

南山村又一个旧改项目即将开工建设,桂庙新村城市更新单元。

5月22日&#xff0c;深圳市南山区城市更新和土地整备局发布关于粤海街道桂庙新村城市更新单元项目实施主体确认的公示。 根据公告&#xff0c;该项目实施方式为与权利主体签订搬迁补偿协议&#xff0c;补偿方式为产权置换和货币补偿相结合&#xff0c;已完成100%签约&#xff0…

【Linux】Nginx编译安装及系统服务添加

提示&#xff1a;文章写完后&#xff0c;目录可以自动生成&#xff0c;如何生成可参考右边的帮助文档 Nginx 一、关于 NginxNginx和Apache的差异Apache的优势Nginx的优势 二. 编译安装Nginx服务1.关闭防火墙2.安装依赖包3.创建运行用户、组4.编译安装Nginx5.检查、启动、重启、…

mac电脑git clone项目时报错证书过期和权限被拒绝

mac电脑使用git clone命令克隆项目时&#xff0c;一开始一直提示证书过期 SSL certificate problem: certificate has expired 执行以下代码关掉验证后&#xff0c;解决了这个问题 找到git目录 Git\git-cmd输入命令跳转到bin目录&#xff0c;cd bin输入命令运行git.exe执行关…

2023网安164道面试题(附答案)

最近有不少小伙伴跑来咨询&#xff1a; 想找网络安全工作&#xff0c;应该要怎么进行技术面试准备&#xff1f;工作不到 2 年&#xff0c;想跳槽看下机会&#xff0c;有没有相关的面试题呢&#xff1f; 为了更好地帮助大家高薪就业&#xff0c;今天就给大家分享两份网络安全工…