CN102798856A - 一种基于小波包的24GHz LFMCW雷达测距方法 - Google Patents

一种基于小波包的24GHz LFMCW雷达测距方法 Download PDF

Info

Publication number
CN102798856A
CN102798856A CN2012103326494A CN201210332649A CN102798856A CN 102798856 A CN102798856 A CN 102798856A CN 2012103326494 A CN2012103326494 A CN 2012103326494A CN 201210332649 A CN201210332649 A CN 201210332649A CN 102798856 A CN102798856 A CN 102798856A
Authority
CN
China
Prior art keywords
signal
frequency
wavelet packet
distance
wavelet
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
CN2012103326494A
Other languages
English (en)
Other versions
CN102798856B (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.)
Chongqing University
Original Assignee
Chongqing 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 Chongqing University filed Critical Chongqing University
Priority to CN201210332649.4A priority Critical patent/CN102798856B/zh
Publication of CN102798856A publication Critical patent/CN102798856A/zh
Application granted granted Critical
Publication of CN102798856B publication Critical patent/CN102798856B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明提出了一种基于小波包的24GHz线性调频连续波LFMCW雷达测距方法，包括以下步骤：将线性调频连续波LFMCW雷达的回波信号与发射信号通过混频和滤波得到差拍信号；利用小波包的分解算法和单节点重构算法对差拍信号进行处理，得到单节点的重构信号；对小波包重构后的信号进行快速傅里叶变换，得到信号的频谱信息；找到频谱的最大幅值所对应的频率，将其代入距离公式计算出目标的距离。本发明提出的测距方法可以在低信噪比下提取出远距离的目标信息，增加雷达的有效作用距离。

Description

一种基于小波包的24GHz LFMCW雷达测距方法
技术领域
本发明涉及一种基于小波包的24GHz LFMCW雷达测距方法，是车辆测距和测速、评估多层的电介质结构、地球表面飞行物的高度和速度的测量以及直升机分类等的一项核心技术，在各类军、民用系统中均可有广泛的应用。
背景技术
毫米波雷达与红外线传感器、超声波传感器以及图像处理传感器设备相比，具有受环境的影响小、高精度和高分辨率等方面的优势。基于线性调频连续波（frequency modulation continue wave，LFMCW）技术的毫米波雷达系统在民用和军事方面都应用得比较广泛。因此，如何进一步提高线性调频连续波LFMCW毫米波雷达的测距精度、增加雷达的有效作用距离仍是一个热门的话题。最常见的方法有：最大似然估计法、曲线拟合法和chirp-z变换法。
影响线性调频连续波LFMCW毫米波雷达测距范围的因素有很多，如调频线性度、收发天线间的隔离度以及信号处理算法等。差拍-傅里叶变换可在频域上等效地完成极大压缩比的脉压功能，因而成为线性调频连续波信号处理的经典方法。差拍信号由发射信号和回波信号混频得到，然而，信号在毫米波段的衰减非常明显，随着目标距离的增加，回波信号的功率会急剧降低，限制了线性调频连续波LFMCW毫米波雷达在远程测距中的应用。
发明内容
本发明的目的在于克服上述线性调频连续波LFMCW毫米波雷达在远程测距中的不足，提供一种基于小波包和傅里叶变换的测距方法，通过小波包分解算法可以实现微弱信号的提取来测量远距离的目标，增加线性调频连续波LFMCW毫米波雷达的有效作用距离。本发明利用小波包的分解算法对差拍信号进行分解，利用单节点的小波包重构算法小波包系数进行重构，利用傅里叶变换对重构信号进行处理，得到目标的频率信息，最终根据距离公式得到目标的距离。
为实现这一目的，本发明的技术方案是以小波包分解算法可以实现微弱信号的提取为切入点。本发明采用以下技术方案：
一种基于小波包的24GHz 线性调频连续波LFMCW雷达测距方法，其特征在于，所述测距方法包括如下步骤：
A、 线性调频连续波LFMCW雷达的发射信号经天线发射出去，遇到目标时形成回波信号，将回波信号与发射信号通过混频和滤波得到差拍信号；
B、 选择合适的小波函数和小波分解层数，对差拍信号进行小波包分解，得到一系列子带上的小波包系数
Figure 136478DEST_PATH_IMAGE001
其中，
Figure 933533DEST_PATH_IMAGE002
，
Figure 227855DEST_PATH_IMAGE003
为差拍信号，
Figure 339030DEST_PATH_IMAGE004
小波包系数，
Figure 349712DEST_PATH_IMAGE005
为小波包函数；
C、 将雷达的有效测距范围划分为若干段，利用距离公式计算出每段距离所对应的频带范围，选择频域与之对应的小波包系数
Figure 504618DEST_PATH_IMAGE006
其中，
Figure 600750DEST_PATH_IMAGE007
为目标距离，
Figure 250037DEST_PATH_IMAGE008
为差拍信号的频率，
Figure 803641DEST_PATH_IMAGE009
为光速，
Figure 942498DEST_PATH_IMAGE010
为调频带宽, 
Figure 463609DEST_PATH_IMAGE011
为调制周期，
Figure 900276DEST_PATH_IMAGE012
为调频斜率；
D、 根据目标所处的距离段来选择小波包系数，对其进行小波包重构，得到单节点的重构信号
Figure 885549DEST_PATH_IMAGE013
；
E、  对小波包重构后的信号进行快速傅里叶变换，得到信号的频谱信息
Figure 132991DEST_PATH_IMAGE014
其中，
Figure 203715DEST_PATH_IMAGE015
；
F、  利用冒泡法找到频谱的最大幅值所对应的频率，将其代入距离公式计算出目标的距离。
附图说明
图1是差拍信号的小波包分解与单节点重构快递算法；
图2是目标距离雷达6.39m时的差拍信号频谱图；
图3是目标距离雷达58.01m时的差拍信号频谱图；
图4是雷达测距方法的流程图。
具体实施方式
本实施方式是以本发明技术方案为前提进行实施，下面结合附图对本发明的具体实施方式做详细说明：
本实施方式是24GHz 线性调频连续波LFMCW雷达对目标进行测距，包括以下步骤：
1）线性调频连续波LFMCW雷达的发射信号经天线发射出去，遇到目标时形成回波信号，将回波信号与发射信号通过混频和滤波得到差拍信号。
2）选择db3小波和小波分解层数
Figure 966045DEST_PATH_IMAGE016
，对差拍信号
Figure 743508DEST_PATH_IMAGE017
进行小波包分解，得到一系列子带上的小波包系数
Figure 224168DEST_PATH_IMAGE018
：
Figure 703560DEST_PATH_IMAGE019
    
式中，
Figure 498341DEST_PATH_IMAGE020
；
Figure 192627DEST_PATH_IMAGE005
为小波包函数，它可以由小波函数
Figure 532604DEST_PATH_IMAGE021
经滤波器和
Figure 910813DEST_PATH_IMAGE023
逐次计算得到：
式中，
Figure 734598DEST_PATH_IMAGE022
、
Figure 939314DEST_PATH_IMAGE023
为小波重构滤波器，
Figure 138215DEST_PATH_IMAGE022
与尺度函数有关，
Figure 423626DEST_PATH_IMAGE023
与小波函数有关。
3）将雷达的有效测距范围划分为
Figure 101918DEST_PATH_IMAGE027
段，利用距离公式计算出每段距离所对应的频带范围，选择频域与之相对应的小波包系数，
Figure 562986DEST_PATH_IMAGE028
其中，为目标距离，为差拍信号的频率，
Figure 135678DEST_PATH_IMAGE009
为光速，
Figure 779149DEST_PATH_IMAGE029
为调频带宽, 
Figure 301266DEST_PATH_IMAGE030
为调制周期，
Figure 30188DEST_PATH_IMAGE031
为调频斜率。
 4）测距时，根据目标所处的距离段来选择小波包系数，对其进行小波包重构，得到单节点的重构信号
Figure 46685DEST_PATH_IMAGE032
，
差拍信号的小波包分解与单节点重构快递算法如图1所示。
5）对小波包重构后的信号
Figure 736216DEST_PATH_IMAGE032
进行快速傅里叶变换，得到信号的频谱信息，
Figure 890117DEST_PATH_IMAGE034
  
式中，
Figure 959573DEST_PATH_IMAGE015
，目标距离雷达6.39m时的差拍信号频谱图如图2所示，目标距离雷达58.01m时的差拍信号频谱图如图3所示。
6）利用冒泡法找出频谱的最大幅值
Figure 312057DEST_PATH_IMAGE035
所对应的频率，即重构信号的频率
Figure 192289DEST_PATH_IMAGE036
，并将其代入距离公式计算出目标的距离，整个算法流程如图4所示。
通过在自制的24GHz 线性调频连续波LFMCW雷达实验平台上验证，与传统的差拍-傅里叶测距方法相比，本方法可以将实验雷达的有效作用范围从1m-30m扩大至1m-60m。

Claims (1)

1.一种基于小波包的24GHz 线性调频连续波LFMCW雷达测距方法，其特征在于，所述检测方法包括如下步骤：
线性调频连续波LFMCW雷达的发射信号经天线发射出去，遇到目标时形成回波信号，将回波信号与发射信号通过混频和滤波得到差拍信号；
选择合适的小波函数和小波分解层数，对差拍信号进行小波包分解，得到一系列子带上的小波包系数
Figure 572367DEST_PATH_IMAGE001
其中，
Figure 2012103326494100001DEST_PATH_IMAGE002
，
Figure 690627DEST_PATH_IMAGE003
为差拍信号，
Figure 2012103326494100001DEST_PATH_IMAGE004
小波包系数，
Figure 869935DEST_PATH_IMAGE005
为小波包函数；
将雷达的有效测距范围划分为若干段，利用距离公式计算出每段距离所对应的频带范围，选择频域与之对应的小波包系数
Figure 2012103326494100001DEST_PATH_IMAGE006
其中，
Figure 269692DEST_PATH_IMAGE007
为目标距离，
Figure 2012103326494100001DEST_PATH_IMAGE008
为差拍信号的频率，
Figure 598649DEST_PATH_IMAGE009
为光速，
Figure 2012103326494100001DEST_PATH_IMAGE010
为调频带宽, 
Figure 871499DEST_PATH_IMAGE011
为调制周期，
Figure 2012103326494100001DEST_PATH_IMAGE012
为调频斜率；
根据目标所处的距离段来选择小波包系数，对其进行小波包重构，得到单节点的重构信号
Figure 2012103326494100001DEST_PATH_IMAGE013
；
对小波包重构后的信号进行快速傅里叶变换，得到信号的频谱信息
Figure 2012103326494100001DEST_PATH_IMAGE014
其中，
Figure 2012103326494100001DEST_PATH_IMAGE015
；
利用冒泡法找到频谱的最大幅值所对应的频率，将其代入距离公式计算出目标的距离。
CN201210332649.4A 2012-09-11 2012-09-11 一种基于小波包的24GHz LFMCW雷达测距方法 Expired - Fee Related CN102798856B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210332649.4A CN102798856B (zh) 2012-09-11 2012-09-11 一种基于小波包的24GHz LFMCW雷达测距方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210332649.4A CN102798856B (zh) 2012-09-11 2012-09-11 一种基于小波包的24GHz LFMCW雷达测距方法

Publications (2)

Publication Number Publication Date
CN102798856A true CN102798856A (zh) 2012-11-28
CN102798856B CN102798856B (zh) 2014-08-20

Family

ID=47198012

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210332649.4A Expired - Fee Related CN102798856B (zh) 2012-09-11 2012-09-11 一种基于小波包的24GHz LFMCW雷达测距方法

Country Status (1)

Country Link
CN (1) CN102798856B (zh)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103914983A (zh) * 2014-04-17 2014-07-09 西安飞达电子科技有限公司 一种交通卡口管理雷达装置及实现方法
CN104155649A (zh) * 2014-09-02 2014-11-19 中国电子科技集团公司第五十四研究所 一种基于三周期调频连续波相参雷达距离速度去耦合方法
CN104793198A (zh) * 2015-02-05 2015-07-22 北京理工大学 测距方法及装置、引信控制系统
CN106908794A (zh) * 2015-12-23 2017-06-30 株式会社万都 用于检测目标物体的方法和设备
CN107424657A (zh) * 2017-07-21 2017-12-01 中国工程物理研究院核物理与化学研究所 一种新型反应堆核噪声频谱分析仪
CN111208084A (zh) * 2020-01-15 2020-05-29 大连理工大学 一种基于相干探测方法的光纤气体浓度遥感探测装置和方法
CN111693955A (zh) * 2020-05-25 2020-09-22 北京理工大学 一种太赫兹调频连续波雷达高分辨率距离向重建方法
CN111929685A (zh) * 2020-08-03 2020-11-13 启航汽车有限公司 基于线性调频连续波的雷达检测方法及装置、雷达装置和机动车
CN112068128A (zh) * 2020-09-19 2020-12-11 重庆大学 一种直道场景线段型雷达数据处理及位姿获取方法
CN112946640A (zh) * 2021-01-05 2021-06-11 重庆大学 一种快速距离-多普勒域的自旋目标isar成像方法
CN113204017A (zh) * 2021-03-16 2021-08-03 西安电子科技大学 一种基于lfmcw的mimo雷达信号处理方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4535417A (en) * 1982-12-23 1985-08-13 Standard Oil Company Method and apparatus for digital time-variant filtering
US20070273567A1 (en) * 2006-05-26 2007-11-29 Lecroy Corporation Adaptive interpolation
CN102508206A (zh) * 2011-10-31 2012-06-20 哈尔滨工程大学 基于小波包去噪和功率谱熵的线性调频信号参数估计方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4535417A (en) * 1982-12-23 1985-08-13 Standard Oil Company Method and apparatus for digital time-variant filtering
US20070273567A1 (en) * 2006-05-26 2007-11-29 Lecroy Corporation Adaptive interpolation
CN102508206A (zh) * 2011-10-31 2012-06-20 哈尔滨工程大学 基于小波包去噪和功率谱熵的线性调频信号参数估计方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
《中国优秀博硕士学位论文全文数据库 (博士) 工程科技Ⅱ辑》 20030630 贺乐厅 智能运输系统--基于毫米波雷达的车辆防撞技术与实验研究 第2003年卷, 第2期 *
徐庆: "信号设计及高分辨雷达目标信息提取技术研究", 《中国优秀博硕士学位论文全文数据库 (博士) 信息科技辑》 *
江登银: "汽车自动防撞雷达系统的研究", 《中国优秀博硕士学位论文全文数据库 (硕士)》 *
贺乐厅: "智能运输系统——基于毫米波雷达的车辆防撞技术与实验研究", 《中国优秀博硕士学位论文全文数据库 (博士) 工程科技Ⅱ辑》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103914983A (zh) * 2014-04-17 2014-07-09 西安飞达电子科技有限公司 一种交通卡口管理雷达装置及实现方法
CN104155649A (zh) * 2014-09-02 2014-11-19 中国电子科技集团公司第五十四研究所 一种基于三周期调频连续波相参雷达距离速度去耦合方法
CN104793198A (zh) * 2015-02-05 2015-07-22 北京理工大学 测距方法及装置、引信控制系统
CN106908794A (zh) * 2015-12-23 2017-06-30 株式会社万都 用于检测目标物体的方法和设备
CN107424657A (zh) * 2017-07-21 2017-12-01 中国工程物理研究院核物理与化学研究所 一种新型反应堆核噪声频谱分析仪
CN111208084A (zh) * 2020-01-15 2020-05-29 大连理工大学 一种基于相干探测方法的光纤气体浓度遥感探测装置和方法
CN111693955A (zh) * 2020-05-25 2020-09-22 北京理工大学 一种太赫兹调频连续波雷达高分辨率距离向重建方法
CN111929685A (zh) * 2020-08-03 2020-11-13 启航汽车有限公司 基于线性调频连续波的雷达检测方法及装置、雷达装置和机动车
CN112068128A (zh) * 2020-09-19 2020-12-11 重庆大学 一种直道场景线段型雷达数据处理及位姿获取方法
CN112946640A (zh) * 2021-01-05 2021-06-11 重庆大学 一种快速距离-多普勒域的自旋目标isar成像方法
CN113204017A (zh) * 2021-03-16 2021-08-03 西安电子科技大学 一种基于lfmcw的mimo雷达信号处理方法

Also Published As

Publication number Publication date
CN102798856B (zh) 2014-08-20

Similar Documents

Publication Publication Date Title
CN102798856B (zh) 一种基于小波包的24GHz LFMCW雷达测距方法
US10228449B2 (en) Method and system for jointly separating noise from signals
CN109031296B (zh) 基于随机间歇采样线性调频信号的宽带雷达成像方法
US9363024B2 (en) Method and system for estimation and extraction of interference noise from signals
CN108132461B (zh) 抑制调频连续波着陆雷达直流泄露的方法
CN101900813B (zh) 基于机动目标距离-瞬时调频的isar成像方法
CN105866749A (zh) 一种距离和速度同步拖引的雷达干扰识别方法
CN103389492B (zh) 多通道随机谐波调制采样雷达接收机及其方法
CN103412295A (zh) 基于回波精确模型的高速机动弱目标检测方法
US20160018512A1 (en) Method for Generating and Compressing Multi-Sweep-Frequency Radar Signals
CN107894586A (zh) 一种基于同步压缩变换的激光雷达回波信号去噪方法
CN106842163B (zh) 一种弹道目标回波信号时频特性估计方法
CN106199602B (zh) 一种Staggered-SAR回波信号重建方法
CN104678395A (zh) 基于循环前缀的mimo-ofdm雷达成像方法
CN106054154B (zh) 一种基于icpf的机动目标的步进频率信号高分辨率成像方法
CN103091665A (zh) 一种压缩感知合成孔径雷达射频干扰抑制处理方法
CN102269808B (zh) 一种风廓线雷达信号探测中抗有源干扰方法
CN106680791A (zh) 基于宽带扫频数据的雷达回波仿真方法
CN109254284A (zh) 用于编码超声波回波检测的方法
CN109085568B (zh) 一种基于二次混频的调频连续波多目标检测方法
CN113848532A (zh) 基于降噪模型的fmcw雷达信号降噪系统及方法
CN103412289A (zh) 基于模拟sinc核采样的雷达高分辨距离成像方法
CN103885044B (zh) 一种基于clean算法的窄带雷达回波杂噪抑制方法
Li et al. Deep learning for interference mitigation in time-frequency maps of fmcw radars
CN107769815A (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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20140820

Termination date: 20160911

CF01 Termination of patent right due to non-payment of annual fee