现有一列震动时域数据data11，使用matlab取其峰值出现之前和之后2秒的数据，绘制其0-50Hz的时频图

在MATLAB中，处理振动时域数据并提取峰值前后2秒的数据，然后绘制0-50Hz的时频图（STFT或短时傅里叶变换），可以按照以下步骤操作：

1. 导入数据：首先，你需要导入时间序列数据`data11`，假设它是时间序列向量。

data11 = csvread('your_data.csv'); % 如果是CSV文件，如果是其他格式，请相应修改读取方式

2. 定义信号长度：确定原始信号的时间长度，以及峰点检测的窗口大小（比如2秒）。

signal_length = length(data11);
window_size = 2 * 60; % 2秒转化为秒，这里假设采样率为1Hz（如果不同请调整）

3. 峰值检测：使用MATLAB的内置函数`findpeaks`来找到峰值位置及其附近的区间。

[peaks,locs] = findpeaks(data11, 'MinPeakHeight', data11(1)*0.8); % 高度设置为原始信号最小值的80%
peak_start = locs - floor(window_size/2);
peak_end = peak_start + window_size;

4. 切割信号：创建子信号，包含峰值前后各2秒的数据。

sub_signals = data11(peak_start:peak_end, :);

5. 计算时频图：使用`spectrogram`或`wavenet`函数计算0-50Hz范围内的时频图，可能需要先将信号进行窗函数处理。

window = hann(window_size); % 使用汉明窗或其他适合的窗函数
fs = 1; % 假设每秒一个样本，实际根据采样率调整
[frequencies, times, spectrogram] = spectrogram(sub_signals, window, window_size, fs, 'yaxis');

6. 显示时频图：最后，用`imagesc`或`contourf`等函数展示结果。

xlim([0 max(frequencies)]); % 设置频率轴范围
ylim([0 window_size]); % 设置时间轴范围
title('0-50Hz时频图 (峰值前后2秒)');
xlabel('Frequency [Hz]');
ylabel('Time [samples]');
imagesc(times, frequencies, abs(spectrogram));
colorbar;

**相关问题--:**

1. `findpeaks`函数的作用是什么？
2. 为什么要使用汉明窗或类似的窗函数？
3. 在绘制时频图时，如何指定频率范围？
4. 如何优化时频图的显示效果？