简介：一组围绕802.11标准OFDM技术的MATLAB仿真与FPGA实现参考资源，适合通信专业学生、无线通信算法工程师及FPGA开发者用于学习多载波调制与物理层设计。压缩包共54个文件，约28KB，以50个.m脚本为主，涵盖发送端的符号映射、交织、卷积编码、IFFT、加循环前缀，以及接收端的定时同步、频率同步、信道估计、解交织、Viterbi译码等功能模块，另有MAT数据文件和说明文档，结构清晰，便于按模块调用与二次开发。已有180人学习该资源。借助这份代码，读者既能通过MATLAB灵活调整子载波数、循环前缀等参数观察系统性能，也能为后续将算法移植到FPGA提供关键实现参考，是理解OFDM收发机原理和硬件化落地的实用素材。
1. 一套能跑通 802.11a/g PHY 的 OFDM 基带链路，值得拆开看
很多人拿到 802.11 OFDM 仿真代码，第一反应是“直接跑通看 BER”。但真正打开这套ofdm_program包后会发现，它不像教材里那种只有 10 个文件的演示 Demo，而是一条完整的发射—信道—接收基带链路：从卷积编码、交织、QPSK/QAM 调制、导频插入、IFFT 加循环前缀，到收端的数据包检测、粗/细时间同步、频偏估计与补偿、信道估计、相位跟踪、解交织和维特比译码，全都有独立的.m文件对应。也就是说，你可以逐级打断链路，观察每个处理块对最终误码率的影响，这对于理解 OFDM 在 802.11 下的工程约束（而非理论公式）非常有价值。这套代码适合三类人：一是刚接触 802.11 PHY 的工程师，需要把“OFDM 原理”和“实际收发处理流程”对应起来；二是在做 FPGA 算法原型验证的人，MATLAB 端行为模型可以作为浮点参考；三是做 Wi-Fi 物理层测试的人，可以用它快速构造特定信道条件下的收发样本。下文会按收发链路顺序，把关键模块的文件映射、参数含义和容易踩的坑逐一拆开。
2. OFDM 收发链路的核心模块映射：从文件清单看清 802.11 处理顺序
2.1 发射端信号链：编码—交织—调制—IFFT—加前缀
802.11a/g 的发射端处理顺序是固定套路：数据比特先经过卷积编码（tx_conv_encoder.m），然后按编码率做打孔（tx_puncture.m，配合get_punc_params.m选择 1/2、2/3、3/4 等码率），打孔后的比特流进入交织器（tx_interleaver.m）。接着根据调制阶数做星座映射（tx_modulate.m，BPSK/QPSK/16QAM/64QAM 由get_bits_per_symbol.m决定），在频域插入导频（tx_add_pilot_syms.m），再补零到 64 点 IFFT 的输入长度（tx_make_int_num_ofdm_syms.m负责凑整符号数，tx_freqd_to_timed.m做 IFFT）。时域符号最后通过tx_add_cyclic_prefix.m加循环前缀。值得注意的是tx_radon_hurwitz.m，这个文件从名字看是发射分集里的旋转矩阵操作，属于可选的空间映射，默认仿真可能不用它。
这个流程里最容易忽略的是tx_power_amplifier.m和phase_noise_model.m。有些教材级仿真根本不建模 PA 非线性和相位噪声，但 802.11a OFDM 对这两项非常敏感：PA 压缩会带来子载波间的互调干扰，相位噪声则会破坏子载波正交性。这套代码把这两个模块独立出来，意味着你可以单独开关它们，观察对星座图和 EVM 的影响。
2.2 接收端信号链：从时域波形到信息比特
接收端对应的文件要更多一些，因为它还包含同步和信道补偿环节。rx_find_packet_edge.m负责检测数据包起始位置，rx_fine_time_sync.m做细定时，rx_timed_to_freqd.m做 FFT。频域之后，rx_estimate_channel.m根据长训练序列估计信道响应，rx_pilot_phase_est.m跟踪残余相位误差，rx_phase_tracker.m做相位补偿。数据部分经过rx_demodulate.m得到软比特或硬比特，再送入rx_deinterleave.m、rx_depuncture.m和rx_viterbi_decode.m。rx_diversity_proc.m和rx_mr_combiner.m对应接收分集合并，属于可选项。rx_radon_hurwitz.m与发射端的旋转矩阵对应，用于恢复 Alamouti 类编码的信号。
这套接收链路的核心特点是：同步模块不是摆设。rx_frequency_sync.m（载波频偏估计）、rx_find_packet_edge.m（包检测）、rx_fine_time_sync.m（精定时）三个模块对 OFDM 性能影响最大。比如频偏如果不校正，FFT 之后子载波间干扰会直接抬高误码平台，这在后面第 4 章会详细展开。
2.3 文件—功能对照速查表
| 文件 | 功能 | 输入/输出 | 所属阶段 |
|---|---|---|---|
tx_conv_encoder.m | 卷积编码（约束长度 7） | 比特流 → 编码比特 | 发射 |
tx_puncture.m | 打孔，匹配码率 | 编码比特 → 打孔后比特 | 发射 |
tx_interleaver.m | 交织，抗突发错误 | 打孔比特 → 交织比特 | 发射 |
tx_modulate.m | 星座映射 | 比特 → 复符号 | 发射 |
tx_add_pilot_syms.m | 插入导频子载波 | 数据符号 → 带导频符号 | 发射 |
tx_freqd_to_timed.m | IFFT，频域转时域 | 频域符号 → 时域样点 | 发射 |
tx_add_cyclic_prefix.m | 加循环前缀 | 时域符号 → 带 CP 符号 | 发射 |
rx_find_packet_edge.m | 包检测 | 时域采样 → 包起始索引 | 接收 |
rx_frequency_sync.m | 载波频偏估计与补偿 | 时域采样 → 去频偏数据 | 接收 |
rx_fine_time_sync.m | 细定时同步 | 去频偏数据 → 对齐符号 | 接收 |
rx_timed_to_freqd.m | FFT，时域转频域 | 时域符号 → 频域符号 | 接收 |
rx_estimate_channel.m | 信道估计 | 长训练序列 → 信道响应 | 接收 |
rx_pilot_phase_est.m | 导频相位估计 | 频域符号 → 相位误差估计 | 接收 |
rx_phase_tracker.m | 残余相位补偿 | 符号序列 → 相位修正数据 | 接收 |
rx_demodulate.m | 星座解映射 | 复符号 → 软比特 | 接收 |
rx_deinterleave.m | 解交织 | 软比特 → 恢复顺序 | 接收 |
rx_depuncture.m | 解打孔（插零） | 软比特 → 恢复编码率 | 接收 |
rx_viterbi_decode.m | 维特比译码 | 软比特 → 信息比特 | 接收 |
这个映射关系是理解整套代码的索引，后续调试任何一个环节，都能直接定位到对应文件。
2.4 一条链路串起来：单包收发的最小流程
以下代码等价于single_packet.m的核心流程。它不做多包统计，专门用来单包调试，观察每个中间变量的维度变化：
% single_packet.m 核心流程（简化） cfg = set_sim_consts(); % 读取仿真参数结构体 bits = randi([0 1], cfg.psdu_bits, 1); % 生成随机 PSDU tx = transmitter(); tx.bits = bits; tx.mod = cfg.modulation; % 'QPSK', '16QAM' 等 tx.code_rate = cfg.code_rate; % 1/2, 2/3, 3/4 tx.n_fft = 64; % 802.11a 固定 64 点 tx.cp_len = 16; % 循环前缀长度（采样点） waveform = tx.tx_chain(); % 经过编码+调制+IFFT+加CP rx = receiver(); rx.waveform = waveform; rx.n_fft = tx.n_fft; rx.cp_len = tx.cp_len; sig = rx.rx_chain(); % 经过同步+FFT+信道估计+译码 errs = sum(bits ~= sig.bits_est); % 统计误码 fprintf('bit errors = %d / %d\n', errs, length(bits));这段代码说明的是收发链路的主干调用关系：tx_chain()内部顺序调用 2.1 小节里的各发射模块，rx_chain()则对应 2.2 小节的接收模块。后面章节分析每个模块时，都建议回到这一段来验证改动效果。需要注意，transmitter()和receiver()是这套包对tx_*.m和rx_*.m函数封装出来的顶层句柄，实际变量名可能不同，但主链路由ui_start_sim.m到runsim.m，再到收发两个核心脚本文件，顺序是一贯的。
3. 关键模块的 MATLAB 实现与参数设置：调制、频偏估计和维特比译码怎么调
3.1 星座映射与解映射：tx_modulate.m / rx_qpsk_demod.m 的参数
tx_modulate.m的核心是根据get_bits_per_symbol.m返回的每符号比特数做查表映射。以 QPSK 为例，每 2 bit 映射一个符号。实际仿真里最影响性能的设置是归一化因子：802.11a 标准要求所有调制阶数的平均符号能量为 1，这样发射功率不随调制方式变化。若你发现 64QAM 链路整体 BER 比理论值差 1 dB 左右，先检查调制函数的归一化常数是否正确。
% QPSK 映射示例（tx_modulate.m 内核心逻辑） bits = reshape(bits, 2, [])'; sym_table = [ 1+1i, -1+1i, -1-1i, 1-1i ] / sqrt(2); % 归一化能量 idx = bi2de(bits, 'left-msb') + 1; symbols = sym_table(idx).';归一化因子sqrt(2)就是让 QPSK 四个星座点的平均能量为 1。如果去掉它，发射功率会变，接收端信噪比计算就会整体偏差。对于 16QAM 是sqrt(10)，64QAM 是sqrt(42)。这个细节在把 MATLAB 浮点模型转到 FPGA 定点实现时尤其重要，因为定点化过程经常因为归一化因子的舍入方式不同，导致定标不一致。
解映射侧，rx_qpsk_demod.m这类函数如果实现的是硬判决，直接找最近星座点即可；如果做软解调，要输出对数似然比。对于后续级联维特比译码，建议至少用 3 bit 量化的软信息，硬判决会让卷积编码的编码增益损失 2 dB 以上。
3.2 载波频偏估计与补偿：rx_frequency_sync.m
OFDM 对载波频偏的容忍度很低，子载波间隔越小越敏感。802.11a 的子载波间隔是 312.5 kHz，频偏是子载波间隔的百分之几时就会显著抬高误码率。rx_frequency_sync.m的常见做法是利用短训练序列的周期性做延迟相关估计：
% 基于短训练序列重复结构的 CFO 估计 L = 16; % 短训练序列重复周期（采样点） N = length(rx_signal) / L; % 重复段数 % 取相邻两段做共轭相关 r1 = rx_signal(1:L); r2 = rx_signal(L+1:2*L); corr = sum(conj(r1) .* r2); % 共轭相关累加 % 频偏估计值（弧度/采样点） delta_phi = angle(corr); cfo_norm = delta_phi / (2 * pi * L); % 归一化到采样率 % 补偿：时域乘指数 t = (0:length(rx_signal)-1).'; rx_corrected = rx_signal .* exp(-1i * 2 * pi * cfo_norm * t / length(rx_signal));这里的核心是angle(corr)的相位差只落在[-pi, pi]区间，所以能估计的频偏范围受限于延迟相关长度 L。L 越大，估计越精细，但可估计范围越小。802.11a 短训练序列用 16 点延迟相关，能覆盖的最大频偏约为采样率的 1/32。如果实际频偏超过这个范围，就需要先做粗频偏校正（通常用短训练序列），再做细频偏校正（用长训练序列）。调试时如果发现频偏估计值总是跳变或者锁定在错误值，先检查接收波形头部的短训练序列是否已经对齐，再看相关窗口是否跨越了两个符号边界。
3.3 卷积编码与维特比译码的参数联动
tx_conv_encoder.m和rx_viterbi_decode.m是误码性能的兜底。802.11a 标准卷积码是约束长度 7、生成多项式 [133, 171]（八进制），这个参数在 MATLAB 中对应poly2trellis(7, [133 171])。rx_viterbi_decode.m的输入如果直接接rx_depuncture.m的输出，注意解打孔是在软比特域完成的——被打掉的比特位置补 0，对应译码器里的“无信息”状态，而不是补一个很大的负值。
打孔参数由get_punc_params.m返回，不同码率的打孔模式不同：
| 码率 | 打孔模式 | 说明 |
|---|---|---|
| 1/2 | 无打孔 | 每个编码比特都发送 |
| 2/3 | 每 4 个编码比特删 1 个 | 对应 P(1,0,1,1) 模式 |
| 3/4 | 每 6 个编码比特删 2 个 | 对应 P(1,0,1,0,1,1) 模式 |
调试维特比译码时有一个常见误区：直接在硬判决比特上做解打孔。软比特域打孔位置补零，等效于该位置对路径度量的贡献为 0，这是正确的；硬判决域则无法表达“不确定”状态，效果会差很多。另外，$1/2$ 码率到 $3/4$ 码率虽然吞吐率提升 50%，但编码增益损失约 0.5 dB，在低信噪比场景下要谨慎选择。
4. 同步、信道估计与相位跟踪的实施要点：OFDM 系统最脆弱的三个环节
4.1 包检测与粗定时：rx_find_packet_edge.m 的门限策略
rx_find_packet_edge.m做的事情是判断“数据包从哪个样本开始”。工程实现常用延迟自相关 + 能量归一化：用短训练序列的周期性，计算滑动窗口内相邻两段波形的相关值，再用该时刻的能量做归一化。门限不能设得太低，否则噪声尖峰就会误触发；太高则会漏检，把包开头的一部分丢掉，导致 FFT 窗口错位，子载波间干扰（ICI）增大。
具体实现时，通常计算一个滑动相关度量：
L = 16; % 短训练序列周期 P = zeros(size(rx)); R = zeros(size(rx)); for n = L+1 : length(rx) P(n) = P(n-1) + conj(rx(n)) * rx(n-L) - conj(rx(n-L)) * rx(n-2*L); R(n) = R(n-1) + abs(rx(n-L))^2 - abs(rx(n-2*L))^2; m(n) = abs(P(n))^2 / (R(n)^2 + eps); end start_idx = find(m > 0.6, 1, 'first'); % 门限 0.6 是经验值这里的门限值取了 0.6，但这个值不是固定的。实际使用中，若信道 SNR 低于 5 dB，建议降到 0.4；若前端有 AGC 增益波动，还需要用两段滑动窗口做归一化。还要注意帧结构里短训练序列之后是长训练序列，中间有保护间隔，定时点应该取在长训练序列的起点附近，而不是短训练序列的峰上，否则 FFT 窗口会向前偏，造成符号间干扰。
4.2 信道估计的可观测性问题：rx_estimate_channel.m 和导频子载波
802.11a 的长训练序列（LTF）在频域是已知序列，接收端在 FFT 之后用收到的 LTF 除以已知 LTF，就能得到每个子载波上的信道响应估计值。代码如下：
% 频域信道估计（rx_estimate_channel.m 核心） H_est = rx_ltf_fft ./ tx_ltf_fft; % 逐子载波复数除法这个估计只有 52 个有效子载波有值，其余子载波置零。估计结果可以直接用于均衡，但若信道是频率选择性的，单个符号上的估计受噪声影响大。常见改进是跨多个符号做平均，或者用插值方法从导频位置恢复数据位置的信道响应。rx_pilot_phase_est.m就是在数据符号内部用 4 个导频子载波估计残余相位误差的。残余相位的来源是发射机和接收机的晶振频率不完全一致，导致每个符号上有累积的公共相位旋转。这个量对高阶层调制（64QAM）特别致命，因为 64QAM 星座点之间的距离小，相位旋转几个度就可能跨到相邻星座点。
4.3 从 MATLAB 浮点模型转向 FPGA 定点时的量化规则
包里的 MATLAB 代码全部是浮点实现的，但既然是面向 FPGA 的 802.11 链路，定点化就是必经之路。从工程实践看，做 FPGA 实现前先在 MATLAB 里做一次“伪定点”仿真能省大量调试时间：把每个关键节点的数据宽度定下来，提前通过仿真确定小数的位宽，而不是到 RTL 阶段再试。较常碰到的量化策略如下：
% 浮点转定点示例（用 fi 对象模拟） % 假设 IFFT 输入数据位宽 16 bit，其中 1 bit 符号，1 bit 整数，14 bit 小数 word_len = 16; frac_len = 14; data_fi = fi(data, true, word_len, frac_len);量化影响最大的三个节点分别是：IFFT/FFT 的蝶形运算中间节点、信道估计的除法结果、维特比译码的软比特量化。FFT 的蝶形运算每级都可能增加位宽，若中间结果不舍入而只在最后截位，资源开销会变大；若每级都截位，噪声会累积。常见做法是每级增加 1 bit 防止溢出，最后统一切到目标位宽。信道估计的除法在 FPGA 上是重点资源消耗项，较常改用查表方式计算倒数，然后用乘法替代除法。软比特到维特比译码的接口一般用 3~4 bit 有符号数，符号位和幅度位要固定好，否则译码器性能衰减明显。
4.4 常见失败模式与排查对照表
| 现象 | 可能原因 | 排查手段 |
|---|---|---|
| BER 在 20 dB 后不再下降 | 频偏未完全补偿或相位跟踪失效 | 查看星座图是否有旋转残差 |
| 包检测频繁误触发 | 门限设置过低 | 打印滑动相关度量值观察底噪 |
| FFT 窗口错位导致地板效应 | 粗定时点偏移到保护间隔内 | 检查符号起点定位逻辑 |
| 64QAM 星座图发散 | AGC 增益未收敛或相位噪声未建模 | 单独开关 phase_noise_model 观察 |
| 维特比输出全 0 | 解打孔位置补错值 | 对比编码前后长度，逐段打印 |
相位噪声的建模在phase_noise_model.m中是用一个随机游走过程叠加到相位上，其功率谱密度由 linewidth 参数控制。在实际 Wi-Fi 芯片里，锁相环（PLL）的相位噪声是性能上限的主要来源之一，这点在 FPGA 原型验证时容易被忽略，因为 FPGA 上的时钟抖动往往比 ASIC 锁相环要差得多。若在 FPGA 上做同频段回环测试发现 EVM 偏高，优先检查时钟抖动和电源噪声。
5. 跑通仿真与调试技巧：从 runsim.m 入口到单包定位
5.1 仿真入口与参数文件组织
runsim.m和ui_start_sim.m是整套代码的入口。仿真参数不散落在各模块里，而是集中在set_sim_consts.m中返回一个结构体。常用参数包括：调制方式、码率、信道模型（channel.m+get_channel_ir.m）、信噪比范围、是否使能频偏/相位噪声/PA 非线性。建议第一次跑通时按以下配置：
cfg.modulation = 'QPSK'; cfg.code_rate = 1/2; cfg.snr_db = [0:2:20]; cfg.freq_offset_hz = 0; % 先关掉频偏 cfg.channel_type = 'awgn'; % 先不用多径信道 cfg.enable_pa = false; % 关闭 PA 非线性 cfg.enable_phase_noise = false;先用最简单的配置确认整条链路无逻辑错误，再逐步打开频偏、多径和相位噪声。一条经验：如果基本配置下 BER 不能和理论曲线吻合，先不要怀疑信道模型或同步模块，往往是某个模块的归一化因子或索引对不齐导致的系统性错误。
5.2 单包调试的操作步骤
单包调试是定位问题最高效的手段。在single_packet.m模式下，把发射端每个中间变量保存下来，接收端每个中间变量再对比一遍。具体做法：
- 在发射链路末尾把
tx_bits、tx_symbols、tx_ifft、tx_with_cp依次保存。 - 在接收端同步完成后，先检查 FFT 输入数据的起点是否落在期望的符号边界上。
- 打印信道估计结果
H_est的幅值和相位，观察是否平滑——如果相邻子载波间相位乱跳，说明定时偏了。
建议在接收端各模块加fprintf输出维度信息，这样能快速定位是哪一级的数据维度与预期不符。尤其在 FFT 前后，时域长度是 $64 + 16 = 80$ 个采样点（含循环前缀），在符号边界对齐后应该是每 80 点一个符号，其中前 16 点是 CP，后 64 点进 FFT。维度对不上的情况基本都出在奇数/偶数索引的偏移上，这是 MATLAB 与硬件实现最常见的差异之一。
5.3 用误码率曲线做端到端验证
链路调通后，跑一组 SNR 扫描画 BER 曲线，与理论值对比，是验证实现是否正确的最终手段。注意事项是，若信道里包含频偏和相位噪声，BER 曲线会在高 SNR 处出现平台（error floor），这是正常现象，不是代码 bug。如果直接对比纯 AWGN 理论曲线，应把所有非理想因素先关掉。多径信道下不要直接和 AWGN 理论曲线比，建议先在同一信道模型下对比不同调制方式和码率的相对增益，看趋势是否符合预期。
最后谈一个实用的排错技巧：在做 OFDM 系统调试时，把接收端rx_timed_to_freqd.m输出的频域数据直接打印成复数数组，和发射端tx_freqd_to_timed.m的输入做逐点对比。如果发现高频子载波上有固定规律的错误图案，优先检查子载波映射表——802.11a 中 64 个 FFT 子载波只有 52 个用于数据，6 个是零填充，4 个是导频，数据子载波的索引映射不是连续排列，中间有 DC 空子载波。这个映射表一旦写错，表现出的现象就是高频子载波数据错乱，而低频部分完全正常，容易误判成信道问题而非映射表问题。对照标准中的子载波编号逐一检查，能节省大量时间。
本文还有配套的精品资源，点击获取
