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

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

% 设置参数
f = 1000; % 频率
A = 1; % 幅度
phi = 0; % 相位
d = 0.1; % 声源间距
c = 343; % 声速
r = 1; % 观测点到声源的距离
t = 0:0.01:1; % 时间序列

% 计算声压级
p1 = A * sin(2*pi*f*t + phi);
p2 = A * sin(2*pi*f*t + phi + pi);
p = p1 + p2;
Lp = 20*log10(abs(p)/2e-5);

% 计算相位差
delta_phi = 2*pi*d*f/c;
phi1 = phi;
phi2 = phi + pi;
phi_diff = phi2 - phi1;

% 计算复合声场
p_r = A * (sin(2*pi*f*t + phi1) + sin(2*pi*f*t + phi2 + delta_phi*r))/r;
Lp_r = 20*log10(abs(p_r)/2e-5);

% 绘制动画
for i = 1:length(t)
    plot(r, Lp_r(i), 'ro', 'MarkerSize', 10);
    hold on;
    plot([0 r], [Lp(i) Lp_r(i)], 'b--');
    hold off;
    xlabel('Distance (m)');
    ylabel('Sound Pressure Level (dB)');
    title('Composite Sound Field');
    ylim([-10 10]);
    pause(0.01);
end

这段代码将生成一个动画，显示两个频率、幅度和相位都相同的点声源的复合声场。在动画中，红色圆圈表示观测点，蓝色虚线表示观测点到声源的距离。您可以根据需要调整参数，例如声源间距、观测点到声源的距离、频率、幅度和相位等。