Mahafza那本《雷达系统分析与设计MATLAB版》第3版我读的时候有个很深的印象前面几章讲雷达方程、检测、波形设计、跟踪内容非常扎实但读起来更像教材到了第五部分“关于雷达的一些特殊主题”画风突然变得务实起来。尤其是第15章《相控阵天线》按目录位置它是靠后的“特殊主题”可仔细读下来就会发现它解决的问题——如何让雷达快速改变波束指向而不机械转动天线——几乎成了现代雷达尤其是毫米波雷达、4D成像雷达等新方向共同的核心底座。这篇笔记就是基于第15章的一次沉淀。我会把这一章里的主要脉络梳理出来把阵因子、栅瓣约束、波束宽度、量化误差这些比较硬核的概念讲得更直白一些同时补上一些我在MATLAB里实际跑方向图仿真时踩过的坑和收获。如果你正在读这本书或者在做阵列信号处理相关的仿真再或者想知道毫米波雷达阵列背后到底是什么原理这篇内容应该能帮你省一点时间。1. 一个挂在“特殊主题”里的章节怎么成了现代雷达的必修课1.1 第15章在整本书里的位置说明了什么这本书的章节安排很清晰基础概念、雷达距离方程、波形设计、检测、MTI与脉冲多普勒、跟踪雷达然后进入特殊主题部分。第一次看书目录的人很容易跳过后面觉得“特殊主题”意味着小众。但第15章相控阵天线这个位置恰恰反映了雷达技术演进史里一个有意思的事实。在相控阵进入工程实用之前雷达要改变波束方向主要靠机械转动天线。这种体制的响应速度慢机械结构复杂而且很难同时跟踪多个目标。相控阵的思路是把单个大天线拆成成百上千个小阵元每个阵元后面接一个可控的移相器更先进的接T/R组件通过控制各阵元之间的相位关系让波束在电学意义上瞬间偏转。这个转变从机械扫描走向电扫描是把相控阵放在“特殊主题”章节的历史背景。Mahafza把这一章放在第五部分还有一个原因相控阵本身牵扯到的数学工具并不算深但它把这些工具综合在了一起。坐标变换、方向图乘积定理、复数运算、信号模型哪一个单拎出来都不难组合到一起就需要一定的熟练度。读这一章之前如果能把书里前几章的窄带信号模型内容翻回去看看会顺畅很多。1.2 读这一章之前需要补足的基础知识清单如果你和我一样是带着一点基础开始读这一章的我建议你先在心里过一遍以下三个概念。第一个是窄带信号的复数表示。雷达信号经过混频、采样之后通常被表示为复数基带形式阵列处理几乎全部工作在这些复数据上。相控阵里的“相位”概念本质上就是复数相乘旋转一个角度。第二个是方向图乘积定理。一个阵列的远场方向图等于阵元方向图与阵因子的乘积。阵元方向图描述单个天线单元本身的辐射特性阵因子描述阵列布局造成的干涉效果。阵因子是这一章的绝对主角很多公式的推导都围绕它展开。第三个是空间频率的概念。波束从一个方向来阵元之间存在波程差这个波程差对应相位差而相位差随角度的变化率就是空间频率。理解了空间频率再去看波束指向、栅瓣、空间混叠这些概念就会觉得它们和时域信号处理里的频率混叠几乎是一回事。这些基础不牢后面看到阵因子公式里那些sinθ和cosθ容易一头雾水基础到位了第15章的推导主线是很清晰的。2. 从波程差到阵因子相控阵“扫描”的几何本质2.1 相邻阵元的相位差是怎么一步步推出的假设一个N元均匀线阵阵元间距为d信号以角度θ入射。由于存在波程差相邻两个阵元接收到信号的时间不同体现在相位上就是Δφ 2πd·sinθ / λ这个式子只有几个变量但含义很深。d和λ确定了阵列的几何尺度sinθ确定了信号方向到阵列法线方向的投影关系。之所以用sinθ而不是θ是因为阵列对方向角的响应不是线性的用sinθ做自变量在很多推导里会大幅简化方向图函数也会保持周期结构。如果对所有阵元的接收信号不施加任何额外相位那么阵列方向图会在法线方向θ0取得最大值因为此时各路信号波程差为零同相叠加。要让波束偏转到某个方向θ₀就需要在各阵元之间人为引入一组相位梯度去抵消来自θ₀方向信号的波程差。也就是说第n个阵元需要补偿的相位是αₙ −n·2πd·sinθ₀ / λ这里n从0到N−1。这组补偿相位就是“相控”二字的来源。它和常规天线不一样的地方在于天线本身没有动是相位梯度在数学上完成了对波程差的“预补偿”让期望方向上的各路信号在合并时重新对齐。2.2 从单个相位补偿到完整阵因子的推导把接收信号按阵元加权求和再取模值就得到阵因子。对均匀加权所有阵元幅度权重都为1的情况阵因子可以写成等比级数求和AF(θ) Σₙ exp[j·n·(2πd/λ)(sinθ − sinθ₀)]这个等比数列求和的解析形式并不复杂AF(θ) sin[N·(πd/λ)(sinθ − sinθ₀)] / sin[(πd/λ)(sinθ − sinθ₀)]方向图则是|AF(θ)|的归一化对数形式通常以dB为单位绘制。如果你第一次在MATLAB里用这个公式画方向图会发现一个非常直观的现象主瓣出现在sinθ sinθ₀的位置旁瓣周期性排列副瓣电平约为−13.26dB均匀加权时第一副瓣。我读这一章时花了挺长时间才真正理解阵因子公式里的相位差分结构。它看起来复杂但本质就是一句话N个等幅等相位梯度信号做相干叠加最大值出现在相位梯度恰好抵消波程差的方向。所有后续的栅瓣分析、波束宽度计算、加窗设计都是在这个公式基础上展开的。2.3 为什么很多人会误以为相控阵真的在“转动”天线我看到不少学习材料里有个说法相控阵通过改变相位等效于“虚拟转动”天线。这个说法方便理解但容易造成概念混淆。孔径投影效应并不是转动阵元的相对位置没有变阵列的物理孔径在任意方向上的投影随扫描角变化了。这个投影效应直接导致两个后果一是波束宽度随扫描角增大而展宽二是阵列有效接收面积随扫描角增大而下降对应增益下降。如果脑子里只有“虚拟转动”的概念这两个后果就很难解释清楚。更准确的理解方式是相控阵的相位梯度改变的是空间中各点信号的干涉条件。哪些方向上各路信号同相叠加哪里就是主瓣那些方向上出现周期性同相叠加哪里就是栅瓣。干涉条件是几何和波长的共同结果和天线是否“转动”没有任何关系。3. 栅瓣的数学边界d/λ约束与60度扫描的完整推导3.1 栅瓣到底什么时候出现从周期性说起栅瓣是相控阵里最需要警惕的现象。简单说除了期望的扫描方向θ₀在其他某些方向上阵元的相位补偿恰好也让各路信号同相叠加于是形成和主瓣幅度接近的额外波瓣。这就是栅瓣。从数学上看阵因子公式里出现的项是(πd/λ)(sinθ − sinθ₀)。等比级数求和的分子sin[Nx]/sin[x]在xkπ时取极值。问题在于当xπk1时对应的方向也可能落在真实空间里即|sinθ|≤1这时栅瓣就出现了。所以避免栅瓣的核心约束是|(πd/λ)(sinθ − sinθ₀)| π对所有满足|sinθ|≤1的θ成立。化简之后就是典型的d/λ约束d/λ 1 / (1 |sinθ₀|)这个公式在几乎每本相控阵书里都有但很多人直接背结论不知道它怎么来的也不知道不同扫描范围对应什么数值。3.2 扫描60度时为什么d要小于0.536λ一个完整计算过程把上面公式代数字进去当最大扫描角θ₀60°时sin60° ≈ 0.8661 0.866 ≈ 1.866d/λ 1/1.866 ≈ 0.536也就是说如果要求波束在±60°范围内扫描阵元间距不能超过0.536λ。工程上取整通常会预留一定余量采用d0.5λ这也是绝大多数均匀线阵仿真的标准配置。d0.5λ还有一个额外好处即使扫描到接近±90°也不会有栅瓣。此时1|sinθ₀|最大值接近2d/λ0.5严格小于1/2刚好卡在安全边界内。这就是为什么大量教材和仿真示例默认用半波长间距不完全是出于方便背后是有数学约束支撑的。如果只工作在很窄的扫描范围比如±20°那么d/λ可以放宽到1/(10.342)≈0.745λ。更大的阵元间距意味着同样的孔径可以用更少的阵元成本更低所以实际工程里不会盲目选择半波长间距而是根据扫描范围反推最优间距。3.3 工程上对栅瓣的额外限制栅瓣可以消失但副瓣还有指标这里必须多说一句。栅瓣约束保证的是“空间中不出现第二个和主瓣等高的极大值”但副瓣电平是另一个维度的问题。均匀加权时第一副瓣是−13.26dB在某些应用场景比如侦察、气象雷达这个副瓣水平是不够的因为副瓣方向进来的强干扰或强杂波会直接污染主瓣信号。要压低副瓣可以在阵元级别做幅度加权也就是给每个阵元乘一个窗函数比如汉明窗、泰勒窗、切比雪夫窗。后面会详细说这里先记住一个结论栅瓣是阵元间距决定的副瓣是幅度加权决定的两者要分开控制。另外栅瓣有一个非常隐蔽的表现形式。有时候主瓣在期望方向但在阵因子周期结构的某个栅瓣位置正好落在低副瓣窗函数的零点上看起来没问题可一旦扫描角变化栅瓣和窗函数零点的相对位置也会变化问题就可能突然冒出来。设计阵列时不能只看静态方向图要看整个扫描范围内的动态方向图变化。4. 波束宽度与扫描损失方向图里的能量守恒和窄带假设4.1 波束宽度公式的由来从sin(x)/x方向图到0.886λ公式继续使用均匀加权均匀线阵来讨论。阵因子的主瓣形状近似为sin(Nx)/x的形式也就是sinc函数型。半功率点对应x满足N·(πd/λ)(sinθ − sinθ₀) ≈ ±0.443π这其实是sinc函数sin(πu)/(πu)在3dB处的u值约为0.443的结论。从这个近似出发可以得到均匀线阵波束宽度以弧度计Δθ₃dB ≈ 0.886λ / (N·d·cosθ₀)在法线方向cos01这个公式退化为常见的Δθ₃dB≈0.886λ/(Nd)。Nd就是阵列总孔径长度所以这个公式的直觉版本是孔径越大、波长越短波束越窄。这个公式在MATLAB里验证非常方便。设N32, d0.5λ法线方向的3dB波束宽度约为0.886/(32×0.5)×(180/π)≈3.17°。我习惯先用理论公式估算再用方向图代码数值测量两者吻合度通常在0.1°以内。4.2 扫描角增大时波束为什么展宽以及阵元增益的cosθ项注意公式里cosθ₀在分母上。当扫描角从0°增大到60°cosθ₀从1降到0.5波束宽度直接变成法线方向的两倍。也就是说扫描到边缘时角度分辨力显著变差。这里面有一个看起来矛盾但其实自洽的物理图景。阵列物理孔径是固定的但在扫描方向上的有效投影孔径是Nd·cosθ₀投影越小等效孔径越小波束自然越宽。同时每个阵元的有效接收面积在偏离法线方向时也按cosθ₀下降两者叠加导致大扫描角时的阵列增益损失接近cos²θ₀的量级。当然这是理想模型。实际工程中阵元方向图本身可能不是各向同性的扫描到一定角度后增益下降更严重。所以很多大型相控阵用“三坐标雷达”的设计思路把扫描限制在两个维度以内或者干脆用多个面阵覆盖不同扇区。4.3 相干积累与时延波束形成窄带假设何时失效这一章里还有一个大概念值得专门提出来相控阵用相位补偿来模拟波程差补偿但这个近似仅在窄带条件下成立。为什么相位补偿本质上是对时延的一阶近似。一个信号从θ方向到达阵列相邻阵元之间的真实时延是Δτ d·sinθ/c。在窄带情形下信号包络变化很慢用相位旋转近似这个时延误差可以忽略。但信号带宽变大后包络沿阵元间的延迟就变得不可忽略直接用相移会导致波束指向偏差和波形失真。量化判断标准是孔径渡越时间bandwidth积B·(D·sinθ/c)是否远小于1。其中B是信号带宽D是阵列孔径。比如一个孔径10米的阵列扫描到30°c约为3×10⁸m/s孔径渡越时间约为16.7ns。如果要保持瞬时带宽100MHz乘积是1.67这早就不能满足窄带假设了。宽带相控阵因此需要用到真时延单元而不是简单的移相器。这一章在教材里也许只是一带而过但在实际工程尤其是宽带雷达和通信一体化系统里是一个核心痛点。5. 副瓣、窗函数与量化误差我在这章里最想分享的MATLAB细节5.1 均匀加权方向图的仿真从阵因子直接画方向图《雷达系统分析与设计》第15章后半部分的重点是仿真实例我自己在MATLAB里重写了一个最基础的均匀线阵方向图函数。核心逻辑只有几行N 32; % 阵元数 d 0.5; % 阵元间距以波长为单位 theta0 30; % 扫描角度 theta -90:0.1:90; % 观测角度范围 thetaRad theta * pi/180; u sin(thetaRad) - sin(theta0*pi/180); % 阵因子等比级数封闭形式 AF sin(N * pi * d * u) ./ (N * sin(pi * d * u)); AF abs(AF); AF_db 20*log10(AF eps); plot(theta, AF_db); xlabel(角度 (deg)); ylabel(归一化方向图 (dB)); grid on; axis([-90 90 -60 0]);有几个实操细节想特别提醒。第一分母里的sin函数在u0时会趋于0代码里必须避免0/0。用eps加在log10参数上可以保证在u0处不出现NaN。第二直接用封闭式求阵因子时方向图在栅瓣位置会出现“数值上的0/0”型不定处理不当会在栅瓣处画出错误的凹坑。我的做法是把角度步长取密一点0.1°以上并用eps卫冕基本可以避开这个问题。第三阵列方向图还应当乘以阵元方向图。前面说过方向图乘积定理如果阵元是理想全向的阵元方向图就是1直接看阵因子即可。实际中阵元方向图通常有一个宽度会让大角度方向的增益进一步下降。5.2 汉明窗与副瓣抑制代码只改一个向量方向图立刻变化均匀加权方向图的副瓣很高第一副瓣−13.26dB。这一章介绍了窗函数压低副瓣MATLAB里做这件事非常简单只需要把原来全1的权向量换成窗函数向量。win hamming(N); % 或者 chebwin(N, 40) 指定副瓣电平 AF_weighted sum(win .* exp(1j * (0:N-1) * 2*pi*d*u), 1); AF_weighted abs(AF_weighted); AF_weighted_db 20*log10(AF_weighted / max(AF_weighted) eps);汉明窗可以把副瓣压到约−43dB但代价是主瓣宽度大约变成均匀加权的两倍。切比雪夫窗则可以指定副瓣电平比如30dB、40dB适合对副瓣有明确指标要求的场景。我做一个提醒设计阵列加权时副瓣和主瓣宽度是跷跷板压制副瓣必然展宽主瓣降低角度分辨力。所以设计时要看系统最优先指标是什么。如果是测角精度和分辨力优先考虑低副瓣但要接受主瓣展宽如果是抗干扰优先可能需要更深的副瓣抑制。5.3 相位量化误差10位还是12位移相器划算算过才知道数字相控阵里移相器不是连续可调的而是量化控制的。B位移相器能把360°分成2^B等份最小相位步进是2π/2^B最大量化误差大约是π/2^B。量化误差会在方向图上产生两类影响一是主瓣指向出现微小偏差二是产生量化副瓣。对均匀线阵做粗略估算B位的量化副瓣电平大约在L_q ≈ 20·log10(π / (2^B · √N))这个公式并不严格但可以用来说明规律位数的提升对副瓣改善是线性dB级别的而阵元数目N增加则摊薄了量化误差的贡献。以N32为例6位移相器量化副瓣大概在−20多dB8位能到−35dB量级10位几乎可以忽略。所以工程里选位移相器要在成本、功耗和副瓣指标之间做权衡。做低成本的车载毫米波雷达6位移相器很常见做高性能相控阵一般选到10位或更高。这里其他博主推荐直接上高位数但我的实际经验是先算清楚系统需要的副瓣水平再选位数别为了纸面上的美观多花钱。6. 从第15章看毫米波雷达相控阵思想在新硬件上的延续6.1 毫米波雷达阵列从“看得见”到“看得清”的演进第15章里的模型一句话概括就是用一组可控相位的阵元合成一个可扫描的窄波束。这个思想落在毫米波频段就是今天车载雷达、机器人雷达和工业雷达的共同底座。77GHz车载雷达的波长大约3.9mm。这个频段的好处是波长小同样的阵元数量可以用很小的物理尺寸形成足够窄的波束。比如一个8阵元线阵d0.5λ总孔径只有1.5cm左右但法线方向的3dB波束宽度已经接近12.7°。配合多片级联和虚拟孔径扩展角度分辨力可以做到1°甚至更小。我在读AWR2243这类毫米波雷达芯片的用户手册时发现芯片内部的BFPbeamforming模式本质上就是这一章的工程实现各接收通道加权求和应用不同的相位梯度实现数字波束扫描。手册里写得虽然复杂底层思想与第15章一脉相承。6.2 4D毫米波雷达与MIMO阵列相控阵思想在新硬件上的延续4D毫米波雷达这几年很热它的“4D”指的是距离、速度、水平角、俯仰角四维信息。要实现俯仰维测角不能只做一维均匀线阵需要在二维平面上布阵。这里用到的就是平面相控阵的推广方向图乘积定理仍然成立只是阵因子变成二维求和。但单纯做二维相控阵阵元数量会很庞大成本过高。于是出现了MIMO阵列的巧妙替代发射阵元之间做正交波形编码接收端通过匹配滤波把虚拟出“发射-接收联合阵列”的方向图其等效孔径远大于实际物理阵列。这可以理解为阵因子概念的扩展——不是所有阵元都在同一时间发射而是在时间和编码维度上获得了更大的自由度。热词里大量出现“awr2243雷达数据读取”也和这个趋势有关。拿到芯片首先把原始ADC数据读出来对多通道数据做距离维FFT、多普勒维FFT再对角度维做数字波束形成本质上就是这本书第15章内容在真实硬件上的复现。数据读出来之后的处理比很多人想象中更接近教材模型。6.3 数字波束形成把模拟移相器变成纯数学运算这一章最后需要延伸到的内容是数字波束形成。传统相控阵用模拟移相器在射频端完成相位控制数字波束形成则把每个阵元的基带数据全部采样到数字域再用矩阵乘法为不同方向同时生成多组波束把一个窄带波束扫描问题变成纯数学问题。这样做有三个明显优势。第一可以同时形成多个波束分别指向不同方向相当于一部雷达同时具备多个“虚拟雷达”。第二可以在数字域实现自适应波束形成用算法把干扰方向的接收信号加权置零。第三没有模拟移相器的量化误差问题相位精度只受ADC位数和计算精度限制。但代价也摆在那里每个阵元都需要独立的接收通道包含低噪声放大、混频、滤波、ADC通道数一多系统成本、功耗和数据率都会成倍上升。所以选择模拟波束成形、数字波束成形还是混合波束成形本质上是在硬件复杂度和系统灵活性之间做取舍。对做仿真的同学来说用MATLAB验证数字波束形成的效果非常简单改几行矩阵运算就能对比不同权重的方向图差异这一步跑通之后对第15章的理解会上一个台阶。我在读这一章时最大的体会是它虽然写在“特殊主题”分类里但几乎把前面所有章节的信号处理思想串起来了。雷达的检测、跟踪、成像最终都依赖精准的波束控制而波束控制的数学核心就是这一章从波程差开始一点点推出的阵因子。把阵因子公式自己在MATLAB里复原一遍比单纯追着书上的公式看要有用得多。后面如果再碰到相控阵相关的新概念比如MIMO虚拟孔径、自适应调零、频控阵回来看这一章的基本模型都能找到熟悉的影子。