做水声装备这些年，我接触最多的需求就是“怎么把水下看明白”。侧扫声呐只能给你一幅声影图，看不出精确深度；单波束测深仪又一针一针地打，效率太低。直到多波束成像声呐出现，这个问题才算真正解决——它用一排换能器同时发射和接收声波，通过波束形成把一次扫掠变成一条覆盖带，同时拿到几十到几百个方向的回波信息，既测深又成像。这篇文章我会从原理讲起，然后把仿真模型的搭建思路和一套可直接运行的Matlab代码完整拆开，适合水声专业的学生、刚进入声呐算法岗的工程师，以及想快速验证多波束信号处理方案的朋友。
1. 多波束成像声呐到底在解决什么问题
1.1 核心价值和典型应用场景
多波束成像声呐的本质，是用相控阵的思路把水下的“单视角”变成“宽覆盖”。传统单波束声呐只能垂直向下打一个点，测完一整片区域就得船来回跑无数遍。多波束声呐在垂直航迹方向排列一个接收基阵，发射声波覆盖一个扇面，接收时对不同方向的回波分别做相位补偿，相当于同时形成了大量指向不同角度的窄波束。这样一来，船只只要沿着测线走一趟，就能得到一条横跨航迹的海底剖面，连续走完就是完整的三维地形。
它的应用场景非常集中：海洋测绘、航道疏浚、水下目标探测、ROV/AUV避障与环境感知。海洋测绘要的是高精度的水深点，多波束能给出带角度和距离的密集点云；水下目标探测则更看重成像分辨率和对比度，多波束成像声呐能够把水雷、沉船、管线这类目标从混响背景里区分出来。所以无论你是做测绘还是做水下安防，这套原理都是绕不开的核心。
1.2 仿真分析要回答哪几个关键问题
我见过不少刚接触声呐仿真的朋友，一上来就急着写代码，结果方向不对，仿真做完了也不知道结果是否可信。其实多波束成像声呐的仿真分析，最终要回答的无非就三个问题。
第一，波束到底能不能把方向接近的目标分开，也就是角度分辨力是否满足指标。第二，回波经过匹配滤波之后，距离向能不能准确还原出目标的相对位置，也就是距离分辨力和测距精度。第三，整个系统在设定信噪比、混响背景下，能不能稳定检测出目标，这涉及到波束形成的增益和信号处理余量。
把这几个问题拆开之后，你会发现仿真建模其实不需要一上来就搞复杂的声场传播模型。对大多数工程验证来说，点目标加高斯白噪声的模型已经足够说明问题，等算法跑通了再逐步加入多径、混响、阵列幅相误差这些现实因素，这才是效率最高的路径。
2. 多波束成像声呐的基本原理拆解
2.1 波束形成的本质：相位差就是方向信息
理解多波束声呐，最先要过的坎就是波束形成。我换个方式讲：假设水下一个目标离阵列很远，远到声波到达不同阵元时幅度近似相等，但路程差不能忽略。这个路程差会导致同一个信号到达不同阵元时出现相位差，而这个相位差恰恰和目标方向有关。
比如均匀直线阵，阵元间距为d，目标方向与法线夹角为θ，相邻阵元间的声程差为d·sinθ，对应的相位差为2πd·sinθ/λ。如果我们给每个阵元的输出人为补偿一个相反的相位延迟，再把所有阵元相加，那么从θ方向来的信号会因同相叠加而增强，从其他方向来的信号则因相位不一致而被削弱，这就是延时求和波束形成的基本原理。
在工程实现上，窄带信号用复加权来实现相移，宽带信号则需要做真正的时延补偿，或者用频域相位补偿。多波束成像声呐通常发射宽带信号来保证距离分辨力，所以仿真中如果只做窄带假设，就需要明确说明适用范围。对于原理验证，窄带复加权已经能很好地展示波束形成的核心逻辑。
2.2 距离分辨力和角度分辨力从哪来
分辨率是声呐最核心的指标之一。角度分辨力由阵列孔径和波束宽度决定，均匀加权时半功率波束宽度近似为0.89λ/L，其中L是阵列孔径长度，λ是波长。想提高角度分辨力，要么降低频率（但会牺牲对目标的细节刻画能力），要么增大阵列尺寸。工程上阵元数、阵元间距和中心频率三者互相牵制，仿真时必须一起考虑。
距离分辨力则主要取决于发射信号的带宽。使用单频脉冲时，距离分辨力为c·τ/2，τ是脉冲宽度，脉冲越短分辨力越高，但能量又不够，探测距离受限。而线性调频（LFM）信号能把能量摊到很宽的频带上，接收端通过匹配滤波把信号压缩成窄脉冲，压缩后的等效脉宽约为1/B，距离分辨力变为c/(2B)。这就是LFM信号几乎统治声呐发射波形的原因——它同时解决了能量和分辨力之间的矛盾。
匹配滤波带来的处理增益约为时间带宽积B·T的平方根，或者用dB表示为10·log10(B·T)。举个例子，带宽8kHz、脉宽2ms的LFM信号，B·T=16，匹配滤波增益大约12dB。这个增益对远距离弱目标检测非常关键，仿真的信噪比设置也应该把这个增益考虑进去。
2.3 声呐方程和系统指标的换算关系
做仿真之前，最好先把声呐方程的账算清楚。主动声呐方程一般写成：
SL - 2TL + TS - (NL - DI) = DT
其中SL是发射声源级，TL是传播损失，TS是目标强度，NL是环境噪声级，DI是接收指向性指数，DT是检测阈。这个式子的意思是：声波从声源传到目标再回来，能量一路衰减，但系统通过波束形成抑制了各向同性噪声，最后能剩下的信噪比刚好要超过检测阈值。
仿真中如果直接把目标回波加到高斯白噪声里，实际上默认了TL、TS、NL这些量已经被简化成某一信噪比的等效值。这条路径完全可行，但有个前提是你清楚自己跳过了什么。等需要评估真实探测距离时，再回到声呐方程逐项代入。
2.4 完整的工作流程一次串起来
多波束成像声呐的一帧处理流程大致是：发射基阵发出一个覆盖扇面的LFM脉冲，声波遇到海底和水中目标后产生散射回波，接收基阵的每个阵元并行采集回波信号，经过放大、滤波和模数转换后进入数字信号处理模块。
数字信号处理模块里，先对每一阵元做匹配滤波或频域脉压，再对不同波束指向做相位补偿和累加，得到一组按角度排列的一维距离像。最后根据波束指向角和回波时延，换算出每个点的斜距和深度，再结合航向姿态数据，生成地理坐标下的点云。后面的成像和显示，都是在这个数据基础上做的。
3. 仿真分析怎么建模才靠谱
3.1 坐标系和阵列方案怎么定
仿真第一步是定义坐标系。我做多波束仿真时习惯用右手系：声呐随平台沿X轴方向运动，Y轴是跨航迹方向，Z轴向下为正，换能器阵面朝Z方向。目标位于阵列前方某个方位角和俯仰角上，声波沿直线路径传播到阵元。这个坐标系既方便生成回波，也方便最后把结果映射成点云。
阵列方案上，最常用的是均匀线阵或均匀平面阵。均匀线阵只能分辨一个维度的角度，适合剖面测量；均匀平面阵可以同时分辨方位和俯仰，适合三维成像。仿真入门阶段，建议先用均匀线阵把波束形成、匹配滤波的逻辑跑通，再扩展到平面阵。
阵元间距一般取半波长，既能保证不出现栅瓣，又能在阵元数固定的情况下获得最大的阵列孔径。中心频率的选择要结合目标尺寸和探测距离：频率高分辨率好但衰减大，频率低传播远但分辨率差。浅水高分辨率成像常用几十千赫到几百千赫，这个范围在仿真中比较容易处理。
3.2 点目标模型的数学表达
回波信号的建模是整个仿真中最关键的环节。假设阵列有N个阵元，目标在方位角θ处，距离为R，那么第n个阵元接收到的目标回波可以写成发射信号的延迟复制，延迟由目标到该阵元的距离除以声速决定。
窄带假设下，不同阵元的回波差异可以简化为相位差，即频域乘一个相位因子exp(-j2πd·n·sinθ/λ)。宽带情况下，延时τ_n = (R + d·n·sinθ)/c必须直接作用在时间轴上，回波为s_tx(t - τ_n)。仿真中如果不做宽带处理而只做窄带近似，会对大孔径阵列造成明显的波束畸变，这一点要特别注意。
目标幅度可以按声呐方程折算，也可以直接设一个信噪比，让回波加噪声后仍能被检测到。多个目标时，总回波是各目标回波的线性叠加，前提是忽略目标之间的二次散射，这对大多数水声仿真场景是合理的。
3.3 仿真流程设计的几个层次
我建议把仿真拆成三个递进的层次。第一层是静态单帧：固定目标和阵列位置，仿真波束形成、匹配滤波、检测和测向测距，主要验证算法逻辑和分辨率指标。第二层是航迹扫描：让阵列沿X轴运动，每个航迹位置生成一帧回波，把所有帧的测点拼成点云，验证覆盖效果。第三层才是复杂环境：加入海底混响、多径、阵列幅相误差、平台姿态扰动等非理想因素，评估系统的鲁棒性。
绝大多数项目的坑都出在第一层没做扎实就急着上第三层。一个点目标下波束形成方向图都对不齐，后面加再多噪声都是白搭。所以这篇文章里重点演示第一层的完整实现，让大家先把基础骨架搭稳。
4. 核心仿真代码实现与结果解读
4.1 一套可直接运行的Matlab仿真
下面这套代码实现了均匀线阵下的多波束成像仿真，包含发射LFM信号、生成多阵元回波、匹配滤波、延时求和波束形成和方向图绘制。代码尽量精简，但每个关键环节都有注释。
% 多波束成像声呐仿真：LFM信号 + 均匀线阵 + 常规波束形成 clear; clc; close all; %% 参数设置 c = 1500; % 声速 m/s f0 = 40e3; % 中心频率 Hz B = 8e3; % 信号带宽 Hz T = 2e-3; % 脉冲宽度 s fs = 200e3; % 采样率 Hz lambda = c / f0; % 波长 m N = 32; % 接收阵元数 d = lambda / 2; % 阵元间距 m theta_target = [-20 8 33]; % 三个目标的方位角 range_target = [50 60 80]; % 三个目标的斜距 m snr_per_element = 5; % 单阵元回波信噪比 dB %% 发射信号：LFM t = 0 : 1/fs : T - 1/fs; s_tx = exp(1j * pi * (B/T) * t.^2); % 复基带LFM K = B / T; %% 生成接收回波 % 计算目标回波时延(参考阵元到目标的往返时延) tau_ref = 2 * range_target / c; % 采样点位置 N_samples = round(T * fs) + 500; r = zeros(N, N_samples); for m = 1:N % 第m个阵元相对参考阵元的声程差 delta_tau = (m - 1) * d * sind(theta_target) / c; for k = 1:length(theta_target) tau = tau_ref(k) + delta_tau(k); % 回波起始采样点 n0 = round(tau * fs) + 1; if n0 < 1 || n0 + length(s_tx) - 1 > N_samples continue; end % 衰减系数：距离平方反比简化模型 amp = 1 / (range_target(k)^2); r(m, n0 : n0 + length(s_tx) - 1) = ... r(m, n0 : n0 + length(s_tx) - 1) + amp * s_tx; end end % 加噪 noise_power = 10^(-snr_per_element / 10) * mean(abs(r(:)).^2); r = r + sqrt(noise_power / 2) * (randn(N, N_samples) + 1j * randn(N, N_samples)); %% 匹配滤波(参考信号共轭反转) mf = conj(fliplr(s_tx)); y = zeros(N, N_samples); for m = 1:N temp = conv(r(m, :), mf); y(m, :) = temp(1:N_samples); end %% 波束形成：对不同方位角扫描 theta_scan = -60 : 1 : 60; beam_out = zeros(length(theta_scan), N_samples); for idx = 1:length(theta_scan) % 计算扫描角度的导向矢量 steer_tau = (0:N-1)' * d * sind(theta_scan(idx)) / c; % 频域相位补偿（窄带近似），乘以共轭相位 w = exp(-1j * 2 * pi * f0 * steer_tau); beam_out(idx, :) = w' * y / N; end %% 提取目标峰值并进行角度-距离映射 distance_axis = c * (0:N_samples-1) / (2 * fs); figure; imagesc(theta_scan, distance_axis, 20*log10(abs(beam_out).' + eps)); set(gca, 'YDir', 'normal'); xlabel('角度 (deg)'); ylabel('斜距 (m)'); title('距离-角度二维成像图'); colorbar; caxis([-40 0]);运行这段代码，你会看到一个典型的距离-角度二维图，三个目标应该出现在对应角度和距离的交叉位置。需要注意，窄带波束形成在宽带宽和大阵元数条件下会有一定的波束指向偏移，这是近似模型导致的，真实系统中通常会用宽带时延补偿替代窄带相位补偿。
4.2 波束方向图和分辨率怎么看
为了让方向图更直观，我通常再画一张静态方向图，看主瓣宽度和旁瓣电平。
%% 画波束方向图 theta_plot = -90 : 0.1 : 90; pat = zeros(size(theta_plot)); w0 = exp(-1j * 2 * pi * f0 * (0:N-1)' * d * sind(0) / c); for idx = 1:length(theta_plot) a_theta = exp(1j * 2 * pi * f0 * (0:N-1)' * d * sind(theta_plot(idx)) / c); pat(idx) = abs(w0' * a_theta) / N; end figure; plot(theta_plot, 20*log10(pat + eps), 'LineWidth', 1.5); grid on; xlabel('方位角 (deg)'); ylabel('归一化幅度 (dB)'); title('32元均匀线阵波束方向图'); ylim([-60 0]);从方向图上能直观看到两个信息：主瓣越窄，角度分辨力越高；旁瓣越低，对旁瓣方向干扰的抑制能力越强。均匀加权时第一旁瓣约-13dB，如果目标强弱悬殊，弱目标很容易被强目标的旁瓣掩盖。这时候引入海明窗或泰勒窗可以压低旁瓣，但主瓣会展宽，对应的角度分辨力会变差，这是典型的取舍问题。
4.3 测角和测距精度怎么评估
仿真得出了峰值位置，接下来要量化误差。测距误差主要由采样率决定，采样率越高，时延量化误差越小。测角误差则受波束宽度和信噪比共同影响，信噪比越高，峰值检测越准。
我在实际项目中通常做蒙特卡洛仿真，就是在同一信噪比下重复几百次，统计测角误差和测距误差的均方根值。比如信噪比从0dB提高到20dB，测角误差往往会从零点几度下降到零点零几度。这类曲线是评估声呐系统性能最直接的证据，也是写方案报告时最有说服力的材料。
5. 常见问题与工程避坑指南
5.1 栅瓣问题：阵元间距是红线
栅瓣是多波束声呐最容易踩的坑。当阵元间距超过半波长时，除了主瓣之外会出现与主瓣幅度几乎相同的栅瓣，目标在栅瓣方向会被误判成主瓣方向的目标。仿真中如果发现成像图里出现对称的假目标，第一件事就检查阵元间距是否满足d不超过λ/2。
工程上有时因为空间限制，阵元间距做不到半波长，这时可以降低中心频率来增大波长，但要牺牲分辨率，或者在算法上使用稀疏阵列构型，让栅瓣变成可控的较高旁瓣。切记栅瓣无法通过加窗消除，只能从阵列设计层面规避。
5.2 幅相误差让波束变差
真实阵列各通道的幅度和相位不可能完全一致，温度漂移、器件公差、阵元老化都会引入误差。波束形成对这种误差比较敏感，误差大会导致主瓣抬升、旁瓣抬高，甚至波束指向偏移。
这个问题我建议在仿真阶段就加进去：给每个阵元回波乘一个随机复增益，幅度误差设为几个百分点到百分之十几，相位误差设为几度到十几度，再跑一遍方向图看恶化程度。如果系统设计需要旁瓣低于-30dB，幅相误差的控制要求会很严格，这对硬件选型和校准方案都有指导意义。
5.3 多径和混响的干扰
多径指声波经海面反射或海底反射后再到达接收阵，会形成时间上拖尾、方向上偏移的虚假回波。混响则是大量随机散射体产生的背景干扰，在近程尤其严重。仿真中混响常用高斯或瑞利分布随机过程模拟，多径则简化成几个离散的镜反射路径。
真实场景里，多波束声呐最容易受多径影响的是浅水区，因为海面和海底距离近，反射路径多。处理手段一般是利用波束形成的空间选择性抑制多径方向，或使用脉冲压缩后的时间窗把多径回波滤除。这些方法在仿真中验证起来非常方便，建议多试几种策略再上实机。
5.4 参数调整的实战心得
仿真中参数调整遵循“先分辨率后增益”的原则。先把带宽设大一些，确认距离维能把目标分开；再把阵元数增大，确认角度维能把目标分开；最后才根据虚警率调整检测门限。如果一上来就把参数调到最优，出了问题很难定位是阵列问题还是信号处理问题。
还有一个细节，低信噪比下峰值检测容易跳变，这时可以在检测前做幅度平滑或使用质心插值算法。质心法在波束域多点加权求峰值，比直接取最大点稳定得多，能够把测角精度提升半个到一个波束宽度级别。
6. 最后说点我自己的体会
从第一次跑通多波束仿真到现在，我最大的感受是：声呐仿真相当考验人对物理模型的判断力。代码只是工具，真正有价值的是你清楚每一个信号的时延从哪来、每一个波束的指向由什么决定、每一个误差项会以什么方式破坏结果。把这个逻辑捋顺，后面不管是做海底测绘点云处理，还是做水下目标自动识别，都会顺手很多。
如果你刚开始接触这套东西，别急着追求复杂的声场模型，先把这篇文章里的点目标仿真吃透，把波束形成和匹配滤波在自己手里跑一遍，再看更高阶的文献就会觉得轻松许多。仿真数据自己会说话，前提是你知道怎么问它问题。
