线性调频连续波雷达测距的信号处理方法技术

编号:201710254182 阅读:43 评论( 0 )

本发明专利技术提供了一种线性调频连续波雷达测距的信号处理方法,包括步骤1,对差拍信号s(t)采样后得到离散差拍信号s(n);步骤2,对s(n)进行离散傅里叶变换得到频谱F(n),通过频谱幅值比较得到分量个数K以及各分量的初始化频率

Signal processing method for linear frequency modulation continuous wave radar ranging

The present invention provides a method of signal processing for linear frequency modulated continuous wave radar ranging, including 1 steps, the beat signal s (T) obtained after sampling discrete beat signal s (n); step 2, the S (n) spectrum F discrete Fu Liye transform (n), through the spectrum amplitude comparison get the number of components and K components initialization frequency

微信分享本专利技术
微信扫一扫分享本专利技术到微信朋友圈,让更多的人了解这个技术。

全部详细技术资料下载

1 技术实现步骤摘要

线性调频连续波雷达测距的信号处理方法
本专利技术涉及雷达测距
,具体地,涉及线性调频连续波雷达测距的信号处理方法
技术介绍
雷达测距是一种非接触式测量技术,可在恶劣环境下稳定工作,不易受雨雪等天气条件的影响,并具有穿透烟雾和灰尘的能力。其中,线性调频连续波(LinearFrequencyModulationContinuousWave,简称LFMCW)雷达则具有辐射功率小、结构紧凑以及成本低廉等优点,在民用及军事领域广泛使用。其原理是:向探测目标发射线性调频的连续波并接收目标散射的电磁回波,将回波与本振信号混频得到差拍信号,进而提取出目标的距离信息。常规地,通过差拍-离散傅里叶变换在频域中寻找频谱峰值所在的频率点,以此计算出目标距雷达的距离。但是离散傅里叶变换存在栅栏效应并受限于频率分辨率,无法精确地提取出差拍信号的真实频率。为了减小栅栏效应带来的误差,可以通过补零或者线性调频Z变换(CZT)等手段来细化频谱,从而找到比较精确的频率值。然而这些方法只能初步提高频率估计精度,而且容易受到噪声的干扰。因此,如何进一步提高线性调频连续波雷达的测距精度仍是一个亟待解决的实际问题。此外,在雷达多目标测距应用中如何自适应的快速提取出多目标的精确距离信息,更是对差拍信号的处理提出了更高的要求。
技术实现思路
针对现有技术中的缺陷,本专利技术的目的是提供一种线性调频连续波雷达测距的信号处理方法。根据本专利技术提供的线性调频连续波雷达测距的信号处理方法,包括如下步骤:步骤1:对差拍信号s(t)进行离散采样,得到离散差拍信号s(n);步骤2:对所述的离散差拍信号s(n)进行离散傅里叶变换得到频谱F(n),通过频谱幅值比较得到分量个数K以及各分量的初始化频率,分别记为步骤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,筛选出的极大值点的序号对应的频率为初始化频率,具体为:式中:为第i个分量的初始化频率,m为筛选出的第i个极大值点在序列{|F(0)|,|F(1)|,…,|F(N-1)|}中的序号;N表示采样点数,fs表示采样频率。优选地,步骤2中所述特定阈值的确定准则如下:根据实际距离探测的先验知识确定雷达发射波覆盖的处于不同距离分布的可识别目标物体的个数,将所述个数作为分量个数K,取极大值点中按值从大到小排列的前K个极大值点,将此K个极大值点对应的极大值中的最小值设置为阈值;当先验知识不明确时,根据实际测距中目标分辨能力的要求,设定阈值为序列{|F(0)|,|F(1)|,…,|F(N-1)|}平均值的多倍值。优选地,所述步骤3中对差拍信号进行信号分解的方法包括:步骤3.1:根据差拍信号的数学模型和K个分量的初始化频率对差拍信号s(t)进行展开表示后得到系数矩阵与系数向量,并将s(t)用矩阵形式表示,其中K个分量的初始化频率分别记为具体地:差拍信号s(t)展开形式如下:式中:Ai为第i个信号分量的幅值,fi为第i个信号分量的频率,θi0为第i个信号分量的初始相位,e(t)为噪声和其他干扰分量,为第i个信号分量的初始化频率;t表示时间分量;对各分量的系数与进行傅里叶级数的拟合,如下式所示:式中:表示傅里叶级数直流分量,表示傅里叶级数m次余弦谐波函数的系数,F0表示傅里叶级数的基频,表示傅里叶级数m次正弦谐波函数的系数;表示傅里叶级数直流分量,表示傅里叶级数m次余弦谐波函数的系数,表示傅里叶级数m次正弦谐波函数的系数,其中F0=fs/(LN),L为正整数,综合考虑测距精度与拟合计算量的要求确定L的取值,P为傅里叶级数的阶次;由傅里叶级数拟合的系数,分别取如下向量:式中:表示拟合的傅里叶级数系数向量,表示拟合的傅里叶级数系数向量,yi表示傅里叶级数拟合的总系数向量;分别取如下矩阵:Hi=[CiDi],其中:t0,…,tN-1分别为差拍信号s(n)对应的时间序列;则差拍信号s(t)用矩阵形式表示为:Hi表示针对第i个分量定义的矩阵,Ci表示Hi的子矩阵,Di表示Hi的子矩阵,i的取值范围为1,2,…K;步骤3.2:取矩阵H=[H1,…HK],向量S=[s(t0)…s(tN-1)]T,则向量y可通过下式计算得出:y=(HTH+λI)-1HTS;式中:I为与HTH同尺寸的单位矩阵,λ为Tikhonov正则化参数;步骤3.3:由求得的向量y得到各子向量y1,…,yK的值,计算得到各个信号分量s1,…,sK,计算公式为:si=Hiyi。优选地,所述步骤4包括:步骤4.1:根据测距精度要求,确定细化频率分辨率,如下式所示:式中:△r为所需的距离分辨精度,c为电磁波在介质中的传播速度,B为雷达信号源扫频带宽,N为采样点数;步骤4.2:选定各分量频率估计的区间,即第i个分量的频率估计区间为:其中fs为采样频率;步骤4.3:利用希尔伯特变换或者结合正交双通道的差拍信号,求得各分量差拍信号的复数解析信号,记为Si(n)。其中,Si(n)表示第i个分量的复数解析信号,i的取值范围为1,…,K;步骤4.4,求取各分量的频率估计值;具体地,通过解调算子e-j2πft寻找取得频谱集中性最好的频率参数值,如下式所示:其中fest_i为第i个分量的频率估计值,表示寻找中取得最大值的参数f,abs(·)表示复数幅值运算,优选地,所述步骤5包括:由雷达差拍信号各个分量的频率估计值,计算目标距雷达的距离,如下式所示:式中:ri表示第i个目标物体距雷达的距离,c为电磁波在介质中的传播速度,Tswp为雷达信号源扫频时间,B为雷达发射波的带宽。与现有技术相比,本专利技术具有如下的有益效果:本专利技术提供的线性调频连续波雷达测距的信号处理方法,能够自适应的自动分离多个目标分量,并显著提高测距精度。基于傅里叶级数幅值估计拟合的信号分解方法和基于降调制的频率参数细化估计方法均具有较高的运算效率和良好的抗噪性能,满足了LFMCW雷达现场测距的实时性和抗干扰能力要求。附图说明通过阅读参照以下附图对非限制性实施例所作的详细描述,本专利技术的其它特征、目的和优点将会变得更明显:图1为本专利技术提出的线性调频连续波雷达测距的信号处理方法的流程图;图2为本专利技术实施例中仿真的差拍信号的时域波形图;图3为本专利技术实施例中仿真的差拍信号的离散频谱图;图4为本专利技术实施例中差拍信号的信号分解结果的时域波形图;图5为本专利技术实施例中基于频谱集中性指标进行频率估计的结果对比图。具体实施方式下面结合具体实施例对本专利技术进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本专利技术,但不以任何形式限制本专利技术。应当指出的是,对本领域的普通技术人员来说,在不脱离本专利技术构思的前提下,还可以做出若干本文档来自技高网...
<a href="http://www.xjishu.com/zhuanli/52/201710254182.html" title="线性调频连续波雷达测距的信号处理方法原文来自X技术">线性调频连续波雷达测距的信号处理方法</a>
【详细说明在详细技术资料中】

全部详细技术资料下载 我是这个专利的主人

2 技术保护点

一种线性调频连续波雷达测距的信号处理方法,其特征在于,包括如下步骤:步骤1:对差拍信号s(t)进行离散采样,得到离散差拍信号s(n);步骤2:对所述的离散差拍信号s(n)进行离散傅里叶变换得到频谱F(n),通过频谱幅值比较得到分量个数K以及各分量的初始化频率,分别记为

3 技术保护范围摘要

1.一种线性调频连续波雷达测距的信号处理方法,其特征在于,包括如下步骤:步骤1:对差拍信号s(t)进行离散采样,得到离散差拍信号s(n);步骤2:对所述的离散差拍信号s(n)进行离散傅里叶变换得到频谱F(n),通过频谱幅值比较得到分量个数K以及各分量的初始化频率,分别记为步骤3:由分量个数K判定是否为单分量信号,若是,则跳到步骤4执行,若否,则对差拍信号进行信号分解,得到K个独立分量,分别记为s1,…,si,…sK;步骤4:对各个分量依次进行精确的频率参数估计;步骤5:由雷达差拍信号各个分量的频率估计值计算目标距雷达的距离。2.根据权利要求1所述的线性调频连续波雷达测距的信号处理方法,其特征在于,所述步骤2中通过频谱幅值比较得到分量个数K以及各分量的初始化频率是指:对离散频谱F(n)进行取模处理,将取模结果记为序列{|F(0)|,|F(1)|,…,|F(N-1)|},判定序列{|F(0)|,|F(1)|,…,|F(N-1)|}的极大值点,并筛选出大于或等于特定阈值的极大值点,将筛选出的极大值点的个数即分量的个数记为K,筛选出的极大值点的序号对应的频率为初始化频率,具体为:式中:为第i个分量的初始化频率,m为筛选出的第i个极大值点在序列{|F(0)|,|F(1)|,…,|F(N-1)|}中的序号;N表示采样点数,fs表示采样频率。3.根据权利要求2所述的线性调频连续波雷达测距的信号处理方法,其特征在于,步骤2中所述特定阈值的确定准则如下:根据实际距离探测的先验知识确定雷达发射波覆盖的处于不同距离分布的可识别目标物体的个数,将所述个数作为分量个数K,取极大值点中按值从大到小排列的前K个极大值点,将此K个极大值点对应的极大值中的最小值设置为阈值;当先验知识不明确时,根据实际测距中目标分辨能力的要求,设定阈值为序列{|F(0)|,|F(1)|,…,|F(N-1)|}平均值的多倍值。4.根据权利要求1所述的线性调频连续波雷达测距的信号处理方法,其特征在于,所述步骤3中对差拍信号进行信号分解的方法包括:步骤3.1:根据差拍信号的数学模型和K个分量的初始化频率对差拍信号s(t)进行展开表示后得到系数矩阵与系数向量,并将s(t)用矩阵形式表示,其中K个分量的初始化频率分别记为具体地:差拍信号s(t)展开形式如下:式中:Ai为第i个信号分量的幅值,fi为第i个信号分量的频率,θi0为第i个信号分量的初始相位,e(t)为噪声和其他干扰分量,为第i个信号分量的初始化频率;t表示时间分量;对各分量的系数与进行傅里叶级数的拟合,如下式所示:

4 专利技术属性

发明(设计)人:彭志科熊玉勇
申请(专利权)人:上海交通大学
专利类型:发明
专利号:201710254182
国别省市:上海,31

5 专利技术项目评估

还没有人评估本专利项目,你来评估一下?

对专利感兴趣?请拨打022-23869559(工作日)13820821600(24小时)咨询

6 相关技术资料

线性调频连续波雷达专利,雷达线性调频信号专利,线性调频雷达专利,信号处理中的线性代数专利,非线性信号处理专利,雷达信号处理基础专利,雷达信号处理专利,mimo雷达信号处理

测距雷达采样信号处理方法线性调频

网友询问留言: 已有 0 条评论

还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1