生成振速时程离散点的时频三维幅值图MATLAB代码

生成振速时程的时频三维幅值图可以使用MATLAB中的STFT（短时傅里叶变换）函数和surf（三维曲面绘制）函数来实现。以下是一个简单的示例代码：

% 生成振速时程离散点的时频三维幅值图
Fs = 1000; % 采样频率
t = 0:1/Fs:2; % 时间向量
x = sin(2*pi*100*t) + sin(2*pi*200*t) + sin(2*pi*300*t); % 生成信号
N = length(x); % 信号长度
window = hamming(128); % 窗函数
noverlap = 120; % 重叠长度
nfft = 256; % FFT长度
[S,F,T] = spectrogram(x,window,noverlap,nfft,Fs); % 计算STFT谱图
figure;
surf(T,F,abs(S)); % 绘制时频三维幅值图
xlabel('Time (s)');
ylabel('Frequency (Hz)');
zlabel('Magnitude');

在这个示例中，我们首先生成一个包含三个正弦波的信号，然后使用MATLAB中的spectrogram函数计算信号的STFT谱图。最后，我们使用surf函数将STFT谱图绘制成时频三维幅值图。你可以根据自己的需要调整窗函数、重叠长度和FFT长度等参数来得到更好的结果。