快速POSP求解LFM的频谱
图中表面相位
ϕ
(
t
)
\phi(t)
ϕ(t)在
d
ϕ
(
t
)
/
d
t
d\phi(t)/dt
dϕ(t)/dt等于0的时刻
t
s
t_s
ts是驻留的,在上图中驻留点出现在
t
=
0
t=0
t=0 处。实虚部幅度变化一样,相位在该点的邻域附近是缓变的,而在其他时间点上是捷变的。当时宽带宽积(TBP)大于100时,POSP是相当准确的。图见黄皮书《合成孔径雷达成像算法与实现》42页。
设线性调频信号为
s
(
t
)
=
r
e
c
t
(
t
T
)
e
x
p
(
j
π
K
t
2
)
s(t)=\mathrm{rect}(\frac{t}T)\mathrm{exp}(j\pi Kt^2)
s(t)=rect(Tt)exp(jπKt2)
s
(
t
)
s(t)
s(t)的频谱为
G
(
f
)
=
∫
−
∞
+
∞
r
e
c
t
(
t
T
)
e
x
p
(
j
π
K
t
2
)
e
x
p
(
−
j
2
π
f
t
)
d
t
G(f)=\int_{-\infty}^{+\infty}\mathrm{rect}(\frac{t}T)\mathrm{exp}(j\pi Kt^2 )\mathrm{exp}(-j2\pi ft)dt
G(f)=∫−∞+∞rect(Tt)exp(jπKt2)exp(−j2πft)dt
其中包络为
w
(
t
)
w(t)
w(t)为矩形,被积相位为
θ
(
t
)
=
π
K
t
2
−
2
π
f
t
\theta(t)=\pi Kt^2-2\pi ft
θ(t)=πKt2−2πft
令
d
θ
d
t
=
0
\frac{d\theta}{dt}=0
dtdθ=0,则线性调频信号的时频关系为
f
=
K
t
,
t
=
f
K
f=Kt, t=\frac{f}K
f=Kt,t=Kf
频域相位为
Θ
(
f
)
=
θ
(
t
=
f
/
K
)
=
π
K
(
f
K
)
2
−
2
π
f
(
f
K
)
=
−
π
f
2
K
\Theta(f)=\theta(t=f/K)=\pi K(\frac{f}K)^2-2\pi f(\frac{f}K)=-\pi\frac{f^2}K
Θ(f)=θ(t=f/K)=πK(Kf)2−2πf(Kf)=−πKf2
频域包络为
W
(
f
)
=
w
(
t
=
f
/
K
)
=
r
e
c
t
(
f
∣
K
∣
T
)
W(f)=w(t=f/K)=\mathrm{rect}(\frac{f}{|K|T})
W(f)=w(t=f/K)=rect(∣K∣Tf)
其中矩形函数的中心频率为0.
直接给出信号频谱的近似
G
(
f
)
≈
C
1
W
(
f
)
e
x
p
{
j
(
Θ
(
f
)
±
π
/
4
)
}
G(f)\approx C_1W(f)\mathrm{exp}\{j(\Theta (f)\pm\pi/4)\}
G(f)≈C1W(f)exp{j(Θ(f)±π/4)}
将算得的
Θ
(
f
)
\Theta(f)
Θ(f)和
W
(
f
)
W(f)
W(f)代入
G
(
f
)
G(f)
G(f)中,并忽略
C
1
C_1
C1和相位
±
π
/
4
\pm\pi/4
±π/4,
可得到线性调频信号信号频谱为
G
(
f
)
=
r
e
c
t
(
f
K
T
)
e
x
p
{
−
j
π
f
2
K
}
G(f)=\mathrm{rect}(\frac{f}{KT})\mathrm{exp}\{-j\pi\frac{f^2}K\}
G(f)=rect(KTf)exp{−jπKf2}