😆
理论知识参考:数字多波束相控阵基础知识
仿真代码
close all; clc;
fc = 1e+6; %载波频率
fs = 128*fc; %采样频率
c = 3e+8; %光速
Wave_length = c/fc; %波长
figure;
for var = 0.1:0.1:1
N = 10; %阵元个数
n = (0:N-1)';%阵元编号
d = var*Wave_length; %阵元间距
theta_0 = 10; %期望角度
tau = d*sind(theta_0)/c; %相邻阵元延时
W = exp(1i*2*pi*fc*n*tau); %权值矩阵
step = 1;
theta = -90:step:90;
for m = 1:length(theta)
V = exp(1i*2*pi*fc*n*d*sind(theta(m))/c); %方向矢量
Beam(m) = abs(W'*V); %方向图
end
plot(theta,Beam);
title(sprintf('d=%2f',var));
pause(0.5);
end
运行结果(改变阵元间距d,阵元个数N=10)
运行结果(改变阵元个数N,阵元间距d=0.5)