CN103823216A - 一种调频连续波雷达系统测距方法 - Google Patents
一种调频连续波雷达系统测距方法 Download PDFInfo
- Publication number
- CN103823216A CN103823216A CN201410084559.7A CN201410084559A CN103823216A CN 103823216 A CN103823216 A CN 103823216A CN 201410084559 A CN201410084559 A CN 201410084559A CN 103823216 A CN103823216 A CN 103823216A
- Authority
- CN
- China
- Prior art keywords
- discrete
- spectrum
- wire size
- frequency
- freuqncy 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/06—Systems determining position data of a target
- G01S13/08—Systems for measuring distance only
- G01S13/32—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/35—Details of non-pulse systems
- G01S7/352—Receivers
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO 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/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/35—Details of non-pulse systems
- G01S7/352—Receivers
- G01S7/356—Receivers involving particularities of FFT processing
Landscapes
- Engineering & Computer Science (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Computer Networks & Wireless Communication (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提供一种调频连续波雷达系统测距方法,该方法包括步骤:步骤S1对调频连续波雷达系统采集的离散中频信号进行离散傅里叶变换,找到离散频谱最大值谱线号km和离散频谱次大值谱线号kc;步骤S2若号kc=km-1进入步骤S4;若kc≠km-1进入步骤S3;步骤S3对离散中频信号的频域向右进行Δf/2的移位,使移位后的谱线号km和kc满足kc=km-1,其中Δf为N点离散傅里叶变换的频率分辨率;步骤S4对离散中频信号减少一个采样点数的离散傅里叶变换,找到离散频谱最大值谱线号k′m;步骤S5对离散频谱最大值谱线号k′m进行判断,若k′m=km重复步骤S4,若k′m≠km,k′m=km-1,进入步骤S6;步骤S6根据离散傅里叶变换的采样点数,再结合距离关系式计算出所测雷达与目标的距离区间(R1,R2),则所测雷达与目标的测距离值为R′=(R1+R2)/2。
Description
技术领域
本发明属于调频连续波(FMCW)雷达系统测距技术领域,特别是用于由于环境因素等影响进行的非接触式测距。该方法可以用于液位计测距计算,也可以用于物位测距计算,具有计算精度高、计算量小、抗噪声性能好等优点。
背景技术
线性调频连续波雷达具有许多其他雷达不具备的优点:无距离盲区、距离分辨力高、辐射功率小。用雷达测量距离是非接触测量,因此不会对目标物体产生损伤,可用于很多特殊环境下,实现很高的精度测量。其原理是:向目标发送线性变换的调频连续波,电磁波到达目标后返回,接收其回波,将回波与本振信号混频得到中频信号,距离信息就包含在这个中频信号中,对中频信号进行处理就可以提取出距离信息。
由于液位和物位测距要求的精度较高,但是离散傅里叶变换具有栅栏效应,选取的峰值点只在测量距离为测量精度整数倍时才是准确的,其他距离时都会出现固有的系统误差。为了减小栅栏效应带来的误差,研究人员提出频谱细化算法:Z域的快速傅里叶变换(ZFFT)、线性调频Z变换(CZT)或者通过补零来细化频谱,从而找到频谱最大值所对应的频率。频谱细化法实则是通过中频回波的时域信号对信号的频域进行插值,在没有噪声的情况下,可以人为调节插值的精度,进而控制算法的精度。在这几种频谱细化的方法中,ZFFT需要移位后进行低通滤波,滤掉高频部分干扰,然后再通过重采样得到所需要的细化频谱。补零频谱细化是根据系统要求的细化倍数补充相应倍数的零,然后进行离散傅里叶变换,可以减小频谱间隔,在所计算的频谱范围内增加更多的谱线,进而提高计算精度。
发明内容
(一)要解决的技术问题
现有技术最主要的技术缺陷:由于ZFFT、CZT和补零频谱细化需要进行离散傅里叶变换,通过插值计算出非整数倍出的频谱幅度,但是在最大值附近的谱线幅度相差较小,由于频谱泄露和噪声对频谱的影响,容易出现对最大值的误判,影响算法的精度,为此本发明提供一种调频连续波雷达系统测距方法。
(二)技术方案
本发明提供一种调频连续波雷达系统测距方法,包括步骤如下:
步骤S1:确定频谱主瓣内的最大值谱线号和次大值谱线号,对调频连续波雷达系统采集的离散中频信号进行离散傅里叶变换(DFT),找到离散频谱最大值谱线号km和离散频谱次大值谱线号kc;
步骤S2:对谱线号km和kc进行判断,若谱线号kc≠km-1,则进入步骤S3;若谱线号kc=km-1,则进入步骤S4;
步骤S3:对离散中频信号的频域向右进行Δf/2的移位,使移位后的离散频谱最大值谱线号km和离散频谱次大值谱线号kc满足kc=km-1,其中Δf为N点离散傅里叶变换的频率分辨率;
步骤S4:对离散中频信号减少一个采样点数的离散傅里叶变换,找到离散频谱最大值谱线号k′m;
步骤S5:对离散频谱最大值谱线号k′m进行判断,若谱线号k′m=km,则重复步骤S4,若谱线号k′m≠km,则谱线号k′m=km-1,进入步骤S6;
步骤S6:根据步骤S4中进行离散傅里叶变换的采样点数,再结合距离关系式计算出所测雷达与目标的距离区间(R1,R2),则所测雷达与目标的测距离值为
(三)有益效果
在实际计算离散傅里叶变换时,本发明可以利用快速傅里叶变换(FFT)代替,快速傅里叶变换是离散傅里叶变换的一种快速运算方法,可以减小运算量,提高运算效率。
a、本发明说计算出的距离区间(R1,R2),其区间长度 与直接进行离散傅里叶变换计算距离相比,精度提高了大约倍,一般N的值较大,而M和km较小,因此精度提高一般大于25倍。
b、由于本发明经过离散傅里叶变换,使噪声主要分布在频谱的高频段,而中频信号的频率属于低频区间,因此可以在存在一定噪声影响的情况下,不至于对计算精度产生太大的影响,可以保证计算结果的可靠性。
附图说明
图1为本发明非均匀频谱细化FMCW测距算法流程图;
图2为本发明离散中频信号频谱图;
图3为本发明离散中频信号非均匀频谱细化频谱图;
图4为本发明的非均匀频谱细化频谱最大值的谱线号变化图;
图5为本发明的仿真误差图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明进一步详细说明。
本发明提出的非均匀频谱细化FMCW测距算法是基于DFT变换后中频信号频率的分布区间来确定频率的一种新型测距方法,该方法通过实验仿真,验证了理论的正确性,理论分析和仿真实验结果证明了该方法的有效性。
请参照图1示出本发明移位能量归一化调频连续波雷达系统测距方法的流程图,该方法包括以下步骤
步骤S1:确定频谱主瓣内的最大值谱线号和次大值谱线号,对调频连续波雷达系统采集的离散中频信号进行离散傅里叶变换(DFT),找到离散频谱最大值谱线号km和离散频谱次大值谱线号kc;
步骤S2:对谱线号km和kc进行判断,若谱线号kc≠km-1,则进入步骤S3;若谱线号kc=km-1,则进入步骤S4;
步骤S3:对离散中频信号的频域向右进行Δf/2的移位,使移位后的离散频谱最大值谱线号km和离散频谱次大值谱线号kc满足kc=km-1,其中Δf为N点离散傅里叶变换的频率分辨率;
步骤S4:对离散中频信号减少一个采样点数的离散傅里叶变换,找到离散频谱最大值谱线号k′m;
步骤S5:对离散频谱最大值谱线号k′m进行判断,若谱线号k′m=km,则重复步骤S4,若谱线号k′m≠km,则谱线号k′m=km-1,进入步骤S6;
步骤S6;根据步骤S4中进行离散傅里叶变换的采样点数,再结合距离关系式计算出所测雷达与目标的距离区间(R1,R2),则所测雷达与目标的距离值为
在步骤S1,所述采集离散中频信号是对中频信号进行离散采样,从而得到离散中频信号其中:t为连续信号时间,π为圆周率,R为雷达与目标的实际距离,c为光速,B为发射信号带宽,T为发射信号扫描周期,f1为起始频率,N为采样点数;n为离散中频信号的时间点且n=0,2…N-1;采样率为中频信号x(t)的频率f与距离R的关系式为
所述找到离散频谱最大值谱线号km和离散频谱次大值谱线号kc的步骤如下:
e为自然对数的底数,j为虚数单位,|·|表示取复数的模值;
步骤S102:根据得到的离散频谱X(k),找到离散频谱X(k)中最大值所对应的谱线号km和次大值所对应的谱线号kc。
在步骤3;根据离散余弦函数的频谱近似为sinc函数,sinc函数关于最大值所在的轴对称,得知频谱次大值kc=km±1,因此当kc=km+1时,则根据离散傅里叶变换的移位公式性质,对离散中频信号频谱x(n)向右进行Δf/2的移位表示如下:其中Δf为N点频率分辨率因此移位后的最大值谱线号仍然为km,移位后的次大值kc满足kc=km-1。
从图2示出本发明的离散中频信号频谱图,表示N点离散傅里叶变换的部分频谱,可以看出离散傅里叶频谱存在栅栏效应,即只有频率为频率分辨率Fs/N的整数倍才进行采样,而在非整数倍频率处不采样,称作栅栏效应。由于实际情况中,离散中频信号频率为分辨率非整数倍居多,因此直接利用离散傅里叶频谱最大值进行计算会产生较大的误差。通过步骤S1、步骤S2和步骤S3的处理后谱线如图2所示,次大值的谱线号kc比最大值的谱线号km小1。
步骤4中,对离散中频信号x(n)减少采样点数的离散傅里叶变换并取模值,根据离散中频信号离散频谱从而找到离散中频信号离散频谱X(k)最大值的谱线号为k′m;由于进行离散傅里叶变换的采样点数在减小,而采样率不变,频率间隔在逐渐增大,因此将这种减小采样点数的离散傅里叶变换称作非均匀频谱细化,其中i为进行减小采样点数离散傅里叶变换时减少的采样点数,WN-i为旋转因子满足欧拉公式: 离散频谱最大值谱线号所对应的频率为
如图3示出本发明离散中频信号的非均匀频谱细化频谱图中作出了在幅度最大值频谱号跳变前一个周期内的非均匀频谱图,从图中可以看出,频谱细化不是均匀的,而是有些地方窄有些地方宽,从而仿真证明了非均匀频谱细化的理论性,但是其中的幅度值最大的点不代表离散中频信号的频率,需要利用其他关系式来计算中频信号的频率。
在步骤S6中利用比值法估计离散中频信号x(n)频率:设离散中频信号为时宽有限信号,等效于无限时宽信号加长度为采样点数N的矩形窗函数表示如下:x(n)=w(n)·x′(n),其中n为离散中频信号的时间点,w(n)为矩形窗函数, x′(n)为无限时长离散中频信号,x′(n),n=1,2,…+∞,离散中频信号x(n)的采样点数为N时的离散傅里叶变换Xw(k)表示为: 其中W(k)为矩形窗函数w(n)的离散傅里叶变换,X′(k)为无限时宽离散中频信号x′(n)的离散傅里叶变换且X′(k)为一脉冲信号,*表示卷积,Xw(k)为加矩形窗后离散中频信号x(n)的离散傅里变换;矩形窗函数w(n)的离散傅里叶变换W(k)表示为:k=0~N-1,得知单目标的回波功率谱按sinc函数变化;因为N>>0,因此在主瓣区域频谱的模为:由此得到矩形窗频谱主瓣函数h(u):h(u)=sin(πu)/(πu),其中u表示连续自变量;
设比值v为: 由此有u·h(u)+(u+1)·h(u+1)=0,说明幅值谱主瓣内两条相邻谱线的重心为主瓣中心,当u=0时为两条谱线的重心,对应的频率为离散中频信号频率Fx;设离散中频信号离散频谱主瓣内最大幅值为yk,对应谱线号为KN,对应的频率为离散中频信号频率Fx所在位置与其离散傅里叶变换幅度最大处对应谱线位置的偏差为:KN的左右两侧相邻谱线幅值分别为和,对应的谱线号分别KN-1和KN+1;若为幅度次大值,则比值离散中频信号频率Fx所在位置与其离散傅里叶变换幅度最大处对应谱线位置的偏差为:从而得到离散中频信号频率估计值F为则有如下关系式:
在步骤4中减少采样点数的傅里叶变换,当减少的采样点数为M时,离散傅里叶变换运算点数为N-M,最大值的谱线号k′m=km,次大值的谱线号为:k′m-1=km-1,则最大值的谱线号k′m对应的频率为:次大值的谱线号k′m-1对应的频率为由比值法得到离散中频信号频率Fx所在位置应该在最大值谱线号和次大值谱线号之间,并且更靠近最大值谱线号,因此得到离散中频信号频率值属于区间ΔF1如下表示:
当减少的采样点数为M+1时离散频谱最大值发生跳变时离散频谱最大值的谱线号为k′m=km-1,此时的离散频谱次大值的谱线号为k′m+1=km,离散傅里叶变换运算的采样点数为N-M-1,则离散频谱最大值谱线号k′m对应的频率为次大值的谱线号k′m+1对应的频率为由比值法得到离散中频信号频率Fx所在位置应该在最大值谱线号和次大值谱线号之间,并且更靠近最大值谱线号,因此得到离散中频信号频率值属于区间ΔF2表示为:
离散中频信号频率Fx属于区间ΔF1和ΔF2的公共部分,且ΔF1∩ΔF2=(f1,f2),则 其中f1表示离散中频信号频率估计下限,f2表示离散中频信号频率估计上限,区间长度为
根据频率与距离的关系:计算距离区间;将频率区间中f1和f2代入距离与频率的关系式中替代f就能计算得到对应的距离区间为 其中R1表示距离估计下限,R2表示距离估计上限,区间长度为: 因此可以取作为距离测量的估计值,因此理想情况下的最大误差为δR/2。
测距正确性判断,根据非均匀频谱细化计算距离的公式对所测距离进行正确性判断。继续做减少点数的离散傅里叶变换,并且找到离散中频信号离散频谱最大值的谱线号为k″m。当减少的点数为M1时,离散频谱最大值的谱线号k″m=km-1,离散傅里叶变换运算的点数为N-M1,则最大值的谱线号k″m对应的频率为:次大值的谱线号k″m-1对应的频率为由比值法得到离散中频信号频率Fx所在位置应该在最大值谱线号和次大值谱线号之间,并且更靠近最大值谱线号,因此得到离散中频信号频率值属于区间ΔF3表示为:
当减少的采样点数为M1+1时离散频谱最大值发生跳变时离散频谱最大值的谱线号为k″m=km-2,此时的离散频谱次大值的谱线号为k″m+1=km-1,离散傅里叶变换运算的采样点数为N-M1-1,则离散频谱最大值谱线号k″m对应的频率为次大值的谱线号k″m+1对应的频率为由比值法得到离散中频信号频率Fx所在位置应该在最大值谱线号和次大值谱线号之间,并且更靠近最大值谱线号,因此得到离散中频信号频率值属于区间ΔF4如下表示:
离散中频信号频率Fx属于区间ΔF3和ΔF4的公共部分,且ΔF3∩ΔF4=(f3,f4),则 其中f3表示离散中频信号频率估计下限,f4表示离散中频信号频率估计上限,区间长度:
由于离散中频信号频率Fx应该落入(f1,f2)和(f3,f4)共同的区间内,若(f1,f2)和(f3,f4)有共同的区间则说明上述计算的距离值R′的正确性;若没有共同区间则说明中频信号受噪声和干涉的影响严重,应该舍弃这次测量,重新采样中频信号进行测量。
如图4示出本发明的非均匀频谱细化频谱最大值的谱线号变化图中作出了进行300次非均匀频谱细化过程中最大值谱线号的变换,从图中可以看出对中频信号做非均匀频谱细化时,最大值谱线号k′m呈现周期性的阶梯减小变化。在每次跳变时可以计算出一个中频信号频率的区间,因为在不同的频率区间内,每次非均匀频谱细化受到噪声影响不同,因此可以有效提高算法的抗造性能。
利用matlab仿真验证了本发明的理论性。构建离散中频信号如下:
测量距离R从5m到6m,步长为1mm,扫描周期T=100us,扫频带宽B=600MHz,光速c=3×108m/s,起始频率f1=24GHz,采样点数N=1024,采样频率Fs=N/T,w(n)为加入20dB信噪比的白噪声。
如图5示出本发明的仿真误差图,是利用非均匀频谱细化法计算存在20dB信噪比的白噪声时目标距离的误差分布图。从图中可以看出误差大部分在4mm以下,满足液位计系统测量精度在mm级误差的要求。
为了减少运算量,可以在第一次出现最大值频谱号跳变时DFT的采样点数为N-M-1,设下一次最大值频谱号出现跳变时的DFT采样点数为N-M1-1,则
(f1,f2)和(f3,f4)应该有共同区间,因此
这样就可以计算出M1的大小,就不用计算整个非均匀频谱细化一个周期内的减少点数的离散傅里叶变换。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种调频连续波雷达系统测距方法,该方法包括步骤如下:
步骤S1:确定频谱主瓣内的最大值谱线号和次大值谱线号,对调频连续波雷达系统采集的离散中频信号进行离散傅里叶变换,找到离散频谱最大值谱线号km和离散频谱次大值谱线号kc;
步骤S2:对谱线号km和kc进行判断,若谱线号kc≠km-1,则进入步骤S3;若谱线号kc=km-1,则进入步骤S4;
步骤S3:对离散中频信号的频域向右进行Δf/2的移位,使移位后的离散频谱最大值谱线号km和离散频谱次大值谱线号kc满足kc=km-1,其中Δf为N点离散傅里叶变换的频率分辨率;
步骤S4:对离散中频信号减少一个采样点数的离散傅里叶变换,找到离散频谱最大值谱线号k′m;
步骤S5:对离散频谱最大值谱线号k′m进行判断,若谱线号k′m=km,则重复步骤S4,若谱线号k′m≠km,则谱线号k′m=km-1,进入步骤S6;
步骤S6:根据步骤S4中进行离散傅里叶变换的采样点数,再结合距离关系式计算出所测雷达与目标的距离区间(R1,R2),则所测雷达与目标的测距离值为
5.根据权利要求1所述的调频连续波雷达系统测距方法,其特征在于,对离散中频信号x(n)减少采样点数的离散傅里叶变换并取模值,根据离散中频信号离散频谱从而找到离散中频信号离散频谱X(k)最大值的谱线号为k′m;由于进行离散傅里叶变换的采样点数在减小,而采样率不变,频率间隔在逐渐增大,因此将这种减小采样点数的离散傅里叶变换称作非均匀频谱细化,其中k为离散频谱谱线号,,N为采样点数,i为进行减小采样点数离散傅里叶变换时减少的采样点数,n为离散中频信号的离散时间点,n=0,1,2,…N-i-1,旋转因子WN-i+1满足欧拉公式:
e为自然对数的底数,j为虚数单位,|·|表示取复数的模值。
6.根据权利要求1所述的调频连续波雷达系统测距方法,其特征在于,所述距离区间(R1,R2)的计算步骤包括:
步骤S51:在减少采样点数的傅里叶变换中,当减少的采样点数为M时,离散傅里叶变换运算点数为N-M,最大值的谱线号k′m=km,次大值的谱线号为:k′m-1=km-1,则最大值的谱线号k′m对应的频率为次大值的谱线号k′m-1对应的频率为由比值法得到离散中频信号频率Fx所在位置应该在最大值谱线号和次大值谱线号之间,并且更靠近最大值谱线号,因此得到离散中频信号频率值属于区间ΔF1如下表示:
步骤S52:当减少的采样点数为M+1时离散频谱最大值发生跳变时离散频谱最大值的谱线号为k′m=k′m-1,此时的离散频谱次大值的谱线号为k′m+1=km,离散傅里叶变换运算的采样点数为N-M-1,则离散频谱最大值谱线号k′m对应的频率为次大值的谱线号k′m+1对应的频率为由比值法得到离散中频信号频率Fx所在位置应该在最大值谱线号和次大值谱线号之间,并且更靠近最大值谱线号,因此得到离散中频信号频率值属于区间ΔF2如下表示:
步骤S53:结合步骤S51和步骤S52得到的区间ΔF1和ΔF2,离散中频信号频率Fx属于区间ΔF1和ΔF2的公共部分,且ΔF1∩ΔF2=(f1,f2),则Fx∈(f1,f2),其中f1表示离散中频信号频率估计下限,f2表示离散中频信号频率估计上限;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410084559.7A CN103823216B (zh) | 2014-03-10 | 2014-03-10 | 一种调频连续波雷达系统测距方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410084559.7A CN103823216B (zh) | 2014-03-10 | 2014-03-10 | 一种调频连续波雷达系统测距方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103823216A true CN103823216A (zh) | 2014-05-28 |
CN103823216B CN103823216B (zh) | 2016-03-02 |
Family
ID=50758367
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410084559.7A Active CN103823216B (zh) | 2014-03-10 | 2014-03-10 | 一种调频连续波雷达系统测距方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103823216B (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105182323A (zh) * | 2015-09-15 | 2015-12-23 | 奇瑞汽车股份有限公司 | 一种fmcw雷达测距系统 |
CN104569961B (zh) * | 2015-01-22 | 2017-04-26 | 中国科学院电子学研究所 | 基于频谱细化的雷达测距方法 |
CN107153189A (zh) * | 2017-04-18 | 2017-09-12 | 上海交通大学 | 线性调频连续波雷达测距的信号处理方法 |
CN107907878A (zh) * | 2017-11-08 | 2018-04-13 | 零八电子集团有限公司 | 高精度获取fmcw雷达距离测量值的方法 |
CN108088406A (zh) * | 2018-01-09 | 2018-05-29 | 中国航空工业集团公司北京长城计量测试技术研究所 | 一种激波测速用压力传感器间有效距离测量方法 |
CN109164441A (zh) * | 2018-09-30 | 2019-01-08 | 厦门大学 | 一种雷达测距的方法 |
CN109471095A (zh) * | 2018-11-06 | 2019-03-15 | 哈尔滨工程大学 | 基于快速迭代插值的fmcw雷达距离估计算法 |
CN109814102A (zh) * | 2019-01-31 | 2019-05-28 | 厦门精益远达智能科技有限公司 | 一种单车道的超高监控方法、装置、设备和存储介质 |
CN113204005A (zh) * | 2021-05-07 | 2021-08-03 | 北京理工大学 | 一种提高调频连续波激光雷达距离解算精度的方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1566985A (zh) * | 2003-06-30 | 2005-01-19 | 武汉大学 | 线性调频连续波体制运动目标参数估计方法 |
CN101295019A (zh) * | 2008-05-23 | 2008-10-29 | 北京航空航天大学 | 一种线性调频连续波合成孔径雷达视频信号的仿真方法 |
WO2012103557A2 (en) * | 2011-01-28 | 2012-08-02 | The Regents Of The University Of Colorado, A Body Corporate | Spectral phase analysis for precision ranging |
-
2014
- 2014-03-10 CN CN201410084559.7A patent/CN103823216B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1566985A (zh) * | 2003-06-30 | 2005-01-19 | 武汉大学 | 线性调频连续波体制运动目标参数估计方法 |
CN101295019A (zh) * | 2008-05-23 | 2008-10-29 | 北京航空航天大学 | 一种线性调频连续波合成孔径雷达视频信号的仿真方法 |
WO2012103557A2 (en) * | 2011-01-28 | 2012-08-02 | The Regents Of The University Of Colorado, A Body Corporate | Spectral phase analysis for precision ranging |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104569961B (zh) * | 2015-01-22 | 2017-04-26 | 中国科学院电子学研究所 | 基于频谱细化的雷达测距方法 |
CN105182323A (zh) * | 2015-09-15 | 2015-12-23 | 奇瑞汽车股份有限公司 | 一种fmcw雷达测距系统 |
CN107153189A (zh) * | 2017-04-18 | 2017-09-12 | 上海交通大学 | 线性调频连续波雷达测距的信号处理方法 |
CN107153189B (zh) * | 2017-04-18 | 2021-08-03 | 上海交通大学 | 线性调频连续波雷达测距的信号处理方法 |
CN107907878A (zh) * | 2017-11-08 | 2018-04-13 | 零八电子集团有限公司 | 高精度获取fmcw雷达距离测量值的方法 |
CN107907878B (zh) * | 2017-11-08 | 2021-05-11 | 零八一电子集团有限公司 | 高精度获取fmcw雷达距离测量值的方法 |
CN108088406B (zh) * | 2018-01-09 | 2019-09-10 | 中国航空工业集团公司北京长城计量测试技术研究所 | 一种激波测速用压力传感器间有效距离测量方法 |
CN108088406A (zh) * | 2018-01-09 | 2018-05-29 | 中国航空工业集团公司北京长城计量测试技术研究所 | 一种激波测速用压力传感器间有效距离测量方法 |
CN109164441B (zh) * | 2018-09-30 | 2020-11-17 | 厦门大学 | 一种雷达测距的方法 |
CN109164441A (zh) * | 2018-09-30 | 2019-01-08 | 厦门大学 | 一种雷达测距的方法 |
CN109471095A (zh) * | 2018-11-06 | 2019-03-15 | 哈尔滨工程大学 | 基于快速迭代插值的fmcw雷达距离估计算法 |
CN109471095B (zh) * | 2018-11-06 | 2023-02-14 | 哈尔滨工程大学 | 一种基于快速迭代插值的fmcw雷达距离估计方法 |
CN109814102A (zh) * | 2019-01-31 | 2019-05-28 | 厦门精益远达智能科技有限公司 | 一种单车道的超高监控方法、装置、设备和存储介质 |
CN113204005A (zh) * | 2021-05-07 | 2021-08-03 | 北京理工大学 | 一种提高调频连续波激光雷达距离解算精度的方法及装置 |
CN113204005B (zh) * | 2021-05-07 | 2023-05-02 | 北京理工大学 | 一种提高调频连续波激光雷达距离解算精度的方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN103823216B (zh) | 2016-03-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103823216B (zh) | 一种调频连续波雷达系统测距方法 | |
CN103823215B (zh) | 线性调频连续波雷达测距方法 | |
CN103777199B (zh) | 一种调频连续波雷达系统的测距方法 | |
CN101825660B (zh) | 欠采样下的正弦信号频率的高效测量方法及实施装置 | |
CN107907878A (zh) | 高精度获取fmcw雷达距离测量值的方法 | |
JP2006516736A (ja) | 速度が曖昧な物標の非コヒーレント積分のための技法 | |
CN103364783B (zh) | 基于单通道sar的运动目标径向速度无模糊估计方法 | |
CN104391282A (zh) | 通过频谱校正来改善成像质量的方法 | |
CN102353946B (zh) | 一种基于x波段雷达图像的海表面流反演方法 | |
CN105229431A (zh) | 具有改进的距离确定的料位测量 | |
US20110239781A1 (en) | Open Channel Meter for Measuring Velocity | |
US11255956B2 (en) | Method for measuring the flow velocity of a medium | |
CN105738889A (zh) | 一种调频连续波测速测距方法 | |
CN103941087A (zh) | 欠采样速率下的高频余弦信号的频率测量方法及其装置 | |
CN103983957B (zh) | 一种多普勒偏移测量方法及其装置 | |
CN104237861A (zh) | 一种未知杂波背景下的cfar检测门限获取方法 | |
CN104076200A (zh) | 一种双子段相位差频率估计方法及其装置 | |
CN105629224A (zh) | 调频连续波雷达高精度测距方法 | |
CN111474524A (zh) | 一种雷达干扰装备干扰效果监测与决策支持系统 | |
CN104678392A (zh) | 一种基于稳定像元加权的地基sar大气延迟改正方法 | |
CN102043091A (zh) | 数字化高精度相位检测器 | |
CN105204011A (zh) | 一种高频地波雷达形成浪场的方法 | |
CN106443614B (zh) | 高超声速目标测加速度方法 | |
CN104977583A (zh) | 一种基于经验正交分解的x波段雷达海浪反演方法 | |
CN108445477A (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 |