CN107153189B - 线性调频连续波雷达测距的信号处理方法 - Google Patents
线性调频连续波雷达测距的信号处理方法 Download PDFInfo
- Publication number
- CN107153189B CN107153189B CN201710254182.9A CN201710254182A CN107153189B CN 107153189 B CN107153189 B CN 107153189B CN 201710254182 A CN201710254182 A CN 201710254182A CN 107153189 B CN107153189 B CN 107153189B
- Authority
- CN
- China
- Prior art keywords
- signal
- frequency
- component
- radar
- fourier series
- 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.)
- Active
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
- G01S13/34—Systems for measuring distance only using transmission of continuous waves, whether amplitude-, frequency-, or phase-modulated, or unmodulated using transmission of continuous, frequency-modulated waves while heterodyning the received signal, or a signal derived therefrom, with a locally-generated signal related to the contemporaneously transmitted signal
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)
- Medicines Containing Antibodies Or Antigens For Use As Internal Diagnostic Agents (AREA)
Abstract
Description
技术领域
本发明涉及雷达测距技术领域,具体地,涉及线性调频连续波雷达测距的信号处理方法。
背景技术
雷达测距是一种非接触式测量技术,可在恶劣环境下稳定工作,不易受雨雪等天气条件的影响,并具有穿透烟雾和灰尘的能力。其中,线性调频连续波(Linear FrequencyModulation Continuous Wave,简称LFMCW)雷达则具有辐射功率小、结构紧凑以及成本低廉等优点,在民用及军事领域广泛使用。其原理是:向探测目标发射线性调频的连续波并接收目标散射的电磁回波,将回波与本振信号混频得到差拍信号,进而提取出目标的距离信息。
常规地,通过差拍-离散傅里叶变换在频域中寻找频谱峰值所在的频率点,以此计算出目标距雷达的距离。但是离散傅里叶变换存在栅栏效应并受限于频率分辨率,无法精确地提取出差拍信号的真实频率。为了减小栅栏效应带来的误差,可以通过补零或者线性调频Z变换(CZT)等手段来细化频谱,从而找到比较精确的频率值。然而这些方法只能初步提高频率估计精度,而且容易受到噪声的干扰。因此,如何进一步提高线性调频连续波雷达的测距精度仍是一个亟待解决的实际问题。此外,在雷达多目标测距应用中如何自适应的快速提取出多目标的精确距离信息,更是对差拍信号的处理提出了更高的要求。
发明内容
针对现有技术中的缺陷,本发明的目的是提供一种线性调频连续波雷达测距的信号处理方法。
根据本发明提供的线性调频连续波雷达测距的信号处理方法,包括如下步骤:
步骤1:对差拍信号s(t)进行离散采样,得到离散差拍信号s(n);
步骤3:由分量个数K判定是否为单分量信号,若是,则跳到步骤4执行,若否,则对差拍信号进行信号分解,得到K个独立分量,分别记为s1,…,si,…sK;
步骤4:对各个分量依次进行精确的频率参数估计;
步骤5:由雷达差拍信号各个分量的频率估计值计算目标距雷达的距离。
优选地,所述步骤2中通过频谱幅值比较得到分量个数K以及各分量的初始化频率是指:
对离散频谱F(n)进行取模处理,将取模结果记为序列 {|F(0)|,|F(1)|,…,|F(N-1)|},判定序列{|F(0)|,|F(1)|,…,|F(N-1)|}的极大值点,并筛选出大于或等于特定阈值的极大值点,将筛选出的极大值点的个数即分量的个数记为K,筛选出的极大值点的序号对应的频率为初始化频率,具体为:
优选地,步骤2中所述特定阈值的确定准则如下:
根据实际距离探测的先验知识确定雷达发射波覆盖的处于不同距离分布的可识别目标物体的个数,将所述个数作为分量个数K,取极大值点中按值从大到小排列的前K 个极大值点,将此K个极大值点对应的极大值中的最小值设置为阈值;
当先验知识不明确时,根据实际测距中目标分辨能力的要求,设定阈值为序列 {|F(0)|,|F(1)|,…,|F(N-1)|}平均值的多倍值。
优选地,所述步骤3中对差拍信号进行信号分解的方法包括:
差拍信号s(t)展开形式如下:
式中:表示傅里叶级数直流分量,表示傅里叶级数m次余弦谐波函数的系数,F0表示傅里叶级数的基频,表示傅里叶级数m次正弦谐波函数的系数;表示傅里叶级数直流分量,表示傅里叶级数m次余弦谐波函数的系数,表示傅里叶级数m次正弦谐波函数的系数,其中F0=fs/(LN),L为正整数,综合考虑测距精度与拟合计算量的要求确定L的取值,P为傅里叶级数的阶次;
由傅里叶级数拟合的系数,分别取如下向量:
分别取如下矩阵:
t0,…,tN-1分别为差拍信号s(n)对应的时间序列;则差拍信号s(t)用矩阵形式表示为:Hi表示针对第i个分量定义的矩阵,Ci表示Hi的子矩阵, Di表示Hi的子矩阵,i的取值范围为1,2,…K;
y=(HTH+λI)-1HTS;
式中:I为与HTH同尺寸的单位矩阵,λ为Tikhonov正则化参数;
步骤3.3:由求得的向量y得到各子向量y1,…,yK的值,计算得到各个信号分量s1,…,sK,计算公式为:si=Hiyi。
优选地,所述步骤4包括:
步骤4.1:根据测距精度要求,确定细化频率分辨率,如下式所示:
式中:△r为所需的距离分辨精度,c为电磁波在介质中的传播速度,B为雷达信号源扫频带宽,N为采样点数;
步骤4.3:利用希尔伯特变换或者结合正交双通道的差拍信号,求得各分量差拍信号的复数解析信号,记为Si(n)。其中,Si(n)表示第i个分量的复数解析信号,i的取值范围为1,…,K;
步骤4.4,求取各分量的频率估计值;具体地,
通过解调算子e-j2πft寻找取得频谱集中性最好的频率参数值,如下式所示:
优选地,所述步骤5包括:由雷达差拍信号各个分量的频率估计值,计算目标距雷达的距离,如下式所示:
式中:ri表示第i个目标物体距雷达的距离,c为电磁波在介质中的传播速度,Tswp为雷达信号源扫频时间,B为雷达发射波的带宽。
与现有技术相比,本发明具有如下的有益效果:
本发明提供的线性调频连续波雷达测距的信号处理方法,能够自适应的自动分离多个目标分量,并显著提高测距精度。基于傅里叶级数幅值估计拟合的信号分解方法和基于降调制的频率参数细化估计方法均具有较高的运算效率和良好的抗噪性能,满足了LFMCW雷达现场测距的实时性和抗干扰能力要求。
附图说明
通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:
图1为本发明提出的线性调频连续波雷达测距的信号处理方法的流程图;
图2为本发明实施例中仿真的差拍信号的时域波形图;
图3为本发明实施例中仿真的差拍信号的离散频谱图;
图4为本发明实施例中差拍信号的信号分解结果的时域波形图;
图5为本发明实施例中基于频谱集中性指标进行频率估计的结果对比图。
具体实施方式
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明,但不以任何形式限制本发明。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变化和改进。这些都属于本发明的保护范围。
如图1所示,本发明提供的线性调频连续波雷达测距的信号处理方法,包括如下步骤:
步骤1:对差拍信号s(t)进行离散采样,得到离散差拍信号s(n);
在实际雷达测距过程中,差拍信号难免会受到扫频调制信号的干扰,因此可采用触发采集的方式,每次均采集扫频单周期信号,而对于连续多周期采集的信号可以采用高通滤波等处理方法滤除扫频调制信号的干扰。如图2所示为仿真的差拍信号时域波形图,其中采样频率fs=40k,采样时间为4ms,差拍信号s(n)如下:
s(n)=3cos(2π*2182n/fs+4π/3)+1.4cos(2π*2580n/fs+π/8)+2cos(2π*4185n/fs+π/3)+e(n/fs),
其中,e(n/fs)为信噪比为10dB的高斯白噪声。
步骤2:对所述的离散差拍信号s(n)进行离散傅里叶变换得到频谱F(n),对离散频谱F(n)进行取模处理,将取模结果记为序列{|F(0)|,|F(1)|,…,|F(N-1)|},判定序列{|F(0)|,|F(1)|,…,|F(N-1)|}的极大值点,并筛选出大于或等于特定阈值的极大值点,将筛选出的极大值点的个数即分量的个数记为K,筛选出的极大值点的序号对应的频率为初始化频率,具体为:
步骤3:由分量个数K判定是否为单分量信号,若是,则跳到步骤4执行;若否,则对差拍信号进行信号分解,得到K个独立分量,分别记为s1,…,sK,信号分解的方法具体包括:
式中:表示傅里叶级数直流分量,表示傅里叶级数m次余弦谐波函数的系数,F0表示傅里叶级数的基频,表示傅里叶级数m次正弦谐波函数的系数;表示傅里叶级数直流分量,表示傅里叶级数m次余弦谐波函数的系数,表示傅里叶级数m次正弦谐波函数的系数,其中F0=fs/(LN),L为正整数,可综合考虑测距精度与拟合计算量的要求确定,一般L的取值应使F0为细化频率分辨率的2到6倍左右, P为傅里叶级数的阶次,一般可取P=L/2或P=L。
由傅里叶级数拟合的系数,取向量:
取矩阵:
其中:
t0,…,tN-1为差拍信号s(n)对应的时间序列。
因此差拍信号s(t)用矩阵形式表示为:
y=(HTH+λI)-1HTS
式中:I为与HTH同尺寸的单位矩阵,λ为Tikhonov正则化参数。
步骤3.3:由求得的向量y得到各子向量y1,…,yK的值,计算得到各个信号分量s1,…,sK,如下式所示:
si=Hiyi
在本实施例中K=3,因此存在多目标分量信号,需要进行信号分解,得到s1,s2和s3。取L=32,P=32,λ=0.2,分解得到的三个信号分量s1,s2和s3如图4所示。
步骤4:对各个分量依次进行精确的频率参数估计,对每一个分量进行频率参数估计的方法具体包括:
步骤4.1:根据测距精度要求,确定细化频率分辨率,如下式所示:
式中:△r为所需的距离分辨精度,c为电磁波在介质中的传播速度,B为雷达信号源扫频带宽,N为采样点数。
步骤4.3:利用希尔伯特变换或者结合正交双通道的差拍信号,求得各分量差拍信号的复数解析信号Si(n),其中Si(n)表示第i个分量的复数解析信号,i的取值范围为 1,…,K。
步骤4.4,求取各分量的频率估计值;
通过解调算子e-j2πft寻找取得降调制性能最好的频率参数值,如下式所示:
在本实施例中取△f=4Hz,如图5所示为针对第二个分量求解得到的abs[X(k)|k=0] 与估计频率f的对应关系图,比较得到fest_2=2567Hz。同样的方法求解得到fest_1=2183 Hz,fest_3=4187Hz,因此可以看出3个分量的频率估计精度均得到了显著提高。
步骤5:由雷达差拍信号各个分量的频率估计值,计算目标距雷达的距离,如下式所示:
式中:ri表示第i个目标物体距雷达的距离,c为电磁波在介质中的传播速度,Tswp为雷达信号源扫频时间,B为雷达发射波带宽。
以上对本发明的具体实施例进行了描述。需要理解的是,本发明并不局限于上述特定实施方式,本领域技术人员可以在权利要求的范围内做出各种变化或修改,这并不影响本发明的实质内容。在不冲突的情况下,本申请的实施例和实施例中的特征可以任意相互组合。
Claims (5)
1.一种线性调频连续波雷达测距的信号处理方法,其特征在于,包括如下步骤:
步骤1:对差拍信号s(t)进行离散采样,得到离散差拍信号s(n);
步骤3:由分量个数K判定是否为单分量信号,若是,则跳到步骤4执行,若否,则对差拍信号进行信号分解,得到K个独立分量,分别记为s1,…,si,…sK;
步骤4:对各个分量依次进行精确的频率参数估计;
步骤5:由雷达差拍信号各个分量的频率估计值计算目标距雷达的距离;
所述步骤4包括:
步骤4.1:根据测距精度要求,确定细化频率分辨率,如下式所示:
式中:Δr为所需的距离分辨精度,c为电磁波在介质中的传播速度,B为雷达信号源扫频带宽,N为采样点数;
步骤4.3:利用希尔伯特变换或者结合正交双通道的差拍信号,求得各分量差拍信号的复数解析信号,记为Si(n);其中,Si(n)表示第i个分量的复数解析信号,i的取值范围为1,…,K;
步骤4.4,求取各分量的频率估计值;具体地,
通过解调算子e-j2πft寻找取得频谱集中性最好的频率参数值,如下式所示:
2.根据权利要求1所述的线性调频连续波雷达测距的信号处理方法,其特征在于,所述步骤2中通过频谱幅值比较得到分量个数K以及各分量的初始化频率是指:
对离散频谱F(n)进行取模处理,将取模结果记为序列{|F(0)|,|F(1)|,…,|F(N-1)|},判定序列{|F(0)|,|F(1)|,…,|F(N-1)|}的极大值点,并筛选出大于或等于特定阈值的极大值点,将筛选出的极大值点的个数即分量的个数记为K,筛选出的极大值点的序号对应的频率为初始化频率,具体为:
3.根据权利要求2所述的线性调频连续波雷达测距的信号处理方法,其特征在于,步骤2中所述特定阈值的确定准则如下:
根据实际距离探测的先验知识确定雷达发射波覆盖的处于不同距离分布的可识别目标物体的个数,将所述个数作为分量个数K,取极大值点中按值从大到小排列的前K个极大值点,将此K个极大值点对应的极大值中的最小值设置为阈值;
当先验知识不明确时,根据实际测距中目标分辨能力的要求,设定阈值为序列{|F(0)|,|F(1)|,…,|F(N-1)|}平均值的多倍值。
4.根据权利要求1所述的线性调频连续波雷达测距的信号处理方法,其特征在于,所述步骤3中对差拍信号进行信号分解的方法包括:
差拍信号s(t)展开形式如下:
式中:表示傅里叶级数直流分量,表示傅里叶级数m次余弦谐波函数的系数,F0表示傅里叶级数的基频,表示傅里叶级数m次正弦谐波函数的系数;表示傅里叶级数直流分量,表示傅里叶级数m次余弦谐波函数的系数,表示傅里叶级数m次正弦谐波函数的系数,其中F0=fs/(LN),fs表示采样频率,L为正整数,N为采样点数,综合考虑测距精度与拟合计算量的要求确定L的取值,P为傅里叶级数的阶次;
由傅里叶级数拟合的系数,分别取如下向量:
分别取如下矩阵:
t0,…,tN-1分别为离散差拍信号s(n)对应的时间序列;则差拍信号s(t)用矩阵形式表示为:Hi表示针对第i个分量定义的矩阵,Ci表示Hi的子矩阵,Di表示Hi的子矩阵,i的取值范围为1,2,…K;
y=(HTH+λI)-1HTS;
式中:I为与HTH同尺寸的单位矩阵,λ为Tikhonov正则化参数;
步骤3.3:由求得的向量y得到各子向量y1,…,yK的值,计算得到各个信号分量s1,…,sK,计算公式为:si=Hiyi。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710254182.9A CN107153189B (zh) | 2017-04-18 | 2017-04-18 | 线性调频连续波雷达测距的信号处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710254182.9A CN107153189B (zh) | 2017-04-18 | 2017-04-18 | 线性调频连续波雷达测距的信号处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107153189A CN107153189A (zh) | 2017-09-12 |
CN107153189B true CN107153189B (zh) | 2021-08-03 |
Family
ID=59793033
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710254182.9A Active CN107153189B (zh) | 2017-04-18 | 2017-04-18 | 线性调频连续波雷达测距的信号处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107153189B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107632292B (zh) * | 2017-09-21 | 2021-07-30 | 北京工业大学 | 一种对雷达信号进行调频傅立叶变换的方法 |
CN110244290A (zh) * | 2019-06-17 | 2019-09-17 | 电子科技大学 | 一种距离扩展目标的检测方法 |
CN113204018B (zh) * | 2021-03-16 | 2022-11-04 | 西安电子科技大学 | 一种基于lfmcw的mimo雷达快速信号处理方法 |
CN113238244B (zh) * | 2021-05-24 | 2022-05-27 | 哈尔滨理工大学 | 一种fmcw激光测距差拍信号频率估计方法及*** |
CN115166681B (zh) * | 2022-09-07 | 2022-12-06 | 武汉新朗光电科技有限公司 | 调频连续波信号体制穿墙雷达目标检测快速方法及*** |
CN115877356B (zh) * | 2023-02-13 | 2023-06-02 | 宁波飞芯电子科技有限公司 | 探测方法、探测设备及计算机可读存储介质 |
CN118050740B (zh) * | 2024-04-16 | 2024-06-14 | 深圳市圳阳精密技术有限公司 | 一种调频连续波激光雷达信号分析方法及*** |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7034743B2 (en) * | 2002-06-04 | 2006-04-25 | Murata Manufacturing Co., Ltd. | Radar |
CN1894595A (zh) * | 2003-12-16 | 2007-01-10 | 株式会社村田制作所 | 雷达 |
CN102707275A (zh) * | 2012-05-25 | 2012-10-03 | 北京航空航天大学 | 一种线性调频连续波雷达高度表的数字化处理方法 |
CN103777199A (zh) * | 2014-02-24 | 2014-05-07 | 中国科学院电子学研究所 | 一种调频连续波雷达***的测距方法 |
CN103823215A (zh) * | 2014-03-03 | 2014-05-28 | 中国科学院电子学研究所 | 线性调频连续波雷达测距方法 |
CN103823216A (zh) * | 2014-03-10 | 2014-05-28 | 中国科学院电子学研究所 | 一种调频连续波雷达***测距方法 |
CN104569961A (zh) * | 2015-01-22 | 2015-04-29 | 中国科学院电子学研究所 | 基于频谱细化的雷达测距方法 |
-
2017
- 2017-04-18 CN CN201710254182.9A patent/CN107153189B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7034743B2 (en) * | 2002-06-04 | 2006-04-25 | Murata Manufacturing Co., Ltd. | Radar |
CN1894595A (zh) * | 2003-12-16 | 2007-01-10 | 株式会社村田制作所 | 雷达 |
CN102707275A (zh) * | 2012-05-25 | 2012-10-03 | 北京航空航天大学 | 一种线性调频连续波雷达高度表的数字化处理方法 |
CN103777199A (zh) * | 2014-02-24 | 2014-05-07 | 中国科学院电子学研究所 | 一种调频连续波雷达***的测距方法 |
CN103823215A (zh) * | 2014-03-03 | 2014-05-28 | 中国科学院电子学研究所 | 线性调频连续波雷达测距方法 |
CN103823216A (zh) * | 2014-03-10 | 2014-05-28 | 中国科学院电子学研究所 | 一种调频连续波雷达***测距方法 |
CN104569961A (zh) * | 2015-01-22 | 2015-04-29 | 中国科学院电子学研究所 | 基于频谱细化的雷达测距方法 |
Non-Patent Citations (1)
Title |
---|
基于周期FRFT的多分量LFMCW雷达信号分离;黄宇等;《航空学报》;20130425;第34卷(第4期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN107153189A (zh) | 2017-09-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107153189B (zh) | 线性调频连续波雷达测距的信号处理方法 | |
CN111352102B (zh) | 一种基于调频连续波雷达的多目标个数检测方法及装置 | |
CN109061589B (zh) | 随机跳频雷达的目标运动参数估计方法 | |
CN107561508B (zh) | 一种用于匀加速运动目标的相参积累检测方法 | |
CN103630888B (zh) | 基于对称三角lfmcw雷达的高精度实时微波测速测距装置 | |
CN107907878B (zh) | 高精度获取fmcw雷达距离测量值的方法 | |
CN103744068B (zh) | 双通道调频连续波sar***的动目标检测成像方法 | |
CN107290589B (zh) | 基于短时分数阶傅里叶变换的非线性信号时频分析方法 | |
CN108761404A (zh) | 一种基于二次相位函数参数估计及补偿的改进算法 | |
CN106814353A (zh) | 一种雷达信号处理*** | |
CN103293521B (zh) | 一种利用x波段雷达探测近海海域水深的方法 | |
CN107255814B (zh) | 一种基于lfmsk波形的雷达目标检测方法 | |
CN110568434B (zh) | 一种多通道匀加速sar动目标二维速度估计方法 | |
CN108196241B (zh) | 一种基于Hough变换的高速动目标速度估计方法 | |
CN110187342B (zh) | 一种基于fmcw移动平台的生命体征检测与成像方法 | |
CN108693523A (zh) | 基于锯齿波线性调频雷达多目标测距测速的方法及*** | |
CN111580061B (zh) | 基于clean算法的电离层电子密度反演方法 | |
CN109521411A (zh) | 一种距离扩展目标的检测方法 | |
Wang et al. | Denoising FMCW ladar signals via EEMD with singular spectrum constraint | |
Molchanov et al. | On micro-Doppler period estimation | |
CN117838083A (zh) | 一种基于毫米波雷达的体征快速精确检测方法 | |
US10782391B2 (en) | Processing received radiation reflected from a target | |
CN109001671B (zh) | 一种跳频信号的目标检测和参数估计方法及装置 | |
CN112906476B (zh) | 一种基于信杂噪比损失的机载雷达训练样本选择方法 | |
CN115494496B (zh) | 单比特雷达成像***、方法及相关设备 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |