熊擇正,袁志超,賀軒,方聲偉,王青,齊向暉,陳徐宗
(北京大學(xué) 電子學(xué)院量子電子學(xué)研究所,北京 100871)
鐘差預(yù)報技術(shù)在守時授時領(lǐng)域具有重要的意義:衛(wèi)星導(dǎo)航系統(tǒng)等精密授時應(yīng)用需要獲取原子鐘與地面基準(zhǔn)或協(xié)調(diào)世界時UTC等參考間的鐘差數(shù)據(jù)以提高精確度,研究鐘差預(yù)報能以相對較小的代價提高數(shù)據(jù)的精度和可靠性[1];同時,鐘差預(yù)報技術(shù)應(yīng)用在原子鐘駕馭時可以明顯提高所馴服原子鐘的長期穩(wěn)定度[2];在建立時間基準(zhǔn)的過程中,也需要對所用多臺原子鐘之間的鐘差數(shù)據(jù)進(jìn)行預(yù)測[3]。目前衛(wèi)星鐘差的預(yù)報算法得到了比較深入的研究,并有各類針對不同數(shù)據(jù)特點(diǎn)的預(yù)報模型被提出。典型的衛(wèi)星鐘差預(yù)報模型包括多項式模型[4]、灰色模型GM(1,1)[5]、卡爾曼模型[6]、神經(jīng)網(wǎng)絡(luò)模型[7]等。
目前被廣泛研究的鐘差數(shù)據(jù)是原子鐘與時間基準(zhǔn)之間或者多臺原子鐘之間的計時差,因此鐘差預(yù)報理論不能直接應(yīng)用于單個原子鐘的性能提升。但是從原子鐘內(nèi)部也可以讀取出一些與鐘差數(shù)據(jù)具有相似的波動特征的信號,例如表示原子鐘內(nèi)訊問信號頻率與參考譜線頻率差值的誤差信號等。使用和鐘差預(yù)報算法類似的方法可以對誤差信號進(jìn)行預(yù)測,預(yù)測結(jié)果可以應(yīng)用于調(diào)節(jié)原子鐘的鎖定方式并提升其性能。此外,銫爐溫度、微波功率、C場電流等參數(shù)的波動會顯著影響原子鐘的穩(wěn)定性[8],對這些參數(shù)進(jìn)行實時的監(jiān)測和預(yù)報也有助于采取更加合適的頻率反饋方式,有望提高原子鐘的頻率穩(wěn)定度。
本文針對原子鐘內(nèi)部誤差信號的預(yù)報展開研究,采用神經(jīng)網(wǎng)絡(luò)模型對誤差信號進(jìn)行預(yù)報,并對所使用的神經(jīng)網(wǎng)絡(luò)模型的預(yù)測性能進(jìn)行了評估和優(yōu)化。論文提出一種基于局域線性小波神經(jīng)網(wǎng)絡(luò)(local linear wavelet neural network,LLWNN)的原子鐘內(nèi)部信號預(yù)測算法,利用北京大學(xué)原子鐘小組光抽運(yùn)小銫鐘內(nèi)的誤差信號數(shù)據(jù)對神經(jīng)網(wǎng)絡(luò)進(jìn)行訓(xùn)練和測試,并構(gòu)建了基于BP(back propagation)神經(jīng)網(wǎng)絡(luò)的信號預(yù)測算法作為對比。論文對基于LLWNN和BP神經(jīng)網(wǎng)絡(luò)的兩種算法進(jìn)行了性能評估和優(yōu)化,在各自優(yōu)化到最佳結(jié)構(gòu)之后,基于LLWNN的預(yù)測算法表現(xiàn)出了更高的預(yù)測精度和更好的預(yù)測穩(wěn)定度,有望在未來應(yīng)用中提升原子鐘的頻率穩(wěn)定度。
銫原子鐘內(nèi)的誤差信號可以近似地表示為
S=-K(vLO-vREF)+δs=-K[vLO-(vhfs+δv)]+δs,
(1)
式(1)中:vLO是微波訊問信號的頻率值,而微波訊問信號由本地晶振經(jīng)過頻率綜合電路產(chǎn)生,其頻率和晶振頻率成正比;vREF=vhfs+δv是銫鐘內(nèi)晶振鎖定參考譜線的中心頻率,受到光功率和微波功率漂移、環(huán)境溫度變化等因素的影響,vREF相對于銫鐘超精細(xì)躍遷頻率的標(biāo)稱值vhfs具有偏離量δv,其波動以長期噪聲為主。原子鐘內(nèi)部通過一系列調(diào)制解調(diào)過程來測量vLO與vREF的頻率差并將其以誤差信號的形式在電路中輸出,這些過程的總增益K在頻率波動很小時可以近似成常數(shù);但是,原子數(shù)散彈噪聲、光功率和微波功率波動、光子散彈噪聲等還會為上述測量帶來不確定度δs[9],這些噪聲會同時帶來短期和長期上的測量誤差,并且在短期產(chǎn)生尤其顯著的影響??傮w而言,銫原子鐘內(nèi)的誤差信號包含白相位噪聲、閃變相位噪聲、白頻率噪聲、閃變頻率噪聲和隨機(jī)行走頻率噪聲5類噪聲的疊加[10],而且短期波動的幅度很大。在原子鐘的反饋機(jī)制調(diào)節(jié)下,誤差信號的長期波動已經(jīng)很小,本文主要對誤差信號進(jìn)行較短時間尺度下的預(yù)報,這與鐘差信號的中長期預(yù)報相比具有很多新的困難。
從原子鐘相位的角度來看,Δt時間內(nèi)原子鐘的輸出相位變化可以表示為
(2)
在鐘差預(yù)報中一般是已知了歷史相位數(shù)據(jù),通過歷史相位數(shù)據(jù)預(yù)測未來的鐘差偏移以提前做出修正,這種方式的前提是通過一個外部的時間基準(zhǔn)作為參考測量出原子鐘的歷史相位。本文討論的場景是在原子鐘的內(nèi)部系統(tǒng)分析參數(shù),此時沒有與外部參考比對而得出的歷史相位數(shù)據(jù),但是分析誤差信號也可以在一定程度上獲取原子鐘的頻率和相位偏移信息;在δs,δv都是小量的中長期范圍內(nèi),預(yù)測誤差信號并對本振頻率做出相應(yīng)的調(diào)整也有利于提高原子鐘的頻率穩(wěn)定度。
在影響銫束原子鐘內(nèi)部誤差信號的諸多參數(shù)中,微波功率、光功率變化等因素導(dǎo)致的信號波動規(guī)律相對簡單一些,對這類信號波動進(jìn)行預(yù)報的難度相對較低;其他類型的信號波動則很難預(yù)報,例如溫度、振動等因素則可能通過多種機(jī)制以更復(fù)雜的方式來擾動信號,原子數(shù)散彈噪聲、光子散彈噪聲等白頻率噪聲則具有相當(dāng)復(fù)雜的波動模式。為了更好地預(yù)報前一類信號波動,本文采用理論上可以逼近任意非線性時間序列的神經(jīng)網(wǎng)絡(luò)[11]來模擬這些波動函數(shù)。
圖1 BP 神經(jīng)網(wǎng)絡(luò)和LLWNN 神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)示意圖
在后期研究中我們希望將神經(jīng)網(wǎng)絡(luò)的訓(xùn)練和預(yù)測過程在原子鐘的嵌入式系統(tǒng)中實現(xiàn),同時該系統(tǒng)也需要控制原子鐘的工作,這導(dǎo)致神經(jīng)網(wǎng)絡(luò)可以使用的計算資源會比較受限。 因此,本文采用由WNN改進(jìn)而來的LLWNN神經(jīng)網(wǎng)絡(luò)[13],這種神經(jīng)網(wǎng)絡(luò)只需要一層隱藏層即可實現(xiàn)比較高效的計算。如圖1(b)所示,LLWNN神經(jīng)網(wǎng)絡(luò)的隱藏層由感知器和一類新的神經(jīng)元(小波神經(jīng)元)組成,每個小波神經(jīng)元都和一個感知器對應(yīng)。小波神經(jīng)元對每個輸入元素都進(jìn)行小波變換,并將變換的結(jié)果相乘作為多元小波函數(shù)Ψj(x1,x2,…,xm)輸出:
(3)
式(3)中,
(4)
式(4)中,A={aj,i},D={dj,i}分別是小波函數(shù)的伸縮因子和平移因子。每個小波神經(jīng)元有一個感知器與其對應(yīng),感知器的激活函數(shù)均為f(x)=x,它們的輸出vj作為對應(yīng)小波神經(jīng)元輸出值的權(quán)重。設(shè)神經(jīng)網(wǎng)絡(luò)中小波神經(jīng)元的數(shù)量為L,則最終輸出的預(yù)測值為
(5)
為了應(yīng)用神經(jīng)網(wǎng)絡(luò)對誤差信號演化的過程進(jìn)行模擬,需要對神經(jīng)網(wǎng)絡(luò)的參數(shù)進(jìn)行訓(xùn)練使建立起來的人工神經(jīng)網(wǎng)絡(luò)模型接近真實的演化機(jī)制。在LLWNN神經(jīng)網(wǎng)絡(luò)中,需要訓(xùn)練的參數(shù)包括小波函數(shù)的伸縮和平移因子A,D以及感知器的輸入權(quán)重與偏置W={wj,i},b={bj}。訓(xùn)練過程采用N組參考信號作為訓(xùn)練集{xn,yR,n},使神經(jīng)網(wǎng)絡(luò)在輸入{xn}時輸出的預(yù)報信號{yn}盡量接近{yR,n}。為此,將訓(xùn)練的評價函數(shù)設(shè)定為
(6)
上述訓(xùn)練過程實際是在求解無約束的最優(yōu)化問題,這種問題有兩類比較成熟的解法:一類是梯度下降法及衍生的各種改進(jìn)算法,另一類是包括遺傳算法、粒子群優(yōu)化算法、退火算法在內(nèi)的啟發(fā)式算法。梯度下降法的普適性較好,但是在處理高維問題時很容易受到局域極值問題的阻礙;各種啟發(fā)式算法可以獲得問題的全局最優(yōu)解,但是需要根據(jù)特定的問題類型來選擇合適的算法并設(shè)定恰當(dāng)?shù)膮?shù)。本文采用動量梯度下降算法(momentum gradient decent,MGD)[14]對LLWNN神經(jīng)網(wǎng)絡(luò)進(jìn)行迭代訓(xùn)練,第t+1次訓(xùn)練后的網(wǎng)絡(luò)參數(shù)為
θt+1=θt+vt+1。
(7)
式(7)中:
vt+1=βvt-αθE。
(8)
式(7)和(8)中,參數(shù)β取為動量梯度下降法應(yīng)用中最常用的值0.9[15],α取為訓(xùn)練效果較好時的值0.001。動量vt是上次迭代參數(shù)的改變量,這一項在梯度較小時加速迭代過程,梯度較大時則可以避免振蕩。此外,我們也采用了粒子群優(yōu)化算法[13](particle swarm optimization,PSO)來驗證網(wǎng)絡(luò)的訓(xùn)練效果。
為了驗證LLWNN神經(jīng)網(wǎng)絡(luò)模型的預(yù)報效果,本文利用北京大學(xué)原子鐘小組研制的光抽運(yùn)小銫鐘[8]進(jìn)行實驗,由原子鐘內(nèi)的單片機(jī)每秒采集一次誤差信號,并通過LabVIEW通信程序傳遞給電腦端的Matlab程序。為便于分析,本文使用的誤差信號都按照增益系數(shù)換算成了相對頻率偏差,這一線性變換不影響信號預(yù)測過程。在訓(xùn)練神經(jīng)網(wǎng)絡(luò)前,本文使用擅長預(yù)測數(shù)據(jù)漂移趨勢的卡爾曼濾波算法[3]分析了一段長度為150 ks的誤差信號:先對這一長段信號進(jìn)行百秒平均處理,再將得到的1 500個數(shù)據(jù)依次輸入卡爾曼濾波模型并與實際值比較。由于卡爾曼模型含有協(xié)方差矩陣等未知參數(shù),也使用PSO算法對這些參數(shù)進(jìn)行了訓(xùn)練,使得預(yù)測誤差的均方差最小。圖2是參數(shù)訓(xùn)練至最佳后的效果,卡爾曼模型輸出的預(yù)測值除了開始一小段之外,主要在誤差信號的平均值附近小幅波動;對未經(jīng)過平均處理的信號用卡爾曼模型處理也可以得到類似的結(jié)果。這驗證了前文的分析:誤差信號的主要成份是短期噪聲,長期波動和漂移很弱,更適合通過神經(jīng)網(wǎng)絡(luò)模型來預(yù)測。
圖2 經(jīng)過卡爾曼濾波算法處理的誤差信號 圖3 MGD和PSO訓(xùn)練算法的迭代過程
根據(jù)歷史信號預(yù)測下一個信號的過程實際是對信號變化量即一次差的預(yù)報[7],為了更清晰地描述預(yù)測結(jié)果,此后的訓(xùn)練和預(yù)測過程都是對誤差信號一次差進(jìn)行處理。在數(shù)據(jù)預(yù)處理中,首先得出上述150 ks誤差信號的一次差,此后從這一段數(shù)據(jù)中抽取連續(xù)的406個數(shù)據(jù)點(diǎn),取最前方7個點(diǎn)為第一組數(shù)據(jù),此后依次后移共得到400組數(shù)據(jù),每組的前6個數(shù)據(jù)作為神經(jīng)網(wǎng)絡(luò)的輸入值,第7個數(shù)據(jù)則作為實際參考值與神經(jīng)網(wǎng)絡(luò)的輸出比較。這400組數(shù)據(jù)的前350組用于訓(xùn)練神經(jīng)網(wǎng)絡(luò),訓(xùn)練效果使用神經(jīng)網(wǎng)絡(luò)的輸出值與實際值相比的均方根誤差來衡量。本文在LLWNN神經(jīng)網(wǎng)絡(luò)中使用的小波基函數(shù)是
(9)
訓(xùn)練LLWNN神經(jīng)網(wǎng)絡(luò)時,先隨機(jī)對參數(shù)θ進(jìn)行初始化,此后再應(yīng)用MGD算法或者PSO算法迭代尋找最優(yōu)參數(shù)直到網(wǎng)絡(luò)收斂。圖3所示是神經(jīng)網(wǎng)絡(luò)訓(xùn)練中一次典型的迭代過程,其中MGD算法在100次迭代以內(nèi)就能穩(wěn)定收斂,且訓(xùn)練結(jié)束后的均方根誤差略小于PSO算法的結(jié)果。由于PSO是全局性算法,這表明小波函數(shù)的應(yīng)用確實幫助解決了局域極值的問題,因此梯度下降算法能找到全局最優(yōu)參數(shù);梯度下降算法擅長求解極值問題,所以MGD算法得出的訓(xùn)練結(jié)果略優(yōu)于PSO算法。
神經(jīng)網(wǎng)絡(luò)訓(xùn)練完成之后,利用后50組數(shù)據(jù)測試網(wǎng)絡(luò)模型的預(yù)測能力:將每組前6個數(shù)據(jù)輸入神經(jīng)網(wǎng)絡(luò),將得到的50個預(yù)測值和對應(yīng)的實際參考值作比較。圖4所示是一次典型的預(yù)測結(jié)果,大部分預(yù)測值和實際值比較符合。計算得到圖4中實際值的標(biāo)準(zhǔn)差為1.81×10-11,而預(yù)測誤差的均方根1.12×10-11相比起來有明顯減小。在理想條件下,如果能精準(zhǔn)預(yù)測下一周期的誤差信號,并通過伺服機(jī)制精準(zhǔn)調(diào)節(jié)原子鐘本振的頻率,可以使下一周期的誤差信號降到零;在不完美預(yù)測的條件下,精準(zhǔn)伺服的方式可以根據(jù)預(yù)測值提前將本振頻率調(diào)節(jié)對應(yīng)的值,則下一周期的誤差信號變成原本的真實值減去預(yù)測值,即預(yù)測誤差值。因此預(yù)測誤差的波動限制了利用神經(jīng)網(wǎng)絡(luò)輸出的預(yù)測值提升原子鐘穩(wěn)定度的上限,例如圖4的數(shù)據(jù)表明利用這一輪預(yù)測結(jié)果最多可以將誤差信號的標(biāo)準(zhǔn)差降低到原來的0.62倍,近似于將原子鐘的短期頻率穩(wěn)定度提升到1.6倍。
圖4 LLWNN神經(jīng)網(wǎng)絡(luò)的預(yù)報結(jié)果 圖5 神經(jīng)網(wǎng)絡(luò)預(yù)測性能隨隱藏層節(jié)點(diǎn)數(shù)的變化情況
作為比較的參考對象,我們也在實驗中對BP神經(jīng)網(wǎng)絡(luò)的預(yù)測效果進(jìn)行了測試。由于神經(jīng)網(wǎng)絡(luò)的訓(xùn)練復(fù)雜度會隨著隱藏層數(shù)的增加而迅速增大,我們選擇單層隱藏層的BP神經(jīng)網(wǎng)絡(luò)進(jìn)行測試,從而控制使兩種神經(jīng)網(wǎng)絡(luò)在訓(xùn)練中使用的計算資源不會相差太多。實際測試了使用不同類型激活函數(shù)的效果后,我們使用Tan-Sigmoid函數(shù):
(10)
作為BP神經(jīng)網(wǎng)絡(luò)隱藏層神經(jīng)元的激活函數(shù),使用線性函數(shù)作為輸出層神經(jīng)元的激活函數(shù)。BP神經(jīng)網(wǎng)絡(luò)的訓(xùn)練算法也采用MGD算法。
為了評估神經(jīng)網(wǎng)絡(luò)的預(yù)測效果,我們在上述15萬個數(shù)據(jù)點(diǎn)中隨機(jī)抽取了100個訓(xùn)練-預(yù)測數(shù)據(jù)段,其中每個訓(xùn)練-預(yù)測數(shù)據(jù)段由前述的400組輸入-輸出數(shù)據(jù)組成,在每個數(shù)據(jù)段中都可以對神經(jīng)網(wǎng)絡(luò)模型進(jìn)行一輪訓(xùn)練-測試的處理過程。每個數(shù)據(jù)段中都完成訓(xùn)練-測試操作之后,我們再將這100輪處理中的預(yù)測誤差結(jié)合起來評估,以盡量減小偶有特殊的規(guī)律性數(shù)據(jù)噪聲對神經(jīng)網(wǎng)絡(luò)預(yù)測性能的影響。在每個數(shù)據(jù)段內(nèi)的測試中,神經(jīng)網(wǎng)絡(luò)會輸出50組預(yù)測值,這些預(yù)測值和實際值之間的均方根誤差體現(xiàn)了這一輪的預(yù)測效果。由于不同數(shù)據(jù)段內(nèi)被預(yù)測數(shù)據(jù)本身的噪聲大小互有差異,這些差異也會影響到均方根誤差的值;而后續(xù)分析中需要綜合考慮神經(jīng)網(wǎng)絡(luò)在不同數(shù)據(jù)段中的預(yù)測性能,因此我們在每個數(shù)據(jù)段內(nèi)的一輪測試后先計算預(yù)測值和實際值之間的均方根誤差,再取其與當(dāng)前數(shù)據(jù)段內(nèi)實際值的標(biāo)準(zhǔn)差這二者之比,將結(jié)果作為當(dāng)前數(shù)據(jù)段內(nèi)測試結(jié)果的評價指標(biāo),本文中我們將這一指標(biāo)簡單地稱為“相對預(yù)測誤差”。在完成所有數(shù)據(jù)段中的訓(xùn)練和測試操作后,我們求出這100輪測試的相對預(yù)測誤差的平均值,將其作為評估神經(jīng)網(wǎng)絡(luò)性能的統(tǒng)一指標(biāo),本文中我們將這一指標(biāo)簡單地稱為“相對預(yù)測誤差平均值”。以相對預(yù)測誤差平均值為依據(jù),我們對不同神經(jīng)網(wǎng)絡(luò)的性能做了比較,同時也對單個神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)進(jìn)行了優(yōu)化,將使得相對預(yù)測誤差平均值最小的神經(jīng)元數(shù)目作為最佳神經(jīng)元數(shù)目。如圖5所示,當(dāng)隱藏層的節(jié)點(diǎn)數(shù)增加時,BP神經(jīng)網(wǎng)絡(luò)和LLWNN神經(jīng)網(wǎng)絡(luò)的相對預(yù)測誤差平均值都有先減小后增大的趨勢,這是因為隨著節(jié)點(diǎn)數(shù)增多,神經(jīng)網(wǎng)絡(luò)的可學(xué)習(xí)參數(shù)更多、結(jié)構(gòu)更加復(fù)雜,可以在更大范圍的函數(shù)空間內(nèi)搜索合適的預(yù)測模型;但在同時,搜索空間的增大也意味著搜索復(fù)雜度的提升,這對訓(xùn)練神經(jīng)網(wǎng)絡(luò)時的計算資源和計算精度都提出了挑戰(zhàn)。對于我們實現(xiàn)的BP神經(jīng)網(wǎng)絡(luò),在隱藏層具有8個神經(jīng)元時取得最佳的預(yù)測性能;引入小波神經(jīng)元的LLWNN神經(jīng)網(wǎng)絡(luò)在隱藏層的小波神經(jīng)元數(shù)小于8時都能實現(xiàn)較好的預(yù)測效果,在小波神經(jīng)元數(shù)目過多時,計算復(fù)雜度的增大會導(dǎo)致預(yù)測誤差變大,其中小波神經(jīng)元數(shù)目為6時的預(yù)測誤差最小。
圖6描述了兩種神經(jīng)網(wǎng)絡(luò)各自設(shè)定到最佳的隱藏層節(jié)點(diǎn)數(shù)時,在每個數(shù)據(jù)段測試中的相對預(yù)測誤差。在少數(shù)數(shù)據(jù)段中BP神經(jīng)網(wǎng)絡(luò)的預(yù)測效果也可以接近甚至超過LLWNN神經(jīng)網(wǎng)絡(luò),但是BP神經(jīng)網(wǎng)絡(luò)的預(yù)測結(jié)果不穩(wěn)定;在一些數(shù)據(jù)段內(nèi),BP神經(jīng)網(wǎng)絡(luò)的相對預(yù)測誤差大于1,即預(yù)測誤差的波動比原始數(shù)據(jù)更大,這些情況下若將神經(jīng)網(wǎng)絡(luò)的輸出反饋到原子鐘反而會惡化穩(wěn)定度。在100個數(shù)據(jù)段的測試中,BP神經(jīng)網(wǎng)絡(luò)相對預(yù)測誤差的平均值為0.887 3,最大值為1.029 1;LLWNN神經(jīng)網(wǎng)絡(luò)相對預(yù)測誤差的平均值為0.789 0,最大值為0.936 9。LLWNN神經(jīng)網(wǎng)絡(luò)的預(yù)測效果與BP神經(jīng)網(wǎng)絡(luò)相比有了較大進(jìn)步,但是這離應(yīng)用到原子鐘的要求還有一段距離。如果希望把預(yù)測結(jié)果應(yīng)用于原子鐘內(nèi)參與頻率調(diào)節(jié)的過程,可能需要使用輸出信號能夠反饋回輸入端的循環(huán)神經(jīng)網(wǎng)絡(luò)(recurrent neural network,RNN)進(jìn)一步分析;此外,也還需要對噪聲中晶振頻率波動的成分與其他波動成分做出更清晰的區(qū)分。
圖6 BP與LLWNN神經(jīng)網(wǎng)絡(luò)預(yù)測效果對比圖
本文采用基于LLWNN神經(jīng)網(wǎng)絡(luò)模型的算法對銫原子鐘的內(nèi)部信號進(jìn)行預(yù)報,并應(yīng)用光抽運(yùn)銫原子鐘的誤差信號進(jìn)行了測試和驗證;神經(jīng)網(wǎng)絡(luò)的參數(shù)訓(xùn)練采用了MGD算法,訓(xùn)練結(jié)果表現(xiàn)出了較好的全局性。論文還評估和比較了BP神經(jīng)網(wǎng)絡(luò)和LLWNN神經(jīng)網(wǎng)絡(luò)的預(yù)報性能,在兩種神經(jīng)網(wǎng)絡(luò)都優(yōu)化至最佳結(jié)構(gòu)時,LLWNN神經(jīng)網(wǎng)絡(luò)表現(xiàn)出了更高的預(yù)報精度,且預(yù)報穩(wěn)定性更好。但是目前對誤差信號的預(yù)報精度還不夠高,也無法區(qū)分組成誤差信號噪聲的各類成份及來源;在未來進(jìn)一步的研究中,還可以考慮建立循環(huán)神經(jīng)網(wǎng)絡(luò)來模擬將預(yù)測結(jié)果加入原子鐘頻率控制環(huán)路的過程,建立起合適的反饋模型并尋找最佳的反饋參數(shù)。