1、我们首先形成一个以0.1为步长从0到10的时间向量，并取所有点的正弦值。
让我们绘制这个基本频率。
命令窗口键入：
t = 0:.1:10;
y = sin(t);
plot(t,y);
2、按“Enter”键。
得图1。
3、现在，将三次谐波添加到基波，并将其绘制出来。
键入：
y = sin(t) + sin(3*t)/3;
plot(t,y);
4、按“Enter”键。
得图2。
5、现在使用第一，第三，第五，第七和第九谐波。
键入：
y = sin(t) + sin(3*t)/3 + sin(5*t)/5 + sin(7*t)/7 + sin(9*t)/9;
plot(t,y);
6、按“Enter”键。
得图3。
7、最后，我们将从基波谐波到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')
8、按“Enter”键。
得图4。
9、这是一个3D曲面，代表正弦波向方波的逐渐转换。
键入：
surf(y);
shading interp
axis off ij
10、按“Enter”键。
得图5。
