简介面向天线工程、通信工程专业学生以及阵列信号处理初学者的技术文档系统讲解阵列天线方向图的波束形成原理与MATLAB仿真方法。文档从阵列天线辐射场的矢量叠加概念出发引入方向图相乘原理说明多元天线阵总方向图函数的简化计算思路为理解相控阵波束扫描与栅瓣抑制提供理论基础。压缩包内包含1个doc文档大小约97KB内容涵盖实验目的、理论原理、三组完整源程序与结果分析结构完整可直接用于课程学习或实验参考。三组源程序分别考察阵列个数n、波长lamda与阵元间距d对方向图的影响并给出对应的方向图曲线与结论可直观看到主瓣宽度、旁瓣电平等参数随变量变化而改变的趋势。文档还结合仿真结果分析了各参数对天线性能的作用帮助读者掌握从理论推导到程序验证再到结果分析的完整研究思路。目前已有345人学习下载对需要掌握阵列天线仿真技巧、完成课程设计或撰写实验报告的读者有较高参考价值。1. 阵列天线方向图及其MATLAB仿真为什么多数仿真结果不能直接用很多做天线方向图仿真的同学第一次跑通代码、看到极坐标图里那条主瓣曲线时心里都会松一口气方向图出来了作业应该能交了。但等到拿仿真方向图去反推波束宽度、副瓣电平或者换一种阵元间距重跑一遍时才发现结果和理论值对不上甚至出现完全多余的栅瓣。我刚接触阵列天线方向图及其MATLAB仿真这个课题时也犯过同样的错——把均匀直线阵的公式一股脑写进循环画出一张“看起来对”的图却没有检查边界条件。真正的问题不在MATLAB能不能算而在于你用什么模型去算、阵元间距和相位差怎么设定、方向图乘积定理在什么前提下才能成立。这篇文章从最常见的线阵出发给你一条能直接复现的仿真路径先立原理再写代码最后用几个典型参数验证结果是不是物理可信的。2. 方向图计算的底层逻辑阵列因子、阵元因子与方向图乘积定理2.1 线阵方向图的解析公式从等幅等相到幅度加权在动手写MATLAB之前我想先花一点篇幅把方向图计算的来龙去脉说清楚。很多人直接拿别人的仿真脚本改参数改出来的图“长得像”但不知道每个量在公式里对应什么。阵列天线的总方向图通常写成阵元因子与阵列因子的乘积这个结论被称为方向图乘积定理。它的物理含义很直白如果你有一排相同的阵元每个阵元本身有方向性那么整个阵列的远场方向图就是这个阵元的方向图“乘上”一组离散点源按同样排列方式形成的阵列因子。前者描述单个辐射单元在空间各方向的增益变化后者描述阵列排布带来的干涉效果。先看最常见的情形沿z轴等间距排列的N个全向点源间距为d。第n个阵元的远场相位差来自两部分一部分是激励电流的初始相位α另一部分是空间程差kdcosθ。把这个相位求和就得到阵列因子的复表达。等幅激励下如果你用等比数列求和公式化简可以得到一个非常经典的闭式解也就是含sin(Nψ/2)/sin(ψ/2)的那种形式。这个形式虽然看起来简单却隐藏着三个关键判断点第一波束最大值出现在ψ0处也就是各阵元辐射同相叠加的方向第二当ψ/2等于π的整数倍时分子分母同时为零需要用极限求值而这一位置的周期性正是栅瓣的来源第三对N取极限后你还能从这个表达式估算第一副瓣相对于主瓣的电平也就是约-13.5dB。在MATLAB里做仿真常见做法是不要直接把闭式解硬编码而是按照定义逐项求和。好处是当你后面引入幅度加权、相位加权时只需要改权值向量不用改公式结构。以俯仰角θ作为自变量阵列因子写成复数累加的形式角度范围扫-90°到90°步长用0.1°。这个做法虽然计算量比闭式解大一点点但对学习阶段来说逻辑更清晰也便于后续扩展到平面阵。2.2 为什么不能忽略阵元因子微带贴片与全向振子的差别很多初学仿真的人会把阵元当成全向点源这样阵列因子就是总方向图。这个假设在线元天线、单极子天线的粗略估算里勉强说得过去但一旦把阵列因子换成微带贴片天线就会发现方向图差距惊人。某个项目中我用8元线阵做波束扫描仿真最初按全向阵元算扫描到60°时增益曲线还是正常的主瓣加副瓣后来把阵元因子改成微带贴片在E面的方向图同样的阵列因子乘上去扫描增益掉了好几个dB第一副瓣也发生了变化。方向图乘积定理在这里的教训是阵元因子不是可选项而是方向图的一部分。具体实现时阵元因子通常写成一个与θ相关的函数可以来自理论公式也可以来自全波仿真软件导出的数据。最省事的做法是用理论公式。例如半波偶极子在某个平面内的方向图因子近似为cos(π/2cosθ)/sinθ而微带贴片则可以用腔模理论的近似公式。两者乘上同一个阵列因子后主瓣宽度、副瓣起伏都会有明显不同。所以你在报告里写明“本次仿真采用全向阵元模型”和“采用贴片阵元模型”结论的适用范围是完全不同的。在代码结构上我一般会把方位角φ也纳入计算因为很多阵列天线仿真需要对整个上半空间画三维方向图。经典做法是用meshgrid生成θ和φ的二维网格然后分别计算每个角度对应的阵元因子与阵列因子。但初学者建议先做一维切片也就是固定φ0°或φ90°只扫描θ。这样能更直观地看到主瓣、副瓣和栅瓣的变化规律也方便和解析公式对照。等一维切片验证通过再把二维网格套进来。3. 用MATLAB算均匀直线阵方向图从一行行代码到能用的函数3.1 最小可复现的主程序参数、循环与极坐标绘图我先把一个能直接跑通的均匀直线阵方向图仿真脚本拆给你看。这个脚本覆盖了最核心的流程设置阵列参数、生成空间角度网格、计算阵列因子、取模归一化、画极坐标图。我习惯把参数全部放在脚本开头方便后续改数。% 均匀直线阵方向图仿真等幅激励 clear; clc; close all; % 基本参数设置 N 8; % 阵元数量 d 0.5; % 阵元间距单位波长lambda theta0 0; % 波束指向角度度0为侧射阵 % 角度扫描范围 theta linspace(-pi/2, pi/2, 1801); % -90°到90°步长0.1° % 波数向量k*d*sin(theta) beta 2 * pi * d * sin(theta); % 空间相位差未含阵元初始相位 % 若波束扫描阵元激励相位为 -k*d*sin(theta0)写为复数权值 % 这里等幅激励各阵元相位均为0权值向量全1 w ones(N, 1); % 计算阵列因子 AF zeros(size(theta)); for n 0:N-1 AF AF w(n1) * exp(1j * n * beta); end % 归一化方向图用最大值归一化取20log10转dB AF_norm abs(AF) / max(abs(AF)); AF_dB 20 * log10(AF_norm eps); % 加eps避免log10(0) % 绘图半空间极坐标 figure; polarplot(theta, AF_norm, LineWidth, 1.5); title(sprintf(均匀直线阵方向图 N%d d/lambda%.2f, N, d));这段代码的逻辑很简洁每个阵元的相位贡献分成n倍的空间相位差权值w控制幅度和初始相位。等幅激励时w全为1所以AF就是等比数列求和。角度网格取了1801个点保证曲线平滑。极坐标图里横坐标是θ纵坐标是归一化后的幅度。注意这里你看到的副瓣高度、零点位置和公式推导完全一致可以直接拿来做解析解对照。参数变更时的行为值得仔细说。阵元间距d用波长作单位这是方向图仿真最核心的尺度参数。d0.5表示间距为半个波长此时-90°到90°范围内不会出现栅瓣是多数阵列设计的基础工况。如果把d改为1.0你会在±90°附近看到幅度与主瓣相等的额外波瓣这在物理上意味着天线会在多个方向同时产生强辐射工程上通常要避免。波束指向theta0参数在等幅激励下先不改后面做波束扫描时它才会真正起作用。3.2 副瓣电平、半功率波束宽度和零点的精确提取画出一张漂亮的极坐标图只是第一步。在实际工程汇报里你需要从方向图数据里提取三个关键指标半功率波束宽度HPBW、第一副瓣电平SLL、主瓣零点位置。用MATLAB做这件事时我倾向于不依赖图形界面手动读数据而是写一个简单的数值提取流程这样参数一改指标自动更新。% 从方向图数据提取关键指标 % 输入theta弧度或度AF_norm归一化幅度0~1 % 输出HPBW、第一副瓣电平、零点位置 % 1. 找到主瓣峰值位置 [maxVal, maxIdx] max(AF_norm); mainlobe_dir theta(maxIdx) * 180 / pi; % 2. 半功率点幅度下降到0.7071即-3dB处的角度范围 half_power 0.7071; cross_idx find(AF_norm half_power); HPBW (theta(max(cross_idx)) - theta(min(cross_idx))) * 180 / pi; % 3. 第一副瓣在主瓣两侧各找局部最大值 % 先把主瓣区域去掉取主瓣两侧区间 left_region AF_norm(1:maxIdx-10); right_region AF_norm(maxIdx10:end); % 用findpeaks找局部极大值取绝对值最大的那个 [pks_left, locs_left] findpeaks(left_region); [pks_right, locs_right] findpeaks(right_region); if ~isempty(pks_left) leftSLL_dB 20 * log10(max(pks_left)); else leftSLL_dB -inf; end if ~isempty(pks_right) rightSLL_dB 20 * log10(max(pks_right)); else rightSLL_dB -inf; end fprintf(主瓣方向: %.2f°, HPBW: %.2f°, 左第一副瓣: %.2f dB, 右第一副瓣: %.2f dB\n, ... mainlobe_dir, HPBW, leftSLL_dB, rightSLL_dB);这段提取代码有一点需要注意findpeaks会找出所有局部极大值包括方向图两边靠近±90°的边界隆起。如果阵列因子在扫描范围内没有出现完整的第一副瓣findpeaks的返回值可能包含边界伪峰判断时要看返回的幅值是否显著低于主瓣。另一方面均匀直线阵的第一副瓣理论值是-13.46dB如果你算出来是-13.0dB左右说明网格步长不够密峰值采样偏低把角度步长从0.1°加密到0.01°能改善这个误差。这是很典型的采样精度问题不是代码逻辑错误。在半功率波束宽度的计算里我用的是“幅度大于0.7071的角度范围”这个判据。这个做法对单主瓣方向图精确有效但如果你把波束扫描到大角度方向图可能只在某一侧有完整主瓣另一侧被扫描盲区截断这时候HPBW会偏大甚至无法计算。遇到这种情况我通常会改用局部寻根的方式在主瓣峰值两侧分别找幅度等于0.7071的角度点再取差值这样即使两侧不对称也能得到稳定的结果。4. 波束扫描与幅度加权从均匀阵到低副瓣设计4.1 相位加权实现波束扫描为什么大角度扫描会变形均匀直线阵的波束指向可以通过改变相邻阵元间的激励相位差来控制。你只需要给第n个阵元乘上一个相位因子使各阵元在目标方向θ0上同相叠加。把前面的脚本稍微改一下权值从全1变成含相位差的复数向量就能看到主瓣跟着θ0移动的效果。% 波束扫描权值叠加相位梯度 N 16; d 0.5; theta0 30; % 期望波束指向度 % 角度网格 theta linspace(-pi/2, pi/2, 1801); beta 2 * pi * d * sin(theta); % 相邻阵元相位差-k*d*sin(theta0) phase_shift -2 * pi * d * sind(theta0); % 权值向量 n (0:N-1); w exp(1j * n * phase_shift); % 阵列因子 AF zeros(size(theta)); for idx 1:length(beta) AF(idx) w. * exp(1j * n * beta(idx)); % 向量化内积 end AF_norm abs(AF) / max(abs(AF)); figure; polarplot(theta, AF_norm, LineWidth, 1.5); title(sprintf(波束扫描 theta0%d°, N%d, d%.1fλ, theta0, N, d));这段代码把原来的循环累加改写成了向量内积这是一个值得养成的习惯权值向量和“相位响应向量”做内积既快又不容易下标出错。这里的关键参数是phase_shift它决定了波束指向。请记住这个公式相邻阵元相位差等于-2π/λ·d·sinθ0。你用度数制sind函数时千万别把0.5波长直接乘进sind单位混淆是这类仿真最常见的低级错误。大角度扫描的物理现象在这段代码里体现得很清楚。当θ0从0°增加到60°时主瓣会变宽峰值增益下降而另一侧的副瓣会抬升。这是阵列口径投影效应造成的波束偏离侧射方向后阵列在视线方向上的有效口径变小方向图被展宽。工程经验是扫描角超过60°后均匀直线阵的方向图质量明显恶化所以很多相控阵设计把扫描范围限制在±60°以内。你仿真时会发现θ060°的主瓣宽度几乎是0°时的两倍这不是脚本写错了而是物理规律本身如此。4.2 切比雪夫与泰勒加权低副瓣的代价是主瓣变宽均匀激励的副瓣电平约-13.5dB很多实际系统要求副瓣低于-30dB那就需要引入幅度加权。切比雪夫加权在给定副瓣电平下能达到最窄的主瓣宽度泰勒加权则是它的近似实现避免切比雪夫加权在阵元数较多时出现两端阵元激励幅度跳变过大的问题。在MATLAB里做这两种加权最常见的路径是调用通信工具箱的函数例如chebwin。如果你没有工具箱自己写切比雪夫加权也不是不行但需要用到切比雪夫多项式的递推关系代码量会大不少。% 切比雪夫加权指定副瓣电平 N 16; side_lobe_db -30; % 希望达到的副瓣电平 w_cheb chebwin(N, abs(side_lobe_db)); % 注意参数是正数 % 泰勒加权需要自行生成省略通常用于差方向图或更大阵列 % 把w_cheb代入前面方向图计算流程比较即可我用这个函数时踩过一次参数坑chebwin第二参数必须是非负的dB值如果你按习惯传-30MATLAB直接报错。此外加权后方向图的主瓣会比均匀阵略宽副瓣被压下去以后换来的代价就是角度分辨率的损失。做仿真对比时一定要在报告里同时画出均匀阵和切比雪夫阵的方向图并标出两者的HPBW和SLL否则别人看不出你有没有理解这个权衡。泰勒加权在MATLAB里没有直接的内置函数常见做法是自己写一个基于泰勒线源分布的近似生成函数。它的思路是用一组“准副瓣”替代切比雪夫的等副瓣条件让靠近主瓣的副瓣逐渐衰减这样阵元激励幅度分布更平滑。实际仿真结果显示泰勒加权的副瓣包络不是完全平坦的而是随角度增大缓慢下降这对天线加工和馈电网络设计比较友好。你只需要记住结论切比雪夫适合小阵列、严格副瓣指标泰勒适合大阵列、工程可实现性优先。5. 平面阵与方向图乘积定理的推广把一维结论延伸到二维5.1 矩形平面阵的阵列因子分解技巧行因子与列因子的相乘当你从线阵转向平面阵时方向图计算的核心技巧是分解。一个M行N列的矩形平面阵如果阵元按矩形栅格排列且行间与列间的激励相位可以分离那么整个阵列因子等于行方向阵列因子和列方向阵列因子的乘积。这个乘积关系的成立条件比较严格阵元在x方向的相位梯度和y方向的相位梯度相互独立且幅度分布可以写成两个一维分布的乘积。在MATLAB中平面阵的典型做法是用二维u-v空间进行计算。定义usinθcosφvsinθsinφ然后让阵列因子成为u和v的二维函数。这样做的好处是方向图在u-v平面上呈周期结构栅瓣条件和一维情况可以直接类比。你只需要分别计算每个方向余弦对应的相位项对全部阵元求和即可。% 矩形平面阵阵列因子矩形栅格 M 8; % x方向阵元数 N 8; % y方向阵元数 dx 0.5; % x方向间距波长 dy 0.5; % y方向间距波长 % 方位角与俯仰角网格 az linspace(-pi/2, pi/2, 361); el linspace(-pi/2, pi/2, 181); [AZ, EL] meshgrid(az, el); % 方向余弦 u sin(EL) .* cos(AZ); v sin(EL) .* sin(AZ); % 阵元位置下标 m (0:M-1); % x方向 n (0:N-1); % y方向 % 阵列因子等幅激励 AF zeros(size(u)); for p 1:numel(u) AF(p) sum(exp(1j*2*pi*(m*dx*u(p) n*dy*v(p))), all); end AF_norm abs(AF) / max(abs(AF)); % 三维方向图显示用surf或mesh figure; surf(u, v, AF_norm); xlabel(u); ylabel(v); zlabel(归一化幅度);这个脚本最需要注意的地方是变量维度。m是列向量n也是列向量在求和时要确保u(p)和v(p)是标量然后利用MATLAB的隐式扩展把mdxu(p)和ndyv(p)叠加成M×N的矩阵。这里如果忘记转置维度对不上要么报错要么得到错误的结果还很难发现。平面阵的方向图在u-v空间看是最直观的主瓣是一个锥形峰栅瓣出现在u方向或v方向的整数倍周期上它们的位置由dx和dy决定。所以平面阵的栅瓣条件是dx和dy都不超过0.5λ和一维线阵的结论一致但需要两个方向同时检查。5.2 栅瓣判据d/λ的极限值别只看公式栅瓣判断是阵列方向图仿真里最容易被忽略的环节。一维均匀直线阵的栅瓣条件由d/λ决定当间距超过一个波长时可见空间内会出现第二个峰值。这个结论的物理本质是空间相位差的周期性与可见范围的有限性之间的矛盾。很多同学把公式背下来但遇到非均匀阵列、波束扫描或者平面阵时就不知道如何用公式判断了。我的排查习惯是先画一张“全空间阵列因子”图也就是不裁剪可见空间直接看阵列因子的周期性栅瓣在哪里。用MATLAB把β从-4π扫到4π画出的阵列因子图像显示出一串等高的峰值相邻峰值间隔为2π。可见空间对应β的取值范围是[-2πd/λ·sinθ, 2πd/λ·sinθ]。只要这个范围里落入第二个峰值方向图里就会出现栅瓣。这个方法比死记dλ的判断准得多特别是在波束扫描时主瓣移到θ0原来的栅瓣位置也会移动某个原来刚好在可见范围边缘的栅瓣可能在扫描后进入可见范围。按上面的判据你能推出一个重要工程结论如果要求波束扫描到±60°阵元间距必须满足d/λ小于约0.577才能保证全扫描范围内不出现栅瓣。这就是很多相控阵实际采用0.5λ~0.55λ间距的原因不是随便定的。我参与过的模拟项目X里最初为了省阵元数量用了0.6λ间距一扫描到45°就看到了栅瓣进入可见区换回0.5λ后问题消失。这个案例说明方向图仿真的价值不只是画一张好看的图而是能在设计阶段提前暴露这类硬约束。6. 阵列方向图仿真避坑指南五个反复出现的翻车点6.1 现象polarplot出现多条放射状直线方向图杂乱无章先说我自己最早翻车的一次。用polarplot画方向图时图形里出现了从圆心向外辐射的直线簇看起来像蜘蛛网完全不是平滑的方向图曲线。检查代码后发现角度数据里有NaN值。原因出在归一化这一步如果某一角度上所有阵元相位相互抵消阵列因子的模恰好为0直接除以最大值后会得到0但如果你先求了log10再归一化log10(0)会给出-inf。绘图函数遇到-inf或NaN时会把曲线画得乱七八糟。解决办法很简单两个位置要处理一是求阵列因子后对abs结果做归一化二是要画dB曲线时给幅度加上一个极小量eps比如用20*log10(AF_normeps)。另外用polarplot时如果角度向量不是严格递增且无缝覆盖整个区间也可能出现跨零点的连线。建议先把角度排序再绘图避免画出一根横穿整个图形的直线。6.2 现象阵元间距大于一个波长方向图“多了一个主瓣”但报告里没解释间距超过一个波长时方向图会出现栅瓣这是理论课上反复强调的内容但仿真时很多人看到多出来一个等高的“主瓣”就直接认为是代码写错了。实际上这是物理规律说明你的阵列已经出现了栅瓣。这个现象的出现条件是d/λ1而且栅瓣的位置可以通过sinθ_g mλ/d来估算m±1时对应第一对栅瓣。正确的处理方式不是删掉栅瓣而是在报告中明确计算并标注栅瓣位置同时说明实际天线系统会通过限制扫描范围或减小阵元间距来避免栅瓣落入感兴趣的区域。如果你在项目里发现“主瓣数量变多”先确认d/λ的值别急着改代码。有些同学为了消除栅瓣把间距调得特别小结果方向图主瓣变得非常宽阵列的增益和分辨率都下降了这就是把栅瓣问题和口径尺寸问题混为一谈。6.3 现象加窗后副瓣没有降下来反而出现不对称的杂散峰这个坑我印象最深。用chebwin给16元阵列加-30dB切比雪夫权值后方向图的左副瓣降到了-30dB左右右副瓣却还有-20dB以上的“杂散峰”怎么看都像代码出了bug。排查了很久才发现问题不在于权值函数而在于角度网格不均匀导致峰值采样偏差。方向图的副瓣位置与角度网格点不完全重合时findpeaks找到的局部最大值比真实副瓣峰值小所以左右副瓣读数不一致。解决方法是加密角度扫描网格把步长从1°改成0.01°让每个副瓣区域有足够的采样点。另一个更稳健的办法是先用插值函数对方向图做样条插值再找峰值。我看到有人用sinc函数插值恢复方向图的连续值效果也不错只是代码量略大。6.4 现象波束扫描方向图在某个角度上主瓣突然消失在仿真波束扫描从0°扫到80°时主瓣在某个角度后幅度骤降有时直接消失在噪声底之下。这通常不是阵列因子算错了而是你选的阵元因子模型在扫描角附近就是没有辐射能力。最典型的例子是阵元因子为cosθ的模型它模拟的是一种理想偶极子的方向性在θ接近±90°方向增益自然趋近于零。总方向图等于阵元因子乘以阵列因子所以在可见区边缘总方向图必然掉下去。此时你需要区分“阵列因子还在增大”和“总方向图在衰减”两个概念报告里如果只画总方向图会让人误以为扫描能力不足。我在实际项目里就遇到过一个类似的困惑某次仿真显示波束扫描到70°时总方向图几乎消失但阵列因子明明还有明显峰值。后来在报告里分开画出阵元因子、阵列因子、总方向图三条曲线问题一目了然。所以仿真时一定要保留中间变量不要只出最终结果。6.5 现象同样参数两次运行结果不一样怀疑随机性方向图本身是确定性计算不该有随机性。但如果你在代码里用了randn生成权值或者误用了未固定的随机数种子结果自然每次不同。更隐蔽的情况出现在使用某些工具箱函数时例如自适应阵列的方向图计算会涉及采样协方差矩阵估计而这个估计通常基于随机快拍数据。排查思路是第一步检查是否用了随机函数第二步检查代码中是否有排序函数对相同数据给出不同顺序第三步检查是否依赖全局变量或并行计算时的工作区分裂。对教学仿真来说建议在脚本开头固定随机种子例如rng(2024)这样别人复现你的结果时数据完全一致。固定随机种子是我在提交仿真作业前必做的一步它不算复杂但能避免掉很多“玄学复现”的争议。7. 方向图仿真的验证闭环独立检验你的结果是否可信最后一章我想给你一个可以立即操作验证流程这也是我自己做方向图仿真时最看重的一个习惯。仿真可以说服自己但要让别人信服你必须提供两条独立路径的比对。我一般会做三件事第一用解析公式手算几个关键角度的方向图数值和MATLAB输出对比第二用全波仿真软件或另一个独立计算方法对同一阵列建模比对主瓣位置和宽度第三把阵元数N和间距d取极限看方向图是否退化为已知的特例。解析对比的做法很轻量以均匀直线阵为例在θ0°时所有阵元同相叠加阵列因子应为N。你把N8代入计算abs(AF)的最大值应该恰好等于8归一化后才是1。如果最大模值不是整数或者明显偏离N说明相位项的计算有误。第二个特例是N1这时阵列因子应该恒等于1方向图是各向同性的圆。任何方向图仿真只要N1时画出方向图不是圆形相位公式肯定有问题。第三个常见的验证方法是把方向图积分看辐射功率是否守恒。虽然方向图仿真不涉及绝对功率但归一化方向图的积分在理想情况下应当与天线辐射的总功率成正比。微带阵列仿真里我经常用方向图积分估算阵列的定向性directivity如果算出的定向性和理论值例如主瓣立体角估算差得超过1dB就要怀疑是不是角度采样太粗或者网格不均匀导致的积分误差。在实际汇报场合中评审人最容易接受的数据呈现方式是一张对比表列出均匀阵、切比雪夫加权阵和扫描状态下的主瓣方向、HPBW、第一副瓣电平、栅瓣位置。这样一张表比十张方向图都更有说服力。我自己在模拟项目X里每次调整参数后都会自动生成这张表然后截图放进实验记录里。多年的习惯养下来我发现写报告时最缺的从来不是“图”而是“图对应的关键参数表”。方向图仿真这个方向入门容易精通难。很多人以为学会了画极坐标图就等于会了阵列天线其实真正的功夫在那些看不见的细节里栅瓣判据、阵元因子的选取、加权方式的权衡、采样密度对峰值提取的影响。我在这些细节上踩过的坑远不止上面五条但每一条都在后来的仿真中变成了下意识的检查项。仿真不是把公式翻译成代码而是用代码验证你对物理的假设是否正确。每次跑完一组结果多问一句“这个图物理上可能的边界条件是什么”你的仿真能力会提升得很明显。希望这篇沿着阵列天线方向图及其MATLAB仿真整理的实战笔记能帮到你也希望你在对照复现之后能画出经得起追问的方向图。本文还有配套的精品资源点击获取 简介面向天线工程、通信工程专业学生以及阵列信号处理初学者的技术文档系统讲解阵列天线方向图的波束形成原理与MATLAB仿真方法。文档从阵列天线辐射场的矢量叠加概念出发引入方向图相乘原理说明多元天线阵总方向图函数的简化计算思路为理解相控阵波束扫描与栅瓣抑制提供理论基础。压缩包内包含1个doc文档大小约97KB内容涵盖实验目的、理论原理、三组完整源程序与结果分析结构完整可直接用于课程学习或实验参考。三组源程序分别考察阵列个数n、波长lamda与阵元间距d对方向图的影响并给出对应的方向图曲线与结论可直观看到主瓣宽度、旁瓣电平等参数随变量变化而改变的趋势。文档还结合仿真结果分析了各参数对天线性能的作用帮助读者掌握从理论推导到程序验证再到结果分析的完整研究思路。目前已有345人学习下载对需要掌握阵列天线仿真技巧、完成课程设计或撰写实验报告的读者有较高参考价值。1. 阵列天线方向图及其MATLAB仿真为什么多数仿真结果不能直接用很多做天线方向图仿真的同学第一次跑通代码、看到极坐标图里那条主瓣曲线时心里都会松一口气方向图出来了作业应该能交了。但等到拿仿真方向图去反推波束宽度、副瓣电平或者换一种阵元间距重跑一遍时才发现结果和理论值对不上甚至出现完全多余的栅瓣。我刚接触阵列天线方向图及其MATLAB仿真这个课题时也犯过同样的错——把均匀直线阵的公式一股脑写进循环画出一张“看起来对”的图却没有检查边界条件。真正的问题不在MATLAB能不能算而在于你用什么模型去算、阵元间距和相位差怎么设定、方向图乘积定理在什么前提下才能成立。这篇文章从最常见的线阵出发给你一条能直接复现的仿真路径先立原理再写代码最后用几个典型参数验证结果是不是物理可信的。2. 方向图计算的底层逻辑阵列因子、阵元因子与方向图乘积定理2.1 线阵方向图的解析公式从等幅等相到幅度加权在动手写MATLAB之前我想先花一点篇幅把方向图计算的来龙去脉说清楚。很多人直接拿别人的仿真脚本改参数改出来的图“长得像”但不知道每个量在公式里对应什么。阵列天线的总方向图通常写成阵元因子与阵列因子的乘积这个结论被称为方向图乘积定理。它的物理含义很直白如果你有一排相同的阵元每个阵元本身有方向性那么整个阵列的远场方向图就是这个阵元的方向图“乘上”一组离散点源按同样排列方式形成的阵列因子。前者描述单个辐射单元在空间各方向的增益变化后者描述阵列排布带来的干涉效果。先看最常见的情形沿z轴等间距排列的N个全向点源间距为d。第n个阵元的远场相位差来自两部分一部分是激励电流的初始相位α另一部分是空间程差kdcosθ。把这个相位求和就得到阵列因子的复表达。等幅激励下如果你用等比数列求和公式化简可以得到一个非常经典的闭式解也就是含sin(Nψ/2)/sin(ψ/2)的那种形式。这个形式虽然看起来简单却隐藏着三个关键判断点第一波束最大值出现在ψ0处也就是各阵元辐射同相叠加的方向第二当ψ/2等于π的整数倍时分子分母同时为零需要用极限求值而这一位置的周期性正是栅瓣的来源第三对N取极限后你还能从这个表达式估算第一副瓣相对于主瓣的电平也就是约-13.5dB。在MATLAB里做仿真常见做法是不要直接把闭式解硬编码而是按照定义逐项求和。好处是当你后面引入幅度加权、相位加权时只需要改权值向量不用改公式结构。以俯仰角θ作为自变量阵列因子写成复数累加的形式角度范围扫-90°到90°步长用0.1°。这个做法虽然计算量比闭式解大一点点但对学习阶段来说逻辑更清晰也便于后续扩展到平面阵。2.2 为什么不能忽略阵元因子微带贴片与全向振子的差别很多初学仿真的人会把阵元当成全向点源这样阵列因子就是总方向图。这个假设在线元天线、单极子天线的粗略估算里勉强说得过去但一旦把阵列因子换成微带贴片天线就会发现方向图差距惊人。某个项目中我用8元线阵做波束扫描仿真最初按全向阵元算扫描到60°时增益曲线还是正常的主瓣加副瓣后来把阵元因子改成微带贴片在E面的方向图同样的阵列因子乘上去扫描增益掉了好几个dB第一副瓣也发生了变化。方向图乘积定理在这里的教训是阵元因子不是可选项而是方向图的一部分。具体实现时阵元因子通常写成一个与θ相关的函数可以来自理论公式也可以来自全波仿真软件导出的数据。最省事的做法是用理论公式。例如半波偶极子在某个平面内的方向图因子近似为cos(π/2cosθ)/sinθ而微带贴片则可以用腔模理论的近似公式。两者乘上同一个阵列因子后主瓣宽度、副瓣起伏都会有明显不同。所以你在报告里写明“本次仿真采用全向阵元模型”和“采用贴片阵元模型”结论的适用范围是完全不同的。在代码结构上我一般会把方位角φ也纳入计算因为很多阵列天线仿真需要对整个上半空间画三维方向图。经典做法是用meshgrid生成θ和φ的二维网格然后分别计算每个角度对应的阵元因子与阵列因子。但初学者建议先做一维切片也就是固定φ0°或φ90°只扫描θ。这样能更直观地看到主瓣、副瓣和栅瓣的变化规律也方便和解析公式对照。等一维切片验证通过再把二维网格套进来。3. 用MATLAB算均匀直线阵方向图从一行行代码到能用的函数3.1 最小可复现的主程序参数、循环与极坐标绘图我先把一个能直接跑通的均匀直线阵方向图仿真脚本拆给你看。这个脚本覆盖了最核心的流程设置阵列参数、生成空间角度网格、计算阵列因子、取模归一化、画极坐标图。我习惯把参数全部放在脚本开头方便后续改数。% 均匀直线阵方向图仿真等幅激励 clear; clc; close all; % 基本参数设置 N 8; % 阵元数量 d 0.5; % 阵元间距单位波长lambda theta0 0; % 波束指向角度度0为侧射阵 % 角度扫描范围 theta linspace(-pi/2, pi/2, 1801); % -90°到90°步长0.1° % 波数向量k*d*sin(theta) beta 2 * pi * d * sin(theta); % 空间相位差未含阵元初始相位 % 若波束扫描阵元激励相位为 -k*d*sin(theta0)写为复数权值 % 这里等幅激励各阵元相位均为0权值向量全1 w ones(N, 1); % 计算阵列因子 AF zeros(size(theta)); for n 0:N-1 AF AF w(n1) * exp(1j * n * beta); end % 归一化方向图用最大值归一化取20log10转dB AF_norm abs(AF) / max(abs(AF)); AF_dB 20 * log10(AF_norm eps); % 加eps避免log10(0) % 绘图半空间极坐标 figure; polarplot(theta, AF_norm, LineWidth, 1.5); title(sprintf(均匀直线阵方向图 N%d d/lambda%.2f, N, d));这段代码的逻辑很简洁每个阵元的相位贡献分成n倍的空间相位差权值w控制幅度和初始相位。等幅激励时w全为1所以AF就是等比数列求和。角度网格取了1801个点保证曲线平滑。极坐标图里横坐标是θ纵坐标是归一化后的幅度。注意这里你看到的副瓣高度、零点位置和公式推导完全一致可以直接拿来做解析解对照。参数变更时的行为值得仔细说。阵元间距d用波长作单位这是方向图仿真最核心的尺度参数。d0.5表示间距为半个波长此时-90°到90°范围内不会出现栅瓣是多数阵列设计的基础工况。如果把d改为1.0你会在±90°附近看到幅度与主瓣相等的额外波瓣这在物理上意味着天线会在多个方向同时产生强辐射工程上通常要避免。波束指向theta0参数在等幅激励下先不改后面做波束扫描时它才会真正起作用。3.2 副瓣电平、半功率波束宽度和零点的精确提取画出一张漂亮的极坐标图只是第一步。在实际工程汇报里你需要从方向图数据里提取三个关键指标半功率波束宽度HPBW、第一副瓣电平SLL、主瓣零点位置。用MATLAB做这件事时我倾向于不依赖图形界面手动读数据而是写一个简单的数值提取流程这样参数一改指标自动更新。% 从方向图数据提取关键指标 % 输入theta弧度或度AF_norm归一化幅度0~1 % 输出HPBW、第一副瓣电平、零点位置 % 1. 找到主瓣峰值位置 [maxVal, maxIdx] max(AF_norm); mainlobe_dir theta(maxIdx) * 180 / pi; % 2. 半功率点幅度下降到0.7071即-3dB处的角度范围 half_power 0.7071; cross_idx find(AF_norm half_power); HPBW (theta(max(cross_idx)) - theta(min(cross_idx))) * 180 / pi; % 3. 第一副瓣在主瓣两侧各找局部最大值 % 先把主瓣区域去掉取主瓣两侧区间 left_region AF_norm(1:maxIdx-10); right_region AF_norm(maxIdx10:end); % 用findpeaks找局部极大值取绝对值最大的那个 [pks_left, locs_left] findpeaks(left_region); [pks_right, locs_right] findpeaks(right_region); if ~isempty(pks_left) leftSLL_dB 20 * log10(max(pks_left)); else leftSLL_dB -inf; end if ~isempty(pks_right) rightSLL_dB 20 * log10(max(pks_right)); else rightSLL_dB -inf; end fprintf(主瓣方向: %.2f°, HPBW: %.2f°, 左第一副瓣: %.2f dB, 右第一副瓣: %.2f dB\n, ... mainlobe_dir, HPBW, leftSLL_dB, rightSLL_dB);这段提取代码有一点需要注意findpeaks会找出所有局部极大值包括方向图两边靠近±90°的边界隆起。如果阵列因子在扫描范围内没有出现完整的第一副瓣findpeaks的返回值可能包含边界伪峰判断时要看返回的幅值是否显著低于主瓣。另一方面均匀直线阵的第一副瓣理论值是-13.46dB如果你算出来是-13.0dB左右说明网格步长不够密峰值采样偏低把角度步长从0.1°加密到0.01°能改善这个误差。这是很典型的采样精度问题不是代码逻辑错误。在半功率波束宽度的计算里我用的是“幅度大于0.7071的角度范围”这个判据。这个做法对单主瓣方向图精确有效但如果你把波束扫描到大角度方向图可能只在某一侧有完整主瓣另一侧被扫描盲区截断这时候HPBW会偏大甚至无法计算。遇到这种情况我通常会改用局部寻根的方式在主瓣峰值两侧分别找幅度等于0.7071的角度点再取差值这样即使两侧不对称也能得到稳定的结果。4. 波束扫描与幅度加权从均匀阵到低副瓣设计4.1 相位加权实现波束扫描为什么大角度扫描会变形均匀直线阵的波束指向可以通过改变相邻阵元间的激励相位差来控制。你只需要给第n个阵元乘上一个相位因子使各阵元在目标方向θ0上同相叠加。把前面的脚本稍微改一下权值从全1变成含相位差的复数向量就能看到主瓣跟着θ0移动的效果。% 波束扫描权值叠加相位梯度 N 16; d 0.5; theta0 30; % 期望波束指向度 % 角度网格 theta linspace(-pi/2, pi/2, 1801); beta 2 * pi * d * sin(theta); % 相邻阵元相位差-k*d*sin(theta0) phase_shift -2 * pi * d * sind(theta0); % 权值向量 n (0:N-1); w exp(1j * n * phase_shift); % 阵列因子 AF zeros(size(theta)); for idx 1:length(beta) AF(idx) w. * exp(1j * n * beta(idx)); % 向量化内积 end AF_norm abs(AF) / max(abs(AF)); figure; polarplot(theta, AF_norm, LineWidth, 1.5); title(sprintf(波束扫描 theta0%d°, N%d, d%.1fλ, theta0, N, d));这段代码把原来的循环累加改写成了向量内积这是一个值得养成的习惯权值向量和“相位响应向量”做内积既快又不容易下标出错。这里的关键参数是phase_shift它决定了波束指向。请记住这个公式相邻阵元相位差等于-2π/λ·d·sinθ0。你用度数制sind函数时千万别把0.5波长直接乘进sind单位混淆是这类仿真最常见的低级错误。大角度扫描的物理现象在这段代码里体现得很清楚。当θ0从0°增加到60°时主瓣会变宽峰值增益下降而另一侧的副瓣会抬升。这是阵列口径投影效应造成的波束偏离侧射方向后阵列在视线方向上的有效口径变小方向图被展宽。工程经验是扫描角超过60°后均匀直线阵的方向图质量明显恶化所以很多相控阵设计把扫描范围限制在±60°以内。你仿真时会发现θ060°的主瓣宽度几乎是0°时的两倍这不是脚本写错了而是物理规律本身如此。4.2 切比雪夫与泰勒加权低副瓣的代价是主瓣变宽均匀激励的副瓣电平约-13.5dB很多实际系统要求副瓣低于-30dB那就需要引入幅度加权。切比雪夫加权在给定副瓣电平下能达到最窄的主瓣宽度泰勒加权则是它的近似实现避免切比雪夫加权在阵元数较多时出现两端阵元激励幅度跳变过大的问题。在MATLAB里做这两种加权最常见的路径是调用通信工具箱的函数例如chebwin。如果你没有工具箱自己写切比雪夫加权也不是不行但需要用到切比雪夫多项式的递推关系代码量会大不少。% 切比雪夫加权指定副瓣电平 N 16; side_lobe_db -30; % 希望达到的副瓣电平 w_cheb chebwin(N, abs(side_lobe_db)); % 注意参数是正数 % 泰勒加权需要自行生成省略通常用于差方向图或更大阵列 % 把w_cheb代入前面方向图计算流程比较即可我用这个函数时踩过一次参数坑chebwin第二参数必须是非负的dB值如果你按习惯传-30MATLAB直接报错。此外加权后方向图的主瓣会比均匀阵略宽副瓣被压下去以后换来的代价就是角度分辨率的损失。做仿真对比时一定要在报告里同时画出均匀阵和切比雪夫阵的方向图并标出两者的HPBW和SLL否则别人看不出你有没有理解这个权衡。泰勒加权在MATLAB里没有直接的内置函数常见做法是自己写一个基于泰勒线源分布的近似生成函数。它的思路是用一组“准副瓣”替代切比雪夫的等副瓣条件让靠近主瓣的副瓣逐渐衰减这样阵元激励幅度分布更平滑。实际仿真结果显示泰勒加权的副瓣包络不是完全平坦的而是随角度增大缓慢下降这对天线加工和馈电网络设计比较友好。你只需要记住结论切比雪夫适合小阵列、严格副瓣指标泰勒适合大阵列、工程可实现性优先。5. 平面阵与方向图乘积定理的推广把一维结论延伸到二维5.1 矩形平面阵的阵列因子分解技巧行因子与列因子的相乘当你从线阵转向平面阵时方向图计算的核心技巧是分解。一个M行N列的矩形平面阵如果阵元按矩形栅格排列且行间与列间的激励相位可以分离那么整个阵列因子等于行方向阵列因子和列方向阵列因子的乘积。这个乘积关系的成立条件比较严格阵元在x方向的相位梯度和y方向的相位梯度相互独立且幅度分布可以写成两个一维分布的乘积。在MATLAB中平面阵的典型做法是用二维u-v空间进行计算。定义usinθcosφvsinθsinφ然后让阵列因子成为u和v的二维函数。这样做的好处是方向图在u-v平面上呈周期结构栅瓣条件和一维情况可以直接类比。你只需要分别计算每个方向余弦对应的相位项对全部阵元求和即可。% 矩形平面阵阵列因子矩形栅格 M 8; % x方向阵元数 N 8; % y方向阵元数 dx 0.5; % x方向间距波长 dy 0.5; % y方向间距波长 % 方位角与俯仰角网格 az linspace(-pi/2, pi/2, 361); el linspace(-pi/2, pi/2, 181); [AZ, EL] meshgrid(az, el); % 方向余弦 u sin(EL) .* cos(AZ); v sin(EL) .* sin(AZ); % 阵元位置下标 m (0:M-1); % x方向 n (0:N-1); % y方向 % 阵列因子等幅激励 AF zeros(size(u)); for p 1:numel(u) AF(p) sum(exp(1j*2*pi*(m*dx*u(p) n*dy*v(p))), all); end AF_norm abs(AF) / max(abs(AF)); % 三维方向图显示用surf或mesh figure; surf(u, v, AF_norm); xlabel(u); ylabel(v); zlabel(归一化幅度);这个脚本最需要注意的地方是变量维度。m是列向量n也是列向量在求和时要确保u(p)和v(p)是标量然后利用MATLAB的隐式扩展把mdxu(p)和ndyv(p)叠加成M×N的矩阵。这里如果忘记转置维度对不上要么报错要么得到错误的结果还很难发现。平面阵的方向图在u-v空间看是最直观的主瓣是一个锥形峰栅瓣出现在u方向或v方向的整数倍周期上它们的位置由dx和dy决定。所以平面阵的栅瓣条件是dx和dy都不超过0.5λ和一维线阵的结论一致但需要两个方向同时检查。5.2 栅瓣判据d/λ的极限值别只看公式栅瓣判断是阵列方向图仿真里最容易被忽略的环节。一维均匀直线阵的栅瓣条件由d/λ决定当间距超过一个波长时可见空间内会出现第二个峰值。这个结论的物理本质是空间相位差的周期性与可见范围的有限性之间的矛盾。很多同学把公式背下来但遇到非均匀阵列、波束扫描或者平面阵时就不知道如何用公式判断了。我的排查习惯是先画一张“全空间阵列因子”图也就是不裁剪可见空间直接看阵列因子的周期性栅瓣在哪里。用MATLAB把β从-4π扫到4π画出的阵列因子图像显示出一串等高的峰值相邻峰值间隔为2π。可见空间对应β的取值范围是[-2πd/λ·sinθ, 2πd/λ·sinθ]。只要这个范围里落入第二个峰值方向图里就会出现栅瓣。这个方法比死记dλ的判断准得多特别是在波束扫描时主瓣移到θ0原来的栅瓣位置也会移动某个原来刚好在可见范围边缘的栅瓣可能在扫描后进入可见范围。按上面的判据你能推出一个重要工程结论如果要求波束扫描到±60°阵元间距必须满足d/λ小于约0.577才能保证全扫描范围内不出现栅瓣。这就是很多相控阵实际采用0.5λ~0.55λ间距的原因不是随便定的。我参与过的模拟项目X里最初为了省阵元数量用了0.6λ间距一扫描到45°就看到了栅瓣进入可见区换回0.5λ后问题消失。这个案例说明方向图仿真的价值不只是画一张好看的图而是能在设计阶段提前暴露这类硬约束。6. 阵列方向图仿真避坑指南五个反复出现的翻车点6.1 现象polarplot出现多条放射状直线方向图杂乱无章先说我自己最早翻车的一次。用polarplot画方向图时图形里出现了从圆心向外辐射的直线簇看起来像蜘蛛网完全不是平滑的方向图曲线。检查代码后发现角度数据里有NaN值。原因出在归一化这一步如果某一角度上所有阵元相位相互抵消阵列因子的模恰好为0直接除以最大值后会得到0但如果你先求了log10再归一化log10(0)会给出-inf。绘图函数遇到-inf或NaN时会把曲线画得乱七八糟。解决办法很简单两个位置要处理一是求阵列因子后对abs结果做归一化二是要画dB曲线时给幅度加上一个极小量eps比如用20*log10(AF_normeps)。另外用polarplot时如果角度向量不是严格递增且无缝覆盖整个区间也可能出现跨零点的连线。建议先把角度排序再绘图避免画出一根横穿整个图形的直线。6.2 现象阵元间距大于一个波长方向图“多了一个主瓣”但报告里没解释间距超过一个波长时方向图会出现栅瓣这是理论课上反复强调的内容但仿真时很多人看到多出来一个等高的“主瓣”就直接认为是代码写错了。实际上这是物理规律说明你的阵列已经出现了栅瓣。这个现象的出现条件是d/λ1而且栅瓣的位置可以通过sinθ_g mλ/d来估算m±1时对应第一对栅瓣。正确的处理方式不是删掉栅瓣而是在报告中明确计算并标注栅瓣位置同时说明实际天线系统会通过限制扫描范围或减小阵元间距来避免栅瓣落入感兴趣的区域。如果你在项目里发现“主瓣数量变多”先确认d/λ的值别急着改代码。有些同学为了消除栅瓣把间距调得特别小结果方向图主瓣变得非常宽阵列的增益和分辨率都下降了这就是把栅瓣问题和口径尺寸问题混为一谈。6.3 现象加窗后副瓣没有降下来反而出现不对称的杂散峰这个坑我印象最深。用chebwin给16元阵列加-30dB切比雪夫权值后方向图的左副瓣降到了-30dB左右右副瓣却还有-20dB以上的“杂散峰”怎么看都像代码出了bug。排查了很久才发现问题不在于权值函数而在于角度网格不均匀导致峰值采样偏差。方向图的副瓣位置与角度网格点不完全重合时findpeaks找到的局部最大值比真实副瓣峰值小所以左右副瓣读数不一致。解决方法是加密角度扫描网格把步长从1°改成0.01°让每个副瓣区域有足够的采样点。另一个更稳健的办法是先用插值函数对方向图做样条插值再找峰值。我看到有人用sinc函数插值恢复方向图的连续值效果也不错只是代码量略大。6.4 现象波束扫描方向图在某个角度上主瓣突然消失在仿真波束扫描从0°扫到80°时主瓣在某个角度后幅度骤降有时直接消失在噪声底之下。这通常不是阵列因子算错了而是你选的阵元因子模型在扫描角附近就是没有辐射能力。最典型的例子是阵元因子为cosθ的模型它模拟的是一种理想偶极子的方向性在θ接近±90°方向增益自然趋近于零。总方向图等于阵元因子乘以阵列因子所以在可见区边缘总方向图必然掉下去。此时你需要区分“阵列因子还在增大”和“总方向图在衰减”两个概念报告里如果只画总方向图会让人误以为扫描能力不足。我在实际项目里就遇到过一个类似的困惑某次仿真显示波束扫描到70°时总方向图几乎消失但阵列因子明明还有明显峰值。后来在报告里分开画出阵元因子、阵列因子、总方向图三条曲线问题一目了然。所以仿真时一定要保留中间变量不要只出最终结果。6.5 现象同样参数两次运行结果不一样怀疑随机性方向图本身是确定性计算不该有随机性。但如果你在代码里用了randn生成权值或者误用了未固定的随机数种子结果自然每次不同。更隐蔽的情况出现在使用某些工具箱函数时例如自适应阵列的方向图计算会涉及采样协方差矩阵估计而这个估计通常基于随机快拍数据。排查思路是第一步检查是否用了随机函数第二步检查代码中是否有排序函数对相同数据给出不同顺序第三步检查是否依赖全局变量或并行计算时的工作区分裂。对教学仿真来说建议在脚本开头固定随机种子例如rng(2024)这样别人复现你的结果时数据完全一致。固定随机种子是我在提交仿真作业前必做的一步它不算复杂但能避免掉很多“玄学复现”的争议。7. 方向图仿真的验证闭环独立检验你的结果是否可信最后一章我想给你一个可以立即操作验证流程这也是我自己做方向图仿真时最看重的一个习惯。仿真可以说服自己但要让别人信服你必须提供两条独立路径的比对。我一般会做三件事第一用解析公式手算几个关键角度的方向图数值和MATLAB输出对比第二用全波仿真软件或另一个独立计算方法对同一阵列建模比对主瓣位置和宽度第三把阵元数N和间距d取极限看方向图是否退化为已知的特例。解析对比的做法很轻量以均匀直线阵为例在θ0°时所有阵元同相叠加阵列因子应为N。你把N8代入计算abs(AF)的最大值应该恰好等于8归一化后才是1。如果最大模值不是整数或者明显偏离N说明相位项的计算有误。第二个特例是N1这时阵列因子应该恒等于1方向图是各向同性的圆。任何方向图仿真只要N1时画出方向图不是圆形相位公式肯定有问题。第三个常见的验证方法是把方向图积分看辐射功率是否守恒。虽然方向图仿真不涉及绝对功率但归一化方向图的积分在理想情况下应当与天线辐射的总功率成正比。微带阵列仿真里我经常用方向图积分估算阵列的定向性directivity如果算出的定向性和理论值例如主瓣立体角估算差得超过1dB就要怀疑是不是角度采样太粗或者网格不均匀导致的积分误差。在实际汇报场合中评审人最容易接受的数据呈现方式是一张对比表列出均匀阵、切比雪夫加权阵和扫描状态下的主瓣方向、HPBW、第一副瓣电平、栅瓣位置。这样一张表比十张方向图都更有说服力。我自己在模拟项目X里每次调整参数后都会自动生成这张表然后截图放进实验记录里。多年的习惯养下来我发现写报告时最缺的从来不是“图”而是“图对应的关键参数表”。方向图仿真这个方向入门容易精通难。很多人以为学会了画极坐标图就等于会了阵列天线其实真正的功夫在那些看不见的细节里栅瓣判据、阵元因子的选取、加权方式的权衡、采样密度对峰值提取的影响。我在这些细节上踩过的坑远不止上面五条但每一条都在后来的仿真中变成了下意识的检查项。仿真不是把公式翻译成代码而是用代码验证你对物理的假设是否正确。每次跑完一组结果多问一句“这个图物理上可能的边界条件是什么”你的仿真能力会提升得很明显。希望这篇沿着阵列天线方向图及其MATLAB仿真整理的实战笔记能帮到你也希望你在对照复现之后能画出经得起追问的方向图。本文还有配套的精品资源点击获取