第四章 均匀线阵的波束形成

各位同学，今天我们来聊聊均匀线阵的波束形成。说实话，这是整个数字波束形成里最基础、也最经典的内容。我当年刚入行时，第一个实战项目就是搞这个，踩了不少坑，今天一并分享给你们。

4.1 均匀线阵的阵列流型

先说说什么是阵列流型。你想想看，一个均匀线阵，就是N个天线等间距排成一条直线。每个天线接收到的信号，因为到达时间不同，会有一个相位差。把这些相位差组合起来，就形成了阵列流型向量。

数学上，对于一个间距为d的N元均匀线阵，假设信号从θ方向入射，阵列流型向量可以写成：

a(θ) = [1, e^(-j·2π·d·sinθ/λ), e^(-j·4π·d·sinθ/λ), ..., e^(-j·2π·(N-1)·d·sinθ/λ)]^T

这里λ是波长。嗯，要注意一点——d通常取半波长，也就是λ/2。为什么？因为这样可以避免栅瓣。我在项目中遇到过有人为了省成本把间距拉大，结果波束图上一堆假目标，排查了三天才发现是阵列设计的问题。

4.2 常规波束形成（CBF）算法

常规波束形成，说白了就是给每个阵元的信号加一个权重，然后求和。这个权重就是阵列流型的共轭。为什么要共轭？因为要补偿相位差，让所有阵元的信号同相相加。

公式很简单：

y(t) = w^H · x(t)

其中w是权重向量，x(t)是接收信号向量。对于常规波束形成，w = a(θ₀)，θ₀是期望的指向方向。

我刚开始做这个的时候，总觉得这算法太简单了，不就是加权求和嘛。但后来发现，简单归简单，里面门道不少。比如权重的归一化问题，如果不做归一化，输出功率会随着阵元数变化，影响后续的检测门限设置。

4.3 波束扫描与指向

波束扫描，就是让主瓣指向不同的方向。怎么做？改变权重向量里的θ₀就行了。你想想看，如果我想让波束指向30°，就用a(30°)作为权重；想指向60°，就用a(60°)。

实际工程中，我们通常预先计算好一组权重，存起来。扫描时直接查表，省得实时计算。我记得有一次做机载雷达项目，要求波束在±60°范围内扫描，步进0.5°。算下来要存241组权重，每组权重N个复数。当时DSP的存储空间有限，我不得不做了一些优化——只存一半的权重，利用对称性推导另一半。

波束指向的精度取决于两个因素：

• 阵元数N：N越大，波束越窄，指向精度越高
• 量化位数：权重系数的量化误差会影响指向精度

4.4 波束图分析

波束图是衡量波束形成性能的重要工具。它描述了波束输出功率随入射角的变化关系。对于均匀线阵，波束图有以下几个特征：

参数	公式	说明
主瓣宽度	BW ≈ 0.886λ/(N·d·cosθ₀)	阵元数越多，波束越窄
第一旁瓣电平	约-13.5dB	均匀加权时的固定值
栅瓣位置	sinθ_g = sinθ₀ ± k·λ/d	d>λ/2时可能出现

这里我要特别强调一下旁瓣电平。均匀线阵的-13.5dB旁瓣，在很多场合是不够的。比如机载雷达，地面杂波从旁瓣进来，可能比主瓣目标信号还强。这时候就需要加窗函数来压低旁瓣，但代价是主瓣会变宽。

嗯，说到加窗，我建议你们先掌握均匀加权，理解它的特性，再去研究切比雪夫、泰勒等加权方法。基础打牢了，后面才能灵活运用。

4.5 代码实现示例

最后，给一个简单的MATLAB代码示例，演示均匀线阵的波束扫描：

% 参数设置
N = 16;          % 阵元数
d = 0.5;         % 阵元间距（波长归一化）
theta_scan = -90:0.5:90;  % 扫描角度范围

% 计算阵列流型
theta_rad = theta_scan * pi / 180;
a = @(theta) exp(-1j * 2 * pi * d * (0:N-1)' * sin(theta));

% 波束扫描
beam_pattern = zeros(size(theta_scan));
for idx = 1:length(theta_scan)
    w = a(theta_rad(idx)) / N;  % 归一化权重
    beam_pattern(idx) = abs(w' * a(theta_rad(idx)))^2;
end

% 绘制波束图
plot(theta_scan, 20*log10(beam_pattern));
xlabel('角度(度)'); ylabel('归一化功率(dB)');
grid on;

这段代码很简单，但包含了波束形成的核心思想。你可以试着改改N和d，看看波束图怎么变化。动手试试，比光看理论强多了。