简介：这份资源聚焦MIMO波束方向图与方向图特性分析，面向无线通信、雷达与5G相关方向的学习者和研究人员，帮助理解波束成形、波束赋形及FDA-MIMO系统的工作机理。内容围绕FDA模式、相控阵模式与FDA-MIMO模式三种工作方式展开，并涉及主瓣、旁瓣、后瓣等方向图关键概念，以及空间复用、抗衰落、能效优化等性能优势，适合作为理论学习与仿真验证的参考素材。资源包为rar压缩格式，仅含1个m文件，整体约3KB，该脚本可用于模拟或分析不同参数（如天线间距、相位偏移）对波束方向图的影响。目前已有408人学习下载，便于读者借助代码快速复现方向图计算流程，加深对MIMO波束特性的直观认识。
1. 从一份 11.m 说起：MIMO 波束方向图到底能拿来干什么
如果你手上正好有一份MIMO波束方向图.rar，解压后只看到一个11.m，第一反应大概率是「就这？」。我当初也一样。但真把这份脚本跑起来、把方向图一张张画出来之后，你会发现它其实是一个很典型的 FDA-MIMO 波束方向图仿真入口——用几十行 MATLAB 代码，把全维 MIMO、相控阵、FDA-MIMO 三种工作模式的方向图特性摆到同一张坐标系里对比。它解决的不是「怎么造天线」，而是「在动手做阵列之前，先在仿真里把波束形状、主瓣旁瓣、栅瓣位置摸清楚」。适合谁？做 5G 大规模阵列预研的、写雷达信号处理论文需要方向图配图的、以及刚接触 MIMO 波束成形想找个能跑的最小例子的从业者。方向图这三个字看着基础，但它决定了你后面所有波束赋形算法的上限。
2. 方向图与 FDA-MIMO 的底层逻辑：为什么先看主瓣旁瓣再谈算法
2.1 方向图三要素：主瓣、旁瓣、后瓣的工程含义
方向图本质上是把天线阵列的辐射强度写成空间角度的函数。对均匀线阵来说，阵因子可以写成
AF(θ) = Σ w_n · exp(j·2π·d/λ·n·sinθ)其中w_n是第 n 个阵元的加权系数，d是阵元间距，λ是波长，θ是扫描角。这个式子里藏了三个你必须在仿真里盯死的量：主瓣宽度、旁瓣电平、栅瓣位置。
主瓣是能量最集中的方向，它的宽度直接决定角分辨率——主瓣越窄，你越能把两个相邻目标分开。旁瓣是主瓣两侧的次强辐射，旁瓣电平高意味着你在接收有用信号的同时也把干扰收进来了，雷达里还会导致虚假目标。后瓣是天线背面的泄漏，通常希望压到 -20dB 以下。这三个量不是孤立的：阵元间距d变大，主瓣变窄，但d > λ/2时栅瓣就会冒出来，方向图出现多个和主瓣差不多高的假峰，这是新手最容易翻车的地方。
在11.m这类脚本里，你通常会看到对d/λ取 0.5、1.0、1.5 几组值做对比，目的就是让你亲眼看到栅瓣是怎么随着间距增大而出现的。我一般会先把d/λ = 0.5跑一遍作为基准，再逐步加大，观察第一个栅瓣出现在哪个角度——这个角度满足sinθ_g = ±λ/d，当d = λ时栅瓣正好在 ±90°，也就是端射方向。
2.2 FDA-MIMO 三种工作模式的差异与选型
FDA-MIMO 的核心思路是在阵元之间除了相位差之外再引入一个微小的频率偏移Δf，让不同阵元发射的信号频率略有不同。这样一来，波束不仅和角度有关，还和距离有关，形成所谓的「距离-角度耦合」。这带来一个传统相控阵没有的能力：可以在同一角度上区分不同距离的目标。
三种模式的区别在于频率偏移和相位加权的组合方式：
| 模式 | 频率偏移 | 相位加权 | 波束特点 | 典型用途 |
|---|---|---|---|---|
| FDA 模式 | 全部阵元启用 Δf | 无 | 距离-角度耦合，波束在距离维展开 | 距离维干扰抑制 |
| 相控阵模式 | Δf = 0 | 有 | 纯角度波束，主瓣窄增益高 | 定向通信、跟踪 |
| FDA-MIMO 模式 | 部分阵元 Δf | 有 | 角度维窄波束 + 距离维可控 | 灵活波束调度 |
选型逻辑很直接：如果你只需要在一个方向上集中能量，相控阵模式就够了，实现简单、增益高。如果你需要对抗来自同一方向但不同距离的干扰，FDA 模式的距离维自由度就是刚需。FDA-MIMO 模式是折中，代价是波束形成算法复杂度上升，因为你要同时优化角度和距离两个维度的加权。
在11.m里，这三种模式通常通过一个mode变量切换，或者写成三个独立的计算分支。你拿到脚本后第一件事就是找到这个切换点，确认每种模式对应的频率偏移向量f = [0, Δf, 2Δf, ...]是怎么构造的。常见做法是f = (0:N-1) * delta_f，其中delta_f一般取载频的千分之一到万分之一量级，太大了会导致波束在距离维散得太开，太小了距离分辨率又不够。
2.3 从阵因子到方向图：仿真里必须对齐的三个参数
跑任何方向图仿真之前，有三个参数必须和你的实际场景对齐，否则图再漂亮也是自嗨。
第一个是阵元数N。N越大主瓣越窄，但计算量和硬件成本也上去。11.m里大概率用的是 8 到 16 元线阵，这是教学和预研的常见规模。第二个是阵元间距d/λ，前面说了，超过 0.5 就有栅瓣风险，但毫米波频段因为波长小，物理间距可以做到很紧凑，这时候反而要关注互耦效应。第三个是扫描角范围，方向图通常只画 -90° 到 90°，但实际扫描到 ±60° 以外时主瓣会明显展宽、增益下降，这个现象叫波束偏斜，在宽角扫描场景里必须提前评估。
我一般会在脚本开头加一段参数检查：
% 参数检查：确保阵元间距不引发栅瓣 N = 16; % 阵元数 d_lambda = 0.5; % 阵元间距与波长之比 theta_scan = 0; % 扫描角（度） if d_lambda > 0.5 warning('阵元间距大于半波长，可能出现栅瓣，请检查方向图'); end % 计算理论栅瓣位置 if d_lambda > 0.5 theta_grating = asind(1/d_lambda); fprintf('理论栅瓣出现在 ±%.1f 度附近\n', theta_grating); end这段代码不复杂，但能帮你在跑图之前就预判问题。d_lambda是核心变量，asind(1/d_lambda)给出的是第一个栅瓣的理论角度，和仿真图对照着看，如果对不上，说明你的阵因子公式或者加权方式有问题。
3. 把 11.m 跑起来：从解压到出图的完整操作链
3.1 环境准备与脚本结构速读
11.m是 MATLAB 脚本，你需要一个能跑 MATLAB 的环境。如果没有正版授权，Octave 也能跑大部分基础绘图和矩阵运算，但涉及相控阵工具箱的函数可能会报错。我建议至少用 MATLAB R2020b 以上版本，因为11.m里可能用到phased命名空间下的函数，老版本没有。
解压MIMO波束方向图.rar后，目录结构通常很简单：一个11.m，可能还有一个data文件夹或者几个.mat文件。先用编辑器打开11.m，不要急着运行，花两分钟扫一遍结构。典型的脚本会分成四段：参数定义、阵因子计算、方向图绘制、结果保存。找到%%分隔符，MATLAB 的代码折叠就是靠它。
如果你看到类似load('xxx.mat')的语句，先确认那个.mat文件在不在当前工作目录。不在的话，要么把文件拷过去，要么改路径。这是最常见的「跑不起来」原因，没有之一。
3.2 逐段拆解：参数区、计算区、绘图区
参数区一般长这样：
%% 参数定义 c = 3e8; % 光速 fc = 10e9; % 载频 10GHz lambda = c/fc; % 波长 N = 16; % 阵元数 d = lambda/2; % 阵元间距 theta = -90:0.1:90; % 角度扫描范围 delta_f = 1e6; % FDA 频率偏移这里fc = 10e9是 X 波段，常见于雷达和卫星通信。delta_f = 1e6是 1MHz 的频率偏移，相对载频是万分之一，属于比较温和的设置。如果你把delta_f改成1e7，方向图在距离维的展开会明显得多，但角度维的主瓣可能会畸变。
计算区是核心，通常是一个循环或者矩阵运算：
%% 阵因子计算 AF = zeros(size(theta)); for n = 1:N % 相控阵模式：仅相位加权 phase_n = 2*pi*d/lambda*(n-1)*sind(theta); AF = AF + exp(1j*phase_n); end AF_db = 20*log10(abs(AF)/max(abs(AF))); % 归一化并转 dB这段代码算的是标准均匀线阵的阵因子。phase_n是第 n 个阵元在角度theta上的相位差，exp(1j*phase_n)是复加权。AF_db做了归一化，所以主瓣峰值是 0dB，旁瓣是负值。注意20*log10而不是10*log10，因为这是场强不是功率。
绘图区一般用plot或polarplot：
%% 方向图绘制 figure; plot(theta, AF_db, 'LineWidth', 1.5); xlabel('角度 (度)'); ylabel('归一化幅度 (dB)'); title('均匀线阵方向图'); grid on; ylim([-40, 0]);ylim([-40, 0])是为了把旁瓣细节看清楚，如果只画到 -20dB，很多旁瓣结构就看不到了。我习惯把纵轴压到 -40dB，这样能同时看到主瓣、第一旁瓣和后瓣的泄漏。
3.3 三种模式的方向图对比与参数扫描
要对比 FDA、相控阵、FDA-MIMO 三种模式，最直接的办法是在同一个figure里画三条曲线。你需要把计算区包一个for循环，每次改变频率偏移向量：
%% 三种模式对比 modes = {'相控阵', 'FDA', 'FDA-MIMO'}; delta_f_list = [0, 1e6, 1e6]; % 相控阵无频偏 phase_weight_list = {ones(1,N), ones(1,N), chebwin(N)'}; % 第三种用切比雪夫加权 figure; hold on; for m = 1:3 AF = zeros(size(theta)); delta_f = delta_f_list(m); w = phase_weight_list{m}; for n = 1:N phase_n = 2*pi*d/lambda*(n-1)*sind(theta); freq_n = 2*pi*delta_f*(n-1)*t; % t 为时间，需在循环外定义 AF = AF + w(n)*exp(1j*(phase_n + freq_n)); end AF_db = 20*log10(abs(AF)/max(abs(AF))); plot(theta, AF_db, 'LineWidth', 1.5, 'DisplayName', modes{m}); end legend; grid on; ylim([-40, 0]);这里有个坑：freq_n里的t是时间变量，在静态方向图里通常取t = 0或者一个固定时刻。如果你直接跑，t未定义会报错。常见做法是在参数区加t = 0;，表示在发射时刻的瞬时方向图。如果要看时间平均效果，就得对t做积分，那是另一套代码了。
参数扫描是这份脚本最值钱的地方。你可以固定其他参数，只改d/lambda从 0.5 到 2.0，步长 0.1，把每次的主瓣宽度和第一旁瓣电平记下来，画成曲线。这样你就能定量地看到：d/lambda = 0.5时主瓣最宽但无栅瓣；d/lambda = 1.0时主瓣窄了但栅瓣出现在 ±90°；d/lambda = 1.5时栅瓣已经进入可见区，方向图彻底乱掉。这个扫描过程比任何教科书上的图都直观。
4. 避坑与排查：方向图仿真里最容易翻车的五个地方
4.1 现象：方向图主瓣不在设定角度上
原因：相位加权符号搞反了。exp(1j*phase_n)和exp(-1j*phase_n)对应的是不同的扫描方向。如果你设了theta_scan = 30°，但主瓣跑到 -30°，就是符号问题。
解决：检查phase_n的表达式。对于扫描角theta_scan，正确的相位加权应该是phase_n = -2*pi*d/lambda*(n-1)*sind(theta_scan)，负号不能丢。在11.m里找到扫描角变量，确认它有没有被正确代入。
4.2 现象：旁瓣电平比理论值高很多
原因：没有做归一化，或者用了10*log10而不是20*log10。均匀线阵的第一旁瓣理论值是 -13.2dB，如果你画出来是 -6.6dB，那就是把场强当功率算了。
解决：统一用20*log10(abs(AF)/max(abs(AF)))。另外检查阵元加权是不是均匀的，如果你不小心加了窗函数（比如汉明窗），旁瓣会降但主瓣会宽，这是正常的。
4.3 现象：FDA 模式下方向图在距离维散不开
原因：delta_f太小，或者距离范围设得太窄。FDA 的距离分辨率取决于delta_f和阵元数N，近似有Δr ≈ c/(2*N*delta_f)。如果delta_f = 1e6、N = 16，Δr大约是 9.4 米，你的距离扫描范围如果只有几米，当然看不出展开。
解决：把delta_f提到1e7量级，或者把距离范围拉到几百米。在11.m里找距离向量r的定义，确认它的步长和范围是否匹配Δr。
4.4 现象：脚本报错「未定义函数或变量」
原因：缺少工具箱或者.mat文件不在路径里。phased.URA、phased.ULA这些需要 Phased Array System Toolbox，Octave 没有。
解决：先which一下报错的函数名，确认它属于哪个工具箱。如果是.mat文件缺失，用dir('*.mat')看看当前目录有什么，把路径补上。实在没有工具箱，就把phased相关调用替换成手写的阵因子公式，前面第 3 章给的代码就是替代方案。
4.5 现象：三种模式曲线完全重合
原因：delta_f没有真正作用到相位上，或者t始终为 0 导致频率偏移项2*pi*delta_f*(n-1)*t恒等于零。
解决：确认t的取值。如果t = 0，FDA 和相控阵在瞬时方向图上确实一样，因为频率偏移的效果需要时间积累。要看 FDA 的特性，要么取t > 0，要么直接看距离-角度二维图。常见做法是令t = 1/delta_f，这样2*pi*delta_f*t = 2*pi，频率偏移项刚好走完一个周期，差异最明显。
5. 进阶：把方向图从「能看」推到「能用」的两个技巧
5.1 用切比雪夫加权压低旁瓣并量化代价
均匀加权的第一旁瓣是 -13.2dB，很多场景不够用。切比雪夫加权可以在给定旁瓣电平的前提下让主瓣宽度最小，是工程上最常用的折中方案。在 MATLAB 里一行就能生成权向量：
% 切比雪夫加权：旁瓣电平设为 -30dB sll = -30; % 目标旁瓣电平 w = chebwin(N, -sll); % 注意 chebwin 的第二个参数是正数 w = w / max(w); % 归一化chebwin(N, -sll)里的-sll是因为 MATLAB 的chebwin要求输入正数表示旁瓣衰减量。生成权向量后，把它乘到阵因子的每一项上，再画方向图。你会看到旁瓣降到 -30dB，但主瓣宽度大约展宽了 1.5 倍。这个代价必须量化：如果你的应用要求角分辨率优于 5°，那 -30dB 旁瓣可能就太奢侈了，-20dB 是更常见的平衡点。
我一般会做一个扫描表，把不同旁瓣电平下的主瓣宽度列出来，选一个刚好满足系统指标的。这个表比任何理论推导都管用，因为它是基于你自己的阵元数和间距算出来的。
5.2 从静态方向图到波束扫描的验证方法
静态方向图只告诉你「波束长什么样」，但实际系统里波束是要扫描的。验证扫描性能的方法很简单：固定加权向量，只改扫描角theta_scan从 -60° 到 60°，步长 10°，每次画一条方向图，叠在一起看。
% 波束扫描验证 theta_scan_list = -60:10:60; figure; hold on; for ts = theta_scan_list AF = zeros(size(theta)); for n = 1:N phase_n = -2*pi*d/lambda*(n-1)*(sind(theta) - sind(ts)); AF = AF + exp(1j*phase_n); end AF_db = 20*log10(abs(AF)/max(abs(AF))); plot(theta, AF_db, 'DisplayName', sprintf('扫描角 %d°', ts)); end legend('Location', 'best'); grid on; ylim([-40, 0]);注意phase_n里用的是sind(theta) - sind(ts)，这是波束扫描的标准形式。跑完这张图，你会看到两个现象：一是主瓣随着ts移动，二是在大扫描角（比如 ±60°）时主瓣明显展宽、增益下降。这个展宽量就是波束偏斜，它和d/lambda成正比。如果你发现 ±60° 时主瓣宽到无法接受，就得减小d/lambda或者接受更小的扫描范围。
从那以后我每次拿到一份新的方向图脚本，都会先跑一遍扫描验证，把主瓣宽度随角度的变化曲线画出来，再决定这个阵列能不能覆盖我需要的空域。这个习惯帮我省掉了至少两次在硬件做完之后才发现扫描范围不够的返工。希望这份拆解能帮到你，把11.m里的每一行都变成你自己阵列设计里的底气。
本文还有配套的精品资源，点击获取
