![ec9afebd0a80f45d2e2e81ed9caaaecf.png](https://i-blog.csdnimg.cn/blog_migrate/f510f28ce0b22387101e2bfd4cec67f4.jpeg)
从正弦波转换为方波(傅里叶验证)
此示例说明方波的傅里叶级数展开式是如何由奇次谐波的和构成的。
1.首先以 0.1 为步长,生成一个从 0 到 10 的时间向量,并求出所有点的正弦。绘制基频图
t = 0:.1:10;
y = sin(t);
plot(t,y);
![63fc35e4b259531b64bf796b4acd7a23.png](https://i-blog.csdnimg.cn/blog_migrate/2b5ef131ab17e47eff7a8406d6ab8fe2.jpeg)
2.接下来,向基频添加第三个谐波,并绘制谐波图
y = sin(t) + sin(3*t)/3;
plot(t,y);
![32872f08d275f464684ca4444824bb93.png](https://i-blog.csdnimg.cn/blog_migrate/0a09ee225c5889f073cbb1b9097ed577.png)
3.接下来使用第一、第三、第五、第七和第九个谐波
y = sin(t) + sin(3*t)/3 + sin(5*t)/5 + sin(7*t)/7 + sin(9*t)/9;
plot(t,y);
![f40bc6d92d36fa2bb80aadffd7380f2b.png](https://i-blog.csdnimg.cn/blog_migrate/ddc71d3c89396c5defaa4556ed07839a.jpeg)
4.最后,从基频开始创建更多连续谐波的向量,一直到第 19 个谐波为止,并将所有中间步长保存为矩阵的行。
在同一个图窗中绘制这些向量,以便显示方波的演变。请注意,吉布斯效应表明它实际上永远不会转换为方波。
t = 0:.02:3.14;
y = zeros(10,length(t));
x = zeros(size(t));
for k = 1:2:19
x = x + sin(k*t)/k;
y((k+1)/2,:) = x;
end
plot(y(1:2:9,:)')
title('The building of a square wave: Gibbs'' effect')
![347d7928447b6ea5df175a562f1eac68.png](https://i-blog.csdnimg.cn/blog_migrate/c631adfb7998dc1f03be8df1aed5e5cf.jpeg)
5.下面提供了一个三维曲面图,该曲面图表示正弦波到方波的逐变过程
surf(y);
shading interp
axis off ij
![6767b9fbc7420a98c81c5e627140a05c.png](https://i-blog.csdnimg.cn/blog_migrate/d35abf4e178021ed86fd1bf22155f896.jpeg)
![984e8e5879f76da6582ba56b67fa3011.png](https://i-blog.csdnimg.cn/blog_migrate/2c06c9274d5131c8482b510550b911bb.jpeg)