[发明专利]一种民航飞机散射信号多普勒频率估计方法

说明书

技术领域

本发明属于无线电信号参数估计技术领域，涉及到多普勒频率估计的方法，特别涉及到民航飞机散射信号中的多普勒频率估计方法。

背景技术

地面信源信号经飞机散射到接收机，由于飞机相对地面信源和接收机有相对位置移动，因此会产生多普勒效应。在散射信号的分段处理过程中，现有方法假设分段长度内的多普勒频率是不变的，不适用于对多普勒频率具有高精度要求的情况。本发明利用线性调频信号对短时间内的多普勒频率进行近似，并提出一种改进的分数阶傅里叶变换对其实现估计。分数阶傅里叶变换对线性调频信号具有最佳的能量聚集特性，被广泛应用于该类信号的参数估计当中，主要包括模值平方最大法、分数阶相关法、分数阶傅里叶高阶矩方法等，但以上算法均对数据进行直接处理，未考虑数据截断带来的分辨率降低以及频谱泄漏等问题。因此，本发明提出一种高斯加权分数阶傅里叶变换，并将之应用于飞机散射信号的多普勒频率估计当中。

发明内容

本发明采用线性调频信号对飞机散射信号中的多普勒频率进行近似，针对数据截断的影响，提出了一种飞机散射信号处理中高精度多普勒频率的估计方法。

一种民航飞机散射信号多普勒频率估计方法，主要包括以下步骤：

A.将采集得到的数据分解成固定长度为T的M段数据；

B.对每段数据采用高斯加权分数阶傅里叶变换进行多普勒频率估计；

C.对多普勒频率进行后处理。

所述步骤B具体包含以下步骤：

B1.采用高斯加权函数进行加权，形式为设置参数c选择加权函数形状，由经验设置参数c的范围为2≤c≤4；

B2.粗遍历，设置分数阶傅里叶变换参数区间p∈[0,2]，遍历步长为0.01；

B3.粗估计，搜索B2结果的模平方最大值对应分数阶傅里叶域位置(p₁，u₁)，设定B4中精细遍历区间[p₂，p₃]，其中p₂＝max(0，p₁-0.1)，p₃＝min(p₁+0.1，2)；

B4.精细遍历，由B3中设定的遍历区间，设置遍历步长为0.001；

B5.参数估计，搜索B4结果的模平方最大值对应分数阶傅里叶域位置(p₄，u₄)，根据下式实现分段时间内多普勒频率的中心频率和多普勒频率变化率的估计m＝1,2,…,M：

所述步骤C具体包含以下步骤：

C1.依据民航飞机的飞行速度等先验信息设置多普勒频率变化率的阈值Th1；

C2.根据多普勒频率变化率对估计结果进行初步修正，具体包含以下三个步骤；

C21.搜索得到B5中小于C1步骤中阈值的第一个时间位置，开始遍历；

C22.如果则

C23.如果则

C3.对C2中结果依据多普勒频率估计进行再次修正，具体包含以下步骤：

C31.对多普勒频率估计进行一阶差分

C32.设定阈值Th3，Th4，搜索C31中差分结果中大于阈值Th3的峰值，如果相邻峰值异号且峰值间距离小于阈值Th4，进行C33中修正，否则保留。

C33.设定C32中两峰值位置为t_i和t_i+1，则求得多普勒频率均值为利用与对多普勒频率及其变化率进行修正，其中m＝t_i+1,…,t_i+1。

C4.根据采样率以及C3中结果得到每个时间点的瞬时多普勒频率。

本发明利用高斯加权分数阶傅里叶变换进行多普勒频率以及频率变化率的估计，最后根据民航飞机速度等先验知识对频率变化率及其对应的多普勒频率实现修正，并对修正后的结果进行中值滤波以得到更加精确的多普勒频率估计。该方法步骤少，估计准确，能够满足实际需求。

附图说明

为了更清楚的说明本发明的实施例或现有技术的技术方案，下面将对实施例或现有技术描述中所需要使用的附图做一简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1是本发明的算法流程图。

图2A是本发明中飞机散射信号的实部时域波形图。

图2B是本发明中飞机散射信号的虚部时域波形图。