如何使用MATLAB编程，针对门信号f(t)=g2(t)，其定义为-1到1范围内的函数，在-1到1的时间间隔内等于1，其余时间等于0，进行均匀采样，并分别展示采样前后的原始信号和对应的频谱图？请包括详细的步骤和关键代码示例。

在MATLAB中，你可以按照以下步骤对门信号f(t)进行均匀采样并绘制原信号和频谱图：

1. **创建门信号函数**:
使用`heaviside`函数表示门信号，因为`g2(t)`可以看作是`heaviside(t+1)-heaviside(t-1)`。这里假设`t`是从-1到1的线性向量。

t = -1:0.001:1; % 时间向量，我们选择0.001作为采样间隔
f = heaviside(t+1) - heaviside(t-1); % 创建门信号函数

2. **采样信号**:
使用`sample`函数对连续信号进行等间隔采样。`f`是离散化前的连续信号，采样率取决于时间向量`t`的步长。

fs = length(t)/duration(t); % 计算采样频率
f_sampled = sample(f, fs); % 等间隔采样

3. **绘制原信号波形**:
使用`plot`函数绘制采样前后的原始信号。

figure;
subplot(2,1,1);
plot(t, f, 'b', 'LineWidth', 2); % 绘制原始连续信号
xlabel('Time (s)');
ylabel('Amplitude');
title('Original Signal');

subplot(2,1,2);
plot(t, f_sampled, 'r', 'LineWidth', 2); % 绘制采样后的信号
xlabel('Time (s)');

4. **计算并显示频谱**:
对采样后的信号，我们可以使用`fft`函数计算傅里叶变换，然后用`abs`获取幅度谱，最后绘制出来。

% 对采样后的信号做FFT
F = fft(f_sampled);

% 将结果转换为频率轴
Fs = fs / length(t); % 总采样频率
frequencies = (0:length(F)-1)*Fs/(length(F)-1);

% 取幅值，得到频谱图
spectrogram = abs(F ./ length(F));

% 绘制频谱
figure;
imagesc(frequencies, t, spectrogram);
xlabel('Frequency (Hz)');
ylabel('Time (s)');
colorbar;
title('Spectrum of Sampled Signal');

注意：为了更好地看到频谱，通常会将`spectrogram`的纵坐标反转，但在本例子中为了直观起见，我们保留了原始的顺序。