现有采样率为1000hz的震动时域数据data，使用matlab取其最大峰值出现之前和之后2秒的数据，绘制其0-50Hz的时频图，请试运行没有错误后再输出代码

在MATLAB中，你可以按照以下步骤操作来完成这个任务：

首先，确保你的时域数据`data`是一个长度为采样点数（1000 Hz * 时间长度）的一维向量。然后，我们需要确定最大峰值的位置，以及提取前后2秒的数据。接着，我们将这部分数据用于绘制0-50 Hz范围内的时频图。这里假设时间长度已知，用`T`表示。

% 假设data是你的原始振动时域信号
Fs = 1000; % 采样率
time_length = length(data) / Fs; % 计算总时间

% 找到最大峰值位置
[~, idx_max] = max(abs(data)); % abs()用于处理可能的负值

% 提取最大峰值前后的2秒数据
t_start = floor(idx_max - 2*Fs); % 前2秒的开始索引，向下取整
t_end = ceil(idx_max + 2*Fs);   % 后2秒的结束索引，向上取整
if t_start < 1 || t_end > length(data)
    error('超出数据边界');
end
data_subsampled = data(t_start:t_end);

% 确保数据长度足够计算50 Hz带宽内的时频图
if t_end - t_start < Fs/2
    warning('提取的数据不足，无法绘制完整的0-50 Hz时频图');
else
    % 使用短时傅立叶变换（STFT）绘制时频图
    win_size = round(Fs/50); % 50 Hz窗口大小
    hop_size = win_size / 2; % 步进大小
    [Y, f] = stft(data_subsampled, win_size, hop_size, Fs);
    
    % 取频率范围0-50 Hz内的部分
    Y_bandlimited = Y(:, find(f <= 50));
    
    % 绘制时频图
    figure;
    imagesc(unwrap(angle(Y_bandlimited)), f(1:end/2), time_length*(t_start/timeshift(win_size,Fs)));
    title('0-50 Hz时频图 (最大峰值附近)');
    xlabel('时间 (s)');
    ylabel('频率 (Hz)');
    colormap(jet);
    colorbar;
end

执行上述代码之前，请确保你的MATLAB环境已经安装了Signal Processing Toolbox。如果遇到问题或错误，根据提示调整参数或检查输入数据是否正确。如果你没有安装该工具箱，可能需要购买或通过其他方法访问该功能。