首页
学习
活动
专区
工具
TVP
发布
精选内容/技术社群/优惠产品,尽在小程序
立即前往
  • 您找到你想要的搜索结果了吗?
    是的
    没有找到

    随机振动 matlab,Matlab内建psd函数在工程随机振动谱分析中的修正方法「建议收藏」

    随机信号的功率谱分析是一种广泛使用的信号处理方法,能够辨识随机信号能量在频率域的分布,同时也是解决多种工程随机振动问题的主要途径之一.Matlab作为大型数学分析软件,得到了广泛应用,目前已推出7.x的版本.Matlab内建了功能强大的信号处理工具箱.psd函数是Matlab信号处理工具箱中自功率谱分析的主要内建函数.Matlab在其帮助文件中阐述psd函数时均将输出结果直接称为powerspectrumdensity,也即我们通常所定义的自功率谱.实际上经分析发现,工程随机振动中功率谱标准定义[1]与Matlab中psd函数算法有所区别,这一点Matlab的帮助文档没有给出清晰解释.因此在使用者如没有详细研究psd函数源程序就直接使用,极易导致概念混淆,得出错误的谱估计.本文详细对比了工程随机振动理论的功率谱定义与Matlab中psd函数计算功率谱的区别,并提出用修正的psd函数计算功率谱的方法,并以一组脉动风压作为随机信号,分别采用原始的psd函数与修正后的psd函数分别对其进行功率谱分析,对比了两者结果的差异,证实了本文提出的修正方法的有效性.1随机振动相关理论1.1傅立叶变换求功率谱理论上,平稳随机过程的自功率谱密度定义为其自相关函数的傅立叶变换:Sxx()=12p+-Rxx(t)eitdt(1)其中,S(xx)()为随机信号x(t)的自功率谱密度,Rxx(t)为x(t)的自相关函数.工程随机振动中的随机过程一般都是平稳各态历经的,且采样信号样本长度是有限的,因此在实用上我们采用更为有效的计算功率谱的方法,即由时域信号x(t)构造一个截尾函数,如式(2)所示:xT(t)=x(t),0tT0,其他(2)其中,t为采样时刻,T为采样时长,x(t)为t时刻的时域信号值.由于xT(t)为有限长,故其傅立叶变换A(f,T)以及对应的逆变换存在,分别如式(3)、(4)所示:A(f,T)=+-xT(t)e-i2pftdt(3)xT(t)=+-A(f,T)ei2pftdt(4)由于所考虑过程是各态历经的,可以证明:Sxx(f)=limT1TA(f,T)2(5)在实际应用中,式(5)是作功率谱计算的常用方法.1.2功率谱分析中的加窗和平滑处理在工程实际中,为了降低工程随机信号的误差,一般对谱估计需要进行平滑处理.具体做法为:将时域信号{x(t)}分为n段:{x1(t)},{x2(t)},…,{xn-1(t)},{xn(t)},对每段按照式(5)求功率谱Sxixi(f),原样本的功率谱可由式(6)求得:Sxx(f)=1nni=1Sxixi(f)(6)如取一样本点为20480的样本进行分析,将样本分割为20段进行分析,每段样本点数为1024.将每段1024个样本点按照式(5)的方法分别计算功率谱后求平均,即可得到经过平滑处理的原样本的功率谱,这样计算出的平滑谱误差比直接计算要降低很多.另一方面,由于实际工程中随机信号的采样长度是有限的,即采样信号相当于原始信号的截断,即相当于用高度为1,长度为T的矩形时间窗函数乘以原信号,导致窗外信息完全丢失,引起信息损失.时域的这种信号损失将会导致频域内增加一些附加频率分量,给傅立叶变换带来泄漏误差.构造一些特殊的窗函数进行信号加窗处理可以弥补这种误差,即构造特殊的窗函数{u(t)},用{u(t)}去乘以原数据,对{x(t)u(t)}作傅立叶变换可以减少泄漏:Aw(f,T)=+-u(t)xT(t)e-i2pftdt(7)其中,Aw(f,T)为加窗后的傅立叶变换.u(t)xT(t)实际上是对数据进行不等加权修改其结果会使计算出

    01

    基于matlab的语音信号频谱分析_声音信号的数字化过程

    随着软硬件技术的发展,仪器的智能化与虚拟化已成为未来实验室及研究机构的发展方向[1]。虚拟仪器技术的优势在于可由用户定义自己的专用仪器系统,且功能灵活,很容易构建,所以应用面极为广泛。基于计算机软硬件平台的虚拟仪器可代替传统的测量仪器,如示波器、逻辑分析仪、信号发生器、频谱分析仪等[2]。从发展史看,电子测量仪器经历了由模拟仪器、智能仪器到虚拟仪器,由于计算机性能的飞速发展,已把传统仪器远远抛到后面,并给虚拟仪器生产厂家不断带来连锅端的技术更新速率。目前已经有许多较成熟的频谱分析软件,如SpectraLAB、RSAVu、dBFA等。

    01

    声音处理之-梅尔频率倒谱系数(MFCC)

    在语音识别(SpeechRecognition)和话者识别(SpeakerRecognition)方面,最常用到的语音特征就是梅尔倒谱系数(Mel-scaleFrequency Cepstral Coefficients,简称MFCC)。根据人耳听觉机理的研究发现,人耳对不同频率的声波有不同的听觉敏感度。从200Hz到5000Hz的语音信号对语音的清晰度影响对大。两个响度不等的声音作用于人耳时,则响度较高的频率成分的存在会影响到对响度较低的频率成分的感受,使其变得不易察觉,这种现象称为掩蔽效应。由于频率较低的声音在内耳蜗基底膜上行波传递的距离大于频率较高的声音,故一般来说,低音容易掩蔽高音,而高音掩蔽低音较困难。在低频处的声音掩蔽的临界带宽较高频要小。所以,人们从低频到高频这一段频带内按临界带宽的大小由密到疏安排一组带通滤波器,对输入信号进行滤波。将每个带通滤波器输出的信号能量作为信号的基本特征,对此特征经过进一步处理后就可以作为语音的输入特征。由于这种特征不依赖于信号的性质,对输入信号不做任何的假设和限制,又利用了听觉模型的研究成果。因此,这种参数比基于声道模型的LPCC相比具有更好的鲁邦性,更符合人耳的听觉特性,而且当信噪比降低时仍然具有较好的识别性能。

    02
    领券