星载面阵SAR多维波形编码技术研究
2022-09-12
来源:爱问旅游网
第32卷第12期 2011年12月 宇航学报 Vo1.32 No12 .Journal of Astronautics December 201 1 星载面阵SAR多维波形编码技术研究 曾祥能 ,何峰 ,张永顺 ,董臻 (1.空军工程大学导弹学院,三原713800;2.国防科技大学电子科学与工程学院,长沙410073) 摘要:为解决传统高分辨和宽测绘带以及高信噪比和宽测绘带之间的矛盾,采用面阵合成孔径雷达(SAR) 多维波形编码(Multi—Dimensional Waveform Coding,MDWC)技术,提出了将面阵分为多行多列的子孔径,将距离维 测绘带分成多个子测绘带进行脉内扫描照射,同时对子测绘带内子脉冲信号进行频分编码,并采用低PRF发射,接 收端则利用方位向多个相位中心解多普勒模糊,获得宽测绘带、高分辨与高信噪比等多项性能,仿真试验表明 方 法的有效性。 关键词:合成孔径雷达;多维波形编码;宽测绘带成像;解多普勒模糊 中图分类号:TN957 文献标识码:A 文章编号:1000—1328(201 1)12-2560-07 DOI:10.3873/j.issn.1000—1328.201 1.12.014 Study on Multi-Dimensional Waveform Coding Technique for Spaceborne Planar Array SAR System ZENG Xiang.neng 一,HE Feng ,ZHANG Yong.shun ,DONG Zhen (1.Missile Institute of Air Force Engineering University,Sanyuan 7 13800,China; 2.School of Electronic Engineering and Science,National Univ.of Defense Tech.,Changsha 410073,China) Abstract:In order to solve contradictions between high azimuth resolution and wide swath or signal to noise ratio (SNR),a muhi—dimensional waveform coding(MDWC)technique for planar array SAR system is proposed,in which the planar array is divided into multiple sub—apertures as rows and columns,meanwhile the imaging swath in range dimension is divided into multiple sub—swath,the sub—swaths are scanned in intra—pulse,and the sub—pulse signal for given sub—swath is coded in frequency domain,the whole system works in low PRF.On the receiver side,multiple phase centers in azimuth is utilized to resolve Doppler ambiguities,the new system can acquire obvious abilities such as wide swath,high resolutions and high signal to noise ratio,and SO on.The simulation results verify the effectiveness of the proposed method. Key words:Synthetic aperture radar;Multi—dimensional waveform coding;Wide swath imaging;Resolving Doppler ambiguities 0引 言 需要PRF尽可能低,二者存在矛盾而不可兼得,针 对高分辨宽测绘带(High Resolution Wide Swath, HRWS)需求,学者们提出了许多技术途径,文献[1] 提出采用高PRF,距离维采用数字波束形成解模糊 来获得宽测绘带,文献[2]采用中PRF,距离维、方 合成孑L径雷达(Synthetic Aperture Radar,SAR) 在对地观测方面得到越来越广泛的应用,一方面,为 了全面准确地获取地面场景信息,要求成像分辨率、 信噪比尽可能高,方位高分辨则要求脉冲重复频率 (Pulse Repetition Frequency,PRF)高而避免方位模 糊;另一方面,为了缩短再访时间,要求宽的测绘带, 收稿日期:2010 10—15; 修回日期:2011-03-19 基金项目:武器装备预研基金项目(9140A21020910KG01) 位维均模糊,接收端在距离维、方位维同时形成多波 束来解模糊,方法计算量大,文献[3]采用低PRF, 保证距离不模糊,方位向利用波束形成来解模糊,文 第12期 曾祥能等:星载面阵SAR多维波形编码技术研究 2561 献[4]利用面阵沿航向子孔径解方位模糊实现 HRWS成像,距离维多个子孔径相干积累来提高回 波信噪比,但其只采用单孔径发射而通道利用率很 低,文献[5]则在此基础上将测绘区划分为多个子 测绘带,运用脉内扫描进一步获得宽测绘带和信噪 比提升。Krieger[6]等提出的多维波形编码技术可 解决高分辨一宽测绘带矛盾,本文通过距离维多波 束多波形技术获得能量的灵活分配,同时获得信噪 比提升,使得重点区域得到充分观测而获得准确全 面的信息,并避免了传统HRWS系统成像回波能量 由近及远递减的趋势,本系统采用频分编码的多波 形技术在接收端能将各子脉冲彻底区分,避免了距 离模糊,而各子测绘带对应的照射信号可灵活应用, 从而实现不同分辨率设置等功能,且能获得一定的 抗干扰能力,接收端采用全阵列接收,并运用DBF 形成零陷,可进一步提高系统抗干扰性能。 定义(・) 为矩阵转置,(・)“为矩阵共轭转 置,(・) 为矩阵的伪逆,rect(・)为窗函数。 发射脉冲 接收窗 豳 r 广1 jr 广 J) 豳圜 _,厂 广1 ・● ———————+ A 图1 面阵SAR距离维波形编码示意图 Fig.1 Sketch of multi—dimensional waveform coding for plane array SAR 1发射信号模型 设发射信号时宽为 ,将其分为 个子脉冲,子 脉冲宽度为 (z=1,…,,J),利用发射阵列加权控 制波束指向,使其在不同子脉冲时间内扫描不同的 子测绘带,也可称作“脉内扫描”。如图1所示,方位 向为 轴,距离向为Y轴,将大型面阵均匀分为N X 个子孔径,方位向中间列的子孔径以低PRF进行 发射,在一个PRI(脉冲重复间隔)内,通过发射阵元 权矢量凋整,控制发射波束在距离向由远及近依次 照射 个子测绘带(图中以L=3为例),每个子测 绘带的子脉冲对应独立的发射波形。 3 显然,发射脉冲总的时长为△丁,=∑ + , 为第2个发射子脉冲后的空歇时间,用以调整各 子脉冲的相对延时,令接收窗时长为△ ,则满足: △7-r+△7R=PRI (1) 这里忽略了发射一接收转换时间。系统采用低PRF 可获取宽测绘带信息,根据场景斜距适当调整 和 ,可使所有回波聚集在小段接收窗口内,即缩小 △r ,在PRI一定的情况下,可相应增长△丁 ,总结 这项技术的优势如下: (1)总的接收窗缩短,从而减少了数据量,利于 星载平台工程实现; (2)缩短了接收窗长度,留给发射端更长的工 作时间,以利于发射信号占空比增加,可提高回波信 噪比或系统作用距离。 (3)各子脉冲的发射时间 可灵活调整,利于 对敏感地区、弱目标所在测绘带增加照射时间,利 于场景成像或运动目标检测。 在图1所示模型中,令距离向的三个子波束对 应的波形为线性调频信号,且其调频斜率为 ,则第 Z个子脉冲波形的调频带宽为B = ,发射中心频 率为: f = +∑B —Bz/2 (2) 式中 为载频,显然各子脉冲信号频谱互不重叠。 设方位向的孔径间距为d ,距离向的孔径间距 为d ,这里以Ⅳ、 为奇数为例,则可等效为N×M个 等效相位中心 的自发自收形式,第n行第,n列子 孔径对应的等效相位中心的坐标(疋, ,Zo)为 : +(m一 ) , Yo+( 一 ) dy:, (n=1,…,Ⅳ,m=l,…, ), ( , , )为发射波束中心。 在第z个子脉冲作用期,发射信号可表示为: ,。、 sin(NTrdy(sin0z—sin0)/(2A)) ’ rect( )唧(j2榔 )+ j ( 一△ ) ) (3) 第12期 曾祥能等:星载面阵SAR多维波形编码技术研究 2563  ̄/(Yo—Y ) +z 为散射点( ,Y , )到航线的垂 直距离;则有 R (£ ) R +( + — 2) (7) 式(8)中最后一个指数项是由于相位中心空间 位置不同形成的附加相位项,可利用此相位解多普 勒模糊,将第z个子脉冲信号对应的 个相位中心 的差异部分写成矢量形式为 对(6)式进行方位傅里叶变换,得到多普勒模 糊信号为: s…2( , ) =)=fexp J ),exp J 唧A(0f)sinc(t 一△ 一丁( ))・ (j )] (9) 现假设有K个多普勒模糊,当K为奇数时,由 脉冲重复频率为PRF得各模糊的频率中心为: 唧-jp 唧),I’. exp/ 一 j 2 "rR.t f.2]‘ /‘ 一 ,= ・[一 ,…,_1’o,1,…, r (10) (j ) =2vf/c。令 (8) 同理当K为偶数时: 式中丁( )为测绘带内某散射点的瞬间延时, ∈ (一PRF/2,PRF/2), s (t , ) =,= ’[一 ,…,一号, ・, ] (11) 则对于以较低PRF采集的多普勒频带内任一 频率点 ,各模糊频率组成如下矢量: =A (0,)sinc(t 一△ 一 ( )). rect(  ̄‘exp(一j 2 ̄rf" xt),・ 722……… 。f+ ・IK 。 (12) 其中, 在[一PRF/2,PRF/2]范围内。 exp/ -rRszf ̄j 47m)‘exp/ a--ot a a ‘一 现利用空域滤波来解模糊的多普勒,假设一组 权向量w =[W W ,…W ,…, r,对于任一 不同的 ,都可通过下列公式得到不同组的权向量 解多普勒模糊。  ̄]J Snml ,…p(j )c (f(1)+ ) z ( )=[Zl(,(1)+ ),…,Zl(,( )+ ),…,Zl(,(K)+ )] (f(2)+ )X (,(2)+厂n) (,( )+ ) exp .2叮T J——‘ (,(1)+ ) (,(K)+ )砭 (f(K)+ ) (13) (f(1)+ ) Zz( )・ :日 , k=1,2,…,K(14) (,(2)+ ) 定义 hf =0。 ==[h h ,…,h ] ,其中h =1, (Zf( )) , =1,2,…, (15) 多普勒信号之和: )= K s ))exp( ) (17) 第z个子脉冲信号对应的 个相位中心的信号 在距离时域方位多普勒域的矢量形式为 f这样,(16)式就可表示为 Sf(t , )=[5¨( , ),s2.f(t , ),…, S删( , )] K = (t , )=[s¨(t , ),s .f( , ),…, (t , )] (16) 各接收阵元接收多普勒频率点信号为多个模糊 ))[eXp( : , 2564 宇航学报 第32卷 唧( )] (18) 测绘带内的目标回波可被压缩在一个200 ̄s长的接 收窗内被接收,则测绘场景宽度为84kin,远远大于 =∑s , ( ))z ( )) 当前的条带式测绘场景宽度,配以发射端的大带宽 利用(15)式得到的权向量对(18)式的信号进 行滤波,得到不模糊的 ( )频率分量: f(t , ( ))= f(t , ) =∑si(t , (p))z ( (p)) =∑s , (p)) (k:1,2,…,K) (19) 最后将K个不模糊的频率分量拼接起来就可 得到解模糊后的信号形式: s;(t , ), ∈[一PRF×K/2,PRF×K/2] (20) 完成多普勒解模糊后,距离时域方位多普勒域 的信号形式为 sf(t , ) =A(0t)sinc(t 一△ 一丁 ))・ e。 xp-p —j : 葡JDt J ̄l,)e。x p(I一 )cJ z2 4仿真实验 本节将主要检验解模糊算法的有效性以及脉内 扫描面阵SAR系统性能分析的正确性,主要验证该 系统的高信噪比性能。 (1)方位解多普勒模糊算法分析 假设大型面阵分为3×3个子孔径,孔径间距 d ×d 为10m×0.8m,各子孔径坐标为( , ,0), X =10(m一2), =0.8(n一2),m,n=1,2,3;将 发射脉冲分为3个子脉冲,子脉冲时宽Tz=100 ̄s(z :1,2,3)。平台速度 =7500m/s,PRF=500Hz,而 根据参数可算得多普勒带宽约为B =2v/d = 1500Hz,因此,方位多普勒谱模糊了3次。利用面 阵在俯仰向形成窄波束在不同子脉冲时间内由远及 近照射子测绘带,三个子测绘带近端与沿航线地面 投影的垂直距离由近及远依次为400km、429km、 457km,场景远端的垂直距离为484km,保证三个子 子脉冲信号,即可获得高分辨宽测绘带成像。设置 3个散射点坐标分别为P.(0,480,一700)km, P2(0,450,一700)km,P3(0,420,一700)km,对点 目标回波信号的距离时域方位多普勒域成像如图4 所示。由于较低的PRF引起了3次多普勒模糊,因 此在图4(a)中距离多普勒域内的谱线不连续,在图 4(b)中场景成像中方位位置两边出现了“模糊点”, 由于在合成孑L径两侧的距离弯曲较大,引起方位向 的模糊点脉压很差,且距离较近的散射点受影响更 大。图4(C)显示通过匹配滤波能彻底将距离向各 散射点回波分开,获得良好的距离向脉压,图4(d) 则显示回波在方位向模糊点严重,与目标散射点幅 度接近。对应的,图5则为利用方位向信息矩阵解 模糊后的成像情况,可见方位模糊已经被消除,且距 离向的脉压输出未受任何影响。 图4接收信号多普勒模糊成像示意图 Fig.4 Receptive signal imaging with Doppler ambiguity (2)系统性能分析 采用面阵多维波形编码技术可以实现高分辨宽 测绘带成像,现讨论其对改善系统信噪比的效果。 由于回波信噪比与天线面积及发射带宽成正比,传 统的HRwS系统采用单孔径发射,全孔径接收,其 发射天线面积d ×d ,接收天线面积Md ×Nd 而 本文采用的中间列在俯仰向扫描发射子脉冲信号, 其发射天线面积d N×d ,接收天线面积Md ×Nd 令二者的PRF相同,对相同幅宽的大场景进行成 第12期 曾祥能等:星载面阵SAR多维波形编码技术研究 2565 ∞ 。 印 ∞ (a)距离多普勒域成像 图5 接收信号解模糊后成像示意图 Fig.5 Receptive signal imaging after removing of Dopp|er ambiguity 像,令场景的回波接收时间为 ,则HRWS的接 收窗近似为 + ,本文中的接收窗近似为 / + / ( 为子脉冲长度和),根据式(1)的分析可 知有PRI=T+ + h。=T+ / + hn/ 近似 成立,则二者SNR比值Ratio=SNRMDwc/SNRHRws :(Ⅳ )/(LT)=(2N(L・PRI— h。))/( ( + 1)(PRI—re。h。)),通常脉冲占空比T/PRI在1%~ 10%,取为10%,则最大场景宽度可达到To = 0.75PRI,图6作出MDWC系统的SNR相对HRWS 的改善情况随子脉冲数 的变化。 图6 Ratio与子测绘带数关系 (N=3—8表示子孔径数) Fig.6 Relationship between ratio and amount of sub—swath(N denotes number of sub—swath) 多维波形编码技术利于采用不同的编码波形信 号对子测绘带进行照射成像,如第1个子测绘带内 的两个散射点q。(0,425.005,一700)km,g:(0, 425.045,一700)km,与第2个子测绘带内的两个散 射点g2(0,455.005,一700)km,g2(0,455.045, 一700)km,第1个子测绘带内照射信号为带宽Br= 3MHz的调频信号,可实现场景内快速探测,进行实 时监视,而对于第2个子测绘带内,令设定的目标在 此范围内,需进行精细识别,采用带宽为10MHz的 宽带信号照射,场景内的距离分辨率约为18m,显 然能将g 、q 分开,如图7(b)所示。进一步,可通 过增长对重点区域的照射时间,以解决弱散射点的 成像问题。 】,/kill (a)第1子测绘带内散射点脉压图 .1 蜘 Y/kin (b)第2子测绘带内散射点脉压图 图7不同子测绘带采用不同波形子脉冲成像 Fig.7 Different sub—swath within different kind of sub—pulse waveform for imaging 5结论 面阵SAR多维波形编码技术能获得高信噪比 下的高分辨宽测绘带成像,而且通过距离向的多个 子测绘带划分进行扫描照射,能有效解决大场景成 像的边缘能量分布不均等问题,同时能使距离分辨 率不随场景的距离变化而变化,其灵活的子测绘带 能量分配可使弱散射区、重点区域得到更清晰成像。 多维波形编码技术在实现星载SAR的多功能成像 方面具有广阔的应用前景,将继续对面阵SAR的方 2566 宇航学报 第32卷 位向多维编码技术、距离一方位维联合多维编码技 术展开研究。 参考文献 Callaghan G D,Long I D.Wide—swath space-borne SAR and range ambiguity[C].The Electronics and Communications Division of the Institution of Electrical Engineers.Radar 97: Edinburgh,1997:248—252. [2] Nathan A,Goodman,Sih C L,Devindran R,et a1.Processing of muhiple-receiver spaceborne arrays for wide—area SAR[J].IEEE Trans.on GRS,2002,40(2):841—852. [3] Li Z F,Bao Z,Wang H Y,et a1.Performance improvement for constellation SAR using signal processing techniques[J].IEEE Trans.on AES.2006,42(2):436—452. [4] Suess M,Grafmueller B,Zahn R.A novel high resolution wide swath SAR system[C].Proceedings of International Geoscience and Remote Sensing Symposium(IGARSS).Sydney,NSW, Australia,2001.1013—1015. [5] 武其松,石洪竺,邢孟道,等.高信噪比、高分辨宽测绘带成 像[J].航空学报,2010,31(5):958~967[Wu Qi—song,Shi Hong—zhu,Xing Meng—dao et a1.High signai to noise ratio and high resolution wide swath imaging[J j Acta Aeronautica ET Astronantica Sinica,2010,31(5):958—967.] [6] Krieger G,Gebert N,Moreira A.Muhidimensioaal wavcfoIt/] encoding:a new digital beamfomfing technique for synthetic aperture radar remote sensing[J].IEEE Trans on Geoscience and Remote Sensing,2008,46(1):31—45. l7] Fischer C,Schaefer C,Heer C.Technology development for the HRWS SAR[c J.International Geoscienee and Remote Sensing Symposium(IGARSS).Barcelona,Spain,2007:1024—1028 作者简介:曾祥能(1983一),男,博士生,主要研究方向为雷达信 号波形编码设计、接收信号处理、波束控制及其最优化设计 通信地址:湖南长沙国防科技大学电子科学与工程学院2所2 室(410073) E—mail:zengxiangneng103@163.eom (编辑:余未)