在通訊和PC行業(yè)，高速串行信號(hào)越來(lái)越普及，在使用示波器測(cè)量和分析這類信號(hào)時(shí)，通常要求測(cè)量總體抖動(dòng)（Total jitter，簡(jiǎn)稱Tj）和固有抖動(dòng)（Deterministic jitter，簡(jiǎn)稱Dj），驗(yàn)證是否滿足相關(guān)規(guī)范的要求。
在力科SDA系列示波器中使用了“Normalized Q-Scale method”（簡(jiǎn)稱NQ-Scale方法）來(lái)求解Tj。而Tj分解為固有抖動(dòng)Dj和隨機(jī)抖動(dòng)Rj時(shí)，力科SDA提供了三種抖動(dòng)分解方法，分別為Conventional、effective、MJSQ，如下圖所示。

圖一：力科SDA的三種抖動(dòng)分解方法

MJSQ方法在Fibre Channel規(guī)范已有定義（MJSQ代表Methodologies for jitter and signal quality specification），這種方法在串行數(shù)據(jù)的抖動(dòng)分析中被廣泛使用。在MJSQ文檔中，Tj是某一測(cè)量樣本數(shù)量下的TIE抖動(dòng)的峰峰值，由Rj和Dj組成，Dj是有邊界的，而Rj是沒(méi)有邊界的，其概率密度函數(shù)滿足高斯分布。Tj的直方圖使用dual-Dirac來(lái)建模。Dual-Dirac模型是由兩個(gè)滿足高斯分布的脈沖組成，左右兩個(gè)脈沖的均值為μL和μR，兩個(gè)脈沖的標(biāo)準(zhǔn)偏差都等于σ，Dj = μR - μL，Rj = σ，Tj@BER-12= 14 * Rj + Dj。如下圖二所示。

圖二：Dual-driac模型與MJSQ方法示意圖

力科SDA中的MJSQ方法直接處理PDF概率密度函數(shù)，使用兩個(gè)高斯分布的曲線分別擬合TIE直方圖的左右兩邊的尾部，調(diào)節(jié)高斯曲線的標(biāo)準(zhǔn)偏差讓曲線能盡量擬合TIE直方圖的尾部。 力科SDA的MJSQ分解方法基于傳統(tǒng)的MJSQ方法進(jìn)行了革新，兩個(gè)高斯分布的均值可以是不以Y軸對(duì)稱的，標(biāo)準(zhǔn)偏差也可以是不相等的。擬合的兩個(gè)高斯曲線的均值之差為Dj，標(biāo)準(zhǔn)偏差的平均值為Rj。
Effective方法是直接對(duì)浴盆曲線（bathtub curve）進(jìn)行分析，將Tj分解為Dj和Rj。這種方法與誤碼率測(cè)試儀（BERT）的方法相同。在使用NQ-Scale方法得到不同BER下的Tj后，由于Tj = Dj + α×Rj，在不同BER下的系數(shù)α是已知，可以推算出不同BER下的Dj和Rj。如右圖所示，BER=10e-12時(shí)α=14.069，BER=10e-10時(shí)α=12.723。
由于effective方法求解Dj和Rj與BERT相似，所以使用這種方法的計(jì)算結(jié)果可以與BERT的進(jìn)行對(duì)比。另外，effective方法求解的Dj和Rj與MJSQ方法的非常相似。

圖三：力科SDA中effective方法分解Tj

Conventional方法先計(jì)算Dj中的周期性抖動(dòng)Pj（periodic jitter）、數(shù)據(jù)相關(guān)性抖動(dòng)DDj（data dependent jitter）和占空比失真DCD（duty cycle distortion），然后Dj = DDj + Pj，Rj = (Tj – Dj) / 14。 在下圖四的流程圖中描繪了DDj、Pj、Tj、Rj的分析流程。

圖四：力科SDA中Conventional方法分解Tj
在conventional方法中，首先計(jì)算出每個(gè)數(shù)據(jù)邊沿和參考時(shí)鐘的偏差，即TIE，把TIE抖動(dòng)隨時(shí)間變化的趨勢(shì)描繪為一條曲線（即TIE track，又稱TIE trend），對(duì)TIE追蹤曲線做快速傅立葉變換（FFT），得到TIE抖動(dòng)的頻譜，通常周期性抖動(dòng)Pj和DDj是TIE頻譜中的峰值部分，隨機(jī)抖動(dòng)Rj是TIE頻譜中的底部，其頻譜范圍非常寬。將代表Pj的抖動(dòng)頻譜部分做反向傅立葉變換IFFT，可以得到各個(gè)頻點(diǎn)貢獻(xiàn)多大的抖動(dòng)值，如下圖五所示，在測(cè)試結(jié)果的左下角對(duì)各個(gè)頻率貢獻(xiàn)的Pj；函數(shù)F1是對(duì)TIE trend作FFT運(yùn)算，即TIE抖動(dòng)的頻譜，在F1中的峰值點(diǎn)即為周期性抖動(dòng)，在圖五中用cursor測(cè)量結(jié)果為350kHz，與左下角的Pj breakdown菜單的分析結(jié)果吻合。

圖五：力科SDA中Pj的分析

DDj的計(jì)算有兩種方法，一種是Synchronous N cycle plot方法，另一種是ISI plot方法。
當(dāng)捕獲的信號(hào)的數(shù)據(jù)碼流是重復(fù)的，使用Synchronous N cycle plot方法，如下圖六所示，
理想情況下，數(shù)據(jù)pattern的每個(gè)邊沿的平均位置μ應(yīng)該在bit interval的整數(shù)倍上，由于通道的非線性以及上升下降時(shí)間的不平衡，該平均位置μ（均值）與理想位置e存在時(shí)間差，這種抖動(dòng)是與數(shù)據(jù)的碼型相關(guān)的，記錄下一個(gè)pattern中每個(gè)邊沿的時(shí)間差offset= e－μ。其中的最大值減去最小值即為DDj抖動(dòng)的峰峰值。
圖六：Synchronous N cycle plot方法計(jì)算DDj

在下圖七中記錄了某偽隨機(jī)碼PRBS7的SnCycle圖，即127個(gè)bit中每個(gè)邊沿的DDj（offset= e－μ）的變化趨勢(shì)圖。相鄰兩個(gè)邊沿之間如果存在連續(xù)的0或者1（即沒(méi)有跳變位），則采用線性插值連接offset與offset。測(cè)量該DDj隨著bit位的變化的曲線的峰峰值即可得到DDj的峰峰值。

圖七：Synchronous N cycle plot方法計(jì)算DDj

當(dāng)捕獲的數(shù)據(jù)信號(hào)的碼流是不重復(fù)時(shí)，使用ISI plot方法來(lái)計(jì)算DDj。由于碼間干擾ISI是由于連接TX到RX的信道對(duì)碼型不同的信號(hào)產(chǎn)生不同的影響，在力科SDA中，首先構(gòu)建選定長(zhǎng)度為N個(gè)bit信號(hào)，幀長(zhǎng)度為N個(gè)bit的信號(hào)可以有2個(gè)組合。從串行數(shù)據(jù)流中找出同種bit組合的信號(hào)，平均運(yùn)算以去除隨機(jī)抖動(dòng)，然后把平均后的各種組合的碼型疊加在一起，可以測(cè)量到碼型相關(guān)抖動(dòng)DDj。在下圖八中，N=5，在左圖中可以看到，從串行數(shù)

圖八：力科SDA中ISI Plot方法分析DDJ

據(jù)碼流中選取“00010”和“11110”，取平均運(yùn)算后，疊加到ISI圖中，可以清晰觀察到兩個(gè)碼型導(dǎo)致的抖動(dòng)。在右上圖為串行數(shù)據(jù)的眼圖，右下圖為六種不同碼型疊加的ISI圖，前者的輪廓與后者完全一致，而后者在加入隨機(jī)抖動(dòng)后與前者很接近。
三種計(jì)算方法的比較：
Conventional方法可以計(jì)算出DDj、Pj、DCD、ISI，可以計(jì)算出Pj的來(lái)自于哪些頻率，對(duì)于串行信號(hào)的分析和調(diào)試非常實(shí)用。Effective和MJSQ方法只能得到Tj、Dj和Rj，不能把Dj進(jìn)一步分解。Effective方法采用與BERT相似的方法來(lái)計(jì)算Dj和Rj，測(cè)量結(jié)果可以與BERT做對(duì)比。
注：在三種方法的示意圖中都可以看到Tj是用TIE直方圖的尾部外插值后推算出來(lái)的，在抖動(dòng)測(cè)試儀器行業(yè)中，對(duì)于直方圖尾部擬合且外插值有幾種算法。 力科SDA使用了NQ-Scale方法來(lái)對(duì)TIE直方圖的尾部進(jìn)行擬合和外插值運(yùn)算，在另一篇文章中將介紹NQ-Scale方法。
參考文獻(xiàn)
1, Understanding the Choices for Jitter Calculation Method, LeCroy Application Brief
2, A Comparison of Methods for Estimating Total Jitter Concerning Precision, Accuracy and Robustness, Martin Miller Ph.D., Michael Schnecker, DesignCon2007.
3, Fibre Channel – Method Jitter and Signal Quality Specification – MJSQ, T11.2/Project 1315-DT/Rev 14.1, June 5, 2005.