简介理解阵列天线的方向图特性是无线通信与雷达系统设计的基础。单元个数、阵元间距和波长共同决定方向图形状为此整理了一套可运行的 MATLAB 仿真脚本通过对比不同配置下的波束与旁瓣变化帮助学习者建立直观认知。压缩包内含4个 .m 文件总大小仅2KB代码风格简洁覆盖阵元数量、阵元间距、波长等关键参数支持在 MATLAB 中直接运行并绘制方向图曲线。资源已有2129人次浏览学习特别适合天线入门者借助脚本快速验证概念将抽象公式转化为可视化结果。运行脚本后可对比不同单元数、不同阵元间距及不同波长条件下方向图的主瓣宽度、副瓣电平与零点分布清晰呈现栅瓣、零深和波束指向等细节脚本参数便于修改可自行扩展不同阵列构型快速生成所需仿真图为课程设计、综合实验或论文写作提供有效参考。 最近一直在做阵列方向图相关的验证工作趁着周末把几种常见构型的对比思路整理出来。做天线阵列的人应该都有体会方向图这东西单独看一副图容易自我感觉良好一旦把线阵、平面阵、圆环阵放在一起对比很多设计上的取舍才真正暴露出来。这篇文章我会从方向图乘积定理出发把均匀线阵、均匀平面阵和圆环阵的阵因子特征拆开讲再结合栅瓣抑制、幅度加权、相位量化这些工程里绕不开的问题最后顺带聊聊MIMO虚拟阵列和毫米波大规模阵列带来的新变化。内容适合刚接触阵列信号处理的学生也适合需要用方向图仿真结果支撑方案选型的系统工程师。1. 方向图比较前必须搞清楚的核心逻辑1.1 方向图到底是什么阵列方向图描述的是阵列对不同来波方向的响应增益。单天线有单元方向图阵列因为多个单元在空间中的位置不同接收或辐射信号时存在路程差这个路程差会带来相位差从而在某个方向上形成相干叠加或相消最终得到一个与角度强相关的响应曲线或响应曲面。在工程比较中我们很少直接盯着三维方向图发呆而是习惯切出方位面或俯仰面的二维切片来看。主瓣宽度、副瓣电平、栅瓣位置、零点深度这些指标才是有实际意义的对比维度。阵列方向图本质上是空间滤波器的频率响应只不过自变量从频率换成了空间角度。1.2 方向图乘积定理的工程意义方向图乘积定理是阵列分析最基础也最好用的工具。它说阵列的总方向图等于单元方向图乘以阵因子。单元方向图就是单个天线单元的方向图阵因子则完全由阵元的空间位置、激励幅度和相位决定与单元类型无关。这个解耦实在太重要了。假设单元方向图是各向同性的总方向图就等于阵因子那么不同阵列构型的差异就浓缩在阵因子这个数学式子里。工程上做方案对比时只要抓住阵因子就抓住了阵列的本质。比如均匀线阵的阵因子是一个Dirichlet核形式均匀平面阵的阵因子在两个维度上分别可分离圆环阵的阵因子则涉及贝塞尔函数的叠加。三者数学形式不同方向图特征自然分化。我第一次用MATLAB写等幅均匀线阵方向图时验证过乘积定理把一个有方向性的贴片单元方向图乘上理想阵因子得到的合成方向图和全波仿真软件里阵列建模的结果很接近只有副瓣区域有小幅偏差。偏差来源是单元的互耦这个后面会提到。1.3 单元方向图与阵因子解耦的好处解耦意味着设计流程可以分层。先定阵元间距和布阵方式计算阵因子看栅瓣和主瓣宽度能不能接受再选单元类型和单元方向图看单元增益能不能覆盖目标扫描范围。两者都满足后再做全波仿真效率会高很多。对比不同阵列构型时也只需要关注阵因子。如果要比较一个8元线阵和一个8元圆环阵单元都是同一个贴片天线那么单元方向图部分完全相同差异完全来自阵因子。这样省去了重复建模的工作量也让比较对象之间的变量更干净。2. 常见阵型的方向图特征对比2.1 均匀线阵一维窄波束的经典之选均匀线阵ULA的阵因子比较简单。对于N个等间距d排列在x轴上的各向同性辐射元等幅同相激励时阵因子绝对值是一个sin(Nψ/2)/sin(ψ/2)的形式其中ψ k·d·sinθ βk是波数β是相邻单元间的激励相位差。通过调控β可以实现波束在θ方向扫描这也是相控阵的基本原理。线阵的方向图有几个基本特征。第一垂直于阵列轴线的方向即侧射方向主瓣最窄波束宽度近似为0.886λ/(N·d)弧度N·d就是阵列总长度所以阵列越大波束越窄。第二随着扫描角从侧射方向往端射方向移动主瓣会逐渐展宽增益下降这个展宽大致与1/cosθ成正比。第三等幅激励时第一副瓣电平固定为-13.26dB左右不随阵元数变化这是很多应用不能接受的所以才有了后面要说的幅度加权。线阵明显的短板是只能在一个维度上控制波束。扫描方向一旦接近端射方向主瓣展宽得厉害栅瓣风险也急剧上升。2.2 均匀平面阵二维扫描与立体波束把线阵在另一维度上复制扩展就得到均匀平面阵UPA。平面阵的阵因子可以看成两个正交方向线阵因子的乘积。如果阵列是矩形的行列方向的阵元间距分别为dx和dy那么总方向图相当于一个沿x方向扫描的线阵和一个沿y方向扫描的线阵相乘的结果。这里需要留意共址均匀平面阵列的概念。共址说的是所有阵元在物理上位于同一块孔径内而非分布在多个分布式节点上。这种阵列的优势在于能同时形成方位和俯仰方向的窄波束对空域干扰抑制能力比线阵强很多因此在雷达和5G毫米波基站里是主流构型。平面阵方向图对比线阵时最大的差异在于波束是“针状波束”而不是“扇形波束”。这个特性决定了平面阵更适合同时需要方位和俯仰分辨力的场景但代价是阵元数量大幅增加。16×16的UPA就需要256个通道每个通道都要独立的移相器和增益控制系统复杂度直接上升一个量级。2.3 圆环阵全方位扫描与对称性优势圆环阵UCA把阵元均匀分布在圆周上这个构型在360度全方位扫描场景里很讨喜。它的阵因子没有可分离的简单解析式通常需要借助贝塞尔函数的级数展开来分析。均匀圆环阵的阵因子在水平面内理论上具有旋转对称性也就是说方向图不随方位角变化这在线阵和平面阵里是不可能实现的。实际用圆环阵时要注意虽然水平面方向图对称但是它的垂直面方向图还是会出现从平面阵到线阵之间的过渡形态。另外圆环阵的主瓣宽度和副瓣电平与阵元数量、圆的直径直接相关。如果环的直径太小阵元之间的耦合会很严重方向图会出现畸变单纯增加阵元数量不一定能改善副瓣可能还会因为互耦增强而变差。我做圆环阵对比时踩过一次坑最初为了追求对称性把阵元从8个加到16个环直径没变结果高频段方向图副瓣不但没降反而因为单元间距太近互耦抬高了副瓣。后来把直径做到一个波长以上问题才缓解。2.4 三种阵型方向图特性速查阵型阵因子特征主瓣宽度扫描范围系统复杂度典型应用均匀线阵Dirichlet核形式窄仅一维一维扫描端射区展宽低雷达测向、波束扫描均匀平面阵可分离的二维线阵因子窄两维均可控方位/俯仰二维扫描高相控阵雷达、毫米波通信圆环阵贝塞尔函数级数随方位角基本不变360度全向覆盖中高电子侦察、卫星通信终端表格只是帮大家快速建立一个认知框架真实工程里还有共形阵列、稀疏阵列、随机阵列等构型。稀疏阵列在相同阵元数量下可以获得更大的物理孔径但副瓣会升高这也是对比时经常遇到的另一个维度。3. 影响方向图的关键参数与工程陷阱3.1 单元间距与栅瓣最不能妥协的红线栅瓣是阵列方向图里最让人头疼的问题。当阵元间距过大方向图在除了主瓣以外的某个方向会出现和主瓣幅度几乎相等的峰值这时阵列就会对那个方向来的信号形成“假响应”在雷达上就是目标角度模糊在通信上就是干扰抑制失效。为了避免栅瓣需要让相邻阵元的相位差在整个可见空间范围内尽量不超过一个周期。对于均匀线阵不出现栅瓣的条件是单元间距d不超过λ/2。这个结论是常规波束形成下的必要条件。如果阵列要扫描到θ_max那条件会更严格最大的允许间距是λ/(1|sinθ_max|)。很多教材里会写成d小于等于λ/(1sinθ_max)前提是只在单侧扫描实际做仿真时最好保留安全余量。我记得有一次仿真一个32元的毫米波阵列间距取了0.6λ想着扫描范围只有正负30度应该不会出问题。结果在扫描到接近30度时方向图里在约-50度位置冒出一个副峰幅度只比主瓣低几个dB。后来算了一下0.6λ在扫描30度时刚好触碰到栅瓣产生的临界条件稍微加一点波纹误差就会越界。从那以后我设计阵列的第一件事就是先用脚本画出“单元间距-最大扫描角-栅瓣位置”的关系曲线再决定间距取值。3.2 幅度加权副瓣与主瓣的零和博弈等幅激励时线阵的第一副瓣约-13.26dB很多场景不够用。降低副瓣最直接的方式是幅度锥削让边缘阵元的激励幅度比中心阵元低。工程中常用的是切比雪夫加权和泰勒加权。切比雪夫加权可以把副瓣压到指定的电平比如-30dB或-40dB代价是主瓣会展宽阵列增益也会有损失。副瓣压得越低主瓣越宽增益损失越大这个trade-off不可回避。泰勒加权是切比雪夫的一个变体它不像切比雪夫那样让所有副瓣都相等而是让远离主瓣的副瓣逐渐衰减这样在方向图远角区域更自然工程上更常用。做方向图对比时光比“等幅激励”下的方向图是不够的。真实系统里一定会有加权加权方式不同阵型之间的相对优劣势也可能反转。比如线阵经过-40dB切比雪夫加权后主瓣已经明显变宽如果还想保持窄波束就必须在物理上加大孔径这时平面阵低副瓣和二维聚焦的优势就更突出。3.3 相位量化误差数字移相器带来的寄生瓣数字移相器的相位不是连续的而是按比特数量化的。一个4位移相器的相位步进是360/1622.5度6位移相器则是5.625度。相位量化会在方向图上产生量化瓣强度与量化位数相关。经验公式是量化副瓣电平大约在-6b dB附近b是移相器比特数。4比特对应约-24dB6比特对应约-36dB这只是一个粗略估计实际值要看具体激励分布。量化瓣通常对称出现在主瓣两侧而且扫描角度越大量化瓣越高。所以做毫米波阵列系统仿真时不能假设理想模拟相移必须把量化误差建模进去。我一般会在MATLAB里先对每个单元的相位做quantize处理然后看方向图变化而不是直接全波仿真这样迭代快得多。3.4 互耦与单元失效实测方向图失真的两大元凶方向图乘积定理里默认阵元是孤立存在的但实际阵列里阵元之间的互耦会让单元方向图不再一致尤其是边缘单元和中心单元的差异非常明显。在阵列扫描大角度时互耦引起的单元方向图畸变可能让实际副瓣比仿真高出2~4dB。单元失效是另一个容易被忽略的问题。功分网络故障、收发通道损坏都会导致某个阵元失去激励。个别阵元失效对主瓣宽度影响不明显但副瓣电平会显著升高而且方向图不对称。做可靠性分析时通常要做“阵元失效”的蒙特卡洛仿真统计不同失效比例下的方向图退化程度。解决互耦影响的工程手段包括在单元周围添加寄生元来缓解耦合、在馈电网络里加入校准通道、使用有源单元方向图替换孤立单元方向图。有源单元方向图的概念是把所有其他单元按实际激励状态接上匹配负载单独激励某一个单元测综合方向图然后再用这个方向图作为阵因子里的“单元方向图”能明显缩小仿真和实测之间的差距。4. 从静态阵列到现代演进方向4.1 MIMO虚拟阵列用少量阵元“拼出”大孔径传统相控阵的每个阵元都要配一路独立的收发通道成本和功耗都很高。MIMO雷达和MIMO通信系统用另一种思路解决这个问题发射端和接收端同时使用多个阵元但发射阵元和接收阵元分开布置利用发射-接收阵元对之间的相位关系可以合成一个比实际物理阵元数量大得多的虚拟阵列。以单发多收为例一个发射阵元和N个接收阵元可以等效成一个N元阵列。如果发射端有M个阵元接收端有N个阵元那么接收信号里能区分出M×N个不同的相位中心这就是虚拟孔径。方向图方面MIMO阵列的等效方向图等于发射方向图和接收方向图的乘积这个关系和方向图乘积定理有异曲同工之妙只不过现在的“单元方向图”被换成了整个发射阵列的阵因子。实际做MIMO虚拟阵列的方向图分析时最怕的是虚拟阵元位置出现重叠或空洞。比如发射阵元间距和接收阵元间距选择不当虚拟阵元会在某些位置重合等效阵列就不是均匀的方向图会出现不期望的高副瓣。设计时通常要保证收发间距满足互质关系这样虚拟阵元填充得最均匀。4.2 毫米波大规模阵列方向图仿真从理想走向OTA实测毫米波频段波长很短同样物理尺寸里可以塞进成百上千个阵元这也让大规模天线阵列建模成为热点。大规模阵列的分层架构很常见先做子阵级波束成形再做子阵间的波束成形这样总方向图就变成了子阵方向图和阵列因子的乘积。这个分层结构有个工程便利阵列方向图的仿真可以先在数学层面完成用方向图乘积定理快速迭代优化子阵划分方式不用每改一次子阵布局就跑全波仿真。但毫米波阵列的实测非常依赖空口测试OTA因为常规的传导测试无法覆盖所有阵元的辐射特性OTA暗室的校准和测量误差对方向图结果的影响比微波频段更敏感。4.3 从超声到麦克风方向图比较方法的跨界复用行列寻址阵列超声成像、麦克风阵列这些场景用的数学工具和天线阵列其实是一样的。天线里叫栅瓣超声里叫伪像天线里叫阵列方向图麦克风里叫波束图。核心都是通过阵元位置和激励权重的设计来控制空间响应。麦克风阵列里方向图乘积定理同样成立只是声波波长远大于电磁波阵元间距通常是厘米到分米级别最大扫描角受限更明显。做系统开发时如果已有天线阵列的仿真代码稍微改一下波长、阵元数和工作频段就能直接用于麦克风阵列的波束图分析。这种跨领域复用省下的时间不是一点半点。4.4 仿真工具建议MATLAB的Phased Array System Toolbox是计算阵因子最顺手的工具内置了多种阵列对象包括ULA、URA和UCA可以直接绘制方向图、计算副瓣电平还能方便地做幅度加权和相位量化仿真。HFSS和CST这类全波软件能精确模拟互耦和边界效应但建模和计算代价高。我的习惯是先用MATLAB或Python快速算阵因子选出两三种有潜力的构型再做全波仿真验证。Python里也可以用NumPy自己写阵因子计算七八十行代码就能实现一个可交互的线阵方向图界面对理解原理帮助很大。5. 常见问题与排查技巧实录5.1 方向图出现意外的“大副瓣”这个最常见的原因就是单元间距超过了无栅瓣条件。先检查扫描角度再看阵列单元间距。如果间距没问题再看激励幅度是否因为加权系数输入错误导致边缘单元异常。还有一种可能是坐标系定义不对导致阵元位置计算错误尤其平面阵的dx和dy混淆。排查步骤一般是先画等幅激励下的方向图确认栅瓣位置再叠加上加权看副瓣变化趋势最后才检查移相量化误差和互耦。如果等幅激励下副瓣就已经异常高那么大概率是阵列几何建模阶段出了问题。5.2 切比雪夫加权设置了-40dB实测只有-33dB这种偏差在工程里非常常见。切比雪夫加权的前提是理想等幅同相单元、各向同性单元归一化阵因子。实际单元的单元方向图在不同角度上增益不同尤其是远离法线的大角度区域单元增益掉得快此时的副瓣相对电平会发生变化。另外功分网络和移相器本身的插损不一致也会打破理想权重的比例关系。优化方向是先用实测的单元方向图或者有源单元方向图替代理想单元方向图在计算加权系数时把单元增益的随角变化考虑进去。这就引出了“方向图加权综合”这一层话题单纯依赖切比雪夫公式套模板是不行的。5.3 平面阵扫描时方位面和俯仰面的波束宽度不对称如果平面阵的行列间距不同比如dx0.5λ、dy0.7λ那么方位和俯仰方向的有效孔径长度不同波束宽度自然不一样这是正常的。如果行列间距相同但波束宽度仍然差异明显就要检查是不是行列方向阵元数不一致或者加权系数在行列方向上用了不同的窗函数。还有一种容易被忽略的情况双极化平面阵中极化方向与扫描面的关系也会影响方向图切片的形状。水平极化和垂直极化在斜切面上看到的增益不同对比方向图切片时最好保证极化匹配。5.4 阵列方向图仿真实测差异速查现象可能原因检查措施副瓣整体偏高2~4dB互耦未建模、单元方向图不一致改用有源单元方向图特定方向出现尖峰栅瓣或量化瓣检查间距和移相器分辨率方向图左右不对称单元失效、馈电网络幅相误差检查各通道幅相校准状态主瓣宽度大于仿真阵元相位误差增大、幅度锥削过深检查波控码和加权系数扫描大角度时增益塌陷有源反射系数恶化、单元匹配失调校准失配、检查功分器一致性5.5 一个值得建立的习惯先画“栅瓣边界图”不论做线阵还是平面阵强烈建议在项目初期画一张单元间距与最大扫描角的关系图。横轴是扫描角纵轴是最大允许间距标注出无栅瓣区域。这张图贴在工位上每次调整阵列参数之前先看一眼能省去大量后发现栅瓣而返工的时间。对平面阵来说还要注意斜切面的等效间距问题。波束指向不是沿主轴方向时在斜切面上看到的阵元投影间距可能会大于主轴间距同样会触发栅瓣。所以设计平面阵的扫描范围时往往以扫描锥角的最差情况来约束间距而不是只看方位或俯仰单维度。6. 最后分享一点个人体会做阵列方向图对比这么久最大的感受是不要把方向图当成一个孤立的仿真结果而要把它和系统指标串起来看。副瓣电平会影响抗干扰能力栅瓣会直接造成测角模糊主瓣宽度决定了空间分辨力这些指标背后都是阵列构型和参数的选择问题。方向图只是一个映射真正要理解的是阵列的物理孔径、单元间距、加权策略、相位量化这些底层设计变量之间的关系。另外一个小建议方向图的对比不必一上来就追求全波仿真精度。先用方向图乘积定理把理论阵因子算清楚确认栅瓣位置、副瓣水平和主瓣宽度是否符合预期再做全波仿真验证互耦影响。这个顺序能节省大量时间也更容易定位问题出在阵列几何设计还是单元设计上。我自己现在做新阵列方案时依然会先打开脚本调理想阵因子看到曲线符合预期了才会去碰全波仿真软件。本文还有配套的精品资源点击获取