现有采样率为1000hz的震动时域数据data1，使用matlab取其峰值出现之前和之后2秒的数据，绘制其0-50Hz的时频图

在MATLAB中，你可以使用以下步骤来处理1000Hz采样的震动时域数据`data1`，提取峰值前后各2秒的数据，并绘制出0-50Hz的时频图（Short-Time Fourier Transform, STFT）：
1. **预处理数据**：
- 首先，确保数据已经正确加载到MATLAB中。如果数据是数组形式，假设它是一个向量。

   data = data1;

2. **确定峰点位置**：
- 使用`findpeaks`函数找到峰值的位置，例如，峰值高度设定为某个阈值。

   [peaks,locs] = findpeaks(data, 'MinPeakHeight', threshold);

3. **计算峰值前后2秒的数据**：
- 根据`locs`找出峰值索引，然后选择每处峰值前后的2秒数据。注意MATLAB的时间尺度是以采样率表示的，所以2秒对应的是1000 * 2个样本。

   start_index = locs - 2000; % 峰值前2秒
   end_index = locs + 2000; % 峰值后2秒
   valid_data = data(start_index:end_index);

4. **窗口ing和STFT**：
- 使用`spectrogram`或`stft`函数，设置窗口大小（如512点），步长（如256点），并限制频率范围为0-50Hz。

   window_size = 512; % 窗口长度
   hop_length = 256; % 步进长度
   fs = 1000; % 采样率
   fmax = 50; % 频率上限

   Y = stft(valid_data, window_size, hop_length, [], fs, [-fmax:fmax]);

5. **绘制时频图**：
- 最后，使用`imagesc`或` spectrogram`函数将结果可视化。

   figure;
   imagesc(unwrap(angle(Y)), logspace(-1, 1, size(Y, 2)), linspace(0, fs-length(valid_data)/hop_length, size(Y, 1)));
   colorbar;
   xlabel('时间 (s)');
   ylabel('频率 (Hz)');
   title('0-50Hz时频分析');

完成以上步骤后，你应该会看到一个显示了峰值前后2秒内0-50Hz振动信号变化的时频图。

**相关问题**:
1. MATLAB中的`findpeaks`函数有哪些参数可以调整？
2. `spectrogram`函数的输入参数如何影响时频图的生成？
3. 如何调整`imagesc`的颜色范围以更好地展示低频和高频信息？