matlab 谱分析函数,科学网—经典谱分析（Power Spectrum Analysis）

频谱分析(或功率谱分析)大家可能都不陌生，然而细究起来，恐怕还是有很多模糊的地方。博主很早就知道并且学会使用各种数学软件来计算功率谱，可是长期以来总是知其然而不知其所以然，一些道理和概念总是含含糊糊过去就算了。前年年末有一段时间比较空闲，博主终于下决心理清这一堆乱麻，于是就有了本帖的学习笔记。笔记中试图弄清楚振幅谱、能谱和功率谱的区别，还有周期性信号、确定性非周期信号和随机信号的功率谱或能谱分析的区别，以及Fourier级数分析和离散Fourier变换的区别等等基本概念。从实际应用的角度出发，博主在笔记中用具体实例介绍了Matlab编写的功率谱分析的程序。

本帖附件完成于2008年1月6日，第一次修改于2009年7月6-7日，内容有：

标题：频谱分析

一、周期函数的Fourier级数展开

1、Fourier级数

2、如何对非周期函数进行Fourier级数展开

3、振幅谱、相位谱和功率谱

二、一般函数的Fourier积分

1、Fourier级数向Fourier积分的过渡

2、能谱和功率谱

3、自相关函数与功率谱及能谱的关系

三、平稳随机过程的Fourier积分

1、平稳随机过程的功率谱密度

2、Wiener-Khintchine公式

四、离散Fourier变换(DFT)

1、Fourier级数的离散化

2、Fourier积分的离散化

3、快速Fourier变化概述(FFT)

五、时间序列的能谱或功率谱分析

1、确定性时间序列的能谱或功率谱分析

2、随机时间序列的功率谱分析

六、直接法(或间接法)的缺陷及其改进方案

1、直接法和间接法总结

2、评价各种谱分析方法的标准

3、“平滑”和“泄漏”现象

4、大方差现象

5、改进方案(Bartlett、Welch、Blackman-Tukey方法)(未完待续)附录

参考文献

笔记中如有任何疏漏、错误，请不吝指教，博主对此表示万分地感谢。另外如果您觉得哪一部分有参考价值而引用到您的文章中，请指明引用的来源：http://www.sciencenet.cn/u/sanshiphy，博主不胜感激。

【2009-7-18】谱分析分为：“非参数化方法”和“参数化方法”，前者直接利用了Fourier变换的离散表达式，将数据代入表达式进行计算。将帖子的标题由“频谱分析”改为“经典谱分析”。此外，对笔记的内容再次做了修正，增加了周期图法的改进方案一节。周期图法(直接法)在计算随机数据时，无论数据量取的多大，其方差可能不会减少至0. 因此，为了减小方差，人们基于周期图法，提出了各种修改方案，这些方案都是以降低分辨率来达到减小方差的目的。

【2013-12-4】订正：p14 Tk 应是单边功率谱，原文双边功率谱是错误的——感谢肖擎曜同学的指正。

【2016-4-19】补充：吕瑞同学讨论了fft补零的问题

转载本文请联系原作者获取授权，同时请注明本文来自刘磊科学网博客。

链接地址：http://blog.sciencenet.cn/blog-200199-242357.html

上一篇：译稿：年度的、Milankovitch以及连续的温度变化之间的联系(Translation of "Links between annual, Milankov

下一篇：日志：《经典谱分析》一文有更新(New Additions to Power Spectral Analysis)

matlab 谱分析函数,科学网—经典谱分析（Power Spectrum Analysis） - 刘磊的博文相关推荐

matlab样条插值拟合,科学网—样条函数插值拟合 - 李继存的博文
样条函数插值拟合 2014–02–11 09:26:49 在拟合势能函数的时候, 除解析式外, 也可以利用样条函数进行拟合. 样条拟合与其插值正好相反: 已知函数在节点上的值求任意位置的值, 做插值; ...
matlab 二维矩形函数,科学网—利用MATLAB对非矩形域实现二维插值 - 张乐乐的博文...
>> load('x1.mat'); >> load('y1.mat') >> load('T.mat'); >> load('Lon.mat'); & ...
单纯性搜索算法 matlab函数,科学网—一种有效的最优化方法——Nelder-Mead单纯形直接搜索算法 - 王福昌的博文...
虽然MATLAB本身自带了fminsearch()函数,可以求解目标函数无梯度的最优化问题,但是感觉下面的程序在很多时候更好用,特别是自变量有边界和非线性约束的时候. 这里是John D'Errico ...
matlab meshgrid例子,科学网—matlab使用笔记-meshgrid函数的使用 - 骆佳良的博文
[X,Y]=meshgrid[x,y] 结果生成两个矩阵 X 与 Y. 矩阵X的每一行都与x相同.矩阵Y的每一列都与y相同. 矩阵X的行数等于矩阵y的行数,矩阵Y的列数等于矩阵x的列数. 例子: 输入 ...
matlab绘制加速度曲线,科学网-Matlab编写由加速度积分得到速度和位移函数-王德才的博文...
研究中经常会用到将加速度数据积分成速度和位移,对前几年编写的程序重新进行了改写.为了测试和验证peer2acc函数.dirFolder函数和acc2vd函数,编写了一个简单的测试程序 acc2vd函数 ...
1071svm函数 r语言,科学网—R中的svm - 吴锐的博文
svm理解: LSSVM: 最小二乘支持向量机(Least squares support vector Maehine,LSSVM)是SVM的一种变体,把问题转化成对一个线性方程求解,所需计算资源 ...
matlab画复变函数,科学网—复数复变函数的Matlab计算与绘图 - 周铁戈的博文
复数复变函数的Matlab计算与绘图周铁戈复数的表示存在两种表示方法,一种是代数式,一种是指数式,在Matlab中的方式如下: >> z=1+2i #代数式,1 ...
matlab聚类算法,科学网—matlab-聚类算法笔记 - 孙月芳的博文
MATLAB提供了两种方法进行聚类分析: 1.利用clusterdata 函数对数据样本进行一次聚类,这个方法简洁方便,其特点是使用范围较窄,不能由用户根据自身需要来设定参数,更改距离计算方法: 2. ...
matlab 水平投影,科学网—Matlab中如何将投影信息写入到shape文件中 - 朱永超的博文...
在Matlab中保存shape格式数据时,没有具体的函数可以将投影信息直接写入到shape文件中,不过可以通过另外一种方式实现.看下shape格式的文件不难发现,shape文件的投影信息是一个单独的文 ...

matlab 谱分析函数,科学网—经典谱分析（Power Spectrum Analysis） - 刘磊的博文

matlab 谱分析函数,科学网—经典谱分析（Power Spectrum Analysis） - 刘磊的博文相关推荐

最新文章

热门文章