第一，sin((2n-1)*t)/(2n-1)累加，将sin函数变成方波，具体公式如下图。

第二，启动MATLAB，新建脚本（Ctrl+N），输入如下代码：

close all; clear all; clc

t=0:.1:2*pi;

y1=sin(t);

y2=sin(t)+sin(3*t)/3;

y3=sin(t)+sin(3*t)/3+sin(5*t)/5;

y4=sin(t)+sin(3*t)/3+sin(5*t)/5+sin(7*t)/7;

figure(Position,,Name,1)

subplot(2,2,1);plot(y1); subplot(2,2,2);plot(y2);

subplot(2,2,3);plot(y3); subplot(2,2,4);plot(y4);

figure(Position,,Name,2)

A=;B=A;

plot(B,LineWidth,2)

axis tight off

其中figure(...Name,1)下面两行代码是将y1，y2，y3，y4分别绘制在图1中，figure(...Name,2)下面两行代码是将y1，y2，y3，y4合成绘制在图2中。

A=是将y1，y2，y3，y4合成一个矩阵，B=A是指矩阵B为矩阵A的转置。

第三，保存和运行上述代码，得到图形figure1，y1，y2，y3，y4分别绘制在图1中，并逐步变成方波。

第四，同时得到图形figur,2，y1，y2，y3，y4合成绘制在图2中。

第五，通过双击工作区（Workspace）中的变量A和B，可以查看A，B的数据情况，可以清楚看到A为4*63的矩阵，B是A的转置，B为64*3的矩阵。

A square wave is made up of a sum of odd harmonics，奇数谐波合成方波。