一种高精度的时频分析方法技术

本发明专利技术提供了一种高精度的时频分析方法，属于数字信号处理领域。所述方法包括以下步骤：(1)输入原始数字信号x(t)，然后对所述原始数字信号x(t)进行EMD分解，将原始数字信号x(t)分解成IMF分量的集合，即其中ci(t)为IMF分量；(2)分别对每个IMF分量ci(t)进行MP分解，得到每个IMF分量的时频谱；(3)将步骤(2)得到的各个IMF分量的时频谱相互叠加，得到原始数据信号x(t)的高时频分辨率的时频谱。本发明专利技术方法避免了直接对复合信号进行MP分解匹配时频原子的不准确性，通过对各IMF分量进行MP分解能够得到准确性更高、时频分辨率更高、信息更丰富的时频谱。

【技术实现步骤摘要】
【专利摘要】本专利技术提供了，属于数字信号处理领域。所述方法包括以下步骤：(1)输入原始数字信号x(t)，然后对所述原始数字信号x(t)进行EMD分解，将原始数字信号x(t)分解成IMF分量的集合，即其中ci(t)为IMF分量；(2)分别对每个IMF分量ci(t)进行MP分解，得到每个IMF分量的时频谱；(3)将步骤(2)得到的各个IMF分量的时频谱相互叠加，得到原始数据信号x(t)的高时频分辨率的时频谱。本专利技术方法避免了直接对复合信号进行MP分解匹配时频原子的不准确性，通过对各IMF分量进行MP分解能够得到准确性更高、时频分辨率更高、信息更丰富的时频谱。【专利说明】
本专利技术属于数字信号处理领域，具体涉及一种针对数字信号的高精度的时频分析方法。
技术介绍
时频分析是分析非平稳信号的传统方法，它能将一维的时域或频域信号变换到二维时频平面上，得到信号的时频谱，可以在时频域中提取并区分信号分量，能够更清楚的了解信号的特征。时频分析方法一般分为线性时频分析方法、双线性时频分析方法和参数化时频分析方法等。典型的线性时频分析方法包含短时傅里叶变换、小波变换、S变换等；由于线性时频分析方法受Heisenberg测不准原理的约束，信号的时间分辨率和频率分辨率为一对矛盾体，不能同时达到较高的时间分辨率和较高的频率分辨率。典型的双线性时频分析方法为Wigner-Ville分布(即WVD);由于WVD中不含任何窗函数，避免了线性时频分析中时间分辨率和频率分辨率矛盾的问题；在分解单分量平稳信号时具有很高的时频分辨率；但是WVD为双线性的，在分解多分量非平稳信号时会出现严重的交叉项。最典型的参数化时频分析方法为匹配追踪法。匹配追踪法为S.Mallat与Z.Zhang于1993年提出的，它是通过对时频原子进行伸缩、时移、频移和相移得到一个原子库，同时根据最大匹配投影原理寻找最佳原子的线性组合，以达到对原始信号进行分解，具有较高的时频分辨率。但是MP算法是一种贪婪算法，它是通过不断地迭代寻找局部最佳原子的方法来实现对原信号的不断逼近过程。为了提高计算效率，YangHua Wang提出了在迭代过程中全局粗粒度预测和局部最优相结合的方法，充分考虑了尺度的重要性，提出了尺度滤波的思想，提高了算法的计算效率。但是对于复合信号，由于波与波相互叠加形成了复合波，直接对复合波进行匹配搜索往往得不到单一的时频原子。如在某一个时间段，彼此间时移较小的几个单波相互叠力口。若直接利用MP算法进行匹配，匹配的时频原子便为一个复合的时频原子，为几个单波相叠加共同作用的波；按照匹配的思想需要继续迭代，需要从原始信号中减去搜索的时频原子，那么剩下的波形便不是一个完整的波形，而是复合波的剩余波，它包含每一个波的部分信息。那么随着不断的搜索，匹配的时频原子越来越失真、不具任何意义。同时，直接利用对非线性的复合信号经希尔伯特变换计算的属性来预测时频原子的时间、频率和相位是没有意义的，不能准确的定位单组份信号的时间和频率等相关信息，导致最终计算的时频谱不准确。瞬时频率为解析信号瞬时相位的导数，为了得到有意义的瞬时频率，E.N.Huang等人对瞬时频率进行了充分的分析，认为只有本征模函数经希尔伯特变换计算的瞬时频率才有意义。在此基础上，提出了一种全新的信号时频分析方法——希尔伯特黄变换(Hilbert-Huang Transform,即 HHT)。HHT 最关键的技术为经验模态分解(Empirical ModeDecomposition,即EMD)，即对时间序列进行经验模态分解，分解成本征模函数组。经验模态分解从本质上讲是对信号进行平稳化的过程，将信号中客观存在的、不同尺度的波动分解出来，得到一系列不同尺度数据序列，即本征模函数。
技术实现思路
本专利技术的目的在于解决目前匹配追踪时频分解中时频原子搜索的准确性问题，提供，充分考虑了 IMF(Intrinsic Mode Function)分量所特有的性质，以及希尔伯特变换对实信号的要求；采取的是对IMF分量进行MP分解(MatchingPursuit Decomposition,即MP分解),能够更准确地定位时频原子的位置、频率、相位及振幅等信息，能够检测出更丰富的频率成分。采用时频分辨率较高的WVD计算时频原子的时频谱，能够得到高时频分辨率的时频谱。本专利技术是通过以下技术方案实现的:所述方法包括以下步骤:(I)输入原始数字信号X (t)，然后对所述原始数字信号X (t)进行EMD分解，将原始数字信号X(t)分解成MF分量的集合，即【权利要求】1.，其特征在于:所述方法包括以下步骤: (1)输入原始数字信号X(t)，然后对所述原始数字信号x(t)进行EMD分解， 将原始数字信号x(t)分解成IMF分量的集合，即 2.根据权利要求1所述的高精度的时频分析方法，其特征在于:所述步骤(1)具体包括以下步骤:(11)初始化，rQ(t)= x(t),令 i = I ; (12)筛选第i个IMF分量，具体包括以下步骤:(12a)初始化，h0(t) = IV1 (t), j = I ； (12b)根据信号V1 (t)的波形，找出IviU)的局部极大值点和局部极小值点的位置和振幅值，并采用边界处理方法分别将所述局部极大值点和局部极小值点向两边延拓1-2个极值点； (12c)利用差值方法对(t)的局部极小值点进行差值，得到IviU)的下包络线；对h^a)的局部极大值点进行差值，得到(t)的上包络线； (12d)计算步骤(12c)得到的上包络线和下包络线的平均值IV1U)； (12e)做减法:hj(t) = hj^ (t) (t)； (12f)如果满足终止条件，则令Ci (t)=比⑴，然后转入步骤(13);否则，j = j+1，转到(12b)； (13)求剩余信号(t) = ri_1(t)-ci(t)； (14)判断ri(t)的极值点的个数是否仍然超过两个，如果是，则i= i+Ι，转到步骤(12)，如果否，则分解结束，ri(t)即为剩余分量。3.根据权利要求2所述的高精度的时频分析方法，其特征在于:所述步骤(2)是这样实现的:首先进行MP分解，即搜索组成一个IMF分量的时频原子，然后计算每个时频原子的时频谱，最后将所有时频原子的时频谱相叠加便得到该IMF分量的时频谱，具体包括以下步骤: (21)初始化，y」(t)= cjt),令 j = I ; (22)对t(t)进行希尔伯特变换，生成复信号，然后分别求解复信号的瞬时振幅、瞬时频率和瞬时相位；其中Un为瞬时振幅最大值对应的时间；fn为该时间对应的瞬时频率；Φη为该时间对应的瞬时相位； 在一组固定的un、fn和Φη情况下，通过计算最优化公式 【文档编号】G01R23/16GK103675444SQ201210316987【公开日】2014年3月26日 申请日期:2012年8月30日 优先权日:2012年8月30日 【专利技术者】陈科, 庞锐, 王鹏燕 申请人:中国石油化工股份有限公司, 中国石油化工股份有限公司石油物探技术研究院本文档来自技高网...

【技术保护点】
一种高精度的时频分析方法，其特征在于：所述方法包括以下步骤：(1)输入原始数字信号x(t)，然后对所述原始数字信号x(t)进行EMD分解，将原始数字信号x(t)分解成IMF分量的集合，即其中ci(t)为IMF分量；(2)分别对每个IMF分量ci(t)进行MP分解，得到每个IMF分量的时频谱；(3)将步骤(2)得到的各个IMF分量的时频谱相互叠加，得到原始数据信号x(t)的高时频分辨率的时频谱。FDA00002080015400011.jpg

【技术特征摘要】

【专利技术属性】
技术研发人员：陈科，庞锐，王鹏燕，
申请(专利权)人：中国石油化工股份有限公司， 中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京;11