1 引言

双基星载合成孔径雷达(Synthetic Aperture Radar, SAR)利用信号收发平台的分置，能够同时获取不同视角的观测数据，在测绘、干涉测量、地面目标识别、自然灾害监测等领域^[1,2]具有重要的应用价值。以Tandem-L为代表的新一代双基星载SAR系统应用多通道、数字波束形成(Digital Beam Forming, DBF)等技术，实现方位向高分辨率和距离向宽测绘带SAR(High Resolution Wide Swath SAR, HRWS-SAR)成像，系统的成像能力得到显著改善^[3,4]。

HRWS-SAR系统在方位向采用多通道采样降低系统的脉冲重复频率(Pulse Recurrence Frequency, PRF)，在不降低方位分辨率的前提下实现系统的宽测绘带成像。HRWS-SAR系统进行成像处理时，若方位向信号为均匀采样可直接采用传统SAR成像处理方法，而受载星平台轨迹约束及应用场景的限制，方位向的非均匀采样更为普遍^[5]，因而对方位向信号进行重构以获得其均匀采样信号或多普勒谱是HRWS-SAR系统成像处理的一项关键技术。单基星载HRWS-SAR系统的方位向信号重构得到了广泛而深入的研究^[6–17]。重排算法^[6]和插值算法^[7]是两种典型的时域重构算法。重排算法依据接收信号的方位向位置将各通道信号重新排列，只能在特定的PRF得到方位向信号的均匀采样。插值算法则依据文献[7]所提出的周期性非均匀采样信号重构公式通过时域插值得到方位向信号的均匀采样，该算法运算复杂度高且精度依赖于插值核的长度。Krieger等依据广义采样定理提出矩阵求逆算法^[8,9]，该算法通过线性方程求解从混叠的多通道信号中重构出无模糊的多普勒谱。该算法不需要协方差矩阵等先验信息且易于实现，但在重叠采样时该算法无法进行信号重构，接近重叠采样时重构性能也较差。文献[10–13]则采用不同的方法对矩阵求逆算法进行改进，实现重叠采样附近的高性能信号重构。文献[14,15]基于统一的信号模型，选取不同代价函数进行优化，提出了多种自适应波束形成(DBF)算法：正交投影算法、信号最大化算法、最大化信号模糊噪声比算法、最小均方误差算法等。与矩阵求逆算法相比，DBF类算法运算复杂度较高，并假定各通道的噪声为高斯白噪声且相互独立。文献[16]则对方位向信号重构算法的性能进行了仿真对比分析。文献[17]采用NUFFT直接重构方位向非均匀采样信号的多普勒频谱。

上述方位向信号重构方法理论上均能推广到双基星载HRWS-SAR系统。考虑到矩阵求逆算法易于实现且研究较为广泛，本文研究该算法在双基HRWS-SAR系统中的实现。本文首先将方位照射时间内时变的发射接收距离比近似为一个常数，推导了双基星载HRWS-SAR系统与某个单基系统在方位向多信道间传递函数上的等效性，从而构建了双基星载HRWS-SAR系统的方位向信号模型。而后，提出了适用于一般双基构型星载HRWS-SAR系统方位向信号重构的矩阵求逆算法，并给出了信噪比缩放因子及方位模糊比这两个重构性能指标的计算公式。最后，通过对几种典型双基构型的星载HRWS-SAR系统进行方位向信号重构仿真，验证了矩阵求逆算法在一般双基构型星载HRWS-SAR系统中的适用性。

2 双基星载HRWS-SAR系统的方位向信号 模型

为保持双基星载HRWS-SAR系统对地面的持续观测，收发平台应置于同一轨道(顺飞模式)或高度相同的平行轨道(平飞模式)。图1给出了一般构型双基星载HRWS-SAR系统的信号收发几何。发射天线和接收天线的最短距离分别为r_T0, r_R0，接收天线共有M个通道，RX_equ为方位向信号重构后的等效接收通道，通道i到等效接收通道RX_equ的方位向距离为 $\Delta {x_i}$ , t_fd为发射天线和接收天线的零多普勒时间差。

不失一般性，假定点目标P在t=0时刻位于接收天线的零多普勒面，容易得到t时刻点P在接收通道i的收发距离和 $R\left( t \right)$ 计算公式为：

$\begin{aligned}R\left( t \right) =& \sqrt {r_{{\rm{T}}0}^2 + {{\left( {vt - v{t_{{\rm{fd}}}}} \right)}^2}} \\& + \sqrt {r_{{\rm{R}}0}^2 + {{\left( {vt - \Delta {x_i}} \right)}^2}} \end{aligned}$

(1)

忽略天线方向图的影响，接收通道i的冲激响应函数可表示为：

$\begin{aligned}{h_i}\left( t \right) = & \exp \left( { - {\rm j}\frac{{2{{π}} }}{\lambda }R\left( t \right)} \right) \\ = & \exp \Biggr( { - {\rm j}\frac{{2{{π}} }}{\lambda }\biggr( {\sqrt {r_{{\rm{T}}0}^2 + {{\left( {vt - v{t_{{\rm{fd}}}}} \right)}^2}}}} \\ & { {\left.+ {\sqrt {r_{{\rm{R}}0}^2 + {{\left( {vt - \Delta {x_i}} \right)}^2}} } \right)} \Biggr)}\end{aligned}$

(2)

2.1 单基星载HRWS-SAR系统的方位向信号模型

单基星载HRWS-SAR系统可视为r_T0=r_R0=r₀, t_fd=0，发射天线为等效接收通道的特殊情况。由式(2)可得此时接收通道i的冲激响应为：

$\begin{aligned}{\tilde h_i}\left( t \right) = & \exp \Biggr( { - {\rm j}\frac{{2{{π}} }}{\lambda }\biggr( {\sqrt {r_0^2 + {{\left( {vt} \right)}^2}}}} \\& + {{\sqrt {r_0^2 + {{\left( {vt - \Delta {x_i}} \right)}^2}} } \biggr)} \Biggr)\end{aligned}$

(3)

令此时单基单通道SAR的冲激响应函数为 ${\tilde h_{{\rm{mono}}}}\left( t \right)$ , ${\tilde h_{{\rm{mono}}}}\left( t \right)$ 可由式(4)表示。

${\tilde h_{{\rm{mono}}}}\left( t \right) = \exp \left( { - {\rm j}\frac{{4{{π}} }}{\lambda }\sqrt {r_0^2 + {{\left( {vt} \right)}^2}} } \right)$

(4)

文献[9]对 ${\tilde h_{{\rm{mono}}}}\left( t \right)$ 及 ${\tilde h_i}\left( t \right)$ 的泰勒级数展开进行2阶相位近似处理，可得

$\left\{ \begin{aligned}&{{\tilde h}_{{\rm{mono}}}}\left( t \right) \approx \exp \left( { - {\rm j}\frac{{4{{π}} }}{\lambda }{r_0}} \right) \cdot \exp \left( { - {\rm j}\frac{{2{{π}} {v^2}{t^2}}}{{\lambda {r_0}}}} \right)\\&{{\tilde h}_i}\left( t \right) \approx \exp \left( { - {\rm j}\frac{{4{{π}} }}{\lambda }{r_0}} \right) \cdot \exp \left( { - {\rm j}\frac{{{{π}} \Delta x_i^2}}{{2\lambda {r_0}}}} \right) \\& \quad\quad\quad\ \cdot \exp \left( { - {\rm j}\frac{{2{{π}} {v^2}{{\left( {t - \displaystyle\frac{{\Delta {x_i}}}{{2v}}} \right)}^2}}}{{\lambda {r_0}}}} \right)\end{aligned} \right.\quad$

(5)

因而 ${\tilde h_{{\rm{mono}}}}\left( t \right)$ 与 ${\tilde h_i}\left( t \right)$ 间存在式(6)所示的关系，其中 $* $ 为卷积运算。

${\tilde h_i}\left( t \right) \!=\! \exp \left( { \!- {\rm j}\!\frac{{{{π}} \Delta x_i^2}}{{2\lambda {r_0}}}} \right) \!\cdot\! \delta \left( {t \!-\! \frac{{\Delta {x_i}}}{{2v}}} \right)\! *\! {\tilde h_{{\rm{mono}}}}\left( t \right)$

(6)

将式(6)进行Fourier变换，可得单基单通道SAR与接收通道i的系统函数关系：

$\begin{aligned}{\tilde H_i}\left( f \right) = & \exp \left( { - {\rm j}\frac{{{{π}} \Delta x_i^2}}{{2\lambda {r_0}}}} \right) \cdot \exp \left( { - {\rm j}2{{π}} \frac{{\Delta {x_i}}}{{2v}}f} \right) \\&\cdot {\tilde H_{{\rm{mono}}}}\left( f \right)\end{aligned}$

(7)

单基单通道SAR到接收通道i的传递函数 ${\tilde G_i}\left( f \right)$ 为：

${\tilde G_i}\left( f \right) = \exp \left( { - {\rm j}\frac{{{{π}} \Delta x_i^2}}{{2\lambda {r_0}}}} \right) \cdot \exp \left( { - {\rm j}2{{π}} \frac{{\Delta {x_i}}}{{2v}}f} \right)$

(8)

基于上述分析，单基星载HRWS-SAR系统的方位向信号生成模型可由图2进行描述。方位向各通道信号在时域为单基单通道SAR信号经对应的相位偏移和时间延迟的结果，而在多普勒域为单基单通道SAR信号与对应的常数相位因子和线性相位因子相乘后的结果。PRF小于多普勒带宽时，方位向各通道信号为时域上的欠采样，将导致多普勒域的频谱混叠。

2.2 双基星载HRWS-SAR系统的方位向信号模型

一般双基构型的星载HRWS-SAR系统中r_T0≠r_R0, t_fd≠0，不能通过对h_mono及h_i的泰勒级数展开进行2阶相位近似处理得到与式(6)、式(7)相似的表达式。本节将对式(2)进一步处理，推导一般双基构型星载HRWS-SAR系统中 ${h_{{\rm{mono}}}}\left( t \right)$ 及 ${h_i}\left( t \right)$ 的关系，进而构建其方位向信号生成模型。

引入点目标P的收发距离比函数 $C\left( t \right) = $ $\displaystyle\frac{{\sqrt {r_{{\rm{T}}0}^2 + {{\left( {vt - v{t_{{\rm{fd}}}}} \right)}^2}} }}{{\sqrt {r_{{\rm{R}}0}^2 + {{\left( {vt} \right)}^2}} }}$ ，式(2)等价于：

$\begin{aligned}{h_i}\left( t \right) = & \exp \Biggr( { - {\rm j}\frac{{2{{π}} }}{\lambda }\biggr( {C\left( t \right)\sqrt {r_{{\rm{R}}0}^2 + {{\left( {vt} \right)}^2}}}} \\& + {{\sqrt {r_{{\rm{R}}0}^2 + {{\left( {vt - \Delta {x_i}} \right)}^2}} } \biggr)} \Biggr)\end{aligned}$

(9)

等效接收通道RX_equ的冲激响应函数为：

${{h}_{\text{mono}}}\left( t \right)=\exp (-\text{j}\frac{2\pi }{\lambda }\left( \left( C\left( t \right)+1 \right)\sqrt{r_{\text{R}0}^{2}+{{\left( vt \right)}^{2}}} \right))$

(10)

对式(9)、式(10)的泰勒级数展开进行2阶相位近似处理，可得：

$\left\{ \begin{aligned}&{h_i}\left( t \right) \approx \exp \left( { - {\rm j}\frac{{2{{π}} }}{\lambda }\left( {C\left( t \right) + 1} \right){r_{{\rm{R}}0}}} \right) \\& \quad\quad\quad\ \cdot \exp \left( { - {\rm j}\frac{{{π}} }{{\lambda {r_{{\rm{R}}0}}}} \cdot \frac{{C\left( t \right)\Delta x_i^2}}{{C\left( t \right) + 1}}} \right) \\& \quad\quad\quad\ \cdot \exp \Biggr( { - {\rm j}\frac{{{π}} }{{\lambda {r_{{\rm{R}}0}}}} \cdot \Big( {C^{} \left( t \right) + 1} \Big) }\\& \quad\quad\quad\ {\cdot {v^2}{{\left( {t - \frac{{\Delta {x_i}}}{{\left( {C\left( t \right) + 1} \right)v}}} \right)}^2}} \Biggr)\\& {h_{{\rm{mono}}}}\left( t \right) \approx \exp \left( { - {\rm j}\frac{{2{{π}} }}{\lambda }\left( {C\left( t \right) + 1} \right){r_{{\rm{R}}0}}} \right) \\& \quad\quad\quad\quad\ \cdot \exp \left( { - {\rm j}\frac{{{π}} }{{\lambda {r_{{\rm{R}}0}}}} \cdot \left( {C\left( t \right) + 1} \right) \cdot {v^2}{t^2}} \right)\end{aligned} \right.$

(11)

在目标照射时间T_a内， $C\left( t \right)$ 的变化极小，用波束中心穿越时刻的取值C₀予以近似。容易得到接收通道i与等效接收通道RX_equ的冲激响应关系为：

$\begin{aligned}{h_i}\left( t \right) = & \exp \left( { - {\rm j}\frac{{{π}} }{{\lambda {r_{{\rm{R}}0}}}} \cdot \frac{{{C_0}\Delta x_i^2}}{{{C_0} + 1}}} \right) \\& \ \cdot \delta \left( {t - \frac{{\Delta {x_i}}}{{\left( {{C_0} + 1} \right)v}}} \right) * {h_{{\rm{mono}}}}\left( t \right)\end{aligned}$

(12)

式(12)进行Fourier变换，可得等效接收通道RX_equ与接收通道i的系统函数关系：

$\begin{aligned}{H_i}\left( f \right) = & \exp \left( { - {\rm j}\frac{{{π}} }{{\lambda {r_{{\rm{R}}0}}}} \cdot \frac{{{C_0}\Delta x_i^2}}{{{C_0} + 1}}} \right) \\& \cdot \exp \left( { - {\rm j}2{{π}} \frac{{\Delta {x_i}}}{{\left( {{C_0} + 1} \right)v}}f} \right) \cdot {H_{{\rm{mono}}}}\left( f \right)\end{aligned}$

(13)

等效接收通道RX_equ到接收通道i的传递函数为：

$\begin{aligned}{G_i}\left( f \right) = & \exp \left( { - {\rm j}\frac{{{π}} }{{\lambda {r_{{\rm{R}}0}}}} \cdot \frac{{{C_0}\Delta x_i^2}}{{{C_0} + 1}}} \right) \\& \cdot \exp \left( { - {\rm j}2{{π}} \frac{{\Delta {x_i}}}{{\left( {{C_0} + 1} \right)v}}f} \right)\end{aligned}$

(14)

考察式(8)、式(14)，传递函数 ${G_i}\left( f \right)$ 与 ${\tilde G_i}\left( f \right)$ 可视为f及 $\Delta {x_i}$ 的二元函数，并存在如下关系：

${G_i}\left( {f,\Delta {x_i}} \right) = {\tilde G_i}\left( {f,\frac{2}{{\left( {{C_0} + 1} \right)}}\Delta {x_i}} \right)$

(15)

式(15)表明双基星载HRWS-SAR系统的方位向信号的生成模型可等效为一单基系统模型，只需将图2中的 $\Delta {x_i}$ 替换为 $2/{{\left( {{C_0} + 1} \right)}}\Delta {x_i}$ 即可予以描述，C₀=1时，双基系统退化为单基系统。

2.3 C(t)近似为C₀的合理性

式(12)、式(13)、式(14)是在 $C\left( t \right)$ 近似为C₀的基础上推导的结果，显然存在近似误差。上述近似误差进行理论上的严格定量分析非常复杂，然而可以从上述表达式中得到简单的定性结果。以式(14)为例，由于在双基星载HRWS-SAR系统中 $\Delta x_i^2 \ll {r_{{\rm{R}}0}}$ ，从而 $\exp \left( { - {\rm j}\displaystyle\frac{{{π}} }{{\lambda {r_{{\rm{R}}0}}}} \cdot \displaystyle\frac{{C\left( t \right)\Delta x_i^2}}{{C\left( t \right) + 1}}} \right)$ 可视为常数1。指数项 $\exp \left( { - {\rm j}2{{π}} \displaystyle\frac{{\Delta {x_i}}}{{\left( {{C_0} + 1} \right)v}}f} \right)$ 的相位近似误差约为 $\Delta \varphi \approx 2{{π}} \left| {\displaystyle\frac{{C\left( t \right) - {C_0}}}{{{C_0} + 1}}} \right| \cdot \left| {\displaystyle\frac{{\Delta {x_i}}}{{\left( {C\left( t \right) + 1} \right)v}}f} \right|$ ，在多普勒带宽B_fd范围内，令接收通道的长度为L_R, $\mathop {\max }\limits_{i = 1,2, \cdot \cdot \cdot ,M} \left| {\Delta {x_i}} \right|\! =\! \alpha {L_{\rm{R}}}$ ，则 $\left| {\displaystyle\frac{{\Delta {x_i}}}{{\left( {C\left( t \right) + 1} \right)v}}f} \right|\!< $ $\left| {\displaystyle\frac{{\alpha {L_a}}}{{\left( {C\left( t \right) + 1} \right)v}}{B_{{\rm{fd}}}}} \right| \approx 0.886\alpha $ 。受天线最小尺寸的限制，L_R不能过小即 $\alpha $ 不会太大，可合理假定 $\alpha $ 在10¹数量级。记 $\varepsilon \left( t \right) \!=\! \left| {\displaystyle\frac{{C\left( t \right) - {C_0}}}{{{C_0} + 1}}} \right|$ ，仿真分析表明 $\mathop {\max }\limits_{t \in \left( { - \frac{{{T_a}}}{2},\frac{{{T_a}}}{2}} \right)} \left( {\varepsilon \left( t \right)} \right)$ 一般不超过10^–3量级，因而 $\Delta \varphi $ 在 $2{{π}} \cdot {10^{ - 2}}$ 量级，为一极小相位。上述分析表明可以忽略式(14)的近似误差，即C₀是 $C\left( t \right)$ 的合理近似。

3 方位向信号重构的矩阵求逆算法

文献[8]给出了单基HRWS-SAR系统方位向信号重构的矩阵求逆算法，由于双基HRWS-SAR系统在方位向信号模型上与单基系统的等效性，矩阵求逆算法在双基系统同样适用。

3.1 双基HRWS-SAR系统的矩阵求逆算法

矩阵求逆算法的理论依据是广义采样定理。在双基HRWS-SAR系统中，通过发射到接收通道i的传递函数，能够得到单基单通道SAR的信号 $S\left( f \right) = U\left( f \right) \cdot H\left( f \right)$ 在接收通道i的表达式 ${S_i}\left( f \right) = $ $ {G_i}\left( f \right) \cdot S\left( f \right)$ 。对于方位向有M个通道的系统，在无重叠采样的情况下，可以得到 $S\left( f \right)$ 的M种独立表达，且采样率均为PRF。根据广义采样定理，能够恢复带宽最高为 $M \cdot {\rm{PRF}}$ 的信号 $S\left( f \right)$ 。

矩阵求逆算法主要包括3个步骤：

首先，构造系统的传递函数矩阵 ${{G}}\left( f \right)$ 及响应矩阵 ${\tilde{{S}}}\left( f \right)$ ：

${{G}}\left( f \right) = \left( {\begin{array}{*{20}{c}}{{G_1}\left( f \right)} & {{G_1}\left( {f + {\rm{PRF}}} \right)} & ·\!·\!· & {{G_1}\left( {f + \left( {M - 1} \right) \cdot {\rm{PRF}}} \right)}\\{{G_2}\left( f \right)} & {{G_2}\left( {f + {\rm{PRF}}} \right)} & ·\!·\!· & {{G_2}\left( {f + \left( {M - 1} \right) \cdot {\rm{PRF}}} \right)}\\ \vdots & \vdots & \ddots & \vdots \\{{G_M}\left( f \right)} & {{G_M}\left( {f + {\rm{PRF}}} \right)} & ·\!·\!· & {{G_M}\left( {f + \left( {M - 1} \right) \cdot {\rm{PRF}}} \right)}\end{array}} \right)$

(16)

${\tilde{{S}}}\left( f \right) = {\Big( {{S_1}\left( f \right),{S_2}\left( f \right), \cdots ,{S_M}\left( f \right)} \Big)^{\rm{T}}}\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\,\, $

(17)

然后，对 ${{G}}\left( f \right)$ 进行求逆运算，得到信号重构矩阵 ${{P}}\left( f \right)$ ：

${{P}}\left( f \right) = \\ \left( {\begin{array}{*{20}{c}}{{P_1}\left( f \right)} & {{P_2}\left( f \right)} & ·\!·\!· & {{P_M}\left( f \right)}\\{{P_1}\left( {f + {\rm{PRF}}} \right)} & {{P_2}\left( {f + {\rm{PRF}}} \right)} & ·\!·\!· & {{P_M}\left( {f + {\rm{PRF}}} \right)}\\ \vdots & \vdots & \ddots & \vdots \\{{P_1}\left( {f + \left( {M - 1} \right) \cdot {\rm{PRF}}} \right)} & {{P_2}\left( {f + \left( {M - 1} \right) \cdot {\rm{PRF}}} \right)} & ·\!·\!· & {{P_M}\left( {f + \left( {M - 1} \right) \cdot {\rm{PRF}}} \right)}\end{array}} \right)$

(18)

最后，利用公式 $P\left( f \right)\cdot S\left( f \right)=\text{ }(\text{ }S\left( f \right),S\left( f+ \right.$ ${ {\left. {{\rm{PRF}}} \right), ·\!·\!· ,S\left( {f + \left( {M - 1} \right) \cdot {\rm{PRF}}} \right)} \big)^{\rm{T}}}$ 即可重构出方位向信号的多普勒频谱。矩阵求逆算法进行方位向重构的原理框图及物理实现原理可参看文献[8]。

3.2 矩阵求逆算法的重构性能

星载HRWS-SAR系统的单通道信号带宽为PRF，采用矩阵求逆算法重构后的方位向信号带宽为 $M \cdot {\rm{PRF}}$ ，因而方位向信号重构必然影响到信噪比和方位模糊比这两个重要的SAR系统性能指标。

3.2.1 信噪比缩放因子

文献[8]定义信噪比缩放因子表征方位向重构对信噪比的影响，其定义式及计算公式为：

${\varPhi _{{\rm{bf}}}} \!=\! \frac{{{\rm{SN}}{{\rm{R}}_{{\rm{in}}}}/{\rm{SN}}{{\rm{R}}_{{\rm{out}}}}}}{{\left( {{\rm{SN}}{{\rm{R}}_{{\rm{in}}}}/{\rm{SN}}{{\rm{R}}_{{\rm{out}}}}} \right)\left| {_{{\rm{PR}}{{\rm{F}}_{{\rm{uni}}}}}} \right.}} \!=\! \frac{{\displaystyle\sum\limits_{j = 1}^M {{\mathop{\rm E}\nolimits} \left( {{{\left| {{P_j}\left( f \right)} \right|}^2}} \right)} }}{N}$

(19)

其中，SNR_in, SNR_out为方位向信号重构前后的信噪比，PRF_uni表示方位向均匀采样时的PRF, ${P_j}\left( f \right)$ 对应于式(18)中的第j行，为作用于接收通道j的重构函数， ${\mathop{\rm E}\nolimits} \left( \cdot \right)$ 为求期望运算。 ${\varPhi _{{\rm{bf}}}}$ 表征了SNR_in与SNR_out的比值随PRF的变化规律。由式(19)可知，越小的 ${\varPhi _{{\rm{bf}}}}$ 意味着更多的信噪比改善。

${\varPhi _{{\rm{bf}}}}$ 还可以通过式(20)进行计算：

${\varPhi _{{\rm{bf}}}} = \frac{{\displaystyle\sum\limits_{j = 1}^M {{\lambda _j}\left( f \right)} }}{N}$

(20)

其中， ${\lambda _j}\left( f \right)$ 为矩阵 ${{P}}\left( f \right) \cdot {{{P}}^{\rm{H}}}\left( f \right)$ 的特征值， ${{{P}}^{\rm{H}}}\left( f \right)$ 为 $\, {{P}}\left( f \right)$ 的共轭转置。

3.2.2 方位模糊比

文献[8]详细推导了单基HRWS-SAR系统的方位模糊比(AASR)计算公式，与之类似可推导出双基HRWS-SAR系统的AASR计算公式，可概括为如下4个计算式：

$\begin{aligned}\\{e_k}\left( f \right) = & {A_k}\left( f \right)\sum\limits_m^{} {\sum\limits_{j = 1}^M {{G_j}\left( {f + k \cdot {\rm{PRF}}} \right)} } \\\text{} & \ \cdot {P_j}\left( {f + m \cdot {\rm{PRF}}} \right)\end{aligned}$

(21)

$\begin{aligned} {e_\Sigma }\left( f \right) = & \sum\limits_{k = - \infty ,k \ne 0}^\infty {{e_k}\left( f \right)} \\= & \sum\limits_{k = - \infty ,k \ne 0}^\infty {\Biggr( {{A_k}\left( f \right)\sum\limits_m^{} {\sum\limits_{j = 1}^M {{G_j}\left( {f + k \cdot {\rm{PRF}}} \right)} } }} \\ \text{ } & {\cdot {P_j}\big( {f + m \cdot {\rm{PRF}}}\big)}\Biggr) \end{aligned}\quad\quad$

(22)

${p_{\rm{s}}} = {\rm{E}}\left( {{{\Big| {A\left( f \right) \cdot {\rm{rect}}\left( {f/{{I}}} \right)} \Big|}^2}} \right)\quad\quad\quad\quad\quad\quad\quad \$

(23)

${\rm{AASR}} = \frac{{{\rm{E}}\left( {{{\left| {{e_\Sigma }\left( f \right)} \right|}^2}} \right)}}{{{p_{\rm{s}}}}}\quad\quad\quad\quad\quad\quad\quad\quad\quad\quad\$

(24)

$A\left( f \right)$ 为发射天线和接收通道合成的天线方向图，双基系统中t_c≠0时 $A\left( f \right) \ne A\left( { - f} \right)$ ，这是导致双基系统AASR计算公式不同于单基系统的根本原因。 ${A_k}\left( f \right)$ 为区间 $\left[ { - {\rm{PRF}}/{2}{\rm{,}}\ {\rm{PRF}}/{2}} \right] + k \cdot {\rm{PRF}}$ 上的 $A\left( f \right)$ , ${e_k}\left( f \right)$ 为区间 $\left[ { - {\rm{PRF}}/{2}{\rm{,}}\ {\rm{PRF}}/{2}} \right] + k \cdot {\rm{PRF}}$ 内的信号产生的方位模糊信号，I为区间 $\left[ - M \right. $ $\left. \cdot {{\rm{PRF}}/{2}}{\rm{,}}{M \cdot {\rm{PRF}}/{2}} \right]$ , p_s为信号能量。上式中m的取值规则为：k>0时， $\max \left( {M - k + 1,1} \right) \le m \le N$ ; k<0时， $1 \le m \le - k$ 。

4 仿真分析

本节对几种典型的双基构型星载HRWS-SAR系统进行方位向信号重构进行仿真，验证矩阵求逆算法的正确性，并分析算法的重构性能。表1列出了系统的方位向系统参数，对于单基系统容易计算出多普勒带宽B_fd=5.61 kHz，照射时间T_a=1.05 s。

表2列出了7种典型的双基构型。表2中采用发射天线和接收天线的零多普勒时间差t_fd与轨道距离L对双基构型进行表征：t_fd=0时系统工作在平飞模式，L=0时系统工作在顺飞模式，L<0时发射天线位于测绘带近端，反之L>0时接收天线位于测绘带近端，构型Ⅰ中t_fd=0且L=0退化为单基系统。

4.1 C₀对C(t)的近似性能

将 $C\left( t \right)$ 近似为C₀是构建双基星载HRWS-SAR系统方位向信号模型及建立矩阵求逆算法的关键步骤，该近似的合理性直接影响到模型及算法的正确性，第3.3小节定性分析了该近似的合理性，本节结合具体的仿真条件予以验证。图3给出了各双基构型下在照射时间范围内的 $C\left( t \right)$ 及 $\varepsilon \!\left( t \right)$ 的变化曲线，图3中对近似误差相对较大的构型Ⅱ、构型Ⅲ采用蓝色线标注，其他构型使用红色线。图3(a)中各双基构型下 $C\left( t \right)$ 曲线近似为直线，直观地的反映出 $C\left( t \right)$ 近似为常数；而图3(b)、图3(c)表明 $\varepsilon \left( t \right)$ 在各双基构型下均小于10^–3，验证了C₀为 $C\left( t \right)$ 的合理近似。图3也表明出目标照射时间T_a及收发天线的零多普勒时间差t_fd是影响 $\varepsilon\! \left( t \right)$ 的主要因素， $\mathop {\max }\limits_t \left( {\varepsilon \left( t \right)} \right)$ 的取值与T_a及t_fd的取值呈正相关性。

4.2 方位向信号重构结果

各接收通道的信噪比设置为20 dB，取PRF=2.0 kHz，对图1所示的点目标P方位向信号重构过程进行仿真，其结果如图4和图5所示。

图4(a)、图4(b)中，由于PRF小于多普勒带宽B_fd，单个接收通道的多普勒频谱有严重的混叠现象，图4(c)则表明采用矩阵求逆算法进行方位向信号重构消除了多普勒频谱中的混叠现象。图4(c)及图5给出了各双基构型下的方位向信号重构后的成像结果，验证了矩阵求逆算法对一般双基构型的适用性，同时也表明重构性能受双基构型的影响，尤其是构型Ⅴ、Ⅶ与其他构型性能差异较为明显。

4.3 方位向信号重构性能

将PRF设置在区间1.4 kHz≤PRF≤2.8 kHz，采用表1的系统参数和表2的双基构型，利用式(20)、式(24)可以得到方位向信号的重构性能曲线，如图6所示。为了便于图6的分析，首先计算各双基构型下的C₀及典型PRF值，计算结果为表3，其中PRF_uni为均匀采样时的PRF值，PRF_rep1及PRF_rep2为PRF范围内出现重叠采样的两个PRF值。构型Ⅰ的PRF_uni, PRF_rep1和PRF_rep2在图6中予以标注。

对比图6和表3，可得到如下结论：

(1) 均匀采样时PRF=PRF_uni，方位向信号的重构性能达到局部最优；重叠采样时PRF=PRF_rep，信噪比缩放因子趋向于无穷大，矩阵求逆算法将不能实现信号重构；PRF位于PRF_rep附近时，矩阵求逆算法的性能急剧下降。从矩阵理论的观点分析，PRF=PRF_rep时，传递函数 ${{G}}\left( f \right)$ 为非满秩矩阵，因而 ${{G}}\left( f \right)$ 及重构矩阵 ${{P}}\left( f \right)$ 的条件数为无穷大，此时信噪比缩放因子 ${\varPhi _{{\rm{bf}}}}$ 也为无穷大，从而无法完成方位向信号的重构。PRF趋于PRF_rep时，重构矩阵 ${{P}}\left( f \right)$ 的条件数急剧增大，方位向信号的重构性能也急剧下降。

(2) C₀值的差异直接导致同一PRF时方位向信号重构性能差异。与其他构型相比构型Ⅴ、Ⅶ的C₀值差异较大因而构性能差异也大，与图5的结论相一致。根本原因在于与双基HRWS-SAR系统等效的单基系统，通道间的方位向距离为 $2/{{\left( {{C_0} + 1} \right)}}$ $ \cdot \Delta {x_i}$ , C₀的取值直接影响到方位向采样的均匀性，进而影响到矩阵求逆算法的重构性能。

5 结束语

本文通过分析单基与双基星载HRWS-SAR系统的方位向信号模型，给出了适用于一般双基构型星载HRWS-SAR系统方位向信号重构的矩阵求逆算法，并使用信噪比缩放因子及方位模糊比两个指标分析了该算法的重构性能。本文的分析方法对其他方位向信号重构算法推广到双基星载HRWS-SAR系统具有借鉴意义。双基星载HRWS-SAR系统方位向信号重构的工程应用中，需改进矩阵求逆算法以改善重叠采样附近的重构性能。