一种高精度频谱校正方法技术

本申请公开了一种高精度频谱校正方法，基于最大旁瓣衰减窗函数的临近频谱间的比例特性，将加权信号频谱实部和虚部拆分，仍能够分别进行三点插值计算。在得到实部和虚部两个归一化频率估计后再采用“数学平均”方式，这里处理方式经验证能够完全消除共轭负频率的泄露干扰，使得算法能够在任意采样窗口长度下实现高精度频偏估计。本发明专利技术提供的算法能够完全消除共轭负频率频谱泄露干扰，使得改进后的三点插值算法在任意采样窗口下实现高精度频谱校正。

A high precision spectrum correction method

【技术实现步骤摘要】
一种高精度频谱校正方法
本申请涉及实时频率估计
，尤其涉及一种高精度频谱校正方法。
技术介绍
在许多工程应用中，如通信、音响系统、雷达、声纳、电力系统、测量和仪表等，经常需要对正弦信号进行实时频率估计。特别对于通信领域，易于实现调制和解调的正弦波多用作其他基带信号的载波，其参数的准确估计和测量是实现准确调制载波通信的前提和关键。目前，基于离散傅里叶变换(DiscreteFourierTransform,DFT)的正弦信号分析算法，因易于实现且具有高效的快速计算方式(FastFourierTransform,FFT)，一直以来是通信工程信号处理应用中优先选择的分析方法。然而，DFT算法对于绝对低频信号及相对低频信号的估计，尤其是极短时数据记录情况下，近似解析插值关系式中所忽略的负频率泄漏干扰量级较大，从而使得基于插值DFT的频率校正估计严重失准，影响极短时数据记录情况下正弦信号频率参数校正估计的精度。
技术实现思路
本申请提供了一种高精度频谱校正方法，以解决极短时数据记录情况下正弦信号频率参数校正估计的精度低的技术问题。为了解决上述技术问题，本申请实施例公开了如下技术方案：本申请实施例公开了一种高精度频谱校正方法，包括：采集短时离散采样信号；对所述短时离散采样信号加窗并进行DFT计算，得到变换后频谱函数；通过局部峰值搜索所述变换后频谱函数，得到峰值；选取峰值所在谱线及相邻的两条谱线作为目标谱线，对所述目标谱线对应的实部和虚部进行拆分，生成实部拆分结果和虚部拆分结果；对所述实部拆分结果进行三点插值计算，得到实部计算结果；对所述虚部拆分结果进行三点插值计算，得到虚部计算结果；对所述实部计算结果和所述虚部计算结果进行数学平均计算，得到高精度频谱校正结果。可选的，对所述短时离散采样信号加窗并进行DFT计算，得到变换后频谱函数，包括：利用窗函数对所述短时离散采样信号进行加窗处理，得到加窗后频谱函数；对所述加窗后频谱函数进行DFT计算，得到变换后频谱函数。可选的，所述窗函数为H项余弦窗函数，表示为：式中，N为采样点数，ah为窗函数系数。可选的，对所述实部计算结果和所述虚部计算结果进行数学平均计算，得到高精度频谱校正结果，包括：对所述实部计算结果进行归一化处理，得到实部归一化结果；对所述虚部计算结果进行归一化处理，得到实部归一化结果；对所述实部归一化结果和所述实部归一化结果进行数学平均计算，得到高精度频谱校正结果。与现有技术相比，本申请的有益效果为：本申请提供了一种高精度频谱校正方法，基于最大旁瓣衰减窗函数的临近频谱间的比例特性，将加权信号频谱实部和虚部拆分，仍能够分别进行三点插值计算。在得到实部和虚部两个归一化频率估计后再采用“数学平均”方式，这里处理方式经验证能够完全消除共轭负频率的泄露干扰，使得算法能够在任意采样窗口长度下实现高精度频偏估计。本专利技术提供的算法能够完全消除共轭负频率频谱泄露干扰，使得改进后的三点插值算法在任意采样窗口下实现高精度频谱校正。应当理解的是，以上的一般描述和后文的细节描述仅是示例性和解释性的，并不能限制本申请。附图说明为了更清楚地说明本申请的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，对于本领域普通技术人员而言，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。图1为本实施例提供的一种高精度频谱校正方法流程图。具体实施方式为了使本
的人员更好地理解本申请中的技术方案，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本申请保护的范围。结合图1所示，本申请提供了高精度频谱校正方法，包括：采集短时离散采样信号。对所述短时离散采样信号加窗并进行DFT计算，得到变换后频谱函数。通过局部峰值搜索所述变换后频谱函数，得到峰值。选取峰值所在谱线及相邻的两条谱线作为目标谱线，对所述目标谱线对应的实部和虚部进行拆分，生成实部拆分结果和虚部拆分结果。对所述实部拆分结果进行三点插值计算，得到实部计算结果。对所述虚部拆分结果进行三点插值计算，得到虚部计算结果。对所述实部计算结果和所述虚部计算结果进行数学平均计算，得到高精度频谱校正结果。具体实施例如下：首先利用窗函数w(n)对单频正弦信号进行加权处理，即xw(n)＝x(n)·w(n)。此处单频正弦信号即为短时离散采样信号。常用的H项余弦窗函数可表示为：式(1)中，N为采样点数，ah为窗函数系数。将xw(n)进行DFT或FFT计算，得到变换后频谱函数。对应sw(n)的DtFT(DiscretetimeFourierTransform,离散时间傅里叶变换)为：SW(γ)＝0.5ARms[ejβW(γ-λ)+e-jβW(γ+λ)](2)对于sw(γ)，采用简单的局部峰值搜索程序，即可得到峰值谱线位置索引l。由于采用DFT计算，则对应坐标为l/bin的频谱表示为：SW(l)＝0.5ARms[ejβW(-v)+e-jβW(2l+v)](3)可将式(3)改写为三部分，分别表示为D(γ)、G(γ)和exp(-jπγ)，即：明显D(γ)保持不变，exp(-jπγ)控制相位变化，非线性函数|G(γ)|则直接控制W(γ)的幅值变化。由此坐标为l/bin的离散频谱可改写为：进一步，将中括号中的两项进行展开，对应实部和虚部分别如式(6)所示。式中，κRe＝cos(β+πν)，κIm＝sin(β+πν)。至此，SW(·)中随谱线位置l变化的部分已然非常明显，进一步将式(6)改写为：同理，可得临近谱线的实部和虚部为：为消掉恒定值Ψ的影响，采用三点插值的处理方式。但区别于前述传统IpDFT算法直接采用幅度谱SW(·)进行插值，此处分别对实部和虚部采用插值计算，由此在两点插值基础上进行三点插值，可得如式(7)所示比例。明显，这种类似传统三点插值处理方式可直接消掉恒定值Ψ，同时能够有效规避掉两点插值中近邻次峰值谱线选择错误所带来的巨大校正估计误差。进一步对上式进行如下处理：将上式进一步简化为：式中Γi参量分别为：名下上述分量为近邻频谱间非线性变化部分的线性插值，为此利用G(·)整数间隔间所存在的比例关系，将式(12)进一步化简为：式中，F(x)＝2(2H-1)/(H2-x2)从上式可明显看出，当负频率分量的非线性部分G(·)可忽略不计时，此时根据v＝H·Γ2/Γ1，即vRe和vIm都能作为高精度频谱本文档来自技高网...

【技术保护点】
1.一种高精度频谱校正方法，其特征在于，包括：/n采集短时离散采样信号；/n对所述短时离散采样信号加窗并进行DFT计算，得到变换后频谱函数；/n通过局部峰值搜索所述变换后频谱函数，得到峰值；/n选取峰值所在谱线及相邻的两条谱线作为目标谱线，对所述目标谱线对应的实部和虚部进行拆分，生成实部拆分结果和虚部拆分结果；/n对所述实部拆分结果进行三点插值计算，得到实部计算结果；/n对所述虚部拆分结果进行三点插值计算，得到虚部计算结果；/n对所述实部计算结果和所述虚部计算结果进行数学平均计算，得到高精度频谱校正结果。/n

【技术特征摘要】
1.一种高精度频谱校正方法，其特征在于，包括：
采集短时离散采样信号；
对所述短时离散采样信号加窗并进行DFT计算，得到变换后频谱函数；
通过局部峰值搜索所述变换后频谱函数，得到峰值；
选取峰值所在谱线及相邻的两条谱线作为目标谱线，对所述目标谱线对应的实部和虚部进行拆分，生成实部拆分结果和虚部拆分结果；
对所述实部拆分结果进行三点插值计算，得到实部计算结果；
对所述虚部拆分结果进行三点插值计算，得到虚部计算结果；
对所述实部计算结果和所述虚部计算结果进行数学平均计算，得到高精度频谱校正结果。


2.根据权利要求1所述的高精度频谱校正方法，其特征在于，对所述短时离散采样信号加窗并进行DFT计算，得到变换后频谱函数，包括：

【专利技术属性】
技术研发人员：李波，罗永睦，张林山，周年荣，曹敏，王浩，杨超，付志红，赖强，朱全聪，利佳，邹京希，
申请(专利权)人：云南电网有限责任公司电力科学研究院，
类型：发明
国别省市：云南;53