%Łukasz Ibron, I03
%PM13P
%Cw 3 - Szereg Furiera
clear all,close all
format short

t=(0:999)/1000;
ft=[1-cos(pi*t(1:500)),zeros(1,500)];


figure(1)
plot(t,ft),grid
xlabel('t')
ylabel('f(t)')
title('zad.3')


f0=pi-2/2*pi
k=(1:20).';
fk=(-1).^k./(pi*(4*k.^2-1))+j*(((4*k.^2-1)*(-1).^k+1)./(2*pi*k.*(4*k.^2-1)));
disp('   k       |Fk|     arg{Fk}   [deg]')
[k, abs(fk),angle(fk)*180/pi]

x=f0;
figure(2)
for l=1:20
   x=x+2*abs(fk(l))*cos(l*2*pi*t+angle(fk(l)));
   plot(t,ft,t,x,'-m'),grid
   xlabel('t')
   ylabel('f(t)')
   title(['ku=u ',num2str(l)])
   pause
end