[发明专利]一种基于对数正态滤波器的局部频率估计方法

摘要

本发明公开了一种基于对数正态滤波器的局部频率估计方法，涉及数字信号处理技术领域，本发明中提出了新的窄带估计结果的表达式：和新的计算宽带估计结果公式原有技术只能在α＝1的情况下进行窄带估计，使得当滤波器带宽设定后，用来覆盖0～π频段的滤波器组中心频率比值关系也就不能改变，使得滤波器估计带宽和滤波器个数相互制约，本发明从理论上突破了这一限制，使得可以保持滤波器估计带宽不变的情况下，采用更多滤波器个数，从而可以最大限度的发挥对数正态滤波器的优越性，从而实现了同时提高基于对数正态滤波器的局部频率估计算法的灵敏度和抗噪声性能的有益效果。

主权项

1.一种基于对数正态滤波器的局部频率估计算法，其特征在于，该方法包括以下步骤：步骤1、获得需要分析局部频率性质的待分析信号SA，根据灵敏度需求设定估计带宽B，将整个0到π频率范围用N个中心频率顺序分割，第i个中心频率为ρi，i＝1,2,..........,N，且使得相邻两个频率范围内的中心频率ρi+1/ρi为定值，确定参数，β＝1/log2(ρi+1/ρi)，β为相邻中心频率的对数倒数比例关系系数；步骤2、将S_A分别通过由中心频率ρ_i确定的传递函数为的对数正态滤波器，获得N个输出结果，其中第i个输出结果为S_Ai；步骤3、将0～360度方向顺序等分为M块，其中第k个块的方向表示为且该方向为表达式为的方向滤波器的方向，将步骤2所得到的每一个S_Ai分别送给M个表达式为的方向滤波器，得到M个输出结果；其中第k个输出结果为k＝1,2,3,...,M，其中是空间向量场u与第k个方向滤波器方向的夹角；步骤4、将步骤3中每个S_Ai对应所得到的M个输出结果组合计算，得到中心频率为ρ_i时的输出结果步骤5、取相邻的ρ_i+1和ρ_i为一组，利用步骤4计算得到的N个f_Ai输出结果计算各组窄带估计结果，共计得到N‑1组窄带估计集结果，其中第h个窄带估计结果的表达式为其中h＝i＝1,2,......N‑1，α是由B和β共同确定的参数，且α＝B²β/8；步骤6、由步骤5所得的窄带估计结果ωh计算宽带估计结果W,其中γ是用来调节加权系数(||f_Aif_A(i+1)||)的参数，且γ＝α。