一种SAR回波历程快速实时计算方法

摘要：      针对合成孔径雷达（Synthetic Aperture Radar， SAR）目标回波模拟计算量过大的问题， 采用两次多阶泰勒级数展开对回波系统函数中的距离历程进行计算， 并结合数据外推的方法完成距离历程的快速实时计算。 该方法使用多次乘法和加法代替计算量较大的开方运算， 可保证在一定计算精度的前提下提高距离计算的实时性， 同时结合距离历程模型的规律性， 采用数据外推代替重复的计算过程， 进一步降低计算量， 提高实时性。 通过仿真对比不同阶数泰勒级数展开在同一距离处的计算误差， 得到了4阶为该仿真条件下的最优泰勒展开阶数，  最后分别对9个分布密集的点目标和真实场景目标进行回波模拟和成像。 点目标仿真结果表明， 该方法在不影响回波信号性能的前提下可以使生成回波信号的时间减少53%， 速度提升2倍； 真实场景仿真结果表明， 该方法可以应用于大场景的回波模拟。

关键词：     回波模拟； 回波系统函数； 实时计算； 距离历程； 数据外推； 点目标； SAR

中图分类号：     TJ760; TN95    文章编号：     1673-5048（2023）04-0098-08

文献标识码：    A    DOI： 10.12132/ISSN.1673-5048.2022.0269

0引言

合成孔径雷达（Synthetic Aperture Radar， SAR）是一种能实现高分辨、 远距离、 广范围的成像雷达［1］， 具有很高的应用价值， 是目前雷达的主要成像方式。

随着真实的雷达图像分辨率不断提高， 对模拟图像分辨率的要求也越来越高［2］， 研究人员普遍开始采用基于电磁计算的雷达成像仿真技术［3-4］， 但这种计算量是相当大的， 一般需要通过数字信号处理器（Digital Signal Processor， DSP）、 现场可编程逻辑门阵列（Field Programmable Gate Array， FPGA）等硬件加速来保证实时性［5-7］。 此外， 图形处理器（Graphics Processing Unit， GPU）由于多核多线程的特性被广泛应用于SAR回波模拟中。 早在2009年， 意大利萨勒诺大学的Clemente等人利用计算统一设备架构（Compute Unified Device Architecture， CUDA）编程平台， 在Nvidia Tesla C1060显卡上实现了SAR回波模拟［8］。 同年， 文献［9］利用中央处理器（Central Processing Unit， CPU）和GPU混合编程实现了SAR的面目标回波模拟。 目前， 国内大部分研究已经在GPU上实现了回波模拟算法的并行化， 但是其采用的回波模拟算法大多是传统的回波模拟算法， 在仿真大场景时依然存在效率较为低下的问题， 有进一步的加速空间。 为了实现实时计算， 文献［10］提出了等距离环的概念， 对成像场景进行等距分割， 并使用多个FPGA并行计算； 文献［11］通过使用宽带数字射频存储器（Digital Radio Frequency Memory， DRFM）技术和并行计算技术实现了较小延迟的SAR回波计算和生成。 即使在硬件平台上进行回波模拟， 仍存在模拟点数过多和系统函数计算量过大的问题， 当前大多数研究围绕在如何减少模拟点数， 针对SAR合成孔径的特性， 越来越多的研究采用子孔径处理［12-17］进行回波模拟， 通过在回波信号方位向划分子孔径， 分段处理， 一定程度上解决了模拟点数过多的问题， 而对系统函数计算量的研究数量不多。

距离历程是回波系统函数中的一个重要参数， 针对回波系统函数计算量过大的问题， 本文提出使用多阶泰勒级数展开对距离历程进行计算， 在保证一定精度的前提下提高距离计算的实时性， 通过对不同阶数泰勒级数展开产生的计算误差进行仿真对比， 找出最优泰勒级数展开阶数， 并从理论上分析了距离历程模型的规律性， 提出存储部分数据来代替重复计算， 进一步减少计算量。

1SAR目标回波信号模型

1.1SAR空间几何模型

SAR的回波信号是雷达平台在各个时刻发射波束覆盖范围内每个点目标回波信号的叠加。 图1是条带式SAR的空间几何模型， 选取直角坐标系为参考坐标系， rOy为地平面， 假设雷达平台高度为H， 沿方位向y正方向以速度V匀速飞行， 机载平台向航向正侧方向地面发射电磁波， 俯仰角度为φ， t时刻平台坐标为（0， y（t）， H）， rOy平面上波束中心坐标为（r0， y（t）， 0）， 波束宽度为θ， 波束照射区域为阴影部分， 测绘带宽度为W。

1.2目标回波信号模型

假设雷达平台发射的脉冲串为

式中： rect［·］为矩形函数； Tp为发射信号的脉冲宽度； fc为载频； k为调频斜率。

假设（ri， yi，0）为探测区域内任一散射点， 由图1可得， 该点与雷达平台的瞬时斜距为

式中： t为慢时间； ta为快时间； A为接收信号的幅度因子， 与目标的后向散射系数和天线方向图有关； c为光速。

解调去掉载波后， 点目标的回波信号可以写为

式中： 第一个指数项为多普勒分量， 是回波信号的方位向相位； 第二个指数项为基带LFM的延迟分量。

实际上， 真实场景的目标是由地面大量散射点组成， 对于真实场景的模拟中， 把每个散射点当成一个单点目标考虑， 整个场景的目标回波可以看作所有散射点回波的叠加， 则真实场景回波信号模型可表示为

式中： L为场景中散射点的总个数； Ak为第k个散射点回波的幅度因子； Rk（t）为第k个散射点在t时刻的瞬时斜距。

2SAR目标回波系统函数的实时计算

2.1SAR目标回波系统函数模型

SAR目标回波是雷达照射波束内全部散射点回波的叠加， 可以看作是雷达发射信号经过一个系统后的输出。 因而SAR目标回波可表示为雷达发射脉冲s（t）与目标回波系统函数h（t）的卷积， 即

sr（t）=s（t）h（t） （7）

式中： 回波系统函数h（t）包含波束照射范围内所有散射点回波的延迟、 幅度以及方位相位等信息， 可表示为

由式（8）可知， 目标回波系统函数是一系列冲激函数的叠加， 这些冲激函数具有不同的延时、 幅度和相位， 其作用就是对雷达发射信号进行延迟、 幅度和相位的调制。 在距离时域脉冲相干法中， 将根据目标回波系统函数的采样周期对散射点的延迟进行近似处理。 此时在某个采样时刻的回波系统函数将是落入该时刻所代表距离门的所有散射点信息的叠加， 即

式中： M为落入该距离门散射点的个数； Ak， Rk和tk分别为该距离门内第k个散射点回波的幅度、 距离和延迟； Ts为系统函数的采样周期。

由于雷达波束在地面照射范围内的散射点数目非常多， 因而在SAR目标回波仿真中， 回波系统函数的计算量将非常大。

2.2高实时性距离历程的计算方法

2.2.1基于多阶泰勒级数展开的距离历程模型

回波系统函数主要是计算每个散射点的距离， 根据距离可以确定回波相位以及该散射点所处的距离门， 最后对处于同一个距离门的目标信息进行累加。 因此， 距离历程是回波系统函数的一个重要参量。

以等距离环和方位向分辨单元对目标场景进行网格化处理， 场景中散射点的瞬时斜距模型， 如图2所示。

设SAR雷达高度为H0， t时刻波束中心的坐标为（x0， y0（t）， 0）， 则雷达发射波束覆盖区域内某个散射点（i， j）的瞬时斜距为

式中： Δx和Δy分别为地面上相邻散射点在距离向和方位向的分辨间隔； i和j分别为目标区域中散射点位于第i个距离环（距离向分辨单元）和位于第j个方位向分辨单元， 即i和j可表示散射点的行列号。

由式（10）可见， 传统的距离计算方法存在两个缺点： （1）使用开方进行计算， 开方运算具有很大的计算量， 且不便于FPGA等硬件平台实现， 不利于成像目标回波的实时计算； （2）通常H0， x0两个参数的数值会很大， 会导致每次计算的距离数值过大。

为了便于目标回波实时计算， 本文提出采用多阶泰勒级数展开对式（10）进行计算， 使用多个乘法和加法来代替开方计算。 考虑到某一行目标数量通常比某一列目标数量多， 首先对于中心列的散射点进行泰勒级数展开。 对于某一列目标， 其行号都相同， 视（i·Δy）为一个整体进行多阶泰勒级数展开， 则式（10）可表示为

R0（t）为雷达平台到波束中心点的斜距， 通过对距离历程模型进行多阶泰勒级数展开， 这个数值最大的参数被分离， 可以作为一个常数来处理， 只需要计算泰勒展开项的各次幂和各次幂的系数即可， 避免了大数值参数的重复计算。

此外， 由于雷达载波信号的波长很短， 因而散射点距离的计算必须保持很高的精度， 以保证回波信号相位的计算精度， 这便为式（11）～（12）中泰勒级数展开所取的阶数提出了一定的要求， 仿真时需要根据具体的参数选取最优泰勒级数展开阶数。

2.2.2距离历程模型的规律性

假设SAR目标模拟器数据参数更新周期为T， T=Δy/v， 表示每当雷达平台移动一个方位分辨单元的距离， 数据更新一次， 则t时刻与经历多个T后某等距离环中心处各等效散射点的距离变化， 如图3所示。

对于任意t+m·T时刻， （i+m，j）处等效散射点的距离如式（13）所示， 其与t时刻（i，j）处等效散射点的距离完全一致：

由于i， j， t， m均是随机取值， 故式（13）具有普遍性， 可以得到

R（i， j， t）=R（i+m， j， t+m·T）（14）

式中： 0≤t<T； 0≤m≤N。

利用R（i，j，t）的规律性， 采用数据外推的方法对距离历程进行实时计算： 首先将t∈［0，T）内所有发射波束覆盖范围内各散射点的距离数据按脉冲编号存储起来， 在后续周期中按照脉冲在该周期中的编号直接调用所需的数据， 以此来代替每个周期都对覆盖范围内各散射点进行距离计算， 可以减少大量的重复操作。 此方法以一定的存储器资源为代价， 大大降低了运算量， 在存储器资源足够的情况下， 对于距离历程以及回波系统函数的实时计算有很大的帮助。

2.3系统函数的实时计算

散射点回波方位向相位表达式为

利用FPGA计算回波相位时， 三角函数的计算需要利用查表的方法来实现， 所以相位计算将转换为计算相位查询表的地址。

文中使用两次多阶泰勒级数展开计算各散射点距离， 两次展开的实现方式分别如下：

（1） 首先计算中心行目标距离及相关信息。 根据泰勒级数展开公式， 可得各目标相对于波束中心点的相对距离表达式为

R1=R0+k1·C1+k2·C2－k3·C3+…（16）

式中： C1， C2， C3， …为泰勒展开项的各次幂； k1， k2， k3， …为展开式中各次幂的系数的分子部分。

各参数表达式为

由于x0， y0， H0， Δy均为定值， 所以式（17）中的所有参数均为定值， 可以直接从上位机中获取。