简介一套基于 MATLAB 的中国地面数字多媒体广播DTMB仿真代码包面向数字电视与通信系统学习者、研究人员用于理解 DTMB 工作原理、完成信号处理算法验证和性能优化实验。压缩包共 12 个文件全部为 .m 脚本整体约 16KB覆盖 TDS-OFDM 时频变换、发射机调制、QAM 调制、循环前缀与信道估计、均衡算法、PN 序列生成、相关性计算及频率扫描测试等模块便于按环节分步调试。已有 886 人学习/下载对入门地面数字电视仿真的开发者有一定参考价值。运行这些代码可模拟 DTMB 信号从发射、传播到接收的全过程观察多径衰落与信道噪声的影响并通过调整编码方式、调制阶数等参数分析误码率与抗干扰能力为工程应用或课题研究提供可修改、可扩展的实验基础。1. DTMB 仿真这套 MATLAB 代码到底能跑通什么上手这套 DTMB 仿真代码前我建议你先放下「找个主函数直接跑」的念头。这份资源对应的是中国地面数字电视标准里最核心的 TDS-OFDM 收发链路不是把完整发射机做成一个大黑匣子而是把 QAM 映射、PN 序列生成、OFDM 调制、信道估计、均衡分开成十几个 .m 文件适合两类人——刚接触地面数字电视想看懂帧结构的学生以及要做信道估计算法对比的工程师。它能帮你从信号生成一路走到星座图恢复而不再对着国标文档里的 3780 这个数字发愣。先跑 test_CR.m再回头看发射机这是我认为最不容易翻车的打开方式。2. 认识这套代码从文件分工看懂收发链路2.1 DTMB 为什么用 TDS-OFDMPN 序列不只是保护间隔传统 OFDM 系统的保护间隔是用循环前缀CP实现的接收端丢掉 CP 就能消除符号间干扰但 CP 里的数据是重复的不携带任何额外信息。DTMB 国标做了一个看起来很「抠频谱」的设计保护间隔里不填重复数据而是填一段 PN 序列这就是时域同步 OFDMTDS-OFDM名字的由来。PN 序列在这个系统里承担三个任务一是作为保护间隔对抗多径时延二是作为同步导频让接收端找到帧起始位置三是作为已知序列做信道估计。也就是说一个 PN 头同时干掉了 CP-OFDM 里三样东西的活。代价是接收端算法更复杂——你不能像 CP-OFDM 那样直接丢掉保护间隔得把 PN 头从接收信号里分离出来做处理这也是这套代码里 correlation.m、CR_loop_base.m 存在的原因。帧结构上DTMB 的一帧由帧头和帧体组成。帧体是 3780 个数据载波帧头 PN420 是 420 个 PN 符号整帧 4200 个符号。3780 正好是 420 的 9 倍所以一个 PN 头能覆盖帧体 1/9 的时延范围这个比例在仿真里直接决定你配置多径信道时最大时延能放到多大。2.2 文件与模块12 个 .m 文件的职责划分拿到压缩包后第一件事不是看代码而是对照下表核对文件缺一个文件后续跑起来报错会非常难定位文件名链路位置主要作用常见坑Tx_Modulator.m发射总入口调用后续模块生成完整基带信号参数在文件头统一改Qammod_CR.m发射端QAM 映射与解映射归一化功率别漏PN_255_Gen.m发射端生成底层 m 序列生成的是 255 长度别直接当帧头用PN420.m发射端由 PN255 扩展成 420 帧头序列帧头模式切换靠它Tds_ofdm_CR.m发射端帧体 IFFT 与帧头拼接输出长度必须是 4200correlation.m / correlation_s.m接收端滑动相关找帧起始慢速版用于验证快速版用于循环CR_loop_base.m接收端基础信道估计与恢复单径无噪声下应先跑通CR_loop_methods.m接收端多种信道估计方法切换方法名改错会静默降级NAFC_CAFC_TO.m接收端频域均衡与定时/频偏校正参数敏感先固定扫频Freq_Sweep_test.m测试频率扫描评估性能跑完看 BER 谷底在哪test_CR.m测试信道恢复算法入口测试建议从它开始跑表里有两类文件值得先读一类是带 _CR 后缀的CR 在这里指 Channel Recovery也就是信道恢复相关模块一类是带 _s 后缀的通常是简化实现。我习惯先读 correlation.m 和 correlation_s.m 的差异能看出作者对性能和速度的取舍也能猜出接收链路的主体流程。2.3 运行前的环境检查先把坑踩在门外这套代码是纯 .m 文件主要依赖 MATLAB 基础函数和信号处理相关函数不额外依赖通信工具箱所以版本要求不高。但有几个环境问题很常见建议每次打开 MATLAB 后先执行一遍% 清空工作区变量防止上次残留的变量名影响本次运行 clear; clc; close all; % 把代码目录加入搜索路径避免脚本里互相调用的文件找不到 addpath(你的代码目录); % 把路径换成实际位置 % 设置随机种子保证每次复现的噪声序列一致 rng(42);clear 和 clc 这两行看起来多余但仿真代码最玄学的报错之一就是工作区里残留的变量覆盖了函数内部同名变量。我见过有人跑 test_CR.m 一直报维度不匹配排查半天发现是上一次调试留下的全局变量没清。addpath 是防止你把代码放在子目录里却忘了把路径加进 MATLAB 搜索路径。rng 设置随机种子则是复现实验的基本素养特别是涉及噪声和多径信道时没有固定种子你很难判断结果是算法改进了还是运气变好了。提示如果运行时报 undefined function先检查文件名大小写是否与调用处一致。MATLAB 在 Windows 下文件名大小写不敏感但在 Linux 或 macOS 下敏感跨平台拷贝代码时最容易踩这个坑。3. 把发射链路跑通从 QAM 映射到 TDS-OFDM 组帧3.1 发射信号流从随机比特到基带帧Tx_Modulator.m 是发射端的入口。我拿到它的第一件事是读文件头部的参数区那里一般会定义调制阶数、载波数、PN 长度这些核心参数。常见做法的信号流是这样随机比特流先做 QAM 映射成复数符号按 3780 个载波排列成帧体对帧体做 IFFT 得到时域波形再在前面拼接 PN420 帧头最后组成完整的 TDS-OFDM 帧。这个流程对应到代码里是四步缺一步帧结构就断了% 发射链路四个核心步骤简化示意 bits randi([0 1], 3780 * 4, 1); % 生成 3780×4 bit对应 QAM16 一帧 qam_symbols Qammod_CR(bits, 16); % QAM16 映射输出 3780 个复数符号 frame_body reshape(qam_symbols, 3780, []); % 按 3780 载波排列成帧体 pn_head PN420(); % 生成 420 点帧头序列 tx_frame [pn_head; frame_body]; % 拼接整帧 4200 点第一行里 3780 × 4 的 4 是 QAM16 的比特数换成 QAM64 就要改成 3780 × 6这是改调制阶数时最容易漏的地方。Qammod_CR(bits, 16) 的第二个参数是调制阶数 M不是 log2(M)传错的话星座点数量全乱。reshape 的作用是把符号流按 3780 一行排成矩阵每一列是一帧的频域数据如果用 reshape(qam_symbols, [], 3780)排列方向就反了后面 IFFT 出来的波形完全是乱的。最后拼帧用 [pn_head; frame_body]PN 头在前、帧体在后顺序反了同步就直接失败。3.2 QAM 映射的手写实现为什么不用工具箱函数Qammod_CR.m 没有用通信工具箱的 qammod这是这类仿真代码的常见选择。原因有两个一是纯函数便于在无工具箱的 MATLAB 环境运行二是手写映射可以控制星座点的排列顺序和功率归一化这对后续信道估计的幅度一致性很重要。QAM 映射的核心是查表16QAM 的星座点按 4bit 分成 I、Q 两路各 2bit映射到 ±1、±3 四个电平上。但注意不同标准里格雷码的排列方式不一样解映射的顺序必须和映射完全一致否则误码率会高得离谱。% 16QAM 映射核心逻辑示意 function sym Qammod_CR(bits, M) k log2(M); % k4每个符号 4bit bits reshape(bits, k, []); % 每列是一个符号的 4bit % 拆分 I/Q 两路各取 2bit按格雷码映射到 ±1、±3 I (2 * bits(1,:) - 1) .* (2 - bits(2,:)); Q (2 * bits(3,:) - 1) .* (2 - bits(4,:)); sym (I 1j * Q) / sqrt(10); % 除以 sqrt(10) 归一化平均功率为 1 end倒数第二行里的 sqrt(10) 是 16QAM 归一化因子16 个星座点的平均功率平方和正好是 10。如果删掉这行星座图整体放大接收端均衡后还要额外乘一个固定增益否则误码率统计会偏。改 QAM64 时归一化因子要重新算64QAM 的平均功率是 42不是想当然的 6 或 8。3.3 Tds_ofdm_CR.m 与 PN420.m 的组帧细节Tds_ofdm_CR.m 做的事情是把 3780 个频域符号做 IFFT 变成时域信号。DTMB 帧体是 3780 个载波这里的关键是 IFFT 点数是 3780 而不是常见的 4096。3780 不是 2 的幂MATLAB 的 ifft 可以处理任意长度但性能会比 4096 慢而且频域索引的计算方式要按 3780 来不能套用 2 的幂次那套 fftshift 习惯。PN420.m 生成帧头序列。从 PN_255_Gen.m 的命名能猜出底层是一个 255 长度的 m 序列PN420 的常见生成方式是把 255 长度 m 序列循环扩展后补位到 420。这个扩展规则在国标文档里有明确定义仿真时直接调用 PN420() 就行不用自己拼。但要提醒的是PN420 的 420 个符号是 BPSK 映射的 ±1 序列不是复数序列拼到 IFFT 输出时要注意数据类型一致性。% TDS-OFDM 组帧验证检查输出结构与星座映射是否一致 tx_frame Tx_Modulator(); % 假设入口无参调用 frame_len 420 3780; % 帧头 帧体 num_frames length(tx_frame) / frame_len; fprintf(帧数: %d, 单帧长度: %d\n, num_frames, frame_len);如果 num_frames 不是整数说明组帧长度算错了常见原因是 PN420 和帧体的拼接顺序或长度不匹配。跑通发射端后我建议先画个星座图看一眼从发射帧里取出帧体位置做 FFT 后解映射回符号如果星座点聚拢在标准 16QAM 位置发射链路就基本没问题后面调试接收端才有个可信的参照物。3.4 首跑建议先用单帧无噪声验证拿到代码别急着直接上多径和噪声。最稳的做法是先把发射链路产出的信号直接丢给接收端中间不加任何信道跑通一个完美信道下的收尾流程。这一步跑不通后面所有信道相关的问题都会混在一起根本分不清是同步的问题还是均衡的问题。% 首跑发射 → 接收不加信道验证全链路基础逻辑 tx_frame Tx_Modulator(); % 这里先跳过信道直接把发射信号喂给接收端 rx_frame tx_frame; % 完美信道 [ber, sym_err] test_CR(rx_frame); % 调用测试入口 fprintf(BER: %.4f\n, ber);test_CR.m 在这套代码里的角色是信道恢复算法的验证入口我一般从它开始而不是从发射机开始原因就在这它能直接告诉你当前算法链路在什么信噪比下能工作省得自己在主流程里翻变量。如果完美信道下 BER 不是 0问题大概率出在 QAM 映射和解映射的对应关系上优先检查星座点顺序和归一化因子。4. 接收端不糊的底气同步、信道估计与均衡4.1 滑动相关怎么从噪声里找到帧头接收端拿到的信号是连续的基带波形你得先找到每一帧从哪里开始。DTMB 的帧头是已知的 PN 序列所以最常见的做法是把接收信号和本地 PN 序列做滑动相关相关峰出现的位置就是帧头位置。这套代码里 correlation.m 和 correlation_s.m 就是干这个的后者通常是前者的简化版用分段相关或隔点采样的方式减少运算量。% 滑动相关找帧起始示意 N length(rx_signal); L length(pn_seq); % 420 corr_out zeros(N - L 1, 1); for idx 1 : N - L 1 % 逐点滑动计算与本地 PN 的内积 corr_out(idx) abs(sum(rx_signal(idx : idx L - 1) .* conj(pn_seq))); end [~, frame_start] max(corr_out); % 峰值即帧起始这个循环逐点算相关速度很慢correlation_s.m 一般会改成用 FFT 或分段累加来加速。峰值位置 frame_start 的含义是从这一点开始的 420 个点是帧头 PN往后 3780 个点是帧体数据。如果多径信道把 PN 头的能量扩散了相关峰可能不是一个尖峰而是一个平坦区域此时取 max 只能覆盖最强径弱径信息要靠信道估计那一步去补。4.2 从帧头提取信道响应CR_loop_base 的核心逻辑找到帧起始后下一个问题是怎么知道信道长什么样。DTMB 的设计巧妙之处在于帧头的 PN 序列是已知的接收端把接收到的 PN 头与本地 PN 做相关就能得到信道的冲激响应。这个过程在代码里就是 CR_loop_base.m 的主循环对每个帧头做相关、截取信道响应、选抽头、把结果送给均衡器。% 信道冲激响应提取示意 rx_pn rx_signal(frame_start : frame_start 419); % 截取接收帧头 cir ifft(fft(rx_pn) ./ fft(pn_seq)); % 频域相除后还原 cir(abs(cir) thr) 0; % 按噪声门限截断抽头第三行的频域相除是信道估计的教科书做法意思是接收 PN 等于发送 PN 卷积信道频域相除等效于去卷积。但直接相除有个问题PN 序列的频谱在某些频点接近 0相除会把噪声放大所以 CR_loop_methods.m 里一般会提供几种替代估计方法比如 DFT 域平滑、最小二乘加阈值截断。阈值 thr 的选取直接影响估计出的抽头数量设大了会丢掉真实多径设小了会留一堆噪声抽头常见做法是先统计噪声底噪再乘一个系数比如 1.5 到 3 倍。4.3 频域均衡与残余误差校正NAFC_CAFC_TO.m 在收什么尾信道响应拿到之后频域均衡本身不复杂把 CIR 补零到 3780 点做 FFT 得到每个子载波上的信道增益帧体数据逐载波除以这个增益就行了。复杂的是残余的定时偏移和频偏——PN 序列相关能把帧起始定到采样点级但符号定时误差会转化为子载波间的相位旋转频偏则会让星座图整体旋转。NAFC_CAFC_TO.m 的名字拆开看NAFC 和 CAFC 是两套均衡策略TO 指定时偏移。常见做法是先用粗同步定位帧起始再做频域均衡最后用导频或判决反馈校正残余相位。% 频域均衡核心三步示意 cfo_est fftshift(...); % 粗略频偏估计 rx_signal rx_signal .* exp(-1j * 2 * pi * cfo_est * (0:length(rx_signal)-1) / fs); cfr fft([cir_est; zeros(3780 - length(cir_est), 1)]); % CIR → CFR eq_symbols fft(frame_body, 3780) ./ cfr; % 帧体逐载波均衡第二步是频偏补偿exp 的指数里 cfo_est 是频偏估计值fs 是采样率如果频偏估计偏差超过子载波间隔的百分之几补偿后星座图还是会转。第三步的 cfr 是信道频率响应注意 CIR 只有几十个抽头补零到 3780 再 FFT 才能和帧体长度对齐。这里最容易翻车的是补零方向CIR 要补在末尾而不是开头否则 FFT 出来的频响相位全错星座图会沿对角线散开。4.4 调试接收链路的顺序先同步、再估计、后均衡接收端三个模块是串联关系调试顺序错了会浪费大量时间。我的习惯是先固定信道为单径无噪声只验证同步模块能不能正确定位帧头然后固定同步正确加一条多径验证信道估计最后再加频偏和噪声调均衡器的参数。% 多径信道参数配置常见做法 chan_taps [1, 0.5, 0.2]; % 三条径主径 两条反射 chan_delays [0, 3, 9]; % 时延单位采样点 rx_noisy filter(chan_taps, 1, tx_frame); % 卷积多径 rx_noisy rx_noisy 0.1 * randn(size(rx_noisy)); % 加噪声时延单位是采样点9 个采样点对应 0.1μs 级别远小于 420 点帧头能覆盖的范围这是合理的首测配置。如果时延超过 420多径会串到下一帧的帧体里导致 PN 头和数据互相污染此时均衡算法基本救不回来需要先调整帧结构或接受误码率上限。5. 避坑与排查这套 DTMB 仿真最常翻车的 5 个地方5.1 星座图转圈频偏没校干净现象均衡后的星座图呈现出同心圆旋转误码率在 0.4 以上看起来像所有符号都在缓慢漂移。原因频偏估计的初值不准确或者直接跳过了频偏补偿直接用均衡器硬扛。OFDM 对频偏非常敏感子载波间隔是符号率的倒数DTMB 帧体 3780 个载波对应约 2kHz 子载波间隔只要频偏超过几十赫兹星座就会开始旋转。解决先用 Freq_Sweep_test.m 对频偏做一次扫描找到误码率最低的频率点作为补偿初值。常见做法是把频偏补偿拆成两步先粗补偿到子载波间隔的 10% 以内再用判决反馈精补偿剩余相位。扫描步长初始可以取 100Hz 量级找到谷底后缩小到 10Hz 量级重扫。5.2 相关峰不突出多径把 PN 能量摊平了现象correlation 输出的峰值和旁瓣差不多高frame_start 每次运行跳来跳去BER 时好时坏。原因多径信道的时延跨度大或者信噪比太低导致相关峰被噪声淹没。PN420 的自相关特性决定了它在低信噪比下的检测能力有限而多径会把相关峰的能量分散到不同时延的峰上。解决不要用单帧找峰用多帧非相干累加后再找峰。也就是把连续几帧的滑动相关结果对应位置相加取绝对值帧头在每一帧都有峰累加后信噪比提升明显。correlation_s.m 如果实现了滑窗累加检查帧数参数是否够大常见做法是 5 到 10 帧累加。5.3 改了调制阶数后误码率飙升归一化因子没跟着换现象QAM16 跑得好好的把参数改成 QAM64 后星座点错位误码率固定在 0.4 附近。原因只改了调制阶数但映射函数里的归一化因子还是 QAM16 的 sqrt(10)。64QAM 的平均功率是 42归一化因子应该是 sqrt(42)用错之后整个星座图的平均幅度变了均衡器输出的判决边界全部错位。解决把归一化因子放进调制阶数的配套参数里不要散落在函数各处。我一般会把 M、k、归一化因子做成一个参数结构体换调制方式只改一处。另外解映射端的判决电平也要按同一次归一化后的星座位置算不要直接用 ±1、±3、±5、±7。5.4 PN420 和 PN255 混用同步和估计全乱现象代码里既看到 PN420.m 又看到 PN_255_Gen.m有人直接拿 PN_255_Gen 的输出当帧头用结果同步相关峰完全消失。原因PN_255_Gen.m 生成的是底层 255 长度的 m 序列PN420.m 是在它基础上做循环扩展和补位得到的两者长度不同、自相关特性不同不能直接互换。255 是帧头模式的底层序列420 才是实际发射的帧头长度。解决发射链路里帧头只能是 PN420PN_255_Gen.m 只作为生成 PN420 的底层工具不直接进入组帧环节。改帧头模式时先确认目标帧头长度与帧体的匹配关系不要跨模式混用。5.5 运行报错 undefined function跨平台拷代码后大小写问题现象在 Windows 上跑得好好的代码拷到 Linux 或 macOS 上运行报某个函数 undefined但文件明明就在目录里。原因MATLAB 在 Windows 下不区分文件名大小写但在 Linux 和 macOS 下严格区分。代码里调用处写成 Correlation.m 而实际文件名是 correlation.m在 Windows 上能跑跨平台就崩。解决拿到代码后先做一个全目录大小写检查直接把调用处和文件名统一成小写是最省事的做法。另一个习惯是每次运行前都用 which 查询一遍关键函数确认 MATLAB 实际加载的是当前目录里的版本而不是搜索路径里其他同名文件。提示如果遇到 which 查出来多个同名文件优先怀疑 MATLAB 搜索路径里残留了旧版本代码目录。清掉多余路径后重新 addpath 当前目录能避免很多「代码没改但结果变了」的诡异问题。6. 进阶验证对齐国标参数测 BER 曲线这套 DTMB.zip 里的仿真代码足够支撑你把一条完整链路跑通但跑通只是起点。真正有价值的是用它对系统性能做量化评估。我建议你做两件事一是把发射端的采样率、符号率等参数对齐到国标标称值二是用 test_CR.m 套一个信噪比循环画出 BER 随 SNR 变化的曲线。对齐参数的技巧是标称符号率 7.56 MSps、帧体 3780 个载波、PN420 帧头仿真里采样率直接用符号率还是过采样会影响频偏单位和多径时延的换算建议统一以采样点为基准写注释别在代码里混用 Hz 和采样点。% 遍历信噪比统计 BER示意 snr_list 0:2:20; % 信噪比扫描范围 ber_list zeros(size(snr_list)); for idx 1:length(snr_list) % 每个信噪比下重发一帧统计误码 ber_list(idx) test_CR(tx_frame, snr, snr_list(idx)); end semilogy(snr_list, max(ber_list, 1e-6)); % 对数坐标画曲线 grid on; xlabel(SNR (dB)); ylabel(BER);max(ber_list, 1e-6) 是为了让零误码点也能在对数坐标里显示不处理的话 semilogy 遇到 0 会画不出点。跑完这条曲线后保存一次工作区变量把不同调制阶数、不同多径配置下的曲线放在同一张图里对比你会直观看到「把 PN 头做得再完美也救不回超过保护间隔的多径」这个结论。从那以后我每次拿到新的仿真包第一件事都是先找入口文件和测试脚本用最小配置跑通一条 BER 曲线基线再开始动算法。这个习惯省下的调试时间远超想象希望帮到你。本文还有配套的精品资源点击获取