heaviside' 需要 Symbolic Math Toolbox。

非常抱歉，我之前提供的方法需要使用 Signal Processing Toolbox 中的 `heaviside` 函数，而非 Symbolic Math Toolbox 中的 `heaviside` 函数。如果您的 MATLAB 没有安装 Signal Processing Toolbox，可以使用以下方法生成阶跃信号：

% 设定采样频率和采样点数
fs = 1000; % 采样频率为1000Hz
N = 1024; % 采样点数为1024

% 生成阶跃信号
t = 0:1/fs:(N-1)/fs;
x = zeros(size(t));
x(t>=0) = 1;

% 计算频谱
X = fft(x,N);
f = (0:N-1)*fs/N;

% 画出频谱图
plot(f,abs(X));
xlabel('频率 (Hz)');
ylabel('幅度');
title('阶跃信号的频谱');

这段代码使用了 `zeros` 函数生成一个全零向量，然后将 $t>=0$ 的部分赋值为 $1$，即得到阶跃信号。接下来计算频谱并画出频谱图。运行上述代码，即可得到阶跃信号的频谱图。