Filtro seno de fourier

A continuación el programa para filtro Fourier modificando los parámetros que se consideraron
%Filtro Fourier
%Funciones seno y coseno de valores enteros de cuatro dígitos
close all;

clearall;
ciclos=input('Introduce el número de ciclos ha analizar: ');
f=20; frecuencia de 20 hz
T=1/f;
N=28; 28 muestras
h=T/N;
t=h:h:ciclos*T;
t1=h:h:ciclos*T-28*h;
Vm=5;amplitude de la señal
w=2*pi*f;
fseno=Vm*sin(w*t);
f=fseno+sin(7*w*t); contaminada con séptimos armónicos
subplot (3,1,1)
plot(t,fseno)
grid;
title('Señal senoidal pura')

subplot(3,1,2)
plot(t,f)
grid;
title('Señal senoidal + septimo armónico ')

for k=1:ciclos*N-28ir(k)=[((f(k)-f(14+k)+f(28+k))*1+(f(1+k)-f(13+k)-f(15+k)+f(27+k))*.97492+(f(2+k)-f(12+k)-f(16+k)+f(26+k))*.90096+(f(3+k)-f(11+k)-f(17+k)+f(25+k))*.7818+(f(4+k)-f(10+k)-f(18+k)+f(24+k))*.6234+(f(5+k)-f(9+k)-f(19+k)+f(23+k))*.4338 +(f(6+k)-f(8+k)-f(20+k)+f(22+k))*.2225+(f(7+k)+f(21))*0)]*.0714;ii(k)=[((f(1+k)+f(13+k)-f(15+k)-f(27+k))*.2225+(f(2+k)+f(12+k)-f(16+k)-f(26+k))*.4338+(f(4+k)+f(11+k)-f(17+k)-f(25+k))*.6243+(f(5+k)+f(10+k)-f(18+k)-f(24+k))*.78183+(f(6+k)+f(9+k)-f(19+k)-f(23+k))*.90096+(f(7+k)+f(8+k)-f(20+k)-f(22+k))*.97492+(f(8+k)-f(21+k))*1)]*.0714;

ip(k)=sqrt((ir(k))^2+(ii(k))^2);
end

subplot (3,1,3)
plot(t1,ip);
plot(ip);
grid;
title('Magnitud de la fundamentalcon Fourier')

Muestras tomadas de matlab
1. 2.1126
2. 2.1694
3. 2.1174
4. 3.9092
5. 5.5048
6. 4.8746
7. 4.0000
8. 4.8746
9. 5.504810. 3.9092

11. 2.1174
12. 2.1694
13. 2.1126
14. 0.0000
15. -2.1126
16. -2.1694
17. -2.1174
18. -3.9092
19. -5.5048
20. -4.874621. -4.0000
22. -4.8746
23. -5.5048
24. -3.9092
25. -2.1174
26. -2.1694
27. -2.1126
28. -0.0000

Estos valores se afectan con el seno y coseno de...