首页
学习
活动
专区
工具
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函数加窗

    数字信号处理的主要数学工具是傅里叶变换.而傅里叶变换研究的是整个时间域和频率域的关系。不过,当运用计算机实现工程测试信号处理时,不可能对无限长的信号进行测量和运算,而是取其有限的时间片段进行分析。做法是从信号中截取一个时间片段,然后用截取的信号时间片段进行周期延拓处理,得到虚拟的无限长的信号,然后就可以对信号进行傅里叶变换、相关分析等数学处理。无限长的信号被截断以后,其频谱发生了畸变,原来集中在f(0)处的能量被分散到两个较宽的频带中去了(这种现象称之为频谱能量泄漏)。 为了减少频谱能量泄漏,可采用不同的截取函数对信号进行截断,截断函数称为窗函数,简称为窗。 窗函数就是时域有限宽的信号。 泄漏与窗函数频谱的两侧旁瓣有关,如果两侧瓣的高度趋于零,而使能量相对集中在主瓣,就可以较为接近于真实的频谱,为此,在时间域中可采用不同的窗函数来截断信号。

    02

    基于MATLAB的语音信号处理

    摘要:语音信号处理是目前发展最为迅速的信息科学研究领域中的一个,是目前极为活跃和热门的研究领域,其研究成果具有重要的学术及应用价值。语音信号处理的研究,对于机器语言、语音识别、语音合成等领域都具有很大的意义。MATLAB软件以其强大的运算能力可以很好的完成对语音信号的处理。通过MATLAB可以对数字化的语音信号进行时频域分析,方便地展现语音信号的时域及频域曲线,并且根据语音的特性对语音进行分析。本文主要研究了基于MATLAB软件对语音信号进行的一系列特性分析及处理,帮助我们更好地发展语音编码、语音识别、语音合成等技术。本文通过应用MATLAB对语音信号进行处理仿真,包括短时能量分析、短时自相关分析等特性分析,以及语音合成等。

    01

    语音信号滤波去噪——使用FLATTOPWIN设计的FIR滤波器

    摘 要 本课程设计主要内容是设计利用窗口设计法选择FLATTOPWIN窗设计一个FIR滤波器,对一段含噪语音信号进行滤波去噪处理并根据滤波前后的波形和频谱分析滤波性能。本课程设计仿真平台为MATLAB7.0,开发工具是M语言编程,通过课程设计了解FIR滤波器设计的原理和步骤,掌握用MATLAB语言设计滤波器的方法,了解FLATTOPWIN对FIR滤波器的设计及编程方法。首先利用windows自带的录音机录制一段语音信号,加入一单频噪声,对信号进行频谱分析以确定所加噪声频率,设计滤波器进行滤波去噪处理,比较滤波前后的波形和频谱并进行分析。由分析结果可知,滤波 后的语音信号与原始信号基本一致,即设计的FIR滤波器能够去除信号中所加单频噪声,达到了设计目的。 关键词 滤波去噪;FIR滤波器;FLATTOPWIN窗;MATLAB

    04

    傅立叶分析和小波分析之间的关系? (通俗讲解)

    从傅里叶变换到小波变换,并不是一个完全抽象的东西,完全可以讲得很形象。小波变换有着明确的物理意义,如果我们从它的提出时所面对的问题看起,可以整理出非常清晰的思路。 下面我就按照傅里叶-->短时傅里叶变换-->小波变换的顺序,讲一下为什么会出现小波这个东西、小波究竟是怎样的思路。(反正题主要求的是通俗形象,没说简短,希望不会太长不看。。) 一、傅里叶变换 关于傅里叶变换的基本概念在此我就不再赘述了,默认大家现在正处在理解了傅里叶但还没理解小波的道路上。(在第三节小波变换的地方我会再形象地讲一下傅里叶变换)

    09

    最强大脑第二场战平听音神童!百度大脑小度声纹识别技术算法解析

    日前,继在江苏卫视《最强大脑》第四季“人机大战”首轮任务跨年龄人脸识别竞赛中击败人类顶级选手后,在上周五晚上,百度的小度机器人再次在声纹识别任务上迎战名人堂选手——11岁的“听音神童”孙亦廷,双方最终以1:1打成平手。被称为“鬼才之眼”的水哥(王昱珩)宣布再度出山,将在下周的第三轮比赛中与“小度”在图像识别方面一决高下。 本轮题目规则为:从“千里眼”到“顺风耳”,节目组将第二场比赛范围划定在“听”的领域,策划出高难度选题《不能说的秘密》,由周杰伦化身出题人,从21位性别相同、年龄相仿、声线极为相似的专业合

    06

    音视频开发之旅(66) - 音频变速不变调的原理

    音频的原始pcm数据是由 采样率、采样通道数以及位宽而定。常见的音频采样率是44100HZ,即一秒内采样44100次,采样通道数 一般为2, 代表双声道,而位宽一般是16bit 即2个字节。 通过改变采样率进行音频的变速,比如音视频播放器中的 2 倍速,0.5 倍速播放。如果想要实现音频的2.0倍速播放,只需要每隔一个样本点丢一个点,即采样率降低一半。如果想要实现0.5倍速播放,只需要每隔一个样本点插入一个值为0的样本点。就可以了,理想很丰满,但是如果仅仅这样做,带来的不止是速度的变化,声音的音调也发生变化了,比如 周杰伦的声音变成了萝莉音,这是我们不期望的。

    02
    领券