MATLAB奇数谐波合成方波

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

2、第二，启动MATLAB，新建脚本（Ctrl+N），输入如下代码：close all; clear all; 艘早祓胂clct=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',[50,50,700,600],'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',[100,100,700,600],'Name','2')A=[y1;y2;y3;y4];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]是将y1，y2，y3，y4合成一个矩阵，B=A'是指矩阵B为矩阵A的转置。

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

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

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