利用matlab软件实现周期矩形脉冲信号的分解，绘制出直流分量及1、3、5、7、9次谐波的波形图。

在MATLAB中，可以使用傅里叶变换（Discrete Fourier Transform, DFT）来分解周期矩形脉冲信号，获取其各次谐波成分。下面是步骤：

1. **生成周期矩形脉冲信号**：首先创建一个周期矩形脉冲信号，比如一个周期为T的矩形波，通常可以通过`rectpulse`函数实现。

T = 1; % 设定周期
n = 0:0.01:T*100; % 时间点
pulsed_signal = rectpulse(n/T, 0.5*T);

这里的`0.5*T`代表矩形波的宽度。

2. **计算DFT**：使用`fft`函数对信号进行离散傅立叶变换，得到频率域的数据。

dft_result = fft(pulsed_signal);

3. **提取直流分量和各次谐波**：直流分量是DFT结果的第一个元素，奇数次谐波对应的是非零偶数次项（从第二个开始），偶数次谐波对应的是非零奇数次项。

dc_component = dft_result(1);
harmonics = abs(dft_result(2:end)); % 提取绝对值，表示幅度大小
odd_harmonics = harmonics(2:2:end); % 奇数次谐波
even_harmonics = harmonics(1:2:end); % 偶数次谐波

4. **绘制波形图**：最后，分别绘制原始信号、直流分量以及指定次数的谐波。

% 绘制原信号
plot(n, pulsed_signal, 'LineWidth', 1.5)
title('Periodic Rectangular Pulse Signal')

% 绘制直流分量
hold on
plot(n, [ones(size(n))/sqrt(size(n)), zeros(size(n)-1)], 'r--', 'LineWidth', 1.5)
text(0, 0.8, 'DC Component', 'HorizontalAlignment', 'center')

% 绘制部分谐波
for i = 1:5
    if mod(i, 2) == 1
        plot(n, harmonic(i)*cos(2*pi*i*n/T), 'LineWidth', 1.5)
    end
end
title(sprintf('Pulse Signal with DC and First 5 Odd Harmonics'))

完成上述步骤后，你将得到周期矩形脉冲信号及其组成部分的图形。