★ FEATURED ARTICLE 简介这是一套面向可见光通信系统正交频分复用OFDM仿真学习的MATLAB源码包覆盖从信号产生、调制解调到性能评估的完整流程适合通信工程、光通信方向的学生、课程设计者及科研人员入门使用。包内共5个m文件压缩包仅5KB以主程序为核心包含16-QAM与QPSK调制解调实现及升余弦窗函数设计结构紧凑、模块划分清晰便于按功能模块逐一研读和二次修改。运行主程序后可直观得到OFDM信号的时域波形图、调制星座图以及误码率BER曲线等关键图形有助于深入理解OFDM在可见光信道下的抗多径衰落能力和频谱效率优势并通过修改子载波数、调制方式等参数对比不同条件下的系统性能差异。目前已有1357人学习下载虽然代码体积小但完整度较高是快速上手可见光通信OFDM仿真的实用参考资料尤其适合用作课程项目或毕业设计的起点。1. 可见光通信OFDM仿真一份能直接跑出波形、星座和误码率的代码包手头这份可见光通信OFDM仿真资源把完整的收发链路压缩在五个.m文件里OFDMdemo.m负责参数设置和流程调度qam16.m和qpsk.m分别实现16-QAM、QPSK两种调制demoduqam16.m做16-QAM解调判决rcoswindow.m处理升余弦窗函数。跑通一次之后波形图、星座图、误码率曲线三张关键图一次性全出来实验报告里最常要的东西都齐了。这套资源真正解决的痛点是理论看了几章但代码落不了地OFDM的IFFT和循环前缀该按什么顺序处理LED信道上双极性信号怎么变成单极性16-QAM解调的判决边界设在哪这些细节不亲手跑一遍很难有手感。它适合刚入门需要一份干净模板把全链路跑通的学生也适合正在调仿真但星座图总是打不开的工程师。文件少反而是优势能一行行看完把每个环节和理论公式一一对上。2. OFDM在可见光链路里的实现逻辑从五个.m文件拆出完整信号链2.1 为什么可见光链路要用OFDM带宽受限与多径反射两个硬约束LED调制带宽通常只有几MHz到几十MHz直接拿高速非归零码去调制高频分量正好落在LED频响的滚降区里波形失真和码间串扰一起来。OFDM的办法是把一路高速数据拆成几十路低速并行数据分别调制到不同子载波上。每路子载波占的带宽很窄在这个窄带区间内信道频响近似平坦接收端只需要对每个子载波配一个复数均衡系数就能恢复信号这是OFDM在可见光通信里最大的价值。另一个约束来自室内环境。可见光链路里除了直视路径还有墙壁、桌面、天花板产生的反射路径多径时延扩展通常在几十纳秒量级单载波系统做高速传输时很难处理这种拖尾。OFDM把符号周期拉到子载波间隔的倒数那么长再在符号前面加一段循环前缀CP反射路径产生的拖尾信号就被CP吸收掉了接收端丢弃CP长度内的样本后多径问题解决了一大半。VLC的OFDM主流实现分成两个方向DCO-OFDM加直流偏置把双极性信号抬成单极性ACO-OFDM利用奇偶子载波的不对称结构只传正半周信息。从这份代码包的文件构成看主链路走的是规整的IFFT加循环前缀结构rcoswindow.m处理时域加窗OFDMdemo.m里有完整的参数区和信道模拟整体更接近DCO-OFDM的收发模型。具体细节打开主程序参数区一眼就能确认。2.2 代码包文件地图五个.m文件各管哪一段信号链文件名在链路中的职责需要关注的关键点OFDMdemo.m主程序设置子载波数、循环前缀长度、调制阶数、信噪比区间三张仿真图、BER统计结果qam16.m16-QAM发射端映射比特流向复数星座点星座映射表、平均符号能量qpsk.mQPSK发射端映射2比特向四个相位状态星座点相位、能量归一化demoduqam16.m16-QAM接收端解调信道均衡加最小距离判决判决边界、误码计数rcoswindow.m发射端时域加窗抑制OFDM频谱泄漏和带外辐射滚降系数、窗长覆盖范围调用关系很清晰OFDMdemo.m在发射端按调制阶数选择qam16.m或qpsk.m在接收端调用demoduqam16.mrcoswindow.m加在IFFT和循环前缀拼接完成之后对时域符号做边缘渐变处理。整条链路的信号走向是随机比特、星座映射、串并转换、IFFT、加循环前缀、加窗、加直流偏置、光信道、光电检测、去直流、去CP、FFT、均衡、解映射、BER统计。我一般会在OFDMdemo.m里先画出信号链每一级的变量尺寸确认IFFT输出长度是N_sub加CP后长度变成N_subcp_len再去逐个函数里看处理细节。这个习惯能避免很多变量对不上号的问题尤其是文件多了之后行列向量方向不一致是最常见的翻车点。2.3 主循环流程从随机比特到误码统计的完整骨架% OFDMdemo.m 主循环节选SNR扫描 每帧收发 误码累计 clear; close all; % 参数区改参数之前先想清楚它影响的是哪张图 N_sub 64; % 子载波总数控制频谱资源划分颗粒度 N_data 48; % 实际承载数据的子载波数其余留给导频和空载波 cp_len 16; % 循环前缀长度必须大于信道时延扩展的样本数 M 16; % 调制阶数M16 对应 16-QAMM4 对应 QPSK N_frm 1000; % 每信噪比点下的仿真帧数帧数越少曲线抖动越大 snr_dB 0:2:20; % 信噪比扫描区间步进太大曲线不光滑 for snr snr_dB total_err 0; total_bit 0; for f 1:N_frm % 发射端随机比特 - 星座映射 bits randi([0 1], N_data*log2(M), 1); sym_tx qam16(bits); % M16 时调用 qam16.m % 频域星座点放进 N_sub 个子载波其余补零 ifft_in zeros(N_sub, 1); ifft_in(1:N_data) sym_tx; % IFFT 转时域MATLAB 的 ifft 自带 1/N 归一化 ifft_out ifft(ifft_in, N_sub); % 加循环前缀拷贝 IFFT 输出末尾 cp_len 个样本到最前面 tx_cp [ifft_out(end-cp_len1:end); ifft_out]; % 升余弦窗平滑符号边缘抑制频谱泄漏 tx_w rcoswindow(tx_cp.).; % 加直流偏置模拟 LED 单极性驱动系数选取见避坑章节 dc_offset 1.2 * abs(min(tx_w(:))); tx_led tx_w dc_offset; % 信道先统计实际信号功率再按当前 SNR 加高斯白噪声 sig_pwr mean(tx_led.^2); n0 sig_pwr / (10^(snr/10)); rx_led tx_led sqrt(n0)*randn(size(tx_led)); % 接收端去掉直流偏置丢 CPFFT 回频域 rx_ac rx_led - dc_offset; rx_cp_off rx_ac(cp_len1:end); sym_rx fft(rx_cp_off, N_sub); % 取数据子载波交给解调函数统计误码 bits_hat demoduqam16(sym_rx(1:N_data).); total_err total_err sum(bits ~ bits_hat.); total_bit total_bit numel(bits); end ber(snr_idx) total_err / total_bit; end % 绘图部分波形图、星座图、BER 曲线见第 5 章链路拆成发射端、信道、接收端三个区块。发射端的关键是ifft_in的构造先给数据符号分配子载波位置剩下的位置补零这一步决定了哪些子载波真正承载信息接收端的关键是去CP和FFT的顺序顺序反了会导致数据全部错位。IFFT输出是双极性复数信号但可见光信道只能传正实数强度所以加窗之后补了一个DC偏置把信号抬到正区间这是DCO-OFDM的标准处理方式。参数方面N_sub64时每个OFDM符号有64个子载波CP开销是16/(6416)20%这个比例在VLC链路里比较常见。cp_len不能乱改它必须大于信道时延扩展对应的采样点数改小星座图尾部会散改大虽然更抗多径但有用功率被浪费在多余样本上BER曲线会引入一个固定的能量损失。N_frm控制统计稳定性帧数少于几百时BER曲线会出现明显抖动高信噪比区间尤其明显。注意MATLAB的ifft自带1/N归一化如果代码里IFFT之后额外乘了sqrt(N_sub)功率就要按比例补回来否则同一个SNR值对应的噪声功率和理论曲线会整体错位。很多BER曲线偏移的根源都在这表面像玄学其实就是归一化约定没对齐。3. 16-QAM与QPSK调制解调实现参数设置、判决边界和适用场景3.1 16-QAM的映射规则与demoduqam16.m的判决逻辑16-QAM每个符号携带4比特信息前两位映射到实部同相分量后两位映射到虚部正交分量实部和虚部各取±1、±3四个电平总共16个星座点。映射排列方式很关键相邻星座点之间只差1个比特这是格雷码的约定。当噪声把星座点推到相邻位置时误码只翻转1比特而不是2到4比特同样的符号错误率下BER能明显改善。demoduqam16.m的硬判决就是按这个网格边界来的。判决的第一步是对接收符号做能量归一化。标准16-QAM星座的平均符号能量是10所以接收端要把符号除以sqrt(10)缩放到±1、±3的标准网格上。然后对实部和虚部分别做四电平判决判决边界固定在-2、0、2三个位置实部大于2判3在0到2之间判1在-2到0之间判-1小于-2判-3虚部同理。% demoduqam16.m 硬判决核心逻辑 function bits_hat demoduqam16(sym_rx) % 前置信道均衡省略仿真场景通常乘一个复数系数补偿频响 sym_norm sym_rx / sqrt(10); % 归一化到标准16-QAM星座能量 I real(sym_norm); Q imag(sym_norm); % 四电平判决判决边界固定在 -2、0、2 I_hat zeros(size(I)); I_hat(I 2) 3; I_hat(I 0 I 2) 1; I_hat(I -2 I 0) -1; I_hat(I -2) -3; % Q 分量判决同 I % 组合星座点索引查格雷码映射表还原 4 比特 end判决逻辑本身不复杂但有两个坑会在实际操作中冒出来。第一归一化系数必须和发射端星座能量对齐发射端如果用qam16.m做了额外能量缩放这里除以sqrt(10)就不够了要先算发射星座的实际平均能量。第二硬判决在AWGN信道下够用但如果链路里加了LED削波模型外层星座点会被压缩向内移动固定的±2、±4边界会误判很多外侧点这时得先做自适应缩放或改用软判决。3.2 QPSK的相位映射与抗干扰特性QPSK每次携带2比特把比特对映射到四个相位状态00、01、10、11对应复平面上的1j、-1j、-1-j、1-j四个点幅度恒为1。调制实现很简洁本质上就是查表选复数点。解调更直接接收符号落在哪个象限就判成对应的比特对幅度方向的畸变基本不影响判决结果因为四个星座点只靠相位区分。% qpsk.m 调制映射逻辑 function sym qpsk(bits) % bits 长度必须是 2 的倍数 table [ 11i, -11i, -1-1i, 1-1i ] / sqrt(2); n_sym length(bits) / 2; sym zeros(1, n_sym); for k 1:n_sym idx bits(2*k-1) * 2 bits(2*k) 1; sym(k) table(idx); end end除以sqrt(2)是为了把平均符号能量归一化到1方便后续加噪时直接用信号功率和SNR的关系。调制方式的切换只要改OFDMdemo.m里的M参数就行M4走qpsk.mM16走qam16.m。两个函数的接口都是输入比特、输出复数符号主程序不需要大改。这种接口设计也适合做对比实验同一个信道模型分别跑两档调制直接对比两条BER曲线。16-QAM和QPSK的性能差距是工程选型时要算清楚的一笔账16-QAM需要区分4个电平幅度LED的非线性增益压缩、直流偏置漂移都会破坏幅度信息QPSK完全依赖相位即使信号幅度被非线性压缩只要星座点还落在正确的象限内判决就不会错。代价是每符号只传2比特同一信噪比下要达到相同误码率16-QAM大约需要多4到5dB的信噪比余量。3.3 两种调制阶数怎么选看信道余量而不是看参数表选16-QAM还是QPSK本质是拿链路预算换吞吐。室内固定位置接收、直视路径占主导、链路余量充足用16-QAM甚至64-QAM把带宽跑满接收端是移动设备、链路里存在强反射或遮挡、LED驱动电路动态范围有限QPSK更稳。我自己接过一个模拟项目天花板高度和灯具间距定了之后接收面靠边角位置的信噪比比正下方低不少16-QAM在这些位置直接翻车切成QPSK后整个覆盖区域才都能用。要快速验证调制阶数的影响把OFDMdemo.m跑两遍就行第一次M4、第二次M16固定N_sub64、cp_len16、SNR从0到20步进2每条曲线各跑1000帧。把两条BER曲线画在同一张图上观察间距和交叉区域。如果16-QAM在高SNR区间性能和理论曲线贴合很好链路余量就是够的如果两条曲线在多数SNR点上都差不多差说明链路的工作点本身就不适合高阶调制。4. OFDM仿真常见问题排查四个容易让星座图翻车的现场4.1 循环前缀加错位置星座图发散的第一嫌疑现象BER曲线在高信噪比下居高不下星座图画出来不是聚类而是呈云团状每个星座点都发散成一片而且帧数增加也不会改善。原因大多数情况是CP的处理顺序出了问题。循环前缀一定是对IFFT输出的时域信号做拼接把符号末尾的cp_len个样本拷贝到最前面而不是在频域数据上做任何操作。有人图省事把数据子载波拼接后直接加CP再一起做IFFT相当于把频域子载波的位置整体打乱重排IFFT输出的子载波正交性被破坏星座点必然发散。另一种常见错误是接收端没有先丢CP就直接FFT符号定时偏移让FFT窗口跨到两个符号边界上同样表现为星座图发散。解决发射端严格按tx_cp[ifft_out(end-cp_len1:end); ifft_out]拼接接收端在FFT之前执行rx_ac(cp_len1:end)丢弃CP。如果用矩阵批量处理务必确认是按行还是按列拼接行向量和列向量混用是MATLAB里最容易翻车的地方。验证方法很简单把CP去掉跑一次如果BER反而变好说明加CP的代码一定写错了。错误写法问题ifft([cp_data; data])对频域数据加CP破坏子载波映射接收端不丢CP直接FFT定时偏移FFT窗口错位行向量与列向量混用拼接维度颠倒符号错乱4.2 Eb/N0与Es/N0混用BER曲线对不齐理论值的根源现象仿真BER曲线形状和理论曲线一致但整条曲线在横轴上平移了2到3dB甚至更多而且平移量随调制阶数变化16-QAM和QPSK的偏移量还不一样。原因Es和Eb被混用了。OFDM系统里一个符号携带log2(M)个比特符号能量Es和比特能量Eb的关系是EsEblog2(M)换算到dB就是10log10(log2(M))。16-QAM是每符号4比特换算差约6dBQPSK是2比特差约3dB。如果发射端加了循环前缀和数据子载波占位实际有效符号能量还要再乘一个占比系数。很多人直接从星座点功率算SNR没把CP开销和空子载波算进去结果噪声功率对不上。解决统一用Es/N0作为仿真变量理论曲线也换算成Es/N0或者反过来全部用Eb/N0。我一般会在发射端IFFT之后加一行信号功率统计以实际时域信号功率为基准加噪声这样CP和补零的功率损失已经自动包含在信号功率里曲线对不齐的问题基本消失。另外要确认噪声生成方式复数高斯噪声用sqrt(N0/2)(randn1irandn)生成时总功率是N0写成sqrt(N0)会在所有SNR点上引入系统性偏差。4.3 升余弦窗滚降过大BER曲线中段的平台期现象低信噪比时曲线还算正常到了中高信噪比突然出现平台期BER停在10^-3或10^-4附近不再下降继续提高SNR也没用。原因rcoswindow.m的滚降系数设得过大。升余弦窗的完整时域长度是(1α)倍的符号周期α是滚降系数多出的部分要靠符号边缘外自然延伸。如果代码把窗强制截断在一个OFDM符号长度内α越大边缘滚降区就越往符号中心挤压中间的平台区缩短甚至消失每个子载波的时域波形被窗函数切边子载波之间的正交性被破坏。这种破坏是确定性的自干扰不是噪声增大SNR自然救不回来。解决滚降系数控制在0到0.25之间且窗长要覆盖整个OFDM符号周期。α0是矩形窗频谱泄漏大但正交性无损α在0.1到0.25之间可以在保留正交性的前提下显著压低带外辐射α超过0.5基本只适合频谱观察不适合做数据传输。如果确实需要强滚降正确做法是扩大FFT点数做过采样在过采样域里加窗而不是在符号长度内硬截断。4.4 LED削波与直流偏置不当星座图偏移加外点内缩现象星座图整体偏向一侧外侧星座点向里收缩BER高信噪比下也压不下去而且发射功率越大症状越明显。原因LED只响应正向电流负半周信号会被硬性削掉。OFDM时域信号近似高斯分布负峰幅度经常远大于均值直接把双极性信号接到LED驱动上负峰削平、信息丢失。加了直流偏置但量不够负峰仍然会碰到零点被削偏置过高则正峰超出LED线性区间产生限幅失真。两种情况的共同结果都是非线性幅度压缩16-QAM对外层星座点幅度最敏感所以外侧四个点最先出错。解决先跑一次不带削波的仿真统计时域信号的最小值和最大值直流偏置设置为1.1到1.2倍的abs(min(x))同时检查偏置加正峰值是否在LED线性区上限内。如果两头顶不住就要从降低PAPR入手预留子载波、削峰或改用ACO-OFDM方案。ACO-OFDM只在奇数子载波上放数据IFFT后信号天然满足反对称结构负半周的信息不重复可以直接裁掉负半周而不损失数据代价是频谱效率减半。5. 三张性能图怎么读波形图看PAPR、星座图看信道损伤、BER曲线看系统边界5.1 波形图先算峰均比再决定系统余量OFDM时域波形最直观的特征是幅度像高斯噪声包络峰谷落差大。这个特征对LED链路意味着两件事第一峰值功率要落在驱动电路线性区之内第二峰均比直接决定直流偏置要抬高多少偏置抬高又会增加光功率消耗。所以拿到波形图先算PAPR再决定后面所有参数怎么调。% 计算 OFDM 时域符号的峰均比 x tx_cp; % 取单个完整 OFDM 符号含循环前缀 papr 10*log10(max(abs(x).^2) / mean(abs(x).^2)); fprintf(PAPR %.2f dB\n, papr);OFDM的PAPR在10到13dB是常见区间超过14dB就要警惕。仿真里PAPR高不一定会立刻毁掉BER因为高斯噪声场景下大峰值出现概率低但真实LED驱动里峰值直接对应削波失真这个概率乘以运行时间就是实实在在的误码来源。降低PAPR的常用手段是削峰和加窗rcoswindow.m对抑制带外辐射有帮助但对PAPR的改善有限想专门压PAPR要找预留子载波或削峰类的技术。5.2 星座图从散点形态判断信道损伤类型星座图是诊断信道损伤最快的工具比看BER曲线直观得多。把接收端均衡后的星座点画出来不同损伤的特征完全不一样。星座图特征可能原因优先排查方向散点呈均匀高斯云团加性噪声主导系统本身正常无需处理散点整体旋转载波频率偏差或采样钟偏差检查有无频偏补偿、采样率是否匹配外侧点向原点收缩非线性压缩或限幅检查LED削波、直流偏置、PAPR单个星座点分裂成双簇符号定时偏移或CP过短检查丢CP位置及时延估计散点偏到某个象限I/Q失衡或直流偏置未去净检查实虚部增益、去直流环节如果星座图只有外侧那几个点散得明显内侧点还保持清晰说明问题多半出在非线性而不是噪声。噪声对所有星座点的影响是均匀的而限幅只削掉大峰值的样本最大幅度的外圈点受害最重。这个判断规则能帮你快速把问题从信道差定位到发射端非线性。5.3 误码率曲线和理论值对比之前先确认三条基础条件BER曲线是最后一道验证但验证之前必须先满足三条基础条件单调下降、无平台期、与理论曲线在高SNR端收敛。如果其中一条不满足说明前面某级信号链还有问题画再多图也掩盖不了。理论值对比用MATLAB自带的berawgn函数最方便在相同Eb/N0轴上叠加即可。% 在 Eb/N0 轴上对比 16-QAM 和 QPSK 的 AWGN 理论误码率 ebno_db 0:0.5:15; ber_16qam berawgn(ebno_db, qam, 16); ber_qpsk berawgn(ebno_db, psk, 4, nondiff); semilogy(ebno_db, ber_16qam, k--, ebno_db, ber_qpsk, b-);berawgn的第一个参数必须是Eb/N0而不是Es/N0画图前先把坐标轴定义确认好。很多项目里BER曲线对不齐问题不在信道仿真而在横坐标定义不统一。仿真低SNR端和理论值存在偏差是正常的帧同步、均衡估计、CP开销都会在低信噪比时扩大损伤高SNR端如果始终差好几个dB就要回头怀疑某个固定损耗没补偿比如导频子载波占据的功率没有扣除。6. 把OFDM仿真参数调出花样从基线记录到单变量扫描拿到代码包先不要急着改参数。我自己的标准流程是先原封不动跑一遍记录默认参数下的波形PAPR、星座图形态、BER曲线形状这三张图就是后续所有改动的基线。然后一次只改一个参数改完再跑对比和基线的差异。如果一次改三个参数曲线坏了根本不知道是哪一步引入的。参数改动方向优先验证哪张图N_sub 子载波数增大后频谱效率提升但PAPR上升波形图cp_len 循环前缀长度减小时多径抗性下降星座图M 调制阶数4改16后吞吐翻倍链路余量要求变高BER曲线滚降系数增大后带外抑制变好超过0.25破坏正交性波形图 BER曲线直流偏置抬高后削波减少光功率开销增加波形图 星座图N_frm 仿真帧数增大后曲线更平滑耗时线性增加BER曲线把OFDMdemo.m改造成函数接口之后可以做参数扫描for循环遍历cp_len或N_sub每组参数算一个BER值最后把几条曲线叠加到同一张图上。这是对比实验报告里最常用也最稳妥的做法。函数化的额外好处是避免工作区变量污染仿真脚本跑多了旧变量残留经常会让新参数看着改了实际没生效函数化之后这个问题自动消失。还有一点如果只是教学演示、不需要精确链路预算N_sub取64、cp_len取16、1000帧就够用了但如果要把结果写进项目报告帧数建议加到5000以上低误码率区间的曲线才有统计意义。我早期做仿真就是图快用300帧高SNR端的BER曲线抖得像心电图后来改成5000帧才总算把平台期的位置看准确。从那以后我每次拿到新的OFDM仿真代码都会强制走一遍基线记录和单变量扫描的流程再开始调参数希望帮到你。本文还有配套的精品资源点击获取
阅读完成 · 觉得有帮助？
