产生指定频率的非周期矩形脉冲，并做频谱分析

要生成指定频率的非周期矩形脉冲（通常称为方波），你可以使用Matlab中的`square`函数或自己编写循环生成信号。这里我会提供一个简单的例子如何用`square`函数生成方波，然后进行快速傅立叶变换（FFT）来进行频谱分析。

首先，你需要确定信号的参数，比如采样率、持续时间、频率等。假设你想生成50Hz的方波，采样率为1kHz（即每秒1000个样本点）。

% 参数设定
fs = 1000;   % 采样率 (Hz)
duration = 1;  % 持续时间 (秒)
freq = 50;     % 频率 (Hz)

% 生成信号
t = 0:1/fs:duration-1/fs;  % 时间轴
signal = square(2*pi*freq*t); % 方波信号

% FFT分析
N = length(signal);
Y = fft(signal, N);  % 傅立叶变换
frequencies = (0:N-1)*(fs/N);  % 频率轴

% 取绝对值并绘制频谱
Pxx = abs(Y/N);  % 平均功率谱密度
Pxx_db = 20*log10(Pxx);  % 对数尺度上的谱密度
plot(frequencies, Pxx_db);
xlabel('Frequency (Hz)');
ylabel('Power spectral density (dB/Hz)');
title(sprintf('Spectrum of a Square Wave at %.1f Hz', freq));

执行这段代码后，你会得到一个显示了方波信号频谱的图表，其中50Hz处应有一个明显的峰值。