B SAR系统的主要优势在于能够获取较高分辨率(5～10 m)、超宽幅(350～400 km)的SAR图像数据，这意味着能够以短重访周期对地球系统的大范围动态变化过程进行较高精度的观测，因此在海洋遥感、环境监测、灾害评估等领域有着重要的应用价值[11,12]。传统的MEB SAR由于受到发射脉冲的干扰，在生成的图像中会存在固定不变的盲区。为了实现对大测绘带场景的连续观测，文献[13,14]提出了改进的成像模式staggered SAR，通过周期性连续变化系统脉冲重复间隔(PRI)来改变测绘带盲区的位置，使得整个测绘带内目标回波均能被接收到，从而消除测绘带内的盲区。此外，方位向多通道技术、多维波形编码技术和MIMO技术的引入还能进一步提升这类系统的成像性能[15–19]。尽管这些星载MEB SAR系统的具体处理方案不同，但是由于对大测绘带场景进行成像，这些系统都需要利用俯仰向的数字波束形成技术对来自不同子测绘带的回波进行分离处理。因此，俯仰向的DBF处理是星载MEB SAR系统实现较高分辨率、超宽幅成像的关键。通常，基于平面相控阵天线的星载MEB SAR系统在发射端产生俯仰向宽波束照射地面场景，在接收端利用俯仰向多通道接收地面回波。基于地球理想球体模型，不同子测绘带目标的回波时间和波达角(AOA)的对应关系可根据确定的成像几何关系计算得到。因此，这类系统可利用零点指向技术[20]生成合适的DBF接收波束，将波束方向图的中心对准需要提取的子测绘带目标回波，同时将波束方向图的零陷位置对准其他子测绘带目标的干扰回波，从而确保不同子测绘带回波能够有效分离。但是，由于存在卫星姿态误差等因素的影响，星载MEBSAR系统的DBF接收波束指向会出现偏差，使得波束方向图的零陷位置不能准确指向干扰回波的方向，这将导致在对具有强散射体的区域(如存在船舶的海面、港口等区域)进行超宽幅成像时会出现明显的鬼影目标，从而严重影响SAR图像的解译。针对这一问题，有一种解决思路是利用自适应数字波束形成(ADBF)进行处理[21,22]。该处理方法采用空间谱估计算法对每个距离门单元的回波信号进行波达角估计，然后利用波达角估计值计算DBF加权矢量来调整接收波束指向。但是这种方法通常需要很多个快拍数据来构造样本协方差矩阵并进行求逆运算，运算量较大，难以满足实时性要求，而且易受到多径效应和相干信号源的干扰，因此能否有效应用于星载SAR系统中还有待进一步研究。考虑矩阵束方法[23]能够在复杂环境中基于单快拍数据准确而且稳健地估计目标的波达角，本文提出了一种基于矩阵束方法的星载MEB SAR俯仰向DBF处理方法。本文详细分析了星载MEB SAR的成像几何模型和俯仰向多通道回波信号模型，并具体推导了主要的处理流程，点目标和分布目标的仿真实验将验证本文方法能够有效抑制鬼影目标的干扰。2    星载MEB SAR的成像几何模型和信号模型s(t) = rect(t=Tp) exp(j2 f0t) exp(j krt2)Tpf0krrect(t)[¡0:5; 0:5]N图1为星载MEB SAR的成像几何模型。星载MEB SAR系统利用高度较小的天线发射脉冲信号照射超宽幅测绘带，发射信号模型的表达式为：。 其 中为脉冲时宽， 为雷达载频， 为调频率，表示范围在 的归一化矩形窗。如图2所示，这类系统的接收天线在俯仰向由 个俯仰向子孔径组成，每个俯仰向子孔径的方向图都覆盖整个成像场景。由于星载MEB SAR的测绘带宽很大，M个子测绘带回波将在同一个雷达接收窗内相互重叠。考虑成像场景中任意位置目标的回波到接收天线各俯仰向子孔径的距离差远小于1个距离门，经  图 1 星载MEB SAR成像几何模型  图 2 星载MEB SAR俯仰向DBF处理示意图2660 电  子  与  信  息  学  报 第 40 卷)过去载频处理后，俯仰向第n个子孔径接收到的目标回波信号 可表示为rn(t) =MXm=1¾m Am rect µt¡t1Tp¶ exp (¡j2 f0tm)¢exp µj2 f0hnsin (µm¡µc)c¶¢exp³j kr(t¡t1)2´+ en(t) (1)¾mAm其中， 表示第m个子测绘带目标的后向散射系数， 表示第m个子测绘带目标受到的双程天线t1tm= t1+ m TTN hn= µn¡N + 12¶ hrhrµmµc c en (t)方向图幅度加权， 表示中间俯仰向子孔径接收到最近端子测绘带目标回波的双程时延。表示中间俯仰向子孔径接收到第m个子测绘带目标回波的双程时延， 表示系统脉冲重复间隔。假设为奇数， 表示第n个俯仰向子孔径相对于中间俯仰向子孔径的距离， 为俯仰向子孔径高度。 表示第m个子测绘带目标对应的下视角， 表示天线法线的下视角， 为光速。表示第n个俯仰向子孔径接收回波受到的噪声干扰。由式(1)可得到回波信号的导向矢量 表达式为=26666666666664exp µj2 f0h1sin (µ1¡µc)c¶¢ ¢ ¢exp µj2 f0h1sin (µM¡µc)c¶:::¢ ¢ ¢:::exp µj2 f0hnsin (µ1¡µc)c¶¢ ¢ ¢exp µj2 f0hnsin (µM¡µc)c¶:::::::::exp µj2 f0hNsin (µ1¡µc)c¶¢ ¢ ¢exp µj2 f0hNsin (µM¡µc)c¶37777777777775(2)m= [!m;1 !m;2¢¢¢!m;N]Tm尽管不同子测绘带目标的回波信号在接收窗内相互重叠，但是这些回波信号的波达角都不相同，因此可利用零点指向技术对不同子测绘带目标回波进行分离提取，例如采用线性约束最小方差(LCMV)的波束形成算法[20]。假定分离提取第m个子测绘带目标回波的加权矢量为 ，采用LCMV算法可得到 的表达式为Hm=Hm¡H¢¡1H¡1(3)mM£M¾2s其中， 表示 单位矩阵的第m列， 表示噪声的协方差矩阵。假设这一噪声为高斯白噪声且功率为 ，则 可表示为= ¾2s N£N (4)Hm=Hm¡H¢¡1H¿ µ将式(4)代入式(3)中，得到： 。基于地球理想球体模型，地球表面目标的回波时刻与波达角 之间的对应关系可表示为µ(¿) = arccosÃ4(Re+ H0)2¡4R2e+ (c¿)24c¿(Re+ H0)!(5)ReH0其中， 表示地球半径， 表示卫星平台轨道高度。在理想情况下，如果星载MEB SAR系统的天线法线指向准确，那么DBF加权矢量可根据式(3)计算得到，并且各个子测绘带回波可有效分离提取。但是由于存在卫星姿态误差等因素的影响，星载SAR系统的天线法线指向会出现偏差。如图1所示，如果仍然按照理想情况计算DBF加权矢量，并形成DBF波束对成像场景进行扫描接收以及分离提取各子测绘带回波，将导致接收波束中心不能准确对准需要提取的子测绘带回波，接收波束方向图的零陷指向不能对准来自其他子测绘带的干扰回波。当对具有强散射体的区域进行超宽幅成像时，比如存在船舶的海面、港口等区域，由于不同子测绘带的强散射体回波信号没有被有效抑制，SAR图像中将出现鬼影目标的干扰，严重影响图像解译。针对这一问题，本文提出了一种基于矩阵束方法的星载MEB SAR俯仰向DBF处理方法，下一节将对该处理方法进行详细阐述和推导。3    基于矩阵束方法的俯仰向DBF处理方法在理想情况下，星载MEB SAR系统根据预先计算好的DBF加权矢量对俯仰向多通道回波信号进行DBF处理。预设的俯仰向DBF加权矢量的主要计算步骤是：根据预先计划照射的成像场景计算每个距离门单元的回波双程时延，利用式(5)计算相应的目标下视角，然后结合已知的预设天线法线指向角计算目标导向矢量，最后根据式(3)和导向矢量计算每个距离门单元的DBF加权矢量。但是，当存在卫星姿态误差等因素的影响时，将导致星载SAR系统的天线法线指向出现偏差，天线方向图出现偏离，这将进一步导致成像场景目标的导向矢量发生改变，从而影响星载MEB SAR的俯仰向DBF处理。µcµ0c如图1所示，在这种情况下，星载MEB SAR系统的预设天线法线指向角 将偏离实际的天线法线指向角 ，这导致了在计算预设DBF加权矢量时第 1 1 期 叶  恺等：基于矩阵束方法的星载MEB SAR俯仰向DBF处理方法 2661 与成像场景的导向矢量 出现偏差，其中 的表达式为=26666666666664exp µj2 f0h1sin (µ1¡µ0c)c¶¢ ¢ ¢exp µj2 f0h1sin (µM¡µ0c)c¶:::¢ ¢ ¢:::exp µj2 f0hnsin (µ1¡µ0c)c¶¢ ¢ ¢exp µj2 f0hnsin (µM¡µ0c)c¶:::::::::exp µj2 f0hNsin (µ1¡µ0c)c¶¢ ¢ ¢exp µj2 f0hNsin (µM¡µ0c)c¶37777777777775(6)此时，如果仍然采用预设的DBF加权矢量对回波信号进行处理，将导致DBF接收波束中心和零陷指向都出现偏差。解决这一问题的关键在于根据实际的天线法线指向角和成像场景的目标下视角重新计算导向矢量，从而确保每个距离门单元的DBF加权矢量计算正确。由于每个距离门单元的距离向时间都是已知的，那么每个距离门单元对应的下视角可以根据式(5)计算得到，这意味着成像场景中目标下视角都是已知的。因此，如果能准确估计成像场景中的任一目标波达角，那么就能结合该目标的下视角计算得到实际的天线法线指向角，从而调整成像场景的导向矢量。根据上述解决思路，同时考虑具有高信噪比的强散射体回波有利于目标波达角的正确估计，下面将详细阐述基于矩阵束方法的俯仰向DBF处理方法。为了便于寻找强散射体的位置以及提升回波信号的信噪比，首先对俯仰向中间子孔径接收的回波信号进行匹配滤波处理，然后根据系统预设的阈值对回波信号脉冲压缩结果进行过滤，预设的阈值主要通过对SAR系统的地面测量和回波仿真计算得到。当脉冲压缩后的回波信号峰值都小于预设阈值时，则直接应用预设DBF加权矢量对所有俯仰向子孔径回波进行处理。当存在超过预设阈值的回波信号脉冲压缩结果时，进一步寻找最大幅度值所在的位置，即对应强散射体目标所在位置。提取最大幅度值位置处的俯仰向多通道回波信号的脉冲压缩结果，得到¹ (t) =h¹r1(t) ¹r2(t)¢¢¢¹rN(t)iT(7)¹rn(t) =MXm=1pjkrjTp¾m Am exp (¡j2 f0tm)¢sinc (krTp(t¡t1))¢exp µj2 f0hnsin (µm¡µc)c¶ + ¹en(t) (8)¹en其中，(t)表示匹配滤波处理后的噪声。接下来利用矩阵束方法对所提取数据进行处理并估计强散射体目标的波达角。矩阵束方法是一种用于参数估(N¡L)£(L + 1)计的直接数据域算法，能够利用单个快拍样本数据准确估计相干信号或非相干信号的波达角，这种算法不仅计算效率高，而且在存在噪声干扰的环境中具有稳健的估计性能。本文采用的矩阵束方法的主要处理步骤跟文献[24]中的描述大体相同，不同点在于考虑了星载MEB SAR系统的实际应用环境来设置数据矩阵 的奇异值筛选条件。根据式(7)和文献[24]，定义 的数据矩阵 ，其表达式为=26664¹r1(t) ¹r2(t)¢ ¢ ¢¹rL+1(t)¹r2(t) ¹r3(t)¢ ¢ ¢¹rL+2(t)::::::::::::¹rN¡L(t) ¹rN¡L+1(t)¢ ¢ ¢¹rN(t)37775(9)L·N3;N2¸为了有效减轻噪声的影响，参数 的取值需要在范围内选取[25]。对数据矩阵 进行奇异值分解，得到=H(10)= [ 1 2 ¢ ¢ ¢ N¡L ](N¡L)£(N¡L)H= [ 1 2¢¢¢L+1](L+1)£(L+1)H(N¡L)£(L + 1)Ns~µ其中， ，由矩 阵 的 特 征 向 量 组 成 。，由矩阵 的特征向量组成。为 的对角矩阵，对角线上的元素对应数据矩阵 的奇异值。接下来对矩阵 中的奇异值进行筛选，筛选的限制条件包含以下两点：由于星载MEB SAR系统的成像场景由M个子测绘带组成，选取的奇异值数目不能超过M；考虑强散射体回波信号具有较高的信噪比，选取的奇异值要求大于最大奇异值的一半。假定选取了 个奇异值，强散射体的波达角估计值 可通过式(11)计算得到：~µ = arcsin µarctan (imag (®) =real (®))2 f0hr¢c¶ (11)® = eig µ³¡¹1¢H´y¡¹2¢H¶ (12)¹= [ 1 2 ¢¢¢ Ns](L+1)£Ns¹1¹ ¹2¹ real (¢) imag (¢)其中， ，矩阵为删去了最后一行元素的矩阵 ，矩阵 为删去了第1行元素的矩阵 。 和 分别表示2662 电  子  与  信  息  学  报 第 40 卷基于矩阵束方法的星载MEB SAR俯仰向DBF处理方法叶   恺*①②      禹卫东①②      王   伟①①(中国科学院电子学研究所   北京   100190)②(中国科学院大学   北京   100049)摘   要：俯仰向数字波束形成(DBF)处理是距离向多波束体制(MEB)星载合成孔径雷达(SAR)系统实现较高分辨率、超宽幅成像的关键。但是，由于存在卫星姿态误差等因素的影响，星载MEB SAR系统的DBF接收波束指向会出现偏差，这将导致在对具有强散射体的区域(如存在船舶的海面、港口等区域)进行成像时会出现鬼影目标。针对这一问题，该文提出一种基于矩阵束方法的俯仰向DBF处理方法。首先对俯仰向通道回波数据进行匹配滤波处理，并根据预设阈值寻找强散射体的峰值位置；然后利用矩阵束方法准确估计强散射体的波达角；最后利用这些信息调整俯仰向DBF加权矢量，确保DBF接收波束指向正确，从而消除鬼影目标的干扰。仿真实验验证了该方法的有效性。关键词：合成孔径雷达；距离向多波束；数字波束形成；矩阵束方法中图分类号：TN959.74 文献标识码：A 文章编号：1009-5896(2018)11-2659-08DOI: 10.11999/JEIT180076Matrix Pencil Method Based Processing Approach for SpaceborneMEB SAR with Digital Beamforming in ElevationYE Kai①②      YU Weidong①②      WANG Wei①①(Institute of Electronics, Chinese Academy of Sciences, Beijing 100190, China)②(University of Chinese Academy of Sciences, Beijing 100049, China)Abstract: Digital Beam Forming (DBF) in elevation plays a crucial role for spaceborne Multiple Elevation Beam(MEB) SAR realizing the High-Resolution Wide-Swath (HRWS) imaging mode. However, due to the influenceof satellite attitude error, the deviation of the DBF receiving beam direction always arises in such system. Thisleads to ghost targets appearing in the SAR image, when mapping the scenes (such as the seaport areas) withstrong scatterers. To address the problem, a matrix pencil method based DBF processing approach in elevationis presented. Firstly, according to the given threshold, the peak position of the strong scatterer is found fromthe range-compressed signals. Then, the direction of arrival angle of the strong scatterer is estimated using thematrix pencil method. Finally, based on the imaging geometry model, the DBF weighting vector is adjusted toensure the receiving beam to precisely point to the signal sources. Thereby, the interferences of ghost targets inSAR image can be eliminated effectively. The theoretical analysis is derived in detail, then it is validated bysimulation experiments.Key words: SAR; Multiple Elevation Beam (MEB); Digital Beam Forming (DBF); Matrix pencil method1    引言传统的单通道星载合成孔径雷达(SAR)系统由于受到最小天线面积的约束[1]，不能同时获取方位向高分辨率和距离向宽测绘带的SAR图像，其有限的成像能力越来越不能满足人们在遥感应用中对SAR图像质量的要求。随着多通道技术的发展和数字波束形成(DBF)技术的应用，配备多个接收通道的新体制星载SAR系统能够有效克服这一限制[2–5]。近年来，距离向多波束(MEB)体制的星载SAR系统[6–8]成为了各大科研机构在下一代星载SAR系统预研过程中的重点关注对象，如德宇航正在研制的双基星载SAR系统Tandem-L就应用了MEB体  收稿日期：2018-01-19；改回日期：2018-05-10；网络出版：2018-08-14*通信作者： 叶恺　yekai_seven@hotmail.com基金项目：国家重点研发计划(2017YFB0502700)，国家自然科学基金(61701479)Foundation Items: The National Key Research and DevelopmentProject of China (2017YFB0502700), The National Natural Sci-ence Foundation of China (61701479)eig (¢取复数实部和虚部的运算符号。)表示计算矩阵特征值的运算符号。µqµ0c= µq¡~µ根据强散射体所在距离门单元的位置和式(5)计算强散射体对应的下视角 ，并进一步计算实际的天线法线指向角： 。然后对DBF加权矢量进行调整，从而确保DBF接收波束中心和零陷位置指向准确。最后，利用调整后的DBF加权矢量对俯仰向多通道回波数据进行处理，分离提取各子测绘带回波并下传地面雷达站进行传统成像处理。N£1(N¡L)£(L + 1)O³(L + 1)3´L [N=3; N=2]L N=2O¡N3=8 + 3N2=4 + 3N=2 + 1¢N£K KN£N¢µ ¢µ¹KeO¡N3+ KN2+¹Ke¡N2+根据上述推导可知，本文方法是对一个维的单快拍回波信号矢量进行处理，计算量主要集中在对 维的数据矩阵 进行奇异值分解的运算上，其计算量约为 次复乘。由于参数 的取值通常在 范围内选取，当 取最大值 时，本文方法的计算量约为次复乘。文献[21]是采用Capon谱估计的方法进行目标波达角估计，需要利用多组快拍样本数据进行处理，因此该方法是对一个 维的信号矩阵进行处理，其中 表示快拍数目。文献[21]方法的计算量主要集中在以下2方面：(1) 对 维样本协方差矩阵的构造和求逆运算；(2) 将发射天线主瓣对应的角度范围以为间隔进行划分( 的选取与精度需求有关)，该方法需要计算每个划分的角度对应的功率谱密度。假定需要计算 个角度对应的功率谱密度，文献[21]方法的计算量约为N¢¢K ¹Ke次复乘。由于 和 的取值通常远大于N，因此本文方法的计算量要明显低于文献[21]方法。4    仿真试验4.1  点目标仿真实验本节将通过一个C波段星载MEB SAR系统来验证本文提出的处理方法。表1为相应的系统参数。首先进行点目标仿真实验。如图1所示，目标P1和目标P2的斜距分别为791.17 km和880.59 km，下视角分别为27.00°和35.13°，其中目标P1位于子测绘带1，目标P2位于子测绘带2。将目标P2的幅度设置为目标P1幅度的3倍，并在回波信号仿真中加入高斯白噪声。假定在仿真实验中预先已知天线法线指向角为26.00°，但是由于存在卫星姿态误差等因素的影响，实际的天线法线指向角为27.00°。图3(a)为没有经过俯仰向DBF处理的直接成像结  图 3 点目标仿真结果表 1  系统仿真参数轨道高度 700 km卫星速度 7504 m/s载频 5.4 GHz信号带宽 60 MHz脉冲宽度 22 μs天线高度 1.5 m俯仰向子孔径数目 23天线长度 10 m脉冲重复频率 1800 Hz第 1 1 期 叶  恺等：基于矩阵束方法的星载MEB SAR俯仰向DBF处理方法 2663测绘带回波在同一个雷达接收窗内相互重叠，当采用子测绘带1对应的成像参数进行处理时，目标P1的成像聚焦良好，目标P2的成像结果在方位向出现散焦。但是，由于目标P2幅度较大，这将导致鬼影目标出现在子测绘带1的成像结果中。采用预设的DBF加权矢量对混合回波进行分离处理，得到子测绘带1和子测绘带2的成像结果分别如图3(b)和图3(c)所示。由于采用的DBF加权矢量出现偏差，DBF接收波束中心和零陷位置指向都发生偏离，导致两个子测绘带成像结果中都出现干扰目标。针对这一问题，采用本文所提的基于矩阵束方法的俯仰向DBF处理方法。图4为利用矩阵束方法通过100次Monte Carlo实验得到的目标P2波达角估计值与实际值的偏差，根据实验数据计算得到目标P2的波达角估计期望值为修改为:8.1218°，均方根估计误差为0.0016°，因此可看出矩阵束方法能够准确地估计强散射体目标的波达角。按照上节阐述的处理方法对回波数据进行处理，最终得到图3(d)和图3(e)所示的子测绘带1和子测绘带2的成像结果，可见本文方法能够准确调整DBF加权矢量，有效消除鬼影目标的干扰。4.2  分布目标仿真实验为了进一步验证本文方法的有效性，本节将采用中国科学院电子学研究所航天微波遥感系统部的海洋机载SAR数据进行分布目标仿真试验。图5(a)和图5(b)分别对应子测绘带1和子测绘带2的SAR图像数据。根据星载MEB SAR成像几何模型生成多通道回波数据。图5(c)为采用子测绘带1成像参数  图 4 目标P2的波达角估计值偏差  图 5 分布目标仿真结果2664 电  子  与  信  息  学  报 第 40 卷得到的直接成像结果。如图中虚线框所示，由于2个子测绘带回波相互混叠，子测绘带1的成像结果中出现了子测绘带2的场景，因此需要利用DBF处理对混合回波数据进行分离提取。图5(d)为采用预设DBF加权矢量处理后的成像结果。由于采用存在偏差的DBF加权矢量对混合回波数据进行处理，在对子测绘带1回波数据进行提取时，DBF波束方向图的零陷指向出现偏离，导致子测绘带2中的强散射体目标回波未能得到有效抑制，因此SAR图像中出现了明显的鬼影目标干扰。采用本文方法对混合回波数据进行处理，得到子测绘带1和子测绘带2的成像结果分别如图5(e)和图5(f)所示。对比图5(c)和图5(d)可见，由于正确调整DBF加权矢量，SAR图像中的鬼影目标被明显消除，这也验证了所提处理方法的有效性。5    结束语针对由卫星姿态误差等影响因素导致的星载MEB SAR系统的DBF波束指向偏差问题，本文提出了一种基于矩阵束方法的星载MEB SAR俯仰向DBF处理方法。该处理方法主要包含以下重要步骤：首先对俯仰向子孔径回波信号进行匹配滤波处理并寻找强散射体目标所在位置；然后利用矩阵束方法准确估计强散射体目标的波达角，并根据星载MEB SAR系统的成像几何模型和回波信号模型调整DBF加权矢量，从而确保DBF接收波束中心和零陷位置指向准确；最后利用该DBF加权矢量分离提取各子测绘带回波，并进行传统成像处理。本文详细阐述了主要处理流程，并给出了点目标和分布目标的仿真实验。仿真结果验证了本文方法能够准确调整DBF加权矢量，有效消除鬼影目标的干扰。参 考 文 献FREEMAN A, JOHNSON W, HUNEYCUTT B, et al. The“myth” of the minimum SAR antenna area constraint[J].IEEE Transactions on Geoscience and Remote Sensing,2000, 38(1): 320–324. doi: 10.1109/36.823926.[1]SUESS M, GRAFMUELLER B, ZAHN R, et al. A novelhigh resolution, wide swath SAR system[C]. Proceedings ofthe International Geoscience and Remote SensingSymposium, Sydney, 2001: 1013–1015.[2]YOUNIS M, ALMEIDA F, LOPEZ-DEKKER P, et al.Techniques and modes for multi-channel SARinstruments[C]. Preceedings of European Conference onSynthetic Aperture Radar, Hamburg, Germany, 2016:812–817.[3]KRIEGER G, HUBER S, VILLANO M, et al. SIMO andMIMO system architectures and modes for high-resolution[4]ultra-wide-swath SAR imaging[C]. Preceedings of EuropeanConference on Synthetic Aperture Radar, Hamburg,Germany, 2016: 187–192.YOUNIS M, HUBER S, PATYUCHENKO A, et al.Performance comparison of reflector- and planar- antennabased digital beam-forming SAR[J]. International Journal ofAntenna and Propagation, 2009, 2009(6): 1–13. doi:10.1155/2009/614931.[5]RINCON R, FATOYINBO T, OSMANOGLU B, et al.Development of NASA’s next generation L-band digitalbeamforming synthetic aperture radar (DBSAR-2)[C].Proceedings of European Conference on Synthetic ApertureRadar, Hamburg, Germany, 2016: 1251–1254.[6]VILLANO M, KRIEGER G , and MOREIRA A. Advancedspaceborne SAR systems with planar antenna[C]. RadarConference, Seattle, USA, 2017: 152–156.[7]YAN Wei, YANG Xu, SUN Jia, et al. An airbornedemonstration for high-resolution wide-swath spacebornereflector SAR systems[C]. International Conference onFrontiers of Sensors Technologies, Shenzhen, China, 2017:350–353.[8]HUBER S, VILLANO M, YOUNIS M, et al. Tandem-L:design concepts for a next-generation Spaceborne SARsystem[C]. Proceedings of European Conference onSynthetic Aperture Radar, Hamburg, Germany, 2016:1237–1241.[9]TRIDON D, BACHMANN M, ZAN F, et al. Tandem-Lobservation concept-contributions and challenges ofsystematic monitoring of earth system dynamics[C]. The18th International Radar Symposium, Prague, CzechRepublic, 2017: 1–9.[10]MOREIRA A, KRIEGER G, HAJNSEK I, et al. Tandem-L:A highly innovative bistatic SAR mission for globalobservation of dynamic processes on the Earth’s surface[J].IEEE Geoscience and Remote Sensing Magazine, 2015, 3(2):8–23. doi: 10.1109/MGRS.2015.2437353.[11]NANNINI M, MARTONE M, RIZZOLI P, et al.Spaceborne demonstration of coherent SAR tomography forfuture companion satellite SAR missions[C]. IEEEInternational Geoscience and Remote Sensing Symposium,Fort Worth, USA, 2017: 129–132.[12]VILLANO M, KRIEGER G, and MOREIRA A. StaggeredSAR: High-resolution wide-swath imaging by continuousPRI variation[J]. IEEE Transactions on Geoscience andRemote Sensing, 2014, 52(7): 4462–4479. doi: 10.1109/TGRS.2013.2282192.[13]VILLANO M, KRIEGER G, JAGER M, et al. StaggeredSAR: Performance analysis and experiments with realdata[J]. IEEE Transactions on Geoscience and RemoteSensing, 2017, 55(11): 6617–6638. doi: 10.1109/TGRS.[14]第 1 1 期 叶  恺等：基于矩阵束方法的星载MEB SAR俯仰向DBF处理方法 26652017.2731047.KRIEGER G, GEBERT N, and MOREIRA A.Multidimensional waveform encoding: A new digitalbeamforming technique for synthetic aperture radar remotesensing[J]. IEEE Transactions on Geoscience and RemoteSensing, 2008, 46(1): 31–46. doi: 10.1109/TGRS.2007.905974.[15]KRIEGER G. MIMO-SAR: Opportunities and pitfalls[J].IEEE Transactions on Geoscience and Remote Sensing,2014, 52(5): 2628–2645. doi: 10.1109/TGRS.2013.2263934.[16]叶恺, 禹卫东, 王伟. 一种基于短偏移正交波形的MIMOSAR处理方案研究[J]. 雷达学报, 2017, 6(4): 376–387. doi:10.12000/JR17048.YE Kai, YU Weidong, and WANG Wei. Investigation onprocessing scheme for MIMO SAR with STSO chirpwaveforms[J]. Journal of Radars, 2017, 6(4): 376–387. doi:10.12000/JR17048.[17]林玉川, 张剑云, 武拥军, 等. 双基星载HRWS-SAR系统方位向信号重构的矩阵求逆算法[J]. 雷达学报, 2017, 6(4):388–396. doi: 10.12000 /JR17060.LIN Yuchuan, ZHANG Jianyun, WU Yongjun, et al. Matrixinversion method for azimuth reconstruction in bistaticspaceborne high-resolution wide-swath SAR system[J].Journal of Radars, 2017, 6(4): 388–396. doi: 10.12000 /JR17060.[18]赵庆超, 张毅, 王宇, 等. 基于多帧超分辨率的方位向多通道星载SAR非均匀采样信号重建方法[J]. 雷达学报, 2017, 6(4):408–419.ZHAO Qingchao, ZHANG Yi, WANG Yu, et al. Signalreconstruction approach for multichannel SAR in azimuthbased on multiframe super resolution[J]. Journal of Radars,[19]2017, 6(4): 408–419.VAN T H L. Optimum Array Processing Part IV ofDetection, Estimation, and Modulation Theory[M]. NewYork: John Wiley & Sons, 2002: 428–669.[20]MAKHOUL V. Adaptive digital beam-forming for high-resolution wide-swath synthetic aperture radar[D]. [Masterdissertation], Polytechnic University of Catalonia, 2009:23–91.[21]BORDONI F, YOUNIS M, MAKHOUL V, et al. Adaptivedigital beamforming algorithm for high-resolution, wide-swath synthetic aperture radar[C]. International RadarSymposium, Hamburg, Germany, 2009: 1–5.[22]KHAN M and TUFAIL M. Comparative analysis of variousmatrix pencil methods for direction of arrival estimation[C].International Conference on Image Analysis and SignalProcessing, Zhejiang, China, 2010: 1–6.[23]SARKAR T and PEREIRA O. Using the matrix pencilmethod to estimate the parameters of a sum of complexexponentials[J]. IEEE Antennas an Propagation Magazine,1995, 37(1): 48–55. doi: 10.1109/74.370583.[24]HUA Y and SARKAR T. On SVD for estimatinggeneralized eigenvalues of singular matrix pencil in noise[J].IEEE Transaction on Signal Processing, 1991, 39(4):892–900. doi: 10.1109/78.80911.[25]叶   恺：男，1988年生，博士生，研究方向为新体制星载SAR系统设计和信号处理.禹卫东：男，1969年生，研究员，博士生导师，研究方向为合成孔径雷达系统设计和信号处理.王   伟：男，1985年生，博士生，研究方向为新体制星载SAR系统设计和信号处理.2666 电  子  与  信  息  学  报 第 40 卷