3. 光纤传输物理效应建模

在光通信链路仿真中，光纤不是一根理想导线，而是一个同时存在衰减、色散和非线性效应的分布参数系统。要准确预测长距离传输后的眼图、误码率或星座图劣化，必须把这些物理效应写成统一的数学模型，再用数值方法求解。本节围绕三个核心问题展开：损耗、色散、非线性效应各自如何用数学描述；它们如何合并为非线性薛定谔方程（NLSE）；以及如何用分步傅里叶法（SSFM）求解并设置典型仿真参数。

一、三类物理效应的数学描述

损耗表现为光功率沿传播方向指数衰减。若以 P(0) 为入射功率，则 P(z) = P(0)·exp(−αz)，其中 α 为功率衰减系数（单位 1/km）。工程上更常用 dB/km 表示：α_dB = 10·α·log₁₀e ≈ 4.343α。例如标准单模光纤在 1550 nm 处约 0.2 dB/km，对应 α ≈ 0.046 /km。在 NLSE 中，损耗通常以复包络幅度衰减项 −α/2·A 出现，因为功率正比于 |A|²。

色散源于不同频率成分的群速度不同。用泰勒展开传播常数 β(ω)，一阶项 β₁ 只带来时延，二阶项 β₂（单位 ps²/km）导致脉冲展宽，三阶项 β₃ 影响高阶色散。工程中更常用色散参数 D = −2πc/λ²·β₂，标准单模光纤在 1550 nm 处 D ≈ 17 ps/(nm·km)，对应 β₂ ≈ −21.7 ps²/km。色散在时域表现为脉冲展宽，在频域表现为相位随频率二次变化。

非线性效应主要来自克尔效应：折射率随光强变化 n = n₀ + n₂·I。由此产生三类典型效应：自相位调制（SPM）——脉冲自身强度调制自身相位，导致频谱展宽；交叉相位调制（XPM）——其他信道的光强改变本信道相位，在 WDM 系统中造成信道间串扰；四波混频（FWM）——多个频率成分相互作用产生新频率，满足 ω₁+ω₂=ω₃+ω₄，在零色散附近尤为严重。三者统一由非线性系数 γ = 2πn₂/(λ·A_eff) 描述，单位 1/(W·km)，典型值约 1.3 /(W·km)。

二、非线性薛定谔方程

把损耗、色散和非线性合并，得到慢变包络 A(z,t) 满足的 NLSE：

∂A/∂z = −(α/2)A − (iβ₂/2)·∂²A/∂t² + (β₃/6)·∂³A/∂t³ + iγ|A|²A

式中各项依次为：损耗项、二阶色散项、三阶色散项、SPM 非线性项。若考虑 WDM 多信道，非线性项需扩展为 iγ(|A|² + 2Σ|A_k|²)A，其中 2 倍因子来自 XPM，而 FWM 则通过多信道耦合项体现。该方程一般没有解析解，必须数值求解。

三、分步傅里叶法求解

SSFM 的核心思想是：把光纤分成若干小段 Δz，在每段内把色散算子和非线性算子交替作用。色散在频域是简单乘法，非线性在时域是简单乘法，因此每步做一次 FFT 和一次 IFFT 即可。对称分步傅里叶法精度更高：

A(z+Δz) ≈ exp(Δz/2·D̂)·exp(Δz·N̂)·exp(Δz/2·D̂)·A(z)

其中 D̂ 为色散算子（频域乘子），N̂ 为非线性算子（时域乘子）。实现时，先做半步色散，再做一步非线性，最后再做半步色散。步长 Δz 需满足：Δz ≪ 1/(γ·P_peak) 且 Δz ≪ L_D（色散长度），通常取 0.1～1 km 量级，或采用自适应步长。

四、典型仿真参数设置与结果判读

以一个 80 km 标准单模光纤、单信道 10 Gb/s NRZ 信号为例。设入纤峰值功率 5 mW，波长 1550 nm，α = 0.2 dB/km，D = 17 ps/(nm·km)，γ = 1.3 /(W·km)。色散长度 L_D = T₀²/|β₂|，若脉冲半宽 T₀ = 50 ps，则 L_D ≈ 115 km；非线性长度 L_NL = 1/(γP) ≈ 154 km。两者与 80 km 同量级，说明色散与非线性均不可忽略。

仿真时，时窗取 12.8 ns（对应 128 个比特），采样点数 2¹⁴ = 16384，频率分辨率约 78 MHz。步长 Δz 先取 0.5 km，共 160 步。运行后观察三项指标：一是输出脉冲时域宽度，若明显大于输入，说明色散主导；二是输出频谱宽度，若比输入宽，说明 SPM 引起频谱展宽；三是若开启多信道，检查是否出现新频率分量，判断 FWM 强度。若把 Δz 减半后输出频谱变化小于 1%，可认为步长已收敛。

综上，损耗、色散和非线性效应分别对应 NLSE 中的衰减项、色散项和非线性项，SSFM 通过频域-时域交替计算实现高效求解。参数设置的关键是让步长、时窗和采样率与信号带宽、色散长度、非线性长度匹配，并通过收敛性检查确认结果可信。