armónicos de una señal cuadrada

n armónicos de una señal cuadrada de periodo 0.2s y amplitud 1.

% Desarrollo en serie de Fourier de una señal cuadrada de periodo 0.2s y amplitud 1
clear;
% frecuencia de la señal cuadrada (=1/T)
f=5;
T=1/f;
% Indice de los coeficientes
n=1:10;
% Generamos la serie de Fourier
t=-1:0.01:1; % vector de tiempos
for i=1:50
for k=1:size(t,2)
s(i,k)=(2*(1-cos(pi*i))/(pi*i))*sin(2*pi*i*f*t(k));
end
end
for k=1:size(t,2)
st(k)=sum(s(:,k));
end
st(1)=st(1)+1;
plot(t,st,'r');
hold on;
% Señal cuadrada original
f_cuadrada=square(2*pi*f*t,50);
plot(t,f_cuadrada);
xlabel('tiempo');
ylabel('Amplitud');