李 嬋，王浩宇，繆海興，韓 森

（上海理工大學 光電信息與計算機工程學院，上海 200093）

引言

Pound-Drever-Hall（PDH）技術是一項用于改善現(xiàn)有激光器頻率穩(wěn)定性的技術，它利用反射光束的相位信息產(chǎn)生誤差信號，通常用于激光穩(wěn)頻。它可以通過將光腔鎖定在頻率穩(wěn)定的激光器上來測量腔長的微小的變化[1-2]，這對引力波探測很重要。引力波探測器光學參考腔的精度一般都很高，其精度要達到10?21m。在實際實驗室中（如LIGO）存在地震噪聲、電子噪聲等環(huán)境擾動，入射至光學參考腔的激光頻率會使腔長發(fā)生微小的變化[3]，進而光腔的共振頻率會發(fā)生抖動，使得PDH 穩(wěn)頻系統(tǒng)鎖定之后激光的頻率不穩(wěn)定。因此，PDH 技術不但可以來測量光腔的長度微小變化，也可以用來穩(wěn)定激光器的頻率[1-2]。

PDH 技術由于其動態(tài)范圍較小，控制回路只有在偏離工作點很小的范圍內(nèi)才能穩(wěn)定，當受環(huán)境震動影響過大時，誤差信號會變得非線性，PDH 穩(wěn)頻系統(tǒng)將失去鎖定。針對上述缺點，有不少團隊進行了研究[4-5]，如Hassen 等演示了使用線性二次高斯控制（LQG）法來鎖定光腔[6-7]；Schütte 等對此進行實驗證明，證明了LQG 方法的穩(wěn)健性，這在一定程度上拓寬了PDH 技術的動態(tài)范圍，有利于PDH 穩(wěn)頻系統(tǒng)的鎖定捕獲，但算法較為繁瑣[8]；Miyoki 等提出了執(zhí)行近Q 相位解調(diào)的解決方案，但降低了PDH誤差信號的靈敏度[9]。

本文以PDH 技術為基礎，建立了一個簡單直觀、抗干擾強的反饋控制系統(tǒng)，新定義了PDH 誤差信號方程，并結合卡爾曼濾波器估計了反射鏡的狀態(tài)。相比于以往的工作，該方法的模型更為簡單直觀。通過仿真可知，本文的方法可以使PDH 穩(wěn)頻系統(tǒng)在更寬的動態(tài)范圍內(nèi)更快實現(xiàn)鎖定，這會對改進引力波探測器中鎖相反饋控制提供新的啟發(fā)。

1 傳統(tǒng)PDH 穩(wěn)頻系統(tǒng)的基本原理

傳統(tǒng)的PDH 技術的裝置圖如圖1 所示。該裝置是通過將光腔鎖定到一個超穩(wěn)的激光器上實現(xiàn)穩(wěn)長。工作原理是：以光腔的共振頻率為基準頻率，對激光頻率進行相位調(diào)制后，可以在激光頻率兩側產(chǎn)生對稱分布在激光載波兩側、幅度相等但相位相反的兩個邊帶[10]；將這三個不同頻率的光入射到光腔，用光電探測器（PD）可以分別檢測光腔的反射功率Pref和透射功率 Ptrans；將反射功率與信號發(fā)生器發(fā)出的正弦信號經(jīng)混頻器混頻，提取低頻分量，得到一個低頻信號，一般稱其為誤差信號；該誤差信號在一定范圍內(nèi)與反射鏡的位移成正比，能用該信號去負反饋，從而補償反射鏡的位移。

2 新PDH 穩(wěn)頻系統(tǒng)的建模

2.1 新PDH 穩(wěn)頻系統(tǒng)

對于傳統(tǒng)的PDH 技術，一般檢測在光腔出射端進行。本文中，我們還考慮了在光腔入射端處對透射功率Ptrans的檢測，并增加了一個狀態(tài)估計器，以估計出射端反射鏡的狀態(tài)，將其反饋到光腔腔體。圖2 顯示了新PDH 穩(wěn)頻系統(tǒng)的基本框架，主要由機械和探測子系統(tǒng)、狀態(tài)估計器和伺服反饋組成。

圖2 新PDH 穩(wěn)頻系統(tǒng)框圖Fig.2 Layout of new PDH frequency stabilization system

此外結合實際物理意義，本文所涉及的各項通用的參數(shù)設計如表1 所示。

表1 新PDH 穩(wěn)頻系統(tǒng)的參數(shù)Tab.1 Parameters of new PDH frequency stabilization system

2.2 機械和探測子系統(tǒng)的設計

我們主要需要控制的是反射鏡的位移。為了對整個光腔進行建模，根據(jù)分離原理[11]，我們將所提出的新PDH 穩(wěn)頻系統(tǒng)分為兩個子系統(tǒng)：機械子系統(tǒng)和探測子系統(tǒng)。機械子系統(tǒng)包括反射鏡及控制其位移的壓電陶瓷（PZT）。探測子系統(tǒng)包括光腔腔體和適當?shù)奶綔y器。機械和探測子系統(tǒng)裝置如圖3 所示。其中，噪聲ng表示地面運動引起的反射鏡位移噪聲，nPDH和nPtrans分別為誤差信號和透射端的電路噪聲。

圖3 機械和探測子系統(tǒng)Fig.3 The mechanical and detection subsystem

2.2.1 機械子系統(tǒng)

對于PDH 穩(wěn)頻系統(tǒng)來說，反射鏡的位移通常使用壓電陶瓷來控制。這里我們用一個四階諧振系統(tǒng)模型來模擬用于驅動反射鏡之一的壓電陶瓷。其傳遞函數(shù)為

式中：G(s)為四階諧振系統(tǒng)的輸出量與輸入量經(jīng)拉普拉斯變換后的函數(shù)；s為復參數(shù)。該模型中包含了兩個約為15 kHz 和30 kHz 的高階諧振頻率。

2.2.2 探測子系統(tǒng)

如圖1 所示，來自激光器的光束由電光調(diào)制器（EOM）調(diào)制。電光調(diào)制器調(diào)制產(chǎn)生兩個頻率為 (ω0±Ω)/2π的一階邊帶，其中 ω0/2π 是載波頻率，Ω/2π 是調(diào)制頻率。當相位調(diào)制邊帶通過諧振腔時，誤差信號與載波諧振產(chǎn)生的信號具有相反的符號。此時，可以從載波和邊帶之間的拍頻信號中提?。ㄅc反射鏡位移成正比的）出相位信息。反射功率PDrefl的同相解調(diào)導致以下（c os 正交）誤差信號，即

式中：P0為載波功率；Ps為邊帶功率；Im 為虛部；ω0為載波角頻率；Ω為調(diào)制角頻率；rc為腔反射率。rc可表示為

式中：r1、r2、t1分別為入射端反射鏡和出射端反射鏡的振幅反射率以及出射端反射鏡的透射率；L為腔體長度（等于波長的整數(shù)）；ω 為激光角頻率；x為用微觀偏移量表示的鏡面位移；c為光速 3×108mm/s 。

光腔透射端處測量的透射功率為

式中：ntrans是噪聲項；tc為腔透射率。tc可表示為

式中：t2為出射端反射鏡的振幅透射率。

2.3 新定義的PDH 誤差信號

在實際PDH 系統(tǒng)中，通常以光腔腔體的反射功率作為判斷激光頻率是否達到共振的方向。光腔的諧振發(fā)生時，反射功率達到最小值。傳統(tǒng)的PDH 系統(tǒng)的線性動態(tài)范圍很窄，光腔附近的PDH 誤差信號變化很快，這給PDH 穩(wěn)頻系統(tǒng)的鎖定捕獲帶來了很大的挑戰(zhàn)。為此，本文提出了一種擴大線性動態(tài)范圍的方法，將透射功率Ptrans除以傳統(tǒng)的PDH 誤差信號 εPDH以獲得一個新PDH誤差信號，即

圖4 所示的是傳統(tǒng)PDH 誤差信號與新誤差信號隨反射鏡偏移量的變化的比較以及在不同反射鏡偏移量下的透射功率。圖中可以看出新方法不僅擴展了PDH 穩(wěn)頻系統(tǒng)的線性范圍，也改善了系統(tǒng)的高度非線性情況。后者在實際應用中，尤其在存在噪聲和大干擾的環(huán)境下尤具優(yōu)勢。

圖4 傳統(tǒng)PDH 誤差信號與新誤差信號的比較以及透射功率Fig.4 Comparison between the traditional PDH error signal and new PDH error signal, and transmitted power

2.4 反射鏡的狀態(tài)估計

在新PDH 穩(wěn)頻系統(tǒng)中，由于除法，當透射端光強接近于零的時候，其電路噪聲被顯著放大，使系統(tǒng)十分容易受到噪聲影響，因此減小噪聲的影響十分必要。一個被廣泛使用的適用于從有噪聲的測量中進行估計的工具是卡爾曼濾波器。卡爾曼濾波器的算法主要分為兩步：預測和更新[12-15]。此外，使用卡爾曼濾波器作為狀態(tài)估計器，也可增強系統(tǒng)的抗干擾能力。

2.5 噪聲時域模型

傳統(tǒng)PDH 誤差信號中的散粒噪聲可表示為

式中：Ps為邊帶功率；? 為約化普朗克常數(shù)。

在本文的模型中，主要考慮了兩個主要部分的噪聲：過程噪聲和測量噪聲。其中，過程噪聲主要來自地面運動，測量噪聲包括常規(guī)誤差信號中的散粒噪聲、透射功率中的散粒噪聲。我們將測量噪聲與光波長進行歸一化，則VnPDH≈10?7,VnPtrans≈10?2。

2.6 控制反饋系統(tǒng)設計

對于傳統(tǒng)的PDH 系統(tǒng)，若系統(tǒng)受到強擾動使其遠離線性工作區(qū)，鎖定獲取會完全失控。因此，控制反饋系統(tǒng)的設計尤為重要。這里，我們考慮了基于LQG 的一種改進的方法：利用卡爾曼濾波器估計光腔反射鏡的狀態(tài)，并尋求合適的控制律。具體是PDH 誤差信號通過狀態(tài)估計器再經(jīng)反饋放大器，隨后反饋給出射端反射鏡，從而確定了新PDH 系統(tǒng)的模型。

3 仿真結果

在引力波探測器中，PDH 技術主要用于降低100 Hz 以內(nèi)的地震噪聲、環(huán)境噪聲以及空氣擾動帶來的諧振腔長度變化。因此，本模型主要是在100 Hz 以內(nèi)的低頻線性范圍內(nèi)進行研究。圖5 為基于新定義的PDH 誤差信號的穩(wěn)頻系統(tǒng)的Simulink 仿真模型，主要包括：用四階諧振系統(tǒng)表示的PZT 模型、基于新PDH 誤差信號觀測光腔反射鏡狀態(tài)的光腔模型和一個恒定大小為20 的負反饋增益。其中，光腔的參數(shù)如表1所示，其頻率響應模型可以看成是一個低通濾波器，截止頻率約為240 kHz，在0～100 Hz范圍內(nèi)的增益等于PDH 誤差信號在零點位置的斜率，大小恒定。另外，PZT 在100 Hz 以內(nèi)的頻率響應也可近似認為大小恒定，故整個PDH 穩(wěn)頻系統(tǒng)回路的頻率特性主要由伺服反饋系統(tǒng)決定。在實際中，伺服反饋系統(tǒng)的頻率特性近似為一個一階低通濾波器，其截止頻率一般為1～10 kHz。

圖5 新PDH 穩(wěn)頻系統(tǒng)的Simulink 模型Fig.5 The Simulink model of the new PDH system

為了驗證新的PDH 誤差信號的穩(wěn)定性，我們在時間t=0 時，施加一個沖擊力，以使反射鏡在x=0 時遠離共振。本文的仿真主要分為兩個方案：（1）對傳統(tǒng)PDH 誤差信號使用卡爾曼濾波；（2）對新的PDH 誤差信號使用卡爾曼濾波。

圖6 和7 分別顯示了在閉環(huán)控制時，傳統(tǒng)PDH 誤差信號經(jīng)卡爾曼濾波和新PDH 誤差信號分別經(jīng)卡爾曼濾波來估計出射端反射鏡位置的兩種方案的模擬結果，同時也分別展示了對應的反饋信號。比較圖6（a）與圖7（a），可以看出用傳統(tǒng)的PDH 誤差信號進行卡爾曼濾波器估計時，系統(tǒng)只能在很窄的線性范圍內(nèi)準確地預測反射鏡的位置，其動態(tài)范圍僅為0.002 λ 。而使用基于新PDH 誤差信號經(jīng)卡爾曼濾波器估計時，可以將反射鏡估計位置的擴展到±0.1 λ 的范圍，系統(tǒng)的動態(tài)范圍擴大了約50 倍。

圖6 基于傳統(tǒng)PDH 誤差信號與卡爾曼濾波器的仿真波形Fig.6 Simulation waveform based on traditional PDH error signal and Kalman filter

此外，我們針對兩種不同方案設置了不同的過程噪聲（來自地面噪聲）和測量噪聲值（來自測量PDH 和透射功率的測量噪聲），并討論了不同模型的噪聲容限。噪聲容限的大小取決于能使得PDH 穩(wěn)頻系統(tǒng)穩(wěn)定的極限情況。表2 總結了兩種方案的動態(tài)范圍，以及其過程噪聲Qk和測量噪聲Rk的噪聲容限。

從圖7 中可以看出，當PDH 穩(wěn)頻系統(tǒng)重新鎖定之后，反饋信號的波動幅度超過10%，這是由于我們在模型中引入了很高的地震噪聲。因此，采用新PDH 控制方案的噪聲容限Qk遠高于傳統(tǒng)PDH 方案。對于足夠高的信號頻率，散粒噪聲是主要的噪聲來源。散粒噪聲的功率譜密度為，其中h為普朗克常數(shù)，Ps、λ的值參考表1。再將Pe做歸一化處理，除以一個因數(shù) λ2，可以得到歸一化的Pe值約為 3×10?8。該噪聲容限的大小符合表2 中Rk的數(shù)量級設定，參考該值，可以驗證本文的設定值符合實際系統(tǒng)標準。

表2 不同方案的動態(tài)范圍以及噪聲容限Tab.2 The broadened dynamic range and noise tolerance of different approaches

圖7 基于新PDH 誤差信號與卡爾曼濾波器的仿真波形Fig.7 Simulation waveform based on new PDH error signal and Kalman filter

4 結論

為了解決傳統(tǒng)PDH 穩(wěn)頻系統(tǒng)線性動態(tài)范圍過小和噪聲容限的問題，本文提出了一種拓展PDH 穩(wěn)頻系統(tǒng)的動態(tài)范圍的方法。將傳統(tǒng)的PDH 誤差信號除以透射功率得到新的PDH 誤差信號以增加線性動態(tài)范圍?；趥鹘y(tǒng)PDH 穩(wěn)頻系統(tǒng)、卡爾曼濾波器以及控制技術的基本原理，通過對比傳統(tǒng)PDH 系統(tǒng)和新的PDH 系統(tǒng)的特點，我們搭建了新的PDH 穩(wěn)頻系統(tǒng)的基本模型。結合卡爾曼濾波器對未知變量的估算能力，對PDH 穩(wěn)頻系統(tǒng)使用卡爾曼濾波器作為狀態(tài)估計器，來估計光腔反射鏡的狀態(tài)，包括光腔反射鏡的位置和動量。并將卡爾曼濾波分別應用于傳統(tǒng)的PDH 誤差信號和新的PDH 誤差信號中，從而可以說明新的PDH 誤差信號的方法可以顯著增加PDH 穩(wěn)頻系統(tǒng)的線性動態(tài)范圍。此外，這種方法可以使得PDH 穩(wěn)頻系統(tǒng)的動態(tài)范圍擴大50 倍，并在系統(tǒng)受到較大擾動而失去鎖定時，實現(xiàn)系統(tǒng)自動鎖定，同時也提高了新的PDH 誤差信號的噪聲容限。本文的研究主要在線性控制部分，針對非線性估計器的選擇以及非線性系統(tǒng)的控制也將是未來研究的一個重點。