自适应时频分析及其研究进展
- 格式:doc
- 大小:161.28 KB
- 文档页数:10
时频分析方法范文时频分析是一种用于分析非平稳信号的方法,它基于时间和频率域的分析技术,能够给出信号在不同时间和频率上的变化规律。
时频分析通常用于处理具有瞬态特征的信号,例如声音、图像、生物信号等。
本文将介绍时频分析的基本原理、常见方法及其在不同领域的应用。
一、基本原理时频分析基于声学和数学等领域的原理,旨在研究信号在时间和频率两个维度上的变化。
传统的傅里叶变换只能提供信号的频域信息,无法描述非定常或非线性信号在时间上的变化。
时频分析通过引入窗函数来实现信号在时间和频率上的分解。
1.窗函数窗函数是时频分析的关键概念,它将信号在时间上切割成多个片段,并将每个片段与一个特定的函数进行乘积。
窗函数通常是时域上的一种窄带滤波器,能够减小信号在时频域的交叉干扰。
常见的窗函数有矩形窗、汉宁窗、高斯窗等。
2.短时傅里叶变换(STFT)短时傅里叶变换是时频分析的最基本方法,它将信号分成多个时间窗口,并对每个窗口进行傅里叶变换。
STFT的窗口长度和重叠率可以根据信号的特性进行调整,从而控制时间和频率分辨率。
STFT分析得到的结果是一个时频矩阵,可以直观地表示信号在不同时间和频率上的能量分布。
3. 维纳-辛钦(Wigner-Ville)分布维纳-辛钦分布是一种时频分析方法,它基于短时傅里叶变换,通过在矩阵的对角线上进行平均来消除交叉干扰。
Wigner-Ville分布能够提供更精确的时频信息,但对噪声和窗口选择比较敏感。
4.小波变换小波变换是一种基于频率域的时频分析方法,它利用小波函数的局部性质,将信号分解成不同频率段的子信号。
小波变换具有良好的时间和频率局部化特性,能够捕捉到信号中的瞬态特征。
常见的小波变换方法有连续小波变换(CWT)和离散小波变换(DWT)。
二、常见方法除了上述方法,时频分析还有一些其他常见的方法,如下所示。
1. 希尔伯特-黄(Hilbert-Huang)变换希尔伯特-黄变换是一种非平稳信号的时频分析方法,它由希尔伯特变换和经验模态分解(EMD)两部分组成。
几种时频分析方法及其工程应用时频分析是一种将时间和频率维度综合起来分析信号的方法,广泛应用于信号处理、通信、音频处理、图像处理等领域。
在实际工程应用中,根据不同的需求和应用场景,可以采用多种不同的时频分析方法。
本文将介绍几种常见的时频分析方法及其工程应用。
短时傅里叶变换是一种将信号分为多个小片段,并对每个小片段进行傅里叶变换的方法。
它在时域上采用滑动窗口的方式将信号分段,然后进行傅里叶变换得到频域信息。
STFT方法具有时间和频率分辨率可调的特点,可用于信号的频域分析、谱估计、声音的频谱显示等。
工程应用:STFT广泛应用于语音处理、音频编解码、信号分析等领域。
例如在音频编解码中,可以利用STFT分析音频信号的频谱特征,进行数据压缩和编码。
2. 小波变换(Wavelet Transform)小波变换是一种时频分析方法,它通过将信号与一系列基函数(小波)进行卷积来分析信号的时间和频率特性。
小波变换具有多分辨率分析的特点,可以在不同尺度上对信号进行分析。
工程应用:小波变换可以用于信号处理、图像压缩等领域。
在图像处理中,小波变换被广泛应用于图像的边缘检测、图像去噪等处理过程中。
3. Wigner-Ville分布(Wigner-Ville Distribution,WVD)Wigner-Ville分布是一种在时间-频率平面上分析信号的方法,它通过在信号的时域和频域上进行傅里叶变换得到瞬时频率谱。
WVD方法可以展现信号在时间和频率上的瞬时变化特性。
工程应用:Wigner-Ville分布在通信领域中被广泛应用于信号的调制识别、通信信号的自适应滤波等方面。
例如在调制识别中,可以利用WVD方法对调制信号的频谱特征进行分析,从而判断信号的调制类型。
4. Cohen类分析(Cohen's class of distributions)Cohen类分析是一种将信号在时间-频率域上进行分析的方法,它结合了瞬时频率和瞬时能量的信息。
时频分析助力精准医疗影像分析一、时频分析与精准医疗影像概述时频分析作为一种强大的信号处理技术,在多个领域展现出了独特的优势。
它能够同时在时间和频率两个维度上对信号进行分析,从而更加全面、细致地揭示信号的内在特征。
在医学影像领域,传统的影像分析方法往往侧重于空间域的信息提取,而时频分析的引入为精准医疗影像分析带来了新的视角和可能性。
精准医疗影像在现代医疗中扮演着至关重要的角色。
它涵盖了多种成像模态,如 X 光、CT、MRI、超声等,这些影像能够为医生提供人体内部结构和生理状态的直观信息。
通过对影像的分析,医生可以诊断疾病、评估病情进展、制定治疗方案以及监测治疗效果。
然而,医学影像数据量巨大、信息复杂,传统分析方法在提取细微特征和准确解读影像方面存在一定的局限性。
二、时频分析在精准医疗影像分析中的应用(一)时频分析技术原理时频分析的核心在于将信号分解为不同时间和频率成分的组合。
常见的时频分析方法包括短时傅里叶变换(STFT)、小波变换(WT)等。
STFT 通过对信号进行加窗处理,在一定时间窗内计算傅里叶变换,从而得到信号在该时间段内的频谱信息。
随着时间窗的滑动,可以获取信号频谱随时间的变化情况。
小波变换则利用小波基函数对信号进行多尺度分解,不同尺度下的小波系数反映了信号在不同频率和时间分辨率下的特征。
这些时频分析方法能够有效地捕捉到信号中的瞬态特征、频率变化以及时变信息,为精准医疗影像分析提供了有力的工具。
(二)不同影像模态中的应用1. CT 影像分析在 CT 影像中,时频分析可用于检测肺部疾病。
例如,在肺部结节的检测中,通过对 CT 影像序列进行时频分析,可以提取结节的纹理特征和动态变化信息。
正常肺部组织和结节组织在时频域上可能呈现出不同的特征模式,时频分析能够增强这些差异,提高结节检测的准确性。
此外,对于肺部疾病的动态监测,如肺炎的发展过程,时频分析可以跟踪肺部影像在不同时间点的频率变化,帮助医生评估病情的演变和治疗效果。
摘要快速傅里叶、Wigner-Ville变换、小波变换在分析非线性非平稳信号时都存在着各自的缺陷与不足。
为了更好地解决这些问题,本文采用了经验模态分解EMD (Empirical Mode Decomposition)来对信号进行分析和滤波。
EMD 分解将复杂信号分解成有限个固有模态函数IMF(Intrinsic Mode Functions)之和,具有很高的频率分辨率和自适应性。
本文对经验模态分解整个理论体系进行了深入的研究,重点研究了EMD 时频分析的应用和基于EMD 的滤波方法。
本文首先引入了瞬时频率的概念,论述了IMF的基本概念和EMD分解算法原理,并且分析了EMD算法的特点,结合IMF给出了边际谱和Hilbert谱的物理含义,对整个EMD 时频分析理论进行了详细的论述。
分析EMD时频分析方法在模态混叠、停止准则等方面存在的不足。
然后通过仿真实验将EMD时频分析方法与传统的时频分析方法进行了比较研究,验证了EMD时频分析方法具有很高的时频分辨率。
最后研究了基于EMD的滤波方法,验证了该方法的有效性及优越性。
关键词:经验模态分解;时频分析;HHT谱ABSTRACTFast Fourier, Wigner-Ville transform and wavelet transform have their own flaws and shortcomings when anglicizing the nonlinear and non-stationary signals.In order to solve these problems, the empirical mode decomposition EMD (Empirical Mode Decomposition) used to signal analysis and filtering in this article. EMD decomposes a complex signal into a finite number of IMF (Intrinsic Mode Functions), and EMD is a high frequency resolution and adaptive method. In this paper, the entire theoretical systems of empirical mode decomposition have been researched deeply, focused on the filtering methods based on EMD and EMD application of time-frequency analysis.In this article, the concept of instantaneous frequency is introduced firstly. Then, the basic concepts of IMF and principles of EMD decomposition algorithm are discussed, and the characteristics of the EMD algorithm are analyzed. Combined with IMF, it is given that the physical meaning of the Hilbert marginal spectrum and the Hilbert spectrum are formulated, and the theories of the whole time-frequency analysis are discussed in detail. Besides, the problem of EMD time-frequency method is also analyzed, including modes mixing, end effect, sifting stop condition and so on. Secondly, compared to traditional time-frequency analysis methods, show that EMD spectrum has perfect time-frequency concentration. Finally, study the EMD-based filtering method, and verified the validity and the superiority of the method.Keywords: EMD; time-frequency analysis; HHT spectrum目录第1章绪论 (1)1.1课题的研究背景和意义 (1)1.2EMD方法的提出 (3)1.3EMD方法的研究发展概况 (4)1.4论文结构与安排 (5)第2章EMD时频分析的基本理论 (7)2.1EMD方法的基本概念 (7)2.2EMD方法的基本原理 (8)2.3基于EMD的HHT时频分析 (12)2.4EMD算法存在的问题及改进 (13)2.5本章小结 (15)第3章EMD 时频分析方法的应用 (16)3.1EMD时频分析仿真及性能比较 (16)3.2EMD在信号趋势提取中的应用 (19)3.3EMD时频分析在奇异信号检测中的应用 (21)3.4本章小结 (24)第4章基于EMD的滤波器设计及仿真 (25)4.1传统的信号滤波和EMD滤波 (25)4.2EMD的二进滤波特性及仿真 (26)4.3基于EMD的滤波方法及仿真 (28)4.4本章总结 (30)结论 (31)参考文献 (33)致谢 (35)附录 (36)第1章绪论1.1 课题的研究背景和意义信号往往包含着许多的重要信息,如时间特征、频率特征等。
地震流动观测数据的时频分析技术研究1. 引言地震是一种自然现象,它在地球内部的断层平面上产生的震动。
由于地震的突发性和破坏性,对于地震的预测和监测是至关重要的。
地震流动观测数据是指在地震事件发生过程中采集到的连续波形数据,在地震研究中起着重要的作用。
本文旨在研究地震流动观测数据的时频分析技术,以提高对地震事件的监测和分析能力。
2. 地震流动观测数据的时频分析概述时频分析是一种将信号在时间和频率上进行分析的方法,可以揭示信号的时变特性。
在地震流动观测数据的时频分析中,常用的方法有短时傅里叶变换(STFT)、连续小波变换(CWT)和经验模态分解(EMD)等。
这些方法可以通过将信号分解成不同频率和时间域分量,从而帮助我们理解地震事件的时间和频率特征。
3. 短时傅里叶变换(STFT)短时傅里叶变换是一种将信号分解为不同频率分量的传统方法。
它将整个信号分成若干个片段,对每个片段进行傅里叶变换。
通过在时域上进行窗函数的滑动,可以得到信号在不同时间和频率上的变化。
然而,STFT存在时间和频率分辨率的局限性,无法兼顾时间和频率的精确描述。
4. 连续小波变换(CWT)连续小波变换是一种基于小波理论的时频分析方法。
与STFT不同,CWT采用可变尺度的小波基函数,可以在不同时间和频率上进行分析。
CWT将信号分解为多个频带,每个频带包含丰富的时频信息。
这使得CWT在处理地震流动观测数据时具有较好的灵活性和适应性。
5. 经验模态分解(EMD)经验模态分解是一种自适应的时频分析方法,可以将信号分解成一组本征模态函数(IMF)。
在地震流动观测数据的时频分析中,EMD可以将信号分解为包含不同时频信息的IMF分量,并能够准确揭示出信号的本质时频特征。
通过分析信号的IMF分量,我们可以获得关于地震事件的重要信息。
6. 地震流动观测数据的时频分析应用时频分析技术在地震研究中具有广泛的应用。
首先,时频分析可以用于地震源过程的研究,揭示地震发生的时间和频率特征。
《EMD时频分析方法在旋转机械耦合故障诊断中的应用研究》篇一一、引言随着工业技术的不断发展,旋转机械在生产过程中的作用日益显著。
然而,由于各种复杂因素的影响,旋转机械常会出现耦合故障,对生产效率和产品质量带来严重影响。
因此,对于旋转机械耦合故障的诊断和识别成为了重要的研究方向。
时频分析作为一种有效的信号处理方法,近年来在旋转机械故障诊断中得到了广泛应用。
其中,EMD(Empirical Mode Decomposition,经验模态分解)时频分析方法以其优秀的自适应性和稳定性在旋转机械耦合故障诊断中显示出其独特的优势。
本文旨在研究EMD时频分析方法在旋转机械耦合故障诊断中的应用,以提高诊断效率和准确性。
二、EMD时频分析方法EMD是一种基于信号本身的自适应时频分析方法,可以有效地处理非线性、非平稳信号。
该方法通过将信号分解为多个本征模态函数(Intrinsic Mode Functions,IMF),进而实现对信号的时频分析。
EMD方法具有自适应性强、抗干扰能力强、无需预设基函数等优点,因此在信号处理领域得到了广泛应用。
三、EMD在旋转机械耦合故障诊断中的应用旋转机械在运行过程中,由于各种因素的影响,可能会出现多种耦合故障。
这些故障往往表现为复杂的振动信号,传统的频域分析方法难以准确诊断。
而EMD时频分析方法可以通过分解振动信号,提取出各个本征模态函数,从而实现对故障的精确诊断。
首先,通过传感器采集旋转机械的振动信号,然后利用EMD 方法对振动信号进行分解,得到多个IMF分量。
每个IMF分量代表了一种振动模式,通过对IMF分量的分析和处理,可以识别出旋转机械的耦合故障类型和程度。
在实际应用中,EMD时频分析方法可以根据故障类型和程度,对IMF分量进行进一步的特征提取和模式识别。
例如,对于轴承故障,可以通过分析IMF分量中的频率成分和幅度变化来诊断故障类型和严重程度;对于齿轮箱故障,可以通过分析IMF分量中的相位关系和能量分布来识别故障源和类型。
时频分析方法在心电信号处理中的应用心电信号是一种记录心脏电活动的生物信号,对于心脏疾病的诊断和治疗具有重要的意义。
随着科技的进步,时频分析方法逐渐被应用于心电信号的处理与分析中。
本文将介绍时频分析方法在心电信号处理中的应用,并探讨其在心脏疾病诊断中的潜力。
一、时频分析方法的原理及常用算法时频分析方法是将信号在时间和频率两个不同的域中进行分析的方法。
其原理是通过对信号进行时频变换,可以获取信号在不同时间和频率上的能量分布特征,进而揭示信号的时频特性。
常用的时频分析算法包括短时傅里叶变换(STFT)、小波变换(Wavelet Transform)和Wigner-Ville变换(Wigner-Ville Distribution)。
这些算法各有特点,可以对信号在时频域上进行不同粒度的分析,从而满足不同应用场景的需求。
二、时频分析方法在心电信号处理中的应用1. 心电信号去噪心电信号常常伴随着各种干扰和噪声,如肌电干扰、基线漂移等。
传统的滤波方法对于非平稳信号可能存在较大的限制。
而时频分析方法能够根据信号时频特性的差异,对信号进行自适应的滤波处理,抑制噪声并保留有用的心电信息。
2. 心律失常检测心律失常是心脏疾病的常见病症之一,对其准确诊断和监测有着重要意义。
时频分析方法可以通过提取心电信号在不同时频域上的特征,实现对心律失常的检测和分类。
例如,可以通过小波变换分析心电信号的R波波形,从而判断是否存在心律失常。
3. 心脏病诊断时频分析方法还可以应用于心脏病的诊断。
通过对心电信号在时频域上的特征进行分析,可以辅助医生判断患者是否存在心脏异常,如心肌梗死、心肌缺血等。
同时,时频分析方法还可以对不同类型的心脏病进行分类,为治疗方案的选择提供依据。
三、时频分析方法在心电信号处理中的局限性与挑战虽然时频分析方法在心电信号处理中具有广泛的应用前景,但仍存在一些局限性和挑战。
1. 算法复杂性常用的时频分析算法,如小波变换和Wigner-Ville变换,对于复杂信号的处理需要较高的计算复杂性,这对于实时应用而言可能存在一定的挑战。
自适应核时频matlab
自适应核时频Matlab是一种用于信号处理和通信系统设计的强大工具。
它可以自动调整参数以适应不同的信号特性和环境条件,以确保系统性能最优化。
本文将介绍自适应核时频Matlab的基本原理和应用。
自适应核时频Matlab基于时频分析的理论,可以更好地捕捉信号的瞬时特征和频率特性。
它通过调整核函数的大小和形状,以适应不同信号的时频结构,从而提高信号处理的准确性和效率。
在Matlab中,可以使用自适应核时频分析工具箱来实现自适应核时频分析。
该工具箱提供了多种核函数和参数选择的方法,用户可以根据具体需求选择合适的核函数和参数,以获得最佳的时频表示结果。
自适应核时频Matlab在许多领域都有广泛的应用。
例如,在语音信号处理中,可以利用自适应核时频Matlab来提取语音信号的时频特征,实现语音识别和分割;在雷达信号处理中,可以利用自适应核时频Matlab来实现目标检测和跟踪等功能。
除了信号处理领域,自适应核时频Matlab还可以应用于图像处理、生物医学信号分析、金融时间序列分析等领域。
它的灵活性和高效性使其成为许多工程和科学领域的首选工具之一。
总的来说,自适应核时频Matlab是一种强大的信号处理工具,它可
以自动调整参数以适应不同信号的时频特性,提高信号处理的准确性和效率。
它在各个领域都有广泛的应用,为工程师和科学家们提供了一个强大的分析工具,帮助他们更好地理解和处理复杂的信号数据。
时频分析方法时频分析是一种用于研究信号在时间和频率两个维度上变化规律的方法。
它在信号处理、通信系统、地震学、生物医学工程等领域都有着广泛的应用。
本文将介绍时频分析方法的基本原理和常见的分析技术,希望能为读者提供一些帮助。
时频分析的基本原理是将信号在时间和频率上进行分解,以揭示信号在不同时间段和频率段的特征。
在时域上,我们可以观察信号的波形和振幅变化;在频域上,我们可以得到信号的频谱信息。
时频分析方法的目的就是将这两个维度结合起来,得到信号在时间和频率上的特性。
常见的时频分析方法包括傅里叶变换、小波变换、时频分布等。
傅里叶变换是一种将信号从时域转换到频域的方法,可以得到信号的频谱信息。
小波变换是一种同时在时域和频域上进行分析的方法,可以更好地捕捉信号的瞬时特性。
时频分布则是一种将信号的时频特性可视化的方法,常用的有Wigner-Ville分布和短时傅里叶变换等。
在实际应用中,选择合适的时频分析方法取决于信号的特性和分析的目的。
如果信号具有明显的频率成分,可以选择傅里叶变换来观察频谱信息;如果信号具有瞬时特性,可以选择小波变换来捕捉信号的瞬时变化;如果需要同时观察信号的时频特性,可以选择时频分布来进行分析。
除了选择合适的时频分析方法,还需要注意信号的预处理和参数的选择。
对于非平稳信号,需要进行平滑处理或者选择适当的小波基函数;对于时频分布方法,需要选择合适的窗口长度和重叠率来得到准确的时频信息。
总之,时频分析是一种重要的信号分析方法,可以帮助我们更好地理解信号的时频特性。
在实际应用中,我们需要根据信号的特性和分析的目的选择合适的时频分析方法,并注意信号的预处理和参数的选择,以得到准确的分析结果。
希望本文能对读者有所帮助,谢谢阅读!。
引力波调制信号的统计时频分析引言:引力波是由于质量的突变或者大质量天体的运动而产生的空间弯曲,它们以波动的形式传播。
近年来,科学家们利用先进的引力波探测器成功探测到了许多引力波信号。
然而,引力波信号往往包含大量的噪声,所以准确地分析和处理这些信号十分具有挑战性。
本文将介绍引力波调制信号的统计时频分析方法,以帮助读者更好地理解和处理这些信号。
1. 引力波调制信号的基本概念引力波调制信号是指由引力波作用于探测器所产生的信号。
这些信号受到许多因素的影响,例如引力波源的位置、质量和运动方式等。
因此,我们需要对这些信号进行统计时频分析,以了解它们的特征和生成机制。
2. 时频分析方法介绍2.1 短时傅里叶变换(Short-time Fourier Transform)短时傅里叶变换是一种经典的时频分析方法,它通过对时域信号的不同时间窗口进行傅里叶变换,得到信号在频域上的分布。
然而,由于引力波信号往往具有较长的持续时间,短时傅里叶变换在提供高频率分辨率的同时,对低频信号的时域分辨率较差。
2.2 小波变换(Wavelet Transform)小波变换是一种基于变尺度分析的时频分析方法,它通过对信号进行连续的尺度变换和平移,得到信号在不同尺度上的频谱信息。
相比于短时傅里叶变换,小波变换可以提供更好的时频分辨率,特别是对于低频信号的分析。
2.3 自适应时频分析方法自适应时频分析方法是一种根据信号的局部特性自动调整时频分辨率的方法。
它结合了短时傅里叶变换和小波变换的优点,在分析引力波调制信号时具有较好的适应性和灵活性。
常见的自适应时频分析方法包括经验模态分解(Empirical Mode Decomposition)和可变窗口谱估计(Variable Window Spectrum Estimation)等。
3. 引力波调制信号的统计特征引力波调制信号具有一定的统计特征,对了解信号的特性和噪声的影响至关重要。
通过对大量引力波数据的分析,科学家们发现了一些重要的统计特征,如信号的功率谱密度、相关性和自相关函数等。
数字信号处理学号:************学生所在学院:测试与光电工程学院学生姓名:XXX********教师所在学院:测试与光电工程学院2013年12月13级4班自适应时频分析及其研究进展XXX(南昌航空大学测试与光电工程学院南昌330063)摘要:通过对自适应时频分析的发展历程的了解,总结分析近几年内,各学科对于自适应时频分析方法的具有创新的研究进展,表现出其在各领域研究中的不可忽视的地位。
对不同的时频分析技术作了简要介绍,并对其优缺点、彼此间的相互关系进行了较为详细的论述。
同时,对时频分析技术所面临的问题及发展方向谈了一点个人看法。
关键词:傅氏变换;适应信号分解;Matching pursuit算法;时频分析引言尽管信号处理的目的不同,如用于数据压缩、降噪、检测、参量估计及模式识别等,但信号处理的基本步骤和方法却是一致的。
那就是首先要获取信号的特征信息。
对待分析信号作各种变换处理的根本目的,就是要通过变换处理使待分析信号的特征信息尽可能地突显出来以利于特征提取。
尽管人们已经提出了形形色色的信号变换方法,但大体上却可分为如下3类:线性变换方法,双线性变换方法及参数化时频表示方法。
本文就现有的各种方法作一综述和比较,并就其中的一些问题谈一点个人的粗浅看法。
1自适应时频分析的现状时频分析方法提供了时间域和频率域的联合分布信息,更清楚地描述了信号的频率是如何随时间变化的关系。
时频分析的方法很多,从短时傅立叶变换到二次型时频分析、Gabor变换、Cohen类时频分布等,各类分布多达几十种,但这些基于传统理论的各种处理方法,有着种种自身难以克服的缺陷。
短时傅里叶分析方法的不足之处是不能同时获得高的时频分辨率。
二次型时频分布具有最高的时频分辨率,但其固有的缺陷是交叉干扰项的存在,使用固定核函数对于某一类信号可能较好地抑制了交叉项,但同一种核函数对另外的信号却常常不能尽如人意。
Gabor 展开的基函数可以体现信号的时频特征,但Gabor 基的频率、带宽是固定的。
同样,小波变换的时移、频移也是固定变化的,只是对时频平面进行了机械的格型分割,其基函数不能自适应选取[1]。
为了克服传统时频分析的种种缺陷,到二十世纪九十年代,自适应时频正式登上了历史舞台。
自适应时频分析是一种应用自适应信号处理手段来完成信号的时频分析以获得更好的时频信息描述与跟踪效果的方法。
在所有的时频分析方法中最受关注的是Cohen 类时频分布[2],包括:Wigner 分布[3],优化平滑Wigner 分布[4],Choi-Willimams 分布[5],锥形核分布[6],交叉项检测的时频表示(CDR)[7],减少交叉顶的分布(RID)[8],具有Besell 核的时频分布[9],具有组合核的时频分布[10]。
基于信号分解的时频分析方法是自适应时频分析方法中的一种,它采用的是信号线性空间扩展的手段,从过完备基函数集合中寻求最佳的信号表示。
如匹配搜索原子分解(MP),自适应高斯基表示(AGR)等。
这类方法实际上参数化和模型化的方法,它结合了线性和二次型时频表示的有点,虽然计算较为复杂,但是当参数模型确实匹配被分析信号时,可以取得比非参数型方法更好的效果。
由于自适应方法潜在的优异性能,引起了人们的广泛关注,形成了时频分析研究的一个新热点。
2 自适应时频分析的各种变换2.1 短时傅氏变换短时傅氏变换[11]定义为*ˆ(,)()()j s t s g t e d ωτωτττ+∞--∞=-⎰ (2.1)其中g (t )为一时间局部化了的窗函数。
显然对应于一定的时刻,短时傅氏变换只对其附近窗口内的信号作分析,能够粗略地反映信号在该时刻附近的局部频谱特征。
而整个变换结果也就能提示信号频谱随时间的演化特性,性能优于傅氏变换。
能够保证该方法有效性的一个隐含假设是,信号在窗函数的有效持续时间内应是平稳的,但此条件通常无法满足或近似满足。
上文中所示信号作傅氏变换等价于在特定时刻对该信号作了一次短时傅氏变换,其所用窗函为一矩形窗,宽度等于信号的有效持续时间。
显然,由于不满足局部平稳性假设,短时傅氏变换给了人们错误的信息。
为满足局部平稳性假设条件,人们又采用自适应的方法对不同的信号段选择不同长度的窗函数[12]。
但如何衡量信号的局部平稳性本身就是一大难题。
再者,同傅氏变换一样,即使信号在窗函数的有效持续时间内是平稳的,它也无法较好地给出信号的局部频谱信息。
2.2 小波变换小波变换定义为*(,)()()s t WT t a s d a ττψτ+∞-∞-= (2.2)其中a > 0 为尺度参数,y(t) 为一个时间、频率均局部化了的带通函数。
对复小波而言 ,母小波y(t)可表示成以下形式()()j t t h t e ωψ= (2.3) 其中w 为小波y (t ) 的中心频率, h(t) 为一低通实函数。
则小波变换又可表示成式(2.4)形式*(,)()()j t j t a as t WT t a e s h e d a ωωτττ-+∞-∞-= (2.4)比较两式可见,小波变换和短时傅氏变换具有很大的相似性。
区别仅在于观察信号的不同频率分量,小波变换使用了不同宽度的窗函数。
同短时傅氏变换相比,小波变换具有多分辨能力[13]。
目前,小波变换是信号处理领域的一个研究热点。
不仅关于它的理论及应用方面的文献浩繁,而且它也被罩上了一层神秘的光环,似乎小波变换无所不能。
其实不然,它只适于分析自相似信号。
在实际应用中人们发现,母小波选择的恰当与否至关重要,几乎是影响小波变换应用成败的决定性因素。
其根本原因就在于此。
3 自适应时频分析的各种分布3.1 径向高斯核时频分布G.B.Richard 和L.J.Douglas 在1993年提出了基于信号的径向高斯核时频分析,这种方法将待求的核函数定义为沿任意径向剖面部都是Gauss 型的二维函数,即22(,)exp 2()τθθτσ⎫⎧+⎪Φ=-⎨⎬ψ⎪⎩⎭(3.1) ()σψ为在径向角arctan()τψθ=上的Gauss 函数的均方差,它控制Gauss 函数的形状,称展形函数[14]。
利用极坐标r =作为半径变量,(3.1)式可以写成2(,)exp 2()r r φσ⎧⎫ψ=-⎨⎬ψ⎩⎭(3.2)且核函数满足约束条件22001(,)2r rdrd πφαπ∞ψψ≤⎰⎰ (3.3) 因为径向Gauss 核函数完全由一维函数()opt σψ所控制,寻找优化的径向Gauss 核函数,等价于寻找优化的展形函数()opt σψ,所以(3.3)式可以等价写为01()2d πσαπψψ≤⎰ (3.4)这种方法将二维核函数的优化问题简化为一维展形函数的优化问题,使优化函数的设计向前推进了一步。
但是它仍基于信号的自项成分集中在模糊平面的原点附近、交叉项远离原点这样的理论基础,所以,当自项成分和交叉项成分重叠时,无论体积参数α取何值,都不能将自项和交叉项分开[15]。
3.2 自适应锥形核时频分布N.C.Richard 提出了自适应锥形核设计,根据信号属性自适应调节锥形核的长度。
锥形核分布是一种Cohen 类[16]二次时频表示,其核函数在时间相关平面呈锥形。
它在提供好的时频分辨率和抑制交叉项的同时保留了时间支持特性,在频率轴上能够使谱峰得到增强。
锥形核分布的离散形式[17-18]为:(,)(,)(,)jkw T kCKD n w n k R n k e -=Φ∑ (3.5)其中*(,)()()k p k R n k x n p k x n p k =-=+++-∑ (3.6) 224(),k T T k e T k T φ-=-≤≤ (3.7)式中x(n)(n=0,1,…,N-1)为待分析信号,CKD 的性质决定于锥形长度T 。
为得到好的时频分布效果,对于变化较快的信号,通常选取较小的T;而对于变化较平缓的信号,常选取较大的T 。
也就是说,对于某一信号,过大或过小的了都会降低CKD 的时间一频率分辨率,使其时频能量分布分散。
因此以时频分布的能量集中程度为准则,确定出最佳的了值。
信号的能量为2212*()()()()k n N T n jk T T m n N k T p k e k x m p k x m p k e d ωωω+--=-=-=-=Φ+++-∑∑∑⎰ 221*2()()()kn N T jk T m n N k Tp k k x m p k x m p k e ω+--=-=-=-=Φ+++-∑∑∑ 2212*()()()()k n N T jk T T k T m n N p k k k x m p k x m p k e ω+--=-=-=-⎡⎤⎢⎥=ΦΦ+++-⎢⎥⎣⎦∑∑∑ (3.8) 如果记221*2()()()k n N jk n m n N p k E k x m p k x m p k e ω+--=-=-=+++-∑∑ (3.9)则22()()()T n T T n k T e k E k =-=Φ∑ (3.10)采用如下优化准则2max ()()()T T T n T T k T C n k E k A =-⎡⎤=Φ⎢⎥⎣⎦∑ (3.11) 其中224T k T T k T A e -=-=∑为核函数能量。
()T C n 表示n 时刻时频分布的能量与核函数能量之比,使()T C n 最大的T 即为最优锥形核长度,记为n opt T 。
从而得到自适应锥形核时频分布为2214()*(,)()()n opt n opt n opt T k k T jk p k k T CKD n e x n p k x n p k e ωω---=-=-⎡⎤=+++-⎢⎥⎢⎥⎣⎦∑∑ (3.12) 3.3 自适应最优核时频分布定义具有时频局部化性质的模糊函数[19-20]如下(,,,)(,)(,)j u jv AF t WVD u v A u t v e dudv θτθτωω-+=•--⎰⎰ (3.13)式中,变量t 和ω给出了二维窗函数的中心位置,将Wigner-Ville 分布与窗函数A(u-t,v-ω)相乘,可以将积分区域限制t 和ω的邻域内。
与加窗傅里叶变换相同,窗函数的选取要保证其既要具有时频域局部化性质,也要具有模糊域局部化性质。
虽然采用加窗的方法降低了模糊域的分辨率,但是对最优化核函数影响并不大。
式中的窗函数要求同时在时频域和模糊域具有较好的分辨率。
在定义局部模糊函数以后,时频自适应最优核函数的设计转换为如下的最优化问题max ()opt g ΦΦ=Φ (3.14) 2()(,;,)(,;,,)g AF t t s d d θτωθτωθτΦ=Φ⎰⎰ 2200(,;,)(,;,,)p AF r t r t s rdrd πϕωϕωϕ∞=Φ⎰⎰ (3.15) 从而得到时频自适应最优核时频表示(,)(,;,)(,;,,)j t j opt p t AF t t s e d d θωτωθτωθτωθτ-=Φ⎰⎰ (3.16)纵观自适应时频表示的发展历程,无一不是以追求完美地刻画复杂多样的信号时频特征为动因,以增加算法的复杂性和降低速度为代价。