董 亮,汪 敏,高亦菲,高冠男

摘 要:首先分析Chirp函數(shù)在頻域上的一般特性,并且分析Altrea 公司提供的數(shù)控振蕩器知識產(chǎn)權(quán)核(NCO IP core)的輸入/輸出特性,通過MegaCore環(huán)境確定其輸入控制字,通過外圍邏輯電路實時向NCO IP core調(diào)入控制頻率控制字以達到改變輸出頻率的目的,并通過在示波器上觀測FPGA的運行情況,驗證了該設(shè)計具有很好的輸出效果。

關(guān)鍵詞:NCO IP core;FPGA;Chirp函數(shù);MegaCore

中圖分類號:TP311文獻標識碼:A

文章編號:1004-373X(2009)20-020-03

Design and Implementation of Chirp Function Based on NCO IP core

DONG Liang1,2,WANG Min2,GAO Yifei2,GAO Guannan2

(1.Graduate School of Chinese Academy of Science,Beijing,100049,China;2.Yunnan Observatory,Chinese Academy of Science,Kunming,650011,China)

Abstract:This design analyses the character of Chirp function in the time and frequency domain zone.And after analyzing the input/output feature of the NCO IP core.The different output frequency by writing different controlling words which can be determined from the MegaCore environment.The FPGA′s running situation by oscillograph to inspect and verify feature of the design is detected.The conclusion that this design can get the Chirp′s function perfectly.

Keywords:NCO IP core;FPGA;Chirp function;MegaCore

0 引 言

IP就是知識產(chǎn)權(quán)核或者知識產(chǎn)權(quán)模塊的意思。在EDA技術(shù)和開發(fā)領(lǐng)域具有十分重要的作用,在半導(dǎo)體產(chǎn)業(yè)中IP定義為用于ASIC或FPGA/CPLD中預(yù)先設(shè)計好的電路功能模塊。IP可以分為軟IP,固IP和硬IP三種。

隨著電子系統(tǒng)的越來越復(fù)雜,PLD設(shè)計的越來越龐大,這就增加了市場對IP核的需求,各大FPGA/CPLD廠商陸續(xù)推出了許多IP核。例如:FIR(有限沖擊響應(yīng))數(shù)字濾波器 core;FFT(快速傅里葉變換)core,NCO(數(shù)控振蕩器)core等,在設(shè)計中如果使用了這些知識產(chǎn)權(quán)核可以大大簡化FPGA/CPLD的設(shè)計,加速設(shè)計速度,縮短研發(fā)周期,并且較之于開發(fā)者自己的設(shè)計程序,這些IP有更好的運算精度、速度、SFDR參數(shù)、SNR參數(shù)等,達到良好的效果!

由于電磁波在傳輸過程中,經(jīng)過色散介質(zhì),如不均勻的波導(dǎo)經(jīng)過高空電離層時會發(fā)生色散現(xiàn)象。Chirp函數(shù)在射電天文信號的消色散處理中發(fā)揮著重要的作用,研究在FPGA中實現(xiàn)Chirp函數(shù)是基于FPGA的射電宇宙信號處理的重要組成部分。如圖1所示。

圖1 基于FPGA的射電宇宙信號處理框圖

該設(shè)計就是通過實時的改變NCO IP core的輸入頻率控制字的辦法,數(shù)控頻率輸出的辦法實現(xiàn)Chirp函數(shù)。

1 系統(tǒng)總體設(shè)計

Chirp函數(shù)根據(jù)輸出頻率的遞變規(guī)律一般分為兩種[1,2]:線性Chirp函數(shù)和非線性Chirp函數(shù),以下是兩種Chirp函數(shù)在頻域上的表現(xiàn)如圖2,圖3所示。

圖2 線性Chirp函數(shù)

圖3 非線性Chirp函數(shù)

從圖2,圖3可以看出Chirp函數(shù)的頻率輸出與時間的f-t關(guān)系可以總結(jié)為:

對于線性Chirp函數(shù):

在連續(xù)域時間域內(nèi)有關(guān)系式:

fout=kt+f0(1)

式中:k為常數(shù);f0為初始輸出頻率;t為連續(xù)時間。

在離散時間域有關(guān)系式:

fout=kn+f0(2)

式中:k為常數(shù);f0為初始輸出頻率;n為采樣點。

對于非線性Chirp函數(shù):在連續(xù)域時間域內(nèi)有關(guān)系式:

fout=f(t)+f0(3)

式中:f(t)為非線性函數(shù);f0為初始輸出頻率;t為連續(xù)時間。

在離散時間域有關(guān)系式:

fout=f(n)+f0(4)

式中:f(n)為非線性函數(shù);f0為初始輸出頻率;n為采樣點。

由上式可以看出Chirp函數(shù)在每一個時刻點具有不同的頻率輸出,而根據(jù)具體的頻率變化的需要在每一個時刻點實時的改變其頻率控制字是實現(xiàn)Chirp函數(shù)的關(guān)鍵。其算法框圖如圖4所示。

圖4 總體設(shè)計框圖

2 NCO IP Core

數(shù)字壓控振蕩器知識產(chǎn)權(quán)核[8](Numerically Controlled Oscillators Intellectual Property Core,NCO IP Core),通過多種算法(相位累加或者CORDIC算法,在此不一一贅述),實現(xiàn)了一個離散幅度和時間的正弦波信號輸出。輸入控制字和輸出頻率之間滿足以下方程:

s(nT)=Asin[2π(f0+fFM)nT+φPM+φDITH)](5)

式中:T為該模塊的工作時鐘;f0是由輸入頻率控制字φINC決定的初始頻率;fFM是由調(diào)制頻率控制字φFM決定的調(diào)制頻率;φPM為該輸出正弦波的調(diào)制相位,φPM=P/2Pwidth,由輸入控制字P的比特位數(shù)(Pwidth)決定了它的精度;φDITH為模塊內(nèi)部自身的不穩(wěn)定而引起的相位雜散(噪聲);幅值量A=2N-1,其中N為幅值精度取值在4~32之間。

該設(shè)計中僅采用通過改變頻率控制字φINC,以實現(xiàn)改變頻率輸出的目的,為此式(5)可以簡化為:

s(nT)=Asin[2π(f0nT+φDITH)](6)

式中,f0由給定的頻率控制字φINC決定,滿足如下方程:

fo=φincfclk/2M(7)

式中:M為累加器精度;fclk為該模塊的輸入時鐘頻率,單位為Hz。例如:在fclk=100 MHz的情況下,如果需要f0=10 MHz的輸出,φINC的計算如下:

φINC=10×106100×106×232=4 294 967 296(8)

通過Altera公司的FPGA編程軟件Quartus Ⅱ提供的MegaWizard Plug-In Manager功能,在NCO IP Core 參數(shù)配置中自動對φINC的計算,很容易得出在輸入頻率的條件下所需輸出頻率的φINC(累加器精度為32 b的情況下)如圖6所示。

圖5 NCO IP內(nèi)部算法框圖

圖6 利用MegaWizard Plug-In Manager計算φINC

值得注意的是:在MegaWizard Plug-In Manage中,φINC的精度只保留到了百位。

3 頻率控制字寄存器及驅(qū)動單元的設(shè)計

頻率控制字寄存器為一個保存有N個輸出頻率所需的相位累加控制字的片上ROM單元,其作用在驅(qū)動單元輸入地址控制字的作用下實時向NCO IP Core調(diào)入所需要的φINC,在該設(shè)計中Chirp函數(shù)的頻率變化規(guī)律是從1 MHz步進1 MHz輸出到16 MHz。在該設(shè)計中選擇的累加器精度為32 b,為此選擇的邏輯單元的規(guī)律為如表1所示。

表1 選擇的邏輯單元的規(guī)律

第N頻點φINCf0/MHz第N頻點φINCf0/MHz

142 949 67319386 547 0579

285 899 346210429 496 73010

3128 849 019311472 446 40311

4171 798 692412515 396 07612

5214 748 365513558 345 74813

6257 698 038614601 295 42114

7300 647 711715644 245 09415

8343 597 384816687 194 76716

為此,建立一個深度為16,每個存儲單元字長32 b位的ROM,將表1內(nèi)所有φINC數(shù)據(jù)保存至nco116.mif文件中,在ROM建立時調(diào)用該mif文件。如圖7所示。

圖7 建立的mif文件內(nèi)部數(shù)據(jù)一覽

兩者連接關(guān)系如圖8所示。

圖8 頻率控制字寄存器與NCO IP Core的連接圖

在設(shè)計中,通過不同時間點向頻率控制字寄存器寫入不同的地址信號驅(qū)動,使存儲器輸出不同的頻率控制字驅(qū)動NCO IP Core,產(chǎn)生不同的頻率信號輸出。

該設(shè)計中采用兩個計數(shù)器級聯(lián)作為驅(qū)動單元,首先第一級計數(shù)器將鐘頻率降至需要的Chirp函數(shù)輸出某頻點的穩(wěn)定時間范圍,將第一級計數(shù)器的進位端作為第二級計數(shù)器的時鐘輸入端;第二級計數(shù)器的作用是,產(chǎn)生地址信號以驅(qū)動頻率控制字存儲器輸出相應(yīng)的控制字,當(dāng)前級進位信號有效時該計數(shù)器輸出加“1”,以達到改變頻率輸出的目的,其連接電路圖如圖9所示。

圖9 驅(qū)動邏輯電路的實現(xiàn)圖

4 仿真與驗證

將該設(shè)計通過將程序下載到Altera公司生產(chǎn)的DSP開發(fā)板[2,3](型號DK-DSP-2C70N)中進行仿真,其核心FPGA(型號為EP2C70F672C6)的資源使用情況如圖10所示。

圖10 該設(shè)計中核心FPGA的資源使用情況

并通過該開發(fā)板上D/A轉(zhuǎn)換器輸出模擬波形(只截取了4個時刻的圖樣)如圖11所示。

圖11 該設(shè)計在示波器(型號DSO3201A)上的顯示波形

通過圖11可以看出該設(shè)計能很好地完成掃頻輸出的功能,并且雜波分量很小,干擾很小。

5 結(jié) 語

該設(shè)計通過采用技術(shù)成熟的NCO IP Core完成,其優(yōu)勢在于:

(1) 利用了成熟的FPGA知識產(chǎn)權(quán)技術(shù),使得設(shè)計更加簡便并易于移植;

(2) 利用NCO IP Core的高穩(wěn)定性,使得Chirp函數(shù)的各項噪聲較之于其他設(shè)計更小,有利于對射電天文這樣微弱信號的處理,減少了處理帶來的各種噪聲。

參考文獻

[1]Altera Corporation.NCO MegaCore Function User Guide[Z].2006.

[2]Altera Corporation.Cyclone II DSP Development Board Re-ference Manual[Z].2006.

[3]Altera Corporation.DSP Development Kit,Cyclone II Edition Getting Started User Guide[Z].2006.

[4] 卞小林,胡生亮,趙佩麗,等.基于Nios Ⅱ軟核處理器的波形發(fā)生器設(shè)計[J].電子測量技術(shù),2008,31(5):181-183.

[5]羅苑棠.CPLD/FPGA常用模塊與綜合系統(tǒng)設(shè)計實例精講[M].北京:電子工業(yè)出版社,2007.

[6]Lyne A G,Joshi B C,Ikin T,et al.Holloway CORBA-A Digital Receiver at Jodrell Bank.

[7]潘松,黃繼業(yè),曾毓.SoPC技術(shù)實用教程[M].北京:清華大學(xué)出版社,2005.

[8]Altera Corporation.Cyclone II Device Handbook(All sections)[Z].2005.

[9]孟憲元,錢偉康.FPGA嵌入式系統(tǒng)設(shè)計[M].北京:電子工業(yè)出版社,2007.