使用matlab，将方波信号的前五次谐波近似合成50HZ，幅值为3的方波

以下是使用matlab实现该操作的代码：

% 定义方波信号
t = 0:0.0001:1; % 时间范围为0到1秒，采样率为10000Hz
f = 50; % 方波频率为50Hz
x = 1.5*square(2*pi*f*t,50)-0.5; % 幅值为3的方波

% 计算前五次谐波的频率和幅值
f_harm = f*(1:5); % 谐波频率
a_harm = [4/pi, 1/3/pi, 4/15/pi, 1/5/pi, 4/35/pi]; % 谐波幅值

% 合成近似信号
x_approx = zeros(size(x)); % 初始化近似信号
for i = 1:5
    x_harm = a_harm(i)*sin(2*pi*f_harm(i)*t); % 计算谐波信号
    x_approx = x_approx + x_harm; % 叠加谐波信号
end
x_approx = 3*x_approx/max(abs(x_approx)); % 幅值调整为3

% 绘制原始信号和近似信号
figure;
subplot(2,1,1);
plot(t,x);
title('原始信号');
xlabel('时间（秒）');
ylabel('幅值');
subplot(2,1,2);
plot(t,x_approx);
title('近似信号');
xlabel('时间（秒）');
ylabel('幅值');

运行上述代码后，将会得到如下图所示的绘图结果：