用Matlab写一段代码，可以用动画表示两个频率，幅度，相位都相同的点声源的复合声场

以下是一个简单的 Matlab 代码，可以用动画表示两个频率、幅度和相位都相同的点声源的复合声场：

% 设置参数
f = 1000; % 频率
A = 1; % 幅度
phi = pi/2; % 相位差
d = 0.1; % 点源间距
c = 343; % 声速
t = 0:0.0001:0.02; % 时间

% 计算声场
x1 = d/2 * cos(2*pi*f*t);
x2 = -d/2 * cos(2*pi*f*t + phi);
r1 = sqrt((x1.^2 + d^2/4));
r2 = sqrt((x2.^2 + d^2/4));
p1 = A ./ r1;
p2 = A ./ r2;
p = p1 + p2;

% 绘制动画
for i = 1:length(t)
    plot([x1(i), x2(i)], [0, 0], 'o', 'MarkerSize', 10, 'MarkerFaceColor', 'r');
    hold on;
    plot([x1(i), x2(i)], [p(i), p(i)], 'k', 'LineWidth', 2);
    hold off;
    xlim([-d/2-0.1, d/2+0.1]);
    ylim([-2*A, 2*A]);
    xlabel('x');
    ylabel('p');
    title(sprintf('t = %.4f s', t(i)));
    drawnow;
end

这段代码会生成一个动画，显示两个点声源在复合声场中的振幅随时间的变化。其中，`f`、`A`、`phi`、`d`、`c`、`t`分别表示声源的频率、幅度、相位差、间距、声速和时间。代码中的计算公式基于点声源的叠加原理，可以用于模拟任意数量、位置和振幅的点声源的复合声场。