MATLAB实战：手把手教你仿真三种天线阵列的波束形成（附完整代码）

MATLAB实战：三种天线阵列波束形成仿真全解析

天线阵列技术在现代无线通信、雷达探测等领域扮演着核心角色。不同于单天线系统，阵列天线通过精确控制多个阵元的信号相位和幅度，能够实现波束的灵活指向和干扰抑制。本文将带您深入理解均匀直线阵列(ULA)、均匀矩形阵列(URA)和均匀圆形阵列(UCA)的波束形成原理，并通过MATLAB代码实现完整的仿真流程。

1. 阵列天线基础与仿真环境搭建

阵列天线的核心思想是利用多个空间分布的阵元协同工作，通过调整各阵元的加权系数，在特定方向形成增益较高的波束。这种技术被称为数字波束形成(DBF)，它使我们能够在不物理移动天线的情况下，实现波束的电子扫描。

仿真环境配置是第一步。我们需要准备MATLAB R2020a或更高版本，并确保安装了以下工具箱：

• Signal Processing Toolbox
• Phased Array System Toolbox（可选，本文采用基础实现）

% 基本参数设置 c = 3e8; % 光速(m/s) f0 = 15e6; % 载波频率15MHz lambda = c/f0; % 波长计算 fs = 4*f0; % 采样频率 N = 1024; % 采样点数 t = (0:N-1)'/fs; % 时间序列

远场假设是阵列分析的基础，它要求目标距离远大于阵列孔径。我们设定信号源位于：

TxLoc = [1000e3, 1000e3, sqrt(2)*1000e3]; % 方位45°,仰角45°

2. 均匀直线阵列(ULA)实现与特性分析

ULA是最简单的阵列形式，所有阵元沿直线等间距排列。其波束形成能力主要体现为方位角分辨，而仰角分辨能力有限。

阵列配置参数：

d = 8; % 阵元间距8米(0.4λ) M = 30; % 阵元数量 RxLoc = [d*(0:M-1)', zeros(M,2)]; % 沿x轴排列

ULA的阵列流形向量计算是关键：

k = 2*pi/lambda; % 波数 phi = linspace(0,2*pi,360); % 方位角扫描范围 theta = linspace(0,pi/2,90); % 仰角扫描范围 a_ULA = exp(1j*k*d*(0:M-1)'*sin(theta0)*cos(phi0)); % 特定方向的流形向量

仿真结果显示几个重要特性：

提示：ULA在实际应用中通常需要与其他传感器配合，才能实现完整的三维空间覆盖。

3. 均匀矩形阵列(URA)设计与波束优化

URA通过在两个正交维度上排列阵元，实现了方位和仰角的二维分辨能力。相比ULA，URA的波束控制更加灵活。

阵列布局代码：

dx = 8; dy = 8; % x,y方向间距 Mx = 30; My = 30; % x,y方向阵元数 RxLoc = zeros(Mx*My,3); for i=1:Mx for j=1:My RxLoc((i-1)*My+j,:) = [(i-1)*dx,(j-1)*dy,0]; end end

URA的阵列流形计算采用Kronecker积：

ax = exp(1j*k*(0:Mx-1)'*dx*sin(theta0)*cos(phi0)); ay = exp(1j*k*(0:My-1)'*dy*sin(theta0)*sin(phi0)); a_URA = kron(ay,ax); % 二维阵列流形

URA的波束特性明显优于ULA：

1. 二维分辨能力：可同时精确控制方位和仰角
2. 无测角模糊：360°范围内唯一确定波束指向
3. 旁瓣控制：通过二维窗函数优化波束形状

% 加窗处理示例 wx = hamming(Mx); wy = hamming(My); W = wy*wx'; % 二维窗函数 a_URA_weighted = a_URA.*W(:); % 加权后的流形

4. 均匀圆形阵列(UCA)实现与全向扫描

UCA因其对称结构，在360°范围内具有一致的波束性能，特别适合全向扫描应用。

圆形阵列生成：

M = 96; % 阵元数量 R = 120; % 阵列半径 alpha = (0:M-1)'*2*pi/M; RxLoc = [R*cos(alpha), R*sin(alpha), zeros(M,1)];

UCA的流形计算需考虑每个阵元的独特位置：

xm = R*cos(alpha); ym = R*sin(alpha); a_UCA = exp(1j*k*(xm*cos(phi0)*sin(theta0) + ym*sin(phi0)*sin(theta0)));

UCA的独特优势包括：

• 全向一致性：各方向波束性能均匀
• 灵活配置：可通过选择激活不同扇区实现重构
• 多模式工作：同时形成多个独立波束

% 多波束形成示例 phi_targets = [30, 150, 270]; % 三个目标方向 theta_target = 45; % 固定仰角 beams = zeros(M,length(phi_targets)); for i=1:length(phi_targets) beams(:,i) = exp(1j*k*(xm*cosd(phi_targets(i))*sind(theta_target) + ... ym*sind(phi_targets(i))*sind(theta_target))); end

5. 三种阵列性能对比与工程选择指南

在实际系统设计中，阵列类型的选择需综合考虑多方面因素。我们通过系统仿真，量化比较三种阵列的关键指标。

性能对比表：

工程选择建议：

1. 移动通信基站：URA适合蜂窝扇区覆盖
2. 雷达预警系统：UCA提供全方位监视能力
3. 车载毫米波雷达：ULA满足前方区域探测需求

% 性能评估代码示例 angles = 0:359; dbi_ULA = zeros(size(angles)); dbi_URA = zeros(size(angles)); dbi_UCA = zeros(size(angles)); for i=1:length(angles) % 各阵列在指定方向的响应计算 % ...（具体计算代码） end plot(angles,dbi_ULA,angles,dbi_URA,angles,dbi_UCA); legend('ULA','URA','UCA'); grid on;

在5G毫米波通信系统的实测中，URA阵列在28GHz频段实现了±60°的扫描范围，波束宽度约5°，完全满足高速移动场景的需求。而UCA在机场雷达系统中，则展现了同时跟踪200个目标的能力。