matlab 旋转机械的阶次分析生成三维瀑布图

以下是一种可能的 MATLAB 代码实现，用于旋转机械的阶次分析并生成三维瀑布图：

% 定义旋转速度、采样率和采样时间
rpm = 1200; % rpm
fs = 2048; % Hz
T = 10; % s

% 生成时间和角度向量
t = linspace(0, T, T*fs);
theta = 2*pi*rpm*t/60;

% 生成测试信号
f1 = 50; % Hz
f2 = 100; % Hz
f3 = 150; % Hz
s = sin(2*pi*f1*theta) + sin(2*pi*f2*theta) + sin(2*pi*f3*theta);

% 进行FFT并计算阶次谱
N = length(s);
f = (0:N-1)*fs/N;
freq = f(f <= fs/2);
Y = abs(fft(s)/N);
amp = Y(1:N/2+1);
amp(2:end-1) = 2*amp(2:end-1);
order = freq/rpm;
order_amp = zeros(size(order));
for i = 1:length(order)
    idx = find(abs(f-order(i)*rpm)<=5, 1);
    order_amp(i) = max(amp(idx-5:idx+5));
end

% 生成三维瀑布图
figure;
waterfall(order, freq, amp);
xlabel('阶次');
ylabel('频率 (Hz)');
zlabel('幅值');
title('旋转机械的阶次分析');

% 生成阶次谱
figure;
plot(order, order_amp);
xlabel('阶次');
ylabel('幅值');
title('旋转机械的阶次谱');

该代码首先定义了旋转速度、采样率和采样时间，并生成了时间和角度向量。然后生成了一个测试信号，包含三个不同频率的正弦波。接下来进行FFT并计算阶次谱，最后生成了三维瀑布图和阶次谱。注意，在生成瀑布图时，需要将频率作为y轴，而不是x轴。