超声声场线性描述与仿真分析工具

资 源 简 介

本工具旨在为超声成像系统的研究与开发提供一套完整的声场线性描述与仿真框架。核心功能基于线性声学理论，通过计算换能器的空间脉冲响应（Spatial Impulse Response），实现对复杂介质中超声波传播、衍射及散射过程的精确模拟。该工具支持自定义多种换能器探头拓扑结构，包括传统的线阵、凸阵以及用于深度成像的相控阵，允许用户灵活配置阵元尺寸、间距、中心频率、激励脉冲带宽及动态聚焦参数。 在仿真设计上，系统能够模拟不同声学特性（如非均匀声速分布、频率相关衰减）的介质环境，并支持导入随机分布的点散射体或精

详 情 说 明

超声声场线性描述与仿真工具

项目介绍

本工具是一个基于MATLAB开发的超声成像系统仿真框架。它建立在物理声学理论的基础上，通过数学建模模拟换能器产生声场、超声波在介质中传播、与散射体相互作用以及最终被接收并处理成图像的全过程。该工具主要用于科研及教学，旨在帮助研究人员深入理解超声成像的物理机制，验证波束合成算法，并评估系统在不同参数配置下的探测性能。

功能特性

1. 物理层建模：支持自定义声速、换能器中心频率、采样率及脉冲带宽等基础声学参数。
2. 换能器阵列配置：具备线阵换能器的几何描述能力，包括阵元数量、尺寸、间距及子元素剖分。
3. 复杂介质模拟：实现了考虑频率相关衰减和功率谱几何扩散的能量损耗模型。
4. 动态发射聚焦：模拟了基于时延的发射波束合成过程，允许设置特定的空间焦点。
5. 全流程信号处理：集成了从原始射频（RF）数据采集、Hilbert变换包络检波到对数压缩成像的完整链路。
6. 性能评估模块：自动计算点扩展函数（PSF）并提供横向分辨率的定性与定量分析。

使用方法

1. 配置环境：确保安装了支持信号处理函数的MATLAB版本。
2. 参数自定义：在代码开头的系统参数设置区域，根据实验需求调整声速、频率、阵元几何等变量。
3. 运行仿真：启动脚本，系统将依次执行脉冲生成、换能器布局、RF信号模拟及图像处理逻辑。
4. 结果查看：程序运行结束后将自动弹出结果窗口，显示阵元布局图、射频脉冲波形图、B模式超声模拟图像以及特定深度的点扩展函数曲线。

系统要求

1. 软件环境：MATLAB R2018b 或更高版本。
2. 工具箱：建议安装 Signal Processing Toolbox（用于Hilbert变换等信号处理操作）。
3. 硬件建议：由于仿真涉及多重嵌套循环计算（每个扫描线、每个散射体以及每个阵元的相互作用），建议配备 8GB 以上内存及多核处理器。

仿真实现逻辑说明

主程序逻辑严格遵循声学物理过程，分为以下关键阶段：

1. 空间脉冲响应近似：通过定义换能器阵元的子剖分结构，将每个物理阵元离散化，为后续基于线性瑞利积分思想的声场仿真打下基础。
2. 激励信号构建：利用高斯调制的正弦波生成满足预设带宽要求的激励脉冲。
3. 线阵扫描模拟：程序通过遍历每一个扫描线（Scan Line），计算换能器各阵元相对于目标的发射时延，以实现声学聚焦模拟。
4. 相干回波合成：这是核心计算模块。对于每一个点散射体，程序计算所有发射-接收对的往返总行程时间，并结合几何扩散带来的 $1/d^2$ 强度衰减和介质吸收的影响，通过相干叠加原理生成各通道的射频回波信号。
5. 信号调优与成像：模拟接收端对原始数据进行Hilbert变换提取瞬时包络，并应用对数增益补偿（Log Compression）和动态范围限制，将弱信号转换为适合视觉显示的灰度图像数据。

关键算法与实现细节分析

1. 延迟叠加（Delay-and-Sum, DAS）原理：代码在发射端实现了发射延时补偿。通过计算焦点到各阵元的几何距离差并将其转化为时间延迟，确保声波在预设焦点处达到同相叠加，增强了系统的侧向分辨率。
2. 衰减模型实现：工具采用了符合医学超声特性的衰减公式，考虑了频率（f0）与传播距离（dist_total）的乘积贡献。通过 $exp(-alpha cdot f cdot d)$ 实现幅度动态调整，使得深度散射体的信号强度表现更符合物理实际。
3. 包络检波算法：利用Hilbert变换将实值的射频（RF）信号转换为复解析信号，其模值即为声压包络。这一步去除了高频载波成分，保留了反应物体反射强度的低频包络。
4. 点扩展函数（PSF）分析：程序通过对图像中特定深度处散射体的横向切面进行采样并进行归一化对数变换，直观地展示了主瓣宽度和旁瓣水平。这是衡量波束合成精度和系统成像质量的核心指标。
5. 坐标系映射：实现了从采样时间轴到空间深度轴的转换（利用 $z = c cdot t / 2$ 关系），确保了生成图像在几何尺寸上的准确性，能够真实反映散射体在物体坐标系中的位置。