4. 频谱分析仪与FFT基础

在量子芯片噪声谱测量中，频谱分析仪和FFT（快速傅里叶变换）是两类最常用的谱估计工具。频谱分析仪通过模拟或数字中频滤波逐点扫描获得频谱，适合观测窄带信号和快速瞬态；FFT则对采样后的时域波形做离散傅里叶变换，适合宽带、相参、可重复采集的噪声谱。无论使用哪种工具，最终得到的功率谱密度（PSD）估计都受四个关键参数支配：分辨率带宽（RBW）、窗函数、平均次数和频率跨度。理解它们之间的耦合关系，是判断谱线可信度和测量不确定度的前提。

分辨率带宽：谱估计的频率刻度

分辨率带宽决定频谱在频率轴上能分辨的最小间隔。对FFT而言，若采样率为 f_s，采样点数为 N，则频率分辨率（bin宽度）为 Δf = f_s/N。例如 f_s = 1 MSa/s、N = 1 048 576 时，Δf ≈ 0.954 Hz。对频谱分析仪，RBW由中频滤波器带宽决定，通常以1-3-10序列设置，如1 Hz、3 Hz、10 Hz、30 Hz、100 Hz。RBW越窄，频率分辨能力越强，但扫描时间越长，且单次扫描内可平均的次数越少。

RBW与噪声本底的关系是：对于白噪声，显示平均噪声电平（DANL）随RBW以10log₁₀(RBW)变化。RBW从100 Hz降到10 Hz，噪声本底下降10 dB。因此比较不同RBW下的谱线绝对幅度没有意义，必须归一化到1 Hz带宽，即得到功率谱密度，单位dBm/Hz或V²/Hz。

窗函数：抑制频谱泄漏

FFT默认假设时域信号在观测窗口内是周期性的。若信号两端不连续，截断处会产生频谱泄漏，使强信号的旁瓣抬高远端噪声本底。窗函数通过在时域两端施加平滑衰减来抑制泄漏，代价是主瓣展宽、频率分辨率下降。

常用窗函数的选择依据：

• 矩形窗：主瓣最窄，频率分辨率最高，但旁瓣仅-13 dB，适合整周期采样或瞬态信号。
• 汉宁窗：旁瓣-31 dB，主瓣宽度为矩形窗2倍，适合连续噪声谱的通用测量。
• 平顶窗：主瓣很宽，但幅度平坦度极好，适合精确测量单频幅度，不适合噪声本底。
• 布莱克曼-哈里斯窗：旁瓣-92 dB，适合大动态范围下观测弱噪声边带。

使用窗函数后，PSD估计需除以窗函数的等效噪声带宽（ENBW）进行归一化。汉宁窗的ENBW约为1.5个bin，矩形窗为1.0个bin。若忽略ENBW修正，PSD幅度会系统性偏高或偏低。

平均次数：降低估计方差

单次FFT得到的周期图是功率谱的随机估计，其标准差与均值同量级，谱线起伏剧烈。通过多次平均可降低方差。对于M次非相干平均，估计的标准差约降为原来的1/√M。例如平均100次，幅度起伏从约100%降到约10%，即约20 dB的平滑改善。

平均方式有两种：线性平均（对功率谱直接平均）和对数平均（对dB值平均）。线性平均是无偏的，适合后续积分计算；对数平均会低估噪声本底约2.5 dB，但视觉上更平滑。量子芯片噪声谱测量推荐使用线性平均后再转dB。

平均次数与测量时间的关系为：T_total = M × N / f_s。若 f_s = 1 MSa/s、N = 1M、M = 100，则总采集时间约100 s。在量子比特退相干时间有限的实验中，需在稳定性和统计精度之间折中。

功率谱密度估计与不确定度

PSD的物理意义是单位带宽内的功率。对电压噪声 v(t)，单边PSD定义为 S_V(f) = 2|V(f)|²/(f_s · N · ENBW)，单位V²/Hz。若用dBm/Hz表示，则 S_dBm = 10log₁₀(S_V/R/1 mW)。

测量不确定度主要来自三部分：统计不确定度（由平均次数决定，约1/√M）、系统不确定度（窗函数ENBW修正、前端增益平坦度、参考电平精度）和本底不确定度（频谱仪自身噪声贡献）。对于量子芯片噪声谱，通常要求统计不确定度小于0.5 dB，对应平均次数不少于20次；若要求0.1 dB，则需约200次。

参数配置示例与结果判读

以下给出一组用于超导量子比特噪声谱测量的FFT参数配置示例，参数为教学假设值，非实测结果：

采样率 fs      = 1 MSa/s
采样点数 N     = 1,048,576
频率分辨率 Δf  = 0.954 Hz
窗函数         = 汉宁窗 (ENBW = 1.5 bin)
平均次数 M     = 100 (线性平均)
总采集时间     ≈ 105 s
频率跨度       = DC ~ 500 kHz (单边谱)
归一化         = S_V(f) = 2|V(f)|² / (fs · N · ENBW)

判读要点：若在100 kHz处观察到1/f噪声拐点，需确认该频率远大于Δf，否则拐点位置会被频率分辨率模糊。若谱线在10 kHz以下出现明显上翘，可能是直流漂移或窗函数泄漏，应检查是否使用了汉宁窗并去除了DC分量。若平均100次后谱线仍起伏超过±1 dB，说明统计不确定度尚未收敛，需增加平均次数或检查触发是否稳定。

测量不确定度的入门估计可按以下方式合成：统计不确定度 u_stat ≈ 4.34/√M dB（对应幅度标准差），系统不确定度 u_sys 取0.2~0.5 dB（取决于校准水平），本底不确定度 u_floor 取0.1~0.3 dB。合成标准不确定度 u_c = √(u_stat² + u_sys² + u_floor²)。当 M = 100 时，u_stat ≈ 0.43 dB，若 u_sys = 0.3 dB、u_floor = 0.2 dB，则 u_c ≈ 0.57 dB。扩展不确定度取 k = 2，约1.1 dB。这一量级在量子芯片噪声谱比较中通常可接受，但若需分辨小于1 dB的谱差异，则需将平均次数提高到1000次以上，或采用交叉谱估计抑制前端噪声。

总结本节的核心关系：RBW决定频率分辨率和噪声本底，窗函数决定泄漏抑制与ENBW修正，平均次数决定统计不确定度，三者共同决定PSD估计的偏差与方差。在实际测量中，应先根据待测噪声的频率范围和幅度动态范围选择RBW和窗函数，再根据所需不确定度确定平均次数，最后统一归一化到1 Hz带宽，才能得到可比较、可复现的量子芯片噪声谱。