声源距离图还是一片空白，重写

好的，下面是重写后的代码，主要修改了计算声源距离的部分：

% 设置声源和接收阵列参数
sourcePos = [0, 0, 0];  % 声源位置
arrayPos = [0, 0, 1; 0, 0, 0; 0, 0, -1];  % 阵列位置
fs = 10000;  % 采样频率
c = 1500;  % 声速
N = size(arrayPos, 1);  % 阵列元素数
f1 = 500;  % 第一个信号源频率
f2 = 550;  % 第二个信号源频率

% 生成信号源
t = 0:1/fs:1-1/fs;
s1 = sin(2*pi*f1*t);
s2 = sin(2*pi*f2*t);

% 传输信号源
d1 = pdist2(sourcePos, arrayPos, 'euclidean');  % 计算声源到各个阵列元素的距离
phi1 = 2*pi*f1/c*d1;  % 计算声源到各个阵列元素的相位差
y1 = s1.' * exp(-1j * phi1);  % 接收到的信号

d2 = pdist2(sourcePos, arrayPos, 'euclidean');  % 计算声源到各个阵列元素的距离
phi2 = 2*pi*f2/c*d2;  % 计算声源到各个阵列元素的相位差
y2 = s2.' * exp(-1j * phi2);  % 接收到的信号

% 计算相位差
dPhi = angle(y1) - angle(y2);

% 计算声源方向和距离
theta = atan2(sin(dPhi), cos(dPhi));  % 方向角
r = d1(1) * sin(theta) ./ sin(theta(1));  % 距离

% 绘制声源方向图
figure;
polar(theta, ones(size(theta)));  % 极坐标图
title('声源方向图');

% 绘制声源距离图
figure;
plot(t, r);
xlabel('时间（秒）');
ylabel('距离（米）');
title('声源距离图');

这样修改后，应该能够正确绘制出声源方向图和距离图。如果还有问题，请随时提出，我会尽力帮助您解决。