clc;close all;clear;
mode = 1; % mode1显示声压谱,mode2显示能量谱
d=1; f=300; nmd=1500/f; N=10;
sita=0:0.001:2*pi;
R=abs( sin(N*pi*d.*sin(sita)./nmd)./N./sin(pi*d.*sin(sita)./nmd) );
% set(gcf,'Position',[100,300,800,300], 'color','w')
figure(1)
if (mode==1)
polar(sita,R)
title( { 'Figure 1', ['d=', num2str(d), ', λ=',num2str(nmd), ', N=', num2str(N) ] },'FontSize',9 )
end
if (mode==2)
shuchu = 20*log10(R/10^(-6))/max(20*log10(R/10^(-6)));
% polar(sita,20*log10(R/10^(-6))) %能量形式
polar(sita,shuchu) %能量形式
title( { 'Figure 1', ['d=', num2str(d), ', λ=',num2str(nmd), ', N=', num2str(N) ] },'FontSize',9 )
end
d=2; f=300; nmd=1500/f; N=10;
sita=0:0.001:2*pi;
R=abs( sin(N*pi*d.*sin(sita)./nmd)./N./sin(pi*d.*sin(sita)./nmd) );
figure(2)
if (mode==1)
polar(sita,R)
title( { 'Figure 2', ['d=', num2str(d), ', λ=',num2str(nmd), ', N=', num2str(N) ] },'FontSize',9 )
end
if (mode==2)
shuchu = 20*log10(R/10^(-6))/max(20*log10(R/10^(-6)));
% polar(sita,20*log10(R/10^(-6))) %能量形式
polar(sita,shuchu) %能量形式
title( { 'Figure 2', ['d=', num2str(d), ', λ=',num2str(nmd), ', N=', num2str(N) ] },'FontSize',9 )
end
d=3; f=300; nmd=1500/f; N=10;
sita=0:0.001:2*pi;
R=abs( sin(N*pi*d.*sin(sita)./nmd)./N./sin(pi*d.*sin(sita)./nmd) );
figure(3)
if (mode==1)
polar(sita,R)
title( { 'Figure 3', ['d=', num2str(d), ', λ=',num2str(nmd), ', N=', num2str(N) ] },'FontSize',9 )
end
if (mode==2)
shuchu = 20*log10(R/10^(-6))/max(20*log10(R/10^(-6)));
% polar(sita,20*log10(R/10^(-6))) %能量形式
polar(sita,shuchu) %能量形式
title( { 'Figure 3', ['d=', num2str(d), ', λ=',num2str(nmd), ', N=', num2str(N) ] },'FontSize',9 )
end