為什麼不能用標準差
假設你手上有一台 10 MHz 頻率標準，接上計數器連續量了一串讀值，想描述它「穩不穩」。直覺的做法是算平均值與標準差——這在量一根鋼尺的長度時完全正確，但用在振盪器上會出事，而且是根本性的出事。
原因在於一般的標準差（樣本變異數）背後有一個假設：待測的隨機過程是平穩的、變異數有限的。對白雜訊而言這成立，取樣愈多，估計愈準，數字會收斂到一個穩定值。但振盪器的雜訊不只有白雜訊，它還帶著閃爍雜訊 (flicker noise，功率密度正比於 1/f) 與隨機漫步 (random walk)，這兩種過程的能量在低頻端是發散的——頻率愈低、能量愈大，而且沒有下界。對這類過程算標準差，結果不會收斂：你量得愈久，讀值裡就混進愈低頻的成分，標準差就愈大。
把數字放進來會更有感覺。同一顆 OCXO，量一分鐘可能得到 2×10⁻¹¹，量一小時得到 5×10⁻¹⁰，量一整天得到 1×10⁻⁸——振盪器並沒有在這段時間裡變差，變的只是你把多少低頻漂移收進了統計裡。這代表「標準差」這個數字同時取決於你量了多久、取了幾點，換一組量測條件就換一個答案。一個會隨量測時間長度改變的數字，不能拿來當規格，也不可能讓兩家廠商公平比較。
1966 年 David W. Allan 提出了雙樣本變異數 (two-sample variance) 來解決這件事，後來成為 IEEE 與 ITU 的標準時域穩定度指標，也就是今天的亞倫變異數 (Allan variance) 與其平方根亞倫偏差 (Allan deviation)。它的巧思其實只有一句話：不要拿每個樣本去和「全體平均」比——因為全體平均本身就在漂——而是拿每個樣本去和它的鄰居比。
定義：相鄰兩段平均的差
先把量測值正規化。振盪器的標稱頻率是 f₀，某一瞬間的實際頻率是 f，定義分數頻率偏差 y = (f − f₀)/f₀。它是無因次的，因為除掉了標稱頻率，所以 10 MHz 標準與 100 MHz 標準的穩定度可以直接放在同一張圖上比較。把 y 在一段長度為 τ 的時間內取平均，得到第 k 段的平均值 ȳₖ；τ 就是平均時間 (averaging time)，也是整件事的自變數。
亞倫變異數的定義是相鄰兩段平均值之差的均方值的一半：σy²(τ) = ½·⟨(ȳₖ₊₁ − ȳₖ)²⟩，開根號就是亞倫偏差 σy(τ)。前面那個 ½ 不是隨便放的——它讓亞倫變異數在純白色頻率雜訊的情況下剛好等於古典變異數，使兩者在最單純的情況下數值一致。因為 y 無因次，σy(τ) 也無因次，所以規格會寫成 5×10⁻¹²、3×10⁻¹¹ 這樣的形式，而且一定要附上 τ：τ = 1 s 與 τ = 100 s 的值可以差一個數量級以上，只給數字不給 τ 的穩定度規格，和只給 dBc/Hz 不給偏移頻率的相位雜訊規格一樣沒有意義。
「相鄰差分」為什麼救得回收斂性？因為對相鄰兩段取差，等效於對訊號做了一次高通濾波：常數的頻率偏移被完全消掉，非常低頻的成分也被大幅壓抑，於是那些在低頻發散的雜訊過程不再把積分撐爆，估計值就能收斂。這裡要順帶澄清一個常見的誤解：亞倫偏差描述的是穩定度，不是準確度。一台頻率偏了 1 Hz 的 10 MHz 標準，只要它偏得很穩定，ADEV 可以非常漂亮——穩定度講的是「會不會變」，準確度講的是「離真值多遠」，是兩個獨立的問題。
實際怎麼量？你需要一個至少和待測物一樣好的參考（或者用三角帽法 (three-cornered hat) 由三台互比拆解出各自的貢獻）。做法通常是用計數器連續記錄待測物相對於參考的相位隨時間變化，再由這串資料算出各個 τ 的 ADEV。這裡有個容易被忽略的陷阱：相鄰兩次量測之間如果有死區時間 (dead time)，被漏掉的那段訊號會讓結果產生偏差，尤其在短 τ 特別明顯。這正是連續零死區時間時間戳記重要的原因——Pendulum CNT-91/91R 的連續時間戳記、CNT-104S 與 CNT-104R 的無間隙 (gap-free) 量測就是為此設計的，這些機種也直接內建 ADEV 統計。計數器本身的解析度與觸發誤差同樣會落在短 τ 的曲線上，相關細節見〈頻率計數器入門：閘門時間與解析度〉(/knowledge/frequency-counter-basics)。
怎麼讀一張 ADEV 圖：斜率就是雜訊的身分證
標準的 ADEV 圖橫軸是平均時間 τ、縱軸是 σy(τ)，兩軸都是對數刻度。橫軸從左往右是「平均得愈久」，縱軸愈往下愈穩定。整條曲線通常由幾段近似直線接成，而每一段的斜率，直接對應到一種物理上的雜訊機制——這是這張圖最有價值的地方：它不只告訴你有多穩，還告訴你不穩的原因是什麼（見圖 1）。
由左至右依序是：斜率 −1，來自白色相位雜訊與閃爍相位雜訊，通常是量測系統本身、緩衝級與放大器的貢獻，拉長平均時間就會等比例改善；斜率 −1/2，是白色頻率雜訊，多半來自被動諧振器的熱雜訊，平均 100 倍的時間可以改善 10 倍；斜率 0 的水平段，是閃爍頻率雜訊，來自諧振器本身，這條平底就是這顆振盪器的穩定度極限，再怎麼平均都不會更好；斜率 +1/2，是隨機漫步頻率雜訊，通常對應環境因素——溫度、電源、機械應力、負載牽引；最後是斜率 +1，那是線性頻率漂移，也就是老化，對漂移率 D 而言 σy(τ) ≈ D·τ/√2。
實務上真正要看的只有三件事。第一，平底在哪個高度——那是這顆振盪器的極限。第二，平底從哪個 τ 開始、到哪個 τ 結束——這決定了你有沒有辦法靠拉長平均時間換到你要的穩定度。第三，你的應用實際平均多久——如果你的量測就是平均 100 秒，那麼 τ = 100 s 那一點以外的曲線再漂亮也與你無關。曲線上的隆起同樣會說話：伺服迴路或鎖相迴路的時間常數常會在對應的 τ 留下一個包；恆溫槽的溫度循環會在循環週期附近留下一個包；GNSS 鎖定的參考源則會在鎖定迴路時間常數附近出現轉折。
最後是一個誠實但常被略過的限制：長 τ 的資料點很不可靠。要估 τ = 1 天的 ADEV，一段三天的量測裡只有寥寥幾個獨立樣本，信賴區間大得驚人。所以看到一條從一天份資料一路畫到 τ = 10⁵ 秒的曲線，合理的反應是懷疑而不是相信。認真的報告會附上誤差棒，也會標明總量測時間；比較兩份 ADEV 資料時，先問對方量了多久，再問 τ。
變形：重疊式 ADEV、MDEV 與 TDEV
最常見也最該預設開啟的變形是重疊式亞倫偏差 (overlapping ADEV)。原始定義是把資料切成互不重疊的區塊，一段 τ 接一段 τ；重疊式則是讓每一個可能的起點都拿來算一次。兩者的期望值相同，但重疊式用同一份資料得到多得多的樣本數，信賴區間明顯收窄，長 τ 的資料點也穩定得多。現代分析軟體多半預設就是重疊式，報告上也應該標明——除非你有特別的理由，重疊式幾乎總是正確的選擇。
第二個變形是修正型亞倫偏差 (modified Allan deviation, MDEV)。前面提過，白色相位雜訊與閃爍相位雜訊在 ADEV 上都呈現斜率 −1，你分不出是哪一種。MDEV 在計算前多做一層相位平均，效果是把這兩者的斜率拉開（白色相位雜訊變成 −3/2，閃爍相位雜訊維持 −1），於是短 τ 區的雜訊型態終於可以辨識。什麼時候需要它？當你關心的正是短 τ 那一段——例如在鑑定分配放大器、緩衝級，或想確認曲線的短 τ 段到底是待測物還是量測系統造成的時候。
第三個是時間偏差 (time deviation, TDEV, σx(τ))，由 MDEV 換算而來：σx(τ) = τ·MDEV(τ)/√3。注意它的單位是秒，不是無因次量——這正是重點。在通訊網路裡，工程師關心的不是分數頻率而是時間誤差有幾奈秒，所以 ITU-T 的同步標準是以 TDEV 與 MTIE 來規範的。Pendulum 的 TimeView™ 3 調變域分析軟體就直接提供 MTIE 與 TDEV 這兩個漂移分析指標。
MTIE 與 TDEV 常被混用，但它們量的是不同的東西。TDEV 是均方根量，描述的是雜訊的統計特性；MTIE (maximum time interval error) 是峰值量，描述的是在任意觀測窗內出現過的最大時間誤差，因此對單一次的相位跳變、切換事件特別敏感。一條鏈路完全可能 TDEV 過關卻 MTIE 不合格——因為統計上很乾淨，但發生過一次跳。要用哪一個，取決於你被稽核的是哪一份標準，兩個都列上通常是最安全的做法。
和相位雜訊是同一件事：L(f) 與 σy(τ)
亞倫偏差與相位雜訊常被當成兩個獨立的規格看待，但它們描述的是同一顆振盪器的同一組雜訊，差別只在座標系：相位雜訊 L(f) 是頻域的描述，橫軸是離載波的偏移頻率；亞倫偏差 σy(τ) 是時域的描述，橫軸是平均時間。兩者可以互相換算，關係式是 σy²(τ) = 2∫Sy(f)·sin⁴(πτf)/(πτf)² df，其中 Sy(f) = (f²/ν₀²)·2L(f)（見圖 2）。那個 sin⁴ 的核就是前面說的「相鄰差分」在頻域裡的長相：一個把直流與極低頻擋掉的濾波器。
換算表本身很值得記起來，而且有一個容易搞混的地方：兩張圖的方向是相反的。頻域的遠偏移對應時域的短平均時間，頻域的近載波對應時域的長平均時間。斜率的對應是：白色相位雜訊 f 的 0 次方對應 τ 的 −1 次方；閃爍相位雜訊 f⁻¹ 對應約 τ⁻¹；白色頻率雜訊 f⁻² 對應斜率 −1/2；閃爍頻率雜訊 f⁻³ 對應水平段；隨機漫步頻率雜訊 f⁻⁴ 對應斜率 +1/2。所以當有人說「這顆的閃爍底很高」，他在兩張圖上指的是同一件事。
既然是同一件事，該用哪一個？答案由時間尺度決定。毫秒以下的行為、以及 RF 系統關心的效應——交互混頻抬高接收機雜訊底、高階 QAM 的 EVM 底限、取樣時脈的抖動——都落在偏移頻率 1 Hz 以上的區域，用相位雜訊描述又自然又量得準（見〈相位雜訊：載波旁邊的那圈裙擺〉/knowledge/phase-noise）。秒級到天級的行為則相反：要用相位雜訊描述 τ = 1 天，等於要量到 10⁻⁵ Hz 的偏移，實務上做不到，這時亞倫偏差才是對的工具。兩者是互補的，不是競爭的。
換算也有它的限制，值得說清楚。那個積分只對某些雜訊型態收斂；線性漂移（老化）在 L(f) 裡根本沒有對應項，所以由相位雜訊換算出來的 ADEV 不會包含漂移，長 τ 那一段會過於樂觀。此外，換算結果會繼承原始量測的所有限制，尤其是量測頻寬的上下限——積分帶不同，換算出的 σy(τ) 就不同。實務上，專用相位雜訊分析儀多半可由相位雜訊資料直接運算亞倫變異數，而計數器則是直接量時域；如果兩邊的數字對不起來，第一個該檢查的就是積分頻寬與量測時間，而不是急著懷疑硬體。
從曲線到採購決定：OCXO、銣鐘、GNSS 鎖定與保持
把前面所有東西收斂成一個可執行的流程，只有兩步：先從應用釘死你關心的 τ，再拿這個 τ 去比規格。τ 從哪裡來？產線上用 1 秒閘門做頻率測試，τ 就是 1 秒；一次衛星通過要維持同調 100 秒，τ 就是 100 秒；校正實驗室以 24 小時平均出具報告，τ 就是 86400 秒；系統要求失去 GNSS 後 8 小時內時間誤差不超過某個值，那要看的是保持規格。τ 沒釘死之前，任何「哪一顆比較好」的討論都不會有答案。
釘死 τ 之後，三種參考源的排序會出現一個讓很多人意外的結果（見圖 3）。以 Pendulum 6688/6689 的實際規格為例：採用 SC 切割 OCXO 的 6688，短期穩定度在 τ = 1 s 與 τ = 10 s 都是 5×10⁻¹²；採用銣原子鐘的 6689，同樣條件下是 3×10⁻¹¹ 與 1×10⁻¹¹。也就是說，在 τ = 1 秒，好的 OCXO 比銣鐘乾淨了將近一個數量級。短期穩定度是晶體的問題，不是原子的問題——原子鐘的價值在中長期：6688 的老化率是 3×10⁻⁹/月、2×10⁻⁸/年，而 6689 是 5×10⁻¹¹/月、10 年累計不超過 1×10⁻⁹。兩條曲線大約在數十秒的 τ 交叉，交叉點之後銣鐘才開始勝出。
再往長 τ 走，老化就成為兩者共同的敵人，這時 GNSS 鎖定登場。GNSS 鎖定的參考源以衛星上的銫原子鐘為長期基準，把本機振盪器持續拉回來，因此長期頻率準確度不再隨時間老化。實際數字：GPS-12R（GPS 鎖定銣鐘）在鎖定狀態下 τ = 1 s 優於 3×10⁻¹¹、τ = 100 s 優於 5×10⁻¹²、τ = 24 h 優於 2×10⁻¹²；GPS-88（本機為 OCXO）鎖定 GPS 時 24 小時平均頻率偏移優於 2×10⁻¹²，並持續比對本機與 GPS 訊號、把偏差存進非揮發性記憶體，可隨時以 GPSView™ 列印可溯源校正報告。注意 GNSS 對短 τ 幾乎沒有幫助——短 τ 仍然由本機振盪器決定，這也是為什麼 GNSS 鎖定機種還要分本機是 OCXO 還是銣鐘。
最後是保持 (holdover)，也就是失去 GNSS 之後的表現，它是另一個獨立的數字，不能用鎖定時的 ADEV 代替。FTR-210R 這類 GNSS 鎖定銣原子參考源標示的是 1 pps 對 UTC 優於 10 ns rms、保持模式漂移典型 1 µs/24 小時、手動保持模式老化率優於 5×10⁻¹¹/月。看保持規格要問三件事：撐多久、允許累積多少時間誤差、以及在什麼溫度條件下——因為保持期間的主導誤差項通常是溫度係數而不是老化，一份沒有標明溫度條件的保持規格幾乎沒有參考價值。對台灣的校正實驗室還有一個現實考量：GNSS 鎖定加上獨立的文件化頻率比較器（例如 FTR-210R 的 Option 220），意味著可以隨時自行產出可溯源報告，而不必每年把標準送出去校；若同時需要多通道分析能力，CNT-104R 把銣鐘與四通道無間隙分析儀合在一台，搭配 GNSS 選項可達 24 小時平均 1×10⁻¹²、時間校準不確定度優於 10 ns rms。
