劉彪，李變，楊劍青，朱張林，王惜康

?

頻率穩(wěn)定度實(shí)時(shí)評(píng)估系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn)

劉彪1,2，李變1,3，楊劍青1,2，朱張林1,2，王惜康1,2

（1. 中國(guó)科學(xué)院國(guó)家授時(shí)中心，西安 710600；2. 中國(guó)科學(xué)院大學(xué)，北京 100049；3. 中國(guó)科學(xué)院時(shí)間頻率基準(zhǔn)重點(diǎn)實(shí)驗(yàn)室，西安 710600）

為了實(shí)時(shí)計(jì)算原子鐘的頻率穩(wěn)定度，設(shè)計(jì)了一種頻率穩(wěn)定度實(shí)時(shí)評(píng)估系統(tǒng)。利用現(xiàn)有理論基礎(chǔ)（穩(wěn)定度測(cè)量原理及其算法等）以及軟件技術(shù)開發(fā)了一種實(shí)時(shí)計(jì)算并以圖形顯示原子鐘穩(wěn)定度的測(cè)量分析軟件。該軟件使用VC++的多媒體定時(shí)器實(shí)現(xiàn)數(shù)據(jù)的實(shí)時(shí)不間斷采集及計(jì)算，同時(shí)采用雙緩存繪圖技術(shù)顯示計(jì)算結(jié)果。測(cè)試結(jié)果表明該軟件工作穩(wěn)定，計(jì)算結(jié)果可靠、可信。

頻率穩(wěn)定度；雙緩存繪圖；Allan方差；雙對(duì)數(shù)坐標(biāo)系

0 引言

對(duì)于一臺(tái)高精度的原子鐘，如果不了解其頻率穩(wěn)定度特性[1-2]，就無法恰當(dāng)?shù)貞?yīng)用它。目前大部分計(jì)算穩(wěn)定度的軟件都屬于事后分析計(jì)算軟件，不能實(shí)時(shí)地發(fā)現(xiàn)問題，實(shí)時(shí)地加以控制和解決，并且事后分析需要大量的人工操作。另外在原子時(shí)實(shí)時(shí)計(jì)算領(lǐng)域，原子時(shí)的計(jì)算需要對(duì)多臺(tái)鐘的數(shù)據(jù)通過一定的算法處理才能實(shí)現(xiàn)。由于各臺(tái)鐘的性能不同，因此在原子時(shí)的計(jì)算中，需要估算它們之間的權(quán)重比值，而它們之間的權(quán)重比值與穩(wěn)定度的評(píng)估結(jié)果有關(guān)，所以為了得到原子時(shí)的實(shí)時(shí)計(jì)算結(jié)果，也需要對(duì)原子鐘的頻率穩(wěn)定度進(jìn)行實(shí)時(shí)評(píng)估。為此，我們研制了一種頻率穩(wěn)定度實(shí)時(shí)評(píng)估系統(tǒng)。

1 頻率穩(wěn)定度實(shí)時(shí)評(píng)估系統(tǒng)的硬件構(gòu)成

圖1是頻率穩(wěn)定度實(shí)時(shí)評(píng)估系統(tǒng)的硬件構(gòu)成圖[1]。

圖1 頻率穩(wěn)定度實(shí)時(shí)評(píng)估系統(tǒng)的硬件構(gòu)成圖

2 軟件功能及其設(shè)計(jì)與實(shí)現(xiàn)

頻率穩(wěn)定度實(shí)時(shí)計(jì)算軟件主要由數(shù)據(jù)采集模塊、數(shù)據(jù)預(yù)處理模塊、數(shù)據(jù)處理模塊、數(shù)據(jù)保存模塊以及圖形顯示模塊等5部分組成。下面主要對(duì)數(shù)據(jù)采集模塊、數(shù)據(jù)處理模塊和圖形顯示模塊作簡(jiǎn)要介紹。

2.1 數(shù)據(jù)采集模塊

數(shù)據(jù)采集模塊是采用串口通信技術(shù)實(shí)現(xiàn)的。目前，在Windows下編程時(shí)，常用的串口通信主要有3種方法：1）用Windows API進(jìn)行編程；2）用MSComm通信控件；3）用第三方提供的一些串口通訊類進(jìn)行編寫[3-5]。

利用Windows API編寫串口程序，特別是編寫復(fù)雜的多線程串口程序時(shí)，對(duì)于程序員的能力要求較高。除了需要程序員熟練掌握和使用眾多的API函數(shù)，能編寫很多底層代碼之外，還必須熟悉線程的編程方法，因此，一般不建議普通程序員使用這種方式。

在只需對(duì)單個(gè)串口操作的情況下，利用MSComm控件會(huì)使編程快捷簡(jiǎn)單。然而，由于做了大量的封裝，降低了編程的可控性和靈活性，因此在多線程多串口編程時(shí)，需要做許多復(fù)雜處理的情況下，最好不使用MSComm控件。

CerialPort類是一個(gè)采用Windows API技術(shù)開發(fā)的串口通信類。這個(gè)類基于多串口、多線程技術(shù)，只需要將這個(gè)類添加到編寫的程序框架中，就可以使用其成員函數(shù)開發(fā)程序。利用CSerialPort類進(jìn)行串口編程時(shí)，既可以使編程效率高，程序可控性強(qiáng)，擴(kuò)展性好，又比Windows API編程簡(jiǎn)單。

比較以上3種串口通信方式，可以發(fā)現(xiàn)使用第三方的CSerialPort類，是實(shí)現(xiàn)Windows下的多線程串口編程的較好的選擇[3-5]。

通常，對(duì)于單個(gè)串口的操作主要包括：打開串口、設(shè)置串口參數(shù)、使用串口（寫串口、讀串口）和關(guān)閉串口等。當(dāng)使用CSerialPort類時(shí)，其串口通信模型框圖如圖2所示。

圖2 CSerialPort類串口通信模型框圖

CSerialPort類工作流程是：對(duì)接收數(shù)據(jù)而言，首先設(shè)置好串口參數(shù)，再開啟串口監(jiān)測(cè)線程，串口監(jiān)測(cè)線程監(jiān)測(cè)到串口接收到的數(shù)據(jù)、流控制事件或其他串口事件后，就以消息方式通知主程序，激發(fā)消息處理函數(shù)。對(duì)發(fā)送數(shù)據(jù)而言，可直接向串口發(fā)送。

2.2 數(shù)據(jù)處理模塊

按觀測(cè)域的不同，頻率穩(wěn)定度分為時(shí)域頻率穩(wěn)定度和頻域頻率穩(wěn)定度。所謂時(shí)域頻率穩(wěn)定度，就是對(duì)受噪聲影響的輸出頻率用一個(gè)時(shí)間函數(shù)來描述，從時(shí)域的角度來分析噪聲對(duì)輸出頻率的影響。對(duì)于頻率穩(wěn)定度的時(shí)域表征，除了常用的阿倫方差（AVAR），還有修正阿倫方差（MVAR）、哈達(dá)瑪方差（HVAR）等[6-12]。在該數(shù)據(jù)處理模塊中，采用了多種方差對(duì)穩(wěn)定度進(jìn)行評(píng)估。在運(yùn)行軟件后，可以對(duì)方差類型進(jìn)行選擇，根據(jù)不同的鐘，去選擇不同的算法對(duì)鐘進(jìn)行分析。

2.3 圖形顯示模塊

對(duì)于圖形顯示模塊，下面主要介紹雙對(duì)數(shù)坐標(biāo)系的繪制和雙緩存繪圖技術(shù)。

2.3.1 雙對(duì)數(shù)坐標(biāo)系的繪制

該雙對(duì)數(shù)坐標(biāo)系可以以任意數(shù)據(jù)為起點(diǎn)、任意數(shù)據(jù)為終點(diǎn)，并能對(duì)坐標(biāo)進(jìn)行縮放。在這里我們定義1×10（=…，-2，-1，0，1，2，3，…）的相應(yīng)數(shù)據(jù)點(diǎn)為標(biāo)準(zhǔn)數(shù)據(jù)。在對(duì)數(shù)坐標(biāo)系中，任意3個(gè)數(shù)據(jù)值1，2，3和在屏幕上的對(duì)數(shù)坐標(biāo)系中的對(duì)應(yīng)點(diǎn)的坐標(biāo)1，2，3是一一對(duì)應(yīng)的（如圖3所示）。它

圖3 數(shù)據(jù)和屏幕上點(diǎn)的坐標(biāo)對(duì)應(yīng)關(guān)系圖

們之間的關(guān)系可表達(dá)為

在繪圖時(shí)，首先要找出起始數(shù)據(jù)（數(shù)據(jù)最小值）和終止數(shù)據(jù)（數(shù)據(jù)最大值），然后確定起始標(biāo)準(zhǔn)數(shù)據(jù)和終止標(biāo)準(zhǔn)數(shù)據(jù)，畫出標(biāo)準(zhǔn)數(shù)據(jù)線，然后根據(jù)標(biāo)準(zhǔn)數(shù)據(jù)線畫出其他數(shù)據(jù)線。起始標(biāo)準(zhǔn)數(shù)據(jù)就是第一個(gè)大于起始數(shù)據(jù)的標(biāo)準(zhǔn)數(shù)據(jù)，終止標(biāo)準(zhǔn)數(shù)據(jù)就是最后一個(gè)小于終止數(shù)據(jù)的標(biāo)準(zhǔn)數(shù)據(jù)。

根據(jù)我們計(jì)算結(jié)果的數(shù)據(jù)，有2種繪圖情況：

1）起始數(shù)據(jù)≤起始標(biāo)準(zhǔn)數(shù)據(jù)<終止標(biāo)準(zhǔn)數(shù)據(jù)≤終止數(shù)據(jù)，如圖4所示。在這種情況下，在確定坐標(biāo)系范圍后，先畫出起始標(biāo)準(zhǔn)數(shù)據(jù)線和終止標(biāo)準(zhǔn)數(shù)據(jù)線。對(duì)于起始數(shù)據(jù)和起始標(biāo)準(zhǔn)數(shù)據(jù)之間的部分，由起始標(biāo)準(zhǔn)數(shù)據(jù)（設(shè)起始標(biāo)準(zhǔn)數(shù)據(jù)為1×10）向起始數(shù)據(jù)方向遞減地畫出表示9×10-1，8×10-1，7×10-1，…的虛線；對(duì)于終止標(biāo)準(zhǔn)數(shù)據(jù)和終止數(shù)據(jù)之間的部分，由終止標(biāo)準(zhǔn)數(shù)據(jù)（設(shè)終止標(biāo)準(zhǔn)數(shù)據(jù)為1×10）向終止數(shù)據(jù)方向遞增地畫出表示2×10+1，3×10+1，4×10+1，…的虛線；對(duì)于起始標(biāo)準(zhǔn)數(shù)據(jù)和終止標(biāo)準(zhǔn)數(shù)據(jù)之間的部分，首先確定二者之間的標(biāo)準(zhǔn)數(shù)據(jù)，然后根據(jù)標(biāo)準(zhǔn)數(shù)據(jù)畫出其間的虛線。

圖4 繪圖情況1

2）起始數(shù)據(jù)<起始標(biāo)準(zhǔn)數(shù)據(jù)=終止標(biāo)準(zhǔn)數(shù)據(jù)<終止數(shù)據(jù)，如圖5所示。這種情況是第1種情況的特例，由于起始標(biāo)準(zhǔn)數(shù)據(jù)=終止標(biāo)準(zhǔn)數(shù)據(jù)，相比于第1種情況，減少了畫出起始標(biāo)準(zhǔn)數(shù)據(jù)與終止標(biāo)準(zhǔn)數(shù)據(jù)之間部分的虛線的過程。

圖5 繪圖情況2

當(dāng)坐標(biāo)系確定后，根據(jù)式（1），可以確定任意數(shù)據(jù)點(diǎn)在屏幕上對(duì)數(shù)坐標(biāo)系中的位置。縱坐標(biāo)的處理同上。

2.3.2 雙緩存繪圖

在GDI的繪圖系統(tǒng)中，每一次調(diào)用區(qū)域繪圖操作，圖形顯示系統(tǒng)就會(huì)在屏幕中對(duì)指定的區(qū)域進(jìn)行刷新操作。如果頻繁地進(jìn)行區(qū)域繪制操作的話，就會(huì)出現(xiàn)屏幕閃爍。產(chǎn)生閃屏的原因是每次在進(jìn)行DC繪圖時(shí)都要把相關(guān)的顯示數(shù)據(jù)發(fā)送到顯卡，經(jīng)顯卡處理后在顯示器上顯示。我們可以通過減少數(shù)據(jù)與顯卡之間的交互次數(shù)來提高繪圖效率。雙緩存繪圖技術(shù)有利于改善閃屏情況，因?yàn)樵陔p緩存繪圖技術(shù)中，首先把更新操作中的所有繪制數(shù)據(jù)寫入內(nèi)存，然后調(diào)用函數(shù)BitBlt一次性地把所有數(shù)據(jù)發(fā)送到顯卡上[13]。

2.4 軟件自動(dòng)實(shí)時(shí)性的實(shí)現(xiàn)

軟件的自動(dòng)實(shí)時(shí)性，主要是通過采用VC++中的多媒體定時(shí)器技術(shù)來實(shí)現(xiàn)的。VC++定時(shí)器的定時(shí)精度可以達(dá)到毫秒級(jí)，因此在本軟件中可以達(dá)到實(shí)時(shí)采集并計(jì)算的要求[14]。具體過程如下：用戶打開軟件后，首先對(duì)參數(shù)進(jìn)行設(shè)置，設(shè)置的參數(shù)包括數(shù)據(jù)的采樣間隔（多長(zhǎng)時(shí)間采集一個(gè)數(shù)據(jù)）和系統(tǒng)的計(jì)算間隔（系統(tǒng)多長(zhǎng)時(shí)間進(jìn)行一次數(shù)據(jù)的保存、計(jì)算并顯示），這2個(gè)參數(shù)的單位都是秒。然后對(duì)下拉列表框中的方差類型進(jìn)行選擇，點(diǎn)擊自動(dòng)計(jì)算按鈕，軟件就開始對(duì)發(fā)送到串口的數(shù)據(jù)進(jìn)行定時(shí)的采集，當(dāng)軟件運(yùn)行的時(shí)間到達(dá)我們所設(shè)置的時(shí)間間隔時(shí)，軟件就會(huì)把顯示在對(duì)話框中的數(shù)據(jù)保存到名為FHASEData00.txt，F(xiàn)HASEData01.txt，F(xiàn)HASEData02.txt…的文件中，在這里我們?cè)O(shè)置了延遲時(shí)間參數(shù)為100ms，經(jīng)過延遲時(shí)間，我們對(duì)保存在FHASEData00.txt，F(xiàn)HASEData01.txt，F(xiàn)HASEData02.txt…的文件中的數(shù)據(jù)進(jìn)行預(yù)處理，把計(jì)算結(jié)果保存到一個(gè)名為FASEDATA.txt的文件中，再調(diào)用數(shù)據(jù)處理模塊中的時(shí)域頻率穩(wěn)定度算法對(duì)數(shù)據(jù)進(jìn)行處理，把計(jì)算結(jié)果保存到一個(gè)名為SATADATA.txt的文件中，最后我們調(diào)用圖形顯示模塊中的算法對(duì)圖形進(jìn)行顯示。這就是軟件的單次運(yùn)行流程，如圖6所示。定時(shí)器的設(shè)置，會(huì)讓軟件重復(fù)執(zhí)行該流程，達(dá)到定時(shí)刷新的效果?？傊?，系統(tǒng)會(huì)不間斷地對(duì)數(shù)據(jù)進(jìn)行采集、保存、運(yùn)算并顯示。

圖6 軟件單次運(yùn)行流程圖

3 軟件測(cè)試及其結(jié)果

該軟件可以實(shí)時(shí)自動(dòng)地采集、保存和計(jì)算，并將結(jié)果以圖形化的方式直觀地展現(xiàn)給用戶。我們對(duì)該軟件進(jìn)行以下2個(gè)方面的測(cè)試：1）軟件的算法的準(zhǔn)確性；2）軟件的實(shí)時(shí)性、自動(dòng)性。

stable32[11]是業(yè)內(nèi)比較權(quán)威的穩(wěn)定度分析軟件，我們用本系統(tǒng)的計(jì)算結(jié)果與stable32的計(jì)算結(jié)果作對(duì)比，驗(yàn)證本文軟件的準(zhǔn)確性。

1）對(duì)軟件的算法的準(zhǔn)確性進(jìn)行驗(yàn)證。

我們對(duì)2013年1月6日09:29:25到2013年1月7日09:30:52商品銣鐘（型號(hào)FS725）相對(duì)于主鐘銫鐘（實(shí)驗(yàn)室編號(hào)Cs2145）的鐘差數(shù)據(jù)進(jìn)行分析，其中數(shù)據(jù)采樣間隔為3s，由于銣鐘存在頻率漂移，用Hadamard方差分析，可以消除線性漂移頻率漂移對(duì)頻率穩(wěn)定度的影響，因此我們采用Hadamard方差進(jìn)行分析，圖7為計(jì)算所得穩(wěn)定度曲線。

圖7 評(píng)估系統(tǒng)與stable32對(duì)銣鐘FS725的分析結(jié)果對(duì)比

對(duì)于Cs鐘，由于不存在頻率漂移，一般采用Allan系列方差進(jìn)行分析，下面我們對(duì)銫鐘（實(shí)驗(yàn)室編號(hào)Cs98）相對(duì)于鐘房主鐘的1d（2013年10月6日10:15:25至10月7日10:15:25）的鐘差數(shù)據(jù)進(jìn)行分析，采樣間隔為1s。圖8為用Allan方差表征的穩(wěn)定度曲線。

圖8 評(píng)估系統(tǒng)與stable32對(duì)Cs98的分析結(jié)果對(duì)比（用1d的鐘差數(shù)據(jù)）

下面我們對(duì)銫鐘（Cs98）相對(duì)于鐘房主鐘的一個(gè)月（2011年2月）的鐘差數(shù)據(jù)進(jìn)行分析，采樣間隔為3600s，圖9為計(jì)算得的Allan方差所表征的穩(wěn)定度曲線。

圖9 評(píng)估系統(tǒng)與stable32對(duì)CS98分析結(jié)果對(duì)比（用一個(gè)月的鐘差數(shù)據(jù)）

由圖7至圖9可知評(píng)估系統(tǒng)得到的穩(wěn)定度結(jié)果與stable32的計(jì)算結(jié)果一致。

2）對(duì)軟件的實(shí)時(shí)性、自動(dòng)性進(jìn)行驗(yàn)證。

在測(cè)量時(shí)，首先運(yùn)行軟件，在界面上進(jìn)行參數(shù)設(shè)置，參數(shù)的設(shè)置界面如圖10所示。參數(shù)設(shè)置包括串口參數(shù)設(shè)置、計(jì)算間隔設(shè)置、采樣間隔設(shè)置以及方差類型（如AVAR、MVAR、HVAR等）的選擇，然后點(diǎn)擊按鈕進(jìn)行自動(dòng)計(jì)算，數(shù)據(jù)采集模塊就會(huì)對(duì)鐘差數(shù)據(jù)進(jìn)行采集，并根據(jù)設(shè)定的計(jì)算間隔進(jìn)行鐘差數(shù)據(jù)的保存并計(jì)算，計(jì)算結(jié)果也會(huì)隨之自動(dòng)保存并生成穩(wěn)定度曲線圖。

圖10 軟件參數(shù)設(shè)置界面圖

下面我們選擇商品銣鐘（FS725）相對(duì)于主鐘（銫鐘Cs2145）的鐘差數(shù)據(jù)進(jìn)行分析，其數(shù)據(jù)采樣間隔為3s，計(jì)算間隔為120s（可以自行設(shè)置），由于銣鐘存在頻率漂移，用Hadamard方差分析，可以消除線性漂移頻率漂移對(duì)頻率穩(wěn)定度的影響，因此我們采用Hadamard方差進(jìn)行分析。圖11是該系統(tǒng)運(yùn)行1d過程中自動(dòng)形成的頻率穩(wěn)定度的部分截圖。

以上測(cè)試表明本軟件計(jì)算結(jié)果準(zhǔn)確、可靠，并且彌補(bǔ)了stable32不能實(shí)時(shí)地對(duì)原子鐘穩(wěn)定度進(jìn)行評(píng)估的不足，完全滿足了自動(dòng)控制的要求，只需人工設(shè)置參數(shù)，就可以自動(dòng)對(duì)原子鐘的穩(wěn)定度進(jìn)行評(píng)估。本軟件還可以圖形顯示鐘差數(shù)據(jù)，自動(dòng)地判別鐘的噪聲類型，這方面內(nèi)容暫不在此詳述。

4 結(jié)語(yǔ)

本文從實(shí)時(shí)性的要求出發(fā)，設(shè)計(jì)并實(shí)現(xiàn)了一種頻率穩(wěn)定度實(shí)時(shí)評(píng)估系統(tǒng)。該系統(tǒng)可以對(duì)原子鐘的穩(wěn)定度進(jìn)行實(shí)時(shí)的評(píng)估和圖形顯示。通過測(cè)試表明：該系統(tǒng)穩(wěn)定、可靠，但仍然存在局限性，例如不能實(shí)時(shí)地對(duì)多臺(tái)原子鐘頻率穩(wěn)定度進(jìn)行評(píng)估。計(jì)劃在以后的工作中予以進(jìn)一步完善。

[1] 胡永輝, 漆貫榮. 時(shí)間測(cè)量原理[M]. 香港: 香港亞太科學(xué)出版社, 2000.

[2] 蔡勇. 頻率穩(wěn)定度測(cè)試儀軟件的設(shè)計(jì)與實(shí)現(xiàn)[D]. 北京: 北京郵電大學(xué), 2011.

[3] 蘇岳龍, 李貽斌, 宋銳. 基于VC++ 6.0的高速串口通信數(shù)據(jù)采集系統(tǒng)[J]. 微計(jì)算機(jī)信息, 2005, 21(5): 147-148.

[4] 胡峪, 劉靜. VC++高級(jí)編程技巧與示例[M]. 西安: 西安電子科技大學(xué)出版社, 2001.

[5] 楊峰, 桂衛(wèi)華, 陳鋒, 等. Win32環(huán)境下基于VC++ 6.0串口通信編程方法[J]. 計(jì)算技術(shù)與自動(dòng)化, 2004, 23(2): 70-73.

[6] 張暉, 金玉雯, 李延. 頻率穩(wěn)定度的時(shí)域表征_阿侖方差[J]. 工業(yè)計(jì)量, 2001(增刊): 192-193.

[7] 郭海榮, 郭樹人, 焦文海, 等. 原子鐘時(shí)域頻率穩(wěn)定性計(jì)算方法[J]. 飛行器測(cè)控學(xué)報(bào), 2007, 26(2): 47-52.

[8] 楊洋, 董緒榮, 李曉宇, 等. 星載原子鐘性能方法比較研究[J]. 科學(xué)技術(shù)與工程, 2012, 12(36): 9801-9804.

[9] GALLEANI L, TAVELLA P. The characterization of clock behavior with the dynamic Allan variance[J]. Proceeding of the 2003 Joint IEEE FCS/EFTF Meeting. 2003: 239-244.

[10] GALLEANI L, TAVELLA P. Tracking nonstationarities in clock noises using the dynamic Allan variance[J]. Proceeding of the 2005 Joint IEEE FCS/PTTI Meeting. 2003: 329-396.

[11] RILEY W J. Stable32 Frequency Stability Analysis USER MANUAL[K]. USA: 2005.

[12] 張敏. 原子鐘噪聲類型和頻率穩(wěn)定度估計(jì)的自由度分析與探討[D]. 西安: 中國(guó)科學(xué)院國(guó)家授時(shí)中心, 2008.

[13] 張磊. 基于VC++的高效繪圖—雙緩沖技術(shù)[J]. 信息科學(xué), 2009, 20: 61.

[14] 卓紅艷, 趙平. 基于VC++的實(shí)時(shí)數(shù)據(jù)采集系統(tǒng)中定時(shí)器的使用與比較[J]. 現(xiàn)代電子技術(shù), 2007, 18(257): 80-82.

Design and implementation of frequency stability real-time evaluation system

LIU Biao1,2, LI Bian1,3, YANG Jian-qing1,2, ZHU Zhang-lin1,2, WANG Xi-kang1,2

(1.National Time Service Center, Chinese Academy of Sciences, Xi′an 710600, China;2. University of Chinese Academy of Sciences, Beijing 100049, China;3. Key Laboratory of Time and Frequency Primary Standard, Chinese Academy of Sciences, Xi′an710600, China)

A frequency stability real-time evaluation system was designed for estimating the stability of atomic clock in real time. On the basis of the existing theory (measurement principle and algorithm of stability) and the software technology, a measuring-analyzing software was developed for estimating and displaying graphically the frequency stability of atomic clocks in real-time. In the software, theVC++multimedia timer is used forcontinuousacquisitionandcalculation of real-time data, and the double buffer memory drawing technique is adopted for displaying the calculationresults. The testresults show that thesoftwareis stable in operation and the calculation is reliable and credible.

frequency stability; double buffer memory drawing; AVAR; double logarithm coordinate system

TM935.12

A

1674-0637(2014)03-0137-08

10.13875/j.issn.1674-0637.2014-03-0137-08

2013-09-10

國(guó)家自然科學(xué)基金青年基金資助項(xiàng)目（11103025）；中國(guó)科學(xué)院“西部之光”人才培養(yǎng)計(jì)劃西部博士資助項(xiàng)目（中國(guó)科學(xué)院人教[2011]180號(hào)）

劉彪，男，碩士，主要從事原子鐘性能分析。