声呐系统中波束形成的MALTAB实现

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