劉浩,劉丹丹,阮軍,王心亮,張首剛
(1.中國科學(xué)院 國家授時中心,西安 710600;2.中國科學(xué)院 時間頻率基準重點實驗室,西安 710600;3.中國科學(xué)院大學(xué),北京 100049)
銫原子噴泉鐘以脈沖方式工作,時序控制系統(tǒng)通過輸出控制信號,同步調(diào)控光學(xué)系統(tǒng)、微波頻率綜合器,實現(xiàn)原子的冷卻—上拋—選態(tài)—微波激勵—探測等周期性運行。通過數(shù)據(jù)采集系統(tǒng)采集銫原子的飛行時間信號,經(jīng)過計算獲得伺服控制系統(tǒng)的參數(shù),實現(xiàn)頻率鎖定[1]。目前運行在NTSC-F1銫原子噴泉鐘上的時序控制系統(tǒng),由于其硬件版本較低,系統(tǒng)整體性能提升空間不足。同時其系統(tǒng)通用性不理想,時序的各項參數(shù)設(shè)置修改比較繁瑣,所以需要研制一套界面友好,操作人性化,性能更優(yōu)的時序控制系統(tǒng)應(yīng)用于NTSC-F2銫原子噴泉鐘上。
本文針對NTSC-F1時序控制系統(tǒng)存在的不足,并依據(jù)NTSC-F2銫原子噴泉鐘對控制時序的要求,基于Labwindows/CVI平臺,通過系統(tǒng)模擬、硬件選型、軟件設(shè)計、性能測試等流程完成了時序控制系統(tǒng)的研制,實現(xiàn)了19路同步時序的控制與輸出。
時序控制系統(tǒng)需要輸出特定的時序,控制光開關(guān)、磁場開關(guān)、聲光調(diào)制器、壓控振蕩器等的關(guān)斷,實現(xiàn)銫原子在高真空、高磁屏蔽的物理環(huán)境中冷卻—上拋—選態(tài)—微波激勵—探測等過程的周期性運行[2]。時序控制系統(tǒng)首先通過輸出TTL信號控制二維磁光阱(2D MOT)冷卻光和磁場開啟,制備水平方向慢速原子束。同時打開三維磁光阱(3D MOT)冷卻光和磁場,等慢速原子束注入3D MOT真空腔后,俘獲并囚禁銫原子束,實現(xiàn)銫原子團裝載。裝載完畢后關(guān)閉磁光阱內(nèi)磁場,通過激光對原子團進行冷卻。隨后改變上下兩組俘獲光失諧,實現(xiàn)銫原子的上拋。上拋約1.5 ms后通過時序控制激光頻率失諧量的增加和激光功率的減小,實現(xiàn)對原子團的進一步冷卻。之后,利用抽運激光將上拋中的銫原子全部抽運至|F=4>態(tài)上。|F=4>態(tài)的銫原子上拋進入選態(tài)腔,經(jīng)由選態(tài)腔中的微波場作用,使得|F=4,mF=0>態(tài)的原子躍遷到|F=3,mF=0>態(tài)上。躍遷后的銫原子飛出選態(tài)腔后,打開選態(tài)推光將|mF≠0>態(tài)的銫原子推走。經(jīng)選態(tài)后的銫原子再進入激勵腔與微波場發(fā)生作用,隨后銫原子繼續(xù)上升達到頂峰后開始回落,在回落過程中再次進入激勵腔,再次與微波場發(fā)生作用,完成整個Ramsey 躍遷過程。銫原子團繼續(xù)下落進入探測區(qū)進行飛行時間信號的探測。利用數(shù)據(jù)采集系統(tǒng)計算銫原子躍遷幾率[3-5]。調(diào)控微波頻率綜合器的頻率獲得Ramsey躍遷條紋,在Ramsey中心頻率附近進行調(diào)制,根據(jù)其兩側(cè)躍遷幾率差值得到誤差控制信號,經(jīng)PID運算后輸出伺服控制晶振,實現(xiàn)頻率鎖定。NTSC-F2時序控制系統(tǒng)工作原理如圖1所示。
圖1 時序控制系統(tǒng)工作圖
NTSC-F2銫原子噴泉鐘的運行由19路時序信號來控制,其時序波形如圖2所示。這19路時序由15路TTL信號和4路模擬信號組成,依次編號為第1路,第2路,……,第19路。其中信號周期為2 s,信號延時為t,用時為T。
圖2 NTSC-F2的控制時序
基于LabWindows/ CVI 為平臺設(shè)計時序控制模擬系統(tǒng)界面,通過C語言進行編程,利用控件觸發(fā)程序運行,實現(xiàn)波形輸出與顯示。具體流程為,通過LabWindows/ CVI搭建系統(tǒng)界面,并設(shè)置各控件初始參數(shù)。再利用GetCtrlVal函數(shù)獲取界面輸入?yún)?shù),包括起始電平、用時、延時、周期和相位等參數(shù)。再由算法將所獲得的波形參數(shù)轉(zhuǎn)化為信號生成SquareWave函數(shù)所需參數(shù),并由SquareWave函數(shù)生成時序,其公式表示為
X[i]=V0×SquareWave (φ0+f×360×i),
(1)
{SquareWave[i]=1.0 0 SquareWave[i]=-1.0 (duty/100.0×360.0) 式(1)中,數(shù)組X[]用于存儲一個周期的信號數(shù)據(jù)點,V0代表起始電平,φ0代表初始相位,f代表頻率,modulo為取模運算,duty代表占空比。 最后利用函數(shù)PlotY將各路波形顯示在Graph控件上。時序控制模擬系統(tǒng)運行效果如圖3所示,信號波形與期望值基本一致,并可以通過參量設(shè)置實現(xiàn)對時序信號的控制。 圖3 時序控制模擬系統(tǒng) 在銫原子噴泉鐘的運行過程中,原子冷卻、上拋,探測過程中激光頻率、功率的變化,以及時序信號控制的各開關(guān)關(guān)斷需要在極短的時間內(nèi)完成,這對硬件性能提出了較高的要求。選擇NI公司基于PXI總線的工控機作為硬件模塊的核心,PXI總線的傳輸速率為132~264 MB/s,各模塊板卡接入10 MHz參考時鐘,能夠有效進行多板卡信號的同步處理[6]。PXIe1082工控機內(nèi)部集成嵌入式控制器PXIe-8840,采用兩塊計數(shù)器卡PXI-6602輸出15路TTL信號,一塊模擬輸出卡PXI-6733輸出4路任意波模擬信號,并通過多功能I/O卡PXIe-6361來實現(xiàn)信號采集。 PXIe-8840是一款基于Intel Core i5-4400E處理器的高性能嵌入式控制器,具有2.7 GHz基頻、3.3 GHz雙核處理器和單通道1 600 MHz DDR3內(nèi)存,極適合用于處理器密集型、模塊化儀器和數(shù)據(jù)采集應(yīng)用。PXI-6602卡有5 V,8通道的計數(shù)器/定時器,用于輸出時序中的TTL電平信號。PXI-6733卡有8路模擬輸出通道,可同時控制8路不同硬件,單通道輸出更新率最高可達1 Ms/s[7-8]。將上述板卡插入PXIe-1082機箱內(nèi),通過配套的線纜接入接線盒,再由BNC接連接線盒輸出端與被控子系統(tǒng),這便構(gòu)成了時序控制系統(tǒng)的硬件模塊,其總體結(jié)構(gòu)如圖4所示。 圖4 硬件模塊框架圖 軟件系統(tǒng)用于實現(xiàn)對硬件模塊的控制,基于時序控制模擬系統(tǒng)進行開發(fā),在工控機平臺上,通過程序算法調(diào)用并控制硬件模塊輸出特定的19路時序信號,實現(xiàn)各開關(guān)和子系統(tǒng)按照特定的方式運行,完成銫原子噴泉鐘的穩(wěn)定工作。 首先借助時序控制模擬系統(tǒng)開發(fā)的系統(tǒng)界面,完善Command Button、Numeric、Text等控件作為時序參量輸入、時序開關(guān)和輔助功能調(diào)用鍵。隨后針對各控件編寫相應(yīng)的回調(diào)函數(shù),由DAQmxCreateAOVoltageChan配置物理通道,再通過GetCtrlVal函數(shù)獲取各時序波形參量,包括時序的起始電平、用時、延時和相位等,并利用算法將這些參量再轉(zhuǎn)化為信號生成函數(shù)所需的各參數(shù)。通過DAQmxCreateCOPulseChanTicks函數(shù)配置計數(shù)器源、TTL電平的極性、高電平時間、低電平時間、延遲時間。同時利用DAQmxCfgSampClkTiming函數(shù)設(shè)置采樣時鐘源和采樣速率,并通過DAQmxCfgDigEdgeStartTrig函數(shù)配置觸發(fā)信號[9]。最后根據(jù)偏振梯度冷卻時間、采樣速率將一個周期的波形數(shù)據(jù)寫入數(shù)組中。當(dāng)接收到觸發(fā)信號時,則將數(shù)據(jù)順序輸出,實現(xiàn)時序輸出。一個周期運行完畢,將其放入循環(huán)系統(tǒng)中便能實現(xiàn)信號周期性輸出。通道采用同一個參考時鐘源,以保證各通道間時間順序上的確定關(guān)系。時序輸出的具體流程圖如圖5所示。 圖5 時序生成流程圖 圖5中4路模擬時序通過GenSquareWave函數(shù)生成,其公式表示為 X[i]=V0×SquareWave (φ0+f×360×i), (2) {SquareWave[i]=(phase_i/360.0<=duty/100.0) ?V0:-V0} 式(2)中,數(shù)組X[]用于存儲一個周期的模擬信號數(shù)據(jù)點,V0代表起始電平,φ0代表初始相位,f代表頻率,duty代表占空比。 4路模擬信號每個周期都有一段用時(高電平向低電平跳變時所需時間),利用程序 A=(yongshi/zhouqi)×numElements; B=(yanshi/zhouqi)×numElements; for(i=k;i {SquareWave[i]=SquareWave[i+1]-DP×(i-A)/numElements+DP;} 和2.2 硬件選型
2.3 軟件設(shè)計