功率谱分析
- 格式:doc
- 大小:222.50 KB
- 文档页数:10
精确测量脑电波频率方法比较脑电波频率是脑部神经活动的表征之一,它可以通过测量脑电图(EEG)来进行分析和研究。
精确测量脑电波频率对于了解脑部活动的功能与疾病状态具有重要意义。
本文将对几种常用的脑电波频率测量方法进行比较,包括波峰法、功率谱分析法和小波变换方法。
1. 波峰法波峰法是一种最常见的脑电波频率测量方法之一。
它通过检测脑电图信号中波峰所出现的时间间隔来计算频率。
尽管在某些情况下,波峰法可以提供可靠的结果,但它有一定的局限性。
首先,波峰法需要一个精确的起始点和终止点,人为地选择这些点可能带来主观偏差。
其次,如果脑电图信号存在噪音或频率不稳定现象,如频率跳变或谐波等,波峰法的测量结果可能会受到严重影响。
因此,波峰法在精确测量脑电波频率方面存在一定的局限性。
2. 功率谱分析法功率谱分析法是一种常用的脑电波频率测量方法。
它通过将脑电图信号转换为频域信号,然后计算信号在不同频率上的功率密度来获得频率信息。
功率谱分析法具有较高的准确性和可靠性。
通过进行窗函数选择、滑动窗口和傅里叶变换等处理,可以有效地分析不同频率带的脑电波。
然而,功率谱分析法也有其限制。
当脑电波信号中存在高噪声干扰时,功率谱分析法可能会在高频和低频端出现伪迹。
此外,频谱估计的分辨率也可能受到样本长度和窗函数选择等因素的影响。
3. 小波变换方法小波变换方法是一种较为新颖的脑电波频率测量方法。
与传统的傅里叶变换方法相比,小波变换方法具有时间-频率局部性的优点。
它可以捕捉到脑电波信号在不同时间尺度上的频率变化。
小波变换方法在研究不同频率带脑电波活动时具有更好的灵敏度和分辨率。
然而,小波变换方法也需要选择适当的小波函数和尺度,所以对于初学者而言,使用小波变换方法进行脑电波频率测量可能需要一定的学习和实践。
综合比较上述三种测量方法,波峰法简单易行,但在脑电波频率分析中存在较大限制。
功率谱分析法具有较高的准确性和可靠性,但可能受到噪声干扰和频谱估计分辨率的限制。
功率谱分析及其运用简答题一、功率谱分析的基本原理功率谱分析的基本思想是将一个连续时间的信号转换为频域上的离散信号,然后对这些离散信号进行傅里叶变换,得到其频谱表示。
频谱表示中的每个峰值代表了一个特定的频率分量,而每个峰值的高度则代表了该频率分量的强度。
通过对频谱表示进行加权平均,可以得到原始信号的能量分布情况。
二、功率谱分析的应用场景1.通信系统:在无线通信系统中,功率谱分析可以用来检测干扰信号或者识别出合法的通信信号。
通过比较接收到的信号与已知的噪声信号之间的功率谱差异,可以判断出是否存在干扰。
此外,功率谱分析还可以用来估计信道容量和误码率等重要参数。
2.音频处理:在音频处理中,功率谱分析可以用来提取音乐中的基音和谐波等信息。
通过对音乐信号进行快速傅里叶变换(FFT),可以得到其频谱表示,然后再通过滤波器等算法提取出所需的信息。
3.雷达系统:在雷达系统中,功率谱分析可以用来检测目标反射回来的信号。
通过对反射回来的信号进行功率谱分析,可以确定目标的位置、速度和形状等信息。
三、实际运用举例下面以一个简单的示例来说明功率谱分析的实际运用过程。
假设我们有一个包含多个正弦波成分的信号x(t),我们需要将其分解成若干个简单的正弦波成分y(i),并计算每个成分的振幅和频率。
具体步骤如下:1.对信号x(t)进行快速傅里叶变换(FFT),得到其频域表示f (k)。
2.对频域表示f(k)进行平滑处理,以减少高频噪声的影响。
常用的平滑方法包括均值滤波和中值滤波等。
3.对平滑后的频域表示f(k)进行平方运算,得到其功率谱密度ρ(f)。
4.根据需要,可以选择不同的窗函数对ρ(f)进行加窗处理,以减少频谱泄漏等问题。
常见的窗函数包括汉宁窗、汉明窗和矩形窗等。
5.最后,根据ρf)的大小和位置等信息,可以确定原始信号中包含的各个正弦波成分以及它们的振幅和频率等特征。
功率谱分析的原理及应用1. 什么是功率谱分析功率谱分析是一种对信号进行频域分析的方法,它可以将信号在频域上表达出来。
通过功率谱分析,我们可以了解信号的频率分布,并从中提取出信号的特征。
功率谱分析广泛应用于信号处理、通信系统、声学分析等领域。
2. 功率谱分析的原理功率谱分析的原理基于傅里叶变换的思想,将时域上的信号转换为频域上的信号。
傅里叶变换可以将一个信号表示为多个不同频率的正弦波的叠加,而功率谱则表示不同频率正弦波的能量分布情况。
功率谱分析的具体步骤如下:- 第一步:将原始信号转换为时域上的离散信号。
- 第二步:对离散信号进行傅里叶变换,得到频域上的信号。
- 第三步:计算频域上信号的幅度谱,得到信号在不同频率上的能量分布。
- 第四步:对幅度谱进行平方处理,得到功率谱。
3. 功率谱分析的应用功率谱分析在许多领域中都有广泛的应用,以下列举了一些常见的应用场景。
3.1 信号处理功率谱分析在信号处理中具有重要的作用。
通过分析信号的功率谱,我们可以了解信号的频率特性,从而帮助我们对信号进行滤波、降噪等处理。
同时,功率谱分析还能够帮助我们检测信号中的周期性成分,并进行信号的识别和分类。
3.2 通信系统在通信系统中,功率谱分析可以用于频谱分析和带宽分配等任务。
通过对信号的功率谱进行分析,可以确定频率段的使用情况,从而辅助我们进行频谱规划和频率资源的分配。
此外,功率谱分析还可以帮助我们评估信道的质量,从而对通信系统进行优化。
3.3 声学分析声学分析是功率谱分析的另一个重要应用领域。
在声学分析中,功率谱分析可以用于声音信号的频谱分析和特征提取。
通过分析声音信号的功率谱,我们可以了解声音的频率成分和能量分布,进而帮助我们进行声音信号的分类、识别和音频处理等任务。
3.4 振动分析功率谱分析在振动分析中也得到了广泛的应用。
通过对振动信号进行功率谱分析,我们可以了解结构物的固有频率和振动模态,从而帮助我们识别结构物中存在的故障和缺陷。
信号处理的功率谱分析(一)信号处理的功率谱分析(一)信号处理的功率谱分析是一种常用的信号处理技术,它可以对信号的频率特征进行分析和研究。
功率谱分析主要用于确定信号在不同频率上的能量分布情况,进而了解信号的频域特性和频谱结构。
在实际应用中,功率谱分析广泛应用于噪声分析、通信系统性能分析、振动信号分析等领域。
功率谱是指信号在不同频率区间上的能量分布情况。
在信号处理中,一般使用离散傅立叶变换(Discrete Fourier Transform,DFT)来计算信号的功率谱。
DFT是傅立叶变换的一种离散形式,将连续时间域信号转换为离散频率域信号。
通过DFT的计算,可以得到信号在不同频率上的幅度和相位信息,进而计算出信号在不同频率区间上的功率谱。
在进行功率谱计算时,首先需要将原始信号进行采样,得到离散时间序列。
然后,对时间序列进行DFT计算,得到信号的频域表达。
最后,通过对频域表达的幅度进行平方运算,得到信号的功率谱。
功率谱分析可以帮助我们了解信号的频率成分和能量分布情况。
通过功率谱分析,我们可以估计信号的主要频率、频率分布范围和功率集中情况,有助于判断信号的特定特征和性质。
例如,在噪声分析中,功率谱分析可以帮助我们确定噪声的频率成分和功率密度,从而判断噪声的类型和影响。
对于实时信号处理和大数据处理,功率谱分析也有着重要的应用价值。
在实时信号处理中,可以通过连续采样和时域滑动窗口的方式,实时计算信号的功率谱,实现对信号的频域特征的实时监测和分析。
在大数据处理中,可以通过对信号进行分块采样和并行计算,从而加快功率谱分析的速度和效率。
此外,功率谱分析还可以与其他信号处理技术相结合,进一步提高信号处理的效果。
例如,可以将功率谱分析与滤波技术相结合,实现对特定频段的信号抑制和增强;还可以将功率谱分析与自适应算法相结合,实现对非平稳信号的频谱跟踪和估计。
综上所述,功率谱分析是一种常用的信号处理技术,它可以对信号的频率特征进行分析和研究,帮助我们了解信号的频域特性和频谱结构。
信号的功率谱分析1、功率谱密度函数的定义对于随机信号)(t x ,由于其任一样本函数都是时间的无限的函数,一般不能满足傅里叶变换的存在条件(即积分⎰∞∞-dt t x )(必须收敛)。
如果将样本函数取在一个有限区间]2,2[T T -内,如图所示,令在该区间以外的0)(=t x ,则积分⎰∞∞-dt t x )(收敛,满足傅里叶变换条件,变换后用功率谱密度函数表示。
2、功率谱密度函数(又称功率谱)的物理意义是在频域中对信号能量或功率分布情况的描述。
功率谱表示振动能量在频率域的分解,其应用十分广泛。
功率谱的横坐标是频率,纵坐标是实部、虚部的模的平方。
功率谱密度函数作为随机信号在频域内描述的函数。
对于随机信号而言,它不存在频谱函数,只存在功率谱密度函数(功率大小在频谱中反映为频谱的面积)。
时域中的相关分析为在噪声背景下提取有用信息提供了途径。
功率谱分析则从频域提供相关技术所能提供的信息,它是研究平稳随机过程的重要方法。
3.功率谱密度函数的应用(1)结构各阶固有频率的测定 工程结构特别是大型结构(如高层楼房、桥梁、高塔和重要机械设备等)要防止共振引起的破坏,需要测定其固有频率。
如果对结构加以激励(或以大地的脉动信号作为激励信号),即可测定结构的响应(振动信号),再对响应信号作自功率谱分析,便可由谱图中谱峰确定结构的各阶固有频率。
(2)利用功率谱的数学特点求取信号传递系统的频率响应函数。
(3)作为工业设备工作状况的分析和故障诊断的依据 根据功率谱图的变化,可以判断机器设备的运转是否正常。
同时.还可根据机器设备正常工作和不正常工作时,振动加速度信号的功率谱的差别,查找不正常工作时,功率谱图中额外谱峰产生的原因以及排除故障的方法。
自功率谱密度函数定义及其物理意义假如)(t x 是零均值的随机过程,即0=x μ(如果原随机过程是非零均值的,可以进行适当处理使其均值为零)又假设)(t x 中没有周期分量,那么当∞→τ,0)(→τx R 。
功率谱分析例要点在进行功率谱分析时,有几个重要的例要点需要注意:1.信号处理前的准备工作:在进行功率谱分析之前,我们需要对信号进行一些预处理,以确保分析的准确性。
这包括去除潜在的噪声、滤波和信号采样等步骤。
这些预处理方法的选择取决于应用的具体要求和信号的特性。
2.快速傅里叶变换(FFT):FFT是计算功率谱的常用方法,它可以在计算上更高效地将信号从时域转换为频域。
FFT通过将信号拆分成不同频率的正弦和余弦函数来实现这种转换。
FFT算法的使用可以大大加快功率谱分析的速度。
3.窗函数的选择:在进行FFT之前,通常需要将信号分成不同的时间窗口。
窗口函数有助于减少谱泄漏(spectral leakage)效应,即当一个窗口函数不匹配信号的特征时,信号能量会泄漏到其他频率上。
常用的窗口函数有矩形窗、汉宁窗、汉明窗等。
选择合适的窗口函数取决于信号的特性以及应用的要求。
4.相对功率谱与绝对功率谱:相对功率谱是指功率谱除以总功率的比例。
它表示不同频率分量的能量在信号中所占的比例。
相对功率谱可以帮助我们了解信号的频率分布情况。
而绝对功率谱表示不同频率分量的能量或功率的绝对值。
绝对功率谱对于分析信号的绝对强度和功率分布很有用。
5.峰值频率和带宽:在功率谱分析中,我们可以通过查找功率谱图中的峰值频率来确定信号中的主要频率分量。
峰值频率表示信号中能量最强的频率。
带宽则表示主要频率分量的频率范围。
对于宽频信号,带宽可能会很大,而对于窄频信号,带宽则较小。
6.平滑功率谱:平滑功率谱可以帮助我们去除谱图中的不稳定和噪声。
平滑功率谱使用低通滤波器对功率谱进行滤波,从而减少高频分量的影响。
平滑功率谱可以提供一个更稳定的频域表示,并突出主要频率分量。
7.谱密度与积分功率谱:谱密度是功率谱密度函数的积分,表示信号的总功率。
通过计算谱密度,我们可以获得信号在整个频谱范围内的功率值。
谱密度是理解信号能量分布的关键指标。
总而言之,功率谱分析是一种重要的信号处理工具,它可以帮助我们理解信号的频率特性、能量分布以及峰值频率等。
单相电动机的电流波形和功率谱分析单相电动机是广泛应用于家庭和工业领域的电动机之一。
了解和分析单相电动机的电流波形和功率谱对于正确运行和故障排除至关重要。
在本文中,我们将深入探讨单相电动机的电流波形和功率谱分析。
一、单相电动机的工作原理单相电动机是利用交变电源生成的交流电流驱动的。
它主要由定子和转子组成。
当交流电源施加在定子上时,定子绕组会产生一个旋转磁场,与转子上的绕组相互作用,从而产生转矩使转子转动。
单相电动机通常具有起动电容器和运行电容器,起动电容器用于提供启动转矩,而运行电容器则用于提高效率和功率因数。
二、单相电动机的电流波形单相电动机的电流波形通常可以分为三个阶段:启动阶段、运行阶段和恢复阶段。
在启动阶段,电流的波形通常呈现出高峰值和大幅度的变化,这是由于起动电容的作用。
一旦电动机达到运行阶段,电流波形将会变得更加稳定,大致呈现类似正弦波的形状。
在恢复阶段,当电动机停止工作或受到额外负载时,电流波形将再次发生变化。
三、单相电动机的功率谱分析功率谱是指信号在频率域的分布情况。
对于单相电动机的功率谱分析,我们主要关注其基波、谐波和杂波等频率成分。
1. 基波:单相电动机的基波频率一般为电源频率,如50Hz或60Hz,它代表了电动机正常运行的主要频率成分。
基波幅值的大小反映了电动机的负载情况。
2. 谐波:谐波是频率是基波频率的整数倍的成分,它们是由于电动机的非线性负载或电源的不纯度引起的。
谐波会增加电动机的功率损耗,导致效率下降并可能引起其他问题,如噪音和振动。
3. 杂波:杂波是指频率不是基波频率的倍数的成分,它们可能是由于电动机运行过程中的故障、电源干扰或其他干扰引起的。
杂波的存在可能导致电动机运行不稳定或产生异常响声。
通过对单相电动机的功率谱分析,我们能够更好地了解电动机的性能和运行状况。
对于基波的监测可以帮助我们评估电动机的负载情况,而谐波和杂波的分析则有助于检测故障和干扰的存在。
四、电流波形和功率谱分析的应用电流波形和功率谱分析可用于以下方面:1. 故障检测:通过分析电流波形和功率谱,我们可以检测到电动机可能存在的故障,如绕组短路、轴承磨损或不平衡负载等。
功率谱密度分析在信号处理中的应用信号是随着时间变化的电压,电流,电磁波等物理量。
信号分析是从信号中获取有用信息的过程。
这种信号常常是含有噪声的,并且要从中提取出所需的信息。
由于信号需要先进行预处理,因此,信号处理是一个复杂的任务。
在信号处理领域,功率谱密度分析是一种常用的技术,被广泛应用于信号处理和系统分析中。
一、功率谱密度分析的基本概念功率谱密度分析的目标是确定一个信号在不同频率下的功率,这是一种分析信号的频域方法。
功率谱密度是指信号在一个频带内的功率的分布,单位是瓦特/赫兹(W/Hz)。
功率谱密度分析的输出结果一般呈现为功率谱密度图,它描述了信号的能量随着频率的变化而变化的情况。
功率谱密度的计算主要基于伯努利-欧拉定理,即将复变量表示为实部与虚部的和。
对于一个实值信号x(t),其傅里叶变换H(f)如下所示:H(f)=∫x(t)exp[-2πi f t]dt然后,对于信号x(t)和其复共轭x* (t),可以计算出它们的积:P(f)=x(t)×x*(t)其中,t 代表时间,f 代表频率。
对于连续时间信号,P(f) 被称为功率谱密度,表示频率 f 的功率。
对于离散时间信号,其内积被替换为求和,并且功率谱密度的单位变为瓦特/赫兹(W/Hz)。
二、功率谱密度分析的应用功率谱密度分析在信号处理中有着广泛的应用,下面我们主要介绍其在音频处理和图像处理中的应用:1. 音频处理中的功率谱密度分析音频信号的功率谱密度是指一段时间内声音量随着频率变化的标志。
在音频处理中,功率谱密度分析可以用于识别音频信号的特定频率成分,并清除噪声。
在使用数字信号处理算法对音频信号进行无噪声处理时,功率谱密度图经常被使用。
通过检测功率谱密度的凸起与波峰,可以识别音频信号的某些特定频率。
功率谱密度分析还可以用于滤波器设计。
具体地说,使用功率谱密度可以确定所需滤波器的特性,例如通带的大小、截止频率等,从而设计出能清除干扰和噪声的专用滤波器。
三、功率谱分析字体[大][中][小]周期信号的功率谱为其双边幅值频谱的平方|c n|2;非周期信号的功率谱为其幅值谱密度的平方|X(ω)|2=X(ω)X*(ω)。
随机信号属于时域无限信号,其频率、幅值和相位为随机变量。
因而,采用具有统计特性的功率谱估计进行谱分析(一)自功率谱密度及其估计各态历经随机信号的功率谱密度S x(ω)与自相关函数R x(τ)为傅里叶变换偶对,即为了方便,也可用在非负频率范围内(ω>0)定义的单边功率谱密度G x(ω)代替双边功率谱密度S x(ω),两者之间的关系为自功率谱估计可分为线性估计法与非线性估计法。
前者以快速变换为基础,应用较早,也称为经典谱分析法; 后者是与时序模型结合的一种新方法,又称为现代谱分析方法。
1. 周期图各态历经随机信号的均方值ψx2为信号能量的时域描述。
巴什瓦定理表明,信号能量的时域计算与频域计算相等,即由此定义自功率谱密度及其估计为:式中表12-45 典型信号的自相关、频谱、概率密度(续)X(ω)为测试数据x(t)的傅里叶变换,X(k)为N个数据x(n)的离散傅里叶变换,由FFT 直接求出。
由于X(k)具有周期函数的性质,所以称由此获得的自功率谱估计为周期图。
自相关估计x′(r)的快速傅里叶变换可作为自功率谱估计的另一计算公式以上两种估计都是自功率谱S x(ω)的有偏估计,只是偏差大小不同。
两种估计在时域对数据或对自相关估计进行截断,相当于加窗处理,致使谱估计成为真实功率谱(或称为真功率谱)与窗谱W(ω)的卷积,即Ŝx(ω)=S x(ω)*W(ω)窗谱旁瓣的泄漏效应和卷积的作用使真功率谱的尖峰数值变化,邻近点的数值变大,造成谱估计的模糊与失真以上两种估计的方差较大; 相距2π/N的各点估计值互不相关,故数据点数N越大,这些点的估计值的随机起伏越严重。
为改善谱估计的估计质量,在增大数据点数的同时,采用平均化处理和窗处理方法减小谱估计的方差。
2. 修正周期图平均化和加窗处理可使改善的周期图估计质量提高,而且还可保留周期图便于应用FFT 计算速度快的优点。
(1) 平均周期图把N点的长数据序列分成k段,每段数据点数为M=N/k; 求得各段周期图Ŝx i(ω)后再用平均法求得平均周期图Ŝx av(ω)。
平均处理使谱估计的方差减小为当N一定时,段数K大则各段数据点数M小,谱估计的偏差大、方差小、谱平滑但频率分辨率低;若K小、M大,则偏差小、方差大。
(2) 加窗谱估计选择适当的窗函数ω(r),对自相关估计序列x′(r)作加窗处理,然后求谱估计加窗谱估计平滑、方差小,故又称为平滑周期图。
常用窗函数有矩形窗、三角窗、汉宁窗、海明窗和高斯窗等。
图12-7表示了分别用矩形窗、三角窗和汉宁窗分析同一数据所得结果。
图a是被分析信号的真谱; 图b和c是用两种时宽的矩形窗分析的结果。
可以看出,时域窗口窄,则频域分辨率低,两相邻谱线无法区分。
图d和e表明三角窗和汉宁窗频率分辨率低,但旁瓣衰减快,谱的分布区域窄而边沿清晰。
显然,应按信号的性质和处理要求适当选择窗函数。
图12-7 窗函数效果比较(3) 平滑平均周期图平滑谱估计计算要求所选用的窗函数保证求出的谱估计非负,但有的窗函数不满足此要求。
如果先将数据x(n)分段,求出分段数据的加窗谱估计,然后再将各分段谱估计作平均化处理,即可满足谱估计非负的要求,又可减小估计偏差和估计方差,使谱估计质量提高。
按此法计算的谱估计称为平滑平均周期图。
N个数据x(n)分成k段,k=N/M。
按照快速傅里叶变换的要求,将段内数据补零使M=2γ(γ为正整数)。
求取各段的加窗谱估计然后用平均法求平滑平均周期图由于P x i(k)计算用的数据经过加窗修正,为使谱估计为真谱的渐近无偏估计,在计算平滑平均谱估计Ŝx(k)时必须作相应的反修正,上式中U就是反修正因子——归一化因子(二)互功率谱密度对于零均值各态历经随机信号x(n)、y(n)的时域和频域描述,与自相关估计、自功率谱估计类似。
但是,由于互相关函数并非偶函数,所以互相关函数的估计量分别给作互功率谱密度简称互谱,与互相关函数为傅里叶变换关系。
互谱估计为互相关估计的离散傅里叶变换,可通快速傅里叶变换FFT求得。
当有限长数据x(n)、y(n)的数据点数为N时,则估计方法和估计质量与自功率谱估计类似。
在互谱分析中,常引用相干函数γy(ω)分析两个各态历经随机信号x(t)、y(t)之间以及它们与系统特性的相互关系。
若在某频率上,γxy2(ω)=0,则表示x(t)与y(t)在此频率上不相干。
如果x(t)与y(t) 是统计独立的,则在所有频率上都有γxy2(ω)=0。
若在所有频率上都有γxy2(ω)=1,则x(t) 与y(t)完全相干。
当γxy2(ω)<1时,则说明:1)测量中有不通过系统的外来干扰; 2)联系x(t)和y(t)的系统不是线性系统; 3)y(t)是x(t)和其他输入的综合输出。
相干函数的估计为其中的自功率谱和互谱的估计都应当是经过频率平均或总体平均后的估计,而且还应保证x(t)、y(t)的测试数据是同时测量、具有相同分析参数(如记录长度T等),以免分析结果出现错误。
(三) 极大熵谱估计伯格(J.P.Burg)提出的极大熵谱分析方法(Maximum Entropy Method,简称MEM)是一种新的功率谱估计方法。
伯格在最大熵谱估计准则的提出和具体算法上有所创新。
由此演变出来的算法有多种,被统称为“现代谱分析”。
1. 极大熵准则传统谱估计方法实际上都是把无限长序列加窗截断后,由有限长序列获得功率谱估计。
不论是对原始数据加窗还是对自相关函数加窗处理,其目的都是减小谱估计的方差、提高频率分辨率。
然而窗处理不可避免地产生频域的“泄漏”,使功率谱失真。
尽管在窗函数和处理方法上进行了许多研究,使得以周期图为基础的谱估计方法广泛应用,但谱估计的频率分辨率并不能令人满意。
此外,传统谱估计方法通过增大数据点数来获得较高的谱估计精度。
这样,不仅数据处理工作量大,而且对短记录数据或缓变信号等显然是无能为力的。
在谱估计计算中,对原始数据或自相关函数加窗处理,假设窗外的数据为零,而且还对窗内部分进行某些修正。
这些人为措施增加了确定性因素,使原来具有的不确定性减少,在一定程度上歪曲了观测得到的信息。
极大熵谱分析方法在谱估计计算中不作加窗处理,而是采用适当的方法把由有限长数据求得的自相关估计外推。
外推的原则是使相应的数据序列在外推点上取值的可能性具有最大的不确定性,亦即不对结果人为地增加任何附加信息。
在数理统计学中,“熵”表征了各种随机试验的统计特性,是随机总体的平均不确定性的量度。
极大熵谱分析法把熵的概念引入谱分析。
上述外推原则就是要求在使随机过程的熵达到最大的条件下,确定未知的自相关函数值,外推原则的最大熵描述就是谱估计的极大熵准则。
2. 极大熵谱估计的基本原理随机试验a有有限个不相容的结果A i,相应的概率为P(A i),且满足则随机试验的熵H(a)为其不确定性的量度。
对于连续随机变量,则熵为式中的对数可以取10或取e为底,在比较熵的大小时并没有影响正态随机变量x的概率密度为则熵为高斯随机过程的样本函数x(n),其数据点数为N,则其熵为式中|R|为N×N阶自相关矩阵[R]的行列式,也可用det[R]表示。
[R]为托布列兹矩阵。
方差为σ的自噪声的熵为Hε=-E[ln(p(ε1)p(ε2)…p(εN))]=Nlnσε由于H随N增长将要发散,故用熵率h代替,熵率定义为由于功率谱密度与自相关函数之间为傅里叶变换关系,则自相关矩阵[R]的行列式|R|的极限与功率谱密度S x(f)的关系是:式中f N=2f c,为Nyquist采样频率;f c——信号的频限,即带宽,由此可以得到熵率为:当由有限长随机信号测试数据计算自相关估计并外推时,为使每步估计的熵达到最大,则应有(k≥m+1,m——数据点数)由此得到:式中,功率谱密度S x(f)应受到熵率导数为零的条件约束,1/S x(f)可用截短的傅里叶级数表示成为2m+1项有限级数之和,于是有和从而求得极大熵谱式中△t——采样时间间隔;a k——自回归模型参数,a k*与a k共轭;P m——外推误差的方差MEM谱与参数a k和P m有关。
典型的算法有Burg、Levinson、Akaike和Marple 算法。
3. 极大熵谱的特点与经典谱估计方法相比,极大熵谱估计方法克服了窗处理带来的一系列缺陷。
而且由于是连续谱,谱线光滑、谱峰陡峭,频率分辨力远高于经典谱估计; 而其谱分辨力与样本函数长度无关,特别适用于短数据样本、缓变过程的谱估计。
图12-8表示短数据分析时FFT谱与MEM谱的对比。
被分析信号是1Hz正弦波上叠加10%的白噪声。
图a样本长度为1s时,两种谱均显示出主峰频率为1Hz;图b表示样本长度为0.57s时,MEM谱主峰频率清晰准确,而FFT谱却无法识别。
图12-9表示缓变信号的MEM谱与FFT谱的比较。
被分析的信号是0.6Hz的正弦波,由图可见,MEM谱较为精确。
图12-8 短数据的MEM谱与FFT谱比较图12-9 缓变数据的MEM谱与FFT谱谱分析与谱估计在生产实践和科学研究中获得了日益广泛的应用。
例如,在声纳系统中,为了寻找海洋水面舰艇或潜艇,需要对噪声信号进行谱分析,以提取有用信息,判断舰艇速度、方位、大小等; 对各种机电产品主体或部件作实际运行的谱分析,可检验设计的效果并提出改进设计的依据,或者进行在线监测与故障诊断,以便及时排除潜在的故障因素,保证安全运行。