3.3 阵列天线、波束合成与MIMO物理

3.3 阵列天线、波束合成与MIMO物理

单个天线的方向图是天生的，阵列的方向图是可编程的。本节把 N 个天线单元排成阵，用干涉原理证明一个关键事实：方向图从"阵列因子"与"单元因子"的乘积中生长出来，而阵列因子完全由各单元的幅相配置决定——改相位就是改波束指向，改幅度就是改旁瓣轮廓。这既是相控阵雷达毫秒级扫描的原理，也是 5G 大规模 MIMO 把波束当资源调度的底牌（旧版 MIMO 物理基础一节并入本节后半）。学完本节，你将能亲手算出阵因子、解释栅瓣从哪来，并理解"空间自由度"这个 6G 时代的高频词在物理上指什么。下一节再看口径与宽带家族如何把方向性推向极致。

一、阵因子：干涉原理的作业本

均匀直线阵（N 元、间距 d）的远场是各单元贡献的相干叠加。若单元全同且等幅激励，总场等于单元因子乘上阵列因子：

AF(θ) = Σ exp(j(n−1)·(k·d·cosθ + φ))

其中 φ 是相邻单元的相位步进。当 k·d·cosθ + φ = 0 时各单元贡献同相相加，出现主瓣——波束指向 θ₀ 由 φ 唯一决定：cosθ₀ = −φ/(k·d)。改一个相位参数就指挥整个波束，这正是"电扫"二字的物理含义；机械转动从此可以退休，扫描速度从秒级跃到微秒级。

方向性随阵长的代价关系同样干脆：N 元等幅阵的主瓣宽度约 2/(N)·波长/间距 弧度量级——阵列孔径越大波束越窄，方向性按阵元数成倍增长。天下没有免费的午餐：阵元数、带宽、功耗在馈电网络里互相较劲。

# 均匀直线阵阵因子与相位扫描
import numpy as np

N, d_over_lam = 8, 0.5            # 八元半波间距
theta = np.radians(np.linspace(-90, 90, 1801))
k_d = 2*np.pi*d_over_lam

def af(phi):
    psi = k_d*np.cos(theta) + phi
    num = np.sin(N*psi/2)
    den = np.sin(psi/2)
    af_val = np.where(np.abs(den) < 1e-9, N, num/np.where(np.abs(den)<1e-9,1,den))
    return np.abs(af_val)/N

for steer in [0, 30]:
    phi = -k_d*np.cos(np.radians(steer))
    pattern_db = 20*np.log10(np.maximum(af(phi), 1e-6))
    peak = theta[np.argmax(af(phi))]
    first_null = theta[np.where(af(phi) < 0.7071)[0][0]]
    print(f"指向 {steer:3d} 度：峰值 {np.degrees(peak):6.1f} 度，"
          f"零点宽度约 {2*abs(np.degrees(first_null)-steer):.1f} 度")
# 修改 d_over_lam 为 1.0 重跑：栅瓣现身——同相叠加方向不止一处
​

把间距改到一个波长再跑：除目标指向外，90 度方向也冒出等强波束——栅瓣。它是空间采样定理的直接后果：单元间距超过半波长，"采样"过疏，波束混叠。相控阵雷达死守半波间距就是为封死栅瓣；而稀疏阵反过来主动放宽间距、用优化算法压制旁瓣，用阵元数换成本。

相位步进指挥波束转向

二、从相控阵到 MIMO：空间自由度的物理账

相控阵用一个波形、靠相位把能量拼成一个窄波束；MIMO 则更进一步——各天线发不同的信号，把空间本身开辟成一条新信道。其物理底牌是散射环境：城市里建筑物把信号打散成多条独立路径，收发两侧的信道矩阵 H 的秩（独立并行子信道数）取决于两端天线数与散射丰富度，容量按 min(Nt, Nr)·log₂(1+SNR) 增长。这些"空间子信道"不是天线造出来的，而是散射环境提供的——天线阵列只是把环境赠送的自由度收进钱包。空旷环境（信道强相关）里 MIMO 容量缩水，正是钱包里没现金的窘境。

互耦是阵列绕不开的暗账：单元挨得近，彼此的近场互相改变阻抗与方向图（3.1 近场概念的阵列版），设计时要用耦合矩阵校准。大规模 MIMO 还有一笔数字账：64 通道的射频链路成本远超 4 通道，于是混合波束赋形用少数射频链路加大量移相器折中——物理原理不变，工程架构在成本与性能间找平衡。

发车平台的调度手册到此凑齐三件套：单元性能（3.1、3.2）、编队控制（阵因子）、空间复用（MIMO）。第8章的相控阵雷达与 5G 基站，全是这三件套的系统级出口。

三、阵因子的数字实验：间距、扫描与加权

阵因子是一条干涉公式，与其背结论，不如跑三个对照实验把规律焊进脑子。八元均匀直线阵，扫描到三十度：半波间距时主瓣之外再无三dB以上的瓣；把间距拉到一个波长，三十度对称的位置冒出一个几乎同高的瓣——栅瓣现身，正是空间采样欠密的混叠。再加一组幅度对照：等幅激励的第一旁瓣约负十三dB；换海明加权，旁瓣直接压到负四十dB 量级，代价是主瓣胖了约七成。

# 八元直线阵方向图：间距、扫描与幅度加权的对照实验
import numpy as np

th = np.radians(np.linspace(-90, 90, 1801))
def af_db(n, d, steer_deg, w=None):
    psi = 2*np.pi*d*(np.sin(th) - np.sin(np.radians(steer_deg)))
    w = np.ones(n) if w is None else w
    af = np.abs(np.exp(1j*np.outer(np.arange(n), psi)).T @ w)
    return 20*np.log10(af/af.max() + 1e-12)

for d, tag in [(0.5, "半波间距"), (1.0, "全波间距")]:
    p = af_db(8, d, 30.0)
    print(f"{tag} 扫描30度：-3dB以上角度样本 {p[p > -3].size} 个（含主瓣）")

sl_zone = np.abs(th) > np.radians(16)          # 避开主瓣的区域
p_uni = af_db(8, 0.5, 0.0)
p_ham = af_db(8, 0.5, 0.0, np.hamming(8))
print(f"等幅激励最高旁瓣 {p_uni[sl_zone].max():.1f} dB")
print(f"海明加权最高旁瓣 {p_ham[sl_zone].max():.1f} dB，主瓣展宽是代价")
​

三个实验对应三条工程军规：间距守半波长红线（栅瓣是混叠，事后无法滤除）；想要低旁瓣就把幅度加权写进馈电网络（面阵同理，行乘列两次加权）；扫描角越大有效口径越短、波束越胖，大角度扫描的性能预算要按最胖时留。这套数字实验也是 3.2 度量衡的自检工具：指标表上的波束宽度、旁瓣电平，都应当能被这条公式粗算复现。

本节要点回顾

• 阵列因子：方向图 = 单元因子 × 阵因子，阵因子由间距与幅相配置决定，波束指向由相位步进唯一指定。
• 电扫原理：改相位不改结构即可微秒级扫描，机械转动在相控阵里成为历史。
• 栅瓣条件：间距超过半波长即空间采样欠密，波束混叠出栅瓣；稀疏阵以算法换阵元成本。
• 旁瓣与主瓣的交易：幅度加权压旁瓣必展宽主瓣，方向性与干扰免疫在同一个旋钮上互相牵制。
• MIMO 底牌：散射环境的独立路径数决定信道矩阵秩，容量随最小天线数对数增长；互耦与信道相关是两笔暗账。

下一节看发车平台的两个进阶家族：把口径做大到米级的反射面天线，与把带宽做到倍频程的非频变天线。