②(华中科技大学电子与信息工程系 武汉 430074)
②(Key Laboratory of Opto-Electronics Information and Technical Science, Ministry of Education, Tianjin 300072, China)
③(School of Electronic Engineering, University of Electronic Science and Technology of China, Chengdu 611731, China)
太赫兹波(Terahertz waves)是指频率为0.1~10 THz (波长3 mm~30 μm)的电磁波[1, 2],该波段在基础研究、生物医学、安全监测等领域具有广泛的应用前景[3, 4, 5, 6, 7]。而近年来国内外的科研机构对太赫兹雷达领域越来越关注[8, 9, 10, 11]。太赫兹雷达研究的蒸蒸日上一方面是因为太赫兹波长比微波短,因此在相同积累角的情况下,太赫兹雷达将比微波雷达具有更高的横向分辨率[12, 13, 14];另一方面是因为太赫兹波能穿过大部分织物和非金属的包装,具有很强的透视效果;另外,太赫兹波光子能量较低,对人体无害,是一种非电离性的、安全的探测技术[15, 16]。目前,大部分太赫兹雷达成像技术都是基于单频或扫频方式完成的,频段大都集中在0.3 THz以下,难以达到标准的太赫兹波段(1 THz及以上),因而没能完全发挥太赫兹短波长所带来的横向分辨率优势;而又由于扫频的频谱宽度一般较窄(10 GHz),极大地限制了成像纵向分辨率的提高,而基于时域光谱技术的太赫兹雷达系统在这两方面显示出特有的优势。时域雷达系统一般采用宽带光导开关发射器和接收器,利用超短激光脉冲作为抽运源,将接收到的太赫兹电场信号转化为微弱电流信号放大,并通过测量电流信号来检测太赫兹脉冲。太赫兹时域光谱系统具有瞬态性(ps级时间分辨率)、宽带性(频谱宽约0.1~5 THz)、低能性(meV量级)和信噪比高(10000:1左右)等特点,其频率多在1 THz附近,频谱宽度大于1.5 THz,因此具有高频率、宽带宽的特征。在该方向,2010年,丹麦技术大学的Jepsen等人利用飞秒激光放大级为光源,以波面倾斜激发的太赫兹时域频谱系统组建了近单站式太赫兹雷达,并利用后向投影算法(BP)对F-16战斗机的1:150金属缩比模型进行了成像,验证了太赫兹时域雷达的成像能力[17, 18, 19]。
本文运用自行搭建的,以钛宝石飞秒激光振荡级为抽运源的太赫兹时域雷达紧缩场,对方形金属体、金属圆球、F-16战机模型、F-22战机模型和辽宁号缩比模型等多种结构和目标通过BP算法进行了成像研究。太赫兹成像的结果显示我们的太赫兹时域雷达系统具有成像能力,而且其理论分辨率极高。本文工作还对后像投影算法所产生的环形中心增强效应进行了分析讨论,并提出了成像背景圆环的产生原因与抑制方法,该研究为太赫兹波段的成像研究提供了新的思路。
2 太赫兹时域雷达系统图 1为搭建的雷达系统示意图,该系统以太赫兹时域光谱系统为基础,利用飞秒激光抽运的光电导天线产生、探测太赫兹波,其激励源为钛宝石飞秒脉冲激光器(脉宽为30 fs,中心波长为800 nm,重复频率为88 MHz)。天线辐射产生的宽带太赫兹波(0~1.3 THz),通过两面抛面镜(焦距10.16 cm)一面球面镜(焦距60.96 cm)扩束准直后形成宽平行光(光束全宽为120 mm)。光束向被测物体传播了1670 mm,散射光由距离被测物1310 mm远处的球面镜(焦距60.96 cm)与抛物面镜(焦距10.16 cm)组成的光学系统收集会聚到接收天线上,入射光束方向与反射光束方向之间的夹角为13°,系统的整体成像范围80 mm×80 mm×80 mm。如图 1所示,被测模型固定在旋转台(角度分辨率0.001°)上,通过旋转转台可测量被测模型不同角度方向的散射信号。
被测模型与旋转台之间由聚苯乙烯塑料泡沫支撑柱连接,在模型测量过程中,该支撑柱对于太赫兹波的散射极小,可忽略不计。图 2展示了雷达系统测量平面反射镜所得的太赫兹雷达信号,平面镜的放置姿态经过调整使得反射信号最大,图中可见其时域电场振幅信噪比可达3392:1。图 2附图展示了用雷达系统单独测量塑料泡沫载物台的时域信号,经实测发现,载物台的姿态和位置对所测得的结果没有影响,由附图可见载物台的散射信号在噪声水平以下(尺寸为12×20 cm),不会对后面测量带来不利影响。
为验证太赫兹雷达系统的探测能力,我们用不同直径的导体金属球的RCS对雷达系统进行校准。所谓RCS是一个用面积单位表征的参数,它定量表征了目标散射到雷达探测端的电磁波强度[20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30],其普遍定义为:
${\rm{RCS}} = \mathop {{\rm{lim}}}\limits_{R \to \infty } 4{{π}}{R^2}\frac{{{{\left| {{{E}_{\rm{s}}}} \right|}^2}}}{{{{\left| {{{E}_{\rm{i}}}} \right|}^2}}}$ | (1) |
其中Ei和Es分别是入射与散射电磁场振幅,而R是雷达与目标物体之间的距离。根据经典Mie理论[31],完美电导体球的散射截面为RCS=${π}r^2$而且 $\rm{RCS} \propto {\left| {{{E}_{\rm{p}}}} \right|^2}$,其中${{E}_{\rm{p}}}$为峰值电场,由此可见Ep~r。我们利用太赫兹雷达系统分别测量了直径0.9 cm,2 cm,3 cm,4 cm,5 cm,6 cm的金属球的峰值电场强度,在图 3中用误差棒表示。可以看出它们是成线性关系的,其比例系数经拟合可得(红线),此结果与Ep~r的关系相符合,证明了系统的有效性。
太赫兹雷达系统可以测量一定体积目标的散射信号,因此可以用于成像,但是时域太赫兹雷达所测得的是宽频时域信号,这与常用窄带雷达非常不同。时域信号的成像处理方法主要有两种:时域相关(Time Domain Correlation,TDC)法和后向投影法(Back-Projection,BP)[32, 33, 34]。时域相关法成像就是通过匹配滤波方法来实现成像,其运算量巨大,而且在大转角成像过程中,散射点到雷达的距离轨迹是弯曲的,而且各点轨迹的弯曲程度均有差异,若采用TDC算法成像,还需要进行大量的插值运算来保证成像精度,更进一步地增加了它的运算量,成像周期会很长,因而时域相关法在实际成像中几乎没有应用。基于时域信号的成像处理方法中最具代表性的就是后向投影算法[24],它的思想来源于计算机层析成像(Computerized Tomography,CT)技术,BP成像算法的基本思想是在回波中沿每个散射点到雷达的距离轨迹对其进行相干累加以获得该散射点的图像,因此对雷达系统的信号带宽和孔径并无太多限制,特别适合用于大转角和全角度的扫描成像。虽然BP算法同样存在计算量大的问题,但是BP算法成像中不存在近似,成像精度较高,适合于超宽带时域太赫兹雷达的成像处理,而且运动补偿相对易于实现,因此本实验采用BP算法对目标全角度回波信号进行处理从而获得目标2维图像。
同时,太赫兹时域雷达系统是目标转动,因此属于逆合成孔径雷达(ISAR)。ISAR是利用目标和雷达之间的相对运动实现对目标的高分辨成像的雷达系统。ISAR的距离向分辨依靠的是发射宽带信号,发射天线发射宽带信号并对回波信号作匹配滤波,通过测量滤波输出的延时就能获得散射点的径向距离信息,距离分辨率$\rho_{_{\rm{\Large R}}}$与发射信号的带宽相关[13]
${\rho_{_{\rm{\Large R}}}} = \frac{{\rm{c}}}{{2{B}}}$ | (2) |
其中c为光速,B为发射信号带宽。雷达系统的带宽为1.2 THz,中心频率为0.6 THz,根据式(5)计算得到的距离向理论分辨率为0.125 mm。
ISAR的方位向分辨则利用雷达和目标之间相对运动所产生的合成孔径来实现。方位向分辨的关键在于多普勒频率的差异,为了形成有效的多普勒差异,目标相对于雷达必须存在姿态视角变化,引起目标姿态变化的运动形式可以是多种多样的,其中雷达固定不动,目标转动的转台模型是一种简单的成像方式[35]。
转台模型如图 4所示,目标围绕聚焦中心匀速转动,转动角频率为$\omega$,雷达到参考点的距离为ra,目标上的某散射点在起始时刻(t=0)相对于聚焦中心的坐标为(r1,$\theta _0$),则在以后的时间里该点到雷达的距离为:
$r\left( t \right) = {\left[{r_0^2 + r_a^2 + 2{r_0}{r_a}\sin \left( {{\theta _0} + \omega t} \right)} \right]^{1/2}}$ | (3) |
实际工作中,雷达与目标之间的距离往往远大于目标尺寸(${r_0} \lt \lt {r_a}$),则式(6)可近似为
$r\left( t \right) = {r_a} + {x_0}\sin \left( {\omega t} \right) + {y_0}\cos \left( {\omega t} \right)$ | (4) |
其中${x_0} = {r_0}\cos {\theta _0}$,${y_0} = {r_0}\sin {\theta _0}$。通过距离关系,对目标区域进行重构,从而完成目标物体的重建。
大转角方位向分辨率公式为[37]
${\delta _{}} = \frac{\lambda }{{4\sin \left( {\theta/2} \right)}}$ | (5) |
我们采用的是360°扫描,步进角度为1°。系统采用宽光束平行光照射物体获得物体的回波信号,得到的是物体的外轮廓信息,成像的结果为压缩的截面图,且与对称性无关,故而,从0到${π}$的旋转和${π}$到2${π}$的转动得到的信号所携带的信息是相同的,因此式(5)中$\theta$取${π}$,通过计算可以得到理论方位向分辨率为0.125 mm。
同时在传统BP成像算法的基础上,对得到的信号进行包络的提取,获得回波信号的能量信息,并以此作为成像数据,从成像结果上来看比直接采用回波信号的强度成像有较大改观。不仅如此,考虑到物体在个别角度下的散射信号比其他角度散射回波信号会有数量级的差异,因此还对成像数据取对数,从而减小差异,便于提取目标的形貌信息。
图 5是我们通过测量磁座(一端带有旋钮,有on和off两种状态决定是否吸附到平台上)来验证成像方法的有效性所成的结果,采样角度为360°,角度采样间隔为1°,目标方柱的长为6 cm,宽为4 cm,从成像的结果来看尺寸几乎相当,但是也有出入,方柱正对一侧的表面有一些凹进去,另一侧则有一个旋钮并且表面有一层稍有脱落的塑料膜,引起了散射上下两侧中间部分比较弱,同时另外两个侧面相对平整,但是其相对粗糙的表面使得漫反射比较严重,导致目标物体左右两个侧面的散射信号较弱。其中上边沿在成像结果中比较明显,这是由于其在实物中为抛光的平面。同时也注意到,图像的中心会有一个亮斑,这是由BP成像算法中通过角度旋转相干叠加的方法来重构目标区域的物体所引起的,中心累加次数较多,故会明显地被加强,可以通过对采集信号滤波和对噪声水平进行补偿的方法来抑制中心加强的现象。
图 6中,我们为了测试系统的实际分辨率,对相距6 mm的直径为3 cm和4 cm的钢球进行了360°全角度测量,角度采样间隔为1°,其中图 6(a)是成像结果,从图中可以看出是两个大小不同的圆球,尺寸上的小差异是因为重建图像时雷达相位中心的选取有偏差,导致重构出来的小球的尺寸被动地展宽了。同时也可以注意到,两圆球靠近的地方信号相对较弱,分析主要由以下的原因造成:首先两个小球中间区域的信号在反射时会在之间的球面多次反射,引起光程上的差别,无法被探到;其次,系统并非是单站系统,双站系统之间有夹角,两个钢球连接处相当于二面角反射,信号出射方向与接收之间偏差会引起信号无法被探到。这些原因都不同程度地引起了中心附近轮廓不是很明显。为了测定系统的分辨能力,我们通过后期重构的方法来鉴定系统是否能将相距6 mm的钢球区分开,如图 6(b)所示,灰色的圆是我们根据成像结果拟合得到的,可以明显看到,系统的分辨能力应高于6 mm。而背景中的圆环是由场景中某些非转动物体的反射信号经过后向投影算法旋转叠加而成,由于这些物体是不动的,所以反射信号大小在一个相对稳定的范围,可通过观察距离向历史选取合适的值来抑制圆环,也可以在下一步的实验中贴附吸波材料减少非被测目标的反射信号强度,以及选取合适窗函数对时域回波进行优化滤波等方式,最大程度减弱成像中环状背景。
除此之外,我们还对飞机、航母等缩比模型进行了成像研究,用以更好地理解物体的散射特性。我们对缩放比为1:2000辽宁号航空母舰模型进行了角采样间隔为1°的全角度测量,采用时域后向投影算法对航母模型进行成像,来观察其散射特性。图 7(a)展示了航母模型的成像结果,并与图 7(b)中的航母模型进行比较。由图 7(a)可见,成像结果与实际模型符合得很好。不仅如此,航母模型的舰桥、甲板等细节都能够分辨[37]。
此外,航母外轮廓上散射较强的特征点,雷达系统能够接收到较强的信号,这些强散射点在成像结果中也有所体现。图中红色的部分表示强信号点,也就是强散射点,如舰岛的结构相比船体更为复杂精细,散射点更多,故而呈现红色。
图 8则是对1:200的F-22战斗机的成像结果与模型,通过比较可以发现图中可以观测到良好的对称性,这与F-22战斗机的对称外形相符合。此外,由文献[36]可知,隐形飞机的特定设计具有整体扁平,表面由一些平面组成的特点,有效降低了其散射效应,使雷达回波主要集中在有限的几个方向(飞机背部,飞机腹部),类似镜面反射。从文献[37]中可得知战斗机进气道形成了类似平面形状物体的散射特性,所以信号较大,从图中也能看出,尾翼与背部也形成了类似的镜面反射,散射信号较强,所以能从成像结果中观测到,而机翼则相对很弱,这与之前文献[36]中的RCS结果符合得很好。
综上所述,本文实现了宽带太赫兹波的时域散射信号测量,并基于此研究了缩比模型的雷达后向投影成像。我们的太赫兹时域雷达成像范围为80 mm×80 mm×80 mm,具有分辨6 mm以上空间结构的分辨能力。基于后向投影算法原理,我们对多种目标模型的太赫兹雷达回波数据进行了图像重构,获得了较好的成像结果。时域太赫兹雷达1.3 THz的带宽,分米量级的模型尺寸和极高的成像分辨率必将会在雷达设计、目标识别、电磁隐身和外形设计等领域发挥独特的作用。
致谢 感谢合肥工业大学尹治平教授对本论文提出的宝贵意见。[1] | Siegel P H. Terahertz technology[J]. IEEE Transactions on Microwave Theory and Techniques, 2002, 50(3): 910–928.(1) |
[2] | Zhang Tong-yi and Cao Jun-cheng. Study of the surface and far fields of terahertz radiation generated by large-aperture photoconductive antennas[J]. Chinese Physics B, 2004, 13(10): 1742–1746.(1) |
[3] | 许景周, 张希成. 太赫兹科学技术和应用[M]. 北京: 北京大学出版社, 2007: 206. Xu Jing-zhou and Zhang Xi-cheng. Applications of Terahertz Science and Technology[M]. Beijing: Peking University Press, 2007: 206.(1) |
[4] | Lee Yun-shik, 崔万照, 等译. 太赫兹科学与技术原理[M]. 北京: 国防工业出版社, 2012: 192–200. Lee Yun-shik. translated by Cui W Z, et al.. Principles of Terahertz Science and Technology[M]. Beijing: National Defence Industry Press, 2012: 192–200.(1) |
[5] | 王瑞君, 王宏强, 庄钊文, 等. 太赫兹雷达技术研究进展[J]. 激光与光电子学进展, 2013, 50(4): 4–20. Wang Rui-jun, Wang Hong-qiang, and Zhuang Zhao-wen, et al.. Research progress of terahertz radar technology[J]. Laser & Optoelectronics Progress, 2013, 50(4): 4–20.(1) |
[6] | Xu Zhuo, Gu Chao, Pei Zhi-bin, et al.. Multiband terahertz metamaterial absorber[J]. Chinese Physics B, 2011, 20(1): 017801.(1) |
[7] | Tian Lu, Zhou Qing-li, Jin Bin, et al.. Consistency-dependent optical properties of lubricating grease studied by terahertz spectroscopy[J]. Science in China, 2009, 39(11): 1589–1593.(1) |
[8] | Goyette T M, Dickinson J C, Waldman J, et al.. A 1.56 THz compact radar range for W-band imagery of scale-model tactical targets[C]. Algorithms for Synthetic Aperture Radar Imagery VII, Orlando, Florida, USA, 2000: 615–622.(1) |
[9] | DeMartinis G B, Coulombe M J, Horgan T M, et al.. A 240 GHz polarimetric compact range for scale model RCS measurements[C]. Antenna Measurements Techniques Association (AMTA), Atlanta, Georgia, USA, 2010: 3–8.(1) |
[10] | Dengler R J, Cooper K B, Chattopadhyay G, et al.. 600 GHz imaging radar with 2 cm range resolution[C]. IEEE/MTT-S International Microwave Symposium, Digest, Honolulu, HL, 2007: 1371–1374.(1) |
[11] | Chattopadhyay G, Cooper K B, Dengler R, et al.. A 600 GHz imaging radar for contraband detection[C]. International Symposium on Space Terahertz Technology 19th, Groningen, Netherlands, 2008: 300–303.(1) |
[12] | Cooper K B, Dengler R J, Llombart N, et al.. Penetrating 3-D imaging at 4- and 25-m range using a submillimeter-wave radar[J]. IEEE Transactions on Microwave Theory and Techniques, 2008, 56(12): 2771–2778.(1) |
[13] | 蔡英武, 杨陈, 曾耿华, 等. 太赫兹极高分辨力雷达成像试验研究[J]. 强激光与粒子束, 2012, 24(1): 7–9. Cai Ying-wu, Yang Chen, Zeng Geng-hua, et al.. Experimental research on high resolution terahertz radar imaging[J]. High Power Laser and Particle Beams, 2012, 24(1): 7–9.(2) |
[14] | 成彬彬, 江舸, 陈鹏, 等. 0.67 THz高分辨力成像雷达[J]. 太赫兹科学与电子信息学报, 2013, 11(1): 7–11. Cheng Bin-bin, Jiang Ge, Chen Peng, et al.. 0.67 THz high resolution imaging radar[J]. Journal of Terahertz Science and Electronic Information Technology, 2013, 11(1): 7–11.(1) |
[15] | Mickan S P and Zhang X C. T-Ray sensing and imaging[J]. International Journal of High Speed Electronics and Systems, 2003, 13(2): 601–676.(1) |
[16] | 冯伟, 张戎, 曹俊诚. 太赫兹雷达技术研究进展[J]. 物理, 2013, 42(12): 1–5. Feng Wei, Zhang Rong, and Cao Jun-cheng. Progress in Terahertz radar technology[J]. Physics, 2012, 42(12): 1–5.(1) |
[17] | Iwaszczuk K, Heiselberg H, and Jepsen P U. Terahertz radar cross section measurements[C]. 2010 35th International Conference on Infrared Millimeter and Terahertz Waves (IRMMW-THz 2010), Rome, Italy, 2010: 26399–26408.(1) |
[18] | Grischkowsky D, Keiding S, Exter M, et al.. Far-infrared time-domain spectroscopy with terahertz beams of dielectrics and semiconductors[J]. Journal of the Optical Society of America B, 1990, 7(10): 2006–2015.(1) |
[19] | McGowan R W, Cheville R A, and Grischkowsky D. Experimental study of the surface waves on a dielectric cylinder via terahertz impulse radar ranging[J]. IEEE Transactions on Microwave Theory and Techniques, 2000, 48(3): 417–422.(1) |
[20] | 江舸, 成彬彬, 张健, 等. 基于0.14 THz成像雷达的RCS测量[J]. 太赫兹科学与电子信息学报, 2014, 11(1): 19–23. Jiang Ge, Cheng Bin-bin, Zhang Jian, et al.. 0.14 THz radar imaging based radar cross section measurement[J]. Information and Electrionic Engineering, 2014, 11(1): 19–23.(1) |
[21] | 蒋彦雯, 邓彬, 王宏强, 等. 基于时域光谱系统的太赫兹圆柱RCS测量[J]. 红外与激光工程, 2014, 43(7): 2223–2227. Jiang Yan-wen, Deng Bin, Wang Hong-qiang, et al.. RCS measurement of cylinders in terahertz band based on the time-domain spectroscopy system[J]. Infrared and Laser Engineering, 2014, 43(7): 2223–2227.(1) |
[22] | 武亚君, 黄欣, 徐秀丽, 等. 太赫兹目标RCS缩比测量技术[J]. 强激光与粒子束, 2013, 25(6): 1541–1544. Wu Ya-jun, Huan Xin, Xu Xiu-li, et al.. Radar cross section measurement technique of scale-model targets at terahertz[J]. High Power Laser and Particle Beams, 2013, 25(6): 1541–1544.(1) |
[23] | Ausherman D A, Kozma A, Walker J L, et al.. Developments in radar imaging[J]. IEEE Transactions on Aerospace and Electronic Systems, 1984, 20(4): 363–400.(2) |
[24] | 王乃卫, 粟毅, 陆仲良, 等. 冲激雷达旋转目标成像方法[J]. 系统工程与电子技术, 1999, 21(4): 15–19. Wang Nai-wei, Su Yi, Lu Zhong-liang, et al.. Impulse radar imaging of rotating objects[J]. Systems Engineering and Electronics, 1999, 21(4): 15–19.(1) |
[25] | Liang M Y, Zhang C L, Zhao R, et al.. Experimental 0.22 THz stepped frequency radar system for ISAR imaging[J]. Journal of Infrared, Millimeter, and Terahertz Waves, 2014, 35(9): 780–790.(1) |
[26] | Jansen C, Krumbholz N, Geisb R, et al.. Alignment and illumination issues in scaled THz RCS measurements[C]. 34th International Conference on Infrared, Millimeter, and Terahertz Waves, Korea, Busan, 2009: 1–2.(1) |
[27] | Knott E F. Radar Cross Section Measurements[M]. New York: SciTech Publishing Inc., 1993: 16–17.(1) |
[28] | Danylov A A, Goyette T M, Waldman J, et al.. Terahertz Inverse Synthetic Aperture Radar (ISAR) imaging with a quantum cascade laser transmitter[J]. Optics Express, 2010, 18(15): 16264–16272.(1) |
[29] | Iwaszczuk K, Heiselberg H, and Jepsen P U. Terahertz radar crosssection measurements[J]. Optics Express, 2010, 18(25): 26399–26408.(1) |
[30] | Gente R, Jansen C, Geise R, et al.. Scaled bistatic radar cross section measurements of aircraft with a fiber-coupled THz time-domain spectrometer[J]. IEEE Transactions on Terahertz Science and Technology, 2012, 2(4): 424–431.(1) |
[31] | Born M and Wolf E. Principles of Optics[M]. Cambridge: Cambridge University Press, 1999: 633–644.(1) |
[32] | 黄培康, 殷红成, 许小剑. 雷达目标特性[M]. 北京: 电子工业出版社, 2005: 71–72. Huang Pei-kang, Yin Hong-cheng, and Xu Xiao-jian. Radar Target Characteristics[M]. Beijing: Publishing House of Electronics Industry, 2005: 71–72.(1) |
[33] | 保铮, 邢孟道, 王彤. 雷达成像技术[M]. 北京: 电子工业出版社, 2005: 6–18. Bao Zheng, Xing Meng-dao, and Wang Tong. Radar Imaging Technology[M]. Beijing: Publishing House of Electronics Industry, 2005: 6–18.(1) |
[34] | 张麟兮, 胡楚峰. 雷达目标散射特性测试与成像诊断[M]. 北京: 中国宇航出版社, 2009: 100–120. Zhang Lin-xi and Hu Chu-feng. Measurement and Imaging of Radar Target Scattering Characteristics[M]. Beijing: China Astronautic Publishing House, 2009: 100–120.(1) |
[35] | Walker J L. Range-Doppler imaging of rotating objects[J]. IEEE Transactions on Aerospace and Electronic Systems, 1980, 16(1): 23–52.(3) |
[36] | Liang Da-chuan, Wei Ming-gui, Gu Jian-qiang, et al.. Broad-band time domain terahertz radar cross-section research in scale models[J]. Acta Physica Sinica, 2014, 63(21): 214102.(1) |
[37] | Mensa D L. High Resolution Radar Cross-Section Imaging[M]. Boston: Artech House Inc., 1981: 85–115.(3) |