一种基于welch的多频点多带宽识别方法技术

本发明专利技术提出的一种基于welch的多频点多带宽识别方法，属于通信技术领域，所述方法包括：采集基带信号进行welch功率谱估计，得到傅立叶变换的点数的功率谱Px；将功率谱Px进行m阶平滑滤波，再进行倒序进行m阶平滑滤波，得到功率谱Pz；将功率谱Pz进行分段，计算噪声值和全带宽阈值；计算出功率谱Pz的最大值max和对应的索引值、3db带宽阈值，判断功率谱Pz的最大值；得到全带宽右边索引值和全带宽左边索引值；求得3db带宽索引差；计算中心频点和3db带宽；进行次峰值信号的中心频点和3db带宽的计算；将接收机接收到的基带信号，往后移nfft，进行下一个welch功率谱估计。本发明专利技术实现了低信噪比情况下，信号中心频点和带宽的识别。信号中心频点和带宽的识别。信号中心频点和带宽的识别。

【技术实现步骤摘要】
一种基于welch的多频点多带宽识别方法


[0001]本专利技术涉及通信
，更具体的说是涉及一种基于welch的多频点多带宽识别方法。

技术介绍

[0002]中心频率估计和带宽估计主要涉及到雷达、声纳、导航、通信、成像、地质勘探、生物医学工程等诸多军事及民用领域。在这些领域中，接收机在接收带宽内同时接收多个不同的信号，而以上领域中需要对该类信号进行分析，进行特征信号的提取。而在进行特征信号的提取之前，首先要对该类信号进行中心频点和带宽的估计。
[0003]关于信号的频率估计，频域估计方法包括周期图法、质心法，其中周期图法是基于最大似然估计的方法，将周期图的最高峰值位置作为载波频率的估计，这种方法适用于存在载波分量的情况，载波被抑制时该方法不适用。质心法适用于频谱对称信号，对非对称信号估计效果不好。因此以上方式不适合接收带宽内的对多个信号进行频点带宽计算的场景。
[0004]常规的带宽估计方法包括均方根法、自相关法、最大熵法和能量集中法等。这几种方法估计出的带宽精度都相对较低，往往不能满足实际应用需求。
[0005]对于干扰比较大，低信噪比的情况，现有的中心频点带宽的估计方法也比较难识别出来。现有频点带宽的计算基于单一信号的居多，不满足实际情况的处理。而且对于mfsk（多进制数字频率调制）这种多峰值的情况，现有频点带宽的方式将它识别为多频点，与实际不符合。

技术实现思路

[0006]针对现有技术中存在的问题，本专利技术的目的在于提供一种基于welch的多频点多带宽识别方法，实现了低信噪比情况下，信号中心频点和带宽的识别。
[0007]本专利技术为实现上述目的，通过以下技术方案实现：一种基于welch的多频点多带宽识别方法，包括如下步骤：步骤1：取长度为N的基带信号进行welch功率谱估计，得到傅立叶变换的点数的功率谱Px；步骤2：将功率谱Px进行m阶平滑滤波，得到功率谱Py；步骤3：将功率谱Py倒序进行m阶平滑滤波，得到功率谱Pz；步骤4：将功率谱Pz进行分段，计算噪声值noise和全带宽阈值FullBwTH；步骤5：计算出功率谱Pz的最大值max和对应的索引值MaxInd、3db带宽阈值ThreeBwTH，其中ThreeBwTH=max/2，判断功率谱Pz的最大值max是否小于预设值，若是则转到步骤10，否则转到下一步；步骤6：根据全带宽阈值FULLBwTH，从索引值MaxInd向右搜索，得到全带宽右边索引值FullBwR，从索引值MaxInd向右搜索，得到全带宽左边索引值FullBwL；
步骤7：当FullBwR和FullBwL同时存在时，根据3db带宽阈值ThreeBwTH得到3db带宽的右边索引值ThreeBwR和3db带宽的左边索引值ThreeBwL，并求得3db带宽索引差；步骤8：根据预设公式计算中心频点和3db带宽；步骤9：将[FullBwL, FullBwR]区间内的Px置为噪声值noise，以进行次峰值信号的中心频点和3db带宽的计算；步骤10：将接收机接收到的基带信号，往后移nfft，并转到步骤1，以进行下一个频谱的计算。
[0008]进一步，步骤1具体为：将接收机接收到的基带信号，每次取长度的数据进行welch功率谱估计，得到点的功率谱；其中，是傅立叶变换的点数，为welch变换时的分段数,每段数据重合一半，每段的长度都是，段数为。
[0009]进一步，步骤2具体为：通过下式对功率谱进行阶的平滑滤波，得到功率谱；。
[0010]进一步，步骤3具体为：通过下式将功率谱倒序进行m阶平滑滤波，得到功率谱；。
[0011]进一步，步骤4具体为：将平滑后得到的功率谱分为段，每个功率谱求均值，取段中最小均值作为噪声值；通过公式计算全带宽阈值，其中，为频谱信号中心频点和带宽计算的阈值。
[0012]进一步，所述判断功率谱Pz的最大值max是否小于预设值具体为：判断功率谱Pz的最大值max是否小于。
[0013]进一步，步骤6包括：
从索引值MaxInd位置向右搜索，找到首次低于全带宽阈值FULLBwTH的索引值的前一个索引值记为全带宽右边索引值FullBwR；从索引值MaxInd位置向左搜索，找到首次低于全带宽阈值FULLBwTH的索引值的后一个索引值记为全带宽左边索引值FullBwL。
[0014]进一步，步骤7包括：如果同时存在FullBwR和FullBwL时计算3db带宽的索引值，从FullBwR到MaxInd往左搜索，找到首个大于ThreeBwTH的索引值，记为ThreeBwR；从FullBwR到MaxInd往右搜索，找到首个大于ThreeBwTH的索引值，记为ThreeBwL；如果FullBwL不存在，则令FullBwL=0，标记为flag=
‑
1；如果FullBwR不存在，则令FullBwR=nfft
‑
1，标记为flag=1；根据公式ThreeBandInd=ThreeBwR
‑
ThreeBwL，求得3db带宽索引差ThreeBandInd；根据公式deltaL= ThreeBwL
‑
FullBwL,求得全带宽3db带宽左边索引差deltaL；根据公式deltaR= FullBwR
ꢀ‑
ThreeBwR,求得全带宽3db带宽右边索引差deltaR；如果deltaL/deltaR＞2，FullBwL= ThreeBwL
‑
deltaR；如果deltaR /deltaL＞2，FullBwR= ThreeBwR+deltaL。
[0015]进一步，步骤8包括：判断flag是否为0，若不为0，则转到步骤5；若为0，则根据下列公式计算出中心频点FreCenter和3db带宽BandWidth；idth；其中，fs为采样率。
[0016]进一步，步骤9包括：将全带宽[FullBwL, FullBwR]索引值内的Px置为噪声值noise，跳转到步骤5，以进行次峰值信号的中心频点和带宽的计算。
[0017]对比现有技术，本专利技术有益效果在于：1、本专利技术算法简单，易于实现，实现了低信噪比情况下，信号中心频点和带宽的识别。
[0018]2、本专利技术可以正确地将mfsk识别出中心频点和它的有效带宽。
[0019]3、本专利技术可以识别多频点多带宽的信号的中心频点与带宽。
[0020]由此可见，本专利技术与现有技术相比，具有突出的实质性特点和显著的进步，其实施的有益效果也是显而易见的。
附图说明
[0021]为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图。
[0022]附图1是本专利技术具体实施方式的方法流程图。
[0023]附图2是本专利技术具体实施方式的fft频谱图。<本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于welch的多频点多带宽识别方法，其特征在于，包括如下步骤：步骤1：取长度为N的基带信号进行welch功率谱估计，得到傅立叶变换的点数的功率谱Px；步骤2：将功率谱Px进行m阶平滑滤波，得到功率谱Py；步骤3：将功率谱Py倒序进行m阶平滑滤波，得到功率谱Pz；步骤4：将功率谱Pz进行分段，计算噪声值noise和全带宽阈值FullBwTH；步骤5：计算出功率谱Pz的最大值max和对应的索引值MaxInd、3db带宽阈值ThreeBwTH，其中ThreeBwTH=max/2，判断功率谱Pz的最大值max是否小于预设值，若是则转到步骤10，否则转到下一步；步骤6：根据全带宽阈值FULLBwTH，从索引值MaxInd向右搜索，得到全带宽右边索引值FullBwR，从索引值MaxInd向右搜索，得到全带宽左边索引值FullBwL；步骤7：当FullBwR和FullBwL同时存在时，根据3db带宽阈值ThreeBwTH得到3db带宽的右边索引值ThreeBwR和3db带宽的左边索引值ThreeBwL，并求得3db带宽索引差；步骤8：根据预设公式计算中心频点和3db带宽；步骤9：将[FullBwL, FullBwR]区间内的Px置为噪声值noise，以进行次峰值信号的中心频点和3db带宽的计算；步骤10：将接收机接收到的基带信号，往后移nfft，并转到步骤1，以进行下一个频谱的计算。2.根据权利要求1所述的基于welch的多频点多带宽识别方法，其特征在于，所述步骤1具体为：将接收机接收到的基带信号，每次取长度的数据进行welch功率谱估计，得到点的功率谱；其中，是傅立叶变换的点数，为welch变换时的分段数,每段数据重合一半，每段的长度都是，段数为。3.根据权利要求2所述的基于welch的多频点多带宽识别方法，其特征在于，所述步骤2具体为：通过下式对功率谱进行阶的平滑滤波，得到功率谱；。4.根据权利要求3所述的基于welch的多频点多带宽识别方法，其特征在于，所述步骤3具体为：
通过下式将功率谱倒序进行m阶平滑滤波，得到功率谱；。5.根据权利要求4所述的基于welch的多频点多带宽识别方法，其特征在于，所述步骤4具体为：将平滑后得到的功率谱分为段，每个功率谱求均值，取段中最小均值作为噪声值；通过公式计算全带宽阈值，其中，为频谱信号中心频点和带宽计算的阈值。6.根据权利要求5所述的基于welch的多频点多带宽识别方法，其特征在于，所述判断功率谱Pz的最...

【专利技术属性】
技术研发人员：郑敏江，谢艳，
申请(专利权)人：中孚信息股份有限公司，
类型：发明
国别省市：