超声系统成像仿真与声场模拟程序

资 源 简 介

本项目提供一套专门用于超声系统成像仿真的MATLAB源代码。该程序的核心实现基于著名的Tupholme-Stepanishen方法，通过计算探头孔径的空间脉冲响应来精确模拟声场分布。为了解决传统解析方法在复杂场景下计算量巨大的问题，程序在实现过程中引入了远场近似技术，在保持物理特性的准确性的同时显著提升了运算速度，使其能够快速生成高质量的仿真成像结果。 该工具支持自定义多种超声换能器的几何参数，能够模拟不同激励信号下的声束形成过程。用户可以利用该源代码进行点散射体成像仿真、囊肿体模成像以及复杂组织结构的超声回波模拟。其应用场景非常广泛，包括医学超声成像算法的验证、新型换能器设计的性能评估、超声信号处理技术的研发以及物理声学原理的教学演示，是超声研究领域重要的仿真实验平台。

详 情 说 明

超声系统仿真程序项目说明

项目介绍

本项目提供一套用于超声系统成像仿真的 MATLAB 源代码。该程序旨在模拟超声波在介质中的传播、散射以及信号接收与处理过程。通过计算换能器孔径的空间分布和声场特性，程序能够准确生成射频（RF）信号，并最终合成高质量的 B 模式超声图像。其核心逻辑兼顾了物理真实感与计算效率，是医学成像研究、算法开发及教学演示的有力工具。

功能特性

1. 灵活的换能器配置：支持自定义线阵换能器的阵元数量、宽度、间距及发射聚焦深度。
2. 复杂体模构建：内置点散射体阵列与随机分布背景，并支持模拟无回声囊肿区域，反映真实组织特征。
3. 声场可视化：计算并展示换能器在空间中的发射声压分布情况。
4. 原始数据仿真：基于声学时延原理生成多通道原始射频信号，包含加性高斯白噪声模拟真实环境。
5. 全流程信号处理：实现从束向合成（延迟叠加）、包络检波到对数压缩的完整成像链条。
6. 标准 B 模式成像：输出符合医学超声标准的纵断面灰度图像。

使用方法

1. 环境配置：确保 MATLAB 环境已安装，建议版本为 R2016b 及以上。
2. 参数自定义：可以在脚本开头的参数设置区域修改声速、频率、采样率以及换能器几何参数。
3. 体模修改：通过调整散射体坐标矩阵或随机点生成逻辑，可以自定义特定的探测目标。
4. 运行仿真：直接在 MATLAB 中执行主程序脚本，程序将自动按顺序完成声场计算、信号生成、波束合成及图像输出。
5. 结果分析：程序运行结束后将弹出可视化窗口，展示声场分布图、RF 信号图及最终的 B 模式图像。

系统要求

1. 软件环境：MATLAB (包含 Signal Processing Toolbox 以支持 Hilbert 变换及相关信号处理功能)。
2. 硬件建议：为了保证快速生成包含大量随机散射体的仿真数据，建议配备 8GB 以上内存。

实现逻辑与算法说明

程序遵循标准的超声成像物理过程进行建模，具体步骤如下：

1. 系统与换能器初始化

设定介质声速为 1540m/s。定义线阵换能器的几何尺寸，计算阵元中心坐标。根据设定的聚焦深度，利用几何路径差计算每个阵元的发射聚焦延迟时延。

1. 激励脉冲生成

生成的激励信号。通过自定义函数实现高斯调制正弦脉冲（Gaussian Pulse），根据中心频率和分数带宽确定其波形，作为声场激发的原始振动信号。

1. 仿真体模设置

构造由离散点散射体和背景散射点组成的体模。通过逻辑判断（距离判定），在特定区域将散射强度设为零，从而模拟低回声或无回声的囊肿结构。

1. 发射声场分布计算

基于相位叠加原理，在远场近似条件下计算空间网格点的复声压。该过程考虑了各个阵元到空间点的传播距离以及预设的发射聚焦时延，展现了声束的聚焦效应和指向性。

1. 射频信号 (RF Data) 仿真

这是程序的核心计算部分。对于每一个接收阵元，程序遍历所有的散射体：

• 计算双程时延：发射波面到达散射体的时间（考虑聚焦延迟后的最短路径）加上散射体返回接收阵元的时间。
• 卷积叠加：根据总时延将激励脉冲波形按比例叠加到对应的时间点。
• 噪声模拟：在合成的信号中加入 30dB 的高斯白噪声。

1. 接收束向合成 (Beamforming)

采用延迟叠加（Delay-and-Sum, DAS）算法实现接收动态聚焦。针对每一条扫描线：

• 计算每个阵元在不同深度处的动态接收延迟。
• 利用线性插值从原始 RF 数据抽取对应时刻的信号并进行累加，从而增强焦点处的信号，抑制背景杂波。

1. 图像后处理

• 包络检波：对波束合成后的数据应用 Hilbert 变换，提取振幅包络，去除载波频率。
• 对数压缩：将提取的包络信号转换为分贝（dB）尺度，并限制动态范围（如 50dB），以适应人眼对亮度的对数感知需求。

1. 可视化输出

程序通过多子图联合展示：

• 发射声压空间分布热力图（dB 刻度）。
• 64 通道原始 RF 信号的时间-阵元编号图。
• 阵列中心位置的单一 RF 接收波形。
• 最终生成的 B 模式超声图像，纵坐标由时间轴映射为实际深度。