简介面向雷达、声纳与无线通信领域的信号处理学习者和研究者这份频率多样性FDA波束形成程序包聚焦多载频体制下的波束形成与干扰抑制通过频率分集和自适应算法解决复杂多径、频率选择性衰落环境下的目标检测与定位问题。压缩包共15个文件包含9个M脚本用于发射/接收方向图仿真、自适应波束形成及信干噪比对比3个FIG图形文件便于直观查看波束方向图3个初始化参数文件用于配置多载频环境整体仅1.95MB结构紧凑。已有714人学习下载。通过运行脚本可复现频率多样性阵列的发射/接收方向图对比不同频率间距与波长间距下的信干噪比表现并观察自适应波束对移动目标的跟踪与旁瓣抑制效果脚本覆盖发射、接收、自适应处理与性能对比等模块便于按需调用修改适合作为课程设计、毕业设计或科研验证的参考代码。1. 先谈清楚FDA 波束形成到底解决了什么问题做阵列信号处理的人大概率都遇到过这个尴尬相控阵能把波束指到某个角度但它分不清这个角度上到底是 10 公里外的目标还是 50 公里外的目标——距离维信息是缺失的。FDAFrequency Diverse Array频控阵就是冲着这个缺口来的它给每个阵元加一个微小的频率偏置让波束方向图同时依赖角度和距离于是一个波束能覆盖的距离-角度二维区域被重新塑形。我最早接触这份资源是因为要给某跨平台系统做近场干扰抑制的预研结果发现网上能直接跑的 FDA 仿真代码并不多大多停留在公式推导阶段。这份资源的意义就在于它把 FDA 发射波束形成、接收端多波束处理完整串起来了适合刚接触频控阵的从业者快速搭起仿真基线也适合做相控阵的人理解距离维波束控制的边界在哪。2. FDA 波束形成的基本原理频率偏置如何撑起距离-角度二维自由度2.1 从相控阵到频控阵那个“微小频偏”到底改了什么相控阵的波束指向只由阵元间的相位差决定公式里没有距离项。FDA 的做法是在每个阵元的工作频率上叠加一个线性递增的频偏 Δf第 n 个阵元的发射频率变成 f0 n·Δf。这个频偏看起来很小通常只有 kHz 到 MHz 量级但它让每个阵元发射的波形在传播过程中积累了不同的相位——注意这个相位差不仅与角度 θ 有关还与目标距离 r 有关。因为频率不同电磁波传播相同距离后相位旋转的圈数不同于是远场某一点r, θ处的合成场强同时受距离和角度控制。这带来的直接后果是波束最大值不再是一条从阵列法线出发的射线而是一条随距离弯曲的曲线传统相控阵的“角度波束”在这里变成了“距离-角度二维响应”。这份资源里给出的 FDA 方向图表达式核心是阵因子项 exp(j·2π·(n·Δf·(r/c) - n·d·sinθ/λ))。注意这个 r/c 项它就是距离维自由度的来源。实际仿真时如果 Δf 选得太小距离维的相位积累不够波束在距离上几乎没有变化退化成相控阵选得太大波束在距离上剧烈摆动可能出现距离模糊。常见的做法是让 Δf 满足 Δf 1/TT 为脉冲宽度保证窄带假设成立。我一般会先固定 Δf 10 kHz、阵元数 N 16、载频 f0 10 GHz 做第一轮扫描观察主瓣随距离的弯曲情况——资源里的参数设置文件也是这个量级。2.2 发射波束与接收波束的差异为什么接收端不能直接套用发射模型这份资源同时给出了发射波束形成和接收波束形成的代码路径但两者处理逻辑完全不同。发射端FDA 通过频偏自然产生距离-角度耦合的波束不需要额外加窗。接收端则要复杂得多因为每个阵元的接收信号已经包含了发射频偏的传播延迟影响所以接收导向矢量必须同时补偿频率差和空间延迟。换句话说接收端做常规延时求和时如果沿用相控阵的导向矢量只含 sinθ 项波束会在距离维上偏掉。资源里有一个关键处理接收端导向矢量写作 s(r, θ) [1, exp(-j·(2π·Δf·r/c 2π·d·sinθ/λ)), ...]注意这个 r 是待聚焦距离。如果处理的是近场目标r 是实际距离如果处理远场r 趋于无穷此时 Δf·r/c 项变化极快导致波束在距离维严重失真——这就是为什么很多论文里强调 FDA 接收波束形成更适合中近场场景。我在复现时发现资源中接收端代码默认对 5 km 以内的目标做了距离匹配超过这个距离输出就开始出现副瓣抬高。所以实验之前要根据自己的场景先改 r 的扫描范围和步进。这个点不搞清楚直接跑接收端多波束模块大概率会得到一团乱的方向图。3. 用 MATLAB 复现 FDA 发射波束形成参数设计与距离-角度扫描3.1 基础参数初始化与阵元相位计算资源里给的第一个可执行脚本是发射波束方向图计算核心逻辑是遍历距离和角度网格累计每个阵元的相位贡献。运行时先初始化常量然后计算每个阵元的相位修正项。下面是我在复现时整理出的最小可运行结构% FDA发射波束形成 - 最小复现示例 % 载频10GHz频偏10kHz16阵元阵元间距半波长 clear; clc; close all; f0 10e9; % 载频 (Hz) delta_f 10e3; % 频偏 (Hz) N 16; % 阵元数 d 0.5 * (3e8 / f0); % 阵元间距半波长 c 3e8; % 光速 % 扫描网格 r linspace(5e3, 50e3, 200); % 距离 5km 到 50km theta linspace(-60, 60, 200) * pi/180; % 角度 -60° 到 60° % 计算阵列因子 AF(r, theta) AF zeros(length(r), length(theta)); for ri 1:length(r) for ti 1:length(theta) sum_val 0; for n 0:N-1 % 相位项1频偏导致的距离-角度耦合相位 phase_freq 2 * pi * n * delta_f * (r(ri)/c - d*sin(theta(ti))/c); % 相位项2常规空间导向相位 phase_space -2 * pi * n * d * sin(theta(ti)) / (c/f0); sum_val sum_val exp(1j * (phase_freq phase_space)); end AF(ri, ti) abs(sum_val); end end % 归一化并绘图 AF_db 20*log10(AF / max(AF(:)) 1e-12); figure; imagesc(theta*180/pi, r/1e3, AF_db); xlabel(角度 (deg)); ylabel(距离 (km)); title(FDA发射波束方向图 (16元, Δf10kHz)); colorbar; caxis([-40 0]);这段代码的循环结构是我刻意保留的因为 FDA 方向图计算的核心就在逐个阵元叠加相位向量化写法虽然快但不利于理解相位累加过程。参数说明delta_f是频偏决定了距离维波束的弯曲速度r(ri)/c是电磁波传播到该距离所需的时间这个量乘上频偏得到额外的相位旋转phase_space是常规相控阵的空间相位项用于控制角度的主瓣指向。注意我特意把角度扫描范围设成 -60° 到 60°因为 FDA 波束在超过 ±60° 后会出现严重的栅瓣周期性重复这是频偏项周期性导致的不是代码错误。3.2 波束扫描与方向图判读主瓣弯曲、周期性和距离模糊复现完这版方向图你会看到主瓣不再是一条竖直的亮线而是从某个角度出发、随距离增加逐渐弯曲的弧线。这个弯曲方向取决于频偏的正负Δf 为正时波束向负角度弯曲Δf 为负时向正角度弯曲。资源里提到了一个常见疑问——“为什么我扫描出来的主瓣在某个距离突然分裂”答案是 FDA 的阵列因子在频率维本身是周期函数当频率偏置导致的相位差累计超过 2π 时会出现第二个主瓣也就是距离维的栅瓣。周期由 c/Δf 决定比如 Δf 10 kHz 时距离周期是 30 km所以在你扫描 5-50 km 范围内大约在 30 km 附近能看到重复主瓣。判断方向图是否合理的三个经验标准一是主瓣峰值是否在预设的参考距离处通常设为 10 km二是主瓣弯曲是否平滑没有突变跳变三是副瓣电平是否低于 -10 dB。如果副瓣过高优先查两个地方一个是阵元间距是否严格半波长另一个是角度网格采样是否过粗导致峰值被截断。我跑通后会把delta_f改成 5 kHz、20 kHz 各跑一遍观察主瓣弯曲程度和副瓣变化这对理解频偏的作用比读十遍理论推导都来得直接。4. 接收波束形成与多波束处理MVDR 自适应在 FDA 下的适配4.1 接收导向矢量的修正频率差补偿是成败关键发射端只需要正向计算相位接收端则需要额外考虑每个阵元收到的信号已经携带了频率差带来的传播相位历史。这意味着匹配滤波或波束形成时导向矢量不再只是空间角度的函数。资源里接收端模块的做法是把接收导向矢量拆成两个部分频偏导致的距离相位项以及空间相位项。关键在于距离项中的 r 是未知的需要扫描或估计。下面这段代码是资源接收模块里最核心的部分——构造修正后的导向矢量% 接收导向矢量构造 - 包含FDA频率差补偿 function a fda_steering_vector(theta, r, f0, delta_f, d, N) c 3e8; lambda c / f0; % 空间相位项第n个阵元相对第0个阵元的相位差 phase_space -2 * pi * d * sin(theta) * (0:N-1). / lambda; % FDA频偏导致的距离相位项 phase_freq 2 * pi * delta_f * (0:N-1). * (r / c - d * sin(theta) / c); % 合成导向矢量 a exp(1j * (phase_space phase_freq)); end为什么不能忽略phase_freq做一个简单的数值实验就明白了当 r 10 kmN 16Δf 10 kHz 时频偏引起的最大相位差是 2π × 10kHz × (10km / c) ≈ 2.09 rad这已经足够让主瓣指向偏移好几度了。如果接收端还用纯空间相位构造导向矢量波束形成器的输出会在错误方向形成峰值。另一个细节是这个修正项里 (0:N-1). 是列向量theta 和 r 是标量保证 a 是 N×1 的列向量。资源里有一种错误版本是把 r 写成向量结果导向矢量变成矩阵运算时报维度不匹配——这就是典型的 MATLAB 行列陷阱。4.2 MVDR 波束形成器在 FDA 下的协方差矩阵构造FDA 接收端最常见的自适应处理是 MVDR最小方差无失真响应但直接套用相控阵 MVDR 代码会翻车因为期望信号的导向矢量不对。正确流程是先构造包含目标距离 r0 和角度 θ0 的修正导向矢量 a再估计采样协方差矩阵 R最后求解权重。注意协方差矩阵一定不能包含期望信号本身否则会导致目标自消——这就是为什么码里要单独留一段“保护间隔”或“不含期望信号的快拍”来估计干扰加噪声协方差。% FDA-MVDR接收波束形成 % 场景期望信号 (10km, 10°)干扰 (8km, -20°) r0 10e3; theta0 10*pi/180; % 期望信号位置 rj 8e3; thetaj -20*pi/180; % 干扰位置 d 0.015; f0 10e9; delta_f 10e3; N 16; L 500; % 快拍数 % 生成接收数据期望干扰噪声 a_target fda_steering_vector(theta0, r0, f0, delta_f, d, N); a_jammer fda_steering_vector(thetaj, rj, f0, delta_f, d, N); snr 10; jnr 20; s_target sqrt(10^(snr/10)) * exp(1j*randn(L,1)); s_jammer sqrt(10^(jnr/10)) * exp(1j*randn(L,1)); X a_target * s_target. a_jammer * s_jammer. (randn(N,L) 1j*randn(N,L))/sqrt(2); % 估计协方差不含期望信号的快拍段 R (X(:, 1:100) * X(:, 1:100)) / 100 1e-6 * eye(N); % MVDR权重 a_desired fda_steering_vector(theta0, r0, f0, delta_f, d, N); w (R \ a_desired) / (a_desired * (R \ a_desired)); % 输出功率扫描验证 r_scan linspace(5e3, 15e3, 100); theta_scan linspace(-40, 40, 100) * pi/180; P zeros(length(r_scan), length(theta_scan)); for ri 1:length(r_scan) for ti 1:length(theta_scan) a_s fda_steering_vector(theta_scan(ti), r_scan(ri), f0, delta_f, d, N); P(ri, ti) abs(w * a_s)^2; end end % 绘图部分略P 的最大值应出现在 (10km, 10°) 附近参数说明L 500是快拍数snr/jnr是信噪比和干噪比dB。X(:, 1:100)取前 100 个快拍估计协方差这部分不含期望信号——这里有个潜在问题MVDR 的协方差估计窗口和信号起点重合如果干扰也在这个窗口内估计是合理的但实际系统里干扰起止时间可能不固定更稳健的做法是取中间段或采用对角加载。代码里 1e-6 * eye(N)就是对角加载项用来防止协方差奇异和低快拍数时的估计误差放大。如果你在实验中发现波束方向图在主瓣附近出现凹陷第一反应应该是调大对角加载系数而不是怀疑导向矢量写错。4.3 多波束输出与距离维“距离相关”的波束现象资源标题里强调“多波束”在接收处理里具体体现在同一个角度上不设置一个波束而是沿距离维设置多个波束。实现方式很直接对一组预先设定的距离 r1, r2, ..., rM 分别构造导向矢量fda_steering_vector(theta0, rm, ...)然后对同一份快拍数据分别做加权求和。多波束的意义是解决 FDA 波束“距离-角度耦合”带来的含糊实际目标距离未知时单波束无法确定距离而多波束同时输出各距离候选的响应强度峰值所在距离即目标估计。我通常设置 3-5 个等间距距离波束比如 5km、10km、15km并画在同一张距离-输出功率图上观察三个波束的输出曲线是否存在明显的“重叠区”或“塌陷区”。塌陷区通常出现在两波束距离的中点原因是 FDA 导向矢量在该距离上对两个波束的相关性都高导致输出功率下降。资源里有一种补救办法增大频偏 Δf 可以拉开波束间的距离选择性但代价是距离周期性变短更容易出现栅瓣。所以多波束的距离间隔不是随意定的经验值是以 c/(2·Δf·N) 为基准比如上面参数下该值约 0.94 km所以 5/10/15 km 的间隔5km远大于栅瓣间隔既能区分距离又不会波束重叠严重。5. 避坑与常见问题FDA 仿真中 5 个最容易翻车的细节5.1 频偏符号导致波束弯曲方向相反现象跑完发射方向图发现主瓣向正角度弯曲但理论推导中预期的是负角度。原因phase_freq里delta_f的正负号定义不统一资源里某一版本的频偏定义为f f0 n*delta_f另一版本定义为f f0 - n*delta_f。解决先固定为正频偏跑一遍观察主瓣弯曲方向再根据需求调整符号。从用户视角来看弯曲方向本身没有绝对对错只要和你后续的接收端导向矢量保持一致即可。我一般会在代码头部写一行注释明确“正向频偏波束向负角度弯曲”。5.2 距离扫描范围超过周期导致栅瓣误判为目标现象方向图在某个非目标距离出现与主瓣等高的亮带。原因FDA 阵列因子的距离周期性导致主瓣重复周期 c/Δf。解决把扫描范围控制在 0 到 c/Δf 之内。如果必须扫描更远距离就在结果图里标注周期边界提醒自己那是重复主瓣而不是真实目标。5.3 MATLAB 维度不匹配导向矢量被写成矩阵现象MVDR 权重计算时报错“维度不匹配”或得到非 1×N 权重向量。原因导向矢量函数里(0:N-1).和theta/r的维度没有对齐最常见的是把 r 定义为向量而不是标量。解决强制要求导向矢量函数第一行把输入标量化——在函数入口加theta theta(1); r r(1);。养成这个习惯后十次有八次维度报错会消失。5.4 快拍数过少导致协方差矩阵奇异现象MVDR 输出方向图出现剧烈抖动主瓣偏离预期位置。原因快拍数 L N 时采样协方差矩阵秩亏求逆结果被噪声主导。解决至少保证 L 2N通常 L 500N 16 足够同时加对角加载系数 1e-6 到 1e-3 之间。从经验来看对角加载加得太大会导致自适应能力下降零陷变浅加得太小低快拍时算法不稳定。5.5 副瓣电平异常高阵元间距违反半波长约束现象方向图除主瓣外在 ±30° 附近出现非对称高副瓣。原因阵元间距 d 超过半波长导致空间栅瓣混入。解决用lambda c/f0; d 0.5*lambda;严格计算不要手输近似值。我遇到过一次把 d 手输成 0.02 m大约 0.67 倍波长副瓣瞬间抬高了 8 dB当时以为是算法问题排查了半天才发现是采样间距写错。6. 进阶技巧把 FDA 波束做成多目标搜索模式到这里你已经能跑通 FDA 发射和接收的基本流程但资源里最有价值的部分其实是如何把 FDA 的“距离相关”能力用到多目标搜索上。传统相控阵要覆盖一个距离-角度区域需要多帧扫描FDA 因为单帧波束本身就覆盖一条距离-角度耦合轨迹所以理论上单帧就能对特定耦合轨迹上的目标进行搜索。实操方法是设定一个期望角度范围对每个角度计算对应的耦合距离 r(θ) r0 R0·sinθR0 为参考距离然后在这个轨迹上设置密集的波束形成“角度-距离联合搜索线”。% 沿FDA天然波束轨迹设置搜索波束 theta_traj linspace(-30, 30, 50) * pi/180; r_ref 15e3; % 参考距离 r_traj r_ref * (1 sin(theta_traj)); % 耦合轨迹近似 search_power zeros(1, length(theta_traj)); for idx 1:length(theta_traj) a fda_steering_vector(theta_traj(idx), r_traj(idx), f0, delta_f, d, N); search_power(idx) abs(w * a)^2; % w 来自上一节的MVDR权重 end % 峰值位置对应目标的角度随后可反查距离 [~, peak_idx] max(search_power); target_theta_est theta_traj(peak_idx) * 180/pi; target_r_est r_traj(peak_idx); fprintf(估计目标: 角度 %.2f°, 距离 %.2f km\n, ... target_theta_est, target_r_est/1e3);这段代码的关键在于r_traj与theta_traj的耦合关系必须与 FDA 实际的波束弯曲方向一致——如果频偏为正导致波束负角度弯曲上面公式的号就要改成-。之前的避坑要点在这里直接生效。我最初做这个多目标搜索时因为忘了调符号估计值在角度上偏移了大约 8°后来统一用一个函数生成轨迹r_traj r_ref * (1 sign(delta_f) * sin(theta_traj))再也没出过偏移问题。多目标场景下可以同时保留两组搜索波束正频偏和负频偏分别搜索正角度区域和负角度区域输出功率叠加后取两组的局部峰值。要注意的是两组波束在距离维的重复周期不同为了避免误检要设置一个距离窗只保留第一周期内的峰值。从那以后我每次做 FDA 仿真都强制走一遍“先定义频偏符号 → 再画轨迹 → 最后搜索”的流程搜索角度结果基本一次到位不会再出现方向对但距离偏的尴尬情况。这份资源对需要做频控阵预研或算法对比的人来说至少省了半个月的前期摸索希望帮到你。本文还有配套的精品资源点击获取 
长河君 · 2026/10/14 7:35:56 · 阅读 1.2k FDA频控阵波束形成：从距离-角度二维自由度到MVDR接收仿真
简介面向雷达、声纳与无线通信领域的信号处理学习者和研究者这份频率多样性FDA波束形成程序包聚焦多载频体制下的波束形成与干扰抑制通过频率分集和自适应算法解决复杂多径、频率选择性衰落环境下的目标检测与定位问题。压缩包共15个文件包含9个M脚本用于发射/接收方向图仿真、自适应波束形成及信干噪比对比3个FIG图形文件便于直观查看波束方向图3个初始化参数文件用于配置多载频环境整体仅1.95MB结构紧凑。已有714人学习下载。通过运行脚本可复现频率多样性阵列的发射/接收方向图对比不同频率间距与波长间距下的信干噪比表现并观察自适应波束对移动目标的跟踪与旁瓣抑制效果脚本覆盖发射、接收、自适应处理与性能对比等模块便于按需调用修改适合作为课程设计、毕业设计或科研验证的参考代码。1. 先谈清楚FDA 波束形成到底解决了什么问题做阵列信号处理的人大概率都遇到过这个尴尬相控阵能把波束指到某个角度但它分不清这个角度上到底是 10 公里外的目标还是 50 公里外的目标——距离维信息是缺失的。FDAFrequency Diverse Array频控阵就是冲着这个缺口来的它给每个阵元加一个微小的频率偏置让波束方向图同时依赖角度和距离于是一个波束能覆盖的距离-角度二维区域被重新塑形。我最早接触这份资源是因为要给某跨平台系统做近场干扰抑制的预研结果发现网上能直接跑的 FDA 仿真代码并不多大多停留在公式推导阶段。这份资源的意义就在于它把 FDA 发射波束形成、接收端多波束处理完整串起来了适合刚接触频控阵的从业者快速搭起仿真基线也适合做相控阵的人理解距离维波束控制的边界在哪。2. FDA 波束形成的基本原理频率偏置如何撑起距离-角度二维自由度2.1 从相控阵到频控阵那个“微小频偏”到底改了什么相控阵的波束指向只由阵元间的相位差决定公式里没有距离项。FDA 的做法是在每个阵元的工作频率上叠加一个线性递增的频偏 Δf第 n 个阵元的发射频率变成 f0 n·Δf。这个频偏看起来很小通常只有 kHz 到 MHz 量级但它让每个阵元发射的波形在传播过程中积累了不同的相位——注意这个相位差不仅与角度 θ 有关还与目标距离 r 有关。因为频率不同电磁波传播相同距离后相位旋转的圈数不同于是远场某一点r, θ处的合成场强同时受距离和角度控制。这带来的直接后果是波束最大值不再是一条从阵列法线出发的射线而是一条随距离弯曲的曲线传统相控阵的“角度波束”在这里变成了“距离-角度二维响应”。这份资源里给出的 FDA 方向图表达式核心是阵因子项 exp(j·2π·(n·Δf·(r/c) - n·d·sinθ/λ))。注意这个 r/c 项它就是距离维自由度的来源。实际仿真时如果 Δf 选得太小距离维的相位积累不够波束在距离上几乎没有变化退化成相控阵选得太大波束在距离上剧烈摆动可能出现距离模糊。常见的做法是让 Δf 满足 Δf 1/TT 为脉冲宽度保证窄带假设成立。我一般会先固定 Δf 10 kHz、阵元数 N 16、载频 f0 10 GHz 做第一轮扫描观察主瓣随距离的弯曲情况——资源里的参数设置文件也是这个量级。2.2 发射波束与接收波束的差异为什么接收端不能直接套用发射模型这份资源同时给出了发射波束形成和接收波束形成的代码路径但两者处理逻辑完全不同。发射端FDA 通过频偏自然产生距离-角度耦合的波束不需要额外加窗。接收端则要复杂得多因为每个阵元的接收信号已经包含了发射频偏的传播延迟影响所以接收导向矢量必须同时补偿频率差和空间延迟。换句话说接收端做常规延时求和时如果沿用相控阵的导向矢量只含 sinθ 项波束会在距离维上偏掉。资源里有一个关键处理接收端导向矢量写作 s(r, θ) [1, exp(-j·(2π·Δf·r/c 2π·d·sinθ/λ)), ...]注意这个 r 是待聚焦距离。如果处理的是近场目标r 是实际距离如果处理远场r 趋于无穷此时 Δf·r/c 项变化极快导致波束在距离维严重失真——这就是为什么很多论文里强调 FDA 接收波束形成更适合中近场场景。我在复现时发现资源中接收端代码默认对 5 km 以内的目标做了距离匹配超过这个距离输出就开始出现副瓣抬高。所以实验之前要根据自己的场景先改 r 的扫描范围和步进。这个点不搞清楚直接跑接收端多波束模块大概率会得到一团乱的方向图。3. 用 MATLAB 复现 FDA 发射波束形成参数设计与距离-角度扫描3.1 基础参数初始化与阵元相位计算资源里给的第一个可执行脚本是发射波束方向图计算核心逻辑是遍历距离和角度网格累计每个阵元的相位贡献。运行时先初始化常量然后计算每个阵元的相位修正项。下面是我在复现时整理出的最小可运行结构% FDA发射波束形成 - 最小复现示例 % 载频10GHz频偏10kHz16阵元阵元间距半波长 clear; clc; close all; f0 10e9; % 载频 (Hz) delta_f 10e3; % 频偏 (Hz) N 16; % 阵元数 d 0.5 * (3e8 / f0); % 阵元间距半波长 c 3e8; % 光速 % 扫描网格 r linspace(5e3, 50e3, 200); % 距离 5km 到 50km theta linspace(-60, 60, 200) * pi/180; % 角度 -60° 到 60° % 计算阵列因子 AF(r, theta) AF zeros(length(r), length(theta)); for ri 1:length(r) for ti 1:length(theta) sum_val 0; for n 0:N-1 % 相位项1频偏导致的距离-角度耦合相位 phase_freq 2 * pi * n * delta_f * (r(ri)/c - d*sin(theta(ti))/c); % 相位项2常规空间导向相位 phase_space -2 * pi * n * d * sin(theta(ti)) / (c/f0); sum_val sum_val exp(1j * (phase_freq phase_space)); end AF(ri, ti) abs(sum_val); end end % 归一化并绘图 AF_db 20*log10(AF / max(AF(:)) 1e-12); figure; imagesc(theta*180/pi, r/1e3, AF_db); xlabel(角度 (deg)); ylabel(距离 (km)); title(FDA发射波束方向图 (16元, Δf10kHz)); colorbar; caxis([-40 0]);这段代码的循环结构是我刻意保留的因为 FDA 方向图计算的核心就在逐个阵元叠加相位向量化写法虽然快但不利于理解相位累加过程。参数说明delta_f是频偏决定了距离维波束的弯曲速度r(ri)/c是电磁波传播到该距离所需的时间这个量乘上频偏得到额外的相位旋转phase_space是常规相控阵的空间相位项用于控制角度的主瓣指向。注意我特意把角度扫描范围设成 -60° 到 60°因为 FDA 波束在超过 ±60° 后会出现严重的栅瓣周期性重复这是频偏项周期性导致的不是代码错误。3.2 波束扫描与方向图判读主瓣弯曲、周期性和距离模糊复现完这版方向图你会看到主瓣不再是一条竖直的亮线而是从某个角度出发、随距离增加逐渐弯曲的弧线。这个弯曲方向取决于频偏的正负Δf 为正时波束向负角度弯曲Δf 为负时向正角度弯曲。资源里提到了一个常见疑问——“为什么我扫描出来的主瓣在某个距离突然分裂”答案是 FDA 的阵列因子在频率维本身是周期函数当频率偏置导致的相位差累计超过 2π 时会出现第二个主瓣也就是距离维的栅瓣。周期由 c/Δf 决定比如 Δf 10 kHz 时距离周期是 30 km所以在你扫描 5-50 km 范围内大约在 30 km 附近能看到重复主瓣。判断方向图是否合理的三个经验标准一是主瓣峰值是否在预设的参考距离处通常设为 10 km二是主瓣弯曲是否平滑没有突变跳变三是副瓣电平是否低于 -10 dB。如果副瓣过高优先查两个地方一个是阵元间距是否严格半波长另一个是角度网格采样是否过粗导致峰值被截断。我跑通后会把delta_f改成 5 kHz、20 kHz 各跑一遍观察主瓣弯曲程度和副瓣变化这对理解频偏的作用比读十遍理论推导都来得直接。4. 接收波束形成与多波束处理MVDR 自适应在 FDA 下的适配4.1 接收导向矢量的修正频率差补偿是成败关键发射端只需要正向计算相位接收端则需要额外考虑每个阵元收到的信号已经携带了频率差带来的传播相位历史。这意味着匹配滤波或波束形成时导向矢量不再只是空间角度的函数。资源里接收端模块的做法是把接收导向矢量拆成两个部分频偏导致的距离相位项以及空间相位项。关键在于距离项中的 r 是未知的需要扫描或估计。下面这段代码是资源接收模块里最核心的部分——构造修正后的导向矢量% 接收导向矢量构造 - 包含FDA频率差补偿 function a fda_steering_vector(theta, r, f0, delta_f, d, N) c 3e8; lambda c / f0; % 空间相位项第n个阵元相对第0个阵元的相位差 phase_space -2 * pi * d * sin(theta) * (0:N-1). / lambda; % FDA频偏导致的距离相位项 phase_freq 2 * pi * delta_f * (0:N-1). * (r / c - d * sin(theta) / c); % 合成导向矢量 a exp(1j * (phase_space phase_freq)); end为什么不能忽略phase_freq做一个简单的数值实验就明白了当 r 10 kmN 16Δf 10 kHz 时频偏引起的最大相位差是 2π × 10kHz × (10km / c) ≈ 2.09 rad这已经足够让主瓣指向偏移好几度了。如果接收端还用纯空间相位构造导向矢量波束形成器的输出会在错误方向形成峰值。另一个细节是这个修正项里 (0:N-1). 是列向量theta 和 r 是标量保证 a 是 N×1 的列向量。资源里有一种错误版本是把 r 写成向量结果导向矢量变成矩阵运算时报维度不匹配——这就是典型的 MATLAB 行列陷阱。4.2 MVDR 波束形成器在 FDA 下的协方差矩阵构造FDA 接收端最常见的自适应处理是 MVDR最小方差无失真响应但直接套用相控阵 MVDR 代码会翻车因为期望信号的导向矢量不对。正确流程是先构造包含目标距离 r0 和角度 θ0 的修正导向矢量 a再估计采样协方差矩阵 R最后求解权重。注意协方差矩阵一定不能包含期望信号本身否则会导致目标自消——这就是为什么码里要单独留一段“保护间隔”或“不含期望信号的快拍”来估计干扰加噪声协方差。% FDA-MVDR接收波束形成 % 场景期望信号 (10km, 10°)干扰 (8km, -20°) r0 10e3; theta0 10*pi/180; % 期望信号位置 rj 8e3; thetaj -20*pi/180; % 干扰位置 d 0.015; f0 10e9; delta_f 10e3; N 16; L 500; % 快拍数 % 生成接收数据期望干扰噪声 a_target fda_steering_vector(theta0, r0, f0, delta_f, d, N); a_jammer fda_steering_vector(thetaj, rj, f0, delta_f, d, N); snr 10; jnr 20; s_target sqrt(10^(snr/10)) * exp(1j*randn(L,1)); s_jammer sqrt(10^(jnr/10)) * exp(1j*randn(L,1)); X a_target * s_target. a_jammer * s_jammer. (randn(N,L) 1j*randn(N,L))/sqrt(2); % 估计协方差不含期望信号的快拍段 R (X(:, 1:100) * X(:, 1:100)) / 100 1e-6 * eye(N); % MVDR权重 a_desired fda_steering_vector(theta0, r0, f0, delta_f, d, N); w (R \ a_desired) / (a_desired * (R \ a_desired)); % 输出功率扫描验证 r_scan linspace(5e3, 15e3, 100); theta_scan linspace(-40, 40, 100) * pi/180; P zeros(length(r_scan), length(theta_scan)); for ri 1:length(r_scan) for ti 1:length(theta_scan) a_s fda_steering_vector(theta_scan(ti), r_scan(ri), f0, delta_f, d, N); P(ri, ti) abs(w * a_s)^2; end end % 绘图部分略P 的最大值应出现在 (10km, 10°) 附近参数说明L 500是快拍数snr/jnr是信噪比和干噪比dB。X(:, 1:100)取前 100 个快拍估计协方差这部分不含期望信号——这里有个潜在问题MVDR 的协方差估计窗口和信号起点重合如果干扰也在这个窗口内估计是合理的但实际系统里干扰起止时间可能不固定更稳健的做法是取中间段或采用对角加载。代码里 1e-6 * eye(N)就是对角加载项用来防止协方差奇异和低快拍数时的估计误差放大。如果你在实验中发现波束方向图在主瓣附近出现凹陷第一反应应该是调大对角加载系数而不是怀疑导向矢量写错。4.3 多波束输出与距离维“距离相关”的波束现象资源标题里强调“多波束”在接收处理里具体体现在同一个角度上不设置一个波束而是沿距离维设置多个波束。实现方式很直接对一组预先设定的距离 r1, r2, ..., rM 分别构造导向矢量fda_steering_vector(theta0, rm, ...)然后对同一份快拍数据分别做加权求和。多波束的意义是解决 FDA 波束“距离-角度耦合”带来的含糊实际目标距离未知时单波束无法确定距离而多波束同时输出各距离候选的响应强度峰值所在距离即目标估计。我通常设置 3-5 个等间距距离波束比如 5km、10km、15km并画在同一张距离-输出功率图上观察三个波束的输出曲线是否存在明显的“重叠区”或“塌陷区”。塌陷区通常出现在两波束距离的中点原因是 FDA 导向矢量在该距离上对两个波束的相关性都高导致输出功率下降。资源里有一种补救办法增大频偏 Δf 可以拉开波束间的距离选择性但代价是距离周期性变短更容易出现栅瓣。所以多波束的距离间隔不是随意定的经验值是以 c/(2·Δf·N) 为基准比如上面参数下该值约 0.94 km所以 5/10/15 km 的间隔5km远大于栅瓣间隔既能区分距离又不会波束重叠严重。5. 避坑与常见问题FDA 仿真中 5 个最容易翻车的细节5.1 频偏符号导致波束弯曲方向相反现象跑完发射方向图发现主瓣向正角度弯曲但理论推导中预期的是负角度。原因phase_freq里delta_f的正负号定义不统一资源里某一版本的频偏定义为f f0 n*delta_f另一版本定义为f f0 - n*delta_f。解决先固定为正频偏跑一遍观察主瓣弯曲方向再根据需求调整符号。从用户视角来看弯曲方向本身没有绝对对错只要和你后续的接收端导向矢量保持一致即可。我一般会在代码头部写一行注释明确“正向频偏波束向负角度弯曲”。5.2 距离扫描范围超过周期导致栅瓣误判为目标现象方向图在某个非目标距离出现与主瓣等高的亮带。原因FDA 阵列因子的距离周期性导致主瓣重复周期 c/Δf。解决把扫描范围控制在 0 到 c/Δf 之内。如果必须扫描更远距离就在结果图里标注周期边界提醒自己那是重复主瓣而不是真实目标。5.3 MATLAB 维度不匹配导向矢量被写成矩阵现象MVDR 权重计算时报错“维度不匹配”或得到非 1×N 权重向量。原因导向矢量函数里(0:N-1).和theta/r的维度没有对齐最常见的是把 r 定义为向量而不是标量。解决强制要求导向矢量函数第一行把输入标量化——在函数入口加theta theta(1); r r(1);。养成这个习惯后十次有八次维度报错会消失。5.4 快拍数过少导致协方差矩阵奇异现象MVDR 输出方向图出现剧烈抖动主瓣偏离预期位置。原因快拍数 L N 时采样协方差矩阵秩亏求逆结果被噪声主导。解决至少保证 L 2N通常 L 500N 16 足够同时加对角加载系数 1e-6 到 1e-3 之间。从经验来看对角加载加得太大会导致自适应能力下降零陷变浅加得太小低快拍时算法不稳定。5.5 副瓣电平异常高阵元间距违反半波长约束现象方向图除主瓣外在 ±30° 附近出现非对称高副瓣。原因阵元间距 d 超过半波长导致空间栅瓣混入。解决用lambda c/f0; d 0.5*lambda;严格计算不要手输近似值。我遇到过一次把 d 手输成 0.02 m大约 0.67 倍波长副瓣瞬间抬高了 8 dB当时以为是算法问题排查了半天才发现是采样间距写错。6. 进阶技巧把 FDA 波束做成多目标搜索模式到这里你已经能跑通 FDA 发射和接收的基本流程但资源里最有价值的部分其实是如何把 FDA 的“距离相关”能力用到多目标搜索上。传统相控阵要覆盖一个距离-角度区域需要多帧扫描FDA 因为单帧波束本身就覆盖一条距离-角度耦合轨迹所以理论上单帧就能对特定耦合轨迹上的目标进行搜索。实操方法是设定一个期望角度范围对每个角度计算对应的耦合距离 r(θ) r0 R0·sinθR0 为参考距离然后在这个轨迹上设置密集的波束形成“角度-距离联合搜索线”。% 沿FDA天然波束轨迹设置搜索波束 theta_traj linspace(-30, 30, 50) * pi/180; r_ref 15e3; % 参考距离 r_traj r_ref * (1 sin(theta_traj)); % 耦合轨迹近似 search_power zeros(1, length(theta_traj)); for idx 1:length(theta_traj) a fda_steering_vector(theta_traj(idx), r_traj(idx), f0, delta_f, d, N); search_power(idx) abs(w * a)^2; % w 来自上一节的MVDR权重 end % 峰值位置对应目标的角度随后可反查距离 [~, peak_idx] max(search_power); target_theta_est theta_traj(peak_idx) * 180/pi; target_r_est r_traj(peak_idx); fprintf(估计目标: 角度 %.2f°, 距离 %.2f km\n, ... target_theta_est, target_r_est/1e3);这段代码的关键在于r_traj与theta_traj的耦合关系必须与 FDA 实际的波束弯曲方向一致——如果频偏为正导致波束负角度弯曲上面公式的号就要改成-。之前的避坑要点在这里直接生效。我最初做这个多目标搜索时因为忘了调符号估计值在角度上偏移了大约 8°后来统一用一个函数生成轨迹r_traj r_ref * (1 sign(delta_f) * sin(theta_traj))再也没出过偏移问题。多目标场景下可以同时保留两组搜索波束正频偏和负频偏分别搜索正角度区域和负角度区域输出功率叠加后取两组的局部峰值。要注意的是两组波束在距离维的重复周期不同为了避免误检要设置一个距离窗只保留第一周期内的峰值。从那以后我每次做 FDA 仿真都强制走一遍“先定义频偏符号 → 再画轨迹 → 最后搜索”的流程搜索角度结果基本一次到位不会再出现方向对但距离偏的尴尬情况。这份资源对需要做频控阵预研或算法对比的人来说至少省了半个月的前期摸索希望帮到你。本文还有配套的精品资源点击获取