ISAC通感一体化论文代码复现：MATLAB仿真避坑与验证指南

简介：基于MATLAB的ISAC通感一体化论文阅读与代码复现毕业设计源码包，项目源自SEU SISE毕设，答辩评审分达到98分。内容围绕毫米波MIMO通感一体化系统的混合波束成形（Hybrid Beamforming）核心实现，并配套2022年相关英文论文PDF，便于读者结合原文逐段对照仿真细节，理解ISAC信道建模与波束设计思路。资源共6个文件，含5个MATLAB脚本和1篇论文文献，压缩包仅355KB，轻量易部署，代码均已调试测试、可直接运行。已有527人学习下载，适合通信、计算机、人工智能、自动化等专业学生及从业者，作为期末课程设计、课程大作业或毕设参考资料。项目从文献阅读到算法复现形成较完整闭环，入门者可按脚本逐步消化，进阶者可在既有代码框架上修改调整参数，以适配不同场景需求。

1. ISAC通感一体化论文复现：为什么照着论文写代码也跑不出结果

ISAC（Integrated Sensing and Communication，通感一体化）是当前无线方向的热门课题，核心思想是让同一份波形在传输数据的同时完成距离、速度的感知。某高校SISE毕业设计里，“ISAC通感一体化论文阅读与代码复现基于Matlab实现”这类题目出现频率很高，表面看就是读论文加写仿真，实际动手才发现论文里的公式、参数和图表经常对不上，照着翻译成代码很容易翻车。这套题的交付物通常是源码加文档说明，而文档能不能把设计逻辑写清楚，往往比代码跑通本身更影响成绩。这篇笔记会重新拆解这套题目的完整落地路径：论文怎么读、MATLAB工程怎么搭、核心估计代码怎么写、哪些参数一改就崩，以及如何验证自己没有“假成功”。适合正在做ISAC仿真课题、手头有论文但没有标准答案的毕设选手。

2. ISAC通感一体化论文怎么读：路线选型与三遍法

ISAC通感一体化论文的读法和普通通信论文不一样：普通论文看懂系统模型就行，ISAC论文必须读懂收发两端怎么共用波形，因为你的MATLAB代码要同时体现通信和感知两条支路。这一步做不好，后面所有代码都是空中楼阁。

2.1 三种一体化路线怎么选：OFDM、FMCW与OTFS对比

常见的ISAC论文里有三种主流波形方案，先看差异，再决定要不要在这个方向上下手。

我一般会建议毕设优先选OFDM方案。理由有三个：第一，OFDM的感知解析模型最成熟，论文里给出的是闭式解而不是纯仿真曲线，方便核对结果；第二，MATLAB对OFDM的支持最好，不需要额外复杂的信道模型就能跑通全链路；第三，它的感知和通信支路相对独立，就算通信代码有bug，距离速度估计还能出结果，分模块调试很舒服。FMCW适合雷达背景更强的题目，OTFS除非论文本身就给了完整代码框架，否则不建议在毕设周期内从零复现。

选好波形方案后，再回来看论文的仿真条件。ISAC论文的仿真参数集中在“系统参数表”或“仿真环境”小节，常见参数组是载频、带宽、子载波数、符号数、循环前缀长度、调制阶数，以及目标距离和速度。这些参数直接决定后续params.m里写什么，所以读论文的第一件事不是看公式，而是先把参数表找出来。

2.2 三遍法读ISAC通感一体化论文：抓什么、抄什么、还原什么

我习惯用三遍法读这类论文，每遍目标不同，缺一遍后期都要补课。

第一遍抓系统框图。只读摘要、引言和系统模型部分，画出收发链路：发射端做什么处理、信号经过什么信道、接收端分成哪两条支路。这一步的目的是建立代码模块清单，比如“发射端需要波形生成模块”“接收端需要同步、去CP、FFT、二维估计、解调模块”。没有模块清单就直接写代码，是后面返工的最大原因。

第二遍抄参数表和公式。把论文里的参数表逐行抄成一份MATLAB脚本或表格，每个参数后面标注它在论文里的符号和单位。抄的时候留意论文有没有给距离分辨率和速度分辨率的公式，这两个公式是后面验证代码正确性的标尺。如果论文只给了参数没给分辨率，需要自己补算。

第三遍把公式翻译成伪代码。重点做发射波形和接收估计两处：发射端看频域符号怎么映射、IFFT和循环前缀在哪个环节插入；接收端看回波经过FFT后哪个维度对应距离、哪个维度对应速度，以及峰值索引怎么换算成物理量。遇到论文公式和参数表对不上的情况，不要硬猜，用参考文献里同类公式交叉验证，确定哪个正确再往下写。

三遍法不是线性走完就结束，我通常在写代码过程中会反复回看第二遍和第三遍的笔记。代码跑不通、结果对不上，先回来查参数和公式，而不是改代码碰运气。

2.3 先列输入输出清单：复现前避免白干半小时

正式建工程之前，花十分钟写一份输入输出清单。例如：

输入：载频fc、带宽bw、子载波数Nsc、OFDM符号数Nofdm、循环前缀长度cp_len、调制阶数M、目标距离R、目标速度v、信噪比SNR 输出：距离估计值range_est、速度估计值speed_est、距离多普勒谱rd_map、误码率BER

这份清单有两个作用。一是让代码模块的接口提前确定，写发射模块时就知道输出要被接收模块怎么消费，写接收模块时也知道输入从哪来。二是在最终验证阶段有一个验收标准：输入一个已知位置和速度的目标，输出估计值和真实值的误差应该在一个距离分辨率和速度分辨率以内。清单不需要写进论文，但它是你自己对“复现成功”的定义。

清单里最容易漏掉的是SNR的定义方式。论文里BER曲线横轴写的可能是Eb/N0、SNR或者Es/N0，三者的换算公式不一样，噪声功率的生成方式也不一样。不管定义直接把SNR输入进去，最后曲线对不上论文是必然的。我通常会在清单里单独标注“本项目SNR按信号功率/噪声功率计算，每个子载波等效”，并且把换算关系写在代码注释里。

3. 从论文到MATLAB工程：目录结构与OFDM波形生成

前面把论文拆清楚了，这一章开始建工程。很多ISAC毕设代码是从一段main.m开始越写越长，最后几百行的变量混在一起，调一次参要滚一遍全文。通感一体化仿真涉及发射、信道、接收、指标统计四个环节，建议一开始就按功能分模块组织工程。

3.1 工程目录这样组织：每个文件只干一件事

下面是一个我常用的MATLAB工程结构：

isac_bachelor/ ├── config/ │ └── params.m # 全部系统参数集中管理 ├── transmitter/ │ ├── data_symbols.m # 生成QAM数据符号与导频符号 │ └── ofdm_waveform.m # IFFT、插循环前缀，输出时域波形 ├── channel/ │ ├── target_channel.m # 目标回波模型：时延+多普勒 │ └── awgn_channel.m # 加噪声 ├── receiver/ │ ├── sync_demod.m # 去CP、FFT、提取频域符号 │ ├── range_doppler.m # 二维FFT，输出距离-多普勒谱 │ └── demap_symbols.m # 通信支路解调 ├── metrics/ │ ├── ber_curve.m # 循环跑SNR点，画误码率曲线 │ └── det_prob.m # 统计检测概率与虚警率 ├── main.m # 入口脚本，串联全流程 └── README.md # 记录参数来源和改动

每个文件对应一个职责，main.m只负责调用和画图。这么做最大的好处是调试时可以单独对receiver模块传入已知回波信号，判断是算法问题还是上游模块问题，不用每次从头跑全链路。另外目录结构本身也方便在论文里写“系统设计与实现”这一章，截图目录树比贴一大段代码更能说明设计思路。

需要强调的是，MATLAB脚本的命名不要带中文和空格，路径里也不要出现中文目录名，否则函数调用在部分旧版本环境下会报找不到文件。这个细节出现过不止一次——不是代码错，是环境不认名字。

3.2 系统参数集中管理：ISAC仿真参数表怎么落进params.m

参数集中放在config/params.m里，用脚本形式而不是函数形式，方便main.m直接执行后所有变量出现在工作区，调试时能直接看到每个参数的值。

% config/params.m % 通感一体化仿真参数，对应论文正文仿真参数表 % 链路参数 fc = 24e9; % 载波频率 24GHz bw = 100e6; % 信号带宽 100MHz M = 16; % 16QAM 调制 num_bits = 10000; % 参与统计的比特数 % OFDM帧结构 Nsc = 64; % 子载波数，决定FFT长度 Nofdm = 32; % 一帧OFDM符号数 cp_len = 16; % 循环前缀长度（采样点数） delta_f = bw / Nsc; % 子载波间隔 1.5625MHz fs = Nsc * delta_f; % 采样率等于带宽 T_sym = 1 / delta_f; % 单个OFDM符号有效时长 T_cp = cp_len / fs; % 循环前缀时长 T_ofdm = T_sym + T_cp; % 完整符号时长（含CP） % 感知目标与评价指标 target_range = 50; % 目标距离 50m target_vel = 10; % 目标速度 10m/s c0 = 3e8; % 光速 R_res = c0 / (2 * bw); % 距离分辨率 c/(2B)，这里是1.5m v_res = c0 / (2 * fc * Nofdm * T_ofdm); % 速度分辨率

这段脚本里，delta_f由带宽除以子载波数得到，采样率fs等于带宽是OFDM仿真的常见设定，好处是距离分辨率公式c0/(2bw)能直接用于验证。T_ofdm计算了带循环前缀的完整符号时长，这个值决定回波脉冲串的等效重复周期，速度分辨率公式里积累时间用的是NofdmT_ofdm，而不是单符号时长。

目标距离50米对应回波时延约0.33微秒，在带宽100MHz下距离分辨率1.5米，目标应该落在第33、34个距离bin附近；目标速度10m/s对应的多普勒频移约1600Hz，在帧长范围内应该落在第3、4个多普勒bin附近。这些预估值是后面核对代码的重要参照。修改参数时只需改这一处，不会出现“这个值在代码里出现过三次只改了两处”的低级错误。参数命名尽量和论文符号对应，写论文时可以直接引用那一章。

3.3 发射波形生成：IFFT与循环前缀的作用

发射模块要回答的问题是：论文里那个频域符号矩阵，怎么变成一串可以送到信道里的时域波形。核心代码如下：

% transmitter/ofdm_waveform.m % 输入: data_sym 为 [Nsc * Nofdm] 的频域QAM符号矩阵 % 输出: tx_wave 为一帧时域发射波形 function tx_wave = ofdm_waveform(data_sym, pilot_idx, params) % 插导频：导频子载波填常数，便于接收端做相位校准 data_sym(pilot_idx, :) = 1 + 1j; % 沿子载波维做IFFT，每列对应一个时域OFDM符号 tx_time = ifft(data_sym, params.Nsc, 1); % 加循环前缀：把每个符号尾部cp_len个点搬到头部 tx_cp = [tx_time(end-params.cp_len+1:end, :); tx_time]; % 串行化输出一维波形 tx_wave = tx_cp(:); end

逻辑上分三步：先按论文指定的导频位置插入已知符号，这是ISAC收发两端共同约定的先验信息；然后对每个OFDM符号沿子载波维度做IFFT，IFFT把频域复符号变成时域采样点；最后把每个符号的尾部样本复制到头部形成循环前缀。循环前缀在这里有两个作用，一个是让发射波形能够对抗多径时延扩展，另一个是让回波到达时依然保持子载波正交性——对ISAC来说后者更关键，因为测距依赖子载波相位，正交性一旦被破坏，距离估计就会产生系统偏差。

导频位置pilot_idx如果设为1:8:Nsc，实际是把所有子载波分成8组，每组第1个做导频。导频密度需要权衡：密度越高相位校准越准，但有效数据率越低。毕设仿真一般取1/8到1/16的密度即可，论文里通常会给出具体导频图样，优先按论文来。导频符号填1+1j表示功率归一化的QPSK导频，如果你用的是平均功率为1的16QAM数据符号，导频功率略高，这在仿真里是允许的，通信估计时再做幅度归一化即可。

生成波形后可以先plot(real(tx_wave(1:256)))看一眼时域波形，确认没有NaN或全是0，再做信道和接收处理。一个常见的低级错误是ifft维度写错：MATLAB的ifft默认沿第一个非单例维度运算，这里显式指定维度参数1，也就是沿子载波维度变换，避免二维矩阵时出现维度歧义。

4. 通感一体化接收与参数估计：MATLAB核心实现

波形发射出去之后，回波经过目标反射和噪声叠加，接收端拿到的是叠加了时延和多普勒频移的发射波形。这一章把接收端的两条核心支路写透：感知支路的距离-多普勒二维估计，通信支路的符号解调与误码率统计。

4.1 通感一体化接收：用二维FFT估计距离和速度

OFDM通感一体化感知的基本原理：回波相对发射波形的时延τ，会让频域每个子载波叠加一个随子载波索引线性变化的相位；目标运动产生的多普勒频移fd，会让相邻OFDM符号之间叠加一个固定相位差。因此，沿子载波维做IFFT可以得到时延信息，沿符号维做FFT可以得到多普勒信息，两者联合就是二维谱。

% receiver/range_doppler.m % 输入: rx_grid 为 [Nsc * Nofdm] 的频域接收符号矩阵（已同步、去CP、FFT） % 输出: rd_map 归一化距离多普勒谱；range_est/speed_est 为估计值 function [rd_map, range_est, speed_est] = range_doppler(rx_grid, params) % 沿子载波维做IFFT，把时延信息搬到距离维 range_ifft = ifft(rx_grid, params.Nsc, 1); % 沿符号维做FFT，把多普勒信息搬到速度维 doppler_fft = fft(range_ifft, params.Nofdm, 2); % 取幅度，转成二维谱并归一化 rd_map = abs(doppler_fft); rd_map = rd_map / max(rd_map(:)); % 峰值索引搜索，注意FFT原点在矩阵的(1,1)位置 [d_idx, v_idx] = find(rd_map == max(rd_map(:)), 1); % 索引换算距离：第一个bin对应距离0 range_est = (d_idx - 1) * params.c0 / (2 * params.bw); % 索引换算速度：先换算多普勒频率，再换算速度 fd_est = (v_idx - 1) / params.Nofdm / params.T_ofdm; speed_est = fd_est * params.c0 / (2 * params.fc); end

关键细节在索引换算。MATLAB的FFT输出从0频开始，所以第一个bin是0距离和0速度，索引必须减去1再乘以分辨率。距离分辨率是c0/(2bw)，速度分辨率是c0/(2fcNofdmT_ofdm)——注意分母里是Nofdm*T_ofdm，也就是一帧的总时长，对应FFT积累时间。有人会把分母写成单符号时长，那算出来的就不是FFT bin对应的速度步长，估计结果会和设定值完全对不上。

按第3.2节的参数跑，距离分辨率1.5米，目标50米应该落在第34个距离bin附近（索引33）。目标速度10m/s对应的多普勒约1600Hz，在多普勒FFT里大概落在第3个bin附近（索引2）。所以看到range_est在50米左右、speed_est在10m/s左右，说明感知支路逻辑基本正确；如果估计值偏了几百米，先把回波模型里有没有正确加入时延和多普勒检查一遍，再回到FFT维度检查。

4.2 通信支路：信道估计与误码率计算

感知做完，通信支路不能只停留在“解调成功”。ISAC论文的通信性能指标最常见的是误码率曲线，实现思路是：对接收频域符号做单抽头均衡，再做QAM解调，和发射比特对比，统计错误比例。

% receiver/demap_symbols.m % 输入: rx_sym 为解调后的频域符号，ch_ref 为参考信道响应 % 输出: ber 当前SNR点下的误码率 % 单抽头均衡：用导频处估计出的信道响应校正幅度和相位 rx_eq = rx_sym ./ ch_ref; % 16QAM解调，'UnitAveragePower'表示符号平均功率为1 rx_sym_idx = qamdemod(rx_eq(:), params.M, 'UnitAveragePower', true); % 十进制符号索引转二进制比特流，注意高低位顺序和发射端一致 rx_bits = de2bi(rx_sym_idx, log2(params.M), 'left-msb'); % 对比发射比特，统计误码率 ber = sum(rx_bits ~= tx_bits) / numel(tx_bits);

ch_ref怎么获得是关键。在纯仿真里，可以直接用真实信道响应当作理想信道估计，这样得到的是理想信道估计下的性能上界，论文里通常会标注CSI known或ideal channel estimation。如果想模拟更真实的情况，可以用导频子载波处的接收符号除以发射导频符号，再插值得到所有子载波的信道估计。两种做法的BER曲线会差1到3dB，对比论文时先确认对方用的是哪一种。

qamdemod需要通信工具箱支持，如果环境工具箱不全，可以在代码里用最小距离判决代替解调函数，不依赖工具箱。我见过因为目标环境缺少工具箱导致qamdemod直接报函数未定义的翻车，所以跑BER曲线前先在命令行执行ver('Communication_Toolbox')确认工具箱可用，不可用就换最小欧氏距离判决的替代实现。

4.3 性能指标怎么对比论文：指标对齐比公式对齐更重要

复现ISAC论文后阶段最容易出现的问题是：自己算的指标和论文对不上，但算法逻辑每一步都对。原因通常不在代码而在指标定义。下表是两类指标的关键对齐项。

对齐时先看坐标轴标签，再看正文对条件的描述。比如论文的检测概率曲线横轴写的是SNR，正文里却可能写“每个子载波的信噪比”，这个前提如果没在代码里体现，曲线差3dB很正常。建议在对比图上保留论文原图作为虚线参考，同一坐标下自己数据用实线，两条线差多少一目了然。这一步也是答辩时老师最爱追问的地方：“你说复现好了，差多少、差在哪？”能答清楚，比单纯把曲线画漂亮更能加分。

5. ISAC通感一体化复现中的5个避坑记录：现象、原因与解决

ISAC通感一体化仿真里的玄学不少。下面五条都是实际调试中出现过的翻车记录，按现象、原因、解决写清楚，排查时可以当成对照表用。

5.1 FFT报错或结果全黑：先查维度和数据类型

现象：range_doppler.m运行报“Matrix dimensions must agree”，或者不报错但rd_map全是一个常数。

原因：rx_grid的维度不是[Nsc, Nofdm]。回波生成时把时域波形reshape错了，或者发射端ifft用了默认维度，二维矩阵在变换时方向搞反，导致接收端每个符号对应错位。另一个原因是没有去CP，rx_grid每个符号的采样点数和发射端不一致。

解决：在range_doppler第一步加一行size(rx_grid)，用第3.2节参数确认输出应该是64×32。不对就先回去检查发射端tx_time的维度：IFFT后是64×32，加入循环前缀后每列变成80×32，串行化后再在接收端按Nsc截断并一列一列恢复成64×32。维度对上了算法才能谈正确性。

5.2 BER曲线比论文差10dB：先查SNR定义再查导频开销

现象：通信模块和论文方法一致，语法上挑不出错，但BER曲线整体比论文右移了10dB以上。

原因：最常见的两个错误。一是噪声功率用了双边谱密度而没有针对实复数基带信号做换算，导致实际噪声功率比预期大一倍；二是仿真统计BER时把导频子载波、循环前缀的开销都算在了有效比特里，但论文的BER横轴只看数据的有效SNR，两者天然差出几个dB。

解决：把SNR定义写成一行注释放在加噪代码旁边：“SNR(dB) = 10log10(信号功率/噪声功率)，信号功率按发射波形平均功率计算”。BER统计时只统计data_sym位置上的比特，导频位置不参与统计。如果论文用Eb/N0，还要根据每比特能量和码率换算，换算公式写在metrics/ber_curve.m顶部。

5.3 仿真跑一次要一小时：先固定随机种子再减小仿真规模

现象：蒙特卡洛循环里每次调用随机函数，运行结果每次都不一样，跑一个SNR点要几分钟，整条BER曲线要好几个小时。更糟的是发现代码有bug，改完又得重跑。

原因：每次运行都生成新的随机信道和噪声，没有固定随机种子；帧数设置太大，没有先做小规模验证。

解决：在main.m开头加rng(2024)固定随机种子，这样同一台机器上的调试结果可复现。验证功能时把Nofdm降到8、蒙特卡洛次数降到20次跑通流程，确认逻辑正确后再把参数恢复到论文设定跑正式数据。另外把SNR循环用parfor并行化，能明显缩短时间，但要保证每个worker的工作区变量显式传递，否则并行循环里引用不到主工作区变量会直接报错。

5.4 距离估计偏移了几百米：索引换算和时延模型二选一

现象：无噪声条件下，目标设在50米，range_est却输出200多米，甚至出现负值。

原因：索引换算公式里忘了减1，或者目标信道模型里计算回波时延用的是采样点数取整，但取整误差在FFT峰值上变成了多bin偏移；还有一种情况是速度高时多普勒相位跨符号累加错误，导致距离维和多普勒维耦合。

解决：先用数字换算确认理论值，第3.2节已经算过目标对应距离bin大约在33附近。如果代码输出在33附近，说明换算无误；不在就检查channel/target_channel.m里时延是不是按R*2/c0再乘以fs换算成采样点，且这个值是否被正确引入发射波形。注意时延不一定是采样点整数倍，回波模型里要做小数时延插值，否则会有恒定偏差，这个偏差在论文的RMSE曲线里会暴露。

5.5 论文数值和代码数值对不上：先对齐仿真条件再看代码

现象：所有参数都和论文参数表一致，但自己的检测概率或BER明显劣于论文，怎么调都补不上那几dB。

原因：论文仿真里可能用了你没有注意到的先验条件，比如理想信道估计、已知目标存在、去掉了导频开销、用了更好的检测器；也可能论文的条件本身就不完整，补全后和你实现的设定存在差异。

解决：不要假定论文全对或全错，把论文里“仿真条件”单独列一节，和论文对比时逐个打勾：天线数、信道模型（AWGN还是多径）、信道信息是否理想、门限怎么设、蒙特卡洛次数。逐项对齐后如果还有1-3dB差距，通常是实现细节差异，写进论文的“对比分析”段落即可；如果差距超过5dB，优先怀疑SNR定义，其次才怀疑算法实现。复现的目的不是一比一复刻论文，而是证明你掌握了该方法的实现路径。

6. 验证与进阶：让ISAC复现结果经得起答辩追问

复现的最后一个环节是验证，很多人跑出一张图就算完了，结果答辩时一问细节就露馅。我现在收尾前习惯做三个自检，耗时不超过半小时，但能把结果从“看起来像”变成“确实对”。

第一个自检：无噪声条件下，距离和速度估计值必须和目标设定值在各自分辨率以内，偏差大于半个分辨率就要回头查代码。第二个自检：SNR从0到30dB跑一遍检测概率和BER，检测概率曲线应呈平滑S型，BER在高SNR区域应逼近理论QAM曲线，如果误差大于2dB先查SNR定义。第三个自检：把rd_map的峰值位置和公式预估值对比，峰值对应的距离和速度应该和设定目标一致，不一致就回头检查回波模型。

如果三个自检都过了，可以往答辩加分方向扩展。最简单的进阶是把单目标换成多目标，峰值查找从max函数换成二维CFAR检测器；想做通信深度的可以把AWGN信道换成多径信道，在信道模块里加一个瑞利衰落模型；想做系统级的可以把参数扩展到MIMO天线，让距离和速度估计同时用到空间维。这些扩展方向都可以在现有工程目录里加一个模块实现，不影响已经完成的部分。

6.1 三个自检信号的具体做法

无噪声自检最简单：把awgn_channel.m里的噪声功率设成0，跑一次全链路，打印range_est和speed_est，和target_range、target_vel对比。蒙特卡洛自检要在metrics/det_prob.m里写循环统计：做200次带噪仿真，统计估计误差小于分辨率一半的次数占比，这个占比就是检测概率的近似值。理论曲线对比则是把仿真BER和qfunc算出的理论误码率画在同一张图上，注意换算成同样的SNR定义。

6.2 答辩级进阶：从单目标到多目标与MIMO

多目标扩展时，把rd_map的峰值搜索改成循环找局部极大值，再配合一个固定门限滤掉旁瓣。MIMO扩展则是在发射端增加天线维度，接收端每个天线都做一次二维FFT，再用波束成形合并各天线的rd_map，信噪比增益会随天线数线性增加。两个方向都有大量公开的同类实现思路可参考，关键是先保住现有代码能跑，再谈扩展。

我现在的习惯是：论文复现题做到“能复现论文一个核心结果加一个自己扩展的对比实验”就够了，前者证明读懂了，后者证明会改了。ISAC通感一体化的仿真链路并不复杂，复杂的是把每一步的物理意义和数值对应关系理清。希望这篇笔记能帮你少走几段弯路，也希望每一个复现卡住的时刻，都能靠这一套排查思路快速定位问题。

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