现有一列震动时域数据data111，采样率为1000hz，使用matlab绘制其0-100Hz频域图

在MATLAB中，你可以使用快速傅立叶变换（FFT）来将时域信号data111转换为频域，然后绘制0-100Hz范围内的频率响应。以下是一般的步骤：

1. **导入数据**:
使用`load`或`readmatrix`函数加载数据文件，假设文件名为'data111.mat'：

   data = load('data111.mat');
   data111 = data.data111;

2. **检查数据长度和采样率**:
确保数据长度是采样点数，并且采样率是1000 Hz：

   fs = 1000; % 采样率
   N = length(data111); % 数据点数量

3. **准备频谱分析**:
对于非整数秒的数据，可能需要对数据进行零填充（`zp = round(N/fs * 100)`），以便得到更精确的频率分辨率。然后创建频率向量 `frequencies`:

   zp = round(N/fs * 100);
   frequencies = (0:zp-1) * fs / zp; % 创建频率范围 [0, 50] Hz with 100 points

4. **计算FFT**:
使用`fft`函数计算离散傅立叶变换（DFT），取前`zp`个结果（因为做了零填充）:

   spectrum = abs(fft(data111, zp));

5. **绘制频域图**:
使用`plot`或`semilogx`（以对数尺度显示低频部分）绘制0-100 Hz范围内的频谱，并标出X轴上的频率刻度：

   subplot(2, 1, 1)
   plot(frequencies(1:zp/2), spectrum(1:zp/2)) % 只画半个周期，消除镜像
   title('0-100 Hz 频域图')
   xlim([0 100]) % 显示0-100 Hz区间

   subplot(2, 1, 2)
   semilogx(frequencies(1:zp/2), spectrum(1:zp/2))
   xlabel('Frequency (Hz)')
   ylabel('Magnitude')

6. **调整图像显示**:
根据需要添加网格线、标签等，并确保图形清晰易读。