信号的频谱、幅度谱、相位谱及能量谱密度、功率谱密度
摘录别人的,因为原始博客公式看不了。下面是原地址。 https://www.cnblogs.com/iliveido/archive/2013/03/22/2976542.html
傅里叶变换一个令人震惊的事实是:Gaussian分布的密度函数 \(e^{-x^2/2}\) 是唯一的一个傅里叶变换不变函数。
泛函分析中,Gaussian密度函数的极限(\(\sigma\to\infty\))是delta-dirac函数 \(\delta(x)\),即脉冲函数。
更简单地,在大学一年级的数学分析课程中,Gaussian密度函数的积分是 \(\sqrt{\pi}\)。
信号经过傅里叶变换之后产生频谱,频谱是一个以频率为自变量的函数。频谱在每一个频率点的取值是一个复数。一个复数由模和辐角唯一地确定,即:
\[
z = r(cos\theta + isin\theta)
\]
所以可将频谱分解为幅度谱(即复数的模关于频率的函数)和相位谱(即复数的辐角关于频率的函数)。
那什么是能量谱密度(energy spectral density)和功率谱密度(power spectral density)?
在英语中,幅度有两个词:amplitude和magnitude,在大多数情况下(包括本文),它们是没有区别的,除了在某个特定的领域(如物理领域),amplitude代表整个信号偏离x轴的最大绝对值,magnitude代表信号上某一点偏离x轴的绝对值。更清晰的阐述如下:
peak amplitude, often shortened to amplitude, is the nonnegative value of the waveform's peak (either positive or negative).
instantaneous amplitude of x is the value of x(t) (either positive or negative) at time t.
instantaneous magnitude, or simply magnitude, of x is nonnegative and is given by |x(t)|.
可见,amplitude是一个全局概念,而magnitude是一个瞬时概念。
在谈及FFT和wavelets的大多数情况下,将amplitude和magnitude认为是同一概念,即瞬时幅度。
信号\(f(t)\)在\(t\)处的瞬时幅度是\(f(t)\)的模,即\(|f(t)|\);
信号\(f(t)\)在\(t\)处的瞬时相位是\(f(t)\)的辐角,即\(Arg f(t)\) 或者 \(\angle f(t)\);
信号\(f(t)\)在\(t\)处的瞬时功率是\(f(t)\)的模的平方,即\(|f(t)|^2\)。
信号的能量是一个全局概念,是瞬时功率的积分值,即
\[ ||f(t)||^2=\int_{-\infty}^{\infty}|f(t)|^2 dt \]
注意\(|f(t)|\)和\(||f(t)||\)的区别,前者是瞬时概念,即信号在某一点的瞬时幅度,后者是全局概念,即整个信号的能量的开方。
需要注意的是,通常所指的能量谱和能量谱密度是一个概念;功率谱和功率谱密度是一个概念,而且功率是指平均功率。
时域上的能量公式:
\[ E(f)=\int_{-\infty}^{\infty}|f(t)|^2 dt \]
其中绝对值号代表取模,当信号是实信号时,显然绝对值号可以去掉,变成
\[
E(f)=\int_{-\infty}^{\infty}f(t)^2 dt
\]
根据Parseval能量恒等式(Parseval’s Identity),能量也可认为是\(f(t)\)的傅里叶变换的模的平方在频域上的积分。
频域上的能量公式:
\[
E(f)=\int_{-\infty}^{\infty}|\hat{f}(\omega)|^2 d\omega
\]
从上述积分可以看出,信号的能量谱密度在某个频率点上的取值就是信号在某个频率上的瞬时功率\(|\hat{f}(\omega)|^2\)。
从上面的公式可以看出,信号的能量可能是无穷。当信号的能量无限时,只能通过平均功率来了解该信号。因为能量\(E(f)\)和时间长度\(\triangle T\)之比就是平均功率\(P(f)\),即:
\[
P(f)=\frac{E(f)}{\triangle T}
\]
易知:当信号在\(t \in (-\infty,\infty)\)的平均功率有限时,能量是无限的;当信号在\(t \in (-\infty,\infty)\)的能量有限时,其平均功率为0。能量有限的信号称为能量信号;平均功率有限的信号称为功率信号。
为方便叙述,记
\[
f_{T}(t)=\left\{\begin{array}{lll} f(t)& , & |t|\leq T\\ 0&,&|t|>T\end{array}\right.
\]
从而平均功率的公式为:
\[
P(f)=\lim_{T\to\infty}\frac{\int_{-T}^{T}|f_{T}(t)|^2 dt}{2T}=\lim_{T\to\infty}\frac{\int_{-T}^{T}|\hat{f_{T}}(\omega)|^2d\omega}{2T}
\]
从上述的积分可以看出,信号的功率谱密度为:
\[
PSD(f) = \lim_{T\to\infty}\frac{|\hat{f_{T}}(\omega)|^2}{2T}
\]
对于比信号更复杂的随机过程\(X(t)\)来说,\(P(f)\)是一个随机变量,所以其平均功率\(P\)必须取加权平均\(E\)(注意这里的\(E\)不是能量):
\[
P=E[P(f)]=E[\lim_{T\to\infty}\frac{\int_{-T}^{T}|\hat{f_{T}}(\omega)|^2d\omega}{2T}]=\lim_{T\to\infty}\frac{\int_{-T}^{T}E[|\hat{f_{T}}(\omega)|^2]d\omega}{2T}
\]
其功率谱密度为:
\[
PSD = \lim_{T\to\infty}\frac{E[|\hat{f_{T}}(\omega)|^2]}{2T}
\]