简介面向本硕博无线通信教研的Matlab代码包围绕宽带瑞丽衰减信道下的误码率对比展开涵盖OFDM、OTFS、C-OFDM、C-OTFS四种方案并统一采用16QAM调制。资源共10个文件含9个Matlab脚本和1个AVI操作录像压缩包约327KB轻量但模块完整。代码按数据生成、调制解调、多径信道、均衡、SFFT/ISFFT变换及绘图等环节组织运行主脚本即可快速复现误码率曲线同时可借助操作录像了解工程路径设置、避免直接运行子函数等关键注意事项。每种方案都有对应对比曲线适合课程设计、课题预研或原理学习目前已有1256人学习浏览是上手OTFS及其改进型误码率仿真的实用参考资料。1. 宽带瑞丽瑞利衰减信道下的误码率仿真为什么 OFDM 和 OTFS 的对比值得亲手跑一遍做无线物理层仿真的人迟早会撞上同一个问题信道一变成频率选择性衰落OFDM 的误码率曲线就开始“掉链子”而教科书里那张漂亮的 16QAM 理论曲线在实测里几乎不存在。宽带瑞利信道下OFDM 靠循环前缀把多径扛过去但子载波间的正交性一旦被时变破坏ICI子载波间干扰会把高信噪比区域的天花板死死压住。OTFS 则把数据从时频域挪到时延-多普勒域等于换了个坐标系跟信道“讲道理”。至于 C-OFDM 和 C-OTFS核心区别在于循环前缀或循环后缀的使用策略两者在同样信道下的误码率表现差异恰好能说明“前缀不是越长约好”这个容易被忽略的结论。这篇文章不打算泛泛讲理论。我会直接给出可复现的 MATLAB 仿真框架用 16QAM 调制在宽带瑞利信道下对比 OFDM、OTFS、C-OFDM、C-OTFS 四条误码率曲线并把每个模块的参数设置、边界条件和踩过的坑写清楚。代码操作视频的录制思路和断点检查点也会放在最后一章单独说。2. 先把四个体制的差异讲透OFDM、OTFS、C-OFDM、C-OTFS 在瑞利信道下各自的“命门”四个体制看着都是“多载波 16QAM”的排列组合但它们对宽带瑞利信道的处理逻辑完全不同。这一章先把差异讲清楚后面给代码时你才知道每个参数在跟什么对抗。2.1 OFDM 与 C-OFDM循环前缀的“补丁”逻辑以及它为什么在高速移动下失效OFDM 的基本思路是把宽带信道切成多个窄带子载波每个子载波上经历近似平坦衰落再用循环前缀CP把多径时延“吸收”掉。CP 的长度必须大于最大时延扩展否则会产生符号间干扰ISI和载波间干扰ICI。在慢变信道里这套逻辑非常可靠16QAM 在低信噪比段的误码率曲线会明显优于单载波体制。C-OFDMCyclic-prefix OFDM是 OFDM 的改进型它的做法是在每个 OFDM 符号末尾也加一段循环后缀CS或者在接收端做额外的循环重构处理。常见做法是发送端用 CP接收端把符号末尾的采样点复制到 FFT 窗口之前人为“制造”一个更长的循环前缀。好处是能用更短的 CP 达到同样的抗多径效果频谱效率更高坏处是如果信道时变太快这种拼接处理会引入额外的相位跳变误码率反而可能比标准 OFDM 还差。仿真里最容易翻车的点是CP 长度和 FFT 点数必须严格匹配信道延迟功率谱的最大抽头位置。很多人只改 CP 长度不改信道参数结果误码率曲线怎么跑都不对就是因为最大时延扩展已经超过 CP 了OFDM 体制下不管怎么调信噪比地板效应都会出现。2.2 OTFS 与 C-OTFS时延-多普勒域的处理逻辑以及“信道的命门”为什么不在频域而在时域OTFS 的出发点是信道在时延-多普勒域里是稀疏的。它把 16QAM 符号放进一个二维网格——一个轴是时延一个轴是多普勒——每个符号占据一个固定的时延-多普勒资源格。发送端做逆辛傅里叶变换ISFFT把数据从多普勒域映射回时频域再通过 OFDM 调制发射接收端做匹配滤波和 SFFT 把信号还原到多普勒域再做均衡。这样做的本质是让每个符号经历的信道响应“几乎恒定”而不是像 OFDM 那样让每个子载波随信道变化而起伏。C-OTFS 是在 OTFS 框架里加入循环前缀/循环后缀的变体。它和 OTFS 的区别不在于多普勒域处理方式而在于保护间隔的插入策略。有些文献里 C-OTFS 也被称为“带循环扩展的 OTFS”前缀长度通常取最大时延扩展的 1.5 到 2 倍以补偿接收端做二维均衡时引入的边界效应。OTFS 的误码率优势主要出现在高多普勒场景。比如 300 km/h 的高铁场景下OFDM 的子载波正交性被多普勒扩展破坏地板效应明显OTFS 因为把时变信道“固化”到多普勒轴上能量泄漏被显著压制。但在低速场景、比如步行或者静止信道里OTFS 因为额外的变换开销和均衡复杂度误码率并不一定优于 OFDM。2.3 16QAM 在四个体制里的映射方式以及为什么调制方式不变、误码率却不同16QAM 在四个体制里都是同样的星座映射格雷码映射、每个符号 4 bit能量归一化也都一样。真正的差异在于比特到资源格的映射顺序以及均衡器拿到的信道估计精度。OFDM 里16QAM 符号直接映射到子载波均衡方式通常是单抽头迫零或 MMSEOTFS 里符号先映射到时延-多普勒网格发射端再做 ISFFT接收端做 SFFT 后还要做二维均衡。C-OFDM 和 C-OTFS 只是在保护间隔上做文章不改变调制映射本身。所以误码率差异的本质是同样的 16QAM 符号在不同资源网格上经历的信道失真程度不同均衡器恢复符号的能力也不同。高信噪比下误差主要来自均衡器的残留干扰——OFDM 的残留干扰来自 ICIOTFS 的残留干扰来自多普勒分辨率的截断效应。这解释了为什么 OTFS 曲线的斜率在 25 dB 后会变缓而不是继续笔直下降。3. 仿真框架搭建从信道模型到四体制收发机一份能直接跑的 MATLAB 代码这一章直接给可复现的代码。我以 MATLAB 为主体因为通信仿真里 MATLAB 的通信工具箱对瑞利信道和 16QAM 的支持最完整。代码按模块拆开你可以一段一段贴进脚本里跑。3.1 宽带瑞利信道模型抽头延迟线怎么设多普勒谱怎么选宽带瑞利信道用抽头延迟线模型TDL描述。每个抽头对应一条多径幅度服从瑞利分布相位服从均匀分布各抽头之间相互独立。我常用的参数是抽头数 12最大时延扩展约 1.76 μs对应 LTE 的 EPA/EVA 扩展信道模型等级。下面这段代码生成信道对象% 宽带瑞利信道12 抽头 TDL 模型 % 采样率 30.72 MHz对应 LTE 20 MHz 带宽 fs 30.72e6; % 采样率 taps 12; % 多径抽头数 maxDelay 1.76e-6; % 最大时延扩展 1.76 us avgPathGains linspace(0, -12, taps); % 功率随抽头线性衰减末抽头-12 dB pathDelays linspace(0, maxDelay, taps); % 抽头时延均匀分布 channel comm.RayleighChannel( ... SampleRate, fs, ... PathDelays, pathDelays, ... AveragePathGains, avgPathGains, ... MaximumDopplerShift, 200, ... % 最大多普勒频移 200 Hz RandomStream, mt19937ar with seed, ... Seed, 42);comm.RayleighChannel是通信工具箱里的标准对象。PathDelays决定多径到来的时刻AveragePathGains决定每条径的平均功率MaximumDopplerShift决定信道随时间变化的快慢。200 Hz 多普勒对应约 3 km/h 的移动速度属于低速场景如果你要跑高速场景把 200 改成 1000 或 2000误码率曲线的差异会非常明显。RandomStream设成带种子的 mt19937ar是为了保证每次仿真结果可重复这一点在对比四个体制时必须做否则曲线差异混入了随机噪声。3.2 OFDM 与 C-OFDM 的收发机代码循环前缀和循环后缀的插入方式差异OFDM 的发射机核心是 IFFT、加 CP、并串转换。C-OFDM 则在 IFFT 后同时拼接循环前缀和循环后缀。下面是两者的对比代码% OFDM 调制 % numSubcarriers 为子载波总数cpLen 为循环前缀长度 function txSignal ofdm_modulate(symbols, numSubcarriers, cpLen) % symbols: 长度为 numSubcarriers 的 16QAM 符号列向量 ifftOut ifft(symbols, numSubcarriers); txSignal [ifftOut(end-cpLen1:end); ifftOut]; % 添加循环前缀 end % C-OFDM 调制多了一个循环后缀 function txSignal cofdm_modulate(symbols, numSubcarriers, cpLen, csLen) ifftOut ifft(symbols, numSubcarriers); cp ifftOut(end-cpLen1:end); cs ifftOut(1:csLen); txSignal [cp; ifftOut; cs]; % 前缀 数据 后缀 endcofdm_modulate里后缀的作用是在接收端给 FFT 窗口提供更多可选的对齐位置。实际接收时C-OFDM 接收端通常会做“滑动窗口搜索”找到误码率最低的 FFT 起始点。这个搜索的开销不算大但代码量会明显增加所以很多简化仿真里干脆省略后缀拼接只保留前缀。省略后 C-OFDM 就退化成 OFDM对比意义就没了这一点务必注意。接收端做 FFT 前要去掉 CP。OFDM 直接去掉前 cpLen 个采样C-OFDM 要去掉前缀加后缀的总长再根据信道估计结果调整 FFT 窗口起点。代码如下% OFDM 解调 function symbols ofdm_demodulate(rxSignal, numSubcarriers, cpLen) rxSignal rxSignal(cpLen1:cpLennumSubcarriers); symbols fft(rxSignal, numSubcarriers); end % C-OFDM 解调移除 CP 和 CS并从信道估计中估计最佳窗口起点 function symbols cofdm_demodulate(rxSignal, numSubcarriers, cpLen, csLen, channelEst) totalGuard cpLen csLen; % 粗略去掉前后保护间隔 rxCore rxSignal(totalGuard1:totalGuardnumSubcarriers); % 依赖信道估计做窗口微调 [~, bestShift] max(abs(channelEst)); rxCore circshift(rxCore, -bestShift); symbols fft(rxCore, numSubcarriers); endC-OFDM 接收端的circshift是根据信道冲激响应最大抽头位置做符号定时同步。这一行是 C-OFDM 能比标准 OFDM 抗多径的关键也是最容易被误写的地方——很多人直接把 OFDM 的接收代码套过来不做窗口微调结果 C-OFDM 的性能和 OFDM 完全一样对比白做。3.3 OTFS 与 C-OTFS 的收发机代码ISFFT / SFFT 变换矩阵与网格映射OTFS 发射机先把 16QAM 符号排列成 M×N 的时延-多普勒网格再做 ISFFT 变换到时频域最后用 OFDM 调制发射。接收端做 SFFT 把信号还原回多普勒域。关键是维度匹配时延轴长度等于子载波数多普勒轴长度等于 OFDM 符号数。% OTFS 发射机 % xDD: M×N 的时延-多普勒域符号矩阵M子载波数N符号数 function txSignal otfs_modulate(xDD) % Step 1: ISFFT把 DD 域映射到 TF 域 xTF fft(ifft(xDD, [], 1), [], 2); % 沿时延轴做 IFFT沿多普勒轴做 FFT % Step 2: 把 TF 域矩阵拉直按列发送每列是一个 OFDM 符号 [M, N] size(xDD); txSignal zeros(M*N, 1); for n 1:N txSymbol ifft(xTF(:, n), M); % 子载波维度的 IFFT txSignal((n-1)*M1:n*M) txSymbol; end end % OTFS 接收机 function xDD_est otfs_demodulate(rxSignal, M, N, channelEst) % Step 1: 按列恢复 TF 域符号 xTF zeros(M, N); for n 1:N rxSymbol rxSignal((n-1)*M1:n*M); xTF(:, n) fft(rxSymbol, M); end % Step 2: SFFT把 TF 域映射回 DD 域 xDD_est fft(ifft(xTF, [], 1), [], 2); % Step 3: 在 DD 域做简单 MMSE 均衡这里用逐点除法近似 xDD_est xDD_est ./ channelEst; end这段代码里fft(ifft(xDD, [], 1), [], 2)的顺序最容易错。ISFFT 的数学定义是先沿时延轴做 IFFT、再沿多普勒轴做 FFT两个轴的顺序不能颠倒。我见过很多版本把第二个fft写成了ifft出来的星座图永远是一团散点根本没法解调。另外xDD ./ channelEst是简化处理——真实系统应该用二维 MMSE 滤波器但那种复杂度对初版对比来说没有必要。C-OTFS 的发射机在 OTFS 的基础上往每个 OFDM 符号前加 CP。它和标准 OTFS 的区别是OTFS 里每个符号的 CP 长度可能不一致取决于多普勒轴索引而 C-OTFS 统一用固定长度 CP。这种“固定化”实现简单代价是多普勒轴上的边界符号会有轻微的能量泄漏。如果仿真的多普勒扩展不大100 Hz 以内C-OTFS 和 OTFS 的误码率曲线几乎重合多普勒拉高后OTFS 的优势才会显现。3.4 四个体制要统一参数和对比基准公平性比“跑得好看”更重要对比误码率最容易出现的伦理问题是只有四条曲线参数却各自偏向对自己有利的配置。比如给 OTFS 用更大的网格、给 OFDM 用更短的 CP出来的对比结论没有意义。我固定以下参数FFT 点数 512有效子载波 128CP 长度 3216QAM 调制蒙特卡洛次数 2000。四个体制都用同一组信道种子、同一帧数据长度。这样跑出来的从 0 到 25 dB 信噪比点才算是“公平决斗”。固定参数之后还有一个必须说的点OTFS 和 C-OTFS 因为引入了多维变换它们的误码率理论上对信噪比的敏感度略低于 OFDM——这不是坏事而是它们的抗时变特性带来的结果。低速信道里 OFDM 在 15 dB 之后会快速下降OTFS 则从 10 dB 开始就拥有更低的绝对误码率。4. 完整仿真流程与参数表从生成数据到画出四条误码率曲线的具体步骤4.1 仿真链路总览每个模块的输入输出、数据降采样和 SNR 换算整体链路包括随机比特生成 → 16QAM 映射 → 四种调制体制 → 信道卷积 → AWGN 加噪 → 信道均衡 → 解调 → 星座解映射 → 误码率统计。这里有一个关键点加噪的信噪比要按“信息比特”计算不是按“调制符号”能量计算。16QAM 每个符号带 4 bit符号能量和比特能量的关系是 Es 4Eb。% 信噪比换算EbNo - SNR EbNo_dB 0:5:25; % 比特信噪比dB EsNo_dB EbNo_dB 10*log10(4); % 符号信噪比 snr EsNo_dB - 10*log10(fftLen/numActiveSubcarriers); % 考虑频谱效率后的实际SNR10*log10(fftLen/numActiveSubcarriers)这一项容易被忽略。FFT 点数 512、有效子载波 128 时这一项等于 6 dB。如果不扣除你会发现自己加噪后的误码率曲线比理论值“好”了约 6 dB而且怎么查都查不出原因。这种现象我见过太多次先在 4.4 的避坑里详细说。4.2 主循环代码蒙特卡洛仿真怎么写帧结构怎么编排主循环结构如下每个信噪比点重复 2000 次每次生成随机比特、调制、过信道、加噪、解调、统计错误比特数。帧结构为每个 OFDM 符号 128 个有效子载波每个子载波映射一个 16QAM 符号每个帧包含 14 个符号模拟一个时隙。OTFS 的网格采用 128×14正好和 OFDM 的总数据量对齐——这也是“公平对比”的另一个维度。numFrames 2000; % 蒙特卡洛次数 ber zeros(4, length(EbNo_dB)); for i 1:length(EbNo_dB) errorBits zeros(4, 1); totalBits zeros(4, 1); for frame 1:numFrames % 生成随机比特16QAM 映射 dataBits randi([0 1], 128*14*4, 1); symbols qammod(dataBits, 16, gray, InputType, bit); % 四种体制发射 txOFDM ofdm_modulate(symbols, 128, 32); txCOFDM cofdm_modulate(symbols, 128, 32, 16); xDD reshape(symbols, [128 14]); txOTFS otfs_modulate(xDD); % OTFS 发射信号补 CP作为 C-OTFS txCOTFS [txOTFS(1:32); txOTFS]; % 简化的 C-OTFS统一加 32 长度 CP % 过信道四种发射信号各自过信道但种子相同 rxOFDM channel(txOFDM); rxCOFDM channel(txCOFDM); rxOTFS channel(txOTFS); rxCOTFS channel(txCOTFS); % 加 AWGN rxOFDM awgn(rxOFDM, snr(i), measured); rxCOFDM awgn(rxCOFDM, snr(i), measured); rxOTFS awgn(rxOTFS, snr(i), measured); rxCOTFS awgn(rxCOTFS, snr(i), measured); % 解调、均衡、误码率统计省略各体制解调内部实现 % ... end ber(:, i) errorBits ./ totalBits; end上面代码里 C-OTFS 的简化是用[txOTFS(1:32); txOTFS]加前缀。严格说来 C-OTFS 应该在每个“OFDM 符号块”前都加 CP而不是在整个帧前加一个 32 长度的 CP。对整个帧加和逐符号加在过信道后的首尾边界表现不同。如果要把 C-OTFS 做成严格版就得把txOTFS按每个符号切分开分别加 32 长度的 CP再拼接。4.3 导频与信道估计最小二乘估计在四个体制里的不同用法误码率仿真的第二个陷阱是信道估计的处理方式不同。很多入门脚本直接假设接收端已知完美信道perfect CSI用信道矩阵直接均衡。这种假设对 OFDM 尚可接受对 OTFS 就不公平了——OTFS 的优势部分来自它对时变信道的估计能力一旦用完美 CSI这个优势就被抹掉了。常规做法是插入导频OFDM 在指定子载波插入已知符号接收端在对应的子载波位置做最小二乘估计再插值得到全信道矩阵。OTFS 通常在时延-多普勒域插入一个“导频符号孤岛”——在一个资源格中放入已知符号周围置零接收端从这个孤岛中估计出稀疏信道响应。下面是最小二乘估计的通用片段% 最小二乘信道估计已知导频处发送符号 txPilot 与接收符号 rxPilot H_ls rxPilot ./ txPilot; % 逐点相除得到信道估计 % 插值OFDM 场景或稀疏重构OTFS 场景 H_interp interp1(pilotIdx, H_ls, allIdx, linear, extrap);如果你用的是带种子的comm.RayleighChannel那么信道估计的误差是可重复的这便于调试。我调试时会先关闭 AWGN、只用信道跑一遍检查各体制均衡后的星座图是否有明显的聚类再逐步加噪声。如果加噪后星座图在低信噪比下出现严重旋转变换优先怀疑信道估计的相位没有正确补偿。4.4 四个必查的仿真参数和容易出现“假性能”的位置建议你把下面这张表打出来放在屏幕旁边每一轮调试都逐个核对。它几乎是“四体制对比”仿真里所有翻车事故的高发区。参数项OFDMC-OFDMOTFSC-OTFS易错点FFT 点数512512512512OTFS 里 FFT 点数等于时延轴长度有效子载波128128128128导频占用后有效数据子载波要少一点CP 长度3232按需32CP 比最大时延短时 OFDM 直接地板多普勒频移200 Hz200 Hz200 Hz200 Hz四个体制必须统一否则对比无效符号映射gray 16QAMgray 16QAMgray 16QAMgray 16QAMOTFS 网格的排列方向影响误码率分布均衡方式单抽头 MMSE滑动窗口 MMSE二维 MMSE二维 MMSE均衡维度写错时曲线会异常平缓信噪比折算扣 6 dB 频谱效率再扣循环后缀开销扣网格内保护开销扣除 CP 开销不扣开销时曲线整体左移每帧符号数14141414OTFS 网格 N14与 OFDM 符号数对齐5. 避坑与常见问题排查跑出的结果不理想先检查这五个位置5.1 现象误码率曲线在高信噪比处突然变成水平线怎么都降不下去原因几乎一定是 ICI 或残余 ISI 主导了误码率。OFDM 体制下CP 长度不足或信道多普勒扩展过大都会导致这种地板效应。OTFS 体制下如果网格维度过小比如时延轴只有 64 格而多径有 12 抽头多普勒分辨率不足也会形成地板。解决先检查信道对象的最大多普勒频移再检查 CP 长度是否大于pathDelays的最后一个抽头时延。我的经验值是 CP 至少取最大时延的 1.25 倍留出裕量。OTFS 的话把时延轴宽度调成子载波数的 2 倍试试。5.2 现象OTFS 的误码率居然比 OFDM 差而且差很多这种情况通常不是算法问题而是网格映射方向和 ISFFT 的定义不一致。fft(ifft(xDD, [], 1), [], 2)与ifft(fft(xDD, [], 1), [], 2)只差一个变换方向但对应完全不同时延-多普勒映射规则。如果你发射端写的是 ISFFT、接收端写的也是 ISFFT那就等于没有完成域变换误码率自然崩坏。解决打印发射端的xTF矩阵检查它是否具有明显的行、列相关性特征——ISFFT 之后的矩阵应该是“能量铺开”的不应该还有明显的星座聚类。5.3 现象四条曲线在 0-5 dB 低信噪比区域几乎重合这其实是正常现象但如果你期望看到明显差异那就要检查是不是 16QAM 在低信噪比下的判决错误已经饱和。低信噪比下噪声主导一切信道体制的差异被噪声淹没。区别拉开的分水岭在 10 dB 以上如果 10 dB 之后四条曲线还是完全重合那就要检查是不是信道对象没有在每次迭代中“推进”——比如你把带种子的信道对象放在循环外面且没有调用reset或advance导致所有帧经历了完全相同的信道实现。解决确认信道对象在每次循环里都被调用或者在不同信噪比点之间用reset(channel)重置。5.4 现象误码率低得不正常好像比理论 AWGN 曲线还低这就是 4.1 里提到的频谱效率扣减问题。你给 AWGN 加的功率是按照符号能量算的但理论曲线是按照比特能量算的两者之间差了log2(M)。16QAM 时差了 4 倍即 6 dB。如果snr里没有扣掉这个数你的曲线就会左移 6 dB。另一个低频翻车点是awgn函数里的measured选项——它会根据输入信号功率自动调整噪声功率如果你输入的是过信道后的信号它会把信道增益也考虑进去导致实际信噪比比设定值高 1-2 dB。我建议自己计算噪声功率不依赖measured。5.5 现象C-OFDM 和 OFDM 曲线完全重合看不出后缀的作用先确认你的 C-OFDM 接收端做了窗口微调也就是 3.2 中的circshift。如果你的 C-OFDM 实现没有这个调整那它本质上还是 OFDM。另外信道时延如果很小比如只有 2-3 个抽头C-OFDM 的后缀带来的增益会被 AWGN 的随机性淹没。解决把最大时延扩展到 2.5 μs 以上并且把多普勒频移调到 500 Hz 以上再看两条曲线是否分离。6. 把仿真过程录成可检查的视频断点设计、画面布局和录屏检查表6.1 录屏分辨率、帧率与关键画面布局代码操作视频不需要华丽的剪辑但画面布局要让人一眼看清“现在跑到哪一步了”。我习惯用 1920×1080 分辨率、30 帧率录制左侧放 MATLAB 编辑器右侧放命令行窗口底部固定显示当前的 SNR 点和误码率累计值。这样的布局能让看视频的人随时判断仿真进度而不必盯着光标猜。录制工具用 OBS Studio 就够了不需要额外付费软件。6.2 每个关键节点设置暂停点而不是一口气跑完视频不是让你表演盲打而是让人能跟着复现。我在视频里会在以下位置停顿一次第一是信道参数设置完成后把延迟功率谱打印出来第二是 OFDM 和 C-OFDM 的收发机代码分别写完后再运行一次显示星座图第三是 OTFS 的 ISFFT 变换完成后显示变换前后的能量分布第四是主循环跑完后显示四条误码率曲线叠加图。停顿期间用鼠标在关键行画圈并用一两句话解释“这里为什么要这样写”。6.3 验证视频结果是否可信的三个自查动作视频录完后我会做三个自查动作。第一个是把视频里的最终误码率曲线截图和代码重新跑一遍得到的结果对比检查有没有因为录制时的环境差异导致结果不同。第二个是检查视频里的代码块是否完整尤其是ofdm_modulate和otfs_modulate这类放在独立段落里的函数——很多录制工具会由于缩进过多把结尾的end吞掉观看者复制后一运行就报语法错误。第三个是检查是否在视频开头和结尾各放了一次完整的clear; close all; clc;否则观众观看中间代码时工作区里的残留变量会严重影响结果。6.4 最后一点习惯把信道种子写进代码注释而不是依赖随机数生成器的默认状态这是我吃过一次亏后养成的习惯。早期做对比仿真时我觉得“反正多条曲线跑完再对比有没有种子无所谓”结果在复现一个 OTFS 论文结果时发现无论怎么调整参数曲线都和论文差 2 dB。排查了两天最后发现是 MATLAB 随机数状态和论文作者的初始状态不一致。把Seed 42写在信道对象里、把rng(42)放在脚本开头并且连同注释一起写进视频要展示的代码里是成本最低的“后悔药”。每一个想在论文或项目里引用这个结果的人都会因为这一行而少走两天弯路。仿真做到这个程度剩下的就是你自己调参和对比的功夫了。我最后想说的是在宽带瑞利信道下对比这四种体制如果你能亲自把代码跑通并且亲手见证 OTFS 在 300 Hz 多普勒下的优势、以及 OFDM 在低多普勒时反超 OTFS 的现象你对“无线通信为什么需要不同的波形设计”的理解会比只读十篇综述都深。希望帮到你。本文还有配套的精品资源点击获取