MATLAB奇数谐波合成方波

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

2、第二，启动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',[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的矩阵。