简介这份资源面向学习雷达信号处理与阵列天线仿真的学生、工程师和研究人员聚焦相控阵波束扫描这一核心问题用Matlab把抽象的天线相位控制过程转化为可直观观察的动态演示。压缩包共2个文件均为m脚本整体约2KB体量轻巧分别承担辐射方向图计算与极坐标分贝显示等功能便于直接运行和二次修改。资源已有3434人学习下载说明其在相控阵入门与教学演示中具有一定认可度。读者可借助脚本理解天线阵列几何建模、单元相位设定、合成波束方向与强度计算等关键环节并通过图形更新观察波束随相位变化扫过空间的过程进而分析波束宽度、旁瓣电平等性能指标。代码结构简洁适合作为课堂演示、课程设计或自学相控阵原理的起点也便于在此基础上扩展阵列形式与扫描策略。1. 从一张“会动的波束图”说起相控阵扫描仿真到底要交付什么做相控阵课题的人大多遇到过同一个尴尬公式推了一黑板阵列流形、导向矢量、阵因子都写出来了可一到答辩现场评委问“波束到底怎么扫的”你只能翻出一张静态方向图指着主瓣位置说“这里会动”。静态图讲不清扫描这是这个项目最直接的痛点。这套 Matlab 仿真要交付的东西很明确一个能看见主瓣随移相量连续偏转的波束扫描动图配套可复现的阵列建模、导向矢量计算、方向图合成与动画导出脚本。它适合正在做相控阵课设、毕设的学生也适合需要快速验证阵元数、阵元间距、扫描角对波束形状影响的技术人员。核心词就三个相控阵、Matlab、波束扫描后面所有设计都围绕这三件事展开。2. 需求拆解与技术选型为什么是均匀线阵加相移扫描2.1 先把“扫描”这件事拆成可计算的量波束扫描的本质不是天线在转而是每个阵元上的激励相位在变。均匀线阵ULA沿一条直线排布 N 个阵元间距 d来波方向与阵列法线夹角为 θ相邻阵元之间的波程差是 d·sinθ对应相位差就是 2πd·sinθ/λ。要让主瓣指向 θ₀就给第 n 个阵元补偿相位 −n·2πd·sinθ₀/λ所有阵元叠加后θ₀ 方向同相相加其他方向相互抵消主瓣就“指”过去了。扫描动图要展示的就是 θ₀ 从负角扫到正角时这一组补偿相位连续变化、主瓣跟着平移的全过程。这里有个容易被忽略的细节扫描角范围不是随便取的。均匀线阵在端射方向±90°附近主瓣会严重展宽甚至畸变实际仿真通常把扫描范围限制在 ±60° 以内动图看起来才“干净”。我一般会在脚本里把扫描角设成 −60° 到 60°、步进 2°一共 61 帧帧率 10 到 15 帧每秒导出的 GIF 体积可控播放也顺。2.2 技术选型Matlab 的三个理由选 Matlab 不是惯性是这个场景下它确实省事。第一矩阵运算天然贴合阵列信号模型导向矢量、阵列流形、方向图合成基本都是一两行矩阵乘法换成别的语言要写不少循环。第二绘图和动画能力开箱即用polarplot、plot 加 drawnow 就能出动态效果导出 GIF 也有现成接口。第三课程和毕设环境普遍装了 Matlab交付脚本别人拿到就能跑不用配一堆依赖。具体到函数层面方向图合成用steervec或手写导向矢量都行我倾向手写因为课设要讲清每一步手写更能体现推导过程。动画用drawnow逐帧刷新或者用getframe抓帧后写 GIF。参数化用脚本顶部的变量区统一管理N、d/λ、扫描范围、步进、帧率全部集中改一个数就能重跑这是交付质量的分水岭。提示d/λ 建议取 0.5这是避免栅瓣的经典取值。取大了会在可见空间出现栅瓣动图里会冒出多余的主瓣答辩时很难解释。2.3 交付物清单与验收口径这套仿真最终要交出的不是一张图而是一组能自洽运行的文件主脚本参数区加扫描循环加动画、方向图计算函数、GIF 导出函数外加一份说明各参数含义的注释。验收口径我通常定三条主瓣指向与设定扫描角的误差在 1° 以内扫描过程中主瓣峰值起伏不超过 1 dB导出的 GIF 能完整播放一个扫描周期且帧数正确。这三条能过基本就说明模型和动画都没跑偏。3. 总体设计阵列模型、方向图合成与动画主循环的三层结构3.1 三层结构划分整个仿真按“模型层—计算层—表现层”三层来搭职责清晰改起来不牵连。模型层负责阵元位置、波长、间距这些静态参数计算层负责给定扫描角时算导向矢量、算阵列因子、归一化表现层负责把每一帧的方向图画出来并串成动画。三层之间只通过参数结构体传递不共享全局变量这是避免脚本越写越乱的关键。用一张表把三层的关键内容列清楚层级核心内容输入输出模型层阵元数 N、间距 d、波长 λ、角度网格设计参数阵元坐标、角度向量计算层导向矢量、阵因子、归一化方向图扫描角 θ₀归一化幅度序列表现层极坐标/直角坐标绘图、帧刷新、GIF 写入每帧幅度动图文件这个划分的好处是想换圆阵或面阵只动模型层想加窗函数压低旁瓣只动计算层想换绘图风格只动表现层。课设里经常被要求“再加一个对比”分层之后加对比就是复制一份计算层调用不用重写整个脚本。3.2 角度网格与分辨率设定方向图是连续函数仿真里要离散成角度网格。网格太粗主瓣峰值位置会偏动图看起来一跳一跳网格太细每帧计算量上去动画会卡。经验值是 0.1° 一格覆盖 −90° 到 90°共 1801 个点单帧计算在毫秒级动画流畅峰值定位也够准。如果只关心 ±60° 扫描范围网格可以缩到 −90° 到 90° 保留完整可见空间因为旁瓣信息对判断栅瓣很重要不建议裁掉。3.3 主循环的节奏控制扫描动图的主循环结构是外层遍历扫描角内层算方向图并绘图每帧之间用drawnow或pause控制节奏。drawnow是立即刷新适合实时看导出 GIF 时改用getframe抓当前 figure 再写入节奏由帧率参数控制。这里有个血泪经验直接在主循环里pause(0.1)再drawnow实时看没问题但导出 GIF 时帧间隔不受控导出的动图会忽快忽慢。正确做法是实时预览和导出分成两个模式用开关变量切换导出时按固定帧率写帧。4. 核心模块实现导向矢量、方向图合成与动图导出4.1 导向矢量与阵因子的代码实现先看计算层最核心的一段给定阵元数、间距波长比、扫描角算出归一化方向图function AF array_factor(N, d_lambda, theta0_deg, theta_grid_deg) % 计算均匀线阵归一化阵因子 % N : 阵元数 % d_lambda : 阵元间距与波长之比 d/λ % theta0_deg : 波束指向角度 % theta_grid_deg : 角度网格度用于绘制方向图 theta0 deg2rad(theta0_deg); theta deg2rad(theta_grid_deg); n (0:N-1).; % 阵元索引列向量 % 相邻阵元波程差对应的空间相位 psi 2*pi*d_lambda*sin(theta); % 1×M 行向量 % 扫描补偿相位使 theta0 方向同相叠加 psi0 2*pi*d_lambda*sin(theta0); % 导向矢量与补偿叠加得到每个角度的复幅度 AF_complex sum(exp(1j*n*(psi - psi0)), 1); % 归一化主瓣峰值归一到 0 dB AF abs(AF_complex) / max(abs(AF_complex)); end这段代码的关键在psi - psi0这一项。psi是来波方向带来的阵元间相位差psi0是扫描角对应的补偿量两者相减后θ 等于 θ₀ 时每一项指数为 0N 个 1 相加得到峰值 N其他角度相位错开、矢量相消。归一化除以最大值保证不同扫描角下主瓣峰值都在 0 dB动图里主瓣高度稳定不会因为扫描角变化而忽高忽低。参数d_lambda取 0.5 时可见空间内不会出现栅瓣取 1.0 时会在 ±90° 附近出现第二个主瓣这正是课上要演示的栅瓣现象可以留作对比实验。4.2 单帧方向图绘制与坐标选择表现层单帧绘制的代码如下直角坐标和极坐标两种视图都保留function draw_pattern(theta_grid_deg, AF, theta0_deg, ax) % 在指定坐标轴绘制单帧方向图 AF_dB 20*log10(AF eps); % 转 dB加 eps 防 log(0) AF_dB max(AF_dB, -40); % 下限截到 -40 dB图形更干净 plot(ax, theta_grid_deg, AF_dB, LineWidth, 1.5); hold(ax, on); xline(ax, theta0_deg, --r, 扫描角); % 标出当前指向 hold(ax, off); xlim(ax, [-90 90]); ylim(ax, [-40 0]); xlabel(ax, 角度 (度)); ylabel(ax, 归一化幅度 (dB)); grid(ax, on); drawnow; end用 dB 而不是线性幅度是因为主瓣和旁瓣差几十 dB线性坐标下旁瓣几乎看不见动图里只能看到一根尖峰在平移说服力不足。截到 −40 dB 是绘图习惯低于这个值的旁瓣对判断波束形状没帮助留着反而让纵轴被拉长。xline标出当前扫描角观众一眼能对上主瓣位置和设定值答辩时这个细节很加分。极坐标视图用polarplot替换plot即可角度和幅度直接传入适合展示波束的空间指向感。4.3 动图导出GIF 写入的完整流程导出 GIF 的代码如下和实时预览分开function export_scan_gif(N, d_lambda, scan_angles, theta_grid, filename, fps) % 扫描并导出 GIF 动图 fig figure(Visible, off, Position, [100 100 800 500]); ax axes(fig); delay 1/fps; for k 1:numel(scan_angles) AF array_factor(N, d_lambda, scan_angles(k), theta_grid); draw_pattern(theta_grid, AF, scan_angles(k), ax); frame getframe(fig); [im, map] rgb2ind(frame2im(frame), 256); % 转索引图控体积 if k 1 imwrite(im, map, filename, gif, ... LoopCount, Inf, DelayTime, delay); else imwrite(im, map, filename, gif, ... WriteMode, append, DelayTime, delay); end end close(fig); endrgb2ind把真彩色帧转成 256 色索引图这是控制 GIF 体积的关键不转的话文件会大好几倍。LoopCount设为Inf让动图循环播放DelayTime用帧率倒数保证播放速度均匀。第一帧用普通写入后续帧用append追加这是 Matlab 写 GIF 的标准套路。Visible设为off是为了导出时不弹窗、不抢焦点批量跑多个参数组合时尤其有用。帧率我一般取 1261 帧的动图大约 5 秒播完一轮节奏刚好。注意getframe抓的是 figure 当前显示内容如果 figure 被遮挡或最小化抓到的帧可能不完整。导出时保持 figure 可见或设成 off-screen 渲染别一边导出一边切窗口。5. 避坑与排查扫描动图最容易翻车的五个点5.1 主瓣指向偏了相位符号搞反现象是设定扫描角 30°主瓣却出现在 −30°。原因几乎都是补偿相位的符号写反psi - psi0写成了psi psi0或者n的索引从 1 开始却按 0 起算。解决方法是拿 θ₀ 等于 0 做基准测试此时主瓣应在 0°再测 θ₀ 等于 30°主瓣应往正方向偏两步都对符号就没问题。我习惯在函数里加一句断言扫描角为 0 时峰值位置偏差超过 1° 就报错早发现早改。5.2 动图里冒出多余主瓣栅瓣没避开现象是扫描到较大角度时方向图里出现第二个和主瓣差不多高的峰。原因是 d/λ 大于 0.5可见空间内出现了栅瓣。解决方法是把间距降到 0.5λ 或更小或者在脚本里加一段栅瓣检查扫描角取最大时除主瓣外是否有峰值超过 −3 dB有就提示间距过大。这个检查写进验收脚本比事后看图靠谱。5.3 动画卡顿或帧率不均刷新方式用错现象是实时预览很流畅导出的 GIF 却忽快忽慢。原因是主循环里混用了pause和drawnow导出时帧间隔由计算耗时决定而非固定值。解决办法是预览和导出分离导出走getframe加固定DelayTime预览走drawnow加短pause。另外每帧重绘时用hold on/off而不是cla全清能减少闪烁。5.4 归一化后主瓣高度不稳每帧各自归一化现象是扫描过程中主瓣峰值在 0 dB 上下跳。原因是每帧都除以当帧最大值而不同扫描角下离散网格可能没正好采到峰值导致归一化基准漂移。解决办法是用解析峰值 N 做归一化基准或者把角度网格加密到 0.05°让峰值采样误差可忽略。前者更稳后者更直观课设里两种都讲一遍学生理解更深。5.5 导出文件过大或打不开编码参数没调现象是 GIF 几十兆或者某些播放器打不开。原因是没做索引色转换或者帧尺寸太大。解决办法是rgb2ind转 256 色、figure 尺寸控制在 800×500 以内、帧数控制在 60 帧上下。如果还要更小可以隔帧抽稀扫描步进从 2° 改成 4°帧数减半动图依然能看出扫描趋势。6. 验收与答辩用三张对照图讲清扫描原理验收这套仿真我一般不看代码写得多漂亮看三张图能不能自圆其说。第一张是固定扫描角 0° 的方向图标出主瓣宽度和第一旁瓣电平证明阵列模型正确第二张是扫描角 −60°、0°、60° 三条曲线叠在一张图上证明主瓣确实随设定角平移且峰值稳定第三张是动图本身证明连续扫描过程没有跳变和异常峰。三张图对应三个验收口径指向误差、峰值起伏、帧完整性答辩时按这个顺序讲逻辑链是完整的。具体验收可以写成一段可执行的检查脚本把关键指标打出来% 验收检查指向误差与峰值起伏 scan_list -60:2:60; peak_err zeros(size(scan_list)); for k 1:numel(scan_list) AF array_factor(16, 0.5, scan_list(k), -90:0.1:90); [~, idx] max(AF); peak_err(k) abs((-90:0.1:90)(idx) - scan_list(k)); end fprintf(最大指向误差: %.2f 度\n, max(peak_err)); assert(max(peak_err) 1, 指向误差超限检查相位符号);这段脚本把每个扫描角的实际峰值位置和设定值比对误差超过 1° 直接报错。答辩前跑一遍指标打印在屏幕上比口头说“我验证过了”有说服力得多。参数说明也简单阵元数 16、间距 0.5λ 是这套仿真的基准配置改这两个数误差指标会跟着变正好用来回答“阵元数对波束有什么影响”这类追问。答辩时还有一个高频问题为什么不用自适应波束形成。回答要落在需求边界上——这套仿真是相移扫描的原理演示自适应算法如 MVDR是另一个课题两者不冲突但目标不同硬塞进来反而让动图讲不清主线。把边界说清楚评委反而觉得你想明白了。从那以后我每次交付仿真类课设都强制先跑一遍验收脚本再导出动图指标不过就不往下走省得答辩前夜才发现符号写反。希望帮到你。本文还有配套的精品资源点击获取