利用FFT在频域进行多目标波束形成

波束形成，可以简单的理解为阵列流形的加权，或者滤波也可以。本质上是对相位关系进行数学计算分析。
基本解释可以从这个图来看：

这里对理论不做过多的阐述，本质上利用波束形成进行DOA估计是判断阵列每个方向上来的能量。
每个角度上的能量累积形成波束方向图。

具体仿真如下：
%% 利用FFT在频域进行多目标波束形成 实际使用过程中需要对相位校准进行预处理 即对exp(-1jphase) phase为通道相位差
clc;clear;close all;
c=3e8;f0=10e6;fs=58.523528e6;%电磁波速度、中频中心频率、采样频率 也可利用带通采样进行波束形成
ome=2pif0;%角频率
d=c/f0/2;%阵元间距
N=2^10;%采样点数
M=16;%阵元个数
t1=(0:2N-1)/fs;
s0=cos(omet1);%参考信号
s0 = awgn(s0,20,10);%加高斯白噪声
theta0=50;%波束形成目标角度0
theta1=40;%波束形成目标角度1
theta2=30;%波束形成目标角度2
theta3=20;%波束形成目标角度3
theta4=10;%波束形成目标角度4
theta5=0;%波束形成目标角度5
theta=-90:0.1:90;%刻画角度坐标
f=(0:N-1)fs/N;%刻画频率坐标
for i=1:M
tao(i)=round((i-1)dsind(theta0)/cfs);%各阵元在目标方向时延差取整，在整个信号上截取时延后的信号
s(i,:)=s0(1,tao(i)+1:tao(i)+N);%生成各阵元时域信号
sf(i,:)=fft(s(i,:),N);%转换到频域
end
for i=1:M
tao1(i)=round((i-1)dsind(theta1)/cfs);%各阵元在目标方向时延差取整，在整个信号上截取时延后的信号
ss(i,:)=s0(1,tao1(i)+1:tao1(i)+N);%生成各阵元时域信号
sfs(i,:)=fft(ss(i,:),N);%转换到频域
end
for i=1:M
tao1(i)=round((i-1)dsind(theta2)/cfs);%各阵元在目标方向时延差取整，在整个信号上截取时延后的信号
ss1(i,:)=s0(1,tao1(i)+1:tao1(i)+N);%生成各阵元时域信号
sfs1(i,:)=fft(ss1(i,:),N);%转换到频域
end
for i=1:M
tao1(i)=round((i-1)dsind(theta3)/cfs);%各阵元在目标方向时延差取整，在整个信号上截取时延后的信号
ss2(i,:)=s0(1,tao1(i)+1:tao1(i)+N);%生成各阵元时域信号
sfs2(i,:)=fft(ss2(i,:),N);%转换到频域
end
for i=1:M
tao1(i)=round((i-1)dsind(theta4)/cfs);%各阵元在目标方向时延差取整，在整个信号上截取时延后的信号
ss3(i,:)=s0(1,tao1(i)+1:tao1(i)+N);%生成各阵元时域信号
sfs3(i,:)=fft(ss3(i,:),N);%转换到频域
end
for i=1:M
tao1(i)=round((i-1)dsind(theta5)/cfs);%各阵元在目标方向时延差取整，在整个信号上截取时延后的信号
ss4(i,:)=s0(1,tao1(i)+1:tao1(i)+N);%生成各阵元时域信号
sfs4(i,:)=fft(ss4(i,:),N);%转换到频域
end
sf =sf+sfs+sfs1+sfs2+sfs3+sfs4;
figure(1)
plot(f,real(sf));legend(‘阵元1’,‘阵元2’,‘阵元3’,‘阵元4’,‘阵元5’);

fl=9e6;fh=11e6;%筛选出信号所在频率谱线进行波束形成
n1=round(fl/fsN);n2=round(fh/fsN);
for i=1:length(theta)
for j=n1:n2
for k=1:M
n11(i)=(k-1)dsind(theta(i))/cfs/N;%对任意角度在该阵元上的时延转换为时域上平移的点
sf1(k)=sf(k,j)exp(-1j2pijn11(i));%时域时延，频域相移
end
s1(j)=sum(sf1);%该谱线上各阵元求和
end
s2(i)=sum(abs(s1(j))^2);%该角度上各谱线求和
end
figure(2)
subplot(121)
plot(theta,s2/max(s2));title(‘频域波束形成波束图’);xlabel(‘角度（°）’);ylabel(‘幅值’);
subplot(122)
polarplot(theta*pi/180,s2/max(s2));title(‘频域波束形成极坐标图’);

波束形成结果极坐标图