线性阵列doa角度估计，适应于均匀线性阵列的情况matlab代码.zip

clear all; close all; element_num =80; %阵元数 d_lamda = 1/2;%阵元间距与波长lamda的关系 theta = -90:0.5:90; theta_in = 40; %来波方向； a_signal =exp(1j*2*pi*d_lamda*sin(theta_in*pi/180)*[0:element_num-1]'); %信号的导向矢量 for m=1:length(theta) a=exp(1j*2*pi*d_lamda*sin(theta(m)*pi/180)*[0:element_num-1]'); singal_out(m)=a_signal'*a; end%接收的导向矢量，从-pi/2到pi/2 patternmag=abs(singal_out); patternmagnorm = patternmag/max(max(patternmag));%归一化 patternmagdbBnorm=20*log10(patternmagnorm);%归一化分贝 figure(1) plot(theta,patternmagdbBnorm); grid on; xlabel('theta/radian'); ylabel('amplitude/dB'); title('阵元均匀线阵方向图');