简介这份资源围绕阵列天线方向图展开面向无线通信、天线设计与电磁场相关课程的学习者及工程入门人员帮助理解单元个数、阵元间距与波长变化对辐射方向图的影响。压缩包共4个文件均为MATLAB脚本.m整体约2KB分别对应不同参数组合下的方向图绘制与对比便于直接运行观察波束宽度、增益与旁瓣变化。资源已有2136人学习下载说明其在阵列天线基础教学中具有一定参考价值。读者可通过修改脚本中的单元数、间距与波长参数复现2单元、4单元、8单元及λ/2、λ、λ/4等典型配置的方向图直观体会增益随单元数增长、间距对旁瓣抑制的作用以及波长与间距比值对主瓣和零点的影响从而为波束赋形、远程通信与空间分辨率优化等场景积累设计直觉。1. 从一张对不齐的波束图说起阵列方向图比较要解决什么做阵列天线方向图比较最怕的不是算不出来而是算出来之后对不上。某同学把均匀线阵、切比雪夫加权、泰勒加权三组方向图叠在一张图上旁瓣电平看着差不多主瓣宽度却差了一截回头查代码发现三组仿真的阵元间距一个用半波长、一个用 0.45 波长横轴角度分辨率也不一致。这类翻车在课程设计和毕设里非常常见根源不在公式而在于没有把“比较”当成一个受控实验来设计。这套方案要解决的就是给定阵元数、间距、加权方式、扫描角这几个输入输出一组可横向对照的方向图曲线并给出主瓣宽度、第一旁瓣电平、零点位置等可量化指标。适合正在做阵列天线课设、毕设或者需要给一份天线选型对比材料的人把“画几张图”升级成“做一次有结论的对比”。2. 需求拆解与参数约定先定比较维度再谈画图方向图比较最容易犯的错是先写仿真再想清楚要比什么。阵列天线的方向图由阵因子和阵元因子共同决定比较时通常固定阵元因子、只比阵因子否则变量太多结论没法归因。所以第一步是把比较维度锁死比的是不同加权方式下的旁瓣抑制能力还是不同阵元数下的波束锐度还是不同间距下的栅瓣出现条件。这三个问题的参数敏感度完全不同混在一起比就是自找麻烦。2.1 把输入参数收敛成一张表我一般会先列一张参数约定表把所有仿真共用的量固定下来只留一个自变量。下面这张表是这套方案里默认的基准配置读者可以按自己的题目改数值但结构建议保留。参数符号基准值说明阵元数N16均匀线阵比较加权时固定阵元间距d0.5λ避免栅瓣的常用取值工作波长λ1归一化角度计算与波长解耦扫描角θ00°比较加权时保持侧射角度扫描范围θ-90°~90°步长 0.1°加权方式w均匀/切比雪夫/泰勒比较的主变量这张表看着朴素但它决定了后面所有曲线能不能叠在一张图上。间距用 0.5λ 是血泪经验一旦超过 0.5λ可见空间里会冒出栅瓣旁瓣电平的比较就失去意义因为最大的那个“旁瓣”其实是栅瓣。扫描角固定为 0° 也是同理扫描会让波束展宽加权比较的结论会被扫描效应污染。2.2 明确输出指标而不只是输出图只给三张方向图答辩时很容易被问“所以呢”。合格的做法是每张图配一组量化指标主瓣宽度通常取 -3dB 波束宽度、第一旁瓣电平dB、零点个数与位置、以及峰值旁瓣电平。这些指标从方向图数据里直接算不靠肉眼读图。比较的结论就落在这些数字上比如“切比雪夫加权把第一旁瓣压到 -30dB代价是主瓣展宽约 1.3 倍”这句话比三张叠图有说服力得多。提示指标计算和画图要用同一份角度采样数据不要画图用一套网格、算指标再插值一套否则主瓣宽度会算出零点几度的偏差评审一眼就能看出来。3. 方向图计算的核心实现阵因子、加权与归一化这一章落到代码。阵列方向图比较的计算量其实很小难点全在细节加权系数怎么生成、方向图要不要归一化、dB 转换时零深怎么处理。下面用 Python 给出一套可直接跑的骨架NumPy 负责向量化Matplotlib 负责出图不依赖任何天线专用库方便读者逐行对照公式。3.1 阵因子与三种加权系数的生成均匀线阵的阵因子在侧射时可写成各阵元相位的累加扫描角引入后加一个线性相位项。加权系数决定了各阵元的幅度分布均匀加权全 1切比雪夫和泰勒加权需要按阶数生成。切比雪夫加权可以用多项式递推实现泰勒加权用理想空间因子近似课程设计里通常给到现成公式这里给出可运行的实现。import numpy as np def array_factor(N, d_lambda, theta_deg, weights, theta0_deg0.0): 计算线阵归一化阵因子dB。 N: 阵元数; d_lambda: 间距/波长; theta_deg: 角度采样(度) weights: 长度 N 的幅度加权; theta0_deg: 扫描角 theta np.deg2rad(theta_deg) theta0 np.deg2rad(theta0_deg) n np.arange(N) # 相邻阵元相位差含扫描相位 psi 2 * np.pi * d_lambda * (np.sin(theta) - np.sin(theta0)) # 各阵元贡献按幅度加权后累加 af np.sum(weights[:, None] * np.exp(1j * n[:, None] * psi[None, :]), axis0) af np.abs(af) af af / np.max(af) # 归一化到主瓣峰值 return 20 * np.log10(np.maximum(af, 1e-6)) # 防止 log(0) def uniform_weights(N): return np.ones(N) def chebyshev_weights(N, sll_db30): 切比雪夫加权sll_db 为设计旁瓣电平(dB)。 R 10 ** (sll_db / 20) x0 np.cosh(np.arccosh(R) / (N - 1)) k np.arange(N) # 用切比雪夫多项式在等效角度上取值作为幅度 u x0 * np.cos(np.pi * k / (N - 1)) w np.cosh((N - 1) * np.arccosh(np.clip(np.abs(u), 1, None))) return w / np.max(w)array_factor里最关键的是psi这一项它把角度映射成相邻阵元的相位差扫描角通过减去sin(theta0)实现波束偏转。np.maximum(af, 1e-6)是防 log 零的后悔药方向图零点处幅度趋近于零直接取对数会得到负无穷画图时曲线会断掉。切比雪夫那段用等效角度上的切比雪夫多项式取值sll_db控制设计旁瓣改这个参数就能得到不同抑制水平的加权比较时建议取 20、30、40 三档曲线差异会很明显。3.2 归一化与 dB 显示的坑方向图比较里最常见的显示问题是归一化基准不一致。有的曲线按各自峰值归一化有的按全局峰值归一化叠在一起主瓣高度看着一样旁瓣的相对关系却错了。正确做法是比较同一组阵元数和间距下的加权差异时按各自峰值归一化是合理的因为关心的是旁瓣相对主瓣的抑制但比较不同阵元数时必须按同一基准归一化否则波束锐度的差异会被归一化吃掉。这个区别在代码里就是一个除法位置的问题但答辩时被问到“为什么两条曲线主瓣一样高”就露馅了。theta np.linspace(-90, 90, 1801) N, d_lambda 16, 0.5 configs { uniform: uniform_weights(N), chebyshev_30dB: chebyshev_weights(N, 30), chebyshev_40dB: chebyshev_weights(N, 40), } for name, w in configs.items(): af_db array_factor(N, d_lambda, theta, w) # 后续画图与指标计算共用 af_db 和 theta这段循环把三组加权跑在同一套角度网格上theta步长 0.1° 对应 1801 个点主瓣宽度算到 0.1° 精度足够课设使用。如果题目要求更细把 1801 改成 3601 即可计算量对 NumPy 来说可以忽略。注意d_lambda传的是间距与波长之比不是绝对长度这样整套代码与频率解耦换频段不用改任何数值。3.3 从曲线到指标主瓣宽度与旁瓣电平的自动提取画完图不算完指标要自动算。主瓣宽度找峰值两侧第一次跌到 -3dB 的角度差第一旁瓣电平找主瓣之外第一个局部极大值。用scipy.signal.find_peaks可以省事但课设环境不一定装了 SciPy手写一个局部极值扫描也不复杂。def beamwidth_3db(theta, af_db): peak_idx np.argmax(af_db) half af_db[peak_idx] - 3.0 left peak_idx while left 0 and af_db[left] half: left - 1 right peak_idx while right len(af_db) - 1 and af_db[right] half: right 1 return theta[right] - theta[left] def first_sidelobe_db(theta, af_db): peak_idx np.argmax(af_db) # 主瓣右侧第一个局部极大 for i in range(peak_idx 1, len(af_db) - 1): if af_db[i] af_db[i-1] and af_db[i] af_db[i1]: return af_db[i] return Nonebeamwidth_3db从峰值向两侧线性扫描遇到低于半功率点就停返回角度差。这个写法对单主瓣方向图够用如果扫描角偏离侧射导致主瓣不对称左右半功率点要分别算再取平均代码里加两行即可。first_sidelobe_db只扫主瓣右侧是因为均匀加权方向图左右对称但加权后仍对称取一侧足够若做了非对称加权两侧都要扫并取较大值。这两个函数返回的数字直接填进比较表比截图读图靠谱得多。4. 比较实验的组织三组对照与结果呈现代码能跑之后真正的设计工作是把比较实验组织清楚。我的习惯是每组对照只改一个变量其余全部锁死并且每组都输出“曲线 指标表 一句结论”。这样即使评审临时要求换参数重跑也能快速定位是哪一组受影响。4.1 三组对照的变量设计第一组比加权方式N16、d0.5λ、扫描角 0°均匀、切比雪夫 30dB、切比雪夫 40dB 三条曲线。预期结论是旁瓣越低主瓣越宽用指标表量化这个 trade-off。第二组比阵元数加权固定为均匀N 取 8、16、32看波束锐度随阵元数的变化主瓣宽度大致与 N 成反比这条规律在答辩时能直接引用。第三组比阵元间距N16、均匀加权d 取 0.4λ、0.5λ、0.6λ重点看 0.6λ 时栅瓣是否进入可见空间这一组是解释“为什么间距不能随便加大”的最好素材。对照组自变量固定量观察指标组一加权方式N16, d0.5λ, θ00°第一旁瓣电平、-3dB 宽度组二阵元数 N均匀加权, d0.5λ主瓣宽度、零点个数组三间距 dN16, 均匀加权栅瓣出现角度、峰值旁瓣这张表建议直接放进设计文档答辩时按组讲逻辑非常顺。每组的三条曲线用同一套坐标轴范围纵轴统一 -60dB 到 0dB横轴统一 -90° 到 90°避免视觉误导。曲线颜色和线型要固定比如均匀加权永远用实线、切比雪夫用虚线读者不看图例也能分辨。4.2 结果呈现的排版细节方向图比较的图最容易犯的毛病是信息过载。一张图里塞五六条曲线、图例挡住主瓣、网格线比曲线还密评审根本看不清。我的做法是一张图最多三条曲线主瓣区域可以加一个局部放大子图把 -3dB 宽度和第一旁瓣的细节放大出来。指标表放在图下方列名与正文一致不要图里叫“旁瓣”、表里叫“SLL”、正文又叫“副瓣电平”术语统一是专业度的直接体现。注意dB 纵轴的下限不要设成 -100dB 以下数值计算的底噪会让曲线在零点附近出现锯齿看起来像“旁瓣”其实是数值噪声。统一截到 -60dB既干净又不影响结论。5. 避坑与排查方向图比较里最容易翻车的五件事这一章是我带课设时踩过的坑合集每条都按“现象、原因、解决”写读者对照自己的代码逐条排查基本能省掉大半调试时间。5.1 曲线叠在一起但旁瓣对不上现象是两组加权曲线主瓣几乎重合旁瓣电平却差了好几个 dB和理论值对不上。原因通常是加权系数没有归一化到峰值 1或者归一化之后又乘了一个常数导致阵因子峰值不在 0dB。解决方法是检查weights的最大值是否为 1并在array_factor里统一做af / np.max(af)不要在两处重复归一化。5.2 间距稍大就冒出莫名其妙的“旁瓣”现象是 d 取 0.6λ 时±60° 附近出现一个和主瓣差不多高的峰误以为是加权出了问题。原因是间距超过半波长后栅瓣进入可见空间这是物理规律不是 bug。解决方法是先确认d_lambda是否大于 0.5若是在结论里明确写出栅瓣角度公式为arcsin(λ/d - sinθ0)把它作为对照组三的核心发现而不是当成错误去“修”。5.3 主瓣宽度算出来是 0 或者跳变现象是beamwidth_3db返回 0 或异常大的值。原因多半是角度网格太粗半功率点落在两个采样点之间线性扫描直接跨过或者方向图峰值不在网格中心左右扫描不对称。解决方法是将步长细化到 0.05° 以内并在扫描前用np.argmax确认峰值索引合理必要时对半功率点做线性插值。5.4 dB 曲线在零点处断裂现象是画出来的曲线在零点附近断开或出现竖直跳变。原因是幅度趋近零时log10得到负无穷Matplotlib 遇到-inf会断线。解决方法是在取对数前加一个下限如np.maximum(af, 1e-6)或者直接用20*np.log10(af 1e-12)两种写法等价选一种统一用。5.5 换了阵元数结论就反了现象是 N8 时切比雪夫加权旁瓣更低N32 时反而均匀加权更好看。原因通常是不同 N 下用了不同的设计旁瓣参数或者归一化基准不一致比较条件没锁死。解决方法是把加权设计参数和归一化方式写进配置字典每组实验打印一份参数快照确保“只改一个变量”这条原则真的落实了。6. 验收与答辩用一张对照表把结论讲死方向图比较这类题目答辩时最容易被追问的不是公式推导而是“你的比较公平吗、结论站得住吗”。我的验收习惯是准备一张总对照表把三组实验的自变量、固定量、关键指标和一句话结论全部列进去评审扫一眼就能判断实验设计是否受控。表里每个指标都要能追溯到代码输出不能有手填的估计值。下面是一个填写示例数值按实际跑出来的结果替换。实验组配置第一旁瓣(dB)-3dB 宽度(°)结论加权对比均匀-13.26.4基准加权对比切比雪夫 30dB-30.18.1旁瓣降 17dB主瓣宽 1.27 倍加权对比切比雪夫 40dB-39.89.6抑制更强展宽代价继续增大阵元数N8-13.112.8宽度约与 N 成反比阵元数N32-13.33.2验证反比规律间距d0.6λ栅瓣 -0.56.4栅瓣出现在约 ±56°答辩陈述按“需求、参数约定、实现要点、对照结果、结论”五段走每段控制在两分钟内。被问到实现细节时直接打开array_factor讲psi的物理含义比背公式更可信。被问到结论局限性时主动说明本方案只比了阵因子、未考虑阵元互耦和馈电误差这反而是加分项说明你知道边界在哪。还有一个实用技巧提前把三组曲线导出成矢量图答辩 PPT 里放大不失真比截图清晰一个档次。从那以后我每次带这类题目都强制要求先交参数约定表和对照表代码可以后补但这两张表不齐就不往下做。方向图比较的功夫不在画图在于把比较条件锁死、把指标算准、把结论说清楚这三件事做到位方案自然立得住。希望帮到你。本文还有配套的精品资源点击获取
