如何在MATLAB中实现阶次分析

Question

我有一个代表马达声音的音频文件。我一直在读到“正常” fft 不能对机器进行有价值的分析，我应该使用 订单分析 来描述机器的“行为”。我在增加它的速度的同时记录了它的声音，我把这些记录放在一起并计算了它的频谱图。这是代码：

%% 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。我的问题是如何计算这个系统的阶次分析！？ 非常感谢您的帮助！

Answer 1

这更像是一个科学问题，而不是关于编程的问题。阶次分析可以使用傅里叶变换来完成。您需要做的就是寻找输出中的对角线，即噪声频率与旋转速度成正比的位置。

在您的示例图像中，有一个从 2000 Hz 时间（速度）0 到 4000 Hz 时间（速度）150 的顺序。之后它似乎向下反映，可能是采样不足的伪影。

那么它是哪个顺序？为此，您需要知道您似乎不知道的电机转速。还有另一个较弱的命令，从 3000 Hz 开始，时间（速度）0。通过两个或更多可见命令，您可能能够推断出电机的旋转速度（通过要求频率是原始速度的整数倍） .

除了更高的频率分辨率之外，代码可能不需要更改。