馬　勇

(民航西南空管局 氣象部,　成都 610225)

?

機場天氣雷達雙通道數(shù)字中頻處理系統(tǒng)的實現(xiàn)

馬勇

(民航西南空管局 氣象部,成都 610225)

機場天氣雷達要求能夠從復(fù)雜的天氣環(huán)境中識別不同的天氣狀況以保障航空飛行安全，其接收機大動態(tài)以及抗噪性能設(shè)計對雷達至關(guān)重要。在分析模數(shù)轉(zhuǎn)換器(ADC)對雷達中頻接收機動態(tài)范圍制約的基礎(chǔ)上，根據(jù)中頻帶通采樣和數(shù)字下變頻的原理，實現(xiàn)了基于現(xiàn)場可編程門陣列的雙通道ADC采樣數(shù)字中頻處理系統(tǒng)，并給出了系統(tǒng)的設(shè)計原理、方法以及測試結(jié)果。通過對雙通道ADC采樣的數(shù)字中頻處理系統(tǒng)的實現(xiàn)，能夠很好地提高天氣雷達接收機的動態(tài)范圍，并應(yīng)用于機場多普勒天氣雷達數(shù)字中頻接收機。

機場天氣雷達；數(shù)字中頻；數(shù)字下變頻；數(shù)字濾波;動態(tài)范圍

0　引　言

數(shù)字中頻接收技術(shù)是近年來迅速發(fā)展的數(shù)字技術(shù)之一，具有瞬時動態(tài)大、失真小、幅相一致性和I/Q正交性高等特點[1-3]，該技術(shù)已經(jīng)廣泛應(yīng)用到機場天氣雷達系統(tǒng)之中。機場天氣雷達的目的是保障終端區(qū)航空飛行氣象安全，因此其監(jiān)測的氣象目標(biāo)較為廣泛，從晴空風(fēng)切變到強降水、雷暴均需要準(zhǔn)確探測，而且監(jiān)測距離范圍達到200 km以上。這些氣象目標(biāo)具有反射能力差別懸殊、變化迅速的特點，其回波幅度變化范圍較大。一般來說，各類氣象目標(biāo)的變化范圍超過100 dB，所以大動態(tài)范圍成為機場天氣雷達設(shè)計追求的目標(biāo)。由于機場天氣雷達需要準(zhǔn)確測量以利于反演降水量及判斷危及航空安全的天氣現(xiàn)象。因此，要求在大動態(tài)范圍中保持較高的線性度。過去雷達采用模擬接收機，中頻則采用自動增益控制技術(shù)來滿足大動態(tài)范圍的要求。近年來，數(shù)字中頻接收機的應(yīng)用提高了雷達性能，簡化了雷達結(jié)構(gòu)，同時，為實現(xiàn)高線性度大動態(tài)范圍提供了更多的方法。

當(dāng)前國內(nèi)機場多普勒天氣雷達普遍采用單通道模數(shù)轉(zhuǎn)換器(ADC)采樣進行數(shù)字化處理，這種方式容易導(dǎo)致系統(tǒng)動態(tài)范圍受限于ADC動態(tài)[4]。本文采用雙通道數(shù)字中頻處理技術(shù)，可將動態(tài)范圍擴展超過100 dB，滿足機場天氣雷達的使用要求。

1　模數(shù)轉(zhuǎn)換器及動態(tài)范圍

采用數(shù)字化中頻接收技術(shù)首先要進行模數(shù)(AD)采樣，雷達的動態(tài)范圍往往受限于ADC的動態(tài)[4]。動態(tài)范圍為最強信號與最弱可檢測信號的功率之比，最強信號功率往往采用ADC輸入的飽和值，最弱信號一般采用靈敏度的功率值?？梢?，對于ADC動態(tài)(最大信噪比)而言，影響其大小的因素主要有兩個：一是ADC的飽和功率大小，二是ADC的量化噪聲功率大小。這就意味著如果回波功率超過ADC的飽和功率，雷達系統(tǒng)就會發(fā)生過載現(xiàn)象，造成目標(biāo)失真；而如果回波功率低于ADC的量化噪聲功率，信號則會淹沒在噪聲當(dāng)中，必須通過數(shù)字信號處理才能把信號從噪聲中提取出來。

ADC的最大輸入范圍通常定義為振幅與ADC的最高電平相匹配的正弦波。如果信號比這個最大電平還大，則輸出波形將被限幅。如果一個輸入信號比該信號小，則不是所有的比特位都能被置位。最大電平通常決定動態(tài)范圍的上限。如果沒有噪聲且輸入電壓與ADC的最大范圍匹配，則最大電壓Vmax為

Vmax=2b-1×q

(1)

式中：b為AD位數(shù)；q為每一量化電平的電壓。式(1)表明正弦波可以達到最大量化電平的頂部和最小量化電平的底部。

幅度與最大電壓相匹配的正弦波的功率為

(2)

式(2)中，輸入阻抗假設(shè)為單位阻抗，在實際計算中需要包含實際阻抗以獲取實際的功率，但不影響我們的分析。

ADC將模擬信號轉(zhuǎn)化成數(shù)字信號的過程是一個非線性過程，在正弦波的真實值和量化值之間存在一個差值，這個差值隨著ADC分辨率(位數(shù))的提高而不斷減小，但一直存在，這就是量化噪聲。研究表明量化噪聲在量化電平q上是均勻分布的，概率密度為1/q，這樣就可以求得量化噪聲功率

(3)

這就是理想ADC的靈敏度電平，將之代入動態(tài)范圍的公式

10lg1.5+20blg2=1.76+6.02b

(4)

采樣后進行信號抽取，則數(shù)字信號處理得益為

D=10lg(fs/2B)

(5)

式中：fs為ADC的采樣時鐘；B為信號帶寬。

由式(5)可見，b和fs兩個指標(biāo)都會對接收機動態(tài)D產(chǎn)生影響。ADC芯片的采樣頻率fs每增加1倍，接收機瞬時動態(tài)D改善3 dB，分辨位數(shù)每增加1位，接收機瞬時動態(tài)增加6 dB。14 bit的ADC動態(tài)范圍理想情況下也僅有84 dB左右。實際上，AD的信噪比還受到其他因素的影響，提高采樣率則時鐘抖動對AD的影響變大，當(dāng)AD位數(shù)提高，導(dǎo)致最小量化電平降低對參考電源要求更高。當(dāng)工藝水平達不到相當(dāng)水平非線性失真和線性失真將大大限制AD的有效信噪比。14 bit的AD變換器目前達到的有效信噪比在76 dB左右，遠低于84 dB的理想水平，這對于100 dB的天氣變化范圍要求，很難完全滿足需求。

2　系統(tǒng)原理和設(shè)計構(gòu)架

本設(shè)計中系統(tǒng)動態(tài)范圍擴展的核心是模擬中頻信號經(jīng)過功分器后變成兩路中頻模擬信號，它們分別進行功率放大和功率衰減后送入ADC中進行數(shù)字化處理，之后對下變頻輸出的從屬于每一路的I/Q信號通過算法進行選擇鏈接，輸出一路I/Q信號。下面結(jié)合天氣雷達接收機對系統(tǒng)原理進行分析。

天氣雷達接收機系統(tǒng)原理圖如圖1所示[5]。雷達接收機的動態(tài)范圍一般是由圖1中的低噪放大器(LNA)限定的(一般大于100 dB)，現(xiàn)有工藝技術(shù)下模擬器件對系統(tǒng)動態(tài)影響很小，可以忽略不計。對于單路中頻采樣的動態(tài)，如果采用14 bit量化的ADC，動態(tài)范圍大約可以達到82dB(1MHz帶寬)，則靈敏度為+6 dBm(最大輸入信號)-82 dBm = -76 dBm。

圖1　雷達接收機框圖

本文設(shè)計方案原理是采用兩個ADC進行采樣，其中，一個ADC處理高端信號，另一個ADC處理低端信號。14 bit的ADC動態(tài)為82 dB，那么兩個ADC協(xié)作的動態(tài)理論上就可以達到2 dB×82 dB=164 dB。因此，充分利用這兩個ADC，使得所需要的動態(tài)范圍通過放大和衰減都能在其中某一個ADC中非飽和情況下得到檢測。

假設(shè)接收系統(tǒng)輸入端等效噪聲就是熱噪聲，模擬前端的噪聲系數(shù)是1 dB，這樣模擬前端的噪聲基底就是

Nn=kTNA=-114dBm/MHz+1dB=

-113dBm/MHz

(6)

式中：k為玻爾茨曼常數(shù)，k=1.38×10-23J/K；T為290 K；NA為模擬前端噪聲系數(shù)。

再假設(shè)模擬前端的增益為G，這樣我們就可以計算出數(shù)字中頻輸入端的噪聲為-113 dBm/MHz+G。為了獲得最佳的動態(tài)范圍和靈敏度需要折中考慮，也就是確定前端增益G。我們用Ni和Nd分別表示輸入到數(shù)字中頻的噪聲功率和數(shù)字中頻自身的噪聲功率。需要說明的是為了使前后的帶寬一樣我們這里討論的均以1 MHz帶寬為基準(zhǔn)，所以這里所說的噪聲功率可以理解為噪聲密度，實際使用中以系統(tǒng)的最終帶寬為準(zhǔn)。這樣就可以計算靈敏度損失和動態(tài)范圍損失，分別記為ΔS和ΔD

(7)

這兩個特性都取決于輸入到數(shù)字中頻的噪聲功率和數(shù)字中頻自身的噪聲功率的比值，即

(8)

從式(7)可以看出接收系統(tǒng)的動態(tài)范圍極限值就是數(shù)字中頻的動態(tài)范圍，同時，靈敏度也將損失殆盡。當(dāng)比值R=1時(輸入噪聲和自身噪聲相同時)，靈敏度和動態(tài)范圍同時損失3 dB，這與傳統(tǒng)計算一致。在靈敏度與動態(tài)范圍的取舍中建議在保證靈敏度的情況下盡量減少動態(tài)范圍的損失。在這里取1 dB的靈敏度損失，動態(tài)范圍將損失6.9 dB，R=5.9 dB，系統(tǒng)的增益G的取值如下(Nd=-76 dBm/MHz)

G=-76dBm/MHz-(-113dBm/MHz)+R=

42.9dB

(9)

假設(shè)LNA的噪聲系數(shù)為1 dB、動態(tài)范圍為105 dB(-113 dB～-8 dB)，這樣輸入到AD端口的信號將達到-70 dBm～+35 dBm，這在小功率放大器中是不可能的，也是不現(xiàn)實的。在實際電路中可以使用不同的增益配置分成兩段：一段采用43 dB的增益，保證小信號能在-70 dBm～+10 dBm(10 dBm為放大器的飽和電平)；一段采用14 dBm的增益，保證大信號不飽和，輸入信號為-94 dBm～+6 dBm。這樣就能保證104 dB動態(tài)范圍。

為驗證系統(tǒng)作如下假設(shè)：

(1)模擬中頻信號經(jīng)過功率分配器產(chǎn)生兩路模擬中頻信號，則每一路中頻信號功率為原來的1/2，兩路模擬中頻信號動態(tài)范圍均為76 dB(-70 dBm～+10 dBm)。

(2)第1路信號加入衰減器，衰減大小4 dB，則第1路模擬中頻信號輸出S1動態(tài)范圍為76 dB(-80 dBm～+6 dBm)，模擬大信號；第2路不處理，模擬信號限幅至10 dBm，第2個通道信號功率為-70 dBm～+10 dBm，模擬小信號。

通過以上分析可見，105 dB 動態(tài)范圍的中頻信號經(jīng)過數(shù)字化處理后都會不失真地存在于第1路或第2路信號中(并有交疊)，我們所需要做的工作就是把這兩個通道中的信號通過數(shù)字信號處理的方式分離并組合出不失真的I/Q信號，進而求出回波信號的幅度和相位信息。

系統(tǒng)設(shè)計構(gòu)架如圖2所示。

圖2　系統(tǒng)設(shè)計框圖

3　雷達中頻處理系統(tǒng)實現(xiàn)

本系統(tǒng)的中頻硬件板設(shè)計構(gòu)成，如圖3所示。圖中，系統(tǒng)的硬件組成主要包括Altera StratixII EP2S30F484C5現(xiàn)場可編程門陣列(FPGA)、AD采樣芯片AD6645、時鐘轉(zhuǎn)換芯片AD9510以及電源管理模塊。同時，為了滿足FPGA硬件配置以及信號差分傳輸?shù)男枰到y(tǒng)中還包括EPCS16配置芯片和AM26LS31差分傳輸芯片。

圖3　雷達中頻處理硬件結(jié)構(gòu)

由于硬件板內(nèi)部芯片對電源要求各不相同，如：FPGA內(nèi)部共需要有1.2 V、3.3 V數(shù)字電壓和1.2 V模擬電壓，AD6645需要5 V模擬電壓和3.3 V數(shù)字電壓，并且其他芯片也需要相應(yīng)指標(biāo)電源供電，所以電源管理模塊需要為硬件系統(tǒng)產(chǎn)生各種所需要的電源。為了減少噪聲，硬件中采用了模擬電和數(shù)字電獨立供電的方式，其中，輸入模擬電壓和數(shù)字電壓均為6 V，其他各種所需要電源電壓分別通過LT1765、LT1763、LT196333等電壓轉(zhuǎn)換芯片完成。在硬件設(shè)計中，由于系統(tǒng)既存在模擬電路也存在數(shù)字電路，并且信號數(shù)據(jù)速率很高，為了減少數(shù)字電路對模擬電路的干擾(主要是數(shù)字電路對AD6645的串?dāng)_)，影響模擬電路的性能，這里對模擬電路和數(shù)字電路采用IT715芯片進行了磁隔離。硬件中經(jīng)過FPGA處理的數(shù)據(jù)(如：定時時鐘、I/Q數(shù)據(jù))對外走線均采用差分走線方式，通過差分轉(zhuǎn)換芯片AM26LS31完成信號的差分數(shù)據(jù)轉(zhuǎn)換，差分轉(zhuǎn)換后的I/Q信號以及相應(yīng)時鐘通過PCI總線送數(shù)據(jù)處理板。

由于 FPGA是使用靜態(tài)存儲器單元存儲配置數(shù)據(jù)的，每次掉電后，配置數(shù)據(jù)必須重新下載到FPGA中。設(shè)計中采用串行配置芯片EPCS16SI16N對FPGA進行配置，設(shè)計中預(yù)留了JTAG接口和主動配置接口，以便于硬件調(diào)試和配置。

系統(tǒng)設(shè)計指標(biāo)如下：

輸入模擬中頻信號：60 MHz

信號模擬中頻帶寬：1 MHz

輸出數(shù)據(jù)率： 1 MHz

采樣頻率：48 MHz

動態(tài)范圍：≥100 dB

I/Q正交一致性： ≤0.2°

I/Q幅度一致性： ≤0.1 dB

4　數(shù)據(jù)驗證與測試

本系統(tǒng)測試平臺如圖4所示，ADC采集的數(shù)據(jù)送入數(shù)字中頻處理系統(tǒng)當(dāng)中，當(dāng)系統(tǒng)工作的時候通過QuartusII 內(nèi)嵌的邏輯分析儀SignalTap功能將處理后的I/Q數(shù)據(jù)導(dǎo)入計算機當(dāng)中，在Matlab中進行時域和頻域分析。在采集數(shù)據(jù)過程中SignalTap數(shù)據(jù)長度設(shè)置為256，即每組數(shù)據(jù)采集256個點。

圖4　測試平臺

當(dāng)輸入信號A(59.9 MHz,-40 dBm)時，輸出I/Q信號時域和頻域變換，如圖5所示。

圖5　測試數(shù)據(jù)圖

由圖5可以看出，系統(tǒng)產(chǎn)生大約0.1 MHz的頻偏。為了測試系統(tǒng)動態(tài)范圍，輸入中頻60 MHz的單頻正弦信號，測出了隨著輸入信號功率變化，輸出信號強度變化的趨勢，測試結(jié)果如表1所示。

表1動態(tài)測試結(jié)果

按照表1中數(shù)據(jù)做出系統(tǒng)動態(tài)范圍曲線，如圖6所示。

圖6　動態(tài)測試曲線

可以看出，在輸入動態(tài)范圍-96 dB～16 dB內(nèi)，輸出都能保持很好的線性，即系統(tǒng)動態(tài)達到112 dB。

5　結(jié)束語

本系統(tǒng)采用雙通道ADC采樣方式解決了天氣雷達數(shù)字中頻接收機動態(tài)范圍受限的重要問題，并且這種解決方式是在不破壞原始信號相位信息的基礎(chǔ)上，適合應(yīng)用于多普勒天氣雷達當(dāng)中。

本文采用的系統(tǒng)硬件主要包括 Altera Striax II FPGA開發(fā)平臺、功分器、衰減器。硬件代碼采用QuartusII6.0集成開發(fā)環(huán)境進行仿真綜合。經(jīng)過綜合后，F(xiàn)PGA程序占用資源為：LC寄存器10 519個、RAM745 688位、鎖相環(huán)1個、DSP單元280個。在雷達中頻60 MHz，采樣頻率48 MHz，帶寬1.2 MHz的工作條件下，該硬件系統(tǒng)完成了對雙通道雷達I/Q信號的解調(diào)，并且工作正常。實驗表明：該設(shè)計方案是可行的，能夠很好地提高雷達接收機的動態(tài)范圍，并應(yīng)用于多普勒天氣雷達數(shù)字中頻接收機當(dāng)中。

[1]弋穩(wěn). 雷達接收機技術(shù)[M]. 北京：電子工業(yè)出版社, 2005.

YI Wen. Radar receiver technology[M]. Beijing: Publishing House of Electronics Industry, 2005.

[2]楊小牛. 軟件無線電原理與應(yīng)用[M]. 北京：電子工業(yè)出版社，2001.

YANG Xiaoniu. Principle and application of software radio[M]. Beijing: Publishing House of Electronics Industry, 2001.

[3]王江濤. 一種數(shù)字非相干接收機的實現(xiàn)[J]. 微波學(xué)報, 2014, 31(增刊): 598-601.

WANG Jiangtao. Realization of a digital noncoherent receiver[J]. Journal of Microwaves, 2014, 31(Supplement): 598-601.

[4]丁家會. 數(shù)字化接收機極限性能指標(biāo)的研究[J]. 現(xiàn)代雷達，2005，27(9)：72-74.

DING Jiahui. A study on limit performance specification of a digital receiver[J]. Modern Radar, 2005, 27(9): 72-74.

[5]肖學(xué)穎, 王世軍. 多通道雷達中頻數(shù)字接收機的工程實現(xiàn)方法研究[J]. 雷達與對抗, 2007(1)：44-47.

XIAO Xueying, WANG Shijun. Implementation of the multi-channel IF digital receiver for radars[J]. Radar and Countermeasure, 2007(1): 44-47.

馬勇男，1983年生，工程師。研究方向為航空氣象裝備管理。

Implementation of Dual-channel Digital IF Processing System for Terminal Weather Radar

MA Yong

(Meteorology Division of South West Air Traffic Management Bureau,Chengdu 610225, China)

Terminal weather radar requires the ability to identify different weather conditions from complex weather conditions for safety of aviation, so the design of large dynamic receiver and noise performance is critical. Based on the constraint analysis of analog digital converter (ADC) to the dynamic range of radar IF receiver, according to band-pass sampling and digital down conversion principle, a digital IF processing system based on field programmable gate array with dual-channel ADC sampling is achieved, and system design principles, methods and test results are given. Through the realization of digital IF processing system with dual-channel ADC sampling, the dynamic range of weather radar receiver is well improved, and this system can be applied to the digital IF receiver of terminal Doppler weather radar.

terminal weather radar; digital IF; digital down conversion; digital filter; dynamic range

馬勇Email：mayong_2002@qq.com

2016-03-18

2016-05-20

TN959.4

A

1004-7859(2016)07-0067-05

·收/發(fā)技術(shù)·

DOI：10.16592/ j.cnki.1004-7859.2016.07.017