微波雷达成像

2R j sr (t ) Ai p (t i )e c i
2 f c Ri c
(2.1)
Ai 和 Ri (tm ) 分别为第 i 个散射点回波的幅度和某时刻的距离； p () 为归一化的回 波包络； f c 为载波频率， c 为光速。 若以单频脉冲发射，脉冲越窄，信号频带越宽。但发射很窄的脉冲，要有很 高的峰值功率，实际困难较大，通常都采用大时宽的宽频带信号，接收后通过处 理得到窄脉冲。为此，我们将(2.1)式的回波信号换到频域来讨论如何处理，这时 有：
(2.3)
实际 P( f ) 的频带虽然较宽，但总是带限信号，所以一种实用距离成像方法
2
是通过匹配滤波，主要将各频率分量的相位校正成一样，为了提高信噪比再按信 号频谱幅度加权，而频谱为零部分是无法恢复的。匹配滤波后的输出为，
* sr M (t ) F(f1) Sr ( f ) P ( f ) 2 ( f c f ) Ri j c F(f1) Ai P( f ) P* ( f )e i
c 。 2
r r 应当指出，如 r 一定，则解线调频后的频率范围为 , ，即信号 c c
最大频宽为
2r 2r r = B= B ，其中 R p 为 Tp 所对应的距离。因此可见，比值 cT p Rp c
6
r 越小，则信号最大频宽比原调频带宽也小得越多，在聚束式 SAR 和 ISAR 里 Rp
2 R ，可得 c
c 1 c 1 ，相应的时间分辨率为 1 B ，这与匹配滤波 2 T p 2 B
脉冲压缩的结果是一致的。
发射频率 ft fc B
场景中心回波
近距离回波
参考信号
远距离回波
Tp
t 0
2r c
(a) 解线频调
t
解线频调后频率 ft
0 0
2 Ra c r c
Tref
Tp
近距离回波
T T p 的时间长度，作傅立叶变换后并作共轭，和接收信号的傅立叶变换相乘
3
后，作傅立叶逆变换，取前 T 时间段的有效数据段。为了便于采用快速傅立叶 变换，可能对匹配函数要补更多的零，对接收信号也要补零。脉压处理过程的如 图 2-2 所示，其中虚框部分可事先计算好，以减小运算量。
Tp
参考信号
FFT
共轭相乘

IFFT
T Tp
FFT
图 2-2 匹配滤波脉压示意图
接收信号
距离匹配滤波压缩后，不管是否补零，其距离分辨率为 c (2 B) ，距离采样率 为 c (2 Fs ) ，其中 Fs 为采样频率， Ts 于距离分辨单元长度。 2.2 线性频调信号和解线频调处理 大时宽宽频带信号可以有许多形式，如脉冲编码等，但用得最多的是线性调 频（LFM）脉冲信号。由于线性调频信号的特殊性质，对它的处理不仅可用一般 的匹配滤波方式，还可用特殊的解线频调（Dechirping）方式来处理。 解线频调脉压方式是针对线性调频信号提出的， 对不同延迟时间信号进行脉 冲压缩，在一些特殊场合，它不仅运算简单，而且可以简化设备，已广泛应用于
即
2 R t ˆ 2 Rt c j 4c ( tˆ 2 Rcref ) R j 4c f c R j 4c 2 e ˆ s if (t , t m ) Arect e e T p
(2.9)
若暂将讨论限制在一个周期里（即 R 为常数） ，则上式为频率与 R 成正比的单 频脉冲。如果所需观测的范围为 [ Rref r 2 , Rref r 2] ，图 2-3 中画出了范围两 侧边缘处的回波。 我们再结合，图 2.3 是解线频调的差频处理示意图作一些说明，图中纵坐标 均为频率，图 2.3(a)中除参考信号外，有远、近的两个回波。参考信号与回波作 其共轭相乘，即作差频处理，回波变成单频信号，且其频率与回波和参考信号的 距离差成正比，因而也叫解线频调处理。由图 2-3(b)可知 f i
S r ( f ) Ai P ( f )e
i
j 2 ( f c f ) Ri c
(2.2)
对理想的几何点目标当然希望重建成冲激脉冲，如果 P( f ) 在所有频率没有零 分量，则冲激脉冲信号可通过逆滤波得到，即 F
1 ( ) 2 f c j Ri Sr ( f ) 2R c e A (t i ) i P( f ) c i
2 2 Ri 1 2 Ri ˆ 2 t c c
(2.8)
解线频调的示意图如图 2-3，若 R Rt Rref ，则其差频输出为
* s if (tˆ, t m ) s r (tˆ, t m ) s ref (tˆ, t m )
(2.7)
式中 Tref 为参考信号的脉宽，它比 T 要大一些(参见图 2-3)。
ˆ, t m ) 为 某点目标到雷达的距离为 Rt ，雷达接收到的该目标信号 s r (t
f t tˆ 2 Rt c j 2 c e ˆ, t m ) Arect s r (t T p
这一比值有时小到几十分之一，甚至几百分之一，以 ISAR 为例，飞机一类目标 的长度一般小于 100 米，对应的时宽为零点几微秒，而大时宽的宽频带信号一般 在几十微秒以上，从而可将信号频带从几百兆赫减小到只有几兆赫，对后续设备
(2.6)
1 其中 rect(u ) 0
u u
1 2 1 2
， f c 为中心频率， Tp 为脉宽， 为调频率， t t mT 为
4
快时间， m 为整数， T 脉冲重复周期， t m mT 为慢时间。 解线频调是用一时间固定， 而频率、 调频率相同的 LFM 信号作为参考信号， 用它和回波作差频处理。设参考距离为 Rref ，则参考信号为
1
3500
3000
2500
2000
1500
1000
500
0
100
110
120
130
140
150
160
图 2-1
飞机回波的一维距离像
不同散射点的径向距离差改变，从而使两者子回波的相位差可能显著变化。以波 长 3 厘米为例，若两散射点的横距为 10 米，当目标转动 0.05°时，两者到雷达 的径向距离差变化为 1 厘米，它们子回波的相位差改变 240°！由此可见，目标 一维距离像中尖峰的位置随视角缓慢变化（由于散射点模型缓变） ，而尖峰的振 幅可能是快变的（当相应距离单元中有多个散射点） 。图 2-1 是 C 波段雷达实测 的飞机一维距离像的例子，图中将视角变化约 3°的回波重合画在一起。一维距 离像随视角变化而具有的峰值位置缓变性和峰值幅度快变性可作为目标特性识 别的基础。 本章将用上述散射点模型对高分辨的一维距离像进行讨论。 2.1 宽带信号的逆滤波、匹配滤波和脉冲压缩 根据散射点模型，设散射点为理想的几何点，若发射信号为 p(t ) ，对不同 距离多个散射点目标，其回波可写成：
2 R 。因此，对 c
解线频调后的信号作傅立叶变换， 便可在频域得到对应的各回波的 sinc 状的窄脉 冲，脉冲宽度为 1 Tp ，而脉冲位置与 R 成正比( 示。
2 R )，如图 2-3(b)的左侧所 c
5
如上所述，变换到频域窄脉冲信号的分辨率为 1 Tp ，利用 f i 相应的距离分辨率为 r =
(2.4)
Ai e
i
j
2 f c Ri c
wk.baidu.com
psf (t
2 Ri ) c
这里 P* () 为 P() 的复共轭，而
2 psf (t ) F(f1) P ( f )
(2.5)
在时域上看，滤波相当于信号与滤波器冲激响应的卷积，对一已知波形的信 号作匹配滤波，其冲激响应为该波形的共轭倒置。当波形的时间长度为 T p ，则 卷积输出信号为 2T p 。实际上，匹配滤波可实现脉冲压缩，输出主瓣的宽度为
2r c
场景中心回波 远距离回波
t
(b)
ft
ft
0
2r c 近距离回波 Tp
场景中心回波 远距离回波
t
(c )
2 Ra c r c
图 2-3 解线频调脉压示意图
由于用解线频调作脉冲压缩的结果表现在频域里， 而不像匹配滤波是在时域 里完成，有些书籍里又把这种方法叫“时频变换脉冲压缩” 。从频率域变换到距 离(相对于参考点的)，应乘以系数
第二章 距离高分辨和一维距离像
雷达采用了宽频带信号后，距离分辨率可大大提高，这时从一般目标（如飞 机等）接收到的已不再是“点”回波，而是沿距离分布开的一维距离像。 雷达回波的性质可以用线性系统来描述， 输入是发射脉冲， 通过系统 （目标） 的作用，输出雷达回波。系统的特性通常用冲激响应（或称分布函数）表示，从 发射波形与冲激响应的卷积可得到雷达回波的波形。 严格分析和计算目标的冲激响应是比较复杂的，要用到较深的电磁场理论， 不属于本书的范围。简单地说，雷达电波作用的目标的一些部件对波前会有后向 散射，当一些平板部分面向雷达时还会有后向镜面反射；这些是雷达回波的主要 部分；此外还有谐振波和爬行波等。因此，目标的冲激响应（分布函数）可以用 散射点模型近似，即目标可用一系列面向雷达的散射点表示，这些散射点位于后 向散射较强的部位。由于谐振波和爬行波的滞后效应，有时也会有少数散射点在 目标本体之外。如上所述，目标的散射点模型显然与雷达的视线向有关，例如当 飞机的平板机身与雷达射线垂直时有很强的后向镜面反射， 而在偏离不大的角度 后，镜向反射射向它方，不为雷达所接收。目标的雷达散射点模型随视角的变化 而缓慢改变，且与雷达波长有关，分析和实验结果表明，在视角变化约 10°的 范围里，可认为散射点在目标上的位置和强度近似不变。顺便提一下，前面曾提 到微波雷达对目标作 ISAR 成像，目标须转动 3°左右，在分析时用散射点模型 是合适的。 虽然目标的散射点模型随视角 作缓慢变化， 但一维距离像的变化要 快得多。可以想像到，一维距离像是 三维分布散射点子回波之和， 在平面 波的条件下， 相当三维子回波以向量 和的方式在雷达射线上的投影， 即相 同距离单元里的子回波作向量相加。 我们知道， 雷达对目标视角的微小变 化， 会使同一距离单元内而横向位置
1 B ( B 为信号的频带宽度，为降低副瓣而作加权，主瓣要展宽一些)，即距离分
辨率为 c (2 B) ， 脉压信号的 B 通常较大( BT 1 )， 输出主瓣是很窄的， 时宽为 2T p 的输出中，绝大部分区域为幅度很低的副瓣。 当反射体是静止的离散点时， 回波为一系列不同延时和复振幅的已知波形之 和，对这样的信号用发射波形作匹配滤波时，由于滤波是线性过程，可分别处理 其相当的时间段为 T (= 2 r c )， 后迭加。 如果目标长度相应的回波距离段为 r ， 考虑到发射信号时宽为 T p ，则目标所对应的回波时间长度为 T T p ，而匹配滤 波后的输出信号长度为 T 2T p 。虽然如此，具有离散点主瓣的时间段仍只有
T ，两端的部分只是副瓣区，没有目标位置信息。
应当指出，通过卷积直接作匹配滤波脉压的运算量相对较大，可以在频率域 通过共轭相乘再作 IFFT 求得。需要注意的是两离散信号频率域相乘相当它们在 时域作圆卷积，为使圆卷积与线性卷积等价，待处理的信号须加零延伸，避免圆 卷积时发生混叠。 实际处理中，为了压低副瓣，通常是将匹配函数加窗，然后加零延伸为
1 为采样周期，距离采样周期要求小于等 Fs
SAR 和 ISAR 中作脉冲压缩。应当指出，解线频调处理和匹配滤波虽然基本原理
相同，但两者还是有些差别的，为了能正确利用解线频调方式作脉冲压缩，我们 对它作一些详细的说明。 假设发射信号为
t j 2 ( f t 1 t 2 ) c 2 s(t, t m ) rect ， T e p
fc t tˆ 2 Rref c j 2 ˆ s ref (t , t m ) rect e Tref 2 Rref 1 2 Rref ˆ t 2 c c
2