企业官网建设流程全解析

简介：Kraken 是声学领域经典的简正波声场计算程序，广泛用于水下声传播建模与海洋环境噪声分析。这份基于 MATLAB 的资源将 Kraken 声场模型的调用、数据读取与后处理流程整合为完整框架，共 11 个文件、约 482KB，包含 9 个功能 M 文件与 2 个已编译的 Windows 可执行程序。krakenc.exe 负责求解简正波本征值，field.exe 完成声场计算，read_shd.m、read_shd_asc.m 等接口可读取二进制/ASCII 格式的声场数据，支持声压、声强、传播损失等参数的提取与可视化；H05.m 与 runtests.m 给出了典型算例测试与结果校验流程，便于用户验证模型配置的正确性。在建模层面，资源覆盖了点源/线源设置、边界条件（硬边界、软边界）以及频域/时域分析等常用环节，用户可根据实际声学环境修改声速剖面、计算域与吸收系数等参数，直接获得声场分布图或传播损失曲线，用于建筑物内部声学设计或水下声呐性能预测。已有 795 人在线学习，适合具备一定 MATLAB 基础、从事水声学/声学工程或教学科研的工程师与研究人员。

1. 这个 kraken00_matlab 声场模型源码包，解决的是浅海声场计算的“最后一公里”

搞水声的同行看到kraken00_matlab_声场模型_声场_源码这种命名，第一反应大概率会联想到那个在水声领域流传多年的 KRAKEN 简正波模型。它解决的不是“声波怎么传播”这个物理问题，而是“你能不能在自己的工作区里快速拿到可靠的声场结果”这个工程问题。源码包把原本需要编译、配环境、命令行交互的简正波求解流程移植到 MATLAB 里，使读环境文件、算模式、输出传播损失全部在脚本内完成。适合两类人：一类是实验前需要快速预估传播损耗的从业者，另一类是想把声场计算嵌进更大仿真链路、不希望频繁切换编程环境的开发者。这篇直接拆解源码包怎么理解、怎么跑、参数怎么调，以及哪些坑值得绕开。

2. 从声场模型到 kraken00：为什么简正波算法是浅海环境的主角

2.1 简正波的基本思路：声场是有限个模式叠加出来的

在浅海这样的波导结构里，声源激发的声场并不是简单的一条直线加一次海底反射。上下两个界面（海面与海底）会迫使声波在垂直方向上形成驻波形态，只有满足特定条件的频率和入射角才能稳定存在，这些稳定的垂直驻波形态就是“简正波”。每一阶简正波对应一个水平传播的波数，携带一部分能量沿水平方向传播，随着距离递增逐渐衰减。接收点的声压是这些模式按幅度和相位叠加的结果。

常见实现里，求解步骤是：先根据声速剖面和环境参数计算每一阶模式对应的垂直波数和特征函数，再把各模式在声源深度上的激励强度、在接收深度上的值、以及沿水平方向的传播衰减组合起来。一个最小化的数值骨架长这样：

% 伪代码：凝固简正波求和的思想，具体函数以包内接口为准 for m = 1:nModes kz = root_m(m); % 第m阶垂直波数（由特征方程解出） psiS = modeShape(m, zSrc); % 模式在声源深度的值 psiR = modeShape(m, zRcv); % 模式在接收深度的值 km = sqrt(k0^2 - kz^2); % 水平波数 p = p + psiS*psiR/sqrt(km) * exp(1i*km*r - alpha_m*r); end

这段伪代码说明的是简正波求和的核心逻辑，不是某个具体函数的实现。实际计算中，特征值求解是最复杂的环节：声速剖面被分层近似，每一层内用解析形式表示场，层间用边界条件衔接，最后得到包含水平波数的特征方程。kraken00 这类源码包的价值，就是把这套求解过程封装成可在 MATLAB 里直接观察到中间变量的函数，而不是只给一个编译好的黑匣子。

从工程角度理解这一点，直接影响后面调参的态度。如果不知道声场是“若干模式叠加”，当传播损失曲线出现异常的起伏时，你很难判断到底是物理现象还是数值问题；知道模式数会随频率升高而增加，你就能预判高频计算为何变慢、为何容易碰到内存问题。

2.2 源码包里的数据流：环境文件、模式求解与声场输出

绝大多数 KRAKEN 移植版都保留了一套固定数据流，kraken00_matlab 也不会跳出这个框架。第一段是环境描述：水深、声速剖面、海底类型、声源和接收深度都写在一个文本文件里，不同实现可能叫.env或.flp，但作用一致。第二段是模式求解：读取环境文件后，程序对声速剖面做分层插值，在每层内得到特征函数的解析形式，再通过迭代搜索找到特征方程的根。第三段是声场合成：把模式信息在指定距离、深度网格上组合，输出复声压或传播损失。

拿到源码包后，我建议按照这个数据流去读文件目录。常见的分法是：环境文件读取函数、模式求解函数、声场合成函数、画图脚本和示例环境文件。读代码时不要从声场合成的入口函数开始翻，那样会被大量参数传递绕晕。先从“环境文件读取”看起，弄清楚它识别哪些字段，再去看模式求解输出的变量结构，最后回到声场合成时，你就能知道矩阵的每一维分别对应距离还是深度。

如果源码包里提供了示例脚本，先原样运行一遍，再逐行打断点看变量变化。这样做的收获比直接改参数大得多：你会亲眼看到环境文件里的几十字节文本，如何变成一个包含全部模式的复振幅矩阵。很多人在这一步跳过去，直接拿真实的海试数据往里塞，一旦结果反常，根本分不清是环境文件格式错了、模式数截断取值不合理，还是声场合成调用的参数顺序有误。

2.3 选型边界：简正波模型不是万能的，但它在水平分层环境中性价比最高

声场模型这十几年被反复比较，无非三类主流思路：射线模型适合高频、近距离、水平变化剧烈也不怕，因为它是完全的几何近似，追每条声线的路径就行，但聚焦区和影区的计算精度不够；抛物方程模型处理水平变化环境最稳，把波动方程在距离方向上一步一步推进，代价是计算量大，网格也要设计得很仔细；简正波模型正好在两者中间，对水平分层（声速剖面只随深度变化）环境最合适，低频、中等距离是它的舒适区。

kraken00 这类源码包被很多人拿来做浅海仿真，正是因为浅海我们可以近似看成水平分层：几十到一百多米的水深，声速剖面沿距离的变化不强，海底参数取一个区域平均值基本可用。这种情况下，简正波方法用较少的模式数就能获得足够精度，计算时间通常是秒级到分钟级。

那什么时候不要选它？如果模拟区域跨度很大，海底地形从几十米急剧变化到几百米，等深线变化剧烈，简正波方法的绝热近似容易失真。这时候强行用简正波模型，得到的场分布可能看起来合理，实际上已经偏离物理。另一个场景是宽带信号仿真需要连续跑上千个频点，模式求解的开销被线性放大，如果还想实时出结果，这类源码包就吃力了。

这个表不是说简正波不能用于那些场景，而是说计算代价和精度代价会迅速上升。选型时把它放在决策链路的早期，能省掉后面很多反复跑不收敛的麻烦。

3. 用 MATLAB 把 kraken00 跑起来：目录整理、环境文件与最小示例

3.1 源码包的目录体检与 MATLAB 路径设置

拿到一份声场模型源码包，先别急着运行，做三个体检动作：看目录里有没有examples文件夹，看有没有README或注释性说明，看入口脚本是m文件还是需要额外编译的mex文件。这三个动作基本能判断源码包是纯 MATLAB 移植版还是混合编译版。纯移植版改起来方便，能单步调试；混合编译版运行快，但断点进去只能看到编译后的函数声明，排查困难一点。

我习惯把源码包放在一个固定目录，然后通过一个独立脚本管理路径：

% setup_path.m：让源码包在任何工作目录下都能被找到 % 注意：把下面的路径改成你本机的实际位置 addpath(genpath('/home/user/acoustics/kraken00_matlab')); % 验证入口是否可见，MATLAB 会返回函数所在路径 which kraken00_run

genpath会把源码包下所有子目录递归加入搜索路径，避免漏掉藏在深层目录里的工具函数。运行后如果返回的是路径字符串，说明路径设置成功；如果返回'kraken00_run' not found，检查两点：路径里是否包含中文或特殊字符，源码包外层是否又多包了一层目录。很多新手在这一步翻车，不是代码问题，是 Windows 下MATLAB的默认工作目录和源码包不在同一个盘符，而addpath时用了相对路径。

路径就位后，先打开示例环境文件，把它和代码里环境读取函数对照着看。示例文件是理解格式的最快途径，比读任何注释都管用。确认示例能跑通之前，先别碰自己的数据。

3.2 环境文件的最小结构：先写一个能跑通的水下波导

环境文件是 kraken00 的输入源头，里面记录水深、声速剖面、海底半空间参数、声源与接收的布放位置。不同版本的移植包字段顺序可能不同，但一定包含这些核心信息。下面是一个最小环境文件的示意结构：

'Demo Shallow Water' 50.0 0.0 ! 水深(m) 与海面粗糙度参数 0.0 ! 声速剖面插值类型，0表示按给定点分段线性 1500.0 0.0 ! 深度0.0m 处声速1500 m/s 1505.0 25.0 ! 深度25.0m 处声速1505 m/s 1510.0 50.0 ! 深度50.0m 处声速1510 m/s 0 ! 海底以下介质类型，0表示流体半空间 1600.0 1.8 0.2 ! 海底声速(m/s)、密度(g/cm^3)、衰减(dB/波长)

这个结构是一个简化记忆，真实环境文件中还会有介质层数、层间界面、可选的上下边界类型等字段。第一次写环境文件时，拿它做骨架能快速验证整个链路；等程序跑通了，再根据源码包内模板把字段补全到完整形式。

写环境文件时最容易埋的雷是单位。有些移植版内部用米、秒、赫兹的 SI 单位制，有的保留了 Fortran 版本的旧习惯，深度和距离可能混用码或者英尺。我建议每次写完环境文件后，先加几百行fprintf或在 MATLAB 里读取一遍，把读入的数值打印出来核对，确认水深是 50 而不是 152.4。这一步只花五分钟，能省下后面几个小时的排障时间。

3.3 最小调用示例：获得第一条传播损失曲线

环境文件就绪后，调用入口函数获取声场。下面给出一个最简调用，函数名以你拿到的包为准，这里用kraken00_run做占位：

% run_demo.m：最小声场计算与传播损失出图 clear; close all; clc; addpath(genpath('/home/user/acoustics/kraken00_matlab')); envFile = 'demo.env'; % 由上个步骤生成的环境文件 freq = 200; % 计算频率 200 Hz zSrc = 10; % 声源深度 10 m zRcv = [5 25 45]; % 三个接收深度 rRcv = 0:10:2000; % 距离 0 到 2000 m，间隔 10 m % p 的结构一般是 nRcvDepth x nRange 的复数矩阵 p = kraken00_run(envFile, freq, zSrc, zRcv, rRcv); % 传播损失换算，注意取某个接收深度画图 TL = -20*log10(abs(p)); figure; plot(rRcv, TL(2, :), 'LineWidth', 1.5); % 取第二个接收深度 xlabel('Distance (m)'); ylabel('TL (dB)'); grid on; title(['f = ' num2str(freq) ' Hz, z_{rcv} = ' num2str(zRcv(2)) ' m']);

kraken00_run返回的矩阵维度由两个输入决定：接收深度个数与距离点数。代码里画的是 25 m 深度上的传播损失，这一层在浅海中通常能看出明显的干涉条纹结构。如果曲线在近场出现剧烈起伏，先看是不是距离采样间隔太大；如果整体趋势没有按距离衰减，检查海底衰减系数的量级是否填小了一个数量级。

传播损失换算成 dB 时注意参考值。很多实现返回的是相对于 1 m 处声压归一化后的复声压，所以直接用上述公式可以。有的移植版返回的不是声压而是强度，那里就要在先取绝对值还是先取平方上做调整。拿示例结果和源码包自带示例图对比一下，如果趋势完全不同，优先检查环境文件读入是否成功，再检查 zRcv 和 rRcv 的定义是否颠倒了。

4. 参数怎么设：声速剖面、海底边界、频率与网格的必调点

4.1 声速剖面的离散方式：分层不是越多越好

KRAKEN 类模型把声速剖面按深度离散成若干层，层内假设声速线性或者恒定，然后用边界条件逐层匹配。离散的层数直接决定特征值求解的精度和速度。层数太少，声速剖面的细节被抹掉，尤其是有温跃层时，特征值的计算结果会偏离真实值；层数太多，每一层的解析表达式都要参与匹配，求解矩阵变大，耗时随之上涨。

我的经验是，声速剖面的离散密度应当跟随声速梯度变化，而不是等深度间隔排列。梯度大的地方，比如温跃层附近，每 1 到 2 m 取一个点；梯度平缓的区域，每 5 到 10 m 取一个点就可以。如果用等间隔 1 m 全深度采样，算出来的结果不会更准，只会更慢。另一个实用技巧是，把实测声速剖面输进去之前，先做一次平滑滤波，剔除明显跳变的野值。野值会让层间边界条件出现伪模式，表现为某一阶模式的特征值反复迭代都不收敛。

源码包里通常会有声速剖面插值选项，默认可能是样条插值，也可能是分段线性。对于大多数浅海场景，分段线性够了；只有当声速剖面存在很浅的混合层、变化非常剧烈时，样条插值才值得考虑。调这个选项时，记得同步查一下模式画出来的垂直特征函数是否光滑，如果有明显折角或突变，说明插值方式选得不合适。

4.2 海底半空间参数：密度和衰减系数的影响远超你预期

海底是声波在浅海中最大的能量汇，海底参数设得不对，传播损失整体就会失真。常见的错误是从文献里抄一个数，却不管那个数的测量条件是否匹配自己模拟的频段。海底密度主要影响反射系数的大小，密度越高，反射损失越大；海底声速决定临界角，声速高，全内反射角度范围大，声能更容易被限制在水层中；衰减系数则决定每次海底反射损失掉的能量。三者共同影响远距离的传播损失斜率。

注意，表里的衰减系数在不同参考书里可能定义不一致，有的是 dB/波长，有的是 dB/m。如果源码包默认用 dB/波长，你在距离方向上按米计算衰减时，频率越高，每米衰减越大，这是符合物理的。反过来，如果把 dB/m 的数当成 dB/波长填进去，低频时衰减会被严重低估。拿到源码包后第一件事，就是确定这几个量的单位换算关系，而不是直接照抄别人的环境文件。

我通常的做法是：先用一组温和参数（细沙底、密度 1.8、衰减 0.5）跑出基准场图，再分别把密度、声速、衰减系数增大一倍，看传播损失曲线的变化趋势是否符合物理直觉。这个敏感性测试做完，你对自己环境参数的信心会完全不同。

4.3 频率、距离网格与模式数的三角关系

频率决定模式数，模式数决定计算量，也决定声场结构。浅海 50 m 水深、200 Hz 条件下，一般只有几个模式起作用；到了 2 kHz，模式数可能增长到几十个，内存和耗时都会上台阶。源码包通常会有一个“最大模式数”或“模式截断”参数，默认值未必匹配你的场景。

我一般按以下顺序设定网格参数：先根据频率算出波长 lambda，距离步长取 lambda/10 到 lambda/20，确保水平方向不欠采样；深度网格取每波长至少 5 到 10 个点，保证特征函数有足够分辨率。代码表示就是：

c0 = 1500; % 近似声速 m/s freq = 500; % 需要计算的频率 lambda = c0 / freq; % 距离步长：半波长以内都算安全，取更密一点避免干涉条纹失真 dr = lambda / 10; % 深度网格：从 0 到水深，按 10 层/波长划分 depth = 50; dz = lambda / 10; zGrid = 0:dz:depth;

这个示例只负责建立网格意识，实际距离向量和深度向量还要与接收点位置匹配。如果计算结果在距离轴上出现等间距的伪振荡，先怀疑距离步长太大；如果模式函数出现锯齿状，再怀疑深度网格太疏。

模式截断是另一个重要旋钮。截断模式数设置太高，会多算那些高掠射角的模式，它们衰减极快，对远场贡献很小，却平白增加计算时间；设置太低，会把真正有能力传播到远距离的低阶模式漏掉。调试办法是：把模式数从 5 逐步增加到 30，观察某几个远距离接收点上的声压幅值是否稳定。一旦稳定，说明已经收敛到足够多的模式。

4.4 用基准环境做参数校验：先让模型说“对”，再让它说“复杂”

参数调试最忌讳一开始就上完整实测剖面加复杂海底，那相当于一个方程组里同时改五个未知量，出了偏差也定位不到来源。我无论用哪个声场模型，都会先构造一个最简单的基准环境：等声速 1500 m/s，平海底，密度 1.5，衰减 0.1，然后对照可解析的传播规律。

等声速环境下，浅海远场传播损失大约按 10log10(r) 到 20log10(r) 之间的规律变化，近场则能看到典型的干涉结构。如果模型在这个环境下都给出离谱结果，那问题基本出在调用方式或环境文件读取上，而不是环境参数本身。基准环境跑通后，再逐步叠加真实声速剖面、海底分层结构、起伏地形，每次只改变一个变量，并记录传播损失曲线的变化。这样做的好处是，你能清楚知道每一条曲线特征是由哪个参数引起的，而不是看到一个复杂结果后凭感觉猜。

5. 跑 kraken00 声场模型的避坑指南：5 个常见问题与排查路径

5.1 现象：运行报错“未定义函数”或“找不到文件”

刚把源码包目录拷到新电脑上运行时，最容易遇到的就是入口函数报错为未定义。原因不一定是代码损坏，而是 MATLAB 路径中没有包含源码包目录，或者入口函数依赖的子函数藏在深层子目录中，addpath只加了根目录。

解决方法是把addpath换成addpath(genpath(...))，让全部子目录递归进入搜索路径。排查时先用which 函数名确认 MATLAB 到底能不能看到目标函数。如果显示file not found，再检查目录是不是被压缩软件解压成了双层嵌套结构，比如kraken00_matlab/kraken00_matlab，这种情况会把源码包的真实根目录藏深一层。还有一种可能：入口函数本身是.mex文件，换机器后需要重新编译。检查本地是否存在对应平台后缀的 mexw64 或 mexa64 文件，没有就找源码包里附带的编译脚本。

5.2 现象：传播损失曲线出现周期性尖峰或剧烈振荡

这种振荡通常不是物理干涉信号，而是数值欠采样。当距离步长过大时，模式之间的干涉图案不能被正确采样，表现为等间距的尖峰，间距正好和步长相关。另一个原因是接收深度太少，导致深度方向上声场结构被完全平均掉，画出来的曲线丢失了应有的细节。

解决方法是先降低频率，看看振荡是否消失。如果低频下曲线平滑、高频下振荡剧烈，说明问题确实出在网格密度上。把距离步长改为lambda/10以下，深度线上增加接收点。改完之后再看曲线，振荡通常会被平滑成连续的干涉起伏。值得注意的是，不是所有振荡都是数值问题，浅海中模式干涉本身就有位置相关性；区分办法是改变频率后振荡位置是否随波长等比移动，数值振荡往往在频率变化后保持同样的距离间隔。

5.3 现象：特征值求解阶段卡死或迭代不收敛

部分源码包在计算模式时会输出迭代过程，如果卡在特征值搜索环节，通常和初始猜测速度设置不当有关。许多实现的特征值求解需要一个初始的水平波数估计值或声速参考值，默认值在深海环境适用，但在浅海低速声速剖面下会落入错误的根区间，导致迭代发散。

解决方法是把初始参考声速改成整个剖面的平均声速附近。如果声速剖面范围是 1480 到 1520 m/s，参考值取 1500 通常是个安全起点。再不行，检查海底声速是否设置得低于水层声速，这种反演情况下模式数会急剧减少，甚至出现无实根，导致直接不收敛。调试时先把环境文件恢复到基准等声速状态，确认特征值求解能收敛，再逐项加入复杂参数。

5.4 现象：改了环境文件参数，结果却没有任何变化

这种情况是参数缓存导致的典型问题。有些源码包在同一个 MATLAB 会话中会缓存环境文件解析结果，或者用持久变量保存上一次计算的模式数据。你改了文本文件里的参数，但缓存没有失效，主函数直接读取的是上一次的解析结果。

解决方法是每次修改环境文件后，先执行clear all或至少清除与模型相关的全局变量和持久变量，再重新运行。更保险的做法是不要把环境文件参数与 MATLAB 工作区变量混用：环境文件里改，工作区里也改，那么代码里很可能有一个变量覆盖了文件里的值。排查时可以打印envFile读入后的结构体，逐字段比对是否与自己输入一致。一旦做过这个比对，就会发现很多“没生效”其实是你自己改了 A 参数，程序读的却是 B 参数。

5.5 现象：高频计算时内存溢出或运行时间不可接受

高频下模式数随频率近似线性增长，存储模式特征函数所需的内存快速增长。如果源码包默认保存全部模式，并且声场合成距离点数又非常大，内存爆掉并不奇怪。

解决思路是分层“断舍离”：先压低距离点数，算完保存结果，再清空中间变量；或者利用模式截断参数限制最高阶模式数。高频场景中，高阶级模式传播损失极大，对远场几乎没贡献。可以先用少数模式跑远距离趋势，近场需要细节时再增加模式数。如果源码包支持在距离轴上分块计算，就分段跑，每段之间把模式输出缓存在磁盘，避免内存同时堆积。

6. 进阶技巧：用批量参数扫描校验模型，而不是靠“跑通一次”就收工

声场模型源码包的价值不在跑通示例，而在你能把它当工具重复使用。我拿到任何声场模型，进阶做的第一件事就是写一个批量频率扫描脚本，用同一环境文件扫一组频率，观察传播损失随频率的变化是否连续、物理上是否自洽。

% scan_freq.m：频率扫描并保存结果，用于校验模型稳定性 clear; close all; clc; addpath(genpath('/home/user/acoustics/kraken00_matlab')); envFile = 'demo.env'; zSrc = 10; zRcv = 25; rRcv = 0:20:3000; freqList = 50:25:500; TL_all = zeros(length(freqList), length(rRcv)); for ii = 1:length(freqList) p = kraken00_run(envFile, freqList(ii), zSrc, [zRcv], rRcv); TL_all(ii, :) = -20*log10(abs(p(1, :))); fprintf('freq %d Hz done, median TL = %.2f dB\n', ... freqList(ii), median(TL_all(ii, :))); end % 把结果存成 mat，方便后续与不同海底参数的结果对比 save('scan_result.mat', 'freqList', 'rRcv', 'TL_all');

批量扫描时最有价值的检查点是低频端。浅海低频传播损失随距离变化的斜率对海底参数最敏感，如果扫描结果在低频段的衰减斜率随频率单调变化，说明参数模型是自洽的；如果出现回弹或剧烈抖动，多半是某一频率的模式截断出了问题。做完这一步，你基本可以判断这套源码包是否值得投入进实际项目。

我个人的习惯是，不管源码包来自哪里，都会先用一组已知参考解的环境做校验，再放心用在自己的场景里。有一回我拿到一套声称“直接可用”的声场模型，没有做基准校验就直接代入实验数据，结果海底衰减系数填错单位，导致整个距离段的传播损失低了快 10 dB，白白消耗了一周时间。后来把“基准验证”固化成流程：换任何环境参数前，先跑等声速浅海，再看曲线必须符合预估趋势。希望这个习惯能帮到你。

本文还有配套的精品资源，点击获取