简介基于MATLAB的OFDM可见光通信系统仿真代码包面向通信工程、电子信息及相关专业的初学者和研究者旨在解决可见光通信链路从信道建模到OFDM收发端仿真的完整实现问题。资源内含基于朗伯模型的可见光无线信道模块以及OFDM发送端与接收端的全部代码所有关键步骤均配有详细中文注释直接运行即可观察结果。压缩包共10个文件其中9个为.m脚本1个为.mat数据文件整体仅10KB文件按信道求解、调制映射、加窗滤波等环节清晰划分便于逐段学习。已有1973人浏览学习特别适合正在入门OFDM或可见光通信仿真、需要快速搭建系统框架的读者。仔细研读HLOS、PRLOS等信道函数以及QAM/QPSK调制解调代码可深入理解朗伯信道特性、子载波映射、循环前缀和信号恢复流程对课程设计或毕业设计具有直接参考价值。 说来也好笑我最早接触到OFDM可见光通信系统仿真是从一个塞满散文件的.zip压缩包开始的。解压之后里面没有文档只有一堆.m脚本和一张画得歪歪扭扭的BER曲线图参数散落在各个文件里连有没有加直流偏置都要自己翻代码猜。我硬着头皮对着论文和标准把整套系统从头跑通了一遍中间踩了不少坑。如果你正准备搭一套可见光通信的仿真环境或者毕设方向正好是OFDM加VLC那这篇东西大概率能帮你省下两周时间。我会从信号链路、参数怎么定、完整代码怎么写到最终怎么调出能看的误码率曲线一次说清楚。1. 先搞清楚这个仿真系统到底在做一件什么事很多人在拿到这类仿真工程后的第一反应是直接读代码但我建议你先退一步把整个系统的逻辑走一遍。可见光通信VLC和传统射频通信有一个本质区别LED发出的光强度只能是非负实数所以通信方式天然是强度调制加直接检测IM/DD。我们没法像射频那样直接发IQ复数信号必须把OFDM信号处理成实数还要叠加一个直流偏置把它抬到LED的工作区间。1.1 为什么可见光通信偏偏要选OFDMOFDM在可见光里的地位基本就是多径信道里的“防弹衣”。室内场景里灯光经过墙壁、天花板、桌面反射后会形成很多条到达路径这些路径叠加起来会产生频率选择性衰落。如果只用单载波调制一个符号受到前后干扰就废掉了必须上复杂的均衡器。而OFDM把整个高速数据流拆成几十上百个窄带子载波每个子载波带宽很窄信道在每个子载波上近似平坦只需要在每个子载波上做一次单抽头均衡就能恢复数据实现复杂度低很多。另外OFDM天然支持自适应调制也就是说信道状态好的子载波发16QAM甚至64QAM信道差的子载波就退到QPSK或者干脆不发数据。在可见光场景里不同子载波受LED非线性响应和信道衰落影响的程度差异很大这个特性特别实用。我做的这套仿真里也保留了按子载波分配比特的接口后面想加自适应调制算法可以直接扩。1.2 从电气接入到光输出一条完整信号链很多人第一次看VLC发射端信号链会懵其实可以想象成一条流水线。原始比特流先做QAM映射变成一个个复数符号接着进行串并变换把这些符号分配到各子载波上然后做IFFT把频域符号变回时域波形。到这里你会发现IFFT输出一般是复数可LED只能发实数光强所以在IFFT输入这一侧要做共轭对称处理把N个频域符号扩成2N个让IFFT输出变成纯实数。随后加上循环前缀抵抗多径时延扩展最后再加一个直流偏置把负半轴的信号也抬到非负区间送入LED的驱动电路。接收端就是这条流水线的逆向过程光电检测器PD把光强变化转成电流经过跨阻放大器放大后由ADC采样接着去掉循环前缀做FFT回到频域再用LS信道估计或简单的一条导频序列估计信道响应对每个子载波做均衡最后QAM解映射还原成比特和原始比特对比就能统计误码率BER。这套链路看起来简单但每个环节都有坑后面我会逐个展开。2. 系统仿真架构设计与关键参数推导仿真不是把代码堆在一起能跑就行参数怎么定直接决定结果有没有参考意义。我给出的参数不一定适用所有场景但可以作为一套比较合理的地基。2.1 发射端从比特到光强信号的三个关键环节QAM调制这一步优先推荐用QPSK验证链路通断等BER曲线趋势对了再换16QAM。仿真中一般用qammod函数注意UnitAveragePower要置为true否则星座点功率没有归一化后面加噪声的SNR算不准。子载波总数N取256其中实际传输数据的子载波设为200多出来的子载波做频谱保护带这样能避免频带边缘受滤波器滚降影响。共轭对称是VLC中OFDM最重要的操作。做法是把频域向量构造成[0 DC, data, 0, conj(flip(data))]的形式。注意索引0的位置是直流子载波不携带数据Nyquist频率对应的那个点也要置零。IFFT出来之后信号就是实数了。如果不做这一步后面的数模转换和LED驱动根本没法送信号。加直流偏置时有个细节偏置大小应该按信号标准差来设置而不是拍脑袋给个固定值。通常设DC_bias 0.35 * PARAM.signal_std意思是直流分量大约是有效信号幅度的0.35倍方差单位。偏置太小会切掉大量负半轴信号造成严重非线性失真偏置太大又会浪费光功率降低能量效率。LED本身的非线性限幅也要建模我用的max(0, min(x, LED_max))看起来简单但已经能模拟LED进削顶区后BER急剧恶化的现象。2.2 信道与接收端光噪声、低通响应与同步恢复室内可见光信道在仿真里可以用一个抽头延迟线模型来表示我常用的冲激响应是h [1.0, 0.25, 0.12, 0.08, 0.05];这组抽头模拟了直射路径和若干条主要反射路径抽头数量4到5个就足够看出OFDM对抗多径的优势。信道总增益要归一化不然SNR会偏。接收端光电检测器输出的噪声用加性高斯白噪声AWGN建模但在可见光系统里噪声并不完全是加性的主要来源是散粒噪声和热噪声二者在高光功率下近似高斯分布所以一般仿真直接叠加AWGN也可以接受。同步问题是你最容易忽略的地方。时域上需要先找到OFDM符号的起始位置仿真里最简单的做法是用一个已知的导频序列做相关。帧同步完成后再去掉循环前缀就可以直接做FFT。我见过很多新人上来就把接收序列直接丢进FFT结果星座图完全是一团浆糊就是同步半步没对齐导致的。3. 完整仿真实现与复现方法这一部分我直接给出可以跑的MATLAB代码框架大家不要盲目复制粘贴要根据自己的参数稍加调整。3.1 模块化代码结构别把所有变量堆进一个主脚本吃过大亏之后我强烈建议把系统拆成参数定义、发射机、信道、接收机、性能统计五个模块。主脚本只保留流程控制和绘图。你是要做参数扫描实验的比如改变循环前缀长度观察BER变化如果代码全部堆在一起每改一次参数就得重新从头跑一遍特别浪费时间。一个比较清晰的主脚本骨架如下%% 主控制脚本ofdm_vlc_sim_main.m clear; clc; close all; % 参数初始化 params.N_subcarriers 256; params.N_used 200; params.CP_len 32; params.M 4; % QPSK params.frame_len 100; % OFDM符号帧数 params.LED_max 1.8; params.DC_bias 0.35; % 扫SNR点统计BER SNR_dB 0:2:20; BER_array zeros(size(SNR_dB)); for idx 1:length(SNR_dB) BER_array(idx) ofdm_vlc_ber(SNR_dB(idx), params); end figure; semilogy(SNR_dB, BER_array, o-, LineWidth, 1.5); grid on; xlabel(SNR (dB)); ylabel(BER); title(OFDM可见光通信系统仿真 BER 曲线);这里ofdm_vlc_ber就是一次完整仿真的入口函数返回当前SNR下的平均误码率。把扫点逻辑和单点仿真分开之后换成16QAM只是改一个params.M换成27dB以下极限场景也方便。3.2 发射端代码映射、串并变换、共轭对称与DCO发射端核心函数我写成了下面这样注意共轭对称那一行的索引逻辑。function tx_serial ofdm_vlc_tx(bits, params) N params.N_subcarriers; N_used params.N_used; M params.M; n_bits N_used * log2(M) * params.frame_len; % 因为我们要构造实数IFFT实际使用子载波数只能是N_used/2个数据子载波 data_sym qammod(bits, M, InputType, bit, UnitAveragePower, true); data_sym reshape(data_sym, N_used/2, params.frame_len); % 构造共轭对称的IF FT输入 ifft_in zeros(N, params.frame_len); ifft_in(2 : N_used/21, :) data_sym; % 对应位置的共轭对称 ifft_in(N - N_used/2 1 : N, :) conj(flipud(data_sym)); % IFFT变换到时域 tx_time ifft(ifft_in, N, 1) * sqrt(N); % 注意这里乘以sqrt(N)是为了保持功率归一化 % 添加循环前缀 tx_cp [tx_time(end-params.CP_len1 : end, :); tx_time]; tx_serial tx_cp(:).; % 加直流偏置 LED限幅非线性 tx_serial tx_serial params.DC_bias * std(tx_serial); tx_serial max(0, min(tx_serial, params.LED_max)); end几个点值得强调。第一data_sym占用的实际子载波是N_used/2个因为另一半要留给共轭对称。第二直流分量所在的第1个子载波索引位置本身是空的这相当于牺牲了一个子载波但换来时域信号是实数。第三ifft后我乘了sqrt(N)做功率归一化很多教程都不写这一步结果稍后加噪声时OFDM信号功率和QAM符号功率对应不上BER曲线整体偏了。3.3 接收端代码同步、解调与BER统计接收端函数里我加入了理想帧同步的假设但保留了一个粗同步的接口方便以后扩展。function ber ofdm_vlc_ber(snr_dB, params) % 生成发送比特 N_bit_per_frame params.N_used * log2(params.M); tx_bits randi([0 1], N_bit_per_frame * params.frame_len, 1); % 发射 tx_serial ofdm_vlc_tx(tx_bits, params); % 可见光多径信道 h [1.0, 0.25, 0.12, 0.08, 0.05]; h h / sum(h); % 归一化 ch_out filter(h, 1, tx_serial); % 加AWGN tx_power mean(ch_out.^2); noise_power tx_power / (10^(snr_dB/10)); rx_serial ch_out sqrt(noise_power) * randn(size(ch_out)); % 去掉循环前缀并转换为矩阵 N params.N_subcarriers; CP params.CP_len; frame_len params.frame_len; rx_reshaped reshape(rx_serial, N CP, frame_len); rx_no_cp rx_reshaped(CP1 : end, :); % FFT回到频域 rx_fft fft(rx_no_cp, N, 1) / sqrt(N); % 最小二乘信道估计用理想信道响应代替实际系统需要导频 H_est fft([h, zeros(1, N - length(h))], N).; H_est H_est(2 : N_used/21); eq_sym rx_fft(2 : N_used/21, :) ./ H_est; % QAM解映射 rx_bits qamdemod(eq_sym(:), params.M, OutputType, bit); % 统计误码 err sum(rx_bits ~ tx_bits); ber err / length(tx_bits); end这版代码里我用了理想信道估计直接把信道频域响应拿来做均衡。真实系统需要一个帧头处的导频OFDM符号通过LS估计得到H_est这部分可以在后续工程扩展时加上不影响理解链路工作方式。接收端最关键的坑就是FFT的缩放因子要和发射端对应发射乘了sqrt(N)接收就要除以sqrt(N)否则星座点大小全是错的。4. 实测结果、调参技巧与问题排查仿真做完肯定要分析结果。这节我把自己跑出来的结果和几个印象深刻的调试经历写出来你们对照着查能少吃点亏。4.1 一组可复现的仿真结果按上面的参数在AWGN信道下的QPSK仿真结果大概是下面这个水平SNR (dB)QPSK 仿真BER16QAM 仿真BER6约 1.8e-2约 1.5e-110约 5.6e-4约 5.3e-214约 8.9e-6约 1.2e-218约 1e-7以下约 9.7e-4QPSK和16QAM之间大概差了8到9dB的功率代价这不奇怪高阶调制星座点距离更近抗噪声能力弱。如果你仿真结果里QPSK在6dB时BER是百分之五以上那多半是共轭对称、功率归一化或者信道估计哪一个环节出了问题不是参数选得不对。另外我记得第一次跑带多径信道版本时SNR12dB下BER一直卡在1e-3下不去。后来把循环前缀从16改成32BER立马降到1e-5以下。这是因为多径时延扩展大约占了6到7个采样点循环前缀只有16时其实够用但我那把发射端滤波器也叠加进去了整体时延扩展变长CP不够就开始产生符号间干扰。这个案例说明CP长度不是越大越好太大浪费频谱效率太小抗不了多径。一般设最大时延扩展的2到4倍比较合适。4.2 我踩过的一些坑这些坑没有一个是写论文能看到的全是实操出来的血泪。第一个星座图旋转且无法收敛。当时我16QAM的BER曲线在10dB左右开始走平降到1e-2就再也不动了。排查很久发现是接收端FFT缩放因子写错了发射端ifft乘了sqrt(N)接收端fft却忘了除结果星座点往外扩了一圈判决边界全乱。这个问题症状很隐蔽因为误码率不是完全不可用而是高到不合常理。第二个把直流偏置加在循环前缀之前还是之后影响了整帧平均光功率的计算。如果先加DC再做CP扩展那么直流分量会被重复复制到循环前缀里等效于直流分量变强但LED限幅是对整个时域信号统一处理的所以大多数人习惯的做法是加CP之后再偏置。为了让标准统一我建议把偏置和限幅都放在加完CP之后否则实际光功率和理论计算偏差很大。第三个统计误码率的窗口太小导致曲线抖动严重。早期我每个SNR点只送10个OFDM符号算下来的BER忽高忽低画出的曲线像心脏骤停的心电图。后来我把帧数提高到100个OFDM符号也就是单个SNR点有20万比特曲线才平滑下来。做仿真实验时尤其是要让曲线和理论值对比统计样本一定给足宁可跑得慢一点也不要拿抖动严重的图糊弄自己。第四个LED限幅参数和直流偏置的关系。如果你把LED_max设得足够大限幅几乎不生效那系统就退化成纯线性DCO-OFDM但你把LED_max设为比如1.2削波会非常严重不仅BER升高还会引入额外的非线性谐波干扰。想要模拟一个真实LED的功率受限特性需要先测波形的峰值因子PAPR再结合驱动摆幅限制去设置限幅门限不然结果没有工程参考价值。最后再说一个我自己的经验做完这套OFDM可见光通信系统仿真之后我最大的体会不是代码怎么写而是整个链路里每一步都环环相扣任何一个细节的参数不匹配都会让最终BER曲线变得不可信。从QAM符号归一化、共轭对称的子载波索引、循环前缀长度选择到接收端FFT缩放和信道估计随便一个地方出错排查起来都非常耗时。如果你也是零基础起步建议先按QPSK把端到端链路跑通确认BER趋势和理论曲线基本贴合再慢慢把16QAM、多径信道、LED非线性限幅一个个加进去。这样每一步都能定位到问题不至于最后曲线错了都不知道在哪一环出的错。仿真里把所有参数和中间变量打出来反复检查比跑一百遍曲线都管用。本文还有配套的精品资源点击获取