k-Wave仿真例程：创建超声换能器并绘制声场分布

k-Wave介绍

k-Wave软件是为了模拟超声波在1D、2D或3D中的传播。
应用示例包括：
- 均匀和非均匀介质中的传播
- 模拟各种类型的传感器
- 模拟多普勒效应
- 衍射、折射和反射
- 光声、超声成像
- 波束合成、成像重建
- 模拟弹性波

安装k-Wave

安装k-Wave需要几个步骤，可能会因您正在使用的操作系统而有所不同。以下是安装k-Wave的一般说明【注意：确保您的系统上安装了MATLAB，因为k-Wave是一个MATLAB工具箱】

1. 下载k-Wave：访问k-Wave官方网站(https://www.k-wave.org/)下载最新版本的软件。
2. 提取ZIP文件：下载后，将ZIP文件提取到您选择的文件夹中
3. MATLAB设置：打开MATLAB并添加k-Wave路径，打开Matlab并包含工具箱路径：＜文件＞＜设置路径＞＜添加文件夹＞＜保存＞
   

基本仿真流程

1. 定义：要计算的空间参数、传播介质参数、超声源、传感器
2. 运行模拟
3. 可视化结果
   

例程：利用k-Wave创建超声换能器

1. 设置 K-wave 网格

% 设置完美匹配层大小 perfectly matched layer (PML)
PML_X_SIZE = 20;            % [grid points]
PML_Y_SIZE = 10;            % [grid points]
PML_Z_SIZE = 10;            % [grid points]
% 设置不包括PML的网格点总数
Nx = 128 - 2*PML_X_SIZE;    % [grid points]
Ny = 128 - 2*PML_Y_SIZE;    % [grid points]
Nz = 64 - 2*PML_Z_SIZE;     % [grid points]
% 在不包括PML的x方向上设置所需的网格大小 ROI大小
x = 40e-3;                  % [m]
% 计算网格点之间的间距
dx = x/Nx;                  % [m]
dy = dx;                    % [m]
dz = dx;                    % [m]
% 创建网格 create the k-space grid
kgrid =