4. 光纤信道损伤建模

相干接收机把光场的幅度与相位都搬到基带电域之后，光纤信道对信号的全部“伤害”也随之显式暴露出来。与强度调制直接检测不同，相干系统不能只用一个实系数衰减来描述信道，而必须建立复基带模型：色散是频率相关的相位弯曲，偏振模色散是随频率变化的偏振旋转，非线性效应是功率相关的相位与偏振耦合，激光器相位噪声与频率偏移则是收发本振引入的随机相位与线性相位斜坡。本节把这些损伤统一写成可仿真的数学形式，并给出工程上可直接落地的离散化方法。

一、统一复基带信道模型

设发射端双偏振复基带信号为 E(t)=[Ex(t), Ey(t)]^T，忽略绝对光频后，接收端基带场可写成：

E_rx(t) = H_CD(ω) · H_PMD(ω) · exp(j·φ_NL(t)) · exp(j·(Δω·t + φ_L(t))) · E(t) + n(t)

其中 H_CD 为色散传递函数，H_PMD 为偏振模色散琼斯矩阵，φ_NL 为非线性相位，Δω 为收发激光器频率偏移，φ_L 为激光器相位噪声，n(t) 为ASE噪声。这个顺序并非唯一，但把线性损伤放在前面、非线性相位放在中间、载波相位放在最后，便于分步验证与DSP模块对应。

二、色散的数学描述

色散源于不同频率分量在光纤中传播速度不同。设色散系数 D（单位 ps/(nm·km)），长度 L，则群速度色散传递函数为：

H_CD(ω) = exp( -j · (D·λ²·L)/(4πc) · ω² )

其中 ω 为相对中心频率的角频率偏移，λ 为波长，c 为光速。以标准单模光纤 D=17 ps/(nm·km)、λ=1550 nm、L=80 km 为例，色散系数约为 1360 ps/nm，在 32 GBaud 信号上会产生显著码间干扰。工程上常用频域相乘实现：对每块信号做FFT，乘以 H_CD(ω)，再IFFT。需要注意块长度要大于色散冲激响应长度，否则循环卷积会引入块间干扰。

三、偏振模色散

PMD 由光纤双折射随机变化引起，表现为两个主偏振态之间的差分群时延 Δτ。一阶PMD可用琼斯矩阵描述：

H_PMD(ω) = R(θ) · diag( exp(+j·ω·Δτ/2), exp(-j·ω·Δτ/2) ) · R(θ)^H

其中 R(θ) 为偏振旋转矩阵，Δτ 为DGD。典型单模光纤 PMD 系数约 0.1 ps/√km，80 km 链路一阶DGD约 0.9 ps，对 32 GBaud 信号影响有限；但长距离或高波特率下需考虑二阶PMD。仿真时可用多个级联的随机旋转与差分延迟段来近似。

四、非线性效应

克尔非线性使折射率随光强变化，导致自相位调制、交叉相位调制和四波混频。对双偏振信号，忽略四波混频时可用 Manakov 方程描述：

∂E/∂z = -α/2·E - j·β2/2·∂²E/∂t² + j·γ·(8/9)·|E|²·E

其中 γ 为非线性系数，典型值约 1.3 /(W·km)，|E|² 为总功率。工程仿真常用分步傅里叶法：把光纤分成若干小段，每段内交替执行色散频域运算和非线性时域相位旋转。步长需足够小，使每段非线性相移小于约 0.05 rad。

五、激光器相位噪声与频率偏移

激光器相位噪声通常建模为维纳过程：

φ_L[n] = φ_L[n-1] + σ_p · w[n],  w[n] ~ N(0,1)

其中 σ_p² = 2π·Δν·T，Δν 为激光器线宽，T 为采样周期。频率偏移则表现为线性相位斜坡：

φ_off[n] = 2π·Δf·n·T

两者叠加后，接收信号星座图会整体旋转并随时间扩散。DSP中通常先用粗频偏估计消除 Δf，再用载波相位恢复跟踪残余相位噪声。

六、离散化仿真实现要点

把上述模型落到代码时，需要统一采样率、块长度和功率归一化。建议按以下顺序搭建仿真链路：

1. 生成双偏振QPSK或16QAM符号，上采样到每符号2或4个采样点。
2. 频域施加 H_CD(ω)，注意 ω 以弧度/秒为单位，且FFT频率轴要按 fftshift 排列。
3. 用两个随机旋转矩阵和差分延迟级联实现PMD，每级 Δτ 按 0.1 ps/√km · √(L/N) 分配。
4. 分步傅里叶法施加非线性：每段先频域色散，再时域乘 exp(j·γ·(8/9)·|E|²·L_eff)。
5. 最后乘上 exp(j·(2π·Δf·nT + φ_L[n]))，其中 φ_L 用累积高斯随机变量生成。

七、典型参数与观察结果

以 32 GBaud 双偏振16QAM、80 km 标准单模光纤为例，设置 D=17 ps/(nm·km)、γ=1.3 /(W·km)、入纤功率 0 dBm、激光器线宽 100 kHz、频偏 200 MHz。仅开色散时，星座图呈明显旋转扩散；加入PMD后，两个偏振态出现不同延迟；加入非线性后，高功率符号的相位被额外旋转；最后叠加相位噪声与频偏，星座图整体旋转且随时间漂移。DSP按“频域色散补偿→自适应均衡→粗频偏→载波相位恢复”的顺序处理，可逐步恢复星座点。

掌握这四类损伤的数学描述后，后续DSP模块的设计就有了明确的“被补偿对象”。色散对应频域二次相位，PMD对应频率相关琼斯矩阵，非线性对应功率相关相位旋转，相位噪声与频偏对应时域相位累积与线性斜坡。把信道模型与DSP模块一一对应，是相干系统仿真与工程调试的基本功。