空间轨道目标的逆合成孔径雷达成像质量分析

周叶剑 张磊 王虹现 邢孟道

引用本文:
Citation:

空间轨道目标的逆合成孔径雷达成像质量分析

    作者简介: 周叶剑(1993–),男,浙江人,2015年获西安电子科技大学工学学士学位,现为西安电子科技大学电子工程学院博士研究生,主要研究方向为ISAR成像与图像解译。E-mail: zhouyejian25@163.com;张 磊(1984–),男,浙江人,2011年获西安电子科技大学工学博士学位,现为西安电子科技大学电子工程学院副教授,主要研究方向为SAR/ISAR成像与运动补偿。E-mail: leizhang@xidian. edu.cn.
    通讯作者: 张磊, leizhang@xidian.edu.cn
  • 基金项目:

    国家自然科学基金(61301280, 61301293)

  • 中图分类号: TN957.52+4

Performance Analysis on ISAR Imaging of Space Targets

    Corresponding author: Zhang Lei, leizhang@xidian.edu.cn ;
  • Fund Project: The National Natural Science Foundation of China (61301280, 61301293)

    CLC number: TN957.52+4

  • 摘要: 空间目标逆合成孔径雷达(ISAR)成像的方位维相干积累增益取决于目标运动状态信息,关于目标运动分析的缺失将直接影响空间目标ISAR观测中成像质量的预估与系统参数的设计。该文提出根据空间目标轨道参数计算其相对雷达视线运动状态,推导一定相干积累角ISAR成像的方位相干积累简化公式,实现基于相干积累角的ISAR成像时间段选择,可有效保证成像方位分辨率;同时详细分析目标轨道高度等参数对ISAR系统回波接收功率、成像质量的影响,为空间目标ISAR成像雷达体制设计提供了设计依据。理论推导和仿真实验均验证了空间目标ISAR成像中随着目标轨道高度升高,单次回波信噪比降低的损失可通过方位相干积累增益的提升实现部分弥补。该文工作可为空间目标ISAR成像体制和处理设计提供了理论基础和指标设计依据。
  • 图 1  空间观测几何模型及成像模型

    Figure 1.  Observation and imaging model of satellite targets

    图 2  雷达斜距变化Δr求解模型

    Figure 2.  Calculation model of change of radar range

    图 3  Δr求解近似模型与实际模型

    Figure 3.  Approximate and actual model of Δr

    图 4  不同轨道高度下目标转角速度计算结果

    Figure 4.  Calculation of target′s rotate speed at different heights

    图 5  归一化回波信号接收功率随目标轨道半径变化曲线

    Figure 5.  Normalization curve of received power changing in different orbit radii

    图 6  目标图像距离、方位剖面图

    Figure 6.  Range and azimuth profiles of target image

    图 7  相近观测视角下,不同轨道高度成像结果对比图

    Figure 7.  Comparison result of target imaging at different heights with similar LOS parameters

    表 1  实验ISAR系统主要参数

    Table 1.  Main parameters of ISAR system

    参数 数值
    载频 16.7 GHz
    带宽 1 GHz
    方位向分辨率 0.18 m
    距离向分辨率 0.15 m
    脉冲重复频率 200 Hz
    下载: 导出CSV

    表 2  实验地基ISAR观测站位置

    Table 2.  Position parameters of radar sites

    地点 经纬度
    库尔勒 41.5°N, 86.8°E
    北京 39.9°N, 116.4°E
    西安 31.1°N, 108.4°E
    下载: 导出CSV

    表 3  不同轨道高度图像质量评价

    Table 3.  Comparison result of imaging quality at different heights

    轨道高度(km) 成像时间(s) 脉冲积累数 TNR 脉冲积累数/TNR
    791 9.94 1988 6.82 291.58
    1200 12.40 2484 8.61 288.58
    1800 16.82 3364 11.45 293.80
    下载: 导出CSV
  • [1] Xing Mengdao and Bao Zheng. High-resolution ISAR imaging of high speed moving targets[C].IEE Proceedings- Radar,Sonar and Navigation, 2005, 152(2): 58–67.
    [2] Bai Xueru, Xing Mengdao, Zhou Feng,et al.. High-resolution 3D imaging of spinning space debris[J].IEEE Transactions on Geoscience and Remote Sensing, 2009, 47(7): 2352–2362. DOI: 10.1109/TGRS.2008.2010854.
    [3] 王琦. 空间目标 ISAR 成像的研究[D]. [博士论文], 西安电子科技大学, 2007.Wang Qi. Study of ISAR imaging for space targets[D]. [Ph.D. dissertation], Xidian University, 2007.
    [4] 王洋, 金胜, 黄璐. 空间目标双基地雷达ISAR成像技术研究[J]. 雷达科学与技术, 2015, 13(5): 485–495.Wang Yang, Jin Sheng, and Huang Lu. Bistatic radar ISAR imaging of space target[J].Radar Science and Technology, 2015, 13(5): 485–495.
    [5] 王珊, 季伟, 陈娟, 等. 基于运动目标检测的多脉冲相干积累算法[J]. 机械与电子, 2010, 9: 490–491.Wang Shan, Ji Wei, Chen Juan,et al.. Coherent integration algorithm of high velocity moving target[J].Science and Technology Information, 2010, 9: 490–491.
    [6] 庞存锁, 张南, 侯慧玲. 高速微弱目标长时间积累对雷达方程的影响[J]. 火控雷达技术, 2015, 44(1): 7–9.Pang Cunsuo, Zhang Nan, and Hou Huilin. Radar equation for high-speed weak target during long time integration[J].Fire Control Radar Technology, 2015, 44(1): 7–9.
    [7] 李亮, 洪峻, 明峰, 等. 电离层时空变化对中高轨SAR成像质量的影响分析[J]. 电子与信息学报, 2014, 36(4): 915–922. DOI: 10.3724/SP.J.1146.2013.00859.Li Liang, Hong Jun, Ming Feng,et al.. Study on ionospheric effects induced by spatio-temporal variability on medium-earth-orbit SAR imaging quality[J].Journal of Electronics &Information Technology, 2014, 36(4): 915–922. DOI: 10.3724/SP.J.1146.2013.00859.
    [8] 黄小红, 文贡坚. L波段雷达电离层高速运动目标ISAR成像补偿方法[J]. 电子与信息学报, 2015, 12(37): 2971–2976. DOI: 10.999/JEIT150646.Huang Xiaohong and Wen Gongjian. Compensating method of L-band radar ISAR imaging for ionosphereic target with high-velocity[J].Journal of Electronics &Information Technology, 2015, 12(37): 2971–2976. DOI: 10.999/ JEIT150646.
    [9] 袁孝康. 合成孔径雷达方程[J]. 上海航天, 2002, 3: 1–5.Yuan Xiaokang. Synthetic aperture radar equation[J].Aerospace Shanghai, 2002, 3: 1–5.
    [10] 范录宏. 逆合成孔径雷达成像与干扰技术研究[D]. [博士论文], 电子科技大学, 2006.Fan Luhong. Research on imaging and jamming of inverse synthetic aperture radar[D]. [Ph.D. dissertation], University of Electronic Science and Technology of China, 2006.
    [11] 许然. 提高雷达成像质量的若干新体制和新方法研究[D]. [博士论文], 西安电子科技大学, 2015.Xu Ran. Study on new systems and techniques for improving radar imaging performances[D]. [Ph.D. dissertation], Xidian University, 2015.
    [12] Cao X, Su F, Sun H,et al.. Space debris observation via space-based ISAR imaging[C]. International Conference on Microwave and Millimeter Wave Technology, Guilin, China, April 18–21, 2007: 1–5.
    [13] Skolnik MI. Radar Handbook (Second edition)[M]. New York: Mcgraw-Hall Publishing Company, 1990.
    [14] 张云彬, 张永生. 近圆轨道遥感卫星星下点轨迹的计算[J]. 测绘学院学报, 2001, 18(4): 257–259.Zhang Yunbin and Zhang Yongsheng. Calculation of sub-satellite track of remote sensing satellite in nearly round orbit[J].Journal of Institute of Surveying and Mapping, 2001, 18(4): 257–259.
    [15] 任锴. 导航卫星精密定轨理论与方法研究[D]. [博士论文], 信息工程大学, 2015.Ren Kai. Research on theory and methodology of precise orbit determination of navigation satellites[D]. [Ph.D. dissertation], PLA Information Engineering University, 2015.
    [16] 黄艳, 张永利, 刘志铭. 一种基于点目标的雷达影像质量评价方法[J]. 测绘工程, 2012, 21(1): 30–34.Huang Yan, Zhang Yongli, and Liu Zhiming. A method of SAR image quality evaluation based on pinpoint target[J].Engineering of Surveying and Mapping, 2012, 21(1): 30–34.
  • [1] 蒋铁珍肖文书李大圣廖同庆 . 基于星载机会源的空间目标外辐射源雷达探测技术可行性研究. 雷达学报, 2014, 3(6): 711-719. doi: 10.12000/JR14080
    [2] 仲利华胡东辉丁赤飚张问一 . 一种稀疏孔径下大尺寸目标的ISAR 成像方法. 雷达学报, 2012, 1(3): 292-300. doi: 10.3724/SP.J.1300.2012.20033
    [3] 黄剑胡卫东 . 基于贝叶斯框架的空间群目标跟踪技术. 雷达学报, 2013, 2(1): 86-96. doi: 10.3724/SP.J.1300.2012.20079
    [4] 谷文堃王党卫马晓岩 . 分布式MIMO-ISAR子图像融合方法. 雷达学报, 2017, 6(1): 90-97. doi: 10.12000/JR16042
    [5] 陈文峰李少东杨军马晓岩 . 基于线性Bregman迭代类的多量测向量ISAR成像算法研究. 雷达学报, 2016, 5(4): 389-401. doi: 10.12000/JR16057
    [6] 邓冬虎张群罗迎李松朱仁飞 . 双基地ISAR 系统中分辨率分析及微多普勒效应研究(英文). 雷达学报, 2013, 2(2): 152-167. doi: 10.3724/SP.J.1300.2013.13039
    [7] 李道京胡 烜 . 合成孔径激光雷达光学系统和作用距离分析. 雷达学报, 2018, 7(2): 263-274. doi: 10.12000/JR18017
    [8] 李海英张珊珊李世强张华春 . 环境一号C 卫星合成孔径雷达相干性分析. 雷达学报, 2014, 3(3): 320-325. doi: 10.3724/SP.J.1300.2014.13060
    [9] 张月婷丁赤飚王宏琦胡东辉 . 分析圆柱油罐目标合成孔径雷达图像的新方法. 雷达学报, 2012, 1(2): 190-195. doi: 10.3724/SP.J.1300.2012.20020
    [10] 吴一戎 . 多维度合成孔径雷达成像概念. 雷达学报, 2013, 2(2): 135-142. doi: 10.3724/SP.J.1300.2013.13047
    [11] 赵雨露张群英李超纪奕才方广有 . 视频合成孔径雷达振动误差分析及补偿方案研究. 雷达学报, 2015, 4(2): 230-239. doi: 10.12000/JR14153
    [12] 詹学丽王岩飞王超李和平 . 一种用于合成孔径雷达的数字去斜方法. 雷达学报, 2015, 4(4): 474-480. doi: 10.12000/JR14117
    [13] 林世斌李悦丽严少石周智敏 . 平地假设对合成孔径雷达时域算法成像质量的影响研究. 雷达学报, 2012, 1(3): 309-313. doi: 10.3724/SP.J.1300.2012.20035
    [14] 王岩飞刘畅詹学丽韩松 . 无人机载合成孔径雷达系统技术与应用. 雷达学报, 2016, 5(4): 333-349. doi: 10.12000/JR16089
    [15] 曾操梁思嘉王威徐青 . 基于频率步进信号的旋转式合成孔径雷达成像方法. 雷达学报, 2014, 3(4): 401-408. doi: 10.3724/SP.J.1300.2014.14043
    [16] 路满宋红军罗运华 . 基于调频连续波信号的圆弧式合成孔径雷达成像方法. 雷达学报, 2016, 5(4): 425-433. doi: 10.12000/JR16007
    [17] 曾涛 . 双基地合成孔径雷达发展现状与趋势分析. 雷达学报, 2012, 1(4): 329-341. doi: 10.3724/SP.J.1300.2012.20093
    [18] 杨汝良戴博伟李海英 . 极化合成孔径雷达极化层次和系统工作方式. 雷达学报, 2016, 5(2): 132-142. doi: 10.12000/JR16013
    [19] 金添 . 叶簇穿透合成孔径雷达增强成像方法. 雷达学报, 2015, 4(5): 503-508. doi: 10.12000/JR15114
    [20] 张珂殊潘洁王然李光祚王宁吴一戎 . 大幅宽激光合成孔径雷达成像技术研究. 雷达学报, 2017, 6(1): 1-10. doi: 10.12000/JR16152
  • 加载中
图(7)表(3)
计量
  • 文章访问数:  865
  • HTML浏览量:  26
  • PDF下载量:  732
  • 被引次数: 0
出版历程
  • 收稿日期:  2016-12-01
  • 录用日期:  2017-03-03
  • 刊出日期:  2017-02-28

空间轨道目标的逆合成孔径雷达成像质量分析

    通讯作者: 张磊, leizhang@xidian.edu.cn
    作者简介: 周叶剑(1993–),男,浙江人,2015年获西安电子科技大学工学学士学位,现为西安电子科技大学电子工程学院博士研究生,主要研究方向为ISAR成像与图像解译。E-mail: zhouyejian25@163.com;张 磊(1984–),男,浙江人,2011年获西安电子科技大学工学博士学位,现为西安电子科技大学电子工程学院副教授,主要研究方向为SAR/ISAR成像与运动补偿。E-mail: leizhang@xidian. edu.cn
  • ①. 西安电子科技大学雷达信号处理国家重点实验室   西安   710071
  • ②. 西安电子科技大学信息感知技术协同创新中心   西安   710071
基金项目:  国家自然科学基金(61301280, 61301293)

摘要: 空间目标逆合成孔径雷达(ISAR)成像的方位维相干积累增益取决于目标运动状态信息,关于目标运动分析的缺失将直接影响空间目标ISAR观测中成像质量的预估与系统参数的设计。该文提出根据空间目标轨道参数计算其相对雷达视线运动状态,推导一定相干积累角ISAR成像的方位相干积累简化公式,实现基于相干积累角的ISAR成像时间段选择,可有效保证成像方位分辨率;同时详细分析目标轨道高度等参数对ISAR系统回波接收功率、成像质量的影响,为空间目标ISAR成像雷达体制设计提供了设计依据。理论推导和仿真实验均验证了空间目标ISAR成像中随着目标轨道高度升高,单次回波信噪比降低的损失可通过方位相干积累增益的提升实现部分弥补。该文工作可为空间目标ISAR成像体制和处理设计提供了理论基础和指标设计依据。

English Abstract

    • 随着雷达空间探测技术的发展,空间目标的高精度成像处理已成为空间探测任务的重要组成部分。对空间目标的雷达观测成像主要采用地基逆合成孔径雷达(Inverse Synthetic Aperture Radar, ISAR)体制实现,目前,针对低轨目标的成像技术比较完善[14],但从空间目标运动特性上进行ISAR成像体制和信号处理方法设计来提升雷达工作性能的研究偏少,需要直接面对若干问题:(1)根据传统雷达方程可知,由于回波信噪比与雷达目标距离4次方成反比关系,目标轨道高度升高将导致单次回波信噪比急剧下降;(2)传统ISAR成像处理通常将目标认为是“完全”非合作目标,ISAR成像体制与信号处理设计均未能紧密结合目标轨道参数。与空中气动目标和海面舰船目标相比,空间轨道目标的运动相对平稳、可预测性强,而且目标姿态通常严格受控,目标相对雷达视线角变化通常可以根据轨道信息精确解算。此外,与合成孔径雷达(SAR)成像类似,ISAR成像过程也是积累多帧脉冲串相干聚焦实现方位成像,方位相干积累增益可有效地提高成像信噪比质量[5]。针对SAR系统中相干积累对成像质量的分析已较为完善[68],考虑方位相干积累增益后,条带SAR雷达方程中成像信噪比将与作用距离的3次方成反比,因此SAR体制参数设计与传统雷达方程是有明显差异的。从成像原理上来看,ISAR成像相干积累角也能带来明显增益,但其成像几何较SAR不同,不能直接应用SAR雷达方程的若干结论[912]。尤其是在空间目标ISAR成像中,轨道参数对成像质量的影响未有深入分析,大角度成像的方位相干积累增益能否弥补轨道升高带来的回波信噪比降低尚无明确结论。

      类比SAR雷达方程,本文结合目标轨道参数推导针对空间目标ISAR成像雷达方程的一般形式,进而分析轨道高度、目标雷达姿态角变化引起的方位相干积累增益变化及其对目标成像质量的影响,得到较为直观的理论来指导ISAR成像体制中发射功率、波形等参数的优化设计。仿真实验就不同轨道高度成像质量进行比较,探究方位向增益与轨道参数间的变化关系及其对成像质量的改善情况,实验结果验证了推导结论的准确性。

      本文结构如下:第2节结合天线理论推导ISAR系统成像雷达方程;第3.1节从地基雷达观测几何出发对空间目标角速度这一关键因素进行推导,并利用其计算结果得到空间轨道目标ISAR系统成像雷达方程;第3.2节分析第3.1节中观测模型引起计算误差;第3.3节对得到的空间轨道目标ISAR系统成像雷达方程作出一些定性结论;第4节结合目标轨道信息,利用空间轨道目标ISAR系统成像雷达方程指导定分辨率成像仿真,并对轨道参数引起的成像质量变化分析,验证了推导公式和结论的正确性。

    • 结合天线理论,对收发共天线的雷达系统,单次回波信号功率与雷达参数和作用距离的关系可由传统雷达方程表示:

      $S = \frac{{P{G\, ^2}{\lambda ^2}\sigma {T\,_{\rm{i}}}}}{{{{\left( {4{{π}} } \right)}^3}{r^4}\eta }}$

      其中,P为雷达发射机功率,G为天线增益, $\lambda $为信号波长, $\sigma $为目标的雷达截面积,Ti为相干积累时间,r为雷达作用距离, $\eta $为系统损耗。

      区别于一般雷达系统,ISAR成像体制具有2维高分辨率,而空间目标尺寸较大。因此,本文针对成像分辨单元信噪比进行成像质量分析,定义目标分辨单元的等效雷达截面积 $\sigma $[13]

      $\sigma = {\sigma _0}{\rho _{\rm{r}}}{\rho _{\rm{a}}}$

      其中, ${\sigma _0}$为目标归一化后向散射系数, ${\rho _{\rm{r}}}$和 ${\rho _{\rm{a}}}$分别为目标ISAR成像纵向和横向分辨率且两者大小相当。一般来说,ISAR成像系统的纵向距离分辨率与斜距r无关而仅和发射信号频率带宽相关[9]

      ${\rho _{\rm{r}}} = \frac{\rm c}{{2B}}$

      其中,c为光速,B为信号带宽。

      横向高分辨是通过在相干测量时间内对多帧回波进行多普勒分析获得,与目标雷达视线的相干积累转角直接相关。空间目标在观测过程中姿态平稳,可采用2维转台模型描述其与雷达间的相对运动。假定雷达视线转速为w,相干测量时间为Ti,则相干积累转角 $\Delta \theta = w{T\, _{\rm{i}}}$,横距分辨率 ${\rho _{\rm{a}}}$可计算为:

      ${\rho _{\rm{a}}} = \frac{\lambda }{{2\Delta \theta }}$

      将式(3)、式(4)带入重写式(1),考虑方位相干积累的ISAR成像目标某分辨单元对应的接收功率可表达为:

      $S = \frac{{P{G \, ^2}{\lambda ^2}{\sigma _0}{\rho _{\rm r}}}}{{{{\left( {4{{π}} } \right)}^3}{r^4}\eta }} \cdot \frac{\lambda }{{2w}}$

      由式(5)可见,与SAR成像雷达方程类似,雷达视角(Light Of Sight, LOS)变化速度w,即下文所述的目标相对转台中心转角速度,将直接影响分辨单元的回波功率,决定分辨单元信噪比质量。下面将结合空间轨道目标轨道参数对以上公式进行扩展分析。

    • 本节将从轨道高度变化引起的空间目标转角速度变化出发对上节得到的空间目标ISAR成像方程进行完善。在分析转角速度变化过程中,还将对关键性的雷达斜距变化进行建模分析,得到定分辨率观测下的雷达方程。

      假定空间目标处于近圆轨道,考虑地球自转后,目标的相对平近角点可表示为[14,15]

      $n = \frac{{\sqrt \mu }}{{{{({{\mathop{R}\nolimits} _{\rm{e}}} + h)}^{3/2}}}}$

      其中,引力常数 $\mu = 3.986 \times {10^{14}}({{\rm m}^3}/{{\rm s}^2})$,Re为地球半径,h为目标轨道高度。

      图1(a)观测几何所示,单次脉冲周期Δtm后,目标从p点运动至p′点,其运行绝对距离可由几何计算得:

      $\mathop {pp'}\limits^{\frown} = ({{\mathop{ R}\nolimits} _{\rm{e}}} + h)n\Delta {t_{\rm{m}}} = \frac{{\sqrt \mu \Delta {t_{\rm{m}}}}}{{{{({{\mathop{R}\nolimits} _{\rm{e}}} + h)}^{1/2}}}}$

      图  1  空间观测几何模型及成像模型

      Figure 1.  Observation and imaging model of satellite targets

      在目标轨道平面内,目标对于观测站点的相对运动如图1(b)所示,雷达站点处于转台原点,运动起点p的雷达斜距为r,运动终点p′的雷达斜距为r′,原点与终点连线上有一斜距为rp0,由几何关系可知:

      $\!\!\! r' = r + \Delta r$

      $\mathop {p{p_0}}\limits^{\frown} = \omega \Delta {t_{\rm{m}}}r$

      考虑空间目标平近点角远小于地球自转、单帧成像时间在秒级的情况下,短弧段 $\mathop {pp'}\limits^{\frown} $可近似为线段pp′,短弧段 $\mathop {p{p_0}}\limits^{\frown} $亦可近似为线段pp0, Δpp0p′可近似为直角三角形满足勾股定理:

      ${\left( {p{p_0}} \right)^2} + {\left( {{p_0}p'} \right)^2} = {\left( {pp'} \right)^2}$

      带入式(7)、式(8)、式(9),可得:

      $w =\frac{{\sqrt { \displaystyle\frac{{\mu \Delta {t_{\rm{m}}}^2}}{{{{\mathop{R}\nolimits} _{\rm{e}}} + h}} - \Delta {r^2}} }}{{r\Delta {t_{\rm{m}}}}}$

      r2=0处泰勒展开,可以得到近似解:

      $w \approx \frac{1}{{r\Delta {t_{\rm{m}}}\sqrt { \displaystyle\frac{{\mu \Delta {t_{\rm{m}}}^2}}{{{{\mathop{R}\nolimits} _{\rm{e}}} + h}}} }}\left( {\frac{{\mu \Delta {t_{\rm{m}}}^2}}{{{{\mathop{R}\nolimits} _{\rm{e}}} + h}} - \frac{{\Delta {r^2}}}{2}} \right)$

      其中,关于雷达斜距变化Δr的求解问题可以简化为研究位于目标轨道平面外一点与轨道上的目标间距离变化关系,如图2所示。

      图  2  雷达斜距变化Δr求解模型

      Figure 2.  Calculation model of change of radar range

      以地心为坐标轴原点,目标所在轨道平面为xOy平面建立直角坐标系xyz,观测站点SxOy平面内投影为 $S'\left( {{x_0},{y_0},0} \right)$,目标运动至 $p\left( {\left( {{{{R}}_{\rm{e}}} + h} \right)} \right.$ $\left. \cdot {\cos \theta \left( t \right),\left( {{{{R}}_{\rm{e}}} + h} \right)\sin \theta \left( t \right),0} \right)$,其中,q(t) = q0+ $n\left( {t - {t_0}} \right)$用以表示目标在轨道上的瞬时位置, ${\theta _0}$为单帧成像中心时刻t0对应的位置参数。对于直角三角形DSSp,由勾股定理知:

      ${r^2} = {d^2} + {l^2}$

      其中,垂直距离d在目标运动过程中视为定值,而水平距离l可由下式计算:

      $\begin{aligned}{l^2} = & {\left( {\left( {{{{R}}_{\rm{e}}} + h} \right) \cos\theta \left( t \right) - {x_0}} \right)^2}\\ & + {\left( {\left( {{{{R}}_{\rm{e}}} + h} \right)\sin \theta \left( t \right) - {y_0}} \right)^2}\end{aligned}$

      将式(14)带入式(13)并对等式两边关于时间求导:

      $\begin{array}{*{20}{c}}\begin{aligned}rr' = & \left[ { - \left( {\left( {{{{R}}_{\rm{e}}} + h} \right) \cos\theta \left( t \right) - {x_0}} \right)\left( {{{{R}}_{\rm{e}}} + h} \right)\sin \theta \left( t \right)} \right.\\& +\left. {\left( {\left( {{{{R}}_{\rm{e}}} + h} \right)\sin \theta \left( t \right) - {y_0}} \right)\left( {{{{R}}_{\rm{e}}} + h} \right)\cos \theta \left( t \right)} \right]n\end{aligned}\\\quad\quad\! { = \left[ {\left( {{{{R}}_{\rm{e}}} + h} \right){x_0}\sin \theta \left( t \right) - \left( {{{{R}}_{\rm{e}}} + h} \right){y_0}\cos \theta \left( t \right)} \right]n}\\\!\!\!\!\!\!\!\! { = \left( {{{{R}}_{\rm{e}}} + h} \right)\sqrt {{x_0}^2 + {y_0}^2} \sin \left( {\theta \left( t \right) - \varphi } \right)n}\end{array}$

      短时Δtm内, $r \approx r' \Delta {t_{\rm{m}}}$,并将式(15)带入可得:

      $\Delta r = \frac{{\left( {{{\mathop{R}\nolimits} _{\rm{e}}} + h} \right)\sqrt {{x_0}^2 + {y_0}^2} \sin \left( {\theta \left( t \right) - \varphi } \right)n\Delta {t_{\rm{m}}}}}{r}$

      其中, $\varphi = \arctan \left( {\displaystyle\frac{{{y_0}}}{{{x_0}}}} \right)$。

      为进一步简化式(16),当以OS′作为x正半轴,即y0=0, $\varphi = 0$时,在雷达可视范围内 $\left| \theta \right| \le {\theta _0}$,式(16)可写为:

      $\Delta r = \frac{{\left( {{{\mathop{R}\nolimits} _{\rm{e}}} + h} \right){x_0}\sin \theta \left( t \right)n\Delta {t_{\rm{m}}}}}{r}$

      重写式(11)、式(12):

      $\!\!\!\!\!\!\! w = \sqrt {\frac{\mu }{{r\left( {{{\mathop{R}\nolimits} _{\rm{e}}} + h} \right)}} - {{\left( {\frac{{\left( {{{\mathop{R}\nolimits} _{\rm{e}}} + h} \right){x_0}\sin \theta \left( t \right)n}}{{{r^3}}}} \right)}^2}} $

      $w \! \approx \! \sqrt {\frac{\mu }{{r\left( {{{\mathop{R}\nolimits} _{\rm{e}}} + h} \right)}}} \!-\! \frac{1}{{2{r^3}\sqrt \mu }}{\left( {{{\mathop{R}\nolimits} _{\rm{e}}} + h} \right)^{5/2}{\left( {{x_0}\sin \theta \left( t \right)n} \right)^2}}$

      为保证方位向分辨率一定,总转角 $\Delta \theta $需固定,那么相干时间Ti将随目标高度变化:

      ${T\,_{\rm{i}}} = \frac{{\Delta \theta }}{w}$

      结合式(6)、式(18),重写空间目标ISAR成像雷达方程:

      $S = \frac{{P{G \, ^2}{\lambda ^3}{\sigma _0}{\rho _{\rm{r}}}}}{{2{{\left( {4{{π}} } \right)}^3}{r^3}\eta }} \cdot \sqrt {\frac{{{r^2}\left( {{{\mathop{R}\nolimits} _{\rm{e}}} + h} \right)}}{{\mu {r^2} - {{\left( {{{\mathop{R}\nolimits} _{\rm{e}}} + h} \right)}^3}{{\left( {{x_0}\sin \theta \left( t \right)n} \right)}^2}}}} $

    • 由于地球自转的影响,空间目标相对于雷达站点的运动轨迹并不是一个闭合的圆,如图3(b)所示。本文采用的近似模型将其轨道近似为圆如图3(a)所示,对于短时间观测任务来说,两者误差可控制在一定程度内。这里分别使用不同高度的轨道参数对LOS转角速度w进行求解并与实际值对比,选择库尔勒作为观测站点,结果如图4所示。其中,真实值为通过STK验证后的实际LOS转角速度,理论值为按式(18)计算得到的结果,近似值为按式(19)计算得到的结果。计算结果与实际值间的误差主要由模型简化引起,可由与轨道高度相关的函数补偿。补偿后的雷达作用距离方程可写成:

      $\begin{aligned}S & = \frac{{P{G \, ^2}{\lambda ^3}{\sigma _0}{\rho _{\rm{r}}}}}{{2{{\left( {4{π} } \right)}^3}{r^3}\eta }}\\& \cdot \!\! {\left(\! {\sqrt {\frac{\mu }{{{{{R}}_{\rm{e}}} \!+\! h}} \!-\! \frac{{{{\left( {\left( {{{{R}}_{\rm{e}}} + h} \right){x_0}\sin \theta \left( t \right)n} \right)}^2}}}{{{r^2}}}} \!-\! C\left( h \right)r} \right)^{ - 1}}\end{aligned}$

      其中,C(h)为补偿函数。需要说明的是,补偿函数与目标轨道参数和观测点坐标均有关,且对于本文关于空间目标ISAR成像雷达方程的影响有限,故未进一步讨论。

      图  3  Δr求解近似模型与实际模型

      Figure 3.  Approximate and actual model of Δr

      图  4  不同轨道高度下目标转角速度计算结果

      Figure 4.  Calculation of target′s rotate speed at different heights

    • 与条带SAR雷达方程类似,式(20)中ISAR系统回波信号的接收功率与横向分辨率 ${\rho _{\rm{a}}}$无关,与LOS转角速度w成反比,而w随轨道目标轨道高度升高而降低。轨道高度升高引起转角速度减小进而导致成像所需相干时间增加。这一变化可部分抵消因目标斜距增大引起的回波能量分散效果,也就是说回波信号接收功率因相干增益不再随作用距离增大呈4次方下降,其下降速度应小于作用距离的4次方,具体数值还与目标与雷达间相对位置有关,一般应介于3次方与4次方之间。采用大角度的ISAR成像处理可弥补轨道升高带来的回波接收功率降低,较普通雷达体制有明显的距离优势。

    • 为验证空间目标ISAR成像雷达方程中,回波信号接收功率下降速度的结论。仿真实验将在成像分辨率固定的情况下,仅改变目标轨道高度引起作用距离变化,对相近姿态下的空间目标进行成像观测,研究得到的RD图像中信号功率的变化以及图像质量的变化。实验中ISAR系统主要参数如表1,目标轨道主要参数倾角为42.8°,升交点赤经(Right Ascension of Ascending Node, RAAN)为34.7°E,观测点分别选取库尔勒、北京、西安,其经纬度信息如表2

      参数 数值
      载频 16.7 GHz
      带宽 1 GHz
      方位向分辨率 0.18 m
      距离向分辨率 0.15 m
      脉冲重复频率 200 Hz

      表 1  实验ISAR系统主要参数

      Table 1.  Main parameters of ISAR system

      地点 经纬度
      库尔勒 41.5°N, 86.8°E
      北京 39.9°N, 116.4°E
      西安 31.1°N, 108.4°E

      表 2  实验地基ISAR观测站位置

      Table 2.  Position parameters of radar sites

    • 对于式(21),为简化计算,本实验中选取 ${\theta _0} = 0$,也就是目标均处于轨道上与观测点最近位置附近,其转角速度w达到该轨道上的最大值,其回波接收功率为同一轨道最小值。

      目标轨道高度变化将引起目标轨道半径的变化,图5为归一化的回波信号接收功率随目标轨道半径的变化曲线,可以看出对于轨道半径在7400 km以上(作用距离在1400 km以上)的目标,其回波信号接收功率随斜距下降速度介于斜距变化的3次方与4次方之间,具体影响因子与目标与观测点相对位置、观测弧段均有关。在雷达位于西安的观测过程中,低轨道观测甚至出现相干增益超过斜距下降3次方影响的现象,这是由于低轨目标其相干增益与其轨道高度直接相关,而斜距变化与轨道高度变化并不是严格的线性相关,也就说低轨观测中,目标与观测点相对位置以及观测弧段的变化也将较大程度影响回波信号的接收功率。

      图  5  归一化回波信号接收功率随目标轨道半径变化曲线

      Figure 5.  Normalization curve of received power changing in different orbit radii

    • 结合ISAR体制下目标所具有的孤立散射特点,参考点目标成像质量评价中的积分旁瓣比[16](Integrated SideLobe Ratio, ISLR),本文计算各图像中心单元的目标与背景噪声像素能量比(Target Noise Ratio, TNR)来反映成像质量。理论上来说,孤立散射点成像后对应像素单元与背景噪声的能量比应与脉冲回波信噪比、脉压长度、脉冲积累数等因素均相关。但为直观反映不同轨高下脉冲积累数变化带来的图像信噪比增益变化,本文将所有单脉冲回波信噪比统一设置为10 dB,仅改变脉冲积累数进行实验。

      ${\rm{TNR = }}\frac{{{{{E}}_{\rm{t}}}}}{{{{{E}}_{\rm{b}}}}}$

      其中,Et为目标像素能量积分,Eb为背景像素能量积分。

      实验选取西安站观测结果进行成像质量分析,其中方位维、距离维无单位,代表像素点位置,图6为某空间目标在791 km轨道某处成像结果,目标中心单元的距离维、方位维剖面如图6所示。选取相近姿态下,目标在3个不同轨道高度的成像结果作为对比,如图7所示,其量化质量评价结果如表3所示。

      表3可以看出,随轨道高度升高脉冲积累数增加,相干处理后图像TNR也相应增大,这与3.3节中轨道高度对大角度ISAR成像体制影响的结论一致;其脉冲积累数与TNR之比直观反映相干积累对图像质量提升的作用,在脉冲回波信噪比、脉压长度、等因素相同的情况下,可近似为一定值,但可以预见的是在实际中将随着高度升高而增大。

      图  6  目标图像距离、方位剖面图

      Figure 6.  Range and azimuth profiles of target image

      图  7  相近观测视角下,不同轨道高度成像结果对比图

      Figure 7.  Comparison result of target imaging at different heights with similar LOS parameters

      轨道高度(km) 成像时间(s) 脉冲积累数 TNR 脉冲积累数/TNR
      791 9.94 1988 6.82 291.58
      1200 12.40 2484 8.61 288.58
      1800 16.82 3364 11.45 293.80

      表 3  不同轨道高度图像质量评价

      Table 3.  Comparison result of imaging quality at different heights

    • 本节实验应用推导的雷达成像方程从回波功率、成像质量两方面进行定分辨成像分析,可总结以下结论。(1)通过有效结合轨道信息,空间目标ISAR成像处理应采用更接近于合作(或半合作)目标的成像处理方式,进一步根据本文方法估计方位分辨性能可在保证横向分辨率的基础上有效指导成像角域优化选择。(2)总体而言,空间目标ISAR观测的回波信噪比受轨道高度升高而下降,另一方面较大的相干积累角ISAR成像,方位相干积累增益可部分补偿目标轨道高度增加引起的信噪比损失,也就是文中所述的回波接收功率下降量级小于斜距4次方,但大于SAR系统中斜距3次方的关系。(3)目标轨道参数、观测几何模型可有效指导空间轨道目标的成像工作功率、波形参数设计,实现ISAR成像信噪比预估计,同时可利用相干积累角的计算进行成像时间段的优化选择,满足高分辨成像任务。

    • 本文从基本雷达方程出发,推导空间轨道目标ISAR成像雷达方程的一般形式,分析空间目标轨道参数对目标转角速度以及成像相干积累增益的影响,定性分析了采用相干体制下的大转角ISAR雷达系统进行空间观测的优势,提出较为简便的成像信噪比估计方法。仿真实验验证ISAR成像的方位相干增益可部分弥补目标轨道高度增加引起的ISAR成像质量下降,为空间目标ISAR成像体制和信号处理设计、成像信噪比估计、成像时间段优化选择提供了理论基础。

参考文献 (16)

目录

    /

    返回文章
    返回