有一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)