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

均匀线阵(ULA)CBF的波束图和方位估计(Direction of Arrival, DOA)空间谱的仿真程序。参考博文：https://blog.csdn.net/weixin_43713224/article/details/87181699

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‘[100300800300] ‘color‘‘w‘）
figure（1）
if （mode==1）
    polar（sitaR）
    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（sita20*log10（R/10^（-6））） %能量形式
    polar（sitashuchu） %能量形式
	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（sitaR）
    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（sita20*log10（R/10^（-6））） %能量形式
    polar（sitashuchu） %能量形式
	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（sitaR）
    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（sita20*log10（R/10^（-6））） %能量形式
    polar（sitashuchu） %能量形式
	title（ { ‘Figure 3‘ [‘d=‘ num2str（d） ‘ λ=‘num2str（nmd） ‘ N=‘ num2str（N） ] }‘FontSize‘9 ）
end

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件      247776  2018-11-02 10:13  Figures_of_Program.docx
     文件        1745  2018-11-02 10:08  Question_1.m
     文件         659  2018-11-02 10:04  Question_2.m