您需要将频谱图提供给MATLAB功能,并且该频谱图应该具有所需的时间分辨率。
请参阅下面的示例。
CLC,清楚
FS = 20480;
f = 0.01 * fs;
ts = 1 / fs;
T = 0:TS:15-TS;
X = COS(2 * PI * F * T);
%计算功率谱图,时间分辨率等于0.02秒
%注意,频谱图窗口之间的重叠设置为零
%如果您在Windows之间引入重叠,则应分辨率应
%改变以保证每0.02秒采样'tp'
[p,fp,tp] = pspectrum(x,fs,'timeresolution',0.02,'重叠的',0,'谱图');
%计算谱熵
[SE,TE] = PENTROPY(P,FP,TP);
fprintf('光谱熵估计每个%g sec \ n',平均值(差异(te)));