clc;
clear all;
close all;
imag=sqrt(-1);
element_num=8;%阵元数为8
d_lamda=1/2;%阵元间距d与波长lamda的关系
theta=linspace(-180,180,200);
theta0=60;%来波方向
w=exp(imag*2*pi*d_lamda*sin(theta0*pi/180)*[0:element_num-1]');
for j=1:length(theta) %(入射角度,即来波方向,计算阵列流形矩阵A)
a=exp(imag*2*pi*d_lamda*sin(theta(j)*pi/180)*[0:element_num-1]');
p(j)=w'*a;% w'为w的共轭转置.方向图的比值约去了信号值
end %对整个空间进行波束扫描就可确定来波方向。
figure;
plot(theta,abs(p)),grid on
xlabel('theta')
ylabel('amplitude')
title('8阵元均匀线阵方向图')