一种FMCW雷达的差频频率估计方法技术

本发明专利技术提供一种FMCW雷达的差频频率估计方法，应用于信号处理技术领域，包括如下步骤：获取待处理实信号，并将所述待处理实信号转换为解析信号；获取所述解析信号的实际相位；对所述实际相位在时域按照预设值进行平均分段，并对每一段进行最小二乘法相位拟合，获取每一段拟合直线的斜率，将所述拟合直线的斜率作为每一段的差频频率；去除因为相位提取以及噪声干扰偏差较大的部分段的差频频率，对剩余部分段的差频频率取均值，并根据该均值计算出FMCW雷达的实际差频频率。本发明专利技术提供的技术方案可以在非线性干扰的情况下，实现对调频连续波雷达的差频频率的精确估计，能够在降低运算量的同时，提高差频的估计精度。提高差频的估计精度。提高差频的估计精度。

【技术实现步骤摘要】
一种FMCW雷达的差频频率估计方法


[0001]本专利技术涉及信号处理
，特别涉及一种FMCW雷达的差频频率估计方法。

技术介绍

[0002]与传统的脉冲体制的雷达相比，调频连续波雷达具有更高的距离分辨率、更强的穿透能力、能够消除距离盲区等优势，因此，调频连续波雷达在高精度测距和成像领域得到广泛应用。FMCW雷达通过发射线性调频信号，并通过将接受信号与发射信号混频得到具有目标位置信息的差频信号，因此高精度位置估计问题就变成了频谱分量(频率和相位)估计问题。影响调频连续波雷达测距成像精度最关键的因素是线性调频信号的线性度，实际上，由于雷达发射机不可能产生理想的线性调频信号，产生的调频信号始终存在一定的非线性，这将会影响测距或者成像精度。调频信号当中的非线性分量，会对差频信号频谱产生一定的影响，例如主瓣展宽、非对称畸变、旁瓣能量过高等问题，这将影响目标差频频率的测量和目标位置的确定。
[0003]大多数的频谱估计方法假设的是一个时间连续信号，并且只受到高斯白噪声的影响，并没有考虑到调频信号非线性对频谱的影响。其中，现有技术中，使用最多的是快速傅里叶变换(FFT)方法，其最主要的缺点是当需要一个较高的频率分辨率时，需要一个较长的窗来进行处理。这可以通过补零实现，但是计算量较大。如果只需要分析信号频谱的一部分即差频频率附近时，可以使用chirp
‑
z变换(CZT)算法，该算法可以在减少计算量的同时提供一个较高的频率分辨率，以便于更加精确的估计差频频率。另外，由于差频信号的相位与时间为线性一次关系，也可以使用最小二乘法对差频信号的相位进行一阶线性拟合，从而获得拟合直线的斜率，即差频频率，但是缺点是计算量大，且抗干扰性差。在涉及到非线性误差时，现有的方法估计的差频频率均会产生一定的估计误差。

技术实现思路

[0004]本专利技术的主要目的是提供一种FMCW雷达的差频频率估计方法，旨在解决现有技术中，在存在非线性误差的情况下，估计的差频频率发误差较大的技术问题
[0005]为实现上述目的，本专利技术提供一种FMCW雷达的差频频率估计方法，包括如下步骤：
[0006]获取待处理实信号，并将所述待处理实信号转换为解析信号；
[0007]获取所述解析信号的实际相位；
[0008]对所述实际相位在时域按照预设值进行平均分段，并对每一段进行最小二乘法相位拟合，获取每一段拟合直线的斜率，将所述拟合直线的斜率作为每一段的差频频率；
[0009]去除因为相位提取以及噪声干扰偏差较大的部分段的差频频率，对剩余部分段的差频频率取均值，并根据该均值计算出FMCW雷达的实际差频频率。
[0010]可选地，在所述获取所述解析信号的实际相位的步骤之前还包括：
[0011]获取解析信号s
b
(t)的包络幅值A(t)，从解析信号中消除信号的包络变化带来的影响，即s
b
(t)/A(t)。
[0012]可选地，所述对所述实际相位在时域按照预设值进行平均分段的步骤包括：
[0013]令预设分段数为k，整个回波数据的长度为N，对所述实际相位数据(t
i
，φ
i
)(i＝1，2，...，N)在时域进行平均分段，每段的数据长度为N
i
＝round(N/k)，k＞N/50，进而得到k组分段数据(N1，N2，...，N
k
)。
[0014]可选地，所述对所述实际相位在时域按照预设值进行平均分段，并对每一段进行最小二乘法相位拟合，获取每一段拟合直线的斜率，将所述拟合直线的斜率作为每一段的差频频率的步骤包括：
[0015]对所述k组分段数据(N1，N2，...，N
k
)做最小二乘法一阶线性拟合，得到k组拟合直线φ＝K
j
t+b
j
，(j＝1，2，...，k)，对每一段的斜率Kj求平均值将所述平均值作为所有数据的拟合直线的斜率；
[0016]其中，K
j
为每一段拟合直线的斜率，b
j
为每一段拟合直线的截距。
[0017]可选地，所述去除因为相位提取以及噪声干扰偏差较大的部分，对剩余部分取均值，并将该均值作为FMCW雷达的实际差频频率的步骤包括：
[0018]如果待拟合数据为理想线性调频信号，则将每一段的拟合直线斜率K
j
直接作为所述平均值
[0019]如果当调频信号中存在非线性误差时，选取更加稳定的部分的斜率取均值，将该均值作为所述平均值
[0020]可选地，FMCW雷达的所述实际差频频率满足如下关系式：
[0021][0022]本专利技术提供的基于分段最小二乘法的FMCW雷达差频频率估计方法，采用希尔伯特变换将待处理实信号转换为解析信号，利用相位解缠技术对解析信号的相位解缠绕，获得实际相位，对解缠绕的实际相位在时域按照预设值进行平均分段，对每一段进行最小二乘法相位拟合，获得每段拟合直线的直线斜率即每段的差频频率，去除因为相位提取以及噪声干扰偏差较大的部分，对剩余部分取均值，将该均值作为FMCW雷达的实际差频频率，仅需要对数据做分段最小二乘法拟合即可得到差频频率估计，简单精确。
附图说明
[0023]为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图示出的结构获得其他的附图。
[0024]图1为本专利技术提供的FMCW雷达的差频频率估计方法一实施例的流程示意图；
[0025]图2为图1中FMCW雷达的差频频率估计方法中数据的分段方法的结构示意图；
[0026]图3为图1中FMCW雷达的差频频率估计方法中使用非理想线性差频信号获取相位的分段拟合直线斜率变化情况示意图；
[0027]图4为图1中FMCW雷达的差频频率估计方法中获取的FMCW雷达的差频信号的实际相位示意图；
[0028]图5为图1中FMCW雷达的差频频率估计方法中获取的FMCW雷达的差频信号的实际相位分段拟合直线的斜率变化情况示意图。
[0029]本专利技术目的的实现、功能特点及优点将结合实施例，参照附图做进一步说明。
具体实施方式
[0030]下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术的一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。
[0031]需要说明，若本专利技术实施例中有涉及方向性指示，则该方向性指示仅用于解释在某一特定姿态下各部件之间的相对位置关系、运动情况等，如果该特定姿态发生改变时，则该方向性指示也相应地随之改变。
[0032]本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种FMCW雷达的差频频率估计方法，其特征在于，包括如下步骤：获取待处理实信号，并将所述待处理实信号转换为解析信号；获取所述解析信号的实际相位；对所述实际相位在时域按照预设值进行平均分段，并对每一段进行最小二乘法相位拟合，获取每一段拟合直线的斜率，将所述拟合直线的斜率作为每一段的差频频率；去除因为相位提取以及噪声干扰偏差较大的部分段的差频频率，对剩余部分段的差频频率取均值，并根据该均值计算出FMCW雷达的实际差频频率。2.根据权利要求1所述的FMCW雷达的差频频率估计方法，其特征在于，在所述获取所述解析信号的实际相位的步骤之前还包括：获取解析信号s
b
(t)的包络幅值A(t)，从解析信号中消除信号的包络变化带来的影响。3.根据权利要求1所述的FMCW雷达的差频频率估计方法，其特征在于，所述对所述实际相位在时域按照预设值进行平均分段的步骤包括：令预设分段数为k，整个回波数据的长度为N，对所述实际相位数据(t
i
，φ
i
)(i＝1，2，...，N)在时域进行平均分段，每段的数据长度为N
i
＝round(N/k)，k＞N/50，进而得到k组分段数据(N1，N2，...，N
k
)。4.根据权利要求3所述的...

【专利技术属性】
技术研发人员：赵博，张曰义，刘小军，方广有，
申请(专利权)人：中国科学院空天信息创新研究院，
类型：发明
国别省市：