请写出order tracking的matlab代码

以下是一个简单的order tracking的MATLAB代码示例：

% 生成测试信号
fs = 1000; % 采样率
f1 = 50; % 基频频率
f2 = 100; % 旋转频率
t = 0:1/fs:1; % 时间向量
x = sin(2*pi*f1*t + 0.5*pi*sin(2*pi*f2*t)); % 信号

% 阶次跟踪
nfft = length(x); % FFT点数
freq = (0:nfft-1)*fs/nfft; % 频率向量
X = fft(x, nfft); % FFT
X = X(1:nfft/2); % 取FFT结果的前半部分
X_mag = abs(X); % 取幅值
[~, max_index] = max(X_mag); % 取最大幅值对应的频率
f = freq(max_index); % 最大幅值对应的频率
N = round(fs/f); % 阶次
X_order = zeros(size(X)); % 初始化阶次谱
for i = 1:N:length(X)-N+1
    X_order(i:i+N-1) = X(i:i+N-1); % 取该阶次范围内的FFT结果
end
x_order = ifft(X_order, nfft); % 逆FFT

% 绘制结果
subplot(2,1,1); plot(t, x); xlabel('时间'); ylabel('幅度');
title('原始信号');
subplot(2,1,2); plot(t, real(x_order)); xlabel('时间'); ylabel('幅度');
title(['阶次跟踪结果（阶次:', num2str(N), ')']);

这段代码首先生成一个包含基频和旋转频率的测试信号，然后使用FFT方法计算信号的频谱，找到最大幅值对应的频率，进而确定信号的阶次。接下来，根据阶次取出相应的FFT结果，再使用逆FFT方法得到阶次跟踪后的信号。最后，将原始信号和阶次跟踪结果绘制在同一张图上进行对比。