方波的合成实验（使用subplot函数） 方波信号可以分解为 用前5项、10项、15项谐波近似合成一频率为50HZ.幅度为3的方波,写出相应MATLAB程序并给出结果，并简述吉布斯现象

以下是MATLAB程序及结果:

t = linspace(0,1,1000);
f50 = 50;
A = 3;
y = zeros(size(t));
for n = 1:5
    y = y + (4*A/pi)/(2*n-1)*sin(2*pi*f50*(2*n-1)*t);
end
subplot(3,1,1)
plot(t,y)
title('Using 5 harmonics')

y = zeros(size(t));
for n = 1:10
    y = y + (4*A/pi)/(2*n-1)*sin(2*pi*f50*(2*n-1)*t);
end
subplot(3,1,2)
plot(t,y)
title('Using 10 harmonics')

y = zeros(size(t));
for n = 1:15
    y = y + (4*A/pi)/(2*n-1)*sin(2*pi*f50*(2*n-1)*t);
end
subplot(3,1,3)
plot(t,y)
title('Using 15 harmonics')

吉布斯现象是指信号经过傅立叶变换后，在频率为信号主频的附近会产生一个明显的振铃现象，振铃的幅度在理论上可以无限大。在方波的情况下，吉布斯现象表现为合成波在方波转折处出现明显的振铃。