在MATLAB中在三维坐标系下设置模拟点声源的声压数据，设置阵列为121个阵列元素数量的正方形平面阵列，利用波束形成算法识别此点声源的代码，并绘制阵列图像、声源识别图像和声压值信息

以下是一个简单的 MATLAB 代码示例，用于模拟点声源的声压数据、设置阵列和进行波束形成：

% 设置模拟点声源的声压数据
f = 1000; % 声源频率
c = 343; % 声速
d = 0.5; % 声源到阵列的距离
theta_s = deg2rad(30); % 声源方向
theta = linspace(-pi/2, pi/2, 121); % 阵列元素方向
P_s = 1; % 声源声压

% 计算声压数据
r = d * cos(theta - theta_s);
P = P_s * exp(1i * 2 * pi * f/c * r) ./ r;

% 设置阵列
d_x = 0.05; % 阵列元素间距
x = -3*d_x : d_x : 3*d_x;
y = -3*d_x : d_x : 3*d_x;
[X, Y] = meshgrid(x, y);
Z = zeros(size(X));
N = numel(X);

% 进行波束形成
theta_scan = linspace(-pi/2, pi/2, 181); % 扫描方向
P_scan = zeros(size(theta_scan));
for i = 1:numel(theta_scan)
    theta_i = theta - theta_scan(i);
    r_i = d * cos(theta_i - theta_s);
    P_i = P_s * exp(1i * 2 * pi * f/c * r_i) ./ r_i;
    w = conj(P_i) ./ sum(abs(P_i).^2);
    P_scan(i) = abs(sum(w .* P))^2;
end
P_scan = P_scan / max(P_scan);

% 绘制阵列图像
figure;
scatter3(X(:), Y(:), Z(:), 'filled');
axis equal;
xlabel('x (m)');
ylabel('y (m)');
title('阵列图像');

% 绘制声源识别图像
figure;
polarplot(theta_scan, P_scan);
title('声源识别图像');

% 绘制声压值信息
figure;
surf(X, Y, real(P));
xlabel('x (m)');
ylabel('y (m)');
zlabel('P (Pa)');
title('声压值信息');

这段代码将产生三个图像：阵列图像、声源识别图像和声压值信息图像。其中，阵列图像显示了阵列的布局，声源识别图像显示了识别出的声源方向和强度，声压值信息图像显示了在阵列上各个点的声压值。