简介该资源面向水声通信、海洋声学方向的研究人员与工程技术人员提供基于Matlab的Kraken水声传播损耗计算与阵列接收波形分析工具可用于声纳系统设计评估、水平阵与垂直阵接收性能分析等场景适合具备一定水声基础的中高级学习者。压缩包共19个文件约713KB以12张jpg结果图、3个m脚本、1个mat数据文件为主辅以docx公式说明、png任务流程图和md使用文档涵盖传播损耗曲线、声场分布、阵元波形等可视化输出与可运行代码。已有1109人学习下载。读者可据此获得从环境参数输入到传播损耗计算、阵列因子与方向图模拟的完整流程结合脚本注释理解算法逻辑并通过结果图直观对照不同条件下的声传播差异为自定义扩展与阵列配置优化提供参考。1. 从一条“对不上”的波形说起Kraken 到底在算什么很多人第一次用 Kraken 做水声传播仿真都会撞上一个反直觉的现象明明声源级设得没问题接收阵元画出来的波形却幅度忽大忽小相位还跟理论上的球面扩展对不上。有人怀疑是网格太粗有人怀疑是边界条件写错折腾半天才发现——问题出在“传播损耗”和“阵元波形”这两件事本来就不是同一个物理量却被当成一回事去校验。Kraken 作为简正波normal mode求解器它输出的不是一条简单的衰减曲线而是一整套模态叠加的场解传播损耗是场强随距离的包络阵元波形是各阶模态在接收点相干叠加后的时域结果。搞懂这两者怎么从同一个本征值问题里长出来你才能解释那些“玄学”般的幅度跳变也才能在 Matlab 里把 Kraken 的输出接对、用对。这篇就拆这条链路从简正波方程到模态系数再到传播损耗曲线和阵元时域波形一层层落到实现细节。2. 简正波求解器Kraken 为什么先解本征值再谈传播2.1 从波动方程到本征值问题模态是怎么被“逼”出来的水声传播在柱对称假设下声压场可以写成距离和深度的分离形式。把简谐点源的波动方程做分离变量深度方向会得到一个标准的 Sturm-Liouville 型本征值问题算子作用在模态函数上等于本征值乘以模态函数本征值的平方根就是水平波数。Kraken 的核心工作就是数值求解这个本征值问题得到一组离散模态和对应的水平波数。这一步决定了后面所有东西——传播损耗的衰减斜率、阵元波形的干涉条纹全部由这组本征值撑起来。常见的实现思路是用有限差分把深度方向离散成网格把微分算子变成三对角或带状矩阵再用移位反幂法迭代求本征值和本征向量。Kraken 的输入文件里那些网格点数、深度范围、边界条件参数本质上都是在描述这个离散化。网格太粗高阶模态的本征值会偏衰减算出来就不对网格太细计算量上去但精度未必线性提升。我一般会这么理解这一步是在给整个场“定调”调定歪了后面波形怎么调都是错的。2.2 模态系数声源和边界条件如何决定“谁被激发”解出本征值和模态函数只是拿到了一套“基”真正决定每个模态被激发多少的是模态系数。它由声源深度处的模态函数值、声源强度以及归一化条件共同决定。换句话说同样的海洋环境声源放浅一点或深一点激发的模态组合就不同传播损耗曲线的形状会跟着变。这也是为什么有人换个声源深度就发现损耗曲线“翻车”——不是求解器坏了是模态系数重新分配了。在 Matlab 里对接 Kraken 输出时模态系数通常需要从输出文件里读出来再和模态函数做加权。这里有个容易忽略的点归一化方式不同系数的量纲就不同直接拿去算损耗会差一个常数因子。常见的做法是确认输出文件里的归一化约定再统一到同一套单位下否则后面阵元波形的绝对幅度永远对不上。2.3 水平波数传播损耗曲线的“斜率”从哪来水平波数是本征值的平方根它的实部决定相位随距离的变化虚部决定衰减。传播损耗随距离的包络本质上就是各阶模态按水平波数虚部衰减、再相干叠加的结果。近距离时高阶模态贡献大损耗曲线可能有起伏远距离高阶模态衰减殆尽只剩低阶模态主导曲线趋于平滑。这就是为什么传播损耗曲线不是一条单调直线而是分段有“拐点”的。理解这一点你在 Matlab 里画损耗曲线时就不会盲目拟合一条指数衰减。正确的做法是按模态逐阶算贡献再叠加而不是拿总场直接套一个球面扩展公式。下面这段伪代码展示了从本征值到损耗的基本链路% 假设已从 Kraken 输出读入k水平波数复数、phi模态函数深度×模态、coef模态系数 r 100:100:10000; % 距离向量单位米 loss zeros(size(r)); for i 1:length(r) % 逐阶模态叠加exp(1i*k*r) 含相位与衰减1/sqrt(k*r) 是柱面扩展因子 field sum(coef .* phi(recv_depth_idx, :) .* ... exp(1i * k * r(i)) ./ sqrt(k * r(i))); loss(i) -20 * log10(abs(field)); % 传播损耗单位 dB end这段代码里k的虚部直接控制衰减速度coef和phi的乘积决定每个模态的权重1/sqrt(k*r)是柱面扩展的几何因子。参数recv_depth_idx是接收深度在网格上的索引取错深度波形和损耗都会偏。实际使用时还要注意k的符号约定不同版本输出可能差一个共轭导致相位反向。3. 传播损耗为什么曲线会“分段”而不是一条直线3.1 模态衰减与干涉损耗曲线的两种成因传播损耗曲线上的起伏来源有两个一是模态自身的指数衰减二是模态之间的干涉。近距离段多个模态的相位差随距离变化叠加后幅度出现周期性起伏这就是所谓的干涉条纹远距离段高阶模态衰减掉只剩少数低阶模态干涉减弱曲线变得平滑。把这两种成因混在一起看就会觉得曲线“不规律”分开看每一段都能对应到具体的模态行为。在 Matlab 里验证这一点可以逐阶画单模态的损耗曲线再画叠加后的总曲线对比起伏出现的位置是否和模态相位差的变化周期吻合。这个对比做一次比看十页公式都直观。3.2 柱面扩展与球面扩展损耗公式里的几何因子很多人习惯用球面扩展的 20log10(r) 去套水声传播但简正波在柱对称假设下几何扩展是柱面性质的表现为 1/sqrt(r)对应 10log10(r) 的衰减斜率。这就是为什么直接用球面扩展公式去校验 Kraken 结果会对不上——几何因子选错了。实际海洋中还有吸收、海底衰减等额外损失这些体现在水平波数的虚部里而不是单独加一个经验项。常见的实现思路是把几何扩展和介质吸收分开处理几何部分由 1/sqrt(k*r) 承担吸收部分由 k 的虚部承担。两者叠加才是完整的传播损耗。如果你的 Matlab 脚本里只写了 20log10(r) 而没处理 k 的虚部远距离损耗会明显偏小波形幅度也会偏大。3.3 从频域到损耗曲线单频与宽带的差别Kraken 一次求解对应一个频率。传播损耗曲线如果是单频的直接按上面的链路算即可如果是宽带需要对每个频点分别求解再合成或者用宽带简正波方法。这里容易踩的坑是拿单频结果去代表宽带损耗曲线的频率依赖被抹掉阵元波形也会丢失色散特征。水声信道的色散正是来自不同频率的模态群速度不同单频求解根本体现不出来。所以在 Matlab 里做宽带仿真时频率循环是绕不开的。每个频点解一次本征值得到该频点的模态和波数再在时域合成。计算量上去了但这是拿到真实波形的必要代价。4. 阵元接收波形模态叠加怎么变成时域信号4.1 频域场到时间序列逆傅里叶变换的那一步阵元接收波形本质上是把频域场解做逆傅里叶变换。每个频点得到的是该频率下的复声压乘以源谱后对所有频点做 IFFT才得到时域波形。这一步里源谱的形状直接决定波形的包络和持续时间源谱设成窄带波形就长而平滑设成宽带波形就短而尖锐。很多人波形画出来“不像”问题往往出在源谱没设对而不是传播计算错了。在 Matlab 里频域场通常存成一个频率×深度的矩阵IFFT 沿频率维做。注意频率采样间隔要满足时域不混叠否则波形会出现周期性重复看起来像“回声”其实是采样问题。4.2 模态群速度与时域到达结构不同模态的群速度不同它们在时域上到达接收点的时间就不同这就是多途结构的来源。阵元波形上看到的多个“包”往往对应不同模态或同一模态的不同群速度分支。把每个模态单独做 IFFT再在时域上标出到达时间就能和群速度曲线对上。这个验证做一次你对波形的理解会从“一条曲线”变成“一组有物理意义的到达”。群速度由水平波数对频率的导数决定数值上可以用相邻频点的波数差来近似。差分步长选得太大群速度会偏选得太小数值噪声会放大。常见的做法是取适中的频率间隔再对结果做平滑。4.3 阵元位置与幅度加权为什么换个阵元波形就变接收阵元在深度方向上的位置不同模态函数在该深度的取值就不同各模态的权重随之改变波形自然不同。这就是为什么同一个传播环境不同阵元的波形幅度和干涉条纹会有差异。阵元靠近模态函数的节点时该模态贡献接近零波形上对应的成分就消失。在 Matlab 里做阵元波形对比时建议固定源和频率只改接收深度索引逐点画出来。这样能直观看到模态函数零点对波形的影响也能帮你判断某个阵元波形“异常”到底是传播问题还是位置选在了节点附近。5. 避坑与排查Kraken 对接 Matlab 时最容易翻车的几处5.1 现象损耗曲线整体偏移一个常数原因归一化约定不一致解决统一量纲后再比较这是最常见的翻车点。Kraken 输出文件的归一化方式和你在 Matlab 里假设的不一致损耗曲线就会整体平移形状对但绝对值错。排查方法是先算一个已知解析解的简单环境对比绝对幅度确认归一化因子后再套到复杂环境。不要一上来就拿实测数据拟合先把量纲对齐。5.2 现象阵元波形出现周期性重复原因频率采样间隔过大导致时域混叠解决加密频率采样或加窗波形上出现不该有的“回声”先查频率采样。频率间隔对应的时域周期如果小于你关心的信号长度就会混叠。解决方法是减小频率间隔或者在 IFFT 前对频域数据加窗抑制截断效应。这个坑在宽带仿真里尤其常见因为频率点多了容易图省事跳着取。5.3 现象远距离损耗偏小原因只用了几何扩展漏了波数虚部的吸收解决确认 k 为复数并正确参与衰减如果 Matlab 脚本里把水平波数当成实数处理衰减就只剩几何扩展远距离损耗会明显偏小。检查方法是看 k 的虚部是否非零以及在 exp(1ikr) 里是否保留了虚部。有些读取脚本为了省事只取实部这就是血泪经验的来源。5.4 现象高阶模态贡献异常大原因网格太粗导致本征值求解偏差解决加密网格并检查收敛性网格分辨率不够时高阶模态的本征值误差大衰减算不准波形上高阶成分会异常突出。排查方法是做网格收敛测试网格加密一倍看损耗曲线和波形是否稳定。如果不稳定说明还没收敛继续加密或调整边界条件。5.5 现象波形相位与理论反向原因波数符号或傅里叶变换约定不一致解决统一正负号约定并做单模态验证相位反向往往来自 exp 的符号约定或 IFFT 的方向。不同工具对正频率的定义可能差一个共轭。排查方法是取单模态、单频点手算相位随距离的变化和脚本输出对比。确认约定后在全流程里统一不要中途换约定。6. 设计取舍为什么简正波方法值得你花时间抠细节简正波方法的核心取舍在于它把三维传播问题降成了本征值问题加模态叠加计算效率高、物理图像清晰但前提是环境近似柱对称、水平变化不剧烈。如果你的场景是深海近场或强水平变化简正波假设会失效这时候耦合模态或抛物方程方法更合适。Kraken 的价值不在于它是万能解而在于它把“模态”这个中间层暴露给你让你能逐阶检查、逐阶验证。这种可解释性是很多黑箱求解器给不了的。在 Matlab 里对接 Kraken真正花时间的不是写循环而是把归一化、符号约定、几何因子这几件事对齐。对齐之后传播损耗和阵元波形就是同一套模态数据的两种视图改一个参数两边都能解释。我自己的习惯是每次换环境或换频段先跑一个单模态、单频点的最小案例确认幅度和相位都对得上再放开到全模态、宽带。这个习惯帮我省掉了大量“波形不对但不知道从哪查”的时间。从那以后我每次接新的传播输出都强制走一遍单模态验证再谈叠加。希望帮到你。本文还有配套的精品资源点击获取
