简介：这份MATLAB仿真程序面向通信工程专业学生、研究人员与工程师，围绕OFDM多载波调制系统的收发链路搭建，帮助读者从理论理解走向可运行的工程模拟。包内共26个文件，以19个txt数据与说明、4个dat波形数据、2个m脚本和1个doc流程说明为主，压缩包约125KB，体量轻便，便于按模块查阅和二次修改。已有790人学习，说明其在OFDM入门、课程实验与算法验证场景中具有一定参考价值。程序覆盖符号生成、IFFT复用、循环前缀加入、信道建模、FFT接收、信道估计与均衡、解调判决及误码率评估等环节，读者可据此调整子载波数量、调制方式、循环前缀长度和信道特性，并比较MMSE、ZF等均衡方法的性能差异。对需要理解OFDM收发流程、完成仿真验证或开展算法对比的读者，这份程序包提供了较完整的实践底稿，尤其适合课程设计、毕业设计与科研预研中的仿真复现。

1. 调 OFDM 仿真，最耗时间的不是写 FFT 本身

调 OFDM 仿真，真正消耗时间的地方通常不是把fft那几行写出来，而是把发射、信道、同步、均衡、解调这几段串成一条能稳定跑出正确误码率曲线的链路。中间任何一处索引错位、归一化漏写、fftshift用反，星座图就会散开，而 MATLAB 不会报出任何错误。一份完整的 OFDM 仿真程序 MATLAB 版本，价值正在这里：它把 IFFT/FFT、循环前缀、导频插入、LS 信道估计、频偏补偿这些环节按可复现的顺序固定下来，参数集中在文件头，改完就能重跑。这份资源适合三类人：通信课程设计阶段需要一条能出 BER 曲线的最小闭环的人；做链路级仿真需要可信基线的人；以及已经会写单点代码、但一直卡在“曲线形状不对”的人。下面按我拆这类程序的顺序讲，先讲链路结构，再讲损伤、同步，最后讲怎么验证你跑出来的东西不是错的。

2. 收发链路拆解：从比特到 OFDM 符号的七道工序

2.1 参数表先固定下来，不要边写边改

绝大多数仿真翻车，起因都是参数没有一次定死。建议开局就把采样率、子载波间距、FFT 点数和 CP 长度写成一组自洽的数，后面所有模块只读这几个变量。我常用的一组是 Δf = 15 kHz、N_fft = 64，于是采样率 fs = N_fft × Δf = 960 kHz，符号周期 1/Δf ≈ 66.7 μs，CP 取 16 个采样点即 16.7 μs，开销 25%。这组数不是随便挑的：CP 长度决定了能吸收多大的时延扩展，1/(2π·CP) 量级以内基本安全，粗算就是十几微秒，室内多径和典型信道模型都落在里面。

数据速率可以直接算出来：16QAM 时 48 个子载波 × 4 bit ÷ 83.3 μs ≈ 2.3 Mbps。这个数后面写报告要用，而且能帮你判断“我加的信道编码到底换来了多少开销”。

子载波分配是最容易出错的一步。我习惯用以 0 Hz 为中心的排布，数组第 33 位对应直流，两侧各留 4 个和 3 个保护带：

% ===== 全局参数：只改这里，全链路跟着变 ===== N_fft = 64; % FFT 点数 N_cp = 16; % 循环前缀长度 N_sym = 2000; % 一次仿真跑多少 OFDM 符号 M = 16; % 4=QPSK, 16=16QAM, 64=64QAM SNR_dB = 0:2:20; % 信噪比扫描范围 rng(2024); % 固定种子，保证曲线可复现 % ===== 子载波分配：数组下标 33 对应 0 Hz ===== guard_idx = [1:4, 62:64]; % 保护带，不承载信息 dc_idx = 33; % 直流置空，避开本振泄漏 pilot_idx = [5 12 19 26 40 47 54 61]; % 8 个导频，频域间隔 7 data_idx = setdiff(5:61, [dc_idx pilot_idx]); % 48 个数据子载波 pilot_sym = ones(numel(pilot_idx),1) * (1+1j)/sqrt(2); % 单位功率导频

这里有个隐式约定要记住：data_idx和pilot_idx是数组下标，不是频率序号。等会儿ifftshift把它们搬到负频率时，全靠这个约定保持一致。凡是“星座图对、频谱却左右反了”的问题，八成都出在这。

2.2 发射机：映射、IFFT、加 CP 三段就够

发射端的核心只有一句话：频域放好数据，IFFT 到时域，再复制尾部到头部。但每一步的归一化都得交代清楚，否则后面算 SNR 会莫名其妙偏 1 dB。

function [tx, bits] = ofdm_tx(N_sym, M, N_fft, N_cp, data_idx, pilot_idx, pilot_sym) k = log2(M); N_data = numel(data_idx); bits = randi([0 1], N_data*N_sym*k, 1); % 串行比特流 sym = qammod(bits, M, 'InputType','bit', 'UnitAveragePower', true); sym = reshape(sym, N_data, N_sym); % 每列一个 OFDM 符号 X = zeros(N_fft, N_sym); X(data_idx, :) = sym; % 数据入位 X(pilot_idx, :) = repmat(pilot_sym, 1, N_sym); % 导频入位 % ifftshift 把直流从第 33 位搬到第 1 位，再 IFFT x_t = ifft(ifftshift(X, 1), N_fft, 1) * sqrt(N_fft); % 复制尾部 N_cp 个样点放到符号前，形成循环前缀 x_cp = [x_t(end-N_cp+1:end, :); x_t]; tx = x_cp(:); % 并串转换 end

UnitAveragePower打开后每个星座点平均功率是 1，所以 IFFT 出来的时域平均功率等于有效子载波占比 56/64 = 0.875，不是 1。这不是 bug，但你在算理论曲线时要记着它。乘sqrt(N_fft)的作用是抵消 MATLABifft自带的 1/N 归一化，让时域和频域功率满足帕塞瓦尔关系，少写这个系数，信噪比会整体偏移 18 dB，曲线看着“特别差”，其实只是标注错了。

CP 那两行也别轻视。x_t(end-N_cp+1:end, :)取的是每个符号的最后 16 个样点，整体放到符号前面。写成x_t(1:N_cp,:)也能跑，但那就不是循环前缀，而是重复了符号开头，频域上不再等效为相位斜坡，均衡之后会残留很强的 ICI。

提示：CP 复制错位置时，低 SNR 下几乎看不出来，因为噪声掩盖了 ICI；等到 SNR 超过 15 dB 才出现误码平台。这类“高位平底”的曲线，第一嫌疑就是 CP。

2.3 接收机：去 CP、FFT、均衡、解映射

接收端是发射端的逆运算，顺序不能乱：先去 CP，再 FFT，再取子载波，再除以信道估计，最后解映射。

function [bits_hat, H_est] = ofdm_rx(rx, N_sym, N_fft, N_cp, M, ... data_idx, pilot_idx, pilot_sym) rx = reshape(rx, N_fft+N_cp, N_sym); rx = rx(N_cp+1:end, :); % 去 CP Y = fftshift(fft(rx, N_fft, 1), 1); % 频域，直流回到第 33 位 % LS 估计：导频位置直接用除法，已知发端导频 H_p = Y(pilot_idx, :) ./ repmat(pilot_sym, 1, N_sym); % 频域线性插值到所有子载波 H_est = interp1(pilot_idx, H_p, (1:N_fft)', 'linear', 'extrap'); % 数据子载波上做迫零均衡 X_eq = Y(data_idx, :) ./ H_est(data_idx, :); bits_hat = qamdemod(X_eq(:), M, 'OutputType','bit', 'UnitAveragePower', true); end

维度顺序是整个接收机最容易埋隐患的地方。qammod是把比特流按每log2(M)个一组切分后映射的，reshape按列优先填充，所以第 j 列天然对应第 j 个符号的第 1 到第 N_data 个子载波。接收端X_eq(:)展开时同样是列优先，比特顺序才对得上。如果发射端用了reshape(sym, N_sym, N_data).'这类转置写法，接收端必须同步转置，否则误码率恒等于 0.5，看起来像“完全没解出来”，实际只是比特打乱了。

interp1的'extrap'会让保护带位置的估计值外插出无意义的结果，但保护带上的数据本来就不用，所以不影响误码统计。真正要注意的是导频间隔：8 个导频、间隔 7 个子载波，对应 105 kHz，按相干带宽 1/τ_max 估算，能撑住的时延扩展大约 9～10 μs。如果你把信道模型的最大时延设成 20 μs，插值出来的 H 在频率选择性深的子载波上会严重失真，这时第一反应应该是加密导频，不是怀疑均衡公式。

把这三段拼起来，一个只含 AWGN 的最小闭环就能跑了：

bits_err = 0; bits_total = 0; for s = 1:numel(SNR_dB) [tx, bits] = ofdm_tx(N_sym, M, N_fft, N_cp, data_idx, pilot_idx, pilot_sym); rx = awgn(tx, SNR_dB(s), 'measured'); % 按实测功率加噪 [bits_hat, ~] = ofdm_rx(rx, N_sym, N_fft, N_cp, M, ... data_idx, pilot_idx, pilot_sym); ber(s) = mean(bits ~= bits_hat); end semilogy(SNR_dB, ber, 'o-'); grid on; xlabel('SNR (dB)'); ylabel('BER');

awgn(..., 'measured')是按输入信号实测功率反推噪声功率的，省掉了手算归一化的麻烦，代价是它把 CP 的能量也算进信号功率里了——这意味着你标注的 SNR 比“每信息比特能量”偏乐观约 1 dB。做定性验证够用，出正式曲线时把这 0.97 dB 扣回去。

3. 信道与估计：多径、导频与均衡怎么接进去

3.1 多径瑞利信道：功率归一化必须先做

AWGN 只能验证链路通不通，验证不了 OFDM 的抗多径能力。加多径时最容易错的一步是功率归一化：每条径的平均功率加起来必须等于 1，否则你标注的 SNR 和实际 SNR 对不上，曲线会整体平移好几 dB，而你还以为是自己算法写错了。

function [rx, h_tap] = chan_multipath(tx, fs, delay_s, pdp_dB, snr_dB) % delay_s : 各径时延（秒），第一个必须为 0 % pdp_dB : 各径平均功率（dB），例如 [0 -3 -8 -12] p = 10.^(pdp_dB(:)/10); p = p / sum(p); % 归一化，总功率为 1 K = numel(p); h_tap = (randn(K,1) + 1j*randn(K,1)) .* sqrt(p/2); % 复高斯，平均功率 p d_samp = round(delay_s(:) * fs); % 时延落到整数采样点 h_imp = zeros(max(d_samp)+1, 1); h_imp(d_samp+1) = h_tap; % 构信道冲激响应 rx_clean = conv(tx, h_imp); rx_clean = rx_clean(1:numel(tx)); % 截断到原长度 rx = awgn(rx_clean, snr_dB, 'measured'); end

round(delay_s * fs)这一步把连续时延量化到采样点，本质是把信道近似成采样间隔的整数倍。当时延不是采样周期整数倍时，会引入额外的能量弥散，短 CP 场景下能明显看出误码率抬高。严谨做法是用分数延时滤波器重采样，课程设计阶段用整数近似足够，但你得知道自己近似掉了什么。

截断的那一行也有讲究。conv之后序列变长，前面 N_cp 个样点里混着“上一个符号的拖尾”，接收端去 CP 正好把这段丢掉——这是 OFDM 用 CP 换来的好处。如果截断位置和发射端 CP 对不齐，比如多截了 3 个点，那 3 个点的偏移就会变成频域线性相位，还能被 LS 估计吸收；多截超过 N_cp，就直接掉进 ISI，曲线在 18 dB 后彻底贴平。

3.2 LS 与 MMSE：什么时候必须换成后者

上一节的接收机用的是迫零（ZF）思路：直接Y/H。它在深衰落子载波上会把噪声一起放大——某个子载波上 |H| 恰好很小，除以它之后噪声被抬上百倍，解调自然出错。这就是为什么很多人的曲线在 20 dB 之后不降反升。

MMSE 均衡多了一项噪声方差，代价是要知道 SNR：

sigma2 = 10^(-SNR_dB(s)/10); % 噪声方差（信号功率归一化为 1） H_d = H_est(data_idx, :); Y_d = Y(data_idx, :); X_mmse = conj(H_d) .* Y_d ./ (abs(H_d).^2 + sigma2); % MMSE

sigma2这里假设每个子载波上的信号功率是 1，这依赖你发射端的UnitAveragePower。如果不小心把导频功率加成 3 dB（pilot_sym = 2*(1+1j)/sqrt(2)），信道估计的幅值就会整体偏大，ZF 均衡后的星座会缩小一圈，16QAM 判决边界全乱。导频功率提升本身是个合理手段（提高估计精度），但提升之后要在接收端除干净。

两种估计的取舍很实际：LS 计算量小、不需要噪声信息，SNR 低于 15 dB 时和 MMSE 差 1～2 dB；MMSE 在低 SNR 和高频选场景下明显更稳，代价是要估噪声方差。我的习惯是先跑 LS 把链路跑通，再切 MMSE 对比，两条曲线如果差距远大于 2 dB，那基本可以确定信道估计的插值或对齐出了问题，而不是均衡器的问题。

3.3 导频图案：间隔定得对不对，看时延扩展

导频不是越多越好。频域间隔 Δk 必须满足 Δk·Δf < 1/τ_max，工程上再留 2 倍余量；时域间隔则由多普勒决定，慢衰落场景每个符号都放导频完全够用，高速场景才需要考虑时频二维导频。

按前面的参数，间隔 7 个子载波即 105 kHz，对应 τ_max ≈ 9.5 μs。你的信道模型如果 PDP 里最后一径时延是 5 μs，那就很安全；如果是 20 μs，那插值误差会主导高 SNR 表现，加多少 dB 都救不回来。这时正确的动作是把导频间隔改成 3 或 4，而不是去调均衡器参数。我见过不少仿真卡在这一步，反复改 MMSE 公式，其实问题在导频密度。

另外，导频序列本身要选低峰均比的，比如 Zadoff-Chu 或者简单的 BPSK 序列。全 1 导频做 LS 没问题，但如果你想顺手做定时或频偏估计，全 1 序列的相关峰会非常宽，定时分辨率会很差。

4. 同步环节：时偏和频偏如何让星座图散开

4.1 时偏和频偏在时域分别是什么样子

时偏分两种。整数倍时偏（偏移量小于 CP 长度）不会产生 ISI，只是在频域上引入一个线性相位斜坡 e^(-j2πkδ/N)，这个斜坡会被 LS 信道估计顺带吸进去，所以 CP 内的定时误差几乎无害——这是 OFDM 一个很反直觉的性质。超出 CP 才会引入前一个符号的拖尾，此时任何均衡都救不了。

频偏就麻烦得多。它来自收发本振的频率差，在时域表现为每个采样点乘一个 e^(j2πΔf·n/fs) 的旋转，在频域表现为所有子载波整体平移。小数倍频偏破坏子载波正交性，产生 ICI，星座图会从“点”散成“圈”，SNR 越高，圈越明显；整数倍频偏则让所有子载波错位，解出来的全是别的子载波上的数据，误码率直接 0.5。这两者的区别是排错的关键线索：圈子说明是小数的，全错说明是整数的或者根本没同步。

4.2 用 CP 自相关做定时和小数倍频偏估计

OFDM 的 CP 是数据尾部的复制，这段天然的冗余正好拿来做盲同步。取长度为 N_cp 的窗，比较它和相隔 N_fft 的那段的相关性，峰值位置就是符号起点，峰值的相位就是频偏。

function [d_hat, f_hat] = cp_sync(rx, N_fft, N_cp) % 返回：符号起点 d_hat，归一化频偏 f_hat（单位：子载波间隔） L = N_cp; P = zeros(N_fft, 1); R = zeros(N_fft, 1); for d = 1:N_fft seg1 = rx(d : d+L-1); % CP 段 seg2 = rx(d+N_fft : d+N_fft+L-1); % 被复制的那段 P(d) = seg1' * seg2; % 共轭相关，相位含频偏 R(d) = seg2' * seg2; % 能量归一化用 end metric = abs(P).^2 ./ (R.^2 + eps); % 归一化定时度量 [~, d_hat] = max(metric); f_hat = -angle(P(d_hat)) / (2*pi); % 小数倍频偏，范围 ±0.5 end

P(d) = seg1' * seg2用的是seg1的共轭乘seg2，所以相位是 +2πε，取负号才得到 ε。这个符号非常容易搞反，搞反的后果是补偿之后频偏变成两倍，星座图散得更厉害。判断方法很直接：故意给信号加一个已知的 0.2 子载波间隔频偏，看估计出来是不是 +0.2，对不上就翻符号。

metric的归一化也别省。不除以R^2的话，度量值会随信号包络起伏，峰值可能落在能量高但位置错的样点上，尤其在多径下特别明显。

4.3 补偿顺序：先小数、再整数、最后按符号重算相位

完整流程是：小数倍频偏估计并补偿 → FFT → 用导频做整数倍频偏搜索 → 频域搬移。顺序错了会互相干扰。

% 1) 小数倍频偏补偿（按全局样点索引，n 从 0 开始） n = (0:numel(rx)-1).'; rx_cfo = rx .* exp(-1j*2*pi*f_hat*n/N_fft); % 2) 整数倍频偏：频域用导频做循环移位搜索 cand = -6:6; metric_i = zeros(size(cand)); for ii = 1:numel(cand) Ys = circshift(Y, cand(ii), 1); metric_i(ii) = abs(sum(sum(conj(repmat(pilot_sym,1,N_sym)) .* Ys(pilot_idx,:)))); end [~, i_best] = max(metric_i); Y = circshift(Y, cand(i_best), 1);

补偿用的n必须是整个帧内的全局索引。如果按每个 OFDM 符号重置成 0 到 N_fft+N_cp-1，那每个符号都会引入一个额外的不连续相位跳变，等效于又加了一次频偏，现象是星座图按符号周期规律地转圈。这个坑我在第一次做同步时踩过，调了一下午公式，最后发现是循环里n的位置写错了。

多径场景下，CP 相关峰会被削弱，因为各径的 CP 到达时间不同，相关窗对不齐。这时可以把多个符号的相关值累加再做峰值搜索，代价是要多缓存一段数据。低速移动场景下这个代价很小，效果却提升明显。

5. 避坑与排查：五条会让你误判“仿真没错”的记录

5.1 星座图与频谱类异常

现象一：星座图整体旋转 45° 或缓慢画圈。原因通常有两个：一是发射端做了ifftshift而接收端忘了fftshift（或者反过来），直流位置错位导致所有子载波搬到错误频点；二是小数倍频偏没补偿，相位随符号线性增长，看上去就是星座在转。解决方法是先用理想信道（H 全 1）跑一遍，如果还转，就是索引或频偏问题；如果理想信道下不转，问题在信道估计的对齐。

现象二：频谱左右对称、出现镜像。原因基本是qammod出来的复数序列在reshape时用了转置，或者 IFFT 之前没有按列放数据，把不同符号的数据混到了同一列。解决办法是把 N_sym 设成 1，打印X(:,1)的非零位置，确认它和data_idx完全一致。

5.2 曲线形态类异常

现象三：低 SNR 正常，高 SNR 出现误码平台（平底）。按可能性排序：CP 长度短于最大时延扩展、信道 PDP 没归一化（总功率大于 1，实际 SNR 比你标称的低）、相位噪声或 IQ 不平衡模型引入的残余。先检查归一化：sum(p)必须等于 1；再把 CP 从 16 加到 32，如果平台消失，就是 ISI 引起的。

现象四：同一 SNR 跑两次，BER 差一个数量级。这在低误码率点特别常见，根因是蒙特卡洛样本不够。BER 要到 1e-4，至少需要 1e6 个比特才有统计意义（期望错误数 100 量级）。N_sym = 2000、48 个数据子载波、16QAM，一共 38.4 万比特，只能勉强看到 1e-3。另外别忘了rng固定种子，否则复现都做不了。

现象五：解调出来误码率恒等于 0.5。三个高频原因：UnitAveragePower在qammod和qamdemod两端设置不一致；比特 reshape 顺序与解调展开顺序不匹配；格雷映射被手动改过而解调还是按默认表查。排查方法很土但有效——把 SNR 设成 30 dB、信道设成理想，如果还是 0.5，那一定不是噪声问题，是映射或顺序问题。

5.3 同步与均衡类异常

现象六：加了频偏补偿，BER 反而更差。除了前面说的符号取反，还有一个隐蔽原因：补偿时用了exp(-1j*2*pi*f_hat*n/N_fft)，但f_hat是从 CP 相关估计出来的，范围只有 ±0.5；如果真实频偏是 1.3 个子载波间隔，估计值会是 0.3，补偿后还剩 1.0 的整数倍频偏，这时需要走整数倍搜索那一步。先估小数、再估整数，顺序不能颠倒。

现象七：ZF 均衡后个别子载波噪声爆炸。这是 ZF 的固有缺陷，不是 bug。判断方法：打印abs(H_est(data_idx,1))，如果最小值比均值低 20 dB 以上，那几个子载波就是噪声放大器。换成 MMSE，或者加一个判决反馈，立刻能看出差别。

6. 把 BER 曲线跑成能写进报告的证据

曲线跑出来不算完，得证明它不是巧合。我的做法是搭三条基准线，逐层加损伤，每加一层记录损失多少 dB，这样任何异常都能定位到具体环节：

A 线是底线，对不上就说明链路本身有问题，后面的都不用看。对标的写法是：

% 理论曲线：注意扣掉 CP 开销带来的信噪比损失 cp_loss = 10*log10((N_fft + N_cp)/N_fft); % 约 0.97 dB ber_ref = berawgn(SNR_dB - cp_loss, 'qam', M); semilogy(SNR_dB, ber_sim, 'o-', SNR_dB, ber_ref, 'k--'); legend('仿真','理论(AWGN)'); grid on;

如果 A 线和理论差了 3 dB 以上，先查sqrt(N_fft)那个归一化系数，再查awgn的'measured'有没有把 CP 能量算进去。这两个地方各能贡献约 1 dB，叠加起来刚好是“看起来像算法错了”的量级。

一个具体技巧：把N_sym分成 10 批跑，每批单独统计误码数，画出误码数的均值和方差。如果 10 批之间的波动超过 30%，说明样本量不够，曲线上的每个点都不可信。我现在的习惯是每个 SNR 点至少跑 100 个错误比特，跑不到就自动加大符号数——这个判断放在脚本里，比事后看图猜靠谱得多。

最后说个血泪经验：任何一次改动参数或公式后，我都会强制重跑 A 线。有一次改了导频功率忘了改回来，A 线偏了 3 dB，我在信道估计那边找了整整两天。从那以后，A 线不通过就不许看别的曲线。希望帮到你。

本文还有配套的精品资源，点击获取