一种基于局部极大值同步提取变换的微动特征与估计方法技术

本发明专利技术涉及一种基于局部极大值同步提取变换的微动特征提取与估计方法，属于雷达信号处理技术领域。本发明专利技术首先通过雷达采集目标回波数据，然后对回波数据进行预处理，再对距离维数据做傅里叶变换得到距离多普勒谱图，然后进行去抖动干扰，提取待分析的信号。对选取目标所在距离单元的回波信号进行短时傅里叶变换得到时频结果，基于局部极大值同步提取变换提取目标的微动特征，估计微动参数。本发明专利技术能够实现微多普勒信号高分辨特征提取，提高参数估计精度。估计精度。估计精度。

【技术实现步骤摘要】
一种基于局部极大值同步提取变换的微动特征与估计方法


[0001]本专利技术涉及一种基于局部极大值同步提取变换的微动特征提取与估计方法，可用于提取目标微动特征，属于雷达信号处理


技术介绍

[0002]近年来，由于无人机飞行器的便捷性和灵活性，其价格也相对合理，因此受到广泛应用和喜爱。无人机可以用于各种目的，例如侦查和监视、航拍、环境检测和灾害救援等。小型无人机带来了许多好处，但同时也存在潜在的危害。随着无人机数据量的增加，空域管理和监管变得更加复杂，可能会被不法分子利用，从而造成社会危害。此外，在机场环境中，小型旋翼无人机在未经授权的情况下飞行会影响民航飞机的起降，并存在与飞机相撞的潜在风险。无人机检测技术的发展一直以来都备受关注，它在各个领域中都发挥着重要作用，无人机、直升机等在雷达目标中是一类有旋转叶片的目标，此类目标的雷达回波信号中会包含旋转叶片的长度、转速等微动参数。分析目标的微多普勒特征可以提取这些参数，准确地提取这些微动参数是识别目标类型的重要依据。因此，对目标旋转叶片参数的提取分析是雷达信号领域的研究热点。
[0003]微动特性是旋翼无人机在飞行过程中微小振动和姿态变化的特性。对于旋翼无人机的微动特性深入研究，对无人机的分类识别具有重要意义。在微动特征分析方面，通常采用时频分析方法来分析目标的微动特性，短时傅里叶变换(Short Time Fourier Transform,STFT)能够分析信号在时间和频率上的变换，适用于分析短时平稳信号和非平稳信号，但是不能兼顾时间分辨率和频率分辨率。Wigner
‑
Ville分布(Wigner
‑
Ville Distribution,WVD)相较于STFT其时频分辨率得到了提升，但对于多分量信号其存在交叉项的影响。这些方法的时频分辨率较低，导致它们无法准确地刻画非平稳信号的非线性行为，近年来，一些学者对传统时频分析方法做了改进，如重新分配方法(RM)、同步压缩变换等。重新分配(RM)是一种时频分辨率较高的时频分析法方法，但它沿着时间和频率方向重新分配时频系数，导致无法重建信号。因此，选择合适的时频方法来提取目标微动特征，对目标的识别分类起着重要作用。

技术实现思路

[0004]本专利技术的目的在于克服传统时频分析方法提取旋翼微动特性时出现的严重时频能量扩散，无法精确估计微动参数问题，在同步提取变换的基础上，提出了局部极大值同步提取变换微动特征提取与估计方法，旨在改善微动时频特征分辨率，提高微动参数估计精度。
[0005]本专利技术的一种基于局部极大值同步提取变换的微动特征与估计方法，其特殊之处在于包括如下步骤：
[0006]步骤1：通过雷达采集目标数据；
[0007]所述步骤1的具体为:采集目标回波数据，分离目标所在距离单元和背景距离单元
信号时域序列，根据设定样本观测时长对雷达信号时间序列截取，得到信号样本序列；
[0008]步骤2：通过信号采集系统实时采集空域中的任意一距离段信号，并对信号进行预处理；
[0009]步骤3：对距离维数据做傅里叶变换得到距离多普勒谱图，提取待分析的信号，通过距离多普勒图像找到目标所在位置，选取目标所在距离单元的信号；
[0010]步骤4：对选取的信号进行去趋势项处理，去除信号中的抖动干扰；
[0011]所述步骤4)中去趋势项处理，去除信号中的抖动干扰，利用原始信号减去均值滤波后的低频部分以得到去除趋势项后的真实信号，其中均值滤波表示为：
[0012][0013]其中x(i)表示某点的信号，k表示某信号点处左右相邻信号点，2m为均值滤波滑窗所包含信号点数，通过计算某信号点处前m个点和后m个点的平均值来代替当前点，实现线性低通滤波。
[0014]步骤5：对数据进行短时傅里叶变换，得到时频结果G(t,ω)；
[0015]所述步骤5)对信号进行短时傅里叶变换表示为：
[0016][0017]其中，G表示短时傅里叶变换结果，t表示时间，ω表示频率，g(u
‑
t)表示高斯函数滑动窗口，s(u)表示回波信号，u表示回波信号的变量；
[0018]步骤6：根据步骤5)中的时频分析结果，加入额外的相位位移e
‑
iωt
，短时傅里叶变化结果可表示为然后对时间进行求导可求取瞬时频频率其中t表示时间，ω表示频率，g'表示高斯窗函数对时间求导。根据ω0(t,ω)得到同步提取变换结果Te(t,ω)＝G(t,ω)δ(ω
‑
ω0(t,ω))；
[0019]步骤7：对于多分量信号短时傅里叶变换后的数据，谱图集中在具有涂抹能量分布的瞬时频率轨迹上，其表示为其中k表示第几个信号分量，A
k
(t)表示第k个分量信号的瞬时幅值，表示相位，表示窗函数进行傅里叶变换后的结果，由此我们提出了一个新的频率提取算子ω
n
其中ω∈|ω
‑
Δ,ω+Δ|，用新的频率提取算子替换同步提取变换中的瞬时频率得到局部极大值同步提取变换LSET(t,ω)＝G(t,ω)δ(ω
‑
ω
n
)；
[0020]步骤8：使用局部极大值同步提取变换对短时傅里叶变换后的数据进行处理，通过时频谱图提取目标微动特征，根据时频谱图中的参数信息来计算微动参数。本专利技术首先通过雷达采集目标回波数据，然后对数据进行预处理，包括解调处理和去直流，再对距离维数据做傅里叶变换得到距离多普勒谱图，对其进行去抖动干扰，提取待分析的信号。将经过预处理和去抖动干扰后的数据进行短时傅里叶变换得到时频结果，之后基于局部极大值同步
提取变换提取目标微动特征，估计目标的微动参数。本专利技术能够实现微多普勒信号高分辨特征提取，提高参数估计精度。
附图说明
[0021]图1为本专利技术流程图；
[0022]图2为无人机信号去趋势前和去趋势后的信号；
[0023]图3为图2对信号进行短时傅里叶变换后的时频谱图；
[0024]图4为图3经本专利技术方法处理后的时频谱图；
[0025]图5为实施例中无人机旋翼模型。
具体实施方式
[0026]下面将结合附图对本专利技术的具体步骤进行详细说明。
[0027]实施例1
[0028]以调频连续波雷达观测旋翼无人机为例，本实施例基于局部极大值同步提取变换的微动特征和估计方法。参照图1，具体实施步骤如下。
[0029](1)步骤1，采集信号。
[0030](2)通过雷达和电脑连接，其中雷达参数设置调制带宽和调制周期为,脉冲周期，采集无人机雷达回波数据，分离目标所在距离单元和背景距离单元信号时域序列，根据设定样本观测时长对雷达信号时间序列截取，得到信号样本序列，对信号进行去趋势项处理，去除信号中的抖动干扰。该算法利用原始信号减去均值滤波后的低频部分以得到去除趋势项后的真实信号。其中均值滤波可以表示为：
[0031][0032]其中x(i)本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于局部极大值同步提取变换的微动特征与估计方法，其特征在于包括如下步骤：步骤1：通过雷达采集目标回波数据；步骤2：信号采集系统实时采集空域中的任意一距离段信号，并对信号进行预处理；步骤3：对距离维数据做傅里叶变换得到距离多普勒谱图，提取待分析的信号，通过距离多普勒图像找到目标所在位置；步骤4：对信号进行去趋势项处理，去除信号中的抖动干扰，选取目标所在距离单元的信号；步骤5：然后对信号进行短时傅里叶变换，得到时频结果G
e
(t,ω)；步骤6：根据步骤5)中的时频分析结果对时间求导可以得到瞬时频频率ω0(t,ω)，则同步提取算子为δ(ω
‑
ω0(t,ω))，进而得到同步提取变换结果Te(t,ω)；步骤7：对于多分量信号短时傅里叶变换后的数据，谱图集中在具有能量分布的瞬时频率轨迹上，其表示为|G(t,ω)|，由此我们定义了一个新的频率提取算子ω
n
(t,ω)，用新的频率提取算子替换同步提取变换中的瞬时频率得到局部极大值同步提取变换LSET(t,ω)；步骤8：使用局部极大值同步提取变换对短时傅里叶变换后的数据进行处理，通过时频谱图提取目标的微动特征，根据时频谱图中的参数信息来计算微动参数。2.按照权利要求1所述的一种基于局部极大值同步提取变换的微动特征与估计方法，其特征在于所述步骤1)具体为：，采集目标回波数据，分离目标所在距离单元和背景距离单元信号时域序列，根据设定样本观测时长对雷达信号时间序列截取，得到信号样本序列。3.按照权利要求1所述的一种基于局部极大值同步提取变换的微动特征与估计方法，其特征在于所述步骤4)中去趋势项处理，去除信号中的抖动干扰，利用原始信号减去均值滤波后的低频部分以得到去除趋势项后的真实信号，其中均值滤波表示为：其中x...

【专利技术属性】
技术研发人员：陈小龙，杜晓林，袁旺，于刚，汪兴海，刘克，关键，苏宁远，
申请(专利权)人：中国人民解放军海军航空大学，
类型：发明
国别省市：