数字射频存储(DRFM)系统是将射频信号保存在数字存储器中,在需要的时候将数字信号恢复为射频信号并输出。它在射频信号模拟、电子对抗等领域有着广泛的应用。随着雷达技术的迅猛发展,瞬时带宽达GHz级的现代宽带微波成像雷达已经开始走向实用。据报道,美国JPL研制中的工作在675 GHz的太赫兹成像雷达带宽可达到28 GHz。针对超宽带雷达,有必要研究瞬时带宽相当的DRFM系统,更高带宽将是DRFM系统的一个重要发展方向。
由于传统的单通道DRFM系统的带宽指标取决于微波前端带宽、基带A/D、D/A的采样速率、存储/信号处理芯片的速率,导致宽带DRFM系统的瞬时带宽指标严重依赖关键器件。以瞬时带宽2GHz的宽带DRFM系统为例,根据奈奎斯特基带采样定理再考虑一定的裕量,必须采用至少5GHz采样率的高速ADC和DAC。这些关键器件目前大多依靠进口,其中某些器件是对中国禁运的,且随着瞬时带宽的提升,甚至国外器件也不能满足系统指标的要求。从物理机理上看,DRFM系统性能的提升根本问题在于材料、工艺领域的技术提升。如目前微波前端采用的半导体工艺制程技术其切换(On-Off)频率在100~120 GHz范围内,该技术能实现的前置放大器最高模拟带宽是30 GHz,新的InP DHBT工艺可使切换频率扩展到350 GHz以上,且电流密度范围更大,可实现最高模拟带宽达70 GHz以上。
在获得高性能器件受限的前提下,提高DRFM系统的瞬时带宽只能在新的架构及突破性方法上进行。传统采用正交双通道[1]、多通道[2]架构来提高DRFM系统的瞬时带宽。缺点是正交双通道DRFM架构只能将瞬时带宽提高一倍,文献[2]提出的多通道DRFM架构不能对通道间的不一致性进行相位级补偿。
多通道架构由于系统较单通道复杂,多硬件带来的系统误差导致系统的底噪增加、信号严重畸变。为实现通道间相位级补偿,通常系统共用一个参考时钟实现通道间相参。该技术已在示波器领域得到广泛利用,如力科、安捷伦等公司通过DBI[3]、realedge[4]等多通道技术成功将基于SiGe工艺的示波器瞬时带宽提升至30 GHz以上,性能与InP工艺的示波器接近。文献[5]提出ATI(asynchronous time interleaving)架构,它利用2片采样率为100GS/s的ADC进行并行采样,实现了从DC到70 GHz的瞬时带宽;文献[6]采用压缩感知提出MWC(modulated wideband converter),该架构利用4片采样率为280 MS/s的ADC实现了对带宽为2 GHz频谱的实时重构,重构信号的瞬时带宽不能超过120 MHz。
对于DRFM系统而言,其架构比示波器更复杂,上述多通道技术不能直接应用到DRFM系统中。此外,新的架构及其失真补偿技术作为各个厂家的绝对机密,外界知之甚少。本文提出了一种基于多通道频域交织的DRFM架构和系统失真补偿算法,实验证明了该系统架构和补偿算法的有效性。
1 多通道频域交织架构与性能分析本文提出的多通道频域交织DRFM架构通过增加微波通道在频域对信号进行交织复用处理,用小带宽及低速器件构建了大带宽的DRFM系统。一个典型的N通道架构如图 1所示,其关键点频谱示意图如图 2所示。实现流程为:将输入信号按频率高低分为N路分别进行处理,在基带通过数字信号处理补偿后,再将N路合为一路输出,获得更大的带宽。
下面对该架构的最大通道数量N及系统极限瞬时带宽Blim进行分析。为方便分析,假设如下:
1) 宽带雷达信号经DRFM转发后,理想重构信号经雷达脉冲压缩信号处理应表现为sinc型函数,当冲激相应宽度(IRW)小于1.5个系统分辨率ρ、且峰值旁瓣比(PSLR)低于-20 dB时,认为重构信号的逼真度达到DRFM系统指标要求;
2) DRFM系统有N个通道,每个通道的带宽均为B MHz,理论上DRFM系统的最大瞬时带宽为BN,对应的系统分辨率ρ为c/(2BN);
3) 每个通道的幅频特性及延时特性均能补偿,补偿后通道及通道间的幅频特性误差不影响DRFM系统对信号的重构,但需要考虑通道间延时不一致误差对信号的重构。按照距离测量噪声误差公式可知,单个通道延时测量误差var{s}为:
$ {\mathop{\rm var}} \{ s\} = \frac{c}{{2B\sqrt {{R_{{\rm{sn}}}}} }} $ | (1) |
考虑重构信号脉冲压缩后的冲激相应宽度(IRW)小于1.5个系统分辨率,即每个通道的延时测量误差var{s}应小于0.75ρ,将其带入式(1),可得通道数N和信噪Rsn比之间的关系为:
$ N < \frac{{3\sqrt {R{\rm{sn}}} }}{2} $ | (2) |
系统极限瞬时带宽Blim和信噪Rsn比之间的关系为:
$ {B_{\lim }} < \frac{{3B\sqrt {R{\rm{sn}}} }}{2} $ | (3) |
由式3可知,系统的理论带宽极限Blim与每个通道的带宽B和通道测量信号功率成正比。
2 补偿算法系统采用多通道架构,多出的硬件带来的系统误差导致系统的底噪增加。且功分器、变频器、功合器等器件本身是宽带的,其幅频特性及群延时特性都做不到理想状态,通道间存在相对时间及幅频不一致。此外,还有频段交界点信号失真、系统时变等误差,必然会加剧系统的总失真,甚至导致系统无法正常工作。因此,必须对各种误差因素分别进行补偿。
2.1 通道内幅频特性及群延时特性误差补偿由于通道内的幅频特性不平坦和群延时不一致等引起的误差,会使经DRFM系统转发信号的带内相频波动增加,导致重构信号质量降低。本文通过校准、补偿相结合的办法[7]在基带对采样数据进行失真修正,分别对每个通道内的幅频特性及群延时特性误差进行抑制,从而保证重构信号与原始信号接近。
当系统在理想的情况下,复检波信号dn是理想的复正弦波,实际得到的I、Q离散样点值cn不是纯正的复正弦波,如何求得一组相应的复数gn,满足:
$ {c_n} \times {g_n} = {d_n} $ | (4) |
将求得的gn与信号处理获得的cn直接相乘,就构成了在信号处理机对系统的补偿。由cn求gn必先求dn。根据样本cn对dn的曲线做拟合运算,求幅度M和初相${\varphi _0}$,圆频率ω=2πkτ,将cn化为极坐标形式${c_n} = {A_n}\exp (j{\varphi _0})$,其最小二乘估计为:
$ \left\{ \begin{array}{c} \hat M = \frac{1}{N}\sum\limits_{n = 1}^N {{A_n}} \\ {{\hat \varphi }_0} = \frac{6}{{N\left( {N - 1} \right)}}\sum\limits_{n = 1}^N {{\varphi _n}\left( {\frac{{2N + 1}}{3} - n} \right)} \\ \hat \omega = \frac{6}{{N\left( {N - 1} \right)\left( {N + 1} \right)}}\sum\limits_{n = 1}^N {{\varphi _n}\left( {2n - N - 1} \right)} \end{array} \right. $ | (5) |
由式(4)和式(5)可得:
$ {g_n} = \frac{{\hat M{\rm{exp}}\left[ {{\rm{j}}\left( {\hat \omega {t_n} + {{\hat \varphi }_0}} \right)} \right]}}{{{c_n}}} $ | (6) |
$\left\{ \begin{array}{c} \left| {{g_n}} \right| = \frac{{\hat M}}{{{A_n}}}\\ \arg ({g_n}) = \hat \omega {t_n} + {{\hat \varphi }_0} - {\varphi _n} \end{array} \right. $ | (7) |
相邻通道间幅频特性不一致,导致合成的波形可能出现相位断裂、波形混乱。通过分析发现,新架构对幅度差不敏感,通道之间的延时差(相位差)是影响重构信号质量的最大因素。因此,通道间的误差补偿重点是对各通道的延时差(相位差)补偿。
各通道的延时量测量原理如下:通过数字信号板板载的数模转换器(DAC)发送一个斜率为k的宽带LFM脉冲信号,该信号通过反馈支路到达射频输入端,经合路器、带通滤波、下变频、低通滤波后,被DAC采集下来。在信号处理器内部,对原始发射信号和被采集信号进行STRETCH去斜处理。由于两个信号之间存在延时差,因此STRETCH会得到一个单频输出Δf,该频率即可以基本代表输入-输出的延时量:Δt=Δf/k。为了得到精确的延时量,可采用插值、递归方法[8, 9]精确估计延时量。
2.3 频段交界点信号失真补偿对于宽带跨通道信号,在频段交界点将不可避免地产生信号失真。文献[10]提出基于分析/综合滤波器组能在一定约束条件下实现全子信道的准确重构(PR)跨道信号。目前最好的结果是Agilent公司的示波器DSAX95004Q[5],该系统实现了频段交界带宽小于30 MHz,通过补偿后重构信号在频段交界处平坦度优于3 dB。本系统通过补偿,频段交界处的平坦度优于6 dB,测试结果见4.1节。
3 系统及关键模块假设每个DRFM通道的带宽为1 GHz,按照本文提出的频域交织DRFM架构,并行10个通道,理论上可以构建瞬时带宽高达10 GHz的DRFM系统。在实际中,2个带宽700 MHz的DRFM通道构建了瞬时带宽1.4 GHz的C波段DRFM系统。该系统由微波模块、数字基带板、CPCI机箱构成,实物如图 3所示。
微波模块功能有信号功分、上下变频、信号功合,并为系统提供时钟信号。射频信号f0±700 MHz分成两个通道完成放大变频滤波,然后输出中频300~1 000 MHz的信号。发射链与接收链相反,通过两个接口输入中频信号,然后滤波变频放大,合路为一个通道后输出,微波模块功能如图 4所示。
数字信号处理集中在一块单板上实现, 2片超高速ADC用于对两路300~1 000 MHz频段内模拟输入信号的采集,2片超高速DAC用于输出两路宽带模拟信号。ADC和DAC参考时钟通过一个SMA射频连接器统一由微波模块提供,在板内进行调理和分配,以保证输入、输出信号的相参性。2片FPGA实现数据采集及存储控制、信号失真补偿等功能。数字信号处理板原理框图如图 5所示。
测试框图如图 6所示。利用泰克任意波形发生器AWG7122B直接产生C波段1.4 GHz带宽线性调频(LFM)信号,该信号输入被测DRFM系统,经多通道频率交织处理后,将重构信号输入泰克时频分析仪RSA6120A测试系统的瞬时带宽及平坦度。测试结果如图 7所示,频谱测试参数中心频率为f0 GHz,span为2 GHz,RBW为1 MHz。原始LFM信号的频域如图 7a所示,原始信号经过被测DRFM系统,重构信号的频域结果如图 7b所示,可见信号带宽依然为1.4 GHz,证明该系统瞬时带宽优于1.4 GHz。但原始信号的平坦度从3 dB恶化为22 dB,采用系统失真补偿后信号平坦度优于6 dB,如图 7c所示。
线性度是衡量DRFM系统的一项重要指标,反映了实际的LFM重构信号偏离理想线性频率的程度,计算式为:
$ \xi = \frac{{{{\left| {{f_e}(t)} \right|}_{\max }}}}{B} $ | (8) |
$ {f_e}(n) = \frac{{F(n)}}{{{\rm{2\pi }}\tau }} - nk{T_{\rm{s}}} + 2{\rm{\pi }}m $ | (9) |
由于通用测量仪器不具备测量超宽带信号线性度的功能,测量仪器采用中物院自行研制的多功能雷达信号处理模拟器测量方法参照文献[11]。测试框图如图 8a所示,实验如图 8b所示。
实验中,LFM信号带宽B为700MHz,时宽为50 μs,采样速率12.5MHz,根据式(1)和式(2)计算可得,补偿前的LFM重构信号的线性度为1.28×10-3,其调频误差曲线和脉冲压缩波形如图 9所示。经过误差补偿后,LFM重构信号的线性度降低为5×10-5,其调频误差曲线和脉冲压缩波形如图 10所示。比较图 9b、图 10b,可以看出补偿后脉压信号的主瓣变窄、旁瓣变低,重构信号的质量明显提高。
采用4.2节的测量方法,并将时间参数转换为距离参数,两通道合并后输出结果如图 11a所示。由于两通道的延时不一致,从图中可以看到两个峰,分别对应两个通道的延时距离,A通道的延时距离为78.08 m,B通道为78.88 m,两者相差0.8 m。经过通道间幅频不一致性误差补偿后结果如图 11b所示。比较图 10b、图 11b,可以发现双通道重构信号的主瓣较单通道的主瓣窄。这是因为单通道的重构信号带宽为700 MHz,双通道重构信号的带宽为1.4 GHz,理论上单通道的脉压信号主瓣比双通道的宽一倍。实验结果与理论分析相一致,表明本文提出的系统架构可行,补偿算法有效。
本文通过DRFM架构上的创新,设计了补偿算法,弥补了系统复杂性带来的失真,运用目前可得到的国产器件(技术指标不需要很高)达到高端产品的系统指标,大大扩展了DRFM的瞬时带宽。此外,实验中发现温度、通道中任何部件的更换、甚至连接器接插的松紧程度等都将导致补偿结果的恶化。下一步工作中需要对系统时变特性进行深入分析,研究自适应实时校正补偿算法。
[1] | 赵国庆. 雷达对抗原理[M]. 西安: 西安电子科技大学出版社, 2003. ZHAO Guo-qing. The principle of radar ECM[M]. Xi'an: Xi'an Electronic and Science University press, 2003. |
[2] | 刘忠, 王国玉, 王雪松. 数字射频存储器(DRFM)结构分析及改进[C]//中国电子学会电子对抗分会第十四届学术年会. 洛阳: 中国电子学会, 2005: 428-432. LIU Zhong, WANG Guo-yu,WANG Xue-song. Digital radio frequency memory (DRFM) analysis and structure improvement[C]//Chinese Electronic Warfare Division of the Fourteenth Annual Meeting. Luoyang: Chinese Institute of Electronics, 2005: 428-432. |
[3] | LECROY Corporation. The interleaving process in digital bandwidth interleaving (DBI) scopes[EB/OL]. [2014-05-05]. http://teledynelecroy.com/files/whitepapers. |
[4] | Agilent Corporation. Agilent infiniium 90000 Q-series oscilloscopes[EB/OL]. [2014-05-05]. http://www. agilent.com/find/9000Q-Series. |
[5] | Tektronix Corporation. The Quest for higher real-time oscilloscope bandwidth[EB/OL]. [2014-05-05]. http:// download.tek.com /secure/. |
[6] | MISHALI M, ELDAR Y C, DOUNEAEVSKY O, et al. Xampling: Analog to digital at sub-Nyquist rates[J]. IET Circuits Devices Syst, 2011, 5(1): 8-20. [7] 万永伦, 姒强, 吕幼新, 等. 基于信号产生器的幅相误差分析与校正[J]. 现代雷达, 2007, 29(12): 101-108. WAN Yong-lun, SI Qiang, LÜ You-xin, et al. Analysis and calibration of amplitude and phase erros in signal generator[J]. Modern Radar, 2007, 29(12): 101-108. |
[8] | ABOUTANIOS E, MULGREW B. Iterative frequency estimation by interpolation on Fourier coefficients[J]. IEEE Trans Signal processing, 2005, 53(4): 1237-1242. |
[9] | 江舸, 陶荣辉, 李钒, 等. 基于FPGA实现的FFT插值正弦波频率估计[J]. 信息与电子工程, 2009, 7(6): 569-572. JIANG Ge, TAO Rong-hui, LI Fan, et al. Frequency estimation of sinusoid waveform by interpolated FFT based on FPGA implementation[J]. Information and Electronic Engineering, 2009, 7(6): 569-572. |
[10] | JIANG J Z, SHUI P L, ZHANG Z J. Design of oversampled DFT-modulated filter banks via modified Newton's method[J]. IET Signal Processing, 2011, 5(3): 271-280. |
[11] | 万永伦, 姒强, 汪学刚. 超宽带线性调频信号线性度的测量方法[J]. 电子测量与仪器学报, 2007, 21(4): 55-58. WAN Yong-lun, SI Qiang, WANG Xue-gang. Linearity measure-ment method for ultra-wideband linear frequency modulated signal[J]. Journal of Electronic Measurement and Instrument, 2007, 21(4): 55-58. |