CN102928834A - 一种基于无相位突变拼接信号的调频连续波测距方法 - Google Patents

一种基于无相位突变拼接信号的调频连续波测距方法 Download PDF

Info

Publication number
CN102928834A
CN102928834A CN2012104829471A CN201210482947A CN102928834A CN 102928834 A CN102928834 A CN 102928834A CN 2012104829471 A CN2012104829471 A CN 2012104829471A CN 201210482947 A CN201210482947 A CN 201210482947A CN 102928834 A CN102928834 A CN 102928834A
Authority
CN
China
Prior art keywords
signal
difference frequency
frequency
ranging method
complex signal
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
CN2012104829471A
Other languages
English (en)
Other versions
CN102928834B (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 Technology BIT
Original Assignee
Beijing Institute of Technology BIT
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 Technology BIT filed Critical Beijing Institute of Technology BIT
Priority to CN201210482947.1A priority Critical patent/CN102928834B/zh
Publication of CN102928834A publication Critical patent/CN102928834A/zh
Application granted granted Critical
Publication of CN102928834B publication Critical patent/CN102928834B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

本发明公开的一种基于无相位突变拼接信号的调频连续波测距方法,涉及一种调频连续波测距方法,属于调频连续波测距领域。本发明具体实现步骤如下:步骤一、将差频实信号变为复信号;步骤二、取规则区内差频复信号;步骤三、用每段信号第一个数据sk(0)除以前一段信号最后一个数据
Figure DDA00002454811800011
1≤k≤N-1;步骤四、通过对步骤三中s'做频谱分析,估计出差频信号频率值fb;步骤五、根据差频信号和目标距离间的关系,计算出目标距离R。本发明的距离分辨率比传统方法的距离分辨率提高了N倍。本发明的测距方法距离分辨率不受限于调制频偏,使得在一些成本和体积受限而无法满足高调制频偏的场合,依然能够使用调频连续波测距得到很高的距离分辨率。

Description

一种基于无相位突变拼接信号的调频连续波测距方法
技术领域
本发明涉及一种调频连续波测距方法,尤其涉及一种基于无相位突变拼接信号的调频连续波测距方法,属于调频连续波测距领域。
背景技术
调频连续波测距系统在高度表、液位测量、防撞雷达、引信等领域有着广泛应用。在此类系统中,发射线性调频信号,发射信号在碰到目标后产生回波,由于回波信号滞后,回波信号频率和发射信号频率之间存在差频,该差频和目标距离之间存在固定关系,准确检测该差频即可准确测距。
目前差频检测方法一般采用快速傅里叶变换(FFT)算法,即选一种窗函数取一段差频信号做FFT,得到其频谱,选取幅值最大的谱线对应频点作为差频估计值,再根据差频和目标距离间的关系,计算出目标距离。但这种方法的距离分辨率受限于调制频偏,二者有关系:ΔRmin=c/2ΔF,其中:ΔRmin是最小距离分辨率,ΔF是调制频偏,c是光速。
在近距离测距时,往往需要很小的距离分辨率,按照目前的方法,就需要很高的调制频偏。而高的调制频偏受限于成本、体积和目前的技术,在很多场合,是难以实现的。
发明内容
本发明的目的是为了解决现有技术中受调制频偏的限制无法实现调频连续波距离测量的问题。本发明公开的一种基于无相位突变拼接信号的调频连续波测距方法,距离分辨率不受限于调制频偏,使得在一些成本和体积受限而无法满足高调制频偏的场合,依然能够使用调频连续波测距得到很高的距离分辨率。
本发明的目的是通过以下技术方案实现的:
在调频连续波测距系统中,发射线性调频信号,发射信号在碰到目标后产生回波,由于回波信号滞后,回波信号频率和发射信号频率之间存在差频,该差频和目标距离之间存在固定关系,准确检测该差频即可准确测距。
本发明的一种基于无相位突变拼接信号的调频连续波测距方法,具体步骤如下:
步骤一、将差频实信号变为复信号。目前已有多种方法,例如在射频电路部分采用正交混频获得复信号,或是在模拟数字转换器ADC采样差频信号后,利用Hilbert变换获得复信号。
步骤二、取规则区内差频复信号。规则区是指回波信号中每个调制周期开始时刻到发射信号每个调制周期结束时刻之间的差频信号。规则区内差频复信号取样可通过下述两种方法实现。
方法1:取步骤一中得到的差频复信号中规则区内的部分,记为
Figure BDA00002454811600021
Figure BDA00002454811600022
T1是每个调制周期规则区的持续时间,Fs是采样频率。将其复制N份,N的值根据期望的距离分辨率和系统的运算能力而定,得到:
s={s0,s1,…,sN-1}
s1,…,sk,…,sN-1是对s0的复制。
方法2:将步骤一中得到的差频复信号中连续的N个调制周期中T1内的部分取出,也记为s={s0,s1,…,sN-1},N的值根据期望的距离分辨率和系统的运算能力而定。
步骤三、用每段信号第一个数据sk(0)除以前一段信号最后一个数据
Figure BDA00002454811600023
1≤k≤N-1。即对N段信号都按公式处理:
s , = { s 0 , s 1 [ s 1 ( 0 ) s 0 ( N T 1 - 1 ) ] - 1 , . . . , s k Π i = 1 i = k [ s i ( 0 ) s i - 1 ( N T 1 - 1 ) ] - 1 , . . . , s N - 1 Π i = 1 i = N - 1 [ s i ( 0 ) s i - 1 ( N T 1 - 1 ) ] - 1 }
因为sk(0)、
Figure BDA00002454811600025
都是复数,所以相位差会包含在
Figure BDA00002454811600026
的e指数上,则s'近似为一段长
Figure BDA00002454811600027
个数据点的正弦信号,其中无相位突变。
步骤四、通过对步骤三中s'做频谱分析,估计出差频信号频率值fb
步骤五、根据差频信号和目标距离间的关系,计算出目标距离R。锯齿波调制时,R=fbTmc/(2ΔF);三角波调制时,R=fbTmc/(4ΔF)。
由于步骤四中求出的频谱只有一个主瓣,且半主瓣宽度为1/(NT1),则当锯齿波调制时,其频率分辨率近似为1/(NT1)≈fm/N,当三角波调制时,其频率分辨率近似为1/(NT1)≈2fm/N。所以该方法的距离分辨率为:
ΔRmin=c/(N2ΔF)                                                 (3)
传统方法的距离分辨率为ΔRmin=c/2ΔF,对比可知,本发明的一种基于无相位突变拼接信号的调频连续波测距方法距离分辨率约为传统方法的1/N。
有益效果
1、传统方法的距离分辨率为ΔRmin=c/2ΔF,本发明的一种基于无相位突变拼接信号的调频连续波测距方法的分辨率为ΔRmin=c/(N2ΔF),本发明的距离分辨率比传统方法的距离分辨率提高了N倍。因此,在成本、体积或现有技术限制而使得调制频偏不能做到很大的情况下,本发明的距离分辨率比传统方法的距离分辨率提高了N倍。
2、本发明的一种基于无相位突变拼接信号的调频连续波测距方法,距离分辨率不受限于调制频偏,使得在一些成本和体积受限而无法满足高调制频偏的场合,依然能够使用调频连续波测距得到很高的距离分辨率。
附图说明
图1是线性调频连续波(以锯齿波为例)发射信号、回波信号和差频信号频率示意图;
图2是线性调频连续波(以三角波为例)发射信号、回波信号和差频信号频率示意图;
其中:f1(t)为发射信号频率,f2(t)为回波信号频率,fb(t)为差频信号频率,ΔF为调频频偏,fc为载波频率,Tm为调制周期,τ为回波延时,T1为差频信号规则区时间,T2为差频信号不规则区时间
具体实施方式
实施例1
调频连续波测距系统,锯齿波调制,调制频偏ΔF=100MHz,调制周期
Figure BDA00002454811600031
按传统方法ΔRmin=c/2ΔF=1.5m。假设该系统测距范围2.5-5m,则回波延时τ=0-3.33×10-8s,差频范围0-666.7KHz,在Fs=40MHz条件下。按照本实施例所述的一种基于无相位突变拼接信号的调频连续波测距方法测距,具体步骤为:
步骤一、将差频实信号变为复信号。目前已有多种方法,这里选用在模拟数字转换器ADC采样差频信号后,利用Hilbert变换获得复信号。即对差频信号做FFT得到频谱,去除频谱的负半轴,再做IFFT,得到复信号。
步骤二、取规则区内的差频复信号。规则区是指回波信号中每个调制周期开始时刻到发射信号每个调制周期结束时刻之间的差频信号。规则区内差频复信号取样可通过下述两种方法实现。
方法1:取步骤一中得到的差频复信号中规则区(图1中T1)内的部分,记为
Figure BDA00002454811600041
每个调制周期对应的数据点数TmFs=200,前后各去除10个数据点,保证去除不规则区(因为τmaxFs=1.33,不规则区大约包含2个数据点),则
Figure BDA00002454811600042
取180。将其复制10份,即s={s0,s1,…,sN-1},一共有1800个数据点。
方法2::将步骤一中得到的差频复信号中连续的10个调制周期中T1内的部分取出,也记为s={s0,s1,…,sN-1}。其余处理同方法1。
步骤三、用每段信号第一个数据sk(0),除以前一段信号最后一个数据sk-1(179),k=1,…,9,对10段信号都按下式处理:
s , = { s 0 , s 1 [ s 1 ( 0 ) s 0 ( 180 - 1 ) ] - 1 , . . . , s 9 Π i = 1 i = 9 [ s i ( 0 ) s i - 1 ( 180 - 1 ) ] - 1 }
因为sk(0)、sk-1(179)都是复数,所以相位差会包含在sk(0)/sk-1(179)的e指数上。则s'近似为一段长1800个数据点的正弦信号,且其中无相位突变。
步骤四、对s'做FFT,估计出差频信号频率值fb
步骤五、根据差频和目标距离间的关系,计算出目标距离。如图1所示,锯齿波调制时,R=fbTmc/(2ΔF);如图2所示,三角波调制时R=fbTmc/(4ΔF)。
本实施例的一种基于无相位突变拼接信号的调频连续波测距方法测距结果与传统方法测距结果比较如下:
实际距离/m 5.00 4.80 4.60 4.40 4.20 4.00 3.80 3.60 3.40 3.20 3.00 2.80
本发明方法/m 4.97 4.85 4.62 4.46 4.18 4.04 3.81 3.60 3.38 3.23 3.06 2.85
由上表可见,在ΔF=100MHz而且Fs=40MHz,只做2048点FFT的条件下,本实施例所述的一种基于无相位突变拼接信号的调频连续波测距方法的距离分辨率至少可达到0.2m。而传统测距方法的距离分辨率根据ΔRmin=c/2ΔF计算,只有1.5m左右。
本实施例的一种基于无相位突变拼接信号的调频连续波测距方法,距离分辨率不受限于调制频偏,使得在一些成本和体积受限而无法满足高调制频偏的场合,依然能够使用调频连续波测距得到很高的距离分辨率。

Claims (5)

1.一种基于无相位突变拼接信号的调频连续波测距方法,其特征在于:具体步骤如下,
步骤一、将差频实信号变为复信号;
步骤二、取规则区内差频复信号;规则区是指回波信号中每个调制周期开始时刻到发射信号每个调制周期结束时刻之间的差频信号;
步骤三、用每段信号第一个数据sk(0)除以前一段信号最后一个数据1≤k≤N-1;即对N段信号都按下式处理:
s , = { s 0 , s 1 [ s 1 ( 0 ) s 0 ( N T 1 - 1 ) ] - 1 , . . . , s k Π i = 1 i = k [ s i ( 0 ) s i - 1 ( N T 1 - 1 ) ] - 1 , . . . , s N - 1 Π i = 1 i = N - 1 [ s i ( 0 ) s i - 1 ( N T 1 - 1 ) ] - 1 }
因为sk(0)、
Figure FDA00002454811500013
都是复数,所以相位差会包含在
Figure FDA00002454811500014
的e指数上,则s'近似为一段长
Figure FDA00002454811500015
个数据点的正弦信号,其中无相位突变;
步骤四、通过对步骤三中s'做频谱分析,估计出差频信号频率值fb
步骤五、根据差频信号和目标距离间的关系,计算出目标距离R。
2.根据权利要求1所述的一种基于无相位突变拼接信号的调频连续波测距方法,其特征在于:步骤二中所述的规则区内差频复信号取样通过下述两种方法中任意一种实现;
方法1:取步骤一中得到的差频复信号中规则区内的部分,记为
Figure FDA00002454811500017
T1是每个调制周期规则区的持续时间,Fs是采样频率;将其复制N份,N的值根据期望的距离分辨率和系统的运算能力而定,得到:
s={s0,s1,…,sN-1}
s1,…,sk,…,sN-1是对s0的复制;
方法2:将步骤一中得到的差频复信号中连续的N个调制周期中T1内的部分取出,也记为s={s0,s1,…,sN-1},N的值根据期望的距离分辨率和系统的运算能力而定。
3.根据权利要求1或2所述的一种基于无相位突变拼接信号的调频连续波测距方法,其特征在于:锯齿波调制时,R=fbTmc/(2ΔF);三角波调制时,R=fbTmc/(4ΔF)。
4.根据权利要求3所述的一种基于无相位突变拼接信号的调频连续波测距方法,其特征在于:步骤一中所述的将差频实信号变为复信号实现方法目前已有多种;例如在射频电路部分采用正交混频获得复信号,或是在模拟数字转换器ADC采样差频信号后,利用Hilbert变换获得复信号。
5.根据权利要求1或2所述的一种基于无相位突变拼接信号的调频连续波测距方法,其特征在于:步骤一中所述的将差频实信号变为复信号实现方法目前已有多种;例如在射频电路部分采用正交混频获得复信号,或是在模拟数字转换器ADC采样差频信号后,利用Hilbert变换获得复信号。
CN201210482947.1A 2012-11-23 2012-11-23 一种基于无相位突变拼接信号的调频连续波测距方法 Expired - Fee Related CN102928834B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210482947.1A CN102928834B (zh) 2012-11-23 2012-11-23 一种基于无相位突变拼接信号的调频连续波测距方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210482947.1A CN102928834B (zh) 2012-11-23 2012-11-23 一种基于无相位突变拼接信号的调频连续波测距方法

Publications (2)

Publication Number Publication Date
CN102928834A true CN102928834A (zh) 2013-02-13
CN102928834B CN102928834B (zh) 2014-05-14

Family

ID=47643677

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210482947.1A Expired - Fee Related CN102928834B (zh) 2012-11-23 2012-11-23 一种基于无相位突变拼接信号的调频连续波测距方法

Country Status (1)

Country Link
CN (1) CN102928834B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103323837A (zh) * 2013-06-27 2013-09-25 北京理工大学 一种基于曲线拟合的调频连续波测距方法
CN104155649A (zh) * 2014-09-02 2014-11-19 中国电子科技集团公司第五十四研究所 一种基于三周期调频连续波相参雷达距离速度去耦合方法
CN105675083A (zh) * 2014-11-20 2016-06-15 中国航空工业集团公司雷华电子技术研究所 一种带有频率细化插值的高精度液位测量方法
CN105738889A (zh) * 2016-02-23 2016-07-06 华域汽车系统股份有限公司 一种调频连续波测速测距方法
CN108445477A (zh) * 2018-01-19 2018-08-24 成都赛英科技有限公司 机场场面异物检测雷达的高精度测距方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1298488A (zh) * 1998-04-28 2001-06-06 丰田自动车株式会社 Fm-cw雷达装置
EP1275979B1 (en) * 2001-07-11 2006-10-04 Fujitsu Ten Limited Radar apparatus
CN101313232A (zh) * 2005-10-24 2008-11-26 三菱电机株式会社 目标测距

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1298488A (zh) * 1998-04-28 2001-06-06 丰田自动车株式会社 Fm-cw雷达装置
EP1275979B1 (en) * 2001-07-11 2006-10-04 Fujitsu Ten Limited Radar apparatus
CN101313232A (zh) * 2005-10-24 2008-11-26 三菱电机株式会社 目标测距

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
孟欣喜等: "LFMCW雷达测距系统及其信号处理算法的设计", 《科学技术与工程》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103323837A (zh) * 2013-06-27 2013-09-25 北京理工大学 一种基于曲线拟合的调频连续波测距方法
CN104155649A (zh) * 2014-09-02 2014-11-19 中国电子科技集团公司第五十四研究所 一种基于三周期调频连续波相参雷达距离速度去耦合方法
CN105675083A (zh) * 2014-11-20 2016-06-15 中国航空工业集团公司雷华电子技术研究所 一种带有频率细化插值的高精度液位测量方法
CN105675083B (zh) * 2014-11-20 2019-03-29 中国航空工业集团公司雷华电子技术研究所 一种带有频率细化插值的高精度液位测量方法
CN105738889A (zh) * 2016-02-23 2016-07-06 华域汽车系统股份有限公司 一种调频连续波测速测距方法
CN105738889B (zh) * 2016-02-23 2018-02-09 华域汽车系统股份有限公司 一种调频连续波测速测距方法
CN108445477A (zh) * 2018-01-19 2018-08-24 成都赛英科技有限公司 机场场面异物检测雷达的高精度测距方法

Also Published As

Publication number Publication date
CN102928834B (zh) 2014-05-14

Similar Documents

Publication Publication Date Title
WO2015190565A1 (ja) チャープ波を使用するレーダ装置及びレーダ装置における信号処理方法
CN103018730B (zh) 分布式子阵波达方向估计方法
CN102928834B (zh) 一种基于无相位突变拼接信号的调频连续波测距方法
CN108415010B (zh) 一种基于梯形lfmcw调制的雷达多目标检测方法
US9746546B2 (en) Method and device for sensing surrounding environment based on frequency modulated continuous wave radar
CN105122079B (zh) 具有间距范围划分的fmcw雷达
CN102508206B (zh) 基于小波包去噪和功率谱熵的线性调频信号参数估计方法
DE102008014918B4 (de) Verfahren zum Erfassen einer Interferenz in einem Radarsystem und Radar, das das gleiche verwendet
DE102012208280B4 (de) Radareinrichtung
CN109375202B (zh) 一种基于车载毫米波雷达的车辆测距测速方法
EP2952926A1 (en) Radar system
CN103675759B (zh) 一种改进的分数阶傅里叶变换机动弱目标检测方法
US9952314B2 (en) Method and device for sensing road environment based on frequency modulated continuous wave radar
WO2014123162A1 (ja) 車両に搭載されるレーダ装置
CN108287335A (zh) 一种利用lfmcw雷达的频率调制信号来对多目标进行测距测速的方法
US7432849B2 (en) Frequency modulation radar apparatus for vehicle use background of the invention
DE102010045980A1 (de) Radarverfahren und Radarsystem
US8884814B2 (en) Processing method for FMCW radar signal with dual pulse repetition frequency
CN104880697B (zh) 基于稀疏约束的线性调频信号参数估计方法
CN109239680B (zh) 一种低截获概率雷达lfm信号的参数估计方法
CN111539323B (zh) 一种循环前缀线性调频信号的频率估计方法与装置
CN103323837A (zh) 一种基于曲线拟合的调频连续波测距方法
CN105675084A (zh) 一种带有迭代插值的高精度液位测量方法
CN105487058B (zh) 能够准确、快速确定雷达系统参数的方法
CN111007473A (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: 20140514

Termination date: 20141123

EXPY Termination of patent right or utility model