matlab图像频谱分析代码_信号频域分析方法的理解（频谱、能量谱、功率谱、倒频谱、小波分析）...

信号的频域分析方法多种多样，这里针对较为常见的（频谱、能量谱、功率谱、倒频谱、小波分析）集中进行说明。这些方法的MATLAB代码实现参见文章频域特征值提取的MATLAB代码实现（频谱、功率谱、倒频谱）

在看这篇文章之前可以参看之前的两篇，其中涉及一些时域特征值介绍和能量、功率信号区别的介绍：

Mr.括号：信号的各种时域分析方法的理解

Mr.括号：能量信号和功率信号的分别

文章如要转载请私信与我联系，并注明来源知乎专栏与信号处理有关的那些东东作者Mr.括号。

一、频谱

一般来说，频谱分析指的是将信号做傅里叶变换从而进行分析。频谱分析是包括幅频谱和相频谱两张图的。不过最常用的是幅频谱。

频谱分析是最常用和最重要也是最基础的频域分析方法。在这里不详细展开说了，不了解概念的可以参考这篇文章：

Heinrich：如果看了这篇文章你还不懂傅里叶变换，那就过来掐死我吧

时常会见到几个概念：DFT，DTFT，DFS，FFT，FT, FS等等。这些概念怎么区分？看这篇文章就能搞懂：一幅图弄清DFT与DTFT,DFS的关系 - BitArt - 博客园

在这里概括一下（不感兴趣的话可以直接跳过）：

FT（Fourier Transformation）：傅里叶变换。就是我们理论上学的概念，但是对于连续的信号无法在计算机上使用。其时域信号和频域信号都是连续的。

DTFT（Discrete-time Fourier Transform）：离散时间傅里叶变换。这里的“离散时间”指的是时域上式离散的，也就是计算机进行了采样。不过傅里叶变换后的结果依然是连续的。

DFT（Discrete Fourier Transform）：离散傅里叶变换。在DTFT之后，将傅里叶变换的结果也进行离散化，就是DFT。

也就是说：FT时域连续、频域连续；DTFT时域离散、频域连续；DFT时域离散、频域离散。

FFT（Fast Fourier Transformation）：快速傅里叶变换。就是DFT的快速算法，一般工程应用时用的都是这种算法。

FS（Fourier Series）：傅里叶级数。是针对时域连续周期信号提出的，结果是离散的频域结果。

DFS（Discrete Fourier Series）：离散傅里叶级数。是针对时域离散周期信号提出的，DFS与DFT的本质是一样的。

另外补充几点相关知识：

• 在实际计算中通常使用快速傅里叶变换（FFT）。它是一种用来计算DFT（离散傅里叶变换）和IDFT（离散傅里叶反变换）的一种快速算法。
• 随机信号是无法做傅里叶变换的（*这里要再补充）

二、能量谱

要理解能量谱和功率谱，首先要弄明白能量有限信号和功率有限信号（参看之前的文章能量信号和功率信号的分别）。

能量谱也叫能量谱密度，能量谱密度描述了信号或时间序列的能量如何随频率分布。能量谱是原信号傅立叶变换的平方。

三、功率谱

功率谱是功率谱密度函数（PSD）的简称，它定义为单位频带内的信号功率。

功率谱是针对功率信号来说的。功率谱的推导公式相对复杂，不过幸运的是维纳-辛钦定理证明了：一段信号的功率谱等于这段信号自相关函数的傅里叶变换。

所以求功率谱就有了两种方法：1.(傅立叶变换的平方)/(区间长度)；2.自相关函数的傅里叶变换。这两种方法分别叫做直接法和相关函数法。

功率谱这里存在着一些问题，整理如下：

1.功率谱密度的单位是什么，看有的写的是dB，还有的说是W/Hz。