从传播损耗到阵元波形：Kraken复数声场提取与Matlab实现

简介这份资源面向水声通信、海洋声学方向的研究人员与工程技术人员提供基于 Matlab 的 Kraken 计算工具用于求解水声传播损耗并模拟水平阵、垂直阵各阵元的接收波形可服务于声纳系统设计、阵列性能评估与信号检测等场景。压缩包共 19 个文件约 713KB包含 3 个 m 脚本分别对应传播损耗与水平阵、垂直阵波形计算、1 个 mat 声速剖面数据、1 个 docx 公式说明以及 12 张 jpg 与 1 张 png 结果图直观呈现传播损耗曲线、声场分布与阵元波形。资源已有 1109 人学习下载配套 README 与任务说明图给出运行流程与参数含义读者可据此复现计算、修改海洋环境参数并对照结果图理解阵列因子与方向图的变化规律为后续自定义扩展打下基础。1. 声场计算里最容易被忽略的一步从传播损耗到阵元波形很多人用 Kraken 算完声场拿到一条 TL 曲线就以为任务结束了。但真正做过阵列信号处理的人会发现传播损耗只是标量幅度它回答不了「每个阵元到底收到什么波形」这个问题。Kraken 输出的场文件里其实藏着复数声压——幅度和相位都在里面只是大多数人只取了abs就丢掉了相位。而阵列波束形成、匹配场处理、时延估计这些下游任务恰恰依赖各阵元之间的相位差和波形细节。换句话说传播损耗是「结果」阵元波形才是「原料」。搞懂怎么从 Kraken 的场输出里把每个阵元的复数波形还原出来以及背后的简正波展开为什么能同时给出幅度和相位是这篇文章要拆开讲清楚的事。Matlab 在这里扮演的角色不是画图工具而是把 Kraken 的二进制场文件翻译成可用的复数矩阵。2. 简正波展开为什么 Kraken 能同时给出幅度和相位2.1 从波动方程到模态求和的形式Kraken 求解的是水平分层介质中的 Helmholtz 方程。在柱坐标下假设介质只随深度变化声压可以写成距离和深度的分离形式。分离变量之后深度方向满足一个本征值问题解出来是一组模态函数 Ψ_m(z) 和对应的水平波数 k_m。最终声压的表达式是一个模态求和p(r, z) (i / (8π r))^(1/2) * exp(-iπ/4) * Σ_m Ψ_m(z_s) * Ψ_m(z) * exp(i k_m r) / sqrt(k_m)这个式子里每一项都带相位exp(i k_m r) 是沿距离的传播相位Ψ_m(z) 是模态函数在接收深度处的值本身可能是复数前面的常数因子也有固定相位。所以 Kraken 输出的场天然是复数传播损耗只是它的模。常见的实现思路是Kraken 先把本征值和模态函数算出来存到.mod文件然后由场计算程序按上面的求和式在指定距离和深度网格上累加结果写进.shd文件。.shd的格式是二进制包含实部和虚部交替存储的声压值Matlab 读进来之后就是一个复数矩阵行对应深度列对应距离。理解这一点很关键传播损耗 TL -20*log10(|p|/|p_ref|)只用了模而阵元波形需要的是 p 本身包括相位。如果你只存了 TL相位信息就永久丢了后面想补都补不回来。这是很多人第一次做阵列处理时踩的坑——回头重新跑一遍 Kraken白白浪费计算时间。2.2 模态函数的物理含义与截断选择模态函数 Ψ_m(z) 描述的是第 m 号简正波在深度方向的驻波结构。它有几个零点零点之间相位翻转。对于阵列接收来说不同深度的阵元落在同一个模态的不同位置上接收到的该模态分量幅度和相位都不同。模态求和时各模态之间还会干涉导致总声压随深度出现起伏——这就是为什么阵列各阵元的波形不是简单的幅度缩放关系。Kraken 在求解本征值时需要指定模态数。模态数不够高频或大距离下会漏掉贡献显著的模态波形失真模态数过多计算量上升而且高阶模态可能已经衰减到数值噪声水平反而引入误差。常见的做法是先跑一次看本征值分布确认最高阶模态的衰减足够大再定截断。我一般会检查.mod文件里各模态的衰减系数确保最后一个模态在最大传播距离上的衰减超过几十 dB这样截断才安全。提示模态数不是越多越好。高阶模态在近距离可能贡献明显但在远距离已经衰减殆尽盲目加模态只会增加计算时间而不改善精度。2.3 场文件里的复数结构实部虚部怎么排.shd文件的二进制布局是理解阵元波形提取的关键。文件头包含频率、深度点数、距离点数、深度网格和距离网格等元信息之后是声压数据区。数据区按距离列优先存储每个距离上沿深度排列每个点存两个 float32——实部和虚部。Matlab 读取时需要先解析文件头拿到维度再按[2*Nz, Nr]的形状读数据然后实部虚部交替拆开合成复数矩阵。如果读错了字节序或者维度顺序得到的波形会完全乱掉但幅度看起来可能还「差不多」这种错误特别隐蔽。血泪经验是读完之后先检查第一个距离上的声压是否随深度平滑变化如果出现剧烈跳变多半是维度或字节序搞错了。另外要注意 Kraken 的距离网格通常从第一个距离点开始不一定包含 r0提取阵元波形时要按实际距离网格插值到阵元位置不能想当然。3. 从场文件到阵元波形Matlab 里的复数矩阵操作3.1 读取 shd 文件并还原复数声压下面这段代码演示如何解析.shd文件头并读出复数声压矩阵。不同版本的 Kraken 文件头字段顺序可能略有差异实际使用时需要对照文档核对字段偏移量。function [pressure, z_grid, r_grid] read_shd(filename) fid fopen(filename, r, ieee-le); % Kraken 默认小端 % 文件头前几个字段依次为频率、Nz、Nr 等 freq fread(fid, 1, float32); Nz fread(fid, 1, int32); Nr fread(fid, 1, int32); % 跳过剩余头字段具体偏移依版本而定 fseek(fid, 16, cof); z_grid fread(fid, Nz, float32); r_grid fread(fid, Nr, float32); % 数据区每个点实部虚部交替按距离列优先 raw fread(fid, 2*Nz*Nr, float32); fclose(fid); raw reshape(raw, 2, Nz, Nr); % 拆实部虚部 pressure squeeze(raw(1,:,:)) 1i * squeeze(raw(2,:,:)); end这段代码的核心在最后三行reshape把一维数据按[2, Nz, Nr]重排第一维是实部虚部后两维是深度和距离squeeze去掉长度为 1 的维度后合成复数。参数说明Nz是深度网格点数Nr是距离网格点数freq是计算频率。注意fseek跳过的字节数需要根据实际文件头结构调整这里只是示意。读出来的pressure是Nz × Nr的复数矩阵每一列是一个距离上的深度剖面每一行是一个深度上的距离剖面。3.2 按阵元深度和距离插值提取波形阵元通常不在 Kraken 的网格点上需要二维插值。深度方向插值到阵元深度距离方向插值到阵元距离。由于声压是复数实部和虚部要分别插值再合成不能对模和相位分别插值——相位在模态干涉处会跳变直接插相位会出错。% 假设阵元深度 z_elem (Ne×1)阵元距离 r_elem (标量或 Ne×1) [Z, R] meshgrid(z_grid, r_grid); % 注意维度对齐 P_real interp2(Z, R, real(pressure)., z_elem, r_elem, spline); P_imag interp2(Z, R, imag(pressure)., z_elem, r_elem, spline); waveform P_real 1i * P_imag; % 各阵元复数波形这里有个容易翻车的地方pressure是Nz × Nr而meshgrid生成的Z、R是Nr × Nz所以要对pressure转置再插值否则深度和距离会对调。插值方法选spline还是linear取决于网格密度网格足够密时linear就够用网格稀疏时spline更平滑但可能过冲。提取出来的waveform是各阵元的复数声压幅度就是该阵元的接收强度相位就是相对于参考点的传播相位直接可以送进波束形成器。3.3 传播损耗与阵元波形的关系验证提取完波形后应该做一次一致性检查取某个阵元的幅度算 TL和 Kraken 直接输出的 TL 曲线对比。如果对不上说明读取或插值环节有问题。TL 的定义是-20*log10(abs(p)/abs(p_ref))参考声压通常取 1 米处的声压。这个检查能帮你快速定位是文件读错了还是插值搞错了比盲目调下游算法高效得多。4. 阵列处理中的常见问题与排查4.1 波形幅度对但相位全错现象各阵元幅度分布合理但波束形成后主瓣指向偏离预期或者匹配场处理的相关峰位置不对。原因读取.shd时实部虚部顺序搞反或者字节序设错导致相位符号翻转。解决用单个距离上的深度剖面做检查正常声压随深度应该是连续变化的如果实部虚部对调剖面会出现不自然的剧烈振荡。确认文件头里的字节序标记Kraken 在多数平台上输出小端但跨平台使用时要注意。4.2 插值后阵元间相位差出现跳变现象相邻阵元的相位差本应平滑变化实际却出现 2π 跳变或突变。原因对相位直接插值而相位在模态干涉零点附近本身不连续。解决始终对实部和虚部分别插值合成后再取相位。如果已经存了相位数据先转回复数再插值不要试图对相位做 unwrap 后插值unwrap 在低信噪比处同样会出错。4.3 模态数不足导致远距离波形失真现象近距离阵元波形和参考解吻合远距离幅度偏小、相位偏差增大。原因高阶模态被截断而远距离上这些模态虽然衰减大但在某些深度上仍有不可忽略的贡献尤其是浅海负梯度环境下。解决增加模态数重新计算检查.mod文件里最高阶模态的衰减是否足够。经验做法是把模态数加到本征值谱出现明显间隙为止。4.4 距离网格不含阵元位置导致外推误差现象阵元距离超出 Kraken 计算的距离范围插值时被外推波形幅度异常。原因interp2对超出网格的点做外推结果不可靠。解决确保 Kraken 的距离网格覆盖所有阵元位置或者在提取前把距离网格加密。不要依赖外推外推的相位误差会直接毁掉阵列处理结果。4.5 频率与网格不匹配现象换频率后波形整体偏移或幅度系统性偏差。原因不同频率下 Kraken 的深度网格可能不同阵元深度插值时用了旧网格。解决每次计算后重新读取网格不要缓存上一次的网格变量。这个坑在多频点批量计算时特别常见。5. 设计取舍为什么值得从 TL 走到复数波形Kraken 把声场算成复数是简正波理论的自然结果不是额外开销。但很多使用流程里人们习惯性地只保留 TL把相位丢掉原因是 TL 直观、好画图、好对比。这个习惯在单点传播损失评估里没问题一旦进入阵列处理就翻车——波束形成靠的就是阵元间的相位差没有相位就没有指向性。从 TL 到复数波形本质上是从「能量视角」切换到「信号视角」前者回答「衰减了多少」后者回答「收到了什么」。这个切换带来的代价是数据量和处理复杂度上升复数矩阵是实数 TL 的两倍存储插值要分实部虚部验证要多做一步一致性检查。但收益是下游所有阵列算法都能直接用不用反复回跑 Kraken。我一般会在第一次计算时就把.shd完整读成复数矩阵存下来后续所有阵元提取都从这个矩阵走避免重复解析二进制文件。另一个取舍是插值方法linear快但网格稀疏时精度不够spline精度高但在模态零点附近可能过冲实际用的时候我会先用两种方法对比差异超过阈值就加密网格而不是换插值方法。验证方法上除了前面说的 TL 一致性检查还可以做一个简单的物理检查把各阵元波形按球面波假设做时延对齐如果对齐后相关性明显提升说明相位提取是对的如果对齐后反而更差相位大概率有问题。这个检查不需要知道真实声源位置只依赖阵列几何适合快速排查。从那以后我每次跑完 Kraken都会强制走一遍「读复数场 → 插值到阵元 → TL 一致性检查」的流程确认相位没问题再往下做。希望帮到你。本文还有配套的精品资源点击获取