matlab 阶次分析代码

以下是一个用MATLAB实现阶次分析的简单示例代码：

% 生成示例振动信号
fs = 1000;  % 采样频率
t = 0:1/fs:1;  % 时间序列
f1 = 50;  % 主频率
f2 = 100;  % 混叠频率
x = sin(2*pi*f1*t) + 0.5*sin(2*pi*f2*t);

% 阶次分析
nfft = 2^12;  % FFT窗口长度，取2的幂次方方便计算
spectrum = fft(x, nfft);  % FFT变换
freq = (0:nfft-1)*(fs/nfft);  % 频率轴
order = freq / (f1/60); % 计算阶次

% 绘制阶次谱
figure;
plot(order, abs(spectrum));
xlabel('Order');
ylabel('Amplitude');
title('Order Spectrum');

这段代码首先生成了一个示例的振动信号，然后进行了阶次分析。使用FFT变换对振动信号进行频谱分析，计算频率轴和阶次，然后绘制阶次谱。你可以根据实际需求修改信号参数和绘图方式。