ISSN 
-
由于相控阵雷达在通道一致性、同时执行多功能、多目标跟踪等方面的要求，迫切需要数字化技术的应用，因而数字阵列雷达(digital array radar, DAR)成为研究热点。DAR是一种采用数字波束形成技术来形成接收和发射波束的全数字阵列电扫描雷达，具有良好的抗干扰能力和多工作模式切换能力[1]。
宽带数字阵列雷达(wide-band digital array radar, WB-DAR)除具有普通数字阵列雷达的优势外，还具有高距离分辨力、获取未知目标距离像、对目标进行分类和识别等特点，因而得到了越来越多的关注[2]。WB-DAR系统要实现大角度范围内的扫描，必须考虑孔径效应和孔径渡越时间的约束问题。因此，针对特定指向，通过模拟或数字方式补偿时延来实现宽带宽角度扫描是WB-DAR系统的必然选择[3-9]。波导或同轴电缆可进行模拟时延补偿，但会带来较大的体积和功耗，实现成本高且易受环境影响。改进的光纤时延补偿无法获得任意时延[5-8]。传统的数字时延补偿有加密采样、时域内插[9]等方法，但也无法获得任意时延且会造成处理数据量的增加。文献[10]针对UWB系统提出数字延迟线和分数时延滤波器相结合的数字时延方式来形成波束。文献[11-12]分别讨论了在WB-DAR系统子阵和阵元层面采用分数时延滤波器来形成宽带波束，这种数字时延波束形成方式适应WB-DAR系统的结构特点，具有良好的可实现性和宽带波束形成性能。基于这种方式的WB-DAR系统对每个阵元收到的回波信号采样进行正交解调、M倍抽取、幅相加权以后，各通道完成不同的整数倍和分数倍时延，合并形成I路和Q路基带信号。
由于目前专用的DDC(digital down convertor)器件尚难以满足WB-DAR输入带宽和数据率的要求，下变频和抽取处理过程通常选择大规模现场可编程逻辑门阵列FPGA(field programmable gate array)器件实现。其中，主要消耗乘法器和加法器资源的模块是由抗混叠滤波器和抽取模块构成的宽带信号抽取器和后继的可变分数时延滤波器。因此，减少接收通道资源消耗的关键在于优化设计方法，降低抗混叠滤波器和分数时延滤波器的阶数来满足指定的性能指标条件。
当前多级抽取框架下常采用积分梳状(cascade integrated comb, CIC)滤波器或半带(half band, HB)滤波器来降低滤波器阶数[13]。对于宽带数字接收通道，在有限的采样率下处理大瞬时带宽的信号，需要在抽取倍数和级数都不高的条件下，在较大的相对带宽范围内保证带内平坦度，提供陡峭的过渡带。但CIC滤波器需要较高级数级联才能提供良好的带外抑制，并且为了保证带内平坦度，还需要设计补偿滤波器[14]，增加了设计复杂度；HB滤波器则不易实现陡峭的过渡带。因此，在工程实现时，宽带数字接收通道常采用具有线性相位特性的多级通用FIR滤波器来提供理想的频响特性。在此基础上，本文提出了一种新的优化设计方法，与传统方法相比，能够在满足同样性能要求的情况下，消耗更少的硬件资源。
-
从当前的高速ADC器件指标和宽带雷达信号带宽参数看，宽带雷达数字接收机的抽取倍数
$M$ 并不高，通常采用两至三级抽取即可实现。在随后分析中，以两级抽取结构为例，相应的结论稍作扩展就适用于更多级数的抽取结构。由两级抽取器和可变分数时延滤波器组成的典型接收通道结构如图1a所示，图中的抗混叠滤波器均为通用线性相位FIR滤波器，其频率响应函数分别为${H_1}({\rm{j}}\omega )$ 和${H_2}({\rm{j}}\omega )$ ，抽取倍数$M = {M_1}{M_2}$ 。可变分数时延滤波器采用Farrow结构[15]实现，其分支滤波器采用通用线性相位FIR滤波器，总体频率响应表示为${H_d}({\rm{j}}\omega )$ 。显然，图1a给出的典型结构中，抽取处理是非线性的，且不同滤波器工作的数据率也不相同。因此，需要利用文献[13,15]给出的等价关系，将典型结构变换为图1b给出的等效结构，在此基础上分析整个通道的总体频率响应特性。
图 1 WB-DAR接收通道结构
基于Farrow结构实现的可变分数时延滤波器结构如图2所示，各分支滤波器系数为
$a(n,k)，n= $ $ 0, 1, 2,\cdots , N-1， k=0,1,2, \cdots ,L$ ，得到并行的分支滤波器和VFD滤波器的频率响应函数分别为：$$ {G_{{k}}}({\rm{j}}\omega ) = \sum\limits_{n = 0}^{N - 1} {a(n,k)} {{\rm{e}}^{ - {\rm{j}}n\omega }} $$ (1) $$ {H_{{d}}}({\rm{j}}\omega ) = \sum\limits_{n = 0}^{N - 1} {\sum\limits_{k = 0}^L {a(n,k){{\rm{e}}^{ - {\rm{j}}n\omega }}{d^k}} } = \sum\limits_{k = 0}^L {{G_{{k}}}} ({\rm{j}}\omega ){d^k} $$ (2) 同样，如果两级抗混叠滤波器系数为
${h_1}(n)$ 和${h_2}(n)$ ，则可得到两级滤波器的频率响应函数分别为：$$ {H_1}({\rm{j}}\omega ) = \sum\limits_{n = 0}^{{N_1} - 1} {{h_1}(n){{\rm{e}}^{ - {\rm{j}}n\omega }}} $$ (3) $$ {H_2}({\rm{j}}\omega ) = \sum\limits_{n = 0}^{{N_2} - 1} {{h_2}(n){{\rm{e}}^{ - {\rm{j}}n\omega }}} $$ (4) 进行等价变换后，整个WB-DAR系统接收通道的总体频率响应函数
$H({\rm{j}}\omega )$ 可以表示为：$$ H({\rm{j}}\omega ) = {H_1}({\rm{j}}\omega ){H_2}({\rm{j}}{M_1}\omega ){H_d}({\rm{j}}M\omega ) $$ (5) 针对
$H({\rm{j}}\omega )$ 的指标要求，设计目标函数和约束条件，求解出两级滤波器系数。
图 2 Farrow结构示意图
-
由于求解对象是级联滤波器的系数，这一类型的最优化求解问题不是凸优化问题，没有全局最优解，但通过合理设置初始化系数，可以求出近似最优解，从而使系统总体频率响应满足指标要求[16]。
在FIR滤波器求解过程中采用的最优化准则一般可以从最小最大化、最小二乘和约束最小二乘等不同准则中选择，最优化准则的选择并不影响求解近似最优解的过程。以下讨论的求解过程都基于最小最大化准则进行。
设总体频率响应函数
$H({\rm{j}}\omega )$ 满足：$$ \left\{ \begin{gathered} \left| {H({\rm{j}}\omega ) - {H_{{\rm{id}}}}({\rm{j}}\omega )} \right| \lt {\delta _{\rm{c}}}\quad \omega \in [0,\;{\omega _{\rm{c}}}] \hfill \\ \left| {H({\rm{j}}\omega ) - {H_{{\rm{id}}}}({\rm{j}}\omega )} \right| \lt {\delta _{\rm{s}}}\quad \omega \in [{\omega _{\rm{s}}},\;{\text π} ] \hfill \\ \end{gathered} \right. $$ (6) 式中，
$$ \left\{ \begin{gathered} {H_{{\rm{id}}}}({\rm{j}}\omega ) = {{\rm{e}}^{ - {\rm{j}}\omega (N/2 + d)}}\quad \omega \in [0,\;{\omega _{\rm{c}}}] \hfill \\ {H_{{\rm{id}}}}({\rm{j}}\omega ) = 0\;\quad \;\;\;\;\;\;\;\quad\quad \omega \in [{\omega _{\rm{s}}},\;{\text π}] \hfill \\ \end{gathered} \right. $$ (7) 且
${\delta _{\rm{c}}}$ 为通带内纹波；${\delta _{\rm{s}}}$ 为阻带内纹波。采取抽取时过渡带不混叠的方式来提高信噪比，将参数设置为：${\omega _{\rm{s}}} = {{\text π} \mathord{\left/ {\vphantom {{\text π} M}} \right. } M}$ 且${\omega _{\rm{c}}} = {{\text π} \mathord{\left/ {\vphantom {{\text π} {M - \Delta }}} \right. } {M - \Delta }}$ ，$\Delta \gt 0$ 。为了得到更好的滤波器性能，定义加权函数：
$$ W({\rm{j}}\omega ) = \left\{ \begin{gathered} 1\quad\quad\quad\quad \omega \in [0,\;{\omega _{\rm{c}}}] \hfill \\ \delta = {{{\delta _{\rm{c}}}} \mathord{\left/ {\vphantom {{{\delta _{\rm{c}}}} {{\delta _{\rm{s}}}}}} \right. } {{\delta _{\rm{s}}}}}\quad \omega \in [{\omega _{\rm{s}}},\;{\text π} ] \hfill \\ \end{gathered} \right. $$ (8) 加权后误差函数的形式如下：
$$ E({\rm{j}}\omega ) = W({\rm{j}}\omega )H({\rm{j}}\omega ) $$ (9) 则求最优解的过程即为寻找未知的滤波器系数
${h_1}(n)\; (n = 0,1,2,\; \cdots ,\;{{{N_1}} \mathord{\left/ {\vphantom {{{N_1}} 2}} \right. } 2})$ 、${h_2}(n)\;(n = 0,1,2,\; \cdots ,\;{{{N_2}} \mathord{\left/ {\vphantom {{{N_2}} 2}} \right. } 2} )$ 、$a(n,k)(n=0,1,2, \cdots, {N}_{m}/2，k=0,1,2, \cdots ,L )$ 以及误差$\delta $ 。实现最小化$\delta $ 的同时满足：$$ \left| {E({\rm{j}}\omega )} \right| \lt \delta $$ (10) 当最小化后的
$\delta $ 满足$\delta \leqslant {\delta _{\rm{c}}}$ 时，得到的总体滤波器频率响应就能够达到式(6)的要求。求解时，使用Matlab软件提供的最优化函数fminimax来完成多次迭代求解。如果在求解时能够给出合适的初始化滤波器系数，则可以降低对滤波器阶数的要求，使求解过程收敛速度加快。求解过程如下。
1)保持
${\delta _{\rm{c}}}$ 和${\delta _{\rm{s}}}$ 不变，设置${\omega _{{\rm{s}}1}} = {M_1}{\omega _{\rm{s}}}$ ，${\omega _{{\rm{c}}1}} = {\omega _{\rm{c}}}$ ，$d = 0$ ，将式(6)～式(9)中的$H({\rm{j}}\omega )$ 替换成${H_1}({\rm{j}}\omega )$ ，利用fminimax求解出幅频响应满足要求的${{{N_1}} \mathord{\left/ {\vphantom {{{N_1}} 2}} \right. } 2} + 1$ 个系数${h_1}(n)$ 。2)保持
${\delta _{\rm{c}}}$ 和${\delta _{\rm{s}}}$ 不变，设置${\omega _{{\rm{s}}2}} = {M_1}{\omega _{\rm{s}}}$ ，${\omega _{{\rm{c}}2}} = $ $ {M_1}{\omega _{\rm{c}}}$ ，$d = 0$ ，将式(6)～式(9)中的$H({\rm{j}}\omega )$ 替换成${H_2}({\rm{j}}\omega )$ ，利用fminimax求解出幅频响应满足要求的${{{N_2}} \mathord{\left/ {\vphantom {{{N_2}} 2}} \right. } 2} + 1$ 个系数${h_2}(n)$ 。3)保持
${\delta _{\rm{c}}}$ 和${\delta _{\rm{s}}}$ 不变，设置${\omega _{{\rm{s}}3}} = M{\omega _{\rm{s}}}$ ，${\omega _{{\rm{c}}3}} = M{\omega _{\rm{c}}}$ 和多个期望的${d_i}$ ，将式(6)～式(9)中的$H({\rm{j}}\omega )$ 替换成${H_{{d}}}({\rm{j}}\omega )$ ，利用fminimax求解出群时延误差满足要求的$({{{N_m}} \mathord{\left/ {\vphantom {{{N_m}} 2}} \right. } 2} + 1)(L + 1)$ 个系数，然后按照式(5)评估得到的总体频率响应函数特性，如果不满足要求，则修改滤波器阶数，重复步骤1)～步骤3)，直到满足要求。4)将得到的
${h_1}(n)$ 、${h_2}(n)$ 和$a(n,k)$ 作为${{{N_1}} \mathord{\left/ {\vphantom {{{N_1}} 2}} \right. } 2} + $ $ {{{N_2}} \mathord{\left/ {\vphantom {{{N_2}} 2}} \right. } 2} + ({{{N_m}} \mathord{\left/ {\vphantom {{{N_m}} 2}} \right. } 2} + 1)(L + 1) + 2$ 个滤波器系数的初始解，根据式(6)～式(9)，将总体频率响应函数给出的${\delta _{\rm{c}}}$ 、${\delta _{\rm{s}}}$ 、${\omega _{\rm{s}}}$ 和${\omega _{\rm{c}}}$ 代入，利用fminimax求解出满足要求的多级滤波器的最优系数。5)在得到第一次的最优解后，减少
${N_1}$ 、${N_2}$ 、${N_m}$ 和$L$ ，重复步骤1)～步骤4)，直到阶数减少后，求解出的滤波器系数对应的频率响应不再满足要求，则上一轮得到的滤波器系数就是阶数最少的近似最优解。为了减少总执行次数，上述过程可以按照通常使用的FIR滤波器阶数估算公式计算出
${N_1}$ 和${N_2}$ 的初始值，按照两分法降阶，可有效地降低计算次数。下面通过一个设计实例来说明本文设计方法的应用性。 -
本文设WB-DAR接收通道抽取倍数
$M = 6$ ，分两级抽取实现${M_1} = 2$ 和${M_2} = 3$ ，通带内纹波${\delta _{\rm{c}}} = $ $ 0.002\;3$ ，阻带纹波${\delta _{\rm{s}}} = 0.003\;0$ ，通带截止频率${\omega _{\rm{c}}} = $ $ {{0.5{\text π} } \mathord{\left/ {\vphantom {{0.5{\text π} } M}} \right. } M}$ ，阻带起始频率${\omega _{\rm{c}}} = {{0.5{\text π}} \mathord{\left/ {\vphantom {{0.5{\text π}} M}} \right. } M}$ ，分数时延间隔$\Delta d = 0.1$ ，时延误差小于0.001。先设置滤波器初始系数为全零，从步骤1)～步骤3)分别利用fminimax函数求解出两级抗混叠滤波器和VFD滤波器系数，评估接收通道总体频率响应函数是否满足要求，当满足要求时，对应的滤波器阶数在表1中给出。随后，以此时的滤波器系数为初始系数，再进行联合优化求解，如果得到的接收通道总体频率响应函数满足要求，则减少各级滤波器阶数，重复步骤1)～步骤4)，经过数次重复后，总体频率响应函数满足要求时的各级滤波器最少阶数也在表1中给出。
表 1 总体频率响应函数满足要求时两种方法所要求的各级滤波器最少阶数对比
方法 第一级抗混
叠滤波器阶
数(${N_1}$)第二级抗混
叠滤波器阶
数(${N_2}$)VFD滤波器
分支滤波器
阶数(${N_m}$)VFD滤波器
分支滤波器
个数($L + 1$)传统方法
(独立优化)12 36 13 4 本文方法
(联合优化)10 10 11 4 由表1中给出的滤波器阶数可知，采用本文提出的优化设计方法可以有效减少各级滤波器的阶数。新旧方法均采用横截型滤波器结构实现，乘法器个数=滤波器阶数+1。可见，本文方法能有效减少硬件乘法器的使用，从而降低整个WB-DAR系统接收通道的资源消耗。数字滤波器使用的乘法器和加法器个数由滤波器阶数决定，采用同样的级联型或其他结构实现时，本文方法亦可降低资源消耗。
采用本文方法得到的接收通道的总体幅频响应和群延时特性曲线如图3所示，图中的曲线对应
${d_i} = 0,0.1,0.2,0.3,0.4,0.5$ 的不同时延情况，仿真结果显示，设计满足指标要求。
图 3 接收通道频率响应特性曲线
-
WB-DAR系统抗干扰能力强，可灵活切换多个工作模式，还能够对目标进行分类和识别，为下一代的软件定义雷达提供了理想的实现平台。数字接收通道是WB-DAR系统的主要组成之一，因此，研究数字接收通道的优化设计方法对降低WB-DAR系统复杂度有重要意义。在WB-DAR系统的数字接收通道中，多级抽取和数字时延模块是主要的资源消耗部分，本文在通用FIR滤波器基础上，基于多级抽取结合可变分数时延滤波器的结构进行研究，提出了一种新的优化设计方法，该方法简单易行，与传统方法相比，在满足同样性能要求的情况下，消耗更少的硬件资源。
Research on Optimal Design of Receiving Channel in Wideband Digital Array Radar
-
摘要: 宽带数字阵列雷达(WB-DAR)具有模式多变、抗干扰能力强的优势，还兼具高距离分辨力，能够获取未知目标距离像，对目标进行分类和识别。大带宽信号的使用引入了孔径渡越时间的约束，因此，宽带数字阵列雷达需采用基于时延的宽带波束形成方法，以数字时延补偿方式取代传统窄带相控阵中的移相器来形成波束。该文在宽带数字阵列基础上，针对多级抽取结合可变分数时延滤波器的典型接收通道结构，提出了一种新的接收通道优化设计方法。与原方法相比，在满足同样性能要求的情况下，该方法能够降低抗混叠滤波器和分数时延滤波器的阶数，有效地节约了系统硬件资源。Abstract: Wide-band digital array radar (WB-DAR) has not only good performance in anti-interference and multifunction but also the advantage of high range resolution and the ability of target classification and recognition. On the other hand, the large bandwidth signal imposes the constraints of aperture transit time, therefore WB-DAR needs variable fractional delay (VFD) filter to obtain digital delay compensation instead of phase shifter in traditional narrowband phased array for accomplishing beamforming. Based on the typical structure of one receiving channel in WB-DAR, which consists of multiple stage anti-aliasing filters and a variable fractional delay (VFD) filter, a new optimal design method is proposed. Compared with the traditional method, the proposed method can apparently reduce the filter order of receiving channel and maintain the same performance.
-
表 1 总体频率响应函数满足要求时两种方法所要求的各级滤波器最少阶数对比
方法 第一级抗混
叠滤波器阶
数(${N_1}$ )第二级抗混
叠滤波器阶
数(${N_2}$ )VFD滤波器
分支滤波器
阶数(${N_m}$ )VFD滤波器
分支滤波器
个数($L + 1$ )传统方法
(独立优化)12 36 13 4 本文方法
(联合优化)10 10 11 4 
下载: 导出CSV
-
[1] TALISA S H, O'HAVER K W, THOMAS M, et al. Benefits of digital phased array radars[J]. Proceedings of the IEEE, 2016, 104(3): 530-543. doi: 10.1109/JPROC.2016.2515842 [2] LU Q, WANG X G. A new wideband digital array radar (WB-DAR) experiment system[C]//2014 IEEE Region 10 Symposium. Malaysia: IEEE, 2014: 440-445. [3] MORGAN L, ANDERSSON H. An efficient beamforming method using a combination of analog true time and digital delay[C]//Proceedings of the 2002 IEEE Radar Conference. Long Beach: IEEE, 2002: 260-265. [4] FISCHMAN M A, LE C. Digital beamforming developments for the Joint NASA-Air force space based radar[C]//2004 IEEE International Geoscience and Remote Sensing. Anchorage: IEEE, 2004: 687-690. [5] 何子述, 金林, 韩蕴洁, 等. 光控相控阵雷达发展动态和实现中的关键技术[J]. 电子学报, 2005, 33(12): 2191-2195. doi: 10.3321/j.issn:0372-2112.2005.12.022 HE Z S, JIN L, HAN Y J, et al. Development and implementation techniques of optically controlled phased array radar[J]. Acta Electronica Sinica, 2005, 33(12): 2191-2195. doi: 10.3321/j.issn:0372-2112.2005.12.022 [6] JUSWARDY B, ALAMEH K, LEE Y T. Integrated microphotonic beamformer for broadband adaptive nulling in smart antennas[C]//2006 IEEE Annual Wireless and Microwave Technology Conference. Clearwater Beach: IEEE, 2006: DOI: 10.1109/WAMICON.2006.351937. [7] TUR M. True time delay photonic beamforming: A review[C]//2009 IEEE International Conference on Microwaves, Communications, Antennas and Electronics Systems. Tel Aviv: IEEE, 2009: 1-2. [8] STEYSKAL H, COX C, ACKERMAN E. Digital and optical beamforming techniques[C]//2003 IEEE International Symposium on Phased Array Systems and Technology. Boston: IEEE, 2003: [9] QUAZI A H. An overview on the time delay estimate in active and passive Systems for target localization[J]. IEEE Transactions on Acoustic, Speech, Signal Processing, 1981, 29(3): 527-533. doi: 10.1109/TASSP.1981.1163618 [10] 王敏, 吴顺君, 杨淑媛. UWB脉冲信号的时域波束形成方法[J]. 电波科学学报, 2006, 21(2): 238-243. doi: 10.3969/j.issn.1005-0388.2006.02.017 WANG M, WU S J, YANG S Y. Time domain beamforming for UWB pulse signal[J]. Chinese Journal of Radio Science, 2006, 21(2): 238-243. doi: 10.3969/j.issn.1005-0388.2006.02.017 [11] 成超, 李会勇, 何子述. 基于子阵时延的数字阵列宽带波束形成[J]. 雷达科学与技术, 2008, 6(6): 459-462. doi: 10.3969/j.issn.1672-2337.2008.06.012 CHENG C, LI H Y, HE Z S. Using fractional delay in subarray to realize WDAR beamforming[J]. Radar Science and Technology, 2008, 6(6): 459-462. doi: 10.3969/j.issn.1672-2337.2008.06.012 [12] CHEUNG C, SHAH R, PARKER M. Time delay digital beamforming for wideband pulsed radar implementation[C]//2013 IEEE International Symposium on Phased Array Systems and Technology. Waltham, MA: IEEE, 2013: 448-455. [13] 陈祝明. 软件无线电技术基础[M]. 北京: 高等教育出版社, 2007. CHEN Z M. The basis of software radio[M]. Beijing: Higher Education Press, 2007. [14] MOLNAR G, VUCIC M. Closed-form design of CIC compensators based on maximally flat error criterion[J]. IEEE Transactions on Circuits and Systems II: Express Briefs, 2011, 58(12): 926-930. doi: 10.1109/TCSII.2011.2172522 [15] LAAKSO T I, VALIMAKI V, KARJALAINEN M, et al. Splitting the unit delay[J]. IEEE Signal Processing Magazine, 1996, 13(1): 30-60. doi: 10.1109/79.482137 [16] JOHANSSON H, GÖCKLER H. Two-stage-based polyphase structures for arbitrary-integer sampling rate conversion[J]. IEEE Transactions on Circuits and Systems II: Express Briefs, 2015, 62(5): 486-490. doi: 10.1109/TCSII.2015.2389291 -
点击查看大图
图(3) / 表(1)
计量
- 文章访问数: 5775
- HTML全文浏览量: 1700
- PDF下载量: 71
- 被引次数: 0
