王裕旗 孙光才* 杨 军 邢孟道 杨小牛 保 铮

①(西安电子科技大学雷达信号处理国家重点实验室 西安 710071)

②(西安科技大学测绘科学与技术学院 西安 710071)

③(中国电子科技集团公司第三十六研究所 嘉兴 314000)

1 引言

在无源电子侦察中，辐射源精确定位具有重要的意义[1-4]。首先，确定通信基站、雷达等辐射源位置后，可以对这些电子设备实施精确打击；其次，确定辐射源位置后，可以降低干扰信号的发射功率对目标辐射源实施更加精确的干扰；最后，信号分选是电子对抗中的基本问题，由于辐射源位置是不能发生捷变的，因此国内已经有学者认为它可以成为信号分选可靠和强有力的依据[4]。根据观测站的数量可以将无源定位技术分为两类：单站无源定位和多站无源定位。

多站无源定位主要通过定位参数获得辐射源的位置，其中主要有多站测向交叉定位、多站时差定位、多站频差定位和时差频差联合定位[5-9]。这类定位方法的优点是其不受信号形式的限制，但需要两个以上的基站，基站间数据的传输和时间的同步等问题增加了系统的设计复杂度。单站定位基本的原理是通过测向获得辐射源的位置[10-12]。传统电子侦察设备主要通过干涉仪[13-17]进行测角定位，这种方法依赖于干涉仪的性能[14,15]，天线的极化误差[16]、基线长度都会影响干涉仪的性能，并且当侦察场景中存在多个同频辐射源时[17]，测角性能会下降，测角定位的方位分辨能力也会随着距离的增加而下降。通过单个基站的运动，可以实现多个位置测角联合定位[18]，也可以通过平台运动过程中测量侦察信号的多普勒和多普勒的变化率完成辐射源的定位[19-22]，这种定位方法以最简单的设备完成辐射源定位，但考虑到对多普勒和其变化率的测量是基于多个时间段的数据单独测量，没有利用长时间数据的相参积累，这种方法可以称为非聚焦类的无源定位方法。受合成孔径雷达信号相参积累的启发，本文提出一种基于长合成孔径的辐射源定位方法，可以将其称为聚焦类的无源定位方法。

合成孔径雷达(Synthetic Aperture Radar,SAR)是一种全天候、全天时的高分辨成像雷达，是一种新的观测手段，因此在军事侦察中应用广泛[23-33]。雷达载体沿方位向飞行过程中，发射和接收信号形成合成阵列。为了获得高的距离分辨率，合成孔径雷达发射的信号频带必须较宽；为了提高方位分辨率，必须加大天线孔径，对于固定的场景，可以通过合成阵列的办法得到长的孔径[24]。目前国内外的文章未有提及合成孔径辐射源侦察的相关概念，本文将长合成孔径的概念引入无源定位中(下文中称为长合成孔径定位)，在侦察机飞行的过程中连续采集一段数据，然后按照合成孔径雷达的工作方式来利用数据，通过在快时间域进行傅里叶变换，可以得到发射信号的频率，然后在慢时间域分析采集到的数据的相位变化历程，对于不同频率，利用驻定相位原理[34]直接在多普勒域构造匹配滤波器，通过聚焦得到辐射源距离位置信息和辐射源方位位置信息。本文针对常见的单频连续波辐射源，首先建立了合成孔径侦察的几何模型和描述了侦察机采集数据的利用方式，给出了侦察机接收到的信号在快时间域和慢时间域的表达式，详细介绍了合成孔径侦察算法，最后通过实测数据的处理对该算法进行了验证。

2 几何模型与采集数据利用方式

长合成孔径无源定位时，雷达不发射信号，只接收外部辐射源发射的信号。图1为长合成孔径定位几何模型，雷达位于侦察机的正侧方，波束宽度为θBW，侦察机以速度v平行于地平面做匀速直线飞行，侦察机从A到C飞行的过程中可以接收到辐射源Pi发射的信号，在其他位置都不能接收到Pi发射过来的信号，B为AC中点。以A点为坐标原点，侦察机飞行的方向为X轴(方位向)，垂直于侦察机飞行的方向为Y轴(距离向)，建立坐标系，辐射源Pi的坐标为(Xi,RB,i),Xi为辐射源Pi的方位坐标，RB,i为辐射源Pi的距离坐标，i=1,2,···,N,N为波束范围内辐射源的总数。

图1 合成孔径定位几何模型Fig.1 Localization geometry model

图2为侦察机采集到数据的利用方式。对于侦察机采集到的单频连续信号，每隔时间T取一段持续时间为Tp的数据。对于单频脉冲信号，T和Tp分别为估计获得的脉冲重复周期和脉冲宽度；对于单频连续信号，如果侦察机的速度为100-200 m/s，则T的大小可以取到毫秒级，Tp可以取到微秒级。这与脉冲体制合成孔径雷达中发射的信号形式一样。用tm=mT表示接收到第m个脉冲的方位慢时间，tm=0时刻侦察机处于A点，则(vtm,0)表示tm时刻侦察机所处的位置，Ri(tm)表示在tm时刻侦察机辐射源Pi(Xi,RB,i)的距离，即瞬时斜距。根据几何关系可得

通信基站发射的调频或调幅信号，频谱范围为3-30 kHz，带宽较窄，当采样速率fs较高时，可以看成是单频信号。若辐射源Pi发射单频连续波信号si(t)

其中Ai为信号的幅度，fi为载频，φi为初始相位。用tr表示距离快时间，tm表示方位慢时间，则侦察机采样时刻可用(tr,tm)表示，侦察机在(tr,tm)时刻接收到的辐射源Pi发射的信号可表示为

其中c为光速，第2项相位中包含辐射源的方位坐标信息和方位调频率特性，方位调频率和辐射源的距离有关，因此这一项可以用来进行方位和距离定位。

3 长合成孔径定位算法

合成孔径侦察主要目的是获得辐射源发射信号的频率，以及辐射源的位置，从而用频率，距离，方位这3维信息来表示辐射源。

图2 采集数据的利用方式Fig.2 Data usage diagram

3.1 单频信号定位原理

从式(3)可以得到，若辐射源发射信号为单频连续波，侦察机接收到的信号在快时间域和慢时间域是2维可分离的。对接收到的信号sr,i(tr,tm)在快时间域tr做傅里叶变换，得到

其中f为距离频率，δ(f)为冲激函数，满足δ(f)=0,

由于采样频率的限制，当侦察信号频带范围较大时，用傅里叶变换估计频率时存在频率模糊的问题，可以通过加一组带通滤波器来解决此问题，对滤波后的信号做傅里叶变换，可以得到发射信号的频率fi。从式(4)可以得到，在快时间域傅里叶变换后，可以得到发射信号频率。对于频率f=fi的数据乘以exp{−j2πfitm}后得到Wi(fi,tm)

其中Ki只和Ai,φi有关且。在不同的距离Rk处，构造匹配滤波器h(fi,tm;Rk)，进行距离位置定位

其中fa为多普勒频率。经过匹配滤波器的输出为

psf(·)称为点散布函数(Point Spread Function,PSF)，它可确定方位分辨率，⊗为卷积。经过匹配滤波处理后，当Rk=RB,i时，信号完全聚焦，通过峰值检测可以得到与航线垂直距离为RB,i，方位位置为Xi，发射频率为fi的辐射源。在不同频率处根据不同距离构造匹配滤波器定位距离位置，可以得到不同频点的辐射源的距离和方位位置，达到侦察目的。

在慢时间域tm采样点数一定的情况下，不同的发射信号频率fi，不同的垂直距离RB,i具有不同的点散布函数，方位分辨率也不同。在慢时间域tm采样点数相同条件下，图3给一组不同发射频率fi，不同垂直距离RB,i时的点散布函数，可以看到，发射信号频率fi大，距离RB,i小，方位分辨性能好。实际上，可以通过增加侦察机数据采集时间来提高方位分辨率，理论上方位维分辨率可以达到vT，其中v为侦察机飞行的速度，T为方位慢时间间隔。

对于单频连续波辐射源，合成孔径侦察算法流程可由图4表示，其中hi(tr),i=1,2,···,M为带通滤波器。带通滤波并进行距离向FFT后，可以获得信号的载频，截取不同载频的信号，可以获得频率为fi的多个矩阵sr(fi,tm)。对每一个截取的信号矩阵sr(fi,tm)，在定位距离RB,i的时候，可以采用分步的方式，初始的间隔可以设的较大，确定辐射源的大致位置，然后采用较小的间隔，对辐射源进行精确定位。每一个距离处可以得到不同频率信号矩阵的匹配滤波结果，将不同频率匹配滤波的结果进行叠加作为此距离处的定位结果，然后把不同距离处的定位结果按距离向重新排列，即可得到辐射源的定位结果图。

3.2 宽波束侦察定位

在采用宽波束侦察时，信号接受过程中，雷达与辐射源的斜距历程可以写为

其中，t=tr+tm为全时间，θi为雷达速度矢量和雷达与辐射源Pi连线的距离矢量夹角，R0,i为t=0时刻雷达到辐射源Pi的斜距，R0,i=RB,i/sinθi。当采样率为kHz级别时，信号排列成2维的形式后，为满足停走停的模型，距离向点数只能取很小，这会导致距离频率的分辨率降低，当距离向取的点数较多时，接收信号在距离向和方位向都存在调制。因此，中频采样的信号模型修改为

图3 不同发射频率，不同垂直距离对应的点散布函数Fig.3 PSF for different frequencies and different ranges

图4 单频连续波辐射源合成孔径侦察算法流程Fig.4 Flow chart of proposed algorithm for single frequency signal

在斜视模式下，辐射源的方位角θi也是需要估计的一个量，可以采用短合成孔径的信号估计方位角。当采样的脉冲数较少时，可以忽略距离变化的二次项

此时通过式(14)进行角度估计的误差为

角度的估计误差会带来多普勒调频率的估计误差，因此载频估计的误差会影响距离定位的结果。

对于多个频点的单个辐射源，辐射源的位置是唯一确定的，调频率和信号的频点是唯一对应的，因此可以在能量最高的频点处进行辐射源的定位。

4 仿真结果

4.1 计算量

该算法的计算量来源于两次FFT、一次IFFT和距离定位的乘法。假设原始数据距离向点数为m，方位向点数为n,m和n为2的整数次幂，搜索的场景宽度为w，粗定位的步长为r1，精定位的步长为r2。距离向m点FFT需要进行次复乘和mlog2m次复加；方位向n点FFT需要进行次复乘和nlog2n次复加；进行距离定位时和匹配滤波器相乘的次数为w/r1+r1/r2，每次匹配滤波相乘需要进行n次复乘，共需要进行(w/r1+r1/r2)n次复乘；方位向n点 IFFT需要进行次复乘和nlog2n次复加。因此，共需要进行Log2m+nlog2n+(w/r1+r1/r2)n次复乘和mlog2m+2nlog2n次复加。

4.2 仿真结果

仿真数据的信号为单频信号，发射信号载频为2.6 GHz，采样率为300 MHz，卫星速度为7500 m/s，卫星的工作模式为正侧视，方位采样频率为4000 Hz，信噪比(Signal-Noise Ratio,SNR)为0 dB。场景中存在3个独立的辐射源，辐射源的位置分别为(0,1200 km),(1 km,1192 km)和(3 km,1202 km)，信噪比为3 dB。辐射源的定位结果如图5所示，3个辐射源的位置分别为(0,1200 km),(0.9994 km,1192 km)和(3 km,1202 km)。图6是通过1000次蒙特卡洛试验获得的不同信噪比下，(0,1200 km)处的辐射源定位的均方根误差(Root Mean Square Error,RMSE)。

5 实测数据处理结果

实测数据的信号为QPSK编码信号，载频为2.1 GHz，采样率为500 KHz，信噪比为3-5 dB。在每个码元为单频信号，不同码元的初相不一致，变换到频域后，信号占据一定的带宽。对下变频后的信号直接进行傅里叶变换可以得到图7所示的频域信号。此时信号并没有完全积累，并且信号在方位向的相干关系被破坏，这给方位向的信号积累带来困难。因此，需要补偿码元间的相位调制，实现信号的全相干积累。可以通过前后两个采样点做相位差分，补偿整数倍π/2相位的跳变，完成码元调制相位的补偿。补偿后的信号傅里叶变换后如图8所示，信号只占一个频点。

图6 不同信噪比下定位的均方根误差Fig.6 RMSE of different SNR

图7 QPSK信号频谱Fig.7 QPSK signal spectrum

将补偿后的信号排列成2维的形式，取距离向点数256点，方位向点数8192点，将重排后的信号进行方位向傅里叶变换可以得到信号的方位谱。在长合成孔径条件下，距离等式包含慢时间的二次项，这相当于信号在慢时间有调制，在方位向上相当于线性调频信号，方位傅里叶变换后信号的多普勒谱占一定的带宽，如图9所示。

在不同的距离构造匹配滤波器进行距离位置定位可以得到图10的结果，当匹配函数中的距离和辐射源的真实距离不相等时，方位脉压后信号会散焦，当匹配函数的距离和信号的真实距离匹配时会形成峰值。从图3的分析可以得到距离越近，方位分辨率越高，同时从多普勒调频率的表达式γm(fi;RB,i)=fiv2/(cRB,i)可以看出，距离越近，调频率对距离变化越敏感，即距离定位步长确定时，RB,i越小，调频率γm的变化越明显。

图8 相位补偿后QPSK信号的频谱Fig.8 Spectrum after phase compensation

图9 2维重排信号的方位谱Fig.9 Azimuth spectrum of two-dimensional signal

图10 距离位置定位的结果Fig.10 Range localization result

图11是距离位置定位结果的局部放大图，红色位置为信号的真实位置，得到的辐射源的距离RB,i为1282.0312 km，方位坐标Xi为229.97 m，误差4.5058 k m。

图12是辐射源在距离频域和方位时域的压缩结果，可以得到辐射源信号的频率和方位坐标，结合距离位置定位的结果，可以获得辐射源距离位置、方位位置和信号频率3维信息，完成辐射源的侦察定位。结合卫星GPS的位置信息，可以获得辐射源在经纬度坐标系下的位置为东经115.73°，北纬40.36°。

传统的单星定位通常使用多个天线进行干涉测角定位，使用四通道天线数据进行测角定位，通过地面划分网格搜索的结果如图13所示，定位误差为15.563 km。取不同时间段的数据进行测角定位和合成孔径定位，结果如表1所示。

定位的误差主要来源于实测数据中，GPS和数据录取设备使用不同的时钟，数据采样的位置和GPS给出的位置不能完全对应。对于测角定位的结果和该算法定位的结果，定位的精度可以提高3-4倍。

图11 距离位置定位结果局部放大Fig.11 Enlarged view of range localization result

图12 辐射源的频率和方位位置Fig.12 Frequency and azimuth location

图13 测角定位的结果Fig.13 Localization result of angle measurement

表1 定位误差对比(km)Tab.1 Comparison of localization errors (km)

6 结束语

获取雷达，通信基站等辐射源发射信号的频率范围及位置信息具有重要意义。本文将合成孔径的概念引入无源侦察中，仿照合成孔径雷达中天线的工作方式来利用侦察机采集到的数据，进行辐射源的距离位置和方位位置定位。实测数据的处理表明，对于单频连续波辐射源，长合成孔径定位算法能够以较高精度获得辐射源的方位位置、距离位置等信息，达到对辐射源侦察目的。相对于传统使用多天线干涉测角的定位方法，该方法只需要使用一个天线，降低了系统设计的复杂度。