上述代码水听器采用二维坐标修改

好的，下面是修改后的代码，这里假设水听器的坐标为二维坐标(x1,y1)和(x2,y2)：

% 生成信号
fs = 1000; % 采样率
t = 0:1/fs:1; % 时间范围
f = 10; % 信号频率
s = sin(2*pi*f*t); % 信号

% 生成两个不同位置的水听器接收到的信号
x1 = 0.1; % 位置1的x坐标
y1 = 0.2; % 位置1的y坐标
x2 = 0.2; % 位置2的x坐标
y2 = 0.3; % 位置2的y坐标
c = 1500; % 声速
d1 = [x1;y1]/norm([x1;y1]); % 位置1到信号源的单位向量
d2 = [x2;y2]/norm([x2;y2]); % 位置2到信号源的单位向量
r1 = norm([x1;y1]); % 位置1到信号源的距离
r2 = norm([x2;y2]); % 位置2到信号源的距离
t1 = r1/c; % 位置1接收到信号的时间
t2 = r2/c; % 位置2接收到信号的时间
s1 = s.*exp(-1i*2*pi*f*(t-t1)); % 位置1接收到的信号
s2 = s.*exp(-1i*2*pi*f*(t-t2)); % 位置2接收到的信号

% 计算相位差
phi = angle(s1./s2);

% 绘制相位差随时间的变化
plot(t,phi);
xlabel('时间（秒）');
ylabel('相位差（弧度）');

这段代码和之前的代码类似，只是将水听器的坐标从一维改为了二维。计算单位向量和距离的方式也做了相应的修改。最后绘制出的曲线仍然是相位差随时间的变化。