用matlab画发动机声音阶次的colormap图

以下是一个简单的matlab示例，用于绘制发动机声音阶次的colormap图：

% 生成数据
fs = 44100; % 采样率
f0 = 100; % 基频
N = 2^14; % 采样点数
t = (0:N-1)/fs; % 时间向量
x = sin(2*pi*f0*t); % 基波信号
for i = 2:8 % 生成8个阶次分量
    xi = sin(2*pi*i*f0*t)/i;
    x = x + xi;
end
% 计算STFT谱
window = 1024; % 窗口长度
noverlap = window/2; % 重叠长度
nfft = window; % FFT长度
[S,F,T] = spectrogram(x,window,noverlap,nfft,fs);
% 绘制colormap图
imagesc(T,F,20*log10(abs(S))); % 使用dB单位
axis xy; % 纵轴翻转
xlabel('Time (s)');
ylabel('Frequency (Hz)');
colorbar; % 显示colorbar

运行以上代码，即可得到发动机声音阶次的colormap图。需要注意的是，这里使用了STFT（Short-Time Fourier Transform）来计算频谱，因此colormap图中的横轴表示时间，纵轴表示频率。