求一个音频的matlab求功率中的功率谱程序,只能用傅里叶变换吗?

小木虫,学术科研互动社区,为中国學术科研免费提供动力

违规贴举报删除请发送邮件至:emuch2018@


}
  • 专业: 信号理论与信号处理

功率谱密度函数对应随机序列相关函数的傅里叶变换一个随机序列的随机性还需要某时间点上随机变量的概率密度函数描述。

    假定你想产生一個功率谱为特性形式如s(f),且满足高斯分布的平稳随机序列可以先用ra ...


非常感谢解答,但生成这个数字滤波器时采样频率应该怎么设定啊……
}

版权声明:本文为博主原创文章未经博主允许不得转载。 /sinat_/article/details/

对于满足狄利克雷条件的周期信号可以进行傅里叶展开对于周期函数f(t),其周期为T1角频率为w1 1/T1,其傅里叶级数展开为:

上式的右端就是所谓的傅里叶级数周期信号按照傅里叶级数展开,可以分解成为直流分量a0和许多正余弦分量an,bn而这些正余弦的汾量必定是基频f1(1/T,T为周期函数的周期)的整数倍)。

将(1)式相同的频率项合并可以写成:

cndnnw1 的关系绘制成下图的曲线,得到的就是周期信号的幅度频谱其中每一条线代表某一频率分量的的幅度。

根据欧拉公式可以将傅里叶级数写成指数形式

而对于非周期信号可以视為周期无限大的周期信号,但当T1无限大时w1 = 2π/T1间隔无限小,因此非周期信号的谱线就由离散变为连续并且此时谱线长度Fn趋于0,则按照周期信号计算的频谱化为乌有但是从物理上来说信号存在,那么必定有相应的能量存在因此为了解决这个问题引入了频谱密度函数从傅裏叶级数导出傅里叶变换。具体推导过程本文就不详讲直接给出傅里叶变换的公式:

傅里叶变换又分离散和非离散,这里给出一个总结

茬上文讲到一个信号可以分解成无数个周期信号的叠加,每个信号都有自己的频率分量描述这种频率分量的图叫做频谱图。频谱图分為两部分一部分是幅度频谱,另一部分是相位谱语音中不怎么关注相位,因此下面会多讲一下幅度谱任意信号经过傅里叶变换(FFT)鈳到的是一个复数序列,每一项记为a+ib以下面代码为例

运行后我们可以看到,其中各个变量的类型与数目如上图所示由于傅里叶变换后嘚点的个数与时域信号采样点个数相同(signal数组数目等于Y数组的数目),但是傅里叶变换是堆成的因此只用显示其中的一半即可(P1一半会哆取一点,这样不用考虑采样点的奇偶情况)时域波形如下图所示

傅里叶变换得到复数的模即为频谱的幅值:


如下图所示,第一张图为单邊频谱第二张图为双边频谱



有时候也会用对数计算后的结果来表示频谱,此时单位为dB


信号的功率谱为信号频谱的平方根据帕斯维尔定悝,实信号的能量等于平均功率谱即


其中为功率谱或者功率谱密度函数

}

我要回帖

更多关于 matlab求功率 的文章

更多推荐

版权声明:文章内容来源于网络,版权归原作者所有,如有侵权请点击这里与我们联系,我们将及时删除。

点击添加站长微信