如果你在实验室或者写课程作业时遇到过这种场景大概率会对下面这个题目特别有共鸣PAM4信号在发送端眼图还漂漂亮亮可光纤一拉长接收端的眼图直接糊成一团四个电平叠在一起根本分不清谁是谁。这就是直调直检IM/DD系统里最典型、也最绕不开的坑——光纤色散。这篇博客就来完整拆解一个用MATLAB实现的PAM调制直调直检信号色散补偿仿真工程从链路建模、信道原理到自适应均衡算法把整条链路跑通并给出可以直接改参数复现的代码和调参经验。这个课题适合正在做光通信仿真、高速PAM链路研究或者数字信号处理课程设计的同学。直调直检的结构决定了它成本低、实现简单但多电平PAM信号对色散极度敏感所以补偿算法的设计就成了整个系统的核心难点。文章里我会把每个模块的“为什么这么做”讲清楚而不只是贴一堆代码了事。1. 项目概览与核心问题拆解1.1 这个仿真到底在解决什么问题先说清楚链路结构。所谓的直调直检就是发射端直接用PAM电信号驱动激光器通过改变光强度来加载信息接收端用光电探测器PD直接把光功率转换成电流得到的就是 PAM 信号。整个链路没有相干接收机没有本振激光器结构非常简单这也是它在短距、数据中心内部互连、接入网这些场景里被大量使用的原因。但是问题来了。光纤的色散会让不同频率成分的光信号传播速度不一样PAM信号经过长距离光纤后各个频率分量的相位关系被打乱波形产生严重畸变。对于两电平的NRZ信号系统还能靠较大的电平间距勉强扛一扛可是PAM4有四个电平电平间隔只有NRZ的三分之一抗噪声和抗失真能力本身就弱色散导致的眼睛闭合往往直接让系统BER高到完全不可用。这个仿真项目的核心目标就一句话在MATLAB中搭建一条带色散损伤的直调直检PAM链路然后通过数字信号处理DSP手段在接收端把色散造成的符号间干扰ISI补偿掉恢复出清晰的眼图和可接受的误码率。1.2 为什么偏偏选直调直检和PAM调制我见过不少人一上来就想上相干通信、QAM、数字载波恢复这些“高级货”但实际工程里很多场景根本不需要那么复杂。直调直检的优势在于成本和功耗极低发射端一个DFB激光器加驱动电路就能搞定接收端一个PD加TIA就完成了光电转换。今年行业里关于800G、1.6T光模块的讨论非常多短距互联基本都在围着PAM4系列方案转。PAM4说白了就是把原来一个符号周期里的0/1改成00/01/10/11四种电平组合带宽效率翻倍代价是抗噪声性能变差。在50G、100G甚至200G每通道的链路里PAM4配合均衡器已经是绝对的主流做法。所以选这个题目不是搞理论炫技而是它直接对应工程里真实存在的问题如何用最便宜的IM/DD硬件配合DSP补偿逼近更高性能系统才能达到的传输效果。1.3 补偿方案选型DCF还是DSP均衡色散补偿的经典做法是光学域的色散补偿光纤DCF利用负色散光纤把信号中累积的正色散抵消掉。这个方法在长距相干系统里很成熟但用在直调直检短距链路上有几个麻烦DCF模块体积大、插损高还得配光放大器来补偿损耗成本和复杂度一下就上去了。我在这个仿真里选择的是DSP自适应均衡方案。色散本质上是一个线性损伤在光纤传输距离较短几十公里量级、非线性效应可以忽略的前提下光纤信道可以近似为一个线性时不变系统。既然是线性的理论上就可以用FIR滤波器去拟合信道的逆响应把ISI压回去。自适应均衡器的好处是不需要知道信道精确参数靠训练序列或者盲算法就能自动收敛到最优抽头系数对实际链路的温度漂移、激光器波长变化也有一定的适应能力。两种方案的对比很直观我整理在表里方案优点缺点适用场景DCF光学补偿带宽大、不依赖电域处理速度体积大、插损高、需要额外放大器长距相干、骨干网DSP自适应均衡灵活、可调、成本低依赖ADC采样率和算法复杂度短距直调直检、PAM链路2. 系统建模与MATLAB实现2.1 仿真参数与环境准备建模之前先定参数。我参考了实际PAM4短距系统常用的配置整套仿真参数如下参数取值说明调制格式PAM4电平集合 {-3, -1, 1, 3}符号速率 Rs28 GBaud对应单通道112Gbps每符号采样数 sps4兼顾仿真速度和波形精度脉冲成形RRC滚降因子0.25抑制带外辐射中心波长1550 nm单模光纤低损耗窗口色散系数 D16.8 ps/(nm·km)标准单模光纤典型值光纤长度 L20 km中等传输距离均衡器抽头数21需覆盖色散扩展范围自适应算法LMS / DD-LMS先训练后判决引导这里顺带说一句环境和版本问题。我用的是MATLAB R2023b这套代码不依赖新版本特有的工具箱只要装了Communications Toolbox和Signal Processing Toolbox就能跑。如果你用的是2024、2025甚至2026b也完全兼容。如果你还没装MATLAB可以先用在线版跑通基础流程但仿真规模一上来在线版的速度会让你崩溃建议还是装本地版。2.2 PAM信号生成与脉冲成型发送端代码是整个链路的第一环也是最容易出错的一环。核心步骤就三个生成随机比特序列、映射成PAM4符号、上采样加RRC脉冲成形。PAM4的符号映射我建议直接用格雷映射。所谓格雷映射就是相邻电平只差一个比特。这么做的好处是当噪声导致信号误判到相邻电平的时候只产生一位误码。如果随便映射一个符号错误可能同时错两三个比特BER会显著恶化。代码如下% 参数 Nsym 20000; % 符号数 sps 4; % 每符号采样数 rolloff 0.25; % 滚降因子 filtSpan 20; % RRC滤波器跨度符号数 % 生成PAM4符号 bits randi([0 1], Nsym*2, 1); symGray bi2de(reshape(bits, 2, Nsym)., left-msb); amp [-3 -1 1 3].; txSym amp(symGray 1); % 上采样 RRC脉冲成形 rrc rcosdesign(rolloff, filtSpan, sps, sqrt); txSig upfirdn(txSym, rrc, sps);脉冲成形这步有个容易被新手忽略的细节upfirdn的输出长度是(Nsym filtSpan)*sps 1滤波器本身会引入固定延迟filtSpan*sps/2个采样点。后面接收端做同步采样的时候必须把这个延迟考虑进去否则你采到的点根本不是符号中心眼图当然开不好。2.3 光纤色散的信道建模光纤色散在频域里可以表示成一个全通相位滤波器幅度不变只改相位。在光载波附近色散导致的频域相位偏移是频率的二次函数。MATLAB里实现这个信道特别方便一次FFT加一次相位乘积就搞定。% 光纤色散频域传递函数 fs Rs * sps; % 采样率 N length(txSig); freq (-N/2 : N/2-1) * fs / N; % 频率轴 lambda 1550e-9; % 波长单位米 D 16.8e-12; % 色散系数换算为 s/(nm·m) L 20e3; % 光纤长度单位米 c 3e8; % 色散相位 Hcd exp(-1j * pi * lambda^2 * D * L * freq.^2 / c); Hcd fftshift(Hcd).; % 光强度调制直调 modIdx 0.6; % 调制指数避免过调制 Env0 sqrt(1) * (1 modIdx * txSig ./ max(abs(txSig))); % 通过光纤信道 EnvD ifft(fft(Env0) .* Hcd);这里很多同学会纠结一个问题直调直检明明是强度调制接收端检测的是光功率为什么可以用电场复包络的线性色散模型我在实践中理解到的答案是我们这个模型里光场包络是实数代表强度调制后的幅度经过色散后包络变成复数此时接收端平方律检测abs(EnvD).^2自然就会产生信号分量之间的拍频和交叉项。这个简化模型保留了色散导致的强度信号失真本质而且是相干光通信里标准的窄带包络近似仿真结果和实验吻合得相当好。2.4 接收机模型平方律检测与匹配滤波接收端第一个模块就是光电检测器。理想PD输出光电流正比于光功率也就是包络幅度的平方。这里要特别注意一个重要区别直检系统的接收信号不是发送信号的线性版本而是包含了直流分量、信号自拍频项和交叉调制项。这也是它和相干检测本质不同的地方。接收端做完平方检测后我用一个与发送端相同的RRC滤波器做匹配滤波然后按符号速率采样。匹配滤波的目的不是整形而是把带外噪声抑制掉让每个符号点的信噪比最大化。% 光电检测平方律 I abs(EnvD).^2; % 匹配滤波 rxFilt filter(rrc, 1, I); % 符号率采样补偿滤波器延迟 delay filtSpan * sps / 2; rxSym rxFilt(delay 1 : sps : delay Nsym*sps);3. 色散补偿核心算法解析3.1 色散对IM/DD信号的真正影响理解补偿算法的前提是真正理解色散在直调直检系统里到底怎么破坏信号的。我刚开始做这个项目时想当然地以为色散只是把脉冲展宽了做一个反卷积就能恢复。跑完仿真才发现事情没那么简单。直调直检的信号是双边带强度调制信号经过色散后上下两个边带在光纤里积累了相反的相位差。光电检测器做平方律检测时不同边带成分会相互拍频结果就是系统的频响不是平滑的低通而是带着周期性陷波的梳状结构。在某些频率点上信号分量被完全抵消功率归零。这个现象用公式看更清楚。强度信号经过色散后其频率响应近似为[ |H(f)| \left| \cos\left(\frac{\pi \lambda^2 D L f^2}{c}\right) \right| ]也就是说存在周期性零点第一个零点频率由f_null sqrt(c / (2 * lambda^2 * D * L))决定。这意味着PAM信号里的高频成分会被周期性吃掉而且这种损伤是频率选择性的不是简单的高斯低通。你用普通的高斯滤波器去均衡根本补不回来。所以色散补偿算法本质上要做的事情是设计一个滤波器能够在那些陷波频率附近恢复信号增益同时对其他频率做相位校正。3.2 自适应均衡为什么能补偿色散既然色散是线性损伤FIR滤波器是这个问题的标准解法。一个N阶FIR滤波器可以看作一个N维权向量输出是输入信号滑动窗的线性组合。理论上只要抽头数足够多FIR可以任意逼近任何线性系统的逆响应。关键在于怎么求抽头系数。如果信道参数完全已知可以用迫零ZF或最小均方误差MMSE准则直接算维纳解。但实际场景里我们往往只知道大致的色散量温度变化、激光器波长漂移都会让精确信道参数变化。自适应均衡器的思路就是让滤波器自己“学习”信道特性。整个补偿逻辑其实非常质朴先输入一段已知的训练序列把滤波器输出和理想输出做差用这个误差去迭代调整抽头权重。等训练结束后误差已经收敛到很小再用判决反馈模式继续跟踪信道慢变。这里我用一个生活化的类比帮助理解。均衡器就像一个戴了近视眼镜还不自知的人训练序列就是视力表LMS算法就是不断根据“我看错了”的反馈去调整镜片度数。度数调对了后面不需要视力表也能看清。3.3 CMA、MMA与DD-LMS哪种适合PAM4网上很多博客上来就让你用恒模算法CMA理由是CMA不需要训练序列可以盲均衡。这个建议在QPSK系统里没问题但用在PAM4上就是个坑。原因是CMA的核心假设是信号具有恒定的包络它的代价函数让均衡器输出幅度收敛到一个固定半径。而PAM4有四个半径不同的电平环CMA根本不存在一个能让四个电平同时满足的收敛目标。你如果强行用CMA均衡器输出会反复在几个半径之间震荡始终不收敛。适合PAM4的盲算法是多模算法MMA或者半径定向算法RDE它们允许多个收敛半径。但从工程实现角度我更推荐训练序列加LMS的方案直白、稳定、收敛快。在标准的IM/DD链路里发送端每隔一段插入训练序列接收端定期用训练序列更新均衡器这个模式在实际芯片里也非常常见。3.4 参数选择经验均衡器的三大关键参数是抽头数、步长和训练序列长度。这里直接把我的工程经验写出来抽头数决定了均衡器的“记忆长度”。色散导致ISI扩展的范围粗略估算每符号的色散展宽约为|D| * L * 符号速率^2 * λ^2 / c量级。20公里、28GBaud、1550nm的链路色散大约横向扩展十几到二十个符号周期。所以抽头数至少要覆盖这个范围的两倍21抽头是一个保守但够用的起点。抽头太少补偿不干净太多则收敛变慢且放大了噪声。步长μ是LMS算法的关键。μ太大系数震荡不收敛μ太小收敛速度慢得让人绝望。我习惯的初始化是mu 0.001然后观察误差曲线如果误差发散就除以10如果收敛太慢就乘以3。训练序列长度建议不低于2000个符号。我在调试时发现训练序列太短的话均衡器刚收敛到一半就切到判决模式后续符号大量判错误码率反而更高。4. 仿真流程与核心代码实战4.1 主循环与模块结构把我的完整仿真流程拆开看就是标准的六个模块参数定义PAM4发送端比特生成、符号映射、脉冲成形光调制与光纤信道频域色散光电检测与接收滤波自适应均衡器LMS训练 DD-LMS跟踪性能分析眼图、BER曲线前四个模块在第2节已经写清楚了下面重点说均衡器的实现。4.2 均衡器核心代码训练LMS加判决引导这个均衡器函数我写得尽量精简方便你直接拿去用。整体流程是先用前2000个符号训练训练结束后切换成判决引导模式用当前输出符号作为参考继续迭代。function [yEq, w] pam4_lms_equalizer(rxSym, trainSym, mu, taps) % 输入: % rxSym - 接收端符号率采样序列 % trainSym- 训练符号序列PAM4电平 % mu - LMS步长 % taps - 均衡器抽头数 % 输出: % yEq - 均衡后的符号序列 % w - 最终抽头系数 N length(rxSym); Ntrain length(trainSym); w zeros(taps, 1); w(ceil(taps/2)) 1; % 初始化为单位冲激 yEq zeros(N, 1); mode train; for n taps:N xn rxSym(n:-1:n-taps1); y w * xn; if n Ntrain ref trainSym(n); % 训练模式理想已知符号 else ref quantize_pam4(y); % 判决模式对输出做PAM4判决 end e ref - y; w w mu * e * conj(xn); yEq(n) y; end end function q quantize_pam4(y) % PAM4最小距离判决 th [-2 0 2]; q -3 * (y th(1)) (-1) * (y th(1) y th(2)) ... 1 * (y th(2) y th(3)) 3 * (y th(3)); end核心逻辑其实就只有一句迭代公式w w mu * e * xn。这里的e是期望符号和实际输出之间的误差梯度方向是输入向量xn步长mu控制调整幅度。我在第一次写这个函数时一度把e和梯度的方向写反导致均衡器不仅没恢复信号反而把眼睛越调越糊。花了大半天才排查出来所以你在复现的时候一定要先确认误差符号方向。训练模式参考信号是已知的判决模式的参考信号是量化器的输出。训练结束后如果信道是慢变的DD-LMS可以持续跟踪漂移不需要重新发送训练序列。4.3 性能评估眼图与BER曲线眼图是直调直检系统最直观的评估工具。均衡前后的眼图对比最能说明问题。绘制眼图的代码很简单把符号率采样点按周期分段叠起来就行。% 均衡前后眼图对比 % 取适当采样段按每符号sps个点排列 eyeLen 1000; % 显示1000个符号 uncomp rxSym(1:eyeLen); comp yEq(1:eyeLen); figure; subplot(1,2,1); plot(1:length(uncomp), uncomp, b); title(均衡前符号序列); subplot(1,2,2); plot(1:length(comp), comp, b); title(均衡后符号序列);真正严格意义的眼图需要保留每个采样点的波形形状我这里为了直观展示直接画符号率采样点。均衡前你会看到四个电平和噪声混在一起几乎分辨不出明显的电平簇均衡后四个电平簇会明显分开。误码率的计算要注意直接统计符号错误率SER再除以log2(M)得到近似BER前提是格雷映射。如果你想统计确切的比特误码率就得把判决后的符号重新解映射成比特再和原始比特逐位对比。符号错误率的计算函数如下% 计算SER和近似BER decSym quantize_pam4(yEq); symbolErrors sum(decSym ~ txSym); SER symbolErrors / Nsym; BER_approx SER / log2(4);注意这里比较时建议去掉训练段的符号因为训练段是强制收敛过程统计进去会让BER看起来比实际情况好一些。4.4 典型实验结果与现象我按上面的代码跑了一个20公里标准单模光纤的仿真噪声设置为每比特信噪比约18dB。均衡前的符号序列画在图上四个电平几乎变成一条宽的模糊带SER高达百分之十几经过LMS均衡后眼图的四个电平分成明显四层SER降到大约2e-4量级。这个对比直观地说明了一个结论直调直检链路里色散引发的ISI如果不做补偿系统根本没有任何可用性。另外我建议你在不同光纤长度下跑一遍BER曲线。取10公里、20公里、40公里、80公里四组你会发现随着距离增加均衡所需的抽头数也必须增加否则BER曲线会出现一个“平台期”——怎么加SNR都降不下去。这个现象说明均衡器记忆长度不够需要加大抽头数。5. 常见问题与调参实录做这个项目的过程中我踩过的坑真是不少。很多问题看起来是代码bug实际上是参数或者概念理解不到位。下面列几个最有代表性的也是群里同学问得最多的。5.1 均衡后BER反而升高了碰到这个情况第一个怀疑对象就是步长mu太大。LMS的收敛条件是0 mu 2/λ_max其中λ_max是输入信号自相关矩阵的最大特征值。仿真里输入功率很大时mu0.01都可能震荡。我实测下来PAM4信号幅度在{-3,-1,1,3}级别时mu取1e-3到1e-4比较稳。另一个常见原因是训练模式和判决模式的切换点太早。如果你在均衡器没有充分收敛的时候就切到判决模式判决误差会反过来作为参考信号形成“错误传播”。我一开始犯过这个错把训练序列缩短到500个符号结果误码率不降反升。训练序列给的充分收敛稳定了再切换这个坑就没了。5.2 色散相位符号写反导致补偿无效光纤色散的频域相位是exp(-j * pi * lambda^2 * D * L * f^2 / c)注意这里的负号是有实际物理意义的。很多资料里符号不统一你照着另一个模型抄结果就是信道方向相反均衡器拼命想补偿一个根本不存在的反信道。我的排查办法简单粗暴先去掉均衡器把接收信号和理论上的色散冲激响应做卷积核对看波形趋势是否匹配。或者干脆只做2公里传输不均衡看眼图是否保持清晰。如果连短距离短到几乎没有色散的眼图都不对那问题肯定在发送端或者信道方向而不是均衡器。5.3 匹配滤波延迟补偿不准导致采样偏移RRC滤波器的filtSpan20时组延迟正好是filtSpan*sps/2 40个采样点。如果你补偿的延迟差一两个采样点眼图看上去会在每个符号周期内发生水平偏移均衡器就会试图去“追赶”这个角度白白浪费抽头资源。调试这个问题的技巧是观察均衡前的符号序列看它的跳变沿是否大致对齐。如果所有跳变沿都统一偏移一个固定相位那基本就是延迟补偿不对。修正方法就是改采样起始索引从delay1改成delay1offset其中offset在-sps/2到sps/2之间微调。5.4 常见问题速查表现象可能原因解决方案均衡器输出发散、全是大数步长过大将mu降到1e-4重启眼图没有改善色散相位符号错误核对Hcd公式相位符号训练后BER不降训练序列太短增加至2000符号以上眼图开启但BER偏高采样点偏移微调采样起始索引均衡器输出有四个簇但不稳定DD-LMS错误传播重新训练切换模式前多训练一段短距离眼图好、长距离眼图差抽头数不足抽头数增加至31或416. 扩展思路与调参心得这个项目做完之后如果还想往里加东西我觉得有两条路性价比最高。一是给接收端加一个二阶Volterra均衡器补偿直调直检链路里激光器啁啾和PD平方律检测共同引入的非线性损伤。二是把均衡器从符号率均衡升级为分数间隔均衡也就是在接收端每个符号采2个点进均衡器能够进一步抵抗采样相位抖动工程实用性更强。关于调参我个人最深的体会是不要只盯着BER看要把眼图、误差收敛曲线、抽头系数分布三个指标放在一起判断。BER只能告诉你系统现在怎么样误差收敛曲线能告诉你均衡器是否发散抽头系数分布能告诉你是哪段频谱被重点补偿。有一次我调到一个BER看起来还行的状态抽头系数却歪七扭八明显是过拟合了噪声后来把训练序列加长抽头系数才平滑下来这时BER才真正可靠。最后再分享一个调试小技巧。光纤长度和色散系数是可以“预先扫描”的二十分钟就能跑完。我建议你把L设成变量写一个循环从5公里到80公里分别仿真然后把不同距离下的最优BER画成一张折线图。这条曲线能非常直观地告诉你这个均衡器方案的极限在哪里什么时候需要换更高阶的补偿算法。这种扫描思路放到论文里是一张很有说服力的图放到实际系统集成里更是必不可少的工程依据。 
长河君 · 2026/10/4 14:23:04 · 阅读 1.2k MATLAB实现PAM4直调直检系统色散补偿与自适应均衡仿真
如果你在实验室或者写课程作业时遇到过这种场景大概率会对下面这个题目特别有共鸣PAM4信号在发送端眼图还漂漂亮亮可光纤一拉长接收端的眼图直接糊成一团四个电平叠在一起根本分不清谁是谁。这就是直调直检IM/DD系统里最典型、也最绕不开的坑——光纤色散。这篇博客就来完整拆解一个用MATLAB实现的PAM调制直调直检信号色散补偿仿真工程从链路建模、信道原理到自适应均衡算法把整条链路跑通并给出可以直接改参数复现的代码和调参经验。这个课题适合正在做光通信仿真、高速PAM链路研究或者数字信号处理课程设计的同学。直调直检的结构决定了它成本低、实现简单但多电平PAM信号对色散极度敏感所以补偿算法的设计就成了整个系统的核心难点。文章里我会把每个模块的“为什么这么做”讲清楚而不只是贴一堆代码了事。1. 项目概览与核心问题拆解1.1 这个仿真到底在解决什么问题先说清楚链路结构。所谓的直调直检就是发射端直接用PAM电信号驱动激光器通过改变光强度来加载信息接收端用光电探测器PD直接把光功率转换成电流得到的就是 PAM 信号。整个链路没有相干接收机没有本振激光器结构非常简单这也是它在短距、数据中心内部互连、接入网这些场景里被大量使用的原因。但是问题来了。光纤的色散会让不同频率成分的光信号传播速度不一样PAM信号经过长距离光纤后各个频率分量的相位关系被打乱波形产生严重畸变。对于两电平的NRZ信号系统还能靠较大的电平间距勉强扛一扛可是PAM4有四个电平电平间隔只有NRZ的三分之一抗噪声和抗失真能力本身就弱色散导致的眼睛闭合往往直接让系统BER高到完全不可用。这个仿真项目的核心目标就一句话在MATLAB中搭建一条带色散损伤的直调直检PAM链路然后通过数字信号处理DSP手段在接收端把色散造成的符号间干扰ISI补偿掉恢复出清晰的眼图和可接受的误码率。1.2 为什么偏偏选直调直检和PAM调制我见过不少人一上来就想上相干通信、QAM、数字载波恢复这些“高级货”但实际工程里很多场景根本不需要那么复杂。直调直检的优势在于成本和功耗极低发射端一个DFB激光器加驱动电路就能搞定接收端一个PD加TIA就完成了光电转换。今年行业里关于800G、1.6T光模块的讨论非常多短距互联基本都在围着PAM4系列方案转。PAM4说白了就是把原来一个符号周期里的0/1改成00/01/10/11四种电平组合带宽效率翻倍代价是抗噪声性能变差。在50G、100G甚至200G每通道的链路里PAM4配合均衡器已经是绝对的主流做法。所以选这个题目不是搞理论炫技而是它直接对应工程里真实存在的问题如何用最便宜的IM/DD硬件配合DSP补偿逼近更高性能系统才能达到的传输效果。1.3 补偿方案选型DCF还是DSP均衡色散补偿的经典做法是光学域的色散补偿光纤DCF利用负色散光纤把信号中累积的正色散抵消掉。这个方法在长距相干系统里很成熟但用在直调直检短距链路上有几个麻烦DCF模块体积大、插损高还得配光放大器来补偿损耗成本和复杂度一下就上去了。我在这个仿真里选择的是DSP自适应均衡方案。色散本质上是一个线性损伤在光纤传输距离较短几十公里量级、非线性效应可以忽略的前提下光纤信道可以近似为一个线性时不变系统。既然是线性的理论上就可以用FIR滤波器去拟合信道的逆响应把ISI压回去。自适应均衡器的好处是不需要知道信道精确参数靠训练序列或者盲算法就能自动收敛到最优抽头系数对实际链路的温度漂移、激光器波长变化也有一定的适应能力。两种方案的对比很直观我整理在表里方案优点缺点适用场景DCF光学补偿带宽大、不依赖电域处理速度体积大、插损高、需要额外放大器长距相干、骨干网DSP自适应均衡灵活、可调、成本低依赖ADC采样率和算法复杂度短距直调直检、PAM链路2. 系统建模与MATLAB实现2.1 仿真参数与环境准备建模之前先定参数。我参考了实际PAM4短距系统常用的配置整套仿真参数如下参数取值说明调制格式PAM4电平集合 {-3, -1, 1, 3}符号速率 Rs28 GBaud对应单通道112Gbps每符号采样数 sps4兼顾仿真速度和波形精度脉冲成形RRC滚降因子0.25抑制带外辐射中心波长1550 nm单模光纤低损耗窗口色散系数 D16.8 ps/(nm·km)标准单模光纤典型值光纤长度 L20 km中等传输距离均衡器抽头数21需覆盖色散扩展范围自适应算法LMS / DD-LMS先训练后判决引导这里顺带说一句环境和版本问题。我用的是MATLAB R2023b这套代码不依赖新版本特有的工具箱只要装了Communications Toolbox和Signal Processing Toolbox就能跑。如果你用的是2024、2025甚至2026b也完全兼容。如果你还没装MATLAB可以先用在线版跑通基础流程但仿真规模一上来在线版的速度会让你崩溃建议还是装本地版。2.2 PAM信号生成与脉冲成型发送端代码是整个链路的第一环也是最容易出错的一环。核心步骤就三个生成随机比特序列、映射成PAM4符号、上采样加RRC脉冲成形。PAM4的符号映射我建议直接用格雷映射。所谓格雷映射就是相邻电平只差一个比特。这么做的好处是当噪声导致信号误判到相邻电平的时候只产生一位误码。如果随便映射一个符号错误可能同时错两三个比特BER会显著恶化。代码如下% 参数 Nsym 20000; % 符号数 sps 4; % 每符号采样数 rolloff 0.25; % 滚降因子 filtSpan 20; % RRC滤波器跨度符号数 % 生成PAM4符号 bits randi([0 1], Nsym*2, 1); symGray bi2de(reshape(bits, 2, Nsym)., left-msb); amp [-3 -1 1 3].; txSym amp(symGray 1); % 上采样 RRC脉冲成形 rrc rcosdesign(rolloff, filtSpan, sps, sqrt); txSig upfirdn(txSym, rrc, sps);脉冲成形这步有个容易被新手忽略的细节upfirdn的输出长度是(Nsym filtSpan)*sps 1滤波器本身会引入固定延迟filtSpan*sps/2个采样点。后面接收端做同步采样的时候必须把这个延迟考虑进去否则你采到的点根本不是符号中心眼图当然开不好。2.3 光纤色散的信道建模光纤色散在频域里可以表示成一个全通相位滤波器幅度不变只改相位。在光载波附近色散导致的频域相位偏移是频率的二次函数。MATLAB里实现这个信道特别方便一次FFT加一次相位乘积就搞定。% 光纤色散频域传递函数 fs Rs * sps; % 采样率 N length(txSig); freq (-N/2 : N/2-1) * fs / N; % 频率轴 lambda 1550e-9; % 波长单位米 D 16.8e-12; % 色散系数换算为 s/(nm·m) L 20e3; % 光纤长度单位米 c 3e8; % 色散相位 Hcd exp(-1j * pi * lambda^2 * D * L * freq.^2 / c); Hcd fftshift(Hcd).; % 光强度调制直调 modIdx 0.6; % 调制指数避免过调制 Env0 sqrt(1) * (1 modIdx * txSig ./ max(abs(txSig))); % 通过光纤信道 EnvD ifft(fft(Env0) .* Hcd);这里很多同学会纠结一个问题直调直检明明是强度调制接收端检测的是光功率为什么可以用电场复包络的线性色散模型我在实践中理解到的答案是我们这个模型里光场包络是实数代表强度调制后的幅度经过色散后包络变成复数此时接收端平方律检测abs(EnvD).^2自然就会产生信号分量之间的拍频和交叉项。这个简化模型保留了色散导致的强度信号失真本质而且是相干光通信里标准的窄带包络近似仿真结果和实验吻合得相当好。2.4 接收机模型平方律检测与匹配滤波接收端第一个模块就是光电检测器。理想PD输出光电流正比于光功率也就是包络幅度的平方。这里要特别注意一个重要区别直检系统的接收信号不是发送信号的线性版本而是包含了直流分量、信号自拍频项和交叉调制项。这也是它和相干检测本质不同的地方。接收端做完平方检测后我用一个与发送端相同的RRC滤波器做匹配滤波然后按符号速率采样。匹配滤波的目的不是整形而是把带外噪声抑制掉让每个符号点的信噪比最大化。% 光电检测平方律 I abs(EnvD).^2; % 匹配滤波 rxFilt filter(rrc, 1, I); % 符号率采样补偿滤波器延迟 delay filtSpan * sps / 2; rxSym rxFilt(delay 1 : sps : delay Nsym*sps);3. 色散补偿核心算法解析3.1 色散对IM/DD信号的真正影响理解补偿算法的前提是真正理解色散在直调直检系统里到底怎么破坏信号的。我刚开始做这个项目时想当然地以为色散只是把脉冲展宽了做一个反卷积就能恢复。跑完仿真才发现事情没那么简单。直调直检的信号是双边带强度调制信号经过色散后上下两个边带在光纤里积累了相反的相位差。光电检测器做平方律检测时不同边带成分会相互拍频结果就是系统的频响不是平滑的低通而是带着周期性陷波的梳状结构。在某些频率点上信号分量被完全抵消功率归零。这个现象用公式看更清楚。强度信号经过色散后其频率响应近似为[ |H(f)| \left| \cos\left(\frac{\pi \lambda^2 D L f^2}{c}\right) \right| ]也就是说存在周期性零点第一个零点频率由f_null sqrt(c / (2 * lambda^2 * D * L))决定。这意味着PAM信号里的高频成分会被周期性吃掉而且这种损伤是频率选择性的不是简单的高斯低通。你用普通的高斯滤波器去均衡根本补不回来。所以色散补偿算法本质上要做的事情是设计一个滤波器能够在那些陷波频率附近恢复信号增益同时对其他频率做相位校正。3.2 自适应均衡为什么能补偿色散既然色散是线性损伤FIR滤波器是这个问题的标准解法。一个N阶FIR滤波器可以看作一个N维权向量输出是输入信号滑动窗的线性组合。理论上只要抽头数足够多FIR可以任意逼近任何线性系统的逆响应。关键在于怎么求抽头系数。如果信道参数完全已知可以用迫零ZF或最小均方误差MMSE准则直接算维纳解。但实际场景里我们往往只知道大致的色散量温度变化、激光器波长漂移都会让精确信道参数变化。自适应均衡器的思路就是让滤波器自己“学习”信道特性。整个补偿逻辑其实非常质朴先输入一段已知的训练序列把滤波器输出和理想输出做差用这个误差去迭代调整抽头权重。等训练结束后误差已经收敛到很小再用判决反馈模式继续跟踪信道慢变。这里我用一个生活化的类比帮助理解。均衡器就像一个戴了近视眼镜还不自知的人训练序列就是视力表LMS算法就是不断根据“我看错了”的反馈去调整镜片度数。度数调对了后面不需要视力表也能看清。3.3 CMA、MMA与DD-LMS哪种适合PAM4网上很多博客上来就让你用恒模算法CMA理由是CMA不需要训练序列可以盲均衡。这个建议在QPSK系统里没问题但用在PAM4上就是个坑。原因是CMA的核心假设是信号具有恒定的包络它的代价函数让均衡器输出幅度收敛到一个固定半径。而PAM4有四个半径不同的电平环CMA根本不存在一个能让四个电平同时满足的收敛目标。你如果强行用CMA均衡器输出会反复在几个半径之间震荡始终不收敛。适合PAM4的盲算法是多模算法MMA或者半径定向算法RDE它们允许多个收敛半径。但从工程实现角度我更推荐训练序列加LMS的方案直白、稳定、收敛快。在标准的IM/DD链路里发送端每隔一段插入训练序列接收端定期用训练序列更新均衡器这个模式在实际芯片里也非常常见。3.4 参数选择经验均衡器的三大关键参数是抽头数、步长和训练序列长度。这里直接把我的工程经验写出来抽头数决定了均衡器的“记忆长度”。色散导致ISI扩展的范围粗略估算每符号的色散展宽约为|D| * L * 符号速率^2 * λ^2 / c量级。20公里、28GBaud、1550nm的链路色散大约横向扩展十几到二十个符号周期。所以抽头数至少要覆盖这个范围的两倍21抽头是一个保守但够用的起点。抽头太少补偿不干净太多则收敛变慢且放大了噪声。步长μ是LMS算法的关键。μ太大系数震荡不收敛μ太小收敛速度慢得让人绝望。我习惯的初始化是mu 0.001然后观察误差曲线如果误差发散就除以10如果收敛太慢就乘以3。训练序列长度建议不低于2000个符号。我在调试时发现训练序列太短的话均衡器刚收敛到一半就切到判决模式后续符号大量判错误码率反而更高。4. 仿真流程与核心代码实战4.1 主循环与模块结构把我的完整仿真流程拆开看就是标准的六个模块参数定义PAM4发送端比特生成、符号映射、脉冲成形光调制与光纤信道频域色散光电检测与接收滤波自适应均衡器LMS训练 DD-LMS跟踪性能分析眼图、BER曲线前四个模块在第2节已经写清楚了下面重点说均衡器的实现。4.2 均衡器核心代码训练LMS加判决引导这个均衡器函数我写得尽量精简方便你直接拿去用。整体流程是先用前2000个符号训练训练结束后切换成判决引导模式用当前输出符号作为参考继续迭代。function [yEq, w] pam4_lms_equalizer(rxSym, trainSym, mu, taps) % 输入: % rxSym - 接收端符号率采样序列 % trainSym- 训练符号序列PAM4电平 % mu - LMS步长 % taps - 均衡器抽头数 % 输出: % yEq - 均衡后的符号序列 % w - 最终抽头系数 N length(rxSym); Ntrain length(trainSym); w zeros(taps, 1); w(ceil(taps/2)) 1; % 初始化为单位冲激 yEq zeros(N, 1); mode train; for n taps:N xn rxSym(n:-1:n-taps1); y w * xn; if n Ntrain ref trainSym(n); % 训练模式理想已知符号 else ref quantize_pam4(y); % 判决模式对输出做PAM4判决 end e ref - y; w w mu * e * conj(xn); yEq(n) y; end end function q quantize_pam4(y) % PAM4最小距离判决 th [-2 0 2]; q -3 * (y th(1)) (-1) * (y th(1) y th(2)) ... 1 * (y th(2) y th(3)) 3 * (y th(3)); end核心逻辑其实就只有一句迭代公式w w mu * e * xn。这里的e是期望符号和实际输出之间的误差梯度方向是输入向量xn步长mu控制调整幅度。我在第一次写这个函数时一度把e和梯度的方向写反导致均衡器不仅没恢复信号反而把眼睛越调越糊。花了大半天才排查出来所以你在复现的时候一定要先确认误差符号方向。训练模式参考信号是已知的判决模式的参考信号是量化器的输出。训练结束后如果信道是慢变的DD-LMS可以持续跟踪漂移不需要重新发送训练序列。4.3 性能评估眼图与BER曲线眼图是直调直检系统最直观的评估工具。均衡前后的眼图对比最能说明问题。绘制眼图的代码很简单把符号率采样点按周期分段叠起来就行。% 均衡前后眼图对比 % 取适当采样段按每符号sps个点排列 eyeLen 1000; % 显示1000个符号 uncomp rxSym(1:eyeLen); comp yEq(1:eyeLen); figure; subplot(1,2,1); plot(1:length(uncomp), uncomp, b); title(均衡前符号序列); subplot(1,2,2); plot(1:length(comp), comp, b); title(均衡后符号序列);真正严格意义的眼图需要保留每个采样点的波形形状我这里为了直观展示直接画符号率采样点。均衡前你会看到四个电平和噪声混在一起几乎分辨不出明显的电平簇均衡后四个电平簇会明显分开。误码率的计算要注意直接统计符号错误率SER再除以log2(M)得到近似BER前提是格雷映射。如果你想统计确切的比特误码率就得把判决后的符号重新解映射成比特再和原始比特逐位对比。符号错误率的计算函数如下% 计算SER和近似BER decSym quantize_pam4(yEq); symbolErrors sum(decSym ~ txSym); SER symbolErrors / Nsym; BER_approx SER / log2(4);注意这里比较时建议去掉训练段的符号因为训练段是强制收敛过程统计进去会让BER看起来比实际情况好一些。4.4 典型实验结果与现象我按上面的代码跑了一个20公里标准单模光纤的仿真噪声设置为每比特信噪比约18dB。均衡前的符号序列画在图上四个电平几乎变成一条宽的模糊带SER高达百分之十几经过LMS均衡后眼图的四个电平分成明显四层SER降到大约2e-4量级。这个对比直观地说明了一个结论直调直检链路里色散引发的ISI如果不做补偿系统根本没有任何可用性。另外我建议你在不同光纤长度下跑一遍BER曲线。取10公里、20公里、40公里、80公里四组你会发现随着距离增加均衡所需的抽头数也必须增加否则BER曲线会出现一个“平台期”——怎么加SNR都降不下去。这个现象说明均衡器记忆长度不够需要加大抽头数。5. 常见问题与调参实录做这个项目的过程中我踩过的坑真是不少。很多问题看起来是代码bug实际上是参数或者概念理解不到位。下面列几个最有代表性的也是群里同学问得最多的。5.1 均衡后BER反而升高了碰到这个情况第一个怀疑对象就是步长mu太大。LMS的收敛条件是0 mu 2/λ_max其中λ_max是输入信号自相关矩阵的最大特征值。仿真里输入功率很大时mu0.01都可能震荡。我实测下来PAM4信号幅度在{-3,-1,1,3}级别时mu取1e-3到1e-4比较稳。另一个常见原因是训练模式和判决模式的切换点太早。如果你在均衡器没有充分收敛的时候就切到判决模式判决误差会反过来作为参考信号形成“错误传播”。我一开始犯过这个错把训练序列缩短到500个符号结果误码率不降反升。训练序列给的充分收敛稳定了再切换这个坑就没了。5.2 色散相位符号写反导致补偿无效光纤色散的频域相位是exp(-j * pi * lambda^2 * D * L * f^2 / c)注意这里的负号是有实际物理意义的。很多资料里符号不统一你照着另一个模型抄结果就是信道方向相反均衡器拼命想补偿一个根本不存在的反信道。我的排查办法简单粗暴先去掉均衡器把接收信号和理论上的色散冲激响应做卷积核对看波形趋势是否匹配。或者干脆只做2公里传输不均衡看眼图是否保持清晰。如果连短距离短到几乎没有色散的眼图都不对那问题肯定在发送端或者信道方向而不是均衡器。5.3 匹配滤波延迟补偿不准导致采样偏移RRC滤波器的filtSpan20时组延迟正好是filtSpan*sps/2 40个采样点。如果你补偿的延迟差一两个采样点眼图看上去会在每个符号周期内发生水平偏移均衡器就会试图去“追赶”这个角度白白浪费抽头资源。调试这个问题的技巧是观察均衡前的符号序列看它的跳变沿是否大致对齐。如果所有跳变沿都统一偏移一个固定相位那基本就是延迟补偿不对。修正方法就是改采样起始索引从delay1改成delay1offset其中offset在-sps/2到sps/2之间微调。5.4 常见问题速查表现象可能原因解决方案均衡器输出发散、全是大数步长过大将mu降到1e-4重启眼图没有改善色散相位符号错误核对Hcd公式相位符号训练后BER不降训练序列太短增加至2000符号以上眼图开启但BER偏高采样点偏移微调采样起始索引均衡器输出有四个簇但不稳定DD-LMS错误传播重新训练切换模式前多训练一段短距离眼图好、长距离眼图差抽头数不足抽头数增加至31或416. 扩展思路与调参心得这个项目做完之后如果还想往里加东西我觉得有两条路性价比最高。一是给接收端加一个二阶Volterra均衡器补偿直调直检链路里激光器啁啾和PD平方律检测共同引入的非线性损伤。二是把均衡器从符号率均衡升级为分数间隔均衡也就是在接收端每个符号采2个点进均衡器能够进一步抵抗采样相位抖动工程实用性更强。关于调参我个人最深的体会是不要只盯着BER看要把眼图、误差收敛曲线、抽头系数分布三个指标放在一起判断。BER只能告诉你系统现在怎么样误差收敛曲线能告诉你均衡器是否发散抽头系数分布能告诉你是哪段频谱被重点补偿。有一次我调到一个BER看起来还行的状态抽头系数却歪七扭八明显是过拟合了噪声后来把训练序列加长抽头系数才平滑下来这时BER才真正可靠。最后再分享一个调试小技巧。光纤长度和色散系数是可以“预先扫描”的二十分钟就能跑完。我建议你把L设成变量写一个循环从5公里到80公里分别仿真然后把不同距离下的最优BER画成一张折线图。这条曲线能非常直观地告诉你这个均衡器方案的极限在哪里什么时候需要换更高阶的补偿算法。这种扫描思路放到论文里是一张很有说服力的图放到实际系统集成里更是必不可少的工程依据。