本发明涉及雷达信号处理，具体涉及一种基于相位估计的fmcw雷达高精度测距方法。
背景技术：
1、调频连续波(fmcw)雷达在工业液位测量中有着非常明显的优势。其非接触性测量特性使其能够避免与液体直接接触，适用于各种液体类型，包括腐蚀性和高粘度液体；其连续测量能力和高分辨率确保了高精度的液位数据，而且能够应对复杂液体环境；fmcw雷达的远程测量能力适用于大型容器或深井等场景，而其优秀的抗干扰能力使其在工业环境中表现出色。实时监测能力和扫频信号的准确性使其成为工业液位监测的重要选择，为生产控制和优化提供了可靠的数据支持。
2、目前的fmcw雷达一般利用发射的扫频信号和回波信号之间频率差进行测距，但是测距精度不高。目前大多数处理差拍信号的方法都是利用fft得到频谱信号，接着对信号进行归一化处理，得到信号的频率谱。通过找到频率谱的最高点来估计差拍信号的频率。在实际应用中，当信号频率不是整数倍的采样频率，频谱泄露会导致主要峰值模糊，从而降低了距离分辨率和测距精度。此外，在信号频率恰好处于fft频率网格的间隙时，会出现栅栏效应，相邻的谱线之间的真实谱线可能丢失，这会造成一定的频率测量误差，进而导致测距误差。与频率测量的方法相比较，相位法测量距离的精度更高，一般高两个数量级，而且稳定性更好，但是相位法存在相位模糊的问题。
技术实现思路
1、为了克服传统方法中的频谱泄露和栅栏效应以及相位模糊的问题。本发明提出了一种基于相位估计的fmcw雷达高精度测距方法，该方法在利用高速ad采集回波信号后与发送的扫频信号进行混频得到差拍信号后，首先利用fft获得频率信息，然后通过相位法来提取相关的相位信息；通过结合这两种信息，可以减小误差，提高系统的测距精度，更准确地确定目标距离，适用于工业高精度测量。
2、本发明采取的技术方案为：
3、基于相位估计的fmcw雷达高精度测距方法，包括以下步骤：
4、s1：获取扫频信号经过目标反射后的第一回波信号，经过目标反射后的第一回波信号和扫频信号混频得到差拍信号sif(t)；
5、s2：基于s1中得到的差拍信号sif(t)，获取目标反射相位分量φr；
6、s3：基于s1中得到的差拍信号sif(t)，在进行ad采集后得到离散信号sif(n)；
7、s4：基于s3得到的离散信号sif(n)，去除波导色散对离散信号sif(n)的相位影响，得到离散差拍信号sif+(n)；
8、s5：基于s4得到的离散差拍信号sif+(n)，对其进行fft变换得到频谱信号，接着对频谱信号进行归一化处理，得到频谱图，找到幅度谱峰值处的频率f0，在f0两侧的预设频率区间内，对频谱进行czt细化，得到细化频率f1；
9、s6：在s5得到的细化频率f1附近寻找两个幅度最接近的频率点，利用这三个点进行多项式拟合，并利用拟合曲线估计第一估计频率
10、s7：估计频谱的最高点处的频率，分析在第一估计频率附近的相位谱，在其附近的相位谱中进行线性插值，接着对估计频谱最高点的相位进行展开并进行插值处理，求得第一估计相位
11、s8：基于s7中得到的第一估计频率和fmcw雷达的参数，根据距离估计公式，求出第一估计距离
12、s9：基于s8中求得的第一估计距离，根据公式求得时间延迟根据求出第二估计相位
13、s10：基于s7中得到的第一估计相位和s9中得到的第二估计相位求得相位间隔nui,根据相位间隔nui，求出第三估计相位
14、s11：根据相位估计和频率估计距离公式，求出第二估计距离
15、所述s1中，经过目标反射后的第一回波信号和扫频信号混频后的差拍信号表达式为：
16、
17、其中：aif为差拍信号的幅值，为发送的扫频信号的带宽与信号频率上升时间之比，τ为时间延迟，ω0为扫频信号的中心频率，φr为目标反射相位分量。
18、所述s2中，目标反射相位分量φr表示发射的扫频信号遇到物体之后反射造成的相位变化。所述s3中，离散信号sif(n)＝sif(t)|t＝nt，其中，t为采样间隔。
19、所述s4中，在电磁波传播过程中，由于波导色散的影响，产生相位速度，造成相位误差，会造成一定程度的传播延迟，利用矩形波导来模拟波导色散造成的影响。
20、雷达系统射频路径相位响应的主要影响是天线馈电波导中的传播延迟色散；在长度为lwg矩形波导中，由于相速度vp(t)引起的额外传播延迟公式如下：
21、
22、其中，ωwg是波导截止频率；
23、波导色散干扰的解析差拍信号为：
24、
25、频谱为：
26、s′if(jω)＝sif(jω)*wrf(jω)
27、为了消去色散的影响，利用如下公式：
28、
29、得到的sif+(n)为消除色散之后的差拍信号。
30、所述s4中，离散差拍信号
31、所述s5中，预设频率区间为f0±f，在预设频率区间内进行频率细化，其中，f为频率细化带宽的一半，细化分辨率为δf，带宽长度为f2-f1，f1为czt细化频率带宽的起始频率，f2为czt细化频率带宽的终点频率，δf＝(f2-f1)/m，m为细化频段的频点数。所述s6中，离散的频率点经过多项式拟合之后近似抛物线，抛物线的表达式为y＝ax2+bx+c；将三个点带入抛物线表达式中求得拟合曲线，对抛物线进行求导y′＝2ax0+b＝0，通过求导公式求得x0的值即为第一估计频率
32、基于频率估计的第一估计距离其中，cm为电磁波传播的速度，为第一估计频率。
33、所述s7中，相位谱可能存在折叠的现象，在相位区间内存在相位跳跃；
34、
35、其中：为展开的信号，为要展开的信号，如果要展开的信号和前一个信号之间的差值大于π，那么前一个信号加上2π；如果要展开的信号与后一个信号的差值大于π，那么后一个信号减去2π；是频率估计的索引值；是向下取整数，假如为5.2，那么为5；是向上取整数，假如为6.2，那么为为6。
36、所述s8中，基于频率估计求得的第二估计距离其中，cm为电磁波传播的速度。
37、所述s9中，第二估计相位其中，ω0是载波信号的角频率；为时间延迟；φr是经过物体得到的相位分量。
38、所述s10中，相位间隔其中，为目标相位的值，
39、为求出第二估计相位在整数采样点上寻找一个值，使得估计相位值与目标相位值的差最小；
40、第三估计相位具体表达式：
41、
42、其中：为第一估计相位；nui为相位间隔。
43、所述s11中，相位估计距离公式为：频率估计公式：
44、
45、其中：为第一估计频率；b为带宽，tr为t为扫频信号扫频时间，cm为光速。
46、本发明一种基于相位估计的fmcw雷达高精度测距方法，技术效果如下：
47、1)本发明所提供的fmcw雷达高精度测距方法，在进行傅里叶变换之后，进行czt细分，接着进行抛物线拟合，利用拟合曲线估计最高点的频率。相比于直接利用czt细分后得到的最高点频率值，该方法减少了栅栏效应的影响，具有更高的精度。
48、2)本发明所提供的fmcw雷达高精度测距方法，在处理差拍信号的时候考虑到了波导色散的影响，在进行频率估计和相位估计之前去除波导色散的影响，进一步地提高了测量的精度；
49、3)本发明提出了一种频率法和相位法结合的方法，在进行频率法测量之后，利用频率法估计的距离求出时间延迟和第二估计相位利用频率法得到的参数进行相位估计，最终求得测量的距离。与单一的频率法测量距离相比较，频率法和相位法相结合的方法具有更高的测量精度。
1.基于相位估计的fmcw雷达高精度测距方法，其特征在于包括以下步骤：
2.根据权利要求1所述基于相位估计的fmcw雷达高精度测距方法，其特征在于：所述s1中，经过目标反射后的第一回波信号和扫频信号混频后的差拍信号表达式为：
3.根据权利要求1所述基于相位估计的fmcw雷达高精度测距方法，其特征在于：所述s4中，在电磁波传播过程中，利用矩形波导来模拟波导色散造成的影响，具体如下：
4.根据权利要求1所述基于相位估计的fmcw雷达高精度测距方法，其特征在于：所述s5中，预设频率区间为f0±f，在预设频率区间内进行频率细化，其中，f为频率细化带宽的一半，细化分辨率为δf，带宽长度为f2-f1，f1为czt细化频率带宽的起始频率，f2为czt细化频率带宽的终点频率，δf＝(f2-f1)/m，m为细化频段的频点数。
5.根据权利要求1所述基于相位估计的fmcw雷达高精度测距方法，其特征在于：
6.根据权利要求1所述基于相位估计的fmcw雷达高精度测距方法，其特征在于：所述s7中，相位谱可能存在折叠的现象，在相位区间内存在相位跳跃；
7.根据权利要求1所述基于相位估计的fmcw雷达高精度测距方法，其特征在于：所述s8中，基于频率估计求得的第二估计距离其中，cm为电磁波传播的速度。
8.根据权利要求1所述基于相位估计的fmcw雷达高精度测距方法，其特征在于：所述s9中，第二估计相位
9.根据权利要求1所述基于相位估计的fmcw雷达高精度测距方法，其特征在于：所述s10中，相位间隔其中，为目标相位的值，为求出第二估计相位在整数采样点上寻找一个值，使得估计相位值与目标相位值的差最小；
10.根据权利要求1所述基于相位估计的fmcw雷达高精度测距方法，其特征在于：
