你能告诉我用matlab代码绘制周期为N，宽度分别为2N1和2N2 的离散时间周期方波序列的时域波形和频谱（频谱主瓣中心位于频率等于0处）吗

当然可以。在MATLAB中，你可以使用`sin()`函数来创建周期性信号，然后利用`plot()`画出时域波形，通过`fft()`计算频谱并进一步分析。这里是一个简单的示例：

% 定义参数
N = 10; % 周期
N1 = N / 4; % 第一个波峰位置
N2 = 3 * N / 4; % 第二个波峰位置

% 创建离散时间周期方波序列
t = 0:N-1; % 时间轴
signal = [zeros(1, floor(N1)); ones(1, round(2*N1)); zeros(1, ceil(N2 - N1) - 1)]; % 方波序列

% 绘制时域波形
figure;
stem(t, signal);
title('离散时间周期方波时域波形');
xlabel('时间 (样本)');
ylabel('幅度');

% 计算频谱并显示主瓣中心
fs = N; % 采样率
yf = fft(signal); % 对信号进行傅立叶变换
Pxy = abs(yf).^2 / length(signal); % 平均功率谱密度
[frequencies, spectrum] = freqz(ones(size(signal))/length(signal), [], fs); % 频谱
figure;
plot(frequencies, Pxy);
xlim([-fs/2 fs/2]); % 窗口显示整个频谱
title('离散时间周期方波频谱');
xlabel('频率 (Hz)');
ylabel('幅度 (dB)');

% 验证主瓣中心
main_peak_index = findPeaks(Pxy, 'MinPeakDistance', fs/(2*N), 'MinPeakHeight', max(Pxy)/2); % 主瓣中心的索引
if ~isempty(main_peak_index)
    main_freq = frequencies(main_peak_index); % 主瓣中心频率
    fprintf('主瓣中心频率位于 %.2f Hz\n', main_freq);
end

在这个例子中，我们首先定义了周期、波峰位置，然后创建了一个离散时间方波序列。接着，分别绘制了时域波形和频谱，并找到了频谱的主瓣中心。