第十章 合成孔径雷达
10.3 合成孔径雷达的另外一种解释
本节从分析回波信号的特性入手,从频谱分析、相关、匹配滤波角度来说明合成孔径原理。
10.3.1 回波信号特性
假设被测目标为一理想点目标p,p点与航线x的垂直斜距为
R0
R
0
,我们把航线x和
R0
R
0
所构成的平面作为坐标平面。设飞机在
t=0
t
=
0
处在坐标原点。在某一瞬时t,飞机的位置在
xa=vat
x
a
=
v
a
t
。点目标p的位置在这个坐标系里是固定的,其坐标为
(xp,R0)
(
x
p
,
R
0
)
。在t时刻,p点与飞机上雷达天线的斜距R为:
一般情况下,
R0>>(xa−xp)
R
0
>>
(
x
a
−
x
p
)
,上式可近似为
天线发出的是周期性的相干等幅高脉冲波,设其频率为
f0
f
0
,振幅为A,脉冲重复频率为
fr
f
r
,重复周期为
Tr=1/fr
T
r
=
1
/
f
r
,脉冲宽度为
τ
τ
。
①假设发射的为一连续波余弦信号,把实际信号看成是对连续信号的抽样,其抽样率即为脉冲重复频率
fr
f
r
;
②假定余弦信号的振幅归一化为1,起始相位为0,则有:
发射信号
S0(t)=Re(eiwt)
S
0
(
t
)
=
R
e
(
e
i
w
t
)
回波信号
Sr(t)=Re[Kσ0ejw(t−τ0)]
S
r
(
t
)
=
R
e
[
K
σ
0
e
j
w
(
t
−
τ
0
)
]
式中,K表示由距离R及其他因素引起的对信号幅度的衰减因子,
τ0
τ
0
为信号往返延迟。
那么
τ0=2R/c
τ
0
=
2
R
/
c
,
τ0=2c[R0+(xa−xp)22R0]=2R0c+(xa−xp)2cR0
τ
0
=
2
c
[
R
0
+
(
x
a
−
x
p
)
2
2
R
0
]
=
2
R
0
c
+
(
x
a
−
x
p
)
2
c
R
0
,
归一化后,
这里, λ=c/f0 λ = c / f 0 ,
取实部后有:
这个信号的相位部分由三项组成:
φ=φ1+φ2+φ3
φ
=
φ
1
+
φ
2
+
φ
3
φ1
φ
1
:原始发射信号的一次相位(线性相位);
φ2
φ
2
:是随
R0
R
0
而变的相位项,但与时间无关。对同一垂直斜距的目标来讲,
R0
R
0
是常数,
φ2
φ
2
是常数相位;
φ3
φ
3
:最重要的相位项,也可写成
φ3=−2π(xa−xp)2λR0
φ
3
=
−
2
π
(
x
a
−
x
p
)
2
λ
R
0
,随时间呈平方律变化的二次相位项。
将相位对时间进行求导,再除以
2π
2
π
,即得回波信号的瞬时频率:
式中,
f0
f
0
是发射信号的载频。第二项是因为天线与目标有相对运动而引起的多普勒频移,通常用
fd
f
d
表示。
它随时间呈线性变化,可见回波信号是一种线性调频信号,其调频斜率为
从上图可以看出,点目标p引起的多普勒频移有一个范围,以
t=t0
t
=
t
0
为中心向正负两方向变化。当
t=t0
t
=
t
0
时,天线位置正好处在p点与航线的垂直斜距点
fd=0
f
d
=
0
;在
t=t0
t
=
t
0
时刻以前,
t−t0<0→fd>0
t
−
t
0
<
0
→
f
d
>
0
,其最大值发生在
t=t0−Ls/2va=t0−Ts2
t
=
t
0
−
L
s
/
2
v
a
=
t
0
−
T
s
2
,此时的多普勒频移为:
Ls L s 为p点所在位置的合成孔径长度, Ts T s 为合成孔径时间。
在
t=t0
t
=
t
0
时刻以后,
t−t0>0→fd<0
t
−
t
0
>
0
→
f
d
<
0
,其最大负值发生在
t=t0+Ls2va=t0+Ts2
t
=
t
0
+
L
s
2
v
a
=
t
0
+
T
s
2
,这时,多普勒频移为:
那么点目标p的回波多普勒频移的带宽为:
由于
Ts
T
s
可以表示成
故多普勒带宽也可以写成:
上式表明,回波信号的多普勒带宽与飞机航速 va v a 、合成孔径长度 Ls L s 或航向波束宽度 θa θ a 成正比。
从线性调频信号的频谱来考察回波信号的特点,即
把调频斜率
ka
k
a
带入得
回波信号的包络为幅度归一化的矩形脉冲:
假定
t0=0→xp=0
t
0
=
0
→
x
p
=
0
,于是有
Sr(t)=a(t)ejπkat2
S
r
(
t
)
=
a
(
t
)
e
j
π
k
a
t
2
,其频谱为:
利用注定相位原理来计算上述积分。被积函数的相位为
φi=wt−πkat2
φ
i
=
w
t
−
π
k
a
t
2
φ,(t)
φ
,
(
t
)
表示回波信号的相位
πkat2
π
k
a
t
2
,驻定相位点的时刻
tk
t
k
为
±△
±
△
表示驻定相位点在
tk
t
k
附近的时刻。
把相位项
wt−πkat2
w
t
−
π
k
a
t
2
在驻定相位点
tk
t
k
展成幂级数,用
φ(t)
φ
(
t
)
表示回波信号的相位
πkat2
π
k
a
t
2
,有
在
t−tk
t
−
t
k
很小的条件下,取前三项即可。
那么,
做变量置换:
式中,积分项为菲涅尔积分的一种形式,有专用的数学表可查。若积分上限较大,则这一积分趋于
ejπ/2
e
j
π
/
2
。
化简后,由于归一化矩形脉冲的振幅a(t)为一常数,则振幅谱也是矩形的,可以写成:
回波信号的相位谱:
10.3.2 从频谱分析、相关、匹配滤波角度解释合成孔径原理
设地面上有两个点目标p1和p2,它们与飞行航向的垂直斜距相同,均为R0,二者所处方位不同,在x方向的坐标分别为x1、x2。根据上一节它们的回波信号都是线性调频信号,这两个线性调频信号的带宽都等于多普勒频移的带宽,其值为
式中,
Ts
T
s
表示合成孔径时间。p1和p2的区别:多普勒变化过程的起始点和终点不同。
设机载雷达天线从时间起点t=0和位置起点x=0开始向前移动,并发射第一个脉冲。天线波束在目标p1处的波束宽度,即合成孔径长度
Ls=θaR0
L
s
=
θ
a
R
0
。若这是波束前沿刚好照射到目标p1,p1回波信号多普勒频率变化的时间起始点t1=0。设第二个目标p2与p1的直线距离为
则p2的回波信号多普勒频率变化的起始点显然为
t2=△x/va
t
2
=△
x
/
v
a
。时间起始点两者的差别为
由于
va
v
a
是常数,如果能分辨△t,也就能分辨△x。从频域角度看:在同一时刻,两个回波信号的瞬时频率不太一样。设在时刻t飞机的位置在
xa=vat
x
a
=
v
a
t
,这时第一个回波的瞬时多普勒频率:
第二个回波的瞬时多普勒频率为:
两者的差为:
如果能够分辨这个频率差,也就能分辨
△x
△
x
。
分辨频率和分辨时间的途径:混频和相关
混频:采用一个具有相同调频斜率的线性调频信号作为本地信号,和两个回波信号
Sr1(t)
S
r
1
(
t
)
和
Sr2(t)
S
r
2
(
t
)
进行混频,即用一个乘法器进行相乘,在乘法器的输出端就会有和频及差频信号,在乘法器后面再接一个低通滤波器,则在滤波器输出端就会得到两个恒定频率(单频)的信号p1(f1)、p2(f2)。
将其振幅归一化后变为矩形振幅的单脉冲信号,脉宽为Ts(合成孔径时间)。
这类脉冲的频谱呈sinc型
其半功率点宽度均为:
△f=1Ts
△
f
=
1
T
s
,带宽与脉冲宽度成反比,即
即
设单个天线方位向孔径为
Dx
D
x
,则
相关技术
第一个回波的线性调频信号可写成,
Sr1(t)=ejπkat2
S
r
1
(
t
)
=
e
j
π
k
a
t
2
,其自相关函数可求出,
其输出自相关函数亦为sinc型,峰值点出现在
τ=0
τ
=
0
,其第一个零点发生在
πkaTsτ1=±π
π
k
a
T
s
τ
1
=
±
π
,即
主瓣宽度为:
半功率点一般取主瓣宽度的一半,即
第二个目标p2的回波相关函数的形状与第一个完全相同,知识峰值出现的时间稍晚,有一个时延
Δt
Δ
t
,
Δt
Δ
t
与目标p1,p2的距离
Δx
Δ
x
有如下关系
两个峰值可分辨的极限是:第二个峰值的延迟
±t
±
t
应不小于主瓣的半功率点宽度。
单个天线方位向孔径为 Dx D x ,则 Ls=λR0Dx→δt=Δt=Dx2va L s = λ R 0 D x → δ t = Δ t = D x 2 v a
于是得到空间分辨力为:
上述这种聚焦型处理中的自相关过程是等效于匹配滤波器的。
有匹配滤波理论:滤波器的冲激响应h(t)应是输入信号
si(t)
s
i
(
t
)
的延时镜像,即