4、频率梳的数学描述：电场表达式、梳齿频率公式、傅里叶级数展开与噪声表示
光学频率梳在时域上表现为一系列等间隔的超短脉冲，在频域上则呈现为一系列离散且等间距的窄线宽纵模。要深入理解其工作原理并应用于精密测量，必须建立严格的数学模型。本章将从电场表达式出发，推导梳齿频率公式，并引入傅里叶级数展开以连接时域与频域，最后讨论相位噪声与幅度噪声的数学表示方法，为后续章节的锁相控制与噪声抑制奠定理论基础。
4.1 电场表达式：从时域脉冲到频域梳齿
光学频率梳的时域电场可以看作一个周期性脉冲序列。设脉冲包络为 A(t)，脉冲重复周期为 T_rep（对应重复频率 f_rep = 1/T_rep），则理想情况下电场表达式为：
E(t) = Σ_{m=-∞}^{∞} A(t - m·T_rep) · e^(i·2π·f_CEO·t) + c.c.其中 f_CEO 为载波包络偏移频率（Carrier-Envelope Offset Frequency），表示载波相位相对于脉冲包络在相邻脉冲间的相位滑移速率。c.c. 表示复共轭项。对 E(t) 进行傅里叶变换，利用傅里叶级数性质可得频域表达式：
E(ω) = Σ_{n=-∞}^{∞} E_n · δ(ω - 2π·f_n)其中 f_n = n·f_rep + f_CEO，E_n 为第 n 根梳齿的复振幅。这一结果揭示了频率梳的本质：时域周期脉冲对应频域等间隔离散谱线，谱线间隔严格等于 f_rep，整体频率偏移由 f_CEO 决定。
4.2 梳齿频率公式 f_n = n·f_rep + f_CEO 的物理内涵
梳齿频率公式是频率梳最核心的数学表达。其中 n 为整数（通常 n 的数量级为 10^5 ~ 10^6），f_rep 由腔长决定（典型值 100 MHz ~ 10 GHz），f_CEO 由脉冲内载波相位与包络相位差决定。在实际系统中，f_rep 和 f_CEO 均可通过锁相环精确锁定到微波参考或光学参考上。
| 参数 | 物理含义 | 典型值 | 测量方法 |
|---|---|---|---|
| f_rep | 脉冲重复频率，由腔长 L 决定：f_rep = v_g / (2L) | 100 MHz ~ 10 GHz | 光电探测 + 频谱分析仪 |
| f_CEO | 载波包络偏移频率，由群速度与相速度差决定 | 0 ~ f_rep | f-2f 自参考干涉法 |
| n | 梳齿序号，正整数 | 10^5 ~ 10^6 | 通过拍频测量确定 |
需要特别注意的是，f_CEO 的测量通常采用 f-2f 自参考技术：将频谱展宽至一个倍频程后，低频端第 n 根梳齿倍频后与高频端第 2n 根梳齿拍频，得到 f_CEO = 2·f_n - f_2n。这一方法使得 f_CEO 可以被直接锁定到射频参考上，从而实现全梳齿频率的稳定。
4.3 傅里叶级数展开与梳齿振幅分布
将时域电场展开为傅里叶级数，可以得到每根梳齿的振幅和相位信息。假设脉冲包络为高斯型 A(t) = A_0 · exp(-t²/2τ²)，则频域包络为：
E_n ∝ A_0 · √(2π)·τ · exp(-(2π·τ·(f_n - f_CEO))²/2)这表明梳齿振幅服从高斯分布，中心频率位于 f_CEO 处，宽度由脉冲宽度 τ 决定。脉冲越短，频域包络越宽，可用的梳齿数量越多。对于飞秒激光器（τ ≈ 100 fs），频域包络宽度可达数十纳米，覆盖数千根梳齿。
4.4 相位噪声与幅度噪声的数学表示
实际频率梳并非理想周期信号，而是受到噪声扰动。引入相位噪声项 δφ(t) 和幅度噪声项 ε(t) 后，电场表达式修正为：
E(t) = [1 + ε(t)] · Σ_m A(t - m·T_rep - δT(t)) · e^(i·2π·f_CEO·t + i·δφ(t)) + c.c.其中 δT(t) 表示脉冲时间位置的抖动，与相位噪声的关系为 δT(t) = δφ(t) / (2π·f_CEO)。在频域中，相位噪声导致梳齿线宽展宽，幅度噪声导致梳齿功率波动。通常用单边带功率谱密度 S_φ(f) 和 S_ε(f) 来表征：
S_φ(f) = |FT{δφ(t)}|² [rad²/Hz]
S_ε(f) = |FT{ε(t)}|² [1/Hz]对于锁模激光器，相位噪声主要来源于泵浦功率波动、腔长抖动和色散波动；幅度噪声主要来源于增益介质的自发辐射和泵浦噪声。在精密测量中，通常需要将 f_rep 的相位噪声抑制到 mHz/√Hz 量级，这要求采用高带宽锁相环和低噪声参考源。
4.5 案例拆解：飞秒激光频率梳的噪声表征与参数提取
假设某课题组搭建了一套基于掺铒光纤锁模激光器的光频梳系统，重复频率 f_rep = 250 MHz，中心波长 1560 nm。实验目标是通过 f-2f 干涉法测量 f_CEO，并评估系统的相位噪声性能。
初始现象：使用频谱分析仪观察光电探测器输出，发现 f_rep 信号清晰，但 f_CEO 信号较弱且不稳定，信噪比仅约 20 dB。
分析判断：f_CEO 信号弱可能原因包括：超连续谱展宽不足导致倍频效率低、探测器响应带宽不足、f_CEO 信号被其他噪声淹没。通过检查光谱，确认频谱已展宽至 1000~2200 nm，满足倍频程要求；探测器带宽为 1 GHz，足以响应 250 MHz 的拍频信号。
操作步骤：
- 优化超连续谱产生效率：调整泵浦功率至 300 mW，微调波片使偏振匹配，光谱展宽范围扩展至 950~2300 nm。
- 使用雪崩光电二极管（APD）替代普通 PIN 管，提高拍频信号幅度约 8 dB。
- 在 APD 后增加窄带带通滤波器（中心频率 250 MHz，带宽 10 MHz），抑制带外噪声。
- 将 f_CEO 信号接入锁相环，与射频参考源（铷原子钟）进行相位锁定。
结果对比：优化后 f_CEO 信噪比提升至 45 dB，锁相环闭环后 f_CEO 残余相位噪声在 1 kHz 偏移处为 -95 dBc/Hz，对应时间抖动约 50 as。同时测量 f_rep 的相位噪声，在 1 kHz 处为 -110 dBc/Hz，表明腔长稳定性良好。
风险与代价：提高泵浦功率增加了系统功耗和热漂移，需要增加温控措施；APD 的增益噪声可能引入额外幅度噪声，需在后续实验中评估。
复盘经验：频率梳噪声表征的关键在于分离 f_rep 与 f_CEO 的贡献。建议同时监测两个参数的功率谱密度，并利用互相关函数判断噪声源是否独立。对于光纤激光器，环境振动和温度漂移是低频噪声的主要来源，应优先做好隔振和温控。
4.6 噪声模型总结与工程应用要点
综合上述分析，频率梳的完整噪声模型可归纳为：
| 噪声类型 | 来源 | 频域特征 | 抑制方法 |
|---|---|---|---|
| f_rep 相位噪声 | 腔长抖动、泵浦功率波动 | 低频为主，1/f 噪声 | 高带宽锁相环、低膨胀系数腔体 |
| f_CEO 相位噪声 | 载波相位随机走动、非线性效应 | 白噪声基底 + 低频漂移 | f-2f 锁定、泵浦电流稳定 |
| 幅度噪声 | 自发辐射、泵浦噪声 | 宽频带，与增益带宽相关 | 平衡探测、功率稳定器 |
| 脉冲时间抖动 | 群速度色散波动 | 与 f_rep 噪声相关 | 色散补偿、被动锁模优化 |
在工程应用中，应根据具体测量需求确定噪声指标。例如，光学频率合成器需要同时锁定 f_rep 和 f_CEO；而低噪声微波产生只需锁定 f_rep，对 f_CEO 要求较低。理解噪声的数学表示是设计锁相环路滤波器的基础，也是评估系统极限性能的必要前提。