CN103529429A - 用于线性调频信号的频域中的脉冲压缩方法 - Google Patents

用于线性调频信号的频域中的脉冲压缩方法 Download PDF

Info

Publication number
CN103529429A
CN103529429A CN201310523053.7A CN201310523053A CN103529429A CN 103529429 A CN103529429 A CN 103529429A CN 201310523053 A CN201310523053 A CN 201310523053A CN 103529429 A CN103529429 A CN 103529429A
Authority
CN
China
Prior art keywords
signal
linear
time domain
domain
frequency domain
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201310523053.7A
Other languages
English (en)
Other versions
CN103529429B (zh
Inventor
李秀华
刘李娟
张峰
李宏宇
高希权
李闯
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Institute of Radio Metrology and Measurement
Original Assignee
Beijing Institute of Radio Metrology and Measurement
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Institute of Radio Metrology and Measurement filed Critical Beijing Institute of Radio Metrology and Measurement
Priority to CN201310523053.7A priority Critical patent/CN103529429B/zh
Publication of CN103529429A publication Critical patent/CN103529429A/zh
Application granted granted Critical
Publication of CN103529429B publication Critical patent/CN103529429B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/02Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
    • G01S7/28Details of pulse systems
    • G01S7/285Receivers

Abstract

本发明公开了一种用于线性调频信号的频域中的脉冲压缩方法,包括如下步骤:对时域中的线性调频信号进行正交解调;对正交解调后的时域中的线性调频信号进行低通滤波;对低通滤波后的时域中的线性调频信号依次进行傅里叶变换和频带搬移得到频域中的线性调频信号;对时域中的匹配滤波信号依次进行傅里叶变换和频带搬移得到频域中的匹配滤波信号;以将时域中的匹配滤波信号转换为频域中的匹配滤波信号;将频域中的线性调频信号与频域中的匹配滤波信号进行复共轭相乘得到频域中的脉冲压缩信号;将频域中的脉冲压缩信号进行傅里叶逆变换得到时域中的脉冲压缩信号。本发明的脉冲压缩方法能够提高距离分辨力,降低信噪比损失;改善图像质量。

Description

用于线性调频信号的频域中的脉冲压缩方法
技术领域
本发明涉及脉冲压缩技术领域,特别涉及一种用于线性调频信号的频域中的脉冲压缩方法。
背景技术
随着遥感技术的发展,对雷达或引信的作用距离、分辨能力和测量精度等性能指标的要求也越来越高。如果要提高雷达系统的分辨能力和测量精度,要求雷达信号具有大的时宽带宽积。但是,由于雷达系统的发射设备和馈电设备的峰值功率的限制,只能通过加大脉宽的方式来获得能量较高的雷达信号。然而,单载频脉冲信号的脉宽即时宽与带宽的乘积接近1,故大的脉宽和大的带宽无法同时获得。
脉冲压缩技术,就是在雷达发射机中发射时宽较宽、峰值功率较低的脉冲,通过对载频进行编码以增加发射波形的带宽;然后在接收机中对回波实施脉冲压缩,从而获得时宽较窄、峰值功率较高的脉冲。
根据现有技术的匹配滤波器理论,用于线性调频信号的脉冲压缩是在宽脉冲内附加线性调频,以扩展信号的频带,使得信号的脉宽与带宽的乘积大于1。经脉冲压缩后的信号称为脉冲压缩信号或大时宽带宽积信号。
脉冲压缩技术能有效解决雷达作用距离与距离分辨力之间的矛盾,在不降低雷达作用距离的前提下提高雷达的距离分辨力,是实现雷达高分辨力,得到聚焦良好图像的有效途径。
现有技术中,用于线性调频信号的脉冲压缩方法存在如下缺点:距离分辨力较差,信噪比损失,图像质量较差。
发明内容
本发明的目的是针对现有技术的上述缺陷,提供一种用于线性调频信号的频域中的脉冲压缩方法。
本发明提供的用于线性调频信号的频域中的脉冲压缩方法包括如下步骤:
对时域中的线性调频信号进行正交解调;
对正交解调后的时域中的线性调频信号进行低通滤波;
对低通滤波后的时域中的线性调频信号依次进行傅里叶变换和频带搬移得到频域中的线性调频信号;
对时域中的匹配滤波信号依次进行傅里叶变换和频带搬移得到频域中的匹配滤波信号;以将时域中的匹配滤波信号转换为频域中的匹配滤波信号;
将频域中的线性调频信号与频域中的匹配滤波信号进行复共轭相乘得到频域中的脉冲压缩信号;
将频域中的脉冲压缩信号进行傅里叶逆变换得到时域中的脉冲压缩信号。
优选地,所述步骤“对时域中的线性调频信号进行正交解调”包括如下子步骤:
将时域中的线性调频信号与余弦信号相乘得到正交解调后的时域中的线性调频信号的实部分量;
将时域中的线性调频信号与正弦信号相乘得到正交解调后的时域中的线性调频信号的虚部分量。
进一步优选地,时域中的线性调频信号x(t)的表达式为:
x ( t ) = rect ( t T ) cos ( 2 π f 0 t + πKt 2 ) ;
其中,t为时间;T为脉冲宽度;f0为载频;K为线性调频率;
余弦信号的表达式为:
Rc=cos(2πf0t);
正弦信号的表达式为:
Rs=-sin(2πf0t)。
进一步优选地,正交解调后的时域中的线性调频信号的实部分量的表达式为:
x c ( t ) = 1 2 rect ( t T ) { cos [ πKt 2 ] + cos [ 4 π f 0 t + πKt 2 ] } ;
正交解调后的时域中的线性调频信号的虚部分量的表达式为:
x s ( t ) = 1 2 rect ( t T ) { sin [ πKt 2 ] - sin [ 4 π f 0 t + πKt 2 ] } .
进一步优选地,所述步骤“对正交解调后的时域中的线性调频信号进行低通滤波”包括如下子步骤:
对正交解调后的时域中的线性调频信号的实部分量进行低通滤波得到低通滤波后的时域中的线性调频信号的实部分量;
对正交解调后的时域中的线性调频信号的虚部分量进行低通滤波得到低通滤波后的时域中的线性调频信号的虚部分量。
进一步优选地,低通滤波后的时域中的线性调频信号的实部分量的表达式为:
x c a ( t ) = 1 2 rect ( t T ) cos [ πKt 2 ] ;
低通滤波后的时域中的线性调频信号的虚部分量的表达式为:
x s a ( t ) = 1 2 rect ( t T ) sin [ πKt 2 ] .
优选地,低通滤波后的时域中的线性调频信号的表达式为:
x a ( t ) = 1 2 rect ( t T ) exp { jπKt 2 } .
优选地,时域中的匹配滤波信号的表达式为:
h ( t ) = rect ( t T ) exp { - jπKt 2 } ;
其中,t为时间;T为脉冲宽度;K为线性调频率。
优选地,进行复共轭相乘的频域中的线性调频信号的数据长度与频域中的匹配滤波信号的数据长度相同。
优选地,时域中的脉冲压缩信号的表达式为:
sout(t)=Tsinc(KTt);
其中,t为时间;T为脉冲宽度;K为线性调频率。
本发明具有如下有益效果:
(1)与现有技术的脉冲压缩方法相比,本发明的脉冲压缩方法能够提高距离分辨力;
(2)与现有技术的脉冲压缩方法相比,本发明的脉冲压缩方法能够降低信噪比损失;
(3)与现有技术的脉冲压缩方法相比,本发明的脉冲压缩方法能够改善图像质量;
(4)本发明的脉冲压缩方法能够用于雷达和地震的探测系统的信号处理。
附图说明
图1为本发明实施例提供的用于线性调频信号的频域中的脉冲压缩方法的流程图。
具体实施方式
下面结合附图及实施例对本发明的发明内容作进一步的描述。
如图1所示,本实施例提供的用于线性调频信号的频域中的脉冲压缩方法包括如下步骤:
S1:对时域中的线性调频信号进行正交解调;
S2:对正交解调后的时域中的线性调频信号进行低通滤波;
S3:对低通滤波后的时域中的线性调频信号依次进行傅里叶变换和频带搬移得到频域中的线性调频信号;
S4:对时域中的匹配滤波信号依次进行傅里叶变换和频带搬移得到频域中的匹配滤波信号;
S5:将频域中的线性调频信号与频域中的匹配滤波信号进行复共轭相乘得到频域中的脉冲压缩信号;
S6:将频域中的脉冲压缩信号进行傅里叶逆变换得到时域中的脉冲压缩信号。
上述步骤S1进一步包括如下子步骤:
S1.1:将时域中的线性调频信号与余弦信号相乘得到正交解调后的时域中的线性调频信号的实部分量;
S1.2:将时域中的线性调频信号与正弦信号相乘得到正交解调后的时域中的线性调频信号的虚部分量。
需要说明的是,上述步骤S1.1与上述步骤S1.2的顺序可以互换。
上述步骤S1.1和步骤S1.2中,时域中的线性调频信号x(t)的表达式为:
x ( t ) = rect ( t T ) cos ( 2 π f 0 t + πKt 2 )    公式(1);
公式(1)中,t为时间;T为脉冲宽度;f0为载频;K为线性调频率;
上述步骤S1.1中,余弦信号的表达式为:
Rc=cos(2πf0t)    公式(2);
上述步骤S1.2中,正弦信号的表达式为:
Rs=-sin(2πf0t)    公式(3)。
上述步骤S1.1得到的正交解调后的时域中的线性调频信号x(t)的实部分量xc(t)的表达式为:
x c ( t ) = 1 2 rect ( t T ) { cos [ πKt 2 ] + cos [ 4 π f 0 t + πKt 2 ] }    公式(4);
上述步骤S1.2得到的正交解调后的时域中的线性调频信号x(t)的虚部分量xs(t)的表达式为:
x s ( t ) = 1 2 rect ( t T ) { sin [ πKt 2 ] - sin [ 4 π f 0 t + πKt 2 ] }    公式(5)。
上述步骤S2进一步包括如下子步骤:
S2.1:对正交解调后的时域中的线性调频信号的实部分量进行低通滤波得到低通滤波后的时域中的线性调频信号的实部分量;
S2.2:对正交解调后的时域中的线性调频信号的虚部分量进行低通滤波得到低通滤波后的时域中的线性调频信号的虚部分量。
需要说明的是,上述步骤S2.1与上述步骤S2.2的顺序可以互换。
上述步骤S2.1得到的低通滤波后的时域中的线性调频信号x(t)的实部分量
Figure BDA0000404440010000056
的表达式为:
x c a ( t ) = 1 2 rect ( t T ) cos [ πKt 2 ]    公式(6);
上述步骤S2.2得到的低通滤波后的时域中的线性调频信号x(t)的虚部分量
Figure BDA0000404440010000055
的表达式为:
x s a ( t ) = 1 2 rect ( t T ) sin [ πKt 2 ]    公式(7)。
根据公式(6)和公式(7),低通滤波后的时域中的线性调频信号的表达式为:
x a ( t ) = 1 2 rect ( t T ) exp { jπKt 2 }    公式(8)。
上述步骤S4中,时域中的匹配滤波信号h(t)的表达式为:
h ( t ) = rect ( t T ) exp { - jπKt 2 }    公式(9);
公式(9)中,t为时间;T为脉冲宽度;K为线性调频率。
优选地,上述步骤S5中,进行复共轭相乘的频域中的线性调频信号的数据长度与频域中的匹配滤波信号的数据长度相同。
上述步骤S6得到的时域中的脉冲压缩信号sout(t)的表达式为:
sout(t)=Tsinc(KTt)    公式(10);
公式(10)中,t为时间;T为脉冲宽度;K为线性调频率。
应当理解,以上借助优选实施例对本发明的技术方案进行的详细说明是示意性的而非限制性的。本领域的普通技术人员在阅读本发明说明书的基础上可以对各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (10)

1.用于线性调频信号的频域中的脉冲压缩方法,其特征在于,该脉冲压缩方法包括如下步骤:
对时域中的线性调频信号进行正交解调;
对正交解调后的时域中的线性调频信号进行低通滤波;
对低通滤波后的时域中的线性调频信号依次进行傅里叶变换和频带搬移得到频域中的线性调频信号;
对时域中的匹配滤波信号依次进行傅里叶变换和频带搬移得到频域中的匹配滤波信号;以将时域中的匹配滤波信号转换为频域中的匹配滤波信号;
将频域中的线性调频信号与频域中的匹配滤波信号进行复共轭相乘得到频域中的脉冲压缩信号;
将频域中的脉冲压缩信号进行傅里叶逆变换得到时域中的脉冲压缩信号。
2.根据权利要求1所述的用于线性调频信号的频域中的脉冲压缩方法,其特征在于,所述步骤“对时域中的线性调频信号进行正交解调”包括如下子步骤:
将时域中的线性调频信号与余弦信号相乘得到正交解调后的时域中的线性调频信号的实部分量;
将时域中的线性调频信号与正弦信号相乘得到正交解调后的时域中的线性调频信号的虚部分量。
3.根据权利要求2所述的用于线性调频信号的频域中的脉冲压缩方法,其特征在于,
时域中的线性调频信号x(t)的表达式为:
x ( t ) = rect ( t T ) cos ( 2 π f 0 t + πKt 2 ) ;
其中,t为时间;T为脉冲宽度;f0为载频;K为线性调频率;
余弦信号的表达式为:
Rc=cos(2πf0t);
正弦信号的表达式为:
Rs=-sin(2πf0t)。
4.根据权利要求2所述的用于线性调频信号的频域中的脉冲压缩方法,其特征在于,
正交解调后的时域中的线性调频信号的实部分量的表达式为:
x c ( t ) = 1 2 rect ( t T ) { cos [ πKt 2 ] + cos [ 4 π f 0 t + πKt 2 ] } ;
正交解调后的时域中的线性调频信号的虚部分量的表达式为:
x s ( t ) = 1 2 rect ( t T ) { sin [ πKt 2 ] - sin [ 4 π f 0 t + πKt 2 ] } .
5.根据权利要求1所述的用于线性调频信号的频域中的脉冲压缩方法,其特征在于,所述步骤“对正交解调后的时域中的线性调频信号进行低通滤波”包括如下子步骤:
对正交解调后的时域中的线性调频信号的实部分量进行低通滤波得到低通滤波后的时域中的线性调频信号的实部分量;
对正交解调后的时域中的线性调频信号的虚部分量进行低通滤波得到低通滤波后的时域中的线性调频信号的虚部分量。
6.根据权利要求5所述的用于线性调频信号的频域中的脉冲压缩方法,其特征在于,
低通滤波后的时域中的线性调频信号的实部分量的表达式为:
x c a ( t ) = 1 2 rect ( t T ) cos [ πKt 2 ] ;
低通滤波后的时域中的线性调频信号的虚部分量的表达式为:
x s a ( t ) = 1 2 rect ( t T ) sin [ πKt 2 ] .
7.根据权利要求1所述的用于线性调频信号的频域中的脉冲压缩方法,其特征在于,低通滤波后的时域中的线性调频信号的表达式为:
x a ( t ) = 1 2 rect ( t T ) exp { jπKt 2 } .
8.根据权利要求1所述的用于线性调频信号的频域中的脉冲压缩方法,其特征在于,时域中的匹配滤波信号的表达式为:
h ( t ) = rect ( t T ) exp { - jπKt 2 } ;
其中,t为时间;T为脉冲宽度;K为线性调频率。
9.根据权利要求1所述的用于线性调频信号的频域中的脉冲压缩方法,其特征在于,进行复共轭相乘的频域中的线性调频信号的数据长度与频域中的匹配滤波信号的数据长度相同。
10.根据权利要求1所述的用于线性调频信号的频域中的脉冲压缩方法,其特征在于,时域中的脉冲压缩信号的表达式为:
sout(t)=Tsinc(KTt);
其中,t为时间;T为脉冲宽度;K为线性调频率。
CN201310523053.7A 2013-10-29 2013-10-29 用于线性调频信号的频域中的脉冲压缩方法 Active CN103529429B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310523053.7A CN103529429B (zh) 2013-10-29 2013-10-29 用于线性调频信号的频域中的脉冲压缩方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310523053.7A CN103529429B (zh) 2013-10-29 2013-10-29 用于线性调频信号的频域中的脉冲压缩方法

Publications (2)

Publication Number Publication Date
CN103529429A true CN103529429A (zh) 2014-01-22
CN103529429B CN103529429B (zh) 2016-08-17

Family

ID=49931585

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310523053.7A Active CN103529429B (zh) 2013-10-29 2013-10-29 用于线性调频信号的频域中的脉冲压缩方法

Country Status (1)

Country Link
CN (1) CN103529429B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104168033A (zh) * 2014-07-11 2014-11-26 天津理工大学 一种时宽带宽积大于10000的脉冲压缩信号产生的方法
CN104811163A (zh) * 2015-05-20 2015-07-29 天津理工大学 一种时宽带宽积超过10000的周期信号产生的方法
CN106908768A (zh) * 2017-02-24 2017-06-30 四川九洲空管科技有限责任公司 一种超低旁瓣的脉冲压缩方法
CN108886509A (zh) * 2016-02-16 2018-11-23 诺基亚通信公司 压缩/解压缩频域信号
CN109490852A (zh) * 2017-09-10 2019-03-19 北京遥感设备研究所 一种线性调频脉冲信号调频斜率极性确定方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101424738A (zh) * 2008-11-13 2009-05-06 国营险峰机器厂 基于fpga的实时数字脉冲压缩系统的处理方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101424738A (zh) * 2008-11-13 2009-05-06 国营险峰机器厂 基于fpga的实时数字脉冲压缩系统的处理方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
邵智超 等: "基于matlab的雷达信号处理仿真系统", 《计算机仿真》, vol. 24, no. 6, 30 June 2007 (2007-06-30), pages 268 - 271 *
陈岳明: "线性调频脉冲压缩相关技术研究", 《中国优秀硕士学位论文全文数据库(电子期刊)》, no. 5, 31 May 2010 (2010-05-31), pages 8 *
高雪迎: "线性调频脉冲压缩系统的设计与研究", 《西安电子科技大学硕士学位论文》, 29 September 2012 (2012-09-29), pages 25 - 44 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104168033A (zh) * 2014-07-11 2014-11-26 天津理工大学 一种时宽带宽积大于10000的脉冲压缩信号产生的方法
CN104811163A (zh) * 2015-05-20 2015-07-29 天津理工大学 一种时宽带宽积超过10000的周期信号产生的方法
CN104811163B (zh) * 2015-05-20 2017-05-17 天津理工大学 一种时宽带宽积超过10000的周期信号产生的方法
CN108886509A (zh) * 2016-02-16 2018-11-23 诺基亚通信公司 压缩/解压缩频域信号
CN106908768A (zh) * 2017-02-24 2017-06-30 四川九洲空管科技有限责任公司 一种超低旁瓣的脉冲压缩方法
CN109490852A (zh) * 2017-09-10 2019-03-19 北京遥感设备研究所 一种线性调频脉冲信号调频斜率极性确定方法

Also Published As

Publication number Publication date
CN103529429B (zh) 2016-08-17

Similar Documents

Publication Publication Date Title
CN103529429A (zh) 用于线性调频信号的频域中的脉冲压缩方法
CN102508200B (zh) 复数扩频信号的到达时间差测量方法
CN100559210C (zh) 移动通信环境中的速度估计装置
CN102998673B (zh) 一种用于合成孔径雷达的压缩感知成像方法
CN109510787A (zh) 水声信道下线性调频信号参数估计方法及装置
WO2015070820A1 (zh) 扩频信号的生成方法、生成装置、接收方法和接收装置
CN104035079B (zh) 一种新型零中频脉冲压缩雷达系统的信号性能补偿方法
CN110071738A (zh) 基于多伪码支路接收机的扩频信号解扩跟踪方法
CN106324598B (zh) 一种多通道子带信号的均衡、同步方法及其系统
CN106842198A (zh) 一种基于凸优化理论的子带融合成像算法
CN110247866A (zh) 一种基于doa估计的dmwc频谱感知相位校准方法
CN102655491A (zh) 相干解调频移键控调制信号的频率偏移估计方法及系统
CN105891600A (zh) 四相相移键控信号频谱估计方法
CN110113273B (zh) 一种基于频谱对称的时频重叠信号盲分离方法
CN106341123B (zh) 一种单音干扰的滤波方法和装置
CN106568421B (zh) 深海走航投弃式剖面仪长距离线性时变传输信道建模的方法
CN101447958B (zh) 一种调制器性能的测量方法和装置
CN105068061B (zh) 基于chirp混频的模拟信息转换方法及系统
CN102780666A (zh) 一种pswf脉冲调制信号的抗定时抖动方法
CN106301464B (zh) 一种针对chirp信号的信噪比估计方法
CN103648160B (zh) 一种信号测距和定位方法
CN106646348A (zh) 适用于短时多频率信号的干涉仪相位差测量电路及测量方法
CN103529430B (zh) 一种用于线性调频信号的匹配压缩装置及匹配压缩方法
CN104199003A (zh) 基于双线性变换的超宽带线性调频信号采样方法
CN201035133Y (zh) 低数据流量到达时间差定位装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant