摘 要： 隨著數(shù)字技術的日益發(fā)展，應用在天文觀測領域的數(shù)字信號處理技術也越來越成熟。本文針對ADC量化對于中性氫（HI）觀測的影響進行研究分析。首先對HI觀測進行了介紹性描述；在此基礎上，對于使用實驗室4.5 m天線進行HI觀測時的信號接收處理系統(tǒng)進行了描述；重點對ADC量化以及射電望遠鏡靈敏度的相關公式進行推導，并結合4.5 m天線實際觀測的數(shù)據(jù)處理結果進行對比說明，指出并證實了ADC量化位數(shù)不足帶來的對于觀測效果的影響可以通過提高積分時間得到補償。
關鍵詞： ADC；量化；中性氫譜線；靈敏度
0 引言
現(xiàn)階段在很多地方都架設有用于天文觀測的射電天文望遠鏡，架設地點通常會經(jīng)過很多的選址對比，最終選擇出一個相對適合的地點。而為了做測試而搭建的天線考慮到使用方便的因素，有些會選擇在城市工作地點附近搭建。此時城市內(nèi)的一些通信信號就成為了射電強干擾信號，給天文觀測者帶來了很大的影響。本文中4.5 m天線是架設在城市內(nèi)的，針對中性氫譜線觀測來進行討論分析，目的是希望能夠?qū)⒁谏潆姀姼蓴_存在的情況下進行中性氫譜線觀測的觀測者提供一些參考，指出并證實ADC量化位數(shù)不足對于靈敏度的影響可以通過提高積分時間得到補償。
1 中性氫（HI）觀測的介紹
荷蘭天文學家Vande Hulst在1944年提出了銀河系中星系氫基態(tài)超精細結構躍遷譜線（對應氫的21 cm譜線）可以被觀測到的預測。中性氫的觀測開始于20世紀50年代Ewen和Purcell對于銀河系中性氫譜線的觀測，在此之后人們進行了大量的中性氫觀測[1]。
用光學手段研究星際區(qū)域是很困難的，但是這些區(qū)域的氫原子卻可以輻射21 cm譜線，并且不會被那些星際物質(zhì)吸收，因此21 cm譜線成了探測宇宙空間的有力武器。通過觀測可以詳細地了解銀河系中性氫的分布及運動規(guī)律，為精確地測量一些天體的距離奠定基礎。高分辨率的中性氫觀測也是研究星系相互作用的一個重要手段。
處于基態(tài)的中性氫原子，其電子自旋所產(chǎn)生的磁矩相對氫核自旋所產(chǎn)生的核磁矩有兩種可能的取向：平行或反平行。前者的能量高于后者，當處于上能級的中性氫原子躍遷到下能級時，發(fā)出頻率為1 420.406 MHz的光子，產(chǎn)生21 cm氫譜線輻射。在沒有外界擾動情況下，它的自發(fā)躍遷概率很低，一個處于基態(tài)高能級的中性氫原子要在長達約1 100萬年中才有一次躍遷到低能態(tài)的機會，導致了21 cm譜線的自然寬度很小[2]。
2 中性氫觀測的信號接收過程以及觀測設備介紹
觀測使用實驗室現(xiàn)有4.5 m口徑的實驗天線，信號經(jīng)由單極性饋源喇叭接收后經(jīng)過放大、濾波、混頻等一系列模擬信號處理之后傳輸?shù)綄嶒炇液蠖诵盘柼幚韰^(qū)域，將信號輸入到由CASPER設計研制的接收機數(shù)字處理平臺，此時信號中包含很強的通信用信號頻段，對于中性氫觀測來說這些通信寬頻信號就是強度很大的射頻干擾信號。接收到的模擬信號經(jīng)過放大混頻之后首先要輸入ADC板卡來實現(xiàn)數(shù)字化，這是十分重要的一個步驟。經(jīng)過ADC量化之后的信號實現(xiàn)了數(shù)字化，數(shù)字化的信號可以完全按照算法要求來進行后續(xù)的數(shù)字信號處理，從而得到HI譜線的觀測結果。圖1對應實驗室4.5 m天線觀測時信號接收處理過程。
基于CASPE硬件平臺的譜線終端的實現(xiàn)是基于軟件無線電的思路。中心思想是構造一個具有開放性、標準化、模塊化的通用硬件平臺，將各種功能（如多相濾波、頻段獲取、調(diào)制解調(diào)類型、數(shù)據(jù)格式、通信協(xié)議等）通過軟件來實現(xiàn)。從基于硬件、面向用途的實現(xiàn)方法中解放出來。功能的軟件化需要減少功能單一、靈活性差的硬件電路，特別是模擬環(huán)節(jié)，讓數(shù)字化過程盡量靠近天線。采用開放性和可編程性的結構來改變硬件的配置結構[3]。
3 ADC量化對于觀測靈敏度的影響分析[4-8]
ADC板卡是整個過程中十分重要的一環(huán)，它是連接著模擬和數(shù)字的環(huán)節(jié)。圖2描述了ADC的結構圖，前置的濾波器可以把輸入信號中高于某一頻率的頻譜去除，以避免高頻信號混疊到ADC的工作頻帶中。隨后是抽樣和保持電路，它可以對輸入的模擬信號采樣并在一個時間點保持不變，其目的是為了讓后級的電路把這個固定的電平轉(zhuǎn)化為所需要的數(shù)字編碼，這個時間也叫做轉(zhuǎn)換時間。轉(zhuǎn)換的過程也就是量化的過程，在這個階段把基準電平分成多個小的區(qū)間（例如N-bit的ADC就有2N個小區(qū)間）。在轉(zhuǎn)換時間內(nèi)找到與輸入抽樣電平相對應的量化區(qū)間，數(shù)字處理器就能夠把它映射到對應的數(shù)字輸出值，進行編碼輸出，得到最后的數(shù)字輸出信號。
電路特性決定ADC的一些關鍵指標，例如采樣和保持電路決定了ADC的精度和線性；量化器決定了ADC的分辨率；輸出緩沖器決定了輸出速率等。數(shù)模轉(zhuǎn)換器的精度（分辨率）定義為模數(shù)轉(zhuǎn)換器中用以描述一個模擬信號值所采用數(shù)字信號的位寬。N-bit的ADC用2N個不同的數(shù)值代表全量程范圍的連續(xù)模擬信號值，通常用LSB（Least Significant Bit）表示。
下面推導N比特量化帶來的影響。設VF為輸入模擬信號的滿量程電壓值，P為量化電平，若ADC的位數(shù)為N bit，則有關系：
ADC實現(xiàn)的是用離散的數(shù)字量來表示連續(xù)的模擬量，所以這個過程中必然會引入一個因量化而帶入的噪聲。在用離散數(shù)字信號表示模擬量的過程中，因為在量化過程中存在著一個類似四舍五入的過程，因此，引入的量化噪聲的最大量化誤差可以達到P/2。若輸入信號確定，并且假設輸入信號在滿量程的范圍均勻分布，也就是在（-VF/2，VF/2）之間均勻分布，則量化噪聲會在（-P/2，P/2）之間均勻分布，此處設定量化噪聲信號為X，則有：
方差是隨機信號在均值上下波動程度的一種統(tǒng)計特征，是用來說明隨機信號各可能值相對于均值的偏離程度的一種運算，若定義X（t）為平穩(wěn)隨機信號，則有：
D[X（t）]=E[X2（t）]-E2[X（t）]（7）
設定信號X（t）為加在單位電阻上的噪聲電流，則 E[X2（t）]表示消耗在單位電阻上的瞬時功率，包括交流和直流功率的一個平均值。E2[X（t）]可看作消耗在單位電阻上的等效直流功率。從而通過式（7）可知D[X（t）]表示消耗在單位電阻上的瞬時功率的交流部分平均值。
靈敏度對于射電望遠鏡來說是一個很重要的指標，通常是指望遠鏡能觀測到的最小流量，用符號?駐T表示。在以下的推導計算中，Tsys是系統(tǒng)噪聲溫度，t是觀測的積分時間，BW是觀測的帶寬，則有：
當觀測射電源沒有進入望遠鏡觀測范圍時，望遠鏡在觀測帶寬BW內(nèi)的輸入功率如下：
Poff=KTsysBW（9）
對應的功率的波動是：
當觀測射電源進入觀測范圍時，由該射電源引起的接收功率的提升如果高于被觀測源進入觀測范圍前的功率波動Poff，那么對于這個源的探測是一個可靠的探測。用Ptotal代表接收處理的總功率，Psource代表目標觀測源帶來的接收處理功率，而此時有：
當式（14）中的A足夠大（大于3）時，則認為被觀測射電源是一個可以被可靠觀測到的源。
對于經(jīng)過ADC處理的電壓信號，ADC器件的位數(shù)越多，則得到的電壓值越精確。上面已經(jīng)計算了量化誤差引入的噪聲影響，據(jù)此可以得到量化噪聲引入功率Pn：
其中，N為ADC量化位數(shù)。可以看出f（N）=2×f（N+1），當量化位數(shù)足夠多時，可以知道Pn十分小，從而Ptotal幾乎不受影響。也就是說在這種情況下當ADC量化位數(shù)足夠多的時候，量化帶來的影響對于觀測的影響很小，幾乎可以忽略。但是在譜線觀測的進行中，如果在觀測頻帶內(nèi)存在強的RFI（Radio Frequency Interference），情況將會發(fā)生很大的改變。在這種情況下，Tsys可以看成是隨頻率變化的函數(shù)。此時，因為ADC量化誤差而引入的Pn帶來的影響則發(fā)生了很大的變化，Ptotal因為強RFI的存在而發(fā)生了大幅度的提高，用P2表示有RFI的情況下的功率，用P1表示無RFI的情況下的功率，則有：
P1，total=P1，off+P1，source+P1，n（17）
P2，total=PRFI+P2，off+P2，source+P2，n（18）
當把整個譜線觀測壓縮到ADC的量化電壓范圍內(nèi)時，各個量的變化如下：
P2，total>P1，total；P2，off<P1，off；P2，source<P1，source；P2，n=P1，n
對比有RFI存在的觀測情況與沒有RFI存在的情況，則有：
也就是說ADC量化噪聲對于觀測的影響變大了，在射電源信號十分微弱的情況下，這種變化可能會導致信號被淹沒在噪聲中。根據(jù)式（13）可知，提高積分時間可以提高系統(tǒng)的靈敏度而得到更好的觀測效果，而ADC量化帶來的影響是式（13）中Ptotal中的一部分，所以如果在有RFI的情況下進行觀測，ADC量化帶來的影響的加劇可以通過增加積分時間來進行優(yōu)化。
4 觀測分析證明
在中性氫信號比較強的時段使用4.5 m天線進行觀測，ADC為8 bit，有效數(shù)據(jù)部分占用為4 bit，將數(shù)據(jù)進行32K個點的FFT，經(jīng)過100 000次累加后得到如圖3中虛線所示的譜線圖，將數(shù)據(jù)進行1 bit量化得到新的數(shù)據(jù)進行同樣的處理得到的譜線圖如圖3中實線圖示，對比可以看出數(shù)據(jù)量化位數(shù)的降低，導致了信噪比明顯的降低。圖4是1 bit量化對應譜線數(shù)據(jù)與4 bit量化對應的譜線數(shù)據(jù)的比值，可以看出比值在對應中性氫信號頻段發(fā)生了很明顯的強于其余頻段的下降，說明較低的量化位數(shù)的量化噪聲帶來的影響更為明顯。
將4 bit量化的數(shù)據(jù)進行32K個點FFT，經(jīng)過了1 000次疊加得到譜線如圖5圖示，100 000次疊加得到的譜線如圖6圖示，實際的數(shù)據(jù)處理結果證實了增加積分時間提高了信噪比，從而可以補償?shù)土炕粩?shù)帶來的對于觀測的影響。
5 結論
在實驗室進行調(diào)試時，窄帶強干擾信號如果不經(jīng)過合適的衰減，則會使ADC飽和，發(fā)生截斷輸出，從而使頻譜輸出波形發(fā)生嚴重的扭曲和變形。此時，必須增大對于輸入觀測信號的衰減，保證強干擾信號不飽和ADC，才能觀測到正常的譜線。但是增加對于輸入觀測信號的衰減使強干擾信號不飽和ADC后，在HI觀測頻段，ADC器件本身器件的噪聲就凸顯出來，所要觀測的HI信號被混在噪聲之中。統(tǒng)計意義上說可以通過增加足夠多的采樣點來進行補償，也可以理解為就是相應的提高積分時間來增加系統(tǒng)的靈敏度，這樣仍然可以得到有效的觀測結果。
參考文獻
[1] EWEN H I， PURCELL E M. Observation of a line in the galactic radio spectrum： radiation from galactic hydrogen at 1，420 Mc./sec[J]. Nature， 1951，168：356.
[2] 錢磊.中性氫研究的前言問題[J].天文學進展，2011，29（4）：408-421.
[3] 楊小牛，樓才義，徐建良.軟件無線電原理與應用[M].北京：電子工業(yè)出版社，2001.
[4] FELLI M， SPENCER R E. Very long baseline interferometry[M]. Netherlands： Kluwer Academic Publishers，1989.
[5] 周濤.量化精度對數(shù)字化接收機性能的影響[J].電子信息對抗技術，2006，21（6）：20-23.
[6] 李迪.高新能sigma-delta_ADC的設計與研究[D].西安：西安電子科技大學，2010.
[7] 史軍，張乃通，王輝宇.UWB數(shù)字接收機量化階數(shù)的研究[J].無線電工程，2008，38（2）：20-22.
[8] 姜碧溈.射電天文工具[M].北京：北京師范大學出版社，2008.