4. 相位噪声的频域与时域表征方法

相位噪声本质上是振荡器输出相位随时间随机起伏的现象。同一个物理过程,既可以在频域用功率谱密度描述,也可以在时域用抖动或方差统计描述。工程上之所以需要两套语言,是因为不同器件、不同链路环节对相位噪声的敏感方式不同:混频器、锁相环环路滤波器更关心某个频偏处的谱密度;ADC 采样时钟、数字接口更关心某个积分带宽内的总抖动;而原子钟、频率源长期稳定度则必须用阿伦方差或哈达玛方差来评估。掌握这些表征量之间的换算与测量原理,是量子测控时钟链路指标分解的基础。

频域表征:从功率谱密度到单边带相位噪声

设振荡器输出为

v(t) = V₀ · [1 + α(t)] · cos[2πf₀t + φ(t)]

其中 α(t) 为幅度起伏,φ(t) 为相位起伏。在多数低噪声振荡器中幅度噪声远小于相位噪声,可近似忽略 α(t)。相位噪声的频域表征有两种常用形式:相位起伏的功率谱密度 Sφ(f),单位为 rad²/Hz;以及单边带相位噪声 ℒ(f),定义为偏离载波 f 处 1 Hz 带宽内的单边带功率与载波功率之比,单位为 dBc/Hz。二者关系为

ℒ(f) = 10·log₁₀[ Sφ(f) / 2 ]   (dBc/Hz)
Sφ(f) = 2 · 10^(ℒ(f)/10)        (rad²/Hz)

注意这里的 f 是频偏(offset frequency),不是绝对频率。ℒ(f) 曲线通常呈现几个特征斜率区:近载波区约 −30 dB/dec(对应 1/f³ 相位噪声,源于谐振器 Q 值有限和闪烁噪声上变频),中段约 −20 dB/dec(1/f²,白频率噪声),远端约 0 dB/dec(白相位噪声,由热噪声基底决定)。测量 ℒ(f) 的经典方法是鉴相法:用参考源与被测源正交鉴相,低通滤波后由 FFT 分析仪读取基带相位噪声谱,再按上式折算。现代仪器(如相位噪声分析仪)多采用互相关双通道结构,通过两路独立鉴相器抑制本底噪声,可把测量本底压到 −180 dBc/Hz 以下。

时域表征:相位抖动与积分带宽

相位抖动是把频域相位噪声在指定带宽内积分后得到的时域量。对同一相位噪声谱,积分带宽不同,抖动值差异可能达数倍,因此工程上必须明确积分区间。常用的 RMS 相位抖动(弧度)为

φ_rms = sqrt( 2 · ∫[f₁→f₂] 10^(ℒ(f)/10) df )   (rad)

换算为时间抖动:

t_rms = φ_rms / (2πf₀)   (s)

例如某 100 MHz 时钟,ℒ(f) 在 1 kHz~10 MHz 积分带宽内积分得到 φrms ≈ 0.35 rad,则 trms ≈ 0.35/(2π×10⁸) ≈ 0.56 ps。若把积分下限改为 10 kHz,抖动会明显下降。因此 ADC 采样时钟选型时,必须核对厂商给出的抖动是在哪个积分带宽下测得的,否则指标不可比。积分带宽的选取原则是:下限取环路带宽或系统关心的最短观测时间对应频率,上限取后级电路实际带宽或采样率的一半。

阿伦方差与哈达玛方差

频域谱密度适合描述平稳随机过程,但振荡器还存在频率漂移、闪烁噪声等非平稳成分,长期稳定度必须用 Allan 方差 σy²(τ) 描述。它通过相邻两个采样区间平均频率的相对偏差来定义:

σ_y²(τ) = (1/2) · ⟨ (ȳ_{k+1} − ȳ_k)² ⟩

其中 ȳk 是第 k 个 τ 区间内的平均相对频率偏差。Allan 方差对白频率噪声、闪烁频率噪声、随机游走频率噪声分别呈现 τ−1/2、τ⁰、τ+1/2 的斜率,因此从 σy(τ) 曲线斜率可以反推噪声类型。但 Allan 方差对频率漂移(线性漂移)不收敛,且无法区分白相位噪声与闪烁相位噪声。Hadamard 方差用三次差分代替二次差分,对线性漂移不敏感,且能区分上述两种相位噪声,适合评估存在明显频率漂移的量子测控时钟源。

实际测量时,常用频率计数器记录频率序列,再按 τ = m·τ₀ 做重叠采样计算 Allan 偏差。重叠采样比非重叠采样能利用更多数据点,在相同测量时长下给出更平滑的曲线。判断时钟源是否满足量子比特门操作要求时,通常要求在某 τ 处 σy(τ) 低于门时间对应的阈值,同时 ℒ(f) 在门操作相关频偏处低于某值,两者需联合评估。

相位噪声频域—时域表征关系图 频域表征 Sφ(f) 相位功率谱密度 rad²/Hz ℒ(f) 单边带相位噪声 dBc/Hz ÷2取对数 特征斜率:−30 / −20 / 0 dB/dec 测量:鉴相法 / 互相关双通道 本底可低至 −180 dBc/Hz f 为频偏,非绝对频率 时域表征 φ_rms / t_rms 相位抖动 / 时间抖动 σ_y(τ) / H_y(τ) Allan / Hadamard 偏差 积分带宽 [f₁, f₂] 必须明确 t_rms = φ_rms / (2πf₀) Allan 斜率判别噪声类型 Hadamard 抗频率漂移 积分 谱估计 工程测量与判断要点 鉴相器 + FFT 分析仪 互相关抑制本底 频率计数器记录序列 重叠采样计算方差 核对积分带宽 统一 f₀ 与 τ 定义 联合评估 ℒ(f) 与 σ_y(τ) 匹配量子门时间要求 示例:100 MHz 时钟抖动换算 ℒ(f) 在 1 kHz~10 MHz 积分 → φ_rms ≈ 0.35 rad → t_rms ≈ 0.56 ps 若积分下限改为 10 kHz,抖动显著下降 → 指标必须标注积分带宽

测量中的关键约定与常见误区

相位噪声测量最容易出错的地方是参照约定不统一。ℒ(f) 的定义中,分母是载波功率,分子是单边带功率,因此同一物理谱用 Sφ(f) 和 ℒ(f) 表示时数值相差 3 dB(即 10log₁₀2)。若把 Sφ 直接当作 ℒ 使用,抖动积分会偏大 √2 倍。另一个常见问题是积分带宽的上下限选取:下限若取 0,积分会因 1/f 噪声发散,必须从某个有限频偏起积;上限若超过测量本底,会把仪器噪声也积进去,导致抖动虚高。工程上应取测量本底至少低于被测谱 10 dB 的频段作为有效积分区间。

Allan 方差与 Hadamard 方差的 τ 序列也必须与测量时长匹配。若总测量时长为 T,最大可用 τ 约为 T/3~T/5,否则统计样本太少,曲线尾部会出现剧烈波动。对于量子测控中常用的 10 MHz~1 GHz 时钟,通常要求同时给出 ℒ(f) 曲线(1 Hz~10 MHz 频偏)和 σy(τ) 曲线(1 ms~1000 s),两者从不同侧面约束时钟源,缺一不可。

频域与时域不是两套独立指标,而是同一随机过程的不同投影。做指标分解时,先确定系统对哪个域敏感,再选择对应表征量,最后用换算关系校验两域指标是否自洽。