维护英语-niel
6.5
动目标显示与动目标检测
引言
1.
目标回波频谱
6.5.1
目标回波和杂波的频谱
2.
杂波频谱
原理
递归
传统
非递归
6.5.2
MTI
滤波器
零点分配算法
滤波器设计
优化
预测误差算法
结语
原理
MTI+FFT
6.5.3
MTD
滤波器
滤波器设计
点最佳
等间隔最佳
结语
6.5.4
改善因子分析
MTI
MTD
6.5
动目标显示与动目标检测
雷达探测的运动目标如飞机,导弹,舰艇,车辆等周围存在各种背景,包括不动的地物和运动着的
云雨,海浪或金属丝干扰等。动目标
显示(
Moving
Target Indicator
:MTI
)与动目标检测(
Moving Target
Detection: MTD
)就是使用各种滤波器,滤去这些背景产生的杂波而取出
运动目标的回波。此外也可以通过把雷达安装在山上、增加雷达天线的倾角、安装防杂波
网来阻止杂波进入天线;或通过调整雷达天线的波
束形式、采用极化技术、降低雷达的分
辨单元、在时域采用
CFAR
检测、自适应门限、杂波图来抑制
杂波。在频域上应用
MTI
与
MTD<
/p>
技术
可以提高信杂比,
改善杂波背景下检
测运动目标的能力。
本节首先分析目标回波和杂波的频谱特性;
然后分别讨论
MTI
与
MTD
原理及滤波
器设计方法;最后分析
MTI
与
MTD
对改善因子的提高。
6.5.1
目标回波和杂波的频谱
运动目标回波
和杂波在频谱结构上有所差别,运动目标检测就是利用这种差别,
从频率上将它们区分,
以达到抑制杂波而显示目标回波
的目的。为此,应首先弄清楚目标和杂波的回波的特性。
(1)
目标回波的频谱
雷达发射相参脉冲串
,
其脉冲宽度为
T
e
< br>,
脉冲重复频率为
f
r
。当天线不扫描而对准目标时,所得脉冲为无限脉冲串。调制信号
u
p>
1
(
t
)
及
其频谱
U
1
(
f
)
分别为
?
t
p>
?
nT
r
u
1
(
t
)
?
A
?
rect
?
?
T
n
?
??
e
?
?
?
?
?
(6.5.1)
?
AT
e<
/p>
sin(
?
f
T
e
)
?
U
p>
1
(
f
)
?
?
?
(
f
?
n
f
< br>r
)
(6.5.2)
T
r
?
f
T
e
n
?
??
A
为信号振幅。而高频载波
u
p>
2
(
t
)
及其频谱
U
2
(
f
)
为
u
p>
2
(
t
)
?
cos
?
0
t
(6.5.3)
1
< br>U
2
(
f
)
?
[
?
(
f
?
f
0
p>
)
?
?
(
f
?
f
0
)]
(6.5.4)
< br>2
发射的相参脉冲串
u
(
t
)
?
u
1
(
t
)
?
u
2
(
t
)
,故其频谱
U
(
f
)
为
p>
U
(
f
)
?
U
1
(
f
)
?
U
< br>2
(
f
)
(6.5.5)
< br>发射信号频谱
f
r
脉冲发射
p>
A
t
回波
(
固定
)
f
0
?
f
0
?
1
?
(a)
固定回波频谱<
/p>
f
0
?
1
?
f
0
f
0
?
1
f
?
?
2R
0
< br>c
(a)
f
r
< br>t
回波
(
运动
< br>)
2(R
0
-
< br>?
R)
c
f
0
?
f
?
1
0
?
(b)
动目标回波频谱
f
0
?
1
?
f
0
f
(b)
f
r
t
u
m
(
t
)
?
(
f<
/p>
0
?
f
d
)
?
(
f
0
?
f
d
)
?
1
(c)
t
(d)
?
天线扫描动目标回
波频谱
U
m
(
f
)
(c)
f
0
?
f
d
1<
/p>
f
0
?
f
d
?
?
f
?
(
f
0
?
f
d
)
1
?
(
f
0
?
f
d
)<
/p>
?
?
(d)
f<
/p>
0
?
f
d
1
f
0
?
f
d
?
?
f
图
6.5.1
发射接收信号时域与频域特性
雷达发
射信号通常是窄带信号
(
如图
6.5.
1(a))
,因而运动目标回波频谱的特征是将发射信号的频谱位置在频率轴上平移一个
多普勒频率
f
d
?
2
v
r
?
(
如图
6.5.1(c))
,
f
d
的符号由目标运动的方向决定,靠近为正,
远离为负。固定回波的频谱与发射信号一致,只是幅度有衰减
(
如
图
6.5.1(b))
。
多普勒频率
f
d
可以直观地解释为:振荡源发射的电磁波以恒速
c
传播,当接收者相对振荡源不动,则它的接收频率等于发射频率。当
接收者与振荡源
之间有相对接近的运动时,则它接收电磁波的频率大于振荡源发射频率,当两者背向运动时,结果相反。运动目标
与雷达有
相对运动,所以存在多普勒频率
f
d
。定义固定回波的波程为
2
R<
/p>
0
,
R
0
为雷达到目标的距离,则回波到达雷达的时间为
标与雷达有相对运动,
其波程为
2(
R
0
?
?
R
)
,时间为
2
R
0
;由于运动目
c
2(
R
0
?
?
R
< br>)
。
c
雷达工作时,天线以各种方式进行扫描,这时收到的回波脉冲为有限数,且其振幅受天线方向图调制。设天线方
向图可用高斯函数来
表示,则收到的回波脉冲串的包络函数可写为
m
(
t
p>
)
?
2
?
?
ex
p
?
?
2
?
2
?
2
t
2
?
(6.5.6)
?
< br>是和天线波瓣宽度及扫描速度有关的参数。
?
减小,表示
观察时间增加。
天线扫描时收到的回波信号,可以用
m
(
t
)
和无限脉冲串
u
r
(
t
)
的乘积表示。
u
r
(
t
)
p>
为天线不扫描时的回波脉冲串,即
u
m
p>
(
t
)
?
m
(
t
)
u
r
(
t
< br>)
(6.5.7)
其包络函数
m
(
t
)
的频谱为
M
(
f
p>
)
?
e
?
f
2
2
?
2
(6.5.8)
m
(
t
)
M
(
f
)
?
s
0
t
0
f
图
6.5.2
天线扫描条件下回波频谱
天线扫描时
回波信号的频谱
U
m
(
f
)
为
p>
U
m
(
f
)
?
M
(
f
)
?
U
< br>r
(
f
)
(6.5.9)
即无限回波脉冲串
频谱
U
r
(
f
)
的每一根谱线均按
M
(
f
)
的形状展宽
(
如图
6.5.1(d))
。谱线展宽的程度反比于天线波束照射目标的时间
T
?
。已
求出当天线方向图为高斯形时谱线展宽的均方值为
?
s
?
p>
0
.
265
0
p>
.
265
f
r
p>
?
T
?
n
(6.5.10)
式中:
f
r
为雷达重复频率
,
n
为在单程天线方向图
3dB
宽度内收到的脉冲数。中频回波信号经过相位检波器后,相当于把中频信号的频谱
搬移到零频率附近,
根据目标多普勒频移
f
d
的不同,
相位检波后谱线
n
f
r
?
f
d
的具体位置也有差异,
每根谱线均
按脉冲串包络的频谱形状
展宽。
⑵
杂波频谱
对于固定点杂波,当天线不扫描时,固定杂波的频谱是位于
n<
/p>
f
r
上的谱线,可以用对消器全部滤去。
当天线扫描时,由于回波数目有
限,谱线将展宽。由于天线扫描引起双程天线方向图对回
波信号调幅,杂波谱展宽可用高斯函数表示为
p>
G
(
f
)
?
G
0
e
?
f
2
2
< br>2
?
s
(6.5.11)
其中
?
s
?
0
.
265
n
,
n
为在单程天线方向图
3d
B
宽度内的脉冲数。设
T
?
为天线照射目标的等效时间,则
n
?
T
?
f
r
< br>,即
?
s
?
0
.
265
T
?
,即
?
s
与目标照射时间成正比。
杂波信号的功率谱的实验公式可近似为
p>
W
(
f
)
?
g
(
f
)
?
g
0
< br>2
2
?
?
f
?
2
?
e
xp
?
?
a
?
?
f
?
?
p>
?
(6.5.12)
?
?
?
0
?
?
< br>?
其中:
W
(
< br>f
)
作为频率函数的杂波功率谱,
g
(
f
)
杂
波的傅立叶变换,
f
0
雷达载波,
p>
a
和杂波相关的参数。
< br>杂波频谱可以用杂波频率分布的均方根值
?
c
(
Hz
)
或速度分布的均
方根值
?
v
(
m
s
)
来表示,式(
< br>2-12
)可写为
?
f
2
?
?
f
2
?
2
< br>?
p>
W
(
f
)
?
W
0
exp
?
?
?
2
?
2
?
?
< br>?
W
0
exp
< br>?
?
?
8
?
2
?
?
(6.5.13)
c
?
v
?
?
?
< br>其中:
W
0
?
< br>g
0
2
c
c
2
2
,
?
c
?
,
?
p>
?
,
可得
a
?
。
?
c
为杂波功率谱方差。
?
v
为杂波内部起伏运动速度的均方根值,和工作波长无
2
f
0
?
8
?
v
2
?
v
关。相同的
?
v
值,对不同的雷达工作波长产生的杂波谱线的宽度也是不同的。工作波长越短,杂波谱的展宽越严重 。图
6.4.1
中的载波频率
f
0
为
1GHz
。
p>
(a)
(b)
0
(c)
f
0
f
0
(d)
f
0
f
0
?
f
p>
d
f
滤波器频率特性
0
f
d
f
f
0
?
f
d
p>
0
f
?
f
2
f
0
f
0
d
0
f
< br>d
f
图
6.5.3
杂波的功率谱
在接收机前端引入发射
信号作为基准电压,可得到收发频率的差频电压,即多普勒电压。图
6.5.4
给出各主要点的频谱图。
(a)
(b)
0
(c)
f
0
f
0
(d)
f
0
f
0
?
f
d
f
< br>滤波器频率特性
0
f
d
f
0
?
f
d
f
0
f
< br>0
?
f
d
2
f
0
f
0
f
d
f
图
p>
6.5.4
主要点频谱图及滤波器特性
C(f)
天线扫描造成
相位检波后的频
谱如图
(d),
固定杂波的频谱在零
于
天线扫描收到的回波脉冲数有限,谱线会有一定程
中还会出现杂乱分量,把它近似看成均
匀谱。本节要
图
6.5.4(d)
所示
,
取出动目标频谱,
滤除杂波频谱。
滤波器凹口和通带的平坦程度是关注的特性。动目标
?
f
0
点,图中给出的
是一条谱线。实际中,由
度的展宽,由于系统不稳等原因,杂波谱
系统不稳造成
设计的滤波器特性如
图
6.5.5
地杂波的功率谱
0
< br>f
0
f
滤波器要满足的要求:<
/p>
(1)
凹口适当扩宽,
与杂波梳状谱宽度
相当。
(2)
杂波有多普勒频移,即不在零频时,滤波器凹口要
对准杂波谱平均多普勒位置。
(3)
凹口深度能使杂波尽量多的
衰减,目标回波能在尽可能大的速度范围内有较大输出。
改善因子是综合评定滤波器性能的参数。
6.5.2
MTI
滤波器
(1)MTI
滤波器原理
当杂波和运动目标回波在雷达显示器上同时显示时,会使目标的观察变得很困难。如果目标处在强杂 波背景内,弱的目标淹没在强杂
波中,
特别是当强杂波使接收机
发生过载时,
将很难发现目标。
目标回波和杂波在时间域上难以
区分,
但由于目标的速度远大于背景的速度,
目标回波的多普勒
频移远大于背景的多普勒频移,从而可在频域上区分目标与杂波。动目标显示滤波器(
M
TI
)利用运动目标回波和杂波在
频谱上的区别,有效地抑制杂
波而提取信号。在雷达上加装
MTI
滤波器,大大的改善了雷达
在强杂波背景中检测运动目标的能力。
MTI
有
多种实现方法,包括传统的相消器和各种优化的
FIR
滤波器。采用重复参差和时变加权的
MTI
体制可以克服盲速。
MTI
和
M
TD
是对多个回波数据进行处理,
因此数据需要进行存储。
p>
其中每一行的数据是沿距离单元采样值,
反映了某一距离单元的信<
/p>
号特征。
每一列的数据为从相同的距离单元,
依次间隔一个脉冲重复周期的采样值,
这些数据的变化反映了在同一距离单元目标的
变化情况。
MTI
和
MTD
都是对同一距离单元的数据,即同一列的数据进行处理。通常
MTI
雷达滤波器如图
6.5.6
所示图,图中
p>
Tr
为雷达重复周期,在
这里作为延迟线的
延迟时间,
w
i
为滤波器权系数值。<
/p>
输入
Tr
Tr
Tr
Tr
w
0
w
1
w
2
p>
求和
输出
w
3
p>
…
w
N
?
1
图
6.5.6
MTI
滤波器的组成
传统的
MTI
相消器可以滤除零频杂波,性能不高,改善因子在
20dB
左右。优化的
MTI
滤波器可以满足上一节提出的动目标滤波器要
求,现在就滤波器的设计讨论如下。
(2)
MTI
滤波器设计
< br>1.
传统
MTI
滤波器设计:<
/p>
在相位检波器输出端,
固定目标的回波
是一串振幅不变的脉冲,
而运动目标的回波是一串振幅调制的脉冲。在把回波信号送到终
端显示
器前,必须先消除固定目标回波。最直接的方法是将相邻重复周期的回波信号相减
,则固定目标回波由于振幅不变而互相抵消,运动目标回
波相减后剩下相邻重复周期振幅
变化的部分。
传统的
MTI
滤波器有两种形式:非递归形和递归形。
(a)
非递归滤波器
不带反馈的滤波器称为
非递归型滤波器。下面以一次对消器为例进行说明。
一次对消
器,
即二脉冲对消。
其结构图如图
6.
5.7(a),
对消器的输入
X(z)
相位检波器的输出信号。
它是一个单零点系统,
零点位置在
p>
z
?
?
1
,
令
s
?
j
?
,即
z
?
e
j
?
T
在
Z
平面上是单位圆。
X
p>
(z)
Z
?
1
p>
?
?
+
Y
(z)
H
(
z
)
?
1
?
z
?
1
(a)
框图
(b)
零点图
图
6.5.7
一次对消滤波器
由相位检波器输出的脉冲包络为
u
(
t
p>
)
?
U
0
cos
?
(
t
)
(6.5.14)
?
为回波与基准电压之间的相位差
p>
?
(
t
)
?
?
?
0
t
r
?
?
< br>?
0
2
(
R
0
?
v
r
t
)
?
?
p>
d
t
?
?
0
(6.5.15)
c
回波信号按重复周期
T
r<
/p>
出现,将回波信号延迟一周期后,其包络是
u
'
(
p>
t
)
?
U
0
cos[
?
d
(
t
?
T
r
)
?
?
0
]
(6.5.16)
p>
?
T
?
?
T
?
?
?
?
u
?
u
< br>'
?
u
?
2
U
0
sin
?
d
r
?
s
in
?
?
d
?
d
r
?
?
p>
0
?
(6.5.17)
2
?
2
?
?
?
< br>输出包络为一多普勒频率的正弦信号,其频率为
2
p>
U
0
sin
?
p>
d
T
r
2
(6.5.18)
为多普勒频率的函数。当
< br>?
d
T
r
2
?
n
?
(
n=1,2,3,
…
)时,输出振幅为
零。这时的目标速度正相当于盲速,盲速是运动目标回波在相位检波
器的输出端与固定目
标回波相同,因而经对消设备后输出为零。
下面从频率域来说明对消器的工作原理。对消器的输出为
u
0
?
p>
u
i
(
1
?
e
?
j
?
T
r
)
< br>
(6.5.19)
对消器的频率响应特性为
j
p>
(
?
?
fT
r
)
u
(6.5.20)
H
(
j
?
)
< br>?
0
?
1
?
e
?
j
?
T
r
?
(
p>
1
?
cos
?
p>
T
r
)
?
j
sin
?
T
?
2
sin
?
f
T
r
e
2
u
i
?
< br>
图
6.5.8
对消器框图及其输出响应
其频率响应特性如图
p>
6.5.8
。
对消器等效于一个梳状滤波器
,
其频率特性在
f
?
< br>nf
r
各点处均为零。
固定目标
频谱的特点是它的谱线位于
nf
r
点<
/p>
上,因而在理想情况下,通过对消器后输出为零。当目标的多普勒频率为重复频率的整数倍
时,其频谱结构也有相同的特点,故通过上述梳
状滤波器后无输出。
表
6.5.1
二项式权系数和增益表
对消器阶数
1
2
3
4
5
增益
G
?
?
j
?
0
p>
n
?
1
w
j
(
线性
)
2
二项式系数
1, -1
1,
–
2, 1
1,
–
3, 3, -1
1,
–
4, 6,
–
4, 1
1, -5,10,
-10, 5, -1
2
2
6
20
70
252
采用级联二项式形的
N
次对消器(
N
?
1
脉冲对消)的权系数为
p>
?
1
?
x
?
展开式的系数,其计算式为
w
i
?
p>
?
?
1
?
i
?
1
N
!
,
i
?
< br>1
,
2
,...
N
?
1
(6.5.21)
?
N
?
i
?
1
?
!
?
i
?
1
?
!<
/p>
非递归滤波器中的权系数变化它的输出也会变化。
MTI
滤波器中的权系数可选二项式系数,也可选其它的系数。在不同的准则下,系
数也不同。在没有先验知识的情况下,二项式系数式一个很好的选择。
如果已知信号和杂波的协方差阵,在最优
MTI
滤
波器的准则是改善因子最大化时,则可以求其最优化的权系数。
如果杂波协方差阵为:
A
?
E
(
c
c
)
,信号协方差阵为:
S
?
E
(
s
s
)<
/p>
。其中:权系数向量为
w
?
(
N
?
1
< br>)
,信号向量为
s
?
(
N
?
1
)
,
杂波向量为
c
?
(
N
?
1
)
,
*
代表共轭,
T
代表转置。则
FIR
滤波器的最优化权系数为
p>
?
1
?
T
?
T
w
i
?
kA
?
1
s
*
(6.5.22)
式中
k
为任意常数,
A
为杂波
协方差阵的逆,
s
为目标回波信号向量。
最优化权系数是相对于一定准则的,准则不同结果也不同。如果没有先验知识或先验
知识不准确,则求得的最优权系数的实现效果可
能会小于二项式权系数实现的结果。
p>
(b)
递归滤波器
带有反馈的滤波器称为递
归型滤波器。递归型一次对消器的传递函数为
Y
(
p>
z
)
1
?
z
?
1
H
(
z
)
?
< br>?
X
(
z
)
1
?
K
1
z
?
1
(6.5.23)
递归滤波器的框图和频响如图
6.4.4
所示。虽然
递归滤波器有较好的频响,可以减小滤波器的阶数,但它的暂态响应差,限制了它在雷达信
号处理中的应用。
输入
?
1
w
1
?
2
w
< br>2
求和
?
3
w
3
?
N
?
1
…
w
N<
/p>
?
1
w
0
输出
(a)
递归滤波器框图
(b)
递归滤波器频率响应
图
6.5.9
波器框图及其输出响应
2.
优化
M
TI
滤波器设计:
在自适应杂波抑制
的应用环境下,
我们需要的理想滤波器是在杂波分布的频率点处,
能使杂波得到最大限度的抑制,而在其他频率点具
有最大平坦幅度。
< br>
(a)
零点分配算法
零点分配算法是在设计带阻滤波器时,在凹口处设置频率响应零点的一种方法,他不需要计算杂
波的相关系数。
对于
FIR
滤波器,如图
6.5.10
,设第
i
阶
输入
的延迟时间是
p>
?
i
。
FIR
p>
滤波器在
t
时刻的输出
是前几个时刻输入信号的加权和。
N
阶滤波器
w
0
?
1
< br>w
1
?
2
w
2
求和
?
3
w
3
?
N<
/p>
?
1
…
的权系数
为
w
i
,
i=
0
,
1
,
2<
/p>
,…,
N
,滤波器
w
N
?
1
图
6.5.10
FIR
滤波器框图
< br>的
传
输
函
数
为
:
输出
H<
/p>
?
f
?
?
?
w
i
e
?
j
2
?
fT
i
i
< br>?
0
N
(6.5.24)
式中,
T
0
=0
且
T
i
?
?
?
k
?
1
i
k
;
i
=0,1,2,
…
, N
。将
H
?
f
?
在
f
=
f
0
< br>处展开成泰勒级数:
H
?
p>
f
?
?
H
?
f
0
?
?
H
'
?
< br>f
0
?
?
(
f
?
f
0
)
?
H
''<
/p>
?
f
0
?
/
2
?
(
f
?
f
0
)
2
?
式中
< br>:
H
?
k
p>
?
(6.5.25)
?
< br>f
0
?
?
?
?
j
2
?
?
?
T
i
p>
k
w
i
e
?
j
2
?
i
?
0
k
< br>N
f
i
T
0
(6.5.26)
k
要在
f
=
f
0
处设计带阻滤波器,为使在
f
=
f
0
处附近幅度最大平坦,则必须
使泰勒级数展开式中
(
f
?
f
0
)
(
k
=0,1,2,
…
,
N
-1
)的系数为零,即
p>
H
?
k
?
?
f
0
?
?
0
,
k
=0<
/p>
,
1
,
2
,…,
N-1
(6.5.27)
这样,
就产生了
N
个关于
w
i
,的齐次线形方程。
?
T
i
?
0
N
k
i
w
i
e
?
j
2
?
f
0
T
i
?
0
k
=0,
1,2,
…
,N-1
(6.5.28)
其中,
T
0
0
?
1
,
w
0
为一个常数。将式
(6.5.28)
写成矩阵形式:
p>
A
?
W
?
?
w
0
U
(6.5.29)
维护英语-niel
维护英语-niel
维护英语-niel
维护英语-niel
维护英语-niel
维护英语-niel
维护英语-niel
维护英语-niel
-
上一篇:齿轮与蜗杆转动汇总
下一篇:《辽宁高职学报》论文格式要求与规范