频谱分析(完整版).docxVIP

Matlab 信号处理工具箱帮助文档谱估计专题翻译：无名网友 Lyra频谱分析Spectral estimation（谱估计）的目标是基于一个有限的数据集合描述一个信号的功率（在频率上的）分布。功率谱估计在很多场合下都是有用的，包括对宽带噪声湮没下的信号的检测。从数学上看，一个平稳随机过程的power spectrum（功率谱）和correlation sequence（相关序列）通过discrete-time Fourier transform（离散时间傅立叶变换）构成联系。从normalized frequency（归一化角频率）角度看，有下式注：，其中。其matlab近似为X=fft(x,N)/sqrt(N)，在下文中就是指matlab fft函数的计算结果了使用关系可以写成物理频率的函数，其中是采样频率相关序列可以从功率谱用IDFT变换求得：序列在整个Nyquist间隔上的平均功率可以表示为上式中的以及被定义为平稳随机信号的power spectral density (PSD)（功率谱密度）一个信号在频带上的平均功率可以通过对PSD在频带上积分求出从上式中可以看出是一个信号在一个无穷小频带上的功率浓度，这也是为什么它叫做功率谱密度。PSD的单位是功率（e.g 瓦特）每单位频率。在的情况下，这是瓦特/弧度/抽或只是瓦特/弧度。在的情况下单位是瓦特/赫兹。PSD对频率的积分得到的单位是瓦特，正如平均功率所期望的那样。对实信号，PSD是关于直流信号对称的，所以的就足够完整的描述PSD了。然而要获得整个Nyquist间隔上的平均功率，有必要引入单边PSD的概念：信号在频带上的平均功率可以用单边PSD求出频谱估计方法Matlab 信号处理工具箱提供了三种方法Nonparametric methods（非参量类方法）PSD直接从信号本身估计出来。最简单的就是periodogram（周期图法），一种改进的周期图法是Welchs method。更现代的一种方法是multitaper method（多椎体法）。Parametric methods （参量类方法）这类方法是假设信号是一个由白噪声驱动的线性系统的输出。这类方法的例子是Yule-Walker autoregressive (AR) method和Burg method。这些方法先估计假设的产生信号的线性系统的参数。这些方法想要对可用数据相对较少的情况产生优于传统非参数方法的结果。Subspace methods（子空间类）又称为high-resolution methods（高分辨率法）或者super-resolution methods（超分辨率方法）基于对自相关矩阵的特征分析或者特征值分解产生信号的频率分量。代表方法有multiple signal classification (MUSIC) method或eigenvector (EV) method。这类方法对线谱（正弦信号的谱）最合适，对检测噪声下的正弦信号很有效，特别是低信噪比的情况。方法描述函数周期图PSD 估计spectrum.periodogram, periodogramWelch重叠，加窗的信号段的平均周期图spectrum.welch, pwelch, cpsd, tfestimate, mscohere多椎体多个正交窗（称为锥）的组合做谱估计spectrum.mtm, pmtmYule-WalkerAR时间序列的估计的自相关函数计算自回归（AR）谱估计spectrum.yulear, pyulearBurg通过最小化线性预测误差计算自回归（AR）谱估计spectrum.burg, pburgCovariance（协方差）通过最小化前向预测误差做时间序列的自回归（AR）谱估计spectrum.cov, pcov修正协方差通过最小化前向及后向预测误差做时间序列的自回归（AR）谱估计spectrum.mcov, pmcovMUSIC多重信号分类spectrum.music, pmusic特征向量法虚谱估计spectrum.eigenvector, peigNonparametric Methods非参数法下面讨论periodogram, modified periodogram, Welch, 和 multitaper法。同时也讨论CPSD函数，传输函数估计和相关函数。Periodogram周期图法一个估计功率谱的简单方法是直接求随机过程抽样的DFT，然后取结果的幅度的平方。这样的方法叫做周期图法。一个长L的信号的PSD的周期图估计是注：这里运用的是matlab里面的fft的定义不带归一化系数，所以要除以L其中实际对的计算可以只在有限的频率点上执行并且使用FFT。实践上大多数周期图法的应用都计算N点PS