使用matlab生成一个矩形脉冲信号，画出其时域图形以及频域图形，使其频域上第一过零点带宽为20Hz；要求时域图和频谱图在一个窗口显示，且都为双边幅度图

在MATLAB中，你可以通过`square`函数生成矩形脉冲信号，然后结合`plot`和`fft`函数来实现时域和频域的绘制。为了限制频域的第一过零点带宽为20Hz，我们需要对信号进行适当的低通滤波处理。以下是步骤：

1. **生成矩形脉冲信号**：

fs = 1000; % 采样频率，例如1kHz
duration = 1; % 信号持续时间，单位秒
t = (0:1/fs:duration)'; % 时间向量
pulse_width = 0.1; % 脉冲宽度，单位秒
rect_pulse = square(2*pi*t/pulse_width); % 矩形脉冲信号

2. **绘制时域图形**：

figure;
plot(t, rect_pulse, 'b', 'LineWidth', 2);
xlabel('Time (s)');
ylabel('Amplitude');
title('Rectangular Pulse Signal - Time Domain');
grid on;

3. **计算并绘制频谱图**：

% 对信号进行FFT并取绝对值，得到双边幅频谱
Y = abs(fft(rect_pulse));
freq = (0:length(Y)-1)*fs/(length(Y)-1); % 频率向量

% 设定截止频率为第一过零点对应的频率减去5Hz（以防边缘效应）
cutoff_freq = fs/2 - pulse_width*fs/2 - 5; % 第一过零点带宽中心减去5Hz
low_pass_filter = freq < cutoff_freq; % 切除高频部分
Y_filtered = Y .* low_pass_filter; % 应用滤波

% 双边幅度谱
Y_bilateral = Y_filtered + conj(flipud(Y_filtered)); % 左右对称

figure;
plot(freq, Y_bilateral, 'r', 'LineWidth', 2);
xlabel('Frequency (Hz)');
ylabel('|Magnitude|');
title('Rectangular Pulse Signal - Frequency Domain (Low-pass filtered)');
xlim([0 fs/2]);
grid on;

4. **在同一窗口显示时域图和频谱图**：

% 将两个图放在同一个窗口中
subplot(2,1,1);
hold on;
plot(t, rect_pulse, 'b');
hold off;

subplot(2,1,2);
plot(freq, Y_bilateral, 'r');
hold off;