clear;
%% 初始化
R=20; %声场半径(m)
d=0.033; %声场采样空间步长(m)
n=16; %扬声器个数(个)
nd=0.1; %扬声器之间的距离(m)
L=-R:d:R; %初始化空间
%% 通过SF函数求解声场
% SF250=SF(R,d,250,n,nd);
% SF500=SF(R,d,500,n,nd);
% SF1000=SF(R,d,1000,n,nd);
SF2000=SF(R,d,2000,n,nd);
% SF4000=SF(R,d,4000,n,nd);
%% 作图
% figure();
% mesh(L,L,SF250);
% colorbar;
% axis([1 20 -10 10]);
% title('阵列扬声器在频率250下的声压分布');
% figure();
% mesh(L,L,SF500);
% colorbar;
% axis([1 20 -10 10]);
% title('阵列扬声器在频率500下的声压分布');
% figure();
% mesh(L,L,SF1000);
% colorbar;
% axis([1 20 -10 10]);
% title('阵列扬声器在频率1000下的声压分布');
figure();
mesh(L,L,SF2000);
colorbar;
axis([1 20 -10 10]);
title('阵列扬声器在频率2000下的声压分布');
% figure();
% mesh(L,L,SF4000);
% colorbar;
% axis([1 20 -10 10]);
% title('阵列扬声器在频率4000下的声压分布');