SAR信号处理的学习总结(三)
前言
关于SAR雷达的基本模型和基础的介绍可以看(一)SAR的基本模型;SAR的相关信号介绍和计算可以看(二)如何得到回波信号;。如果遇到未说明的参数,请看(一)、(二)中的说明。本篇主要对SAR信号的相关的操作进行介绍。
快时间与慢时间
在(二)中得到了SAR的回波信号:
S
0
(
τ
,
η
)
=
A
0
⋅
w
r
(
τ
−
2
R
(
η
)
/
c
)
⋅
w
α
(
η
−
η
c
)
⋅
exp
{
j
π
K
r
(
τ
−
2
R
(
η
)
/
c
)
2
−
j
4
π
f
0
R
(
η
)
/
c
}
S_0\left( \tau ,\eta \right) =A_0\cdot w_r\left( \tau -2R\left( \eta \right) /c \right) \cdot w_{\alpha}\left( \eta -\eta _c \right) \cdot \exp \left\{ j\pi K_r\left( \tau -2R\left( \eta \right) /c \right) ^2-j4\pi f_0R\left( \eta \right) /c \right\}
S0(τ,η)=A0⋅wr(τ−2R(η)/c)⋅wα(η−ηc)⋅exp{jπKr(τ−2R(η)/c)2−j4πf0R(η)/c}
(二)中仅仅从信号处理的角度得到了该信号,但就从信号接收的角度来说,现实中接收信号来说只有一个时间,又是如何分出了快时间还是慢时间。
快时间和慢时间实际就是将一个随现实时间流动的信号根据多次发射的不同分为了两个维度,一个脉冲回波的幅值变化(快时间)以及多个脉冲的发射时间(慢时间)。
脉冲压缩
脉冲压缩的目的
获得高信噪比的途径通常有两条,增大峰值功率与增加发射信号的持续时间。
高峰值功率较难实现,通常的办法为增加持续时间,但是由于高分辨率(分辨率值小)与发射信号持续时间之间为负相关;有一个很显然的例子:雷达信号中距离向的分辨率
ρ
r
≈
c
2
B
=
c
⋅
T
p
2
\rho _r\approx \frac{c}{2B}=\frac{c \cdot T_p}{2}
ρr≈2Bc=2c⋅Tp其中
B
B
B为信号的带宽,
T
p
T_p
Tp为脉冲的持续时间;脉冲持续时间越短,分辨率越小。
而脉冲压缩的目的正是来解决信号持续时间与分辨率之间的矛盾。即在发射一个时宽脉冲信号的同时获得较高分辨率。
脉冲压缩的实现方法
-
如何在有限条件下获得一个短脉冲?
在实际信号发射的过程中,尽量去延展信号的带宽是一个不错的选择,即将能量平均于每一个频率分量。因此频域上均匀的信号 X X X
X ( f ) = r e c t ( f T ) = { 1 , i f ∣ f T ∣ ≤ 1 2 0 , e l s e X\left( f \right) =rect\left( fT \right) =\left\{ \begin{array}{l} 1,\ if\left| fT \right|\le \dfrac{1}{2}\\ 0,\ else\\ \end{array} \right. X(f)=rect(fT)={1, if∣fT∣≤210, else得到最短持续时间的脉冲
x = I D F T ( X ) = sin c ( t T ) = sin ( π t T ) π ⋅ t T x=IDFT\left( X \right) =\sin c\left( \frac{t}{T} \right) =\frac{\sin \left( \pi \frac{t}{T} \right)}{\pi \cdot \frac{t}{T}} x=IDFT(X)=sinc(Tt)=π⋅Ttsin(πTt) s i n c sinc sinc函数在时域能量分布非常的集中,其时间带宽积近似为1,具有良好的性质特点。
由此,对于一般的脉冲压缩过程,其目的在于使得信号在频域上的幅值近似于矩形函数,信号的相位线性变化(仅包含常量和线性分量)。 -
脉冲压缩方法
频域上,脉冲压缩本质上是消除信号频谱中含有二次相位的部分,即信号与含有二次共轭相位的滤波器在频域相乘。该过程与信号滤波一致,消除高次分量,保留低次分量,因此脉冲压缩也叫匹配滤波。
时域上,匹配滤波可以通过相关操作实现, S = ∫ − ∞ ∞ s r ( u ) g ∗ ( u − t ) d u S=\int_{-\infty}^{\infty}{s_r\left( u \right) g^*\left( u-t \right)}du S=∫−∞∞sr(u)g∗(u−t)du其中 S r S_r Sr为接收信号, g ∗ g^* g∗为预期信号(发射信号)的复共轭。
对应的匹配滤波器为 h ( t ) = g ∗ ( − t ) h\left( t \right) =g^*\left( -t \right) h(t)=g∗(−t) -
线性调频信号的压缩
发射信号: s p u l = r e c t ( t T ) ⋅ exp ( j π K t 2 ) s_{pul}=rect\left( \frac{t}{T} \right) \cdot \exp \left( j\pi K t ^2 \right) spul=rect(Tt)⋅exp(jπKt2)
延时 t 0 t_0 t0后的目标回波:
s r = r e c t ( t − t 0 T ) ⋅ exp ( j π K ( t − t 0 ) 2 ) s_r=rect\left( \frac{t-t_0}{T} \right) \cdot \exp \left( j\pi K\left( t-t_0 \right) ^2 \right) sr=rect(Tt−t0)⋅exp(jπK(t−t0)2)
则对应匹配滤波器 h ( t ) h\left(t \right) h(t)为 t 0 = 0 t_0=0 t0=0时的时间反褶后 s r s_r sr的复共轭
h ( t ) = s r ∗ ( − t ) = r e c t ( t T ) exp ( − j π K t 2 ) h\left( t \right) =s_{r}^{*}\left( -t \right) =rect\left( \frac{t}{T} \right) \exp \left( -j\pi Kt ^2 \right) h(t)=sr∗(−t)=rect(Tt)exp(−jπKt2)
压缩后的输出
s o u t = T ⋅ s i n c ( K T ( t − t 0 ) ) s_{out}=T\cdot sinc\left( KT\left( t-t_0 \right) \right) sout=T⋅sinc(KT(t−t0))