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) 在门操作相关频偏处低于某值,两者需联合评估。
测量中的关键约定与常见误区
相位噪声测量最容易出错的地方是参照约定不统一。ℒ(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),两者从不同侧面约束时钟源,缺一不可。