CN107783129B - 一种旋翼无人机防撞毫米波雷达信号处理方法 - Google Patents

Abstract

旋翼无人机防撞毫米波雷达信号处理方法，属于信号处理领域，用于解决旋翼无人机低空飞行时易发生与障碍物之间的碰撞，导致旋翼无人机的损坏的问题，技术要点是：S1.对通道1和通道2中的A/D采集到的IQ数据去直流；S2.对去直流后的通道1和通道2中A/D采集到的IQ数据，进行FFT变换，将时域数据转换成频率数据；S3.将FFT变换后的复数模值做CFAR门限检测，输出过门限的第一个峰值点，得到通道1中上扫频频率值和下扫频对应的频率值，以及通道2中的上扫频频率值，并计算通道1和通道2中，根据各自上扫频分别计算得到相位。

Description

一种旋翼无人机防撞毫米波雷达信号处理方法

技术领域

本发明属于信号处理领域，涉及一种雷达信号处理方法。

背景技术

近几年，随着技术的不断发展，民用小型旋翼无人机价格越来越低，被广泛用于航拍、电影拍摄、农药喷洒、现场救援、大地遥感测绘、高压线电网巡视等领域。但是因为旋翼无人机低空飞行时易发生与障碍物之间的碰撞，导致旋翼无人机的损坏。目前威胁旋翼无人机室外低空飞行安全的物体主要有树木等自然物体以及电力线实施例、电线杆、建筑物等人造物体。

由于毫米波雷达工作波长介于1mm～10mm之间，与其他的探测方式相比，主要有探测性能稳定、环境适应良好、尺寸小、价格低，可以在相对恶劣的雨雪天气使用等优点。因此，本发明重点介绍基于毫米波雷达的无人机避障功能方法的实现。

发明内容

为了解决旋翼无人机低空飞行时易发生与障碍物之间的碰撞，导致旋翼无人机的损坏的问题，本发明提出了一种旋翼无人机防撞毫米波雷达信号处理方法，以解算得到障碍物的速度、距离和角度，从而可以进行障碍物规避。

为了实现上述目的，本发明的技术方案是：

一种旋翼无人机防撞毫米波雷达信号处理方法，包括如下步骤：

S1.对通道1和通道2中的A/D采集到的IQ数据去直流；

S2.对去直流后的通道1和通道2中A/D采集到的IQ数据，进行FFT变换，将时域数据转换成频率数据；

S3.将FFT变换后的复数模值做CFAR门限检测，输出过门限的第一个峰值点，得到通道1中上扫频频率值和下扫频对应的频率值，以及通道2中的上扫频频率值，并计算通道1和通道2中，根据各自上扫频分别计算得到相位；

S4.使用步骤S3中得到的通道1中上扫频频率值和下扫频对应的频率值，计算得到无人机前向障碍物目标的距离；

S5.将步骤S3中得到的通道1和通道2中，根据各自上扫频分别计算得到的相位计算方位角。

进一步的，所述步骤S3中，

设通道1中上扫频段第一个过门限点的峰值坐标为p1_up，则该点对应的频率值为f1_up，对应的FFT变换后的数据为a_p1_up+1j*b_p1_up，相位为

通道2中上扫频段第一个过门限点的峰值坐标为p2_up，则该点对应的频率值为f2_up，对应的FFT变换后的数据为a_p2_up+1j*b_p2_up，相位为

设通道1中下扫频段第一个过门限点的峰值坐标为p1_down，则该点对应的频率值为f1_down；

其中：a表示I路的数据值，b表示Q路的数据值，a_p1表示在a+j*b组成的数组中，过门限的峰值点对应的坐标为p1，a_p2表示在a+j*b组成的数组中，过门限的峰值点对应的坐标为p2，b_p1表示在a+j*b组成的数组中，过门限的峰值点对应的坐标为p1，b_p2表示在a+j*b组成的数组中，过门限的峰值点对应的坐标为p2。

进一步的，所述步骤S4，将步骤S3中得到的通道1中上扫频频率值f1_up和下扫频对应的频率值f1_down，根据公式

计算得到无人机前向障碍物目标的距离，其中T为三角波周期，B为调频带宽，c为光速；

根据公式

计算得到无人机前向障碍物目标的速度，其中f₀为中心频率。

进一步的，所述步骤S5，将步骤3中得到的通道1和通道2中，根据各自上扫频分别计算得到的相位

和

根据计算公式

得到相位差为Δψ；根据公式

计算方位角，其中d为天线间距，λ为雷达波长。

进一步的，所述步骤S1去直流的方法是：计算通道1中AD采集到的上、下扫频IQ两路数据的均值，IQ两路各数据点分别减掉计算得到的均值；计算通道2中AD采集到的上扫频IQ两路数据的均值，IQ两路各数据点分别减去计算得到的均值。

进一步的，所述步骤S2中，还包括加窗的步骤，位于去直流的步骤之后。

进一步的，所述方法还包括步骤S6.滤波跟踪，并预测下一测量时刻的距离以及速度值。

进一步的，所述滤波，使用α-β滤波器，其常增益滤波器的预测方程为

X(k+1/k)＝ΦX(k/k)；

其滤波方程为：

X(k+1/k+1)＝X(k+1/k)+K[Z(k+1)-H(k+1/k)]；

其中，X(k/k)为k时刻的滤波值，X(k+1/k)为k时刻对下一个时刻的预测值，Z(k)是k时刻的观测值；

当目标运动方程采用常速度模型时，常数增益矩阵K＝[α,β/T]^T,其状态转移矩阵

此模型的量测矩阵为H＝[1，0]；

其中：0＜α＜1，0＜β＜1。

有益效果：

1、本发明首次给出了一种基于线性调频三角波实现旋翼无人机防撞毫米波雷达系统的波形设计；

2、本发明给出基于线性调频三角波实现的旋翼无人机防撞毫米波雷达信号处理分系统的处理过程，该分系统可以实现对前方障碍物的相对距离以及相对速度的检测，同时可以实现目标方向角的检测功能。

附图说明

图1线性调频三角波FMCW在一个扫频周期内的频率变化图；

图2实施例旋翼无人机短距离防撞系统信号处理流程图。

具体实施方式

实施例1：一种旋翼无人机防撞毫米波雷达信号处理方法，包括如下步骤：

S1.对通道1和通道2中的A/D采集到的IQ数据去直流；所述步骤S1去直流的方法是：计算通道1中AD采集到的上、下扫频IQ两路数据的均值，IQ两路各数据点分别减掉计算得到的均值；计算通道2中AD采集到的上扫频IQ两路数据的均值，IQ两路各数据点分别减去计算得到的均值；该步骤主要是起到去直流的目的。

S2.对去直流后的通道1和通道2中A/D采集到的IQ数据，进行FFT变换，将时域数据转换成频率数据；所述步骤S2中，还包括加窗的步骤，位于去直流的步骤之后。去直流后的I、Q数据，合并成I+jQ的数据形式，然后进行加窗处理，将对通道1中的上扫频段和下扫频段、通道2中的上扫频段，各自数据进行加窗处理。可以选择汉宁窗或是海明窗等，降低旁瓣，从而提高目标的检测性能；汉宁窗会导致主瓣加宽并降低，但是旁瓣会显著减小。

汉宁窗计算公式为

S3.将FFT变换后的复数模值做CFAR门限检测，输出过门限的第一个峰值点，得到通道1中上扫频频率值和下扫频对应的频率值，以及通道2中的上扫频频率值，并计算通道1和通道2中，根据各自上扫频分别计算得到相位；作为一种实施例，所述步骤S3中，主要是考虑到对无人机飞机危险程度最大的为距离无人机最近的物体，所以并不是找所有过门限的最大值，而是选择第一个过门限的峰值。设通道1中上扫频段第一个过门限点的峰值坐标为p1_up，则该点对应的频率值为f1_up，对应的FFT变换后的数据为a_p1_up+1j*b_p1_up，相位为

通道2中上扫频段第一个过门限点的峰值坐标为p2_up，则该点对应的频率值为f2_up，对应的FFT变换后的数据为a_p2_up+1j*b_p2_up，相位为

设通道1中下扫频段第一个过门限点的峰值坐标为p1_down，则该点对应的频率值为f1_down；其中：a表示I路的数据值，b表示Q路的数据值，a_p1表示在a+j*b组成的数组中，过门限的峰值点对应的坐标为p1，a_p2表示在a+j*b组成的数组中，过门限的峰值点对应的坐标为p2，b_p1表示在a+j*b组成的数组中，过门限的峰值点对应的坐标为p1，b_p2表示在a+j*b组成的数组中，过门限的峰值点对应的坐标为p2。

S4.使用步骤S3中得到的通道1中上扫频频率值和下扫频对应的频率值，计算得到无人机前向障碍物目标的距离；

作为一种实施例：所述步骤S4，将步骤S3中得到的通道1中上扫频频率值f1_up和下扫频对应的频率值f1_down，根据公式

计算得到无人机前向障碍物目标的距离，其中，T为三角波周期，T＝20ms，B为调频带宽，B＝200MHz，c为光速，c＝3.0×10⁸；

根据公式

计算得到无人机前向障碍物目标的速度，其中f₀为中心频率，f₀＝24.125GHz；

S5.将步骤S3中得到的通道1和通道2中，根据各自上扫频分别计算得到的相位计算方位角。

作为一种实施例：所述步骤S5，将步骤3中得到的通道1和通道2中，根据各自上扫频分别计算得到的相位

和

根据计算公式

得到相位差为Δψ；根据公式

计算方位角，其中，d为天线间距。

作为一种实施例，还包括步骤：S6.滤波跟踪，并预测下一测量时刻的距离以及速度值。

进一步的，由上述步骤则完成基于锯齿波和恒频波的组和波形实现的无人机短距离防撞毫米波雷达系统的信号处理，对于单一目标的相对速度、相对距离以及相应的方位角的解算过程。

当无人机短距离防撞毫米波雷达系统完成对目标的相对速度、相对距离以及相应的方位角的解算过程后，需要进行滤波跟踪模块。由于系统输出数据刷新率快，在很短时间内距离以及速度等变化量很小，可以近似看作为匀速运动，可以通过一定的算法来估计高度的变化率，预测下一测量时刻的距离以及速度值等。跟踪与预测方法是自适应跟踪跟踪滤波的前提和基础。主要方法目前有，线性自回归滤波，维纳滤波，加权最小二乘滤波，α-β与α-β-γ滤波，卡尔曼滤波和简化的卡尔曼滤波等。

本发明推荐采用α-β滤波器。由于α-β滤波器适用于跟踪误差变化速率相对较均匀的情况，所以基本适合无人机飞行场景。

α-β滤波器中，其常增益滤波器的预测方程为X(k+1/k)＝ΦX(k/k)，其滤波方程为X(k+1/k+1)＝X(k+1/k)+K[Z(k+1)-H(k+1/k)]，其中，X(k/k)为k时刻的滤波值，X(k+1/k)为k时刻对下一个时刻的预测值，Z(k)是k时刻的观测值。

当目标运动方程采用常速度模型时，常数增益矩阵K＝[α,β/T]^T,其状态转移矩阵

此模型的量测矩阵为H＝[1，0]。α-β滤波器就是满足长增益矩阵K，状态转移矩阵Φ以及量测矩阵H分别为上述表达式描述的常增益滤波器，即

α-β滤波器中的参数α和β的选择对于跟踪的响应、收敛速度以及跟踪稳定性都有关系。一般要求0＜α＜1，0＜β＜1。工程中，可以根据公式计算得到α和β的取值，即

和

其中k为次数，α和β随着k的变化取不同的值，实际在使用中，这两个参数会趋于定值。

经过α-β滤波器对单次结算的目标速度、距离就可以实现滤波跟踪以及预测。可以更好地实现目标的跟踪，同时使得输出的数据更加平滑，减少异常值的出现，有效的提高系统的稳定性。

实施例2：作为实施例1的补充，本实施例主要是完成旋翼无人机对其飞行前方环境障碍物的距离、速度以及方位进行测量。前方障碍物主要针对人，树、墙、网以及高压线等目标。

本实施例所设计的毫米波雷达的工作频率在24GHz或77GHz，采用FMCW连续波体制，采用线性调频其距离分辨率高。波形采用线性调频三角波FMCW，主要是因为本实施例要实现对目标距离以及速度的计算。通过三角波的上扫频和下扫频可以实现目标距离以及速度解算。本实施例设计的旋翼无人机的最大飞行速度为40km/h，无人机防撞的最大测距为50m，比目前市面上的无人机防撞距离高出3倍以上。

本实施例主要是给出无人机防撞毫米波雷达信号处理部分的设计以及信号处理方法。

本实施例设计的雷达中心频率f为24.125GHz。发射波形选择三角波，周期为20ms，带宽为200MHz。发射波形如图1所示。

本实施例通过单路IQ数据实现对目标距离速度的解算，由于本实施例实现目标方位角的计算，所以本实施例采用双接收天线的方式，即双通道IQ数据，通过对双通道各自上扫频段的计算实现该目标的测角功能。

旋翼无人机防撞毫米波雷达信号处理流程图，如图2所示，具体实现步骤如下：

1.计算通道1中AD采集到的上、下扫频IQ两路数据的均值，IQ两路各数据点分别减掉计算得到的均值；计算通道2中AD采集到的上扫频IQ两路数据的均值，IQ两路各数据点分别减去计算得到的均值。该步骤主要是起到去直流的目的。

2.对去直流后的通到1和通道2中A/D采集到的IQ数据，进行FFT变换，将时域数据转换成频率数据。

3.实施例将FFT变换后的复数模值做CFAR门限检测，输出过门限的第一个峰值点，主要是考虑到对无人机飞机危险程度最大的为距离无人机最近的物体，所以并不是找所有过门限的最大值，而是选择第一个过门限的峰值。

设通道1中上扫频段第一个过门限点的峰值坐标为p1_up，则该点对应的频率值为f1_up，对应的FFT后的数据为a_p1_up+1j*b_p1_up，相位

通道2中上扫频段第一个过门限点的峰值坐标为p2_up，则该点对应的频率值为f2_up，对应的FFT后的数据为a_p2_up+1j*b_p2_up，相位

设通道1中下扫频段第一个过门限点的峰值坐标为p1_down，则该点对应的频率值为f1_down。

4.将步骤三中得到的通道一中上扫频频率值f1_up和下扫频对应的频率值f1_down，根据公式

其中，T为三角波周期，T＝20ms，B为调频带宽，B＝200MHz，c为光速，c＝3.0×10⁸；根据公式

其中f₀为中心频率，f₀＝24.125GHz。根据这两个公式，得到无人机前向障碍物目标的距离和速度。

5.将步骤3中得到的通道1和通道2中，根据各自上扫频分别计算得到的相位

和

计算根据计算公式

得到相位差为Δψ。

根据公式

计算方位角，其中，d为天线间距。

由上述步骤则完成旋翼无人机防撞毫米波雷达对无人机运行前方障碍物距离、速度以及方位角等信息的的解算功能。

实施例3：对于上述各方案中，峰值处理，本实施例提供一种应用于无人机信号的峰值处理方法：

设置一个峰值点阈值因子α，其用于限制检测出的过门限最大峰值点与上一周期出现的最大峰值点的差值绝对值，使得该差值绝对值不得大于该峰值点阈值因子α：

表达式如下：

|L_max(k)-L_max(k-1)|≤α；

其中：L_max(k)为k周期的过门限最大峰值点坐标，L_max(k-1)为上一周期的最大峰值点坐标，k表示第k时刻；v_max为无人机最大飞行速度，λ为毫米波雷达波长，fs为采样率，N为FFT的点数；

如果k时刻，过门限最大峰值点与k-1时刻过门限最大峰值点的绝对值差值在所设置的峰值点阈值因子α范围内，则认为第k周期的峰值点有效；如果k时刻，过门限最大峰值点超过所设置的峰值点阈值因子α，则k时刻输出的峰值点用k-1时刻的峰值点进行替换。

作为上述技术手段的解释，在相邻周期的一个时间单元内，当前周期解算出的峰值点，与上个周期的峰值点，如果在相邻周期内，速度没有发生变化，则峰值点在相邻周期内也会保持不变，但是如果在相邻周期时间内，无人机水平飞行速度发生变化，会导致当前周期的峰值点在上一周期的峰值点发生一定的变化，如果是无人机靠近目标，则当前周期的点数会小于上一周期的点数，如果无人机远离目标，则当前周期的点数会大于上一周期的点数，该峰值点的变化范围即是所设计的峰值点阈值因子α，该因子选取的取值范围，主要取决于在相邻周期内，无人机的最大飞行速度，即公式

其中v_max为无人机最大飞行速度，λ为毫米波雷达波长，fs为采样率，N为FFT的点数。

但是如果旋翼无人机飞行环境发生突变后，对应的过门限的峰值点数也可能会连续发生超出所设计的阈值因子。如果不进行修正，发生突变后，每个周期检测到的过门限最大峰值点都会超过设置的阈值因子，每次过门限最大峰值点坐标都会被修正为上一时刻的峰值点坐标，即同理值也会保持突变前的值，不能适应突变后的值。为了提高无人机对各种环境的适应能力，为此引入一个峰值点突变累计因子φ。

设置一个峰值点突变累计因子φ，该峰值点突变累计因子φ的定义为，如果从k时刻开始，连续b个周期，b的取值范围为5～10，过门限最大峰值点与前一周期的过门限最大峰值点相比，都超过阈值门限因子a，则第k+b时刻，将当前时刻解算出的过门限最大峰值点作为当前时刻的过门限最大峰值点。为了保证跟踪的实时性，建议b的取值为5～10个。

通过上一步得出过门限最大峰值点后，为了提高表系统值测量的精度，提出提高测距精度的谱最大估计算法。

理想情况下，回波差频信号的频谱只有一个谱线，但是实际在使用过程中，由于采样存在栅栏效应，离散频谱最大幅值谱线必然会发生偏移谱峰位置，从而通过峰值点计算出的距离值与实际距离将会存在一定的误差。当谱峰发生偏移的时候，相对于主瓣峰值所对应的中央谱线将会两种情况，即左偏或是右偏。如果过门限最大值峰值点的左右峰值中，左边峰值大于右边峰值，则中央谱线所在的位置，在最大峰值点与左边峰值点之间，反之，则在最大峰值点与右边峰值点之间。

由于FFT计算得到的频谱对连续距离普等间距采样，其频谱幅值最大点必定位于其曲线的主瓣内，主瓣内有且仅有两个采样点。设过门限最大峰值点A1的坐标为(a1,k1)，其中，a1表示过门限最大峰值点的值，k1表示过门限峰值点对应的幅度值；最大峰值点左右两边，次峰值点坐标为A3(a3，k3)，设所求的中央峰值点A为(amax，kmax)，则e＝amax-a1，则A1点，关于A点对称点A2坐标为(a2,k1)＝(a1+2e,k1)，复包络的零点A4为(a4,k1)＝(a3+e,0)；

其中：a2、a3、a4是对应点的过门限最大峰值点的值，k3、k4是对应点的过门限峰值点对应的幅度值；

A2、A3和A4近似为一条直线，其线性关系为：

令

则

设定误差E与偏差e进行比对，如果|e|<E，则此时的过门限峰值点的值则为所要求的中央峰值点的值，如果偏差e大于所设定的误差E时，

β为修正因子，取值范围为1.5～1.9，该修正因子的选取理由是：由于初始的时候

A点对称点A2坐标为(a2,k1)＝(a1+2e,k1)，初始条件时A点横轴坐标点与A2横轴坐标是关于最大峰值点对称的，即A2的坐标点是a1+2e，如果偏差e大于所设定的误差E时，说明A2的坐标选取过大，也即是最大峰值点在a1+2e之间，2倍的偏差e需要进行取小，本发明采用的修正方法是，通过改变修正因子β的大小从而改变l值，然后进行e的不断迭代，直到e小于设定的误差E为止。修正因子β的取值原则可以根据所需求达到的E值进行选取，如果E需求精度不高，修正因子β可以选择1.9进行修正，如果E需求精度很高，可能需要多次迭代达到要求，则需要修正因子β尽量选择小一点，可以选择1.5进行修正，本发明给出了一个快速解算出最大峰值点的修正因子的区间范围值，即修正因子β＝1.5～1.9。改变修正因子计算出e的值，以计算得到中央峰值点的值amax＝a1+e。

作为另一种实施例，还包括步骤：距离跟踪：设置一个阈值因子ε，其用于限制当前距离数据H(k)与上一周期出现的距离数据H(k-1)的差值绝对值，使得该差值绝对值不得大于该阈值因子ε；

表达式如下：

|H(k)-H(k-1)|≤ε，ε取值范围为0.8～1.3；

如果k时刻的数据与k-1时刻的绝对值差值，在所设置的阈值因子ε范围内，则认为第k周期的峰值点有效；如果k时刻，数据超过所设置的阈值因子ε，则k时刻输出的数据用k-1时刻的数据进行替换。

设置一个突变累计因子θ，该突变累计因子θ的定义为，如果从k时刻开始，连续b个周期，数据与前一周期的数据相比，都超过阈值门限因子θ，则第k+b时刻，将当前时刻解算出的数据作为当前时刻的数据。

作为一种实施例，具体到本实施例中，对于上述未执行距离跟踪或执行了距离跟踪的，输出时，对于单次输出的距离数据，采用滑窗算法进行距离值的输出；

第k时刻的数据等于滑窗中的N_c个值去掉最大值和最小值后的均值，作为最后的数据输出，其计算公式为

其中N_c表示滑窗所采用的数据点数。

采用峰值跟踪算法和跟踪算法，可以有效避免由于单次或是多次峰值搜索的错误而导致一次或是多次数据解算的异常现象，如在单次峰值搜索过程中，发生峰值跳变，相邻周期之间的峰值差值很大，同时由与峰值的跳变，而引起的发生很大的跳变，即该周期内，峰值跳变引起的的跳变范围，已经远远大于由无人机速度引起的一个周期所产生的距离变化范围。由此峰值跟踪以及跟踪可以有效避免这种异常峰值导致的异常值，从而有效地的提高跟踪的数据的稳定度。

以上所述，仅为本发明创造较佳的具体实施方式，但本发明创造的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明创造披露的技术范围内，根据本发明创造的技术方案及其发明构思加以等同替换或改变，都应涵盖在本发明创造的保护范围之内。

Claims (8)

1.一种旋翼无人机防撞毫米波雷达信号处理方法，其特征在于，包括如下步骤：

S1.对通道1和通道2中的A/D采集到的IQ数据去直流；

S2.对去直流后的通道1和通道2中A/D采集到的IQ数据，进行FFT变换，将时域数据转换成频率数据；

S3.将FFT变换后的复数模值做CFAR门限检测，输出过门限的第一个峰值点，得到通道1中上扫频频率值和下扫频对应的频率值，以及通道2中的上扫频频率值，并计算通道1和通道2中，根据各自上扫频分别计算得到相位；

CFAR门限检测的过门限峰值点的处理的方法，其特征在于：设置一个峰值点阈值因子α，其用于限制检测出的过门限最大峰值点与上一周期出现的最大峰值点的差值绝对值，使得该差值绝对值不得大于该峰值点阈值因子α：

表达式如下：

|L_max(k)-L_max(k-1)|≤α；

其中：L_max(k)为k周期的过门限最大峰值点坐标，L_max(k-1)为上一周期的最大峰值点坐标，k表示第k时刻；v_max为无人机最大飞行速度，λ为毫米波雷达波长，fs为采样率，N为FFT变换的点数，FFT变换的对象是加窗后的锯齿波数据；

如果k时刻，过门限最大峰值点与k-1时刻过门限最大峰值点的绝对值差值在所设置的峰值点阈值因子α范围内，则认为第k周期的峰值点有效；如果k时刻，过门限最大峰值点超过所设置的峰值点阈值因子α，则k时刻输出的峰值点用k-1时刻的峰值点进行替换；

S4.使用步骤S3中得到的通道1中上扫频频率值和下扫频对应的频率值，计算得到无人机前向障碍物目标的距离；

S5.将步骤S3中得到的通道1和通道2中，根据各自上扫频分别计算得到的相位计算方位角。

2.如权利要求1所述的一种旋翼无人机防撞毫米波雷达信号处理方法，其特征在于，所述步骤S3中

设通道1中上扫频段第一个过门限点的峰值坐标为p1_up，则该点对应的频率值为f1_up，对应的FFT变换后的数据为a_p1_up+1j*b_p1_up，相位为

通道2中上扫频段第一个过门限点的峰值坐标为p2_up，则该点对应的频率值为f2_up，对应的FFT变换后的数据为a_p2_up+1j*b_p2_up，相位为

设通道1中下扫频段第一个过门限点的峰值坐标为p1_down，则该点对应的频率值为f1_down；

其中：a表示I路的数据值，b表示Q路的数据值，a_p1表示在a+j*b组成的数组中，过门限的峰值点对应的坐标为p1，a_p2表示在a+j*b组成的数组中，过门限的峰值点对应的坐标为p2，b_p1表示在a+j*b组成的数组中，过门限的峰值点对应的坐标为p1，b_p2表示在a+j*b组成的数组中，过门限的峰值点对应的坐标为p2。

3.如权利要求1所述的一种旋翼无人机防撞毫米波雷达信号处理方法，其特征在于，所述步骤S4，将步骤S3中得到的通道1中上扫频频率值f1_up和下扫频对应的频率值f1_down，根据公式

计算得到无人机前向障碍物目标的距离，其中T为三角波周期，B为调频带宽，c为光速；

根据公式

计算得到无人机前向障碍物目标的速度，其中f₀为中心频率。

4.如权利要求1所述的一种旋翼无人机防撞毫米波雷达信号处理方法，其特征在于，所述步骤S5，将步骤3中得到的通道1和通道2中，根据各自上扫频分别计算得到的相位

和

根据计算公式

得到相位差为Δψ；根据公式

计算方位角，其中d为天线间距，λ为雷达波长。

5.如权利要求1所述的旋翼无人机防撞毫米波雷达信号处理方法，其特征在于，所述步骤S1去直流的方法是：计算通道1中AD采集到的上、下扫频IQ两路数据的均值，IQ两路各数据点分别减掉计算得到的均值；计算通道2中AD采集到的上扫频IQ两路数据的均值，IQ两路各数据点分别减去计算得到的均值。

6.如权利要求1所述的旋翼无人机防撞毫米波雷达信号处理方法，其特征在于，所述步骤S2中，还包括加窗的步骤，位于去直流的步骤之后。

7.如权利要求1所述的旋翼无人机防撞毫米波雷达信号处理方法，其特征在于，还包括步骤S6.滤波跟踪，并预测下一测量时刻的距离以及速度值。

8.如权利要求7所述的旋翼无人机防撞毫米波雷达信号处理方法，其特征在于，所述滤波，使用α-β滤波器，其常增益滤波器的预测方程为

X(k+1/k)＝ΦX(k/k)；

其滤波方程为：

X(k+1/k+1)＝X(k+1/k)+K[Z(k+1)-H(k+1/k)]；

其中，X(k/k)为k时刻的滤波值，X(k+1/k)为k时刻对下一个时刻的预测值，Z(k)是k时刻的观测值；

当目标运动方程采用常速度模型时，常数增益矩阵K＝[α,β/T]^T,其状态转移矩阵

此模型的量测矩阵为H＝[1，0]；

其中：0＜α＜1，0＜β＜1。

Images