做超声阵列信号处理的人应该都绕不开 Field II 这个工具包。它是丹麦技术大学 Jørgen Arendt Jensen 教授团队开发和维护的开源超声声场仿真软件运行在 MATLAB 和 C 语言环境中用来计算换能器产生的声场分布、脉冲回波信号、波束特性等。我最早接触 Field II 是研究生阶段做阵列成像算法验证那时候实验室设备排期紧张全靠它先把换能器布局、聚焦参数、发射接收流程“跑”一遍把算法里的坑提前踩掉。后来工作了在换能器预研项目里做声场评估Field II 依然是我最常用的验证工具。可以说凡是接触过超声成像和阵列信号处理的人迟早都要跟它打交道。这篇博文是我完整捋过一遍 Field II 之后的整理笔记会从学习者的视角讲清楚三件事Field II 到底能做什么、它的核心仿真原理是什么、怎么从零开始搭一个可运行的线性阵列仿真以及我在实际调试中碰到过的问题和优化经验。适合刚接触超声仿真、还没跑通第一个例程的同学也适合已经能跑基础仿真、但觉得计算太慢或者结果看不明白的工程师。我尽量把每个关键步骤都写透包括参数为什么这么选、代码里每行在干什么、结果怎么解读。先纠正一个小误区Field II 里面的 II 是罗马数字 2正确读法是 Field Two不是 Field 二型 或者 FSAAII。这个命名是为了纪念早期研究基础同时也和产品名区分开大家交流的时候念 Field Two 就行。接下来我们一步步往下走。1. Field II是什么超声仿真领域绕不开的核心工具1.1 它的来头和定位Field II 的源头可以追溯到上世纪 80 年代末到 90 年代初Jørgen Arendt Jensen 当时在丹麦技术大学DTU做超声成像方面的研究写这套工具的核心目标是提供一个灵活的、基于线性系统理论的计算框架让研究人员不需要自己从底层推导每个换能器几何的声场公式而是靠通用的空间冲激响应法就能算出任意形状、任意激励下的声场和回波信号。它的底层实现是 C 语言通过 MATLAB 的 mex 接口暴露给用户所以既有 C 的计算速度又有 MATLAB 的矩阵操作便利性。从最早发布到现在Field II 保持了比较稳定的 API 设计你翻十年前的一篇论文里面的脚本放到当前版本上基本还能跑通。这一点对科研场景非常重要因为学术复现最怕的就是工具版本跨度太大导致结果对不上。它的学术认可度也相当高。在超声成像、波束合成、换能器设计这些方向上大量论文里的仿真结果都是用 Field II 跑的引用率非常可观。如果你准备写这个方向的论文审稿人看到仿真工具是 Field II基本不会质疑仿真结果的可靠性。1.2 它能做什么解决什么问题我个人的理解Field II 能解决的核心问题可以分成四块第一声场计算。给定一个换能器几何、激励信号、聚焦参数可以算出空间任意一点的声压时间波形。这个功能可以用于观察波束形状、焦点位置、旁瓣水平也可以用来验证相控阵的聚焦延迟算得对不对。第二回波信号仿真。在声场里放置若干散射点每个点都会产生一个回波Field II 会根据发射和接收两个过程的声场来计算最终的射频回波信号。这就是模拟超声成像的基础。你可以把散射点布置成血管截面、组织边界也可以布置成均匀分布的散斑用来模拟组织的背向散射。第三阵列设计与参数验证。阵元数量、阵元宽度、阵元间距、阵元高度、孔径大小、有无 apodization幅度加权这些参数对波束质量的影响都可以在仿真中定量地观察。实测工作中换能器一旦做出来就很难改参数仿真阶段多验证几次能避免很多返工。第四血流与运动估计。Field II 可以模拟血流方向上的运动散射体并产生一系列回波信号用于验证血流速度估计算法。这也是它在超声多普勒研究里被大量使用的原因。当然Field II 不是万能的。它基于线性声学假设本身不处理非线性传播、也不模拟组织中的衰减和多次散射这些复杂效应。真要研究谐波成像或者高强度聚焦超声的非线性声场还要结合其他工具。但就“超声阵列系统级仿真”这个定位来说Field II 在精度和效率之间的平衡确实做得最好。1.3 它和热词里其他仿真工具的区别很多同学在搜索时会把 Field II 和 k-Wave、COMSOL、ANSYS 这些仿真工具混在一起其实定位完全不同。k-Wave 更偏重声波在介质中的传播可以模拟非线性效应、不均匀介质中的散射但它的输入通常是一个已经定义好的初始声场源而不是一个精细的换能器阵列。Field II 更像“换能器机制造商 成像系统前端”k-Wave 更像“介质传播实验室”。实际使用中有人会把 Field II 算出来的声场作为 k-Wave 的初始源接力做非线性传播仿真这个组合在超声研究里非常流行。COMSOL 和 ANSYS 这类有限元工具胜在能算非常复杂的换能器结构力学、压电耦合、热效应但代价是计算量大、建模周期长而且很难直接输出成像意义上的射频回波信号。Field II 的优势是快一个 64 阵元的线性阵列聚焦声场普通电脑上跑到毫秒到秒级就能出结果适合做参数扫描和算法迭代。明白了这些区别你就能理解为什么 Field II 在超声阵列信号处理这个细分领域几乎是“标配”。它不管压电材料怎么振动不管非线性怎么演化它只管从换能器表面振动到接收信号这一条线性链路上的事而这一条链路恰恰是阵列成像算法验证最关心的地方。2. 核心原理空间冲激响应为什么能“算”出声场2.1 线性系统视角理解超声声场很多人第一次用 Field II直接把函数一跑、图一画看起来挺好但完全不知道内部发生了什么。这样用久了容易出问题因为一旦结果异常根本不知道是参数设置错误还是物理上本来如此。我比较推荐先用“线性系统”这个角度去理解。想象一下你在空旷房间里拍一下巴掌房间各个位置听到的声音其实是“拍巴掌这个信号”和“房间声学冲激响应”卷积的结果。房间不同位置的冲激响应不同所以你听到的回声和混响也不一样。超声声场本质上一个道理换能器表面振动产生的声压可以看作换能器表面振动速度与空间某点的声学冲激响应的卷积。Field II 使用的正是这个思路。它把换能器看作一个发射面源表面每个点都被同一个电激励信号驱动但根据几何位置不同对空间某一点的声压贡献有时间差和幅值差。把这些贡献做叠加就得到该点的总声压时间波形。2.2 空间冲激响应法的数学直觉在数学上对于线性、均匀、无耗散的声学介质换能器表面振动速度 ( u(t) ) 在空间某点产生的声压可以写成[ p(\mathbf{r}, t) \rho \cdot u(t) * h(\mathbf{r}, t) ]其中 ( \rho ) 是介质密度( h(\mathbf{r}, t) ) 就是空间冲激响应。这个冲激响应描述了如果换能器表面被一个理想的脉冲速度激励空间点 ( \mathbf{r} ) 处会接收到什么样的压力变化。一旦我们得到 ( h(\mathbf{r}, t) )再和任意激励信号卷积就可以算出任意激励下的声压波形。这就像线性电路里得到冲激响应后任意输入信号都能通过卷积得到输出原理是同一套。Field II 计算 ( h(\mathbf{r}, t) ) 的方法是基于 Tupholme 和 Stepanishen 等人的经典工作——把矩形换能器表面分成若干个小单元对每个小单元计算它在空间点产生的冲激响应再按几何时间延迟叠加起来。这与惠更斯原理的精神一脉相承每个小面积元都是一个新的次级声源空间某一点听到的声压是所有小面积元在当前时刻贡献的总和。2.3 为什么换能器会被切分成很多小元Field II 里的换能器并不是一个整体矩形就完事它会自动把每个阵元在长度和宽度方向划分成更小的子单元内部管这叫 mathematical elements。为什么这么切因为空间冲激响应的解析表达式对复杂形状换能器很难直接写出来但切成足够小的矩形小元之后每个小元可以近似为点源或者简单矩形源贡献叠加起来就有很高的精度。这里有一个核心代价权衡子单元划分得越细声场计算精度越高但计算量会显著上升。Field II 提供了一些参数来控制这种划分比如xdc_set里的no_ele_w和no_ele_h分别控制每个阵元在宽度方向和高度方向的划分倍数。如果你仿真的阵列频率很高、孔径很大默认划分可能不够需要手动调大但如果随便调得太大计算时间可能直接翻几倍甚至几十倍。用生活化一点的类比这就像画一幅马赛克画砖块越小图像细节越清晰但拼起来的工作量越大。你得在“够看”和“时间成本”之间找平衡。这也是我在后面优化部分会重点讲的内容。3. 环境搭建从下载到验证的一整套流程3.1 获取软件与目录准备Field II 的官方下载地址在 DTU 的主页上搜索 Field II 就能找到不需要注册下载下来是一个压缩包解压后里面有Field_II_Matlab_...这样的目录里面包含全部 m 文件和编译好的 mex 文件源码。C 语言版本也有但大多数用户用 MATLAB 版本就够了。解压之后我的习惯是在 MATLAB 里单独建一个工作目录比如D:\ultrasound_sim\把 Field II 目录放进去仿真脚本也放进去保持目录整洁。这一步看似简单但非常重要因为 Field II 运行时会在当前目录生成一些临时文件比如fieldII.bin目录混乱容易导致初始化失败。注意千万不要把 Field II 目录放在中文路径下。MATLAB 对中文路径的兼容性虽然一直在改善但 mex 文件加载和临时文件写入层面有时候还是会踩坑出现莫名其妙的初始化失败。用全英文路径最稳妥。3.2 在 MATLAB 中设置路径与编译环境启动 MATLAB 后第一步不是直接敲代码而是先把 Field II 目录加到 MATLAB 搜索路径里。可以用addpath命令也可以右键文件夹添加到路径。具体如下addpath(D:\ultrasound_sim\Field_II_MATLAB_3_30);这里要注意Field II的版本号会变目录名以你解压出来的实际名字为准。添加完路径后还要确认 MATLAB 的 mex 编译环境可用。现代 MATLAB 自带编译器配置但偶尔也会出问题。在命令行里跑一句mex -setup会弹出一个窗口让你选择编译器Windows 上通常选择 minGW 或者微软的 C/C 编译器。如果没有可用的编译器需要先安装 MATLAB 支持的编译器。这一步不能跳过因为 Field II 的底层 mex 文件如果平台不匹配就需要用你本地的编译器现场编译。配置完成后在 Field II 目录下运行make_field或对应的编译脚本它会重新编译所有 mex 源文件。官方发布的压缩包里会带各个平台预编译好的 mex 文件正常情况直接能用但如果 MATLAB 版本太新或者平台匹配不上重新编译就是必要的。3.3 验证安装是否成功路径和编译都没问题后可以跑一个最简初始化流程验证一下field_init(-1); field_info;field_init表示初始化 Field II 环境参数 -1 表示关闭写盘功能只做内存计算。正常情况会输出类似Field II simulation initialized的信息。field_info会打印 Field II 的版本和基本信息。然后跑一个最朴素的例子做一个单阵元换能器计算轴上一点的脉冲声场画出波形。如果这一步能出图说明环境完全可用。fs 100e6; field_init(-1); Th xdc_linear_array(1, 0.5e-3, 10e-3, 0); t_exc 0 : 1/fs : 2e-6; xdc_excitation(Th, sin(2*pi*1e6*t_exc)); points [0; 0; 50e-3]; [hp, tstart] calc_hp(Th, points); t_axis tstart (0:length(hp{1})-1)/fs; plot(t_axis*1e6, hp{1}); xlabel(时间 (μs)); ylabel(声压 (任意单位));这个脚本如果顺利运行会看到一条类似衰减振荡的波形。至此环境搭建结束可以进入下一步正儿八经的建模。4. 核心函数与建模方法4.1 换能器怎么描述xdc 系列命令Field II 里换能器建模用的主要是xdc_开头的一系列函数常见的几种xdc_linear_array创建线性阵列xdc_focused_array创建球面聚焦阵列xdc_concave创建凹面换能器xdc_rectangles用任意矩形单元拼接换能器xdc_apodization设置幅度加权以最常用的xdc_linear_array为例它的完整调用格式是Th xdc_linear_array(no_elements, width, height, kerf, no_ele_x, no_ele_y, no_sub_x, no_sub_y, focus);no_elements是阵元数量width是阵元宽度方向尺寸沿 x 轴height是阵元高度沿 y 轴kerf是相邻阵元间的间隙。阵元中心间距等于width kerf。后面四个参数是物理划分和数学划分的倍数物理划分和数学划分是用于提高计算精度的内部设定最后的focus是几何聚焦点可以暂时为空。在这个函数里有一个高频踩坑点阵元数量和阵元间距直接决定阵列的总孔径。比如 64 阵元、间距等于一个波长那么总孔径大约是 64 个波长。孔径越大波束越窄但旁瓣也越多越复杂。很多人一上来就抄论文里的 128 阵元结果算出来栅瓣一堆就是因为间距超过了半波长太多。4.2 激励信号与脉冲响应换能器的振动不是瞬间完成的它由两个时间信号决定第一个是电激励信号用xdc_excitation设置。这是你给换能器加的电压或驱动脉冲在 Field II 里表示为一个离散时间序列。比如中心频率 5 MHz 的换能器激励通常用 1 到 2 个周期的正弦波采样率设 100 MHz就是 1.5 个周期内采样 30 个点。第二个是换能器自身的脉冲响应用xdc_impulse设置。它描述的是换能器把电信号转换为声信号的响应特性以及接收模式下声信号转电信号的响应。简单建模时可以直接设为一个单位幅度、长度很短的脉冲相当于理想宽带换能器想模拟实际换能器带宽时可以设为中心频率附近的带通脉冲。这两个信号在物理上是卷积关系。实际发射声场的时间波形是激励信号、换能器脉冲响应和空间冲激响应三者共同卷积的结果。理解这一点对后续看波形结果很有帮助如果你仿真出来的波形包络比预期宽很多不要慌先看激励信号和脉冲响应是不是设置得太宽频带。4.3 聚焦、偏转与延迟阵列换能器的核心能力之一就是电子聚焦和波束偏转。Field II 中控制聚焦和偏转的方式有三种xdc_focus简单指定空间焦点Field II 自动计算各阵元延迟xdc_focus_times手动设置每个阵元的延迟时间灵活性更高xdc_dynamic_focus接收模式下动态聚焦使用xdc_focus是最省事的。比如想聚焦到 (0, 0, 50mm)只需xdc_focus(Th, 0, [0 0 50e-3]);第二个参数 0 表示从当前时刻开始聚焦。Field II 内部会计算每个阵元到焦点的传播时间差自动加上对应的延迟。这个时间差就是“聚焦延迟”它的物理意义是较远阵元先激励较近阵元后激励让所有阵元发出的声波能同时到达焦点。如果你想手动确认这个延迟算得对不对可以用xdc_focus_times自己算一遍pitch width kerf; x_i ((0:no_elements-1) - (no_elements-1)/2) * pitch; dist sqrt(x_i.^2 focus_z^2); delay (max(dist) - dist) / c; xdc_focus_times(Th, 0, delay);这个手动计算的过程对理解阵列波束合成非常有帮助我建议就算用xdc_focus也抽出十分钟手动验证一次。4.4 声场计算三兄弟场计算有几个函数新手最容易弄混calc_hp计算脉冲声场也就是在给定激励下的声压时间波形calc_h计算空间冲激响应即前面原理部分提到的 ( h(\mathbf{r}, t) )不包含激励信号calc_scat计算散射点回波是发射和接收双向过程的总和实际使用中calc_hp用得最多。它的输入除了换能器句柄还要给一个 3×N 的坐标矩阵每一列是一个空间点的 ( (x, y, z) ) 坐标。输出hp是 N×1 的 cell 数组每个 cell 是该点的声压时间序列另一个输出tstart是该点的时间起点。这里有一个非常影响效率的点calc_hp的坐标矩阵是逐列处理的如果你把几万个点一次性传进去内部会循环计算内存占用和耗时都会直线上升。我一般会把问题拆成多个小网格分别计算比如先算轴向再算横向最后用cat拼起来。这比一次性把所有网格点塞进去要稳得多。4.5 控制计算精度与速度的取舍Field II 的计算精度主要由两个因素控制阵元内部的子单元划分数量以及换能器表面边缘的采样密度。默认参数对低频、小孔径的换能器够用但如果你仿真高频阵列、大尺寸阵元声场结果可能会出现微小的伪波纹或数值噪声。这时候就要调高no_sub_x和no_sub_y也就是每个阵元在宽度和高度方向的划分数量。我个人的经验法则是先保持默认参数跑一版结果然后把划分数量翻倍再跑一次对比声场曲线如果两次结果差异很小说明当前划分精度够了如果曲线明显变化说明精度不足需要继续加细。这种方法比凭借经验乱调参数高效得多也能帮你对仿真结果的可靠性建立信心。5. 实操线性阵列聚焦声场仿真全流程5.1 场景设定与参数选择这一节我们来复现一个非常典型的场景一个 64 阵元的线性阵列中心频率 5 MHz阵元间距为一个波长焦点设在轴线上 50 mm 处观察焦点深度上的横向波束分布。参数为什么这么选我逐一解释中心频率 5 MHz医学超声里很常用的频段兼具分辨率和穿透深度实验数据也比较好找。64 阵元常用相控阵/线性阵列规模不多不少计算量也比较适中。阵元间距等于波长这是波束合成的经典约束。间距小于半波长可以避免栅瓣但阵元数不变时总孔径会变小波束变宽等于一个波长是折中方案在很多线性阵列应用中是可行的。焦点 50 mm中等深度既能体现聚焦增益又不会让近场效应干扰主瓣判断。声速取 1540 m/s这是软组织中声速的常用近似值。波长计算出来大约是 0.308 mm64 个阵元的阵列总孔径大约是 19.7 mm这个因径比孔径与焦点深度的比值约 0.4在超声成像中非常典型。5.2 完整代码与逐行注释下面这段代码可以直接复制到 MATLAB 里运行我已经把每一行的作用都写在注释里。%% Field II 线性阵列聚焦声场仿真 % 功能计算64阵元线性阵列在焦点深度处的横向声场 % 清理环境 clear; clc; % 初始化 Field II-1 表示不写盘 field_init(-1); % 基本参数 f0 5e6; % 中心频率 5 MHz fs 100e6; % 采样率 100 MHz c 1540; % 声速 m/s lambda c / f0; % 波长约 0.308 mm % 阵列参数 num_ele 64; % 阵元数量 width 0.95 * lambda; % 阵元宽度留 5% 的间隙 kerf 0.05 * lambda; % 阵元间隙 height 10 * lambda; % 阵元高度远大于宽度接近线源 focus_z 50e-3; % 焦点深度 50 mm % 创建发射换能器 th_emit xdc_linear_array(num_ele, width, height, kerf); % 设置激励信号1.5个周期的正弦波 t_exc 0 : 1/fs : 1.5/f0; exc sin(2 * pi * f0 * t_exc); xdc_excitation(th_emit, exc); % 设置理想宽带换能器脉冲响应 xdc_impulse(th_emit, 1); % 聚焦到 (0, 0, focus_z) xdc_focus(th_emit, 0, [0 0 focus_z]); %% 计算横向声场z 50 mmx 从 -15 mm 到 15 mm x linspace(-15e-3, 15e-3, 121); z focus_z * ones(size(x)); points [x; zeros(size(x)); z]; [hp, tstart] calc_hp(th_emit, points); % 提取每个点的声压包络峰值 hp_peak zeros(size(x)); for i 1:length(x) hp_peak(i) max(abs(hp{i})); end % 归一化并转成分贝 hp_db 20 * log10(hp_peak / max(hp_peak)); % 绘制横向声场 figure; plot(x * 1e3, hp_db, LineWidth, 1.5); xlabel(横向位置 x (mm)); ylabel(声压幅度 (dB)); title([焦点深度 , num2str(focus_z*1e3), mm 处的横向波束]); grid on; ylim([-60 0]); % 标注 -6dB 波束宽度 idx find(hp_db -6); half_width x(idx(end)) - x(idx(1)); text(-18, -10, [-6dB 宽度 , num2str(half_width*1e3, 3), mm]);运行这段代码你会看到一条类似高斯形状的波束曲线。主瓣在 x0 处最高两侧逐渐下降大约在 -6 dB 处的主瓣宽度值得记录下来这就是这个阵列在 50 mm 深度上的横向分辨率的近似值。5.3 轴向声场与焦点位置的验证只看横向声场还不够我还喜欢同时看轴向声场验证焦点是不是真的在 50 mm 附近。%% 计算轴向声场x 0z 从 10 mm 到 80 mm z_axis linspace(10e-3, 80e-3, 141); points_axis [zeros(size(z_axis)); zeros(size(z_axis)); z_axis]; [hp_axis, tstart_axis] calc_hp(th_emit, points_axis); hp_axis_peak zeros(size(z_axis)); for i 1:length(z_axis) hp_axis_peak(i) max(abs(hp_axis{i})); end % 归一化并转分贝 hp_axis_db 20 * log10(hp_axis_peak / max(hp_axis_peak)); figure; plot(z_axis * 1e3, hp_axis_db, LineWidth, 1.5); xlabel(轴向深度 z (mm)); ylabel(声压幅度 (dB)); title(轴向声场分布); grid on; ylim([-40 0]);轴向声场曲线通常在焦点附近出现一个明显的峰值峰值对应的深度应接近 50 mm。如果峰值位置偏了首先要检查声速设置是否正确其次检查xdc_focus的焦点坐标是否书写正确。这两个都是多月工作中高频出错的地方。5.4 参数变化对波束的影响跑通基础仿真之后我强烈建议做一个小实验改变阵元数量观察波束宽度和旁瓣水平的变化。方法是把num_ele从 32 改成 128重新运行脚本对比三条曲线。你会很直观地看到阵元越多主瓣越窄但旁瓣结构也可能更复杂。阵列孔径决定了波束聚焦能力孔径越大理论波束宽度越窄这就是相控阵成像提高分辨率的基本原理。再做一个实验把xdc_focus改成不聚焦即让所有阵元同时发射看看声场变成什么样。不聚焦时没有几何增益焦点深度处的横向声场会明显变宽峰值也低很多。这个对比能帮你直观体验“聚焦”在超声成像里的核心作用。6. 常见问题与调优经验6.1 初始化失败报 “Unable to open fieldII.bin”这个问题我遇到过很多次原因基本都是 MATLAB 当前工作目录不是 Field II 目录而 Field II 初始化时需要在当前目录写临时文件。解决办法有两种一种是每次运行前cd到 Field II 目录另一种是在脚本开头把工作目录临时切过去再切回来。如果用了field_init(-1)参数不写盘跑内纯计算问题通常也会消失。6.2 编译不通过mex 报错如果在make_field阶段报编译错误八成是 MATLAB 找不到合适的编译器。用mex -setup重新选择编译器Windows 上推荐 minGW-w64MAC 上注意 Xcode Command Line Tools 要装完整。还有一个小概率情况是 MATLAB 版本太新Field II 的 mex 源码里某个头文件写法与新版不兼容这时需要对照官方更新说明看看有没有提供针对新版 MATLAB 的补丁包或源码修复。6.3 计算慢到怀疑人生这是最影响使用体验的问题。我个人的优化顺序是降低采样率。很多场景不需要 100 MHz 采样如果换能器中心频率只有 3.5 MHz50 MHz 采样已经完全满足奈奎斯特条件计算量直接减半。减少计算网格点。先跑一个粗网格看趋势比如 41 个点再在波束峰值附近加细不要一开始就 201×201 个点全算。检查阵元划分参数。no_sub_x和no_sub_y如果从默认调大了很多乘以阵元数之后计算量会爆炸。尽量在精度验证之后再固化这一项。分支计算。把横向和轴向分开算比一次性算一个平面快很多因为每个方向的点数少。关闭不必要的写盘。使用field_init(-1)而不是默认的写盘模式。另外Field II 本身是单线程的MATLAB 启动并行池对calc_hp的帮助不大。我试过用parfor循环并行计算多个点的声场但xdc对象句柄在并行 worker 之间的传递经常出问题所以不太推荐除非你把计算任务拆成完全独立的子进程分别跑。6.4 声场结果出现 NaN 或 INFNaN 和 INF 几乎都是延迟时间设置异常导致的。比如xdc_focus_times里延迟时间向量长度不等于阵元数量或者焦点距离为 0都会导致除以零或者根号下负数的运算。排查顺序先检查阵元数量是否一致再检查焦点坐标是否有 0 或负值最后检查声速单位是否统一为 m/s。Field II 所有长度单位都用米如果漏写了e-3结果会离谱到无法直视。6.5 常见问题速查表现象可能原因排查与解决初始化失败找不到 fieldII.bin工作目录错误使用field_init(-1)或切换目录mex 编译报错编译器未配置mex -setup配置合适的 C 编译器计算结果全为 NaN延迟时间向量异常检查焦点坐标、阵元数量一致性波形包络宽于预期激励信号周期太多减少激励信号周期到 1~2 个计算耗时过长采样率高 / 网格密降采样、粗网预处理、分支计算旁瓣异常偏高阵元间距过大检查widthkerf是否接近或超过波长焦点峰值位置偏后声速设置偏小将声速改为 1540 m/s 或实际介质声速用实际项目经验收个尾我花在这份笔记上的时间其实很大一块不是在学习 API而是在反复调参和验证结果。Field II 是一个上手容易精通难的工具它的文档相对零散很多细节藏在一个个示例脚本里。我的建议是一定不要只盯着官网给的 example 看完就关而是改一个参数、看一个结果亲手把每个函数的行为摸清楚。我第一次用xdc_focus_times的时候因为没搞明白延迟时间的参考基准整整浪费了半天。还有一个小技巧值得分享Field II 仿真结果最好和理论值对比验证一下。比如对于均匀激励的矩形孔径横向波束应该符合 sinc 函数的形式主瓣宽度和旁瓣位置都有理论公式。如果仿真结果和理论预期明显不一致就先不要急着往下做系统仿真回头查参数设置直到两个结果吻合。这一步能帮你把对 Field II 的“黑盒感”彻底消除后续用起来心里才踏实。最后Field II 本身的生态很丰富学会了基础建模之后可以继续拓展到动态聚焦接收、多角度平面波成像、血流速度仿真这些更高级的方向。每一步拓展其实都是在现有知识树上加新枝。希望这份笔记能帮你把第一根枝条扎稳后面的事情就顺理成章了。