拿到这套10G QPSK光纤通信系统仿真工程也就是Matlab源码编号15103的那套时我的第一反应不是去看代码而是先问自己一个问题光纤通信做到10G这个速率为什么放着最简单的OOK强度调制不用非要绕一圈上QPSK这种相位调制这个问题想明白了后面看链路模型、看损伤补偿、看误码率曲线都会顺很多。这篇东西我打算从工程复现的角度拆开讲覆盖调制选型逻辑、光纤损伤建模方式、发射-信道-接收三段链路的实现细节、不同传输条件下的性能差异以及我自己调参过程中踩过的一些坑。对刚接触光通信仿真的朋友来说这是一条比较完整的上手路径已经跑过类似仿真的也可以重点看看第5节的参数敏感性问题。1. 为什么光纤通信到了10G反而爱用QPSK选型逻辑先讲清楚1.1 10G速率下OOK的瓶颈在哪里先回忆一下最朴素的强度调制OOK用光的有无表示0和1接收端直接探测光功率就行。10Gbps的OOK意味着符号速率就是10GBaud单极性归零或不归零脉冲的物理带宽大约要占10~20GHz这对1550nm窗口来说已经是一个相当宽的信号了。更麻烦的是光纤色散。标准单模光纤在1550nm附近的色散系数大约17ps/(nm·km)10Gbps信号如果用的是普通外调制器光源等效谱宽按0.1nm估算那么在80km的标准单模光纤上累积色散大约有1360ps/nm。这是什么概念5GBaud的QPSK符号周期是200ps1360ps意味着脉冲展宽已经超过符号周期的6倍会产生非常严重的码间干扰。OOK到这一步眼睛图基本就关没了哪怕加色散补偿光纤代价也不小。QPSK的优势说到底就是两个字降速。它一个符号带两个比特10Gbps只需要5GBaud的符号速率信号带宽直接减半。带宽减半带来两个连锁好处一是同样距离下的色散代价大幅下降二是同等带宽内能塞下更多波长通道频谱效率上去了。这也是经典WDM系统里10G时代后期和40G初期大量采用QPSK及衍生格式的根本原因。1.2 QPSK的四象限表示与复包络仿真从数学上看QPSK就是把两路比特分别调制到I路和Q路每个符号点落在复平面的四个象限上。在Matlab仿真里我们不需要真的在光频率上做运算而是用复基带包络来表示光场这样计算量小得多物理含义也足够清晰。% QPSK映射00-0, 01-1, 11-2, 10-3 bits randi([0 1], Nbits, 1); symbols reshape(bits, 2, []).; idx bi2de(symbols, left-msb); % 四个星座点等间隔分布在单位圆上 qpsk_sym exp(1j * (pi/4 idx * pi/2));这段代码有一个关键点星座点默认能量归一化每个符号的幅度是1。仿真里后面如果要计算信噪比、OSNR必须把能量基准搞清楚不然误码率曲线会对不上号。1.3 为什么从“强度”转向“相位”后系统鲁棒性更好除了带宽优势QPSK还有一个容易被忽略的特性理想情况下它是一个恒包络调制。什么意思每个符号的幅度都是1光功率不会在符号切换时大起大落。光纤非线性中最常见的SPM自相位调制作用强度直接和瞬时功率相关。OOK那种功率从0跳到满幅的信号在每个上升沿和下降沿都会产生强烈的相位啁啾而QPSK幅度恒定非线性相位扰动会被平均掉很大一部分。这一点在入纤功率稍高时体现得特别明显后面说功率扫描的部分我会再展开。2. 真实光纤损伤的建模方法不是把效应全堆上去就叫真实2.1 衰减与ASE噪声OSNR才是第一指标光纤衰减本身不难建模0.2dB/km100km下来就是20dB损耗。难的是光放大器带来的ASE噪声。每经过一个EDFA信号里都会叠加自发辐射噪声它是宽带的加性高斯白噪声。工程上判断系统好坏最常用的还不是单纯的信噪比而是OSNR光信噪比因为它跟测量带宽有关。ASE噪声功率在单偏振下的典型公式是P_ase 2 * n_sp * h * f * B_ref其中n_sp是自发辐射因子h是普朗克常数f是光频率B_ref通常取0.1nm对应的12.5GHz。仿真里最简单且实用的做法是直接在接收端加复高斯白噪声用OSNR反推噪声方差。这样能快速扫出一条BER-OSNR曲线用于验证DSP算法。如果要模拟真实级联链路就需要每80km放一个放大器每段加一次ASE复杂度会明显上去。2.2 色度色散频域传递函数的写法色散是线性效应在频域里写传递函数非常方便。单模光纤的色散传递函数可以写成H(ω) exp(j * D * λ² * ω² * L / (4πc))其中D是色散系数λ是中心波长ω是角频率L是光纤长度。在Matlab里实现就是在频域乘上这个相位因子freq (-Nfft/2 : Nfft/2-1) / Nfft * fs; omega 2 * pi * freq; Hcd exp(1j * D * lambda^2 * omega.^2 * L / (4 * pi * c)); E_out ifft(fft(E_in) .* Hcd);写这段代码时有个容易踩的坑fft之后频谱的顺序和freq数组的排列必须一一对应否则相位因子会加反。我建议先画出Hcd的实部虚部确认一下趋势再丢进链路里跑。2.3 非线性效应什么时候非做不可如果只是验证QPSK的线性性能不加非线性完全没问题。但如果要分析系统在不同入纤功率下的表现SPM等非线性效应就必须进模型。SPM的本质是光场的瞬时功率通过光纤的三阶非线性效应改变自身相位表现为一个和功率成正比的相位旋转。精确模拟非线性需要分步傅里叶法把光纤分割成很多小段每段先走一半色散再走一步非线性再走一半色散。简化写法如下% 非线性步长 E E .* exp(1j * gamma * dz * abs(E).^2); % 色散步长 E ifft(fft(E) .* Hcd(dz));gamma通常取1.2~1.3 W^-1·km^-1。关键参数是光纤截面积对应的有效面积决定非线性强度。仿真里如果你把步长设置得太粗非线性相位估计就会很不准设置得太细计算量又受不了。一般功率在5dBm以下时步长取1km足够如果扫到10dBm以上建议把步长缩到0.2km以内再对比一次结果。2.4 偏振相关损伤单偏振仿真的局限性这套标题里的系统如果只建单偏振QPSK偏振模色散PMD其实没有完整建模。PMD会让两个正交偏振分量以不同速度传播对单偏振信号的影响体现在偏振旋转和波形畸变仿真里可以用随机DGD来近似。但如果是单偏振仿真我建议做法是先用理想单偏振把核心算法跑通再在扩展阶段加入PMD一上来就把所有效应都堆上去出了问题根本不知道是谁导致的。3. 仿真链路三段式实现从比特流到星座图的完整路径3.1 发射端映射、脉冲成型与过采样的设置发射端做的事情很清晰随机比特 - QPSK符号 - 过采样 - 脉冲成型。过采样率每个符号的采样点数建议设16这是为了后面色散和非线性分步仿真留足够的频率余量。脉冲成型用根升余弦滤波器滚降因子选0.2~0.35之间比较常见。滚降因子越小频谱效率越高但对定时误差越敏感实际系统里0.25是个比较平衡的取值。sps 16; % 每符号16个采样点 upSym upsample(qpsk_sym, sps); % 插零过采样 rrcFilter rcosdesign(0.25, 10, sps, sqrt); % 根升余弦 txSignal filter(rrcFilter, 1, upSym);这里有一个重要的配合关系发射端用根升余弦接收端也必须用同样的根升余弦做匹配滤波两者级联才是完整的升余弦特性才能做到无码间干扰采样。只发端滤波、收端不做匹配误码率会很难看。3.2 信道部分把损伤按“确定性/随机性”分层加我的推荐顺序是先加衰减和色散这类确定性损伤再加非线性最后加ASE噪声。也就是说先把信号走一遍光纤的线性、非线性过程得到纯波形损伤后的结果再叠加噪声。这个顺序和真实物理过程有一定差异但对于大多数算法验证场景是够用的。如果你要做级联放大器的精确模拟那就是另一套逻辑每跨段结束先放大再加该段的ASE噪声然后进下一段。这个完整模型计算量大一般放在最后精调阶段用前期用末端加噪声的方式扫参效率更高。3.3 接收端的DSP链CD补偿、匹配滤波、相位恢复接收端我习惯拆成四个模块色散补偿、匹配滤波、采样同步、相位估计。色散补偿可以做一个频域逆滤波直接把第2节里Hcd的相位取共轭乘回去就能把80km的色散全部补回来。在Matlab里这属于“理想DSP”实际系统还有信道估计误差但仿真里先做到理想再往里面加扰动是合理的推进方式。QPSK相位恢复最常用的是四次方估计法也叫Viterbi-Viterbi算法% 对接收符号做四次方去除调制信息 phEst angle(mean(rxSym .^ 4)) / 4; rxSym rxSym .* exp(-1j * phEst);这里有一个绕不开的工程细节四次方把相位映射到了0~π/2的范围内所以恢复出来的相位有π/2的模糊直接判决会得到随机旋转的星座。解决办法一般是发端做差分编码收端判决后进行差分译码或者用已知训练序列校准这四种可能。4. 不同条件下的性能对比照着仿真结果说人话4.1 背靠背与80km传输OSNR代价从哪里来先说背靠背也就是发射机直接连接收机不加光纤只加ASE噪声。这种情况下损伤只有噪声误码率曲线接近理论值。5GBaud的QPSK理论上达到1e-3误码率需要的OSNR大约在8~9dB附近实际仿真加上滤波器非理想性会在10dB左右。加上80km标准单模光纤后如果不做色散补偿纯色散导致的码间干扰会让误码率即使开大OSNR也降不下去。这时候跑出来的BER曲线可能停在1e-2级别。加上数字色散补偿后80km的性能可以非常接近背靠背代价只有零点几dB主要来自色散补偿对带外噪声的放大。场景OSNRBER1e-3 近似值主要代价来源背靠背10 dB滤波器和匹配误差80km无CD补偿不可达码间干扰完全主导80kmDSP色散补偿10.5 dB噪声增强、残余CD4.2 入纤功率扫描你能找到最优功率点这是我最喜欢做的一个实验。把发射功率从-5dBm扫到10dBm固定传输距离和其他参数记录BER变化你会发现曲线是个U形。低功率端信号太弱ASE噪声主导BER很差随着功率升高OSNR变好BER逐步下降但功率超过某个点后SPM引起的非线性相位噪声开始增加BER又开始往上涨。这个U形的最低点就是该系统在当前配置下的最佳入纤功率。常见QPSK单通道系统里这个点大概在0~5dBm之间具体和色散累积量、跨段长度都有关系。入纤功率表现-5 dBmASE主导BER偏高2 dBm接近最佳工作点10 dBm非线性损伤明显BER回升这个扫描的意义在于它直观地说明了为什么光纤通信系统不能单纯靠加大功率来改善性能。同时它也是检验非线性建模是否正确的好方法——如果功率升高后BER纹丝不动说明你的非线性步长没起作用或者根本没有实现。4.3 残余色散对系统的影响补偿后也要留余量数字色散补偿不是万能的实际中会有残余色散。仿真里可以直接控制残余色散量比如名义上补1000ps/nm但实际链路是1100ps/nm就有100ps/nm的失配。从仿真数据看QPSK在5GBaud下对残余色散有一定容忍度200ps/nm以内基本无感知500ps/nm开始出现0.5dB左右的OSNR代价超过1000ps/nm代价就明显了。这说明在系统设计时色散补偿精度不用追求到零但要控制在符号周期的1/4以内。4.4 光滤波带宽的影响真实系统里接收端会有光滤波器限制ASE噪声。仿真里如果直接忽略光滤波器只用电源滤波器去滤噪声得到的BER会比实际系统乐观一点。把光滤波器带宽设成25GHz、30GHz、50GHz分别跑一遍你会发现窄带滤波器能压低ASE噪声但太窄会损伤信号高频分量。5GBaud QPSK加0.25滚降理论带宽是6.25GHz单个信号放25GHz光滤波里绰绰有余选30GHz是兼顾未来扩展和噪声抑制的常规做法。5. 调参最容易翻车的五个细节我踩过的坑都在这5.1 比特数不足误码率曲线在“抖”这是新手最容易忽略的问题。BER1e-4意味着平均每1万比特错1个你如果只仿真2万个比特运气好可能一个误码都没碰到BER算出来是0运气差点可能碰到两个得到1e-4。但这些数字波动很大根本不可信。工程上一般要求至少要观察到20个左右的误码再统计BER才比较稳。反推需要的比特数大约是20/BER。要测1e-4就得跑至少20万到40万个比特。对应到仿真里5GBaud、每符号16个采样点信号长度是40万×16640万个点这个量级Matlab跑起来并不慢但如果你把每个采样点的非线性分步步长设得很小就会开始吃时间。5.2 采样率和滤波器滚降别随便给过采样率太低会造成频域混叠尤其在色散补偿阶段高频部分被污染后补偿效果会很差。我建议每个符号至少8个点稳定起见16个点。滤波器滚降因子也不要乱用标准值以外的数滚降为0虽然频谱效率最高但RRC滤波器对定时误差的敏感度会急剧上升仿真里只要抽样时刻偏一点误码率立刻恶化。5.3 相位恢复的块长不能一刀切Viterbi-Viterbi算法需要把接收符号分成若干个块每个块单独估计平均相位。块长太短ASE噪声在四次方后放大4倍估计出来的相位抖动很大块长太长相位中的快速变化分量被平均掉对于有激光相位噪声的系统会低估损伤。经验上5GBaud QPSK取32~64个符号一块比较合适。你可以在仿真里跑一个块长扫描看哪一段BER最低把这个图存档方便以后参考。5.4 频域补偿注意循环边界用fft做色散补偿时如果直接把整帧数据丢进fft再ifft回来会引入循环卷积效应具体表现就是帧的首尾会有异常的误码集中点。解决办法是数据分段后做重叠保留法overlap-save或者干脆在帧头和帧尾各留一段保护比特统计误码时把保护段排除掉。提示我在第一次做80km色散补偿时总BER看起来合理但画误码位置图发现首尾各有一撮错误这就是循环边界问题。切掉保护段后数值立刻稳定了。5.5 不要把“加了ASE噪声”简单等同于“真实OSNR”ASE噪声是宽带均匀叠加的但接收机里的滤波器会切掉带外噪声所以OSNR测量点是影响最终BER的直接因素。仿真中如果只关心误码率可以直接在接收符号上按目标SNR加噪声但如果要和真实系统对标必须把OSNR换算清楚。0.1nm光带宽对应约12.5GHz信号带宽只有10GHz左右两者之间有一个约0.97dB的差异把这个折算进噪声方差结果才严谨。6. 从这套QPSK往上升级的三个方向6.1 PDM-QPSK多一个偏振维度信息速率直接翻倍在QPSK基础上把两路独立信号分别放在两个正交偏振态上传输就得到PDM-QPSK。同样5GBaud符号速率信息速率变成20Gbps。代价是接收端需要做偏振解复用通信算法里加一个自适应均衡器比如经典的恒模算法去分离两路信号。这套仿真如果要扩展成100G系统PDM-QPSK其实就是最经典的初代方案模板。6.2 16QAM同带宽传更多比特但要接受OSNR代价从QPSK换到16QAM每个符号带4个比特5GBaud符号速率可以做到20Gbps。星座点间距变小同样的OSNR下误码率会明显上升。粗略估算16QAM比QPSK需要多约7dB的OSNR才能达到同样的BER。这就引出一个现实问题调制阶数不能无限往上走高阶调制必须配合更低的噪声系数、更强的FEC纠错能力以及更完善的相位噪声抑制算法。6.3 迭代开发每次都只改一个变量我给所有想上手这套源码的朋友一个建议拿到代码后先别急着把所有功能打开。第一步跑背靠背确认误码率曲线和理论曲线吻合第二步加色散确认补偿算法有效第三步加ASE确认OSNR关系第四步再开非线性。每一步只改一个变量记录一张结果图。这样出了问题你能立刻锁定是哪一层加的损伤导致性能恶化。如果第一次就跑全链路数据一出问题你会陷入“是色散算错了还是非线性步长不对还是相位恢复失效”的三连猜排查成本高得多。这也正好呼应了本节的主题——这套QPSK仿真的价值不只是复现一个通信系统它本身就是一个调试光通信DSP算法的良好试验台。