常规波束成形算法仿真.doc

有一8个阵元的均匀线阵，假设一功率为1的单频信号，假设信号从10度方向入射到基阵。不考虑噪声和其他干扰信号。则其波束图和方位图的仿真为：
%1.beampattern
clear all
N=8;
theta=10*pi/180;
As=exp(-j*pi*[0:N-1]'*sin(theta));
wc=As/N;
thetas=-90:90;
tm=thetas*pi/180;
Am=exp(-j*pi*[0:N-1]'*sin(tm));
p=wc'*Am;
P=20*log(abs(p));
plot(thetas,P);

%2．DOA
clear all
N=8;
d=1/2; %相对于波长的阵元间距
theta=10*pi/180; %信号从10度方向入射到基阵
As=exp(-j*pi*[0:N-1]'*sin(theta));%导向向量a(10)
thetas=-90:90;
tm=thetas*pi/180;
Rx=As*As'; %数据协方差矩阵
for i=1:length(tm)
Am=exp(-j*pi*[0:N-1]'*sin(tm(i)));
w=Am/N; %常规波数形成加权向量
b=w'*Rx*w;
B=10*log(abs(b));
C(i)=B;
i=i+1;
end
plot(thetas,C)