CN101776747B - 基于自适应重采样的雷达信号压缩方法 - Google Patents

基于自适应重采样的雷达信号压缩方法 Download PDF

Info

Publication number
CN101776747B
CN101776747B CN2010100118791A CN201010011879A CN101776747B CN 101776747 B CN101776747 B CN 101776747B CN 2010100118791 A CN2010100118791 A CN 2010100118791A CN 201010011879 A CN201010011879 A CN 201010011879A CN 101776747 B CN101776747 B CN 101776747B
Authority
CN
China
Prior art keywords
signal
radar signal
compression
radar
resampling
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.)
Expired - Fee Related
Application number
CN2010100118791A
Other languages
English (en)
Other versions
CN101776747A (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.)
Shandong University
Original Assignee
Shandong University
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 Shandong University filed Critical Shandong University
Priority to CN2010100118791A priority Critical patent/CN101776747B/zh
Publication of CN101776747A publication Critical patent/CN101776747A/zh
Application granted granted Critical
Publication of CN101776747B publication Critical patent/CN101776747B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明提供了一种基于自适应重采样的雷达信号压缩方法,具体包括以下步骤:(1)对雷达信号进行分段,根据对雷达信号压缩比和处理速度的要求合理确定每段雷达信号的长度;(2)对每段雷达信号进行频谱分析,通过分析确定每段雷达信号的重采样率,并确定雷达信号的压缩比;(3)根据新的采样频率对雷达信号重采样,实现对雷达信号的压缩;(4)在需要的情况下对雷达信号进行重构,首先对雷达信号进行零插值,将压缩后的信号恢复到原始长度,并通过截止频率为信号重采样频率一半的低通滤波器滤波实现信号的重构。本发明通过对不同时段的信号分别进行压缩,降低了采样频率,解决了雷达信号采集时采样频率过高的问题。

Description

基于自适应重采样的雷达信号压缩方法
技术领域
本发明涉及一种用于雷达辐射源信号识别、分析等处理的雷达信号压缩方法。
技术背景
雷达已广泛应用于军事和民用的各个领域。为了获得有效信息,需要对采集到的雷达信号进行识别和分析。但是在实际应用中,在对雷达信号进行识别和分析之前,通常需要先对雷达信号进行存储或者传输。而现实应用中采集到的雷达信号数据量极大,这给数据实时存储和传输带来很大的困难。在同样的通信容量或者存储容量下,如果先对雷达信号进行压缩,然后再进行存储或传输,就可以提高信息的存储和传输能力,因此,对雷达信号进行有效的压缩是十分有意义的。
随着雷达技术的发展,雷达信号的种类越来越多,复杂度越来越高。不同种类不同参数的雷达信号对采样率的要求是不一样的,但是雷达信号在采集时通常不能根据信号特点灵活变化采样频率,而且为了保证信号信息的完整性,雷达信号采集时会采用很高的采样频率。
发明内容
本发明针对雷达信号由于采集时采样频率远高于所需采样频率而造成的数据量大的问题,提供一种通过对不同时段的信号分别进行压缩而降低采样频率的基于自适应重采样的雷达信号压缩方法。
本发明的基于自适应重采样的雷达信号压缩方法,是在雷达信号采集设备采集到雷达信号后,根据信号特点及实际需要对雷达信号分段,针对每段信号的频谱特点确定其重采样频率,然后通过对雷达信号重新采样实现对信号的压缩,具体包括以下步骤:
(1)对雷达信号进行分段
根据对雷达信号压缩比和处理速度的要求合理确定每段雷达信号的长度,根据雷达信号特点和处理设备的要求进行调整,分段加长会提高处理速度,但是分段增长后会因为频率分量丰富而造成压缩效果的降低,同时,分段长度也会受到硬件处理设备处理能力的限制;而分段减短后压缩比会提高但是会因为计算量的增加而影响处理速度;
在本发明中,每段雷达信号的长度是以每10ms采集的数据为一个分段。
(2)对每段雷达信号进行频谱分析,通过分析确定每段雷达信号的重采样率,并确定雷达信号的压缩比;计算雷达信号的奈奎斯特率(即最小采样频率),并根据雷达信号的实际情况确定其压缩比M;
在计算奈奎斯特率时,将雷达信号的实际采样率Fs归一化,然后计算其奈奎斯特率Fn,得到Fn后,计算Fs与Fn的比值M,即:M=Fs/Fn,M即为雷达信号的压缩比;
具体计算方法如下:
对信号x[n]作N点FFT,得到离散傅里叶变换序列X[k],对|X[k]|2做中值滤波,选定阈值Tf,以|X[k]|2中所有大于Tf的值之中距离零频率最远的数值所对应的频率作为信号的最高频率分量fm,从而得到Fn=2fm,其中将阈值Tf确定为|X[k]|2均值,根据fm可以得到雷达信号压缩比M的值为:
Figure G2010100118791D00021
即M取Fs/(2fm)结果的整数部分。此时,信号的重采样率为:Fn=Fs/M。
(3)根据新的采样频率对雷达信号重采样,实现对雷达信号的压缩;
压缩时,为了方便信息处理和信号重构,在传输和存储的编码过程中需要对每段信号添加附加信息,用来保留信号的压缩比,以便在处理时获取信号的重采样频率(信号接收方已知信号原始采样频率),或者在重构时用于恢复信号;
(4)在需要的情况下对雷达信号进行重构,首先对雷达信号进行零插值,将压缩后的信号恢复到原始长度,并通过截止频率为信号重采样频率一半的低通滤波器滤波实现信号的重构。
本发明通过对不同时段的信号分别进行压缩,降低了采样频率,解决了雷达信号采集时采样频率过高的问题。
附图说明
图1是本发明雷达信号压缩模型示意图。
图2是本发明雷达信号重构模型示意图。
具体实施方式
在雷达接收机(或接收端)采集到雷达信号后,将得到的数据信号按本发明的方法压缩,详细步骤如下:
(1)对雷达信号进行分段
根据对雷达信号压缩比和处理速度的要求合理确定每段雷达信号的长度;在本发明中,以每10ms采集的数据为一个分段。由于雷达信号采样率的不同需要对雷达信号进行分段,分段的长度需要根据实际需要灵活确定,分段过长会因为频率分量过于丰富造成无法到达理想的压缩效果,而分段过短会因为计算量的增加而影响处理速度。雷达信号压缩模型如图1所示。
(2)雷达信号的频谱估计和采样率的确定
对每段雷达信号进行频谱分析,通过分析确定其重采样率,并确定雷达信号的压缩比。本发明中对雷达信号的压缩是在实际采样频率高于奈奎斯特率(即最小采样频率)的情况,最后取得的压缩效果也与采样频率的冗余度有关,因此需要计算雷达信号的奈奎斯特率,并根据雷达信号的实际情况确定其压缩比M。
在计算奈奎斯特率时,将信号的实际采样率Fs归一化,然后计算其奈奎斯特率Fn,得到Fn后,计算Fs与Fn的比值M,即M=Fs/Fn,M即为雷达信号的压缩比。
具体计算方法如下:
对信号x[n]作N点FFT,得到离散傅里叶变换序列X[k],对|X[k]|2做中值滤波;
选定阈值Tf,以|X[k]|2中所有大于Tf的值之中距离零频率最远的数值所对应的频率作为信号的最高频率分量fm,从而得到Fn=2fm,其中将阈值Tf确定为|X[k]|2均值;
根据fm可以得到雷达信号压缩比M的值为:
即M取Fs/(2fm)结果的整数部分。此时,信号的重采样率为:Fn=Fs/M。
(3)通过重采样实现对每段雷达信号进行压缩
压缩时,为了方便信息处理和信号重构,在传输和存储的编码过程中需要对每段信号添加附加信息,用来保留信号的压缩比,以便在处理时获取信号的重采样频率(信号接收方已知信号原始采样频率),或者在重构时用于恢复信号。
确定了雷达信号的压缩比后,通过对信号的抽取实现对信号的重采样和压缩;假设原始雷达信号为x[n],压缩后得到的信号为:xd[n]=Mx[nM];
这样就实现了对信号的压缩。之所以将抽取得到的序列增加M倍是为了保证信号频谱幅度保持不变。由于在重采样率的计算过程中保证了信号的频谱不会发生混叠,因此压缩后的信号完整保留了信号的频谱信息。
(4)在需要的情况下对雷达信号进行重构,首先对雷达信号进行零插值,将压缩后的信号恢复到原始长度,并通过截止频率为信号重采样频率一半的低通滤波器滤波实现信号的重构。
在本发明的雷达信号压缩过程中,由于保留了信号的频谱信息,因此在对雷达信号的分析和识别过程中可以直接对压缩后的信号进行处理,但有时也会需要对雷达信号进行还原,因此本发明设计了雷达信号的重构模型,如图2所示。
首先对雷达信号进行零插值,将压缩后的信号xd[n]恢复到原始长度,得到
将xe[n]通过一个低通滤波器进行滤波,就实现了对雷达信号的重构,其中低通滤波器的增益为1,截止频率为Fs/2M,即重采样率的一半。

Claims (3)

1.一种基于自适应重采样的雷达信号压缩方法,其特征是,在雷达信号采集设备采集到雷达信号后,根据信号特点及实际需要对雷达信号分段,针对每段信号的频谱特点确定其重采样频率,然后通过对雷达信号重新采样实现对信号的压缩,具体包括以下步骤:
(1)对雷达信号进行分段,根据对雷达信号压缩比和处理速度的要求合理确定每段雷达信号的长度;
(2)对每段雷达信号进行频谱分析,通过分析确定每段雷达信号的重采样率,并确定雷达信号的压缩比;计算雷达信号的奈奎斯特率,并根据雷达信号的实际情况确定其压缩比M;在计算奈奎斯特率时,将雷达信号的实际采样率Fs归一化,然后计算其奈奎斯特率Fn,得到Fn后,计算Fs与Fn的比值M,即:M=Fs/Fn,M即为雷达信号的压缩比;
(3)根据新的采样频率对雷达信号重采样,实现对雷达信号的压缩;
压缩时,为了方便信息处理和信号重构,在传输和存储的编码过程中需要对每段信号添加附加信息,用来保留信号的压缩比,以便在处理时获取信号的重采样频率,或者在重构时用于恢复信号;
(4)在需要的情况下对雷达信号进行重构,首先对雷达信号进行零插值,将压缩后的信号恢复到原始长度,并通过截止频率为信号重采样频率一半的低通滤波器滤波实现信号的重构。
2.根据权利要求1所述的基于自适应重采样的雷达信号压缩方法,其特征是,所述步骤(1)中的每段雷达信号的长度是以每10ms采集的数据为一个分段。
3.根据权利要求1所述的基于自适应重采样的雷达信号压缩方法,其特征是,所述步骤(2)中雷达信号压缩比的具体计算方法如下:
对信号x[n]作N点FFT,得到离散傅里叶变换序列X[k],对|X[k]|2做中值滤波,选定阈值Tf,以|X[k]|2中所有大于Tf的值之中距离零频率最远的数值所对应的频率作为信号的最高频率分量fm,从而得到Fn=2fm,其中将阈值Tf确定为|X[k]|2均值,根据fm可以得到雷达信号压缩比M的值为:
Figure FSB00000763343700011
即M取Fs/(2fm)结果的整数部分,此时,信号的重采样率为:Fn=Fs/M。
CN2010100118791A 2010-01-15 2010-01-15 基于自适应重采样的雷达信号压缩方法 Expired - Fee Related CN101776747B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010100118791A CN101776747B (zh) 2010-01-15 2010-01-15 基于自适应重采样的雷达信号压缩方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010100118791A CN101776747B (zh) 2010-01-15 2010-01-15 基于自适应重采样的雷达信号压缩方法

Publications (2)

Publication Number Publication Date
CN101776747A CN101776747A (zh) 2010-07-14
CN101776747B true CN101776747B (zh) 2012-07-18

Family

ID=42513255

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010100118791A Expired - Fee Related CN101776747B (zh) 2010-01-15 2010-01-15 基于自适应重采样的雷达信号压缩方法

Country Status (1)

Country Link
CN (1) CN101776747B (zh)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103093608B (zh) * 2011-10-31 2014-12-24 北京掌城科技有限公司 高频gps数据的过滤方法及系统
CN103518356A (zh) * 2013-03-07 2014-01-15 华为技术有限公司 数据压缩方法、还原方法,装置及系统
CN103475376A (zh) * 2013-09-24 2013-12-25 上海无线电设备研究所 一种雷达多普勒回波信号数据压缩方法
CN105279193B (zh) * 2014-07-22 2020-05-01 腾讯科技(深圳)有限公司 文件处理方法及装置
CN107870330B (zh) * 2017-09-27 2019-11-29 西安空间无线电技术研究所 分段变重频模式采样均匀化数字处理方法及装置
CN108362939B (zh) * 2018-01-31 2020-06-23 成都泰格微波技术股份有限公司 一种线性调频信号的频域参数测量方法
CN110932753B (zh) * 2019-12-11 2021-04-09 电子科技大学 一种基于智能决策的变换域自适应通信传输方法
CN111323794B (zh) * 2020-03-21 2023-05-05 哈尔滨工程大学 一种基于周期重采样的周期调频干扰消除方法
CN116660612A (zh) * 2023-07-31 2023-08-29 青岛鼎信通讯科技有限公司 一种特征电流检测方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4053840A (en) * 1975-08-01 1977-10-11 Thomson-Csf Band compression device
US4673941A (en) * 1983-12-07 1987-06-16 Hollandse Signaalapparaten B.V. Digital pulse compression filter
US5850202A (en) * 1997-07-29 1998-12-15 Erim International, Inc. Method for improving SAR system sensitivity in the presence of RF interference
CN1664860A (zh) * 2004-03-05 2005-09-07 中国科学院计算技术研究所 合成孔径雷达复数图像数据的实时自动压缩方法
CN1996042A (zh) * 2006-12-31 2007-07-11 武汉大学 大时间带宽积线性调频信号的数字脉冲压缩方法
CN101165510A (zh) * 2006-10-18 2008-04-23 中国科学院电子学研究所 星载合成孔径雷达可变位数baq压缩系统及方法
CN101178434A (zh) * 2007-02-02 2008-05-14 上海易罗信息科技有限公司 雷达数据采集系统及其采集方法
JP2008170287A (ja) * 2007-01-12 2008-07-24 Japan Radio Co Ltd レーダ装置

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4053840A (en) * 1975-08-01 1977-10-11 Thomson-Csf Band compression device
US4673941A (en) * 1983-12-07 1987-06-16 Hollandse Signaalapparaten B.V. Digital pulse compression filter
US5850202A (en) * 1997-07-29 1998-12-15 Erim International, Inc. Method for improving SAR system sensitivity in the presence of RF interference
CN1664860A (zh) * 2004-03-05 2005-09-07 中国科学院计算技术研究所 合成孔径雷达复数图像数据的实时自动压缩方法
CN101165510A (zh) * 2006-10-18 2008-04-23 中国科学院电子学研究所 星载合成孔径雷达可变位数baq压缩系统及方法
CN1996042A (zh) * 2006-12-31 2007-07-11 武汉大学 大时间带宽积线性调频信号的数字脉冲压缩方法
JP2008170287A (ja) * 2007-01-12 2008-07-24 Japan Radio Co Ltd レーダ装置
CN101178434A (zh) * 2007-02-02 2008-05-14 上海易罗信息科技有限公司 雷达数据采集系统及其采集方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
丁留贯
丁留贯;顾松山;牛生杰;陈小兰.气象雷达的脉冲压缩技术研究.《大气与环境光学学报》.2007,第2卷(第4期), *
牛生杰
陈小兰.气象雷达的脉冲压缩技术研究.《大气与环境光学学报》.2007,第2卷(第4期),
顾松山

Also Published As

Publication number Publication date
CN101776747A (zh) 2010-07-14

Similar Documents

Publication Publication Date Title
CN101776747B (zh) 基于自适应重采样的雷达信号压缩方法
CN104468427B (zh) 基于fri时频域综合分析的信号高效采样及信号重构方法
CN112466290B (zh) 异常声音检测模型的训练方法、装置和计算机存储介质
CN115883670A (zh) 一种医疗数据分析采集方法及装置
CN117828285B (zh) 一种生产数据管理系统及方法
CN117749253B (zh) 一种遥测数据压缩方法、卫星及地面设备
CN109102818B (zh) 一种基于信号频率概率密度函数分布的去噪音频采样算法
CN112819082A (zh) 一种基于深度学习的卫星频谱感知数据重构方法
CN113923723B (zh) 流量重构方法、装置、设备及存储介质
CN114822578A (zh) 语音降噪方法、装置、设备及存储介质
CN102129073A (zh) 一种合成孔径雷达方位预处理的两级滤波的方法
CN105898211A (zh) 一种处理多媒体信息的方法及装置
CN111049684B (zh) 数据分析方法、装置、设备和存储介质
CN106656201B (zh) 一种基于采样数据幅频特性的压缩方法
CN102565773A (zh) 一种合成孔径雷达高速数据抽取的两级滤波实现方法
CN117278048A (zh) 多导联生物电信号的深度压缩感知方法、设备和介质
CN116904569A (zh) 信号处理方法、装置、电子设备、介质和产品
CN106686029B (zh) 一种生理信号压缩传输、重建方法及装置
CN113506065B (zh) 配送轨迹修正方法及装置
EP2953266B1 (en) Data compression device, data compression method, and program
CN110109803A (zh) 一种用户行为上报方法及系统
CN112116917A (zh) 基于相位跃变度的电抗器本体与风机声信号分离方法
CN113625971B (zh) 数据存储方法、装置、设备以及计算机存储介质
CN215420783U (zh) 胎心检测装置及系统
CN117078781B (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
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20120718

Termination date: 20140115