引　言

光学频率梳（optical frequency comb, OFC）因具有高精度和高光谱分辨率，已经成为宽带光谱学的一种新兴技术。光学频率梳在时域是等间距的激光脉冲序列，在频域中展示为等频率间距分布的梳齿线。这些梳齿具有间距相等、数量众多、光谱范围宽、梳齿宽度窄等特点，这些特点为获得高光谱分辨率奠定了基础。一直以来，由于光波段的频率远大于电子仪器工作带宽，导致没有一种直接的方法可以准确地在电子电路中获取光波段的频率和相位信息，而光学频率梳的发明^[1]在射频电子设备和光频段之间提供了一种新型的连接方式^[2]，从而解决了这一难题。

飞秒激光频率梳近年来在高精度光谱测量中贡献了重要力量，提供了许多新型光谱测量方法^[3-4]。其中基于光频梳的相干性，使用两个具有微小重复频率差的光频梳组成的双梳光谱技术进一步展示了光频梳在宽带光谱学的发展潜力。双梳光谱技术也因其高光谱采集速度、高频率分辨率和良好的稳定性等特点^[5]得到了广泛的应用^[6-13]。2002年，Schiller^[14]首次提出了双梳光谱的概念。2004 年，Keilmann等^[15]首次设计了实验，验证了其有效性和可行性。2005年，Janke等^[16]首次在太赫兹时域光谱系统中引入了两个重复频率差异较小的同步飞秒激光器，采用异步光学采样技术，使用电控延迟线代替传统机械延迟线完成光谱的采样，结果表明使用异步光学采样方法显著提高了太赫兹时域光谱仪的性能。自此以后，异步光学采样技术被广泛地应用于太赫兹时域光谱系统的研究中^[17-19]。2006年，Yasui等基于光学频率梳提出了太赫兹双梳光谱系统，并基于异步光学采样技术证明了太赫兹双梳光谱系统与太赫兹时域光谱系统之间的等价性^[20-21]。在过去的20年里，太赫兹双梳光谱系统取得了长足的进步，具有代表性的突破有：扩展时间延迟窗口，开发多个太赫兹脉冲的采集方法，获得了梳状模式分辨光谱，克服了重复频率限制的频率分辨率问题^[22]；为了不受泵浦和检测飞秒激光器的重复频率稳定和同步的限制，在太赫兹双梳光谱系统中采用了两个独立的自由运行飞秒激光器或一个具有两个输出端口的双向自由运行飞秒激光器^[23-24]；为了消除干涉图中有关激光定时抖动的失真，引入了参考通道来检测激光定时抖动，并开发了后/实时校正^[25]和自适应采样^[26-28]技术等。

在太赫兹双梳光谱系统中，包含测量信息的泵浦光频梳（频率间隔 $ ={f}_{\mathrm{r}} $ ）被另一个频率有微小偏移的探测光频梳（频率间隔 $ ={f}_{\mathrm{r}}+\Delta f $ ）采样，产生了射频梳齿（频率间隔 $ =\Delta f $ ）。由于射频梳本质上是泵浦光梳通过压缩因子（ $ {f}_{\mathrm{r}}/\Delta f $ ）降低频率后的线性变换，所以可以通过直接利用射频梳齿来分辨光谱。因此在泵浦光频梳重复频率确定的情况下，选择合适的重复频率差 $ \Delta f $ 决定了太赫兹双梳光谱系统的光谱质量。基于两台重复频率稳定的飞秒激光器搭建了一套太赫兹双梳光谱采样系统，通过微小改变两飞秒激光激光器的重复频率，来控制重复频率差 $ \Delta f $ 的变化，通过多频外差的方法直接使用频谱分析仪在不同的重复频率差下测量了系统的太赫兹频域信号，获得太赫兹频域谱。通过对比不同的重复频率差 $ \Delta f $ 得到的频域结果图，验证和分析了不同的重复频率差对太赫兹双梳光谱的影响，选取了合适的重复频率差，为后续太赫兹双梳光谱系统的重复频率差选定提供了方法。

1 实验原理及装置 1.1 光学频率梳原理

光学频率梳由光纤锁模激光器产生，在时域上的表现为等间隔的光学脉冲序列，在对飞秒光学频率梳的重复频率 $ f $ 和载波包络相位 $ \varphi $ 进行同步控制后，通过傅里叶变换，在频域上，光学频率梳则表现为离散且等间隔的频率梳齿线。

理想飞秒锁模激光器谐振腔输出的脉冲电场强度 $ \tilde{E}\left(t\right)$ 在时域中表达示为^[29]：

$ \tilde{E}(t)=A(t) \exp \left[{\rm{i}} \omega_{\mathrm{c}} t+{\rm{i}}\left(\frac{\varphi_{\rm {ceo }} t}{T_{\mathrm{r}}}\right)\right] \otimes \sum\nolimits_{m=-\infty}^{+\infty} \delta\left(t-m T_{\mathrm{r}}\right)$

(1)

式中： $ t $ 表示时间； $ A\left(t\right) $ 表示脉冲高斯包络函数； $ {\omega }_{\mathrm{c}} $ 表示载波频率； $ {\varphi }_{\rm{c}\rm{e}\rm{o}} $ 表示载波包络偏移； $ \delta \left(t- \right. \left. m{T}_{\mathrm{r}}\right) $ 表示狄拉克函数； $ m $ 为整数； $\otimes $ 表示卷积运算。在频域中可表示为等间隔的梳齿，即：

$ \tilde E\left(\omega \right)=\tilde{A}\left(\omega -{\omega }_{\mathrm{c}}\right){\sum }_{n=-\infty }^{+\infty }\delta \left(\omega -2 {\text{π}}n{f}_{\mathrm{r}}-2\mathrm{ {\text{π}} }{f}_{\mathrm{c}\mathrm{e}\mathrm{o}}\right) $

(2)

式中： $ n $ 为整数； $ \omega $ 表示角频率； $ {f}_{\mathrm{r}} $ 表示激光脉冲的重复频率； $ {f}_{\mathrm{c}\mathrm{e}\mathrm{o}} $ 表示载波包络偏移频率； $ \tilde{A}\left(\omega -{\omega }_{\mathrm{c}}\right) $ 表示频域电场包络。设第 $ n $ 个梳齿频率为 $ {f}_{n} $ ，若想让狄拉克函数不等于0，则 $ {f}_{n} $ 必须满足梳齿方程：

$ {f}_{n}=n{f}_{r}+{f}_{\mathrm{c}\mathrm{e}\mathrm{o}} $

(3)

式中： $ {f}_{\mathrm{c}\mathrm{e}\mathrm{o}} $ 表示载波包络偏移频率。式（3）也表示了通过光学频率梳可以连接光频波段（ $ {f}_{n} $ ）和射频波段（ $ {f}_{\mathrm{r}} $ ）。

1.2 太赫兹双梳光谱原理

太赫兹双梳光谱系统的测量原理与以光学频率梳作为光源的傅里叶变换光谱系统的测量原理十分相似^[30]，区别在于傅里叶变换光谱系统是将一束光分成两束，这两束光经过反射后合束，再经过样品被探测器采集，通过移动反射镜实现光谱的采集。而在太赫兹双梳光谱系统中，使用两个独立的脉冲激光器作为光源，泵浦光频率梳（重复频率 $ {=f}_{\mathrm{r}} $ ）和探测光频率梳（重复频率= $ {f}_{\mathrm{r}}+\Delta f $ ）之间存在一个微小的重复频率差 $ \Delta f $ ，重复频率差导致泵浦脉冲光和探测脉冲光之间的相对时间延迟增加，最大延迟时间是泵浦光频率梳的重复频率的倒数（ $ 1/{f}_{\mathrm{r}} $ ），因此在不需要机械延迟线的情况下泵浦脉冲激光激发的太赫兹时域信号可以被探测激脉冲激光所采集。理想上两频率梳的电场强度 $ E\left(t\right) $ 分别表示为：

$ {E}_{1}\left(t\right)=\sum\nolimits_{n}{A}_{1,n}\mathrm{exp}\left[-{\rm{i}}2 {\text{π}}\left\{{f}_{\mathrm{1,0}}+n\left(f+\Delta f\right)\right\}t\right] $

(4)

$ {E}_{2}\left(t\right)=\sum\nolimits _{n}{A}_{2,n}\mathrm{e}\mathrm{x}\mathrm{p}[-{\rm{i}}2 {\text{π}} \{{f}_{\mathrm{2,0}}+nf\left\}t\right] $

(5)

式中： $ {A}_{i,n} $ 和 $ {f}_{i,0} $ 分别表示脉冲包络的傅里叶分量和每个梳齿的载波包络偏移量频率(i = 1,2)。这两个脉冲序列在一个光电探测器上被检测，探测到的信号强度可以表示为：

$ I\left(t\right)\propto {\left|{E}_{1}\right(t)+{E}_{2}(t\left)\right|}^{2} $

(6)

这个信号包含许多频率分量，但在这里，只考虑频率梳1和梳2最近的梳齿之间的组合。在本实验系统中，使用低通电子滤波器来完成滤波。由式（4）、（5）和（6）可知，探测后的时域信号可表示为：

$ S\left(t\right)\propto \sum\nolimits _{n}{A}_{1,n}{A}_{2,n}^{*}\mathrm{e}\mathrm{x}\mathrm{p}\{-{\rm{i}}2 {\text{π}}(\Delta {f}_{0}+n\Delta f\left)t\right\} $

(7)

式中： $ \Delta f $ 表示两梳的重复频率差。时域经过傅里叶变换后的频域谱表示为：

$ B\left(\nu \right)\propto \sum\nolimits _{n}{A}_{1,n}{A}_{2,n}^{*}\delta \{\nu -(\Delta {f}_{0}+n\Delta f\left)\right\} $

(8)

式（8）表明待测的高频频率梳被下转到一个低频率区域的过程，其下转的压缩因子为 $ {f}_{\mathrm{r}}/\Delta f $ 。该压缩因子可以通过重复频率差的变化进行改变。

太赫兹双梳光谱的采样过程在时域上呈现为异步光学采样，产生相关干涉图，通过压缩因子（ $ {f}_{\mathrm{r}}/\Delta f $ ）的变换即可得到原太赫兹脉冲。异步光学采样技术即为使用两脉冲激光产生的电延迟线代替传统时域光谱中的机械延迟线实现采样。与使用机械延迟线的等效采样相比，异步光学采样技术具有采样速度快、误差小、可连续多信号连续采样和自由缩短或延长采样时间窗口等优点。在频域上两双梳光谱交差的区域产生外差干涉，在射频区域产生拍频信号，拍频后的梳齿间距为两光频梳的重复频率差 $ \Delta f $ ，探测光频梳对泵浦光频梳进行多频外差采样，通过两光频梳之间拍频作用将难以探测的光频段下换到电子器件可探测的射频段，通过对采集到射频段的电信号乘以压缩因子（ $ {f}_{\mathrm{r}}/\Delta f $ ）的计算，即可完成对太赫兹频段信号的采样还原。本文中的太赫兹双梳光谱系统与传统的太赫兹时域光谱系统之间的优劣势如表1所示。

太赫兹双梳光谱的采样原理如图1所示。图1（a）展示了异步光学采样的原理（太赫兹双梳光谱时域信号采集原理）^[31]，重复频率为 ${f}_{{\rm{r}}}$ 的泵浦脉冲激光照射进发射光电导天线的电极间隙处激发出电子空穴对，此时给天线外加偏压后外加电场使电子空穴对加速运动，并辐射出重复频率为 $ {f}_{r} $ 太赫兹脉冲（蓝）。重复频率为 $ {f}_{\mathrm{r}}+\Delta f $ 探测脉冲激光（红）照射进探测天线中，激发出电子空穴对，这些电子空穴对在太赫兹信号产生的电场作用下进行加速，产生电信号被探测仪器接收，对接收到的电信号进行处理后得到原始的太赫兹信号从而完成对太赫兹时域脉冲信号的采集。

设在实验室时间 $ t=0 $ 时，发射光路的太赫兹脉冲和探测光路的飞秒脉冲同时到达探测器，那么在 $ t=1/{f}_{\mathrm{r}} $ 时，太赫兹脉冲与探测脉冲之间存在时间间隔 $ \tau $ 为：

$ \tau =\frac{1}{{f}_{\mathrm{r}}}-\frac{1}{{f}_{\mathrm{r}}+\Delta f}=\frac{\Delta f}{(\Delta f+{f}_{\mathrm{r}}){f}_{\mathrm{r}}}\approx \frac{\Delta f}{{{f}_{\mathrm{r}}}^{2}} $

(9)

采集一个完整太赫兹脉冲所用时间为：

$ t=\frac{1}{{f}_{\mathrm{r}}+\Delta f}\times \frac{{f}_{\mathrm{r}}}{\Delta f}\approx \frac{1}{\Delta f} $

(10)

所以当实验室时间 $ t=1/\Delta f $ 时，一个完整太赫兹脉冲被采样完成，且总采样时间 $ \tau =1/{f}_{\mathrm{r}} $ 。实验室时间 $ t $ 和采样时间 $ \tau $ 之间的关系表达式为：

$ t=\tau \times \frac{{f}_{\mathrm{r}}}{\Delta f} $

(11)

当太赫兹脉冲被完全采样后，即可对时域信号进行傅里叶变换从而获得相对应的信号频域谱。

如图1（b）展示了多频外差采样的原理（太赫兹双梳光谱频域信号采集原理）。当重复频率为 $ {f}_{r} $ 的泵浦光频梳（红）入射到太赫兹发射光电导天线时，给天线外加偏压后辐射出太赫兹梳，由于光电导天线发射器通过二次检测，使光梳的电场平方，将光强度转换为载流子，并通过偶极子辐射发射电磁太赫兹波的频率梳，因此它向下转换至太赫兹频段，而不改变梳齿频率间距。泵浦太赫兹梳（红）是一个没有频率偏移的谐波频率梳，其频率 $={f}_{{\rm{r}}}$ , 2 ${f}_{{\rm{r}}}$ , 3 ${f}_{\mathrm{r}},\;\cdots, \;n{f}_{\mathrm{r}}$ 。重复频率为 $ {{f}_{\mathrm{r}2}=f}_{\mathrm{r}}+\Delta f $ 的探测光频梳（蓝）照射进探测天线中，探测光频梳的频率 $ ={f}_{\mathrm{r}2} $ , 2 $ {f}_{\mathrm{r}2} $ , 3 ${f}_{\mathrm{r}2} , \; \cdots ,\; n{f}_{\mathrm{r}2}= {f}_{\mathrm{r}}+\Delta f$ , 2 $ {(f}_{\mathrm{r}}+\Delta f) $ , 3 ${(f}_{\mathrm{r}}+\Delta f),\;\cdots ,\; n({f}_{\mathrm{r}}+\Delta f)$ ，此时在光电导探测器中探测梳与泵浦梳发生多频外差光导混合，产生一个重复频率为 $ \Delta f $ 的射频梳，使用频谱分析仪即可直接对射频梳进行采集。对采集到的射频梳进行压缩因子 $ {f}_{\mathrm{r}}/\Delta f $ 的转换即可回转得到太赫兹梳的信息，完成对太赫兹脉冲的频域信号采集^[32]。

异步光学采样和多频外差采样的方法都将高频的太赫兹信号下转到低频的射频信号后进行采集。时域测量完成后需要完成傅里叶变换后得到频域信息，用频谱仪可以直接在频域中进行全带宽的测量和采集，同时也可对某一特定频段的信号进行直接采集和分析，这使太赫兹信号的采集和分析变得简易且方便。该采样方法在分辨梳齿线的精确测量和超快速测量方面都具有巨大的潜力。

1.3 实验装置

为了验证测量重复频率差对太赫兹双梳光谱的影响，设计了图2（a）的光路原理图，并根据光路原理图搭建了图2（b）实验系统光路图。实验开始前，我们需要开启两台飞秒脉冲激光器（Menlo, C-Fiber 780）预热90 min后打开稳频同步系统（Menlo Systems, ASOPS, EM100），设置两台激光器的的重复频率后开启同步控制系统。激光器1的重复频率为100 MHz，激光器2的重复频率为100 MHz– $\Delta f$ ，两激光器的输出波长均为780 nm，平均功率均大于100 mW，脉冲宽度均小于100 fs。为了防止激光的光功率过高对天线造成损伤，在光路系统中加入光学衰减片来控制光功率，使照射进天线的光功率符合天线输入功率标准防，止损伤天线电极。激光器1射出的泵浦脉冲激光（红）经过反射镜调整光路后，照射进发射光电导天线（Teravil，ETM-8）中，通过电源表（KEITHLEY，2400）对发射天线外加30 V电压后辐射出太赫兹脉冲。另一激光器2射出的探测脉冲激光（蓝）经过反射镜后射进探测天线（Batop，PCA-44-06-10-800-h）中，将探测天线输出的电信号经过前置放大器（FEMTO，DHPCA-100 I/V）的放大（放大增益为10⁶ V/A）后接入探测仪器中，可以实现对太赫兹脉冲信号的探测。为了探测时域信号，将放大后的电信号接入示波器（KEYSIGHT，DSAZ504A，采样率为1 G sample/s）中，通过调整示波器的时间轴标度，来放大和缩小延迟时间范围，从而可以探测一个或者多个周期的太赫兹脉冲信号，经过傅里叶变换后，得到信号的频谱信息。

由于太赫兹双梳光谱具有同时探测时域和频域信号的优点，也可以将探放大后的电信号接入频谱信号分析仪（KEYSIGHT，N9040B）中来直接探测太赫兹频域信号。通过调整频谱信号分析仪的分辨率带宽、扫宽等参数即可直接获得射频波段的频谱数据，再将射频数据经过压缩因子 $ {f}_{\mathrm{r}}/\Delta f $ 的变换后得到太赫兹段的频谱数据从而完成太赫兹信号的频域测量。

在太赫兹双梳光谱系统中，希望所有的数据采样值都是真实有效的，所以此系统需要一个合适的重复频率差值。假设需要测试的频谱宽度为∆F，由奈奎斯特-香农采样定理可知，光学采样的自由光谱范围为0到 $ {f}_{\mathrm{r}}/2 $ 。由于下转换光谱的带宽为 $\Delta f \times \Delta F/{f}_{\mathrm{r}}$ ，由此关系式可知改变重复频率差会影响光谱跨度范围。而如果被下转换的太赫兹双梳光谱不在这个自由谱的范围内，采样过程就会产生混叠效应。为了避免出现这种混叠影响信号质量，重复频率差需要满足以下条件：

$ \mathrm{\Delta }F\times \frac{\Delta f}{f}\leqslant \frac{f}{2} $

(12)

$ \Delta f\leqslant \frac{{f}^{2}}{2\mathrm{\Delta }F} $

(13)

因此需要式（13）所示的范围内选择合适的重复频率差值。由压缩因子 $ {f}_{\mathrm{r}}/\Delta f $ 可知，改变两光频梳的重复频率差 $ \Delta f $ 会直接影响低频到高频的变换系数从而影响到光谱的质量，同时由式（11）可知改变 $ \Delta f $ 也会影响太赫兹脉冲的采集的速度。由式（10）可知两光梳间的重复频率差越小，光采样时间间隔越小，得到的频谱就越精细，因此理论上减小两光频梳的重复频率差有利于测量精度的提升。

2 结果与分析讨论

以下所有测量均在环境温度为25 ℃，相对湿度为50%条件下进行。首先开启两台飞秒激光器预热90 min后，开启重复频率稳定和同步控制系统，将激光器1的重复频率设为100 MHz，激光器2的重复频率分别设置为100 MHz～50 Hz、100 MHz～40 Hz、100 MHz～25 Hz、100 MHz～10 Hz、100 MHz～5 Hz、100 MHz～2 Hz。每次改变激光器2的重复频率后，待稳定控制系统90 min后测量。频域信号的结果如图3（a）～（f）所示。

由于频谱仪直接测得的信号位于kHz的射频波段，需要通过压缩因子 $ {f}_{\mathrm{r}}/\Delta f $ 的变换上转到太赫兹波段，得到原始的太赫兹频谱。如图3所示，下横轴坐标为测试所得射频段，上横轴坐标为上转至所在的太赫兹波段。由6次测量的上横坐标可知，本次实验所测高频太赫兹信号波段为0～1.5 THz，而当 $ \Delta f $ 改变时，所需测量的射频波段的范围也会随之发生改变。例如当 $ \Delta f $ = 50 Hz时，此时的压缩因子 $ {f}_{\mathrm{r}}/\Delta f= $ 100 MHz/50 Hz = 2 $ \text{×} $ 10⁶，那么高频0～1.5 THz范围下转到射频段范围为0～750 kHz，所以当两光梳重复频率差为50 Hz时，实验所需测量射频波段为0～750 kHz。同理，当 $ \Delta f $ 分别为40 Hz、25 Hz、10 Hz、5 Hz和2 Hz时，频谱仪所需测得的射频范围分别为0～600 kHz、0～375 kHz、0～150 kHz、0～75 kHz、0～30 kHz。由计算和测量结果可知测量同一波段的高频信号时，两光梳间的重复频率差越小，需要测得射频波段范围越小。

为了验证理论计算与实际测量的准确性，在6次测量的每个结果中都选取了1.097 THz和1.163 THz位置的水汽吸收峰^[33]，分别标识为水汽吸收峰1和水汽吸收峰2，将6个测试结果图叠加后如图4所示。由图4的吸收峰放大叠加图可以清晰地看出，6次测量的频谱结果图在1.097 THz和1.163 THz的水汽吸收峰位置严格对应，没有偏移。证实了此次测量的频谱结果的正确性和一致性。由图3所示结果可知，尽管都是对原太赫兹频谱信号的还原，但是由于重复频率差 $ \Delta f $ 的不同，压缩因子 $ {f}_{\mathrm{r}}/\Delta f $ 的值会变化，导致实际所测得射频段范围也不同，进而对原信号的还原效果也有显著差异。所以选择合适的重复频率差，也是影响太赫兹双梳光谱测量质量的重要因素。

当重复频率差 $ \Delta f $ 大于25 Hz时，结果显示在射频段150 kHz位置底噪有明显升高。这是由于信号在输入频谱分析仪前接入了直流低通滤波器，导致低频段的噪声通过，高频段部分被滤除。当重复频率差 $ \Delta f $ 变大时，压缩因子 $ {f}_{\mathrm{r}}/\Delta f $ 的值变小，则需要测量的射频波段增长，高频段引入的偏移噪声增加，因此影响了高频信号的频谱质量，在结果图中显示为1.3 THz后信号梳齿显著降低。但是当重复频率差 $ \Delta f $ 过小时，梳齿会出现断层，采样仪器无法采集到完整的梳齿以及底噪，造成频谱信息的缺失。为了直观对比结果，选取了0.4888～0.490 THz高频位置重复频率差为10 Hz和2 Hz的梳齿对比如图5所示。

图5（a）为重复频率差 $ \Delta f $ = 10 Hz时太赫兹梳齿的放大图，由频谱结果图看出梳齿完整稳定，底噪位置清晰确定。图5（b）为重复频率差 $ \Delta f $ = 2 Hz时太赫兹梳齿的放大图，可以看出梳齿清晰但不稳定，出现上下浮动的现象，梳齿不完整且无法观测到信号底噪，严重影响了频谱的质量。因此相比于其他结果，重复频率差 $ \Delta f $ = 5 Hz和10 Hz时，输出频谱质量最好。而相比于 $ \Delta f $ = 5 Hz， $ \Delta f $ = 10 Hz的梳齿信号的动态范围更大，约为47 dB（ $ \Delta f $ = 5Hz时约为40 dB），所以选取 $ \Delta f $ = 10 Hz为最佳重复频率差。因此可知并非重复频率差 $ \Delta f $ 越小越好，当重复频率差过小时，频谱信息失真，影响信号的质量。因此选取合适的重复频率差可获得最佳的频谱信号。

3 结　论

本文从光学频率梳和太赫兹双梳光谱的原理出发，基于可稳定重复频率的两台飞秒激光器和两个光电导天线搭建了一套完整的太赫兹双梳光谱采样系统。通过稳定泵浦激光器的重复频率并轻微改变探测激光器的重复频率来控制两激光器的重复频率差。在不同的重复频率差下测量了系统的太赫兹频域信号，通过对比实验结果分析了改变重复频率差对太赫兹双梳光谱的影响。结果表明，重复频率差变小时，压缩因子随之变大，光谱系统测得频谱越精细，频谱质量越好。但当重复频率差 $ \Delta f $ 低于10 Hz时，光谱质量下降，频谱信息失真，影响测试结果。因此在环境温度为25 ℃、相对湿度为50 %和测量太赫兹频谱范围在0～1.5 THz等测量条件的情况下选取 $ \Delta f $ = 10 Hz作为最佳重复频率差条件。由此，通过原理分析和实验结合的方法，选取了最佳的重复频率差测试条件，为基于飞秒激光器的太赫兹双梳光谱提供了选取最佳重复频率差的方法。