matlab拟合简单方波傅里叶级数(附代码和计算过程)

如图所示是一个周期为2π的方波，分析如下：

1. 该函数是奇函数，故an余弦项全部为0。

2.该函数是奇谐函数，故只含有奇次谐波。

下面进行计算：

t = -5:0.001:5;         %创建时间轴从-5到5，步长为0.001
y1 = square(t);         %创建方波，默认周期为2π
subplot(2,1,1);         
plot(t,y1);
ylim([-2 2]);


%用for循环实现累加,n越大谐波次数越多
y2 = 0;
for n=1:200
    y2 = y2 + (1/n) * (1-(-1).^n) * sin(n*t);
end

subplot(2,1,2);
plot(t,2/pi*y2);
ylim([-2 2]);

n=200时运行结果如下所示：