均匀线阵方向图.m

来自「matlab,线阵」· M 代码 · 共 44 行

% 绘制均匀线阵方向图 
% 《空间谱估计理论与算法》 P31 
clc
clear

M = 16;             % 阵元数
k = 0.5;             % k = d/lamada = 1/2

% 形成波束指向doa方向的权值
doa = 30;    % doa=0，为静态方向图
m = [0 : M-1];
w = exp(-j*2*pi*k*m'*sin(doa*pi/180)); 

% 波束扫描，绘制方向图
theta = linspace(-90,90,360);                        % 方位角
a = exp(-j*2*pi*k*m'*sin(theta*pi/180));             % 导向适量
Y = abs(w'*a);
Y = 20*log10(Y/max(Y));
subplot(211)
plot(theta, Y); title('未加权方向图');
axis([-90,90,-50,0]);
grid on; xlabel('方位角/(\circ)');ylabel( 'P/(dB)');

% =================================
% 加切比雪夫权后方向图

M = 16;             % 阵元数
k = 0.5             % k = d/lamada = 1/2

% 形成波束指向doa方向的权值
doa = 30;    % doa=0，为静态方向图
m = [0 : M-1];
w = exp(-j*2*pi*k*m'*sin(doa*pi/180)); 
w = chebwin(M, 30) .* w;    % 加切比雪夫权

% 波束扫描，绘制方向图
theta = linspace(-90,90,360);                        % 方位角
a = exp(-j*2*pi*k*m'*sin(theta*pi/180));             % 导向适量
Y = abs(w'*a);
Y = 20*log10(Y/max(Y));
subplot(212)
plot(theta, Y); title('加-30dB切比雪夫权方向图');
axis([-90,90,-50,0]);
grid on; xlabel('方位角/(\circ)');ylabel( 'P/(dB)');