基于FPGA的等效時間采樣

劉建博，郭文秀，張 捷，伍守豪

（1.西北工業大學 電子信息學院，陜西 西安 710072；2.深圳清華大學研究院數字電視系統重點實驗室，廣東 深圳 518057）

在現代電子測量、通訊系統以及生物醫學等領域，經常涉及對寬帶模擬信號進行數據采集和存儲，以便計算機進一步進行數據處理。為了對高速模擬信號進行不失真采集，根據奈奎斯特定理，采樣頻率必須為信號頻率的2倍以上，但在電阻抗多頻及參數成像技術中正交序列數字解調法的抗噪性能對信號每周期的采樣點數決定，采樣點數越多，抗噪性能越高[1]。當采樣信號頻率很高時，為了在被采樣信號的一周期內多采樣，就需要提高采樣時鐘的頻率，但是由于系統的ADC器件時鐘速率并不能達到要求的高頻速率或者存儲處理速度等不能滿足要求因此我們可以采用低速ADC器件通過等效時間采樣來對寬帶模擬信號進行數據采集從而使系統易于實現。

1 等效時間采樣原理

等效時間采樣技術是把周期性或準周期性的高頻、快速信號變換為低頻的慢速信號。在電路上只對取樣前的電路具有高頻的要求，大大降低采樣變換后的信號處理、顯示電路對速度的要求，簡化了整個系統的設計難度[3]。等效時間采樣分為順序采樣 （sequential equivalent sampling）、 隨機采樣（random equivalent sampling）以及結合這兩種方式的混合等效采樣（compound equivalent sampling）[2]。在文獻[3]、[4]中分別介紹了兩種硬件實現的等效時間采樣中的順序采樣。

下面我將介紹等效時間采樣中的混合時間采樣，對于周期性信號的等效時間采樣如圖1（a）所示。

在第一周期中的橫軸（時間）的第2與第6處的時鐘上升沿對模擬信號進行采樣，圖中的箭頭表示采樣時刻。在一個周期中可以采集兩個點，緊接著在第二個周期橫軸的第11與第15處的時鐘上升沿對模擬信號進行采樣。為了方便觀察在此將第一至第五周期的波形縱向排列。可以看到第二周期比第一周的采樣點距離各自周期起始點的時間晚了一個時鐘周期。第三周期比第二周的采樣點距離第三周期起始點的時間晚了一個時鐘周期。在第四周期進行采樣時我們可以發現第二個采樣點已經進入第五周期。如果我們在第五周期周試圖繼續用以上方式進行采樣即第五周期比第四周的采樣點距離起始點的時間晚一個時鐘周期，那么我們會發現在第五周期的采樣起始點采樣到的值重復了第一周期采樣到的數值。所以此時我們可以終止采樣那么我們就得到了如圖1中的第6個波形示意圖所表示的在一個周期的正弦波形中采到的8個數據點。

在文獻[5]中給出了等效時間采樣中每個周期可以采集多個點時的理論依據，在文獻[6]中給出了等效時間采樣中每個周期可以采集單個點時的理論依據。

圖1 等效時間采樣示意圖Fig.1 Schematic diagram of ETS

我們通過將高頻時鐘進行分頻已達到或者接近滿足處理速度時鐘要求。在圖1（b）中幅度最小的時鐘信號為采樣時鐘。由圖1（b）可以很清楚的看到分頻后的時鐘波形，分頻后的時鐘波形在時鐘的上升沿對信號進行采樣，那么就會得到如圖1（a）中所表示的等效時間采樣。

2 基于FPGA的等效時間采樣實現

2.1 系統硬件實現框圖

系統的總體框圖如圖2，FPGA控制的等效采樣時鐘連接到ADC器件的時鐘部分，ADC器件在時鐘的控制下對寬帶模擬信號進行采樣，采集到的數據傳送到FPGA中的FIFO，FPGA再將FPGA中 FIFO的數據傳遞到USB中的FIFO，然后USB將USB中FIFO數據推送到計算機，計算機對接收到的數據進行重構處理。對于信號周期的獲取，在電阻抗多頻及參數成像技術中采集信號的周期是由發送信號的周期決定，而對于其他復雜周期信號的周期獲得可以通過[5]所采用的方法獲得。

圖2 系統方案框圖Fig.2 System block diagram

2.2 等效時間采樣時鐘的程序實現

圖3 等效時間采樣模塊圖Fig.3 ETSmodule diagram

圖3 展示了基于FPGA生成的等效時間采樣模塊的輸入端口與輸出端口。其中CLK表示高頻時鐘的輸入，RESET表示的是復位輸入端，FREN_CON表示的是分頻控制輸入用于控制高頻時鐘的分頻數，SANM_CONT表示的是模擬信號的周期包含多少個高頻時鐘信號的波形，CLK_ADC_OUT表示的是輸出時鐘端口，此端口連接到模數轉換器件（ADC）的時鐘輸入端口。

以下是實現等效時間采樣所需時鐘的代碼：

SIGNAL SAMP_CONTS:STD_LOGIC_VECTOR （11 DOWNTO 0）:=（OTHERS=>'0'）；

SIGNAL ADC_CLK_BANK:STD_LOGIC_VECTOR （11 DOWNTO 0）:=（OTHERS=>'0'）；

SIGNAL CLK_CNT:INTEGER RANGE 0 TO 5000:=0；

SIGNAL CLK_TANK:STD_LOGIC:='0'；

SIGNAL EN :STD_LOGIC:='1'；

BEGIN

PROCESS（CLK，RESET）

BEGIN

IF RESET='1'THEN EN<='1'；

ELSE

IF CLK'EVENT AND CLK='1'THEN

IF SAMP_CONTS

SAMP_CONTS<=SAMP_CONTS+'1'；

IF ADC_CLK_BANK

ADC_CLK_BANK<=ADC_CLK_BANK+'1'；

ELSE EN<='0'；

END IF；

ELSE

ADC_CLK_BANK<=（OTHERS=>'0'）；

SAMP_CONTS<=（OTHERS=>'0'）；

EN<='1'；

END IF；

END IF；

END IF；

END PROCESS；

---分頻控制部分程序，作ADC的時鐘使用

PROCESS（CLK，RESET，EN）

BEGIN

IF RESET='1'THEN

CLK_TANK<='0'；

CLK_CNT<=0；

ELSE

IF EN='1'THEN

IF CLK'EVENT AND CLK='1'THEN

IF （CLK_CNT=（CONV_INTEGER （FREN_CON）/2）-1）THEN

CLK_TANK<=NOT CLK_TANK；

CLK_CNT<=0；

ELSE CLK_CNT<=CLK_CNT+1；

END IF；

END IF；

ELSE

CLK_TANK<='0'；

CLK_CNT<=0；

END IF；

END IF；

END PROCESS；

CLK_ADC_OUT<=CLK_TANK；

END；

3 波形仿真

圖4中的波形仿真是以模擬信號的一周期等于8個CLK時鐘周期，CLK_ADC_OUT是對CLK進行4分頻且分頻后的時鐘占空比為50%為假設的。1號箭頭指向的時鐘上升沿標志著第一周期結束，上升沿之后進入第二周期。同理，2號箭頭所指時鐘的上升沿標志著第二周期的結束，上升沿之后標志著進入第三周期。

圖4 波形仿真Fig.4 Waveform simulation

在第一個周期中從CLK的第一個上升沿開始計時同時對CLK進行分頻可以得到CLK_ADC_OUT時鐘信號，在第一周期中在CLK的第二個上升沿CLK_ADC_OUT電平翻轉（存在延時），在第二周期中在第三個上升沿CLK_ADC_OUT電平翻轉，在第三個周期中在CLK的第四個上升沿CLK_ADC_OUT電平翻轉。可以看出波形仿真圖是對圖1（a）、（b）兩圖表達時鐘的實現。在這里應該注意到，在第一周期中雖然也有8個CLK的上升沿，但是并沒有表示出如1號箭頭所指CLK時鐘上升沿之后與第二周期第一個CLK時鐘上升沿之間的波形。

4 結 論

本文介紹了等效時間采樣的基本原理、系統實現的具體方案。等效時間采樣技術實現了利用低速的ADC器件對寬帶模擬信號的采集，降低了系統對ADC器件的要求以及系統實現的復雜度。本文介紹的等效時間采樣技術由于使用了FPGA采樣技術，使得在被采樣信號的一個周期中相較于一個周期僅能采集一個點的順序等效時間采樣有很大的提高，并且可以控制被采集信號一個周期中的采集點數從而可以根據后續器件處理速度實現變頻控制采樣。通過FPGA實現等效采樣時間，降低了系統實現的復雜度，同時可以十分方便的對代碼進行修改使系統的調試更加簡便。

[1]史學濤.用于電阻抗多頻及參數成像的多頻數據采集技術研究[D].西安:第四軍醫大學，2001.

[2]葉旭.導波式雷達物位計信號處理方法研究與實驗[D].合肥:合肥工業大學，2012.

[3]張杰，曹思樟，祝懷標，等.基于游標法的等效時間采樣模塊的設計[J].計量與測試技術2012，39（4）:59-60，63.ZHANG Jie，CAO Si-zhang，ZHU Huai-biao，et al.Design of equivalent time sanding module based on vernier method[J].Measurement and Test Technology，2012，39（4）:59-60，63.

[4]陳小橋，劉愛榮，亓長軍.基于等效時間采樣的高速數據采集技術[J].電測與儀表，2002，39（440）:23-26.CHEN Xiao-qiao，LIU Ai-rong，QI Chang-jun.High speed data acquisition system based on equivalent time sampling[J].Electrical Measurement&Instrumentation，2002，39 （440）:23-26.

[5]任斌，余成，陳衛，等.基于EDA技術的等效采樣的設計實現[J].微計算機信息，2007，23（5-2）:293-295.REN Bin，YU Cheng，CHEN Wei，et al.Design and realization of equivalent-time sampling based on EDA[J].Micro Computer Information，2007，23（5-2）:293-295.

[6]周維，王赤，田茂，等.基于等效時間采樣的探地雷達回波信號采樣方法研究[J].雷達科學與技術，2004（2）:43-47.ZHOU Wei，WANG Chi，TIAN Mao，et al.Research on a sample-collecting method of echo data of ground penetrating radar based on equivalent time sampling[J].Radar Science and Technology，2004（2）:43-47.