matlab 案例 订单分配,如何在MATLAB中实现订单分析

我有一个代表电机声音的音频文件.我一直在读“正常”fft没有对机器进行有价值的分析,我应该使用Order分析来描述机器的“行为”.我在增加速度的同时记录了它的声音,我将记录放在一起并计算出它的频谱.这是代码：

%% read the Audio file and plot its

clear, clc , clf;

M = 512;

k =1

data= 0 ;

for i =100:5:180

fileName =['A10_uSp0_Mic100k_2nd_V_',int2str(i),'.wav']

[x(:,k),Fs] = audioread(fileName);

k = k+1 ;

end

for i = 1:length(x(1,:))

data = [data(:);x(:,i)];

end

k = length(x(1,:));

while k ~= 2

data = [data(:);x(:,k)];

k = k -1;

end

spectrogram(data,hamming(M),M/2,0:20:4000,Fs,'yaxis');

结果如下：

正如您可能看到的那样,我为每个时间点将机器的速度提高了5.我的问题是如何计算该系统的订单分析！？

非常感谢您的帮助！