注意:这个函数现在可以从IoSR iosr.dsp.ltas Matlab工具箱。
- - - - - - - - - - - - - - - - - - - - - - - - -
S =节目(X, FS)计算信号的长期平均频谱(节目)X,采样在FS赫兹。的频谱计算平均功率谱密度(PSD)获得一系列重叠的fft算法;FFT长度是4096,跳的大小是2048。Hann-windowed X的片段。然后平均PSD Gaussian-smoothed 1/3-octave决议。
X是一个向量,矩阵,或多维数组;节目将在第一个non-signleton维度,并返回节目每个对应的行/列/ etc。
S =节目(X, FS,“参数”,“价值”)允许大量的参数被指定。这些参数有:-
“暗”:{找到(大小(X) > 1, 1, '第一次')}|标量
指定操作的维数(默认为第一个单体尺寸)。
“图”:{假}|真的
选择是否策划节目的图。
“跳”:{NFFT / 2} |标量
指定步长通过X用来计算每一部分。NFFT是由“赢得”参数。
“noct”:{3} |标量
应用1 / noct-octave平滑的频谱。“noct”设置为0的结果不平滑。
“缩放”:{‘没有’}|“max0”
指定任何扩展适用于美国在默认情况下,没有扩展应用。如果缩放设置为“max0”,将比例有一个最大的值为0。
“单位”:{dB} |“没有”
指定输出单元。默认情况下,在dB PSD计算。否则PSD直接返回。
“赢得”:{4096}|标量|向量
指定了窗口或FFT长度NFFT用来计算频谱。如果“赢”是一个标量,它指定了FFT长度和损害窗口是应用于每一段。如果“赢”是一个矢量,NFFT向量的长度,向量乘以每段。
例子
%绘制1/6th-octave-smoothed节目韩德尔的例子
负载handel.mat
图
节目(y, Fs,“noct”6“图”,真正的);
参见STFT, SMOOTH_SPECTRUM。
引用作为
克里斯托弗Hummersone (2022)。长期平均频谱GitHub (https://github.com/IoSR-Surrey/MatlabToolbox)。检索。
