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;endplot(y(1:2:9,:)')title('The building of a square wave: Gibbs'' effect')
8、按“Enter”键。得图4。
9、这是一个3D曲面，代表正弦波向方波的逐渐转换。键入：surf(y);shading interpaxis off ij
10、按“Enter”键。得图5。
