一种民航飞机散射信号多普勒频率估计方法技术

本发明专利技术属于无线电信号处理技术领域，提供一种飞机散射信号多普勒频率提取方法。其特征是首先将采集信号进行分段处理，其次采用两级高斯加权分数阶傅里叶变换进行多普勒中心频率及多普勒频率变化率的估计；然后根据估计得到的多普勒频率变化率进行初步修正，再根据估计得到的多普勒中心频率进行二次修正；最后根据采样频率等信息实现瞬时多普勒频率的估计。实验证明本发明专利技术算法性能良好，在低信噪比情况下能够获得多普勒频率的有效估计。

【技术实现步骤摘要】
一种民航飞机散射信号多普勒频率估计方法
本专利技术属于无线电信号参数估计
，涉及到多普勒频率估计的方法，特别涉及到民航飞机散射信号中的多普勒频率估计方法。
技术介绍
地面信源信号经飞机散射到接收机，由于飞机相对地面信源和接收机有相对位置移动，因此会产生多普勒效应。在散射信号的分段处理过程中，现有方法假设分段长度内的多普勒频率是不变的，不适用于对多普勒频率具有高精度要求的情况。本专利技术利用线性调频信号对短时间内的多普勒频率进行近似，并提出一种改进的分数阶傅里叶变换对其实现估计。分数阶傅里叶变换对线性调频信号具有最佳的能量聚集特性，被广泛应用于该类信号的参数估计当中，主要包括模值平方最大法、分数阶相关法、分数阶傅里叶高阶矩方法等，但以上算法均对数据进行直接处理，未考虑数据截断带来的分辨率降低以及频谱泄漏等问题。因此，本专利技术提出一种高斯加权分数阶傅里叶变换，并将之应用于飞机散射信号的多普勒频率估计当中。
技术实现思路
本专利技术采用线性调频信号对飞机散射信号中的多普勒频率进行近似，针对数据截断的影响，提出了一种飞机散射信号处理中高精度多普勒频率的估计方法。一种民航飞机散射信号多普勒频率估计方法，主要包括以下步骤：A.将采集得到的数据分解成固定长度为T的M段数据；B.对每段数据采用高斯加权分数阶傅里叶变换进行多普勒频率估计；C.对多普勒频率进行后处理。所述步骤B具体包含以下步骤：B1.采用高斯加权函数进行加权，形式为设置参数c选择加权函数形状，由经验设置参数c的范围为2≤c≤4；B2.粗遍历，设置分数阶傅里叶变换参数区间p∈[0,2]，遍历步长为0.01；B3.粗估计，搜索B2结果的模平方最大值对应分数阶傅里叶域位置(p1，u1)，设定B4中精细遍历区间[p2，p3]，其中p2＝max(0，p1-0.1)，p3＝min(p1+0.1，2)；B4.精细遍历，由B3中设定的遍历区间，设置遍历步长为0.001；B5.参数估计，搜索B4结果的模平方最大值对应分数阶傅里叶域位置(p4，u4)，根据下式实现分段时间内多普勒频率的中心频率和多普勒频率变化率的估计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中两峰值位置为ti和ti+1，则求得多普勒频率均值为利用与对多普勒频率及其变化率进行修正，其中m＝ti+1,…,ti+1。C4.根据采样率以及C3中结果得到每个时间点的瞬时多普勒频率。本专利技术利用高斯加权分数阶傅里叶变换进行多普勒频率以及频率变化率的估计，最后根据民航飞机速度等先验知识对频率变化率及其对应的多普勒频率实现修正，并对修正后的结果进行中值滤波以得到更加精确的多普勒频率估计。该方法步骤少，估计准确，能够满足实际需求。附图说明为了更清楚的说明本专利技术的实施例或现有技术的技术方案，下面将对实施例或现有技术描述中所需要使用的附图做一简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1是本专利技术的算法流程图。图2A是本专利技术中飞机散射信号的实部时域波形图。图2B是本专利技术中飞机散射信号的虚部时域波形图。图3是本专利技术飞机散射信号的时频分析图。图4A是本专利技术初步估计得到多普勒频率图。图4B是本专利技术初步估计得到多普勒频率变化率图。图5是本专利技术多普勒频率及变化率修正流程图。图6A是本专利技术修正后得到的多普勒频率图。图6B是本专利技术修正后得到的多普勒频率变化率图。图7是本专利技术多普勒频率及变化率二次修正流程图。图8A是本专利技术二次修正后得到的多普勒频率图。图8B是本专利技术二次修正后得到的多普勒频率变化率图。图9是本专利技术得到的多普勒频率估计图。具体实施方式为使本专利技术实施例的目的、技术方案及其优点更加清楚，下面结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚完整的描述，整体算法流程图如图1所示：A.将采集得到的数据分解成固定长度为T的M段数据；采集得到的数据是IQ复信号，由感兴趣的多普勒频率及背景噪声组成，波形如图2A和图2B所示，其时频图如图3所示。可以看出，接收数据信噪较低，感兴趣信号完全被噪声所淹没，为更高精度地估计多普勒频率，利用线性调频信号对多普勒频率进行近似，并采用对线性调频信号具有最佳能量聚集性能的分数阶傅里叶变换进行处理。将接收数据分成固定长度为2049的多段数据，为克服分段处理过程中数据截断的影响，提出利用高斯加权分数阶傅里叶变换进行处理。所述步骤B具体包含以下步骤：B1.采用高斯加权函数进行加权，形式为设置参数c＝2；B2.粗遍历，根据民航飞机速度范围、接收机位置等先验知识，设定分数阶傅里叶变换参数区间p∈[0.9,1.1]，遍历步长为0.01；B3.粗估计，搜索B2结果的模平方最大值对应分数阶傅里叶域位置(p1，u1)，设定B4中精细遍历区间[p2,p3]，其中p2＝max(0，p1-0.1)，α3＝min(p1+0.1,2)；B4.精细遍历，由B3中设定的遍历区间，设置遍历步长为0.001；B5.参数估计，搜索B4结果的模平方最大值对应分数阶傅里叶域位置(p4，u4)，根据下式实现该分段时间内多普勒频率的中心频率和多普勒频率变化率的估计m＝1,2,…,M：初步估计得到的结果如图4A和4B所示。所述步骤C具体包含以下步骤：C1.依据民航飞机的飞行速度等先验信息设置多普勒频率变化率的阈值Th1＝5，即认为分段长度为2049时，该段时间内的多普勒频率变化率不大于5Hz/s；C2.根据多普勒频率变换率对多普勒频率和频率变化率进行修正，具体包含以下三个步骤，如图5所示。C21.搜索得到B5中小于C1阈值的第一个时间位置；C22.如果则C23.设定Th2＝3即两段数据间的多普勒频率变化率之差的绝对值不超过3Hz/s，如果则根据多普勒频率变化率的初步修正结果如图6A和图6B所示。C3.对C2中结果依据多普勒频率估计进行再次修正，具体包含以下步骤，如图7所示：C31.对多普勒频率估计进行一阶差分C32.设定阈值Th3＝10即相邻数据段多普勒频率的中心频率之差不超过10Hz，如果大于阈值则该段数据估计值可能为野点，搜索C31中差分结果中大于阈值Th3的峰值，设定Th4＝5若相邻峰值异号且峰值间距离小于阈值Th4，进行C33中修正，否则保留。C33.设定C32中两峰值位置为ti和ti+1，则求得多普勒频率均值为利用与对多普勒频率及其变化率进行修正，其中m＝ti+1，…，ti+1。根据多普勒频率修正得到的估计结果如图8A和图8B所示。C4.根据采样率以及C3中结果得到每个时间点的瞬时多普勒频率，如图9所示。本文档来自技高网...

【技术保护点】
一种飞机散射信号处理中高精度多普勒频率的估计方法，其特征在于具有如下步骤：A.将采集得到的数据分解成固定长度为T的M段数据；B.对每段数据采用高斯加权分数阶傅里叶变换进行多普勒频率估计；具体步骤如下：B1.采用高斯加权函数进行加权，形式为设置参数c选择加权函数形状，设置参数c的范围为2≤c≤4；B2.粗遍历，设置分数阶傅里叶变换参数区间p∈[0,2]，遍历步长为0.01；B3.粗估计，搜索B2结果的模平方最大值对应分数阶傅里叶域位置(p1,u1)，设定B4中精细遍历区间[p2,p3]，其中p2＝max(0,p1‑0.1)，α3＝min(p1+0.1,2)；B4.精细遍历，由B3中设定的遍历区间，设置遍历步长为0.001；B5.参数估计，搜索B4结果的模平方最大值对应分数阶傅里叶域位置(p4,u4)，根据下式实现该分段时间内多普勒频率的中心频率和多普勒频率变化率的估计m＝1,2,…,M：μ^m=-12πcot(12p4π)f^m=12πu4csc(12p4π)]]>C.对多普勒频率进行后处理；具体步骤如下：C1.依据民航飞机的先验信息设置多普勒频率变化率的阈值Th1；C2.根据多普勒频率变化率对估计结果进行初步修正；包含以下步骤：C21.搜索得到B5中小于C1阈值的第一个时间位置，开始遍历；C22.如果μ^m=0,]]>则μ^m=μ^m-1,f^m-f^m-1+μ^m-1·T;]]>C23.如果|μ^m-μ^m-1|>Th2,]]>则μ^m=μ^m-1,μ^m-f^m-1+μ^m-1·T;]]>C3.对C2中结果依据多普勒频率估计进行再次修正；包含以下步骤：C31.对多普勒频率估计进行一阶差分C32.设定阈值Th3，Th4，搜索C31中差分结果中大于阈值Th3的峰值，如果相邻峰值异号且峰值间距离小于阈值Th4，进行C33中修正，否则保留；C33.设定C32中两峰值位置为ti和ti+1，则求得多普勒频率均值为mf=(f^ti+1+1-f^ti)/(ti+1-ti+2),]]>利用f^m=f^ti+mf×(m-ti)]]>与μ^m=mf/T]]>对多普勒频率及其变化率进行修正，其中m＝ti+1,…,ti+1；C4.根据采样率以及C3中结果得到每个时间点的瞬时多普勒频率。...

【技术特征摘要】
1.一种飞机散射信号处理中高精度多普勒频率的估计方法，其特征在于具有如下步骤：A.将采集得到的数据分解成固定长度为T的M段数据；B.对每段数据采用高斯加权分数阶傅里叶变换进行多普勒频率估计；具体步骤如下：B1.采用高斯加权函数进行加权，形式为设置参数c选择加权函数形状，设置参数c的范围为2≤c≤4；B2.粗遍历，设置分数阶傅里叶变换参数区间p∈[0,2]，遍历步长为0.01；B3.粗估计，搜索B2结果的模平方最大值对应分数阶傅里叶域位置(p1,u1)，设定B4中精细遍历区间[p2,p3]，其中p2＝max(0,p1-0.1)，p3＝min(p1+0.1,2)；B4.精细遍历，由B3中设定的遍历区间，设置遍历步长为0.001；B5.参数估计，搜索B4结果的模平方最大值对应分数阶傅里叶域位置(p4,u4)，根据下式实现分段时间内多普勒频率的中心频率和多普勒频率变化率的估计

【专利技术属性】
技术研发人员：邱天爽，王鹏，栾声扬，史益新，马济通，
申请(专利权)人：大连理工大学，
类型：发明
国别省市：辽宁;21