线性调频信号数字脉冲压缩技术分析
2011年1月1日第34卷第1期
M odern Electro nics T echnique
现代电子技术
Jan. 2011V ol. 34N o. 1
线性调频信号数字脉冲压缩技术分析
郑力文, 孙晓乐
(中国空空导弹研究院, 河南洛阳 471009)
摘 要:在线性调频信号脉冲压缩原理的基础上, 利用M atlab 对数字脉冲压缩算法进行仿真, 得到了雷达目标回波信号经过脉冲压缩后的仿真结果。运用数字脉冲压缩处理中的中频采样技术与匹配滤波算法, 对中频采样滤波器进行了优化, 降低了实现复杂度, 减少了运算量与存储量。最后总结了匹配滤波的时域与频域实现方法, 得出在频域实现数字脉冲压缩方便, 运算量小, 更适合线性调频信号。
关键词:线性调频信号; 脉冲压缩; 中频采样; 匹配滤波
中图分类号:T N911 34 文献标识码:A 文章编号:1004 373X(2011) 01 0039 04
Digital Pulse C ompression Technology of Linear Frequency Modulation Signal
ZH ENG L i w en, SU N X iao le
(Chi na Airborne Missi le Academy, L uo yang 471009, China)
Abstract :Based o n the pr inciple of pulse com pr essio n techno lo gy o f linear fr equency mo dulat ion signal, the simulatio n r e sult of radar echo sig nal co mpressed by the pulse can be ga ined by using M atlab to simulate the dig ital pulse com pr essio n algo r ithm. Co mbining the techno log y o f IF sampling with the matching filt er alg or ithm in the digit al pulse compression processing and optimazing the I F sampling filter, which can remarkably reduce the complex ity and decr ease t he mult iplier operation and the memo ry. Finally , the implementation methods of matching filter algo rithm in time domain and fr equency doma in are summar ized, the dig ital pulse compression can be im plemented on frequency do main.
Keywords :linear frequency modulatio n signal; pulse com pr essio n; IF sampling ; matching f ilter
为了提高雷达系统的发现能力, 以及测量精度和分
辨能力, 要求雷达信号具有大的时宽带宽积[1 2]。但是, 在系统的发射和馈电设备峰值功率受限制的情况下, 大的信号能量只能通过加大信号的时宽来得到。然而单载频脉冲信号的时宽和带宽乘积接近1, 故大的时宽和带宽不可兼得。因此, 对这种信号来说, 测距精度和距离分辨力同测速精度和速度分辨力以及作用距离之间存在着不可调和的矛盾。在匹配滤波器理论的指导下, 提出了线性调频脉冲压缩的概念, 即在宽脉冲内附加线性调频, 以扩展信号的频带, 提供了一类信号, 其时宽带宽乘积大于1, 称之为脉冲压缩信号或大时宽带宽积信号。线性调频信号是应用最广泛的脉冲压缩信号, 因此线性调频信号的特性、脉冲压缩的原理及其实现技术都是比较受人关注的[3 5]。
1 线性调频信号脉冲压缩基本原理1. 1 线性调频信号简介
线性调频信号是通过非线性相位调制或线性频率调制(LFM ) 来获得大的时宽带宽积[6 7], 这种信号又称
为chirp 信号, 它是研究得最早而且应用最广泛的一种脉冲压缩信号。线性调频信号的时域波形如图1所示,
其频谱如图2所示。
线性调频信号可以表示为:
x (t) =A r ect exp j 2 f 0t +2
1, 0,
t/ 1/2t/
2
(1)
式中:A 为信号幅度; rect (t/ ) 为矩形函数, 即: rect (t/ ) =
(2)
线性调频信号的瞬时角频率 i 为:
i ==2 f 0+ t
d t
(3)
图1 线性调频信号的时域波形
在脉冲宽度 内, 信号的角频率由2 f 0- /2变
40
化到2 f 0+ /2, 调频带宽B = , 调频斜率为: =B/ 线性调频信号的时宽带宽积为:
D =B = 2
信号的复频谱为:X (f ) =
现代电子技术
2011年第34卷
速率与回波脉内速率的增加速度一致。于是就使得回
(4) (5)
波中的低频先到部分比高频后到部分通过滤波器的时
间滞后要长, 所以脉冲内的各频率分量在时域被积叠(或者说压缩) 在一起, 形成了幅度增大、宽度变窄的滤波器输出信号, 其理想包络如图3(e) 所示。
设回波信号的频率特性为:
X (f ) =|X (f ) |e j (f )
则匹配滤波器的频率特性H (f ) 应满足:
H (f ) =K X (f ) e
根据式(8) , 若令K =
(7)
-j (f )
在D 1的情况下
, 对式(6) 积分可得:X (f ) =
A 0,
j[-2 (f -f 0) /2 + /4]
e , f -f 0 B/2f -f 0>B/2
A
e
- /2- /2
+
t 2/2)
A rect e j(2 f 0t+ e -j2 f t d t =
j[2 (f -f )+ t /2]
2
(8)
-j2 ft
d0
d t (6)
e
(9)
, 则可得:
A
2/2 - /4-2 f
0d0
H (f ) =e j[2 (f -f 0)
特性) t d (f ) 为:
t ]
(10)
在式(10) 中, 匹配滤波器的群时延特性(频率延时=0+t d0, |f -f 0| (11) t d (f ) =d f B 2
式中:t d0是一个与滤波器物理实现有关的附加时延。
线性调频脉冲压缩匹配滤波器的输出信号为:
图2 线性调频信号的频谱
U(t) =
1. 2 脉冲压缩的原理
脉冲压缩的过程其实就是匹配滤波, 脉冲压缩的基本原理如图3所示。
A
- -
X (f ) H (f ) e j2 f t d f =
e j2 f 0(t-t d0) d f =d 0j2 f 0(t-t d0)
e
d0(12)
A 式(12) 为信号的复表示, 实际信号为实数, 故取其实部:U(t) =A
d 0
cos 2 f 0(t -t d0) (13)
B(t-t d 0)
故脉冲压缩匹配滤波器的输出包络为:
U 0(t) =A B(t -t d0) d0(14)
可以看出它具有辛克函数的形式, 如图4所示。
图3 脉冲压缩的基本原理示意图
图3(a) ~图3(c) 表示脉冲宽度为T 的线性调频信号, 也即回波信号。其中, (a) 为输入信号的波形; (b) 为输入信号的包络; (c) 为信号的载频调制特性; (d) 表示压缩滤波器(也即匹配滤波器) , 为压缩滤波器的延时频率特性; (e) 为压缩滤波器输出信号的包络。
如图3(a) 所示, 假定其载频在脉冲内按恒速(线性) 增加, 它通过脉冲压缩滤波器, 该滤波器具有如图3(d) 所示的时延频率特性, 即延时t d 随频率线性减小, 且减小
图4 线性调频信号的匹配滤波器输出
在图4中, 输出脉冲幅度下降到-4dB 处的脉冲宽度为T 0, 近似等于发射信号频谱宽度B 的倒数:T 0=1/B =T /D (15)
输出脉冲宽度T 0比输入脉冲宽度T 缩小了D 倍, 输出脉冲幅度为输入幅度A 的倍:
第1期
郑力文等:线性调频信号数字脉冲压缩技术分析
41
A 0=A (16)
在图4中, 主副瓣比约为13. 2dB (第一副瓣) , 第二副瓣再降低约4dB, 以后依次下降。太大的副瓣会影响对邻近弱目标的检测, 所以通常需要采取措施来降低副瓣。最常用的方法是窗函数加权, 加权函数可以选取海明加权, 余弦平方加权等。如果采用了窗函数加权, 副瓣电平将大大降低, 但同时也会使主瓣展宽和产生一定的信噪比损失, 信噪比损失约为1~1. 5dB 。图5
所示为加海明窗的脉冲压缩结果。
拉定理, 可以将其写为(e j 0t -e -j 0t ) /(2j ) , 显然, 该信号的频谱分布在 0和- 0处, 且符号相反, 如图6(a) 所示。对于带通信号情况也类似, 即实信号的谱总是共轭对称的两部分。
前面讲的采样定理都是针对实信号而言的, 对于复信号来说, 要进行I , Q 正交双路采样, 采样率可以降低50%, 仍以单频信号x (t) =cos 0t +jsin 0t 为例。实部的频谱如上所述, 虚部cos t =(e 0+e 0) /2在 0和- 0处有两个等幅相同的谱。将二者叠加就可得到一个单边谱, 即镜频抵消, 可以形象地用图6表示。图6(a) 和图6(b) 叠加后得到图6(c) , 镜频抵消, 正频能量增加1倍。
j t
-j t
图6 单频信号频谱
图5 加海明窗的脉冲压缩结果
把以上理论搬到对带限信号的分析, 也是如此。正交复信号的带宽相当于实信号带宽的50%, 因此, 对于实信号所要求的最低采样率f s =2B, 对复信号就可降低为f s =B 。另外, 由于信号的单边谱包含了原信号的所有信息, 故只需要单边谱就可以进行信号处理。
以常见的低通滤波法为例来说明正交采样的实现过程及原理, 如图7所示。
(17)
2 中频数字脉冲压缩信号处理实现技术2. 1 中频采样技术
将中频信号表示为:
X (t) =Re {a(t) ex p [j (2 f 0+ (t) ) ]}
=x (t) exp [2 f 0t+ (t) ]
式中:f 0为中频频率; a(t) 和x (t) 分别表示信号的幅度
和相位; x (t) =a(t) exp (j (t) ) =I +j Q, 称为X (t) 的复包络, 它包含了带通信号X (t) 的所有信息。令:
I =X I =a(t) co s (t) Q =X Q (t) =a(t) sin (t)
则:
a(t) =
X I (t) +X Q (t) =
图7 低通滤波法
在该例中, 根据带通信号采样定理选取f 0=B/2, f s =2B , 则A/D 变换后频谱示意如图6(b) 。之后将
(18)
信号分别乘以cos 0t n 和-sin 0t n , 由于t n =n/f s , 而
f 0=f s /4, 所以实际上cos 0t n 可以化简为cos n/2, 如图6所示, 其取值为1, 0, -1, 0的交替变换, -sin 0t n 也是一样。这样相乘以后相当于频谱中心左移 /2, 即将正频谱的中心移到了零频, 时域信号也相应分成了实部和虚部。然后通过低通滤波器, 频域上就是滤除掉等价于高频的镜频分量, 而在时域上是将分离出的实虚部通过滤波插值得到I , Q 两路同一时刻的采样值。滤波后数据率仍为2B, 进行1/2抽取是为了降低数据率, 在频域上等效为降低频谱间的间隔, 提高了频带的利用率, 最后得到的就是所需信号的复包络。
可见, 通过数字正交采样, 不仅得到了信号的复包
+Q (t) =tan -1[X I (t) /X Q (t) ]=tan -1(I /Q) (19) 式中:X I (t) 和X Q (t) 分别为X (t) 的同相分量和正交分量, 可分别用I 和Q 表示。在雷达信号处理中, X (t) 的同相分量和正交分量应保持式(18) 和式(19) 所表示的严格的幅度和相位关系。可以看出, 通过正交采样所得到的I , Q 两路正交信号, 可以很方便地得到信号的幅度和相位信息。
实信号的频谱是共轭对称的两部分, 通常将其中在频率正半轴的部分称为正频, 而将对称的在频率负半轴的部分称为负频或镜频。通过一个简单的例子就可以
, u(t) =0t,
42
现代电子技术
则输出y (n) 为:
2011年第34卷
络, 可以方便地获取信号的幅度及相位信息, 同时还降低了数据率, 大大减轻了对后续信号处理的运算负荷。2. 2 数字脉冲压缩
用数字技术实现脉冲压缩可采用时域方法和频域方法。一般而言, 对于大时宽带宽积信号, 用频域脉压较好; 对于小时宽带宽积信号, 用时域脉压较好2. 2. 1 时域卷积法
[8 10]
y (n) =IFFT (Y( ) ) =IFFT (S ( ) ) (24) 根据式(24) 可以画出频域快速卷积法实现脉冲压缩的原理框图, 如图9
所示。
2
。
时域脉冲压缩的过程是通过对接收信号s(t) 与匹配滤波器脉冲响应h(t) 求卷积的方法实现的。根据匹配滤波理论h(t) =s(t 0-t) , 即匹配滤波器是输入信号的共轭镜像, 并有相应的时移t 0。
用数字方法实现时, 输入离散信号为s(n) , 其匹配滤波器为h(n) , 匹配滤波器的输出为输入离散信号s(n) 与其匹配滤波器h(n) 的卷积:
N-1
图9 频域快速卷积法数字脉压原理图
图9是频域脉冲压缩的原理性框图, 在具体工程实现时常常将匹配滤波器权系数存于只读存储器中, 当需要加权以降低距离旁瓣时, 存于只读存储器中的权系数应当是匹配滤波器频率响应与加权函数的乘积。
在时域横向滤波器实现数字脉压, 对于N 点长度的输入信号, 需要进行N 次复数乘法, 而采用频域快速卷积法实现数字脉压, 需要进行N +N log 2N 次复数乘法。3 结 语
线性调频信号是应用广泛的雷达信号, 线性调频信号处理算法中比较重要的是中频采样技术与匹配滤波算法。本文从理论上分析了线性调频信号的特性以及脉冲压缩的基本原理, 并对正交变换技术的原理以及匹配滤波算法的两种实现方法进行了分析与比较。
参 考 文 献
[1]丁鹭飞, 耿富录. 雷达原理[M ]. 西安:西安电子科技大学出
版社, 2002.
[2]林茂庸, 柯有安. 雷达信号理论[M ].北京:国防工业出版
2
Y(n) =
s(k) *
k=0
N-1
h(n-k) =
k=0
h(k) *
s(n -k) (20)
式中:N 为信号采样点数。
按式(20) 构成的滤波器如图8所示。这是一种非递归的横向滤波器。由式(20) 可以排列出其他的计算方法, 图8所示的横向滤波器结构还有许多其他的等效网络结构。应当指出, 图8所示仅是原理性的, 在实际应用中往往在复数域进行滤波处理, 因此实际应用中应采用正交双通道滤波器, 按式(20)
完成复卷积运算。
图8 经典横向滤波器
社, 1984.
[3]马晓岩, 向家彬, 朱裕生, 等. 雷达信号处理[M ]. 长沙:湖南
科学技术出版社, 1998.
[4]杨小牛, 楼才义, 徐建良, 等. 软件无线电原理与应用[M ].
北京:电子工业出版社, 2001.
[5]车俐, 蒋留兵. 脉冲压缩在雷达中的应用及抗干扰性能研究
[J].甘肃科技, 2009(9) :49 51.
[6]孟泽, 孙合敏, 董礼, 等. 线性调频信号仿真及其特性分析
[J].舰船电子工程, 2009(8) :114 117.
[7]魏选平, 姚敏立, 张周生, 等. 脉冲压缩雷达原理及其M atlab
仿真[J].电子产品可靠性与环境试验, 2008(4) :36 38. [8]范红旗, 翟庆林, 王胜, 等. 基于IP 核的锐截止中频采样滤
波器优化设计[J].电子测量与仪器学报, 2008(1) :99 103. [9]赵树杰, 史林. 数字信号处理[M ]. 西安:西安电子科技大学
出版社, 1997.
[10]苏涛, 强生斌, 吴顺君, 等. 数字正交采样与脉压的高效算
法实现[J]. 现代雷达, 2001, 23(1) :39 48.
2. 2. 2 频域卷积实现方法
脉冲压缩过程是对输入信号s(n) 与匹配滤波器的脉冲响应h(n) 求卷积的过程。由傅里叶变化的性质可
知, 时域卷积相当于频域相乘。这个过程可以表示为:设输入离散信号为s(n) , 其傅里叶变换为S ( ) ; 匹配滤波器脉冲响应为h(n) , 其傅里叶变换为H ( ) ; 匹配滤波器输出为y (n) , 其傅里叶变换为Y( ) 。
对公式y (n) =s(n) *h(n) 两边同时进行傅里叶变换可得:
Y( ) =S( ) H ( )
又因为:
H ( ) =S ( )
代入式(21) 可得:
Y( ) =S( ) S *( )
*
(21) (22) (23)
作者简介:郑力文 男, 1981年出生, 河南郑州人, 硕士, 助理工程师。主要研究方向为雷达信号处理。
, 人,