基于FPGA的超宽带数字波束形成技术

doi:10.3969/j.issn.1003-3114.2022.01.022
引用格式:顾明超,李春晓,边疆,等.基于FPGA的超宽带数字波束形成技术[J].无线电通信技术,2022,48(1):173-179.[GUMingchao,LIChunxiao,BIANJiang,etal.TechniqueofUWBDigitalBeam⁃FormingforFPGA[J].RadioCommunicationsTech⁃nology,2022,48(1):173-179.]
基于FPGA的超宽带数字波束形成技术
顾明超1,2,李春晓1,2,边㊀疆1,2,张汉卿1,2
(1.中国电子科技集团公司第五十四研究所,河北石家庄050081;
2.河北省电磁频谱认知与管控重点实验室,河北石家庄050081)摘㊀要:随着高速ADC和高性能FPGA芯片工艺技术的飞速发展,芯片的小型化和集成化促进了电子战系统瞬时覆盖带宽由几十兆赫兹上升到几百兆赫兹甚至到吉赫兹量级,属于超宽带范畴㊂为了获得更远的探测距离性能,具备同时对多个辐射源侦测的能力,需要瞬时形成多个数字波束对目标区域密集覆盖㊂针对宽带时域校准与合成处理资源占用多㊁形成波束少的缺点,在FPGA中采用整数点移位与相位校准滤波方法实现了通道
校准,通过整数点移位与分数时延滤波处理降低了合成滤波器的阶数㊂对瞬时带宽内多个等间隔单音信号的波束指向进行仿真,通过波束图验证了宽带合成各频率指向的一致性㊂提出了一种处理实时㊁资源占用少的频域处理方法,结合FPGA芯片特性进行了硬件化实现,对合成前后频谱进行了对比,验证了合成增益㊂
停车场闸机
关键词:阵列信号处理;FPGA;超宽带数字多波束形成;电子侦察
中图分类号:TN911㊀㊀㊀文献标志码:A㊀㊀㊀开放科学(资源服务)标识码(OSID):
文章编号:1003-3114(2022)01-0173-07
TechniqueofUWBDigitalBeam⁃FormingforFPGA
GUMingchao1,2,LIChunxiao1,2,BIANJiang1,2,ZHANGHanqing1,2
(1.The54thResearchInstituteofCETC,Shijiazhuang050081,China;
2.HebeiKeyLaboratoryofElectromagneticSpectrumCognitionandControl,Shijiazhuang050081,China)Abstract:Withthefastdevelopmentofhigh⁃speedADCandhigh⁃performanceFPGA,miniaturizationandsystemintegrationforchipsmakeinstantaneousbandwidthofelectronicwarfaresystemincreasedfromtensofMHztohundredsoreventhousandsofMHz,whi
chcanbeclassifiedtoultra⁃widebandcategory.Comparedwithtraditionalelectronicreconnaissancesystem,inordertoachievedetection/identificationofmultipleradiationsources,multipledigitalbeamsneedtobeformedtocovertargetarea.High⁃performanceFPGAhasabundantlogicoperationsandstorageresourceson⁃chip,whichishighlysuitbleforrealizationofmulti⁃beamforminginarraysignalprocessing.Inordertoreducetheresourceoccupancyrateofwidebandtime⁃domaincalibrationandsynthesisprocessingandformmorebeams,integerpointshiftandphasecalibrationfilteringmethodareusedtorealizechannelcalibration,andintegerpointshiftandfractionaldelayfilterprocessingareadoptedtoreducesynthesisfiltercoefficientsinFPGA.Throughsimulatingthedirectivityofmultipleequally⁃spacedtonesignalswithintheinstantaneou
sbandwidth,theconsistencyofeachfrequencydirectionofthewidebandsynthesisisverifiedbythegeneratedbeampattern.Acalibrationandsynthesisintegratedfrequencydomainprocessingmethodwithreal⁃timepro⁃cessingandlessresourceoccupancyisproposed,andrealizedbyhardwareincombinationwiththecharacteristicsoftheFPGAchip.ThenthespectrumbeforeandaftertheDBFiscompared.Andthecombininggaincanbeverified.
Keywords:arraysignalprocessing;FPGA;UWBmultidigitalbeam⁃forming;electronicreconnaissance
收稿日期:2021-10-20
基金项目:国家科技计划专项经费(2016QY10W1802)
FoundationItem:SpecialFundsofNationalScienceandTechnologyPlan(2016QY10W1802)
0㊀引言
相控阵技术是近年发展迅速的前沿热点技术之
一,早先应用于雷达探测技术,多为窄带阵列处理体
制[1],典型带宽值为几兆赫兹至几十兆赫兹㊂随着技术推广,目前已经逐步应用到电子对抗系统中,瞬
时带宽也要求到几百兆赫兹甚至更高㊂2020年12月中国工程院发布的信息与电子工程领域技术前沿,相控阵技术在该领域10项前沿工程技术中占据两席,足以说明相控阵体制的技术先进性,对其关键技术工程化研究应用具有一定紧迫性,本文主要介绍其关键技术之一的超宽带数字波束形成技术的原理与其在FPGA中实现过程㊂
相控阵技术从有源放大和波控单元的位置区
分,可分为无源和有源相控阵;从波束形成方式区
分,可分为模拟波束合成㊁数字波束合成以及模数混
合波束合成;从技术发展角度看,相控阵从无源发展
到有源,从模拟发展到数字架构㊂数字波束形成是
相控阵领域中的一个关键技术,其本质是一种数字
化空域滤波,增强特定方位入射到阵列的信号,提高
信噪比,抑制其他方向的干扰与噪声,为之后信号处
理提供必要条件㊂在电子对抗领域,若能形成多个
指向可变的宽带数字波束密集覆盖敏感区域,通过
多波束并行侦察方式,势必能提高电子侦察系统的
信号截获概率㊂
电子对抗中对宽带没有明确定义,一般认为信
号带宽与中心频率之比小于1%时,为窄带信号;大
于1%小于25%时,为宽带信号;若带宽与中心频率
的比值大于25%,为超宽带(UWB)信号[2]㊂在电子对抗领域,由于为非合作接收方式,不明确信号调制类型和带宽以及持续时间,需通过高速ADC采集瞬时覆盖目标频段,再通过宽带侦察和信道化检测等措施对信号参数进行测量㊂目前瞬时带宽典型应用范围为300MHz及以下,瞬时覆盖500MHz到GHz带宽的需求也将日益增多[3]㊂
Xilinx公司Virtex⁃7系列的大容量FPGA是一类性价比高且应用较广的芯片㊂这个系列FPGA有数量众多的逻辑/存储资源以及较为丰富的乘法器和数量众多的高速串行总线接口,非常适合于超宽带数字多波束的技术验证㊂
1㊀数字波束形成实现原理
1.1㊀数字波束合成原理
本文主要讨论信号接收波束形成,不涉及发射波束㊂早期的波束合成往往采用模拟移相的方式对
多路模拟信号进行移相后合路,由于模拟移相器在
宽频带中具有 非散 特性,对覆盖频带中各频
点移相值一样,应用于窄带信号的合成效果较为
理想㊂但是对于电子侦察而言,瞬时覆盖带宽高
达几百MHz甚至GHz,由于天线孔径效应的存在,
导致宽带波束方向与理想方向产生偏差,瞬时带
宽越大,天线的孔径效应越明显[4]㊂模拟信号移相合成方式已经不能满足宽带合成需求,必须采
用宽带数字校准合成技术实现对多阵元宽带数据
流的处理㊂
在阵列信号接收过程中,辐射源发出的目标信
号在空间传输后到达接收阵列天线的各阵元,由于
天线阵元的位置差异导致其接收信号的波程差异,
因此各阵元的输出信号存在一定的相对时间延迟㊂
对各阵元接收信号的时延做出补偿,让经过补偿后
的各个输出信号在所期望的空间方向上幅度为同向
相加,则此时可以最大化指定方向的波束输出信号
幅度,并且使其他方向上的波束幅值相应变小,由此
对非期望方向上的无用信号起到了抑制作用,用时
域滤波处理方式达到了空域滤波的效果,这也是数三方通话>实心锥形喷嘴
字波束合成的一般性原理㊂
下面以均匀线阵㊁单音信号输入为例进行说明,
相邻两个阵元间隔为d,如图1所示
图1㊀信号入射均匀线阵示意图
Fig.1㊀Schematicdiagramofsignalincidentuniformlineararray
假设阵元1为相位参考,来波方向与法线方向
夹角为θ,c为光速,各阵元的时延表达式为[5]:
τi=d(i-1)sinθ
c,i=1,2,...N㊂(1)则各阵元相位差如式(2)所示:
Δφi=2πfcτi=2πfcd(i-1)sinθ
c=
2πd(i-1)sinθ
λ
,i=1,2,...N㊂(2)若采样率为fs,则整数时延就近取整后,整数延迟值计算如式(3)所示:
Lint=round(τi㊃fs)㊂
(3)
分数时延计算如式(4)所示:
Lfraction=τi㊃fs-round(τi㊃fs)㊂(4)式(1) (4)的推导隐含了一个重要前提,即各阵元接收到数字化过程的通道一致性是理想的,但实际情况并非如此,需要校准源产生扫频或梳状谱信号作为阵列接收模块的输入激励,对通道一致性进行校准,再按上述分析过程进行合成处理㊂
通道间不一致包括在法线方向时,各阵元接收信号时延不一致和通道间非线性相位差异,时延不一致可分解为整数采样点和分数点延迟,整数采样点可通过采样点校准在FPGA中对数据缓存或移位寄存器实现,分数点延迟和非线性相位差异可以与波束合成功能的分数点延迟功能进行合并,通过频域实现,这种组合实现方式有利于降低硬件实现的复杂度和资源消耗㊂
1.2㊀宽带数字波束合成实现过程
一般阵列处理系统包含阵列天线㊁模拟信道㊁数字处理和软件控制席位四部分㊂波束控制由软件席位控制模拟信道和数字部分电路完成㊂简化框图如图2所示
图2㊀宽带数字波束形成原理图
Fig.2㊀Blockdiagramofwide⁃bandDBF在图2中,若ADC采样率足够高,采用射频直接采样体制接收覆盖目标频段,各阵元接收的多通道射频信号经过滤波放大后进入高速ADC完成模数转换,通过实采样校准和数字延迟滤波完成波束合成功能,或者实采样信号通过宽带变频滤波后形成多相结构的零中频数据流,再经过宽带复校准和数字复合成滤波计算后,完成宽
带数字波束合成功能㊂
宽带侦察系统中,除了要求形成宽带波束,往往还采用宽带侦察引导窄带测向的方式,需要形成多个窄带波束进行控守,只需按照阵列模型和频点计算合成系数,再与按中心频点计算的通道间相位校准系数相乘,通过窄带滤波后的乘加运算完成信号合成功能㊂
2㊀宽带数字波束合成实现过程
波束合成的实现可分为时域法和频域法,无论哪种方法,预处理包含的采样同步㊁采样点校准都是必须的,都是DSP通过回读FPGA内同步存储的多路数据,经过校准和合成算法,将计算后系数写入
FPGA中[6]进行后续运算㊂时域法和频域法均可实现宽带数字波束,处理方法的选择主要取决于波束数量㊁硬件资源消耗与系统成本㊂
2.1㊀时域合成法
宽带波束形成时域实现常见的有FIR滤波器实现㊁基于Farrow的分数点延时法和微波光子技术的时延方法㊂Farrow结构[7]可实现任意分数点延迟,延迟改变只需更改输入参数即可,延迟精度取决于相数和阶数,对宽带多相数据而言,乘加器使用数量较大,工程化优势不明显㊂微波
光子技术[8]通过光延迟实现真延时功能,延迟精度受环境因素影响较大,目前技术成熟度暂不满足工程化要求㊂
宽带时域波束合成最常见做法是通过FIR滤波器组实现的,利用FIR滤波器实现各阵元不同延时后的加权㊂其设计思路是:若需在期望方向形成指向波束,选择带宽内一定数量频点来进行波束设计,得出在这些特定频点上的加权值,也就是得到每个阵元的幅度权与相位权㊂设计一组滤波器,使每个滤波器的幅相响应分别在这些频点上,与各阵元的幅度加权和相位加权近似相同㊂换言之,就是设计一组FIR滤波器,用其幅频响应和相频响应分别拟合各阵元的幅度权与相位权㊂时域合成框图如图3所示
图3㊀时域合成框图
Fig.3㊀Blockdiagramoftimedomain
DBF
模拟信号数字化后,通过FIR滤波器加权后相加,形成波束输出㊂在不考虑通道一致的情况下,为降低FIR滤波器阶数,只需要滤波器实现分数点延迟,整数点延迟由FPGA内部缓存或移位寄存器完成㊂
理想的分数时延滤波器的冲激响应可表示为:hid(n)=sinc(n-D)㊂
(5)
当延迟不是正整数时,式(5)表示的滤波器是非因果的,若直接使用截断后的sinc函数来设计分数时延滤波器,其性能往往是不可接受的㊂为降低吉布斯效应的影响,时域加窗是常用的方法㊂
加窗后的冲激响应如式(6)所示:
h(n)=
W(n-D)sinc(n-D),0ɤnɤN
0,其他
{
,(6)
其中,理想冲激响应hid(n)被窗函数截断,窗长L=N+1㊂窗函数可选择海明窗㊁汉宁窗以及切比雪夫窗等㊂窗函数法计算量小㊁实时性高,但硬件资源受限导致滤波器长度较短时,该方法难以控制幅度响应误差㊂若波束指向变化时,需不断加载更新滤波器系数,造成波束数据间断,影响对目标信号接收效果㊂
设计举例:24阵元线阵,阵元间距0.1m,采样率1600Msample/s,瞬时覆盖500MHz带宽,要求形成3个波束,波束覆盖范围ʃ30ʎ㊂
校准滤波器系数64阶,时延滤波器系数32阶,满足设计要求㊂单片FPGA接收2路采样数据,8相200Msample/s结构,若采用实信号校准合成方法,校准滤波消耗乘法器为64ˑ8ˑ2=1024,合成滤波消耗乘法器为32ˑ8ˑ2ˑ3=1536,合计2560个乘法器,采用XC7VX690T可实现预期功能,但只能形成1个波束㊂由此可见,宽带时域处理方法资源消耗多,尤其是数据为多相结构时,资源按相数线性增长㊂2.2㊀频域合成法
时域合成处理占用FPGA内部资源较多,因此必须转换思路,选择一种处理实时资源消耗少的校准合成方法支撑宽带阵列处理数字多波束应用需求㊂
时域校准合成实质为数据流与时域系数卷积运算,需通过并行乘加实现㊂若将运算变换至频域处理,校准合成系数先进行卷积运算,补零后通过FFT运算变换至频域,频域系数与频域数据串行相乘,各路相加后再IFFT变换至时域,形成DBF数据流㊂FFT和IFFT消耗乘法器和存储资源较少,乘法器的使用数量较时域大为减少㊂频域合成FPGA实现如图4所示
终端准入系统㊂图4㊀频域合成FPGA实现框图
Fig.4㊀FPGAimplementationoffrequencydomainDBF
图4中,通过JESD204B总线接收的采样数据经过多相滤波形成四相零中频数据流,对多相数据流按N点等间隔进行划分,数据块顺序与FPGA中4路处理资源对应,经过数据缓存并串转换后进行FFT运算,并与DSP下发的频域系数相乘后与其他数据通道的对应支路求和,IFFT变换至时域,串并转换后,按照时间对应关系重新排列数据,形成连续的波束合成时域数据流㊂
3㊀FPGA实现验证
在超宽带波束合成FPGA实现过程中,遇到了多个技术实现问题㊂例如频率分辨率与自校源步进不匹配,FPGA处理时钟高造成时序收敛困难,数据传输量过大导致形成波束数量减少,频域合成后的时域数据断数,校准合成一体设计降低资源量,提升波束扫描连续性等FPGA实现问题㊂3.1㊀采样率变换提升处理与传输
在工程应用中,校准时理想状态是让校准频点精确落在整数谱线上㊂自校源产生的频率步进最小为500kHz,若按采样率进行计算,采样
800Msample/s,FFT点数256,此时频率步进为3.125MHz,不能与自校源进行适配㊂
在XC7VX690T中,当FPGA资源占用较多时,时钟速率运行至200MHz会造成时序收敛的困难㊂采样率为800Msample/s,一个波束的数据量为25.6Gbit/s,再加之64/66bit编码效率以及95%的传输效率,数据量为27.79Gbit/s,需通过一组4XAurora总线才能将一个波束传出,造成FPGA高速串行接口使用过多,减少了波束数据的传输,从而限制了波束形成数量㊂
为了解决以上三点问题,设计了基于多相结构的采样率变换模块,采样率由800Msample/s变换至640Msample/s,校准的频率步进由3.125MHz变为640/256=2.5MHz,适配了自校源频率步进特性㊂通过采样率变换后,FPGA内部逻辑处理时钟为640/4=160MHz,有效降低了FPGA电路时序收敛的难度,降低了FPGA技术开发难度㊂
通过变采样率运算后,FPGA内部形成一个波束数据量为22.232Gbit/s,采用2XAurora(线速率为12.5Gbit/Lane)便可完成波束数据传输,提高了产生极传输的数量㊂
3.2㊀IFFT变换后波束数据不连续
在频域处理时,通过仿真发现IFFT后的波束数据存在不连续现象,即当进行IFFT变换时,存在瞬态信息丢失问题㊂因此频域处理必须消除两次IFFT交界处相位不连续对合成造成的影响[9]㊂结合FPGA实现特点,采用1/2交叠运算的办法对相邻两次IFFT后时域数据交叠处进行去重处理㊂以单音信号为输入,通过仿真对比,频域FFT处理不作1/2交叠与作1/2交叠,输出时域数据的对比㊂光滑曲线是经过1/2交叠后的曲线,三角符
号曲线代表未经过1/2交叠的曲线,如图5所示
视频处理
图5㊀FFT计算1/2交叠对比图
Fig.5㊀1/2overlapcomparisonchartofFFTcalculation由图5不难看出,未经过1/2交叠处理的数据波形存在明显数据点周期性跳变现象,跳变周期与FFT点数一致㊂与不交叠处理相比,交叠处理会使FPGA波束形成模块的资源消耗翻倍,但与时域校准合成相比,频域处理资源消耗小,仍具明显优势㊂3.3㊀频域校准合成一体设计
传统的宽带时域数字波束形成多采用多相分解滤波器的结构,对通道校准滤波器和时延滤波器独立设计,通过滤波器级联的方式硬件实现,导致乘加滤波资源消耗过大,宽带波束形成个数较少,工程化使用受限㊂
如果将时域信号变换至频域处理,可在DSP处理器中,对校准系数与合成系数首先进行时域卷积运算,达到校准系数与合成系数一体设计的效果,再转化成频域系数进行乘加运算后转换为时域数据,通过这样的转化运算能有效降低FPGA内乘法器资源消耗,提升波束形成个数,工程化优势明显㊂3.4㊀提升波束扫描连续性
传统的数字波束合成采用延时滤波器实现,在FPGA中通过对FIR滤波器在线配置实现㊂当波束
指向发生变化时,需要实时对时延滤波器系数进行加载更新和复位,在更新系数时,输出的DBF存在时间上间断现象,影响对目标信号侦察控守㊂
频域处理时,合成系数是按块运算的,只要提前将指向系数写入FPGA缓存内,在FPGA的时序控制下,将系数写入波束合成模块,在不复位电路的情况下,对波束指向进行 捷变 ,切换速率为FPGA内部一个时钟周期㊂自动线喷涂
3.5㊀频率指向一致性验证
为了验证瞬时带宽内DBF后各频率指向一致性,500MHz内产生26个等间隔20MHz的单音信号进行仿真,在30ʎ指向时,得到各频率的波束图㊂由图6可以看出,所有频率分量的单音信号都指向了30ʎ,没有其他方向的波束峰值出现,从而验证了宽带下各频率指向一致性
图6㊀多音信号波束指向仿真
Fig.6㊀Beampointingsimulationofmultitone
signal

本文发布于:2024-09-24 05:26:19,感谢您对本站的认可!

本文链接:https://www.17tex.com/tex/3/214576.html

版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系,我们将在24小时内删除。

上一篇:机载天线综述
下一篇:天线测量报告
标签:波束   合成   信号   形成   数字   处理   实现
留言与评论(共有 0 条评论)
   
验证码:
Copyright ©2019-2024 Comsenz Inc.Powered by © 易纺专利技术学习网 豫ICP备2022007602号 豫公网安备41160202000603 站长QQ:729038198 关于我们 投诉建议