简介无线通信物理层波形设计是6G技术演进的核心方向之一。传统OFDM在高铁、低轨卫星等高速移动场景下因多普勒频移导致子载波间干扰加剧性能急剧下降。AFDM仿射频率分复用通过离散仿射傅里叶变换DAFT将信号映射到时延-多普勒域从根本上将多普勒效应纳入信号设计成为6G候选波形的有力竞争者。本文从OFDM局限性引入DAFT变换原理拆解了完整的MATLAB仿真链路发射机调制、时变多径信道建模、接收机均衡与误码率统计。通过BER/PEP仿真与参数调优可直观验证AFDM在高移动性场景下的性能增益为无线通信从业者及研究生提供从公式到可运行代码的工程参考。1. AFDM波形是什么高速移动场景下替代OFDM的6G候选做无线通信物理层的人这几年应该都听过AFDM全称是仿射频率分复用Affine Frequency Division Multiplexing。它本质上是把OFDM的傅里叶变换基底换成了离散仿射傅里叶变换DAFT在时延-多普勒域里重新排布信号能量从而在高移动性场景下比OFDM更抗多普勒干扰。如果你正在做6G波形预研或者被高铁、低轨卫星这类高速场景下的子载波间干扰ICI折磨过这套AFDM波形理论解析源码就是拿来拆的它给出了从发射机调制、时变多径信道到接收机解调的完整链路并带误码率BER和成对错误概率PEP仿真脚本。适合无线通信方向的从业者、研究生以及任何想快速把AFDM从论文公式落成可运行代码的人。2. 从OFDM到AFDM离散仿射傅里叶变换到底改了哪里2.1 为什么OFDM在时变信道下会翻车OFDM之所以成为4G和5G的物理层基石靠的是把宽带信道切成一组窄带子载波每个子载波经历近似平坦衰落。这个设计隐含了一个假设在一个OFDM符号持续时间内信道是时不变的。但在高铁每小时350公里、或者低轨卫星以7.6公里每秒划过地面的时候多普勒频移可以轻松达到千赫兹量级子载波之间的正交性被破坏能量从一个子载波泄漏到相邻子载波这就是ICI。传统的解决办法是缩短符号时长、增大子载波间隔但代价是循环前缀CP开销占比上升频谱效率掉得很难看。AFDM的思路不是去“对抗”多普勒而是把多普勒放进信号设计里。它把时延维度和多普勒维度同时映射到一组带调频率的仿射基底上让信号在时延-多普勒域里呈现稀疏且可分离的结构接收端用一个矩阵变换就能把时延和多普勒同时估计出来。2.2 离散仿射傅里叶变换DFT加上两个chirp参数离散仿射傅里叶变换DAFT是AFDM的核心运算单元。DFT用的是复指数基底 e^{-j2\pi kn/N}DAFT在此基础上额外乘了两个二次相位项一个作用在时域采样点 n 上一个作用在变换域索引 k 上分别由参数 c1 和 c2 控制。这两个参数就是chirp调频率它们决定了基底在时间-频率平面上的倾斜方向。下面是构造 N 点 DAFT 矩阵的 MATLAB 实现这段代码也是这套源码里最底层的基础函数function F daft_matrix(N, c1, c2) % 构造 N 点离散仿射傅里叶变换矩阵 % F(k1, n1) exp(-j*2*pi*(k*n/N c1*n^2 c2*k^2)) / sqrt(N) % c1: 时域 chirp 调频率 % c2: 变换域 chirp 调频率 [n, k] meshgrid(0:N-1, 0:N-1); F exp(-1i * 2 * pi * (k .* n ./ N c1 * n.^2 c2 * k.^2)) / sqrt(N); end逻辑说明meshgrid 生成两组网格索引k .* n ./ N 是传统 DFT 项c1 * n.^2 相当于给时域信号叠了一个二次相位包络c2 * k.^2 相当于给变换域加相位旋转。整个矩阵除以 sqrt(N) 是为了保持能量归一化保证 DAFT 是酉变换。当 c1 c2 0 时这个矩阵退化为标准 DFT 矩阵——这一点后面排错时非常关键。参数说明N 一般取子载波数量通常是 128 到 1024 之间的2的幂c1 和 c2 的取值决定了基底对时延-多普勒平面的覆盖能力典型取值是 N 的整数倍分之一或者直接用 0 到 1 之间的实数搜索得到。c1 负责时延分辨c2 负责多普勒分辨两者需要联合设计。2.3 AFDM信号模型从矩阵到发射链路有了DAFT矩阵AFDM的收发模型就清晰了。设 N 为子载波数M 为时隙内符号块长度实际就是一次处理的数据块发射端把 M×N 个QPSK或QAM符号 d[m, n] 排成 M×N 矩阵先沿着 N 维度做逆DAFTIDAFT再加CP得到时域发送信号。接收端做正DAFT把时域信号映射回DD域做均衡和判决。% 发射机核心QPSK符号 - IDAFT - 加循环前缀 function x afdm_tx(bits, M, N, c1, c2, cp_len) % bits: 0/1比特流 % M: 数据块长度每块处理 M 个符号 % N: DAFT 变换点数即子载波数 % c1, c2: chirp 调频率 % cp_len: 循环前缀长度 sym qpsk_mod(bits, M, N); % QPSK 映射输出 M×N 复符号矩阵 F daft_matrix(N, c1, c2); % 正变换矩阵 x_tf sym * F.; % 对每行做 IDAFT用 F 的转置 x_cp [x_tf(:, end-cp_len1:end), x_tf]; % 加循环前缀 x x_cp(:).; % 串行化输出 end逻辑说明sym * F. 是对符号矩阵的每一行做时域变换。用 F. 而不是 F 是因为这里只需要共轭对称的逆变换不取共轭保持与正变换方向一致。加 CP 时把每行末尾 cp_len 个样本复制到行首和 OFDM 的做法一致为的是把线性卷积变成循环卷积让接收端能做频域/DD域单抽头均衡。参数说明M 是每帧处理的符号块数影响处理时延和频谱效率cp_len 必须大于信道的最大时延扩展典型取 N/4 或 N/8。3. 源码拆解发射机、信道与接收机怎么协同工作3.1 工程文件结构与运行入口这套源码的目录结构按功能拆得比较清楚拿到手先看入口脚本顺着调用关系往下追很快就能理清。建议用 MATLAB R2020a 及以上版本打开不需要额外工具箱纯脚本加几个自写函数。文件名职责关键接口main_afdm_ber.m主仿真入口跑BER扫描配置结构体 cfgafdm_tx.m / afdm_rx.m发射/接收链路(bits, M, N, c1, c2, cp_len)daft_matrix.m生成DAFT变换矩阵(N, c1, c2)channel_twv.m时变多径信道模型(x, fd, delays, gains)qpsk_mod.m / qpsk_demod.mQPSK调制解调(bits / syms)utility_plot.m星座图/DD域可视化(data, fig_idx)3.2 发射机从比特流到DAFT域信号主脚本里按顺序执行生成随机比特、QPSK调制、AFDM发射、过信道、接收、解调、对比误码率。发射模块是最值得读的部分它决定了信号在DAFT域的形状。% main_afdm_ber.m 中的配置与循环节选 N 128; % 子载波数 M 64; % 每帧符号块数 c1 0.1 / N; % 时域 chirp 系数 c2 0.1 / N; % 变换域 chirp 系数 cp_len 32; % 循环前缀长度 snr_list 0:2:20; % 信噪比扫描范围 dB for idx 1:length(snr_list) bits randi([0 1], M * N * 2, 1); % QPSK: 每符号2比特 x afdm_tx(bits, M, N, c1, c2, cp_len); [y, h] channel_twv(x, 500, [0 3 7], [1 0.5 0.3]); y_noisy awgn(y, snr_list(idx), measured); bits_hat afdm_rx(y_noisy, h, M, N, c1, c2, cp_len); ber(idx) mean(bits ~ bits_hat); end逻辑说明每帧的符号总数是 M×NQPSK再乘以2得到比特数。信道函数返回经过时变多径信道后的信号 y 和信道冲激响应 h接收端把 h 一起传进去做均衡。awgn 的 measured 选项会让噪声功率按信号实际功率计算避免因发射功率归一化误差导致信噪比偏移。参数说明c1 和 c2 这里取了 0.1/N是让调频率周期远大于符号块长度保证信号在DD域不扩散多径延时用采样点表示[0 3 7] 对应三条径分别落后0、3、7个采样点幅度依次衰减。3.3 信道时变多径怎么建模才接近真实时变信道是AFDM仿真里最容易写错的部分。常见的错误是把信道写成时不变抽头延迟线那样多普勒频移完全体现不出来。正确的做法是给每条径叠加一个随采样点变化的相位旋转。function [y, h] channel_twv(x, fd, delays, gains) % 时变多径信道每径带独立多普勒频移 % x: 时域发射信号串行 % fd: 最大多普勒频移 Hz % delays: 各径时延采样点 % gains: 各径复数增益 fs 10e6; % 采样率 10 MHz t (0:length(x)-1) / fs; % 时间轴 y zeros(size(x)); h zeros(length(delays), length(x)); for p 1:length(delays) phase exp(1j * 2 * pi * fd * t); % 多普勒相位旋转 h(p, :) gains(p) .* phase; if delays(p) 0 y y [zeros(1, delays(p)), gains(p) .* phase(1:end-delays(p)) .* x(1:end-delays(p))]; else y y gains(p) .* phase .* x; end end end逻辑说明每个抽头延迟线分支乘以一个随时间线性增长的复数相位等效于在频域搬移该径的频谱这就是多普勒频移的实现。最直接的验证方法是把 fd 设成 0信道退化为普通多径把 fd 设到 1kHz 以上接收端星座图会开始旋转。参数说明fd 是最大多普勒频移典型值按 v*f_c/c 计算比如 300 km/h、3.5GHz 载波下大约 1000 Hzdelays 和 gains 的长度一致分别对应径的时延和平均功率。3.4 接收机DAFT域均衡与符号判决接收端做的是发射端的逆过程去掉CP、做正DAFT、在DD域做均衡、QPSK解调。均衡这一步用最小二乘LS还是最小均方误差MMSE会直接影响误码率曲线源码里给的是LS版本简单直接。function bits_hat afdm_rx(y, h, M, N, c1, c2, cp_len) % AFDM 接收机去CP - DAFT - LS均衡 - QPSK解调 L M cp_len; y_mat reshape(y, L, M).; % 串行转矩阵每行是一个带CP的块 y_mat y_mat(:, cp_len1:end); % 去掉循环前缀 F daft_matrix(N, c1, c2); Y y_mat * F.; % 正DAFT变换到DD域 % 基于导频估计的信道频响做LS均衡 H_est fft(h(1,:), N); % 以第一条径做参考信道估计 X_eq Y ./ max(H_est, 1e-6); % 频域单抽头均衡 bits_hat qpsk_demod(X_eq(:)); end逻辑说明这里用了简化处理把信道等效成DD域对角矩阵做均衡。严格做法是构造完整的时域信道矩阵再投影到DAFT域但仿真里LS单抽头在低阶调制下已经能反映误码率趋势。H_est 用的是第一条径的FFT代表参考信道响应除以它相当于频域均衡。参数说明max(H_est, 1e-6) 是为了避免除零实际项目中如果导频密度足够应该用插值后的完整信道估计代替 H_est否则在高多普勒下均衡会残留较大误差。4. 参数怎么调从仿真配置到系统级权衡4.1 可调参数与默认值上手跑通主脚本之后下一步就是改参数。参数之间是强耦合的调一个不改另一个仿真的结论就可能失真。先把核心参数列出来参数符号默认值调节范围影响子载波数N12864 ~ 1024DAFT粒度越大频域分辨率越高符号块数M6416 ~ 256块越长抗多普勒越好但时延变大时域chirpc10.1/N0 ~ 1/N控制时延域的倾斜程度变换域chirpc20.1/N0 ~ 1/N控制多普勒域的倾斜程度循环前缀cp_len32N/16 ~ N/4必须大于最大时延扩展多普勒频移fd500 Hz0 ~ 2000 Hz超过阈值误码率急剧上升这些参数没有一组“万能组合”需要根据目标场景折中。低速场景下c1、c2设小一点接近OFDM高速场景下加大c1、c2让基底更好地对齐多普勒方向。4.2 调频率的选择逻辑与搜索策略c1 和 c2 是AFDM区别于其他波形的核心自由度它们的作用相当于在时延-多普勒平面里“旋转”基底。选得好时延和多普勒会自然分离选得不好信号在DD域变成一坨接收端均衡器怎么调都救不回来。常见的做法是网格搜索加误码率验证。因为 c1、c2 与 N 的乘积有周期性搜索范围限定在 [0, 1) 即可。先做粗搜步长 0.05/N找到BER最低的区间后再以 0.005/N 的步长细搜。这个过程大概要跑几十轮仿真建议把 BER 扫描的 SNR 固定在 10 dB 或 12 dB而不是全信噪比扫描节省时间。% 调频率网格搜索示例固定 SNR10dB扫 c1 和 c2 best_ber Inf; for c1 0:0.02:0.2 for c2 0:0.02:0.2 ber run_single_sim(c1/N, c2/N, 10); if ber best_ber best_ber ber; best_c [c1/N, c2/N]; end end end fprintf(最优调频率: c1%.5f, c2%.5f, BER%.4f\n, best_c(1), best_c(2));逻辑说明run_single_sim 封装了完整收发链路返回给定信噪比下的误码率。粗搜阶段不要追求精确最优值只要锁定量级。细搜时可以把c1、c2绑定成同一个值减少搜索维度因为时延和多普勒在实际信道中往往同时存在。4.3 与OFDM对比仿真的控制变量方法如果是为了写报告或论文OFDM作为基准波形是绕不开的对照。最常犯的错误是OFDM和AFDM用不同的CP长度或不同的信道估计方式导致对比不公平。对比时必须保证带宽、子载波数N、CP长度、信道模型、调制方式全部一致。% OFDM基准链路的参数必须与AFDM严格对齐 ofdm_cfg.N 128; % 和AFDM的N一致 ofdm_cfg.cp 32; % 和AFDM的cp_len一致 ofdm_cfg.mod QPSK; % 调制方式一致 % 唯一不同的是OFDM用FFT/IFFTAFDM用DAFT/IDAFT逻辑说明OFDM本身就是 c1c20 时AFDM的特例所以这套源码里把OFDM仿真单独写了一个分支复用信道和均衡模块。对比的结果一般是在低速场景两者BER几乎重合随fd增大OFDM劣化加速AFDM保持平稳。5. 避开这五个坑AFDM仿真的常见翻车点5.1 调频率设成0跑了半天发现和OFDM一模一样现象c1 和 c2 设成 0 后误码率曲线和OFDM完全重合怀疑代码没有生效。原因c1c20 时 DAFT 退化为 DFTAFDM 在数学上就是 OFDM仿真结果重合是正确行为而不是Bug。解决先确认这是预期行为再检查 daft_matrix 里非零调频率下的矩阵元素是否与DFT有明显差异——打印 F 的相位矩阵非零 chirp 下应该能看到对角线上的二次相位条纹。5.2 DAFT矩阵的转置方向用反现象接收端星座图完全散开BER 接近 0.5。原因发射用 F 做 IDAFT接收又用 F 做 DAFT两次用同一个方向等于做了两次逆变换信号在DD域没有正确对齐。解决统一约定 F 是正变换矩阵发射端用 F.转置接收端也用 F.因为复指数项的对称性转置等效于逆变换方向前后保持一致。5.3 功率没归一化信噪比越高误码率越差现象SNR从10dB加到20dBBER不降反升。原因DAFT矩阵如果漏了 sqrt(N) 缩放因子发射信号功率随N变化awgn 用 measured 计算噪声功率时会把功率偏移算进去导致实际SNR偏移。解决检查 daft_matrix 里是否有 1/sqrt(N)以及发射链路里 x 的能量是否约等于符号能量——直接跑一行 sum(abs(x).^2) / length(x)应该接近1。5.4 时不变信道代替时变信道现象fd 从100Hz加到2000HzBER曲线纹丝不动。原因channel_twv 里相位旋转写成了常量或者信道系数在符号块内被平均了。解决把信道响应的实部画出来时变信道下应该是振荡曲线时不变信道是水平直线。多普勒相位必须按采样点逐一更新不能按块更新。5.5 低阶调制掩蔽了多普勒损伤现象用QPSK仿真AFDM优势不明显误码率差距不到一个数量级。原因QPSK的抗噪能力太强很多ICI损伤被噪声容限掩盖。解决至少在16QAM下做一次对比高阶调制对残余ICI更敏感AFDM和OFDM的差距会更直观——如果高阶调制下差距仍然很小大概率是信道多普勒没有打进信号里。6. 拿什么证明这套AFDM仿真没白做三项自查与DD域可视化先做零参数退化测试把 c1、c2 全部设为0BER曲线应该与OFDM基准完全重合如果不重合说明DAFT矩阵的某个方向写反了。再做能量守恒测试发射信号经过完整收发链路后无噪声场景下解调星座点应该和发射符号完全重合误差小于 1e-8 量级误差大了就回去检查矩阵归一化。第二项自查是画DD域稀疏图。取一个无噪声的发射块在接收端做正DAFT后把 M×N 的变换域矩阵取模画成图像。正确的AFDM信号应该在 DD 域表现为少量亮斑而不是整个平面均匀分布亮斑越集中说明调频率对齐了信道的时延-多普勒响应。如果图上糊成一片优先怀疑 c1、c2 的搜索范围没选对。第三项自查是追一次PEP成对错误概率事件。源码里PEP统计脚本会记录出错时发射符号和判决符号的差观察这些差向量的结构如果错误集中在相邻DD网格点说明时延-多普勒分辨率刚好卡在极限如果错误随机分布可能是噪声主导。这一步能帮你判断系统是受限于信道弥散还是受限于信噪比。验证通过之后可以往两个方向扩展一是换成16QAM/64QAM把调制阶数参数提上去看高多普勒下的增益是否保持二是把仿真里的时变信道换成 5G CDL 信道模型加入空间相关性和极化损耗逼近真实链路。我给自己定了个习惯每次拿到新的波形源码第一件事永远是跑零参数退化和能量守恒这两项自查过了才敢去调 chirp 参数不然任何实验结论都可能是代码Bug造成的假象。希望帮到你。本文还有配套的精品资源点击获取