一个周期矩阵脉冲信号的matlab,频率是步进频的

下面是一个周期矩阵脉冲信号的 Matlab 代码，其中频率是步进频率：

% 设置参数
N = 200; % 采样点数
f0 = 100; % 基频
n = 5; % 周期数
step = 5; % 步进频率

% 生成时间序列
t = linspace(0, (n-1)*1/f0, N*n);

% 生成信号
x = zeros(1, N*n);
for i = 1:n
    x((i-1)*N+1:i*N) = rectpuls(t((i-1)*N+1:i*N));
end

% 绘制时域波形
subplot(2,1,1);
plot(t, x);
title('时域波形');
xlabel('时间 / s');
ylabel('幅值');

% 进行傅里叶变换
X = fft(x);
X = fftshift(X);

% 生成频率序列
fs = N*(n-1)/t(end);
f = linspace(-fs/2, fs/2, N*n);

% 绘制频域波形
subplot(2,1,2);
plot(f, abs(X));
title('频域波形');
xlabel('频率 / Hz');
ylabel('幅值');
axis([-fs/2, fs/2, 0, max(abs(X))*1.1]);

% 添加步进频率
f_step = f0:step:f0*(n-1);
for i = 1:length(f_step)
    line([f_step(i), f_step(i)], [0, max(abs(X))*1.1], 'Color', 'r');
end

该代码将生成一个周期矩阵脉冲信号，并绘制出时域波形和频域波形。在频域波形中，还将添加步进频率的红色竖线。你可以根据需要更改参数来自定义你的信号。