使用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
    harmonics = harmonics + (2 * E / i / pi) * sin(2 * pi * i * t / T);
end

% 合成方波信号
square_wave = sign(harmonics);

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

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

subplot(3, 1, 3);
plot(t, square_wave);
title('合成方波信号');
xlabel('时间');
ylabel('幅值');

运行这段代码后，你将会看到三个子图，分别为矩形脉冲信号、谐波分量和合成的方波信号。