使用matlab编写程序，对T=4，τ=2，E=1的周期矩形脉冲信号进行分解

好的，下面是用MATLAB编写的程序，用于对T=4，τ=2，E=1的周期矩形脉冲信号进行分解：

T = 4;  % 周期
tau = 2;  % 脉宽
E = 1;  % 振幅

% 定义时间范围
t = 0:0.01:T;

% 计算矩形脉冲信号
rect_pulse = E * rectpuls(mod(t, T), tau);

% 傅里叶级数展开
n = 1:20;  % 谐波次数范围
harmonics = zeros(size(t));
for i = n
    coefficient = (2 * E / i / pi) * sin(2 * pi * i * tau / T);
    harmonics = harmonics + coefficient * cos(2 * pi * i * t / T);
end

% 绘制结果
subplot(2, 1, 1);
plot(t, rect_pulse);
title('矩形脉冲信号');
xlabel('时间');
ylabel('幅值');

subplot(2, 1, 2);
plot(t, harmonics);
title('谐波分量');
xlabel('时间');
ylabel('幅值');

运行这段代码后，你将会看到两个子图，第一个子图显示的是矩形脉冲信号，第二个子图显示的是谐波分量。