一些关键的时间和频率测量概念（抖动，漂移，MTIE和TDEV）
1 与时频相关的技术指标说明
1.1 概述
在卫星导航中，我们常遇到秒脉冲的概念，其精度以峰峰值方式表示，一般峰峰值可以理解为秒脉冲误差分布的 3
σ
\sigma
σ。这里的3
σ
\sigma
σ是样本标准差，是方差算术平方根，克服了误差均值
μ
\mu
μ的影响。因此，可以看到，精度是抖动的表征，而准确度也就是偏差，则用误差的均值
μ
\mu
μ来表征。
在时频传递系统中，一方面我们关心重建时间的特征，主要是指重建的1PPS的上升沿，其与原始参考时间（可能是某一个原子钟生成的秒脉冲，也可能是一个UTC整秒）之间的对齐程度称作零值偏差。对于宽带卫星通信系统，与一般卫星导航系统一样，其核心网时间是一个“纸面时间”，由多个信关站的原子钟联合维护，该时间与UTC时间的差需要一个专门的校准系统进行持续维护，其与UTC时钟存在一个零值偏差。另一方面，我们关心驯服后的频率的特征。重建的秒脉冲本身也是一个特殊的1Hz时钟，除了1PPS外，一般通信系统都还另有至少一个频率时钟通常不低于10MHz的重建时钟。显然，重建的秒脉冲和时钟与原始秒脉冲和时钟之间是存在频率和相位误差的，如何评价这些误差呢？
1.2 频率指标
1.2.1 时钟频率
通常所说的原子钟频率数据，实际上指的是瞬时相对频率偏差数据，它的物理意义是相位变化对时间的微分，其表达式为
y
i
=
x
i
+
1
−
x
i
τ
0
{y_{_i}} = \frac{{{x_{i + 1}} - {x_i}}}{{{\tau _0}}}
yi=τ0xi+1−xi
式中
x
i
+
1
{x_{i + 1}}
xi+1和
x
i
{x_i}
xi分别表示第
i
+
1
i+1
i+1和第
i
i
i历元的钟差（相位）值，
τ
0
{\tau _0}
τ0表示两个历元间的采样间隔。
####1.2.2 频率准确度
准确度是测量值或计算值与其定义的符合程度，其表达式为：
A
=
f
−
f
0
f
0
A = \frac{{f - {f_0}}}{{{f_0}}}
A=f0f−f0
式中
f
f
f为被测信号的实际频率，
f
0
f_0
f0 为测量信号的标称频率。实际实验中一般都要进行多次测量取平均求得对应的准确度指标，这是因为受测频标的实际输出频率并不是固定不变的。此外，现实情况下实际频率和标称频率之间的偏差无法直接测量，一般把参考信号的实际频率当做标准进行测量，且参考频率准确度应该比被测频标的准确度高出至少一个数量级。
1.2.3 频率漂移率
原子钟在连续运行的过程中，由于硬件设备的老化和周围环境变化等因素的影响，导致其输出的频率值经常会随运行时间呈现出线性递增或递减的现象。这种频率随时间单调变化的速率称为频率漂移率，也叫频率老化率，其表达式为：
z
i
=
y
i
+
1
−
y
i
τ
0
{z_{_i}} = \frac{{{y_{i + 1}} - {y_i}}}{{{\tau _0}}}
zi=τ0yi+1−yi
式中
y
i
+
1
{y_{i + 1}}
yi+1 和
y
i
{y_{i}}
yi分别表示第
i
+
1
i+1
i+1和第
i
i
i历元的频率值，
τ
0
{\tau _0}
τ0为相邻历元间的采样间隔。实际应用中，通常取多个历元瞬时频率漂移率的平均值，或利用最小二乘法获取最优频率漂移率。
1.2.4 频率稳定度
频率稳定度是表征时钟在一定时间内产生同样频率的能力，评估稳定度主要采用的指标是Allan方差，其表达式为：
σ
y
2
(
τ
)
=
⟨
(
y
i
+
1
−
y
i
)
2
2
⟩
=
1
n
∑
i
=
1
n
(
y
i
+
1
−
y
i
)
2
2
=
1
2
n
τ
0
2
∑
i
=
1
n
(
x
i
+
2
−
2
x
i
+
1
+
x
i
)
2
\sigma _y^2(\tau ) = \left\langle {\frac{{{{({y_{i + 1}} - {y_i})}^2}}}{2}} \right\rangle = \frac{1}{n}\sum\limits_{i = 1}^n {\frac{{{{({y_{i + 1}} - {y_i})}^2}}}{2}} = \frac{1}{{2n{\tau _0}^2}}\sum\limits_{i = 1}^n {{{({x_{i + 2}} - 2{x_{i + 1}} + {x_i})}^2}}
σy2(τ)=⟨2(yi+1−yi)2⟩=n1i=1∑n2(yi+1−yi)2=2nτ021i=1∑n(xi+2−2xi+1+xi)2
式中
n
n
n为采样点的个数，
x
i
x_i
xi表示第
i
i
i历元的相位，
y
i
y_i
yi表示第i历元的频率值， 为相邻历元间的采样间隔。随时间间隔
n
τ
0
n{\tau _0}
nτ0的取值从秒到数天的跨度，分别对应短期、中期和长期稳定度。下图直观的反映了频率准确度和稳定度之间的关系。
图1-1 频率稳定度和准确度示例
1.3 时钟抖动和游走指标
1.3.1 时钟抖动
生成的时钟相对于原始参考时钟，会存在抖动，如下图。
图1-2 时钟抖动
根据ITU-T G.810，低于10Hz的时钟抖动（Jitter）称作游走（Wander）。
时钟的抖动常用的简单方法通过峰峰值表征，比如GNSS接收机输出秒脉冲的抖动是±50ns，有时候也会用RMS来表示，这里的RMS是抖动对应的“高斯”分布的1σ。显然，这里认为的抖动对应的长期零偏为0的。
显然，这样的表征不能有效的表示出“游走”的特征。为此，《IEEE 2414-2020 Standard Jitter and Phase Noise》专门定义了TIE、MTIE、TDEV、ADEV等技术参数，来表征“游走”， MTIE和TDEV通常使用的较多。
1.3.2 TIE、PPTIE、MTIE
时间间隔误差（Time Interval Error，TIE）是指在特定时间点上，一个给定信号相对于参考时钟的相位差，测量的间隔为
τ
0
{\tau _0}
τ0 ，下图中 。TIE的测量周期为
τ
0
=
1
s
{\tau _0}=1s
τ0=1s ，其起始时刻的初值通常设置为0，如下图的例子，
T
=
100
s
，
N
=
100
T=100s，N=100
T=100s，N=100。
图1-3 TIE游走测量示意
下一篇：
彻底理解晶振的工作原理和参数
显然，给定一个窗口
τ
=
n
τ
0
\tau {\rm{ = n}}{\tau _0}
τ=nτ0 ，上图中为20s，都能够得到一个最大的TIE值，称作PPTIE，而最大时间间隔误差（Maximum Time Interval Error，MTIE）是指在一个测量周期
T
=
100
s
T=100s
T=100s内，所有给定滑动窗口 内的最大相位变化，也就是这段时间的PPTIE的最大值。其表达式和示意图如下：
M
T
I
E
(
n
τ
)
=
max
1
≤
k
≤
N
−
n
+
1
(
max
k
≤
i
≤
k
+
n
−
1
(
x
i
)
−
min
k
≤
i
≤
k
+
n
−
1
(
x
i
)
)
MTIE(n\tau ) = \mathop {\max }\limits_{1 \le k \le N - n + 1} (\mathop {\max }\limits_{k \le i \le k + n - 1} ({x_i}) - \mathop {\min }\limits_{k \le i \le k + n - 1} ({x_i}))
MTIE(nτ)=