郭業(yè)才， 萬逸儒， 吳 星， 胡 崢

(1.南京信息工程大學 電子與信息工程學院，南京 210044； 2.江蘇省大氣環(huán)境與裝備技術(shù)協(xié)同創(chuàng)新中心，南京 210044)

0 引 言

針對現(xiàn)代衛(wèi)星通信盲均衡系統(tǒng)產(chǎn)業(yè)化實現(xiàn)的需求，搭建一個完整的電子測量系統(tǒng)顯得十分必要。然而，只有在能夠采集到某類信號時，這些工具才能進行測量，絕大多數(shù)情況下，這些信號是沒有的，所以研制可以穩(wěn)定精確地產(chǎn)生調(diào)制信號的發(fā)生器能彌補電子測量系統(tǒng)激勵不足的問題，同時為了方便工程師探測盲均衡系統(tǒng)性能的極限，信號發(fā)生器必須還具有高階調(diào)制及各個參數(shù)靈活可變的功能，使余量測試和檢定成為可能。信號發(fā)生器的性能指標和考慮因素主要有：存儲深度、采樣速率、帶寬、幅度分辨率、定時分辨率、區(qū)域位移功能、輸出通道數(shù)量、數(shù)字輸出、排序、集成編輯器、數(shù)據(jù)導入功能、無雜散動態(tài)范圍、相位噪聲密度等。目前，國內(nèi)的數(shù)字合成信號發(fā)生器在技術(shù)指標上與國外相比還存在很大的差距[1]，例如，比較有名的生產(chǎn)信號發(fā)生器的大公司有安捷倫、美國國家儀器NI、羅德與施瓦茨、泰克科技、臺灣百科精密儀器、利達電子株式會社等，其中安捷倫儀器的PSG系列信號發(fā)生器已經(jīng)可達毫米波范圍，最高可以產(chǎn)生70GHz的射頻信號，高品質(zhì)的射頻調(diào)制帶寬支持為雷達、電子戰(zhàn)、和衛(wèi)星通信測試生成更真實的寬帶信號，最大限度提高多普勒雷達或接收機阻塞測試的接收機靈敏度，或在相位噪聲測試系統(tǒng)中用作理想的基準：100 Hz 頻偏時為-91 dBc/Hz；10 kHz 頻偏(67 GHz)時為-126 dBc/Hz，低雜散和諧波(分別為 10 GHz 處-70 dBc 和 55 dBc)的出色性能可以改善信號質(zhì)量并提高測試精度。

1 系統(tǒng)設計

本設計方案使用硬件邏輯可編程的FPGA設計基帶硬件結(jié)構(gòu)及基于軟件無線電思想的通用、標準化、模塊的硬件設計平臺，通過Altera FPGA內(nèi)部集成的NiosⅡ軟核[2]設置參數(shù)值。系統(tǒng)總體框架，如圖1所示。圖中，AD9957上變頻芯片工作在1GHz頻率，內(nèi)部可通過8-252可編程插值輸出同步時鐘到FPGA內(nèi)部可變參數(shù)基帶處理模塊，AD9957通過T1-1T射頻變壓器將差分型號轉(zhuǎn)換為單端信號，并經(jīng)過160 MHz截止頻率的7階橢圓模擬重構(gòu)濾波器輸出到BNC接口端。同時，F(xiàn)PGA內(nèi)部用QSYS搭建SOPC系統(tǒng)，該系統(tǒng)集成了NiosⅡ/fast軟核、高速緩存單元、內(nèi)部緊密耦合寄存器、SDRAM控制模塊、EPCS控制模塊、JTAG_UART仿真器控制模塊、系統(tǒng)時鐘單元、時間戳單元、總線(8080總線、SPI總線、I2C總線)控制模塊。同時該設計可以通過上位機和基于圖形庫支持脫機操作的觸控屏來完成人機友好界面的交互。

圖1 系統(tǒng)總體框架圖

2 數(shù)字脈沖成型原理

2.1 升余弦濾波器

由于矩形脈沖的頻譜形狀是sinc響應，在碼元速率為f0的整數(shù)倍處周期性的出現(xiàn)零點，造成基帶信號在頻譜上無限制的拓展，在大部分數(shù)據(jù)傳輸應用中，所傳輸?shù)男盘柋仨毾拗圃谝欢ǖ膸拑?nèi)[3-5].文獻[6]中提供了在衛(wèi)星對地數(shù)據(jù)傳輸系統(tǒng)領(lǐng)域，國際空間數(shù)據(jù)系統(tǒng)咨詢委員會(Consultative Committee for Space Data Systems，CCSDS)建議中給出的固態(tài)功放(Solid-state Power Amplifier，SSPA)下各種調(diào)制方式的頻譜密度示意圖。由圖可見，限制矩形脈沖帶寬的低通濾波器是必須的，這種將矩形脈沖形狀從門形狀變?yōu)闊o陡峭邊沿的平滑輪廓線的低通濾波器稱之為脈沖整形濾波器。然而，理想低通濾波器在頻域表現(xiàn)為門形狀，因此在時域必定表現(xiàn)為阻尼衰減的特性，碼元經(jīng)過阻尼衰減的拖尾會產(chǎn)生對鄰近碼元有干擾的紋波，紋波會導致數(shù)據(jù)解碼的錯誤。綜上所述，必須選擇既能限制頻域脈沖帶寬又能保證時域波形無碼間干擾的濾波器。滿足這兩個條件的著名濾波器就是升余弦濾波器，其頻率響應為：

(1)

式中：T是碼元周期；α是余弦滾降系數(shù)；ω是角頻率。

相應的時域響應形式為：

(2)

2.2 余弦滾降系數(shù)

滾降系數(shù)的概念是建立在升余弦特性頻譜帶寬相對于理想低通特性帶寬的基礎(chǔ)上提出來的，余弦滾降系數(shù)決定了兩個比較重要的頻率點，一個是阻帶頻率點fstop，另一個是通帶頻率點fpass，它們相對于滾降系數(shù)的關(guān)系式為：

(3)

式中：f0是脈沖速率；f0/2是奈奎斯特頻率點，也稱理想低通濾波器的截止頻率。

余弦滾降系數(shù)與濾波器帶寬成正相關(guān)關(guān)系，在插值倍數(shù)為8、雙邊截斷符號碼元數(shù)為8、碼元速率為10 MB、采樣速率為80 MB條件下，圖2和圖3分別仿真了滾降系數(shù)為1,0.5,0.1時的濾波器時域和頻域特性。

圖2 滾降系數(shù)為1,0.5,0.1時的時域特性圖

圖3 滾降系數(shù)為1,0.5,0.1時的頻域特性圖

由圖可見，在滾降系數(shù)為0.1的情況下，頻域特性急劇惡化，圖中產(chǎn)生了最大約0.1左右的帶內(nèi)波紋，這里可以通過增加截取符號碼元數(shù)來彌補波紋的惡化，但會增加濾波器階數(shù)從而提高了硬件系統(tǒng)復雜度和資源耗用。從時域角度來說，滾降系數(shù)越小時域拖尾也會越大，接收端符號碼同步不可能完全沒有誤差，只要有一點誤差，就會因為較大的時域拖尾而造成很大的碼間干擾，惡化系統(tǒng)的信噪比，影響系統(tǒng)的誤比特率。當然滾降系數(shù)也不能過大，過大會造成系統(tǒng)頻帶利用率的下降，同時會引入更寬帶寬的帶內(nèi)噪聲，對數(shù)模轉(zhuǎn)換后的低通濾波器和接收端模數(shù)轉(zhuǎn)換前的抗混疊濾波來說，在過寬的帶寬內(nèi)保證幅度一致和線性相位是非常困難的，文獻[3]中指出了系統(tǒng)誤碼率在滾降系數(shù)取0.5～1之間時區(qū)別不大，因此設計了可調(diào)滾降系數(shù)范圍為0.3～0.5。

2.3 雙邊截斷碼元數(shù)

雙邊截斷碼元數(shù)與濾波器階數(shù)復雜度成正比，因此要使用滿足濾波要求的最小雙邊截斷碼元值，文獻[7]中給出了基于浮點運算模擬生成的通帶邊緣頻率處誤差、奈奎斯特頻率處誤差、阻帶頻率處誤差隨余弦滾降系數(shù)和雙邊截斷碼元數(shù)變化的三維曲面圖，說明了可接受的最小截斷碼元值是由頻域特性的帶外第一個衰減旁瓣峰值決定的，尤其是在滾降系數(shù)極小時，更需要考慮最低阻帶衰減特性。在插值倍數(shù)為8、雙邊截斷碼元數(shù)為8、碼元速率為10MB、采樣速率為80MB條件下，圖4分別仿真了余弦滾降系數(shù)為0.5，0.3，0.2時的頻域特性曲線圖。

圖4 雙邊截斷碼元數(shù)為8情況下不同滾降系數(shù)頻域特性曲線圖

可以看出，滾降系數(shù)為0.5時阻帶第一個衰減旁瓣峰值在-40 dB以下，滾降系數(shù)在0.3時第一個衰減旁瓣峰值正好在-40 dB上，滾降系數(shù)在0.2時第一個衰減旁瓣峰值在-30 dB左右，此時會引入更多的旁瓣噪聲，導致信噪比的惡化，因此余弦滾降系數(shù)下限設計在0.3 ，雙邊截斷碼元數(shù)為固定值8。

截斷碼元數(shù)之所以取固定值是由兩方面因素決定的，① 截斷碼元數(shù)在多相時變系數(shù)插值濾波器中直接映射為每組子濾波器的系數(shù)個數(shù)，也就是FPGA中乘法器使用單元數(shù)，過多的截斷碼元會導致乘法器資源的消耗。② 濾波器抽頭階數(shù)一般是由Matlab產(chǎn)生的浮點數(shù)，需要進行一定比例的放大并取整，放大倍數(shù)為

(4)

2.4 插值倍數(shù)

由奈奎斯特采樣定理可知，采樣頻率必須大于或等于帶內(nèi)信號的最高頻率，插值倍數(shù)是采樣頻率相對于碼元傳輸頻率的比值，根據(jù)數(shù)字采樣定理可知，雙邊截斷碼元是對升余弦濾波器時域響應的有限長度截取，而在頻域則表現(xiàn)為按采樣頻率進行周期的幅頻特性的搬移。原始的周期脈沖頻域形狀是眾所周知的sinc響應，零點總是出現(xiàn)在碼元速率f0的整數(shù)倍處，多進制波形可看做是多個二進制振幅波形的疊加，因此頻域特性也是多個信號的頻域特性之和，盡管疊加后頻譜特性較復雜，但就零點出現(xiàn)特性而言，在f0的整數(shù)倍上也總是表現(xiàn)為0幅值，當插值倍數(shù)為整數(shù)時，濾波器按采樣頻率搬移的通帶正好出現(xiàn)在f0的整數(shù)倍上，可以有效地降低頻率搬移對濾波特性的影響。在雙邊截斷碼元數(shù)為8、余弦滾降系數(shù)為0.5、碼元速率為10 MB條件下，圖5給出了插值倍數(shù)分別為4,8,16即采樣速率為40、80、160 MB的幅頻特性圖。

圖5 不同插值倍數(shù)下的幅頻特性圖

從圖中可以看出，不同插值倍數(shù)雖然不會對第一個衰減旁瓣峰值有所影響，但會影響后級旁瓣峰值大小，采樣率越高，時域成型效果越好，頻域后級旁瓣衰減越大，但會增加濾波器的階數(shù)，提高FIR復雜度，也會使輸出碼元速率增大，增大了系統(tǒng)實現(xiàn)難度。

3 可變參數(shù)基帶處理模塊設計

3.1 串行IQ鏈路多相時變系數(shù)濾波器的結(jié)構(gòu)

AD9957、AD9857、AD9854等具有正交調(diào)制功能的上變頻芯片，它們在并行端口上交替出現(xiàn)I和Q路數(shù)據(jù)字，并在時鐘有效沿進行捕捉，因此在基帶處理部分，也可以設計成串行IQ兩路串行的串并轉(zhuǎn)換、星座映射、成型濾波，這樣可以不用進行時鐘分層以及節(jié)省下一路的濾波資源和映射編碼邏輯，傳統(tǒng)的FPGA插值濾波器采用的是全并行分布式算法原理[5],這種基于查找表結(jié)構(gòu)的濾波器可以大大提高電路的執(zhí)行速度及不需要乘法資源，目前大多數(shù)FPGA/CPLD都集成了DSP處理單元，內(nèi)部嵌入了細顆粒單周期乘法器，用這些乘法器來設計卷積濾波器將大大降低系統(tǒng)復雜度。

假設插值濾波器抽頭為

H=[h(1),h(2),h(3),…,h(N)]

式中：N=D×M；D是雙邊截斷碼元數(shù)；M為插值倍數(shù)。將H分解為M個子濾波器：

(5)

這種結(jié)構(gòu)的插值濾波器只需要D個乘法器資源，相對于采用FIR的IP核要節(jié)省M倍的乘法器，圖6給出了串行IQ鏈路多相時變系數(shù)濾波器[4]的結(jié)構(gòu)圖。

在每個PDCLK時鐘周期，計數(shù)控制模塊都會控制子濾波器系數(shù)組都會進行更新，當M組子濾波器都遍歷更新過一遍之后的第一個PDCLK周期，計數(shù)控制模塊控制I路用戶數(shù)據(jù)更新入I路移位寄存器，第二個PDCLK周期，Q路用戶數(shù)據(jù)更新入Q路移位寄存器，并在每一個PDCLK周期由計數(shù)控制的MUX選擇器控制IQ兩路鏈路串行進入并行乘加單元進行濾波[8-15]，并行乘加單元輸出18位(1位符號位+17位數(shù)據(jù)位)數(shù)據(jù)給AD9957，因為最高使用256QAM，也就是16ASK正交調(diào)制，共16個電平位，電平位間隔為2，需要1位符號位+4位數(shù)據(jù)位，所以用戶數(shù)據(jù)設置為5位，AD9957并行數(shù)據(jù)輸入口為18位，留給串行IQ鏈路多相時變系數(shù)濾波器的有效數(shù)據(jù)位數(shù)是1位符號位+13位數(shù)據(jù)位，相乘后得到的數(shù)據(jù)進行3級流水線加法器得到累加輸出值。

圖6 串行IQ鏈路多相時變系數(shù)濾波器的結(jié)構(gòu)

3.2 可變參數(shù)基帶處理模塊數(shù)字邏輯設計

可變參數(shù)基帶處理模塊主要實現(xiàn)的是信號調(diào)制方式、調(diào)制階數(shù)、成型濾波器插值倍數(shù)和余弦滾降系數(shù)的可變調(diào)度，具體結(jié)構(gòu)圖，如圖7所示。

NIOSⅡ軟核通過AVALON-MM總線傳輸?shù)刂?、?shù)據(jù)、寫使能給從端的基帶控制模塊IP核從而改變IP核內(nèi)部對參數(shù)控制的數(shù)據(jù)，需要一個地址位映射到兩個寄存器：控制數(shù)據(jù)寄存器和FIR抽頭系數(shù)寄存器，為了節(jié)約FPGA內(nèi)部存儲資源，F(xiàn)IR抽頭系數(shù)查詢表放在SDRAM中，通過Nios的讀取并傳輸?shù)交鶐KIP

圖7 可變參數(shù)基帶處理模塊結(jié)構(gòu)圖

的FIR抽頭系數(shù)寄存器中。數(shù)據(jù)總線設計為16位，在控制數(shù)據(jù)寄存器中[1:0]位傳輸插值倍數(shù)控制參數(shù)，[3:2]位傳輸調(diào)制階數(shù)參數(shù)，第4位傳輸調(diào)制模式控制參數(shù)，第5位傳輸DAC和正交調(diào)制切換選擇參數(shù)，[15:6]傳輸余弦滾降系數(shù)參數(shù)；在FIR抽頭系數(shù)寄存器中，數(shù)據(jù)總線只需要14位，每次傳輸一個抽頭系數(shù)，因為抽頭系數(shù)是對稱的，因此直到抽頭系數(shù)通過總線傳輸個數(shù)達到65個即停止。每次參數(shù)更改都需要軟復位清零一下計數(shù)控制單元的計數(shù)值，調(diào)制方式可選值為QAM和ASK，使用QAM調(diào)制時Q路移位寄存器處于工作狀態(tài)，使用ASK調(diào)制時，Q路移位寄存器被旁路，Q路相對I路間隔輸出零值，即關(guān)閉AD9957的DDS余弦輸出功能。調(diào)制階數(shù)可選為4QAM/2ASK、16QAM/4ASK、64QAM/8ASK、256QAM /16ASK，插值倍數(shù)可調(diào)階數(shù)為2、4、8、16，選擇這個階數(shù)的原因是這些抽頭系數(shù)的公倍數(shù)是最高可插值倍數(shù)16，這樣可以節(jié)省抽頭階數(shù)的存儲，16倍插值的成型濾波器抽頭系數(shù)包含了2、4、8倍插值的抽頭系數(shù)，插值倍數(shù)控制模塊通過控制計數(shù)單元的步進梯度從而間隔選擇不同的子濾波器。余弦滾降系數(shù)可選范圍為0.3～0.5，更新余弦滾降系數(shù)FIR抽頭系數(shù)查詢表切會換到相對應滾降系數(shù)的FIR抽頭系數(shù)RAM，并導入到多相時變系數(shù)插值濾波器組中。計數(shù)控制單元的最低位映射到IQ鏈路的切換，高位映射到插值子濾波器的移位更新及用戶數(shù)據(jù)星座映射的更新。

4 調(diào)制信號發(fā)生器仿真和實驗結(jié)果

為了提高基帶處理模塊的最高運行頻率，需要對AD9957并行輸入端口進行時序約束，根據(jù)AltiumDesigner畫出的PCB線路統(tǒng)計出18根并行總線和1根時鐘總線的PCB傳輸時延，在介電常數(shù)為4.0的條件下，單位走線時延[5]為0.18 ns/1 000 mil，在TimeQuest Timing Analyzer中約束輸出最大和最小延時：

Output Max Delay=PCB Max Delay+

AD9957 Setup Time+Clock Max Delay+

(Latch Time-Launch Time)=

0.35 ns+1.75 ns+0.24 ns-5 ns=-2.66 ns

Output Min Delay=PCB Min Delay-

AD9957 Hold Time+Clock Min Delay+

(Latch Time-Launch Time)=0.18 ns+0 ns+

0.23 ns-5 ns=-4.59 ns

(6)

仿真后得到PDCLK最高時鐘118.96 MHz，這個結(jié)果已經(jīng)可以滿足AD9957內(nèi)部絕大部分的可編程插值倍數(shù)值，其中PDCLK需要配置為極性反轉(zhuǎn)，上升沿驅(qū)動，下降沿鎖存。

圖8和圖9分別是用戶數(shù)據(jù)經(jīng)過脈沖成型后的Modelsim基帶波形圖和SignalTapⅡ基帶波形圖，仿真條件：余弦滾降系數(shù)0.5，插值倍數(shù)8，雙邊碼元截取長

圖8 Modelsim基帶波形仿真圖

圖9 SignalTapⅡ基帶波形仿真圖

找到波形，用戶輸入數(shù)據(jù)采用了步進遞增符號碼。

圖10和圖11是硬件實驗結(jié)果圖，基帶配置條件不變，AD9957配置為內(nèi)部時鐘1GB，CCI+HB插值倍數(shù)為20倍，反sinc濾波器和反CCI濾波器使能，圖10是AD9957配置為DAC模式時輸出基帶波形，PDCLK時鐘為50MB，圖11是配置為QDUC模式時輸出的16ASK調(diào)制波形，PDCLK時鐘為100 MB，載波設置為50 MB。

圖10 基帶波形實驗結(jié)果

圖11 調(diào)制波形實驗結(jié)果

由圖可見，基帶波形實驗結(jié)果和仿真結(jié)果基本相同，調(diào)制波形的包絡很好地還原了基帶波形的形狀。

5 結(jié) 語

本文闡述了升余弦濾波器的時域頻域特性，并充分分析了FIR升余弦濾波器的余弦滾降系數(shù)、雙邊截斷碼元數(shù)、插值倍數(shù)的取值關(guān)鍵要點，設計了串行IQ鏈路多相時變系數(shù)濾波器，提高了系統(tǒng)的時鐘資源、邏輯資源及DSP資源的利用率。最后利用FPGA可編程邏輯完成了可變參數(shù)基帶處理模塊，提升了系統(tǒng)的靈活性，并集成為AVALON總線規(guī)范的IP核，增加了該模塊的可移植性，該設計使用了較低的資源實現(xiàn)了軟件可控的高階調(diào)制信號發(fā)生器。

參考文獻(References)：

[1] 吳彬彬，郭業(yè)才. 高階調(diào)制信號發(fā)生器設計與實現(xiàn)[D]. 南京：南京信息工程大學碩士學位論文, 2015.

[2] Pong P. CHU，基于NiosⅡ的嵌入式SoPC系統(tǒng)設計與Verilog開發(fā)實例[M]. 北京:電子工業(yè)出版社,2015.

[3] 梁婷婷，遲永剛，劉軍軍，等.一種具有任意倍升采樣的高速率脈沖成型算法[J].哈爾濱：哈爾濱工業(yè)大學學報，2015,47(11):53-58.

[4] 付文君，姜景山，王拴榮，等. 一種可變插值比的脈沖成型數(shù)字濾波器FPGA設計[J]. 遙感技術(shù)與應用, 2006,21(5):445-450.

[5] 吳厚航. 愛上FPGA開發(fā)特權(quán)同學和你一起學NIOSⅡ[M].北京:北京航空航天大學出版社,2011.

[6] 胡俊杰.衛(wèi)星數(shù)據(jù)傳輸高速調(diào)制器關(guān)鍵技術(shù)研究[D]. 北京:中國科學院研究生院博士學位論文，2011.

[7] Gentile K. The Care and Feeding of Digital Pulse-Shaping Filters[DB/OL].(2002-04-30)[2017-11-10].HTTP://www.rfdesign.com.

[8] 杜 勇. 數(shù)字調(diào)制解調(diào)技術(shù)的Matlab與FPGA實現(xiàn)[M].北京:電子工業(yè)出版社, 2015.

[9] 郭業(yè)才，章 濤，袁 濤. 衛(wèi)星信道模擬器數(shù)字下變頻系統(tǒng)設計[J]. 科學技術(shù)與工程, 2014，24(1):54-59.

[10] 劉華僑，林靜然.64QAM調(diào)制解調(diào)關(guān)鍵技術(shù)研究與實現(xiàn)[D].成都: 電子科技大學碩士學位論文,2013.

[11] 于偉龍，宗成閣.基于AD9957的信號發(fā)生器設計[D]. 哈爾濱: 哈爾濱工業(yè)大學碩士學位論文,2014.

[12] 譚龍龍，宗成閣.列車防撞雷達信號源開發(fā)設計[D]. 哈爾濱: 哈爾濱工業(yè)大學碩士學位論文, 2013.

[13] 任 哲. 嵌入式實時操作系統(tǒng)uC/OS-Ⅱ原理及應用[M]. 2版.北京航空航天大學出版社,2013.

[14] 李 敏,王 孜, 邵奕強,等.自適應低通濾波器的設計與實現(xiàn)[J].實驗室研究與探索,2014,33(8):89-93.

[15] 樊開陽,杜小峰, 楊紅兵. 基于FPGA的FIR濾波器設計方案優(yōu)化[J].實驗室研究與探索,2014,33(5):91-95.