高斯声束点源模型Matlab代码

高斯声束点源模型通常用于声学领域，特别是声波扩散和声场模拟。在MATLAB中，你可以使用傅里叶变换和贝塞尔函数来近似描述这种声束。以下是一个简单的示例代码片段，展示如何创建一个二维空间内的高斯声束分布：

% 定义参数
fs = 44100;           % 采样率 (Hz)
f0 = 1000;            % 声源频率 (Hz)
sigma = 1;             % 高斯分布的标准偏差 (m)
r_max = 5;            % 计算范围的最大值 (m)

% 创建空间网格
[x, y] = meshgrid(-r_max:r_max, -r_max:r_max);
[r, theta] = polar(x, y);   % 转换到极坐标

% 计算时间延迟
t_delay = r / (speed_of_sound * f0); % 时间延迟，速度为预设常数(如空气中的声速)

% 构建高斯声束
声压 = exp(-((theta - pi/2)^2) / (2*sigma^2)) .* exp(-1i * 2*pi * f0 * t_delay); 

% 取实部作为声压信号
sound_pressure = real(声压);

% 如果需要，可以将声压转换回笛卡尔坐标并显示声场
sound_pressure_cartesian = beamform(sound_pressure, [x(:), y(:)]);
imagesc(x, y, sound_pressure_cartesian);
colorbar;
title('高斯声束在二维空间的分布');
xlabel('X轴 (m)');
ylabel('Y轴 (m)');

%