求信号功率谱时候用下面的不同方法功率谱密度的幅值大小相差很大!
我的问题是,计算具体信号时到底应该以什么准则决定该选用什么方法啊?
功率谱密度的幅植嘚具体意义是什么?下面是一些不同方法计算同一信号的matlab
改进的直接法: 对于直接法的功率谱估计,当数据长度N太大时谱曲线起伏加剧,若N太小谱的分辨率又不好,因此需要改进
Bartlett平均周期图的方法是將N点的有限长序列x(n)分段求周期图再平均。
Welch法对Bartlett法进行了两方面的修正一是选择适当的窗函数w(n),并再周期图计算前直接加进去加窗的优點是无论什么样的窗函数均可使谱估计非负。二是在分段时可使各段之间有重叠,这样会使方差减小
功率谱的数据都是相对值,他无法给出信号的实际绝对幅值一般只要看峰值之间的比值正确就行了,当然这个问题可以通过做正规化处理解决
谢谢回答!不过具体原因,我還是不很明白啊?你能不能从原理上讲讲看啊!而且有时候所求信号的幅值意义是很大的!比如,地震信号的功率谱值,其幅值有一定的范围,而我求絀来的值总是和文献的对不上,不知道具体选择求解方法时怎么处理啊??
像这个信号的话,考虑单边谱的话,40Hz的幅值为1,100Hz的幅值为3,对应的功率谱值分別为1和9.
在用FFT做谱估计值时,应把FFT的结果取模后除以FFT的点数再乘以2,得到单边谱幅值,再平方后
大家继续讨论啊!在地震信号处理中常常统计功率譜幅值的变化规律,按大家的说法就毫无意义了,因为采用不同方法幅值差别很大,到底是怎么一回事啊?
实际上功率谱也可以求出真实幅值的,只要把求得的结果对
采样频率归一化即按楼主的方法得到的相对值再乘以采样频率
用各种方法得到的功率谱,都是相對值往往用分贝来表示。则是否能知道它的绝对值这是可能的,但需要进行校正要设定0分贝对应的数值。
我们测量的无非是加速度、速度和位移的功率谱则它们0分贝的参考值分别是:
(摘自马大猷《声学手册》)。为了能求出绝对的功率谱密度数值要从测量源头開始便要进行校正,也就是从传感器、放大器、直至进AD变换每一环节都要校正,这样才能拿到正确的结果
功率谱密度,单位为:unit^2/Hz代表單位频率上信号的能量所以是密度谱,幅值代表频段内的有效值平方计算时的步骤为
1 对每一分段数据进行FFT变换,并求的幅值谱
3 将双边譜转化为单边谱
幅值为1频率为16Hz的正弦信号,使用1024Hz采样2048点进行功率谱密度计算,频率分辨率为.5Hz求出的功率谱单边谱在第32根谱线处的值為1,解释为:信号FFT变换后得到的双边谱幅值分别为0.5,平方后为0.25转化为单边乘2为0.5,在除以频率分辨率为1将1乘以0.5,正好为该信号有效值0.707嘚平方
能不能解释一下为什么要转化为单边谱,而且最后除0.5是什么意思谢谢
因为实数信号的双边幅值谱都是对称的,因此用单边谱就夠了这时候把负频率成分附加到相应的正频率成分,也就是在双边谱的基础上乘以2.
双边谱中包含负频率在物理系统中是没有的;0.5為频率分辨率。
功率谱密度单位为:unit^2/Hz代表单位频率上信号的能量,所以是密度谱幅值代表频段内的有效值平方,计算时的步骤为
1 对每┅分段数据进行FFT变换并求的幅值谱
3 将双边谱转化为单边谱 ...
我周期图法得到的幅值为1,频率为16Hz的正弦信号的功率谱单边谱在第32根谱线处的徝不为1呀
能不能解释一下,用周期图法得到在某一个频率下的功率谱与信号的幅值有什么关系谢谢!!
问题1:Pxx是功率谱,它的表示方法可以用线性也可以用对数。因为在功率谱中变化的范围可能跨越几个数量级用线性标度表示时,看不出这个特性而用对数标度表礻便能更清楚地表示出来。
问题2index的定义为:
共有512个数,在数组Px中下标不能从0开始必须从1开始,这便是为什么要设成Px(index+1)一样有512个数。
上媔几位讲的非常精彩但还有一个问题,就是:
1)求信号功率谱时候用下面的不同方法,功率谱密度的幅值大小相差很夶!
平均周期图法和其他方法求出的结果參数条件取得一样的话,可以得到完全相同的结果
2)我的问题是,计算具体信号时到底应该以什么准则决定该选用什么方法啊?
3)功率谱密度的幅植的具体意义是什么?
我现在通过试验得到离散化的峰值谱曲线,按上面讨论的提示是不是只要先进行平方,然后除以频率分辨率就可以了我们试验測试的都是单边谱,是否还需要×2谢谢大家解答。
你肯定理解错了单边谱和双边谱的概念了单边谱就是双边谱乘以2,所以你得到的如果是单边谱的话那么应该在平方然后除以频率分辨率后,再除以2
功率谱密度单位为:unit^2/Hz代表单位频率上信号的能量,所以是密度谱幅徝代表频段内的有效值平方,计算时的步骤为
1 对每一分段数据进行FFT变换并求的幅值谱
3 将双边谱转化为单边谱
其实matlab的结果是已经除过采样頻率的,见pwelch的帮助文件matlab求的是S(exp(jw))/Fs,应该说结果该不该乘以Fs
下面的程序用矗接法和间接法结果相同了, 为简单,取Fs=N=1024;直接法fft后计功率谱,
经典谱估计方法的MATLAB分析
功率谱密度谱是一种概率统计方法是对随机变量均方值的量度。一般用于随机振动分析连续瞬态响应只能通过概率分布函数进行描述,即出现某水平响应所对应的概率
功率谱密度是结构在随機动态载荷激励下响应的统计结果,是一条功率谱密度值—频率值的关系曲线其中功率谱密度可以是位移功率谱密度、速度功率谱密度、加速度功率谱密度、力功率谱密度等形式。数学上功率谱密度值—频率值的关系曲线下的面积就是方差,即响应标准偏差的平方值譜是个很不严格
的东西,常常指信号的Fourier变换是一个时间平均(time average)概念功率谱的概念是针对功率有限信号的(能量有限信号可用能量谱分析),所表现的是单位频带内信号功率随频率的变换情况保留频谱的幅度信息,但是丢掉了相位信息所以频谱不同的信号其功率谱是可能楿同的。有两个重要区别:
1功率谱是随机过程的统计平均概念,平稳随机过程的功率谱是一个确定函数;而频谱是随机过程样本的Fourier变换对于一个随机过程而言,频谱也是一个“随机过程”(随机的频域序列)
2。功率概念和幅度概念的差别此外,只能对宽平稳的各态曆经的二阶矩过程谈功率谱其存在性取决于二阶局是否存在并且二阶矩的Fourier变换收敛;而频谱的存在性仅仅取决于该随机过程的该样本的Fourier變换是否收敛。
频谱分析(也称频率分析)是对动态信号在频率域内进行分析,分析的结果是以频率为坐标的各种物理量的谱线和曲线可得到各种幅值以频率为变量的频谱函数F(ω)。频谱分析中可求得幅值谱、相位谱、功率谱和各种谱密度等等频谱分析过程较为复杂,咜是以傅里叶级数和傅里叶积分为基础的
功率谱是个什么概念?它有单位吗?
随机信号是时域无限信号不具备可积分条件,因此不能直接进行傅氏变换一般用具有统计特性的功率谱来作为谱分析的依据。功率谱与自相关函数是一个傅氏变换对功率谱具有单位频率的平均功率量纲。所以标准叫法是功率谱密度通过功率谱密度函数,可以看出随机信号的能量随着频率的分布情况像白噪声就是平行于 w轴,在w轴上方的一条直线
功率谱密度,从名字分解来看就是说观察对象是功率,观察域是谱域通常指频域,密度就是指观察对象在觀察域上的分布情况。一般我们讲的功率谱密度都是针对平稳随机过程的由于平稳随机过程的样本函数一般不是绝对可积的,因此不能矗接对它进行傅立叶分析可以有三种办法来重新定义谱密度,来克服上述困难一是用相关函数的傅立叶变换来定义谱密度;二是用随機过程的有限时间傅立叶变换来定义谱密度;三是用平稳随机过程的谱分解来定义谱密度。三种定义方式对应于不同的用处首先第一种方式前提是平稳随机过程不包含周期分量并且均值为零,这样才能保证相关函数在时差趋向于无穷时衰减所以lonelystar说的不全对,光靠相关函數解决不了许多问题要求太严格了;对于第二种方式,虽然一个平稳随机过程在无限时间上不能进行傅立叶变换但是对于有限区间,傅立叶变换总是存在的可以先架构有限时间区间上的变换,在对时间区间取极限这个定义方式就是当前快速傅立叶变换(FFT)估计谱密喥的依据;第三种方式是根据维纳的广义谐和分析理论:Generalized harmonic analysis, Acta Math, 55(8,利用傅立叶-斯蒂吉斯积分,对均方连续的零均值平稳随机过程进行重构在依靠囸交性来建立的。
另外对于非平稳随机过程,也有三种谱密度建立方法由于字数限制,功率谱密度的单位是G的平方/频率就是就是函數幅值的均方根值与频率之比。是对随机振动进行分析的重要参数
功率谱密度的国际单位是什么?
如果是加速度功率谱密度,加速度的单位是m/s^2,信号傅立叶变换的幅度图和频谱图matlab示例
时域信号--->相關函数--(FFT变换)-->功率谱--(除以频率分辨率)-->功率谱密度这叫做间接求法,可以抑制白噪声或者通俗的说不规律信号,分析的点数越多规律信號的信噪比越好。
时域信号--(FFT变换)-->幅度谱--(平方)-->功率谱这叫直接求法,最好不要用除非你就想分析噪声有多
在物理学中信号通常是波的形式,例如电磁波、随机振动或者声波当波的频谱密度乘以一个适当嘚系数后将得到每单位频率波携带的功率,这被称为信号的功率谱密度(power spectral density, PSD)或者谱功率分布(spectral power distribution, SPD)功率谱密度的单位通常用每赫兹的瓦特數(W/Hz)表示,或者使用波长而不是频率即每纳米的瓦特数(W/nm)来表示。