CN102928834B - 一种基于无相位突变拼接信号的调频连续波测距方法 - Google Patents
一种基于无相位突变拼接信号的调频连续波测距方法 Download PDFInfo
- Publication number
- CN102928834B CN102928834B CN201210482947.1A CN201210482947A CN102928834B CN 102928834 B CN102928834 B CN 102928834B CN 201210482947 A CN201210482947 A CN 201210482947A CN 102928834 B CN102928834 B CN 102928834B
- Authority
- CN
- China
- Prior art keywords
- signal
- difference frequency
- frequency
- complex signal
- ranging method
- 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
Links
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
- Optical Radar Systems And Details Thereof (AREA)
Abstract
本发明公开的一种基于无相位突变拼接信号的调频连续波测距方法,涉及一种调频连续波测距方法,属于调频连续波测距领域。本发明具体实现步骤如下:步骤一、将差频实信号变为复信号;步骤二、取规则区内差频复信号;步骤三、用每段信号第一个数据sk(0)除以前一段信号最后一个数据1≤k≤N-1;步骤四、通过对步骤三中s'做频谱分析,估计出差频信号频率值fb;步骤五、根据差频信号和目标距离间的关系,计算出目标距离R。本发明的距离分辨率比传统方法的距离分辨率提高了N倍。本发明的测距方法距离分辨率不受限于调制频偏,使得在一些成本和体积受限而无法满足高调制频偏的场合,依然能够使用调频连续波测距得到很高的距离分辨率。
Description
技术领域
本发明涉及一种调频连续波测距方法,尤其涉及一种基于无相位突变拼接信号的调频连续波测距方法,属于调频连续波测距领域。
背景技术
调频连续波测距系统在高度表、液位测量、防撞雷达、引信等领域有着广泛应用。在此类系统中,发射线性调频信号,发射信号在碰到目标后产生回波,由于回波信号滞后,回波信号频率和发射信号频率之间存在差频,该差频和目标距离之间存在固定关系,准确检测该差频即可准确测距。
目前差频检测方法一般采用快速傅里叶变换(FFT)算法,即选一种窗函数取一段差频信号做FFT,得到其频谱,选取幅值最大的谱线对应频点作为差频估计值,再根据差频和目标距离间的关系,计算出目标距离。但这种方法的距离分辨率受限于调制频偏,二者有关系:ΔRmin=c/2ΔF,其中:ΔRmin是最小距离分辨率,ΔF是调制频偏,c是光速。
在近距离测距时,往往需要很小的距离分辨率,按照目前的方法,就需要很高的调制频偏。而高的调制频偏受限于成本、体积和目前的技术,在很多场合,是难以实现的。
发明内容
本发明的目的是为了解决现有技术中受调制频偏的限制无法实现调频连续波距离测量的问题。本发明公开的一种基于无相位突变拼接信号的调频连续波测距方法,距离分辨率不受限于调制频偏,使得在一些成本和体积受限而无法满足高调制频偏的场合,依然能够使用调频连续波测距得到很高的距离分辨率。
本发明的目的是通过以下技术方案实现的:
在调频连续波测距系统中,发射线性调频信号,发射信号在碰到目标后产生回波,由于回波信号滞后,回波信号频率和发射信号频率之间存在差频,该差频和目标距离之间存在固定关系,准确检测该差频即可准确测距。
本发明的一种基于无相位突变拼接信号的调频连续波测距方法,具体步骤如下:
步骤一、将差频实信号变为复信号。目前已有多种方法,例如在射频电路部分采用正交混频获得复信号,或是在模拟数字转换器ADC采样差频信号后,利用Hilbert变换获得复信号。
步骤二、取规则区内差频复信号。规则区是指回波信号中每个调制周期开始时刻到发射信号每个调制周期结束时刻之间的差频信号。规则区内差频复信号取样可通过下述两种方法实现。
s={s0,s1,…,sN-1}
s1,…,sk,…,sN-1是对s0的复制。
方法2:将步骤一中得到的差频复信号中连续的N个调制周期中T1内的部分取出,也记为s={s0,s1,…,sN-1},N的值根据期望的距离分辨率和系统的运算能力而定。
步骤四、通过对步骤三中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,调制周期按传统方法Δ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)内的部分,记为每个调制周期对应的数据点数TmFs=200,前后各去除10个数据点,保证去除不规则区(因为τmaxFs=1.33,不规则区大约包含2个数据点),则取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段信号都按下式处理:
因为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 (4)
1.一种基于无相位突变拼接信号的调频连续波测距方法,其特征在于:具体步骤如下,
步骤一、将差频实信号变为复信号;
步骤二、取规则区内差频复信号;规则区是指回波信号中每个调制周期开始时刻到发射信号每个调制周期结束时刻之间的差频信号;
步骤四、通过对步骤三中s'做频谱分析,估计出差频信号频率值fb;
步骤五、根据差频信号和目标距离间的关系,计算出目标距离R;
步骤二中所述的取规则区内差频复信号通过下述两种方法中任意一种实现;
方法1:取步骤一中得到的差频复信号中规则区内的部分,记为 T1是每个调制周期规则区的持续时间,Fs是采样频率;将其复制N份,N的值根据期望的距离分辨率和系统的运算能力而定,得到:
s={s0,s1,…,sN-1}
s1,…,sk,…,sN-1是对s0的复制;
方法2:将步骤一中得到的差频复信号中连续的N个调制周期中T1内的部分取出,也记为s={s0,s1,…,sN-1},N的值根据期望的距离分辨率和系统的运算能力而定。
2.根据权利要求1所述的一种基于无相位突变拼接信号的调频连续波测距方法,其特征在于:锯齿波调制时,R=fbTmc/(2ΔF);三角波调制时,R=fbTmc/(4ΔF);其中,Tm为调制周期,ΔF为调频频偏。
3.根据权利要求1所述的一种基于无相位突变拼接信号的调频连续波测距方法,其特征在于:步骤一中所述的将差频实信号变为复信号实现方法为:在射频电路部分采用正交混频获得复信号,或是在模拟数字转换器ADC采样差频信号后,利用Hilbert变换获得复信号。
4.根据权利要求1或2所述的一种基于无相位突变拼接信号的调频连续波测距方法,其特征在于:步骤一中所述的将差频实信号变为复信号实现方法为:在射频电路部分采用正交混频获得复信号,或是在模拟数字转换器ADC采样差频信号后,利用Hilbert变换获得复信号。
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 CN102928834A (zh) | 2013-02-13 |
CN102928834B true 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) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103323837A (zh) * | 2013-06-27 | 2013-09-25 | 北京理工大学 | 一种基于曲线拟合的调频连续波测距方法 |
CN104155649A (zh) * | 2014-09-02 | 2014-11-19 | 中国电子科技集团公司第五十四研究所 | 一种基于三周期调频连续波相参雷达距离速度去耦合方法 |
CN105675083B (zh) * | 2014-11-20 | 2019-03-29 | 中国航空工业集团公司雷华电子技术研究所 | 一种带有频率细化插值的高精度液位测量方法 |
CN105738889B (zh) * | 2016-02-23 | 2018-02-09 | 华域汽车系统股份有限公司 | 一种调频连续波测速测距方法 |
CN108445477B (zh) * | 2018-01-19 | 2022-02-01 | 成都赛英科技有限公司 | 机场场面异物检测雷达的高精度测距方法 |
CN116045752A (zh) * | 2022-12-30 | 2023-05-02 | 北京理工大学 | 一种三角波调频引信的差频信号的模拟方法和系统 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3534164B2 (ja) * | 1998-04-28 | 2004-06-07 | トヨタ自動車株式会社 | Fm−cwレーダ装置 |
JP2003028951A (ja) * | 2001-07-11 | 2003-01-29 | Fujitsu Ten Ltd | レーダ装置 |
EP1777549B1 (en) * | 2005-10-24 | 2012-10-03 | Mitsubishi Electric Information Technology Centre Europe B.V. | Object ranging |
-
2012
- 2012-11-23 CN CN201210482947.1A patent/CN102928834B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN102928834A (zh) | 2013-02-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102928834B (zh) | 一种基于无相位突变拼接信号的调频连续波测距方法 | |
CN108415010B (zh) | 一种基于梯形lfmcw调制的雷达多目标检测方法 | |
CN103018730B (zh) | 分布式子阵波达方向估计方法 | |
US9746546B2 (en) | Method and device for sensing surrounding environment based on frequency modulated continuous wave radar | |
JP6364986B2 (ja) | レーダ装置 | |
CN105122079B (zh) | 具有间距范围划分的fmcw雷达 | |
CN109375202B (zh) | 一种基于车载毫米波雷达的车辆测距测速方法 | |
CN103675759B (zh) | 一种改进的分数阶傅里叶变换机动弱目标检测方法 | |
US20170299711A1 (en) | Radar device and position-determination method | |
DE102014009648A1 (de) | Verfahren und Vorrichtung zum Abtasten einer Straßenumgebung basierend auf einem frequenzmodulierten Dauerstrichradar | |
CN102508206A (zh) | 基于小波包去噪和功率谱熵的线性调频信号参数估计方法 | |
US20200379081A1 (en) | Detection of interference-induced perturbations in fmcw radar systems | |
US11994576B2 (en) | Radar systems and methods | |
CN104215952B (zh) | 基于微动特性的车载目标识别系统及其识别方法 | |
CN102226839A (zh) | 一种低采样率线扫频脉冲时延估计方法 | |
Enggar et al. | Performance comparison of various windowing On FMCW radar signal processing | |
US8884814B2 (en) | Processing method for FMCW radar signal with dual pulse repetition frequency | |
CN109239680B (zh) | 一种低截获概率雷达lfm信号的参数估计方法 | |
CN104880697A (zh) | 基于稀疏约束的线性调频信号参数估计方法 | |
CN111539323B (zh) | 一种循环前缀线性调频信号的频率估计方法与装置 | |
CN117008130A (zh) | 一种基于毫米波雷达的车辆目标速度变化测量方法及系统 | |
CN103323837A (zh) | 一种基于曲线拟合的调频连续波测距方法 | |
CN109085568B (zh) | 一种基于二次混频的调频连续波多目标检测方法 | |
CN106814350A (zh) | 基于压缩感知的外辐射源雷达参考信号信杂比估计方法 | |
CN105675084A (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 |