一种基于散射变换的无人机目标旋翼参数估计方法

摘要

本发明公开了一种基于线性调频连续波雷达的无人机旋翼参数估计方法。步骤1：初始化散射变换尺度，生成小波滤波器组；步骤2：对回波信号进行一阶散射变换；步骤3：对散射系数设置能量阈值，并进行第二阶散射变换，估计频谱带宽；步骤4：提取微多普勒分量所在的第二阶散射系数，提取周期分量进行FFT得到频谱，通过频率分量数目对旋翼数目进行判断；步骤5：若无人机旋翼数目为0，则不存在旋翼，跳转至8，否则进入6；步骤6：根据旋翼数目N和峰值间隔频率对转动频率进行估计；步骤7：结合微动频谱带宽和转动频率反推叶片尺寸；步骤8：输出旋翼参数估计结果。本发明能实现线性调频连续波雷达对无人机目标旋翼转速、数目和尺寸的估计。

说明书

技术领域

本发明属于的技术领域；具体涉及一种基于散射变换的无人机目标旋翼参数估计方法。

背景技术

随着商用无人机使用规模的增长，无人机滥用带来的安全问题正越来越受到关注。目前有很多手段对无人机目标进行检测，相比于其他方法，雷达系统具有全天候、全天时，探测距离远的优势，线性调频连续波雷达凭借结构简单、分辨率高、无距离盲区等优点取代脉冲体制雷达，可以作为对中近程低空无人机探测的重要手段。但是由于无人机目标一般尺寸较小，运动速度较慢，飞行高度较低，容易与飞鸟或其他杂波混淆，给检测和识别造成困难。

由于无人机一般具有一个或多个旋翼，这些转动的旋翼会在雷达回波中产生微多普勒效应。微多普勒效应由目标自身结构的微运动产生，可以有效反映目标的运动状态、尺寸和形状等信息，且这种特征具有“唯一性”，对于无人机目标而言，通过微多普勒频率对旋翼的尺寸、转速和数目等特征参数进行估计，对无人机目标进行分类和威胁度评估具有重要作用。

但是通常情况下，小型无人机的旋翼长度相对较小，一般在亚米级，导致微多普勒分量能量较低，而转速相对较快，一般为每分钟几千转，会形成较大的多普勒频率展宽，给参数提取造成困难。且旋翼叶片间的相互作用会在回波中产生转动频率的高次谐波分量，对转动频率的估计带来干扰。

现有的方法线性时频分析为主，如短时傅里叶变换，其时间分辨率和频率分辨率相互制约导致，限制了能量聚集效果，降低了转动频率的估计效果，某些双线性时频分析如Cohen类方法虽然能量积累效果有所提高，但是对多分量信号存在交叉项问题，且现有方法会受到主体回波的影响，很难对旋翼参数进行估计。

由于上述原因，传统非平稳信号分析方法并不适用于对无人机旋翼参数的估计。

发明内容

本发明提供了一种基于散射变换的无人机目标旋翼参数估计方法，能实现线性调频连续波雷达对无人机目标旋翼转速、数目和尺寸的估计。

本发明通过以下技术方案实现：

一种基于线性调频连续波雷达的无人机旋翼参数估计方法，所述无人机旋翼参数估计方法包括以下步骤：

步骤1：对线性调频连续波雷达的回波信号进行散射变换尺度的初始化，进而生成小波滤波器组；

步骤2：对回波信号进行一阶散射变换；

步骤3：对无人机旋翼分量的散射系数设置能量阈值，选择目标所在的散射系数进行第二阶散射变换，通过对应的协变分量根据Carson准则估计频谱带宽f

步骤4：对无人机旋翼分量的散射系数设置能量阈值，提取微多普勒分量所在的第二阶散射系数，提取周期分量进行FFT得到频谱，通过频率分量数目对无人机旋翼数目N进行判断；

步骤5：若无人机旋翼数目N为0，则不存在旋翼，跳转至步骤8输出结果，否则进入步骤6；

步骤6：根据旋翼数目N和峰值间隔频率Δf对转动频率进行估计

步骤7：结合微动频谱带宽f

步骤8：输出旋翼参数估计结果。

进一步的，所述步骤2具体为，对于时间信号x，其一阶散射变换为，

其中，

进一步的，所述步骤3具体为散射系数S

随着散射变换阶数的增多，小波尺度函数时宽减小，频率分辨率更高；

对于信号，一阶散射系数提取基频成分，而二阶散射系数提取信号中的调制特征。

进一步的，所述步骤4具体为，线性调频连续波雷达的无人机目标回波差频信号为，

其中A

进一步的，所述无人机目标回波由一个固定频率信号和一组调频信号组成。

进一步的，中频由瞬时距离R

本发明的有益效果是：

本发明有效的对小型无人机的旋翼参数进行估计，在信噪比高于-16dB时，对于旋翼参数估计的可靠性在90％以上，对旋翼转动频率和叶片长度的估计误差在10％以内。

附图说明

附图1本发明算法流程图。

附图2一阶散射系数结果示意图。

附图3二阶散射系数示意图。

附图4微多普勒分量频谱展宽示意图。

附图5周期分量频谱示意图。

附图6不同信噪比条件下的参数估计准确率示意图。

具体实施方式

下面将结合本发明实施例中的附图对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

实施例1

目前需要对旋翼无人机的参数估计方法进行研究。小波变换通过改变小波尺度函数实现多分辨分析，得到更高的时频分辨率，在此基础上散射变换通过对小波模算子的迭代平均保持信号的局部形变稳定性，保留了局部细节信息，防止信号能量发散，也有助于提取时变的频率调制信息。因此可以通过散射变换实现无人机旋翼参数的估计。

一种基于线性调频连续波雷达的无人机旋翼参数估计方法，所述无人机旋翼参数估计方法包括以下步骤：

步骤1：对线性调频连续波雷达的回波信号进行散射变换尺度的初始化，进而生成小波滤波器组；

步骤2：对回波信号进行一阶散射变换；

步骤3：对无人机旋翼分量的散射系数设置能量阈值，选择目标所在的散射系数进行第二阶散射变换，通过对应的协变分量根据Carson准则估计频谱带宽f

步骤4：对无人机旋翼分量的散射系数设置能量阈值，提取微多普勒分量所在的第二阶散射系数，提取周期分量进行FFT得到频谱，通过频率分量数目对无人机旋翼数目N进行判断；

步骤5：若无人机旋翼数目N为0，则不存在旋翼，跳转至步骤8输出结果，否则进入步骤6；

步骤6：考虑到每个旋翼一般具有两个叶片，根据旋翼数目N和峰值间隔频率Δf对转动频率进行估计

步骤7：结合微动频谱带宽f

步骤8：输出旋翼参数估计结果。

进一步的，所述步骤2具体为，对于时间信号x，其一阶散射变换为，

其中，

进一步的，所述步骤3具体为散射系数S

随着散射变换阶数的增多，小波尺度函数时宽减小，频率分辨率更高；

对于信号，一阶散射系数提取基频成分，而二阶散射系数提取信号中的调制特征。

进一步的，所述步骤4具体为，线性调频连续波雷达的无人机目标回波差频信号为，

其中A

进一步的，所述无人机目标回波由一个固定频率信号(主体回波)和一组调频信号(旋翼回波)组成。

进一步的，中频由瞬时距离R

第一阶散射系数中以目标主体位置产生的频率为中心频率，对主体运动状态进行提取，选出相应的第一阶散射系数，对其进行第二阶散射变换，此时回波旋翼产生的微多普勒频率体现为二阶散射系数强度的变化，通过相应的二阶散射系数直接对转动频率进行提取，据此来对转动频率进行估计；

一阶散射系数对应的协变部分U

实施例2

本发明实例中，设置线性调频连续波雷达信号参数如下：

无人机目标旋翼转动频率为30Hz，叶片长度为12cm，每个旋翼有两个叶片。设置目标距离雷达天线300m，运动速度为1m/s。

选取128个调频周期进行积累，则信号长度为0.128s，由于积累时间较短，在此过程中，旋翼的转速可以近似看作不会发生突变，当无人机状态稳定时，每个旋翼的转速可以视为大致相同。

S1、对回波进行第一阶散射变换，结果如图2所示，通过阈值检测可以得到目标信号分量处于第五个尺度的散射系数中；

S2、对其进行第二阶散射系数，其结果如图3所示，通过阈值检测可以得到旋翼微多普勒信号分量位于第6个尺度的二阶散射系数中；

S3、在完成转动速度估计后，通过微动分量分离结果对微多普勒带宽进行估计如图5所示，根据Carson准则，可以得到单边微多普勒展宽频率约为f

S4、提取该散射系数的周期变化结果，对其进行FFT得到频谱如图4所示，由于只有一个频率峰值，可以判断目标的旋翼数目为1；

S5、判断峰值数目大于0，存在旋翼，继续执行；

S6、且可以得到变化频率为62.5Hz，由于旋翼上有两个叶片，因此可以估计出旋翼转动频率为f

S7、此时雷达信号波长为

S8、输出旋翼参数估计结果。

通过该方法可以对旋翼参数进行有效估计，结果与设置参数相符。

通过蒙特卡洛实验对本文所提方法在不同信噪比和不同旋翼数目条件下的估计性能，结果如图6所示，设置信噪比以5dB为步长从-20dB变化到20dB，旋翼数目分别设置为1，2和4，进行400次蒙特卡洛实验，通过图6可以看出，该方法在信噪比高于-16dB的情况下可以得到较好的估计效果。