谱密度

的谱密度是光波长的函数，箭头所指的原子跃迁处会出现峰值。]]

時間序列 x(t) 的功率谱 S_{xx}(f) 描述了信号功率在频域的分布状况。根据傅里叶分析，任何物理信号都可以分解成一些离散频率或连续范围的频谱。对特定信号或特定种类信号（包括噪声）频率内容的分析的统计平均，称作其频谱。

当信号的能量集中在一个有限时间区间的时候，尤其是总能量是有限的，就可以计算能量谱密度。更常用的是应用于在所有时间或很长一段时间都存在的信号的功率谱密度。由于此种持续存在的信号的总能量是无穷大，功率谱密度（英語：Power Spectral Density，縮寫PSD）则是指单位时间的频谱能量分布。频谱分量的求和或积分会得到（物理过程的）总功率或（统计过程的）方差，这与帕塞瓦尔定理描述的将 x^2(t) 在时间域积分所得相同。

物理过程 x(t) 的频谱通常包含与 x 的性质相关的必要信息。比如，可以从频谱分析直接确定乐器的音高和音色。电磁波电场 E(t) 的频谱可以确定光源的颜色。从这些时间序列中得到频谱就涉及到傅里叶变换以及基于傅里叶分析的推广。许多情况下时间域不会具体用在实践中，比如在攝譜儀用散射棱镜来得到光谱，或在声音通过内耳的听觉感受器上的效应来感知的过程，所有这些都是对特定频率敏感的。

不过本文关注的是时间序列（至少在统计意义上）已知，或可以直接测量（如经麦克风采集再由电脑抽样）的情形。功率谱在与随机过程的统计研究以及物理和工程中的许多其他领域中都很重要。通常情况下，该过程是时间的函数，但也同样可以讨论空间域的数据按空間頻率分解。

解释
在物理学中，信号通常是波的形式，例如电磁波、随机振动或者声波。当波的频谱密度乘以一个适当的系数后将得到每单位频率波携带的功率，这被称为信号的功率谱密度（power spectral density, PSD）或者谱功率分布（spectral power distribution, SPD）。功率谱密度的单位通常用每赫兹的瓦特数（W/Hz）表示，或者使用波长而不是频率，即每纳米的瓦特数（W/nm）来表示。

尽管并非一定要为信号或者它的变量赋予一定的物理量纲，下面的讨论中假设信号在时域内变化。

定义
能量谱密度
能量谱密度描述的是信号或者时间序列的能量如何随频率分布。这里，能量这个术语是用作信号处理中的推广含义； 也就是说，信号 x(t) 的能量 E 为
: E =\int\limits_{-\infty}^\infty |x(t)|^2\, dt.
能量谱密度对总能量有限的瞬变信号（也就是类似于脉冲信号的）最为适用。在这种情况下，帕塞瓦尔定理给出了用傅里叶变换 \hat{x}(f)=\int\limits_{-\infty}^\infty e^{-2\pi ift}x(t) dt 表示信号能量的形式。

:\int\limits_{-\infty}^\infty |x(t)|^2\, dt = \int\limits_{-\infty}^\infty |\hat{x}(f)|^2\, df.

这里频率 f 单位为Hz，即每秒周期数。经常使用角频率 \omega=2\pi f。由于右边的积分是信号的能量，被积函数 |\hat{x}(f)|^2 可以理解为频率为 f 的信号中单位频率包含的能量的密度函数。鉴于此，信号 x(t) 的能量谱密度定义为
）仍旧使用非归一化的傅里叶变换来定义功率谱密度

: \langle \hat x(\omega) \hat x^\ast(\omega') \rangle = 2\pi\,f(\omega)\,\delta(\omega-\omega'),

其中 \delta(\omega-\omega') 为狄拉克δ函数。这些形式陈述有时对引导直觉会比较有用，但要极其谨慎地使用。 由于这种复杂性，可以用仅仅在有限区间 [0, T] 把信号积分的截短傅里叶变换 \hat{x}_T(\omega)：

: \hat{x}_T(\omega) = \frac{1}{\sqrt{T}} \int_0^T x(t) e^{-i\omega t}\, dt.

因此功率谱密度可以被定义为

: S_{xx}(\omega) = \lim_{T \rightarrow \infty} \mathbf{E} \left[ | \hat{x}_T(\omega) | ^ 2 \right].

这里 E 表示期望值；明确地，我们有

给定频带 [f_1, f_2]（或[\omega_1,\omega_2]）中信号的功率可以通过对频率积分计算。由于 S_{xx}(-\omega) = S_{xx}(\omega)，正、负频率的功率相同，因而下面形式中的因子为2（这种因子取决于使用的惯例）：
: P_\mathsf{bandlimited} =
2 \int_{f_1}^{f_2}\,S_{xx}(2\pi \! f) \, df =
\frac{1}{\pi} \int_{\omega_1}^{\omega_2}\,S_{xx}(\omega) d\omega

更一般地，类似的技术可以被用来估计一个随时间变化的光谱密度。更一般地，可以使用类似的方法来估计时变谱密度。在这种情况下上面定义的 (0, T) 上的截短傅里叶变换不是通过 T 趋近于无穷的极限计算的。这导致光谱覆盖率和分辨率降低，因为不会采样小于 1/T 的频率，而 1/T 的整数倍频率的结果不是独立的。

性质

• f(t) 的谱密度和 f(t) 的自相关组成一个傅里叶变换对（对于功率谱密度和能量谱密度来说，使用着不同的自相关函数定义）。

• 通常使用傅里叶变换技术估计谱密度，但是也可以使用如和最大熵这样的技术。

• 傅里叶分析的结果之一就是帕塞瓦爾定理，这个定理表明能量谱密度曲线下的面积等于信号幅度平方下的面积，总的能量是：

::\int_{-\infty}^\infty \left| f(t) \right|^2 dt = \int_{-\infty}^\infty \Phi(\omega) \,d\omega.

：上面的定理在离散情况下也是成立的。另外的一个结论是功率谱密度下总的功率与对应的总的平均信号功率相等，它是逐渐趋近于零的自相关函数。

相关概念

• 大多数“频率”图实际上仅仅表示了谱密度。有时完整的频率要用两部分来表示，一部分是对应于频率的“幅度”（它就是谱密度），另外一部分是对应于频率的“相位”（它包含了频谱中剩余的其它信息）。信号 f(t) 可以从一个完整的频谱进行恢复。需要注意的是 f(t) 不能仅仅从谱密度这一部分进行恢复——它丢失了“临时信息”。

• 信号的 是谱密度函数的中点，也就是说将整个分布切分成两个相等部分的点。

• 谱密度是频率的函数，而不是时间的函数。但是，也可以计算一个较长信号上一小段“窗口”的谱密度，并且根据与事件相关的窗口进行绘图，这样的图形称为频谱图。这是短时傅里叶变换和小波等许多谱分析技术的基础。

应用
电子工程
信号功率谱的概念和应用是电子工程的基础，尤其是在电子通信系统中，例如无线电和微波通信、雷达以及相关系统。人们已经花费了很大的精力和大量的金钱投入到开发、生产“频谱分析仪”这种电子设备，用来帮助电子工程师、技术人员、技工观察、测量电子信号的功率谱。频谱分析仪的价格根据带宽和精度的不同而不同，质量最好的仪器的价格超过 100,000 美元。

色度学
光源的频谱是每个频率携带的功率或者光源中“颜色”的度量。光谱通常是沿着可见光在波长空间而不是频率空间测量的不同点（通常是 31 个点）进行测量，它不是严格意义上的谱密度。一些分光光度计能够分辨高达 1 到 2 纳米 的增量精度，测量值用来计算其它的规格然后绘制出来显示光源的频谱属性。这对于分析特定光源的颜色特性来说是一个非常有用的工具。

参见

• 噪声的颜色
• 谱泄漏
• 窗函数
• 频域
• 频谱
• 双谱

参考资料
外部链接

• [https://web.archive.org/web/20060716045107/http://documents.wolfram.com/applications/timeseries/UsersGuidetoTimeSeries/1.8.1.html Power Spectral Density Function] - from Wolfram Research's "Time Series Documentation"